1.进入网站www.52ocean.cn& b, _) X4 G1 e |
/ s* [! S% J) O- f% E @2.点击进入数据下载界面
9 n/ f+ z5 @ @. l# C5 s0 H; V+ b8 r9 C M( g1 J
/ Y, Y6 D# V8 r* u5 K
3.设置需要提取的地形数据经纬度范围，由于数据精度比较大， 大范围的数据需要分块进行下载，设置好之后，点击“get data”8 h) F2 {' D$ E0 b0 g
* z( E+ M0 i! J) _3 t/ v: ^/ [0 |) t I. F8 \3 T
* [% [( u; k. d/ g' ~$ j z4. 将数据拷贝到文本中，例如‘.txt’
/ W! i' W( E* s【附】Matlab对上述数据进行网格化处理的程序：
' S9 k* v3 S( n" p2 B/ Cclear all;clc, z: p: e) n9 Y5 p5 }* X F
load('E:\Dataprocesstopodatatopo_raw.txt')
: R3 C" {* T) e f2 H- mLon=sort(unique(topo_raw(:,1)),'ascend');! }, P# g1 F& p0 g0 P2 X
Lat=sort(unique(topo_raw(:,2)),'ascend');
/ ]; \& f6 p- Xlon=repmat(Lon',[length(Lat),1]);
* y. T! H4 s3 Ylat=repmat(Lat,[1,length(Lon)]);. G" l7 X2 L# A) F
for i=1:length(Lat)+ L2 |& X: t- x1 `2 V5 T, p! D
for j=1:length(Lon)) M2 Q1 ^* M' G* M6 P% y- Y
idx = find(abs(topo_raw(:,1)-lon(i,j))>=0 & abs(topo_raw(:,1)-lon(i,j))<0.009...
* A" ]+ x* ^+ i m/ I3 q& Q &abs(topo_raw(:,2)-lat(i,j))>=0 & abs(topo_raw(:,2)-lat(i,j))<0.009);
$ n6 o2 N, A( x; Q. n$ B: L topo(i,j)=topo_raw(idx,3); , U3 y% G0 z, X8 j5 x; D
end; d2 y6 [. ^7 B! D
end5 ]8 s7 u0 y, j
最终用到的数据变量为：lon, lat, topo.: D7 Z/ a- c' H) h- ^
9 ]) e& ^) [: f- y
) m% N$ c. l! w+ w) L
9 F# s: a6 Z! o# Z6 a( t
5 U4 t( d% e2 Q4 X1 r3 w/ @ 转载本文请联系原作者获取授权，同时请注明本文来自叶瑞杰科学网博客。 |
$ R, o" ?; O a, Z A+ T2 d
8 ]3 g t6 v/ I ^5 O6 p