如何使用 rmapshaper 简化空间多边形？
Posted
tags:
篇首语：本文由小常识网(cha138.com)小编为大家整理，主要介绍了如何使用 rmapshaper 简化空间多边形？相关的知识，希望对你有一定的参考价值。
我试图用下面的代码简化一个空间多边形数据框。
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对象？
答案
我运行了这段代码，如果能成功的话，最好更新一下 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 功能
以上是关于如何使用 rmapshaper 简化空间多边形？的主要内容，如果未能解决你的问题，请参考以下文章