r - ESRIシェープファイルを穴でラスタライズしますが、穴スロットはFALSEです

Question

シェープファイル (北海のさまざまな堆積物クラスを示す) を で読み込んでいreadOGR()ます。多くのポリゴンに「あるべき」穴がたくさんありますが、穴スロットrasterize()としてマークされていないため、使用するとすべての穴がなくなります。成功せずにTRUE使用されました。rasterize(...,fun='first')それにもかかわらず、QGIS はすべての穴をうまく表示します。また、over()穴などのフィールド値を正しく評価します。おそらく、スロットの「プロット順序」を利用しているため、次のようなものを思いつきました。

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
}

ラスタライズをバイパスしますが、完了するまでに 50 日かかります。

だからここに私の質問：

誰かが同様のことを経験し、その回避策を見つけましたか?

(例のデータを含めたかったdput()のですが、SpatialPolygons で失敗しますか?!?)