1.进入网站www.52ocean.cn
9 s E( Q2 M( F! S7 H2 _" k4 U- _9 A' T$ J
2.点击进入数据下载界面
& l# |' \# c d r
$ b# f/ K4 ?, K. _% A) T
; C/ M+ o) H* \7 t3.设置需要提取的地形数据经纬度范围，由于数据精度比较大， 大范围的数据需要分块进行下载，设置好之后，点击“get data”" }$ y% X' K- v! M3 [! Z! \
3 n! n6 X9 I/ [3 v% L" }
% @4 W) Z: A0 R4 ^# s2 q' v
5 I! w I7 o& J* S9 j B t
4. 将数据拷贝到文本中，例如‘.txt’
. O, F f: ~5 G- R- ~, z4 G【附】Matlab对上述数据进行网格化处理的程序：% z# i* V3 ?& t
clear all;clc
) x2 q0 h/ [7 o3 {& Qload('E:\Dataprocesstopodatatopo_raw.txt'). a+ @8 W# h7 h; m7 j
Lon=sort(unique(topo_raw(:,1)),'ascend');. v ?3 \" k: S
Lat=sort(unique(topo_raw(:,2)),'ascend');
" E8 R4 \& A/ f- w5 [* G& I `lon=repmat(Lon',[length(Lat),1]);
) d: k% ?4 M' g i" v b! a! l' Plat=repmat(Lat,[1,length(Lon)]);
' V8 e+ C. X1 B7 s! L6 W3 Cfor i=1:length(Lat)
0 J' P. y5 q5 Z: i for j=1:length(Lon)# y1 i) `1 t% {7 o) q
idx = find(abs(topo_raw(:,1)-lon(i,j))>=0 & abs(topo_raw(:,1)-lon(i,j))<0.009...$ q% ]; h) x! R$ V$ \( i3 ?4 m
&abs(topo_raw(:,2)-lat(i,j))>=0 & abs(topo_raw(:,2)-lat(i,j))<0.009);" ]; ^4 _1 q7 R1 g7 q
topo(i,j)=topo_raw(idx,3); ! G0 ~4 j2 {6 j; w
end: m- o0 w- g$ ~
end* m! g% R5 }! h0 l
最终用到的数据变量为：lon, lat, topo.
. v8 x! @: B5 [' v* M2 a5 D
3 x3 y" p* g! c& T Z. h2 P; J6 C% ]
4 R$ ?4 g% \0 r( S0 Z' d; Y
@& n- L3 g1 e8 ^
1 s7 j. n$ [. C# @) X! K 转载本文请联系原作者获取授权，同时请注明本文来自叶瑞杰科学网博客。 |
+ {# H2 n# `8 ^) T9 s
8 q, K$ x/ ?. h0 [* ^+ R. i! f4 s