У меня есть Large SpatialPolygons объект rtu со списком из 161 полигона. В первом списке 146 Polygons. Я перемещаю один многоугольник в этот список (пытаясь добавить) из того же объекта: список индекса 92. В конце я удаляю эту строку.
r = rtu@polygons[[92]]@Polygons
rp <- Polygons(r, 1)
ps = c(rtu@polygons[[1]]@Polygons, rp@Polygons)
rtu@polygons[[1]]@Polygons <- ps
rtu@polygons[[92]] <- NULL
На данный момент все выглядит хорошо. Структура теперь имеет 147 Polygons и сюжет в порядке. Но пока я иду дальше и пробую:
over(rtu[1], rtu[2])
это выдает ошибку:
> ov = over(rtu[1], rtu[2])
Error in RGEOSBinPredFunc(spgeom1, spgeom2, byid, func) :
lengths of comment and Polygons slot differ
Я полагаю, что моя трансформация была неправильной в какой-то момент, но где? Я думаю, что есть некоторые функции, чтобы быть лучше с такими преобразованиями. Пожалуйста, помогите мне найти эту ошибку.