去UHSLC網(wǎng)站上批量下載驗(yàn)潮站數(shù)據(jù)
首先打開(kāi)UHSLC網(wǎng)頁(yè)http://uhslc.soest.hawaii.edu/opendap/rqds/global/hourly/,如下圖所示

可以看到同一個(gè)站點(diǎn)會(huì)有多個(gè)數(shù)據(jù)(以字母a,b,c等來(lái)區(qū)分)。這是因?yàn)樵诤车哪硞€(gè)站點(diǎn)的位置在某個(gè)時(shí)間可能有小范圍的移動(dòng),這種移動(dòng)可能只有幾百米或者幾十米,所以編號(hào)還是一樣的。在大部分情況下,這類(lèi)站點(diǎn)經(jīng)過(guò)校正后可以合并。
我們?cè)谑謩?dòng)下載這些站點(diǎn)數(shù)據(jù)的時(shí)候,打開(kāi)F12開(kāi)發(fā)者工具,就可以捕獲下載鏈接:
url='http://uhslc.soest.hawaii.edu/opendap/rqds/global/hourly/h135a.nc.nc?sea_level%5B0:1:0%5D%5B0:1:15863%5D,time%5B0:1:15863%5D,lat%5B0:1:0%5D,lon%5B0:1:0%5D,station_name%5B0:1:0%5D,station_country%5B0:1:0%5D';
上面的鏈接下載了代號(hào)135a站nc格式數(shù)據(jù)下的sea level, time, lat, lon, station name, station country.
使用如下命令實(shí)現(xiàn)下載135a.nc文件
urlwrite(url,'h135a.nc')
我們只要修改url里的135a,把它變成001a,就可以下載編號(hào)為001a的站點(diǎn)的數(shù)據(jù)文件了。再寫(xiě)一個(gè)循環(huán)即可實(shí)現(xiàn)批量下載!