Utwórz wielokąty zasięgu i granicy danego rastra w R.

Chcę utworzyć dwa wielokąty.

Jeden jest z prostokątnych rozmiarów mojego rastra. Wiem, że mogę uzyskać zasięg użycia rastra, r <- raster("band5.tif") e <- extent(r) plot(e)daje mi to:
Jak mogę utworzyć JEDEN wielokąt granicy rastra, jak pokazano poniżej?

Tak wygląda mój raster:

raster r rgdal

— csheth
źródło

Oto przykład.

library(raster)
# example data
x <- raster(system.file("external/test.grd", package="raster"))

Aby uzyskać zasięg prostokątny

e <- extent(x)
# coerce to a SpatialPolygons object
p <- as(e, 'SpatialPolygons')

Aby uzyskać wielokąt otaczający komórki, które nie są NA

# make all values the same. Either do
r <- x > -Inf
# or alternatively
# r <- reclassify(x, cbind(-Inf, Inf, 1))

# convert to polygons (you need to have package 'rgeos' installed for this to work)
pp <- rasterToPolygons(r, dissolve=TRUE)

# look at the results
plot(x)
plot(p, lwd=5, border='red', add=TRUE)
plot(pp, lwd=3, border='blue', add=TRUE)

— Robert Hijmans
źródło

Nie mogę znaleźć dokumentacji z tej r <- r > -Infstrony. Co to właściwie robi? I jak różni się od values(r)[!is.na(values(r))] <- 1(co ustawia wszystkie lokalizacje, które nie są NA w r na 1).

— csheth

r > -Infjest podstawowy R. Sprawdź, c(1,NA,3,NA) > -Infjak to działa. Dodałem alternatywę ( reclassify). Twoja alternatywa działa, ale nie jest dobra dla dużych obiektów.

— Robert Hijmans,

Czekam tylko, aż komputer przetworzy tę rasterToPolygons()część, i wtedy ocenię twoją odpowiedź. Dzięki za -Infkawałek będzie bardzo przydatny!

— csheth

W przypadku bardzo dużego rastra może to chwilę potrwać lub nawet zawieść. W takim przypadku możesz rozważyć użycie najpierw, aggregateponieważ i tak nie zobaczysz różnicy.

— Robert Hijmans,

Nawet z zaledwie 100 komórkami po uruchomieniu aggregate()on r, nie przetworzył (komputer wyposażony w 16 GB pamięci RAM). Problem polega na tym, że muszę zachować granicę tak natywną, jak to tylko możliwe, ponieważ DEM musi zostać wyodrębniony z powstałego wielokąta. Jeśli spróbuję dalej pobrać próbkę rastra, utracę tę granicę. Masz jakieś prace?

— csheth,