一种基于可见光谱的空间碎片姿态分析方法
【专利摘要】一种基于可见光谱的空间碎片姿态分析方法，包括如下步骤：步骤1，空间碎片光谱数据采集，记为观测谱；步骤2，光谱数据预处理；步骤3，建立空间碎片的三维模型，确定各面元的法线方向及材料谱；步骤4，利用在轨测量，确定空间碎片质心到测站的方向矢量和空间碎片质心到太阳的方向矢量；步骤5，记空间碎片的偏航角、俯仰角、滚动角，计算旋转矩阵R；步骤6，计算在本体参考系下的空间碎片指向测站的方向矢量和空间碎片指向太阳的方向矢量；步骤7，由朗伯散射模型计算预测谱；步骤8，计算预测谱与观测谱之差的欧氏范数(残余)作为姿态指向判别度量；步骤9，估计目标指向，以最小残余所对应的偏航角、俯仰角、滚动角作为空间碎片姿态的空间指向。本发明有效适用于深空环境、准确性更高。
【专利说明】
-种基于可见光谱的空间碎片姿态分析方法
技术领域
[0001] 本发明设及载人航天实施空间监测技术领域，尤其是一种基于可见光谱的空间碎 片姿态分析方法。
【背景技术】
[0002] 自1957年10月前苏联发射首颗人造地球卫星W来，美国、法国、日本也相继发射了 人造地球卫星，中国也于1970年发射了第一颗自己的人造地球卫星"东方红一号"；据2013 年统计，截止2013年5月31号，在轨运行的人造地球卫星总数为1071颗;而目前具体有多少 颗卫星，并不确定，因为该数字几乎每天都会在变，有新的卫星发射，也有老卫星失效后被 废弃;并且随着小(纳、皮)卫星的发射，其数量变化更快。
[0003] 人类几十年的空间探索活动产生了大量的空间碎片，对人类航天活动的安全造成 严重威胁，成为空间环境的主要污染源，并在一定程度上对航天活动的正常开展产生了影 响。
[0004] 传统空间碎片测量W位置信息测量为主，包括空间碎片的Ξ维位置坐标、速度、加 速度等参数，可衍生出各类地球轨道参数。为了提高空间碎片的监测预警能力，对测量系统 除了要求获得其位置信息之外，更需要获得空间碎片的特征信息，如空间碎片姿态等特征 信息，它对空间监测、预警尤为重要。
[0005] 光谱测量技术是空间碎片状态分析的一种重要方法，空间碎片一般本身并不发 光，其亮度来自太阳光的反射，仅由亮度变化不足W分辨空间碎片的材料组成及空间指向， 分析其光谱特征成为我们辨认空间碎片姿态的主要手段。
[0006] 在地基(地面观测站)条件下，光谱仪所获取光谱由空间碎片材料成分、太阳光谱、 地球大气吸收谱，W及测量过程中产生的噪声等因素综合决定，因而由光谱推断空间碎片 指向并没有一个解析解。本发明采用穷举捜索方法，对空间碎片姿态指向作判别。
[0007] 空间碎片的测量有雷达测量和光学测量，雷达测量可W克服天气、太阳和及地影 的影响，能全天候全天时工作，但由于雷达测量时，其反射回波信号的强度与距离的四次方 成反比，因而雷达测量比较适合于低地球轨道的空间碎片；而对于无源光学测量，信号反射 强度与物体的距离的平方成反比，因而能探测高轨道空间碎片，另外它仅接收空间碎片对 太阳光的反射，属于被动探测形式，可较好的适合于某些特殊场合。
【发明内容】
[000引为了克服已有空间碎片姿态指向分析方法无法适用于深空环境、准确性较差的不 足，本发明提供一种有效适用于深空环境、准确性更高的基于可见光谱的空间碎片姿态分 析方法。
[0009] 本发明解决其技术问题所采用的技术方案是：
[0010] -种基于可见光谱的空间碎片姿态分析方法，所述方法包括如下步骤：
[0011] 步骤1，空间碎片光谱数据采集:利用光谱仪，采集空间碎片光谱数据，记为观测谱 Ytest;
[0012] 步骤2，光谱数据预处理:进行去噪处理；
[0013] 步骤3,建立空间碎片的3维模型，确定各面元的法线方向及材料谱；
[0014] 步骤4,利用在轨测量，确定空间碎片质屯、到测站的方向矢量0勺日空间碎片质屯、到 太阳的方向矢量S%
[0015] 步骤5,记空间碎片的偏航角Yaw、俯仰角Pitch、滚动角Roll,计算旋转矩阵R，有：
[0016] R = Rx(R〇ll)*Ry(Pitch)Rz(化 W)
[0017] 其中；
[0020] 步骤6，计算在本体参考系下的空间碎片指向测站的方向矢量0和空间碎片指向太 阳的方向矢量S
[0021]
[0022] 步骤7，由朗伯散射模型计算预测谱：
[0023]
[0024] 其中，A为k个材料谱所组成矩阵，ill表示第i面元的面积，m表示第i面元的法线矢 量，其中，点积运算符〈·，·〉还要求满足，对任意实数yER，y表示由点积运算结果值：
[0025]
[0026] 步骤8,计算预测谱Ypred与观测谱Ytest之差的欧氏范数(残余)作为姿态指向判别度 量，记残余：
[0027] 残余=I lYpred-YtestN2
[0028] 步骤9,估计空间碎片姿态指向，W最小残余所对应的偏航角Yaw、俯仰角Pitch、滚 动角Roll作为空间碎片姿态的空间指向。
[0029] 进一步，所述步骤2中，去噪处理的过程如下：
[0030] 2.1)阔值化；
[0031]
[0032] 选取阔值化，若光谱数据X大于该阔值，则设置该项数据为所选择的阔值化;若小 于ο，则设置该值为ο;其它情形，保持不变；
[0033] 2.2)将可见光范围数据规整到[400-700]醒，^消除二级谱的混入。
[0034] 本发明的有益效果主要表现在:有效适用于深空环境、准确性更高。
【附图说明】
[0035] 图1是一种基于可见光谱的空间碎片姿态分析方法的流程图。
【具体实施方式】
[0036] 下面结合附图对本发明作进一步描述。
[0037] 参照图1，一种基于可见光谱的空间碎片姿态分析方法，包括如下步骤：
[0038] 步骤1，空间碎片光谱数据采集:利用光谱仪，采集空间碎片光谱数据；
[0039] 步骤2,光谱数据预处理；
[0040] 由于测量仪器精度及大气端流等因素影响，所获取光谱数据中存在噪声，为此采 用如下方法作数据预处理：
[0041] 2.1)阔值化；
[0042]
[0043] 选取阔值化，若光谱数据X大于该阔值，则设置该项数据为所选择的阔值化;若小 于0，则设置该值为0;其它情形，保持不变。
[0044] 2.2)将可见光范围数据规整到[400-700]11111，^消除二级谱的混入；
[0045] 步骤3,建立空间碎片的3维模型，确定各面元的法线方向及材料谱；
[0046] 步骤4,利用在轨测量，确定空间碎片质屯、到测站的方向矢量0勺日空间碎片质屯、到 太阳的方向矢量S%
[0047] 步骤5,记空间碎片的偏航角(Yaw)、俯仰角(Pitch)、滚动角(Roll),计算旋转矩阵 R，有：
[004引 R = Rx(R0ll)*Ry(Pitch)Rz(化 W)
[0049]其中；
[0052]步骤6,计算在本体参考系下的空间碎片指向测站的方向矢量0和空间碎片指向太 阳的方向矢量S
[0化3]
[0054]步骤7，由朗伯散射模型计算预测谱：
[0化5]
[0056]其中A为k个材料谱所组成矩阵，日1表示第i面元的面积，m表示第i面元的法线矢 量，其中，点积运算符〈·，·〉还要求满足，对任意实数yER，y表示由点积运算结果值：
[0化7]
[0058]步骤8,计算预测谱Ypred与测试谱Ytest之差的欧氏范数(残余)作为姿态指向判别度 量，记残余：
[0059]残余=I |Ypred-Ytest| |2
[0060] 步骤9,估计空间碎片指向，W最小残余所对应的偏航角Yaw、俯仰角Pitch、滚动角 Roll作为空间碎片姿态的空间指向。
[0061] 实例:一种基于可见光谱的空间碎片姿态分析方法，包括如下步骤
[0062] (1)构造空间碎片3D模型;利用STK的V0模块设置验证空间碎片的3D模型；
[0063] (2)选择空间碎片表面覆盖材料谱，基于公开报道文献获取材料光谱数据；
[0064] (3)计算W空间碎片本体为参考系下的面元法线矢量，如单位立方体，则六面的法 线矢量表示为
[00 化]
[0066] (4)由仿真观测谱验证空间碎片姿态指向，取Ξ个姿态角的捜索范围为[-180， 180]，捜索步长为5°，设置的Ξ个姿态角为[13.6，-37.9，167.0];在前10个最小残余所对应 的姿态角中有两个满足我们的要求(估计姿态角与实际姿态角之差小于10°): [15.0,-40， 160]、[20，-35，165]。
【主权项】
1. 一种基于可见光谱的空间碎片姿态分析方法，其特征在于:所述方法包括如下步骤： 步骤1，空间碎片光谱数据采集:利用光谱仪，采集空间碎片光谱数据，记为观测谱Ytest; 步骤2，光谱数据预处理:进行去噪处理； 步骤3，建立空间碎片的3维模型，确定各面元的法线方向及材料谱； 步骤4,利用在轨测量，确定空间碎片质心到测站的方向矢量一和空间碎片质心到太阳 的方向矢量s% 步骤5,记空间碎片的偏航角Yaw、俯仰角Pitch、滚动角Roll，计算旋转矩阵R，有： R=Rx(Roll)*Ry(Pitch)Rz(Yaw) 其中：步骤6，计算在本体参考系下的空间碎片指向测站的方向矢量〇和空间碎片指向太阳的 方向矢量s步骤7，由朗伯散射模型计算预测谱：其中，A为k个材料谱所组成矩阵，ai表示第i面元的面积，m表示第i面元的法线矢量，其 中，点积运算符〈·，· >还要求满足，对任意实数y￡R，y表示由点积运算结果值：步骤8,计算预测谱Y-与观测谱Yt(3St之差的欧氏范数(残余)作为姿态指向判别度量， 记残余： 残余=| | Ypred-Ytest I I 2 步骤9,估计空间碎片指向，以最小残余所对应的偏航角Yaw、俯仰角Pitch、滚动角Roll 作为空间碎片姿态的空间指向。2. 如权利要求1所述的一种基于可见光谱的空间碎片姿态分析方法，其特征在于:所述 步骤2中，去噪处理的过程如下： 2.1)阈值化；选取阈值Th，若光谱数据X大于该阈值，则设置该项数据为所选择的阈值Th;若小于0， 则设置该值为〇;其它情形，保持不变； 2.2)将可见光范围数据规整到[400-700]nm，以消除二级谱的混入。
【文档编号】G06F19/00GK105825053SQ201610146136
【公开日】2016年8月3日
【申请日】2016年3月15日
【发明人】庄德文, 唐轶峻, 康湘华, 夏晨阳
【申请人】浙江工业大学