matalb基于shp文件对栅格文件进行裁剪并保存

要基于shapefile文件对栅格文件进行裁剪并保存，您可以使用以下代码：

% 加载shapefile文件
shp = shaperead('your_shapefile.shp');

% 加载栅格文件
raster = imread('your_raster_file.tif');

% 获取shapefile的边界框
bbox = [min([shp.X]), min([shp.Y]), max([shp.X]), max([shp.Y])];

% 在栅格文件中裁剪边界框
raster_cropped = mapcrop(raster, bbox);

% 保存裁剪后的栅格文件
imwrite(raster_cropped, 'your_cropped_raster_file.tif');

请注意，这将裁剪栅格文件并将结果保存为新文件。如果您希望在原始文件上进行裁剪，请确保对原始文件进行备份，并使用新文件名保存结果。

相关问题

matalb基于shp文件对栅格文件进行裁剪

您可以使用MATLAB中的Mapping Toolbox来裁剪栅格文件。以下是一个简单的示例代码：

% 加载shapefile文件
shp = shaperead('your_shapefile.shp');

% 加载栅格文件
raster = imread('your_raster_file.tif');

% 获取shapefile的边界框
bbox = [min([shp.X]), min([shp.Y]), max([shp.X]), max([shp.Y])];

% 在栅格文件中裁剪边界框
raster_cropped = mapcrop(raster, bbox);

% 显示裁剪后的栅格文件
imshow(raster_cropped);

请注意，此示例代码仅适用于裁剪矩形区域的栅格文件。如果您需要裁剪更复杂的形状，请使用Mapping Toolbox中的其他函数来实现。

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`函数进行裁剪。最后，将裁剪结果保存为栅格数据。