1.进入网站www.52ocean.cn9 r! r; ]! R1 s' c% u6 \
. @2 r, i; L. W( r: X9 u2.点击进入数据下载界面9 Q9 L3 N0 A5 T/ ]+ `' B/ o- k
0 r Z, U& i7 \& J8 b
+ [9 K) J, [; h6 Q7 L: i3.设置需要提取的地形数据经纬度范围,由于数据精度比较大, 大范围的数据需要分块进行下载,设置好之后,点击“get data”
$ {. ^3 J$ Q) R3 J; j9 L6 C; d$ y( _- `4 R) ^
, i t7 M1 X8 w# t1 M2 \
9 V; P/ a$ `" U2 p# k8 s4. 将数据拷贝到文本中,例如‘.txt’ p9 F7 ^. \4 e
【附】Matlab对上述数据进行网格化处理的程序:" V: n. a6 L, u- G. v W
clear all;clc
* h& t+ Q7 p+ E7 @; Eload('E:\Dataprocesstopodatatopo_raw.txt')( A3 F& I+ `" |
Lon=sort(unique(topo_raw(:,1)),'ascend');
7 H, C0 b' n/ F/ `Lat=sort(unique(topo_raw(:,2)),'ascend');/ ]/ F) ~4 [7 a6 B8 h( J/ d8 r
lon=repmat(Lon',[length(Lat),1]);
/ O( L8 \, d+ H1 t: g& W' E0 {lat=repmat(Lat,[1,length(Lon)]);
0 v/ Y& `. W3 D4 D b6 ]5 bfor i=1:length(Lat)2 ~5 p. O/ _' ~: [) N" N* q+ E
for j=1:length(Lon)
, o5 T5 q# e) Q3 d. m, w idx = find(abs(topo_raw(:,1)-lon(i,j))>=0 & abs(topo_raw(:,1)-lon(i,j))<0.009...
/ d! D4 |0 K. A) ^ &abs(topo_raw(:,2)-lat(i,j))>=0 & abs(topo_raw(:,2)-lat(i,j))<0.009);, b9 S" G7 W6 U( F2 {5 w; i
topo(i,j)=topo_raw(idx,3); # t6 P4 J) V- x% ~, `9 m
end
8 }) C: c6 c# U9 [) i% vend, R. J5 ]" _ X' w
最终用到的数据变量为:lon, lat, topo.
8 h9 u2 |* M1 }0 E) U5 U5 u% _* z0 \/ c' q: R: j
9 L* B) t% n/ E0 M
3 c- [) ^" T) d, i, E0 f& E6 }7 z& n' J8 A9 C( [# n
转载本文请联系原作者获取授权,同时请注明本文来自叶瑞杰科学网博客。 |