将R绘图中的多边形导出为shape文件

Question

我一直试图将绘图（线条/多边形）的内容导出为可在ArcMap中打开的图层/ shapefile。这些都是一些，我一直在使用的库，

library(shapefiles) 
library(PBSmapping) 
library(adehabitatHR) 
library(maptools) 
library(maps) 
library(rgdal) 
library(igraph) 

我的纬度/经度数据是这样的：

tagdata<-read.table(text="meanlat meanlong 
-18.63327 147.0248 
-18.6368 147.0238 
-18.62068 147.294 
-18.62953 147.2942 
-18.62953 147.2942 
-18.62091 147.2938 
-18.62953 147.2942 
-18.62466 147.2926 
-18.73393 147.2816 
-18.73393 147.2816 
-18.75383 147.2541 
-18.75383 147.2541 
-18.75383 147.2541 
-18.75383 147.2541 
-18.6368 147.0238 
-18.63063 147.0256 
-18.63063 147.0256 
-18.68133 147.1164 
-18.6368 147.0238 
-18.63063 147.0256 
-18.63063 147.0256 
-18.75383 147.2541 
-18.61273 147.0682 
-18.69655 147.09 
-18.6368 147.0238 
-18.63063 147.0256 
-18.63063 147.0256 
-18.63217 147.0251 
-18.75383 147.2541 
-18.75383 147.2541 
-18.75383 147.2541 
-18.63063 147.0256 
-18.68133 147.1164 
-18.68133 147.1164 
-18.63217 147.0251 
-18.69922 147.0909 
-18.73393 147.2816 
-18.63632 147.0792 
-18.69522 147.0896 
-18.6368 147.0238 
-18.75383 147.2541 
-18.75383 147.2541 
-18.75383 147.2541",header=TRUE) 

我绘制的位置和计算的使用AdehabitatHR软件包的最小凸多边形（MCP）。

plot(tagdata$meanlong,tagdata$meanlat, col="red",pch=1)     

loc<-tagdata[ ,c("meanlong","meanlat")] 
coord<-SpatialPoints(loc) 
poly<-mcp(coord,percent=100) 
plot(poly,add=TRUE) 

我知道如何导出/写点作为shape文件，我可以在ArcMap或类似的软件打开，

例如：

loc<-SpatialPoints(loc) # #convert loc to spatial points 
rem<-tagdata[c(-1:-2)] 
SpatialPointsDataFrame(coords=loc,data=rem) 
obj<-SpatialPointsDataFrame(coords=loc,data=rem) 
writePointsShape(obj,"myshape.shp") 

然而，我还没有找到使用多边形或多边形对象的好方法。我希望能够将MCP的poly对象作为shapefile导出/写入。有什么建议么？

Answer 1

rgdal对于这类事情来说非常棒。 http://www.gdal.org/拥有大部分信息，您可以根据需要支持哪些格式。

在这种情况下你想要的writeOGR功能

# this will create a shapefile called poly within the working directory 
library(rgdal) 
writeOGR(poly, dsn = '.', layer = 'poly', driver = "ESRI Shapefile") 

你可以很方便地用它来写任何形状文件（点，多边形等），但它们必须​​对象

coorddf <- SpatialPointsDataFrame(coord, data = data.frame(dummy = rep(1,nrow(coord@coords)))) 
writeOGR(coorddf, dsn = '.', layer = 'mypoints', driver = "ESRI Shapefile")