code: def interpolate_bilinear_skycoord(self, skycoord, values):
          
          if self.frame is None:
              raise NoFrameError("interpolate_bilinear_skycoord")
          skycoord = skycoord.transform_to(self.frame)
          representation = skycoord.represent_as(UnitSphericalRepresentation)
          lon, lat = representation.lon, representation.lat
          return self.interpolate_bilinear_lonlat(lon, lat, values)
docs: Interpolate values at specific celestial coordinates using bilinear interpolation.
  
          If a position does not have four neighbours, this currently returns NaN.
  
          Note that this method requires that a celestial frame was specified when
          initializing HEALPix. If you don't know or need the celestial frame, you
          can instead use :meth:`~astropy_healpix.HEALPix.interpolate_bilinear_lonlat`.
  
          Parameters
          ----------
          skycoord : :class:`~astropy.coordinates.SkyCoord`
              The celestial coordinates at which to interpolate
          values : `~numpy.ndarray`
              1-D array with the values in each HEALPix pixel. This must have a
              length of the form 12 * nside ** 2 (and nside is determined
              automatically from this).
  
          Returns
          -------
          result : `~numpy.ndarray`
              1-D array of interpolated values