SAR成像系列:【15】合成孔径雷达(SAR)运动补偿

2023-11-17

       不同于光学图像,SAR图像的获取的本质是方位信号的累积结果,也就是说是在合成孔径时间内的电磁波照射结果。类似于光学成像中的延时拍照。通常使用手机拍照时,若拍照的手臂出现抖动,那么得到的照片就会变模糊。同样的现象也会出现在SAR图像中。SAR的成像时,雷达天线是安装在各种运动平台上的,如卫星、飞机、无人机、飞艇、导弹等,这些平台在实际的飞行过程中,受环境及动力影响,通常无法保持匀速直线运动,这将导致在进行方位压缩时存在相位误差,该误差会影响成像质量,使成像效果模糊(模糊的本质是方位匹配滤波函数与实际的方位信号产生失配)。因此运动补偿是SAR成像不可缺少的环节。

        其实,由平台运动引起的运动误差是SAR“相位误差”中的一个主要部分。在相位误差中除了运动误差(需要运动补偿进行消除),还有大气传播误差、地面曲率误差等。这些误差在不同的场景中要具体考虑。这里主要讨论由平台运动引起的相位误差补偿,即运动补偿。

        运动补偿主要分为两类:基于平台运动数据的粗补偿和基于回波数据的精细补偿。由于飞行平台上的运动数据主要来源于导航仪(如 GPS、北斗 等),它的精度不高,只能补偿掉较大的运动误差。要获得高质量的SAR图像就必须从回波数据中估计多普勒参数,进而实现SAR图像的精确聚焦。

        在多普勒参数中,多普勒中心影响图像的位置,而多普勒调频率影响SAR图像的模糊 程度。在运动补偿中需要对多普勒调频率进行精确计算从而消除图像模糊现象。主要方法由:图像偏置法,最大对比度自聚焦算法,最大熵聚焦算法和相位梯度法等。下面介绍几种常用的运动补偿方法。

1.相位梯度自聚焦算法-PGA

        PGA算法是各种运动平台最常用的运动补偿算法。它不是基于相位误差估计模型,而是基于 SAR 图像中的强散射点,直接进行运动误差估计,可以对任意阶数的多普勒误差进行估计。

        PGA 算法主要包括四步,圆周移位、加窗、相位梯度估计和迭代相位校正。

①圆周位移

        从每个距离单元中选取最强散射点(SAR图像较大时,可只选取一部分距离单元),并将其圆周移位到图像中心,同时也去除了散射点的多普勒中心频率。圆周移位不仅仅是对最强散射点进行排列,同时也去除了多普勒频率,只保留下了相位误差信息。
②加窗

        对圆周移位后的图像域数据进行加窗,提高信噪比,减少背景杂波的影响,提高相位误差估计的准确度。加窗的宽度可根据图像的模糊程度来确定。

③相位梯度估计

        相位梯度估计是 PGA 算法的核心部分,它是在距离多普勒域完成的。因此首先应该对加窗的图像数据进行IFFT,再计算两个相位点的相位误差梯度:

         Sr为加窗内的IFFT信号。d(n)为

         将所有的相位点的相位误差梯度相加可得到总的相位误差。

④迭代相位校正

        利用相位误差对 SAR 距离多普勒数据进行补偿。并将补偿过的数据FFT转换到图像域,用作下一次迭代。反复迭代几次(一般4-6次即可),图像会越来越聚焦。

        下图为PGA聚焦的一幅星载SAR图像。

2、最大对比度算法(与最大熵算法类似)

        当平台以理想状态下运动时,进行方位维压缩之后,所得到的目标聚焦函数应该是 sinc 函数,匹配后的图像能量集中,目标与背景的对比度高;当飞行平台不能按照理想状况运动时,此时匹配滤波函数与实际的方位信号存在存在多普勒匹配误差,导致信号主瓣展宽,旁瓣升高,目标与背景部分的对比度下降。

        所以,以信号的对比度位指标可以对 SAR 成像质量进行评估,当图像的对比度最大时,SAR 成像的方位匹配压缩效果应为最好。因此,SAR图像的运动补偿问题就转换为求取图像对比度最大值的问题,就是在匹配函数中搜索最优的多普勒调频率。

      关于怎样搜索调频率的值是一类搜索策略问题。以进退法为例,其运动补偿过程如下:

①根据平台运动数据,设计方位匹配滤波函数并对回波数据进行聚焦处理;

②计算所得聚焦图像的对比度;

③确定调频率迭代步长,使用迭代调频率聚焦的图像计算对比度,并与②比较;

④根据比较结果修正迭代步长,如果新计算对比度大于原对比度,则迭代值不变,继续迭代,否则令迭代步长变号,并将迭代步长绝对值减半,继续迭代。同时计算图像对比度;

⑤重复③、④直到迭代步长小于某值,则认为找到了正确的调频率。

本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)

SAR成像系列:【15】合成孔径雷达(SAR)运动补偿 的相关文章

随机推荐

  • 3、backbone中的model实例

    关于backbone 最基础的一个东西就是model 这个东西就像是后端开发中的数据库映射那个model一样 也是数据对象的模型 并且应该是和后端的model有相同的属性 仅是需要通过前端来操作的属性 下面就从实例来一步一步的带大家来了解b
  • 文心一言、讯飞星火与GPT-4/3.5在回答中文历史问题的表现

    最近 随着备受关注的文心一言正式免费向全社会开放 再次引起了社会层面对国产大模型的兴趣 以文心一言为代表的国产大模型性能究竟如何 如果将它们相互比较 并且和GPT系列模型等国际前沿水平的LLM进行比较 会得到什么样的结果呢 笔者对此非常好奇
  • Python所有方向的学习路线图

    在放这个路线图之前 我先做一个简单的介绍 避免新手小白有什么不理解的 这个学习路线上面写的是某个方向建议学习和掌握的知识带你 这样学习下来之后会更加容易掌握 知识体系会比较全面 比自己单纯的自学效果好很多 不至于看到什么就学什么 容易走弯路
  • 二,ES6中新增const关键字的使用方法

    之前用var声明变量 变量想怎么改就怎么改 这里const关键字也是声明变量的 不过声明的是常量 常量就是固定的一个值 不能改变 例如 const name 唐僧 name 老沙 报错 因为它要更改常量name 只在块级作用于起作用 和le
  • Windows10 java环境变量的配置详细教程(Windows10 和Windows11)

    java环境变量的配置详细教程 1 首先要区分一下Windows10 2021年之前的版本和Windows10 2021年之后的版本 Windows10 2021年之后的版本和Windows11 系统在配置java上差不多 故不作区分 1
  • Intellij IDEA svn的使用记录

    这里的忽略一直灰色的 可以进入 这里的版本控制里进行忽略选择 或者 这里进行添加 这里有三个选择 按照顺序 1 忽略指定的文件 2 忽略文件夹下所有文件 3 忽略符合匹配规则的文件 到Commit Changes 这里有几个选项需要了解的
  • Vue实例挂载的过程

    一 思考与分析 我们都听过知其然知其所以然这句话 那么不知道是否思考过new Vue 这个过程中究竟做了些什么 过程中是如何完成数据的绑定 又是如何将数据渲染到视图的等等 首先找到vue的构造函数 源码位置 src core instanc
  • 因为强行关机, 而导致的fedora23 不能重新启动, 卡在开机logo那里的 修复 解决方案...

    其实 fedora23的U盘live 也很好用 很流畅 主要还是 要用一个比较好的 快的 U盘 这样live U盘在4GB 3 75GiB 的内存中还是较快的 原来的U盘live系统用得很卡 可能是因为 U盘太烂的原因 要方便的使用live
  • SAXParserFactoryImpl cast SAXParserFactory异常

    Caused by java lang ClassCastException com sun org apache xerces internal jaxp SAXParserFactoryImpl cannot be cast to ja
  • web攻击日志分析之新手指

    0x00 前言 现实中可能会经常出现web日志当中出现一些被攻击的迹象 比如针对你的一个站点的URL进行SQL注入测试等等 这时候需要你从日志当中分析到底是个什么情况 如果非常严重的话 可能需要调查取证谁来做的这个事情 攻击流程是什么样子的
  • AlibabaProtect 卸载,不使用其他软件

    背景 发现系统中存在AlibabaProtect服务 停止不掉 文件夹也删除不掉 还占用内存 CPU 在网上也搜了很多其他的步骤 发现不太容易 这是整理的比较简单的 不需要装其他软件 步骤 1 删除注册表 AlibabaProtect搜索之
  • 华为OD机试-最长连续方波信号

    Online C compiler to run C program online include
  • docker-engine安装

    最近一直在使用docker 做一些试验 每个新机器都需要部署docker的环境 环境信息如下 RedHat 7 2 安装 docker官方的安装 docker engine 1 sudo rpm import https sks keyse
  • Webpack5优化之提高代码运行性能(Preload、Network Cache、Core-js、PWA)

    文章目录 一 Preload Prefetch 1 1 为什么 1 2 是什么 1 2 1 共同点 1 2 2 区别 1 2 3 问题 1 2 4 总结 1 3 怎么样 1 3 1 安装依赖 1 3 2 配置 1 3 3 测试 二 Netw
  • python获取微信群消息_python-itchat 统计微信群、好友数量,及原始消息数据的实例...

    coding utf 8 import itchat from itchat content import TEXT from itchat content import import sys import time import re r
  • LeetCode Week 4

    LeetCode Week 4 练腿是最虐的项目 没有之一 问题集合 1 Reverse Words in a String III Easy 557 Given a string you need to reverse the order
  • 如何高效安装MindSpore的GPU版本

    作者 王磊 更多精彩分享 欢迎访问和关注 https www zhihu com people wldandan MindSpore的GPU版本以前的安装指南 只写清楚了安装依赖 但没有明确指出安装具体执行的命令 缺乏实操性 比较依赖开发者
  • 整流七 - 三相PWM整流器—公式推导篇

    此篇文章为了进一步理解三相pwm整流器 前期的四象限产品 以及仿真模型都没有彻底理解三相pwm整流器的数学模型 于是现在开始一步步推到整流器各个环节的数学公式 三相PWM整流器拓扑结构 补充 三相 VSR 中 交流侧电感的设计尤为重要 起到
  • 面积积分_A-level数学:必考题型之积分求面积解题技巧汇总!!!

    对于A level 数学的pure Math考试部分 总有那么一道积分求面积的题 很多同学在做这种题的时候总是觉得即使自己充分调动学过的所有公式都无法做出来 导致失分 那么今天潘老师带大家一起总结一下 1积分的物理意义 我们知道积分其实是微
  • SAR成像系列:【15】合成孔径雷达(SAR)运动补偿

    不同于光学图像 SAR图像的获取的本质是方位信号的累积结果 也就是说是在合成孔径时间内的电磁波照射结果 类似于光学成像中的延时拍照 通常使用手机拍照时 若拍照的手臂出现抖动 那么得到的照片就会变模糊 同样的现象也会出现在SAR图像中 SAR