CN101612047A - 图像显示装置、图像显示方法、存储介质以及程序 - Google Patents
图像显示装置、图像显示方法、存储介质以及程序 Download PDFInfo
- Publication number
- CN101612047A CN101612047A CN200910160229A CN200910160229A CN101612047A CN 101612047 A CN101612047 A CN 101612047A CN 200910160229 A CN200910160229 A CN 200910160229A CN 200910160229 A CN200910160229 A CN 200910160229A CN 101612047 A CN101612047 A CN 101612047A
- Authority
- CN
- China
- Prior art keywords
- image
- breast
- zone
- forcing
- mammary gland
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本发明的目的在于提供一种去掉乳腺的重叠,显示容易掌握乳腺构造的图像的图像显示装置、图像显示方法、存储介质以及程序。图像显示装置包括:存储与包含被检体的乳房的摄像区域有关的数据的存储单元(10)、根据数据生成与根据压迫板的位置信息分割的摄像区域内的多个分割区域分别对应的多个图像的处理单元(13);显示已生成的图像的显示单元(16)。
Description
本申请是申请号为200510091028.1、申请日为2005年8月3日、发明名称为“图像显示装置、图像显示方法、存储介质以及程序”的发明专利申请的分案申请。
技术领域
本发明涉及显示与在压迫板之间被压成扁平状的乳房有关的图像的图形显示装置、图像显示方法、存储介质以及程序。
背景技术
作为拍摄乳房的以往技术,有压迫乳房扩大乳腺,从1个方向进行2维图像收集的乳房专用X射线摄像。在乳房中的乳腺是由如葡萄那样的小室组成的立体性的构造,只原样拍摄时各室的重叠很多,不能进行详细的评价。为此通过夹着乳房使其变薄延伸,在减少乳腺的重叠的同时,用更少的X射线拍摄。
但是,难以判别和诊断在摄像时管球跟前的乳腺和检测器侧的乳腺重叠的肿瘤。
[非专利文献1]Tao Wu和另外10人著「Medical Physics,Vol.30No.3“Tomographic mammography using a limited number oflow-dose conebeam projection images”」American Association ofPhysicists in Medicine发行,2003年3月,p.365-380
发明内容
本发明的目的在于提供一种显示去掉乳腺的重叠且容易掌握乳腺构造的图像的图像显示装置、图像显示方法、存储介质以及程序。
本发明的第1方面,图像显示装置包括:存储与包含被检体的乳房的摄像区域有关的数据的部件;由上述数据生成与根据压迫板的位置信息分割的上述摄像区域内的多个分割区域分别对应的多个图像的部件;显示上述生成的图像的部件。
本发明的第2方面,在生成乳房的图像的图像显示装置中,包括:多个生成上述图像的局部区域图像的部件;在把上述生成的部分区域图像的全部显示在1个画面上的同时,显示与乳房的构造相关联的部分区域的记号的显示部件。
本发明的第3方面,在生成乳房的图像的图像显示装置中,包括:算出区域分割位置的部件;生成用上述区域分割位置分割的分割区域的图像的部件;显示上述生成的图像的部件,算出上述区域分割位置的部件以上述图像为基础,检测上述乳房的乳头的位置信息,根据上述被检测出的乳头的位置信息算出上述区域分割位置。
本发明的第4方面,在生成乳房的图像的图像显示装置中,包括:算出区域分割位置的部件;生成用上述算出的区域分割位置分割的分割区域的图像的部件,显示上述生成的图像的部件,上述算出装置以上述图像为基础,提取乳房内部的概略构造,根据该概略构造的位置算出分割位置。
本发明的第5方面,图像显示方法包括:从与包含被检体的乳房的摄像区域有关的数据中,生成与根据压迫板的位置信息分割的上述摄像区域内的多个分割区域分别对应的多个图像的步骤;显示上述生成的图像的步骤。
本发明的第6方面,是记录有用于让算出机实现以下部件功能的程序的计算机可以读取的存储介质,即,从与包含在被检体的乳房的摄像区域的数据中,生成与根据压迫板的位置信息分割的上述摄像区域内的多个分割区域分别对应的多个图像的部件;显示上述生成的图像的部件。
本发明的第7方面,是让计算机实现以下部件功能的程序,即,从与包含在被检体的乳房的摄像区域有关的数据中,生成与根据压迫板的位置信息分割的上述摄像区域内的多个分割区域分别对应的多个图像的部件;显示上述生成的图像的部件。
如果采用本发明,则可以提供显示去掉乳腺的重叠且容易掌握乳腺构造的图像的图像显示装置、图像显示方法、存储介质以及程序。
附图说明
图1是展示本发明的实施方式的X射线乳房摄像装置的结构图。
图2是展示被显示在图1的显示单元上的显示画面例子的图。
图3是展示图1的X射线乳房摄像装置的X射线摄像台的主要部分结构的斜视图。
图4是展示在本实施例中X射线摄像时的X射线管球的动作的图。
图5是展示图1的成像处理单元的可成像区域(摄像区域)的图。
图6是用于补足图1的投影处理单元的部分投影图像生成的图。
图7是展示图6的区域分割面(中央面)的图。
图8是用于补足图1的区域分割面设定处理单元的图像处理产生的区域分割面(中央面)的确定处理的图。
图9是展示图6的区域分割面(规定面)的图。
图10是展示图6的区域分割面(乳头面)的图。
图11是用于补足由图1的区域分割面设定处理单元的图像处理进行的区域分割面(乳头面)的确定处理的图。
图12是示例在本实施方式中形成在压迫板的表面上的沟的图。
图13是在图1的区域分割面设定处理单元12中,乳头提取处理的补足图。
图14是在图1的区域分割面设定处理单元12中,另一乳头提取处理的补足图。
图15是在图1的区域分割面设定处理单元12中,分割面设定处理的补足图。
图16是在图1的区域分割面设定处理单元12中,阈值确定处理的补足图。
图17是在图1的区域分割面设定处理单元12中,用于补足乳腺区域的提取处理的图。
图18是展示与乳腺层被置换为脂肪层的状态的乳房的体积数据有关的直方图。
图19是图18的乳腺区域的提取处理内的直方图的曲线近似处理的补足图。
图20是展示与不均匀高浓度的乳房的体积数据有关的直方图。
图21是图20的乳腺区域的提取处理内的直方图的曲线近似处理的补足图。
图22是展示与接收了边缘强调处理的体积数据有关的直方图的图。
图23是展示与乳腺密度高的状态的体积数据有关的直方图的图。
图24是展示在本实施方式中,用超声波扫描乳房区域的装置的概要的图。
图25是在图1的区域分割面设定处理单元12中,乳头和标识器(marker)定位处理的补足图。
图26是在图1的区域分割面设定处理单元12中,乳管的跟踪处理的补足图。
图27是在本实施方式中,展示相对摄像法等的分割方法的对应关系的图。
具体实施方式
以下参照附图说明本发明的实施方式。如果采用本实施方式,则本发明可以提供作为X射线诊断装置的X射线立体摄像装置、被装入X射线立体摄像装置内的图像显示装置、图像显示方法、用计算机实现图像显示部件的组合的程序、记录该程序的计算机可以读取的存储介质。
在本实施方式中采用X射线立体摄像法。X射线立体摄像法根据由多方向摄像的多个X射线透过图像再构成断面图像(为了和X射线CT的再构成区别,以下称为“成像”)。作为实际的处理,根据由多个方向拍摄的多个X射线透过图像,在指定的断面上的各点(象素)上选择通过该点的多条放射线,选择与这些放射线对应的象素值加算。通过该处理原理上指定断面上的部位鲜明映现,离开断面的位置的部位模糊。在这样的X射线立体摄像法中,在全摄像方向的X射线锥重叠的区域中,可以进行断面图像的成像(参照图4,图5)。
图1展示本发明的实施方式的X射线乳房摄像装置的结构。X射线乳房摄像装置具有X射线摄像台1。X射线摄像台1如图3所示,在C型臂17的一方的端部上装备X射线管球2。X射线管球2接收来自高电压发生器3的高电压的施加发生X射线。在C型臂17的另一方的端部上,在与X射线管球2相对的方向上安装平面检测器(还称为平板检测器)5。平面检测器5具有把入射X射线直接变换为电气信号的直接变换型,或者把入射X射线用荧光体变换为光,把该光变换为电气信号的间接变换型的多个半导体检测元件。多个半导体检测元件被排列成2维格子状。在多个半导体检测元件中伴随X射线入射发生的信号电荷经由数据收集单元6作为数字信号读出。平面检测器5可以代替图像增强器和光学性照相机的组合。
平面检测器5被收容在固定于C型臂17上的下侧压迫板18中。下侧压迫板18和被支撑在C型臂17上可以在箭头α方向移动的上侧压迫板19一同,被设置成用于压迫被检体的乳房。C型臂17由被安装在支柱20内的臂旋转机构4支撑为可以在箭头β方向上摇头旋转。例如,该旋转中心轴被设定在下侧压迫板18的附近,通过在和C型臂17的旋转相反的方向上旋转下侧压迫板18,如图4(a)所示,在多方向摄像时,放置被检体的乳房的下侧压迫板18停止,X射线管球2可以进行摇头运动。
本实施方式的X射线乳房摄像装置在具备X射线摄像台1的同时,还具备:系统控制单元7、摄像控制单元8、操作单元9、图像存储单元10、成像处理单元11、区域分割面设定处理单元12、投影处理单元13、断面变换处理单元14、画面构成单元15、显示单元16。摄像控制单元8通过根据事前设定的摄像条件控制高电压发生单元3、臂旋转机构4、数据收集单元6,执行把在压迫板18、19上被压成扁平状的被检体的乳房作为对象的多方向摄像(立体摄像)(参照图4(a)、图4(b))。对右侧的乳房和左侧的乳房各自进行立体摄像,把和摄像方向不同的右侧乳房有关的多个X射线透过图像的数据以及和摄像方向不同的左侧乳房有关的多个X射线透过图像的数据存储在图像存储单元10上。成像处理单元11如图6所示,对于用X射束以及一对压迫板18、19规定的大致正方形状的可成像区域(简称为摄像区域)内的等间隔的多个断面,成像分别对应的多个断面图像。把由多个断面图像组成的集合称为体积图像或者只称为体积。多个断面图像的数据,即体积图像的数据被存储在图像存储单元10中。
区域分割面设定处理单元12把摄像区域设定为多个,典型的是对于厚度方向(Y轴方向),设定用于分割成2个部分区域的区域分割面。区域分割面可以操纵操作单元9来手动地任意设定,而该初始位置由区域分割面设定处理单元12自动地设定。其初始的区域分割面从3种中选择。选择操作在事前或者处理过程中进行。
第1区域分割面如图7所示,是夹着乳房的上下2块压迫板18、19的中央面。具体地说,把下侧压迫板18的上面和上侧压迫板19的下面之间的距离(被压迫的乳房的厚度)设置为T,压迫板18、19平行,在从下侧压迫板18的上面到位于T/2的距离的面(中央面)上设定第1区域分割面。因为下侧压迫板18被固定在C型臂17上,所以如果上侧压迫板19的下面位置被确定,则可以设定中央面。上侧压迫板19的下面位置可以通过旋转编码器等的位置传感器来检测。而且通过图像处理可以确定上侧压迫板19的下面位置。上侧压迫板19的下面位置因为是被压迫的乳房和上侧压迫板19的边界,所以如图8所示,可以通过从其附近的多个断面中提取断面图像从软组织像变化到压迫板19的均匀像的断面来确定。
这时,如图12所示,如果在压迫板18、19的各个表面上把非常潜的细沟21形成为格子状,则从已生成的统计图像数据中,容易确定压迫板18、19的表面位置。因为知道设置的沟21的图形,所以通过图像匹配可以在2维图像中算出该图形包含什么程度的指标。对于成像图像的各分层图像,调查其指标,指标最高的分层是压迫板19的表面位置。如果在上下压迫板18、19上设置不同图像的沟,则可以独立确定上下压迫板表面。也可以不用沟21而在压迫板18、19上设置不透过X射线的标识器。
第2区域分割面如图9所示,被设置在从夹着乳房的上下2块压迫板18、19的下侧压迫板18的上面开始只离开规定的距离D的位置的面上。如上所述,因为选择压迫板18被固定在C型臂17上,所以不管上侧压迫板19的下面的位置如何,都始终被设定在一定位置上。
下侧或者上侧压迫板18、19的位置如以下那样确定。体积图像的各分层假设是和压迫板18、19大致平行。压迫板的方向因为由成像处理的条件确定,所以容易生成上述那样的分层。
把体积图像的分层层数设置为N,把最初的分层的位置设置为z1,把最后的分层位置设置为z2,把下侧或者上侧的压迫板18、19的表面设置为z3。z1、z2是作为成像时的条件指示的值。既有z3是装置设计上的常数的情况,如上所述,也有是用传感器测量的测量值的情况。此外,也有如上所述是根据图像确定的值的情况。
区域分割面例如在被指定为从压迫板18、19的表面开始在z轴正的方向上离开距离L的位置上的情况下,区域分割面的位置成为第a=1+N(z3+L-z1)/(z2-z1)层的分层位置。一般该值不是整数,表示在不超过a的最大的整数和比a大的最小的整数间存在区域分割面。
如图10、图11所示,从根据断面图像利用阈值处理以及图形匹配处理可以高精度地提取的、通过乳头且与压迫板18、19平行或者倾斜的多个面中,选择各面内象素值的积分值最高的断面,将该被选择的面设定为第3区域分割面。进而,代替通过乳头的面,也可以把通过摄像区域内的象素值的重心的面作为第3区域分割面设定。
进而,有关乳头提取处理的细节以后叙述。
摄像区域由区域分割面分割成上侧压迫板19侧的大致上半部分区域(A)和下侧压迫板18侧的大致下半部分区域(B)。投影处理单元13从大致上半部分区域内的多个断面图像中,通过最大值投影处理(MIP),生成与大致上半部分区域对应的投影图像。此外,投影处理单元13从大致下半部分区域内的多个断面图像中通过最大值投影处理,生成与大致下半部分区域对应的投影图像。MIP处理如公知的那样,作为投影面在此设定与压迫板18、19平行的面,把与该投影面垂直的轴(放射线)上的象素值的最大值作为该放射线和投影面的交点象素的值。即,对于被压迫的乳房的厚度方向对每个部分区域独立实施投影处理。作为投影处理,并不限于MIP,也可以是积分放射线上的象素值的处理。投影处理单元13在对左右乳房各自生成大致上半部分区域的部分投影图像和大致下半部分区域的部分投影图像的同时,对于左右乳房的各自生成全部投影区域的投影图像。把左右乳房各自的大致上半部分区域的部分投影图像,和左右乳房各自的大致下半部分区域的部分投影图像,和左右乳房各自的全摄像区域的投影图像的数据存储在存储单元10中。
断面变换单元14根据摄像区域内的多个断面图像,在和该断面正交的同时相对左右乳房的各自生成通过该断面的中央的图6所示的纵断面的断面图像。右乳房的纵断面图像和左乳房的纵断面图像的数据被存储在存储单元10中。
画面构成单元15构成显示单元17的显示画面。如图2所示,夹着用虚线表示的画面的中央线,把有关右乳房的大致上半部分的部分投影图像,和有关左乳房的大致下半部分的部分投影图像配置成对称。在大致上半部分的部分投影图像的下侧上,把有关右乳房的大致下半部分的部分投影图像和有关左乳房的大致下半部分的部分投影图像夹着画面中央线配置成对称。在大致上半部分的部分投影图像的外侧上把有关右乳房的全区域的投影图像和有关左乳房的全区域的投影图像夹着画面中央线对称配置。进而,有关左右乳房的全区域的投影图像,与操作单元9的适当的操作相应,可以切换到通过断面变换处理单元14生成的与区域分割面有关的断面图像。在大致下半部分的部分投影图像的外侧上,把有关右乳房的纵断面图像和有关左乳房的纵断面图像夹着画面中央线大致对称配置。在各纵断面图像中,分别与区域分割面对应的线光标彼此重合。线光标的显示位置与操作单元9的适当的操作相应可以上下移动。在线光标移动时,跟随该移动更新部分投影图像。
进而,在图27中,展示乳房构造、与之对应的摄像方法以及显示方法。在图27中,“外侧”、“内侧”是与乳房的构造对应的记号。根据分割面把体积区域分割为2个,生成各自区域的图像,而此时,把生成分割面的表面侧区域的方法作为方法A,把生成分割面的背面侧区域的方法作为方法B。与方法A、方法B对应的各个记号可以用图27那样的表来确定。
如上所述,在本实施方式中,制成分割了压迫乳房的体积图像(多层断面图像)的图像,显示表面和背面的图像。作为分割的区域分割面的确定方法,使用把压迫板的中央断面作为区域分割面的方法,和通过图像处理提取脂肪层作为区域分割面的方法等。在通常的X射线乳腺摄影法中因为多个乳腺重叠,所以有可能错误地诊断为肿瘤,当在本发明中通过制成乳房分割图像,可以提供没有乳腺重叠的图像,诊断精度提高。
进而,上述的第3区域分割面是通过乳头的面。以下说明区域分割面设定处理单元12的乳头提取处理。
如图13所示,把体积数据的3个轴设置为X,Y,Z。体积数据假设在X,Y,Z方向上分别具有n、m、l[象素]的大小。首先,从体积数据中通过滤波处理除去散布的图像干扰(斑点噪声等)。以下,被滤波后的体积数据通过阈值处理被2值化。阈值为了区别空气区域,和此外的以乳房以及乳头为主的区域被设定在适宜的值。对于具有超过阈值的值的体素(voxel),即主要包含乳房以及乳头的区域内的体素例如给出“1”,对于具有阈值以下的值的体素,即空气区域内的体素例如给出“0”。
在此,假设向量f是表示乳头存在的概率高的方向的向量。计算该向量f和给予了“1”的体素各自的内积。从给予了“1”的多个体素中,选择内积最大的体素。把选择出的体素的中心坐标值确定为乳头位置。
作为其他的乳头提取方法有半自动方法。首先,在根据体积数据生成的投影显示、MIP显示(最大值投影)等的2维图像上操作者手动指定乳头位置。在该指定阶段中,只是确定了在2维画面上的位置,进而,需要确定进深方向(投影方向)的位置。
操作者沿着通过了在2维图像上指定的点的投影直线(典型的是与画面正交的方向)检索象素值,用以下的(a)、(b)之一的方法确定乳头位置。
(a)从投影直线上的多个体素中选择具有最大象素值的体素,把该选择出的体素的位置作为乳头位置。
(b)从投影直线上的多个体素中,提取具有超过规定阈值的象素值的体素。把提取的体素连续的区域(体素序列)的中央位置设置成乳头位置。当体素序列存在多个的情况下,选择最长的体素序列。进而,阈值例如可以根据MIP的显示条件(灰电平和窗口宽度)来确定。
进而在其他的乳头提取方法中,在从体积数据由MPR处理(断面变换处理:multiplanar reconstruction processing)生成的任意断面上,操作者指定乳头位置,把指定的位置作为乳头位置。
进而在其他的乳头提取方法中,在由体积构图(VolumeRendering)显示的图像上,操作者指定乳头位置。从该位置开始,如以下那样确定乳头的3维坐标。沿着在体积构图图像中通过点击的点的投影直线,检索在体积构图制作中使用的不透明度。如果该不透明度大于等于第1阈值,则累计不透明度,把该累计值达到第2阈值的点作为乳头位置。
进而其他的乳头提取方法是对于MPR图像显示的各断面实施图像处理,自动地算出乳头位置的图像处理方法。在该图像处理方法中主要包含2种图像处理方法。第1种方法是铺上判别乳头形状的2维掩模,进行图形匹配(Pattern Matching Algorithm)的方法。第2图像处理方法如下。首先,从体积数据中,通过滤波处理除去散布的图像噪声(斑点噪声等)。以下,根据象素值2值化经滤波处理的体积数据。用于2值化的阈值如上所述和用于分离空气区域和此外的乳房等的区域用的阈值相同。对于被2值化的体积数据实施一般性的一次微分、2次微分等的滤波处理,提取轮廓。整理被提取的轮廓线。
如图14示例的那样对轮廓线上的多点各自求切线的倾斜,算出倾斜是零的位置,即典型的是使用与被检体的身体轴大致平行的切线B的接点位置算出乳头位置。或者以微分值最大的位置作为乳头位置。
进而在其他的乳头提取方法中,首先,对乳头的体积数据适用高频除去滤波器处理,除去包含在体积数据中的空间性的高频成分。通过阈值处理把体积数据2值化为空气区域和此外的包含乳房和乳头的区域。如果需要,则作为补缺处理,对包含空气区域以外的乳房以及乳头的区域(“1”区域),在小距离的膨胀处理(dilation)后,进行同样距离的收缩处理(erotion)。
只以距离n1对接收了补缺处理的体积数据A的“1”区域进行收缩处理,对于该结果只以距离n2进行膨胀处理。此时,n2假设大于等于n1。把接收了该处理的体积数据称为B。该运算有平滑化“1”区域的表面的作用,其结果,体积数据B与除去了突起部分的体积数据A对应。通过从具有突起部分的体积数据A中减去除去了突起部分的体积数据B,生成残留多个突起部分区域的体积数据C。体积数据C内的多个突起部分区域之一是乳头区域。
对体积数据C内的联结区域进行标记(称为:connectedcomponent labeling)。求各个区域的体积和重心坐标值。预先以对每种摄像方法(压迫乳房的标准的几种方法)确定的选择条件为基础选择1个区域。例如,选择体积在V1以上V2以下的区域中重心坐标值的向量f的内积最大的区域。向量f是表示乳头存在的概率高的方向的向量。求选择出的区域的重心,把它作为乳头位置。
以下,在分割面设定处理单元12中,还具有根据乳腺区域的分布变化确定分割面的功能。参照图5,把分离2个乳腺区域的中央面作为分割面。首先,说明从体积数据中提取乳腺的方法。乳腺的区域在X射线断层图像上,是比乳房的其他许多构造更大的象素值。因而,如果确定某一阈值,则可以通过2值化处理大致求得乳腺的区域。阈值的确定方法以后叙述。在摄像区域内设定与压迫板18、19大致水平的多个分层面。换句话说,在被压迫的乳房厚度方向上以一定的间隔设置分层面(XZ面)。求对各分层面的每面提取的乳腺区域的断面积。如果对每个分层位置求断面积,则如图15(a)所示,可以得到表示与厚度方向有关的乳腺区域断面积的变化的曲线图。在该曲线图上,求具有比某一定面积(虚线)还大的区域的范围(从T点到B点的范围),把该范围的中点M1作为确定的分层位置。
作为其他方法,如图15(b)所示,求从中点M1到下方断面积出现峰值的分层位置P1。同样,从M1开始求上方的峰值,确定该分层位置P2。设置为P2。把P1和P2的中点作为M2。
把这样求得的分层位置(M1或者M2)作为分割面,在本例子中,也可以确定T、P1、M1、P2、B这多个分层,把它作为分割面制成分割在多个区域上的投影像来显示。
进而,阈值可以从操作者设定的MIP和MPR等的图像显示的显示条件中求得。这种情况下,希望构成为在MIP和MPR等的图像中,(用例如红色等)重叠显示在已确定的阈值中提取的阈值,当显示条件变更的情况下,根据该操作者操作,再次确定阈值,重新提取区域,更新重叠图像。此外,已知有最佳阈值、大津法和百分比拼贴法等确定阈值的几种自动的方法,这些方法都可以使用。
下面说明处理单元12中的阈值的另一确定方法,处理单元12首先从体积数据中制成象素值的直方图(参照图16)。根据直方图判定乳房以外的区域,把该区域作为α。根据直方图把乳房内的区域设置为β、γ。在此把图16的黑色部分的值设置为φ,把白色部分的值设置为ψ(φ<ψ)。
区域α:相当于乳房的脂肪层。
区域β:不与乳房的主要构造对应。
区域γ:相当于乳房的乳腺区域等。
背景(空气)区域的象素值比φ还小,在图16的直方图中未展示。
在区域α以及区域β上和区域γ上设定不同的值,生成2值图像,提取乳腺构造。
以下,确定阈值。乳房一般有3种乳腺构造状态。状态依赖于年龄、哺乳经历等。
状态a:乳腺层全部被置换为脂肪层的乳房(参照图17(a))
状态b:不均匀高浓度乳房的状态(参照图17(b))
状态c:乳腺密度高的乳房(参照图17(c))
在此,表示乳房具有3个状态。图18示例乳腺层被全部置换为脂肪层的乳房时的直方图。在乳腺层全部被置换为脂肪层的乳房中,乳腺构造区域γ从直方图中容易分成乳腺区域和非乳腺区域这2个区域。
如果把象素值设置为X,把象素总数设置为Y,则在用近似曲线f(X)表现直方图时可以用图19表现。
操作者指定区域α和β的边界位置φ,以及区域β和γ的边界位置ψ。假设把Y=f(X)的微分值从负变为正的X设置为X1(n是整数)。把X1作为中心,把体素宽度η的ε%设置为区域β。进而,有ε由操作者设定的情况,和使用模型的直方图设定默认值的情况。
在不均匀高浓度的乳房时的直方图示例于图20。如果把象素值设置为X,把象素总数设置为Y,则在用近似曲线f(X)表现直方图时可以用图21表示。但是,看不出显著的变化。
因而,如果经由傅立叶变换把图像变换为频率区域的信息,在频率空间上降低到低频下,则边缘被强调。在反傅立叶变换中返回实际空间。在频率处理后的直方图中,如图22所示,边缘被强调。
如果把象素值设置为X,把象素总数设置为Y,用近似曲线f(X)近似直方图,则可以得到边缘被强调的直方图。在边缘被强调的直方图上,操作者手动指定区域α和β的边界位置φ、以及ψ。
在乳腺密度高的乳房中,如图23示例的那样,几乎都是乳腺区域。操作者指定区域α和γ的边界位置φ。
对于各乳房,也有预先用默认值确定象素值的阈值,选择乳房的种类,设定直方图的阈值的方法。操作者选择被提示的图像是乳腺层全部被置换为脂肪层的乳房、不均匀高浓度的乳房、乳腺密度高的乳房的哪个。预先确定与各个乳房图像对应的阈值,使用该阈值进行阈值处理。
此外,有自动地确定乳腺层全部被置换为脂肪层的乳房、不均匀高浓度的乳房、乳腺密度高的乳房的哪个的方法和选择后的阈值的设定方法。首先被提示的图像把乳腺层全部被置换为脂肪层的乳房、不均匀高浓度的乳房、乳腺密度高的乳房的哪个和预先模型化的直方图取相关,从相关值中选择是哪类的乳房。
事前确定各个乳房图像的阈值。把相对被选择出的乳房的种类的阈值作为默认值(参考值)进行阈值处理。
以下,说明从体积数据中设定分割面的方法。本方法如图24所示,从由超声波扫描求得的乳房的体积数据中,半自动地设定分割面。首先,摄像者指定通过乳头或者想分割显示的断面(可以是多个)。用超声波3D扫描包含乳房的区域来收集体积数据。把设定的平面作为分割面进行图像显示。
此外,如以下所示,也可以根据体积数据设定分割面。是一种如知道乳头的位置和图像的位置关系那样使用可以定位的卡具,收集包含乳房的区域的体积数据,机械地设定断面的方法。如图25所示,可以使标识器和乳头的一方与另一方对齐。
收集图像设定通过乳头的分割面。例如,使用距离压迫板18、19的距离设定分割面。把乳头(标识器)作为旋转中心,使平面旋转,操作者设定分割面。
下面说明从体积数据中提取乳腺构造的另一方法。本方法是操作者选择跟踪开始点,从该开始点使用一般的跟踪方法提取乳腺构造的方法。2维图像显示(投影,MIP,MPR)体积数据。如图26所示,操作者点击乳管的位置。使用公知的血管跟踪方法跟踪乳管。
进而,本发明并不限于上述实施方式本身,在实施阶段在不脱离其主旨的范围中可以变形构成要素而具体化。此外,通过在上述实施方式中公开的多个构成要素的适宜组合,可以形成各种发明。例如,可以从本实施方式展示的全部构成要素中删除几个构成要素。进而,也可以适宜组合不同的实施方式的构成要素。
Claims (1)
1.一种图像显示装置,用于生成乳房图像,包括:
算出区域分割位置的部件;
生成用上述区域分割位置分割的分割区域的图像的部件;
显示上述已生成的图像的部件,
算出上述区域分割位置的部件以上述图像为基础,检测上述乳房的乳头的位置信息,根据上述被检测出的乳头的位置信息算出上述区域分割位置。
Applications Claiming Priority (6)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2004226942 | 2004-08-03 | ||
JP2004226942 | 2004-08-03 | ||
JP2004-226942 | 2004-08-03 | ||
JP2005-216587 | 2005-07-27 | ||
JP2005216587 | 2005-07-27 | ||
JP2005216587A JP4709604B2 (ja) | 2004-08-03 | 2005-07-27 | 画像表示装置、画像表示方法、記憶媒体及びプログラム |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200510091028A Division CN100594843C (zh) | 2004-08-03 | 2005-08-03 | 图像显示装置、图像显示方法、存储介质以及程序 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101612047A true CN101612047A (zh) | 2009-12-30 |
CN101612047B CN101612047B (zh) | 2013-03-06 |
Family
ID=36092084
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 200910160229 Active CN101612047B (zh) | 2004-08-03 | 2005-08-03 | 图像显示装置 |
CN200510091028A Active CN100594843C (zh) | 2004-08-03 | 2005-08-03 | 图像显示装置、图像显示方法、存储介质以及程序 |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200510091028A Active CN100594843C (zh) | 2004-08-03 | 2005-08-03 | 图像显示装置、图像显示方法、存储介质以及程序 |
Country Status (1)
Country | Link |
---|---|
CN (2) | CN101612047B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TWI511085B (zh) * | 2013-08-08 | 2015-12-01 | Univ Chien Hsin Sci & Tech | Nipple image segmentation method |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5528518B2 (ja) * | 2012-09-28 | 2014-06-25 | 富士フイルム株式会社 | 放射線画像生成装置および方法 |
CN105228525B (zh) * | 2013-03-29 | 2019-10-18 | 通用电气公司 | 乳房造影装置以及使用该装置对待检目标照相的方法 |
CN104127202B (zh) * | 2013-05-03 | 2016-12-28 | 上海西门子医疗器械有限公司 | 压迫器的定位方法和定位装置、以及胃肠x光机 |
CN104207797B (zh) * | 2013-06-03 | 2018-06-19 | 上海联影医疗科技有限公司 | 基于ct定位片确定病人尺寸、位置及x射线经过病人衰减的方法 |
CN104182965B (zh) * | 2014-01-22 | 2016-03-30 | 上海联影医疗科技有限公司 | 一种乳腺图像中分割胸肌的方法 |
WO2017054775A1 (en) * | 2015-09-30 | 2017-04-06 | Shanghai United Imaging Healthcare Co., Ltd. | System and method for determining a breast region in a medical image |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FI80996C (fi) * | 1988-05-26 | 1991-10-25 | Automed Oy | Mammografiamenetelmä ja -laite |
US5579360A (en) * | 1994-12-30 | 1996-11-26 | Philips Electronics North America Corporation | Mass detection by computer using digital mammograms of the same breast taken from different viewing directions |
US5572565A (en) * | 1994-12-30 | 1996-11-05 | Philips Electronics North America Corporation | Automatic segmentation, skinline and nipple detection in digital mammograms |
US5640956A (en) * | 1995-06-07 | 1997-06-24 | Neovision Corporation | Methods and apparatus for correlating ultrasonic image data and radiographic image data |
-
2005
- 2005-08-03 CN CN 200910160229 patent/CN101612047B/zh active Active
- 2005-08-03 CN CN200510091028A patent/CN100594843C/zh active Active
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TWI511085B (zh) * | 2013-08-08 | 2015-12-01 | Univ Chien Hsin Sci & Tech | Nipple image segmentation method |
Also Published As
Publication number | Publication date |
---|---|
CN101612047B (zh) | 2013-03-06 |
CN100594843C (zh) | 2010-03-24 |
CN1739452A (zh) | 2006-03-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4709604B2 (ja) | 画像表示装置、画像表示方法、記憶媒体及びプログラム | |
CN105705096B (zh) | 用于显示医学图像的方法和设备 | |
KR102109588B1 (ko) | 유방 조직 이미지를 프로세싱하고, 디스플레잉하고, 네비게이팅하기 위한 방법 | |
JP6392309B2 (ja) | 自動合焦を含む、トモシンセシススタックをナビゲートするためのシステム | |
US8165365B2 (en) | Method and system for displaying tomosynthesis images | |
US9401019B2 (en) | Imaging tomosynthesis system, in particular mammography system | |
US6487432B2 (en) | Method and system for selecting and displaying medical image data | |
EP2212859B1 (en) | Method and apparatus for volume rendering of data sets | |
CN107405126B (zh) | 检索成对的医学图像的对应结构 | |
US8184890B2 (en) | Computer-aided diagnosis and visualization of tomosynthesis mammography data | |
CN102763137B (zh) | 用于由层析数据集生成2d图像的系统和方法 | |
CN100594843C (zh) | 图像显示装置、图像显示方法、存储介质以及程序 | |
US9842415B2 (en) | Method for processing tomosynthesis acquisitions in order to obtain a representation of the contents of an organ | |
US20120157819A1 (en) | Imaging method and imaging device for displaying decompressed views of a tissue region | |
US8611492B2 (en) | Imaging method for rotating a tissue region | |
CN101443815A (zh) | 重建图像的方法和装置 | |
CN101005803A (zh) | 灵活融合3dra-ct的方法 | |
US10973485B1 (en) | Enhanced volume viewing | |
US11158045B2 (en) | Method and apparatus for performing 3D imaging examinations of a structure under differing configurations and analyzing morphologic changes | |
KR20160117807A (ko) | 유방 병변 진단을 위한 영상 판독 보조 방법 및 장치 | |
Kita et al. | Correspondence between different view breast X-rays using a simulation of breast deformation | |
CN100583161C (zh) | 用于显示以立体数据组成像的对象的方法 | |
EP3503020A1 (en) | Clinical task-based processing of images | |
JPH07236634A (ja) | 診断支援装置 | |
Costa | Automated Deformable Registration of Breast Images: towards a software-assisted multimodal breast image reading |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C41 | Transfer of patent application or patent right or utility model | ||
TR01 | Transfer of patent right |
Effective date of registration: 20160720 Address after: Japan Tochigi Patentee after: Toshiba Medical System Co., Ltd. Address before: Tokyo, Japan, Japan Patentee before: Toshiba Corp Patentee before: Toshiba Medical System Co., Ltd. |