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

1.进入网站www.52ocean.cn

2 U- z& q- E9 I) \# F) J                               
登录/注册后可看大图

2.点击进入数据下载界面
- \% q, `/ M$ O

4 b! z8 P6 q# j, o* G7 B2 F' Z+ u                               
登录/注册后可看大图

3.设置需要提取的地形数据经纬度范围，由于数据精度比较大， 大范围的数据需要分块进行下载，设置好之后，点击“get data”
2 w0 |& M2 x% G. `( }# _

                               
登录/注册后可看大图

& @2 Y! x! y+ @% j, j$ k$ t) F! J
4. 将数据拷贝到文本中，例如‘.txt’
0 j7 z+ i1 b: S2 [【附】Matlab对上述数据进行网格化处理的程序：
1 p8 \! J8 v" [; s! L$ ?clear all;clc
' P$ \) X! E9 ]5 v0 Hload('E:\Dataprocesstopodatatopo_raw.txt')
Lon=sort(unique(topo_raw(:,1)),'ascend');
4 |" D9 z( j- f; x; v1 A2 ?+ vLat=sort(unique(topo_raw(:,2)),'ascend');
lon=repmat(Lon',[length(Lat),1]);
lat=repmat(Lat,[1,length(Lon)]);
5 c6 G# G1 ~1 |+ f$ hfor i=1:length(Lat)
; q1 h( s- l8 U0 `6 ^' d  w   for j=1:length(Lon)
; w: E5 Y6 ~& L) V+ t( W2 F       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);  
  ^" R4 e* E$ u5 K  s; S   end
; A# ^) g% A, Z7 f: Q# eend
最终用到的数据变量为：lon, lat, topo.

2 }* `& E2 u; G8 h
                    
2 h! F" P5 P& R1 f* J
. X1 @, n1 y# b' h' F- c7 J, e/ \9 o                                        转载本文请联系原作者获取授权，同时请注明本文来自叶瑞杰科学网博客。