rasterio

Read and write georeferenced raster data with a Pythonic GDAL interface.

19|2|Updated Jan 31, 2026
One-click install
npx skills add https://github.com/tondevrel/scientific-agent-skills --skill rasterio
Or copy as Structured Prompt for Agent
Please help me install this Agent Skill.
Skill: rasterio
Source: https://github.com/tondevrel/scientific-agent-skills/tree/main/skills/rasterio
Command: npx skills add https://github.com/tondevrel/scientific-agent-skills --skill rasterio

SYSTEM DOCUMENTATION & REQUIREMENTS

💡 This Skill requires numpy, matplotlib, geopandas, shapely, fiona, and includes scripts (resource) and references (resource) components.

What problem does it solve?

This Skill simplifies the complex task of reading, writing, and manipulating georeferenced raster data, such as satellite imagery and elevation models, making geospatial analysis accessible in Python.

Core Features & Use Cases

  • Read/Write Georeferenced Rasters: Handles formats like GeoTIFF, NetCDF, etc.
  • Coordinate Reference System (CRS) Management: Works with CRS and performs reprojections.
  • Band Math & Spectral Indices: Computes indices like NDVI, NDWI for land analysis.
  • Vector Clipping: Clips raster data to polygon boundaries.
  • Memory-Efficient I/O: Processes large files without loading them entirely into memory.
  • Use Case: Analyze changes in vegetation health over time by calculating NDVI from satellite imagery for specific regions, clipping the results to administrative boundaries.

Quick Start

Use the rasterio skill to read band 1 from 'satellite_image.tif' and display it.

Frequently Asked Questions about rasterio

High-intent search queries and answers about installing and using this skill.

FAQPage Schema
How do I read and write georeferenced raster data like GeoTIFF in Python?

To read and write georeferenced raster data, use a Pythonic interface to GDAL that handles formats like GeoTIFF and NetCDF. It simplifies complex geospatial file I/O for satellite imagery and elevation models.

What's the best way to compute spectral indices like NDVI from satellite imagery?

Computing spectral indices like NDVI from satellite imagery is done through band math. This calculates indices using multi-band raster data to analyze land characteristics and vegetation health.

How do I clip raster data to vector polygon boundaries for a specific region?

Clipping raster data to vector polygon boundaries involves using vector geometries to mask the raster. This extracts and limits satellite imagery or elevation models to specific administrative regions.

Can I process large raster files without loading them entirely into memory?

Yes, you can process large raster files without loading them entirely into memory using windowed I/O. This memory-efficient approach reads specific portions of georeferenced rasters as needed for analysis.

How do I handle coordinate reference systems and reproject geospatial rasters?

Handling coordinate reference systems and reprojecting geospatial rasters is supported natively. You can manage CRS transformations to align satellite imagery and elevation models from different sources.

Do I need numpy and geopandas to perform geospatial raster analysis?

Yes, numpy and geopandas are required dependencies for performing geospatial raster analysis. Numpy supports band math array operations, while geopandas handles the vector data needed for clipping.