1.进入网站www.52ocean.cn u$ |, |- @3 V
2 l c* m% q9 d I1 l9 \: w, N2.点击进入数据下载界面) t, R( Q8 b8 K- h; s+ D
& E8 M8 J( C9 p: D
+ ]# w# Q$ t l# |/ b
3.设置需要提取的地形数据经纬度范围,由于数据精度比较大, 大范围的数据需要分块进行下载,设置好之后,点击“get data”0 U+ O1 r+ O- o( G, ]! ^
" |7 j9 |3 P0 t
8 v ?3 h8 Y3 P2 w1 d& k, i
, P! Q& c* M, ?6 S8 Z5 W4. 将数据拷贝到文本中,例如‘.txt’
1 q( s7 S( Z [) e【附】Matlab对上述数据进行网格化处理的程序:
, H4 M+ j x; T& I. _1 y5 K+ {# vclear all;clc5 _$ l$ t) N; _, o3 l# P" G
load('E:\Dataprocesstopodatatopo_raw.txt')1 E) b$ K5 T4 {* ]2 o7 Z% F
Lon=sort(unique(topo_raw(:,1)),'ascend');4 ~0 o3 [' C! q S( m4 z/ C
Lat=sort(unique(topo_raw(:,2)),'ascend');
! q1 l! M# u5 I; E* [: llon=repmat(Lon',[length(Lat),1]); y+ x6 o6 d3 _
lat=repmat(Lat,[1,length(Lon)]);
8 S8 S, d" G; y& S9 u- }; @for i=1:length(Lat)
3 N+ G5 g5 p( j for j=1:length(Lon)$ N" M% Z; U! O, c8 [: [2 @8 J
idx = find(abs(topo_raw(:,1)-lon(i,j))>=0 & abs(topo_raw(:,1)-lon(i,j))<0.009...
- P8 g# T8 [6 i/ x4 ] &abs(topo_raw(:,2)-lat(i,j))>=0 & abs(topo_raw(:,2)-lat(i,j))<0.009);8 h* x5 m6 r$ M# b/ h |
topo(i,j)=topo_raw(idx,3); ) p; |, b. z2 A9 I- u
end7 E( W6 w( h3 |* r5 o0 J, D' Q, w; d
end
+ r. d5 o# E/ O$ f最终用到的数据变量为:lon, lat, topo.
, d3 x" Q+ n: i2 ]0 B$ D
5 H, r8 Q; ]8 C4 v7 i) Z7 J2 P( O# o$ C9 k4 j
& K, [' R" \& n
+ h6 A& \5 Q/ a. P5 l 转载本文请联系原作者获取授权,同时请注明本文来自叶瑞杰科学网博客。 |