1.进入网站www.52ocean.cn& Z, P% _ W6 p& L* v* Q, H; W
' z, ?3 Z% m7 x/ g2 V/ [0 }2.点击进入数据下载界面. _, G8 G# T8 K1 g, I: {1 [
( u; M* \. H. x" r* z4 D1 t k/ Q% A" f) L2 D- C
3.设置需要提取的地形数据经纬度范围,由于数据精度比较大, 大范围的数据需要分块进行下载,设置好之后,点击“get data”
. N' U" _3 h/ h9 Q% u, \3 E. x
+ X' ^. u- h$ C6 s+ ?
) z$ b4 }/ V/ P+ L$ O) |6 a2 P2 S3 F) E9 P6 r/ a
4. 将数据拷贝到文本中,例如‘.txt’
& \8 U. h* A# D, A【附】Matlab对上述数据进行网格化处理的程序:
3 w- R) a! B& g% v4 ]clear all;clc
0 W1 _* h% U3 { g' zload('E:\Dataprocesstopodatatopo_raw.txt')
( G; t7 F; x A7 }- t$ Q# M1 y% YLon=sort(unique(topo_raw(:,1)),'ascend');
8 j! p& k" a- ]0 ~8 qLat=sort(unique(topo_raw(:,2)),'ascend');" m0 X3 a# } c, j
lon=repmat(Lon',[length(Lat),1]);
$ C9 U) `# ]) {) alat=repmat(Lat,[1,length(Lon)]);+ j% ]1 K6 ?- P7 I- h7 h
for i=1:length(Lat)
' H" [2 o8 e4 Z9 L for j=1:length(Lon)! K% x2 a9 _' j# {3 G. P6 V2 E
idx = find(abs(topo_raw(:,1)-lon(i,j))>=0 & abs(topo_raw(:,1)-lon(i,j))<0.009...3 S! }$ Q2 P# d8 g% G; _( y
&abs(topo_raw(:,2)-lat(i,j))>=0 & abs(topo_raw(:,2)-lat(i,j))<0.009);
0 j8 t* u& ?, A3 `* K# E+ M5 F$ s/ N topo(i,j)=topo_raw(idx,3);
v; p: d/ A, Y1 x$ z end
5 ^! f i I5 p2 O! L! Eend* R! ~8 l. b) o! b
最终用到的数据变量为:lon, lat, topo.9 E4 V) \/ ~7 H* `/ o
5 u: f: N# }$ j+ T! i& U6 T$ } J6 d9 D
# p: }# B, i \; D4 O8 f
/ ?0 [1 G6 h0 r0 C 转载本文请联系原作者获取授权,同时请注明本文来自叶瑞杰科学网博客。 |