1.进入网站www.52ocean.cn
- m0 d1 z' ^% [" {5 j3 b# S& X& e( J: x/ M4 X" {8 b
2.点击进入数据下载界面( d: [: k% [4 B
9 {6 H+ M& K; ^) y( B7 S! f
. T7 w* Q) {* e5 C8 t0 A( }) L/ [) _3.设置需要提取的地形数据经纬度范围,由于数据精度比较大, 大范围的数据需要分块进行下载,设置好之后,点击“get data”3 ]# a1 U/ \! @: J# K! q6 ~
3 q! v. H2 J2 u% ]$ i; s
$ o0 i6 n. k) x$ H7 T, y2 r' v, a" u- [) m
4. 将数据拷贝到文本中,例如‘.txt’! ] ^% E, \9 Q4 ]
【附】Matlab对上述数据进行网格化处理的程序:
. \/ r1 w, {. T% M$ G3 a4 Y$ vclear all;clc! o7 `/ S% k9 T+ a5 f! c
load('E:\Dataprocesstopodatatopo_raw.txt')
( O0 ], a* J# J7 X) gLon=sort(unique(topo_raw(:,1)),'ascend');; t. [4 m* ]9 M2 @& G! L0 J* d# w7 P
Lat=sort(unique(topo_raw(:,2)),'ascend');
+ `* P5 Z" c' c8 \# t3 H7 Y- F7 R: ]( Mlon=repmat(Lon',[length(Lat),1]);+ q1 n& n7 e/ l+ @) r5 h0 {$ v e
lat=repmat(Lat,[1,length(Lon)]);$ e. }( c$ s( s. X0 D( R
for i=1:length(Lat)
- E$ c" }( a. A9 ? for j=1:length(Lon) S2 W+ k1 U% Y! T7 l
idx = find(abs(topo_raw(:,1)-lon(i,j))>=0 & abs(topo_raw(:,1)-lon(i,j))<0.009...
2 m" p9 Z: r* i1 g, C- M &abs(topo_raw(:,2)-lat(i,j))>=0 & abs(topo_raw(:,2)-lat(i,j))<0.009);8 [+ c/ ^- }( @, I' d3 F
topo(i,j)=topo_raw(idx,3);
* ?% t$ m; `. W( N end, t* \2 H1 f6 H( l$ q7 r( K
end# Y/ g3 Y4 q7 `* ^+ q9 s
最终用到的数据变量为:lon, lat, topo.% [$ ?8 ~0 x/ }: v# {
: E$ Z5 o2 }: f; f0 X' M2 g% i- t+ I- N0 ^, K9 s
* Y5 z0 u" [" j5 R. t0 e% B
% }! Q5 F! c! t, o 转载本文请联系原作者获取授权,同时请注明本文来自叶瑞杰科学网博客。 |