|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
用Matlab画三维地形图6 [ Y- j2 |0 Y
2 o( Y9 H6 \( f
代码如下:
. M' F5 k* e( o$ V/ E0 q# x. I%%%%—————Code to draw 3D bathymetry—————————— |# ~# o& k9 P" Q% b# Q
%-------Created by 猛牛,May 14,2012--------------------3 J7 J- m$ W9 F
clear;clf;clc;4 g6 _1 T1 ? ^1 M% {: O( |
ncload data_bathymetry.cdf
) m4 ~! b6 S5 g+ ?d_elev=elev(709:901,1441:1597); % Luzon Strait
" V+ a% _0 d7 H0 [ ]: R0 a! olon=X(1441:1597);
* ]& A6 F1 @; C Olat=Y(709:901);
3 q( f2 |) Q8 _" I; y; u) n8 i[xx,yy]=meshgrid(lon,lat); O3 k! v: b4 }) L6 u' ?" M0 G
%---------------------------------------------
- w, n! A9 T! E. _view([10 -10 55])" k1 J' a; }& x+ U" ?
hold on
7 h' o5 u7 B2 z9 H4 psuRF(xx,yy,d_elev);
& S: Y" k" @1 U# G" {$ [# Tshading interp4 q n6 w3 I7 K9 v5 u( N/ r
axis([120 135 14 31 -7000 900])9 o/ p% a" M- f) m/ g& ^
xlabel('Longitude','FontSize',16);' q8 m& [( e. L
ylabel('Latitude','FontSize',16);
; l1 y* Z/ f# r0 ezlabel('Elevation(m)','FontSize',16);- J9 Q" _- E7 C3 l
h=colorbar('v');
+ A2 E. A9 }; t5 f, ^/ q, Kcaxis([-7000 3000]);/ |6 o: w8 @5 H' D2 t% M
set(get(h,'ylabel'),'string','Elevation(m)','fontsize',21);9 Y. j- W% K3 U9 L5 C8 d1 c4 p
text(118,22,1000,'Taiwan Island','color','k','fontsize',16,'fontweight','bold');
, y: G0 C! q7 gtext(118,19,0,'Luzon Strait','color','k','fontsize',16,'fontweight','bold'); }: t4 f h& C
text(118,16,100,'Luzon Island','color','k','fontsize',16,'fontweight','bold');
" z% w1 r" t8 ^title('Bathymetry of Luzon Strait','FontSize',31,'fontweight','bold');/ j+ P5 `. U- a: D# ~; V: m
|
|