Earth Science Data Queries¶

This notebook demonstrates GeoAgent's expanded data query capabilities across many different Earth observation collections — fire detection, surface water, snow cover, temperature, vegetation indices, and more.

These examples mirror the types of queries supported by Microsoft Earth Copilot.

In [ ]:

# %pip install geoagent

# %pip install geoagent

In [ ]:

from geoagent import GeoAgent

agent = GeoAgent()

from geoagent import GeoAgent agent = GeoAgent()

1. Fire Detection / Thermal Anomalies¶

Search for MODIS thermal anomaly data to detect active fires and burned areas.

In [ ]:

result = agent.chat("Show me thermal anomalies in Australia in January 2024")
print(
    f"Success: {result.success} | Items: {result.data.total_items if result.data else 0}"
)
if result.data and result.data.items:
    print(f"Collection: {result.data.items[0].get('collection')}")
result.map

result = agent.chat("Show me thermal anomalies in Australia in January 2024") print( f"Success: {result.success} | Items: {result.data.total_items if result.data else 0}" ) if result.data and result.data.items: print(f"Collection: {result.data.items[0].get('collection')}") result.map

In [ ]:

result = agent.chat("Show me burned area mapping for Montana wildfire regions in 2023")
print(
    f"Success: {result.success} | Items: {result.data.total_items if result.data else 0}"
)
result.map

result = agent.chat("Show me burned area mapping for Montana wildfire regions in 2023") print( f"Success: {result.success} | Items: {result.data.total_items if result.data else 0}" ) result.map

In [ ]:

print(result.code)

print(result.code)

2. Surface Water¶

Visualize surface water extent and changes using the JRC Global Surface Water dataset.

In [ ]:

result = agent.chat("Show me surface water in Bangladesh")
print(
    f"Success: {result.success} | Items: {result.data.total_items if result.data else 0}"
)
if result.data and result.data.items:
    print(f"Collection: {result.data.items[0].get('collection')}")
result.map

result = agent.chat("Show me surface water in Bangladesh") print( f"Success: {result.success} | Items: {result.data.total_items if result.data else 0}" ) if result.data and result.data.items: print(f"Collection: {result.data.items[0].get('collection')}") result.map

In [ ]:

print(result.code)

print(result.code)

3. Snow Cover¶

Visualize snow cover extent using MODIS snow cover data.

In [ ]:

result = agent.chat("Show me snow cover in Quebec in January 2024")
print(
    f"Success: {result.success} | Items: {result.data.total_items if result.data else 0}"
)
if result.data and result.data.items:
    print(f"Collection: {result.data.items[0].get('collection')}")
result.map

result = agent.chat("Show me snow cover in Quebec in January 2024") print( f"Success: {result.success} | Items: {result.data.total_items if result.data else 0}" ) if result.data and result.data.items: print(f"Collection: {result.data.items[0].get('collection')}") result.map

In [ ]:

print(result.code)

print(result.code)

4. Surface Temperature¶

Visualize land surface temperature (LST) and sea surface temperature (SST) from MODIS.

In [ ]:

result = agent.chat("Show me sea surface temperature near Madagascar")
print(
    f"Success: {result.success} | Items: {result.data.total_items if result.data else 0}"
)
if result.data and result.data.items:
    print(f"Collection: {result.data.items[0].get('collection')}")
result.map

result = agent.chat("Show me sea surface temperature near Madagascar") print( f"Success: {result.success} | Items: {result.data.total_items if result.data else 0}" ) if result.data and result.data.items: print(f"Collection: {result.data.items[0].get('collection')}") result.map

In [ ]:

result = agent.chat("Show me land surface temperature for Phoenix Arizona in July 2024")
print(
    f"Success: {result.success} | Items: {result.data.total_items if result.data else 0}"
)
result.map

result = agent.chat("Show me land surface temperature for Phoenix Arizona in July 2024") print( f"Success: {result.success} | Items: {result.data.total_items if result.data else 0}" ) result.map

In [ ]:

print(result.code)

print(result.code)

5. Vegetation Indices (MODIS)¶

Visualize MODIS vegetation index products (NDVI/EVI at 250m resolution).

In [ ]:

result = agent.chat("Show me vegetation indices for Ukraine in June 2024")
print(
    f"Success: {result.success} | Items: {result.data.total_items if result.data else 0}"
)
if result.data and result.data.items:
    print(f"Collection: {result.data.items[0].get('collection')}")
result.map

result = agent.chat("Show me vegetation indices for Ukraine in June 2024") print( f"Success: {result.success} | Items: {result.data.total_items if result.data else 0}" ) if result.data and result.data.items: print(f"Collection: {result.data.items[0].get('collection')}") result.map

In [ ]:

print(result.code)

print(result.code)

6. Cropland¶

Visualize USDA Cropland Data Layer (CDL) for agricultural mapping.

In [ ]:

result = agent.chat("Show me cropland data for Florida")
print(
    f"Success: {result.success} | Items: {result.data.total_items if result.data else 0}"
)
if result.data and result.data.items:
    print(f"Collection: {result.data.items[0].get('collection')}")
result.map

result = agent.chat("Show me cropland data for Florida") print( f"Success: {result.success} | Items: {result.data.total_items if result.data else 0}" ) if result.data and result.data.items: print(f"Collection: {result.data.items[0].get('collection')}") result.map

In [ ]:

print(result.code)

print(result.code)

7. Radar / SAR Imagery¶

Search for Sentinel-1 SAR (Synthetic Aperture Radar) imagery — useful for all-weather, day/night monitoring.

In [ ]:

result = agent.chat("Show me radar imagery of Houston Texas during August 2017")
print(
    f"Success: {result.success} | Items: {result.data.total_items if result.data else 0}"
)
if result.data and result.data.items:
    print(f"Collection: {result.data.items[0].get('collection')}")
result.map

result = agent.chat("Show me radar imagery of Houston Texas during August 2017") print( f"Success: {result.success} | Items: {result.data.total_items if result.data else 0}" ) if result.data and result.data.items: print(f"Collection: {result.data.items[0].get('collection')}") result.map

In [ ]:

print(result.code)

print(result.code)

8. LiDAR Height Data¶

Visualize 3DEP LiDAR-derived height above ground data.

In [ ]:

result = agent.chat("Show me LIDAR height data for Denver Colorado")
print(
    f"Success: {result.success} | Items: {result.data.total_items if result.data else 0}"
)
if result.data and result.data.items:
    print(f"Collection: {result.data.items[0].get('collection')}")
result.map

result = agent.chat("Show me LIDAR height data for Denver Colorado") print( f"Success: {result.success} | Items: {result.data.total_items if result.data else 0}" ) if result.data and result.data.items: print(f"Collection: {result.data.items[0].get('collection')}") result.map

In [ ]:

print(result.code)

print(result.code)

9. HLS (Harmonized Landsat Sentinel)¶

Search for Harmonized Landsat Sentinel (HLS) imagery — combined Landsat and Sentinel-2 data.

In [ ]:

result = agent.chat("Show me HLS Landsat imagery for Greece in June 2024")
print(
    f"Success: {result.success} | Items: {result.data.total_items if result.data else 0}"
)
if result.data and result.data.items:
    print(f"Collection: {result.data.items[0].get('collection')}")
result.map

result = agent.chat("Show me HLS Landsat imagery for Greece in June 2024") print( f"Success: {result.success} | Items: {result.data.total_items if result.data else 0}" ) if result.data and result.data.items: print(f"Collection: {result.data.items[0].get('collection')}") result.map

In [ ]:

print(result.code)

print(result.code)

10. Elevation (Multiple DEMs)¶

GeoAgent can search different DEM sources: Copernicus 30m, NASADEM, ALOS World.

In [ ]:

result = agent.chat("Show me elevation profile for the Grand Canyon")
print(
    f"Success: {result.success} | Items: {result.data.total_items if result.data else 0}"
)
result.map

result = agent.chat("Show me elevation profile for the Grand Canyon") print( f"Success: {result.success} | Items: {result.data.total_items if result.data else 0}" ) result.map

In [ ]:

result = agent.chat("Show me ALOS World DEM for Berlanga Spain")
print(
    f"Success: {result.success} | Items: {result.data.total_items if result.data else 0}"
)
result.map

result = agent.chat("Show me ALOS World DEM for Berlanga Spain") print( f"Success: {result.success} | Items: {result.data.total_items if result.data else 0}" ) result.map

11. Net Ecosystem Production¶

Visualize MODIS net primary production / gross primary production data.

In [ ]:

result = agent.chat("Show me net production for San Jose California")
print(
    f"Success: {result.success} | Items: {result.data.total_items if result.data else 0}"
)
if result.data and result.data.items:
    print(f"Collection: {result.data.items[0].get('collection')}")
result.map

result = agent.chat("Show me net production for San Jose California") print( f"Success: {result.success} | Items: {result.data.total_items if result.data else 0}" ) if result.data and result.data.items: print(f"Collection: {result.data.items[0].get('collection')}") result.map

In [ ]:

print(result.code)

print(result.code)

12. NADIR BRDF Reflectance¶

Visualize MODIS BRDF-adjusted nadir reflectance.

In [ ]:

result = agent.chat("Show me nadir BRDF reflectance for Mexico")
print(
    f"Success: {result.success} | Items: {result.data.total_items if result.data else 0}"
)
if result.data and result.data.items:
    print(f"Collection: {result.data.items[0].get('collection')}")
result.map

result = agent.chat("Show me nadir BRDF reflectance for Mexico") print( f"Success: {result.success} | Items: {result.data.total_items if result.data else 0}" ) if result.data and result.data.items: print(f"Collection: {result.data.items[0].get('collection')}") result.map

In [ ]:

print(result.code)

print(result.code)

---

Summary¶

Category	Collection	Example Query
Fire / Thermal	modis-14A1-061	"Show thermal anomalies in Australia"
Surface Water	jrc-gsw	"Show surface water in Bangladesh"
Snow Cover	modis-10A1-061	"Show snow cover in Quebec"
Surface Temperature	modis-11A1-061	"Show sea surface temperature near Madagascar"
Vegetation Indices	modis-13Q1-061	"Show vegetation indices for Ukraine"
Cropland	usda-cdl	"Show cropland data for Florida"
Radar / SAR	sentinel-1-grd	"Show radar imagery of Houston"
LiDAR Height	3dep-lidar-hag	"Show LIDAR height for Denver"
HLS	hls-l30	"Show HLS imagery for Greece"
Elevation	cop-dem-glo-30	"Show elevation for Grand Canyon"
Net Production	modis-17A2H-061	"Show net production for San Jose"
BRDF Reflectance	modis-43A4-061	"Show nadir BRDF for Mexico"