计算共形几何-USTC

计算共形几何暑期讲习班 - 大连理工大学数学科学学院

计算共形几何暑期讲习班时间：2012年7月2日-7月6日上午9：00地点：大连理工大学研教楼301报告人：顾险峰教授纽约州立大学石溪分校课程摘要：Computational conformal geometry is an emerging interdisciplinary field, which combines modern geometry and computer science, and is applied in a broad range of fields in engineering, medicine and fundamental sciences.This course will cover the following fundamental concepts and theories: homology, homotopy and covering space in algebraic topology; Riemann mapping, uniformization and conformal welding; Riemann surface, Hodge theory, Abel differential, Riemann-Roch theorem, conformal module, Tiechmuller space, quasi-conformal mapping theory; differential geometry of surfaces, exterior calculus, Gauss-Bonnet theorem, characteristic class; hyperbolic geometry, surface Ricci flow etc.The course will cover the following computational algorithms: representation of discrete surface, discrete exterior calculus, homology groups, homotopy groups, covering space, harmonic maps, harmonic differential forms, holomorphic differential forms, conformal structure, discrete conformal mapping, conformal module, discrete Ricci flow, uniformization metric, Teichmuller coordinates, quasi-conformal mapping, conformal welding, persistence homology and its stability, and so on.The course will introduce the following engineering applications: in computer graphics, mesh parameterization, texture mapping, vector field design; in computer vision, surface tracking and matching with large deformations, geometric database indexing, shape analysis; in geometric modeling, manifold spline theory, mesh-spline conversion; in medical imaging, conformal brain mapping, virtual colonoscopy; in computational geometry, homotopy detection, graph embedding; in wireless sensor network, routing based on metric design, load balancing; in computational topology, computing homology。

中科院课程

主要参考书： 1. J.柯歇尔，邹异明，《辛几何引论》，科学出版社，北京，1999。 2. Kobayashi, S. and Nomizu, K. ， Foundations of Differential
Geometry, Vol.1,2， Interscience Publishers， New York， 1969. 撰写人：肖良(中国科学院研究生院) 撰写日期： 2001 年 10 日
同构定理；共轭定理。第五章存在定理
通用包络代数；PBW 定理；生成元与定理关系。第六章表示理论
有限维表示；基础表示与初等表示；旋表示；表示的 Freudeuthal 公式；特征标理论；Weyl 公式；Kostant 公式和 Steinberg 公式。第七章李群与李代数
指数映射；伴随表示；李群与李代数。
本课程为数学学科各专业博士、硕士研究生的学科基础课。同时也可作为物理学、力学等专业研究生的选修课。微分流形己成为现代数学研究的基本对象。本课程讲授微分流形与李群的基本知识。通过本课程的学习，希望学生能初步掌握微分流形的基本概念、方法和技巧。为进一步学习微分几何、微分拓扑、几何分析等相关课程打下坚实基础。内容提要: 第一章欧氏空间
单纯同调群；奇异同调群；一般系数同调群；长正合同调列； Mayer-Vietoris 序列；球面同调群及几何应用；Lefschetz 不动点定理； CW 复形及其同调群。第四章上同调与对偶定理
上同调群；正合上同调列；上同调环；Poincare 对偶定理； Alexander 对偶定理；Lefschetz 对偶定理。
主要参考书： 1．Maunder, C.R.F.，Algebraic Topology，Cambridge University Press，

中科大考博辅导班：2019中科大数学科学学院考博难度解析及经验分享

中科大考博辅导班：2019中科大数学科学学院考博难度解析及经验分享中国科学院大学2019年博士研究生招生统一实行网上报名。

报考者须符合《中国科学院大学2019年招收攻读博士学位研究生简章》规定的报考条件。

考生在报考前请联系所报考的研究所(指招收博士生的中科院各研究院、所、中心、园、台、站)或校部相关院系，了解具体的报考规定。

下面是启道考博辅导班整理的关于中国科学技术大学数学科学学院考博相关内容。

一、院系简介数学科学学院的前身数学系于1958年由著名数学家华罗庚教授亲自主持创办并任首任系主任，关肇直、吴文俊、冯康、龚昇、王元、万哲先、陆启铿、石钟慈、林群、张景中、陈希孺等一大批知名专家曾在此任教。

2011年5月，数学科学学院正式挂牌成立，首任院长为马志明院士。

本院为首批全国理科人才培养基地、中国科学院博士生重点培养基地、长江学者特聘岗位设置学科，并获得首批数学一级学科博士学位授予权（涵盖数学所有博士点）,2007年获首批一级重点学科，是教育部985、211工程、中科院知识创新工程建设学科。

为吸引高水平的学者来我院讲学，学校为本院设立了“华罗庚大师讲席”及“吴文俊大师讲席”。

二、招生信息中国科学技术大学数学科学学院博士招生专业有1个：070100数学研究方向：随机分析与数理金融．计算机辅助几何设计．计算机图形学．应用逼近论、并行计算．组合优化．李代数及相关理论．微分几何．可积系统与子流形．几何分析．Ads/ds 空间的几何．可积系统．代数表示论．微分几何．非线性演化方程．可积系统．一维动力系统．材料科学与结构分析的计算方法研究．计算机辅助几何设计．计算机图形学．应用逼近．密码学．李代数及相关理．组合数学．信息安全．编码理论．无穷维系统控制．复杂系统控制及系统可靠．几何拓扑．拓扑量子场．动力系统．遍历理论．拓扑．图论．代数组合．偏微分方程．几何分析．亚纯函数值分布相关理论．一维动力系统．计算机辅助几何设计．计算机图形学．应用逼近论．生物数学．抛物方程动力学．应用分析．计算机图形学．图像处理．微分几何．离散几何分析．偏微分方程．几何分析．大范围分析．极值组合．图论．概率方法．组合优化．偏微分方程．算术代数几何．代数数论．函数空间与算子理论．Clifford分析．哈密尔顿系统．动力系统．遍历理论．拓扑学．一维动力系统．代数几何．微分方程动力系统．生物数学．几何分析．微分几何．代数群与量子群．算术代数几何．概率论与随机分析．偏微分方程数值方法．数据科学．最优化计算方．计算几何（三维打印的优化设计）．黎曼几何．动力系统．遍历理论．拓扑学．代数表示论．量子群与张量范畴．弦理论．共形场．最优化计算方法．随机优化．偏微分方程数值方法理论及应用．李代数与量子群．随机分析．微分几何．几何分析．组合设计与编码．宇宙学．偏微分方程数值解．可积系统三、报考条件（1）中华人民共和国公民；拥护中国共产党的领导，愿意为祖国社会主义现代化建设服务；品德良好，遵纪守法，学风端正，无任何考试作弊、学术剽窃及其它违法违纪行为；（2）身体健康状况符合我校规定的体检要求，心理正常；（3）申请者原则上应来自国内重点院校或所在高校学习专业为重点学科；（4）专业基础好、科研能力强，在某一领域或某些方面有特殊学术专长及突出学术成果；（5）对学术研究有浓厚的兴趣，有较强的创新意识、创新能力和专业能力；（6）申请者的学位必须符合下述条件之一：应届硕士毕业生须在博士入学前取得硕士学位；或已获得硕士或博士学位；在境外获得学位的考生，须凭教育部留学服务中心的认证书报名；（7）具有较强的语言能力，外语（限本单位招生专业目录中公布的语种）水平较高。

中国科学技术大学2017年线性代数与解析几何考研试题及解答

2.
设直线
l:
1−x 3
=y+1=
3−z 2
在平面
x−y+z
=2
上的投影为
l1,
则
l1
的方程为
,l
绕 l1 旋转所得的曲面方程是 .
101
3. 矩阵 1 1 = −1 1
3 + a1b1 a1b2
,
行列式
det
a2b1
3 + a2b2
a1b3
a2b3
=
.
a3b1
a3b2 3 + a3b3
3 −2 1
2. (15分) 考虑二阶复方阵 M (C) 组成的复线性空间, 方阵 A = 7 2 以及线性变换 B : 37
M2(C) → M2(C) 满足 B(X) = AX − XA, 其中 X 为任意 2 阶方阵, 试证明: B 是可对角化的线性变换.
3. (20分) 设 V 是由次数不超过 3 的实系数多项式组成的线性空间. 对于任意的 f (x), g(x) ∈
解得
a
=
3 5
.
2. 设 l 与平面的交点为 (1 − 3t, t − 1, 3 − 2t), 由交点在平面上得 1 − 3t − t + 1 + 3 − 2t = 2, 解得
t
=
1 2
,
于是交点为
−
1 2
,
−
1 2
,
2
. l1 的一个方向向量为 (−3, 1, −2) × (1, −1, 1) × (1, −1, 1) =
就马上得到结论. 至于上面例题的证明可以翻书查阅, 书上给了两种证明, 第二种证明与证明惯性定理类似.

计算共形几何讲义

计算共形几何讲义共形几何是数学中的一个重要分支，研究的是保持角度不变的变换。

它在几何学、物理学和工程学等领域都有广泛的应用。

本文将介绍共形几何的基本概念和计算方法。

共形几何的基本概念是保持角度不变的几何变换。

在平面几何中，共形变换包括平移、旋转和缩放三种基本变换。

平移是将图形沿着平行线移动，旋转是将图形绕着某一点旋转一定角度，缩放是将图形按照比例因子放大或缩小。

这三种变换都是保持角度不变的，因此它们都是共形变换。

在计算共形几何中，我们经常使用的工具是复变函数和复平面。

复变函数是将复数映射到复数的函数，它具有两个实部和虚部。

复平面是由复数构成的平面，其中实部和虚部分别对应于平面上的横轴和纵轴。

复变函数在共形几何中有着重要的应用，可以用来描述共形变换。

共形几何的计算方法主要包括两个方面：共形映射和共形度量。

共形映射是指将一个复变函数应用到一个图形上，从而得到一个经过共形变换的新图形。

共形映射可以通过复变函数的特定形式来实现，例如使用指数函数、对数函数或三角函数等。

共形度量是指计算两个图形之间的共形距离，即它们之间的角度差异。

共形度量可以通过计算两个图形的曲率来实现，曲率是指曲线在某一点处的弯曲程度。

共形几何在实际应用中有着广泛的应用。

在物理学中，共形几何被用来描述弯曲时空中的物理现象，例如黑洞的形成和宇宙膨胀等。

在工程学中，共形几何被用来设计无人机的飞行路径和机器人的运动轨迹等。

共形几何还被应用于图像处理和计算机图形学等领域，用来实现图像的变换和渲染。

总结起来，共形几何是研究保持角度不变的几何变换的数学分支。

它使用复变函数和复平面来描述共形变换，并通过共形映射和共形度量来计算共形距离。

共形几何在几何学、物理学和工程学等领域有着广泛的应用。

它不仅可以用来描述自然界中的现象，还可以用来解决工程和计算机图形学中的实际问题。

因此，深入了解和掌握共形几何的基本概念和计算方法对于学习和应用共形几何是非常重要的。

非线性问题初步-USTC

连续模型：微分方程（组）
线性算子和非线性算子
线性算子的定义
1） L(φ + ϕ ) = Lφ + Lϕ
2） L(cφ) = cLφ
例如：
L = d + ⎜⎛ d ⎟⎞2 dx ⎝ dx ⎠
离散模型：代数方程（组）
[K({a}]{a} = {f {a}}
(二) 固体力学平面问题的基本方程
１）几何关系（应变-位移方程）非线性
非线性问题初步
第一节基本概念
(一) 线性算子和非线性算子 (二) 固体力学平面问题的基本
方程 (三) 几种情况的示意图
(一) 线性算子和非线性算子
自然现象和工程问题的数学模型数学模型的分类离散模型和连续模型
静态问题动态问题特征值问题线性问题和非线性问题
L(φcφ+)ϕ=)c=LφLφ + Lϕ
根据几何关系，有: sin β = Lsinα −u L−∆
L cosα = ( L − ∆) cos β L sin α = ( L − ∆) sin β + u
由上式得到
Lα = L − ∆ = L2 − 2Lu sin α + u 2 sin β = L sinα − u
L−∆
平衡方程是（变形后形态的平衡方程）
了，而又抓住了问题本质。
回到线性问题的有限元法，我们得到的是一组线性代数方程
组 [K]{a} = { f }
{ f } = ∫De [ N ]T {b}dD ：外载荷向量，对变形前区域积
分，没有区分变形前后形状的变化，变形以后的区域是未知
的。{a}
：节点位移向量
{ } [K] = ∫De [B]T[D][B]dD ：刚度矩阵，与 a 无关，

计算共形几何

Note: This is mostly copied or translated from the papers ofProf. Xianfeng Gu (/~gu/).计算共形几何张威应用数学110060591引言共形几何是纯数学中很多学科的交叉领域,比如黎曼曲面理论、微分几何、代数曲线、代数拓扑、偏微分方程、复分析等等.它有很长的历史,至今在现代几何与现代物理中仍然非常活跃.比如超弦理论中的共形场和理论物理中的模空间理论都是当今快速发展的研究领域.近些年来,随着三维数字扫描仪、计算机辅助几何设计、生物信息和医学成像的快速发展,出现了越来越多的三维数字模型.因此迫切需要有效的算法来表示、处理和使用这些模型.计算共形几何在数字几何处理中扮演了一个重要角色.它已经应用在很多重要的领域,比如曲面修复、光顺、去噪、分片、特征提取、注册、重新网格化、网格样条转换、动画和纹理合成.特别地,共形几何奠定了曲面参数化的理论基础,同时也提供了严格的算法.计算共形几何还应用于计算机视觉中的人脸跟踪、识别和表情转换,医学成像中的脑电图、虚拟结肠镜和数据融合,几何建模中的具有任意拓扑流形上的样条构造.共形几何之所以如此有用是基于以下一些事实:•共形几何研究的是共形结构.日常生活中的所有曲面都有一个自然的共形结构,因此共形几何算法非常普遍.•共形结构比黎曼度量结构更灵活、比拓扑结构更具有刚性.它能处理大量黎曼几何不能有效处理的变换,这些变换还能保持很多拓扑方法会丢失的几何信息.•共形映射比较容易控制.比如,两个单连通封闭曲面之间的共形映射构成一个6维空间,因此只要固定3个点,这个映射就是唯一的.这个事实使得共形几何方法在曲面匹配和比较中非常有价值.•共形映射保持局部形状,因此在可视化方面有很好的应用.•所有的曲面都可以根据共形结构进行分类,而且所有的共形等价类形成一个有限维流形.这个流形有丰富的几何结构,容易对其分析和研究.与之相反,曲面的等距类形成一个难以分析处理的无穷维流形.计算共形几何•计算共形几何算法是以椭圆偏微分方程为基础的,而椭圆偏微分方程又是容易求解而且稳定的,因此计算共形几何方法对于实际工程应用非常有用.•共形几何中,所有单连通曲面都能共形变换成某种标准空间:球面、平面、双曲空间.也就是说,任何曲面都具有三种标准几何(球几何、欧式几何、双曲几何)中的一种.这样大部分三维数字几何处理任务都能转化成二维标准空间中的任务.历史上,计算共形几何方法已经广泛应用于许多工程领域,然而绝大部分应用都基于平面区域的共形映射.最近,随着数学理论的发展和计算能力的提升,计算共形几何方法已经从平面区域推广到具有任意拓扑的曲面.1.1共形变换和共形结构图1:共形映射保持角度(a)Circle packing(b)Checkboard图2:共形映射根据Fleix Klein的Erlangen纲领,几何就是研究在特定的变换群下保持不变的空间性质.共形几何就是研究保角变换群下的不变量.它介乎于拓扑和黎曼几何之间.计算共形几何共形映射就是保角映射,如图1所示.在无穷小邻域,共形映射就是放缩变换.它保持局部形状,比如它将无穷小圆周映成无穷小圆周.如图2所示,这个bunny曲面通过一个共形映射映到平面.如果平面有一个circle packing,则通过拉回得到bunny曲面上的一个circle packing.如果给平面铺上棋盘格,则同样得到bunny曲面的棋盘格修饰,其中直角和正方形都是保持的.曲面上的两个黎曼度量是共形的,如果它们定义的角度是相同的.共形结构就是指曲面上度量的共形等价类,而黎曼曲面就是带有共形结构的光滑曲面.因此在黎曼曲面上,我们可以度量角度,但不能度量长度.每一个带有度量的曲面都自动成为一个黎曼曲面.如果两个黎曼曲面之间存在共形映射,则称它们是共形等价的.显然,共形等价是黎曼曲面间的一个自然的等价关系.共形几何的目的就是在共形等价意义下对黎曼曲面进行分类,这就是所谓的模空间问题.给定一张光滑曲面,考察它上面的所有共形结构在共形等价下的模,这个集合被称为曲面的模空间.对于具有正亏格的封闭曲面,模空间是正维数的有限维空间.1.2基本任务下面的问题是计算共形几何最基本的一些任务.这些问题是相互依赖的:1.共形结构给定一张带有黎曼度量的曲面,计算它的内蕴共形结构的不同表示.一种方法是计算它的Abelian微分群,另一种方法是计算标准的黎曼度量.2.共形模完全共形不变量称为黎曼曲面的共形模.正如前面所讲的,理论上存在一组有限的数完全决定了黎曼曲面,这些称为黎曼曲面的共形模.一个比较难的问题是显式计算任意给定曲面的共形模.3.标准黎曼度量黎曼曲面的uniformization定理揭示了每一个黎曼度量都共形等价于一个常Gauss 曲率度量.除了球面和环面外,这个度量是唯一的.计算这个度量在计算共形几何中具有基本的重要性.4.共形映射计算共形几何计算两个共形等价的曲面之间的共形映射可以简化为计算它们到标准形状空间(球面、平面、双曲空间中的圆域)之间的共形映射.5.拟共形映射大部分微分同胚都不是共形的,它们将无穷小圆周映成无穷小椭圆.如果这些椭圆的长短轴比一致有界,那么就称为拟共形映射.拟共形映射的微分是由所谓的Beltrami微分刻画的,Beltrami微分记录了长轴方向和长短轴比.有一个基本定理是说通过Beltrami微分可以恢复拟共形映射.而至于怎样通过Beltrami微分计算拟共形映射,则是一个具有很多应用价值的重要问题.6.共形镶嵌粘合带边黎曼曲面并研究缝曲线形状和粘合样式之间的关系.这与拟共形映射问题紧密相关.1.3共形几何在工程应用中的优点计算共形几何已经被证实在许多工程领域中有重要应用.下面是一些主要理由:1.标准区域所有带度量的曲面都能共形地映成球面、平面或双曲圆盘中的标准区域.这可以帮助我们将三维几何处理问题转化为二维问题.2.通过曲率设计度量每一个共形结构都有一个常Gauss曲率的标准度量,这个度量在很多几何应用中非常有价值.例如在双曲度量下,每一个非平凡同伦类都有一个闭测地线代表元.此外,我们可以根据预先给定的曲率设计黎曼度量,这在几何建模中非常有用.3.一般几何结构共形几何方法能够用来构造其他几何结构,比如仿射结构、射影结构等等.这些结构在几何建模应用中是关键的.4.微分同胚的构造共形映照和拟共形映照可以被用来构造曲面间的微分同胚.可以应用于曲面注册和比较这些计算机视觉和医学图像中最基本的问题.计算共形几何5.等温坐标共形结构可以被当作曲面上的等温坐标图册.在这种坐标下,度量的表达式最简单,因此所有的微分算子,比如Laplace-Beltrami算子,具有很简洁的表达式.这可以用来简化偏微分方程.等温坐标保持局部形状,对于可视化和纹理映射是非常完美的.2已有工作图3:基本群的多边形表示图4:四种类型的图.(a)A cut graph但不是system of loops.(b)A system of loops.(c)基本群的基但不是system of loops.(d)同调基,但既不是同伦基也不是cut graph.计算共形几何已经有比较长的历史了,起初大部分共形几何方法都是针对平面区域或拓扑圆盘的.直到近些年来才推广到具有复杂拓扑的曲面上,这与当前研究热点计算拓扑的发展是紧密相关的.有相当一部分计算共形算法是依赖于计算拓扑算法的.计算共形几何图5:从左到右:输入模型,handle(绿)和tunnel(红)环,handle特征2.1计算拓扑计算拓扑的研究对象主要是单纯复形,比如三角网格曲面就是一种二维单纯复形.计算拓扑的一些基本问题就是计算曲面单纯复形的cut graph,同调群和基本群.知道了cut graph,就能将曲面剪开成拓扑圆盘,这对于参数化和纹理映照是非常有用的.知道了基本群,就能完全确定封闭曲面的拓扑分类.通常基本群是用多边形表示的,其中有约化多边形表示、标准多边形表示等一些特殊的多边形表示,见图3.Vegter和Yap[1]给出了计算标准多边形表示的算法,Lazarus等[2]简化了Vegter和Yap的算法,复杂度为O(gn),其中n,g分别为曲面的大小和亏格.标准多边形表示由2g 个具有公共点的环组成,令人不满意的是它们可能具有公共边.在大部分应用中,标准多边形表示不是必须的,只需要知道cut graph就足够了.cut graph是指三角网格曲面上的一族边集使得曲面去掉这些边后变成拓扑圆盘,如图4所示.Dey和Schipper[3]使用cut graph得到了检验曲线是否可缩或两曲线是否同伦的快速算法.2002年,Erickson等[4]第一次提出了最优cut graph的问题:给每一条边赋一个权,求权最小的cut graph.这篇文章证明了最小cut graph问题是NP难的,并给出了计算近似最小cut graph的贪婪算法,该算法复杂度为O(g2n log n),精度为O(log2g).de Verdière和Lazarus[5]研究了一种称为"systerm of loops"的特殊cut graph,每一个systerm of loops都由具有公共基点的2g个环组成,如图4所示.Erickson等[6]给出了一个简单的贪婪算法来计算基本群和第一同调群的生成元,这同时也解决了de Verdière和Lazarus提出了最短systerm of loops问题.计算共形几何曲面上一些与特征紧密相关的环,比如handles和tunnels,在诸如拓扑修复、曲面参数化和特征识别等应用中非常有用,如图5所示.Dey等[7--9]给出了计算handles和tunnels的算法.Yin等[10]提出了使用universal covering space来计算最短非平凡环的算法.de Verdière[11]使用最优pants分解来找同伦类中的最短环.Xin等[12]研究了高亏格曲面上的测地环.关于曲面上的非平凡环,还有许多文章[13--19],目前仍是一个热点问题.(a)调和1形式(b)共轭的调和1形式图6:调和1形式2.2平面区域传统的数值复分析方法集中在平面区域的共形映射.[20--25]作了详尽总结. Schwarz-Christoffel映射已经被广泛应用于计算共形映射[26,27].[28--30]研究了多连通区域的Schwarz-Christoffel映射.[31]给出了一种基于交比和Delaunay三角化的鲁棒算法.最近,Marshall等[32]引进了基于迭代简单映射的geodesic zippper算法,Bishop[33]提出了基于双曲几何的线性复杂度的共形映射算法.2.3亏格为0的曲面[34]利用余切公式构造了离散调和映射,Lévy等[35]引进了Cauchy-Riemann方程的一阶有限元逼近.2002年,Desbrun等[36]最小化Dirichlet能量得到了离散内蕴参数化. Floater[37]引进平均值坐标来计算广义调和映射.[38,39]对拓扑球面的共形映射进行了研究.在CG领域,[40,41]对拓扑圆盘和球面上的共形映射进行了详尽地总结.计算共形几何2.4高亏格曲面高亏格曲面上共形结构的计算有两个主要方法:全纯微分方法和离散曲率流方法.图7:全纯1形式2.4.1全纯微分离散全纯微分方法是由Gu和Yau[42,43]提出的,用来计算高亏格曲面上的共形结构.这个方法基于Hodge理论,使用热扩散方法来计算每一个上同调类中的调和形式,见图6(a)和6(b).然后利用Hodge星算子构造全纯形式,见图7.所有计算都是在离散多边形曲面上进行的.Pinkall等[34]定义了一种不同的离散Hodge星算子来计算极小曲面. Mercat[44]通过离散Cauchy-Riemann方程来构造离散全纯映射,不过这个方法要求曲面是四边形网格.2003年,Hirani在他的博士论文[45]中详细介绍了离散外微分方法.Gortler等[46,47]使用离散1形式来对亏格为1的曲面进行参数化.Tong等[48]推广到带有锥奇点的1形式方法,并用来进行remeshing和tiling.2009年,Zeng等[49,50]将全纯微分方法应用于计算带多个边界的亏格为0的曲面上的共形映射,以及拟共形映射.2.4.2Ricci FlowRicci流是由Hamilton[51]在Princeton的讨论班上引进的.Ricci流对曲面和三维流形的几何研究产生了革命性的影响,是当前几何研究中极为活跃的方向.特别地,计算共形几何它导致了三维Poincaré猜测的证明.Hamilton[52]使用2维Ricci流给出了正亏格曲面uniformization定理的证明.这表明了它在CG领域中有极大的潜在应用价值.图8:Circle packing存在很多方法离散光滑曲面,其中与共形映射离散特别相关的是由Thurston[53]引进的circle packing度量.circle packing这个概念首先是由Koebe提出的,Thurston猜测平面Jordan区域的circle packing离散序列收敛于Riemann映射,见图8.这个猜测被Rodin 和Sullivan[54]所证明.Colin de Verdière[55]建立了circle packing的第一变分原理,然后证明了Thurston的circle packing度量存在性.这为实现计算circle packing度量的快速算法铺平了道路,比如Collins和Stephenson[56]就给出了一个算法.Chow和Luo[57]推广了Colin de Verdière 的工作,并引进了曲面上的离散Ricci流和离散Ricci能量.它们证明了离散Ricci流的一般存在性和收敛性定理,同时证明了Ricci能量的凸的.基于此,Jin[58]在2008年给出了离散Ricci流的算法实现.另外一个相关的离散方法是circle pattern,它同时考虑网格曲面的组合与几何信息,可以看成是circle packing的变种.Circle pattern是由Bowers和Hurdal[59]提出的,已经被证明是Bobenko和Springborn[60]提出的凸能量的极小元.2006年,Kharevych等[61]给出了一个有效的circle pattern算法.2.4.3Yamabe FlowYamabe问题目的是为紧Riemann流形找到具有常标量曲率的共形度量.Yamabe[62]给出了第一个带有漏洞的证明,之后被一些包括Trudinger[63],Aubin[64]和Schoen[65]在内的研究者所更正.关于这个问题的详细总结可以参看Lee和Parker的文章[66].2004年,Luo[67]研究了曲面上的离散Yamabe流.他引进了多边形度量的离散共形计算共形几何变换概念,这在研究离散Yamabe流和相应的变分原理方面具有关键性的作用.基于离散共形类和几何方面的考虑,Luo定义离散Yamabe能量为微分1形式的积分,然后证明了这个能量是一个局部凸函数.他还从这得出Yamabe流下的曲率发展是一个热方程.最近,Springborn等[68]的一个非常漂亮的工作表明可以将Yamabe能量和Milnor-Lobachevsky函数等同起来,将曲率发展的热方程和余切Laplace方程等同起来.近期Gu 等[69]使用离散Yamabe能量构造了一个计算离散共形度量的有效算法.Bobenko等[70]讨论了离散双曲Yamabe流,它被Zeng等[71]用来计算双曲结构和标准同伦群生成元.图9:全局共形的纹理映照图10:使用特殊的平坦度量设计向量场2.5应用2.5.1图形学共形几何方法广泛应用于CG中.等温坐标对于全局共形参数化来说是很自然的[42].因为共形映照不会扭曲局部形状,所以它对于纹理映照来说是非常理想的.图9显示了一个使用全纯微分进行纹理映照的例子.特殊的平坦度量对于设计曲面上的向量场是非常有价值的,图10显示了使用曲率流方法设计的向量场.计算共形几何图11:流形样条框架2.5.2几何建模几何建模中的一个最基本的问题就是将传统的欧式区域上的样条系统地推广到流形区域上,这与曲面上的一般几何结构关系紧密,见图11和12.传统的样条方案是建立在仿射不变量上的,如果流形具有仿射结构,那么传统的样条也能直接推广到流形上.然而由于拓扑障碍,一般的流形不具有仿射结构.但是删去一些奇点后,曲面上允许存在仿射结构.详细的讨论请参考Gu等[72].仿射结构可以通过共形几何方法显式计算.比如可以通过曲率流和全纯微分方法来计算仿射结构,图13显示的是用曲率流方法构造的仿射结构.与其他的方法比较起来,比如基于平凡联络的方法[73],这些方法的优点是它能给出样条曲面的全局共形参数化,即等温坐标.在这种坐标下,微分算子,比如梯度算子和Laplace-Beltrami算子,具有最简单的形式.2.5.3医学成像共形几何在医学成像的很多领域都有应用.比如在脑电图中,不同部分的大脑皮层曲面片的注册是关键性的.大脑表面高度卷曲,而且不同的人的大脑有不同的解剖结构,计算共形几何图12:流形样条图13:仿射结构因此在大脑皮层曲面片之间找到好的匹配是非常有挑战性的任务.图14演示了一个解决方案[39],先通过用标准方式将大脑表面映到球面,然后找到球面的自同构,这样曲面片的注册就很容易建立了.2.5.4视觉曲面匹配在计算机视觉中是一个基本问题.图15演示了曲面匹配的基本框架.图16演示了同一个人的不同表情的人脸匹配的例子,细节可参考[74--76].Teichmüller理论可以应用于曲面分类[77,78].通过使用Ricci曲率流,我们可以计算双曲uniformization度量.然后我们使用测地线进行pants分解并计算Fenchel-Nielsen坐标.图17演示了万有覆盖空间在双曲空间中的一个有限部分的计算,图18演示了计算Teichmüller坐标的流程.计算共形几何图14:大脑表面的共形映照图15:曲面匹配框架2.5.5计算几何在计算几何中,同伦检测是一个重要的问题.可以使用Ricci流计算双曲uniformiza-tion度量[71].根据Gauss-Bonnet定理,每一个同伦类中有唯一的闭测地线.给定一个环,我们可以计算这个环的同伦类对应的Möbius变换,这个变换的轴就是双曲度量下的闭测地线.如图19所示,同伦的环的标准代表元是相同的.计算共形几何图16:不同表情的人脸匹配图17:万有覆盖空间在双曲空间中的有限部分3总结计算共形几何是一门数学和计算机科学的交叉学科,它对数学理论的掌握要求比较高.很多算法都依赖于拓扑算法,与计算拓扑紧密相关.本文介绍了这个方向的一些基本问题、算法以及应用.目前仍然有很多基本问题没有解决,它们需要更深刻的理解和更强大精确的计算方法.下面列举一些对理论和应用都非常重要的问题.1.Teichmüller映射给定两个度量曲面和映射同伦类,计算唯一的angle distortion最小的映射,也就是Teichmüller映射.计算共形几何图18:Teichmüller空间中的Fenchel-Nielsen坐标图19:使用双曲度量进行同伦检测2.Abel微分计算不同类型的Abel微分群,特别是全纯的二次微分.3.组合结构和共形结构之间的关系给定一个拓扑曲面,每一个三角化都有一个由circle packing决定的自然的共形结构,讨论这两个结构之间的关系.4.逼近理论尽管计算共形不变量的算法已经有了,但是逼近的理论结果仍没有建立.对于平面区域的共形映射,不同离散方法的收敛性都已经建立了.但是对于一般曲面而言,收敛性分析仍然是开放的.5.精度和稳定性双曲几何计算对于数值误差是非常敏感的.要提高计算精度是非常具有挑战性的.计算几何中的精确计算方法给出了解决这个问题的方向.在逆距离circle packing计算共形几何方法和组合Yamabe流方法中,可接受曲率空间的非凸性造成了算法的不稳定.因此,曲面需要高精度的三角化.在实际应用中,提高三角化精度对这些算法是非常重要的.交角为锐角的circle packing算法更加稳定,全纯微分方法最稳定.计算共形几何参考文献[1]Vegter,G.,Yap,C.K.,putational complexity of combinatorial surfaces.In:Proceedingsof the sixth annual symposium on Computational geometry.New York:ACM,102-111.[2]Lazarus,F.,Pocchiola,M.,Vegter,G.,et al,puting a canonical polygonal schema of anorientable triangulated surface.In:SCG'01:Proceedings of the seventeenth annual symposium on Computational geometry.New York:ACM,80-89.[3]Dey,T.K.,1994.A new technique to compute polygonal schema for2-manifolds with applica-tion to null-homotopy detection.In:Proceedings of the tenth annual symposium on Computational geometry.New York:ACM,277-284.[4]Erickson,J.,Har-Peled,S.,2002.Optimally cutting a surface into a disk.In:SCG'02:Proceedingsof the eighteenth annual symposium on Computational geometry.New York:ACM,244-253.[5]de Verdière,E.C.,Lazarus,F.,2002.Optimal system of loops on an orientable surface.In:Proc.43rd Annual IEEE Symp.Foundations of Computer Science.Springer-Verlag627-636.[6]Erickson,J.,Whittlesey,K.,2005.Greedy optimal homotopy and homology generators.In:SODA'05:Proceedings of the sixteenth annual ACM-SIAM symposium on Discrete algorithms.Philadel-phia,PA,USA:Society for Industrial and Applied Mathematics,1038-1046.[7]Dey,T.K.,Li,K.,Sun,J.,2007.On computing handle and tunnel loops.In:CW'07:Proceedingsof the2007International Conference on Cyberworlds.Washington,DC,USA:IEEE Computer Society,357-366.[8]Dey,T.K.,Li,K.,Sun,J.,et al,puting geometry-aware handle and tunnel loops in3dmodels.In:SIGGRAPH'08:ACM SIGGRAPH2008papers.New York:ACM,1-9.[9]Dey,T.K.,Li,K.,Sun,J.,puting handle and tunnel loops with knot put.Aided Des.,41(10):730-738.[10]Yin,X.,Jin,M.,Gu,X.,puting shortest cycles using universal covering space.Vis.Comput.,23(12):999-1004.[11]Verdière,E.C.D.,Lazarus,F.,2007.Optimal pants decompositions and shortest homotopic cycleson an orientable surface.J.ACM,54(4):18.[12]Xin,S.Q.,He,Y.,Fu,C.W.,et al,2011.Euclidean geodesic loops on high-genus surfaces appliedto the morphometry of vestibular systems.In:MICCAI'11:Proceedings of the14th international conference on Medical image computing and computer-assisted intervention.Berlin,Heidelberg: Springer-Verlag,384-392.[13]Chambers,E.W.,Erickson,J.,Worah,P.,2008.Testing contractibility in planar rips complexes.In:SCG'08:Proceedings of the twenty-fourth annual symposium on Computational geometry.New York:ACM,251-259.[14]Chambers,E.W.,Erickson,J.,Nayyeri,A.,2009.Minimum cuts and shortest homologous cycles.In:SCG'09:Proceedings of the25th annual symposium on Computational geometry.New York: ACM,377-385.[15]Chambers,E.W.,Erickson,J.,Nayyeri,A.,2009.Homology flows,cohomology cuts.In:STOC计算共形几何'09:Proceedings of the41st annual ACM symposium on Theory of computing.New York:ACM, 273-282.[16]Erickson,J.,2010.Maximum flows and parametric shortest paths in planar graphs.In:SODA'10:Proceedings of the Twenty-First Annual ACM-SIAM Symposium on Discrete Algorithms.Philadelphia,PA,USA:Society for Industrial and Applied Mathematics,794-804.[17]Erickson,J.,Worah,P.,puting the shortest essential cycle.Discrete Comput.Geom.,44(4):912-930.[18]Dey,T.K.,Hirani,A.N.,Krishnamoorthy,B.,2010.Optimal homologous cycles,total unimodular-ity,and linear programming.In:STOC'10:Proceedings of the42nd ACM symposium on Theory of computing.New York:ACM,221-230.[19]Erickson,J.,2011.Shortest non-trivial cycles in directed surface graphs.In:SoCG'11:Proceedingsof the27th annual ACM symposium on Computational geometry.New York:ACM,236-243.[20]Delillo,T.K.,1994.The accuracy of numerical conformal mapping methods:A survey of examplesand results.SIAM Journal on Numerical Analysis,31(3):pp.788-812.[21]Ivanov,M.,V.T,1995.Handbook of conformal mapping with computer-aided visualization.BocaRaton:CRC Press.[22]Binder,I.,Braverman,M.,Yampolsky,M.,2007.On the computational complexity of the riemannmapping.ARKIV FOR MATEMATIK,45(2):221-239.[23]Trefethen,L.,1986.Numerical conformal mapping.Amsterdam:North-Holland Publishing Co.[24]Wegmann,R.,2005.Handbook of Complex Analysis:Geometric Function Theory.vol2.Ams-terdam:Elsevier,351-477.[25]Henrici,,1993.Discrete Fourier Analysis,Cauchy Integrals,Construction of Conformal Maps,Univalent Functions.vol3.London:Wiley.[26]Driscoll,T.A.,Trefethen,L.N.2002.Schwarz-Christoffel Mapping.Cambridge Press.[27]Banjai,L.,Trefethen,L.N.,2003.A multipole method for schwarz--christoffel mapping of poly-gons with thousands of sides.SIAM put.,25(3):1042-1065.[28]DeLillo,T.K.,Elcrat,A.R.,Pfaltzgraff,J.A.,2001.Schwarz-christoffel mapping of the annulus.SIAM Review,43(3):pp.469-477.[29]Delillo,T.,Elcrat,A.,Pfaltzgraff,J.,2004.Schwarz-christoffel mapping of multiply connecteddomains.Journal d'Analyse Mathématique,94:17-47.[30]Crowdy,D.,2005.The schwarz-christoffel mapping to bounded multiply connected polygonaldomains.Proceedings:Mathematical,Physical and Engineering Sciences,461(2061):pp.2653-2678.[31]Driscoll,T.A.,Vavasis,S.A.,1998.Numerical conformal mapping using cross-ratios and delaunaytriangulation.SIAM put.,19:1783-1803.[32]Marshall,D.E.,Rohde,S.,2007.Convergence of a variant of the zipper algorithm for conformalmapping.SIAM J.Numer.Anal.,45:2577-2609.计算共形几何[33]Bishop,C.J.,2010.Conformal mapping in linear time.Discrete Comput.Geom.,44:330-428.[34]Pinkall,U.,Polthier,K.,puting discrete minimal surfaces and their conjugates.Exper-imental Mathematics,2:15-36.[35]Lévy,B.,Petitjean,S.,Ray,N.,et al,2002.Least squares conformal maps for automatic texture atlasgeneration.In:Proceedings of the29th annual conference on Computer graphics and interactive techniques.New York:ACM,362-371.[36]Desbrun,M.,Meyer,M.,Alliez,P.,2002.Intrinsic parameterizations of surface puterGraphics Forum,21(3):209-218.[37]Floater,M.S.,2003.Mean value put.Aided Geom.Des.,20:19-27.[38]Gotsman,C.,Gu,X.,Sheffer,A.,2003.Fundamentals of spherical parameterization for3d meshes.In:SIGGRAPH'03:ACM SIGGRAPH2003Papers.New York:ACM,358-363.[39]Gu,X.,Wang,Y.,Chan,T.F.,et al,2004.Genus zero surface conformal mapping and its applicationto brain surface mapping.IEEE Transactions on Medical Imaging,23:949-958.[40]Floater,M.S.,Hormann,K.,2005.Surface parameterization:a tutorial and survey.In:Dodgson,N.A.,Floater,M.S.,Sabin,M.A..Advances in multiresolution for geometric modelling.Springer Verlag,157-186.[41]Sheffer,A.,Praun,E.,Rose,K.,2006.Mesh parameterization methods and their applications.put.Graph.Vis.,2:105-171.[42]Gu,X.,Yau,S.T.,2003.Global conformal surface parameterization.In:SGP'03:Proceedingsof the2003Eurographics/ACM SIGGRAPH symposium on Geometry processing.Aire-la-Ville, Switzerland,Switzerland:Eurographics Association,127-137.[43]Jin,M.,Wang,Y.,Yau,S.T.,et al,2004.Optimal global conformal surface parameterization.In:VIS'04:Proceedings of the conference on Visualization'04.Washington,DC,USA:IEEE Computer Society,267-274.[44]Mercat,C.,2001.Discrete riemann surfaces and the ising munications in MathematicalPhysics,218:177-216.[45]Hirani,A.N.,2003.Discrete exterior calculus.PhD thesis.Pasadena,CA,USA:California Instituteof Technology.[46]Gortler,S.J.,Gotsman,C.,Thurston,D.,2006.Discrete one-forms on meshes and applications to3d mesh put.Aided Geom.Des.,23:83-112.[47]Tewari,G.,Gotsman,C.,Gortler,S.J.,2006.Meshing genus-1point clouds using discrete one-puters and Graphics,30:917-926.[48]Tong,Y.,Alliez,P.,Cohen-Steiner,D.,et al,2006.Designing quadrangulations with discreteharmonic forms.In:Proceedings of the fourth Eurographics symposium on Geometry processing.Aire-la-Ville,Switzerland,Switzerland:Eurographics Association,201-210.[49]Zeng,W.,Yin,X.,Zhang,M.,et al,2009.Generalized koebe's method for conformal mappingmultiply connected domains.In:SPM'09:2009SIAM/ACM Joint Conference on Geometric and计算共形几何Physical Modeling.New York:ACM,89-100.[50]Zeng,W.,Luo,F.,Yau,S.T.,et al,2009.Surface quasi-conformal mapping by solving beltramiequations.In:Proceedings of the13th IMA International Conference on Mathematics of Surfaces XIII.Berlin,Heidelberg:Springer-Verlag,391-408.[51]Hamilton,R.,1982.Three manifolds with positive ricci curvature.J.Differ.Geom.,17:255-306.[52]Hamilton,R.,1988.Mathematics and General Relativity.vol71.American Mathematical Society.[53]Thurston,W.,1980.Geometry and Topology of Three-Manifolds.Princeton University Press.[54]Rodin,B.,Sullivan,D.,1987.The convergence of circle packings to the riemann mapping.J.Differ.Geom.,26(2):349-360.[55]Colin de Verdière,Y.,1991.Un principle variationnel pour les empilements de cercles.Invent.Math.,104(3):655-669.[56]Collins,C.R.,Stephenson,K.,2003.A circle packing put.Geom.Theory Appl.,25:233-256.[57]Chow,B.,Luo,F.,binatorial ricci flows on surfaces.J.Differ.Geom.,63(1):97-129.[58]Jin,M.,Kim,J.,Luo,F.,et al,2008.Discrete surface ricci flow.IEEE Transactions on Visualizationand Computer Graphics,14(5):1030-1043.[59]Bowers,P.L.,Hurdal,M.K.,2003.Planar conformal mappings of piecewise flat surfaces.In:In Vi-sualization and Mathematics III,H.-C.Hege and K.Polthier,Eds.Mathematics and Visualization.Springer-Verlag3-34.[60]Bobenko,A.I.,Springborn,B.A.,2004.Variantional principles for circle patterns and koebe'stheorem.Trans.Am.Math.Soc.,356:659-689.[61]Kharevych,L.,Springborn,B.,Schröder,P.,2006.Discrete conformal mappings via circle patterns.ACM Trans.Graph.,25:412-438.[62]Yamabe,H.,1960.The yamabe problem.Osaka Math.J.,12(1):21-37.[63]Trudinger,N.,1968.Remarks concerning the conformal deformation of riemannian structures oncompact manifolds.Ann.Scuola Norm.Sup.Pisa,22(2):265-274.[64]Aubin,T.,1976.équations diffréntielles non linéaires et problème de yamabe concernant la cour-bure scalaire.J.Math.Pure Appl.,55(3):269-296.[65]Schoen,R.,1984.Conformal deformation of a riemann metric to constant scalar curvature.J.Differ.Geom.,20(2):479-495.[66]Lee,J.M.,Parker,T.H.,1987.The yamabe problem.Bull.Am.Math.Soc.,17(1):39-71.[67]Luo,F.,Khan,L.,Bastani,F.,et al,2004.A dynamically growing self-organizing tree(dgsot)forhierarchical clustering gene expression profiles.Bioinformatics,20(16):2605-2617.[68]Springborn,B.,Schröder,P.,Pinkall,U.,2008.Conformal equivalence of triangle meshes.In:SIGGRAPH Asia'08:ACM SIGGRAPH ASIA2008courses.New York:ACM,1-11.[69]Yin,X.,Jin,M.,Luo,F.,et al,2009.Emerging trends in visual computing.In:.Berlin,Heidelberg:。

USTC中国科技大学《计算机图形学》暑期课程_计算机图形学进展_刘利刚

计算机图形学进展
刘利刚中国科学技术大学
我们的显示设备是二维的
2D图像
0.6 R + 0.3 G + 0.1 B
我们生活的世界是三维的
如何表达三维物体的表面？
陶瓷兔子实物的照片
3D模型：三维坐标点的表达方式
采用不同密度的点来表达的3D兔子模型
3D物体的数字表达：点和三角面
三角面片
线点
80
85
90
95
00
05
10
图像
几何
Visual Computing
• Visual media
– Computer graphics to visual information
Image
Video
Geometry
Animated geometry
计算机图形学的主要内容？
1. 建模 Modeling
3D设计（3D建模）
将想象中的3D物体构建出来
3D Modeling from X
Sketches
Hand drawing
Point Clouds
3D Model
Images
3D Scanner
Taking Photos
1.1 3D扫描建模
3D数字采集设备
基于深度相机的3D扫描
通过深度相机设备测量物体表面的点坐标数据，然后通过数学算法构建3D模型
3D打印
3D扫描 3D动画 3D设计广泛的行业应用 3D娱乐 3D交互
3D科技
3D科技：以3D数据为表达信息以3D处理为技术手段
城市规划
3D交互：手势、身体等
3D交互：虚拟头盔+全息影像

黎曼几何简介-USTC

λµeν
∇fα fβ = Γ˜γ αβ fγ
=
Γ˜ γ
∂xν αβ ∂yγ eν
上面两个结果相对照，得到：
Γ˜ γ
αβ
=
∂xλ ∂yα
∂xµ ∂yβ
∂yγ Γν ∂xν
λµ
+
∂2xν ∂yα∂yβ
∂yγ ∂xν
3 度规相容联络
由于联络系数的具体形式在上面并未赋予，我们进一步从物理角度提出要求：在 V 上平行移动的矢量 X, Y 的矢量积 g(X, Y ) 应不变：
3 度规相容联络
4
也有莱布尼茨法则：
∇X (f Y ) = X[f ]Y + f ∇X Y
∇X (T1 ⊗ T2) = (∇X T1) ⊗ T2 + T1 ⊗ (∇X T2)
协变导数在坐标变换下，联络系数的变换：取
2
组坐标
{eµ}
=
{
∂
∂ xµ
},
{fα}
=
{
∂
∂ yα
},
对后者的协变导数：
∇fα fβ = Γ˜γ αβ fγ
特征值的数目，其中 j = 1 是 Lorentz 度规，总可坐标伸缩为闵氏度规 η = diag(−1, 1, · · · , 1)
以后称呼流形和其上的度规组 (M, g) 是相应的黎曼或洛伦兹型的；在洛伦兹型的流形上，将 TpM 分成三个子空间：g(U, U ) > 0 是 spacelike 的， = 0 是 lightlike 的，< 0 是 timelike 的。
6. 嵌入度规 (Induced metric) 是子流形嵌入到维度更大流形上，获得的度规结构：若 M 是 n 维 Riemann 流形 N 的 m 维子空间，有映射 f : M → N 将 M 嵌入到 N 中，（N 的度规已知），则：

计算机图形学-USTC

• J. D. Foley et al. Computer Graphics - Principles and Practice, 3nd Ed., 2013.(影印版: 计算机图形学原理及实践-C语言描述(第二版), 机械工业出版社, 2002)
• D. F. Rogers. Procedural Elements for Computer Graphics, 2nd Ed., 1译9版98).(影印版: 计算机图形学的算法基础, 机械工业出版社, 2002. 有中
述的图形等等
• 构成图形的要素
– 几何要素：刻画对象的轮廓、形状的点、线、面、体等
– 非几何要素：刻画对象的颜色、材质、纹理、线型、线宽等
11
表示图形的方法
• 参数表示
– 由图形的形状参数(方程或分析表达式的系数，线段的端点坐标等)+属性参数(颜色、线型等) 来表示图形
– 称为参数图或图形（graphics）
• 中译版: OpenGL编程指南 (原书第6版), 机械工业出版社, 2008.8
• 在线OpenGL 1.1版本 http://www.glprogramming.c om/red/
5
OpenGL
• R. S. Wright, et al, OpenGL SuperBible: Comprehensive Tutorial and Reference (4th Edition), 2007.
• 点阵表示
– 枚举出图形中所有的点(强调图形由点构成) – 称为像素图或图像（image, 数字图像）
12
图形与图像
• 图形含有几何属性，或者说更强调场景的几何表示，是由场景的几何模型和物理属性共同描述

计算共形几何

书中还介绍了一些有趣的共形几何现象。例如，通过共形变换可以将一个平面分成有限个区域，每个区域都是一个凸多边形，这一过程被称为凸包映射。这一现象在地图制作、网络路由等领域有着广泛的应用。
《计算共形几何》这本书的精彩摘录展示了共形几何的广泛应用和重要价值。通过深入了解共形几何的基本概念和原理，我们可以更好地理解数学和物理中的一些基本概念和规律，并更好地应用这些知识来解决实际问题。
计算共形几何
读书笔记
01 思维导图
03 精彩摘录 05 目录分析
目录
02 内容摘要 04 阅读感受 06 作者简介
思维导图
本书关键字分析思维导图
变换
这些
几何
读者
数学
理解
保持
计算
共形
几何重要
几何量
共形
深入
基础
不变
书籍
性质
习题
内容摘要
内容摘要
《计算共形几何》是一本深入探讨共形几何理论的书籍，其内容涵盖了从基础概念到最新进展的全面知识。共形几何作为数学的一个重要分支，主要研究图形在连续变换下的不变性质，这种变换可以想象为把纸张扭曲、折叠、拉伸或收缩，但保持其边界形状不变。这本书首先从共形几何的基础知识开始，解释了共形变换和几何量的定义、性质以及基本的定理。它详细阐述了共形变换如何在不同的维度和曲面上保持几何量不变，以及这些性质在几何学、拓扑学和物理学中的重要应用。然后，这本书深入探索了共形几何的高级主题，包括黎曼面、复分析、非线性微分方程和流形上的几何。这些主题的讨论为读者提供了理解共形几何在现代数学和物理中的重要性的深入视角。这本书还包含了大量的习题和示例，以帮助读者理解和掌握共形几何的概念。
内容摘要

Note: This is mostly copied or translated from the papers ofProf. Xianfeng Gu (/~gu/).计算共形几何张威应用数学110060591引言共形几何是纯数学中很多学科的交叉领域,比如黎曼曲面理论、微分几何、代数曲线、代数拓扑、偏微分方程、复分析等等.它有很长的历史,至今在现代几何与现代物理中仍然非常活跃.比如超弦理论中的共形场和理论物理中的模空间理论都是当今快速发展的研究领域.近些年来,随着三维数字扫描仪、计算机辅助几何设计、生物信息和医学成像的快速发展,出现了越来越多的三维数字模型.因此迫切需要有效的算法来表示、处理和使用这些模型.计算共形几何在数字几何处理中扮演了一个重要角色.它已经应用在很多重要的领域,比如曲面修复、光顺、去噪、分片、特征提取、注册、重新网格化、网格样条转换、动画和纹理合成.特别地,共形几何奠定了曲面参数化的理论基础,同时也提供了严格的算法.计算共形几何还应用于计算机视觉中的人脸跟踪、识别和表情转换,医学成像中的脑电图、虚拟结肠镜和数据融合,几何建模中的具有任意拓扑流形上的样条构造.共形几何之所以如此有用是基于以下一些事实:•共形几何研究的是共形结构.日常生活中的所有曲面都有一个自然的共形结构,因此共形几何算法非常普遍.•共形结构比黎曼度量结构更灵活、比拓扑结构更具有刚性.它能处理大量黎曼几何不能有效处理的变换,这些变换还能保持很多拓扑方法会丢失的几何信息.•共形映射比较容易控制.比如,两个单连通封闭曲面之间的共形映射构成一个6维空间,因此只要固定3个点,这个映射就是唯一的.这个事实使得共形几何方法在曲面匹配和比较中非常有价值.•共形映射保持局部形状,因此在可视化方面有很好的应用.•所有的曲面都可以根据共形结构进行分类,而且所有的共形等价类形成一个有限维流形.这个流形有丰富的几何结构,容易对其分析和研究.与之相反,曲面的等距类形成一个难以分析处理的无穷维流形.计算共形几何•计算共形几何算法是以椭圆偏微分方程为基础的,而椭圆偏微分方程又是容易求解而且稳定的,因此计算共形几何方法对于实际工程应用非常有用.•共形几何中,所有单连通曲面都能共形变换成某种标准空间:球面、平面、双曲空间.也就是说,任何曲面都具有三种标准几何(球几何、欧式几何、双曲几何)中的一种.这样大部分三维数字几何处理任务都能转化成二维标准空间中的任务.历史上,计算共形几何方法已经广泛应用于许多工程领域,然而绝大部分应用都基于平面区域的共形映射.最近,随着数学理论的发展和计算能力的提升,计算共形几何方法已经从平面区域推广到具有任意拓扑的曲面.1.1共形变换和共形结构图1:共形映射保持角度(a)Circle packing(b)Checkboard图2:共形映射根据Fleix Klein的Erlangen纲领,几何就是研究在特定的变换群下保持不变的空间性质.共形几何就是研究保角变换群下的不变量.它介乎于拓扑和黎曼几何之间.计算共形几何共形映射就是保角映射,如图1所示.在无穷小邻域,共形映射就是放缩变换.它保持局部形状,比如它将无穷小圆周映成无穷小圆周.如图2所示,这个bunny曲面通过一个共形映射映到平面.如果平面有一个circle packing,则通过拉回得到bunny曲面上的一个circle packing.如果给平面铺上棋盘格,则同样得到bunny曲面的棋盘格修饰,其中直角和正方形都是保持的.曲面上的两个黎曼度量是共形的,如果它们定义的角度是相同的.共形结构就是指曲面上度量的共形等价类,而黎曼曲面就是带有共形结构的光滑曲面.因此在黎曼曲面上,我们可以度量角度,但不能度量长度.每一个带有度量的曲面都自动成为一个黎曼曲面.如果两个黎曼曲面之间存在共形映射,则称它们是共形等价的.显然,共形等价是黎曼曲面间的一个自然的等价关系.共形几何的目的就是在共形等价意义下对黎曼曲面进行分类,这就是所谓的模空间问题.给定一张光滑曲面,考察它上面的所有共形结构在共形等价下的模,这个集合被称为曲面的模空间.对于具有正亏格的封闭曲面,模空间是正维数的有限维空间.1.2基本任务下面的问题是计算共形几何最基本的一些任务.这些问题是相互依赖的:1.共形结构给定一张带有黎曼度量的曲面,计算它的内蕴共形结构的不同表示.一种方法是计算它的Abelian微分群,另一种方法是计算标准的黎曼度量.2.共形模完全共形不变量称为黎曼曲面的共形模.正如前面所讲的,理论上存在一组有限的数完全决定了黎曼曲面,这些称为黎曼曲面的共形模.一个比较难的问题是显式计算任意给定曲面的共形模.3.标准黎曼度量黎曼曲面的uniformization定理揭示了每一个黎曼度量都共形等价于一个常Gauss 曲率度量.除了球面和环面外,这个度量是唯一的.计算这个度量在计算共形几何中具有基本的重要性.4.共形映射计算共形几何计算两个共形等价的曲面之间的共形映射可以简化为计算它们到标准形状空间(球面、平面、双曲空间中的圆域)之间的共形映射.5.拟共形映射大部分微分同胚都不是共形的,它们将无穷小圆周映成无穷小椭圆.如果这些椭圆的长短轴比一致有界,那么就称为拟共形映射.拟共形映射的微分是由所谓的Beltrami微分刻画的,Beltrami微分记录了长轴方向和长短轴比.有一个基本定理是说通过Beltrami微分可以恢复拟共形映射.而至于怎样通过Beltrami微分计算拟共形映射,则是一个具有很多应用价值的重要问题.6.共形镶嵌粘合带边黎曼曲面并研究缝曲线形状和粘合样式之间的关系.这与拟共形映射问题紧密相关.1.3共形几何在工程应用中的优点计算共形几何已经被证实在许多工程领域中有重要应用.下面是一些主要理由:1.标准区域所有带度量的曲面都能共形地映成球面、平面或双曲圆盘中的标准区域.这可以帮助我们将三维几何处理问题转化为二维问题.2.通过曲率设计度量每一个共形结构都有一个常Gauss曲率的标准度量,这个度量在很多几何应用中非常有价值.例如在双曲度量下,每一个非平凡同伦类都有一个闭测地线代表元.此外,我们可以根据预先给定的曲率设计黎曼度量,这在几何建模中非常有用.3.一般几何结构共形几何方法能够用来构造其他几何结构,比如仿射结构、射影结构等等.这些结构在几何建模应用中是关键的.4.微分同胚的构造共形映照和拟共形映照可以被用来构造曲面间的微分同胚.可以应用于曲面注册和比较这些计算机视觉和医学图像中最基本的问题.计算共形几何5.等温坐标共形结构可以被当作曲面上的等温坐标图册.在这种坐标下,度量的表达式最简单,因此所有的微分算子,比如Laplace-Beltrami算子,具有很简洁的表达式.这可以用来简化偏微分方程.等温坐标保持局部形状,对于可视化和纹理映射是非常完美的.2已有工作图3:基本群的多边形表示图4:四种类型的图.(a)A cut graph但不是system of loops.(b)A system of loops.(c)基本群的基但不是system of loops.(d)同调基,但既不是同伦基也不是cut graph.计算共形几何已经有比较长的历史了,起初大部分共形几何方法都是针对平面区域或拓扑圆盘的.直到近些年来才推广到具有复杂拓扑的曲面上,这与当前研究热点计算拓扑的发展是紧密相关的.有相当一部分计算共形算法是依赖于计算拓扑算法的.计算共形几何图5:从左到右:输入模型,handle(绿)和tunnel(红)环,handle特征2.1计算拓扑计算拓扑的研究对象主要是单纯复形,比如三角网格曲面就是一种二维单纯复形.计算拓扑的一些基本问题就是计算曲面单纯复形的cut graph,同调群和基本群.知道了cut graph,就能将曲面剪开成拓扑圆盘,这对于参数化和纹理映照是非常有用的.知道了基本群,就能完全确定封闭曲面的拓扑分类.通常基本群是用多边形表示的,其中有约化多边形表示、标准多边形表示等一些特殊的多边形表示,见图3.Vegter和Yap[1]给出了计算标准多边形表示的算法,Lazarus等[2]简化了Vegter和Yap的算法,复杂度为O(gn),其中n,g分别为曲面的大小和亏格.标准多边形表示由2g 个具有公共点的环组成,令人不满意的是它们可能具有公共边.在大部分应用中,标准多边形表示不是必须的,只需要知道cut graph就足够了.cut graph是指三角网格曲面上的一族边集使得曲面去掉这些边后变成拓扑圆盘,如图4所示.Dey和Schipper[3]使用cut graph得到了检验曲线是否可缩或两曲线是否同伦的快速算法.2002年,Erickson等[4]第一次提出了最优cut graph的问题:给每一条边赋一个权,求权最小的cut graph.这篇文章证明了最小cut graph问题是NP难的,并给出了计算近似最小cut graph的贪婪算法,该算法复杂度为O(g2n log n),精度为O(log2g).de Verdière和Lazarus[5]研究了一种称为"systerm of loops"的特殊cut graph,每一个systerm of loops都由具有公共基点的2g个环组成,如图4所示.Erickson等[6]给出了一个简单的贪婪算法来计算基本群和第一同调群的生成元,这同时也解决了de Verdière和Lazarus提出了最短systerm of loops问题.计算共形几何曲面上一些与特征紧密相关的环,比如handles和tunnels,在诸如拓扑修复、曲面参数化和特征识别等应用中非常有用,如图5所示.Dey等[7--9]给出了计算handles和tunnels的算法.Yin等[10]提出了使用universal covering space来计算最短非平凡环的算法.de Verdière[11]使用最优pants分解来找同伦类中的最短环.Xin等[12]研究了高亏格曲面上的测地环.关于曲面上的非平凡环,还有许多文章[13--19],目前仍是一个热点问题.(a)调和1形式(b)共轭的调和1形式图6:调和1形式2.2平面区域传统的数值复分析方法集中在平面区域的共形映射.[20--25]作了详尽总结. Schwarz-Christoffel映射已经被广泛应用于计算共形映射[26,27].[28--30]研究了多连通区域的Schwarz-Christoffel映射.[31]给出了一种基于交比和Delaunay三角化的鲁棒算法.最近,Marshall等[32]引进了基于迭代简单映射的geodesic zippper算法,Bishop[33]提出了基于双曲几何的线性复杂度的共形映射算法.2.3亏格为0的曲面[34]利用余切公式构造了离散调和映射,Lévy等[35]引进了Cauchy-Riemann方程的一阶有限元逼近.2002年,Desbrun等[36]最小化Dirichlet能量得到了离散内蕴参数化. Floater[37]引进平均值坐标来计算广义调和映射.[38,39]对拓扑球面的共形映射进行了研究.在CG领域,[40,41]对拓扑圆盘和球面上的共形映射进行了详尽地总结.计算共形几何2.4高亏格曲面高亏格曲面上共形结构的计算有两个主要方法:全纯微分方法和离散曲率流方法.图7:全纯1形式2.4.1全纯微分离散全纯微分方法是由Gu和Yau[42,43]提出的,用来计算高亏格曲面上的共形结构.这个方法基于Hodge理论,使用热扩散方法来计算每一个上同调类中的调和形式,见图6(a)和6(b).然后利用Hodge星算子构造全纯形式,见图7.所有计算都是在离散多边形曲面上进行的.Pinkall等[34]定义了一种不同的离散Hodge星算子来计算极小曲面. Mercat[44]通过离散Cauchy-Riemann方程来构造离散全纯映射,不过这个方法要求曲面是四边形网格.2003年,Hirani在他的博士论文[45]中详细介绍了离散外微分方法.Gortler等[46,47]使用离散1形式来对亏格为1的曲面进行参数化.Tong等[48]推广到带有锥奇点的1形式方法,并用来进行remeshing和tiling.2009年,Zeng等[49,50]将全纯微分方法应用于计算带多个边界的亏格为0的曲面上的共形映射,以及拟共形映射.2.4.2Ricci FlowRicci流是由Hamilton[51]在Princeton的讨论班上引进的.Ricci流对曲面和三维流形的几何研究产生了革命性的影响,是当前几何研究中极为活跃的方向.特别地,计算共形几何它导致了三维Poincaré猜测的证明.Hamilton[52]使用2维Ricci流给出了正亏格曲面uniformization定理的证明.这表明了它在CG领域中有极大的潜在应用价值.图8:Circle packing存在很多方法离散光滑曲面,其中与共形映射离散特别相关的是由Thurston[53]引进的circle packing度量.circle packing这个概念首先是由Koebe提出的,Thurston猜测平面Jordan区域的circle packing离散序列收敛于Riemann映射,见图8.这个猜测被Rodin 和Sullivan[54]所证明.Colin de Verdière[55]建立了circle packing的第一变分原理,然后证明了Thurston的circle packing度量存在性.这为实现计算circle packing度量的快速算法铺平了道路,比如Collins和Stephenson[56]就给出了一个算法.Chow和Luo[57]推广了Colin de Verdière 的工作,并引进了曲面上的离散Ricci流和离散Ricci能量.它们证明了离散Ricci流的一般存在性和收敛性定理,同时证明了Ricci能量的凸的.基于此,Jin[58]在2008年给出了离散Ricci流的算法实现.另外一个相关的离散方法是circle pattern,它同时考虑网格曲面的组合与几何信息,可以看成是circle packing的变种.Circle pattern是由Bowers和Hurdal[59]提出的,已经被证明是Bobenko和Springborn[60]提出的凸能量的极小元.2006年,Kharevych等[61]给出了一个有效的circle pattern算法.2.4.3Yamabe FlowYamabe问题目的是为紧Riemann流形找到具有常标量曲率的共形度量.Yamabe[62]给出了第一个带有漏洞的证明,之后被一些包括Trudinger[63],Aubin[64]和Schoen[65]在内的研究者所更正.关于这个问题的详细总结可以参看Lee和Parker的文章[66].2004年,Luo[67]研究了曲面上的离散Yamabe流.他引进了多边形度量的离散共形计算共形几何变换概念,这在研究离散Yamabe流和相应的变分原理方面具有关键性的作用.基于离散共形类和几何方面的考虑,Luo定义离散Yamabe能量为微分1形式的积分,然后证明了这个能量是一个局部凸函数.他还从这得出Yamabe流下的曲率发展是一个热方程.最近,Springborn等[68]的一个非常漂亮的工作表明可以将Yamabe能量和Milnor-Lobachevsky函数等同起来,将曲率发展的热方程和余切Laplace方程等同起来.近期Gu 等[69]使用离散Yamabe能量构造了一个计算离散共形度量的有效算法.Bobenko等[70]讨论了离散双曲Yamabe流,它被Zeng等[71]用来计算双曲结构和标准同伦群生成元.图9:全局共形的纹理映照图10:使用特殊的平坦度量设计向量场2.5应用2.5.1图形学共形几何方法广泛应用于CG中.等温坐标对于全局共形参数化来说是很自然的[42].因为共形映照不会扭曲局部形状,所以它对于纹理映照来说是非常理想的.图9显示了一个使用全纯微分进行纹理映照的例子.特殊的平坦度量对于设计曲面上的向量场是非常有价值的,图10显示了使用曲率流方法设计的向量场.计算共形几何图11:流形样条框架2.5.2几何建模几何建模中的一个最基本的问题就是将传统的欧式区域上的样条系统地推广到流形区域上,这与曲面上的一般几何结构关系紧密,见图11和12.传统的样条方案是建立在仿射不变量上的,如果流形具有仿射结构,那么传统的样条也能直接推广到流形上.然而由于拓扑障碍,一般的流形不具有仿射结构.但是删去一些奇点后,曲面上允许存在仿射结构.详细的讨论请参考Gu等[72].仿射结构可以通过共形几何方法显式计算.比如可以通过曲率流和全纯微分方法来计算仿射结构,图13显示的是用曲率流方法构造的仿射结构.与其他的方法比较起来,比如基于平凡联络的方法[73],这些方法的优点是它能给出样条曲面的全局共形参数化,即等温坐标.在这种坐标下,微分算子,比如梯度算子和Laplace-Beltrami算子,具有最简单的形式.2.5.3医学成像共形几何在医学成像的很多领域都有应用.比如在脑电图中,不同部分的大脑皮层曲面片的注册是关键性的.大脑表面高度卷曲,而且不同的人的大脑有不同的解剖结构,计算共形几何图12:流形样条图13:仿射结构因此在大脑皮层曲面片之间找到好的匹配是非常有挑战性的任务.图14演示了一个解决方案[39],先通过用标准方式将大脑表面映到球面,然后找到球面的自同构,这样曲面片的注册就很容易建立了.2.5.4视觉曲面匹配在计算机视觉中是一个基本问题.图15演示了曲面匹配的基本框架.图16演示了同一个人的不同表情的人脸匹配的例子,细节可参考[74--76].Teichmüller理论可以应用于曲面分类[77,78].通过使用Ricci曲率流,我们可以计算双曲uniformization度量.然后我们使用测地线进行pants分解并计算Fenchel-Nielsen坐标.图17演示了万有覆盖空间在双曲空间中的一个有限部分的计算,图18演示了计算Teichmüller坐标的流程.计算共形几何图14:大脑表面的共形映照图15:曲面匹配框架2.5.5计算几何在计算几何中,同伦检测是一个重要的问题.可以使用Ricci流计算双曲uniformiza-tion度量[71].根据Gauss-Bonnet定理,每一个同伦类中有唯一的闭测地线.给定一个环,我们可以计算这个环的同伦类对应的Möbius变换,这个变换的轴就是双曲度量下的闭测地线.如图19所示,同伦的环的标准代表元是相同的.计算共形几何图16:不同表情的人脸匹配图17:万有覆盖空间在双曲空间中的有限部分3总结计算共形几何是一门数学和计算机科学的交叉学科,它对数学理论的掌握要求比较高.很多算法都依赖于拓扑算法,与计算拓扑紧密相关.本文介绍了这个方向的一些基本问题、算法以及应用.目前仍然有很多基本问题没有解决,它们需要更深刻的理解和更强大精确的计算方法.下面列举一些对理论和应用都非常重要的问题.1.Teichmüller映射给定两个度量曲面和映射同伦类,计算唯一的angle distortion最小的映射,也就是Teichmüller映射.计算共形几何图18:Teichmüller空间中的Fenchel-Nielsen坐标图19:使用双曲度量进行同伦检测2.Abel微分计算不同类型的Abel微分群,特别是全纯的二次微分.3.组合结构和共形结构之间的关系给定一个拓扑曲面,每一个三角化都有一个由circle packing决定的自然的共形结构,讨论这两个结构之间的关系.4.逼近理论尽管计算共形不变量的算法已经有了,但是逼近的理论结果仍没有建立.对于平面区域的共形映射,不同离散方法的收敛性都已经建立了.但是对于一般曲面而言,收敛性分析仍然是开放的.5.精度和稳定性双曲几何计算对于数值误差是非常敏感的.要提高计算精度是非常具有挑战性的.计算几何中的精确计算方法给出了解决这个问题的方向.在逆距离circle packing计算共形几何方法和组合Yamabe流方法中,可接受曲率空间的非凸性造成了算法的不稳定.因此,曲面需要高精度的三角化.在实际应用中,提高三角化精度对这些算法是非常重要的.交角为锐角的circle packing算法更加稳定,全纯微分方法最稳定.计算共形几何参考文献[1]Vegter,G.,Yap,C.K.,putational complexity of combinatorial surfaces.In:Proceedingsof the sixth annual symposium on Computational geometry.New York:ACM,102-111.[2]Lazarus,F.,Pocchiola,M.,Vegter,G.,et al,puting a canonical polygonal schema of anorientable triangulated surface.In:SCG'01:Proceedings of the seventeenth annual symposium on Computational geometry.New York:ACM,80-89.[3]Dey,T.K.,1994.A new technique to compute polygonal schema for2-manifolds with applica-tion to null-homotopy detection.In:Proceedings of the tenth annual symposium on Computational geometry.New York:ACM,277-284.[4]Erickson,J.,Har-Peled,S.,2002.Optimally cutting a surface into a disk.In:SCG'02:Proceedingsof the eighteenth annual symposium on Computational geometry.New York:ACM,244-253.[5]de Verdière,E.C.,Lazarus,F.,2002.Optimal system of loops on an orientable surface.In:Proc.43rd Annual IEEE Symp.Foundations of Computer Science.Springer-Verlag627-636.[6]Erickson,J.,Whittlesey,K.,2005.Greedy optimal homotopy and homology generators.In:SODA'05:Proceedings of the sixteenth annual ACM-SIAM symposium on Discrete algorithms.Philadel-phia,PA,USA:Society for Industrial and Applied Mathematics,1038-1046.[7]Dey,T.K.,Li,K.,Sun,J.,2007.On computing handle and tunnel loops.In:CW'07:Proceedingsof the2007International Conference on Cyberworlds.Washington,DC,USA:IEEE Computer Society,357-366.[8]Dey,T.K.,Li,K.,Sun,J.,et al,puting geometry-aware handle and tunnel loops in3dmodels.In:SIGGRAPH'08:ACM SIGGRAPH2008papers.New York:ACM,1-9.[9]Dey,T.K.,Li,K.,Sun,J.,puting handle and tunnel loops with knot put.Aided Des.,41(10):730-738.[10]Yin,X.,Jin,M.,Gu,X.,puting shortest cycles using universal covering space.Vis.Comput.,23(12):999-1004.[11]Verdière,E.C.D.,Lazarus,F.,2007.Optimal pants decompositions and shortest homotopic cycleson an orientable surface.J.ACM,54(4):18.[12]Xin,S.Q.,He,Y.,Fu,C.W.,et al,2011.Euclidean geodesic loops on high-genus surfaces appliedto the morphometry of vestibular systems.In:MICCAI'11:Proceedings of the14th international conference on Medical image computing and computer-assisted intervention.Berlin,Heidelberg: Springer-Verlag,384-392.[13]Chambers,E.W.,Erickson,J.,Worah,P.,2008.Testing contractibility in planar rips complexes.In:SCG'08:Proceedings of the twenty-fourth annual symposium on Computational geometry.New York:ACM,251-259.[14]Chambers,E.W.,Erickson,J.,Nayyeri,A.,2009.Minimum cuts and shortest homologous cycles.In:SCG'09:Proceedings of the25th annual symposium on Computational geometry.New York: ACM,377-385.[15]Chambers,E.W.,Erickson,J.,Nayyeri,A.,2009.Homology flows,cohomology cuts.In:STOC计算共形几何'09:Proceedings of the41st annual ACM symposium on Theory of computing.New York:ACM, 273-282.[16]Erickson,J.,2010.Maximum flows and parametric shortest paths in planar graphs.In:SODA'10:Proceedings of the Twenty-First Annual ACM-SIAM Symposium on Discrete Algorithms.Philadelphia,PA,USA:Society for Industrial and Applied Mathematics,794-804.[17]Erickson,J.,Worah,P.,puting the shortest essential cycle.Discrete Comput.Geom.,44(4):912-930.[18]Dey,T.K.,Hirani,A.N.,Krishnamoorthy,B.,2010.Optimal homologous cycles,total unimodular-ity,and linear programming.In:STOC'10:Proceedings of the42nd ACM symposium on Theory of computing.New York:ACM,221-230.[19]Erickson,J.,2011.Shortest non-trivial cycles in directed surface graphs.In:SoCG'11:Proceedingsof the27th annual ACM symposium on Computational geometry.New York:ACM,236-243.[20]Delillo,T.K.,1994.The accuracy of numerical conformal mapping methods:A survey of examplesand results.SIAM Journal on Numerical Analysis,31(3):pp.788-812.[21]Ivanov,M.,V.T,1995.Handbook of conformal mapping with computer-aided visualization.BocaRaton:CRC Press.[22]Binder,I.,Braverman,M.,Yampolsky,M.,2007.On the computational complexity of the riemannmapping.ARKIV FOR MATEMATIK,45(2):221-239.[23]Trefethen,L.,1986.Numerical conformal mapping.Amsterdam:North-Holland Publishing Co.[24]Wegmann,R.,2005.Handbook of Complex Analysis:Geometric Function Theory.vol2.Ams-terdam:Elsevier,351-477.[25]Henrici,,1993.Discrete Fourier Analysis,Cauchy Integrals,Construction of Conformal Maps,Univalent Functions.vol3.London:Wiley.[26]Driscoll,T.A.,Trefethen,L.N.2002.Schwarz-Christoffel Mapping.Cambridge Press.[27]Banjai,L.,Trefethen,L.N.,2003.A multipole method for schwarz--christoffel mapping of poly-gons with thousands of sides.SIAM put.,25(3):1042-1065.[28]DeLillo,T.K.,Elcrat,A.R.,Pfaltzgraff,J.A.,2001.Schwarz-christoffel mapping of the annulus.SIAM Review,43(3):pp.469-477.[29]Delillo,T.,Elcrat,A.,Pfaltzgraff,J.,2004.Schwarz-christoffel mapping of multiply connecteddomains.Journal d'Analyse Mathématique,94:17-47.[30]Crowdy,D.,2005.The schwarz-christoffel mapping to bounded multiply connected polygonaldomains.Proceedings:Mathematical,Physical and Engineering Sciences,461(2061):pp.2653-2678.[31]Driscoll,T.A.,Vavasis,S.A.,1998.Numerical conformal mapping using cross-ratios and delaunaytriangulation.SIAM put.,19:1783-1803.[32]Marshall,D.E.,Rohde,S.,2007.Convergence of a variant of the zipper algorithm for conformalmapping.SIAM J.Numer.Anal.,45:2577-2609.计算共形几何[33]Bishop,C.J.,2010.Conformal mapping in linear time.Discrete Comput.Geom.,44:330-428.[34]Pinkall,U.,Polthier,K.,puting discrete minimal surfaces and their conjugates.Exper-imental Mathematics,2:15-36.[35]Lévy,B.,Petitjean,S.,Ray,N.,et al,2002.Least squares conformal maps for automatic texture atlasgeneration.In:Proceedings of the29th annual conference on Computer graphics and interactive techniques.New York:ACM,362-371.[36]Desbrun,M.,Meyer,M.,Alliez,P.,2002.Intrinsic parameterizations of surface puterGraphics Forum,21(3):209-218.[37]Floater,M.S.,2003.Mean value put.Aided Geom.Des.,20:19-27.[38]Gotsman,C.,Gu,X.,Sheffer,A.,2003.Fundamentals of spherical parameterization for3d meshes.In:SIGGRAPH'03:ACM SIGGRAPH2003Papers.New York:ACM,358-363.[39]Gu,X.,Wang,Y.,Chan,T.F.,et al,2004.Genus zero surface conformal mapping and its applicationto brain surface mapping.IEEE Transactions on Medical Imaging,23:949-958.[40]Floater,M.S.,Hormann,K.,2005.Surface parameterization:a tutorial and survey.In:Dodgson,N.A.,Floater,M.S.,Sabin,M.A..Advances in multiresolution for geometric modelling.Springer Verlag,157-186.[41]Sheffer,A.,Praun,E.,Rose,K.,2006.Mesh parameterization methods and their applications.put.Graph.Vis.,2:105-171.[42]Gu,X.,Yau,S.T.,2003.Global conformal surface parameterization.In:SGP'03:Proceedingsof the2003Eurographics/ACM SIGGRAPH symposium on Geometry processing.Aire-la-Ville, Switzerland,Switzerland:Eurographics Association,127-137.[43]Jin,M.,Wang,Y.,Yau,S.T.,et al,2004.Optimal global conformal surface parameterization.In:VIS'04:Proceedings of the conference on Visualization'04.Washington,DC,USA:IEEE Computer Society,267-274.[44]Mercat,C.,2001.Discrete riemann surfaces and the ising munications in MathematicalPhysics,218:177-216.[45]Hirani,A.N.,2003.Discrete exterior calculus.PhD thesis.Pasadena,CA,USA:California Instituteof Technology.[46]Gortler,S.J.,Gotsman,C.,Thurston,D.,2006.Discrete one-forms on meshes and applications to3d mesh put.Aided Geom.Des.,23:83-112.[47]Tewari,G.,Gotsman,C.,Gortler,S.J.,2006.Meshing genus-1point clouds using discrete one-puters and Graphics,30:917-926.[48]Tong,Y.,Alliez,P.,Cohen-Steiner,D.,et al,2006.Designing quadrangulations with discreteharmonic forms.In:Proceedings of the fourth Eurographics symposium on Geometry processing.Aire-la-Ville,Switzerland,Switzerland:Eurographics Association,201-210.[49]Zeng,W.,Yin,X.,Zhang,M.,et al,2009.Generalized koebe's method for conformal mappingmultiply connected domains.In:SPM'09:2009SIAM/ACM Joint Conference on Geometric and计算共形几何Physical Modeling.New York:ACM,89-100.[50]Zeng,W.,Luo,F.,Yau,S.T.,et al,2009.Surface quasi-conformal mapping by solving beltramiequations.In:Proceedings of the13th IMA International Conference on Mathematics of Surfaces XIII.Berlin,Heidelberg:Springer-Verlag,391-408.[51]Hamilton,R.,1982.Three manifolds with positive ricci curvature.J.Differ.Geom.,17:255-306.[52]Hamilton,R.,1988.Mathematics and General Relativity.vol71.American Mathematical Society.[53]Thurston,W.,1980.Geometry and Topology of Three-Manifolds.Princeton University Press.[54]Rodin,B.,Sullivan,D.,1987.The convergence of circle packings to the riemann mapping.J.Differ.Geom.,26(2):349-360.[55]Colin de Verdière,Y.,1991.Un principle variationnel pour les empilements de cercles.Invent.Math.,104(3):655-669.[56]Collins,C.R.,Stephenson,K.,2003.A circle packing put.Geom.Theory Appl.,25:233-256.[57]Chow,B.,Luo,F.,binatorial ricci flows on surfaces.J.Differ.Geom.,63(1):97-129.[58]Jin,M.,Kim,J.,Luo,F.,et al,2008.Discrete surface ricci flow.IEEE Transactions on Visualizationand Computer Graphics,14(5):1030-1043.[59]Bowers,P.L.,Hurdal,M.K.,2003.Planar conformal mappings of piecewise flat surfaces.In:In Vi-sualization and Mathematics III,H.-C.Hege and K.Polthier,Eds.Mathematics and Visualization.Springer-Verlag3-34.[60]Bobenko,A.I.,Springborn,B.A.,2004.Variantional principles for circle patterns and koebe'stheorem.Trans.Am.Math.Soc.,356:659-689.[61]Kharevych,L.,Springborn,B.,Schröder,P.,2006.Discrete conformal mappings via circle patterns.ACM Trans.Graph.,25:412-438.[62]Yamabe,H.,1960.The yamabe problem.Osaka Math.J.,12(1):21-37.[63]Trudinger,N.,1968.Remarks concerning the conformal deformation of riemannian structures oncompact manifolds.Ann.Scuola Norm.Sup.Pisa,22(2):265-274.[64]Aubin,T.,1976.équations diffréntielles non linéaires et problème de yamabe concernant la cour-bure scalaire.J.Math.Pure Appl.,55(3):269-296.[65]Schoen,R.,1984.Conformal deformation of a riemann metric to constant scalar curvature.J.Differ.Geom.,20(2):479-495.[66]Lee,J.M.,Parker,T.H.,1987.The yamabe problem.Bull.Am.Math.Soc.,17(1):39-71.[67]Luo,F.,Khan,L.,Bastani,F.,et al,2004.A dynamically growing self-organizing tree(dgsot)forhierarchical clustering gene expression profiles.Bioinformatics,20(16):2605-2617.[68]Springborn,B.,Schröder,P.,Pinkall,U.,2008.Conformal equivalence of triangle meshes.In:SIGGRAPH Asia'08:ACM SIGGRAPH ASIA2008courses.New York:ACM,1-11.[69]Yin,X.,Jin,M.,Luo,F.,et al,2009.Emerging trends in visual computing.In:.Berlin,Heidelberg:。

数学分析视频教程-全套220讲-史济怀-中国科技大学

数学分析视频教程全套220讲史济怀中国科技大学国家精品课程-中国科技大学数学分析视频222讲中科大数学分析史济怀8DVD赠pdf格式课件和部分期末考试试卷一、所用教材《数学分析教程》(上、下册),常庚哲,史济怀编，高等教育(2003年)二、章节容数学分析一77讲数学分析二88讲数学分析三55讲目前,本课程使用的教材是由我校数学系常庚哲和史济怀两位教授编著的《数学分析教程》上下册（高等教育，2003年5月，第一版）。

该教材是普通高等教育“十五”国家级规划教材，是在1998年教育出版的《数学分析教程》的基础上写成的，原书融合了20多年来数学系讲授数学分析课程的教师的教学经验，同时也参考了国外同类书籍中的许多名著，在全国同类教材中有非常积极的影响。

该教材已经在本校数学系使用了5年，教学效果很好。

该教材的第二版正在修订中。

参考书：1．《数学分析》，何琛，史济怀，徐森林编，高等教育(1985年)。

2．《数学分析新讲》，筑生编，大学(1991年)。

第一学期：主要讲授单变量函数的微积分学。

主要容有：实数理论，极限理论，单变量函数的微分学和积分学。

教学重点：极限理论，导数的概念和运算，Taylor公式，可积性理论和积分的计算。

教学难点：实数理论，极限理论，上、下极限，Taylor公式，可积性理论。

教材：《数学分析教程》(上册)，常庚哲，史济怀编，高等教育(2003年)。

参考书：《数学分析新讲》，筑生编，大学(1991年)。

第一章实数15学时§1 无尽小数1学时§2 收敛数列及其性质5学时§3 收敛原理和上下确界5学时§4 上、下极限和Stolz定理4学时第二章函数的连续性19学时§1 集合的映射和势2学时§2 函数的极限6学时§3 连续函数7学时§4 混沌现象4学时第三章函数的导数15学时§1 导数的定义和计算5学时§2 微分学中值定理及其应用5学时§3 凸函数及函数作图5学时第四章Taylor定理6学时第五章插值与逼近初步5学时第六章求导的逆运算5学时说明：讲课共用80学时，余下的学时用作习题课和期中考试。

清华笔记：计算共形几何讲义（15）拓扑圆盘的调和映照

清华笔记：计算共形几何讲义（15）拓扑圆盘的调和映照图1. 从三维人脸曲面到平面圆盘的调和映照。

我们前面的课程介绍了亏格为0带边界曲面的典范共形映射，包括古典理论和基于全纯微分的计算方法。

这里，我们着重介绍拓扑球面到标准球面的典范共形映射，主要是基于调和映照理论。

简介物理上，调和映射极小化弹性形变势能，因而物理意义明确；偏微分方程理论证明了调和映照的存在性，唯一性，正则性，稳定性和光滑性；有限元方法保证了离散解到连续解的收敛性；数值计算方法的共轭梯度法保证了调和映照计算的高效性；微分几何保证了调和映射的共形不变性和微分同胚性。

因此，调和映照简单直观，理论完备，在工程实践中被广泛应用。

同时，调和映照是一个极好的例子，从理论到实践，横跨物理学，偏微分方程理论，有限元理论，数值计算，微分几何等诸多领域，使我们能够融会贯通，体会到这些领域各有侧重，同时相辅相成的关系。

如图1所示，人脸曲面到平面区域的调和映射，可以表示成坐标映射，这里是相互独立的调和函数。

如果坐标分量函数彼此非独立而共轭，则调和映照称为共形映照。

下面，我们首先考察调和函数。

物理解释黎曼考察了这样一个物理问题：（Dirichlet问题）假设是平面中的有界区域，由某种电阻率处处相同的材料制成。

我们在的边缘处设置电压，问内部电压函数是多少？根据物理定律，内的电场诱导电流，电流发热做功，那么真实可能的电压函数必定使得发热功率最小。

电流强度正比于电压梯度，电阻率处处相同，因此电流发热功率可以表示成所谓的调和能量：，如果函数极小化调和能量，则我们称其为调和函数。

我们进一步考察调和函数应该满足的条件。

令试探函数为在边界上取值为0的无穷阶光滑函数。

假设，则对一切，的调和能量在为0的时候取到极值，因此，由关系式，我们得到由斯托克斯定理和，我们有，因为h任意，因此我们得到调和函数的欧拉-拉格朗日方程，即所谓的Laplace方程。

这里Laplace算子的物理意义是梯度的散度。

基于共形几何代数的C-球壳聚类方法及其实现

基于共形几何代数的C-球壳聚类方法及其实现李茂宽;刘超【摘要】提出一种新的基于共形几何代数的C-球壳聚类算法.首先C-球壳聚类分析中的模式和球壳原型,在共形几何代数中可以表示为矢量,而共形空间中模式与球壳的内积,等价于模式与球壳的距离,从而能够有效地表示球壳聚类分析中的目标函数,并且推导给出求解该目标函数最小化的新方法;介绍该算法基于VB和MatrixVB 的实现,该方法将Matlab的强大计算功能与VB的windows用户界面开发方面优势结合起来,充分发挥了各自优势,缩短了软件的开发周期.实验结果表明,该算法能够有效地对球壳形数据进行聚类,计算方法正确、计算速度快、系统资源消耗少,能满足数据聚类的要求.【期刊名称】《现代电子技术》【年(卷),期】2010(033)010【总页数】4页(P102-104,121)【关键词】共形几何代数;球壳聚类;圆检测;MatrixVB【作者】李茂宽;刘超【作者单位】海军航空工程学院,山东,烟台,264001;总参气象水文局,北京,100081【正文语种】中文【中图分类】TP391.4基于目标函数的聚类算法在模式识别、图像处理和计算机视觉中有着广泛的应用。

这种聚类方法是通过极小化一个目标函数得到数据的一种最优分类，目标函数聚类算法的研究始于20世纪60年代中期，最早的目标函数聚类算法是1966年G.H.Ball和D.J.Hall提出的ISODATA算法，即C均值聚类算法。

C均值聚类算法以点为原型，能够实现球形数据集的聚类。

随着聚类算法研究的不断深入，聚类原型也逐步由点扩展到线、面、球壳等。

Dave[1]最早开始壳聚类算法研究，提出了球壳状分布数据的聚类问题，Dave采用球壳作为聚类原型，定义目标函数，推导出球壳聚类的迭代算法，证明算法能有效地聚类球壳状分布的数据。

本文将介绍一种基于共形几何代数理论的C-球壳聚类算法，利用几何代数理论能有效地表示聚类分析中目标函数及其最小化求解问题。

计算共形几何

计算共形几何
共形几何就是把两个空间中的点，变为同一个空间中的点，因此它的定义是两点确定一条直线。

计算共形几何，我们可以用同一种方法把其他形式几何也统一起来。

“共形”一词来源于几何学，是指把空间分割为多个互相平行的空间，使这些空间之间的连接处变为平面。

在计算共形几何时，有很多情况需要把几何图形或几何量看成多个共形几何图形，由于共形几何与其他几何形式的相似性，所以当需要确定某些几何元素时，人们往往选择几何元素作为共形几何图形的模板，然后进行拼合，或者通过替换的方法将各种几何元素融入共形几何图形中去，从而得到最终的结果。

无论用什么方法计算都不可能得到光滑曲线，这是因为“不共形”（不光滑）的原因。

如果每一条曲线都是光滑的，那么一般方法也无法计算，这又是由于“不共形”的原因。

既然几何对象在外观上看上去不可能是光滑的，因此也就没有什么计算可言。

“光滑性”是一种计算上的技巧，但并不是真正的原因。

无论是“不共形”还是“光滑性”，它们都是在外观上看起来不可能是光滑的，也就是说，只要在计算上想办法使它们不是光滑的，就可以把计算过程简化了。

曲面的光滑性有时还要求把曲面切成许多小块，把这些小块看成共形几何图形。

但是，这个过程又太复杂，计算量太大。

因此人们在计算光滑曲线和曲面时，先假设这样的曲面被剪裁后会变得光
滑，再按照需要制作出曲面的剪裁模板，最后把剪裁下来的部分分别作为曲面的共形几何图形加以研究。

这样的做法对于那些虽然外观看起来不可能光滑，但是在计算上有可能被简化的曲面或曲线是十分有用的。

黎曼几何简介-USTC

3 度规相容联络
7
3. 李希张量 Ricci Tensor：
标曲率定义为：
Ricµν = Rλ µλν
R = gµν Ricµν
3.2 Levi-Civita Connection
{}
无挠率的联络显然是 Γλ µν =
λ , 满足 Γλ µν = Γλ νµ，只有对称 µν
的贡献。物理上经常计算这个联络，即假设流行上无挠率。本质上是因为这
Hodge Star 可以保证构造出满足上面条件的 m-form，即能写成 f ΩM 的形式。
5.1 Hodge Star
若 m 维流形 M 上有度规 g，可以定义 Ωr, Ωm−r 之间的同构映射，就是 Hodge Star 操作 (∗ : Ωr → Ωm−r)
√
∗(dxµ1
∧ dxµ2
∧ · · · ∧ dxµr ) =
−
1 gµν R) 2
=
∇µGµν
=
0
其中 G 是爱因斯坦张量，若这个张量和能动量张量成正比，就给出了
爱因斯坦场方程。
由以上的这些关系，m 维流形上的黎曼曲率张量的自由度有 F (m) =
1 12
m2(m2
−
1)
个，例如
2
维流形上只有一个自由度
R1212，且分量统一表示
成：
R Rkλµν = 2 (gkµgλν − gkν gλµ)
∇V [g(X, Y )] = V k[(∇kg)(X, Y ) + g(∇kX, Y ) + g(X, ∇kY )] = 0
3 度规相容联络
5
于是有条件 (∇kg)µν = 0，可化简为：
∂λgµν − Γk λµgkν − Γk λν gkµ = 0 满足这个条件的联络称为度规相容联络（metric connection）。对上式轮换指标并加和得到：

计算共形几何简介

计算共形几何简介
顾险峰;雷娜
【期刊名称】《大学数学》
【年(卷),期】2016(032)003
【摘要】计算共形几何是数学与计算机科学的交叉领域.本文介绍了计算共形几何中理论基础,计算方法及其在工程领域中的应用等.基本理论概念包括共形映射,单值化定理,共形模,拟共形映射,模空间等.基本计算方法包括调和映照方法,基于黎曼面理论的全纯微分方法,基于几何分析的曲面Ricci流方法.基本应用涵盖了医学图像,机械制造,数字媒体和虚拟/增强现实等领域.
【总页数】13页(P1-13)
【作者】顾险峰;雷娜
【作者单位】纽约州立大学石溪分校计算机系,纽约11794;大连理工大学软件学院,辽宁大连116620
【正文语种】中文
【中图分类】O186
【相关文献】
1.共形几何代数与几何不变量的代数运算 [J], 李洪波
2.共形几何代数--几何代数的新理论和计算框架 [J], 李洪波
3.五峰山长江特大桥北锚碇沉井前两次下沉几何姿态监测及计算简介 [J], 周海生;周亚军
4.HP^3中共形极小曲面的几何 [J], 高紫娟;焦晓祥;
5.基于共形几何代数的空间机构位置解方法研究 [J], 亢书华;李瑞琴
因版权原因，仅展示原文概要，查看原文内容请购买。

中国科学技术大学线性代数课程讲义8

定义 8.1. 设 V 是一个非空集合，V 的元素称为向量，F 是一个数域，F 的元素称为数．在 V 上定义了加法运算 α + β ∈ V 和数乘运算 λ · α ∈ V ，其中 α, β ∈ V ，λ ∈ F．若这两种运算满足下列运算律
(A1) (A2) (A3) (A4) (M1) (M2) (D1) (D2)
比较以上关于 “线性空间” 和 “数域” 的定义，两者有许多相似之处．我们可以认为 “向量” 是 “数” 的推广，向量运算降低了对 “乘法” 的要求，从而有更广泛的应用．
例 8.3. 容易验证，下面是一些常见的线性空间的例子．
1. 在自身的加法和乘法运算下，F 是 F 上的线性空间．
2. 在复数的运算下，复数域 C 是实数域 R 上的线性空间，也是有理数域 Q 上的线性空间．
(A1) (A2) (A3) (A4) (M1) (M2) (M3) (M4) (D1)
加法交换律加法结合律加法有单位元加法有逆元乘法交换律乘法结合律乘法有单位元乘法有逆元加乘分配律
∀a, b ∈ F，a + b = b + a ∀a, b, c ∈ F，(a + b) + c = a + (b + c) 存在零元素 0 ∈ F 使得 0 + a = a，∀a ∈ F ∀a ∈ F，存在 b ∈ F 使得 a + b = 0 ∀a, b ∈ F，a · b = b · a． ∀a, b, c ∈ F，(a · b) · c = a · (b · c) 存在幺元素 1 ∈ F 使得 1 · a = a，∀a ∈ F ∀a ∈ F，若 a ̸= 0，则存在 b ∈ F 使得 a · b = 1 ∀a, b, c ∈ F，(a + b) · c = (a · c) + (b · c)