Inverted axis in numpy


gabriel@...
 

Hi everyone. I started using rasterio a few days ago to read some DSM/DTM and process them using numba.

I was trying to debug some weird behaviour in my numba function and i realized that the "read" function in Rasterio actually flips the axis of my dataset.
The file I am opening is a Geotiff, with EPSG:3003 projection.

This how i reproduced this behaviour:

dataset = rasterio.open("raster.tif")
band = dataset.read(1)
print(dataset.width)    -> 84
print(dataset.heigth)  -> 44
print(band.shape)   -> (44,84)


Is this behaviour normal? I don't understand why my data are being flipped in this way.

Thanks, Gabriel

Join main@rasterio.groups.io to automatically receive all group messages.