发明名称 一种基于成像传感器视差的卫星姿态颤振探测方法
摘要 本发明涉及一种基于成像传感器视差的卫星姿态颤振探测方法,包括以下步骤:1)获取视差影像,并对其进行预处理,所述的视差影像包括相互间具有观测视差的参考影像和待匹配影像;2)对视差影像进行精确密集匹配处理,获得两幅影像的视差偏差值;3)根据视差偏差值计算卫星姿态颤振值。与现有技术相比,本发明可在不依赖地面控制点和附加硬件设备的基础上,对星体的姿态颤振进行探测和分析,得到真实的星体在垂轨和沿轨方向的姿态颤振值。
申请公布号 CN103791899B 申请公布日期 2016.02.10
申请号 CN201410051078.6 申请日期 2014.02.14
申请人 同济大学 发明人 童小华;徐聿升;叶真;刘世杰;李凌云;李天鹏;王凤香
分类号 G01C21/00(2006.01)I 主分类号 G01C21/00(2006.01)I
代理机构 上海科盛知识产权代理有限公司 31225 代理人 赵继明
主权项 一种基于成像传感器视差的卫星姿态颤振探测方法,其特征在于,包括以下步骤:1)获取视差影像,并对其进行预处理,所述的视差影像包括相互间具有观测视差的参考影像和待匹配影像;2)对视差影像进行精确密集匹配处理,获得两幅影像的视差偏差值;3)根据视差偏差值计算卫星姿态颤振值,具体为:301)使用频谱分析的方法由步骤2)获得的视差偏差值得到卫星姿态颤振的频率和幅值,当颤振只含有一种频率分量时,执行步骤302),当颤振含有多种频率分量时,执行步骤303);302)通过频谱分析获得视差偏差值的幅值<img file="FDA0000872867190000011.GIF" wi="67" he="77" />和频率ω,并由相位谱获取其初相位<img file="FDA0000872867190000012.GIF" wi="86" he="63" />将t时刻的视差偏差值与影像几何偏移值表示为:<img file="FDA0000872867190000013.GIF" wi="454" he="78" />f(t)=A<sub>f</sub>cos(ωt+f<sub>0</sub>)f(t+Δt)=A<sub>f</sub>cos(ωt+f<sub>0</sub>+Δf)式中,A<sub>f</sub>、f<sub>0</sub>分别为影像几何偏移值的幅值和初相位,Δt为视差影像的成像时间间隔,Δf=ωΔt;根据<img file="FDA0000872867190000014.GIF" wi="444" he="63" />计算得:<img file="FDA0000872867190000015.GIF" wi="509" he="106" /><img file="FDA0000872867190000016.GIF" wi="661" he="198" />从而获得影像几何偏移值f(t);303)通过频谱分析获得视差偏差值各频率分量的幅值<img file="FDA0000872867190000017.GIF" wi="102" he="77" />频率ω<sub>i</sub>和初相位<img file="FDA0000872867190000018.GIF" wi="95" he="63" />影像几何偏移值的各频率等于视差偏差值的频率;对式<img file="FDA0000872867190000019.GIF" wi="414" he="63" />进行离散傅里叶变换,获得:<img file="FDA00008728671900000110.GIF" wi="990" he="71" />式中,<img file="FDA00008728671900000111.GIF" wi="173" he="75" />为虚数单位,<img file="FDA00008728671900000112.GIF" wi="214" he="127" />k为离散傅里叶变换后频域下对应t时刻的虚数,N为参与离散傅里叶变换的数据总数;则<img file="FDA0000872867190000021.GIF" wi="466" he="134" />对f(t)进行频谱分析,获取影像几何偏移值各个频率分量的幅值A<sub>fi</sub>和频率ω<sub>i</sub>,由相位谱获取各个频率分量的初相位f<sub>0i</sub>;从而获得影像几何偏移值f(t):<img file="FDA0000872867190000022.GIF" wi="556" he="135" />M为频率分量个数;304)根据影像几何位移值f(t)计算星体在垂轨方向和沿轨方向的姿态颤振值<img file="FDA0000872867190000023.GIF" wi="180" he="78" />和<img file="FDA0000872867190000024.GIF" wi="204" he="71" /><maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mo>&part;</mo><msub><mi>g</mi><mrow><mi>c</mi><mi>r</mi><mi>o</mi><mi>s</mi><mi>s</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>f</mi><mi>x</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><msub><mi>f</mi><mi>c</mi></msub></mfrac></mrow>]]></math><img file="FDA0000872867190000025.GIF" wi="338" he="142" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mo>&part;</mo><msub><mi>g</mi><mrow><mi>a</mi><mi>l</mi><mi>o</mi><mi>n</mi><mi>g</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>f</mi><mi>y</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><msub><mi>f</mi><mi>c</mi></msub></mfrac></mrow>]]></math><img file="FDA0000872867190000026.GIF" wi="342" he="142" /></maths>式中,f<sub>x</sub>(t)和f<sub>y</sub>(t)分别为f(t)在像平面行方向和列方向的分量,f<sub>c</sub>为相机主距。
地址 200092 上海市杨浦区四平路1239号