matlab利用shp图层边缘对栅格进行裁剪代码

Matlab可以使用Mapping Toolbox中的函数对栅格数据进行裁剪。具体代码如下：

% 读取栅格数据
grid = imread('grid.tif');

% 读取shp图层
S = shaperead('boundary.shp');

% 转换shp图层为地理坐标系
boundary = mapshape(S,'XY',true);

% 创建地理参考对象
R = georasterref('RasterSize',size(grid),'LatitudeLimits',[-90 90],'LongitudeLimits',[-180 180]);

% 将栅格数据转换为地理坐标系
grid_geo = geotiff2mat('grid.tif',R);

% 对栅格数据进行裁剪
grid_clip = mapcrop(grid_geo,R,boundary);

% 将裁剪结果保存为栅格数据
geotiffwrite('grid_clip.tif',grid_clip,R);

其中，`grid.tif`为需要裁剪的栅格数据，`boundary.shp`为shp图层，`grid_clip.tif`为裁剪后的栅格数据输出路径。这段代码首先将shp图层转换为地理坐标系，然后将栅格数据也转换为地理坐标系，并使用`mapcrop`函数进行裁剪。最后，将裁剪结果保存为栅格数据。