Saltar al contenido principal

23  Práctica 2 · Sistemas de coordenadas y SIRGAS-Chile 2025

Duración: 3 horas · Software: QGIS 4.4 (o 3.44 LTR), Python (PyProj) · Lectura previa: Capítulo 6 y Capítulo 7 Datos: data/puntos_control_santiago.csv

23.1 Objetivos

  1. Transformar coordenadas entre PSAD56, SAD69, SIRGAS-Chile 2002 y SIRGAS-Chile 2025.
  2. Configurar y verificar los SRC UTM oficiales de SIRGAS-Chile 2025 (EPSG:11475 y 11476).
  3. Cuantificar las diferencias entre datums y explicar su origen.

23.2 Parte A · Conversión de formatos de coordenadas

Convierte a grados decimales y luego a UTM 19S las coordenadas de los puntos del archivo data/puntos_control_santiago.csv. Usa la función gms_a_decimal del Capítulo 7.

23.3 Parte B · SRC oficiales de SIRGAS-Chile 2025 en QGIS

  1. Proyecto → Propiedades → SRC y busca 11476 (SIRGAS-Chile 2025 / UTM zone 19S).
  2. Revisa su definición (WKT) y compárala con SIRGAS-Chile 2021 / UTM zone 19S (EPSG:20049): mismo elipsoide y proyección, distinto datum y época.
  3. Repite para la zona 18S (EPSG:11475) y anota el código de Rapa Nui (EPSG:11474, UTM 12S).
Nota¿Y si QGIS no encuentra el código?

Los códigos 11474–11476 están en la base EPSG desde 2026. Si tu instalación no los reconoce, actualiza QGIS/PROJ. No definas un SRC personalizado con +towgs84=0,...: ocultaría la época y el datum reales.

23.4 Parte C · Transformaciones con PyProj

Ejecuta el script scripts/p02_transformaciones.py:

import pandas as pd
from pyproj import Transformer

df = pd.read_csv("data/puntos_control_santiago.csv")

# 1) Geográficas SIRGAS-Chile 2002 -> UTM 19S SIRGAS-Chile 2002
t1 = Transformer.from_crs("EPSG:5360", "EPSG:5361", always_xy=True)
df["E_sirgas2002"], df["N_sirgas2002"] = t1.transform(df.lon.values, df.lat.values)

# 2) UTM 19S SIRGAS-Chile 2002 -> UTM 19S PSAD56 (dato heredado)
t2 = Transformer.from_crs("EPSG:5361", "EPSG:24879", always_xy=True)
df["E_psad56"], df["N_psad56"] = t2.transform(df.E_sirgas2002.values, df.N_sirgas2002.values)

df["dE"] = df.E_psad56 - df.E_sirgas2002
df["dN"] = df.N_psad56 - df.N_sirgas2002
print(df[["punto", "dE", "dN"]].round(1))

Observa el orden de magnitud de las diferencias entre PSAD56 y SIRGAS-Chile: en Santiago son cercanas a 180 m en Este y 370 m en Norte (unos 415 m en total) con la transformación por defecto de PROJ.

23.5 Parte D · Transformaciones disponibles

Lista las transformaciones que PROJ conoce entre dos sistemas y su precisión:

from pyproj.transformer import TransformerGroup
tg = TransformerGroup("EPSG:4248", "EPSG:5360")   # PSAD56 -> SIRGAS-Chile 2002
for t in tg.transformers:
    print(f"{t.description} | precisión: {t.accuracy} m")

23.6 Preguntas

  1. ¿Por qué las diferencias PSAD56–SIRGAS son de cientos de metros y las diferencias SIRGAS-Chile 2021–2025 de centímetros a decímetros?
  2. ¿Qué transformación elegirías para migrar una carta IGM 1:50.000 en PSAD56 y por qué su precisión es suficiente para esa escala?

23.7 Entregable

Tabla con las coordenadas en los cuatro sistemas, las diferencias calculadas y una respuesta argumentada a las preguntas.