Estoy tratando de transformar puntos tomados de los mapas de Google (como latitud y longitud) a Irish Grid TM65 utilizando el spTransform en R, y estoy confundido en algo.
Todo lo que encuentro en la red me dice que google usando una proyección mercator como tal:
"+proj=merc +lon_0=0 +k=1 +x_0=0 +y_0=0 +a=6378137 +b=6378137 +towgs84=0,0,0,0,0,0 +units=m +no_defs"
Para la cuadrícula irlandesa estoy utilizando una cadena proj4 extraída de los shapefiles del censo irlandés
"+proj=tmerc +lat_0=53.5 +lon_0=-8 +k=1.000035 +x_0=200000 +y_0=250000 +date=ire65 +units=m +no_defs +ellps=mod_airy +towgs84=482.530,-130.596,564.557,-1.042,-0.214,-0.631,8.15"
Así que eso está bien - he enchufado todo eso en R, pero estaba recibiendo respuestas extrañas. Así que me puse a jugar con diferentes CRS en QGIS para ver lo que estaba pasando. De todos modos encontré que en lugar de utilizar la proyección anterior para google, en lugar de utilizar esta proyección para el CRS para los puntos tomados de googlemaps entonces me dio proyecciones correctas en mi CRS Irish Grid: "+proj=longlat +datum=WGS84 +no_defs"
¿Cómo es que la proyección que todo Internet me dijo que usara me dio resultados extraños y la que elegí casi al azar funcionó?
Por cierto, esto funcionó tanto en QGIS como en R. Estoy razonablemente feliz con los resultados - pero claramente me estoy perdiendo algo y tengo que entender esto correctamente para el trabajo que estoy haciendo
