|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
用Matlab画三维地形图, N8 W% F6 p$ [* e/ z
7 N! _: B" t/ z; y代码如下：
H% N9 [- b" N* V7 |. q. H) r# k7 i%%%%—————Code to draw 3D bathymetry——————————
0 q/ b# l* E$ N# C! ?+ }% J/ ?$ A%-------Created by 猛牛,May 14,2012--------------------
: n9 d) L; O0 ]6 N3 P* [( R$ eclear;clf;clc;6 n2 N) W& V$ i( n$ i8 d
ncload data_bathymetry.cdf. ^: `1 n: Z O7 r7 Y3 R& J# r/ m
d_elev=elev(709:901,1441:1597); % Luzon Strait
8 G: w+ v/ R9 d) }lon=X(1441:1597);0 l8 C7 ~/ d) Z4 M8 u
lat=Y(709:901);
+ U) x _; g' X+ ^, o[xx,yy]=meshgrid(lon,lat);0 G, y$ r' R+ D" W1 [& S5 [
%---------------------------------------------
* e7 e7 ^- c4 i! {- D0 @view([10 -10 55])
/ W1 \& H4 t& ^8 h+ J6 Thold on. g. l0 [3 Y l: }5 u! Z0 l- d' I; q
suRF(xx,yy,d_elev);
/ [$ ]5 f* K7 ushading interp
9 c2 h0 u$ A: o4 Qaxis([120 135 14 31 -7000 900])3 v, e$ d" ^0 n: Q% F" _" W
xlabel('Longitude','FontSize',16);9 K' `# a) J9 Q8 s/ W
ylabel('Latitude','FontSize',16);
& @- W0 {8 N9 d7 }zlabel('Elevation(m)','FontSize',16);
v, s# |! X$ z7 E \$ [# Yh=colorbar('v');
' s) a: p2 c' ^caxis([-7000 3000]);
. f& X1 A2 E- i. x7 ?: Z& sset(get(h,'ylabel'),'string','Elevation(m)','fontsize',21);
4 D8 `5 V0 b d! H# M0 ptext(118,22,1000,'Taiwan Island','color','k','fontsize',16,'fontweight','bold');
" w! [( S: M. r, C# U6 {) F; ^text(118,19,0,'Luzon Strait','color','k','fontsize',16,'fontweight','bold');: Y/ l8 b/ `3 e- B, y
text(118,16,100,'Luzon Island','color','k','fontsize',16,'fontweight','bold');( w! a1 R1 H2 m+ W1 k) o. q4 H1 s
title('Bathymetry of Luzon Strait','FontSize',31,'fontweight','bold');
1 {& i- K. W9 z: N; k+ u: k |
|
/1 
发表于 2020-4-10 10:08
收藏
淘帖
支持！
反对！