這是我對你要求的嘗試。
站IDS必須有資料集與冒號的前綴,如此GHCND:USC00200230
代替USC00200230
你不可錯過datacategoryid到ncdc
功能。請參閱該文檔的API here
例如,獲得降水數據:
out <- ncdc(datasetid='GHCND', datatypeid = 'PRCP', stationid='GHCND:USC00200230', startdate = "2013-09-03", enddate = "2013-09-30", limit=30)
head(out$data)
station value attributes datatype date
1 GHCND:USC00200230 0 ,,7,1800 PRCP 2013-09-03T00:00:00
2 GHCND:USC00200230 0 ,,7,1800 PRCP 2013-09-04T00:00:00
3 GHCND:USC00200230 0 ,,7,1800 PRCP 2013-09-05T00:00:00
4 GHCND:USC00200230 0 ,,7,1800 PRCP 2013-09-06T00:00:00
5 GHCND:USC00200230 8 ,,7,1800 PRCP 2013-09-07T00:00:00
6 GHCND:USC00200230 0 ,,7,1800 PRCP 2013-09-08T00:00:00
另一個例子,拿最低溫度
out <- ncdc(datasetid='GHCND', datatypeid = 'TMIN', stationid='GHCND:USC00200230', startdate = "2013-09-03", enddate = "2013-09-30", limit=30)
head(out$data)
station value attributes datatype date
1 GHCND:USC00200230 139 ,,7,1800 TMIN 2013-09-03T00:00:00
2 GHCND:USC00200230 128 ,,7,1800 TMIN 2013-09-04T00:00:00
3 GHCND:USC00200230 111 ,,7,1800 TMIN 2013-09-05T00:00:00
4 GHCND:USC00200230 83 ,,7,1800 TMIN 2013-09-06T00:00:00
5 GHCND:USC00200230 139 ,,7,1800 TMIN 2013-09-07T00:00:00
6 GHCND:USC00200230 183 ,,7,1800 TMIN 2013-09-08T00:00:00
如果你離開了datatypeid
你得到的數據對所有數據類型。
請顯示您迄今爲止嘗試創建[可重現的示例]的代碼(http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example)。 – Thomas
請參閱上面的我編輯的編輯,非常感謝 – kellis13
hi @ kellis13我是包的創建者。我們正在努力完成一個更好的版本,以更好的文檔到達CRAN。現在看看你的問題。 – sckott