Le portail francophone de la géomatique

La largeur du fleuve, c'est à dire, la distance qu'il faut parcourir à la nage pour atteindre l'autre rive ? À vue de nez, c'est très compliqué.

Le dessin de bassin versant dépend déjà d'un certain nombre de paramètres, mais la largeur du fleuve, il faudrait à vue de nez en plus de ce que tu as énoncé :

Un MNT assez fin du lit du fleuve.
La période de l'année.
Les précipitations.
L'occupation du sol du BV.
Le type de sol du BV, sans parler des couches géologiques.
Les prélèvements artificiels...
Pour la proximité de l'embouchure, il faut aussi les coefficients de marrée !

Et alors il faudrait un modèle en béton pour avoir en sortie ce que tu demandes. Un joli sujet d'hydrologie. Mais je crois qu'il va falloir prendre un bon nombre d'hypothèse et donc s'y connaitre vraiment en hydrologie pour faire tout ça.


Je pense que le plus simple est la télédétection et non pas le calcul.

Bonjour à vous,

Merci Jérome C et Robin pour les conseils, désolé je  ne trouve pas d'outils en tapant hydrologie ? Ou dois je chercher svp ?

je fais aussi un peu de Télédétection, je possède des images JERS Basses Eaux et Hautes eaux de ma zone d'étude que je compte exploiter pour déterminer une zone en eaux BE et HE. Ces images sont de plusieurs dates mais peu importe cela donnera largeur min et largeur max.
Pour ce qui est des autres données je n'ai rien à part les données Hydrosheds dont le SRTM3 (90m) sur la zone.
Le fleuve c'est l'Amazone de Manaus à Obidos pour les connaisseurs (700 000 km²).... et c'est du gros, donc en premier temps du travail un peu grossier sans complexes.
Le but est d'obtenir un profil avec l'altitude en fonction de la distance amont-aval, et sur ce même graphique une 'largeur de fleuve' à intervalles réguliers en fonction de la distance.

Pour la première partie c'est faisable avec l'outil profil ou un équivalent. Maintenant la difficulté est d'avoir une seule ligne matérialisant le cours d'eau.
Dans un même temps estimer la largeur en fonction de la distance parcouru de manière régulière ...c'est plutôt compliqué !
Le but est d'intégrer ces 'mailles' dans un modèle hydrologique et la découverte d'une méthode nous serait très utile !

Jérome qu'entends tu par : prélèvements artificiels ?

Mes idées de calculs :
1> pour extraire le cours d'eau matérialisé par une ligne au milieu des 2 berges. Passer par une classification (grossière) de 2 classes Eau/pas Eau. Transformer résultat en vecteur 'polygone' ne garder que le long du fleuve, nettoyage/union des objets si besoin. Convertir en ligne ou scinder en 2 ??= cours d'eau (ligne) si cela est possible ?

2> Convertir ligne en points (un points tout les 1km ou 5Km au début) ?...faire des tampons autour de mes points, intersections des tampons et des polygones, calcul des surfaces / rayons pour chaque objet ?? = une largeur tout les '5km' de distance.....

.....c'est toujours pareil on se perd avec tout ces outils et fonctions.....

J'entends bien Jérome ce sont beaucoup d'hypothèses et d'approximations mais ce n'est qu'un début de travail.
Merci de vos retours.

BBR

Ok je cherchais vraiment ailleurs et je l'avais sous les yeux.....
Merci Robin

bbr

Hello,

J'ai bien une idée mais il me faut comme ingrédients :
- une couche de polygone qui représente le fleuve en plein (donc toute la largeur que l'on veut mesurer)
- une polyligne qui passe au milieu de ce fleuve
- l'intervalle entre les mesures (tout les combien de mètres ?)

Est ce que tu peux avoir ça ?
Robin.

Bonjour à tous,

j'ai ma couche polygone matérialisant le cours d'eau (1 seul objet) = cours d'eau principal + les plus gros affluents.
Cette couche est issue d'une classification, puis convertie en polygones, puis simplifiée/agrégée afin d'obtenir un polygone unique.

Pour obtenir l'image de la classification j'ai utilisée ArcHydro (nombreux outils équivalents à Spatial Analyst/Hydrologie) j'ai ainsi pu déterminé ma couche cours d'eau à partir de données SRTM.

Mes problèmes :
1> ce polygone n'est pas "plein" il a des trous, comment puis procéder si je veux obtenir seulement les contours ?
2> je ne sais pas comment obtenir une ligne au milieu de ce polygone, savez vous comment je pourrais procéder pour convertir ce polygone en ligne ? Ou le scinder en 2 parties égales ?

Merci de votre aide.

BBR

Bonjour suite à ces questions/réponses,

Pour mon polygone "complexe" avec de nombreux trous je n'ai pas procédé à la suppression des "donuts". J'ai convertit mon polygone en lignes puis j'ai sélectionné à la main les lignes de contour pour créer un nouveau polygone sans trous appelé masque.

Pour le calcul de la ligne centrale il existe un outil appelé "collapse dual lines to centerline" mais pour l'instant je n'arrive pas à obtenir quelque chose de correct.

A+
bbr

Normal, ce n'est que pour ArcInfo

Determining proximity to points: Near, Point Distance, Point Node(**ArcInfo only**)

Bebert,

Si tu te sens prêt à programmer un petit peu, tu peux toujours utiliser les ArcObjects et l'outils qui fait la même chose que l'outils de la toolBox, mais en plus élémentaire : http://edndoc.esri.com/arcobjects/9.2/C … stance.htm je ne vois pas de restriction de license. Le seul problème, c'est la barrière de la programmation.

Sinon, ce problème de topologie, il est grâve à quel point. Ces dépassement de limite, elles sont de quelle longueur à peu près ? Est-ce que ces dépassements sont beaucoup plus grands que la valeur de tolérance que tu peux entrer en option dans la commande de découpage ?

Bebert, pour diagnostiquer, pas pour envisager une méthode de traitement par petits lots, est-ce que tu as essayé de faire les traitements juste sur la zone que tu nous montres. ArcGIS a toujours été un peu frileux avec les très gros lots de données. Il faudrait valider pour ton cas l'hypothèse de frilosité fasse la zone d'étude géante.

Comprend pas, pourquoi ne pas mettre cette info dans la table attributaire des segments ? C'est un calcul simple ("calculer la géométrie" dans la table attrib)

EDIT pour Daniel : vaut mieux 2 fois qu'une

Bonsoir,

et merci de ces réponses. Oui  c'est une très bonne idée de calculer dans la table attributaire la taille des segments.
Mais n'est-il pas possible d'utiliser cette même fonction en VBA à la suite de la création des segments, ou retourner directement les distances à l'issue d'un Clip ?
Merci
Bonne Soirée.

Bebert

SAlut la Berte !, et salut tout le monde.


Toujours en train de batailler !...

Je te site : "Deuxièmement je suis toujours à la recherche d'une fonction en batch pour mesurer la longueurs de mes segments?"

Et je site daniel :" Je me permet d'ajouter une peut-être réponse au dernier problème de Bébert.  Je suis les l'histoire du début.  Pouquoi ne par se créer une colonne longueur dans la couche de segment et de faire le calcul de longueur?"


Effectivement, la calculette c'est le must ! Mais encore mieux,c'est de l'inclure en automatique dans ton prog la calculette OUI,  mais en vba ... voili voilou...


Ps : pour ta question du bug bizarre... tres bizarre en effet...KML format googleearth...keske ça vient faire la !

A plus jeune!
Longue vie aux casseurs de cailloux

'Fonction qui renvoie un layer a partir de son nom
Function FindLayerByName(pMap As IMap, sName As String) As Ilayer
    Dim i As Integer
    For i = 0 To pMap.LayerCount - 1
            If pMap.layer(i).name = sName Then
            Set FindLayerByName = pMap.layer(i)
        End If
    Next
End Function


Sub launch()
    
    Dim pApp As IApplication
    Set pApp = Application
    Dim pMxDoc As IMxDocument
    Set pMxDoc = pApp.Document
    Dim Mylayer As IFeatureLayer
    Set Mylayer = FindLayerByName(pMxDoc.FocusMap, "MACOUCHE")
    Dim MyFeatureClass As IFeatureClass
    Dim MyTab As ITable
    Set MyTab = Mylayer
    Dim pField As IField
    Dim pFieldEdit As IFieldEdit
    Set MyFeatureClass = Mylayer.FeatureClass
    
    
    'insertion du champ devant recevoir lA LONGUEUR
    Set pField = New Field
    Set pFieldEdit = pField
    With pFieldEdit
        .name = "longueur"
        .Type = esriFieldTypeDouble
    End With
    'Si un nom de champ est deja present -> on passe sur L' exception
    On Error Resume Next
    MyFeatureClass.AddField pField
    On Error GoTo 0
    
    Call CalculatorLong(MyTab, "longueur")
    
End Sub


Private Sub CalculatorLong(ptab As ITable, nom_champ As String)
    ' Declarations
    Dim pCalc As ICalculator
    Dim pCursor As ICursor
    Set pCalc = New Calculator
    Set pCursor = ptab.Update(Nothing, True)
    
    'Expression vba calculette de valeur
    With pCalc
       Set .Cursor = pCursor
                .PreExpression = "dim c as icurve" & vbNewLine & _
                                    "Set c = [Shape]" & vbNewLine & _
                                    "Dim l As Double" & vbNewLine & _
                                    "l = c.length"
                .Expression = "l"
            'on rempli le champ tata
           .Field = nom_champ
    End With
    
  'pEditor.StartOperation
  pCalc.Calculate
End Sub

Et me dit pas qu'il ne marche po ! Certfié, testé sans erreur.

On fire ! Yeah....

A votre service les geomaticiens copains !



PS : pour ton autre probleme, essaye un truc :

Tu vas sur ton editeur vba,
tu click sur le normal.mxt dans l'explorateur de projet
Puis tu choisi module/ arcID,
Enfin, tu essayes de reperer les lignes qui posent probleme, et tu les vires purement et simplement .....

Sinon, essaye aussi de supprimer le normal.mxt pour le regenerer...
il se trouve la(par contre tu vas perdre tes configs par defaut Arcmap..mais bon, vite refait) :
C:\Documents and Settings\lima\Application Data\ESRI\ArcMap\Templates\normal.mxt

Si ça marche tu me dois encore un apéro !
a plouch

Dernière modification par lima (Thu 04 September 2008 20:25)

lol...
Desolé, mais c'était un peu private joke sur les bords... Qd je disais re-apéro, c'est vraiment le cas ...je connais l'annimal...! Mais un pot à chaque resolution pourrai etre sympa aussi moiup ... A la votre!

Hello,

Juste pour info, si ça peut dépanner, je viens de tomber sur un script dans la base arcscript US, pour créer des segments perpendiculaires (si ca peut dépanner) :
http://arcscripts.esri.com/details.asp?dbid=15756

Pas testé mais ca me semblait intéressant dans le cadre de ce sujet.

Trad rapide : Cet outil prend un shp de lignes et génère des lignes perpendiculaires à chaque enregistrement.. Les nouvelles lignes perpendiculaire intersectent les lignes originales au milieu de chaque segment. La géométrie résultante se présente sous forme d'un fichier texte. Pour reconstituer le fichier, on peut utiliser "Créer un shp depuis un fichier txt" de la Toolbox.

Cet outil se présente sous forme d'un script python dont il faut changer les paramètres.

Le fichier entrant original devrait être constitué uniquement d'entités unique (non multiples) idéalement découpé avec un script de découpage aux sommets (vertices)

Robin.