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

1.进入网站www.52ocean.cn
. X; I- [0 a: R

" ?+ Q8 s9 J3 `                               
登录/注册后可看大图

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

" w1 M4 W: h4 v

/ Q% T7 h( t3 [7 `( @: b" C                               
登录/注册后可看大图

) {* ]7 Z$ k: D* v3.设置需要提取的地形数据经纬度范围，由于数据精度比较大， 大范围的数据需要分块进行下载，设置好之后，点击“get data”
' ^0 L2 \; E+ a$ X. w
6 r# B$ G$ G  k: U' _) X; @

# X7 `' E: [% V  Y2 A                               
登录/注册后可看大图

4. 将数据拷贝到文本中，例如‘.txt’
: Q0 o0 O' K8 H  l$ {: @8 @【附】Matlab对上述数据进行网格化处理的程序：
clear all;clc
( @( ~( f$ s6 {7 p+ d, _8 ~8 C5 Y0 iload('E:\Dataprocesstopodatatopo_raw.txt')
$ u; r- g" @9 p, p" u  ^Lon=sort(unique(topo_raw(:,1)),'ascend');
Lat=sort(unique(topo_raw(:,2)),'ascend');
- h2 c# ?8 a8 R" [& B; F" ulon=repmat(Lon',[length(Lat),1]);
1 z2 ^7 J) }# e' j, Xlat=repmat(Lat,[1,length(Lon)]);
) ~" C5 G5 S' i' J& k- \, }* Wfor 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);
- B6 i' h3 n' D. X+ p4 G2 T, k       topo(i,j)=topo_raw(idx,3);  
. z& s" H5 I, c: H& g2 ^; s1 i   end
6 q6 _# r3 W. d& Gend
最终用到的数据变量为：lon, lat, topo.


$ f  f- l* G) ?) z9 B5 j/ i                    

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