本发明专利技术公开了一种基于轨迹编目模型的空间点目标天基探测方法，本发明专利技术对经过预处理的图像进行阈值分割和连通域计算，并对每个连通域进行轨迹编目，根据每个轨迹的运动特征和判断概率来识别并跟踪目标。本发明专利技术适用于大视场可见光空间点目标天基探测，可以扩大搜索范围，实现尽早探测、尽快预警。可以同时探测多个不同运动方式的目标，并对多个目标进行识别和跟踪。本发明专利技术对于探测空间点目标可以实现快速探测识别和稳定跟踪。
【技术实现步骤摘要】
本专利技术涉及目标探测技术，具体涉及一种利用轨道编目模型对在复杂星空背景下 运动的空间点目标进行天基探测的方法，主要用于对空间目标的探测、识别、跟踪和测量。
技术介绍
天基空间目标探测系统是指对光学传感器获取的观测图像序列进行处理，从中提 取目标和背景恒星的位置、亮度等信息，并对空间目标进行定位跟踪、编目维护以及识别分 类。漆长松提出了一种星点间部最小距离特征概念，构建局部最小距离特征，根据点 群间的相对位置关系，采用递归搜索的方法对运动目标进行检测。见漆长松等，“天基平台 深空运动小目标检测”，光电工程，第35卷第5期。李骏提出对帧集内目标运动条痕线性拟合提取帧集测量信息的空间目标被动跟 踪方法。见李骏等，“基于天基光学监视的空间目标跟踪方法研究”，空间科学学报，2009， 29(3) :3洸-331。张健提出了基于序列图像多帧累加背景掩模帧生成的目标提取算法。见张健等， “空间目标运动轨迹提取算法研究”，红外技术，第四卷第8期。上述探测方法都没能很好的解决在无任何初始化条件下，在复杂的星空背景中， 对快速机动的多目标进行实时探测并长时间稳定跟踪的要求。
技术实现思路
本专利技术的目的在于提供一种在复杂星空背景中，对多个运动的空间点目标进行实 时探测并可长时间稳定跟踪的方法。该方法可以探测并跟踪多个不同运动方式的空间目 标，算法简单易用实现，可以减轻星载硬件设备的计算量。本专利技术提出利用轨迹编目的方法对不同的提取点进行编目，根据每个轨迹的运动 特征和判断概率来辨别哪些是目标。具体步骤包括(1)读取实时图像数据；(2)对图像进行图像预处理；(3)自适应确定图像阈值；(4)对图像进行阈值分割；(5)对图像进行连通域计算；(6)确定每个连通域的质心；(7)进行轨迹编目；(8)利用轨迹编目进行目标识别。本专利技术是利用轨道编目模型对复杂星空背景下的空间点目标进行天基探测。探 测波段选取在可见光波段，即380nm-780nm，探测系统为大视场可见光探测系统，视场为 20° X15°，探测器分辨率为10MX768，探测系统瞬时视场为1.2'。恒星最大张角为0.05〃，可以计算ImXlm的目标在1.5km处的张角为1.15'。根据以上分析可知，在可见 光范围内即能看到目标又能看到恒星，在15°视场内比5等星亮的星平均有6. 2颗，恒星在 视场内始终是点目标，ImX Im的目标在1. 5km以外也是点目标。从上面的分析可知，要探测 到目标，必须在恒星点和目标点中找到它们的区别。本专利技术采用对所有提取点全部进行轨 迹编目和轨迹维护，在维护的轨迹中进行目标探测和识别。识别的方法为(1)在所有轨迹中统计灰度变化，如果灰度变化超过阈值，则此轨迹编目被判为目 标的概率加1 ；(2)在所有轨迹中统计速度变化，首先统计出大众速度，如果轨迹的速度与大众速 度差值超过阈值，则此轨迹编目被判为目标的概率加1 ；(3)在所有轨迹中统计尺寸变化，如果尺寸变化超过阈值，则此轨迹编目被判为目 标的概率加1 ；最后，依据统计的概率来判断是否为目标，如果在整个轨迹编目中，被连续判为目 标的数目大于整个轨迹编目的数目的一半，即被判为目标。本专利技术的优点(1)本专利技术适用于大视场可见光空间点目标天基探测。探测系统采用大视场，可以 扩大搜索范围，实现尽早探测、尽快预警。但同时也使背景星点成倍增加，加大探测识别的 难度。本专利技术采用轨迹编目的方法，一旦进入稳定的轨迹维护，对视场大小不敏感，可以实 现大视场探测识别。(2)本专利技术适用于多目标同时探测、识别和跟踪。由于采用对视场内所有可疑点全 部进行轨迹编目，可以同时探测多个目标，并对多个目标进行识别和跟踪。(3)本专利技术对于探测空间点目标可以实现高稳定性和快速探测识别。由于对所有 轨迹编目的特征进行概率统计，使得判别目标的虚警率极大降低，增强判别目标的稳定性。 通常情况下，在10帧之内即可判别出目标，可以实现对空间目标的快速探测识别。附图说明图1为流程图。图2为连通域计算用的九宫图。图3为轨迹编目算法流程图。具体实施例方式下面结合附图和实例对本专利技术作进一步详细的说明。使用本专利技术对复杂星空背景中运动的空间点目标进行探测、识别和跟踪，其流程 图如图1所示，具体步骤为1)首先，对采集到的实时图像数据进行坏点剔除、畸变校正和图像滤波预处理；2)接着，对经过预处理的图像数据进行自适应阈值的确定，阈值D计算公式为D = aX (M+2)其中，M为整幅图像灰度均值，a常数，选取范围在1.8 2. 4之间；3)对采集到的图像进行阈值分割，阈值分割公式为权利要求1. 一种，其特征在于包括以下步骤1)首先，对采集到的实时图像数据进行坏点剔除、畸变校正和图像滤波预处理；2)接着，对经过预处理的图像数据进行自适应阈值的确定，阈值D计算公式为 D = aX (M+2)其中，M为整幅图像灰度均值，a常数，选取范围在1. 8 2. 4之间；3)对采集到的图像进行阈值分割，阈值分割公式为全文摘要本专利技术公开了一种，本专利技术对经过预处理的图像进行阈值分割和连通域计算，并对每个连通域进行轨迹编目，根据每个轨迹的运动特征和判断概率来识别并跟踪目标。本专利技术适用于大视场可见光空间点目标天基探测，可以扩大搜索范围，实现尽早探测、尽快预警。可以同时探测多个不同运动方式的目标，并对多个目标进行识别和跟踪。本专利技术对于探测空间点目标可以实现快速探测识别和稳定跟踪。文档编号G01V8/10GK102116876SQ20111000856公开日2011年7月6日 申请日期2011年1月14日 优先权日2011年1月14日专利技术者凌丽, 崔文楠, 张涛, 杨波, 王磊, 童广辉 申请人:中国科学院上海技术物理研究所本文档来自技高网...
【技术保护点】
１．一种基于轨迹编目模型的空间点目标天基探测方法，其特征在于包括以下步骤：１）首先，对采集到的实时图像数据进行坏点剔除、畸变校正和图像滤波预处理；２）接着，对经过预处理的图像数据进行自适应阈值的确定，阈值Ｄ计算公式为：Ｄ＝ａ×（Ｍ＋２）其中，Ｍ为整幅图像灰度均值，ａ常数，选取范围在１．８～２．４之间；３）对采集到的图像进行阈值分割，阈值分割公式为：（ｍａｔｈ）??（ｍｆｅｎｃｅｄｏｐｅｎ＝＇｛＇ｃｌｏｓｅ＝＇＇）?（ｍｔａｂｌｅ）?（ｍｔｒ）?（ｍｔｄ）?（ｍｉ）ｉｆ（／ｍｉ）?（／ｍｔｄ）?（ｍｔｄ）?（ｍｉ）ｆ（／ｍｉ）?（ｍｒｏｗ）?（ｍｏ）（（／ｍｏ）?（ｍｉ）ｘ（／ｍｉ）?（ｍｏ），（／ｍｏ）?（ｍｉ）ｙ（／ｍｉ）?（ｍｏ））（／ｍｏ）?（／ｍｒｏｗ）?（ｍｏ）＆ＧｒｅａｔｅｒＥｑｕａｌ；（／ｍｏ）?（ｍｉ）Ｄ（／ｍｉ）?（ｍｏ），（／ｍｏ）?（ｍｉ）ｆ（／ｍｉ）?（ｍｒｏｗ）?（ｍｏ）（（／ｍｏ）?（ｍｉ）ｘ（／ｍｉ）?（ｍｏ），（／ｍｏ）?（ｍｉ）ｙ（／ｍｉ）?（ｍｏ））（／ｍｏ）?（／ｍｒｏｗ）?（ｍｏ）＝（／ｍｏ）?（ｍｉ）ｆ（／ｍｉ）?（ｍｒｏｗ）?（ｍｏ）（（／ｍｏ）?（ｍｉ）ｘ（／ｍｉ）?（ｍｏ），（／ｍｏ）?（ｍｉ）ｙ（／ｍｉ）?（ｍｏ））（／ｍｏ）?（／ｍｒｏｗ）?（／ｍｔｄ）?（／ｍｔｒ）?（ｍｔｒ）?（ｍｔｄ）?（ｍｉ）ｉｆ（／ｍｉ）?（／ｍｔｄ）?（ｍｔｄ）?（ｍｉ）ｆ（／ｍｉ）?（ｍｒｏｗ）?（ｍｏ）（（／ｍｏ）?（ｍｉ）ｘ（／ｍｉ）?（ｍｏ），（／ｍｏ）?（ｍｉ）ｙ（／ｍｉ）?（ｍｏ））（／ｍｏ）?（／ｍｒｏｗ）?（ｍｏ）＆ｌｔ；（／ｍｏ）?（ｍｉ）Ｄ（／ｍｉ）?（ｍｏ），（／ｍｏ）?（ｍｉ）ｆ（／ｍｉ）?（ｍｒｏｗ）?（ｍｏ）（（／ｍｏ）?（ｍｉ）ｘ（／ｍｉ）?（ｍｏ），（／ｍｏ）?（ｍｉ）ｙ（／ｍｉ）?（ｍｏ））（／ｍｏ）?（／ｍｒｏｗ）?（ｍｏ）＝（／ｍｏ）?（ｍｎ）０（／ｍｎ）?（／ｍｔｄ）?（／ｍｔｒ）?（／ｍｔａｂｌｅ）?（／ｍｆｅｎｃｅｄ）?（／ｍａｔｈ）其中，ｆ（ｘ，ｙ）为图像坐标（ｘ，ｙ）处的灰度值；４）对经过阈值分割的图像进行连通域计算，计算连通域的方法如下：（１）从像元（１，１）开始，按照自上而下、从左到右的顺序逐行处理每一个像元，如果ｆ（ｘ，ｙ）≥Ｄ，则执行步骤（２），否则执行步骤（３）；（２）在已经处理的相邻４个像元（ｘ－１，ｙ－１）、（ｘ－１，ｙ）、（ｘ－１，ｙ＋１）和（ｘ，ｙ－１）中，统计不同连通域编号的个数，如果不同连通域编号的个数为０，待处理的像元（ｘ，ｙ）属于新的连通域，在已经编号的连通域中取最大值并且加１赋予像元（ｘ，ｙ）新的连通域编号；如果不同连通域编号的个数为１，则待处理的像元（ｘ，ｙ）和一个已知连通域连通，将此连通域的编号赋予像元（ｘ，ｙ）；如果不同连通域编号的个数为２，则待处理的像元（ｘ，ｙ）和两个已知连通域连通，将两个连通域中较小的编号赋予像元（ｘ，ｙ），并用较小的编号替代较大编号的连通域；（３）如果全部像元处理完毕，执行步骤（４），否则，像元指针加１，执行步骤（２）；（４）在不更改原有排列顺序的前提下，按照编号不间断的原则更新连通域的编号；５）计算每个连通域的质心，质心坐标ｘｃ，ｙｃ计算公式为：（ｍａｔｈ）??（ｍｒｏｗ）?（ｍｓｕｂ）?（ｍｉ）ｘ（／ｍｉ）?（ｍｉ）ｃ（／ｍｉ）?（／ｍｓｕｂ）?（ｍｏ）＝（／ｍｏ）?（ｍｆｒａｃ）?（ｍｒｏｗ）?（ｍｉ）＆Ｓｉｇｍａ；ｘ（／ｍｉ）?（ｍｏ）＆ＣｅｎｔｅｒＤｏｔ；（／ｍｏ）?（ｍｉ）ｆ（／ｍｉ）?（ｍｒｏｗ）?（ｍｏ）（（／ｍｏ）?（ｍｉ）ｘ（／ｍｉ）?（ｍｏ），（／ｍｏ）?（ｍｉ）ｙ（／ｍｉ）?（ｍｏ））（／ｍｏ）?（／ｍｒｏｗ）?（／ｍｒｏｗ）?（ｍｒｏｗ）?（ｍｉ）＆Ｓｉｇｍａ；ｆ（／ｍｉ）?（ｍｒｏｗ）?（ｍｏ）（（／ｍｏ）?（ｍｉ）ｘ（／ｍｉ）?（ｍｏ），（／ｍｏ）?（ｍｉ）ｙ（／ｍｉ）?（ｍｏ））（／ｍｏ）?（／ｍｒｏｗ）?（／ｍｒｏｗ）?（／ｍｆｒａｃ）?（／ｍｒｏｗ）?（／ｍａｔｈ）（ｍａｔｈ）??（ｍｒｏｗ）?（ｍｓｕｂ）?（ｍｉ）ｙ（／ｍｉ）?（ｍｉ）ｃ（／ｍｉ）?（／ｍｓｕｂ）?（ｍｏ）＝（／ｍｏ）?（ｍｆｒａｃ）?（ｍｒｏｗ）?（ｍｉ）＆Ｓｉｇｍａ；ｙ（／ｍｉ）?（ｍｏ）＆ＣｅｎｔｅｒＤｏｔ；（／ｍｏ）?（ｍｉ）ｆ（／ｍｉ）?（ｍｒｏｗ）?（ｍｏ）（（／ｍｏ）?（ｍｉ）ｘ（／ｍｉ）?（ｍｏ），（／ｍｏ）?（ｍｉ）ｙ（／ｍｉ）?（ｍｏ））（／ｍｏ）?（／ｍｒｏｗ）?（／ｍｒｏｗ）?（ｍｒｏｗ）?（ｍｉ）＆Ｓｉｇｍａ；ｆ（／ｍｉ）?（ｍｒｏｗ）?（ｍｏ）（（／ｍｏ）?（ｍｉ）ｘ（／ｍｉ）?（ｍｏ），（／ｍｏ）?（ｍｉ）ｙ（／ｍｉ）?（ｍｏ））（／ｍｏ）?（／ｍｒｏｗ）?（...
【技术特征摘要】
【专利技术属性】
技术研发人员：崔文楠，张涛，王磊，凌丽，童广辉，杨波，
申请(专利权)人：中国科学院上海技术物理研究所，
类型：发明
国别省市：31
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。