根据shp文件创建Mask文件_mask.shp 廊坊

1. GIS中生成渔网

可以参考渔网生成教程。

2. 利用选择工具确定shp内部的点

3. 直接根据导出的点shp文件生成mask

% Global scale
for k = 1:length(data)
    x = data(k).X + 360;
    y = data(k).Y;
    for i = 1:360
        for j = 1:180
            if LON(i,j) == x && LAT(i,j) == y
                maskorig(i,j) = 1;
            end
        end
    end
end
1
2
3
4
5
6
7
8
9
10
11
12

% China scale
data = shaperead('SemiArid_Points.shp');
rlon = [70.25 : 0.5 : 139.75];
rlat = [15.25 : 0.5 : 54.75];
[lat, lon] = meshgrid(rlat, rlon);

maskorig = nan(140, 80);
for k = 1:length(data)
    x = data(k).X;
    y = data(k).Y;
    for i = 1:140
        for j = 1:80
            if lon(i,j) == x && lat(i,j) == y
                maskorig(i,j) = 1;
            end
        end
    end
end

h = pcolor(maskorig');
h.LineStyle = 'none';
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21