[QGIS-it-user] QGIS line_interpolate_point con EPSG 4326

Ciao lista,
ho notato che la funzione line_interpolate_point non lavora come mi
aspetterei se la geometry è in EPSG 4326; con geometry proiettate funziona.

allego [0] uno shapefile linestring in EPSG 4326 e vi chiedo di fare la
seguente prova:
aggiungere un campo testuale -lungo 80 e popolatelo con la seguente
espressione:

geom_to_wkt( line_interpolate_point( $geometry, $length/2) )

poi esportate in csv per verificare i punti.

In QGIS 2.18.23 LTR il campo viene compilato ma le coordinate sono del primo
vertice e non quello atteso del centroide lungo la linea;
In QGIS 3.2.2 il campo viene compilato ma le coordinate sono del primo
vertice e non quello atteso del centroide lungo la linea;
In QGIS 3.3 dev il campo NON viene compilato, o meglio il valore è NULL

in tutte e tre le versioni l'espressione funziona bene se utilizzata per
tematizzare con geometry generator

domanda: è normale questo comportamento con EPSG 4326 oppure è un bug?

saluti

[0] https://mega.nz/#!sAwlyTLI!GmVZP4k2kiC914CutDD1r9D2VV327y3WryUGkL3FEkM

-----
https://pigrecoinfinito.wordpress.com/
--
Sent from: http://osgeo-org.1560.x6.nabble.com/QGIS-Italian-User-f5250612.html

Immagino che sia legato al fatto che non vengono usate sistemi di
riferimento metrici, ma in gradi.

direi che e' un bug, ma non ho fatto prove, meglio se qualcuni conferma
provando con numeri piccoli.

Saluti, e grazie.

On 08/27/2018 07:19 PM, Totò wrote:

Ciao lista,
ho notato che la funzione line_interpolate_point non lavora come mi
aspetterei se la geometry è in EPSG 4326; con geometry proiettate funziona.

allego [0] uno shapefile linestring in EPSG 4326 e vi chiedo di fare la
seguente prova:
aggiungere un campo testuale -lungo 80 e popolatelo con la seguente
espressione:

geom_to_wkt( line_interpolate_point( $geometry, $length/2) )

poi esportate in csv per verificare i punti.

In QGIS 2.18.23 LTR il campo viene compilato ma le coordinate sono del primo
vertice e non quello atteso del centroide lungo la linea;
In QGIS 3.2.2 il campo viene compilato ma le coordinate sono del primo
vertice e non quello atteso del centroide lungo la linea;
In QGIS 3.3 dev il campo NON viene compilato, o meglio il valore è NULL

in tutte e tre le versioni l'espressione funziona bene se utilizzata per
tematizzare con geometry generator

domanda: è normale questo comportamento con EPSG 4326 oppure è un bug?

saluti

[0] https://mega.nz/#!sAwlyTLI!GmVZP4k2kiC914CutDD1r9D2VV327y3WryUGkL3FEkM

-----
https://pigrecoinfinito.wordpress.com/
--
Sent from: http://osgeo-org.1560.x6.nabble.com/QGIS-Italian-User-f5250612.html
_______________________________________________
QGIS-it-user mailing list
QGIS-it-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/qgis-it-user

--
Paolo Cavallini - www.faunalia.eu
QGIS & PostGIS courses: http://www.faunalia.eu/training.html
https://www.google.com/trends/explore?date=all&geo=IT&q=qgis,arcgis

pcav wrote

direi che e' un bug, ma non ho fatto prove, meglio se qualcuni conferma
provando con numeri piccoli.

Saluti, e grazie.

Ho aperto ticket
https://issues.qgis.org/issues/19756

saluti

PS. cosi non andiamo da nessuna parte, nessuno collabora; non credo sia una
richiesta cosi difficile la mia.
pazienza.
saluti

-----
https://pigrecoinfinito.wordpress.com/
--
Sent from: http://osgeo-org.1560.x6.nabble.com/QGIS-Italian-User-f5250612.html

On 09/03/2018 02:49 PM, Totò wrote:

Ho aperto ticket
https://issues.qgis.org/issues/19756

saluti

PS. cosi non andiamo da nessuna parte, nessuno collabora; non credo sia una
richiesta cosi difficile la mia.

in che senso?

--
Paolo Cavallini - www.faunalia.eu
QGIS & PostGIS courses: http://www.faunalia.eu/training.html
https://www.google.com/trends/explore?date=all&geo=IT&q=qgis,arcgis

pcav wrote

On 09/03/2018 02:49 PM, Totò wrote:

Ho aperto ticket
https://issues.qgis.org/issues/19756

saluti

PS. cosi non andiamo da nessuna parte, nessuno collabora; non credo sia
una
richiesta cosi difficile la mia.

in che senso?

nel senso che ho chiesto di fare una semplice verifica, ho condiviso
shapefile e descritto cosa fare;
ma nessuno ha tempo o voglia... pazienza
saluti

-----
https://pigrecoinfinito.wordpress.com/
--
Sent from: http://osgeo-org.1560.x6.nabble.com/QGIS-Italian-User-f5250612.html

On 09/03/2018 04:42 PM, Totò wrote:

nel senso che ho chiesto di fare una semplice verifica, ho condiviso
shapefile e descritto cosa fare;
ma nessuno ha tempo o voglia... pazienza

è l regola: ascolatano in tanti, parlano in pochi, fanno in pochissimi.
per questo è giusto stimoleare, come fai tu, ma non usiamo toni
disfattisti :slight_smile:
Saluti, e grazie per quel che fai.

--
Paolo Cavallini - www.faunalia.eu
QGIS & PostGIS courses: http://www.faunalia.eu/training.html
https://www.google.com/trends/explore?date=all&geo=IT&q=qgis,arcgis

Totò wrote

Ciao lista,
ho notato che la funzione line_interpolate_point non lavora come mi
aspetterei se la geometry è in EPSG 4326; con geometry proiettate
funziona.

allego [0] uno shapefile linestring in EPSG 4326 e vi chiedo di fare la
seguente prova:
aggiungere un campo testuale -lungo 80 e popolatelo con la seguente
espressione:

geom_to_wkt( line_interpolate_point( $geometry, $length/2) )

poi esportate in csv per verificare i punti.

Nonostante tutto descrivo la soluzione del problema.

Nyall ha osservato, giustamente, che la funzione $length dipende dalle
impostazioni di progetto, se tali impostazioni vengono lasciate di default
la $length restituirà sempre un valore in metri (nel caso mio converte la
lunghezza in metri); questo è il motivo per cui la mia espressione non
funzionava.

Soluzione proposta da Nyall usare la seguente espressione

line_interpolate_point( $geometry, length($geometry)/2);

oppure, aggiungo io, usare sempre la prima espressione
line_interpolate_point( $geometry, $length/2)
e modificare le impostazioni di progetto, nel caso specifico:
1. usare None/Planimetrico nell'ellissoide
2. gradi nell'unità per il calcolo delle lunghezze

saluti

-----
https://pigrecoinfinito.wordpress.com/
--
Sent from: http://osgeo-org.1560.x6.nabble.com/QGIS-Italian-User-f5250612.html