基于面圖標的三維對稱張量場可視化研究

宋偉杰， 崔俊芝， 葉正麟

（1. 西北工業大學理學院應用數學系，陜西 西安 710072；2. 中國科學院數學與系統科學研究院, 北京 100080）

基于面圖標的三維對稱張量場可視化研究

宋偉杰1， 崔俊芝2， 葉正麟1

（1. 西北工業大學理學院應用數學系，陜西 西安 710072；2. 中國科學院數學與系統科學研究院, 北京 100080）

現有的三維對稱張量場可視化的面圖標方法，僅表現了張量場的某一特征向量場在面圖標上的方向信息，而沒有表現出特征向量場的大小，基于此，提出了在面圖標上通過繪制三維等值線及利用顏色屬性來表現特征向量場的大小信息，這樣，面圖標就將特征向量場的方向和大小信息同時表現出來，從而進一步完善了三維對稱張量場的面圖標方法，更全面地揭示了張量場位于空間曲面上的信息。

工程圖學；可視化；面圖標；等值線；對稱；張量場

張量場可視化是科學計算可視化的一個極其重要的研究領域，這是因為張量場在物理學、力學和醫學等諸多領域中有著非常廣泛的應用。例如，固體力學中的應力和應變，流體力學中的應力、粘性應力、雷諾應力、應變率等等都是二階張量。在醫學領域，由核磁共振圖像重構得到的擴散張量(DT-MRI)也為二階張量。因此，深入研究張量場的可視化對于幫助人們理解復雜的物理和醫學現象以及揭示其中的科學規律具有非常重要的科學意義。同時，與標量場和向量場相比，張量場更為復雜、抽象，且難以利用已有成熟的計算機圖形學技術進行可視化。因此，張量場的可視化一直是科學計算可視化中最具挑戰的研究領域，目前研究剛剛起步。二階張量場中應用最為廣泛的一類是二階對稱張量場，而非對稱張量場往往可以通過分解轉化為對稱張量場，因此二階對稱張量場的可視化是目前張量場可視化的研究熱點。

目前二階對稱張量場的可視化方法可以劃分為兩類。一類是采用標量場的可視化方法來實現張量場的可視化，如描繪出張量場的每個分量場的等高線或等值面。這類方法存在的一個致命的缺點是，可視化的結果嚴重依賴于坐標系的選擇，即隨著坐標軸的旋轉顯示的圖形也隨之發生改變。第二類方法是通過對與之等價的特征向量場進行可視化來實現張量場的可視化。由于特征向量場關于坐標軸的旋轉具有不變性，因此可視化的結果不依賴于坐標系的選擇。這類方法自提出之后受到了人們的普遍關注，出現了許多經典的可視化方法。這些方法包括基本圖標方法[1-7]、基于紋理的方法[8-9]、基于拓撲的方法[10-15]和基于變形的方法[16-17]等。其中，基本圖標方法包括點圖標方法[1-2,6]、線圖標方法[3,5]和面圖標方法[4,7]，這3種方法分別在“點”、“線”和“面”的層面上來表現張量場的信息。相對于點圖標和線圖標，面圖標可以在更大的程度上體現張量場的信息和張量場自身的連續性，因而是目前三維對稱張量場的一種很有前景的可視化方法。Zhang Song[4]和宋偉杰[7]提出的面圖標方法通過構造空間曲面來展現張量場，其中，張量場在空間曲面上的某一特征向量場垂直于該曲面。這樣，面圖標方法通過空間曲面表現了位于曲面上的某一特征向量場的方向（即垂直于曲面）。由于特征向量場不僅有方向，而且有大小，上述面圖標方法僅僅在空間曲面上表現了特征向量場的方向，而沒有表現其大小；另一方面，等值線作為科學計算可視化中一種重要的可視化方法，可以直觀地反映標量場的分布情況及變化趨勢，同時，顏色屬性也是一種非常有效的表征標量場大小的手段。因此，本文提出在原有面圖標即空間曲面上通過繪制等值線和利用顏色屬性來表現特征向量場的大小，這樣，在利用面圖標表現特征向量場的方向的同時，其大小信息也表現了出來。本文的工作進一步完善了對稱張量場的面圖標方法，表現了張量場位于空間曲面上的更多的信息。

1 面圖標上等值線的生成及顏色屬性的設置

面圖標方法[4,7]構造的面圖標即空間曲面均是采用三角網格來近似表示的，因此，在空間曲面上生成等值線及設置顏色屬性的問題就歸結為在空間三角網格上的相應問題。

1.1 空間三角網格上等值線的生成

當三角網格頂點的密度較高，即當每個三角片都很小時，可以假定場量函數沿三角片的邊呈線性變化。基于這樣的假設，在每個三角片內，可以通過對它每條邊的兩端點處的場量值進行線性插值求出等值線與該邊的交點，并將這些交點用直線段連接起來，形成等值線的一部分。遍歷所有的三角片，最終生成整條等值線。下面具體給出等值線的生成過程。

1.1.1 判斷等值線與三角片的邊是否相交

即，當且僅當邊 )(IJ的其中一個頂點的場量值大于等于0f，而另一個頂點的場量值小于等于0f時，場量值為0f的等值線與該邊相交。當公式（1）中的“=”成立時，等值線通過該邊的某個頂點（一個或兩個），為統一處理等值線與邊相交的各種情況和便于編程實現，將該頂點的場量值暫時增加一個微小的增量ε，這樣，等值線通過頂點的情況就可視為與邊（內部）相交的情況來處理，同時這種處理方式對生成的等值線的精度影響不大。由于三角片 )(IJKT= 的每個頂點的場量值現在只有大于0f和小于0f兩種可能的情況，因此，等值線與三角片的位置關系只有圖1所示的八種情況。

圖1 等值線與三角片的八種位置關系。其中，頂點處的“+”和“-”分別表示該頂點的場量值大于f0和小于f0

如果將等值線的值和三個頂點的場量值的大小關系顛倒過來，那么等值線與三角片三個頂點之間的拓撲關系不會改變，生成的等值線是相同的，這樣上面的八種情況就可以視為四種情況；再利用三角片的三個頂點存在的旋轉對稱性，可以進一步將四種情況視為兩種情況，即等值線與三角片的兩條邊相交和等值線不與三角片任何一條邊相交。

1.1.2 計算等值線與三角片邊的交點

1.1.3 生成等值線

在場量函數沿三角片的邊呈線性變化的假設下，公式(2)可以精確地計算出等值線在三角形兩條邊上的等值點。對于等值線在兩個等值點之間的部分，如果三角網格的每個三角片都很小，則可以近似地用直線段連接起來。遍歷三角網格的每一個三角片，求出所有的等值線段，這樣就生成了整條等值線。

1.2 空間三角網格上顏色屬性的設置

等值線可以反映出場量值的分布情況，且等值線之間的疏密程度還反映出場量變化的梯度。如果輔以顏色屬性，則可以連續的表現出場量值的大小及分布規律。本文采用通常的藍色紅色顏色變化，藍色表示較小的場量值，紅色表示較大的場量值。將三角網格的每個頂點根據其場量值的大小賦予相應的顏色，每個三角片內各處的顏色則通過對其頂點處已設置的顏色作線性插值得到。

2 實例結果

本例中的三維對稱張量場為兩豎直向下的外力作用在一正立方體的上表面所形成的應力場。其中，立方體棱長為30mm，材料為鋼，彈性模量E=1.93e11Pa，泊松比V=0.29，邊界條件為下表面固定，承受載荷為兩外力豎直向下作用于上表面中心的兩點上，均為600N。

本例采用宋偉杰[7]的面圖標方法構造出上述應力場的最小主應力的一個面圖標（見圖2），然后在該面圖標上繪制最小主應力的等值線（見圖3）和設置顏色屬性（見圖4）。為表現最小主應力的實際大小，這里取最小主應力的模的絕對值。

圖2 應力場的最小主應力的一個面圖標

圖3 在面圖標上繪制場量值分別為1.1、2.1和3.1的三條等值線

圖4 在面圖標上設置顏色屬性

圖5 在面圖標上同時采用等值線和顏色屬性表現最小主應力的大小

由圖3～圖5看出，面圖標作為空間的一個曲面，在表現應力場位于其上的最小主應力的方向（垂直于曲面）的同時，其上的等值線和顏色屬性則可以表現最小主應力在面圖標上的大小信息。

3 結 論

已有的三維對稱張量場的面圖標方法采用空間曲面來表現某一特征向量場，與點圖標和線圖標相比，面圖標可以在更大的程度上揭示張量場的內在結構和連續性，但是，現有面圖標方法僅僅表現出了特征向量場的方向信息，而特征向量場除了具有方向屬性之外，其大小也是一個很重要的信息，因此現有面圖標方法所揭示的特征向量場的信息是不完整的。本文提出在已構造的面圖標上通過繪制等值線以及設置顏色屬性來表現特征向量場的大小，這樣，在用面圖標表現特征向量場的方向的同時，大小信息也表達了出來。本文工作完善了三維對稱張量場的面圖標方法，更為全面的揭示了張量場位于空間曲面上的信息。

盡管本文工作完善了面圖標方法，使得面圖標可以同時表現特征向量場的方向和大小，但是卻只能表現其中的一個特征向量場，而對于三維對稱張量場來說，有三個特征向量場，因此，研究如何利用面圖標同時表現與之相切的另外兩個特征向量場是一件非常有意義的工作，這可能需要采用van Wijk的IBFVS方法[19]，至于如何利用該方法來實現將是作者今后研究工作的一個目標。

[1]Haber R B. Visualization techniques for engineering mechanics [J]. Computing Systems in Engineering,1990, 1(1)：37-50.

[2]Kindlmann G. Superquadric tensor glyphs [C]//Proceedings of IEEE TVCG/EG Symposium on Visualization, 2004：147-154.

[3]Delmarcelle T, Hesselink L. Visualizing second-order tensor fields with hyperstream-lines [J]. IEEE Computer Graphics and Applications, 1993, 13(4)：25-33.

[4]Zhang S, Demiralp C, Laidlaw D. Visualizing diffusion tensor MR images using streamtubes and streamsurfaces [J]. IEEE TVCG, 2003, 9(4)：454-462.

[5]陳 莉, 竹島由里子, 藤代一成, 等. 大規模數據場的并行可視化[J]. 浙江大學學報(理學版), 2001,28(2)：222-226.

[6]Chen W, Zhang S, Correia S. Visualizing diffusion tensor imaging data with merging ellipsoids [C]//IEEE Pacific Visualization Symposium, 2009：145-151.

[7]Song W J, Cui J Z, Ye Z L. Visualizing 3D symmetric tensor fields using a type of surface icons [C]//IEEE CAD/Graphics’09, 2009：627-631.

[8]Cabral B, Leedom L C. Imaging vector fields using line integral convolution [C]//SIGGRAPH, Computer Graphics Proceedings, 1993：263-270.

[9] Zheng X Q, Pang A. HyperLIC [C]//IEEE Visualization,2003：249-256.

[10]Delmarcelle T. The visualization of second-order tensor fields [D]. Stanford University, 1995.

[11]Zheng X Q, Pang A. Topological lines in 3D tensor fields [C]//IEEE Visualization, 2004：313-320.

[12]Zheng X Q, Parlett B, Pang A. Topological lines in 3D tensor fields and discriminant hessian factorization [J]. IEEE Transaction on Visualization and Computer Graphics, 2004, 11(4)：395-407.

[13]Zheng X Q, Tricoche X, Pang A. Degenerate 3D tensors [C]//Visualization and Processing of Tensor Fields, 2006：247-256.

[14]Zheng X Q, Parlett B, Pang A. Topological structures of 3D tensor fields [C]//IEEE Visualization, 2005：551-558.

[15]Schultz T, Theisel H, Seidel H P, Topological visualization of brain diffusion MRI data [J]. IEEE Transactions on Visualization and Computer Graphics, 2007, 13(6)：1496-1503.

[16]Zheng X Q, Pang A. Volume deformation for tensor visualization [C]//IEEE Visualization, 2002：379-386.

[17]Zheng X Q, Pang A. Interaction of light and tensor fields [C]//Proceedings of The Symposium on Data Visualization, 2003：157-166.

[18]Lorenson W E, Cline H E. Marching cubes：a high resolution 3D surface construction algorithm [J].Computer Graphics, 1987, 21(3)：163-169.

[19]Jarke J, van Wijk. Image based flow visualization for curved surfaces [C]//Proceedings of the IEEE Visualization Conference, 2003：123-130.

Research on the Visualization of 3D Symmetric Tensor Fields Based on Surface Icons

SONG Wei-jie1, CUI Jun-zhi2, YE Zheng-lin1
( 1. Department of Applied Mathematics , School of Science, Northwestern Polytechnical University, Xi’an Shaanxi 710072, China;2. Academy of Mathematics and System Sciences, Chinese Academy of Sciences, Beijing 100080, China )

The current surface icon method to visualize a 3D symmetric tensor field just displays the direction of certain eigenvector field of the tensor field, but not its magnitude. The paper presents that rendering 3D isolines and setting color on the surface icons can display the magnitude of the eigenvector field, thus the direction and magnitude of the eigenvector field are displayed simultaneously on the surface icons, which improves the surface icon method of 3D symmetric tensor fields and reveals more completely the information of tensor fields on surfaces.

engineering graphics; visualization; surface icon; isoline; symmetric; tensor field

TP 391

A

1003-0158(2010)06-0146-05

2010-07-13

國家自然科學基金資助項目（10802068）；西北工業大學基礎研究基金資助項目（JC200949）

宋偉杰（1977-），男，山西長治人，講師，在讀博士，主要研究方向為科學與工程計算、計算機圖形學。

book=150,ebook=151