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

Я хотел бы разрезать полигональный слой в соответствии с высотой на две части (верхнюю и нижнюю). Многоугольник может быть выпуклым или вогнутым, а положение разреза может отличаться друг от друга. Линия контура имеет интервал 5 м, что означает, что мне может потребоваться создать контур с сильно сжатыми линиями контура, например, с интервалом 1 м. Есть идеи, как это сделать, лучше в ArcGIS или R? Многоугольник может быть выпуклым или вогнутым, а положение разреза может отличаться друг от друга. The ...

Я хотел бы разрезать полигональный слой в соответствии с высотой на две части (верхнюю и нижнюю). Многоугольник может быть выпуклым или вогнутым, а положение разреза может отличаться друг от друга. Линия контура имеет интервал 5 м, что означает, что мне может потребоваться создать контур с сильно сжатыми линиями контура, например, с интервалом 1 м. Есть идеи, как это сделать, лучше в ArcGIS или R? Многоугольник может быть выпуклым или вогнутым, а положение разреза может отличаться друг от друга. The ...

Я хотел бы разрезать полигональный слой в соответствии с высотой на две части (верхнюю и нижнюю). Многоугольник может быть выпуклым или вогнутым, а положение разреза может отличаться друг от друга. Линия контура имеет интервал 5 м, что означает, что мне может потребоваться создать контур с сильно сжатыми линиями контура, например, с интервалом 1 м. Есть идеи, как это сделать, лучше в ArcGIS или R? Это означает, что мне может потребоваться создать контур с сильно сжатыми линиями контура, например, с интервалом 1 м. Есть идеи, как это сделать, лучше в ArcGIS или R? Это означает, что мне может потребоваться создать контур с сильно сжатыми линиями контура, например, с интервалом 1 м. Есть идеи, как это сделать, лучше в ArcGIS или R? Ниже приведен рабочий пример для Q:

library(sp)
library(raster)
r<-raster(ncol=100,nrow=100)
values(r)<-rep(1:100,100)
plot(r)   ### I have no idea why half of the value is negative...
p1<-cbind(c(-100,-90,-50,-100),c(60,70,30,30,60))
p2<-cbind(c(0,50,100,0),c(0,-25,10,0))
p1p<-Polygons(list(Polygon(p1,hole=T)),"p1")
p2p<-Polygons(list(Polygon(p2,hole=T)),"p2")
p<-SpatialPolygons(list(p1p,p2p),1:2)
plot(p,add=T)
segments(-90,80,-90,20)  ##where the polygon could be devided
segments(50,20,50,-30)  ##

Заранее спасибо ~

Марко