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

Touuqu

data=load('data.txt');
lat=data(:,1);
# {) \. x' K7 I! z2 Nlon=data(:,2);
rain=data(:,3);
& h, g9 u% c4 n3 }# p4 Z: G$ z. Vclear data
. V6 [! G1 R% N' L9 `( V# Rlevel0=find(rain==0);
( f8 [- D- z4 i  B9 Klevel1=find(rain>0&rain<=0.1);
level2=find(rain>0.1&rain<=0.5);
level3=find(rain>0.5&rain<=1);
2 _9 `7 H! n0 e% H' J4 xlevel4=find(rain>1&rain<=5);
7 k0 q7 J% N" E" m0 ?! F1 I! ]level5=find(rain>5&rain<=10);
level6=find(rain>10);
* J: P/ j$ m& H3 n1 c  V  q6 Frain(level0)=nan;
% y: c+ s/ n. E; t" A  lrain(level1)=1;
& u  Z7 [+ v) i4 ^1 `rain(level2)=2;
7 Q6 o, j: V2 }6 Z% R: Vrain(level3)=3;
' d1 H) I; ^0 Y9 L  j2 Train(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;
0 U' j0 s$ P2 A, o  q# [6 ufigure
/ }; H8 O6 P% Lset(gcf,'position',[100 100 800 650])
& `7 u$ C0 f( b0 D7 Y6 @8 jscatter(lon,lat,10,rain,'fill')
/ K  `% B: l& u& J" X4 yset(gca,'fontsize',18,'linewidth',1.5,'box','on')
- C8 |8 w2 }' w& u0 H: G1 z8 Q' axlim([floor(min(lon))-1 floor(max(lon))+1])
ylim([floor(min(lat))-1 floor(max(lat))+1])
: T7 {" B5 A7 H" Iylabel('Latitude')
xlabel('Longitude')
. l0 d1 e2 z/ c) ?- S! Dcolormap(cmap)
cb=colorbar;
% ]$ R8 q5 b: V4 z* ccaxis([-0.5 6.5])
; L6 t, X9 {4 @- m7 oset(cb,'ytick',-0.5:1:6.5,'yticklabel',{'','0','0.1','0.5','1','5','10','15'})

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

yxlk

路过