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