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

10

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: entrez la description de l'image ici

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.

multigoodverse
la source
1
L'Antarctique est généralement en projection stéréographique polaire, donc tout va du nord au pôle sud. nsidc.org/data/polar_stereo/ps_grids.html
Mapperz
1
Bizarre, ça ressemble à 0-360, mais pourquoi complote-t-il dans l'hémisphère nord?
mkennedy
1
Si les extensions signalées dans l'en-tête NetCDF ressemblent à lat / lon, mais sont lat: 0 à?, Lon: 0 à 360. Pouvez-vous les pirater à lat: -90 à (? - 90) et lon: -180 à + 180?
mkennedy
@Mapperz & mkennedy, j'ai essayé de lui attribuer une projection stéréographique polaire, mais cela n'a pas fonctionné. Je peux probablement le refléter verticalement vers le bas, mais sa moitié droite sera toujours au-delà des longitudes 180E. J'ai ouvert le netCDF avec le bloc-notes, mais je n'ai trouvé aucune information liée à lat, long. J'attache le raster à la question si vous êtes curieux de l'ouvrir de votre côté.
multigoodverse
Ce lien ne fonctionne pas - dropbox indique "le propriétaire n'a pas autorisé l'accès".
BradHards

Réponses:

6

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:

entrez la description de l'image ici

AndreJ
la source
1

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

djq
la source
+ 1 et oui, le géoréférencement pourrait être la dernière puce, mais le processus générerait des erreurs car il est difficile de faire correspondre les points de contrôle. En outre, le raster netCDF devrait être divisé car il n'a pas le même méridien principal. Par exemple, la queue de l'Antarctique est à droite, tandis que le raster orange la montre à gauche.
multigoodverse
0

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.

entrez la description de l'image ici

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.

entrez la description de l'image ici

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

entrez la description de l'image ici

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.

radouxju
la source