CN105678713B - 基于线性扩散板与反射分量分离理论的高光去除方法 - Google Patents
基于线性扩散板与反射分量分离理论的高光去除方法 Download PDFInfo
- Publication number
- CN105678713B CN105678713B CN201610076239.6A CN201610076239A CN105678713B CN 105678713 B CN105678713 B CN 105678713B CN 201610076239 A CN201610076239 A CN 201610076239A CN 105678713 B CN105678713 B CN 105678713B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- pixel
- mover
- mfrac
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 24
- 238000000926 separation method Methods 0.000 title claims abstract description 20
- 238000004737 colorimetric analysis Methods 0.000 claims abstract description 6
- 238000005070 sampling Methods 0.000 claims description 15
- 239000003795 chemical substances by application Substances 0.000 claims description 12
- FGUUSXIOTUKUDN-IBGZPJMESA-N C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 Chemical compound C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 FGUUSXIOTUKUDN-IBGZPJMESA-N 0.000 claims description 8
- GNFTZDOKVXKIBK-UHFFFAOYSA-N 3-(2-methoxyethoxy)benzohydrazide Chemical compound COCCOC1=CC=CC(C(=O)NN)=C1 GNFTZDOKVXKIBK-UHFFFAOYSA-N 0.000 claims description 7
- 238000001228 spectrum Methods 0.000 claims description 6
- 238000005286 illumination Methods 0.000 claims description 3
- 230000035945 sensitivity Effects 0.000 claims description 3
- 239000013598 vector Substances 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 3
- 230000001629 suppression Effects 0.000 abstract description 3
- 240000002234 Allium sativum Species 0.000 description 9
- 235000004611 garlic Nutrition 0.000 description 9
- 239000004579 marble Substances 0.000 description 9
- 239000000919 ceramic Substances 0.000 description 6
- 230000005484 gravity Effects 0.000 description 6
- 238000005259 measurement Methods 0.000 description 6
- 239000000203 mixture Substances 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 239000002184 metal Substances 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 206010068052 Mosaicism Diseases 0.000 description 1
- 238000012356 Product development Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000011195 cermet Substances 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000000615 nonconductor Substances 0.000 description 1
- 230000010287 polarization Effects 0.000 description 1
- 238000003908 quality control method Methods 0.000 description 1
- 238000001028 reflection method Methods 0.000 description 1
- 210000003765 sex chromosome Anatomy 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T15/00—3D [Three Dimensional] image rendering
- G06T15/10—Geometric effects
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0004—Industrial image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30108—Industrial image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30108—Industrial image inspection
- G06T2207/30136—Metal
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30168—Image quality inspection
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Quality & Reliability (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Geometry (AREA)
- Computer Graphics (AREA)
- Image Processing (AREA)
- Spectrometry And Color Measurement (AREA)
Abstract
本发明属于三维结构光测量领域，具体涉及一种基于线性扩散板与反射分量分离理论的高光去除方法；该方法应用了线性扩散板与反射分量分离理论，其中，线性扩散板放置在投影仪与被测物之间，反射分量分离理论首先建立待处理图像信息模型，然后建立待处理图像漫反射与强反射色度模型，再依次建立归一化图像模型，建立非强反射图像模型，确定强反射像素点，最后处理强反射像素区域；本发明将线性扩散板与反射分量分离理论进行软硬有机结合，使二者彼此支持，实现线性扩散板或反射分量分离理论单独存在时，都不具有的抑制高光的作用，从而降低高光部分对三维重构的影响。
Description
技术领域
本发明属于三维结构光测量领域，具体涉及一种基于线性扩散板与反射分量分离理论的高光去除方法。
背景技术
上世纪70年代以来，编码光三维测量技术在高速检测、产品开发、质量控制、反求工程等领域得到广泛的应用和发展，其具有高精度、高效率、非接触等优点。在编码光三维测量中，需要在物体表面投射条纹，条纹携带了物体表面的三维信息，条纹的灰度变化反映了物体表面轮廓的变化。然而在实际工业测量中，条纹截面的灰度变化会受到很多因素的干扰，其中高光是影响最为强烈的一个因素。高光的存在不仅可能会使相机饱和，丢失条纹灰度变化信息，而且还将改变原有漫反射条纹的灰度分布，从而影响条纹中心提取的准确性。去高光问题目前已是利用光学方法测量金属陶瓷等强反射物体所面临的一个难点与共性问题。
现有的基于漫反射与镜面反射的分离高光抑制方法中，双色模型只适用于非导体材质，不适用于金属表面。极化方法在入射角接近90度时,金属表面的菲涅耳系数才有少许变化，易造成相机的饱和。而在避免和减少镜面反射方法中，单独使用多光源、多角度、多曝光、喷雾法、光度学等方法仍然存在着高光区域重叠部分不能完全去除、不能满足测量的精度需求、对金属表面造成腐蚀等问题。总之,在编码光三维测量领域，现有的方法仍然无法克服高光成分对三维重构的影响。
发明内容
针对上述问题，本发明公开了一种基于线性扩散板与反射分量分离理论的高光去除方法，该方法将线性扩散板与反射分量分离理论进行软硬有机结合，使二者彼此支持，实现线性扩散板或反射分量分离理论单独存在时，都不具有的抑制高光的作用，从而降低高光部分对三维重构的影响。
本发明的目的是这样实现的：
基于线性扩散板与反射分量分离理论的高光去除方法，包括以下步骤：
步骤a、在投影仪与被测物之间加入线性扩散板；
步骤b、建立待处理图像信息模型；
步骤c、建立待处理图像漫反射与强反射色度模型；
步骤d、建立归一化图像模型；
步骤e、建立非强反射图像模型；
步骤f、确定强反射像素点；
步骤g、处理强反射像素区域。
上述基于线性扩散板与反射分量分离理论的高光去除方法，步骤b到步骤g具体为：
步骤b、建立待处理图像信息模型
待处理图像表示为：
AR(X)＝α(X)∫ΩT(λ,X)B(λ)QR(λ)dλ+β(X)∫ΩB(λ)QR(λ)dλ
AG(X)＝α(X)∫ΩT(λ,X)B(λ)QG(λ)dλ+β(X)∫ΩB(λ)QG(λ)dλ
AB(X)＝α(X)∫ΩT(λ,X)B(λ)QB(λ)dλ+β(X)∫ΩB(λ)QB(λ)dλ
其中，A(X)为图像强度，α(X)为漫反射加权因子，β(X)为强反射加权因子，X＝{x,y}为像素点坐标，T(λ,X)为漫反射功率谱，B(λ)为强反射功率谱，Q(λ)为传感器灵敏度，下角标R、G、B分别表示图像的红色通道、绿色通道、以及蓝色通道；
令：
DR(X)＝∫ΩT(λ,X)B(λ)QR(λ)dλ
DG(X)＝∫ΩT(λ,X)B(λ)QG(λ)dλ
DB(X)＝∫ΩT(λ,X)B(λ)QB(λ)dλ
SR＝∫ΩB(λ)QR(λ)dλ
SG＝∫ΩB(λ)QG(λ)dλ
SB＝∫ΩB(λ)QB(λ)dλ
有：
AR(X)＝α(X)DR(X)+β(X)SR
AG(X)＝α(X)DG(X)+β(X)SG
AB(X)＝α(X)DB(X)+β(X)SB
步骤c、建立待处理图像漫反射与强反射色度模型
定义图像色度为：
在像素中只有漫反射的条件下，β(X)＝0，此像素的色度表达式为：
在像素中只有强反射的条件下，α(X)＝0，此像素的色度表达式为：
此时有：
AR(X)＝md(X)MR(X)+ms(X)KR
AG(X)＝md(X)MG(X)+ms(X)KG
AB(X)＝md(X)MB(X)+ms(X)KB
其中，md(X)MR(X)、md(X)MG(X)、md(X)MB(X)为漫反射分量，ms(X)KR、ms(X)KG、ms(X)KB为强反射分量，并且有：
md(X)＝α(X)[DR(X)+DG(X)+DB(X)]
ms(X)＝β(X)(SR+SG+SB)
步骤d、建立归一化图像模型
归一化的图像可以表示为：
AR′(X)＝AR(X)/KR′
A′G(X)＝AG(X)/K′G
AB′(X)＝AB(X)/KB′
其中，A′(X)为归一化的图像强度，K′为强反射分量的色度估计值；
步骤e、建立非强反射图像模型
定义经过颜色补偿之后的图像为非强反射图像，定义式如下：
其中，为非强反射图像的强度，为非强反射图像的漫反射色度，为非强反射图像的漫反射加权因子；
步骤f、确定强反射像素点
步骤f1、确定非漫反射像素
的像素为非漫反射像素，其中，d()为导数微分运算，log()为取对数运算；
步骤f2、判断以下条件是否同时满足：
如果满足，步骤f结束；
如果不满足，进入步骤f3；
步骤f3、判断以下条件是否满足：
max(KR(X)、KG(X)、KB(X))＝max(KR(X-1)、KG(X-1)、KB(X-1))
如果满足，则像素X为噪声像素；
如果不满足，则像素X为强反射像素；
按照此步骤将所有像素处理完毕，即可标记出被测图像中所有强反射像素区域Ω；
步骤g、处理强反射像素区域
步骤g1、根据强反射像素区域Ω，按照以下公式确定源区域Φ：
Φ＝A(X)-Ω
进而确定强反射像素区域Ω与源区域Φ的边界σΩ；
步骤g2、边界σΩ上的每个像素点o对应一个采样模块ψo，计算边界σΩ上每一个像素点o的优先级P(o)，对于区域ψo：其优先级可以表示为：
P(o)＝C(o)×D(o)
其中，C(o)为置信度，D(o)为数据项，并且有：
其中，|ψo|是采样窗口ψo的面积，α为规范化因子，⊥代表正交操作；no是点o正交于边界σΩ的单位正交向量，是o点处的照度；
步骤g3、寻找满足以下公式的像素：
其中，q均为像素点，为优先级最高的采样模块，表示优先级最高的采样模块与采样模块ψq的像素的平方差的和，为满足条件的源模块；在找到符合条件的源模块后，用源模块中的像素点按照对应关系代替模块的像素点；反复迭代，直到强反射像素区域Ω都处理完毕。
有益效果：本发明将线性扩散板与反射分量分离理论进行软硬有机结合，使二者彼此支持，实现线性扩散板或反射分量分离理论单独存在时，都不具有的抑制高光的作用，从而降低高光部分对三维重构的影响。
附图说明
图1是在投影仪与被测物之间加入线性扩散板的原理图。
图2是在投影仪与被测物之间加入线性扩散板的实物图。
图3是大理石捣蒜罐照片。
图4是未加入线性扩散板条件下的大理石捣蒜罐三维重构图。
图5是加入线性扩散板条件下的大理石捣蒜罐三维重构图。
图6是本发明方法得到的大理石捣蒜罐三维重构图。
具体实施方式
下面结合附图对本发明具体实施方式作进一步详细描述。
本实施例的基于线性扩散板与反射分量分离理论的高光去除方法，包括以下步骤：
步骤a、在投影仪与被测物之间加入线性扩散板，原理图如图1所示，实物图如图2所示；
步骤b、建立待处理图像信息模型
待处理图像表示为：
AR(X)＝α(X)∫ΩT(λ,X)B(λ)QR(λ)dλ+β(X)∫ΩB(λ)QR(λ)dλ
AG(X)＝α(X)∫ΩT(λ,X)B(λ)QG(λ)dλ+β(X)∫ΩB(λ)QG(λ)dλ
AB(X)＝α(X)∫ΩT(λ,X)B(λ)QB(λ)dλ+β(X)∫ΩB(λ)QB(λ)dλ
其中，A(X)为图像强度，α(X)为漫反射加权因子，β(X)为强反射加权因子，X＝{x,y}为像素点坐标，T(λ,X)为漫反射功率谱，B(λ)为强反射功率谱，Q(λ)为传感器灵敏度，下角标R、G、B分别表示图像的红色通道、绿色通道、以及蓝色通道；
令：
DR(X)＝∫ΩT(λ,X)B(λ)QR(λ)dλ
DG(X)＝∫ΩT(λ,X)B(λ)QG(λ)dλ
DB(X)＝∫ΩT(λ,X)B(λ)QB(λ)dλ
SR＝∫ΩB(λ)QR(λ)dλ
SG＝∫ΩB(λ)QG(λ)dλ
SB＝∫ΩB(λ)QB(λ)dλ
有：
AR(X)＝α(X)DR(X)+β(X)SR
AG(X)＝α(X)DG(X)+β(X)SG
AB(X)＝α(X)DB(X)+β(X)SB
步骤c、建立待处理图像漫反射与强反射色度模型
定义图像色度为：
在像素中只有漫反射的条件下，β(X)＝0，此像素的色度表达式为：
在像素中只有强反射的条件下，α(X)＝0，此像素的色度表达式为：
此时有：
AR(X)＝md(X)MR(X)+ms(X)KR
AG(X)＝md(X)MG(X)+ms(X)KG
AB(X)＝md(X)MB(X)+ms(X)KB
其中，md(X)MR(X)、md(X)MG(X)、md(X)MB(X)为漫反射分量，ms(X)KR、ms(X)KG、ms(X)KB为强反射分量，并且有：
md(X)＝α(X)[DR(X)+DG(X)+DB(X)]
ms(X)＝β(X)(SR+SG+SB)
步骤d、建立归一化图像模型
归一化的图像可以表示为：
AR′(X)＝AR(X)/KR′
A′G(X)＝AG(X)/K′G
AB′(X)＝AB(X)/KB′
其中，A′(X)为归一化的图像强度，K′为强反射分量的色度估计值；在本实施例中：
步骤e、建立非强反射图像模型
定义经过颜色补偿之后的图像为非强反射图像，定义式如下：
其中，为非强反射图像的强度，为非强反射图像的漫反射色度，为非强反射图像的漫反射加权因子；
步骤f、确定强反射像素点
步骤f1、确定非漫反射像素
的像素为非漫反射像素，其中，d()为导数微分运算，log()为取对数运算；
步骤f2、判断以下条件是否同时满足：
如果满足，步骤f结束；
如果不满足，进入步骤f3；
步骤f3、判断以下条件是否满足：
max(KR(X)、KG(X)、KB(X))＝max(KR(X-1)、KG(X-1)、KB(X-1))
如果满足，则像素X为噪声像素；
如果不满足，则像素X为强反射像素；
按照此步骤将所有像素处理完毕，即可标记出被测图像中所有强反射像素区域Ω；
步骤g、处理强反射像素区域
步骤g1、根据强反射像素区域Ω，按照以下公式确定源区域Φ：
Φ＝A(X)-Ω
进而确定强反射像素区域Ω与源区域Φ的边界σΩ；
步骤g2、边界σΩ上的每个像素点o对应一个采样模块ψo，计算边界σΩ上每一个像素点o的优先级P(o)，对于区域ψo：其优先级可以表示为：
P(o)＝C(o)×D(o)
其中，C(o)为置信度，D(o)为数据项，并且有：
其中，|ψo|是采样窗口ψo的面积，α为规范化因子，⊥代表正交操作；no是点o正交于边界σΩ的单位正交向量，是o点处的照度；
步骤g3、寻找满足以下公式的像素：
其中，q均为像素点，为优先级最高的采样模块，表示优先级最高的采样模块与采样模块ψq的像素的平方差的和，为满足条件的源模块；在找到符合条件的源模块后，用源模块中的像素点按照对应关系代替模块的像素点；反复迭代，直到强反射像素区域Ω都处理完毕。
为了验证本发明方法的有效性，在不同条件下对大理石捣蒜罐进行三维重构。其中，图3是大理石捣蒜罐照片，图4是未加入线性扩散板条件下的大理石捣蒜罐三维重构图，图5是加入线性扩散板条件下的大理石捣蒜罐三维重构图，图6是本发明方法得到的大理石捣蒜罐三维重构图。对比这几幅图可以发现，随着线性扩散板与反射分量分离理论的逐步施行，高光区域逐渐减小，最终实施本发明基于线性扩散板与反射分量分离理论的高光去除方法后，已经彻底消除高光对于三维重构的影响。
为了进一步验证本专利方法的有效性，分别对陶瓷瓶、陶瓷盘、陶瓷罐进行高光去除实验，实验结果下表所示：
从表中的结果可以看出，本发明基于线性扩散板与反射分量分离理论的高光去除方法实施后，图片中的高光像素数目大幅降低。对于陶瓷瓶、陶瓷盘、以及陶瓷罐，高光像素分别降低了约44倍、41倍、以及33倍，充分验证了本发明方法的有效性。
Claims (1)
1.基于线性扩散板与反射分量分离理论的高光去除方法，包括以下步骤：
步骤a、在投影仪与被测物之间加入线性扩散板；
步骤b、建立待处理图像信息模型；
步骤c、建立待处理图像漫反射与强反射色度模型；
步骤d、建立归一化图像模型；
步骤e、建立非强反射图像模型；
步骤f、确定强反射像素点；
步骤g、处理强反射像素区域；
其特征在于，步骤b到步骤g具体为：
步骤b、建立待处理图像信息模型
待处理图像表示为：
AR(X)＝α(X)∫ΩT(λ,X)B(λ)QR(λ)dλ+β(X)∫ΩB(λ)QR(λ)dλ
AG(X)＝α(X)∫ΩT(λ,X)B(λ)QG(λ)dλ+β(X)∫ΩB(λ)QG(λ)dλ
AB(X)＝α(X)∫ΩT(λ,X)B(λ)QB(λ)dλ+β(X)∫ΩB(λ)QB(λ)dλ
其中，A(X)为图像强度，α(X)为漫反射加权因子，β(X)为强反射加权因子，X＝{x,y}为像素点坐标，T(λ,X)为漫反射功率谱，B(λ)为强反射功率谱，Q(λ)为传感器灵敏度，下角标R、G、B分别表示图像的红色通道、绿色通道、以及蓝色通道；
令：
DR(X)＝∫ΩT(λ,X)B(λ)QR(λ)dλ
DG(X)＝∫ΩT(λ,X)B(λ)QG(λ)dλ
DB(X)＝∫ΩT(λ,X)B(λ)QB(λ)dλ
SR＝∫ΩB(λ)QR(λ)dλ
SG＝∫ΩB(λ)QG(λ)dλ
SB＝∫ΩB(λ)QB(λ)dλ
有：
AR(X)＝α(X)DR(X)+β(X)SR
AG(X)＝α(X)DG(X)+β(X)SG
AB(X)＝α(X)DB(X)+β(X)SB
步骤c、建立待处理图像漫反射与强反射色度模型
定义图像色度为：
<mrow>
<msub>
<mi>I</mi>
<mi>R</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>A</mi>
<mi>R</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>A</mi>
<mi>R</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>A</mi>
<mi>G</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>A</mi>
<mi>B</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>I</mi>
<mi>G</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>A</mi>
<mi>G</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>A</mi>
<mi>R</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>A</mi>
<mi>G</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>A</mi>
<mi>B</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>I</mi>
<mi>B</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>A</mi>
<mi>B</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>A</mi>
<mi>R</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>A</mi>
<mi>G</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>A</mi>
<mi>B</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
在像素中只有漫反射的条件下，β(X)＝0，此像素的色度表达式为：
<mrow>
<msub>
<mi>M</mi>
<mi>R</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>D</mi>
<mi>R</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>D</mi>
<mi>R</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>D</mi>
<mi>G</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>D</mi>
<mi>B</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>M</mi>
<mi>G</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>D</mi>
<mi>G</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>D</mi>
<mi>R</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>D</mi>
<mi>G</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>D</mi>
<mi>B</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>M</mi>
<mi>B</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>D</mi>
<mi>B</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>D</mi>
<mi>R</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>D</mi>
<mi>G</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>D</mi>
<mi>B</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
在像素中只有强反射的条件下，α(X)＝0，此像素的色度表达式为：
<mrow>
<msub>
<mi>K</mi>
<mi>R</mi>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>S</mi>
<mi>R</mi>
</msub>
<mrow>
<msub>
<mi>S</mi>
<mi>R</mi>
</msub>
<mo>+</mo>
<msub>
<mi>S</mi>
<mi>G</mi>
</msub>
<mo>+</mo>
<msub>
<mi>S</mi>
<mi>B</mi>
</msub>
</mrow>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>K</mi>
<mi>G</mi>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>S</mi>
<mi>G</mi>
</msub>
<mrow>
<msub>
<mi>S</mi>
<mi>R</mi>
</msub>
<mo>+</mo>
<msub>
<mi>S</mi>
<mi>G</mi>
</msub>
<mo>+</mo>
<msub>
<mi>S</mi>
<mi>B</mi>
</msub>
</mrow>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>K</mi>
<mi>B</mi>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>S</mi>
<mi>B</mi>
</msub>
<mrow>
<msub>
<mi>S</mi>
<mi>R</mi>
</msub>
<mo>+</mo>
<msub>
<mi>S</mi>
<mi>G</mi>
</msub>
<mo>+</mo>
<msub>
<mi>S</mi>
<mi>B</mi>
</msub>
</mrow>
</mfrac>
</mrow>
此时有：
AR(X)＝md(X)MR(X)+ms(X)KR
AG(X)＝md(X)MG(X)+ms(X)KG
AB(X)＝md(X)MB(X)+ms(X)KB
其中，md(X)MR(X)、md(X)MG(X)、md(X)MB(X)为漫反射分量，ms(X)KR、ms(X)KG、ms(X)KB为强反射分量，并且有：
md(X)＝α(X)[DR(X)+DG(X)+DB(X)]
ms(X)＝β(X)(SR+SG+SB)
步骤d、建立归一化图像模型
归一化的图像可以表示为：
A′R(X)＝AR(X)/K′R
A′G(X)＝AG(X)/K′G
A′B(X)＝AB(X)/K′B
其中，A′(X)为归一化的图像强度，K′为强反射分量的色度估计值；
步骤e、建立非强反射图像模型
定义经过颜色补偿之后的图像为非强反射图像，定义式如下：
<mrow>
<msub>
<mover>
<mi>A</mi>
<mo>~</mo>
</mover>
<mi>R</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mover>
<mi>m</mi>
<mo>~</mo>
</mover>
<mi>d</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
<msub>
<mover>
<mi>M</mi>
<mo>~</mo>
</mover>
<mi>R</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mover>
<mi>A</mi>
<mo>~</mo>
</mover>
<mi>G</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mover>
<mi>m</mi>
<mo>~</mo>
</mover>
<mi>d</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
<msub>
<mover>
<mi>M</mi>
<mo>~</mo>
</mover>
<mi>G</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mover>
<mi>A</mi>
<mo>~</mo>
</mover>
<mi>B</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mover>
<mi>m</mi>
<mo>~</mo>
</mover>
<mi>d</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
<msub>
<mover>
<mi>M</mi>
<mo>~</mo>
</mover>
<mi>B</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
</mrow>
其中，为非强反射图像的强度，为非强反射图像的漫反射色度，为非强反射图像的漫反射加权因子；
步骤f、确定强反射像素点
步骤f1、确定非漫反射像素
的像素为非漫反射像素，其中，d()为导数微分运算，log()为取对数运算；
步骤f2、判断以下条件是否同时满足：
<mrow>
<msub>
<mi>&Delta;</mi>
<mi>R</mi>
</msub>
<mo>=</mo>
<msubsup>
<mi>I</mi>
<mi>R</mi>
<mo>&prime;</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msubsup>
<mi>I</mi>
<mi>R</mi>
<mo>&prime;</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>></mo>
<mover>
<mi>R</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
<mrow>
<msub>
<mi>&Delta;</mi>
<mi>G</mi>
</msub>
<mo>=</mo>
<msubsup>
<mi>I</mi>
<mi>G</mi>
<mo>&prime;</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msubsup>
<mi>I</mi>
<mi>G</mi>
<mo>&prime;</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>></mo>
<mover>
<mi>G</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
<mrow>
<msub>
<mi>&Delta;</mi>
<mi>B</mi>
</msub>
<mo>=</mo>
<msubsup>
<mi>I</mi>
<mi>B</mi>
<mo>&prime;</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msubsup>
<mi>I</mi>
<mi>B</mi>
<mo>&prime;</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>X</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>></mo>
<mover>
<mi>B</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
如果满足，步骤f结束；
如果不满足，进入步骤f3；
步骤f3、判断以下条件是否满足：
max(KR(X)、KG(X)、KB(X))＝max(KR(X-1)、KG(X-1)、KB(X-1))
如果满足，则像素X为噪声像素；
如果不满足，则像素X为强反射像素；
按照此步骤将所有像素处理完毕，即可标记出被测图像中所有强反射像素区域Ω；
步骤g、处理强反射像素区域
步骤g1、根据强反射像素区域Ω，按照以下公式确定源区域Φ：
Φ＝A(X)-Ω
进而确定强反射像素区域Ω与源区域Φ的边界σΩ；
步骤g2、边界σΩ上的每个像素点o对应一个采样模块ψo，计算边界σΩ上每一个像素点o的优先级P(o)，对于区域ψo：其优先级可以表示为：
P(o)＝C(o)×D(o)
其中，C(o)为置信度，D(o)为数据项，并且有：
<mrow>
<mi>C</mi>
<mrow>
<mo>(</mo>
<mi>o</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&Sigma;</mi>
<mrow>
<mi>o</mi>
<mo>&Element;</mo>
<msub>
<mi>&psi;</mi>
<mi>o</mi>
</msub>
<mo>&cap;</mo>
<mi>&Omega;</mi>
</mrow>
</msub>
<mi>C</mi>
<mrow>
<mo>(</mo>
<mi>o</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>|</mo>
<msub>
<mi>&psi;</mi>
<mi>o</mi>
</msub>
<mo>|</mo>
</mrow>
</mfrac>
</mrow>
<mrow>
<mi>D</mi>
<mrow>
<mo>(</mo>
<mi>o</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mo>|</mo>
<mo>&dtri;</mo>
<msubsup>
<mi>I</mi>
<mi>o</mi>
<mo>&perp;</mo>
</msubsup>
<mo>&CenterDot;</mo>
<msub>
<mi>n</mi>
<mi>o</mi>
</msub>
<mo>|</mo>
</mrow>
<mi>&alpha;</mi>
</mfrac>
</mrow>
其中，|ψo|是采样窗口ψo的面积，α为规范化因子，⊥代表正交操作；no是点o正交于边界σΩ的单位正交向量，是o点处的照度；
步骤g3、寻找满足以下公式的像素：
<mrow>
<msub>
<mi>&psi;</mi>
<mover>
<mi>q</mi>
<mo>^</mo>
</mover>
</msub>
<mo>=</mo>
<mi>arg</mi>
<munder>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
<mrow>
<msub>
<mi>&psi;</mi>
<mi>q</mi>
</msub>
<mo>&Element;</mo>
<mi>&Phi;</mi>
</mrow>
</munder>
<mi>d</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&psi;</mi>
<mover>
<mi>p</mi>
<mo>^</mo>
</mover>
</msub>
<mo>,</mo>
<msub>
<mi>&psi;</mi>
<mi>q</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
其中，q均为像素点，为优先级最高的采样模块，表示优先级最高的采样模块与采样模块ψq的像素的平方差的和，为满足条件的源模块；在找到符合条件的源模块后，用源模块中的像素点按照对应关系代替模块的像素点；反复迭代，直到强反射像素区域Ω都处理完毕。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610076239.6A CN105678713B (zh) | 2016-02-03 | 2016-02-03 | 基于线性扩散板与反射分量分离理论的高光去除方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610076239.6A CN105678713B (zh) | 2016-02-03 | 2016-02-03 | 基于线性扩散板与反射分量分离理论的高光去除方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105678713A CN105678713A (zh) | 2016-06-15 |
| CN105678713B true CN105678713B (zh) | 2018-04-17 |
Family
ID=56303525
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610076239.6A Expired - Fee Related CN105678713B (zh) | 2016-02-03 | 2016-02-03 | 基于线性扩散板与反射分量分离理论的高光去除方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105678713B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108665421B (zh) * | 2017-03-31 | 2021-05-14 | 北京旷视科技有限公司 | 人脸图像高光分量去除装置和方法、存储介质产品 |
| CN110889810A (zh) * | 2019-11-29 | 2020-03-17 | 合肥视展光电科技有限公司 | 一种基于偏振的滤光薄膜提取图像方法及系统 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2015042939A1 (en) * | 2013-09-30 | 2015-04-02 | Carestream Health, Inc. | Method and system for intra-oral imaging using hdr imagine and highlight removal |
| CN104677309A (zh) * | 2015-03-24 | 2015-06-03 | 哈尔滨理工大学 | 基于稀疏表达的强反射表面编码光测量的高光抑制方法 |
| CN105023249A (zh) * | 2015-06-26 | 2015-11-04 | 清华大学深圳研究生院 | 基于光场的高光图像修复方法及装置 |
-
2016
- 2016-02-03 CN CN201610076239.6A patent/CN105678713B/zh not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2015042939A1 (en) * | 2013-09-30 | 2015-04-02 | Carestream Health, Inc. | Method and system for intra-oral imaging using hdr imagine and highlight removal |
| CN104677309A (zh) * | 2015-03-24 | 2015-06-03 | 哈尔滨理工大学 | 基于稀疏表达的强反射表面编码光测量的高光抑制方法 |
| CN105023249A (zh) * | 2015-06-26 | 2015-11-04 | 清华大学深圳研究生院 | 基于光场的高光图像修复方法及装置 |
Non-Patent Citations (3)
| Title |
|---|
| Real-time Specular Highlight Removal Using Bilateral Filtering;Qingxiong Yang 等;《ECCV"10 Proceedings of the 11th European conference on Computer vision:Part IV》;20100911;87-100 * |
| 基于单幅图像的高光分量可靠去除方法;郑芝寰 等;《杭州电子科技大学学报》;20110831;第31卷(第4期);引言，正文第1.1，1.3节 * |
| 自然光图像反射分量分离;张潇云 等;《湖南大学学报（自然科学版）》;20150228;第42卷(第2期);126-133 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105678713A (zh) | 2016-06-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105741249B (zh) | 一种强反射表面高光去除方法 | |
| Narasimhan et al. | Structured light in scattering media | |
| WO2017164971A3 (en) | Spatial relationships for integration of visual images of physical environment into virtual reality | |
| CN105678713B (zh) | 基于线性扩散板与反射分量分离理论的高光去除方法 | |
| CN106936964B (zh) | 一种基于霍夫变换模板匹配的手机屏幕角点检测方法 | |
| CN102063706A (zh) | 一种快速去雾方法 | |
| CN101017567A (zh) | 用于将数字照片转换成人像素描画的方法和设备 | |
| CN110288561A (zh) | 基于频率域滤波增强的耐火砖表面划痕识别方法 | |
| US10643376B2 (en) | Computer system and method for improved gloss representation in digital images | |
| CN106408526B (zh) | 一种基于多层矢量图的能见度检测方法 | |
| JP2017527863A5 (zh) | ||
| CN106504230B (zh) | 基于相位一致性的全参考彩色屏幕图像质量评估方法 | |
| TW201335892A (zh) | 特殊邊界量測系統及方法 | |
| CN107462182B (zh) | 一种基于机器视觉和红线激光器的截面轮廓形变检测方法 | |
| Mahiddine et al. | Underwater image preprocessing for automated photogrammetry in high turbidity water: An application on the Arles-Rhone XIII roman wreck in the Rhodano river, France | |
| Pintus et al. | State‐of‐the‐art in Multi‐Light Image Collections for surface visualization and analysis | |
| Valente et al. | Methods for ancient wall graffiti documentation: overview and applications | |
| CN107481201A (zh) | 一种基于多视角图像特征匹配的高光去除方法 | |
| CN105488475A (zh) | 手机中人脸检测方法 | |
| CN105488784A (zh) | 一种人像自动抠图方法 | |
| TW201814240A (zh) | 用於平面拍攝的長度量測方法 | |
| CN101930600B (zh) | 一种基于复合二阶分数阶次信号处理的边缘检测方法 | |
| Kotoula | Reflectance transformation imaging beyond the visible: ultraviolet reflected and ultraviolet induced visible fluorescence | |
| CN104392417A (zh) | 基于像素暗通道和各向异性扩散滤波的图像去雾方法 | |
| CN104966273B (zh) | 适用于光学遥感影像的dcm-htm去雾霾方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180417 Termination date: 20200203 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |