У меня есть шейп-файл (показывающий различные классы отложений в северном море), прочитанный сreadOGR(), Оно имеет много & quot; что должно быть & quot; дыры во многих полигонах, но с использованием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 дней.

Так вот & apos; мой вопрос:

Кто-нибудь испытывал нечто подобное и нашел обходной путь для этого?

(Я хотел бы включить пример данных, ноdput()терпит неудачу на SpatialPolygons?!?)