发明名称 一种基于机器视觉的空中受油引导控制方法
摘要 本发明属于飞行控制技术领域,具体涉及到一种基于机器视觉引导的空中受油引导控制方法。本发明是针对为拓宽作战范围为提高成功回收的几率,具有加受油功能的飞机。其特征在于,实现空中加油的关键在于准确测量加油锥管和受油机之间的相对位置,以及对受油机实施有效的引导和控制。对于视觉引导系统来说,可以分为视觉导航与引导控制两个主要组成部分。假定加油机进行定高匀速直线平飞,在加油机的加油锥管上安装辅助光学标记,受油机上安装摄像机,采用视觉测量技术实时测量加油锥管的位置和姿态,设计适当的飞行控制率对受油机进行引导,完成加油对接。本发明适应受油任务特点、降低受油任务过程中飞行员操纵负担,大大提高受油任务成功率。
申请公布号 CN104133480B 申请公布日期 2017.02.15
申请号 CN201410154224.8 申请日期 2014.04.17
申请人 中国航空工业集团公司沈阳飞机设计研究所 发明人 范彦铭;张杨;杨宁;邵铮
分类号 G05D1/10(2006.01)I 主分类号 G05D1/10(2006.01)I
代理机构 北京航信高科知识产权代理事务所(普通合伙) 11526 代理人 周良玉
主权项 一种基于机器视觉的空中受油引导控制方法,其特征在于,包括如下步骤:第一,定义加油机、受油机、加油锥、摄像头所在的坐标系,通过各个坐标系之间的变换来描述几何位姿关系,求解加油锥坐标系的中心到加油探头的距离;第二,对摄像头采集图像进行彩色空间到色相、饱和度和明度化处理,H、S、V三个通道的图像,彩色空间简称RGB,色相、饱和度和明度简称HSV,选取其中的H、S通道图像进行给定阈值分割,此时,得到二值化的图像;二值化图像存在噪声,用腐蚀、膨胀形态处理,将噪声去掉,对两个通道的二值图像进行融合,然后检测到LED光点,并为最后加油锥与摄像机相对位姿测量做准备;由于检测对象是加油锥上的红色LED点,通过RGB彩色空间到HSV空间的变换,得到三个通道H、S、V的值,然后对H、S通道的图像进行阈值分割,可以检测到目标点;第三,进行对应点匹配,即匹配标记连续帧图像对应特征点的位置;通过摄像机采集到加油锥特征点P<sub>j</sub>在图像上的投影点坐标,假定投影点集为<img file="FDA0001154138430000011.GIF" wi="171" he="69" />从摄像机采集图像提取特征点得到的观测点集为{p<sub>1</sub>,p<sub>2</sub>,…p<sub>n</sub>},其中p<sub>j</sub>=(u<sub>j</sub>,v<sub>j</sub>)是图像坐标,而对应的投影点集为<img file="FDA0001154138430000012.GIF" wi="309" he="69" />其中<img file="FDA0001154138430000013.GIF" wi="244" he="63" />是利用针孔投影成像模型计算得到的图像坐标;对于点集p<sub>j</sub>,<img file="FDA0001154138430000014.GIF" wi="83" he="63" />计算两者之间的欧式距离矩阵Err,矩阵大小为m×n:<maths num="0001"><math><![CDATA[<mrow><mi>E</mi><mi>r</mi><mi>r</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><mi>d</mi><mrow><mo>(</mo><msub><mover><mi>p</mi><mo>^</mo></mover><mn>1</mn></msub><mo>,</mo><msub><mi>p</mi><mn>1</mn></msub><mo>)</mo></mrow></mrow></mtd><mtd><mrow><mi>d</mi><mrow><mo>(</mo><msub><mover><mi>p</mi><mo>^</mo></mover><mn>1</mn></msub><mo>,</mo><msub><mi>p</mi><mn>2</mn></msub><mo>)</mo></mrow></mrow></mtd><mtd><mn>...</mn></mtd><mtd><mrow><mi>d</mi><mrow><mo>(</mo><msub><mover><mi>p</mi><mo>^</mo></mover><mn>1</mn></msub><mo>,</mo><msub><mi>p</mi><mi>n</mi></msub><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mi>d</mi><mrow><mo>(</mo><msub><mover><mi>p</mi><mo>^</mo></mover><mn>2</mn></msub><mo>,</mo><msub><mi>p</mi><mn>1</mn></msub><mo>)</mo></mrow></mrow></mtd><mtd><mrow><mi>d</mi><mrow><mo>(</mo><msub><mover><mi>p</mi><mo>^</mo></mover><mn>2</mn></msub><mo>,</mo><msub><mi>p</mi><mn>2</mn></msub><mo>)</mo></mrow></mrow></mtd><mtd><mn>...</mn></mtd><mtd><mrow><mi>d</mi><mrow><mo>(</mo><msub><mover><mi>p</mi><mo>^</mo></mover><mn>2</mn></msub><mo>,</mo><msub><mi>p</mi><mi>n</mi></msub><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mrow><mi>d</mi><mrow><mo>(</mo><msub><mover><mi>p</mi><mo>^</mo></mover><mi>m</mi></msub><mo>,</mo><msub><mi>p</mi><mn>1</mn></msub><mo>)</mo></mrow></mrow></mtd><mtd><mrow><mi>d</mi><mrow><mo>(</mo><msub><mover><mi>p</mi><mo>^</mo></mover><mi>m</mi></msub><mo>,</mo><msub><mi>p</mi><mn>2</mn></msub><mo>)</mo></mrow></mrow></mtd><mtd><mn>...</mn></mtd><mtd><mrow><mi>d</mi><mrow><mo>(</mo><msub><mover><mi>p</mi><mo>^</mo></mover><mi>m</mi></msub><mo>,</mo><msub><mi>p</mi><mi>n</mi></msub><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0001154138430000015.GIF" wi="1124" he="356" /></maths>点集对应数学模型如下:假设有两个点集A,B,s<sub>ik</sub>表示集合A中的点i与集合B中的点k之间的欧式距离,要建立A、B两个点集之间的一一对应的关系,即转化为如下公式:<maths num="0002"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mi>m</mi><mi>i</mi><mi>n</mi><msubsup><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></msubsup><msubsup><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></msubsup><msub><mi>s</mi><mrow><mi>i</mi><mi>k</mi></mrow></msub><msub><mi>x</mi><mrow><mi>i</mi><mi>k</mi></mrow></msub></mrow></mtd></mtr><mtr><mtd><mtable><mtr><mtd><mrow><msubsup><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></msubsup><msub><mi>x</mi><mrow><mi>i</mi><mi>k</mi></mrow></msub><mo>=</mo><mn>1</mn></mrow></mtd><mtd><mrow><mi>k</mi><mo>=</mo><mn>1</mn><mo>...</mo><mi>M</mi></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></msubsup><msub><mi>x</mi><mrow><mi>i</mi><mi>k</mi></mrow></msub><mo>=</mo><mn>1</mn></mrow></mtd><mtd><mrow><mi>i</mi><mo>=</mo><mn>1</mn><mo>...</mo><mi>N</mi></mrow></mtd></mtr></mtable></mtd></mtr><mtr><mtd><mrow><msub><mi>s</mi><mrow><mi>i</mi><mi>k</mi></mrow></msub><mo>&gt;</mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>x</mi><mrow><mi>i</mi><mi>k</mi></mrow></msub><mo>&Element;</mo><mo>{</mo><mn>0</mn><mo>,</mo><mn>1</mn><mo>}</mo></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0001154138430000021.GIF" wi="630" he="534" /></maths>其中S={s<sub>ik</sub>|s<sub>ik</sub>=dist(a<sub>i</sub>,b<sub>k</sub>),a<sub>i</sub>∈A,b<sub>k</sub>∈B,i=1...N,k=1...M};测量得到LED点在图像上的位置信息,然后与加油锥上的标记点对应起来,对应点标记可转化为数学指派问题来描述;P表示LED点在加油锥上的位置:P={p<sub>1</sub>,…,p<sub>N</sub>}A表示点集P通过旋转、缩放矩阵R和转移向量t以及摄像机投影变换到图像坐标系的点集:A={a<sub>1</sub>,...,a<sub>N</sub>}<img file="FDA0001154138430000022.GIF" wi="253" he="70" />B表示图像检测得到的特征点集:B={b<sub>1</sub>,...,b<sub>M</sub>}C表示点集A与点集B对应点匹配,定义如下:C={(a<sub>i</sub>,b<sub>j</sub>)|a<sub>i</sub>∈A,b<sub>j</sub>∈B}第四,在特征点匹配后,进行位姿计算,加油锥图像特征点2D到3D的转化,通过旋转、缩放矩阵R和平移向量T的几何变化,获得加油锥到摄像机坐标系的相对距离信息,然后转化得到相对于受油机坐标系的位置信息[x<sub>drogue</sub> y<sub>drogue</sub> z<sub>drogue</sub>];已知受油机加油管在受油机坐标系的位置信息[x<sub>probe</sub> y<sub>probe</sub> z<sub>probe</sub>],那么空中加油的控制输入为:<maths num="0003"><math><![CDATA[<mrow><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>e</mi><mi>x</mi></msub></mtd></mtr><mtr><mtd><msub><mi>e</mi><mi>y</mi></msub></mtd></mtr><mtr><mtd><msub><mi>e</mi><mi>z</mi></msub></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><msub><mi>x</mi><mrow><mi>d</mi><mi>r</mi><mi>o</mi><mi>g</mi><mi>u</mi><mi>e</mi></mrow></msub><mo>-</mo><msub><mi>x</mi><mrow><mi>p</mi><mi>r</mi><mi>o</mi><mi>b</mi><mi>e</mi></mrow></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>y</mi><mrow><mi>d</mi><mi>r</mi><mi>o</mi><mi>g</mi><mi>u</mi><mi>e</mi></mrow></msub><mo>-</mo><msub><mi>y</mi><mrow><mi>p</mi><mi>r</mi><mi>o</mi><mi>b</mi><mi>e</mi></mrow></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>z</mi><mrow><mi>d</mi><mi>r</mi><mi>o</mi><mi>g</mi><mi>u</mi><mi>e</mi></mrow></msub><mo>-</mo><msub><mi>z</mi><mrow><mi>p</mi><mi>r</mi><mi>o</mi><mi>b</mi><mi>e</mi></mrow></msub></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0001154138430000031.GIF" wi="670" he="339" /></maths>有GLSDC和LHM两种算法来估计摄像机坐标系下加油锥的位置<img file="FDA0001154138430000032.GIF" wi="200" he="81" />和Euler姿态角<img file="FDA0001154138430000033.GIF" wi="318" he="79" />在机器视觉测量系统获取加油机上加油锥的姿态和位置信息后传给受油机的控制系统;控制系统根据受油机自身的姿态、位置信息与从机器视觉测量系统获取的加油锥的信息进行结算,得到受油机与加油锥之间的相对距离矢量r<sub>d</sub>;r<sub>d</sub>与期望相对距离r相减得到误差信号e;e通过控制器模块产生控制指令u以控制受油机。
地址 110035 辽宁省沈阳市皇姑区塔湾街40号
您可能感兴趣的专利