EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
data=load('data.txt'); _% A9 K( a f5 P! i6 A
lat=data(:,1);. O! l0 f. j. b$ ^/ c7 k. m3 G7 v2 @
lon=data(:,2);3 u5 T6 M" g, U9 B; c4 V
rain=data(:,3);( i! t) f$ ]' e4 ^
clear data
% c5 B( Z8 F: y) L: `- Vlevel0=find(rain==0);
* N# H% Z$ z' L) `3 f; |level1=find(rain>0&rain<=0.1);
: d+ u7 v( k4 q" `level2=find(rain>0.1&rain<=0.5);
% U. ?% L- i; W d; v; Qlevel3=find(rain>0.5&rain<=1);
# S7 e5 {! e) a' Q( I+ ?% {/ a* d) ?level4=find(rain>1&rain<=5);
1 q, |5 {( P: Q qlevel5=find(rain>5&rain<=10);; @* A, X+ m( r' j
level6=find(rain>10);
+ y6 j4 e: ?7 T) u! @rain(level0)=nan;6 K- ^$ k; V: R1 V# x
rain(level1)=1;
( ~8 V7 p6 S; x: P3 N" mrain(level2)=2;
9 e0 ?* M a. u6 f# R" irain(level3)=3;- |7 ?! q! O* \
rain(level4)=4;1 r Y- _. d& ]2 p+ g7 q
rain(level5)=5;" j. `. b1 S; d
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 a n, S+ \7 _. V, Ffigure
1 ?, }( D0 y2 i* n3 s5 bset(gcf,'position',[100 100 800 650])8 b: R# I; C5 I% z9 A! J7 [
scatter(lon,lat,10,rain,'fill')
8 g! q! S6 q8 R; @4 mset(gca,'fontsize',18,'linewidth',1.5,'box','on')
" Y* m0 B0 [# f1 L" b$ W4 j. Cxlim([floor(min(lon))-1 floor(max(lon))+1])8 Y! t2 w6 s% A8 E( X7 R
ylim([floor(min(lat))-1 floor(max(lat))+1])
4 ^6 M3 A/ ~4 tylabel('Latitude')6 M$ S3 A7 ]. ?& }
xlabel('Longitude')
6 s4 z9 X5 B* Q8 xcolormap(cmap)
U2 `$ R7 [9 e s$ @cb=colorbar;- @2 m% G9 u5 U+ W8 H8 \
caxis([-0.5 6.5])
3 u7 p$ F! Y8 eset(cb,'ytick',-0.5:1:6.5,'yticklabel',{'','0','0.1','0.5','1','5','10','15'}) 用matlanb画了一张散点图，现在要根据.shp文件，加上地图，并且白化地图区域外，求大神帮帮忙 3 f- ~' `; q/ j4 c, J, O
; G8 w; E1 W6 p/ _1 H: x$ E: A | 
/1 
发表于 2019-7-26 15:27
收藏
淘帖
支持！
反对！