EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
data=load('data.txt');4 H- f. y) e7 y! g
lat=data(:,1);5 c6 ?9 Y: W* G1 S
lon=data(:,2);% N" c! I, P# U5 w
rain=data(:,3);
& S: z" T* i% T/ G% g7 `clear data
# e5 j- p* G( S, T clevel0=find(rain==0);, [0 r4 n2 g3 s: P: Y
level1=find(rain>0&rain<=0.1);( l$ I7 d. f. x9 A4 a S* J/ y
level2=find(rain>0.1&rain<=0.5);( _8 n) H/ |& Q
level3=find(rain>0.5&rain<=1);
% @* ~, Y8 U2 F6 e; alevel4=find(rain>1&rain<=5);
$ v9 @6 ^$ x) a, Elevel5=find(rain>5&rain<=10);5 a4 y, r) N4 L' g, { Z6 W
level6=find(rain>10);; N, s% ^1 T. v0 g# y
rain(level0)=nan;0 k2 B5 C8 u6 C- K: o1 ~
rain(level1)=1;
& f4 t/ p* T& y/ G+ Mrain(level2)=2;
2 |& Q j8 P8 G" Crain(level3)=3;
) `, Y( z8 B: m" i: Crain(level4)=4;
' M. a# D$ M9 l% p/ B7 [rain(level5)=5;
) I2 A# [& S) o) m; b; O$ u r* `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;7 D& A* l. J d4 Q1 w, p* N! R$ A
figure
! J9 B& q+ n$ Y Jset(gcf,'position',[100 100 800 650])' A3 N& o! R/ s/ L% Q' H
scatter(lon,lat,10,rain,'fill')2 v$ X- U: b! ?$ T
set(gca,'fontsize',18,'linewidth',1.5,'box','on') V ]! F$ J1 K1 a# \
xlim([floor(min(lon))-1 floor(max(lon))+1])# j$ o1 b: J4 ]0 y8 u) g2 L' o
ylim([floor(min(lat))-1 floor(max(lat))+1])) J: t+ d; q/ R5 `, D0 q
ylabel('Latitude')
, k+ z4 C5 L* i) q+ {% O2 Xxlabel('Longitude')1 y/ I& R0 ?5 ~0 G2 z6 w
colormap(cmap)
9 {3 z, [9 G! q& L. t- scb=colorbar;
5 ` c( J1 z1 v0 Pcaxis([-0.5 6.5])
* {# W- i7 g& }/ t6 Kset(cb,'ytick',-0.5:1:6.5,'yticklabel',{'','0','0.1','0.5','1','5','10','15'}) 用matlanb画了一张散点图,现在要根据.shp文件,加上地图,并且白化地图区域外,求大神帮帮忙 ; W3 P3 b9 }# w1 ^) J6 e
7 K3 w* C) t4 m; g% t. h |