보존 동형성을 이용한 보로노이 보간 중단 기준

본 연구는 보로노이 보간(Voronoi interpolation) 과정에서 전역 토폴로지 정보를 고려하지 못하는 기존 방법들의 한계를 극복하고자, 지속 동형성(persistent homology)을 중단 기준(stopping‑criterion)으로 도입한다. 논문은 먼저 보로노이 테셀레이션이 듈러네 삼각분할(Delaunay triangulation)과 이중 관계를 가지며, 이 듈러네 복합체를 이용하면 점 집합 X에 대한 필터링된 복합체를 효율적으로 구성할 수 있음을 설명한다. 1. **수학적 배경** - **단순 복합체와 필터링**: 점 집합 X⊂ℝᵈ에서 정의되는 Čech 복합체, Vietoris‑Rips 복합체, 듈러네 복합체, Witness 복합체를 차례로 소개하고, 각각의 정의와 포함 관계를 정리한다. 특히 Čech 복합체는 Nerve 보조정리를 통해 듈러네 복합체와 동형동치임을 보이며, 이는 복합체의 위상적 특성을 정확히 반영한다는 점을 강조한다. - **지속 동형성 이론**: 필터링 파라미터 r에 따라 복합체 K_r이 증가하는 과정을 통해 동형성 군 H_k(K_r)와 그 변화를 추적한다. 영속성 다이어그램(Persistence diagram)은 (birth, death) 쌍으로 나타내며, 각 차원 k에 대해 β_k(=rank H_k)와 멀티플리시티 μ_{n,r}^k를 정의한다. 2. **알고리즘 설계** - **보간 단계**: 초기 점 집합 X에 새로운 점을 하나씩 삽입하면서 듈러네 복합체 Del(X, r) 를 재구성한다. 각 단계마다 복합체에 대한 지속 동형성을 계산하고, 두 연속 단계의 영속성 다이어그램 P₁, P₂ 사이의 Bottleneck 거리 d_B와 p‑Wasserstein 거리 d_W를 구한다. - **중단 기준**: 두 거리의 L₂ 노름 ‖(d_B, d_W)‖₂가 사전에 설정한 임계값 τ를 초과하면, 더 이상의 점 삽입이 위상적 변화를 크게 일으키지 않는다고 판단하고 알고리즘을 종료한다. τ는 실험적 경험법칙에 의해 결정되며, 데이터의 노이즈 수준과 목표 정확도에 따라 조정 가능하다. 3. **이론적 정당성** - **Nerve Lemma 적용**: 듈러네 복합체가 Čech 복합체의 Nerve와 동형동치임을 이용해, 복합체를 통해 얻은 영속성 다이어그램이 실제 점 집합의 볼록합(Union of balls) 위상의 변화를 정확히 반영함을 증명한다. - **복합체 선택과 복잡도**: Čech 복합체는 계산 비용이 O(2ⁿ)으로 비현실적이지만, 듈러네 복합체와 Vietoris‑Rips 복합체는 O(n²) 혹은 O(n log n) 수준으로 구현 가능함을 제시한다. 또한 Witness 복합체를 이용하면 대규모 데이터에 대한 스케일링이 가능함을 논한다. 4. **실험 및 결과** - **데이터셋**: 서명(필기) 이미지에서 추출한 점 구름을 사용해 H₀, H₁, H₂ 차원의 영속성을 분석한다. 서명 획의 교차점에서는 H₁(터널)과 H₂(볼록체)가 나타나며, 이는 기존 보간 방법이 놓치기 쉬운 고차 위상 특성이다. - **거리 변화 관찰**: 초기 몇 번의 삽입 단계에서 d_B와 d_W가 급격히 상승하고, 이후 단계에서는 완만하게 감소하거나 거의 변하지 않는다. 이는 알고리즘이 “위상적 안정화” 시점을 자동으로 탐지함을 의미한다. - **중단 효과**: 제안된 τ=0.05(예시) 기준으로 평균 12~15번의 삽입 후 중단되며, 최종 보간 결과는 원본 점 집합과 시각적으로 거의 구분되지 않는다. 또한, 과잉 삽입으로 인한 계산 비용과 메모리 사용량이 크게 절감된다. 5. **결론 및 향후 연구** - 지속 동형성을 중단 기준으로 활용함으로써 보로노이 보간이 전역 위상 구조를 보존하면서도 불필요한 과잉 샘플링을 방지한다는 점을 입증하였다. 향후 연구에서는 비유클리드 거리, 고차원 데이터, 그리고 실시간 스트리밍 상황에서의 동적 τ 조정 메커니즘을 탐색할 예정이다.