matlab根据shp文件裁剪tif数据代码

在MATLAB中，你可以使用`geotiffread`函数读取TIFF图像数据，然后利用`shaperead`和`imcrop`函数配合地理空间信息从Shapefile（`.shp`）裁剪数据。这里是一个简单的步骤示例：

% 加载Shapefile数据
shapes = shaperead('your_shapefile.shp'); % 替换为你的shapefile路径

% 获取Shapefile中的几何信息（例如，边界）
exteriorPoly = shapes(1).Part; % 假设第一个形状是外边界
bbox = bbox(shapes); % 获取所有形状的边界盒

% 选择需要裁剪的TIFF文件
tifData = geotiffread('your_tif_file.tif'); % 替换为你的TIFF文件路径

% 转换边界框到相同的坐标参考系统（如果你的TIFF和Shapefile使用不同的CRS）
crs TIFF.GeoCoordTransform(tifData.GeoCoordTransform, bbox.CRS);

% 使用边界框裁剪TIFF数据
croppedData = imcrop(tifData, [bbox.Xmin, bbox.Ymin, bbox.Width, bbox.Height], 'OutputView', crs);

% 将结果保存为新的TIFF文件
geotiffwrite('output_cropped.tif', croppedData, tifData.GeoTransform);

注意，这个例子假设你的Shapefile和TIFF文件都有地理元数据，并且可以相互转换坐标系。如果数据格式或位置有误，你需要相应地调整代码。