1.进入网站www.52ocean.cn4 X, c0 K" n6 ^7 m' G" z2 l: P8 ~
5 w% e3 o" _' ^8 W7 k2.点击进入数据下载界面* m! l2 Z W3 \5 B' q& }! @
, [8 f/ f, Y7 R
- L7 b" w! w# Y1 s, m: u) Y; @3.设置需要提取的地形数据经纬度范围,由于数据精度比较大, 大范围的数据需要分块进行下载,设置好之后,点击“get data”
" R0 ^+ p% m0 l n& U$ `$ @ U ]7 A" `& @8 S9 m" { {
0 l8 P& ?3 K9 o" \$ Y
* @' w; f1 n$ O: Z6 m4. 将数据拷贝到文本中,例如‘.txt’1 T( Z! `0 }& E( k
【附】Matlab对上述数据进行网格化处理的程序:
- L6 T5 K# w% }3 x# }clear all;clc
7 P( B' C! i4 t/ `5 ?+ Gload('E:\Dataprocesstopodatatopo_raw.txt')
0 j0 a: @" D! X( T/ G G( d4 KLon=sort(unique(topo_raw(:,1)),'ascend'); \) ?+ [8 J+ L/ |+ X' f6 }
Lat=sort(unique(topo_raw(:,2)),'ascend');3 t1 s* ~1 X5 `% k! Z& R8 V
lon=repmat(Lon',[length(Lat),1]);
2 v6 X7 h6 Y4 [9 u" W2 ?lat=repmat(Lat,[1,length(Lon)]);% E2 ^* b( g' T. L) c, o2 T
for i=1:length(Lat)5 @. s1 t( j) G9 Z! v1 Z
for j=1:length(Lon)
7 s* {: ?, T. X8 S. w! U; l idx = find(abs(topo_raw(:,1)-lon(i,j))>=0 & abs(topo_raw(:,1)-lon(i,j))<0.009... \: m5 W/ b% D' ]
&abs(topo_raw(:,2)-lat(i,j))>=0 & abs(topo_raw(:,2)-lat(i,j))<0.009);
! L+ F9 ?/ y" l1 L8 T: S+ D; ]$ B topo(i,j)=topo_raw(idx,3); 2 U2 } n8 T; G5 o! C2 y
end
7 b& M4 n* ~% x5 Q2 J2 O! Jend5 R! K4 a, Q; U6 U
最终用到的数据变量为:lon, lat, topo.: q; Q4 |( f- F* N2 u9 O5 W2 @; N3 _
2 s5 M. k& x4 `4 V# q- ~& M; f, k
4 G$ g8 o8 G: _) h; d
$ @: B- `! D6 A# J9 B/ f2 q- M
% h; W5 |' V0 {5 S5 y 转载本文请联系原作者获取授权,同时请注明本文来自叶瑞杰科学网博客。 |