Как прикрепить простой data.frame к SpatialPolygonDataFrame в R?

Пусть df = фрейм данных, sp = объект пространственного полигона и by = имя или номер столбца общего столбца. Затем вы можете объединить фрейм данных с объектом sp, используя следующую строку кода

sp@data = data.frame(sp@data, df[match(sp@data[,by], df[,by]),])

Вот как работает код. Функция match внутри выравнивает столбцы так, чтобы порядок сохранялся. Поэтому, когда мы объединяем его с sp@data, порядок сохраняется правильно. Быстрая проверка того, сработал ли код, заключается в проверке двух столбцов, соответствующих общему столбцу, и проверке их идентичности (общие столбцы дублируются, и их легко удалить, но я сохраняю их такими, какие они есть). хороший чек)

Это так же просто:

require(sp) # the trick is that this package must be loaded!

oo <- merge(SPDF,DF, by="QDGC")

Я испытал на себе. Но это работает, только если вы используете слияние из пакета sp. Это значение по умолчанию при загрузке пакета sp. Затем функция merge перегружается, а sp::merge используется, если первым аргументом является пространственная структура.

слияние может создать фрейм данных с большим количеством строк, чем оригиналы, если нет простого сопоставления 1-1 двух фреймов данных. В этом случае пришлось бы копировать всю геометрию и создавать несколько полигонов, что, вероятно, не очень хорошо.

Если у вас есть фрейм данных с тем же количеством строк, что и SpatialPointsDataFrame, вы можете просто напрямую заменить слот @data.

library(sp)
example(overlay) # to get the srdf object
srdf@data
spplot(srdf)
srdf@data=data.frame(x=runif(3),xx=rep(0,3))
spplot(srdf)

если вы ошиблись с количеством строк:

srdf@data=data.frame(x=runif(2),xx=rep(0,2))
spplot(srdf)
Error in data.frame(..., check.names = FALSE) : 
  arguments imply differing number of rows: 3, 2

Возможно, функция joinCountryData2Map в пакете rworldmap может дать вдохновение. (Но я могу ошибаться, как и в прошлый раз.)

Еще одно решение — использовать функцию append_data из пакета tmaptools . Он вызывается с такими аргументами:

append_data(shp, data, key.shp = NULL, key.data = NULL,
  ignore.duplicates = FALSE, ignore.na = FALSE,
  fixed.order = is.null(key.data) && is.null(key.shp))

Немного прискорбно, что это называется append, так как я понимаю добавление больше в смысле rbind, и мы хотим, чтобы здесь было что-то вроде join или merge.

Игнорируя этот факт, функция действительно полезна для проверки правильности ваших соединений и наличия некоторых строк только на одной стороне соединения. Из документов:

Недостаточное покрытие (элементы формы, которые не соответствуют записям данных), избыточное покрытие (записи данных, которые не соответствуют элементам формы соответственно), а также наличие повторяющихся значений ключей автоматически проверяются и сообщаются с помощью сообщений консоли. С помощью under_coverage и over_coverage можно получить значения ключа недостаточного и избыточного покрытия из последнего вызова append_data,

Если необходимо объединить два шейп-файла в один объект, просто используйте rbind().

При использовании rbind() просто убедитесь, что оба аргумента, которые вы используете, равны SpatialDataFrames. Вы можете проверить это с помощью class(sf). Если это не кадр данных, используйте st_as_sf(), чтобы преобразовать их в SpatialDataFrame, прежде чем привязывать их.

Примечание. Вы также можете использовать это для добавления к NULLs, особенно когда вы используете результат из цикла и хотите суммировать результаты.