EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
data=load('data.txt');3 X7 C. G3 H, n' I0 [. f& c6 a
lat=data(:,1);' Z+ k8 R4 K4 V+ o) L
lon=data(:,2);5 E) N) B# _' }" W4 l( }* _# B
rain=data(:,3);/ P/ M2 x. t: {" s V* m* n
clear data1 G' u! m6 X* g* K0 |
level0=find(rain==0);
. I) \$ k5 u+ z ^' T7 }% {level1=find(rain>0&rain<=0.1);
" e8 ]% A3 S9 }& dlevel2=find(rain>0.1&rain<=0.5);
) e3 m) q: o0 C) e" u. N# Y+ clevel3=find(rain>0.5&rain<=1);* B2 C8 O# i M5 h
level4=find(rain>1&rain<=5);; x3 c' t- s# [) ?# k! Y
level5=find(rain>5&rain<=10);
% a' h3 i5 y% G: w, d1 O% ?2 M/ ?/ Mlevel6=find(rain>10);
' V' s, G- ]6 Erain(level0)=nan;' H% [' V) R; W( r- G7 A
rain(level1)=1;
' ], \1 b9 X& E' Q arain(level2)=2;
5 P, X6 w8 n9 S1 urain(level3)=3;
: |$ ~) ~- W8 u, frain(level4)=4;# O) G, H! ~6 G
rain(level5)=5;2 Z8 }0 q5 n% V5 Y6 p
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;
- v0 s \, R% J5 `- d( {% ~, @figure
9 a6 Y. X) n7 h8 T5 Yset(gcf,'position',[100 100 800 650])2 n" l/ n3 F4 S6 M' n( a
scatter(lon,lat,10,rain,'fill')
7 P1 q) {6 ~0 d8 Zset(gca,'fontsize',18,'linewidth',1.5,'box','on')
1 c7 L% @. O+ Sxlim([floor(min(lon))-1 floor(max(lon))+1])& w8 K; W% [) t: T1 S7 u
ylim([floor(min(lat))-1 floor(max(lat))+1])6 M, ]/ j8 E8 S6 S7 ?
ylabel('Latitude')
* v6 c g3 w2 G+ s- o0 ~xlabel('Longitude')
, g; E! T1 ^" [: T5 P! ycolormap(cmap)" i- Z! T' F& h0 f' J* h% `
cb=colorbar;
0 \/ p, _" r4 N* u7 i$ _, ocaxis([-0.5 6.5])
' M2 k2 \. m t$ N' dset(cb,'ytick',-0.5:1:6.5,'yticklabel',{'','0','0.1','0.5','1','5','10','15'}) 用matlanb画了一张散点图，现在要根据.shp文件，加上地图，并且白化地图区域外，求大神帮帮忙 4 _, T+ F: k- T) w: z
! Q: g+ V! p+ }' ^. h | 
/1 
发表于 2019-7-26 15:27
收藏
淘帖
支持！
反对！