1.进入网站www.52ocean.cn/ r- B* X8 p* ]
8 G% t. @. i5 D( a) i2.点击进入数据下载界面+ ] G! \* }6 T
+ y% a, |3 a( j6 N4 N G }) I( c; k4 `. y) A* b8 i! X6 e
3.设置需要提取的地形数据经纬度范围，由于数据精度比较大， 大范围的数据需要分块进行下载，设置好之后，点击“get data”
9 A" C$ a( _' M4 c2 q' s
/ z! s8 a% t2 t: A( z, z: q4 \' u- }2 h% X0 T3 F
M4 r1 x3 |0 W# l+ j
4. 将数据拷贝到文本中，例如‘.txt’7 v+ e# f/ d0 F: F- q; x) }
【附】Matlab对上述数据进行网格化处理的程序：
% ?& k h9 `. k- C, tclear all;clc0 T: q2 Q; X+ q5 C# R r' A3 |- c
load('E:\Dataprocesstopodatatopo_raw.txt')
. K+ J4 f8 j0 j3 q9 s2 h1 M9 `Lon=sort(unique(topo_raw(:,1)),'ascend');$ X7 R- Y" b3 _% {) B
Lat=sort(unique(topo_raw(:,2)),'ascend');' o8 [1 }5 z) @
lon=repmat(Lon',[length(Lat),1]);8 Y/ x# L5 t% W9 g4 l5 ^
lat=repmat(Lat,[1,length(Lon)]);5 w& @# m& z* N/ |* N! d
for i=1:length(Lat)
% v8 q8 H4 c( m9 f% B0 s$ [ for j=1:length(Lon)
( Z" O D% a) i1 w) _: f idx = find(abs(topo_raw(:,1)-lon(i,j))>=0 & abs(topo_raw(:,1)-lon(i,j))<0.009...- G5 Z1 ?0 _! ~, I# {! W! W
&abs(topo_raw(:,2)-lat(i,j))>=0 & abs(topo_raw(:,2)-lat(i,j))<0.009); ~, I! C0 }8 b3 E) q$ e: r
topo(i,j)=topo_raw(idx,3); . D4 s2 Z- Q8 W& h7 s$ T
end
8 t( O: Y+ S! l6 ~9 {7 h0 rend
7 ?2 r& _5 C; O0 G8 B0 G" x最终用到的数据变量为：lon, lat, topo.: a) b! t$ o& s8 q/ ^9 W- _
4 k" M% @% u, m2 ^ @
) Q. w6 ?; O5 ]3 A& A9 ]' g # B1 R) s+ J3 X0 A7 G. T: m/ B
. {1 c( j+ V! |. y! q7 D$ l6 g
转载本文请联系原作者获取授权，同时请注明本文来自叶瑞杰科学网博客。 |
1 F1 b; Q# O/ M' N
