Преобразование SpatialPointsDataframe в SpatialPolygons в R

У меня есть SpatialPointsDataFrame в R, который выглядит так:

              coordinates id order  hole piece group box_id
    326  (-94.4, 27.6586) 47     1 FALSE     1  47.1      1
    327 (-93.64, 27.6232) 47     2 FALSE     1  47.1      1
    328 (-92.04, 27.7649) 47     3 FALSE     1  47.1      1
    329 (-90.36, 27.0903) 47     4 FALSE     1  47.1      1
    330 (-91.12, 25.6929) 47     5 FALSE     1  47.1      1
    331 (-92.92, 25.6569) 47     6 FALSE     1  47.1      1
    332 (-93.44, 26.0169) 47     7 FALSE     1  47.1      1
    333  (-94.4, 25.9809) 47     8 FALSE     1  47.1      1
    334  (-94.4, 27.6586) 47     9 FALSE     1  47.1      1
    335 (-92.04, 27.7649) 48     1 FALSE     1  48.1      2
    336 (-93.64, 27.6232) 48     2 FALSE     1  48.1      2
    337  (-94.4, 27.6586) 48     3 FALSE     1  48.1      2
    338  (-94.4, 27.8356) 48     4 FALSE     1  48.1      2
    339 (-93.64, 27.7649) 48     5 FALSE     1  48.1      2
    340 (-90.28, 28.1182) 48     6 FALSE     1  48.1      2
    341 (-90.56, 27.9417) 48     7 FALSE     1  48.1      2
    342 (-92.04, 27.7649) 48     8 FALSE     1  48.1      2
    100  (-94.4, 27.8356) 20     1 FALSE     1  20.1      3
    101  (-94.4, 28.0829) 20     2 FALSE     1  20.1      3
    102 (-90.28, 28.1182) 20     3 FALSE     1  20.1      3
    103 (-93.64, 27.7649) 20     4 FALSE     1  20.1      3
    104  (-94.4, 27.8356) 20     5 FALSE     1  20.1      3

(имена/номера строк не в порядке, потому что я отсортировал по столбцу box_id)

Эти точки являются узлами полигонов (обозначаются box_id). Я хочу написать это как файл формы многоугольника для чтения в программах ГИС, но у меня возникли проблемы с преобразованием его в SpatialPolygonDataFrame (чтобы затем использовать writeOGR) или напрямую записать его в файл shp. Любая помощь будет оценена. Я новичок в ГИС с R, поэтому прошу прощения, если упустил что-то очевидное.


r gis spatial polygons
person user3306720    schedule 13.02.2014    source источник

Ответы (2)


arrow_upward
5
arrow_downward

Итак, вот еще одно решение, концептуально похожее на @JoshO'Brien. Этот вмещает подполигоны (несколько групп для заданного идентификатора) и создает объект класса SpatialPolygonsDataFrame с включенным разделом данных (например, таблица атрибутов). Это полностью экспортируется в шейп-файл ESRI.

Во-первых, в качестве примера мы создаем SpatialPointsDataFrame с таким же форматом, как у вас. Столбец id определяет полигон ID, а столбец group определяет подполигоны для данного полигона. В этом примере у нас есть 3 полигона, два из которых имеют 1 субполигон, а третий имеет 3 субполигона. Мы также создаем фрейм данных data (таблица атрибутов), который в этом случае связывает ID полигона с box_id.

# this code just creates a sample SpatialPointsDataFrame 
x <- c(-10,-10,10,10,-10)
y <- c(-10,10,10,-10,-10)
df.1 <- data.frame(x,y,id=47, order=1:5,hole=F,piece=1,group=47.1,box_id=1)
coordinates(df.1)=c("x","y")
x <- c(+15+3*cos(2*pi/5*(0:5)))
y <- c(-15+3*sin(2*pi/5*(0:5)))
df.2 <- data.frame(x,y,id=48, order=1:6,hole=F,piece=1,group=48.1,box_id=2)
coordinates(df.2)=c("x","y")
x <- c(-15+2*cos(2*pi/8*(0:8)))
y <- c(+15+2*sin(2*pi/8*(0:8)))
df.3.1 <- data.frame(x,y,id=20, order=1:9,hole=F,piece=1,group=20.1,box_id=3)
coordinates(df.3.1)=c("x","y")
x <- c(0+2*cos(2*pi/8*(0:8)))
y <- c(+15+2*sin(2*pi/8*(0:8)))
df.3.2 <- data.frame(x,y,id=20, order=1:9,hole=F,piece=1,group=20.2,box_id=3)
coordinates(df.3.2)=c("x","y")
x <- c(+15+2*cos(2*pi/8*(0:8)))
y <- c(+15+2*sin(2*pi/8*(0:8)))
df.3.3 <- data.frame(x,y,id=20, order=1:9,hole=F,piece=1,group=20.3,box_id=3)
coordinates(df.3.3)=c("x","y")
df <- rbind(df.1,df.2,df.3.1,df.3.2,df.3.3)
df
#              coordinates id order  hole piece group box_id
# 1             (-10, -10) 47     1 FALSE     1  47.1      1
# 2              (-10, 10) 47     2 FALSE     1  47.1      1
# 3               (10, 10) 47     3 FALSE     1  47.1      1
# 4              (10, -10) 47     4 FALSE     1  47.1      1
# 5             (-10, -10) 47     5 FALSE     1  47.1      1
# 6              (18, -15) 48     1 FALSE     1  48.1      2
# 7  (15.92705, -12.14683) 48     2 FALSE     1  48.1      2
# 8  (12.57295, -13.23664) 48     3 FALSE     1  48.1      2
# ...
data <- data.frame(box_id=unique(df$box_id),row.names=unique(df$id))

Теперь это код для points2polygons(...). С вашим существующим SpatialPointsDataFrame это все, что вам нужно.

points2polygons <- function(df,data) {
  get.grpPoly <- function(group,ID,df) {
         Polygon(coordinates(df[df$id==ID & df$group==group,]))
         }
  get.spPoly  <- function(ID,df) {
         Polygons(lapply(unique(df[df$id==ID,]$group),get.grpPoly,ID,df),ID)
         }
  spPolygons  <- SpatialPolygons(lapply(unique(df$id),get.spPoly,df))
  SpatialPolygonsDataFrame(spPolygons,match.ID=T,data=data)
}
spDF <- points2polygons(df,data)
plot(spDF,col=spDF$box_id+1)

library(rgdal)
writeOGR(spDF,dsn=".",layer="myShapefile", driver="ESRI Shapefile")

person jlhoward    schedule 13.02.2014
comment
Спасибо за отзыв. - person user3306720; 14.02.2014
comment
Я выбрал этот, чтобы включить полный процесс написания с использованием writeOGR. Я думаю, что ( в конце строки 2 должен быть { но в остальном он работает отлично - person user3306720; 14.02.2014

arrow_upward
0
arrow_downward

Проверьте это:

Из: SpatialPointsDataFrame --› SpatialPolygons --› SpatialPolygonsDataFrame

https://stackoverflow.com/a/66565985/14361772

person Yuri Gelsleichter    schedule 10.03.2021