尝试将空间多边形导出为 R 中的 shapefile 时 writeOGR 出错

问题描述

我正在尝试导出我为家庭范围数据创建的多边形,但是当我导出它们时收到此错误消息:

Error in writeOGR(colkud95[1,],".","colonyhomerange",driver = "ESRI Shapefile") : 
  Spatial polygons objects should have valid ordering comments
use rgeos::createSPComment()
dput(HRtracks[1:10,])
new("SpatialPointsDataFrame",data = structure(list(TripID = c(1,1,1),`Breeding stage` = c(1,1)),row.names = c(NA,-10L),class = c("tbl_df","tbl","data.frame")),coords.nrs = 1:2,coords = structure(c(174.78327,174.78326,174.78332,174.78323,174.78325,174.78277,174.78238,174.78244,174.78278,-41.09482,-41.09491,-41.09484,-41.09496,-41.09564,-41.09637,-41.09685,-41.09747,-41.09771),.Dim = c(10L,2L),.Dimnames = list(NULL,c("x","y"))),bBox = structure(c(174.78238,-41.09771,-41.09482
),.Dim = c(2L,.Dimnames = list(c("x","y"),c("min","max"
))),proj4string = new("CRS",projargs = "+proj=merc +a=6378137 +b=6378137 +lat_ts=0 +lon_0=0 +x_0=0 +y_0=0 +k=1 +units=m +nadgrids=@null +wktext +no_defs"))

#main package for this analysis
library(adehabitatHR)  
library(sp)
library(readxl)
library(rgdal)

HRtracks<-read_excel("allinterpolatedtracks.xlsx") 
coordinates(HRtracks)<-c("x","y") 
proj4string(HRtracks)<- CRS("+init=epsg:3857")
#Organising variable 1 refers to Trip ID
colkud<-kernelUD(tracks.utm[,1],h="href",grid=1000,same4all=T)
#Organising variable 2 refers to Breeding stage
colkudb<-kernelUD(tracks.utm[,2],same4all=T)
kernel.area(colkud,percent=50)
kernel.area(colkudb,percent=50)
kernel.area(colkud,percent=95)
kernel.area(colkudb,percent=95)
colkud95<- getverticeshr(colkud,percent=95)
colkud95b<- getverticeshr(colkudb,percent=95)
colkud50b<- getverticeshr(colkudb,percent=50)
colkud50<- getverticeshr(colkud,percent=50)
#This is where the error occurs
writeOGR(colkud95,layer="colkud95",driver="ESRI Shapefile")

我没有找到答案的运气,我的目标是为家庭范围创建多边形,然后将它们导出为要在 QGIS 中使用的 shapefile。

解决方法

暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!

如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@)