CN104360891A - 可见光图像制导导弹仿真简易目标模拟系统及其模拟方法 - Google Patents

可见光图像制导导弹仿真简易目标模拟系统及其模拟方法 Download PDF

Info

Publication number
CN104360891A
CN104360891A CN201410531204.8A CN201410531204A CN104360891A CN 104360891 A CN104360891 A CN 104360891A CN 201410531204 A CN201410531204 A CN 201410531204A CN 104360891 A CN104360891 A CN 104360891A
Authority
CN
China
Prior art keywords
angle
target
axle
straight line
axle table
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.)
Granted
Application number
CN201410531204.8A
Other languages
English (en)
Other versions
CN104360891B (zh
Inventor
王懿
崔建鹏
张辉
包源
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
No 41 Institute Of Fourth Academy China Aerospace Science & Industry Corp
Original Assignee
No 41 Institute Of Fourth Academy China Aerospace Science & Industry Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by No 41 Institute Of Fourth Academy China Aerospace Science & Industry Corp filed Critical No 41 Institute Of Fourth Academy China Aerospace Science & Industry Corp
Priority to CN201410531204.8A priority Critical patent/CN104360891B/zh
Publication of CN104360891A publication Critical patent/CN104360891A/zh
Application granted granted Critical
Publication of CN104360891B publication Critical patent/CN104360891B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Image Analysis (AREA)

Abstract

本发明提供了一种可见光图像制导导弹仿真简易目标模拟系统,采用投影仪向投影幕布投射目标十字;采用控制计算机控制三轴转台的偏航和俯仰方向角;采用可见光导引头探测目标十字;采用目标显示计算机进行目标模拟:接收外部输入的模拟目标相对于三轴转台中心的偏航方向角和俯仰方向角,依据内部存有的二维映射表计算获得十字目标在投影图像中的位置坐标,依据该位置坐标控制投影仪向投影幕布投射十字目标;采用标定计算机进行标定:建立一个三轴转台的偏航方向角和俯仰方向角与十字目标位置一一对应的二维映射表,二维映射表建立完成之后发送至目标显示计算机。该系统能够进行简易图像生成,并能够保证一定的仿真精度。

Description

可见光图像制导导弹仿真简易目标模拟系统及其模拟方法
技术领域
本发明属于半实物仿真技术领域,可用于替代五轴转台及目标模拟器完成图像制导武器的半实物仿真实验。
背景技术
半实物仿真是精确武器系统研制过程中不可缺少的一个重要的环节,对于采用可见光图像制导的导弹,目前半实物仿真常规的做法是采用五轴转台,利用五轴转台的内三轴模拟弹体姿态,利用外两轴模拟目标的相对位置,并在外两轴上安装可见光目标模拟器,模拟可见光目标及其背景。
虽然采用五轴转台方案可以得到较高的仿真精度高,但由于五轴转台和专用目标模拟器的投资大、建设周期长(需12个月以上),导致很多半实物仿真实验室在型号研制初期无法具备末制导仿真的条件,从而延误了型号的研制进度。
发明内容
有鉴于此,本发明提出了可见光图像制导导弹仿真简易目标模拟系统,能够进行简易图像生成,投资成本低、建设周期短,并且能够保证一定的仿真精度。
为达到上述目的,本发明的技术方案为:该系统由投影仪、投影幕布、摄像头、可见光导引头、三轴转台及其控制计算机、目标显示计算机和标定计算机组成。
投影仪在目标显示计算机的控制下向投影幕布投射投影图像,该投影图像为一个位置确定的目标十字。
控制计算机接收转台控制信号,并根据该转台控制信号控制三轴转台的偏航方向角和俯仰方向角。
摄像头安装于三轴转台的滚转轴上,用于在标定过程中对投影图像进行采集。
可见光导引头安装于三轴转台的滚转轴上,用于对投影幕布上的目标十字进行探测。
目标显示计算机在进行目标模拟时,接收外部输入的模拟目标相对于三轴转台中心的偏航方向角和俯仰方向角,依据内部存有的二维映射表计算获得十字目标在投影图像中的位置坐标,依据该位置坐标控制投影仪向投影幕布投射十字目标。
标定计算机采用设定的标定过程建立一个三轴转台的偏航方向角和俯仰方向角与十字目标位置一一对应的二维映射表,二维映射表建立完成之后发送至目标显示计算机。
进一步地,设定的标定过程具体为:
第一步、将摄像头中心轴固定于三轴转台的滚转轴上;在摄像头镜头中心处设定一个标定十字,将投影图像通过投影仪投射于投影幕布上,投影图像中具有一个目标十字;采用摄像头采集投影图像,并在显示计算机显示摄像头所采集的图像。
第二步、确定二维映射表中三轴转台偏航角和俯仰角度的范围:
A、获取粗略的偏航角度范围[αmin,αmax]以及俯仰角度范围[βmin,βmax],在该范围内,投影幕布处于摄像头的视场中;其中α为偏航角度、β为俯仰角度,以下角标区分不同的偏航角度和俯仰角度。
B、控制三轴转台初始角度为(αminmin)。
C、控制三轴转台的滚转角度分别为0°、90°、180°和270°时,标定计算机发送一个目标控制指令至目标显示计算机,从而控制投影仪向投影幕布投射目标十字,该目标控制指令满足目标十字与标定十字中心重合,则依次对应滚转角度分别为0°、90°、180°和270°时,目标十字在投影图像上的坐标分别为(x1,y1)、(x2,y2)、(x3,y3)、(x4,y4)。
D、获得当前步骤下的(x1,y1)和(x3,y3)组成的直线L1与(x2,y2)和(x4,y4)组成的直线L2;判断直线L1和L2之间有无交点:
若无,则控制三轴转台的偏航角增加角度Kα*Xr,俯仰角增加角度Kβ*Yr,重复本步骤C和步骤D,其中Kα和Kβ为投影图像像素到转台角度的横向变化系数;
若有,获取其交点(Xm1,Ym1)。
E、获取此时(Xm1,Ym1)对应的三轴转台的位置为(αmin1min1)。
F、控制三轴转台的初始角度分别为(αminmax)、(αmaxmin)和(αmaxmax)时,执行步骤C~E,获取对应的三轴转台的位置分别为(αmin2max1)、(αmax1min2)和(αmax2max2)。
G、取(αmin1,βmin1)、(αmin2,βmax1)、(αmax1,βmin2)、(αmax2,βmax2)的最小子集,获得精确的三轴转台偏航角和俯仰角的范围[αminf,αmaxf],[βminf,βmaxf]。
第三步、生成标定数据:
I、将[αminf,αmaxf]平均分成m+1个点αi,其中i=0,1,2,3…m;将俯仰角度范围[βminf,βmaxf]平均分成n+1个点βj(j=0,1,2,3…n);以αi为行,βj为列建立二维映射表。
II、以(αij)作为三轴转台的偏航角和俯仰角,执行步骤C,并计算(x1,y1)和(x3,y3)组成的直线与(x2,y2)和(x4,y4)组成的直线的交点,将该交点写入二维映射表中对应第i行第j列交汇的表项中,其中i取遍0,1,2,3…m,j取遍0,1,2,3…n,最终获得完整的二维映射表。
优选地,Kα和Kβ的计算方法为:
a、将标定十字的中心点对准投影图像的中心点进行采集,此时三轴转台的位置为(α1,β1)。
b、执行步骤C,计算(x1,y1)和(x3,y3)组成的直线与(x2,y2)和(x4,y4)组成直线的交点(Xp1,Yp1)。
c、计算上述交点(Xp1,Yp1)与(x1,y1)的偏差(Xr,Yr)。
d、将标定十字的中心点对准投影图像的中心点进行采集,此时三轴转台的位置为(α22);其中α1与α2之间偏差大于1°并在一定范围之内,β1与β2之间的偏差大于1°并在一定范围之内。
e、执行步骤C,并计算(x1,y1)和(x3,y3)组成的直线与(x2,y2)和(x4,y4)组成的直线的交点(Xp2,Yp2)。
f、计算投影图像像素到转台角度的横向变化系数Kα和Kβ,公式为:Kα=(α21)/(Xp2-Xp1),Kβ=(β21)/(Yp2-Yp1)。
优选地,步骤D中获取交点(Xm1,Ym1)之后,对(Xm1,Ym1)进行如下步骤的优化:
D1、判断(Xm1,Ym1)与投影图像中偏航方向边缘的像素距离是否2*Xr1,若是,则控制三轴转台的偏航角减少Kα*Xr1,执行步骤C,获得由(x1,y1)和(x3,y3)组成的直线与(x2,y2)和(x4,y4)组成的直线的交点(Xm2,Ym2),以(Xm2,Ym2)替代(Xm1,Ym1),执行步骤D2;否则,直接执行步骤D2。
D2、判断像素偏移(Xm1,Ym1)在俯仰方向与相对应的边缘的像素距离是否超过2*Yr1,若是,则控制三轴转台俯仰角减少Kβ*Yr1,执行步骤C,获得由(x1,y1)和(x3,y3)组成的直线与(x2,y2)和(x4,y4)组成的直线的交点(Xm3,Ym3),以(Xm3,Ym3)替代(Xm1,Ym1),结束步骤D;否则直接结束步骤D。
进一步地,目标显示计算机接收外部输入的模拟目标相对于三轴转台中心的偏航角αx和俯仰角βy,依据内部存有的二维映射表计算获得十字目标在投影图像中的位置坐标(x,y)的具体方法为:判断(αxy)是否为二维映射表所记载的角度,若存在i',j';i'=0,1,2,3…m,j'=0,1,2,3…n,使得αx=αi',且βx=βj',则直接在二维映射表中定位获取第i'行第j'列交汇的表项作为(x,y),若不存在符合上述条件的i',j',则采用如下方法计算(αxy)对应的(x,y):
S1、建立坐标系,以偏航方向定义为X坐标轴,向右为正,俯仰方向定义为Y坐标轴,向下为正,以二维表内均匀分布的偏航角度值为X坐标轴刻度值,以二维表内均匀分布的俯仰角度值为Y坐标轴刻度值。
S2、在X轴上取距αx最近的右侧三个相邻已知刻度值αk1、αk2、αk3,取βy在Y轴上距其最近的下侧已知刻度值βk,在二维映射表中查找与(αk1k)、(αk2k)、(αk3k)一一对应的表项的X坐标:Xmk1k、Xmk2k、Xmk3k
S3、以Xmk1k、Xmk2k、Xmk3k三点所在直线Lk为底边,以三轴转台回转中心为顶点构建直角三角形Δk,俯仰角为βk,偏航角为0°时三轴转台回转中心线与直线Lk的交点为Xm
S4、直角三角形Δk的直角边与底边的交点坐标x0以及直角边长度zx通过如下公式计算获得:
x 0 = X mk 1 k + X mk 3 k - X mk 2 k - 2 X mk 2 k - X mk 1 k - X mk 3 k tan ( α k 2 - α 0 ) tan ( α k 1 - α 0 + α k 3 - α k 2 ) - 1
z x = X mk 1 k - x 0 tan ( α k 1 - α 0 )
其中 α 0 = α k 1 - arctan ( ( X mk 2 k - X mk 1 k ) tan ( α k 3 - α k 2 ) - ( X mk 3 k - X mk 2 k ) tan ( α k 2 - α k 1 ) tan ( α k 2 - α k 1 ) tan ( α k 3 - α k 2 ) ( 2 X mk 2 k - X mk 1 k - X mk 3 k ) ) .
S5、则偏航为αx时的像素点横坐标x=zxtan(αx0)+x0
S6、取βy在Y轴上距其最近的下侧三个相邻已知刻度值βp1、βp2、βp3,取αx在X轴上距其最近右侧的已知刻度值αp,在二维映射表中查找与(αpp1)、(αpp2)、(αpp3)一一对应的表项的Y坐标Ympp1、Ympp2和Ympp3
S7、以Ympp1、Ympp2和Ympp3三点所在直线Lp为底边,以三轴转台回转中心为顶点构建直角三角形Δp,偏航角为αp,俯仰角为0°时三轴转台回转中心线与直线Lp的交点为Ym
S8、直角三角形Δp的直角边与底边的交点坐标y0以及直角边长度zy通过如下公式计算获得:
y 0 = Y mpp 1 + Y mpp 3 - Y mpp 2 - 2 Y mpp 2 - Y mpp 1 - Y mpp 3 tan ( β p 2 - β 0 ) tan ( β p 1 - β 0 + β p 3 - β p 2 ) - 1
z y = Y mpp 1 - y 0 tan ( β p 1 - β 0 )
其中 β 0 = β p 1 - arctan ( ( Y mpp 2 - Y mpp 1 ) tan ( β p 3 - β p 2 ) - ( Y mpp 3 - Y mpp 2 ) tan ( β p 2 - β p 1 ) tan ( β p 2 - β p 1 ) tan ( β p 3 - β p 2 ) ( 2 Y mpp 2 - Y mpp 1 - Y mpp 3 ) )
S9、则俯仰角为βy时的像素点纵坐标y为:y=zytan(βy0)+y0
有益效果:
1、本发明所提供的可见光图像制导导弹仿真简易目标模拟系统,能够进行简易图像生成,投资成本低、建设周期短,并且能够保证一定的仿真精度。
2、该系统利用摄像头获取图像的位置信息,可实现标定过程的自动化,大幅度的提高了标定的进度并降低了操作人员的工作强度。标定的结果是一个三轴转台偏航、俯仰角度与投影图像像素偏移之间的二维映射表,通过这个映射表可以得到图像显示位置与角度之间的映射关系,该二维映射表为后续进行目标模拟提供了方便。
3、该系统目标显示计算机接收到外部输入的模拟目标相对于三轴转台中心在偏航、俯仰方向的角度,查询标定建立的二维映射表,得到目标图像在投影幕布上显示的位置,对于二维映射表中没有直接记载的角度,则通过一定的方法进行计算,该方法能够精确找到目标图像在投影幕布上显示的位置,相比传统的插值法更具备较高的精度。
附图说明
图1是本发明的系统组成图;
图2直角三角形Δk
具体实施方式
下面结合附图并举实施例,对本发明进行详细描述。
本发明提供了可见光图像制导导弹仿真简易目标模拟系统,该系统如图1所示,由投影仪、投影幕布、摄像头、三轴转台及其控制计算机、目标显示计算机和标定计算机组成。
投影仪在目标显示计算机的控制下向投影幕布投射目标十字。
控制计算机接收转台控制信号,并根据该信号控制三轴转台的位置,该位置可以用偏航方向角和俯仰方向角来表示。
目标显示计算机在进行目标模拟时,接收外部输入的模拟目标相对于三轴转台中心的偏航方向角和俯仰方向角,依据内部存有的二维映射表计算获得十字目标在投影图像中的位置坐标,依据该位置坐标控制投影仪向投影幕布投射十字目标。
标定计算机采用设定的标定方法建立一个三轴转台位置与十字目标位置一一对应的二维映射表,二维映射表建立完成之后发送至目标显示计算机。
本实施例中的标定计算机的功能是为系统标定以获得二维映射表;系统的标定过程是在三轴转台与投影幕布的相对位置固定后进行的。标定过程只需一次,只有当三轴转台与投影幕布的相对位置发生变化后,才需再次标定。
上述设定的标定过程具体为:
第一步、将摄像头中心轴固定于三轴转台的滚转轴上;在摄像头镜头中心处设定一个标定十字,将投影图像通过投影仪投射于投影幕布上,投影图像中具有一个目标十字;采用摄像头采集投影图像,并在显示计算机显示摄像头所采集的图像;投影图像的像素分辨率为Xmax,Ymax;在实际过程中,可以将目标十字设定为红色、将标定十字设定为白色以示区别。
本实施例中摄像头的放大倍率尽可能设置为高放大倍率,使得通过摄像头观察的图像尺寸是原投影图像尺寸的1/6至1/5。
第二步、确定二维映射表中三轴转台偏航角和俯仰角度的范围。
A、依据三轴转台的具体参数获取粗略的偏航角度范围[αmin,αmax]以及俯仰角度范围[βmin,βmax],在该范围内,投影幕布处于摄像头的视场中。
B、控制三轴转台初始位置为(αminmin),该位置为右下边界位置,其中αmin为偏航角度、βmin为俯仰角度;
C、控制三轴转台的滚转角度分别为0°、90°、180°和270°时,标定计算机发送一个目标控制指令至所述目标显示计算机,从而控制所述投影仪向投影幕布投射目标十字,该目标控制指令满足所述目标十字与标定十字中心重合,则依次对应滚转角度分别为0°、90°、180°和270°时,目标十字在投影图像上的坐标分别为(x1,y1)、(x2,y2)、(x3,y3)、(x4,y4);
D、执行步骤C,获得当前步骤下的(x1,y1)和(x3,y3)组成的直线L1与(x2,y2)和(x4,y4)组成的直线L2;判断直线L1和L2之间有无交点,
若无,则控制三轴转台的偏航角增加角度Kα*Xr,俯仰角增加角度Kβ*Yr,重复本步骤C和步骤D,其中Kα和Kβ为投影图像像素到转台角度的横向变化系数。
若有,获取其交点(Xm1,Ym1)。
计算上述交点(Xm1,Ym1)与(x1,y1)的偏差(Xr1,Yr1)。
本实施例中可以采用如下具体步骤进行变化系数Kα和Kβ的计算:
a、将标定十字的中心点对准投影图像的中心点进行采集,此时三轴转台的位置为(α1,β1),其中α为偏航角度、β为俯仰角度,以下角标区分不同的偏航角度和俯仰角度;
b、控制三轴转台的滚转角度分别为0°、90°、180°和270°时,标定计算机发送一个目标控制指令至目标显示计算机,从而控制投影仪向投影幕布投射目标十字,该目标控制指令满足目标十字与标定十字中心重合,即红色十字与白色十字重合,则依次对应滚转角度分别为0°、90°、180°和270°时,目标十字在投影图像上的坐标分别为(x1,y1)、(x2,y2)、(x3,y3)、(x4,y4);计算(x1,y1)和(x3,y3)组成的直线与(x2,y2)和(x4,y4)组成直线的交点(Xp1,Yp1)该交点就是滚转轴与投影图像交点的像素偏移,从而便于后续消除摄像安装误差带来的影响。
c、计算上述交点(Xp1,Yp1)与(x1,y1)的偏差(Xr,Yr),该偏差即为摄像头中心与旋转轴中心的偏差。
d、将标定十字的中心点对准投影图像的中心点进行采集,此时三轴转台的位置为(α22),其中α2为偏航方向角、β2为俯仰方向角;其中α1与α2之间偏差大于1°并且在一定范围之内,β1与β2之间的偏差大于1°并且在一定范围之内,其中在本实施例中,只要α1与α2之间偏差和β1与β2之间的偏差稍微大于1°即可,可以控制在1°~3°的范围之内。
g、重复步骤b,并计算(x1,y1)和(x3,y3)组成的直线与(x2,y2)和(x4,y4)组成的直线的交点(Xp2,Yp2)。
在获得了(Xm1,Ym1)之后,本实施例中又对(Xm1,Ym1)进行了优化,使(Xm1,Ym1)在偏航方向和俯仰方向更靠近投影图像边缘;
优化方法具体为:
D1、判断(Xm1,Ym1)与投影图像中偏航方向边缘的像素距离是否超过2*Xr1,若是,则控制三轴转台的偏航角减少Kα*Xr1,执行步骤C,获得由(x1,y1)和(x3,y3)组成的直线与(x2,y2)和(x4,y4)组成的直线的交点(Xm2,Ym2),以(Xm2,Ym2)替代(Xm1,Ym1)执行步骤D2;否则,直接执行步骤D2。
此步骤的作用是使(Xm1,Ym1)在偏航方向更靠近边缘。
D2、判断像素偏移(Xm1,Ym1)在俯仰方向与相对应的边缘的像素距离是否超过2*Yr1,若是,则控制三轴转台俯仰角减少Kβ*Yr1,执行步骤C,获得由(x1,y1)和(x3,y3)组成的直线与(x2,y2)和(x4,y4)组成的直线的交点(Xm3,Ym3),以(Xm3,Ym3)替代(Xm1,Ym1),结束步骤D;否则直接结束步骤D。
此优化的作用是使(Xm1,Ym1)在俯仰方向更靠近边缘。
E、获取此时(Xm1,Ym1)对应的三轴转台的位置为(αmin1min1);
F、分别控制三轴转台转动到粗略计算的其余三个边界位置,包括右上边界位置(αminmax)、左下边界位置(αmaxmin)和左上边界位置(αmaxmax),执行步骤C~E,获取对应的三轴转台的位置分别为(αmin2max1)、(αmax1min2)和(αmax2max2)。
G、取(αmin1,βmin1)、(αmin2,βmax1)、(αmax1,βmin2)、(αmax2,βmax2)的最小子集,获得精确的三轴转台偏航角和俯仰角的范围[αminf,αmaxf],[βminf,βmaxf];
第三步、生成标定数据;
I、将[αminf,αmaxf]平均分成m+1个点αi(i=0,1,2,3…m);将俯仰角度范围[βminf,βmaxf]平均分成n+1个点βj(j=0,1,2,3…n);以αi为行,βj为列建立二维映射表;
II、以(αij)作为三轴转台的偏航角和俯仰角,执行步骤b,并计算(x1,y1)和(x3,y3)组成的直线与(x2,y2)和(x4,y4)组成的直线的交点(Xmij,Ymij),将(Xmij,Ymij)写入二维映射表中对应第i行第j列的表项中,最终获得完整的二维映射表。
本实施例中,目标显示计算机接收外部输入的模拟目标相对于三轴转台中心的偏航角αx和俯仰角βy,依据内部存有的二维映射表计算获得十字目标在投影图像中的位置坐标(x,y),首先判断(αxy)是否为二维映射表所记载的角度,若存在i',j';(i'=0,1,2,3…m,j'=0,1,2,3…n),使得αx=αi',且βx=βj',则直接在二维映射表中定位获取第i'行第j'列交汇的表项(Xmi'j',Ymi'j')作为(x,y),若不存在符合上述条件的i',j',则采用如下方法计算x和y的值:
S1、建立坐标系,以偏航方向定义为X坐标轴,向右为正,俯仰方向定义为Y坐标轴,向下为正,以二维表内均匀分布的偏航角度值为X坐标轴刻度值,以二维表内均匀分布的俯仰角度值为Y坐标轴刻度值;
S2、在X轴上取距αx最近的右侧三相邻已知刻度值αk1,αk2,αk3(k2-k1=k3-k2=1),取βy在Y轴上距其最近的下侧已知刻度值βk,在二维映射表中查找与(αk1k)、(αk2k)、(αk3k)一一对应的表项的X坐标:Xmk1k、Xmk2k、Xmk3k
S3、以Xmk1k、Xmk2k、Xmk3k三点所在直线Lk为底边,以三轴转台回转中心为顶点构建直角三角形Δk,俯仰角为βk,偏航角为0°时三轴转台回转中心线与直线Lk的交点为Xm
S4、直角三角形Δk如图2所示,的直角边与底边的交点坐标x0以及直角边长度zx通过如下公式计算获得:
x 0 = X mk 1 k + X mk 3 k - X mk 2 k - 2 X mk 2 k - X mk 1 k - X mk 3 k tan ( α k 2 - α 0 ) tan ( α k 1 - α 0 + α k 3 - α k 2 ) - 1
z x = X mk 1 k - x 0 tan ( α k 1 - α 0 )
其中 α 0 = α k 1 - arctan ( ( X mk 2 k - X mk 1 k ) tan ( α k 3 - α k 2 ) - ( X mk 3 k - X mk 2 k ) tan ( α k 2 - α k 1 ) tan ( α k 2 - α k 1 ) tan ( α k 3 - α k 2 ) ( 2 X mk 2 k - X mk 1 k - X mk 3 k ) ) ;
S5、则偏航为αx时的像素点横坐标x=zxtan(αx0)+x0
S6、取βy在Y轴上距其最近的下侧三相邻已知刻度值βp1、βp2、βp3(p2-p1=p3-p2=1),取αx在X轴上距其最近右侧的已知刻度值αp,在二维映射表中查找与(αpp1)、(αpp2)、(αpp3)一一对应的表项的Y坐标Ympp1、Ympp2和Ympp3
S7、以Ympp1、Ympp2和Ympp3三点所在直线Lp为底边,以三轴转台回转中心为顶点构建直角三角形Δp,偏航角为αp,俯仰角为0°时三轴转台回转中心线与直线Lp的交点为Ym
S8、直角三角形Δp的直角边与底边的交点坐标y0以及直角边长度zy通过如下公式计算获得:
y 0 = Y mpp 1 + Y mpp 3 - Y mpp 2 - 2 Y mpp 2 - Y mpp 1 - Y mpp 3 tan ( β p 2 - β 0 ) tan ( β p 1 - β 0 + β p 3 - β p 2 ) - 1
z y = Y mpp 1 - y 0 tan ( β p 1 - β 0 )
其中 β 0 = β p 1 - arctan ( ( Y mpp 2 - Y mpp 1 ) tan ( β p 3 - β p 2 ) - ( Y mpp 3 - Y mpp 2 ) tan ( β p 2 - β p 1 ) tan ( β p 2 - β p 1 ) tan ( β p 3 - β p 2 ) ( 2 Y mpp 2 - Y mpp 1 - Y mpp 3 ) )
S9、则俯仰角为βy时的像素点纵坐标y为:y=zytan(βy0)+y0
综上所述,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (5)

1.可见光图像制导导弹仿真简易目标模拟系统,其特征在于,该系统由投影仪、投影幕布、摄像头、可见光导引头、三轴转台及其控制计算机、目标显示计算机和标定计算机组成; 
所述投影仪在所述目标显示计算机的控制下向投影幕布投射投影图像,该投影图像为一个位置确定的目标十字; 
所述控制计算机接收转台控制信号,并根据该转台控制信号控制所述三轴转台的偏航方向角和俯仰方向角; 
所述摄像头安装于所述三轴转台的滚转轴上,用于在标定过程中对投影图像进行采集; 
所述可见光导引头安装于所述三轴转台的滚转轴上,用于对投影幕布上的目标十字进行探测; 
所述目标显示计算机在进行目标模拟时,接收外部输入的模拟目标相对于所述三轴转台中心的偏航方向角和俯仰方向角,依据内部存有的二维映射表计算获得十字目标在所述投影图像中的位置坐标,依据该位置坐标控制所述投影仪向投影幕布投射十字目标; 
所述标定计算机采用设定的标定过程建立一个三轴转台的偏航方向角和俯仰方向角与十字目标位置一一对应的二维映射表,所述二维映射表建立完成之后发送至所述目标显示计算机。 
2.如权利要求1所述的可见光图像制导导弹仿真简易目标模拟系统,其特征在于,所述设定的标定过程具体为: 
第一步、将摄像头中心轴固定于三轴转台的滚转轴上;在摄像头镜头中心处设定一个标定十字,将投影图像通过投影仪投射于投影幕布上,所述投影图像中具有一个目标十字;采用摄像头采集投影图像,并在显示计算机显示摄像 头所采集的图像; 
第二步、确定二维映射表中三轴转台偏航角和俯仰角度的范围 
A、获取粗略的偏航角度范围[αmin,αmax]以及俯仰角度范围[βmin,βmax],在该范围内,投影幕布处于摄像头的视场中;其中α为偏航角度、β为俯仰角度,以下角标区分不同的偏航角度和俯仰角度; 
B、控制三轴转台初始角度为(αminmin); 
C、控制三轴转台的滚转角度分别为0°、90°、180°和270°时,标定计算机发送一个目标控制指令至所述目标显示计算机,从而控制所述投影仪向投影幕布投射目标十字,该目标控制指令满足所述目标十字与标定十字中心重合,则依次对应滚转角度分别为0°、90°、180°和270°时,目标十字在投影图像上的坐标分别为(x1,y1)、(x2,y2)、(x3,y3)、(x4,y4); 
D、获得当前步骤下的(x1,y1)和(x3,y3)组成的直线L1与(x2,y2)和(x4,y4)组成的直线L2;判断直线L1和L2之间有无交点: 
若无,则控制三轴转台的偏航角增加角度Kα*Xr,俯仰角增加角度Kβ*Yr,重复本步骤C和步骤D,其中Kα和Kβ为投影图像像素到转台角度的横向变化系数; 
若有,获取其交点(Xm1,Ym1); 
E、获取此时(Xm1,Ym1)对应的三轴转台的位置为(αmin1min1); 
F、控制三轴转台的初始角度分别为(αminmax)、(αmaxmin)和(αmaxmax)时,执行步骤C~E,获取对应的三轴转台的位置分别为(αmin2max1)、(αmax1min2)和(αmax2max2); 
G、取(αmin1,βmin1)、(αmin2,βmax1)、(αmax1,βmin2)、(αmax2,βmax2)的最小子集,获得精确的三轴转台偏航角和俯仰角的范围[αminf,αmaxf],[βminf,βmaxf]; 
第三步、生成标定数据; 
I、将[αminf,αmaxf]平均分成m+1个点αi,其中i=0,1,2,3…m;将俯仰角度范围[βminf,βmaxf]平均分成n+1个点βj(j=0,1,2,3…n);以αi为行,βj为列建立二维映射表; 
II、以(αij)作为三轴转台的偏航角和俯仰角,执行步骤C,并计算(x1,y1)和(x3,y3)组成的直线与(x2,y2)和(x4,y4)组成的直线的交点,将该交点写入所述二维映射表中对应第i行第j列交汇的表项中,其中i取遍0,1,2,3…m,j取遍0,1,2,3…n,最终获得完整的二维映射表。 
3.如权利要求2所述的可见光图像制导导弹仿真简易目标模拟系统,其特征在于,Kα和Kβ的计算方法为: 
a、将所述标定十字的中心点对准所述投影图像的中心点进行采集,此时三轴转台的位置为(α1,β1); 
b、执行步骤C,计算(x1,y1)和(x3,y3)组成的直线与(x2,y2)和(x4,y4)组成直线的交点(Xp1,Yp1); 
c、计算上述交点(Xp1,Yp1)与(x1,y1)的偏差(Xr,Yr); 
d、将所述标定十字的中心点对准所述投影图像的中心点进行采集,此时三轴转台的位置为(α22);其中α1与α2之间偏差大于1°并在一定范围之内,β1与β2之间的偏差大于1°并在一定范围之内; 
e、执行步骤C,并计算(x1,y1)和(x3,y3)组成的直线与(x2,y2)和(x4,y4)组成的直线的交点(Xp2,Yp2); 
f、计算投影图像像素到转台角度的横向变化系数Kα和Kβ,公式为:Kα=(α21)/(Xp2-Xp1),Kβ=(β21)/(Yp2-Yp1)。 
4.如权利要求2所述的可见光图像制导导弹仿真简易目标模拟系统,其特征在于,所述步骤D中获取交点(Xm1,Ym1)之后,对(Xm1,Ym1)进行如下步骤的优化: 
D1、判断(Xm1,Ym1)与投影图像中偏航方向边缘的像素距离是否2*Xr1,若是, 则控制三轴转台的偏航角减少Kα*Xr1,执行步骤C,获得由(x1,y1)和(x3,y3)组成的直线与(x2,y2)和(x4,y4)组成的直线的交点(Xm2,Ym2),以(Xm2,Ym2)替代(Xm1,Ym1),执行步骤D2;否则,直接执行步骤D2; 
D2、判断像素偏移(Xm1,Ym1)在俯仰方向与相对应的边缘的像素距离是否超过2*Yr1,若是,则控制三轴转台俯仰角减少Kβ*Yr1,执行步骤C,获得由(x1,y1)和(x3,y3)组成的直线与(x2,y2)和(x4,y4)组成的直线的交点(Xm3,Ym3),以(Xm3,Ym3)替代(Xm1,Ym1),结束步骤D;否则直接结束步骤D。 
5.如权利要求2所述的可见光图像制导导弹仿真简易目标模拟系统,其特征在于,所述目标显示计算机接收外部输入的模拟目标相对于所述三轴转台中心的偏航角αx和俯仰角βy,依据内部存有的二维映射表计算获得十字目标在所述投影图像中的位置坐标(x,y)的具体方法为:判断(αxy)是否为二维映射表所记载的角度,若存在i’,j’,i’=0,1,2,3…m,j’=0,1,2,3…n,使得αx=αi',且βx=βj',则直接在二维映射表中定位获取第i'行第j'列交汇的表项作为(x,y),若不存在符合上述条件的i',j',则采用如下方法计算(αxy)对应的(x,y): 
S1、建立坐标系,以偏航方向定义为X坐标轴,向右为正,俯仰方向定义为Y坐标轴,向下为正,以二维表内均匀分布的偏航角度值为X坐标轴刻度值,以二维表内均匀分布的俯仰角度值为Y坐标轴刻度值; 
S2、在X轴上取距αx最近的右侧三个相邻已知刻度值αk1、αk2、αk3,取βy在Y轴上距其最近的下侧已知刻度值βk,在二维映射表中查找与(αk1k)、(αk2k)、(αk3k)一一对应的表项的X坐标:Xmk1k、Xmk2k、Xmk3k; 
S3、以Xmk1k、Xmk2k、Xmk3k三点所在直线Lk为底边,以三轴转台回转中心为顶点构建直角三角形Δk,俯仰角为βk,偏航角为0°时三轴转台回转中心线与直线Lk的交点为Xm; 
S4、所述直角三角形Δk的直角边与底边的交点坐标x0以及直角边长度zx通过如下公式计算获得: 
其中
S5、则偏航为αx时的像素点横坐标x=zxtan(αx0)+x0; 
S6、取βy在Y轴上距其最近的下侧三个相邻已知刻度值βp1、βp2、βp3,取αx在X轴上距其最近右侧的已知刻度值αp,在二维映射表中查找与(αpp1)、(αpp2)、(αpp3)一一对应的表项的Y坐标Ympp1、Ympp2和Ympp3; 
S7、以Ympp1、Ympp2和Ympp3三点所在直线Lp为底边,以三轴转台回转中心为顶点构建直角三角形Δp,偏航角为αp,俯仰角为0°时三轴转台回转中心线与直线Lp的交点为Ym; 
S8、所述直角三角形Δp的直角边与底边的交点坐标y0以及直角边长度zy通过如下公式计算获得: 
其中
S9、则俯仰角为βy时的像素点纵坐标y为:y=zytan(βy0)+y0。 
CN201410531204.8A 2014-10-10 2014-10-10 可见光图像制导导弹仿真简易目标模拟系统及其模拟方法 Active CN104360891B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410531204.8A CN104360891B (zh) 2014-10-10 2014-10-10 可见光图像制导导弹仿真简易目标模拟系统及其模拟方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410531204.8A CN104360891B (zh) 2014-10-10 2014-10-10 可见光图像制导导弹仿真简易目标模拟系统及其模拟方法

Publications (2)

Publication Number Publication Date
CN104360891A true CN104360891A (zh) 2015-02-18
CN104360891B CN104360891B (zh) 2018-01-09

Family

ID=52528155

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410531204.8A Active CN104360891B (zh) 2014-10-10 2014-10-10 可见光图像制导导弹仿真简易目标模拟系统及其模拟方法

Country Status (1)

Country Link
CN (1) CN104360891B (zh)

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105371701A (zh) * 2015-12-01 2016-03-02 中国科学院长春光学精密机械与物理研究所 电视制导导引头动态模拟器
CN106767161A (zh) * 2016-12-09 2017-05-31 北京仿真中心 一种红外与可见光复合制导半实物仿真方法
CN106839877A (zh) * 2015-12-07 2017-06-13 上海机电工程研究所 一种光学多波段复合目标模拟系统的误差补偿方法
CN108008642A (zh) * 2017-11-29 2018-05-08 北京航空航天大学 一种导弹制导系统半实物仿真方法
CN108843912A (zh) * 2018-06-20 2018-11-20 石家庄硕华电子科技有限公司 一种导引头伺服系统的低反射传动机构
CN109932925A (zh) * 2017-12-15 2019-06-25 北京机电工程研究所 多侦察点仿真方法
CN110285711A (zh) * 2019-05-14 2019-09-27 中国人民解放军海军工程大学 基于信息物理系统的制导弹药飞行姿态修正模拟系统及方法
CN111025935A (zh) * 2019-12-23 2020-04-17 北京世纪联信科技有限公司 一种基于机器视觉的抛射物落点实景仿真系统及方法
CN111080708A (zh) * 2019-11-07 2020-04-28 中国人民解放军国防科技大学 捷联光学成像导引头标定方法
CN111459046A (zh) * 2020-02-20 2020-07-28 南京理工大学 一种图像导引头用目标及场景实时动态生成系统与方法
CN112629342A (zh) * 2020-10-30 2021-04-09 西北工业大学 弹体姿态角测量方法
CN114894047A (zh) * 2021-09-17 2022-08-12 中国人民解放军63875部队 一种靶场关键段姿态交会测量优化布站分析方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103954179A (zh) * 2014-04-30 2014-07-30 北京理工大学 捷联红外导引头隔离度寄生回路评估系统

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103954179A (zh) * 2014-04-30 2014-07-30 北京理工大学 捷联红外导引头隔离度寄生回路评估系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
袁志勇: "《新型光学目标模拟器的原理与结构》", 《红外与激光工程》 *

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105371701B (zh) * 2015-12-01 2017-02-01 中国科学院长春光学精密机械与物理研究所 电视制导导引头动态模拟器
CN105371701A (zh) * 2015-12-01 2016-03-02 中国科学院长春光学精密机械与物理研究所 电视制导导引头动态模拟器
CN106839877A (zh) * 2015-12-07 2017-06-13 上海机电工程研究所 一种光学多波段复合目标模拟系统的误差补偿方法
CN106767161A (zh) * 2016-12-09 2017-05-31 北京仿真中心 一种红外与可见光复合制导半实物仿真方法
CN106767161B (zh) * 2016-12-09 2019-05-03 北京仿真中心 一种红外与可见光复合制导半实物仿真方法
CN108008642A (zh) * 2017-11-29 2018-05-08 北京航空航天大学 一种导弹制导系统半实物仿真方法
CN109932925B (zh) * 2017-12-15 2021-12-10 北京机电工程研究所 多侦察点仿真方法
CN109932925A (zh) * 2017-12-15 2019-06-25 北京机电工程研究所 多侦察点仿真方法
CN108843912A (zh) * 2018-06-20 2018-11-20 石家庄硕华电子科技有限公司 一种导引头伺服系统的低反射传动机构
CN110285711A (zh) * 2019-05-14 2019-09-27 中国人民解放军海军工程大学 基于信息物理系统的制导弹药飞行姿态修正模拟系统及方法
CN111080708A (zh) * 2019-11-07 2020-04-28 中国人民解放军国防科技大学 捷联光学成像导引头标定方法
CN111025935A (zh) * 2019-12-23 2020-04-17 北京世纪联信科技有限公司 一种基于机器视觉的抛射物落点实景仿真系统及方法
CN111459046A (zh) * 2020-02-20 2020-07-28 南京理工大学 一种图像导引头用目标及场景实时动态生成系统与方法
CN112629342A (zh) * 2020-10-30 2021-04-09 西北工业大学 弹体姿态角测量方法
CN112629342B (zh) * 2020-10-30 2023-10-03 西北工业大学 弹体姿态角测量方法
CN114894047A (zh) * 2021-09-17 2022-08-12 中国人民解放军63875部队 一种靶场关键段姿态交会测量优化布站分析方法
CN114894047B (zh) * 2021-09-17 2023-10-03 中国人民解放军63875部队 一种靶场关键段姿态交会测量优化布站分析方法

Also Published As

Publication number Publication date
CN104360891B (zh) 2018-01-09

Similar Documents

Publication Publication Date Title
CN104360891A (zh) 可见光图像制导导弹仿真简易目标模拟系统及其模拟方法
CN102914262B (zh) 一种基于附加视距的非合作目标贴近测量方法
CN105067011A (zh) 一种基于视觉标定及坐标转换的测量系统整体校准方法
CN101539397B (zh) 物体三维姿态的精密光学测量方法
CN103604411B (zh) 基于图像识别的经纬仪自动准直方法
CN106681170B (zh) 一种半实物制导仿真方法及仿真系统
CN101566476B (zh) 基于6自由度机械臂的景像匹配半实物仿真系统
CN108507403B (zh) 基于智能化光电标定技术的自行高炮多轴线一致性检测装置
CN113340277B (zh) 一种基于无人机倾斜摄影的高精度定位方法
CN103376097B (zh) 自动布局和点转换系统
CN103017686A (zh) 一种利用激光跟踪仪调整主镜光轴与水平轴垂直的方法
CN104990533B (zh) 卫星地面物理仿真系统超高精度姿态测量方法及装置
CN104932535A (zh) 利用仿真系统对机载前视红外搜索设备闭环测试的方法
CN104914877A (zh) 一种轻型无人机自主感知与规避方法
CN103345073A (zh) 单透镜光学光轴定位装置及方法
CN104880200A (zh) 复合制导系统初始姿态现场校准系统及方法
CN203857967U (zh) 复合制导系统初始姿态现场校准系统
CN104535078A (zh) 一种基于标志点的光电设备对飞行目标的测量方法
CN105424060B (zh) 一种飞行器星敏感器与捷联惯组安装误差的测量方法
CN202885828U (zh) 一种激光影像定位导航装置
CN101858712B (zh) 一种应用于飞机总装配的航炮数字化校靶方法
CN203216389U (zh) 应用于火炮射弹初速测量装置的对准装置
CN103791837B (zh) 一种跨区域高精度空间坐标基准体系
CN104260908B (zh) 一种跨平台对地观测卫星联合演示验证系统
CN102798391B (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