Разделите многоугольник неправильной формы на 3 равные части - PullRequest
Новая
       93

Разделите многоугольник неправильной формы на 3 равные части

0 голосов
/ 10 ноября 2018

У меня есть многоугольник объекта линии разреза, который я хотел бы разделить на 3 равные части, расположенные один за другим вдоль линии разреза. До сих пор мне удавалось разрезать многоугольник только по вертикали, что приводит к образованию не непрерывных частей.

Вот воспроизводимый пример с использованием данных Гарвардского леса, доступных по следующей ссылке

После того, как вы скачали и распаковали данные, вы увидите папку с именем HARV. В этой папке находится шейп-файл, HARV_roads.shp

Import HARV_roads.shp 

lines_HARV <- readOGR("/Your file path/NEON-DS-Site-Layout-Files/HARV/HARV_roads.shp")

Получить только линии пешеходной дорожки 

footpath_HARV <- lines_HARV[lines_HARV$TYPE == "footpath",]

Участок тропинки 

plot(footpath_HARV,
     lwd=6,
     main="NEON Harvard Forest Field Site\n Footpath")

footpath

Создание более толстого многоугольного объекта из линии пешеходной дорожки с использованием буфера для имитации линии разреза с шириной X 

require(raster)
footpath_buffer <- buffer(footpath_HARV,width=40)

Участок footpath_buffer 

plot(footpath_buffer)

buffer

Разделите буфер на 3 равные части, используя код Барри Роулингсона (ниже) из этого поста : 

makeVchopper <- function(pol){
  bb = bbox(pol)
  delta = (bb[2,2] - bb[2,1])/10
  xmin = bb[1,1]-delta
  ymin = bb[2,1]-delta
  ymax = bb[2,2]+delta

  choppoly = function(xmax){
    readWKT(sprintf("POLYGON((%s %s, %s %s, %s %s, %s %s, %s %s))",
                    xmin,ymin, xmin,ymax, xmax,ymax, xmax,ymin,
                    xmin,ymin))
  }
  choppoly
}

slicer <- function(pol, xmin, xmax){
  bb = bbox(pol)
  delta = (bb[2,2] - bb[2,1])/10
  ymax = bb[2,2] + delta
  ymin = bb[2,1] - delta
  r = readWKT(sprintf("POLYGON((%s %s, %s %s, %s %s, %s %s, %s %s))",
                      xmin,ymin, xmin,ymax, xmax,ymax, xmax,ymin, xmin,ymin))
  gIntersection(pol,r)
}

chop_thirds <- function(pol, fractions=c(1/3, 2/3)){
  chopper = makeVchopper(pol)
  bb = bbox(pol)
  xmin = bb[1,1]
  xmax = bb[1,2]

  totalArea = gArea(pol)

  chopped_area = function(x){
    gArea(gIntersection(chopper(x),pol))
  }

  edges = lapply(fractions, function(fraction){
    target = totalArea * fraction
    target_function = function(x){
      chopped_area(x) - target
    }
    uniroot(target_function, lower=xmin, upper=xmax)$root
  })

  xdelta = (xmax-xmin)/10
  chops = matrix(c(xmin-xdelta, rep(edges,rep(2,length(edges))),
                   xmax+xdelta), ncol=2, byrow=TRUE)
  apply(chops, 1, function(edges){
    slicer(pol, edges[1], edges[2])
  })

}

Разделите footpath_buffer на 3 равные части: 

parts.footpath <- chop_thirds(footpath_buffer)

Участок footpath_buffer и 3 равные части 

plot(footpath_buffer)
plot(parts.footpath[[1]], add=TRUE, col=1,border=F)
plot(parts.footpath[[2]], add=TRUE, col=2,border=F)
plot(parts.footpath[[3]], add=TRUE, col=3,border=F)

3 parts

footpath_buffer теперь разделен на 3 равные части, но части разбиты. Я надеюсь достичь чего-то подобного, где каждая часть непрерывна вдоль разреза:

3 continuous parts transect

...