1.进入网站www.52ocean.cn
0 C7 l. {: {* i. |$ {0 x, y" P' C$ T
8 I; r- B b+ i1 `2 L0 ]- e2.点击进入数据下载界面4 ` n2 J8 w/ a- n J' K( y1 D
: h, \; R8 G- j4 k! I# N2 z
1 }; G# w! G# ?2 ]. l
3.设置需要提取的地形数据经纬度范围,由于数据精度比较大, 大范围的数据需要分块进行下载,设置好之后,点击“get data”1 ^6 {+ o0 @9 k
- `2 L o. H) [) l0 F% a+ Z2 k
) \; R: O# k& u7 \+ K
s6 q' h1 U# i. X6 \* F4 x4. 将数据拷贝到文本中,例如‘.txt’0 } a3 p! y; w
【附】Matlab对上述数据进行网格化处理的程序:
8 q6 v5 d+ o2 L- b0 jclear all;clc
# u: n# c. _$ k" }' yload('E:\Dataprocesstopodatatopo_raw.txt')
/ I4 D0 w' v! z! u: _% oLon=sort(unique(topo_raw(:,1)),'ascend');1 g. Q$ [* v% q, y' ~5 e
Lat=sort(unique(topo_raw(:,2)),'ascend');
; \; ~/ F0 j2 c2 l0 \lon=repmat(Lon',[length(Lat),1]);
* P6 }" ^$ E k6 V# Clat=repmat(Lat,[1,length(Lon)]);9 V: ^. l9 I7 h- {6 I5 x- X
for i=1:length(Lat)
8 a2 a. p! U a( Q/ ]9 K for j=1:length(Lon): N: C/ h+ k U w- i6 @% v T6 x
idx = find(abs(topo_raw(:,1)-lon(i,j))>=0 & abs(topo_raw(:,1)-lon(i,j))<0.009...
( F6 f6 O6 m3 l$ E6 Q &abs(topo_raw(:,2)-lat(i,j))>=0 & abs(topo_raw(:,2)-lat(i,j))<0.009);
. X0 [: h- M' m0 g topo(i,j)=topo_raw(idx,3); 2 Y+ h- e& j6 ^ S/ [/ s
end
4 M( L! A* H. s' t1 L- m! oend" e) \' y# Q2 f: n
最终用到的数据变量为:lon, lat, topo. q- H# r% I4 K( g
! V! I9 H/ V: M' p% y- n6 l+ P! M1 c& W
& f: Q8 ~6 ]* y5 Q; `
9 h9 x/ e( @7 ]9 T. }5 }* b
转载本文请联系原作者获取授权,同时请注明本文来自叶瑞杰科学网博客。 |