[논문 리뷰] Fast OPED algorithm for reconstruction of images from Radon data

Yuan Xu, Oleg Tischenko|ArXiv.org|2007. 03. 21.

Medical Imaging Techniques and Applications참고 문헌 8인용 수 7

한 줄 요약

이 논문은 FFT 기반 이산 sine 변환과 선형 보간을 사용하여 Radon 데이터로부터 영상을 복원하기 위한 OPED(Orthogonal Polynomial Expansion on the Disk) 알고리즘의 빠른 구현을 제시한다. 이 방법은 원래 OPED 대비 26배 빠른 속도 향상을 이룩하면서도 복원 정확도를 유지하며, 수치 실험에서는 오차 차이가 거의 없고, 약간의 부드러운 함수에 대해 이론적으로 수렴성이 입증된다.

ABSTRACT

A fast implementation of the OPED algorithm, a reconstruction algorithm for Radon data introduced recently, is proposed and tested. The new implementation uses FFT for discrete sine transform and an interpolation step. The convergence of the fast implementation is proved under the condition that the function is mildly smooth. The numerical test shows that the accuracy of the OPED algorithm changes little when the fast implementation is used.

연구 동기 및 목표

• Radon 데이터로부터 영상 복원을 위한 OPED 알고리즘의 계산 효율적인 구현을 개발하기 위해.
• FFT와 보간 기법을 활용하여 OPED가 FBP와 동등한 속도로 경쟁할 수 있도록 하기 위해.
• 이미지 함수의 미약한 부드러움 조건 하에서 빠른 OPED 알고리즘의 수렴성을 증명하기 위해.
• 기본 편형체를 사용한 수치 실험을 통해 빠른 OPED의 정확성과 성능을 검증하기 위해.

제안 방법

• 빠른 OPED 알고리즘은 직접적인 수직 다항식 전개 계산을 대체하기 위해 FFT를 사용하여 이산 sine 변환을 계산한다.
• 비균일한 샘플링을 처리하기 위해 선형 보간 단계를 통합하여 효율적인 계산을 가능하게 한다.
• 알고리즘은 체비셰프 노드 위의 가우스 구적법에 기반하며, 재구성 공식은 각도 시야와 반경 구적 점에 대한 합으로 표현된다.
• 알고리즘은 제2종 체비셰프 다항식의 성질을 활용하여 FFT 가속을 가능하게 한다.
• 빠른 구현은 원래 OPED 공식의 연속 적분을 이산 구적과 보간으로 근사함으로써 유도된다.
• 이론적으로 재구성 연산자의 연산자 노름을 유계함으로써 수렴성을 확립하였으며, 보간 하에 있어서도 O(m log m)로 유지됨을 보였다.

실험 결과

연구 질문

• RQ1 OPED 알고리즘은 계산 속도를 FBP 수준으로 높일 수 있을까, 동시에 복원 정확도를 유지할 수 있을까?
• RQ2 OPED에 FFT와 선형 보간을 사용할 경우 원래 알고리즘의 수렴 성질이 유지되는가?
• RQ3 부드럽고 비부드러운 이미지에 대해 보간에 의해 유도된 오차는 원래 OPED와 비교해 얼마나 되는가?
• RQ4 빠른 OPED는 원래 OPED 대비 정량적으로 얼마나 빠르며, 오차는 어떤 식으로 상쇄되는가?
• RQ5 이차 또는 사차 도함수가 연속인 함수에 대해 빠른 OPED 알고리즘이 균일하게 수렴할 수 있는가?

주요 결과

• 빠른 OPED 알고리즘은 원래 OPED 대비 26배 빠른 속도 향상을 이룩하였으며, m=512일 때 복원 시간을 344초에서 13초로 단축시켰다.
• 원래 Shepp-Logan 편형체와 빠른 OPED 재구성 간의 상대 제곱 오차(RSE)는 0.00249574로, 원래 OPED의 0.00239702와 거의 동일한 수준이다.
• 빠른 OPED의 평균 오차(ME)는 0.00981329로, 원래 OPED의 0.0129175보다 낮아 평균 오차가 향상됨을 나타낸다.
• 두 재구성 간의 차이 이미지에서는 시각적·정량적 차이가 거의 없으며, 차이 RSE는 오직 0.000515499에 불과하다.
• 이론적 분석을 통해 빠른 OPED가 보간을 포함하더라도 C⁴(B)에 속하는 함수에 대해 O(m log m)로 수렴함을 확인하였으며, 원래 알고리즘의 수렴 속도를 유지함을 입증하였다.
• 수치 결과는 빠른 OPED가 Shepp-Logan 편형체와 같이 급격한 불연속성이 있는 이미지에 대해서도 원래 OPED와 거의 동일한 성능을 보임을 보여주었다.