发明名称 一种矢量几何模型压缩方法
摘要 本发明公开了一种矢量几何模型压缩方法。方法中采用随机矩阵对网格几何进行抽样并完成压缩,之后通过抽样信号恢复出原始信号的稀疏表达方式,再通过逆拉普拉斯算子恢复出原始的几何信号,完成解压缩。该压缩方法可以快速完成对矢量几何模型的压缩,并可以按照任意精度设定阈值,控制恢复出模型的效果。
申请公布号 CN102665027A 申请公布日期 2012.09.12
申请号 CN201210118488.9 申请日期 2012.04.20
申请人 西北大学 发明人 周明全;耿国华;杜卓明;李康;贺毅岳;史哲;郭竞;张婧
分类号 H04N1/41(2006.01)I 主分类号 H04N1/41(2006.01)I
代理机构 西安恒泰知识产权代理事务所 61216 代理人 史玫
主权项 1.一种矢量几何模型压缩方法,该方法中二维矢量几何模型的几何信息由几何信号x<sub>2</sub>和几何信号y<sub>2</sub>构成,三维矢量几何模型的几何信息由几何信号x<sub>3</sub>、几何信号y<sub>3</sub>和几何信号z<sub>3</sub>构成,方法具体通过下列步骤实现:(1)对于二维矢量几何模型:其拉普拉斯算子<maths num="0001"><![CDATA[<math><mrow><msub><mi>L</mi><mrow><msub><mi>n</mi><mn>1</mn></msub><mo>&times;</mo><msub><mi>n</mi><mn>1</mn></msub></mrow></msub><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac></mtd><mtd><mn>0</mn></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mn>0</mn></mtd><mtd><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac></mtd></mtr><mtr><mtd><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac></mtd><mtd><mn>1</mn></mtd><mtd><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac></mtd><mtd><mn>0</mn></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mn>0</mn></mtd><mtd><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac></mtd><mtd><mn>1</mn></mtd><mtd><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac></mtd></mtr><mtr><mtd><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac></mtd><mtd><mn>0</mn></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mn>0</mn></mtd><mtd><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths>n<sub>1</sub>为二维矢量几何模型的顶点总数,n<sub>1</sub>取正整数;对于三维矢量几何模型:其拉普拉斯算子<img file="FDA0000155546850000012.GIF" wi="428" he="149" />其中:A为三维矢量几何模型的邻接矩阵,D为三维矢量几何模型的顶点度矩阵,且<img file="FDA0000155546850000013.GIF" wi="831" he="558" />其中:d<sub>i</sub>为三维矢量几何模型的第i个顶点的度,n<sub>2</sub>为三维矢量几何模型的顶点总数,n<sub>2</sub>和i均取正整数;(2)对于二维矢量几何模型:将二维矢量几何模型的拉普拉斯算子<img file="FDA0000155546850000014.GIF" wi="114" he="65" />作用到二维矢量几何模型的几何信号x<sub>2</sub>,得到向量λ′<sub>1</sub>:<maths num="0002"><![CDATA[<math><mrow><msubsup><mi>&lambda;</mi><mn>1</mn><mo>&prime;</mo></msubsup><mo>=</mo><msub><mi>L</mi><mrow><msub><mi>n</mi><mn>1</mn></msub><mo>&times;</mo><msub><mi>n</mi><mn>1</mn></msub></mrow></msub><msub><mi>x</mi><mn>2</mn></msub><mo>;</mo></mrow></math>]]></maths>根据设定的阈值ε<sub>1</sub>,将向量λ′<sub>1</sub>中的绝对值小于ε<sub>1</sub>的元素赋值为0,得到几何信号x<sub>2</sub>的稀疏几何信号λ<sub>1</sub>;将二维矢量几何模型的拉普拉斯算子<img file="FDA0000155546850000021.GIF" wi="135" he="79" />作用到二维矢量几何模型的几何信号y<sub>2</sub>,得到向量λ′<sub>2</sub>:<maths num="0003"><![CDATA[<math><mrow><msubsup><mi>&lambda;</mi><mn>2</mn><mo>&prime;</mo></msubsup><mo>=</mo><msub><mi>L</mi><mrow><msub><mi>n</mi><mn>1</mn></msub><mo>&times;</mo><msub><mi>n</mi><mn>1</mn></msub></mrow></msub><msub><mi>y</mi><mn>2</mn></msub><mo>;</mo></mrow></math>]]></maths>根据设定的阈值ε<sub>2</sub>,将向量λ′<sub>2</sub>中的绝对值小于ε<sub>2</sub>的元素赋值为0,得到几何信号y<sub>2</sub>的稀疏几何信号λ<sub>2</sub>;其中:λ<sub>1</sub>和λ<sub>2</sub>的维数均为n<sub>1</sub>,ε<sub>1</sub>和ε<sub>2</sub>满足:λ<sub>1</sub>和λ<sub>2</sub>中的非0元素的个数相等;对于三维矢量几何模型:将三维矢量几何模型的拉普拉斯算子<img file="FDA0000155546850000023.GIF" wi="143" he="81" />作用到三维矢量几何模型的几何信号x<sub>3</sub>,得到向量λ′<sub>3</sub>:<maths num="0004"><![CDATA[<math><mrow><msubsup><mi>&lambda;</mi><mn>3</mn><mo>&prime;</mo></msubsup><mo>=</mo><msub><mi>L</mi><mrow><msub><mi>n</mi><mn>2</mn></msub><mo>&times;</mo><msub><mi>n</mi><mn>2</mn></msub></mrow></msub><msub><mi>x</mi><mn>3</mn></msub><mo>;</mo></mrow></math>]]></maths>根据设定的阈值ε<sub>3</sub>,将向量λ′<sub>3</sub>中的绝对值小于ε<sub>3</sub>的元素赋值为0,得到几何信号x<sub>3</sub>的稀疏几何信号λ<sub>3</sub>;将三维矢量几何模型的拉普拉斯算子<img file="FDA0000155546850000025.GIF" wi="143" he="81" />作用到三维矢量几何模型的几何信号y<sub>3</sub>,得到向量λ′<sub>4</sub>:<maths num="0005"><![CDATA[<math><mrow><msubsup><mi>&lambda;</mi><mn>4</mn><mo>&prime;</mo></msubsup><mo>=</mo><msub><mi>L</mi><mrow><msub><mi>n</mi><mn>2</mn></msub><mo>&times;</mo><msub><mi>n</mi><mn>2</mn></msub></mrow></msub><msub><mi>y</mi><mn>3</mn></msub><mo>;</mo></mrow></math>]]></maths>根据设定的阈值ε<sub>4</sub>,将向量λ′<sub>4</sub>中的绝对值小于ε<sub>4</sub>的元素赋值为0,得到几何信号y<sub>3</sub>的稀疏几何信号λ<sub>4</sub>;将三维矢量几何模型的拉普拉斯算子<img file="FDA0000155546850000027.GIF" wi="143" he="81" />作用到三维矢量几何模型的几何信号z<sub>3</sub>,得到向量λ′<sub>5</sub>:<maths num="0006"><![CDATA[<math><mrow><msubsup><mi>&lambda;</mi><mn>5</mn><mo>&prime;</mo></msubsup><mo>=</mo><msub><mi>L</mi><mrow><msub><mi>n</mi><mn>2</mn></msub><mo>&times;</mo><msub><mi>n</mi><mn>2</mn></msub></mrow></msub><msub><mi>z</mi><mn>3</mn></msub><mo>;</mo></mrow></math>]]></maths>根据设定的阈值ε<sub>5</sub>,将向量λ′<sub>5</sub>中的绝对值小于ε<sub>5</sub>的元素赋值为0,得到几何信号z<sub>3</sub>的稀疏几何信号λ<sub>5</sub>;其中:λ<sub>3</sub>、λ<sub>4</sub>和λ<sub>5</sub>的维数均为n<sub>2</sub>,ε<sub>3</sub>、ε<sub>4</sub>和ε<sub>5</sub>满足:λ<sub>3</sub>、λ<sub>4</sub>和λ<sub>5</sub>中的非0元素的个数相等;(3)对于二维矢量几何模型:记录稀疏几何信号λ<sub>1</sub>或稀疏几何信号λ<sub>2</sub>中的非0元素的个数r<sub>1</sub>,r<sub>1</sub><<n<sub>1</sub>;对于三维矢量几何模型:记录稀疏几何信号λ<sub>3</sub>、稀疏几何信号λ<sub>4</sub>或稀疏几何信号λ<sub>5</sub>中的非0元素的个数r<sub>2</sub>,r<sub>2</sub><<n<sub>2</sub>;步骤二,生成随机矩阵对矢量几何模型的几何信息进行抽样:对于二维矢量几何模型:生成随机抽样矩阵<img file="FDA0000155546850000031.GIF" wi="211" he="94" />利用<img file="FDA0000155546850000032.GIF" wi="172" he="94" />对二维矢量几何模型的几何信号x<sub>2</sub>进行抽样,得到抽样后的信号θ<sub>1</sub>:<maths num="0007"><![CDATA[<math><mrow><msub><mi>&theta;</mi><mn>1</mn></msub><mo>=</mo><msub><mi>S</mi><mrow><msub><mrow><mn>4</mn><mi>r</mi></mrow><mn>1</mn></msub><mo>&times;</mo><msub><mi>n</mi><mn>1</mn></msub></mrow></msub><msub><mi>x</mi><mn>2</mn></msub><mo>;</mo></mrow></math>]]></maths>利用<img file="FDA0000155546850000034.GIF" wi="172" he="94" />对二维矢量几何模型的几何信号y<sub>2</sub>进行抽样,得到抽样后的信号θ<sub>2</sub>:<maths num="0008"><![CDATA[<math><mrow><msub><mi>&theta;</mi><mn>2</mn></msub><mo>=</mo><msub><mi>S</mi><mrow><msub><mrow><mn>4</mn><mi>r</mi></mrow><mn>1</mn></msub><mo>&times;</mo><msub><mi>n</mi><mn>1</mn></msub></mrow></msub><msub><mi>y</mi><mn>2</mn></msub><mo>;</mo></mrow></math>]]></maths>其中:θ<sub>1</sub>和θ<sub>2</sub>长度均为4r<sub>1</sub>,4r<sub>1</sub><<n<sub>1</sub>,由此完成压缩;对于三维矢量几何模型:生成随机抽样矩阵<img file="FDA0000155546850000036.GIF" wi="219" he="93" />利用<img file="FDA0000155546850000037.GIF" wi="186" he="94" />对三维矢量几何模型的几何信号x<sub>3</sub>进行抽样,得到抽样后的信号θ<sub>3</sub>,<maths num="0009"><![CDATA[<math><mrow><msub><mi>&theta;</mi><mn>3</mn></msub><mo>=</mo><msub><mi>S</mi><mrow><msub><mrow><mn>4</mn><mi>r</mi></mrow><mn>2</mn></msub><mo>&times;</mo><msub><mi>n</mi><mn>2</mn></msub></mrow></msub><msub><mi>x</mi><mn>3</mn></msub><mo>;</mo></mrow></math>]]></maths>利用<img file="FDA0000155546850000039.GIF" wi="186" he="94" />对三维矢量几何模型的几何信号y<sub>3</sub>进行抽样,得到抽样后的信号θ<sub>4</sub>,<maths num="0010"><![CDATA[<math><mrow><msub><mi>&theta;</mi><mn>4</mn></msub><mo>=</mo><msub><mi>S</mi><mrow><msub><mrow><mn>4</mn><mi>r</mi></mrow><mn>2</mn></msub><mo>&times;</mo><msub><mi>n</mi><mn>2</mn></msub></mrow></msub><msub><mi>y</mi><mn>3</mn></msub><mo>;</mo></mrow></math>]]></maths>利用<img file="FDA0000155546850000041.GIF" wi="186" he="93" />对三维矢量几何模型的几何信号z<sub>3</sub>进行抽样,得到抽样后的信号θ<sub>5</sub>,<img file="FDA0000155546850000042.GIF" wi="415" he="88" />其中:θ<sub>3</sub>、θ<sub>4</sub>和θ<sub>5</sub>长度均为4r<sub>2</sub>,4r<sub>2</sub><<n<sub>2</sub>;进而完成压缩;步骤三,进行模型的恢复:对于二维矢量几何模型:利用最小化稀疏信号的1范数,通过抽样后的信号θ<sub>1</sub>恢复出几何信号x<sub>2</sub>的稀疏几何信号λ<sub>1</sub>:<maths num="0011"><![CDATA[<math><mrow><mi>min</mi><msub><mrow><mo>|</mo><mo>|</mo><msub><mi>&lambda;</mi><mn>1</mn></msub><mo>|</mo><mo>|</mo></mrow><mn>1</mn></msub></mrow></math>]]></maths><maths num="0012"><![CDATA[<math><mrow><mi>s</mi><mo>.</mo><mi>t</mi><msub><mi>&theta;</mi><mn>1</mn></msub><mo>=</mo><msup><mrow><msub><mi>S</mi><mrow><msub><mrow><mn>4</mn><mi>r</mi></mrow><mn>1</mn></msub><mo>&times;</mo><msub><mi>n</mi><mn>1</mn></msub></mrow></msub><msub><mi>L</mi><mrow><msub><mi>n</mi><mn>1</mn></msub><mo>&times;</mo><msub><mi>n</mi><mn>1</mn></msub><mo></mo></mrow></msub></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><msub><mi>&lambda;</mi><mn>1</mn></msub><mo>;</mo></mrow></math>]]></maths>利用最小化稀疏信号的1范数,通过抽样后的信号θ<sub>2</sub>恢复出几何信号y<sub>2</sub>的稀疏几何信号λ<sub>2</sub>:<maths num="0013"><![CDATA[<math><mrow><mi>min</mi><msub><mrow><mo>|</mo><mo>|</mo><msub><mi>&lambda;</mi><mn>2</mn></msub><mo>|</mo><mo>|</mo></mrow><mn>1</mn></msub></mrow></math>]]></maths><maths num="0014"><![CDATA[<math><mrow><mi>s</mi><mo>.</mo><mi>t</mi><msub><mi>&theta;</mi><mn>2</mn></msub><mo>=</mo><msup><mrow><msub><mi>S</mi><mrow><msub><mrow><mn>4</mn><mi>r</mi></mrow><mn>1</mn></msub><mo>&times;</mo><msub><mi>n</mi><mn>1</mn></msub></mrow></msub><msub><mi>L</mi><mrow><msub><mi>n</mi><mn>1</mn></msub><mo>&times;</mo><msub><mi>n</mi><mn>1</mn></msub><mo></mo></mrow></msub></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><msub><mi>&lambda;</mi><mn>2</mn></msub><mo>;</mo></mrow></math>]]></maths>然后,利用二维矢量几何模型的逆拉普拉斯算子<img file="FDA0000155546850000047.GIF" wi="115" he="97" />和λ<sub>1</sub>恢复出原始的几何信号x<sub>2</sub>:<maths num="0015"><![CDATA[<math><mrow><msub><mi>x</mi><mn>2</mn></msub><mo>=</mo><msup><msub><mi>L</mi><mrow><msub><mi>n</mi><mn>1</mn></msub><mo>&times;</mo><msub><mi>n</mi><mn>1</mn></msub></mrow></msub><mrow><mo>-</mo><mn>1</mn></mrow></msup><msub><mi>&lambda;</mi><mn>1</mn></msub><mo>;</mo></mrow></math>]]></maths>利用二维矢量几何模型的逆拉普拉斯算子<img file="FDA0000155546850000049.GIF" wi="115" he="97" />和λ<sub>2</sub>恢复出原始的几何信号y<sub>2</sub>:<maths num="0016"><![CDATA[<math><mrow><msub><mi>y</mi><mn>2</mn></msub><mo>=</mo><msup><msub><mi>L</mi><mrow><msub><mi>n</mi><mn>1</mn></msub><mo>&times;</mo><msub><mi>n</mi><mn>1</mn></msub></mrow></msub><mrow><mo>-</mo><mn>1</mn></mrow></msup><msub><mi>&lambda;</mi><mn>2</mn></msub><mo>;</mo></mrow></math>]]></maths>对于三维矢量几何模型:利用最小化稀疏信号的1范数,通过抽样后的信号θ<sub>3</sub>恢复出几何信号x<sub>3</sub>的稀疏几何信号λ<sub>3</sub>:<maths num="0017"><![CDATA[<math><mrow><mi>min</mi><msub><mrow><mo>|</mo><mo>|</mo><msub><mi>&lambda;</mi><mn>3</mn></msub><mo>|</mo><mo>|</mo></mrow><mn>1</mn></msub></mrow></math>]]></maths><maths num="0018"><![CDATA[<math><mrow><mi>s</mi><mo>.</mo><mi>t</mi><msub><mi>&theta;</mi><mn>3</mn></msub><mo>=</mo><msup><mrow><msub><mi>S</mi><mrow><msub><mrow><mn>4</mn><mi>r</mi></mrow><mn>2</mn></msub><mo>&times;</mo><msub><mi>n</mi><mn>2</mn></msub></mrow></msub><msub><mi>L</mi><mrow><msub><mi>n</mi><mn>2</mn></msub><mo>&times;</mo><msub><mi>n</mi><mn>2</mn></msub><mo></mo></mrow></msub></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><msub><mi>&lambda;</mi><mn>3</mn></msub><mo>;</mo></mrow></math>]]></maths>利用最小化稀疏信号的1范数,通过抽样后的信号θ<sub>4</sub>恢复出几何信号y<sub>3</sub>的稀疏几何信号λ<sub>4</sub>:<maths num="0019"><![CDATA[<math><mrow><mi>min</mi><msub><mrow><mo>|</mo><mo>|</mo><msub><mi>&lambda;</mi><mn>4</mn></msub><mo>|</mo><mo>|</mo></mrow><mn>1</mn></msub></mrow></math>]]></maths><maths num="0020"><![CDATA[<math><mrow><mi>s</mi><mo>.</mo><mi>t</mi><msub><mi>&theta;</mi><mn>4</mn></msub><mo>=</mo><msup><mrow><msub><mi>S</mi><mrow><msub><mrow><mn>4</mn><mi>r</mi></mrow><mn>2</mn></msub><mo>&times;</mo><msub><mi>n</mi><mn>2</mn></msub></mrow></msub><msub><mi>L</mi><mrow><msub><mi>n</mi><mn>2</mn></msub><mo>&times;</mo><msub><mi>n</mi><mn>2</mn></msub><mo></mo></mrow></msub></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><msub><mi>&lambda;</mi><mn>4</mn></msub><mo>;</mo></mrow></math>]]></maths>利用最小化稀疏信号的1范数,通过抽样后的信号θ<sub>5</sub>恢复出几何信号z<sub>3</sub>的稀疏几何信号λ<sub>5</sub>:<maths num="0021"><![CDATA[<math><mrow><mi>min</mi><msub><mrow><mo>|</mo><mo>|</mo><msub><mi>&lambda;</mi><mn>5</mn></msub><mo>|</mo><mo>|</mo></mrow><mn>1</mn></msub></mrow></math>]]></maths><maths num="0022"><![CDATA[<math><mrow><mi>s</mi><mo>.</mo><mi>t</mi><msub><mi>&theta;</mi><mn>5</mn></msub><mo>=</mo><msup><mrow><msub><mi>S</mi><mrow><msub><mrow><mn>4</mn><mi>r</mi></mrow><mn>2</mn></msub><mo>&times;</mo><msub><mi>n</mi><mn>2</mn></msub></mrow></msub><msub><mi>L</mi><mrow><msub><mi>n</mi><mn>2</mn></msub><mo>&times;</mo><msub><mi>n</mi><mn>2</mn></msub><mo></mo></mrow></msub></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><msub><mi>&lambda;</mi><mn>5</mn></msub><mo>;</mo></mrow></math>]]></maths>然后,利用三维矢量几何模型的逆拉普拉斯算子<img file="FDA0000155546850000055.GIF" wi="127" he="101" />和稀疏几何信号λ<sub>3</sub>恢复出原始的几何信号x<sub>3</sub>:<maths num="0023"><![CDATA[<math><mrow><msub><mi>x</mi><mn>3</mn></msub><mo>=</mo><msup><msub><mi>L</mi><mrow><msub><mi>n</mi><mn>2</mn></msub><mo>&times;</mo><msub><mi>n</mi><mn>2</mn></msub></mrow></msub><mrow><mo>-</mo><mn>1</mn></mrow></msup><msub><mi>&lambda;</mi><mn>3</mn></msub><mo>;</mo></mrow></math>]]></maths>利用三维矢量几何模型的逆拉普拉斯算子<img file="FDA0000155546850000057.GIF" wi="126" he="101" />和稀疏几何信号λ<sub>4</sub>恢复出原始的几何信号y<sub>3</sub>:<maths num="0024"><![CDATA[<math><mrow><msub><mi>y</mi><mn>3</mn></msub><mo>=</mo><msup><msub><mi>L</mi><mrow><msub><mi>n</mi><mn>2</mn></msub><mo>&times;</mo><msub><mi>n</mi><mn>2</mn></msub></mrow></msub><mrow><mo>-</mo><mn>1</mn></mrow></msup><msub><mi>&lambda;</mi><mn>4</mn></msub><mo>;</mo></mrow></math>]]></maths>利用三维矢量几何模型的逆拉普拉斯算子<img file="FDA0000155546850000059.GIF" wi="126" he="101" />和稀疏几何信号λ<sub>5</sub>恢复出原始的几何信号z<sub>3</sub>:<maths num="0025"><![CDATA[<math><mrow><msub><mi>z</mi><mn>3</mn></msub><mo>=</mo><msup><msub><mi>L</mi><mrow><msub><mi>n</mi><mn>2</mn></msub><mo>&times;</mo><msub><mi>n</mi><mn>2</mn></msub></mrow></msub><mrow><mo>-</mo><mn>1</mn></mrow></msup><msub><mi>&lambda;</mi><mn>5</mn></msub><mo>.</mo></mrow></math>]]></maths>
地址 710069 陕西省西安市太白北路229号