Quel système de coordonnées est ce raster netCDF?

J'ai obtenu un fichier raster netCDF, mais je n'ai pas pu obtenir de métadonnées pour obtenir le nom du système de coordonnées sur lequel le raster a été construit. Le raster lui-même n'a pas de système de coordonnées intégré. Je pensais que ce serait juste un WGS84, et cela ressemblait à cela à première vue, mais avec une enquête plus approfondie dans ArcMap, j'ai vu que c'était un système plutôt pas commun. Voici comment il s'affiche:

Le raster orange est un raster normal dans WGS84 que j'ai inséré ici à des fins de comparaison. Le violet est le raster avec le système de coordonnées inconnu. Avez-vous une idée de ce que cela pourrait être?

Quelques mises à jour: Voici le raster netCDF: https://www.dropbox.com/s/nottbl9yt6dwss6/sic_average_nclimate.nc?dl=0 J'ai également pu obtenir des métadonnées auprès du fournisseur d'images:

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

Ils montrent la frontière lat, mais apparemment aucune information concernant le système de coordonnées.

Vous pouvez déplacer les données dans la position correcte à l'aide de ces commandes 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

(avec un peu d'aide de Frank Warmerdam: Comment reprojeter un raster de 0 360 à -180 180 avec une coupe de 180 méridiens )

et le résultat ressemble à ceci dans une projection stéréographique:

(Ce n'est pas une réponse, mais plutôt une suggestion de processus trop longue pour être un commentaire)

Vous pouvez le géoréférencer dans QGIS en utilisant l'outil GeoReferencer car vous semblez être en mesure d'identifier des points communs sur les deux rasters. Cela générerait un script GDAL avec un ensemble de points GCP et des paramètres de transformation. Vous pouvez ajouter ce script avec un échantillon de ces points à votre question, car il pourrait expliquer la différence entre les deux projections.

Dans le pire des cas, vous auriez une estimation approximative des valeurs basée sur le raster repositionné qui pourrait vous aider à atteindre votre objectif.

Le coin inférieur gauche est aux coordonnées (0,0), mais pour moi (0,0) est la coordonnée du coin supérieur droit dans votre cas et je ne connais aucun système de coordonnées géographiques avec son origine à -90.

Lorsque vous travaillez dans ArcGIS, vous pouvez utiliser la barre d' outils de géoréférencement pour déplacer vos données au bon endroit. Déplacez le coin supérieur droit vers (0,0) et le coin inférieur gauche vers (-360, -90) et cela devrait produire une transformation exacte mais vous auriez toujours un problème avec la ligne de données.

Vous pouvez ainsi essayer ce lien qui devrait vous aider à créer un système de coordonnées géographiques personnalisé. Pour résumer, dans votre système de référence, créez un WGS 84 personnalisé avec le méridien principal à -180, puis une transformation personnalisée en utilisant la rotation en longitude.

Le géoréférencement dans ce nouveau système consisterait à fixer l'origine à (-180, -90). Un point de contrôle suffit car vous avez juste besoin d'une traduction. Sélectionnez soigneusement le pixel inférieur gauche, faites un clic droit et ajoutez les valeurs XY exactes). Ensuite, cela devrait ressembler à ceci

Enfin, il n'y a aucune certitude que le système de coordonnées est basé sur WGS 84. Mais c'est à coup sûr un système de coordonnées géographiques en degrés et la résolution de votre raster est probablement plus grossière que l'erreur potentielle due à une mauvaise sélection de référence.