基于频率域处理的脉诊曲线图像网格消除方法
技术领域
[0001] 本
发明属于脉诊曲线图像数字化提取之前的背景网格线去除方法，特别涉及到一种在频率域处理的脉诊曲线图像网格消除方法。
背景技术
[0002] 对医学曲线的研究一直是生理学、病理学以及
生物医学工程等学科领域的重要工作。由于以前的医学曲线多是以纸张为介质并采用“档案式”管理模式，这种保存方式有几个主要的缺点，例如资料容易破损、色彩褪变，且存在管理繁琐等弊端，难以查询和维护，不便于与医院信息化
数据库连接和集成。为了解决这个问题，有必要对医学曲线进行数字化研究以为医院数字化、信息化
进程的发展奠定
基础。
[0003] 在进行脉诊曲线数字化时，我们需要对现有的纸质的脉诊曲线采用扫描或者
数码相机拍摄的方式获得脉诊曲线图像，图像背景通常为规整的网格，目标脉诊曲线与网格交叠重合，如果不把坐标网格去掉，将严重影响曲线提取的效果。所以对目标脉诊曲线进行检测与提取之前，首先有必要进行脉诊图片背景坐标网格的消除。背景坐标网格去除通常可以采用
空域去除和频域去除两种方法。空域去除法通常采用如均值滤波法、中值滤波法、维纳滤波法、十字模板匹配法、Hough曲变换法、投影法、直线模板匹配法、形态学滤波等方法。与空域处理不同的是，频域滤波主要针对在空域有一定变换规律的
信号。频域去除法是指通过傅里叶变换、余弦变换、
小波变换等变换方法将图像从空域变换到频域，在频域的
角度分析图像的特征并进行处理。频域法的优点是对背景网格依赖性不强，一般不会损失曲线信息。本发明提出了一种基于频率域处理的脉诊曲线图像网格消除方法。
发明内容
[0004] 为了在频率域消除脉诊图片的网格，有必要先对仿真网格的频域图像特点进行分析。本发明根据自制的仿真网格，分别对无较大倾斜的仿真网格和存在较大倾斜的仿真网格进行了实验，分析了网格在频域的分布特点，并得出了在频率域消除网格的具体步骤,为实现本发明之目的，采用以下技术方案予以实现：
[0005] (1)首先对脉诊曲线图像进行傅里叶变换并移频取对数后得到频率域图像并保存；
[0006] (2)对频率域图像进行二值化；
[0007] (3)在原始
频谱图像中，将与二值化后的对数频谱图的高亮部分所对应的
位置上的值
修改为0；
[0008] (4)最后对修改后的频率域图像进行傅里叶反变换并二值化，得到消除网格后的脉诊曲线图像。
[0009] 上述步骤1)、2)、3)、4)分别按如下方式进行：
[0010] 1)对原始脉诊曲线图像变换得到频谱图：对原始脉诊曲线图像进行快速傅里叶变换并移频取对数后得到对数频谱图；
[0011] 2)对频谱图进行二值化：对对数频谱图进行二值化；
[0012] 3)消除网格在频域中的分量：将对数频谱图按照二值化后的对数频谱图中的高亮部分所对应的位置进行修改，将对应位置上的值改为0，得到去网格后的频谱图；
[0013] 4)对去网格后的频率域图像进行反变换及后处理：对去网格后的频率域图像进行傅里叶反变换，对反变换后的图像进行二值化，然后用SUSAN滤波法对去网格后的二值图像中的点状噪声进行处理，得到消除网格后的脉诊曲线图像。
附图说明
[0015] 图2为快速傅里叶变换前后的空域脉诊图和频域脉诊图，图2a)为原始脉诊图，图2b)为对数频谱图
[0016] 图3为有无倾斜的仿真网格在空域和频域的一系列图像，图3a)为仿真网格原图，图3b)为仿真网格对数频谱图，图3c)为二值化的仿真网格对数频谱图，图3d)为倾斜仿真网格原图，图3e)为倾斜仿真网格对数频谱图，图3f)为二值化的倾斜仿真网格对数频谱图[0017] 图4为处理前后的对数频谱图，图4a)为二值频谱图，图4b)为处理后的频谱图[0018] 图5为有无二值化及后处理的去网格脉诊图，图5a)为初步去网格后的图像，图5b)为二值化并进行SUSAN滤波后的图像
具体实施方式
[0019] 为了更好地理解本发明的技术方案，以下结合图1-5作进一步说明并着重对其中所涉及的几个关键步骤进行较为详细的描述。
[0020] 1)频域脉诊曲线图像特点分析
[0021] 为了消除网格在频域的分量，我们必须找出网格在频域的分布规律。为此，通过自制的仿真网格实验，我们对无较大倾斜的网格和存在较大倾斜的网格分别进行了快速傅里叶变换(图2所示为快速傅里叶变换前后的空域脉诊图和频域脉诊图。其中，图2a)为原始脉诊图，图2b)为对数频谱图)。
[0022] 为了更清晰地观察到网格在频域的分布特点，我们对频谱图进行了二值化。由图3可以看出，不论仿真网格是否有倾斜，网格在频域中的强
能量点在频谱图中一个十字线上分布明显，只是十字叉存在角度的差异，但不管怎样我们都可以通过消除频域中对应规整网格的这些高能量点来达到消除网格的目的。
[0023] 2)二值化
[0024] 本发明对原始脉诊曲线图像进行快速傅里叶变换得到频谱图，再通过移频及对数变换得到对数频谱图，增加了频谱图的可视细节。为了更准确的得到频域中网格所对应的这些高能量点的位置，我们要对对数频谱图进行二值化。图4为处理前后的对数频谱图，图4a)为二值频谱图，图4b)为处理后的频谱图。
[0025] 然后将对数频谱图按照二值化后的对数频谱图中的高亮部分所对应的位置进行修改，将对应位置上的值改为0，得到处理后的频谱图。在消除了频谱图上的这些强能量点之后，再进行傅里叶反变换，就得到初步去网格后的图像。
[0026] 3)后处理
[0027] 为了得到更好的去网格效果，我们要对反变换后的去网格图像进行后处理。首先我们要进行二值化，然后用SUSAN滤波法对去网格后的二值图像中的点状噪声进行处理，得到很好的去网格的效果。如图5所示为所得到的去网格脉诊图。其中，图5a)为去网格后的图像，图5b)为二值化与SUSAN滤波后的图像。
