1.进入网站www.52ocean.cn" t1 V! J( g V) r
7 D3 p. s" K, R$ ?- l9 `
2.点击进入数据下载界面
% r; |& f1 O% A) t$ ~; f
3 u! k. e0 p; F1 i3 K4 k, n
1 b; f8 F# B1 N3.设置需要提取的地形数据经纬度范围,由于数据精度比较大, 大范围的数据需要分块进行下载,设置好之后,点击“get data” n, K* h2 v: N1 } K
9 I- G% i: _2 U
4 x1 T& Y- T C8 n
0 t+ {8 _3 t* i; s2 i4. 将数据拷贝到文本中,例如‘.txt’( Y' u$ D) ~% E/ X+ T+ E
【附】Matlab对上述数据进行网格化处理的程序:( W I1 Z1 d) m/ H' I# S, Z
clear all;clc: }3 h4 Q" u. E2 `( i
load('E:\Dataprocesstopodatatopo_raw.txt')5 s7 z5 l7 D0 m* B( ~3 H6 k1 ?- b% L
Lon=sort(unique(topo_raw(:,1)),'ascend');
2 c: n# V9 B" V3 c) gLat=sort(unique(topo_raw(:,2)),'ascend');! ] c) w: T2 H4 G# k1 I
lon=repmat(Lon',[length(Lat),1]);
1 i: f: L& u' D$ C' j0 s) R* Vlat=repmat(Lat,[1,length(Lon)]);! x) ~+ L7 q8 _0 J$ @
for i=1:length(Lat)
5 a% ~* m4 l4 q" h/ _. }2 R for j=1:length(Lon)
7 Z, Q1 R; J! j idx = find(abs(topo_raw(:,1)-lon(i,j))>=0 & abs(topo_raw(:,1)-lon(i,j))<0.009.../ N' ]& i9 K( h3 T
&abs(topo_raw(:,2)-lat(i,j))>=0 & abs(topo_raw(:,2)-lat(i,j))<0.009);
3 }1 k. o0 G; H3 N$ \1 ~, h; Z topo(i,j)=topo_raw(idx,3);
& X s% }8 Z5 e+ P end3 P2 U3 H- G8 ~( V% }7 F4 v6 r" O
end/ W/ c0 F( ^2 p% G8 n% I
最终用到的数据变量为:lon, lat, topo.- M$ e# Q+ v g. E
8 M+ v9 ?1 F+ S+ f; p( w$ e$ A: s
; }, a( Z$ f4 d' G! T9 k% n
- v" o3 N& C& ?7 G2 w$ ^8 T3 W/ @# m$ C
转载本文请联系原作者获取授权,同时请注明本文来自叶瑞杰科学网博客。 |