1.进入网站www.52ocean.cn
& G% R! o! [0 l2 t: P1 p+ A$ _' O$ [+ Y# S% U
2.点击进入数据下载界面) c7 `" y/ \ j! b, _2 o& a
! i1 ?0 i& h/ ~" r; B( J4 M
! a- X, D! d3 L* I* p. f; O0 j; W& M3.设置需要提取的地形数据经纬度范围,由于数据精度比较大, 大范围的数据需要分块进行下载,设置好之后,点击“get data”
( V- k0 u6 I1 x; V/ c1 J( T- ^: @4 ^- `; q- W
, D# B$ h5 ~% n3 L7 f8 u$ t4 j3 y1 h9 U
4. 将数据拷贝到文本中,例如‘.txt’
0 h0 }. T5 u: R/ y/ G$ X( y3 I【附】Matlab对上述数据进行网格化处理的程序:
" @1 N# p, V/ I/ j9 C' mclear all;clc0 F d6 n) R7 J% _" j
load('E:\Dataprocesstopodatatopo_raw.txt')( Y, _2 y9 R, W' S8 Y
Lon=sort(unique(topo_raw(:,1)),'ascend');
) H( @- K$ D! C1 L5 |9 r8 BLat=sort(unique(topo_raw(:,2)),'ascend');! `2 {2 X5 ~& e9 C1 V
lon=repmat(Lon',[length(Lat),1]);3 j/ c& R/ k5 i: k7 Q
lat=repmat(Lat,[1,length(Lon)]);
) h) W3 g" Z/ f+ ~$ p. k7 Cfor i=1:length(Lat) `/ N" g' W9 \8 G. J! p- Z) |
for j=1:length(Lon)
7 L+ ^/ Q9 x" f4 | y$ H8 ~ idx = find(abs(topo_raw(:,1)-lon(i,j))>=0 & abs(topo_raw(:,1)-lon(i,j))<0.009...5 }1 s6 q4 S- I+ |
&abs(topo_raw(:,2)-lat(i,j))>=0 & abs(topo_raw(:,2)-lat(i,j))<0.009);5 n3 P& t0 ^( R( I
topo(i,j)=topo_raw(idx,3); , w3 f7 b% \/ w2 D w5 Y- H" y
end. a# d& Y2 d( z; U2 r2 g
end
& }0 B' P$ R$ }: w" {6 z9 P最终用到的数据变量为:lon, lat, topo.
% x+ }/ j) z* C
5 h2 U; `7 m" a0 A
( M$ B& u3 k9 p* s& u9 y
' F2 _3 H/ G9 F3 @, u0 J8 P+ L6 G# J. t- g
转载本文请联系原作者获取授权,同时请注明本文来自叶瑞杰科学网博客。 |