EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
data=load('data.txt');( |* S0 ]8 C" a/ E# F4 Q
lat=data(:,1);7 N+ B; Q' D1 t2 [/ X' \+ S: D
lon=data(:,2);$ `- W/ m, A# L+ {( [: w
rain=data(:,3);, J3 |7 V+ W6 p1 |9 ?+ B0 s5 q
clear data6 l6 P! F4 P4 e( D) c( _
level0=find(rain==0);
]% G. g n) V2 b7 W* ?level1=find(rain>0&rain<=0.1);
& d+ A, n' V4 w$ o4 u1 B/ klevel2=find(rain>0.1&rain<=0.5);
7 S( _7 p, [% d( v6 L3 o; mlevel3=find(rain>0.5&rain<=1);
% y2 O1 ~0 l" `6 w) q/ {' D/ O, vlevel4=find(rain>1&rain<=5);) d- L$ S* l( q# y
level5=find(rain>5&rain<=10);
2 g% _# k8 v, klevel6=find(rain>10);
& H" o2 [7 h7 h1 J/ ^rain(level0)=nan;) t6 j7 o' ^2 Z# n N
rain(level1)=1;
* T: B3 h9 [# urain(level2)=2;
5 K6 @$ B: T, P; \3 srain(level3)=3;0 g i1 ~% J* d. y% ?. f$ A k
rain(level4)=4;# C4 |6 U: d- T2 B; R! H
rain(level5)=5;8 M9 y) P8 A- J3 w2 X+ |2 ^. Q& l
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;3 X( r! q9 i- e, C) b6 t
figure
- ?7 ?& w* L5 ~. y0 d' i6 iset(gcf,'position',[100 100 800 650])
# P- e6 V* G6 Z% ?4 \* fscatter(lon,lat,10,rain,'fill')/ z0 G% F, m2 Y7 w k. p. y
set(gca,'fontsize',18,'linewidth',1.5,'box','on')4 U, L6 P& ~ o
xlim([floor(min(lon))-1 floor(max(lon))+1])0 ]& q/ c- h* z" y. D
ylim([floor(min(lat))-1 floor(max(lat))+1])
3 G5 i* @4 G$ wylabel('Latitude')! L1 L8 k7 L7 |8 c1 `# v4 k2 W. x
xlabel('Longitude')9 Z( k% W1 N) i& b6 z+ e
colormap(cmap)7 U0 h( d! x i+ {# _
cb=colorbar;
* l4 w; R4 f! l8 b1 V. m9 O$ icaxis([-0.5 6.5])
( j1 C( a4 w( V: T: r8 q* ` h0 mset(cb,'ytick',-0.5:1:6.5,'yticklabel',{'','0','0.1','0.5','1','5','10','15'}) 用matlanb画了一张散点图,现在要根据.shp文件,加上地图,并且白化地图区域外,求大神帮帮忙 : f6 ?, V# N$ G: i
, I7 T2 q( t2 u- x" q
|