经验-matlab制作陆地掩膜

gershon

利用matlab自带landareas制作海陆掩膜。
代码如下：

1. res = 0.05; % grid size of land mask
   
2. Latlim = [-90 90];
   
3. Lonlim = [-180 180];
   
4. 
   
5. [S,~] = shaperead('landareas');
   
6. 
   
7. landmask = false(diff(Latlim)/res,diff(Lonlim)/res);
   
8. ref = makerefmat('RasterSize', size(landmask), ...
   
9.    'Latlim', Latlim, 'Lonlim',Lonlim ,...
   
10.    'ColumnsStartFrom','north','RowsStartFrom','west');
    
11. 
    
12. for i=1:length(S)
    
13.     n = find(isnan(S(i).X), 1, 'first') - 1;
    
14.     [r,c] = map2pix(ref, S(i).X(1:n), S(i).Y(1:n));
    
15.     mask = roipoly(landmask,c,r);
    
16.     landmask(mask)=true;
    
17. end
    
18. 
    
19. % geotiffwrite('landmask',uint8(landmask),ref)

复制代码

zdcgl

谢谢楼主分享！

aizaizai1989

制作好了掩模咋用呢，附带来点啥数据演示下呗？帅哥