主权项 |
一种基于机器视觉的空中受油引导控制方法,其特征在于,包括如下步骤:第一,定义加油机、受油机、加油锥、摄像头所在的坐标系,通过各个坐标系之间的变换来描述几何位姿关系,求解加油锥坐标系的中心到加油探头的距离;第二,对摄像头采集图像进行彩色空间到色相、饱和度和明度化处理,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>Σ</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></msubsup><msubsup><mi>Σ</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>Σ</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>Σ</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>></mo><mn>0</mn></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>x</mi><mrow><mi>i</mi><mi>k</mi></mrow></msub><mo>∈</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以控制受油机。 |