发明名称 TDI成像的像移自配准方法
摘要 TDI成像的像移自配准方法,涉及光电探测成像领域,解决现有像移补偿方法增加成像系统的设计难度、复杂度及成本,并且存在补偿范围小、精度低以及实时性差等问题,本发明基于像移速度得到图像传感器每帧像素的位置偏差,在数字域TDI过程中,根据偏差值进行行列配准和对位累加,从而实现在TDI的同时完成像素级像移补偿的目的。其中像移速度矢的角度偏差通过列向配准来调整;像移速度矢的速度大小偏差通过行向配准来调整。本发明采用设计的像移自配准方法在扫描成像的时间延迟积分过程中实现配准像素的对位累加,依赖于灵活的数字域操作,从而简单有效地实现不同像移方向和像移速度大小情形下的TDI相机扫描清晰成像。
申请公布号 CN103985129B 申请公布日期 2017.01.18
申请号 CN201410224831.7 申请日期 2014.05.26
申请人 中国科学院长春光学精密机械与物理研究所 发明人 陶淑苹;曲宏松;郑晓云;金光
分类号 G06T7/00(2006.01)I;H04N5/232(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 长春菁华专利商标代理事务所 22210 代理人 陶尊新
主权项 TDI成像的像移自配准方法,构建像移自配准算法理论模型,其特征是,该方法由以下步骤实现:步骤一、实时采集当前帧图像P<sub>i+1</sub>(m,n),并计算得到当前帧图像的像移量S<sub>P</sub>(i)对应的像素个数Δn(i),用公式表示为:<maths num="0001"><math><![CDATA[<mrow><mi>&Delta;</mi><mi>n</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>S</mi><mi>P</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></mrow><mi>a</mi></mfrac></mrow>]]></math><img file="FDA0001050559620000011.GIF" wi="253" he="109" /></maths>其中,a为像元尺寸;步骤二、判断步骤一中所述的像素个数Δn(i)是否为整数,如果是,则执行步骤三,如果否,则进行图像插值,用公式表示为:<maths num="0002"><math><![CDATA[<mrow><msubsup><mi>P</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow><mo>&prime;</mo></msubsup><mrow><mo>(</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mtable><mtr><mtd><mrow><mi>&alpha;</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><msubsup><mi>P</mi><mi>i</mi><mo>&prime;</mo></msubsup><mrow><mo>(</mo><mi>m</mi><mo>+</mo><mi>&Delta;</mi><mi>&beta;</mi><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>,</mo><mi>n</mi><mo>)</mo><mo>+</mo><mo>...</mo></mrow></mtd></mtr><mtr><mtd><mrow><mo>...</mo><mo>+</mo><mo>&lsqb;</mo><mn>1</mn><mo>-</mo><mi>&alpha;</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>&rsqb;</mo><msub><mi>P</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow></mrow></mtd></mtr></mtable><mo>,</mo></mrow></mtd><mtd><mrow><mi>m</mi><mo>&Element;</mo><mi>R</mi><mo>&lsqb;</mo><mn>1</mn><mo>,</mo><msub><mi>n</mi><mrow><mi>l</mi><mi>i</mi><mi>n</mi><mi>e</mi><mi>s</mi></mrow></msub><mo>-</mo><mi>&Delta;</mi><mi>&beta;</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>&rsqb;</mo></mrow></mtd></mtr><mtr><mtd><mrow><mtable><mtr><mtd><mrow><mi>&alpha;</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><msub><mi>P</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>m</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>+</mo><mo>...</mo></mrow></mtd></mtr><mtr><mtd><mrow><mo>...</mo><mo>+</mo><mo>&lsqb;</mo><mn>1</mn><mo>-</mo><mi>&alpha;</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>&rsqb;</mo><msub><mi>P</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow></mrow></mtd></mtr></mtable><mo>,</mo></mrow></mtd><mtd><mrow><mi>m</mi><mo>&Element;</mo><mi>R</mi><mo>&lsqb;</mo><msub><mi>n</mi><mrow><mi>l</mi><mi>i</mi><mi>n</mi><mi>e</mi><mi>s</mi></mrow></msub><mo>-</mo><mi>&Delta;</mi><mi>&beta;</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>+</mo><mn>1</mn><mo>,</mo><msub><mi>n</mi><mrow><mi>l</mi><mi>i</mi><mi>n</mi><mi>e</mi><mi>s</mi></mrow></msub><mo>&rsqb;</mo></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0001050559620000012.GIF" wi="1326" he="382" /></maths>式中:α(i)=rem(Δn(i))β(i)=int(Δn(i))Δβ(i)=β(i)‑β(i‑1);步骤三、采用向上追溯l个邻近帧的方法插值计算缺失像素值,其表达式如下:式中,x和k为中间变量,M表示积分级数;<maths num="0003"><math><![CDATA[<mrow><msubsup><mi>P</mi><mrow><mi>i</mi><mo>+</mo><mi>j</mi><mo>+</mo><mn>1</mn></mrow><mo>&prime;</mo></msubsup><mrow><mo>(</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mn>3</mn></munderover><msub><mi>P</mi><mrow><mi>i</mi><mo>+</mo><mi>j</mi><mo>-</mo><mi>k</mi></mrow></msub><mo>&lsqb;</mo><mn>1</mn><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>x</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></munderover><mi>&beta;</mi><mrow><mo>(</mo><mi>i</mi><mo>+</mo><mi>j</mi><mo>-</mo><mi>x</mi><mo>)</mo></mrow><mo>,</mo><mi>n</mi><mo>&rsqb;</mo></mrow><mi>l</mi></mfrac><mo>,</mo></mrow></mtd><mtd><mrow><mi>M</mi><mo>&GreaterEqual;</mo><mi>l</mi></mrow></mtd></mtr><mtr><mtd><mrow><mfrac><mrow><msub><mi>P</mi><mrow><mi>i</mi><mo>+</mo><mi>j</mi><mo>-</mo><mi>k</mi></mrow></msub><mo>&lsqb;</mo><mn>1</mn><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>x</mi><mo>=</mo><mn>2</mn></mrow><mi>M</mi></munderover><mi>&beta;</mi><mrow><mo>(</mo><mi>i</mi><mo>+</mo><mi>j</mi><mo>-</mo><mi>x</mi><mo>)</mo></mrow><mo>,</mo><mi>n</mi><mo>&rsqb;</mo></mrow><mi>M</mi></mfrac><mo>,</mo></mrow></mtd><mtd><mrow><mi>M</mi><mo>&lt;</mo><mi>l</mi></mrow></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0001050559620000013.GIF" wi="1438" he="502" /></maths>步骤四、插值完成后与上一帧像素矩阵叠加,得到时间延迟积分后的第i帧像素矩阵P″<sub>i+1</sub>:<maths num="0004"><math><![CDATA[<mrow><msubsup><mi>P</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msubsup><mrow><mo>(</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msubsup><mi>P</mi><mi>i</mi><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msubsup><mrow><mo>(</mo><mi>m</mi><mo>+</mo><mi>&Delta;</mi><mi>&beta;</mi><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>,</mo><mi>n</mi><mo>)</mo><mo>+</mo><msubsup><mi>P</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow><mo>&prime;</mo></msubsup><mrow><mo>(</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>,</mo></mrow></mtd><mtd><mrow><mi>m</mi><mo>&Element;</mo><mi>R</mi><mo>&lsqb;</mo><mn>1</mn><mo>,</mo><msub><mi>n</mi><mrow><mi>l</mi><mi>i</mi><mi>n</mi><mi>e</mi><mi>s</mi></mrow></msub><mo>-</mo><mi>&Delta;</mi><mi>&beta;</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>&rsqb;</mo></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mi>P</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow><mo>&prime;</mo></msubsup><mrow><mo>(</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>,</mo></mrow></mtd><mtd><mrow><mi>m</mi><mo>&Element;</mo><mi>R</mi><mo>&lsqb;</mo><msub><mi>n</mi><mrow><mi>l</mi><mi>i</mi><mi>n</mi><mi>e</mi><mi>s</mi></mrow></msub><mo>-</mo><mi>&Delta;</mi><mi>&beta;</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>+</mo><mn>1</mn><mo>,</mo><msub><mi>n</mi><mrow><mi>l</mi><mi>i</mi><mi>n</mi><mi>e</mi><mi>s</mi></mrow></msub><mo>&rsqb;</mo></mrow></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0001050559620000021.GIF" wi="1534" he="183" /></maths>步骤五、每个帧周期将第一行累加数据输出作为像移自配准的M级积分值,得到第i个帧周期的M级积分图像,用公式表示为:<maths num="0005"><math><![CDATA[<mrow><msub><mi>O</mi><mrow><mi>M</mi><mo>,</mo><mi>i</mi></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>&beta;</mi><mrow><mo>(</mo><mi>i</mi><mo>-</mo><mi>M</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow></munderover><mfenced open = "{" close = "}"><mtable><mtr><mtd><mrow><msub><mi>P</mi><mrow><mi>i</mi><mo>-</mo><mi>M</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>&lsqb;</mo><mi>M</mi><mi>&beta;</mi><mrow><mo>(</mo><mi>i</mi><mo>-</mo><mi>M</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><mi>j</mi><mo>,</mo><mi>n</mi><mo>&rsqb;</mo><mo>+</mo></mrow></mtd></mtr><mtr><mtd><mrow><munderover><mo>&Sigma;</mo><mrow><mi>x</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>M</mi><mo>-</mo><mn>1</mn></mrow></munderover><msubsup><mi>P</mi><mrow><mi>i</mi><mo>-</mo><mi>M</mi><mo>+</mo><mn>1</mn><mo>+</mo><mi>x</mi></mrow><mo>&prime;</mo></msubsup><mo>&lsqb;</mo><mrow><mo>(</mo><mi>M</mi><mo>-</mo><mi>x</mi><mo>)</mo></mrow><mi>&beta;</mi><mrow><mo>(</mo><mi>i</mi><mo>-</mo><mi>M</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><mi>j</mi><mo>,</mo><mi>n</mi><mo>&rsqb;</mo></mrow></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0001050559620000022.GIF" wi="1470" he="278" /></maths>实现在TDI的同时完成像移的自配准。
地址 130033 吉林省长春市东南湖大路3888号