在matlab中，我們有時需要創(chuàng)建180*360,360*720,720*1440等不同分辨率的柵格掩膜，用于提取特定研究區(qū)域的數(shù)據(jù)。

貼上代碼：

lon = -179.875:0.25:179.875; %分辨率0.25

lat = -89.875:0.25:89.875;

% lon = 0:1:360;

% lat = -90:1:90;

[lon,lat] = meshgrid(lon,lat);

nc_plain_boundary1=load('taiwan.txt');??

nc_plain_boundary2=load('hainan.txt');??

nc_plain_boundary3=load('mainland.txt');??

% data

ind1=inpolygon(lon,lat,nc_plain_boundary1(:,1),nc_plain_boundary1(:,2));

ind2=inpolygon(lon,lat,nc_plain_boundary2(:,1),nc_plain_boundary2(:,2));

ind3=inpolygon(lon,lat,nc_plain_boundary3(:,1),nc_plain_boundary3(:,2));

t = ind2+ind1+ind3;

% t = ind1;

u.lon = lon;

u.lat = lat;

u.rg = t;

rg_plot(u); %繪圖函數(shù)

【其中的taiwan.txt hainan.txt mainland.txt是邊界數(shù)據(jù)，是經(jīng)緯度兩列的數(shù)據(jù)】

得到如下的結(jié)果：