CN104573561B - 基于稀疏双随机相位加密和qr码的认证方法 - Google Patents

基于稀疏双随机相位加密和qr码的认证方法 Download PDF

Info

Publication number
CN104573561B
CN104573561B CN201510050296.2A CN201510050296A CN104573561B CN 104573561 B CN104573561 B CN 104573561B CN 201510050296 A CN201510050296 A CN 201510050296A CN 104573561 B CN104573561 B CN 104573561B
Authority
CN
China
Prior art keywords
mrow
msub
msup
prime
phase
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201510050296.2A
Other languages
English (en)
Other versions
CN104573561A (zh
Inventor
汪小刚
周国泉
戴朝卿
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang A&F University ZAFU
Original Assignee
Zhejiang A&F University ZAFU
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Zhejiang A&F University ZAFU filed Critical Zhejiang A&F University ZAFU
Priority to CN201510050296.2A priority Critical patent/CN104573561B/zh
Publication of CN104573561A publication Critical patent/CN104573561A/zh
Application granted granted Critical
Publication of CN104573561B publication Critical patent/CN104573561B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F21/00—Security arrangements for protecting computers, components thereof, programs or data against unauthorised activity
    • G06F21/60—Protecting data
    • G06F21/602—Providing cryptographic facilities or services
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F21/00—Security arrangements for protecting computers, components thereof, programs or data against unauthorised activity
    • G06F21/60—Protecting data
    • G06F21/62—Protecting access to data via a platform, e.g. using keys or access control rules
    • G06F21/6209—Protecting access to data via a platform, e.g. using keys or access control rules to a single file or object, e.g. in a secure envelope, encrypted and accessed using a key, or with access control rules appended to the object itself
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2221/00—Indexing scheme relating to security arrangements for protecting computers, components thereof, programs or data against unauthorised activity
    • G06F2221/21—Indexing scheme relating to G06F21/00 and subgroups addressing additional information or applications relating to security arrangements for protecting computers, components thereof, programs or data against unauthorised activity
    • G06F2221/2107—File encryption

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Software Systems (AREA)
  • General Health & Medical Sciences (AREA)
  • Computer Hardware Design (AREA)
  • Computer Security & Cryptography (AREA)
  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Bioethics (AREA)
  • Two-Way Televisions, Distribution Of Moving Picture Or The Like (AREA)
  • Compression Or Coding Systems Of Tv Signals (AREA)

Abstract

一种基于稀疏双随机相位加密和QR码的认证方法。按如下两大步骤进行：一是加密：利用QR码和两块随机相位板将待加密图像加密成一振幅图像，通过相位恢复算法得到输出面相位的近似分布，最终通过随机抽取元素的方法生成稀疏双随机相位加密图像；二是解密与认证：首先利用解密密钥对稀疏双随机相位加密图像进行解密，解密后得到的图像与原始图像和QR码的乘积进行对比认证。本发明提出的安全认证方法结合了光学加密和QR码认证两种方式，因此具有很高的安全性。

Description

基于稀疏双随机相位加密和QR码的认证方法
【技术领域】
本发明涉及信息安全技术领域，特别是图像的安全认证方法，适用于图像的加密和认证。
【背景技术】
随着信息技术的不断发展，信息安全问题也日显突出。如何确保信息系统的安全已成为全社会关注的问题，它也成为信息科学的热点课题。基于光学理论和方法的信息加密技术是近些年逐步发展起来的新一代的信息安全处理技术。1995年，美国Connecticut大学的学者P.Refregier和B.Javidi提出了一种基于“双随机相位编码”的光学图像加密技术，它的提出引起了世界各国光学信息安全领域科研人员的极大关注。利用双随机相位加密技术对图像进行光学加密后，通过光子成像或者计算机数字技术对光学加密结果进行稀疏化处理，生成的稀疏加密图像可用于信息的安全认证。这种基于稀疏双随机相位加密的认证方法结合了光学加密和图像认证两种方式，安全性极高。
QR码的“QR”是Quick Response的缩写，即快速反应的意思。QR码能够对文字、URL地址和其他类型的数据进行加密。它具有信息容量大、可靠性高、保密防伪性强以及快速读取等优点。相位恢复算法则是一种通过可测量的光场强度确定光场相位分布的方法，它已被广泛应用在电子显微镜、波前再现、面形检测等领域。QR码和相位恢复算法目前都已经成功地运用于图像加密领域，它们在基于稀疏双随机相位加密的认证方法中的应用将降低光学加密系统的复杂性并提升认证过程的安全性。
【发明内容】
本发明要解决的技术问题是提供基于稀疏双随机相位加密和QR码的认证方法。
解决上述技术问题采用如下技术措施：基于稀疏双随机相位加密和QR码的认证方法按如下步骤进行：
(1)加密：
(i)Q(x，y)代表光学加密时作为加密密钥和认证时作为认证密钥的QR码，f(x，y)代表待加密并用于认证的原始图像，利用菲涅耳域双随机相位加密系统对f(x，y)进行加密，r1(x，y)和r2(x′，y′)是双随机相位加密系统中的两块随机相位板，分别表示成exp[2πa(x，y)]和exp[2πb(x′，y′)]，其中a(x，y)、b(x′，y′)代表两个统计无关且在区间[0，1]上具有均匀概率分布的随机矩阵，(x，y)和(x′，y′)分别表示输入平面和菲涅耳衍射平面的坐标，加密时图像f(x，y)和加密密钥Q(x，y)相乘，它们的乘积与r1(x，y)相乘后作一次波长为λ，距离为d1的菲涅耳变换，得到的结果与r2(x′，y′)相乘后再作一次波长为λ，距离为d2的菲涅耳变换，对变换后的结果取振幅后得到密文C(x″，y″)，即：
其中PT{ }代表取振幅操作，即除去复振幅的相位信息，只保留振幅部分的信息，FrT{ }代表菲涅耳变换，(x″，y″)表示菲涅耳域双随机相位加密系统输出平面的坐标；
(ii)利用相位恢复算法计算菲涅耳域双随机相位加密系统输出面上相位的近似分布，设总的迭代次数为某一正整数K，在第k次(k≤K)迭代运算过程中，系统输入面上的图像为fk(x，y)，则其所对应的输出面上的加密结果为：
其中，当k＝1时，初始输入信号f1(x，y)是一元素值均为1的矩阵，由式(2)得到的加密结果Ck(x″，y″)的相位信息可表示为：
qk(x″，y″)＝PR{Ck(x″，y″)} (3)
其中PR{ }表示相位保留操作，即除去复振幅的振幅信息，只保留相位部分的信息，接着，式(3)得到的qk(x″，y″)与密文C(x″，y″)相乘，其乘积作为双随机相位加密系统解密过程中的输入信号，解密结果的振幅分布具体表示为：
其中“*”表示复共轭，IFrT{ }代表菲涅耳变换的逆运算，即逆菲涅耳变换，由φk(x，y)计算得到第k+1次迭代运算中输入面上的图像，即：
fk+1(x，y)＝PT{φk(x，y)} (5)
随后进入下轮迭代过程，即第k+1次迭代过程；
(iii)重复迭代运算过程直至迭代次数大于K时，迭代终止，由式(3)得到qK(x″，y″)，接着对相位qK(x″，y″)进行随机抽取操作，生成稀疏双随机相位加密图像qsp(x″，y″)，即：
qsp(x″，y″)＝SP{qK(x″，y″)} (6)
其中SP{ }表示随机抽取元素操作，即通过保留qK(x″，y″)中被抽取的元素的值，而未被抽取的元素值则以零值代替，从而生成稀疏双随机相位加密图像qsp(x″，y″)；
(2)解密与认证：
(i)稀疏双随机相位加密图像qsp(x″，y″)输入解密与认证系统后经过一次波长为λ，距离为d2的逆菲涅耳变换，变换后得到的结果与相位相乘后再经过一次波长为λ，距离为d1的逆菲涅耳变换，对变换后得到的结果进行取振幅操作得到最终的解密结果，即：
从式(7)可以看出，解密过程需要的密钥包括波长λ，衍射距离d1和d2，以及相位
(ii)将原图f(x，y)与Q(x，y)的乘积与上一步骤中得到的fsp(x，y)进行对比认证，认证采用的非线性相关方法的计算表达式包括：
c(μ，v)＝FT[fsp(x，y)]×{FT[f(x，y)Q(x，y)]}* (8)
NC(x，y)＝|IFT[c(μ，v)|c(μ，v)|ω-1]|2 (9)
其中(μ，v)表示傅立叶频域的坐标，| |表示取模，FT[ ]和IFT[ ]分别表示傅立叶变换和逆傅立叶变换，ω表示非线性的强度，当认证成功时，函数NC(x，y)的分布图将出现尖锐的相关峰。
本发明的有益效果在于：首先，加密过程不需要运用全息技术进行相位的记录，降低了光学双随机相位加密系统的复杂性；其次，输入解密与认证系统的信息是一稀疏相位分布，数据量小，其解密后的结果为噪声图，避免了原图信息的泄露；最后，在认证之前，需要对输入信息进行解密，认证过程中则需要输入QR码，系统安全性高。
【附图说明】
图1为加密过程流程图。
图2为相位恢复算法流程图。
图3为解密与认证过程流程图。
图4为本发明的光电加密装置图。
图5(a)原始图像“Lena”(256×256)，即f(x，y)；(b)用于生成QR码的文字；(c)QR码Q(x，y)；(d)振幅图像C(x″，y″)。
图6(a)原图f(x，y)和K次迭代运算后得到的图像fK+1(x，y)之间的CC值与迭代运算次数K之间的关系图；(b)迭代次数K＝50对应的相位分布q50(x″，y″)。
图7(a)稀疏度为7％的稀疏双随机相位加密图像qsp(x″，y″)的相位分布图；(b)解密结果fsp(x，y)。
图8不同ω对应的解密图像fsp(x，y)与f(x，y)Q(x，y)(即图5(a)所示的原始图像与图5(c)所示的QR码的乘积)进行对比认证的结果：(a)ω＝0.3；(b)ω＝0.4。
图9(a)解密过程中波长参数为λ＝642nm所对应的解密结果；(b)解密过程中衍射距离为d1＝41cm所对应的解密结果；(c)解密过程中衍射距离为d2＝31cm所对应的解密结果；(d)图(a)与f(x，y)Q(x，y)进行对比认证的结果；(e)图(b)与f(x，y)Q(x，y)进行对比认证的结果；(f)图(c)与f(x，y)Q(x，y)进行对比认证的结果。
图10(a)解密过程中解密密钥发生错误时得到的解密结果；(b)图(a)与f(x，y)Q(x，y)进行对比认证的结果。
图11(a)用于生成QR码的文字；(b)对应于图(a)的QR码；(c)f(x，y)[图5(a)]和图(b)所示QR码的乘积与fsp(x，y)[图7(b)]进行对比认证的结果(ω＝0.3)；(d)f(x，y)和fsp(x，y)直接进行对比认证的结果(ω＝0.3)。
图12(a)图像“Cameraman”(256×256)；(b)用于生成对图(a)进行加密的QR码的文字；(c)对应于图(b)的QR码。
图13(a)对应于图像“Cameraman”的稀疏双随机相位加密图像的相位分布(稀疏度为6％)；(b)正确的解密结果；(c)图(a)和图12(c)所示QR码的乘积与图(b)进行对比认证后的结果(ω＝0.3)。
图14(a)图13(b)与f(x，y)Q(x，y)(即图5(a)所示的原始图像与图5(c)所示QR码的乘积)对比认证的结果(ω＝0.3)；(b)f(x，y)[图5(a)]和图12(c)所示QR码的乘积与解密结果[图13(b)]进行对比认证后的结果(ω＝0.3)。
【具体实施方式】
下面本发明结合实施例并参照附图予以详述：本发明所述方法的加密过程可以通过图4所示的光电混合系统实现。空间光调制器(spatial light modulator，SLM)具有调制复振幅信号的能力。加密过程分为两步：
(1)在加密的过程中，利用计算机可控的SLM1和SLM2分别显示信息f(x，y)×Q(x，y)×r1(x，y)和r2(x′，y′)，在波长为λ的单位振幅平面光波的照射下，f(x，y)、Q(x，y)和加密密钥r1(x，y)的乘积作一次距离为d1的菲涅耳变换，变换得到的信号经r2(x′，y′)调制后再作一次距离为d2的菲涅耳变换，输出信号的强度由输出面上放置的光强探测器CCD记录并存储于计算机中，由计算机处理得到f(x，y)经菲涅耳域双随机相位加密系统后输出的振幅图像，即其中PT{ }代表取振幅操作，即除去复振幅的相位信息，只保留振幅部分的信息，FrT{ }代表菲涅耳变换。
(2)利用相位恢复算法计算菲涅耳域双随机相位加密系统输出面上相位的近似分布，设总的迭代次数为某一正整数K，在第k次(k≤K)迭代运算过程中，系统输入面上的图像为fk(x，y)，则其所对应的输出面上的加密结果为：其中，当k＝1时，初始输入信号f1(x，y)是一元素值均为1的矩阵，加密结果Ck(x″，y″)的相位信息可表示为：qk(x″，y″)＝PR{Ck(x″，y″)}，其中PR{ }表示相位保留操作，即除去复振幅的振幅信息，只保留相位部分的信息，接着，qk(x″，y″)与密文C(x″，y″)相乘，其乘积作为双随机相位加密系统解密过程中的输入信号，解密结果的振幅分布具体表示为：其中“*”表示复共轭，IFrT{ }代表菲涅耳变换的逆运算，即逆菲涅耳变换，由φk(x，y)计算得到第k+1次迭代运算中输入面上的图像，即：fk+1(x，y)＝PT{φk(x，y)}，随后进入下轮迭代过程，即第k+1次迭代过程，重复迭代运算过程直至迭代次数大于K时，迭代终止，由式(5)得到qK(x″，y″)，接着对相位qK(x″，y″)进行随机抽取操作，生成稀疏双随机相位加密图像qsp(x″，y″)，即：qsp(x″，y″)＝SP{qK(x″，y″)}，其中SP{ }表示随机抽取元素操作。
本发明提出的方法中运用相关系数(the correlation coefficient，CC)衡量两幅图像的相似度，已知f(x，y)是原始图像，fk+1(x，y)是相位恢复算法中第k次迭代运算得到的加密系统输入面上的输入图像，它们两者间的CC值表示为：
其中E[]代表数学期望运算符，式(10)中函数的坐标已经略去，通过CC值可以反映出本方法中迭代运算过程的收敛性。
本发明提出的安全认证方法，每一幅需要加密的图像都会分配一个不同的QR码，因此在解密与认证过程中，用户除了需要输入波长和衍射距离等这些解密密钥外，还需要输入与原图相对应的QR码。另外，由于稀疏双随机相位加密图像及其解密结果都不会泄露原图的信息，这样就保证了原始图像的安全。
解密与认证过程中，稀疏双随机相位加密图像qsp(x″，y″)输入解密与认证系统后作一次波长为λ，距离为d2的逆菲涅耳变换，变换后得到的结果与相位相乘后再作一次波长为λ，距离为d1的逆菲涅耳变换，对变换后得到的结果进行取振幅操作得到最终的解密结果，即：可以看出，解密过程需要的密钥包括波长λ，衍射距离d1和d2，以及相位将原图f(x，y)与Q(x，y)的乘积与fsp(x，y)进行对比认证，认证采用的非线性相关方法的计算表达式包括：c(μ，v)＝FT[fsp(x，y)]×{FT[f(x，y)Q(x，y)]}*，NC(x，y)＝|IFT[c(μ，v)|c(μ，v)|ω-1]|2，其中(μ，v)表示傅立叶频域的坐标，| |表示取模，FT[ ]和IFT[ ]分别表示傅立叶变换和逆傅立叶变换，ω表示非线性的强度，当认证成功时，函数NC(x，y)的分布图将出现尖锐的相关峰，从中可以看出，Q(x，y)是认证所需的密钥。
下面结合实施例和附图对本发明的内容进行进一步的解释。
选择如图5(a)所示大小为256×256的归一化灰度图“Lena”作为待加密图像和认证中的原图，图5(b)显示的是用来生成QR码的文字，即“Zhejiang A&FUniversity”，通过二维码软件加密后生成的QR码则如图5(c)所示。仿真中采用的入射光波波长为λ＝632nm，衍射距离分别为d1＝40cm和d2＝30cm，根据加密流程图图1所示，在图5(c)所示的QR码和计算机生成的两块随机相位板的作用下，图像“Lena”最终被加密成振幅图像C(x″，y″)，其分布如图5(d)所示。
接着，根据相位恢复算法流程图图2所示，通过相位恢复算法计算得到相位分布qK(x″，y″)。图6(a)给出了迭代运算过程中原图f(x，y)和K次迭代运算后得到的图像fK+1(x，y)之间的CC值与不同的迭代运算次数K之间的关系图，从中可以看出，该相位恢复算法具有收敛速度快的特点。在迭代运算5次之后，CC值就达到了0.9865，迭代次数K为10时，CC值为0.9953，当迭代次数K的大小超过40时，CC值达到理论最大值，即1。迭代运算次数为K＝50时计算得到的相位分布q50(x″，y″)则如图6(b)所示。
根据qsp(x″，y″)＝SP{qK(x″，y″)}，其中SP{ }表示随机抽取元素操作，通过保留q50(x″，y″)中被抽取的元素的值，而未被抽取的元素值则以零值代替，从而生成稀疏双随机相位加密图像qsp(x″，y″)。将稀疏双随机相位加密图像中非零元素占所有元素的比例定义为稀疏度，则图7(a)给出了稀疏度为7％的稀疏双随机相位加密图像qsp(x″，y″)的相位分布图。
解密与认证的过程参照图3，稀疏双随机相位加密图像qsp(x″，y″)进行一次波长为λ，距离为d2的逆菲涅耳变换，变换得到的结果与解密密钥相乘后再作一次波长为λ，距离为d1的逆菲涅耳变换，对变换后得到的结果进行取振幅操作，则有图7(b)是对应于图7(a)所示的qsp(x″，y″)的解密图像，即fsp(x，y)。从图7可以看出，稀疏双随机相位加密图像及其对应的解密结果是充满噪声的乱码，无法从中观察到原图f(x，y)的信息。接着，将原图f(x，y)和Q(x，y)的乘积与fsp(x，y)进行对比认证，对比过程中非线性强度ω取值的不同将会影响函数NC(x，y)值的大小，从而影响相关峰的高度和尖锐程度，但不会改变认证是否获得通过这一最终结果。非线性的强度ω的取值为0.3和0.4所对应的函数NC(x，y)的分布图分别如图8(a)和8(b)所示，两幅图中均出现了相关峰，这表明在使用不同的非线性强度情况下，信息的认证均获得通过。
下面考察解密与认证过程中光学参数的作用。当对图7(a)进行解密时使用的光学参数(即波长λ、衍射距离d1以及衍射距离d2)中的其中一个发生错误时，都将得到与图7(b)不同的解密结果。解密中若采用的波长为λ＝642nm，得到的解密结果则如图9(a)所示。解密时使用的衍射距离为d1为41cm时得到的解密结果如图9(b)所示，当输入的衍射距离d2为31cm时所对应的解密结果如图9(c)所示。将这些解密结果依次与原图f(x，y)和QR码Q(x，y)的乘积进行对比认证，认证结果分别如图9(d)-9(f)所示。接着考察解密密钥的作用，当发生错误时，得到的解密结果也将与图7(b)不同，因此认证的结果也将不同。仿真中用一计算机生成的随机相位分布代替进行解密，得到的解密结果如图10(a)所示，将原图f(x，y)和Q(x，y)的乘积与图10(a)进行对比认证后，得到的结果如图10(b)所示。从图9和图10可以看出，图像的解密过程中有一个密钥发生错误时，即使其余密钥全部正确，信息都将无法通过认证，即在认证结果中不会出现尖锐的相关峰。
接着考察认证密钥即QR码在认证过程中的作用。图11(a)是用来生成不同于图5(c)所示QR码的文字，即“double random phase encoding”，生成的QR码如图11(b)所示。将图11(b)所示的QR码用于图7(b)所示解密图像fsp(x，y)的认证，即原图f(x，y)和图11(b)所示的QR码相乘后再与图像fsp(x，y)进行对比认证，对比过程中非线性的强度ω的取值为0.3，得到的认证结果如图11(c)所示。当认证过程中不使用任何QR码，即原图f(x，y)和图7(b)所示的图像fsp(x，y)直接进行对比认证(对比过程中非线性的强度ω的取值为0.3)，得到的认证结果如图11(d)所示。从图11(c)和图11(d)可以看出，认证过程中若使用错误的QR码或者不使用QR码，认证结果中都不会出现尖锐的相关峰，即信息无法通过认证。
按照加密流程图图1所示，对图12(a)所示的归一化灰度图像“Cameraman”进行加密，加密过程使用的QR码如图12(c)所示，图12(b)是用来生成图12(c)所示QR码的文字，即“information authentication”。加密过程最终生成的稀疏度为6％的稀疏双随机相位加密图像的相位分布如图13(a)所示。根据图3所示的解密与认证流程得到的解密图像如图13(b)所示。接着，图像“Cameraman”和图12(c)所示的QR相乘后与解密结果图13(b)进行对比认证，当ω取值为0.3时，认证结果如图13(c)所示，图中出现的相关峰表明在所有解密密钥和认证密钥都正确的情况下，解密与认证系统的输入信息即图13(a)顺利通过了认证。
下面考察本发明提出的认证方法的识别能力。将图5(a)所示的图像“Lena”作为认证过程的原图，图13(a)所示的稀疏双随机相位加密图像(对应于图像“Cameraman”)和图5(c)所示的QR码(对应于图像“Lena”)分别作为解密与认证系统的输入信息和认证密钥，由此得到的认证结果如图14(a)所示。仍将图5(a)所示的图像“Lena”作为认证过程使用的原图，而图13(a)所示的稀疏双随机相位加密图像(对应于图像“Cameraman”)和图12(b)所所示的QR码(对应于图像“Cameraman”)分别作为解密与认证系统的输入信息和认证密钥，由此得到的认证结果则如图14(b)所示，以上认证过程中ω取值均为0.3。从图14可以看出，使用不同于原图的图像所生成的稀疏双随机相位加密图像进行解密与认证时，不论输入的认证密钥是否正确，认证得到的图像中都不会出现尖锐的相关峰，即该信息无法通过认证。

Claims (1)

1.一种基于稀疏双随机相位加密和QR码的认证方法，其特征是按如下步骤进行：
(1)加密：
(i)Q(x，y)代表光学加密时作为加密密钥和认证时作为认证密钥的QR码，f(x，y)代表待加密并用于认证的原始图像，利用菲涅耳域双随机相位加密系统对f(x，y)进行加密，r1(x，y)和r2(x′，y′)是双随机相位加密系统中的两块随机相位板，分别表示成exp[2πa(x，y)]和exp[2πb(x′，y′)]，其中a(x，y)、b(x′，y′)代表两个统计无关且在区间[0，1]上具有均匀概率分布的随机矩阵，(x，y)和(x′，y′)分别表示输入平面和菲涅耳衍射平面的坐标，加密时图像f(x，y)和加密密钥Q(x，y)相乘，它们的乘积与r1(x，y)相乘后作一次波长为λ，距离为d1的菲涅耳变换，得到的结果与r2(x′，y′)相乘后再作一次波长为λ，距离为d2的菲涅耳变换，对变换后的结果取振幅后得到密文C(x″，y″)，即：
<mrow> <mi>C</mi> <mrow> <mo>(</mo> <mrow> <msup> <mi>x</mi> <mrow> <mo>&amp;prime;</mo> <mo>&amp;prime;</mo> </mrow> </msup> <mo>,</mo> <msup> <mi>y</mi> <mrow> <mo>&amp;prime;</mo> <mo>&amp;prime;</mo> </mrow> </msup> </mrow> <mo>)</mo> </mrow> <mo>=</mo> <mi>P</mi> <mi>T</mi> <mrow> <mo>{</mo> <mrow> <msub> <mi>FrT</mi> <mrow> <msub> <mi>d</mi> <mn>2</mn> </msub> <mo>,</mo> <mi>&amp;lambda;</mi> </mrow> </msub> <mrow> <mo>{</mo> <mrow> <msub> <mi>FrT</mi> <mrow> <msub> <mi>d</mi> <mn>1</mn> </msub> <mo>,</mo> <mi>&amp;lambda;</mi> </mrow> </msub> <mrow> <mo>{</mo> <mrow> <mi>f</mi> <mrow> <mo>(</mo> <mrow> <mi>x</mi> <mo>,</mo> <mi>y</mi> </mrow> <mo>)</mo> </mrow> <mi>Q</mi> <mrow> <mo>(</mo> <mrow> <mi>x</mi> <mo>,</mo> <mi>y</mi> </mrow> <mo>)</mo> </mrow> <msub> <mi>r</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mrow> <mi>x</mi> <mo>,</mo> <mi>y</mi> </mrow> <mo>)</mo> </mrow> </mrow> <mo>}</mo> </mrow> <mo>&amp;times;</mo> <msub> <mi>r</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <mrow> <msup> <mi>x</mi> <mo>&amp;prime;</mo> </msup> <mo>,</mo> <msup> <mi>y</mi> <mo>&amp;prime;</mo> </msup> </mrow> <mo>)</mo> </mrow> </mrow> <mo>}</mo> </mrow> </mrow> <mo>}</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
其中PT{}代表取振幅操作，即除去复振幅的相位信息，只保留振幅部分的信息，FrT{}代表菲涅耳变换，(x″，y″)表示菲涅耳域双随机相位加密系统输出平面的坐标；
(ii)利用相位恢复算法计算菲涅耳域双随机相位加密系统输出面上相位的近似分布，设总的迭代次数为某一正整数K，在第k次(k≤K)迭代运算过程中，系统输入面上的图像为fk(x，y)，则其所对应的输出面上的加密结果为：
<mrow> <msub> <mi>C</mi> <mi>k</mi> </msub> <mrow> <mo>(</mo> <mrow> <msup> <mi>x</mi> <mrow> <mo>&amp;prime;</mo> <mo>&amp;prime;</mo> </mrow> </msup> <mo>,</mo> <msup> <mi>y</mi> <mrow> <mo>&amp;prime;</mo> <mo>&amp;prime;</mo> </mrow> </msup> </mrow> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>FrT</mi> <mrow> <msub> <mi>d</mi> <mn>2</mn> </msub> <mo>,</mo> <mi>&amp;lambda;</mi> </mrow> </msub> <mrow> <mo>{</mo> <mrow> <msub> <mi>FrT</mi> <mrow> <msub> <mi>d</mi> <mn>1</mn> </msub> <mo>,</mo> <mi>&amp;lambda;</mi> </mrow> </msub> <mrow> <mo>{</mo> <mrow> <msub> <mi>f</mi> <mi>k</mi> </msub> <mrow> <mo>(</mo> <mrow> <mi>x</mi> <mo>,</mo> <mi>y</mi> </mrow> <mo>)</mo> </mrow> <mi>Q</mi> <mrow> <mo>(</mo> <mrow> <mi>x</mi> <mo>,</mo> <mi>y</mi> </mrow> <mo>)</mo> </mrow> <msub> <mi>r</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mrow> <mi>x</mi> <mo>,</mo> <mi>y</mi> </mrow> <mo>)</mo> </mrow> </mrow> <mo>}</mo> </mrow> <mo>&amp;times;</mo> <msub> <mi>r</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <mrow> <msup> <mi>x</mi> <mo>&amp;prime;</mo> </msup> <mo>,</mo> <msup> <mi>y</mi> <mo>&amp;prime;</mo> </msup> </mrow> <mo>)</mo> </mrow> </mrow> <mo>}</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>
其中，当k＝1时，初始输入信号f1(x，y)是一元素值均为1的矩阵，由式(2)得到的加密结果Ck(x″，y″)的相位信息可表示为：
qk(x″，y″)＝PR{Ck(x″，y″)} (3)
其中PR{}表示相位保留操作，即除去复振幅的振幅信息，只保留相位部分的信息，接着，式(3)得到的qk(x″，y″)与密文C(x″，y″)相乘，其乘积作为双随机相位加密系统解密过程中的输入信号，解密结果的振幅分布具体表示为：
<mrow> <msub> <mi>&amp;phi;</mi> <mi>k</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>P</mi> <mi>T</mi> <mo>{</mo> <msub> <mi>IFrT</mi> <mrow> <msub> <mi>d</mi> <mn>1</mn> </msub> <mo>,</mo> <mi>&amp;lambda;</mi> </mrow> </msub> <mo>{</mo> <msub> <mi>IFrT</mi> <mrow> <msub> <mi>d</mi> <mn>2</mn> </msub> <mo>,</mo> <mi>&amp;lambda;</mi> </mrow> </msub> <mo>{</mo> <mi>C</mi> <mrow> <mo>(</mo> <msup> <mi>x</mi> <mrow> <mo>&amp;prime;</mo> <mo>&amp;prime;</mo> </mrow> </msup> <mo>,</mo> <msup> <mi>y</mi> <mrow> <mo>&amp;prime;</mo> <mo>&amp;prime;</mo> </mrow> </msup> <mo>)</mo> </mrow> <msub> <mi>q</mi> <mi>k</mi> </msub> <mrow> <mo>(</mo> <msup> <mi>x</mi> <mrow> <mo>&amp;prime;</mo> <mo>&amp;prime;</mo> </mrow> </msup> <mo>,</mo> <msup> <mi>y</mi> <mrow> <mo>&amp;prime;</mo> <mo>&amp;prime;</mo> </mrow> </msup> <mo>)</mo> </mrow> <mo>}</mo> <mo>&amp;times;</mo> <msubsup> <mi>R</mi> <mn>2</mn> <mo>*</mo> </msubsup> <mrow> <mo>(</mo> <msup> <mi>x</mi> <mo>&amp;prime;</mo> </msup> <mo>,</mo> <msup> <mi>y</mi> <mo>&amp;prime;</mo> </msup> <mo>)</mo> </mrow> <mo>}</mo> <mo>}</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>
其中“*”表示复共轭，IFrT{}代表菲涅耳变换的逆运算，即逆菲涅耳变换，由φk(x，y)计算得到第k+1次迭代运算中输入面上的图像，即：
fk+1(x，y)＝PT{φk(x，y)} (5)
随后进入下轮迭代过程，即第k+1次迭代过程；
(iii)重复迭代运算过程直至迭代次数大于K时，迭代终止，由式(3)得到qK(x″，y″)，接着对相位qK(x″，y″)进行随机抽取操作，生成稀疏双随机相位加密图像qsp(x″，y″)，即：
qsp(x″，y″)＝SP{qK(x″，y″)} (6)
其中SP{}表示随机抽取元素操作，即通过保留qK(x″，y″)中被抽取的元素的值，而未被抽取的元素值则以零值代替，从而生成稀疏双随机相位加密图像qsp(x″，y″)；
(2)解密与认证：
(i)稀疏双随机相位加密图像qsp(x″，y″)输入解密与认证系统后经过一次波长为λ，距离为d2的逆菲涅耳变换，变换后得到的结果与相位相乘后再经过一次波长为λ，距离为d1的逆菲涅耳变换，对变换后得到的结果进行取振幅操作得到最终的解密结果，即：
<mrow> <msub> <mi>f</mi> <mrow> <mi>s</mi> <mi>p</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>P</mi> <mi>T</mi> <mo>{</mo> <msub> <mi>IFrT</mi> <mrow> <msub> <mi>d</mi> <mn>1</mn> </msub> <mo>,</mo> <mi>&amp;lambda;</mi> </mrow> </msub> <mo>{</mo> <msub> <mi>IFrT</mi> <mrow> <msub> <mi>d</mi> <mn>2</mn> </msub> <mo>,</mo> <mi>&amp;lambda;</mi> </mrow> </msub> <mo>{</mo> <msub> <mi>q</mi> <mrow> <mi>s</mi> <mi>p</mi> </mrow> </msub> <mrow> <mo>(</mo> <msup> <mi>x</mi> <mrow> <mo>&amp;prime;</mo> <mo>&amp;prime;</mo> </mrow> </msup> <mo>,</mo> <msup> <mi>y</mi> <mrow> <mo>&amp;prime;</mo> <mo>&amp;prime;</mo> </mrow> </msup> <mo>)</mo> </mrow> <mo>}</mo> <mo>&amp;times;</mo> <msubsup> <mi>r</mi> <mn>2</mn> <mo>*</mo> </msubsup> <mrow> <mo>(</mo> <msup> <mi>x</mi> <mo>&amp;prime;</mo> </msup> <mo>,</mo> <msup> <mi>y</mi> <mo>&amp;prime;</mo> </msup> <mo>)</mo> </mrow> <mo>}</mo> <mo>}</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>
从式(7)可以看出，解密过程需要的密钥包括波长λ，衍射距离d1和d2，以及相位
(ii)将原图f(x，y)与Q(x，y)的乘积与上一步骤中得到的fsp(x，y)进行对比认证，认证采用的非线性相关方法的计算表达式包括：
c(μ，v)＝FT[fsp(x，y)]×{FT[f(x，y)Q(x，y)]}* (8)
NC(x，y)＝|IFT[c(μ，v)|c(μ，v)|ω-1]|2 (9)
其中(μ，v)表示傅立叶频域的坐标，||表示取模，FT[]和IFT[]分别表示傅立叶变换和逆傅立叶变换，ω表示非线性的强度，当认证成功时，函数NC(x，y)的分布图将出现尖锐的相关峰。
CN201510050296.2A 2015-01-30 2015-01-30 基于稀疏双随机相位加密和qr码的认证方法 Expired - Fee Related CN104573561B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510050296.2A CN104573561B (zh) 2015-01-30 2015-01-30 基于稀疏双随机相位加密和qr码的认证方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510050296.2A CN104573561B (zh) 2015-01-30 2015-01-30 基于稀疏双随机相位加密和qr码的认证方法

Publications (2)

Publication Number Publication Date
CN104573561A CN104573561A (zh) 2015-04-29
CN104573561B true CN104573561B (zh) 2017-12-12

Family

ID=53089594

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510050296.2A Expired - Fee Related CN104573561B (zh) 2015-01-30 2015-01-30 基于稀疏双随机相位加密和qr码的认证方法

Country Status (1)

Country Link
CN (1) CN104573561B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106533667B (zh) * 2016-11-08 2019-07-19 深圳大学 基于双光束干涉的多级密钥生成方法和用户分级认证方法
CN107104790A (zh) * 2017-03-20 2017-08-29 上海师范大学 基于qr码和阵列指纹密钥的光学加密方法
CN112765623B (zh) * 2021-01-15 2022-08-02 浙江科技学院 基于相位恢复算法的光学多图像认证和加密方法
CN112765624B (zh) * 2021-01-19 2022-05-27 浙江科技学院 基于相位优化和稀疏约束的可认证纯相位全息图生成方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103279916A (zh) * 2013-05-24 2013-09-04 浙江农林大学 基于双随机相位编码和干涉原理的图像加密方法
CN104159073A (zh) * 2014-07-21 2014-11-19 中国人民武装警察部队工程大学 一种基于双随机相位编码的视频光学加解密系统及方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103279916A (zh) * 2013-05-24 2013-09-04 浙江农林大学 基于双随机相位编码和干涉原理的图像加密方法
CN104159073A (zh) * 2014-07-21 2014-11-19 中国人民武装警察部队工程大学 一种基于双随机相位编码的视频光学加解密系统及方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
A special attack on the asymmetric cryptosistem based on phase-truncated Fourier transforms;Xiaogang Wang 等;《Optics Communications》;20121231;第1078至1081页 *
一种基于迭代振幅-相位恢复算法和非线性双随机相位编码的图像加密方法;陈翼翔 等;《光学学报》;20140831;第34卷(第8期);第0810003-1至0810003-6页 *
利用QR码在光学干涉多图像加密系统中实现信息高质量恢复;王红娟 等;《光学学报》;20140930;第34卷(第9期);第0907001-1至0907001-7页 *
基于双随机相位编码的非线性双图像加密方法;陈翼翔 等;《光学学报》;20140731;第34卷(第7期);第0710001-1至0710001-5页 *
基于菲涅耳域光学图像加密系统的解密研究;肖永亮 等;《激光技术》;20090831;第33卷(第4期);第433至436页 *

Also Published As

Publication number Publication date
CN104573561A (zh) 2015-04-29

Similar Documents

Publication Publication Date Title
CN104517261B (zh) 基于相位恢复和稀疏双随机相位加密的安全认证方法
Verma et al. An optical asymmetric encryption scheme with biometric keys
Li et al. Asymmetric multiple-image encryption based on the cascaded fractional Fourier transform
CN102567943B (zh) 基于联合傅立叶变换和相位切除的非对称双图像加密方法
CN103279916B (zh) 基于双随机相位编码和干涉原理的图像加密方法
Singh Watermarking image encryption using deterministic phase mask and singular value decomposition in fractional Mellin transform domain
CN104573561B (zh) 基于稀疏双随机相位加密和qr码的认证方法
EP2089794B1 (en) A method of generating arbitrary numbers given a seed
Wang et al. Multiple-image encryption and authentication in interference-based scheme by aid of space multiplexing
Wei et al. Optical image encryption using QR code and multilevel fingerprints in gyrator transform domains
US20200177380A1 (en) Seed key expansion method and its uses
Musanna et al. A novel chaos-based approach in conjunction with MR-SVD and pairing function for generating visually meaningful cipher images
CN103023633A (zh) 基于混沌随机相位和相干叠加原理的数字图像隐藏方法
CN103258315B (zh) 基于切相分数傅立叶变换的双图像加密方法
Lang et al. The discrete multiple-parameter fractional Fourier transform
Kumar et al. Multiuser optical image authentication platform based on sparse constraint and polar decomposition in Fresnel domain
CN104376524B (zh) 基于光阑加密和相位恢复算法的二值图像加密方法
Farwa et al. An efficient image encryption scheme using Fresnelet transform and elliptic curve based scrambling
Chen et al. Optical information authentication using optical encryption and sparsity constraint
CN104134184B (zh) 基于迭代切相分数傅立叶变换的图像加密方法
Mohan et al. Secure visual cryptography scheme with meaningful shares
Liansheng et al. Multiple-image authentication based on the single-pixel correlated imaging and multiple-level wavelet transform
Chen et al. Information authentication using sparse representation of double random phase encoding in fractional Fourier transform domain
CN104376525B (zh) 基于迭代非线性双随机相位编码的图像加密方法
Liu et al. Optical asymmetric JTC cryptosystem based on multiplication-division operation and RSA algorithm

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20171212

Termination date: 20190130

CF01 Termination of patent right due to non-payment of annual fee