1.进入网站www.52ocean.cn# H9 E2 u }) [
2 F6 X H S: d) m$ N& x
2.点击进入数据下载界面
* `& C% y$ |( l) t: m/ [
3 K/ I- `; R1 z( ?. M% h" a P* J0 m
" a5 [( a# O/ b6 O- D7 D+ i9 A0 w0 S( q3.设置需要提取的地形数据经纬度范围,由于数据精度比较大, 大范围的数据需要分块进行下载,设置好之后,点击“get data”
8 w; A9 V7 Z; R' O- k8 q B1 V8 E" |/ l
Y1 h# A9 H8 O) g8 a+ Z
; n2 B! J! U+ @( _4. 将数据拷贝到文本中,例如‘.txt’2 m) w f0 p8 n% ^! i' U6 {
【附】Matlab对上述数据进行网格化处理的程序:
( a# l4 M: P+ g" x( j' Iclear all;clc4 ~+ z& o# J/ K# Z: @: S! B& ? R6 J
load('E:\Dataprocesstopodatatopo_raw.txt')
* B) {: R8 r9 J! f4 }Lon=sort(unique(topo_raw(:,1)),'ascend');8 M! g* |; R1 Y/ v- X0 S/ o# C
Lat=sort(unique(topo_raw(:,2)),'ascend');
% V3 Z4 K; }2 |9 U& A6 J# klon=repmat(Lon',[length(Lat),1]);4 d0 n) `- \7 b. j# D
lat=repmat(Lat,[1,length(Lon)]);
% Y4 i5 P0 z8 i+ Lfor i=1:length(Lat)
9 x4 V% c0 {4 @+ c for j=1:length(Lon)6 [5 x1 Z T1 D# q2 Q+ i2 `
idx = find(abs(topo_raw(:,1)-lon(i,j))>=0 & abs(topo_raw(:,1)-lon(i,j))<0.009...
9 @0 ~+ C# F Y1 l9 _; m, @ &abs(topo_raw(:,2)-lat(i,j))>=0 & abs(topo_raw(:,2)-lat(i,j))<0.009);
! v) F3 F$ t' }6 r topo(i,j)=topo_raw(idx,3);
6 T8 o) T m" y$ v8 ]4 j end0 Q/ }6 R* s, R& @0 ~; H7 R
end
7 m0 @+ y1 h; G2 K2 s; M最终用到的数据变量为:lon, lat, topo.
' I; l6 _- \; E1 d$ U1 D/ R7 Z2 S! J- k% p4 l
' F6 e N0 P( g- Y
4 q, {! {9 i( s
R8 c* c, ?% B i5 }; x 转载本文请联系原作者获取授权,同时请注明本文来自叶瑞杰科学网博客。 |