GeoRaster2.get_tile messes up masked values in some cases

[astrojuanlu]

This code:

import telluric as tl
rs2 = tl.GeoRaster2.open("../tests/data/raster/overlap2.tif", lazy_load=False)
rs2.get_tile(4377, 3038, 13)

produces this image in the notebook:

This is what @guydou said about it in #108:

about the mask, the mask will not have an effect on a raster that nodata value set, the paramter tells rasterio not to read the mask from the raster, we are doing it, because when you have no-data value set up rasterio is not using the overview masks

[drnextgis]

This is consequence of having nodata=0 in GeoRaster2 constructor. If we change this row:

telluric/telluric/georaster.py

Line 1598 in 9cabe2a

raster = self.copy_with(image=array, affine=affine)

to:

raster = self.copy_with(image=array, affine=affine, nodata=None)

then image looks as expected:

[guydou]

it looks a projection issue and not masking

crop with out projection

and then projection