Gewusst wie: extrahieren von Werten aus rasterstack mit xy-Koordinaten?

Ich habe eine rasterstack (5 raster-Layer), die tatsächlich eine Zeit-Serie raster. 

r <- raster(nrow=20, ncol=200)
s <- stack( sapply(1:5, function(i) setValues(r, rnorm(ncell(r), i, 3) )) )
s

class       : RasterStack
dimensions  : 20, 200, 4000, 5  (nrow, ncol, ncell, nlayers)
resolution  : 1.8, 9  (x, y)
extent      : -180, 180, -90, 90  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
names       :   layer.1,   layer.2,   layer.3,   layer.4,   layer.5
min values  : -9.012146, -9.165947, -9.707269, -7.829763, -5.332007
max values  :  11.32811,  11.97328,  15.99459,  15.66769,  16.72236

Mein Ziel ist es, zum zeichnen jedes pixel und erforschen Ihr Verhalten im Laufe der Zeit.

Wie konnte ich das extrahieren der einzelnen Pixel, die zusammen mit Ihren x -, y-Koordinaten und plot-Zeitreihen-Kurve?

InformationsquelleAutor Eddie | 2013-11-07
  1. 9

    Können Sie extract bestehen und ein Vektor, der die Anzahl der Zellen, die Sie extrahieren wollen, um wieder eine matrix der Werte in jedem pixel. Jede Zeile repräsentiert ein pixel, die Spalten sind Schichten...

    mat <- extract( s , 1:ncell(s) )
head( mat )
#        layer.1  layer.2  layer.3    layer.4   layer.5
#[1,] -0.2138718 3.114061 3.670945  1.2560295  2.881104
#[2,]  3.3580783 5.008205 2.315353  2.3247236 11.539837
#[3,]  3.2173875 2.958985 1.055389  3.1016730  4.064339
#[4,]  4.1113162 4.469828 3.113790  8.5329679  8.771459
#[5,] -2.4011283 4.747527 4.299707  2.2111643  9.457012
#[6,] -2.6159294 5.659211 1.926900 -0.3886837  5.661419

    Jedoch extract ist mehr nützlich, wenn Sie versuchen, um bestimmte Pixel. Um alle Pixel mit der x /y Koordinaten können Sie einfach rasterToPoints...

    head( rasterToPoints( s ) )
#          x    y    layer.1  layer.2  layer.3    layer.4   layer.5
#[1,] -179.1 85.5 -0.2138718 3.114061 3.670945  1.2560295  2.881104
#[2,] -177.3 85.5  3.3580783 5.008205 2.315353  2.3247236 11.539837
#[3,] -175.5 85.5  3.2173875 2.958985 1.055389  3.1016730  4.064339
#[4,] -173.7 85.5  4.1113162 4.469828 3.113790  8.5329679  8.771459
#[5,] -171.9 85.5 -2.4011283 4.747527 4.299707  2.2111643  9.457012
#[6,] -170.1 85.5 -2.6159294 5.659211 1.926900 -0.3886837  5.661419
    • Vielen Dank @SimonO101.
    • Ohhh @SimonO101...gerade festgestellt, dass rasterToPoints nicht extrahieren NA Werte, falls Sie existieren in jeder raster-layer. Dies wäre ein problem für mich, da muss ich produzieren zurück, das raster mit der gleichen dimension. Beispiel: r <- raster(nrow=10, ncol=10) s < - - stack( sapply(1:2, function(i) setValues(r, rnorm(ncell(r), i, 3) )) ) s plot(s) s[1:3]<-NA-plot(s) pix.Werte<-rasterToPoints(s) pix.Werte #nur 97 Zeilen generiert wurden. 3 NA-Werte wurden ausgelassen. Jemand bereit zu werfen etwas Licht? Vielen Dank im Voraus
    • verwenden Sie die extract Methode mit na.rm = FALSE und cbind dass xyFromCell( 1:ncell(s) )
    InformationsquelleAutor Simon O'Hanlon
  2. 6

    Vielen Dank @SimonO101 !

    Dieser code funktioniert.

    r <- raster(nrow=10, ncol=10)
s <- stack( sapply(1:5, function(i) setValues(r, rnorm(ncell(r), i, 3) )) )
s[1:3]<-NA

vals <- values(s)

coord <- xyFromCell(s,1:ncell(s))

combine <- cbind(coord,vals)

write.table(combine,"xyvalues.txt")
    InformationsquelleAutor Eddie
Schreibe einen Kommentar