Geospatial Data Essentials

import geoai
import rasterio

raster_url = "https://data.source.coop/opengeos/geoai/naip-train.tif"
raster_path = geoai.download_file(raster_url)

with rasterio.open(raster_path) as src:
    print(f"Shape: {src.height} x {src.width}")
    print(f"Bands: {src.count}")
    print(f"CRS: {src.crs}")
    print(f"Resolution: {src.res}")
    band1 = src.read(1)  # Read the first band as a NumPy array
    print(f"Band 1 dtype: {band1.dtype}")

geoai.print_raster_info(raster_path, figsize=(18, 10))

clip_raster_path = "naip_clip.tif"
geoai.clip_raster_by_bbox(
    raster_path,
    clip_raster_path,
    bbox=(0, 0, 500, 500),
    bands=[1, 2, 3],
    bbox_type="pixel",
)
geoai.view_image(clip_raster_path)

import geopandas as gpd

vector_url = "https://data.source.coop/opengeos/geoai/naip-train-buildings.geojson"
vector_path = geoai.download_file(vector_url)
gdf = gpd.read_file(vector_path)
print(f"Features: {len(gdf)}")
print(f"CRS: {gdf.crs}")
gdf.head()

geoai.print_vector_info(vector_path, figsize=(18, 10))

geoai.view_vector(vector_path, basemap=True)

geoai.view_vector(
    vector_path,
    raster_path=raster_path,
    outline_only=True,
    edge_color="red",
    figsize=(18, 10),
)

gdf_utm = gdf.to_crs(crs="EPSG:26911")  # Reproject to UTM Zone 11N
print(f"Original CRS: {gdf.crs}")
print(f"New CRS: {gdf_utm.crs}")
print(f"Sample coordinates before: {gdf.geometry.iloc[0].centroid}")
print(f"Sample coordinates after:  {gdf_utm.geometry.iloc[0].centroid}")

import leafmap

leafmap.reproject(
    image=raster_path,
    output="naip_reprojected.tif",
    dst_crs="EPSG:3857",
    resampling="bilinear",
)

leafmap.split_raster(
    raster_path,
    out_dir="naip_tiles",
    tile_size=512,
    overlap=128,
)