Как создать случайные точки вне полигонов?

Инструмент Создать случайную точку способен генерировать определенное количество точек внутри полигонов. Мне интересно, учитывая ограничивающий прямоугольник, есть ли способ, которым я могу генерировать случайные точки за пределами этого многоугольника?

Лично мне не нравится алгоритм случайных точек в ArcGIS. В качестве альтернативы используйте функцию Geospatial Modeling Environment (GME) genrandompnts. Вы сможете идентифицировать определенные полигоны, где случайные точки будут исключены (см. Выделенную область в прикрепленном .jpg). Лучше всего это программное обеспечение бесплатно.

GME предоставляет вам набор инструментов для анализа и моделирования, начиная от небольших «строительных блоков», которые вы можете использовать для создания сложного рабочего процесса, до полностью автономных программ анализа. Он также использует чрезвычайно мощное программное обеспечение с открытым исходным кодом R в качестве статистического механизма для управления некоторыми инструментами анализа. Одной из сильных сторон R является то, что она имеет открытый исходный код, полностью прозрачна и хорошо документирована: важные характеристики для любого научно-аналитического программного обеспечения.

Вам нужно будет создать кольцевой многоугольник с дыркой от бублика, представляющей внутреннее неточечное пространство, и некоторый пространственный экстент, представляющий ограничивающую область многоугольной области.

Извините, я просто не могу устоять. Поскольку всегда полезно знать, что происходит под капотом с чем-то вроде GME, вот решение в реальном коде R.

require(sp)
require(rgeos)

# Create example polygon data
x <- readWKT("POLYGON ((0 0, 0 10, 10 10, 10 0, 0 0))")
y <- readWKT("POLYGON ((3 3, 7 3, 7 7, 3 7, 3 3))")

# Calculate difference in polygon geometries to create null polygon
d <- gDifference(x,y)

# Create random sample in non-null polygon
rs <- spsample(d, 20, type="random")  

# Plot results  
plot(d, col="red")
  plot(rs,pch=19,col="black",add=TRUE)

Этот подход, вероятно, сильно отличается от того, как это делает GME, но использует собственные пространственные классы R sp и довольно новую библиотеку топологий, что делает код очень эффективным. Это также дает пример, который можно легко обернуть в цикл for.