Suite

Le raster en sortie de l'outil d'interpolation IDW va de 3.40282e+038 à -3.40282e+038 dans ArcPy


Je ne connais pas Python dans ArcGIS, j'ai donc d'abord utilisé ModelBuilder pour définir l'environnement et y mettre l'outil. Le modèle fonctionnait correctement et le raster en sortie semblait bon.

Ensuite, j'ai exporté le modèle vers un script Python. La sortie raster est générée, mais la valeur de raster est comprise entre 3.40282e+038 et -3.40282e+038. Pourquoi cela s'est-il produit et comment puis-je y remédier?

Mon script est comme ceci :

# Importer le module arcpy import arcpy # Vérifier toutes les licences nécessaires arcpy.CheckOutExtension("3D") # Définir les environnements de géotraitement arcpy.env.scratchWorkspace = "C:DocumentsArcGISiomass.gdb" arcpy.env.outputCoordinateSystem = "PROJCS['WGS_1984_Plate_Carree',GEOGCS['GCS_WGS_1984',DATUM['D_WGS_1984',SPHEROID['WGS_1984',6378137.0,298.257223563]],PRIMEM['Greenwich',0.0],UNIT,['5329291994' ,PROJECTION['Plate_Carree'],PARAMETER['False_Easting',0.0],PARAMETER['False_Northing',0.0],PARAMETER['Central_Meridian',0.0],UNIT['Meter',1.0]]" arcpy.env.extent = "60.8785969898536 -11.0048500049116 148.898540012983 53.560734" arcpy.env.cellSize = "MAXOF" arcpy.env.workspace = "C:Desktopshape" arcpy.env.mask = "newasia.shp" # Variables locales : xyz_shp = "xyz.shp" b1_ald22222 = "C:DocumentsArcGISmass.gdb1_ald22222" # Processus : IDW arcpy.Idw_3d(xyz_shp, "ALDX", b1_ald22222, "27000", "2", " VARIABLE 8", "")

La raison est probablement dans ArcMap incapable de reconnaître une valeur NoData définie par l'outil IDW comme 3.40282e+038 ou -3.40282e+038.

je recommanderais d'utiliserSetNulloutil pour ces deux valeurs sur le résultat et voyez si cela élimine toutes les valeurs insensées (devrait les définir correctement sur NoData). Sinon, une faille dans les données d'entrée peut être présente (je ne sais pas si IDW gère correctement les valeurs NULL, par exemple).


J'ai résolu ce problème et il y a une petite erreur ici. La portée de l'interpolation est de 27 000 mètres mais l'étendue de l'environnement arcpy est en degré. La seule chose à faire est donc d'unifier les unités.


Voir la vidéo: Arcgis: Merge the raster datasets fusionner raster arcgis Mosaic raster dataset ArcGIS (Octobre 2021).