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

Touuqu

data=load('data.txt');
lat=data(:,1);
lon=data(:,2);
# y' i! K3 b' i' P6 X* J5 Lrain=data(:,3);
4 U: _( j, S# O1 lclear data
; q  r& s  Q9 Q- W5 O% X% ]level0=find(rain==0);
: p* d' c  p+ R. j+ d( e3 C, C- glevel1=find(rain>0&rain<=0.1);
level2=find(rain>0.1&rain<=0.5);
* R# \$ i, P& x1 alevel3=find(rain>0.5&rain<=1);
$ {1 ?& w. h; O2 ]1 A/ Plevel4=find(rain>1&rain<=5);
  B( C, X' _9 |8 t4 j) \! d& vlevel5=find(rain>5&rain<=10);
+ N) K6 c& ]$ g" L: q- ylevel6=find(rain>10);
rain(level0)=nan;
$ ^1 {' U* `, ^rain(level1)=1;
6 l. c- y/ f9 f, X5 s5 S$ }rain(level2)=2;
* }. c6 Q3 a# vrain(level3)=3;
rain(level4)=4;
rain(level5)=5;
rain(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])
6 \. J- w, s- l$ L; Sscatter(lon,lat,10,rain,'fill')
set(gca,'fontsize',18,'linewidth',1.5,'box','on')
xlim([floor(min(lon))-1 floor(max(lon))+1])
ylim([floor(min(lat))-1 floor(max(lat))+1])
' W9 c9 F% `9 |ylabel('Latitude')
xlabel('Longitude')
" H/ _4 V: m5 A1 {" fcolormap(cmap)
cb=colorbar;
caxis([-0.5 6.5])
0 n  O- S) _0 @8 \6 d! ]# yset(cb,'ytick',-0.5:1:6.5,'yticklabel',{'','0','0.1','0.5','1','5','10','15'})

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

* t$ g' \0 K$ |) J: A3 }& G

yxlk

路过