Utilities for processing raster data
input_image = "../data/phl_ppp_2020_constrained.tif"
raster = rio.open(input_image)
raster.meta
show(raster.read(1), cmap="plasma", transform=raster.transform)
bbox = (120.888062, 14.394778, 121.199112, 14.705822)
circle_gdf
print(circle_gdf.crs)
fig, ax = plt.subplots(1, 1, figsize=(4, 8))
show(raster.read(1), cmap="viridis", ax=ax, transform=raster.transform)
circle_gdf.plot(ax=ax, facecolor="none", edgecolor="yellow")
ax
output_folder = Path("../data")
query_window_by_gdf(input_image, output_folder, circle_gdf, mask=False)
with rio.open(output_folder / "output_0.tif") as dst:
fig, ax = plt.subplots(1, 1, figsize=(4, 4))
show(dst.read(1), cmap="viridis", ax=ax, transform=dst.transform)
circle_gdf.plot(facecolor="none", edgecolor="yellow", ax=ax)
print(dst.read(1))
ax
grid_generator = grids.SquareGridGenerator(100_000)
grid_gdf = grid_generator.generate_grid(circle_gdf)
grid_gdf["name"] = (
"gridxy-" + grid_gdf["x"].astype(str) + "-" + grid_gdf["y"].astype(str)
)
grid_gdf
fig, ax = plt.subplots(1, 1, figsize=(4, 8))
show(raster.read(1), cmap="viridis", ax=ax, transform=raster.transform)
grid_gdf.plot(ax=ax, facecolor="none", edgecolor="yellow")
ax
query_window_by_gdf(input_image, output_folder, grid_gdf, name_col="name", mask=False)
for name in grid_gdf["name"]:
image_path = output_folder / (name + ".tif")
with rio.open(image_path) as dst:
fig, ax = plt.subplots(1, 1, figsize=(4, 4))
ax.set_title(image_path)
show(dst.read(1), cmap="viridis", ax=ax, transform=dst.transform)