CN109239384B - 一种二维三维融合的非合作目标转速转轴测量方法 - Google Patents
一种二维三维融合的非合作目标转速转轴测量方法 Download PDFInfo
- Publication number
- CN109239384B CN109239384B CN201811232918.3A CN201811232918A CN109239384B CN 109239384 B CN109239384 B CN 109239384B CN 201811232918 A CN201811232918 A CN 201811232918A CN 109239384 B CN109239384 B CN 109239384B
- Authority
- CN
- China
- Prior art keywords
- frequency
- dimensional
- value
- rotating speed
- target
- 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
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000003384 imaging method Methods 0.000 claims abstract description 9
- 230000000737 periodic effect Effects 0.000 claims description 22
- 238000004364 calculation method Methods 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 9
- 238000009616 inductively coupled plasma Methods 0.000 claims description 6
- 238000004458 analytical method Methods 0.000 claims description 3
- 238000010183 spectrum analysis Methods 0.000 claims description 3
- 238000012360 testing method Methods 0.000 claims description 2
- 238000012545 processing Methods 0.000 abstract description 17
- 238000001914 filtration Methods 0.000 abstract description 2
- 238000005259 measurement Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 230000004927 fusion Effects 0.000 description 3
- 238000005096 rolling process Methods 0.000 description 3
- 238000007499 fusion processing Methods 0.000 description 2
- 238000004140 cleaning Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000011045 prefiltration Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P3/00—Measuring linear or angular speed; Measuring differences of linear or angular speeds
- G01P3/36—Devices characterised by the use of optical means, e.g. using infrared, visible, or ultraviolet light
- G01P3/38—Devices characterised by the use of optical means, e.g. using infrared, visible, or ultraviolet light using photographic means
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Electromagnetism (AREA)
- Power Engineering (AREA)
- General Physics & Mathematics (AREA)
- Image Analysis (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
一种二维三维融合的非合作目标转速转轴测量方法,通过对序列二维图像进行处理得到目标的转速信息,再对序列的三维点云进行处理,同时得到目标的转速和转轴信息,并以基于图像信息得到的目标转速为依据,对基于三维点云信息得到的目标转轴转速信息进行滤波,得到目标的旋转轴信息,达到综合利用二维和三维成像敏感器观察旋转目标获取的二维图像和三维点云数据,并进行目标的转速和转轴的测量与计算的目的。
Description
技术领域
本发明涉及一种二维三维融合的非合作目标转速转轴测量方法,属于空间旋转非合作目标测量、图像处理和三维点云数据处理领域。
背景技术
面向空间非合作目标的在轨服务技术可实现对目标航天器的在轨维护、空间碎片清理、抓捕或破坏敌方军事卫星等任务,具有很高的社会效益和军事价值。为执行各种在轨任务中,服务航天器在捕获、跟踪、抵近和抓捕目标时,都需对非合作目标的运动参数进行测量。因此,面向空间非合作目标的测量技术和方法是空间技术领域的热点研究方向。
“非合作”体现在三个方面:无信息交换、无预先设计标识和姿态不可控,因此为数众多的空间失效卫星处于翻滚旋转状态。在空间目标失控翻滚/旋转状态下,对其实施抓捕、维修、加注燃料等操作将会十分困难。因此,需要再抵近翻滚旋转目标的过程中对其进行转速转轴等运动参数的测量。
面向空间旋转非合作目标的转速转轴测量方法包括基于三维成像敏感器获取的序列三维点云的匹配结算、基于二维成像相机序列图像的特征匹配算法,以及基于双目视觉技术的特征点追踪方法。以上几种方式均能从观测数据中结算出目标的转轴和转速信息。但是,对基于序列点云处理获取目标转速转轴信息的方法受点云噪声以及匹配精度的影响,其结果精度一般,波动较大;而基于序列图像的方法需要先提取特征点,然后再利用图像特征匹配的方式进行相邻帧之间的相对姿态和目标转轴的求解,运算量大,效率较低;而基于双目视觉技术的特征点追踪方法,需要进行左、右图同名点特征点匹配,效率不高。本专利将二维图像数据处理结果和三维点云数据处理结果进行融合处理,利用灰度区间划分方式计算目标转速,再以此结果作为依据,对三维点云的处理结果进行滤波,提高了序列点云转轴、转速结算结果的精度,获取了旋转目标的高精度转速、转轴信息,实现了图像处理结果和点云处理结果的融合输出。
发明内容
本发明解决的技术问题是:针对目前现有技术中,单独的图像处理匹配算法容易受到外界干扰、图像波动较大,计算精度一般的问题,提出了一种二维三维融合的非合作目标转速转轴测量方法。
本发明解决上述技术问题是通过如下技术方案予以实现的:
一种二维三维融合的非合作目标转速转轴测量方法,步骤如下:
(1)利用二维成像敏感器对旋转的非合作目标进行观测,得到二维序列图像,对任一图像进行直方图分析,并获取该图像的灰度值分布;
(2)根据步骤(1)所得图像的灰度值分布进行灰度区间划分,得到n个灰度区间;
(3)实时统计步骤(2)所得n个灰度区间内的像素个数并进行记录,统计特定时长的像素个数后,获取n个像素个数关于统计时长的周期性信号;
(4)对步骤(3)所得n个周期性信号进行频谱分析,于每个周期性信号中取幅值最大的频率作为该周期性信号的频率,记为f1,f2…fn;
(5)根据最大幅值对应频率和最大幅值一半所对应的频率之间的差值,设定相对频率阈值,并对步骤(4)所得周期性信号的频率序列f1,f2…fn进行检索分类,对分类后数量较多的一组频率数据取频率平均值作为输出频率,再根据输出频率计算目标转速;
(6)利用三维成像敏感器对旋转的非合作目标进行观测并获取三维点云数据,利用ICP算法对相邻帧点云进行匹配并获取相邻帧点云的旋转矩阵,计算该非合作目标的转轴、旋转角度、转速。
所述步骤(2)中,划分灰度区间的方法为:
选取n个图像灰度的局部极大值H1,H2…Hn,取每个局部极大值波动±M范围内的区间作为一个灰度区间,得到分别为 [H1-M,H1+M],[H2-M,H2+M]…[Hn-M,Hn+M]的n个灰度区间,其中,n为正整数,M 为图像灰度饱和值的1/20。
所述步骤(5)中,选取频率值的具体方法为:
(5a)于频率序列中任选一个周期性信号的频率fb作为A类频率值,计算其他频率值相对于该频率值的相对频率差dF,计算公式如下:
式中,fi为除去fb后频率序列中任意一个周期性信号频率;
(5b)根据最大幅值对应频率和最大幅值一半所对应的频率之间的差值设定相对频率阈值Δf,同时根据步骤(5a)所得相对频率差dF与Δf进行比较,若dF小于Δf,则fi与fb均为A类频率值,否则fi为B类频率值,对频率序列内全部频率值进行判断后得到A类频率值序列及B类频率值序列;
(5c)根据步骤(5b)所得的A类、B类频率值序列内频率值个数进行比较,取数值较多的一组频率值序列求频率平均值。
所述步骤(6)中,计算非合作目标的转轴、旋转角度、转速的公式为:
式中,[q0,q1,q2,q3]为转轴四元数,R为目标旋转矩阵,mij为目标旋转矩阵R中,第i行第j列的元素,θ为相邻帧之间的夹角,t为相邻帧之间的时间间隔,ωpointcloud为角速度,[nx,ny,nz]为转轴向量。
本发明与现有技术相比的优点在于:
(1)本发明提供的一种二维三维融合的非合作目标转速转轴测量方法,基于图像数据灰度区间内像素数统计的方法能够准确第求解出旋转目标的转速,将二维图像数据处理结果和三维点云数据处理结果进行融合处理,利用灰度区间划分方式计算目标转速,再以此结果作为依据,对三维点云的处理结果进行滤波,提高了序列点云转轴、转速结算结果的精度,获取了旋转目标的高精度转速、转轴信息,实现了图像处理结果和点云处理结果的融合输出;
(2)本发明依据基于序列图像测量得到的高精度转速信息为以及,对基于序列点云数据测量得到的转速转轴信息进行滤波,可有效剔除波动较大的转轴测量数据,大大提高转轴的测量精度,同时提高了计算效率。
附图说明
图1为发明提供的计算方法流程图;
图2为发明提供的图像灰度值分布示意图;
图3为发明提供的周期性信号示意图;
图4为发明提供的滤波前转轴计算结果图;
图5为发明提供的滤波后转轴计算结果图
具体实施方式
一种二维三维融合的非合作目标转速转轴测量方法,如图1所示,具体步骤如下:
(1)利用二维成像敏感器对旋转的非合作目标进行观测,得到一系列的二维序列图像,对任一图像进行直方图分析,并获取该图像的灰度值分布;
(2)根据步骤(1)所得图像的灰度值分布进行灰度区间划分,得到n个灰度区间,其中,选取灰度空间的具体方法为:
如图2所示,选取n个图像灰度的局部极大值H1,H2…Hn,取每个局部极大值波动±M范围内的区间作为一个灰度区间,得到分别为 [H1-M,H1+M],[H2-M,H2+M]…[Hn-M,Hn+M]的n个灰度区间,其中,n为正整数,M 为图像灰度饱和值的1/20;
(3)实时统计步骤(2)所得n个灰度区间内的像素个数并进行记录,统计特定时长的像素个数后,获取n个像素个数关于统计时长的周期性信号,如图3所示;
(4)对步骤(3)所得n个周期性信号进行频谱分析,于每个周期性信号中取幅值最大的频率作为该周期性信号的频率,记为f1,f2…fn;
(5)根据最大幅值对应频率和最大幅值一半所对应的频率之间的差值,设定相对频率阈值,并对步骤(4)所得周期性信号的频率序列f1,f2…fn进行检索分类,对分类后数量较多的一组频率数据取频率平均值作为输出频率,再根据输出频率计算目标转速,其中,选取频率值的具体方法为:
(5a)于频率序列中任选一个周期性信号的频率fb作为A类频率值,计算其他频率值相对于该频率值的相对频率差dF,计算公式如下:
式中,fi为除去fb后频率序列中任意一个周期性信号频率;
(5b)根据最大幅值对应频率和最大幅值一半所对应的频率之间的差值设定相对频率阈值Δf,同时根据步骤(5a)所得相对频率差dF与Δf进行比较,若dF小于Δf,则fi与fb均为A类频率值,否则fi为B类频率值,对频率序列内全部频率值进行判断后得到A类频率值序列及B类频率值序列;
(5c)根据步骤(5b)所得的A类、B类频率值序列内频率值个数进行比较,取数值较多的一组频率值序列求频率平均值;
(6)利用三维成像敏感器对旋转的非合作目标进行观测并获取三维点云数据,利用ICP算法对相邻帧点云进行匹配并获取相邻帧点云的旋转矩阵,计算该非合作目标的转轴、旋转角度、转速;
计算非合作目标的转轴、旋转角度、转速的公式为:
式中,[q0,q1,q2,q3]为转轴四元数,R为目标旋转矩阵,mij为目标旋转矩阵R中,第i行第j列的元素,θ为相邻帧之间的夹角,t为相邻帧之间的时间间隔,ωpointcloud为角速度,[nx,ny,nz]为转轴向量。
(7)在试验当中,还需要根据步骤(5)所得平均频率作为转速阈值,并对步骤(6)所得的转速数值进行筛选,选取转速数值大于转速阈值的部分并获取其他转轴信息。
下面结合具体实施例进行进一步说明:
如图4、图5所示,为依照本发明提供的转速转轴测量方法得出的测量结果。
本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。
Claims (4)
1.一种二维三维融合的非合作目标转速转轴测量方法,其特征在于步骤如下:
(1)利用二维成像敏感器对旋转的非合作目标进行观测,得到二维序列图像,对任一图像进行直方图分析,并获取该图像的灰度值分布;
(2)根据步骤(1)所得图像的灰度值分布进行灰度区间划分,得到n个灰度区间;
(3)实时统计步骤(2)所得n个灰度区间内的像素个数并进行记录,统计特定时长的像素个数后,获取n个像素个数关于统计时长的周期性信号;
(4)对步骤(3)所得n个周期性信号进行频谱分析,于每个周期性信号中取幅值最大的频率作为该周期性信号的频率,记为f1,f2…fn;
(5)根据最大幅值对应频率和最大幅值一半所对应的频率之间的差值,设定相对频率阈值,并对步骤(4)所得周期性信号的频率序列f1,f2…fn进行检索分类,对分类后数量较多的一组频率数据取频率平均值作为输出频率,再根据输出频率计算目标转速;
(6)利用三维成像敏感器对旋转的非合作目标进行观测并获取三维点云数据,利用ICP算法对相邻帧点云进行匹配并获取相邻帧点云的旋转矩阵,计算该非合作目标的转轴、旋转角度、转速;
(7)在试验当中,根据步骤(5)所得平均频率作为转速阈值,并对步骤(6)所得的转速数值进行筛选,选取转速数值大于转速阈值的部分并获取其他转轴信息。
2.根据权利要求1所述的一种二维三维融合的非合作目标转速转轴测量方法,其特征在于:所述步骤(2)中,划分灰度区间的方法为:
选取n个图像灰度的局部极大值H1,H2…Hn,取每个局部极大值波动±M范围内的区间作为一个灰度区间,得到分别为[H1-M,H1+M],[H2-M,H2+M]…[Hn-M,Hn+M]的n个灰度区间,其中,n为正整数,M为图像灰度饱和值的1/20。
3.根据权利要求1或2所述的一种二维三维融合的非合作目标转速转轴测量方法,其特征在于:所述步骤(5)中,选取频率值的具体方法为:
(5a)于频率序列中任选一个周期性信号的频率fb作为A类频率值,计算其他频率值相对于该频率值的相对频率差dF,计算公式如下:
式中,fi为除去fb后频率序列中任意一个周期性信号频率;
(5b)根据最大幅值对应频率和最大幅值一半所对应的频率之间的差值设定相对频率阈值Δf,同时根据步骤(5a)所得相对频率差dF与Δf进行比较,若dF小于Δf,则fi与fb均为A类频率值,否则fi为B类频率值,对频率序列内全部频率值进行判断后得到A类频率值序列及B类频率值序列;
(5c)根据步骤(5b)所得的A类、B类频率值序列内频率值个数进行比较,取数值较多的一组频率值序列求频率平均值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811232918.3A CN109239384B (zh) | 2018-10-23 | 2018-10-23 | 一种二维三维融合的非合作目标转速转轴测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811232918.3A CN109239384B (zh) | 2018-10-23 | 2018-10-23 | 一种二维三维融合的非合作目标转速转轴测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109239384A CN109239384A (zh) | 2019-01-18 |
CN109239384B true CN109239384B (zh) | 2021-02-09 |
Family
ID=65081119
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811232918.3A Active CN109239384B (zh) | 2018-10-23 | 2018-10-23 | 一种二维三维融合的非合作目标转速转轴测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109239384B (zh) |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AT399222B (de) * | 1992-10-19 | 1995-04-25 | Tabarelli Werner | Interferometrische einrichtung zur messung der lage eines reflektierenden objektes |
FR2919731A1 (fr) * | 2007-08-03 | 2009-02-06 | Thales Sa | Architecture radar modulaire |
CN101493317B (zh) * | 2009-03-06 | 2010-11-10 | 北京理工大学 | 抗径向跳动的非合作目标激光转角及速率测量方法与装置 |
CN102981011B (zh) * | 2012-12-10 | 2014-10-29 | 西北工业大学 | 一种基于双目视觉的空间目标旋转状态辨识方法 |
CN107644433B (zh) * | 2017-08-30 | 2020-07-14 | 北京控制工程研究所 | 改进的最近点迭代点云配准方法 |
CN108519075B (zh) * | 2018-03-22 | 2021-04-13 | 北京控制工程研究所 | 一种空间多目标位姿测量方法 |
-
2018
- 2018-10-23 CN CN201811232918.3A patent/CN109239384B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN109239384A (zh) | 2019-01-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US5801970A (en) | Model-based feature tracking system | |
CN108734103B (zh) | 卫星视频中运动目标的检测与跟踪方法 | |
CN102496015B (zh) | 二维高斯分布光斑图像中心快速高精度定位方法 | |
CN105976353A (zh) | 基于模型和点云全局匹配的空间非合作目标位姿估计方法 | |
CN108960190B (zh) | 基于fcn图像序列模型的sar视频目标检测方法 | |
CN111292369B (zh) | 激光雷达的伪点云数据生成方法 | |
CN113627473A (zh) | 基于多模态传感器的水面无人艇环境信息融合感知方法 | |
Qian et al. | Robust visual-lidar simultaneous localization and mapping system for UAV | |
CN109724586A (zh) | 一种融合深度图和点云的航天器相对位姿测量方法 | |
CN112990128A (zh) | 一种基于视频跟踪的多车辆测速方法 | |
CN113744315A (zh) | 一种基于双目视觉的半直接视觉里程计 | |
CN115588144A (zh) | 基于高斯动态阈值筛选的实时姿态捕捉方法及装置、设备 | |
CN114689030A (zh) | 一种基于机载视觉的无人机辅助定位方法及系统 | |
CN112580683B (zh) | 一种基于互相关的多传感器数据时间对齐系统及其方法 | |
CN109239384B (zh) | 一种二维三维融合的非合作目标转速转轴测量方法 | |
CN104156977A (zh) | 基于多线列时差扫描扩展采样的点目标运动速度探测方法 | |
CN117173215A (zh) | 一种跨摄像头的内河航船全程轨迹识别方法及系统 | |
CN102123234B (zh) | 无人机侦察视频分级运动补偿方法 | |
CN111031258A (zh) | 月球车导航相机曝光参数确定方法及装置 | |
CN114170376B (zh) | 一种针对室外大场景的多源信息融合分组式运动恢复结构方法 | |
CN116091804A (zh) | 一种基于相邻帧构型匹配的恒星抑制方法 | |
JP5465001B2 (ja) | 対象推定装置 | |
CN111899289B (zh) | 一种基于图像特征信息的红外图像与可见光图像的配准方法 | |
US5116118A (en) | Geometric fiedlity of imaging systems employing sensor arrays | |
CN114111602A (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 |