Получение сообщения об ошибке при использовании команды over ()

Question

У меня есть два файла формы, shp1 и shp2. После вычисления центроида первого я хотел бы увидеть, где центроид находится во втором.

Моя настройка:

shp1@data$centroid <- gCentroid(shp1, byid = TRUE)
foo <- over(shp1$centroid, shp2)

После выполнения второй строки я получаюследующее сообщение об ошибке:

Error in xj[i] : NAs not permitted in row index

Любая идея, что я должен искать, чтобы избежать этой ошибки? (примечание: grep("NA", row.names(shp1)) возвращает integer(0), также для другого файла формы.)

Answer 1

Просто выяснил, почему возникает эта ошибка:

В дополнение к shp1@data$centroid <- gCentroid(shp1, byid = TRUE) я также сделал shp2@data$centroid <- gCentroid(shp2, byid = TRUE) перед вызовом foo <- over(shp1$centroid, shp2), что, очевидно, вызвало ошибку.