기술 요약: 고유 교정 섭동 하에서의 척추 DRR에 대한 투영 영역 민감도 분석

문제 정의

정확한 기하학적 교정은 투시 가이드 척추 영상, 골격 디지털 재구성 방사선(DRR) 생성, 그리고 2D–3D 척추 등록을 위한 전제 조건이다. 전통적으로 교정 품질은 재구성 오차나 체적 재구성 정확도와 같은 재구성 영역 메트릭을 사용하여 평가된다. 그러나 최근의 많은 투시 가이드 등록 알고리즘은 재구성 정확도 자체보다는 획득된 방사선 영상과 합성된 DRR 사이의 일관성에 의존하며, 주로 투영 영역에서 작동한다.

최근 연구들은 교정 불확실성이 DRR 생성에 사용되는 영상 기하학에 영향을 미친다는 점을 시사하고 있으나, 고유 교정 섭동이 투영 영역의 일관성에 미치는 구체적인 영향은 충분히 이해되지 않았다. 이는 특히 복잡한 3D 구조, 겹쳐진 골격 특징, 제한된 시점을 포함하는 척추 해부학적 구조를 가진 척추 영상에서 매우 중요하다. 영상 기하학의 작은 변화는 척추의 윤곽, 추궁판(pedicle)의 외형, 종판(endplate) 정렬을 변화시켜, 이미지 유사도 메트릭과 등록 수렴 성능을 저하시킬 가능성이 있다. 본 연구는 해부학적 구조와 획득 포즈가 고정된 상태에서도 고유 교정 섭동이 어떻게 척추 투영 일관성과 다운스트림 2D–3D 등록 정확도에 영향을 미치는지 조사함으로써, 재구성 기반 교정 평가와 투영 영역 등록 동작 사이의 간극을 해결하고자 한다.

방법론

저자들은 척추 투시 영상의 고유 교정 섭동에 대한 투영 민감도를 분석하기 위해 통제된 합성 프레임워크를 제시한다.

1. 합성 데이터셋 및 기하학:

• 모델: 실험에는 격리된 척추 및 다중 척추 요추 분절(L1–L5)을 포함하여 CT 유래 척추 분할 마스크 및 메쉬 모델(VerSe, SpineWeb, TotalSegmentor 등의 데이터셋 활용)을 사용한다.
• 영상 설정: 고정된 선원-검출기 간 거리(SDD) 및 선원-등심점 간 거리(SID)를 가진 비등심(non-isocentric) C-arm 기하학을 시뮬레이션한다. 두 가지 표준 뷰인 전후방(AP) 및 측방(LAT) 뷰를 생성한다. 교정 변화의 영향을 해부학적 움직임으로부터 격리하기 위해 뷰 전반에 걸쳐 척추 모델의 포즈는 고정된다.
• DRR 생성: 선적분 감쇠 누적을 기반으로 한 레이 드라이븐(ray-driven) 투영 렌더링을 사용하여 1024 × 1024 해상도의 DRR을 생성한다.

2. 섭동 모델:

• 고유 파라미터: 외적 파라미터($R, t$)는 일정하게 유지하면서, 고유 교정 행렬($K$)에 대해서만 통제된 섭동을 적용한다.
• 특정 섭동:
  • 초점 스케일(Focal Scale): 확대율 불확실성을 시뮬레이션하기 위해 유효 초점 거리($f_x, f_y$)에 곱셈적 스케일링을 적용한다.
  • 관통점 변위(Piercing Point Displacement): 검출기 평면상의 주점($c_x, c_y$)에 대한 2D 이동을 적용한다.
• 비교: 그라운드 트루스 파라미터로 생성된 DRR을 섭동된 파라미터로 생성된 DRR과 비교한다.

3. 분석 메트릭:

• 투영 영역 일관성: 해부학적 랜드마크 변위, 윤곽 기반 거리(Chamfer, Hausdorff), 실루엣 중첩(Dice 계수, IoU), 그리고 원본 DRR 차이 히트맵을 사용하여 평가한다.
• 다운스트림 등록: 랜드마크 기반 2D–3D 강체 등록 실험을 수행하여 투영 불일치가 포즈 추정 오차로 어떻게 전파되는지 정량화한다. 메트릭으로는 병진 오차(mm), 회전 오차(degrees), 잔여 재투영 오차를 포함한다.

주요 기여

1. 합성 프레임워크: 해부학적 움직성으로부터 교정 효과를 격리하여 척추 투시 영상의 고유 교정 섭동에 대한 투영 영역 민감도를 분석하는 통제된 방법론을 도입하였다.
2. 체계적 정량화: 랜드마크 변위, 윤곽 변형, 실루엣 중첩, 이미지 유사성을 포함하는 상호 보완적인 메트릭을 사용하여 투영 일관성을 종합적으로 평가하였다.
3. 뷰 의존적 민감도 특성화: 측방(LAT) 척추 투영이 전후방(AP) 투영보다 고유 교정 변화에 훨씬 더 민감하다는 것을 입증하였다.
4. 등록에 대한 정량적 연결: 고유 교정 섭동, 투영 영역 불일치, 그리고 다운스트림 2D–3D 등록 오차 사이의 정량적 관계를 확립하였다.
5. 상호 보완적 평가: 투영 영역 분석이 기존의 재구성 기반 교정 평가와 구별되면서도 이를 보완하는 정보를 제공함을 입증하였다.

실험 결과

1. 정성적 및 전역적 일관성:

• 뷰 의존성: 동일한 섭동 크기(예: 100–500 픽셀 초점 거리 이동) 하에서 LAT 투영은 AP 투영에 비해 원본 DRR 차이 히트맵에서 현저히 더 강하고 공간적으로 확장된 불일치를 보인다.
• 실루엣 저하: 섭동 크기가 증가함에 따라 LAT 투영의 실루엣 일관성이 눈에 띄게 감소한다. 초점 길이 섭동이 100에서 500 픽셀로 증가할 때, LAT Dice 계수는 약 0.995에서 약 0.973으로 떨어지고, IoU는 약 0.991에서 약 0.948로 떨어진다. 반면, AP 투영은 높은 안정성을 유지한다(Dice > 0.997, IoU ~0.995).

2. 랜드마크 변위:

• 이동 크기: LAT 투영은 실질적으로 더 큰 해부학적 랜드마크 변위를 나타낸다. 예를 들어, 특정 섭동 하에서 LAT 뷰의 하부 종판과 극돌기 끝부분은 약 15–18 픽셀 이동한 반면, AP 랜드마크는 2 픽셀 미만이었다.
• 기하학적 원인: 이러한 LAT 뷰에서의 높은 민감도는 척추가 검출기에 더 가까운 비등심 설정에 기인하며, 이로 인해 더 큰 투영 스케일이 발생한다. 결과적으로 고유 교정 부정확성이 겉보기 랜드마크 변위를 증폭시킨다.

3. 2D–3D 등록 영향:

• 오차 전파: 교정 섭동에 의해 유도된 투영 불일치는 다운스트림 등록 오차로 전파된다.
• 회전 민감도: 회전 정렬 정확도가 특히 민한성을 보인다. 직접 등록 실험에서 초점 스케일 섭동이 증가함에 따라 회전 오차는 약 0.5°에서 약 3.8°로 증가하였다. Kabsch 기반 등록 또한 섭동에 따라 회전 오차가 증가(0.05°에서 0.37°)하는 것을 보여주었다.
• 병진 안정성: 병진 오차는 상대적으로 안정적이었으나, 섭동 크기에 따라 단조 증가하는 경향을 보였다.

의의 및 주장

본 논문은 척추 영상의 교정 품질을 평가할 때 재구성 기반 메트릭(예: 재투영 오차)에만 의존하는 것은 불충분하다고 주장한다. 본 연구는 다음을 입증한다:

• 작은 섭동의 중요성: 비교적 작은 고유 교정 변화라도 척추 투영 기하학, 윤곽 형태 및 랜드마크 국소화에 측정 가능한 변화를 일으킬 수 있다.
• 뷰 의존성의 중요성: 영상 시스템의 민감도는 균일하지 않다. 측방 뷰는 깊이 압축 및 해부학적 중첩으로 인해 본질적으로 교정 오류에 더 취чески 취약하다.
• 다운스트림 영향: 투영 영역의 불일치는 특히 회전 정렬에 영향을 미쳐 2D–3D 등록 정확도를 직접적으로 저하시킨다.

저자들은 본 연구를 단순한 체적 재구성 충실도를 넘어 투영 영역 작업(DRR 생성 및 등록)의 특정 요구 사항을 고려하는 보다 총체적인 교정 평가를 향한 단계로 위치시킨다. 제안된 프레임워크는 기하학적 민감도를 특성화하기 위한 통제된 환경을 제공하며, 저자들은 이것이 향후 교정 분석 및 DRR 생성 평가에 유용할 수 있다고 제안한다.

한

저자들은 다음과 같은 한계점을 명시적으로 인정한다:

• 합성 데이터의 특성: 실험은 산란(scatter), 검출기 노이즈, 빔 경화(beam hardening) 또는 비균일한 검출기 응답과 같은 물리적 효과를 모델링하지 않은 CT 유래 메쉬 모델을 사용한다.
• 강체 해부학: 연구는 연부 조직의 상호작용, 변형 또는 환자 움직임을 제외한 강체 척추 해부학을 가정한다.
• 섭동 범위: 분석은 분석적으로 정의된 고유 섭동에 국한되어 있으며, 실제 시스템은 고유, 외적 및 시간적 변화의 복잡한 조합을 포함할 수 있다.
• 등록 범위: 등록 평가는 전체 강도 기반 등록 파이프라인이나 실제 임상 영상의 열화가 아닌, 사전 정의된 랜드마크를 가진 통제된 합성 환경으로 제한된다.

향ol 향후 연구에서는 실제 투시 데이터셋을 사용하여 이러한 발견을 검증하고, 현실적인 검출기 물리학을 통합하며, 변형 가능한 해부학적 시뮬레이션을 평가할 것을 제안한다.