Stability of 0-dimensional persistent homology in enriched and sparsified point clouds
๐ Abstract
**
์ฐ๋ฆฌ๋ ์ ํ ์งํฉ ์ฌ์ด์ ํ๋ถํ(์ ์์ ์ถ๊ฐ), ํฌ์ํ(์์ ์ ๊ฑฐ), ๊ฒฉ์ ์ ๋ ฌ(์ขํ๋ฅผ ๊ณ ์ ๋ ์คํ
ํฌ๊ธฐ๋ก ๋๋์ด ๋ชซ์ ์ทจํจ) ๊ด๊ณ์ ์ํด ์ฐ๊ฒฐ๋ VietorisโRips, ฮฑ, ๊ทธ๋ฆฌ๊ณ ํ๋น
๋ณตํฉ์ฒด ํํฐ๋ง์์ 0์ฐจ์ ์ง์ ๋ํํ๊ณผ ์ฝ๋ค๋ฉ์
1 ๋ํํ์ ๋ํ ๊ฒฝ๊ณ๋ฅผ ์ ์ํ๋ค. ํ๋ถํ๋ ์ค์ฌ๋ถํ (barycentric subdivision)์, ํฌ์ํ๋ ์ต์ ๋ถ๋ฆฌ ๊ฑฐ๋ฆฌ(minimum separating distance)๋ฅผ, ๊ฒฉ์ ์ ๋ ฌ์ ๊ณ ์ ์คํ
ํฌ๊ธฐ๋ก ์ขํ๋ฅผ ๋๋ ๋ค ๋ชซ์ ์ทจํ๋ ๋ฐฉ๋ฒ์ ์ฌ์ฉํ๋ค. ๋๊ท๋ชจยท๋ถ๊ท์น ๋ฐ์ดํฐ์
์ ๋ค๋ฃจ๋ ์ค์ ์์ฉ์ ๋๊ธฐ๋ก ์ผ์ผ๋ฉฐ, ํนํ ํ๊ฒฝ ๋ณ์๋ค์ ์ฐจ์์ผ๋ก ๊ฐ๋ ๊ณ ์ฐจ์ โํ์ดํผ๋ณผ๋ฅจโ์ ํตํด ๊ด์ฐฐ๋ ์ข
์ ์ํ๋ฅผ ์ถ์ ํ๋ ์ํํ์ ์ฌ๋ก๋ฅผ ๋ค๋ฃฌ๋ค. ์ด ํ์ดํผ๋ณผ๋ฅจ์ ๋ถํผยท๋ณผ๋ก์ฑ ๊ฐ์ ๊ธฐํํ์ ํน์ฑ๊ณผ ์ฐ๊ฒฐ์ฑยท๋ํํ ๊ฐ์ ์์์ ํน์ฑ์ ๋์์ ๊ฐ์ผ๋ฉฐ, ์ด๋ ์ข
์ ํ์ฌยท๋ฏธ๋ ์ํ์ ์ฐ๊ด๋๋ค. ์ฐ๋ฆฌ๋ VietorisโRips์ ฮฑ ๋ณตํฉ์ฒด์ ์ง์ ๋ค์ด์ด๊ทธ๋จ ์ฌ์ด์ ์ ํํ ๊ฒฝ๊ณ๋ฅผ ์ ๊ณตํ๊ณ , ํ๋น
๋ณตํฉ์ฒด์ ๋ํ ์ด์ค์ฑ ์ ์ฒด์ฑ์ ์ ์ํ๋ค. ์ ์๋ ๋ฐฉ๋ฒ๋ก (TopoAware)์ C++, Python, R ์ธํฐํ์ด์ค๋ฅผ ํตํด GUDHI ๋ผ์ด๋ธ๋ฌ๋ฆฌ๋ฅผ ๊ธฐ๋ฐ์ผ๋ก ๊ตฌํ๋์ด ๊ณต๊ฐ๋๋ค.
**
๐ก Deep Analysis
**
1. ์ฐ๊ตฌ ๋ฐฐ๊ฒฝ ๋ฐ ๋๊ธฐ
- ๋๊ท๋ชจยท๋ถ๊ท์น ๋ฐ์ดํฐ: ํ๋ ๋ฐ์ดํฐ ๊ณผํ์์ ์ ๊ตฌ๋ฆ์ด ์์ฒยท์๋ง ๊ฐ์ ๋ฌํ๊ณ , ์ธก์ ์ค์ฐจยท๊ฒฐ์ธก์นยท๋ถ๊ท ํ ๋ถํฌ๊ฐ ํํ ๋ฐ์ํ๋ค. ๊ธฐ์กด ์ง์ ๋ํํ(Persistent Homology) ์๊ณ ๋ฆฌ์ฆ์ ์ด๋ฌํ ์ํฉ์์ ๊ณ์ฐ ๋น์ฉ๊ณผ ์์ ์ฑ ๋ฌธ์ ๊ฐ ์ ๊ธฐ๋๋ค.
- ์ํํ์ ํ์ดํผ๋ณผ๋ฅจ: ํ๊ฒฝ ๋ณ์(๊ธฐํ, ํ ์, ํ ์ pH ๋ฑ)๋ฅผ ์ฐจ์์ผ๋ก ํ๋ ๊ณ ์ฐจ์ ๊ณต๊ฐ์์ ์ข ์ ์์ ๊ฐ๋ฅ ์์ญ์ ์ ์ํ๋ค. ์ด ์์ญ์ ์์(์: ์ฐ๊ฒฐ์ฑ ํํธ)๊ณผ ๊ธฐํ(์: ๋ถํผ ๊ฐ์)๋ ๋ฉธ์ข ์ํ์ ์์ธกํ๋ ๋ฐ ์ค์ํ ์ ํธ๊ฐ ๋๋ค.
2. ์ฃผ์ ๊ธฐ์ฌ
| ๊ตฌ๋ถ | ๋ด์ฉ | ์์ |
|---|---|---|
| ์ด๋ก ์ ๊ฒฝ๊ณ | ํ๋ถํยทํฌ์ํยท๊ฒฉ์ ์ ๋ ฌ์ ๋ํด VietorisโRips, ฮฑ, ํ๋น ๋ณตํฉ์ฒด์ 0โ์ฐจ ์ง์ ๋ฐ์ฝ๋์ ์ฝ๋ค๋ฉ์ โ1 ๋ํํ์ ๋ํ ๋ช ์์ ์ํ/ํํ์ ์ ๊ณต | ๋ฐ์ดํฐ ์ ์ฒ๋ฆฌ(์ถ๊ฐยท์ ๊ฑฐยท๊ฒฉ์ํ) ํ์๋ ์์์ ์ ๋ณด๊ฐ ์ผ๋ง๋ ๋ณํ๋์ง ์ ๋ํ |
| ํ๋ถํ ๊ธฐ๋ฒ | ์ค์ฌ๋ถํ (barycentric subdivision) ์ฌ์ฉ โ ๊ธฐ์กด ์ ์ ์๋ก์ด ์ ์ ์ฝ์ ํ๋ฉด์ ๋ณตํฉ์ฒด ๊ตฌ์กฐ๋ฅผ ๋ณด์กด | ๋ฐ์ดํฐ ๋ณด๊ฐ(์: ์๋ฎฌ๋ ์ด์ ๋ฐ์ดํฐ ์ถ๊ฐ) ์ ์์ ์๊ณก ์ต์ํ |
| ํฌ์ํ ๊ธฐ๋ฒ | ์ต์ ๋ถ๋ฆฌ ๊ฑฐ๋ฆฌ(minimum separating distance) ๊ธฐ๋ฐ ์ํ๋ง โ ์ ๊ฐ ๊ฑฐ๋ฆฌ๊ฐ ์ถฉ๋ถํ ํฐ ๊ฒฝ์ฐ์๋ง ์ ๊ฑฐ | ๋ฐ์ดํฐ ์์ถยท์ฐ์ฐ๋ ๊ฐ์์ ๋์์ ์์ ์์ค์ ์ ์ด |
| ๊ฒฉ์ ์ ๋ ฌ ๊ธฐ๋ฒ | ๊ณ ์ ์คํ ํฌ๊ธฐ๋ก ์ขํ๋ฅผ ์์ํํ๊ณ ๋ชซ์ ์ทจํจ โ ๊ฒฉ์ ์์์ ๋๋ฑ ํด๋์คํ | ๋ฉ๋ชจ๋ฆฌ ํจ์จ ํฅ์ ๋ฐ GPU/๋ถ์ฐ ํ๊ฒฝ์์์ ๊ตฌํ ์ฉ์ด |
| ์ด์ค์ฑ ์ ์ฒด์ฑ | ํ๋น ๋ณตํฉ์ฒด์ ๋ํด โdualโ ๊ด๊ณ๋ฅผ ์์ํ โ ฮฑ ๋ณตํฉ์ฒด์์ ๋น๊ต ๊ฐ๋ฅ์ฑ ํ๋ณด | ๋ค์ํ ๋ณตํฉ์ฒด ์ ํ ์ ์ผ๊ด๋ ์์ ํด์ ์ ๊ณต |
| ์ํํธ์จ์ด | TopoAware ํจํค์ง (C++, Python, R) ๊ณต๊ฐ, GUDHI ๊ธฐ๋ฐ | ์ค๋ฌดยท์ฐ๊ตฌ์๋ค์ด ๋ฐ๋ก ์ ์ฉ ๊ฐ๋ฅํ๋๋ก ๊ตฌํ ์ฅ๋ฒฝ ๋ฎ์ถค |
3. ๋ฐฉ๋ฒ๋ก ์์ธ
- ์ ๊ตฌ๋ฆ ๋ชจ๋ธ๋ง
- ์ ๋ ฅ: ์ ํ ์งํฉ (X \subset \mathbb{R}^d).
- ํ๋ถํ: ๊ฐ ๋จ์์ฒด์ ์ค์ฌ์ ์ถ๊ฐ โ ์๋ก์ด ์ ์งํฉ (X^{+}).
- ํฌ์ํ: ๊ฑฐ๋ฆฌ ํจ์ (d_{\min}) ์ ์ ํ, (d($x_i$, $x_j$) < d_{\min}) ์ธ ์ ์ ํ๋๋ง ๋จ๊น โ (X^{-}).
- ๊ฒฉ์ ์ ๋ ฌ: ์คํ (\delta >0) ์ ํ, (\lfloor x/\delta \rfloor) ๋ก ๋งคํ โ ๊ฒฉ์ ํด๋์ค (
๐ Full Content
์ฐ๋ฆฌ๋ ์ ํ ์งํฉ๋ค ์ฌ์ด์์ ํ๋ถํ(์ ์์ ์ถ๊ฐ), ํฌ์ํ(์์ ์ ๊ฑฐ), ๊ฒฉ์ ์ ๋ ฌ(์์๋ฅผ ๊ท ์ผํ๊ฒ ์ด์ฐํ)๊ณผ ๊ฐ์ ์ฐ์ฐ์ ์ํด ๊ด๋ จ๋ VietorisโRips, ์ํ, ๊ทธ๋ฆฌ๊ณ ํ๋น์ปฌ ๋ณตํฉ์ฒด ํํฐ๋ง์ ์ฐจ์โฏ0 ์ง์ ๋ํ๋ก (persistent homology) ๋ฐ ์ฝ์ฐจ์โฏ1 ํธ๋ชฐ๋ก์ง์ ๋ํ ๊ฒฝ๊ณ๋ฅผ ์ ์ํ๋ค. ํ๋ถํ์ ๋ํด์๋ ์ค์ฌ๋ถํ (barycentric subdivision)์ ์ฌ์ฉํ๊ณ , ํฌ์ํ์ ๋ํด์๋ ์ต์ ๋ถ๋ฆฌ ๊ฑฐ๋ฆฌ(minimum separating distance)๋ฅผ ์ด์ฉํ๋ฉฐ, ๊ฒฉ์ ์ ๋ ฌ์ ๋ํด์๋ ๊ฐ ์ขํ ๊ฐ์ ๊ณ ์ ๋ ์คํ ํฌ๊ธฐ๋ก ๋๋ ๋ค ๋ชซ์ ์ทจํ๋ ๋ฐฉ๋ฒ์ ์ ์ฉํ๋ค.
์ด๋ฌํ ๋ฐฉ๋ฒ๋ก ์ ํฌ๊ธฐ๊ฐ ํฌ๊ณ ๋ถ๊ท์นํ ๋ฐ์ดํฐ์ ์ ๋ค๋ฃจ๋ ์ค์ ์์ฉ์ ๋๊ธฐ๋ฅผ ๋ถ์ฌ๋ฐ์์ผ๋ฉฐ, ์ง์ ๋ํ๋ก ์ ์ด๋ฌํ ๋ฐ์ดํฐ์ ๋ณด๋ค ํจ๊ณผ์ ์ผ๋ก ์ ์ฉํ๊ธฐ ์ํ ์ฐ๊ตฌ ๊ฐ๋ฐ์ ์ผํ์ด๋ค. ํนํ, ์ฐ๋ฆฌ๋ ์ํํ์ ์์ฉ์ ๊ณ ๋ คํ๋ค. ์ฌ๊ธฐ์ ๊ด์ฐฐ๋ ์ข ์ ์ํ๋ ํ๊ฒฝ ๋ณ์๋ค์ ์ฐจ์์ผ๋ก ๊ฐ๋ ๊ณ ์ฐจ์ ๊ณต๊ฐ์ ํตํด ์ถ๋ก ๋๋ค. ์ด๋ฌํ โํ์ดํผ๋ณผ๋ฅจโ(hypervolume)์ ๋ถํผ์ ๋ณผ๋ก์ฑ ๊ฐ์ ๊ธฐํํ์ ํน์ฑ(geometry)๊ณผ ์ฐ๊ฒฐ์ฑ ๋ฐ ํธ๋ชฐ๋ก์ง์ ๊ฐ์ ์์ํ์ ํน์ฑ(topology)์ ๋์์ ๊ฐ์ง๊ณ ์์ผ๋ฉฐ, ์ด๋ค ํน์ฑ์ ์ข ์ ํ์ฌ ์ํ์ ์ ์ฌ์ ์ธ ๋ฏธ๋ ์ํ์ ์ฐ๊ด์ด ์์์ด ์๋ ค์ ธ ์๋ค.
์ฐ๋ฆฌ๋ ๊ธฐ์กด์ ํต๊ณ์ ๋ฐฉ๋ฒ์ ๋ณด์ํ๋ ํํ๋ก, ์์ํ์ ๋ณด์ฅ์ ์ ๊ณตํ๋ ํ์ดํผ๋ณผ๋ฅจ ๋ถ์ ์ ๊ทผ๋ฒ์ ์ ์ํ๋ค. ๊ตฌ์ฒด์ ์ผ๋ก, VietorisโRips ๋ณตํฉ์ฒด์ ์ํ ๋ณตํฉ์ฒด์ ์ง์์ฑ ๋ค์ด์ด๊ทธ๋จ(persistence diagram) ์ฌ์ด์ ์ ํํ ๊ฒฝ๊ณ๋ฅผ ์ ๊ณตํ๊ณ , ํ๋น์ปฌ ๋ณตํฉ์ฒด์ ๋ํด์๋ ์ด์ค์ฑ(duality) ์ ์ฒด์ฑ์ ์ ์ํ๋ค.
์ฐ๋ฆฌ์ ๋ฐฉ๋ฒ๋ก ์ TopoAware๋ผ๋ ์ด๋ฆ์ผ๋ก ๊ตฌํ๋์์ผ๋ฉฐ, C++, Python, R ์ธํฐํ์ด์ค๋ฅผ ํตํด ์ ๊ณต๋๋ค. ๊ตฌํ์ GUDHI ๋ผ์ด๋ธ๋ฌ๋ฆฌ๋ฅผ ๊ธฐ๋ฐ์ผ๋ก ํ์ฌ, ์ฌ์ฉ์๊ฐ ๋๊ท๋ชจ ๋ฐ์ดํฐ์ ๋ํด ํจ์จ์ ์ผ๋ก ์ง์ ๋ํ๋ก ์ ๊ณ์ฐํ๊ณ , ์์์ ์ ์ํ ์ด๋ก ์ ๊ฒฝ๊ณ์ ์ด์ค์ฑ ๊ฒฐ๊ณผ๋ฅผ ์คํ์ ์ผ๋ก ๊ฒ์ฆํ ์ ์๋๋ก ์ค๊ณ๋์๋ค. TopoAware๋ ๋ฐ์ดํฐ ์ ์ฒ๋ฆฌ ๋จ๊ณ์์ ํ๋ถํ, ํฌ์ํ, ๊ฒฉ์ ์ ๋ ฌ์ ์ ํ์ ์ผ๋ก ์ ์ฉํ ์ ์๋ ํ์ดํ๋ผ์ธ์ ์ ๊ณตํ๋ฉฐ, ๊ฐ ๋จ๊ณ์์ ์ฌ์ฉ์๊ฐ ์ค์ ํ ์ ์๋ ํ๋ผ๋ฏธํฐ(์: ์ค์ฌ๋ถํ ์ ๋ฐ๋ณต ํ์, ์ต์ ๋ถ๋ฆฌ ๊ฑฐ๋ฆฌ์ ์๊ณ๊ฐ, ๊ฒฉ์ ์คํ ํฌ๊ธฐ ๋ฑ)๋ฅผ ํตํด ๋ถ์ ๊ฒฐ๊ณผ์ ์ ๋ฐ๋์ ๊ณ์ฐ ๋น์ฉ ์ฌ์ด์ ๊ท ํ์ ์กฐ์ ํ ์ ์๋ค.
๋ํ, ๊ฒฐ๊ณผ๋ก ์ป์ด์ง๋ ์ง์์ฑ ๋ค์ด์ด๊ทธ๋จ์ ๊ธฐ์กด์ ํต๊ณ์ ์์ฝ๊ฐ(์: ํ๊ท , ๋ถ์ฐ)๊ณผ ํจ๊ป ์๊ฐํ ๋ฐ ํด๋ฌ์คํฐ๋ง์ ํ์ฉ๋ ์ ์์ด, ์ฐ๊ตฌ์๊ฐ ์ข ์ ์์์ง ๋ณํ๋ ํ๊ฒฝ ์คํธ๋ ์ค ์์ธ์ ์ํฅ์ ์ ๋์ ์ผ๋ก ํ๊ฐํ๋ ๋ฐ ์๋ก์ด ํต์ฐฐ์ ์ ๊ณตํ๋ค. ๊ถ๊ทน์ ์ผ๋ก ๋ณธ ์ฐ๊ตฌ๋ ๋ณตํฉ์ ์ธ ๋ฐ์ดํฐ ๊ตฌ์กฐ์ ์์ํ์ ๋ถ์์ ์ฐ๊ฒฐํจ์ผ๋ก์จ, ๋๊ท๋ชจ ๋ถ๊ท์น ๋ฐ์ดํฐ์ ๋ํ ์ง์ ๋ํ๋ก ์ ์ ์ฉ ๊ฐ๋ฅ์ฑ์ ํฌ๊ฒ ํ์ฅํ๊ณ , ์ํํ์ ๋ฐ์ดํฐ ํด์๋ฟ๋ง ์๋๋ผ ๋ค์ํ ๋ถ์ผ์ ๊ณ ์ฐจ์ ๋ฐ์ดํฐ ๋ถ์์ ์ ์ฉํ ๋๊ตฌ๋ฅผ ์ ๊ณตํ๋ค.
ํฅํ ์ฐ๊ตฌ์์๋ TopoAware์ ํ์ฅ์ฑ์ ๊ฒ์ฆํ๊ธฐ ์ํด ๋ค์ํ ๋๋ฉ์ธ(์: ์๋ฃ ์์, ๋ฌผ๋ฆฌํ ์๋ฎฌ๋ ์ด์ , ์ฌํ ๋คํธ์ํฌ)์์์ ์ ์ฉ ์ฌ๋ก๋ฅผ ์กฐ์ฌํ ๊ณํ์ด๋ฉฐ, ํนํ ๊ณ ์ฐจ์ ๋ฐ์ดํฐ์ ๋ ธ์ด์ฆ์ ๋ํ ๊ฐ์ธ์ฑ์ ํ๊ฐํ๊ธฐ ์ํด ํฉ์ฑ ๋ฐ์ดํฐ์ ์ค์ ์ธก์ ๋ฐ์ดํฐ๋ฅผ ๋์์ ํ์ฉํ ์์ ์ด๋ค. ๋ํ, GUDHI ๊ธฐ๋ฐ ๊ตฌํ์ ๋ณ๋ ฌํ์ GPU ๊ฐ์์ ํตํด ์ด๋๊ท๋ชจ ๋ฐ์ดํฐ์ (์์ต ๊ฐ์ ํฌ์ธํธ)์์๋ ์ค์๊ฐ์ ๊ฐ๊น์ด ๋ถ์์ด ๊ฐ๋ฅํ๋๋ก ์ต์ ํ ์์ ์ ์งํ ์ค์ด๋ค. ์ด๋ฌํ ๊ธฐ์ ์ ์ง๋ณด๋ ์ง์ ๋ํ๋ก ์ ์ค์ฉ์ ์ธ ๋ฐ์ดํฐ ๊ณผํ ๋๊ตฌ๋ก ์๋ฆฌ๋งค๊นํ๊ฒ ํ ๊ฒ์ด๋ฉฐ, ๊ถ๊ทน์ ์ผ๋ก๋ ๋ฐ์ดํฐ ๊ธฐ๋ฐ ์์ฌ๊ฒฐ์ ๊ณผ์ ์์ ์์ํ์ ์ธ์ฌ์ดํธ๋ฅผ ์ ๋์ ์ผ๋ก ํ์ฉํ ์ ์๋ ๊ธฐ๋ฐ์ ๋ง๋ จํ ๊ฒ์ด๋ค.