KR20050084991A - 등각 구조에 의한 기하학적 표면의 분석 방법 - Google Patents

등각 구조에 의한 기하학적 표면의 분석 방법 Download PDF

Info

Publication number
KR20050084991A
KR20050084991A KR1020057008156A KR20057008156A KR20050084991A KR 20050084991 A KR20050084991 A KR 20050084991A KR 1020057008156 A KR1020057008156 A KR 1020057008156A KR 20057008156 A KR20057008156 A KR 20057008156A KR 20050084991 A KR20050084991 A KR 20050084991A
Authority
KR
South Korea
Prior art keywords
conformal
mapping
harmonic
mesh
geometric
Prior art date
Application number
KR1020057008156A
Other languages
English (en)
Other versions
KR100819960B1 (ko
Inventor
싱-텅 야우
시안펭 구
얄린 왕
Original Assignee
지오메트릭 인포매틱스 인코퍼레이티드
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 지오메트릭 인포매틱스 인코퍼레이티드 filed Critical 지오메트릭 인포매틱스 인코퍼레이티드
Publication of KR20050084991A publication Critical patent/KR20050084991A/ko
Application granted granted Critical
Publication of KR100819960B1 publication Critical patent/KR100819960B1/ko

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00Scenes; Scene-specific elements
    • G06V20/60Type of objects
    • G06V20/64Three-dimensional objects
    • G06V20/653Three-dimensional objects by matching three-dimensional models, e.g. conformal mapping of Riemann surfaces
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T15/003D [Three Dimensional] image rendering
    • G06T15/10Geometric effects
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/40Extraction of image or video features
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00Scenes; Scene-specific elements
    • G06V20/60Type of objects
    • G06V20/64Three-dimensional objects

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Multimedia (AREA)
  • Software Systems (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Evolutionary Computation (AREA)
  • Data Mining & Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Artificial Intelligence (AREA)
  • Geometry (AREA)
  • Computer Graphics (AREA)
  • Processing Or Creating Images (AREA)
  • Image Generation (AREA)

Abstract

기하학적 표면을 분석하고, 분류하며, 인식하는 방법이 개시된다. 기하학적 표면은 리만 다면체로 취급되며, 표면에 상응하는 등각 구조가 계산된다. 표면의 등각 구조는 표면에 대한 내부의 기하학 정보를 포함하지만, 다른 표면들과 비교할 때 훨씬 더 많이 치밀하다. 표면을 원반, 구, 또는 평면과 같은 기준 파라미터 도메인에 등각으로 매핑시키는 것은 표면의 기하학적 정보를 보유하며, 등각 구조의 계산을 훨씬 더 쉽게 한다. 그러한 등각 표현에 의해 가능해진 다양한 애플리케이션은 표면 매칭, 표면 카탈로깅, 표면 인식, 표면들간의 애니메이션 및 모핑, 및 다른 수학적 분석을 포함한다.

Description

등각 구조에 의한 기하학적 표면의 분석 방법{ANALYSIS OF GEOMETRIC SURFACES BY CONFORMAL STRUCTURE}
본 출원은 표면의 분석에 관한 것이며 특히, 컴팩트한 리만 표면이론을 컴퓨터로 계산가능한 알고리즘으로 변환하여 표면 분석을 위한 기초적인 기하 도구를 제공함으로써 표면의 등각 구조를 계산하여 표면을 분석하는 것에 관한 것이다.
기하 표면의 분류 및 식별은 컴퓨터 그래픽 및 컴퓨터 캐드(CAD) 분야에서 기본적인 문제이다. 스캐닝 및 이미징 기술이 발전하면서, 다수의 채색된 메쉬가 데이터베이스, WWW, 그리고 인터넷에서 이용 가능해지고 있다. 또한, MRI나 PET 이미징과 같은 의료 이미징 기술이 신체 내부구조의 3차원 모델을 생성할 수 있다. 예를 들면, 뇌 이미징에서의 최근의 발전은 뇌 지도의 데이터베이스에서 이와 같은 이미지의 수집과 저장을 가속시켜 왔다. 마찬가지로, 바이오메트릭 보안 응용에서, 안면 인식은 이미징, 저장, 및 미리 저장된 안면들에 대한 3차원 안면 특징들의 매칭을 포함한다. 또한, 3차원 웹페이지를 사용하는 엔터테인먼트 시스템은 그 수가 증가하고 있으며, 모핑(morphing)과 텍스처(texture) 매핑과 같은 컴퓨터 애니메이션 기술 또한 3차원 표면의 생성과 조작을 포함한다.
이 모든 응용들에 있어서, 기하적인 데이터는 미분구조(differential structure) 대신에 조합구조(combinatorial structure)를 갖는 삼각 메쉬로서 표현된다. 따라서, 이들 구조들을 미분 기하 기법을 사용하여 처리하는 것은 어렵다. 현재의 분석 방법은 2개 표면 사이의 Hausdorff 거리를 측정하나; 상기 표면들 사이에 일치를 발견하는 일반적인 접근방법은 없고, 또한 조합 검색은 비효율적이다. 또한, 현재의 표면 분석 방법은 표면의 삼각법(triangulation)과 분해능(resolution)에 크게 의존하고 있다. 그러나, 다른 삼각법과 분해능은 크게 변화하는 결과를 가져올 수 있다. 결국, 기하 표면 데이터는 극히 크다. 하나의 표면이 수백만의 꼭지점과 면들을 가질 수 있으므로 현재의 시스템에서 필요한 순수한 계산 횟수는 효과적이고 효율적인 알고리즘을 만드는 것을 극히 어렵게 한다. 또한, 현재로서는 토폴로지 불변량을 사용하여 표면을 분류하는 효과적이고 일반적인 방법이 없는데 이는 상기 분류가 보통 너무 조악하기 때문이며, 또한 상기 분류가 너무 엄격하기 때문에 유클리드 기학 불변량을 사용하는 방법도 없다.
따라서, 표면의 형상에만 의존하는 점에서 본질적인 시스템이고; 표면들을 효율적으로 분류하여, 동일 분류에서 2개의 표면 사이에 일치를 발견하는 일반적인 방법을 제공하고; 수적으로 안정하고 정확한 효율적이면서 달성가능한 계산(computation)을 제공하는 기하적 분석 방법을 제공하는 것이 유용하다.
본 발명은 첨부된 도면과 함께 다음의 상세한 설명으로부터 더욱 완전하게 이해될 수 있다:
도 1A는 사람의 안면과 구체 사이의 등각 매핑을 도시하고,
도 1B는 도 1A의 사람 안면으로부터 도 1B의 평면에 매핑된 체커 보드 텍스처를 도시하고,
도 2A 내지 2D는 2홀 토루스(torus)의 홀로모픽 1-형의 다양한 성분들을 도시하고,
도 3은 구체에서 가고일 모델(gargoyle model)에 대한 구형의 등각 임베딩을 도시하고,
도 4는 구체에 등각으로 매핑된 인간의 뇌 모델을 도시하고,
도 5는 상기 단위 구체에 매핑된 토끼 모델을 도시하고,
도 6A 및 6B는 파라미터화의 영점을 도시하고,
도 7A 내지 7D는 지너스 2 및 3인 토루스에 대한 글로벌한 등각 환추(atlas)를 도시하고,
도 8A 내지 8D는 2개의 지너스-1 토루스의 등각 등가물(conformal equivalence)이 아닌 토폴로지 등가물(topological equivalence)을 도시하고,
도 9A 내지 9D는 상이한 등각 구조를 갖는 지너스-1 표면을 도시하고,
도 10A 내지 10D는 글로벌한 등각 파라미터화의 균일성에서의 향상을 도시하고,
도 11A 내지 11D는 상이한 등각 구조를 갖는 다양한 지너스-2 표면을 도시하고,
도 12A 및 12B는 토끼 표면의 삼각법의 규칙화의 사용을 도시하고,
도 12C 및 12D는 등각 기하 이미지로부터 토끼 표면의 재구성을 도시하고,
도 13A는 뇌 표면 모델을 도시하고,
도 13B는 구체로 등각 매핑된 도 13A의 뇌 표면 모델을 도시하고,
도 13C는 도 13A의 뇌 표면 모델의 구형 기하 이미지를 도시하고,
도 13D는 도 13C가 256배 압축된 후 재구성된 뇌 표면을 도시하고,
도 14는 등각 구조를 사용하여 여성 안면으로부터 남성 안면으로의 기하 모핑을 도시하고,
도 15A 및 15B는 삼각법의 원시 레벨에서 차주전자 모델의 전체적인 파라미터화를 도시하고,
도 15C 및 15D는 삼각법의 단순화된 레벨에서 차주전자의 전체적인 파라미터화를 도시하고, 및
도 16A 내지 16D는 4개의 높은 지너스 표면에 대한 전체적인 파라미터화 결과를 도시하는 도면이다.
기하 표면을 분석, 분류 및 인식하는 방법이 개시된다. 기하 표면들은 리만 다양체(Riemann manifolds)로서 취급되고 상기 표면에 대응하는 등각 구조가 계산된다. 상기 표면의 등각 구조는 표면에 관한 본질적인 기하 정보를 포함하지만, 훨씬 더 콤팩트한 포맷이다. 그러나 일반적으로 표면들은 상당히 큰 메쉬 데이터 포인트 수를 갖는 복수의 메쉬 데이터로서 표현된다. 이와 같은 메쉬 표면의 등각 구조를 계산하는 것은 많은 메쉬 데이터 포인트의 갯수와 필요한 훨씬 더 많은 계산 횟수 때문에 어려운 일이다. 상기 표면을 디스크, 구체, 또는 평면과 같은 정규 파라미터 도메인으로 등각으로 매핑하는 것은 상기 표면의 기하 정보를 계속 유지하고, 등각 구조의 계산을 훨씬 용이하게 한다.
특히, 일 실시예에서, 제 1 및 제 2 표면은 제 1 및 제 2 매핑 표면을 형성하는 정규 파라미터 도메인으로 등각으로 매핑된다. 각각의 매핑 표면에 대한 등각 파라미터화는 상기 표면들이 매칭하는지를 판단하기 위해 계산되고 서로 비교된다.
다른 실시예에서, 표면을 분류하는 방법이 개시되며, 이 방법에서은 표면이 상기 등각 파라미터화에 따라 분류된다. 구체적으로는, 상기 표면에 대응하는 주기 매트릭스(R)가 결정되어 저장된다. 그 후, 특정 표면에 대한 검색이 미리 저장된 주기 매트릭스(R)를 조사하고 이 매트릭스를 원하는 표면에 대응하는 제 2 주기 매트릭스(R')에 비교함으로써 수행될 수 있다.
또 다른 실시예에서, 표면 인식에 대한 방법이 제공된다. 구체적으로는, 표면을 나타내는 메쉬가 제공되고 하나 이상의 피처 포인트가 순차적으로 제거된다. 제거된 각 피처 포인트에 대해 대응하는 주기 매트릭스(R)가 계산된다. 얻어진 주기 매트릭스의 순서를 미리 계산된 주기 매트릭스의 순서와 비교함으로써, 표면이 인식된다.
또는, 모든 피처 포인트가 한번에 제거될 수 있고 한 포인트가 상기 표면 내에서 선택된다. 이 포인트가 표면을 소정의 궤도로 이동하면서, 주기 매트릭스의 순서가 계산되고 주기 매트릭스의 미리 계산된 순서와 비교된다.
또 다른 실시예에서, 이미지 압축 방법이 개시된다. 표면을 나타내는 메쉬가 제공되고 상기 메쉬에 대한 등각 파라미터화가 계산된다. 상기 등각 파라미터화를 사용하여, 평균 곡률이 계산될 수 있으며, 이들 2개의 파라미터를 가지고 원래 표면이 유일하게 결정된다.
다른 실시예에서, 의료 이미징에 대한 응용이 개시된다. 예를 들어 뇌나 다른 기관의 의료 이미지는 전형적인 지너스-0 표면(genus-zero surface)이다. 상기 지너스-제로 표면을 구체로 등각 매핑하는 것은 상기 표면을 분석하는 것을 가능하게 한다.
다른 실시예에서, 표면에 애니메이션하는 방법이 개시된다. 2개의 유사한 구체가 주어진 경우, 각 표면으로부터 피처 포인트들이 제거되고 각 표면의 배증(doubling)이 계산된다. 각 표면은 하나 이상의 패치로 분해되고 각 패치는 평면으로 매핑된다. 한 평면으로부터 다른 평면으로의 등각 매핑이 결정되고 제어 포인트를 선택한 후 BSpline곡선이나 다른 부드러운 곡선 함수가 상기 2개의 평면간 부드러운 전이를 발생시키기 위해 사용된다.
다른 실시예에서, 주어진 표면을 커버하기 위한 텍스처를 발생시키는 방법이 제공된다. 상기 표면은 등각 파라미터화를 사용하여 평면 표면과 같은 정규 파라미터 표면으로 매핑되며, 상기 텍스처는 상기 파라미터 표면을 위해 계산된다. 상기 텍스처 패치들을 전체적으로 부드럽게 하기 위해, 텍스처 패치간의 경계를 확산시키는데 디히클레(Dirichlete) 방법이 사용된다. 이런 방식으로, 텍스처 패치는 성장하고 함께 꿰매진 다음 상기 파라미터 표면으로 매핑된다.
상술한 방법 및 시스템의 다른 형태, 피처, 및 측면들이 다음에서 보다 상세하게 설명된다.
다음 실시예에서 2차원 표면들은 리만 표면으로 취급되며 상기 표면에 대응하는 등각 구조가 계산된다. 방위를 맞출 수 있는(orientable) 모든 면들은 리만 표면이고, 등각 변형하에서 불변인 고유한 등각 구조를 갖는다. 일반적으로, 상기 등각 구조는 토폴로지 구조보다 더욱 정교하고 미터법 구조보다 덜 정확하다. 지너스-1 표면에 대해서, 모든 등각 구조의 공간은 2차원이다. 그리하여, 2개의 파라미터를 사용함으로써, 모든 지너스-1 표면들이 분류될 수 있다. 일반적으로, 지너스 g 표면에 대해서, 모든 가능한 등각 구조의 공간은 6g-6 차원이다. 그리하여, 모든 지너스 g 표면들은 g x g 복소 매트릭스(complex matrix)를 사용하여 분류될 수 있다.
2개의 표면 사이에 등각 등가물을 시스템적으로 계산하는 방법이 제공된다. 구체적으로, 동일한 등가 구조를 갖는 임의의 2개 표면에 대해, 상기 2개 표면 사이에 등각 1 대 1 매핑을 시스템적으로 계산하는 방법이 제공된다. 지너스-0 표면에 대해, 그와 같은 매핑의 그룹은 6차원적이다. 지너스-1 표면에 대해, 그와 같은 그룹은 2차원적이다. 지너스-1 이상의 표면에 대해, 그와 같은 매핑의 그와 같은 그룹은 단지 1차원만을 포함한다. 그리하여, 유리한 것은, 이하에서 설명되는 방법들이 최적의 매핑을 발견하고 동일한 등각 구조를 갖는 임의의 2개 표면 사이의 Hausdorff 거리를 측정하는 효율적인 방법을 제공한다.
표면의 등가 구조가 결정되는 이하에서 설명된 방법들에 대해, 상기 등가 구조는 상기 표면의 기하의 함수일 뿐이다. 그것은 삼각법과 분해능 어느것에 의해서도 영향을 받지 않으며 또한 등각 매핑은 상기 표면의 형상을 보존한다.
모든 면들이 리만 표면이라는 것은 잘 알려져 있다. 임의의 리만 표면은 등각 좌표 환추(conformal coordinate atlas), 또는 등각 구조를 갖는다. 등각 변형은 하나의 등각 구조를 하나의 등각 구조에 매핑한다. 2개의 리만 표면 사이의 등각 변형에 의해 어디에서든지 각들이 유지된다. 알려진 바와 같이, 1차원 접속된 복소 다양체(complex manifold)는 리만 표면으로 알려져 있다. 리만 균일화 정리에 의해, 모든 표면들은 정규의 공간(canonical space)에 전체적으로 등각으로 매립될 수 있다. 상기 정규 공간은 보통 디스크, 평면, 또는 구체이며, 그 선택은 표면의 고유의 기하구조에 의해 결정된다. 상기 등각 매립된 표면은 상기 정규 공간 상에 매립된 원시 기하 정보의 많은 부분을 포함한다. 등각 매립을 통해서, 3차원 매칭 문제가 이들 3개의 정규 공간에서 2차원 매칭 문제로 변환될 수 있다. 이하에서 보다 상세히 설명된 바와 같이, 이 방법은 비-강체의 기형의 표면 매칭을 위한 잠재성을 갖는다.
표면을 상기 정규 공간에 매립하는 방법은 상기 표면의 등각 구조를 반영한다. 구체적으로는, 하나의 표면으로부터 정규 공간으로의 모든 전체적인 등각 매립은 특수 그룹을 형성한다. 만일 2개의 표면이 서로에게 등가으로 매핑될 수 있다면, 그것들은 동일한 그룹 구조를 공유한다. 다시 말해, 그와 같은 그룹 구조는 완전 등각 불변량이다. 그러므로, 등각 불변량을 사용하여 모든 구조를 분류할 수 있다. 토폴로지상의 등가 분류 각각에 대해, 무수한 등각 등가 분류가 존재한다. 이것은 표면 분류화 문제를 위해 가치있는 것이다.
S1과 S2를 2개의 정규 표면으로 하고, (x1, x2)에 의해 파라미터화 되는 것으로 하자. 로컬 좌표에서 로서 표현되는 매핑 φ: S1 -> S2를 정의한다.
S1과 S2의 제 1 기초 형태(리만 메트릭스(Riemann metrics))를:
로 하자.
에 의해 유도된 S1에 대한 풀 백 메트릭(pull back metric)은
이다.
만약 다음과 같은 함수 λ(x1, x2)가 존재한다면,
,
φ는 S1과 S2 사이에서 등각 매핑이라고 할 수 있다. 구체적으로, S1으로부터 로컬 좌표 평면(x1, x2)으로의 매핑이 등각이라면, (x1, x2)는 S1의 등각 좌표이고, 또한 등온 좌표로서 지칭된다. 도 1A는 사람의 안면과 평면 상의 정사각형 사이의 등각 매핑을 도시한다. 도 1B는 체커보드를 표면에 텍스처 매핑함으로써 상기 매핑의 등각 성질을 도시한다. 도 1A 와 1B를 정밀 검사는 상기 체커보드 상의 모든 직각은 상기 표면의 텍스처 상에 유지된다는 것을 설명해준다. 도 16은 높은 지너스를 갖는 4개의 표면들 즉 지너스 > 1 인 표면의 글로벌 파라미터화 결과를 도시한다. 도시된 바와 같이 상기 체커보드 패턴 상의 모든 각들은 상기 매핑의 등각 성질을 나타내는 직각이다.
복소 다양체에 대해, U ⊂ C 는 개방 셋이라고 가정하고 f는 복소함수 f: U →C 라고 하자. 그러면 임의의 z0 ∈ U 에 대해 ε> 0 라면, f는 홀로모픽(holomorphic)하다고 한다. 그러므로 디스크에서
이면,
f는 다음과 같은 수렴 파워 시리즈(power series)로서 표현된다.
U ⊂ C와 V ⊂ C를 C의 개방 셋이라고 하자. f가 1 대 1 이고 홀로모픽하다면 매핑 f: U →V는 바이홀로모픽(biholomorphic) 하며, f-1: V →U 역시 홀로모픽하다.
S를 다음과 같은 3개의 조건을 만족하는 {(Uj, zj)}jεJ 군(family)과 연결된 Hausdorff 공간이라 하자:
1. 모든 Uj는 S의 개방 서브셋이고, S = ∪jεJUj
2. 모든 zj는 상기 복소 평면에서 개방 서브셋 Dj로의 Uj의 유사형(homeomorphism)이다
3. 만일 Uj ∩Uk ≠φ라면, 전이 매핑
은 바이홀로모픽 매핑이며, 역시 홀로모픽 유사형이다.
이와 같이, {(Uj, zj)}jεJ 은 S 상의 좌표 이웃들(coordinate neighborhoods)의 시스템이고 S상의 1차원 복소 구조를 정의한다. 리만 표면의 좌표 이웃 (U, z)는 S에서 한 쌍의 개방 셋(U)이고 상기 복소 평면으로의 U의 유사형 z이다. U는 S의 좌표 이웃으로서 지칭되며, 유사형 z는 로컬 좌표 또는 로컬 파라미터로서 지칭된다.
일반적으로, 리만 표면 R로의 S의 매핑 f는, w·f·z-1이 f(U) ⊂ V에 의해 S의 (U, z)과 R의 (V, w)의 모든 좌표 이웃들에 대해 홀로모픽하다면, 홀로모픽 매핑이라고 한다. 바이홀로모픽 매핑 f:S →R은 R상으로의 S의 홀로모픽 매핑 f가 홀로미픽 역매핑 f-1: R →S를 갖는 것을 의미한다.
이와 같이, 2개의 리만 표면 S 와 R은 둘 사이에 바이홀로모픽 매핑이 존재한다면 바이홀로모픽 등가이다. 이와 같은 매핑이 존재한다면, S와 R은 같은 리만 표면으로 간주되고 S와 R은 같은 등각 구조를 갖는다. 일반적으로, 복소 구조, 바이홀로모픽 매핑 및 바이홀로모픽 등가는 각각 등각 구조, 등각 매핑 및 등각 등가라고도 한다.
표면 S가 와 같은 리만 메트릭을 갖는 것으로 하면, 상기 메트릭은 등각 구조 {(Ui, zi)}를 유일하게 결정하는데 사용될 수 있으며, 따라서 좌표이웃 (Ui, zi) 에서 ds2의 로컬 표현은
이다.
여기서 λ(zi)는 양의 실수 함수이다.
상기 등각 구조를 결정하기 위해 S에에서 모든 홀로모픽 미분 형태들이 발견되어야 한다. S를 리만 표면이라고 하면, S에서 홀로모픽 미분형태 ω는 다음과 같은 2개의 조건을 만족하는 {(Ui, zi, ωi)}의 군으로 주어진다:
1. {(Ui, zi)}를 등각 구조라고 가정하면, ω는 다음과 같은 로컬 표현 ω1을 갖는다.
여기서, fi는 Ui에 대한 홀로모픽(holomorphic) 함수이다.
2. 만약 에 대한 좌표 이행이라면,
그리하여, 로컬 표현은 연쇄 규칙을 만족한다.
S에 대한 모든 홀로모픽 미분의 집합은 로 나타내며, 는 S의 동상동(cohomology) 그룹에 대해 등정형이다. 그리하여, 을 계산하기 위해서는 S의 상동 그룹이 계산되어야만 한다.
S를 메트릭 g, 조밀한 2차원 다면체를 갖는 2차원 리만(Riemann) 다면체라고 하자. C1에 대해,
에 대한 로컬 좌표 에서의 에너지 밀도라고 하자. φ의 C1 이형(variation)은 파라미터 에 부드럽게 의존하는 C1의 패밀리 이며, 그리하여 이다. φ의 이형 은 만약 모든 에 대한 S/ω에 대해서 인 콤팩트 집합 이 존재한다면 조밀하게 지원된다.
C1 에 대한 하모닉(harmonic) 맵은 조밀하게 지원된 이형에 관한 디리치렛(Dirichlet)의 에너지에 대해 정적인 맵 이며, 다음으로 주어진다.
로컬 좌표에서 이며, 여기서 이다. 맵 은 다음을 만족할 때에만 하모닉하다.
여기서, λ는 S에 대해 전역적으로 정의된 함수이며, 는 N에 대한 이미지 점에서 수직이다. 유-영(genus-zero) 표면에 대해서, 하모닉 매핑은 등각 매핑이다. 만약 N이 R이라면, 은 하모닉 함수로 불린다. 모든 등각 맵은 하모닉하지만, 모든 하모닉 맵이 등각은 아니라는 것을 주의하라.
S에 대한 실수 미분 1-형태 τ는 만약 S에서의 임의의 점에 대해, 다음의 식을 만족하는 열린 집합 가 존재한다면, 하모닉하다.
여기서, f는 S에 대한 하모닉 함수이며, d는 외부-미분 연산자이다.
모든 하모닉 미분은 동상동 그룹 에 등정형인 특수 그룹 H를 형성한다. 호지(Hodge) 이론에 따르면, 각각의 동상동 클래스에는, 고유한 하모닉 미분 형태가 있다.
홀로모픽 1-형태 ω는 가 되는 두 개의 실수 미분 1-형태 로 분해될 수 있으며, 여기서 는 모두 하모닉이다. 하모닉 1-형태를 평면에 대해서 적분함으로써, 평면은 복소수 평면으로 등각으로 매핑될 수 있다.
모든 하모닉 1-형태는 상동 그룹 H1(S,Z)에 듀얼(dual)이다. 유(genus) g 평면 S에 대해서, H1(S,Z)의 2g 생성기가 있다. 각각의 핸들에 대응해서, 다음의 식이 되는 두 개의 생성기가 있다.
여기서, 은 두 개의 닫힌 곡선의 대수적 교차 수를 나타낸다. 은 기준 상동 기초라 불린다. 만약 가 H1(S,Z)의 기초라면, 듀얼 홀로모픽 1-형태 기초는 이며, 다음을 만족한다.
도 2A-2D는 4개의 닫힌 곡선을 구성하는 도 2A에서의 2-구멍 토러스(torus)의 상동 기초를 묘사한다. 도 2B는 빗금친 곡선이 ω의 적분 선들인 e1에 대한 하모닉 1-형태 ω듀얼을 묘사한다. 도 2C는 도 2B에 묘사된 하모닉 1-형태에 직교인 공액 하모닉 1-형태 ω*를 묘사한다.
등각 등가에 대한 완전 불변식은 복소수 행렬에 의해 제공된다. 가 기준 상동 기초이고 의 기초라고 가정하면, 행렬 P=(pij)은 S의 주기 행렬이라 불리며, 여기서
이다.
P1과 P2로 주어지는 2개 평면의 주기 행렬에 대한 검사는 2개 평면 사이의 등각 매핑을 계산할 필요없이, 2개 평면이 서로에 대해 등각으로 동등한지를 판정할 수 있다.
일반적으로, 평면은 삼각형 메시(meshes)에 의해 표현된다. 모든 심플리셜(simplicial) 표면은 자연스런, 아래에 놓이는 복소수 구조를 갖는다. K를 심플리셜 복소수라고 하고, 매핑 이 R3를 끼워넣는다면, M(K,f)는 삼각형 메시라 불리며, Kn(여기서 n=0,1,2)은 n-심플리시의 집합들이다. 은 n-복소수를 나타내며, 이고, 이다.
연쇄 공간은 심플리시들의 선형 결합이며, 다음의 식으로 주어진다.
Cn(n=0,1,2)에서의 요소들은 n-연쇄라고 불린다. 또한, 모든 페이스(faces) 의 합은 C2에 있으며, M은 또한 이러한 2-연쇄를 나타내기 위해 사용된다.
연쇄 공간들 중의 경계 연산자 는 선형 연산자이다. 라고 하면,
그다음 Cn 의 n-체인에 대해, 바운더리 오퍼레이터는
로서 정의된다.
1의 영공간(null space)을 나타내기 위해서, ker∂1⊂C1 는 M 상의 모든 닫혀진 루프를 나타낸다. 유사하게, img∂2⊂C1 는 모든 표면 패치 바운더리를 나타내는 ∂2의 이미지 공간을 나타낸다. ∂1ㆍ∂2 = 0 이므로, img∂2 ⊂ ker∂1이다. 그러므로, M의 호몰로지 그룹(homology group), Hn(M, Z)는 다음과 같이 주어진다.
H 1(M, Z)는 M 상의 임의의 표면 패치의 바운더리가 아닌 모든 닫혀진 루프를 나타낸다. M의 위상은 H1(M, Z)에 의해 판정된다.
M이 지너스(genus) g의 닫혀진 메쉬가 되게 하고, B = { 1, 2,..., 2g} 그 호몰로지 그룹의 임의의 기저가 되게 하자. 그러면 B의 교차 행렬 C는 다음과 같이 주어진다.
여기서 ㆍ은 교차의 개수를 나타내며, 교차점에서의 ei 및 ej 의 탄젠트 벡터의 크로스 프로덕트(cross product)의 방향이 정상적인 방향과 일치하는 경우 +1 을 하고, 그렇지 않으면 -1 을 한다.
코-체인 공간은 R로의 체인 공간 사이의 위상동형(homeomorphism)의 세트이고 다음과 같이 주어진다.
여기서 Hom(Cn, R)은 R로의 Cn 사이의 모든 위상동형의 세트를 나타낸다. Cn 의 요소는 n-코체인 또는 n-폼으로 불린다. 코바운더리 오퍼레이터는 δn : Cn →Cn+1로서 정의된다. ωn∈Cn 가 n-폼이 되게 하고, 가 n+1 체인이 되게 하자, 그러면
,
δ1ㆍδ2 = 0.
코호몰로지 그룹 H n (M,R)은 아래와 같이 정의된다.
kerδ1에서 1-폼은 닫혀진 1-폼으로 불려지고 imgδ0에서의 1-폼은 이그젝트(exact) 1-폼으로 불려진다. 두개의 닫혀진 1-폼은 그들이 이그젝트 1-폼과 다르다면 코호몰로지한 것으로 불려진다. 코호몰로지 그룹 H 1 (M,R)은 호몰로지 그룹 H 1 (M,Z)와 동형이다.
n-체인을 따라 n-폼의 집적은 cn ∈ Cn 및 ωn ∈ Cn 인 경우 아래와 같이 정의된다.
바운더리 및 코바운더리 오퍼레이터는 Stokes 공식
에 의해 관련된다.
웨지 프로덕트는(Wedge product)는 쌍일차(bilinear) 오퍼레이터 ∧: C1xC1 →C0 이다. f∈K2 가 M상의 페이스, δ2f=e0+e1+e2, ω,∈ C1 되게 하면,
이다.
쌍일차 오퍼레이터 스타 웨지 프로덕트 ∧*: C1xC1 →C2 는 유사하게 정의된다. f∈K2 이고, 세개의 에지의 길이를 l0, l1, l2, 및 f의 영역이 A라고 하면,
여기서,
,
,
및 2차 폼 G는 폼
을 갖는다.
닫혀진 1-폼의 조화 에너지 ω는
여기서,
및 e가 바운더리 에지, e∈δ2M 라면, e는 하나의 페이스 f0에 부착되고 그다음 we 는 다음과 같이 주어진다.
조화 에너지를 최소화시킨다면, 즉
로서 한정된 Laplacian 오퍼레이터가 0과 같다면, 닫혀진 1-폼은 조화 1-폼으로 불려진다. 따라서, Laplacian 오퍼레이터가 영인 경우에 한해 닫혀진 1-폼은 조화된다. M이 호몰로지 기저 {r1, r2,...,r2g} 및 조화 1-폼 기저 {ω1, ω2,...,ω2g}를 갖게 하고,
여기서, - i ㆍ j i j 대수학적인 교차 개수이면, 호몰로지 기저 및 조화 1-폼 기저는 서로 듀얼(dual)이 되는 것으로 알려져 있다.
M이 메쉬가 되게 하면 N은 R3 의 부드러운 표면이다. 구분적인 선형 맵 u:M →N ⊂ R3
으로서 M ~ N 의 모든 정점을 도시한다.
u = (u1, u2, u3)의 조화 에너지는 다음과 같이 주어진다.
여기서, E(δuα)는 1-폼 (δuα)에 대해 정의된 조화 에너지이다. 만약 u가 조화 에너지, E(u)를 최소화시키면, u는 조화 맵이고 다음의 조건을 만족한다.
여기서, n은 N 상의 노말(nomal) 필드이다.
조화 1-폼 ω가 제공되면, 고유의 켤레 조화 1-폼 ω* 이 존재한다. 홀로모로픽(holomorophic) 1-폼은
으로서 정의된다.
모든 홀로모로픽 1-폼은 H1(M,R)와 동형인 그룹 1(M)을 형성한다. 1(M)의 기저는 조화 1-폼 그룹의 기저로부터 직접 구성될 수 있다. {ω1, ω2,...,ω2g} 의 기저를 갖는 조화 1-폼 그룹이 제공되면, 1(M)의 기저는 {ω1+ ω1 *2+ ω2 *,...,ω2g+ ω2g *,} 에 의해 제공된다.
제공된 B = { 1, 2,..., 2g-1, 2g}가 H1(M,Z)의 기저이고, B* = {ω1, ω2,...,ω2g-12g}가 1(M)의 듀얼 기저이면, 행렬 C2gx2g=(cij) 및 행렬 S2gx2g=(sij)는
,
로서 정의된다.
그런 다음 M의 주기 행렬 R은
CR=SI
과 같이 정의된다.
여기서 R은 R2=-I를 만족한다. 상기 행렬(C, R)은 M의 등각동등한 클래스를 결정한다. 특히, (R1, C1)과 (R2, C2)을 갖는 임의의 두 표면 M1, M2 각각에 대해, M1과 M2 는 정수 행렬 N이 존재한다면 서로 등각 동등하여,
이 된다.
종수 g>인 메시의 등각 구조는 {(Ui, zi)}의 패밀리여서,
1. Ui는 단순히 연결되고, M의 면에 의해 형성된다.
2. M ⊂∩Ui
2. zi는 모두 선형이고, 홀로모픽 1 폼 가 존재하여 이다.
종수가 0인 메시에 대해, 홀로모픽(holomorphic) 1 폼은 존재하지 않는다. 이경우, 종수가 0인 표면은 단위구 S2의 표면으로 등각으로 매핑될 수 있고, S2의 등각 구조는 M의 등각 구조를 정의하는 데에 사용될 수 있다. 따라서, 불연속 조화맵 u:M→ S2는 M의 등각구조를 정의한다. 임의의 표면에 대해, c ∈ C1를 따라 M을 자름으로써, D위상적인 디스크(DM)이 형성될 수 있고, 그와 함께 특수한 1-체인이 형성될 수 있다. 이 c를 따라 자는 것은 궤적 또는 커트 그래프라고 하고, DM은 M의 기본 영역이 된다. c의 선택은 고유하지 아니하며, 따라서 기본 영역이 없다.
등각 맵 u:DM→ C는 홀로모픽 1 폼 을 사용하여 발견될 수 있다. 이 선택되고 임의의 꼭지점(vertex) 에 대해 임의의 경로 이 선택되어, 이고,
이다.
전술한 바와 같이, 모든 종수 0인 표면은 등각 동등하다. S2에서 그자신으로의 모든 등각 맵은 6차원 뫼비우스 변환 그룹을 형성한다. 복소 평면으로 상기 구면을 매핑하기 위한 스테레오 그래픽 사영을 사용하여, 모든 뫼비우스 변환은
의 폼이 된다.
그러나, 뫼비우스 변환의 폼에 따라 고유 해답을 형성하기 위해, 구면으로 종수 0인 표면을 매핑하는 등각 맵을 연산하기 위한, 뫼비우스 변환에 대한 추가적인 제약들이 필요하다.
다른 어려움은 맵의 이미지가 S2 상에 있고 R3 상에는 없는 것이다. 따라서, 맵이 갱신될 때, 이미지는 있고 R3가 아닌 S2의 탄젠트 공간에서 이동되어야 한다.
전술한 것들이 구축되면, 전술한 등각 구조를 연산하기 위해 여러 알고리즘이 하기에 제공된다. 컴퓨터 그래픽, 컴퓨터 비전, 메디컬 이미징 영역에서의 적용이 기술된다.
다음의 알고리즘인 알고리즘(1)에서, 임의의 종수가 0인 표면과 구면 사이의 등각 맵이 연산된다. 우선 이미지의 질량 중심이 연산되어야하고, 다음의 폼이다.
불연속인 경우, 다음의 근사치가 사용된다.
여기서, A[u, v, w]는 면[u, v, w]의 영역이다.
알고리즘(1)은 현재 S2에 대한 종수가 0인 메시의 등각 맵을 연산하는 데에 사용될 수 있다.
알고리즘(1): 종수가 0인 메시의 등각 패러미터화
입력 : 닫힌 종수가 0인 메시 M
출력 : 전체 등각 지도 φ: M→S2
1. M을 S2로 매핑한 가우스 맵을 연산.
2. M, 각각의 꼭지점에서의 라플라시안을 연산.
3. ∈S2의 탄젠트 공간으로의 의 사영.
4. 네거티브하게 사영된 를 따라서 를 갱신.
5. 의 질량중심 , S2의 중심에 대한 질향 중심을 연산하고, S2 상에 있는 를 재정규화(renormalize).
6. 사영된 라플라시안이 0이 될때까지 모든 꼭지점에 대한 꼭지점 스텝(2-5)을 반복.
도 3, 4, 및 5는 3 개의 다른 종수가 0인 표면에 대한 구체 등각 매핑을 도시한다. 특히, 도 3은 S2로 등각으로 매핑된 가고일 모양의 모델을 도시하고, 도 4는 S2로 등각으로 매핑된 뇌를 도시하며, 도 5는 S2로 등각으로 매핑된 토끼의 모델을 도시한다.
임의의 2 개의 위상의 디스크 사이의 등각 맵을 연산하는 것을 고려하여, 모든 그러한 매핑은 전술한 뫼비우스 그룹의 서브 그룹이고 다음의 식으로 나타내는 3차원 그룹을 형성한다.
위상의 디스크와 단위 디스크 사이의 등각 맵을 연산하기 위해, 더블링이라고 하는 기술이 사용된다.
더블링은 경계가 있는 표면을 닫힌 대칭형 표면으로 변환한다. 경계 을 갖는 표면 M이 주어지면, 대칭의 닫힌 면 이 구축되어 이 M을 2번 덮는다. 즉, 면 을 면 에 대해 등축으로 매핑하는 등축 사영(isometric projection) 이 존재한다. 각 면 에 대해, 에 2개의 원상(preimage)이 있다. 알고리즘(2)은 일반 메시 M의 더블링을 연산한다.
알고리즘(2). 열린 메시의 더블링 연산
입력 : 경계를 갖는, 메시 M
출력 : M, 의 더블링
1. -M으로 표시되는 M의 사본을 만든다.
2. -M의 배향을 역전한다.
3. 임의의 경계 꼭지점 에 대해, 고유하게 대응하는 경계 꼭지점이 존재하고, 의 에지에 대해, 고유한 경계 에지 가 존재한다. 대응하는 모든 꼭지점과 에지를 찾아라.
4. M과 을 결합하여 대응하는 꼭지점과 에지가 일치하도록하라. 그 결과인 메시는 을 더블링하는 것이다.
알고리듬 2에서 설명된 더블링 기술을 사용하면, S2에 대한 위상 디스크의 등각 맵핑은 직접 계산될 수 있다. 더블링 표면은 대칭이기 때문에, M 및 -M 은 분리된 반구로 맵핑되게 될 것이며, 입체 투영(stereographic projection) π를 사용하면, 구의 반구는 유닛 디스크로 맵핑될 수 있다. 이 방식에서, 위상 디스크와 유닛 디스크 D2 사이에 맵을 만드는 등각 맵핑은 계산된다. 수학식 52 에서의 뫼비우스 변환을 적용함으로써, 모든 가능한 등각 맵핑은 계산될 수 있다.
알고리듬 3. 위상 디스크로부터 D2 까지의 글로벌 등각 맵 계산.
입력 : 위상 디스크 M.
출력 : M으로부터 유닛 디스크 D2 까지의 글로벌 등각 맵 Φ.
1. M의 더블링 계산.
2. 대칭을 유지하면서, 글로벌 등각 맵 Φ: →S2 계산.
3. Φ(∂M)이 적도가 되게 Φ()회전.
4. 유닛 디스크에 대한 상위 반구를 맵핑하기 위해 스테레오-그래픽 투영 π사용.
5. °Φ출력.
논-제로 지너스를 갖는 표면에 대해, 표면의 위상에 의해 판정된 홀로모로픽 1-폼 그룹 1(M)은 이러한 표면에 대한 글로벌 등각 파라미터화를 계산하는데 중요하다. 이 그룹을 계산하기 위해서, 호몰로지 기저는 먼저 계산되고, 듀얼 조화 1-폼 기저는 그다음에 계산되며, 그다음 조화 1-폼은 기저 홀로모로픽 1-폼으로 변환된다.
호몰로지 및 조화 1-폼을 계산하는 대수적인 알고리듬이 도입된다. 메쉬 M이 제공되면, 상응하는 호몰로지 기저는 대수 위상 방법을 사용하여 계산된다. δi n∈Kn 및 δi n-1∈Kn-1 이 되게 하면,
을 정의한다.
그러면 n-차원 바운더리 행렬은
로서 정의된다.
그다음 호몰로지 기저는 다음의 오퍼레이터
의 제로 고유치에 상응하는 고유벡터로부터 형성된다.
알고리듬 4. 메쉬 M에 대한 호몰로지 기저 계산.
입력 : 메쉬 M.
출력 : 호몰로지 기저 { 1, 2,..., 2g}
1. ∂1, ∂2 에 대한 바운더리 행렬 계산.
2. 행렬 D = ∂1 T1 + ∂22 T 의 Smith 노말 폼 계산.
3. { 1, 2,..., 2g}를 형성하기 위해서, 제로 고유치에 상응하는 D의 고유벡터 구하기.
모든 조화 1-폼은 호몰로지 그룹 H1(M,Z)의 듀얼인 코호몰로지 그룹을 형성한다. 조화 1-폼은 닫혀지고 조화된다. Hodge 이론에 따르면 모든 조화 1-폼은 호몰로지 그룹의 듀얼 공간인 선형 공간을 형성한다. 또한, 각 코호몰로지 등급은 고유의 조화 1-폼을 갖는다.
알고리듬 5. 일 세트의 조화 1-폼 기저 계산.
입력 : M의 호몰로지 기저 { 1, 2,..., 2g}.
출력 : 조화 1-폼 기저 {ω1, ω2,...,ω2g}.
1. ci j = - i j, i,j = 1, 2,...,2g 의 값 설정.
2. ωi 에 대한 다음의 선형 시스템 해결.
δωi = 0
△ωi = 0
< ωi , j>= - i j
3. {ω1, ω2,...,ω2g} 출력
상기에서 사용된 대수적인 접근에 대한 대안으로서, 호몰로지, 코호몰로지, 및 조화 1-폼은 다음과 같은 알고리듬의 조합을 사용하여 계산될 수 있다.
알고리듬 6. 메쉬 M의 기초 도메인 계산.
입력 : 메쉬 M.
출력 : M의 기초 도메인 DM.
1. DM=f0, ∂DM = ∂f0 로 하고 임의의 페이스 f0 ∈ M 선택, 큐 Q로 에지와 f0 를 공유하는 f0 의 모든 이웃하는 페이스를 출력.
2. Q가 비어있지 않는 동안
a. ∂f = e0 + e1 + e2 로 하고, Q 내의 제 1 페이스 f 제거.
b. DM= DM∪f.
c. -ei ∈ ∂DM 이 그 차수를 유지하면서 ∂DM 내의 -ei 를 {ei+1 , ei+2}로 대체하는, 최초의 ei ∈ ∂f를 구한다.
d. 에지를 f 및 Q로의 DM 또는 Q 가 아닌 것을 공유하는 모든 이웃하는 페이스 입력.
3. 서로에 대해 대향하는 ∂DM 의 모든 이웃하는 방향의 에지를 제거, 즉 ∂DM 으로부터 모든 쌍 {ek , -ek} 제거.
결과적인 기초 도메인 DM 은 그들의 삽입 순서에 따라 분류되는 M의 모든 페이스를 포함한다. DM 의 최종 바운더리의 향하는 방향이 없는 에지 및 정점은 컷(cut) 그래프로서 참조되는 그래프 G를 형성한다.
컷 그래프에 대해, 알고리듬 7은 또한 M의 호몰로지 기저인 상응하는 호몰로지 제너레이터를 계산한다.
알고리듬 7. M의 호몰로지 기저 계산.
입력 : 메쉬 M.
출력 : 호몰로지 기저 .
1. M의 기초 도메인 DM 계산 및 상응하는 커트 그래프 G 판정.
2. G의 확장 트리 T 계산, G/T = { e1, e2,..., e2g}라 하자.
3. 루트(root) 정점 r∈T, 깊은 제 1 트래버스(traverse) T 선택.
4. ∂ei= ti - si 로 하고, 루트 r로부터 [r, ti] 및 [r, si] 로서 나타나는 ti 및 si로의 경로가 존재하며, 그다음 그들을 루프 i = [r, ti]-[r, si] 에 연결.
5. H1(G,Z) 및 H1(M,Z)의 기저로서 { 1, 2,..., 2g} 출력.
M, H1(M,Z)의 코호몰로지 그룹에 대한 기저를 명백하게 계산하기 위해서, 일 세트의 닫혀진 1-폼 {ω1, ω2,...,ω2g}는 다음과 같이 구해진다.
여기서, δi j 는 Kronecker 델타 이고, I 는 호몰로지 기저임.
알고리듬 8. M의 코호몰리지 기저 계산.
입력 : 메쉬 M.
출력 : 코호몰로지 기저 .
1. 기초 도메인 DM 및 메쉬 M의 커트 그래프 G 계산, 및 확장 트리 T, G/T = { e1, e2,..., e2g} 계산.
2. 임의의 에지 e∈T 에 대해 ωi(ei)=1 및 ωi(e)=0 로 함.
3. DM은 DM={f1, f2,...,fn}인 방식으로 정렬된다면, DM의 차수를 {fn, fn-1,...,f1}으로 반대로 전환.
4. DM이 비어있지 않는 동안,
a. DM의 제 1 페이스 유도, DM이, ∂f = e0 + e1 + e2로부터 f 제거.
b. {ek}를 두개의 세트, = {e ∈ ∂f | -e ∈ ∂DM}, = {e ∈ ∂f | -e ∂DM}로 나눔.
c. ωi(ek),ek가 비었다면, 우변이 영과 동일한, ωi(e) = ωi(e) 인, 의 값을 임의로 선택.
d. ∂DM = ∂DM + ∂f 라 하고, DM의 바운더리 업데이트.
일단 코호몰로지 기저 {ω1, ω2,...,ω2g} 가 계산되면, 호몰로지 기저 { 1, 2,..., 2g}의 듀얼은
와 같이 선형 변환 {ω1, ω2,...,ω2g} 에 의해 구해질 수 있다.
알고리듬 9. 닫혀진 1-폼을 조화 1-폼으로 확산.
입력 : 메쉬 M, 닫혀진 1-폼 ω.
출력 : ω에 대해 코호몰로지한 조화 1-폼,
1. △(ω+ ∂f) ≡0 과 같이 f ∈ C0(M)선택.
2. f에 대한 상기 희소 선형 시스템(sparse linear system) 해결.
3. ω+ ∂f 출력.
여기서,
조화 1-폼 {ω1, ω2,...,ω2g}이 제공되면, 켤레 조화 1-폼 ω* 다음의 선형 시스템을 해결함으로써 구해질 수 있다.
일단 기초 도메인이 계산되면, 등각 맵핑은 홀로모로픽 1-폼 ω을 적분함으로써 직접 계산될 수 있다. 먼저, 루트 정점 v0∈DM 선택하고 그다음 DM 을 트래버스하기 위해서 깊은 제 1 탐색 방법을 사용하자. 각 정점 u∈DM 은 v0로부터 u까지의 고유한 경로를 갖고, 그다음 (u) =< ω,> 를 정의한다.
알고리듬 10. 메쉬 M의 글로벌 등각 파라미터화.
입력 : 메쉬 M, 조화 1-폼 ω.
출력 : 맵 , 또는 글로벌 등각 파라미터화.
1. 기초 도메인 M의 DM 계산.
2. 정점 u∈DM 를 트래버스하기 위해서 깊은 제 1 탐색 방법을 사용하고, u 로서 표시되는 루트 정점 v0로부터 u까지의 경로 기록.
3. 적분 (u) =< ω, u > 계산.
4. u의 등각 좌표로서 (u) 출력.
알고리듬 11. 메쉬 M의 등각 구조.
입력 : 메쉬 M.
출력 : M{(Ui, zi)}의 등각 구조.
1. 홀로모로픽 1-폼 기저 계산.
2. M⊂Ui, Ui가 단순하게 연결되는, 파티션 {Ui} 계산.
3. 각각의 Ui에 대해 홀로모로픽 베이스 선택, Ui 상의 홀로모로픽 1-폼 집적, zi로서 맵핑 표시. 제로 포인트가 존재한다면, Ui 세분 및 단계 3 반복.
4. {(Ui, zi)} 출력.
기초 도메인 상의 홀로모로픽 1-폼을 집적함으로서 얻어지는 글로벌 등각 파라미터화는 메쉬의 정준 분해(canonical decomposition)를 위해 사용될 수 있고, 메쉬를 텐서 프로덕트 스플라인 표면(tensor product spline surface), 표면 매칭 및 인식, 및 다른 유용한 이미지 프로세싱 적용으로 변환시킨다.
Poincare-Hopf 이론에 따르면, 홀로모로픽 1-폼 ω는 토러스(torus)로의 위상동형(homeomorphic)이 아니라면 제로 포인트를 가져야만 한다. ω의 제로 포인트는 등각 인자가 제로인 포인트이다. 지너스-g 표면은 2g-2 제로 포인트를 갖는다. 등각 맵핑은 각 포인트의 이웃을 두배 및 이중으로 싸곰 복소 평면 상의 이미지 p의 이웃을 커버한다. 맵에 국부적으로, : C →C는 이웃에서
와 유사하다.
도 6A 및 6B 는 개방 티폿(teapot) 모델 및 복소 평면 각각에 대해 글로벌 등각 파라미터화 상의 제로 포인트를 묘사한다.
우리는 표면 M으로부터 단위 원 S1 까지의 맵핑으로서 조화 1-폼 ω를 취급할 수 있다. 그러면 홀로모로픽 1-폼에 대해, 실질적인 부분의 조화 1-폼은 원 평가된(circle valued) 맵핑이다. 가상 부분의 조화 1-폼은 그레디언트 필드(grident field)이다. 제로 포인트를 통하는 적분 곡선은 표면을 균등한 패치로 세분할 것이다. 특히, 메쉬 M 및 홀로모로픽 1-폼 ω=+ *에 대해, 또는 *를 따라 및 제로 포인트를 통하는 적분 곡선은 표면을 위상 디스크 또는 실린더로 분할한다.
M이 C에 대해 등각으로 맵핑된 위상 토러스 M이 되게 하자. 보편적인 커버링 공간 상에서 홀로모로픽 1-폼 ω을 집적함으로써, 주기적인 등각 맵이 결과로서 산출된다. 베이스 포인트 u0를 선택하면, 베이스 포인트의 이미지 세트는 다음과 같다.
이러한 맵핑은 주기적이거나 또는 모듈러(modular)하다. 전체 토러스는 M의 주기로서 참조되는, < 1, ω>, < 2, ω> 에 의해 확장된 평행사면형인 일 주기로 맵핑된다. M의 지너스-g 가 1 보다 크다면, 다른 핸들들은 다른 주기를 갖을 수 있다. 전체 표면은 g 가 덮여진 모듈러 평행사면형에 맵핑된다. 평행사면형은 부착될 수 있고 제로 포인트의 이미지를 통해 서로 교차할 수 있다.
도 7A-7D는 이러한 현상을 나타낸다. 도 7A 및 7B에서, 2-홀 토러스는 두개의 핸들로 분리되고 각 핸들은 모듈러 공간으로 등각으로 맵핑된다. 도 7C 및 7D는 지너스-3 토러스 및 모듈러 공간으로의 등각 맵핑을 나타낸다.
여기에서 설명된 방법을 일반화하기 위해서, 바운더리를 갖는 메쉬는 현재 고려될 것이다. 바운더리를 갖는 메쉬 M이 주어지면, M의 더블링 은 계산된다. 각 내부 정점 u∈M에 대해, u1 및 u2 로서 나타나는 에서의 u의 두개의 카피가 존재한다. u1 및 u2 는 아래와 같은 것으로서 서로에게 듀얼이다.
각 바운더리 정점 u∈∂M에 대해, 에서 단지 하나의 카피가 존재하여, u는 그 자체에 대해 듀얼이된다.
M 상의 조화 1-폼을 계산하기 위해서, 의 모든 대칭적인 조화 1-폼은 또한 M 상의 조화 1-폼이다.
각 조화 1-폼 ω에 대한 듀얼 오퍼레이터을 아래로 정의하자.
임의의 ω는 대칭 부분 및 아래와 같은 비대칭 부분으로 분해될 수 있다.
여기서, 는 대칭 부분이고, 는 비대칭 부분임.
알고리듬 12. 바운더리를 갖는 메쉬에 대해 일 세트의 홀로모로픽 1-폼 기저 계산.
입력 : 바운더리를 갖는 메쉬 M.
출력 : 폼 { 1 + 1 *, 2 + 2 *,..., k + k *}의 메쉬 M에 대한 홀로모로픽 1-폼.
1. M의 더블링, 계산.
2. 의 조화 1-폼 기저 계산.
3. i = 할당, 중복되는(redundant) 것들 제거.
4. i *로서 표시되는 i 의 켤레 조화 1-폼 계산.
5. 홀로모로픽 기저 { 1 + 1 *, 2 + 2 *,..., k + k *} 출력.
도 8A 및 8C는 비록 그들은 위상학적으로 동일한, 즉 두개의 표면이 등각적으로 동일하지 않은 두개의 지너스-1 표면이지만, 두개의 지너스-1 표면을 도시한다. 각 토러스는 개방되게 커트될 수 있고 도 8B 및 8D에 각각 나타난 것과 같은 2차원 평행사변형에 등각으로 맵핑될 수 있다. 각 평행사변형의 모양은 등각의 동일한 종류를 나타낸다. 등각의 동일한 종류들은 평행사변형의 날카로운 각, 이 두개의 경우에서의 우측 각, 및 등각의 불변을 나타내는 두개의 이웃하는 에지 사이의 길이 비율, 또는 이러한 두개의 지너스-1 표면의 모양 인자들에 의해 판정된다. 도 8B 및 8D에 나타난 것과 같이, 두개의 토러스는 다른 모양 인자를 갖고 등각적으로 동등하지 않다.
표 1은 이하에서 도 9A-9D에 나타난 지너스-1 표면의 등각의 불변을 포함한다. 도 9A-9D에 의해 나타난 표면의 어떠한 것도 등각적으로 동일하지 않음은 명백하다.
지너스-1 표면의 등각의 동일한 것들
메쉬 모양 인자 정점 페이스
토러스 1.0-1.142i 1089 2048
노트(knot) 1.0-.272i 5808 11616
노트2 1.0+0.128i 2050 3672
락커 1.0-3.509i 3750 7500
티폿 1.0-.112i 17024 34048
알고리듬 13. M1 및 M2가 등각적으로 동등함을 증명
입력 : 두개의 메쉬 M1 및 M2.
출력 : M1 및 M2의 등각적으로 동등하거나 또는 동등하지 않은 표시
1. 각각 M1 및 M2 에 상응하는 주기 행렬 (R1, C1) 및 (R2, C2) 계산.
2. R1 = P1 1P1 -1 및 R1 = P2 2P2 -1의 Jordan 노말(normal) 폼 계산.
3. 1 2 라면 거짓으로 복귀.
4. N = P1P2 -1 로하고, N이 역이 가능한 정수 행렬 및 NC1Nt = C2 라면 참으로 복귀, 그렇지 않으면 거짓으로 복귀.
등각 인자 λ(u, v)는 표면 S의 제 1 기초 폼을 나타낸다. 만약 λ가 상수라면 표면의 Gaussian 곡률은 제로다. 표면 상에서 선택적으로 커팅함으로써, 새로운 바운더리는 도입되고, 따라서 등각의 구조는 변경될 수 있다. 실제로, 파라미터화의 일치를 개선하는 것이 도움이 되고, 일반적으로 이러한 커트는 높은 Gaussian 곡률을 갖는 표면의 영역 상에서 이루어진다. 도 10A-10D는 일치의 개선을 나타낸다. 도 10A에 나타난 구의 파라미터화에서, 이어(ear) 부분은 크게 아래로 샘플링된다. 이어 팁(ear tip)에서 위상 커트를 도입함으로써, 파라미터화는 더욱 많이 일치하게 된다.
일반적으로 계산의 안정성은 삼각형 분할의 품질에 크게 의존한다. 삼각형 분할의 모든 각이 날카로운 각도라면, 계산 알고리듬은 안정되고 수렴되도록 보장된다. 도 15는 표면 모델 복잡성의 두가지 다른 레벨에서의 티폿 모델의 글로벌 파라미터화를 나타낸다. 도 15A-15B에 나타난 바와 같이, 많은 복잡한 본래의 티폿에 대해 글로벌 파라미터화는 모든 각도가 날카로운 각도 특히 정확한 각도로 귀착된다. 도 15C-15D는 모든 각도가 날카로운 각도 특히 정확한 각도인 단순화된 티폿의 글로벌 등각 파라미터화를 나타낸다. 양자의 경우에, 모델의 복잡성과 관계없이, 계산 알고리듬은 수렴되고 안정적이다. 다음의 알고리듬은 모두 날카로운 각도를 갖는 삼각형 분할에 근접한다.
알고리듬 14. 모두 날카로운 각도를 갖는 표면의 삼각형 분할.
입력 : 메쉬 M.
출력 : 모두 날카로운 각도를 갖는 리메쉬(remesh) M
1. 루프 세분 방법을 사용하여 메쉬 세분.
2. 최소 에지 길이 표준을 사용하는 선분병합(Edge Collapse)을 사용하여 메쉬 단순화.
3. M 상의 모든 각도가 날카로울때 까지 단계 1 및 2 반복.
4. 리메쉬 M 출력.
등각 패러미터(conformal parameter)에 기초한 표면 매칭과 평균 곡선 매칭
인간의 표정(human expression)이나 피부의 변형과 같이, 하나의 표면이 너무 많은 왜곡없이 다른 하나로 변형될 수 있다면, 상기 변형은 글로벌 등각 매핑에 의해 정확하게 시늉내어질 수 있다. 등각 패러미터화는 표면의 제 1 기본 폼에 의존하고, 특히 등각 구조는 리만 계량 텐서(Riemannian metric tensor)가 너무 많이 변화하지 않고 상기 등각 구조가 유사한 한 연속한 리만거리(Riemannian metric)에 의존한다. 따라서, 정준 패러미터(canonical parameter) 영역으로의 두 표면의 매핑과 상기 패러미터 영역에서의 표면의 매핑은 3-D 매핑 문제가 보다 효율적으로 해결되도록 한다.
패러미터 영역에 등각 팩터 λ(u,v)와 노멀 n(u,v)를 저장함으로써, 원표면이 R3에서 로테이션과 전환에 따라 고유하게 재구성될 수 있다. λ(u,v)는 제 1 기본 폼을 정의하고, n(u,v)는 제 3 기본 폼을 정의하여 그에 따라 제 2 기본 폼을 정의하여, 즉, R3에 임베딩한 것이 연산될 수 있다. 따라서, 상기 표면은 유클리드 변환에 따라 고유하게 구성될 수 있다.
보다 효율적인 방법은 상기 등각 패러미터 영역 상의 평균 곡률을 사용하는 것이다. 경계를 갖지 않는 임의의 닫힌 표면에 대해, 상기 표면은 등각 팩터 λ(u,v)와 평균 곡률 H에 의해 고유하게 결정된다. 경계를 갖는 임의의 열린 표면에 대해, 표면은 등각 팩터 λ(u,v)와 평균 곡률 H, 및 경계 상의 제 2 기본 폼에 의해 고유하게 결정된다.
가우시안 곡률과 평균 곡률에 기초한 표면을 매치하기 위해, 매치될 표면은 정준 패러미터 영역에 포함된다. 예를 들면, 인간의 얼굴은 단위 디스크로 매핑될 수 있다(mapped). 가우시안 곡률과 평균 곡률은 등각 패러미터화(parameterization)를 사용하여 연산된다. 가우시안 곡률과 평균 곡률의 레벨의 세트는 패러미터 영역 상의 평면 곡선의 패밀리이다. 이 곡선의 레벨의 세트는 상기 표면들을 매칭하는 데에 사용된다.
특별한 특성을 가진 표면을 매칭하기 위해, 상기 특성의 포인트는 우선 제거되고 상기 표면의 더블링(doubling)이 연산된다. 다음으로, 상기 맵의 호모토피(homotopy) 형은 상기 제 1 표면에서의 특성이 제 2 표면의 대응하는 특성에 매치가 되는것을 보장하도록 제약을 가한다. 상기 등각 구조는 그런다음 전술한 바와 같이 매칭을 수행하도록 연산된다. 예를 들면, 인간의 얼굴을 매칭하기 위해, 눈, 코끝, 및 입과 같은 특성이 등각 구조의 연산에 앞서서 제거된다.
표면 분류(Surface Classification)
효과적인 데이터베이스화와 검색을 할 수 있는 표면 분류를 위해, 각 표면에 대한 주기 행렬(period matrix)의 형태로 상기 등각 구조가 연산되고 저장된다. 도 11A 내지 11D는 여러 종의 2 개의 표면을 도시한다. 하기에서 볼수 있는 바와 같이, 상기 주기 행렬 R이 등적이 아닌것 처럼 도 11A 내지 11D에 도시된 표면 중 어떤 것도 등적이 아니다.
도 11A의 2 홀의 원환체(torus)는 861 개의 꼭지점과 1,536 개의 면을 포함하고,
의 주기 행렬 R을 갖는다.
도 11B에 도시된 꽃병의 모델은 1,582의 꼭지점과 2,956의 면과,
의 주기 행렬 R을 갖는다.
도 11C에 도시된 꽃의 모델은 5,112의 꼭지점과 10,000의 면과,
의 주기 행렬 R을 갖는다.
도 11D에 도시된 옹이가 있는 병은 15,000의 꼭지점과 30,000의 면과,
의 주기 행렬 R을 갖는다.
표면 인식(Surface Recognition)
서로 직접적인 매칭 없이 표면이 인식될 수 있는 것이 바람직하다. 정준 방식으로 표면의 등각 구조를 변형하고 각 변형에 대해 주기 행렬을 연산하는 것은 상기 표면의 고유한 기하학적인 특성을 나타내는 일련의 주기 행렬을 제공한다.
예를 들면, 인간의 얼굴을 인식하기 위해, 왼쪽 눈의 중심, 오른쪽 눈의 중심, 코끝, 입의 중심과 같은 특성 포인트가 제거된다. 얼굴 표면에 대한 각 변형에 대해, 표면의 더블링과 주기 행렬이 연산된다. 일련의 주기 행렬을 비교함으로써, 우리는 얼굴과 같은 기하학적인 표면을 인식할 수 있다.
또는, 모든 주된 특성 포인트가 제거되고, 상기 표면 내의 다른 포인트가 선택되고 이동되어, 상기 더블링의 주기 행렬이 상기 선택된 포인트의 각 이동에 대해 연산된다. 예를 들면 인간의 얼굴을 인식하기 위해 눈의 중심, 코끝, 및 입의 중심에서의 포인트가 제거되고, 얼굴의 다른 포인트가 미리 정해진 궤도를 따라 이동된한다. 각 스텝에서, 현재 위치의 포인트가 제거되고 주기 행렬이 연산된다. 미리 정해진 궤도를 따라 각 포인트에 대해 일련의 주기 행렬이 연산된다. 표면을 인식하기 위해 이 주기 행렬들이 사용된다.
조화 스펙트럼 분석
또는 전술한 라플라스 연산자는 부정 고유값과 고유함수를 갖는다. 모든 고유값의 스펙트럼은 상기 표면의 고유한 기하학적 형태의 많은 것을 나타낸다. 또한, 고유함수는 상기 표면을 재구성하는 데에 사용될 수 있다. 리(Rhe) 표면은 오직 표면의 서명(signature)으로서 표면의 스펙트럼을 사용해서만이 인식될 수 있다. 예를 들면, 의학 영역에서, 내장의 형태의 스펙트럼을 분석함으로써 일련의 병이 검지될 수 있다.
원하는 고유값과 고유함수가 라플라스 행렬의 고유값과 고유벡터를 찾아냄으로써 삼각형 메시에 의해 나타내는 표면에 대해 연산될 수 있다.
조화 고유함수를 사용하는 표면 데이터의 압축
종수가 0인 표면은 단위 구에 대해 등각적으로 매핑되고, 상기 표면의 위치 벡터는 상기 구면 상에 정의된 벡터값을 가진 함수로서 나타낸다. 상기 구면 상의 라플라스 연산자의 고유함수들은 상기 구면의 함수 공간에 대한 기저를 형성하는 구면 조화함수이다. 상기 위치 벡터는 그런 다음 상기 함수의 기저에 대해 분석되고 스펙트럼이 얻어진다. 고주파 성분을 필터링 함으로써, 표면 데이터가 압축된다. 전술한 뫼비우스 변환을 사용하는 동안, 영역은 추가적인 조사를 위해 "확대(zoomed into)" 될 수 있다. 일반적인 표면에 대해, 라플라스 연산자의 고유함수를 사용하여 그의 등각 동등한 클래스에서의 정준 형태로 상기 표면을 등각적으로 매핑하고 표면 위치 벡터를 분석하는 것은 고주파수 성분이 저장에 앞서 제거될 수 있는 원하는 함수의 기저를 제공한다.
또한, 상기 등각 죄표계에서 정의된 등각 팩터와 평균 곡률은 유클리드 변환에 대해 고유하게 상기 표면을 결정하는데에 사용될 수 있다. 이러한 방법으로, 평면 상에 정의된 2 개의 함수, 즉 등각 팩터와 평균 곡률이 표면을 나타내는 데에 사용된다. 따라서, 3분의 1의 저장공간의 절약이 실현된다. 또한 압축은 전술한 고유함수의 기술을 사용하거나 다른 알려진 압축 기술을 사용함으로써 얻어질 수 있다.
리메슁(remeshing)과 하드웨어 설계
등각 구조를 사용함으로써, 우리는 표면이 패러미터 영역으로 등각적으로 매핑된 후에 상기 표면을 리메시(remesh)할 수 있다. 이러한 방법으로, 불규칙한 접속이 정삼각형의 분할(regular triangulation)로 변경될 수 있다. 이론적으로, 재구성된 노멀(normal)들은 정확하다. 이것은 기하학적 데이터의 표시를 단순화시키고, 그래픽 하드웨어 아키텍처를 단순화 시킨다. 현재, 일반적인 그래픽 하드웨어에 대해, 접속정보를 저장하기 위한 메모리 버퍼가 있다. CPU와 이 접속정보를 특정짓기 위해 필요한 그래픽 카드 사이의 통신은 극도로 시간 소모적이다. 그래픽 카드 상의 메모리 내에 저장된 데이터의 접속이 정상적이고, 상기 그래픽 카드가 그 스스로 그것을 예측할 수 있다면, 그때에는 접속 정보를 위한 추가적인 메모리가 필요가 없다. 따라서, 프로세서와 그래픽 카드 사이의 통신의 필수 레벨을 감소 시킨다. 그래픽 카드의 아키텍처에 대하여, 현재 기하학적 배열을 처리하는 파이프라인과 표면 텍스처를 처리하는 파이프라인이 개별적이다. 정상적인 접속을 사용하여, 기하학적 배열이 텍스처로서 또한 나타낼 수 있고, 이 두개의 개별적인 파이프라인들은 결합될 수 있다. 이런 방식으로 그래픽 카드 아키텍처의 복잡성이 감소될 수 있다.
또한 리메시에 의해, 기하학적인 이미지가 구성될 수 있고, 상기 이미지 포맷이 표면의 기하학적 배열을 나타내기 위해 사용될 수 있다. 이러한 방식으로, 다른 것들 중에서 압축, 다중 분해능(multi-resolution), 및 필터링과 같은 기하학적 배열상에서 동작하는 많은 이미지 처리 기술이 사용될 수 있다.
도 12A는 원 메시의 비정상적인 연결을 갖는 토끼 모델을 도시한다. 등각 구조를 사용하여 리메시를 한 후, 도 12B에 도시된 것과 같이 접속이 매우 정상적으로 되고, 재구성된 노멀들은 매우 정확해진다. 상기 등각 기하학적 이미지가 도 12C에 도시되고, 재구성된 형태가 도 12D에 도시된다.
패러매트릭 표면과 메시 컨버전(Parametric Surface and Mesh Conversion)
CAGD 필드에서, 비스플라인(BSpline) 표면과 Bezier 표면과 같은 패러메트릭 표면이 자주 사용된다. 제조업에서, 이러한 종류의 패러메트릭 표면을 사용하는 컨트롤러가 자주 처리 머신을 가이드한다. 그러나, 기하학적 데이터는 자주 삼각형 메시로 나타낸다. 현재의 기하학적 데이터 획득 디바이스는 짙은 포인트 구름으로 기하학적인 데이터를 출력한다. 이 스캔된 포인트 구름을 메시로 변환하는 것이 보다 쉽고, 메시로부터 패러메트릭 표면을 변환하고 메시로 패러메트릭 표면을 변환하는 것은 매우 중요하다. 현재, 메시를 스플라인 표면으로 자동으로 변환하는 방법은 없다.
본문에 기술된 등각 기하학적인 기술을 사용하여, 이 문제는 해결될 수 있다. 전술된 것과 같이, 표면의 글로벌 등각 패러미터화가 연산되고, 상기 표면은 제로 포인트를 통해 그래디언트 필드를 따라 집적 라인을 사용하여 정준 패치로 분해된다. 각 정준 패치는 평면상의 사각형으로 매핑되고, 텐서 프로덕트 스플라인 표면이 그 위에 구축된다. 상기 결과로서 나타나는 패러미터화는 바운드리 상에 제어 포인트를 매칭함으로서 전체적으로 부드럽게 될 수 있다. 따라서, 상기 메시를 임의의 바람직한 연속체를 가진 패러메트릭 표면으로 변환하는 것이 편리하다. 추가로, 이 구성은 정확한 노멀 정보를 유지한다.
표면에서의 수치계산
등각 구조는 표면에서의 공변미분(covariant differentiation)을 연산하기에 좋은 패러미터화이다. 공변미분은 표면 기하학적 배열에 대해 고유한 것이어서, 유클리드 표면에서의 삽입은 무관련하다. 등각 구조 분석은 왜곡된 형태의 표면에서의 자연적 물리적 프로세스를 연산하는 포텐셜을 갖는다.
등각좌표계를 사용하여, 미분연산자는 매우 간단한 포맷을 갖는다. 예를 들면, 라플라스 연산자는
이다.
이 기술은 나비어-스토크스 방정식 및 맥스웰 방정식과 같은 표면의 편미분 방정식의 보다 쉬운 해답을 갖도록 한다. 전술한 등각 구조를 사용하여, 표면의 가우시안 곡률이 쉽게 결정된다.
메디컬 이미징
전술한 등각 구조는 또한 뇌의 매핑, 뇌의 표시(brain registration), 심장표면의 매칭, 및 혈관의 표면 분석과 같은 메디컬 이미징 영역에 적용될 수 있다.
예를 들면, 단위구로 뇌의 표면을 매핑함으로써, 2 개의 뇌를 비교하고 그 특성들을 비교하는 것이 편리하다. 뇌의 기하학적인 구조를 분석함으로써, 뇌의 오버타임의 변화와 잠재적인 질병을 보다 쉽게 발견할 수 있다.
뇌의 표면에서 구까지 등각 맵은 삼각측량과 분해능에 독립적이다. 등각 매핑은 우리가 2 개의 뇌의 표면을 표시하고 비교할 수 있는 훌륭한 표준 공간을 제공한다. 뇌의 표면은 매우 복잡하기 때문에, 다른 방법으로 꼭지점(vertex)의 흐름의 전개를 추적하는 것은 매우 어렵다. 본문에 기술된 방법은 정확한 각도 정보를 관리하면서 복잡한 표면 구조를 제어한다. 뇌는 일반적으로 종수가 0인(genus-zero) 표면이기 때문에, 전술한 알고리즘(1)이 단위구에 대해 뇌표면을 매핑하는데에 사용된다. 도 14는 뇌 매핑에 대한 예를 도시한다.
애니메이션
등각 기하학적 배열은 또한 컴퓨터 그래픽 애니메이션에 적용될 수 있다. 현재의 데이터 획득 기술을 사용하여, 배우의 3D 형태를 다른 제스처와 표현으로 스캔할 수 있다. 전술한 등각 분석 기술을 사용하면, 이 키 제스처와 표현들이 서로 매핑될 수 있다. 스플라인 삽입 기술을 사용하여, 제스처와 표현 사이의 부드러운 변화가 그것들 사이에 생성될 수 있다. 따라서, 부드러운 형태와 왜곡된 모델을 포함한 임의의 형태가 동화상화 될 수 있고, 이것은 현재 방법들을 사용하여 동화상을 만드는 것을 극도로 어렵게 한다.
우리에게 2개의 유사한 형태가 주어졌다고 가정하자. 제 1 특성 포인트가 위치지정되고 그런다음 제거된다. 상기 표면의 더블링은 연산되고, 매핑의 호모토피 유형이 결정된다. 각 표면의 홀로모픽(holomorphic) 1 폼이 선택되어, 2 개 표면의 코호모로지 유형이 호모토미 유형의 매핑에 의해 결정된다. 제로 포인트가 위치지정되고, 표면들이 제로 포인트를 지나는 그래디언트 선을 사용하여 패치되도록 분해된다. 각 패치들은 패러미터 영역의 사각형에 등각으로 매핑된다. 표면 사이의 맵을 얻기 위해, 평면 상의 이 패치들이 매칭된다.
키 형태 사이의 매핑이 알려지면, 키 형태상의 포인트들이 제어 포인트로서 작용하도록 선택된다. 비-스플라인이 키 형태들 사이에서의 부드러운 변화를 생성하도록 사용된다. 이것은 등각 구조를 사용하여 인간 여성의 얼굴이 인간 남성의 얼굴로 모핑되는(morphed) 도 15에 도시되어 있다. 이 방법으로, 우리는 임의의 형태를 동화상화 할 수 있다. 이것은 특히 인간 배우에게 유용하다. 다른 연령대의 배우의 얼굴의 표현, 제스처, 및 변형이 데이터 베이스에 저장될 수 있다. 이 저장된 기하학적 데이터는 가상의 배우를 형성해내도록 동화상화 될 수 있다.
왜곡없는 텍스처 매핑
표면의 텍스처 매핑은 컴퓨터 게임 산업과 영화 산업에 있어서 매우 중요하다. 표면의 렌더링 속도는 여러 다른 팩터 중에서 표시되는 기하학적 모델의 복잡도에 의해 결정된다. 컴퓨터 게임과 같이 리얼타임 적용을 위해, 단순 모델이 일반적으로 바람직하다. 이미지의 인식 화질을 개선하기 위해, 이미지들이 텍스처 매핑이라고 하는 프로세스를 사용하여 기하학적 표면 상에 붙여진다.
곡선 표면에 대해, 텍스처 매핑은 표시되는 이미지에서의 일종의 왜곡을 도입한다. 텍스처 왜곡의 가장 힘든 일은 평면과 곡선 표면의 텍스처 사이에서의 왜곡을 피하는 것이다. 산업에서, 기하학적 모델러와 텍스처 디자이너들은 일반적으로 다른 전문적 기술을 가진 다른 전문가이다. 텍스처 매핑은 기하학적 배열과 텍스처 모두를 변형할 필요가 있기 때문에, 이 두 다른 기술 세트사이의 협력은 대개 어렵고 시간 소모적이다.
전술한 바와 같이, 등각 패러미터화는 국부적인 왜곡을 갖지 않는다. 전술한 기술을 사용하여, 기하학적인 모델러와 텍스처 디자이너는 이전보다 쉽고 더 효율적으로 그들의 기술을 통합할 수 있다.
딜리크리트(Dirichlete) 방법을 사용한 텍스처 합성
텍스처 합성은 작은 텍스처 샘플로부터 주어진 표면을 덮을 수 있는 텍스처를 생성하는 것을 목적으로 한다. 이것은 그래픽 디자인, 영화 산업, 및 컴퓨터 게임 산업에서 중요하게 고려된다.
등각 패러미터화를 사용하여, 기하학적 표면 상의 텍스처 합성의 난제가 평면 상의 텍스처라는 보다 쉬운 문제로 변환될 수 있다. 전술한 등각 팩터 분석과 기술을 사용하여, 표면 상에서 텍스처를 늘리는 것이 제어될 수 있고, 표면 상의 텍스처의 기하학적 특성이 정확하게 예측될 수 있다.
전체적으로 부드럽게 합성된 텍스처를 만들기 위해, 텍스처의 패치의 경계를 확산시키는 딜리크리트 방법을 사용한다. 이것은 텍스처를 보다 부드럽게 그리고 부드럽게 만든다. 우선 제어된 스트레칭 효과를 가진 패러미터 평면 상의 개별적으로 떨어진 텍스처가 결정된다. 이 패치들은 그들 각각의 경계가 만나지만 서로 겹치지 않을때까지 성장된다. 이 패치들의 경계는 고정되고 딜리크리트 문제는 표면 상의 덮여지지 않은 영역에서 해결된다. 각 컬러 채널들은 함수로써 취급되고, 전술한 기술들이 표면상의 전체적인 부드러운 텍스처를 제공하는 해결책을 제공한다.
부피 조화 매핑(Volumetric Harmonic Mapping)
주어진 3D 복사본, M, f: M →R3은 조화 에너지를 최소화하는 것이 바람직하다. 이러한 방식으로, 원 3D 사본의 부피 매핑이 정준 공간에서 고려될 수 있다. f: M →R3에 대한 조화에너지는 다음의
과 같이 정의된다.
불연속 계에서, 조화 에너지는
로 정의된다.
여기서, 에서의 는 주어진 에지에 대해 반대의 2 면각이고, 1이 에지의 길이이다.
켤레 그래디언트 방법이 그런 다음 조화 매핑을 얻기 위해 상기 조화 에너지를 최소화하는 데에 사용될 수 있다. 부피 조화 맵은 구상으로 종수가 0인 3D 개체를 매핑하는 것이 발견될 수 있다. 구상의 정준 원에 대해, 종수 0인 개체에서의 닫혀진 단순한 곡선이 발견될 수 있다. 곡선 상의 플래토우(Plateau) 문제가 등각의 변형된 거리에 대해 해결될 수 있다. 이러한 방법으로, 표면에 의해 닫혀진 정준 작도(canonical description)가 얻어질 수 있다.
조화 매핑은 또한 외과수술의 시뮬레이션과 계획을 세우는데에 유용한 도구이다. 외과의사는 인체의 관심영역의 하나 이상의 MRI 이미지로부터 3D 뇌 부피 모델을 구성할 수 있다. 상기 외과 의사는 인체의 관심영역의 3D 지도를 구축할 수 있고, 새로운 환자의 인체의 관심영역의 3D 부피 데이터를 존재하는 지도 데이터와 비교할 수 있다. 조화 매핑은 고유하기 때문에, 이 기술은 뇌의 부피 데이터를 표시하는 유용한 방법이고, 외과 수술 시뮬레이션을 발달시키는 데에 유용할 것이다.
당업자는 또한 전술한 방법에 대한 변형과 변경이 본문에 기술된 발명의 개념에서 벗어남없이 이루어질 수 있다는 것을 이해할 것이다. 따라서, 본 발명은 첨부된 청구범위에 의해서 그 영역과 취지가 제한된다.
관련 출원에 대한 상호 참조
본 출원은 미국특허법(35 U.S.C.) 제119조(e)에 의거 2002년 11월 6일자 출원된 가출원번호 제60/424,141호에 대해 우선권을 주장하며, 거기에 기재된 사항은 본 명세서에 참조를 위해 병합되는 것으로 한다.

Claims (1)

  1. 제 1 및 제 2 표면을 매칭시키는 방법에 있어서,
    상기 제 1 및 제 2 표면의 제 1 및 제 2 메시(mesh) 표현을 각각 수신하는 단계;
    상기 제 1 및 제 2 메시 표현을 제 1 및 제 2 매핑된 표면을 형성하는 제 1 기준(canonical) 파라미터 도메인에 등각으로(conformally) 매핑시키는 단계;
    상기 제 1 및 제 2 매핑된 표면의 제 1 및 제 2 등각(conformal) 파라미터 표현(parameterizations)을 각각 계산하는 단계;
    상기 제 1 및 제 2 매핑된 표면에 대한 가우스 곡률과 평균 곡률의 제 1 및 제 2 레벨 세트들을 각각 계산하는 단계 - 여기서, 상기 가우스 곡률 및 평균 곡률의 제 1 및 제 2 레벨 세트들은 각각 상기 제 1 및 제 2 등각 파라미터 표현의 함수임 - ; 및
    상기 가우스 곡률 및 평균 곡률의 제 1 및 제 2 레벨 세트들을 비교하고, 상기 비교가 소정의 문턱값을 초과하는 경우에는 상기 제 1 및 제 2 표면간의 일치를 선언하며, 그렇지 않은 경우에는 상기 제 1 및 제 2 표면간의 불일치를 선언하는 단계를 포함하는 것을 특징으로 하는 방법.
KR1020057008156A 2002-11-06 2003-11-06 등각 구조에 의한 기하학적 표면의 분석 방법 KR100819960B1 (ko)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US42414102P 2002-11-06 2002-11-06
US60/424,141 2002-11-06

Publications (2)

Publication Number Publication Date
KR20050084991A true KR20050084991A (ko) 2005-08-29
KR100819960B1 KR100819960B1 (ko) 2008-04-07

Family

ID=32312758

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020057008156A KR100819960B1 (ko) 2002-11-06 2003-11-06 등각 구조에 의한 기하학적 표면의 분석 방법

Country Status (7)

Country Link
US (1) US20060013505A1 (ko)
EP (1) EP1559060A4 (ko)
JP (1) JP2006518886A (ko)
KR (1) KR100819960B1 (ko)
CN (1) CN1781111A (ko)
AU (1) AU2003286922A1 (ko)
WO (1) WO2004044689A2 (ko)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101401417B1 (ko) * 2012-10-22 2014-05-30 한국과학기술연구원 변형체 모델의 조각 시뮬레이션 방법 및 장치

Families Citing this family (78)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7421098B2 (en) 2002-10-07 2008-09-02 Technion Research & Development Foundation Ltd. Facial recognition and the open mouth problem
US6947579B2 (en) 2002-10-07 2005-09-20 Technion Research & Development Foundation Ltd. Three-dimensional face recognition
US7224356B2 (en) * 2004-06-08 2007-05-29 Microsoft Corporation Stretch-driven mesh parameterization using spectral analysis
US7127095B2 (en) 2004-10-15 2006-10-24 The Brigham And Women's Hospital, Inc. Factor analysis in medical imaging
US7671860B2 (en) * 2004-10-28 2010-03-02 Microsoft Corporation Representing implicit curves of procedural geometric surfaces
US7425959B2 (en) 2004-10-28 2008-09-16 Microsoft Corporation Representation of implicit curves for procedural surfaces
US7408548B2 (en) 2005-06-30 2008-08-05 Microsoft Corporation Triangulating procedural geometric objects
CN101268667B (zh) * 2005-09-23 2011-06-15 皇家飞利浦电子股份有限公司 估计发送数据向量的方法、码元确定设备和系统
US8786613B2 (en) 2006-04-08 2014-07-22 Alan Millman Method and system for interactive simulation of materials and models
EP2005342B1 (en) * 2006-04-08 2019-06-05 Allan Millman Method and system for interactive simulation of materials
WO2007132746A1 (ja) 2006-05-11 2007-11-22 Panasonic Corporation 処理装置
US8155207B2 (en) * 2008-01-09 2012-04-10 Cisco Technology, Inc. Processing and managing pictures at the concatenation of two video streams
US8243066B2 (en) * 2006-12-11 2012-08-14 Richard Garfinkle System and method for model creation and computer animation
US8538183B1 (en) 2007-03-08 2013-09-17 Nvidia Corporation System and method for approximating a diffusion profile utilizing gathered lighting information associated with an occluded portion of an object
US8064726B1 (en) * 2007-03-08 2011-11-22 Nvidia Corporation Apparatus and method for approximating a convolution function utilizing a sum of gaussian functions
US8773423B2 (en) * 2007-05-07 2014-07-08 Microsoft Corporation Creating optimized gradient mesh of a vector-based image from a raster-based image
US7868885B2 (en) * 2007-06-22 2011-01-11 Microsoft Corporation Direct manipulation of subdivision surfaces using a graphics processing unit
US8218844B2 (en) * 2008-02-21 2012-07-10 Siemens Aktiengesellschaft Method for generating analytical tool for analyzing abnormalities in human organs
US8368712B2 (en) * 2008-08-28 2013-02-05 Pixar Mesh transfer in n-D space
US20090213138A1 (en) * 2008-02-22 2009-08-27 Pixar Mesh transfer for shape blending
US8482569B2 (en) * 2008-08-28 2013-07-09 Pixar Mesh transfer using UV-space
EP2306891A1 (en) * 2008-07-08 2011-04-13 IT University of Copenhagen Eye gaze tracking
KR20100034843A (ko) * 2008-09-25 2010-04-02 (주)엠엑스알커뮤니케이션즈 3차원 얼굴 인식 보안 장치 및 그 방법
US8405659B2 (en) * 2009-06-24 2013-03-26 International Business Machines Corporation System and method for establishing correspondence, matching and repairing three dimensional surfaces of arbitrary genus and arbitrary topology in two dimensions using global parameterization
US8624893B2 (en) * 2009-07-31 2014-01-07 Adobe Systems Incorporated System and method for generating 2D texture coordinates for 3D meshed surfaces
AU2011207453A1 (en) * 2010-01-22 2012-08-02 The Research Foundation Of The State University Of New York System and method for prostate visualization and cancer detection
JP5914470B2 (ja) * 2010-06-21 2016-05-11 ゲニキャップ ビヘール ビー.ヴィー. コンピュータ実装ツールボックスシステム及び方法
US8861800B2 (en) * 2010-07-19 2014-10-14 Carnegie Mellon University Rapid 3D face reconstruction from a 2D image and methods using such rapid 3D face reconstruction
US20130170726A1 (en) * 2010-09-24 2013-07-04 The Research Foundation Of State University Of New York Registration of scanned objects obtained from different orientations
US8494237B2 (en) * 2010-11-08 2013-07-23 Cranial Technologies, Inc Method and apparatus for processing digital image representations of a head shape
US8442288B2 (en) * 2010-11-08 2013-05-14 Cranial Technologies, Inc. Method and apparatus for processing three-dimensional digital mesh image representative data of three-dimensional subjects
CN102542527B (zh) * 2011-12-27 2015-08-19 深圳万兴信息科技股份有限公司 一种变换图像形状的方法及系统
US9152743B2 (en) 2012-02-02 2015-10-06 Gehry Technologies, Inc. Computer process for determining best-fitting materials for constructing architectural surfaces
WO2013116859A1 (en) * 2012-02-02 2013-08-08 Gehry Technologies, Inc. Computer process for determining best-fitting materials for constructing architectural surfaces
US9298853B2 (en) * 2012-11-07 2016-03-29 Nec Laboratories America, Inc. System and methods for prioritizing queries under imprecise query execution time
EP2936288A1 (en) * 2012-12-18 2015-10-28 Van Der Westhuizen, Willem Morkel Event and causality-based human-computer interaction
CN104570928A (zh) * 2013-10-29 2015-04-29 中国科学院沈阳自动化研究所 基于共形参数化的网格曲面上数控加工轨迹规划方法
CN105719310B (zh) * 2014-08-13 2019-12-20 北京毕普创新科技有限公司 碰撞检测方法及装置
US10115035B2 (en) * 2015-01-08 2018-10-30 Sungkyunkwan University Foundation For Corporation Collaboration Vision system and analytical method for planar surface segmentation
JP2016142595A (ja) 2015-01-30 2016-08-08 富士通株式会社 移動体端末、位置特定方法、位置特定プログラムおよび位置特定装置
US9892506B2 (en) * 2015-05-28 2018-02-13 The Florida International University Board Of Trustees Systems and methods for shape analysis using landmark-driven quasiconformal mapping
US10262414B2 (en) * 2015-07-29 2019-04-16 University Of Louisville Research Foundation, Inc. Computer aided diagnostic system for mapping of brain images
CN105160358B (zh) * 2015-09-07 2018-08-17 苏州大学张家港工业技术研究院 一种图像分类方法及系统
US9940747B2 (en) * 2015-12-02 2018-04-10 Biosense Webster (Israel) Ltd. Mapping 3D to 2D images
EP3179389B1 (en) * 2015-12-07 2019-05-22 Dassault Systèmes Detecting cut-outs
US20170186208A1 (en) * 2015-12-28 2017-06-29 Shing-Tung Yau 3d surface morphing method based on conformal parameterization
KR101978316B1 (ko) * 2016-02-22 2019-05-14 연세대학교 산학협력단 메쉬모핑을 이용한 동맥의 혈류역학 3차원 볼륨 격자 생성 방법
WO2017175888A1 (ko) * 2016-04-05 2017-10-12 삼성전자 주식회사 이미지 처리 방법 및 장치
KR101812001B1 (ko) 2016-08-10 2017-12-27 주식회사 고영테크놀러지 3차원 데이터 정합 장치 및 방법
US9755661B1 (en) * 2016-08-16 2017-09-05 Macau University Of Science And Technology Method for coding data
CN108764140B (zh) * 2018-05-29 2021-11-30 中北大学 一种基于Ricci曲率流的三维人脸识别方法
CN108805063A (zh) * 2018-05-31 2018-11-13 王红军 一种对物体与环境的多分辨率视觉感知方法
US20210110602A1 (en) * 2018-06-25 2021-04-15 Hewlett-Packard Development Company, L.P. Three-dimensional registration procedures
CN108981717A (zh) * 2018-08-23 2018-12-11 王红军 一种基于双曲度规的路径规划方法
CN109190231B (zh) * 2018-08-24 2023-04-28 中建中原建筑设计院有限公司 基于grasshopper的折板网壳参数化建模方法
CN109377561B (zh) * 2018-10-26 2024-01-19 北京工商大学 一种基于共形几何的数模表面网格生成方法
US10622518B1 (en) * 2018-12-12 2020-04-14 Bolb Inc. Light-emitting diode with a mesa constructed from a unit cell
CN110176079B (zh) * 2019-05-26 2023-03-24 中北大学 一种基于准共形映射的三维模型变形算法
CN110188318B (zh) * 2019-05-28 2023-03-28 武汉轻工大学 曲面积分的自动计算方法、设备、存储介质及装置
CN110544310A (zh) * 2019-08-23 2019-12-06 太原师范学院 一种双曲共形映射下三维点云的特征分析方法
CN110488752B (zh) * 2019-08-23 2020-08-18 武汉数字化设计与制造创新中心有限公司 一种大型复杂曲面机器人自动化表面加工自主分片方法
CN111353187B (zh) * 2019-09-09 2022-12-09 中国铁路设计集团有限公司 一种获取空间线位平曲线要素的方法
CN110654565B (zh) * 2019-09-26 2021-10-08 东方航空技术有限公司 一种飞机外表复杂曲面蒙皮快速修理方法及系统
TWI812874B (zh) 2019-10-01 2023-08-21 美商杜拜研究特許公司 張量乘積之b平滑曲線預測子
CN110766808B (zh) 2019-11-05 2020-03-27 北京智拓视界科技有限责任公司 人脸数据的处理方法、设备和计算机可读存储介质
WO2021092329A1 (en) * 2019-11-06 2021-05-14 Nanotronics Imaging, Inc. Systems, methods, and media for manufacturing processes
CN115210744A (zh) * 2020-03-09 2022-10-18 纳米电子成像有限公司 用于制造过程的系统、方法和介质
CN111127492B (zh) * 2020-03-27 2020-07-17 北京智拓视界科技有限责任公司 用于脑图像注册的方法、装置和计算机可读存储介质
CN113240788A (zh) * 2020-07-08 2021-08-10 北京智拓视界科技有限责任公司 三维数据的传输和接收方法、设备和计算机可读存储介质
CN114117833B (zh) * 2020-08-28 2024-05-14 中国科学院沈阳自动化研究所 基于增材制造后处理的谢尔宾斯基曲线路径规划方法
CN112652047A (zh) * 2020-10-23 2021-04-13 成都完美时空网络技术有限公司 扭曲效果生成方法、装置、设备和存储介质
WO2022093216A1 (en) * 2020-10-28 2022-05-05 Hewlett-Packard Development Company, L.P. Generating conformal structures for 3d object models
CN113268789B (zh) * 2021-05-31 2024-01-12 西安交通大学 基于微分同胚的曲面注册方法、系统、设备及可读存储介质
CN113538219A (zh) * 2021-07-15 2021-10-22 北京智拓视界科技有限责任公司 三维人脸数据的传输和接收方法、设备和计算机可读存储介质
US11967025B2 (en) * 2021-11-05 2024-04-23 Tencent America LLC Detection of boundary loops in non-manifold meshes
CN115115729B (zh) * 2022-06-21 2024-04-30 浙江理工大学 一种基于球谐函数的纺织品图案单花型生成方法
CN116597109B (zh) * 2023-01-13 2023-12-22 东莘电磁科技(成都)有限公司 一种复杂三维曲面共型网格生成方法
CN117574470A (zh) * 2023-11-28 2024-02-20 清华大学 工业结构的智能化模型重构方法和装置

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6285372B1 (en) * 1998-05-08 2001-09-04 Lawrence C. Cowsar Multiresolution adaptive parameterization of surfaces
US6424745B1 (en) * 1998-05-19 2002-07-23 Lucent Technologies Inc. Method and apparatus for object recognition
US6256039B1 (en) * 1998-08-14 2001-07-03 The Board Of The Leland Stanford Junior University Methods for manipulating curves constrained to unparameterized surfaces
US6961055B2 (en) * 2001-05-09 2005-11-01 Free Radical Design Limited Methods and apparatus for constructing virtual environments

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101401417B1 (ko) * 2012-10-22 2014-05-30 한국과학기술연구원 변형체 모델의 조각 시뮬레이션 방법 및 장치

Also Published As

Publication number Publication date
JP2006518886A (ja) 2006-08-17
AU2003286922A8 (en) 2004-06-03
US20060013505A1 (en) 2006-01-19
WO2004044689A3 (en) 2004-09-02
WO2004044689A2 (en) 2004-05-27
CN1781111A (zh) 2006-05-31
KR100819960B1 (ko) 2008-04-07
EP1559060A2 (en) 2005-08-03
EP1559060A4 (en) 2007-06-13
AU2003286922A1 (en) 2004-06-03

Similar Documents

Publication Publication Date Title
KR100819960B1 (ko) 등각 구조에 의한 기하학적 표면의 분석 방법
Näf et al. 3D Voronoi skeletons and their usage for the characterization and recognition of 3D organ shape
Li et al. On computing mapping of 3d objects: A survey
Lim et al. Planar arrangement of high-dimensional biomedical data sets by isomap coordinates
Li et al. The wrinkle generation method for facial reconstruction based on extraction of partition wrinkle line features and fractal interpolation
Biasotti et al. Mathematical tools for shape analysis and description
Jin et al. Conformal geometry
Reina Molina et al. Cell complexes and membrane computing for thinning 2D and 3D images
Sivignon et al. Decomposition of a three-dimensional discrete object surface into discrete plane pieces
Hua et al. Spectral geometry of shapes
Mhamdi et al. A local approach for 3D object recognition through a set of size functions
Imielinska et al. Technical challenges of 3D visualization of large color data sets
Dai et al. IsoExplorer: an isosurface-driven framework for 3D shape analysis of biomedical volume data
Jia et al. 3D personalized human modeling and deformation technology for garment CAD
Krahnstoever et al. Computing curvature-adaptive surface triangulations of three-dimensional image data
Zhou et al. Adaptive algorithm for three-dimensional mesh generation based on constrained delaunay
Meyer Dynamic particles for adaptive sampling of implicit surfaces
Ge et al. Expansion of 3D face sample set based on genetic algorithm
CN113421326B (zh) 一种基于坐姿点检测的助行设备设计方法及系统
Anoshkina et al. Ridges and ravines on a surface and related geometry of skeletons, caustics, and wavefronts
Lin et al. Orthogonal integral transform for 3D shape recognition with few examples
Patane et al. Surface-and volume-based techniques for shape modeling and analysis
Ranjan A union of spheres representation for 3D objects
Jadhav et al. Volume Exploration Using Multidimensional Bhattacharyya Flow
Sun Deep learning on curved surfaces: manifold-formulation of convolutional neural networks and its operations

Legal Events

Date Code Title Description
A201 Request for examination
E902 Notification of reason for refusal
E90F Notification of reason for final refusal
E701 Decision to grant or registration of patent right
GRNT Written decision to grant
FPAY Annual fee payment

Payment date: 20110216

Year of fee payment: 4

LAPS Lapse due to unpaid annual fee