1.进入网站www.52ocean.cn) D" j+ E. e, G8 k% F
- L; N8 k& k6 ^3 ?8 Y
2.点击进入数据下载界面* \( G6 @- z/ F# t" Z
& Y# H$ \ Q- u7 h p
% d' @# ^. t+ K0 D, z0 s
3.设置需要提取的地形数据经纬度范围，由于数据精度比较大， 大范围的数据需要分块进行下载，设置好之后，点击“get data”
3 R( z7 M# d& ~8 @3 j* W0 @5 p
1 V6 r# k. B9 V! i& I& s0 z( U9 {/ U! B3 } {8 }# q, L
5 v; L0 v+ f! b6 {4. 将数据拷贝到文本中，例如‘.txt’( m; j/ }/ |% d/ M! Z! J; g
【附】Matlab对上述数据进行网格化处理的程序：
: Q" z% \* a& B8 `- _clear all;clc$ d) X* `! x; L' ]& x
load('E:\Dataprocesstopodatatopo_raw.txt'); Q- M) s, e t" K: S
Lon=sort(unique(topo_raw(:,1)),'ascend');9 G! Z! N. P+ i& |6 M' e( q1 x
Lat=sort(unique(topo_raw(:,2)),'ascend');- s& K. ]4 y9 [; F6 v
lon=repmat(Lon',[length(Lat),1]);# t: x! C$ Z9 d% g
lat=repmat(Lat,[1,length(Lon)]);/ M& v$ u" R2 S5 b
for i=1:length(Lat)/ j1 ~& Z& n3 x9 L% r8 v q; l
for j=1:length(Lon)
6 J4 J8 }* x! J0 c! |4 L# a- w idx = find(abs(topo_raw(:,1)-lon(i,j))>=0 & abs(topo_raw(:,1)-lon(i,j))<0.009...; V' d/ w% n3 f$ t0 s2 a
&abs(topo_raw(:,2)-lat(i,j))>=0 & abs(topo_raw(:,2)-lat(i,j))<0.009);
3 t: |. W) V0 Q% _7 @0 p topo(i,j)=topo_raw(idx,3);
; i" K W( N$ ~: @, S) @* o! D8 F1 J end+ N: h% |! t* C% c- k9 |) P
end1 `3 A, X. Z; S% y% I: z
最终用到的数据变量为：lon, lat, topo.
3 L2 \/ u( ^3 Z! Y- C6 F3 S
/ T5 v0 l( m/ D4 Q, P
- i- b+ W/ [# T- a( A , k+ T" t$ F4 X4 p4 j) \
4 Q& q% G# d1 ` \5 c: J6 I2 N
转载本文请联系原作者获取授权，同时请注明本文来自叶瑞杰科学网博客。 |
4 A) G# Y6 S9 i, \5 i