[논문 리뷰] Sinogram constrained TV-minimization for metal artifact reduction in CT

Clemens Schiffer, Kristian Bredies|arXiv (Cornell University)|2014. 04. 26.

Advanced X-ray and CT Imaging참고 문헌 2인용 수 3

한 줄 요약

이 논문은 금속 임플랜트로 인해 손상된 시노그램 데이터에 부등식 제약 조건을 적용하면서도 영상 재구성에 총변동(TV) 정규화를 사용하여 CT 영상에서 금속 아티팩트를 감소시키는 새로운 시노그램 제약 TV 최소화 방법을 제안한다. 챔볼-포크 알고리즘을 사용하여, 인공 쉬프-로건 페리티브에서 금속 임플랜트 근처의 해부학적 세부 정보를 성공적으로 유지하였으며, 딱딱한 제약 조건 하에서는 PSNR 47.6 dB, 노이즈가 있는 데이터에서는 40.1 dB의 성능을 달성하였다.

ABSTRACT

A new method for reducing metal artifacts in X-ray computed tomography (CT) images is presented. It bases on the solution of a convex optimization problem with inequality constraints on the sinogram, and total variation regularization for the reconstructed image. The Chambolle-Pock algorithm is used to numerically solve the discretized version of the optimization problem. As proof of concept we present and discuss numerical results for synthetic data.

연구 동기 및 목표

• 고밀도 임플랜트로 인한 비균일성과 절단으로 인한 CT 영상에서의 금속 아티팩트를 해결하기 위해.
• 금속 물체로 인해 일부 또는 손상된 시노그램 데이터가 존재할 경우의 영상 재구성 품질 향상.
• 총변동 정규화와 시노그램 수준의 제약 조건을 결합한 변분 접근법을 개발하여 금속 근처의 구조적 세부 정보를 유지하기 위해.
• 반복 최적화를 사용하여 시뮬레이션 데이터를 바탕으로 제안된 방법의 타당성을 입증하기 위해.

제안 방법

• 손상되지 않은 시노그램 영역에 대한 데이터 일치 항과 영상에 대한 총변동(TV) 정규화 항을 포함하는 볼록 최적화 문제를 수립한다.
• 시노그램에 대한 지점별 부등식 제약 조건을 도입하여, 금속 영향을 받는 영역에서 재구성된 영상 투과값이 최소 감지 가능한 강도를 나타내는 임계값 C 이하로 떨어지지 않도록 보장한다.
• 비부드러운 TV와 지표 함수를 효과적으로 처리할 수 있는 챔볼-포크 원-이중 알고리즘을 사용하여 최적화 문제를 해결한다.
• 디지털화된 라돈 변환과 그 쌍대 연산자를 사용하며, 수렴을 보장하기 위해 전처리 인자(Preconditioning factor)를 적용한다.
• 금속 영향을 받는 시노그램 영역의 값이 최소한 C 이상이 되도록 보장하기 위해 투영 단계를 통해 제약 조건을 구현한다.
• 근접 연산자 기반의 이중 형식을 통해 효율적인 구현이 가능하며, TV에 대한 소프트 스레시홀딩과 시노그램 제약 조건에 대한 클리핑 연산이 포함된다.

실험 결과

연구 질문

• RQ1 표준 TV 최소화 방법과 비교해 볼 때, 시노그램에 대한 부등식 제약 조건이 금속 영향을 받는 CT 스캔에서 영상 재구성 품질을 향상시키는가?
• RQ2 시노그램 수준의 제약 조건이 금속 임플랜트 근처의 미세한 해부학적 구조를 유지하는 데 어떤 영향을 미치는가?
• RQ3 다양한 데이터 일치 항과 정규화 파라미터의 조합이 아티팩트 감소와 영상 품질에 어떤 영향을 미치는가?
• RQ4 제안된 방법이 노이즈가 있거나 잘린 시노그램 데이터를 효과적으로 처리하면서도 진단 가능한 영상 품질을 유지할 수 있는가?

주요 결과

• 딱딱한 제약 조건을 적용한 제안된 방법은 10×10 크기의 금속 삽입부가 있는 시뮬레이션된 쉬프-로건 페리티브에서 PSNR 47.6 dB를 달성하였으며, 잘린 시노그램으로부터의 표준 재구성 방식보다 유의미하게 뛰어난 성능을 보였다.
• 5% 가우시안 노이즈가 존재하는 조건에서, 80,000회 반복 후 PSNR 40.1 dB를 달성하여 노이즈에 대한 강건성을 입증하였다.
• 표준 재구성 방식에서 완전히 가려졌던 금속 임플랜트 근처의 미세한 조직 구조를 성공적으로 유지하였다.
• 챔볼-포크 알고리즘이 비부드러운 제약 조건이 있는 최적화 문제를 안정적이고 수렴 가능하게 해결하는 데 기여하였다.
• 높은 영상 품질에도 불구하고, 80,000회 반복에 대해 GPU에서 약 1시간의 계산 시간이 소요되어 향후 가속화가 필요함을 시사하였다.
• 결과적으로, 이 방법은 금속 아티팩트 외에도 다른 형태의 손상되거나 누락된 투과 데이터에 일반화 가능할 것으로 보인다.