一种利用单个光流传感器获取球形电机旋转速度的方法与流程

1.本发明涉及速度传感器与计算机图像处理领域，具体是一种利用单个光流传感器获取球形电机旋转速度的方法。


背景技术：

2.在航天姿态控制、机器人以及数控机床等新兴高科技领域，球形电机凭借其相较于传统电机的优势——可以输出任意轴的力矩、质量更低、体积更小而获得更广泛的应用前景。为了实现对电机高精度的控制，就需要对球形电机进行闭环控制，所以获取球形电机的转速是必要的。
3.球形电机的转速不同于传统电机，测量球形电机的转速不仅需要测量转速大小，还需要测量转速的方向。现在市场上多采用多个摄像头配合获取球形电机的转速大小与方向。主要包括霍尔传感器、光电编码器等等方法，这些方法为了测得多个自由度的球形电机转速，需要多个传感器之间的配合。


技术实现要素：

4.本发明的目的在于提供一种利用单个光流传感器获取球形电机旋转速度的方法，通过确定传感器与球形转子的相对位置，可以实现利用单个摄像头对球形电机的转速测量。
5.本发明的目的可以通过以下技术方案实现：
6.一种利用单个光流传感器获取球形电机旋转速度的方法，包括以下步骤：
7.s1：确定光流传感器与球形电机转子的位置关系；
8.s2：光流摄像头获取球形电机转子表面图像；
9.s3：对球形电机转子图像进行数据处理；
10.s4：球形电机转子图像特征点选取；
11.s5：对每一个特征点进行l
‑
k光流运算；
12.s6：由每个特征点的光流计算完整图像的光流；
13.s7：利用图像的光流获取图像投影速度；
14.s8：利用球形电机与传感器之间的位置关系求解球形转子的速度。
15.进一步的，所述s1中光流传感器安装于球形电机定子外壳上，传感器摄像头参考平面与球形电机转子表面的垂直距离在摄像头工作范围之内。
16.进一步的，所述s2中获取照片的摄像头分辨率在500以上，拍摄时外部光照条件保持不变或变化缓慢。
17.进一步的，所述s3具体包括以下步骤：
18.s31：对rgb图像数据进行灰度计算；
19.s32：将灰度值转化为双精度值，便于计算机处理。
20.进一步的，所述s4具体包括以下步骤：
21.s41：计算图像中每个像素点的梯度矩阵g及其最小特征值λ
m
；
[0022][0023]
式中p
x
，p
y
为计算像素点的位置，w
x
，w
y
为综合选取窗口大小，选取范围在2
‑
7像素范围之内；i
x
，i
y
为像素点灰度值在x，y方向的偏导；
[0024]
s42：寻找图像中所有最小特征值的最大值λ
max
；
[0025]
s43：保留最小特征值大于10％λ
max
的像素点；
[0026]
s44：保留λ
m
局部最大的像素点，保留的像素点的λ
m
要比邻域3*3像素范围内的任意像素点大。
[0027]
进一步的，所述s5具体包括以下步骤：
[0028]
s51：建立图像分层金字塔：
[0029][0030]
式中l为层数，i为某位置像素点的灰度值大小；将获取得到的图像作为金字塔的第0层，从第0层计算得到第1层，再从第1层计算得到第2层；以这样的递归方式得到图像金字塔；式中x、y为图像像素点的位置，图像左上角为(0，0)；
[0031]
s52：最上层特征点的光流值无可靠初始估计值：
[0032][0033]
式中为最大层的光流；
[0034]
s53：迭代计算下一层特征点光流值：
[0035][0036]
式中d
l
为第l层残差位移矢量；
[0037]
s54：若最终收敛，得到最终光流值：
[0038][0039]
式中为该特征像素点光流的最终解，每个特征点都会有计算得到的光流。
[0040]
进一步的，所述s6由每个特征点的光流计算完整图像的光流：
[0041][0042]
式中为第i个特征点的光流，n为特征点的个数，为图像的光流。
[0043]
进一步的，所述s7由光流传感器摄像头的分辨率计算图像光流所代表的投影速
度：
[0044][0045]
式中为投影速度，ppi为摄像头分辨率，δt为前后两帧图像的间隔时间，计算结果单位为cm/s。
[0046]
进一步的，所述s8具体包括以下步骤：
[0047]
s81：摄像头到球形转子摄像面的距离由图像大小计算得到：
[0048][0049][0050]
式中d为摄像头到球形转子摄像面的距离，f为摄像头的焦距，r为球形转子的半径，d
′
为摄像头到球形转子表面的距离，a、b为图像大小；
[0051]
s82：由传感器摄像头与球形转子的位置关系：
[0052][0053]
为实际的转速，单位为cm/s，方向与光流方向相反。
[0054]
本发明的有益效果：
[0055]
本发明的球形转子转速识别方法与传统的球形转子转速识别方法相比，具有无接触的优点，降低损耗；只利用单个摄像头，依靠单个光流摄像头与空间位置关系直接计算速度大小与方向，而其他无接触的识别方法一般都用到了两个或多个摄像头协同工作。
附图说明
[0056]
图1是本发明光流传感器获取球形电机转子照片示意图；
[0057]
图2是本发明球形转子空间坐标系的示意图；
[0058]
图3是本发明转子图像金字塔分层与光流计算示意图；
[0059]
图4是本发明球形电机转子与摄像头相对位置示意图。
具体实施方式
[0060]
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其它实施例，都属于本发明保护的范围。
[0061]
一种利用光流传感器获取球形电机旋转速度的方法，包括以下步骤：
[0062]
s1：球形电机转子照片获取依靠光流传感器上的摄像头，如图1所示，球形电机应用场合不要求高转速输出，转速不高保证了前后两帧图像差异性不大，电机的工作环境一般不发生改变，保证了外部光照环境变换是缓慢均匀的。
[0063]
获取的前后两帧2d图像分别为图像i与图像j。
[0064]
s2：图像数据处理的灰度处理依靠灰度函数，将rgb图像转化为灰度图；将灰度图转化为双精度数据，方便计算机处理。灰度函数：
[0065]
e
′
y＝0.299r 0.587g 0.114b
[0066]
式中r，g，b为rgb图像像素点的值，e
′
y为像素点的灰度值；
[0067]
两幅图像i，j的灰度值分别为i(x，y)与j(x，y)，图像左上顶点的像素坐标为[0 0]
t
。n
x
、n
y
分别为图像的宽度与高度。
[0068]
s3：图像特征点选取，包括：
[0069]
s31：计算图像1中每个像素点的梯度矩阵g及其最小特征值λ
m
；
[0070][0071]
式中p
x
，p
y
为计算像素点的位置，w
x
，w
y
为综合选取窗口大小，选取范围在2～7像素范围之内；i
x
，i
y
为像素点灰度值在x、
y
方向的偏导。
[0072]
s32：寻找图像1中所有最小特征值的最大值λ
max
；
[0073]
s33：保留最小特征值大于10％λ
max
的像素点；
[0074]
s34：保留λ
m
局部最大的像素点，保留像素点的λ
m
要比邻域(3*3大小)中任意像素点的大；
[0075]
得到第一幅图像i的特征点为μ＝[μ
x μ
y
]
t
，在第二幅图像j中对应的特征点为v＝[v
x v
y
]
t
。
[0076]
s4：光流运算基于灰度一致性假设，即：
[0077][0078]
光流d使得上述残差公式最小。
[0079]
对特征点进行l
‑
k光流运算，具体包括：
[0080]
s41：如附图3所示，建立图像分层金字塔：
[0081][0082]
为了满足上式，在图像周围定义一些虚拟元素：
[0083]
i
l
‑1(
‑
1，y)＝i
l
‑1(0，y)
[0084]
i
l
‑1(x，
‑
1)＝i
l
‑1(x，0)
[0085][0086][0087][0088]
图像i
l
的宽度与高度的最大整数满足下列两个条件：
[0089][0090][0091]
由此建立起图像金字塔，对于两帧图像就有两个图像金字塔序列，如对于图像大小为640
×
480的图像，第一、二、三层图像大小分别为320
×
240、160
×
120、80
×
60。
[0092]
s42：第l层的灰度一致性假设：
[0093][0094]
式中g
l
为第l层的初始光流，d
l
为残差位移矢量。
[0095]
最上层特征点的光流值无可靠初始估计值：
[0096][0097]
式中为最大层的光流值；
[0098]
迭代计算下一层特征点光流值：
[0099]
g
l
‑1＝2(g
l
d
l
)
[0100]
式中d
l
为第l层残差位移矢量；
[0101]
若最终收敛，得到最终光流值：
[0102]
d＝g0 d0[0103]
式中d为光流的最终解。
[0104]
s43：l
‑
k光流迭代算法是每一层金字塔都要运行一遍的算法，利用中心差分法计算图像梯度：
[0105][0106][0107]
得到光流矢量：
[0108][0109]
式中为使得误差函数最小的最优光流矢量，g为梯度矩阵，为：
[0110][0111]
式中k代表迭代次数。
[0112]
图像光流迭代计算是基于金字塔分层的l
‑
k光流迭代计算。图像金字塔分层一般为3
‑
5层，每层金字塔每个窗口的l
‑
k光流迭代次数一般在5次以内结束。
[0113]
s5：利用位置关系求解球形转子的速度，包括：
[0114]
s51：由每个特征点的光流计算整个图像的光流：
[0115][0116]
式中为第i个特征点的光流，n为特征点的个数，为图像的光流。
[0117]
s52：由光流传感器摄像头的分辨率计算光流所代表的投影速度：
[0118][0119]
式中为投影速度，ppi为摄像头分辨率，δt为前后两帧图像的间隔时间，计算结果单位为cm/s。
[0120]
s53：如附图2所示，摄像头到球形转子摄像面的距离可以由图像大小计算得到：
[0121][0122][0123]
式中d为摄像头到球形转子摄像面的距离，f为摄像头的焦距，r为球形转子的半径，d
′
为摄像头到球形转子表面的距离，a、b为图像大小。
[0124]
s54：如附图4所示，由传感器摄像头与球形转子的位置关系：
[0125][0126]
为实际的转速，单位为cm/s，方向与光流方向相反。
[0127]
在本说明书的描述中，参考术语“一个实施例”、“示例”、“具体示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
[0128]
以上显示和描述了本发明的基本原理、主要特征和本发明的优点。本行业的技术人员应该了解，本发明不受上述实施例的限制，上述实施例和说明书中描述的只是说明本发明的原理，在不脱离本发明精神和范围的前提下，本发明还会有各种变化和改进，这些变化和改进都落入要求保护的本发明范围内。