Building Footprint Extraction for Africa¶

Install package¶

To use the geoai-py package, ensure it is installed in your environment. Uncomment the command below if needed.

In [ ]:

# %pip install geoai-py

# %pip install geoai-py

Import libraries¶

In [ ]:

import geoai

import geoai

Download sample data¶

In [ ]:

raster_url = "https://huggingface.co/datasets/giswqs/geospatial/resolve/main/buildings_africa.tif"

raster_url = "https://huggingface.co/datasets/giswqs/geospatial/resolve/main/buildings_africa.tif"

In [ ]:

raster_path = geoai.download_file(raster_url)

raster_path = geoai.download_file(raster_url)

Initialize the model¶

In [ ]:

extractor = geoai.BuildingFootprintExtractor(model_path="building_footprints_usa.pth")

extractor = geoai.BuildingFootprintExtractor(model_path="building_footprints_usa.pth")

Extract building footprints¶

In [ ]:

masks_path = extractor.generate_masks(
    raster_path,
    output_dir="building_masks.tif",
    min_object_area=1000,
    confidence_threshold=0.5,
    threshold=0.5,
)

masks_path = extractor.generate_masks( raster_path, output_dir="building_masks.tif", min_object_area=1000, confidence_threshold=0.5, threshold=0.5, )

In [ ]:

geoai.view_raster(masks_path, opacity=0.7, colormap="tab20", basemap=raster_url)

geoai.view_raster(masks_path, opacity=0.7, colormap="tab20", basemap=raster_url)

Vectorize masks¶

In [ ]:

gdf = geoai.orthogonalize(
    input_path=masks_path, output_path="building_footprints.geojson", epsilon=1.0
)

gdf = geoai.orthogonalize( input_path=masks_path, output_path="building_footprints.geojson", epsilon=1.0 )

Add geometric attributes¶

In [ ]:

gdf = geoai.add_geometric_properties(gdf)

gdf = geoai.add_geometric_properties(gdf)

Visualize results¶

In [ ]:

geoai.view_vector_interactive(
    gdf, style_kwds={"color": "red", "fillOpacity": 0.2}, tiles=raster_url
)

geoai.view_vector_interactive( gdf, style_kwds={"color": "red", "fillOpacity": 0.2}, tiles=raster_url )