r - gSimplify 不简化 R 中的 shapefile

Question

我无法在 R 中简化 shapefile

来自此处的形状文件:https://geoportal.statistics.gov.uk/Docs/Boundaries/Local_authority_district_(GB)_2014_Boundaries_(Generalised_Clipped).zip

library(tmap)
library(maptools)
library(ggmap)

England <- readOGR(dsn = "...")

#works fine
print(qtm(England, "LAD14CD", borders = NA, fill.title = "A-Level" )) 

# simplify the polygons
England<-gSimplify(England,tol=0.01, topologyPreserve=TRUE)

print(qtm(England, "LAD14CD", borders = NA, fill.title = "A-Level" )) 

给出一个错误:

Error in process_fill(data, g$tm_fill, gborders, gt, gf, z = z + which(plot.order ==  : 
Fill argument neither colors nor valid variable name(s)

如果您查看 UK 数据对象，您会发现它已从大型空间多边形数据帧更改为大型空间多边形并删除了 @data

相反，如果您尝试仅简化 Shapefile 中的多边形:

England@polygons<-gSimplify(England@polygons,tol=0.01, topologyPreserve=TRUE)

上面写着:

Error in gSimplify(England@polygons, tol = 0.01, topologyPreserve = TRUE) : 
cannot get a slot ("proj4string") from an object of type "list"

如何简化 shapefile 中的多边形？

Answer 1

gSimplify 的返回只是几何图形，而不是属性，因此您必须使用简化的几何图形和原始属性数据构造一个新的 SpatialPolygonsDataFrame:

> England2 <-gSimplify(England,tol=0.01, topologyPreserve=TRUE)
> England3 = SpatialPolygonsDataFrame(England2, data=England@data)

我认为多边形可以保证是相同的顺序，除非有任何东西被简化了。检查 length(England2) 的行数是否与 England 相同，或匹配 ID 上的行。