数字X线成像系统的二维调制传递函数测量方法技术方案

本发明专利技术属于医学工程技术领域，尤其涉及一种数字X线成像系统的二维调制传递函数测量方法和系统。本发明专利技术数字X线成像系统的二维调制传递函数测量方法，包括获取不同倾斜角度的刃边图像，计算其相对应位置的调制传递函数，然后通过插值拟合获得数字X线成像系统二维调制传递函数的曲面图。本发明专利技术的数字X线成像系统的二维调制传递函数测量方法和系统摆脱了现有技术过度依赖探测器对称性的缺点；通过引入角度测量仪获得刃边的位置，与通过线性拟合和霍夫变换等通过图像处理算法得到的刃边位置相比，本发明专利技术所用的方法具有更高的精度，从而降低了由于刃边位置的不精确性带来的噪声。

【技术实现步骤摘要】
数字X线成像系统的二维调制传递函数测量方法
本专利技术属于医学工程
，尤其涉及数字X线成像系统的二维调制传递函数测量方法。
技术介绍
随着X线成像系统的发展，成像质量的评价也在逐渐发展，数字X线成像系统的调制传递函数(modulationtransferfunction,MTF)不仅是反映其空间分辨力的一个重要参数，而且是获得系统的量子检出效率(detectivequantumefficiency,DQE)的必要条件。调制传递函数是用空间频率表示输出调制和输入调制之比的函数，可对成像系统的分辨力参数进行定量描述，并可以精确地描述成像系统和其组成部分的信息再现率。MTF是衡量成像系统性能的重要指标参数，作为空间频率因素，调制传递函数用来测量分辨率在整个成像系统中的传递情况。现有测量MTF的方法中主要有测试卡法、狭缝法和刃边法，但是这些方法只能获得系统在某一个方向上的MTF，即只能获得水平或垂直方向的一维MTF，无法对成像系统在整个平面的分辨力传递特性进行描述，具有很大的局限性。因此若要获得系统的完整描述，必须计算成像系统的二维调制传递函数。现有的数字X线成像系统的二维调制传递函数测量方法包括：现有技术一：在文献<A.Kuhls-Gilcrist,A.Jain,D.R.Bednarek,K.R.Hoffmann,andS.Rudin,“AccurateMTFmeasurementindigitalradiographyusingnoiseresponse,”Med.Phys.37,724-735(2010).>中根据噪声功率谱和调制传递函数之间的关系NPS(u,v)＝ηMTFdet(u,v)+NPSelectronic(u,v)，通过测量成像系统的噪声功率谱，和特定空间频率下的调制传递函数值计算两者之间的相关系数，然后获得系统的调制传递函数，但是现有技术一的缺点在于：要获得二维MTF的函数形式，必须进行深入的线性级联系统分析；现有技术二：在文献<K.A.Fetterly,N.J.Hangiandreou,B.A.Schueler,andE.R.Ritenour,“Measurementofthepresampledtwo-dimensionalmodulationtransferfunctionofdigitalimagingsystems,”Med.Phys.29,913-921(2002)>中，利用精加工的针孔设备来产生X线的近似点光源，获得成像系统的二维点扩散函数，然后通过傅里叶变换得到系统的二维MTF；但是现有技术二的缺点在于：若要获取高精度的MTF，需要处理大量的图像数据，另外，测量装置精度要求太高，不适合大面积推广应用现有技术三：在文献<M.Bath，P.Sund，L.G.Mansson.“Methodfordeterminingthetwo-dimensionalpre-samplingmodulationtransferfunctionindigitalradiography,”Proc.SPIEVOL.4320>中通过对上面构造有N×N个圆孔的不透光圆盘进行成像，得到圆盘扩散函数，然后经过重新采样，傅里叶变换及校正等计算步骤得到成像系统的二维MTF；但是现有技术三的缺点在于：由于高频混叠现象使整体的测量精度尤其是高频精度大大降低。
技术实现思路
本专利技术提供了一种数字X线成像系统的二维调制传递函数测量方法，旨在解决现有的数字X线成像系统的二维调制传递函数测量方法对探测器对称性依赖高，且测量精度尤其是高频精度低的技术问题。本专利技术提供的技术方案为：一种数字X线成像系统的二维调制传递函数测量方法，包括：步骤a：将刃边装置放置在探测器表面，调整刃边与探测器采样方向成初始角度θ0，设置曝光参数，在一定辐射质量下对对刃边装置进行多次曝光；步骤b：旋转刃边装置，调整刃边装置使刃边与探测器的角度均匀变化，对刃边装置进行曝光，每个位置曝光多次，获取多幅图像；步骤c：对图像进行线性化处理，将二维调制传递函数的分析区域内的图像数据沿着刃边的方向进行投影得到亚像素数组，通过投影变换得到该方向的边缘响应函数；步骤d：对得到的边缘响应函数ERF降噪处理，对ERF曲线进行拟合，处理后的ERF进行差分运算得到线扩散函数LSF(xi)；步骤e：对LSF(xi)进行傅里叶变换并取模，利用零频率处MTF值进行归一化处理，得到刃边位置为θi时系统的归一化调制传递函数MTFi(u,v)；步骤f：对得到的所有方向的一维调制传递函数MTFk(x)通过双线性差值进行曲面拟合得到系统的二维调制传递函数MTF(u,v)。本专利技术的技术方案还包括：所述步骤a前进一步包括：构建刃边装置，所述刃边装置包括环形导轨、角度测量仪、铅板和钨板，且钨板的边缘经过抛光处理，所述铅板固定在环形导轨上，所述钨板固定在铅板中央切割出的矩形区域。本专利技术的技术方案还包括：在所述步骤b中，调整刃边装置使刃边与探测器的角度在0和180之间均匀变化，刃边与探测器采样方向的夹角为θk，且满足(θkmod45)≠0对刃边装置进行曝光。本专利技术的技术方案还包括：在所述步骤c中，所述对图像进行线性化处理步骤后还包括：通过线性拟合得到刃边与探测器采样方向的夹角θk'，比较θk和θk'的大小关系。本专利技术的技术方案还包括：在所述步骤c中，所述通过投影变换得到该方向的边缘响应函数步骤包括：将二维调制传递函数的分析区域中的所有像素点(i,j)投影到s轴为：s(i,j)＝ipcosθk-jpsinθk，设采样间隔为Δs＝0.1p，采样公式为：其中：p为像素大小，i、j分别表示ROI区域的第i行和第j列，s(i,j)为像素(i,j)到刃边的距离，nm为到刃边的距离在(m-0.5)Δs和(m+0.5)Δs之间的像素个数；对落在同一间隔内的所有数据点的几何平均作为该间距的数据，得到系统的过采样边缘响应函数ERF。本专利技术的技术方案还包括：在所述步骤d中，所述对ERF曲线进行拟合步骤具体为：通过改进的费米函数对ERF曲线进行拟合，拟合公式为：本专利技术的技术方案还包括：所述步骤e包括：方向θk的边缘响应函数经过差分计算得到该方向的线扩散函数：对其做傅里叶变换得到方向θk的光学传递函数，对其取模得到调制传递函数MTFk'(x)；根据零频率处的MTF值进行归一化，得到归一化的调制传递函数MTFk(x)。本专利技术的技术方案还包括：所述步骤f包括：将得到的所有方向的一维调制传递函数MTFk(x)公式转换为二维平面中的MTF值：通过双线性差值，经过曲面拟合得到数字X线成像系统的曲面图。本专利技术的技术方案具有如下优点或有益效果：本专利技术数字X线成像系统的二维调制传递函数测量方法获取不同倾斜角度的刃边图像，计算其相对应位置的调制传递函数，然后通过插值拟合获得数字X线成像系统二维调制传递函数的曲面图，本专利技术可本文档来自技高网...

【技术保护点】
一种数字X线成像系统的二维调制传递函数测量方法，包括：步骤a：将刃边装置放置在探测器表面，调整刃边与探测器采样方向成初始角度θ0，设置曝光参数，在一定辐射质量下对对刃边装置进行多次曝光；步骤b：旋转刃边装置，调整刃边装置使刃边与探测器的角度均匀变化，对刃边装置进行曝光，每个位置曝光多次，获取多幅图像；步骤c：对图像进行线性化处理，将二维调制传递函数的分析区域内的图像数据沿着刃边的方向进行投影得到亚像素数组，通过投影变换得到该方向的边缘响应函数；步骤d：对得到的边缘响应函数ERF降噪处理，对ERF曲线进行拟合，处理后的ERF进行差分运算得到线扩散函数LSF(xi)；步骤e：对LSF(xi)进行傅里叶变换并取模，利用零频率处MTF值进行归一化处理，得到刃边位置为θi时系统的归一化调制传递函数MTFi(u,v)；步骤f：对得到的所有方向的一维调制传递函数MTFk(x)通过双线性差值进行曲面拟合得到系统的二维调制传递函数MTF(u,v)。

【技术特征摘要】
1.一种数字X线成像系统的二维调制传递函数测量方法，包括：
步骤a：将刃边装置放置在探测器表面，调整刃边与探测器采样方向成初始角度θ0，设置曝光参数，在一定辐射质量下对刃边装置进行多次曝光；
步骤b：旋转刃边装置，调整刃边装置使刃边与探测器的角度均匀变化，对刃边装置进行曝光，每个位置曝光多次，获取多幅图像；
步骤c：对图像进行线性化处理，将二维调制传递函数的分析区域内的图像数据沿着刃边的方向进行投影得到亚像素数组，通过投影变换得到该方向的边缘响应函数；
步骤d：对得到的边缘响应函数ERF降噪处理，对ERF曲线进行拟合，处理后的ERF进行差分运算得到线扩散函数；
步骤e：对所述线扩散函数进行傅里叶变换并取模，利用零频率处MTF值进行归一化处理，得到刃边与探测器采样方向的夹角为θi时系统的归一化调制传递函数；
步骤f：对得到的所有方向的一维调制传递函数通过双线性差值进行曲面拟合得到系统的二维调制传递函数。


2.根据权利要求1所述的数字X线成像系统的二维调制传递函数测量方法，其特征在于，所述步骤a前进一步包括：构建刃边装置，所述刃边装置包括环形导轨、角度测量仪、铅板和钨板，且钨板的边缘经过抛光处理，所述铅板固定在环形导轨上，所述钨板固定在铅板中央切割出的矩形区域。
...

【专利技术属性】
技术研发人员：孟书先，熊璟，李志成，谢耀钦，刘勇，李生广，陈鸣闽，王丽艳，
申请(专利权)人：中国科学院深圳先进技术研究院，
类型：发明
国别省市：广东;44