¿Usa el Sistema de coordenadas personalizado en ArcGIS Desktop?

16

No sé mucho sobre sistemas de coordenadas ...

En mi oficina utilizamos para tratar datos espaciales procedentes de sitios arqueológicos. Cada sitio tiene su propio sistema de coordenadas xyz (GCS). Tres ejes cartesianos ortogonales simples . En los últimos años hemos estado administrando estos datos espaciales a través del software GIS (ArcGIS), sin utilizar un sistema de coordenadas específico (solo déjelo como "indefinido")

Me gustaría saber si existe algún GCS diseñado para tratar dichos conjuntos de datos utilizando un eje ortogonal cartesiano simple, sin distorsiones de la cuadrícula del GCS típico. Además, me gustaría saber si este sistema es adecuado para usarlo en una aplicación de mapeo en línea.

Por cierto, gestionamos entornos 2D (ArcMap) y 3D (ArcScene) y trabajamos con "mm" como unidad base de longitud.

Si tal cosa no existe, tal vez alguien sepa cómo crearla.

Maik
fuente

Respuestas:

17

Suponiendo que sus áreas de interés son bastante pequeñas en comparación con el mundo, puede configurar una proyección mercator transversal personalizada.

Debe conocer las coordenadas geográficas lat_0 y lon_0 del origen de su CRS, y la dirección de los ejes x e y:

+proj=tmerc +lat_0=51.4 +lon_0=7 +k=1 +x_0=0 +y_0=0 +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs

Para usar con Arcgis, el archivo .prj debe contener:

PROJCS["Transverse_Mercator",GEOGCS["GCS_WGS_1984",DATUM["D_unknown",SPHEROID["WGS84",6378137,298.257223563]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",51.4],PARAMETER["central_meridian",7],PARAMETER["scale_factor",1],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["Meter",1]]

Si no están en paralelo con las direcciones Norte y Este, sería necesaria una rotación utilizando una proyección oblicua de mercator:

+proj=omerc +lat_0=51.4 +lonc=7 +alpha=-10 +k=1 +x_0=0 +y_0=0 +gamma=0 +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs

El archivo .prj para esta proyección es:

PROJCS["Hotine_Oblique_Mercator_Azimuth_Center",GEOGCS["GCS_WGS_1984",DATUM["D_unknown",SPHEROID["WGS84",6378137,298.257223563]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Hotine_Oblique_Mercator_Azimuth_Center"],PARAMETER["latitude_of_center",51.4],PARAMETER["longitude_of_center",7],PARAMETER["azimuth",-10],PARAMETER["scale_factor",1],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["Meter",1]]

Aquí están ambos juntos (tmerc en azul y omerc en rojo):

ingrese la descripción de la imagen aquí

Lat_0 y lon_0 (o lonc) son las coordenadas de su origen local de CRS. Tienes que explorarlos con un GPS o desde otros mapas. El ángulo que debe calcular a partir de las coordenadas de dos puntos conocidos en su eje X local.

Lo único que no funciona como se espera es el milímetro como unidades. Debería funcionar, pero no funciona en QGIS actual :-(

http://hub.qgis.org/issues/9414


Actualizar

Las versiones actuales de QGIS ahora pueden trabajar con proyecciones locales utilizando mm como unidades.

AndreJ
fuente