matlab使用shp文件裁剪TIF影像

可以使用MATLAB中的Mapping Toolbox中的函数来实现shp文件对TIF影像的裁剪。

以下是一个简单的例子：

% 读取shp文件
s = shaperead('example.shp');

% 读取TIF文件
[I, R] = geotiffread('example.tif');

% 将TIF文件的地理参考信息与shp文件匹配
refmat = worldfileread('example.tfw', 'planar', size(I));
RA = maprefcells(R, size(I), 'ColumnsStartFrom', 'north');

% 将shp文件转换成地理坐标系
x = [s.X];
y = [s.Y];
[x, y] = projinv(s(1).BoundingBox.CRS, x, y);

% 将shp文件转换成像素坐标系
[r, c] = map2pix(RA, x, y);

% 创建一个mask矩阵
mask = poly2mask(r, c, size(I, 1), size(I, 2));

% 裁剪TIF影像
I = I .* uint16(mask);

% 显示裁剪后的TIF影像
imshow(I);

这个例子假设你有一个名为`example.shp`的shp文件和一个名为`example.tif`的TIF影像。你需要将`example.tfw`文件放在与`example.tif`相同的文件夹中，并且该文件必须包含TIF文件的地理参考信息。你可以根据自己的实际情况修改这个例子，以适应你的数据。