update transformation metadata when cropto used in read_geotif

danforthcenter · Sep 20, 2024 · b878289 · b878289
1 parent 0648992
commit b878289
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/plantcv/geospatial/read_geotif.py b/plantcv/geospatial/read_geotif.py
@@ -90,8 +90,9 @@ def _read_geotif_and_shapefile(filename, cropto):
                 shapes = [mapping(convex_hull)]
         # rasterio does the cropping within open
         with rasterio.open(filename, 'r') as src:
-            img_data, _ = mask(src, shapes, crop=True)
+            img_data, trans_metadata = mask(src, shapes, crop=True)
             metadata = src.meta.copy()
+            metadata.update({"transform": trans_metadata})
             d_type = src.dtypes[0]
 
     else: