비정상 심근 변형 추정을 위한 일관된 모델

초록

본 논문은 심근경색 등 비정상적인 심근 상태에서의 변형을 정확히 추정하기 위해, 심내·심외막 윤곽점을 중심점 기준으로 정렬하고 변위 벡터를 추출한다. 추출된 변위에 기반한 경계조건을 적용한 타원형 방정식을 유한요소법(FEM)으로 풀어 변형장과 응변 지도를 얻는다. ANSYS로 만든 비균질 링 모델을 이용해 알고리즘을 검증하였다.

상세 분석

이 연구는 기존의 전역적인 변형 모델이 심근경색 부위와 같이 물성비가 급격히 변하는 상황에서 한계를 보이는 점을 정확히 지적한다. 저자는 먼저 심내막(endocardial)과 심외막(epicardial) 윤곽을 이미지에서 추출한 뒤, 내막의 무게중심을 기준으로 각 점을 각도 순서대로 정렬한다는 ‘점‑별(point‑wise)’ 접근법을 채택한다. 이 과정은 복잡한 3차원 형태를 2차원 평면에 투사하면서도 각 점의 상대적 위치 정보를 보존하므로, 변위 벡터를 정확히 계산할 수 있다. 변위 벡터는 시간에 따라 변화하는 내부 압력에 의해 발생하는 실제 물리적 변형을 반영하도록 설계되었으며, 이를 경계조건으로 사용해 타원형(elliptic) 편미분 방정식, 즉 선형 탄성학의 기본 방정식을 정의한다.

경계조건 적용 후, 저자는 유한요소법(FEM)을 통해 전체 심근 영역에 대한 변형장을 수치적으로 해석한다. FEM 구현은 MATLAB 환경에서 이루어졌으며, 요소(mesh) 설계 시 비균질성을 반영하기 위해 물성 파라미터를 지역별로 다르게 할당한다. 특히, 경색 부위는 탄성계수가 현저히 낮은 재료로 모델링함으로써 실제 병변의 ‘연성(softness)’을 재현한다. 변형장이 구해지면, 변형 텐서로부터 1차원 응변(strain) 지도를 추출한다. 이 응변 지도는 정상 부위와 병변 부위 간의 강성 차이를 시각적으로 드러내어, 임상적 진단 및 치료 계획 수립에 직접 활용될 수 있다.

알고리즘 검증을 위해 저자는 ANSYS를 이용해 비균질 원형 링 모델을 구축하였다. 이 모델은 내부 압력이 시간에 따라 선형적으로 증가하는 상황을 시뮬레이션하며, 실제 심실 수축을 근사한다. 링의 일부 구역에만 낮은 탄성계수를 부여해 ‘인프라크트(infarcted)’ 영역을 재현하고, MATLAB 기반 변형 추정 알고리즘을 동일한 입력 조건에 적용한다. 결과적으로 두 시뮬레이션 간 변위 및 응변 분포가 높은 일치도를 보였으며, 이는 제안된 방법이 비균질성(heterogeneity)을 효과적으로 포착함을 증명한다.

기술적 강점으로는 (1) 점‑별 정렬을 통한 변위 추출의 정밀성, (2) 경계조건을 물리적 변위에 직접 연결함으로써 모델의 물리적 일관성 확보, (3) FEM을 활용한 연속적인 변형장 해석으로 전역적인 응변 지도를 제공한다는 점을 들 수 있다. 반면, 제한점으로는 2차원 평면 가정이 실제 3차원 심실 형태를 완전히 대변하지 못한다는 점, 그리고 물성 파라미터를 사전에 정확히 측정하기 어려워 모델 민감도가 높아질 수 있다는 점이 있다. 향후 연구에서는 3D 이미지(예: MRI, CT) 기반의 점‑별 정렬을 확장하고, 환자별 조직강도 측정을 위한 초음파 탄성성 영상(Elastography)과 연계함으로써 모델의 임상 적용성을 크게 향상시킬 수 있을 것이다.