发明名称 成像光电系统光谱响应非均匀性测量方法
摘要 本发明公开了一种成像光电系统光谱响应非均匀性测量方法,属于光学计量测试领域。该方法利用可调谐激光光源、激光稳功率控制器结合积分球形成比较稳定的均匀的激光束,均匀的激光束经过准直透镜准直后由均匀性补偿器对准直的激光束进行均匀性修正,再由带有二维扫描机构的类针孔辐射计对修正后的均匀激光束进行均匀性测量,然后根据被测成像光电系统的像元数通过算法对激光束的均匀性数组进行匀滑处理,得到被测成像光电系统每个像元对应的修正系数,从而实现成像光电系统不同波长处的光谱响应非均匀性测量。该方法解决了目前成像光电系统光谱响应非均匀性的测量难题,具有测量准确度高,应用前景广的特点。
申请公布号 CN103105286B 申请公布日期 2015.09.23
申请号 CN201310026239.1 申请日期 2013.01.24
申请人 中国兵器工业第二0五研究所 发明人 范纪红;俞兵;侯西旗;袁良;李琪;卢飞;秦艳
分类号 G01M11/02(2006.01)I 主分类号 G01M11/02(2006.01)I
代理机构 代理人
主权项 一种成像光电系统光谱响应非均匀性测量方法,其特征在于:该方法包括以下步骤:第一步:搭建成像光电系统光谱响应非均匀性测量系统将可调谐激光器、激光稳功率控制器、带有监视探测器的积分球、准直透镜、均匀性补偿器、安装在二维扫描机构上的类针孔辐射计依次放置在光学平台上,其中,积分球的出口处于准直透镜的物方焦面上,二维扫描机构受计算机控制,类针孔辐射计的输出送入所述计算机,可调谐激光器输出的激光由激光稳功率控制器的中心入射,由激光稳功率控制器输出的功率稳定的激光垂直入射在所述积分球入口中心,积分球上的监视探测器将探测信号反馈给激光稳功率控制器,由积分球出射的均匀激光由准直透镜形成准直激光束,准直激光束经均匀性补偿器补偿后照射到类针孔辐射计的输入端;第二步:计算机控制测量2.1将可调谐激光器调整到波长λ<sub>1</sub>,计算机首先控制类针孔辐射计测量波长λ<sub>1</sub>处的背景信号,然后控制二维扫描机构带动类针孔辐射计根据设置的扫描间隔对类针孔辐射计输入端的准直激光束截面进行二维扫描,计算机将类针孔辐射计输出的激光信号扣除背景信号得到波长为λ<sub>1</sub>且对应二维扫描坐标为(x<sub>p</sub>,y<sub>q</sub>)的实测信号V(x<sub>p</sub>,y<sub>q</sub>,λ<sub>1</sub>),x<sub>p</sub>表示行方向x的扫描坐标点且p=0、1、2、…、P‑1,P为行方向扫描距离与扫描间距的比值,y<sub>q</sub>表示列方向y的扫描坐标点且q=0、1、2、…、Q‑1,Q为列方向扫描距离与扫描间距的比值,然后根据公式(1)计算激光波长为λ<sub>1</sub>时准直激光束截面的均匀性U(x<sub>p</sub>,y<sub>q</sub>,λ<sub>1</sub>)并保存该组数据:<img file="FSB0000139515870000011.GIF" wi="1237" he="135" />式中,V<sub>最大</sub>(x<sub>p</sub>,y<sub>q</sub>,λ<sub>1</sub>)表示类针孔辐射计在波长λ<sub>1</sub>处所有实测信号值中的最大值;2.2计算机根据被测成像光电系统即CCD的像元数对所述均匀性U(x<sub>p</sub>,y<sub>q</sub>,λ<sub>1</sub>)数组进行匀滑处理,具体操作步骤如下:2.2.1根据公式(2)求出在p=0时q=0和q=1两个扫描点对应到CCD第一列两个像元点之间且含该两个像元点的n个像元的修正系数<img file="FSB0000139515870000026.GIF" wi="154" he="81" /><maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>C</mi><mrow><mn>0</mn><mo>,</mo><mi>k</mi><mo>,</mo><msub><mi>&lambda;</mi><mn>1</mn></msub></mrow></msub><mo>=</mo><mi>U</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>0</mn></msub><mo>,</mo><msub><mi>y</mi><mn>0</mn></msub><mo>,</mo><msub><mi>&lambda;</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>+</mo><mfrac><mrow><mi>U</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>0</mn></msub><mo>,</mo><msub><mi>y</mi><mn>1</mn></msub><mo>,</mo><msub><mi>&lambda;</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>-</mo><mi>U</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>0</mn></msub><mo>,</mo><msub><mi>y</mi><mn>0</mn></msub><mo>,</mo><msub><mi>&lambda;</mi><mn>1</mn></msub><mo>)</mo></mrow></mrow><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></mfrac><mo>&times;</mo><mi>k</mi></mrow>]]></math><img file="FSB0000139515870000021.GIF" wi="1118" he="112" /></maths>(k=0,1,…,n‑1且k<J)    (2)式中,U(x<sub>0</sub>,y<sub>0</sub>,λ<sub>1</sub>)为p=0,q=0扫描点的均匀性且(x<sub>0</sub>,y<sub>0</sub>)为该扫描点的坐标;U(x<sub>0</sub>,y<sub>1</sub>,λ<sub>1</sub>)为p=0,q=1扫描点的均匀性且(x<sub>0</sub>,y<sub>1</sub>)为该扫描点的坐标,J为CCD的像元行数;2.2.2根据公式(3)求出在p=1时q=0和q=1两个扫描点对应到CCD第m列两个像元点之间且含该两个像元点的n个像元的修正系数<img file="FSB0000139515870000027.GIF" wi="177" he="78" />其中m为p=0和p=1两个扫描列对应到CCD两列之间同一行中的像元数:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>C</mi><mrow><mi>m</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>k</mi><mn>1</mn><msub><mi>&lambda;</mi><mn>1</mn></msub></mrow></msub><mo>=</mo><mi>U</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>1</mn></msub><mo>,</mo><msub><mi>y</mi><mn>0</mn></msub><mo>,</mo><msub><mi>&lambda;</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>+</mo><mfrac><mrow><mi>U</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>1</mn></msub><mo>,</mo><msub><mi>y</mi><mn>1</mn></msub><mo>,</mo><msub><mi>&lambda;</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>-</mo><mi>U</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>1</mn></msub><mo>,</mo><msub><mi>y</mi><mn>0</mn></msub><mo>,</mo><msub><mi>&lambda;</mi><mn>1</mn></msub><mo>)</mo></mrow></mrow><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></mfrac><mo>&times;</mo><mi>k</mi></mrow>]]></math><img file="FSB0000139515870000022.GIF" wi="1150" he="113" /></maths>(k=0,1,…,n‑1且k<J)  (3)式中,U(x<sub>1</sub>,y<sub>0</sub>,λ<sub>1</sub>)为p=1,q=0扫描点的均匀性且(x<sub>1</sub>y<sub>0</sub>)为该扫描点的坐标;U(x<sub>1</sub>,y<sub>1</sub>,λ<sub>1</sub>)为p=1,q=1扫描点的均匀性且(x<sub>1</sub>,y<sub>1</sub>)为该扫描点的坐标;2.2.3根据公式(4)求出p=0和p=1时q=0和q=1这四个相邻扫描点对应到CCD四个点所含区域中第一行至第n行所有像元的修正系数<img file="FSB0000139515870000023.GIF" wi="130" he="62" /><img file="FSB0000139515870000024.GIF" wi="740" he="125" />(l=0,1,…,m‑1;k=0,1,…,n‑1且l<I)   (4)式中,I为CCD的像元列数;2.2.4按照2.2.1到2.2.3步骤中的公式逐行求出P×Q个扫描点中每四个相邻扫描点对应到CCD中m×n个像元的修正系数;至此,获得CCD中每个像元对应的修正系数<img file="FSB0000139515870000025.GIF" wi="138" he="58" />i代表CCD像元的行坐标且i=0、1、…、I,j代表CCD像元的列坐标且j=0、1、…、J;2.3将被测成像光电系统安装在光学平台上并使其处于准直激光光路中,计算机控制二维扫描机构带动类针孔辐射计移出光路,用被测成像光电系统分别测量波长λ<sub>1</sub>处的背景和激光并将两种测量值依次送入计算机,计算机将被测成像光电系统每个像元的激光信号值扣除背景信号值得到波长λ<sub>1</sub>处每个像元的实测信号<img file="FSB0000139515870000031.GIF" wi="128" he="54" />计算机根据公式(5)计算每个像元对应的修正值<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><msup><mi>V</mi><mo>&prime;</mo></msup><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><msub><mi>&lambda;</mi><mn>1</mn></msub></mrow></msub><mo>:</mo></mrow>]]></math><img file="FSB0000139515870000032.GIF" wi="142" he="56" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msubsup><mi>V</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><msub><mi>&lambda;</mi><mn>1</mn></msub></mrow><mo>&prime;</mo></msubsup><mo>=</mo><mfrac><msub><mi>V</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><msub><mi>&lambda;</mi><mn>1</mn></msub></mrow></msub><msub><mi>C</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><msub><mi>&lambda;</mi><mn>1</mn></msub></mrow></msub></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FSB0000139515870000033.GIF" wi="1115" he="152" /></maths>2.4计算机根据公式(6)计算被测成像光电系统在第一波长λ<sub>1</sub>处光谱响应非均匀性NU(λ<sub>1</sub>):<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>NU</mi><mrow><mo>(</mo><msub><mi>&lambda;</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>=</mo><msqrt><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mi>I</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>0</mn></mrow><mi>J</mi></munderover><msup><mrow><mo>(</mo><msub><msup><mi>V</mi><mo>&prime;</mo></msup><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><msub><mi>&lambda;</mi><mn>1</mn></msub></mrow></msub><mo>-</mo><mover><msub><msup><mi>V</mi><mo>&prime;</mo></msup><msub><mi>&lambda;</mi><mn>1</mn></msub></msub><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mn>2</mn></msup></mrow><mrow><msup><mover><msubsup><mi>V</mi><msub><mi>&lambda;</mi><mn>1</mn></msub><mo>&prime;</mo></msubsup><mo>&OverBar;</mo></mover><mn>2</mn></msup><mrow><mo>(</mo><mi>M</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow></mfrac></msqrt><mo>&times;</mo><mn>100</mn><mo>%</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FSB0000139515870000034.GIF" wi="1349" he="243" /></maths>式中,<img file="FSB0000139515870000035.GIF" wi="85" he="71" />表示在波长λ<sub>1</sub>处被测成像光电系统各像元修正值的平均值,M表示被测成像光电系统像元总数;2.5将可调谐激光器调整到波长λ<sub>d</sub>且d=2、3…、D,在每个波长点下重复第2.1步至第2.4步,最终保存非均匀性NU(λ<sub>d</sub>)数据且d=1、2、…、D。
地址 710065 陕西省西安市雁塔区电子三路西段九号