Processing math: 100%

4 votos

Cómo definir un CRS usando una proyección girada sobre un punto dado

Tengo un conjunto de datos que utiliza una proyección normal, pero tiene su propia falsa coordenada y coordenada y, y se gira 315 grados. Necesito ser capaz de trabajar en este CRS y traducir fácilmente de un CRS a otro, por lo que el uso de hacks como herramientas CAD no es práctico. Yo tengo una definición de ella para MapInfo, que funciona bien. Sin embargo, la necesito para trabajar también en QGIS y otros programas de software, los cuales requieren .proj4 formato. He definido un CRS que da los resultados que quiero, pero al cambiar de una proyección a otra, hay un desplazamiento que genera es diferente dependiendo de la proyección que se utiliza.

El MapInfo definición es la siguiente: "GRDMF", 3018, 31, 7, 173, -41, 2510000, 6023150, 7, 0.686785, 0.726871, -5539726.41, -0.726876, 0.686796, -2108424.34, -1000000, -1000000, 1000000, 1000000

El .proj4 definición que tengo es: +proj=nzmg +lat_0=-41 +lon_0=173 +x_0=2510000 +y_0=6023150 +no_defs +a=6378388 +rf=297 +ellps=GRS80 +datum=nzgd49 +nadgrids=nzgd2kgrid0005.gsb +towgs84=84,-22,209,-0,-0,-0,0 +to_meter=1

La proyección es de Nueva Zelanda Mapa de Cuadrícula con NZGD1949. El CRS necesito se basa en NZMG de proyección, pero girada 315 acerca de la coordenada 2310655.80 5535100.29, que se traduce en 70501.32 13497.54 en la costumbre de CRS.

Una prueba de coordenadas son:

NZMG coords            Custom CRS
2292000 5539000        60522.871 29736.79
2296000 5552000        72720.155 35757.861
2330000 5527000        77898.623 -6126.353
2339000 5521000        79718.234 -16788.504
2333000 5511000        68329.450 -19294.963
2321000 5520000        66629.725 -4392.060

Cualquier ayuda se agradece. Gracias

2voto

Perdón por la demora en responder. Los parámetros para la transformación de su costumbre de CRS a EPSG:27200 parecen ser de la siguiente manera:

Escala x: 0.9999708285575817

Escala y: 0.9999269459420771

Rotación: 46.62502359527837

DX: 2223279.845

DY: 5516242.139

El problema es que no puedo hacer una custom de QGIS/Proj4 CRS que utiliza estos (es posible, pero no sé cómo hacerlo, la rotación me tiene perplejo).

He obtenido los parámetros de uso (gratis) OpenJump y el uso de esto sería posible afín a su cuadrícula personalizada shapefiles a EPSG:27200 bastante fácilmente. Hay un ejemplo de cómo hacer esto en el antiguo QGIS foro, consulte:

Enlace al hilo del foro

cerca de la parte inferior de la primera página.

Sería incluso mejor si alguien que es más brillante que me podría venir para arriba con una costumbre Proj4 cadena de parámetro :)

Saludos, Nick.

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X