一种激光反射断层成像多角度激光回波配准方法

技术领域

本发明属于激光雷达成像技术领域，尤其涉及一种激光反射断层成像多角度激光回波配准方法。

背景技术

激光反射断层成像是一种兼顾远距离和高分辨率成像的新型激光雷达成像技术，最早由Parker J K于1988年提出(Parker J K , Craig E B , Klick D I , et al.Reflective tomography: images from rangeresolved laser radar measurements[J].Appl Opt, 1988, 27(13): 2642-2643)。这项技术是通过激光雷达探测目标的多个角度并收集回波信号，获取目标多角度的深度信息，利用成像算法来计算重构目标的断层面轮廓图像。但在实际工作过程中，由于目标或探测器的抖动，激光脉冲发射和探测器接收的随机抖动等，会在时间维度上出现信号漂移，影响重构图像的准确性，因此需要将采集到的回波信号按照统一的时间基准重新排列，即对各个角度的激光回波进行配准处理。

2、与本专利最接近的现有技术

现有的配准方法主要包括参考屏法、特征点跟踪法、分光法和相位恢复算法。参考屏法是在目标前后找到可以参考的时间基准，接收的回波按参考峰统一配准，其要求参考屏始终相对光束不动且最好满足垂直入射条件，反射率要尽量大，不能被目标遮挡或者遮挡住目标。参考屏法简单易操作，但在实际空间目标附近很难有理想的参考屏，因此这种方法主要适用于地面验证实验；特征点跟踪法（Jin X, Sun J, Yan Y, et al. Featuretracking for projection registration in laboratory- scale reflectivetomography laser radar imaging[J]. Optics Communications, 2010, 283(18):3475–3480）是基于目标有反射系数突变的点，其在投影回波上反映为尖峰或低谷，据此进行配准。但是目标表面往往缺少特征点，且要想将特征点分辨出来，需要发射激光的脉宽足够小，这都对此种方法的应用提出了限制；分光法是将发射激光通过分束镜一分为二，一部分照射目标，另一部分返回到探测器，作为配准的基准。这种方法比较准确，但是由于远距离探测回波间隔时间相对投影回波脉宽较长，对探测器的采集时间长度要求较高。相位恢复算法（赵楠翔, 胡以华. 激光反射层析成像相位恢复算法研究[J]. 红外与激光工程,2019, 48(10)）是通过已知衍射强度分布与目标振幅分布函数，在频域和空间进行反复傅里叶正逆变换，逐步恢复目标的相位信息。通过反复光强迭代消除目标随机运动造成的相位误差，达到减少重构误差、恢复目标图像的目的。但实际重构结果表明其主要消除了重构图像的伪迹，从而改善了重构图像质量，并不能完全消除多角度激光回波数据与投影中心之间存在的相对位移。

图1给出了一种现有的反射投影成像投影图中心的配准方法。

该方法是先利用特征点对反射投影中心进行粗对准，再使用相位恢复算法对投影中心进行进一步对准的方法，这种方法原理可靠，易于实现投影中心的对准，可以有效提高成像图像质量，减少外界因素的影响，更加有效及快速的实现反投影成像。是目前实现激光反射断层成像多角度激光回波配准后重建图像质量提升较为显著的方法。

2.2现有技术的缺点

实际空间目标表面往往缺少特征点，且要想从空间目标中将特征点分辨出来，尤其是例如厘米级空间碎片这类尺度10cm内的微小空间目标，这需要发射激光的脉冲宽度足够小，而能够实现天基远距离探测的高功率脉冲激光器难以达到这样的数量级，这些因素都对利用特征点对反射投影中心进行粗对准的应用提出了限制，且滤波反投影和相位恢复算法实际重构结果对比如图2所示，可以看到相位恢复算法主要消除了重构图像的伪迹，从而改善了重构图像质量，并不能完全消除多角度激光回波数据与投影中心之间存在的相对位移，且重建过程中会损失一部分的目标信息。

因此当前对基于激光反射断层成像空间背景下远距离的空间目标探测，尚缺乏原理切实可行、适用范围广、简单易操作的多角度激光回波配准方法。

发明内容

为解决上述技术问题，本发明提出了一种激光反射断层成像多角度激光回波配准方法，该方法用于激光雷达成像中对激光回波执行配准处理，所述方法包括以下步骤：

步骤1，从所述激光雷达提取该雷达接收的多个角度的雷达回波信号，转到步骤2；

步骤2，对接收的所述雷达回波信号使用滤波反投影算法复原目标图像得到重构图像，并对重构图像进行阈值化处理；转到步骤3；

步骤3，根据阈值化处理结果给出直线初始估计，由这几条估计直线围成闭合目标轮廓，构建估计图像；转到步骤4；

步骤4，计算重建图像与估计图像的误差函数，改变直线参数直到误差函数最小；转到步骤5；

步骤5，确定各条直线的参数，提取出其围成封闭图形的轮廓，得到各角度估计投影的峰值位置坐标；转到步骤6；

步骤6，按照每个投影峰值相对拟合结果的偏差，进行配准；转到步骤7

步骤7，用配准后的投影数据再次重建目标图像。

进一步的，在步骤2中所述的滤波反投影算法采用R-L函数滤波反投影算法。

进一步的， 在步骤2中所述的重构图像进行阈值化处理过程是采用迭代阈值法确定的最佳阈值

步骤2.1，先求出图像的最大灰度值T

步骤2.2，以平均值

步骤2.3，重复该流程直到满足：

步骤2.4，以此时的

进一步的，步骤3中所述的直线初始估计包括以下子步骤：

步骤3.1,计算图像中所有点的梯度大小和方向，然后将梯度方向变化小且相邻的点作为一个连通域，根据每一个连通域的矩形度判断是否需要按照规则将其断开以形成多个矩形度较大的连通域，最后对生成的所有的连通域做改善和筛选，保留其中满足条件的连通域，即为最后的直线检测结果；

步骤3.2，由直线段L

步骤3.21，求出各线段所在直线的交点坐标，然后将相距小于距离分辨率的交点坐标合并，取其坐标平均值形成新的交点坐标作为封闭图形的顶点；

步骤3.22，依次连接各新的交点坐标形成首位连接的直线段

进一步的，步骤4中所述计算重建图像与估计图像的误差函数，其包括以下子步骤：

选取归一化平均绝对值距离测量值r来评价图像重建误差大小，r的表达式为

其中t

进一步的，在步骤6的配准过程包括：

步骤6.1，确定直线段的参数，利用校正后的估计图像做radon变换取第一个峰值点作为该角度下的投影峰值位置；

步骤6.2，通过修正各角度下每个投影峰值相对估计结果的偏差，对各角度下激光回波波形数据进行配准；

步骤6.3，用配准后的投影数据再次重建目标图像。

采用本发明的激光反射断层成像多角度激光回波配准方法，该方法的原理可靠，易于实现旋转中心的自配准，可以有效减少目标或探测器的抖动，激光脉冲发射和探测器接收的随机抖动等因素的干扰，解决成像图像模糊的问题，提高成像图像的质量。这种利用阈值分割图像和轮廓提取的多角度激光回波配准方法，是一种激光反射断层成像多角度激光回波配准的全新思路。。

附图说明

图1为现有技术的反射投影成像投影图中心的配准流程图；

图2为现有技术中滤波反投影和相位恢复算法重构图像对比；

图3是激光反射断层成像多角度激光回波配准流程图；

图4是本发明实施方式的流程图。

具体实施方式

激光反射断层成像系统在实际工作过程中，由于目标或探测器的抖动，激光脉冲发射和探测器接收的随机抖动等，会在时间维度上出现信号漂移，因此会出现成像图像模糊的问题，影响重构图像的准确性。本申请提供了一种针对非合作空间目标的激光反射断层成像多角度激光回波配准方法，该方法的原理可靠，易于实现例如厘米级空间碎片等尺度微小而外形相对简单的空间目标的多角度激光回波的自配准，可以有效减少目标或探测器的抖动，激光脉冲发射和探测器接收的随机抖动等因素的干扰，解决成像图像模糊的问题，提高成像图像的质量。

以下结合附图对本发明的具体实施方式作出详细说明。

本发明所述的激光反射断层成像多角度激光回波配准方法，其流程图如图3所示，这种方法是利用滤波反投影算法得到目标重建图像，然后对图像进行自动阈值分割处理得到阈值分割图像，根据阈值化处理图像利用线段提取算法给出直线的初始估计，然后由多条估计直线围成的闭合目标轮廓来构建估计图像，通过对估计图像与阈值分割图像的误差函数比较来调制直线参数直到误差函数达到最小值，确定各条直线的参数，提取出其围成封闭图形的轮廓，得到各角度估计投影的峰值位置坐标，最后按照每个投影峰值相对估计结果的偏差进行配准，将各角度投影峰值移动至分段线性估计结果处，从而实现激光反射断层成像多角度激光回波配准。

2、结合上述附图，描述总体技术方案的实现过程

基于阈值分割图像和轮廓提取的激光反射断层成像多角度激光回波自配准方法包括以下六个步骤：

步骤一：利用滤波反投影算法复原目标图像，并对重构图像进行阈值化处理；

步骤二：根据阈值化处理图像给出直线初始估计，由这几条估计直线围成闭合目标轮廓，构建估计图像；

步骤三：计算重建图像与估计图像的误差函数，改变直线参数直到误差函数最小；

步骤四：确定各条直线的参数，提取出其围成封闭图形的轮廓，得到各角度估计投影的峰值位置坐标；

步骤五：按照每个投影峰值相对拟合结果的偏差，进行配准；

步骤六：用配准后的投影数据再次重建目标图像。

实施例1。

图4是实施例1的流程图。

实施例1技术方案详细的实现过程。

①利用R-L(Ram-Lak)函数滤波反投影算法复原目标图像G；

②依据迭代阈值法确定的最佳阈值

③根据阈值化处理图像Gʹ使用LSD(Line Segment Detector)直线段检测算法，给出直线初始估计L

④计算重建图像Gʹ与估计图像G

其中t

⑤确定直线段

最后应说明的是，以上实施方式仅用以说明本发明实施例的技术方案而非限制，尽管参照以上较佳实施方式对本发明实施例进行了详细说明，本领域的普通技术人员应当理解，可以对本发明实施例的技术方案进行修改或等同替换都不应脱离本发明实施例的技术方案的精神和范围。