一种量子逆向传播神经网络的数字全息压缩传输方法
【专利摘要】本发明公布了一种量子逆向传播神经网络的数字全息压缩传输方法，将原始图像通过菲涅尔离轴全息计算方法记录为数字全息图、构建量子逆向传播(QBP)神经网络、使用数字全息图对QBP神经网络进行训练、使用训练完毕的QBP神经网络对数字全息图进行压缩、传输和解压缩，得到重构后的全息图，以及对重构后的全息图进行再现，得到再现像。本发明提供方法具有更快的并行处理速度和更强的存储数据的能力，适用于大数据量的数字全息图的计算处理；能够使用更少的训练次数完成压缩传输网络结构的训练，提高数字全息图的压缩传输速度；图像压缩率可调，再现图像质量较好。
【专利说明】
一种量子逆向传播神经网络的数字全息压缩传输方法
技术领域
[0001] 本发明属于数字图像处理领域，涉及一种采用量子逆向传播（Quantum Back-Pr〇pagati〇n，QBP)神经网络的数字全息压缩传输方法，具体涉及数字全息图像的制作与再 现、QBP神经网络图像压缩等方法。
【背景技术】
[0002] 1947年英国物理学家丹尼斯?伽博(Dennis Gabor)在研究如何提高电子显微镜 的分辨率的过程中发明了全息术。全息术是一种能够记录和再现真实世界的三维(Three-Dimensional ，3D) 物体的技术。物体反射的光波与参考光波相干叠加产生干涉条纹，被记录 的这些干涉条纹称为全息图。全息图在一定的条件下再现，便可重现原物体逼真的三维像。 1967年美国科学家顾德门（J.W.Goodman)首次提出了数字全息。数字全息是一种全新的获 取光学信息的方法，它是传统的全息术和数字技术相结合的产物。随着计算机技术和高分 辨率图像传感器的飞速发展，数字全息技术的优势正在越来越明显地显示出来，其应用范 围已涉及三维形貌测量、形变测量、粒子场测试、显微和防伪等许多领域。
[0003] 数字全息图记录了物体的振幅与相位信息，数据量非常庞大，为了进行传输必须 进行压缩处理。数字全息图的信息是一种具有高度非线性强度分布的干涉条纹图案，这种 强度分布包含着物体的信息。人工神经网络已经被广泛的应用在非线性领域，它可以自适 应地和智能地调整来处理不同的图像信息分布，因此非常适合处理全息信息分布问题。据 此，张超等人在文南犬一 "Chao Zhang ,Guangl in Yang and Haiyan Xie , ^Information Compression of Computer-Generated Hologram Using BP Neural Network,〃in Biomedical Optics and 3-D Imaging，0SA Technical Digest(CD)(Optical Society of America)，paper JMA2,2010."中第一次验证了使用人工神经网络传输全息信息的可行性， 提出了使用传统逆向传播(Back-Pr〇pagati 〇n，BP)神经网络进行数字全息图压缩传输的技 术，并且使用计算机仿真得到了实验结果。然而，由于传统BP神经网络具有在信息量大的情 况下处理速度过慢、记忆容量有限等缺陷，这些缺陷限制了这种技术应用于数字全息图像 压缩传输的效果与应用价值。
【发明内容】
[0004] 为了克服上述现有技术的不足，本发明提供一种量子逆向传播神经网络的数字全 息压缩传输方法，基于QBP神经网络对数字全息图进行压缩传输，能够使用更少的训练次 数完成压缩传输网络结构的训练，提高数字全息图的压缩传输速度，同时保证了图像的恢 复质量。
[0005] 本发明的原理是:本发明将QBP神经网络技术用于数字全息图的压缩传输，由于数 字全息图是干涉条纹图像，每一像素都携带了其它像素的部分信息，数据量很大;而量子理 论的态叠加原理使QBP神经网络具有比传统BP神经网络更快的并行处理速度和更强的存储 数据的能力，因而非常适合应用于大数据量的数字全息图的计算处理。根据量子计算原理， 一个η位量子寄存器可以同时保存2"个11位二进制数(0到2n-l)，它们各以一定的概率存在。 量子计算系统以这种方式指数地增加存储能力并能并行处理一个η位量子寄存器的所有2n 个数，它的一次运算可产生2n个运算结果，相当于常规计算2n次操作。因此，通过将QBP神经 网络应用于数字全息信息处理，可以得到比使用传统BP神经网络更好的效果。
[0006] 本发明提供的技术方案是：
[0007] -种量子逆向传播神经网络的数字全息压缩传输方法，包括数字全息图像的制作 与再现、量子逆向传播(QBP)神经网络图像压缩等过程，具有更快的并行处理速度和更强的 存储数据的能力，迭代次数少，压缩传输速度快，适用于大数据量的数字全息图的计算处 理;能够使用更少的训练次数完成压缩传输网络结构的训练，提高数字全息图的压缩传输 速度;图像压缩率可调，再现图像质量较好;包括如下步骤：
[0008] 1)将原始图像通过菲涅尔离轴全息计算方法记录为数字全息图；
[0009] 2)构建量子逆向传播神经网络，所述量子逆向传播神经网络包含一个输入层、一 个隐含层和一个输出层，所述量子逆向传播神经网络每一层包含多个量子神经元;使用步 骤1)得到的数字全息图对所述量子逆向传播神经网络进行训练，得到训练完毕的量子逆向 传播神经网络；
[0010] 3)使用步骤2)得到的训练完毕的量子逆向传播神经网络对步骤1)得到的数字全 息图进行压缩、传输和解压缩，得到重构后的全息图；
[0011] 4)对步骤3)得到的重构后的全息图进行再现，得到再现像。
[0012] 针对上述量子逆向传播神经网络的数字全息压缩传输方法，进一步地，步骤1)所 述记录包括菲涅尔衍射计算方法和离轴参考光干涉计算方法；
[0013] 所述菲涅尔衍射计算方法由菲涅耳衍射公式计算，所述菲涅耳衍射公式为式2:
[0014]
[0015] 式2中，d为衍射距离;k为波数;j为虚数单位；（XQ，yQ)为初始物波平面；（x，y)为衍 射平面;〇( XQ，yQ)为初始物波平面信息，U(x，y)为经过衍射距离为d的物光波前复振幅函数；
[0016] 所述离轴参考光干涉计算方法利用干涉记录物光衍射波前复振幅函数形成全息 图 lH(x，y)，如式 1:
[0017]
[0018]式1中，lH(x，y)为全息图记录的光强;R(x，y)为参考光波前函数;U(x，y)为原始物 体波前函数;Rlxd)与υ?χ,γ)分别为R(x，y)与U(x，y)的共辄项。
[0019]针对上述量子逆向传播神经网络的数字全息压缩传输方法，更进一步地，通过式3 进行快速傅里叶变换，对菲涅尔衍射进行计算：
[0020]
[0021] 式3中，F表示傅立叶变换；d为衍射距离;k为波数；j为虚数单位;Lo为初始物波平 面采样范围；初始物波平面的采样点数为NXN，即初始图像的像素个数;初始物波平面采样 点的间隔为Δχ〇= Ay〇 = L()/N; Δχ= Ay是初始物波平面离散傅立叶变换后采样间隔，在空 间域内，与之对应的为衍射平面的采样间隔。
[0022] 针对上述量子逆向传播神经网络的数字全息压缩传输方法，进一步地，步骤2)具 体包括如下步骤：
[0023] 21)基于量子门建立新的量子神经元模型，用于实现量子计算;所述新的量子神经 元模型中，权值与输入相乘的结果通过相移门对神经元态相移实现：
[0024] 设]^(1 = 1，2，>"，]\〇为第;[个输入到神经元的量子态11 = ;^10(1 = 1，2，>"，]\〇的相 位，则该量子神经元输出表示为如下量子神经元模型：
[0025]
[0026]
[0027] 〇 = f(y)(式 9)
[0028] 其中，0,(1 = 1，2，···，Μ)为权值的相位转移系数;λ为阈值系数;δ为相位控制因子； 〇为输出态;3找(11)是对复数11提取相位，8卩3坪(11)=3(^38(1111(11)/1^(11))，其中，1111(11)为求 复数U的虚部;Re(u)为求复数u的实部;g(5)是sigmoid函数;函数f( ·)的定义如式5所示；
[0029] f (Θ) =e10 = cos9+isin9 (式 5)
[0030] 其中，:i .=：_λΡ?表示虚数单元;Θ表示量子态的相位；
[0031] 22)构建用于图像压缩传输的量子逆向传播神经网络，所述量子逆向传播神经网 络每一层的量子神经兀的输入和输出是多个量子态的叠加；
[0032] 23)定义神经网络的逆向传播训练方法，训练步骤22)构建的量子逆向传播神经网 络;所述神经网络的逆向传播训练方法采用近似均方误差的最速下降算法，调整所述量子 逆向传播神经网络中相位转动系数Θ、阈值系数λ和相位控制因子δ，使训练误差的均值小于 期望目标时结束训练，得到训练完毕的量子逆向传播神经网络。
[0033] 针对上述量子逆向传播神经网络的数字全息压缩传输方法，更进一步地，步骤21) 所述量子门由一位相移门和两位受控非门组成，作为神经网络的激活函数，用复数形式表 示量子态及通用量子逻辑门组。
[0034] 步骤22)中，所述量子神经元的输入和输出具体是：
[0035] 设定{INl}(l = l，2，…，L)、{HIDEk}(k=l，2，…，K)和{0UTn}(n = l，2，…，N)分别表 示所述量子逆向传播神经网络的输入层、隐含层和输出层的量子神经元集合;L，K和Ν分别 对应各层的神经元个数;当归一化的输入值inputi(1 = 1，2，…，L)被输入到输入层{ΙΝι}，输 入层各神经元将使介于[0，1]之间的输入值转换为介于[0，V2]之间量子态xi(l = l，2,…， L)的相位值Ii(l = 1，2，…，L)，然后将其量子态输出至隐含层;表达式为：
[0036]
[0037] xi = f (Ιι)(式 11)
[0038] 其中，函数f( ·)的定义如式5所示；
[0039] 所述隐含层{HIDEk}和所述输出层{0UTn}的处理过程按照式7~式9;
[0040] 定义任意一个量子神经元的量子态为激活态和抑制态的叠加态，所述神经元的最 终输出值putputn (η=1，2，···，Ν)为激活态11 >时的概率;所述量子逆向传播神经网络输出 层第η个神经元的最终输出值表示为式12:
[0041] outputn= | Im(On) |2 (式12)
[0042] 式12中，0为输出态。
[0043] 针对上述量子逆向传播神经网络的数字全息压缩传输方法，更进一步地，步骤23) 中，所述神经网络的逆向传播训练方法使用步骤1)得到的数字全息图对所述量子逆向传播 神经网络进行训练，具体包括如下步骤：
[0044] 231)初始化所述量子逆向传播神经网络，将步骤1)得到的大小为XX Y的数字全息 图归一化处理，然后划分为xP X yi^图像块，每个图像块变为Μ X 1 (M = xP X yP)维的向量，作 为所述量子逆向传播神经网络的训练样本；
[0045] 232)将一个训练样本输入到所述量子逆向传播神经网络中计算网络输出；
[0046] 233)记录网络实际输出与目标输出的误差，根据最速下降法，逆向逐层调整各层 参数；
[0047] 234)重复步骤232)~233)，直至所有训练样本都已输入；
[0048] 235)设定训练次数，设定累计总误差为输入数字全息图与网络处理后输出的数字 全息图之间的误差;计算所有训练样本的实际输出与目标输出的累计总误差，训练次数加 1;当所述累计总误差小于设定的累计总误差或所述训练次数大于设定的训练次数时，训练 结束;否则重新进入步骤232)。
[0049] 针对上述量子逆向传播神经网络的数字全息压缩传输方法，进一步地，步骤3)具 体包括:将步骤1)得到的数字全息图归一化处理后输入到所述训练完毕的量子逆向传播神 经网络的输入层，通过输入层到隐含层实现图像的编码/压缩，通过隐含层到输出层实现图 像的解码/解压缩，得到重构后的归一化全息图；将所述重构后的归一化全息图再经过反归 一化处理，得到重构后的全息图。
[0050] 针对上述量子逆向传播神经网络的数字全息压缩传输方法，更进一步地，所述量 子逆向传播神经网络的输入层和输出层的量子神经元数量相等，记作L = N;每个神经元对 应一个像素;所述隐含层的量子神经元数量记为K，K小于N;通过调整所述隐含层的神经元 数量Κ的大小，达到调节图像压缩率的目的。
[0051] 针对上述量子逆向传播神经网络的数字全息压缩传输方法，进一步地，步骤4)所 述再现具体通过离轴参考光照射过程和菲涅尔衍射过程实现菲涅尔离轴全息图再现。 [0052]与现有技术相比，本发明的有益效果是：
[0053]本发明提供一种量子逆向传播神经网络的数字全息压缩传输方法，使用QBP神经 网络针对数字全息图进行压缩传输。
[0054]本发明的优点主要体现在以下几方面：
[0055](一)本发明所使用的QBP神经网络基于其量子理论的态叠加原理，具有比传统BP 神经网络更快的并行处理速度和更强的存储数据的能力，非常适合应用于大数据量的数字 全息图的计算处理。
[0056](二)对于一幅相同的数字全息图，训练QBP神经网络所需的训练次数远小于传统 BP神经网络，缩短了计算时间，提高了传输全息信息的速度；
[0057](三)本发明将QBP神经网络算法引入到数字全息信息的压缩传输中，图像压缩率 可调，并且可以保证较好的再现图像质量。
【附图说明】
[0058] 图1是本发明提供的量子逆向传播神经网络的数字全息压缩传输方法的流程框 图。
[0059] 图2是量子神经元模型的示意图。
[0060] 图3是三层QBP神经网络的示意图。
[0061]图4是QBP网络学习和训练计算的流程图。
[0062]图5为原始lena图像的数字全息图及其再现像；
[0063]其中，（a)为原始lena图像（128 X 128像素）；（b)为lena图像的数字全息图（256 X 256像素）；（c)为lena图像的再现像。
[0064]图6为QBP和BP神经网络在四个不同压缩率下的训练次数图；
[0065]其中，（&)(13)对应1? = 0.5000 3^狀=0.0005的情况；（(3)((1)对应1? = 0.2500， Errmax = 0.0015的情况；（e)(f)对应R = 0.1250，Err臆= 0.0023的情况；（g)(h)对应R = 0 · 0625，Errmax = 0 · 0030 的情况。
[0066] 图7为在最佳学习速率_ρ = 0.40和t1bp = 0.20时，得到的QBP和BP的处理后的数字 全息图和再现图像的对比图；
[0067]其中，图7(a)(b)为基于BP神经网络处理后的数字全息图和再现图像；图7(c)(d) 为基于QBP神经网络处理后的数字全息图和再现图像。
【具体实施方式】
[0068] 下面结合附图，通过实施例进一步描述本发明，但不以任何方式限制本发明的范 围。
[0069] 本发明提供的量子逆向传播神经网络的数字全息压缩传输方法的流程框图如图1 所示，该方法适合处理菲涅尔离轴全息图，包括由计算机仿真产生的或由光学系统记录的， 并不限于数字全息图。在本发明实施例中，本发明提供方法具体包括如下步骤：
[0070] 1)将原始图像通过菲涅尔离轴全息计算方法记录为数字全息图；
[0071] 根据菲涅尔离轴全息算法，将原始图像(物体)记录为数字全息图。菲涅尔离轴全 息图记录由菲涅尔衍射和离轴参考光干涉两部分组成。
[0072]其中，菲涅尔衍射(Fresnel diffraction)指的是光波在近场区域的衍射。菲涅尔 衍射积分式可以用来计算光波在近场区域的传播。计算由物体每一点所出射光波至全息记 录平面的菲涅尔衍射并叠加，即可计算物体在全息记录平面的物光波前。
[0073] 离轴参考光干涉指的是，全息图再现时同时存在有再现像与孪生像及零级衍射亮 斑，相互重叠影响再现图像质量。在全息图记录时引入与物光波前存在一定夹角的参考光， 可以将不同级衍射像分离，得到单独清晰的再现图像。
[0074] 根据离轴全息公式，利用干涉记录物光衍射波前复振幅函数形成全息图IH(x，y):
[0075]
[0076] 式1中，IH(x，y)为全息图记录的光强，R(x，y)为参考光波前函数，U(x，y)为原始物 体波前函数。Rlxd)与υ?χ,γ)分别为R(x，y)与U(x，y)的共辄项。
[0077] 物体光波的衍射过程由菲涅耳衍射公式(式2)计算：
[0078]
[0079] 式2中，d为衍射距离，k为波数，j为虚数单位；（XQ，yQ)为初始物波平面，（x，y)为衍 射平面;〇( XQ，yo)为初始物波平面信息，U(x，y)为经过衍射距离为d的物光波前复振幅函数。 使用计算机对该衍射积分进行计算，可以通过式3进行快速傅里叶变换：
[0080]
[0081]式3中，F表示傅立叶变换，d为衍射距离，k为波数，j为虚数单位;Lo为初始物波平 面采样范围；初始物波平面的采样点数为NXN，即初始图像的像素个数;初始物波平面采样 点的间隔为Δχ〇= Ay〇 = L()/N; Δχ= Ay是初始物波平面离散傅立叶变换后采样间隔，在空 间域内，与之对应的为衍射平面的采样间隔。
[0082] 2)构建QBP神经网络，使用步骤1)得到的数字全息图对其进行训练；
[0083] 此过程主要包括建立量子神经元模型，构建和训练QBP神经网络。
[0084] 2.1建立量子神经元模型
[0085] 建立量子神经元模型是构建QBP神经网络的基础，量子神经元是QBP神经网络的基 本组成单元。基于量子门（具体涉及一位相移门和两位受控非门）建立新的量子神经元模 型，实现量子计算。
[0086] 量子门是物理实现量子计算的基础，可以由量子门构成任意的量子门组网络，本 发明采用的就是由一位相移门和两位受控非门组成基本的计算单元，作为神经网络的激活 函数来构成新的量子神经元模型。为了便于应用，用复数形式表示量子态及通用量子逻辑 门组。
[0087] 量子比特与经典比特的不同在于:一个量子比特态| Φ>是处于状态|〇>和11>的相 干叠加态上，用式4表示：
[0088] Φ> = α|〇>+β|?> (式 4)
[0089]其中，α和β是复数表示概率幅，满足归一化要求I α 12+ I β 12 = 1。用一个复数函数来 描述量子态的状态，复数函数的表不形式为式5:
[0090] f (Θ) =e10 = cos9+isin9 (式 5)
[0091 ]其中，i = 表示虚数单元；Θ表示量子态的相位。| 〇>的概率幅用这个复数函数 的实部表不，11>的概率幅用它的虚部表不，则一个量子态可以描述为式6:
[0092] i])> = cos0 | 0>+sin9 11> (式6)
[0093] 根据式5的定义，基于1位相移门和2位受控非门构成的量子神经元模型如图2所 不，其输入使用多个量子态叠加的形式，通过该三位量子门对输入量子态的幅度和相位分 别进行处理，来得到多个量子态叠加的输出。
[0094] 在图2中，Xi(i = l，2，···，M)表示第i个输入到神经元的量子态;0i(i = l，2，···，M)为 权值的相位转移系数;λ为阈值系数;δ为相位控制因子;〇为输出态;arg( u)是对复数u提取 相位，即arg(u) =actag(Im(u)/Re(u))，其中，Im(u)为求复数u的虚部;Re(u)为求复数u的 实部;函数f( ·)的定义如(式5)所示;g( ·)是sigmoid函数。
[0095] 本文设Ii(i = 1，2，…，M)为第i个输入到神经元的量子态Xi = f (Ii) (i = 1，2，…，M) 的相位，则该量子神经元输出表示为如下量子神经元模型：
[0096]
[0097]
[0098] 〇 = f(y)(式 9)
[0099] 在这个量子神经元模型中，存在着两种类型的参数形式，一种是对应于相移门的 相位的权值参数Θ和阈值参数λ;另一种翻转控制参数δ对应于控制非门。与传统神经元不同 的是，量子神经元中权值fOO与输入 Xl = f(L·)相乘的结果，是通过相移门对神经元态的相 移来实现的。
[0100] 2.2构建〇8?神经网络
[0101]根据上面提出的量子神经元模型，借助传统BP神经网络用于图像压缩传输的原理 和网络结构，建立用于图像压缩传输的QBP神经网络。为了减少系统的运行时间和提高网络 的效率，所建立的用于图像压缩的QBP网络隐含层只有一层，其网络结构如图3所示。可以看 出其网络结构与传统的三层BP网络相同，包含一个输入层、一个隐含层和一个输出层，同层 神经元之间无连接，层与层之间各个神经元则实现全连接。不同的在于QBP神经网络中每一 层的神经元（图3中网络的黑色结点）为量子神经元，并且其输入和输出是多个量子态的叠 加。
[0102]在图3所示的三层QBP网络中，{ΙΝι}(1 = 1，2,…，L)、{HIDEk}(k=l，2，."，KWP {〇υΤη}(η = 1，2，···，Ν)分别表示输入层、隐含层和输出层的量子神经元集合;L，K和N分别对 应各层的神经元个数。
[0103]当归一化的输入值inputi(l = l，2，···，L)被输入到输入层{ΙΝι}，输入层各神经元 将使介于[0，1 ]之间的输入值转换为介于[0，V2 ]之间量子态XI (1 = 1，2，…，L)的相位值11 (1 = 1，2，…，L)，然后将其量子态输出至隐含层。表达式为：
[0104]
)
[0105] xi = f(Ii)(式11)隐含层{HIDEk}和输出层{0UTn}的处理过程依据式7-9进行。
[0106] 在量子神经元里，量子态|1>相当于为神经元的激活态，量子态|0>相当于神经元 的抑制态，那么任意一个神经元的量子态被定义为激活态和抑制态的叠加态，其神经元的 最终输出值〇utput n(n = 1，2，…，N)为激活态11>时的概率。所以，三层QBP神经网络输出层 第η个神经元的最终输出值为：
[0107] outputn= | Im(On) |2 (式12)
[0108] 2.3训练QBP神经网络
[0109]为训练QBP神经网络，借助复数BP算法，定义网络中的逆向传播为量子逆向传播算 法(QBP算法），采用近似均方误差的最速下降算法，来调整网络中相位转动系数Θ、阈值系数 λ和相位控制因子δ，使训练误差的均值小于期望目标为止，具体表达为：
[0110]
[0111]
[0112]
[0113] 式12-14中，η为学习速率。Etotai为均方误差函数，其表达式为：
[0114]
[0115]其中，B代表训练样本的总个数，第b个训练样本由第b个图像块产生。targetn,b和 outputn,b分别代表输入第b个训练样本时，输出层第η个神经元的目标输出和实际输出。
[0116] QBP神经网络的逆向传播训练过程如图4:
[0117] (1)初始化QBP神经网络，将步骤1)得到的大小为X X Υ的数字全息图归一化处理， 然后划分为XP X yP的图像块，每个图像块变为Μ X 1 (Μ = χΡ X yP)维的向量，作为此QBP神经网 络的训练(学习）样本。
[0118] (2)将一个训练样本输入到QBP神经网络中计算网络输出。
[0119] (3)记录网络实际输出与目标输出的误差，根据最速下降法，逆向逐层调整各层参 数。
[0120] (4)重复第(2)、（3)步，直至所有训练样本都已输入。
[0121] (5)计算所有训练样本的实际输出与目标输出的累计总误差，即输入数字全息图 与网络处理后输出的数字全息图之间的误差，训练次数加1，若此总误差小于设定总误差或 训练次数大于设定训练次数，网络训练结束，否则重新进入步骤(2)。网络总误差定义如下：
[0122]
[0123] 其中，g(b，m)为归一化处理后的输入神经网络的数字全息图，沉δ,m)是经过QBP 神经网络压缩解压缩处理后的归一化数字全息图。
[0124] 3)使用步骤2)得到的训练完毕的QBP神经网络对步骤1)得到的数字全息图进行压 缩、传输和解压缩，得到重构后的全息图；
[0125] 将步骤1)得到的数字全息图归一化处理后输入到训练完毕的QBP神经网络的输入 层，通过细腰性的隐含层压缩传输数字全息图。其中，通过输入层到隐含层实现图像的编 码/压缩，通过隐含层到输出层实现图像的解码/解压缩，得到重构后的归一化全息图（对应 式12的最终输出值）；将此全息图再经过反归一化处理后得到重构后的全息图，用于步骤 4) 〇
[0126] 图像数据输入到输入层后，强迫其通过细腰性的隐含层达到压缩传输的目的，然 后由隐含层到输出层实现图像的解码重建。因此，在图3所示的3层QBP网络中，网络的输入 层和输出层通常取相等的神经元数，即L = N，其大小由原图像在压缩的过程中根据具体情 况来确定，每个神经元对应一个像素;隐含层神经元数为K，且K〈N，通过调整K的大小，实现 调节图像压缩率的功能。
[0127] 4)对步骤3)得到的重构后的全息图进行再现，得到再现像。
[0128] 对重构全息图进行再现，具体根据菲涅尔离轴全息算法，将数字全息图再现为原 始图像，需要通过离轴参考光照射与菲涅尔衍射两个过程实现菲涅尔离轴全息图再现。
[0129] 根据离轴全息再现公式，利用与参考光一致的再现光R(x，y)照射全息图IH(x，y) 进行再现：
[0130] IH(x,y) XR(x,y) = ( |R(x,y) |2+|U(x,y) |2)R(x,y)+|R(x,y) |2U(x,y)+R2(x,y)U* (x,y) (式 16)
[0131]
[0132] 式16中，IH(x，y)为全息图记录的光强，R(x，y)为参考光与再现光波前函数，U(x， y)为原始物体再现波前函数，#( X，y)为再现像的共辄项。
[0133] 经过再现光照射的数字全息图光波的衍射过程由菲涅尔衍射公式(式2~式3)计 算，在衍射平面可以获得原始图像的再现像。
[0134] 为了测试本发明的压缩传输速度（即网络训练速度)与再现图像质量，采用与BP神 经网络对比训练次数的方法衡量压缩传输速度，采用压缩率(R)、均方误差(MSE)与峰值信 噪比(PSNR)等评价方法衡量再现图像的质量。
[0135] R，MSE，PSNR 分别定义为：
[0139] 其中，S。代表压缩后的数字全息图数据大小，S。代表原始数字全息图数据大小。f (x，y)代表原始全息图的再现图像函数，代表经过神经网络压缩和解压缩处理后的
[0136]
[0137]
[0138] 全息图的再现函数。Xpeak代表图像数据的峰峰值。
[0140] 图5为原始lena图像的数字全息图及其再现像。图5(a)为原始lena图像（128X128 像素）；图5(b)为lena图像的数字全息图（256X256像素）；图5(c)为lena图像的再现像。在 本实施例中，我们使用图5(b)中的数字全息图经过归一化，然后作为QBP神经网络的输入。 设定QBP和BP神经网络的x P = yP = 8,因此该数字全息图被划分为1024个训练样本。神经网络 的输入输出层神经元个数为L = N = p2 = 64，隐含层神经元个数K取值在{32，16，8，4}变化， 分别对应不同的压缩率{〇 · 5000，0 · 2500，0 · 1250，0 · 0625}。
[0141] 图6为QBP和BP神经网络在上述四个不同压缩率下的训练次数。图6(a) (b)中R = 0.5000 3^狀=0.0005;图6(。）（(1)中1? = 0.2500 3^狀=0.0015;图6(6)(〇中1?=0.1250， 已^狀=0.0023;图6(8)(11)中1? = 0.0625 3^狀=0.0030。其中训练次数上限设为5000。网 络总误差的最大值Errmax根据不同的压缩率设定。学习速率η取值范围为[0.05,0.70]。
[0142] 表1QBP与ΒΡ不同压缩率下的实验结果数据对比
[0143]
[0144] 表1是不同压缩率下，本发明提供的QBP压缩方法与ΒΡ方法的实验结果数据对比。 可以看出，在最佳学习速率下，QBP神经网络的训练次数远小于相同压缩率下的ΒΡ神经网 络，因此本发明的数字全息图压缩传输方法的速度更快。
[0145] 在本实施例中，我们对压缩率为0.5000的情况下QBP和ΒΡ神经网络压缩处理图像 的结果进行了质量评估。图7为在最佳学习速率nQBP = 0.40和ηΒΡ = 0.20时，得到的QBP和ΒΡ的 处理后的数字全息图和再现图像的对比。其中，图7(a) (b)为基于ΒΡ神经网络处理后的数字 全息图和再现图像；图7(c)(d)为基于QBP神经网络处理后的数字全息图和再现图像。表2为 QBP与BP的再现图像质量评估参数。可以看出，在此压缩率下，本发明采用的QBP方法比BP方 法能够得到更好的图像恢复质量。
[0146] 表2QBP与BP再现图像质量评估参数(R = 0 · 5000)
[0147]
[0148] 综上所述，本方案在三维全息信息压缩传输中具有良好的应用前景。本方案思路 对于数字全息图压缩传输技术在压缩传输速度和图像恢复质量上均能达到较优水平。
[0149] 需要注意的是，公布实施例的目的在于帮助进一步理解本发明，但是本领域的技 术人员可以理解:在不脱离本发明及所附权利要求的精神和范围内，各种替换和修改都是 可能的。因此，本发明不应局限于实施例所公开的内容，本发明要求保护的范围以权利要求 书界定的范围为准。
【主权项】
1. 一种量子逆向传播神经网络的数字全息压缩传输方法，包括如下步骤： 1) 将原始图像通过菲涅尔离轴全息计算方法记录为数字全息图； 2) 构建量子逆向传播神经网络，所述量子逆向传播神经网络包含一个输入层、一个隐 含层和一个输出层，所述量子逆向传播神经网络每一层包含多个量子神经元;使用步骤1) 得到的数字全息图对所述量子逆向传播神经网络进行训练，得到训练完毕的量子逆向传播 神经网络. 3) 使用步骤2)得到的训练完毕的量子逆向传播神经网络对步骤1)得到的数字全息图 进行压缩、传输和解压缩，得到重构后的全息图； 4) 对步骤3)得到的重构后的全息图进行再现，得到再现像。2. 如权利要求1所述量子逆向传播神经网络的数字全息压缩传输方法，其特征是，步骤 1) 所述记录包括菲涅尔衍射计算方法和离轴参考光干设计算方法； 所述菲涅尔衍射计算方法由菲涅耳衍射公式计算，所述菲涅耳衍射公式为式2:巧2) 式2中，d为衍射距离;k为波数;j为虚数单位；（xo，yo)为初始物波平面；（x，y)为衍射平 面;O(x〇，y〇)为初始物波平面信息，U(x，y)为经过衍射距离为d的物光波前复振幅函数； 所述离轴参考光干设计算方法利用干设记录物光衍射波前复振幅函数形成全息图Ih (x，y)，如式 1:式1中，lH(x，y)为全息图记录的光强;R(x，y)为参考光波前函数;U(x，y)为原始物体波 前函数;R*(x，y)与U*(x，y)分别为R(x，y)与U(x，y)的共辆项。3. 如权利要求2所述量子逆向传播神经网络的数字全息压缩传输方法，其特征是，通过 式3进行快速傅里叶变换，对菲涅尔衍射进行计算：式3中，F表示傅立叶变换;d为衍射距离;k为波数;j为虚数单位;Lo为初始物波平面采样 范围；初始物波平面的采样点数为NXN，即初始图像的像素个数;初始物波平面采样点的间 隔为A χ〇= Δ y〇 = Lo/N; Δ χ= Δ y是初始物波平面离散傅立叶变换后采样间隔，在空间域 内，与之对应的为衍射平面的采样间隔。 4 .如权利要求1所述量子逆向传播神经网络的数字全息压缩传输方法，其特征是，步骤 2) 具体包括如下步骤： 21)基于量子口建立新的量子神经元模型，用于实现量子计算;所述新的量子神经元模 型中，权值与输入相乘的结果通过相移口对神经元态相移实现： 设Ii(i = l，2,…，M)为第i个输入到神经元的量子态xi = f (Ii)(i = l，2,…，M)的相位， 该量子神经元输出表示为如下量子神经元模型：其中，θια = 1，2，···，Μ)为权值的相位转移系数;λ为阔值系数;δ为相位控制因子;0为输 出态;日巧(11)是对复数11提取相位，即日'旨(11)=日(31:日旨（1111(11)/1^(11))，其中，1111(11)为求复数11 的虚部;Re(u)为求复数U的实部;g(S)是sigmoid函数;函数f( ·)的定义如式5所示； f (白）= e" = cos 目+isin 目（式 5) 其中，? = 表示虚数单元;Θ表示量子态的相位； 22) 构建用于图像压缩传输的量子逆向传播神经网络，所述量子逆向传播神经网络每 一层的量子神经元的输入和输出是多个量子态的叠加； 23) 定义神经网络的逆向传播训练方法，训练步骤22)构建的量子逆向传播神经网络； 所述神经网络的逆向传播训练方法采用近似均方误差的最速下降算法，调整所述量子逆向 传播神经网络中相位转动系数Θ、阔值系数λ和相位控制因子δ，使训练误差的均值小于期望 目标时结束训练，得到训练完毕的量子逆向传播神经网络。5. 如权利要求4所述量子逆向传播神经网络的数字全息压缩传输方法，其特征是，步骤 21) 所述量子口由一位相移口和两位受控非口组成，作为神经网络的激活函数，用复数形式 表示量子态和通用量子逻辑口组。6. 如权利要求4所述量子逆向传播神经网络的数字全息压缩传输方法，其特征是，步骤 22) 中，所述量子神经元的输入和输出具体是： 设定{I化}α = l，2，...，L)、化IDEkKk=l，2，...，K)和{0UTn}(n=l，2，...，N)分别表示所 述量子逆向传播神经网络的输入层、隐含层和输出层的量子神经元集合;L，K和Ν分别对应 各层的神经元个数;当归一化的输入值?ηριι?ι(1 = 1，2，···，υ被输入到输入层{I化}，输入层 各神经元将使介于[〇，1]之间的输入值转换为介于[〇，V2]之间量子态χι(1 = 1，2，···，υ的 相位值11(1 = 1，2,…，L)，然后将其量子态输出至隐含层;表达式为式10和式11:xi = f (Ii)试 11) 其中，函数f( ·)的定义如式5所示； 所述隐含层化IDEk}和所述输出层{ΟυΤη}的处理过程按照式7~式9; 定义任意一个量子神经元的量子态为激活态和抑制态的叠加态，所述神经元的最终输 出值OU化ιι?η(η=1，2，···，Ν)为激活态|1〉时的概率;所述量子逆向传播神经网络输出层第η 个神经元的最终输出值表示为式12: OUtpUtn= I Im(On) I 2 (式 12) 式12中，0为输出态。7. 如权利要求4所述量子逆向传播神经网络的数字全息压缩传输方法，其特征是，步骤 23)中，所述神经网络的逆向传播训练方法使用步骤1)得到的数字全息图对所述量子逆向 传播神经网络进行训练，具体包括如下步骤： 231) 初始化所述量子逆向传播神经网络，将步骤1)得到的大小为XXY的数字全息图归 一化处理，然后划分为XpXyp的图像块，每个图像块变为MXl(M=XpXyp)维的向量，作为所 述量子逆向传播神经网络的训练样本； 232) 将一个训练样本输入到所述量子逆向传播神经网络中计算网络输出； 233) 记录网络实际输出与目标输出的误差，根据最速下降法，逆向逐层调整各层参数； 234) 重复步骤232)~233 )，直至所有训练样本都已输入； 235) 设定训练次数，设定累计总误差为输入数字全息图与网络处理后输出的数字全息 图之间的误差;计算所有训练样本的实际输出与目标输出的累计总误差，训练次数加1;当 所述累计总误差小于设定的累计总误差或所述训练次数大于设定的训练次数时，训练结 束;否则重新进入步骤232)。8. 如权利要求1所述量子逆向传播神经网络的数字全息压缩传输方法，其特征是，步骤 3)具体包括:将步骤1)得到的数字全息图归一化处理后输入到所述训练完毕的量子逆向传 播神经网络的输入层，通过输入层到隐含层实现图像的编码/压缩，通过隐含层到输出层实 现图像的解码/解压缩，得到重构后的归一化全息图；将所述重构后的归一化全息图再经过 反归一化处理，得到重构后的全息图。9. 如权利要求8所述量子逆向传播神经网络的数字全息压缩传输方法，其特征是，所述 量子逆向传播神经网络的输入层和输出层的量子神经元数量相等，记作L = N;每个神经元 对应一个像素;所述隐含层的量子神经元数量记为Κ，Κ小于N;通过调整所述隐含层的神经 元数量Κ的大小，达到调节图像压缩率的目的。10. 如权利要求1所述量子逆向传播神经网络的数字全息压缩传输方法，其特征是，步 骤4)所述再现具体通过离轴参考光照射过程和菲涅尔衍射过程实现菲涅尔离轴全息图再 现。
【文档编号】G06T9/00GK105976408SQ201610273648
【公开日】2016年9月28日
【申请日】2016年4月28日
【发明人】杨光临, 刘梦佳
【申请人】北京大学