1.进入网站www.52ocean.cn
( s: a4 V% r4 }0 ~* T5 e8 p7 N+ L0 l! y9 ?8 ~+ `6 ^
2.点击进入数据下载界面. t1 O5 b: l, ~& F" ?) r
3 \1 i6 } I9 O4 x4 Y$ r5 h& `" [
) G# ?- @( S5 m; I3 N& }3.设置需要提取的地形数据经纬度范围，由于数据精度比较大， 大范围的数据需要分块进行下载，设置好之后，点击“get data”
2 y" `9 Q, K: m% u, f
9 W6 P8 `; m4 {# X
9 T( k# @: v2 y, K* T% V+ s
; S9 h. \6 C" C3 U4. 将数据拷贝到文本中，例如‘.txt’
7 }$ i1 W% D5 Z. k【附】Matlab对上述数据进行网格化处理的程序：
2 j4 [/ I! S; c+ `clear all;clc
% n* T! |5 |# Kload('E:\Dataprocesstopodatatopo_raw.txt')
$ {8 z! x; q( a" GLon=sort(unique(topo_raw(:,1)),'ascend');, J" K B* M' v1 N! `9 {
Lat=sort(unique(topo_raw(:,2)),'ascend'); Z2 e/ o" Q. f
lon=repmat(Lon',[length(Lat),1]);
4 h' s0 w0 G( E! F/ z1 tlat=repmat(Lat,[1,length(Lon)]);
' {9 ~0 U& X, }) u/ y* `$ t6 nfor i=1:length(Lat)4 H' Z) w" c; y) r4 Y) C% y! g: B
for j=1:length(Lon)' ?# T( K# L+ z8 k' z2 {
idx = find(abs(topo_raw(:,1)-lon(i,j))>=0 & abs(topo_raw(:,1)-lon(i,j))<0.009...
7 R* G5 f( b* d$ n &abs(topo_raw(:,2)-lat(i,j))>=0 & abs(topo_raw(:,2)-lat(i,j))<0.009);8 I+ V) t# O6 @
topo(i,j)=topo_raw(idx,3); : u' B2 S1 [6 ^ T# m* u. }, q
end
" O n) l) I5 U2 U6 Q8 J( send
7 o. I/ x- z8 R2 a: ^3 i最终用到的数据变量为：lon, lat, topo.
0 f R" |4 u! `% K2 Y+ U5 ^: {4 w4 d$ W9 W6 B
6 F ?' X' d- s' J; }3 F: f8 h
0 r. [7 ]8 s. Q8 y
& j5 u) v. `4 c' f) R/ x4 V1 F0 f6 H 转载本文请联系原作者获取授权，同时请注明本文来自叶瑞杰科学网博客。 |
' r# o$ f T8 ]* D
0 n( F% I4 L2 ^; Y3 z6 ?/ E' _
- A) r+ S4 n/ X