1.进入网站www.52ocean.cn+ i( q* d, D) c) X9 T
4 n1 ?3 p! _( C
2.点击进入数据下载界面
9 u" T" h; d& B5 M/ F! ^/ H5 n, F( Y; m
5 `" l/ z5 e- z' k0 |/ L6 X. M3.设置需要提取的地形数据经纬度范围，由于数据精度比较大， 大范围的数据需要分块进行下载，设置好之后，点击“get data”1 X* ^) y- n7 \! h
6 }; C( ^; ^& B/ b# K8 M0 H5 h
5 T; K+ q8 T. [" B# n6 R8 {1 P
/ a; N, y7 \, q% i
4. 将数据拷贝到文本中，例如‘.txt’
, [$ |3 s, B* J【附】Matlab对上述数据进行网格化处理的程序：7 x; k$ a% f1 g8 _& h0 L
clear all;clc7 [0 h& H8 x& L7 j- ^
load('E:\Dataprocesstopodatatopo_raw.txt')& q5 N! R0 A F% d4 v6 m+ I
Lon=sort(unique(topo_raw(:,1)),'ascend');- r% j2 b) l& I! f P. d9 k/ z
Lat=sort(unique(topo_raw(:,2)),'ascend');
# |' X. R1 e0 c$ w4 i5 Tlon=repmat(Lon',[length(Lat),1]);9 S4 r/ T c1 v) s& X' P
lat=repmat(Lat,[1,length(Lon)]);) b* r- z' Q0 j' M3 Q5 D1 v
for i=1:length(Lat)
1 C" A$ P+ e- s3 F" S for j=1:length(Lon)
5 i# d- C( H1 ?; d idx = find(abs(topo_raw(:,1)-lon(i,j))>=0 & abs(topo_raw(:,1)-lon(i,j))<0.009...
: m/ o9 }* A- P. l: b; U" O8 [ &abs(topo_raw(:,2)-lat(i,j))>=0 & abs(topo_raw(:,2)-lat(i,j))<0.009);
( { b; n* R8 ^) r- c+ L topo(i,j)=topo_raw(idx,3);
2 j- A' s2 J" _ end2 P- |) \4 [6 a- }
end
) e1 i, X/ [# i" @; K最终用到的数据变量为：lon, lat, topo. Y+ {& x4 w! q4 `0 Z5 o
T$ Z6 }8 |% }. r$ p# `- q+ }4 _
' Q9 S5 @. [, p' C( D# j# w
# J# e+ n9 S: J! R1 Q) _0 \9 o2 P
转载本文请联系原作者获取授权，同时请注明本文来自叶瑞杰科学网博客。 |
' p6 l4 _) d+ i4 b& N3 j
