[논문 리뷰] The fixed-point iteration method for IMRT optimization with truncated dose deposition coefficient matrix

Zhen Tian, Masoud Zarepisheh외 2명|arXiv (Cornell University)|2013. 03. 14.

Advanced Radiotherapy Techniques참고 문헌 11인용 수 7

한 줄 요약

이 논문은 메모리 및 계산 비용을 줄이기 위해 흔히 사용되는 잘라낸 복사선 부여 계수(Dose Deposition Coefficient, DDC) 행렬을 사용할 때 강도 조절 방사선 치료(IMRT) 치료 계획 품질을 향상시키기 위해 고정점 반복 방법을 제안한다. 이 방법은 전체 DDC 해와 정확히 일치하지는 않지만, 잘라낸 문제에 대한 난이한 해보다 훨씬 더 정확한 해로 수렴하므로, 전체 행렬 저장이 필요 없이도 계획 정확도를 향상시킨다.

ABSTRACT

In the treatment plan optimization for intensity modulated radiation therapy (IMRT), dose-deposition coefficient (DDC) matrix is often pre-computed to parameterize the dose contribution to each voxel in the volume of interest from each beamlet of unit intensity. However, due to the limitation of computer memory and the requirement on computational efficiency, in practice matrix elements of small values are usually truncated, which inevitably compromises the quality of the resulting plan. A fixed-point iteration scheme has been applied in IMRT optimization to solve this problem, which has been reported to be effective and efficient based on the observations of the numerical experiments. In this paper, we aim to point out the mathematics behind this scheme and to answer the following three questions: 1) whether the fixed-point iteration algorithm converges or not? 2) when it converges, whether the fixed point solution is same as the original solution obtained with the complete DDC matrix? 3) if not the same, whether the fixed point solution is more accurate than the naive solution of the truncated problem obtained without the fixed-point iteration? To answer these questions, we first performed mathematical analysis and deductions using a simplified fluence map optimization (FMO) model. Then we conducted numerical experiments on a head-and-neck patient case using both the simplified and the original FMO model. Both our mathematical analysis and numerical experiments demonstrate that with proper DDC matrix truncation, the fixed-point iteration can converge. Even though the converged solution is not the one that we obtain with the complete DDC matrix, the fixed-point iteration scheme could significantly improve the plan accuracy compared with the solution to the truncated problem obtained without the fixed-point iteration.

연구 동기 및 목표

• 메모리 및 계산 제약으로 인해 복사선 부여 계수(DDC) 행렬의 작은 요소를 잘라내는 것에 의해 발생하는 IMRT 치료 계획 품질 저하 문제를 해결하기 위해.
• IMRT 최적화의 맥락에서 잘라낸 DDC 행렬에 대해 고정점 반복 방법의 수학적 행동을 분석하기 위해.
• 고정점 반복 해가 잘라낸 DDC 행렬에서 직접 구한 표준 해보다 우수한지 평가하기 위해.
• 고정점 반복 해가 수렴하는지, 그리고 완전한 DDC 행렬에서의 최적 해에 근접하는지 확인하기 위해.

제안 방법

• 고정점 반복 스킴의 수학적 분석을 위해 단순화된 복사선 조절 지도 최적화(FMO) 모델을 사용한다.
• 고정점 반복 방법을 반복적으로 적용하여 잘라낸 DDC 행렬을 사용해 복사선 강도 값을 업데이트함으로써 잘라내기로 인한 근사 오차를 보정한다.
• 이 방법은 고정점 갱신으로서 표현된다: x^{k+1} = x^k + α(D^T D x^k - D^T d), 여기서 D는 잘라낸 DDC 행렬이다.
• 수치 실험을 통해 단순화된 모델과 원래 FMO 모델을 모두 사용하여 환자 한 명의 두부·경부 부위에 대해 방법을 검증한다.
• 반복 간 복사선 벡터의 변화를 모니터링하여 수렴성을 평가한다.
• 선형 복사선 분포도(DVH) 및 감마 분석 지표를 사용하여 고정점 해, 난이한 잘라낸 해, 전체 DDC 해 간의 계획 품질을 비교한다.

실험 결과

연구 질문

• RQ1 IMRT 최적화에서 잘라낸 DDC 행렬에 고정점 반복 알고리즘이 적용되었을 때 수렴하는가?
• RQ2 고정점 반복이 수렴할 경우, 그 결과 해가 완전한 DDC 행렬에서 얻은 해와 동일한가?
• RQ3 동일하지 않다면, 고정점 해가 잘라낸 DDC 행렬에서 직접 유도된 표준 해보다 더 정확한가?
• RQ4 DDC 행렬 잘라내기가 계획 품질에 어떤 영향을 미치며, 고정점 반복은 이 영향을 완화할 수 있는가?

주요 결과

• 적절한 DDC 행렬 잘라내기 조건 하에서 고정점 반복 방법은 수학적 분석과 수치 실험을 통해 수렴하는 것으로 확인되었다.
• 수렴한 고정점 해는 완전한 DDC 행렬에서 얻은 해와 정확히 일치하지 않으며, 잘라내기로 인해 다른 고정점을 가지게 됨을 시사한다.
• 완전한 DDC 해와 일치하지는 않지만, 고정점 해는 잘라낸 행렬에서 직접 유도된 난이한 해보다 훨씬 더 정확하다.
• 두부·경부 환자 사례에 대한 수치 결과는 고정점 방법이 계획 품질을 향상시킴을 보여주며, 더 나은 복사선 분포와 목표 부위의 냉전역 감소로 나타난다.
• 이 방법은 DDC 행렬 잘라내기로 인한 오차를 효과적으로 줄여주며, 전체 행렬 저장이 필요 없이도 계획 정확도를 향상시킨다.