CN107220644A - 一种心电图扫描图像倾斜度校正方法 - Google Patents
一种心电图扫描图像倾斜度校正方法 Download PDFInfo
- Publication number
- CN107220644A CN107220644A CN201710251629.7A CN201710251629A CN107220644A CN 107220644 A CN107220644 A CN 107220644A CN 201710251629 A CN201710251629 A CN 201710251629A CN 107220644 A CN107220644 A CN 107220644A
- Authority
- CN
- China
- Prior art keywords
- image
- point
- ecg
- represent
- value
- 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
- 238000000034 method Methods 0.000 claims abstract description 27
- 238000012545 processing Methods 0.000 claims abstract description 19
- 230000002146 bilateral effect Effects 0.000 claims abstract description 6
- 238000001914 filtration Methods 0.000 claims abstract description 6
- 230000010339 dilation Effects 0.000 claims abstract description 4
- 230000000877 morphologic effect Effects 0.000 claims abstract description 4
- 239000004576 sand Substances 0.000 claims description 2
- 238000010276 construction Methods 0.000 abstract description 4
- 230000002708 enhancing effect Effects 0.000 abstract description 2
- 238000001514 detection method Methods 0.000 description 6
- 238000004422 calculation algorithm Methods 0.000 description 5
- 238000006243 chemical reaction Methods 0.000 description 4
- 238000012937 correction Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000007639 printing Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 238000003759 clinical diagnosis Methods 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 238000002405 diagnostic procedure Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- 238000004171 remote diagnosis Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Classifications
-
- 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/24—Aligning, centring, orientation detection or correction of the image
- G06V10/245—Aligning, centring, orientation detection or correction of the image by locating a pattern; Special marks for positioning
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Multimedia (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明涉及一种心电图扫描图像倾斜度校正方法,包括:双边滤波处理;从RGB颜色空间转换到HSV颜色空间;将饱和度分量图FS进行增强处理;使用Otsu法进行二值化构造二值图;形态学膨胀处理;计算二值图中取值为1的各点的梯度值和方向;构造梯度方向直方图HOG;得到局部极值点序列{HSTm},并选取接近0°的第一主方向角和接近90°的第二主方向角;计算心电图扫描图像的倾斜角;得到校正后的图像。
Description
技术领域
本发明涉及心电图扫描图像的数字化处理技术,尤其是针对心电图扫描图像的倾斜度校正方法。
背景技术
在临床诊断过程中,医院一直依赖心电图作为监测病人心脏电活动的仪器。由于纸张易破碎且热敏纸字迹不稳定、易消退,多数纸质心电图会遭到了不同程度的破坏。而这些纸质心电图是一笔宝贵的财富,可以建立丰富的病例数据库、有助于学者进行科研分析、增加医生临床诊断经验以及支持在线远程诊断等。因此纸质心电图数字化信息的提取成为一个亟待解决的问题。
纸质心电图是由纵横交错的背景坐标网格、打印在网格上的心电波形曲线、导联字符以及定标电压组成。纸质心电图数字化信息提取,不是单纯指将纸质心电图通过扫描仪等设备转换成数字图片保存,还包括通过一系列的智能处理技术将纸面上的波形数据提取并保存的数字化过程。使用扫描仪将纸质心电图转换为数字图像,是完成纸质心电图数字化过程的第一步。但是由于人工操作或者扫描设备等原因,得到的扫描心电图像通常存在不同程度的倾斜,这对于后续的数字化过程,特别是心电图信号定标和量化过程不利。因此有必要在心电图信号数字化实现过程中,检测心电图扫描图像的倾斜角度,并做相应校正。
目前针对扫描图像的倾斜角度检测的方法有:基于投影的方法、基于变换的方法(如常用Hough变换或Rodon变换)和基于特征的方法。基于投影的方法简单易行,处理速度快,但对于密布网格线的心电图图像效果不佳。基于变换的方法通过检测直线段来判断倾斜角,但无论是Hough变换还是Rodon变换,所需的计算量较大,如果输入的扫描图像分辨率较高,将导致处理速度不能满足实时性要求。相对于前两种方法,基于特征的方法,可以快速准确地检测扫描图像的倾斜角。
发明内容
本发明的目的是提供一种针对心电图扫描图像的基于特征的倾斜度校正方法,本发明能够快速准确检测心电图扫描图像的倾斜角。技术方案如下:
一种心电图扫描图像倾斜度校正方法,包括下列步骤:
步骤1:将输入彩色扫描ECG图像进行双边滤波处理,处理结果用F表示;
步骤2:将F从RGB颜色空间转换到HSV颜色空间,用FH、FS和FV分别表示色调、饱和度和强度分量图;
步骤3:为了加大背景与网格点的差别,将饱和度分量图FS进行增强处理,包括对比度拉伸和归一化处理,处理结果用FSE表示;
步骤4:使用Otsu法计算FSE的全局阈值TH1,构造二值图BW1,其中BW1中取值为1对应于FSE中取值大于TH1的点,其它点取值为0;
步骤5:使用半径为2的碟状算子,对BW1进行形态学膨胀处理,处理结果二值图用BW2表示;
步骤6:使用强度分量FV,计算BW2中取值为1的各点的梯度值和方向,分别用GRDi和DIRi表示,其中下标i表示序号;
步骤7:计算梯度值的最大值,用GRDmax表示,设定一个小于0.5的α值,只选取满足条件GRDi>α×GRDmax的点对应的梯度方向值来构造梯度方向直方图HOG,;
步骤8:用HSTj(j=0,1,...,180)表示直方图数值,选取其中的局部极值点,并用{HSTm}表示局部极值点序列;
步骤9:从局部极值点序列{HSTm}中选取接近0°的第一主方向角θ1和接近90°的第二主方向角θ2;
步骤10:计算心电图扫描图像的倾斜角φ=0.5×(|θ1-0|+|θ2-90|);
步骤11:根据倾斜角φ旋转图像,得到校正后的图像。
计算机仿真结果表明,使用本发明所提方法能够快速准确检测心电图扫描图像的倾斜角,能够满足实时处理的要求。
附图说明
图1是心电图扫描图像构成示意图
图2所提方法流程图
图3所提方法处理结果示意图,(a)为原图(b)为处理结果图
具体实施方式
目前国内大多数医院打印心电图所用的纸张通常是热敏纸。该纸张使用纵横交错的红色网格来定标,心电波形则以黑色波形打印在纸上。扫描彩色心电图图像中主要包括三类颜色点(如图1所示):(1)黑色点,包括心电图波形和标注字符;(2)红色点,主要是背景网格点,(3)白色背景点。三类颜色点中,白色背景点数目最多,数目较少的黑色点主要用于提取心电图波形信息,而红色网格点可用于帮助检测图像旋转角度。红色网格点由许多相互垂直的水平线和垂直线构成,其中水平线的旋转角与扫描图像的旋转角基本一致。本发明所提方法利用水平网格线的梯度方向确定心电图扫描图像的倾斜角。
根据心电图图像背景分布着密集、规则的红色网格的特点,本发明提出一种心电图扫描图像倾斜度校正方法。首先,使用双边滤波对心电图扫描图像做平滑处理,将转换到HSV颜色空间;然后对S分量进行增强处理,并提取网格区域;再使用V分量计算网格区域内各点梯度的方向值,构造网格点区域对应的梯度方向直方图(HOG);最后,根据得到的HOG计算倾斜度角,完成心电图扫描图像的倾斜度校正。图2所示为本发明所提方法的框图。具体流程如下;
1、双边滤波
纸质心电图图像在打印、存储、扫描等过程中都有可能引入噪声,这些噪声对于倾斜度校正和后续的数字化过程都会造成干扰,因此本发明选用双边滤波技术,对输入的彩色心电图扫描图像进行滤波处理,处理结果用F表示。
2、颜色通道分离
为了提取位于背景中的红色网格点,将F由RGB空间转换到HSV空间,分别用FH、FS和FV表示色调、饱和度和强度三个分量图像。
3、网格点提取
在心电图扫描图像中,黑色或者深灰色的心电图波形,以及白色的背景区域,两类区域像素点的饱和度值都较小,相比较而言,红色网格线各点的饱和度数值较高。因此,可以使用饱和度值来区分网格点和背景点。具体的做法是:
算法1:网格点提取算法
第一步:对FS进行对比度拉伸和归一化处理,目的是加大背景与网格点的差别,处理结果用FSE表示。
第二步:使用Otsu法计算FSE的全局阈值TH1,构造二值图BW1,其中BW1中取值为1对应于FSE中取值大于TH1的点,其它点取值为0。
第三步:使用半径为2的碟状算子,对BW1进行形态学膨胀处理,处理结果二值图用BW2表示。
4、主梯度方向检测
使用强度分量FV,计算BW2中取值为1的点的梯度方向值(0°到180°之间),构造梯度方向直方图,用HOG表示。由于网格线是由相互垂直的线构成,因此对于不存在倾斜的心电图扫描图像,其HOG中0°和90°的取值应明显大于其他角度,形成两个极值,我们称为两个主方向。
当扫描图像发生倾斜时,两个主方向角仍然存在,且与0°和90°的偏差反映了倾斜程度。基于上述分析,本算法使用以下算法确定倾斜度角:
算法2:倾斜角检测算法
第一步:计算BW2中取值为1的各点的梯度数值和方向,分别用GRDi和DIRi表示,其中下标i表示序号。
第二步:计算梯度值的最大值,用GRDmax表示,只选取满足条件GRDi>α×GRDmax的点对应的梯度方向值来计算HOG,直方图中角度的取值范围是[0°,180°],间隔为1°,令α=0.2。
第三步:用HSTj(j=0,1,...,180)表示直方图数值,下标j表示对应的角度。对于某一个梯度方向值HSTk,如果满足HSTk的值是{HSTk-l,HSTk-l+1,...,HSTk-1,HSTk,HSTk+1,...,HSTk+l-1,HSTk+l}(l=15)中的最大值时,才将HSTk判为局部极值点。用{HSTm}表示局部极值点序列。
第四步:从局部极值点{HSTm}中选取最大值点和次大值点,计算它们对应的角度值,分别用θ1和θ2表示,如果满足85°<|θ1-θ2|<95°,则将θ1和θ2判为主方向角1和主方向角2,并用θ1表示接近0°的主方向角,用θ2表示接近90°的主方向角。
第五步:使用下式计算心电图扫描图像的倾斜角φ(假定顺时针旋转为正方向),
φ=0.5×(|θ1-0|+|θ2-90|) (1)
5、倾斜角校正
根据倾斜角φ,逆时针方向旋转φ,得到倾斜度校正后的图像,并用白色填充旋转后图像中数值没有定义的区域。
采用Windows10系统下的matlab2015b作为实验仿真平台。选用专利申请人从项目合作医院得到50幅ECG扫描图像作为测试集。采用本发明提出的方法对测试图像进行处理,得到了良好的处理效果。对于1750×1275大小的图像,采用所提方法的处理速度平均为35ms,处理速度非常迅速。图2给出了部分处理结果图像,其中左侧为输入的倾斜图像,右侧为处理后的图像。由实验结果可见,采用本发明所提方法,可以有效检测并校正ECG扫描图像的倾斜情况,并且能够满足实时处理的要求。
Claims (1)
1.一种心电图扫描图像倾斜度校正方法,包括下列步骤:
步骤1:将输入彩色扫描ECG图像进行双边滤波处理,处理结果用F表示;
步骤2:将F从RGB颜色空间转换到HSV颜色空间,用FH、FS和FV分别表示色调、饱和度和强度分量图;
步骤3:为了加大背景与网格点的差别,将饱和度分量图FS进行增强处理,包括对比度拉伸和归一化处理,处理结果用FSE表示;
步骤4:使用Otsu法计算FSE的全局阈值TH1,构造二值图BW1,其中BW1中取值为1对应于FSE中取值大于TH1的点,其它点取值为0;
步骤5:使用半径为2的碟状算子,对BW1进行形态学膨胀处理,处理结果二值图用BW2表示;
步骤6:使用强度分量FV,计算BW2中取值为1的各点的梯度值和方向,分别用GRDi和DIRi表示,其中下标i表示序号;
步骤7:计算梯度值的最大值,用GRDmax表示,设定一个小于0.5的α值,只选取满足条件GRDi>α×GRDmax的点对应的梯度方向值来构造梯度方向直方图HOG;
步骤8:用HSTj(j=0,1,...,180)表示直方图数值,选取其中的局部极值点,并用{HSTm}表示局部极值点序列;
步骤9:从局部极值点序列{HSTm}中选取接近0°的第一主方向角θ1和接近90°的第二主方向角θ2;
步骤10:计算心电图扫描图像的倾斜角φ=0.5×(|θ1-0|+|θ2-90|);
步骤11:根据倾斜角φ旋转图像,得到校正后的图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710251629.7A CN107220644B (zh) | 2017-04-18 | 2017-04-18 | 一种心电图扫描图像倾斜度校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710251629.7A CN107220644B (zh) | 2017-04-18 | 2017-04-18 | 一种心电图扫描图像倾斜度校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107220644A true CN107220644A (zh) | 2017-09-29 |
CN107220644B CN107220644B (zh) | 2020-04-24 |
Family
ID=59927892
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710251629.7A Expired - Fee Related CN107220644B (zh) | 2017-04-18 | 2017-04-18 | 一种心电图扫描图像倾斜度校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107220644B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109636744A (zh) * | 2018-11-26 | 2019-04-16 | 山东航天电子技术研究所 | 一种水下距离选通成像的自适应图像处理方法 |
CN110507311A (zh) * | 2019-08-27 | 2019-11-29 | 中科麦迪人工智能研究院(苏州)有限公司 | 一种基于图像信息的心电分析方法、装置、设备及介质 |
CN114663443A (zh) * | 2022-02-24 | 2022-06-24 | 清华大学 | 12导联纸质心电图数字化方法与装置 |
CN115579109A (zh) * | 2022-11-24 | 2023-01-06 | 合肥心之声健康科技有限公司 | 医疗环境下心电图图像分析方法、装置和终端设备 |
CN117078913A (zh) * | 2023-10-16 | 2023-11-17 | 第六镜科技(成都)有限公司 | 对象倾斜矫正方法、装置、电子设备和存储介质 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101447017A (zh) * | 2008-11-27 | 2009-06-03 | 浙江工业大学 | 一种基于版面分析的选票快速识别统计方法及系统 |
CN102038498A (zh) * | 2011-01-20 | 2011-05-04 | 天津理工大学 | 一种纸质心电图数字化方法 |
CN104103047A (zh) * | 2014-07-25 | 2014-10-15 | 上海理工大学 | 一种心电图像倾斜度校正方法 |
US20140323850A1 (en) * | 2013-04-24 | 2014-10-30 | Siemens Aktiengesellschaft | Method and ekg trigger device for correcting an ekg signal in magnetic resonance image acquisition |
CN106097254A (zh) * | 2016-06-07 | 2016-11-09 | 天津大学 | 一种扫描文档图像纠偏方法 |
CN106447615A (zh) * | 2016-08-02 | 2017-02-22 | 浪潮软件股份有限公司 | 一种扫描文档图像的倾斜校正方法 |
-
2017
- 2017-04-18 CN CN201710251629.7A patent/CN107220644B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101447017A (zh) * | 2008-11-27 | 2009-06-03 | 浙江工业大学 | 一种基于版面分析的选票快速识别统计方法及系统 |
CN102038498A (zh) * | 2011-01-20 | 2011-05-04 | 天津理工大学 | 一种纸质心电图数字化方法 |
US20140323850A1 (en) * | 2013-04-24 | 2014-10-30 | Siemens Aktiengesellschaft | Method and ekg trigger device for correcting an ekg signal in magnetic resonance image acquisition |
CN104103047A (zh) * | 2014-07-25 | 2014-10-15 | 上海理工大学 | 一种心电图像倾斜度校正方法 |
CN106097254A (zh) * | 2016-06-07 | 2016-11-09 | 天津大学 | 一种扫描文档图像纠偏方法 |
CN106447615A (zh) * | 2016-08-02 | 2017-02-22 | 浪潮软件股份有限公司 | 一种扫描文档图像的倾斜校正方法 |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109636744A (zh) * | 2018-11-26 | 2019-04-16 | 山东航天电子技术研究所 | 一种水下距离选通成像的自适应图像处理方法 |
CN110507311A (zh) * | 2019-08-27 | 2019-11-29 | 中科麦迪人工智能研究院(苏州)有限公司 | 一种基于图像信息的心电分析方法、装置、设备及介质 |
CN110507311B (zh) * | 2019-08-27 | 2022-07-19 | 中科麦迪人工智能研究院(苏州)有限公司 | 一种基于图像信息的心电分析方法、装置、设备及介质 |
CN114663443A (zh) * | 2022-02-24 | 2022-06-24 | 清华大学 | 12导联纸质心电图数字化方法与装置 |
CN115579109A (zh) * | 2022-11-24 | 2023-01-06 | 合肥心之声健康科技有限公司 | 医疗环境下心电图图像分析方法、装置和终端设备 |
CN117078913A (zh) * | 2023-10-16 | 2023-11-17 | 第六镜科技(成都)有限公司 | 对象倾斜矫正方法、装置、电子设备和存储介质 |
CN117078913B (zh) * | 2023-10-16 | 2024-02-02 | 第六镜科技(成都)有限公司 | 对象倾斜矫正方法、装置、电子设备和存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN107220644B (zh) | 2020-04-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107220644A (zh) | 一种心电图扫描图像倾斜度校正方法 | |
Emre Celebi et al. | Border detection in dermoscopy images using statistical region merging | |
JP2021518966A (ja) | 医用画像認識方法及びシステム、並びに、モデルトレーニング方法、コンピュータ装置、及びプログラム | |
Chanwimaluang et al. | An efficient blood vessel detection algorithm for retinal images using local entropy thresholding | |
CN109816742B (zh) | 基于全连接卷积神经网络的锥束ct几何伪影去除方法 | |
CN107174232B (zh) | 一种心电图波形提取方法 | |
JP4303598B2 (ja) | 画素の符号化法、画像処理法と、1以上の画素により再現されるオブジェクトの定性的な認識を目的とする画像処理法 | |
CN108830149B (zh) | 一种目标细菌的检测方法及终端设备 | |
CN110110723B (zh) | 一种图像中目标区域自动提取的方法及装置 | |
CN107133933B (zh) | 基于卷积神经网络的乳腺x线图像增强方法 | |
JP5709216B2 (ja) | 画像処理プログラム、方法及び装置 | |
CN110751605A (zh) | 一种图像处理方法、装置和电子设备及可读存储介质 | |
CN107564048A (zh) | 基于分叉点特征配准方法 | |
Cueva et al. | Detection of skin cancer” Melanoma” through computer vision | |
Garg et al. | ECG paper records digitization through image processing techniques | |
CN113889238B (zh) | 一种图像识别方法、装置、电子设备及存储介质 | |
Kumari et al. | Blood vessel extraction using wiener filter and morphological operation | |
CN114066846A (zh) | 一种基于深度学习的ctp非急性闭塞缺血评估方法和系统 | |
Cai et al. | Perception preserving decolorization | |
Allegra et al. | Virtual unrolling using X-ray computed tomography | |
Ruslau et al. | Edge detection of digital image with different edge types | |
Pearl Mary et al. | Unified adaptive framework for contrast enhancement of blood vessels | |
Poostchi et al. | Diabetic retinopathy dark lesion detection: preprocessing phase | |
CN114663443A (zh) | 12导联纸质心电图数字化方法与装置 | |
CN110390339B (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200424 |