CN108921171A - 一种骨关节x线片自动识别分级方法 - Google Patents
一种骨关节x线片自动识别分级方法 Download PDFInfo
- Publication number
- CN108921171A CN108921171A CN201810648749.5A CN201810648749A CN108921171A CN 108921171 A CN108921171 A CN 108921171A CN 201810648749 A CN201810648749 A CN 201810648749A CN 108921171 A CN108921171 A CN 108921171A
- Authority
- CN
- China
- Prior art keywords
- bones
- joints
- ray film
- contour line
- key position
- 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
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/22—Matching criteria, e.g. proximity measures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/20—Image preprocessing
- G06V10/30—Noise filtering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/40—Extraction of image or video features
- G06V10/44—Local feature extraction by analysis of parts of the pattern, e.g. by detecting edges, contours, loops, corners, strokes or intersections; Connectivity analysis, e.g. of connected components
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H30/00—ICT specially adapted for the handling or processing of medical images
- G16H30/20—ICT specially adapted for the handling or processing of medical images for handling medical images, e.g. DICOM, HL7 or PACS
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V2201/00—Indexing scheme relating to image or video recognition or understanding
- G06V2201/03—Recognition of patterns in medical or anatomical images
- G06V2201/033—Recognition of patterns in medical or anatomical images of skeletal patterns
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Data Mining & Analysis (AREA)
- Health & Medical Sciences (AREA)
- Multimedia (AREA)
- Artificial Intelligence (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- General Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Evolutionary Computation (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Radiology & Medical Imaging (AREA)
- Epidemiology (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Primary Health Care (AREA)
- Public Health (AREA)
- Image Processing (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明提供一种骨关节X线片自动识别分级方法,该方法根据骨关节X线片图像模板库,对骨关节X线片图像样本分级,判断骨关节状况,其中提取骨关节X线片图像模板库和图像样本的轮廓线形状向量方法相同。本发明首先对骨关节X线片图像预处理,运用小波去噪、中值滤波、Log边缘检测算子法提取骨关节关键部位轮廓线。接着利用邻域特征算法截取关键部位轮廓线,并且采用连通区域标记法滤波得到优化的关键部位轮廓线。最后采用改进的傅里叶描述子算法,提取骨关节X线片图像样本的轮廓线形状向量,由分类算法匹配骨关节X线片图像模板库的形状向量,判断图像样本的骨关节状况。本发明的骨关节X线片自动识别分级方法可以取代人工识别,快速且精度高。
Description
技术领域
本发明涉及一种骨关节X线片自动识别分级方法,涉及骨科图像处理技术。
背景技术
骨关节X光片是临床检测骨关节状况的主要方式,对X光片目前是人工识别,根据现代医学的发展,采用自动方法进行骨关节分级,辅助医生监视骨关节状况十分有必要。
现代骨科医疗自动判别领域是一个很有应用前景的方向,但是通过医师观看骨关节X线片来进行人工分级,这样往往过于主观,不同人的分级结果难免会存在差异,且工作量很大。骨科X光片自动识别目前采用的是对图像二值化处理,进行相关识别。但是针对对象不能直接二值化的情况,由于图像失真度非常大,骨关节的轮廓提取精度要求高,这样导致X光片自别识别困难。并且,由于划分的级别多,而级与级之间差别又不是很大,有些级别之间甚至肉眼都无法区分,必须开发新算法对骨关节这个关键部位进行截取放大才能保证识别精度,因此,目前来说实现骨关节X光片自动分级仍然是难点。
发明内容
鉴于以上所述针对现有技术的空缺与不足,本发明的目的在于提供一种高精度高效率骨关节X线片的自动识别分级方法。
本发明提供一种骨关节X线片自动识别分级方法,特征在于,根据骨关节X线片图像模板库,对骨关节X线片图像样本分级,判断骨关节状况,其中提取骨关节X线片图像模板库和图像样本的轮廓线形状向量方法相同,本发明包括如下步骤:
1)运用小波去噪、中值滤波算法对骨关节X线片进行滤波处理;
2)使用Log边缘检测算子法提取骨关节X线片轮廓线;
3)用邻域特征算法截取关键部位轮廓线,并且采用连通区域标记法优化;
4)使用改进的傅里叶描述子算法,提取关键部位轮廓线形状向量;
5)匹配骨关节X线片图像模板库和图像样本的形状向量,对骨关节X线片图像样本分级。
优选的,步骤3)采用邻域特征算法截取关键部位轮廓线,其具体包括如下步骤:
1)确定骨关节X线片轮廓线数据矩阵中每行白像素个数h(i),其中i表示行数;i=1,2,3,L;
2)确定i,使得h(i)满足:h(i)≥H(H为阈值,可根据实际确定),且h(i-1)<h(i)<h(i+1);
3)定义截取高度k1,将以外的行中的白像素替换为黑像素;
4)将所有不含白像素的行和列截去,提取出关节间隙部分。
优选的,步骤4)使用改进的傅里叶描述子算法,对关键部位轮廓线进行特征向量提取,其具体包括如下步骤:
1)设定描述子的大小;
2)设置骨关节X线片关键部位轮廓线的偏移中心;
3)将提取的骨关节X线片关键部位轮廓线傅里叶展开,二维数据转化成一维特征向量数据。
提取骨关节部位十分重要,采用上述的算法可以精确定位骨关节部分并将其截取出来,这步操作直接决定骨关节形状特征识别精度,十分关键。傅里描述子只能用于封曲线,当用于非封闭曲线其精确度十分低下,这里改进的傅里叶描述子算法可以用于非封闭曲线,由于经过关键部分提取,图像基本在一个很小的区域内,只要设定取点中心和提取的描述子大小范围,结果表明可行度很高。因此,采用上述方法后,识别精度和效率将远远高于人工方法,使用本发明算法就可大大降低误判率。
附图说明
图1为本发明自动识别分级方法流程图。
图2为本发明实施例骨关节X线片图像模板库。(KL表示骨关节状况等级,KL=0表示等级为0;KL=1表示等级为1;KL=2表示等级为2;KL=3表示等级为3;KL=4表示等级为4)
图3为本发明实施例预处理图。
图4为本发明实施例关键部位优化轮廓线。
具体实施方式
以下通过特定的具体实例说明本发明的实施方式,本领域技术人员可由本说明书所揭露的内容轻易地实施。
对于骨关节状况的分级问题,只需考虑关节部位所能体现的信息。自动识别分级方法流程如图1所示,骨关节X线片图像模板库如图2所示,本实施例采用5个模板。注意:提取骨关节X线片图像模板库和图像样本的轮廓线形状向量方法相同。对X线片自动分级主要分成以下5步。
步骤1运用小波去噪、中值滤波滤波算法对骨关节X线片进行滤波处理;
在这一步当中,主要有2个部分,即小波去噪和中值滤波。
小波去噪方法是:通过小波变换将目标图像实现小波级数展开,分成基波和谐波信号;对分解后的谐波系数进行阈值处理;最后将阈值处理后谐波系数和未处理的基波系数逆变回图像即重构。其中主要是阈值获取和阈值去噪两方面,本发明实施例的阈值获取采用ddencmp函数,阈值去噪采用wdencmp函数。
MATLAB对应程序段如下:
[thr,sorh,keepapp]=ddencmp(a,b,c);
S2=wdencmp(a,b,‘wname’,N,THR,sorh,keepapp);
ddencmp是阈值获取函数,c为要求处理的输入图像矩阵数据,a分为两种情况,当其为den时表示实现去噪,当其为cmp时表示实现阈值返回,sorh是软硬阈值参数选择,其余参数由于本发明不使用。
wdencmp函数是在获取阈值后的基础上实现去噪,a可取gbl和lvd,前者表示小波分解的不同层之间使用同一个阈值,后者表示不同层之间采用不同的阈值,wname是对应的小波变换算法,THR是阈值向量。sorh是从ddencmp函数获取的处理参数,keepapp是决定是否阈值量化的参数,只为值为1时进行,值为0时不进行。S2为经过小波去噪后的输出图像矩阵。
中值滤波是把数字图像或数字序列中一点的值用该点的一个邻域中各点值的中值代替。在MATLAB中使用medfilt2函数实现中值滤波操作:
S2=medfilt2(S1,[a,b]);
其中S2为中值滤波输出图像矩阵数据,S1为要处理的输入图像矩阵数据,a、b为指定输出图像大小。
步骤2使用Log边缘检测算子法提取骨关节X线片轮廓线;
运用Log边缘检测算子法等算法提取骨关节轮廓线,在MATLAB中实现Log边缘检测算法只需要调用edge函数即可,程序如下所示:
S2=edge(S1,'log');
其中S1是要求处理的图像矩阵数据,log表示采用拉普拉斯算子进行边缘提取,S2是经过拉普拉斯算子处理后的输出矩阵数据,对骨关节X线片图像样本轮廓线处理结果如图3所示。
步骤3用邻域特征算法截取关键部位轮廓线,并且采用连通区域标记法优化;
采用邻域特征算法截取关键部位轮廓线,其具体包括如下步骤:
1)确定骨关节X线片轮廓线数据矩阵中每行白像素个数h(i),其中i表示行数;i=1,2,3,L;
2)确定i,使得h(i)满足:h(i)≥H(H为阈值,可根据实际确定),且h(i-1)<h(i)<h(i+1);
3)定义截取高度k1,将以外的行中的白像素替换为黑像素;
4)将所有不含白像素的行和列截去,提取出关节间隙部分。
采用连通区域标记法优化关节部位轮廓线,本发明采用八邻域来作为连通区域标记法识别,过程如下:
1)判断此点八邻域中的最左,左上,最上,上右点的情况。如果都没有点,则表示一个新的区域的开始。
2)如果此点八邻域中的最左有点,上右都有点,则标记此点为这两个中的最小的标记点,并修改大标记为小标记。
3)如果此点八邻域中的左上有点,上右都有点,则标记此点为这两个中的最小的标记点,并修改大标记为小标记。
4)否则按照最左,左上,最上,上右的顺序,标记此点为四个中的一个。
MATLAB使用bwlabel函数实现连通区域标记法,程序段如下:
[S2,n]=bwlabel(S1,N);
其中N的取值可为8,表示采用八邻域来作为连通区域标记法识别核心。S1为要求处理图像矩阵数据,一般是经过滤波操作甚至边缘预处理(边缘提取)后的图像矩阵数据,S2为经过连通区域标记法处理后的输出矩阵。n为返回连通数,本发明没有使用。
根据骨关节间隙的位置特征截取关节部位,同时采用连通区域标记法就能得到理想的处理后的图像。用邻域特征算法截取关键部位优化轮廓线,并且采用连通区域标记法优化结果如图4所示。
步骤4使用改进的傅里叶描述子算法,对关键部位轮廓线进行特征向量提取;
1)使用getSampleFD从五个图像样本中获取傅里叶描述子,设定描述子的大小。其中提取傅立叶描述符,对傅立叶描述符进行规范化,调整傅立叶描述符的大小,以及计算区别分别采用extractFD(),normaliseFD(),resizeFD(),compareFD()函数实现。
2)设置骨关节X线片关键部位轮廓线的偏移中心;
3)将提取的骨关节X线片关键部位轮廓线傅里叶展开,二维数据转化成一维特征向量数据。
步骤5匹配骨关节X线片图像模板库和样本的形状向量,对骨关节X线片图像样本分级。
1)提取模板库5个图像特征向量矩阵;
2)提取图像样本图像特征向量矩阵;
3)利用abs()函数计算以上二者元素之间的差的绝对值之和;
4)求出图像样本与模板库五者每个元素之间差的绝对值之和最小值;
5)记下最小值的模板标号,匹配出最合适的模板,对骨关节X线片图像样本分级。
将图3中的5幅图分别命名为O,A,B,C,D,对应等级0,1,2,3,4。数据结果如表1所示。
表1实验数据统计表
从结果来看,本发明骨关节X线片自动识别分级方法是准确的,而且识别精度和效率将远远高于人工方法,可大大降低误判率。本发明有效克服了现有技术中的种种缺点而具高度实用价值。
Claims (3)
1.一种骨关节X线片自动识别分级方法,特征在于,根据骨关节X线片图像模板库,对骨关节X线片图像样本分级,判断骨关节状况,其中提取骨关节X线片图像模板库和图像样本的轮廓线形状向量方法相同,本发明包括如下步骤:
1)运用小波去噪、中值滤波算法对骨关节X线片进行滤波处理;
2)使用Log边缘检测算子法提取骨关节X线片轮廓线;
3)用邻域特征算法截取关键部位轮廓线,并且采用连通区域标记法优化;
4)使用改进的傅里叶描述子算法,提取关键部位轮廓线形状向量;
5)匹配骨关节X线片图像模板库和图像样本的形状向量,对骨关节X线片图像样本分级。
2.根据权利要求1所述的骨关节X线片自动识别分级方法,其特征在于:步骤3)采用邻域特征算法截取关键部位轮廓线,其具体包括如下步骤:
1)确定骨关节X线片轮廓线数据矩阵中每行白像素个数h(i),其中i表示行数;i=1,2,3,L;
2)确定i,使得h(i)满足:h(i)≥H(H为阈值,可根据实际确定),且h(i-1)<h(i)<h(i+1);
3)定义截取高度k1,将以外的行中的白像素替换为黑像素;
4)将所有不含白像素的行和列截去,提取出关节间隙部分。
3.根据权利要求1所述的骨关节X线片自动识别分级方法,其特征在于:步骤4)使用改进的傅里叶描述子算法,对关键部位轮廓线进行特征向量提取,其具体包括如下步骤:
1)设定描述子的大小;
2)设置骨关节X线片关键部位轮廓线的偏移中心;
3)将提取的骨关节X线片关键部位轮廓线傅里叶展开,二维数据转化成一维特征向量数据。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810648749.5A CN108921171B (zh) | 2018-06-22 | 2018-06-22 | 一种骨关节x线片自动识别分级方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810648749.5A CN108921171B (zh) | 2018-06-22 | 2018-06-22 | 一种骨关节x线片自动识别分级方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108921171A true CN108921171A (zh) | 2018-11-30 |
CN108921171B CN108921171B (zh) | 2022-02-08 |
Family
ID=64421013
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810648749.5A Active CN108921171B (zh) | 2018-06-22 | 2018-06-22 | 一种骨关节x线片自动识别分级方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108921171B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109886320A (zh) * | 2019-01-30 | 2019-06-14 | 姜姿君 | 一种人体股骨x光片智能识读方法及系统 |
CN113763340A (zh) * | 2021-08-25 | 2021-12-07 | 中国人民解放军空军军医大学 | 基于多任务深度学习强直性脊柱炎的自动分级方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100177946A1 (en) * | 2007-05-18 | 2010-07-15 | Marleen De Bruijne | Semi-automatic contour detection |
US20110091089A1 (en) * | 2008-04-08 | 2011-04-21 | Yamaguchi University | X-ray Image Processing Method and X-ray Image Processing System Using Open Curve Fourier Descriptors |
CN103247056A (zh) * | 2013-05-29 | 2013-08-14 | 中国人民解放军第三军医大学第一附属医院 | 人体骨关节系统三维模型-二维影像空间配准方法 |
CN105389802A (zh) * | 2015-10-20 | 2016-03-09 | 无锡日联科技股份有限公司 | 一种基于X-Ray图像的IC元器件缺陷检测方法 |
CN105512629A (zh) * | 2015-12-04 | 2016-04-20 | 南方医科大学 | 一种人体手指静脉图像骨关节定位方法 |
CN107368792A (zh) * | 2017-06-29 | 2017-11-21 | 广州智慧城市发展研究院 | 一种基于滤波器以及主干边缘的指静脉识别方法及系统 |
-
2018
- 2018-06-22 CN CN201810648749.5A patent/CN108921171B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100177946A1 (en) * | 2007-05-18 | 2010-07-15 | Marleen De Bruijne | Semi-automatic contour detection |
US20110091089A1 (en) * | 2008-04-08 | 2011-04-21 | Yamaguchi University | X-ray Image Processing Method and X-ray Image Processing System Using Open Curve Fourier Descriptors |
CN103247056A (zh) * | 2013-05-29 | 2013-08-14 | 中国人民解放军第三军医大学第一附属医院 | 人体骨关节系统三维模型-二维影像空间配准方法 |
CN105389802A (zh) * | 2015-10-20 | 2016-03-09 | 无锡日联科技股份有限公司 | 一种基于X-Ray图像的IC元器件缺陷检测方法 |
CN105512629A (zh) * | 2015-12-04 | 2016-04-20 | 南方医科大学 | 一种人体手指静脉图像骨关节定位方法 |
CN107368792A (zh) * | 2017-06-29 | 2017-11-21 | 广州智慧城市发展研究院 | 一种基于滤波器以及主干边缘的指静脉识别方法及系统 |
Non-Patent Citations (1)
Title |
---|
韦巍 等: "骨科X线图像的边缘检测", 《中国影像学技术》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109886320A (zh) * | 2019-01-30 | 2019-06-14 | 姜姿君 | 一种人体股骨x光片智能识读方法及系统 |
CN109886320B (zh) * | 2019-01-30 | 2020-04-21 | 姜姿君 | 一种人体股骨x光片智能识读方法及系统 |
CN113763340A (zh) * | 2021-08-25 | 2021-12-07 | 中国人民解放军空军军医大学 | 基于多任务深度学习强直性脊柱炎的自动分级方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108921171B (zh) | 2022-02-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP3176751B1 (en) | Information processing device, information processing method, computer-readable recording medium, and inspection system | |
CN109919960B (zh) | 一种基于多尺度Gabor滤波器的图像连续边缘检测方法 | |
CN108154519A (zh) | 眼底图像中血管的分割方法、装置及存储介质 | |
WO2013099772A1 (ja) | 細胞輪郭線形成装置及びその方法、コンピュータにより処理可能な細胞輪郭線形成プログラムを記憶する記憶媒体 | |
Valliammal et al. | Plant leaf segmentation using non linear K means clustering | |
CN110111387B (zh) | 一种基于表盘特征的指针表定位及读数方法 | |
CN112308854B (zh) | 一种芯片表面瑕疵的自动检测方法、系统及电子设备 | |
CN113177915B (zh) | 一种胫骨平台后倾角测量方法及装置、存储介质 | |
Wang et al. | Leaf recognition based on elliptical half Gabor and maximum gap local line direction pattern | |
CN115797356A (zh) | 一种核磁共振肿瘤区域提取方法 | |
CN108921171A (zh) | 一种骨关节x线片自动识别分级方法 | |
CN111105427A (zh) | 一种基于连通区域分析的肺部图像分割方法及系统 | |
Parthasarathy et al. | Predicting source and age of brain tumor using canny edge detection algorithm and threshold technique | |
Hazarika et al. | A new breast border extraction and contrast enhancement technique with digital mammogram images for improved detection of breast cancer | |
US12106334B2 (en) | Artificial intelligence-based system and method for grading collectible trading cards | |
CN112801031A (zh) | 静脉图像识别方法、装置、电子设备及可读存储介质 | |
Soni et al. | CT scan based brain tumor recognition and extraction using Prewitt and morphological dilation | |
CN117314940B (zh) | 基于人工智能的激光切割零件轮廓快速分割方法 | |
CN112258532B (zh) | 一种超声图像中胼胝体的定位及分割方法 | |
CN110647889A (zh) | 医学图像识别方法、医学图像识别装置、终端设备及介质 | |
EP3764319B1 (en) | Image comparison device | |
CN107729863B (zh) | 人体指静脉识别方法 | |
MUSTAFA et al. | A comparative study of automated segmentation methods for cell nucleus detection | |
Vyshnavi et al. | Breast density classification in mammogram images | |
CN106295478A (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 |