找回密码
 注册
关于网站域名变更的通知
查看: 629|回复: 1
打印 上一主题 下一主题

Matlab 画图,根据shp地图文件白化???

[复制链接]
  • TA的每日心情

    2019-11-19 15:34
  • 签到天数: 1 天

    [LV.1]初来乍到

    跳转到指定楼层
    1#
    发表于 2019-7-26 15:27 | 只看该作者 |只看大图 回帖奖励 |倒序浏览 |阅读模式

    EDA365欢迎您登录!

    您需要 登录 才可以下载或查看,没有帐号?注册

    x
    data=load('data.txt');( |* S0 ]8 C" a/ E# F4 Q
    lat=data(:,1);7 N+ B; Q' D1 t2 [/ X' \+ S: D
    lon=data(:,2);$ `- W/ m, A# L+ {( [: w
    rain=data(:,3);, J3 |7 V+ W6 p1 |9 ?+ B0 s5 q
    clear data6 l6 P! F4 P4 e( D) c( _
    level0=find(rain==0);
      ]% G. g  n) V2 b7 W* ?level1=find(rain>0&rain<=0.1);
    & d+ A, n' V4 w$ o4 u1 B/ klevel2=find(rain>0.1&rain<=0.5);
    7 S( _7 p, [% d( v6 L3 o; mlevel3=find(rain>0.5&rain<=1);
    % y2 O1 ~0 l" `6 w) q/ {' D/ O, vlevel4=find(rain>1&rain<=5);) d- L$ S* l( q# y
    level5=find(rain>5&rain<=10);
    2 g% _# k8 v, klevel6=find(rain>10);
    & H" o2 [7 h7 h1 J/ ^rain(level0)=nan;) t6 j7 o' ^2 Z# n  N
    rain(level1)=1;
    * T: B3 h9 [# urain(level2)=2;
    5 K6 @$ B: T, P; \3 srain(level3)=3;0 g  i1 ~% J* d. y% ?. f$ A  k
    rain(level4)=4;# C4 |6 U: d- T2 B; R! H
    rain(level5)=5;8 M9 y) P8 A- J3 w2 X+ |2 ^. Q& l
    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;3 X( r! q9 i- e, C) b6 t
    figure
    - ?7 ?& w* L5 ~. y0 d' i6 iset(gcf,'position',[100 100 800 650])
    # P- e6 V* G6 Z% ?4 \* fscatter(lon,lat,10,rain,'fill')/ z0 G% F, m2 Y7 w  k. p. y
    set(gca,'fontsize',18,'linewidth',1.5,'box','on')4 U, L6 P& ~  o
    xlim([floor(min(lon))-1 floor(max(lon))+1])0 ]& q/ c- h* z" y. D
    ylim([floor(min(lat))-1 floor(max(lat))+1])
    3 G5 i* @4 G$ wylabel('Latitude')! L1 L8 k7 L7 |8 c1 `# v4 k2 W. x
    xlabel('Longitude')9 Z( k% W1 N) i& b6 z+ e
    colormap(cmap)7 U0 h( d! x  i+ {# _
    cb=colorbar;
    * l4 w; R4 f! l8 b1 V. m9 O$ icaxis([-0.5 6.5])
    ( j1 C( a4 w( V: T: r8 q* `  h0 mset(cb,'ytick',-0.5:1:6.5,'yticklabel',{'','0','0.1','0.5','1','5','10','15'})
    用matlanb画了一张散点图,现在要根据.shp文件,加上地图,并且白化地图区域外,求大神帮帮忙
    : f6 ?, V# N$ G: i
    , I7 T2 q( t2 u- x" q

    该用户从未签到

    2#
    发表于 2019-7-27 10:30 | 只看该作者
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

    推荐内容上一条 /1 下一条

    EDA365公众号

    关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

    GMT+8, 2025-10-5 01:07 , Processed in 0.125000 second(s), 27 queries , Gzip On.

    深圳市墨知创新科技有限公司

    地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

    快速回复 返回顶部 返回列表