1.进入网站www.52ocean.cn/ u0 ?6 q: m( }9 B E# F- Y
8 S% d. D+ @" d
2.点击进入数据下载界面
/ Y7 p f! N2 _. O: K5 m" Z t" _8 o8 X' O# P& I
3 T. @: e2 O4 o! I
3.设置需要提取的地形数据经纬度范围,由于数据精度比较大, 大范围的数据需要分块进行下载,设置好之后,点击“get data”) F( {6 w+ ]5 f: |6 b, c) K
+ h& t3 i0 ~; Y9 ^9 q) z, h
1 q! }; e: s* ^+ d! t
9 j3 H8 s6 O- x: L1 E4. 将数据拷贝到文本中,例如‘.txt’
& z3 r# Z4 h8 j. m( l8 s【附】Matlab对上述数据进行网格化处理的程序:
) M9 J3 [. I8 ]" H! k- Y& R0 Hclear all;clc& h, H. B$ y* c" u O" e
load('E:\Dataprocesstopodatatopo_raw.txt')
/ i/ k+ ~& O! N+ p$ n8 DLon=sort(unique(topo_raw(:,1)),'ascend'); e& j8 ]/ |$ K( U$ i9 j6 k
Lat=sort(unique(topo_raw(:,2)),'ascend');: X& Z* J6 L( E, p8 l
lon=repmat(Lon',[length(Lat),1]);" D; I I- z- K8 ^
lat=repmat(Lat,[1,length(Lon)]);
" }1 t1 P; V% [$ d k# S: Gfor i=1:length(Lat)+ P' E% a0 ?( N7 |" h
for j=1:length(Lon)9 ]! u# } r7 e% T# r. X
idx = find(abs(topo_raw(:,1)-lon(i,j))>=0 & abs(topo_raw(:,1)-lon(i,j))<0.009...% O( J$ m$ ? q3 x! P3 {
&abs(topo_raw(:,2)-lat(i,j))>=0 & abs(topo_raw(:,2)-lat(i,j))<0.009);# \+ T& H9 x: t+ p l& T
topo(i,j)=topo_raw(idx,3);
* F. S& k3 o# O4 d1 p3 O* z. e end4 G, `; N# N, w8 N* H: W
end
: e1 H: m9 E, ?" P最终用到的数据变量为:lon, lat, topo.
( J2 { h0 i0 w- _% O/ J+ m, w
9 ~9 o9 X" T# z& z3 G: l
}- T4 t8 W0 {% H% ^
% @+ Y3 h3 `% Q% h( |- i9 l3 y! u v; w3 u3 Q" E D6 H
转载本文请联系原作者获取授权,同时请注明本文来自叶瑞杰科学网博客。 |