r.stream.extractgrass - Online sa Cloud

PROGRAMA:

NAME

r.stream.extract - Nagsasagawa ng stream network extraction.

KEYWORDS

raster, hydrology, stream network

SINOPSIS

r.stream.extract
r.stream.extract - Tumulong
r.stream.extract elevation=pangalan [akumulasyon=pangalan] [depresyon=pangalan] threshold=lumutang
[d8cut=lumutang] [mexp=lumutang] [stream_length=kabuuan] [memorya=kabuuan]
[stream_raster=pangalan] [stream_vector=pangalan] [utos=pangalan] [--patungan] [--Tulungan]
[--pandiwang] [--tahimik] [--ui]

Mga Bandila:
--patungan
Pahintulutan ang mga output file na i-overwrite ang mga kasalukuyang file

- Tumulong
I-print ang buod ng paggamit

--verbose
Verbose na output ng module

--tahimik
Tahimik na output ng module

--ui
Piliting ilunsad ang dialog ng GUI

parameter:
elevation=pangalan [kailangan]
Pangalan ng input elevation raster map

akumulasyon=pangalan
Pangalan ng input accumulation raster map
Gagamitin ng stream extraction ang ibinigay na akumulasyon sa halip na kalkulahin itong muli

depresyon=pangalan
Pangalan ng input raster na mapa na may mga totoong depression
Ang mga stream ay hindi iruruta sa mga totoong depression

threshold=lumutang [kailangan]
Minimum na akumulasyon ng daloy para sa mga stream
Dapat ay > 0

d8cut=lumutang
Gamitin ang SFD sa itaas ng threshold na ito
Kung ang akumulasyon ay mas malaki kaysa sa d8cut, SFD ang ginagamit sa halip na MFD. Nalalapat lamang kung hindi
ibinibigay ang accumulation map.
Default: kawalang-hanggan

mexp=lumutang
Montgomery exponent para sa slope, hindi pinagana na may 0
Montgomery: ang akumulasyon ay pinarami ng pow(slope,mexp) at pagkatapos ay inihambing sa
threshold
Default: 0

stream_length=kabuuan
Tanggalin ang mga segment ng stream na mas maikli kaysa sa mga cell ng stream_length
Nalalapat lang sa mga first-order na segment ng stream (mga spring/stream head)
Default: 0

memorya=kabuuan
Pinakamataas na memorya na gagamitin (sa MB)
Laki ng cache para sa mga raster row
Default: 300

stream_raster=pangalan
Pangalan para sa output raster map na may mga natatanging stream id

stream_vector=pangalan
Pangalan para sa output vector map na may mga natatanging stream id

utos=pangalan
Pangalan para sa output raster map na may direksyon ng daloy

DESCRIPTION

r.stream.extract kinukuha ang mga stream sa parehong format ng raster at vector mula sa isang kinakailangang input
elevation mapa at opsyonal na input akumulasyon mapa.

NOTA

NULL (nodata) na mga cell sa input elevation ang mapa ay binabalewala, ang zero at ang mga negatibong halaga ay
wastong data ng elevation. Mga gaps sa elevation map na matatagpuan sa loob ng lugar ng
dapat punan muna ang interes, hal r.fillnulls, upang maiwasan ang mga pagbaluktot.

Lahat ng non-NULL at non-zero na mga cell ng depresyon ang mapa ay ituturing na tunay na mga depresyon.
Ang mga stream ay hindi iruruta mula sa mga depresyon. Kung ang isang lugar ay minarkahan bilang depresyon ngunit ang
Ang modelo ng elevation ay walang depresyon sa lokasyong ito, ang mga batis ay hindi titigil doon. Kung isang daloy
akumulasyon mapa at isang mapa na may tunay na depressions ay ibinigay, ang daloy ng akumulasyon mapa
dapat tumugma sa mapa ng depresyon upang ang daloy ay hindi naipamahagi sa labas ng ipinahiwatig
mga depresyon. Inirerekomenda na gumamit ng internally computed flow accumulation kung a
Ang mapa ng depresyon ay ibinigay.

Opsyon threshold tumutukoy sa minimum (opsyonal na binago) na halaga ng akumulasyon ng daloy
magsisimula ng bagong stream. Kung gagamitin ang pamamaraan ni Montgomery para sa pagsisimula ng channel, ang
cell value ng accumulation input map ay pinarami ng (tan(lokal na slope))mexp at pagkatapos
kumpara sa threshold. Kung mexp ay ibinigay kaysa sa paraan ng Montgomery at
Foufoula-Georgiou (1993) upang simulan ang isang stream na may ganitong halaga. Ang halaga ng cell ng
ang accumulation input map ay pinarami ng (tan(lokal na slope))mexp at pagkatapos ay inihambing sa
threshold. Kung naabot o nalampasan ang threshold, magsisimula ang isang bagong stream. Ang default
hindi pinapagana ng value 0 ang Montgomery. Karaniwang inirerekomenda ni Montgomery at Foufoula-Georgiou (1993).
gamitin ang 2.0 bilang exponent. mexp ang mga value na mas malapit sa 0 ay gagawa ng mga stream na mas katulad ng
mga stream na nakuha nang hindi pinagana ang Montgomery. Mas malaki mexp binabawasan ng mga halaga ang bilang ng
batis sa mga patag na lugar at dagdagan ang bilang ng mga batis sa matarik na lugar. Kung timbang is
ibinigay, ang timbang ay inilapat muna.

Opsyon d8cut tumutukoy sa pinakamababang halaga ng daloy sa lupa (akumulasyon) kung kailan magiging SFD (D8).
ginamit sa halip na MFD (FD8) upang kalkulahin ang akumulasyon ng daloy. Nalalapat lamang kung walang akumulasyon
ibinigay ang mapa. Ang pagtatakda sa 0 ay ganap na hindi pinapagana ang MFD.

Opsyon stream_length tumutukoy sa pinakamababang haba ng stream sa bilang ng mga cell para sa first-order
(ulo/tagsibol) mga segment ng stream. Lahat ng first-order na mga segment ng stream na mas maikli kaysa stream_length
tatanggalin

Pagbubuhos utos Ang mapa ng raster ay naglalaman ng direksyon ng daloy para sa lahat ng mga non-NULL na cell sa input
elevation. Ang direksyon ng daloy ay uri ng D8 na may hanay na 1 hanggang 8. Pagpaparami ng mga halaga sa
45 ay nagbibigay ng digri CCW mula sa Silangan. Ang direksyon ng daloy ay nababagay sa panahon ng pagnipis, pagkuha
mga shortcut at paglaktaw ng mga cell na inalis ng pamamaraan ng pagnipis.

sapa pagkuha
Kung hindi akumulasyon ang input map ay ibinigay, ang daloy ng accumulation ay tinutukoy gamit ang a
hydrological analysis katulad ng r.watershed. Ang algorithm ay MFD (FD8) pagkatapos ng Holmgren
1994, para sa r.watershed. ang threshold Tinutukoy ng opsyon ang bilang ng mga stream at detalye
ng mga stream network. Sa tuwing umabot ang akumulasyon ng daloy threshold, nagsimula ang isang bagong stream
at natunton sa ibaba ng agos hanggang sa labasan nito. Tungkol naman sa r.watershed, ang akumulasyon ng daloy ay
kinakalkula bilang ang bilang ng mga cell na dumadaloy sa isang cell.

If akumulasyon ay ibinigay kaysa sa mga halaga ng akumulasyon ng ibinigay akumulasyon mapa ay
ginamit at hindi kinakalkula mula sa input elevation mapa. Sa kasong ito ang elevation mapa dapat
maging eksakto ang parehong mapa na ginamit sa pagkalkula akumulasyon. Kung akumulasyon ay kalkulado
sa r.terraflow, ang punong elevation na output ng r.terraflow dapat gamitin. Karagdagan pa, ang
kasalukuyang rehiyon ay dapat na nakahanay sa akumulasyon mapa. Ang direksyon ng daloy ay una
kinakalkula mula sa elevation at saka nag-adjust sa akumulasyon. Hindi naman kailangan
magbigay akumulasyon bilang bilang ng mga cell, maaari rin itong maging opsyonal na i-adjust o
tinimbang ang kabuuang lugar na nag-aambag sa metro kuwadrado o anumang iba pang yunit. Kapag ang isang orihinal na daloy
akumulasyon mapa ay nababagay o weighed, ang pagsasaayos o pagtimbang ay hindi dapat i-convert
wastong mga halaga ng akumulasyon sa mga halaga ng NULL (nodata).

Tinimbang pag-agos akumulasyon
Maaaring kalkulahin muna ang akumulasyon ng daloy, hal r.watershed, at pagkatapos ay binago bago
gamit ito bilang input para sa r.stream.extract. Sa pangkalahatang anyo nito, ang isang weighed accumulation map ay
nabuo sa pamamagitan ng unang paggawa ng weighing map at pagkatapos ay pagpaparami ng accumulation map sa
ang weighing map gamit r.mapcalc. Lubhang inirerekomenda na suriin ang natimbang na daloy
accumulation map muna, bago ito gamitin bilang input para sa r.stream.extract.

Nagbibigay-daan ito hal. upang bawasan ang bilang ng mga sapa sa mga tuyong lugar at dagdagan ang bilang ng
batis sa mga basang lugar sa pamamagitan ng pagtatakda timbang sa mas maliit sa 1 sa mga tuyong lugar at mas malaki sa 1 in
mga basang lugar.

Ang isa pang posibilidad ay paghigpitan ang pagsisimula ng channel sa mga lambak na tinutukoy mula sa lupain
morpolohiya. Ang mga lambak ay maaaring matukoy gamit ang r.param.scale param=cross (cross-sectional o
tangential curvature). Ang mga halaga ng curvature < 0 ay nagpapahiwatig ng mga malukong tampok, ibig sabihin, mga lambak. Ang
ang laki ng window sa pagpoproseso ay tumutukoy kung magiging makitid o malalawak na lambak
nakilala (Tingnan ang halimbawa sa ibaba).

Pagtukoy a rehiyon of interes
Ang pamamaraan ng pagkuha ng stream ay maaaring paghigpitan sa isang partikular na rehiyon ng interes, hal a
subbasin, sa pamamagitan ng pagtatakda ng computational region na may g.rehiyon at/o paggawa ng MASK. ganyan
rehiyon ng interes ay dapat na isang kumpletong catchment area, kumpleto sa kahulugan na ang
Ang kumpletong lugar sa itaas ng agos ng isang outlet point ay kasama at na-buffer na may hindi bababa sa isang cell.

sapa output
Ang output raster at vector ay naglalaman ng mga segment ng stream na may mga natatanging ID. Tandaan na ang mga ID na ito
ay iba sa mga ID na itinalaga ni r.watershed. Ang output ng vector ay naglalaman din ng mga puntos
sa lokasyon ng pagsisimula ng isang stream segment, sa mga confluence at sa stream network
mga lokasyon ng outlet.

Pagbubuhos stream_raster Ang raster map ay nag-iimbak ng mga nakuhang stream. Ang mga cell value ay nag-encode ng isang natatanging ID
para sa bawat segment ng stream.

Pagbubuhos stream_vector Ang vector map ay nag-iimbak ng mga na-extract na mga segment at puntos ng stream. Ang mga puntos ay
nakasulat sa panimulang lokasyon ng bawat stream segment at sa labasan ng isang stream
network. Sa layer 1, ang mga kategorya ay mga natatanging ID, na kapareho ng cell value ng raster
output. Ang talahanayan ng katangian para sa layer 1 ay naglalaman ng impormasyon tungkol sa uri ng stream
segment: simulang segment, o intermediate na segment na may mga tributaries. Ang mga column ay cat int,
stream_type varchar(), type_code int. Ang pag-encode para sa type_code ay 0 = simula, 1 =
nasa pagitan. Sa layer 2, ang mga kategorya ay kapareho ng type_code sa layer 1 na may karagdagang
kategorya 2 = outlet para sa mga outlet point. Mga puntos na may kategorya 1 = intermediate sa layer 2
ay nasa lokasyon ng mga confluence.

Halimbawa

Ang halimbawang ito ay batay sa elevation map na "elev_ned_30m" sa sample ng North Carolina
dataset at gumagamit ng mga lambak na tinutukoy sa r.param.scale upang timbangin ang isang mapa ng akumulasyon
ginawa ng r.watershed.
# itakda ang rehiyon
g.region -p raster=elev_ned_30m@PERMANENT
# kalkulahin ang akumulasyon ng daloy
r.watershed ele=elev_ned_30m@PERMANENT acc=elevation.10m.acc
# curvature para makakuha ng makikitid na lambak
r.param.scale input=elev_ned_30m@PERMANENT na output=tangential_curv_5 size=5 param=cross
# curvature upang makakuha ng medyo mas malawak na mga lambak
r.param.scale input=elev_ned_30m@PERMANENT na output=tangential_curv_7 size=7 param=cross
# curvature para makakuha ng malalawak na lambak
r.param.scale input=elev_ned_30m@PERMANENT na output=tangential_curv_11 size=11 param=cross
# lumikha ng mapa ng timbang
r.mapcalc "weight = 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)))"
# timbangin ang akumulasyon na mapa
r.mapcalc expr="elev_ned_30m.acc.weighed = elev_ned_30m.acc * weight"
# kopyahin ang talahanayan ng kulay mula sa orihinal na mapa ng akumulasyon
r.colors map=elev_ned_30m.acc.weighed raster=elev_ned_30m.acc
Ipakita ang parehong orihinal at ang weighed accumulation map. Ihambing ang mga ito at magpatuloy kung
ang weighed accumulation map ay may katuturan.
# extract stream
r.stream.extract elevation=elev_ned_30m@PERMANENT \
accumulation=elev_ned_30m.acc.weighed \
threshold=1000 \
stream_rast=elev_ned_30m.streams
# extract stream gamit ang orihinal na accumulation map
r.stream.extract elevation=elev_ned_30m@PERMANENT \
accumulation=elev_ned_30m.acc \
threshold=1000 \
stream_rast=elev_ned_30m.streams.noweight
Ipakita ngayon ang parehong mga stream ng mapa at magpasya kung alin ang mas makatotohanan.

Mga sanggunian

· Ehlschlaeger, C. (1989). paggamit ang AT Maghanap Algorithm sa Magbuo Hydrologic
Modellen mula Digital Kataasan data, Mga pamamaraan of Internasyonal Heograpiko
Impormasyon Systems (IGIS) panayam '89, pp 275-281 (Baltimore, MD, 18-19 Marso
1989). URL: http://faculty.wiu.edu/CR-Ehlschlaeger2/older/IGIS/paper.html

· Holmgren, P. (1994). Maramihang pag-agos utos algorithm para patakbuhan pagmomolde in
parilya batay elevation mga modelo: An empirical pagsusuri. Hydrological Mga Proseso Vol
8(4), pp 327-334. DOI: 10.1002/hyp.3360080405

· Montgomery, DR, Foufoula-Georgiou, E. (1993). channel network pinagmulan
pagkatawan paggamit digital elevation na mga modelo. tubig Mga mapagkukunan Pananaliksik Vol
29(12), pp 3925-3934.

Gumamit ng r.stream.extractgrass online gamit ang mga serbisyo ng onworks.net