def read(self, **keys):
if not self.has_data():<EOL><INDENT>return None<EOL><DEDENT>dtype, shape = self._get_dtype_and_shape()<EOL>array = numpy.zeros(shape, dtype=dtype)<EOL>self._FITS.read_image(self._ext+<NUM_LIT:1>, array)<EOL>return array<EOL>
Read the image.

If the HDU is an IMAGE_HDU, read the corresponding image.  Compression
and scaling are dealt with properly.