除“href”之外,R 中内核 UD() 函数的合适“h”参数,以包含数据中的所有空间点

Suitable "h" parameter for the Kernel UD() function in R other than "href" to encompass all spatial points in the data

提问人:PrabhaV 提问时间:11/11/2023 更新时间:11/11/2023 访问量:9

问:

此脚本中使用的数据的示例图像#

这是我的 R 脚本,它导致了 50% 和 95% 占用率等值线的附加输出。

library(sp)  
library(adehabitatHR)  
library(raster)  
library(ggplot2)  

#read 包含数据点的纬度和经度的.csv中

Bat1 <- read.csv("D:/Bats/Bat1.csv", header = TRUE)  

#create 空间点对象

datapoints<-SpatialPoints(cbind(Bat1$Lat, Bat1$Lon))   

kud <- kernelUD(datapoints, h= "href", grid = 3000, extent = 5.0) 

image(kud)

hr50<-getverticeshr(kud, percent= 50, unout = "m2") 

hr95<-getverticeshr(kud, percent= 95, unout = "m2") 

设置坐标并写入可在 ArcGIS 中投影的 shapefile

WGScoor1<- hr50
proj4string(WGScoor1)<- CRS("+proj=longlat +ellips=WGS84")
raster::shapefile(WGScoor1, "D:/Bats/Bat1HR_50_new.shp")

WGScoor2<- hr95
proj4string(WGScoor2)<- CRS("+proj=longlat +ellips=WGS84")
raster::shapefile(WGScoor2, "D:/Bats/Bat1HR_95_new.shp")

我注意到我调整了多少“网格”、“范围”和“h”值,它并没有包含大多数空间点。我发现平滑参数“h”是要操作的参数,但给它“LSCV”或数字都无济于事。

我还能尝试什么?产生 50% 和 95% 的轮廓 shapefile

参数 kernel-density

评论


答: 暂无答案