创建矩阵，显示许多多边形之间重叠的比例

Question

我有一个多边形列表，我需要知道多边形组合的重叠的比例。换句话此计算。我该怎么做？另外，我该提供什么使这篇文章有帮助？

这是我用来计算单个多边形列表重叠的代码的一个示例（处理）。是否可以调整它来创建每个poly计算的矩阵？

###################################################
#read in and create list of polygons
home.names <- list.files("", 
                         pattern="*.shp", recursive=T, full.names=T)
home.names


#####################################################
#create list to hold results and loop home-range overlap calculations
all.home.ranges <- list()

#loop
for (i in 1:length(home.names)){
  home.range <- st_read(home.names[i])
  home.range.t <- st_transform(home.range,'+proj=utm +zone=14 
        +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs')
  home.range$Area.ha <- st_area(home.range.t)*0.0001
  pi <- st_intersection(treatment, home.range.t) # treatment is a polygon
  if(nrow(pi) != 0){
    home.range$Treat.Overlap.per <- as.numeric(st_area(pi) / st_area(home.range.t) *100)
  } else {
    home.range$Treat.Overlap.per <- 0
  }

Answer 1

一种方法是构建具有2列的矩阵。例如，如果您有5个polys：

i <- rep(1:5, each=5)
j <- rep(1:5, 5)
comp <- cbind(i, j)[i < j, ]
comp
#       i j
#  [1,] 1 2
#  [2,] 1 3
#  [3,] 1 4
#  [4,] 1 5
#  [5,] 2 3
#  [6,] 2 4
#  [7,] 2 5
#  [8,] 3 4
#  [9,] 3 5
# [10,] 4 5

loop thun comp，并使用i和j选择两个多粒子进行比较。没有可重现的数据，我将无法进行调试，但是类似：

for (k in 1:nrow(comp){
  home.range <- st_read(home.names[comp[k, 1]])
  home.range.t <- st_transform(home.range,'+proj=utm +zone=14 
        +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs')
  home.range$Area.ha <- st_area(home.range.t)*0.0001

  home.range2 <- st_read(home.names[comp[k, 2]])
  home.range.t2 <- st_transform(home.range2,'+proj=utm +zone=14 
        +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs')
  home.range2$Area.ha <- st_area(home.range.t2)*0.0001

  pi <- st_intersection(home.range.t, home.range.t2)
  if(nrow(pi) != 0){
    home.range$Treat.Overlap.per <- as.numeric(st_area(pi) / st_area(home.range.t) *100)
  } else {
    home.range$Treat.Overlap.per <- 0
  }

主要区别是重叠的百分比可能基于home.range和home.range2。您可以根据计划对结果的计划进行计算。