发明名称 一种基于图像对比度的调频斜率估计方法
摘要 本发明属于估计方法,具体涉及一种基于图像对比度的调频斜率估计方法。它包括:步骤一:预处理,计算初始的调频斜率Ka1,Ka2,步骤二:计算调频斜率,由黄金分割法计算得到调频斜率Ka3,并利用这个调频斜率对原始SAR数据进行距离向压缩和方位向压缩,得到二维图像,同时计算所得二维图像的对比度C3;步骤三:曲线拟合,利用上述计算出的三个对应点(Ka1,C1)、(Ka2,C2)和(Ka3,C3)进行二次抛物线拟合,得到拟合曲线,步骤四:计算斜率,步骤五:成像。使用本发明的效果是:该方法是根据图像对比度与方位向调频斜率的关系,使用抛物线法来估计多普勒调频斜率,避免了传统算法中的反复迭代,可以大大提高计算效率。
申请公布号 CN104656067A 申请公布日期 2015.05.27
申请号 CN201310606379.6 申请日期 2013.11.25
申请人 北京环境特性研究所 发明人 徐志明;卢永革;姚京萍
分类号 G01S7/40(2006.01)I 主分类号 G01S7/40(2006.01)I
代理机构 核工业专利中心 11007 代理人 高尚梅;刘昕宇
主权项 一种基于图像对比度的调频斜率估计方法,其特征在于,包括下述步骤:步骤一:预处理根据雷达惯导数据计算初始的调频斜率<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原始数据进行距离向压缩和方位向压缩,得到合成孔径雷达图像。
地址 100076 北京市海淀区永定路50号