Mam wczytany shapefile (pokazujący różne klasy osadów na Morzu Północnym)readOGR(). Ma wiele „co powinno być” dziur w wielu wielokątach, ale używarasterize() eliminuje wszystkie dziury, ponieważ nie są oznaczone jakoTRUE w ich otworach. Używanyrasterize(...,fun='first') bez powodzenia. Niemniej jednak QGIS ładnie pokazuje dziury. Również,over() poprawnie ocenia wartości pola, np. w otworze, prawdopodobnie wykorzystując slot „kolejność wydruku”, dlatego wymyśliłem coś takiego:

for (i in 1:ncell(raster)){
    coo<-xyFromCell(raster,i,spatial=T)
    col<-colFromX(ra,coo@coords[1,1])
    row<-rowFromY(ra,coo@coords[1,2])
    proj4string(coo)<-proj4string(shape)
    n<-over(coo,shape)
    raster[col,row]<-n$Prime_FOLK
}

aby pominąć rasteryzację, ale zajmie to 50 dni.

Więc tutaj moje pytanie:

Czy ktoś doświadczył czegoś podobnego i znalazł dla niego obejście?

(Chciałbym dodać dane przykładowe, aledput()nie działa na SpatialPolygons?!?)