Eu tenho um shapefile (mostrando diferentes classes de sedimentos no mar do Norte) lendo comreadOGR(). Tem muitos buracos "o que deve ser" em muitos polígonos, mas usandorasterize() elimina todos os furos já que eles não estão marcados comoTRUE nas fendas dos buracos. Usavarasterize(...,fun='first') sem sucesso. No entanto, o QGIS mostra todos os buracos bem. Além disso,over() corretamente avalia os valores do campo, por exemplo, em um buraco, provavelmente aproveitando a "ordem de plotagem" do slot, e é por isso que eu criei algo como:

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
}

para ignorar o rasterize, mas vai demorar 50 dias para ser feito.

Então aqui 'minha pergunta:

Alguém já experimentou algo semelhante e encontrou uma solução para isso?

(Eu gostaria de incluir dados de exemplo, masdput()falha em SpatialPolygons?!?)