code: def numpy_array_2d_from_fits(file_path, hdu):
      
      hdu_list = fits.open(file_path)
      return np.flipud(np.array(hdu_list[hdu].data))
docs: Read a 2D NumPy array to a .fits file.
  
      After loading the NumPy array, the array is flipped upside-down using np.flipud. This is so that the arrays \
      appear the same orientation as .fits files loaded in DS9.
  
      Parameters
      ----------
      file_path : str
          The full path of the file that is loaded, including the file name and '.fits' extension.
      hdu : int
          The HDU extension of the array that is loaded from the .fits file.
  
      Returns
      -------
      ndarray
          The NumPy array that is loaded from the .fits file.
  
      Examples
      --------
      array_2d = numpy_array_from_fits(file_path='/path/to/file/filename.fits', hdu=0)