0

所以我正在绘制我国的邮政区域，这是一个巨大的 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)
4

0 回答 0