Georreferenciación CRS y QGIS para la proyección Cassini

Question

Tengo escaneados mapas de la serie OS 6" que entiendo son de proyección Cassini. Según un artículo de Mugnier, Photographic & remote sensing, octubre de 2003, tenían el típico factor de escala igual a la unidad, ningún origen falso y una única estación de triangulación como origen de la proyección. De esto deduzco que la cadena proj4 debería ser algo parecido a:

 +proj=cass +lat_0=53.9619966389 +lon_0=-1.0804806111 +x_0=0 +y_0=0 +units=m +no_defs

donde los valores lat_0 y lon_0 son las coordenadas del origen, en este caso York Minster. Si luego utilizo el georreferenciador QGIS con una transformación lineal con puntos de control tomados de las coordenadas impresas de las esquinas de la entrada y la salida establecidos en este CRS espero que genere un archivo de mundo para la exploración con las coordenadas son la parte superior izquierda de la hoja, x & y rotaciones de cero y x e y los valores de píxel que son iguales, excepto que y es negativo. De hecho, los valores de los píxeles son muy diferentes, de modo que cuando QGIS abre el archivo las proporciones están muy distorsionadas. Es evidente que algo en mi comprensión de esto está mal, pero ¿qué?

Se trata de puntos muy cercanos a las esquinas del mapa impreso. El archivo escaneado real es algo mayor, ya que incluye toda la hoja. Las proyecciones no se imprimen. La serie de condados de 6" tenía diferentes proyecciones para cada condado o, creo, en algunos casos para pequeños grupos de condados. En el caso de Yorkshire, se centraba en York Minster. Esta hoja en concreto es la 272.

-1.72916666666666652,53.72888888888888914
-1.87361111111111112,53.52027777777777828
-1.72833333333333328,53.52111111111111086

Answer 1

Las coordenadas de tus esquinas están en grados. Si desea el resultado en la proyección Cassini, tiene que introducir las coordenadas Cassini (en metros, relativas al punto de anclaje). Puede calcularlas con cs2cs.

Tomé la imagen de http://maps.nls.uk/view/102345139 y leer esas coordenadas de las esquinas (convertidas a grados decimales):

-1.8747 53.5781
-1.7290 53.5789
-1.8737 53.5202
-1.7285 53.5211

Póngalos en un archivo llamado Yorkdeg.txt, y escriba un archivo por lotes con este contenido:

cs2cs +init=epsg:4277 +to +proj=cass +ellps=airy +lat_0=53.9619966389 +lon_0=-1.0804806111 +x_0=0 +y_0=0 +towgs84=446.448,-125.157,542.06,0.15,0.247,0.842,-20.489 +units=m +no_defs -f "%%.2f" <Yorkdeg.txt >Yorkm.txt

Las coordenadas resultantes son:

-52600.52   -42430.79 -0.00
-42950.41   -42439.52 -0.00
-52606.04   -48874.84 -0.00
-42975.81   -48872.07 -0.00

Utilizando esos valores en el georreferenciador QGIS con la misma proyección que el CRS personalizado para el CRS de destino, el resultado se ajusta bien al fondo moderno de Openstreetmap alrededor de Holmfirth:

Como alternativa, puede importar las coordenadas de las esquinas en grados a QGIS como texto delimitado utilizando EPSG:4277 (que es lo más parecido a lo que utilizaban los topógrafos como lat/lon), y establecer la opción proyecto CRS al CRS personalizado mencionado anteriormente. A continuación, utilice from map canvas en el georreferenciador para obtener las coordenadas cassini de los puntos de esquina.