专利名称：一种检测图像中条状伪影的系统和方法
技术领域：
本发明涉及医疗成像领域，尤其涉及一种检测图像中条状伪影的系统和方法。
背景技术：
在X射线计算机断层成像(X-ray Computed Tomography,CT)的图像重建中，当待检对象对X射线的衰减系数呈非连续性的跳跃变化，如一种物质的密度远远大于另一种物质的密度时，CT图像中高密度物质的周围就会产生条状伪影。条状伪影在CT成像中是一种常见的伪影，其影响成像结果，使图像不能正确反映组织的病变特性，干扰医生对病变组织的客观判断，容易发生误诊或漏诊。因此，同提高信噪比和分辨率一样，识别和消除条状伪影也是提高图像质量的重要环节。目前，大多数与条状伪影相关的文章和专利都关注在用统计学的方法来减少条状伪影°如 Patrick J. La Rivigre 等人的 Reduction of noise-induced streak artifacts in χ-ray CT throughspline-based penalized-likelihood sinogram smoothing(通过基于样条的惩罚似然正弦图平滑滤波以降低χ射线计算机断层扫描设备中由噪音引起的条状伪影)(Conf Record IEEENSS-MIC，2003)，他们在文章中提出一种统计学上的原则性的正弦图平滑方法，通过使用惩罚泊松似然函数以减少条状伪影引起的噪声，其中，惩罚泊松似然函数在CT中的应用已经很成熟。但是在该文章中却没有指出条状伪影的位置。这样就无法检测到设备在出厂前由于硬件故障(如探测器通道损坏)所引起的条状伪影；同时这些条状伪影都是通过人眼进行判断的，存在人为误差。
发明内容
有鉴于此，本发明提出一种检测图像中条状伪影的系统和方法，以检测条状伪影及其准确位置，并以此作为反馈来评估CT扫描仪的运行状况。本发明提出一种检测图像中条状伪影的系统，包括梯度图像获取模块、感兴趣区选择模块、雷登变换模块、比较模块和条状伪影标记模块，其中梯度图像获取模块，用于获取所述图像在水平方向和竖直方向的两幅梯度图像，并将所述梯度图像传送给感兴趣区选择模块；感兴趣区选择模块，用于在所述梯度图像中选择一感兴趣区，并将所述感兴趣区传送给雷登变换模块；雷登变换模块，用于计算所述感兴趣区的雷登变换，并将雷登变换结果传送给比较模块；比较模块，用于将所述雷登变换结果与门限值进行比较，并将比较结果传送给条状伪影标记模块；条状伪影标记模块，用于根据所述比较结果在所述感兴趣区中标记出条状伪影。所述感兴趣区选择模块包括区域选择单元和比较单元，其中区域选择单元，用于在所述水平方向和竖直方向的两幅梯度图像中分别选取一个区域，并将这两个区域传送给比较单元；比较单元，用于比较所述两个区域中的最大灰度梯度值，将较大者所在的区域作为所述感兴趣区，并将所述感兴趣区传送给雷登变换模块。所述比较模块包括角度计算单元、角度门限值比较单元和概率门限值比较单元，其中角度计算单元，用于获取所述雷登变换结果中最大的η个值在所述感兴趣区中所对应的η条直线的角度，计算这η条直线的平均角度，其中η是自然数，并将所述η条直线的角度和所述平均角度传送给角度门限值比较单元；角度门限值比较单元，用于计算所述η条直线的角度分别与所述平均角度之差的绝对值，得到η个绝对值，将这η个绝对值与一角度门限值进行比较，得到小于该角度门限值的s条直线，其中s为不大于η的自然数，并将s 与η传送给概率门限值比较单元；概率门限值比较单元，计算s与11的比值，将该比值与一概率门限值进行比较，并将比较结果传送给条状伪影标记模块。所述条状伪影标记模块进一步用于，根据所述比较结果得到所述条状伪影的准确位置。本发明还提出一种检测图像中条状伪影的方法，包括如下步骤获取所述图像在水平方向和竖直方向的两幅梯度图像；在所述梯度图像中选择一感兴趣区；计算所述感兴趣区的雷登变换；将所述雷登变换的结果与门限值进行比较；根据所述比较的结果在所述感兴趣区中标记出条状伪影。所述选择感兴趣区包括分别选择所述两幅梯度图像的一个区域，取这两个区域中的最大灰度梯度值进行比较，将较大者所在的区域作为所述感兴趣区。所述比较包括获取所述雷登变换结果中最大的η个值在所述感兴趣区中所对应的η条直线的角度，并计算这η条直线的平均角度，其中η是自然数；计算所述η条直线的角度分别与所述平均角度之差的绝对值，得到η个绝对值，将这η个绝对值与一角度门限值进行比较，得到小于该角度门限值的s条直线，其中s为不大于η的自然数；计算s与11的比值，将该比值与一概率门限值进行比较，若所述比值大于所述概率门限值，则这s条直线为条状伪影。优选地，所述角度门限值小于等于30度。优选地，所述角度门限值小于等于5度。优选地，所述概率门限值小于等于1。进一步，根据所述条状伪影得到所述条状伪影在所述图像中的准确位置。从上述方案中可以看出，由于本发明能得到条状伪影所在直线的角度，所以能在原始图像中标记出条状伪影及其准确位置，进一步以此作为反馈来评估CT扫描仪或X射线管的运行状况。另一方面，由于传统的条状伪影检测方法是通过人的肉眼进行判断，存在人为误差，而本发明则是完全自动化的，避免了人为判断造成的误差，提高了判断条状伪影的准确率。
图1是原始的CT图像；图2是图1的梯度图像的示意图，其中图加是图1在X方向的梯度图像FX，图2b 是图1在Y方向的梯度图像FY的示意图；图3是图2的梯度图像的一个方形区域的示意图，其中图3a是图加的一个方形区域crop_FX的示意图，图北是图2b的一个方形区域crop_FY的示意图；图4是感兴趣区的Radon变换结果R的示意图；图5是图像中条状伪影的检测结果的示意图6是本发明检测图像中条状伪影系统的组成示意图；图7是本发明感兴趣区选择模块的组成示意图；图8是本发明比较模块的组成示意图；图9是本发明检测图像中条状伪影方法的流程图。
具体实施例方式为使本发明的目的、技术方案和优点更加清楚，以下参照实施例对本发明作进一步详细说明。本发明中水平向右方向为X方向的正方向，竖直向下方向为Y方向的正方向。本发明提出一种检测图像中条状伪影的系统，如图6所示，是本发明检测图像中条状伪影系统的组成示意图，包括梯度图像获取模块10、感兴趣区选择模块20、雷登变换模块30、比较模块40和条状伪影标记模块50，其中梯度图像获取模块10，用于获取所述图像在水平方向和竖直方向的两幅梯度图像，并将所述梯度图像传送给感兴趣区选择模块20。图1是原始的CT图像。按下式(1)来计算图1所示图像分别在X方向和Y方向上的灰度梯度值。 FX (X，y) = (I (x+1，y) -I (x_l，y)) /2FY (x, y) = (I (χ, y+1)-I (χ, y-1))/2 (1)其中(χ，y)是图像中的任一点(即像素)，I (X，y)是在(x，y)处像素的灰度值， FX(X，y)和FY(X，y)分别是像素(x，y)在X方向和Y方向上的灰度梯度值(x，y)。灰度值的范围从0到255不等。对于所述图像中的边界像素，按下式( 来计算灰度梯度值。FX (X，y) = (I (x+1，y) -I (x, y)) /2FY (χ, y) = (I (χ, y+1) -I (χ, y)) /2FX (x,y) = (I (χ, y) -I (χ-l，y)) /2FY(x,y) = (I (χ, y)-I (χ, y-1))/2 (2)通过式(1)- )得到如图2所示的两幅梯度图像，其中图加是图1所示图像在X 方向上的梯度图像FX，图2b是图1所示图像在Y方向上的梯度图像FY，这两幅图的横坐标为像素在X方向的位置，纵坐标为像素在Y方向的位置。感兴趣区选择模块20，用于根据所述梯度图像来选择一感兴趣区(Region of interest, R0I)P,以下用R0I_P表示，并将所述感兴趣区R0I_P传送给雷登变换模块30。雷登变换模块30，用于计算所述感兴趣区的雷登变换，并将雷登变换结果传送给比较模块40。雷登变换(即Radon变换)的原理是将原始图象(如XY平面内的二维图像)内的每一点通过线积分的形式变换到另一对参数域(如投影值，角度)内，若原始图像内存在直线，则会导致Radon变化的结果在该直线对应的另一对参数域内的函数值出现极值，所以可用Radon变换来检测原始图像中是否含有直线成分。在本发明的实施例中，进行Radon变换的做法如下将感兴趣区R0I_P分别沿着0 度到179度的投影方向，对每一条投影线计算R0I_P的线积分，得到如图4所示的Radon变换结果R，其中投影方向每次旋转1度，横坐标为旋转角度，纵坐标为R0I_P矩阵的对角线长度。在本发明的实施例中，由于设置ROI是90*90的方形矩阵，所以图4中纵坐标的最大值为90。图4中线条的明暗表示灰度梯度值的累积值的大小，越明亮的区域表明累积值越大，则该区域可能是极值所在的区域，也就说明该角度附近可能存在条状伪影，图中可见横坐标在150度附近较其他角度要亮很多，这反映了图像原始CT图像在150度附近可能存在条状伪影。比较模块40，用于将所述雷登变换结果R与门限值进行比较，并将比较结果传送给条状伪影标记模块50。条状伪影标记模块50，用于根据所述比较结果在R0I_P中标记出条状伪影。图7是本发明感兴趣区选择模块20的组成示意图。图7中感兴趣区选择模块20 包括区域选择单元201和比较单元202，其中区域选择单元201，用于在所述水平方向和竖直方向的两幅梯度图像中分别选取一个区域，并将这两个区域传送给比较单元202。通常在图2中将梯度图像分成若干个方形区域(或称为方形矩阵)，每次按顺序来检测一个方形区域(如90*90)中是否存在条状伪影。本发明中选择图2中的两幅梯度图像FX和FY的一个方形区域，得到如图3所示的示意图，其中图3a是图加的一个方形区域 crop_FX,图3b是图2b的一个方形区域crop_FX。通常crop_FX与crop_FX在梯度图像中的位置一致，以便确定该方形区域在原始图像中是在X方向还是Y方向的灰度梯度值变化较大。比较单元202，用于比较所述两个区域中的最大灰度梯度值，将较大者所在的区域作为所述感兴趣区，并将所述感兴趣区传送给雷登变换模块。这里，比较单元202按下式(3)来计算crop_FX的最大灰度梯度值max_FX和crop_ FY的最大灰度梯度值max_FY。max_FX = max(crop_FX)max_FY = max (crop_FY)(3)其中，max是找出矩阵中最大值的函数。之后，比较max_FX和max_FY的大小，将二者中的较大者所在的方形区域作为所述感兴趣区R0I_P。例如，若max_FX > max_FY,则R0I_P对应crop_FX。若二者相等，则R0I_ P可以是这二者中的任一者。通常认为，灰度梯度值变化大的区域存在条状伪影的可能性大。图8是本发明比较模块40的组成示意图。图8中比较模块40包括角度计算单元 401、角度门限值比较单元402和概率门限值比较单元403，其中角度计算单元401，用于获取所述雷登变换结果中最大的η个值在所述感兴趣区中所对应的η条直线的角度，计算这η条直线的平均角度，其中η是自然数，并将所述η条直线的角度和所述平均角度传送给角度门限值比较单元402 ；这里，角度计算单元401要获取Radon变换结果R中亮度最高的η个值，其与R0I_ P中的η条直线对应，接着计算这η条直线的角度(或斜率)gl，g2，...，gn。其中η可依据需要来设置，如η为10。角度门限值比较单元402，用于计算所述η条直线的角度分别与所述平均角度之差的绝对值，得到η个绝对值，将这η个绝对值与一角度门限值进行比较，得到小于该角度门限值的s条直线，其中s为不大于η的自然数，并将s与η传送给概率门限值比较单元 403 ；首先按下式(4)来计算gl，g2，· · ·，gn的平均角度g。接着按下式(5)来分别计算角度gl，g2,...,gn与平均角度g之间的差的绝对值
其中abs是计算绝对值的函数。最后得到al，a2，. .，an中小于预设角度门限t_g的个数S。其中s为不大于η的自然数，t_g的值依经验设定。优选地，可以设定角度门限t_ g为5度，则绝对值al，a2，. . .，an中小于5度的可认为这些直线具有非常接近的方向，而大于5度的则认为是肌1_ 图像中杂乱无章的直线。当然也可以设定角度门限值其他数值， 如小于等于30度或更优选地为小于5度。概率门限值比较单元403，计算s与η的比值r，将该比值与一概率门限值进行比较，并将比较结果传送给条状伪影标记模块。首先按下式(6)来计算比值r。
S” = - (6)
η其中r为斜率较接近的直线个数s占总的直线个数η的比例。接着将r与概率门限值t_r进行比较。其中t_r的值依经验设定，数值范围在(0，1]。若！· > t_r，则说明R0I_P中存在条状伪影，否则，就说明R0I_P中没有条状伪影。若r> t_r,则条状伪影标记模块50在R0I_ P中标记出这s条条状伪影。进一步，条状伪影标记模块50还用于，根据比较模块40的比较结果得到所述条状伪影的准确位置。若概率门限值比较单元403的比较结果为r > t_r，则说明R0I_P中存在条状伪影，也就是说角度门限值比较单元402中的s条直线都是条状伪影。由于原始图像和梯度图像的大小一致，且各个点的位置坐标也一致，所以可以直接将这些条状伪影在原始图像中标记出来。这s条直线斜率分布在平均角度g附近，参照图4，此平均角度g在150度附近。同时由于R0I_P是crop_FX或crop_FY，实际上也是一梯度图像，R0I_P中s条直线的斜率和位置坐标与这s条直线在原始CT图像中的斜率和位置坐标一致，所以能根据R0I_P 中的条状伪影得到这些条状伪影在原始CT图像中的具体位置。在检测出R0I_P在图1所示的原始图像中对应的条状伪影后，继续在图2中选取另一方形区域进行检测，直到检测出原始CT图像中所有的条状伪影及其准确位置，得到如图5所示的结果，图中的白线即是检测到的条状伪影。由于Y方向的正方向是竖直向下，由此这些条状伪影的斜率大概是150度。同时由于原始CT图像是各个方形区域存在的条状伪影的角度不同，同一方形矩阵中条状伪影的角度也只是接近，而不是完全一致，所以在图 4可见条状伪影的角度不尽相同。
8
本发明还提供一种检测图像中条状伪影的方法，如图9所示，包括如下步骤步骤101，获取所述图像在水平方向和竖直方向的两幅梯度图像。可按式(1)-( 来计算图1所示原始图像分别在X方向和Y方向上的两幅梯度图像。步骤102，分别选择所述两幅梯度图像的一个方形区域crop_FX和crop_FX。步骤103，取这两个区域中的最大灰度梯度值进行比较，将较大者所在的区域作为感兴趣区。按式(3)来计算crop_FX的最大灰度梯度值max_FX和crop_FY的最大灰度梯度值maX_FY，将较大者所在的方形区域作为感兴趣区R0I_P。步骤104，计算所述感兴趣区的雷登变换，得到如图4所示的Radon变换结果R。步骤105，获取所述雷登变换结果中最大的η个值在所述感兴趣区中所对应的η条直线的角度，并计算这η条直线的平均角度，其中η是自然数。按式(4)来计算这η条直线的平均角度。步骤106，计算所述η条直线的角度分别与所述平均角度之差的绝对值，得到η个绝对值，将这η个绝对值与一角度门限值进行比较，得到小于该角度门限值的s条直线，其中s为不大于η的自然数；优选地，所述角度门限值小于等于30度。更优选地，所述角度门限值小于等于5度。步骤107，计算s与η的比值，将该比值与一概率门限值进行比较，若所述比值大于所述概率门限值，则这s条直线为条状伪影。按式(6)来计算所述比值r。优选地，所述概率门限值小于等于1。步骤108，根据所述比较的结果在所述感兴趣区中标记出条状伪影。进一步，本发明的检测图像中条状伪影的方法还包括步骤109，根据所述条状伪影得到所述条状伪影在所述图像中的准确位置。这点在条状伪影标记模块50中已有具体说明，这里不再赘述。本发明的检测图像中条状伪影的系统和方法能用来检查出厂前的医学成像系统中扫描仪或X射线管的运行状况。若图像中存在条状伪影，则说明扫描仪的探测器通道损坏或X射线管的球管打火，需要检查以对损坏的部件进行维修或更换。这样就可以在出厂前检测到由于硬件损坏所引起的条状伪影。另外，本发明中检测条状伪影的过程是完全自动化的，不需要肉眼进行人为判断，避免了人为误差。理论上，本发明可以用于所有的医学成像系统，包括CT，磁共振成像系统和X射线摄影系统等。以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
权利要求
1.一种检测图像中条状伪影的系统，包括梯度图像获取模块、感兴趣区选择模块、雷登变换模块、比较模块和条状伪影标记模块，其中梯度图像获取模块，用于获取所述图像在水平方向和竖直方向的两幅梯度图像，并将所述梯度图像传送给感兴趣区选择模块；感兴趣区选择模块，用于在所述梯度图像中选择一感兴趣区，并将所述感兴趣区传送给雷登变换模块；雷登变换模块，用于计算所述感兴趣区的雷登变换，并将雷登变换结果传送给比较模块；比较模块，用于将所述雷登变换结果与门限值进行比较，并将比较结果传送给条状伪影标记模块；条状伪影标记模块，用于根据所述比较结果在所述感兴趣区中标记出条状伪影。
2.根据权利要求1所述的系统，其特征在于，所述感兴趣区选择模块包括区域选择单元和比较单元，其中区域选择单元，用于在所述水平方向和竖直方向的两幅梯度图像中分别选取一个区域，并将这两个区域传送给比较单元；比较单元，用于比较所述两个区域中的最大灰度梯度值，将较大者所在的区域作为所述感兴趣区，并将所述感兴趣区传送给雷登变换模块。
3.根据权利要求1所述的系统，其特征在于，所述比较模块包括角度计算单元、角度门限值比较单元和概率门限值比较单元，其中角度计算单元，用于获取所述雷登变换结果中最大的η个值在所述感兴趣区中所对应的η条直线的角度，计算这η条直线的平均角度，其中η是自然数，并将所述η条直线的角度和所述平均角度传送给角度门限值比较单元；角度门限值比较单元，用于计算所述η条直线的角度分别与所述平均角度之差的绝对值，得到η个绝对值，将这η个绝对值与一角度门限值进行比较，得到小于该角度门限值的 s条直线，其中s为不大于η的自然数，并将s与η传送给概率门限值比较单元；概率门限值比较单元，计算s与η的比值，将该比值与一概率门限值进行比较，并将比较结果传送给条状伪影标记模块。
4.根据权利要求1所述的系统，其特征在于，所述条状伪影标记模块进一步用于，根据所述比较结果得到所述条状伪影的准确位置。
5.一种检测图像中条状伪影的方法，包括如下步骤 获取所述图像在水平方向和竖直方向的两幅梯度图像； 在所述梯度图像中选择一感兴趣区；计算所述感兴趣区的雷登变换；将所述雷登变换的结果与门限值进行比较；根据所述比较的结果在所述感兴趣区中标记出条状伪影。
6.根据权利要求5所述的方法，其特征在于，所述选择感兴趣区包括分别选择所述两幅梯度图像的一个区域，取这两个区域中的最大灰度梯度值进行比较，将较大者所在的区域作为所述感兴趣区。
7.根据权利要求5所述的方法，其特征在于，所述比较包括获取所述雷登变换结果中最大的η个值在所述感兴趣区中所对应的η条直线的角度， 并计算这η条直线的平均角度，其中η是自然数；计算所述η条直线的角度分别与所述平均角度之差的绝对值，得到η个绝对值，将这η 个绝对值与一角度门限值进行比较，得到小于该角度门限值的s条直线，其中s为不大于η 的自然数；计算s与η的比值，将该比值与一概率门限值进行比较，若所述比值大于所述概率门限值，则这s条直线为条状伪影。
8.根据权利要求7所述的方法，其特征在于，所述角度门限值小于等于30度。
9.根据权利要求7所述的方法，其特征在于，所述角度门限值小于等于5度。
10.根据权利要求7所述的方法，其特征在于，所述概率门限值小于等于1。
11.根据权利要求5所述的方法，其特征在于，进一步，根据所述条状伪影得到所述条状伪影在所述图像中的准确位置。
全文摘要
本发明公开一种检测图像中条状伪影的系统，包括梯度图像获取模块用于获取所述图像在水平方向和竖直方向的两幅梯度图像，并将所述梯度图像传送给感兴趣区选择模块；感兴趣区选择模块用于在所述梯度图像中选择一感兴趣区，并将所述感兴趣区传送给雷登变换模块；雷登变换模块用于计算所述感兴趣区的雷登变换，并将雷登变换结果传送给比较模块；比较模块用于将所述雷登变换结果与门限值进行比较，并将比较结果传送给条状伪影标记模块；条状伪影标记模块用于根据所述比较结果在所述感兴趣区中标记出条状伪影。本发明还提供一种检测图像中条状伪影的方法。采用本发明的系统和方法，能检测条状伪影及其准确位置，并评估扫描仪的各通道的运行状况。
文档编号G06T5/00GK102402787SQ201010286368
公开日2012年4月4日 申请日期2010年9月19日 优先权日2010年9月19日
发明者杜春华 申请人:上海西门子医疗器械有限公司