¿Utilizar transform() en el editor de funciones de QGIS?

Question

Estoy intentando transformar las coordenadas de una característica como parte de un script para rellenar un campo de atributo calculado.

Esto funciona en el cuadro de expresión:

x(transform($geometry,'EPSG:4326','EPSG:3857'))

Así es como intenté implementarlo en el editor de funciones:

def gridLocation(value1, feature, parent):

    place = feature.geometry()
    tr = QgsCoordinateTransform('EPSG:4326','EPSG:3857')
    place.transform(tr)

    x_val = place.asPoint().x()
    # ...

pero x_val sigue devolviendo el valor x original.

¿Cuál es la forma correcta de utilizar transform?

Answer 1

Pocas cosas podrías cambiar:

• Añadir usesgeometry=True como otro parámetro de su función.
• Puede eliminar value1 ya que en realidad no estás usando esto.
• En su tr = QgsCoordinateTransform('EPSG:4326','EPSG:3857') línea, sustituya EPSG:4326 con QgsCoordinateReferenceSystem(4326) y repetir con EPSG:3857 .

---

A mí me funciona lo siguiente:

from qgis.core import *
from qgis.gui import *

@qgsfunction(args='auto', group='Custom', usesgeometry=True)
def gridLocation(feature, parent):
    place = feature.geometry()
    tr = QgsCoordinateTransform(QgsCoordinateReferenceSystem(4326),QgsCoordinateReferenceSystem(3857))
    place.transform(tr)
    x_val = place.asPoint().x()

    #Add return statement to return values to the attribute table
    return x_val