A Novel Differential Pathlength Factor Model for Near-Infrared Diffuse Optical Imaging

기술 요약: 근적외선 확산 광학 이미징을 위한 새로운 차등 경로 길이 계수(Differential Pathlength Factor) 모델

문제 제기
연속파 근적외선(CW-NIR) 확산 광학 이미징은 조직의 산소 포화도와 혈역학을 모니터링하기 위해 널리 사용되는 비용 효율적인 방식이다. 그러나 이 방식의 정량적 정확도는 측정된 광학적 감쇠를 크로모포어(chromophore) 농도로 변환하기 위해 정확한 차등 경로 길이 계수(DPF)를 필요로 하는 수정된 비어-람베르트 법칙(Modified Beer–Lambert Law, MBLL)에 크게 의존한다. 기존의 DPF 정의는 공식에 따른 한계점을 가지고 있어, 특히 반사 모드(reflectance mode)의 짧은 광원-검출기 간격에서 상당한 오류를 유발한다.

현재 사용되는 접근 방식은 다음과 같다:

1. 평균 경로 길이 DPF (Mean-pathlength DPF): 평균 광자 경로 길이를 광원-검출기 간격으로 나눈 값으로 정의된다. 이는 흡수 가중치 광자 통계(absorption-weighted photon statistics)를 무시하므로, 흡수율이 높은 조직에서 과대평가되는 경향이 있다.
2. 기울기 기반 DPF (Slope-based DPF): 광학 밀도 변화를 흡수 계수에 대해 미분한 기울기로 정의된다. 이론적으로는 근거가 있으나, 여전히 특정 농도 추정이라는 구체적인 목표보다는 평균적인 궤적에 초점을 맞추고 있다.
3. 상수/반무한 DPF (Constant/Semi-infinite DPF): 점근적 확산 이론에 따라 간격이 2.5 cm보다 큰 경우 상수로 근사되기도 한다. 이 가정은 DPF의 거리 의존성과 조직의 광학적 특성 변화를 무시하며, 실제 상황에서 100%를 초과하는 추정 오류를 빈번하게 발생시킨다.

이러한 부정확성은 신호 강도가 가장 높은 짧은 간격 영역에서 오류가 두드라지기 때문에, 넓은 영역의 검출이나 신호 깊이가 변하는 상황에서 CW-NIR 이미징의 적용 가능성을 제한한다.

방법론
저자들은 기존 모델을 평가하고 새로운 모델을 제안하기 위해 컴퓨터 시뮬레이션과 물리적 실험을 결합한 이중 접근 방식을 채택하였다.

• 컴퓨터 시뮬레이션 프레임워크: 균일한 조직 슬래브 기하 구조(29 × 29 × 5 cm³)에서 MCmatlab을 사용하여 몬테카를로 광자 전달 시뮬레이션을 수행하였다. 연구 범위는 흡수 계수($\mu_a$) 0 ~ 0.5 cm⁻¹ 및 산란 계수($\mu_s$) 0 ~ 500 cm⁻¹의 광범위한 광학적 특성을 포함한다. 각 조건에 대해 "참값(true)" 광학 밀도를 계산하였고, 6가지 서로 다른 DPF 공식을 적용하여 흡수 계수를 추정하였다. 성능 지표로는 추정된 흡수 계수와 참값 사이의 상대 오차를 사용하였다.
• 실험적 검증: CW-NIR 설정(750 nm LED 및 실리콘 포토다이오드 사용)을 사용하여 제어된 팬텀 실험을 수행하였다. 서로 다른 산란 영역을 생성하기 위해 이산화티타늄(TiO₂)을 농도별(0.25 mg/mL 및 1.00 mg/mL)로 도핑한 두 종류의 균일한 조직 모사 팬텀을 초흡수성 고분자 하이드로젤 매트릭스를 사용하여 제작하였다. 측정은 광원-검출기 간격을 1 cm에서 8.5 cm까지 변화시키며 확산 반사 기하 구조에서 수행되었다.

핵심 기여: 새로운 DPF 모델
본 논문은 시뮬레이션 데이터를 통해 도출된 두 가지 새로운 D형 DPF 모델을 소개한다:

1. 참 DPF ($DPF_{true}$): 알려진 참값 흡수 계수와 시뮬레이션된 광학 밀도를 사용하여 DPF를 구하도록 MBLL 식을 재배열하여 직접 정의된다. 이 모델은 흡수에 의한 지수 가중치를 포함한 광자 전달의 완전한 통계적 특성을 포착하는 이론적 골드 스탠다드 역할을 한다.
2. 역거리 DPF ($DPF_{inv}$): 참 DPF를 근사하기 위해 제안된 실용적인 해석적 모델이다. 광학 밀도가 광원-검출기 간격($d$)과 $\log_{10}(d)$의 선형 결합(0이 아닌 절편 포함)으로 정확하게 모델링될 수 있다는 관찰에 기초하여, 저자들은 다음과 같은 공식을 도출하였다:
   $$DPF_{inv}(d) = A + \frac{B}{d} + \frac{C\log_{10}(d)}{d}$$
   여기서 $A, B, C$는 조직의 흡수 및 산란 특성에 따라 달라지는 회귀 계수이다. 저자들은 매 측정마다 새로운 시뮬레이션을 수행할 필요 없이 완전히 경험적으로 구현할 수 있도록 이 계수들에 대한 경험적 관계식을 제공하였다.

결과

• 시뮬레이션 성능: 참 DPF는 테스트된 모든 조건에서 무시할 만한 오차를 보였다. 제안된 역거리 DPF는 대부분의 광학적 특성 공간과 간격 범위에서 10% 미만의 상대 오차를 달성하였다. 반면, 기존 모델들(상수, 평균 경로 길이, 기울기 기반, 반무한 모델)은 특히 광원-검출기 간격이 8 cm 미만인 구간에서 지속적으로 100%를 초과하는 오차를 나타냈다.
• 실험적 검증: 팬텀 실험은 시뮬레이션의 경향성을 확인해주었다. 광학 밀도 측정값은 몬테카를로 예측값과 일치하여 실험 설정을 검증하였다. 실험 데이터로부터 흡수 계수를 추정할 때, 참 모델과 역거리 모델은 측정된 간격 전반에서 20% 미만의 상대 오차를 유지하였다. 반대로 기존 모델들은 특히 짧은 간격에서 빈번하게 100%를 초과하는 오차를 보였다.
• 매개변수 거동: 역거리 DPF 모델의 계수($A, B, C$)는 광학적 특성에 따라 매끄럽고 예측 가능하게 변하는 것으로 나타나 견고한 경험적 보정이 가능함을 보여주었다.

의의 및 주장
본 논문은 제안된 역거리 DPF 모델이 물리적 엄밀함과 계산 효율성 사이의 최적의 균형을 제공한다고 주장한다. 단순한 상수 DPF나 단순 평균 경로 길이 비율에서 벗어남으로써, 새 모델은 특히 기존 모델들이 실패하는 짧은 간격 영역에서 CW-NIR 이미징의 정량적 정확성을 크게 향상시킨다.

저자들은 참 DPF가 정확도의 이론적 한계를 나타내는 반면, 역거리 DPF는 일상적인 임상 및 연구 응용 분야에 적합한 실용적인 대안임을 강조한다. 본 연구는 정확한 DPF 모델링이, 특히 반사 모드에서 크로모포어 농도 추정의 신뢰성을 높이는 데 결정적임을 입증한다. 이러한 연구 결과는 정교한 모델을 채택함으로써 CW-NIR 시스템의 단순함과 견고한 기능적/임상적 응용에 필요한 정량적 엄밀성 사이의 간극을 메울 수 있음을 시사한다. 논문은 향후 연구에서 불균일한 조직 구조를 다루고 이러한 모델링 전략을 시간 영역(time-domain) 및 주파수 영역(frequency-domain) 방식으로 확장해야 함을 언급하며 마무리된다.