-
摘要: 研究陶瓷晶粒尺寸分布对估计陶瓷样品的物理属性具有重要意义, 当前主要依赖人工方法测量晶粒尺寸, 由于晶粒形状不规则且大小不一, 因此人工方法测量效率低、误差大. 针对该问题, 提出一种数据与模型联合驱动的陶瓷材料晶粒分割算法. 该算法首先通过图像预处理解决材料表面反光导致的灰度不均匀问题; 其次利用本文提出的鲁棒分水岭变换实现图像中晶粒的预分割, 解决传统分水岭算法存在的过分割以及分割区域个数与轮廓精度难以平衡的问题; 最后提出轻量级富卷积特征网络输出晶粒轮廓, 并利用该轮廓对预分割结果进行优化. 与主流图像分割算法相比, 该算法一方面利用鲁棒分水岭变换实现了更为准确的晶粒区域定位, 另一方面利用图像的低层与高层特征融合获取了更为精准的晶粒轮廓. 实验结果表明, 该算法不仅能够实现陶瓷材料晶粒尺寸的精准计算, 而且具有较高的计算效率, 为分析陶瓷材料物理属性提供了客观准确的数据.Abstract: The research of the distribution of ceramic grain sizes is significant for estimating physical properties of ceramic materials. At present, researchers mainly use manual methods to measure ceramic grains, which leads to low efficiency and large error due to the irregular shape and different sizes of grains in ceramic images. To solve this issue, a new algorithm used for grain segmentation of ceramic materials is proposed in this paper. This algorithm first performs image pre-processing to solve the problem of intensity inhomogeneity caused by surface reflection of materials in the imaging process. Secondly, a robust watershed transform (RWT) is presented and used for pre-segmentation of grains in scanning electron microscope images, which solves two problems existing in traditional watershed transform, one is over-segmentation and the other is the difficulty of balancing the number of segmentation regions and the accuracy of boundaries. Finally, we present a lightweight and richer convolutional features network (LRCF) used for grain contour prediction, and use the results generated by LRCF to optimize the pre-segmentation result. Compared with popular image segmentation algorithms, the proposed algorithm shows two advantages. On one hand, it achieves more accurate grain region segmentation due to the employment of RWT; on the other hand, it provides more accurate grain boundaries due to the optimization of LRCF. Experimental results show that the proposed algorithm not only provides accurate grain sizes of ceramic materials, but also shows high calculation efficiency, which is significant for property analysis of ceramic materials.
-
机器人从自由度上可分为非冗余机器人和冗余机器人[1]. 非冗余机器人是指关节自由度小于或等于任务空间维度的机器人. 不同于非冗余机器人, 冗余机器人的关节自由度大于任务空间维度, 其额外的关节自由度提升了控制过程中的灵活性, 使机器人能够实现障碍物躲避等其他附加任务的能力[2-4]. 现存的大多数冗余机器人控制方案都需要准确的机器人模型信息, 如机器人出厂标定的DH (Denavit-Hartenberg)参数[5]. 它通过4个参数确定关节与连杆间的转化关系, 辅助科研人员实现机器人的精确建模. 文献[6]提出一种复杂工作空间下的机器人实时规划方法, 使机器人在跟踪目标运动的同时, 实现了躲避障碍物的功能. 为提高移动机器人搬运物品的效率, 文献[7]提出一种基于双向采样的方案, 用于为给定的移动机器人移动拾取和运输线路生成期望轨迹. 实验表明该方案能够降低35%的移动时间并提升16倍的计算速度. 一种正交投影机器人重复运动方案在文献[8]中提出, 解决了现有的重复运动方案存在的因关节角度误差而导致的位置误差问题. 上述冗余机器人的控制方案都是在假定机器人模型已知的情况下, 设计并构建控制方案驱动机器人完成目标任务. 然而, 机器人的模型信息通常会因实际使用过程中配备不同类型的末端执行器而发生改变. 这会导致此类方案控制精度降低, 甚至导致任务的失败.
数据驱动技术自提出以来就受到了国内外的高度重视[9-10]. 其主要思想是通过收集并利用系统运行过程中大量的在线、离线数据, 实现对系统有效、鲁棒、稳定的建模与控制[11-12]. 数据驱动技术的提出有效地解决了模型不确定或者模型发生改变时的系统控制问题. 已有学者将数据驱动的思想应用于机器人控制领域, 并取得了卓越的成果. 针对机器人在实际使用过程中出现的模型不确定问题, 文献[13]提出一种基于数据驱动的位置/力控制方案, 从运动学的角度解决了模型未知的机器人位置/力控制问题. 文献[14]提出一种加速度层数据驱动的重复运动规划方案, 并构造了相应的递归神经网络, 对模型未知的冗余机器人提供了一种新的加速度层解决方案. 大量研究表明数据驱动技术能够对模型未知的机器人进行有效控制, 具有较强的应用价值.
冗余机器人的轨迹跟踪控制同样是一个重要的研究方向. 机器人面临的大多数任务需要控制机器人的末端执行器沿预定的轨迹运动. 基于速度层、加速度层的算法被提出用于求解冗余机器人的最优解. 值得指出的是, 目前大多数从运动学和优化角度建模的方法存在两方面的不足[15-17]. 一方面, 这些方法基于机器人结构信息已知的情况实现对机器人的控制. 例如, 文献[18]基于冗余机器人的模型构建了一个动态神经网络方法, 实现了良好的控制性能. 然而, 结构信息不准确或未知的机器人是可能存在的. 例如, 物流机器人可能会受到货物形状、重量等因素的影响, 导致其结构信息不准确. 同样地, 自研机器人的设计和制造过程中也可能因人为疏忽或技术限制等原因导致其结构信息错误. 因此, 面向结构信息不准确或未知的机器人, 上述基于准确结构信息的方法将无法被应用. 另一方面, 这些方法仅考虑机器人末端执行器的位置跟踪, 忽略了其姿态控制的重要性. 例如, 文献[14]构建了一个可应用于结构信息未知机器人的重复运动控制方案, 但其仅考虑了末端执行器的位置跟踪, 在一定程度上限制了其应用场景. 事实上, 机器人末端执行器的姿态控制在大多数应用场景下都是不可或缺的, 在机器人执行任务的同时, 还需对其末端执行器姿态加以控制. 例如医疗手术、工业制造等领域[19-20]. 为弥补现有方法以上两方面的不足, 本文构建了基于数据驱动的冗余机器人末端执行器位姿控制方案. 该方案从机器人的运动学出发, 用四元数表示了末端执行器的姿态, 实现了对末端执行器位姿的控制. 此外, 一种神经动力学求解器被设计用于方案求解. 相关的理论分析与仿真实验均验证了方案的有效性与可行性, 与现有冗余机器人位姿控制方案的比较结果体现了所提方案的新颖性.
1. 问题描述
本节从冗余机器人的运动学出发, 给出机器人末端执行器位姿与关节角间的关系, 并推导了四元数的表示方法以控制末端执行器的姿态. 此外, 设计了估计雅可比矩阵的数据驱动算法. 值得注意的是, 所提机器人末端执行器位姿控制方案所涉及到的符号含义在表1 作具体解释.
表 1 所提冗余机器人控制方案的符号含义Table 1 Definitions of variables of the proposed scheme for redundant manipulators符号 含义 $ {{\boldsymbol{\theta}}} \in {\bf{R}}^a $ 机器人关节角向量 $ \dot{\boldsymbol{\theta}}\in {\bf{R}}^a $ 机器人关节角速度向量 $ \dot{\boldsymbol{\theta}}^{-}(\dot{\boldsymbol{\theta}}^{+}) $ 关节角速度的下界(上界) $ {\boldsymbol r}\in {\bf{R}}^b $ 末端执行器的位置向量 $ \boldsymbol{r}^{d}\in {\bf{R}}^b $ 末端执行器的期望位置向量 $ \dot{\boldsymbol r}\in {\bf{R}}^b $ 末端执行器的速度向量 $ \dot{\hat{\boldsymbol r}}\in {\bf{R}}^b $ 末端执行器的估计速度向量 $ f(\cdot): {\bf{R}}^a \rightarrow {\bf{R}}^b $ 机器人非线性前向运动学映射 $ J=\dfrac{\partial f({{\boldsymbol{\theta}}})}{\partial {{\boldsymbol{\theta}}}}\in {\bf{R}}^{b\times a} $ 机器人雅可比矩阵 $ \hat{J}\in {\bf{R}}^{b\times a} $ 机器人估计雅可比矩阵 $ {\dot{\hat{J}}}\in {\bf{R}}^{b\times a} $ 机器人估计雅可比矩阵的导数 $ M(\boldsymbol \theta)\in {\bf{R}}^{3\times 3} $ 末端执行器的方向旋转矩阵 $ {\boldsymbol q}_{E}(\boldsymbol \theta)\in {\bf{R}}^{4} $ 末端执行器的方向四元数 $ \boldsymbol{\overline{o}}(\boldsymbol \theta)\in {\bf{R}}^{5} $ 末端执行器的方向向量 $ \tilde{\boldsymbol q}\in {\bf{R}}^{5} $ 末端执行器的期望方向向量 $ H({\boldsymbol \theta})=\dfrac{\partial{\boldsymbol q}_{E}(\boldsymbol \theta)}{\partial{\boldsymbol \theta}}\in {\bf{R}}^{4\times a} $ $ {\boldsymbol q}_{E} $ 的雅可比矩阵 $ G({\boldsymbol{\theta}})=\dfrac{\partial{\boldsymbol{\overline{o}}({\boldsymbol{\theta}}})}{\partial{{\boldsymbol{\theta}}}}\in {\bf{R}}^{5\times a} $ $ \boldsymbol{\overline{o}}({\boldsymbol{\theta}}) $的雅可比矩阵 $ \kappa(\boldsymbol q)=\dfrac{\partial{{\tilde{\boldsymbol q}}}}{\partial{\boldsymbol q}}\in {\bf{R}}^{5\times 4} $ $ \tilde{\boldsymbol q} $ 的雅可比矩阵 $ \boldsymbol{u}\in {\bf{R}}^a $ 方差为极小值的独立同分布零均值随机噪声 ${\boldsymbol{u} }_{0}\in {\bf{R} }^a$ $ \boldsymbol{u} $的上界 $ \hat{\dot{{\boldsymbol{\theta}}}}\in {\bf{R}}^a $ 受噪声驱动的关节角速度 $ \Vert \cdot \Vert_2 $ 向量的二范数 $ \mathrm{tr(\cdot)} $ 矩阵的迹 1.1 基础知识
冗余机器人的冗余度解析在其控制过程中十分重要. 冗余机器人的$ a $维自由度能够通过非线性映射函数转化到$ b $维的末端笛卡尔坐标系空间, 具体的数学表达式为
$$ \begin{array}{*{20}{l}} f({{\boldsymbol{\theta}}})={{\boldsymbol{r}}}\end{array} $$ (1) 进一步, 对式(1)关于时间求导, 可得
$$ \begin{array}{*{20}{l}} J\dot{{\boldsymbol{\theta}}}=\dot{{\boldsymbol{r}}}\end{array} $$ (2) 冗余机器人的末端执行器在任务空间中通常需要保持特定角度, 因此需要对末端执行器的位姿加以控制. 机器人末端执行器的姿态信息通常用旋转矩阵$ M $表示, 即
$$ \begin{array}{*{20}{l}} M({\boldsymbol{\theta}})= \begin{bmatrix} a_{11} & a_{12} & a_{13} \\ a_{21} & a_{22} & a_{23} \\ a_{31} & a_{32} & a_{33} \end{bmatrix} \end{array} $$ (3) 其中, $a_{ij}\;(i,j\in{1,2,3})$表示矩阵$ M $中的第$ i $行$ j $列个元素. 四元数由实数加上三个虚数单位组成, 主要用于在三维空间中表示旋转操作. 末端执行器的四元数可以看作是将基坐标系绕向量旋转得到的一个新的坐标系. 因此, 本文使用四元数表示末端执行器的姿态. 定义一个归一化四元数为
$$ {{\boldsymbol{q}}(\phi,\boldsymbol{n})=\left[\cos\frac{\phi}{2},\sin\frac{\phi}{2}{\boldsymbol{n}}\right]=[\phi,v_{x},v_{y},v_{z}]^\text{T}} $$ (4) 其中, $ \phi $表示旋转角度; $ {\boldsymbol{n}}$表示旋转轴向量. 由式(4)可得
$$ \left\{ \begin{split} &\phi=\cos\frac{\phi}{2}, \quad v_{x}=\sin\frac{\phi}{2}{n}_{x} \\ &v_{y}=\sin\frac{\phi}{2}{n}_{y}, \quad v_{z}=\sin\frac{\phi}{2}{n}_{z} \end{split} \right.$$ (5) 式(5)中四元数的复数形式可以表示为${\boldsymbol{q}}=\phi \;+ v_{x}{{\rm{i}}}+ v_{y}{\rm{j}} + v_{z}{\rm{k}}$. 符号${{\rm{i}}}$, ${{\rm{j}}}$, ${{\rm{k}}}$分别表示x轴, y轴, z轴方向的虚数单位. 根据四元数的定义, 可得
$$ \begin{array}{*{20}{l}} ||{\boldsymbol{q}}||_{2}=v_{x}^{2}+v_{y}^{2}+v_{z}^{2}+\phi^{2}=1 \end{array}$$ (6) 1.2 四元数与姿态控制
对于基坐标系的向量$ \boldsymbol s=[s_{x},s_{y},s_{z}]^\mathrm{T} $, 可以通过旋转运算写为四元数的形式$ {\tilde{{\boldsymbol{s}}}=[0,s_{x},s_{y},s_{z}]^\mathrm{T}}$. 在此基础上, 变换后的坐标可以表示为${ \boldsymbol s}_{w}={\boldsymbol q}\;\cdot {\boldsymbol s} \cdot {\boldsymbol q}^{\ast} .$ 其中“$ \cdot $”表示四元数的乘法, $ \boldsymbol q^{\ast} $为$ \boldsymbol q $的共轭形式. 给定四元数$ \boldsymbol q $ 的旋转矩阵排列可以表示为
$$ \begin{split} & R({{\boldsymbol{q}}})=\\ &\quad\begin{bmatrix} 1-2v_{y}^{2}-2v_{z}^{2} \quad & 2v_{x}v_{y}-2v_{z}\phi \quad & 2v_{x}v_{z}+2v_{y}\phi\\ 2v_{x}v_{y}+2v_{z}\phi \quad & 1-2v_{x}^{2}-2v_{z}^{2} \quad & 2v_{y}v_{z}-2v_{x}\phi \\ 2v_{x}v_{z}+2v_{y}\phi \quad & 2v_{y}v_{z}+2v_{x}\phi \quad & 1-2v_{x}^{2}-2v_{y}^{2} \end{bmatrix} \end{split} $$ (7) 联立式(3), (6), (7), 机器人末端执行器的方向可以用四元数$ {{\boldsymbol{q}}_{E}(\boldsymbol \theta)=[\phi^{e},v_{x}^{e},v_{y}^{e},v_{z}^{e}]^\text{T}}$表示为
$$ \left\{ \begin{aligned} &\phi^{e}({\boldsymbol{\theta}}) =\frac{\sqrt{{\text{tr}}(\tilde{{\boldsymbol{q}}})+1}}{2}\\ &v_{x}^{e}({\boldsymbol{\theta}}) = \frac{a_{32}-a_{23}}{4\phi} \\ & v_{y}^{e}({\boldsymbol{\theta}}) = \frac{a_{13}-a_{31}}{4\phi} \\ & v_{z}^{e}({\boldsymbol{\theta}}) = \frac{a_{21}-a_{12}}{4\phi} \end{aligned} \right. $$ (8) 对式(8)关于时间求导, 可得运动学控制式为
$$ {\dot{{\boldsymbol{q}}}^{e}=\frac{\partial{\boldsymbol{q}}_{E}(\boldsymbol \theta)}{\partial{\boldsymbol{\theta}}}\frac{\partial{\boldsymbol{\theta}}}{\partial t}=H({\boldsymbol{\theta}}){\dot{\boldsymbol{\theta}}}} $$ (9) 在此基础上, 机器人末端执行器的方向由旋转矩阵$ M(\boldsymbol \theta) $控制, 进一步可以推广到四元数的实时跟踪即旋转矩阵$ R({\boldsymbol{q}}) $的收敛问题. 由于旋转矩阵$ M(\boldsymbol \theta) $为正交矩阵, 其行、列均为正交的单位向量, 旋转矩阵$ M(\boldsymbol \theta) $中任意2 阶子矩阵即可唯一确定旋转矩阵$ M(\boldsymbol \theta) $. 此外需要一个额外元素来确认正负, 方向向量$ \boldsymbol{\overline{o}}(\boldsymbol \theta) $可以表示为
$$ \begin{array}{*{20}{l}} {\boldsymbol{\overline{o}}}(\boldsymbol \theta)=[a_{11}, a_{12}, a_{21}, a_{22},a_{33}]^\mathrm{T}\in {\bf{R}}^5 \end{array}$$ (10) 对于给定的四元数$ \boldsymbol q $, 其旋转矩阵$ R(\boldsymbol q) $对应的方向向量期望值$ \tilde{{\boldsymbol{q}}}$可以表示为
$$ \begin{array}{*{20}{l}} {\tilde {{\boldsymbol{q}}}}=\begin{bmatrix} 1-2v_{y}^{2}-2v_{z}^{2} \\ 2v_{x}v_{y}-2v_{z}\phi \\ 2v_{x}v_{y}+2v_{z}\phi \\ 1-2v_{x}^{2}-2v_{z}^{2} \\ 1-2v_{x}^{2}-2v_{y}^{2} \end{bmatrix}\in {\bf{R}}^5 \end{array} $$ (11) 为跟踪所需的方向向量, 设计误差函数为${\boldsymbol{e}}= {\overline{{\boldsymbol{o}}}}({\boldsymbol{\theta}})-{\tilde{\boldsymbol{q}}}$, 对该式求导可得
$$ G({\boldsymbol{\theta}}){\dot{{\boldsymbol{\theta}}}=\kappa(\boldsymbol q){\dot{\boldsymbol{q}}}} $$ (12) 其中, $ \dot{{\boldsymbol{q}}}$为向量$ \boldsymbol q $关于时间的导数. 针对不同的任务需求, 冗余机器人的末端执行器方向控制可以通过式(12)实现.
1.3 数据驱动
冗余机器人的控制方案需要使用机器人的雅可比矩阵. 机器人制造商生产机器人时, 为机器人标定了标准的结构参数, 通过这些参数可以计算出机器人的雅可比矩阵. 然而, 机器人的结构信息通常会因实际使用过程中装配不同类型的末端执行器而发生改变, 并且机器人结构参数的测量是复杂且费时的.
数据驱动技术能够将系统参数作为输入量, 动态地更新系统的未知量. 通过数据驱动技术, 机器人能够学习自身结构信息, 并利用学习到的信息协助完成给定的任务. 本文结合数据驱动技术, 提出一种雅可比矩阵实时更新的学习方案以控制模型未知的冗余机器人. 设计误差函数$ \Vert \dot{\hat{{\boldsymbol{r}}}}-\dot{\boldsymbol{r}}\Vert_2^2$, 其中$\dot{\hat{\boldsymbol{r}}}= \hat{J}\dot{{\boldsymbol{\theta}}}$. 通过梯度下降法可以得到
$$ {\dot{\hat{J}}}=-\sigma\frac{\partial{\Vert {\dot{\hat{{\boldsymbol{r}}}}-{\dot{\boldsymbol{r}}}\Vert_2^2}}}{\partial{\hat{J}}}=-\sigma(\hat{J}{\dot{{\boldsymbol{\theta}}}}-{\dot{\boldsymbol{r}}}){\dot{{\boldsymbol{\theta}}}}^\mathrm{T} $$ (13) 其中, $ \sigma>0 $代表控制收敛速率. 值得指出的是, 式(13)中的关节角速度$ \dot{{\boldsymbol{\theta}}}$和末端执行器速度$ \dot{{\boldsymbol{r}}}$为真实的测量数据, 本文以$ \dot{{\boldsymbol{\theta}}}$和$ \dot{{\boldsymbol{r}}}$为核心进行雅可比矩阵的学习. $ \dot{{\boldsymbol{\theta}}}$和$ \dot{{\boldsymbol{r}}}$可通过RGB-D传感器、惯性运动单元、动作捕捉系统等途径获得. 在此基础上, 冗余机器人估计的雅可比矩阵$ \hat J $能够收敛到真实的雅可比矩阵$ J $, 从而实现对机器人的控制.
2. 控制方案的构建与求解
本节首先构建一种基于数据驱动技术的冗余机器人末端执行器位姿控制方案, 其考虑了速度层的控制误差反馈. 随后设计了相应的神经动力学求解器对所提方案进行求解, 能有效控制模型未知的机器人进行末端执行器位置与姿态的同步跟踪. 此外, 相关的理论分析证明了方案的收敛性.
2.1 方案构建
冗余机器人末端执行器的位姿精确度是评价方案的重要因素. 为保证机器人轨迹跟踪任务的顺利执行, 末端执行器通常需要与任务平面保持特定的角度. 本文提出一种基于数据驱动的冗余机器人末端执行器位姿控制方案, 其表述为
$$ \min \;\; \frac{1}{2}\dot{{\boldsymbol{\theta}}}^\mathrm{T}\dot{{\boldsymbol{\theta}}}\;\;\;\tag{14a} $$ $$ \; {\rm{s.t.}} \quad {\hat{J}}{\dot{\hat{\boldsymbol{\theta}}}}={\dot{\boldsymbol{r}}}^{d}\tag{14b} $$ $$\qquad\;\; G({\boldsymbol{\theta}})\dot{\hat{{\boldsymbol{\theta}}}}=\kappa(\boldsymbol q)\dot{{\boldsymbol{q}}}\tag{14c} $$ $$\qquad\;\; {\dot{\hat{J}}}=-\sigma(\hat{J}\hat{\dot{{\boldsymbol{\theta}}}}-{\dot{{\boldsymbol{r}}})\hat{\dot{{\boldsymbol{\theta}}}}^\mathrm{T}}\tag{14d} $$ $$ \qquad\; \; \hat{\dot{{\boldsymbol{\theta}}}}=\dot{{\boldsymbol{\theta}}}+{\boldsymbol{u},\quad 0 \le \left\|\boldsymbol{u}\right\| \le {\boldsymbol{u}}_{0}} \tag{14e\,} $$ $$ \qquad\;\; {\dot{{\boldsymbol{\theta}}}}^-\le \dot{{\boldsymbol{\theta}}}\le {\dot{{\boldsymbol{\theta}}}}^+ \tag{\,14f\,} $$ 值得注意的是, $ \hat{J} $或$ \dot{{\boldsymbol{\theta}}}$的初始值为零时, 系统的输入、输出始终为零, 机器人无法准确地执行轨迹跟踪任务. 为避免该情况的发生, 引入受噪声驱动的关节角速度$ \hat{\dot{{\boldsymbol{\theta}}}}=\dot{{\boldsymbol{\theta}}}+{\boldsymbol{u}}$. 为便于下文的推导, 合并$ \hat{J}\dot{\hat{{\boldsymbol{\theta}}}}=\dot{{\boldsymbol{r}}}$, $G({\boldsymbol{\theta}}){\dot{\hat{{\boldsymbol{\theta}}}}=\kappa(\boldsymbol q){\dot{\boldsymbol{q}}}}$, 并改写为$ A\dot{\hat{{\boldsymbol{\theta}}}}={\dot{{{\boldsymbol{v}}}}}^{d}$, 其中, $A=\Bigg[\begin{aligned}\hat{J}\;\;\;\\G({\boldsymbol{\theta}})\end{aligned}\Bigg] ,$ ${{\boldsymbol{v}}}^{d}=\Bigg[\begin{aligned}\boldsymbol{r}^{d}\\{\tilde {\boldsymbol{q}}}\;\end{aligned}\Bigg] , {\boldsymbol{v}}^{d}$ 表示 ${\boldsymbol{v}}=\Bigg[\begin{aligned}\boldsymbol{r}\;\;\,\\ \boldsymbol{\overline{o}}({\boldsymbol{\theta}})\end{aligned}\Bigg]$的期望值.
2.2 神经动力学求解器设计
设计误差函数为$ E=\Vert {\boldsymbol{v}}-{\boldsymbol{v}}^{d}\Vert^2 $. 由梯度下降法可得, 其沿负梯度下降方向可以表示为
$$ \dot{\hat{{\boldsymbol{\theta}}}}=-\eta\frac{\partial{E}}{\partial{\hat{{\boldsymbol{\theta}}}}}=\eta A^\mathrm{T}({{\boldsymbol{v}}}^{d}-{{\boldsymbol{v}}}) $$ (15) 其中, $ \eta>0 $为设计参数. 为修正式(15)求解过程中的滞后误差问题, 添加补偿项$ {\boldsymbol{\varphi}}$可得
$$ \begin{array}{*{20}{l}} \dot{\hat{{\boldsymbol{\theta}}}}=\eta A^\mathrm{T}({{\boldsymbol{v}}^{d}-{\boldsymbol{v}})+\boldsymbol{\varphi}}\end{array} $$ (16) 对式(16)两边同乘$ A $, 可得
$$ \begin{array}{*{20}{l}} A\dot{\hat{{\boldsymbol{\theta}}}}=\eta AA^\mathrm{T}({{\boldsymbol{v}}^{d}-{\boldsymbol{v}})+A\boldsymbol{\varphi}}\end{array} $$ (17) 当系统趋于稳态时, 应有$ {\boldsymbol{v}}^{d}-{\boldsymbol{v}}=0 $成立. 则式(17)可重写为
$$ \begin{array}{*{20}{l}} A\dot{\hat{{\boldsymbol{\theta}}}}=A{\boldsymbol{\varphi}=\dot{{\boldsymbol{v}}}^{d}}\end{array} $$ (18) 对于式(18), 考虑最小化函数(14a), 并应用最小二乘法可以得到$ {\boldsymbol{\varphi}=A^{\dagger}\dot{{\boldsymbol{v}}}^{d}}$, 其中, $ A^{\dagger}=A^{\mathrm{T}}(AA^{\mathrm{T}})^{-1} $. 为避免接近奇异时求矩阵伪逆结果极大的问题, 设计具有辅助参数的递归过程为
$$\left\{ \begin{aligned} &{\boldsymbol{\varphi}} =A^\mathrm{T}{\boldsymbol{\beta}}\\ &{\dot{{\boldsymbol{\beta}}}} =\mu(\dot{{{\boldsymbol{v}}}}^{d}-AA^{\mathrm{T}}\boldsymbol{\beta}) \end{aligned} \right. $$ (19) 其中, $ \mu>0 $为收敛系数, $ {\boldsymbol{\beta}}$为设计参数. 联立式(16)和式(19), 可得
$$\qquad\qquad\;\; \left\{ \begin{aligned} &\dot{\hat{{\boldsymbol{\theta}}}} =P_{\boldsymbol{\ell}}(\eta A^\mathrm{T}({\boldsymbol{v}}^{d}-{\boldsymbol{v}})+A^{\mathrm{T}}\boldsymbol{\beta}) &&\qquad (20{\rm{a}})\\ &\dot{\boldsymbol{\beta}} =\mu(\dot{{\boldsymbol{v}}}^{d}-AA^{\mathrm{T}}\boldsymbol{\beta}) &&\qquad (20{\rm{b}})\end{aligned} \right. $$ 其中, $ P_{\boldsymbol{\ell}}(\cdot) $表示对机器人施加物理约束的投影函数, 具体地, $ \boldsymbol{\ell}=\{\dot{\hat{{\boldsymbol{\theta}}}} \in {\bf{R}}^{a}, \dot{{\boldsymbol{\theta}}}^{-} \le {\dot{\hat{{\boldsymbol{\theta}}}}} \le \dot{{\boldsymbol{\theta}}}^{+}\} $. 至此, 面向基于数据驱动的冗余机器人末端执行器位姿控制方案的神经动力学求解器为
$$\qquad\qquad\;\;\left\{ \begin{aligned} &{\dot{\hat{{\boldsymbol{\theta}}}}} =P_{\boldsymbol{\ell}}(\eta A^\mathrm{T}({{\boldsymbol{v}}}^{d}-{\boldsymbol{v}})+A^{\mathrm{T}}\boldsymbol{\beta}) &&\qquad (21{\rm{a}})\\ &\dot{\boldsymbol{\beta}}=\mu(\dot{{{\boldsymbol{v}}}}^{d}-AA^{\mathrm{T}}\boldsymbol{\beta}) &&\qquad (21{\rm{b}}) \\ &{\dot{\hat{J}}} =-\sigma(\hat{J}{\hat{\dot{{\boldsymbol{\theta}}}}}-{\dot{{\boldsymbol{r}}}){\hat{\dot{{\boldsymbol{\theta}}}}}^\mathrm{T}}&&\qquad (21{\rm{c}})\\ &{\hat{\dot{{\boldsymbol{\theta}}}}} ={\dot{{\boldsymbol{\theta}}}}+{\boldsymbol{u}},\quad 0 \le \left\|\boldsymbol{u}\right\| \le {\boldsymbol{u}}_{0}&&\qquad (21{\rm{d}}) \end{aligned} \right. $$ 通过神经动力学求解器(21)能够得到机器人执行任务时各关节角速度, 从而实现机器人末端执行器姿态的有效控制.
2.3 理论分析与证明
本节对所提的基于数据驱动的冗余机器人末端执行器位姿控制方案的收敛性进行分析证明.
定理 1. 当$ A^{\mathrm{T}}{\boldsymbol{\beta}}\in \boldsymbol{\ell} $时, 面向基于数据驱动的冗余机器人末端执行器位姿控制方案(式(14))的神经动力学求解器(21)的控制误差$ {\boldsymbol{\epsilon}}={\boldsymbol{v}}-{\boldsymbol{v}}^{d}$全局收敛至$ {{\boldsymbol{0}}} $.
证明. 首先对含有存储状态信息$ {\boldsymbol{\beta}}$的式(21b)进行收敛性分析, 设计李雅普诺夫函数如下:
$$ N_{1}=(\dot{{\boldsymbol{v}}}^{d}-AA^{\mathrm{T}}\boldsymbol{\beta})^{\mathrm{T}}(\dot{{\boldsymbol{v}}}^{d}-AA^{\mathrm{T}}\boldsymbol{\beta})/2 \geq 0 $$ (22) 以参数$ {\boldsymbol{\beta}}$为决策变量, 式(22)关于时间求导, 可得
$$ \begin{split} \dot{N}_{1}=\; &-(\dot{{\boldsymbol{v}}}^{d}-AA^{\mathrm{T}}\boldsymbol{\beta})^{\mathrm{T}}AA^{\mathrm{T}}\dot{\boldsymbol{\beta}}= \\ &\ \mu(\dot{{\boldsymbol{v}}}^{d}-AA^{\mathrm{T}}\boldsymbol{\beta})^{\mathrm{T}}AA^{\mathrm{T}}(\dot{{\boldsymbol{v}}}^{d}-AA^{\mathrm{T}}\boldsymbol{\beta}) \leq \\ &-\eta \boldsymbol{\ell}_{\min}(\dot{{\boldsymbol{v}}}^{d}-AA^{\mathrm{T}}\boldsymbol{\beta})^{\mathrm{T}}(\dot{{\boldsymbol{v}}}^{d}-AA^{\mathrm{T}}\boldsymbol{\beta})=\\ &-\eta\boldsymbol{\ell}_{\min}N_{1}\leq 0 \end{split}$$ (23) 其中, $ {\boldsymbol{\ell}_{\min}}$表示$ AA^{\mathrm{T}} $的最小特征值. 根据拉萨尔不变集原理[21], 当$ \dot{N}_{1} $= 0 时, 参数$ {\boldsymbol{\beta}}$指数收敛至$ (AA^{\mathrm{T}})^{-1}\dot{{{\boldsymbol{v}}}^{d}}$. 将参数$ {\boldsymbol{\beta}}$的稳态值代入(20a), 可得
$$ \begin{split} \dot{\hat{{\boldsymbol{\theta}}}} =\; &P_{\boldsymbol{\ell}}(\eta A^\mathrm{T}({\boldsymbol{v}}^{d}-{\boldsymbol{v}})+A^{\mathrm{T}}(AA^{\mathrm{T}})^{-1}\dot{{\boldsymbol{v}}}^{d})= \\ &P_{\boldsymbol{\ell}}(\eta A^\mathrm{T}({\boldsymbol{v}}^{d}-{\boldsymbol{v}})+A^{\dagger}\dot{{\boldsymbol{v}}}^{d}) \end{split}$$ (24) 控制误差的导数$ \dot{{{\boldsymbol{\epsilon}}}}$可以表示为
$$ \begin{split} \dot{{{\boldsymbol{\epsilon}}}} =\; &A\dot{\hat{{\boldsymbol{\theta}}}}-\dot{{\boldsymbol{v}}}^{d}=\\ & AP_{\boldsymbol{\ell}}(-\eta A^\mathrm{T}{\boldsymbol{\epsilon}}+A^{\dagger}\dot{{\boldsymbol{v}}}^{d})-\dot{{\boldsymbol{v}}}^{d}=\\ &A(P_{\boldsymbol{\ell}}(-\eta A^\mathrm{T}{\boldsymbol{\epsilon}}+A^{\dagger}\dot{{\boldsymbol{v}}}^{d})-A^{\dagger}\dot{{\boldsymbol{v}}}^{d}) \end{split}$$ (25) 为分析控制误差$ {{\boldsymbol{\epsilon}}}$的稳定性, 设计李雅普诺夫函数为$ N_{2}={\boldsymbol{\epsilon}}^{\mathrm{T}}{\boldsymbol{\epsilon}}/2 $, 其时间导数可以表示为
$$ \begin{split} \dot{N}_{2} =\; & {{\boldsymbol{\epsilon}}}^{\mathrm{T}}\dot{{\boldsymbol{\epsilon}}}= \\ & {\boldsymbol{\epsilon}}^{\mathrm{T}}A(P_{\boldsymbol{\ell}}(-\eta A^\mathrm{T}{\boldsymbol{\epsilon}}+A^{\dagger}\dot{{\boldsymbol{v}}}^{d})-A^{\dagger}\dot{{\boldsymbol{v}}}^{d})=\\ &-\frac{1}{\eta}((-\eta A^\mathrm{T}{\boldsymbol{\epsilon}}+A^{\dagger}\dot{{\boldsymbol{v}}}^{d})-A^{\dagger}\dot{{\boldsymbol{v}}}^{d})^{\mathrm{T}}\times \\ &(P_{\boldsymbol{\ell}}(-\eta A^\mathrm{T}{{\boldsymbol{\epsilon}}}+A^{\dagger}\dot{{\boldsymbol{v}}}^{d})-A^{\dagger}\dot{{\boldsymbol{v}}}^{d}) \end{split}$$ (26) 对于投影函数$ P_{\boldsymbol{\ell}} $, 假定其区间下(上)界为${c^{-}}\; ({c^{+}})$, 则当${c^{+}}\in{\boldsymbol{\ell}}$时, 应有$\Vert P_{\boldsymbol{\ell}}({c^{-}})-{c^{-}} \Vert \leq \Vert {c^{+}}-{c^{-}} \Vert$恒成立. 令${c^{-}}=-\eta A^\mathrm{T}{{\boldsymbol{\epsilon}}+A^{\dagger}\dot{{\boldsymbol{v}}}^{d},{c^{+}}=A^{\dagger}\dot{{\boldsymbol{v}}}^{d}}$, 则有以下关系成立
$$ \begin{split} &\Vert P_{\boldsymbol{\ell}}(-\eta A^\mathrm{T}{{\boldsymbol{\epsilon}}+A^{\dagger}\dot{{\boldsymbol{v}}}^{d})+\eta A^\mathrm{T}{\boldsymbol{\epsilon}}-A^{\dagger}\dot{{\boldsymbol{v}}}^{d}}\Vert^{2}\leq \Vert \eta A^\mathrm{T}{{\boldsymbol{\epsilon}}}\Vert^{2} \end{split} $$ (27) 不等式(27)左侧展开, 可得
$$ \begin{split} &\Vert P_{\boldsymbol{\ell}}(-\eta A^\mathrm{T}{{\boldsymbol{\epsilon}}+A^{\dagger}\dot{{\boldsymbol{v}}}^{d})+\eta A^\mathrm{T}{\boldsymbol{\epsilon}}-A^{\dagger}\dot{{\boldsymbol{v}}}^{d}}\Vert^{2}=\\ &\qquad\Vert P_{\boldsymbol{\ell}}(-\eta A^\mathrm{T}{{\boldsymbol{\epsilon}}+A^{\dagger}\dot{{\boldsymbol{v}}}^{d})-A^{\dagger}\dot{{\boldsymbol{v}}}^{d}}\Vert^{2}\;+\\ &\qquad\Vert \eta A^\mathrm{T}{{\boldsymbol{\epsilon}}}\Vert^{2}-2((-\eta A^\mathrm{T}{{\boldsymbol{\epsilon}}}+A^{\dagger}\dot{{\boldsymbol{v}}}^{d})-A^{\dagger}\dot{{\boldsymbol{v}}}^{d})^\mathrm{T}\;\times \\ &\qquad(P_{\boldsymbol{\ell}}(-\eta A^\mathrm{T}{\boldsymbol{\epsilon}}+A^{\dagger}\dot{{\boldsymbol{v}}}^{d})-A^{\dagger}\dot{{\boldsymbol{v}}}^{d}) \\[-10pt] \end{split}$$ (28) 进一步地, 对式(27)和式(28)微分, 可得
$$ \begin{split} &\Vert P_{\boldsymbol{\ell}}(-\eta A^\mathrm{T}{{\boldsymbol{\epsilon}}+A^{\dagger}\dot{{\boldsymbol{v}}}^{d})-A^{\dagger}\dot{{\boldsymbol{v}}}^{d}}\Vert^{2}\leq\\ &\qquad2((-\eta A^\mathrm{T}{{\boldsymbol{\epsilon}}}+A^{\dagger}\dot{{\boldsymbol{v}}}^{d})-A^{\dagger}\dot{{\boldsymbol{v}}}^{d})^\mathrm{T}\;\times \\ &\qquad(P_{\boldsymbol{\ell}}(-\eta A^\mathrm{T}{\boldsymbol{\epsilon}}+A^{\dagger}\dot{{\boldsymbol{v}}}^{d})-A^{\dagger}\dot{{\boldsymbol{v}}}^{d}) \end{split}$$ (29) 联立式(26)和式(29), 可得
$$ \begin{array}{*{20}{l}} \dot{N}_{2} \leq -\frac{2}{\eta} \Vert P_{\boldsymbol{\ell}}(-\eta A^\mathrm{T}{{\boldsymbol{\epsilon}}+A^{\dagger}\dot{{\boldsymbol{v}}}^{d})-A^{\dagger}\dot{{\boldsymbol{v}}}^{d}}\Vert^{2} \leq 0 \end{array} $$ (30) 根据拉萨尔不变集原理, 当$ \dot{N}_{2} $ = 0时, 可以得到
$$ \begin{array}{*{20}{l}} \begin{split} P_{\boldsymbol{\ell}}(-\eta A^\mathrm{T}{{\boldsymbol{\epsilon}}+A^{\dagger}\dot{{\boldsymbol{v}}}^{d})=A^{\dagger}\dot{{\boldsymbol{v}}}^{d}} \end{split} \end{array} $$ (31) 当$ A^{\mathrm{T}}{\boldsymbol{\beta}}\in \boldsymbol{\ell} $时, 可得$ -\eta A^\mathrm{T}{\boldsymbol{\epsilon}}={\boldsymbol{0}}$, 因此$ {{\boldsymbol{\epsilon}}}={\boldsymbol{0}}$. 综上, 本文所提方案(式(14))的控制误差全局收敛至$ {{\boldsymbol{0}}}$.
□ 3. 仿真与实验
为验证本文所提方案对机器人雅可比矩阵估计的准确性及末端执行器位姿控制的有效性, 本节基于德国7自由度Franka Emika Panda 冗余机器人分别给出末端执行器姿态保持和姿态跟踪的两组仿真结果.
值得指出的是, 仿真过程中的相关参数设置如下: 冗余机器人的初始关节角度被设置为$[0,-\pi/4, 0,-3\pi/4,0,\pi/2,3\pi/4]^\mathrm{T}$ rad; 收敛系数$ \mu=500 $; 设计参数$ \eta=500 $; 控制学习的收敛参数$\sigma=1\, 000$; 估计雅可比矩阵初始值为$ [-1, 1] $区间内的随机数; 设计参数$ {\boldsymbol{\beta}=[0, 0, 0, 0, 0, 0, 0, 0]^\mathrm{T}}$. 值得注意的是, 很难通过分析或者直观的设计步骤选择这些参数, 本文则是通过多次仿真确定它们的值.
图1是一组机器人末端执行器展示姿态保持功能的仿真结果. 具体而言, 给定机器人末端执行器跟踪四叶草任务. 图1(a)描绘了机器人执行跟踪任务的整体运动轨迹, 绿色曲线为任务过程中末端执行器的姿态变化, 红色曲线为末端执行器跟踪预定的四叶草图案. 图1(b)和图1(c)分别展示了机器人在执行任务过程中关节角度、关节角速度的变化曲线. 可以看出曲线变化平滑, 机器人关节角度、关节角速度均在合理的范围内. 图1(d)展示了机器人在执行任务过程中估计的雅可比矩阵与实际的雅可比矩阵误差. 随着任务时间的持续, 其误差迅速小于$ 10^{-3} $, 证明本文所提方案(14)能够有效地估计雅可比矩阵, 从而实现机器人的准确控制; 图1(e)给出了机器人任务过程中的四元数误差变化趋势, 可以看出四元数误差小于$ 10^{-5} $, 证明机器人实现了末端执行器姿态保持的功能; 图1(f)给出了机器人在执行任务过程中末端执行器的空间位置误差, 沿三个坐标轴的误差均小于$ 10^{-4} $ m.
图 1 采用所提方案(14)实现冗余机器人末端执行器位置跟踪与姿态保持的仿真结果Fig. 1 Simulation results of the redundant manipulator using the proposed scheme (14) to achieve position tracking and orientation maintenance图2为一组机器人末端执行器展示姿态跟踪功能的仿真结果. 具体而言, 机器人末端执行器沿着圆锥体表面进行打磨、焊接等任务. 值得指出的是, 在任务执行期间, 机器人末端执行器的方向需要一直保持在圆锥体表面的法线方向. 不同于图1中姿态保持的仿真结果, 图2(a)展示了机器人执行跟踪任务的整体运动轨迹, 绿色曲线描绘了机器人末端执行器位姿, 可见末端执行器位姿是不断变化的. 图2(b)和图2(c)中关节角度和关节角速度的变化曲线表明了机器人执行过程中关节运行平滑. 图2(d)、图2(e)和图2(f)展示了机器人任务过程中雅可比矩阵误差、冗余机器人末端执行器的姿态与位置误差变化趋势, 证明了由所提方案(14)驱动的机器人在未知机器人结构信息情况下成功地执行了位姿跟踪任务. 此外, 为更好地说明本方案对末端执行器位姿的控制过程, 考虑图2的应用场景, 图3给出了基于CoppeliaSim仿真平台所提方案(14)与文献[22]方案的对比仿真实验. CoppeliaSim是一款功能强大的动力学仿真软件, 它可以根据物体的物理特性计算运动、旋转和碰撞, 因此许多研究应用CoppeliaSim以快速验证算法. 对比仿真实验中, 由所提方案(14)驱动的机器人在执行任务的过程中, 其末端执行器的方向始终保持在圆锥体表面的法线方向. 相比之下, 由文献[22]方案驱动的机器人在执行过程中因无法控制其末端执行器方向而导致任务失败. 由此可见, 机器人末端执行器的姿态跟踪在一些要求姿态不断变化的应用场景中具备重要的作用.
图 2 采用所提方案(14)实现冗余机器人位置与姿态跟踪的仿真结果Fig. 2 Simulation results of the redundant manipulator using the proposed scheme (14) to achieve position and orientation tracking
图 3 基于CoppeliaSim平台冗余机器人实现位置与姿态跟踪的对比结果Fig. 3 Comparison results of the redundant manipulator achieving position and orientation tracking based on CoppeliaSim platform此外, 表2给出了所提方案(14)与其他5种方案的对比. 具体而言, 所提方案(14)与对比方案使用同样的机器人执行相同的任务, 即图2中的应用场景, 同时选取方案层级、末端执行器控制模式、结构信息是否已知以及跟踪误差等多个指标对方案进行了评价. 值得指出的是, 每个方案参数的选择都基于其对应的文献. 从表2可看出, 文献[13]为基于数据驱动的速度层控制方案, 该方案能够实现末端执行器姿态保持的功能; 文献[22]和文献[23]为加速度层的控制方案, 需要已知机器人的模型信息; 文献[24]和文献[25]为速度层的控制方案, 同样需要已知机器人的模型信息. 此外, 从位置跟踪误差方面来看, 每个方案的跟踪误差均维持在$ 10^{-3} $ m数量级. 在姿态跟踪误差方面, 多数方案不支持姿态跟踪, 文献[13]和文献[25]的方案的跟踪精度高于所提方案(14). 然而, 文献[13]方案仅具备简单的姿态保持功能, 文献[25]方案则无法应用于结构信息不准确或未知的机器人. 综上, 本文所提方案(14)能够有效地控制结构信息未知的机器人进行轨迹跟踪任务, 在对末端执行器姿态有约束的情况下具有一定优势.
表 2 冗余机器人不同轨迹跟踪控制方案对比Table 2 Comparison of different trajectory tracking control schemes for redundant manipulators4. 结束语
本文针对模型未知的冗余机器人提出了一种基于数据驱动的末端执行器位姿控制方案. 该方案从机器人运动学出发, 用四元数表示末端执行器的姿态, 进而实现了对末端执行器位姿控制. 通过引入数据驱动技术, 本文所提方案能够有效地控制模型未知的冗余机器人, 且估计的雅可比矩阵能够收敛至真实值. 此外, 设计了一种神经动力学方案用于本方案的求解. 理论分析证明了所提求解器的收敛性, 方案对比和仿真实验体现了所提方案对机器人末端执行器控制的优越性和有效性.
-
图 3 不同参数的MGR-WT对图像的分割结果对比
Fig. 3 Segmentation results comparison using MRG-WT with different values of r
图 4 RWT与MGR-WT对陶瓷材料晶粒的分割结果对比
Fig. 4 Comparison of the segmentation results of ceramic grains between RWT and MGR-WT
图 9 基于LRCF与分水岭变换的图像分割
Fig. 9 Image segmentation using the combination of LRCF and watershed transform
图 11 第一组分割结果对比(未镀金图像)
Fig. 11 Comparison of the first group of segmentation results (unplated image)
图 12 第二组分割结果对比(镀金图像)
Fig. 12 Comparison of the first group of segmentation results (gilded image)
表 1 不同方法对陶瓷晶粒分割的性能指标对比(第1组实验, 未镀金的图像)
Table 1 Performance comparison of different approaches for ceramic grain segmentation (the first group of experiments for unplated image)
Methods CV↑ VI↓ GCE↓ BDE↓ Liu＇s-MGR[38] 0.2889 3.4270 0.4742 7.3230 Random Walker[39] 0.3556 2.9003 0.1407 13.2147 SLIC[14] 0.3547 3.0524 0.4396 10.1678 LSC[15] 0.3455 2.8820 0.3563 7.5911 Banerjee＇s[30] 0.5959 2.1992 0.2031 3.9182 SE-MGR-WT[32] 0.4680 2.3887 0.1364 5.0346 SE-AMR-WT[40] 0.8287 1.1280 0.1122 1.6261 RCF-MGR-WT[23] 0.6636 1.4952 0.0955 3.5651 LRCF-RWT 0.8697 0.8710 0.0763 1.6262 
下载: 导出CSV
表 2 不同方法对陶瓷晶粒分割的性能指标对比(第2组实验, 镀金的图像)
Table 2 Performance comparison of different approaches for ceramic grain segmentation (the second group of experiments for gilded image)
Methods CV↑ VI↓ GCE↓ BDE↓ Liu＇s-MGR[38] 0.2622 3.8053 0.3565 6.9440 Random Walker[39] 0.3823 2.9517 0.2202 16.4378 SLIC[14] 0.3279 3.0962 0.4070 11.3350 LSC[15] 0.3347 2.8418 0.3265 8.0651 Banerjee＇s[30] 0.7035 1.7175 0.1052 2.7484 SE-MGR-WT[32] 0.7979 1.2031 0.1033 2.0565 SE-AMR-WT[40] 0.8757 0.9909 0.1110 1.2623 RCF-MGR-WT[23] 0.5771 1.7691 0.0895 4.8813 LRCF-RWT 0.9217 0.6699 0.0628 1.0201
下载: 导出CSV
表 3 人工测量晶粒尺寸结果(像素)
Table 3 Grain sizes using manual method (pixels)
测量者 1 测量者 2 测量者 3 测量者 4 测量者 5 1 94.55 89.17 93.39 94.22 88.51 2 90.92 100.33 105.38 91.48 99.91 3 107.50 100.91 102.09 96.49 89.91 4 101.61 89.91 92.08 94.42 93.38 5 108.31 103.88 95.16 102.45 93.52 6 112.51 108.21 112.34 109.70 107.84 7 101.85 104.13 102.80 94.40 89.73
下载: 导出CSV
表 4 不同方法对陶瓷晶粒尺寸的计算结果对比(像素)
Table 4 Comparison of ceramic grain sizes using different approaches (pixels)
人工测量 Ground Truth Liu＇s-MGR[38] RW[39] SLI[14] LSC[15] [30] SE-MGR-WT[32] SE-AMR-WT[40] RCF-MGR-WT[23] LRCF-RWT 1 92.26 97.80 88.00 195.16 74.33 63.95 92.58 48.88 83.73 63.07 98.56 2 97.24 98.00 85.60 161.54 74.48 63.66 86.59 55.09 94.34 75.08 99.15 3 99.83 92.33 82.81 175.15 76.66 62.39 105.29 50.92 90.52 63.08 92.47 4 93.29 93.34 65.97 206.96 75.72 62.73 86.45 53.17 87.70 65.21 92.48 5 100.50 96.09 74.38 192.80 75.99 68.04 102.02 67.25 93.87 59.95 96.76 6 110.08 98.93 69.83 177.56 76.48 70.01 104.08 76.38 96.00 59.31 100.65 7 99.68 96.61 78.18 183.03 75.50 71.71 114.28 85.29 93.98 53.59 97.67
下载: 导出CSV
表 5 不同方法计算陶瓷晶粒尺寸结果的误差(像素)
Table 5 Error comparison of different approaches on ceramic grain size computation (pixels)
Liu＇s-MGR[38] RW[39] SLIC[14] LSC[15] [30] SE-MGR-WT[32] SE-AMR-WT[40] RCF-MGR-WT[23] LRCF-RWT 1 −9.80 +97.36 −23.47 −33.85 −5.22 −48.92 −14.07 −34.73 −0.76 2 −12.40 +63.54 −23.52 −34.34 −11.41 −42.91 −3.66 −22.92 +1.15 3 −9.52 +82.82 −15.67 −29.94 −12.96 −41.41 −1.81 −29.25 −0.14 4 −27.37 +113.62 −17.62 −30.61 −6.89 −40.17 −5.64 −28.13 −0.86 5 −21.71 +96.71 −20.1 −28.05 +6.07 −28.84 −2.22 −36.14 −0.67 6 −29.10 +18.63 −22.45 −28.92 +5.15 −19.55 −2.93 −39.62 +1.72 7 −18.43 +86.42 −21.11 −24.90 +17.67 −11.32 −2.63 −43.02 −1.06
下载: 导出CSV
-
[1] Viola G, Chong K B, Eriksson M, Shen Z, Zeng J, Yin Q, et al. Effect of grain size on domain structures, dielectric and thermal depoling of Nd-substituted bismuth titanate ceramics. Applied Physics Letters, 2013, 103(18): 182903. [2] Yuan Q, Yao F, Wang Y, Ma R, Wang H. Relaxor ferroelectric 0.9 BaTiO3–0.1 Bi (Zn0.5Zr0.5) O3 ceramic capacitors with high energy density and temperature stable energy storage properties. Journal of Materials Chemistry C, 2017, 5(37): 9552-9558. [3] Yang Z, Gao F, Du H, Jin L, Yan L, Hu Q, et al. Grain size engineered lead-free ceramics with both large energy storage density and ultrahigh mechanical properties. Nano Energy, 2019, 58: 768-777. [4] 邹耀斌, 雷帮军, 臧兆祥, 王俊英, 胡泽海, 董方敏. 归一化互信息量最大化导向的自动阈值选择方法. 自动化学报, 2019, 45(7): 1373-1385.Zou Yao-Bin, Lei Bang-Jun, Zang Zhao-Xiang, Wang Jun-Ying, Hu Ze-Hai, Dong Fang-Min. Automatic threshold selection guided by maximizing normalized mutual information. Acta Electronica Sinica, 2019, 45(7): 1373-1385. [5] Yang M, Liu S, Chen K, Zhang H, Zhao E, Zhao T. A hierarchical clustering approach to fuzzy semantic representation of rare words in neural machine translation. IEEE Transactions on Fuzzy Systems, 2020, 28(5): 992-1002. [6] Lei T, Jia X, Zhang Y, He L, Meng H, Nandi A K. Significantly fast and robust fuzzy c-means clustering algorithms based on morphological reconstruction and membership filtering. IEEE Transactions on Fuzzy Systems, 2018, 26(5): 3027-3041. [7] 陶新民, 王若彤, 常瑞, 李晨曦, 刘艳超. 基于低密度分割密度敏感距离的谱聚类算法. 自动化学报, 2020, 46(7): 1479-1495.Tao Xin-Min, Wang Ruo-Tong, Chang Rui, Li Chen-Xi, Liu Yan-Chao. Low density separation density sensitive distance-based spectral clustering algorithm. Acta Electronica Sinica, 2020, 46(7): 1479-1495. [8] 雷涛, 张肖, 加小红, 刘侍刚, 张艳宁. 基于模糊聚类的图像分割研究进展. 电子学报, 2019, 47(08): 1776-1791. doi: 10.3969/j.issn.0372-2112.2019.08.023Lei Tao, Zhang Xiao, Jia Xiao-Hong, Liu Shi-Gang, Zhang Yan-Ning. Research progress on image segmentation based on fuzzy clustering. Acta Electronica Sinica, 2019, 47(8): 1776-1791. doi: 10.3969/j.issn.0372-2112.2019.08.023 [9] 马超, 刘亚淑, 骆功宁, 王宽全. 基于级联随机森林与活动轮廓的3D MR图像分割. 自动化学报, 2019, 45(5): 1004-1014.Ma Chao, Liu Ya-Shu, Luo Gong-Ning, Wang Kuan-Quan. Combining concatenated random forests and active contour for the 3D MR images segmentation. Acta Electronica Sinica, 2019, 45(5): 1004-1014. [10] Zhang L, Wu K T, Li P. An improved adaptive level set method for image segmentation. International Journal of Pattern Recognition & Artificial Intelligence, 2017, 32(5): 1854013. [11] Liu X, Guo S, Yang B, Ma S, Fu Y. Automatic organ segmentation for CT scans based on super-pixel and convolutional neural networks. Journal of Digital Imaging, 2018, 31(6): 748-760. [12] Zhang Z, Pan X, Cheng L, Zhan S, Zhou H, Chen R, et al. SAR image segmentation using hierarchical region merging with orientated edge strength weighted Kuiper’s distance. IEEE Access, 2019, 7: 84479-84496. doi: 10.1109/ACCESS.2019.2910735 [13] Ghosh P, Mali K, Das S K. Use of spectral clustering combined with normalized cuts (N-Cuts) in an iterative k-means clustering framework (NKSC) for superpixel segmentation with contour adherence. Pattern Recognition and Image Analysis, 2018, 28(3): 400-409. [14] Achanta R, Shaji A, Smith K, Lucchi A, Fua P, Süsstrunk, S. SLIC superpixels compared to state-of-the-art superpixel methods. IEEE Transactions on Pattern Analysis and Machine Intelligence, 2012, 34(11): 2274-2282. [15] Li Z, Chen J. Superpixel segmentation using linear spectral clustering. IEEE Transactions on Image Processing, 2017, 26(7): 3317-3330. doi: 10.1109/TIP.2017.2651389 [16] Liu M Y, Tuzel O, Ramalingam S, Chellappa R. Entropy rate superpixel segmentation. In: Proceedings of IEEE Computer Vision and Pattern Recognition, CO, USA: IEEE, 2011. 2097-2104 [17] Zhao Z Q, Zheng P, Xu S, Wu X. Object detection with deep learning: A review. IEEE transactions on neural networks and learning systems, 2019, 30(11): 3212-3232. doi: 10.1109/TNNLS.2018.2876865 [18] Wang K, Liew J H, Zou Y, Zhou D, Feng J. Panet: Few-shot image semantic segmentation with prototype alignment. In: Proceedings of International Conference on Computer Vision. Seoul, Korea: IEEE, 2019. 9197−9206 [19] Long J, Shelhamer E, Darrell T. Fully convolutional networks for semantic segmentation. IEEE Transactions on Pattern Analysis and Machine Intelligence, 2014, 39(4): 640-651. [20] He K, Zhang X, Ren S, Sun J. Spatial pyramid pooling in deep convolutional networks for visual recognition. IEEE Transactions on Pattern Analysis and Machine Intelligence, 2015, 37(9): 1904-1916. doi: 10.1109/TPAMI.2015.2389824 [21] Zhao H, Shi J, Qi X, Wang X, Jia J. Pyramid scene parsing network. In: Proceedings of IEEE Computer Vision and Pattern Recognition, Honolulu, USA: IEEE, 2017. 6230−6239 [22] Chen L C, Papandreou G, Kokkinos I, Murphy K, Yuille A L. DeepLab: Semantic image segmentation with deep convolutional nets, atrous convolution, and fully connected CRFs. IEEE Transactions on Pattern Analysis and Machine Intelligence, 2018, 40(4): 834-848.. doi: 10.1109/TPAMI.2017.2699184 [23] Liu Y, Cheng M M, Hu X, Wang K, Bai X. Richer convolutional features for edge detection. IEEE Transactions on Pattern Analysis and Machine Intelligence, 2019, 41(8): 1939-1946. doi: 10.1109/TPAMI.2018.2878849 [24] Gao S, Cheng M M, Zhao K, Zhang X Y, Yang M H, Philip H S T. Res2net: A new multi-scale backbone architecture. IEEE transactions on pattern analysis and machine intelligence, 2019. doi: 10.1109/TPAMI.2019.2938758. [25] Li H, Xiong P, Fan H, Sun J. DFANet: Deep feature aggregation for real-time semantic segmentation. In: Proceedings of IEEE Computer Vision and Pattern Recognition, LA, USA: IEEE, 2019. 9522−9531 [26] Ding H, Jiang X, Shuai B, Liu A Q, Wang G. Semantic segmentation with context encoding and multi-path decoding. IEEE Transactions on Image Processing, 2020, 29: 3520-3533. doi: 10.1109/TIP.2019.2962685 [27] Andrew M. A quantified study of segmentation techniques on synthetic geological XRM and FIB-SEM images. Computational Geosciences, 2018, 22(6): 1503-1512. doi: 10.1007/s10596-018-9768-y [28] 薛维华, 王浩, 刘国权, 冯敏楠, 齐琳, 陈展鹏. 基于图像序列间相似性的晶粒组织图像分割方法. 中国体视学与图像分析, 2015, 20(04): 400-407.Xue Wei-Hua, Wang Hao, Liu Guo-Quan, Feng Min-Nan, Qi Lin, Chen Zhan-Peng. Segmentation method of grain microstructure based on similarity among image sequence. Chinese Journal of Stereology and Image Analysis, 2015, 20(04): 400-407. [29] Jiang F, Gu Q, Hao H, Li N, Wang B, Hu X. A method for automatic grain segmentation of multi-angle cross-polarized microscopic images of sandstone. Computers & Geosciences, 2018, 115: 143-153. [30] Banerjee S, Chakraborti P C, Saha S K. An automated methodology for grain segmentation and grain size measurement from optical micrographs. Measurement, 2019, 140: 142-150. doi: 10.1016/j.measurement.2019.03.046 [31] Vincent L, Soille P. Watersheds in digital spaces: An efficient algorithm based on immersion simulations. IEEE Transactions on Pattern Analysis and Machine Intelligence, 1991, 13(6): 583-598. doi: 10.1109/34.87344 [32] Dollár P, Zitnick C L. Fast edge detection using structured forests. IEEE Transactions on Pattern Analysis and Machine Intelligence, 2015, 37(8): 1558-1570. doi: 10.1109/TPAMI.2014.2377715 [33] Jobson D J, Rahman Z, Woodell G A. Properties and performance of a center/surround retinex. IEEE transactions on image processing, 1997, 6(3): 451-462. doi: 10.1109/83.557356 [34] Jobson D J, Rahman Z, Woodell G A. A multiscale retinex for bridging the gap between color images and the human observation of scenes. IEEE Transactions on Image Processing, 1997, 6(7): 965-976. doi: 10.1109/83.597272 [35] Jackway P T. Gradient watersheds in morphological scale-space. IEEE Transactions on Image Processing, 1996, 5(6): 913-921. doi: 10.1109/83.503908 [36] Yu Z, Feng C, Liu M Y, Ramalingam S. Casenet: Deep category-aware semantic edge detection. In: Proceedings of IEEE Computer Vision and Pattern Recognition, Honolulu, USA: IEEE, 2017. 5964−5973 [37] Liu Y, Lew M S. Learning relaxed deep supervision for better edge detection. In: Proceedings of IEEE Computer Vision and Pattern Recognition, NV, USA: IEEE, 2016. 231−240 [38] Liu G, Zhang Y, Wang A. Incorporating adaptive local onformation onto fuzzy clustering for image segmentation. IEEE Transactions on Image Processing, 2015, 24(11): 3990-4000. doi: 10.1109/TIP.2015.2456505 [39] Grady L. Random walks for image segmentation. IEEE Transactions on Pattern Analysis & Machine Intelligence, 2006, 28(11): 1768-1783. [40] Lei T, Jia X, Liu T, Liu S, Meng H, Nandi A K. Adaptive morphological reconstruction for seeded image segmentation. IEEE Transactions on Image Processing, 2019, 28(11): 5510-5523. doi: 10.1109/TIP.2019.2920514 [41] Arbelaez P, Maire M, Fowlkes C, Malik J. Contour detection and hierarchical image segmentation. IEEE Transactions on Pattern Analysis and Machine Intelligence, 2011, 33(5): 898-916. doi: 10.1109/TPAMI.2010.161 [42] Meila M. Comparing clusterings: An axiomatic view. In: Proceedings of the International Conference on Machine Learning, Bonn, Germany: IMLS, 2005. 577−584 [43] Martin D, Fowlkes C, Tal D, Malik J. A database of human segmented natural images and its application to evaluating segmentation algorithms and measuring ecological statistics. In: Proceedings of International Conference on Computer Vision. BC, Canada: IEEE, 2001. 416−423 [44] Freixenet J, Muñoz X, Raba D, Martí J, Cufí X. Yet another survey on image segmentation: Region and boundary information integration. In: Proceedings of the European Conference on Computer Vision, Copenhagen, Denmark: IEEE, 2002. 408−422 期刊类型引用(2)
1. 马靖雯，张廷丰，陆明超. 基于数据驱动的跨声速风洞控制方法研究. 辽宁工业大学学报(自然科学版). 2024(05): 298-302+309 . 
百度学术2. 曾旭翔，孔颖. 求解时变二次规划的自适应参数归零神经网络. 浙江科技大学学报. 2024(05): 384-393 . 百度学术其他类型引用(3)
-
下载:
计量
- 文章访问数: 1998
- HTML全文浏览量: 752
- PDF下载量: 334
- 被引次数: 5
