1.进入网站www.52ocean.cn
7 i4 A( {9 k/ N6 f. z
3 s i9 H& s C8 L, U0 t2.点击进入数据下载界面
/ R+ \8 v3 W) ]6 V7 B4 ]. t% u: W) @: M
7 @$ g7 R9 ]' n I$ Y! A3.设置需要提取的地形数据经纬度范围，由于数据精度比较大， 大范围的数据需要分块进行下载，设置好之后，点击“get data”3 A- H. k2 J* n$ ?! T" i
9 z! H- h- e0 q, x2 p- \5 G
8 D) s1 p9 W5 z7 M/ h$ D2 s5 @5 j, M. L
- {. z) N1 W) c6 n4. 将数据拷贝到文本中，例如‘.txt’2 Z+ M7 o5 a! v# z# u5 F
【附】Matlab对上述数据进行网格化处理的程序：
; |( @" c5 R8 W& T; `, J8 Mclear all;clc
) x# a" t# U/ Fload('E:\Dataprocesstopodatatopo_raw.txt'), K+ w- ?. q/ \, B9 q
Lon=sort(unique(topo_raw(:,1)),'ascend');9 c, ^% L. [6 n& {# H/ ?5 m5 i
Lat=sort(unique(topo_raw(:,2)),'ascend');& R! E L8 d$ k5 I8 ?8 [/ e
lon=repmat(Lon',[length(Lat),1]);; d1 e) h& T4 e
lat=repmat(Lat,[1,length(Lon)]);
1 N3 t) X$ d' bfor i=1:length(Lat)4 E1 ^' R+ W3 L( x* } w4 c( d% K
for j=1:length(Lon)/ ~5 \% J! Q! `, j
idx = find(abs(topo_raw(:,1)-lon(i,j))>=0 & abs(topo_raw(:,1)-lon(i,j))<0.009...' Z* V. m" ]9 A0 }9 V* j3 i
&abs(topo_raw(:,2)-lat(i,j))>=0 & abs(topo_raw(:,2)-lat(i,j))<0.009);
& [$ x" C9 C1 P2 H topo(i,j)=topo_raw(idx,3); % o Y. Z2 K s7 a$ c/ O4 o: q( ^; X# X
end/ s7 Y' h, s( K( T, }* L" B
end
6 d9 v2 e. k, X最终用到的数据变量为：lon, lat, topo.
( Y4 u7 W' x! Z# Y# f# m/ o
: q0 G9 Y# x2 F0 w( R" j" |
+ M. @ } N3 k5 O3 Q( `
! N8 w! ]0 J( L8 p( ] h1 d8 ]4 n* O
转载本文请联系原作者获取授权，同时请注明本文来自叶瑞杰科学网博客。 |
: N. N* P5 ^8 U% S" f: S4 @
% M$ ?) n" h) `