2016-04-09 3 views
0

Ich versuche, das Beispiel geben here unter „einem polaren Beispiel“ folgenNSIDC Meereis .bin-Datei Raster()

Dieses Beispiel verwendet Meereis Daten in .bin-Format als Raster zu zeichnen. Ich versuche das gleiche mit einer anderen Datei, die von der original ftp server des Nationalen Schnee und Eis Datenzentrums verfügbar ist. Daher nehme ich an, dass es kein Problem geben sollte. Aber wenn ich versuche, das folgende Skript

# from NSIDC sea ice concentration file 
baseurl <- "ftp://sidads.colorado.edu/pub/DATASETS/" 

f2 <- paste(baseurl, 
     "nsidc0051_gsfc_nasateam_seaice/final-gsfc/north/daily/2013/nt_20130111_f17_v1.1_n.bin", 
     sep='') 

if (!file.exists(basename(f2))) download.file(f2, basename(f2), mode = "wb") 
ice2 <- raster(basename(f2)) 

Fehler in .rasterObjectFromFile (x, Band = Band, Objecttype = "Rasterebene" aufzufordern R laufen,.: Es kann keine Rasterebene Objekt aus dieser Datei erstellen

Wo bin ich falsch gehen? ist die .bin-Datei beschädigt? Jede Hilfe dankbar!

Dank!

Antwort

0

ok, eine Lösung auf github gefunden, die funktioniert wirklich nett.

https://github.com/cran/raster/blob/master/R/nsidcICE.R

Ersetzen Sie einfach in Zeile 14 des Skripts

hemi <- tolower(substr(bx, 21L, 21L)) 

von

hemi <- tolower(substr(bx, 22L, 22L)), 

als die neue Struktur Name um eine Stelle namly ein Punkt etwas von der ursprünglichen unterscheidet sich in die v1.1 Sequenz!

vergleichen:

"nt_19781119_f07_v01_s.bin" 

auf die Version, die ich in

interessiert war
"nt_20130111_f17_v1.1_n.bin"