发明名称 一种投影仪伽马校正方法
摘要 本发明公开了一种投影仪伽马校正方法。首先搭建硬件校正系统,利用计算机生成两幅不同的单一灰度值图像,经投影仪投射到表面平整的平面,同时相机采集平面上的图像并将其传送到计算机中进行图像处理;然后在频域分析图像的分布特征,进而提取特征点求得投影仪伽马,最后将投影仪伽马预编码到正弦光栅条纹中,经投影仪投射即可得到保持良好正弦性的光栅条纹。该方法操作简单,由于只处理两幅图像,在保证精度的前提下,可快速实现伽马校正,极大提高了系统的整体测量速度。
申请公布号 CN105491315A 申请公布日期 2016.04.13
申请号 CN201610109658.5 申请日期 2016.02.29
申请人 河北工业大学 发明人 戴士杰;付津昇;周国霞;宋涛
分类号 H04N5/74(2006.01)I;H04N17/00(2006.01)I 主分类号 H04N5/74(2006.01)I
代理机构 天津翰林知识产权代理事务所(普通合伙) 12210 代理人 李济群;王瑞
主权项 一种投影仪伽马校正方法,包括以下步骤:1)搭建硬件校正系统:相机、投影仪和平面形成三角关系;投影仪的视场覆盖整个平面;平面占相机的60%‑100%视场;2)计算机生成两幅不同的单一灰度值图像:同一幅图像上每个像素的灰度值都相等,而且两幅图像间的像素灰度值不同,将两幅图像的灰度值分别用I<sub>1</sub>和I<sub>2</sub>表示,所述灰度值不小于0.3;3)投影图像与采集图像:将计算机生成的两幅灰度图像经投影仪投射到表面平整的平面上;同时相机采集平面上的图像并将其传送到计算机中进行图像处理;相机采集的图像可以表示为<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msubsup><mi>I</mi><mn>1</mn><mi>C</mi></msubsup><mo>=</mo><mi>&alpha;</mi><msup><mrow><mo>&lsqb;</mo><msub><mi>I</mi><mn>1</mn></msub><mo>&rsqb;</mo></mrow><mi>&gamma;</mi></msup><mo>&CircleTimes;</mo><mi>g</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000931212420000011.GIF" wi="1509" he="95" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msubsup><mi>I</mi><mn>2</mn><mi>C</mi></msubsup><mo>=</mo><mi>&alpha;</mi><msup><mrow><mo>&lsqb;</mo><msub><mi>I</mi><mn>2</mn></msub><mo>&rsqb;</mo></mrow><mi>&gamma;</mi></msup><mo>&CircleTimes;</mo><mi>g</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000931212420000012.GIF" wi="1510" he="87" /></maths>式中:<img file="FDA0000931212420000013.GIF" wi="69" he="78" />表示相机采集的第一幅图像,<img file="FDA0000931212420000014.GIF" wi="67" he="86" />为相机采集的第二幅图像,α∈[0,1]为待测物体表面反射率,γ为投影仪伽马,<img file="FDA0000931212420000015.GIF" wi="68" he="63" />表示卷积运算,[I<sub>1</sub>]<sup>γ</sup>表示I<sub>1</sub>的γ次幂,[I<sub>2</sub>]<sup>γ</sup>表示I<sub>2</sub>的γ次幂,g(x,y)表示系统散焦作用的点分布函数,即g(x,y)=(1/2πσ<sup>2</sup>)exp[‑(x<sup>2</sup>+y<sup>2</sup>)/2σ<sup>2</sup>]   (3)式中:σ为系统的散焦系数;g(x,y)的傅里叶变换式表示系统的光学传递函数,即<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>G</mi><mrow><mo>(</mo><msub><mi>f</mi><mi>x</mi></msub><mo>,</mo><msub><mi>f</mi><mi>y</mi></msub><mo>)</mo></mrow><mo>=</mo><mi>F</mi><mo>{</mo><mi>g</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>}</mo><mo>=</mo><mi>exp</mi><mo>&lsqb;</mo><mo>-</mo><mn>2</mn><msup><mi>&pi;</mi><mn>2</mn></msup><msup><mi>&sigma;</mi><mn>2</mn></msup><mrow><mo>(</mo><msubsup><mi>f</mi><mi>x</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>f</mi><mi>y</mi><mn>2</mn></msubsup><mo>)</mo></mrow><mo>&rsqb;</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000931212420000016.GIF" wi="1542" he="95" /></maths>式中F{}为傅里叶变换,f<sub>x</sub>和f<sub>y</sub>分别对应于变量x和y的频域变量,即为沿x轴和y轴的空间频率;4)频域变换:将等式(1)、(2)两边同时进行傅里叶变换<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msubsup><mi>F</mi><mn>1</mn><mi>C</mi></msubsup><mrow><mo>(</mo><msub><mi>f</mi><mi>x</mi></msub><mo>,</mo><msub><mi>f</mi><mi>y</mi></msub><mo>)</mo></mrow><mo>=</mo><msub><mi>F</mi><mn>1</mn></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>x</mi></msub><mo>,</mo><msub><mi>f</mi><mi>y</mi></msub><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>G</mi><mrow><mo>(</mo><msub><mi>f</mi><mi>x</mi></msub><mo>,</mo><msub><mi>f</mi><mi>y</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000931212420000021.GIF" wi="1542" he="87" /></maths><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msubsup><mi>F</mi><mn>2</mn><mi>C</mi></msubsup><mrow><mo>(</mo><msub><mi>f</mi><mi>x</mi></msub><mo>,</mo><msub><mi>f</mi><mi>y</mi></msub><mo>)</mo></mrow><mo>=</mo><msub><mi>F</mi><mn>2</mn></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>x</mi></msub><mo>,</mo><msub><mi>f</mi><mi>y</mi></msub><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>G</mi><mrow><mo>(</mo><msub><mi>f</mi><mi>x</mi></msub><mo>,</mo><msub><mi>f</mi><mi>y</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000931212420000022.GIF" wi="1550" he="85" /></maths>式中<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>F</mi><mn>1</mn></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>x</mi></msub><mo>,</mo><msub><mi>f</mi><mi>y</mi></msub><mo>)</mo></mrow><mo>=</mo><mi>F</mi><mo>{</mo><msubsup><mi>&alpha;I</mi><mn>1</mn><mi>&gamma;</mi></msubsup><mo>}</mo><mo>,</mo><msub><mi>F</mi><mn>2</mn></msub><mrow><mo>(</mo><msub><mi>f</mi><mi>x</mi></msub><mo>,</mo><msub><mi>f</mi><mi>y</mi></msub><mo>)</mo></mrow><mo>=</mo><mi>F</mi><mo>{</mo><msubsup><mi>&alpha;I</mi><mn>2</mn><mi>&gamma;</mi></msubsup><mo>}</mo></mrow>]]></math><img file="FDA0000931212420000023.GIF" wi="895" he="94" /></maths>5)特征点提取:通过对采集图像的频域分析,提取图像在频域内的特征点,求取投影仪伽马<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mi>&gamma;</mi><mo>=</mo><msubsup><mi>log</mi><mrow><mo>(</mo><mfrac><msub><mi>I</mi><mn>1</mn></msub><msub><mi>T</mi><mn>2</mn></msub></mfrac><mo>)</mo></mrow><mrow><mo>(</mo><mfrac><mrow><msubsup><mi>F</mi><mn>1</mn><mi>C</mi></msubsup><mrow><mo>(</mo><mn>0</mn><mo>,</mo><mn>0</mn><mo>)</mo></mrow></mrow><mrow><msubsup><mi>F</mi><mn>2</mn><mi>C</mi></msubsup><mrow><mo>(</mo><mn>0</mn><mo>,</mo><mn>0</mn><mo>)</mo></mrow></mrow></mfrac><mo>)</mo></mrow></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000931212420000024.GIF" wi="1438" he="247" /></maths>6)预编码图像:将求得的投影仪伽马预编码到正弦条纹图像中,经投影仪投射,即可得到保持良好正弦性的正弦光栅条纹;预编码后的正弦条纹图像表示为I<sup>Y</sup>=[a+b cos(2πfx+δ)]<sup>1/γ</sup>   (8)式中:a为正弦函数的基本偏移量,b为正弦函数的幅值,f为正弦函数的数字频率,x为沿x轴的变化量,δ为相移变化量。
地址 300130 天津市红桥区丁字沽光荣道8号河北工业大学东院330#