[发明专利]基于多线列时差扫描扩展采样的点目标运动速度探测方法

权利要求书

1.一种基于多线列时差扫描扩展采样的点目标运动速度探测方法，其特征在于步骤如下：

(1)构造多线列时差扫描探测装置，该装置包括光学系统、扫描机构和多线列探测器；所述的扫描机构包括摆镜及其驱动转轴；所述的多线列探测器包含N_d个线列探测器，每个线列探测器平行排列，相邻两个线列探测器之间的距离为d_i，其中i＝1,2,…,N_d-1；光学系统和扫描机构一起将视场内场景成像于焦平面，驱动转轴驱动摆镜旋转，线视场内的场景所成像以一定的速率扫过焦平面上前后排列的多线列探测器，扫描角速率为多个线列探测器对视场内同一位置场景先后成像，相邻两个线列成像时间间隔其中d_min为相邻两个线列探测器之间的最小排列距离，v_min可探测的目标最小运动速度，f为光学系统焦距，GSD为线列探测器的地面采样距离；所述N_d取值范围为N_d≥2；

(2)每个线列探测器采用扩展采样方式进行点目标探测，即线列探测器采用N_t个探测阵列组成，像元对应的瞬时视场为IFOV，相邻两个探测阵列平行排列，在垂直扫描方向依次错开1/N_t个像元，并设置探测阵列在扫描方向、在一个采样长度内采样S_t次；所述采样长度为像元对应的瞬时视场；所述的N_t大于等于2；所述S_t取值范围为S_t≥2；

(3)分别将每个线列探测器所采集的N_t组图像数据进行处理，得到亚像元图像；

(4)分别对相邻两个线列探测器处理后的两幅亚像元进行非均匀性校正、亚像素匹配后进行差分计算，完成背景消除；

(5)分别对步骤(4)中相邻两个线列探测器处理后的差分图像进行阈值滤波，并采用邻域约束准则提取其中的正负点对，完成一次扫描过程中运动点目标探测提取；根据相邻线列探测器成像时差和目标运动速度形成邻域约束准则为其中正负区域距离D_i，探测的目标最大运动速度v_max；

(6)对步骤(5)中分别提取的成对的正负标记区域，根据正负标记区域的位置关系，计算目标的运动速度v_i和方向θ_i；

(7)对步骤(6)得到的运动速度v_i和方向θ_i求取均值，得到所探测目标的速度和方向，即

2.根据权利要求1所述的一种基于多线列时差扫描扩展采样的点目标运动速度探测方法，其特征在于：所述步骤(3)中的将每个线列探测器所采集的N_t组图像数据进行处理，得到亚像元图像的方式，具体过程如下：

(3.1)每个线列探测器中的N_t个探测阵列同时成像，得到N_t组图像数据；每个线列探测器按照步骤(2)中设置采样次数对应的采样间隔在扫描方向进行扫描成像，每次成像分别得到N_t组图像数据，每次得到N_t组图像数据后都立即转入步骤(3.2)；

(3.2)分别将每个线列探测器的N_t组图像数据进行对齐拼接处理后形成一帧探测图像；

(3.3)在扫描方向上完成预设的采样次数后，将每个线列探测器对应得到的帧探测图像按照时间进行拼接得到对应的亚像元图像。

3.根据权利要求1所述的一种基于多线列时差扫描扩展采样的点目标运动速度探测方法，其特征在于：在所述步骤(4)中的两幅亚像元图像进行亚像素匹配之前，以在先成像的线列探测器对应的亚像元图像为基准，将后成像的线列探测器对应的亚像元图像在扫描方向向前移动行，L_pi取整数行，线列探测器阵列像元尺寸为a×a。

4.根据权利要求1所述的一种基于多线列时差扫描扩展采样的点目标运动速度探测方法，其特征在于：所述步骤(5)中的滤波采用阈值滤波的方式，具体过程如下：

(4.1)对步骤(4)处理后得到的差分图像中的每个像元进行增强处理，即待处理像元处理后的像素值＝待处理像元的像素值*α+四邻域像元的像素值，α取2～4，处理后得到图像

(4.2)图像取绝对值，记为图像

(4.3)将相邻的数个图像列分为一组，分为N_sub个子图像；子图像的列数L_sub＝H/N_sub,H为图像的总列数；同时对图像采用与图像相同的方式进行分组，图像与的子图像一一对应；

(4.4)从当前子图像的第一行记为当前行开始，设定一个阈值初值T，T取第一个像素的像素值；

(4.5)将当前子图像当前行内的L_sub个像元的像素值依次与阈值初值T相对比，如果像素值超过阈值初值T，则按照式(a)对T提高一定幅度：

T′＝T+0.5^k(V_j-T) (a)

T＝T′

其中k为系数，根据图像信噪比进行调整，V_j代表当前子图像当前行的第j个像元的像素值；

在当前行的最后一个像元处理之后得到的T，即为当前子图像当前行r的最终滤波阈值T_r；

(4.6)利用上述得到的最终阈值T_r对中对应子图像的对应行进行阈值滤波，即对于同一行图像内，将大于滤波阈值的像素标记为1；如果像素值为负，则将小于滤波阈值负值的像素标记为-1；其余的像素均标记为0；

(4.7)将步骤(4.5)最终得到的阈值T_r按公式(b)进行比例降低，做为当前子图像下一行的阈值初值：

T＝(1-0.5^m)T_r (b)

其中m为比例系数；

将当前子图像下一行作为当前行，从步骤(4.5)开始执行，直至当前子图像所有行处理完毕；

(4.8)选择下一子图像作为当前子图像，重复步骤(4.4)～(4.7)，直至完成对图像所有子图像所有行的阈值滤波处理，记为图像

(4.9)遍历的每一列，如果连续2个以上像素被标为1，则保留原标不变，否则将该像素标记为0；如果连续2个以下像素被标记为-1，则保留原标记，否则将该像素标记为0；完成所有像素遍历后，记为图像