Fix .grid call that escaped replacing to .get_grid()

Author: tmillenaar

gridkit-py/examples/cell_centric_operations/downslope_path.py

@@ -58,7 +58,7 @@
 def downslope_path(cell_id):
     visited_cells.append(cell_id)
     # add new neigbours and their values to their respective lists
-    neighbours = dem.grid.neighbours(cell_id, connect_corners=True)
+    neighbours = dem.get_grid().neighbours(cell_id, connect_corners=True)
     for cell, value in zip(neighbours, dem.value(neighbours)):
         cell = tuple(cell.index)
         if not cell in visited_cells and not cell in all_neighbours:
@@ -80,7 +80,7 @@ def downslope_path(cell_id):
 
 # %%
 # Determine the starting cell and call the function
-start_cell_id = dem.grid.cell_at_point((29330, 167848))
+start_cell_id = dem.get_grid().cell_at_point((29330, 167848))
 downslope_path(tuple(start_cell_id.index))
 
 

gridkit-py/examples/raster_operations/ndvi.py

@@ -61,6 +61,6 @@
 fig.colorbar(im_ndvi, ax=ax, fraction=0.022, pad=0.01)
 ax.set_xlabel("lon")
 ax.set_ylabel("lat")
-ax.set_title(f"NDVI of scene in the alps \n EPSG:{ndvi.grid.crs.to_epsg()}")
+ax.set_title(f"NDVI of scene in the alps \n EPSG:{ndvi.get_grid().crs.to_epsg()}")
 
 plt.show()

gridkit-py/examples/raster_operations/resampling.py

@@ -25,7 +25,7 @@
 dem = raster_to_data_tile(
     "../../tests/data/alps_dem.tiff", bounds=(28300, 167300, 28700, 167700)
 )
-print(f"Original resolution in (dx, dy): ({dem.grid.dx}, {dem.grid.dy})")
+print(f"Original resolution in (dx, dy): ({dem.get_grid().dx}, {dem.get_grid().dy})")
 
 
 # %%
@@ -49,8 +49,10 @@
 # Let's define a rectangular grid with a cell size of 10x10 degrees.
 # I am calling it rectdem for later on I will define hexagonal ones as well.
 # To make sure it worked we can print out the cellsize after resampling
-rectdem = dem.resample(RectGrid(dx=10, dy=10, crs=dem.grid.crs), method="linear")
-print(f"Downsampled resolution in (dx, dy): ({rectdem.grid.dx}, {rectdem.grid.dy})")
+rectdem = dem.resample(RectGrid(dx=10, dy=10, crs=dem.get_grid().crs), method="linear")
+print(
+    f"Downsampled resolution in (dx, dy): ({rectdem.get_grid().dx}, {rectdem.get_grid().dy})"
+)
 
 # %%
 #
@@ -71,10 +73,12 @@
 # for a fair visual comparison.
 #
 hexdem_flat = dem.resample(
-    HexGrid(area=rectdem.grid.area, shape="flat", crs=dem.grid.crs), method="linear"
+    HexGrid(area=rectdem.get_grid().area, shape="flat", crs=dem.get_grid().crs),
+    method="linear",
 )
 hexdem_pointy = dem.resample(
-    HexGrid(area=rectdem.grid.area, shape="pointy", crs=dem.grid.crs), method="linear"
+    HexGrid(area=rectdem.get_grid().area, shape="pointy", crs=dem.get_grid().crs),
+    method="linear",
 )
 
 

gridkit-py/gridkit/tile.py

@@ -995,7 +995,7 @@ def resample(self, alignment_grid, method="nearest", **interp_kwargs):
         tile_is_given = isinstance(alignment_grid, Tile)
         if tile_is_given:
             tile = alignment_grid
-            alignment_grid = alignment_grid.grid
+            alignment_grid = alignment_grid.get_grid()
 
         if self.get_grid().crs is None or alignment_grid.crs is None:
             warnings.warn(