Obtuve un archivo ráster netCDF, pero no pude obtener ningún metadato para obtener el nombre del sistema de coordenadas en el que se ha construido el ráster. El ráster en sí no tiene ningún sistema de coordenadas incrustado. Pensé que sería solo un WGS84, y parecía a primera vista, pero con más investigación en ArcMap, vi que es un sistema bastante poco común. Así es como se muestra:
La trama naranja es una trama normal en WGS84 que he insertado aquí para fines de comparación. El morado es el ráster con el sistema de coordenadas desconocido. ¿Tienes alguna idea de lo que podría ser?
Algunas actualizaciones: Aquí está el ráster netCDF: https://www.dropbox.com/s/nottbl9yt6dwss6/sic_average_nclimate.nc?dl=0 También pude obtener algunos metadatos del proveedor de imágenes:
netcdf sic_average_nclimate {
dimensions:
nlon = 361 ;
nlat = 90 ;
nseas = 4 ;
variables:
float SIC_Change(nlat, nlon) ;
SIC_Change:Title = "Gridded Multi-Model Ensemble Mean Annual Mean Change in Ice Concentration 21C-20C" ;
float SIC_Season_Change(nseas, nlat, nlon) ;
SIC_Season_Change:Title = "Gridded Multi-Model Ensemble Mean Seasonal Mean Change in Ice Concentration 21C-20C" ;
float SIC_Change_STD(nlat, nlon) ;
SIC_Change_STD:Title = "Gridded Multi-Model Standard Deviation of the Annual Mean Change in Ice Concentration 21C-20C" ;
float SIC_Season_Change_STD(nseas, nlat, nlon) ;
SIC_Season_Change_STD:Title = "Gridded Multi-Model Standard Deviation of the Seasonal Mean Change in Ice Concentration 21C-20C" ;
float LAT(nlat) ;
LAT:Title = "Latitude" ;
float LON(nlon) ;
LON:Title = "Longitude" ;
// global attributes:
:Title = "Ice Concentration metrics for Model subset as in Figure 1 of NCLIMATE paper" ;
Muestran el límite de largo, pero aparentemente no hay ninguna información sobre el sistema de coordenadas.
fuente
Respuestas:
Puede cambiar los datos a la posición correcta utilizando estos comandos GDAL:
(con un poco de ayuda de Frank Warmerdam: Cómo reproyectar el ráster de 0 360 a -180 180 con un meridiano de corte 180 )
y el resultado se ve así en una proyección estereográfica:
fuente
(Esto no es una respuesta, sino una sugerencia de proceso que es demasiado larga para ser un comentario)
Puede geo-referenciarlo en QGIS utilizando la herramienta GeoReferencer, ya que parece ser capaz de identificar puntos comunes en ambos rásteres. Esto generaría un script GDAL con un conjunto de puntos GCP y parámetros de transformación. Puede agregar este script con una muestra de estos puntos a su pregunta, ya que podría explicar cuál es la diferencia entre las dos proyecciones.
En el peor de los casos, tendría una estimación aproximada de los valores en función del ráster reposicionado que podría ayudarlo con su objetivo.
fuente
La esquina inferior izquierda está en coordenadas (0,0), pero para mí (0,0) es la coordenada de la esquina superior derecha en su caso y no conozco ningún sistema de coordenadas geográficas con su origen en -90.
Mientras trabaja en ArcGIS, puede usar la barra de herramientas de georreferenciación para mover sus datos a la ubicación correcta. Mueva la esquina superior derecha a (0,0) y la esquina inferior izquierda a (-360, -90) y esto debería producir una transformación exacta, pero aún tendría un problema con la línea de fecha.
Por lo tanto, puede probar este enlace que debería ayudarlo a crear un sistema de coordenadas geográficas personalizado. Para resumir, en su sistema de referencia, cree un WGS 84 personalizado con el meridiano principal en -180, luego una transformación personalizada utilizando la rotación de longitud.
La georreferenciación en este nuevo sistema consistiría en establecer el origen en (-180, -90). Un punto de control es suficiente porque solo necesita una traducción. Seleccione cuidadosamente el píxel inferior izquierdo, haga clic derecho y agregue los valores XY exactos). Entonces debería verse así
Como observación final, no hay certeza de que el sistema de coordenadas se base en WGS 84. Pero es seguro un sistema de coordenadas geográficas en grados y la resolución de su ráster es probablemente más gruesa que el error potencial debido a una selección de datos incorrecta.
fuente