CN104392014B - 一种动态星模拟器及其星图修正方法 - Google Patents
一种动态星模拟器及其星图修正方法 Download PDFInfo
- Publication number
- CN104392014B CN104392014B CN201410370367.2A CN201410370367A CN104392014B CN 104392014 B CN104392014 B CN 104392014B CN 201410370367 A CN201410370367 A CN 201410370367A CN 104392014 B CN104392014 B CN 104392014B
- Authority
- CN
- China
- Prior art keywords
- star
- star chart
- chart
- grid
- axis
- 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
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公布了一种动态星模拟器及其星图修正方法,动态星模拟器包括星图控制计算机、星图显示***、星图修正软件、星图显示器件、光学***和卧式调整机构;星图修正方法包括星图测试、星图修正和星图复测三个步骤。如果星图复测结果不满足模拟精度要求,则再次进行星图修正;如果星图复测结果满足模拟精度要求,则生成星图修正软件写入星图显示***,配合动态星模拟器,实现无穷远处的高精度星图模拟。
Description
技术领域
本发明涉及航天器标定技术,尤其涉及一种动态星模拟器及其星图修正方法。
背景技术
星模拟器是航天技术领域中星敏感器的地面标定装置。星模拟器按星图显示方式不同,分为动态星模拟器和静态星模拟器。按模拟精度分为甚高精度星模拟器、高精度星模拟器和中等精度星模拟器。
动态星模拟器根据姿轨控制仿真动力学计算机提供的星敏感器坐标轴在惯性坐标系下的指向,由星表数据生成当前时刻星敏感器所能观测到的星图数据,再通过接口及驱动电路在星图显示器件上产生星图。由模拟星点发出的光线经准直光学***汇聚后形成平行光,可在室内有限距离上模拟对真实恒星的观测效果。
由于受到动态星模拟器光学***像差和显示器件安装位置等因素的影响,现有动态星模拟器在显示星图之前,需要一套星图修正软件,以提高星图模拟精度。
发明内容
本发明提供一种动态星模拟器星图修正方法,以提高动态星模拟器的星图模拟精度。
本发明提高一种动态星模拟器星图修正方法,其特征在于,包括:
星图测试部分,包括利用经纬仪测试网格并建立测量模型;
所述星图测试部分,所述网格由星图显示***控制星图显示器件显示,所述网格由64个特征点组成,所述经纬仪测试所述64个特征点的位置信息,利用所述64个特征点的所述位置信息,计算所述动态星模拟器主距f主和实测所述64个特征点之间的星间角距γ测(i,j),并根据国家天文台测试天空中各星点的位置信息,得到所述64个特征点的理论星间角距γ理(i,j),计算所述网格上所述64个特征点的星间角距模拟误差Δγ测(i,j)=γ测(i,j)-γ理(i,j);
星图修正部分,包括根据经纬仪对所述64个特征点位置信息的测量结果,计算所述64个特征点的修正系数、在星图显示区域内对所述修正系数进行分区并拟合修正曲线;
所述星图修正部分,计算所述64个特征点的所述修正系数,所述修正系数分为对所述64个特征点的X轴的修正系数Kx(m)和Y轴的修正系数Ky(m);其中,m为所述网格上任一所述特征点,m∈[1,64];
X轴的所述修正系数Kx(m),由式(1)计算:
式中,X理(m)为网格上任一特征点m在X轴的理论位置;
dpix为星图显示器件的像素尺寸;
f主为光学***的实测主距;
αm为经纬仪实测任一特征点m的方位角度值;
Y轴的所述修正系数Ky(m),由式(2)计算:
式中,Y理(m)为网格上任一特征点m在Y轴的理论位置;
dpix为星图显示器件的像素尺寸;
f主为光学***的实测主距;
βm为经纬仪实测任一特征点m的俯仰角度值。
所述星图显示区域内对所述修正系数分区,是根据所述修正系数Kx(m)和Ky(m)的分布情况,在所述星图显示范围内进行所述Kx(m)和所述Ky(m)的分区;
根据分区情况拟合所述修正系数Kx(m)的修正曲线和所述修正系数Ky(m)的修正曲线,利用MATLAB软件,拟合所述X轴修正系数Kx(m)的所述修正曲线f[K(x)](u),其中u表示X轴的修正系数的分区个数,取u∈[1,∞];拟合所述Y轴修正系数Ky(m)的所述修正曲线f[K(y)](v),v表示Y轴的修正系数的分区个数,取v∈[1,∞];
所述Ky(m)和所述Ky(m)的分区结果,将所述64个特征点在X轴的所述修正系数Kx(m)分区拟合的所述修正曲线f[K(x)](u)和在Y轴的所述修正系数Ky(m)分区拟合的所述修正曲线f[K(y)](v)分别写入所述星图显示***,则所述星图显示区域被分为N个区域,其中,N表示星图显示区域同时被X轴的修正分区和Y轴的修正分区共同分割后,形成的分区数量,取N∈[1,∞];
所述X轴的所述修正曲线f[K(x)](u)和所述Y轴的所述修正曲线f[K(y)](v)分别写入所述星图显示***。
星图复测部分,包括测试修正后的网格和测量显示区域内的任意网格,比较复测星间角距γ复测(i,j)与理论星间角距γ理(i,j)的模拟误差Δγ测(i,j)。
所述测试修正后的所述网格,所述修正后的网格由写入了所述星图修正软 件的所述星图显示***生成,所述修正后的网格上所述64个特征点经过所述修正曲线的调制后,所述位置信息发生变化,根据复测的位置信息,计算所述64个特征点之间的复测星间角距γ复测(i,j),所述复测星间角距γ复测(i,j)与所述理论星间角距γ理(i,j)对比,计算所述复测网格上所述64个特征点的星间角距模拟误差Δγ复测(i,j);
所述复测网格上所述64个特征点的所述星间角距模拟误差Δγ复测(i,j)=γ复测(i,j)-γ理(i,j),先计算所述任意网格上特征点间的理论星间角距γ理(i,j),再测量所述任意网格上特征点间的实测星间角距γ复测(i,j),并与动态星模拟器要求的星图模拟精度Δγ(i,j)对比;当Δγ复测(i,j)>Δγ(i,j)时,应重新计算网格上所述64个特征点的所述修正位置、所述修正系数Kx(m)和Ky(m)、所述修正分区、所述修正曲线f[K(x)](u)和f[K(y)](v);当Δγ复测(i,j)<Δγ(i,j)时,则所述星图修正工作完成,并输出;
本发明提供一种动态星模拟***,包括星图控制计算机、星图显示***、星图修正软件、星图显示器件、光学***和卧式调整机构,其特征在于:还包括权利要求1-5任一所述的动态星模拟器,所述星图控制计算机、所述星图显示***、所述星图修正软件、所述星图显示器件、所述光学***设置在所述卧式调整机构上。
由上述技术方案可知,本发明提供的动态星模拟器及其星图修正方法,动态星模拟器包括星图控制计算机、星图显示***、星图修正软件、星图显示器件、光学***和卧式调整机构;星图修正方法包括星图测试、星图修正和星图复测三个部分,根据星图修正方法生成的星图修正软件写入星图显示***,配合星图控制计算机、星图显示器件、光学***和卧式调整机构,完成动态星模拟器的无穷远处星图模拟,星图修正方法提高了星图的模拟精度。
附图说明
图1为本发明实施例提供的一种动态星模拟器的星图修正流程图。
图2为星图显示***显示的网格Ⅰ。
图3为星图显示***显示的网格Ⅱ。
图4为本发明实施例提供的一种动态星模拟器的组成结构示意图。
图5为星图显示***显示的网格Ⅲ。
图6为网格Ⅲ上各特征点在X轴的修正系数分区图。
图7为网格Ⅲ上各特征点在Y轴的修正系数分区图。
图8为网格Ⅲ上各特征点在X轴和Y轴的组合分区形式。
图9为写入星图修正软件的星图显示***显示的网格Ⅳ。
图10为实施例三提供的动态星模拟器。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例,对本发明实施例中的技术方案进行清楚、完整地描述。需要说明的是,在附图或说明书中,相似或相同的元件皆使用相同的附图标记。
实施例一
图1为本发明实施例提供的一种动态星模拟器的星图修正流程图。图2为星图显示***显示的网格Ⅰ。图3为星图显示***显示的网格Ⅱ。如图1、图2和图3所示,本发明实施例提供的星图修正方法具体可分为三个部分:星图测试部分、星图修正部分和星图复测部分。
星图测试部分,首先,利用经纬仪测试动态星模拟器的星图显示视场和星图中心点;其次,利用星图显示软件显示网格Ⅰ,结合星图中心点和网格Ⅰ,测试动态星模拟器主距;然后,根据网格Ⅰ上各特征点的像素位置,计算网格 Ⅰ上各特征点之间的理论星间角距,利用经纬仪测试网格Ⅰ上各特征点的实测星间角距;最后,计算网格Ⅰ上各特征点的星间角距模拟误差。
星图修正部分,首先,根据网格Ⅰ上各特征点的星间角距模拟误差,修正网格Ⅰ上各特征点在X轴和Y轴的像素位置;其次,计算网格Ⅰ上各特征点在X轴和Y轴的修正系数,并根据网格Ⅰ上各特征点在X轴和Y轴的修正系数进行分区;然后,依据网格Ⅰ上各特征点在X轴和Y轴的修正系数的分区结果,拟合各特征点在X轴和Y轴的修正曲线;最后,将拟合的修正曲线写入星图显示软件。
星图复测部分,首先,利用已经写入了星图修正曲线的星图显示软件显示网格Ⅰ,其中,各特征点已经被修正曲线修正;其次,利用经纬仪测试网格Ⅰ上被修正后的各特征点之间的实测星间角距,与其理论星间角距相比,判断是否满足精度要求。如果不能满足精度要求,则返回第二部分,重新计算网格Ⅰ上各特征点的修正位置、修正系数、修正分区和修正曲线;如果满足图精度要求,则利用写入了修正曲线的星图显示***显示网格Ⅱ;然后,根据网格Ⅱ上各特征点的像素位置,计算网格Ⅱ上各特征点之间的理论星间角距,利用经纬仪测试网格Ⅱ上各特征点之间的实测星间角距;最后,比较网格Ⅱ上各特征点之间的实测星间角距与理论星间角距的误差,判断是否满足动态星模拟器的星图模拟精度。如果不能满足星图模拟精度,则返回第二部分,重新计算网格Ⅰ上各特征点的修正位置、修正系数、修正分区和修正曲线;如果满足星图模拟精度,则星图修正完成。
实施例二
图4为本发明实施例提供的一种动态星模拟器的组成结构示意图。图5为星图显示***显示的网格Ⅲ。图6为网格Ⅲ上各特征点在X轴的修正系数分区 图。图7为网格Ⅲ上各特征点在Y轴的修正系数分区图。图8为网格Ⅲ上各特征点在X轴和Y轴的组合分区形式。图9为写入星图修正软件的星图显示***显示的网格Ⅳ。
本发明实施例提供的动态星模拟器主要由星图控制计算机、星图显示***、星图显示器件、光学***和星图修正软件组成。动态星模拟器的工作原理与星图修正过程如下:
首先,星图控制计算机发出指令控制星图显示***,星图显示***控制星图显示器件显示网格Ⅲ,由于星图显示器件位于光学***焦面位置,故网格Ⅲ经过光学***后,以平行光的形式出射,经纬仪接收并测得网格Ⅲ上各特征点的星间角距γ测(i,j),与理论星间角距γ理(i,j)比较,得到星间角距模拟误差Δγ(i,j),其中,i和j分别代表经纬仪测试的网格Ⅲ上任意两个特征点的序号,i∈[1,64],j∈[1,64];
然后,根据星间角距模拟误差Δγ(i,j)=γ测(i,j)-γ理(i,j),利用式1计算网格Ⅲ各特征点X轴的修正系数Kx(m),利用式2计算网格Ⅲ各特征点Y轴的修正系数Ky(m);
其次,根据修正系数Kx(m)和Ky(m)的分布情况,在星图显示范围内进行Kx(m)和Ky(m)的分区;根据分区情况拟合修正系数Kx(m)的修正曲线和修正系数Ky(m)的修正曲线,并写入星图显示***。利用经纬仪测试修正后的像素点,如果测试结果不满足精度要求,则应重新计算修正系数并拟合修正曲线写入星图显示***。如果测试结果满足精度要求,则修正工作完成。
最后,利用写入了星图修正软件的星图显示***显示网格Ⅳ,实测网格Ⅳ上各特征点的星间角距γ测(p,q),计算网格Ⅳ上各特征点的理论星间角距γ理(p,q),得到网格Ⅳ上各特征点的星间角距模拟精度Δγ测(p,q)=γ测(p,q)-γ理(p,q),满足星敏感 器的使用要求。
具体的,星图显示器件的显示像素范围为1920×1920个像素,每个像素尺寸为8μm;网格Ⅲ设有64个特征点,各特征点的像素位置如表1;实测光学***主距f主=39.52mm,计算64个特征点的理论星间角距γ理(i,j)如表2所示;利用经纬仪测试网格Ⅲ上64个特征点的实测星间角距γ实(i,j),其中,取i∈[1,64],取j为显示像素范围的中心点,即0点,并得到模拟误差Δγ(i,j),如表3所示。
表1网格Ⅲ的64个特征点的像素位置
表2网格Ⅲ的64个特征点的理论星间角距值
表3网格Ⅲ的64个特征点的理论星间角距值
具体的,计算网格Ⅲ上各特征点X轴的修正系数Kx(m)的公式,如式(1)所示;计算网格Ⅲ各特征点Y轴的修正系数Ky(m)的公式,如式(2)所示;Kx(m)和Ky(m)的计算结果,如表4所示。其中,m为网格Ⅲ上任意特征点的序号,取m∈[1,64]。
式中,X理(m)为网格Ⅲ上任意特征点m在X轴的理论像素坐标位置;
dpix为星图显示器件的像素尺寸,dpix=8μm;
f主为光学***的实测主距,f主=39.52mm;
αm为经纬仪实测任意像素点m的方位角度值。
式中,Y理(m)为网格Ⅲ上任意特征点m在Y轴的理论像素坐标位置;
dpix为星图显示器件的像素尺寸,dpix=8μm;
f主为光学***的实测主距,f主=39.52mm;
βm为经纬仪实测任意像素点m的俯仰角度值。
表4网格Ⅲ的64个特征点的X轴修正系数Kx(m)和Y轴修正系数Ky(m)
特征点序号m | Kx(m) | Ky(m) | 特征点序号m | Kx(m) | Ky(m) |
1 | 0.001775 | 0.008547 | 33 | 0.003834 | 0.000037 |
2 | 0.001715 | -0.000482 | 34 | 0.001969 | -0.001131 |
3 | 0.000843 | -0.000142 | 35 | 0.001176 | -0.000311 |
4 | 0.001601 | 0.001876 | 36 | 0.000881 | -0.000886 |
5 | 0.000866 | 0.000344 | 37 | -0.000114 | 0.001535 |
6 | 0.001148 | -0.001431 | 38 | -0.001085 | -0.000168 |
7 | 0.001482 | 0.002227 | 39 | -0.002639 | 0.000223 |
8 | 0.002307 | 0.000347 | 40 | -0.003043 | 0.000555 |
9 | -0.000437 | -0.000373 | 41 | 0.000509 | -0.002165 |
10 | 0.000052 | -0.000013 | 42 | -0.000686 | -0.001788 |
11 | 0.000551 | 0.001772 | 43 | -0.000391 | -0.000568 |
12 | 0.000263 | 0.002408 | 44 | -0.000353 | -0.000984 |
13 | 0.000725 | 0.000680 | 45 | -0.000678 | 0.000350 |
14 | 0.001238 | -0.000862 | 46 | -0.000625 | -0.000585 |
15 | 0.001577 | -0.000368 | 47 | -0.000625 | 0.000036 |
16 | 0.001811 | 0.001024 | 48 | -0.000625 | 0.000068 |
17 | -0.002826 | -0.001147 | 49 | -0.000368 | -0.002262 |
18 | -0.000852 | -0.000335 | 50 | -0.001066 | -0.001751 |
19 | 0.000297 | 0.001118 | 51 | -0.001468 | -0.001227 |
20 | -0.000391 | 0.001259 | 52 | -0.000879 | -0.001383 |
21 | 0.000343 | -0.000031 | 53 | -0.000711 | -0.000655 |
22 | 0.014183 | -0.000048 | 54 | -0.000606 | -0.000549 |
23 | 0.000797 | 0.000638 | 55 | -0.000423 | -0.000492 |
24 | 0.000751 | 0.001177 | 56 | -0.000414 | -0.000216 |
25 | -0.005576 | -0.001395 | 57 | -0.001508 | -0.001772 |
26 | -0.003896 | 0.001441 | 58 | -0.001722 | -0.002332 |
27 | 0.000474 | 0.000128 | 59 | -0.001021 | -0.000110 |
28 | -0.001085 | -0.001067 | 60 | -0.000975 | -0.003093 |
29 | -0.000091 | -0.000547 | 61 | -0.001150 | -0.000932 |
30 | 0.001357 | -0.000084 | 62 | -0.000725 | -0.000565 |
31 | 0.001856 | 0.000298 | 63 | -0.001029 | 0.003066 |
32 | 0.002151 | 0.000455 | 64 | -0.001018 | -0.000695 |
具体的,根据修正系数Kx(m)对星图显示范围进行分区,如图5所示;根据修正系数Ky(m)对星图显示范围进行分区,如图6所示。
具体的,对X轴的修正系数Kx(m)分区后,利用Matlab软件拟合X轴修正系数Kx(m)的修正曲线f[K(x)](u),如表5所示,根据分区结果,u表示X轴的修正系数的分区个数,取u∈[1,5]。
表5 X轴修正系数Kx(m)的修正曲线f【K(x)】(u)
X轴分区序号u | 修正方程 |
1 | f[K(x)](1)=-2.204×10-11y3+2.403×10-8y2-3.4316×10-6y-1.78×10-3 |
2 | f[K(x)](2)=-2.099×10-11y3-3.56×10-8y2-1.4387×10-5y-2.95×10-3 |
3 | f[K(x)](3)=6.9968×10-12y3+1.0012×10-8y2+5.5818×10-6y-2.0125×10-3 |
4 | f[K(x)](4)=1.3994×10-10y3+2.1694×10-7y2+1.022×10-4y+6.0625×10-3 |
5 | f[K(x)](5)=7.8714×10-12y3-1.3072×10-8y2+7.9304×10-6y-2.031×10-3 |
对Y轴的修正系数Ky(m)分区后,利用Matlab软件拟合Y轴修正系数Ky(m)的修正曲线f[K(y)](v),如表6所示,根据分区结果,v表示Y轴的修正系数的分区个数,取v∈[1,9]。
表6 Y轴修正系数Ky(m)的修正曲线f【K(y)】(v)
依据Ky(m)和Ky(m)的分区结果,将X轴的修正曲线f[K(x)](u)和Y轴的修正曲线f[K(y)](v)分别写入星图显示***,则星图显示区域被分为12个区域,每个星图显示区域写入不同的修正曲线,如表7所示。
表7星图显示区域分区与修正曲线
具体的,写入了星图修正软件的星图显示***显示的网格Ⅳ如图7所示,网格Ⅳ上设有100个特征点,特征点序号为s,s∈[65,164];计算100个特征点的理论星间角距、实测星间角距和模拟误差,结果如表8所示,可以看出,写入了星图修正软件的星图显示***显示的网格可以满足星敏感器的使用要求。
表8网格Ⅳ上100个特征点的理论星间角距、实测星间角距和模拟误差
实施例三
图10为实施例三提供的动态星模拟器。如图10所示,动态星模拟器包括,卧式调整机构,设置于卧式调整结构上的星图显示器件,设置于星图显示器件后为其提供星点位置信号和星等模拟信号的星图显示***,设置于星图显示***后为其提供星点位置修正方程的星图修正软件,设置于星图显示器件前用于将星图显示器件模拟的星图成像到无穷远的光学***。
具体的,星图修正软件用于根据实测星点位置的星间角距与恒星理论星间角距,计算星点位置信息,利用MATLAB拟合星点位置修正方程,并将星点位置修正方程写入星图显示***。星图控制软件接收星敏感器提供的当前时刻星图的星点位置信号和星等模拟信号,配合星图修正软件写入的星点修正方程,控制星图显示器件显示当前时刻下的模拟星图,星图模拟精度优于18″,星图模拟视场为22°,星等模拟范围为-2Mv~+6Mv;光学***是一种由八片透镜组成的小畸变、波像差和倍率色差小、弥散斑均匀的高成像质量投影光学***,完成对星图显示器件模拟的星图的高质量无穷远成像。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (2)
1.一种动态星模拟器星图修正方法,其特征在于,包括:
星图测试部分,包括利用经纬仪测试网格并建立测量模型;星图修正部分,包括根据经纬仪对64个特征点位置信息的测量结果,计算所述64个特征点的修正系数、在星图显示区域内对所述修正系数进行分区并拟合修正曲线;星图复测部分,包括测试修正后的网格和测量显示区域内的任意网格,比较复测星间角距γ复测(i,j)与理论星间角距γ理(i,j)的模拟误差Δγ复测(i,j);
所述星图测试部分,所述网格由星图显示***控制星图显示器件显示,所述网格由64个特征点组成,所述经纬仪测试所述64个特征点的位置信息,利用所述64个特征点的所述位置信息,计算所述动态星模拟器主距f主和实测所述64个特征点之间的星间角距γ测(i,j),并根据国家天文台测试天空中各星点的位置信息,得到所述64个特征点的理论星间角距γ理(i,j),计算所述网格上所述64个特征点的星间角距模拟误差Δγ测(i,j)=γ测(i,j)-γ理(i,j);
所述星图修正部分,计算所述64个特征点的所述修正系数,所述修正系数分为对所述64个特征点的X轴的修正系数Kx(m)和Y轴的修正系数Ky(m);其中,m为所述网格上任一所述特征点,m∈[1,64];
所述星图显示区域内对所述修正系数分区,是根据所述修正系数Kx(m)和Ky(m)的分布情况,在所述星图显示范围内进行所述Kx(m)和所述Ky(m)的分区;
根据分区情况拟合所述修正系数Kx(m)的修正曲线和所述修正系数Ky(m)的修正曲线,利用MATLAB软件,拟合所述X轴修正系数Kx(m)的所述修正曲线f[K(x)](u),其中u表示X轴的修正系数的分区个数,取u∈[1,∞);拟合所述Y轴修正系数Ky(m)的所述修正曲线f[K(y)](v),v表示Y轴的修正系数的分区个数,取v∈[1,∞);
所述Ky(m)和所述Ky(m)的分区结果,将所述64个特征点在X轴的所述修正系数Kx(m)分区拟合的所述修正曲线f[K(x)](u)和在Y轴的所述修正系数Ky(m)分区拟合的所述修正曲线f[K(y)](v)分别写入所述星图显示***,则所述星图显示区域被分为N个区域,其中,N表示星图显示区域同时被X轴的修正分区和Y轴的修正分区共同分割后,形成的分区数量,取N∈[1,∞);
所述X轴的所述修正曲线f[K(x)](u)和所述Y轴的所述修正曲线f[K(y)](v)分别写入所述星图显示***;
所述测试修正后的所述网格,所述修正后的网格由写入了所述星图修正软件的所述星图显示***生成,所述修正后的网格上所述64个特征点经过所述修正曲线的调制后,所述位置信息发生变化,复测所述任意网格上特征点间的实测星间角距γ复测(i,j),即所述64个特征点之间的复测星间角距γ复测(i,j),所述复测星间角距γ复测(i,j)与所述理论星间角距γ理(i,j)对比,计算所述复测网格上所述64个特征点的星间角距模拟误差Δγ复测(i,j);
所述复测网格上所述64个特征点的所述星间角距模拟误差Δγ复测(i,j)=γ复测(i,j)-γ理(i,j),与动态星模拟器要求的星图模拟精度Δγ(i,j)对比;当Δγ复测(i,j)>Δγ(i,j)时,应重新计算网格上所述64个特征点的所述修正位置、所述修正系数Kx(m)和Ky(m)、所述修正分区、所述修正曲线f[K(x)](u)和f[K(y)](v);当Δγ复测(i,j)<Δγ(i,j)时,则所述星图修正工作完成,并输出;
X轴的所述修正系数Kx(m),由式(1)计算:
式中,X理(m)为网格上任一特征点m在X轴的理论位置;
dpix为星图显示器件的像素尺寸;
f主为光学***的实测主距;
αm为经纬仪实测任一特征点m的方位角度值;
Y轴的所述修正系数Ky(m),由式(2)计算:
式中,Y理(m)为网格上任一特征点m在Y轴的理论位置;
dpix为星图显示器件的像素尺寸;
f主为光学***的实测主距;
βm为经纬仪实测任一特征点m的俯仰角度值。
2.一种动态星模拟***,包括星图控制计算机、星图显示***、星图修正软件、星图显示器件、光学***和卧式调整机构,其特征在于:所述动态星模拟***采用权利要求1所述的动态星模拟器星图修正方法,所述星图控制计算机、所述星图显示***、所述星图修正软件、所述星图显示器件、所述光学***设置在所述卧式调整机构上。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410370367.2A CN104392014B (zh) | 2014-07-31 | 2014-07-31 | 一种动态星模拟器及其星图修正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410370367.2A CN104392014B (zh) | 2014-07-31 | 2014-07-31 | 一种动态星模拟器及其星图修正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104392014A CN104392014A (zh) | 2015-03-04 |
CN104392014B true CN104392014B (zh) | 2017-07-25 |
Family
ID=52609918
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410370367.2A Active CN104392014B (zh) | 2014-07-31 | 2014-07-31 | 一种动态星模拟器及其星图修正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104392014B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106123925B (zh) * | 2016-08-24 | 2019-01-18 | 长春理工大学 | 一种用于动态星模拟器星图显示器件的焦面调整机构 |
CN110926501B (zh) * | 2019-11-08 | 2022-03-22 | 中国科学院长春光学精密机械与物理研究所 | 一种光测设备的自动标校方法、***及终端设备 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5473746A (en) * | 1993-04-01 | 1995-12-05 | Loral Federal Systems, Company | Interactive graphics computer system for planning star-sensor-based satellite attitude maneuvers |
CN1912547A (zh) * | 2006-08-23 | 2007-02-14 | 北京航空航天大学 | 一种高精度低成本星光模拟器 |
-
2014
- 2014-07-31 CN CN201410370367.2A patent/CN104392014B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5473746A (en) * | 1993-04-01 | 1995-12-05 | Loral Federal Systems, Company | Interactive graphics computer system for planning star-sensor-based satellite attitude maneuvers |
CN1912547A (zh) * | 2006-08-23 | 2007-02-14 | 北京航空航天大学 | 一种高精度低成本星光模拟器 |
Non-Patent Citations (2)
Title |
---|
《一种静态星模拟器的设计与星点位置修正方法》;孙高飞,张国玉,姜会林 等;《激光与光电子学进展》;20111231;第1-5页 * |
《基于TFT-LCD的动态星模拟器星点位置修正方法》;王俣,张国玉,高玉军 等;《中国光学》;20110630;第4卷(第3期);第248-250页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104392014A (zh) | 2015-03-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105528500B (zh) | 一种分米级星载tdi ccd立体测绘相机成像仿真方法和*** | |
CN104154928B (zh) | 一种适用于惯性平台内置式星敏感器的安装误差标定方法 | |
CN103822615B (zh) | 一种多控制点自动提取与聚合的无人机地面目标实时定位方法 | |
CN108221603A (zh) | 一种道路的路面三维信息检测装置、方法及*** | |
CN103344258B (zh) | 光电经纬仪伺服***性能测试装置及测试方法 | |
CN108072319A (zh) | 一种运动平台的快速标定***及标定方法 | |
CN110285827A (zh) | 一种距离约束的摄影测量高精度目标定位方法 | |
CN103868528B (zh) | 靶场光学测量设备姿态测量精度的测量方法 | |
CN108759798A (zh) | 一种高精度航天器精度测量的实现方法 | |
CN108510546B (zh) | 一种适用于图谱及结构信息同步探测***的相机标定方法 | |
CN105222806A (zh) | 一种运载火箭双捷联惯组方位偏差标定装置及方法 | |
CN111537003A (zh) | 一种基于折射面共线的星光大气折射测量校正方法 | |
CN110262048A (zh) | 集成成像光场显示***透镜阵列的校准方法 | |
CN106871929A (zh) | 一种测量星敏感器测量坐标系与棱镜坐标系关系的方法 | |
CN108733711B (zh) | 基于三维gis技术的配电线路空间距离获取方法 | |
CN110501026B (zh) | 基于阵列星点的相机内方位元素标定装置及方法 | |
CN108154535B (zh) | 基于平行光管的摄像机标定方法 | |
CN109887041A (zh) | 一种机械臂控制数字相机摄影中心位置和姿态的方法 | |
CN104392014B (zh) | 一种动态星模拟器及其星图修正方法 | |
CN110751601A (zh) | 一种基于rc光学***的畸变校正方法 | |
CN105627916B (zh) | 一种建立跟踪仪地理坐标系及六自由度测量的方法 | |
CN106289317A (zh) | 一种单镜头数字航测相机的单机检校方法及装置 | |
CN110068313B (zh) | 一种基于投影变换的数字天顶仪定向方法 | |
CN109242909A (zh) | 一种面向高精度二维尺寸测量的线阵相机标定算法 | |
CN109682395B (zh) | 星敏感器动态噪声等效角评估方法及*** |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |