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

1.进入网站www.52ocean.cn
0 x! u, n7 q% K! d

0 B1 k9 G4 u; S; K& {8 e# g- m                               
登录/注册后可看大图

2.点击进入数据下载界面

, J& r& d( {# R/ G. B+ [

                               
登录/注册后可看大图

7 S8 S1 q' Z6 M8 k. @0 r% {3.设置需要提取的地形数据经纬度范围，由于数据精度比较大， 大范围的数据需要分块进行下载，设置好之后，点击“get data”

$ |" E6 g0 {( W. w, T

6 X) L2 n3 q2 d; O                               
登录/注册后可看大图

7 j/ g( X! h0 \  a0 B# i6 @! k
4. 将数据拷贝到文本中，例如‘.txt’
* a  T* A1 M: l3 k/ D$ `【附】Matlab对上述数据进行网格化处理的程序：
clear all;clc
load('E:\Dataprocesstopodatatopo_raw.txt')
( r2 P# F$ {# Z- M( x+ {2 cLon=sort(unique(topo_raw(:,1)),'ascend');
Lat=sort(unique(topo_raw(:,2)),'ascend');
lon=repmat(Lon',[length(Lat),1]);
( r4 Q/ p# T2 m/ _& j4 C7 ~& K- _lat=repmat(Lat,[1,length(Lon)]);
' G9 Y5 O  g4 a7 N5 yfor i=1:length(Lat)
   for j=1:length(Lon)
       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);
       topo(i,j)=topo_raw(idx,3);  
- y; ~/ p- [: D" y. P   end
end
* p. c% K) z7 O$ I最终用到的数据变量为：lon, lat, topo.
) C6 N$ W4 U! N+ j) k8 \: u# j" v

; v- Q% Z, `! j' G8 E% _                    

2 J0 o- G8 S: Z7 z                                        转载本文请联系原作者获取授权，同时请注明本文来自叶瑞杰科学网博客。