【数据】全球高精度地形数据下载(Sandwell and Smith地形数据)

1.进入网站www.52ocean.cn
% E$ b  ?1 _+ E0 B

7 A) K9 b  R! A- u- z                               
登录/注册后可看大图

0 j9 k, N4 G, B  K% O7 H" x/ \2.点击进入数据下载界面

# S9 l! l0 N) J& i+ k: w9 j, G

+ }) d2 v' n( a( Y9 w                               
登录/注册后可看大图

2 j- I2 Q1 a2 v: Q* c3.设置需要提取的地形数据经纬度范围，由于数据精度比较大， 大范围的数据需要分块进行下载，设置好之后，点击“get data”
  G2 ?$ q% I' q

3 f9 I) t: @, W6 \4 f# J                               
登录/注册后可看大图

. z% z& C8 V! B
) k/ F- K. @4 f4. 将数据拷贝到文本中，例如‘.txt’
1 O5 j4 r' v0 S, O' L; t【附】Matlab对上述数据进行网格化处理的程序：
clear all;clc
load('E:\Dataprocesstopodatatopo_raw.txt')
Lon=sort(unique(topo_raw(:,1)),'ascend');
; ^! I. P" H8 p+ BLat=sort(unique(topo_raw(:,2)),'ascend');
lon=repmat(Lon',[length(Lat),1]);
lat=repmat(Lat,[1,length(Lon)]);
for i=1:length(Lat)
) n3 o! \# M7 i! w/ @& [   for j=1:length(Lon)
5 T& V) w/ p6 u( g       idx = find(abs(topo_raw(:,1)-lon(i,j))>=0 & abs(topo_raw(:,1)-lon(i,j))<0.009...
                  &abs(topo_raw(:,2)-lat(i,j))>=0 & abs(topo_raw(:,2)-lat(i,j))<0.009);
. g$ ?1 D  K( O8 u( @+ ?0 A       topo(i,j)=topo_raw(idx,3);  
  S, b7 q, M4 D& i' S7 y   end
end
7 o. z7 J. O% t6 d最终用到的数据变量为：lon, lat, topo.
7 n; w$ r# L8 A. X

# j- H3 q' S  T/ L7 I                    

9 s) V' M  k  G8 n  [" @7 s1 e                                        转载本文请联系原作者获取授权，同时请注明本文来自叶瑞杰科学网博客。