Le portail francophone de la géomatique

Bonjour,

j'aimerai pouvoir utiliser l'indice IDPR (Indice de Développement et de Persistance des Réseaux) du BRGM pour mettre en évidence mes zones d'infiltrations.
Quelqu'un a-t-il déjà travaillé avec cet indice et connaît-il la méthodologie SIG ?

"L'IDPR permet de quantifier la perméabilité du bassin en comparant un réseau théorique établi selon l'hypothèse d'un milieu parfaitement homogène, au réseau naturel mis en place sous le contrôle d'un contexte géologique hétérogène. Il permet ainsi de rendre compte si le sol laisse infiltrer ou ruisseler les eaux de surface[...] L’IDPR résulte du rapport, pour tout point de l'espace étudié, entre la distance au plus proche cours d'eau calculé et la distance au plus proche cours d'eau réel. Par soucis de simplification du critère, L’IDPR est arbitrairement multiplié par 1000."

j'ai bien essayer d'établir une suite de traitement mais le résultat final ne donne pas ce que je veux:
- Etablissement du MNA à partir de mes données altimétriques (ArcToolbox > 3D Analyst Tools > Raster Interpolation > Topo to raster)
- Création de la carte de direction de flux, qui est nécessaire pour la suite du traitement (ArcToolbox > Spatial Analyt Tools > Hydrology > Flow Direction)
- Création de la carte d'accumulation de flux, qui me permet d'obtenir mon réseau hydrographique théorique (ArcToolbox > Spatial Analyt Tools > Hydrology > Flow accumulation)
- Vectorisation du réseau hydrographique théorique :
           Application d’un seuil dans le raster (Boîte à outils Spatial Analyst > Raster Calculator) : streamnet50 = con (flowacc > 50, 1)
           puis Ecoulement vers entité (ArcToolbox > Spatial Analyt Tools > Hydrology > Stream to Feature)
- Lissage des lignes : ArcToolbox > Data Management Tools > Generalization > Smooth line
- Calcul de la distance, pour tout point de l’espace, au plus proche cours d’eau théorique, sur la base du réseau hydro théorique obtenu précédemment (ArcToolBox > Spatial Analyst Tools > Distance > Euclidean Distance) --> [Carte distance euclidienne réseau hydro théorique]
- Calcul de la distance, pour tout point de l’espace, au plus proche cours d’eau réel, sur la base de mon réseau hydro réel (ArcToolBox > Spatial Analyst Tools > Distance > Euclidean Distance) --> [Carte distance euclidienne réseau hydro réel]
- Calcul de l’IDPR (Boîte à outils Spatial Analyst > Raster Calculator) : et c'est là que je dois avoir un soucis puisque mes résultats ne sont pas cohérents (notamment au niveau du réseau où les valeurs des pixels sont égales à zéro, mais aussi ailleurs où mes zones infiltrantes et mes zones de ruissellement sont inversées). J'utilise la formule suivante :
         IDPR = ([Carte distance euclidienne réseau hydro théorique] / ([Carte distance euclidienne réseau hydro réel]) * 1000


Si quelqu'un a déjà utilisé cet indice, merci de m'indiquer où je me trompe dans le traitement.

Merci d'avance pour vos réponses