CN101964293B - 一种金相显微组织图像处理方法 - Google Patents

一种金相显微组织图像处理方法 Download PDF

Info

Publication number
CN101964293B
CN101964293B CN201010262281XA CN201010262281A CN101964293B CN 101964293 B CN101964293 B CN 101964293B CN 201010262281X A CN201010262281X A CN 201010262281XA CN 201010262281 A CN201010262281 A CN 201010262281A CN 101964293 B CN101964293 B CN 101964293B
Authority
CN
China
Prior art keywords
image
pixel
value
pattern
images
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.)
Active
Application number
CN201010262281XA
Other languages
English (en)
Other versions
CN101964293A (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.)
AECC Aviation Power Co Ltd
AVIC Aviation Engine Corp PLC
Original Assignee
Xian Aviation Power Co Ltd
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 Xian Aviation Power Co Ltd filed Critical Xian Aviation Power Co Ltd
Priority to CN201010262281XA priority Critical patent/CN101964293B/zh
Publication of CN101964293A publication Critical patent/CN101964293A/zh
Application granted granted Critical
Publication of CN101964293B publication Critical patent/CN101964293B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

本发明属于理化检测技术,涉及对金相显微组织图像处理方法的改进。处理的步骤如下:通过电子显微镜采集金相图片,规定金相图片的图像模式,对变换到HIS模式的图像进行图像滤波,对滤波后的图像进行Hotelling变换,采用人工神经网络聚类算法处理图像,图像特征提取,计算提取后的所有晶粒元素中的像素点的个数,得到晶粒的面积百分比。本发明大大提高了处理效率和分析精度。

Description

一种金相显微组织图像处理方法
技术领域
本发明属于理化检测技术,涉及对金相显微组织图像处理方法的改进。
背景技术
目前,金相电子显微镜采用传统的数字图像分析方法是依赖试验人员经验进行金相图像分析,具体分析方法是:将同一金相显微图片打印出数十份称重,然后将其中的晶粒部分抠出再次称重,最后通过计算晶粒部分的重量占总重量百分比的方法计算晶粒的面积百分比。这种人工处理方法的缺点是:处理效率和分析精度不高。另外,还有采用数字图像处理软件对金相显微组织图像进行处理的方法,其缺点是:现有的数字图像处理软件在实际应用中无法自动识别图像中的干扰项,无法自动剔除图像干扰,处理精度低。
发明内容
本发明的目的是:提出一种处理效率和分析精度高的金相显微组织图像处理方法。
本发明的技术方案是:一种金相显微组织图像处理方法,处理的步骤如下:
1、通过电子显微镜采集金相图片,图片分辨率不低于800×600;
2、规定金相图片的图像模式:将电子显微采集的金相图片的图像模式由RGB模式转换为HIS模式;
2.1、对RGB三色信号做归一化处理:归一化处理的具体方法是:将图像中任意像素点RGB三色信号色度值的除以该像素点的Bit值,得到图像所有像素点数值在0~1之间的归一值;
2.2、对进行RGB模式到HIS模式的转换,转换公式如下:
I=1/3(R+G+B)............................................................[1]
S=1-3[min{R,G,B}]/(R+G+B)..........................................[2]
H=arccos{1[(R-G)+(R+B)]/2[(R-G)2+(R-B)(G-B)]1/2}...............[3]
上面三式中:I表示亮度,S表示饱和度,H表示色度;
3、对变换到HIS模式的图像进行图像滤波:采用Epanechnikov小波函数均值平移向量的方法,对HIS模式的图像进行滤波,具体方法是:以任意像素点为中心像素点,计算中心像素点周围8个像素点与该中心像素点的距离,然后,将距离小于设定阈值的像素点的平均像素值作为该中心像素点的色度值;
4、对滤波后的图像进行Hotelling变换;
5、采用人工神经网络聚类算法处理图像:具体方法是:设定一个色度阈值,将滤波后图像中的所有像素点分成两类,大于色度阈值的像素点为黑色像素点,小于色度阈值的像素点为白色像素点,将图像变为二值黑白图像;
6、图像特征提取:利用数学形态算法进行图像特征提取,去除比晶粒结构元素小的图像细节;
7、计算提取后的所有晶粒元素中的像素点的个数,然后将晶粒元素的像素点个数除以图像的总像素点个数,即可得到晶粒的面积百分比。
本发明的优点是:大大提高了处理效率和分析精度。与人工处理方法相比,处理效率提高了100倍以上,分析精度比人工处理方法和数字图像处理软件处理方法提高了20%以上。
具体实施方式
下面对本发明做进一步详细说明。一种金相显微组织图像处理方法,处理的步骤如下:
1、通过电子显微镜采集金相图片,图片分辨率不低于800×600;
2、规定金相图片的图像模式:将电子显微采集的金相图片的图像模式由RGB模式转换为HIS模式;
现有的电子显微采集金相图片为RGB模式图片(红、绿、蓝三色组成),由于RGB三个分量的相关性很强,例如当光线发生变化时,RGB三个分量将同时发生变化。因此,检测结果与转换为灰度图像后的检测结果相比变化不大。为了克服RGB颜色空间的不均匀和不直观的缺点,在彩色图像处理中可以采用更加符合颜色视觉特性的颜色空间。HSI颜色空间便是其中的一种,HSI颜色空间是从人的视觉感知角度建立的。H为色度,是描述纯色的属性;S为饱和度,给出一种纯色被白光稀释的程度的度量;I为亮度,是视觉对光刺激的亮度响应。HSI颜色空间在彩色图像中从携带的彩色信息(色度和饱和度)里消去强度分量的影响,因此这种颜色空间对人来说是自然的、直观的。
2.1、对RGB三色信号做归一化处理:归一化处理的具体方法是:将图像中任意像素点RGB三色信号色度值的除以该像素点的Bit值,得到图像所有像素点数值在0~1之间的归一值;
2.2、对进行RGB模式到HIS模式的转换,转换公式如下:
I=1/3(R+G+B)............................................................[1]
S=1-3[min{R,G,B}]/(R+G+B)..........................................[2]
H=arccos{1[(R-G)+(R+B)]/2[(R-G)2+(R-B)(G-B)]1/2}...............[3]
上面三式中:I表示亮度,S表示饱和度,H表示色度;
3、对变换到HIS模式的图像进行图像滤波:采用Epanechnikov小波函数均值平移向量的方法,对HIS模式的图像进行滤波,具体方法是:以任意像素点为中心像素点,计算中心像素点周围8个像素点与该中心像素点的距离,然后,将距离小于设定阈值的像素点的平均像素值作为该中心像素点的色度值;参见《数字图像处理与分析实例教程》(赵书兰主编化学工业出版社2009年版222-234页)
4、对滤波后的图像进行Hotelling变换;
将Hotelling变换应用到金相显微图像各像素数据值组成的数据集合,将颜色进行降维处理即将色彩空间将三维降到二维。应用此变换降维后,晶粒的图形与周围背景的可区分度并没有降低,同时带来了数据运算量的大幅度减少。详细算法参见《数字图像处理》(朱秀昌主编北京邮电大学出版社2002年版51-52页)。
5、采用人工神经网络聚类算法处理图像:具体方法是:设定一个色度阈值,将滤波后图像中的所有像素点分成两类,大于色度阈值的像素点为黑色像素点,小于色度阈值的像素点为白色像素点,将图像变为二值黑白图像;
6、图像特征提取:利用数学形态算法进行图像特征提取,数学形态学算法是应用在模式识别领域的算法,用独立的变换和数学概念描述图像基本特征,去除比晶粒结构元素小的图像细节;详细算法参见《数字图像处理》(朱秀昌主编北京邮电大学出版社2002年版144-145页)。
7、计算提取后的所有晶粒元素中的像素点的个数,然后将晶粒元素的像素点个数除以图像的总像素点个数,即可得到晶粒的面积百分比。
实施例
通过对一组金相显微晶粒图片进行分析得到如下数据:
使用方法           分析用时    分析结果误差;
人工分析方法       2小时       约为20%;
现有软件分析方法   30分钟      约为10%;
本发明方法         10分钟      约为5%。

Claims (1)

1.一种金相显微组织图像处理方法,处理的步骤如下:
1.1、通过电子显微镜采集金相图片,图片分辨率不低于800×600;
1.2、规定金相图片的图像模式:将电子显微采集的金相图片的图像模式由RGB模式转换为HIS模式;
1.2.1、对RGB三色信号做归一化处理:归一化处理的具体方法是:将图像中任意像素点RGB三色信号色度值的除以该像素点的Bit值,得到图像所有像素点数值在0~1之间的归一值;
1.2.2、对进行RGB模式到HIS模式的转换,转换公式如下:
I=1/3(R+G+B)............................................................[1]
S=1-3[min{R,G,B}]/(R+G+B)..........................................[2]
H=arccos{1[(R-G)+(R+B)]/2[(R-G)2+(R-B)(G-B)]1/2}............[3]
上面三式中:I表示亮度,S表示饱和度,H表示色度;
1.3、对变换到HIS模式的图像进行图像滤波:采用Epanechnikov小波函数均值平移向量的方法,对HIS模式的图像进行滤波,具体方法是:以任意像素点为中心像素点,计算中心像素点周围8个像素点与该中心像素点的距离,然后,将距离小于设定阈值的像素点的平均像素值作为该中心像素点的色度值;
1.4、对滤波后的图像进行Hotelling变换;
1.5、采用人工神经网络聚类算法处理图像:具体方法是:设定一个色度阈值,将滤波后图像中的所有像素点分成两类,大于色度阈值的像素点为黑色像素点,小于色度阈值的像素点为白色像素点,将图像变为二值黑白图像;
1.6、图像特征提取:利用数学形态算法进行图像特征提取,去除比晶粒结构元素小的图像细节;
1.7、计算提取后的所有晶粒元素中的像素点的个数,然后将晶粒元素的像素点个数除以图像的总像素点个数,即可得到晶粒的面积百分比。
CN201010262281XA 2010-08-23 2010-08-23 一种金相显微组织图像处理方法 Active CN101964293B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201010262281XA CN101964293B (zh) 2010-08-23 2010-08-23 一种金相显微组织图像处理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201010262281XA CN101964293B (zh) 2010-08-23 2010-08-23 一种金相显微组织图像处理方法

Publications (2)

Publication Number Publication Date
CN101964293A CN101964293A (zh) 2011-02-02
CN101964293B true CN101964293B (zh) 2012-01-18

Family

ID=43517133

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201010262281XA Active CN101964293B (zh) 2010-08-23 2010-08-23 一种金相显微组织图像处理方法

Country Status (1)

Country Link
CN (1) CN101964293B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103824268A (zh) * 2014-02-08 2014-05-28 江西赛维Ldk太阳能高科技有限公司 晶粒图像边缘连接方法及其装置
CN108446706A (zh) * 2018-02-27 2018-08-24 西安交通大学 一种基于颜色主分量提取的磨粒材质自动识别方法

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102331391B (zh) * 2011-06-16 2012-11-14 东方电气集团东方锅炉股份有限公司 一种测算耐热钢中多相颗粒聚集粗化程度的方法
CN102494976B (zh) * 2011-11-18 2014-04-09 江苏大学 一种超细晶粒钢晶粒的自动测量及其形态分类方法
CN102494987A (zh) * 2011-12-13 2012-06-13 天津卓朗科技发展有限公司 球墨铸铁中显微粒子的自动分类评级方法
CN103018248B (zh) * 2012-12-11 2016-03-02 国电锅炉压力容器检验中心 基于轮廓跟踪对珠光体球化的分级方法
CN106353326A (zh) * 2016-08-12 2017-01-25 京东方科技集团股份有限公司 金属层中小丘的检测方法和装置
CN111325706A (zh) * 2018-12-13 2020-06-23 阿里巴巴集团控股有限公司 一种晶界提取、晶界提取模型生成方法及装置
CN111487272B (zh) * 2020-04-21 2023-06-02 中国航发沈阳发动机研究所 一种航空发动机涡轮叶片表面产物层分析方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4488245A (en) * 1982-04-06 1984-12-11 Loge/Interpretation Systems Inc. Method and means for color detection and modification
CN1258058A (zh) * 1998-12-14 2000-06-28 佳能株式会社 图像处理方法和装置,图像处理系统以及存储介质
CN101620060A (zh) * 2009-08-13 2010-01-06 上海交通大学 颗粒大小分布自动检测方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2194534A1 (en) * 1997-01-07 1998-07-07 Maztech Microvision Ltd. Method and apparatus for quantifying particle components

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4488245A (en) * 1982-04-06 1984-12-11 Loge/Interpretation Systems Inc. Method and means for color detection and modification
CN1258058A (zh) * 1998-12-14 2000-06-28 佳能株式会社 图像处理方法和装置,图像处理系统以及存储介质
CN101620060A (zh) * 2009-08-13 2010-01-06 上海交通大学 颗粒大小分布自动检测方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
朱秀昌等.离散K-L变换.《数字图像处理与图像通信》.2002,51-52. *
王桂棠等.基于数字图像处理技术的金相组织定量分析.《金属热处理》.2006,第31卷(第2期), *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103824268A (zh) * 2014-02-08 2014-05-28 江西赛维Ldk太阳能高科技有限公司 晶粒图像边缘连接方法及其装置
CN103824268B (zh) * 2014-02-08 2017-02-15 江西赛维Ldk太阳能高科技有限公司 晶粒图像边缘连接方法及其装置
CN108446706A (zh) * 2018-02-27 2018-08-24 西安交通大学 一种基于颜色主分量提取的磨粒材质自动识别方法
CN108446706B (zh) * 2018-02-27 2021-01-19 西安交通大学 一种基于颜色主分量提取的磨粒材质自动识别方法

Also Published As

Publication number Publication date
CN101964293A (zh) 2011-02-02

Similar Documents

Publication Publication Date Title
CN101964293B (zh) 一种金相显微组织图像处理方法
CN101193317B (zh) 对图像进行自动白平衡处理的方法和装置
CN103345617B (zh) 中药识别的方法及其系统
CN103034838B (zh) 一种基于图像特征的特种车辆仪表类型识别与标定方法
CN106294705B (zh) 一种批量遥感影像预处理方法
CN105572541A (zh) 一种基于视觉注意机制的高压巡线故障检测方法及系统
CN109118548A (zh) 一种综合智能水质识别方法
CN107144353A (zh) 一种基于数码相机的纺织品色差测量方法
CN106650606A (zh) 人脸图像的匹配及处理方法、人脸图像模型构建系统
CN110705639A (zh) 一种基于深度学习的医学精子图像识别系统
CN102300100A (zh) 一种白平衡控制装置及方法
CN101976338A (zh) 一种基于梯度方向直方图的判决式视觉显著性检测方法
CN114067222B (zh) 一种城市水体遥感分类方法及装置
Jin et al. Perceptual Gradient Similarity Deviation for Full Reference Image Quality Assessment.
CN103559486B (zh) 一种视频检测地板正反面的方法
CN112884795A (zh) 一种基于多特征显著性融合的输电线路巡检前景与背景分割方法
Crichton et al. Skin chromaticity gamuts for illumination recovery
CN103871084B (zh) 蓝印花布图案识别方法
CN116642827A (zh) 基于高光谱成像技术的大批量织物沾色评级方法及系统
CN111652055B (zh) 一种基于两阶段定位的开关仪表智能识别方法
CN111178229A (zh) 一种基于深度学习的静脉显像方法和装置
CN115508288A (zh) 一种基于掩码正态拟合的茶汤混浊度测定方法
CN116052041A (zh) 基于深度网络的指示灯状态识别方法
CN111007020B (zh) 双帧四光谱成像方法及应用
Ding et al. Rice lodging area extraction based on YCbCr spatial and texture features

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CP01 Change in the name or title of a patent holder

Address after: Weiyang Xu Jia Wan 710021 Shaanxi city of Xi'an Province

Patentee after: AECC AVIATION POWER CO,LTD.

Address before: Weiyang Xu Jia Wan 710021 Shaanxi city of Xi'an Province

Patentee before: AVIC AVIATION ENGINE Corp.,PLC

Address after: Weiyang Xu Jia Wan 710021 Shaanxi city of Xi'an Province

Patentee after: AVIC AVIATION ENGINE Corp.,PLC

Address before: Weiyang Xu Jia Wan 710021 Shaanxi city of Xi'an Province

Patentee before: XI'AN AVIATION POWER Co.,Ltd.

CP01 Change in the name or title of a patent holder