Efficient Beam-Hardening Correction Algorithm for Metal Artifact Reduction in Computed Tomography

Abstract

빔경화는 다색 X선을 사용하고 에너지 의존적인 물질 감쇠 계수를 이용하는 CT 시스템의 특성상 불가피한 현상이며, 이는 특히 금속 영역을 포함하는 프로젝션 상의 값을 오측정하여 결과적으로 CT 영상에 허상물을 유발한다. 금속 허상물 저감화는 CT 영상에 존재하는 이러한 허상물을 제거하고 가려진 실제 정보를 복원하는 과정이다. 영상을 통한 진단과 방사선치료를 위한 계획 수립에 있어서 정확한 CT 영상을 획득하기 위해 금속 허상물의 제거는 필수적이다. 반복적인 재구성에 의한 수치적 방법에 기반을 둔 효과적인 금속 허상물 제거에 관한 최신 연구들이 발표되었으나 무거운 계산량으로 인해 임상 실습에 적용이 어려운 상황이다.

본 논문에서는 이러한 계산적인 이슈를 해결하기 위한 효율적인 빔 경화 추정 모델과 이를 이용한 금속 허상물 저감화 방법을 제안한다. 제안한 모델은 금속 물체의 기하정보와 다색 X선이 물체를 통과하면서 발생하는 빔경화의 물리적인 특성을 반영한다. 모델에 필요한 대부분의 매개변수들은 수치학적인 방법으로 교정 전의 CT 영상과 CT 시스템으로부터 추가적인 최적화 과정 없이 획득한다. 빔경화 허상물과 관련된 매개 변수 중 단 하나만 재구성 이후의 영상 단계에서 선형 최적화를 통해 탐색된다. 또한 제안한 방법으로 교정된 결과 영상에 잔존하는 허상물들을 제거하기 위한 추가적인 두가지 개선 방법을 제시한다.

다수의 시뮬레이션 데이터와 실제 데이터를 사용하여 정성적 및 정량적 비교를 통해 제안 기법의 유효성이 체계적으로 평가되었다. 제안 알고리즘은 정확성 및 견고성 측면에서 유의미한 결과를 보여주었고, 기존의 기법들에 비해 향상된 결과 영상의 품질 뿐만 아니라 임상적으로 적용할만한 빠른 수행 시간을 보여주었다. 이 연구는 CT 영상을 통한 진단과 방사능 치료의 계획 수립을 위한 정확성 향상에 유의미한 의미를 갖는다.

Beam hardening in X-ray computed tomography (CT) is an inevitable problem due to the characteristics of CT system that uses polychromatic X-rays and energy-dependent attenuation coefficients of materials. It causes artifacts in CT images as the result of underestimation on the projection data, especially on metal regions. Metal artifact reduction is the process of reducing the artifacts in CT and restoring the actual information hidden by the artifacts. In order to obtain exact CT images for more accurate diagnosis and treatment planning on radiotherapy in clinical fields, it is essential to reduce metal artifacts. State-of-the-art approaches on effectively reducing metal artifact based on numerical methods by iterative reconstruction have been presented. However, it is difficult to be applied in clinical practice due to a heavy computational burden.

In this dissertation, we proposes an efficient beam-hardening estimation model and a metal artifact reduction method using this model to address this computational issue. The proposed model reflects the geometric information of metal objects and physical characteristics of beam hardening during the transmission of polychromatic X-ray through a material. Most of the associated parameters are numerically obtained from an initial uncorrected CT image and CT system without additional optimization. Only the unknown parameter related to beam-hardening artifact is fine-tuned by linear optimization, which is performed only in the reconstruction image domain. Two additional refinement methods are presented to reduce residual artifacts in the result image corrected by the proposed metal artifact reduction method.

The effectiveness of the proposed method was systematically assessed through qualitative and quantitative comparisons using numerical simulations and real data. The proposed algorithm showed significant results in the aspects of accuracy and robustness. Compared to existing methods, it showed improved image quality as well as fast execution time that is clinically applicable. This work may have significant implications in improving the accuracy of diagnosis and treatment planning for radiotheraphy through CT imaging.