|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
用Matlab画三维地形图' ]3 b0 E/ F |' _( c
- E$ k% |2 X; N9 g# S代码如下:. `9 U- b. j& R
%%%%—————Code to draw 3D bathymetry——————————
& O1 U& p& ^" U t7 i%-------Created by 猛牛,May 14,2012--------------------
$ K! W3 h4 B" o; F/ F! b" jclear;clf;clc;7 `! ?# j# q1 s+ e% N1 `% |2 ^) r2 }' d
ncload data_bathymetry.cdf
( {+ c& }, [4 `& V) hd_elev=elev(709:901,1441:1597); % Luzon Strait
- m: }' O- J% P0 @7 A" \lon=X(1441:1597);1 H' W' T/ O& }
lat=Y(709:901);
* P8 u- Q: v6 g9 _7 Y- Z[xx,yy]=meshgrid(lon,lat);
5 o% \( j8 o9 b8 n& ~1 c%---------------------------------------------
# D/ m e" N1 E- {# t" m; {view([10 -10 55]); K$ H4 H3 ]" l0 j
hold on- [) T7 `0 }) }) z/ c4 ]% y
suRF(xx,yy,d_elev); V2 a3 H4 r0 o4 M
shading interp
! Q, t( a& h5 g* Jaxis([120 135 14 31 -7000 900])0 `5 [0 l2 I& O' U1 g5 y2 g9 H) z# O( \
xlabel('Longitude','FontSize',16);: r1 U& Q& a. Y" Q
ylabel('Latitude','FontSize',16);
1 K$ I! z4 Z% y s/ m/ O+ czlabel('Elevation(m)','FontSize',16);
8 j8 P2 o6 _3 u* {1 `5 Fh=colorbar('v');
, M; S. h6 d) [" Tcaxis([-7000 3000]);4 `4 Z, V( h& K& j7 L
set(get(h,'ylabel'),'string','Elevation(m)','fontsize',21);, C1 ? o; k. n# l" q: k
text(118,22,1000,'Taiwan Island','color','k','fontsize',16,'fontweight','bold');3 F7 @! a- j' x
text(118,19,0,'Luzon Strait','color','k','fontsize',16,'fontweight','bold');8 {5 ?" m6 S: C q
text(118,16,100,'Luzon Island','color','k','fontsize',16,'fontweight','bold');8 G Y7 W5 Y' ~, X2 E
title('Bathymetry of Luzon Strait','FontSize',31,'fontweight','bold');, s( d& o$ A% D' f$ j! |9 w% w
|
|