Message posté par : Geraud DBDA (geraud.debouchard(a)gmail.com)
----------------------------------------
Bonjour à tous,
Je me retrouve face à une difficulté sans avoir trouver de réponse, peut-être que vous
pouvez m'aider.
Tout d'abord, j'ai 2 jeux de données, sur une grande partie dans les Alpes :
- Une couche vecteur représentant un maillage (fishnet) où chaque maille mesure 25x25 m.
Ce maillage est calé au km Lambert pour le besoin de la base de données.
- Un MNT de 5 m de résolution qui, une fois dérivé, servira à alimenter les données de ces
fameuses mailles.
Ce MNT n'est (évidemment) pas calé au km Lambert et j'ai un décalage de +2,5 m en
X et Y entre les pixels de 5m du raster et mon maillage vecteur.
Exemple :
Coin inférieur gauche d'une maille vecteur : 962875.00 , 6427200,00
Coin inférieur gauche d'un pixel du MNT : 962872.5 , 6427197.5
Mon problème :
Je voudrais ajouter qu'une seule information provenant du MNT dans mes mailles.
Exemple :
Maille n : Altitude = 2000
Dans mon cas, je me retrouve donc avec plusieurs pixels au sein de ma maille 25m (au total
25 pixels, ceux sur les cotés sont coupés en deux, ceux dans les coins sont coupés en
quatre par rapport aux mailles) et, j'aimerai avoir une sorte de centroïde qui, au
centre de chaque maille, renseignerai une valeur lissée du MNT. Une maille, une valeur du
MNT.
J'avais imaginé le cheminement suivant :
1. Décaler tous les pixels du MNT de 2,5m en X,Y. Mais de ne pas faire une simple
transposition, sinon les nouvelles valeurs n'auront plus aucuns sens. Donc en même
temps de recalculer chaque pixel pour que leur valeurs correspondent bien à leur nouvelle
position.
2. Redimensionner le MNT de 5m en 25m pour avoir une belle superposition entre les mailles
vecteurs et les pixels du MNT. Ce qui suppose de faire une sorte de moyenne des valeurs
des 25 pixels de 5m.
3. Extraire les données du MNT dans mes mailles (ça je sais faire).
Déjà pensez-vous que ce serait la bonne façon de faire ? N'y-aurait-il pas de
raccourci ou de formule magique ?
Du coup je ne sais vraiment pas comment m'y prendre pour les points 1 et 2...
Logiciel : QGis, ArcGis, même SagaGis peu importe.
Une illustration d'une maille de 25 m et du MNT derrière (ici Aspect) :
https://zupimages.net/up/21/22/p28t.jpg
(
https://zupimages.net/viewer.php?id=21/22/p28t.jpg)
Merci d'avance de vos réponses.
----------------------------------------
Le message est situé
https://georezo.net/forum/viewtopic.php?pid=343817#p343817
Pour y répondre : geomatique(a)ml.georezo.net ou reply de votre messagerie
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net