답변:
@jul은 "추세 표면을 계산 한 다음 초기 DEM에서 빼서 '추세를 제거한'것을 얻는 것이 맞지만이 경우에는"딥을 유지 "하기 위해 더 간단한 절차가 필요한 것 같습니다. "추세 표면"이 원래 DEM을 너무 가깝게 따르는 경우 잔차는 표면의 로컬 특성을 유지하지 않습니다. 따라서 피해야 할 기술 중에는 모든 로컬 (스플라인, 필터 및 특히 크릭)이 있으며 전체적으로 선호되는 기술이 있습니다.
간단하고 강력하며 직접적인 접근 방식은 강 근처의 비행기를 DEM에 맞추는 것입니다. (유클리드에 따르면) 비행기는 공간에서 3 개의 (비-선형) 점에 의해 결정되기 때문에 멋진 기술이나 무거운 계산은 필요하지 않습니다. 따라서 하천 끝에서 한 지점 (x1 ', y1', z1 ') = (좌표, 고도)를, 하류 끝에서 다른 지점 (x2', y2 ', z2')을 선택하십시오. (x0 ', y0', z0 ')을 선택하면 처음 두 점을 연결하는 선분에서 멀어집니다. (이 좌표는 곧 변경 될 것이므로 소수로 표시됩니다.)이 마지막 점은 그렇지 않습니다지면 위 또는 근처의 지점에 해당해야합니다! 실제로 초기 선택은 표고를 상류 및 하류의 평균 인 z0 '= (z1'+ z2 ') / 2로 설정하는 것이 좋습니다.
점 (x0 ', y0', z0 ')을 로컬 좌표계의 원점으로 채택하면 계산이 쉬워 집니다. 이 좌표에서 다른 두 점은
(x1,y1,z1) = (x1'-x0',y1'-y0',z1'-z0')
(x2,y2,z2) = (x2'-x0',y2'-y0',z2'-z0').
원래 좌표계의 (x ', y')에있는 임의의 위치는이 새 시스템의 좌표 (x, y) = (x'-x0 ', y'-y0')입니다. 원점 (0,0,0)을 통과하는 평면은 z = a * x + b * y 형식의 방정식을 가져야하므로 문제가 다음과 같이 줄어 듭니다.
점 (0,0,0), (x1, y1, z1) 및 (x2, y2, z2)를 통과하는 평면에 대해 z = a * x + b * y 형식의 방정식을 찾으십시오.
독특한 솔루션은 계산하는 것입니다
u = z1 y2 - z2 y1
v = x1 z2 - x2 z1
w = x1 y2 - x2 y1
어떤 용어로
a = u/w, b = v/w.
이 두 개의 숫자 a와 b를 발견하고 두 개의 원래 좌표 x0 '와 y0'을 회상 하여 양식의 래스터 계산
[DEM] - a * ([X'] - x0') - b * ([Y'] - y0')
DEM에서 "틸트"를 제거합니다 . 이 식에서 [X ']는 원래 좌표 의 x 좌표 그리드 를 나타내고 [Y']는 원래 y 좌표 그리드를 나타냅니다. 결과 DEM은 원래 선택한 세 점 각각에서 동일한 높이 (즉, z0 ')를 갖습니다. 다른 곳에서하는 일은 DEM 자체에 달려 있습니다!
(일부 독자들은이 접근법이 삼각법 또는 최소 제곱 기계에 대한 모든 참조를 피하는 방법을 높이 평가하기를 바랍니다. :-)
강물에 대한 경향 분석이 불만족 스러웠습니다. 중서부와 같은 지역에서는 작동하지만 서구에서는 단일 지형 또는 다항식 표면이 복잡한 지형 때문에 강 구배와 거의 비슷한 것으로 보이지 않습니다. 대신에 사용한 것은 가중 평균입니다. 여기서 커널 밀도를 사용하여 고지대 픽셀의 높이를 강의 가중 평균과 비교합니다 (Winter 2010 ArcUser의 기사 참조). 보다 "현지화 된"결과가 생성됩니다. 음수 값을 얻는 것이 여전히 가능하지만 그 값이 강 수로에 국한된 것으로 보이며 가파른 경사 지역에서 가장 일반적입니다. 더 최근에는 "flood height"를 매핑하기위한 경로 거리 접근법을 실험했지만이 방법은 그렇지 않습니다.
트렌드 표면을 계산 한 다음 초기 DEM에서 빼면 "추세가있는"것을 얻을 수 있습니다.
트렌드 표면을 계산하기 위해, 임의의 표면 단순화 기능 (필터링, 스플라인 스무딩, 크 래깅 등)이 사용될 수 있습니다. 이 선택은이 추세 표면의 모양에 따라 다릅니다.
다음은 arcGIS 10.0에서 크릭을 사용하는 방법에 대한 예 입니다.