本发明涉及一种颜色补偿映射域的期望可视图像归约方法，其特征在于：预投影图像的像素分辨率为，其中：预投影图像的横向像素分辨率、预投影图像的纵向像素分辨率，能够准确地得到应用于预投影图像的颜色变换参数，有效地优化了期望可视图像的全局表现和投影显示画面的量化调控；能够根据预投影图像的颜色分布情况不同，动态地对预投影图像进行颜色变换而生成最佳的期望可视图像。
1.颜色补偿映射域的期望可视图像归约方法，其特征在于：预投影图像GP的像素分辨率为Q×P，其中：预投影图像GP的横向像素分辨率Q＝1024、预投影图像GP的纵向像素分辨率P＝768，该方法的具体实现步骤如下：步骤S1、设置投影仪-摄像机系统的图像采集像素分辨率为N×M，此处：图像采集的横向像素分辨率N＝5616、图像采集的纵向像素分辨率M＝3744，使用投影仪-摄像机系统拍摄一幅投影仪工作状态时的工作图像GW，工作图像GW的有效投影区域记为PZ，工作图像GW中所有像素点的亮度映射曲线构成了离散映射集合MS；
步骤S2、以摄像机采集的有效投影区域PZ为底版，绘制有效投影区域PZ的单通道掩模图像GM，单通道掩模图像GM的像素分辨率为N×M，单通道掩模图像GM中第i行第j列像素点的取值GM(i,j)如表达式EQ1所示；
步骤S3、根据单通道掩模图像GM中第i行第j列像素点的取值GM(i,j)，将离散映射集合MS中无用像素点对应的亮度映射曲线过滤掉，得到过滤后的离散映射集合MF，离散映射集合MF中第i行第j列像素点颜色通道c对应的亮度映射曲线MF(i,j,c)如表达式EQ2所示；
其中：aijc、bijc、dijc、eijc均为第i行第j列像素点颜色通道c对应的亮度映射曲线MF(i,j,c)表达式的已知计算系数，Sijc为投影仪投射图像中第i行第j列像素点颜色通道c对应的亮度值，Rijc(Sijc)为摄像机拍摄图像中第i行第j列像素点颜色通道c对应的亮度值，并且i∈[0,M-1]，j∈[0,N-1]，c∈{red green blue}，Sijc∈[0,255]，Rijc(Sijc)∈[0,255]；
步骤S4、设定离散映射集合MF中亮度映射曲线MF(i,j,c)的响应区间阈值ε为180，以摄像机拍摄图像中第i行第j列像素点的颜色通道c对应的亮度差值变量DRijc＝Rijc(255)-Rijc(0)为限定条件，对离散映射集合MF中存在异常的亮度映射曲线进行清洗，则可得到异常数据清洗后的离散映射集合MW，离散映射集合MW中第i行第j列像素点颜色通道c对应的亮度映射曲线MW(i,j,c)如表达式EQ3所示；
表达式EQ3：WRijc(WSijc)＝aijc×arctan(bijc×WSijc+dijc)+eijc DRijc≥ε其中：aijc、bijc、dijc、eijc均为第i行第j列像素点颜色通道c对应的亮度映射曲线MF(i,j,c)表达式的已知计算系数，WSijc为投影仪投射图像中第i行第j列像素点颜色通道c对应的亮度值，WRijc(WSijc)为摄像机拍摄图像中第i行第j列像素点颜色通道c对应的亮度值，i∈[0,M-1]，j∈[0,N-1]，c∈{red green blue}，WSijc∈[0,255]，WRijc(WSijc)∈[0,255]；
步骤S5、根据亮度映射曲线MW(i,j,c)的表达式EQ3，利用WRijc(0)计算得到离散映射集合MW中第i行第j列像素点颜色通道c对应的最小响应值NSijc，利用WRijc(255)计算得到离散映射集合MW中第i行第j列像素点颜色通道c对应的最大响应值XSijc，其中：i∈[0,M-1]，j∈[0,N-1]，c∈{red green blue}，所有最小响应值NSijc构成最小响应值序列NS，所有最大响应值XSijc构成最大响应值序列XS；
步骤S6、将预投影图像GP以填充的方式叠加于有效投影区域PZ之上，则形成一幅模拟图像GS，模拟图像GS的像素分辨率为N×M，模拟图像GS中第i行第j列像素点颜色通道c的取值GS(i,j,c)如表达式EQ4所示，其中：i∈[0,M-1]，j∈[0,N-1]，c∈{red green blue}；
步骤S7、令最小响应值序列NS的最小值fsmin＝min(NS)，令最小响应值序列NS的最大值fsmax＝max(NS)，最大响应值序列XS的最小值csmin＝min(XS)，最大响应值序列XS的最大值csmax＝max(XS)，定义变量：数据量tfn、数据量tcx、全局偏移量tdf、全局压缩比trc，并且tfn＝fsmin、tcx＝csmax、trc＝1.0；
步骤S8、令全局偏移量tdf＝fsmax-tfn、临时数据量ptcx＝tcx×trc，则可对数据量tcx进行值更新，使得tcx＝ptcx+tdf，并利用表达式EQ5对全局压缩比trc的值进行更新赋值；
步骤S9、利用表达式EQ5对模拟图像GS的每一个像素点进行全局映射偏移与压缩处理，模拟图像GS中第i行第j列像素点颜色通道c的取值GS(i,j,c)如表达式EQ6所示，其中：i∈[0,M-1]，j∈[0,N-1]，c∈{red green blue}；
表达式EQ6：GS(i,j,c)＝(GS(i,j,c)+tdf)×trc
步骤S10、令临时数据量ptfn＝(tfn+tdf)×trc，则可对全局偏移量tdf和数据量tfn进行更新赋值，使得tdf＝fsmax-ptfn、tfn＝ptfn，如果全局偏移量tdf＞5则回退至步骤S8继续处理，否则结束全局映射偏移与压缩处理，最后得到符合颜色补偿映射域的期望可视图像。
颜色补偿映射域的期望可视图像归约方法
技术领域
[0001] 本发明涉及一种颜色补偿映射域的期望可视图像归约方法，属于增强现实技术和投影显示技术领域。
背景技术
[0002] 传统普适性的投影接收表面具有漫反射或近似漫反射的物理特性，它是具有统一反射率的标准化理想表面。顺应数据可视化、增强现实、互动影视等领域的突进式发展，为了呈现出新奇、炫动、震撼、融入感的视觉画面和特殊显示效果，多种构成形式的智能型投影系统被使用的越来越具针对性。由于投影表面材质、纹理、形状和外界环境因素等物理属性会产生变化，因此必将引起投影显示画面的几何畸变和颜色失真，二者共同决定了人眼对投影显示画面的视觉感受。
[0003] 关于投影显示画面几何畸变方面的校正技术研究，已经存在多种有效的解决方案能够保障投影显示画面几何校正的可用性，这方面取得了相当数量的重要成果并显见成效，其中：既存在依靠大量人工交互调节的几何校正实现方式，又存在无需人工交互调节的自主几何校正实现方式。相对于投影显示画面的几何畸变而言，人类视觉特性决定了人眼对投影显示画面颜色失真的敏感性更高，使之成为投影显示技术领域的重点研究问题之一。关于投影显示画面颜色失真方面的补偿技术研究，虽然领域内开展了大量的、细致的、极具针对性的补偿技术研究工作，但是由于颜色失真的补偿过程严重受限于投影系统所处的现实环境和摄像系统的颜色复原能力，因此，不同工作环境中投影显示画面补偿后的显示结果也会产生色度、亮度、对比度等方面的视觉差异，甚至某些情况会导致颜色补偿过程完全失效而无法修正投影显示画面的颜色失真。投影仪-摄像机(ProCam)系统的颜色转换与传递过程是复杂的，它不仅与投影仪和摄像机的光学特性相关，更与现实环境的光照、硬件工作参数的设定、光线路径的中心指向等息息相关。
[0004] 预投影图像颜色补偿的投影显示画面是提供给观察者的，而观察者对投影显示画面的主观视觉感受却无法直接获得，故需借助摄像机间接地模拟或替代观察者的主观视觉响应，并以此为依据具体量化如何补偿预投影图像。利用摄像机拍摄的多幅亮度标定图像，可以拟合得到投影仪与摄像机之间的亮度映射曲线，由每个像素点亮度响应曲线构成的离散映射集合仅反映了预投影图像与投影显示画面的间接变换关系，而无法保证任意预投影图像都能得到有效的补偿。
发明内容
[0005] 本发明的目的在于提供一种颜色补偿映射域的期望可视图像归约方法，旨在解决期望可视图像所对应预投影图像的颜色补偿失效问题。通过对预投影图像进行颜色补偿映射域的压缩和剔除处理，确保颜色补偿过程中所有像素点的预补偿颜色值都能够得到有效计算，而不会产生颜色溢出或存在无效颜色补偿值等情况。
[0006] 本发明的技术方案是这样实现的：颜色补偿映射域的期望可视图像归约方法，其特征在于：预投影图像GP的像素分辨率为Q×P，其中：预投影图像GP的横向像素分辨率Q＝1024、预投影图像GP的纵向像素分辨率P＝768，该方法的具体实现步骤如下：
[0007] 步骤S1、设置投影仪-摄像机系统的图像采集像素分辨率为N×M，此处：图像采集的横向像素分辨率N＝5616、图像采集的纵向像素分辨率M＝3744，使用投影仪-摄像机系统拍摄一幅投影仪工作状态时的工作图像GW，工作图像GW的有效投影区域记为PZ，工作图像GW中所有像素点的亮度映射曲线构成了离散映射集合MS；
[0008] 步骤S2、以摄像机采集的有效投影区域PZ为底版，绘制有效投影区域PZ的单通道掩模图像GM，单通道掩模图像GM的像素分辨率为N×M，单通道掩模图像GM中第i行第j列像素点的取值GM(i,j)如表达式EQ1所示；
[0009] 表达式EQ1：
[0010] 步骤S3、根据单通道掩模图像GM中第i行第j列像素点的取值GM(i,j)，将离散映射集合MS中无用像素点对应的亮度映射曲线过滤掉，得到过滤后的离散映射集合MF，离散映射集合MF中第i行第j列像素点颜色通道c对应的亮度映射曲线MF(i,j,c)如表达式EQ2所示；
[0011] 表达式EQ2：
[0012] 其中：aijc、bijc、dijc、eijc均为第i行第j列像素点颜色通道c对应的亮度映射曲线MF(i,j,c)表达式的已知计算系数，Sijc为投影仪投射图像中第i行第j列像素点颜色通道c对应的亮度值，Rijc(Sijc)为摄像机拍摄图像中第i行第j列像素点颜色通道c对应的亮度值，并且i∈[0,M-1]，j∈[0,N-1]，c∈{red green blue}，Sijc∈[0,255]，Rijc(Sijc)∈[0,255]；
[0013] 步骤S4、设定离散映射集合MF中亮度映射曲线MF(i,j,c)的响应区间阈值ε为180，以DRijc＝Rijc(255)-Rijc(0)为限定条件，对离散映射集合MF中存在异常的亮度映射曲线进行清洗，则可得到异常数据清洗后的离散映射集合MW，离散映射集合MW中第i行第j列像素点颜色通道c对应的亮度映射曲线MW(i,j,c)如表达式EQ3所示；
[0014] 表达式EQ3：WRijc(WSijc)＝aijc×arctan(bijc×WSijc+dijc)+eijc DRijc≥ε[0015] 其中：aijc、bijc、dijc、eijc均为第i行第j列像素点颜色通道c对应的亮度映射曲线MF(i,j,c)表达式的已知计算系数，WSijc为投影仪投射图像中第i行第j列像素点颜色通道c对应的亮度值，WRijc(WSijc)为摄像机拍摄图像中第i行第j列像素点颜色通道c对应的亮度值，i∈[0,M-1]，j∈[0,N-1]，c∈{red green blue}，WSijc∈[0,255]，WRijc(WSijc)∈[0,255]；
[0016] 步骤S5、根据亮度映射曲线MW(i,j,c)的表达式EQ3，利用WRijc(0)计算得到离散映射集合MW中第i行第j列像素点颜色通道c对应的最小响应值NSijc，利用WRijc(255)计算得到离散映射集合MW中第i行第j列像素点颜色通道c对应的最大响应值XSijc，其中：i∈[0,M-1]，j∈[0,N-1]，c∈{red green blue}，所有最小响应值NSijc构成最小响应值序列NS，所有最大响应值XSijc构成最大响应值序列XS；
[0017] 步骤S6、将预投影图像GP以填充的方式叠加于有效投影区域PZ之上，则形成一幅模拟图像GS，模拟图像GS的像素分辨率为N×M，模拟图像GS中第i行第j列像素点颜色通道c的取值GS(i,j,c)如表达式EQ4所示，其中：i∈[0,M-1]，j∈[0,N-1]，c∈{red green blue}；
[0018] 表达式EQ4：
[0019] 步骤S7、令最小响应值序列NS的最小值fsmin＝min(NS)，令最小响应值序列NS的最大值fsmax＝max(NS)，最大响应值序列XS的最小值csmin＝min(XS)，最大响应值序列XS的最大值csmax＝max(XS)，定义变量：数据量tfn、数据量tcx、全局偏移量tdf、全局压缩比trc，并且tfn＝fmsin、tcx＝csmax、trc＝1.0；
[0020] 步骤S8、令全局偏移量tdf＝fsmax-tfn、数据量tcx＝tcx×trc，使得数据量tcx＝tcx+tdf，则得到如表达式EQ5所示的全局压缩比trc值；
[0021] 表达式EQ5：
[0022] 步骤S9、利用表达式EQ5对模拟图像GS的每一个像素点进行全局映射偏移与压缩处理，模拟图像GS中第i行第j列像素点颜色通道c的取值GS(i,j,c)如表达式EQ6所示，其中：i∈[0,M-1]，j∈[0,N-1]，c∈{red green blue}；
[0023] 表达式EQ6：GS(i，j，c)＝(GS(i，j，c)+tdf)×trc
[0024] 步骤S10、令数据量tfn＝(tfn+tdf)×trc，则全局偏移量tdf＝fsmax-tfn，如果全局偏移量tdf＞5则回退至步骤S8继续处理，否则结束全局映射偏移与压缩处理，最后得到符合颜色补偿映射域的期望可视图像。
[0025] 通过使用以上步骤可以得出离散映射集合的多个统计特性，根据这些统计特性不仅可以明确判断预投影图像是否符合颜色补偿映射域的取值限定要求，还能够利用这些统计结果对不符合颜色补偿映射域限定要求的预投影图像进行颜色变换，进而得到视觉感官最为理想的期望可视图像。
[0026] 本发明的积极效果在于提供了一种颜色补偿映射域的期望可视图像归约方法，该方法通过对离散映射集合进行统计与分析，能够准确地得到应用于预投影图像的颜色变换参数，利用这些具有指导意义的辅助数据完成预投影图像的颜色变换，有效地优化了期望可视图像的全局表现和投影显示画面的量化调控；能够根据预投影图像的颜色分布情况不同，动态地对预投影图像进行颜色变换而生成最佳的期望可视图像。
附图说明
[0027] 图1为本发明中期望可视图像归约方法的工作流程图，其中：“异常数据清洗”对应于具体实现步骤中的S2～S4；“全局映射偏移与压缩”对应于具体实现步骤中的S5～S10。
[0028] 图2为本发明中投影仪-摄像机系统工作状态时的工作图像和有效投影区域示意图。
具体实施方式
[0029] 下面结合附图和实施例对本发明做进一步的描述；图1为期望可视图像归约方法的工作流程图，其中：“异常数据清洗”对应于具体实现步骤中的S2～S4；“全局映射偏移与压缩”对应于具体实现步骤中的S5～S10；图2为投影仪-摄像机系统工作状态时的工作图像和有效投影区域示意图。颜色补偿映射域的期望可视图像归约方法，其特征在于:预投影图像GP的像素分辨率为Q×P，其中：预投影图像GP的横向像素分辨率Q＝1024、预投影图像GP的纵向像素分辨率P＝768，该方法的具体实现步骤如下：
[0030] 步骤S1、设置投影仪-摄像机系统的图像采集像素分辨率为N×M，此处：图像采集的横向像素分辨率N＝5616、图像采集的纵向像素分辨率M＝3744，使用投影仪-摄像机系统拍摄一幅投影仪工作状态时的工作图像GW，工作图像GW的有效投影区域记为PZ，工作图像GW中所有像素点的亮度映射曲线构成了离散映射集合MS；
[0031] 步骤S2、以摄像机采集的有效投影区域PZ为底版，绘制有效投影区域PZ的单通道掩模图像GM，单通道掩模图像GM的像素分辨率为N×M，单通道掩模图像GM中第i行第j列像素点的取值GM(i,j)如表达式EQ1所示；
[0032] 表达式EQ1：
[0033] 步骤S3、根据单通道掩模图像GM中第i行第j列像素点的取值GM(i,j)，将离散映射集合MS中无用像素点对应的亮度映射曲线过滤掉，得到过滤后的离散映射集合MF，离散映射集合MF中第i行第j列像素点颜色通道c对应的亮度映射曲线MF(i,j,c)如表达式EQ2所示；
[0034] 表达式EQ2：
[0035] 其中：aijc、bijc、dijc、eijc均为第i行第j列像素点颜色通道c对应的亮度映射曲线MF(i,j,c)表达式的已知计算系数，Sijc为投影仪投射图像中第i行第j列像素点颜色通道c对应的亮度值，Rijc(Sijc)为摄像机拍摄图像中第i行第j列像素点颜色通道c对应的亮度值，并且i∈[0,M-1]，j∈[0,N-1]，c∈{red green blue}，Sijc∈[0,255]，Rijc(Sijc)∈[0,255]；
[0036] 步骤S4、设定离散映射集合MF中亮度映射曲线MF(i,j,c)的响应区间阈值ε为180，以DRijc＝Rijc(255)-Rijc(0)为限定条件，对离散映射集合MF中存在异常的亮度映射曲线进行清洗，则可得到异常数据清洗后的离散映射集合MW，离散映射集合MW中第i行第j列像素点颜色通道c对应的亮度映射曲线MW(i,j,c)如表达式EQ3所示；
[0037] 表达式EQ3：WRijc(WSijc)＝aijc×arctan(bijc×WSijc+dijc)+eijc DRijc≥ε[0038] 其中：aijc、bijc、dijc、eijc均为第i行第j列像素点颜色通道c对应的亮度映射曲线MF(i,j,c)表达式的已知计算系数，WSijc为投影仪投射图像中第i行第j列像素点颜色通道c对应的亮度值，WRijc(WSijc)为摄像机拍摄图像中第i行第j列像素点颜色通道c对应的亮度值，i∈[0,M-1]，j∈[0,N-1]，c∈{red green blue}，WSijc∈[0,255]，WRijc(WSijc)∈[0,255]；
[0039] 步骤S5、根据亮度映射曲线MW(i,j,c)的表达式EQ3，利用WRijc(0)计算得到离散映射集合MW中第i行第j列像素点颜色通道c对应的最小响应值NSijc，利用WRijc(255)计算得到离散映射集合MW中第i行第j列像素点颜色通道c对应的最大响应值XSijc，其中：i∈[0,M-1]，j∈[0,N-1]，c∈{red green blue}，所有最小响应值NSijc构成最小响应值序列NS，所有最大响应值XSijc构成最大响应值序列XS；
[0040] 步骤S6、将预投影图像GP以填充的方式叠加于有效投影区域PZ之上，则形成一幅模拟图像GS，模拟图像GS的像素分辨率为N×M，模拟图像GS中第i行第j列像素点颜色通道c的取值GS(i,j,c)如表达式EQ4所示，其中：i∈[0,M-1]，j∈[0,N-1]，c∈{red green blue}；
[0041] 表达式EQ4：
[0042] 步骤S7、令最小响应值序列NS的最小值fsmin＝min(NS)，令最小响应值序列NS的最大值fsmax＝max(NS)，最大响应值序列XS的最小值csmin＝min(XS)，最大响应值序列XS的最大值csmax＝max(XS)，定义变量：数据量tfn、数据量tcx、全局偏移量tdf、全局压缩比trc，并且tfn＝fmsin、tcx＝csmax、trc＝1.0；
[0043] 步骤S8、令全局偏移量tdf＝fsmax-tfn、数据量tcx＝tcx×trc，使得数据量tcx＝tcx+tdf，则得到如表达式EQ5所示的全局压缩比trc值；
[0044] 表达式EQ5：
[0045] 步骤S9、利用表达式EQ5对模拟图像GS的每一个像素点进行全局映射偏移与压缩处理，模拟图像GS中第i行第j列像素点颜色通道c的取值GS(i,j,c)如表达式EQ6所示，其中：i∈[0,M-1]，j∈[0,N-1]，c∈{red green blue}；
[0046] 表达式EQ6：GS(i,j,c)＝(GS(i,j,c)+tdf)×trc
[0047] 步骤S10、令数据量tfn＝(tfn+tdf)×trc，则全局偏移量tdf＝fsmax-tfn，如果全局偏移量tdf＞5则回退至步骤S8继续处理，否则结束全局映射偏移与压缩处理，最后得到符合颜色补偿映射域的期望可视图像。
[0048] 通过使用以上步骤可以得出离散映射集合的多个统计特性，根据这些统计特性不仅可以明确判断预投影图像是否符合颜色补偿映射域的取值限定要求，还能够利用这些统计结果对不符合颜色补偿映射域限定要求的预投影图像进行颜色变换，进而得到视觉感官最为理想的期望可视图像。
