Kann jemand TRMM 3B42 Zeitreihendaten herunterladen?

Question

Ich versuche TRMM 3B42 3-Stunden-Binärdaten für eine bestimmte Zeitspanne von diesem NASA FTP server herunterladen.

Es ist ein ausgezeichneter Code gemacht von Florian Detsch das tägliche Produkt zum Download (hier ist der Link: https://github.com/environmentalinformatics-marburg/Rsenal/blob/master/R/downloadTRMM.R) im GitHub-only Rsenal Paket. Leider funktioniert es nicht für die 3-Stunden-Daten.

änderte ich den Code:

downloadTRMM <- function(begin, end, dsn = ".", format = "%Y-%m-%d.%H") { 

    ## transform 'begin' and 'end' to 'Date' object if necessary 
    if (!class(begin) == "Date") 
    begin <- as.Date(begin, format = format) 

    if (!class(end) == "Date") 
    end <- as.Date(end, format = format) 

    ## trmm ftp server 
    ch_url <-"ftp://disc2.nascom.nasa.gov/data/TRMM/Gridded/3B42_V7/" 

    ## loop over daily sequence 
    ls_fls_out <- lapply(seq(begin, end, 1), function(i) { 

    # year and julian day (name of the corresponding folder) 
    tmp_ch_yr <- strftime(i, format = "%Y%m") 
    #tmp_ch_dy <- strftime(i, format = "%j") 

    # trmm date format 
    tmp_dt <- strftime(i+1, format = "%Y%m%d.%H") 

    # list files available on server 
    tmp_ch_url <- paste(ch_url, tmp_ch_yr, "", sep = "/") 

    tmp_ch_fls <- tmp_ch_fls_out <- character(2L) 
    for (j in 1:2) { 
     tmp_ch_fls[j] <- paste0("3B42.", tmp_dt, "z.7.precipitation", 
           ifelse(j == 1, ".bin")) 

     tmp_ch_fls[j] <- paste(tmp_ch_url, tmp_ch_fls[j], sep = "/") 
     tmp_ch_fls_out[j] <- paste(dsn, basename(tmp_ch_fls[j]), sep = "/") 

     download.file(tmp_ch_fls[j], tmp_ch_fls_out[j], mode = "wb") 
    } 

    # return data frame with *.bin and *.xml filenames 
    tmp_id_xml <- grep("xml", tmp_ch_fls_out) 
    data.frame(bin = tmp_ch_fls_out[-tmp_id_xml], 
       xml = tmp_ch_fls_out[tmp_id_xml], 
       stringsAsFactors = FALSE) 

    }) 

    ## join and return names of processed files 
    ch_fls_out <- do.call("rbind",ls_fls_out) 
    return(ch_fls_out) 
} 

getwd() 

setwd("C:/Users/joaoreis/Documents/Bases_Geograficas/trmm_3h/") 

fls_trmm <- downloadTRMM(begin = "2008-01-01.00", end = "2008-01-05.00") 
fls_trmm 

Aber ich erhalte den folgenden Fehler:

trying URL ' ftp://disc2.nascom.nasa.gov/data/TRMM/Gridded/3B42_V7//200801//3B42.20080102.00z.7.precipitation.bin ' Error in download.file(tmp_ch_fls[j], tmp_ch_fls_out[j], mode = "wb") : cannot open URL ' ftp://disc2.nascom.nasa.gov/data/TRMM/Gridded/3B42_V7//200801//3B42.20080102.00z.7.precipitation.bin ' In addition: Warning message: In download.file(tmp_ch_fls[j], tmp_ch_fls_out[j], mode = "wb") : InternetOpenUrl failed: '' Called from: download.file(tmp_ch_fls[j], tmp_ch_fls_out[j], mode = "wb")

Wer weiß, wie es zu beheben R verwenden?

Danke!

Answer 1

Ab dem Commit 909f98a habe ich den automatisierten Abruf von 3-stündlichen Daten von ftp://disc3.nascom.nasa.gov/data/s4pa/TRMM_L3 aktiviert. Stellen Sie sicher, dass Sie die neueste Version von Rsenal installiert mit

devtools::install_github("environmentalinformatics-marburg/Rsenal") 

und dann bei den Beispielen in ?downloadTRMM einen Blick. Für den Moment unterstützt die Funktion sowohl character (erfordert 'Format' Argument übergeben an strptime) und POSIXlt Eingabe. Zum Beispiel, so etwas wie

downloadTRMM(begin = "2015-01-01 12:00", end = "2015-01-03 12:00", 
      type = "3-hourly", format = "%Y-%m-%d %H:%M") 

3-stündliche Daten von 1 bis 3 Januar 2015 (Mittag bis Mittag) zum Download soll jetzt gut funktionieren.

Beachten Sie, dass die Daten im Gegensatz zu dem von Ihnen genannten FTP-Server im .HDF-Format vorliegen und eine -Methode bisher nicht implementiert wurde, so dass Sie die Container-Dateien selbst bearbeiten müssen. Ich werde versuchen, bald etwas bequemer in Bezug auf die automatisierte Rasterung der Daten herauszufinden.