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

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

[复制链接]
1.进入网站www.52ocean.cn- Y% X+ D7 [0 `$ j: \- v6 B
( R3 _1 n4 ^- P$ e- b( G/ \
                               
登录/注册后可看大图
* _' O" c* \: `2 G
2.点击进入数据下载界面
5 X# @2 s" D2 O: v
  z  t+ ^' `: m$ ~

2 k4 ^9 [& m5 ~0 M                               
登录/注册后可看大图

# z% b7 L: i" i- p2 ~: i6 p3.设置需要提取的地形数据经纬度范围,由于数据精度比较大, 大范围的数据需要分块进行下载,设置好之后,点击“get data”
+ C) c; b0 E/ Y' ?0 V/ d2 ?1 b8 D1 a- {) Z+ |

  }& k4 m8 U# j# h* o( g. P5 I                               
登录/注册后可看大图
( W" f2 K2 D2 C0 r

$ i" \7 h6 |9 w- Y! E4 h4. 将数据拷贝到文本中,例如‘.txt’
# m7 W$ ^+ n4 T6 l3 `【附】Matlab对上述数据进行网格化处理的程序:+ ^6 X0 c& d6 w, L
clear all;clc) l0 K# l' @: m" l1 L$ _( r5 j  W+ d
load('E:\Dataprocesstopodatatopo_raw.txt')- n8 F4 K7 K; f5 P! _9 Y8 o! p2 H
Lon=sort(unique(topo_raw(:,1)),'ascend');
) g- m, g5 H, D- C4 Q& z8 N3 ELat=sort(unique(topo_raw(:,2)),'ascend');# `8 |2 V9 e5 Z7 z$ w
lon=repmat(Lon',[length(Lat),1]);. F! _! p" l9 \/ R: j, N
lat=repmat(Lat,[1,length(Lon)]);$ I' g, s5 J6 ]" N& p0 [
for i=1:length(Lat)
+ e' v5 S5 Y0 n  V   for j=1:length(Lon): d* e$ h  `: g% k' I' d
       idx = find(abs(topo_raw(:,1)-lon(i,j))>=0 & abs(topo_raw(:,1)-lon(i,j))<0.009...
" A2 {7 f& r+ L2 ~3 T. q4 ^5 i                  &abs(topo_raw(:,2)-lat(i,j))>=0 & abs(topo_raw(:,2)-lat(i,j))<0.009);
8 y/ z% E; I; Y' {+ a7 o- o1 q7 y       topo(i,j)=topo_raw(idx,3);  * t+ k% X+ A8 F  W7 n
   end0 ?2 V% ~% }, \5 Z
end
# E3 _) J* d$ K3 t% s最终用到的数据变量为:lon, lat, topo.
" m# P9 W" J' d) A2 P: \& c! {6 _' p- g6 ?6 }. S+ u

9 _2 a) Z) U& w# A. I                    1 A5 c( V5 G4 P

6 m2 p1 g4 x- I& x1 S- j  v                                        转载本文请联系原作者获取授权,同时请注明本文来自叶瑞杰科学网博客。
回复

举报 使用道具

相关帖子

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