EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
data=load('data.txt');
4 V( w( e; Y B( mlat=data(:,1);
6 B/ ?6 W0 P3 d( J6 b1 ?, ]- ylon=data(:,2);
g+ \4 s% \0 ~+ n# R. train=data(:,3);. Q0 ?; W0 R* g1 b( J
clear data& n0 S- N" E/ H- j* d- h) V1 K, b& x
level0=find(rain==0);. v' G, Q! e9 c6 ?' ^
level1=find(rain>0&rain<=0.1);
* _2 j: b! C& l, ~0 q6 q' }level2=find(rain>0.1&rain<=0.5);. j; p7 h6 r( _ l. S9 g
level3=find(rain>0.5&rain<=1);
2 r* v8 c, T$ z7 E' Mlevel4=find(rain>1&rain<=5);
& @+ I+ }0 l1 N, P$ Alevel5=find(rain>5&rain<=10);7 k) P" m. }) D
level6=find(rain>10);
1 ? T1 K% @) [rain(level0)=nan;
1 ^) T, h1 a/ G# Irain(level1)=1; m; c, Q2 N; ?# ?
rain(level2)=2;
- L, Y0 Y4 m1 X" Irain(level3)=3;5 D! o) P0 V9 i9 V6 E! c
rain(level4)=4;- @( L: F/ J4 i. M
rain(level5)=5;3 _; E' M# J0 _, h' ~8 f7 q
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;
2 _1 P7 f" q" w$ y5 _figure
$ \ _8 f/ V" y0 t; t8 z' O0 rset(gcf,'position',[100 100 800 650])
0 a I- R( w( I: J. ^; C6 `# i; mscatter(lon,lat,10,rain,'fill')1 w1 V" s: }; S& w
set(gca,'fontsize',18,'linewidth',1.5,'box','on')
5 m$ i: c, E: _" Txlim([floor(min(lon))-1 floor(max(lon))+1])
7 E p/ V* B R( C; O9 `" I8 _. I' Jylim([floor(min(lat))-1 floor(max(lat))+1]) e( m" I! K4 U' \1 z8 x$ E; D4 f
ylabel('Latitude')" ^( B- }: N) a( p; f
xlabel('Longitude')/ A) K3 b b& `* F
colormap(cmap)
8 H( j' k4 Q4 x' h. K/ g7 u$ pcb=colorbar;. t$ d& L/ d4 F8 u/ L5 i
caxis([-0.5 6.5])
6 [* Z$ l, `- v+ h" o6 s% rset(cb,'ytick',-0.5:1:6.5,'yticklabel',{'','0','0.1','0.5','1','5','10','15'}) 用matlanb画了一张散点图，现在要根据.shp文件，加上地图，并且白化地图区域外，求大神帮帮忙 ; K a y; T" T9 `9 U. Z/ S
) Z6 g+ h1 @; h
| 
/1 
发表于 2019-7-26 15:27
收藏
淘帖
支持！
反对！