1.进入网站www.52ocean.cn- H: C4 T; n Q) ]
/ _5 C9 f/ _. X) ?8 ?) U0 l
2.点击进入数据下载界面2 b2 P' E- t) k7 C/ H: l; z+ B$ s) T
4 e# v! b' v. s4 n8 `
& k( a( r' \4 [. P- F. o4 ]
3.设置需要提取的地形数据经纬度范围，由于数据精度比较大， 大范围的数据需要分块进行下载，设置好之后，点击“get data”
1 R1 l& [, K2 g
& F8 D7 O' r( n5 `/ U
; l2 l9 C7 S8 r7 K( J! g0 n6 L3 }# f1 y* ~, p; f4 ~% g2 ]. Z$ c$ i
4. 将数据拷贝到文本中，例如‘.txt’
# ^ m% v6 S- F【附】Matlab对上述数据进行网格化处理的程序：
- Q, Q ?7 R$ `" \. r- z1 w& Oclear all;clc
' p; u* T/ q: \ V+ kload('E:\Dataprocesstopodatatopo_raw.txt')
, z# y, I) }( f: Z- {) C6 [Lon=sort(unique(topo_raw(:,1)),'ascend');
^9 D2 ^2 q4 h* V& l* `3 P) [Lat=sort(unique(topo_raw(:,2)),'ascend');6 O* w9 A o( v5 [5 v
lon=repmat(Lon',[length(Lat),1]);
* C, D+ H% c0 dlat=repmat(Lat,[1,length(Lon)]);0 _" [' t; g, B
for i=1:length(Lat)
1 \' Y8 y. e4 J0 U; B9 u! T for j=1:length(Lon)3 p3 ~3 I- ?2 o5 \# j, E- G0 F
idx = find(abs(topo_raw(:,1)-lon(i,j))>=0 & abs(topo_raw(:,1)-lon(i,j))<0.009..." ~# N- z5 f! |, F9 w* ?
&abs(topo_raw(:,2)-lat(i,j))>=0 & abs(topo_raw(:,2)-lat(i,j))<0.009);8 |; s0 f4 ]) G* y9 R0 B
topo(i,j)=topo_raw(idx,3);
( f- m/ o) o! Z6 Q! H/ V( c end S: K+ t: u, c5 F
end+ n5 g& `# `0 }1 T0 f( ^2 y
最终用到的数据变量为：lon, lat, topo.- L; |9 ^! z. o
# p- ~/ v1 }- F, v4 C# o3 L
) u& U% _2 A V5 t; t
& W; K; \. _ I; U! `+ M$ Z' K' S4 {/ p/ V' J
转载本文请联系原作者获取授权，同时请注明本文来自叶瑞杰科学网博客。 |
