本質(zhì)上是用shp文件對(duì)范圍以外的數(shù)據(jù)取空值,既可以用于繪圖數(shù)據(jù),也可用于河流地形等數(shù)據(jù),僅需要提供一組經(jīng)緯度
Example
% 安徽shp
region_ah = shaperead('S:\5.模式與編程學(xué)習(xí)\重要的畫圖數(shù)據(jù)\全國(guó)各省矢量圖(全)-by 朱燕\anhui.shp');
% 河流文件
hyd2 = shaperead('S:\5.模式與編程學(xué)習(xí)\重要的畫圖數(shù)據(jù)\國(guó)家基礎(chǔ)地理數(shù)據(jù)\hyd2_4m\hyd2_4l.shp');
hyd2_4lx=[hyd2(:).X];hyd2_4ly=[hyd2(:).Y];
% 利用安徽shp白化河流文件
inshp = inpolygon(hyd2_4lx,hyd2_4ly,region_ah.X,region_ah.Y);
hyd2_4lx = hyd2_4lx.*inshp;hyd2_4ly = hyd2_4ly.*inshp;
hyd2_4lx(hyd2_4lx==0) = nan;
hyd2_4ly(hyd2_4ly==0) = nan;
% 繪制地圖
m_plot(region_ah.X,region_ah.Y,'Color', 'k' ,'LineWidth',2);hold on; %安徽省行政邊界?
m_plot(hyd2_4lx,hyd2_4ly,'color',[0.2 0.4 0.9],'LineWidth',1.5);hold on; %繪制一級(jí)河流