发明名称 四通道大气偏振信息检测传感器信号处理与补偿方法
摘要 四通道大气偏振信息检测传感器信号处理与补偿方法,其特征是采用逐级信号处理的方法,对于来自四通道大气偏振信息检测传感器的实时动态检测数据进行逐级处理,在逐级处理过程中,通过硬补偿和软补偿两种分级的补偿方法,对传感器逐级信号的误差和噪声进行信号补偿与匹配,使得由于传感器各通道逐级物理器件的信号响应与传输特性不同所造成的传感器各通道信号的漂移误差分量和幅值误差,以及由于传感器的物理噪声和环境光噪声,所引起的随机噪声,得到有效的一致性补偿和抑制,最终得到有效的传感器输出,通过对获得的数字量电信号进行大气偏振模式信息处理,得到准确、可靠的大气偏振模式信息。
申请公布号 CN101441171A 申请公布日期 2009.05.27
申请号 CN200810246233.4 申请日期 2008.12.31
申请人 合肥工业大学 发明人 高 隽;范之国;胡 翊;潘登凯;范宁生
分类号 G01N21/21(2006.01)I 主分类号 G01N21/21(2006.01)I
代理机构 安徽省合肥新安专利代理有限责任公司 代理人 何梅生
主权项 1、四通道大气偏振信息检测传感器信号处理与补偿方法,其特征是按如下过程进行:采用逐级信号处理的方法,首先,利用四通道大气偏振信息检测传感器,获得四个通道的偏振光信号分别是P<sub>11</sub>、P<sub>12</sub>、P<sub>21</sub>和P<sub>22</sub>,有:P<sub>11</sub>=n<sub>11</sub>I(1+d cos 2(φ-φ<sub>11</sub>))          (1)P<sub>12</sub>=n<sub>12</sub>I(1+d cos 2(φ-φ<sub>12</sub>))          (2)P<sub>21</sub>=n<sub>21</sub>I(1+d cos 2(φ-φ<sub>21</sub>))          (3)P<sub>22</sub>=n<sub>22</sub>I(1+d cos 2(φ-φ<sub>22</sub>))          (4)式(1)、式(2)、式(3)和式(4)中的n<sub>11</sub>、n<sub>12</sub>、n<sub>21</sub>和n<sub>22</sub>分别为所述传感器中四个通道的偏振感知器对特定波段光的传输响应系数,I为归一化光强系数,d为大气偏振模式中偏振度信息,φ为大气偏振模式中偏振化方向信息,即检测时刻的太阳方位角,φ<sub>11</sub>、φ<sub>12</sub>、φ<sub>21</sub>和φ<sub>22</sub>分别为四个通道的偏振感知器的起始偏振感知方向,且有:φ<sub>12</sub>-φ<sub>11</sub>=π/2、φ<sub>22</sub>-φ<sub>21</sub>=π/2、0&lt;φ<sub>21</sub>-φ<sub>11</sub>&lt;π/2;所述四个通道的偏振感知器输出的偏振光信号P<sub>11</sub>、P<sub>12</sub>、P<sub>21</sub>和P<sub>22</sub>分别经光电转换器后输出的模拟量电信号分别是I<sub>11</sub>、I<sub>12</sub>、I<sub>21</sub>和I<sub>22</sub>,有:I<sub>11</sub>=m<sub>11</sub>P<sub>11</sub>=l<sub>11</sub>I(1+d cos 2(φ-φ<sub>11</sub>))     (5)I<sub>12</sub>=m<sub>12</sub>P<sub>12</sub>=l<sub>12</sub>I(1+d cos 2(φ-φ<sub>12</sub>))     (6)I<sub>21</sub>=m<sub>21</sub>P<sub>21</sub>=l<sub>21</sub>I(1+d cos 2(φ-φ<sub>21</sub>))     (7)I<sub>22</sub>=m<sub>22</sub>P<sub>22</sub>=l<sub>22</sub>I(1+d cos 2(φ-φ<sub>22</sub>))     (8)式(5)、式(6)、式(7)和式(8)中的m<sub>11</sub>、m<sub>12</sub>、m<sub>21</sub>和m<sub>22</sub>分别为四个通道的光电转换器的传输响应系数,l<sub>11</sub>=m<sub>11</sub>×n<sub>11</sub>、l<sub>12</sub>=m<sub>12</sub>×n<sub>12</sub>、l<sub>21</sub>=m<sub>21</sub>×n<sub>21</sub>和l<sub>22</sub>=m<sub>22</sub>×n<sub>22</sub>为偏振感知器与光电转换器的合传输系数;所述四个通道的光电转换器输出的模拟量电信号I<sub>11</sub>、I<sub>12</sub>、I<sub>21</sub>和I<sub>22</sub>,以两个通道为一组,分别经对数运算单元,输出的模拟量电信号分别是G<sub>1</sub>和G<sub>2</sub>,有:<maths num="0001"><![CDATA[<math><mrow><msub><mi>G</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>&phi;</mi><mo>)</mo></mrow><mo>=</mo><mi>log</mi><mrow><mo>(</mo><mfrac><msub><mi>l</mi><mn>11</mn></msub><msub><mi>l</mi><mn>12</mn></msub></mfrac><mo>)</mo></mrow><mo>+</mo><mi>log</mi><mrow><mo>(</mo><mfrac><mrow><mn>1</mn><mo>+</mo><mi>d</mi><mi>cos</mi><mrow><mo>(</mo><mn>2</mn><mi>&phi;</mi><mo>)</mo></mrow></mrow><mrow><mn>1</mn><mo>-</mo><mi>d</mi><mi>cos</mi><mrow><mo>(</mo><mn>2</mn><mi>&phi;</mi><mo>)</mo></mrow></mrow></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><msub><mi>G</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>&phi;</mi><mo>)</mo></mrow><mo>=</mo><mi>log</mi><mrow><mo>(</mo><mfrac><msub><mi>l</mi><mn>21</mn></msub><msub><mi>l</mi><mn>22</mn></msub></mfrac><mo>)</mo></mrow><mo>+</mo><mi>log</mi><mrow><mo>(</mo><mfrac><mrow><mn>1</mn><mo>+</mo><mi>d</mi><mi>cos</mi><mrow><mo>(</mo><mn>2</mn><mi>&phi;</mi><mo>-</mo><msub><mrow><mn>2</mn><mi>&phi;</mi></mrow><mn>0</mn></msub><mo>)</mo></mrow></mrow><mrow><mn>1</mn><mo>-</mo><mi>d</mi><mi>cos</mi><mrow><mo>(</mo><mn>2</mn><mi>&phi;</mi><mo>-</mo><msub><mrow><mn>2</mn><mi>&phi;</mi></mrow><mn>0</mn></msub><mo>)</mo></mrow></mrow></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow></math>]]></maths>式(9)和式(10)中,0&lt;φ<sub>0</sub>=φ<sub>21</sub>-φ<sub>11</sub>&lt;π/2,<img file="A200810246233C00033.GIF" wi="169" he="136" />和<img file="A200810246233C00034.GIF" wi="174" he="136" />即是由于传感器中四个通道的偏振光感知器与光电转换器的信号响应特性不一致所造成的加性漂移误差分量;所述两路对数运算单元输出的模拟量电信号G<sub>1</sub>和G<sub>2</sub>,分别经两路运算与放大单元,输出的模拟量电信号为<img file="A200810246233C00035.GIF" wi="44" he="48" />和<img file="A200810246233C00036.GIF" wi="67" he="49" />有:<maths num="0003"><![CDATA[<math><mrow><msub><mi>F</mi><mn>1</mn></msub><mo>&prime;</mo><mrow><mo>(</mo><mi>&phi;</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>k</mi><mn>1</mn></msub><mi>log</mi><mrow><mo>(</mo><mfrac><msub><mi>l</mi><mn>11</mn></msub><msub><mi>l</mi><mn>12</mn></msub></mfrac><mo>)</mo></mrow><mo>+</mo><msub><mi>k</mi><mn>1</mn></msub><mi>log</mi><mrow><mo>(</mo><mfrac><mrow><mn>1</mn><mo>+</mo><mi>d</mi><mi>cos</mi><mrow><mo>(</mo><mn>2</mn><mi>&phi;</mi><mo>)</mo></mrow></mrow><mrow><mn>1</mn><mo>-</mo><mi>d</mi><mi>cos</mi><mrow><mo>(</mo><mn>2</mn><mi>&phi;</mi><mo>)</mo></mrow></mrow></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><msubsup><mi>F</mi><mn>2</mn><mo>&prime;</mo></msubsup><mrow><mo>(</mo><mi>&phi;</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>k</mi><mn>2</mn></msub><mi>log</mi><mrow><mo>(</mo><mfrac><msub><mi>l</mi><mn>21</mn></msub><msub><mi>l</mi><mn>22</mn></msub></mfrac><mo>)</mo></mrow><mo>+</mo><msub><mi>k</mi><mn>2</mn></msub><mi>log</mi><mrow><mo>(</mo><mfrac><mrow><mn>1</mn><mo>+</mo><mi>d</mi><mi>cos</mi><mrow><mo>(</mo><mn>2</mn><mi>&phi;</mi><mo>-</mo><msub><mrow><mn>2</mn><mi>&phi;</mi></mrow><mn>0</mn></msub><mo>)</mo></mrow></mrow><mrow><mn>1</mn><mo>-</mo><mi>d</mi><mi>cos</mi><mrow><mo>(</mo><mn>2</mn><mi>&phi;</mi><mo>-</mo><msub><mrow><mn>2</mn><mi>&phi;</mi></mrow><mn>0</mn></msub><mo>)</mo></mrow></mrow></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>12</mn><mo>)</mo></mrow></mrow></math>]]></maths>式(11)和式(12)中,k1和k2分别为两路运算与放大单元的传输响应系数,k<sub>1</sub>和k<sub>2</sub>的不同造成了运算与放大单元后输出的模拟量电信号存在乘性幅值误差,<img file="A200810246233C00039.GIF" wi="213" he="138" />和<img file="A200810246233C000310.GIF" wi="219" he="140" />就是运算与放大单元后输出模拟量电信号中的被放大了的加性漂移误差分量;针对加性漂移误差分量和乘性幅值误差,按以下方式进行硬补偿:所述硬补偿是在运算与放大过程中,调节两路运算与放大单元中的放大反馈参数,用以补偿两路运算与放大单元中由于k<sub>1</sub>和k<sub>2</sub>的不同所造成的乘性幅值误差,使得经过两路运算与放大单元后的模拟量电信号的幅值相等;在运算与放大过程中,调节两路运算与放大单元中的基准电平信号,用以补偿<img file="A200810246233C000311.GIF" wi="213" he="139" />和<img file="A200810246233C000312.GIF" wi="219" he="138" />所造成的加性漂移误差分量,使得经过两路运算与放大单元后的模拟量电信号最大最小值调谐至零对称,经过硬补偿后的两路运算与放大单元输出的模拟量电信号分别为F<sub>1</sub>和F<sub>2</sub>,有:<maths num="0005"><![CDATA[<math><mrow><msub><mi>F</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>&phi;</mi><mo>)</mo></mrow><mo>=</mo><mi>log</mi><mrow><mo>(</mo><mfrac><mrow><mn>1</mn><mo>+</mo><mi>d</mi><mi>cos</mi><mrow><mo>(</mo><mn>2</mn><mi>&phi;</mi><mo>)</mo></mrow></mrow><mrow><mn>1</mn><mo>-</mo><mi>d</mi><mi>cos</mi><mrow><mo>(</mo><mn>2</mn><mi>&phi;</mi><mo>)</mo></mrow></mrow></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>13</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0006"><![CDATA[<math><mrow><msub><mi>F</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>&phi;</mi><mo>)</mo></mrow><mo>=</mo><mi>log</mi><mrow><mo>(</mo><mfrac><mrow><mn>1</mn><mo>+</mo><mi>d</mi><mi>cos</mi><mrow><mo>(</mo><mn>2</mn><mi>&phi;</mi><mo>-</mo><msub><mrow><mn>2</mn><mi>&phi;</mi></mrow><mn>0</mn></msub><mo>)</mo></mrow></mrow><mrow><mn>1</mn><mo>-</mo><mi>d</mi><mi>cos</mi><mrow><mo>(</mo><mn>2</mn><mi>&phi;</mi><mo>-</mo><msub><mrow><mn>2</mn><mi>&phi;</mi></mrow><mn>0</mn></msub><mo>)</mo></mrow></mrow></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>14</mn><mo>)</mo></mrow></mrow></math>]]></maths>硬补偿后的模拟量电信号F<sub>1</sub>和F<sub>2</sub>,经过模数转换单元,转换成数字量电信号;对于所述数字量电信号,若误差超过规定的阈值,再按以下方法进行软补偿:所述软补偿是当检测数据中的误差超过规定的阈值时,则启动软补偿,所述软补偿为自适应实时校准过程,首先,传感器以自体中心为轴,匀速旋转不少于360度,在旋转过程中等间隔采样,获取两路模数转换单元的采样点数据,并对所有采样点数据进行统计与计算,分别计算出两路模数转换单元之间的信号总体误差影响因子,总体误差影响因子包括均值误差比和峰值误差比,在实际检测过程中利用总体误差影响因子对两路模数转换单元输出的数字量电信号进行自适应的实时补偿与调整,用以抵消由于传感器各通道逐级物理器件的信号响应与传输特性不同所造成的传感器各通道信号的漂移误差分量和幅值误差,以及由于传感器的物理噪声和环境光噪声所引起的随机噪声,提高模数转换单元所获得的数字量电信号的检测精度与稳定性;最后,获得的数字量电信号进行大气偏振模式信息处理,得到大气偏振模式的偏振度信息d和偏振化方向信息φ,即:<maths num="0007"><![CDATA[<math><mrow><mi>&phi;</mi><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mi>arctan</mi><mrow><mo>(</mo><mfrac><mrow><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mn>2</mn><msub><mover><mi>F</mi><mo>&OverBar;</mo></mover><mn>2</mn></msub><mrow><mo>(</mo><mi>&phi;</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mn>2</mn><msub><mover><mi>F</mi><mo>&OverBar;</mo></mover><mn>1</mn></msub><mrow><mo>(</mo><mi>&phi;</mi><mo>)</mo></mrow><mo>)</mo></mrow><mi>cos</mi><mrow><mo>(</mo><msub><mrow><mn>2</mn><mi>&phi;</mi></mrow><mn>0</mn></msub><mo>)</mo></mrow></mrow><mrow><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mn>2</mn><msub><mover><mi>F</mi><mo>&OverBar;</mo></mover><mn>1</mn></msub><mrow><mo>(</mo><mi>&phi;</mi><mo>)</mo></mrow><mo>)</mo></mrow><mi>sin</mi><mrow><mo>(</mo><msub><mrow><mn>2</mn><mi>&phi;</mi></mrow><mn>0</mn></msub><mo>)</mo></mrow></mrow></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>15</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0008"><![CDATA[<math><mrow><mi>d</mi><mo>=</mo><msqrt><msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mn>2</mn><msub><mover><mi>F</mi><mo>&OverBar;</mo></mover><mn>1</mn></msub><mrow><mo>(</mo><mi>&phi;</mi><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mn>2</mn><msub><mover><mi>F</mi><mo>&OverBar;</mo></mover><mn>2</mn></msub><mrow><mo>(</mo><mi>&phi;</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mn>2</mn><msub><mover><mi>F</mi><mo>&OverBar;</mo></mover><mn>1</mn></msub><mrow><mo>(</mo><mi>&phi;</mi><mo>)</mo></mrow><mo>)</mo></mrow><mi>cos</mi><mrow><mo>(</mo><msub><mrow><mn>2</mn><mi>&phi;</mi></mrow><mn>0</mn></msub><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup><mo>/</mo><msup><mi>sin</mi><mn>2</mn></msup><mrow><mo>(</mo><msub><mrow><mn>2</mn><mi>&phi;</mi></mrow><mn>0</mn></msub><mo>)</mo></mrow></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>16</mn><mo>)</mo></mrow></mrow></math>]]></maths>式(15)和式(16)中的F<sub>1</sub>(φ)和F<sub>2</sub>(φ)是式(13)和式(14)中的F<sub>1</sub>(φ)和F<sub>2</sub>(φ)由S型函数<maths num="0009"><![CDATA[<math><mrow><mover><mi>F</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>&phi;</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><msup><mn>10</mn><mrow><mi>F</mi><mrow><mo>(</mo><mi>&phi;</mi><mo>)</mo></mrow></mrow></msup><mo>+</mo><mn>1</mn></mrow></mfrac></mrow></math>]]></maths>进行去对数处理后而得到。
地址 230009安徽省合肥市屯溪路193号