
This course covers agentic coding workflows specializing in geospatial data science. You will get hands-on experience setting up your agent and understand how to use Skills and MCP servers to augment domain-specific capabilities. We cover best practices for using coding agents to rapidly solve complex problems - grounded in established science and human-in-the-loop validation. We cover hands-on examples using Claude Code - including crime hotspot mapping, route optimization, and object detection from aerial imagery. You will also learn how to use agents to manage a personal knowledge base and create a website that automatically updates with new information.
Visit GitHub.com and create a free account. If you already have an account, skip this step.
Follow our Git and GitHub CLI Installation Guide to install and configure Git and GitHub CLI.
We will use conda to install the required Python packages and manage local development environments.
Follow our step-by-step Conda Installation Guide to install Miniconda for your operating system. Once you have a working conda installation, follow the steps below.
conda-forge.Your local development environment is now ready.
We recommend Visual Studio Code (VS Code) editor for this workshop.
Follow our step-by-step Visual Studio Code Installation Guide to install and configure VS Code on your system. Make sure to complete the steps to set VS Code as your default editor.
Follow our step-by-step Claude Code Installation Guide to install and configure Claude Code on your system.
Obsidian is a free app for creating and managing notes. We will be using Obsidian for setting up your personal knowledge base using Claude. Visit Obsidian Download page and install Obsidian for your platform.
OpenRouteService (ORS) provides a free API for routing, distance matrix, geocoding, route optimization etc. using OpenStreetMap data. To obtain an API key, visit the HeiGIT Sign Up page and create an account. Once your account is activated, you can visit the Dashboard and copy the long string displayed under Basic Key.
The code examples and reference materials used in this workshop are
supplied to you in the agentic_coding_geospatial.zip file.
Unzip this file to a directory - preferably to the
<home folder>/Desktop/agentic_coding_geospatial/
folder.
Download agentic_coding_geospatial.zip.
Note: Certification and Support are only available for participants in our paid instructor-led classes.
A good practice is to always setup a CLAUDE.md file for
your project. This file contains instructions and rules to be followed
when working inside the project directory. If you already have a
directory with some code or have cloned a Git repository, Claude Code
can automatically create the CLAUDE.md file using the
/init command. Here we are starting with a blank project,
so we will create a new CLAUDE.md file with
instructions.
<home folder>/Desktop/agentic_coding_geospatial/
where you have unzipped the data package.

CLAUDE.md. Note
that the filename is case sensitive and needs to have the
.md extension.
# CLAUDE.md
## Geospatial Stack Preferences
- Prefer Python-based approaches
- Use existing packages instead of building solutions from scratch
- My preference for geospatial packages are below
- Pandas for tabular data
- GeoPandas for vector geospatial data
- XArray ecosystem (rioxarray, xarray-spatial etc) for raster geospatial data
- Scikit-learn for Machine Learning
- GeoAI (geoai-py) for Deep Learning
- Web Apps
- Streamlit for data driven apps
- Leaflet for interactive mapping apps
- Self contained HTML for small apps
## Other Preferences
- Write code that is simple to understand and explain
- Always install python packages in a conda environment. Never install anything in the base environment. Ask the user for confirmation on their preferred conda environment before installing anything.
- Do not use emojis

cd command to change the
current working directory to the project directory.Windows
Mac/Linux


/help to see all available
commands. Let’s run a few commands. First enable Plan
Mode by entering the /plan command. Next choose the
model using the /model command. The default model choice is
good for now, and you can change the model at any time.
Now we are ready to use Claude Code for our first data analysis task.
We will now use Claude Code to do some data exploration, cleaning and analysis. This will give you a feel of how agents work and how to guide them. We will take the Police.uk data on street-level crime for the City of London Police and identify crime hotspots.
I would like to do an exploratory data analysis of crime data in the folder london_crime_2024.
* Give me a summary of the rows, columns and any missing data.
* Flag any data quality issues.
* Create charts showing different types of crime and monthly patterns.

CLAUSE.md states, all the work
needs to be done in a conda virtual environment. Claude will prompt you
to choose an environment. Choose the
claude_code_workshop environment created for this class
and press Enter.

Depending on your operating system, Claude Code version, and selected model, you may see a slightly different set of questions.




output/ subfolder.
Now that we have a good understanding of the data structure, we can ask Claude to help us with analyzing the data. In this section, we will go through the workflow of identifying hotspots off thefts. Starting with data cleaning, we will iteratively improve our analysis and learn how we can use existing references (research papers, manuals etc.) to ground the behavior of the model.
I want to analyze thefts data. Suggest which categories should be merged.




Filter to data to only for the 7 theft-family categories and
write a single merged CSV with those records.
Remove missing data and unnecessary fields.


I want to map the theft hotspots.
Suggest what would be the best technique for the analysis.


Create a static visualiztion using matplotlib using gaussian_kde.



I have some reference material for crime mapping in documents/hotspots.pdf.
Incorporate suggestions from that and give me an updated plan.



Save the updated script and chart with the _v2 suffix.



Use the criteria "Greater than 5x mean" to extract hotspots and convert then to polygons.
Save the results as a GeoJSON london_theft_hotspots.geojson.



It’s a good practice to keep an eye on your token usage and costs. You can use the following command to know the utilization in the current session and how much quota is left in your plan.
/usage

In this section, we will learn how to use Claude Code to build Jupyter notebooks for data analysis and interactively guide the agent to reach the desired solution and validate the results. We will be implementing a real-world scenario of optimizing delivery routes based on a route optimization algorithm.
In your data package, you will find the following files in the
data/route_optimization/ folder.
Delivery_Locations.geojson: A set of 100 addresses in
Washington DC area representing delivery locations.Grocery_Store_Locations.geojson: A set of 14 grocery
store locations which needs to fulfill the delivery to these
locations.We will use Claude Code to build a Jupyter Notebook that builds an end-to-end workflow that loads the data, finds optimal delivery routes and writes PDF manifests with delivery schedule for each store.
Let’s install the Claude Code for VS Code extension.



Now that we have Claude Code running inside VS Code, we can prompt it to build a Jupyter notebook for us.
agentic_coding_geospatial is
open in VS Code. Enter the following prompt.I want to build a notebook to find optimal delivery routes.
Create a new notebook scripts/route_optimization.ipynb
Read the .geojson files in data/route_optimization/ folder and
display the data on an interactive map.

CLAUDE.md already has preferences listed for how to setup
the environment and what packages to use. Claude Code will start
creating the notebook and will pause periodically for confirmation. You
can choose Yes to proceed.
route_optimization.ipynb will be created in the
scripts/ folder as instructed.

Now, let’s build the core of our data analysis to find the optimal delivery routes. This is a complex problem that requires a lot more context and direction. We will use the Plan mode and guide Claude Code towards building the code that solves the problem.
route_optimization.ipynb is
open. Open Claude Code from the Toggle Chat button.

I want to use OpenRouteService API.
Give me options on different approaches to plan for the most optimal delivery routes from each store.
Here are the constrants:
- A single trip will consist of 10 deliveries at maximum.
- The delivery will originate at a grocery store and end at the same delivery store.
- All stores are capable of fulfilling the order but can do a maximum of 5 delivery trips each day.
I want to know the most optimal route to minimize the distance traveled along with a schedule of the trips.


<your-api-key> with your actual API key. If you do
have your key handy, you can login to your HeiGIT Account and copy
the long string displayed under Basic Key.Yes, I have a key. Save the API key in a separate cell in the notebook.
ORS_API_KEY='<your-api-key>'

claude_code_workshop conda environment that we have setup.
You can then click Submit answers.



openrouteservice Python package. We can instead ask it to
use the API directly instead of a client library by adding a comment
like below.Instead, use the REST API directly.






When building with AI - you need to spend most of your time validating the results and guiding the model. Let’s learn some techniques for validation.
I see that some stores are not assigned any deliveries even if they are quite close to an address.
The assignment tries to assign 10 deliveries in each route -
making some of them quite long while nearby stores end up without any deliveries.
Update the plan to address this problem.



Save the routes as a single GeoJSON file in the outputs/route_optimization/
folder so I can validate the results.

routes.geojson will be created. Locate it in
the VS Code Explorer. Right-click and select Reveal in File
Explorer to open the folder on your computer.
routes.geojson, along with the source data
Delivery_Locations.geojson and
Grocery_Store_Locations.geojson and verify the
results.
Generate PDF manifests in the outputs/route_optimization/ folder containing
delivery order and addresses for each trip.
1 manifest per store.



Now you are ready to put your agentic coding skills to test by solving a spatial analysis problem from scratch.
Problem Statement
You are a data scientist for a real-estate company who wants to enhance their property listings by assigning a livability score to each of their properties to help buyers pick properties that suit their lifestyle.
Your task is to calculate such a score for all the Apartment Buildings in your city. Consider the factors that are important to home buyers. You can design the score using metrics such as:
These are some examples and you can pick the criteria that make sense to you based on the knowledge of your city. You can use OpenStreetMap to obtain the data for apartment building polygons and the required amenities.
Required output
Tips
osmnx which can easily extract the required data from
OpenStreetMap.network distance/isochrone
instead of simple circular buffer for a more realistic distances.Skills are a set of instructions that teach an AI agent how to
perform specific, multi-step workflows. Skills in Claude
Code are defined by a SKILL.md containing
natural-language instructions for the agent along with specific examples
and patterns. You can also publish your skill on GitHub and make it
available for others to download and use with their own agent.
We will install the Humanizer skill that improves the text generated by AI models and makes them sound more natural.
cd command
to change the current working directory to the project directory. Once
you are inside the project folder, launch Claude Code by entering the
following command.
npx skills
command provides an easy to to install skills for a wide range of coding
agents. However this requires installation of Node.js on your system.
For Claude Code, you can install Skills as plugins. Run the following
command to add the GitHub repository.
humanizer
skills as a plugin.


Once installed, Skills are triggered automatically when your prompt matches the task handled by the skill. We will use the Humanizer skills to review and rewrite a paragraph.
Give me a 100-word summary of this article
https://cloudnativegeo.org/blog/2026/02/the-technical-debt-of-earth-embedding-products/

humanizer skill by asking Claude
Code to review and update the text. Enter the following prompt and press
Enter.Review the summary and remove signs of AI writing

humanizer skill will be used to rewrite the text by
removing common AI-writing patterns.
The real power of Skills lies in automating custom workflows. We will
learn how to create a new skill for your custom workflow and use it. We
will create a new skill named create-cog for converting
raster data into a Cloud Optimized GeoTIFF (COG) format.
Create a new skill from the following text: and paste the
content below.---
name: create-cog
description: Convert a raster to a Cloud-Optimized GeoTIFF (COG)
---
Use this skill when you need to convert any geospatial raster data to a
Cloud-Optimized GeoTIFF.
## Workflow
1. Activate the conda environment that has GDAL. Never use the base environment.
2. Check the GDAL version, since the conversion command differs by version.
3. Check if the input is already a valid COG. If it is, stop and report
that. Do not re-convert.
4. Convert, using the compression settings and naming convention below.
5. Validate the output and report the result.
## Required Tools
* `gdal` Conda package. Install using `conda install -c conda-forge gdal`
* `rio-cogeo` Python package for validation. Install using `pip install
rio-cogeo`
Both must come from a conda environment, not base. Check what is available
before installing:
conda env list
conda activate <env>
gdal --version
If neither is installed, ask the user which environment to install into
before proceeding.
## Validating COGs
rio cogeo validate <file>
Always validate the input before converting and the output after converting.
## Converting Raster to a COG
Check `gdal --version` first, then use the matching command.
GDAL version >= 3.11 (the `gdal` subcommand interface)
gdal raster convert -f COG <input file> <output file>
GDAL version < 3.11 (the legacy utility)
gdal_translate -of COG <input file> <output file>
Note the flag difference. The new `gdal raster convert` takes `-f` (or
`--format`). It does NOT accept `-of`, which is the legacy `gdal_translate`
flag, and fails with "Option 'o' is not a boolean option."
Creation options are passed with `--co KEY=VALUE` on the new interface,
`-co KEY=VALUE` on the legacy one.
If the output file already exists, the command errors out. Add
`--overwrite` only when you intend to replace it.
## Compression
The COG driver already applies LZW compression by default, so the output is
compressed even if you pass nothing. But the default uses no predictor,
which wastes a lot of space on continuous data. Set the compression and
predictor explicitly. Use DEFLATE as the default, since it is readable by
every GeoTIFF client.
gdal raster convert -f COG --co COMPRESS=DEFLATE --co PREDICTOR=YES <input> <output>
Only exception - If the input data is JPEG compressed - typically used in
aerial/drone imagery), retain the same with COMPRESS=JPEG
## File Naming Convention
Append the text `_cog` to the converted filename. If the input is
`data.tif`, the output should be named `data_cog.tif`
## Remote Files
If the given file is a remote file, do not download it. Use GDAL Virtual
File Systems by prefixing the URL with `/vsicurl/`
rio cogeo validate /vsicurl/https://example.com/data.tif
gdal raster convert -f COG /vsicurl/https://example.com/data.tif data_cog.tif
A remote input gives no obvious output location. Write to the current
working directory unless a location is specified by the user.
In Windows Powershell, set the following environment variable `MSYS_NO_PATHCONV=1`
to avoid errors when using paths starting with `/vsicurl`

.claude directory in the current
folder.
/reload-skills

Convert the file data/chirps/chirps-v3.0.2025.tif to a Cloud Optimized GeoTIFF


Convert the file https://data.chc.ucsb.edu/products/CHIRPS/v3.0/annual/global/tifs/chirps-v3.0.2024.tif to a COG

/vsicurl virtual file system to access the file instead of
downloading it.
/create-cog Convert this file https://storage.googleapis.com/spatialthoughts-public-data/ntl/viirs/viirs_ntl_2021_global.tif


skills CLI.Coming soon
Coming soon
Coming soon
MCP (Model Context Protocol) is an open-source standard for connecting agents to external systems. MCP servers are similar to APIs - but are designed specifically for use by AI Agents. MCP servers allow your agents to access specialized data sources and applications. We will setup and configure a connection to the Google Colab MCP server which will allow Claude Code to create, modify and run Google Colab notebooks.
cd command
to change the current working directory to the project directory.
Setup an MCP server according to the instructions at
https://github.com/googlecolab/colab-mcp


.claude.json with the required
information. Restart Claude Code for the changes to take effect. You can
exit the session using Ctrl + C and start it again by
typing claude.
/mcp

colab-mcp
server connected.
Now that we have configured the MCP server, we can prompt Claude Code to create cloud-hosted notebooks using Google Colab.
Create a new colab notebook

colab-mcp server to open a
connection. Approve the execution requests.


Add a cell to print "hello world" and run it


Colab MCP server can be a bit flaky and you may see errors when establishing the connection. Here are a few things you can try.
/mcp and check whether the colab-mcp
server is connected. If it is disconnected, select it and press
Enter to reconnect.


Create a new Colab Notebook named "building_detection".
This notebook will use a pre-trained model for building footprint extraction from GeoAI and run it on the provided image.
Structure the notebook following this example https://opengeoai.org/examples/building_footprints_africa/
The provided image can be large. So use the GeoAI package to split the input image into tiles
Reproject the tiles to a suitable local UTM crs.
The model is trained on images around 50cm-1m resolution images. Check the resolution of the provided image and if it is higher resolution (i.e. 5cm, 10cm) - resample the tiles to 50cm resolution before running the inference.
Save the results as a GeoJSON file in EPSG:4326 CRS.
Use this image from OpenAerialMap to test the notebook https://oin-hotosm-temp.s3.us-east-1.amazonaws.com/69493c8084a859b011c94266/0/69493c8084a859b011c94267.tif




Run the notebook



building_footprint.geojson file.


https://oin-hotosm-temp.s3.us-east-1.amazonaws.com/69493c8084a859b011c94266/0/69493c8084a859b011c94267.tif




If you want to catch up to this step and use a sample notebook generated in this step, click the button below.
You will notice that since the footprints were vectorized from a raster layer, they have jagged edges. We can make the output polygons much better by Regularizing the polygons. GeoAI has several algorithms for Building Regularization. Prompt Claude Code to add a step at the end to regularize the polygons and save them.

This course material is licensed under a Creative Commons Attribution 4.0 International (CC BY 4.0). You are free to re-use and adapt the material but are required to give appropriate credit to the original author as below:
Agentic Coding for Geospatial course by Ujaval Gandhi www.spatialthoughts.com
© 2026 Spatial Thoughts www.spatialthoughts.com
If you want to report any issues with this page, please comment below.