|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
用Matlab画三维地形图
' w3 s- y1 ^/ Q0 b8 G: p7 S$ l5 I+ @; I' ~
代码如下：
6 Q3 M' ?- v5 d! j%%%%—————Code to draw 3D bathymetry——————————
: P% U; x9 |0 P0 W: Y9 \- S$ a j! V%-------Created by 猛牛,May 14,2012--------------------
3 V" K8 @% g) k; J% P {) K, Zclear;clf;clc;
) |: G7 s7 @/ Oncload data_bathymetry.cdf* l8 G2 z1 e6 X Z, E/ n8 |0 z
d_elev=elev(709:901,1441:1597); % Luzon Strait/ O- M& L% f- P! V
lon=X(1441:1597);9 ~2 d5 A: X l
lat=Y(709:901);
4 O$ {( j0 {& `8 r[xx,yy]=meshgrid(lon,lat);
& Z4 M5 e0 D& d3 ?%---------------------------------------------. e9 L! e5 O2 y! U7 o( K" r
view([10 -10 55]): v' D' W/ r) ~% L/ h3 v8 C
hold on
+ |% h( s" X% K6 p6 AsuRF(xx,yy,d_elev);
) ~- ~0 `. U! F0 Ishading interp/ S0 l* S9 {0 Y
axis([120 135 14 31 -7000 900])/ e" ]) @7 P o8 j1 ]; H
xlabel('Longitude','FontSize',16);6 T2 Y% v3 \9 Z! A, @, v
ylabel('Latitude','FontSize',16);7 O% p0 ~/ F M& y/ D
zlabel('Elevation(m)','FontSize',16);
: S# F, U$ i2 U- H5 ]% |h=colorbar('v');6 `2 x3 M; ^6 O F3 @6 E& g
caxis([-7000 3000]);8 V0 c% {* ^" F; U! B* o9 h
set(get(h,'ylabel'),'string','Elevation(m)','fontsize',21);9 U! y/ ^4 K& [, t
text(118,22,1000,'Taiwan Island','color','k','fontsize',16,'fontweight','bold');
1 s; u* {- A% Y1 O8 dtext(118,19,0,'Luzon Strait','color','k','fontsize',16,'fontweight','bold');+ A5 H) @$ f+ |9 f: u1 w; G
text(118,16,100,'Luzon Island','color','k','fontsize',16,'fontweight','bold');; e+ B: l: F; ~5 p
title('Bathymetry of Luzon Strait','FontSize',31,'fontweight','bold');7 }* Y3 ~7 F# f: L4 I( s3 i
|
|
/1 
发表于 2020-4-10 10:08
收藏
淘帖
支持！
反对！