1.进入网站www.52ocean.cn
" i; U/ S. v6 @3 c
1 W9 e( b) ]& Q8 T# p& P3 f" D2.点击进入数据下载界面
% |2 ?- B# w4 v3 Y; L- U' c, {$ P# B8 X$ L9 `" C# i! I
! H- e) S# | I. ^4 y
3.设置需要提取的地形数据经纬度范围，由于数据精度比较大， 大范围的数据需要分块进行下载，设置好之后，点击“get data”
7 N+ \, O2 J# J9 n w, j N) a% r5 s! _7 y1 x8 q5 Q3 F
8 T7 }% M1 d" b- N' r. ?' n! E6 B
, \' H( b8 S) d! Q4 T( f4 O6 n4. 将数据拷贝到文本中，例如‘.txt’6 i; Y( U1 d1 s: f* N
【附】Matlab对上述数据进行网格化处理的程序：
% i b+ K& h. Y. @; uclear all;clc& p R4 h# P8 n, O. S) L% u) R1 O
load('E:\Dataprocesstopodatatopo_raw.txt') o" |! Y' x9 [5 E
Lon=sort(unique(topo_raw(:,1)),'ascend');
/ p( T. r. u. ^4 HLat=sort(unique(topo_raw(:,2)),'ascend');
. o( o0 L: K8 l1 V0 Wlon=repmat(Lon',[length(Lat),1]);1 n3 t7 P7 M. `# i# l2 r$ b% T! ?; }
lat=repmat(Lat,[1,length(Lon)]);2 B- U& t- M( K3 c7 Q) N& I. X5 L
for i=1:length(Lat)
5 g* u6 `0 C, z/ r2 ]9 I for j=1:length(Lon)3 r0 w2 V9 ^1 \- w B5 j( J7 n! b
idx = find(abs(topo_raw(:,1)-lon(i,j))>=0 & abs(topo_raw(:,1)-lon(i,j))<0.009...
- i/ D6 J! X2 L7 T. Z$ b2 t f! W Q &abs(topo_raw(:,2)-lat(i,j))>=0 & abs(topo_raw(:,2)-lat(i,j))<0.009);
6 Z& }0 \' o' Y" d; y topo(i,j)=topo_raw(idx,3); 8 w" t. z1 \9 @0 z0 c! F9 u' ?/ x
end
. o1 J( U4 ~( f5 Yend
9 V1 @5 _" R7 a/ v) `( U最终用到的数据变量为：lon, lat, topo.
( V) g6 o# k) N1 T8 f5 d. Q7 d7 m% n8 S3 Q3 g
7 I) a4 `1 u# _4 i6 ]7 {& f& _0 n
: F3 ~% l- @0 k/ z6 x0 }# W2 K% R7 L) Z" ]7 o
转载本文请联系原作者获取授权，同时请注明本文来自叶瑞杰科学网博客。 |
: _, v2 W: W) W, Y, [. V. K/ V