技术领域本发明涉及计算机图像修复领域，具体涉及一种老电影褪色闪烁的修复方法。
背景技术：
老电影是历史文化的缩影，也是一份极为珍贵的可视化资料。然而由于我国早期电影保存条件有限，再加上转印、播放次数过多等人为因素，现存的老电影基本受到了一定程度的损伤，直接影响了观看效果。随着数字图像处理技术的发展，修复老电影的损伤成为了可能。闪烁是老电影损伤中一种最为常见的问题，现有的研究主要集中在亮度闪烁上，主要的修复方法都需要对图像序列中的闪烁进行建模，然后估计该模型的闪烁参数。周磊等设计了一种基于块的局部区域参数估计和补偿方法，通过对图像序列进行基于块的闪烁模型参数估计和当前帧补偿模型的设计，实现了参照模型下的受损图像亮度修复。谌安军等通过对图像序列中的局部闪烁进行基于余弦基的建模和逐次超松弛迭代的方法，实现了对严重局部亮度闪烁的电影的修复。JulieDelon基于尺度空间理论使用尺度时间均衡的方法实现了全局亮度闪烁的修复。然而，经典老电影不仅有黑白老电影，还有诸多彩色老电影，此类电影同样存在闪烁现象，统称为褪色闪烁。在老电影修复中，相对于亮度闪烁修复，褪色闪烁的研究较少。
技术实现要素：
针对现有技术中存在的问题，本发明提出了一种老电影的褪色闪烁修复方法。本方法基于颜色传递算法和色彩校正算法，首先对序列采用均值比较的方法进行闪烁判断，对判定闪烁的序列进行参考帧选取，若闪烁序列无参考帧，则采用色彩校正来产生参考帧，最后通过颜色传递进行全局褪色闪烁修复。为达到上述目的，本发明才有如下技术方案：一种老电影褪色闪烁的修复方法，包括如下步骤：101：记整个待处理的序列为P，对序列P进行分镜处理，记分镜个数为n，分镜序列为Pi；由于后期处理闪烁使用的颜色传递算法对图像的结构有一定要求，只有在参考图像帧与待处理帧之间的布局结构相似时，处理结果最为理想，故分镜是必须的步骤；102：对分镜序列Pi进行RGB均值统计，判断是否存在闪烁；闪烁具有单帧不易观测的特点，只有连续放映时，才能明显发现，利用RGB均值统计，能客观判断Pi是否存在闪烁；103：若判定Pi序列存在闪烁，则在Pi序列中筛选出一帧适合的参考帧；若无参考帧，则使用色彩校正算法产生参考帧；对于参考帧的选取具有主观因素影响，主要是通过修复人员判定是否在Pi序列存在褪色不明显的图像帧，在无法确定合适参考帧的情况下，在Pi序列中选择褪色现象较轻的一帧进行色彩校正；104：使用颜色传递算法，对Pi序列的每一帧待修复图像进行处理；105：在Pi序列处理完成后，对下一组分镜序列P(i+1)重复步骤102-104处理过程，直至Pi＝Pn。所述步骤103中的色彩校正算法的主要过程为：1)将待处理帧的原始像素值范围映射到[0,1]，并使用公式(1)将其从RGB空间转换到lαβ空间；1αβ=0.57740.57740.57740.40820.4082-0.81650.7071-0.707100.38110.57830.04020.19670.72440.07820.02410.12880.8444RGB---(1)]]>2)使用公式(2)在lαβ空间进行色彩校正处理：{α′=α-mα-CαNβ′=β-mβ-CβN---(2)]]>其中，α，β表示lαβ颜色空间通道内所有像素的原始值；α′，β′表示处理好之后的所有像素值；mα，mβ表示α与β通道均值；Cα，Cβ表示实际需求的通道均值，若图像仅需要色彩校正，则Cα，Cβ均为0；3)使用公式(3)将处理好的参考帧转回RGB空间，并将其像素值映射到原始值范围：RGB=4.4679-3.58730.1193-1.21862.3809-0.16240.0497-0.24391.20450.57740.40820.70710.57740.4082-0.70710.5774-0.816511αβ---(3)]]>在进行映射前，利用公式(4)修正参考帧像素值的溢出，I(s){s=1,s>s=0,s<0---(4)]]>其中，I表示参考帧图像，s表示其各通道像素值。所述步骤104中的颜色传递算法，具体如下：1)将参考帧与待修复帧转到lαβ颜色空间，计算参考帧的平均值参考帧的方差待修复帧的平均值待修复帧的方差2)使用公式(5)弱化待修复帧的整体色彩信息；lS′=lS-mSlαS′=αS-mSαβS′=βS-mSβ---(5)]]>其中，lS，αS，βS是待修复帧三通道的所有像素值，lS′，α′S，βS′是待修复帧三通道弱化后的像素值；3)使用公式(6)将待修复帧与参考帧标准差之比作为通道值偏移量的系数，将参考帧的细节信息映射到待修复帧中；l′=σRlσSl×lS′α′=σRασSα×αS′β′=σRβσSβ×βS′---(6)]]>4)使用公式(7)将参考帧的整体信息添加到处理后的待修复帧中；l=l′+mRlα=α′+mRαβ=β′+mRβ---(7)]]>5)将处理好的图像帧转回RGB颜色空间。本发明与现有技术相比的优点在于：(1)、本发明提供了一种彩色老电影闪烁的修复方案，其交互简单，公式简洁，效率较高，为修复褪色闪烁提供了可能；(2)、采用了色彩校正的方案，在没有参考帧的情况下仍然可以进行修复，扩宽了可修复范围；(3)、在有合适参考帧的前提下，忽略图像颜色空间的变换，本发明同样适用于全局亮度闪烁的修复。附图说明图1为一种彩色老电影闪烁修复方案的流程图。图2为色彩校正对比图，其中(a)为色彩校正处理前图像，(b)为处理后图像。图3为动画老电影序列《没头脑和不高兴》处理效果对比图，其中(a)为原始图像序列帧，(b)为修复后序列帧。图4为动画老电影序列《没头脑和不高兴》的50帧RGB通道均值统计对比图，其中(a)为处理前统计图，(b)为处理后统计图。具体实施方式下面结合附图以及具体实施方式进一步说明本发明。如图1所示，一种老电影褪色闪烁的修复方法，包括如下步骤：101：记整个待处理的序列为P，对序列P进行分镜处理，记分镜个数为n，分镜序列为Pi；由于后期处理闪烁使用的颜色传递算法对图像的结构有一定要求，只有在参考图像帧与待处理帧之间的布局结构相似时，处理结果最为理想，故分镜是必须的步骤；102：对分镜序列Pi进行RGB均值统计，判断是否存在闪烁；闪烁具有单帧不易观测的特点，只有连续放映时，才能明显发现，利用RGB均值统计，能客观判断Pi是否存在闪烁；103：若判定Pi序列存在闪烁，则在Pi序列中筛选出一帧适合的参考帧；若无参考帧，则使用色彩校正算法产生参考帧；对于参考帧的选取具有主观因素影响，主要是通过修复人员判定是否在Pi序列存在褪色不明显的图像帧，在无法确定合适参考帧的情况下，在Pi序列中选择褪色现象较轻的一帧进行色彩校正；具体步骤为：1)将待处理帧的原始像素值范围映射到[0,1]，并使用公式(1)将其从RGB空间转换到lαβ空间；1αβ=0.57740.57740.57740.40820.4082-0.81650.7071-0.707100.38110.57830.04020.19670.72440.07820.02410.12880.8444RGB---(1)]]>2)使用公式(2)在lαβ空间进行色彩校正处理：{α′=α-mα-CαNβ′=β-mβ-CβN---(2)]]>其中，α，β表示lαβ颜色空间通道内所有像素的原始值；α′，β′表示处理好之后的所有像素值；mα，mβ表示α与β通道均值；Cα，Cβ表示实际需求的通道均值，若图像仅需要色彩校正，则Cα，Cβ均为0；3)使用公式(3)将处理好的参考帧转回RGB空间，并将其像素值映射到原始值范围：RGB=4.4679-3.58730.1193-1.21862.3809-0.16240.0497-0.24391.20450.57740.40820.70710.57740.4082-0.70710.5774-0.816501αβ---(3)]]>在进行映射前，利用公式(4)修正参考帧像素值的溢出，I(s){s=1,s>1s=0,s<0---(4)]]>其中，I表示参考帧图像，s表示其各通道像素值。104：使用颜色传递算法，对Pi序列的每一帧待修复图像进行处理；具体如下：1)将参考帧与待修复帧转到lαβ颜色空间，计算参考帧的平均值参考帧的方差待修复帧的平均值待修复帧的方差2)使用公式(5)弱化待修复帧的整体色彩信息；lS′=lS-mSlαS′=αS-mSαβS′=βS-mSβ---(5)]]>其中，lS，αS，βS是待修复帧三通道的所有像素值，lS′，α′S，βS′是待修复帧三通道弱化后的像素值；3)使用公式(6)将待修复帧与参考帧标准差之比作为通道值偏移量的系数，将参考帧的细节信息映射到待修复帧中；l′=σRlσSl×lS′α′=σRασSα×αS′β′=σRβσSβ×βS′---(6)]]>4)使用公式(7)将参考帧的整体信息添加到处理后的待修复帧中；l=l′+mRlα=α′+mRαβ=β′+mRβ---(7)]]>5)将处理好的图像帧转回RGB颜色空间。105：在Pi序列处理完成后，对下一组分镜序列P(i+1)重复步骤102-104处理过程，直至Pi＝Pn。下面结合具体的实例证明本发明提供的一种老电影闪烁修复方法的可行性，详见下文描述：本发明截取动画老电影《没头脑和不高兴》部分序列进行实验，每帧图像大小为2048×1556，分为4组镜头共578张图片。图2(a)取自第二组分镜序列第33帧，该序列共50帧图像，由于闪烁严重无法筛选出恰当参考帧，利用色彩校正算法对其进行校正，处理结果如图2(b)所示。对比原始图像，图2(a)的偏色现象得到明显改善。图3(a)取自第二组分镜序列第30帧至35帧，经过颜色传递算法处理后，效果如图3(b)所示，图3(a)中序列主色调逐渐变红的情况得到了极大改善，基本保持在同一色调上。图4是第二组序列的RGB均值统计，从处理前后对比来看，图4(b)的统计曲线趋近平滑，客观说明本发明的有效性。