reclassify_raster
Reclassify raster values by mapping old values to new values using a dictionary, then save the reclassified raster to a specified output path.
Instructions
Reclassify raster values using a mapping dictionary. Args: raster_path: Path to the input raster. reclass_map: Dictionary mapping old values to new values (e.g., {1: 10, 2: 20}). output_path: Path to save the reclassified raster. Returns: Dictionary with status and message.
Input Schema
TableJSON Schema
| Name | Required | Description | Default |
|---|---|---|---|
| raster_path | Yes | ||
| reclass_map | Yes | ||
| output_path | Yes |
Output Schema
TableJSON Schema
| Name | Required | Description | Default |
|---|---|---|---|
No arguments | |||
Implementation Reference
- src/gis_mcp/rasterio_functions.py:87-117 (handler)The core handler function for the 'reclassify_raster' tool. It reads the input raster, applies the reclassification mapping using numpy, and writes the output raster using rasterio.
def reclassify_raster(raster_path: str, reclass_map: dict, output_path: str) -> Dict[str, Any]: """ Reclassify raster values using a mapping dictionary. Args: raster_path: Path to the input raster. reclass_map: Dictionary mapping old values to new values (e.g., {1: 10, 2: 20}). output_path: Path to save the reclassified raster. Returns: Dictionary with status and message. """ try: import rasterio import numpy as np with rasterio.open(raster_path) as src: data = src.read(1) profile = src.profile.copy() reclass_data = np.copy(data) for old, new in reclass_map.items(): reclass_data[data == old] = new output_path_resolved = resolve_path(output_path, relative_to_storage=True) output_path_resolved.parent.mkdir(parents=True, exist_ok=True) with rasterio.open(str(output_path_resolved), "w", **profile) as dst: dst.write(reclass_data, 1) return { "status": "success", "message": f"Raster reclassified and saved to '{output_path_resolved}'.", "output_path": str(output_path_resolved) } except Exception as e: logger.error(f"Error in reclassify_raster: {str(e)}") return {"status": "error", "message": str(e)} - src/gis_mcp/main.py:66-72 (registration)Imports the rasterio_functions module (line 69), which triggers the @gis_mcp.tool() decorators to register the reclassify_raster tool with the MCP server.
from . import ( geopandas_functions, shapely_functions, rasterio_functions, pyproj_functions, pysal_functions, ) - src/gis_mcp/rasterio_functions.py:11-35 (registration)Resource function listing 'reclassify_raster' among available rasterio operations, aiding tool discovery.
@gis_mcp.resource("gis://operation/rasterio") def get_rasterio_operations() -> Dict[str, List[str]]: """List available rasterio operations.""" return { "operations": [ "metadata_raster", "get_raster_crs", "clip_raster_with_shapefile", "resample_raster", "reproject_raster", "weighted_band_sum", "concat_bands", "raster_algebra", "compute_ndvi", "raster_histogram", "tile_raster", "raster_band_statistics", "extract_band", "zonal_statistics", "reclassify_raster", "focal_statistics", "hillshade", "write_raster" ] } - src/gis_mcp/__init__.py:7-7 (registration)Wildcard import of rasterio_functions.py in __init__.py, ensuring tools are registered when the package is imported.
from .rasterio_functions import * - Function signature and docstring define the input schema (raster_path: str, reclass_map: dict, output_path: str) and output format.
def reclassify_raster(raster_path: str, reclass_map: dict, output_path: str) -> Dict[str, Any]: """ Reclassify raster values using a mapping dictionary. Args: raster_path: Path to the input raster. reclass_map: Dictionary mapping old values to new values (e.g., {1: 10, 2: 20}). output_path: Path to save the reclassified raster. Returns: Dictionary with status and message. """