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

1.进入网站www.52ocean.cn
0 C7 l. {: {* i. |$ {0 x, y" P' C$ T

                               
登录/注册后可看大图

8 I; r- B  b+ i1 `2 L0 ]- e2.点击进入数据下载界面

2 x' ^9 h) y; N9 C4 k2 N( N! N* \                               
登录/注册后可看大图

3.设置需要提取的地形数据经纬度范围，由于数据精度比较大， 大范围的数据需要分块进行下载，设置好之后，点击“get data”

2 N8 M& j# E/ b. L                               
登录/注册后可看大图

) \; R: O# k& u7 \+ K
  s6 q' h1 U# i. X6 \* F4 x4. 将数据拷贝到文本中，例如‘.txt’
【附】Matlab对上述数据进行网格化处理的程序：
8 q6 v5 d+ o2 L- b0 jclear all;clc
# u: n# c. _$ k" }' yload('E:\Dataprocesstopodatatopo_raw.txt')
/ I4 D0 w' v! z! u: _% oLon=sort(unique(topo_raw(:,1)),'ascend');
Lat=sort(unique(topo_raw(:,2)),'ascend');
; \; ~/ F0 j2 c2 l0 \lon=repmat(Lon',[length(Lat),1]);
* P6 }" ^$ E  k6 V# Clat=repmat(Lat,[1,length(Lon)]);
for i=1:length(Lat)
8 a2 a. p! U  a( Q/ ]9 K   for j=1:length(Lon)
       idx = find(abs(topo_raw(:,1)-lon(i,j))>=0 & abs(topo_raw(:,1)-lon(i,j))<0.009...
( F6 f6 O6 m3 l$ E6 Q                  &abs(topo_raw(:,2)-lat(i,j))>=0 & abs(topo_raw(:,2)-lat(i,j))<0.009);
. X0 [: h- M' m0 g       topo(i,j)=topo_raw(idx,3);  
   end
4 M( L! A* H. s' t1 L- m! oend
最终用到的数据变量为：lon, lat, topo.

! V! I9 H/ V: M' p% y- n
                    

                                        转载本文请联系原作者获取授权，同时请注明本文来自叶瑞杰科学网博客。