[논문 리뷰] Inverse Potential Problems for Divergence of Measures with Total Variation Regularization

Laurent Baratchart, Cristóbal Villalobos Guillén외 3명|arXiv (Cornell University)|2018. 09. 21.

Numerical methods in inverse problems참고 문헌 43인용 수 8

한 줄 요약

이 논문은 3차원에서 벡터 값 측도의 발산에 대한 역전위 문제를 다루며, 자기장 측정값에서 자화 분포를 복원하기 위해 총 변동(regularization)을 사용한다. 정규화 파rameter와 노이즈가 0으로 수렴할 때 희박하거나 단일 방향의 자화가 유일하게 복원될 수 있는 조건을 확립하였으며, 수치 결과는 진정한 듀플레로 수렴함을 확인한다.

ABSTRACT

We study inverse problems for the Poisson equation with source term the divergence of an $\\mathbf{R}^3$-valued measure, that is, the potential $\\Phi$ satisfies $$ \\Delta \\Phi= \ ext{div} \\boldsymbol{\\mu}, $$ and $\\boldsymbol{\\mu}$ is to be reconstructed knowing (a component of) the field grad $\\Phi$ on a set disjoint from the support of $\\boldsymbol{\\mu}$. Such problems arise in several electro-magnetic contexts in the quasi-static regime, for instance when recovering a remanent magnetization from measurements of its magnetic field. We develop methods for recovering $\\boldsymbol{\\mu}$ based on total variation regularization. We provide sufficient conditions for the unique recovery of $\\boldsymbol{\\mu}$, asymptotically when the regularization parameter and the noise tend to zero in a combined fashion, when it is uni-directional or when the magnetization has a support which is sparse in the sense that it is purely 1-unrectifiable. Numerical examples are provided to illustrate the main theoretical results.

연구 동기 및 목표

• 측도의 외부에서 자기장을 측정함으로써 3차원 자화 분포를 복원하기 위한 프레임워크를 개발한다.
• 벡터 측도의 발산을 포함하는 역전위 문제의 불안정성 문제를 다룬다.
• 측도의 총 변동을 통한 정규화를 적용할 때 자화가 유일하게 복원될 수 있는 충분한 조건을 확립한다.
• 복원 정밀도 향상을 위해 순수 1-비직선성 및 단일 방향성 지지 집합과 같은 기하학적 및 측도론적 희박성 개념을 탐색한다.
• FISTA 및 In-Crowd 알고리즘을 사용한 이산화된 듀플레 소스 모델에서의 수치 예제를 통해 이론적 결과를 검증한다.

제안 방법

• 측도 μ ∈ M(S)³에 대해 ΔΦ = div μ 를 만족하는 역문제를 수식화하며, S와 분리된 집합 Q에서 grad Φ 의 측정값으로부터 μ 를 복원한다.
• 역문제의 안정성을 확보하기 위해 총 변동 정규화를 사용하여 데이터 적합성과 μ 의 총 변동을 포함하는 기능을 최소화한다.
• 자기장이 코일-자유 및 발산-자유 성분의 합으로 표현될 수 있도록 포아송-헤지 분해를 적용하며, 이는 순환 측도와 연결된다.
• S-침묵 자화(외부에서 필드가 0이 되는 자화) 개념을 적용하고, 3차원 르베그 측도가 0이고 보완 집합이 무한한 구성 요소를 가지는 미세한 집합(slim set) 성질을 활용하여, 발산-자유 조건 하에서의 유일성 보장을 한다.
• 디지털 격자상의 듀플레 소스에 대해 FISTA 및 In-Crowd 알고리즘을 사용하여 수치적 복원을 구현하며, 정규화 파rameter λ 가 희박성 제어에 기여한다.
• 전진 연산자 A: μ ↦ v·b(μ)|Q 가 S-침묵 클래스에서 유계이자 단사임을 가정하여 해의 식별 가능성을 보장한다.

실험 결과

연구 질문

• RQ1 자화 μ 가 외부에서 자기장 측정값으로부터 유일하게 복원될 수 있는 기하학적 및 측도론적 조건(S의 지지 집합에 대해)은 무엇인가?
• RQ2 무한 차원 함수 공간에서 총 변동 정규화가 희박하거나 단일 방향의 자화 분포를 복원하는 데 어떻게 기여하는가?
• RQ3 미세한 집합 성질(3차원 르베그 측도가 0이고 보완 집합의 구성 요소가 무한함)이 S-침묵 자화가 반드시 발산-자유임을 보장하는 데 어떤 역할을 하는가?
• RQ4 정규화 파arameter λ 와 노이즈 수준 ε 이 동시에 0으로 수렴할 때 복원이 진정한 듀플레 구성으로 수렴할 수 있는가?
• RQ5 자화의 구조적 성질(예: 순수 1-비직선성 또는 조각별 단일 방향성)은 복원의 수렴성과 정확성에 어떤 영향을 미치는가?

주요 결과

• 유한한 수의 점 듀플레가 포함된 컴팩트하고 미세한 집합 S 상에서, 정규화된 측도 μλ,e 는 λ → 0 이고 노이즈 → 0 일 때 총 변동에서 진짜 듀플레 강도로 수렴한다.
• 진짜 듀플레 위치 주변의 공 위에서 정규화된 측도 μλ,e 의 총 변동은 진짜 듀플레 모멘트로 수렴한다: |μλ,e|(B(x(j),r)) → |vj| as λ, ε → 0.
• 진짜 듀플레 위치의 이웃 영역 외부에서 정규화된 측도 μλ,e 의 총 변동은 0으로 수렴한다: |μλ,e|(S∖∪jB(x(j),r)) → 0 as λ, ε → 0.
• 정규화된 측도 μλ,e 의 지지 집합은 유한 개의 해석적 표면, 곡선 및 점들의 합집합에 포함되며, 이는 해의 구조적 단순성을 나타낸다.
• 20개의 듀플레에 대한 수치적 복원 결과, λ = 10−12 에서 진짜 μ0 와의 상대적 총 변동 거리는 0.004로 나타나 수렴을 확인한다.
• 네 개의 분리된 구성 요소에 대해 조각별 단일 방향 자화가 적용된 경우, λ = 10−14 에서 상대적 거리는 0.015로 감소하여 구조적 가정 하에서도 강건함을 보여준다.