[논문 리뷰] Cine cone beam CT reconstruction using low-rank matrix factorization: algorithm and a proof-of-princple study

Jian‐Feng Cai, Xun Jia외 4명|arXiv (Cornell University)|2012. 04. 16.

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

한 줄 요약

이 논문은 시간 해상도가 높은 콘빔CT(CT) 영상의 재구성에 대해, 시간적 일관성과 주기적인 호흡 패턴을 활용하여 희박한 시간 해상도의 콘빔CT 투과 영상에서 심플한 콘빔CT 영상(시네-4DCBCT)을 재구성하기 위한 새로운 저질서 행렬 분해 알고리즘을 제안한다. 이 방법은 영상 시퀀스를 U = LR로 모델링하며, L은 영상 특징을 나타내고 R은 시간 동적 특성을 나타낸다. L과 R의 웨이블릿 및 푸리에 변환에서의 희박성을 유도함으로써 고해상도 재구성을 달성하였으며, 상대 오차는 노이즈가 없는 경우 3.97%, 노이즈가 있는 경우 6.81%를 기록하여 단일 회전 시네-4DCBCT의 실현 가능성을 입증하였다. 이는 기존의 4DCBCT에 비해 더 높은 시간 해상도를 확보한 것이다.

ABSTRACT

Respiration-correlated CBCT, commonly called 4DCBCT, provide respiratory phase-resolved CBCT images. In many clinical applications, it is more preferable to reconstruct true 4DCBCT with the 4th dimension being time, i.e., each CBCT image is reconstructed based on the corresponding instantaneous projection. We propose in this work a novel algorithm for the reconstruction of this truly time-resolved CBCT, called cine-CBCT, by effectively utilizing the underlying temporal coherence, such as periodicity or repetition, in those cine-CBCT images. Assuming each column of the matrix $\bm{U}$ represents a CBCT image to be reconstructed and the total number of columns is the same as the number of projections, the central idea of our algorithm is that the rank of $\bm{U}$ is much smaller than the number of projections and we can use a matrix factorization form $\bm{U}=\bm{L}\bm{R}$ for $\bm{U}$. The number of columns for the matrix $\bm{L}$ constraints the rank of $\bm{U}$ and hence implicitly imposing a temporal coherence condition among all the images in cine-CBCT. The desired image properties in $\bm{L}$ and the periodicity of the breathing pattern are achieved by penalizing the sparsity of the tight wavelet frame transform of $\bm{L}$ and that of the Fourier transform of $\bm{R}$, respectively. A split Bregman method is used to solve the problem. In this paper we focus on presenting this new algorithm and showing the proof of principle using simulation studies on an NCAT phantom.

연구 동기 및 목표

• 기존의 단계별 분류 기반 4DCBCT가 저해상도 또는 높은 방사선 선량을 야기하는 바, 진정으로 시간 해상도가 높은 시네-4DCBCT 재구성 방법을 개발한다.
• 호흡 운동의 내재된 시간적 일관성—특히 주기성—을 활용하여 희박한 투과 영상에서 영상 재구성 품질을 향상시킨다.
• 다중 회전 4DCBCT 프로토콜에 비해 스캔 시간과 선량을 줄이기 위해 단일 기구 이동 회전만으로도 부피형 CBCT 영상을 재구성할 수 있도록 한다.
• 외부 호흡 모니터를 사용하지 않고도 투과 데이터에서 시간 행렬 R을 통해 직접 호흡 운동 신호를 추출할 수 있도록 한다.

제안 방법

• 이 방법은 시네-4DCBCT 영상 시퀀스를 U = LR로 모델링하며, U는 열로 구성된 재구성 영상, L은 영상 특징, R은 시간 동적 특성을 나타낸다.
• U의 질서는 L의 열 수에 의해 제약되며, 이는 이미지 간 시간 유사성을 암묵적으로 강제한다.
• L의 타이트 웨이블릿 프레임 변환에서의 희박성은 영상 특징 일관성과 노이즈 감소를 촉진한다.
• R의 푸리에 변환에서의 희박성은 호흡 패턴의 주기성을 강제로 유도하며, 호흡의 규칙성을 활용한다.
• 데이터 적합성과 정규화 항을 균형 잡는 최적화 문제를 해결하기 위해 스플릿 브레그먼 알고리즘을 사용한다.
• 실현 가능성을 검증하기 위해 노이즈가 없는 경우와 노이즈가 있는 경우 모두에서 NCAT 페트 페르몬 시뮬레이션 데이터에 적용하였다.

실험 결과

연구 질문

• RQ1 저질서 행렬 분해가 단일 회전 콘빔CT 투과 영상에서 시네-4DCBCT 영상을 효과적으로 재구성할 수 있는가?
• RQ2 웨이블릿 및 푸리에 변환에서의 희박성 제약을 적용할 경우 재구성 정확도와 시간 해상도는 어느 정도 향상되는가?
• RQ3 이러한 방법은 기존의 4DCBCT에 비해 영상 품질과 선량 효율성 측면에서 어떻게 비교되는가?
• RQ4 분해 과정에서 추출된 행렬 R은 외부 대체 측정 장치 없이도 신뢰성 있게 호흡 운동 신호를 인코딩할 수 있는가?

주요 결과

• 제안된 시네-4DCBCT 재구성 방법은 노이즈가 없는 시뮬레이션 케이스에서 상대 오차 3.97%를 기록하여 높은 정확도를 입증하였다.
• 노이즈가 있는 조건에서는 재구성 오차가 6.81%로 증가하였으며, 현실 조건에서도 안정성을 보였다.
• 이 방법은 다중 회전이나 느린 스캔이 필요 없이 단일 기구 이동 회전만으로도 고시간 해상도 재구성을 가능하게 하였다.
• 분해 과정에서 추출된 행렬 R은 호흡 운동을 성공적으로 인코딩하여 외부 대체 장치 없이도 투과 데이터에서 직접 호흡 신호를 추출할 수 있었다.
• 알고리즘은 시간적 일관성과 주기성을 효과적으로 활용하여, 기존의 4DCBCT에 비해 시간 해상도와 선량 효율성 측면에서 뛰어난 성능을 보였다.
• 임상용 시네-4DCBCT에 대한 잠재력은 있으나, 전체 3D 부피 재구성에서 계산적 과제가 여전히 남아 있다.