[논문 리뷰] Inverse diffusion problems with redundant internal information

François Monard, Guillaume Bal|arXiv (Cornell University)|2011. 06. 21.

Numerical methods in inverse problems참고 문헌 10인용 수 5

한 줄 요약

이 논문은 타원형 방정식의 해 uᵢ에 대한 내부 측정값 Hᵢⱼ(x) = σ²ᵃ(x)|∇uᵢ|² 형태의 중복 측정값으로부터 스칼라 확산 계수 σ(x)를 복원하기 위한 두 가지 새로운 복원 방법을 제시한다. 첫 번째 방법은 과잉 결정된 상미분 방정식계를 풀어 2차원 및 3차원 결과를 임의의 차원으로 일반화하고, 거의 모든 α ∈ ℝ에 대해 적용 가능하다. 두 번째 방법은 특정 조건 하에서 유일하게 해를 갖는 선형 타원형 방정식계를 풀며, 기저 기반 경계 조건과 호환성 조건을 통해 안정적인 복원을 가능하게 한다.

ABSTRACT

This paper concerns the reconstruction of a scalar diffusion coefficient $σ(x)$ from redundant functionals of the form $H_i(x)=σ^{2α}(x)| abla u_i|^2(x)$ where $α\in\Rm$ and $u_i$ is a solution of the elliptic problem $ abla\cdot σ abla u_i=0$ for $1\leq i\leq I$. The case $α=\frac12$ is used to model measurements obtained from modulating a domain of interest by ultrasound and finds applications in ultrasound modulated electrical impedance tomography (UMEIT) as well as ultrasound modulated optical tomography (UMOT). The case $α=1$ finds applications in Magnetic Resonance Electrical Impedance Tomography (MREIT). We present two explicit reconstruction procedures of $σ$ for appropriate choices of $I$ and of traces of $u_i$ at the boundary of a domain of interest. The first procedure involves the solution of an over-determined system of ordinary differential equations and generalizes to the multi-dimensional case and to (almost) arbitrary values of $α$ the results obtained in two and three dimensions in \cite{CFGK} and \cite{BBMT}, respectively, in the case $α=\frac12$. The second procedure consists of solving a system of linear elliptic equations, which we can prove admits a unique solution in specific situations.

연구 동기 및 목표

• 중복 내부 측정값을 이용한 확산 계수 σ(x)의 안정적 복원을 2차원 및 3차원에서 임의의 공간 차원으로 일반화하고, 거의 모든 α ∈ ℝ에 대해 적용 가능하게 하는 것.
• UMEIT, UMOT, MREIT, ImpACT와 같은 하이브리드 영상 기법을 모델링하는 내부 기능성 Hᵢⱼ(x) = σ²ᵃ(x)∇uᵢ·∇uⱼ를 활용한 복원 프레임워크를 개발하는 것.
• 영역 내 각 점에서 n개의 벡터 Sᵢ = σᵃ∇uᵢ가 프레임을 이룰 수 있도록 보장함으로써 σ 복원의 유일성과 안정성을 확립하는 것.
• 중복 측정값으로 인해 발생하는 과잉 결정된 시스템의 호환성 조건을 분석하고 강제 적용하는 것.
• 측정값 Hᵢⱼ에 의해 유도된 계수를 갖는 새로운 타원형 방정식계를 제안하고, σ에 의존하지 않으며 프레드홀름 대체 조건 하에서 유일하게 해를 갖는다는 것을 증명하는 것.

제안 방법

• 원래의 역문제를 σ가 직접적으로 제거된 Sᵢ = σᵃ∇uᵢ에 대한 방정식계로 변환하는 것.
• ∇uᵢ가 컬-free임과 타원형 방정식 ∇·(σ∇uᵢ) = 0임을 이용해 Hᵢⱼ = Sᵢ·Sⱼ를 포함하는 비선형 과잉 결정된 시스템을 유도하는 것.
• 벡터 Sᵢ가 영역 내 모든 점에서 ℝⁿ에서 프레임을 이룰 수 있도록 uᵢ에 대한 경계 조건을 구성하는 것. 이는 역행성과 안정성 보장을 위해 필수적이다.
• 첫 번째 방법에서 과잉 결정된 상미분 방정식계를 풀어 이전의 2차원 및 3차원 결과를 고차원으로 일반화하는 것.
• 측정값 Hᵢⱼ에만 의존하는 계수를 갖는 선형 타원형 방정식계를 uᵢ에 대해 구성하고, 프레드홀름 대체 조건 하에서 유일한 해 존재성을 증명하는 것.
• 미분기하학적 제약 조건, 예를 들어 F = ∇logσ에 대해 dF^♭ = 0를 이용해 호환성 조건을 강제 적용하며, 이는 ∑ₚ,ₚ Mₚₚⁱʲ : ∇uₚ ⊗ ∇uₚ = 0로 표현된다.

실험 결과

연구 질문

• RQ1 2차원 및 3차원에서의 σ(x) 복원을 중복 내부 측정값을 이용해 임의의 공간 차원 n ≥ 2 및 거의 모든 α ∈ ℝ로 일반화할 수 있는가?
• RQ2 Sᵢ = σᵃ∇uᵢ가 영역 내 모든 점에서 ℝⁿ에서 프레임을 이룰 수 있도록 보장하기 위해 경계 데이터에 어떤 조건이 필요한가?
• RQ3 과잉 결정된 상미분 방정식계 및 타원형 방정식계의 호환성 조건을 어떻게 유도하고 강제 적용하여 일관성과 경로 독립성 보장을 확보할 수 있는가?
• RQ4 제안된 uᵢ에 대한 선형 타원형 방정식계가 σ에 독립적으로 유일하게 해를 갖는 조건은 무엇인가?
• RQ5 측정 기능성 Hᵢⱼ = σ²ᵃ|∇uᵢ|²를 사용하여 비중복 측정값에 의존하지 않고도 σ를 복원할 수 있는가? 비중복 측정값은 기존에 안정성이 떨어짐이 알려져 있음.

주요 결과

• 첫 번째 방법은 Sᵢ 벡터가 프레임을 이룰 경우, 임의의 차원 n ≥ 2 및 거의 모든 α ∈ ℝ에 대해 σ 복원의 전역적 유일성과 안정성을 보장한다.
• 프레임 조건과 내부 측정값 Hᵢⱼ를 바탕으로 유도된 과잉 결정된 상미분 방정식계는 해가 존재하며, σ의 안정적 복원을 가능하게 한다.
• 두 번째 방법은 측정값 Hᵢⱼ에만 의존하는 계수를 갖는 선형 타원형 방정식계를 구성하며, (n, α) 값이 이산적인 집합을 제외한 전 범위에서 유일한 해 존재성을 증명한다.
• 복원의 호환성 조건은 F = ∇logσ가 그라디언트여야 한다는 조건에 의해 표현되며, 이는 데이터가 ∑ₚ,ₚ Mₚₚⁱʲ : ∇uₚ ⊗ ∇uₚ = 0를 충족해야 한다는 것을 의미한다.
• 복소 기하학적 옹도(CGO) 해를 사용함으로써 n ≥ 3 차원에서 프레임에 독립적인 성질을 확보하여 적절한 경계 조건을 구성할 수 있다.
• 이 프레임워크는 이전의 2차원 및 3차원 결과(α = 1/2 및 α = 1)를 일반화하며, UMEIT, UMOT, ImpACT, MREIT 응용 분야를 모두 포함한다.