收藏本站 劰载中...网站公告 | 吾爱海洋论坛交流QQ群:835383472

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

[复制链接]
1.进入网站www.52ocean.cn) D" j+ E. e, G8 k% F

& Y0 {. Q( Q1 n& }! @$ `                               
登录/注册后可看大图
- L; N8 k& k6 ^3 ?8 Y
2.点击进入数据下载界面* \( G6 @- z/ F# t" Z
& Y# H$ \  Q- u7 h  p

( Y7 N4 X; ~9 {% M/ y7 ?) e6 P                               
登录/注册后可看大图
% d' @# ^. t+ K0 D, z0 s
3.设置需要提取的地形数据经纬度范围,由于数据精度比较大, 大范围的数据需要分块进行下载,设置好之后,点击“get data”
3 R( z7 M# d& ~8 @3 j* W0 @5 p
1 V6 r# k. B9 V! i& I& s0 z
4 A) G# Y6 S9 i, \5 i
                               
登录/注册后可看大图
( U9 {/ U! B3 }  {8 }# q, L

5 v; L0 v+ f! b6 {4. 将数据拷贝到文本中,例如‘.txt’( m; j/ }/ |% d/ M! Z! J; g
【附】Matlab对上述数据进行网格化处理的程序:
: Q" z% \* a& B8 `- _clear all;clc$ d) X* `! x; L' ]& x
load('E:\Dataprocesstopodatatopo_raw.txt'); Q- M) s, e  t" K: S
Lon=sort(unique(topo_raw(:,1)),'ascend');9 G! Z! N. P+ i& |6 M' e( q1 x
Lat=sort(unique(topo_raw(:,2)),'ascend');- s& K. ]4 y9 [; F6 v
lon=repmat(Lon',[length(Lat),1]);# t: x! C$ Z9 d% g
lat=repmat(Lat,[1,length(Lon)]);/ M& v$ u" R2 S5 b
for i=1:length(Lat)/ j1 ~& Z& n3 x9 L% r8 v  q; l
   for j=1:length(Lon)
6 J4 J8 }* x! J0 c! |4 L# a- w       idx = find(abs(topo_raw(:,1)-lon(i,j))>=0 & abs(topo_raw(:,1)-lon(i,j))<0.009...; V' d/ w% n3 f$ t0 s2 a
                  &abs(topo_raw(:,2)-lat(i,j))>=0 & abs(topo_raw(:,2)-lat(i,j))<0.009);
3 t: |. W) V0 Q% _7 @0 p       topo(i,j)=topo_raw(idx,3);  
; i" K  W( N$ ~: @, S) @* o! D8 F1 J   end+ N: h% |! t* C% c- k9 |) P
end1 `3 A, X. Z; S% y% I: z
最终用到的数据变量为:lon, lat, topo.
3 L2 \/ u( ^3 Z! Y- C6 F3 S
/ T5 v0 l( m/ D4 Q, P
- i- b+ W/ [# T- a( A                    , k+ T" t$ F4 X4 p4 j) \
4 Q& q% G# d1 `  \5 c: J6 I2 N
                                        转载本文请联系原作者获取授权,同时请注明本文来自叶瑞杰科学网博客。
回复

举报 使用道具

相关帖子

全部回帖
暂无回帖,快来参与回复吧
懒得打字?点击右侧快捷回复 【吾爱海洋论坛发文有奖】
您需要登录后才可以回帖 登录 | 立即注册
董永部落
活跃在2021-12-4
快速回复 返回顶部 返回列表