发明名称 面向手机服务的数字图像选择加密方法
摘要 一种面向手机服务的数字图像选择加密方法。该方法包括图像小波子带的选择、混沌序列初始值的确定、图像加密及加密质量评价四部分过程。采用计算原始图像各尺度分数小波低频子带的重要性值来选择重要的分数小波低频子带,将该子带作为待加密的图像低频子带,而不是对所有的图像子带进行加密;采用四个参数的集合作为图像加密的密钥,其中对于被选择的待加密分数小波低频子带,其混沌加密序列初始值的确定完全依赖于所选择的重要低频子带;采用Chebyshev映射而不是Logistic映射对所选择的低频子带进行加密;采用峰值信噪比的倒数与图像的边缘差异比之和作为图像加密质量的量化评价,而不是仅仅通过肉眼观察来判断图像加密的质量。
申请公布号 CN102831569A 申请公布日期 2012.12.19
申请号 CN201210294355.7 申请日期 2012.08.17
申请人 华中师范大学 发明人 金聪;金枢炜
分类号 G06T1/00(2006.01)I 主分类号 G06T1/00(2006.01)I
代理机构 湖北武汉永嘉专利代理有限公司 42102 代理人 张安国;伍见
主权项 1.一种面向手机服务的数字图像选择加密方法,其特征在于:包括待加密分数小波低频子带的选择、图像混沌加密密钥的确定、图像的加密操作及加密质量评价四部分;(一)待加密分数小波低频子带的选择步骤1.对于给定的大小为M<sub>1</sub>×M<sub>2</sub>的原始图像I,进行β尺度分数小波变换;步骤2.对于所获得的β个分数小波低频子带,按照下列公式分别计算每个分数小波低频子带的重要性:<maths num="0001"><![CDATA[<math><mrow><msubsup><mi>SIG</mi><mi>&beta;</mi><mrow><mi>&theta;</mi><mo>=</mo><mn>4</mn></mrow></msubsup><mo>=</mo><mfrac><mn>1</mn><mrow><mrow><mo>(</mo><msub><mi>M</mi><mn>1</mn></msub><mo>/</mo><msup><mn>2</mn><mi>&beta;</mi></msup><mo>)</mo></mrow><mo>&times;</mo><mrow><mo>(</mo><msub><mi>M</mi><mn>2</mn></msub><mo>/</mo><msup><mn>2</mn><mi>&beta;</mi></msup><mo>)</mo></mrow></mrow></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><msub><mi>M</mi><mn>1</mn></msub><mo>/</mo><msup><mn>2</mn><mi>&beta;</mi></msup></mrow></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mrow><msub><mi>M</mi><mn>2</mn></msub><mo>/</mo><msup><mn>2</mn><mi>&beta;</mi></msup></mrow></munderover><msup><mrow><mo>{</mo><msubsup><mi>I</mi><mi>&beta;</mi><mrow><mi>&theta;</mi><mo>=</mo><mn>4</mn></mrow></msubsup><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>}</mo></mrow><mn>2</mn></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,θ∈{1,2,3,4}分别表示水平子带、垂直子带、对角线子带和低频子带,(M<sub>1</sub>/2<sup>β</sup>)×(M<sub>2</sub>/2<sup>β</sup>)是小波低频子带<img file="FDA00002024208800012.GIF" wi="94" he="70" />的大小,<img file="FDA00002024208800013.GIF" wi="153" he="70" />是分数小波低频子带<img file="FDA00002024208800014.GIF" wi="94" he="70" />的重要性;步骤3.取选择阈值<img file="FDA00002024208800015.GIF" wi="335" he="63" />其中0&lt;Q&lt;1是调节参数;步骤4.将每个小波低频子带的重要性与选择阈值α进行比较,进行如下操作:若满足不等式<img file="FDA00002024208800016.GIF" wi="277" he="71" />则判定此小波子带<img file="FDA00002024208800017.GIF" wi="94" he="71" />是重要的;如果不满足不等式<img file="FDA00002024208800018.GIF" wi="268" he="70" />则认为是不重要的;步骤5.若判断是重要的分数小波低频子带多于一个,则选择重要性最大的分数小波低频子带为最终选择的待加密小波低频子带;若具有最大重要性的小波低频子带不止一个,则任选其一作为待加密小波低频子带;若满足不等式<img file="FDA00002024208800019.GIF" wi="248" he="70" />的小波低频子带不存在,则调整调节参数Q的值,调整量由用户决定;(二)图像混沌加密密钥的确定要实现分数小波变换和混沌加密需要确定以下参数:(1)混沌序列的控制参数α,由用户给出;(2)图像小波变换的尺度β,由用户给出;(3)分数小波变换的分数阶数γ,且满足0&lt;γ&lt;2,由用户给出;(4)混沌序列初始值X(0),它的确定过程如下:步骤1.设原始图像是m比特的,从被选择的子带中任选5个像素点P<sub>1</sub>,P<sub>2</sub>,P<sub>3</sub>,P<sub>4</sub>,P<sub>5</sub>,由这5个点的像素值来获得该子带的混沌序列初始值;步骤2.按照下列公式计算所选小波子带的混沌序列初始值:<maths num="0002"><![CDATA[<math><mrow><mi>X</mi><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>P</mi><mn>1</mn></msub><mo>&times;</mo><msup><mn>10</mn><mn>4</mn></msup><mo>+</mo><msub><mi>P</mi><mn>2</mn></msub><mo>&times;</mo><msup><mn>10</mn><mn>3</mn></msup><mo>+</mo><msub><mi>P</mi><mn>3</mn></msub><mo>&times;</mo><msup><mn>10</mn><mn>2</mn></msup><mo>+</mo><msub><mi>P</mi><mn>4</mn></msub><mo>&times;</mo><msup><mn>10</mn><mn>1</mn></msup><mo>+</mo><msub><mi>P</mi><mn>5</mn></msub><mo>&times;</mo><msup><mn>10</mn><mn>0</mn></msup></mrow><msup><mn>10</mn><mi>m</mi></msup></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths>上述四个参数α,β,γ,X(0)构成所选子带图像加密的密钥;(三)图像的加密步骤1.对于所选择的小波低频子带,设该子带的大小是N<sub>1</sub>×N<sub>2</sub>,利用这个分数小波低频子带生成的系列初始值X(0),按照如下映射生成混沌序列:X(n+1)=cos(α·arccos(X(n))    (4)其中,α是控制参数,且X(n)∈[-1,1],此时所生成的序列是混沌的,由于这个混沌序列是双极性的,按照下述变换将其转化成[0,1]内的单极性序列:<maths num="0003"><![CDATA[<math><mrow><mi>Y</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mi>X</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths>这个序列的长度取为N<sub>1</sub>×N<sub>2</sub>,并将这N<sub>1</sub>×N<sub>2</sub>个数值转换成N<sub>1</sub>行N<sub>2</sub>列的混沌矩阵,记为U;步骤2.对所选择的分数小波低频子带进行加密,加密按如下公式计算:<maths num="0004"><![CDATA[<math><mrow><msub><mi>Value</mi><mi>New</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mi>U</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>&CircleTimes;</mo><msub><mi>Value</mi><mi>Old</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,<img file="FDA00002024208800032.GIF" wi="40" he="39" />是XOR运算,Value<sub>New</sub>(i,j)、Value<sub>Old</sub>(i,j)分别是所选择分数小波低频子带加密后与加密前在第i行第j列的低频子带分数小波系数;步骤3.用加密后的分数小波低频子带替换加密前的分数小波低频子带之后,进行小波逆变换,获得初始加密图像;步骤4.将所获得的初始加密图像分割成4幅大小相同的子图像,之后对这4幅子图像依照顺时针方向做轮换处理后,所获得的图像就是最终的加密图像;(四)图像选择加密方法的质量评价步骤1.按照下列公式计算峰值信噪比:<maths num="0005"><![CDATA[<math><mrow><mi>PSNR</mi><mo>=</mo><mn>10</mn><msub><mi>log</mi><mn>10</mn></msub><mfrac><mrow><mrow><mo>(</mo><msub><mi>M</mi><mn>1</mn></msub><mo>&times;</mo><msub><mi>M</mi><mn>2</mn></msub><mo>)</mo></mrow><mi>max</mi><mrow><mo>(</mo><mi>max</mi><mrow><mo>(</mo><msup><mi>I</mi><mn>2</mn></msup><mo>)</mo></mrow><mo>)</mo></mrow></mrow><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>M</mi><mn>1</mn></msub></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>M</mi><mn>2</mn></msub></munderover><msup><mrow><mo>(</mo><mi>I</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><msup><mi>I</mi><mo>*</mo></msup><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,I<sup>*</sup>是加密后的图像;步骤2.按照Canny算子计算原始图像和加密图像的边缘,所获得的原始图像I的二值边缘图像记为B,加密图像I<sup>*</sup>的二值边缘图像记为B<sup>*</sup>;步骤3.按照下列公式计算图像的边缘差异比:<maths num="0006"><![CDATA[<math><mrow><mi>EDR</mi><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>M</mi><mn>1</mn></msub></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>M</mi><mn>2</mn></msub></munderover><mo>|</mo><mi>B</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><msup><mi>B</mi><mo>*</mo></msup><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>|</mo></mrow><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>M</mi><mn>1</mn></msub></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>M</mi><mn>2</mn></msub></munderover><mrow><mo>(</mo><mi>B</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>+</mo><msup><mi>B</mi><mo>*</mo></msup><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow></math>]]></maths>步骤4.图像加密质量的评价值是1/PNSR+EDR,此值越大,则加密质量越好。
地址 430079 湖北省武汉市洪山区珞瑜路152号