matlab利用shp文件对HDF5文件裁剪

你好，你可以使用MATLAB中的Mapping Toolbox来裁剪HDF5文件。首先，你需要读入shp文件并创建一个地理参考对象。然后，你需要读入HDF5文件并将其投影到与shp文件相同的地图投影中。最后，你可以使用shaperead函数读取shp文件的多边形，并使用geoshow函数在地图上显示它们。使用inpolygon函数可以查找HDF5文件中每个像素是否位于多边形内部，并将它们裁剪掉。以下是一个可能的代码示例：

% 读入shp文件并创建地理参考对象
S = shaperead('filename.shp');
R = georasterref('RasterSize', [height width], ...
    'LatitudeLimits', [minlat maxlat], 'LongitudeLimits', [minlon maxlon]);

% 读入HDF5文件并投影到地图投影中
hdf5data = h5read('filename.h5', '/dataset');
hdf5ref = imref2d(size(hdf5data), [minlon maxlon], [minlat maxlat]);
[X, Y] = meshgrid(hdf5ref.XWorldLimits(1):hdf5ref.PixelExtentInWorldX:hdf5ref.XWorldLimits(2), ...
    hdf5ref.YWorldLimits(1):hdf5ref.PixelExtentInWorldY:hdf5ref.YWorldLimits(2));
[lat, lon] = projinv(proj, X, Y);
hdf5data = interp2(hdf5data, lon, lat);

% 绘制多边形并裁剪HDF5数据
figure
geoshow(S, 'FaceColor', 'none')
hold on
for k = 1:numel(S)
    [x, y] = projfwd(proj, S(k).X, S(k).Y);
    h = plot(x, y, 'r');
    in = inpolygon(lon, lat, S(k).X, S(k).Y);
    hdf5data(in) = NaN;
end
hold off

% 显示裁剪后的数据
figure
mapshow(hdf5data, hdf5ref, 'DisplayType', 'surface')

这个例子中，你需要设置`proj`变量为shp文件使用的地图投影。如果不知道地图投影是什么，可以使用`shaperead`函数读取shp文件中的空间参考信息。