主权项 |
一种基于图像对比度的调频斜率估计方法,其特征在于,包括下述步骤:步骤一:预处理根据雷达惯导数据计算初始的调频斜率<img file="FDA0000422023300000011.GIF" wi="322" he="137" /><img file="FDA0000422023300000012.GIF" wi="307" he="137" />其中,V表示雷达载体的速度,R表示雷达到目标的距离,λ表示雷达波长,并利用这两个调频斜率对原始SAR数据分别进行距离向压缩和方位向压缩,得到二维图像,同时分别计算所得二维图像的对比度<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>C</mi><mn>1</mn><mo>=</mo><mfrac><msqrt><mi>A</mi><mo>{</mo><msup><mrow><mo>[</mo><msup><mi>I</mi><mn>2</mn></msup><mrow><mo>(</mo><mi>n</mi><mo>,</mo><mi>Ka</mi><mn>1</mn><mo>)</mo></mrow><mo>-</mo><mi>A</mi><mo>[</mo><msup><mi>I</mi><mn>2</mn></msup><mrow><mo>(</mo><mi>n</mi><mo>,</mo><mi>Ka</mi><mn>1</mn><mo>)</mo></mrow><mo>]</mo><mo>]</mo></mrow><mn>2</mn></msup><mo>}</mo></msqrt><mrow><mi>A</mi><mo>[</mo><msup><mi>I</mi><mn>2</mn></msup><mrow><mo>(</mo><mi>n</mi><mo>,</mo><mi>Ka</mi><mn>1</mn><mo>)</mo></mrow><mo>]</mo></mrow></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000422023300000013.GIF" wi="795" he="151" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>C</mi><mn>2</mn><mo>=</mo><mfrac><msqrt><mi>A</mi><mo>{</mo><msup><mrow><mo>[</mo><msup><mi>I</mi><mn>2</mn></msup><mrow><mo>(</mo><mi>n</mi><mo>,</mo><mi>Ka</mi><mn>2</mn><mo>)</mo></mrow><mo>-</mo><mi>A</mi><mo>[</mo><msup><mi>I</mi><mn>2</mn></msup><mrow><mo>(</mo><mi>n</mi><mo>,</mo><mi>Ka</mi><mn>2</mn><mo>)</mo></mrow><mo>]</mo><mo>]</mo></mrow><mn>2</mn></msup><mo>}</mo></msqrt><mrow><mi>A</mi><mo>[</mo><msup><mi>I</mi><mn>2</mn></msup><mrow><mo>(</mo><mi>n</mi><mo>,</mo><mi>Ka</mi><mn>2</mn><mo>)</mo></mrow><mo>]</mo></mrow></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000422023300000014.GIF" wi="821" he="155" /></maths>其中为I(n,Ka),(n=1,...,N)方位向一维图像的幅度,A()表示空间求平均运算;步骤二:计算调频斜率由黄金分割法计算得到调频斜率Ka3=0.618(Ka2‑Ka1)+Ka1,并利用这个调频斜率对原始SAR数据进行距离向压缩和方位向压缩,得到二维图像,同时计算所得二维图像的对比度<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>C</mi><mn>3</mn><mo>=</mo><mfrac><msqrt><mi>A</mi><mo>{</mo><msup><mrow><mo>[</mo><msup><mi>I</mi><mn>2</mn></msup><mrow><mo>(</mo><mi>n</mi><mo>,</mo><mi>Ka</mi><mn>3</mn><mo>)</mo></mrow><mo>-</mo><mi>A</mi><mo>[</mo><msup><mi>I</mi><mn>2</mn></msup><mrow><mo>(</mo><mi>n</mi><mo>,</mo><mi>Ka</mi><mn>3</mn><mo>)</mo></mrow><mo>]</mo><mo>]</mo></mrow><mn>2</mn></msup><mo>}</mo></msqrt><mrow><mi>A</mi><mo>[</mo><msup><mi>I</mi><mn>2</mn></msup><mrow><mo>(</mo><mi>n</mi><mo>,</mo><mi>Ka</mi><mn>3</mn><mo>)</mo></mrow><mo>]</mo></mrow></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000422023300000015.GIF" wi="849" he="163" /></maths>步骤三:曲线拟合利用上述计算出的三个对应点(Ka1,C1)、(Ka2,C2)和(Ka3,C3)进行二次抛物线拟合,得到拟合曲线,步骤四:计算斜率计算拟合得到的二次曲线的最大值及所对应的调频斜率值Ka0,该调频斜率值就是调频斜率的估计值,步骤五:成像利用Ka0对SAR原始数据进行距离向压缩和方位向压缩,得到合成孔径雷达图像。 |