[논문 리뷰] Reconstruction Algorithms for Positron Emission Tomography and Single Photon Emission Computed Tomography and their Numerical Implementation

A. S. Fokas, Arieh Iserles외 1명|ArXiv.org|2004. 12. 04.

Medical Imaging Techniques and Applications참고 문헌 33인용 수 3

한 줄 요약

이 논문은 비선형 적분 가능 체계의 수학적 기법을 활용하여 PET 및 SPECT에 대한 새로운 분석적 재구성 알고리즘을 제시하며, 라돈 변환과 감쇠 라돈 변환의 정확한 역변환을 가능하게 한다. 삼차 스퍼인 보간을 통한 수치 구현은 시준형 필터링을 통해 기존의 기브스 아티팩트를 감소시키고, 임상 영상 응용 분야에서의 방법의 강건성을 검증함으로써 표준 필터 모형(예: 쉐프-로그란)에 대해 정확한 재구성을 달성한다.

ABSTRACT

The modern imaging techniques of Positron Emission Tomography and of Single Photon Emission Computed Tomography are not only two of the most important tools for studying the functional characteristics of the brain, but they now also play a vital role in several areas of clinical medicine, including neurology, oncology and cardiology. The basic mathematical problems associated with these techniques are the construction of the inverse of the Radon transform and of the inverse of the so called attenuated Radon transform respectively. We first show that, by employing mathematical techniques developed in the theory of nonlinear integrable equations, it is possible to obtain analytic formulas for these two inverse transforms. We then present algorithms for the numerical implementation of these analytic formulas, based on approximating the given data in terms of cubic splines. Several numerical tests are presented which suggest that our algorithms are capable of producing accurate reconstruction for realistic phantoms such as the well known Shepp--Logan phantom.

연구 동기 및 목표

• 비선형 적분 가능 체계의 고급 수학적 기법을 활용하여 PET 및 SPECT에서 방사성약물 분포를 정확한 분석 공식으로 재구성하는 것.
• 실용적 재구성을 위해 삼차 스퍼인 근사법을 활용한 이러한 분석 공식의 수치적 구현을 위한 것.
• 실제 영상 조건에서 실제 의료 필터 모형(예: 쉐프-로그란 필터 모형 및 흉부 모델 포함)에 대한 알고리즘의 유효성 검증을 위한 것.
• 기본적으로 기브스–윌브라함 현상과 같은 수치적 아티팩트를 타겟팅한 필터링 절차를 통해 감소시키기 위한 것.
• 중간 데이터의 굵은 샘플링(예: 200개 대비 10개)이 정확한 재구성에 충분한지 확인하여 계산 비용을 감소시키기 위한 것.

제안 방법

• PET의 역라돈 변환은 비선형 적분 가능 방정식의 수학적 방법을 활용하여 유도되며, 방사성약물 분포 g(x₁,x₂)의 정확한 분석 공식을 도출한다.
• SPECT의 감쇠 라돈 변환 역시 분석 기법을 활용하여 역변환되며, CT로 측정된 감쇠 f(x₁,x₂)를 고려한다.
• 수치적 구현은 측정된 데이터를 삼차 스퍼인으로 근사화하여 역변환의 안정적이고 정확한 계산을 가능하게 한다.
• 알고리즘은 스퍼인 보간 및 적분 관계를 통해 중간 함수 h(ρ,θ), fᶜᵐᵉ(ρ,θ), fˢᵐᵉ(ρ,θ), hᶜ(ρ,θ), hˢ(ρ,θ), I(τ,ρ,θ)를 계산한다.
• 2차원 격자에서 공식 (3.10)을 활용한 역변환을 수행하며, 최종 이미지 g(x₁,x₂)는 계산된 함수 r(τ,ρ,θ)로부터 재구성된다.
• 후처리 단계로 필터링을 수행한다: PET에는 평균 필터(a=0.005)를 5번 적용하고 임계값 처리(최댓값의 1/20 이하 값은 0으로 설정)를 적용하며, SPECT에는 중앙값 필터링과 임계값 처리를 적용한다.

실험 결과

연구 질문

• RQ1 비선형 적분 가능 체계에서 유도된 분석적 역변환 공식이 PET 및 SPECT의 역라돈 및 감쇠 라돈 변환에 효과적으로 적용될 수 있는가?
• RQ2 삼차 스퍼인 기반 수치적 구현이 쉐프-로그란 필터 모형과 같은 복잡한 의료 필터 모형을 고정밀도로 재구성할 수 있는가?
• RQ3 라돈 변환 데이터의 굵은 샘플링(예: 200개 대비 10개)이 재구성 정확도를 유지하고 계산 부담을 줄이는 데 얼마나 효과적인가?
• RQ4 평균 및 중앙값 필터링과 임계값 처리를 결합한 필터링 기법이 재구성된 이미지에서 기브스–윌브라함 아티팩트를 억제하는 데 얼마나 효과적인가?
• RQ5 제안된 알고리즘이 인간 흉부 및 심장근육을 모의한 임상적으로 관련성이 있는 필터 모형에 대해 정확한 재구성을 수행할 수 있는가?

주요 결과

• 제안된 분석적 재구성 알고리즘은 시각적 비교를 통해 쉐프-로그란 필터 모형과 임상 필터 모형을 고정밀도로 재구성함을 입증하였다.
• 삼차 스퍼인 근사법을 활용한 수치적 구현은 제한된 데이터 샘플링 조건에서도 역변환의 안정성과 정확성을 보장한다.
• 특히 PET에 대해 평균 필터링과 임계값 처리(a=0.005)를, SPECT에 대해 중앙값 필터링과 임계값 처리를 적용한 필터링 절차는 재구성된 이미지에서 기브스–윌브라함 아티팩트를 크게 감소시켰다.
• 중간 데이터를 t에 대해 10개의 등간격 점으로 샘플링하더라도 재구성 정확도가 유지됨을 확인하여 계산 비용을 크게 감소시키면서도 정밀도 손실 없이 작동함을 입증하였다.
• PET에 대해 500×500 격자, SPECT에 대해 140×140 격자에서 정확한 재구성을 달성하여 표준 영상 격자에서의 확장성과 강건성을 입증하였다.
• 비선형 적분 가능 체계 기반 수학적 프레임워크가 임상용 PET 및 SPECT 재구성에 실현 가능하고 정확한 길임을 확인하였다.