基于多融合態(tài)的低密度三維模型信息隱藏算法

任帥 徐振超 王震 賀媛 張弢 蘇東旭 慕德俊

摘 要：針對(duì)現(xiàn)有三維模型信息隱藏算法無法有效抵抗不均勻壓縮的問題，提出一種基于多融合態(tài)的多載體低密度的信息隱藏算法。首先通過平移和縮放對(duì)多個(gè)三維模型進(jìn)行定位、定向及定型;其次對(duì)三維模型進(jìn)行不同角度的旋轉(zhuǎn)，以中心點(diǎn)作為融合點(diǎn)進(jìn)行融合，得到多個(gè)融合態(tài);再次，利用局部高度和Mean Shift聚類分析算法對(duì)融合態(tài)模型的頂點(diǎn)進(jìn)行能量劃分，得到不同能量的頂點(diǎn);最后，通過修改頂點(diǎn)坐標(biāo)的方法將經(jīng)過Arnold置亂變化的秘密信息快速隱藏于多個(gè)融合態(tài)和三維模型中。實(shí)驗(yàn)結(jié)果表明，該算法對(duì)抵御不均勻壓縮的攻擊有很好的魯棒性且具有很高的不可見性。

關(guān)鍵詞：信息隱藏;三維模型;融合態(tài);局部高度理論;Mean Shift聚類分析

中圖分類號(hào)：TP309.2;TP301.6

文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：1001-9081（2019）04-1100-06

Abstract： Aiming at the problem that the existing 3D model information hiding algorithms cannot effectively resist uneven compression， a multi-carrier low-density information hiding algorithm based on multiple fusion states was proposed. Firstly， multiple 3D models were positioned， oriented and stereotyped by translation and scaling. Secondly， the 3D models were rotated at different angles and merged by using the center point as merging point to obtain multiple fusion states. Thirdly， local height and Mean Shift clustering analysis were used to divide the energy of the vertices of the fusion state model， obtaining the vertices with different energies. Finally， by changing the vertex coordinates， the secret information changed by Arnold scrambling was quickly hidden in multiple fusion states and 3D models. Experimental results show that the proposed algorithm is robust against uneven compression attacks and has high invisibility.

Key words： information hiding; 3D model; fusion state; local height theory; Mean Shift clustering analysis

0?引言

三維模型的信息隱藏技術(shù)是當(dāng)今信息隱藏安全領(lǐng)域中的一個(gè)研究熱點(diǎn)。目前，基于三維模型的信息隱藏技術(shù)的研究主要是基于空間域和變換域的。文獻(xiàn)[1]首先以協(xié)方差分析得到三維模型的初始特征點(diǎn)為核心，在它的K近鄰域中構(gòu)建不跨越區(qū)域最小三角形為嵌入單元的底面;其次將剩余頂點(diǎn)按照升序排序，尋找合適的嵌入頂點(diǎn);最后通過改變頂點(diǎn)信息來嵌入水印。文獻(xiàn)[2]通過計(jì)算模型的重心，利用頂點(diǎn)和重心連線與模型求交的方法計(jì)算出重心交點(diǎn)距離比，并通過修改三角形面片3個(gè)頂點(diǎn)對(duì)應(yīng)的重心點(diǎn)距離比（Ratio of Barycenter and Crosspoint， RBC）來嵌入信息。文獻(xiàn)[3]通過修改某些頂點(diǎn)的積分不變量來進(jìn)行秘密信息的嵌入，在一定程度上能抵抗網(wǎng)格編輯和裁剪等攻擊。文獻(xiàn)[4]通過計(jì)算頂點(diǎn)的局部高度，利用Mean Shift聚類分析把模型的頂點(diǎn)分成了魯棒點(diǎn)、亞魯棒點(diǎn)、脆弱點(diǎn)三類，通過實(shí)現(xiàn)不同信息的嵌入提高了隱藏的魯棒性。

文獻(xiàn)[5]利用3D模型自身的結(jié)構(gòu)，選擇和提取幾何特征來進(jìn)行信息隱藏。文獻(xiàn)[6]通過對(duì)特征點(diǎn)和非特征點(diǎn)構(gòu)成的二維矩陣進(jìn)行小波變換來嵌入信息。文獻(xiàn)[7]提出了一種新的三維多邊形網(wǎng)格自適應(yīng)隱寫算法，該算法具有自適應(yīng)性、高效性等特點(diǎn)。由于三維模型在變換域進(jìn)行信息隱藏的難度較大，三維模型信息隱藏算法大部分都是在三維網(wǎng)格的情況下進(jìn)行秘密信息隱寫，可提取的特征值不是很多，對(duì)抵御不均勻壓縮的攻擊魯棒性也不是很好。

本文在上述研究成果的基礎(chǔ)上，提出一種基于多融合態(tài)的多載體低密度信息隱藏算法。該算法通過對(duì)多個(gè)三維模型進(jìn)行定位、定向及定型，以三維模型的中心為融合點(diǎn)，通過沿X軸、Y軸旋轉(zhuǎn)不同角度來得到多個(gè)融合態(tài);然后，利用局部高度和Mean Shift聚類分析算法對(duì)融合態(tài)模型的頂點(diǎn)進(jìn)行能量劃分得出魯棒點(diǎn)、亞魯棒點(diǎn)和脆弱點(diǎn);最后，通過修改頂點(diǎn)坐標(biāo)的值來完成信息隱藏。

1?本文的信息隱藏算法

本文的信息隱藏算法主要分為三個(gè)階段：首先，對(duì)三維模型進(jìn)行預(yù)處理，通過一定的方法構(gòu)建多個(gè)融合態(tài);其次利用局部高度和Mean Shift聚類分析算法對(duì)融合態(tài)模型的頂點(diǎn)進(jìn)行能量劃分，標(biāo)記出融合態(tài)的魯棒點(diǎn)、亞魯棒點(diǎn)和脆弱點(diǎn);最后，通過修改頂點(diǎn)坐標(biāo)的方法對(duì)經(jīng)過Arnold置亂變化的秘密信息進(jìn)行快速信息隱藏。

1.1?以三維模型中心點(diǎn)融合的多融合態(tài)算法

融合態(tài)載體是對(duì)多個(gè)單載體進(jìn)行有效標(biāo)識(shí)和預(yù)處理，然后通過一定的方法組合在一起形成一個(gè)融合態(tài)載體。用多個(gè)載體進(jìn)行信息隱藏，可以降低信息隱藏的密度，形成一種大容量、低密度的隱藏環(huán)境。圖1表示的就是一個(gè)簡單的融合態(tài)模型。

對(duì)三維模型預(yù)處理后，得到三維模型的中心點(diǎn)Cmid、方向向量u和中間頂點(diǎn)到中心的平均距離Lmid后，可以進(jìn)行融合態(tài)的構(gòu)建， 本文算法通過以三維模型的中心Cmid為匯聚點(diǎn)，對(duì)三維模型進(jìn)行平移、旋轉(zhuǎn)，得到三個(gè)不同的融合態(tài)S1、S2、S3。

融合態(tài)S1：把三維模型M′1，M′2，…，M′n依次以同方向向量u，中心坐標(biāo)Cmid，平移到三維坐標(biāo)原點(diǎn)O，得到融合態(tài)S1，結(jié)果如圖3所示。

融合態(tài)S2：在計(jì)算其他融合態(tài)時(shí)，為了在旋轉(zhuǎn)的過程中使每個(gè)三維模型露出來的定點(diǎn)數(shù)量基本均勻，本文算法采用90°除以三維模型個(gè)數(shù)n來計(jì)算旋轉(zhuǎn)的角度α（即α=90°/n）。

把三維模型M′1平移到坐標(biāo)的原點(diǎn)，三維模型M′2以中心點(diǎn)融合，使用式（2）計(jì)算出的方向特征向量u對(duì)三維模型沿著X軸正方向旋轉(zhuǎn)α角度，三維模型M′3沿著X軸正方向旋轉(zhuǎn)2α角度，依次進(jìn)行旋轉(zhuǎn)，三維模型M′n旋轉(zhuǎn)（n-1）α角度，得到融合態(tài)S2，融合態(tài)S2如圖4所示。

融合態(tài)S3：計(jì)算模型需要旋轉(zhuǎn)的角度β（β=90°/n），三維模型M′1的中心移到原點(diǎn)O，使用式（2）計(jì)算出的方向特征向量u對(duì)三維模型M′2沿著Y軸正方向旋轉(zhuǎn)β角度，模型M′3繼續(xù)沿著Y軸正方向旋轉(zhuǎn)2β，模型M′n沿Y軸正方向旋轉(zhuǎn)（n-1）β，這樣就得到了融合態(tài)S3，得到的融合態(tài)S3如圖5所示。

本文算法利用多融合態(tài)理論，通過對(duì)多個(gè)三維模型進(jìn)行縮放、平移、旋轉(zhuǎn)得到多個(gè)融合態(tài)，該算法可以使用三維模型盡量多的點(diǎn)。在后期，通過對(duì)每個(gè)融合態(tài)進(jìn)行秘密信息隱藏，對(duì)提取的秘密信息進(jìn)行校驗(yàn)，提高信息隱藏的魯棒性。

1.2?局部高度的頂點(diǎn)標(biāo)注與篩選

局部高度是南京大學(xué)Lin等[9]提出的一種三維模型的顯著性度量方式，可以計(jì)算三維模型的凸起程度。設(shè)三維模型的頂點(diǎn)v的R-鄰居點(diǎn)集合為NR（v），縮寫為NR，則頂點(diǎn)v的局部高度由式（5）計(jì)算得出：

計(jì)算出頂點(diǎn)的局部高度后，利用Mean Shift算法[10]這種非參數(shù)化的概率密度估計(jì)方法對(duì)載體表面的局部高度進(jìn)行聚類分析，按照局部高度值的大小得出載體的特征，通過測量三維模型頂點(diǎn)的凸起程度（如圖6），可以計(jì)算得到特征點(diǎn)（Feature Vertice， FV）、亞特征點(diǎn)（Sub-feature Vertice，SV）和背景點(diǎn)（Background Vertice， BV）這三類特征點(diǎn)，根據(jù)信息隱藏能量分析理論，F(xiàn)V、SV、BV這3類的能量依次降低，因此把其命名為魯棒點(diǎn)、亞魯棒點(diǎn)和脆弱點(diǎn)。

本文算法中對(duì)融合態(tài)S1、S2、S3的頂點(diǎn)進(jìn)行能量信息標(biāo)注，標(biāo)記三維模型的亞魯棒點(diǎn)為Vsen，此為信息隱藏的嵌入點(diǎn)。

1.3?修改頂點(diǎn)坐標(biāo)的快速信息隱藏算法

本文的三維模型信息隱藏算法是通過對(duì)三維模型頂點(diǎn)坐標(biāo)值的個(gè)位數(shù)字進(jìn)行修改，來實(shí)現(xiàn)信息的隱藏。偶數(shù)代表二進(jìn)制數(shù)字0，奇數(shù)代表二進(jìn)制數(shù)字1，通過對(duì)數(shù)字進(jìn)行加1的方式來對(duì)需要修改的點(diǎn)進(jìn)行修改。例如：頂點(diǎn)V1、V2、V3、V4的X軸坐標(biāo)個(gè)位表示的二進(jìn)制信息為0011（如圖7（a）所示），如果要存儲(chǔ)0010的二進(jìn)制信息，則對(duì)頂點(diǎn)的X坐標(biāo)值修改為圖7（b）的結(jié)果，結(jié)果存儲(chǔ)的二進(jìn)制信息為0010。

該算法是對(duì)三維的頂點(diǎn)坐標(biāo)進(jìn)行修改，因此該方法對(duì)秘密信息具有快速存儲(chǔ)和讀取的性能。圖8表示的是在3dMax仿真平臺(tái)中對(duì)一個(gè)茶壺頂點(diǎn)的坐標(biāo)值進(jìn)行修改。

1.4?信息隱藏的詳細(xì)嵌入步驟

基于多融合態(tài)的多載體低密度信息隱藏算法分為以下8個(gè)步驟，本文算法流程如圖9所示。

步驟4?利用式（3）計(jì)算出三維模型中心到頂點(diǎn)的平均距離（L1，L2，…，Ln），然后對(duì)距離進(jìn)行排序，選取中間的平均距離Lmid為三維模型縮放的標(biāo)準(zhǔn);用式（4）計(jì)算出縮放比例，對(duì)其他三維模型進(jìn)行縮放處理。目的是為了保證三維模型在實(shí)現(xiàn)融合態(tài)時(shí)的大小能基本一致，同時(shí)抵抗均勻壓縮。

步驟5?在3dMax2012平臺(tái)中將縮放好的三維模型以中心點(diǎn)為匯聚點(diǎn)依次平移到坐標(biāo)原點(diǎn)，得到融合態(tài)S1。利用1.2節(jié)中描述的局部高度計(jì)算公式和Mean Shift聚類分析算法，標(biāo)記融合態(tài)S1的亞魯棒點(diǎn)為Vsen-S1并獲取其頂點(diǎn)坐標(biāo)VS1，利用1.3節(jié)中的隱藏算法，通過修改頂點(diǎn)中Z軸坐標(biāo)值的個(gè)位來實(shí)現(xiàn)隱藏，奇數(shù)代表1，偶數(shù)代表0，依次對(duì)頂點(diǎn)坐標(biāo)對(duì)比和修改來嵌入比特序列Cre。

步驟6?把融合態(tài)S1的三維模型拆分開，在以中心點(diǎn)為匯聚點(diǎn)的基礎(chǔ)上，使用方向特征向量u把三維模型依次沿著X軸正方向旋轉(zhuǎn)α角度，得到融合態(tài)S2，標(biāo)記三維模型的亞魯棒點(diǎn)為Vsen-S2。獲取融合態(tài)S2的頂點(diǎn)坐標(biāo)VS2，通過修改頂點(diǎn)坐標(biāo)X軸坐標(biāo)值的個(gè)位數(shù)字來對(duì)秘密信息Cre進(jìn)行嵌入。

步驟7?把融合態(tài)S2的三維模型拆分開，讓三維模型依次沿著Y軸正方向旋轉(zhuǎn)β角度，得到融合態(tài)S3，標(biāo)記三維模型的亞魯棒點(diǎn)為Vsen-S3。獲取融合態(tài)S3的頂點(diǎn)坐標(biāo)VS3，通過對(duì)頂點(diǎn)坐標(biāo)Y軸坐標(biāo)值的個(gè)位數(shù)字進(jìn)行修改嵌入秘密信息Cre。

步驟8?把融合態(tài)S3分開后得到含有秘密信息的三維模型M*1，M*2，…，M*n。

1.5?信息提取步驟

秘密信息的提取是信息嵌入的逆過程，主要分為以下幾個(gè)過程：第1步?按照信息隱藏嵌入步驟中的步驟1～5對(duì)含有秘密信息的三維模型M*1，M*2，…，M*n進(jìn)行平移、旋轉(zhuǎn)和縮放處理，使三維模型恢復(fù)到組合融合態(tài)之前的狀態(tài)，然后恢復(fù)3個(gè)融合態(tài)，并從中提取秘密信息。

第2步?按照步驟8恢復(fù)融合態(tài)S3，對(duì)三維頂點(diǎn)的Y軸值進(jìn)行提取，得到比特序列Cre3。

第3步?按照步驟7恢復(fù)融合態(tài)S2，對(duì)三維頂點(diǎn)的X軸值進(jìn)行提取，得到比特序列Cre2。

第4步?按照步驟6恢復(fù)融合態(tài)S1，對(duì)三維頂點(diǎn)的Z軸值進(jìn)行提取，得到比特序列Cre1。

第5步?對(duì)得到的Cre1、Cre2、Cre3三個(gè)比特序列進(jìn)行對(duì)比、校驗(yàn)得出最終的比特序列Ce，即秘密信息。

2?實(shí)驗(yàn)與結(jié)果分析

對(duì)本文算法進(jìn)行仿真實(shí)驗(yàn)，實(shí)驗(yàn)環(huán)境為Matlab2014，3dMax2012。

所選載體是由12個(gè)三維模型、隱藏圖像為64×64的baboon二值圖像和含有秘密信息的三維模型組成（如圖10所示）。

2.1?不可見性

2.1.1?不可見性分析

已有基于三維模型的信息隱藏算法的信息嵌入為修改頂點(diǎn)到中心的距離，而本文的隱藏算法是利用多個(gè)三維模型進(jìn)行融合得出一個(gè)新的三維模型，需要隱藏的信息被分配到多個(gè)三維模型上面，這樣使得每個(gè)模型中都有隱藏的信息存在，提高了算法的不可見性。

2.1.2?不可見性實(shí)驗(yàn)

1）Hausdorff距離。

Hausdorff距離[11]是度量兩個(gè)點(diǎn)集間的最大不匹配程度，距離越小，則表示匹配程度越高。

圖11為算法基于Hausdorff距離的不可見性實(shí)驗(yàn)對(duì)比結(jié)果。

文獻(xiàn)[12]算法是基于球型分割（Spherical Segmentation， SS）的算法，該算法首先對(duì)三維模型進(jìn)行主元分析、球形分割等預(yù)處理，以立體分區(qū)中法向量變化較大的點(diǎn)為特征點(diǎn)，根據(jù)待嵌入秘密信息量對(duì)特征點(diǎn)進(jìn)行小波變換，最后將秘密信息嵌入到載體中。文獻(xiàn)[13]算法是通過處理非重疊的2×2像素塊，利用修改方向技術(shù)將秘密數(shù)據(jù)嵌入到每個(gè)四像素塊中，不僅可以動(dòng)態(tài)修改像素值，還可以改變嵌入的數(shù)字，提高了秘密信息嵌入的魯棒性。本文實(shí)驗(yàn)使用的三維模型共含有2341874個(gè)頂點(diǎn)，取各個(gè)頂點(diǎn)3個(gè)坐標(biāo)值十進(jìn)制的尾數(shù)進(jìn)行信息隱藏，即可容納性約為2341874×3=223bit。以文中選定的局部高度為聚類參數(shù)對(duì)頂點(diǎn)進(jìn)行魯棒性劃分，所得到亞魯棒點(diǎn)約為1159017個(gè)，即用于隱藏秘密信息的容量約為1159017×3=222bit。由圖11可知，本文算法的Hausdorff距離在k>10時(shí)，除去k=20時(shí)的數(shù)據(jù)噪點(diǎn)，始終低于文獻(xiàn)[12-13]算法，表明本文算法相比基于文獻(xiàn)[12-13]算法，在嵌入量較大時(shí)，不可見性較好。

2）骨架相識(shí)度匹配。

骨架相似度En是另一種描述的是兩個(gè)三維模型的骨架相似程度，En越大，表明三維模型越相似，算法中對(duì)模型的修改程度越小。圖12為基于骨架相似度En的不可見性實(shí)驗(yàn)對(duì)比結(jié)果，從中可看出，在開始信息隱藏時(shí)，本文算法的骨架相似度匹配曲線位置除了部分坐標(biāo)值相同，其他始終高于文獻(xiàn)[12-13]算法，表明本文算法相比文獻(xiàn)[12-13]算法，不可見性較好。

2.2?魯棒性

2.2.1?魯棒性分析

魯棒性是衡量對(duì)三維模型經(jīng)過處理后的改變程度，反映信息的修改程度。三維模型可以抵抗旋轉(zhuǎn)、平移、等比例縮放等變換，在本文中使用多個(gè)三維模型以及融合態(tài)對(duì)需要隱藏的信息進(jìn)行隱藏，增加了對(duì)隱藏信息的校驗(yàn)，使算法能抵抗一定程度的非均勻縮放、隨機(jī)加噪攻擊等攻擊。

2.2.2?魯棒性實(shí)驗(yàn)

本文采用式（7）正確信息的比特?cái)?shù)比率（Bit Correct Ration， BCR）作為算法魯棒性的客觀評(píng)判標(biāo)準(zhǔn)：

BCR=正確的秘密信息比特?cái)?shù)秘密信息總比特?cái)?shù)（7）

本文實(shí)驗(yàn)中首先對(duì)含密的三維模型進(jìn)行旋轉(zhuǎn)、均勻縮放、隨機(jī)加噪、不均勻壓縮等仿真實(shí)驗(yàn)，然后對(duì)進(jìn)行干擾處理的三維模型進(jìn)行秘密信息提取，計(jì)算BCR。BCR越高，則算法的魯棒性越強(qiáng);否則算法的魯棒性越弱。

表1給出模型在遭受不同程度的攻擊之后，從模型中提取秘密信息并計(jì)算BCR的結(jié)果。

從表1可看出，模型對(duì)旋轉(zhuǎn)、不均勻壓縮、隨機(jī)加噪等攻擊有較強(qiáng)的魯棒性。

根據(jù)魯棒性檢驗(yàn)算法，將本文算法與文獻(xiàn)[13-14]算法進(jìn)行比較，本文重點(diǎn)對(duì)非均勻縮放攻擊和隨機(jī)加噪攻擊進(jìn)行了對(duì)比，圖13～14分別給出在非均勻壓縮、隨機(jī)加噪等攻擊下與文獻(xiàn)[13-14]算法的BCR對(duì)比。

由圖13～14可看出本文算法魯棒性在非均勻壓縮、隨機(jī)加噪等攻擊下的曲線的坐標(biāo)值，即BCR數(shù)值均高于文獻(xiàn)[13-14]算法，表明本文算法對(duì)非均勻壓縮、隨機(jī)加噪攻擊具有很強(qiáng)的魯棒性。

3?結(jié)語

本文基于融合態(tài)的思想對(duì)三維模型進(jìn)行預(yù)處理和旋轉(zhuǎn)操作得到多個(gè)融合態(tài)，再利用局部高度理論和Mean Shift聚類技術(shù)對(duì)融合態(tài)頂點(diǎn)進(jìn)行能量劃分，最后通過修改頂點(diǎn)的方法實(shí)現(xiàn)信息隱藏。由于本文算法使用的是多個(gè)三維模型，而且是對(duì)頂點(diǎn)的個(gè)位數(shù)值進(jìn)行修改，所以能保證算法的不可見性和魯棒性。通過實(shí)驗(yàn)對(duì)比可知，在相同攻擊下，

本文算法比文獻(xiàn)[13-14]算法具有更強(qiáng)的魯棒性，尤其是抵抗非均勻壓縮，達(dá)到了對(duì)信息隱藏的要求。

參考文獻(xiàn)（References）

[1] 齊向明， 史雙宇， 楊曉陶.基于特征點(diǎn)提取的三維點(diǎn)云模型水印算法[J]. 計(jì)算機(jī)應(yīng)用， 2014， 34（5）： 1309-1312， 1353. （QI X M， SHI S Y， YANG X T. Watermarking algorithm of 3D point cloud model based on feature point extraction[J]. Journal of Computer Applications， 2014， 34（5）： 1309-1312， 1353.）

[2] 廖學(xué)良， 王瑀屏.一種新的三維模型水印嵌入空域算法[J]. 計(jì)算機(jī)學(xué)報(bào)， 2008， 31（10）： 1848-1856. （LIAO X L， WANG Y P. A new 3D model watermark embedding spatial algorithm[J]. Chinese Journal of Computers， 2008， 31（10）： 1848-1856.）

[3] WANG Y P， HU S M. A new watermarking method for 3D models based on integral invariants[J]. IEEE Transactions on Visualization & Computer Graphics， 2009， 15（2）： 285.

[4] 任帥， 趙祥模， 張弢， 等.基于局部高度與Mean Shift的三維模型信息隱藏算法[J]. 計(jì)算機(jī)科學(xué)， 2017， 44（3）： 187-191. （REN S， ZHAO X M， ZHANG T， et al. Information hiding algorithm for 3D model based on local height and Mean Shift[J]. Computer Science， 2017， 44（3）： 187-191.）

[5] OHBUCHI R， MASUDA H， AONO M. Watermarking three-dimensional polygonal models[C]// Proceedings of the Fifth ACM International Conference on Multimedia. New York： ACM， 1997： 261-272.

[6] 張建明， 周小梅， 王新宇， 等.局部特征點(diǎn)的3維模型變換域水印算法[J]. 中國圖象圖形學(xué)報(bào)， 2014， 19（4）： 613-621. （ZHANG J M， ZHOU X M， WANG X Y， et al. Three-dimensional model transform domain watermarking algorithm for local feature points[J]. Journal of Image and Graphics， 2014， 19（4）： 613-621.）

[7] CHENG Y M， WANG C M. An adaptive steganographic algorithm for 3D polygonal meshes[J]. The Visual Computer， 2007， 23（9/10/11）： 721-732.

[8] 翟浩， 楊有， 萬輝.基于頂點(diǎn)排序的三維模型數(shù)字水印算法[J]. 重慶文理學(xué)院學(xué)報(bào)（社會(huì)科學(xué)版）， 2014， 33（5）： 128-133. （ZHAI H， YANG Y， WAN H. Digital watermarking algorithm based on vertex sorting for 3D models[J]. Journal of Chongqing University of Arts and Sciences （Social Sciences Edition）， 2014， 33（5）： 128-133.）

[9] 林金杰，朱代輝，楊育彬.3維模型局部高度研究[J].中國圖象圖形學(xué)報(bào)， 2011， 16（10）： 1841-1849.（LIN J J， ZHU D H， YANG Y B. Three-dimensional model study on local height[J]. Journal of Image & Graphics， 2011， 16（10）： 1841-1849.）

[10] COMANICIU D， MEER P. Mean shift analysis and applications[C]// Proceedings of the Seventh IEEE International Conference on Computer Vision. Piscataway， NJ： IEEE， 2002： 1197.

[11] YI X， CAMPS O I. Line-based recognition using a multidimensional hausdorff distance[J]. IEEE Transactions on Pattern Analysis and Machine Intelligence， 1999， 21（9）： 901-916.

[12] 任帥， 張弢， 楊濤， 等.基于三維模型球型分割的信息隱藏算法[J]. 計(jì)算機(jī)應(yīng)用， 2017， 37（9）： 2576-2580. （REN S， ZHANG T， YANG T， et al. Information hiding algorithm based on three-dimensional model spherical segmentation[J]. Journal of Computer Applications， 2017， 37（9）： 2576-2580.）

[13] LIAO X， DING L. Data hidng in digital image using four pixel-value-differencing and multiple-base notational[C]// Proceedings of the 2015 International Conference on Intelligent Information Hiding and Multimedia Signal Processing. Piscataway， NJ： IEEE， 2016： 76-80.

[14] 張弢， 慕德俊， 任帥， 等.利用內(nèi)切球解析的三維模型信息隱藏算法[J]. 西安電子科技大學(xué)學(xué)報(bào)（自然科學(xué)版）， 2014， 41（2）： 185-190. （ZHANG T， MU D J， REN S， et al. Three-dimensional model information hiding algorithm using inner sphere parsing[J]. Journal of Xidian University， 2014， 41（2）： 185-190.）