基于压缩感知高反光成像技术研究

范剑英; 马明阳; 赵首博

doi:10.11999/JEIT190512

基于压缩感知高反光成像技术研究

doi: 10.11999/JEIT190512

1.
哈尔滨理工大学测控技术与通信工程学院哈尔滨 150080
2.
哈尔滨理工大学测控技术与仪器黑龙江省高校重点实验室哈尔滨 150080

基金项目: 国家自然科学基金(61801148, 61803128)，黑龙江省自然科学基金(QC2016067)

详细信息

作者简介:
范剑英：男，1963年生，教授，硕士生导师，研究方向为光电检测、数字图像与重建

马明阳：男，1993年生，硕士生，研究方向为压缩感知与数字信号处理

赵首博：男，1985年生，副教授，硕士生导师，研究方向为精密光电测量、计算视觉成像

通讯作者:
赵首博　shoubozh@126.com

中图分类号: TN911.73; TN957.52
计量
- 文章访问数: 2595
- HTML全文浏览量: 1023
- PDF下载量: 63
- 被引次数: 0
出版历程
- 收稿日期: 2019-07-09
- 修回日期: 2020-01-17
- 网络出版日期: 2020-02-17
- 刊出日期: 2020-06-04

Research on High Reflective Imaging Technology Based on Compressed Sensing

1.
School of Measurement and Control Technology and Communication Engineering, Harbin University of Science and Technology, Harbin 150080, China
2.
Measurement and Control Technology and Instrument Key Laboratory of Heilongjiang Province, Harbin University of Science and Technology, Harbin 150080, China

Funds: The National Natural Science Foundation of China (61801148, 61803128), The Scientific Research Foundation of Heilongjiang Province (QC2016067)

摘要

摘要:
高反光物体成像时反射的光强容易超出传感器接收光强的最大量化值，使得采集图像部分区域图像失真，严重影响信息传递。为了改善高反光成像饱和区域中数据丢失的状况，该文结合压缩感知这一新的采样理论提出基于压缩感知高反光成像方法，利用特定测量矩阵对目标图像进行线性采样，将CCD图像传感器的单个光强采样值与测量矩阵中的分布数据对应结合，对整合后的数据用算法进行恢复重建实现被测目标在高光环境中成像。以峰值信噪比和灰度直方图作为客观评定标准。实验表明，该成像方法鲁棒性较强、可行性较高，直方图检测饱和像素占比为0%，峰值信噪比为58.37 dB实现了在高光环境下不含饱和光成像，为压缩感知在成像应用中提供了新的方向。
- 高反光成像 /
- 压缩感知 /
- 数据整合 /
- 直方图
Abstract:
When imaging a highly reflective object, the light intensity reflected easily exceeds the maximum quantized value of the light intensity received by the sensor, which causes image distortion of the captured image in the saturated region of light intensity and seriously affects the quality of information transmission. In order to improve the data loss in the high-reflection imaging saturation region, a compression-sensing of high-reflection imaging method based on the new sampling theory of compressed sensing is proposed. A specific measurement matrix is used to conduct linear sampling of the target image, and the single light intensity sampling value of the CCD image sensor is combined with the distribution data in the measurement matrix, and the integrated data is restored and reconstructed with the algorithm to achieve the imaging of the measured target in the high-light environment. The peak signal to noise ratio and gray histogram are used as objective evaluation criteria. Experiments show that this imaging method is robust and feasible, with the proportion of saturated pixels in histogram detection 0% and the peak signal to noise ratio 58.37 dB, realizing the imaging without saturated light in the high-light environment, providing a new direction for the application of compressed sensing in imaging.
- High reflected light imaging /
- Compressed sensing /
- Data integration /
- Histogram

HTML全文

1. 引言

随着浅层陆地资源的日益匮乏与开发难度的增大，世界各国对海洋资源的开发与海域权益的争夺日趋激烈。面对部分邻国蓄意挑起的海洋争端及来自境外军事侦查、非法入侵和资源掠夺的敌方潜艇及无人潜航器等对我国的临海主权，尤其是对港口、码头及军事设施等的安全构成严重威胁^[1,2]。磁探测技术是一种被动探测技术，因其能够有效弥补声学近海、浅海探测的不足受到国内外的广泛关注^[3]。

磁性目标的特性识别主要包括优化定位与实时滤波跟踪两大类。其中，优化定位算法通常依赖大量观测数据，比如传感器网络、磁梯度系统或航空测线数据等，通过最优化的方法获取目标位置、磁矩等参量。优化定位算法研究的典型代表有，以色列的Sheinker等人^[4,5]在2005～2008年提出依次基于遗传算法、模拟退火算法及种群增量学习算法实现磁偶极子目标的定位。2010年Yang等人^[6]将粒子群与(Levenberg-Marquardt, LM)优化算法结合实现了3个磁偶极子目标的定位与跟踪，2018年Gao等人^[7]将LM算法用于运动磁性目标的定位。此外，利用磁梯度张量系统对UXO目标的定位与识别UXO同样吸引了大量研究人员的关注^[8-10]，其通常将搭载在移动平台上，如无人机等。在目标的实时滤波跟踪方面，通常利用静止平台的磁传感器阵列测量目标磁异常，实现对目标体的实时跟踪。滤波跟踪算法研究的典型代表有，Birsan等人^[11,12]在2003～2006年利用扩展/无迹卡尔曼滤波方法和无迹粒子滤波算法等迭代贝叶斯算法对磁偶极子目标进行实时跟踪；2008年，Kozick等人^[13]提出磁传感器阵列目标跟踪的矢量与标量模型，并进行仿真实验研究；2012年日本学者Alhmiedat等人^[14]提出将无线磁传感器网络用于实现军事目标的识别与跟踪应用。国内研究学者在近几年针对磁性目标的跟踪已开展大量工作，主要的研究成果有：2010年，张朝阳等人^[15]对磁性目标跟踪中磁偶极子模型适用性进行研究；2012年，于振涛等人^[16]利用海底磁力仪阵列对水中磁目标定位与跟踪方法进行研究；2013年，吴志东等人^[17,18]提出用于实现磁偶极子目标的非线性跟踪算法；2014年，张宏欣等人^[19]提出改进的粗糙化粒子滤波的算法改善磁偶极子目标跟踪中的粒子贫化；2015年，周穗华等人^[20]提出混合卡尔曼滤波算法解决经典粒子滤波在磁偶极子跟踪中高维模型估计时维数灾问题；2017年吴垣甫等人^[21]提出一种递推更新的卡尔曼滤波方法对磁偶极子目标进行跟踪，有效抑制磁偶极子跟踪中初始误差过大导致的滤波发散问题；张宏欣等人^[22]提出渐近贝叶斯估计算法和多模型自适应滤波算法对磁偶极子目标进行跟踪，有效抑制目标跟踪中初始误差导致的性能下降和发散等问题；在2018年张宏欣^[23]提出多模型自适应滤波算法对磁偶极子目标进行跟踪，算法可在目标源及位置等先验信息完全未知条件下准确估计出真实目标状态，提升真实目标跟踪中的适应力。

分析当前国内外研究现状，磁目标跟踪的主要研究对象在于非线性跟踪模型的优化与不同体制跟踪系统的研究，其中很少考虑到复杂跟踪环境对目标跟踪性能的影响，比如海水中目标在运动过程中存在的涡流磁场、尾流磁场^[18]，以及目标跟踪过程中存在地磁背景干扰等，都会对磁目标的跟踪性能产生影响。本文主要针对磁目标跟踪过程中存在的地磁背景干扰的问题开展相关的研究，并在传统磁跟踪模型的基础上提出利用差量磁场对磁偶极子进行跟踪的模型，并考虑到差量磁场模型的高度非线性度，提出基于蒙特卡洛卡尔曼滤波(Monte Carlo Kalman Filter, MCKF)跟踪方法提升跟踪性能。

2. 磁偶极子运动模型

铁磁性目标周围产生磁场叠加原背景磁场导致周围的磁场发生畸变，称为磁异常。当目标探测距离较近时，磁性目标可等效为若干个偶极子源叠加；当探测距离为目标尺寸3倍以上时，可将目标视为磁偶极子，其矢量磁场如式(1)，式中 ${\mu _0}$ 为自由空间磁导率， ${{m}}$ 为目标磁矩矢量， ${{r}}$ 为距离矢量，其幅度为 $r$ 。

$B\left( {{{m}},{{r}}} \right) = \frac{{{\mu _0}}}{{4\pi }}\frac{{3\left( {{{r}} \cdot {{m}}} \right){{r}} - {r^2}{{m}}}}{{{r^5}}}$

(1)

确立直角坐标系，式(1)可以写为如式(2)形式。式中 $\left( {{x_0},{y_0},{z_0}} \right)$ 代表了在3维坐标系中矢量磁力仪的坐标， $\left( {x,y,z} \right)$ 代表坐标系中磁偶极子的坐标； $\left( {{B_x},{B_y},{B_z}} \right)$ 代表了磁场矢量在各坐标轴上的投影分量，而 $\left( {{m_x},{m_y},{m_z}} \right)$ 代表了磁矩矢量在3维坐标系统的投影分量。

$\left[ {\begin{array}{*{20}{c}} {{B_x}} \\ {{B_y}} \\ {{B_z}} \end{array}} \right] = \frac{{{\mu _0}}}{{4\pi {r^5}}} \cdot \left[ {\begin{array}{*{20}{c}} {3{{\left( {x - {x_0}} \right)}^2} - {r^2}}&{3\left( {x - {x_0}} \right)\left( {y - {y_0}} \right)}&{3\left( {x - {x_0}} \right)\left( {z - {z_0}} \right)} \\ {3\left( {y - {y_0}} \right)\left( {x - {x_0}} \right)}&{3{{\left( {y - {y_0}} \right)}^2} - {r^2}}&{3\left( {y - {y_0}} \right)\left( {z - {z_0}} \right)} \\ {3\left( {z - {z_0}} \right)\left( {x - {x_0}} \right)}&{3\left( {z - {z_0}} \right)\left( {y - {y_0}} \right)}&{3{{\left( {z - {z_0}} \right)}^2} - {r^2}} \end{array}} \right] \cdot \left[ {\begin{array}{*{20}{c}} {{m_x}} \\ {{m_y}} \\ {{m_z}} \end{array}} \right]$

(2)

由式(2)可见，任意磁偶极子目标产生的磁场仅依赖于目标与传感器的相对位置及磁偶极子的磁矩。因此，通过对目标位置与磁矩参量的估计，即可完成对运动目标的实时跟踪。

磁偶极子目标跟踪过程中状态变量由两部分组成：在t时刻，描述目标状态的位置参数 $\left[ {x\left( t \right),y\left( t \right),z\left( t \right)} \right]$ 和速度参数 $\left[ {{v_x}\left( t \right),{v_y}\left( t \right),{v_z}\left( t \right)} \right]$ 以及目标的磁矩参数 $\left[ {{m_x}\left( t \right),{m_y}\left( t \right),{m_z}\left( t \right)} \right]$ 。因此，目标状态向量 ${{{x}}_{{t}}}$ 如式(3)。

$\begin{split} {{{x}}_{{t}}} =\,& \left[ x\left( t \right),y\left( t \right),z\left( t \right),{v_x}\left( t \right),{v_y}\left( t \right),{v_z}\left( t \right),\right. \\ & \left.{m_x}\left( t \right),{m_y}\left( t \right),{m_z}\left( t \right) \right] \end{split}$

(3)

考虑到水面舰船或航行器等磁性目标运动状态相对稳定，目标运动状态可以用恒速模型和恒磁矩模型描述，因此目标运动状态方程如式(4)。

$\left. \begin{array}{l} x\left( t \right) = x\left( {t - \Delta t} \right) + \Delta t \cdot {v_x}\left( {t - \Delta t} \right) \\ y\left( t \right) = y\left( {t - \Delta t} \right) + \Delta t \cdot {v_y}\left( {t - \Delta t} \right) \\ z\left( t \right) = z\left( {t - \Delta t} \right) + \Delta t \cdot {v_z}\left( {t - \Delta t} \right) \\ {v_x}\left( t \right) = {v_x}\left( {t - \Delta t} \right) \\ {v_y}\left( t \right) = {v_y}\left( {t - \Delta t} \right) \\ {v_z}\left( t \right) = {v_z}\left( {t - \Delta t} \right) \\ {m_x}\left( t \right) = {m_x}\left( {t - \Delta t} \right) \\ {m_y}\left( t \right) = {m_y}\left( {t - \Delta t} \right) \\ {m_z}\left( t \right) = {m_z}\left( {t - \Delta t} \right) \\ \end{array} \right\}$

(4)

根据目标的位置参量与磁矩参量，通过式(2)，可以获得不同位置传感器的磁异常，当使用3个矢量传感器阵列时，各传感器测量场表示为 $\left[ b_x^i\left( t \right),\right.$ $\left.b_y^i\left( t \right),b_z^i\left( t \right),i = 1,2,3 \right]$ 。因此，t时刻的观测场 ${{{z}}_{{t}}}$ 可表示如式(5)。

$\begin{split} {{{z}}_{{t}}} =\,& \left[ b_x^1\left( t \right),b_y^1\left( t \right),b_z^1\left( t \right),b_x^2\left( t \right),b_y^2\left( t \right),b_z^2\left( t \right),\right.\\ & \left.b_x^3\left( t \right),b_y^3\left( t \right),b_z^3\left( t \right) \right] \end{split}$

(5)

综上，偶极子目标跟踪过程中的状态方程式(4)为线性的；而量测方程式(2)、式(5)为非线性方程。

3. 基于差量磁场的MCKF跟踪算法

3.1 差量磁场信号

传感器测量磁场中同时包含目标磁异常与背景地磁场。通常，采用远参考传感器来消除背景地磁场的影响，但基于参考传感器的方法存在两个明显缺陷：其一，参考传感器与传感器阵列的距离很难把握，距离过近则容易受目标磁异常影响；距离过远，则由于相关性的下降致使背景场消除不彻底。其二，参考传感器设置处需要可靠的背景环境，要保证工作中不受外界磁干扰影响。本文针对上述缺陷，提出传感器阵列测量磁场差量实现磁偶极子目标跟踪，在消除背景地磁场干扰的同时，又避免参考传感器的设置难题。

通过仿真示例，阐述磁偶极子差量磁场与测量场及梯度磁场的特性。仿真磁偶极子运动轨迹如图1所示，其轨迹中类似波浪式的运动用于模拟真实目标运动过程中的随机性，“*”代表测量磁传感器，仿真磁偶极子磁矩为 $\left[ {10,\;20,\;30} \right]{\rm{A}} \cdot {{\rm{m}}^2}$ 。

图 1 仿真磁偶极子运动轨迹与传感器分布位置

下载: 全尺寸图片幻灯片

传感器测量磁场如图2所示，其中图2(a)代表传感器1测量的磁场3分量；图2(b)代表传感器1处的磁梯度场中3个典型量，依次是x分量场在x方向梯度、y分量场在y方向梯度及z分量场在z方向的梯度；图2(c)代表传感器1, 2的差量场。

图 2 传感器测量结果

下载: 全尺寸图片幻灯片

从图2可知，差量场与传感器原始测量磁场在磁异常幅度值上相当，而梯度磁场最小。因此，与磁梯度场跟踪相比，采用差量磁场进行目标跟踪时可以获得更远的探测距离。

3.2 基于MCKF的跟踪算法

常用的非线性滤波跟踪算法有扩展卡尔曼滤波(Extended Kalman Filter, EKF)和无迹卡尔曼滤波(Unscented Kalman Filter, UKF)等。EKF算法利用泰勒级数展开，对非线性的系统方程作1阶线性化近似处理，仅适用于弱非线性系统。UKF以确定性的方式选择一组近似高斯分布的样本点，通过非线性方程变换，获得均值和方差的估计，是最常用的非线性卡尔曼滤波算法。本文提出的蒙特卡洛卡尔曼滤波方法，与UKF原理相似，区别在于通过蒙特卡洛的方法生成随机sigma样本点，通过非线性方程变换，获取非线性系统均值与方差估计，完成卡尔曼滤波增益计算与状态的更新，算法流程如图3所示。

图 3 基于MCKF跟踪算法的流程图

下载: 全尺寸图片幻灯片

基于图3的算法流程，目标状态变量 ${{x}}$ 及其状态方程与观测方程分别用 $f\left( {{x}}\right)$ 和 $h\left( {{x}} \right)$ 表示，则磁偶极子的蒙特卡洛卡尔曼滤波跟踪算法的步骤可归结如下：

(1) 滤波器初始化：定义初始状态变量为 $9 \times 1$ 维列向量，记为x₀；初始化协方差矩阵为 $9 \times 9$ 矩阵，标记为P₀^xx。

(2) 状态变量预测：根据蒙特卡洛规则随机生成N个样本点，记 ${{{x}}^{(i)}}_{n - 1|n - 1} = {{{x}}_{n - 1|n - 1}} +$ ${{{D}}_{n - 1|n - 1}}c(i),i =1,2, ··· ,N$ ，样本点权重记为 ${w_i} = 1/N,i = 1,2, ··· ,N$ 。根据随机样本点获取状态量均值预测记为式(6)，其协方差矩阵预测如式(7)，式中 $c\left( i \right)$ ～ $N\left( {0,{{I}}} \right),i = 1,2, ··· ,N$ ; ${{{D}}_{n - 1|n - 1}} =$ ${\left[ {{{P}}_{n - 1|n - 1}^{{xx}}} \right]^{ - 1}},n = 2,3,4, ···$ ; Q为变量预测协方差矩阵。

${{{x}}_{n|n - 1}} = \mathop \sum \limits_{i = 1}^N {w_i}f\left( {{{x}}_{n - 1|n - 1}^{\left( i \right)}} \right) \hspace{58pt}$

(6)

$\begin{split} {{P}}_{n|n - 1}^{{xx}} =\,& \mathop \sum \limits_{i = 1}^N {w_i}\left[ {f\left( {{{x}}_{n - 1|n - 1}^{\left( i \right)}} \right) - {{{x}}_{n|n - 1}}} \right]\\ &\cdot {\left[ {f\left( {{{X}}_{n - 1|n - 1}^{\left( i \right)}} \right) - {{{x}}_{n|n - 1}}} \right]^{\rm{T}}} + {{Q}} \end{split}$

(7)

(3) 观测值及其协方差矩阵估计：根据蒙特卡洛规则生成样本点 ${{x}}_{n|n - 1}^{\left( i \right)} = {{{x}}_{n|n - 1}} + {{{D}}_{n|n - 1}}c\left( i \right),$ $i = 1,2, ··· ,N$ ，通过量测方程估计观测均值如式(8)。利用输出样本点及输入样本点获取观测值协方差及观测量与状态变量互协方差矩阵如式(9)，式中 ${{{D}}_{n|n - 1}} = {\left[ {{{P}}_{n|n - 1}^{{xx}}} \right]^{ - 1}},n = 2,3,4, ···$ ，R为观测噪声协方差矩阵。

${{{z}}_{n|n - 1}} = \mathop \sum \limits_{i = 1}^N {w_i}h\left( {{{X}}_{n|n - 1}^{\left( i \right)}} \right)\hspace{73pt}$

(8)

$\left. \begin{split} {{P}}_{n|n - 1}^{{zz}} =\,& \mathop \sum \limits_{i = 1}^N {w_i}\left[ {h\left( {{{X}}_{n - 1|n - 1}^{\left( i \right)}} \right) - {{{z}}_{n|n - 1}}} \right] \\ &\cdot {\left[ {h\left( {{{X}}_{n - 1|n - 1}^{\left( i \right)}} \right) - {{{z}}_{n|n - 1}}} \right]^{\rm{T}}} + {{R}} \\ {{P}}_{n|n - 1}^{{xz}} =\, & \mathop \sum \limits_{i = 1}^N {w_i}\left[ {f\left( {{{X}}_{n - 1|n - 1}^{\left( i \right)}} \right) - {{{x}}_{n|n - 1}}} \right]\\ & \cdot {\left[ {h\left( {{{X}}_{n - 1|n - 1}^{\left( i \right)}} \right) - {{{z}}_{n|n - 1}}} \right]^{\rm{T}}} \end{split} \right\}$

(9)

(4) 增益计算与状态变量更新：利用第(3)步得到的协方差矩阵求得滤波的增益估计，如式(10)所示。利用卡尔曼增益与实时观测值，更新状态变量及其协方差矩阵，如式(11)和式(12)。

${K_n} = {{P}}_{n|n - 1}^{{{xz}}}{\left( {{{P}}_{n|n - 1}^{{{zz}}}} \right)^{ - 1}}\hspace{40pt}$

(10)

${{{x}}_{n|n}} = {{{x}}_{n|n - 1}} + {K_n}\left( {{{{z}}_n} - {{{z}}_{n|n - 1}}} \right)$

(11)

${{P}}_{n|n}^{{{xx}}} = {{P}}_{n|n - 1}^{{{xx}}} - {K_n}{{P}}_{n|n - 1}^{{{zz}}}K_n^{\rm{T}} \hspace{20pt}$

(12)

(5) 存储与迭代：将该时刻状态变量及其协方差存储至跟踪文件，返回第(2)步进行下一时刻迭代。

4. 实验与算法验证

4.1 仿真实验

仿真实验如3.1节中描述，磁偶极子目标运动轨迹如图1。用于实现目标跟踪的传感器阵列由3个矢量磁力仪构成，坐标分别为(–2 m, 0, 0), (0, 0, 0), (2 m, 0, 0)，依次标记为传感器1、传感器2、传感器3。通过传感器1、传感器2和传感器2、传感器3分别获得差量磁场，并在原始差量场上叠加高斯白噪声模拟工作的噪声。利用观测磁场，分别采用EKF, UKF与本文提出的基于MCKF跟踪算法，实现目标实时跟踪。为评估不同算法的跟踪性能，采用100次随机实验，将跟踪误差的统计平均值作为评判指标。不同仿真算法误差结果如图4所示，其中图4(a)为x轴投影轨迹的误差，图4(b)为y轴投影轨迹的误差，图4(c)为z轴投影轨迹的误差。

图 4 仿真磁偶极子跟踪结果对比

下载: 全尺寸图片幻灯片

从图4可以看出，EKF与UKF具有类似的跟踪效果，MCKF的初期跟踪效果较差，而后期跟踪效果明显优于EKF/UKF的跟踪结果。定性来讲，MCKF在收敛速度上较EKF/UKF更慢，但稳定后在跟踪效果上跟踪精度最高。为定量分析不同算法跟踪性能，计算3种跟踪方法跟踪误差结果如表1。

表 1 不同跟踪算法各方向投影轨迹跟踪误差(m)

	跟踪算法	时间区间
	跟踪算法	1～40点	41～80点	81～120点	120～160点	160～200点
x方向	EKF	0.0651	0.0127	0.0023	0.0057	0.0409
	UKF	0.0689	0.0127	0.0024	0.0057	0.0410
	MCKF	0.2143	0.0153	0.0024	0.0039	0.0218
y方向	EKF	0.0491	0.0113	0.0059	0.0100	0.0319
	UKF	0.0512	0.0112	0.0061	0.0096	0.0318
	MCKF	0.1237	0.0115	0.0042	0.0061	0.0158
z方向	EKF	0.0430	0.0089	0.0021	0.0042	0.0163
	UKF	0.0436	0.0088	0.0022	0.0041	0.0164
	MCKF	0.0456	0.0086	0.0022	0.0043	0.0198

下载: 导出CSV

| 显示表格

表1中列举不同时间区间，x, y, z等3个方向上投影轨迹跟踪误差。从表1中看出MCKF算法与EKF/UKF算法相比，在1～40点内跟踪误差较大，随后在40～80点内跟踪误差相当，而在稳定的80～200点内，MCKF的跟踪精度明显优于其余两种跟踪算法；表1结果表明，在z方向3种算法跟踪误差相当。

依据MCKF算法原理，蒙特卡洛样本点数会影响算法跟踪性能。为评估不同蒙特卡洛样本点数的影响，依次采用50, 100, 200, 400样本点进行仿真，结果如表2。从表2中得出随粒子数由50增至200，算法的跟踪误差逐步减小，说明增加粒子数可有效提升算法跟踪性能；而当粒子数由200增至400时，跟踪性能的提升并不明显。由于粒子数多少决定计算量，因此目标跟踪过程中应在满足性能时尽量减少粒子数。

表 2 不同蒙特卡洛样本点数的算法各方向投影轨迹跟踪误差(m)

表2(a) x方向
蒙特卡洛样本点数	时间区间
蒙特卡洛样本点数	1～40	41～80	81～120	120～160	160～200
50	0.2817	0.0242	0.0035	0.0052	0.0264
100	0.2481	0.0194	0.0030	0.0046	0.0231
200	0.2143	0.0153	0.0024	0.0039	0.0218
400	0.2316	0.0159	0.0023	0.0039	0.0212
表2(b) y方向
蒙特卡洛样本点数	时间区间
蒙特卡洛样本点数	1～40	41～80	81～120	120～160	160～200
50	0.1658	0.0168	0.0049	0.0069	0.0192
100	0.1442	0.0141	0.0047	0.0063	0.0166
200	0.1237	0.0115	0.0042	0.0061	0.0158
400	0.1346	0.0118	0.0043	0.0059	0.0151
表2(c) z方向
蒙特卡洛样本点数	时间区间
蒙特卡洛样本点数	1～40	41～80	81～120	120～160	160～200
50	0.0569	0.0124	0.0031	0.0053	0.0242
100	0.0496	0.0105	0.0026	0.0046	0.0202
200	0.0456	0.0086	0.0022	0.0043	0.0198
400	0.0466	0.0090	0.0021	0.0042	0.0195

下载: 导出CSV

| 显示表格

4.2 地磁背景跟踪实验

地磁背景跟踪实验在张家口张北县某地草原进行，选用Mag-03磁通门传感器，传感器分布与磁性目标运动轨迹示意如图5，其中模拟跟踪磁性目标体为通电螺旋圈如图6，磁矩约为60 A·m²。根据图5(b)中运动轨迹，磁性目标在运动过程中的拐点位置由于磁性源的姿态发生变化，此时目标分量磁矩会产生突变。

图 5 真实目标跟踪实验

下载: 全尺寸图片幻灯片

图 6 跟踪实验模拟磁性目标体

下载: 全尺寸图片幻灯片

传感器1、传感器2、传感器3测量的磁异常如图7(a)、图7(b)、图7(c)所示，传感器1、传感器2的差量磁场如图7(d)，传感器2、传感器3的差量磁场如图7(e)。差量磁异常信号与原始磁异常信号在幅度上相近，甚至大于原始测量磁异常。利用本节中提出的基于MCKF跟踪方法与传统的基于EKF/UKF的滤波跟踪方法对目标进行跟踪，目标位置跟踪结果如图8所示，其中图8(a)、图8(b)、图8(c)分别代表x, y, z的跟踪结果，黑色实线为参考轨迹；目标磁矩的跟踪结果如图9所示，其中图9(a)、图9(b)、图9(c)分别代表磁矩x, y, z分量磁矩的跟踪结果，黑色实线为参考磁矩。

图 7 阵列传感器测量磁场与差量磁场

下载: 全尺寸图片幻灯片

图 8 模拟磁性目标的位置跟踪结果

下载: 全尺寸图片幻灯片

图 9 模拟磁性目标的磁矩跟踪结果

下载: 全尺寸图片幻灯片

由图8和图9的跟踪结果可以看出：(1)在跟踪初始阶段及拐点位置，EKF与UKF算法跟踪响应速度要优于MCKF跟踪方法，与仿真结果一致；(2)当目标运动状态稳定后(图中3500点后)跟踪过程中，MCKF的跟踪性能明显优于其余两种算法。其中，EKF与UKF算法在4300点后，由于信噪比下降算法的跟踪轨迹发生明显偏离，而MCKF直到4800点后轨迹才逐渐发生偏移。综上所述，真实目标的跟踪结果与仿真结果表现一致，说明本文提出的算法具有较好的抗地磁干扰性能，在低信噪比时跟踪性能优于传统滤波方法。

5. 结束语

本文分析了地磁场背景下磁偶极子目标跟踪过程中存在背景干扰问题，提出利用传感器差量磁场信号实现目标跟踪的模型。新模型的优势有：(1)传感器阵列测量磁场差量与原始测量场相比磁场强度相似，因此可以保证足够远的跟踪距离；(2)测量磁场差量不易受到背景地磁场的干扰。同时，考虑差量磁场模型的高度非线性，本文提出基于MCKF的滤波跟踪算法，从而有效解决减小跟踪模型的非线性误差。通过仿真实验，结果表明基于新模型的跟踪算法较经典的UKF、EKF算法具有更高的跟踪精度；地磁背景跟踪实验，同样证实算法的有效性。

图 1 压缩感知框架图

下载: 全尺寸图片幻灯片

图 2 不同环境下成像状态

下载: 全尺寸图片幻灯片

图 3 图像分块

下载: 全尺寸图片幻灯片

图 4 CCD所采集含高反光图像

下载: 全尺寸图片幻灯片

图 5 去除成像中高亮区域效果

下载: 全尺寸图片幻灯片

图 6 压缩感知不同恢复算法去除饱和光成像

下载: 全尺寸图片幻灯片

图 7 压缩感知去饱和光成像光路图

下载: 全尺寸图片幻灯片

图 8 高亮光环境下采集的被测目标信息

下载: 全尺寸图片幻灯片

图 9 原始图像和压缩感知高反光成像直方图

下载: 全尺寸图片幻灯片

表 1 不同采样率下两种恢复算法的MSE值与PSNR值

CS采样率 OMP算法 SAMP算法
MSE PSNR MSE PSNR

0.30 0.0158 66.1481 0.0166 65.9403
0.35 0.0150 66.3739 0.0154 66.2512
0.40 0.0143 66.5625 0.0142 66.6061
0.45 0.0139 66.7149 0.0135 66.8149
0.50 0.0136 66.7928 0.0133 66.9078

下载: 导出CSV

参考文献(23)

赵首博, 曲兴华, 冯维, 等. 成像前光学调制系统的眩光测量[J]. 光电工程, 2016, 43(1): 13–17. doi: 10.3969/j.issn.1003-501X.2016.01.003

ZHAO Shoubo, QU Xinghua, FENG Wei, et al. A measurement method for glare based on the Optical modulation system before image formation[J]. Opto-Electronic Engineering, 2016, 43(1): 13–17. doi: 10.3969/j.issn.1003-501X.2016.01.003

范剑英, 刘力源, 赵首博. 电机铜排表面毛刺缺陷检测技术研究[J]. 仪器仪表学报, 2019, 40(3): 14–22. doi: 10.19650/j.cnki.cjsi.J1804559

FAN Jianying, LIU Liyuan, and ZHAO Shoubo. Research on detection technology of burr defects in motor copper[J]. Chinese Journal of Scientific Instrument, 2019, 40(3): 14–22. doi: 10.19650/j.cnki.cjsi.J1804559

ZHAO Shoubo, MA Mingyang, and GUO Cong. Accurate Pixel-to-Pixel alignment method with Six-Axis adjustment for computational photography[J]. IEEE Photonics Journal, 2018, 10(3): 6802416. doi: 10.1109/jphot.2018.2839093

JIANG Hongzhi, ZHAO Huijie, and LI Xudong. High dynamic range fringe acquisition: A novel 3-D scanning technique for high-reflective surfaces[J]. Optics and Lasers in Engineering, 2012, 50(10): 1484–1493. doi: 10.1016/j.optlaseng.2011.11.021

ZHAO Shoubo, ZHANG Fumin, QU Xinghua, et al. Removal of parasitic image due to metal specularity based on digital micromirror device camera[J]. Optical Engineering, 2014, 53(6): 063105. doi: 10.1117/1.oe.53.6.063105

YAN Qingsen, ZHU Yu, and ZHANG Yanning. Robust artifact-free high dynamic range imaging of dynamic scenes[J]. Multimedia Tools and Applications, 2019, 78(9): 11487–11505. doi: 10.1007/s11042-018-6625-x

KALANTARI N K, SHECHTMAN E, BARNES C, et al. Patch-based high dynamic range video[J]. ACM Transactions on Graphics, 2013, 32(6): No. 202. doi: 10.1145/2508363.2508402

ZHAO Shoubo, LIU Liyuan, and MA Mingyang. Adaptive high-dynamic range three-dimensional shape measurement using DMD camera[J]. IEEE Access, 2019, 7: 67934–67943. doi: 10.1109/access.2019.2918843

DONOHO D L. Compressed sensing[J]. IEEE Transactions on Information Theory, 2006, 52(4): 1289–1306. doi: 10.1109/TIT.2006.871582

ELDAR Y C, KUTYNIOK G. Compressed Sensing: Theory and Applications[M]. Cambridge: Cambridge University Press, 2012: 1289–1306. doi: 10.1017/CBO9780511794308.

CANDES E J, ROMBERG J, and TAO T. Robust uncertainty principles: Exact signal reconstruction from highly incomplete frequency information[J]. IEEE Transactions on Information Theory, 2006, 52(2): 489–509. doi: 10.1109/TIT.2005.862083

RISTOBAL-HUERTA A, POOT D H J, VOGEL M W, et al. Compressed sensing 3D-GRASE for faster high-resolution MRI[J]. Magnetic Resonance in Medicine, 2019, 82(3): 984–999. doi: 10.1002/mrm.27789

ZHU Xiaoxiang and BAMLER R. Superresolving SAR tomography for multidimensional imaging of urban areas: Compressive sensing-based tomoSAR inversion[J]. IEEE Signal Processing Magazine, 2014, 31(4): 51–58. doi: 10.1109/MSP.2014.2312098

王伟, 胡子英, 龚琳舒. MIMO雷达三维成像自适应Off-grid校正方法[J]. 电子与信息学报, 2019, 41(6): 1294–1301. doi: 10.11999/JEIT180145

WANG Wei, HU Ziying, and GONG Linshu. Adaptive off-grid calibration method for MIMO radar 3D imaging[J]. Journal of Electronics &Information Technology, 2019, 41(6): 1294–1301. doi: 10.11999/JEIT180145

金艳, 周磊, 姬红兵. 基于稀疏时频分布的跳频信号参数估计[J]. 电子与信息学报, 2018, 40(3): 663–669. doi: 10.11999/JEIT170525

JIN Yan, ZHOU Lei, and JI Hongbing. Parameter estimation of frequency-hopping signals based on sparse time-frequency distribution[J]. Journal of Electronics &Information Technology, 2018, 40(3): 663–669. doi: 10.11999/JEIT170525

FAZEL F, FAZEL M, and STOJANOVIC M. Random access compressed sensing for energy-efficient underwater sensor networks[J]. IEEE Journal on Selected Areas in Communications, 2011, 29(8): 1660–1670. doi: 10.1109/jsac.2011.110915

孙玉宝, 李欢, 吴敏, 等. 基于图稀疏正则化多测量向量模型的高光谱压缩感知重建[J]. 电子与信息学报, 2014, 36(12): 2942–2948. doi: 10.3724/SP.J.1146.2014.00566

SUN Yubao, LI Huan, WU Min, et al. Compressed sensing reconstruction of hyperspectral image using the graph sparsity regularized multiple measurement vector model[J]. Journal of Electronics &Information Technology, 2014, 36(12): 2942–2948. doi: 10.3724/SP.J.1146.2014.00566

LIAO Wenchao, HSIEH J, WANG Chengming, et al. Compressed sensing spectral domain optical coherence tomography with a hardware sparse-sampled camera[J]. Optics Letters, 2019, 44(12): 2955–2958. doi: 10.1364/OL.44.002955

DUARTE M F, DAVENPORT M A, TAKHAR D, et al. Single-pixel imaging via compressive sampling[J]. IEEE Signal Processing Magazine, 2008, 25(2): 83–91. doi: 10.1109/msp.2007.914730

余慧敏, 方广有. 压缩感知理论在探地雷达三维成像中的应用[J]. 电子与信息学报, 2010, 32(1): 12–16. doi: 10.3724/SP.J.1146.2009.00040

YU Huimin and FANG Guangyou. Research on compressive sensing based 3D imaging method applied to ground penetrating radar[J]. Journal of Electronics &Information Technology, 2010, 32(1): 12–16. doi: 10.3724/SP.J.1146.2009.00040

庄佳衍, 陈钱, 何伟基, 等. 基于压缩感知的动态散射成像[J]. 物理学报, 2016, 65(4): 040501. doi: 10.7498/aps.65.040501

ZHUANG Jiayan, CHEN Qian, HE Weiji, et al. Imaging through dynamic scattering media with compressed sensing[J]. Acta Physica Sinica, 2016, 65(4): 040501. doi: 10.7498/aps.65.040501

LI Bo, LIU Falin, ZHOU Chongbin, et al. Mixed sparse representation for approximated observation-based compressed sensing radar imaging[J]. Journal of Applied Remote Sensing, 2018, 12(3): 035015. doi: 10.1117/1.JRS.12.035015

LI Yunhui, WANG Xiaodong, WANG Zhi, et al. Modeling and image motion analysis of parallel complementary compressive sensing imaging system[J]. Optics Communications, 2018, 423: 100–110. doi: 10.1016/j.optcom.2018.04.018

施引文献

资源附件(0)

访问统计

图(9) / 表(1)

计量

文章访问数: 2595
HTML全文浏览量: 1023
PDF下载量: 63
被引次数: 0

1. 引言
2. 磁偶极子运动模型
3. 基于差量磁场的MCKF跟踪算法
3.1 差量磁场信号
3.2 基于MCKF的跟踪算法
4. 实验与算法验证
4.1 仿真实验
4.2 地磁背景跟踪实验
5. 结束语

姓名
邮箱
手机号码
标题
留言内容
验证码

留言板

基于压缩感知高反光成像技术研究

doi: 10.11999/JEIT190512

通讯作者:
赵首博　shoubozh@126.com

计量

Research on High Reflective Imaging Technology Based on Compressed Sensing

1. 引言

2. 磁偶极子运动模型

3. 基于差量磁场的MCKF跟踪算法

3.1 差量磁场信号

3.2 基于MCKF的跟踪算法

4. 实验与算法验证

4.1 仿真实验

4.2 地磁背景跟踪实验

5. 结束语

计量

目录

1. 引言

2. 磁偶极子运动模型

3. 基于差量磁场的MCKF跟踪算法

3.1 差量磁场信号

3.2 基于MCKF的跟踪算法

4. 实验与算法验证

4.1 仿真实验

4.2 地磁背景跟踪实验

5. 结束语

CS采样率	OMP算法		SAMP算法
CS采样率	MSE	PSNR	MSE	PSNR
0.30	0.0158	66.1481	0.0166	65.9403
0.35	0.0150	66.3739	0.0154	66.2512
0.40	0.0143	66.5625	0.0142	66.6061
0.45	0.0139	66.7149	0.0135	66.8149
0.50	0.0136	66.7928	0.0133	66.9078

留言板

基于压缩感知高反光成像技术研究

doi: 10.11999/JEIT190512

通讯作者: 赵首博 shoubozh@126.com

计量

出版历程

Research on High Reflective Imaging Technology Based on Compressed Sensing

1. 引言

2. 磁偶极子运动模型

3. 基于差量磁场的MCKF跟踪算法

3.1 差量磁场信号

3.2 基于MCKF的跟踪算法

4. 实验与算法验证

4.1 仿真实验

4.2 地磁背景跟踪实验

5. 结束语

计量

出版历程

目录

1. 引言

2. 磁偶极子运动模型

3. 基于差量磁场的MCKF跟踪算法

3.1 差量磁场信号

3.2 基于MCKF的跟踪算法

4. 实验与算法验证

4.1 仿真实验

4.2 地磁背景跟踪实验

5. 结束语

通讯作者:
赵首博　shoubozh@126.com