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/patterns/game_of_life.py

@@ -13,10 +13,10 @@
 
 The rules as taken from https://en.wikipedia.org/wiki/Conway%27s_Game_of_Life:
 
- 1. Any live cell with fewer than two live neighbours dies, as if by underpopulation.
- 2. Any live cell with two or three live neighbours lives on to the next generation.
- 3. Any live cell with more than three live neighbours dies, as if by overpopulation.
- 4. Any dead cell with exactly three live neighbours becomes a live cell, as if by reproduction.
+1. Any live cell with fewer than two live neighbours dies, as if by underpopulation.
+2. Any live cell with two or three live neighbours lives on to the next generation.
+3. Any live cell with more than three live neighbours dies, as if by overpopulation.
+4. Any dead cell with exactly three live neighbours becomes a live cell, as if by reproduction.
 
 Initial conditions
 ------------------

gridkit-py/examples/raster_operations/coordinate_transformations.py

@@ -56,7 +56,7 @@
 
 utm_epsg = 32629
 transformer = Transformer.from_crs(
-    data_tile_wgs84.grid.crs, CRS.from_user_input(utm_epsg), always_xy=True
+    data_tile_wgs84.get_grid().crs, CRS.from_user_input(utm_epsg), always_xy=True
 )  # UTM zone 28N
 
 wgs84_centroids_utm = transformer.transform(

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(