1.进入网站www.52ocean.cn8 J$ E/ X2 Y. v1 D- w! f
( ^! M/ G7 Y& ]& N) i
2.点击进入数据下载界面
- \3 ?+ L6 k# d( e' _0 Z; X* n0 q% ?% b8 H$ g2 w1 f; [5 d) b
6 K7 ^; ^/ _ p4 b( w3.设置需要提取的地形数据经纬度范围,由于数据精度比较大, 大范围的数据需要分块进行下载,设置好之后,点击“get data”
8 U( |/ M9 z8 T3 V5 S5 h/ l7 i5 Z1 p' E
6 i) V% W! x+ S9 G6 w- _4 b
+ H* ^2 x6 n1 C: F
4. 将数据拷贝到文本中,例如‘.txt’
% f% V) `/ y: M- h* C* b【附】Matlab对上述数据进行网格化处理的程序:
4 c6 a7 s$ u# k8 B! j1 j2 r+ ]9 dclear all;clc) t6 o! E4 [" Z. }9 u( J% u
load('E:\Dataprocesstopodatatopo_raw.txt')
1 \0 l5 I( U& ]( D e {Lon=sort(unique(topo_raw(:,1)),'ascend');
$ W0 {2 d; o5 J6 E u$ ~, A: CLat=sort(unique(topo_raw(:,2)),'ascend');
! z9 i3 u# o! e" i* V5 Llon=repmat(Lon',[length(Lat),1]); Q! \) X% {- @8 n. A
lat=repmat(Lat,[1,length(Lon)]);% g" i7 ]; X. g0 H# ?9 H! T: p
for i=1:length(Lat)3 q; n' O7 ]; o5 ?7 H! y
for j=1:length(Lon)
+ `# }% C, j4 i# I- s: K idx = find(abs(topo_raw(:,1)-lon(i,j))>=0 & abs(topo_raw(:,1)-lon(i,j))<0.009...* a( C1 w5 [2 r- B# a0 r2 T# `
&abs(topo_raw(:,2)-lat(i,j))>=0 & abs(topo_raw(:,2)-lat(i,j))<0.009);- n. u$ j/ N# C' F$ K" v
topo(i,j)=topo_raw(idx,3); + e, Q* Q4 W! ^5 h: c
end
/ v0 }" @1 C! rend
# d+ Q: A+ ~3 A& S3 I最终用到的数据变量为:lon, lat, topo.
& p, @% ^! y0 @/ ~; G$ ?$ I% T$ o6 B7 o( k" M }! @0 U. @
5 D2 V" C$ H$ Y, {4 u$ H $ j2 x- n7 _4 h6 b, l
( ~9 H$ N( O/ j
转载本文请联系原作者获取授权,同时请注明本文来自叶瑞杰科学网博客。 |