code: def datapt_to_wcspt(self, datapt, coords='data', naxispath=None):
          
          # We provide a list comprehension version for WCS packages that
          # don't support array operations.
          if naxispath:
              raise NotImplementedError
  
          return np.asarray([self.pixtoradec((pt[0], pt[1]), coords=coords)
                             for pt in datapt])
docs: Convert multiple data points to WCS.
  
          Parameters
          ----------
          datapt : array-like
              Pixel coordinates in the format of
              ``[[x0, y0, ...], [x1, y1, ...], ..., [xn, yn, ...]]``.
  
          coords : 'data' or None, optional, default to 'data'
              Expresses whether the data coordinate is indexed from zero.
  
          naxispath : list-like or None, optional, defaults to None
              A sequence defining the pixel indexes > 2D, if any.
  
          Returns
          -------
          wcspt : array-like
              WCS coordinates in the format of
              ``[[ra0, dec0], [ra1, dec1], ..., [ran, decn]]``.