Suite

Existe-t-il un moyen de modifier les valeurs d'élévation des pixels d'un MNT en fonction de l'élévation du point de BM et du point de contrôle à partir du relevé de terrain réel ?


Nous avons un DEM rapide d'oiseau de notre site. Après avoir obtenu les valeurs d'altitude réelles au sol avec un appareil GPS différentiel, nous avons constaté qu'il existe une différence moyenne de 12 pieds entre le DEM et les points au sol réels. nw, je veux modifier les élévations de ce DEM en fonction des valeurs réelles.


Deux manières de procéder à l'aide d'ArcGIS :

Interpoler - a plus de sens : utilisez l'outil Raster to Point sur votre DEM et définissez l'option Field sur Value. Vous obtiendrez une couche de points pour chaque centre de pixel, avec la valeur d'élévation du pixel stockée pour chaque point. Utilisez ensuite l'outil Créer un TIN et utilisez les deux couches comme entrées - la couche des valeurs d'altitude du GPS et la nouvelle couche de l'outil Raster à Point. après cela, utilisez TIN to Raster. Vous pouvez à la place copier les points d'altitude GPS dans la nouvelle couche de points DEM et effectuer une interpolation telle que IDW ou Natural Neighbor. Ceci ou le TIN peut également être effectué après la suppression des points DEM adjacents aux points GPS.

Remplacer des pixels spécifiques - un peu délicat : utilisez l'outil Pointer vers raster sur vos mesures d'altitude GPS, définissez la taille de cellule de la même manière que la taille de pixel de votre DEM et dans l'outil Environnements, définissez l'étendue du traitement comme votre DEM d'origine. Vous obtiendrez un raster avec de nombreuses valeurs NoData et quelques valeurs d'altitude, en fonction du nombre d'entre elles. Ouvrez maintenant l'outil Field Calculator, utilisez une expression qui remplace les pixels de votre DEM par des pixels du raster GPS qui ne sont pas des NoData :

Con(IsNull("GPS_Raster") == Faux, "GPS_Raster","DEM")

J'ajouterai que j'utilise l'étendue du DEM d'origine et que j'ignore ensuite les NoData car sinon, fonctionnant sur deux rasters, il semble que le calculateur de champ utilise leur intersection, et non l'étendue du plus grand.

Et il serait peut-être plus logique d'utiliser la calculatrice raster pour ajouter ou soustraire une valeur constante de déviation (vous avez mentionné 12 pieds) de votre raster.


Voir la vidéo: Resampling Raster ArcGis. changing the cell size of Raster dataset in ArcGis (Octobre 2021).