[논문 리뷰] Weighted Diffeomorphic Density Matching with Applications to Thoracic Image Registration

Caleb Rottman, Martin Bauer외 2명|arXiv (Cornell University)|2015. 08. 11.

Medical Imaging Techniques and Applications참고 문헌 7인용 수 3

한 줄 요약

이 논문은 확장성과 공간에 따라 변하는 조직 압축성에 따라 질량 보존을 강제하는 계산적으로 효율적인 디퍼모르픽 밀도 매칭 프레임워크를 제안한다. 이는 디퍼모르피즘의 밀도 작용과 압축성에 따라 변화하는 자코비안 행렬에 대한 페널티를 사용한다. 이 방법은 표준 LDDMM보다 질량을 보존하고 부피 변화를 탄성 있는 조직에 국한함으로써 성능을 뛰어넘으며, 런타임 속도가 10배 빠르다 (분당 400회 대비 45회).

ABSTRACT

In this article we study the problem of thoracic image registration, in particular the estimation of complex anatomical deformations associated with the breathing cycle. Using the intimate link between the Riemannian geometry of the space of diffeomorphisms and the space of densities, we develop an image registration framework that incorporates both the fundamental law of conservation of mass as well as spatially varying tissue compressibility properties. By exploiting the geometrical structure, the resulting algorithm is computationally efficient, yet widely general.

연구 동기 및 목표

• 표준 LDDMM이 질량 보존에 실패하고 조직 특이적 압축성을 모델링하지 못하는 한계를 해결한다.
• 숨쉬는 주기 동안 폐는 매우 탄성 있고 연조직/골격은 거의 비압축성인 흉부解부학의 물리적 현실을 모델링한다.
• 질량 보존과 해부학적 조직 특성을 존중하는 계산적으로 효율적인 정렬 프레임워크를 개발한다.
• 디퍼모르피즘의 리만 기하학과 밀도를 통합하여 안정적이고 기하학적으로 일致하는 변환을 가능하게 한다.
• 동적 흉부 영상에서의 변형 추정을 향상시켜 방사선 옹념에서 정확한 4차원 치료 계획 수립을 가능하게 한다.

제안 방법

• 디퍼모르피즘의 밀도 작용(즉, 밀도 작용)을 사용하여 체적 형식에 대한 당겨짐 작용을 통해 이미지 정렬을 밀도 매칭 문제로 공식화함으로써 질량 보존을 확보한다.
• 디퍼모르피즘 기하학과 밀도 변환 간의 연결을 정의하기 위해 밀도 공간에서 피셔-레이 메트릭을 사용하여 리만 기하학적 구조를 정의한다.
• 디퍼모르피즘의 자코비안 행렬에 공간에 따라 변하는 페널티 함수 $ f(x) $ 를 적용하며, $ f(x) $ 는 초기 CT 강도를 기반으로 로지스틱 함수를 통해 조직의 압축성을 반영한다.
• 부드러움과 일치성의 균형을 이루는 에너지 기능을 최소화하기 위해 소볼레프 $ H^1 $ 메트릭을 사용하여 디퍼모르피즘 공간에서 지오데식 경로를 구성한다.
• 지오데식 방정식의 통합이 필요 없도록 PyCA 소프트웨어 패키지를 사용하여 GPU 기반으로 알고리즘을 구현함으로써 상당한 속도 향상을 달성한다.
• 추정된 디퍼모르피즘을 통해 내쉬기 이미지를 푸시포워드하여 들이쉬기 이미지를 합성함으로써 참값과 직접 비교할 수 있도록 한다.

실험 결과

연구 질문

• RQ1 대규모 폐 부피 변화가 있는 흉부 CT 이미지 정렬에서 질량을 보존하는 디퍼모르픽 정렬 프레임워크를 설계할 수 있는가?
• RQ2 폐에서의 높은 압축성과 뼈에서의 비압축성과 같은 공간에 따라 변하는 조직 압축성을 이미지 정렬에서 효과적으로 모델링할 수 있는가?
• RQ3 표준 $ L^2 $ 이미지 작용 대신 디퍼모르피즘의 밀도 작용을 사용하면 더 정확하고 물리적으로 타당한 변형 장을 얻을 수 있는가?
• RQ4 자코비안 행렬에 공간에 따라 변하는 페널티를 적용하면 정확도는 높이고 계산 효율성은 유지할 수 있는가?
• RQ5 동적 흉부 영상 데이터에서 제안된 방법이 표준 LDDMM에 비해 속도, 정확도, 해부학적 타당성 측면에서 얼마나 뛰어나게 성능을 냈는가?

주요 결과

• 제안된 방법은 표준 LDDMM 대비 10배 빠른 속도를 기록하여 분당 400회, LDDMM는 분당 45회로 실행된다.
• 공간에 따라 변하는 페널티를 적용한 밀도 작용은 부피 변화를 주로 폐와 체외로 국한시키며, 비압축성 조직에서 비물리적인 압축/팽창을 방지한다.
• 합성된 들이쉬기 이미지와 진짜 들이쉬기 이미지 간의 절대 차이는 LDDMM에 비해 크게 감소하였으며, 폐 부위에서 오차가 최소화되었다.
• 공간에 따라 변하는 페널티가 없을 경우에도 LDDMM보다 개선되었지만, 폐 외부에서 비물리적인 부피 변화가 나타나, 적응형 페널티의 필요성을 확인한다.
• 자코비안 행렬 지도는 제안된 방법이 고밀도 구조물(예: 뼈)에서 비압축성을 유지하면서도 저밀도 영역(예: 폐)에서는 압축성을 허용함을 확인한다.
• 에너지 플롯과 로지스틱 페널티 함수는 안정적인 수렴을 보이며, 초기 CT 강도를 기반으로 압축성과 비압축성을 효과적으로 구분함을 확인한다.