[논문 리뷰] H\"older stability of quantitative photoacoustic tomography based on partial data

Faouzi Triki, Qi Xue|arXiv (Cornell University)|2021. 03. 30.

Photoacoustic and Ultrasonic Imaging참고 문헌 35인용 수 3

한 줄 요약

이 논문은 경계 근처의 부분 내부 데이터로부터 양자화 광음향 단층촬영(PAT)에서 산란 계수 D(x)와 흡수 계수 µ(x)를 복원하는 데 있어 허더 안정성(Hölder stability)을 확립한다. 이는 √D에 대한 코시 문제의 새로운 안정성 분석을 통해 이루어지며, 하위영역이 경계 쪽으로 수축함에 따라 안정성이 향상되며, 허더 지수는 하위영역에 관계없이 일정한 양의 상수로 수렴한다. 수치적 복원 결과는 매끄럽고 불연속적인 매질 모두에서 상대 오차가 낮은(3–8.5%) 높은 정확도를 보이며 검증되었다.

ABSTRACT

We consider the reconstruction of the diffusion and absorption coefficients of the diffusion equation from the internal information of the solution obtained from the first step of the inverse photoacoustic tomography (PAT). In practice, the internal information is only partially provided near the boundary due to the high absorption property of the medium and the limitation of the equipment. Our main contribution is to prove a H\"older stability of the inverse problem in a subregion where the internal information is reliably provided based on the stability estimation of a Cauchy problem satisfied by the diffusion coefficient. The exponent of the H\"older stability converges to a positive constant independent of the subregion as the subregion contracts towards the boundary. Numerical experiments demonstrates that it is possible to locally reconstruct the diffusion and absorption coefficients for smooth and even discontinuous media.

연구 동기 및 목표

• 組織 흡수 및 측정 한계로 인해 내부 데이터 H(x)가 경계 근처에서만 이용 가능한 경우 양자화 PAT에 대한 안정성 분석의 부족을 해결한다.
• 측정 경계 근처의 하위영역에서 부분 내부 데이터를 사용하여 산란 계수 D(x)와 흡수 계수 µ(x)를 복원하는 데 있어 허더 안정성을 입증한다.
• 기존의 코시 문제 안정성 추정치가 하위영역이 경계에 가까워질수록 발산하는 문제를 해결하기 위해 하위영역 크기와 무관한 bound를 유도한다.
• 실제 부분 데이터 조건 하에서 매끄럽고 불연속적인 매질 모두에 대해 국소 복원의 가능성을 입증한다.

제안 방법

• 내부 데이터 H(x) = µ(x)u(x)를 바탕으로 √D에 대한 코시 문제를 수립하며, Γ ≡ 1을 가정한다.
• 해의 비율 uⱼ/u₁을 이용해 ∇lnσ(σ = D|u₁|²)를 복원하기 위해 Hⱼ/H₁로부터 ∇lnσ를 포함하는 방정식계를 유도한다.
• 노이즈 증폭을 완화하기 위해 수치 미분 및 필터링 기법을 사용하여 타원형 방정식의 계수 및 소스 항을 복원한다.
• 알려진 영역 근처에 영향을 최소화하기 위해 배경 또는 평균 값으로 코시 문제의 누락된 데이터를 외삽한다.
• H(x)가 노이즈에 의존하는 임계값을 초과하는 영역에서만 √D에 대한 타원형 PDE를 안정적인 수치적 방법으로 풀도록 한다.
• 직접 σ를 미분하는 것을 피하기 위해 식 Δ√σ / √σ = ½|∇lnσ|² + ½div(∇lnσ)를 이용해 계수 항을 계산한다.

실험 결과

연구 질문

• RQ1 내부 데이터 H(x)가 경계 근처에서만 이용 가능한 경우, 광음향 단층촬영(PAT)에서 D와 µ를 복원하는 역문제에 대해 허더 안정성을 확립할 수 있는가?
• RQ2 복원 하위영역이 측정 경계에 가까워질수록 안정성 추정치가 여전히 유계이면서 비퇴화되는가?
• RQ3 이 복원 방법은 부분 데이터 조건 하에서 매끄럽고 불연속적인 매질 모두를 다룰 수 있는가?
• RQ4 하위영역이 경계 쪽으로 수축함에 따라 안정성 추정의 허더 지수는 어떻게 행동하는가?
• RQ5 실제로 수치적 노이즈와 도함수 증폭이 복원 품질에 미치는 영향은 어느 정도인가?

주요 결과

• 하위영역이 경계 쪽으로 수축함에 따라 D와 µ를 복원하는 데 있어 허더 안정성 지수는 하위영역에 관계없이 일정한 양의 상수로 수렴하며, 조명원 근처에서 안정성이 향상됨을 나타낸다.
• 매끄러운 매질의 수치적 복원 결과, y > 0.2 영역에서 D의 상대 오차는 3.42%, µ의 상대 오차는 3.19%로 매우 높은 정확도를 보였다.
• 불연속적인 매질의 경우, 도메인 내 포함 영역의 상대 오차가 D에 대해 8.56%, µ에 대해 4.70%로 유의미하게 높지만, 불연속성 근처의 도함수 증폭에도 불구하고 성공적으로 복원되었다.
• 혈관 모델 복원 결과, D의 상대 오차는 4.74%, µ의 상대 오차는 2.47%로, 실제 생물학적 구조에서의 강건성을 확인하였다.
• 불연속성 경계에서 진짜 극단적 값이 필터링 기법으로 감소하더라도 복원이 안정적이고 정확하게 유지됨을 확인하여 수치적 노이즈에 대한 내성과 견고성을 입증하였다.
• H(x)가 노이즈 임계값을 초과하는 영역에서만 복원이 국소적으로 효과적으로 이루어지며, 이는 부분 데이터 조건 하에서의 실용적 타당성을 확인한다.