三維服裝仿真中人體左右剖面的構造

張妙維， 修 毅

(北京服裝學院 計算機信息中心， 北京 100029)

在三維服裝仿真人體變形過程中，人體軀干部分形態復雜，變化豐富，研究方法多樣。隨變形技術的更新和研究的深入，變形效果和交互效率得到提高。

目前，研究較多的一類變形方法是通過掃描大量的人體信息建立人體模型庫[1]，將新導入的人體關鍵信息與模型庫中的人體相匹配，將匹配度高的人體經再次合成得到的人體用于新人體[2]。為使精確度和體型覆蓋度得到提升，這類研究方法所需人體模型量大，隨著研究的展開，模型數量持續增多，過程繁瑣。另一類研究較多的是特征驅動變形：根據體型特征將人體分成多個特征區域，針對每個特征區域制定不同的約束和驅動關系[3]；多種變形方法共同驅動復雜的人體變形，增加局部變形調整的靈活性，但變形效果常與實際人體體型規律不符[4]。

孫守遷等[5]提出的參數化人體建模方法中，對腰腹部的變形是基于計算機圖形學軸變形技術，將人體模板剖成左右2個部分，并且使用徑向變形權重曲線來增加調整的靈活性。該方法建模速度快，但使用的3個人體模板均為標準人體模板，且變形之后的人體模板體型仍在標準體型內。耿玉磊等[4]提出了多因素驅動人體模型變形的方法，可對已有人體進行快速變形。其中對腰圍的變型采用的是基于截面環的尺寸驅動算法，這種算法使腰圍前后統一變形，變形之后的人體腰腹部突兀，不符合實際人體體型規律。

為得到符合實際人體體型規律的變形，引入一種左右剖面模型算法，該算法沿著肩線和體側線，橫跨人體左右將人體分為前后2部分剖面，可用于對前后2個部分的人體分別進行變形。在人體解剖學里，將人體分成前后2個部分的平面，稱為冠狀面，與鉛垂軸平行。由于本文引入的這種剖面是與人體相匹配的傾斜曲面，不滿足冠狀面的定義和性質，因此，本文給引入剖面新的定義，根據其橫貫人體左右的性質，稱其為左右剖面。獲得三維服裝仿真人體左右剖面的主要流程包括：掃描并處理人體；獲取仿真人體用于構造左右剖面的特征點；確定用于構造左右剖面的二次曲面方程；利用獲取的特征點的坐標數據和確定的二次曲面方程，構成非齊次線性方程組求解，獲得二次曲面系數；求取左右剖面。

1 左右剖面特征點獲取

三維人體模型是數字化三維服裝造型與設計技術的基礎，模擬出的人體要包含真實人體信息，對其進一步處理時獲得的數據才有意義。本文采用三維激光儀，通過非接觸人體掃描技術掃描人體，獲得stl.格式的點云數據人體模型[6]。

通過掃描獲取的人體模型，原始數據處理量大，有干擾信息，因此，需要使用Geomagic Studio逆向建模工程軟件對掃描人體數據進行拼合、簡化、過濾、三角化等，進行反復修剪、刪除、排除體外孤點、降噪、修補以及平滑處理得到可用于仿真的人體模型。

由于obj.文件格式簡單、通用，對這種以純文本形式存儲了頂點、法線、紋理坐標和材質信息的模型使用起來方便。通過Geomagic Studio將處理過的掃描得到的stl.格式人體模型數據以obj.格式文件導出。在Visual Studio中，利用Csharp OpenGL圖形處理接口，搭建para 3D model仿真系統，導入用于仿真的obj.格式的人體數據，得到仿真人體如圖1所示。

圖1 仿真人體Fig.1 Simulation body

在該仿真系統中，利用射線求交法手動標識仿真人體上的點，選擇肩線和體測線上明顯、固定、易測部位的一組點為特征點(包括側頸點，肩端點，前腋窩點，胸圍、腰圍、臀圍分別與體側線的交點)進行標識，并獲取特征點三維坐標數據，特征點選取情況如圖2所示。

圖2 特征點的標注Fig.2 Mark of feature points

2 二次曲面方程確定

在獲取了左右剖面特征點的三維坐標數據x,y,z后，需對構成左右剖面的曲面方程類型進行確定。由于獲取的特征點三維坐標數據是確定的測量值，構成左右剖面的曲面曲率變化平緩，引入二次曲面構成人體左右剖面的模型為：

Ax2+By2+Cz2+Dxy+Exz+Fyz+

Gx+Hy+Iz+J=0

(1)

式中，A～J為從人體左右剖面二次曲面議程的各項系數。

從解析二次曲面方程[7]分析，如果將1組已知三維坐標的點代入式(1)中，將得到一個多元一次齊次方程組。該方程組只有零解，無法求取二次曲面方程系數。

從人體結構上分析，人體軀干傾斜度整體趨于平緩，仿真人體左右剖面上的三維坐標x、y、z中，z坐標對剖面影響最小，所以去掉z的二階項，只保留一階項，經簡化后得到二次曲面方程[8]為：

Ax2+By2+Cxy+Dx+Ey+F=z

(2)

將該方程作為左右剖面的模型。

3 非齊次線性方程組求解

使用特征點的三維坐標數據和確定的二階曲面方程，構成方程組求解方程系數[9]。將標識的特征點(xi,yi,zi)(i=1,2,…,n)三維坐標數據代入式(2)中，由于x、y、z值已知，多元二次曲面方程組構成多元一次方程組，如式(3)所示：

(3)

為求解特征點列和方程組系數[10]，將式(3)(n為特征點個數)轉換為矩陣PX=S進行求解，將每行的已知x、y、z值轉化為系數矩陣P,S為右邊值矩陣,X為反求系數的值，具體表示式為：

(4)

代入不同數量的特征點會導致系數值與方程個數出現以下幾種情況：

1)n=6,此時為定解方程組，利用矩陣除法X=S/P可得方程組有唯一解；

2)n<6,此時為奇異方程組，方程組有無數解；

3)n>6,此時為超定方程組，可用最小二乘法求得近似解。

4 左右剖面求取

在標識特征點的過程中，獲取的三維坐標值成對出現，由于人體左右面對稱使得左右對稱特征點的x值互為相反數，y、z相等。將左右對稱的特征點代入二次曲面方程，并組成方程組時，所有包含x奇次冪的項都將被抵消。此時式(3)的系數發生變化，求解的不確定因素增加，因此，選擇使用仿真人體一側的特征點，即頸圍線的中心點和腰圍線的中心點作為新的特征點進行標識，并用這些點的數據組成新的線性方程組求解。在仿真人體上重新取點結果如圖3所示。

圖3 特征點的重新標注Fig.3 Re-mark of feature points.

選擇跨平臺、跨編程語言支持的OpenGL作為3-D繪圖的API進行3-D仿真，使用面向對象語言的C#語言在Visual Studio、Para 3D Model仿真系統中添加左右剖面的功能。

對指定的特征點進行標識，將指定特征點的三維坐標值中包含x、y項的值存入系數矩陣P中；z的值存放在等號右側P中；調用數學庫中的QR矩陣方法求解PX=S；求出系數A～F的值。

根據非齊次線性方程組的求解原理，求解二階曲面方程組成的方程組時，當n≥6時方程組有可用解。求解定解方程組、超定方程組獲得的二階曲面方程，呈現的是效果不同的曲面，因此，把求取左右剖面分成標識的特征點數等于6和大于6的2種情況進行分析。

對包含左右剖面的仿真人體進行前、后、左、右觀察，當特征點數為6時，僅取最關鍵的特征點(側頸點，肩端點，胸圍、腰圍與體側線的交點，頸圍線、腰圍線的中心點)，剖面精確地經過這些特征點，效果良好，其結果如圖4所示。

圖4 標識點為6的左右剖面Fig.4 Six marked feature points of left-right section. (a) Left view; (b) Right view

當特征點個數為9時，此時標識的點包含人體關鍵部位的所有特征點(側頸點，肩端點，前腋窩點，胸圍、胸下圍、腰圍、臀圍與體側線的交點，頸圍線、臀圍線的中心點)，其結果如圖5所示。

圖5 標識點數為9的左右剖面Fig.5 Nine marked feature points of left-right section. (a)Left view; (b)Right view

圖6 局部細節對比Fig.6 Local detail contrasts.(a) Six points; (b) Nine points

標識不同數量特征點生成的剖面局部細節對比如圖6所示。

由圖6可知：當標識點個數為6時，剖面全部經特征點；標識點個數為9時剖面大部分經特征點，部分靠近特征點，滿足將人體分成前后2個部分的要求。

除對能表征人體測量基準點和測量部位的關鍵點的標識，即對具有典型代表性的6個和9個特征點的標識外，選取特征點個數為8、10、12、15進行標識，結果如圖7所示。在圖7(a)～(d)所示4幅剖面展示圖中，隨標識點的增多非特征點標識的個數也隨之增加，剖面對肩線和體側線上一般點的標識也具有適用性，可滿足將仿真人體很好地分成前后2個部分的要求。

圖7 標識實例Fig.7 Mark of examples. (a) Eight points; (b) Ten points ; (c) Twelve points; (d) Fifteen points

5 結束語

在構造空間曲面時，經常采用插值算法和擬合算法，但需要具有一組分布在曲面空間上的數據點作為這些算法的輸入數據。在構造人體左右剖面時，由于通過三維掃描只能獲取仿真人體表面點的數據，僅已知仿真人體表面肩線和體測線上相關的特征點，無法標識體內點并獲取點的三維坐標信息進行利用，因此，本文引入空間二次曲面方程作為左右剖面模型，再利用非齊次方程組求解，最終得到了二次曲面方程，構造出左右剖面。

針對不同時期人體胖瘦脂肪分布規律，將標準人體變形為特殊體型時，人體前后需要進行不同規律變形的目的，生成人體左右變形的基準剖面，通過對肩線和體側線上點的標識以及形成非齊次線性方程組來求取曲面，該過程簡捷、高效。

監測輸液正常的關鍵指標之一: 輸液液滴的狀態,液滴的快慢代表著輸液的不同狀態。圖2為液滴信號監測示意圖。紅外對管是檢測輸液液滴的關鍵器件。圖3是實現輸液液滴監測的具體實現電路。

在后續研究中發現，應用左右剖面技術不僅對仿真人體前后區別變形具有良好可用性，還適用于作為服裝造型前后部分的基準，對前后相連的統一圍線等的分別獲取，三維服裝到二維平面展平中的中間媒介，三維仿真服裝通過軸旋轉方式映射到左右剖面上等研究與應用。這些應用為三維服裝CAD仿真人體的變形，以及三維服裝到二維平面的展平的研究提供了新的思路。

FZXB

參考文獻：

[1] WANG C C L.Parameterization and parametric design of mannequins [J]. Computer-Aided Design, 2005, 37(1)：83-98.

[2] CHU C H, TSAI Y T, WANG C C L, et al. Exemplar-based statistical model for semantic parametric design of human body[J]. Computers in Industry, 2010, 61(6):541-549.

[3] 路寧,寧濤,唐榮錫,等.自由變形方法及其在建立尺寸驅動人體模型中的應用[J]. 工程圖形學報，2004(1)：76-83.

LU Ning, NING Tao, TANG Rongxi, et al. Dirichlet free-form deformations and their applicationtion dimension driven body model[J]. Journal of Engineering Graphics, 2004(1): 76-83.

[4] 耿玉磊,陸國棟,鄧衛燕,等. 面向服裝CAD的多因素驅動人體模型變形技術研究[J]. 中國圖象圖形學報，2009，14(11):2349-2355.

GENG Yulei, LU Guodong, DENG Weiyan, et al.

[5] 孫守遷,徐愛國. 三維服裝仿真中的參數化人體建模技術[J]. 紡織學報,2007,28(12):89-93.

SUN Shouqian, XU Aiguo.Study on parameterized mannequin modeling in 3-D garment simulation[J]. Journal of Textile Research, 2007, 28(12): 89-93.

[6] 陳亮,趙曙光,張麗娟,等. 虛擬三維服裝展示的發展歷史與研究熱點[J]. 紡織學報,2011,32(10):153-160.

CHEN Liang, ZHAO Shuguang, ZHANG Lijuan, et al. History and Research hot spot of three-dimensional virtual garment display[J]. Journal of Textile Research, 2011, 32(10): 153-160.

[7] 李二濤,張國煊,曾虹. 基于最小二乘的曲面擬合算法研究[J]. 杭州電子科技大學學報,2009,29(2):48-51.

LI Ertao, ZHANG Guoxuan, ZENG Hong.Algorithm of surface fitting research based on least-squares me-thods[J]. Journal of Hangzhou Dianzi University, 2009, 29(2):48-51.

[8] 霍曉程,李小平. 用最小二乘法擬合曲面方程[J]. 赤峰學院學報(自然科學版),2009,25(6):11-13.

HUO Xiaocheng, LI Xiaoping.Using least squares fitting surface equation[J].Journal of Chifeng Univer-sity (Natural Science Edition),2009, 25(6):11-13.

[9] 張文信,黃陽,潘衛國. 反求工程的一種曲面插值算法[J]. 機電產品開發與創新,2008,21(5):25-27.

ZHANG Wenxin, HUANG Yang, PAN Weiguo. Palletizing robot control system design and implementation[J].Development & Innovation of Machinery & Electrical Products,2008, 21(5):25-27.

[10] 戴中林. 多元線性不定方程的矩陣解法[J]. 四川師范學院學報(自然科學版),2001,22(3):270-273.

DAI Zhonglin.Solution of the multivariate linear indeterminate equation by a matrix[J].Journal of Sichuan Teachers College (Natural Science Edition),2001,22(3):270-273.