Añadir polilíneas de forma programada a partir de puntos existentes en QGIS

Question

Estoy usando QGIS 2.18 y necesito ayuda para encontrar un python script o una herramienta que me permita dibujar programáticamente una línea desde los puntos GPS existentes (a veces hasta 500 o más) que tenga una longitud de 250 pies en un rumbo listado en un campo de atributo de acimut para cada punto.

Answer 1

Si no necesita realmente una geometría y sólo necesita dibujar las líneas como se describe, puede utilizar el constructor de geometría de QGIS.

Si necesitas las geometrías, tengo otra respuesta en esta misma pregunta que muestra cómo generar las geometrías usando un python script aquí .

Dentro del editor de expresiones necesitarías algo así

make_line($geometry, make_point($x + 250*cos(radians("azimuth")), $y + 250*sin(radians("azimuth"))))

Desglosando...

make_line toma dos puntos y crea una línea

$geometry la geometría actual. Si tiene una capa de puntos, este es el punto actual

make_point crea un punto con una x y una y dadas

$x el valor x de la geometría actual.

$y el valor y de la geometría actual.

250*cos(radians("azimuth")) Desplazamiento x calculado a partir del acimut. Esto supone que está utilizando un sistema de coordenadas proyectadas, la distancia es de 250 unidades (cualquiera que sean las unidades de su sistema de coordenadas. Para la mayoría serían metros), y el acimut se almacena en un atributo llamado acimut y se almacena en grados, y la forma en que está escrito es utilizando un sistema de coordenadas matemáticas, por lo que 0 grados es el Este aumentando en sentido contrario a las agujas del reloj, por lo que el Norte es 90, el Oeste es 180, etc.

A continuación se muestra un ejemplo de salida.

Answer 2

Pude modificar mi script para crear las geometrías que buscabas hacer. La distancia debe estar en las unidades de tu sistema de coordenadas. Yo sugeriría transformar a un sistema de referencia de coordenadas proyectado apropiado que utilice pies o metros.

Para utilizar este código, abra la caja de herramientas de Processing, luego expanda Scripts y luego expanda Herramientas. Seleccione "Crear nuevo script" y copie y pegue el código de abajo en la ventana script (tenga cuidado al copiar y pegar código python ya que los espacios en blanco son sintácticamente significativos. Si tienes problemas pon el código en un editor de texto que muestre los espacios en blanco y asegúrate de que se ha copiado correctamente). Guárdalo donde quieras y hay un botón de ejecutar script en la parte superior de la ventana. Después de guardarlo puedes "Añadir script desde el archivo" y tener permanentemente el script bajo "Scripts de usuario".

##input_layer=vector
##distance=number
##output_layer=output vector

from qgis.core import QgsWKBTypes, QgsField, QgsVectorFileWriter, QgsFeature, QgsGeometry, QgsPoint
from PyQt4.QtCore import QVariant

def addLines( geometry, writer, inFeature ):
    coordinateSequence = geometry.coordinateSequence()
    for rings in coordinateSequence:
        for points in rings:
            for point in points:
                feature = QgsFeature( fields )
                azimuth = inFeature.attribute('azimuth')
                simplePoint = QgsPoint(point.x(), point.y())
                pointAtAngle = simplePoint.project(distance, azimuth)
                lineGeometry = QgsGeometry.fromPolyline([simplePoint, pointAtAngle])
                feature.setGeometry( QgsGeometry( lineGeometry ) )
                type = lineGeometry.wkbType()
                attributes = inFeature.attributes()
                feature.setAttributes( attributes )
                writer.addFeature( feature )
    return

inlayer = processing.getObject( input_layer )
provider = inlayer.dataProvider()
fields = provider.fields()
geomType = QgsWKBTypes.Type(inlayer.wkbType())
outputGeomType = QgsWKBTypes.MultiLineString

writer = QgsVectorFileWriter( output_layer, 'UTF-8', fields,  outputGeomType , inlayer.crs() )

features = inlayer.getFeatures()
featureCount = inlayer.featureCount()
featureIndex = 0

for f in features:
    percent = ( featureIndex/float( featureCount ) ) * 100
    progress.setPercentage( percent )
    g = f.geometry().geometry()
    addLines( g, writer, f )
    featureIndex +=1

del writer

Answer 3

Busque el uso de https://plugins.qgis.org/plugins/qchainage/

Una vez que instales el plugin deberías poder crear esas líneas.

Este plugin toma las características de las líneas y crea una nueva capa de puntos en distancias proporcionadas sobre estas Líneas. Opcionalmente se puede establecer el punto inicial (distancia desde el inicio de la línea) y el punto final (parada antes de llegar al final de la línea) y etiquetar automáticamente los nuevos puntos con las distancias.

Answer 4

Si no quieres usar expresiones puedes usar el renderizador de campos vectoriales para tus propósitos. He creado un conjunto de datos de muestra de cinco puntos con un bearing atributo...

...y con la calculadora de campos añadió un atributo length que define como 2500m para cada característica (puede utilizar diferentes valores de longitud para cada uno o utilizar la velocidad del GPS para la longitud):

Luego en Layer properties | Point Style añadí una capa de símbolos de tipo 'vector field renderer' en la que se pueden definir atributos para la longitud y el ángulo:

lo que resulta en