Topics

Windowed reads and coordinates

Koshy Thomas
 

Hi all

I'm trying to create vector features out of raster files using the `features.shapes` module and I'm running into a weird issue where when I run a windowed read on a `.tif`, the Y coordinates of the polygons are thrown off while the X coordinates seem correct. If I run the same logic without a windowed read (against the whole tif), it seems to work correctly. Any ideas? I'm attaching the .tif and the code screenshot is below:

Koshy Thomas
 

Turns out the transform object needs to actually be

transform = rasterio.windows.transform(window, dataset.transform)

for windowed reads.