Как упростить пространственный многоугольник, используя rmapshaper? - PullRequest
Новая
       40

Как упростить пространственный многоугольник, используя rmapshaper?

0 голосов
/ 20 апреля 2020

Я пытаюсь упростить структуру данных пространственных полигонов с помощью следующего кода: 

library(rmapshaper)
library(rgdal)
countyshapes <- readOGR("https://raw.githubusercontent.com/kjhealy/us-county/master/data/geojson/gz_2010_us_050_00_500k.json")
countyshapes <- ms_simplify(countyshapes, keep = 0.1)

Но это возвращает следующую ошибку: 

Error in sp::SpatialPolygonsDataFrame(polys, data = input@data) : 
  row.names of data and Polygons IDs do not match

Есть ли способ, которым я могу обойти эта ошибка для упрощения моего объекта SP?

1 Ответ

0 голосов
/ 22 апреля 2020

Я запустил код и, если он сработал, желательно обновить пакет rmapshaper: 

packageVersion("rmapshaper") 
#[1] ‘0.4.4’

countyshapes[1:5,]
#An object of class "SpatialPolygonsDataFrame"
#Slot "data":
#  GEO_ID STATE COUNTY      NAME   LSAD CENSUSAREA
#1 0500000US01029    01    029  Cleburne County    560.100
#2 0500000US01031    01    031    Coffee County    678.972
#3 0500000US01037    01    037     Coosa County    650.926
#4 0500000US01039    01    039 Covington County   1030.456
#5 0500000US01041    01    041  Crenshaw County    608.840

Другая альтернатива - использование функции sf::st_transform

...