Добавление полигона в SpatialPolygonsDataFrame

Question

Как добавить многоугольник в SpatialPolygonsDataFrame?

Пример: скрипт внизу создаст SpatialPolygonsDataFrame. Я хотел бы добавить многоугольник, который состоит из большого квадрата вокруг существующих многоугольников.

   library(rgdal)
   dsn <- system.file("vectors", package = "rgdal")[1]
   Scotland <- readOGR(dsn=dsn , layer="scot_BNG")
   plot(Scotland)

Предпочтительный конечный результат:

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

Спасибо!

Answer 1

Следующий код создает прямоугольник, который охватывает исходные пространственные многоугольники, и добавляет его как пространственный многоугольник к исходной форме.

library(rgdal)
library(rgeos)

dsn <- system.file("vectors", package = "rgdal")[1]
Scotland <- readOGR(dsn=dsn , layer="scot_BNG")

# change the width parameter to make the rectangle the desired size
# this results in an extent object that surrounds the original shape
eScotland <- extent(gBuffer(Scotland, width = 50000))
# turn the extent into a Spatial Polygon
pScotland <- as(eScotland, 'SpatialPolygons')
crs(pScotland) <- crs(Scotland)
newScotland <- bind(pScotland, Scotland)
plot(newScotland)