MATLAB shp裁剪栅格

MATLAB是一种强大的数学软件，它也支持地理信息系统（GIS）功能。如果你想要在MATLAB中裁剪一个栅格数据（通常是.shp格式的Shapefile），可以使用`shaperead`函数读取.shp文件，然后结合`imcrop`或`geotiffread`函数来裁剪栅格数据。以下是基本步骤：

1. **读取.shp数据**：

   shpData = shaperead('your_shapefile.shp');

2. **获取栅格数据**：
如果你的栅格数据是以GeoTIFF或其他地理格式存在，例如`rasterData.tif`：

   raster = geotiffread('rasterData.tif');

3. **裁剪栅格基于.shp边界**：

   % 获取Shapefile的几何信息
   geometry = poly2mask(shpData.X, shpData.Y, size(raster, 1), size(raster, 2));
   
   % 裁剪栅格
   croppedRaster = raster .* geometry;

4. **保存结果**：
要保存裁剪后的栅格数据，可以使用`imwrite`或`geotiffwrite`：

   imwrite(croppedRaster, 'cropped_raster_data.tif', 'GeoKeyDirectory', 'rasterData');

注意：这里假设你的栅格数据和Shapefile有相同的投影坐标系统。如果需要转换坐标，你需要先进行相应的坐标变换。