专利名称：医用图像处理装置的制作方法
技术领域：
本发明涉及一种对利用MRI (Magnetic Resonance Imaging:磁共振成像)等得到的脑图像进行图像处理后进行脑疾病的诊断支援的技术，尤其涉及一种将利用MRI等得到的脑图像处理成适于诊断支援的状态的技术。
背景技术：
由于高龄化社会的到来，认知症疾病的患者逐年增加。认知症疾病有各种种类，必需在诊断中加以区别，以实施对应于疾病的适当处置。另一方面，为了应对这种要求，近年来能通过SPECT(Single Photon EmissionComputed Tomography:单光子发射计算机断层成像术)或PET (Positron EmissionTomography:正电子发射型计算机断层显像)等核医学检查、或者CT(ComputerizedTomography:计算机断层扫描)或MRI来取得与脑的状态有关的信息。其结果是，可以了解脑的确定部位的血流或代谢低下或组织萎缩的现象因疾病不同而不同，需要与之对应的定量的评价方法。例如，脑的局部部位的血流或代谢低下能通过利用SPECT或PET图像进行比较来检定。另外，对于组织的萎缩，能利用MRI图像来求出确定部位的容积，比较其相对大小后判别有无异常。作为使用这种脑图像来评价脑的萎缩的方法，已知以作为3维像素的体素为单位图像处理对被检查者 的头部进行摄像而取得的脑图像来进行的VBM(Voxel BasedMorphometry:基于体素的形态测量学)(例如参照专利文献I)。该VBM方法是有效识别老年痴呆病的评价方法，对于健康者与老年痴呆病的识另U，报告称能诊断出87.8% (参照非专利文献I)。专利文献1:特开2005-237441号公报非专利文献1:Yoko Hirata, Hiroshi Matsuda, Kiyotaka Nemoto, TakashiOhnishi,Kentaro Hirao，Fumio Yamashita，Takashi Asada，Satoshi Iwabuchi,HirotsuguSamejima， Voxel-based morphometry to discriminate early Alzheimer' s diseasefrom controls.Neurosci Lett 382:269-274,2005非专利文献2:Bookstein FL “Voxel-based morphometry, shouId not be usedwith imperfectly registered images.Neuroimage.2001 ；14(6):1454-62.
非专利文献3:J.Ashburner and K.J.Friston.Unified segmentation.Neuroimage.2005 ；26:839-851.
非专利文献 4 ；Ashburner J，A fast diffeomorphic image registrationalgorithm.Neuroimage.2007 Oct 15 ；38 (I):95-113.
非专利文献5:松田博史:SPECT O統計学的画像解析(SPECT的统计学图像解析)。r ^ ^ 一型痴呆ο画像診断(老年痴呆型痴呆的图像诊断)，> '7'力 > 匕工一社:pp.76-86(2001).
但是，在基于上述现有方法的脑图像处理中，存在因作为处理对象的脑图像而不能顺利进行组织分离处理，从而无法正确确定组织的问题。
发明内容
因此，本发明的课题在于提供一种医用图像处理装置，该医用图像处理装置能判定取得的脑图像的状态，并调整为适于进行组织分离处理的脑图像。为了解决上述课题，在本发明第I方式中，提供一种医用图像处理装置，该医用图像处理装置具有:对象切片选择单元，其从由多个切片图像构成的脑图像中，选择作为处理对象的切片图像作为对象切片；脑实质测量单元，其进行脑实质的像素的信号值的有效最大值、即脑实质有效最大值的测量处理；脑图像测量单元，其进行脑图像整体的像素的信号值的有效最大值的测量处理；头盖测量单元，其测量处理成为头盖周边区域的峰值的像素的信号值的平均作为头盖周围的峰值平均值；抑制处理判定单元，其根据测量出的脑实质的有效最大值、脑图像整体的有效最大值、头盖周围的峰值平均值，判定是否需要高信号值抑制处理；以及高信号值抑制处理单元，其在由所述抑制处理判定单元判定为需要的情况下，对所述脑图像执行高信号值抑制处理。根据本发明第I方式，选择脑图像的对象切片，并对选择到的切片，测量脑实质的像素的信号值的有效最大值、脑图像整体的像素的信号值的有效最大值、成为头盖周边区域峰值的像素的信号值的平均即头盖周围的峰值平均值，根据测量结果，在必要的情况下，对脑图像执行高信号值抑制处理，因此能判定取得的脑图像的状态，并调整为适于组织分离处理的脑图像。另外，在本发明第2方式中，对于本发明第I方式的医用图像处理装置，其特征在于，所述抑制处理判定单元执行头盖区域判定以及图像整体判定，该头盖区域判定是头盖周围的峰值平均值是否比所述脑实质有效最大值高一定比例的判定，该图像整体判定是图像整体的信号值是否比所述脑实质有效最大值高一定比例的判定，在所述头盖区域判定、图像整体判定的任一方满足条件的情况下，判定为需要高信号值抑制处理。
根据本发明第2方式，进行头盖周边区域的峰值平均值、图像整体的信号值是否比脑实质有效最大值高一定比例的判定，在头盖周边区域的峰值平均值、图像整体的信号值的任一个比脑实质有效最大值高一定比例的情况下，判定为需要高信号值抑制处理，因此，即使在头盖周边区域的信号值比脑实质明显高的情况、头盖周边区域、脑实质以外部分中的信号值比脑实质明显高的情况的任一情况下，都能准确进行需要高信号值抑制处理的判定。另外，在本发明第3方式中，对于本发明第I或第2方式的医用图像处理装置，其特征在于，所述对象切片选择单元对各切片图像在图像内设定预定的线段，求出位于该线段上且信号值为预定值以上的像素中间隔最远的像素间的长度，将该像素间的长度与全部切片图像中最长的长度相比较，将成为预定比率以上的切片图像选为对象切片。根据本发明第3方式，在各切片图像中，将具有预定值以上信号值的像素与全部切片图像中最长的像素相比较，将成为预定比率以上的切片图像选为对象切片，因此，能选择适于判定是否需要高信号值抑制处理的、除头盖周边区域与脑实质以外不包含高信号值的切片图像。另外，在本发明第4方式中，对于本发明第I至第3方式的任一方式的医用图像处理装置，其特征在于，所述脑实质测量单元对选择的所述各对象切片设定预定数量的通过脑中心部的线段，按各线段的每一个，将线段上像素的信号值分布分割成预定数量的区域，将中央的预定数量的区域作为脑实质区域，求出该脑实质区域中的信号值的直方图，求出去除高值的预定像素后的信号最大值作为脑实质的有效最大值。根据本发明第4方式，对选择的各对象切片，将通过脑中心部的多个线段各自的信号值分布分割成规定数量的区域，将中央的预定数量的区域作为脑实质区域，求出脑实质区域中的信号值的直方图，求出去除高值的预定像素后的信号最大值，作为脑实质的有效最大值，因此，能准确求出作为代表脑实质的值的脑实质的有效最大值。另外，在本发明第5方式中，对于本发明第I至第4方式的任一方式的医用图像处理装置，其特征在于，所述脑图像测量单元求出对象切片整体的像素的信号值的直方图，求出去除高值的预定像素后的信号的最大值作为图像整体的有效最大值。根据本发明第5方式，求出对象切片整体的像素的信号值的直方图，求出去除高值的预定像素后的信号的最大值作为图像整体的有效最大值，因此，能准确求出作为代表图像整体的值的图像整体的有效最大值。另外，在本发明第6方式中，对于本发明第I至第5方式的任一方式的医用图像处理装置，其特征在于，所述头盖 测量单元对选择到的所述各对象切片设定预定数量的通过脑中心部的线段，按各线段的每一个，将线段上的像素的信号值分布分割成预定数量的区域，将两端的预定数量的区域作为头盖周边区域，求出该头盖周边区域中的信号最大值在全部对象切片的全部线段中的平均值作为头盖周围的峰值平均值。根据本发明第6方式，对各对象切片，将通过脑中心部的多个线段各自的信号值分布分割成预定数量的区域，将两端的预定数量的区域作为头盖周边区域，求出该头盖周边区域中的信号最大值在全部对象切片的全部线段中的平均值作为头盖周围的峰值平均值，因此，能准确求出作为代表头盖周边区域的值的头盖周围的峰值平均值。另外，在本发明第7方式中，对于本发明第I至第6方式的任一方式的医用图像处理装置，其特征在于，所述高信号值抑制处理单元通过执行对信号值不足预定值的像素不改变其信号值、对信号值为所述预定值以上的像素将信号值设为预定值的处理，由此，相对地实现高信号值抑制处理。根据本发明第7方式，执行对信号值不足预定值的像素不改变其信号值、对信号值为预定值以上的像素将信号值设为预定值的处理，因此，能抑制输入图像的高信号值。另外，在本发明第8方式中，对于本发明第I至第6方式的任一方式的医用图像处理装置，其特征在于，所述高信号值抑制处理单元通过执行对信号值不足预定值的像素不改变其信号值、对信号值为所述预定值以上的像素将信号值变换成比原值小的处理，由此，相对地实现高信号值抑制处理。根据本发明第8方式，执行对信号值不足预定值的像素不改变其信号值、对信号值为所述预定值以上的像素将信号值变换成比原值小的处理，因此，能抑制输入图像的高
信号值。另外，在本发明第9方式中，对于本发明第8方式的医用图像处理装置，其特征在于，所述高信号值抑制处理单元通过执行对信号值不足预定值的像素不改变其信号值、对信号值为所述预定值以上的像素利用斜率大于等于零小于I的一次函数来变换信号值的处理，由此，相对地实现高信号值抑制处理。根据本发明第9方式，执行对信号值不足预定值的像素不改变其信号值、对信号值为预定值以上的像素利用斜率大于等于零小于I的一次函数来变换信号值的处理，因此，能抑制输入图像的高信号值。另外，组织分离处理结果构成自然的图像。发明效果根据本发明，实现如下效果，即能判定取得的脑图像的状态，并调整为适于组织分离处理的脑图像。
图1是表示本发明所涉及的一实施方式的医用图像处理装置的概要的框图。图2是表示基于本实施方式的医用图像处理方法的基本处理步骤的流程图。图3是表示基于本实施方式的医用图像处理方法的事先处理步骤的流程图。图4是示意性地表示脑的切片图像与体素的特征的概念图。图5是表示准备处理的细节的流程图。图6是表示对象切片的选择处理的状态的图。图7是表示脑实质与头盖周边区域的信号值分布状态的图。图8是表示脑实质的有效最大值、头盖周围的峰值平均值的测量处理的状态的图。图9是表不灰白质、白质的不同年龄模板的概念图。图10是表示灰白质、白质的性别模板的概念图。图11是示意性地表示空间标准化的特征的概念图。图12是表示非线性变换的特征的概念图。图13是表示空间标准化与平滑化的结果的概念图。图14是表示每个体素的比较检定的特征的概念图。图15是表示基于ROI的解析特征的概念图。图16是表示制作ROI时的特征的概念图。图17是表不检查的阴阳性与有无疾病的关系的图表。图18是表示ROC曲线的一例的线图。图19是表示基于I个参数的识别方法的线图。
具体实施例方式下面，参照附图来详细说明本发明的优选实施方式。图1是表示本发明一实施方式中的医用图像处理装置的框图。图1所示的本实施方式的医用图像处理装置具有用户界面10、图像/统计处理部20、数据库部30。用户界面10具有如下功能:输入MRI图像作为输入图像的图像输入功能12、以及显示由图像/统计处理部20处理的结果的结果显示功能14，图像/统计处理部20具有如下功能:处理从用户界面10输入的MRI图像的图像处理功能22、进行各种统计运算等的统计处理功能24、以及对输入图像判定是否需要高信号值抑制的是否需要高信号值抑制的判定功能26。另外，数据库部30中保存有图像/统计处理部20在后述的处理中使用的白质脑图像模板32、灰白质脑图像模板34、健康者图像数据库36、疾病特异ROI38等。图2是表示涉及本实施方式的医用图像处理装置的处理概要的流程图，图3是表示基于本实施方式的医用图像处理方法的事先处理步骤的流程图。图2中，首先，当输入被检查者的MRI脑图像时，对该脑图像进行预定的处理，判定是否应进行高信号值抑制处理，并根据判定结果，进行高信号值抑制处理(SI)。在SI中，有时也因判定结果不同而不进行高信号值抑制处理。该SI中的处理是成为本发明特征的部分。对于SI中的处理细节将在后文进行说明。接着，对高信号值抑制处理后的脑图像(在不进行高信号值抑制处理的情况下为输入的MRI脑图像)进行用于校正空间错位的定位(S2)。接下来，根据定位后的脑图像，制作利用组织分离处理抽取出白质的白质图像与同时抽取出灰白质的灰白质图像(S3)。接着，对制作出的两个图像进行适用后述的DARTEL算法的空间标准化(S4)，并且，对该标准化后的白质及灰白质图像分别进行平滑化(S5)。之后，在平滑化后的两个图像与健康者的白质及灰白质图像之间分别进行统计比较(S6)，再进行基于ROI的解析(S7)，输出解析结果作为诊断结果，用于诊断的支援。在本实施方式中，上述SI S7的各处理能在由计算机构成的所述图像/统计处理部20中通过程序来实施，并且，下面说明的Sll S14的各处理的白质及灰白质模板的制作也同样能通过程序来实施。在执行上述SI S7的基本处理流程之前，在图3的Sll S14中进行制作在S4的空间标准化中使用的模板的事先处理。首先，从尽可能多的健康者输入Tl强调MRI脑图像(图中为被检查者图像)。对从各被检查者取得的MRI脑图像进行预处理。具体地说，如图4中示出脑整体与切出其一部分而得的切片图像的影像那样，以包含被检查者的脑整体的方式摄像成预定厚度的切片状的例如100 200张Tl强调MRI图像作为输入。另外，进行切片图像的重新取样，以使各切片图像中体素(voxel)的各边长度事先相等。这里，体素是具有“厚度”的图像的坐标单位，相当于2维图像中的像素。在输入进行这种预处理的MRI脑图像之后，检查其切片图像的撮像方向或分辩率是否符合事先设定于系统中的条件。如上所述，在确认MRI脑图像符合设定条件地被输入的情况下，进行准备处理(Sll)。所谓准备处理是在模板制作中，为了准确进行组织分离处理，而判定是否应对输入图像进行高信号值抑制处理，在判定为需要的情况下，进行高信号值抑制处理。图5是表示准备处理的细节的流程图。首先，进行从输入图像中选择用于判定对象的切片图像作为对象切片的处理(S21)。图6中示出对象切片的选择处理的状态。作为对象切片，期望选择不包含头盖周边区域与脑实质以外的高信号值的部分。具体地说，对于各切片图像，在图像内设定预定的线段，求出位于该线段上且信号值为预定值以上的像素中间隔最远的像素间的长度。设定通过撮影到的头盖周边区域与脑实质的大部分的线段 作为在图像内设定的线段。通常，因为以头盖周边区域与脑实质为中心摄影四边形的图像，所以只要使用图像的对角线作为上述线段即可。作为与信号值进行比较的上述预定值，可适当设定，例如可以设为向图像整体的信号值的平均乘以预定常数而得的值。信号值为预定以上的像素表示头盖周边区域，信号值为预定以上的像素中间隔最远的像素间的长度表示线段与头部重合的长度。如上所述，期望选择不包含头盖周边区域与脑实质以外的高信号值的切片图像作为对象切片，该部分存在于长度为头部最大部分的80 95%的位置。因此，在本实施方式中，确定上述像素间的长度最大的切片图像，选择该像素间长度为确定的切片图像的80 95%长度的切片图像作为对象切片。一旦选择了对象切片，则接着进行脑实质的有效最大值的测量处理(S22)。当进行脑实质的有效最大值的测量时，首先取得脑实质的信号值分布。在取得脑实质的信号值分布时，同时还取得头盖周边区域的信号值分布。图7中示出脑实质与头盖周边区域的信号值分布的状态。作为信号值分布的取得处理，首先对选择的各对象切片设定预定数量的通过脑中心部的线段。在本实施方式中，如图7(b)所示，设定有10条线段。之后，对各线段描绘线段上像素的信号值。再根据线段将信号值分布分割成预定数量的区域，将中央的预定数量的区域设为脑实质区域，将两端的预定数量的区域设为头盖周边区域。在本实施方式中，如图7(a)右侧的信号值分布所示，分割成10个区域，设中央的4区域为脑实质区域，设各两端分别I个区域(共2个区域)为头盖周边区域。信号值分布中的A、B对应于图7(a)左侧所示的脑图像上的A、B。接着，求出全部对象切片的脑实质区域的信号值的直方图，求出其中去除受到高信号伪影(MRI的噪声)影响的部分时的最大值作为脑实质有效最大值Bmax。对于将哪部分作为受到高信号 伪影的影响的部分而去除，可以适当进行设定，具体地说，去除上述直方图中高值访他％的像素。在本实施方式中，设定为thHa = 3。图8(a)中示出脑实质的有效最大值的测量处理的状态。如图8(a)的下部直方图所示，去除画有影线的高值3%像素而得的箭头所示部位的值为脑实质有效最大值Bmax。接着，进行图像整体的有效最大值的测量处理(S23)。具体地说，求出四边形的对象切片整体的像素的信号值的直方图，求出其中去除受到高信号伪影影响的部分时的最大值作为图像整体的有效最大值Imax。对于将哪部分作为受到高信号伪影影响的部分而去除，可以适当进行设定，具体地说，去除上述直方图中高值thHb%的像素。在本实施方式中，设定为thHb = I。接着，进行头盖周围的峰值平均值的测量处理(S24)。具体地说，从上述S22中设定的头盖周边区域中抽取最大的信号值。头盖周边区域中的最大的信号值为图7(a)右侧的信号值分布中由〇包围的两个部位。对全部对象切片的全部线段进行该抽出处理，算出抽取的信号值的平均值。图8(b)中示出头盖周围的峰值平均值的测量处理的状态。在本实施方式的情况下，如图7(a) (b)所示，对I个线段抽取2个信号值，因此，对I个对象切片抽取20个信号值。因此，算出对象切片数X20个信号值的平均作为头盖周围的峰值平均值 Pave。接着，进行是否需要高信号值抑制的判定处理(S25)。具体地说，进行头盖区域判定与图像整体判定这2个判定，在满足该任一判定条件的情况下，判定为需要高信号值抑制，在任一判定条件均不满足的情况下，判定为不需要高信号值抑制。头盖区域判定进行头盖周边区域的信号值是否比脑实质区域高一定比例的判定。具体地说，执行根据以下〔数式I〕的处理。〔数式I〕thl X (脑实质有效最大值Bmax) < (峰值平均值Pave)上述〔数式I〕中，thl是确定允许头盖周边区域比脑实质高至何程度比例的系数，可适当进行设定，在本实施方式中为thl = 1.8。判定是否满足上述〔数式I〕所示的条件，在满足的情况下，判定为需要高信号值抑制。这是因为在头盖周边区域的峰值平均值比脑实质有效最大值高一定比例的情况下，需要抑制头盖周边区域的高信号值。另一方面，图像整体判定进行图像整体的信号值是否比脑实质区域高一定比例的判定。具体地说，执行根据以下〔数式2〕的处理。〔数式2〕th2 X (脑实质有效最大值Bmax) < (图像整体的有效最大值Imax)上述〔数式2〕中，th2是确定允许还包含脑实质以外的图像整体比脑实质高至何程度比例的系数，可以适当进行设定，在本实施方式中为th2 = 1.7。判定是否满足上述〔数式2〕所示的条件，在满足的情况下，判定为需要高信号值抑制。这是因为在图像整体的有效最大值比脑实质有效最大值高一定比例的情况下，需要抑制脑实质以外的高信号值。在满足头盖区域判定与图像整体判定任一方的条件的情况下，即满足上述〔数式I〕与〔数式2〕任一方的情况下，前进到高信号值抑制处理(S26)。相反，在头盖区域判定与图像整体判定任一条件均不满足的情况下，即，在上述〔数式I〕与〔数式2〕任一均不满足的情况下，不进行高信号值抑制处理，终止准备处理。这是因为在头盖区域判定与图像整体判定的任一条件均不满足的情况下，脑实质以外的高信号值不会高至妨碍组织分离处理，不需要高信号值抑制。S25中头盖区域判定与图像整体判定哪个在先执行均可，在在先执行一方满足判定条件的情 况下，进入到高信号值抑制处理(S26)而不对另一方判定条件进行判断。对高信号值抑制处理(S26)进行说明。高信号值抑制处理是将输入图像的像素中信号值高的像素值抑制为相对不高的处理。因此，只要是将信号值高的像素值抑制为与其他像素相比相对不高的处理，就能使用各种方法。在本实施方式中，作为第I方法，通过执行对高至某一定值以上的像素设其信号值固定、对不足某一定值的像素不改变其信号值的处理，来相对实现高信号值抑制处理。在本实施方式中，通过执行根据以下〔数式3〕的处理，将输入图像的各像素(x，y)的信号值s (X, y)校正为信号值s’ (X, y)。〔数式3〕S，(X，y) = th3不过，限于s (x, y) ^ th3的情况。上述〔数式3〕中，th3是成为抑制对象的像素的信号值的阈值，在本实施方式中，设定为th3 =(脑实质有效最大值Bmax)。另外，在本实施方式中，作为第2方法，也可通过执行对高至某一定值以上的像素利用斜率a大于零小于I的一次函数变换其信号值、对不足某一定值的像素不改变其信号值的处理，由此，相对地实现高信号值抑制处理。由此，能将某一定值以上的信号值的变化率抑制成比不足某一定值的信号值的变化率小，将信号值高的像素值抑制成相对不高。在第2方法中，因为使某一定值以上的信号值平缓上升，所以组织分离处理的结果为较自然的图像。在本实施方式中，作为第2手法，通过执行根据以下〔数式4〕的处理，将输入图像的各像素(X, y)中的信号值s(x, y)校正为信号值s’ (x, y)。〔数式4〕S，(X，y) = th3+a (s (x、y) _th3)不过，限于s (x, y) ^ th3的情况。上述〔数式4〕中，th3是成为抑制对象的像素的信号值的阈值，在本实施方式中，设定为th3=(脑实质有效最大值Bmax)。另外，在本实施方式中，为a = 0.2的常数。此外，在本实施方式中，利用根据上述〔数式3〕或〔数式4〕的处理来实现高信号值抑制处理，但只要是对某一定值以上的信号值像素将其信号值变换成比原值小的处理即可，也可使用Y校正等公知的其他方法来实现。当准备处理结束时，接着，返回图3的流程图，执行以后的处理。在以后的S12 S14中，对于进行了高信号值抑制的图像，将校正后的图像作为对象，对于未进行高信号值抑制的图像，将输入图像作为对象进行处理。首先，进行空间定位处理(S12)。这相当于为了提高检查对象的脑图像与标准的脑图像进行比较时的精度而利用线性变换(仿射变换)进行空间位置与角度的校正。在以上空间定位结束后，进行组织分离处理(S13)，分别制作抽取白质与灰白质而得的白质图像与灰白质图像。
·
因为输入的Tl强调MRI脑图像中包含对应于神经纤维的呈现高信号值的白质、对应于神经细胞的呈现中间信号值的灰白质、呈现低信号值的脑脊液这3种组织，所以着眼于该差别，进行分别抽取白质与灰白质的处理。该处理在所述专利文献I中进行了说明。并且，本发明中进行抽取精度比所述专利文献I中使用的方法更高的综合型组织分离处理。综合型组织分离处理是将标准化、组织分离、信号不均一的校正编入I个模块中的组织分离方法。具体内容记载在非专利文献3中。在综合型组织分离处理中，还具有如下特征:除制作白质图像和灰白质图像外，还制作表不MRI图像坐标与标准脑坐标的对应关系的变换场。变换场在后述的标准化中使用。这样，通过事先从多个健康者的MRI脑图像进行组织分离，分别取得3维抽取白质、灰白质而得的白质图像、灰白质图像作为多个采样。如上所述，分别组织分离许多(多个)健康者的MRI脑图像后制作抽取白质而得的白质图像作为采样，通过在制作的全部采样间进行空间标准化来制作白质模板(步骤14)。同样，分别组织分离多个健康者的MRI脑图像后制作抽取灰白质而得的灰白质图像作为采样，通过在制作的全部采样间进行空间标准化来制作灰白质模板。在这里执行的空间标准化中适用DARTEL算法。为了解决所述现有的VBM方法的精度差，JohnAshbumer开发了新的VBM方法、即DARTEL(Diffeomorphic Anatomical Registration using Exponentiated Lie algebra:通过指数化李代数的微分同胚的解剖登记)(参照非专利文献4)。DARTEL与现有的VBM方法相比，空间标准化精度好，作为提高老年痴呆病的图像统计解析的诊断性能的技术而受到期待。另外，因为进行DARTEL的空间标准化比以往的方法更精密，所以不仅以前作为评价对象的灰白质，还能将白质作为评价对象。在所述S14的模板制作处理中，对白质与灰白质分别制作与年龄、性别等被检查者属性对应的不同层次的模板，作为所述白质脑图像模板32、灰白质脑图像模板34保存在所述数据库部30中。如图9所示，所述白质模板与灰白质模板按不同年龄来制作。灰白质(A)、白质(C)根据54岁 69岁的各健康者的图像制作而成，灰白质(B)、白质⑶根据70岁 84岁的各健康者的图像制作而成。如该图所示，可知模板因年龄不同而不同，尤其是白质的情况下，在(C)与(D)中从左起第6个模板中箭头所示的脑室大小大不相同。另外，如图10所示，所述白质模板与灰白质模板按性别来制作。若按男女来比较该图中以四边形包围表示的位置的模板，则尤其在脑室形状上可以发现因性别的不同而存在差异。如上所述那样制作的白质与灰白质的模板以按不同年龄、性别进行准备为前提，实施根据所述SI 7的诊断支援信息制作处理。此外，下面，将上述白质与灰白质的模板称为DARTEL模板。首先，作为被检查者图像，输入摄像成预定厚度的切片状的Tl强调MRI图像，进行切片图像的重新取样，使各切片图像中体素(voxel)的各边长度事先相等。接着，与所述事先处理的Sll的情况一样进行所述SI的高信号值抑制处理。接着，与所述事先处理的S12的情况一样进行所述S2的空间定位处理。在以上空间定位结束之后，进行所述S3的组织分离的处理。该组织分离与所述S13的情况一样，抽取白质与灰白质之后，分别制作被检查者的白质图像与灰白质图像。对以上制作的被检查者的白质图像与灰白质脑图像进行所述S4的空间标准化处理。在这里进行的空间标准化中，与所述S14的情况一样，适用DARTEL算法。该空间标准化处理是为了吸收个人之间存在的脑图像的解剖学上的差异，而对脑整体大小进行大局上的校正，并对部分大小进行局部校正的处理。下面，为了方便，主要围绕灰白质进行说明，但在白质的情况下也进行实质上相同的处理。S4的DARTEL中的空间标准化处理由以下3个步骤的处理构成。(S4-1)初始位置决定处理(S4-2)向DARTEL模板的变换处理(S4-3)向标准脑模板的变换处理在S4-1的初始位置决定处理中，使用向上述综合型组织分离处理得到的标准脑的变换场，对灰白质图像、白质图像进行决定初始位置的处理。该处理中，其特征在于，因为进行刚体(rigid body)变换，所以图像形状不发生变化。在S4-2的向DARTEL模板的变换处理中，对实施S4-1处理的图像使用DARTEL算法，使形状与DARTEL模板吻合。在S4-3的向标准脑模板的变换处理中，进行使S4-2中得到的与DARTEL模板吻合的图像与标准脑模板匹配的处理。事先求出从DARTEL模板向标准脑模板的变换场，使用该变换场来进行向标准脑坐标系的变换。在S4-2与 S4-3的处理中，通过保持各体素的信号值的合计不变地进行标准化来保持体积信息，因此，能在标准化后测量体积。S4-1中进行线性变换，在S4-2与S4-3中进行线性变换与非线性变换。若以步骤3-2为例进行说明，则如图11中概括性地示出处理特征那样，使用线性变换与非线性变换进行图像处理以使与从所述数据库部30读出的、所述S14中制作的平均灰白质脑图像模板34的误差的平方和最小。在该空间标准化处理中，首先进行基于线性变换的位置或大小、角度的大局校正，接着利用非线性变换来进行局部的凹凸等形状的校正。这里进行的线性变换是与所述S2的定位一样的仿射变换。另外，非线性变换如图12中示出处理的影像那样，就X方向、y方向分别推定表示局部位移的变形场，根据该变形场来进行原图像的变换。S4-2的处理是以所述S14制作的模板为雏形、匹配处理对象图像的处理，使用的模板适用DARTEL算法并以高精度制作，所以其形状尖锐。因此，由于空间标准化，各个处理对象匹配为无个体差的近似形状，所以个体间的形状接近相同，但因为萎缩等反映于局部密度，因此能够提高空间标准化的精度。如上所述，对实施了空间标准化的白质图像与灰白质脑图像(以下也称为标准化脑图像)进行所述S5的图像平滑化处理。该平滑化处理是用于提高上述标准化脑图像的S/N比并且使后面进行比较时作为标准使用的健康者的图像群与图像的smoothness相等的处理，使用3维高斯核来进行。该平滑化中使用的滤波器的FWHM(半值宽度)为8mm左右。具体地说，如所述专利文献I中所述，进行3维脑图像与3维高斯函数的3维卷积(convolution)。这通过依次在x、y、z各方向上进行I维卷积来实现。通过这样进行平滑化，能使空间标准化处理中不完全一致的个体差降低。图13中示出对在S3中分离出的白质图像与灰白质图像进行S4的空间标准化处理与S5的平滑化处理后的结果图像。在S4的处理中，对脑的体积信息进行了保存。因此，也可在实施下一浓度值校正之前，对白质和灰白质的处理结果图像，测量整体或后述的感兴趣区域(regions ofinterest:R0I)的积分值作为体积并活用作诊断支援信息。对如上所述那样进行了图像平滑化而得的标准化脑图像省略了图示，但为了与之后进行比较时用作标准的健康者的图像群中的体素值分布匹配，进行校正脑整体的体素值的浓度值校正。之后，进行所述S6的统计比较。这里，进行通过以上SI 5各处理进行标准化而得的被检查者的灰白质(白质)的MRI脑图像、与事先收集后作为健康者图像数据库36保存在所述数据库部30中的健康者的同一 MRI脑图像群的比较检定。期望使用的健康者图像群由接近被检查者年龄的图像群构成。具体地说，如图14中示出影像那样，以体素为单位与这种健康者图像群进行1:N(N为健康者图像的总数)的比较检定，检测在统计上发现有意差(推定为异常)的体素。首先，对全部体素算出分别由下式表示的Z分数。数I
权利要求
1.一种医用图像处理装置，其特征在于，具有: 对象切片选择单元，从由多个切片图像构成的脑图像中，选择作为处理对象的切片图像作为对象切片； 脑实质测量单元，进行脑实质的像素的信号值的有效最大值、即脑实质有效最大值的测量处理； 脑图像测量单元，进行脑图像整体的像素的信号值的有效最大值的测量处理； 头盖测量单元，测量处理成为头盖周边区域的峰值的像素的信号值的平均作为头盖周围的峰值平均值； 抑制处理判定单元，根据测量出的脑实质的有效最大值、脑图像整体的有效最大值、头盖周围的峰值平均值，判定是否需要高信号值抑制处理；以及 高信号值抑制处理单元，在由所述抑制处理判定单元判定为需要的情况下，对所述脑图像执行高信号值抑制处理。
2.如权利要求1所述的医用图像处理装置，其特征在于: 所述抑制处理判定单元执行头盖区域判定以及图像整体判定，该头盖区域判定是头盖周围的峰值平均值是否比所述脑实质有效最大值高一定比例的判定，该图像整体判定是图像整体的信号值是否比所述脑实质有效最大值高一定比例的判定， 在所述头盖区域判定、图像整体判定的任一方满足条件的情况下，判定为需要高信号值抑制处理。
3.如权利要求1或2所述的医用图像处理装置，其特征在于: 所述对象切片选择单元对各切片图像在图像内设定预定的线段，求出位于该线段上且信号值为预定值以上的像素中间隔最远的像素间的长度，将该像素间的长度与全部切片图像中最长的长度相比较，将成为预定比率以上的切片图像选为对象切片。
4.如权利要求1至3的任一项所述的医用图像处理装置，其特征在于: 所述脑实质测量单元对选择到的所述各对象切片设定预定数量的通过脑中心部的线段，按各线段的每一个，将线段上像素的信号值分布分割成预定数量的区域，将中央的预定数量的区域作为脑实质区域，求出该脑实质区域中的信号值的直方图，求出去除高值的预定像素后的信号最大值作为脑实质的有效最大值。
5.如权利要求1至4的任一项所述的医用图像处理装置，其特征在于: 所述脑图像测量单元求出对象切片整体的像素的信号值的直方图，求出去除高值的预定像素后的信号最大值作为图像整体的有效最大值。
6.如权利要求1至5的任一项所述的医用图像处理装置，其特征在于: 所述头盖测量单元对选择到的所述各对象切片设定预定数量的通过脑中心部的线段，按各线段的每一个，将线段上的像素的信号值分布分割成预定数量的区域，将两端的预定数量的区域作为头盖周边区域，求出该头盖周边区域中的信号最大值在全部对象切片的全部线段中的平均值作为头盖周围的峰值平均值。
7.如权利要求1至6的任一项所述的医用图像处理装置，其特征在于: 所述高信号值抑制处理单元通过执行对信号值不足预定值的像素不改变其信号值、对信号值为所述预定值以上的像 素将信号值设为预定值的处理，由此，相对地实现高信号值抑制处理。
8.如权利要求1至6的任一项所述的医用图像处理装置，其特征在于: 所述高信号值抑制处理单元通过执行对信号值不足预定值的像素不改变其信号值、对信号值为所述预定值以上的像素将信号值变换成比原值小的处理，由此，相对地实现高信号值抑制处理。
9.如权利要求8所述的医用图像处理装置，其特征在于: 所述高信号值抑制处理单元通过执行对信号值不足预定值的像素不改变其信号值、对信号值为所述预定值以上的像素利用斜率大于等于零小于I的一次函数来变换信号值的处理，由此，相对地实现高信号值抑制处理。
10.一种用于使计算机作为权利要求1至9的任一项所述的医用图像处理装置发挥功能 的程序。
全文摘要
本发明提供一种医用图像处理装置，其能判定取得的脑图像的状态，并调整为适于组织分离处理的脑图像。医用图像处理装置从由多个切片图像构成的脑图像中选择作为处理对象的切片图像作为对象切片，根据选择到的对象切片，执行脑实质的有效最大值的测量处理、脑图像整体的有效最大值的测量处理、头盖周围的峰值平均值的测量处理，根据测量到的脑实质的有效最大值、脑图像整体的有效最大值、头盖周围的峰值平均值，判定是否需要高信号值抑制处理，在判定为需要高信号值抑制处理的情况下，对脑图像执行高信号值抑制处理。
文档编号A61B5/055GK103249358SQ20118005811
公开日2013年8月14日 申请日期2011年12月1日 优先权日2010年12月2日
发明者后藤智章, 小野彻太郎 申请人:大日本印刷株式会社