基于MIMO毫米波雷达的多模式形变监测系统及方法

基于mimo毫米波雷达的多模式形变监测系统及方法
技术领域
1.本发明涉及形变监测技术领域，尤其涉及一种基于mimo毫米波雷达的多模式形变监测系统及方法。


背景技术：

2.近些年，随着科学技术的不断发展，基建设施得到了快速发展，但相对的，也发生了很多老化建筑坍塌崩裂的悲惨事故，造成了不可轻视的损失和影响。无论是在基础设施的建设或者完工投入使用的过程中，都可能因为各种各样的原因而发生一定程度的形变，而当该形变值超出了该设施的最大可承受范围之后，就可能发生上述坍塌崩裂的事故。因此，就需要对建筑物的微弱形变进行持续有效的监测。
3.但是目前的形变监测方法或只能够用于微弱振动的监测，不能够对大规模坍塌场景进行有效直观的监测，或是能够针对大幅度形变场景进行监测，但不能够获取准确有效的局部形变信息，且成像精度低，难以识别。
4.因此，亟需一种能够适用于多场景且成像精度高的形变监测系统。


技术实现要素：

5.基于此，有必要针对上述技术问题，提供一种基于mimo毫米波雷达的多模式形变监测系统及方法。
6.一种基于mimo毫米波雷达的多模式形变监测系统，包括：精密滑轨、第一mimo雷达、第二mimo雷达和主控计算机；所述第一mimo雷达用于进行sar模式的目标形变监测，所述第一mimo雷达由所述精密滑轨带动对目标发射调频连续波信号，并接收对应的回波信号，将所述回波信号传输至所述主控计算机；所述主控计算机接收所述回波信号进行处理，获取天线通道数*方位位置*高度位置*距离维的四维立方体离散数据，采用ram成像算法对所述四维立方体离散数据进行处理，获取不同时刻的sar主副图像，对所述sar主副图像进行干涉、配准及解缠绕处理，计算出目标的形变值；所述第二mimo雷达固定在系统最高处，用于进行rar模式的目标微形变监测，通过所述第二mimo雷达对目标连续发射射频连续波，并接收回波，通过所述主控计算机处理所述回波，提取回波的相位信息进行解缠绕，反演得到形变-时序图，根据所述形变-时序图计算出目标的形变值。
7.在其中一个实施例中，所述第一mimo雷达和第二mimo雷达的前端模块均包括有发射天线阵列、接收天线阵列、信号处理存储模块、射频控制电路、功率放大器、混频器、低噪声放大器、中频滤波器和模数转换器；所述发射天线阵列与所述功率放大器连接，所述射频电路连接所述混频器、功率放大器和信号处理存储模块；所述接收天线阵列与所述低噪声放大器连接，所述低噪声放大器与所述混频器连接，所述混频器与所述中频滤波器连接，所述中频滤波器与所述模数转换器连接，所述模数转换器与所述信号处理存储模块连接。
8.一种基于mimo毫米波雷达的多模式形变监测方法，采用如上所述的基于mimo毫米波雷达的多模式形变监测系统，能够实现sar模式和/或rar模式的目标形变监测，其中，sar
模式包括：根据第一mimo雷达的天线布局，接收回波数据并进行处理，获取通道、方位向、高度向和距离向的四维立方体离散数据，所述四维立方体离散数据包括待测目标的方位、高度和距离信息；根据通道将对应的四维立方体离散数据放置在随滑轨移动时发射调频连续波的位置处，根据等效相位中心原理进行相位补偿，并进行距离向的傅里叶变换，通过聚焦成像获取不同时刻sar图像形成的多幅时序复数图像；从所述多幅时序复数图像中选取一幅图像作为主图像，未选中的作为副图像，采用相关系数法对所述主图像和副图像进行图像配准处理，将配准后的主图像和副图像共轭相乘，获得若干幅干涉图像；采用振幅均值法和振幅离差法在所述若干幅干涉图像选取永久散射点，根据所述永久散射点构建三角网络；采用基于fft的最小二乘法对所述三角网络进行相位解缠绕操作，获取每个所述永久散射点的真实相位，根据所述真实相位反演获取待测目标的形变值。
9.在其中一个实施例中，所述根据通道将对应的四维立方体离散数据放置在随滑轨移动时发射调频连续波的位置处，根据等效相位中心原理进行相位补偿，并进行距离向的傅里叶变换，通过成像聚焦获取不同时刻sar图像形成的多幅时序复数图像，具体包括：在慢时间t时刻，接收阵元n接收来自发射阵元m的静态场景目标回波信号的实际多普勒信号相位为：
[0010][0011]
所述接收阵元和发射阵元的中间位置为等效相位中心，等效多普勒信号相位为等效阵元到目标点距离re(t)的双程延时得到的相位，为：
[0012][0013]
将所述实际多普勒信号相位与等效多普勒信号相位的差值进行等效中心相位补偿，将四维立方体离散数据根据等效相位中心原理进行相位补偿，并进行距离向的傅里叶变换，则雷达回波表示为：
[0014][0015]
将相位补偿后的四维立方体离散数据送入rma快速层析算法中，通过成像聚焦形成高精度sar图像，为：
[0016][0017]
其中，σ(x,y,z0)为雷达散射截面积，是目标在第一mimo雷达接收方向上反射第一mimo雷达信号能力的度量，k＝2πf/c为对应发射频率的波数，z0＝z+r0为目标距离第一mimo雷达孔径的距离；根据不同时刻的高精度sar图像，形成时序复数图像。
[0018]
在其中一个实施例中，所述rma快速层析算法具体包括：将任意散射点的散射系数用卷积表示为：
[0019]
[0020]
根据卷积定理将上式改写为：
[0021][0022]
结合驻定相位原理求解上式，则成像区域任意散射点强度分布为：
[0023][0024]
其中，k＝2πf/c为对应发射频率的波数，z0＝z+r0为目标距离雷达孔径的距离。
[0025]
在其中一个实施例中，所述从所述多幅时序复数图像中选取一幅图像作为主图像，未选中的作为副图像，采用相关系数法对所述主图像和副图像进行图像配准处理，将配准后的主图像和副图像共轭相乘，获得若干幅干涉图像，具体包括：在所述主图像中选取主图像中心点，并确定所述主图像中心点在所述副图像上的同名点；根据所述主图像中心点和同名点，计算副图像相对主图像在行、列方向上的坐标偏移量；以所述主图像中心点为中心选取匹配窗口，在所述副图像对应的位置选取搜索框，并以整像素为步进搜索计算不同偏移时的相关系数，获取最大相关系数，所述最大相关系数对应的点为配准点，其中，所述相关系数的计算公式为：
[0026][0027]
式中，m1表示主图像，m2表示副图像，m、n为相关计算的窗口大小，u、v为窗口的偏移量，/表示复共轭；根据所述配准点对所述副图像进行配准，所述副图像逐一配准完成后，将所述副图像分别与所述主图像共轭相乘，获取干涉图像。
[0028]
在其中一个实施例中，所述采用振幅均值法和振幅离差法在所述若干幅干涉图像选取永久散射点，根据所述永久散射点构建三角网络，具体包括：计算所述干涉图像的振幅均值，公式为：
[0029][0030]
式中，m和n分别表示图像中行和列的像元数目；遍历所述干涉图像上所有像素点的幅值，保留所述干涉图像上幅值大于所述振幅均值的像素点，获取若干幅新干涉图像；在所述若干幅新干涉图像中，用(i,j)表示像素点的坐标，ma(i,j)表示振幅均值，ma(i,j)表示方差，其中，振幅均值以及方差的计算公式分别为：
[0031][0032]
式中，ak(i,j)表示第k幅干涉图像的振幅图，n表示振幅图个数；设置阈值ε，计算振幅离差指数da(i,j)，公式为：
[0033][0034]
若所述振幅离差指数大于阈值ε，则认定对应的像素点为永久散射点；基于delaunay三角网的构造方法，根据所述永久散射点构建三角网络。
[0035]
在其中一个实施例中，所述采用基于fft的最小二乘法对所述三角网络进行相位解缠绕操作，获取每个所述永久散射点的真实相位，根据所述真实相位反演获取待测目标的形变值，具体包括：在所述三角网络中选取一个永久散射点p
i,j
作为参考点，计算p
i,j
的值，公式为：
[0036][0037]
设定干涉相位矩阵大小为m
×
n，根据周期函数对每一行的p
i,j
进行镜像对称操作，并按列执行相同操作得到其中，周期函数为：
[0038][0039]
对作二维傅里叶变换，得到p
k,l
，根据下式计算
[0040][0041]
对作二维傅里叶逆变换，得到解缠绕函数的最小二乘估算值，完成解缠绕操作，得到真实相位；基于所述真实相位采用下式反演出待测目标的形变值：
[0042][0043]
其中，λ为电磁波波长，为干涉相位差。
[0044]
在其中一个实施例中，所述rar模式包括：通过第二mimo雷达向待测目标发射调频连续波信号，并接收待测目标的回波信号；对所述回波信号进行距离维的快速傅里叶变换，获取待测目标的距离信息，并根据所述距离信息，确定待测目标所在的距离门并提取对应回波信号的相位值；重复获取所有回波信号的距离信息和相位值，得到相位-时序数据；对所述相位-时序数据进行相位解缠绕处理，获取真实相位；对所述真实相位进行相位差分操
作，并根据差分后的真实相位计算获取待测目标的形变值。
[0045]
相比于现有技术，本发明的优点及有益效果在于：能够实现sar模式和/或rar模式的目标形变监测，从而兼顾了微形变监测和大幅度形变监测，使得系统能够适用于多种场景下的目标形变监测，此外，还能够提供对待测目标的实时、长时间和全方位的有效监测，且提升了雷达图像的精确度，使得形变测量能够达到毫米级，监测精度高。
附图说明
[0046]
图1为一个实施例中一种基于mimo毫米波雷达的多模式形变监测系统的结构示意图；
[0047]
图2为一个实施例中一种基于mimo毫米波雷达的多模式形变监测系统的外观示意图；
[0048]
图3为一个实施例中第一雷达前端模块的结构示意图；
[0049]
图4为一个实施例中一种基于mimo毫米波雷达的多模式形变监测方法的流程示意图；
[0050]
图5为一个实施例中sar模式监测方法的流程示意图；
[0051]
图6为一个实施例中sar模式的运行示意图；
[0052]
图7为一个实施例中rar模式监测方法的流程示意图。
具体实施方式
[0053]
为了使本发明的目的、技术方案及优点更加清楚明白，下面通过具体实施方式结合附图对本发明做进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
[0054]
为了方便理解，下面对具体实施方式中使用的相关专业术语进行说明。
[0055]
ground-based synthetic aperture radar(gbsar)：地基合成孔径雷达，应用在地表微形变监测领域。其不需要布点作业、监测效率高、携带方便、安装灵活、监测范围广，可以满足基础设施建设、建筑结构监测的应用需求。
[0056]
multiple input-multiple output(mimo)：多输入多输出，mimo雷达指在发射端和接收端都使用多阵元的天线阵。mimo雷达系统可以构成远多于实际收发阵元数目的观测信道，多观测信道使得mimo雷达能够采集携带有目标不同幅度、时延和相位的回波信息，实现更多的信息自由度。
[0057]
rar(real aperture radar，真实孔径雷达)模式，具有便于监测，持续时间长的优点，能够用于器械等的微弱振动监测。
[0058]
sar(synthetic aperture radar，合成孔径雷达)模式，能够适用于山体滑坡、大坝坍塌等大幅度形变场景中。
[0059]
如图1至图3所示，提供了一种基于mimo毫米波雷达的多模式形变监测系统10，包括：主控计算机11、精密滑轨12、第一雷达13和第二雷达14；第一mimo雷达13用于进行sar模式的目标形变监测，第一mimo雷达13由精密滑轨12带动对目标发射调频连续波信号，并接收对应的回波信号，将回波信号传输至主控计算机11；主控计算机11接收回波信号进行处理，获取天线通道数*方位位置*高度位置*距离维的四维立方体离散数据，采用ram成像算
法对四维立方体离散数据进行处理，获取不同时刻的sar主副图像，对sar主副图像进行干涉、配准及解缠绕处理，计算出目标的形变值；第二mimo雷达14固定在系统最高处，用于进行rar模式的目标微形变监测，通过第二mimo雷达14对目标发射射频连续波，并接收回波，主控计算机处理回波，提取回波的相位信息进行解缠绕，反演得到形变-时序图，根据形变-时序图计算出目标的形变值。
[0060]
在一个实施例中，系统的外观如图2所示，在使用前，初始化数据采集的各个参数，雷达移动速度、雷达脉冲发射频率等，通过该系统进行sar模式的形变监测时，设置系统x轴(方位向)和y轴(高度向)每次移动距离分别为dx(mm)、dy(mm)，y轴上升高度；雷达采用一发四收的收发方式，并将天线阵元上接收的回波数据全部存储起来，用一个小天线沿着方位向长线阵的轨迹等速移动并辐射相参信号，就能获得一个相当于很长线阵的方位向高分辨力，本系统按照“蛇”型依次按行接收并存储回波信号，将上位机接收到的数据进行预处理，提取天线通道数
×
方向维
×
距离维
×
高度维的四维立方体离散数据，便于后续的数据处理。
[0061]
在本实施例中，系统在sar模式的工作方法为：首先，第一mimo雷达13通过其mimo天线随着精密滑轨12在1m*0.5m的矩形空间里的每一个位置上发射调频连续波信号并接收对应位置的回波信号；其次，将接收到的雷达回波传送到主控计算机11进行数据的处理，形成天线通道数*方位位置*高度位置*距离维的四维矩阵；第三，将获取的四维矩阵送入到rma成像算法中形成高精度sar图像，通过长时间的采集，获得不同时刻的sar主副图像集；第四，进行sar主副图像的配准并共轭相乘获得干涉图，并对干涉图进行相位解缠获取真实的相位差；最后，反演出形变信息在主控计算机11上进行显示。
[0062]
在本实施例中，系统在rar模式的工作方法为：在将第二mimo雷达14调整至最优角度后，对待测目标进行连续的发射调频连续波，接收回波，通过主控计算机11提取每个回波的相位信息来进行解缠，最后反演出形变值组成形变-时序图，在主控计算机11进行显示。其中，第二mimo雷达14固定在系统最高处，可进行监测角度的调整，便于将第二mimo雷达14调整到最优角度对待测目标进行监测。
[0063]
其中，第一mimo雷达13和第二mimo雷达14均包括其对应的雷达子系统，通过雷达子系统能够向待测目标发射调频连续波信号，并接收每个位置的回波信号。
[0064]
其中，精密滑轨12为二维滑轨，能够带动第一mimo雷达13滑动至矩形空间内的任一位置，增大了电磁波的辐射范围，便于获取对应的方位距离图像，并根据方位距离图像获取待测目标的方位和距离信息，从而提升对待测目标的形变监测精确度。
[0065]
在一个实施例中，如图3所示，第一mimo雷达13和第二mimo雷达14的前端模块包括有发射天线阵列、接收天线阵列、信号处理存储模块、射频控制电路、功率放大器、混频器、低噪声放大器、中频滤波器和模数转换器；发射天线阵列与功率放大器连接，射频电路连接混频器、功率放大器和信号处理存储模块；接收天线阵列与低噪声放大器连接，低噪声放大器与混频器连接，混频器与中频滤波器连接，中频滤波器与模数转换器连接，模数转换器与信号处理存储模块连接。
[0066]
具体地，通过发射天线阵列向待测目标发射调频连续波信号，并通过接收天线阵列接收对应的回波信号，回波信号依次进行低噪声放大、混频、中频滤波和模数转换，最后输入信号处理存储模块进行处理或存储，结合主控计算机11，得到待测目标的形变值，从而
实现对待测目标的形变监测。
[0067]
在本实施例中，结合了sar技术和rar技术，形成了一种基于mimo毫米波雷达的多模式形变监测系统，从而兼顾了微形变监测和大幅度形变监测，使得系统能够适用于多种场景下的目标形变监测，此外，能够提供对待测目标的实时、长时间和全方位的有效监测，且提升了雷达图像的精确度，使得形变测量能够达到毫米级，监测精度高。
[0068]
在一个实施例中，如图4所示，提供了一种基于mimo毫米波雷达的多模式形变监测方法，基于上述基于mimo毫米波雷达的多模式形变监测系统，能够实现sar模式和/或rar模式的目标形变监测。在采用sar模式对目标进行形变监测时，将雷达接收到的雷达原始数据进行处理和加窗滤波，并输入rma算法成像，获取二维雷达的高精度图像，多次获取高精度图像并进行共轭相乘，得到干涉图，将干涉图进行相位解缠和误差校正，获取真实相位，并根据真实相位反演得出目标的形变信息。
[0069]
在采用rar模式对目标进行形变监测时，将雷达接收到的雷达原始数据进行数据处理和加窗滤波，计算每一帧各chirp(一个脉冲的时间)对应距离的相位值，将每个chirp的相位值组成相位-时序图，并将相位时序进行相位解缠，从而能够实时计算出目标的位移形变信息。
[0070]
在一个实施例中，如图5至6所示，提供了一种基于mimo毫米波雷达的多模式形变监测方法，采用sar模式实现目标的形变监测，包括以下步骤：
[0071]
步骤s501，根据第一mimo雷达的天线布局，接收回波数据并进行处理，获取通道、方位向、高度向和距离向的四维立方体离散数据，四维立方体离散数据包括待测目标的方位、高度和距离信息。
[0072]
具体地，在进行监测之前，初始化数据采集的各个参数，第一mimo雷达的移动速度和脉冲发射频率等，设置gbsar系统x轴(方位向)和y轴(高度向)每次移动距离分别为dx(mm)和dy(mm)，表示x轴和y轴的上升高度。第一mimo雷达采用一发四收的收发方式，并将天线阵元上接收的回波数据全部进行存储，用一个小天线沿着方位向长线阵的轨迹等速移动并辐射相参信号，就能获得一个相当于长线阵的方位向高分辨能力，系统按照“蛇”型依次按行接收并存储回波信号，其按sar模式运行如附图6所示。将接收到的回波数据进行预处理，提取通道
×
方位维
×
距离维
×
高度维的四维立方体离散数据，便于后续的数据处理。
[0073]
步骤s502，根据通道将对应的四维立方体离散数据放置在随滑轨移动时发射调频连续波的位置处，根据等效相位中心原理进行相位补偿，并进行距离向的傅里叶变换，通过成像聚焦获取不同时刻sar图像形成的多幅时序复数图像。
[0074]
具体地，根据通道将对应的四维立方体离散数据放置在随滑轨移动时发射调频连续波的位置处，并根据等效相位中心原理进行相位补偿，在某一时刻t，距雷达r
t
处的目标回波为：
[0075][0076]
式中，a为第一mimo雷达接收信号的幅度。
[0077]
在对四维立方体离散数据进行相位补偿后，进行距离向的傅里叶变换，使用匹配滤波函数与二维波数谱相乘，通过重建图像函数生成高精度sar图像，将不同时刻形成sar图像组成时序复数图像，从而能够获取较大场景下的待测目标的监测图像，便于对大幅度
形变场景进行监测。
[0078]
其中，在慢时间t时刻，接收阵元n接收来自发射阵元m的静态场景目标回波信号的实际多普勒信号相位为：
[0079][0080]
接收阵元和发射阵元的中间位置为等效相位中心，等效多普勒信号相位为等效阵元到目标点距离re(t)的双程延时得到的相位，为：
[0081][0082]
将所述实际多普勒信号相位与等效多普勒信号相位的差值进行等效中心相位补偿，将四维立方体离散数据根据等效相位中心原理进行相位补偿，并进行距离向的傅里叶变换，则雷达回波表示为：
[0083][0084]
将相位补偿后的四维立方体离散数据送入rma快速层析算法中，通过成像聚焦形成高精度sar图像，为：
[0085][0086]
其中，σ(x,y,z0)为雷达散射截面积，是目标在第一mimo雷达接收方向上反射第一mimo雷达信号能力的度量，k＝2πf/c为对应发射频率的波数，z0＝z+r0为目标距离第一mimo雷达孔径的距离；根据不同时刻的高精度sar图像，形成时序复数图像。
[0087]
具体地，将实际多普勒信号相位与等效多普勒信号相位的差值进行等效中心相位补偿，将四维立方体离散数据根据等效相位中心原理进行相位补偿，再进行距离向的傅里叶变换，将相位补偿后的数据送入rma快速层析算法中，形成高精度sar图像，得到n幅不同时刻的高精度sar时序复数图像，从而提升了sar雷达图像的精确度。
[0088]
其中，雷达散射截面积的计算公式为：
[0089][0090]
其中，rma快速层析算法具体包括：将任意散射点的散射系数用卷积表示为：
[0091][0092]
根据卷积定理将上式改写为：
[0093]
[0094]
结合驻定相位原理求解上式，则成像区域任意散射点强度分布为：
[0095][0096]
其中，k＝2πf/c为对应发射频率的波数，z0＝z+r0为目标距离雷达孔径的距离。
[0097]
具体地，为了便于对回波信号进行频域处理，根据卷积定理将散射系数进处理，利用驻定相位原理(principle of stationary phase，posp)求解，由于posp对于幅度变化慢、相位变化快的信号有很好的计算精度，其实现原理是在相位变化很快的地方，幅度在一个完整相位周期内近似为常数，由于相位周期的正负部分相互抵消，故其对于傅里叶变换中的积分运算贡献几乎为零，根据posp求解，得到成像区域任意散射点散射强度分布。使用匹配滤波函数与二维波数谱相乘，通过重建图像函数生成高精度sar图像，将不同时刻形成sar图像组成时序复数图像。
[0098]
步骤s503，从多幅时序复数图像中选取一幅图像作为主图像，未选中的作为副图像，采用相关系数法对主图像和副图像进行图像配准处理，将配准后的主图像和副图像共轭相乘，获得若干幅干涉图像。
[0099]
具体地，获取多幅时序复数图像后，选取一幅图像作为主图像，未选中的作为副图像，将主图像与任意副图像进行配准，首先需要在主图像和副图像之间找到足够多的同名点，即两幅图像之间的相同特征点，通过同名点之间的对应关系，确定主图像与副图像的几何变换模型，根据几何变换关系对副图像进行重采样，获取配准后的图像，根据前述方法逐一对所有副图像进行配准，将配准后的主图像和副图像共轭相乘，获取对应的干涉图像。
[0100]
其中，在主图像中选取主图像中心点，并确定主图像中心点在副图像上的同名点；根据主图像中心点和同名点，计算副图像相对主图像在行、列上的坐标偏移量；以主图像中心点为中心选取匹配窗口，在副图像对应的位置选取搜索框，并以整像素为步进搜索计算不同偏移量时的相关系数，获取最大相关系数，最大相关系数对应的点为配准点，其中，相关系数的计算公式为：
[0101][0102]
式中，m1表示主图像，m2表示副图像，m、n为相关计算的窗口大小，u、v为窗口的偏移量，*表示复共轭；根据配准点对副图像进行配准，副图像逐一配准完成后，将副图像分别与主图像共轭相乘，获取干涉图像。
[0103]
具体地，选取主图像中心点pm(um,vm)，计算点pm在副图像中的像素坐标ps(us,vs)，该点即主图像中心点pm在副图像上的同名点；则副影像相对主影像在行、列方向上的坐标偏移量分别为u＝u
s-um和v＝v
s-vm，以主图像中心点为控制点，确定主图像的控制点后，以该控制点为中心选取匹配窗口，在辅图像的对应位置选取搜索框，并以整像素为步进搜索计算不同偏移时的相干系数，相关系数最大的点即为配准点。图像配准完毕后，副图像分别与主图像共轭相乘，获取对应的干涉图像。
[0104]
步骤s504，采用振幅均值法和振幅离差法在若干幅干涉图像中选取永久散射点，根据永久散射点构建三角网络。
[0105]
具体地，由于永久散射点通常具有较强的散射特性，在采用振幅离差法进行处理之前，可以采用振幅均值阈值法，将振幅较小的点去除，然后采用振幅离差法获取干涉图像中的永久散射点，从而提升筛选永久散射点的精确度，并根据干涉图像的永久散射点构造三角网络，便于获取永久散射点的相位。
[0106]
其中，计算所述干涉图像的振幅均值，公式为：
[0107][0108]
式中，m和n分别表示图像中行和列的像元数目；遍历干涉图像上所有像素点的幅值，保留干涉图像上幅值大于振幅均值的像素点，获取若干幅新干涉图像；在若干幅新干涉图像中，用(i,j)表示像素点的坐标，ma(i,j)表示振幅均值，ma(i,j)表示方差，其中，振幅均值以及方差的计算公式分别为：
[0109][0110]
式中，ak(i,j)表示第k幅干涉图像的振幅图，n表示振幅图个数；设置阈值ε，计算振幅离差指数da(i,j)，公式为：
[0111][0112]
若振幅离差指数大于阈值ε，则认定对应的像素点为永久散射点；基于delaunay三角网的构造方法，根据永久散射点构建三角网络。
[0113]
具体地，根据获取的干涉图像，采用振幅离差法和振幅均值法，选取干涉图像中相位稳定性强、变化小的强散射点，作为永久散射点。由于永久散射点通常具有较强的散射特性，因此，振幅过小的点不会是永久散射点，因此，在采用振幅离差法进行处理之前，可以采用振幅均值阈值法，将振幅较小的点去除，防止误判的发生，提升最终形变监测的精确度。
[0114]
获取干涉图像中的永久散射点后，需要根据永久散射点建立ps基线网络，采用基于delaunay三角网的构造方法，即其中每两个相邻的三角形构成的凸四边形的对角线，在相互交换后，六个内角的最小角不再增大，且任意一个三角形的外接圆不含面内其他点，该三角网络具有结构好、数据冗余小及存储简单等优点，从而能够便于永久散射点的相位获取。
[0115]
步骤s505，采用基于fft的最小二乘法对三角网络进行相位解缠绕操作，获取每个永久散射点的真实相位，根据真实相位反演获取待测目标的形变值。
[0116]
具体地，在根据永久散射点组成三角网络后，采用基于fft的最小二乘法对三角网络进行相位解缠绕，获取每个永久散射点的真实相位，并根据真实相位反演获取各个永久散射点的形变值，从而获取待测目标的形变大小和形变速率，通过sar模式的形变监测，使
得该方法能够适用于大幅度形变场景中。
[0117]
其中，在三角网络中选取一个永久散射点p
i，j
作为参考点，计算p
i，j
的值，公式为：
[0118][0119]
设定干涉相位矩阵大小为m
×
n，根据周期函数对每一行的p
i,j
进行镜像对称操作，并按列执行相同操作得到其中，周期函数为：
[0120][0121]
对作二维傅里叶变换，得到p
k,l
，根据下式计算
[0122][0123]
对作二维傅里叶逆变换，得到解缠绕函数的最小二乘估算值，完成解缠绕操作，得到真实相位；基于真实相位采用下式反演出待测目标的形变值：
[0124][0125]
其中，λ为电磁波波长，为干涉相位差。
[0126]
具体地，ps点组成三角网络后，采用相位解缠绕的方式获得各个ps点的形变大小和形变速率。在网络中，选择一个稳定的点为参考点，利用最小二乘法即可求得每个ps点的差分相位。主要算法思想为相位解缠前的相位数据导数和解缠后的相位数据导数的差最小化，从而得到其解缠相位的值，其相位解缠的数学表达式为：
[0127][0128]
其中，k
i,j
为整数，-π≤ψ(i,j)≤π，i∈[0,m-1],j∈[0,m-1]，某一像元点(i，j)解缠相位的梯度和缠绕相位的梯度之差为:
[0129][0130]
其中，要使得j最小等价于其平
方和最小，对其求平方后对求偏导解得：
[0131][0132]
其中令：
[0133][0134][0135]
求极值问题中，要使j最小，必然有δj＝0，故需满足b
i,j-b
i-1,j
+a
i,j-a
i-1,j
＝0，代入有：
[0136][0137]
ρ(i,j)＝δ
x
(i,j)-δ
x
(i-1,j)-δy(i,j)-δy(i,j-1)
[0138]
其边界条件为：
[0139][0140]
该式就是离散形式的具有纽曼边界的泊松方程，即相位解缠问题就简化为求离散泊松方程的解为：
[0141][0142]
最后，利用解缠的真实相位反演出目标微形变值，从而实现对大幅度形变场景的形变监测，通过采用fft的最小二乘法对三角网络进行相位解缠绕操作，从而获取到每个永久散射点的真实相位，并能够根据真实相位反演获取待测目标的形变值，从而实现对待测目标大幅度形变的精准监测，提升了对待测目标形变监测的精确度。
[0143]
在本实施例中，根据第一mimo雷达的天线布局，接收回波数据并进行处理，得到四维立方体离散数据，基于通道将对应的四维立方体离散数据放置在随滑轨移动时发射射频连续波的位置处，根据等效相位中心原理进行相位补偿，并进行距离向为傅里叶变换，得到不同时刻的sar图像组成的多幅时序复数图像；在多幅时序复数图像中确定主图像和副图像，采用相关系数法对主图像和副图像进行图像配准，并在配准后进行共轭相乘，得到若干幅干涉图像；采用振幅均值法和振幅离差法在若干幅干涉图像中选取永久散射点，并构建三角网络，采用基于fft的最小二乘法对三角网络进行相位解缠绕操作，获取所有永久散射点的真实相位，并反演得到待测目标的形变值，实现了sar模式的目标形变监测，从而能够实现对大幅度形变场景的监测，且提升了获取的sar图像的精度，从而能够得到更加清晰的成像图，便于精准识别形变信息。
[0144]
在一个实施例中，如图7所示，提供了一种基于mimo毫米波雷达的多模式形变监测方法，采用rar模式实现目标的形变监测，包括以下步骤：
[0145]
步骤s701，通过第二mimo雷达向待测目标发射调频连续波信号，并接收待测目标的回波信号。
[0146]
具体地，确定待测目标后，将系统中的相关参数初始化，启动系统，将第二mimo雷达调整到监测待测目标的最优角度，通过设置在系统最高处的第二mimo雷达向待测目标发射调频连续波信号，并接收对应的回波信号，获取原始回波数据。
[0147]
步骤s702，对回波信号进行距离维的快速傅里叶变换，获取待测目标的距离信息，并根据距离信息，确定待测目标所在的距离门并提取对应回波信号的相位值。
[0148]
具体地，对回波信号进行预处理后，再进行距离维的快速傅里叶变换，在某时刻t，距雷达rt处的目标回波为：
[0149][0150]
式中，a为雷达接收信号的幅值；此时，根据第二mimo雷达与待测目标的实际距离关系，确定待测目标所在的距离门信息和在距离维的频点数，对待测目标所在的频点位置处的复数值进行相位求解操作，得到该点回波信号的相位值，为：
[0151][0152]
步骤s703，重复获取所有回波信号的距离信息和相位值，得到相位-时序数据。
[0153]
具体地，通过去斜处理消除上述相位中最后一项非线性残余相位，同时由于电磁波的传输速度快，待测目标移动前后，信号在回波和第二mimo雷达之间的传输时间近似异质，所以电磁波传输时间差值近似为0，此时相位简化为：
[0154][0155]
对接收的每一个回波信号重复进行步骤s702的操作，得到相位-时序图，从而能够根据相位-时序图提取对应的相位变化关系。
[0156]
步骤s704，对相位-时序数据进行相位解缠绕处理，获取真实相位。
[0157]
具体地，在获取相位-时序数据后，对其进行解缠绕处理，解缠绕的步骤与sar模式中步骤s505的解缠绕方法相同，此处不再赘述。由于相位值只能在[-π，π]之间，因此每当连续值之间的相位差大于或者小于
±
π时，将相位减去2π获得真实相位。
[0158]
步骤s705，对真实相位进行相位差分操作，并根据差分后的真实相位计算获取待测目标的形变值。
[0159]
具体地，在真实相位后，对其进行相位差分操作，当r
t
不同时，求得两次回波之间的相位差与距离差之间的关系为：
[0160][0161]
其中，λ为电磁波波长，为干涉相位差，从而求得目标的微弱形变值d
los
。通过上
述步骤实现了rar模式的形变监测，从而能够对生产器械的微弱振动进行实时持续的监测，且形变测量能够达到毫米级，监测精度高。
[0162]
在本实施例中，通过第二mimo雷达向待测目标发射调频连续波信号，接收对应的回波信号，对回波信号进行fft处理，获取待测目标的距离信息，从而确定待测目标所在的距离门，提取对应的相位，重复获取回波信号的距离信息和相位值，形成相位-时序数据，并进行相位解缠绕处理，获取真实相位，通过对真实相位进行相位差分操作，计算得到待测目标的形变值，实现了通过rar模式对目标的形变监测，能够对微弱振动进行实时持续监测，且监测精度高。
[0163]
在一个实施例中，提供了一种基于mimo毫米波雷达的多模式形变监测方法，能够通过同时启动sar模式和rar模式，对目标进行形变监测，使得该方法能够适用于多种场景的目标形变监测，能够同时对大幅度形变和微弱形变进行监测，便于提取待测目标的有效形变信息，且能够得到更加清晰的成像图，提升了待测目标的监测精确度。
[0164]
显然，本领域的技术人员应该明白，上述本发明的各模块或各步骤可以用通用的计算装置来实现，它们可以集中在单个的计算装置上，或者分布在多个计算装置所组成的网络上，可选地，它们可以用计算装置可执行的程序代码来实现，从而，可以将它们存储在计算机存储介质(rom/ram、磁碟、光盘)中由计算装置来执行，并且在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤，或者将它们分别制作成各个集成电路模块，或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。所以，本发明不限制于任何特定的硬件和软件结合。
[0165]
以上内容是结合具体的实施方式对本发明所做的进一步详细说明，不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干简单推演或替换，都应当视为属于本发明的保护范围。