CN108052957A - 一种航天器目标快速识别方法 - Google Patents
一种航天器目标快速识别方法 Download PDFInfo
- Publication number
- CN108052957A CN108052957A CN201711083754.8A CN201711083754A CN108052957A CN 108052957 A CN108052957 A CN 108052957A CN 201711083754 A CN201711083754 A CN 201711083754A CN 108052957 A CN108052957 A CN 108052957A
- Authority
- CN
- China
- Prior art keywords
- mrow
- spacecraft
- feature
- image
- point
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 52
- 238000012545 processing Methods 0.000 claims abstract description 27
- 230000000007 visual effect Effects 0.000 claims abstract description 26
- 238000000605 extraction Methods 0.000 claims abstract description 22
- 238000004364 calculation method Methods 0.000 claims abstract description 20
- 238000011084 recovery Methods 0.000 claims abstract description 16
- 238000004458 analytical method Methods 0.000 claims abstract description 9
- 239000011159 matrix material Substances 0.000 claims description 68
- 230000006870 function Effects 0.000 claims description 31
- 230000009466 transformation Effects 0.000 claims description 31
- 238000013139 quantization Methods 0.000 claims description 12
- 238000005316 response function Methods 0.000 claims description 11
- 230000001629 suppression Effects 0.000 claims description 9
- 238000003384 imaging method Methods 0.000 claims description 8
- 239000013598 vector Substances 0.000 claims description 8
- 238000003064 k means clustering Methods 0.000 claims description 7
- 238000013507 mapping Methods 0.000 claims description 7
- 230000011218 segmentation Effects 0.000 claims description 6
- 230000003287 optical effect Effects 0.000 claims description 5
- 238000002187 spin decoupling employing ultra-broadband-inversion sequences generated via simulated annealing Methods 0.000 claims description 5
- 238000007781 pre-processing Methods 0.000 claims description 3
- 238000000354 decomposition reaction Methods 0.000 claims description 2
- 239000000203 mixture Substances 0.000 claims 1
- 239000002994 raw material Substances 0.000 claims 1
- 108020001568 subdomains Proteins 0.000 claims 1
- 230000008569 process Effects 0.000 abstract description 11
- 230000004438 eyesight Effects 0.000 abstract description 7
- 230000002708 enhancing effect Effects 0.000 abstract description 3
- 230000004313 glare Effects 0.000 abstract description 3
- 230000009286 beneficial effect Effects 0.000 description 14
- 230000008901 benefit Effects 0.000 description 7
- 238000005457 optimization Methods 0.000 description 7
- 238000005286 illumination Methods 0.000 description 6
- 238000005070 sampling Methods 0.000 description 5
- 230000004075 alteration Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000009792 diffusion process Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 230000001965 increasing effect Effects 0.000 description 4
- 230000005764 inhibitory process Effects 0.000 description 4
- 230000015556 catabolic process Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000006731 degradation reaction Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000010606 normalization Methods 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000003702 image correction Methods 0.000 description 2
- 238000012634 optical imaging Methods 0.000 description 2
- 238000013519 translation Methods 0.000 description 2
- PCTMTFRHKVHKIS-BMFZQQSSSA-N (1s,3r,4e,6e,8e,10e,12e,14e,16e,18s,19r,20r,21s,25r,27r,30r,31r,33s,35r,37s,38r)-3-[(2r,3s,4s,5s,6r)-4-amino-3,5-dihydroxy-6-methyloxan-2-yl]oxy-19,25,27,30,31,33,35,37-octahydroxy-18,20,21-trimethyl-23-oxo-22,39-dioxabicyclo[33.3.1]nonatriaconta-4,6,8,10 Chemical compound C1C=C2C[C@@H](OS(O)(=O)=O)CC[C@]2(C)[C@@H]2[C@@H]1[C@@H]1CC[C@H]([C@H](C)CCCC(C)C)[C@@]1(C)CC2.O[C@H]1[C@@H](N)[C@H](O)[C@@H](C)O[C@H]1O[C@H]1/C=C/C=C/C=C/C=C/C=C/C=C/C=C/[C@H](C)[C@@H](O)[C@@H](C)[C@H](C)OC(=O)C[C@H](O)C[C@H](O)CC[C@@H](O)[C@H](O)C[C@H](O)C[C@](O)(C[C@H](O)[C@H]2C(O)=O)O[C@H]2C1 PCTMTFRHKVHKIS-BMFZQQSSSA-N 0.000 description 1
- 102000002274 Matrix Metalloproteinases Human genes 0.000 description 1
- 108010000684 Matrix Metalloproteinases Proteins 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006386 neutralization reaction Methods 0.000 description 1
- 101150050759 outI gene Proteins 0.000 description 1
- 230000003071 parasitic effect Effects 0.000 description 1
- 210000001747 pupil Anatomy 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
- G06F18/232—Non-hierarchical techniques
- G06F18/2321—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
- G06F18/23213—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions with fixed number of clusters, e.g. K-means clustering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/20—Image preprocessing
- G06V10/30—Noise filtering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/40—Extraction of image or video features
- G06V10/42—Global feature extraction by analysis of the whole pattern, e.g. using frequency domain transformations or autocorrelation
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Multimedia (AREA)
- Probability & Statistics with Applications (AREA)
- Image Analysis (AREA)
Abstract
本发明涉及一种航天器目标快速识别方法及装置,该方法包括:对获取的航天器观测图像进行杂光抑制增强处理,得到所述航天器观测图像的恢复图像;对所述恢复图像进行特征点及特征信息的提取;对所述特征信息进行优化处理,并通过语义分析和判别进行航天器目标的识别。通过杂光抑制及增强处理,能提高图像信噪比,增强图像质量。对观测图像的特征点以及不变矩特征信息和局部灰度特征信息的提取,能更有利于鉴别航天器型号。信息优化过程减轻计算量,有助于提高处理速度。根据图像中局部显著特征的特性,不同型号航天器的视觉特征和其语义概率值,然后依据后验概率值和经验阈值来判断语义信息,能够最大程度的保证视觉类似的航天器目标的识别。
Description
技术领域
本发明涉及图像识别技术领域,具体涉及一种航天器目标快速识别方法。
背景技术
随着天基空间观测探测技术的快速发展,对空间场景中观测探测技术要求也越来越高。现有识别方法依靠观测的椭圆度、直线度、周长比等特征向量作为后续识别的基础和前提输入,采用聚类或者二分类的方法对航天器进行识别。忽略了航天器成像环境杂光的干扰,杂光形成原理如图1所示,同时由于航天器自身的姿态变化使得特征信息获取有很大的不确定性和不完整性,导致最后的识别结果精度下降。
发明内容
本发明针对现有技术中存在的技术问题,提供航天器目标快速识别方法及装置,首先对航天器图像进行杂光纠正,增强图像质量,进而提取更为精细的判别性好的局部特征进行快速语义识别,来实现对卫星目标进行稳定、可靠的语义识别。
本发明采用的技术方案如下:
本发明一方面提供一种航天器目标快速识别方法,包括,
步骤1,对获取的航天器观测图像进行杂光抑制增强处理,得到所述航天器观测图像的恢复图像;
步骤2,对所述航天器观测图像的恢复图像进行特征点及特征信息的提取;所述特征信息包括所述特征点的不变矩特征信息和局部灰度特征信息;
步骤3,对所述特征信息进行优化处理,并通过语义分析和判别进行航天器目标的识别。
本发明的有益效果是:航天器在空中被动成像,由于光照原因会形成杂光,使得局部特征信息受到干扰,导致后续特征提取不准确,影响识别效果。本方法通过杂光抑制及增强处理,能提高图像信噪比,增强图像质量,利于后续目标快速识别;对观测图像的特征点以及包含特征点的不变矩特征和局部灰度特征的特征信息的提取,能更有利于鉴别航天器型号;信息优化过程减轻计算量,有助于提高处理速度。根据图像中局部显著特征的特性,不同型号航天器的视觉特征和其语义概率值,然后依据后验概率值和经验阈值来判断语义信息,能够最大程度的保证视觉类似的航天器目标的识别。
进一步,所述步骤1包括:
步骤101,根据下式对获取的所述航天器观测图像进行图像恢复,
式中y为航天器观测图像,为恢复图像,cos2γ是阴影系数,γ表示光轴和成像系统的夹角,SSTR是纯杂散光部分的点扩散函数矩阵,β为其权重系数;
步骤102,对纯杂散光部分的点扩散函数SSTR进行量化,并采用Frobenius误差度量量化形变引起的误差,然后根据下式得到矩阵
式中,E为单位正交阵,Λy为对航天器观测图像进行特征分解得到的对角阵;
步骤103,采用K-L变换和小波变换对所述矩阵的行和列进行去相关运算;
步骤104,对去相关后的矩阵进行量化,并利用稀疏矩阵变换对量化后的矩阵的行进行去相关,通过下式得到经过杂光抑制增强处理的所述航天器观测图像的恢复图像
式中,W为二维正交小波变换矩阵,为量化后的空变点扩散函数,Tk为根据吉文斯旋转定理得到的稀疏矩阵。
采用上述进一步方案的有益效果是根据航天器目标图像成像时杂散光不同的成因,建立了归一化的空变点扩散函数模型,并且对模型增加了光照的阴影系数,使模型更加符合光学成像系统的实际情况;Van Cittert迭代方法是线性反卷积过程,最终收敛为一个简单的逆滤波,适用于杂散光降质图像校正时的图像恢复问题,在初始迭代情况下原始清晰图像受到移变杂散光点扩散函数的影响最小;稀疏矩阵是指大多数元素为零元素,只有极少非零元素的矩阵。和传统的变换方式不同,稀疏矩阵变换(Sparse Matrix Transform,SMT)只用其有限数目点坐标的旋转,即吉文斯旋转的乘积就可以实现对既定目标矩阵的精确估计,且其运算和存储只需要针对非零元即可,在实现数据有效压缩的同时,又可对原始变换进行精确估计,作为估计正交变换的有效工具,吉文斯旋转可以有效降低反矩阵求解问题的计算复杂度。通过杂光抑制,增加航天器目标图像的可见度,同时纠正图像中的色差问题,提高图像质量。
进一步,所述步骤2,包括:
步骤201,对所述航天器观测图像的恢复图像进行候选特征点提取;
步骤202,利用SUSAN算子,对所述候选特征点进行区域分割,并对分割获得的候选特征点区域进行分析,去掉伪特征点,获得所述航天器观测图像的特征点;
步骤203,将所述特征点周围预设区域作为感兴趣区域,采用Sobel算子提取感兴趣区域的灰度特征信息,即为所述特征点的灰度特征信息;
步骤204,以所述特征点为中心,求取所述恢复图像的归一化中心矩,并依据二阶和三阶的归一化中心矩,构造多个不变矩,所述多个不变矩构成一组特征量,所述特征量即为所述特征点的不变矩特征信息。
采用上述进一步方案的有益效果是通过特征点的二次提取,去除伪特征点信息,提高特征点特征信息提取准确度。不变矩可以有效描述空间目标的统计特性,且对于图像的平移、缩放、旋转都不敏感,因此通过不变矩特征信息提取,增强对航天器观测图像处理时的几何不变性。
进一步,所述步骤201包括,
采用预设大小n*n的窗口(优选为5*5大小)对每个像元逐一探测,以所述像元为中心,分别计算所述像元四个方向(0°、45°、90°、135°)和相邻像素的灰度差的平方和,
其中k=[n/2]。
取四个所述平方和V1,V2,V3,V4中的最小值作为该像元的兴趣值,同时计算四个所述平方和的中值:将所述兴趣值大于中值的像元作为侯选特征点。
进一步,所述步骤202包括,
以候选特征点所在的圆形区域作为SUSAN算子的模板,选取所述模板中构成候选特征点的最大的子域作为所述候选特征点的第一区域,而将其他子域视为单一的背景区域;将所述模板分为多个同心的圆,第一区域的特征点位置为圆中心;定义Si是第i个同心圆上第一区域的对应圆弧,对于每一个圆,则有:
用计算信息量的方法获得候选特征点响应函数,将其定义为模板内像元n与平均信息量差的平方和,以模板内像元幅度值的累加值除以模板面积作为该模板的平均信息量gr(n),其计算方法如下:
其中,fk,l(n)表示内坐标(k,l)点的幅值,r表示圆形模板的半径,Ω表示模板的覆盖范围,特征点响应函数为:
选取函数值最大的点作为特征点,其余的候选特征点作为伪特征点被除去。
采用上述进一步方案的有益效果是去除候选特征点中的伪特征点,降低特征提取时的计算复杂度。
进一步,所述步骤203包括,
将特征点周围预设区域作为目标图像的感兴趣区域,采用Sobel算子计算感兴趣区域内每个像素点水平和垂直两个方向上的灰度梯度分量,获取感兴趣区域内每个像素点灰度梯度向量的幅值和幅角,即所述特征点的灰度特征信息。
采用上述进一步方案的有益效果是不变矩可以有效描述空间目标的统计特性,且对于图像的平移、缩放、旋转都不敏感,因此通过不变矩特征信息提取,增强对航天器观测图像处理时的几何不变性。
进一步,所述步骤203,还包括:采用高斯滤波器对以特征点为中心的感兴趣区域中的每个方向属性平面进行亚采样。
采用上述进一步方案的有益效果是增强局部灰度特征信息的鲁棒性。
进一步,所述步骤3,包括:
步骤310,对所述特征信息进行优化,将所有高维特征数据映射到低维空间,获得优化后的特征信息;
步骤320,对优化后的特征信息进行K-means聚类得到k个聚类中心,每个聚类中心对应一个视觉单词,获取视觉单词w、所述航天器观测图像的恢复图像及其潜在语义主题z所对应的联合概率分布,进而根据所述联合概率分布计算概率值
步骤330,采用EM学习方法获取参数阈值,根据参数阈值来判别所述航天器观测图像的恢复图像所属的语义信息;根据不同类型航天器目标图像,选择概率值大于经验阈值的主题,识别航天器目标的类别。
采用上述进一步方案的有益效果是信息优化过程减轻计算量,有助于提高处理速度;根据图像中局部显著特征的特性,不同型号航天器的视觉特征和其语义概率值,依据后验概率值和经验阈值来判断语义信息,能够最大程度的保证视觉类似的航天器目标的识别。
进一步,所述步骤310包括:
步骤311,利用K近邻算法求取每个所述特征点的k个近邻特征点,采用Dijkstra距离进行度量来保持各个特征点之间的几何特性;
步骤312,根据下式计算特征信息的局部重建权值矩阵:
其中,xi表示第个特征点,xij为特征点xi的满足Dijkstra距离的第个近邻点,Wi是用特征点xi的近邻点对其进行线性重构时的表示坐标;
采用拉格朗日乘子法求局部最优化重建权值矩阵:
其中,Qi表示依据特征点xi和特征点xj构造的Gram矩阵;
步骤313,在条件下将所有高维特征数据映射到低维空间,获得优化后的特征信息;式中yi表示在低维空间中重构的特征点。
采用上述进一步方案的有益效果是所有高维特征数据映射到低维空间能够有效减轻计算量,有助于提高处理速度。
进一步,所述步骤320,包括:
对优化后的特征信息进行K-means聚类得到k个聚类中心,每个聚类中心对应一个视觉单词w;k个所述视觉单词w构成视觉单词表W;
基于pLSA语义模型,视觉单词w、增强后的航天器观测图像及潜在语义主题z所对应的联合概率分布为:而不包含未观测到的潜在变量z时的概率分布为:
根据联合概率分布关系由上式可得的概率表示为:进而可以获得的概率值。
采用上述进一步方案的有益效果是采用K-means聚类算法对提取的特征信息进行聚类,速度快,效率高,时间复杂度近于线性;基于pLSA语义模型,获取视觉单词w、增强后的航天器观测图像及潜在语义主题z所对应的联合概率分布,进而求取潜在语义主题z的概率值,利用语义识别能够最大程度的保证视觉类似的航天器目标的识别。
本发明另一方面提供一种航天器目标快速识别装置,包括:
图像预处理模块,用于对获取的航天器观测图像的进行杂光抑制及增强处理;
特征提取模块,用于对处理后的航天器观测图像进行特征点及特征信息提取;所述特征信息包括所述特征点的局部灰度特征信息提取和不变矩特征信息;
识别模块,用于对所述特征信息进行优化处理,并通过语义分析和判别进行航天器目标的识别。
本发明的有益效果是:航天器在空中被动成像,由于光照原因会形成杂光,使得局部特征信息受到干扰,导致后续特征提取不准确,影响识别效果。本装置通过杂光抑制及增强处理,能提高图像信噪比,增强图像质量,利于后续目标快速识别。对观测图像的特征点以及包含特征点的不变矩特征和局部灰度特征的特征信息的提取,能更有利于鉴别航天器型号。信息优化过程减轻计算量,有助于提高处理速度。根据图像中局部显著特征的特性,不同型号航天器的视觉特征和其语义概率值,然后依据后验概率值和经验阈值来判断语义信息,能够最大程度的保证视觉类似的航天器目标图像的目标识别。
附图说明
图1为杂光形成示意图;
图2为本发明实施例提供的一种航天器目标快速识别方法流程图;
图3为阴影系数的几何示意图;
图4吉文斯旋转示意图;
图5为由四个区域组成的X型特征点示意图;
图6为感兴趣区域内灰度梯度投影示意图;
图7为本发明实施例提供的一种航天器目标快速识别装置结构示意图。
具体实施方式
以下结合实例对本发明的原理和特征进行描述,所举实例只用于解释本发明,并非用于限定本发明的范围。
航天器在空间中成像会受到环境杂光的干扰杂光形成原理如图1所示。
图2为本发明实施例提供的一种航天器目标快速识别方法流程图;如图2所示,本发明提供一种航天器目标快速识别方法,包括,
步骤1,对获取的航天器观测图像进行杂光抑制增强处理,得到所述航天器观测图像的恢复图像;
步骤2,对所述航天器观测图像的恢复图像进行特征点及特征信息的提取;所述特征信息包括所述特征点的不变矩特征信息和局部灰度特征信息;
步骤3,对所述特征信息进行优化处理,并通过语义分析和判别进行航天器目标的识别。
该实施例中,航天器在空中被动成像,由于光照原因会形成杂光,使得局部特征信息受到干扰,导致后续特征提取不准确,影响识别效果。本方法通过杂光抑制及增强处理,能提高图像信噪比,增强图像质量,利于后续目标快速识别;对观测图像的特征点以及包含特征点的不变矩特征和局部灰度特征的特征信息的提取,能更有利于鉴别航天器型号;信息优化过程减轻计算量,有助于提高处理速度。根据图像中局部显著特征的特性,不同型号航天器的视觉特征和其语义概率值,然后依据后验概率值和经验阈值来判断语义信息,能够最大程度的保证视觉类似的航天器目标的识别。
可选地,作为本发明一个实施例,步骤1,对获取的航天器观测图像的进行杂光抑制及增强处理,获得所述航天器观测图像的恢复图像包括以下内容:
步骤101,根据式(1-1)对航天器观测图像进行基于Van Cittert迭代法的图像恢复,
式(1-1)中y为观测图像,为在第t次迭代后得到的观测图像,A是杂散光降质图像的点扩散函数;
根据场景图像成像时杂散光的不同成因,可将杂散光模型表示为:A=((1-β)GDA+βSSTR),其中SSTR是纯杂散光部分的点扩散函数矩阵,β为其权重系数,GDA为光学成像系统自身所引起的衍射和色差部分对图像的干扰所形成的杂散光;
对于一幅受杂散光干扰的观测图像而言,图像中某一像素点(ip,jp)的衍射和色差部分的点扩散函数模型为:
其中,K是归一化因子,能够保证整个杂散光空变点扩散函数模型在二维的x-y平面上g(iq,jq,ip,jp;σ)的和为1,且有助于增加图像中细节内容的高频信息,同时使得图像中总的能量得到保持。在光学系统的成像过程中,来自外部杂光对于图像造成干扰部分的纯杂散光部分的点扩散函数模型表达式为:
其中,z为归一化因子,其分析表达式为:使得在图像中都满足同时,有研究结果表明,在视角很广的情况下,非正常光路进入的光线其衰减符合余弦四次方定律,根据图3中相机在焦平面上成像的几何关系示意图可以得到cos2γ的计算方法为:D代表相机出瞳到图像平面的距离,而γ表示光轴和成像系统的夹角,可将cos2γ作为杂散光空变点扩散函数模型的阴影系数来模拟成像时照度的衰减。
设在相机光学系统成像过程中,增加阴影系数后的杂散光空变点扩散模型S表达式可写为:
其中,衍射和色差部分的点扩散函数GDA在实际情况中可用脉冲响应函数来对其近似表达,所以可设GDA=I(I是单位阵)。则(1-1)可改写为:
由于在初始迭代情况下原始清晰图像受到移变杂散光点扩散函数的影响最小,即因此式(1-2)可简化为:
在公式(1-3)中,场景图像成像时杂散光的点扩散函数模型不仅是全局的,而且其本质也是空变的,模型中相关的参数为(a,b,c,α,β,σ),利用在杂散光切面方向构造函数来对模型的参数进行估计。对于一幅受杂散光影响而降质的场景图像,将其作为观测图像,利用(1-3)式进行图像恢复得到最终的清晰图像,而杂散光的点对点移变点扩散函数SSTR与观测图像y的卷积计算在图像恢复算法中计算复杂度和存贮代价较高,如何对其进行加速优化是杂散光降质图像校正算法实现和应用的关键问题。
将(1-3)式中点对点杂散光空变点扩散函数SSTR与观测图像y卷积的部分表示为:
X=SSTRy
在数字图像中,将SSTR量化表示为[S],用δSSTR表示由量化引起的误差,则SSTR=[S]+δSSTR,同时由于SSTR量化误差造成的X误差可表示为:
δX=δSSTRy
在有损信源编码技术中,编码质量的优劣和度量其形变的方法密切相关。由于Frobenius范数的优良特性,本文采用Frobenius误差来度量由有损信源编码机制导致的形变,可记为||δX||2,在实际算法中,可以根据需求来调整量化步长。假设观测到的杂散光降质图像y和量化误差δX彼此之间相互独立,则SSTR的量化导致的X的形变可表示为:
其中,Ry=E[yyt],如果观测图像y是白的,Ry则必定是一个单位阵。
步骤102,对y进行分解并使其满足约束条件Ry=EΛyEt,E是特征向量的正交阵,而Λy是关于特征值的对角阵。由上述可得:
对于估计图像矩阵发生形变时,其均方差的最小值等于形变误差的期望最小值。由此,矩阵SSTR可近似表示为有效最小化同时在最小均方误差(Normalized RootMean Squared Error,NRMSE)条件下进行图像恢复;
步骤103,采用KL变换和小波变换对矩阵的行和列进行去相关运算,其表达式如下:
通过相关行列变换后,只有小部分非零元素的矩阵能量被集中到矩阵中,基于小波变换的算法在图像压缩编码中已得到稳健的应用,因此可采用2D的正交小波变换Wt为的列进行变换。同时设ψ是协方差矩阵的特征向量矩阵,M、N为点扩散函数矩阵SSTR的行与列个数,采用ψ来对去的行相关性,便于图像的恢复重建。同时,选择ψ作为特征向量正交阵,ΛS作为SSTR的特征向量对角阵,则有:
RS=ψΛSψt
由上述可得杂散光点扩散函数与观测图像可分别表示为如下形式:
其中,而杂散光空变点扩散函数SSTR与观测图像y的卷积可近似表示为:
式中的稀疏表示已经过量化和编码阶段,而和观测图像y紧密相关的变换T依然不合符稀疏表示条件,需要作进一步处理。类似KL变换或者小波变换的方法已经被用于矩阵T的变换中,由于杂散光的空变点扩散函数已经被稀疏表示,当T也满足稀疏表示条件时,其计算速度就会显著提高。
步骤104,稀疏矩阵变换
稀疏矩阵是指大多数元素为零元素,只有极少非零元素的矩阵。和传统的变换方式不同,稀疏矩阵变换(Sparse Matrix Transform,SMT)只用其有限数目点坐标的旋转,即吉文斯旋转的乘积就可以实现对既定目标矩阵的精确估计,且其运算和存储只需要针对非零元即可。
对于式(1-5)中的采用K对非零元素坐标的旋转的乘积对其进行矩阵行去相关处理,对观测图像y的变换T采用K对坐标的吉文斯旋转的乘积以代替精确的变换近似得到;稀疏矩阵变换T的K对坐标的吉文斯旋转形式表示为:
其中,每个矩阵Ti(i=K-1,K-2,…,0)为对应坐标点的吉文斯旋转;每次只需执行不超过4次乘法,图3给出了点坐标为(xk,yk)时吉文斯旋转变换的示意图。实际应用中,目标变换矩阵中非零元素K的数目远小于矩阵规模,即满足假设K<<MN。对于观测图像而言,采用稀疏矩阵变换对其进行稀疏化表达将加快计算速度,且可获得与精确矩阵变换近似的表达。
根据吉文斯旋转的定义,Ti表示为:
Ti=BiΛiAi (1-7)
其中,Ai和Bi为吉文斯旋转,Λi是对角标准化矩阵。每个吉文斯旋转可由坐标点在其自身张成的平面中正交旋转得到,有限数目个吉文斯旋转即可获得有效的稀疏矩阵变换。
为了降低计算复杂度,式(1-7)中可设Λi=I,则Ai和Bi可由非零元素点(xi,yj)以及其旋转角度φk和θk表示如下:
Ai=I+Y(xi,yj,φk) (1-8)
Bi=I+Y(xi,yj,θk) (1-9)
其中,函数Y(m,n,θ)的取值可作定义如下:
其中,对应于不同的坐标点及旋转角度,由函数Y(m,n,θ)可获得矩阵中非零元素Ai和Bi相对应的取值,每次的乘法运算,不超过4次。作为估计正交变换的有效工具,吉文斯旋转可以有效降低反矩阵求解问题的计算复杂度。吉文斯旋转的个数仅限于目标矩阵中有限个非零元素,因此对应的稀疏变换矩阵只采用较少的简单矩阵乘法来近似,且能保证算法的有效性。实验中,吉文斯旋转的个数K等于目标矩阵中非零元素的个数,为了简化计算设旋转的角度为45°。
步骤105,将稀疏矩阵变换T的K个分量及其对应位置归位,其他位置置零,进行解码得到代入式(1-3),最终得到
可选地,作为本发明一个实施例,步骤2,对所述航天器观测图像的恢复图像进行特征点及特征信息的提取,包括以下内容:
步骤201,候选特征点提取
采用预设大小n*n的窗口(优选为5*5大小)对每个像元逐一探测,以所述像元为中心,分别计算所述像元四个方向(0°、45°、90°、135°)和相邻像素的灰度差的平方和,
其中k=[n/2]。
取四个所述平方和V1,V2,V3,V4中的最小值作为该像元的兴趣值,同时计算四个所述平方和的中值:将所述兴趣值大于中值的像元作为候选特征点。
步骤202,特征点提取
利用SUSAN算子,对所述候选特征点进行区域分割,并对分割获得的候选特征点区域进行分析,去掉伪特征点信息,获得所述航天器观测图像的特征点;
图4为由ABCD四个区域构成的“X”连接型候选特征点,选取构成候选特征点的区域中的最大的区域作为候选特征点的第一区域,而将其他区域视为单一的背景区域;假设整个候选特征点区域足够大,能把这个区域分为许多同心的圆,第一区域的特征点位置为圆中心;同样地,使Si是第i个同心圆上第一区域的对应圆弧,对于每一个圆Ci,则有:
用计算信息量的方法获得所述候选特征点响应函数,将其定义为模板内像元n与平均信息量差的平方和,以模板内像元幅度值的累加值除以模板面积作为该模板的平均信息量gr(n),其计算方法如下:
其中,fk,l(n)表示内坐标(k,l)点的幅值,r表示圆形模板的半径,Ω表示模板的覆盖范围,特征点响应函数为:
选取函数值最大的点作为特征点,其余的候选特征点作为伪特征点被除去。
步骤203,局部灰度特征信息提取
在步骤202的基础上,将特征点周围20*20的区域作为目标图像的感兴趣区域,采用Sobel算子计算感兴趣区域内每个像素点的灰度梯度,水平和垂直两个方向上的灰度梯度分量分别为:
gx(x,y)=f(x+1,y-1)+2f(x+1,y)+f(x+1,y+1)-f(x-1,y-1)-2f(x-1,y)-f(x-1,y+1)
gy(x,y)=f(x-1,y+1)+2f(x,y+1)+f(x+1,y+1)-f(x-1,y-1)-2f(x,y-1)-f(x+1,y-1)
则区域内每个像素点灰度梯度向量g=gxx+gyy的幅值和幅角如下:
θ=arcsin(gy(x,y)/||g||),θ∈[0,2π]
其中,幅角θ在取值范围[0,2π]中被均匀量化为多个(优选为8个)区间,并按照量化后各个幅角所属区间的不同,将灰度梯度的幅值分解与所述多个区间对应的8个方向的属性平面;最后判定梯度向量幅角所属的范围,并将梯度向量按平行四边形法则投影到如图6所示的nθ0和(n+1)θ0两个方向上,记录投影长度gn和gn+1:
同时采用高斯滤波器对以特征点为中心的感兴趣区域中的每个方向属性平面进行亚采样,增强局部灰度特征信息的鲁棒性,高斯滤波器冲激响应函数为:
其中,高斯滤波器的尺度参数tx是采样的间隔,高斯滤波函数为:
其中,(x0,y0)为以特征点为中心的感兴趣区域中的像素点。
以特征点为中心的20*20的感兴趣区域以间隔tx=4进行二维亚采样,即可得到25个采样值,并分别在8个方向属性平面上将得到的采样值进行合并就可得到200维的灰度梯度特征信息。
步骤204,不变矩特征信息提取
对于大小为M*N的图像f(x,y),其(p+q)阶矩定义为:
(p+q)阶中心矩定义为:
其中和是图像的矩心,中心矩心为:
当航天器目标的观测图像发生变化时,mpq也发生变化,而μpq则具有平移不变性,归一化中心矩为:
在步骤202的基础上,以特征点为中心,依据二阶和三阶的归一化中心矩,构造以下7个不变矩:
φ1=y20+y02
φ3=(y30+3y12)2+(3y21-y03)2
φ4=(y30+y12)2+(y21+y03)2
φ5=(y30-3y12)(y30+y12)[(y30+y12)2-3(y21+y03)2]+(3y21-y03)(y21+y03)[3(y30+y12)2-(y21+y03)2]
φ6=(y20-y02)[(y30+y12)2-(y21+y30)2]+4y11(y30+y21)(y21+y03)
φ7=(3y21+y03)(y30+y12)[(y30+y12)2-3(y21+y03)2]+(y30-3y12)(y21+y30)[3(y30+y12)2-(y21+y03)2]
这7个不变矩构成一组特征向量,根据待识别分析的空间目标图像的实际情况,对不变矩特征变量φi取绝对值的对数调整其取值范围为φi=lg|φi|,i=1,2,…,7。由于不变矩需要对目标区域中的每一个像素进行运算,为了减小计算量,仅选用对分类影响最大的4个分量φ1,φ2,φ3和φ4,作为目标的矩特征。
可选地,作为本发明一个实施例,步骤3,对所述特征信息进行优化处理,并通过语义分析和判别进行航天器目标的识别,包括以下内容:
步骤310,对所述特征信息进行优化,将所有高维特征数据映射到低维空间,获得优化后的特征信息;
具体包括:
步骤311,利用K近邻算法求取每个所述特征点的k个近邻特征点,采用Dijkstra距离进行度量来保持各个特征点之间的几何特性;
步骤312,根据下式计算特征信息的局部重建权值矩阵:
其中,xi表示第个特征点,xij为特征点xi的满足Dijkstra距离的第个近邻点,Wi是用特征点xi的近邻点对其进行线性重构时的表示坐标;
采用拉格朗日乘子法求局部最优化重建权值矩阵:
其中,Qi表示依据特征点xi和特征点xj构造的Gram矩阵;
步骤313,在条件下将所有高维特征数据映射到低维空间,获得优化后的特征信息;式中yi表示在低维空间中重构的特征点。
步骤320,对优化后的特征信息进行K-means聚类得到k个聚类中心,每个聚类中心对应一个视觉单词,获取视觉单词w、所述航天器观测图像的恢复图像及其潜在语义主题z所对应的联合概率分布,进而根据所述联合概率分布计算概率值
具体的,对提取的信息进行K-means聚类得到k个聚类中心,每个聚类中心对应一个视觉单词w;k个所述视觉单词w构成视觉单词表W;
基于pLSA语义模型,视觉单词w、增强后的航天器观测图像及潜在语义主题z所对应的联合概率分布为:而不包含未观测到的潜在变量z时的概率分布为:
根据联合概率分布关系由上式可得的概率表示为:进而可以获得的概率值。
步骤330,采用EM学习方法获取参数阈值,根据参数阈值来判别所述航天器观测图像的恢复图像所属的语义信息;根据不同类型航天器目标图像,选择概率值大于经验阈值的主题,识别航天器目标的类别。
本发明另一方面提供一种航天器目标快速识别装置,如图7所示,包括:
图像预处理模块,用于对获取的航天器观测图像的进行杂光抑制及增强处理;
特征提取模块,用于对处理后的航天器观测图像进行特征点及特征信息提取;所述特征信息包括所述特征点的局部灰度特征信息提取和不变矩特征信息;
识别模块,用于对所述特征信息进行优化处理,并通过语义分析和判别进行航天器目标的识别。
航天器在空中被动成像,由于光照原因会形成杂光,使得局部特征信息受到干扰,导致后续特征提取不准确,影响识别效果。本装置通过杂光抑制及增强处理,能提高图像信噪比,增强图像质量,利于后续目标快速识别。对观测图像的特征点以及包含特征点的不变矩特征和局部灰度特征的特征信息的提取,能更有利于鉴别航天器型号。信息优化过程减轻计算量,有助于提高处理速度。根据图像中局部显著特征的特性,不同型号航天器的视觉特征和其语义概率值,然后依据后验概率值和经验阈值来判断语义信息,能够最大程度的保证视觉类似的航天器目标图像的目标识别。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种航天器目标快速识别方法,其特征在于:包括,
步骤1,对获取的航天器观测图像进行杂光抑制增强处理,得到所述航天器观测图像的恢复图像;
步骤2,对所述航天器观测图像的恢复图像进行特征点及特征信息的提取;所述特征信息包括所述特征点的不变矩特征信息和局部灰度特征信息;
步骤3,对所述特征信息进行优化处理,并通过语义分析和判别进行航天器目标的识别。
2.根据权利要求1所述一种航天器目标快速识别方法,其特征在于,所述步骤1包括:
步骤101,根据下式对获取的所述航天器观测图像进行图像恢复,
<mrow>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mo>=</mo>
<mn>2</mn>
<mi>y</mi>
<mo>-</mo>
<msup>
<mi>cos</mi>
<mn>2</mn>
</msup>
<mi>&gamma;</mi>
<mrow>
<mo>(</mo>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mi>&beta;</mi>
</mrow>
<mo>)</mo>
<mo>+</mo>
<msub>
<mi>&beta;S</mi>
<mrow>
<mi>S</mi>
<mi>T</mi>
<mi>R</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mi>y</mi>
</mrow>
式中,y为所述航天器观测图像,为恢复图像,cos2γ是阴影系数,γ表示光轴和成像系统的夹角,SSTR是纯杂散光部分的点扩散函数,β为其权重系数;
步骤102,对纯杂散光部分的点扩散函数SSTR进行量化,并采用Frobenius误差度量量化形变引起的误差,然后根据下式得到矩阵
<mrow>
<mover>
<mi>S</mi>
<mo>~</mo>
</mover>
<mo>=</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>S</mi>
<mi>T</mi>
<mi>R</mi>
</mrow>
</msub>
<msubsup>
<mi>E&Lambda;</mi>
<mi>y</mi>
<mrow>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msubsup>
</mrow>
式中,E为单位正交阵,Λy为对航天器观测图像进行特征分解得到的对角阵;
步骤103,采用K-L变换和小波变换对所述矩阵的行和列进行去相关运算;
步骤104,对去相关后的矩阵进行量化,并利用稀疏矩阵变换对量化后的矩阵的行进行去相关,通过下式得到经过杂光抑制增强处理的所述航天器观测图像的恢复图像
<mrow>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mo>=</mo>
<mn>2</mn>
<mi>y</mi>
<mo>-</mo>
<msup>
<mi>cos</mi>
<mn>2</mn>
</msup>
<mi>&gamma;</mi>
<mrow>
<mo>(</mo>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mi>&beta;</mi>
</mrow>
<mo>)</mo>
<mo>+</mo>
<mi>&beta;</mi>
<mi>W</mi>
<mo>&lsqb;</mo>
<mover>
<mi>S</mi>
<mo>&OverBar;</mo>
</mover>
<mo>&rsqb;</mo>
<msub>
<mi>T</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mi>y</mi>
</mrow>
式中,W为二维正交小波变换矩阵,为量化后的空变点扩散函数,Tk为根据吉文斯旋转定理得到的稀疏矩阵。
3.根据权利要求1所述一种航天器目标快速识别方法,其特征在于,所述步骤2包括:
步骤201,对所述航天器观测图像的恢复图像进行候选特征点提取;
步骤202,利用SUSAN算子,对所述候选特征点进行区域分割,并对分割获得的候选特征点区域进行分析,去掉伪特征点,获得所述航天器观测图像的特征点;
步骤203,将所述特征点周围预设区域作为感兴趣区域,采用Sobe l算子提取感兴趣区域的灰度特征信息,即为所述特征点的灰度特征信息;
步骤204,以所述特征点为中心,求取所述恢复图像的归一化中心矩,并依据二阶和三阶的归一化中心矩,构造多个不变矩,所述多个不变矩构成一组特征量,所述特征量即为所述特征点的不变矩特征信息。
4.根据权利要求3所述一种航天器目标快速识别方法,其特征在于,所述步骤201包括:
采用预设大小的窗口对每个像元逐一探测,分别以每个所述像元为中心,分别计算每个所述像元四个方向和相邻像素的灰度差的平方和,取四个所述平方和中的最小值作为该像元的兴趣值,同时计算四个所述平方和的中值,将所述兴趣值大于所述中值的像元作为侯选特征点。
5.根据权利要求3所述一种航天器目标快速识别方法,其特征在于,所述步骤202包括:
以候选特征点所在的圆形区域作为SUSAN算子的模板,选取所述模板中构成候选特征点的最大的子域作为所述候选特征点的第一区域,而将其他子域视为单一的背景区域;将所述模板分为多个同心的圆,第一区域的特征点位置为圆中心;定义Si是第i个同心圆上第一区域的对应圆弧,对于每一个圆,则有:
用计算信息量的方法获得候选特征点响应函数,将其定义为模板内像元n与平均信息量差的平方和,以模板内像元幅度值的累加值除以模板面积作为该模板的平均信息量gr(n),其计算方法如下:
<mrow>
<msub>
<mi>g</mi>
<mi>r</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msup>
<mi>&pi;r</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<munder>
<mo>&Sigma;</mo>
<mi>&Omega;</mi>
</munder>
<msub>
<mi>f</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>l</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>,</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
</mrow>
其中,fk,l(n)表示内坐标(k,l)点的幅值,r表示圆形模板的半径,Ω表示模板的覆盖范围,特征点响应函数为:
<mrow>
<msub>
<mi>w</mi>
<mi>r</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munder>
<mo>&Sigma;</mo>
<mi>&Omega;</mi>
</munder>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>l</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>g</mi>
<mi>r</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
选取函数值最大的点作为特征点,其余的候选特征点作为伪特征点被除去。
6.根据权利要求3所述一种航天器目标快速识别方法,其特征在于,所述步骤203包括:
将特征点周围预设区域作为目标图像的感兴趣区域,采用Sobe l算子计算感兴趣区域内每个像素点水平和垂直两个方向上的灰度梯度分量,获取感兴趣区域内每个像素点灰度梯度向量的幅值和幅角,即所述特征点的灰度特征信息。
7.根据权利要求1所述一种航天器目标快速识别方法,其特征在于,所述步骤3包括:
步骤310,对所述特征信息进行优化,将所有高维特征数据映射到低维空间,获得优化后的特征信息;
步骤320,对优化后的特征信息进行K-means聚类得到k个聚类中心,每个聚类中心对应一个视觉单词,获取视觉单词w、所述航天器观测图像的恢复图像及其潜在语义主题z所对应的联合概率分布,进而根据所述联合概率分布计算概率值
步骤330,采用EM学习方法获取参数阈值,根据参数阈值来判别所述航天器观测图像的恢复图像所属的语义信息;根据不同类型航天器目标图像,选择概率值大于经验阈值的主题,识别航天器目标的类别。
8.根据权利要求7所述一种航天器目标快速识别方法,其特征在于,所述步骤310,包括:
步骤311,利用K近邻算法求取每个所述特征点的k个近邻特征点,采用D i jkstra距离进行度量来保持各个特征点之间的几何特性;
步骤312,根据下式计算所述特征信息的局部重建权值矩阵:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&epsiv;</mi>
<mrow>
<mo>(</mo>
<mi>W</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<mo>|</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<msubsup>
<mi>w</mi>
<mi>j</mi>
<mi>i</mi>
</msubsup>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msubsup>
<mi>W</mi>
<mi>j</mi>
<mi>i</mi>
</msubsup>
<mo>=</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,xi表示第个特征点,xij为特征点xi的满足Di jkstra距离的第个近邻点,Wi是用特征点xi的近邻点对其进行线性重构时的表示坐标;
采用拉格朗日乘子法求局部最优化重建权值矩阵:
<mrow>
<msubsup>
<mi>w</mi>
<mi>j</mi>
<mi>i</mi>
</msubsup>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<msubsup>
<mrow>
<mo>(</mo>
<msup>
<mi>Q</mi>
<mi>i</mi>
</msup>
<mo>)</mo>
</mrow>
<mrow>
<mi>j</mi>
<mi>m</mi>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>p</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>q</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<msubsup>
<mrow>
<mo>(</mo>
<msup>
<mi>Q</mi>
<mi>i</mi>
</msup>
<mo>)</mo>
</mrow>
<mrow>
<mi>p</mi>
<mi>q</mi>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
</mrow>
</mfrac>
</mrow>
其中,Qi表示依据特征点xi和特征点xj构造的Gram矩阵;
步骤313,在条件下,将所有高维特征数据映射到低维空间,获得优化后的特征信息;式中yi表示在低维空间中重构的特征点。
9.根据权利要求7所述一种航天器目标快速识别方法,其特征在于,所述步骤320,包括:
对优化后的特征信息进行K-means聚类得到k个聚类中心,每个聚类中心对应一个视觉单词w;k个所述视觉单词w构成视觉单词表W;
基于pLSA语义模型,视觉单词w、所述航天器观测图像的恢复图像及潜在语义主题z所对应的联合概率分布为:而不包含未观测到的潜在变量z时的概率分布为:
<mrow>
<mi>P</mi>
<mrow>
<mo>(</mo>
<mi>w</mi>
<mo>,</mo>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>z</mi>
<mo>&Element;</mo>
<mi>Z</mi>
</mrow>
</munder>
<mi>P</mi>
<mrow>
<mo>(</mo>
<mi>w</mi>
<mo>,</mo>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mo>,</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>P</mi>
<mrow>
<mo>(</mo>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mo>)</mo>
</mrow>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>z</mi>
<mo>&Element;</mo>
<mi>Z</mi>
</mrow>
</munder>
<mi>P</mi>
<mrow>
<mo>(</mo>
<mi>w</mi>
<mo>|</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mi>P</mi>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>|</mo>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mo>)</mo>
</mrow>
</mrow>
根据联合概率分布关系由上式可得的概率表示为:进而可以获得的概率值。
10.一种航天器目标快速识别装置,其特征在于,包括:
图像预处理模块,用于对获取的航天器观测图像的进行杂光抑制及增强处理;
特征提取模块,用于对处理后的航天器观测图像进行特征点及特征信息提取;所述特征信息包括所述特征点的局部灰度特征信息提取和不变矩特征信息;
识别模块,用于对所述特征信息进行优化处理,并通过语义分析和判别进行航天器目标的识别。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711083754.8A CN108052957B (zh) | 2017-11-07 | 2017-11-07 | 一种航天器目标快速识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711083754.8A CN108052957B (zh) | 2017-11-07 | 2017-11-07 | 一种航天器目标快速识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108052957A true CN108052957A (zh) | 2018-05-18 |
CN108052957B CN108052957B (zh) | 2021-09-14 |
Family
ID=62119022
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711083754.8A Active CN108052957B (zh) | 2017-11-07 | 2017-11-07 | 一种航天器目标快速识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108052957B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109115229A (zh) * | 2018-09-17 | 2019-01-01 | 中国人民解放军国防科技大学 | 利用低频测姿传感器对航天器高频姿态进行测量的方法 |
CN109166122A (zh) * | 2018-09-18 | 2019-01-08 | 广西电网有限责任公司电力科学研究院 | 基于图像处理的断路器操作机构弹簧伸缩特性测试方法 |
CN109631912A (zh) * | 2019-01-10 | 2019-04-16 | 中国科学院光电技术研究所 | 一种深空球形目标被动测距方法 |
CN109816650A (zh) * | 2019-01-24 | 2019-05-28 | 强联智创(北京)科技有限公司 | 一种基于二维dsa图像的目标区域识别方法及其系统 |
CN111563517A (zh) * | 2020-04-20 | 2020-08-21 | 腾讯科技(深圳)有限公司 | 图像处理方法、装置、电子设备及存储介质 |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101540043A (zh) * | 2009-04-30 | 2009-09-23 | 南京理工大学 | 单幅图像复原的解析迭代快速频谱外推方法 |
CN101733746A (zh) * | 2009-12-22 | 2010-06-16 | 哈尔滨工业大学 | 空间机器人非合作目标自主识别与捕获方法 |
CN103020904A (zh) * | 2011-09-23 | 2013-04-03 | 波音公司 | 反射消除系统 |
CN103049889A (zh) * | 2012-12-06 | 2013-04-17 | 华中科技大学 | 一种强噪声气动光学效应退化图像的预处理方法 |
CN103473551A (zh) * | 2013-09-16 | 2013-12-25 | 中国传媒大学 | 基于sift算子的台标识别方法及系统 |
CN103617328A (zh) * | 2013-12-08 | 2014-03-05 | 中国科学院光电技术研究所 | 一种飞机三维姿态解算方法 |
CN103886559A (zh) * | 2014-04-04 | 2014-06-25 | 北京航空航天大学 | 一种光谱图像处理方法 |
CN103996027A (zh) * | 2014-05-19 | 2014-08-20 | 上海微小卫星工程中心 | 一种天基空间目标识别方法 |
US20160171727A1 (en) * | 2014-12-16 | 2016-06-16 | The Regents Of The University Of California | Feature-preserving noise removal |
US20160228000A1 (en) * | 2015-02-06 | 2016-08-11 | Richard F. Spaide | Volume analysis and display of information in optical coherence tomography angiography |
US20160291343A1 (en) * | 2013-03-11 | 2016-10-06 | Sudhakar Prasad | Rotating point-spread function (psf) design for three-dimensional imaging |
CN106650629A (zh) * | 2016-11-21 | 2017-05-10 | 河海大学 | 一种基于核稀疏表示的快速遥感目标检测识别方法 |
CN106679634A (zh) * | 2016-06-20 | 2017-05-17 | 山东航天电子技术研究所 | 一种基于立体视觉的空间非合作目标位姿测量方法 |
-
2017
- 2017-11-07 CN CN201711083754.8A patent/CN108052957B/zh active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101540043A (zh) * | 2009-04-30 | 2009-09-23 | 南京理工大学 | 单幅图像复原的解析迭代快速频谱外推方法 |
CN101733746A (zh) * | 2009-12-22 | 2010-06-16 | 哈尔滨工业大学 | 空间机器人非合作目标自主识别与捕获方法 |
CN103020904A (zh) * | 2011-09-23 | 2013-04-03 | 波音公司 | 反射消除系统 |
CN103049889A (zh) * | 2012-12-06 | 2013-04-17 | 华中科技大学 | 一种强噪声气动光学效应退化图像的预处理方法 |
US20160291343A1 (en) * | 2013-03-11 | 2016-10-06 | Sudhakar Prasad | Rotating point-spread function (psf) design for three-dimensional imaging |
CN103473551A (zh) * | 2013-09-16 | 2013-12-25 | 中国传媒大学 | 基于sift算子的台标识别方法及系统 |
CN103617328A (zh) * | 2013-12-08 | 2014-03-05 | 中国科学院光电技术研究所 | 一种飞机三维姿态解算方法 |
CN103886559A (zh) * | 2014-04-04 | 2014-06-25 | 北京航空航天大学 | 一种光谱图像处理方法 |
CN103996027A (zh) * | 2014-05-19 | 2014-08-20 | 上海微小卫星工程中心 | 一种天基空间目标识别方法 |
US20160171727A1 (en) * | 2014-12-16 | 2016-06-16 | The Regents Of The University Of California | Feature-preserving noise removal |
US20160228000A1 (en) * | 2015-02-06 | 2016-08-11 | Richard F. Spaide | Volume analysis and display of information in optical coherence tomography angiography |
CN106679634A (zh) * | 2016-06-20 | 2017-05-17 | 山东航天电子技术研究所 | 一种基于立体视觉的空间非合作目标位姿测量方法 |
CN106650629A (zh) * | 2016-11-21 | 2017-05-10 | 河海大学 | 一种基于核稀疏表示的快速遥感目标检测识别方法 |
Non-Patent Citations (5)
Title |
---|
C. KWAN等: "Mastcam image enhancement using estimated point spread functions", 《2017 IEEE 8TH ANNUAL UBIQUITOUS COMPUTING, ELECTRONICS AND MOBILE COMMUNICATION CONFERENCE (UEMCON)》 * |
DONG-HAN LEE等: "image restoration of the asymmetric point spread function of a high-resolution remote sensing satellite with time-delayed integration", 《ADVANCES IN SPACE RESEARCH》 * |
崔毅: "推扫式宽视场CCD成像光谱技术研究", 《中国博士学位论文全文数据库 工程科技II辑》 * |
徐贵力等: "基于光照模糊相似融合不变矩的航天器目标识别", 《航空学报》 * |
曹雷: "低对比度目标探测跟踪技术研究", 《中国博士学位论文全文数据库 信息科技辑》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109115229A (zh) * | 2018-09-17 | 2019-01-01 | 中国人民解放军国防科技大学 | 利用低频测姿传感器对航天器高频姿态进行测量的方法 |
CN109115229B (zh) * | 2018-09-17 | 2019-05-31 | 中国人民解放军国防科技大学 | 利用低频测姿传感器对航天器高频姿态进行测量的方法 |
CN109166122A (zh) * | 2018-09-18 | 2019-01-08 | 广西电网有限责任公司电力科学研究院 | 基于图像处理的断路器操作机构弹簧伸缩特性测试方法 |
CN109631912A (zh) * | 2019-01-10 | 2019-04-16 | 中国科学院光电技术研究所 | 一种深空球形目标被动测距方法 |
CN109631912B (zh) * | 2019-01-10 | 2022-08-23 | 中国科学院光电技术研究所 | 一种深空球形目标被动测距方法 |
CN109816650A (zh) * | 2019-01-24 | 2019-05-28 | 强联智创(北京)科技有限公司 | 一种基于二维dsa图像的目标区域识别方法及其系统 |
CN111563517A (zh) * | 2020-04-20 | 2020-08-21 | 腾讯科技(深圳)有限公司 | 图像处理方法、装置、电子设备及存储介质 |
CN111563517B (zh) * | 2020-04-20 | 2023-07-04 | 腾讯科技(深圳)有限公司 | 图像处理方法、装置、电子设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN108052957B (zh) | 2021-09-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108052957B (zh) | 一种航天器目标快速识别方法 | |
Wei et al. | Spatiotemporal fusion of MODIS and Landsat-7 reflectance images via compressed sensing | |
CN108182449A (zh) | 一种高光谱图像分类方法 | |
Nascimento et al. | Comparing edge detection methods based on stochastic entropies and distances for PolSAR imagery | |
CN103632155B (zh) | 基于慢特征分析的遥感影像变化检测方法 | |
CN109766858A (zh) | 结合双边滤波的三维卷积神经网络高光谱影像分类方法 | |
CN101526995B (zh) | 基于对角子类判决分析的合成孔径雷达目标识别方法 | |
JP5063776B2 (ja) | 幾何学変換に基づく一般化統計的テンプレートマッチング | |
Uss et al. | Multimodal remote sensing image registration with accuracy estimation at local and global scales | |
Paillassa et al. | MAXIMASK and MAXITRACK: Two new tools for identifying contaminants in astronomical images using convolutional neural networks | |
Bakken et al. | The effect of dimensionality reduction on signature-based target detection for hyperspectral remote sensing | |
Gelvez-Barrera et al. | Joint nonlocal, spectral, and similarity low-rank priors for hyperspectral–multispectral image fusion | |
Su et al. | Superpixel-based weighted collaborative sparse regression and reweighted low-rank representation for hyperspectral image unmixing | |
CN112036235A (zh) | 一种高光谱图像目标检测方法和系统 | |
CN111105363A (zh) | 一种含噪高光谱图像快速解混方法 | |
CN111199251B (zh) | 一种基于加权邻域的多尺度高光谱图像分类方法 | |
CN104463245A (zh) | 一种目标识别方法 | |
CN111046844B (zh) | 一种基于邻域选取约束的高光谱图像分类方法 | |
Xiong et al. | Small target detection for infrared image based on optimal infrared patch-image model by solving modified adaptive RPCA problem | |
CN110263777B (zh) | 基于空谱结合的局部保持投影算法的目标检测方法及系统 | |
CN105139353B (zh) | 一种置换混叠图像的盲分离方法 | |
Prasad et al. | Rotation invariance through structured sparsity for robust hyperspectral image classification | |
Kim et al. | FMPR-Net: False Matching Point Removal Network for Very-High-Resolution Satellite Image Registration | |
Sireesha et al. | Hyperspectral Image Segmentation and Classification using Hybrid IB-CDA with Least Square Unmixing | |
Cao et al. | Sparse representation denoising framework for 3-D building reconstruction from airborne LiDAR data |
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 |