5。域

Question

对不起，关于我的问题。我仍在学习如何在R中进行空间分析。我知道以前曾有人问过这个问题（ here ）。

目标：我无法在专门针对我的经度（X）值的参数内使用模拟数据运行此代码（请参见第24行）。我想用模拟数据绘制家庭范围（如下）。

错误：“ getverticeshr.estUD（x [[i]]中的错误，百分比，ida =名称（x）[i]，unin ，：网格太小，无法进行估算您应该使用更大的范围参数重新运行kernelUD。

# 1. Packages
library(adehabitatHR)         # Package for spatal analysis

# 2. Empty Dataframe
points <- data.frame(ID = double())
XY_cor <- data.frame(X = double(),
                 Y = double())
# 3. Assigning values (this will be our spatial coordinates)
set.seed(17)
for(i in c(1:100)){
  if(i >= 50){points[i, 1] <- 1}
  else {points[i, 1] <- 2}

  XY_cor[i, 1] <- runif(1, -78.86887, -78.86440) ## error is here!
  XY_cor[i, 2] <- runif(1, 0.958533, 0.960777)} 

# 4. Transform to SpatialDataframe
coordinates(points) <- XY_cor[, c("X", "Y")]
class(points)

# 5. Domain
x <- seq(-80.0, -77.0, by=1.) # resolution is the pixel size you desire 
y <- seq(-200, 200, by=1.)
xy <- expand.grid(x=x,y=y)
coordinates(xy) <- ~x+y
gridded(xy) <- TRUE
class(xy)

# 6. Kernel Density
kud_points <- kernelUD(points, h = "href", grid = xy)
image(kud_points)

# 7. Get the Volum
 vud_points <- getvolumeUD(kud_points)

# 8. Get contour
levels <- c(50, 75, 95)
list <- vector(mode="list", length = 2)

list[[1]] <- as.image.SpatialGridDataFrame(vud_points[[1]])
list[[2]] <- as.image.SpatialGridDataFrame(vud_points[[2]])

# 9. Plot
par(mfrow = c(2, 1))
image(vud_points[[1]])
contour(list[[1]], add=TRUE, levels=levels)
image(vud_points[[2]])
contour(list[[2]], add=TRUE, levels=levels)


# 10. Get vertices 
vkde_points <- getverticeshr(kud_points, percent = 50,
                         unin = 'm', unout='m2')
plot(vkde_points)

Answer 1

我可以使用以下代码解决它。希望这可以帮助其他人解决同样的问题！我必须手动更改边界框的分辨率。

5。域

x <- seq(737326.0, 737639.6, by=5)        
y <-seq(106071.4, 106259.9, by=5)
xy <- expand.grid(x=x,y=y)
coordinates(xy) <- ~x+y    
gridded(xy) <- TRUE

R代码adehabitatHR-网格对于内核UD / getverticeshr / adehabitatHR主范围估计而言太小

1 个答案:

5。域