r.stream.extractgrass - En ligne dans le Cloud

PROGRAMME:

Nom

r.stream.extrait - Effectue l'extraction du réseau de flux.

MOTS-CLÉS

raster, hydrologie, réseau hydrographique

SYNOPSIS

r.stream.extrait
r.stream.extrait --Aidez-moi
r.stream.extrait élévation=prénom [accumulation=prénom] [la Dépression=prénom] порог=flotter
[d8couper=flotter] [mexicain=flotter] [longueur_flux=entier] [Mémoire=entier]
[flux_raster=prénom] [vecteur_flux=prénom] [direction=prénom] [--écraser] [--vous aider]
[--verbeux] [--calme] [--ui]

Drapeaux:
--écraser
Autoriser les fichiers de sortie à écraser les fichiers existants

--Aidez-moi
Imprimer le récapitulatif d'utilisation

--verbeux
Sortie du module verbeux

--silencieux
Sortie module silencieuse

--interface utilisateur
Forcer le lancement de la boîte de dialogue GUI

Paramètres:
élévation=prénom [obligatoire]
Nom de la carte raster d'altitude en entrée

accumulation=prénom
Nom de la carte raster d'accumulation en entrée
L'extraction de flux utilisera l'accumulation fournie au lieu de la recalculer

la Dépression=prénom
Nom de la carte raster en entrée avec de vraies dépressions
Les flux ne seront pas acheminés hors des véritables dépressions

порог=flotter [obligatoire]
Accumulation de débit minimum pour les cours d'eau
Doit être > 0

d8couper=flotter
Utiliser SFD au-dessus de ce seuil
Si l'accumulation est supérieure à d8cut, SFD est utilisé à la place de MFD. S'applique seulement si non
carte d'accumulation est donnée.
Valeur par défaut: infini

mexicain=flotter
Exposant de Montgomery pour la pente, désactivé avec 0
Montgomery : l'accumulation est multipliée par pow(slope,mexp) puis comparée à
порог
Valeur par défaut: 0

longueur_flux=entier
Supprimer les segments de flux plus courts que les cellules stream_length
S'applique uniquement aux segments de cours d'eau de premier ordre (sources/têtes de cours d'eau)
Valeur par défaut: 0

Mémoire=entier
Mémoire maximale à utiliser (en Mo)
Taille du cache pour les lignes raster
Valeur par défaut: 300

flux_raster=prénom
Nom de la carte raster en sortie avec des identifiants de flux uniques

vecteur_flux=prénom
Nom de la carte vectorielle de sortie avec des identifiants de flux uniques

direction=prénom
Nom de la carte raster en sortie avec la direction du flux

DESCRIPTION

r.stream.extrait extrait les flux au format raster et vectoriel à partir d'une entrée requise
élévation carte et entrée facultative accumulation carte

NOTES

Cellules NULL (nodata) dans l'entrée élévation map sont ignorés, les valeurs zéro et négatives sont
données d'altitude valides. Les lacunes dans la carte d'élévation qui sont situées dans la zone de
les intérêts doivent être remplis au préalable, par exemple avec r.fillnulls, pour éviter les distorsions.

Toutes les cellules non NULL et non nulles de la Dépression carte seront considérés comme de véritables dépressions.
Les cours d'eau ne seront pas acheminés hors des dépressions. Si une zone est marquée comme dépression mais que le
modèle d'élévation n'a pas de dépression à cet endroit, les ruisseaux ne s'arrêteront pas là. Si un flux
une carte d'accumulation et une carte avec des dépressions réelles sont fournies, la carte d'accumulation de flux
doit correspondre à la carte de dépression de telle sorte que le flux ne soit pas distribué en dehors de la zone indiquée
dépressions. Il est recommandé d'utiliser l'accumulation de débit calculée en interne si un
une carte de dépression est fournie.

Option порог définit la valeur d'accumulation de débit minimale (éventuellement modifiée) qui
lancera un nouveau flux. Si la méthode de Montgomery pour l'initiation du canal est utilisée, le
la valeur de la cellule de la carte d'entrée d'accumulation est multipliée par (tan(pente locale))mexp, puis
par rapport à порог. Si mexicain est donnée que la méthode de Montgomery et
Foufoula-Georgiou (1993) pour initier un flux avec cette valeur. La valeur de la cellule du
la carte d'entrée d'accumulation est multipliée par (tan(pente locale))mexp puis comparée à
порог. Si le seuil est atteint ou dépassé, un nouveau flux est lancé. Le défaut
la valeur 0 désactive Montgomery. Montgomery et Foufoula-Georgiou (1993) recommandent généralement
d'utiliser 2.0 comme exposant. mexicain des valeurs plus proches de 0 produiront des flux plus similaires à
flux extraits avec Montgomery désactivé. Plus grand mexicain les valeurs diminuent le nombre de
cours d'eau dans les zones plates et augmenter le nombre de cours d'eau dans les zones escarpées. Si poids is
donné, le poids est appliqué en premier.

Option d8couper définit la quantité minimale d'écoulement de surface (accumulation) lorsque le SFD (D8) sera
utilisé à la place du MFD (FD8) pour calculer l'accumulation de débit. S'applique uniquement si aucun cumul
la carte est fournie. Le réglage sur 0 désactive complètement le MFD.

Option longueur_flux définit la longueur minimale du flux en nombre de cellules pour le premier ordre
segments de cours d'eau (tête/source). Tous les segments de flux de premier ordre plus courts que longueur_flux
sera supprimé.

Sortie direction la carte raster contient la direction du flux pour toutes les cellules non NULL en entrée
élévation. Le sens d'écoulement est de type D8 avec une plage de 1 à 8. Multiplication des valeurs par
45 donne des degrés CCW de l'Est. La direction du flux a été ajustée pendant l'éclaircissage, en prenant
raccourcis et sauts de cellules qui ont été éliminés par la procédure d'amincissement.

Discussions extraction
Sinon accumulation une carte d'entrée est fournie, l'accumulation de flux est déterminée avec un
analyse hydrologique similaire à bassin versant. L'algorithme est MFD (FD8) d'après Holmgren
1994, comme pour bassin versantL’ порог l'option détermine le nombre de flux et le détail
des réseaux de flux. Chaque fois que l'accumulation de débit atteint порог, un nouveau flux est lancé
et tracée en aval jusqu'à son point de sortie. Pour ce qui est de bassin versant, l'accumulation de débit est
calculé comme le nombre de cellules s'écoulant à travers une cellule.

If accumulation est donné que les valeurs d'accumulation des accumulation carte sont
utilisé et non calculé à partir de l'entrée élévation carte. Dans ce cas le élévation la carte doit
être exactement la même carte utilisée pour calculer accumulation. Si accumulation a été calculé
avec r.terraflow, la sortie d'élévation remplie de r.terraflow doit être utilisé. Plus loin, le
la région actuelle doit être alignée sur la accumulation carte. La direction du flux est la première
calculé à partir de élévation puis ajusté à accumulation. Il n'est pas nécessaire de
fournir accumulation comme le nombre de cellules, il peut également être le éventuellement ajusté ou
superficie contributive totale pondérée en mètres carrés ou toute autre unité. Quand un flux original
carte d'accumulation est ajustée ou pesée, l'ajustement ou la pesée ne doit pas convertir
valeurs d'accumulation valides aux valeurs NULL (nodata).

Pesé flux accumulation
L'accumulation de débit peut être calculée en premier, par ex. bassin versant, puis modifié avant
l'utiliser comme entrée pour r.stream.extrait. Dans sa forme générale, une carte d'accumulation pondérée est
généré en créant d'abord une carte de pesée puis en multipliant la carte d'accumulation avec
la carte de pesée à l'aide r.mapcalc. Il est fortement recommandé d'évaluer le débit pesé
carte d'accumulation d'abord, avant de l'utiliser comme entrée pour r.stream.extrait.

Cela permet par exemple de diminuer le nombre de cours d'eau dans les zones sèches et d'augmenter le nombre de
cours d'eau dans les zones humides en fixant poids à plus petit que 1 dans les zones sèches et plus grand que 1 dans
zones humides.

Une autre possibilité consiste à restreindre l'initiation du canal aux vallées déterminées à partir du terrain
morphologie. Les vallées peuvent être déterminées avec r.param.échelle param=crosc (transversal ou
courbure tangentielle). Les valeurs de courbure < 0 indiquent des caractéristiques concaves, c'est-à-dire des vallées. Les
la taille de la fenêtre de traitement détermine si des vallées étroites ou larges seront
identifié (voir exemple ci-dessous).

Définir a région of intérêt
La procédure d'extraction de flux peut être limitée à une certaine région d'intérêt, par exemple un
sous-bassin, en définissant la région de calcul avec g.région et/ou créer un MASQUE. Tel
région d'intérêt doit être un bassin versant complet, complet en ce sens que la
une zone complète en amont d'un point de sortie est incluse et tamponnée avec au moins une cellule.

Discussions sortie
Le raster et le vecteur en sortie contiennent des segments de flux avec des ID uniques. Notez que ces identifiants
sont différents des identifiants attribués par bassin versant. La sortie vectorielle contient également des points
à l'emplacement du début d'un segment de cours d'eau, aux confluents et au réseau de cours d'eau
emplacements des points de vente.

Sortie flux_raster la carte raster stocke les flux extraits. Les valeurs des cellules encodent un identifiant unique
pour chaque segment de flux.

Sortie vecteur_flux la carte vectorielle stocke les segments et points de flux extraits. Les points sont
écrit à l'emplacement de départ de chaque segment de flux et à la sortie d'un flux
réseau. Dans la couche 1, les catégories sont des identifiants uniques, identiques à la valeur de cellule du raster
sortir. La table attributaire de la couche 1 contient des informations sur le type de flux
segment : segment de départ, ou segment intermédiaire avec affluents. Les colonnes sont cat int,
stream_type varchar(), type_code int. L'encodage pour type_code est 0 = début, 1 =
intermédiaire. Dans la couche 2, les catégories sont identiques à type_code dans la couche 1 avec des
catégorie 2 = sortie pour points de sortie. Points de catégorie 1 = intermédiaire dans la couche 2
sont à l'endroit des confluents.

EXEMPLE

Cet exemple est basé sur la carte d'altitude "elev_ned_30m" dans l'échantillon de Caroline du Nord
jeu de données et utilise les vallées déterminées avec r.param.échelle peser une carte d'accumulation
produit avec bassin versant.
# région définie
g.region -p raster=elev_ned_30m@PERMANENT
# calculer l'accumulation de flux
r.bassin versant ele=elev_ned_30m@PERMANENT acc=elevation.10m.acc
# courbure pour obtenir des vallées étroites
r.param.scale input=elev_ned_30m@PERMANENT output=tangential_curv_5 size=5 param=crosc
# courbure pour obtenir des vallées un peu plus larges
r.param.scale input=elev_ned_30m@PERMANENT output=tangential_curv_7 size=7 param=crosc
# courbure pour obtenir de larges vallées
r.param.scale input=elev_ned_30m@PERMANENT output=tangential_curv_11 size=11 param=crosc
# créer une carte de poids
r.mapcalc "poids = if(tangential_curv_5 < 0, -100 * tangential_curv_5, \
if(tangential_curv_7 < 0, -100 * tangential_curv_7, \
if(tangential_curv_11 < 0, -100 * tangential_curv_11, 0.000001)))"
# carte d'accumulation de poids
r.mapcalc expr="elev_ned_30m.acc.weighed = elev_ned_30m.acc * poids"
# copie la table des couleurs de la carte d'accumulation d'origine
r.colors map=elev_ned_30m.acc.weighed raster=elev_ned_30m.acc
Affichez à la fois l'original et la carte d'accumulation pesée. Comparez-les et continuez si
la carte d'accumulation pesée a du sens.
# extraire les flux
r.stream.extract altitude=elev_ned_30m@PERMANENT \
accumulation=elev_ned_30m.acc.pesé \
seuil=1000 \
stream_rast=elev_ned_30m.streams
# extraire les flux à l'aide de la carte d'accumulation d'origine
r.stream.extract altitude=elev_ned_30m@PERMANENT \
cumul=elev_ned_30m.acc \
seuil=1000 \
stream_rast=elev_ned_30m.streams.noweight
Affichez maintenant les deux cartes de flux et décidez laquelle est la plus réaliste.

Références

· Ehlschläger, C. (1989). En utilisant le AT Rechercher Algorithme à Développer Hydrologique
Modèles photo de Ressources Élévation Sauvegarde de, Procédures of INTERNATIONAL Portée
Info placement (IGIS) Symposium «89, pp 275-281 (Baltimore, MD, 18-19 mars
1989). URL : http://faculty.wiu.edu/CR-Ehlschlaeger2/older/IGIS/paper.html

· Holmgren, P. (1994). Multiple flux direction algorithmes pour ruissellement la modélisation in
grille basé élévation des modèles: An empirique évaluation. Hydrologique Processus Vol
8(4), pages 327-334. DOI : 10.1002/hyp.3360080405

· Montgomery, DR, Foufoula-Georgiou, E. (1993). Développement réseau source
représentation en utilisant numérique élévation . L'eau Ressources Recherche Vol
29(12), pages 3925-3934.

Utilisez r.stream.extractgrass en ligne en utilisant les services onworks.net