将R图中的多边形导出为shapefile

Question

我一直试图将地块(行/多边形)的内容导出为可以在ArcMap中打开的图层/shapefile。这些是我一直在使用的一些图书馆，

库(Shapefiles)库(PBSmapping)库(AdehabitatHR)库(Maptools)库(映射)库(Rgdal)库

我的lat/长数据如下：

tagdata<-read.table(text=“均数”均值为长-18.63327 147.0248 -18.63327 147.0238 -18.63327 147.294 -text= 147.2942 - 147.2942 -旋即147.2938 -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.63063 147.1164 -18.63063 147.0238 -.0256 147.0256 - 147.0256 -18.7538147.2541-18.753147.2541-18.7538147.2541-18.7538147.2541 -18.63063 147.0256 -18.63063 147.0256 -18.63063 147.0256 -18.63063 147.0256-18.63063 147.0256 -18.75383 147.2541 -18.75383 147.2541 -18.63063 147.0256 -18.68133 147。1164 -18.68133 147.1164 -18.68133 147.0251 -18.68133 147.0909 -1164 147.2816 -488- 147.0792 -专属147.0896 - 147.0238 -18.75383 147.2541 -18.75383 147.2541 -18.75383 147.2541"，header=TRUE)

我绘制了位置图，并使用AdehabitatHR包计算了最小凸多边形(MCP)。

Col=<-mcp(coord，percent=100)图( add=TRUE)

我知道如何以shapefile的形式导出/写点，我可以在ArcMap或类似的软件中打开，

例如：

loc <-空间点(Loc)#将loc转换为空间点rem<-tagdatac(-1:-2) SpatialPointsDataFrame(coords=loc，data=rem) obj<-SpatialPointsDataFrame(coords=loc，data=rem) writePointsShape(obj，"myshape.shp")

但是，我还没有找到一种很好的方法来处理多边形o多边形对象。我希望能够导出/编写带有MCP作为shapefile的poly。有什么建议吗？

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")

您可以很容易地使用它来编写任何形状文件(点、多边形等)，但它们必须是SpatialxxxDataFrame对象。

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