EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
data=load('data.txt');
! e& z8 |: _/ l; W M- nlat=data(:,1);1 u: D1 A- s9 C
lon=data(:,2);
; a4 p. A5 `4 E. o( rrain=data(:,3);4 t6 v# x4 N. [, _
clear data
/ E& Z8 R* n( ^level0=find(rain==0);
+ P, ]( @: B6 l9 s+ S' l3 klevel1=find(rain>0&rain<=0.1);
, ^: L& P0 l" q- r( n* e+ Ulevel2=find(rain>0.1&rain<=0.5);
$ J+ O; g6 H; xlevel3=find(rain>0.5&rain<=1);$ y1 h: O$ L% }3 C3 |
level4=find(rain>1&rain<=5);( Y- a) m a T; e X
level5=find(rain>5&rain<=10);
% y4 N& A1 a5 T/ ?- m8 x; ulevel6=find(rain>10);
2 m/ V$ {) W% Q7 xrain(level0)=nan;: F6 i* p- ?" j, j- c* g
rain(level1)=1;: T3 v# o ^+ x4 }( b
rain(level2)=2;9 q7 g. L+ _: c" |9 o3 r2 [: P1 Y+ K
rain(level3)=3;
* O; y, \) E" Xrain(level4)=4; Q6 h2 |, L# P7 g
rain(level5)=5;3 D, B0 p8 \# k. i
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;. g0 k& l' P4 b
figure
+ k4 }/ Z9 b: N9 {8 N Z2 vset(gcf,'position',[100 100 800 650])% U0 t- a" g: X2 d9 m4 _
scatter(lon,lat,10,rain,'fill')
8 _3 w9 l4 D2 {1 @9 N+ Q7 [5 t+ Jset(gca,'fontsize',18,'linewidth',1.5,'box','on'), l6 J& R. b0 I7 k1 o" b
xlim([floor(min(lon))-1 floor(max(lon))+1]); G' C7 j+ T8 x9 c5 }6 N9 B# i* K$ v
ylim([floor(min(lat))-1 floor(max(lat))+1])
$ o4 L& ^& @* D+ I0 L5 aylabel('Latitude')
$ j$ s' k5 G* D7 P6 Xxlabel('Longitude')
7 M* U; W; g( @: z B2 X: F. ?+ `colormap(cmap)- ]1 }- G8 y) T: l
cb=colorbar;
% J. Q+ W2 Q! w: K' N# B+ jcaxis([-0.5 6.5])& v8 t# X* x3 O/ y( g1 v% l
set(cb,'ytick',-0.5:1:6.5,'yticklabel',{'','0','0.1','0.5','1','5','10','15'}) 用matlanb画了一张散点图，现在要根据.shp文件，加上地图，并且白化地图区域外，求大神帮帮忙
' D9 |# Y3 ~! O
5 K% Y( s, W; n | 
/1 
发表于 2019-7-26 15:27
收藏
淘帖
支持！
反对！