发明名称 一种基于可见光谱的空间碎片姿态分析方法
摘要 一种基于可见光谱的空间碎片姿态分析方法,包括如下步骤:步骤1,空间碎片光谱数据采集,记为观测谱;步骤2,光谱数据预处理;步骤3,建立空间碎片的三维模型,确定各面元的法线方向及材料谱;步骤4,利用在轨测量,确定空间碎片质心到测站的方向矢量和空间碎片质心到太阳的方向矢量;步骤5,记空间碎片的偏航角、俯仰角、滚动角,计算旋转矩阵R;步骤6,计算在本体参考系下的空间碎片指向测站的方向矢量和空间碎片指向太阳的方向矢量;步骤7,由朗伯散射模型计算预测谱;步骤8,计算预测谱与观测谱之差的欧氏范数(残余)作为姿态指向判别度量;步骤9,估计目标指向,以最小残余所对应的偏航角、俯仰角、滚动角作为空间碎片姿态的空间指向。本发明有效适用于深空环境、准确性更高。
申请公布号 CN105825053A 申请公布日期 2016.08.03
申请号 CN201610146136.2 申请日期 2016.03.15
申请人 浙江工业大学 发明人 庄德文;唐轶峻;康湘华;夏晨阳
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 杭州斯可睿专利事务所有限公司 33241 代理人 王利强
主权项 一种基于可见光谱的空间碎片姿态分析方法,其特征在于:所述方法包括如下步骤:步骤1,空间碎片光谱数据采集:利用光谱仪,采集空间碎片光谱数据,记为观测谱Y<sub>test</sub>;步骤2,光谱数据预处理:进行去噪处理;步骤3,建立空间碎片的3维模型,确定各面元的法线方向及材料谱;步骤4,利用在轨测量,确定空间碎片质心到测站的方向矢量o<sup>*</sup>和空间碎片质心到太阳的方向矢量s<sup>*</sup>;步骤5,记空间碎片的偏航角Yaw、俯仰角Pitch、滚动角Roll,计算旋转矩阵R,有:R=R<sub>x</sub>(Roll)*R<sub>y</sub>(Pitch)R<sub>z</sub>(Yaw)其中:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>R</mi><mi>x</mi></msub><mrow><mo>(</mo><mrow><mi>R</mi><mi>o</mi><mi>l</mi><mi>l</mi></mrow><mo>)</mo></mrow><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mi>cos</mi><mrow><mo>(</mo><mrow><mi>R</mi><mi>o</mi><mi>l</mi><mi>l</mi></mrow><mo>)</mo></mrow></mrow></mtd><mtd><mrow><mo>-</mo><mi>sin</mi><mrow><mo>(</mo><mrow><mi>R</mi><mi>o</mi><mi>l</mi><mi>l</mi></mrow><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mi>sin</mi><mrow><mo>(</mo><mrow><mi>R</mi><mi>o</mi><mi>l</mi><mi>l</mi></mrow><mo>)</mo></mrow></mrow></mtd><mtd><mrow><mi>cos</mi><mrow><mo>(</mo><mrow><mi>R</mi><mi>o</mi><mi>l</mi><mi>l</mi></mrow><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced><mo>,</mo><msub><mi>R</mi><mi>y</mi></msub><mrow><mo>(</mo><mrow><mi>P</mi><mi>i</mi><mi>t</mi><mi>c</mi><mi>h</mi></mrow><mo>)</mo></mrow><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><mi>cos</mi><mrow><mo>(</mo><mrow><mi>P</mi><mi>i</mi><mi>t</mi><mi>c</mi><mi>h</mi></mrow><mo>)</mo></mrow></mrow></mtd><mtd><mn>0</mn></mtd><mtd><mrow><mi>sin</mi><mrow><mo>(</mo><mrow><mi>P</mi><mi>i</mi><mi>t</mi><mi>c</mi><mi>h</mi></mrow><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mrow><mo>-</mo><mi>sin</mi><mrow><mo>(</mo><mrow><mi>P</mi><mi>i</mi><mi>t</mi><mi>c</mi><mi>h</mi></mrow><mo>)</mo></mrow></mrow></mtd><mtd><mn>0</mn></mtd><mtd><mrow><mi>cos</mi><mrow><mo>(</mo><mrow><mi>P</mi><mi>i</mi><mi>t</mi><mi>c</mi><mi>h</mi></mrow><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000941430160000011.GIF" wi="1526" he="220" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>R</mi><mi>z</mi></msub><mrow><mo>(</mo><mi>Y</mi><mi>a</mi><mi>w</mi><mo>)</mo></mrow><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><mi>c</mi><mi>o</mi><mi>s</mi><mrow><mo>(</mo><mi>Y</mi><mi>a</mi><mi>w</mi><mo>)</mo></mrow></mrow></mtd><mtd><mrow><mo>-</mo><mi>s</mi><mi>i</mi><mi>n</mi><mrow><mo>(</mo><mi>Y</mi><mi>a</mi><mi>w</mi><mo>)</mo></mrow></mrow></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mrow><mi>s</mi><mi>i</mi><mi>n</mi><mrow><mo>(</mo><mi>Y</mi><mi>a</mi><mi>w</mi><mo>)</mo></mrow></mrow></mtd><mtd><mrow><mi>cos</mi><mrow><mo>(</mo><mi>Y</mi><mi>a</mi><mi>w</mi><mo>)</mo></mrow></mrow></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000941430160000012.GIF" wi="694" he="198" /></maths>步骤6,计算在本体参考系下的空间碎片指向测站的方向矢量o和空间碎片指向太阳的方向矢量s<maths num="0003" id="cmaths0003"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mi>s</mi><mo>=</mo><mi>R</mi><mo>*</mo><msup><mi>s</mi><mo>*</mo></msup></mrow></mtd></mtr><mtr><mtd><mrow><mi>o</mi><mo>=</mo><mi>R</mi><mo>*</mo><msup><mi>o</mi><mo>*</mo></msup></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000941430160000013.GIF" wi="221" he="158" /></maths>步骤7,由朗伯散射模型计算预测谱:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>Y</mi><mrow><mi>p</mi><mi>r</mi><mi>e</mi><mi>d</mi></mrow></msub><mo>=</mo><mi>A</mi><mo>*</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>a</mi><mn>1</mn></msub><mo>&lt;</mo><msub><mi>n</mi><mn>1</mn></msub><mo>,</mo><mi>o</mi><mo>&gt;</mo><mo>&lt;</mo><msub><mi>n</mi><mn>1</mn></msub><mo>,</mo><mi>s</mi><mo>&gt;</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>a</mi><mi>k</mi></msub><mo>&lt;</mo><msub><mi>n</mi><mi>k</mi></msub><mo>,</mo><mi>o</mi><mo>&gt;</mo><mo>&lt;</mo><msub><mi>n</mi><mi>k</mi></msub><mo>,</mo><mi>s</mi><mo>&gt;</mo></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000941430160000014.GIF" wi="651" he="223" /></maths>其中,A为k个材料谱所组成矩阵,a<sub>i</sub>表示第i面元的面积,n<sub>i</sub>表示第i面元的法线矢量,其中,点积运算符&lt;·,·&gt;还要求满足,对任意实数y∈R,y表示由点积运算结果值:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mo>&lt;</mo><mi>y</mi><mo>&gt;</mo><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mi>y</mi></mtd><mtd><mrow><mi>y</mi><mo>&gt;</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mi>y</mi><mo>&le;</mo><mn>0</mn></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000941430160000021.GIF" wi="356" he="143" /></maths>步骤8,计算预测谱Y<sub>pred</sub>与观测谱Y<sub>test</sub>之差的欧氏范数(残余)作为姿态指向判别度量,记残余:残余=||Y<sub>pred</sub>‑Y<sub>test</sub>||<sub>2</sub>步骤9,估计空间碎片指向,以最小残余所对应的偏航角Yaw、俯仰角Pitch、滚动角Roll作为空间碎片姿态的空间指向。
地址 310014 浙江省杭州市下城区朝晖六区潮王路18号浙江工业大学