|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
用Matlab画三维地形图( ~! n9 o0 v2 i" E d$ ?
5 D! z* w4 O7 {( z$ Q9 S5 t% I代码如下:* M+ ?: I( c0 U( U1 X. }% w7 X
%%%%—————Code to draw 3D bathymetry——————————- {% ?8 f8 y. d* T0 ~
%-------Created by 猛牛,May 14,2012--------------------
C; A2 h% Z8 X ^clear;clf;clc;# K" z1 }$ h3 J: q
ncload data_bathymetry.cdf; a1 R) Z3 D: ?* G$ ^) C" ~
d_elev=elev(709:901,1441:1597); % Luzon Strait! s' [1 @- _; j2 T$ F! l
lon=X(1441:1597);; q# |! M& N7 P a4 z& L4 k, S
lat=Y(709:901);
9 C9 W: Z! i6 o" a# W[xx,yy]=meshgrid(lon,lat);+ s3 h7 W# T. `& T
%---------------------------------------------8 D0 a% h) A" R1 b9 O
view([10 -10 55])6 R$ |, z/ ?7 B6 a* _7 c
hold on
: ]2 h6 f0 P2 qsuRF(xx,yy,d_elev);" w. _' l. [' j* m* F) ?! I" Y
shading interp
8 m: p& ^7 R$ R; gaxis([120 135 14 31 -7000 900])" ?3 d3 N% i0 {8 I
xlabel('Longitude','FontSize',16);
' ]+ Y8 Z B$ Z( eylabel('Latitude','FontSize',16);
1 \2 d0 u: W$ ?* y# c5 F2 _9 Rzlabel('Elevation(m)','FontSize',16);
! J( r7 I! M0 I- Y0 G& U2 Ch=colorbar('v'); @) F* F, b7 k, ?/ m# O, o
caxis([-7000 3000]);
0 {9 }3 Q: s. S! U+ Wset(get(h,'ylabel'),'string','Elevation(m)','fontsize',21);
4 R! w6 x! j5 g! P0 `/ vtext(118,22,1000,'Taiwan Island','color','k','fontsize',16,'fontweight','bold');
5 q m: N' D' t0 F! P* g) ?text(118,19,0,'Luzon Strait','color','k','fontsize',16,'fontweight','bold');
3 X, P7 Y: T$ |$ G o/ E, [8 }text(118,16,100,'Luzon Island','color','k','fontsize',16,'fontweight','bold');
+ `, [9 H* n! p3 M5 _title('Bathymetry of Luzon Strait','FontSize',31,'fontweight','bold');
3 F5 Q" \* D5 v |
|