基于压缩感知的高动态范围图像去伪影融合方法
专利汇可以提供基于压缩感知的高动态范围图像去伪影融合方法专利检索，专利查询，专利分析的服务。并且本 发明 公开了一种基于 压缩 感知 的高动态范围图像去伪影融合方法。首先，对输入的多曝光 图像序列 进行压缩 采样 ；接着，使用重构方法进行重构得到压缩感知后的多曝光图像序列；然后归一化经过压缩感知的图像集，对图像集使用基于PatchMatch和秩最小化 算法 进行多曝光图像去伪影融合得到目标的高动态（high dynamic range,HDR）图像。本发明利用K-SVD字典学习、压缩感知和去伪影融合的最新研究成果，能够有效降低采样率、存储空间和计算复杂度，得到去除伪影和模糊的HDR图像。，下面是基于压缩感知的高动态范围图像去伪影融合方法专利的具体信息内容。
1.一种基于压缩感知的高动态范围图像去伪影融合方法，其特征在于包括有如下步骤：
1)对输入多曝光图像序列进行分块压缩采样；
2)对压缩采样后的图像块进行LDR图像序列的重构；
3)对经过压缩感知后的多曝光图像序列进行高动态范围图像的去伪影融合。
2.根据权利要求1所述的基于压缩感知的高动态范围图像去伪影融合方法，其特征在于，所述步骤1)的具体实现过程为：
11)以当前输入的LDR图像为参考，进行K-SVD字典学习，得到用于稀疏表示的超完备字典ψ；
12)对分块后的图像信号矩阵F进行稀疏表示，得到F＝ψθ，设计观测矩阵Φ，则压缩采样的过程重写为：
y＝ΦF＝Φψθ＝Θθ\*MERGEFORMAT(1)。
3.根据权利要求1所述的基于压缩感知的高动态范围图像去伪影融合方法，其特征在于，所述步骤11)还包括如下内容：
111)判断当前的LDR图像是否为灰度图像，若不是，则转换为灰度图像，并进一步转换为double类型的图像；
112)选择所需的图像分块大小，将111)中获取的灰度图像进行分块，并重排图像块为矩阵列，形成字典学习的数据集矩阵，也即待稀疏表示的信号矩阵Y，进一步将矩阵Y表示为其中N为矩阵Y的列数；
113)K-SVD算法参数初始化，包括：
numAtom：需要训练的字典元素个数；
numIteration：迭代次数
errorFlag：等于0则表示每个信号的稀疏系数个数固定，需要配置参数L；否则表示信号的稀疏系数不固定，需要配置参数errorGoal；
errorGoal：最大允许表示误差
preserveDCAtom：等于1则表示字典的第一个原子为常量；
InitializationMethod：若为“DataElements”则用信号本身作为初始字典，若为“GivenMatrix”则用给定的矩阵作为初始字典，并对初始字典进行归一化处理，初始字典表示为ψ(0)∈RN×K；
初始化算法迭代次数J＝1；
114)稀疏编码阶段：使用正交匹配追踪算法OMP来计算Y中每一个列向量在当前字典上的稀疏表示所对应的系数向量xi，通过求解下式中的目标函数：
其中 为由系数向量所构成的矩阵，i表示列向量的索引，|| ||F为Frobenius范数，|| ||0为l0范数，T0为稀疏表示稀疏中非零分量的数目的上限，即系数向量中的最大差异度；
115)字典更新阶段：初始字典往往不是最优的，满足稀疏性的稀疏矩阵表示的数据和原数据可能会有比较大的误差，因此需在满足稀疏的条件下逐行逐列进行字典的更新优化，减小整体误差，逼近最优字典，对于ψ(J-1)中的每一列k＝1,2,...K，定义使用到字典中第k个原子ψk的所有信号集合{yi}的索引所构成的集合为：
其中 为系数矩阵X中原子ψk所对应的第k行，不同于X的第k列xk的转置，则目标函数式(2)重写为：
上式中，字典与系数矩阵的乘积ψX被分解为K个秩为1的矩阵的和的形式，当前更新的矩阵列为第k列，则其余K-1个项固定，剥离第k个原子的贡献，得到的矩阵Ek为去除原子ψk的成分其余原子所造成的误差；
116)由于 中0的存在，用SVD得到的 更新向量中的非零元素的数量和位置和原 中的都不相同，即产生了“发散”，因此需要去除 中所有0元素，定义矩阵Ωk为N×ωk的矩阵，它在(ωk(i),i)处的值为1，其他位置上的值都为0，则我们定义 Y、Ek去除零项后的收缩结果为 的长度为|ωk|， 是用到原子ψk的信号的集
合， 为剥离原子ψk的误差，两者都为n×|ωk|大小的矩阵，则式(4)进一步转换为：
117)对 进行SVD分解，即 矩阵U的第一列表示为 即为列ψk的更新项，
更新的 为矩阵V的第一列乘上Δ(1,1)，J＝J+1进入下一次迭代直至满足迭代条件为止，迭代条件根据设定，选择为选定的迭代次数或者是误差值，迭代完成后，即得到最优的超完备字典ψ。
4.根据权利要求3所述的基于压缩感知的高动态范围图像去伪影融合方法，其特征在于，所述112)图像块大小选为8×8，16×16，32×32。
5.根据权利要求1所述的基于压缩感知的高动态范围图像去伪影融合方法，其特征在于，所述步骤12)还包括：
121)观测矩阵Φ的行数为选择的图像块大小的平方乘上采样率，用M表示，列数定义为图像块的大小的平方，用N表示，此外传感矩阵Θ＝Φψ应满足RIP性质(有限等距性质)：保证观测矩阵不会把两个不同的稀疏信号映射到同一个集合中，保证原空间到稀疏空间的一一映射关系，定义观测矩阵Φ的RIP常数δf为满足下式的最小值：
其中F即为分块后的图像信号矩阵，也即稀疏表示后的信号，假设为f稀疏的，即含有f个非零元素，若δf≤1，测成矩阵满则f阶RIP条件。
6.根据权利要求5所述的基于压缩感知的高动态范围图像去伪影融合方法，其特征在于，所述步骤2)的具体实现过程为：
21)初始化：令标签索引集 残差向量r0＝y，y即为步骤12)中所获得的观测向量，迭代次数t＝1；
22)辨识：在传感矩阵Θ中求出与残差向量rt-1最相关的所对应的列向量λt，即内积值最大时所对应的索引(列序号)：
其中N为传感矩阵的列数，Θj表示传感矩阵的第j列；
23)更新：更新标签索引集Λt＝Λt-1∪{λt}，将找到的原子所对应的列向量 添加到集合 其中，t＞0，Λt表示t次迭代的索引，Θt表示按索引Λt选出的矩阵Θ的列集合(大小为M×t的矩阵)；
24)估计：求y＝Θθ的最小二乘解：
25)更新残差向量
26)令t＝t+1；并不断重复步骤22)至步骤26)，若满足某个迭代停止条件，则停止迭代并进入第b7)步；
27)求得输出系数向量
7.根据权利要求6所述的基于压缩感知的高动态范围图像去伪影融合方法，其特征在于，所述步骤26)还包括：
在上述算法的步骤26)中，其迭代停止条件有以下常见的3种情况：
261)当运行到t＞s时，迭代停止，其中，s表示固定的迭代步数，如M/4；
262)残差向量的能量值小于某个预先给定的常数ε，例如取1e-6；
||rt||2≤ε\*MERGEFORMAT(11)
263)当传感矩阵Θ的任何一列都没有残差向量rt的明显能量时：
||ΘTrt||∞≤ε\*MERGEFORMAT(12)。
8.根据权利要求1至7任一项所述的基于压缩感知的高动态范围图像去伪影融合方法，其特征在于，所述步骤3)的具体实现过程为：
31)定义输入的经过压缩感知后的图像中任意一点处的像素值为B，对灰度值B归一化处理之后该点的灰度值为I：
I＝B/255\*MERGEFORMAT(13)
32)现假设源图像为S，参考图像R，源图像和参考图像合成的图像为L，PatchMatch算法就是一个以参考图像为模板，配准源图像生成图像L的过程；由于PatchMatch算法是处理一对图像，现假设输入图像为I1...IN，以N＝5为例，首先令I3为参考图像R，则I3和I4作为其源图像S，然后令I2和I4作为参考图像R，I1和I5分别作为I2和I4对应的源图像S；
33)现定义PatchMatch算法合成图像L的二次函数：
其中τ为灰度映射函数，Ω为图像R和图像S的图像域，i为图像域上的任意一个像素点，n(i)为以i中心p×p的邻域，其中p为邻域的大小，故j是邻域n(i)上的像素点，R(i)是图像R上的第i个像素点，S(i+u(j))是图像S上的第(i+u(j))个像素点，其中u(j)表示从图像L上的像素点j映射到图像S的偏移量，α为一个归一化的因子 其
中wτ和wu为一对权重函数，wτ(i)表示图像R中像素点i映射到图像L的比重，wu(j)表示图像S中像素点j加上偏移量u(j)映射到图像L的比重；
34)灰度映射函数τ定义如下：
其中灰度映射函数的导数τ′≥0，τ()∈[0,1]，i为图像域Ω上的像素点，故L(i)表示为图像L上的第i个像素点，使用迭代重加权最小二乘法算法求解灰度映射函数，则将灰度映射函数τ的目标函数重写为：
其中求解目标函数τ过程中，τ和权重因子ω更新为：
其中n表示迭代次数，δ为一个很小的正常数，令δ＝10-10；
35)权重函数wτ定义如下：
当图像R中图像域上亮度太暗或太亮(即图像域的像素点的灰度值小于3/255，或者大于252/255)时，像素点将会被clipped，否则，就不clipped；
36)权重函数wμ定义如下：
其中d(·,·)表示输入参数之间的空间距离；υ1,υ2为两个归一化参数，分别取对应空间距离的75百分位数；
37)对于参数x和y，d(x,y)＝||x-y||2，而对于 表示为：在图像R和图像S
的图像域上的任意一个像素点i，取以i为中心，大小为p×p的邻域，得到图像块 和 然后图像块 经过灰度映射函数得到 而图像块 相对于i平移了u(i)得到 最
后求两者的空间距离；同理得 其中τ-1()灰度映射函数的逆函数；
38)通过上面定义的函数知，PatchMatch算法实际上就是求解二次函数的过程，输入图像R和图像S，并分别对两幅图像进行向下采样，分别得到图像R和图像S的金字塔图像集，从金字塔顶端向下迭代，求得在对应每层金字塔图像下合成的图像L和灰度映射函数τ，将此结果作为下一次的迭代的初始值，当迭代完成后即得到最终的配准图像L，依照此方法，即得到输入图像I1 IN配准后的图像L1 LN；
39)利用伽马曲线对已经配准好的图像集进行辐射校准，消除潜在的相机的移动而产生的噪音.定义伽马(gamma)函数为：
gamma＝crγ\*MERGEFORMAT(21)
其中c和γ为常数，c取为1，γ取为2.2；这里将经PatchMatch配准过的图像集用gamma曲线进行辐射校准；
310)对配准后的图像集使用秩最小化算法得到批量的对齐图像，首先列向量化所有输入的矩阵得到矩阵D，并初始化低秩矩阵A和噪声矩阵E，根据增广拉格朗日乘子法经过内迭代和外迭代得到最佳低秩矩阵A，则噪声矩阵E＝D-A，最后对得到的低秩矩阵和噪声矩阵进行调整m×n大小的图像，即得到输入图像L对应的低秩图像和噪声图像；
311)输入对齐后的图像集A，将图像A合成目标HDR图像：
其中nImg表示为输入图像的数量，x∈{r,g,b}，r，g，b为彩色图像的三个通道；A(x)和H(x)分别输入图像和HDR图像的x通道图像，最后融合H(x)即得到HDR图像H。
9.根据权利要求8所述的基于压缩感知的高动态范围图像去伪影融合方法，其特征在于p＝7。
基于压缩感知的高动态范围图像去伪影融合方法
技术领域
背景技术
发明内容
数；
具体实施方式
到图像L的比重，wu(j)表示图像S中像素点j加上偏移量u(j)映射到图像L的比重；
最后求两者的空间距离；同理可得 其中τ ()灰度映射函数的逆
函数；
| 标题 | 发布/更新时间 | 阅读量 |
|---|---|---|
| 一种通过照片直方图调节相机参数的方法 | 2020-05-25 | 512 |
| 使用长寿命放射性同位素的实时伽马相机死区时间确定 | 2020-05-17 | 971 |
| 伽马相机和利用伽马相机检测辐射射线的方法 | 2020-05-12 | 184 |
| 一种多伽马光子符合成像系统及方法 | 2020-05-19 | 1013 |
| 具有轨迹优化的心脏SPECT系统 | 2020-05-26 | 123 |
| 适用于在放射性环境中工作的相机 | 2020-05-24 | 511 |
| 人员携带放射性物质的实时检测系统及方法 | 2020-05-24 | 811 |
| 基于活体小动物伽马相机成像的图像显示及优化方法 | 2020-05-17 | 625 |
| 用于在使用时间编码的图案目标的情况下校准由至少一个传感器相机构成的图像采集传感器的方法 | 2020-05-23 | 600 |
| 相机伽马校正曲线提取方法 | 2020-05-13 | 228 |
高效检索全球专利专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台，是提供专利分析，专利查询，专利检索等数据服务功能的知识产权数据服务商。
我们的产品包含105个国家的1.26亿组数据，免费查、免费专利分析。
专利汇分析报告产品可以对行业情报数据进行梳理分析，涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度，系统通过AI智能系统对图表进行解读，只需1分钟，一键生成行业专利分析报告。
