主权项 |
1.一种基于双随机相位编码和干涉原理的图像加密方法,其特征是按如下步骤进行:(1)加密:(i)f(x,y)和g(x,y)是待加密的两幅原始图像,R(x,y)是一块随机相位板,可以具体表示成exp[2πm(x,y)],其中m(x,y)代表在区间[0,1]上具有均匀概率分布的随机矩阵,将图像g(x,y)分解成两块相位板,一块为R(x,y),另一块相位板为P<sub>0</sub>(x,y),即有如下关系:g(x,y)=PT{R(x,y)+P<sub>0</sub>(x,y)} (1)其中PT{}表示切相操作,(x,y)代表空间域的坐标,切相操作的结果是除去复振幅的相位信息,只保留振幅信息,将P<sub>0</sub>(x,y)作为g(x,y)的加密结果,用函数<img file="FSA00000902794400011.GIF" wi="373" he="60" />具体表述,则R(x,y)为解密密钥,从式(1)可得到<img file="FSA00000902794400012.GIF" wi="174" he="59" />的一个解:<img file="FSA00000902794400013.GIF" wi="1467" he="151" />(ii)f(x,y)与随机相位板P<sub>0</sub>(x,y)相乘,进行一次阶数为α的分数傅立叶变换,切相操作后的结果为f′(u,υ)=PT{F<sup>α</sup>[f(x,y)·P<sub>0</sub>(x,y)]} (3)切除的相位可以表示为P<sub>1</sub>(u,υ)=PR{F<sup>α</sup>[f(x,y)·P<sub>0</sub>(x,y)]} (4)其中PR{}表示相位保留运算,即除去复振幅的振幅信息,只保留相位信息,(u,υ)代表频域坐标,F<sup>α</sup>[]代表阶数为α的分数傅立叶变换,函数f(x,y)·P<sub>0</sub>(x,y)的α阶分数傅立叶变换定义为:<maths num="0001"><![CDATA[<math><mrow><msup><mi>F</mi><mi>α</mi></msup><mo>[</mo><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>·</mo><msub><mi>P</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>]</mo><mrow><mo>(</mo><mi>u</mi><mo>,</mo><mi>υ</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mo>∫</mo><mrow><mo>-</mo><mo>∞</mo></mrow><mrow><mo>+</mo><mo>∞</mo></mrow></msubsup><msub><mi>K</mi><mi>α</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>;</mo><mi>u</mi><mo>,</mo><mi>υ</mi><mo>)</mo></mrow><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><msub><mi>P</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mi>dxdy</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中K<sub>α</sub>(x,y;u,υ)是二维分数傅立叶的变换核,即<maths num="0002"><![CDATA[<math><mrow><msub><mi>K</mi><mi>α</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>u</mi><mo>)</mo></mrow><mo>=</mo><mi>Aexp</mi><mrow><mo>(</mo><mi>iπ</mi><mfrac><mrow><msubsup><mi>x</mi><mn>1</mn><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>y</mi><mn>1</mn><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>x</mi><mn>2</mn><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>y</mi><mn>2</mn><mn>2</mn></msubsup></mrow><mrow><mi>λ</mi><mi>f</mi><mi>tan</mi><mi>φ</mi></mrow></mfrac><mo>-</mo><mn>2</mn><mi>iπ</mi><mfrac><mrow><msub><mi>x</mi><mn>1</mn></msub><msub><mi>y</mi><mn>1</mn></msub><msub><mi>x</mi><mn>2</mn></msub><msub><mi>y</mi><mn>2</mn></msub></mrow><mrow><msup><mi>λ</mi><mn>2</mn></msup><msup><mi>f</mi><mn>2</mn></msup><mi>sin</mi><mi>φ</mi></mrow></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中<img file="FSA00000902794400022.GIF" wi="747" he="153" />且φ=απ/2,α是分数形式的阶数;(iii)生成一块相位板P<sub>2</sub>(u,υ),其值为P<sub>2</sub>(u,υ)=P<sub>1</sub>(u,υ)·R<sup>*</sup>(u,υ) (7)其中“*”表示共轭,f′(u,υ)与相位板P<sub>2</sub>(u,υ)相乘,进行一次阶数为α的分数傅立叶变换后再进行切相操作,得到加密结果E(x,y)=PT{F<sup>α</sup>[f′(u,υ)·P<sub>2</sub>(u,υ)]} (8)切除的相位可以作为解密过程中的密钥,表示为P<sub>3</sub>(u,υ)=PR{F<sup>α</sup>[f′(u,υ)·P<sub>2</sub>(u,υ)]} (9)(2)解密:(i)将加密结果E(x,y)与作为解密密钥的相位板P<sub>3</sub>(x,y)相乘后进行一次阶数为-α的分数傅立叶变换,则由式(8)、(9)可知变换后的结果为F<sup>-α</sup>[E(x,y)·P<sub>3</sub>(x,y)]=f′(u,υ)·P<sub>2</sub>(u,υ);(ii)f′(u,υ)·P<sub>2</sub>(u,υ)与解密密钥R(u,υ)相乘后实行一次阶数为-α的傅立叶变换,由式(3),(4),(7)可得变换后的结果为f(x,y)·P<sub>0</sub>(x,y);(iii)对f(x,y)·P<sub>0</sub>(x,y)进行切相操作后得到第一幅解密结果f(x,y),作相位保留操作则得到P<sub>0</sub>(x,y),将P<sub>0</sub>(x,y)与解密密钥R(x,y)相加后作切相操作,由式(1)可知,其结果为g(x,y),即解密得到另一幅原始图像。综合以上各过程可以看出,两幅原始图像f(x,y)和g(x,y)最终被加密成一振幅分布E(x,y),解密过程中需要使用的解密密钥为R(u,υ)、P<sub>3</sub>(x,y)和分数傅立叶阶数-α。 |