1.进入网站www.52ocean.cn- Y% X+ D7 [0 `$ j: \- v6 B
* _' O" c* \: `2 G
2.点击进入数据下载界面
5 X# @2 s" D2 O: v
z t+ ^' `: m$ ~
# z% b7 L: i" i- p2 ~: i6 p3.设置需要提取的地形数据经纬度范围,由于数据精度比较大, 大范围的数据需要分块进行下载,设置好之后,点击“get data”
+ C) c; b0 E/ Y' ?0 V/ d2 ?1 b8 D1 a- {) Z+ |
( 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 转载本文请联系原作者获取授权,同时请注明本文来自叶瑞杰科学网博客。 |