光电跟踪系统高精度模板匹配跟踪算法

吴浩, 张勇, 李欣, 佀明华, 王伟明

吴浩, 张勇, 李欣, 佀明华, 王伟明. 光电跟踪系统高精度模板匹配跟踪算法[J]. 红外技术, 2022, 44(12): 1301-1308.
引用本文: 吴浩, 张勇, 李欣, 佀明华, 王伟明. 光电跟踪系统高精度模板匹配跟踪算法[J]. 红外技术, 2022, 44(12): 1301-1308.
WU Hao, ZHANG Yong, LI Xin, SI Minghua, WANG Weiming. High-precision Template Matching Tracking Algorithm for Optoelectronic Tracking System[J]. Infrared Technology , 2022, 44(12): 1301-1308.
Citation: WU Hao, ZHANG Yong, LI Xin, SI Minghua, WANG Weiming. High-precision Template Matching Tracking Algorithm for Optoelectronic Tracking System[J]. Infrared Technology , 2022, 44(12): 1301-1308.

光电跟踪系统高精度模板匹配跟踪算法

基金项目: 

国家自然科学基金 51305455

陆军装备部资助课题 ZS2014070140A12002

详细信息
    作者简介:

    吴浩(1993-),男,河北辛集人,硕士研究生,研究方向为光电检测与图像处理

    通讯作者:

    张勇(1974-),男,河北邢台人,博士,高级工程师,研究方向为光电检测与图形图像处理。E-mail:bit10701159@163.com

  • 中图分类号: O439;TP391

High-precision Template Matching Tracking Algorithm for Optoelectronic Tracking System

  • 摘要: 为实现在光电跟踪系统条件下的高精度测量并且满足复杂环境下高精度目标匹配,本文选用去均值归一化互相关匹配算法。为提高匹配速度以及跟踪实时性,利用和表法计算公式中图像求和、平方和图像匹配互相关来简化计算复杂度;采用小波分层金字塔法作为搜索策略,并将模板质心作为参考点进行十字形搜索,引入终止阈值减少误匹配点进一步提高搜索速度。为验证该算法的有效性,实验中将光电跟踪系统放置在二维转台上,调整转台利用该算法跟踪目标靶板。实验结果表明,目标脱靶量控制在3个像素以内,该算法在光电跟踪系统上可实现高精度稳定跟踪。
    Abstract: To achieve high-precision measurements under the operating conditions of optoelectronic tracking systems and satisfy high-precision target matching in complex environments, in this study we adopted the average normalized cross-correlation algorithm. To improve the matching speed and real-time tracking, the computational complexity was simplified by using the sum table method to correlate the sum of images, squares, and the correlation of images. The wavelet pyramid method was used as the search strategy, and the center of the template was used as the reference point for cross-shaped search. A termination threshold was introduced, which reduced the number of mismatched points to increase the search speed. To verify the effectiveness of the algorithm, an optoelectronic tracking system was placed on a two-dimensional turntable in an experiment that used the algorithm to track a target. The experimental results show that the missed target was controlled within 3 pixels. The proposed algorithm can realize high-precision and stable tracking in optoelectronic tracking systems.
  • 基于液晶可调谐滤波器的高光谱成像设备具有成本低、体积小、可靠性高、环境适应性好等特点,且能够在获得目标的空间信息同时,也获得对应的光谱信息,在军事侦察、化学分析、医疗检测等领域有一定的应用研究[1-3]。近年来受到了许多领域学者的关注,张海丹[4]等构建了一套基于LCTF的高光谱火焰成像系统,并基于该系统建立一个碳氢火焰三维温度场和烟黑体积浓度场重建模型,并验证了该模型的合理性,为火场情形科学性判断及火场救灾指挥提供了理论帮助;刘逸飞[5]等利用基于LCTF的凝视型高光谱成像仪实现人脸活体检测,验证了使用高光谱成像技术解决人脸识别系统无法甄别人脸真假的问题,有效减小该系统受到三维面具、视频或人脸照片等欺骗的可能性。朱思祁[6]等利用了基于LCTF的高光谱成像设备结构简单、谱段范围宽、适合探测静止目标的特性,构建了一套高光谱显微系统,实现了生物医学检测应用,该系统能够发现花粉中的特殊结构,并能对胃癌细胞进行精确鉴别,其相关研究成果对高光谱显微成像技术在生物检测应用领域的发展有着积极的意义。

    该类设备利用LCTF实现连续可调谐滤波,在每一波长条件下,使用面阵探测器对目标窗口成像,最后融合为一组高光谱数据立方体,实现高光谱成像。但受成像机理和光学器件的限制,该组合设备成像时间较长,完成典型近红外(900~1700 nm)波段滤波成像往往需要数十秒甚至数分钟时间,大大限制了该类设备在一些有快速、实时检测需求场景下的应用,无法满足动态目标检测和运动检测的需求。

    超分辨光谱重建可实现由低分辨率的光谱数据中得到高分辨率的光谱数据,这使得在使用LCTF进行滤波时,可以在低采样率情况下,超分辨率得到所需的高分辨率光谱信息。在保证超分辨率重建光谱质量的前提下,可以减少LCTF滤波通道数量,降低采样时间,从而有效提高该类设备的光谱采集速度。

    在信号处理领域,信号的采集需要遵循Nyquist采样定理,即采样频率应当大于信号中最高频率的2倍。对于基于LCTF滤波的光谱采集系统来讲,可近似认为,传统采样方法需要至少2倍的最小特征峰半峰宽的滤波步长才能得到满足可用于有效识别物质的光谱信息。2006年,D. Donoho、E. Candès及Tao[7-8]等总结提出了压缩感知理论(Compressed Sensing,CS),该理论认为:如果信号是稀疏的或可压缩的,则可以直接采集压缩后的M个有效测量值,用以实现无失真或较低失真地重建原始信号,即得到满足Nyquist采样定理的N个采样值。根据该理论特点,在该理论框架内,可实现高光谱数据的压缩采集和重建,从而为在保证光谱准确性的前提下提高信息采集速度提供了思路。中国科学院大学汪琪[9]等仿真验证了该思路在珍珠石目标光谱细节识别中的可应用性。

    本文使用CS理论进行光谱超分辨率重建,优化构建了CS理论框架内的要素矩阵,实验验证了光谱重建的可能性,比较了重建结果的识别性能,并找到一组重建效果较好的重建参数,提高了LCTF高光谱数据的采集速度。

    压缩感知理论中信号的采样、压缩编码发生在同一个步骤,利用信号的稀疏性,以远低于Nyquist采样率的速率对信号进行非自适应的测量编码。测量值并非信号本身,而是从高维到低维的投影值,从数学角度看,每个测量值是传统理论下的每个样本信号的组合函数,即一个测量值已经包含了所有样本信号的少量信息。解码过程不是编码的简单逆过程,而是在盲源分离中的求逆思想下,利用信号稀疏分解中已有的重建方法在概率意义上实现信号的精确重建或者一定误差下的近似重建,解码所需测量值的数目远小于传统理论下的样本数。

    根据压缩感知理论原理的数学表示,针对稀疏信号x,给出观测模型yΦ×x,Tao[8]等证明了如果式中观测矩阵Φ满足约束等距理论(Restricted Isometry Property, RIP),则y可由x准确重建。又因重建目标x的维数远高于欠采样信号y,因此重建即为求欠定方程组yΦ×x的最优解。由上可知,压缩感知理论应用的重点在于两个条件的满足:信号x的稀疏性(sparsity)和观测矩阵Φ符合RIP理论要求[10]。另外,信号的稀疏性要求是近似的,若x在原信号域内不是稀疏的,则可以对xΨ的稀疏基上进行稀疏表示为:xΨ×s。将信号投影到稀疏域内,称Ψ为稀疏基矩阵,s为稀疏系数,观测模型则表示为yΦ×Ψ×s。令ΘΦ×Ψ,则通用表达式可写为yΘ×s,称Θ为传感矩阵,该数学表达如图 12所示。

    图  1  观测模型的数学表达示意图
    Figure  1.  Schematic diagram of the observation model's mathematical representation
    图  2  观测模型的数学表达
    Figure  2.  The mathematical representation of the observation model

    至此,依据压缩感知理论框架,就可以用一个与稀疏基矩阵Ψ不相关的观测矩阵Φ将变换所得高维信号投影到一个低维空间上,然后通过求解一个优化问题就可以从这些少量的投影中以高概率重建出原信号。该理论为光谱超分辨率重建提供了理论支撑。

    对于基于LCTF的高光谱成像系统,将目标在波长λ处的反射光能量记为I0(λ),镜头在波长λ处的透过率记为C(λ),LCTF在波长λ处的透过率记为T(λ),相机传感器在波长λ处典型量子效率为Q(λ),则探测器的测量值S(λ)满足:

    $$ S(\lambda ) \propto C(\lambda )T(\lambda )Q(\lambda ){I_0}(\lambda ) $$ (1)

    又因探测器实际探测能量值为其响应范围内所有波长能量值的总和,即对式(1)中λ进行积分,由于LCTF在实际应用中阻断波段处仍无法避免地有极小的透过率,故我们以探测器的响应波长范围作为积分波段。假定探测器的响应范围为λ1λ2,则LCTF调谐波段为i时,探测器接收到的能量值为:

    $$ {S_i} = K\int_{{\lambda _1}}^{{\lambda _2}} C (\lambda )T(\lambda )Q(\lambda ){I_0}(\lambda ){\rm{d}}\lambda $$ (2)

    为便于分析,设定标系数K=1,被积函数C(λ)T(λ)Q(λ)I0(λ)=fi(λ),积分原函数为Fi(λ),将式(2)按牛顿-莱布尼茨公式写为:

    $$ \begin{array}{l} {S_i} = \sum\limits_{k = 1}^n {\int_{{\lambda _1} + (k - 1)\Delta \lambda }^{{\lambda _1} + k\Delta \lambda } {{f_i}} } (\lambda ){\rm{d}}\lambda = \\ \;\;\;\;\;\;\sum\limits_{k = 1}^n {\left( {{F_i}\left( {{\lambda _1} + k\Delta \lambda } \right) - {F_i}\left( {{\lambda _1} + (k - 1)\Delta \lambda } \right)} \right)} \end{array} $$ (3)

    式中:n=(λ2λ1)/∆λ为求和的段数。对式(3)中原函数Fiλ1+kλ处的值进行Taylor展开,展开点取λ1+(k-1)∆λ,则有:

    $$ \begin{array}{*{20}{l}} {{F_i}\left( {{\lambda _1} + k\Delta \lambda } \right) = {F_i}\left[ {{\lambda _1} + (k - 1)\Delta \lambda } \right] + F_i^\prime \left[ {{\lambda _1} + (k - 1)\Delta \lambda } \right]\Delta \lambda + }\\ {F_i^{\prime \prime }\left[ {{\lambda _1} + (k - 1)\Delta \lambda } \right]{{(\Delta \lambda )}^2} + \cdots + }\\ {F_i^{(n)}\left[ {{\lambda _1} + (k - 1)\Delta \lambda } \right]{{(\Delta \lambda )}^n} + {R_n}\left( {{\lambda _1} + k\Delta \lambda } \right)} \end{array} $$ (4)

    式中:Rn(λ1+kλ)为n阶Taylor余项。由原函数的定义Fi'(λ)=fi(λ),并舍去二阶以上展开余项和Taylor余项,则式(4)可写为:

    $$\begin{array}{l} {F_i}({\lambda _1} + k\Delta \lambda ) - {F_i}[{\lambda _1} + (k - 1)\Delta \lambda ] \approx \\ \quad \quad {f_i}[{\lambda _1} + (k - 1)\Delta \lambda ]\Delta \lambda \end{array}$$ (5)

    代入式(3)中,得到离散求和形式:

    $$\begin{array}{l} {S_i} = \sum\limits_{k = 1}^n {{f_i}} [{\lambda _i} + (k - 1)\Delta \lambda ] + o(\Delta \lambda ) = \\ \quad \quad \Delta \lambda \sum\limits_{k = 1}^n {{P_i}[{\lambda _1} + (k - 1)\Delta \lambda ] \times {I_0}[{\lambda _1} + (k - 1)\Delta \lambda ] + o(\Delta \lambda )} \end{array}$$ (6)

    式中:函数Pi(λ)为镜头透过率函数C(λ)、LCTF透过率函数Ti(λ)与探测器响应函数Q(λ)的乘积,可以看作是修正后的透过率函数o(∆λ)为Taylor展开近似的舍入误差。

    从式(6)可以看出,在特定调谐谱段i下得到的测量值与待测目标光谱间的数学关系。其中未知量即为目标反射的光能量I0(λ),其线性系数Pi(λ)是已知的,事先可以通过对镜头、LCTF和探测器的定标测量得出,o(∆λ)忽略不计,等式左边的Si即为探测器在LCTF对应滤波波段处的测量值。

    为求解目标反射光能量I0n个不同波长处的值,在不同LCTF调谐波长处获取传感器的测量值,假设测量次数为mm < < n),则可实现通过m个测量值,计算得到分辨率足够高的n个目标光谱反射值,从而实现目标反射光谱的完整采集。根据这一方法,将式(6)写成矩阵形式:

    $$\begin{array}{l} \left( {\begin{array}{*{20}{c}} {{S_1}}\\ {{S_2}}\\ \vdots \\ {{S_m}} \end{array}} \right) = \Delta \lambda \left( {\begin{array}{*{20}{c}} {{P_1}({\lambda _1})}&{{P_1}({\lambda _1} + \Delta \lambda )}& \cdots &{{P_1}({\lambda _1} + (n - 1)\Delta \lambda )}\\ {{P_2}({\lambda _1})}&{{P_2}({\lambda _1} + \Delta \lambda )}& \cdots &{{P_2}({\lambda _1} + (n - 1)\Delta \lambda )}\\ \vdots & \vdots & \ddots & \vdots \\ {{P_m}({\lambda _1})}&{{P_m}({\lambda _1} + \Delta \lambda )}& \cdots &{{P_m}({\lambda _1} + (n - 1)\Delta \lambda )} \end{array}} \right) \cdot \\ \quad \quad \quad \,\left( {\begin{array}{*{20}{c}} {{I_0}({\lambda _1})}\\ {{I_0}({\lambda _1} + \Delta \lambda )}\\ \vdots \\ {{I_0}({\lambda _1} + (n - 1)\Delta \lambda )} \end{array}} \right) + o(\Delta \lambda ) \end{array}$$

    将上式整理后得:

    $$\mathit{\boldsymbol{S}} = \mathit{\boldsymbol{P}}{\mathit{\boldsymbol{I}}_0} + o(\Delta \lambda )$$ (7)

    式中:Sm×1的测量值向量;I0n×1的目标反射光谱向量;o(∆λ)为误差向量;Pm×n的系数矩阵。

    至此,线性方程组(7)中,测量值m远小于未知光谱构成点数n,即该方程组中未知数个数多于方程个数,故引入压缩感知理论框架,求解该欠定方程组。

    根据CS理论框架可知,使用CS方法实现光谱超分辨率重建,需设计构建适合目标物光谱特征的观测矩阵Φ,稀疏矩阵Ψ,以及寻找一种准确且快速的欠定方程组求解方法。在1.2观测模型的构建理论推导中,可知矩阵P可作为CS理论框架内的观测矩阵Φ。在CS理论应用中常用的稀疏化方法有离散余弦变换(Discrete Cosine Transform,DCT),离散傅里叶变换(Discrete Fourier Transform,DFT),离散小波变换(Discrete Wavelet Transformation,DWT),快速沃尔什变换(Fast Walsh Transform,FWT)以及哈达玛变换(Hadamard Transform,HT)等。使用通用性较强的DCT将目标光谱信号稀疏化。即产生一个n×n的DCT稀疏基矩阵D使得:

    $${\mathit{\boldsymbol{I}}_0} = \mathit{\boldsymbol{D}}x$$ (8)

    因此,式(7)可表示为:

    $$\mathit{\boldsymbol{S}} = \mathit{\boldsymbol{PD}}x + o(\Delta \lambda )$$ (9)

    在稀疏约束下求解x,即求解最优化问题:

    $$\hat x = \arg \min {\left\| x \right\|_0}\quad {\rm{ s}}{\rm{.t}}{\rm{. }}\quad {\left\| {\mathit{\boldsymbol{P}}{\mathit{\boldsymbol{D}}_x} - \mathit{\boldsymbol{S}}} \right\|_2} \le o(\Delta \lambda )$$ (10)

    式中:l0范数为非0值的个数,表示信号的稀疏程度。求解该方程组是一个NP问题,常用贪婪迭代的方法求解。另一种方式,由Candès等证明上述l0范数最小化问题在求解时近似等价于l1范数最小化问题,l1范数最小化问题求解可使用凸优化的方法,文中选用凸优化的方法求解欠定方程组的最优解x。得到方程组的解x即为目标反射光谱经过DCT变换后的投影,根据式(8)得,Dx即为重建后的光谱I0

    搭建基于面阵式近红外相机、LCTF、光学镜头和卤素补光灯的高光谱数据采集平台。平台选用凌云TB-S640-GigE型近红外相机,中达瑞和LCTFSWIR-15-20-S型LCTF以及Schneider Lens Componon-S 4.0/80型号镜头,为弥补实验中红外环境弱导致的信噪比低的问题,设置了两盏50 W的卤素灯对样品进行红外补光,采集平台及环境如图 3所示。

    图  3  高光谱数据采集平台及实验环境
    Figure  3.  Hyperspectral data acquisition platform and experimental environment

    选用化学战剂沙林模拟剂甲基膦酸二甲酯(Dimethyl Methylphosphonate,DMMP)作为检测目标物。用移液枪取600 μL DMMP放入700 μL微量比色皿中,将该比色皿置于标准白板前,根据镜头参数,将样品及背景置于距高光谱成像系统1.5 m处,并设置两盏50 W卤素灯在距离样品10 cm处照射补光。设置LCTF滤波范围为1040~1700 nm,滤波步长为5 nm,相机曝光时间设置为1 ms,使用相机默认参数对LCTF调谐各波长下目标图像进行采集,共得到133幅不同滤波波长下的目标灰度图像,图 4中所示的是滤波波长在1050 nm和1170 nm处,半高宽15 nm条件下,红外相机采集到的背景和样品灰度图像,可以看出在吸收波长附近,目标物有明显的红外吸收。

    图  4  高光谱平台采集特征图像
    Figure  4.  Feature image collected by hyperspectral platform

    将上述133幅灰度图像数据读取并进行预处理,可得到中心指定样品区域内,每个像素点1040~1700 nm的光谱数据。为减少样品区域数据奇点对重建结果的影响,将该区域像素进行了17行4列的栅格划分,共得到68块区域,并以每块区域内所有像素点灰度值的平均值作为该块区域的灰度值,再以行为单位,取每行各块灰度值的中位数作为该行像素点的典型灰度值代表,并以此灰度值进行重建计算。

    由于LCTF透过率的限制,使得反射光能量未能达到探测器满量程,故使用标准白板作为实验背景,视其最大灰度值处反射率为1,以此将探测器的灰度值读数转换为光的反射率,即做了归一化处理。

    信号在某种表示方式下的稀疏性,是压缩感知应用的理论基础。本实验中采用稀疏化能力较强且无需先验信息的DCT离散余弦变换对目标反射光谱进行稀疏化。假设目标重建光谱数据点个数(分辨率)为n,则随机生成一匹配信号数量级的n×n的DCT系数矩阵D图 5为5 nm步长下全采样反射光谱及其在DCT变换后的信号谱图,可以看出原始信号向量左乘稀疏矩阵D后,得到的信号是稀疏的。

    图  5  DCT变换前后信号对比
    Figure  5.  Comparison of signals before and after DCT transform

    本实验中观测矩阵的设计同时也代表着相应采集方法的实现,根据前文理论分析可知,由镜头透过率曲线、LCTF的透过率曲线和探测器量子响应效率的乘积得到的矩阵P可作为CS理论框架中的观测矩阵。在本实验中,根据欠采样时的采样波段和重建目标分辨率n,选取对应波段行中目标分辨率列的向量组成参加计算的观测矩阵P

    根据厂家提供的镜头、LCTF各波段的透过率以及探测器对各波段的量子响应效率得到完备的测量矩阵。本实验使用了凸优化算法中的基于1范数最小的求解算法,使用Cvx凸优化工具箱[11]对方程进行求解。

    使用自主搭建的LCTF高光谱成像系统进行1 nm步长的样品数据采样,得到DMMP目标样品在1040~1700 nm波长范围内的光谱曲线。然后以5 nm间隔从15~100 nm范围内共选择了18种不同的等间距采样步长,基于上述采样方式得到的样本数据进行目标分辨率分别为1.60 nm、3.21 nm和4.81 nm三种符合识别分辨率要求的光谱重建,并对不同参数组合的重建性能进行了评价。找到一组相对较好的重建参数,构建一套具有可行性的LCTF高光谱快速成像方法。

    对于重建结果,从识别性能和误差参数两方面对其进行评价。

    识别性能方面,使用全局相关系数对其进行评价,类比现有使用测量光谱与数据库标准光谱间相关系数划定阈值的识别方法,计算了完整波段范围内重建光谱与传统全采样光谱之间的相关系数,该系数取值范围在0~1之间,相关系数越大代表两光谱越相似,有利于提高判断阈值,从而提高识别的准确性。

    误差参数方面,引用了图像超分辨率重建评价工作中通常使用的峰值信噪比作为评价标准。PSNR是一种常用的客观指标,通常用来衡量有损变换的图像质量。类比图像重建领域的PSNR定义,在光谱重建中,通过反射光谱的最大反射系数L和均方误差(Mean Squared Error,MSE)定义,PSNR与重建光谱的MSE的对数成反比。假设全采样光谱离散数据I和重建光谱数据$\hat I$,两者的光谱采集点个数(即分辨率)都是N,MSE和PSNR(dB)定义如下:

    $${\rm{MSE}} = \frac{1}{N}\sum\limits_{i = 1}^N {{{(I(i) - \hat I(i))}^2}} $$ (11)
    $${\rm{PSNR}} = 10 \times {\rm{lg}}(\frac{{{L^2}}}{{{\rm{MSE}}}})$$ (12)

    式中:对于同一款探测器L是固定值,为探测器的光谱取值范围,本实验选用的近红外相机对光谱表示为14 bit,取值范围为0~16383,则L值为16383,从公式中可看出,L一定,PSNR只与重建光谱的MSE有关,且MSE越小,PSNR越大,说明光谱重建质量越好。

    使用近红外光谱进行目标物识别,一般要求光谱分辨率不大于5 nm,以观察分辨物质种类的特征峰。受到定标仪器分辨率的限制,符合使用需求的先验矩阵分辨率只有1.60 nm,3.21 nm和4.81 nm三种。如图 6所示,将不同采样步长得到的样本,对应不同目标分辨率的重建效果及时间绘制成图。根据图 6中(a)、(b)两图的对比分析,首先找到图(a)中PSNR较高的重建参数点,参考该参数条件下图(b)中的重建时间,在重建时间相差不大的前提下,尽可能选择重建光谱PSNR值高的参数组合作为待选采样方法。据此,选择了图中标注的6组采样参数进行进一步比较和分析。

    图  6  对比选择合适的重建分辨率和采样分辨率
    Figure  6.  Comparison and selection the appropriate reconstruction resolution and sampling resolution

    实际上,需要在保证光谱重建质量的前提下,尽可能提高采样速度,故优先评价重建光谱质量。由图 7所示的重建光谱识别性能评价,在采样步长30 nm,重建分辨率4.81 nm时,重建光谱与全采样光谱相关系数最高,达到0.91,且与其他参数条件下的重建光谱相比,特征峰位相对更准确,故认为该组为较好的重建参数。重建光谱误差参数评价,在该组重建参数下,重建光谱PSNR值第二高,达到99.63 dB。故在优先考虑光谱重建质量的前提下,均衡主客观评价结果,我们认为在采样步长30 nm,重建分辨率4.81 nm参数组合情况下,有着较好的光谱重建质量。

    图  7  不同参数下的原始光谱与重建光谱对比
    Figure  7.  Comparison of original and reconstructed spectrum under different parameters

    根据LCTF厂家提供的数据,其调谐波长时所需液晶稳定时间最少为200 ms,设置相机曝光时间为1 ms,可计算出按4.8 nm步长采集1040~1700 nm波长范围的高光谱数据时,完成一次采集所需时间约为27.54 s。使用该方法重建得到4.81 nm分辨率的高光谱数据时,所需采集步长为30 nm,完成一次采集所需时间约为4.42 s,再需0.56 s完成光谱重建,总耗时约为4.98 s。因此,使该套系统采样时间为传统方式达到同样分辨率所需采样时间的18.08%,速度约为传统采集方法的5.53倍。由此可见,此参数采集条件下,该方法能够在保证光谱质量的前提下实现快速采集。

    本文利用压缩感知原理构建了一种光谱超分辨率重建模型,提出了一种适用于LCTF光谱成像系统的快速采集方法,并通过实验验证了该方法的可行性。鉴于该方法可在保证光谱质量的前提下,实现光谱数据的快速、轻量化采集,为LCTF光谱成像设备在测量动态目标和快速检测等方向的应用提供了可行的应用手段,对LCTF光谱成像设备在更多领域的应用提供了可行性。

  • 图  1   金字塔分层图像的互相关匹配结果

    Figure  1.   Cross-correlation matching results of pyramid images

    图  2   质心法

    Figure  2.   Centroid method

    图  3   直线检测法

    Figure  3.   Line detection method

    图  4   待测图像模板匹配相关度

    Figure  4.   Template matching correlation of image

    图  5   待测图像匹配相关度大于0.8

    Figure  5.   Image matching correlation is greater than 0.8

    图  6   跟踪算法流程

    Figure  6.   Flowchart of tracing algorithm

    图  7   实验装置示意图

    Figure  7.   Experimental setup

    图  8   跟踪结果图像

    Figure  8.   Tracking result of image

    图  9   图像质心坐标序列

    Figure  9.   Sequence of image centroid coordinate

    图  10   非合作目标跟踪图像

    Figure  10.   Tracking image of non-cooperative target

    图  11   非合作目标质心坐标序列

    Figure  11.   Sequence of non-cooperative target centroid coordinate

    表  1   计算复杂度对比

    Table  1   Comparison of computational complexity

    Image size: M×N
    Window size: m×n
    The traditional method The proposed method
    Addition/subtraction Multiplication Addition/subtraction Multiplication
    Construction of sum-tables S1 0 0 M×N 0
    S2 0 0 M×N M×N
    S3 0 0 M×N M×N
    Calculation of NCC ∑∑ M×N×m×n 0 M×N 0
    ∑∑f2 M×N×m×n M×N×m×n M×N 0
    ∑∑f·r M×N×m×n M×N×m×n M×N 0
    Total M×N×m×n M×N×m×n 18×M×N M×N
    下载: 导出CSV
  • [1] 朱泓谕. 模板匹配技术在图像识别中的运用[J]. 电子技术与软件工程, 2021(5): 122-123. https://www.cnki.com.cn/Article/CJFDTOTAL-DZRU202105062.htm

    ZHU Hongyu. Application of template matching technique in image recognition[J]. Electronic Technology & Software Engineering, 2021(5): 122-123. https://www.cnki.com.cn/Article/CJFDTOTAL-DZRU202105062.htm

    [2]

    Lewis J P. Fast normalized cross correlation[C]//Proceeding of Vision Interface, 1995: 120-123.

    [3] 陈翔, 陈鹏. 基于改进模板匹配的目标跟踪算法[J]. 计算机应用, 2011, 31(z2): 127-128. https://www.cnki.com.cn/Article/CJFDTOTAL-JSJY2011S2040.htm

    CHEN Xiang, CHEN Peng. Object tracking algorithm based on improved template matching[J]. Journal of Computer Applications, 2011, 31(z2): 127-128. https://www.cnki.com.cn/Article/CJFDTOTAL-JSJY2011S2040.htm

    [4] 谢维达, 周宇恒, 寇若岚. 一种改进的快速归一化互相关算法[J]. 同济大学学报(自然科学版), 2011, 39(8): 1233-1237. DOI: 10.3969/j.issn.0253-374x.2011.08.025

    XIE Weida, ZHOU Yuheng, KOU Ruolan. An improved fast normalized cross correlation algorithm[J]. Journal of Tongji University (Natural Science), 2011, 39(8): 1233-1237. DOI: 10.3969/j.issn.0253-374x.2011.08.025

    [5] 王斌, 何中市, 伍星, 等. 基于高斯金字塔的图像运动估计算法[J]. 计算机工程与应用, 2015, 51(7): 174-178. DOI: 10.3778/j.issn.1002-8331.1305-0110

    WANG Bin, HE Zhongshi, WU Xing, et al. Image motion estimation algorithm based on Gaussian pyramid[J]. Computer Engineering and Applications, 2015, 51(7): 174-178. DOI: 10.3778/j.issn.1002-8331.1305-0110

    [6] 胡敏, 贺晓佳, 王晓华. 快速区域质心图像匹配算法[J]. 电子测量与仪器学报, 2011, 25(5): 455-462. https://www.cnki.com.cn/Article/CJFDTOTAL-DZIY201105013.htm

    HU Min, HE Xiaojia, WANG Xiaohua. Fast image matching algorithm with area centroid[J]. Journal of Electronic Measurement and Instrumentation, 2011, 25(5): 455-462. https://www.cnki.com.cn/Article/CJFDTOTAL-DZIY201105013.htm

    [7]

    YANG Zhuo. Fast template matching based on normalized cross correlation with centroid bounding[C]//International Conference on Measuring Technology and Mechatronics Automation of IEEE, 2010: 224-227.

    [8] 穆欣侃, 罗海波. 一种对亮度变化鲁棒的相关跟踪方法[J]. 红外与激光工程, 2012, 40(1): 255-260. DOI: 10.3969/j.issn.1007-2276.2012.01.048

    MU Xinkan, LUO Haibo. Correlation tracking approach robust to the variation of image brightness[J]. Infrared and Laser Engineering, 2012, 40(1): 255-260. DOI: 10.3969/j.issn.1007-2276.2012.01.048

    [9] 徐一鸣, 顾菊平, 袁媛, 等. 基于改进归一化积相关算法的目标跟踪方法研究[J]. 南通大学学报(自然科学版), 2013, 12(2): 11-15. DOI: 10.3969/j.issn.1673-2340.2013.02.003

    XU Yiming, GU Juping, YUAN Yuan, et al. Research on target tracking method based on an improved normalized product correlation algorithm[J]. Journal of Nantong University(Natural Science Edition), 2013, 12(2): 11-15. DOI: 10.3969/j.issn.1673-2340.2013.02.003

    [10]

    WEI Shouder, LAI Shanghong. Fast template matching based on normalized cross correlation with adaptive multilevel winner update[J]. IEEE Transactions on Image Processing, 2008, 17(11): 2227-2235. DOI: 10.1109/TIP.2008.2004615

    [11]

    Tsai D M, LIN C T. Fast normalized cross correlation for defect detection[J]. Pattern Recognition Letters, 2003, 24: 2625-2631. DOI: 10.1016/S0167-8655(03)00106-5

    [12]

    Gharavi-Alkhansari M. A fast globally optimal algorithm for template matching using low-resolution pruning[J]. IEEE Transactions on Image Processing, 2001, 10(4): 526-533. DOI: 10.1109/83.913587

    [13] 吴强, 任琳, 张杰. 快速归一化互相关算法及DSP优化实现[J]. 电子测量与仪器学报, 2011, 25(6): 495-499. https://www.cnki.com.cn/Article/CJFDTOTAL-DZIY201106002.htm

    WU Qiang, REN Lin, ZHANG Jie. Fast algorithm of normalized cross correlation and optimized implementation on DSP[J]. Journal of Electronic Measurement and Instrument, 2011, 25(6): 495-499. https://www.cnki.com.cn/Article/CJFDTOTAL-DZIY201106002.htm

    [14] 陈岳军, 孙广玲, 姚恒. 结合小波金字塔的空频域亚像素图像配准[J]. 上海大学学报(自然科学版), 2012, 18(4): 342-348. DOI: 10.3969/j.issn.1007-2861.2012.04.003

    CHEN Yuejun, SUN Guangling, YAO Heng. Hybrid spatial-frequency domain sub-pixel image registration with wavelet pyramid[J]. Journal of Shanghai University(Natural Science Edition), 2012, 18(4): 342-348. DOI: 10.3969/j.issn.1007-2861.2012.04.003

    [15] 吴鹏, 徐洪玲, 宋文龙. 结合小波金字塔的快速NCC图像匹配算法[J]. 哈尔滨工程大学学报, 2017, 38(5): 791-796. https://www.cnki.com.cn/Article/CJFDTOTAL-HEBG201705022.htm

    WU Peng, XU Hongling, SONG Wenlong. A fast NCC image matching algorithm based on wavelet pyramid search strategy[J]. Journal of Harbin Engineering University, 2017, 38(5): 791-796. https://www.cnki.com.cn/Article/CJFDTOTAL-HEBG201705022.htm

    [16] 刘元琳, 宋春凤, 王玲玲. 基于金字塔的渐进分辨率匹配算法研究[J]. 电子制作, 2020(20): 27-29. https://www.cnki.com.cn/Article/CJFDTOTAL-DZZZ202020012.htm

    LIU Yuanlin, SONG Chunfeng, WANG Lingling. Research on progressive resolution matching algorithm based on pyramid[J]. Electronic Production, 2020(20): 27-29. https://www.cnki.com.cn/Article/CJFDTOTAL-DZZZ202020012.htm

    [17] 张万绪, 吴佳丽, 赵丽平, 等. 改进的十字菱形搜索算法INCDS[J]. 西北大学学报(自然科学版), 2011, 41(2): 226-230. https://www.cnki.com.cn/Article/CJFDTOTAL-XBDZ201102009.htm

    ZHANG Wanxu, WU Jiali, ZHAO Liping, et al. Improved new cross-diamond search algorithm INCDS[J]. Journal of Northwest University(Natural Science Edition), 2011, 41(2): 226-230. https://www.cnki.com.cn/Article/CJFDTOTAL-XBDZ201102009.htm

  • 期刊类型引用(0)

    其他类型引用(2)

图(11)  /  表(1)
计量
  • 文章访问数:  164
  • HTML全文浏览量:  37
  • PDF下载量:  58
  • 被引次数: 2
出版历程
  • 收稿日期:  2018-06-19
  • 修回日期:  2022-11-17
  • 刊出日期:  2022-12-19

目录

/

返回文章
返回