EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
data=load('data.txt');- L! s* L! k9 L. s( a; E
lat=data(:,1);. Y3 L6 n/ r, j; N: M
lon=data(:,2);0 N- x4 E2 n# W; p' M
rain=data(:,3);8 Q6 c# y7 F& }9 o
clear data' R8 p/ _( l- I- r% j: t/ | H& f$ n
level0=find(rain==0);
X6 N4 Q2 M; blevel1=find(rain>0&rain<=0.1);* _5 b( y2 ^2 V' {- M* p: h4 h
level2=find(rain>0.1&rain<=0.5);: K0 o, ^ p0 p! E
level3=find(rain>0.5&rain<=1);4 n p$ E6 g0 G R% ]2 I) R
level4=find(rain>1&rain<=5);
. t2 y; m- J8 |! zlevel5=find(rain>5&rain<=10);* G) f% t0 I' ?9 I1 Q5 \; I7 X; [
level6=find(rain>10);
. ]1 i7 a( F3 Nrain(level0)=nan;
6 ], |% t8 c7 [. Q3 Srain(level1)=1;$ ?( X/ p5 P5 s5 u$ }
rain(level2)=2;: T: C2 S, S! R" @. ]9 ]
rain(level3)=3;
$ u4 `/ B9 {3 \4 D- u' a- p- |rain(level4)=4;# m. K+ Z9 p' c: m4 Q
rain(level5)=5;2 l7 }. U9 O7 S- z. b% k4 c% y" F
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 ]3 X. {/ p7 \! |% m- _) X
figure$ I9 X/ H& b! Y# Z) U e$ b
set(gcf,'position',[100 100 800 650]): H9 p w$ P3 S, ?+ Q* t1 Q& `
scatter(lon,lat,10,rain,'fill') c. F/ o$ C3 X& S* o7 P- b
set(gca,'fontsize',18,'linewidth',1.5,'box','on') |+ f" M) `& q$ G1 k4 J: ~. N
xlim([floor(min(lon))-1 floor(max(lon))+1]), `' w3 ]5 r0 H) S
ylim([floor(min(lat))-1 floor(max(lat))+1])2 q6 s- t! e- Q5 U. c
ylabel('Latitude')
" H O& H; C: a6 W axlabel('Longitude')' [5 Q+ Y6 {" e# B
colormap(cmap)2 f; [* V, s2 {# T/ B7 U/ C% l$ P
cb=colorbar;; \7 B0 \" q" }2 N7 v7 u
caxis([-0.5 6.5]) g! ]. N" |# b8 M4 l9 x
set(cb,'ytick',-0.5:1:6.5,'yticklabel',{'','0','0.1','0.5','1','5','10','15'}) 用matlanb画了一张散点图，现在要根据.shp文件，加上地图，并且白化地图区域外，求大神帮帮忙
, r7 l0 G- a' c# N5 T
2 V& F0 E, K: F) U7 z, d
| 
/1 
发表于 2019-7-26 15:27
收藏
淘帖
支持！
反对！