1.进入网站www.52ocean.cn+ P3 p3 L* ?- p
; m& {5 t+ H Q& F# ~
2.点击进入数据下载界面9 u d9 A# A8 V* o3 F$ v
4 v- l; @" s `% N0 P% P! w2 \. ?; u5 d$ V; k2 ?
3.设置需要提取的地形数据经纬度范围，由于数据精度比较大， 大范围的数据需要分块进行下载，设置好之后，点击“get data”) q# o7 d! l+ r4 Y- d
0 Z5 `, K# u1 M. D: A4 ?
, d4 Q4 i' z" h6 b6 g) q6 }9 p& @; v- r
4. 将数据拷贝到文本中，例如‘.txt’
6 i7 A1 S; ^# Q- Q【附】Matlab对上述数据进行网格化处理的程序：
; C. z# }$ C! e& K( a3 ?clear all;clc
- x. b# l% M3 `. a' Cload('E:\Dataprocesstopodatatopo_raw.txt')
* g, A1 _$ |3 j' D8 Z1 n% }2 XLon=sort(unique(topo_raw(:,1)),'ascend');
! X5 z' d) I7 ^9 W- x" N5 CLat=sort(unique(topo_raw(:,2)),'ascend');
! q! n9 w$ n8 a7 l2 L* wlon=repmat(Lon',[length(Lat),1]);
0 _; H3 D& r T9 Q" y1 _lat=repmat(Lat,[1,length(Lon)]);9 f2 K- G N5 A. H; e' p
for i=1:length(Lat)
( B+ q, V8 |6 U$ Y3 |9 v, [ for j=1:length(Lon)$ O& _2 \5 | ^# o+ _8 Z
idx = find(abs(topo_raw(:,1)-lon(i,j))>=0 & abs(topo_raw(:,1)-lon(i,j))<0.009...
/ ?- W- b+ m% D; g A$ I t% D &abs(topo_raw(:,2)-lat(i,j))>=0 & abs(topo_raw(:,2)-lat(i,j))<0.009);
" c2 D, V6 F) W# E& G2 l! K topo(i,j)=topo_raw(idx,3);
' a7 J, z- S- W. H) W end
' _, x" ^0 N5 Pend
$ G* u4 |$ f& d9 E$ e* s最终用到的数据变量为：lon, lat, topo.
! B t+ B) s& z3 d# l4 M) P2 }4 B% K: a7 m6 O( l$ d
% ~3 y8 r8 E6 j7 a k
4 k3 N) Q+ [6 u6 R
8 F% V0 p( {9 e' w" r5 O4 Z, `' I! A; {2 E 转载本文请联系原作者获取授权，同时请注明本文来自叶瑞杰科学网博客。 |
0 D) W" M2 ~3 \* C& { J
