主权项 |
TDI成像的像移自配准方法,构建像移自配准算法理论模型,其特征是,该方法由以下步骤实现:步骤一、实时采集当前帧图像P<sub>i+1</sub>(m,n),并计算得到当前帧图像的像移量S<sub>P</sub>(i)对应的像素个数Δn(i),用公式表示为:<maths num="0001"><math><![CDATA[<mrow><mi>Δ</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>′</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>α</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><msubsup><mi>P</mi><mi>i</mi><mo>′</mo></msubsup><mrow><mo>(</mo><mi>m</mi><mo>+</mo><mi>Δ</mi><mi>β</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>[</mo><mn>1</mn><mo>-</mo><mi>α</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>]</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>∈</mo><mi>R</mi><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>-</mo><mi>Δ</mi><mi>β</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>]</mo></mrow></mtd></mtr><mtr><mtd><mrow><mtable><mtr><mtd><mrow><mi>α</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>[</mo><mn>1</mn><mo>-</mo><mi>α</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>]</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>∈</mo><mi>R</mi><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>Δ</mi><mi>β</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>]</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>′</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>Σ</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>[</mo><mn>1</mn><mo>+</mo><munderover><mi>Σ</mi><mrow><mi>x</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></munderover><mi>β</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>]</mo></mrow><mi>l</mi></mfrac><mo>,</mo></mrow></mtd><mtd><mrow><mi>M</mi><mo>≥</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>[</mo><mn>1</mn><mo>+</mo><munderover><mi>Σ</mi><mrow><mi>x</mi><mo>=</mo><mn>2</mn></mrow><mi>M</mi></munderover><mi>β</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>]</mo></mrow><mi>M</mi></mfrac><mo>,</mo></mrow></mtd><mtd><mrow><mi>M</mi><mo><</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>′</mo><mo>′</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>′</mo><mo>′</mo></mrow></msubsup><mrow><mo>(</mo><mi>m</mi><mo>+</mo><mi>Δ</mi><mi>β</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>′</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>∈</mo><mi>R</mi><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>-</mo><mi>Δ</mi><mi>β</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>]</mo></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mi>P</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow><mo>′</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>∈</mo><mi>R</mi><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>Δ</mi><mi>β</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>]</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>Σ</mo><mrow><mi>j</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>β</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>[</mo><mi>M</mi><mi>β</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>]</mo><mo>+</mo></mrow></mtd></mtr><mtr><mtd><mrow><munderover><mo>Σ</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>′</mo></msubsup><mo>[</mo><mrow><mo>(</mo><mi>M</mi><mo>-</mo><mi>x</mi><mo>)</mo></mrow><mi>β</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>]</mo></mrow></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0001050559620000022.GIF" wi="1470" he="278" /></maths>实现在TDI的同时完成像移的自配准。 |