Read image unit

The unit is read from the ``BUNIT`` keyword if present. If the ``unit`` argument to ``read`` is provided **and** there is a a unit in the FITS header, the ``unit`` argument is used and an INFO-level message is issued.