1.进入网站www.52ocean.cn+ v4 c, L( t2 y# W) ^8 U
. Y1 s1 i8 j# B; Q2.点击进入数据下载界面. U: Q1 O2 d5 j1 x4 ]; `# t1 f7 [
0 V0 ?: s* V* ~! W7 R( D; s
. Q4 @/ u) T6 [# X8 b3.设置需要提取的地形数据经纬度范围,由于数据精度比较大, 大范围的数据需要分块进行下载,设置好之后,点击“get data”: c! U6 x, q0 x% H
, b2 B5 {' D( {( V1 j* r
0 F- O) W/ S2 q* O' W: v1 p& ]1 ~2 c! b. f/ w
4. 将数据拷贝到文本中,例如‘.txt’
6 {9 `- Q8 {1 h8 x# `! i$ w【附】Matlab对上述数据进行网格化处理的程序:9 m" n' g# o% d% d
clear all;clc# `6 V0 F: m% q" e8 P' L& T8 t* d
load('E:\Dataprocesstopodatatopo_raw.txt')
: ^+ Z: O/ j- R! ^5 Z0 ^Lon=sort(unique(topo_raw(:,1)),'ascend');
" |' z1 U8 E6 b6 W }Lat=sort(unique(topo_raw(:,2)),'ascend');
) j: Z' `' p1 e9 Mlon=repmat(Lon',[length(Lat),1]);" p: A1 z1 k$ x! B! Q
lat=repmat(Lat,[1,length(Lon)]);
) y# T( D9 f1 [$ y. d; N: wfor i=1:length(Lat)
, z& a, |( C E6 E$ o; o for j=1:length(Lon)
; u* J' x8 ~+ k+ k5 m l" { idx = find(abs(topo_raw(:,1)-lon(i,j))>=0 & abs(topo_raw(:,1)-lon(i,j))<0.009...# Y" g; Z. ^/ [/ k2 Y* x! U
&abs(topo_raw(:,2)-lat(i,j))>=0 & abs(topo_raw(:,2)-lat(i,j))<0.009);4 f' j3 T U* V" q$ Q3 b$ J
topo(i,j)=topo_raw(idx,3);
, p# O7 [# m3 d2 g9 U( ^ end# o+ A; L6 I0 `9 S6 c7 }
end
2 g _6 F9 ]/ t1 u) d2 m最终用到的数据变量为:lon, lat, topo.
; L7 o. T; f* ]
J" h7 |0 K; l% ~+ V d# C. D- J1 x" I
( r4 N4 a' [& q% }7 r& m/ Z7 W- W0 }8 s
转载本文请联系原作者获取授权,同时请注明本文来自叶瑞杰科学网博客。 |