extract edinburgh v1

import%20marimo%0A%0A__generated_with%20%3D%20%220.18.4%22%0Aapp%20%3D%20marimo.App(width%3D%22medium%22)%0A%0A%0A%40app.cell%0Adef%20_()%3A%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_()%3A%0A%20%20%20%20import%20geopandas%20as%20gpd%0A%20%20%20%20import%20pandas%20as%20pd%0A%20%20%20%20import%20numpy%20as%20np%0A%20%20%20%20import%20duckdb%0A%20%20%20%20from%20shapely%20import%20wkt%2C%20LineString%0A%20%20%20%20import%20matplotlib.pyplot%20as%20plt%0A%20%20%20%20import%20h3%0A%20%20%20%20from%20shapely.geometry%20import%20Polygon%2C%20shape%2C%20mapping%2C%20box%0A%20%20%20%20import%20rasterio%0A%20%20%20%20from%20rasterio.mask%20import%20mask%0A%20%20%20%20from%20typing%20import%20List%2C%20Tuple%0A%20%20%20%20from%20osgb%20import%20format_grid%0A%20%20%20%20import%20pyarrow%20as%20pa%0A%20%20%20%20from%20pyarrow%20import%20feather%0A%20%20%20%20return%20(%0A%20%20%20%20%20%20%20%20List%2C%0A%20%20%20%20%20%20%20%20Polygon%2C%0A%20%20%20%20%20%20%20%20box%2C%0A%20%20%20%20%20%20%20%20duckdb%2C%0A%20%20%20%20%20%20%20%20format_grid%2C%0A%20%20%20%20%20%20%20%20gpd%2C%0A%20%20%20%20%20%20%20%20h3%2C%0A%20%20%20%20%20%20%20%20mapping%2C%0A%20%20%20%20%20%20%20%20mask%2C%0A%20%20%20%20%20%20%20%20np%2C%0A%20%20%20%20%20%20%20%20pd%2C%0A%20%20%20%20%20%20%20%20rasterio%2C%0A%20%20%20%20%20%20%20%20wkt%2C%0A%20%20%20%20)%0A%0A%0A%40app.cell%0Adef%20_(duckdb)%3A%0A%20%20%20%20def%20duckdb_connection()%3A%0A%20%20%20%20%20%20%20%20conn%20%3D%20duckdb.connect()%0A%20%20%20%20%20%20%20%20conn.install_extension(%22spatial%22%2C%20force_install%3DFalse)%0A%20%20%20%20%20%20%20%20conn.load_extension(%22spatial%22)%0A%20%20%20%20%20%20%20%20return%20conn%0A%20%20%20%20return%20(duckdb_connection%2C)%0A%0A%0A%40app.cell%0Adef%20_(gpd%2C%20wkt)%3A%0A%20%20%20%20def%20query_to_dataframe(conn%2C%20query)%3A%0A%20%20%20%20%20%20%20%20arrow_table%20%3D%20conn.execute(query).fetch_arrow_table()%0A%20%20%20%20%20%20%20%20df%20%3D%20arrow_table.to_pandas()%0A%20%20%20%20%20%20%20%20gdf%20%3D%20gpd.GeoDataFrame(%0A%20%20%20%20%20%20%20%20%20%20%20%20df%2C%20geometry%3Ddf%5B%22geometry%22%5D.apply(wkt.loads)%2C%20crs%3D%22EPSG%3A4326%22%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20return%20gdf%0A%20%20%20%20return%20(query_to_dataframe%2C)%0A%0A%0A%40app.cell%0Adef%20_(query_to_dataframe)%3A%0A%20%20%20%20def%20load_city_boundary(conn%2C%20gers_id%2C%20release%3D%222026-02-18.0%22)%3A%0A%20%20%20%20%20%20%20%20query%20%3D%20f%22%22%22%0A%20%20%20%20%20%20%20%20%20%20%20%20SELECT%20*%20EXCLUDE(geometry)%2C%20ST_AsText(geometry)%20AS%20geometry%20%0A%20%20%20%20%20%20%20%20%20%20%20%20FROM%20read_parquet('s3%3A%2F%2Foverturemaps-us-west-2%2Frelease%2F%7Brelease%7D%2Ftheme%3Ddivisions%2Ftype%3Ddivision_area%2F*'%2C%20filename%3Dtrue%2C%20hive_partitioning%3D1)%0A%20%20%20%20%20%20%20%20%20%20%20%20WHERE%20id%20%3D%20'%7Bgers_id%7D'%0A%20%20%20%20%20%20%20%20%20%20%20%20%22%22%22%0A%20%20%20%20%20%20%20%20return%20query_to_dataframe(conn%2C%20query)%0A%20%20%20%20return%20(load_city_boundary%2C)%0A%0A%0A%40app.cell%0Adef%20_(duckdb_connection)%3A%0A%20%20%20%20conn%20%3D%20duckdb_connection()%0A%20%20%20%20return%20(conn%2C)%0A%0A%0A%40app.cell%0Adef%20_(conn%2C%20load_city_boundary)%3A%0A%20%20%20%20edinburgh_gers%20%3D%20%2258a34fa4-bc76-476e-81a8-1ed8a5cd693f%22%0A%20%20%20%20edinburgh_gdf%20%3D%20load_city_boundary(conn%2C%20edinburgh_gers)%0A%20%20%20%20return%20(edinburgh_gdf%2C)%0A%0A%0A%40app.cell%0Adef%20_(edinburgh_gdf)%3A%0A%20%20%20%20edinburgh_gdf.explore()%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(edinburgh_gdf%2C%20gpd)%3A%0A%20%20%20%20%23%20the%20boundary%20looks%20a%20bit%20over%20simplified%20and%20is%20missing-out%20bits%20of%20the%20coast.%20so%2C%20bloat%20a%20bit%20to%20recover%20the%20coast%0A%20%20%20%20def%20buffer(gdf%3A%20gpd.GeoDataFrame%2C%20buffer_metres%3A%20float)%20-%3E%20gpd.GeoDataFrame%3A%0A%20%20%20%20%20%20%20%20original_crs%20%3D%20gdf.crs%0A%20%20%20%20%20%20%20%20gdf_buffered%20%3D%20gdf.to_crs(epsg%3D27700).copy()%20%23%20EPSG%2027700%20has%20a%20unit%20of%20metres%20so%20is%20valid%20to%20use%20here%0A%20%20%20%20%20%20%20%20gdf_buffered%5B%22geometry%22%5D%20%3D%20gdf_buffered.geometry.buffer(buffer_metres)%0A%20%20%20%20%20%20%20%20return%20gdf_buffered.to_crs(original_crs)%0A%0A%20%20%20%20edinburgh_buffered_gdf%20%3D%20buffer(edinburgh_gdf%2C%20buffer_metres%3D400)%0A%20%20%20%20return%20(edinburgh_buffered_gdf%2C)%0A%0A%0A%40app.cell%0Adef%20_(edinburgh_buffered_gdf)%3A%0A%20%20%20%20edinburgh_buffered_gdf.explore()%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(Polygon%2C%20gpd%2C%20h3)%3A%0A%20%20%20%20def%20get_h3_cells_for_gdf(gdf%2C%20resolution)%3A%0A%20%20%20%20%20%20%20%20%22%22%22%0A%20%20%20%20%20%20%20%20Get%20all%20H3%20cells%20at%20a%20given%20resolution%20that%20cover%20the%20geometry%20in%20a%20GeoDataFrame.%0A%0A%20%20%20%20%20%20%20%20Parameters%3A%0A%20%20%20%20%20%20%20%20-%20gdf%3A%20GeoDataFrame%20with%20a%20single%20row%0A%20%20%20%20%20%20%20%20-%20resolution%3A%20H3%20resolution%20(0-15)%0A%0A%20%20%20%20%20%20%20%20Returns%3A%0A%20%20%20%20%20%20%20%20-%20GeoDataFrame%20with%20H3%20cell%20polygons%0A%20%20%20%20%20%20%20%20%22%22%22%0A%20%20%20%20%20%20%20%20%23%20Ensure%20we're%20working%20in%20WGS84%20(required%20by%20H3)%0A%20%20%20%20%20%20%20%20gdf_wgs84%20%3D%20gdf.to_crs(epsg%3D4326)%0A%0A%20%20%20%20%20%20%20%20%23%20Get%20the%20geometry%0A%20%20%20%20%20%20%20%20geom%20%3D%20gdf_wgs84.geometry.iloc%5B0%5D%0A%0A%20%20%20%20%20%20%20%20%23%20Get%20H3%20cells%20covering%20the%20polygon%0A%20%20%20%20%20%20%20%20h3_cells%20%3D%20h3.geo_to_cells(geom%2C%20res%3Dresolution)%0A%0A%20%20%20%20%20%20%20%20%23%20Convert%20H3%20cells%20to%20polygons%0A%20%20%20%20%20%20%20%20polygons%20%3D%20%5B%5D%0A%20%20%20%20%20%20%20%20for%20cell%20in%20h3_cells%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20boundary%20%3D%20h3.cell_to_boundary(cell)%0A%20%20%20%20%20%20%20%20%20%20%20%20%23%20h3%20returns%20(lat%2C%20lng)%20tuples%2C%20need%20to%20flip%20to%20(lng%2C%20lat)%20for%20Shapely%0A%20%20%20%20%20%20%20%20%20%20%20%20polygon%20%3D%20Polygon(%5B(lng%2C%20lat)%20for%20lat%2C%20lng%20in%20boundary%5D)%0A%20%20%20%20%20%20%20%20%20%20%20%20polygons.append(polygon)%0A%0A%20%20%20%20%20%20%20%20%23%20Create%20new%20GeoDataFrame%0A%20%20%20%20%20%20%20%20result_gdf%20%3D%20gpd.GeoDataFrame(%0A%20%20%20%20%20%20%20%20%20%20%20%20%7B%22h3_index%22%3A%20list(h3_cells)%7D%2C%20geometry%3Dpolygons%2C%20crs%3D%22EPSG%3A4326%22%0A%20%20%20%20%20%20%20%20)%0A%0A%20%20%20%20%20%20%20%20return%20result_gdf%0A%0A%0A%20%20%20%20%23%20Usage%3A%0A%20%20%20%20%23%20h3_gdf%20%3D%20get_h3_cells_for_gdf(my_gdf%2C%20resolution%3D9)%0A%20%20%20%20return%20(get_h3_cells_for_gdf%2C)%0A%0A%0A%40app.cell%0Adef%20_(edinburgh_buffered_gdf%2C%20get_h3_cells_for_gdf)%3A%0A%20%20%20%20h3_resolution%20%3D%2010%0A%20%20%20%20edinburgh_h3_cells_gdf%20%3D%20get_h3_cells_for_gdf(edinburgh_buffered_gdf%2C%20h3_resolution)%0A%20%20%20%20edinburgh_h3_cells_gdf.head()%0A%20%20%20%20return%20edinburgh_h3_cells_gdf%2C%20h3_resolution%0A%0A%0A%40app.cell%0Adef%20_(edinburgh_h3_cells_gdf)%3A%0A%20%20%20%20edinburgh_h3_cells_gdf.explore()%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(box%2C%20format_grid%2C%20gpd)%3A%0A%20%20%20%20def%20get_bng_quadrants(gdf)%3A%0A%20%20%20%20%20%20%20%20%22%22%22Find%20all%20BNG%20quadrant%20cells%20(e.g.%20NT16NE)%20intersecting%20a%20GeoDataFrame.%22%22%22%0A%20%20%20%20%20%20%20%20gdf%20%3D%20gdf.to_crs(epsg%3D27700)%0A%20%20%20%20%20%20%20%20geom%20%3D%20gdf.union_all()%0A%0A%20%20%20%20%20%20%20%20minx%2C%20miny%2C%20maxx%2C%20maxy%20%3D%20geom.bounds%0A%20%20%20%20%20%20%20%20minx%20%3D%20(minx%20%2F%2F%205000)%20*%205000%0A%20%20%20%20%20%20%20%20miny%20%3D%20(miny%20%2F%2F%205000)%20*%205000%0A%0A%20%20%20%20%20%20%20%20cells%20%3D%20%5B%5D%0A%20%20%20%20%20%20%20%20x%20%3D%20minx%0A%20%20%20%20%20%20%20%20while%20x%20%3C%3D%20maxx%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20y%20%3D%20miny%0A%20%20%20%20%20%20%20%20%20%20%20%20while%20y%20%3C%3D%20maxy%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20cell_geom%20%3D%20box(x%2C%20y%2C%20x%20%2B%205000%2C%20y%20%2B%205000)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20if%20cell_geom.intersects(geom)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20km_ref%20%3D%20format_grid(x%2C%20y%2C%20form%3D%22SS%20E%20N%22).replace(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22%20%22%2C%20%22%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20)%20%20%23%20NT16%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20qy%20%3D%20%22N%22%20if%20(y%20%25%2010000)%20%3E%3D%205000%20else%20%22S%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20qx%20%3D%20%22E%22%20if%20(x%20%25%2010000)%20%3E%3D%205000%20else%20%22W%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20cells.append(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%7B%22grid_ref%22%3A%20km_ref%20%2B%20qy%20%2B%20qx%2C%20%22geometry%22%3A%20cell_geom%7D%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20y%20%2B%3D%205000%0A%20%20%20%20%20%20%20%20%20%20%20%20x%20%2B%3D%205000%0A%0A%20%20%20%20%20%20%20%20return%20gpd.GeoDataFrame(cells%2C%20crs%3D%22EPSG%3A27700%22)%0A%20%20%20%20return%20(get_bng_quadrants%2C)%0A%0A%0A%40app.cell%0Adef%20_(edinburgh_gdf%2C%20get_bng_quadrants)%3A%0A%20%20%20%20edinburgh_bng_quadrants_gdf%20%3D%20get_bng_quadrants(edinburgh_gdf)%0A%20%20%20%20edinburgh_bng_quadrants_gdf.head()%0A%20%20%20%20return%20(edinburgh_bng_quadrants_gdf%2C)%0A%0A%0A%40app.cell%0Adef%20_(edinburgh_bng_quadrants_gdf)%3A%0A%20%20%20%20edinburgh_bng_quadrants_gdf.explore()%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(edinburgh_bng_quadrants_gdf)%3A%0A%20%20%20%20edinburgh_geotiff_urls%20%3D%20%5B%0A%20%20%20%20%20%20%20%20f%22https%3A%2F%2Fsrsp-open-data.s3-eu-west-2.amazonaws.com%2Flidar%2Fphase-5%2Fdtm%2F27700%2Fgridded%2F%7Bref%7D_50CM_DTM_PHASE5.tif%22%0A%20%20%20%20%20%20%20%20for%20ref%20in%20edinburgh_bng_quadrants_gdf%5B%22grid_ref%22%5D%0A%20%20%20%20%5D%0A%20%20%20%20edinburgh_geotiff_urls%0A%20%20%20%20return%20(edinburgh_geotiff_urls%2C)%0A%0A%0A%40app.cell%0Adef%20_(box%2C%20gpd%2C%20mapping%2C%20mask%2C%20np%2C%20pd%2C%20rasterio)%3A%0A%20%20%20%20def%20process_geotiff(%0A%20%20%20%20%20%20%20%20tiff_url%3A%20str%2C%20h3_gdf%3A%20gpd.GeoDataFrame%0A%20%20%20%20)%20-%3E%20gpd.GeoDataFrame%3A%0A%20%20%20%20%20%20%20%20%22%22%22%0A%20%20%20%20%20%20%20%20Process%20a%20single%20GeoTIFF%20and%20compute%20minimum%20height%0A%20%20%20%20%20%20%20%20for%20each%20H3%20cell%20that%20intersects%20with%20the%20raster.%0A%0A%20%20%20%20%20%20%20%20Args%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20tiff_url%3A%20URL%20or%20path%20to%20the%20GeoTIFF%20file%0A%20%20%20%20%20%20%20%20%20%20%20%20h3_gdf%3A%20GeoDataFrame%20with%20H3%20cell%20polygons%2C%20indexed%20by%20h3_index%20(EPSG%3A4326)%0A%0A%20%20%20%20%20%20%20%20Returns%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20GeoDataFrame%20with%20h3_index%20index%2C%20geometry%2C%20and%20columns%3A%20min_height%0A%20%20%20%20%20%20%20%20%20%20%20%20CRS%20is%20EPSG%3A4326%20(H3%20native%20CRS)%0A%20%20%20%20%20%20%20%20%22%22%22%0A%20%20%20%20%20%20%20%20print(f%22Opening%20raster%3A%20%7Btiff_url%7D%22)%0A%0A%20%20%20%20%20%20%20%20with%20rasterio.open(tiff_url)%20as%20src%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20raster_crs%20%3D%20src.crs%0A%20%20%20%20%20%20%20%20%20%20%20%20raster_bounds%20%3D%20src.bounds%0A%20%20%20%20%20%20%20%20%20%20%20%20nodata_value%20%3D%20src.nodata%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20print(f%22%20%20CRS%3A%20%7Braster_crs%7D%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20print(f%22%20%20Bounds%3A%20%7Braster_bounds%7D%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20print(f%22%20%20Nodata%20value%3A%20%7Bnodata_value%7D%22)%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20%23%20Reproject%20H3%20cells%20to%20raster%20CRS%0A%20%20%20%20%20%20%20%20%20%20%20%20print(f%22%20%20Reprojecting%20H3%20cells%20from%20%7Bh3_gdf.crs%7D%20to%20%7Braster_crs%7D%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20h3_reprojected%20%3D%20h3_gdf.to_crs(raster_crs)%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20%23%20Filter%20to%20H3%20cells%20that%20intersect%20the%20raster%20bounds%0A%20%20%20%20%20%20%20%20%20%20%20%20raster_bbox%20%3D%20box(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20raster_bounds.left%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20raster_bounds.bottom%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20raster_bounds.right%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20raster_bounds.top%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20%20%20%20%20intersecting_mask%20%3D%20h3_reprojected.intersects(raster_bbox)%0A%20%20%20%20%20%20%20%20%20%20%20%20h3_subset%20%3D%20h3_reprojected%5Bintersecting_mask%5D%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20print(f%22%20%20Found%20%7Blen(h3_subset)%7D%20H3%20cells%20intersecting%20raster%22)%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20len(h3_subset)%20%3D%3D%200%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20return%20gpd.GeoDataFrame(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20columns%3D%5B%22geometry%22%2C%20%22total_height%22%2C%20%22total_samples%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20crs%3D%22EPSG%3A4326%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20)%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20%23%20Process%20each%20H3%20cell%0A%20%20%20%20%20%20%20%20%20%20%20%20results%20%3D%20%5B%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20for%20h3_index%2C%20row%20in%20h3_subset.iterrows()%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20geom%20%3D%20%5Bmapping(row.geometry)%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20out_image%2C%20out_transform%20%3D%20mask(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20src%2C%20geom%2C%20crop%3DTrue%2C%20all_touched%3DTrue%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20data%20%3D%20out_image%5B0%5D%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%23%20Mask%20out%20nodata%20values%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20if%20nodata_value%20is%20not%20None%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20valid_mask%20%3D%20data%20!%3D%20nodata_value%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20valid_mask%20%3D%20~np.isnan(data)%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20valid_data%20%3D%20data%5Bvalid_mask%5D%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20if%20len(valid_data)%20%3E%200%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20results.append(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22h3_index%22%3A%20h3_index%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22min_height%22%3A%20float(np.min(valid_data))%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%7D%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20)%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20print(f%22%20%20Processed%20%7Blen(results)%7D%20H3%20cells%20with%20valid%20data%22)%0A%0A%20%20%20%20%20%20%20%20if%20not%20results%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20return%20gpd.GeoDataFrame(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20columns%3D%5B%22geometry%22%2C%20%22min_height%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20crs%3D%22EPSG%3A4326%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20)%0A%0A%20%20%20%20%20%20%20%20%23%20Create%20DataFrame%20and%20join%20back%20to%20original%20h3_gdf%20geometry%0A%20%20%20%20%20%20%20%20df%20%3D%20pd.DataFrame(results).set_index(%22h3_index%22)%0A%20%20%20%20%20%20%20%20gdf%20%3D%20h3_gdf.join(df%2C%20how%3D%22inner%22)%0A%0A%20%20%20%20%20%20%20%20return%20gdf%0A%20%20%20%20return%20(process_geotiff%2C)%0A%0A%0A%40app.cell%0Adef%20_(edinburgh_geotiff_urls%2C%20edinburgh_h3_cells_gdf%2C%20process_geotiff)%3A%0A%20%20%20%20edinburgh_height_gdfs%20%3D%20%5B%5D%0A%20%20%20%20for%20i%2C%20url%20in%20enumerate(edinburgh_geotiff_urls)%3A%0A%20%20%20%20%20%20%20%20print(f%22Processing%20%7Bi%20%2B%201%7D%2F%7Blen(edinburgh_geotiff_urls)%7D%3A%20%7Burl%5B%3A80%5D%7D...%22)%0A%20%20%20%20%20%20%20%20gdf%20%3D%20process_geotiff(url%2C%20edinburgh_h3_cells_gdf)%0A%20%20%20%20%20%20%20%20if%20len(gdf)%20%3E%200%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20edinburgh_height_gdfs.append(gdf)%0A%20%20%20%20%20%20%20%20%20%20%20%20print(f%22Found%20data%20for%20%7Blen(gdf)%7D%20H3%20cells%22)%0A%20%20%20%20%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20print(f%22No%20intersecting%20H3%20cells%22)%0A%20%20%20%20return%20(edinburgh_height_gdfs%2C)%0A%0A%0A%40app.cell%0Adef%20_(edinburgh_height_gdfs)%3A%0A%20%20%20%20edinburgh_height_gdfs%5B20%5D.explore(%22min_height%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(List%2C%20gpd%2C%20pd)%3A%0A%20%20%20%20def%20merge_and_compute_min(%0A%20%20%20%20%20%20%20%20gdfs%3A%20List%5Bgpd.GeoDataFrame%5D%2C%0A%20%20%20%20)%20-%3E%20gpd.GeoDataFrame%3A%0A%20%20%20%20%20%20%20%20%22%22%22%0A%20%20%20%20%20%20%20%20Merge%20multiple%20height%20GeoDataFrames%20by%20finding%20min%20of%20min_height%0A%20%20%20%20%20%20%20%20for%20each%20H3%20cell.%0A%0A%20%20%20%20%20%20%20%20Args%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20gdfs%3A%20List%20of%20GeoDataFrames%20with%20h3_index%20index%2C%20geometry%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20and%20columns%3A%20min_height%0A%0A%20%20%20%20%20%20%20%20Returns%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20GeoDataFrame%20with%20h3_index%20index%2C%20geometry%2C%20and%20columns%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20min_height%20(CRS%3A%20EPSG%3A4326)%0A%20%20%20%20%20%20%20%20%22%22%22%0A%20%20%20%20%20%20%20%20%23%20Filter%20out%20empty%20GeoDataFrames%0A%20%20%20%20%20%20%20%20non_empty_gdfs%20%3D%20%5Bgdf%20for%20gdf%20in%20gdfs%20if%20len(gdf)%20%3E%200%5D%0A%0A%20%20%20%20%20%20%20%20print(f%22Merging%20%7Blen(non_empty_gdfs)%7D%20non-empty%20GeoDataFrames%22)%0A%0A%20%20%20%20%20%20%20%20if%20not%20non_empty_gdfs%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20return%20gpd.GeoDataFrame(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20columns%3D%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22geometry%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22min_height%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20crs%3D%22EPSG%3A4326%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20)%0A%0A%20%20%20%20%20%20%20%20%23%20Extract%20geometry%20from%20first%20GeoDataFrame%20for%20each%20h3_index%0A%20%20%20%20%20%20%20%20%23%20(all%20GeoDataFrames%20should%20have%20the%20same%20geometry%20for%20the%20same%20h3_index)%0A%20%20%20%20%20%20%20%20all_geometries%20%3D%20pd.concat(%5Bgdf%5B%5B%22geometry%22%5D%5D%20for%20gdf%20in%20non_empty_gdfs%5D)%0A%20%20%20%20%20%20%20%20geometries%20%3D%20all_geometries.groupby(level%3D0).first()%0A%0A%20%20%20%20%20%20%20%20%23%20Concatenate%20and%20min%20the%20numeric%20columns%0A%20%20%20%20%20%20%20%20all_stats%20%3D%20pd.concat(%0A%20%20%20%20%20%20%20%20%20%20%20%20%5Bgdf%5B%5B%22min_height%22%5D%5D%20for%20gdf%20in%20non_empty_gdfs%5D%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20merged_stats%20%3D%20all_stats.groupby(level%3D0).agg(%0A%20%20%20%20%20%20%20%20%20%20%20%20%7B%22min_height%22%3A%20%22min%22%7D%0A%20%20%20%20%20%20%20%20)%0A%0A%20%20%20%20%20%20%20%20%23%20Combine%20geometry%20and%20stats%20into%20GeoDataFrame%0A%20%20%20%20%20%20%20%20result%20%3D%20gpd.GeoDataFrame(merged_stats.join(geometries)%2C%20crs%3D%22EPSG%3A4326%22)%0A%0A%20%20%20%20%20%20%20%20print(f%22%20%20Merged%20into%20%7Blen(result)%7D%20H3%20cells%22)%0A%20%20%20%20%20%20%20%20print(%0A%20%20%20%20%20%20%20%20%20%20%20%20f%22%20%20Height%20range%3A%20%7Bresult%5B'min_height'%5D.min()%3A.2f%7D%20to%20%7Bresult%5B'min_height'%5D.max()%3A.2f%7D%22%0A%20%20%20%20%20%20%20%20)%0A%0A%20%20%20%20%20%20%20%20return%20result%0A%20%20%20%20return%20(merge_and_compute_min%2C)%0A%0A%0A%40app.cell%0Adef%20_(edinburgh_height_gdfs%2C%20merge_and_compute_min)%3A%0A%20%20%20%20edinburgh_height_gdf%20%3D%20merge_and_compute_min(edinburgh_height_gdfs)%0A%20%20%20%20edinburgh_height_gdf.head()%0A%20%20%20%20return%20(edinburgh_height_gdf%2C)%0A%0A%0A%40app.cell%0Adef%20_()%3A%0A%20%20%20%20%23%20edinburgh_height_gdf.explore(%22avg_height%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(edinburgh_gdf)%3A%0A%20%20%20%20edinburgh_gdf.to_feather(%0A%20%20%20%20%20%20%20%20f%22edinburgh.arrow%22%2C%20compression%3D%22uncompressed%22%0A%20%20%20%20)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(edinburgh_buffered_gdf)%3A%0A%20%20%20%20edinburgh_buffered_gdf.to_feather(%0A%20%20%20%20%20%20%20%20f%22edinburgh_buffered.arrow%22%2C%20compression%3D%22uncompressed%22%0A%20%20%20%20)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(edinburgh_height_gdf%2C%20h3_resolution)%3A%0A%20%20%20%20edinburgh_height_gdf.to_parquet(f%22edinburgh_min_%7Bh3_resolution%7D.geoparquet%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(edinburgh_height_gdf%2C%20h3_resolution)%3A%0A%20%20%20%20edinburgh_height_gdf.to_feather(%0A%20%20%20%20%20%20%20%20f%22edinburgh_min_%7Bh3_resolution%7D.arrow%22%2C%20compression%3D%22uncompressed%22%0A%20%20%20%20)%0A%20%20%20%20return%0A%0A%0Aif%20__name__%20%3D%3D%20%22__main__%22%3A%0A%20%20%20%20app.run()%0A