Suite

Les données d'altitude de « nivellement » coïncidant avec le polygone ?

Les données d'altitude de « nivellement » coïncidant avec le polygone ?


J'ai un fichier de formes contenant des géométries de construction (polygones) et un jeu de données d'altitude. Les données d'élévation ne correspondent pas parfaitement aux bâtiments.

Je veux "aplatir" les données d'altitude de sorte que partout où il y a un bâtiment, l'altitude soit complètement plate.

En bref : pour tous les pixels qui coupent un polygone, modifiez leur valeur à la moyenne des pixels qui se coupent.

Comment puis-je accomplir cela?

Cette question n'est plus ou moins liée à un logiciel spécifique, bien que je préfère utiliser des logiciels gratuits comme QGIS, GDAL/OGR, etc.


Pour les valeurs moyennes d'altitude dans les polygones de construction…

Dans QGIS :

  1. Activer Statistiques zonales brancher
    • Exécutez-le - "Statistiques raster/zonales", il ajoutera de nouvelles colonnes à la couche de polygones
    • celui qu'il faut est "signifier".
  2. Dupliquez votre couche d'altitude (raster/grille) - vous l'utiliserez pour le fichier final
  3. Utilisation "Raster / Conversion / Rasteriser"
    • sélectionnez votre couche de polygone (bâtiment), champ "moyenne"
    • comme cible, utilisez votre copie de raster/grille d'altitude.
    • Conserver la taille et la résolution existantes
    • (pour une raison quelconque, ma version QGIS 2.2 écrit "Le fichier sélectionné n'est pas un format OGR pris en charge" pour tous les fichiers vectoriels que j'ai essayés, j'ai donc dû utiliser QGIS dev version 2.3)

Le résultat sera une combinaison des valeurs d'altitude d'origine à l'extérieur des polygones et de la moyenne à l'intérieur des polygones.


Cette réponse ne fonctionne que si vous voulez une valeur constante - dans l'exemple 0.

Dans QGIS :

  1. Découper l'élévation par polygones (bâtiments) vers un nouveau raster (grille) - buildings_grd

    • Raster/Extraction/Clipper
    • (Utilisez le mode Découpage "Couche de masque" et votre fichier de formes de bâtiment)
  2. Exécutez la calculatrice raster et calculez

    • Raster/Calculatrice raster…
    • Expression : ("élévation" = "buildings_grd")*0 + ("élévation" != "buildings_grd")*"élévation"

Voir la vidéo: nivellement par cheminement encadrée