Suite

Erreur de projection, géotransformation non valide

Erreur de projection, géotransformation non valide


J'essaie d'extraire des données d'ensembles de données météorologiques maillées tels que GLDAS à partir de pixels dans lesquels se trouvent les stations météorologiques. J'ai créé un fichier de formes ArcGIS de points de station météo. Ces coordonnées de point de station météo étaient dans GCS_North_American_1983.

Mes données raster GLDAS ont un système de coordonnées géographiques étrange (voir ci-dessous). Pour extraire correctement les valeurs de pixels de mon raster sous les points de mon fichier de formes, je dois reprojeter le fichier de formes pour qu'il soit dans le même GCS que les rasters GLDAS. Lorsque j'essaie de projeter les points de la station météo vers le GLDAS GCS, j'obtiens l'erreur suivante : "000151 : La ou les géotransformations ne sont pas valides."

J'utilise l'outil Projet dans ArcMap. Je sélectionne le système de coordonnées en entrée pour correspondre aux données à reprojeter (GCS_North_American_1983) et je sélectionne le système de coordonnées en sortie pour correspondre à mon raster GLDAS. La référence spatiale des données raster est la suivante :

Autorité GCS_Coordinate_System_imported_from_GRIB_file : Personnalisé

Unité angulaire : Degré (0,0174532925199433) Premier méridien : Greenwich (0,0) Datum : D_unknown Sphéroïde : Sphère Axe semi-majeur : 6367470,0 Axe semi-mineur : 6367470,0 Aplatissement inverse : 0,0

Je n'ai pas créé ce GCS, c'est ce que le raster est venu avec. La projection vers le GCS ci-dessus échoue, mais j'ai pu projeter les mêmes points de station météo vers ce système de coordonnées de sortie subtilement différent à partir d'un ensemble de données météo quadrillé similaire appelé NLDAS :

Autorité GCS_Coordinate_System_imported_from_GRIB_file : Personnalisé

Unité angulaire : Degré (0,0174532925199433) Méridien premier : Greenwich (0,0) Datum : D_unknown Sphéroïde : Sphère Axe semi-majeur : 6371200,0 Axe semi-mineur : 6371200,0 Aplatissement inverse : 0,0

Je n'arrive pas à comprendre pourquoi j'obtiens une erreur pour la première projection GCS et pas la seconde.


Je ne sais pas pourquoi l'un fonctionne mais pas l'autre. Ni l'un ni l'autre ne devrait travailler.

Vous avez deux choix. L'une consiste à utiliser l'outil Créer une transformation géographique personnalisée. Configurez-le comme l'outil Projet et importez le système de référence de coordonnées en sortie (CRS) à partir des données raster. Étant donné que vous obtenez le même nom CRS géographique à partir des deux jeux de données, faites en sorte que le nom de la transformation reflète la sphère ou le jeu de données avec lequel il doit fonctionner.

Essayez d'utiliser la méthode Null. S'il existe un moyen de le vérifier et que vous pensez que les données ne se superposent pas correctement, réessayez en utilisant la méthode de traduction géocentrique. Dans le premier cas, les valeurs de longitude et de latitude sont inchangées. Dans le second cas, les valeurs de latitude sont converties entre NAD83 et la sphère.

Une fois cela fait, vous devriez pouvoir l'utiliser dans l'outil Projet ou dans ArcMap lui-même.

Vous pouvez également ajouter d'abord la couche raster à un nouvel ArcMap afin qu'il définisse le système de coordonnées du bloc de données, puis ajouter la couche de station. Si vous ne souhaitez pas effectuer de transformation (AKA Null fonctionne), cliquez avec le bouton droit sur la couche de la station et sélectionnez les données, exportez les données. Dans cette boîte de dialogue, choisissez d'utiliser le système de coordonnées du bloc de données.