Continuous Dynamic Time Warping (CDTW) measures the similarity of polygonal curves robustly to outliers and to sampling rates, but the design and analysis of CDTW algorithms face multiple challenges. We show that CDTW cannot be computed exactly under the Euclidean 2-norm using only algebraic operations, and we give an exact algorithm for CDTW under norms approximating the 2-norm. The latter result relies on technical fundamentals that we establish, and which generalise to any norm and to related measures such as the partial Fréchet similarity.
논문 ID : 2511.20420제목 : Fundamentals of Computing Continuous Dynamic Time Warping in 2D under Different Norms저자 : Kevin Buchin (TU Dortmund), Maike Buchin (Ruhr University Bochum), Jan Erik Swiadek (Ruhr University Bochum), Sampson Wong (University of Copenhagen)분류 : cs.CG (계산 기하학)발표 시간/학회 : WALCOM 2026 (사전인쇄본, 2025년 11월 제출)논문 링크 : https://arxiv.org/abs/2511.20420 연속 동적 시간 정렬(CDTW)은 다각형 곡선의 유사성을 견고하게 측정할 수 있으며, 이상치 및 샘플링 속도에 대한 견고성을 가지고 있습니다. 그러나 CDTW 알고리즘의 설계 및 분석은 다중 과제에 직면해 있습니다. 본 논문은 유클리드 2-노름 하에서 대수 연산만을 사용하여 CDTW를 정확하게 계산할 수 없음을 증명하고, 근사 2-노름 하에서 CDTW를 계산하는 정확한 알고리즘을 제시합니다. 후자는 저자들이 확립한 기술적 기초에 의존하며, 이는 임의의 노름 및 관련 메트릭(예: 부분 Fréchet 유사성)으로 일반화될 수 있습니다.
본 논문은 2차원 공간에서 다각형 곡선 사이의 연속 동적 시간 정렬(CDTW) 거리를 효율적이고 정확하게 계산하는 방법을 연구합니다.
광범위한 실제 응용 : CDTW는 서명 검증, 지도 매칭, 궤적 클러스터링 등의 분야에서 중요한 응용을 가집니다.이산 방법의 한계 극복 : 전통적인 이산 DTW는 곡선의 연속성을 무시하며, 샘플링 속도가 다르거나 충분하지 않을 때 좋지 않은 결과를 생성합니다.견고성 요구사항 : Fréchet 거리의 최댓값 메트릭과 비교하여 이상치에 민감한 반면, CDTW는 경로 적분을 사용하여 이상치를 더 잘 처리할 수 있습니다.근사 및 휴리스틱 방법 : 많은 기존 방법은 입력 곡선을 이산화하여 적분을 처리하므로, 해의 품질 또는 실행 시간이 이산화 해상도에 의존합니다.차원 제한 : 이전의 정확한 알고리즘은 주로 1차원 경우에 국한되거나, 2D 유클리드 2-노름에서는 의사다항식 시간의 (1+ε)-근사 알고리즘만 존재합니다.이론적 이해 부족 : CDTW의 기본 성질, 특히 2D 및 서로 다른 노름 하에서의 성질은 아직 충분히 이해되지 않았습니다.저자들은 다음을 목표로 합니다:
CDTW의 계산 복잡성, 특히 2-노름 하에서의 비계산 가능성을 깊이 있게 이해 임의의 노름에 적용 가능한 기술적 기초 확립 곡선 이산화를 피하는 정확/근사 알고리즘 설계 국소 최적성 이론 (Section 2): 경로 적분 기반 CDTW 정의가 알고리즘 관점에서 유리한 국소 최적 매칭을 허용함을 증명하고, 계곡(valley)의 존재성 및 계산 방법을 확립하며, 이는 임의의 노름에 적용됩니다.비계산 가능성 결과 (Section 3): 유클리드 2-노름 하에서 CDTW와 관련된 수치가 초월수(transcendental numbers)일 수 있으므로, 대수 연산(ACMQ 모델)만을 사용하여 정확하게 계산할 수 없음을 증명합니다.다각형 노름 알고리즘 (Section 4): 다각형 등값 집합을 가진 노름 하에서 CDTW를 계산하는 정확한 알고리즘을 제시하며, 이 알고리즘은:입력 곡선의 이산화가 필요하지 않음 2-노름 하에서 (1+ε)-근사를 얻는 데 사용 가능 k ∈ O(ε^(-1/2))의 정다각형 노름 사용 기술적 성질 : 최적 함수 연속성, 전파 순서 등을 포함한 여러 기술적 성질을 확립하여 복잡도 분석의 기초를 마련합니다.입력 :
두 개의 2차원 다각형 곡선 P = ⟨p₀, ..., pₙ⟩ 및 Q = ⟨q₀, ..., qₘ⟩ 노름 ‖·‖ 출력 :
CDTW 정의 (Definition 1):
cdtw ∥ ⋅ ∥ ( P , Q ) : = inf ( f , g ) ∈ Π ( P ) × Π ( Q ) ∫ 0 1 d ( f ( t ) , g ( t ) ) ⋅ ∥ ( ∥ f ′ ( t ) ∥ ∥ g ′ ( t ) ∥ ) ∥ 1 d t \text{cdtw}_{\|\cdot\|}(P,Q) := \inf_{(f,g)\in\Pi(P)\times\Pi(Q)} \int_0^1 d(f(t), g(t)) \cdot \left\|\begin{pmatrix}\|f'(t)\| \\ \|g'(t)\|\end{pmatrix}\right\|_1 dt cdtw ∥ ⋅ ∥ ( P , Q ) := inf ( f , g ) ∈ Π ( P ) × Π ( Q ) ∫ 0 1 d ( f ( t ) , g ( t )) ⋅ ( ∥ f ′ ( t ) ∥ ∥ g ′ ( t ) ∥ ) 1 d t
여기서:
Π(P)는 P의 모든 단조이고 구간별 연속 미분 가능한 매개변수화를 포함 d(·,·)는 거리 함수(본 논문에서는 d(p,q) = ‖p-q‖ 사용) 1-노름을 사용하여 속도를 정규화하므로 경로 호장은 σ = ‖P‖ + ‖Q‖ 매개변수 공간 (Definition 2): 0, ‖P‖ × 0, ‖Q‖ , n×m 셀로 분할
계곡(Valley)의 개념 (Definition 4):
계곡 ℓ은 매개변수 공간의 직선(기울기 ≠ -1) 각 점 z ∈ ℓ은 싱크(sink): 기울기 -1인 직선을 따라 함수가 z에서 최솟값에 도달 핵심 정리 (Theorem 8):
임의의 노름 ‖·‖ 및 다각형 세그먼트 P, Q에 대해, 양의 기울기를 가진 계곡이 존재합니다. 이는 쌍대성 및 기하학적 분석을 통해 증명됩니다:
Lemma 7을 사용하여 직선 위의 게이지 함수 최소화 최대화 점 v*의 양의 성분 존재 증명 다각형 노름의 경우, O(1) 시간에 계곡 계산 가능 (Corollary 9) 최적 경로 특성 (Theorem 5):
계곡 ℓ이 주어졌을 때, 최적 (x,y)-경로는 다음과 같이 추적됩니다:
ℓ이 직사각형 Ξ = x₁,y₁ ×x₂,y₂ 과 교차하면, 경로는 x에서 x̂(ℓ 위)로, ŷ(ℓ 위)로, y로 이동 그렇지 않으면, 경로는 x에서 ξ(Ξ에서 ℓ에 가장 가까운 점)로 y로 이동 주요 정리 (Theorem 11):
정수 꼭짓점을 가진 곡선 P, Q를 구성하여:
(a) cdtw‖·‖₂(P,Q)는 초월수 (b) 각 최적 경로의 전환점 좌표는 모두 초월수 증명 개요 :
특정 2-세그먼트 곡선 P 및 3-세그먼트 곡선 Q 구성 적분 계산을 통해 로그 항을 포함하는 CDTW 값 도출 Baker 정리(초월수론의 결과, Lemma 10)를 이용하여 관련 수가 대수수가 아님을 증명 (b)의 경우, 도함수를 최소화하는 점도 초월수임을 증명 의의 : 이는 2-노름 하에서 CDTW가 근호로 표현될 수 없을 뿐만 아니라 대수수도 아니므로, 대수 연산을 사용하는 정확한 알고리즘이 존재할 수 없음을 나타냅니다.
알고리즘 프레임워크 : 매개변수 공간 셀을 통해 최적 경로 비용을 전파하는 동적 프로그래밍
노름 설정 :
1-부분수준 집합이 ψ(Rₖ)인 노름 Gψ(Rₖ) 사용 Rₖ는 정k각형(k는 짝수), 꼭짓점은 vᵣ = (cos(θᵣ), sin(θᵣ))ᵀ, θᵣ = r·2π/k ψ: ℝ² → ℝ²는 전단사 선형 매핑 핵심 성질 (Lemma 14):
(a) Gψ(Rₖ)는 O(1) 시간에 평가 가능하며, 각 원뿔에서 선형 (b) 전파 공간 ΣA,B는 O(k)개 직선의 배열로 표현 가능하며, optA,B는 각 면에서 이차 (c) 최적 함수 opt₀,B는 구간별 이차 전파 프로세스 (Algorithm Propagate):
입력: 곡선 세그먼트 P,Q, 경계 B, 인접 및 대변의 최적 함수
출력: B의 최적 함수(구간별 이차)
1. 계곡 ℓ 계산(O(1) 시간, Corollary 9)
2. 각 경계 A ∈ {adj(B), opp(B)}에 대해:
- O(k)개 직선의 배열 구성
- opt₀,A의 단절점에서 수직선과 오버레이
- 적절한 순서로 구간 처리(Lemma 15)
- 각 구간 I에 대해:
* 모서리 및 면의 후보 함수 수집
* 하부 포락선 계산(O(Xᵢlog(Xᵢ)α(Xᵢ)) 시간)
* 스택을 사용하여 결과 업데이트
3. 구간별 이차 최적 함수 반환
전파 순서 (Lemma 15):
해당 경로의 접미사가 교차하지 않음을 증명하여 올바른 전파 순서 결정:
A와 B가 같은 방향이면(예: A = opp(B)), s < s' A와 B가 직교하면(예: A = adj(B)), s > s' 근사 보장 (Corollary 17):
정k각형 노름 Gψ(Rₖ)를 사용하여 CDTW를 정확하게 계산 가능 k ∈ O(ε^(-1/2))을 통해 2-노름 하에서 (1+ε)-근사 획득 증명: 1/cos(π/k)² ≤ 1+ε 기하학적 쌍대 방법 : 게이지 함수의 쌍대 성질 및 볼록 기하학을 사용하여 계곡의 존재성 및 양의 기울기 성질 증명초월수론 응용 : CDTW에 Baker 정리를 처음 적용하여 이전의 "근호로 표현 불가"보다 더 강한 결과 증명이산화 회피 : 다각형 노름의 구간별 선형성을 활용하여 연속 곡선에서 직접 작업하며 이산화 불필요스택 기반 동적 프로그래밍 : 전파 순서 성질(Lemma 15)을 활용하여 스택 데이터 구조로 하부 포락선 계산 가속화통일된 프레임워크 : 확립된 기술 기초는 임의의 노름에 적용 가능하며 관련 메트릭(예: 부분 Fréchet 유사성)으로 일반화 가능주의 : 본 논문은 이론 논문이며, 주요 기여는 알고리즘 및 복잡도 분석이며 전통적 의미의 실험 부분을 포함하지 않습니다. 논문의 초점은:
초월성 검증 (Section 3):Theorem 11을 검증하기 위해 구체적 예시 구성 예시 (a): P = ⟨(1,2)ᵀ, (1,-4)ᵀ⟩, Q = ⟨(0,0)ᵀ, (5,0)ᵀ⟩ 계산 결과: cdtw = ½·ln(α₁) - 1/√2 - ½·ln(α₂) + √5 + 17√2 여기서 α₁ = √2-1, α₂ = √5-2 Lemma 10을 통해 이것이 초월수임을 증명 알고리즘 정확성 :Theorem 16이 Propagate 알고리즘의 정확성 증명 Theorem 20이 최적 함수의 연속성 증명 Lemma 19가 경로 수열 수렴성 증명 실행 시간 (Proposition 18):
총 시간: O(N·k²log(k)α(k)) N: 모든 최적 함수의 이차 세그먼트 총 개수 α: 역 Ackermann 함수 미해결 문제 :
1D 경우에는 N ∈ O(n⁵)임이 증명됨 2D 경우의 N의 다항식 경계는 아직 미확립 주요 어려움: 배열의 음의 기울기 직선이 이중화 행동을 초래할 수 있음(Figure 5) 비계산 가능성 :2-노름 하에서 CDTW가 초월수를 포함함을 명확히 증명 순수 대수 알고리즘의 가능성 배제 근사 알고리즘에 이론적 지원 제공 알고리즘 유효성 :다각형 노름 하에서 정확하게 계산 가능 2-노름의 (1+ε)-근사 획득, k = O(ε^(-1/2)) 입력 곡선의 이산화 불필요 통용성 :계곡 존재성은 모든 노름에 적용 가능 기술 프레임워크는 다른 메트릭으로 일반화 가능 예시 1 (Figure 4, Theorem 11a):
간단한 2-세그먼트 및 1-세그먼트 곡선 단일 매개변수 공간 셀 최적 경로는 3개 세그먼트: 계곡 위의 두 세그먼트, 수평 세그먼트 하나 CDTW 값은 로그 항을 포함하며, 초월수임이 증명됨 예시 2 (Figure 4, Theorem 11b):
3-세그먼트 및 1-세그먼트 곡선 두 개의 매개변수 공간 셀 최적 경로 후보 γₛ는 s ∈ 0,10 으로 매개변수화 C'(s) = 0의 해를 분석하여 최소화 점 s*가 초월수임을 증명 수치 검증: s* ≈ 2.08, 유일한 대수 후보 s₀* ≈ 4.31 표준 DTW : O(n²) 시간 Vintsyuk 1968 Fréchet 거리 : O(n²log n) 시간 Alt & Godau 1995 개선된 경계 : 더 나은 상한 존재 Gold & Sharir 2018; Buchin et al. 2017; Cheng & Huang 2025 Serra & Berthod 1994 : 첫 번째 연속 버전, 연속 매칭이지만 유한 합Munich & Perona 1999 : 동등한 정의 및 분석Serra & Berthod 1995 : 거리 벡터 변화 기반 평행이동 불변 버전Efrat et al. 2007 : 더 일반적인 적분 버전Buchin 2007 : 본 논문에서 채택한 정의Klaren 2020 , Buchin et al. 2022 : 1D 다항식 시간 정확 알고리즘Maheshwari et al. 2018 : 2D 유클리드 2-노름 하에서 의사다항식 시간 (1+ε)-근사Brankovic 2022 : 2D 1-노름 하에서 알고리즘De Carufel et al. 2014 : 부분 Fréchet 유사성은 근호로 계산 불가Bajaj 1988 , De Carufel et al. 2014 : 관련 기하학적 최적화 문제의 대수 차수본 논문 : 더 강한 초월성 결과사전식 Fréchet 거리 Rote 2014 부분 Fréchet 유사성 Buchin et al. 2009 이러한 메트릭은 CDTW와 국소 최적성 성질 공유 이론적 기여 :2-노름 하에서 CDTW의 정확한 계산은 초월수를 필요로 하며, 대수 계산 모델을 초과 임의의 노름 하에서 양의 기울기 계곡이 존재하여 알고리즘 설계의 가능성 보장 임의의 노름에 적용 가능한 기술 기초 확립 알고리즘 기여 :다각형 노름 하에서의 정확 알고리즘 제공 k = O(ε^(-1/2))의 정k각형을 통해 2-노름의 (1+ε)-근사 획득 입력 곡선의 이산화 회피 개방 문제 :2D 경우의 다항식 시간 경계 미확립 주요 과제는 배열의 음의 기울기 직선이 초래할 수 있는 이중화 행동 복잡도 분석 불완전 :O(N·k²log(k)α(k)) 경계를 제시했지만 N의 다항식 경계 미확립 1D의 O(n⁵) 분석 기법이 2D로 직접 일반화되지 않음 실제 효율성 미검증 :논문은 순수 이론이며 구현 및 실험 평가 없음 실제 실행 시간은 k와 N의 영향을 크게 받을 수 있음 근사 인자 의존성 :k = O(ε^(-1/2))는 작은 ε이 큰 k를 필요로 함을 의미 예를 들어 ε = 0.01은 k ≈ 314 필요 수치 안정성 :초월수의 정확한 계산을 회피했지만, 누적 오차 문제는 여전히 존재(Section 3 언급) 복잡도 분석 :2D 경우의 N의 다항식 경계 문제 해결 특히 Figure 5의 이중화 행동 처리 실제 구현 :알고리즘 구현 및 실험 평가 수행 기존 이산화 방법과 비교 추가 응용 :기술을 부분 Fréchet 유사성 등 관련 메트릭으로 일반화 다른 응용 분야 탐색 근사 개선 :동일한 근사 보장으로 더 작은 k 사용 가능 여부 연구 다른 근사 전략 탐색 이론적 깊이 :초월성 결과는 해당 분야의 중요한 이론적 기여이며, 이전의 "근호로 표현 불가"보다 더 강함 Baker 정리 사용의 증명 기법은 새롭고 엄격 계곡 존재성의 기하학적 증명은 우아하고 통용적 기술적 엄밀성 :모든 정리는 완전한 증명을 가짐(본문 또는 부록) 수학적 유도는 세밀하며, 특히 Appendix B의 초월성 상세 증명 다양한 경계 경우 고려 통용성 :확립된 프레임워크는 임의의 노름에 적용 가능 관련 메트릭(사전식 Fréchet 거리, 부분 Fréchet 유사성)으로 일반화 가능 Theorem 8 및 Lemma 15 등의 결과는 독립적 가치 보유 알고리즘 설계 :이산화 회피는 중요한 방법론적 기여 스택 기반 전파는 문제의 기하학적 구조 활용 Propagate 알고리즘은 명확하고 이해하기 쉬움 작성 품질 :구조가 명확하며 동기에서 이론을 거쳐 알고리즘으로 단계적 진행 그림(Figures 1-9)은 이해에 도움 관련 연구 종합이 포괄적 실험 부재 :알고리즘 논문으로서 실제 구현 및 실험 평가 부재 알고리즘의 실제 효율성 및 실용성 평가 불가 기존 방법과의 실제 성능 비교 부재 복잡도 미완성 :N의 다항식 경계는 핵심 개방 문제 이중화 행동 해결을 위한 명확한 방향 미제시 이는 알고리즘의 이론적 완전성 제한 근사 매개변수 :k = O(ε^(-1/2))의 의존 관계는 상대적으로 약함 작은 ε은 큰 k 필요로 하여 실제 효율성에 영향 가능 k의 실제 값이 성능에 미치는 영향 미논의 수치 문제 :초월수 계산 회피했지만 Section 3 언급의 누적 오차 문제 미충분 논의 구간별 이차 함수의 수치 안정성 미분석 응용 논의 :실제 응용 시나리오에 대한 논의 부족 CDTW를 DTW 또는 Fréchet 거리 대신 사용해야 할 시기 미논의 이론적 영향 :초월성 결과는 곡선 유사성 메트릭 계산 복잡성의 중요한 진전 근사 알고리즘의 필요성에 견고한 이론적 기초 제공 기술 프레임워크는 관련 문제 연구에 영감 가능 실용적 가치 :(1+ε)-근사 알고리즘은 실제 응용에 가치 있음 이산화 회피는 해의 품질 향상 가능 그러나 실제 효율성은 실험 검증 필요 재현성 :알고리즘 설명이 상세하여 이론적으로 재현 가능 그러나 구현 세부사항 및 코드 부재 부록의 상세 증명은 이해에 도움 후속 연구 :개방된 복잡도 문제는 후속 연구에 명확한 방향 제공 기술 기초는 다른 메트릭 및 응용으로 일반화 가능 새로운 알고리즘 설계 사고 자극 가능 이론 연구 :곡선 유사성 메트릭의 계산 복잡성 연구 기하학적 최적화 문제의 알고리즘 설계 초월수의 계산 기하학 응용 실제 응용 (잠재적):궤적 유사성 분석(샘플링 속도 차이가 클 때) 서명 검증(이상치에 견고성 필요) 지도 매칭(GPS 궤적 매칭) 시계열 클러스터링 부적합 시나리오 :실시간 계산이 필요한 응용(복잡도 높음) 고차원 데이터(현재 2D만 제한) 정확도 요구사항이 낮은 응용(DTW로 충분할 수 있음) Alt & Godau 1995 : Fréchet 거리의 고전 알고리즘Vintsyuk 1968 : DTW의 원래 정의Baker 1990 : 초월수론 기초(Lemma 10의 출처)Buchin 2007 : CDTW 정의의 출처Buchin, Nusser & Wong 2022 : 1D CDTW의 정확 알고리즘Maheshwari, Sack & Scheffer 2018 : 2D CDTW의 근사 알고리즘Brankovic 2022 : 2D 1-노름 하에서의 알고리즘Boyd & Vandenberghe 2009 : 볼록 최적화(게이지 함수)Schaefer & Wolff 1999 : 위상 벡터 공간(노름 이론)De Carufel et al. 2014 : 부분 Fréchet 유사성의 비계산 가능성종합 평가 : 이것은 CDTW의 계산 복잡성 및 알고리즘 설계 측면에서 중요한 기여를 한 고품질의 이론 계산 기하학 논문입니다. 초월성 결과는 해당 분야의 획기적인 진전이며, 기술 프레임워크는 우수한 통용성을 가집니다. 주요 부족점은 실험 평가 및 완전한 복잡도 분석의 부재입니다. 논문은 계산 기하학 최상위 학회(예: WALCOM, SoCG)에 발표하기에 적합하며, 이론 연구자 및 곡선 유사성 메트릭에 관심 있는 연구자에게 중요한 참고 가치를 가집니다.