|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
用Matlab画三维地形图
2 x5 X% e% [6 {+ g% d& }- D6 M: ]
3 ]/ f% @( p0 M: q& ]代码如下：$ i6 p5 r3 h% O* g
%%%%—————Code to draw 3D bathymetry——————————* s: R8 B. S, a4 w) ~
%-------Created by 猛牛,May 14,2012--------------------% ^' r O# g" }6 k0 M
clear;clf;clc;
4 U" [& f9 _% C: W4 @/ Xncload data_bathymetry.cdf
3 Z- o* P$ E1 |0 T: r3 Q/ T& bd_elev=elev(709:901,1441:1597); % Luzon Strait
9 U/ u4 Q4 C: r0 O. S" {6 Zlon=X(1441:1597);" C8 F9 z2 j! G) H
lat=Y(709:901);; G3 L' K6 C/ k1 G [* y
[xx,yy]=meshgrid(lon,lat);, t- u2 t7 W" ~7 y4 B& X" c
%---------------------------------------------4 G \7 b+ S0 k1 \' Y6 p
view([10 -10 55])
# o) z$ t* k* x, d$ N) Rhold on
) }/ H; ?- f: K7 ]suRF(xx,yy,d_elev);
! s& V& N6 a6 ~* a3 j9 nshading interp8 m0 n4 [/ t1 J, E3 H9 a e7 J
axis([120 135 14 31 -7000 900])3 B! F2 ` F- L2 q8 Y2 C' ]
xlabel('Longitude','FontSize',16);5 u8 O$ t- \5 L. R2 H7 N
ylabel('Latitude','FontSize',16);& p9 r3 M2 k. R9 p6 @# o; x
zlabel('Elevation(m)','FontSize',16);* k4 J6 o$ R2 H; z9 _' l4 P1 t
h=colorbar('v');
9 r9 N' l( x/ |6 P4 }5 P+ ^caxis([-7000 3000]);
" R6 {$ t+ r: q* gset(get(h,'ylabel'),'string','Elevation(m)','fontsize',21);) ]4 o' [- _/ ~0 A9 Z3 U1 d, C; L
text(118,22,1000,'Taiwan Island','color','k','fontsize',16,'fontweight','bold');0 y4 I; a; S5 L( D
text(118,19,0,'Luzon Strait','color','k','fontsize',16,'fontweight','bold');: _/ M& q1 y0 g; e+ h
text(118,16,100,'Luzon Island','color','k','fontsize',16,'fontweight','bold');1 u4 m: n' Y; Z6 e
title('Bathymetry of Luzon Strait','FontSize',31,'fontweight','bold');( u" _$ B5 X; I# V* e( S1 C6 O
|
|
/1 
发表于 2020-4-10 10:08
收藏
淘帖
支持！
反对！