以編程方式提取澳大利亞BOM天氣數據與R

Question

這裏http://www.bom.gov.au/climate/data/我可以輸入一個變電站號碼，比如009572;選擇變量（比如溫度）及其類型（比如說最大值）。點擊「獲取數據」會將我帶到一個鏈接爲「所有年份的數據」的頁面。點擊它，你有一個zip文件。我知道this questions，但在這裏我沒有直接鏈接到一個zip文件。可以採用R來自澳大利亞氣象局網站的天氣數據自動提取？

Answer 1

下面是我做即時下載代碼，同時也解決了您的P_C問題。如果你想和發佈，你可以改進功能。

#daily code = 136 
#monthy code = 139 

bomdata<- function(station,code){ 
for(i in 1: length(station)){ 
p.url<-paste("http://www.bom.gov.au/jsp/ncc/cdio/weatherData/av?p_stn_num=",station[i],"&p_display_type=availableYears&p_nccObsCode=",code,sep ="") 
download.file(p.url,"test.txt") 
filelist <- list.files(pattern = ".txt") 
foo<- file(filelist,"r") 
text<- suppressWarnings(readLines(foo)) 
close(foo) 
l<- regexpr(":",text[1]) 
m<- unlist(gregexpr(",", text[1], perl = TRUE)) 
pc<- substr(text[1],l[[1]]+1,l[[1]]+(m[2]-(l[[1]]+1))) 
url<-paste("http://www.bom.gov.au/jsp/ncc/cdio/weatherData/av?p_display_type=dailyZippedDataFile&p_stn_num=",station[i],"&p_c=",pc,"&p_nccObsCode=",code,"&p_startYear=2013", sep ="") 
suppressWarnings(download.file(url,paste(station[i],".zip",sep= ""), mode = "wb")) 
unlink("test.txt") 
} 
} 

例

bomdata(073137,136) 

Answer 2

雖然我仍然無法看到如何使用download.file（）來做到這一點，但下面幾乎提供了Chrome的「詢問下載之前保存每個文件的位置」這項工作。

system(paste('"C:/Documents and Settings/UserName/Local Settings/Application Data/Google/Chrome/Application/chrome.exe"', 
     '-url http://www.bom.gov.au/jsp/ncc/cdio/weatherData/av?p_display_type=dailyZippedDataFile&p_stn_num=009572&p_c=-18465084&p_nccObsCode=136'), wait = FALSE) 

然後，我可以通過不同的站號使用paste0（）和循環，如果我知道P_C = -18465084手段，它從站如何變化，以車站。

Answer 3

你可以試試這個，它是metvurst包使用的代碼序列。 metvurst

## SET URL FOR DATA DOWNLOAD 
url <- "http://www.bom.gov.au/ntc/IDO70004/IDO70004_" 

## YEARS TO BE DOWNLOADED 
yr <- 1993:2012 

## READ DATA FOR ALL YEARS FROM URL INTO LIST 
fijilst <- lapply(seq(yr), function(i) { 
read.csv(paste(url, yr[i], ".csv", sep = ""), na.strings = c(-9999, 999)) 
})