发明名称 一种NMOS线阵图像传感器系统固定误差动态标定与校正方法
摘要 一种NMOS线阵图像传感器系统固定误差动态标定与校正方法属于线阵图像传感器系统的误差校正方法领域,该方法将导致NMOS线阵图像传感器及其外部读出电路的各种误差分为与积分时间成正比的变量部分和一个常量部分,并在此基础上建立系统固定误差的数理模型,设计图像传感器的误差实时测量实验,根据实测结果最终给出基于该模型的固定误差的动态标定与校正方法。该方法摒弃了传统的针对各种误差源,分别从系统硬件设计或图像传感器工作环境条件上采取措施,以求减小测量结果中固定误差的长周期、高成本处理方式。
申请公布号 CN104864960A 申请公布日期 2015.08.26
申请号 CN201510217126.9 申请日期 2015.04.30
申请人 中国科学院长春光学精密机械与物理研究所 发明人 张佩杰;宋克非
分类号 G01J3/28(2006.01)I 主分类号 G01J3/28(2006.01)I
代理机构 长春菁华专利商标代理事务所 22210 代理人 田春梅
主权项 一种NMOS线阵图像传感器系统固定误差动态标定与校正方法,其特征在于,该方法包括如下步骤:步骤一、对NMOS线阵图像传感器及其外部读出电路的系统固定误差进行整体数学建模,具体方法为:根据系统固定误差与积分时间的关系类型,将NMOS线阵图像传感器及其外部读出电路的系统固定误差的数学模型表示为:E=E<sub>R</sub>+T×E<sub>D</sub>   (1)式(1)中,E表示NMOS线阵图像传感器及其外部读出电路的系统固定误差的总量,T为积分时间(单位为秒),E<sub>D</sub>表示随积分时间变化的固定误差在1秒的单位时间内的累积量,E<sub>R</sub>表示不随积分时间变化的误差的总量;即图像传感器的m个像元中,任意一个像元的固定误差可以表示为:E<sub>k</sub>=E<sub>Rk</sub>+T×E<sub>Dk</sub>   (2)式(2)中,E<sub>k</sub>是第k个像元的固定误差总量;E<sub>Rk</sub>是第k个像元不随积分时间变化的误差的总量;E<sub>Dk</sub>表示第k个像元随积分时间变化的固定误差在单位时间内的累积量;步骤二、在无光照条件下进行N(N≥2,N为自然数)次误差测量,获得N组误差测量值,其具体包括如下子步骤:步骤2.1:在无光照条件下,使NMOS线阵图像传感器进行第一次测量,积分时间为T<sub>1</sub>,得到由m个像元各自在第一次测量中的输出信号值组成的第一组误差测量值,A<sub>s1</sub>={E<sub>11</sub>,E<sub>21</sub>,……E<sub>k1</sub>,……E<sub>m1</sub>}   (3)其中,E<sub>k1</sub>表示第k个像元在第1次误差测量中的输出信号值;步骤2.2:使用与步骤2.1相同的方法,在无光照条件下,使NMOS线阵图像传感器再进行N‑1次测量,并且使N‑1次测量中每次测量的积分时间按照如下要求进行取值:T<sub>1</sub>&lt;T<sub>2</sub>&lt;…&lt;T<sub>i</sub>&lt;…&lt;T<sub>N</sub>   (4)式(4)中,i(i≤N,i取自然数)为测量的次数序号;T<sub>i</sub>为第i次测量中积分时间的取值;T<sub>1</sub>为第1次测量中积分时间的取值;T<sub>N</sub>为第N次测量中积分时间的取值;N‑1次测量后,得到N‑1组误差测量值A<sub>s2</sub>,A<sub>s3</sub>,……A<sub>si</sub>,……A<sub>sN</sub>,其中,A<sub>si</sub>是由m个像元各自在第i次测量中输出信号值E<sub>1i</sub>,E<sub>2i</sub>,……E<sub>ki</sub>,……,E<sub>mi</sub>组成的第i组误差测量值;步骤二中,共获得N组误差测量值A<sub>s1</sub>,A<sub>s2</sub>,A<sub>s3</sub>,……A<sub>si</sub>,……A<sub>sN</sub>;步骤三、利用步骤二所述N组误差测量值,对步骤一所述像元固定误差数学模型中的未知参数进行求解,m个像元中每个像元都进行一次求解,得到对应的一组误差参数E<sub>Dk</sub>、E<sub>Rk</sub>;其具体包括如下子步骤:步骤3.1:在步骤二的N组误差测量值A<sub>s1</sub>,A<sub>s2</sub>,……A<sub>si</sub>,……A<sub>sN</sub>中,将第一个像元在N次误差测量中的误差输出值E<sub>11</sub>、E<sub>12</sub>、……、E<sub>1N</sub>取出,将N次误差测量中使用的N个积分时间T<sub>1</sub>、T<sub>2</sub>、……、T<sub>N</sub>取出,根据步骤1所述的误差模型,可以得到m个方程:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfenced open='' close=''><mtable><mtr><mtd><msub><mi>E</mi><mn>11</mn></msub><mo>=</mo><msub><mover><mi>E</mi><mo>^</mo></mover><mrow><mi>R</mi><mn>1</mn></mrow></msub><mo>+</mo><msub><mi>T</mi><mn>1</mn></msub><mo>&times;</mo><msub><mover><mi>E</mi><mo>^</mo></mover><mrow><mi>D</mi><mn>1</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>E</mi><mn>12</mn></msub><mo>=</mo><msub><mover><mi>E</mi><mo>^</mo></mover><mrow><mi>R</mi><mn>1</mn></mrow></msub><mo>+</mo><msub><mi>T</mi><mn>2</mn></msub><mo>&times;</mo><msub><mover><mi>E</mi><mo>^</mo></mover><mrow><mi>D</mi><mn>1</mn></mrow></msub></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>E</mi><mrow><mn>1</mn><mi>N</mi></mrow></msub><mo>=</mo><msub><mover><mi>E</mi><mo>^</mo></mover><mrow><mi>R</mi><mn>1</mn></mrow></msub><mo>+</mo><msub><mi>T</mi><mi>N</mi></msub><mo>&times;</mo><msub><mover><mi>E</mi><mo>^</mo></mover><mrow><mi>D</mi><mn>1</mn></mrow></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000709273370000021.GIF" wi="1024" he="409" /></maths>式(5)中,<img file="FDA0000709273370000022.GIF" wi="92" he="84" />是第一个像元输出信号值中不随积分时间变化的误差E<sub>R1</sub>的估计值;<img file="FDA0000709273370000023.GIF" wi="82" he="83" />是第一个像元输出信号值中随积分时间变化的固定误差在单位时间内累积量E<sub>D1</sub>的估计值;步骤3.2:应用最小二乘法求解式(5)中<img file="FDA0000709273370000031.GIF" wi="80" he="84" />和<img file="FDA0000709273370000032.GIF" wi="105" he="84" /><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mfenced open='' close=''><mtable><mtr><mtd><msub><mover><mi>E</mi><mo>^</mo></mover><mrow><mi>R</mi><mn>1</mn></mrow></msub><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msup><mrow><mo>(</mo><msub><mi>T</mi><mi>i</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>E</mi><mi>i</mi></msub><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>T</mi><mi>i</mi></msub><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mrow><mo>(</mo><msub><mi>T</mi><mi>i</mi></msub><msub><mi>E</mi><mi>i</mi></msub><mo>)</mo></mrow></mrow><mrow><mi>N</mi><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msup><mrow><mo>(</mo><msub><mi>T</mi><mi>i</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><msup><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>T</mi><mi>i</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></mfrac></mtd></mtr><mtr><mtd><msub><mover><mi>E</mi><mo>^</mo></mover><mrow><mi>D</mi><mn>1</mn></mrow></msub><mo>=</mo><mfrac><mrow><mi>N</mi><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mrow><mo>(</mo><msub><mi>T</mi><mi>i</mi></msub><msub><mi>E</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>T</mi><mi>i</mi></msub><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>E</mi><mi>i</mi></msub></mrow><mrow><mi>N</mi><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msup><mrow><mo>(</mo><msub><mi>T</mi><mi>i</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><msup><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>T</mi><mi>i</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></mfrac></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000709273370000033.GIF" wi="1226" he="628" /></maths>式(6)中i取自然数;将<img file="FDA0000709273370000034.GIF" wi="80" he="86" />和<img file="FDA0000709273370000035.GIF" wi="76" he="86" />作为第1个像元E<sub>D1</sub>和E<sub>R1</sub>的数值大小,即:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>E</mi><mrow><mi>D</mi><mn>1</mn></mrow></msub><mo>=</mo><msub><mover><mi>E</mi><mo>^</mo></mover><mrow><mi>D</mi><mn>1</mn></mrow></msub><mo>,</mo><msub><mi>E</mi><mrow><mi>R</mi><mn>1</mn></mrow></msub><mo>=</mo><msub><mover><mi>E</mi><mo>^</mo></mover><mrow><mi>R</mi><mn>1</mn></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000709273370000036.GIF" wi="896" he="92" /></maths>步骤3.3:对于第2、3、…、m个像元中的每个像元,重复步骤3.1、3.2,求出第2、3、…、m个像元中每个像元所对应的E<sub>Dk</sub>、E<sub>Rk</sub>,即分别求得E<sub>D2</sub>和E<sub>R2</sub>;E<sub>D3</sub>和E<sub>R3</sub>……E<sub>Dm</sub>和E<sub>Rm</sub>;步骤三中,对每个像元都会求解出一组误差参数E<sub>Dk</sub>、E<sub>Rk</sub>,m个像元总共会有m组误差参数,即E<sub>D1</sub>和E<sub>R1</sub>;E<sub>D2</sub>和E<sub>R2</sub>;E<sub>D3</sub>和E<sub>R3</sub>……E<sub>Dm</sub>和E<sub>Rm</sub>;步骤四、在用NMOS线阵图像传感器接对外部入射光信号进行测量时,利用步骤二、三得到误差模型与参数,并对测量结果进行校正,其具体包括如下子步骤:步骤4.1:在用NMOS线阵图像传感器每次开始工作时或工作过程中每隔一段时间,遮挡入射到传感器的光信号,按照步骤二、三进行固定误差的实时测量与标定,得到图像传感器系统中m个像元的固定误差参数E<sub>D1</sub>和E<sub>R1</sub>;E<sub>D2</sub>和E<sub>R2</sub>;E<sub>D3</sub>和E<sub>R3</sub>……E<sub>Dm</sub>和E<sub>Rm</sub>;步骤4.2:固定误差测量标定完成后,NMOS线阵图像传感器正常工作,进行光谱测量,并在每次测量读出时,针对NMOS线阵图像传感器的每一个像元,利用该次测量的积分时间和固定误差标定得到的模型与参数对像元的输出进行校正;某次测量的积分时间为t,m个像元的输出分别为y<sub>1</sub>,y<sub>2</sub>,…,y<sub>m</sub>;m个像元中每个像元固定误差的大小计算为:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mfenced open='' close=''><mtable><mtr><mtd><msub><mi>e</mi><mn>1</mn></msub><mo>=</mo><msub><mi>E</mi><mrow><mi>R</mi><mn>1</mn></mrow></msub><mo>+</mo><mi>t</mi><mo>&times;</mo><msub><mi>E</mi><mrow><mi>D</mi><mn>1</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>e</mi><mn>2</mn></msub><mo>=</mo><msub><mi>E</mi><mrow><mi>R</mi><mn>2</mn></mrow></msub><mo>+</mo><mi>t</mi><mo>&times;</mo><msub><mi>E</mi><mrow><mi>D</mi><mn>2</mn></mrow></msub></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>e</mi><mi>m</mi></msub><mo>=</mo><msub><mi>E</mi><mi>Rm</mi></msub><mo>+</mo><mi>t</mi><mo>&times;</mo><msub><mi>E</mi><mi>Dm</mi></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000709273370000041.GIF" wi="935" he="308" /></maths>步骤4.3:利用公式(8)的计算结果,对m个像元的输出分别进行校正:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mfenced open='' close=''><mtable><mtr><mtd><msub><mover><mi>y</mi><mo>^</mo></mover><mn>1</mn></msub><mo>=</mo><msub><mi>y</mi><mn>1</mn></msub><mo>-</mo><msub><mi>e</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><msub><mover><mi>y</mi><mo>^</mo></mover><mn>2</mn></msub><mo>=</mo><msub><mi>y</mi><mn>2</mn></msub><mo>-</mo><msub><mi>e</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><msub><mover><mi>y</mi><mo>^</mo></mover><mi>m</mi></msub><mo>=</mo><msub><mi>y</mi><mi>m</mi></msub><mo>-</mo><msub><mi>e</mi><mi>m</mi></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000709273370000042.GIF" wi="891" he="308" /></maths>当用步骤一至四所述方法分别求得<img file="FDA0000709273370000043.GIF" wi="344" he="77" />时,即完成了对NMOS线阵图像传感器及其外部读出电路的系统固定误差的校正过程。
地址 130033 吉林省长春市东南湖大路3888号