Matlab 画图，根据shp地图文件白化？？？

Touuqu

data=load('data.txt');
6 P  `0 P, S2 `/ p5 c6 vlat=data(:,1);
& Z9 |0 I* X; |" a3 x& ylon=data(:,2);
0 ]2 t8 w' b& z9 x+ x9 V0 prain=data(:,3);
. d: g/ M! K6 T3 w3 {clear data
1 q. }+ b6 X# {1 ~level0=find(rain==0);
  S1 `4 _2 g$ T9 j2 C# Qlevel1=find(rain>0&rain<=0.1);
level2=find(rain>0.1&rain<=0.5);
) H6 b* D. c0 I; Rlevel3=find(rain>0.5&rain<=1);
level4=find(rain>1&rain<=5);
level5=find(rain>5&rain<=10);
  j( }3 l" A6 o. m" o4 t7 i$ U  u& ylevel6=find(rain>10);
. P% f+ z+ g  j/ ?% g% a2 Z- arain(level0)=nan;
0 ]) @7 q4 i0 g& H5 w& wrain(level1)=1;
# X; Q, I. }: S- B1 \+ k+ ~rain(level2)=2;
* m! L7 |6 Q! W2 @6 D% J( `rain(level3)=3;
; ~$ S% ~4 X. V+ |% d/ Crain(level4)=4;
rain(level5)=5;
: i. _; p" e; L6 O4 R1 Drain(level6)=6;

cmap=[255,255,255;0,255,3;0,255,253;0,0,254;255,253,6;249,2,255;254,0,2]/255;
figure
set(gcf,'position',[100 100 800 650])
scatter(lon,lat,10,rain,'fill')
0 M6 ^4 N/ \, _+ kset(gca,'fontsize',18,'linewidth',1.5,'box','on')
xlim([floor(min(lon))-1 floor(max(lon))+1])
- r0 B. t1 m% T& _2 S* j0 {ylim([floor(min(lat))-1 floor(max(lat))+1])
ylabel('Latitude')
0 c. l0 P8 X- wxlabel('Longitude')
colormap(cmap)
cb=colorbar;
- \& K! X6 g1 [) M3 Ecaxis([-0.5 6.5])
. O. i  {/ [, W! ~set(cb,'ytick',-0.5:1:6.5,'yticklabel',{'','0','0.1','0.5','1','5','10','15'})

用matlanb画了一张散点图，现在要根据.shp文件，加上地图，并且白化地图区域外，求大神帮帮忙

& q" I: K' Q2 Z! G! @8 X! Z0 t

yxlk

路过