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

Touuqu

data=load('data.txt');
/ v& N& n7 g1 V! A* Q9 clat=data(:,1);
! I& z& R6 v; n" l$ ]lon=data(:,2);
% g8 F/ s$ v% D& r8 rrain=data(:,3);
% ], {+ K! g, D' X. l* jclear data
& c3 d' v6 `9 ~( _7 rlevel0=find(rain==0);
level1=find(rain>0&rain<=0.1);
level2=find(rain>0.1&rain<=0.5);
$ _# N$ }  e8 N5 Zlevel3=find(rain>0.5&rain<=1);
- o7 W6 y5 r: alevel4=find(rain>1&rain<=5);
level5=find(rain>5&rain<=10);
level6=find(rain>10);
rain(level0)=nan;
+ `0 b- R4 m. b0 [8 O7 S1 D" G: Arain(level1)=1;
rain(level2)=2;
rain(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;
- @' ?+ z  @/ U/ j4 k! Ifigure
, m* \2 }7 y( i- `1 r6 \' wset(gcf,'position',[100 100 800 650])
# }2 R+ C$ Y  H* G8 Sscatter(lon,lat,10,rain,'fill')
set(gca,'fontsize',18,'linewidth',1.5,'box','on')
. Y/ m7 r) ?# T& ?  O( Zxlim([floor(min(lon))-1 floor(max(lon))+1])
. z% ~3 i0 M% d* C3 F4 \9 `ylim([floor(min(lat))-1 floor(max(lat))+1])
# Q7 Y( [! U" o5 ]# I7 Cylabel('Latitude')
xlabel('Longitude')
4 }6 J; G6 `$ n. T- Ecolormap(cmap)
cb=colorbar;
( ]1 A9 N: |$ zcaxis([-0.5 6.5])
set(cb,'ytick',-0.5:1:6.5,'yticklabel',{'','0','0.1','0.5','1','5','10','15'})

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

! e7 b, v4 D) z" @" u- Q* k' n) r( C& _

yxlk

路过