code: def read(self):
          
          data = xr.open_rasterio(self.finfo["filename"],
                                  chunks=(1, CHUNK_SIZE, CHUNK_SIZE))
          attrs = data.attrs.copy()
          # Create area definition
          if hasattr(data, 'crs'):
              self.area = self.get_geotiff_area_def(data.crs)
  
          # Rename to Satpy convention
          data = data.rename({'band': 'bands'})
  
          # Rename bands to [R, G, B, A], or a subset of those
          data['bands'] = BANDS[data.bands.size]
  
          # Mask data if alpha channel is present
          try:
              data = mask_image_data(data)
          except ValueError as err:
              logger.warning(err)
  
          data.attrs = attrs
          self.file_content['image'] = data
docs: Read the image