CN114758219B - 基于光谱数据和红外温度数据融合的尾迹识别方法 - Google Patents
基于光谱数据和红外温度数据融合的尾迹识别方法 Download PDFInfo
- Publication number
- CN114758219B CN114758219B CN202210658792.6A CN202210658792A CN114758219B CN 114758219 B CN114758219 B CN 114758219B CN 202210658792 A CN202210658792 A CN 202210658792A CN 114758219 B CN114758219 B CN 114758219B
- Authority
- CN
- China
- Prior art keywords
- trail
- data
- infrared
- spectral
- position information
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 230000003595 spectral effect Effects 0.000 title claims abstract description 85
- 238000000034 method Methods 0.000 title claims abstract description 38
- 230000004927 fusion Effects 0.000 title claims abstract description 19
- 238000003672 processing method Methods 0.000 claims abstract description 34
- 238000001514 detection method Methods 0.000 claims abstract description 21
- 230000000694 effects Effects 0.000 claims abstract description 10
- 238000004364 calculation method Methods 0.000 claims abstract description 9
- 239000011159 matrix material Substances 0.000 claims description 30
- 238000001228 spectrum Methods 0.000 claims description 24
- 239000013598 vector Substances 0.000 claims description 10
- 230000009467 reduction Effects 0.000 claims description 8
- 238000007781 pre-processing Methods 0.000 claims description 6
- 238000001914 filtration Methods 0.000 claims description 5
- 238000009499 grossing Methods 0.000 claims description 5
- 238000000513 principal component analysis Methods 0.000 claims description 5
- 238000012847 principal component analysis method Methods 0.000 claims description 4
- 230000011218 segmentation Effects 0.000 claims description 4
- 238000007619 statistical method Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 238000012163 sequencing technique Methods 0.000 claims description 2
- 239000000126 substance Substances 0.000 claims 1
- 230000009286 beneficial effect Effects 0.000 abstract 1
- 230000010354 integration Effects 0.000 abstract 1
- 239000013535 sea water Substances 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 4
- 230000002159 abnormal effect Effects 0.000 description 2
- 239000000498 cooling water Substances 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000003331 infrared imaging Methods 0.000 description 1
- 238000002329 infrared spectrum Methods 0.000 description 1
- 238000002620 method output Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000011144 upstream manufacturing Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
- G01J5/0022—Radiation pyrometry, e.g. infrared or optical thermometry for sensing the radiation of moving bodies
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/22—Matching criteria, e.g. proximity measures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/25—Fusion techniques
- G06F18/251—Fusion techniques of input or preprocessed data
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Evolutionary Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Artificial Intelligence (AREA)
- Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Radiation Pyrometers (AREA)
Abstract
本发明属于水下目标探测技术领域，公开了一种基于光谱数据和红外温度数据融合的尾迹识别方法，其包括如下步骤：基于已采集的尾迹红外温度数据样本，用红外热像仪获得的基于时间序列的红外图像提取尾迹目标温度数据，得到尾迹红外位置信息数据；基于已采集的尾迹光谱数据样本，用光谱仪采集的光谱数据进行处理，得到尾迹光谱位置信息数据；将尾迹红外位置信息数据与尾迹光谱位置信息数据进行余弦相似度计算；若计算得到的匹配度值大于或等于预设匹配度阈值，表明面向尾迹识别的光谱数据处理方法以及红外图像处理方法融合效果达到最优；否则不断地调节参数，直到二者融合效果达到最优，最后对目标尾迹识别。本发明利于提高水下目标探测的探测精度。
Description
技术领域
本发明属于水下目标探测技术领域，特别涉及一种基于光谱数据和红外温度数据融合的尾迹识别方法。
背景技术
水下目标探测根据其检测的信号源可以分为两类：直接探测和间接探测。声呐探测一直以来是最有效的直接探测手段，但随着水下目标降噪技术的不断发展，极大地增加了声纳探测的难度。非声探测手段应运而生，这种方法使用灵活，拥有与航空载机平台性能相关的动态探测能力，因而具有广阔的应用前景。非声的直接探测手段主要是通过检测水下目标本身的红外特征、光学特征、雷达特征和磁场特征，探测发现水下目标，利用红外成像技术进行尾流探测，已成为世界上重点发展的一项侦察技术。
水下目标高温冷却水排放入海使航行区域呈现“热尾流”特性，热尾流成为暴露水下目标位置、航向等重要信息的重要因素。目前红外探测系统的最低温度分辨率可以达到0.001℃，足以探测到海面异常的温度，获取或判断导航信息，识别水下目标的位置、速度和方向。
水下目标冷却水的排放使周围海水迅速升温，由于温度升高，密度相对减小，在扩散的过程中会逐渐向上浮升。若距离较近，会随时间浮升到自由水表面，由于海面恢复到原态需要较长时间，因此，水下目标运动过的海域由于热尾迹的存在会在海表面形成热尾迹。当热尾迹浮升不受任何抑制时，海表会形成大面积热尾迹，其特征原则上会被探测到。通过利用热尾迹目标和海洋背景的温差，可以实现水下航行器的红外探测。
水下目标运动过程中，海水冷却核装置排放出大量的温热尾迹，这些轨迹的特征与海水的性质相关。水下目标扰动温度密度分层，海水浮升形成冷尾迹温差信号，与温度密度均匀海水相比，海表温变区域显著增大、尾迹温差由6.13mK增大到 84 mK，可以通过海表上游冷尾迹特征判断是否存在水下目标。红外图像是通过“测量”物体向外辐射的热量而获得的，虽然可以探测到海洋背景与尾迹的温度差异，但是具有分辨率差、信噪比低、视觉效果模糊等缺点，因而，影响了水下目标尾迹的探测效率以及识别精度。
发明内容
本发明的目的在于提出一种基于光谱数据和红外温度数据融合的尾迹识别方法，在利用红外数据识别尾流的基础上同时综合了光谱数据，通过对光谱数据以及红外温度数据的处理与对比，利于提高水下目标的探测与识别准确度。
本发明为了实现上述目的，采用如下技术方案：
基于光谱数据和红外温度数据融合的尾迹识别方法，包括如下步骤：
步骤1. 基于已采集的尾迹红外温度数据样本，用红外热像仪获得的基于时间序列的原始红外图像提取尾迹目标温度数据，进一步得到尾迹红外位置信息数据；
步骤2. 基于已采集的尾迹光谱数据样本，用光谱仪采集的光谱数据进行处理，得到尾迹光谱位置信息数据；
步骤3. 将步骤1得到的尾迹红外位置信息数据与步骤2得到的尾迹光谱位置信息数据进行余弦相似度计算，进而得到它们之间的匹配度；
判断步骤3得到的匹配度值与预设匹配度阈值的大小关系，若该匹配度值小于预设匹配度阈值，则转到步骤4；若该匹配度值大于或等于预设匹配度阈值，则转到步骤5；
步骤4. 比较尾迹红外位置信息数据与尾迹光谱位置信息数据的数量；
经过比较，若尾迹红外位置信息数据的数量小于尾迹光谱位置信息数据的数量，则需要进一步改进步骤1中红外图像处理方法的帧间差分法的背景选取；
经过比较，若尾迹红外位置信息数据的数量大于或等于尾迹光谱位置信息数据的数量，则需要进一步改进步骤2中光谱数据处理方法的主成分分析n值选取；
参数调节完成后，并重复执行上述步骤1至步骤3；
步骤5. 判定为此目标区域为尾迹目标区域；此时，面向尾迹识别的光谱数据处理方法以及红外图像处理方法融合效果达到最优；
步骤6. 利用融合后的红外图像处理方法以及光谱数据处理方法，对目标尾迹进行识别。
本发明具有如下优点：
如上所述，本发明述及了一种基于光谱数据和红外温度数据融合的尾迹识别方法，该尾迹识别方法通过红外相机输出基于时间序列的红外图像温度数据对尾迹进行目标识别，通过红外光谱传感器获取长波红外波段的光谱数据，可以反演出海水的发射率，进行谱段选择和特征变换；本发明通过融合时间序列、空间维度的红外数据和光谱数据，能够使得海洋背景与尾迹之间的差别加大，进而提高水下目标尾迹的识别效率及精度。
附图说明
图1为本发明实施例中基于光谱和红外温度数据融合识别水下目标尾迹的方法的流程图。
图2为本发明实施例中利用红外热像仪采集的原始红外图像示意图。
图3为本发明实施例中红外背景图像示意图。
图4为本发明实施例中红外图像帧间差分后结果示意图。
图5为本发明实施例中红外图像阈值分割结果示意图。
图6为本发明实施例中利用光谱仪采集一个像元的光谱数据示意图。
图7为本发明实施例中利用光谱仪采集多个像元不同波段数据示意图。
具体实施方式
下面结合附图以及具体实施方式对本发明作进一步详细说明：
如图1所示，基于光谱数据和红外温度数据融合的尾迹识别方法，包括如下步骤：
步骤1. 基于已采集的尾迹红外温度数据样本，用红外热像仪获得的基于时间序列的原始红外图像提取尾迹目标温度数据，进一步得到尾迹红外位置信息数据。
该步骤1具体为：
步骤1.1. 用红外热像仪采集基于时间序列的原始红外图像温度数据Fi(x,y)；其中，i表示第i帧图像，第i帧图像的大小为M×N，(x,y)表示尾迹红外位置信息。
步骤1.2. 用小波降噪方法对步骤1.1得到的基于时间序列的原始红外图像温度数据Fi(x,y)进行预处理，得到处理后的红外图像温度数据fi(x,y)。
步骤1.3. 将步骤1.2中经过降噪处理后的红外图像温度数据fi(x,y)，进行背景选取，得到背景数据fk(x,y) ，即第k帧红外图像温度数据。
将步骤1.2中经过降噪处理后的红外图像温度数据fi(x,y)与步骤1.3中的背景数据fk(x,y)相减，突出尾迹目标信息，得到红外图像D(x,y)。
步骤1.3中，将降噪处理后的红外图像温度数据fi(x，y)进行帧间差分法，得到红外图像D(x,y)，帧间差分法包括以下步骤：
选取基于时间序列的第k帧红外图像作为背景数据fk(x,y)；选取基于时间序列的第i帧红外图像作为尾迹检测的红外图像温度数据fi(x,y)。
帧间差分法的公式如下：
当| fk(x,y)- fi(x,y)|＞T时，D(x,y)=1，否则，D(x,y)=0，阈值T由阈值分割法确定。
步骤1.4. 根据红外图像D(x,y)识别尾迹，提取出包含尾迹信息的二值图像d(x,y)即矩阵D，得到尾迹红外位置信息数据{(x1,y1) ,(x2 ,y2) ,…,(xd ,yd)}。
其中，(x1,y1)、(x2 ,y2) ,…, (xd ,yd)分别表示第1,2,…,d个尾迹红外位置信息。
步骤2. 基于已采集的尾迹光谱数据样本，用光谱仪采集的光谱数据进行处理，得到尾迹光谱位置信息数据。该步骤2具体为：
步骤2.1. 用光谱仪采集的光谱数据X是一个 P 维随机变量。
X=( x1,x2,…, xN) =( X1, X2,…, XP)T；
其中， N 是光谱数据的像元数目, P是波段数目，xi表示光谱数据的第i个像元，i=1, 2,…, N, Xk表示图像的第 k 个波段，k=1, 2,…, P。
步骤2.2. 将光谱数据X进行预处理，使用Savitzky-Golay卷积平滑滤波, 提高光
谱的平滑性，并降低噪声干扰，公式如下：
。
其中，
为Savitzky-Golay卷积平滑滤波所得结果；
表示k点处经平
滑的光谱值，w表示滤波窗口宽度，xk+i表示第k+i段波段光谱数据。
步骤2.3. 通过主成分分析法分析光谱数据特征贡献率，通过对贡献率的排序选取前n个主成分作为参考波段；主成分分析法包括以下步骤：
步骤2.3.1. 计算数据协方差矩阵，并获得协方差矩阵∑x的特征值。其中，λ1,λ2,…,λP为协方差矩阵∑x的特征值，且λ1≥λ2≥…≥λP≥0。
定义m表示X 的均值向量，m=∑N i=1Xi/N。
其中，Xi表示第i个像元的光谱数据，∑x表示协方差矩阵；∑x =E{(X-m) (X-m)T}。
步骤2.3.2. 将对应的特征值按照大小排序，最大的特征值是第一主成分，第二大的特征值是第二主成分，Vn为选取前n个主成分作为参考波段。
Y=( Y1, Y2, …, YP)T=AX；A=( t1, t2, …, tP)；Vn=∑n i=1λi/∑P i=1λi。
其中，A表示主成分变换矩阵，Y1, Y2, …, YP 表示X的第1主成分，第2主成分，…,第p主成分，t1, t2, …, tP为特征值λ1,λ2, …,λP对应的单位特征向量。
步骤2.4. 设δi1是∑x对角线上的值, 根据ρ(Yk,Xi)计算主成分Yk与原来第i个波段 Xi的相关系数，Vi是进行波段选择的主要依据，选择最大Vi值对应的波段作为最优波段。
ρ(Yk,Xi)= 
tik/ 
k,i=1,2,…,P；Vi=∑m k=1ρ2(Yk,Xi)= ∑m k=1λkt2 ik/δi1。
其中，λk表示∑x的特征值，tik表示单位特征向量ti第k个值。
步骤2.5. 基于步骤2.4得到的最优波段的光谱数据构建光谱曲线。
步骤2.6. 在步骤2.5构建的光谱曲线中基于统计的方法识别与全局背景光谱存在较大差异的像元，提取尾迹光谱位置信息数据{(x1’,y1’) ,(x2’ ,y2’),…,(xd’,yd’)}；其中，(x1’,y1’)、(x2’ ,y2’) ,…, (xd’ ,yd’)分别表示第1,2,…,d个尾迹光谱位置信息。
步骤3. 将步骤1得到的尾迹红外位置信息数据与步骤2得到的尾迹光谱位置信息数据进行余弦相似度计算，进而得到它们之间的匹配度。
判断经过计算得到的匹配度值与预设匹配度阈值的大小关系，若匹配度值小于预设匹配度阈值，则转到步骤4；若匹配度值大于或等于预设匹配度阈值，则转到步骤5。
本实施例中预设匹配度阈值例如设定为0.9，当然，该值只是示例性的，例如还可以根据需要，将预设匹配度阈值设定为0.85、0.95等值。
该步骤3具体为：
步骤3.1. 将提取的尾迹光谱位置信息数据，构成一个大小为M×N的矩阵，其中，尾迹位置信息置1，其他位置置0，构成一个光谱位置矩阵G。
步骤3.2. 将光谱位置矩阵G与红外位置矩阵D进行余弦相似度计算。
将光谱位置矩阵G表示为向量G=(g1,g2，…，gd)，将红外位置矩阵D表示为向量D=(d1,d2，…，dd)，则光谱位置矩阵G和红外位置矩阵D的余弦相似度为：
其中，g1=(x1’,y1’)，g2=(x2’,y2’) ，…，gd= (xd’,yd’)。
d1=(x1 ,y1),d1=(x2 ,y2)，…，dd=(xd ,yd)。
如果两个矩阵的余弦相似度，即匹配度值大于或等于预设匹配度阈值，则判定尾迹识别成功，此时面向尾迹识别的光谱数据处理方法以及红外图像处理方法融合效果达到最优。
步骤4. 比较尾迹红外位置信息数据与尾迹光谱位置信息数据的数量。
经过比较，若尾迹红外位置信息数据的数量小于尾迹光谱位置信息数据的数量，则需要进一步改进步骤1中红外图像处理方法的帧间差分法的背景选取。
改进帧间差分法第k帧背景的选取，每次参数调整在选取背景帧数k值基础上减50。
经过比较，若尾迹红外位置信息数据的数量大于或等于尾迹光谱位置信息数据的数量，则需要进一步改进步骤2中光谱数据处理方法的主成分分析n值选取。
每次参数调整在n值基础上加1，重新进行计算参考波段选取以及最优波段选取。
参数调节完成后，重复执行上述步骤1至步骤3；通过不断地调节参数，使得面向尾迹识别的光谱数据处理方法以及红外图像处理方法融合效果达到最优。
步骤5. 判定为此目标区域为尾迹目标区域；此时，面向尾迹识别的光谱数据处理方法以及红外图像处理方法融合效果达到最优，进而精确、实时识别尾迹。
步骤6. 利用融合后的红外图像处理方法以及光谱数据处理方法，对目标尾迹进行识别。
该步骤6具体为：利用参数调整好的红外图像处理方法和光谱数据处理方法，分别获取尾迹红外位置信息数据G’和尾迹光谱位置信息数据D’。
将数据G’和D’合并，组成一张尾迹位置信息二值化图片，得到尾迹识别结果。
下面根据本发明方法工作流程，对尾迹进行目标识别，如图2至图7所示。
用红外热像仪进行数据收集，图2为依据红外热像仪采集温度数据显示的原始红外图像，将原始红外图像进行预处理，采用小波降噪后得到的图像，用帧间差分法对红外原始温度数据进行处理，减去图3背景图，得到图4的结果。由图4中的结果看出，尾迹目标得以明显凸显。将帧间差分法处理后的图像采用OTSU算法进行阈值分割，如图5所示，可以凸显尾迹目标，根据阈值分割后的二值化图像可进一步得到尾迹红外位置信息数据。
在红外热像仪采集尾迹红外温度数据的同时，用光谱仪进行光谱数据采集，得到三维的光谱数据，将光谱数据采用Savitzky-Golay卷积平滑滤波进行预处理，去除噪声等影响因素。采用主成分分析选取n个参考波段，如图6所示，将所有波段与参考波段对比，得出包含大部分参考波段信息的最优波段，如图7所示，选取了最优波段并构建光谱曲线。用统计学方法对选取的谱线进行处理，分析出相对于海面背景异常的区域，提取该目标像元的光谱特征。
将红外位置信息数据与光谱位置信息数据进行余弦相似度计算，若余弦相似度大于或等于0.9，将红外位置信息矩阵和光谱位置信息矩阵合并，此时识别的目标区域即为尾迹区域，此时，面向尾迹识别的光谱数据处理方法以及红外图像处理方法融合效果达到最优。若余弦相似度小于0.9，则通过不断调整参数，使得上述两种方法融合效果达到最优。
当然，以上说明仅仅为本发明的较佳实施例，本发明并不限于列举上述实施例，应当说明的是，任何熟悉本领域的技术人员在本说明书的教导下，所做出的所有等同替代、明显变形形式，均落在本说明书的实质范围之内，理应受到本发明的保护。
Claims (7)
1.基于光谱数据和红外温度数据融合的尾迹识别方法，其特征在于，
包括如下步骤：
步骤1. 基于已采集的尾迹红外温度数据样本，用红外热像仪获得的基于时间序列的原始红外图像提取尾迹目标温度数据，得到尾迹红外位置信息数据；
步骤2. 基于已采集的尾迹光谱数据样本，用光谱仪采集的光谱数据进行处理，得到尾迹光谱位置信息数据；
步骤3. 将步骤1得到的尾迹红外位置信息数据与步骤2得到的尾迹光谱位置信息数据进行余弦相似度计算，得到它们之间的匹配度；
判断所述匹配度值与预设匹配度阈值的大小关系，若该匹配度值小于预设匹配度阈值，则转到步骤4；若该匹配度值大于或等于预设匹配度阈值，则转到步骤5；
步骤4. 比较尾迹红外位置信息数据与尾迹光谱位置信息数据的数量；
经过比较，若尾迹红外位置信息数据的数量小于尾迹光谱位置信息数据的数量，则进一步改进步骤1中红外图像处理方法的帧间差分法的背景选取；
经过比较，若尾迹红外位置信息数据的数量大于或等于尾迹光谱位置信息数据的数量，则进一步改进步骤2中光谱数据处理方法的主成分分析n值选取；
参数调节完成后，并重复执行上述步骤1至步骤3；
步骤5. 判定为此目标区域为尾迹目标区域；此时，面向尾迹识别的光谱数据处理方法以及红外图像处理方法融合效果达到最优；
步骤6. 利用融合后的红外图像处理方法以及光谱数据处理方法，对目标尾迹进行识别。
2.根据权利要求1所述的尾迹识别方法，其特征在于，
所述步骤1具体为：
步骤1.1. 用红外热像仪采集基于时间序列的原始红外图像温度数据Fi(x,y)；其中，i表示第i帧图像，第i帧图像的大小为M×N，(x,y)表示尾迹红外位置信息；
步骤1.2. 用小波降噪方法对步骤1.1得到的基于时间序列的原始红外图像温度数据Fi(x,y)进行预处理，得到处理后的红外图像温度数据fi(x,y)；
步骤1.3. 将步骤1.2中经过降噪处理后的红外图像温度数据fi(x,y)，进行背景选取，得到背景数据fk(x,y) ，即第k帧红外图像温度数据；
将步骤1.2中经过降噪处理后的红外图像温度数据fi(x,y)与步骤1.3中的背景数据fk(x,y)相减，突出尾迹目标信息，得到红外图像D(x,y)；
步骤1.4. 根据红外图像D(x,y)识别尾迹，提取出包含尾迹信息的二值图像d(x,y)即矩阵D，得到尾迹红外位置信息数据{(x1,y1) ,(x2 ,y2) ,...,(xd ,yd)}；
其中，(x1,y1)、(x2 ,y2) ,..., (xd ,yd)分别表示第1,2,...,d个尾迹红外位置信息。
3.根据权利要求2所述的尾迹识别方法，其特征在于，
所述步骤1.3中，将降噪处理后的红外图像温度数据fi(x，y)进行帧间差分法，得到红外图像D(x,y)，帧间差分法包括以下步骤：
选取基于时间序列的第k帧红外图像作为背景数据fk(x,y)；选取基于时间序列的第i帧红外图像作为尾迹检测的红外图像温度数据fi(x,y)；
帧间差分法的公式如下：
当| fk(x,y)- fi(x,y)|＞T时，D(x,y)=1，否则，D(x,y)=0，阈值T由阈值分割法确定。
4.根据权利要求3所述的尾迹识别方法，其特征在于，
所述步骤2具体为：
步骤2.1. 用光谱仪采集的光谱数据X是一个 P 维随机变量；
X=( x1,x2, …, xN) =( X1, X2, …, XP)T；
其中， N 是光谱数据的像元数目, P是波段数目，xi表示光谱数据的第i个像元，i=1,2, …, N, Xk表示图像的第 k 个波段，k=1, 2, …, P；
步骤2.2. 将光谱数据X进行预处理，使用Savitzky-Golay卷积平滑滤波,公式如下：
；
其中，
为Savitzky-Golay卷积平滑滤波所得结果；
表示k点处经平滑的光谱值，w表示滤波窗口宽度，xk+i表示第k+i段波段光谱数据；
步骤2.3. 通过主成分分析法分析光谱数据特征贡献率，通过对贡献率的排序选取前n个主成分作为参考波段；主成分分析法包括以下步骤：
步骤2.3.1. 计算数据协方差矩阵，并获得协方差矩阵∑x的特征值；其中，λ1,λ2, …,λP为协方差矩阵∑x的特征值，且λ1≥λ2≥…≥λP≥0；
定义m表示X 的均值向量，m=∑N i=1Xi/N；
其中，Xi表示第i个像元的光谱数据，∑x表示协方差矩阵；∑x =E{(X-m) (X-m)T}；
步骤2.3.2. 将对应的特征值按照大小排序，最大的特征值是第一主成分，第二大的特征值是第二主成分，Vn为选取前n个主成分作为参考波段；
Y=( Y1, Y2, …, YP)T=AX；A=( t1, t2, …, tP)；Vn=∑n i=1λi/∑P i=1λi；
其中，A表示主成分变换矩阵，Y1, Y2, …, YP 表示X的第1主成分，第2主成分，…,第p主成分，t1, t2, …, tP为特征值λ1,λ2, …,λP对应的单位特征向量；
步骤2.4. 设δi1是∑x对角线上的值, 根据ρ(Yk,Xi)计算主成分Yk与原来第i个波段 Xi的相关系数，Vi是进行波段选择的主要依据，选择最大Vi值对应的波段作为最优波段；
ρ(Yk,Xi)= 
tik/ 
k,i=1,2,…,P；Vi=∑m k=1ρ2(Yk,Xi)= ∑m k=1λkt2 ik/δi1；
其中，λk表示∑x的特征值，tik表示单位特征向量ti第k个值；
步骤2.5. 基于步骤2.4得到的最优波段的光谱数据构建光谱曲线；
步骤2.6. 在构建的光谱曲线中基于统计的方法识别与全局背景光谱存在较大差异的像元，提取尾迹光谱位置信息数据{(x1’,y1’) ,(x2’,y2’),…,(xd’,yd’)}； 其中，(x1’,y1’)、(x2’,y2’) ,..., (xd’ ,yd’)分别表示第1,2,...,d个尾迹光谱位置信息。
5.根据权利要求4所述的尾迹识别方法，其特征在于，
所述步骤3具体为：
步骤3.1. 将提取的尾迹光谱位置信息数据，构成一个大小为M×N的矩阵，其中，尾迹位置信息置1，其他位置置0，构成一个光谱位置矩阵G；
步骤3.2. 将光谱位置矩阵G与红外位置矩阵D进行余弦相似度计算；
将光谱位置矩阵G表示为向量G=(g1,g2，…，gd)，将红外位置矩阵D表示为向量D=(d1,d2，…，dd)，则光谱位置矩阵G和红外位置矩阵D的余弦相似度为：
其中，g1=(x1’,y1’)，g1=(x2’,y2’) ，…，gd= (xd’,yd’)；
D1=(x1 ,y1),d2=(x2 ,y2)，…，dd= (xd ,yd)。
6.根据权利要求5所述的尾迹识别方法，其特征在于，
所述步骤4具体为：
若尾迹红外位置信息数据的数量小于尾迹光谱位置信息数据的数量，则红外数据改进帧间差分法第k帧背景的选取，每次参数调整在选取背景帧数k值基础上减50；
否则，则光谱数据改进主成分分析中n值参数选取，每次参数调整在n值的基础上加1，重新进行计算参考波段选取以及最优波段选取；
通过不断地调节光谱数据处理方法和红外图像处理方法的参数，使得面向尾迹识别的光谱数据处理方法以及面向尾迹识别的红外图像处理方法融合效果达到最优。
7.根据权利要求6所述的尾迹识别方法，其特征在于，
所述步骤6具体为：利用参数调整好的红外图像处理方法和光谱数据处理方法，分别获取尾迹红外位置信息数据G’和尾迹光谱位置信息数据D’；
将数据G’和D’合并，组成一张尾迹位置信息二值化图片，得到尾迹识别结果。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210658792.6A CN114758219B (zh) | 2022-06-13 | 2022-06-13 | 基于光谱数据和红外温度数据融合的尾迹识别方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210658792.6A CN114758219B (zh) | 2022-06-13 | 2022-06-13 | 基于光谱数据和红外温度数据融合的尾迹识别方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114758219A CN114758219A (zh) | 2022-07-15 |
| CN114758219B true CN114758219B (zh) | 2022-09-02 |
Family
ID=82336463
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210658792.6A Active CN114758219B (zh) | 2022-06-13 | 2022-06-13 | 基于光谱数据和红外温度数据融合的尾迹识别方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114758219B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115327657B (zh) * | 2022-10-13 | 2022-12-20 | 维飒科技(西安)有限公司 | 探测目标定位方法及装置 |
| CN116258932B (zh) * | 2023-03-09 | 2024-06-14 | 中国人民解放军海军潜艇学院 | 一种用于水下运动目标尾迹的红外融合探测方法及系统 |
| CN116735527B (zh) * | 2023-06-09 | 2024-01-05 | 湖北经济学院 | 一种近红外光谱优化方法、装置、系统以及存储介质 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE10359868B3 (de) * | 2003-12-18 | 2005-06-16 | Deutsches Zentrum für Luft- und Raumfahrt e.V. | Verfahren zur Erkennung von sich hinter einem Flugzeug bildendem/n Kondensstreifen sowie Einrichtung zur Durchführung des Verfahrens |
| CN113673385B (zh) * | 2021-08-06 | 2024-06-21 | 南京理工大学 | 基于红外图像的海面舰船检测方法 |
| CN114169264A (zh) * | 2021-11-30 | 2022-03-11 | 哈尔滨工程大学 | 一种舰船尾迹红外仿真方法 |
| CN113920283B (zh) * | 2021-12-13 | 2022-03-08 | 中国海洋大学 | 基于聚类分析和特征过滤的红外图像尾迹检测及提取方法 |
-
2022
- 2022-06-13 CN CN202210658792.6A patent/CN114758219B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN114758219A (zh) | 2022-07-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114758219B (zh) | 基于光谱数据和红外温度数据融合的尾迹识别方法 | |
| CN107301661B (zh) | 基于边缘点特征的高分辨率遥感图像配准方法 | |
| Nilsson et al. | Localization of corresponding points in fingerprints by complex filtering | |
| CN111667470B (zh) | 一种基于数字图像的工业管道探伤内壁检测方法 | |
| CN106815819B (zh) | 多策略粮虫视觉检测方法 | |
| CN107220628B (zh) | 红外干扰源检测的方法 | |
| CN111208479B (zh) | 一种降低深度网络检测中虚警概率的方法 | |
| CN111353496B (zh) | 一种红外弱小目标实时检测方法 | |
| CN116152115B (zh) | 基于计算机视觉的垃圾图像去噪处理方法 | |
| CN109063669B (zh) | 一种基于图像识别的桥区船舶航行态势分析方法及装置 | |
| CN116665095B (zh) | 一种运动舰船检测方法、系统、存储介质和电子设备 | |
| CN112308873A (zh) | 多尺度Gabor小波PCA融合图像边缘检测方法 | |
| He et al. | Multiscale local gray dynamic range method for infrared small-target detection | |
| Ma et al. | A method for infrared sea-sky condition judgment and search system: Robust target detection via PLS and CEDoG | |
| CN115861226A (zh) | 基于特征值梯度变化的深度神经网络智能识别表面缺陷方法 | |
| CN115359258A (zh) | 一种成分不确定度测量的弱小目标检测方法及系统 | |
| CN113222898B (zh) | 基于多元统计量和深度学习的双航过sar图像痕迹检测方法 | |
| Wu et al. | Research on crack detection algorithm of asphalt pavement | |
| CN116934808A (zh) | 一种基于水面漂浮物目标跟踪的河流表面流速测量方法 | |
| CN116243289A (zh) | 一种基于成像声呐的无人艇水下目标智能识别方法 | |
| CN115797366A (zh) | 一种基于霍夫变换的移动目标提取图像处理方法和系统 | |
| CN115345845A (zh) | 基于方向梯度直方图的特征融合烟幕干扰效能评估及处理方法及电子设备 | |
| CN112686222B (zh) | 星载可见光探测器对船舶目标检测的方法和系统 | |
| CN114429593A (zh) | 基于快速导向滤波的红外小目标检测方法及其应用 | |
| CN111401195B (zh) | 一种基于多波段红外图像的海面目标检测方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |