[논문 리뷰] Stability estimates in stationary inverse transport

Guillaume Bal, Alexandre Jollivet|ArXiv.org|2008. 04. 08.

Numerical methods in inverse problems참고 문헌 7인용 수 4

한 줄 요약

이 논문은 경계 데이터를 속도 벡터에 수직인 평면으로 매핑하는 기하학적 변환을 사용하여 정적 선형 운반 방정식에서 공간에 따라 변하는 산란 및 흡수 계수의 재구성에 대한 $L^1$-안정성 추정을 수립한다. 주요 기여는 이전 연구가 공간에 따라 변하지 않는 산란 계수에 국한되어 있었던 바, 차원 $n \geq 3$ 에서 엄밀한 안정성 결과를 확장한 것이다.

ABSTRACT

We study the stability of the reconstruction of the scattering and absorption coefficients in a stationary linear transport equation from knowledge of the full albedo operator in dimension $n\geq3$. The albedo operator is defined as the mapping from the incoming boundary conditions to the outgoing transport solution at the boundary of a compact and convex domain. The uniqueness of the reconstruction was proved in [M. Choulli-P. Stefanov, 1996 and 1999] and partial stability estimates were obtained in [J.-N. Wang, 1999] for spatially independent scattering coefficients. We generalize these results and prove an $L^1$-stability estimate for spatially dependent scattering coefficients.

연구 동기 및 목표

• 정적 선형 운반 방정식에서 공간에 따라 변하는 산란 및 흡수 계수의 재구성에 대한 $L^1$-안정성 추정을 증명하는 것.
• 이전 연구에서 공간에 따라 변하지 않는 산란 계수에 대해서만 부분적 안정성을 확립한 바를 확장하는 것.
• 입사 조건이 한 평면에서 출사 측정값으로 매핑되는 기하학적 프레임워크를 개발하여 분석을 단순화하는 것.
• 부하 조건 하에서 $L^1$ 공간에서 운반 연산자의 잘 정의됨과 유계성을 확립하는 것.
• 실제 영상 기하학에서 경계 측정값으로부터 안정적인 산란 및 흡수 계수 재구성의 기초를 마련하는 것.

제안 방법

• 원래 경계 $\Gamma_\pm$ 를 속도 벡터에 수직인 새로운 평면 $F_\pm$ 로 매핑하기 위해 반경 방향 연장 변환을 사용하여 단일 산란 기여의 분석을 단순화한다.
• 입사 데이터에서 출사 데이터로의 전진 운반 해를 나타내는 연산자 $\mathcal{A}$ 를 $L^1(F_-)$ 에서 $L^1(F_+)$ 로의 사상으로 정의한다.
• 함수해석적 접근을 통해 알보도 연산자를 특이 성분으로 분해하고, 지지 집합이 수축하는 테스트 함수를 통해 단일 산란 기여를 분리한다.
• 자유 운반 연산자 $\mathbf{T}_1$ 과 산란을 나타내는 $A_2$ 를 포함하는 연산자 $\mathbf{T} = \mathbf{T}_1 + A_2$ 를 기반으로 변형을 적용하여 유계 역가능성을 증명한다.
• 프레드홀름 대안과 노이만 급수를 적용하여 $\mathbf{T}^{-1}$ 이 $L^1(O, |v|\,dxdv)$ 에서 존재하고 유계임을 보이며, 운반 방정식의 잘 정의됨을 보장한다.
• 알보도 연산자의 연산자 노름 $\|\mathcal{A}\|_{\mathcal{L}(L^1(F_-), L^1(F_+))}$ 를 통해 계수의 차이를 유계로 제한함으로써 안정성을 확립한다.

실험 결과

연구 질문

• RQ1 정적 선형 운반 방정식에서 공간에 따라 변하는 산란 및 흡수 계수의 재구성에 대해 $L^1$-안정성 추정을 수립할 수 있는가?
• RQ2 측정 설정의 기하학적 구조가 역운동반 방정식 재구성의 안정성에 어떤 영향을 미치는가?
• RQ3 단일 산란 기여가 안정성 분석에서 어떤 역할을 하는가? 그리고 지지 집합이 수축하는 테스트 함수를 통해 어떻게 분리할 수 있는가?
• RQ4 운동반 연산자 $\mathbf{T}$ 가 어떤 조건에서 $L^1$-기반 함수 공간에서 역가능한가?
• RQ5 입사 및 출사 데이터가 속도 벡터에 수직인 평면에 정의된 변형된 기하학에서 알보도 연산자를 분석할 수 있는가?

주요 결과

• 공간에 따라 변하는 산란 및 흡수 계수의 재구성에 대해 $L^1$-안정성 추정이 수립되었으며, 이는 이전 연구에서 공간에 따라 변하지 않는 산란 계수에 국한된 결과를 일반화한 것이다.
• 적절한 부하 조건 하에서 계수에 대해 알보도 연산자 $\mathcal{A}$ 는 $L^1(F_-)$ 에서 $L^1(F_+)$ 로의 유계 연산자로 나타난다.
• 운동반 방정식의 해는 $L^1(O, |v|\,dxdv)$ 에서 유일하게 해가 존재하며, 해 연산자의 노름은 $\|\sigma\|_{\infty}$ 와 $\|\sigma_p\|_{\infty}$ 에 의해 제어된다.
• 운동반 연산자 $\mathbf{T}$ 의 역연산자는 $L^1(O, |v|\,dxdv)$ 에서 유계이므로 전진 문제의 안정성 보장과 함께 안정성 추정 유도가 가능하다.
• 안정성 추정은 $\mathbf{T}^{-1} = \mathbf{T}_1^{-1}(I + A_2\mathbf{T}_1^{-1})^{-1}$ 의 노이만 급수 전개를 통해 유도되었으며, $A_2\mathbf{T}_1^{-1}$ 의 연산자 노름은 $1 - e^{-2R\mathrm{v}_0^{-1}\|\sigma_p\|_{\infty}} < 1$ 으로 유계이므로 수렴이 보장된다.
• 최종 안정성 상한은 $\|f\|_{L^1(O,|v|dxdv)} \leq C_0 \|f_-–\|_{F_-}$ 로 표현되며, 여기서 $C_0$ 는 $R$, $\mathrm{v}_0^{-1}$, $\|\sigma\|_{\infty}$ 및 변형 연산자의 역연산자 노름에 의존한다.