[논문 리뷰] Fast Backprojection Techniques for High Resolution Tomography

Nikolay Koshev, Elias S. Helou외 1명|arXiv (Cornell University)|2016. 08. 11.

Medical Imaging Techniques and Applications참고 문헌 19인용 수 6

한 줄 요약

이 논문은 고해상도 단층촬영을 위한 새로운 빠른 역백프로젝션 기법을 제안한다. 로그-극좌표계와 반극좌표계를 사용하여 계산 복잡도를 감소시킨다. 역백프로젝션 연산자를 원형 적분으로 재표현하고, 로그-극좌표계에서 델타 함수 성질을 활용함으로써, O(N² log N)의 복잡도를 달성한다. 이는 기존의 O(N³) 방법에 비해 대규모 싱크로트론 단층촬영 데이터의 재구성 속도를 크게 향상시킨다.

ABSTRACT

Fast image reconstruction techniques are becoming important with the increasing number of scientific cases in high resolution micro and nano tomography. The processing of the large scale three-dimensional data demands new mathematical tools for the tomographic reconstruction task because of the big computational complexity of most current algorithms as the sizes of tomographic data grow with the development of more powerful acquisition hardware and more refined scientific needs. In the present paper we propose a new fast back-projection operator for the processing of tomographic data and compare it against other fast reconstruction techniques.

연구 동기 및 목표

• 싱크로트론 영상에서 고해상도 단층촬영 재구성의 증가하는 계산 부담을 해결한다.
• 기존의 O(N³) 복잡도를 가지는 표준 역백프로젝션을 개선하기 위해, 다른 좌표계에서 연산자를 재정의한다.
• 대규모 3차원 단층촬영 데이터에 적합한 수학적으로 엄밀하고 계산 효율성이 높은 역백프로젝션 방법을 개발한다.
• 이미지 품질 손실 없이 마이크로 및 나노 스케일 단층촬영에 대해 더 빠른 재구성 가능하게 한다.

제안 방법

• 기하학적 대칭성을 활용하기 위해 반극좌표계와 로그-극좌표계를 사용하여 역백프로젝션 연산자를 재구성한다.
• 역라돈 변환 커널의 역수를 사용하여, 원점과 점 x를 지나는 원을 따라 적분하는 방식으로 역백프로젝션을 표현한다.
• 로그-극좌표계에서 델타 함수 항등식을 활용하여 적분을 단순화하고, 문제를 일차원 근의 찾기 및 평가 단계로 축소한다.
• 디랙 델타 함수의 걸러내기 성질을 통해 로그-극좌표계에서의 역백프로젝션 이미지에 대한 닫힌 형태의 표현식을 유도한다.
• 수치적 적분을 효율적으로 수행하기 위해, 싱그램 데이터를 로그-극좌표 표현으로 매핑하기 위한 좌표 변환을 적용한다.
• 표준 카르테시안 좌표계에서의 역백프로젝션과의 성능 비교를 위해 시펄-로그런 페phant를 사용하여 방법을 검증한다.

실험 결과

연구 질문

• RQ1 역백프로젝션 연산자는 로그-극좌표계로 재표현될 수 있는가? 이는 계산 복잡도 감소에 기여하는가?
• RQ2 원형 적분 경로의 사용은 정확도를 유지하면서 재구성 속도를 향상시키는가?
• RQ3 원점과 이미지 점을 지나는 원을 따라 적분하는 역백프로젝션의 분석적 형태는 무엇인가? (점 소스에 대해)
• RQ4 고해상도 단층촬영에서 기존의 O(N³) 역백프로젝션과 비교해 본다면, 제안된 방법은 속도와 정확도 측면에서 어떻게 다른가?
• RQ5 역백프로젝션을 원점과 이미지 점을 지나는 원을 따라 적분하는 것으로 기하학적으로 해석할 수 있는가? 이는 알고리즘 최적화에 활용될 수 있는가?

주요 결과

• 제안된 역백프로젝션 방법은 O(N² log N)의 계산 복잡도를 달성하여 기존의 표준 O(N³) 복잡도에 비해 크게 향상되었다.
• 로그-극좌표계에서의 분석적 해는 역백프로젝션 이미지에 대한 닫힌 형태의 표현식을 제공한다: $ b(e^{ ho}m{ u}_{ heta}) = \frac{1}{\sqrt{(\cos(\theta-\phi)e^{\rho})^{2} + (e^{A} - \sin(\theta-\phi)e^{\rho})^{2}}} $.
• 역백프로젝션 기여가 원점과 이미지 점을 지나는 원에서 유래한다는 기하학적 성질을 활용하여, 효율적인 수치적 적분이 가능해졌다.
• 반극좌표계와 로그-극좌표계의 사용은 고해상도 환경에서 라돈 변환의 역변환을 더 자연스럽고 효율적으로 표현할 수 있도록 했다.
• 유도된 표현식을 통해 최적화된 좌표 변환과 델타 함수 적분을 통해 수행해야 할 연산 수를 줄여 더 빠른 재구성 가능해졌다.
• 시펄-로그런 페phant을 사용한 검증을 통해, 계산 부담 감소와 함께 정확한 재구성 결과를 도출하였다.