¿Qué sistema de coordenadas es este ráster netCDF?


10

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: ingrese la descripción de la imagen aquí

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.


1
La Antártida suele estar en proyección estereográfica polar, por lo que todo va hacia el norte desde el polo sur. nsidc.org/data/polar_stereo/ps_grids.html
Mapperz

1
Extraño, parece 0-360, pero ¿por qué está tramando en el hemisferio norte?
mkennedy

1
Si las extensiones informadas en el encabezado NetCDF se ven como lat / lon, pero son lat: 0 a?, Lon: 0 a 360. ¿Puede hackearlas a lat: -90 a (? - 90) y lon: -180 a + 180?
mkennedy

@Mapperz & mkennedy, intenté asignarle una proyección estereográfica polar, pero no funcionó. Probablemente pueda reflejarlo verticalmente hacia abajo, pero su mitad derecha seguirá estando más allá de las longitudes de 180E. Abrí el netCDF con el bloc de notas, pero no pude encontrar ninguna información relacionada con lat, por mucho tiempo. Adjunto el ráster en la pregunta si tiene curiosidad por abrirlo de su lado.
multigoodverse

Ese enlace no funciona: Dropbox dice "el propietario no ha autorizado el acceso".
BradHards

Respuestas:


6

Puede cambiar los datos a la posición correcta utilizando estos comandos GDAL:

gdal_translate -a_srs EPSG:4326 -a_ullr 0 0 360 -90 NETCDF:"sic_average_nclimate.nc":SIC_Change change360.tif
gdalwarp -t_srs WGS84 change360.tif change180.tif  -wo SOURCE_EXTRA=1000 --config CENTER_LONG 0

(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:

ingrese la descripción de la imagen aquí


1

(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.


+ 1 y sí, la georreferenciación podría ser la última viñeta, pero el proceso generaría algunos errores porque es difícil hacer coincidir los puntos de control. Además, el ráster netCDF tendría que dividirse porque no tiene el mismo meridiano principal. Por ejemplo, la cola de la Antártida está a la derecha, mientras que la trama naranja la muestra a la izquierda.
multigoodverse

0

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.

ingrese la descripción de la imagen aquí

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.

ingrese la descripción de la imagen aquí

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í

ingrese la descripción de la imagen aquí

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.

Al usar nuestro sitio, usted reconoce que ha leído y comprende nuestra Política de Cookies y Política de Privacidad.
Licensed under cc by-sa 3.0 with attribution required.