1.进入网站www.52ocean.cn
' R" E! n' n2 S/ Q8 m
# J& G0 c5 C7 k, A2.点击进入数据下载界面9 y: }$ Q- V2 c( Q& H3 T! `2 l/ |
9 `9 q5 I) J0 K! B/ }4 q& b
* h( w8 a" |: C& K8 D) p3.设置需要提取的地形数据经纬度范围,由于数据精度比较大, 大范围的数据需要分块进行下载,设置好之后,点击“get data”
3 e2 F% y% K; g: h. }5 w1 g% k9 {, d1 r% T/ |
" K- m+ i5 N: m! F
. L+ R- F R+ ?: b2 i6 G$ l1 r4. 将数据拷贝到文本中,例如‘.txt’
2 ]0 W3 u& {3 u6 B1 G8 k【附】Matlab对上述数据进行网格化处理的程序:) {- q. s8 r; {# L8 q6 I
clear all;clc
5 S# Y) G, {( |# K0 o& dload('E:\Dataprocesstopodatatopo_raw.txt')
# ^) M0 I! `+ D4 V& lLon=sort(unique(topo_raw(:,1)),'ascend');$ h$ g( I! C' u+ D( I
Lat=sort(unique(topo_raw(:,2)),'ascend');
6 T' s/ n$ r- k: Y5 j0 G: mlon=repmat(Lon',[length(Lat),1]);
* r1 A8 U0 z, ^; D* ~lat=repmat(Lat,[1,length(Lon)]);! Q; U S6 r3 h; P- ]% n
for i=1:length(Lat)" S" @: }: A0 Q4 S
for j=1:length(Lon)0 A: [7 r( }0 o! D; O( s3 M& u0 u, j
idx = find(abs(topo_raw(:,1)-lon(i,j))>=0 & abs(topo_raw(:,1)-lon(i,j))<0.009...
, z# x- z, o* ~2 E. r; R* S &abs(topo_raw(:,2)-lat(i,j))>=0 & abs(topo_raw(:,2)-lat(i,j))<0.009);& h2 L% a+ U* Q, E6 x
topo(i,j)=topo_raw(idx,3);
6 u3 c* H8 }' x: Z: }) X end& ~. _ a4 g9 k
end
* v3 E K. A( R1 l$ X& e G9 Q最终用到的数据变量为:lon, lat, topo.
2 {1 ]+ j( `: l% v( ?1 U) A$ P
' r3 ~- T9 n" J3 K
: U. M1 a/ b* V$ F: R( W8 ^ ' k& y9 V! F2 |- V) E
% @0 z; {, v: J# a- h6 Y# W5 T' t 转载本文请联系原作者获取授权,同时请注明本文来自叶瑞杰科学网博客。 |