Definir Proj4 cadena para CANGRD conjunto de datos

Question

Estoy tratando de georreferencia de la .grd archivos proporcionados aquí http://open.canada.ca/data/en/dataset/3d4b68a5-13bc-48bb-ad10-801128aa6604.

La documentación indica que "El CANGRD red de la proyección estereográfica polar con un 50 km de resolución espacial. La red es un 125 (columnas) por 95 (filas) de la matriz, donde la esquina suroeste (0,0) está en 40.0451°de latitud N y 129.8530°W de longitud. La proyección es cierto a 60.0°N y centrada en 110.0°W."

El encabezado de cada uno .grd archivo contiene esto:

DSAA
125 95
0.0 124
0.0 94
-3.051277 3.268792 

He intentado con el siguiente Proj4 cadena, pero todo el conjunto de datos aparece inclinada.

+proj=stere +lat_0=40.0451 +lat_ts=60 +lon_0=-129.8530 +x_0=0 +y_0=0 +datum=WGS84 +k=1 +to_meter=50000

Alguna idea de cómo correctamente la georreferencia de estos archivos?

Answer 1

Así que hice lo que AndreJ sugerido y poner lat_0=90 para generar este sistema temporal

+proj=stere +lat_0=90 +lat_ts=60 +lon_0=-110 +x_0=0 +y_0=0 +datum=WGS84 +to_meter=50000

Entonces me transformó la latitud, la longitud de las esquinas a este nuevo sistema de referencia y utiliza esas coordenadas para calcular el x_0 y y_0.

+proj=stere +lat_0=90 +lat_ts=60 +lon_0=-110 +x_0=1884770 +y_0=5220000 +datum=WGS84 +to_meter=50000

Esto parece dar una definición perfecta para el CANGRID archivos. Tenga en cuenta que en este sistema de referencia, las unidades son "las celdas de la cuadrícula", o metros/50000

Creo que el equivalente WKT .prj contenido del archivo es:

PROJCS["CANGRD_Polar_Stereographic",GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Stereographic_North_Pole"],PARAMETER["false_easting",37.6954],PARAMETER["false_northing",104.4],PARAMETER["central_meridian",-110.0],PARAMETER["standard_parallel_1",60.0],UNIT["50_Kilometers",50000.0]]