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

Touuqu

data=load('data.txt');
! B$ M; |- d$ B6 E2 Nlat=data(:,1);
% I0 ?( \) d4 L- ^lon=data(:,2);
2 k  Z* K: K2 I; f( k! Q& s9 Srain=data(:,3);
# r  Y7 [, B4 m- p* Jclear data
level0=find(rain==0);
( O; v$ i& \! t6 ]# _% P0 ?( Klevel1=find(rain>0&rain<=0.1);
4 |, d0 ^6 r+ v( i8 g5 z. blevel2=find(rain>0.1&rain<=0.5);
) n% Z$ u9 q& x- p* Alevel3=find(rain>0.5&rain<=1);
level4=find(rain>1&rain<=5);
; e3 K% r: h6 Q6 slevel5=find(rain>5&rain<=10);
level6=find(rain>10);
3 \+ n9 {" q% J3 zrain(level0)=nan;
* |: H1 S! x. D4 d/ b  srain(level1)=1;
8 [5 Q& R1 @6 Z! `1 {rain(level2)=2;
# g$ p$ |: s4 d( rrain(level3)=3;
$ K+ e6 o( r( D/ F8 d, xrain(level4)=4;
rain(level5)=5;
9 j1 ^* P# ~2 [& mrain(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])
4 t0 j7 U/ ]. ]0 p& _) G! l; O) uscatter(lon,lat,10,rain,'fill')
set(gca,'fontsize',18,'linewidth',1.5,'box','on')
2 M/ O7 i; ?1 Dxlim([floor(min(lon))-1 floor(max(lon))+1])
# G' A+ x$ O! ?+ qylim([floor(min(lat))-1 floor(max(lat))+1])
8 P* o6 X/ y! V% gylabel('Latitude')
, j$ G7 u- \- D9 I* y6 nxlabel('Longitude')
) s7 h! Z& C+ Y. [. ncolormap(cmap)
4 [/ c& J4 N: Acb=colorbar;
caxis([-0.5 6.5])
set(cb,'ytick',-0.5:1:6.5,'yticklabel',{'','0','0.1','0.5','1','5','10','15'})

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

yxlk

路过