code: def read(self, **keys):
          
          if not self.has_data():
              return None
  
          dtype, shape = self._get_dtype_and_shape()
          array = numpy.zeros(shape, dtype=dtype)
          self._FITS.read_image(self._ext+1, array)
          return array
docs: Read the image.
  
          If the HDU is an IMAGE_HDU, read the corresponding image.  Compression
          and scaling are dealt with properly.