Suite

Preuve de faibles écarts dans le DEM par rapport à l'ensemble de données de points


Je ne sais pas comment obtenir les résultats suivants de manière cohérente :

J'ai un jeu de données de points distribués réguliers (élévations de surface). J'en ai tiré plusieursDEMutilisantInterpolation QGIS IDWméthode. Maintenant, j'obtiens des résultats différents en fonction de la valeur du coefficient de distance que j'entre. J'ai examiné les valeurs entre1.0et3.5, où des valeurs plus faibles lissent davantage la surface et des valeurs plus élevées augmentent la rugosité de la surface et ainsi les caractéristiques du terrain apparaissent plus clairement.

Dans l'étape suivante, j'ai comparé les valeurs individuelles de monDEMrasters avec la valeur à la même position dans le jeu de données de points d'origine (que je suppose être les valeurs "correctes" du monde réel, je veux également être représenté dans le DEM). J'ai trouvé des écarts plus importants dans les rasters avec un coefficient de distance plus faible, en raison du lissage de la surface.

Maintenant, je veux m'assurer quel coefficient de distance convient le mieux à mon calcul, essayant ainsi de calculer l'écart moyen entre un raster et le jeu de données de points sur toute la zone. Je pensais que je pouvais juste calculerraster - jeu de données de points, puis calcul de la valeur moyenne sur le jeu de données résultant, mais cela ne fonctionnera pas avec la calculatrice raster.

Je ne suis malheureusement pas non plus un expert en QGIS.

Existe-t-il une fonction réalisant ce qui précède ou une solution de contournement ?

En cherchant sur Google, je n'ai rien trouvé qui puisse faire ça pour moi.


Sur la base des informations ajoutées dans votre commentaire, je suggère d'utiliser le module GRASS r.in.xyz. Il y a eu quelques problèmes récemment avec le plugin QGIS GRASS, donc je travaillerais directement dans GRASS, de manière autonome. Tout d'abord, lorsque vous démarrez GRASS, vous devez configurer la base de données GRASS et LOCATION. L'EMPLACEMENT est défini par les paramètres de projection des coordonnées. Assurez-vous qu'il est le même que le système de coordonnées de vos points de données d'entrée.

Une fois que vous avez fait cela, utilisez le paramètre -s dur.in/xyzmodule pour rechercher l'étendue complète de votre fichier de données de points d'entrée. Utilisez ces étendues pour définir la région (module d'herbeg.région). De plus, réglez la résolution sur 20 m, en utilisant à nouveau le mêmeg.régionmodule.

Maintenant tu peux juste courirr.in.xyz, et un raster GRASS sera créé avec la valeur de chaque pixel égale à la valeur du point qu'il contient. Vous avez donc votre DEM. Et puis vous pouvez exporter vers un GTiff, ou un autre format raster pour d'autres utilisations avecr.out.gdal.