¿Cómo obtener la latitud/longitud en grados/minuto segundo (DMS) usando la calculadora de campo QGIS?

Question

He ideado la siguiente expresión para convertir 5.1234 a 5° 7`24.24"

toint ( $x) || '° ' || toint((($ x) - toint ( $x)) * 60) ||'`' || substr( (tostring(((($ x) - toint ( $x)) * 60) - toint((($ x) - toint ($x)) * 60)) * 60),1,5) || '"'

Este es el problema: para algunos puntos como el 5.1234, funciona. pero para otros, no. Sospecho que el problema es la conversión de números enteros que redondean las cifras decimales en lugar de truncarlas.

¿Hay alguna otra opción? Gracias.

Answer 1

Actualmente existe una función nativa to_dms() para esta intención incluida en la calculadora de expresiones.

Ahora, puedes simplemente ejecutar to_dms(5.1234, 'x', 2) . A continuación se puede ver una captura del documento actual

Answer 2

Me gusta la expresión que has puesto - probablemente no hay solución en QGIS 1.8, pero en QGIS 1.9-dev hay un floor() función de redondeo que redondea abajo . Para D°M'S'':

(CASE WHEN $x < 0 THEN '-' ELSE '' END) || floor (abs($x)) || '° ' || floor(((abs($x)) - floor (abs($x))) * 60) ||'\'' || substr( (tostring((((abs($x)) - floor (abs($x))) * 60) - floor(((abs($x)) - floor (abs($x))) * 60)) * 60),1,5) || '"'

Tenga en cuenta el apóstrofe escapado ( \' ).

Para D°M.MMM':

(CASE WHEN $x < 0 THEN '-' ELSE '' END) || floor (abs($x)) || '° ' || ((abs($x) - floor(abs($x))) * 60) ||'\''

Para limitar el número de minutos decimales mostrados, sustituya # en la expresión siguiente por el número de dígitos:

(CASE WHEN $x < 0 THEN '-' ELSE '' END) || floor (abs($x)) || '° ' ||  format_number( (abs($x) - floor(abs($x))) * 60, #) ||'\'

Answer 3

Se puede utilizar el operador módulo para hacer el truncamiento, pero la expresión resultante sería muy fea. Es más bonito usar la sustitución de cadenas, pero desafortunadamente QGIS no expone ninguna función strpos o similar. Utilice regexp_replace($x, '\\..*', '') para obtener toda la parte y regexp_replace($x, '^[0-9]*\\.', '') para obtener la parte decimal. Utilice toreal en lugar de toint para los cálculos con la segunda expresión para asegurarse de que no habrá redondeo.

Answer 4

He hecho algunos cambios en el código para obtener las coordenadas en DMS, si es necesario. Sólo se necesitan ajustes en la información EPSG.

Latitud

(CASE WHEN $y < 0 THEN '-' ELSE '' END) || floor (abs(Y( transform( $geometry, 'EPSG:31983', 'EPSG:4326'))
)) || '°' || floor(((abs(Y( transform( $geometry, 'EPSG:31983', 'EPSG:4326'))
)) - floor (abs(Y( transform( $geometry, 'EPSG:31983', 'EPSG:4326'))
))) * 60) ||'\'' || substr( (tostring((((abs(Y( transform( $geometry, 'EPSG:31983', 'EPSG:4326'))
)) - floor (abs(Y( transform( $geometry, 'EPSG:31983', 'EPSG:4326'))
))) * 60) - floor(((abs(Y( transform( $geometry, 'EPSG:31983', 'EPSG:4326'))
)) - floor (abs(Y( transform( $geometry, 'EPSG:31983', 'EPSG:4326'))
))) * 60)) * 60),1,5) || '"'

Longitud

(CASE WHEN $x < 0 THEN '-' ELSE '' END) || floor (abs(X( transform( $geometry, 'EPSG:31983', 'EPSG:4326'))
)) || '°' || floor(((abs(X( transform( $geometry, 'EPSG:31983', 'EPSG:4326'))
)) - floor (abs(X( transform( $geometry, 'EPSG:31983', 'EPSG:4326'))
))) * 60) ||'\'' || substr( (tostring((((abs(X( transform( $geometry, 'EPSG:31983', 'EPSG:4326'))
)) - floor (abs(X( transform( $geometry, 'EPSG:31983', 'EPSG:4326'))
))) * 60) - floor(((abs(X( transform( $geometry, 'EPSG:31983', 'EPSG:4326'))
)) - floor (abs(X( transform( $geometry, 'EPSG:31983', 'EPSG:4326'))
))) * 60)) * 60),1,5) || '"'