因此,我正在绘制我国的邮政地区,它是一个巨大的,540万行的数据帧。我不需要每一个小小的扭曲和转弯的边界,所以我非常兴奋地发现mapshaper.org。看似很容易,我只是把我的.shp文件扔进去了,它吐出了10%大小的文件!优秀!
没有..现在它没有玩readShapeSpatial ...希望我不需要提供shapefile(我可以跟踪它,如果这些部分真的有必要的话)
原始代码,作品:
canada.shp <- readShapeSpatial("canada", verbose=TRUE, delete_null_obj=TRUE)
新地图,不是:
canada.shp <- readShapeSpatial("simplecanada", verbose=TRUE, delete_null_obj=TRUE)
Shapefile type: Polygon, (5), # of Shapes: 1621
Error in `rownames<-`(`*tmp*`, value = c("0", "1", "2", "3", "4", "5", :
attempt to set 'rownames' on an object with no dimensions
In addition: Warning message:
In .Map2PolyDF(Map, IDs = IDvar, proj4string = proj4string, force_ring = force_ring, :
Null objects with the following indices deleted: 895
traceback()显示:
4: stop("attempt to set 'rownames' on an object with no dimensions")
3: `rownames<-`(`*tmp*`, value = c("0", "1", "2", "3", "4", "5",
"6", "7", "8", "9", "10", "11", "12", "13", "14", "15", "16",
"17", "18", "19", "20", "21", "22", "23", "24", "25", "26", "27",
"28", "29", "30", "31", "32", "33", "34", "35", "36", "37", "38",
"39", "40", "41", "42", "43", "44", "45", "46", "47", "48", "49",
"50", "51", "52", "53", "54", "55", "56", "57", "58", "59", "60", (ETC!!!!)
2: .Map2PolyDF(Map, IDs = IDvar, proj4string = proj4string, force_ring = force_ring,
delete_null_obj = delete_null_obj, retrieve_ABS_null = retrieve_ABS_null)
1: readShapePoly("simplecanada",
verbose = TRUE, delete_null_obj = TRUE)