Nous avons un DEM LIDAR de 1 mètre d'une ville.
Un petit sous-ensemble peut être téléchargé à partir de ce lien :
Cette capture d'écran montre un DEM brut avec une palette de gris (les ceintures plus sombres sont des rues et les rectangles grisâtres et blanchâtres sont des blocs):
Cela correspond à un endroit dans la ville de Santo Domingo, qui peut être vu dans cette capture d'écran de Google:
En moyenne, les blocs sont "élevés" ca. 2 mètres des rues , ce qui n'est pas vrai. Nous voulons avoir un MNT propre pour générer un réseau de cours d'eau et un indice d'humidité topographique (TWI). Avec le DEM fourni (nous n'avons pas les bandes originales du scanner laser), le réseau hydrgraphique semblait suivre une disposition rectangulaire, et TWI a abouti à un motif de blocs . Ces images montrent les résultats:
Ceci est le résultat de flux réseau, généré avec r.watershed
en Grass GIS
:
Et voici le résultat TWI, généré avec SAGA
:
Nous avons essayé quelques procédures pour résoudre cette inexactitude sans succès:
1) Outil de débruitage . Nous avons appliqué l' r.denoise
outil Grass GIS
, mais nous avons eu quelques problèmes avec l'installation du module. Nous l'avons à nouveau exécuté avec un shell dans Windows et nous avons reçu un message de mémoire insuffisante.
2) Filtres . Nous avons effectué différents types de filtres (filtre passe-bas, la médiane, moyenne, etc.), avec des fenêtres différentes tailles, et d' essayer de mettre du poids dans la direction des rues ( Grass GIS
, SAGA
, QGIS
).
3) Géostatistique . Nous avons généré des nuages de points strictement sur les rues (essayé 1000 et 2000 points), généré un modèle de variogramme et ensuite exécuté un krigeage ordinaire pour remplir les blocs. La modélisation des variogrammes et le krigeage ordinaire ont été effectués en R
utilisant différents packages. Nous avons obtenu un variogramme linéaire, donc nous ne comptions pas sur les résultats du krigeage.
4) Autres outils . ALDPAT
Outil installé , mais impossible de le faire fonctionner car le programme n'a pas pu lire le DEM.
Dans tous les cas, les résultats en termes de réseau de drainage n'étaient pas bons , car nous ne pouvions pas éviter le réseau de cours d'eau rectangulaire; aussi, TWI a toujours abouti à un modèle de blocs .
En particulier, avec un résultat interpolé OK, nous avons obtenu un modèle DEM de type point qui a affecté le résultat du réseau. Cependant, l'effet du motif des blocs a été diminué.
De plus, nous avons jeté un œil à cette question et à ses réponses ...
Filtrer les auvents et les bâtiments du DSM pour avoir une élévation de terre nue
... qui nous a redirigés vers Whitebox Geospatial Analysis Tools
, mais n'a pas pu convertir notre DEM en LAS format
. De plus, nous n'étions pas sûrs de l'efficacité de Bare-Earth DEM tool
pour nous, car il est conçu pour supprimer des objets semi-transparents, et non des blocs incorrectement "soulevés", ce qui est notre cas.
Nous voulons toujours générer un DEM de haute qualité pour effectuer notre analyse hydrographique, mais nous ne savons pas quoi d'autre pouvons-nous essayer.