[논문 리뷰] A class of second-order geometric quasilinear hyperbolic PDEs and their application in imaging science

Guozhi Dong, Michael Hintermüller외 1명|arXiv (Cornell University)|2019. 01. 01.

Numerical methods in inverse problems참고 문헌 31인용 수 1

한 줄 요약

이 논문은 영상 처리를 위한 두 번째 차수 기하학적 비선형 쌍곡형 편미분방정식의 클래스를 제안하며, 감쇠된 두 번째 차수 총변동성 및 평균 곡률 흐름에 초점을 맞춘다. 해의 존재성과 유일성을 확립하고, 단순한 초기 자료에 대해 해석적 해를 제공하며, 특히 노이즈 제거 및 재편성 오류 수정(예: dejittering) 응용에서 일阶 흐름보다 우수한 수렴성과 경계 유지 성능을 수치적 비교를 통해 입증한다.

ABSTRACT

In this paper, we study damped second-order dynamics, which are quasilinear hyperbolic partial differential equations (PDEs). This is inspired by the recent development of second-order damping systems for accelerating energy decay of gradient flows. We concentrate on two equations: one is a damped second-order total variation flow, which is primarily motivated by the application of image denoising; the other is a damped second-order mean curvature flow for level sets of scalar functions, which is related to a non-convex variational model capable of correcting displacement errors in image data (e.g. dejittering). For the former equation, we prove the existence and uniqueness of the solution. For the latter, we draw a connection between the equation and some second-order geometric PDEs evolving the hypersurfaces which are described by level sets of scalar functions, and show the existence and uniqueness of the solution for a regularized version of the equation. The latter is used in our algorithmic development. A general algorithm for numerical discretization of the two nonlinear PDEs is proposed and analyzed. Its efficiency is demonstrated by various numerical examples, where simulations on the behavior of solutions of the new equations and comparisons with first-order flows are also documented.

연구 동기 및 목표

• 최적화에서 두 번째 차수 역학의 뛰어난 수치 성능에 기반하여, 개선된 영상 복원을 위한 두 번째 차수 기하학적 비선형 쌍곡형 PDE를 개발하는 것.
• 영상 노이즈 제거 및 이동 오차 보정(예: dejittering)을 위해 감쇠된 두 번째 차수 흐름을 총변동성 및 평균 곡률 흐름으로 확장하는 것.
• 이러한 새로운 PDE의 해석적 및 수치적 기초를 확립하는 것, 특히 해의 존재성, 유일성 및 수렴성에 대해 논의하는 것.
• 수치 실험을 통해 두 번째 차수 흐름이 일계 경사 흐름보다 경계 유지 능력과 진동 감소 성능에서 뛰어나다는 것을 입증하는 것.

제안 방법

• 감쇠된 두 번째 차수 총변동성 흐름(TVF)을 쌍곡형 PDE로 제안: $\ddot{w} + \eta(t)\dot{w} = -\partial\Phi(w)$, 여기서 $\Phi(w)$는 총변동성 함수기능이다.
• 레벨 집합을 위한 감쇠된 두 번째 차수 평균 곡률 흐름을 도입하여, 곡면의 표면적 최소화 및 이동 오차 보정을 모델링한다.
• 에너지 안정성과 수치 수렴성을 보장하기 위해 시간 이산화에 심플렉틱 오일러 스킴을 사용한다.
• 행렬 형태의 반연속 시스템 $\mathbf{z}_{k+1etro{A}_k \mathbf{z}_k$를 유도하며, $\mathbf{A}_k$는 시간 스텝 연산자이다.
• 거스터긴 원판 정리와 고유값 분석을 활용하여, CFL 유사 조건 하에서 이산화 스킴이 수축성임을 증명한다.
• 잘 정의된 문제와 수치 안정성을 확보하기 위해 평균 곡률 흐름 PDE의 정규화된 버전을 구현한다.

실험 결과

연구 질문

• RQ1 노이즈 제거 및 dejittering과 같은 영상 복원 작업을 위한 두 번째 차수 기하학적 PDE를 제작하고 분석할 수 있는가?
• RQ2 이러한 새로운 두 번째 차수 쌍곡형 PDE의 해에 대해 존재성과 유일성은 어떻게 보장되는가?
• RQ3 이러한 두 번째 차수 흐름의 해가 일계 경사 흐름과 비교해 수렴 속도와 경계 유지 능력 측면에서 어떻게 다를 수 있는가?
• RQ4 두 번째 차수 TVF 프레임워크에서 단순한 초기 자료에 대해 해석적 해를 유도할 수 있는가?
• RQ5 반연속 수치 스킴의 안정성과 수렴성을 보장하기 위한 시간 스텝에 대한 조건은 무엇인가?

주요 결과

• 감쇠된 두 번째 차수 총변동성 흐름은 유일한 해를 가지며, 단순한 초기 자료에 대해서는 명시적인 해석적 해가 유도된다.
• 정규화된 두 번째 차수 평균 곡률 흐름에 대해서는 해의 존재성과 유일성이 엄밀히 증명된다.
• 심플렉틱 오일러 방법에 기반한 수치 스킴은 감쇠 계수 $\eta$와 시스템 행렬의 최대 고유값을 포함하는 시간 스텝 조건 하에서 수렴성이 입증된다.
• 이산 반복 행렬 $\mathbf{A}_k$의 고유값이 $|\mu_{i,\pm}^k| \leq 1$을 만족함을 보여, 수축성과 수렴성을 보장한다.
• 수치 비교 결과, 두 번째 차수 흐름이 일계 흐름보다 더 빠른 수렴 속도를 보이며, 특히 경계가 풍부한 영역에서 더 효과적으로 경계를 유지함을 입증한다.
• 레벨 집합을 평균 곡률 흐름을 통해 진화시킴으로써 이미지의 원래 구조를 복원하는 방식으로, 이동 오차(예: dejittering)를 성공적으로 보정한다.