고해상도 LiDAR DEM에서 정확한 배수 네트워크 (및 집수구)를 생성하는 방법론?


28

내가이 문제를 우연히 발견 한 것은 이번이 처음이 아닙니다. 전체 배수 (1m 셀) LiDAR 데이터에서 올바른 배수 네트워크 모델과 결과물을 생성 할 수없는 것 같습니다.

LiDAR 데이터 세트를 일반화하고 정수 DEM 및 채우기 싱크로 변환하면 모든 것이 좋으며 매우 일반화 된 모델로 보이는 것을 쉽게 만들 수 있습니다. 그러나 대규모지도에 대한 자세한 사이트 모델을 만들고 싶습니다. 이곳에서 문제가 발생합니다.

나는 대부분의 문제가 더 평평한 지역에서 발생한다는 것을 지적해야한다.

배수 네트워크가 지형을 정확하게 따르기를 원하지만 정수 DEM 입력 에서 배수 네트워크를 만들 때 발생하는 스트림은 매우 일반적이며 자주 연결되지 않아야하는 영역에서 "분리"됩니다. 개울은 지형의 자연 능선을 밀접하게 따르지 않습니다. "고아"또는 "아무데도"세그먼트가 많이 있습니다. 내가 사용하는 경우 부동 소수점 DEM 입력 , 결과 배수 네트워크는 상세하고 정확하지만, 아주, 연결 클러스터와 고아 스트림을 "산재"입니다.

내 문제가 데이터 준비 어딘가에 있다고 생각합니다. 정수 대 부동 소수점 래스터 DEM 입력, 싱크가 올바르게 채워지는 등. 또는 어떻게하면 먼저 "수 문학적으로 올바른"입력 DEM을 생성하기 위해 표면 데이터를 처리해야합니까?

고해상도 LiDAR을 사용하여 연속 배수 네트워크 및 집수구를 작성하는 올바른 방법을 설명 할 수 있습니까?

정수 DEM 입력에서 모델을 만드는 데 더 성공했습니다. 그러나 자세한 대규모 분석에는 적합하지 않습니다.

첫 번째 첨부 이미지는 정수 DEM 입력에서 생성 된 모델입니다. 몇 가지 명백한 문제 영역이 있습니다. 실제로 주 배수 채널로 보이는 곳에 스트림이 있습니다. 매우 일반적인 버전의 스트림을 추가했습니다. 여기에 이미지 설명을 입력하십시오

편집 : 이미 언급했듯이 정수 DEM 입력에서 모델을 만드는 데 더 많은 성공을 거두었습니다. 다음 화면 캡처는 그 이유를 보여줍니다. 정수 DEM 입력은 위에서 볼 수있는 바와 같이 많은 문제를 가지고 있지만, 지형 특성에 맞지 않더라도 여전히 연결이 끊긴 배수 네트워크를 생성합니다. 부동 소수점 DEM 입력을 사용하여 바로 아래 이미지에서 볼 수 있듯이 작은 고아 세그먼트로 가득한 매우 분리되고 클러스터 된 네트워크가 생성됩니다.

부동 소수점 DEM에서 생성 된 유량 누적 래스터 여기에 이미지 설명을 입력하십시오

정수 DEM으로 생성 된 Flow Accumulation 래스터 여기에 이미지 설명을 입력하십시오

제가 공제 할 수있는 한, 두 방법 모두 극적으로 다른 결과를 산출하며, 두 방법 모두 상세 모델에는 사용할 수 없습니다.

편집 : 나는이 게시물을 더 길고 길게 만드는 것에 대해 사과드립니다 (아마도 영어로 명확하게 표현하지 않을 것입니다) 부동 소수점 DEM을 입력에 사용하는 문제를 추가로 설명하기 위해 결과 Stream Link 출력과 결과 유역을 첨부합니다. 내가 기대하는 것은 연속 스트림 네트워크와 모든 영역이 서로 흘러 들어가는 분지로 덮여 있습니다.

부동 소수점 입력 DEM에서 생성 된 스트림 링크 : 여기에 이미지 설명을 입력하십시오

부동 소수점 입력 DEM에서 생성 된 유역 유역 : 여기에 이미지 설명을 입력하십시오

정수 DEM 입력을 사용하여 유역의 전체 흐름 방향이 변경되는 예 (가까운 영역, 동일한 데이터)는 다음과 같습니다. 빨간색 화살표는 모델의 흐름 방향이고 파란색 화살표는 실제 흐름의 방향을 나타냅니다. . (파란색 선-실제 스트림, 빨간색 네트워크는 LiDAR 파생 스트림 네트워크 Strahler 주문) 여기에 이미지 설명을 입력하십시오

데이터 링크 : https://www.yousendit.com/download/MEtSOGNVNXZvQnRFQlE9PQ(2011 년 5 월 13 일 만료)



언덕 음영은 어디에서 왔습니까? (검은 색) 유량 누적 결과는 언덕이있는 고도에서 도출 된 것으로 보이지 않습니다. 아마도 당신은 우리에게 같은 맵을 보여줄 수 있지만 유량 누적 값을 얻는 데 사용되는 그리드의 언덕 음영 표현이 있습니다.
whuber

권리. 나는 그것을 언급 했어야했다. 언덕 모양은 동일한 그리드에서 파생됩니다. (그리고 검은 색 스트림 네트워크는 스트림 링크 래스터에서 파생 된 스트림 순서 (Strahler)입니다) 스트림 위치 (파란색)를 제외한이 맵의 모든 항목은 동일한 그리드에서 생성됩니다.
Jakub Sisak GeoGraphics

2
복잡한 문제를 좁히기위한 나의 충고는 간단한 테스트 사례를 사용하는 것입니다. 원시 래스터 소스에서 작은 조각을 자르고 원하는대로 단계를 시도하십시오 (예 : 플로트로 유지). 항상 싱크대를 채우십시오. 각 단계의 출력을주의 깊게 검사하여 "올바로"보이는지 확인하십시오.
Mike T

2
야콥 나는 똑같은 문제를 겪는다. 당신은 혼자가 아닙니다! 내가 이전에받은 대답은 배수 네트워크를 만들기 위해 LiDAR 데이터를 사용하지 않는 것이 었습니다.
Jacques Tardie

답변:


11

GRASS GIS 분석을 사용해 보셨습니까? GRASS 알고리즘이 수 문학 분석에 매우 적합하다는 기대감을 가지고 있습니다. 예를 들어, 해상도 5x5m의 DTM에서 배수 네트워크와 같은 것을 생성하고 싶습니다. ArcMap (ArcHydro Tools 포함)의 도구를 비교했으며 첫 번째 그림 (빨간 선)에서 결과를 볼 수 있습니다. 그런 다음 GRASS GIS 기능 'r.stream.extract'를 사용하려고했는데 그림 2 (빨간 선)에 결과가 표시되었습니다. 두 배수 라인 모두 3 헥타르의 카테 멘트 영역으로 생성됩니다.

그것은 실제로 다르며 실제 스트림과 비교하여 꽤 정확성이 있습니다 (그림 3, 실제 스트림은 파란색입니다). GRASS GIS에는 많은 수 문학적 도구가 있습니다.

ArcMap을 이용한 배수관] GRASS GIS를 사용한 배수 라인 GRASS GIS 배수 라인과 실제 스트림 비교


1
매우 흥미로운! ESRI 도구를 사용하는 것과 동일한 오류가 발생할 수 있습니다. 이것은 ESRI 알고리즘이 고해상도 데이터를 처리 할 수 ​​없다고 믿게합니다. 이것은 거의 질문에 대한 답변입니다. 영상 주셔서 감사합니다-놀라운! 유역 / 배수 분석을 위해 GRASS 도구를 사용한 경험이 없습니다. 기본 "방법"자습서를 알려 주시면 감사하겠습니다.
Jakub Sisak GeoGraphics

1
이것이 훌륭하다고 말하고 싶었습니다! Lidar 데이터 세트에서 동료와 함께 예비 테스트를 실행하면 초기 결과가 매우 유망 해 보입니다. 많은 기능과 매개 변수 및 일부지도를 추가하는 기능도 훌륭합니다. 결과는 실제 스트림과 일치합니다. 또한 80 년대 중반 이후로 ESRI 알고리즘이 얼마나 오래된 지 알아 내십시오. 많은 설명이 있습니다. 고맙습니다!
Jakub Sisak GeoGraphics

도와 주셔서 감사합니다! 나는 많은 수 문학적 분석과 그것이 제공하는 아주 좋은 결과를 위해 GRASS GIS를 좋아합니다. 말했듯이 ESRI는 정말 구식입니다. 나는 오래된 것이 얼마나 오래된지도 모른다. 더 많은 수 문학적 분석을 시도하려면 grasswiki.osgeo.org/wiki/Hydrological_Sciencesgrass.osgeo.org/grass70/manuals/topic_hydrology.html 페이지를 확인하십시오 .
david_p

GRASS GIS는 내가 생각하는 단 하나의 단점이 있으며 이는 레이어의 특별한 기본 환경입니다. GRASS를 모르는 사람들에게는 약간 성가신 일입니다. 그러나 잠시 후, 당신은 그것에 사용했습니다.
david_p

8

수 문학적으로 올바른 고도 모델 (배수 강제, ANUDEM 이라고도 함)을 생성하는 것과 관련하여 제 지식으로는 최고의 품종입니다. 캐나다 국가 별 고도 데이터 세트 (CDED, 아이런 식으로 정수 미터로 저장 됨) 를 생성하는 데 사용되는 프로그램 입니다. 또한 ArcGIS 의 TopoToRaster 도구는 후드 아래에서 Anudem을 사용합니다 (현재 개정 또는 3 배).

USGS는 AverStar의 미국 모델 Delta3D에 대해 다른 프로그램을 사용했지만 10 년 전에 문의했을 때이 프로그램은 맞춤형 프로그램이었으며 기성품으로 사용할 수 없었습니다 (100k 정도는 우리의 요구에 맞게 조정했지만) ).

배수 강제 입면 모델을 생성하기위한 다른 도구는 알지 못하지만 그에 대해 듣고 싶습니다.


나는 실제로 이것을 시도했지만 toold는 많이 충돌합니다. LiDAR 파생 윤곽선 (2K x 2K 하위 집합)을 사용한 다음 작은 중요하지 않은 윤곽선을 제거하여 표면을 더 단순하게 만들고 TopoToRaster를 시도했지만 계속 죽어 가고 있습니다. (형상 폴리 라인 오류에서 너무 많은 점) 대신 점 입면을 시도해야합니까?
Jakub Sisak GeoGraphics

그리고 CDED에 관해서는 정수 반올림과 그로 인한 "terrace anomaly"문제와 관련된 모든 종류의 문제 (아직 해결되지 않음)가 있습니다.
Jakub Sisak GeoGraphics

LiDAR 점을 점 (점) 입력으로 사용하여 TopoToRaster 도구를 사용하여 "수 문학적으로 올바른"표면을 성공적으로 만들 수있었습니다. 출력 셀 크기가 서로 다른 2 개의 표면을 생성했습니다. 2 및 4. 결과 흐름 누적 래스터도 동일한 문제를 겪습니다. ArcGIS에서이 작업을 수행 할 수없는 것으로 의심되기 시작했습니다. 또한 TopoToRaster를 실행하는 데 시간이 오래 걸린다는 점도 지적하고 싶습니다.
Jakub Sisak GeoGraphics

5

대학으로 돌아가서 나는 이것을 아주 잘한 프로젝트에서 일했다. 나는 수 문학자가 아니고 프로젝트를 마치지 않았지만 (점검), 이것을 확인하고 싶을 수도 있습니다.

TauDEM 5.0

내가 기억하는 것에서 그것은 꽤 잘 작동했습니다. 무료 도구이며 필요한 것일 수도 있습니다.

편집 : 질문을보다 신중하게 읽은 후에는 이것이 정확히 필요한 도구라고 생각합니다. 설명에 따라 단절이 없으며 모든 흐름이 다운 스트림으로 계속됩니다 (예 : 고아 스트림 없음). 대부분의 DEM은 N, E, S, W 및 NE, SE, SW, NW의 가능한 방향으로 만 흐름 방향을 계산합니다. 이것은 부 자연스러운 흐름을 초래합니다. TauDEM의 무게 방향은 360 도입니다. 더 자연스러운 흐름을 가질 것이고 더 정확한 흐름을 가정합니다.

또한 코어가 여러 개인 경우이를 사용합니다. TauDEM은 고해상도 LiDAR을 사용하여 필요한 것을 신속하게 처리해야합니다.


2
나는 이것을 두 번째로 할 것이다! D8 흐름 방향은 TauDEM이 D 무한 흐름 방향을 사용할 수있는 고해상도 데이터에서 바람직하지 않은 결과를 생성합니다. 또한 수 문학적 흐름 모델의 의도를 명심하십시오. 클수록 항상 더 좋은 것은 아닙니다 (현명한 해결책). 초 고해상도 DEM은 모델보다 더 큰 문제입니다. Lidar 파생 DEM은 본질적으로 흐름 모델에 사용하도록 의도되지 않은 실제 "소음"을 가지고 있습니다. DEM을 축소하는 것이 좋습니다.
Jeffrey Evans

SAGA GIS 소프트웨어도 확인하십시오. 정보 (즉, x, y, z)가 각기 다른 유량 누적 방법 (D8, Dinfinity 등)에서 항상 동일하기 때문에 이것은 관련없는 문제가 아님을 지적하고 싶습니다. . SAGA GIS에서 발견 된 병렬 처리는 라이더 데이터의 빠른 처리를 허용합니다. 나는이 방법을 상당히 큰 계산에 사용했으며 잘 작동했습니다. 문제는 데이터를 올바르게 전처리하는 것입니다. 즉, 배수 구조 (구름대, 다리)를 태우고 채우면 유량 누적 계산이 이루어집니다!
reima

타우 뎀은 다중 프로세서 기능도 가지고있었습니다.
모르는 경우 GIS

4

기부 해 주셔서 감사합니다. 전 해상도의 LiDAR 표면이 이러한 유형의 분석에 적합하지 않다는 결론을 내 렸습니다.


이 기사 를 사용하는 10 가지 이유Terrain Datasets 는 DEM 래스터 표면이 귀하의 경우에 사용하기에 잘못된 데이터 모델이라고 생각했습니다. 패싯이 실험에서 너무 많은 인공물을 생성함에 따라 고도 모델에 대한 TIN을 거부했습니다. 그러나 소스 데이터는 윤곽이며 Lidar와 같은 밀집된 스팟 높이 필드가 아닙니다.
matt wilkie 16:17에

4

정수 또는 부동 소수점 사용에 대한 질문에 특히 : 정수는 속도, 저장에 가장 적합하며 반올림 오류로 인해 어떤 종류의 드리프트를 피합니다. 그러나 정수를 사용할 때 Z (고도) 값에 미터를 사용하지 마십시오! 세로 단위를 센티미터 또는 밀리미터로 변경하거나 미터로 유지하고 동일한 효과를 갖는 값을 100 또는 1000으로 곱합니다. 이것이 가능하지 않으면 부동 소수점을 사용하십시오.

경사 및 종횡비 분석 및 기타 2 차 및 3 차 미분은 미터 기반 정수 고도의 조도에 특히 민감합니다. 실제로 나쁜 연습이지만 표준 연습이기도합니다.

참조 지형 분석 : 원리와 응용 프로그램 특정 섹션에서 (존 피터 윌슨 & 존 C. 호협) 2.7.2 고도 단위 및 수직 정밀도디지털 고도 모델의 지형 학적 특성 ( 조 우드 ), "정수 반올림"를 검색합니다. 이 두 문서는 모두 무겁습니다. ANUDEM 소프트웨어를 사용하여 호주 최초의 대륙 별 고도 모델 (2000 년경)을 구축하는 문서에서 문제에 대한 간결하고 이해하기 쉬운 설명을 통해 문제를 처음 알게 되었지만 지금은 찾을 수 없습니다.


1
고마워 매트. 좋은 물건. 나는 이것을 시도하고 다시보고합니다. 매우 흥미로운 정보가 많이 있습니다. 노력해 주셔서 감사합니다.
Jakub Sisak GeoGraphics

1
입력 래스터에 1000을 곱하면 이전과 동일한 결과를 얻습니다. 정수와 부동 소수점을 모두 시도했습니다. 결과 흐름 누적 래스터는 두 경우 모두 거의 동일합니다. TopoToRaster 기술을 사용하려고합니다.
Jakub Sisak GeoGraphics

3

이것이 도움이 될지 모르지만 나는 1cm LIDAR DEM을 위해 하이드로 네트워크에 다시 블로그 게시물을 썼습니다. 너겟 몇 개가 있겠지

http://www.thadwester.com/1/post/2011/03/hydrologic-networks.html


감사. 특히 평평한 지역에서 사용 가능한 배수 네트워크를 도출 할 수있는 연속 흐름 방향 래스터를 얻는 것이 나의 주요 문제인 것 같습니다. 흐름 방향 래스터를 생성하기 위해 ArcGIS에서 D8 방법을 사용하는 방법을 설명해 주시겠습니까?
Jakub Sisak GeoGraphics 2016 년

이것에 추가하십시오. 메모리 문제가있을 수 있다고 판단되면 플로우 업스트림을 추적 할 수없는 컷오프 매개 변수가있는 것처럼 기존 알고리즘이 무한하지 않다고 생각합니다.
Jakub Sisak GeoGraphics 2016 년

Arcgis에서 흐름 방향 래스터를 만들 수 있습니다. 당신이 능력이 없다면 당신을 위해 그것을 실행할 수 있습니다.
Thad

죄송합니다. 흐름 방향이 아니라 위의 주석에서 흐름 축적을 의미했습니다. 이것은이 질문에서 설명한 초기 문제입니다. 거짓말 방향이 낮은 곳에서 밀집된 Lidar 데이터를 실행할 때 흐름 방향 도구 도구는 유용한 결과를 얻지 못합니다. 실제로 부동 소수점 래스터를 사용하면 돌이킬 수없는 오류가 발생하지만 정수 래스터를 사용하면 데이터가 너무 일반화됩니다. 그대로 ArcGIS 도구 만 사용하여 LiDAR 데이터에서 정확한 배수 모델을 도출하는 것은 불가능합니다.
Jakub Sisak GeoGraphics

데이터 밀도가 어떻게 중요한지 이해하지 못합니다. 내가 사용하는 데이터는 1 cm ^ 2입니다. 더 조밀 한. 데이터를 다운로드하겠습니다. 시도하겠습니다.
Thad

1

여기에 대해 더 생각할 것이 있다고 생각했습니다. 유역 분지 묘사 프로세스가 작동하는지 여부를 묻습니다. 수동으로 편집 한 모델이 있으며 계속 잘못된 영역을 찾아오고 있습니다. ArcGIS 컴퓨터 생성 모델에 전혀 의존 할 수 없다고 생각합니다.

스크린 샷

당사 사이트를 사용함과 동시에 당사의 쿠키 정책개인정보 보호정책을 읽고 이해하였음을 인정하는 것으로 간주합니다.
Licensed under cc by-sa 3.0 with attribution required.