我在R中有一个SpatialPointsDataFrame对象,它的结构如下:
coordinates id order hole piece group box_id
326 (-94.4, 27.6586) 47 1 FALSE 1 47.1 1
327 (-93.64, 27.6232) 47 2 FALSE 1 47.1 1
328 (-92.04, 27.7649) 47 3 FALSE 1 47.1 1
329 (-90.36, 27.0903) 47 4 FALSE 1 47.1 1
330 (-91.12, 25.6929) 47 5 FALSE 1 47.1 1
331 (-92.92, 25.6569) 47 6 FALSE 1 47.1 1
332 (-93.44, 26.0169) 47 7 FALSE 1 47.1 1
333 (-94.4, 25.9809) 47 8 FALSE 1 47.1 1
334 (-94.4, 27.6586) 47 9 FALSE 1 47.1 1
335 (-92.04, 27.7649) 48 1 FALSE 1 48.1 2
336 (-93.64, 27.6232) 48 2 FALSE 1 48.1 2
337 (-94.4, 27.6586) 48 3 FALSE 1 48.1 2
338 (-94.4, 27.8356) 48 4 FALSE 1 48.1 2
339 (-93.64, 27.7649) 48 5 FALSE 1 48.1 2
340 (-90.28, 28.1182) 48 6 FALSE 1 48.1 2
341 (-90.56, 27.9417) 48 7 FALSE 1 48.1 2
342 (-92.04, 27.7649) 48 8 FALSE 1 48.1 2
100 (-94.4, 27.8356) 20 1 FALSE 1 20.1 3
101 (-94.4, 28.0829) 20 2 FALSE 1 20.1 3
102 (-90.28, 28.1182) 20 3 FALSE 1 20.1 3
103 (-93.64, 27.7649) 20 4 FALSE 1 20.1 3
104 (-94.4, 27.8356) 20 5 FALSE 1 20.1 3
(行名 / 行号不按顺序,因为我按 box_id 列排序)
这些点是多边形的节点(由 box_id 确定)。 我想将其写成多边形形状文件以便读入 GIS 程序,但我无法将其转换为 SpatialPolygonDataFrame(然后使用 writeOGR)或直接将其写入 shp 文件。如有帮助会非常感激。我是 R 中的 GIS 新手,如果我漏掉了一些显而易见的东西,请谅解。