Generación de puntos dentro de un polígono con QGIS

Question

Estoy utilizando QGIS.

Tengo una característica poligonal y quiero poder generar puntos dentro de ella. Necesito esto para una tarea de clasificación.

Generar puntos aleatorios hasta que uno esté dentro del polígono no funcionaría porque es realmente impredecible el tiempo que tarda.

Answer 1

Comience por descomponer el polígono en triángulos y, a continuación generar puntos dentro de esos . (Para una distribución uniforme, ponderar cada triángulo por su área).

Answer 2

La herramienta Puntos aleatorios puede utilizarse con una capa límite.

Si está buscando código, el código fuente del plugin subyacente debería serle de ayuda.

Answer 3

Hay algunas buenas bibliotecas que hacen la mayor parte del trabajo pesado por ti.

Ejemplo de uso de [shapely][1] en python.

import random
from shapely.geometry import Polygon, Point

def get_random_point_in_polygon(poly):
     minx, miny, maxx, maxy = poly.bounds
     while True:
         p = Point(random.uniform(minx, maxx), random.uniform(miny, maxy))
         if poly.contains(p):
             return p

p = Polygon([(0, 0), (0, 2), (1, 1), (2, 2), (2, 0), (1, 1), (0, 0)])
point_in_poly = get_random_point_in_polygon(mypoly)

O utilice .representative_point() para obtener un punto dentro del objeto (como menciona dain):

Devuelve un punto calculado a bajo coste que se garantiza que está dentro del objeto geométrico.

poly.representative_point().wkt
'POINT (-1.5000000000000000 0.0000000000000000)'

  [1]: https://shapely.readthedocs.io

Answer 4

Podría determinar la extensión del polígono y luego restringir la generación de números aleatorios para los valores X e Y dentro de esas extensiones.

Proceso básico: 1) Determinar maxx, maxy, minx, miny de los vértices del polígono, 2) Generar puntos aleatorios utilizando estos valores como límites 3) Comprobar la intersección de cada punto con el polígono, 4) Deje de generar cuando usted tiene suficientes puntos que satisfagan la prueba de intersección

He aquí un algoritmo (C#) para la prueba de intersección:

bool PointIsInGeometry(PointCollection points, MapPoint point)
{
int i;
int j = points.Count - 1;
bool output = false;

for (i = 0; i < points.Count; i++)
{
    if (points[i].X < point.X && points[j].X >= point.X || points[j].X < point.X && points[i].X >= point.X)
    {
        if (points[i].Y + (point.X - points[i].X) / (points[j].X - points[i].X) * (points[j].Y - points[i].Y) < point.Y)
        {
            output = !output;
        }
    }
    j = i;
}
return output;
}

Answer 5

Si R es una opción, véase ?spsample en el sp paquete. Los polígonos se pueden leer desde cualquier formato compatible con GDAL incorporado en el paquete rgdal, y luego spsample trabaja directamente sobre el objeto importado con diversas opciones de muestreo.