1.进入网站www.52ocean.cn
* i& c; c, r5 m5 W+ J+ n- @ ~8 k& ~) ]) X3 |1 v9 _4 c
2.点击进入数据下载界面
" V E2 O( g9 O6 k
( O$ V& n) c+ L) j" \$ H" R8 E! A
" u/ S# a3 \2 y: L4 n3.设置需要提取的地形数据经纬度范围,由于数据精度比较大, 大范围的数据需要分块进行下载,设置好之后,点击“get data”+ s. n$ h1 ~+ V7 A" K
6 ?0 R1 P% h1 z6 d' b$ T) A. G9 u$ Q0 D2 L2 Y7 [) m6 N4 {% X- a4 {
" b" I6 I: ~/ V% M6 O: t
4. 将数据拷贝到文本中,例如‘.txt’6 G4 i! N5 @& h; U) Q
【附】Matlab对上述数据进行网格化处理的程序:% _, g# h, x; P
clear all;clc9 w. X/ q9 w7 M7 l, e
load('E:\Dataprocesstopodatatopo_raw.txt')5 j: w* A! _1 K/ B; F$ [+ m* d# O& E! e: C
Lon=sort(unique(topo_raw(:,1)),'ascend');
. n3 U5 x9 ]- B7 {4 j1 G' ILat=sort(unique(topo_raw(:,2)),'ascend');
4 r. S, b9 [6 z. Dlon=repmat(Lon',[length(Lat),1]);2 L9 O; E. @6 q8 Q" V6 {4 P9 F, \
lat=repmat(Lat,[1,length(Lon)]);) x( m3 [- v: p
for i=1:length(Lat)* F: @0 K( T; G! ?7 v
for j=1:length(Lon)
) u0 D; g% g% K% ~1 T6 [ idx = find(abs(topo_raw(:,1)-lon(i,j))>=0 & abs(topo_raw(:,1)-lon(i,j))<0.009...) A" C8 Y% ]& a5 q1 @% c
&abs(topo_raw(:,2)-lat(i,j))>=0 & abs(topo_raw(:,2)-lat(i,j))<0.009);
) U' O8 O3 \! \9 W) P# R w topo(i,j)=topo_raw(idx,3);
$ U% R, C5 ^. D# }5 C6 n% ~ end
$ M0 W" Z' r- t+ e6 D2 r/ Iend( L3 k: L6 z. g0 G% ?
最终用到的数据变量为:lon, lat, topo.
$ Y- I0 v% T; c0 Y" j1 ?2 u _5 ?4 `& l- c
% V Z* a8 e3 K% j/ ?9 L4 g0 u
2 G' n# o# E8 K; T1 I9 _
5 y. \0 O+ f8 `% L" u7 O 转载本文请联系原作者获取授权,同时请注明本文来自叶瑞杰科学网博客。 |