[논문 리뷰] Approximate Cech Complexes in Low and High Dimensions

Michael Kerber, R. Sharathkumar|arXiv (Cornell University)|2013. 07. 11.

Topological and Geometric Data Analysis참고 문헌 17인용 수 3

한 줄 요약

이 논문은 위상적 데이터 분석에서 Čech 필터링에 대한 두 가지 새로운 근사 방법을 제안한다: 고정된 차원에서 잘 분리된 단체 분해(WSSD)를 사용한 선형 크기 근사, 그리고 임의의 환경 차원에 관계없이 $(1+\varepsilon)$-근사 Čech 필터링을 가능하게 하는 Rips 보조정리의 일반화. 주요 기여는 환경 차원에 의존하지 않는 근사 프레임워크로, 스케일러비리티를 향상시키면서도 위상적 정확성을 유지한다.

ABSTRACT

Čech complexes reveal valuable topological information about point sets at a certain scale in arbitrary dimensions, but the sheer size of these complexes limits their practical impact. While recent work introduced approximation techniques for filtrations of (Vietoris-)Rips complexes, a coarser version of Čech complexes, we propose the approximation of Čech filtrations directly. For fixed dimensional point set $S$, we present an approximation of the Čech filtration of $S$ by a sequence of complexes of size linear in the number of points. We generalize well-separated pair decompositions (WSPD) to well-separated simplicial decomposition (WSSD) in which every simplex defined on $S$ is covered by some element of WSSD. We give an efficient algorithm to compute a linear-sized WSSD in fixed dimensional spaces. Using a WSSD, we then present a linear-sized approximation of the filtration of Čech complex of $S$. We also present a generalization of the known fact that the Rips complex approximates the Čech complex by a factor of $\sqrt{2}$. We define a class of complexes that interpolate between Čech and Rips complexes and that, given any parameter $\eps > 0$, approximate the Čech complex by a factor $(1+\eps)$. Our complex can be represented by roughly $O(n^{\lceil 1/2\eps ceil})$ simplices without any hidden dependence on the ambient dimension of the point set. Our results are based on an interesting link between Čech complex and coresets for minimum enclosing ball of high-dimensional point sets. As a consequence of our analysis, we show improved bounds on coresets that approximate the radius of the minimum enclosing ball.

연구 동기 및 목표

• 정확한 Čech 복합체의 높은 계산 비용을 해결하기 위해, $k$-스켈레톤에서 $O(n^k)$개의 단체를 가질 수 있는 문제점을 해결한다.
• 상수 인자 $(1+\varepsilon)$ 이내의 위상 정확성을 유지하면서도 확장 가능한 Čech 필터링 근사법을 개발한다.
• 고정된 차원의 점 집합에 대해 잘 분리된 단체 분해(WSSD)를 통해 잘 분리된 쌍 분해(WSPD)를 고차원 구조로 일반화한다.
• 완전 복합체를 이용한 차원에 독립적인 Čech 복합체 근사법을 수립하여 환경 차원에 대한 지수적 의존성을 피한다.
• MEB-코어셋과 Čech 복합체 근사 간의 연관성을 활용하여 최소 둘레 구의 코어셋 경계를 향상시킨다.

제안 방법

• 잘 분리된 단체 분해(WSSD)를 도입하여 WSPD의 일반화로, 점 집합을 $O(n/\varepsilon^d)$개의 클러스터로 분해하여 임의의 단체가 클러스터의 튜플에 의해 커버되도록 한다.
• WSSD를 사용해 각 척도에서 필수적인 위상적 구조만을 표현함으로써 Čech 필터링의 선형 크기 근사 복합체를 구성한다.
• $k$-완성 복합체 $\mathcal{M}_k(K)$를 정의하여 $K$와 같은 $k$-스켈레톤을 가지는 최대 복합체로, 표현 크기와 근사 품질 사이의 트레이드오프를 가능하게 한다.
• $\delta = \lceil 1/(2\varepsilon + \varepsilon^2) + 1 \rceil$일 때 $\mathcal{M}_{\delta-1}(\mathcal{C}_\alpha) \subseteq \mathcal{C}_{(1+\varepsilon)\alpha}$임을 증명하여 완성 복합체와 $(1+\varepsilon)$-근사 간의 연결 고리를 확립한다.
• $\varepsilon$-코어셋의 크기가 $\lceil 1/\varepsilon \rceil$ 이하인 최소 둘레 구에 대해, 근사에서 반경 확장을 유한하게 제한하기 위해 이를 활용한다.
• Helly의 정리와 코어셋 성질을 이용하여, 단체의 경계 정점들이 작은 코어셋을 가지면 그 단체가 확장된 Čech 복합체에 포함됨을 보인다.

실험 결과

연구 질문

• RQ1 고정된 차원의 점 집합에 대해, 점의 수에 대해 선형 크기의 복합체로 Čech 필터링을 근사할 수 있는가?
• RQ2 환경 차원에 의존하지 않고 $(1+\varepsilon)$-근사 Čech 복합체를 가능하게 하는 Rips 보조정리의 일반화가 존재하는가?
• RQ3 잘 분리된 단체 분해(WSSD)는 효율적으로 계산되고 Čech 필터링의 확장 가능한 근사에 사용될 수 있는가?
• RQ4 최소 둘레 구 반경의 $\varepsilon$-근사 보장을 보장하는 최소 코어셋 크기는 얼마이며, 이는 위상 근사와 어떻게 관련되는가?
• RQ5 단체 스켈레톤에 기반한 완성 복합체는 확장 가능하고 차원에 독립적인 Čech 복합체 근사법을 제공할 수 있는가?

주요 결과

• 고정된 차원의 점 집합에 대해, WSSD를 사용하여 $O(n)$개의 단체로 제한되는 선형 크기 근사 복합체를 구성한다.
• WSSD 구축은 $n$에 대해 다항식 시간 내에 수행되며, 고정된 차원에 대해서만 지수적 의존성을 가지며 환경 차원에 의존하지 않는다.
• 논문은 $\delta$-완성 복합체가 스케일 $\alpha$에서 $\delta = \lceil 1/(2\varepsilon + \varepsilon^2) + 1 \rceil$일 때 $(1+\varepsilon)\alpha$ 스케일의 Čech 복합체에 포함됨을 증명하여 $(1+\varepsilon)$-근사가 성립함을 보였다.
• 완성 복합체의 크기는 환경 차원에 관계없이 $O(n^{\lceil 1/2\varepsilon \rceil})$이며, 이는 $(1+\varepsilon)$-근접도를 달성하기 위해 최적임이 입증되었다.
• 근사 필터링과 정확한 Čech 필터링의 영속성 다이어그램은 최대 $(1+\varepsilon)$의 요소로 다를 뿐이며, 이는 위상적 정확성을 보장한다.
• 개선된 코어셋 경계를 도출: 최소 둘레 구의 $\varepsilon$-코어셋으로 크기가 $\lceil 1/\varepsilon \rceil$ 이하이면 충분하며, 이 경계는 점의 수와 환경 차원에 독립적이다.