DUNEuro 바이오전기자기학을 위한 고성능 전방 모델링 툴박스

초록

**
DUNEuro는 DUNE 프레임워크 위에 구축된 오픈소스 C++ 툴박스로, EEG·MEG 전방 문제를 해결하기 위해 적합(fitted) 및 비적합(unfitted) 유한요소법을 제공한다. 다양한 전류원 모델과 Python·Matlab 인터페이스를 지원하며, 6구획 현실 두뇌 모델을 이용한 체감각 유발 전위 분석 예제로 실용성을 입증한다.

**

상세 분석

**
본 논문은 뇌전도(EEG)와 뇌자기(MEG) 전방 문제를 고해상도 현실 두뇌 해부학에 적용하기 위한 수치적 기반을 제시한다. 핵심은 DUNE(Distributed and Unified Numerics Environment) 라이브러리를 활용해 모듈화된 구조를 구현한 점이다. DUNEuro는 ‘적합(fitted)’과 ‘비적합(unfitted)’ 두 가지 FEM 접근법을 모두 제공한다. 적합 방식은 전두부·뇌실·두개골·피질 등 실제 조직 경계에 격자를 정확히 맞추어 CG‑FEM(연속 Galerkin)과 DG‑FEM(불연속 Galerkin)을 구현한다. 비적합 방식은 레벨셋(level‑set)으로 조직 경계를 정의하고, CutFEM·UDG(불연속 Galerkin)와 같은 기법으로 격자와 경계의 불일치를 약하게 처리한다. 이러한 설계는 복잡한 6구획(뇌척수액, 뇌실, 회백질, 백질, 두개골, 피부) 모델을 손쉽게 다룰 수 있게 한다.

전류원 모델링 측면에서는 부분적분, St. Venant, 전위 차감(full/subtraction), Whitney 등 기존에 제안된 여러 방법을 소스 모델 팩토리(SourceModelFactory) 형태로 캡슐화한다. 이는 전류쌍극자와 같은 특이점 소스를 FEM 격자에 삽입할 때 발생하는 수치적 발산을 효과적으로 억제한다. 또한 전방 해를 직접 구하지 않고 전극·센서 위치에서만 전위·자기장을 얻는 전이 행렬(transfer matrix) 개념을 도입해, 다수의 소스에 대해 선형 연산만으로 빠르게 결과를 산출한다.

소프트웨어 구조는 MEEGDriverInterface라는 고수준 인터페이스를 중심으로, FittedMEEGDriver와 UnfittedMEEGDriver가 각각 적합·비적합 방식을 구현한다. Solver와 SourceModelFactory는 템플릿 파라미터화되어 새로운 선형 솔버(예: 직접법, 멀티그리드)나 새로운 소스 모델을 추가하기 용이하도록 설계되었다. DUNE‑uggrid, dune‑subgrid, dune‑pdelab, dune‑istl 등 기존 DUNE 모듈을 재사용함으로써 메모리 효율과 코드 재사용성을 극대화하였다.

실험에서는 6구획 현실 두뇌 모델을 이용해 체감각 유발 전위(SSEP)를 분석하였다. 전극 위치에 대한 전위와 자기장을 계산한 뒤, 전위 역전파(source localization) 알고리즘을 적용해 실제 뇌 활동 위치를 복원하였다. 결과는 기존 SimBio·FieldTrip 구현과 비교해 정확도와 계산 속도 모두에서 경쟁력을 보였으며, Python·Matlab 바인딩을 통해 기존 파이프라인에 손쉽게 통합할 수 있음을 시연하였다.

전반적으로 DUNEuro는 고급 FEM 기법을 모듈화·표준화함으로써 연구자와 임상 엔지니어가 복잡한 전두부 모델을 손쉽게 다루고, 다양한 소스 모델과 해석 파이프라인에 유연하게 연결할 수 있는 플랫폼을 제공한다. 향후 확장성 있는 설계 덕분에 전기·자기 자극(E‑/TMS), 전위 지도화, 뇌-컴퓨터 인터페이스 등 다양한 바이오전기자기 응용 분야로의 파생이 기대된다.

**