|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
用Matlab画三维地形图
$ M- C& w d7 k$ L# G
: S' d9 G$ i& \" Z- E5 z3 n' \8 H2 G代码如下:
! D' N2 [: B2 \%%%%—————Code to draw 3D bathymetry——————————3 c$ [* W9 }2 `6 R0 q/ A u
%-------Created by 猛牛,May 14,2012--------------------/ Y' M' P, v3 s2 R6 {* J
clear;clf;clc;% ?' d7 _$ p# c+ K) B: a4 v% P' @& A
ncload data_bathymetry.cdf* X6 {* Y9 |/ K1 R3 l# _
d_elev=elev(709:901,1441:1597); % Luzon Strait
8 H" f1 u$ `. w- Z) r) G7 _lon=X(1441:1597);
: g4 ^9 Y( s$ R2 X. D) a7 D9 ~lat=Y(709:901);
: t1 z. z( T% n8 r- P: T: ~+ f. P[xx,yy]=meshgrid(lon,lat);4 U9 P, P c% w) A! Y2 P
%---------------------------------------------! s! e' n" F) k
view([10 -10 55])
! @. h' C7 ]& o; D2 ^- K# ]hold on
Z1 U8 F1 r5 vsuRF(xx,yy,d_elev);% P: B: U7 W& m8 Y; c
shading interp7 T) _ L- h- l" H
axis([120 135 14 31 -7000 900])
- s3 y4 `" g& c) s% Mxlabel('Longitude','FontSize',16);
* e' v. O' _. Z! b( P) Yylabel('Latitude','FontSize',16);- ~- \, {+ ^$ f# v
zlabel('Elevation(m)','FontSize',16);
; x! N. b% ?9 _8 q0 b4 Y4 E' Ah=colorbar('v');
1 i4 L8 ~+ y" L9 Rcaxis([-7000 3000]);5 @) G: C! R2 ^. l6 L
set(get(h,'ylabel'),'string','Elevation(m)','fontsize',21);
4 N! Q. t! K [text(118,22,1000,'Taiwan Island','color','k','fontsize',16,'fontweight','bold');
' d* ~( E; ]7 ttext(118,19,0,'Luzon Strait','color','k','fontsize',16,'fontweight','bold');2 H7 g, f- Z# _% b4 a
text(118,16,100,'Luzon Island','color','k','fontsize',16,'fontweight','bold');
" V4 s0 K& Y1 Ytitle('Bathymetry of Luzon Strait','FontSize',31,'fontweight','bold');
* b- D/ o! I. E, J |
|