Quale sistema di coordinate è questo raster netCDF?

Ho ottenuto un file raster netCDF, ma non sono riuscito a ottenere metadati per ottenere il nome del sistema di coordinate su cui è stato costruito il raster. Il raster stesso non ha alcun sistema di coordinate incorporato. Ho pensato che sarebbe stato solo un WGS84, e a prima vista sembrava così, ma con ulteriori ricerche in ArcMap, ho visto che si trattava di un sistema piuttosto non comune. Ecco come viene visualizzato:

Il raster arancione è un normale raster in WGS84 che ho inserito qui a scopo di confronto. Quello viola è il raster con il sistema di coordinate sconosciuto. Hai idea di cosa potrebbe essere?

Alcuni aggiornamenti: Ecco il netCDF raster: https://www.dropbox.com/s/nottbl9yt6dwss6/sic_average_nclimate.nc?dl=0 Sono stato anche in grado di ottenere alcuni metadati dal provider di immagini:

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" ;

Mostrano il confine lat-long, ma apparentemente nessuna informazione riguardo al sistema di coordinate.

Puoi spostare i dati nella posizione corretta usando questi comandi 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 piccolo aiuto di Frank Warmerdam: come riproiettare il raster da 0 360 a -180 180 con il taglio di 180 meridiani )

e il risultato si presenta così in una proiezione stereografica:

(Questa non è una risposta, ma piuttosto un suggerimento di processo che è troppo lungo per essere un commento)

È possibile geo-referenziarlo in QGIS utilizzando lo strumento GeoReferencer in quanto sembra essere in grado di identificare punti comuni su entrambi i raster. Ciò genererebbe uno script GDAL con un set di punti GCP e parametri di trasformazione. È possibile aggiungere questo script con un campione di questi punti alla domanda, poiché potrebbe spiegare qual è la differenza tra le due proiezioni.

Nel peggiore dei casi, avresti una stima approssimativa dei valori basata sul raster riposizionato che potrebbe aiutarti con il tuo obiettivo.

L'angolo in basso a sinistra è a (0,0) coordinate, ma per me (0,0) è la coordinata dell'angolo in alto a destra nel tuo caso e non conosco alcun sistema di coordinate geografiche con la sua origine a -90.

Mentre lavori in ArcGIS, puoi utilizzare la barra degli strumenti di georeferenziazione per spostare i dati nella posizione corretta. Sposta in alto a destra su (0,0) e in basso a sinistra su (-360, -90) e questo dovrebbe produrre una trasformazione esatta ma avresti comunque un problema con la linea di dati.

Puoi quindi provare questo link che dovrebbe aiutarti a creare un sistema di coordinate geografiche personalizzato. Per riassumere, nel tuo sistema di riferimento, crea un WGS 84 personalizzato con il meridiano primo a -180, quindi una trasformazione personalizzata usando la rotazione della longitudine.

La georeferenziazione in questo nuovo sistema consisterebbe nel impostare l'origine su (-180, -90). Un punto di controllo è sufficiente perché hai solo bisogno di una traduzione. Seleziona con attenzione il pixel in basso a sinistra, fai clic con il pulsante destro del mouse e aggiungi i valori XY esatti). Quindi dovrebbe apparire così

Come osservazione finale, non vi è alcuna certezza che il sistema di coordinate sia basato su WGS 84. Ma è sicuramente un sistema di coordinate geografiche in gradi e la risoluzione del raster è probabilmente più grossolana del potenziale errore a causa di una scelta errata del riferimento.