code: def from_fits_with_pixel_scale(cls, file_path, hdu, pixel_scales, origin=(0.0, 0.0)):
          
          return cls(array_util.numpy_array_2d_from_fits(file_path, hdu), pixel_scales, origin)
docs: Loads the image from a .fits file.
  
          Parameters
          ----------
          file_path : str
              The full path of the fits file.
          hdu : int
              The HDU number in the fits file containing the hyper.
          pixel_scales: (float, float)
              The arc-second to pixel conversion factor of each pixel.