CN114049310A - 一种磁悬浮轴承转子与其保护轴承相对位置图像分析方法 - Google Patents

一种磁悬浮轴承转子与其保护轴承相对位置图像分析方法 Download PDF

Info

Publication number
CN114049310A
CN114049310A CN202111248602.5A CN202111248602A CN114049310A CN 114049310 A CN114049310 A CN 114049310A CN 202111248602 A CN202111248602 A CN 202111248602A CN 114049310 A CN114049310 A CN 114049310A
Authority
CN
China
Prior art keywords
rotor
bearing
matrix
circle
image
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
CN202111248602.5A
Other languages
English (en)
Other versions
CN114049310B (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.)
Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical University
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 Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN202111248602.5A priority Critical patent/CN114049310B/zh
Publication of CN114049310A publication Critical patent/CN114049310A/zh
Application granted granted Critical
Publication of CN114049310B publication Critical patent/CN114049310B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/0002—Inspection of images, e.g. flaw detection
    • G06T7/0004—Industrial image inspection
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/60—Analysis of geometric attributes
    • G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30108—Industrial image inspection
    • G06T2207/30164—Workpiece; Machine component

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • Quality & Reliability (AREA)
  • Length Measuring Devices By Optical Means (AREA)
  • Magnetic Bearings And Hydrostatic Bearings (AREA)

Abstract

本发明公开了一种磁悬浮轴承转子与其保护轴承相对位置图像分析方法，首先基于转子和保护轴承的具体尺寸以及摄像机的分辨率建立数学图像模型和相应的数据库；然后使用相机拍摄转子和其保护轴承的图像，并对图像进行处理，使用图像增强的方法将转子和保护轴承区域灰度值置成0，两者中间间隔区域置成1；最后处理图像对应矩阵，统计出在图像的右侧和上侧区域值为1的个数，根据1的个数，来通过对数据库数据的检索，再由单个像素与实际尺寸的关系，确定转子和其保护轴承的相对位置。本发明方法简单，计算便捷，有利于实现转子的高精度位置测量。

Description

一种磁悬浮轴承转子与其保护轴承相对位置图像分析方法
技术领域
本发明属于图像处理技术领域，具体涉及一种磁悬浮轴承转子与其保护轴承位置分析方法。
背景技术
磁悬浮轴承稳定运转时无机械接触，转子转速上限高，且具备机械磨损小、能耗低、噪声小、寿命长、无需润滑、无油污染等优点，特别适用于高速、真空、超净等特殊环境中。而且磁悬浮轴承的刚度、阻尼可调节，因此在旋转机械领域的应用前景良好。磁悬浮转子的位置是反馈控制中的反馈量，其精度直接影响到磁悬浮系统的稳定性和稳态精度，所以位置测量技术在磁悬浮控制系统中至关重要。当前磁悬浮控制轴承位置测量手段主要基于电涡流位移传感器、光电传感器、磁敏传感器、电感式传感器和电容式传感器等，但这些传感器的精度都存在局限性。
发明内容
为了克服现有技术的不足，本发明提供了一种磁悬浮轴承转子与其保护轴承相对位置图像分析方法，首先基于转子和保护轴承的具体尺寸以及摄像机的分辨率建立数学图像模型和相应的数据库；然后使用相机拍摄转子和其保护轴承的图像，并对图像进行处理，使用图像增强的方法将转子和保护轴承区域灰度值置成0，两者中间间隔区域置成1；最后处理图像对应矩阵，统计出在图像的右侧和上侧区域值为1的个数，根据1的个数，来通过对数据库数据的检索，再由单个像素与实际尺寸的关系，确定转子和其保护轴承的相对位置。本发明方法简单，计算便捷，有利于实现转子的高精度位置测量。
本发明解决其技术问题所采用的技术方案包括如下步骤：
步骤1：基于磁悬浮轴承转子与其保护轴承的尺寸以及摄像机的分辨率建立数学图像模型和相应的数据库；
步骤1-1：根据相机所拍摄照片的分辨率m×n定义矩阵E，矩阵E的大小为m×n；矩阵E中每个元素对应一个像素点；
步骤1-2：设定磁悬浮轴承转子的半径为r1个像素，保护轴承的半径为r2个像素；
步骤1-3：以矩阵E的中心O点建立x-y直角坐标系，每一个像素的长度为1；
步骤1-4：以O点为中心、r2为半径建立一个轴承圆，将轴承圆区域以外的矩阵E元素值全部置成0；在轴承圆内随机选取一点，以该点为圆心、r1为半径建立一个转子圆，将转子圆区域内的矩阵E元素值全部置成0；矩阵E中剩余元素值置成1；转子圆的圆心在以O为圆心、r2-r1为半径的圆D内移动；
步骤1-5：以过矩阵E的中心O点的水平线为界，将矩阵E分为上半侧和下半侧；以过矩阵E的中心O点的竖直线为界，将矩阵E分为左半侧和右半侧；
步骤1-6：将转子圆的圆心遍历圆D中每一个点；转子圆圆心的每一个横坐标x1都对应一个矩阵E右半侧值为1的元素总数，由此对应关系建立一个横坐标数组；转子圆圆心的每一个纵坐标y1都对应一个矩阵E上半侧值为1的元素总数，由此对应关系建立一个纵坐标数组；横坐标数组和纵坐标数组构成数据库；
步骤2：使用相机拍摄磁悬浮轴承转子与其保护轴承的图像建立与步骤1中矩阵对应的图像矩阵；
步骤2-1：使用相机拍摄磁悬浮轴承转子与其保护轴承的图像，使图像中心为保护轴承的中心点；
步骤2-2：相机拍摄的图像的实际尺寸l毫米×k毫米，则图像中每一个像素值对应的实际长度c＝l/m；
步骤2-3：将图像中保护轴承以外区域和磁悬浮轴承转子区域的像素灰度值全部置位0，磁悬浮轴承转子和保护轴承之间的区域像素灰度值全部置位1；生成一个大小m×n的矩阵F；
步骤3：确定磁悬浮轴承转子实际位置；
步骤3-1：分别统计矩阵F中右半侧和上半侧值为1的元素个数；
步骤3-2：使用矩阵F中右半侧和上半侧值为1的元素个数在步骤1建立的数据库中进行查询，得到磁悬浮轴承转子的圆心像素坐标值(x，y)；
步骤3-3：计算磁悬浮轴承转子的圆心的实际坐标位置X＝x×c，Y＝y×c，得到磁悬浮轴承转子相对于其保护轴承的相对位置。
优选地，所述磁悬浮轴承转子上涂抹高亮度涂料，方便后期图像处理。
优选地，所述统计矩阵F中右半侧和上半侧值为1的元素个数采用软件自动统计。
本发明的有益效果如下：
本发明方法相比传统采用传感器测量的方法，使用图像处理，实施简单，计算便捷，有利于实现转子的高精度位置测量。
附图说明
图1为本发明实施例经过步骤2处理后相机拍摄图像。
具体实施方式
下面结合附图和实施例对本发明进一步说明。
一种磁悬浮轴承转子与其保护轴承相对位置图像分析方法，包括如下步骤：
步骤1：基于磁悬浮轴承转子与其保护轴承的尺寸以及摄像机的分辨率建立数学图像模型和相应的数据库；
步骤1-1：根据相机所拍摄照片的分辨率m×n定义矩阵E，矩阵E的大小为m×n；矩阵E中每个元素对应一个像素点；
步骤1-2：设定磁悬浮轴承转子的半径为r1个像素，保护轴承的半径为r2个像素；
步骤1-3：以矩阵E的中心O点建立x-y直角坐标系，每一个像素的长度为1；
步骤1-4：以O点为中心、r2为半径建立一个轴承圆，将轴承圆区域以外的矩阵E元素值全部置成0；在轴承圆内随机选取一点，以该点为圆心、r1为半径建立一个转子圆，将转子圆区域内的矩阵E元素值全部置成0；矩阵E中剩余元素值置成1；转子圆的圆心在以O为圆心、r2-r1为半径的圆D内移动；
步骤1-5：以过矩阵E的中心O点的水平线为界，将矩阵E分为上半侧和下半侧；以过矩阵E的中心O点的竖直线为界，将矩阵E分为左半侧和右半侧；
步骤1-6：将转子圆的圆心遍历圆D中每一个点；转子圆圆心的每一个横坐标x1都对应一个矩阵E右半侧值为1的元素总数，由此对应关系建立一个横坐标数组；转子圆圆心的每一个纵坐标y1都对应一个矩阵E上半侧值为1的元素总数，由此对应关系建立一个纵坐标数组；横坐标数组和纵坐标数组构成数据库；
步骤2：使用相机拍摄磁悬浮轴承转子与其保护轴承的图像建立与步骤1中矩阵对应的图像矩阵；
步骤2-1：使用相机拍摄磁悬浮轴承转子与其保护轴承的图像，使图像中心为保护轴承的中心点；
步骤2-2：相机拍摄的图像的实际尺寸l毫米×k毫米，则图像中每一个像素值对应的实际长度c＝l/m；
步骤2-3：使用图像增强的方法将图像中保护轴承以外区域和磁悬浮轴承转子区域的像素灰度值全部置位0，磁悬浮轴承转子和保护轴承之间的区域像素灰度值全部置位1；生成一个大小m×n的矩阵F；
步骤3：确定磁悬浮轴承转子实际位置；
步骤3-1：使用软件分别统计矩阵F中右半侧和上半侧值为1的元素个数；
步骤3-2：使用矩阵F中右半侧和上半侧值为1的元素个数在步骤1建立的数据库中进行查询，得到磁悬浮轴承转子的圆心像素坐标值(x，y)；
步骤3-3：计算磁悬浮轴承转子的圆心的实际坐标位置X＝x×c，Y＝y×c，得到磁悬浮轴承转子相对于其保护轴承的相对位置。
具体实施例：
假设m＝1024，n＝1024，l＝40cm，k＝40cm，计算得到c＝l/m＝0.39mm。基于m和n的数值，生成1024×1024矩阵，然后按照步骤1生成该范例模型下的数据库。其中设定r2长度对应100像素，r1长度对应80像素。
用相机拍摄实际图像，经过步骤2方法处理后得到附图1；然后记录其右半侧和上半侧元素值为1的个数，得到上半侧值为1的个数为4533，右半侧值为1的个数为4096，然后在检索数据库，可以找到在数据库中，上半侧数据中为4533的数据对应纵坐标为-7，右半侧数据中为4096的数据对应横坐标为10,。得到输出转子圆心的坐标(10，-7)，再经过计算从而得到具体的坐标为(3.90，2.73)，单位为毫米。

Claims (3)

1.一种磁悬浮轴承转子与其保护轴承相对位置图像分析方法，其特征在于，包括如下步骤：
步骤1：基于磁悬浮轴承转子与其保护轴承的尺寸以及摄像机的分辨率建立数学图像模型和相应的数据库；
步骤1-1：根据相机所拍摄照片的分辨率m×n定义矩阵E，矩阵E的大小为m×n；矩阵E中每个元素对应一个像素点；
步骤1-2：设定磁悬浮轴承转子的半径为r1个像素，保护轴承的半径为r2个像素；
步骤1-3：以矩阵E的中心O点建立x-y直角坐标系，每一个像素的长度为1；
步骤1-4：以O点为中心、r2为半径建立一个轴承圆，将轴承圆区域以外的矩阵E元素值全部置成0；在轴承圆内随机选取一点，以该点为圆心、r1为半径建立一个转子圆，将转子圆区域内的矩阵E元素值全部置成0；矩阵E中剩余元素值置成1；转子圆的圆心在以O为圆心、r2-r1为半径的圆D内移动；
步骤1-5：以过矩阵E的中心O点的水平线为界，将矩阵E分为上半侧和下半侧；以过矩阵E的中心O点的竖直线为界，将矩阵E分为左半侧和右半侧；
步骤1-6：将转子圆的圆心遍历圆D中每一个点；转子圆圆心的每一个横坐标x1都对应一个矩阵E右半侧值为1的元素总数，由此对应关系建立一个横坐标数组；转子圆圆心的每一个纵坐标y1都对应一个矩阵E上半侧值为1的元素总数，由此对应关系建立一个纵坐标数组；横坐标数组和纵坐标数组构成数据库；
步骤2：使用相机拍摄磁悬浮轴承转子与其保护轴承的图像建立与步骤1中矩阵对应的图像矩阵；
步骤2-1：使用相机拍摄磁悬浮轴承转子与其保护轴承的图像，使图像中心为保护轴承的中心点；
步骤2-2：相机拍摄的图像的实际尺寸l毫米×k毫米，则图像中每一个像素值对应的实际长度c＝l/m；
步骤2-3：将图像中保护轴承以外区域和磁悬浮轴承转子区域的像素灰度值全部置位0，磁悬浮轴承转子和保护轴承之间的区域像素灰度值全部置位1；生成一个大小m×n的矩阵F；
步骤3：确定磁悬浮轴承转子实际位置；
步骤3-1：分别统计矩阵F中右半侧和上半侧值为1的元素个数；
步骤3-2：使用矩阵F中右半侧和上半侧值为1的元素个数在步骤1建立的数据库中进行查询，得到磁悬浮轴承转子的圆心像素坐标值(x，y)；
步骤3-3：计算磁悬浮轴承转子的圆心的实际坐标位置X＝x×c，Y＝y×c，得到磁悬浮轴承转子相对于其保护轴承的相对位置。
2.根据权利要求1所述的一种磁悬浮轴承转子与其保护轴承相对位置图像分析方法，其特征在于，所述磁悬浮轴承转子上涂抹高亮度涂料，方便后期图像处理。
3.根据权利要求1所述的一种磁悬浮轴承转子与其保护轴承相对位置图像分析方法，其特征在于，所述统计矩阵F中右半侧和上半侧值为1的元素个数采用软件自动统计。
CN202111248602.5A 2021-10-26 2021-10-26 一种磁悬浮轴承转子与其保护轴承相对位置图像分析方法 Active CN114049310B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111248602.5A CN114049310B (zh) 2021-10-26 2021-10-26 一种磁悬浮轴承转子与其保护轴承相对位置图像分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111248602.5A CN114049310B (zh) 2021-10-26 2021-10-26 一种磁悬浮轴承转子与其保护轴承相对位置图像分析方法

Publications (2)

Publication Number Publication Date
CN114049310A true CN114049310A (zh) 2022-02-15
CN114049310B CN114049310B (zh) 2024-03-08

Family

ID=80206279

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111248602.5A Active CN114049310B (zh) 2021-10-26 2021-10-26 一种磁悬浮轴承转子与其保护轴承相对位置图像分析方法

Country Status (1)

Country Link
CN (1) CN114049310B (zh)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE19623172C1 (de) * 1996-06-10 1997-10-23 Univ Magdeburg Tech Verfahren zur dreidimensionalen optischen Vermessung von Objektoberflächen
JP2005227929A (ja) * 2004-02-12 2005-08-25 Nippon Telegr & Teleph Corp <Ntt> 対象物の撮影画像の処理方法、画像表示システム、プログラム、および、記録媒体
WO2013004788A1 (de) * 2011-07-07 2013-01-10 Siemens Aktiengesellschaft Verfahren zum betreiben einer rotierbaren vorrichtung, insbesondere medizinisches bildgebendes system, sowie eine solche vorrichtung
RU2013121462A (ru) * 2013-05-07 2014-11-20 Открытое акционерное общество "ТВЭЛ" Устройство для точного позиционирования на цилиндрической поверхности ротора
CN106791331A (zh) * 2017-01-13 2017-05-31 成都微晶景泰科技有限公司 基于透镜阵列成像的图像处理方法、装置及成像系统
CN107037305A (zh) * 2017-03-15 2017-08-11 江苏大学 一种无轴承电机悬浮绕组故障检测方法
CN109730683A (zh) * 2018-12-21 2019-05-10 重庆金山医疗器械有限公司 内窥镜目标物大小计算方法及分析系统
CN110084127A (zh) * 2019-03-29 2019-08-02 南京航空航天大学 一种基于视觉的磁悬浮转子振动测量方法
CN111307073A (zh) * 2020-03-16 2020-06-19 湖南米艾西测控技术有限公司 一种测量旋变定子与转子同轴度偏差的装置

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE19623172C1 (de) * 1996-06-10 1997-10-23 Univ Magdeburg Tech Verfahren zur dreidimensionalen optischen Vermessung von Objektoberflächen
JP2005227929A (ja) * 2004-02-12 2005-08-25 Nippon Telegr & Teleph Corp <Ntt> 対象物の撮影画像の処理方法、画像表示システム、プログラム、および、記録媒体
WO2013004788A1 (de) * 2011-07-07 2013-01-10 Siemens Aktiengesellschaft Verfahren zum betreiben einer rotierbaren vorrichtung, insbesondere medizinisches bildgebendes system, sowie eine solche vorrichtung
RU2013121462A (ru) * 2013-05-07 2014-11-20 Открытое акционерное общество "ТВЭЛ" Устройство для точного позиционирования на цилиндрической поверхности ротора
CN106791331A (zh) * 2017-01-13 2017-05-31 成都微晶景泰科技有限公司 基于透镜阵列成像的图像处理方法、装置及成像系统
CN107037305A (zh) * 2017-03-15 2017-08-11 江苏大学 一种无轴承电机悬浮绕组故障检测方法
CN109730683A (zh) * 2018-12-21 2019-05-10 重庆金山医疗器械有限公司 内窥镜目标物大小计算方法及分析系统
CN110084127A (zh) * 2019-03-29 2019-08-02 南京航空航天大学 一种基于视觉的磁悬浮转子振动测量方法
CN111307073A (zh) * 2020-03-16 2020-06-19 湖南米艾西测控技术有限公司 一种测量旋变定子与转子同轴度偏差的装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
张少军, 苟中魁, 李庆利, 李忠富, 金剑: "利用数字图像处理技术测量直齿圆柱齿轮几何尺寸", 光学精密工程, no. 06, 25 December 2004 (2004-12-25) *
熊文豪;何毅斌;陈宇晨;戴乔森;涂盼盼;: "基于图像识别的齿轮尺寸测量研究", 现代制造工程, no. 03, 18 March 2020 (2020-03-18) *

Also Published As

Publication number Publication date
CN114049310B (zh) 2024-03-08

Similar Documents

Publication Publication Date Title
CN106643545B (zh) 激光位移技术测量钢轨轮廓的校准方法
CN110906875B (zh) 一种孔径测量的视觉处理方法
CN106643965B (zh) 一种利用模板匹配精确识别液位的方法
CN111288967A (zh) 一种基于机器视觉的远距离高精度位移检测方法
CN110763136A (zh) 一种高精度三维工件尺寸检测系统及方法
El-Hakim et al. Multicamera vision-based approach to flexible feature measurement for inspection and reverse engineering
CN108627119B (zh) 一种基于线扫激光的板式家具轮廓尺寸检测方法
Wang et al. Error analysis and improved calibration algorithm for LED chip localization system based on visual feedback
CN107238727A (zh) 基于动态视觉传感器芯片的光电式转速传感器及探测方法
CN109373896A (zh) 一种单方向测量位移的机器视觉系统及其实现方法
CN110458785B (zh) 一种基于图像传感的磁悬浮球悬浮间隙检测方法
CN114049310B (zh) 一种磁悬浮轴承转子与其保护轴承相对位置图像分析方法
Du et al. Roundness detection of end face for shaft workpiece based on Canny-Zernike sub pixel edge detection and improved hough transform
CN105674893B (zh) 基于cmos图像传感器的绝对式光栅尺及其测量方法
CN109000566B (zh) 激光扫描三维成像和ccd二维成像组合测量方法及装置
CN110544276A (zh) 最小二乘法椭圆拟合活塞裙部最大点处尺寸测量方法
CN111521204B (zh) 一种基于绝对位置旋转型编码器的角位移视觉测量方法
CN105606623A (zh) 一种基于线阵相机的钢板表面缺陷检测系统
CN111981983A (zh) 一种复杂形貌物体三维测量的扫描仪精度评价方法
CN110842930A (zh) 基于dlp和相机标定的用于机器人的视觉装置及测量方法
CN112184790B (zh) 基于深度相机的物体尺寸高精度测量方法
CN113012143B (zh) 一种基于二维数字图像相关法的试件质量检测方法
CN210154538U (zh) 一种基于机器视觉的金属结构变形的测量装置
CN209134558U (zh) 一种大视野成像装置
CN114708318B (zh) 一种基于深度相机的未知待测面曲率测定方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant