def read(self, indexes=None, **kwargs):
band_indexes = self._get_band_indexes(indexes)<EOL>arr = self.process.get_raw_output(self.tile)<EOL>if len(band_indexes) == <NUM_LIT:1>:<EOL><INDENT>return arr[band_indexes[<NUM_LIT:0>] - <NUM_LIT:1>]<EOL><DEDENT>else:<EOL><INDENT>return ma.concatenate([ma.expand_dims(arr[i - <NUM_LIT:1>], <NUM_LIT:0>) for i in band_indexes])<EOL><DEDENT>
Read reprojected & resampled input data.

Parameters
----------
indexes : integer or list
    band number or list of band numbers

Returns
-------
data : array