技术领域
[0001] 本
发明属于航天遥感技术领域中的光学卫星
图像处理的一种量子信息处理方法，特别是一种航天光学遥感多路图像量子滤波与合成方法。
背景技术
[0002] 目前光学遥感器在轨成像阶段，因空间环境变化等多种复杂因素造成探测元器件性能退化，再加上探测器件固有特性，在焦面出现大量噪声，造成成像
质量退化。目前，卫星在轨阶段，难以有效提升成像质量，需要通过数字图像处理的手段进一步提升像质。
[0003] 由于探测器像元规模受限，造成光学成像卫星在轨成像
分辨率难以提高，
采样到信息量较低；另外因空间环境变化，很多复杂因素造成探测元器件性能退化，或因探测器件固有特性，退化图像质量，影响图像判读，采用传统方法往往难以得到理想效果。
发明内容
[0004] 本发明技术解决问题：克服
现有技术的不足，提供一种航天光学遥感高密度量子化信息处理方法，在量子空间对
信号进行处理，同时采用多路量子信息合成的方法进一步提升数据的信息量，能够更有效的去除噪声，大大提升了图像的有效信息量。
[0005] 本发明技术解决方案：一种航天光学遥感图像量子滤波与合成方 法，包括如下步骤：
[0006] 步骤一、获取光学遥感卫星对目标成像的多路观测信息，并结合光学遥感成像系统信息，对成像系统建模，得到多路光学遥感图像；
[0007] 步骤二、将多路光学遥感图像转换到量子空间，分解到个
量子比特面上，得到多路光学遥感图像对应的量子比特面中量子态向量；
[0008] 步骤三、将多路光学遥感图像对应的量子比特面中量子态向量进行滤波，通过对多路成像信息的量子比特面信息进行范数优化，实现对量子信息的滤波，得到滤波后的多路量子信息；
[0009] 步骤四、将滤波后的多路量子信息合成，即在量子比特面内对多路分解后的量子信息进行重采样，进行量子信息合成，得到合成后的量子信息；
[0010] 步骤五、信息重构，即将合成后量子信息重构转换回图像信息，得到处理后的高清晰图像。
[0011] 所述步骤一中，成像系统建模公式如下
[0012]
[0013] 其中:(m,n)为获取到光学遥感图像数据的尺寸，其中m为图像在横轴方向的尺寸，n为图像在纵轴方向的尺寸；H为PSF的二维高斯状的s×s矩阵表示，由测量得到，s为点扩展函数PSF在一个方向的尺寸；I(m,n)为理想场景信息，N为噪声项，f(m,n)为光学遥感图像。
[0014] 所述s取值为：
[0015] 所述步骤二中，将多路光学遥感图像转换到量子空间表示如下：
[0016]
[0017] 上式是光学遥感图像转换到量子空间转化表示，其中t代表多路成像信息，b代表成像信息的量化位数，p(t,b)代表第t路信息内第b个比特面f的概率的均方，|ψt,b(μ,ν)>为图像f对应第t路信息内第b个量子比特面中量子态向量，＜ψt,b(μ,ν)|代表图像对应量子态向量的左矢，为|ψt,b(μ,ν)>的共轭转置，(μ,ν)代表量子态向量在量子比特面上的
位置。
[0018] 得到多路光学遥感图像对应的量子比特面中量子态向量为：
[0019]
[0020] 上式中θ和φ分别表示量子比特空间上对应向量与z轴夹
角和投影在xy面与x轴之间的夹角。
[0021] 所述步骤三中，多路光学遥感图像对应的量子比特面中量子态向量进行滤波，表示如下
[0022]
[0023] 其中：λt,b为每个比特面的调整系数，取0～1之间； 表示求2范数最小化，||·||1表示1范数；量子范畴的范数空间近似于图像域的范数空间，即在不同比特面内求得的范数表达，其展开运算遵守量子运算规则，通过优化使得 趋近于理想值，最终得到优化后的光学遥感图像对应的量子态向量；
[0024] 为光学遥感成像系统PSF对应的矩阵H在第t路第b个量子比特面量子态向量；
[0025] 为理想光学遥感图像I在第t路第b个量子比特面量子态向量；
[0026] 为光学遥感图像第t路第b个量子比特面上的量子态向量；
[0027] 表示光学遥感成像系统噪声N对应第t路第b个量子比特面上的量子态向量。
[0028] 所述步骤四中，在各个量子比特面内对多路分解后的量子信息进行重采样，信息合成得到
[0029]
[0030] |ψ't,b(μ,ν)>为得到的合成后的量子信息，上式重采样采用插值的方法；t为多路成像信息中的第t路信息(本发明取t＝4)，b为量子比特面数(本发
明取b＝8)。
[0031] 所述步骤五中，信息从量子空间转换回
图像空间，表示为
[0032]
[0033] 其中，|ψ't,b(μ,ν)>为滤波后的量子态向量；＜ψ't,b(μ,ν)|为滤波后量子态向量的左矢，为|ψ't,b(μ,ν)>的共轭转置；f'(m,n)为转换后的图像。
[0034] 所述H由两种测量方法得到，第一种是实验室通过点源靶标，实现对系统点扩展函数PSF的精确测量；第二种是在卫星过境星下点铺设地面靶标，实现对卫星动态点扩展函数PSF的测量。
[0035] 所述步骤(三)中对多路光学遥感图像对应的量子比特面中量子态向量进行滤波时，量子向量之间的运算满足如下法则：
[0036] 且 则有 和 且有 和
[0037] 上式中代表任意两个量子态向量之间的运算关系，其中量子态向量|ψj>的矩阵表示为 而|ψk>的矩阵表示为 因此|ψj>与|ψk>之间的运算关系均满足展开式 与之间运算。这样量子态向量之间的运算关系就可以由矩阵之间运算来实现。
[0038] 所述插值方法采用双线性插值、三次插值等插值方法。
[0039] 本发明与现有技术相比的优点在于：
[0040] (1)本发明采用量子比特的概念，将
空域信息转换到量子空间分析。在量子比特面对图像进行滤波、信息合成，实现良好的滤波效果。由于信息是通过量子比特的形式表达的，能够表示出基态之间的
叠加态，这种存在于叠加态的噪声信息是传统图像处理方法无法处理的，同时采用多路量子信息合成一路量子化信息的方法。由于合成后的信息结合了多路信息，可提高图像的有效信息量，如峰值
信噪比、模糊度等指标。
[0041] (2)本发明对多路成像信息采用量子比特转换的方法将图像转换到量子空间，考虑到各种叠加态信息，因而滤波更明显。
[0042] (3)本发明适于在星上分辨率受限，存在大量噪声难以处理的图像，适于在星上无法处理或地面难以滤除的图像，或分辨率难以提 升的图像，在有效降低计算复杂度后，能够实现星上实时处理。
附图说明
[0043] 图1为本发明的量子图像处理的主要过程示意图；
[0044] 图2为原图中量子比特面下的分解示意图，其中a.第1层分解图，b.第2层分解图，c.第3层分解图，d.第4层分解图，e.第5层分解图，f.第6层分解图，g.第7层分解图，h.第8层分解图；
[0045] 图3为原始图及本发明处理结果比较，其中a.遥感器获取图像，b.传统处理方法(滤波+插值)，c.为本发明量子空间处理方法。
具体实施方式
[0046] 如图1所示，本发明方法具体实现如下：
[0047] 步骤(一)，获取光学遥感卫星对目标成像的多路观测信息，并结合光学遥感成像系统信息，对成像系统建模，得到多路光学遥感图像。数值经过归一化后，在[0,1]区间。光学遥感成像系统可视为线性系统，表示如下：
[0048]
[0049] 其中(m,n)为获取到图像数据的尺寸，s为PSF在一个方向的尺寸，I(m,n)为理想场景信息，N为噪声项；
[0050] 上式中H代表光学遥感成像系统中的点扩展函数PSF，可由测量得到。测量方法有两种，一种是实验室通过点源靶标，实现对系统点扩展函数PSF的精确测量。第二种是在卫星过境星下点铺设地面靶标，实现对卫星动态点扩展函数PSF的测量。总之，测得的PSF以二维矩阵H形式表达，表示为高斯状的s×s矩阵，为运算方便取
[0051] 步骤(二)，将多路光学遥感图像转换到量子空间，分解到个量子比特面上，得到多路光学遥感图像对应的量子比特面中量子态向量。
[0052]
[0053] 上式子(2)是多路信息从图像空间到量子空间的转化表达，其中t代表多路成像信息，b代表成像信息的量化位数，p(t,b)代表第t路信息内第b个比特面f的概率的均方，|ψt,b(μ,ν)>为图像f对应第t路信息内第b个量子比特面中量子态向量；＜ψt,b(μ,ν)|代表图像对应量子态向量的左矢，为|ψt,b(μ,ν)>的共轭转置。
[0054] 信息量子化建模，将图像f归一化，即f(m,n)满足f(m,n)∈[0,1],此时图像f(m,n)的量化比特可表达为
[0055]
[0056] 这里，|0>和|1>分别表示为图像量子比特中的两个基态。如果α>和β>可理解为光2 2
入射后形成不同的偏振态，且满足α+β＝1。
[0057] 进一步推导为
[0058]
[0059] 式子(4)中因子eiξ对信息没有明显的效果，因而可忽略掉。θ和φ分别表示量子比特空间上对应向量与z轴夹角和投影在xy面与x轴之间的夹角；
[0060] 本发明
实施例中t＝4，即有4路信息，每个通道内有8个比特面，得到32个量子态向量。
[0061] 步骤(三)，将多路光学遥感图像对应的量子比特面中量子态向量进行滤波，通过对多路成像信息的量子比特面信息进行范数优化，实现对量子信息的滤波，得到滤波后的多路量子信息；量子信息滤波过程：对每一路成像信息进行滤波。通过对每一路成像信息量子态向量|ψt,b(μ,ν)>进行范数最小凸优化，实现对量子信息的滤波，即
[0062]
[0063] 这里λt,b为每个比特面的调整系数，这里取0～1之间； 表示求2范数最小化，||·||1表示1范数。通过优化使得 趋近于理想值，最终得到优化后的光学遥感图像对应的量子态向量。
[0064] 为光学遥感成像系统PSF对应的矩阵H在第t路、第b个量子比特面量子态向量；
[0065] 为理想光学遥感图像在第t路、第b个量子比特面量子态向量；
[0066] 为光学遥感图像在第t路、第b个量子比特面上的量子态向量；
[0067] 表示光学遥感成像系统噪声N对应第t路、第b个量子比特面上的量子态向量；
[0068] 量子范畴的范数空间近似于图像域的范数空间，即在不同比特面内求得的范数表达，其量子向量运算遵守量子运算规则：
[0069] 定理 ：若 且 则有 和且有 和
[0070] 上式中代表任意两个量子态向量之间的运算关系，其中量子态向量|ψj>的矩阵表示为 而|ψk>的矩阵表示为 因此|ψj>与|ψk>之间的运算关系均满足展开式与 之间运算。这样量子态向量之间的运算关系就可以由矩阵之间运算来实现。
[0071] 步骤(四)，将滤波后的多路量子信息合成，即在量子比特面内对多路分解后的量子信息进行重采样，进行量子信息合成，得到合成后的量子信息；量子信息合成过程，即在第b个比特面内多路分解后的量子信息进行重采样，合成高密度的量子信息，即[0072]
[0073] 重采样采用插值的方法，具体各类插值方法为业内公知(如双线性插值、三次插值等插值方法)，这里取 t为多路成像信息中的第t路信息，b为量子比特面数；
[0074] 步骤(五)，信息重构，即将合成后量子信息重构转换回图像信息，得到处理后的高清晰图像，表示如下
[0075]
[0076] 其中，|ψ't,b(μ,ν)>为滤波后的量子态向量；＜ψ't,b(μ,ν)|为滤波后量子态向量的左矢，为|ψ't,b(μ,ν)>的共轭转置；f'(m,n)为转换后的图像。最终，得到处理后的高清晰图像。
[0077] 实施例1
[0078] 本发明实施例中选用全色谱段成像，图像选取航母作为仿真图，如图3中的(a)，同时包含稀疏的海洋背景信息。实施例中t＝4，每路信息分解8个比特面，共32个图像的量子态向量。如图1中所示，通过获取上述4路成像信息进行结合系统先验信息点扩展函数PSF的量子滤波方法，然后合成高密度的量子信息，最后重建滤波后的高信息量图像。图2为不同量子面内的分解图，层数越大信息分解的越精细。通过逐层分解可以将信号和噪声在量子空间精细分解。
[0079] 图3中的(b)可见，在图像空间采用传统图像滤波方法(均值滤波+线性插值)合成一幅图像，图像中的去除了部分噪声，同时丢失很多细节信息。而图3中的(c)量子滤波方法考虑到更多叠加态特性并且在量子空间合成多路量子信息，能够更有效去除噪声，增强图像的细节信息，得到的图像像质更清晰。
[0080] 如表1为采用传统图像滤波方法与本发明量子方法处理后图像质量的评价结果表。由表可见，本发明处理后图像的细节信息评价指标与传统图像滤波方法相比更高，说明处理后的图像质量优于在传统图像滤波处理结果。
[0081] 表1传统处理方法与本发明量子化信息处理方法后对比(海面舰船图)
[0082]
[0083] 以上所述，仅为本发明最佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉
本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。
[0084] 本发明
说明书中未作详细描述的内容属于本领域专业技术人员的公知技术。
