1.进入网站www.52ocean.cn
T4 j/ o, n* S6 d9 {# g, c, m5 C! y$ ?. y" {
2.点击进入数据下载界面' ]0 t; B4 u2 {9 Q
: d) J$ ]2 Y I5 e4 S: C# ^% k
Y4 |4 x4 @3 I- i P7 s: ]3.设置需要提取的地形数据经纬度范围,由于数据精度比较大, 大范围的数据需要分块进行下载,设置好之后,点击“get data” g+ D. g4 ~2 { n9 R5 P1 x
. K- i: e& X6 Y8 |$ i. @. Y% f
$ a* N$ i% a; ]; A# C
+ f7 k, L5 I% A5 N: f$ P4. 将数据拷贝到文本中,例如‘.txt’
% S, `: ~6 U, U【附】Matlab对上述数据进行网格化处理的程序:
0 A4 x( \, Y1 e+ o7 O; O4 Sclear all;clc( A9 S0 t6 H( ]+ O" J
load('E:\Dataprocesstopodatatopo_raw.txt')
$ Y& R: }# U8 L. P/ w2 T qLon=sort(unique(topo_raw(:,1)),'ascend');+ Y2 S& f. z$ \6 @1 j
Lat=sort(unique(topo_raw(:,2)),'ascend');
; ]$ D i& b9 Z2 }4 | Blon=repmat(Lon',[length(Lat),1]);
9 J6 S! g$ Z$ x% n5 glat=repmat(Lat,[1,length(Lon)]);3 l! H$ L/ M- J/ F {
for i=1:length(Lat); { g1 \8 S/ \* | v$ T Q g
for j=1:length(Lon)
5 U8 Y7 Z9 O$ r) f idx = find(abs(topo_raw(:,1)-lon(i,j))>=0 & abs(topo_raw(:,1)-lon(i,j))<0.009...
2 w, x% z8 d. Y2 Z) c9 X" r &abs(topo_raw(:,2)-lat(i,j))>=0 & abs(topo_raw(:,2)-lat(i,j))<0.009);% j; b* X( H7 m, y; z! M# X
topo(i,j)=topo_raw(idx,3);
. m' \6 ?/ O4 ^4 ?& N* f* Q0 l. P end% C6 I" t: ?: `# G) T# S% D
end
) c- p8 O- q" ?7 S8 J. y最终用到的数据变量为:lon, lat, topo.8 M. |7 l$ v* `" o/ R: G# E
; ^1 E7 ^/ D7 F: Y
; t$ v; ~+ z9 d+ f( m; r% G% j
" }8 L" ~( G: L+ M- J
4 o1 d. ~! ~" f+ g4 o5 j# H 转载本文请联系原作者获取授权,同时请注明本文来自叶瑞杰科学网博客。 |