一种利用线对卡自动计算成像平板探测器线对比度的方法技术

本发明专利技术涉及图像处理技术领域，具体公开了一种利用线对卡自动计算成像平板探测器线对比度的方法，包括：通过成像平板探测器采集到线对卡的图像；对所述线对卡的图像进行处理，得到所述线对卡的调制传递函数曲线；根据所述线对卡的调制传递函数曲线，得到所述成像平板探测器的线对比度值。本发明专利技术提供的利用线对卡自动计算成像平板探测器线对比度的方法，通过计算调制传递函数MTF值，得到MTF曲线，通过设置MTF的阈值为0.4，取不低于MTF阈值的最大线对比度值，作为最终的成像平板探测器的线对比度值；本发明专利技术能够更加精确地计算出平板探测器的线对比度，而且提高对成像平板探测器线对比度的检验效率。

【技术实现步骤摘要】
一种利用线对卡自动计算成像平板探测器线对比度的方法
本专利技术涉及图像处理
，更具体地，涉及一种利用线对卡自动计算成像平板探测器线对比度的方法。
技术介绍
通常测试平板探测器线对比度的方法是把线对卡贴在平板探测器上，然后利用X射线球管进行曝光，读取平板探测器采集的图像，然后人工去查看每个线对下的成像质量，当不能清晰地分辨三条线的时候，则取上一个线对值作为最终的线对比度值。通过上述方式，判定时间较长，会受到个人主观臆断的影响，不能准确的得到成像平板探测器的线对比度。
技术实现思路
本专利技术的目的在于克服现有技术中存在的不足，提供了一种利用线对卡自动计算成像平板探测器线对比度的方法，能够更加精确地计算出平板探测器的线对比度，而且提高对成像平板探测器线对比度的检验效率。作为本专利技术的第一个方面，提供一种利用线对卡自动计算成像平板探测器线对比度的方法，包括：通过成像平板探测器采集到线对卡的图像；对所述线对卡的图像进行处理，得到所述线对卡的调制传递函数曲线；根据所述线对卡的调制传递函数曲线，得到所述成像平板探测器的线对比度值。进一步地，所述对所述线对卡的图像进行处理，得到所述线对卡的调制传递函数曲线中，还包括：读取所述线对卡的原始灰度图像；对所述线对卡的原始灰度图像进行图像二值化处理，并去掉所述原始灰度图像边缘处灰度值为0的黑色边框，得到所述线对卡的二值化图像；对所述线对卡的二值化图像进行边缘检测，得到所述线对卡的边缘图像，其中，所述边缘图像中黑色像素灰度值为0，白色像素灰度值为255；对所述线对卡的边缘图像进行孔洞填充，得到所述线对卡的孔洞图像；对所述线对卡的孔洞图像进行边缘检测，得到所述线对卡的轮廓图像；对所述线对卡的轮廓图像进行直线检测，得到所述线对卡的四条边所在的直线集合；对所述直线集合中的直线进行交点的计算，得到四条直线的交点集合；在所述四条直线的交点集合中，选择所述线对卡的轮廓图像上最左侧的交点和最下侧的交点，及通过计算两交点的向量方向，获得出所述线对卡的旋转角度；在所述线对卡的边缘图像上，利用得到的所述旋转角度，得到所述线对卡的调制传递函数曲线。进一步地，所述在所述线对卡的边缘图像上，利用得到的所述旋转角度，得到所述线对卡的调制传递函数曲线中，还包括：在所述线对卡的边缘图像上，利用得到的所述旋转角度，计算出所述线对卡上每个线对区的矩形感兴趣区域位置，其中，所述线对卡上共有20个线对区，每个线对区相对所述线对卡的位置是不变的；在每个线对区的矩形感兴趣区域中，依次按列进行扫描，计算所述线对卡的边缘图像对应列中，灰度值为255的最小行值点和最大行值点；根据计算出的每一列的最小行值点和最大行值点，在所述线对卡的原始灰度图像中，统计出每个线对区的矩形感兴趣区域中每一列的最小行值点和最大行值点之间的所有像素灰度值；分别计算出20个线对区的矩形感兴趣区域内的所有像素灰度值中的最大值和最小值；选择线对值为0.6LP/mm的线对区的矩形感兴趣区域内的最大像素灰度值MaxValue和最小像素灰度值MinValue，设置所述最大像素灰度值MaxValue与最小像素灰度值MinValue之间的差值作为实际对比度Δμ；根据所述实际对比度，计算其他的每个线对区相对0.6LP/mm线对区的调制传递函数MTF值；根据计算出的其他每个线对区的调制传递函数MTF值以及对应的线对值，得到所述线对卡的调制传递函数曲线。进一步地，所述实际对比度Δμ的计算公式如下：Δμ＝MaxValue–MinValue其中，MaxValue代表最大像素灰度值，MinValue代表最小像素灰度值。进一步地，所述调制传递函数MTF值的计算公式如下：MTF＝(MaxValue(n)–MinValue(n))/Δμ＝Δμn/Δμ其中，n为2,3……20，设定0.6LP/mm线对区作为第1个线对区；其中，当n为2时，计算第2个线对区相对于0.6LP/mm线对区的调制传递函数MTF值，即将第2个线对区的矩形感兴趣区域内的最大像素灰度值MaxValue(2)与最小像素灰度值MinValue(2)之间的差值除以所述实际对比度Δμ后得到的值作为该第2个线对区的MTF值；以此类推，直至计算出第20个线对区相对于0.6LP/mm线对区的调制传递函数MTF值。进一步地，所述根据所述线对卡的调制传递函数曲线，得到所述成像平板探测器的线对比度值中，还包括：设置一调制度阈值，根据所述线对卡的调制传递函数曲线，得到满足所述调制度阈值的所述成像平板探测器的线对比度值。进一步地，所述调制度阈值为0.4。进一步地，所述通过成像平板探测器采集到线对卡的图像中，还包括：将所述线对卡放置并固定在所述成像平板探测器上；通过X射线球管正对所述成像平板探测器进行曝光，所述成像平板探测器采集到线对卡的原始灰度图像。本专利技术提供的利用线对卡自动计算成像平板探测器线对比度的方法具有以下优点：采用图像处理技术，通过计算调制传递函数(MTF)值，得到MTF曲线，通过设置MTF的阈值为0.4，取不低于MTF阈值的最大线对比度值，作为最终的成像平板探测器的线对比度值；本专利技术能够更加精确地计算出平板探测器的线对比度，而且提高对成像平板探测器线对比度的检验效率。附图说明附图是用来提供对本专利技术的进一步理解，并且构成说明书的一部分，与下面的具体实施方式一起用于解释本专利技术，但并不构成对本专利技术的限制。图1为本专利技术提供的利用线对卡自动计算成像平板探测器线对比度的方法的流程图。图2为本专利技术提供的利用线对卡自动计算成像平板探测器线对比度的方法的具体实施方式流程图。图3为本专利技术提供的线对卡的结构示意图。图4为本专利技术提供的线对卡的原始灰度图像示意图。图5为本专利技术提供的线对卡的二值化图像示意图。图6为本专利技术提供的线对卡的原始边缘图像示意图。图7为本专利技术提供的线对卡的边缘图像局部放大示意图。图8为本专利技术提供的线对卡的孔洞图像示意图。图9为本专利技术提供的线对卡的轮廓图像直线检测图。图10为本专利技术提供的线对卡的轮廓图像直线检测局部放大图。图11为本专利技术提供的线对卡的四条直线的交点集合示意图。图12为本专利技术提供的线对卡上每个线对区的感兴趣区域位置示意图。图13为本专利技术提供的线对卡的边缘图像对应列中，灰度值为255的最小行值点和最大行值点示意图。图14为本专利技术提供的线对卡的MTF曲线示意图。图15为本专利技术提供的采集线对卡图像的装置结构图。具体实施方式为更进一步阐述本专利技术为达成预定专利技术目的所采取的技术手段及功效,以下结合附图及较佳实施例，对依据本专利技术提出的利用线对卡自动计算成像平板探测器线对比度的方法其具体实施方式、结构、特征及其本文档来自技高网...

【技术保护点】
1.一种利用线对卡自动计算成像平板探测器线对比度的方法，其特征在于，包括：/n通过成像平板探测器采集到线对卡的图像；/n对所述线对卡的图像进行处理，得到所述线对卡的调制传递函数曲线；/n根据所述线对卡的调制传递函数曲线，得到所述成像平板探测器的线对比度值。/n

【技术特征摘要】
1.一种利用线对卡自动计算成像平板探测器线对比度的方法，其特征在于，包括：
通过成像平板探测器采集到线对卡的图像；
对所述线对卡的图像进行处理，得到所述线对卡的调制传递函数曲线；
根据所述线对卡的调制传递函数曲线，得到所述成像平板探测器的线对比度值。


2.根据权利要求1所述的利用线对卡自动计算成像平板探测器线对比度的方法，其特征在于，所述对所述线对卡的图像进行处理，得到所述线对卡的调制传递函数曲线中，还包括：
读取所述线对卡的原始灰度图像；
对所述线对卡的原始灰度图像进行图像二值化处理，并去掉所述原始灰度图像边缘处灰度值为0的黑色边框，得到所述线对卡的二值化图像；
对所述线对卡的二值化图像进行边缘检测，得到所述线对卡的边缘图像，其中，所述边缘图像中黑色像素灰度值为0，白色像素灰度值为255；
对所述线对卡的边缘图像进行孔洞填充，得到所述线对卡的孔洞图像；
对所述线对卡的孔洞图像进行边缘检测，得到所述线对卡的轮廓图像；
对所述线对卡的轮廓图像进行直线检测，得到所述线对卡的四条边所在的直线集合；
对所述直线集合中的直线进行交点的计算，得到四条直线的交点集合；
在所述四条直线的交点集合中，选择所述线对卡的轮廓图像上最左侧的交点和最下侧的交点，及通过计算两交点的向量方向，获得出所述线对卡的旋转角度；
在所述线对卡的边缘图像上，利用得到的所述旋转角度，得到所述线对卡的调制传递函数曲线。


3.根据权利要求2所述的利用线对卡自动计算成像平板探测器线对比度的方法，其特征在于，所述在所述线对卡的边缘图像上，利用得到的所述旋转角度，得到所述线对卡的调制传递函数曲线中，还包括：
在所述线对卡的边缘图像上，利用得到的所述旋转角度，计算出所述线对卡上每个线对区的矩形感兴趣区域位置，其中，所述线对卡上共有20个线对区，每个线对区相对所述线对卡的位置是不变的；
在每个线对区的矩形感兴趣区域中，依次按列进行扫描，计算所述线对卡的边缘图像对应列中，灰度值为255的最小行值点和最大行值点；
根据计算出的每一列的最小行值点和最大行值点，在所述线对卡的原始灰度图像中，统计出每个线对区的矩形感兴趣区域中每一列的最小行值点和最大行值点之间的所有像素灰度值；
分别计算出20个线对区的矩形感兴趣区域内的所有像素灰度值中的最大值和最小值；
选择线对值为0.6LP/mm的...

【专利技术属性】
技术研发人员：王嘉军，梅建辉，马善达，凡小会，
申请(专利权)人：江苏瑞尔医疗科技有限公司，
类型：发明
国别省市：江苏;32