具体实施方式
图1为表示本发明的实施中优选的超声波诊断装置的整体结构的图。探头10为,相对于包含诊断对象在内的区域发送接收超声波的超声波探头。探头10具备发送接收超声波的多个振动元件,且多个振动元件通过发送接收部12而被实施发送控制从而形成发送信号束。此外,多个振动元件对从包含诊断对象在内的区域内获得的超声波进行接收,由此向发送接收部12输出所得到的信号,从而发送接收部12形成接收信号束而沿着接收信号束收集回声数据。
探头10优选为三维探头,所述三维探头在三维空间内对超声波束(发送信号束和接收信号束)进行扫描从而以立体的方式来收集回声数据。例如,通过使由一维排列的多个振动元件(1D阵列振子)以电的方式所形成的扫描面机械地进行移动,从而立体地扫描超声波束。此外,也可以对二维排列的多个振动元件(2D阵列振子)进行电控制从而立体地扫描超声波束。当然,也可以利用在断层面内扫描超声波束的二维的超声波探头。
当在三维空间内扫描超声波束以收集回声数据时,关于构成与该三维空间相对应的三维数据空间的多个体素的回声数据(体素数据),被存储于未图示的存储器等中。而且,对于构成三维数据空间的多个体素,在对象特定部20以后的各部中执行各种处理。在此,对上述各种处理进行说明。另外,关于图1所示的部分(结构),在以下的说明中利用图1中的符号。
图2为用于对对象特定部20中的处理进行说明的图。图2(A)图示了二值化处理。对象特定部20为,通过对构成三维数据空间的多个体素实施二值化处理,从而形成图2(A)所示的二值化处理后的图像数据。作为本实施方式中的诊断对象,优选为生物体内的卵泡。对象特定部20通过对用于二值化的阈值和各个体素的体素值(回声数据的大小)进行比较,从而对与卵泡F相对应的体素与除此之外的体素进行识别。而且,例如,将与卵泡F相对应的体素的体素值设为“1”,并将除此之外的体素的体素值设为“0”。图2(A)中用白色来描绘与卵泡F相对应的体素群,用黑色来描绘除此之外的相当于背景的体素群。
在生物体内,多个卵泡以互相紧贴的方式而密集在一起。因此,在超声波图像内,如图2(A)所示,多个卵泡F以互相连结的状态而被图像化,从而难以对各个卵泡F的大小或形状等单独地进行确认。在此,在本实施方式中,通过下文说明的各种处理而使多个卵泡F相互分割。另外,虽然在图2中以二维的方式绘制了各个图像数据,但各种处理是在三维数据空间内以三维的方式被执行的。
图2(B)图示了收缩分离处理。对象特定部20在构成三维数据空间的二值化处理后的体素数据内,即,在图2(A)所示的二值化图像数据内,对多个卵泡F实施收缩处理,从而如图2(B)所示那样,将多个卵泡F分离为卵泡F1~F3。对象特定部20反复执行n次(n为自然数)阶段性地使卵泡F收缩的收缩处理。在该各个阶段的收缩处理中,利用收缩处理用的过滤器,并且该过滤器在三维数据空间内的整个区域内被扫描。
图3为用于对三维数据空间100内的过滤器120的扫描进行说明的图。在图3中,三维数据空间100由xyz直角坐标系表示。此外,过滤器120成为如下结构,即,x轴方向、y轴方向和z轴方向上的各个长度相当于三个体素,从而相当于合计27个体素的立体结构。位于过滤器120的中心的体素为关注体素,而包围该关注体素的26个体素为周边体素。而且,以三维数据空间100内的全部体素成为关注体素的方式,使过滤器120在x轴方向、y轴方向和z轴方向上进行移动以在三维数据空间100内的整个区域内进行扫描。
在收缩分离处理中,在各扫描位置处,过滤器120内的26个周边体素中只要存在一个体素值为“0”的体素,则位于过滤器120的中心的关注体素的体素值将被设为“0”。例如,关注体素为“1(卵泡)”,在周边体素中只要存在一个体素值为“0(背景)”的体素,则该关注体素将被变更为“0(背景)”。而且,通过在三维数据空间100内的整个区域内粗略地扫描过滤器120,且在各扫描位置处实施过滤处理,从而结束一阶段的收缩处理。另外,与关注体素相关的体素值的变更,是在过滤器120于三维数据空间100内的整个区域内粗略地被扫描之后执行的。即,在过滤器120被扫描的中途,不实施体素值的变更,并且在任一扫描位置处均将变更前的体素值作为对象来执行过滤处理。
如此,当结束一阶段的收缩处理,并根据其结果来实施体素值的变更时,对由该变更后的体素值构成的三维数据空间100,执行第二阶段的收缩处理。在第二阶段的收缩处理中,也执行与第一阶段相同的过滤处理。即,在各个扫描位置,过滤器120内的26个周边体素中只要存在一个体素值为“0”的体素,则位于过滤器120的中心的关注体素的体素值将被设为“0”,并且在过滤器120于三维数据空间100内的整个区域内粗略地被扫描之后,变更体素值。
对象特定部20反复执行n次(n为自然数)该阶段性的收缩处理。该反复次数n根据各个体素的大小或过滤器的大小等而被适当地确定,例如被设定为10以内的程度。当然,也可以采用如下方式,即,用户能够对次数n进行调节。
另外,在二维的情况下,只需代替图3所示的过滤器120,而利用纵向与横向的长度分别相当于三个体素从而与合计九个体素相对应的二维过滤器,并且将位于中心的体素设为关注体素,将包围所述关注体素的八个体素设为周边体素即可。
返回图2,当通过收缩处理而如图2(B)所示那样被分离为多个卵泡F1~F3时,对象特定部20将在构成三维数据空间的收缩处理后的体素数据内,即,在图2(B)所示的收缩处理后的图像数据内,实施分类处理,并对多个卵泡F1~F3分配互不相同的标签。作为分类处理,可以使用公知的方法,例如检测出三维数据空间内具有相同的体素值的多个体素的群,并对每个群赋予标签号。例如,如图2(C)所示,对作为体素值为“0”的群的背景部分分配标签0,对作为体素值为“1”的群的卵泡F1~F3,分别分配标签1~3。
当分类处理被实施时,对象特定部20将在构成三维数据空间的分类处理后的体素数据内,即,在图2(C)所示的分类处理后的图像数据内,对多个卵泡中的各个卵泡实施膨胀处理。在该膨胀处理中,在从各个卵泡所得到的膨胀部分中,分配有该卵泡的标签,而且,在由于被实施膨胀处理而相互重叠的膨胀部分(膨胀的卵泡)彼此的重复部分处形成边界,且多个卵泡的大小被复原。如此,如图2(D)所示,在与互不相同的标签相对应的卵泡之间形成有边界(背景像素),且各个卵泡的大小被复原至收缩处理前(刚实施二值化处理后)的大小。
对象特定部20反复执行n次(n与收缩处理的次数相同)阶段性地使卵泡F膨胀的膨胀处理。在该各个阶段的膨胀处理中,利用膨胀处理用的过滤器,并且该过滤器在三维数据空间内的整个区域内被扫描。在膨胀处理中,也利用图3所示的相当于合计27个体素的立体过滤器120,并且位于过滤器120的中心的体素被设为关注体素,包围所述关注体素的26个体素被设为周边体素。而且,以三维数据空间100内的全部体素成为关注体素的方式,使过滤器120在x轴方向、y轴方向和z轴方向上进行移动以在三维数据空间100内的整个区域内进行扫描。但是,在收缩处理和膨胀处理中,过滤处理有所不同。
图4为用于对膨胀处理中的过滤处理进行说明的图,在图4中图示了与形成边界且进行膨胀的处理(膨胀边界处理)中的体素值的变更相关的条件表。在膨胀边界处理中,参照各个体素的标签值。
在位于过滤器120(图3)的中心的关注体素为标签0(背景)的情况下,当26个周边体素全部为标签0(背景)时,关注体素将被设为标签0。此外,在关注体素为标签0(背景)的情况下,当26个周边体素中存在标签0以外(卵泡)的标签,且这些标签全部为相同的标签N(同一卵泡)时,关注体素将被变更为标签N。即,标签N的卵泡被膨胀。
而且,在关注体素为标签0(背景)的情况下,当26个周边体素中存在标签0以外(卵泡)的标签,且这些标签中包含互不相同的标签号(互不相同的卵泡)时,关注体素将被设为标签0。即,关注体素被维持为标签0,而成为互不相同的卵泡的边界。
另一方面,在位于过滤器120的中心的关注体素为标签M(卵泡)的情况下,无论周边体素的状况如何,关注体素均被维持为标签M。
而且,通过在三维数据空间100内的整个区域内粗略地扫描图3所示的过滤器120,且在各扫描位置,根据图4所示的条件而实施过滤处理,从而结束一阶段的膨胀边界处理。另外,与关注体素相关的标签值的变更,是在过滤器120于三维数据空间100内的整个区域内粗略地被扫描之后执行的。即,在过滤器120被扫描的中途,不实施标签值的变更,并且在任一扫描位置处均将变更前的标签值作为对象来执行过滤处理。
如此,当结束一阶段的膨胀边界处理,并根据其结果而实施标签值的变更时,对由该变更后的标签值构成的三维数据空间100,执行第二阶段的膨胀边界处理。在第二阶段的膨胀边界处理中,也执行与第一阶段相同的过滤处理。即,在各扫描位置处,根据图4所示的条件来实施过滤处理,并且在过滤器120于三维数据空间100内的整个区域内粗略地被扫描之后,变更标签值。
对象特定部20反复执行n次该阶段性的膨胀边界处理。优选为,该反复次数n与收缩处理的反复次数n相同。如此,如图2(D)所示,在与互不相同的标签相对应的卵泡之间形成边界(背景像素),并且各个卵泡的大小被复原至收缩处理前的大小。
另外,在膨胀边界处理中,在二维的情况下,只需代替图3所示的过滤器120,而利用纵向和横向的长度分别相当于三个体素从而与合计九个体素相对应的二维过滤器,并且将位于中心的体素设为关注体素,且将包围所述关注体素的八个体素设为周边体素即可。
在本实施方式中,如图2所示,以相互紧贴的方式而密集在一起的多个卵泡被相互分离并被特定。而且,如图2(D)所示,当对各个卵泡对应附加单独的标签时,将能够通过各个标签而对所对应的卵泡进行特定,从而能够针对于每个标签而计算与对应于该标签的各个卵泡的大小或形状相关的测量值,例如,可以针对于每个标签而计算对应于该标签的各个卵泡的体积、长轴的长度或短轴的长度等。
另外,对各个卵泡进行特定时,例如,通过用户对所需的标签进行指定,从而能够对与该标签相对应的卵泡进行特定。此外,由于各个卵泡被分离,因此在显示有多个卵泡的图像上,用户可以对光标等显示方式进行操作以指定所需的卵泡,由此仅显示所特定的卵泡的图像。
在本实施方式中,在对由用户所指定的卵泡进行显示时,形成与该卵泡的形态相对应的显示图像。而且,作为卵泡的形态所涉及的特征量,图1中的三轴运算部30对卵泡的长轴和两个短轴进行特定。
图5为用于对卵泡的长轴和两个短轴进行说明的图。在本实施方式中,对于所特定的卵泡F,如图5所示那样通过考虑外接于该卵泡F的最小的长方体,从而将该长方体的边的长度设为卵泡F的三轴长度。例如,图5所示的最长的边D1被设为卵泡F的长轴,与边D1正交的边D2和边D3被设为卵泡F的两个短轴。
返回图1,三轴运算部30在对卵泡的三轴进行特定时,利用主成分分析的方法。在被称为主成分分析的公知的方法中,将最佳地表示数据的波动的方向,即,数据的分散成为最大的方向设为第一主成分。而且,在本实施方式中,在进行主成分分析时,例如,利用在以下进行说明的公知的协方差矩阵。
为了得到协方差矩阵,通过数学式1来计算平均位置m。另外,在数学式1中,Pi为,关于构成卵泡的第i个像素(体素)的、三维数据空间(参照图3)内的坐标值,通过利用构成所特定的卵泡的i=1~N的全部像素(体素)的坐标值,来计算平均位置(重心位置)m。
【数学式1】
Pi=<xi,yi,zi>
而且,利用数学式1的平均位置m,来计算数学式2所表示的协方差矩阵C。另外,数学式2所表示的协方差矩阵C为3×3的矩阵,并成为具有数学式3所表示的六个独立成分的对称矩阵。
【数学式2】
【数学式3】
在使用了协方差矩阵C的主成分分析中,计算通过数学式2、数学式3式而得到的协方差矩阵C的固有矢量,并且与最大的固有值相对应的固有矢量被设为第一主成分。在本实施方式中,将利用协方差矩阵C而得到的第一主成分的方向设为卵泡的长轴。由此,穿过卵泡的重心并沿着卵泡的长边方向的长轴被特定。而且,将利用协方差矩阵C而得到的第二主成分和第三主成分各自的方向分别设为卵泡的两个短轴。例如,第二主成分的方向被设为短轴1,第三主成分的方向被设为短轴2。如此,作为卵泡的三轴,而特定了长轴和与之正交的两个短轴。
另外,虽然也可以在卵泡的图像数据内,沿着连结重心和距重心最远的像素的直线而设定长轴,但由于还考虑到最远的像素为噪声等的情况,因此优选为通过主成分分析而进行的长轴的设定。
当通过三轴运算部30而特定了三轴时,诊断坐标系设定部40将对以卵泡的形态为基准的诊断坐标系进行设定。诊断坐标系设定部40对将卵泡的三轴设为坐标轴的诊断坐标系进行设定。
图6为表示以卵泡为基准的诊断坐标系的图。诊断坐标系设定部40将图6所示的直角坐标系设为诊断坐标系,在所述直角坐标系中,将卵泡F的重心G的位置设为坐标原点,且将第一轴、第二轴和第三轴设为坐标轴,所述第一轴为第一主成分的方向,即,卵泡F的长轴的方向,所述第二轴为第二主成分的方向,即,卵泡F的一个短轴的方向,所述第三轴为第三主成分的方向,即,卵泡F的另一个短轴的方向。
返回图1,当通过诊断坐标系设定部40而设定了诊断坐标系时,坐标系匹配部50将使成为显示图像的基准的显示坐标系与诊断坐标系相互匹配,从而将卵泡的图像数据配置于显示坐标系上。
图7为用于对显示坐标系与诊断坐标系的匹配进行说明的图。在图7中,显示坐标系由XYZ直角坐标系来表示。显示坐标系为,成为形成显示图像时的基准的坐标系,且为使相对于三维数据空间(参照图3)的相对位置关系较为明确的坐标系。在本实施方式中,将三维数据空间(参照图3)的XYZ直角坐标系以原状态设为显示坐标系。
此外,在图7中,诊断坐标系为,由第一轴、第二轴和第三轴(参照图6)而被特定的坐标系。诊断坐标系的第一轴至第三轴为,根据三维数据空间内的各个像素(体素)的坐标,通过主成分分析并利用例如数学式1至数学式3而得到的轴,从而三维数据空间内的位置和方向被特定。因此,当将三维数据空间的XYZ直角坐标系设为显示坐标系时,相对于该显示坐标系的诊断坐标系的位置和方向将被特定。
图7的(A)图示了与被特定在显示坐标系(XYZ轴)上的卵泡F相关的诊断坐标系(第一轴至第三轴)的一个示例。由于诊断坐标系为以卵泡F的长轴和短轴为基准的坐标系,因此与显示坐标系内卵泡F的位置和方向相对应。
因此,坐标系匹配部50首先使诊断坐标系相对于显示坐标系而平行移动以使显示坐标系的原点与诊断坐标系的原点相一致。当然,与卵泡F相关的体素数据(图像数据)也与诊断坐标系一起平行移动。
图7的(B)图示了诊断坐标系被进行了平行移动的状态,诊断坐标系的原点被移动至显示坐标系的原点的位置处,由此,作为诊断坐标系的原点的、卵泡F的重心位置被移动至显示坐标系的原点。
接下来,坐标系匹配部50对与卵泡F的长轴相对应的轴,即,诊断坐标系的第一轴与显示坐标系的XYZ轴的各个轴进行比较,从而对XYZ轴中与第一轴之间的角度最小的轴进行特定。例如,对第一轴与XYZ轴的各个轴的内积进行比较,以对与第一轴之间的角度最小的轴进行特定。而且,以所特定的轴与第一轴重合的方式旋转移动诊断坐标系。例如,在特定了X轴以作为与第一轴之间的角度最小的轴的情况下,如图7的(C)所示那样,以第一轴重合在X轴上的方式旋转诊断坐标系,从而卵泡F的图像数据也被旋转。
而且,坐标系匹配部50对与卵泡F的短轴相对应的诊断坐标系的第二轴、与显示坐标系中的其余的轴进行比较,从而对与第二轴之间的角度最小的轴进行特定。例如,在第一轴与X轴重合的情况下,对其余的Y轴和Z轴中与第二轴之间的角度最小的轴进行特定。而且,以所特定的轴与第二轴重合的方式旋转移动诊断坐标系。例如,在特定了Z轴以作为与第二轴之间的角度最小的轴的情况下,如图7的(D)所示那样,以第二轴重合在Z轴上的方式旋转诊断坐标系,从而卵泡F的图像数据也被旋转。
另外,当诊断坐标系为直角坐标系时,通过使第一轴和第二轴与X轴和Z轴重合,从而第三轴将沿着Y轴而配置。在图7的(D)中,第三轴与Y轴以相互相同的方向被重合在一起。
返回图1,当通过坐标系匹配部50而使显示坐标系与诊断坐标系被匹配从而使卵泡的图像数据被配置在显示坐标系上时,显示图像形成部60将根据被配置于显示坐标系上的卵泡的图像数据,而形成该卵泡的显示图像,并使所形成的显示图像显示在显示部70上。在该显示图像的形成中,形成由以显示坐标系为基准的截面所表示的卵泡的断层图像。
图8为用于对以显示坐标系为基准的截面进行说明的图。图8中图示了,通过图7(D)所示的显示坐标系和诊断坐标系的匹配而被配置于显示坐标系上的卵泡F的图像数据。在图8中,截面A为包含显示坐标系的Z轴和X轴的平面,截面B为包含显示坐标系的Y轴和Z轴的平面,截面C为包含显示坐标系的X轴和Y轴的平面。
由于如图7(D)所示那样使显示坐标系与诊断坐标系匹配,因此与第一轴相对应的卵泡F的长轴被配置于X轴上,与第二轴相对应的卵泡F的短轴1被配置于Z轴上,与第三轴相对应的卵泡F的短轴2被配置于Y轴上。因此在图8中,截面A成为包含卵泡F的长轴和短轴1的截面,截面B成为包含卵泡F的短轴1和短轴2的截面,截面C成为包含卵泡F的长轴和短轴2的截面。
图9为表示显示图像62的具体示例的图。在构成图9所示的显示图像62的图像中,<3D>为与多个卵泡相关的三维图像。该三维图像是根据从三维空间内收集的回声数据(体素数据),并通过例如体绘制处理而形成的。通过采用用户能够变更该体绘制处理中的视点的位置的结构,从而能够得到从所需的方向放映出多个卵泡的三维图像。
例如,在图9所示的三维图像内,用户通过从所需的方向放映出多个卵泡,而且,通过使光标等的显示方式进行移动,以使该显示方式与所需的卵泡的位置相一致,从而能够使用户选择欲进行诊断的卵泡。如此,例如在图9的三维图像内卵泡F1被特定。另外,在三维图像内,也可以通过颜色等显示方式或标记等,以能够从视觉上对所特定的卵泡与其他的卵泡进行识别的方式而进行显示。
例如,当利用三维图像而由用户特定了卵泡F1时,通过三轴运算部30而特定卵泡F1的三轴(参照图5、数学式1至数学式3),且通过诊断坐标系设定部40而设定与卵泡F1的三轴相对应的诊断坐标系(参照图6),而且,通过坐标系匹配部50而使显示坐标系与诊断坐标系(参照图7)匹配。然后,通过显示图像形成部60而形成截面A至C(参照图8)中的卵泡F1的断层图像。
在图9的构成显示图像62的图像中,<截面A>为关于卵泡F1的截面A中的断层图像,<截面B>为关于卵泡F1的截面B中的断层图像,<截面C>为关于卵泡F1的截面C中的断层图像。由于显示坐标系与诊断坐标系被匹配,因此截面A显示包含卵泡F1的长轴和短轴1的截面,截面B显示包含卵泡F1的短轴1和短轴2的截面,截面C显示包含卵泡F1的长轴和短轴2的截面。
如此,在本实施方式中,通过由用户从多个卵泡中对所需的卵泡进行特定,从而形成包含该指定的卵泡的三轴的断层图像。因此,例如,能够减少用于用户对切面的设定等的烦杂的操作,优选能够省略用于切面的设定的操作。
此外,在显示坐标系与诊断坐标系的匹配中,通过使交差角度成为最小的坐标轴彼此重合,从而能够尽可能地减少诊断坐标系的旋转移动,由此能够尽可能地抑制由于该旋转移动而使用户感觉到的视觉上的不协调。
另外,关于所特定的卵泡F1,长轴的长度和两个短轴的长度或体积等测量值也可以作为显示图像62的一部分而被显示。此外,由于多个卵泡以相互分离的方式而被特定(参照图2),因此可以针对各个卵泡而对长轴的长度和两个短轴的长度或体积等测量值进行计算,并显示与多个卵泡相关的测量值的一览。而且,也可以由用户从这些测量值的一览中特定所需的卵泡,由此表示所特定的卵泡的截面。
以上,虽然对作为本发明的优选的实施方式的超声波诊断装置进行了说明,但例如可以采用如下方式,即,通过计算机来实现图1所示的对象特定部20、三轴运算部30、诊断坐标系设定部40、坐标系匹配部50和显示图像形成部60中的至少一个,并使该计算机作为超声波图像处理装置而发挥作用。
符号说明
10探头;20对象特定部;30三轴运算部;40诊断坐标系设定部;50坐标系匹配部;60显示图像形成部;70显示部。