` `
하천망 분석¶
D8 오르막 극값¶
설명¶
입력 그리드에서 D8 유수 모델을 바탕으로 오르막 극값(최대 또는 최소 가운데 하나)을 판정합니다. 이 알고리즘은 원래 (낙차 분석에 따른) 최적의 하천망을 산출하는 경사와 면적의 곱의 한계값을 식별하기 위한 하천 래스터를 생성하기 위해 개발됐습니다.
부가적인 배출 포인트 shapefile을 사용하는 경우, 배출 셀과 해당 셀의 (D8 유수 모델에 따른) 오르막 셀들이 있는 영역만 평가합니다.
기본적으로, 이 도구는 경계 오염을 확인합니다. 경계 오염이란 영역 바깥에 있는 그리드 셀을 계산에 넣지 않기 때문에 결과가 저평가될 수도 있다는 가능성으로 정의됩니다. 경계 오염은 배수가 경계로부터 안 쪽으로 흐르거나 또는 표고값이 NODATA인 영역에서 발생합니다. 이 알고리즘은 이를 인식하고 이 그리드 셀들의 결과에 대해 NODATA를 보고합니다. 경계선에서 영역으로 들어오는 유수 경로를 따라 NODATA 값들이 쭉 나열되는 것을 흔히 볼 수 있습니다. 이는 기대대로의 효과로서, 해당 그리드 셀들의 결과가 사용할 수 있는 데이터 영역 외부에 있는 지형에 의존하기 때문에 그 기여 영역을 알 수 없다는 사실을 나타냅니다. 사용자가 경계 오염이 문제가 되지 않는다는 사실을 알고 있거나 또는 이런 문제를 무시하고자 하는 경우, 예를 들면 유역 외곽선을 따라 DEM을 잘라낸 경우라면 경계 오염 확인 옵션을 끌 수도 있습니다.
파라미터¶
- D8 Flow Directions Grid [raster]
각 셀 별로 그 인접 또는 대각선 방향의 이웃 8개 셀 가운데 가장 심한 내리막 경사를 가지고 있는 셀의 방향으로 정의한 D8 유수 방향 그리드입니다. 이 그리드는 D8 유수 방향 알고리즘의 산출물로 얻을 수 있습니다.
- Upslope Values Grid [raster]
오르막 최대값 또는 최소값을 선택한 값의 그리드입니다. 가장 흔히 쓰이는 값은 낙차 분석에 따라 하천 래스터를 생성하는 경우 필요한 경사와 면적의 곱입니다.
- Outlets Shapefile [vector: point]
Optional.
관심 배출 포인트를 정의하는 포인트 shapefile입니다. 이 파일을 입력하는 경우, 알고리즘이 해당 배출 셀의 주변 오르막 영역만 평가할 것입니다.
- Check for edge contamination [boolean]
알고리즘이 경계 오염을 확인할지 여부를 나타내는 옵션입니다.
기본값: True
- Use max upslope value [boolean]
오르막 경사의 최대값 또는 최소값 가운데 어느 쪽을 계산할지 나타내는 옵션입니다.
기본값: True
산출물¶
- Extreme Upslope Values Grid [raster]
오르막 경사의 최대값/최소값 그리드입니다.
콘솔 사용법¶
processing.runalg('taudem:d8extremeupslopevalue', -p, -sa, -o, -nc, -min, -ssa)
참고¶
길이 면적 하천원¶
설명¶
오르막 경로 길이, D8 기여 영역 그리드 입력, 그리고 파라미터 M 과 y 를 바탕으로 A >= (M)(Ly) 를 평가하는 지표(indicator) 그리드 (1, 0)을 생성합니다. 이 그리드는 하천원(stream source)일 가능성이 높은 그리드 셀들을 표시합니다. 이 알고리즘은 실험적인 메소드로, 하천 L ~ A 0.6 에 대해 앞의 내용을 서술하고 있는 핵의 법칙 을 그 이론적 기반으로 삼고 있습니다. 하지만 산사면(hillslope)의 경우 평행 유수 L ~ A 를 대상으로 합니다. 따라서 산사면에서 하천으로의 전이를 L ~ A 0.8 로 표현할 수도 있는데, 그렇다면 그리드 셀이 A > M (L (1/0.8)) 인 경우 하천 셀이라고 식별할 수 있을 것입니다.
파라미터¶
- Length Grid [raster]
각 셀 별 최장 오르막 길이를 담고 있는 그리드입니다. 최장 오르막 길이는 각 셀로 흘러드는 셀 가운데 가장 멀리 있는 셀에서 시작되는 유수 경로의 길이로 계산합니다. 이 길이는 셀 크기 및 방향의 인접/대각선 여부를 고려해서 셀 중심점 간의 거리로 측정합니다. 어떤 셀을 하천 셀로 간주할 것인지 결정하기 위한 A >(M)(Ly) 공식에 사용되는 길이 L 이 바로 이 최장 오르막 길이입니다. 이 그리드는 그리드 망 알고리즘의 산출물입니다.
- Contributing Area Grid [raster]
D8 알고리즘을 사용해서 계산한 각 셀 별 기여 영역 값 그리드입니다. 어느 셀의 기여 영역 값은 해당 셀 고유의 기여와 해당 셀로 흘러드는 모든 오르막 이웃에서 나오는 기여를 더한 값으로, 셀 개수로 측정됩니다. 일반적으로 이 그리드는 D8 기여 영역 알고리즘의 산출물입니다. 길이 면적 하천원 알고리즘에서 하천으로의 전이 여부를 결정하는 A > (M)(Ly) 공식 안의 비교 대상인 기여 영역 A 가 바로 이 기여 영역 값입니다.
- Threshold [number]
하천의 시작을 식별하는 공식 A > (M)(Ly) 에 사용되는 곱수 한계값 M 파라미터입니다.
기본값: 0.03
- Exponent [number]
하천의 시작을 식별하는 공식 A > (M)(Ly) 에 사용되는 지수(exponent) y 파라미터입니다. 분기 체계에서, 핵의 법칙은 1/y = 0.6 (또는 0.56)인 L = 1/M A(1/y) 를 제안하고 있습니다. (y 가 약 1.7입니다.) 평행 유수 체계에서 L 은 A 에 비례합니다. (y 가 약 1입니다.) 이 메소드는 그 사이 어딘가에 있는 지수 y 를 사용해서 이 두 패러다임 사이의 전이를 식별하고자 합니다. (y 가 약 1.3입니다.)
기본값: 1.3
산출물¶
- Stream Source Grid [raster]
최장 오르막 경로 길이, D8 기여 영역 그리드 입력, 그리고 파라미터 M 과 y 를 바탕으로 A >= (M)(L^y) 를 평가하는 지표(indicator) 그리드 (1, 0)입니다. 이 그리드는 하천원(stream source)일 가능성이 높은 그리드 셀들을 표시합니다.
콘솔 사용법¶
processing.runalg('taudem:lengthareastreamsource', length_grid, contrib_area_grid, threshold, exponent, stream_source_grid)
참고¶
배출 포인트를 하천으로 이동¶
설명¶
하천 래스터 그리드에서 하천 셀과 정렬되지 않은 배출 포인트를 내리막 D8 유수 방향을 따라 하천 래스터 셀과 만날 때까지, 또는 그리드 셀들의 max_dist 개수가 판정될 때까지, 또는 유수 경로가 영역을 벗어날 때까지 (예를 들어 D8 유수 방향의 NODATA 값이 나올 때까지) 이동시킵니다. 이 알고리즘의 산출 파일은 새로운 배출 shapefile로, 가능한 경우 각 배출 포인트를 하천 래스터 그리드와 일치하도록 이동시킨 shapefile입니다. 이 새로운 shapefile에는 각 배출 포인트가 변경된 내용을 나타내는 dist_moved 필드가 추가돼 있습니다. 이미 하천 셀 위에 있는 포인트들은 이동되지 않으며, 해당 포인트의 dist_moved 필드엔 0 값이 할당됩니다. 초기에 하천 셀 위에 없던 포인트들은 내리막 D8 유수 방향을 따라 다음 사항 가운데 하나가 일어날 때까지 이동됩니다: 1) 그리드 셀의 max_dist 개수를 넘기 전에 하천 그리드 셀과 만나는 경우, 포인트를 이동시키고 dist_moved 필드에 포인트가 몇 개의 그리드 셀들을 이동했는지 나타내는 값을 할당합니다. 2) 그리드 셀의 max_dist 개수를 넘거나, 또는 3) 영역 밖으로 나가는 경우 (예를 들어 D8 유수 방향의 NODATA 값을 만나는 경우) 배출 포인트를 이동시키지 않고 dist_moved 필드에 -1 값을 할당합니다.
파라미터¶
- D8 Flow Direction Grid [raster]
각 셀 별로 그 인접 또는 대각선 방향의 이웃 8개 셀 가운데 가장 심한 내리막 경사를 가지고 있는 셀의 방향으로 정의한 D8 유수 방향 그리드입니다. 이 그리드는 D8 유수 방향 알고리즘의 산출물로 얻을 수 있습니다.
- Stream Raster Grid [raster]
이 입력물은 하천인 경우 그리드 셀 값을 1, 나머지 셀 값을 0으로 해서 하천의 위치를 나타내는 지표(indicator) 그리드 (1, 0)입니다. 그리드입니다. 하천망 분석 도구 모음 가운데 하나로 이런 그리드를 생성할 수 있습니다.
- Outlets Shapefile [vector: point]
관심 포인트 또는 배출 포인트를 정의하는 포인트 shapefile입니다. 이상적으로는 이런 포인트들이 하천 위에 위치해 있어야 하지만, shapefile 포인트 위치가 하천 래스터 그리드에 대해 정확히 등록돼 있지 않을 수도 있기 때문에 하천 위에 정확히 위치하지 않을 수도 있습니다.
- Maximum Number of Grid Cells to traverse [number]
이 입력 파라미터는 입력 배출 shapefile 안에 있는 포인트들이 산출 배출 shapefile로 저장되지 전에 이동할 수 있는 최대 그리드 셀 개수입니다.
기본값: 50
산출물¶
- Output Outlet Shapefile [vector]
관심 포인트 또는 배출 포인트를 정의하는 포인트 shapefile입니다. 이 파일이 보유한 포인트 개수는 입력 배출 shapefile 이 보유한 포인트 개수와 동일합니다. 입력 포인트가 하천 위에 존재하는 경우, 포인트는 이동하지 않습니다. 입력 포인트가 하천 위에 존재하지 않는 경우, 포인트를 내리막 D8 유수 방향을 따라 하천을 만날 때까지 또는 최장 거리까지 이동시킵니다. dist_moved 필드가 이 파일에 추가되는데, 포인트가 이동한 셀 개수를 담고 있습니다. 포인트가 원래 하천 위에 위치한 경우 이 필드의 값은 0이며, 하천이 최장 거리 내에 없기 때문에 이동하지 않았을 경우 -1이고 포인트가 이동한 경우 이동한 셀 개수의 값입니다.
콘솔 사용법¶
processing.runalg('taudem:moveoutletstostreams', -p, -src, -o, -md, -om)
참고¶
더글러스-패커¶
설명¶
더글러스-패커 알고리즘 이 생성한 오르막 만곡 그리드 셀들의 지표(indicator) 그리드 (1, 0)을 생성합니다.
이 알고리즘은 먼저 DEM을 중심, 테두리, 그리고 대각선 방향에 가중치를 부여하는 커널을 통해 평활화합니다. 그 다음 Peuker & Douglas (1975) 방법을 사용해서 오르막 만곡 그리드 셀을 식별합니다. (이 방법은 Band (1986)에서도 설명하고 있습니다.) 이 기술은 전체 그리드를 표시한 다음, 그리드 셀 4개로 이루어진 각 사분면을 한번에(single pass) 평가해서 가장 높은 셀의 표시를 지웁니다. 나머지 표시된 셀들을 “오르막 만곡”으로 간주하는데, 이를 살펴보면 수로망과 닮아 있습니다. 이 원시 수로망은 일반적으로 연결성이 부족해서 솎아내기(thinning)를 해야 합니다. Band (1986)에서 이런 문제들을 자세하게 논의하고 있습니다.
파라미터¶
- Elevation Grid [raster]
표고 값을 보유한 그리드를 지정합니다. 일반적으로 이 그리드는 구덩이 제거 도구의 산출물로, 구덩이를 제거한 표고입니다.
- Center Smoothing Weight [number]
이 알고리즘이 오르막 만곡 그리드 셀을 식별하기 전에 커널이 DEM을 평활화하기 위해 사용하는 중심 가중치 파라미터입니다.
기본값: 0.4
- Side Smoothing Weight [number]
이 알고리즘이 오르막 만곡 그리드 셀을 식별하기 전에 커널이 DEM을 평활화하기 위해 사용하는 테두리 가중치 파라미터입니다.
기본값: 0.1
- Diagonal Smoothing Weight [number]
이 알고리즘이 오르막 만곡 그리드 셀을 식별하기 전에 커널이 DEM을 평활화하기 위해 사용하는 대각선 방향 가중치 파라미터입니다.
기본값: 0.05
산출물¶
- Stream Source Grid [raster]
더글러스-패커 알고리즘이 생성한 오르막 만곡 그리드 셀들의 지표(indicator) 그리드 (1, 0)로, 이를 살펴보면 수로망과 닮아 있습니다. 이 원시 수로망은 일반적으로 연결성이 부족해서 솎아내기(thinning)를 해야 합니다. Band (1986)에서 이런 문제들을 자세하게 논의하고 있습니다.
콘솔 사용법¶
processing.runalg('taudem:peukerdouglas', elevation_grid, center_weight, side_weight, diagonal_weight, stream_source_grid)
참고¶
- Band, L. E., (1986), “Topographic partition of watersheds with digital elevation models”, Water Resources Research, 22(1): 15-24.
- Peuker, T. K. and D. H. Douglas, (1975), “Detection of surface-specific points by local parallel processing of discrete terrain elevation data”, Comput. Graphics Image Process., 4: 375-387.
경사-면적 조합¶
설명¶
경사, 특정 집수 영역 그리드 입력물, 그리고 m 과 n 파라미터를 바탕으로 하는 경사-면적 값 (Sm) (An) 을 담고 있는 그리드를 생성합니다. 이 알고리즘은 경사-면적 하천 래스터 묘사 방법의 일부로 사용하기 위해 개발됐습니다.
파라미터¶
- Slope Grid [raster]
경사 값을 가진 입력 그리드입니다. 이 그리드는 고저차를 거리로 나눈 값으로 경사를 측정하며, D-Infinity 유수 방향 알고리즘의 산출물입니다.
- Contributing Area Grid [raster]
각 셀의 고유 기여에 해당 셀로 흘러드는 오르막 이웃들의 기여를 더한 각 셀 별 집수 영역 값을 제공하는 그리드입니다. 이 그리드는 D-Infinity 기여 영역 알고리즘을 사용해서 얻을 수 있습니다.
- Slope Exponent [number]
경사-면적 그리드를 생성하는 (Sm)(An) 공식에 사용되는 경사 지수(exponent) m 파라미터입니다.
기본값: 2
- Area Exponent [number]
경사-면적 그리드를 생성하는 (Sm)(An) 공식에 사용되는 면적 지수(exponent) n 파라미터입니다.
기본값: 1
산출물¶
- Slope Area Grid [raster]
경사 그리드, 특정 집수 영역 그리드, 경사 지수 파라미터 m, 그리고 면적 지수 파라미터 n 으로 계산한 경사-면적 값 (Sm) (An) 을 담고 있는 그리드입니다.
콘솔 사용법¶
processing.runalg('taudem:slopeareacombination', slope_grid, area_grid, slope_exponent, area_exponent, slope_area_grid)
참고¶
한계값으로 하천 정의¶
설명¶
어떤 그리드를 입력받아 입력값이 한계값 이상인 셀들을 식별하는 지표(indicator) 그리드 (1, 0)을 산출합니다. 일반적으로 집적원(accumulated source) 영역 그리드를 입력받아 하천 래스터 그리드를 생성합니다. 부가적인 입력 마스크 그리드를 사용하는 경우, 마스크 값이 0 이상인 셀 영역만 평가하도록 제한합니다. D-infinity 기여 영역 그리드(*sca)를 마스크 그리드로 사용하는 경우, 경계 오염 마스크로 기능합니다. 한계값 논리는 다음과 같습니다:
src = ((ssa >= thresh) & (mask >= s0)) ? 1:0
파라미터¶
- Accumulated Stream Source Grid [raster]
이 그리드는 일반적으로 유역(watershed)의 일부 특성 또는 특성들의 조합을 집적합니다. 정확한 특성(들)은 사용되는 하천망 래스터 알고리즘에 따라 달라집니다. 이 그리드는 산출되는 하천망이 연속하도록 그리드 셀 값이 내리막 D8 유수 방향을 따라 점증하는 속성을 보유해야 합니다. 집적 알고리즘이 이 그리드를 주로 산출하긴 해도, 최장 오르막 함수 같은 다른 알고리즘도 이런 그리드를 생성할 수 있습니다.
- Threshold [number]
어느 셀을 하천 셀로 간주해야 하는지 여부를 결정하기 위해 이 파라미터와 집적 하천원(Accumulated Stream Source) 그리드(*ssa) 내부의 값을 비교합니다. *ssa 값이 이 한계값 이상인 그리드 셀을 하천으로 식별합니다.
기본값: 100
- Mask Grid [raster]
Optional.
이 부가적인 입력물은 관심 영역을 마스크하기 위한 그리드로, 산출물은 이 그리드에서 0 이상의 값을 가진 영역만 커버합니다. 일반적으로 D-Infinity 기여 영역 그리드를 마스크로 사용해서 하천망을 D-Infinity 기여 영역 안으로 제한해서 묘사합니다. 이는 경계 오염 마스크의 기능과 동일합니다.
산출물¶
- Stream Raster Grid [raster]
이 입력물은 하천인 경우 그리드 셀 값을 1, 나머지 셀 값을 0으로 해서 하천의 위치를 나타내는 지표(indicator) 그리드 (1, 0)입니다.
콘솔 사용법¶
processing.runalg('taudem:streamdefinitionbythreshold', -ssa, -thresh, -mask, -src)
참고¶
하천 낙차 분석¶
설명¶
입력 집적 하천원 그리드(*ssa)에 (입력 파라미터로 결정된) 일련의 한계값을 적용해서 그 결과를 하천 낙차 통계 테이블 *drp.txt 파일로 산출합니다. 이 알고리즘은 하천을 묘사하는 데 사용되는 지형학적으로 객관적인 한계값을 결정하기 위해 개발됐습니다. 낙차 분석은 한계값들의 범위에서 하천망을 평가하고, 산출되는 스트랄러 하천의 일관적인 낙차 속성을 조사해서 올바른 한계값을 자동적으로 선택하려 합니다. 근본적으로 낙차 분석은 1차 하천의 평균 하천 낙차가 고차 하천의 평균 하천 낙차와 통계적으로 다른지를 T-검증을 사용해서 판정합니다. T-검증으로 두드러진 차이가 드러나는 경우 하천망이 이 “법칙”을 따르지 않는 것이기 때문에 더 큰 한계값을 선택해야 합니다. T-검증으로 두드러진 차이가 드러나지 않는 최소 한계값을 선택하면 지형학적으로 일관적인 하천 낙차 “법칙”을 따르는 최고 해상도의 하천망을 산출하며, 또한 이 최소 한계값은 DEM에서 하천의 “객관적인” 또는 자동적인 매핑을 위해 선택되는 한계값이기도 합니다. 하천망 래스터를 구축하는 데 이 알고리즘을 사용할 수 있습니다. 이때 집적 하천망 그리드에서 집적된 정확한 유역 특성(들)은 하천망 래스터 구축에 사용되는 메소드에 따라 달라집니다.
Broscoe (1959)에서 일관적인 하천 낙차 “법칙”을 식별했습니다. 하천 묘사 한계값을 결정하기 위해 이 법칙을 사용하는 데 쓰이는 과학에 대해 알고 싶다면, Tarboton et al. (1991, 1992), Tarboton & Ames (2001)을 참조하십시오.
파라미터¶
- D8 Contributing Area Grid [raster]
D8 알고리즘을 사용해서 계산한 각 셀 별 기여 영역 값 그리드입니다. 어느 셀의 기여 영역 값은 해당 셀 고유의 기여와 해당 셀로 흘러드는 모든 오르막 이웃에서 나오는 기여를 더한 값으로, 셀 개수로 측정됩니다. 일반적으로 이 그리드는 D8 기여 영역 알고리즘의 산출물입니다. 하천 낙차 테이블에 보고된 배수 농도를 평가하는 데 이 그리드를 사용합니다.
- D8 Flow Direction Grid [raster]
각 셀 별로 그 인접 또는 대각선 방향의 이웃 8개 셀 가운데 가장 심한 내리막 경사를 가지고 있는 셀의 방향으로 정의한 D8 유수 방향 그리드입니다. 이 그리드는 D8 유수 방향 알고리즘의 산출물로 얻을 수 있습니다.
- Pit Filled Elevation Grid [raster]
표고 값을 보유한 그리드를 지정합니다. 일반적으로 이 그리드는 구덩이 제거 도구의 산출물로, 구덩이를 제거한 표고입니다.
- Accumulated Stream Source Grid [raster]
이 그리드는 내리막 D8 유수 방향을 따라 점증해야만 합니다. 하천의 시작을 결정하기 위해 이 그리드와 일련의 한계값을 비교합니다. 이 그리드는 D8 기여 영역 알고리즘으로, 또는 D8 유수 경로 극값 알고리즘의 최대값 옵션을 사용해서 유역의 일부 특성들 또는 특성들의 조합을 집적해서 생성하는 경우가 많습니다. 정확한 방법은 사용하는 알고리즘에 따라 달라집니다.
- Outlets Shapefile [vector: point]
낙차 분석 수행 대상인 배출 상류를 정의하는 포인트 shapefile입니다.
- Minimum Threshold [number]
이 파라미터는 낙차 분석으로 가능한 한계값을 탐색하는 범위의 가장 낮은 값입니다. 이 기술은 t-통계의 절대값이 2 미만인 범위에 있는 최소 한계값을 찾습니다. 낙차 분석에 쓰이는 과학에 대해 알고 싶다면 Tarboton et al. (1991, 1992), Tarboton & Ames (2001)을 참조하십시오.
기본값: 5
- Maximum Threshold [number]
이 파라미터는 낙차 분석으로 가능한 한계값을 탐색하는 범위의 가장 높은 값입니다. 이 기술은 t-통계의 절대값이 2 미만인 범위에 있는 최대 한계값을 찾습니다. 낙차 분석에 쓰이는 과학에 대해 알고 싶다면 Tarboton et al. (1991, 1992), Tarboton & Ames (2001)을 참조하십시오.
기본값: 500
- Number of Threshold Values [number]
이 파라미터는 낙차 분석으로 가능한 한계값을 탐색하는 범위를 나누기 위한 단계의 개수입니다. 이 기술은 t-통계의 절대값이 2 미만인 범위에 있는 최소 한계값을 찾습니다. 낙차 분석에 쓰이는 과학에 대해 알고 싶다면 Tarboton et al. (1991, 1992), Tarboton & Ames (2001)을 참조하십시오.
기본값: 10
- Spacing for Threshold Values [selection]
이 파라미터는 낙차 분석으로 가능한 한계값을 탐색할 때 로그 간격(logarithmic spacing)을 사용해야 할지 또는 선형 간격을 사용해야 할지를 나타내는 파라미터입니다.
옵션:
0 – 로그
1 – 선형
기본값: 0
산출물¶
- D-Infinity Drop to Stream Grid [file]
다음 헤더 줄을 보유한 쉼표 구분 텍스트 파일입니다:
- ::
- Threshold,DrainDen,NoFirstOrd,NoHighOrd,MeanDFirstOrd,MeanDHighOrd,StdDevFirstOrd,StdDevHighOrd,T
이 파일은 조사한 각 한계값 별 데이터 한 줄, 그 다음으로 최적 한계값을 나타내는 요약 줄을 담고 있습니다. 이 기술은 t-통계의 절대값이 2 미만인 범위에 있는 최소 한계값을 찾습니다. 낙차 분석에 쓰이는 과학에 대해 알고 싶다면 Tarboton et al. (1991, 1992), Tarboton & Ames (2001)을 참조하십시오.
콘솔 사용법¶
processing.runalg('taudem:streamdropanalysis', d8_contrib_area_grid, d8_flow_dir_grid, pit_filled_grid, accum_stream_source_grid, outlets_shape, min_treshold, max_threshold, treshold_num, step_type, drop_analysis_file)
참고¶
- Broscoe, A. J., (1959), “Quantitative analysis of longitudinal stream profiles of small watersheds”, Office of Naval Research, Project NR 389-042, Technical Report No. 18, Department of Geology, Columbia University, New York.
- Tarboton, D. G., R. L. Bras and I. Rodriguez-Iturbe, (1991), “On the Extraction of Channel Networks from Digital Elevation Data”, Hydrologic Processes, 5(1): 81-100.
- Tarboton, D. G., R. L. Bras and I. Rodriguez-Iturbe, (1992), “A Physical Basis for Drainage Density”, Geomorphology, 5(1/2): 59-76.
- Tarboton, D. G. and D. P. Ames, (2001), “Advances in the mapping of flow networks from digital elevation data”, World Water and Environmental Resources Congress, Orlando, Florida, May 20-24, ASCE, http://www.engineering.usu.edu/dtarb/asce2001.pdf.
하천 범위 및 유역¶
설명¶
이 알고리즘은 하천 래스터 그리드로부터 벡터망 및 shapefile을 생성합니다. 하천 래스터를 따라 유수 경로를 연결하기 위해 유수 방향 그리드를 사용합니다. 각 하천 구간의 스트랄러 순서를 계산합니다. 각 하천 구간(직선 유역)으로 흘러드는 소유역(subwatershed)도 묘사하며, 하천 직선 유역(Stream Reach) shapefile 안에 있는 WSNO(유역 번호)에 대응하는 값 식별자를 라벨로 표시합니다.
이 알고리즘은 하천망을 스트랄러 순서 체계에 따라 순서를 부여합니다. 다른 어떤 하천도 흘러들지 않는 하천이 1차(order 1)입니다. 서로 다른 순서의 하천 구간 2개가 합쳐질 경우 그 하류 구간의 순서는 해당 하천 구간들 가운데 더 높은 순서를 따릅니다. 동일한 순서의 하천 구간 2개가 합쳐질 경우 그 하류 구간의 순서는 1만큼 증가합니다. 하천 구간 2개 이상이 합쳐질 경우 유수 경로 순서는 가장 높은 구간 순서의 최대값, 또는 두 번째로 높은 구간 순서에 1을 더한 값으로 계산됩니다. 어느 한 지점에서 2개 이상의 구간이 합쳐지는 경우 스트랄러 순서는 일반적으로 이렇게 정의됩니다. 하천망의 지형적 연결성은 하천망 트리 파일에 저장되며, 하천망을 따라 각 그리드 셀의 좌표 및 속성은 하천망 좌표 파일에 저장됩니다.
이 알고리즘은 하천 래스터 그리드를 하천망의 소스로 사용하며, 유수 방향 그리드를 하천망 내부의 연결점을 추적하는 데 사용합니다. 하천망 좌표 파일 안에 표고 및 기여 영역 속성을 결정하는 데 표고 그리드 및 기여 영역 그리드를 사용합니다. 관측 지점의 상류 및 하류 유역을 쉽게 표현하기 위해 하천 구간을 논리적으로 나누는 데 배출 shapefile에 있는 포인트들을 사용합니다. 이 알고리즘은 배출 shapefile에 있는 id 속성 필드를 하천망 트리 파일의 식별자로 사용합니다. 그 다음 하천망 트리 및 좌표 파일에 있는 벡터망 표현 텍스트를 shapefile로 변환합니다. 더 많은 속성들도 평가합니다. 이 알고리즘은 산출 유역 그리드에 하천망으로 흘러드는 전체 영역을 단일 값으로 표현해서 단일 유역을 묘사하는 옵션도 가지고 있습니다.
파라미터¶
- Pit Filled Elevation Grid [raster]
표고 값을 보유한 그리드를 지정합니다. 일반적으로 이 그리드는 구덩이 제거 도구의 산출물로, 구덩이를 제거한 표고입니다.
- D8 Flow Direction Grid [raster]
각 셀 별로 그 인접 또는 대각선 방향의 이웃 8개 셀 가운데 가장 심한 내리막 경사를 가지고 있는 셀의 방향으로 정의한 D8 유수 방향 그리드입니다. 이 그리드는 D8 유수 방향 알고리즘의 산출물로 얻을 수 있습니다.
- D8 Drainage Area [raster]
D8 알고리즘을 사용해서 해당 셀 고유의 기여와 해당 셀로 흘러드는 모든 오르막 이웃에서 나오는 기여를 더한 값을 각 셀 별 그리드 셀 개수(또는 가중치의 총체)라는 측면에서 계산한 기여 영역 값을 제공하는 그리드입니다. 일반적으로 이 그리드는 D8 기여 영역 알고리즘의 산출물입니다. 하천망 좌표 파일의 기여 영역 속성을 결정하기 위해 이 그리드를 사용합니다.
- Stream Raster Grid [raster]
하천인 경우 그리드 셀 값을 1, 나머지 셀 값을 0으로 해서 하천을 나타내는 지표(indicator) 그리드입니다. 하천망 분석 도구 모음 가운데 하나로 이런 그리드를 생성할 수 있습니다. 하천 래스터 그리드를 하천망의 소스로 사용합니다.
- Outlets Shapefile as Network Nodes [vector: point]
Optional.
관심 지점을 정의하는 포인트 shapefile입니다. 이 파일을 사용하는 경우, 알고리즘이 이 배출 지점 상류의 하천망만을 묘사할 것입니다. 또한, 관측 지점의 상류 및 하류 유역을 쉽게 표현하기 위해 하천 구간을 논리적으로 나누는 데 배출 shapefile에 있는 포인트들을 사용합니다. 이 알고리즘은 배출 shapefile이 정수형 속성 필드 id 를 가지고 있을 것을 요구합니다. 하천망 트리 파일이 id 값을 식별자로 사용하기 때문입니다.
- Delineate Single Watershed [boolean]
산출 유역 그리드에 하천망으로 흘러드는 전체 영역을 단일 값으로 표현해서 단일 유역을 묘사하도록 하는 옵션입니다. 이 옵션을 거짓으로 설정하면 각 하천 구간 별로 개별 유역을 묘사합니다. 기본값은 False (개별 유역)입니다.
기본값: False
산출물¶
- Stream Order Grid [raster]
하천 순서 그리드는 스트랄러 순서 체계에 따라 순서를 부여한 하천 셀 값을 담고 있습니다. 스트랄러 순서 체계는 다른 어떤 하천도 흘러들지 않는 하천 구간을 1차(order 1)로 정의합니다. 서로 다른 순서의 하천 구간 2개가 합쳐질 경우 그 하류 구간의 순서는 해당 하천 구간들 가운데 더 높은 순서를 따릅니다. 동일한 순서의 하천 구간 2개가 합쳐질 경우 그 하류 구간의 순서는 1만큼 증가합니다. 하천 구간 2개 이상이 합쳐질 경우 유수 경로 순서는 가장 높은 구간 순서의 최대값, 또는 두 번째로 높은 구간 순서에 1을 더한 값으로 계산됩니다. 어느 한 지점에서 2개 이상의 구간이 합쳐지는 경우 스트랄러 순서는 일반적으로 이렇게 정의됩니다.
- Watershed Grid [raster]
이 산출 그리드는 각 구간 유역을 유일 ID 번호로 식별하거나, 또는 단일 유역 옵션을 참으로 설정한 경우 하천망으로 흘러드는 전체 영역을 단일 ID로 식별합니다.
- Stream Reach Shapefile [vector]
이 산출물은 하천망 내부의 링크를 지정하는 폴리라인 shapefile입니다. 속성 테이블에 있는 열들은 다음과 같습니다:
LINKNO – 링크 번호입니다. 각 링크(합류점 사이의 수로 구간)에 관련된 유일한 숫자입니다. 랜덤한 숫자로, 공간 처리 횟수에 따라 달라집니다.
DSLINKNO – 하류 링크의 링크 번호입니다. -1 은 하류 링크가 존재하지 않는다는 뜻입니다.
USLINKNO1 – 첫 번째 상류 링크의 링크 번호입니다. (예를 들어 소스 링크의 경우처럼, -1 은 상류 링크가 존재하지 않는다는 뜻입니다.)
USLINKNO2 – 두 번째 상류 링크의 링크 번호입니다. (예를 들어 소스 링크의 경우 또는 하천 구간이 논리적으로 나누어지지만 하천망은 두 갈래로 나뉘지 않는 내부 관측 지점의 경우처럼, -1 은 두 번째 상류 링크가 존재하지 않는다는 뜻입니다.)
DSNODEID – 하천 구간의 하류 종단점에 있는 노드의 노드 식별자입니다. 이 식별자는 배출 shapefile에서 노드를 지정하기 위해 사용되는 id 속성 필드와 대응합니다.
Order – 스트랄러 하천 순서
Length – 링크의 길이입니다. 기저 DEM 그리드의 수평 맵 단위를 사용합니다.
Magnitude – 링크의 슈리브 규모(Shreve Magnitude)입니다. 상류 소스의 총 개수를 뜻합니다.
DS_Cont_Ar – 링크의 하류 종단점에 있는 배수 영역입니다. 일반적으로 이 영역은 하류 종단점의 상류 그리드 셀 1개인데, 왜냐하면 하류 종단점 그리드 셀에 있는 배수 영역이 하천이 합류하는 영역을 포함하기 때문입니다.
Drop – 링크의 시작점에서 종단점 사이의 표고 낙차입니다.
Slope – 링크의 평균 경사입니다. (낙차를 길이로 나눈 값입니다.)
Straight_L – 링크의 시작점에서 종단점까지의 직선 거리입니다.
US_Cont_Ar – 링크의 상류 종단점에 있는 배수 영역입니다.
WSNO – 유역 번호입니다. 링크로 직접 흘러드는 유역의 식별 번호를 지정하는 *w.shp 파일 및 *w 그리드 파일을 교차 참조합니다.
DOUT_END – 링크의 하류 종단점에서 (예를 들어 하천망에서 가장 하류 지점인) 최종 배출 지점까지의 거리입니다.
DOUT_START – 링크의 상류 종단점에서 최종 배출 지점까지의 거리입니다.
DOUT_MID – 링크의 중앙 지점에서 최종 배출 지점까지의 거리입니다.
- Network Connectivity Tree [file]
이 산출물은 하천망 트리 파일에 저장된 하천망의 지형적 연결성을 상세히 설명하는 텍스트 파일입니다. 다음과 같은 열들을 담고 있습니다:
링크 번호 (랜덤 – 공간 처리 횟수에 따라 달라집니다.)
하천망 좌표 파일(*coord.dat)에서 시작 지점 번호 (인덱스가 0부터 시작합니다.)
하천망 좌표 파일(*coord.dat)에서 종단 지점 번호 (인덱스가 0부터 시작합니다.)
다음 (하류) 링크 번호 (링크 번호를 가리킵니다. 예를 들어 최종 링크의 경우처럼, -1 은 상류 링크가 존재하지 않는다는 뜻입니다.)
첫 번째 이전 (상류) 링크 번호. (링크 번호를 가리킵니다. -1 은 상류 링크가 존재하지 않는다는 뜻입니다.)
두 번째 이전 (상류) 링크 번호. (링크 번호를 가리킵니다. -1 은 상류 링크가 존재하지 않는다는 뜻입니다. 이전 링크 하나만 -1 인 경우, 하천 구간이 논리적으로 나누어지지만 하천망은 두 갈래로 나뉘지 않는다는 뜻입니다.)
링크의 스트랄러 순서
링크의 하류 종단점에 있는 관측 지점의 식별자 (-1 은 하류 종단점이 관측 지점이 아니라는 뜻입니다.)
링크의 하천망 규모 (슈리브의 정의에 따라 상류 소스의 개수로 계산합니다.)
- Network Coordinates [file]
이 산출물은 하천망을 따라 각 포인트의 좌표 및 속성을 담고 있는 텍스트 파일입니다. 다음과 같은 열들을 담고 있습니다:
X 좌표
Y 좌표
최종 링크의 하류 종단점까지의 수로 거리
표고
기여 영역
콘솔 사용법¶
processing.runalg('taudem:streamreachandwatershed', -fel, -p, -ad8, -src, -o, -sw, -ord, -w, -net, -tree, -coord)
참고¶
측정기 유역¶
설명¶
측정기 유역(Gage Watershed) 그리드를 계산합니다. 각 그리드 셀에 해당 셀이 다른 어떤 측정기도 통과하지 않고 직접 흘러드는 측정기의 (id 열에서 나온) 식별자 라벨을 부여합니다.
파라미터¶
- D8 Flow Directions Grid [raster]
각 셀 별로 그 인접 또는 대각선 방향의 이웃 8개 셀 가운데 가장 심한 내리막 경사를 가지고 있는 셀의 방향으로 정의한 D8 유수 방향 그리드입니다. 이 그리드는 D8 유수 방향 알고리즘의 산출물로 얻을 수 있습니다.
- Gages Shapefile [vector: point]
유역을 묘사할 측정기들을 정의하는 포인트 shapefile입니다. 이 shapefile은 id 열을 가지고 있어야 합니다. 이 shapefile에서 각 포인트로 직접 흘러드는 그리드 셀들에 이 id 라벨을 부여할 것입니다.
산출물¶
- Gage Watershed Grid [raster]
각 측정기 유역을 식별하는 그리드입니다. 각 그리드 셀에 해당 셀이 다른 어떤 측정기도 통과하지 않고 직접 흘러드는 측정기의 (id 열에서 나온) 식별자 라벨을 부여합니다.
- Downstream Identifiers File [file]
유역의 내리막 연결성을 지정하는 텍스트 파일입니다.
콘솔 사용법¶
processing.runalg('taudem:gagewatershed2', -p, -o, -gw, -id)