La incertidumbre de una línea entre dos coordenadas con las incertidumbres en QGIS

Question

Estoy tratando de trazar la incertidumbre de una línea entre dos coordenadas con las incertidumbres. He creado polígonos alrededor de las 2 de coordenadas para representar a la región de incertidumbre alrededor de cada punto. Yo, a continuación, crear puntos aleatorios en cada polígono, y luego unirse a ellos mediante el uso de puntos de las rutas de la función. Se ve algo como esto:

Quiero trazar una región que representa la probabilidad de las líneas, y también elimina la incertidumbre de la región de la original extremos en cuenta - debe tener este aspecto;

La región de incertidumbre, debe estar en el 95% de nivel de confianza. Soy nuevo en QGIS, pensé que tal vez mediante la creación de una cuadrícula y rellenar las celdas dentro de la cuadrícula de alguna manera podría venir para arriba con un polígono que representa una densidad de probabilidad de la zona? Estoy usando QGIS Escritorio 3.0.1 con HIERBA 7.4.0. Creo que se necesita un gran número de puntos aleatorios/líneas para lograr lo que quieres, que puede hacer el procesamiento de cualquier algoritmos lento?

Answer 1

Que fue un divertido problema! He resuelto usando QGIS 2.18 pero no creo que las herramientas han cambiado drásticamente en 3.0.

1. Generar las líneas

No voy a dar mucho detalle aquí, ya que usted ya tiene las líneas que usted necesita. He escrito un script para generar estas líneas, pero funciona sólo hasta QGIS 2.18. Huelga decir que, las más de las líneas, el mejor de su estimación.

2. Crear una cuadrícula

Ir a Vector > Herramientas de Investigación > Vector de Cuadrícula (QGIS 2.18). En QGIS 3.0, parece que tiene que abrir el procesamiento (caja de herramientas de Procesamiento > caja de herramientas), a continuación, busque la herramienta de cuadrícula.

Para la cuadrícula de medida, seleccione "Usar capa/lienzo medida" (3 puntos a la derecha). Seleccione cualquier celda del tamaño que quieras, pero estar preparado para enfrentarse a los largos tiempos de cómputo si su red es demasiado delgada. Seleccione "Salida de la cuadrícula como polígonos" en "tipo de cuadrícula".

3. El conteo de líneas por celda de la cuadrícula

Esto es un poco complicado. Primero tienes que cortar sus líneas de manera que cada línea se divide en las células. Para ello, utilice la intersección de la herramienta (Vector > Herramientas de Geoprocesamiento > Intersección en tanto QGIS 2.18 y 3.0). Seleccione la capa de línea como "capa de Entrada", y la red como "Intersección de la capa".

Para preparar el siguiente paso, tenemos que agregar un campo a la capa que acabamos de obtener. Seleccione esta capa y abrir la calculadora de campo; marque la casilla "crear un nuevo campo de opción", establece el 'newid' como nombre, mantener 'numérico' como tipo. En la "Expresión" del cuadro de, escriba $rownum, a continuación, haga clic en ACEPTAR. Usted recibirá un nuevo campo con un identificador único para cada segmento de línea.

A continuación, vamos a utilizar la media de las coordenadas de la herramienta (Vector > herramientas de Análisis > coordenadas en QGIS 2.18, de nuevo en la caja de herramientas de QGIS 3.0). Establecer la interceptó capa (de la operación anterior) en la "capa de entrada", y "newid' como "campo IDENTIFICADOR Exclusivo". Esto le da una capa de puntos, cada punto corresponde a un segmento de línea.

Por último, el uso de la "Contar los puntos en la herramienta polígono" (en el Vector > herramientas de Análisis para ambos QGIS 2.18 y 3.0). Conjunto de la red en los "Polígonos" campo y el recién creado capa de punto en los "Puntos" del campo. Esto creará un nuevo campo en la cuadrícula con el número de segmentos de línea en cada celda de la cuadrícula. Se puede dividir este número por el número total de líneas en el campo de la calculadora si lo desea, si usted necesita una estimación de la línea de la probabilidad.

4. Ejemplo de resultados de

Obtenidos con más de 2000 líneas y un uniforme de punto de volver a crear particiones en el interior de la incertidumbre de las regiones (es posible que desee probar una gaussiana repartition demasiado).

La misma incertidumbre de los parámetros, con 5000 líneas y un diluyente de la cuadrícula (tomó varias horas para que se ejecute):