1.进入网站www.52ocean.cn: _% K5 {) n7 n5 H/ ~/ _
5 ? M& K1 C( ^% r2 T
2.点击进入数据下载界面
4 P4 w O7 T3 G- T9 Y% L; m7 w- [0 X+ r7 e' _8 p; ~- ~
3 U' A+ B" L1 I S
3.设置需要提取的地形数据经纬度范围，由于数据精度比较大， 大范围的数据需要分块进行下载，设置好之后，点击“get data”
. p' r' _+ j4 a
0 t2 g( u% I. ?5 n" c. r1 V& l0 y+ Y u
# @, _" m( F' Q1 M9 q4. 将数据拷贝到文本中，例如‘.txt’# l5 B+ O" e. X6 Q3 D; Q
【附】Matlab对上述数据进行网格化处理的程序：
; {; W, ~: Y0 t* b4 qclear all;clc, {( M; l* b% z8 E6 @( g
load('E:\Dataprocesstopodatatopo_raw.txt')& }; t* D% C! F! Y
Lon=sort(unique(topo_raw(:,1)),'ascend');" ^9 j8 k2 ?. h1 M
Lat=sort(unique(topo_raw(:,2)),'ascend');
8 }. C: t7 j, v( a# R! e$ jlon=repmat(Lon',[length(Lat),1]);
( r9 u' R* _; f' v& klat=repmat(Lat,[1,length(Lon)]);
6 p8 H, _7 R- {1 t3 Y+ dfor i=1:length(Lat)
2 i2 N1 Q! ]: K' h3 R& x, w for j=1:length(Lon)
8 q. i1 ]1 x; y9 I$ c: ^+ V8 V$ t idx = find(abs(topo_raw(:,1)-lon(i,j))>=0 & abs(topo_raw(:,1)-lon(i,j))<0.009...
+ t* r+ o- j- A &abs(topo_raw(:,2)-lat(i,j))>=0 & abs(topo_raw(:,2)-lat(i,j))<0.009);
, H- B- ` F; V" x topo(i,j)=topo_raw(idx,3); . f- D# R; L& m1 k- W7 B
end
& x, V. G) r; e+ e6 M$ D; u3 rend- k8 ^/ p+ S6 n; B! j# k* K s/ |
最终用到的数据变量为：lon, lat, topo.- C3 S" ?* j6 d7 Z
) A$ C. s, \1 ] r" E% l
, R9 D" i. j0 F! T# P/ s
7 K" f; i+ F! H8 d& U. _8 M/ N
8 @/ b4 N1 W4 ^8 K- @/ C 转载本文请联系原作者获取授权，同时请注明本文来自叶瑞杰科学网博客。 |
! \2 F5 s/ q4 H7 X* E