동적 스키닝 기반 관상동맥 변형 모델링: 실시간 외과 시뮬레이션을 위한 새로운 프레임워크

초록

본 연구는 관상동맥을 해부학적으로 정확하게 변형시키기 위해 이중조화 에너지 최소화를 이용한 스키닝 가중치 계산과 3차원 테트라헤드론 메쉬를 결합한 동적 모델링 프레임워크를 제안한다. 시간 샘플링·보간을 통해 심장 주기 전체에 걸친 연속 변형을 구현하고, 부피 보존 및 기계적 제약을 적용한다. 5명의 임상 데이터에 대한 정량적 평가에서 평균 Hausdorff 거리 4.96 mm, 평균 표면 거리 1.78 mm 등 실시간 성능과 충분한 기하학적 정확성을 입증하였다.

상세 분석

이 논문은 관상동맥 시뮬레이션에서 가장 큰 난제 중 하나인 고해상도 해부학적 변형을 실시간으로 구현하는 방법을 제시한다. 핵심 아이디어는 스키닝(skinning) 기법을 의료 영상에 적용하는 것으로, 기존 그래픽스 분야에서 사용되는 선형 블렌드 스키닝(LBS)이나 듀얼 쿼터니언 스키닝(DQS)과 달리, 여기서는 해부학적 구조와 물리적 제약을 동시에 만족하도록 가중치를 설계한다. 가중치 계산은 biharmonic energy minimization을 기반으로 하며, 이는 스무스하면서도 변형에 대한 국소적 강성을 보존한다. 구체적으로, 관상동맥 중심선과 분기점을 제어 포인트로 설정하고, 전체 혈관을 tetrahedral mesh로 볼륨화한다. 그런 다음 각 정점에 대해 제어 포인트와의 거리와 곡률 정보를 포함한 에너지 함수를 정의하고, 라그랑주 승수를 이용해 부피 보존 제약을 추가한다. 최적화는 선형 시스템을 풀어 빠르게 수렴하도록 설계돼, 실시간 인터랙션에 적합하다.

시간적 측면에서는 심장 주기를 여러 정밀 프레임(예: 30 Hz)으로 샘플링하고, 각 프레임에 대해 위에서 구한 스키닝 가중치를 적용한다. 두 인접 프레임 사이의 변형은 선형 보간이 아니라, 가중치 자체를 시간에 따라 스무스하게 변화시키는 방식으로 수행한다. 이를 통해 급격한 수축·이완 구간에서도 비물리적인 진동이나 겹침을 최소화한다. 또한, 가이드와이어와 같은 도구와 혈관 벽 사이의 충돌 검출을 위해 가중치 기반 변형된 메쉬에 대한 BVH(Bounding Volume Hierarchy)를 실시간으로 업데이트한다.

검증에서는 5명의 환자로부터 CT 또는 MRI 데이터를 이용해 4개의 주요 심장 단계(수축기, 이완기 등)를 각각 세분화하고, 이를 ground truth로 삼았다. 제안된 프레임워크는 각 단계 사이를 보간한 결과와 실제 segmentation 사이의 Hausdorff 거리 평균 4.96 mm(표준편차 1.78 mm), 평균 표면 거리 1.78 mm(표준편차 0.75 mm)를 기록했다. Branch Completeness Ratio가 1.82 ± 0.46으로 높은 편이며, Branch Continuity Score는 0.84 ± 0.06으로 거의 완전한 연속성을 유지한다. 이러한 정량적 결과는 기존 스프링-매시 기반 방법보다 우수함을 시사한다.

한계점으로는 현재 가중치 최적화가 정적 해부학적 모델에 국한되어 있어, 병변(플라그)이나 급격한 혈관 경화와 같은 비선형 변형을 직접 모델링하기 어렵다는 점을 들 수 있다. 또한, 5명의 제한된 데이터셋으로 일반화 가능성을 완전히 검증하지 못했다. 향후 연구에서는 비선형 물성 모델을 가중치에 통합하고, 더 다양한 환자군을 대상으로 임상 시험을 진행할 계획이다.