上次課程，我們講解了拓?fù)鋱A盤間的調(diào)和映照。這次，我們討論拓?fù)淝蛎骈g的調(diào)和映照。

圖1. 共形腦圖。

拓?fù)淝蛎骈g的調(diào)和映照在醫(yī)學(xué)圖像領(lǐng)域被頻繁使用，特別是構(gòu)造大腦皮層曲面到單位球面間的映射，我們稱之為“共形腦圖”，如圖1所示。

圖2. 拓?fù)淝蛎骈g的調(diào)和映射。

外蘊(yùn)方法 我們將目標(biāo)曲面單位球面嵌在三維歐式空間中，這樣曲面間的映射可以被表示成從源曲面到三維歐式空間的映射，并且像集被限制在單位球面上，

進(jìn)一步，這一映射由三個(gè)坐標(biāo)函數(shù)來表示，

映射的拉普拉斯由坐標(biāo)函數(shù)的拉普拉斯給出，

傳統(tǒng)的熱流方法將一個(gè)函數(shù)經(jīng)過“熱力學(xué)擴(kuò)散”成一個(gè)調(diào)和函數(shù)，使得其調(diào)和能量隨時(shí)間單調(diào)下降：

但是，在我們目前的情況下，映射的像被限制在單位球面上，因此每一點(diǎn)的像只能沿著單位球面的切方向移動(dòng)，而無法沿著球面的法方向移動(dòng)。因此，我們需要將傳統(tǒng)的熱流方法修改成非線性熱流：

這里等式右側(cè)是映射拉普拉斯的切向分量，

這樣我們可以保證對(duì)任意時(shí)間和任意點(diǎn)，映射像一直在單位球面上，

由此可見，所謂的“非線性”是指在每一點(diǎn)，每一步我們都需要向切空間投影，因此能量泛函的歐拉-拉格朗日方程不是線性偏微分方程，其離散化的近似不是線性方程組：

非線性熱流方法可以確保映射的調(diào)和能量單調(diào)遞降，但是由于調(diào)和映射不唯一，計(jì)算并不穩(wěn)定，我們需要添加更多的條件來確保解的唯一性。事實(shí)上，如果存在兩個(gè)映射都是調(diào)和映射，則它們彼此相差一個(gè)球面到自身的共形變換。首先，我們用球極投影將球面映到擴(kuò)展復(fù)平面上，

圖4. 球極投影。

如圖4所示，我們?cè)趩挝磺虻谋睒O放置一盞燈，從北極發(fā)射的光線將球面上所有點(diǎn)映到過南極的切平面上，同時(shí)將北極點(diǎn)映成無窮遠(yuǎn)點(diǎn)。直接計(jì)算表明，球極投影是共形映射。擴(kuò)展復(fù)平面到自身的所有共形變換構(gòu)成所謂的莫比烏斯變換群，每一個(gè)自映射都具有如下形式：

如果存在兩個(gè)調(diào)和映射，則它們相差一個(gè)莫比烏斯變換：

為了去掉莫比烏斯變換群帶來的歧義性，我們添加一些歸一化條件，例如我們要求映射滿足

這個(gè)限制去掉了3個(gè)自由度，余下的莫比烏斯變換是單位球面到自身的旋轉(zhuǎn)。非線性熱流方法是穩(wěn)定的，不會(huì)誘發(fā)旋轉(zhuǎn)，因此這一歸一化條件足夠確保調(diào)和映射解的唯一性。

非線性熱流方法將一個(gè)初始映射“擴(kuò)散”成一個(gè)調(diào)和映射，那么我們?nèi)绾芜x取初始映射？理論上，任何一個(gè)映射度為1的映射都可以，例如最為常見的高斯映射。但在實(shí)踐中，光滑曲面被離散的多面體網(wǎng)格所逼近，優(yōu)化過程有可能落入局部最優(yōu)的陷阱。一種有效的避免局部最優(yōu)陷阱的方法如下：我們將源曲面一分為二，將每一半用拓?fù)鋱A盤的調(diào)和映射映到單位圓盤，并且邊界映射相互一致。然后我們用球極投影將兩個(gè)單位圓盤映到上半球面和下半球面，如此得到曲面到單位球面的初始映射。然后，我們?cè)儆梅蔷€性熱流方法將初始映射進(jìn)一步擴(kuò)散。

內(nèi)蘊(yùn)法 內(nèi)蘊(yùn)方法只需要曲面的黎曼度量，不需要曲面在歐式空間中的等距嵌入。外蘊(yùn)法引導(dǎo)我們得到計(jì)算方法，內(nèi)蘊(yùn)法使我們能夠洞察調(diào)和映射更為深刻的性質(zhì)。假設(shè)曲面間的映射給定，我們?nèi)∏娴牡葴刈鴺?biāo)：

我們記

映射的調(diào)和能量密度為

映射的調(diào)和能量為

由此，我們可以看到調(diào)和能量只和源曲面的共形結(jié)構(gòu)有關(guān)，和具體的共形黎曼度量無關(guān)。

進(jìn)一步，我們得到調(diào)和能量的歐拉-拉格朗日為

從而，內(nèi)蘊(yùn)的非線性熱流方程為

圖5. 調(diào)和映射

圖6. 共形映射

通常情況下，共形映射一定是調(diào)和的，調(diào)和映射不一定是共形的。在上圖情形，當(dāng)我們固定邊界映射，優(yōu)化調(diào)和能量，我們得到調(diào)和映射。如果，我們放開邊界，令邊界的像在單位圓上可以自由滑動(dòng)，從而進(jìn)一步減小調(diào)和能量，則我們得到共形映射。換言之，共形映射是所有調(diào)和映射中調(diào)和能量最小者。

但是，對(duì)于球面而言，拓?fù)淝蛎骈g的調(diào)和映射一定是共形的。直觀來說，拓?fù)淝蛎骈g的映射可以在球面上自由滑動(dòng)，從而使調(diào)和能量無障礙地達(dá)到最優(yōu)，從而得到共形映射。下面，我們將這一直覺述諸嚴(yán)格的證明。我們定義曲面間映射所誘導(dǎo)的二次微分，霍普夫微分：

我們可以證明，如果霍普夫微分全純，則映射必為調(diào)和；如果霍普夫微分為0，則映射必為共形。

曲面上所有的全純二次微分構(gòu)成一個(gè)群，根據(jù)黎曼-羅赫定理，虧格為g>1的曲面，這個(gè)群的維數(shù)為6g-6。黎曼-羅赫定理是指標(biāo)定理的特殊形式，本質(zhì)上是說流形上橢圓型偏微分方程解的空間維數(shù)被流形拓?fù)錄Q定。

拓?fù)淝蛎嫔纤械娜兌挝⒎直厝粸?。假如球面上存在一個(gè)全純二次微分，那么它誘導(dǎo)一個(gè)帶有奇異點(diǎn)的平直度量，奇異點(diǎn)對(duì)應(yīng)著全純二次微分的零點(diǎn)，因此奇異點(diǎn)處的曲率測(cè)度為-180度。根據(jù)高斯博內(nèi)定理，總曲率等于球面的歐拉示性數(shù)乘以360度，從而奇異點(diǎn)的個(gè)數(shù)為0，全純二次微分沒有零點(diǎn)。我們能夠取其平方根的一個(gè)分支，得到一個(gè)全純1-形式。球面上的全純1-形式的實(shí)部為調(diào)和1-形式。因?yàn)榍蛎娴囊浑A上同調(diào)群為0，所以調(diào)和1-形式必為0，進(jìn)而全純二次微分必為0。

因此，拓?fù)淝蛎骈g的所有的調(diào)和映射，其霍普夫微分為全純二次微分，必然為0。我們得到拓?fù)淝蛎骈g的所有調(diào)和映射必為共形映射。

我們直接計(jì)算調(diào)和能量：

這里雅克比矩陣

因此我們得到

這里

是映射誘導(dǎo)的Beltrami系數(shù)， 

是映射的伸縮商。

圖2. Beltrami系數(shù)和伸縮商的幾何意義。

我們考察從人臉曲面到平面單位圓盤的映射，如圖2所示，映射將曲面上的無窮小橢圓映到平面的無窮小圓。無窮小橢圓的偏心率，就是伸縮商K，偏心率和橢圓長(zhǎng)軸方向構(gòu)成了Beltrami系數(shù)。從上面的推導(dǎo)我們可以看出，調(diào)和能量的下界等于目標(biāo)曲面的面積；當(dāng)所有的伸縮商K都是1時(shí)，調(diào)和能量的下界能夠被達(dá)到，這時(shí)，映射為共形映射。

換言之，在所有可能的調(diào)和映射中，如果存在共形映射，則共形映射的調(diào)和能量達(dá)到全局最優(yōu)。

當(dāng)映射為調(diào)和映射時(shí)，其像為閉測(cè)地線。如果目標(biāo)曲面的曲率有正有負(fù)，那么在同一同倫類中，測(cè)地線并不唯一。如果目標(biāo)曲面的曲率處處為負(fù)，那么同倫的測(cè)地線一定重合。直觀上講，假如曲面上有兩條彼此同倫的閉測(cè)地線，它們之間的同倫在曲面上掃過一個(gè)拓?fù)鋱A柱面。根據(jù)高斯-博內(nèi)定理，我們有

等式左邊，因?yàn)閳A柱面的邊界為測(cè)地線，所有第二項(xiàng)為0，第一項(xiàng)非正。等式右邊為0，所以第一項(xiàng)為0，圓柱面面積為0。兩條測(cè)地線彼此重合。

封閉曲面間的調(diào)和映射，如果目標(biāo)曲面上的高斯曲率處處為負(fù)，且M在N上的像不是一條閉測(cè)地線，則同倫的調(diào)和映射必然重合。一種想法是基于調(diào)和能量的凸性，設(shè)定二階光滑同倫連接著兩個(gè)調(diào)和映射，

我們計(jì)算調(diào)和能量的二階導(dǎo)數(shù)

這里切矢量場(chǎng)定義為

拉回聯(lián)絡(luò)算子，源曲面上的標(biāo)準(zhǔn)正交基底和高斯曲率分別為

我們選取測(cè)地同倫，亦即

為測(cè)地線，那么最后一項(xiàng)恒為0，因?yàn)?img doc360img-src='http://image108.360doc.com/DownloadImg/2017/08/0605/107624720_37' data-ratio='0.32432432432432434' data-type='png' data-w='74' src='http://pubimage.360doc.com/wz/default.gif'>

因此調(diào)和能量的二階導(dǎo)數(shù)恒正，調(diào)和能量為嚴(yán)格凸。 同時(shí)在時(shí)間為0和1點(diǎn)處，映射為調(diào)和映射，調(diào)和能量關(guān)于時(shí)間的一階導(dǎo)數(shù)為0，由此我們有并進(jìn)一步我們得到V必然處處為0，因此起始和終點(diǎn)處的調(diào)和映射重合。

我們首先定義兩個(gè)輔助函數(shù)

那么，映射的雅克比行列式為

對(duì)于調(diào)和映射，我們得到Bochner公式

調(diào)和映射誘導(dǎo)的輔助函數(shù)的零點(diǎn)是孤立零點(diǎn)，我們可以定義零點(diǎn)的階數(shù)。

如果輔助函數(shù)不恒為0，那么輔助函數(shù)零點(diǎn)的總階數(shù)和曲面的拓?fù)浼坝成涞挠成涠戎g有著整體的關(guān)系

由此，我們能夠在特定條件下，討論調(diào)和映射的存在性。例如，不存在從拓?fù)洵h(huán)面到拓?fù)淝蛎?，度為一的調(diào)和映射。

如果調(diào)和映射存在，并且非共形，那么零點(diǎn)的總階數(shù)應(yīng)該非負(fù)，

矛盾。因此調(diào)和映射必然為共形映射，進(jìn)一步映射必為分支覆蓋映射。因?yàn)橛成涞亩葹?，所以映射為同胚。但是我們知道拓?fù)洵h(huán)面和拓?fù)淝蛎嬷g不存在同胚，矛盾。由此，拓?fù)洵h(huán)面和拓?fù)淝蛎嬷g不存在度為1的調(diào)和映射。

我們用反證法給出簡(jiǎn)略證明。假設(shè)在某一內(nèi)點(diǎn)，雅克比行列式為負(fù)，因此區(qū)域

非空。研究函數(shù)

則它在D的邊界上為0， 在D的內(nèi)部處處為負(fù)，因此函數(shù)

為上調(diào)和函數(shù)(super harmonic), 在D的內(nèi)部處處為正，那么雅克比行列式

在D的內(nèi)部處處為正，這和D的定義相矛盾。

References

1. Xianfeng Gu, Yalin Wang, Tony F. Chan, Paul M. Thompson and Shing-Tung Yau. Genus Zero Surface Conformal Mapping and Its Application to Brain Surface Mapping. IEEE Transaction on Medical Imaging (TMI), 23(8):949-958, August 2004.

2. Craig Gotsman, Xianfeng Gu and Alla Sheffer. Fundamentals of Spherical Parameterization for 3D Meshes. ACM Transaction on Graphics (TOG), 22(3):358-363,2003.