专利名称：用于形成3d超声图像的设备和方法
技术领城本发明一般涉及一种三维(3D)超声诊断系统，更具体地说涉及一种用于从目标对象的二维(2D)超声图像中自动检测出轮廓并利用该轮廓内的体积数据形成3D超声图像的设备和方法。
背景技术：
三维(3D)超声诊断系统是用于提供诸如空间信息、解剖信息等临床信息，这些临床信息不能从传统的二维图像中提供。3D超声诊断系统从通过探针从目标对象接收的信号中获得体积数据，并且针对所获得的体积数据执行扫描转换程序。通过根据从扫描转换后的数据中获得的图像执行再现程序而在例如监视器、屏幕等显示设备上显示目标对象的3D超声图像。因此用户可以获得目标对象的临床信息。
如在现有技术中众所周知的，探针一般具有多个传感器，其中分别将脉冲信号输入每个传感器的定时被适当地延迟。因此聚焦超声束沿着发射扫描线被发射到目标对象中。每个传感器在不同的接收时间接收从发射扫描线上的焦点反射的回声信号。接收信号被发射给束成形器。接收信号被适当地延迟，其中在波束成形器中对延迟的接收信号求和。因此，从发射扫描线上的焦点反射的表示能级的接收聚焦束被输出。以上程序被反复执行，直到产生由多个扫描线的接收聚焦束形成的目标对象的2D片层图像。
体积数据获取单元通过综合2D超声图像而输出体积数据，该2D超声图像表示从束成形器输入的目标对象的截面。体积数据从由存在于3D空间内并用圆环坐标系定义的目标对象所反射的信号中产生。因此，为了在具有笛卡尔坐标系的显示设备(例如监视器、平面等)中执行体积数据的再现程序，需要用于执行体积数据的坐标转换的扫描转换。扫描转换在扫描转换器中实现。
扫描转换器中的扫描转换后的体积数据通过典型的容积再现程序被这样再现，使得3D超声图像得以显示。用户通过显示在显示设备上的3D超声图像获得目标对象的临床信息。
3D超声诊断系统主要被用于在产科和妇科领域中用3D超声图像显示胎儿的形状。在通过扫描孕妇的腹部区域而获取了体积数据之后，根据所获取的体积数据执行容积再现程序。因此可以用3D超声图像显示胎儿的形状。然而由于体积数据包括子宫组织、脂肪组织、羊水、浮动物质和胎儿的混合数据，如果再现程序被直接应用于体积数据，则很难用3D超声图像清楚地显示胎儿的形状。因此，为了用3D超声图像显示胎儿的形状，需要从例如羊水等相邻区域中分割出胎儿区域。
因此，通过使用与3D超声诊断系统连接的外部接口设备(例如鼠标、键盘等)，可以产生包围显示设备上显示的2D超声图像中胎儿形状的感兴趣区域(ROI)方框，如图1A所示。此后，如图1B所示，通过细致地操作外部接口产生最后的ROI方框。存在于从该ROI方框内的图像中检测出的轮廓内的体积数据被这样再现，使得胎儿的3D超声图像能够被显示。
然而，由于ROI方框的产生以及目标对象图像的轮廓的检测都是由用户在3D超声诊断系统中手动地进行操作，所以最后显示的3D超声图像的质量取决于用户的专门技能。换句话说，根据产生ROI方框的用户，ROI方框的尺寸是不一致的。同样地，由于不能精确地显示所期望的目标对象的3D超声图像，经常存在问题。
而且即使用户是专家，也存在问题，即产生ROI方框以及从2D超声图像中检测出目标对象图像的轮廓需花费很长的工作时间。这是因为用户直接在2D超声图像上产生ROI方框。此外，如果对于所期望的目标对象的3D超声图像而言ROI方框的尺寸不精确，则存在以下问题，即在体积数据再现程序或ROI方框中的目标对象的轮廓检测程序中可能会产生误差。
发明内容
因此，本发明的一个目的是提供用于在减小可能在体积数据再现过程和目标对象的轮廓检测过程中产生的误差的同时形成目标对象的精确的3D超声图像的设备和方法。本发明的另一个目的是通过自动地产生ROI方框并且检测存在于ROI方框内的目标对象的轮廓来减少在产生感兴趣(ROI)方框的过程以及检测目标对象的轮廓的过程中耗费的时间。
根据本发明的一个方面，提供一种用于形成三维(3D)超声图像的设备，包括第一单元，用于在2D超声图像上产生感兴趣区域(ROI)方框；第二单元，用于检测ROI方框中的目标对象的轮廓；以及第三单元，用于通过再现存在于检测到的轮廓中的体积数据来形成3D超声图像。
根据本发明的另一方面，提供一种用于形成三维(3D)超声图像的方法，包括以下步骤a)在2D超声图像上产生感兴趣区域(ROI)方框；b)检测ROI方框中的目标对象的轮廓；以及c)通过再现存在于检测到的轮廓中的体积数据来形成3D超声图像。
根据以下结合附图给出的优选实施例的说明，本发明的上述以及其他目的和特征将变得显而易见，其中图1A和1B是描述在2D超声图像上产生的感兴趣区域(ROI)方框的图；图2是说明用于形成3D超声图像的设备的示意性框图，该设备根据本发明的优选实施例构成；图3是说明图2的ROI方框产生单元的详细的框图；图4A至4C是说明对2D超声图像进行图像分割的结果的实例；图5A至5F示出了自动调节用于胎儿面部图像的ROI方框的尺寸的过程；图6A示出了用于确定目标对象图像的轮廓检测能力的掩码；图6B是说明应用掩码的2D超声图像的图；图7A和7B示出了对ROI方框中的目标对象图像执行平滑和亮度对比度的实例；图8A和8B示出了对ROI方框中的目标对象图像执行二值化的实例；以及图9A和9B示出了检测胎儿图像的上表面上的轮廓的实例。
具体实施例方式
图2是说明3D超声诊断设备中的三维(3D)超声图像形成设备3的示意性框图，该3D超声图像形成设备3根据本发明的优选实施例构成。该3D超声图像形成设备3包括感兴趣区域(ROI)方框产生单元31、第一确定单元33、第二确定单元35、轮廓检测单元37和3D图像处理单元39。当安装在该3D超声图像形成设备3中的控制面板(未示出)的3D按钮被激活时，该3D超声图像形成设备3开始工作。
如果用户激活了该3D按钮，通过3D超声诊断设备中的探针和束成形器获取的2D超声图像中的一个则以亮度模式(B-mode)显示在显示设备(未示出)上。根据本发明，显示在显示设备上的2D超声图像是表示目标对象的一个中心片层的2D超声图像。B-mode表示从目标对象反射的信号能量以亮度水平进行显示。ROI方框产生单元31在显示在显示设备上的2D超声图像上自动产生ROI方框。
第一确定单元33确定自动产生的ROI方框的尺寸是否适合目标对象的2D超声图像的尺寸。第二确定单元35确定存在于ROI方框内的目标对象的轮廓是否能够被检测。轮廓检测单元37检测存在于ROI方框内的目标对象的轮廓。3D图像处理单元39在存储于体积数据获取单元内的体积数据中选择存在于目标对象轮廓中的体积数据，并且通过再现这些体积数据来形成3D超声图像。
在下文中将通过参考图3至5F来详细讲述用于在2D超声图像上自动产生ROI方框的ROI方框产生单元31，该2D超声图像显示在3D超声诊断设备的显示设备上。
图3是说明图2的ROI方框产生单元31的详细的框图。如图3所示，ROI方框产生单元31包括图像分割单元311、ROI方框设置单元313和ROI方框调节单元315。图像分割单元311将2D超声图像分割成目标对象图像和与该目标对象相邻的背景图像。ROI方框设置单元313在2D超声图像上设置具有预定尺寸的ROI方框，而ROI方框调节单元315调节该ROI方框的尺寸。
由于在用于形成目标对象的体积数据的2D超声图像中存在一些使清楚地显示3D超声图像变得困难的因素(例如斑点噪声等)，所以应首先执行用于消除这种因素的程序。图像处理单元311通过滤波程序来消除存在于2D超声图像中的斑点噪声。根据本发明的优选实施例，通过使用一个Lee滤波器来执行滤波程序，以便只消除斑点噪声，同时保留2D超声图像中目标对象的边缘信息。
图像分割单元311为已消除斑点噪声的2D超声图像的二值化设置一个阈值。通过参考该阈值这样对2D超声图像进行二值化，以致该2D超声图像能够被分割成两类。由于该2D超声图像以B-mode显示，所以根据该2D超声图像的亮度来设置该阈值。例如，如果具有像素值0至255的2D超声图像的阈值为“t”，则该2D超声图像被分割成具有像素值{0，1，2…，t}的第一类以及具有像素值{t，t+1，t+2，…，255}的第二类。
图像分割单元311根据第一二值化阈值t1对2D超声图像执行第一分割程序。因此，该2D超声图像可以被分割成目标对象图像和背景图像。通过参考第二阈值t2执行针对通过第一分割程序被分割出来的背景图像的第二分割程序。因此该背景图像可以被分割成目标对象图像和背景图像。此后，通过参考第三阈值t3执行针对通过第二分割程序被分割出来的背景图像的第三分割程序，以致于背景图像最后被分割成目标对象图像和背景图像。由于反复执行这些分割程序，所以图像分割单元311将2D超声图像更清楚地分割成目标对象图像和与该目标对象图像相邻的背景图像。在图像分割单元311中产生的分割图像的值被输出到ROI方框调节单元315。
图4A示出了原始的2D超声图像。图4B示出了执行了一次图像分割程序的2D超声图像。图4C示出了执行了两次图像分割程序的2D超声图像。
ROI方框设置单元313设置配置有四个边界的ROI方框的尺寸，并且将所设置的ROI方框输出到ROI方框调节单元315，该边界在最初显示在显示设备上的2D超声图像上产生。一般根据需要以3D超声图像显示的目标对象的类型来预先设置ROI方框的尺寸。
ROI方框调节单元315调节由ROI方框设置单元313设置的ROI方框的尺寸，以适合存在于ROI方框内的目标对象的尺寸。在这段时间期间，目标对象图像的特征对于调节ROI方框的尺寸而言变成一个重要的因素。胎儿是目标对象的一个例子。胎儿的超声图像在头部和身体之间的边界处具有一个凹部，其中头部的前表面比身体的前表面更弯曲，该头部的前表面为面部的表面。还存在面部比身体长的特征。
参考图5来描述在ROI方框调节单元315中将针对整个胎儿的ROI方框的尺寸自动调节成针对面部的ROI方框的尺寸的一个实例。
首先，如图4C所示，ROI方框调节单元315从经图像分割后的2D超声图像中选择存在于由ROI方框设置单元313设置的ROI方框内的图像。如图5A所示，为所选择的胎儿图像设置一个二值化阈值，其中通过参考该阈值而将二值化处理应用于所选择的胎儿图像。接着，ROI方框调节单元315从胎儿图像的二值化区域中去除噪声区域。因此如图5B所示检测到利用一部分胎儿图像确定的二值化区域。更具体地说，ROI方框调节单元315检查每个二值化区域的亮度并且计算亮度平均值。当二值化区域的亮度低于亮度平均值时，该二值化区域则被认为是噪声区域并因此被去除。而且即使亮度大于亮度平均值，噪声区域也可能存在。因此，ROI方框调节单元315在二值化区域上设置一个阈值。同样地，当二值化区域的一个像素值低于该阈值时，该二值化区域则被认为是噪声区域并因此被去除。
在从二值化区域中去除了噪声区域之后，ROI方框调节单元315为全部像素分配像素值“255”，该全部像素存在于从位于所检测到的二值化区域的上表面上的像素到位于由ROI方框设置单元313所设置的ROI方框下边界上的像素中。因此，可以产生如图5C所示的胎儿的二值化图像。接着，ROI方框调节单元315如图5D所示在该二值化图像上产生表面曲线。然后，ROI方框调节单元315搜索该表面曲线上彼此相连的像素中相对最大的点和相对最小的点，并且指定对应于胎儿头部和身体之间凹部的相对最小的点。
如图5D所示，在对应于胎儿头部和身体之间凹部的相对最小的点被指定为N2之后，ROI方框调节单元315指定表面曲线的起始点N1和端点N3。因此，该表面曲线被分割成两个区域[N1，N2]和[N2，N3]。ROI方框调节单元315计算在第一表面曲线区域[N1，N2]和第二表面曲线区域[N2，N3]上存在的相对最大的点的数量。此外，具有较大数量的相对最大的点的区域被确定为胎儿的面部表面。由于胎儿面部表面弯曲度大比胎儿身体表面弯曲度，所以胎儿面部上的相对最大的点的数量大于胎儿身体上的相对最大的点的数量。
在确定了胎儿的面部区域之后，ROI方框调节单元315将ROI方框的左边界移动到表面曲线的相对最大的点N2。因此，ROI方框调节单元315自动地调节由ROI方框设置单元313所设置的ROI方框的左/右边界的位置以适应胎儿图像的面部大小。图5E示出了一个实例，其中设置在显示于显示设备上的2D超声图像上的ROI方框的左/右边界由ROI方框调节单元315自动地进行调节。
在为胎儿面部图像调节了ROI方框的左/右边界之后，ROI方框调节单元315为存在于调节后的ROI方框内的胎儿面部图像执行二值化、噪声区域去除、二值化图像产生和表面曲线产生程序。为了从胎儿图像中产生二值化图像，R0I方框调节单元315给由从胎儿面部图像中检测到的二值化区域的最左边的表面到调节后的ROI方框的右边界而存在的全部像素分配像素值“255”。
ROI方框调节单元315将ROI方框的上边界朝着ROI方框的下边界移动，直到上边界与表面曲线接触。然后其将上边界与表面曲线接触的位置确定为ROI方框的上边界位置。同样地，ROI方框调节单元315将ROI方框的下边界朝着ROI方框的上边界移动，直到下边界到达表面曲线。然后其将ROI方框的下边界与表面曲线接触的位置确定为ROI方框的下边界位置。图5F示出了一个实例，其中在2D超声图像上所设置的ROI方框最后由ROI方框调节单元315针对胎儿面部图像而调节。
为了简便起见，虽然针对胎儿图像描述了ROI方框产生单元31中的ROI方框调节程序，但是显而易见的是，通过根据上述ROI方框调节程序自动调节ROI方框的边界，可以调节在2D超声图像的任意目标对象上所设置的ROI方框的尺寸以适合任意目标对象的尺寸。
图2所示的第一确定单元33确定在ROI方框产生单元31中产生的ROI方框的尺寸对目标对象图像而言的适宜性。一般对相同类型的目标对象图像而言从ROI方框产生单元31中产生的ROI方框具有相似的尺寸。因此，第一确定单元33从构造在3D超声诊断系统中的存储器(没有示出)中选择标准直方图中的一个，该标准直方图预先为每种类型的目标对象图像所设置。然后计算为目标对象图像所选择的一个直方图和多个直方图的均方误差值，该目标对象图像存在于从ROI方框产生单元31输入的ROI方框内。
结果，如果所计算的均方误差值大于以前预定的值，第一确定单元33则确定从ROI方框产生单元31中产生的ROI方框的尺寸是不适当的，于是停止3D超声诊断系统的操作。此后，其通知用户所产生的ROI方框的尺寸是不适当的。另一方面，如果所计算的均方误差值小于预定的值，第一确定单元33则确定从ROI方框产生单元31中产生的ROI方框的尺寸是适当的。然后，其将从ROI方框产生单元31输入的2D超声图像输出给第二确定单元35。从ROI方框产生单元31输入的ROI方框中的目标对象图像的直方图被发送给存储器以更新目标对象图像的标准直方图。
图2所示的第二确定单元35确定是否能检测到ROI方框中目标对象图像的轮廓。也即，当第二确定单元35确定ROI方框中目标对象图像的轮廓检测的可能性时，其提高用于在3D超声诊断系统中显示3D超声图像的功率。因为超声的吸收和反射根据目标对象(例如胎儿以及包围胎儿的羊水)的每部分而不用，所以由于清楚地表示亮度差别的边缘效应而在2D超声图像的目标对象中形成边界。第二确定单元35基于边缘效应确定目标对象的轮廓检测的可能性。
图6A示出了用于确定是否能检测2D超声图像的轮廓的掩码。图6B示出了应用图6A的掩码的超声图像。通常使用配置有3×3的像素的掩码以检测2D超声图像中目标对象的边界。在设计2D超声图像的掩码时，(1)对应于黑暗部分的像素值被设置为“-1”，(2)对应于明亮部分的像素值被设置为“1”，以及(3)形成黑暗部分和明亮部分之间的边界的像素值被设置为“0”。根据本发明的优选实施例，使用如图6A所示的、配置有扩展了形成边界的像素区域的5×5像素的掩码。构成掩码的像素的大小与构成2D超声图像的像素的大小相同。
第二确定单元35通过使位于掩码(图6A所示)中心的像素与在与从ROI方框的上边界到胎儿面部图像的一半的一部分对应的区域“A”中存在的所有像素一一匹配来检测形成边界的像素。只对(从ROI方框的上边界到胎儿面部图像的一半)这样的部分执行形成边界的像素的检测程序是为了提高检测速度并且防止在不期望的区域产生边缘区域。
更具体地说，5×5的掩码的中心像素首先与区域“A”中存在的一个任意像素匹配。接着，第二确定单元35在区域“A”中选择分布在该任意像素附近并且与构成5×5的掩码的像素一一匹配的像素。此后，第二确定单元35将所选择的像素的像素值分别与构成5×5的掩码的像素的像素值、即“-1”、“1”或“0”相乘。然后其计算产生的值的总和，以便第二确定单元35最后将总值确定为该任意像素的新像素值，该任意像素存在于区域“A”中并且与5×5的掩码的中心像素匹配。第二确定单元35将上述程序应用于区域“A”中存在的所以像素并且确定形成目标对象的边界的像素，该目标对象具有像素值超过最新确定的像素值中预定的像素值的像素。根据本发明的优选实施例，对应于最新确定的像素值中高等级的20％的像素被确定为表示边缘像素的参考像素。
在检测出形成边界的像素之后，第二确定单元35计算边界像素的数量与构成区域“A”的所有像素的数量的比率以及形成边界的像素的方差。此后，该比率以及方差被应用于下列等式1，以便第二确定单元35确定目标对象的轮廓检测的可能性。
Di＝αR1i+βR2iEq.1其中，i是被输入到第二确定单元35中的相同类型目标对象的2D超声图像的数量，Di是表示轮廓检测可能性的确定数值，R1i是边界像素的数量与构成区域“A”的所有像素的数量的比率，以及R2i是形成边界的像素的方差。α和β是等式、例如等式1的系统，其通过相同类型目标对象的各种2D超声图像的轮廓检测试验获得。
特别地，第二确定单元从构造在3D超声诊断系统中的存储器中选择相同类型目标对象的各种2D超声图像。第二确定单元35为其轮廓能够被检测的2D超声图像赋予一个确定值“1”。其还为所选择的2D超声图像中其轮廓不能被检测的2D超声图像赋予一个确定值“0”。当在下列等式2中定义的均方误差值通过使用从存储器中选择的每个2D超声图像的确定值以及比率R1i和R2i而被最小化时，第二确定单元35计算α和β。
ϵ=Σi(Di-(αR1i+βR2i))2---Eq.2]]>为了使均方误差值最小化，第二确定单元35对α和β的等式2执行偏微分。在执行偏微分的等式的产生值变成“0”的情况下，α和β被确定为等式1的系数。第二确定单元35将R1i、R2I、α和β应用于等式1，以便通过第一确定单元33输入的ROI方框中的目标对象图像的确定值能够被计算。最后，第二确定单元35确定当所计算的确定值大于预定的设置值时能够检测ROI方框中的目标对象图像的轮廓。然后，所输入的2D超声图像被输出到(图2中所示的)轮廓检测单元37。另一方面，当所计算的确定值小于预定的设置值时，第二确定单元35确定不能检测ROI方框中的目标对象图像的轮廓并且停止3D超声诊断系统的操作。此后，第二确定单元35通知用户不能检测目标对象图像的轮廓。
为了简便起见，描述了通过检测在ROI方框的上边界和胎儿面部图像之间存在的边界像素来确定胎儿面部图像上表面的轮廓检测可能性的程序。然而，根据上述程序，应检测在ROI方框的底、左和右边界和胎儿面部图像之间存在的边界像素，以便确定胎儿面部图像的轮廓检测可能性。针对这样的程序，应使用转换图6A的5×5掩码的像素值的掩码。也即，针对ROI方框的底边界，应使用一个掩码，其中位于(图6A中所示的)5×5掩码的第一行42的像素值为“-1”并且位于第五行44的像素值为“1”。针对ROI方框的左边界，应使用一个掩码，其中位于(图6A中所示的)5×5掩码的第一列46的像素值为“-1”并且位于第五列48的像素值为“1”。同样地，针对ROI方框的右边界，应使用一个掩码，其中位于(图6A中所示的)5×5掩码的第一列46的像素值为“1”并且位于第五列48的像素值为“-1”。
由于根据上述程序执行用于确定胎儿面部图像的底、左和右表面的轮廓检测可能性的程序，所以在此将省略详细的说明。而且上述用于确定轮廓检测可能性的确定程序也可应用于任意的目标对象。
(图2中所示的)轮廓检测单元37检测从第二确定单元35输入的、ROI方框中存在的目标对象的轮廓。首先，为了补偿不具有足够亮度对比度的目标对象的缺陷，轮廓检测单元37经由通过使用Lee滤波器等去除存在于ROI方框内的目标对象中的斑点噪声来执行目标对象的平滑。此后，通过给构成ROI方框内的目标对象图像的像素提供在等式3中定义的像素值来执行目标对象的亮度对比度缩放。
F(x)=0forx≤low255×(x-low)(high-low)forlow<x<high255forhigh≤x---Eq.3]]>其中，F(x)是由轮廓检测单元37提供的新像素值，x是构成从第二确定单元35输入的目标对象图像的像素的老像素值，“low”是从第二确定单元35输入的目标对象的直方图低范围中的临界值，以及“high”是从第二确定单元35输入的目标对象的直方图高范围中的临界值。
图7A和7B表出说明在轮廓检测单元37中对ROI方框中的目标对象执行平滑和亮度对比度缩放的结果的实例，该ROI方框是以前设置的。
下面将描述用于检测胎儿图像的轮廓的程序。
轮廓检测单元37为ROI方框中的胎儿图像设置一个二值化阈值，并且通过作为参考的所设置的二值化阈值而执行胎儿图像的二值化。接着，轮廓检测单元37通过从胎儿图像的二值化区域中去除噪声区域来检测被确定为一部分胎儿图像的区域。也即，轮廓检测单元37检查二值化区域的亮度值并计算其平均值。在二值化区域中，如果二值化区域的亮度值小于平均亮度值，则该二值化区域被视为噪声区域并因此被去除。而且即使亮度值大于平均亮度值，由于噪声区域也可能存在，所以轮廓检测单元37为二值化区域的尺寸设置一个阈值。因此，如果二值化区域的尺寸小于该阈值，则该二值化区域被视为噪声区域并因此被去除。图8A和8B表示说明在轮廓检测单元37中对ROI方框中的胎儿图像执行二值化并从二值化区域中去除噪声区域的结果的实例。
下面考虑图9来说明在轮廓检测单元37中执行的轮廓检测程序。如果如图8B所示产生胎儿图像的二值化区域，轮廓检测单元37则将(没有示出的)ROI方框的上边界朝着(没有示出的)下边界移动，以便检测出胎儿图像的轮廓。接着，轮廓检测单元37提取对应于与ROI方框上边界接触的二值化区域的上表面的像素，并将所提取的像素与相邻像素相互连接。由于在胎儿图像的头部和身体之间存在一个深的凹部，所以获得如图9A所示的胎儿图像的上表面的轮廓。如果检测到如图9A所示的分离的轮廓，轮廓检测单元37则连接这些轮廓以由此获得最终的轮廓。图9B示出了一个实例，其将检测到的轮廓应用于在显示设备上显示的2D超声图像。
为了简便起见，只描述了胎儿图像上表面的轮廓检测程序。同样地，通过将上述程序应用于底、左和右边界可以检测底、左和后表面的轮廓。最后，轮廓检测单元37输出2D超声图像，在该2D超声图像中从胎儿图像中检测到的轮廓被显示给(如图2所示的)3D图像处理单元39。
3D图像处理单元39通过只再现目标对象的体积数据来形成3D超声图像。针对这种程序，3D图像处理单元39从存储在(没有示出的)电压数据获取单元中的电压数据中选择对应于由轮廓所包围的目标对象的体积数据，该轮廓由轮廓检测单元37输入。接着，3D图像处理单元39对所选择的体积数据执行扫描转换，然后应用典型的容积再现程序，以便3D超声图像能够被更精确地显示。
如上所述，由于ROI方框是自动产生的并且该ROI方框中存在的目标对象图像的轮廓是自动进行检测的，所以ROI方框产生和目标对象图像的轮廓检测的时间消耗能够被减少。而且由于自动产生的ROI方框的尺寸的适宜性以及从ROI方框内的目标对象图像中检测轮廓的可能性是确定的，所以由体积数据再现和目标对象图像的轮廓检测产生的误差能够被减少。同样地，更精确的目标对象的3D超声图像可提供给3D超声诊断系统的用户。
虽然已根据本发明的优选实施例描述和解释了本发明，但对本领域的技术人员而言显而易见的是在不偏离本发明的主要原理和技术启示的情况下可以进行改变和修改，本发明的主要原理和技术启示应只由所附的权利要求的范围所限定。
权利要求
1.一种用于形成三维(3D)超声图像的设备，包括第一装置，用于在2D超声图像上产生感兴趣区域(ROI)方框；第二装置，用于检测该RIO方框中的目标对象的轮廓；以及第三装置，用于通过再现存在于检测到的轮廓中的体积数据来形成3D超声图像。
2.如权利要求1所述的设备，还包括第一确定装置，用于确定ROI方框的尺寸是否适合目标对象的2D超声图像的尺寸。
3.如权利要求2所述的设备，还包括第二确定装置，用于确定目标对象的2D超声图像的轮廓是否可检测。
4.如权利要求1所述的设备，其中第一装置包括第一单元，用于将2D超声图像分割成目标对象图像和背景图像；第二单元，用于根据2D超声图像上的目标对象的类型来设置具有预定尺寸的ROI方框；以及第三单元，用于调节该ROI方框的尺寸以适合目标对象的尺寸。
5.如权利要求4所述的设备，其中第三单元执行以下功能在执行了ROI方框中目标对象图像的二值化之后，去除目标对象图像的二值化区域中的噪声区域；通过为二值化区域提供预定的像素值而产生目标对象的二值化图像以及二值化图像的表面曲线；以及通过将每个边界移动至最靠近该表面曲线的像素来调节ROI方框的边界。
6.如权利要求2所述的设备，其中第一确定装置通过计算针对目标对象的每种类型的2D目标对象图像的、预定的标准直方图和针对ROI方框中目标对象图像的直方图的均方误差值来确定ROI方框的尺寸是否合适。
7.如权利要求3所述的设备，其中第二确定装置执行以下功能通过使一种预定类型的掩码的像素与ROI方框中目标对象图像的像素一一匹配来提取形成目标对象边界的像素；通过根据以下等式应用边界像素的数量与和掩码一一匹配的所有像素的数量的比率以及形成边界的像素的方差来确定目标对象的轮廓检测的可能性Di＝αR1i+βR2i其中，i是先前存储在存储器中的相同类型目标对象的2D超声图像的数量，Di是表示轮廓检测可能性的确定数值，R1i是边界像素的数量与和掩码一一匹配的所有像素的数量的比率，以及R2i是形成边界的像素的方差，以及其中通过相同类型目标对象的各种2D超声图像的轮廓检测试验获得的系数α和β是当在以下等式中定义的均方误差值最小化时的值ϵ=Σi(Di-(αR1i+βR2i))2]]>其中为了使均方误差值最小化，针对每个α和β执行均方误差值的偏微分，以及其中当产生的偏微分值为“0”时α和β被确定为系数。
8.如权利要求1所述的设备，其中第二装置执行以下功能执行目标对象图像的亮度对比度缩放，以便使RIO方框中的目标对象有足够的亮度对比度；在执行了RIO方框中目标对象的二值化之后，从目标对象的二值化区域中去除噪声区域；提取在二值化区域的表面内存在的像素，该二值化区域的表面通过调节RIO方框边界而与RIO方框边界接触；以及通过连接所提取的像素来检测目标对象的轮廓。
9.一种用于形成三维(3D)超声图像的方法，包括以下步骤a)在2D超声图像上产生感兴趣区域(ROI)方框；b)检测该RIO方框中的目标对象的轮廓；以及c)通过再现存在于检测到的轮廓中的体积数据来形成3D超声图像。
10.如权利要求9所述的方法，还包括以下步骤，即确定ROI方框的尺寸是否适合目标对象的2D超声图像的尺寸。
11.如权利要求10所述的方法，还包括以下步骤，即确定目标对象的2D超声图像的轮廓是否可检测。
12.如权利要求9所述的方法，其中步骤a)包括以下步骤a1)将2D超声图像分割成目标对象图像以及与目标对象相邻的背景图像；a2)在2D超声图像上设置具有预定尺寸的ROI方框；以及a3)调节该ROI方框的尺寸以适合目标对象的2D超声图像的尺寸。
13.如权利要求12所述的方法，其中步骤a3)包括以下步骤a3-1)在执行了ROI方框中目标对象图像的二值化之后，去除目标对象图像的二值化区域中的噪声区域；a3-2)通过为二值化区域提供预定的像素值而产生目标对象的二值化图像以及二值化图像的表面曲线；以及a3-3)通过将每个边界移动至最靠近该表面曲线的像素来调节ROI方框的边界。
14.如权利要求10所述的方法，其中步骤d)包括以下步骤，即通过计算针对目标对象的每种类型的2D目标对象图像的、预定的标准直方图和针对ROI方框中目标对象图像的直方图的均方误差值来确定ROI方框的尺寸。
15.如权利要求11所述的方法，其中步骤e)包括以下步骤e1)通过使一种预定类型的掩码的像素与ROI方框中目标对象图像的像素一一匹配来提取形成目标对象边界的像素；e1)通过根据以下等式应用边界像素的数量与和掩码一一匹配的所有像素的数量的比率以及形成边界的像素的方差来确定目标对象的轮廓检测的可能性Di＝αR1i+βR2i其中，i是先前存储在存储器中的相同类型目标对象的2D超声图像的数量，Di是表示轮廓检测可能性的确定数值，R1i是边界像素的数量与和掩码一一匹配的所有像素的数量的比率，以及R2i是形成边界的像素的方差，以及其中通过相同类型目标对象的各种2D超声图像的轮廓检测试验获得的系数α和β是当在以下等式中定义的均方误差值最小化时的值ϵ=Σi(Di-(αR1i+βR2i))2]]>其中为了使均方误差值最小化，针对每个α和β执行均方误差值的偏微分，以及其中当所产生的偏微分值为“0”时α和β被确定为系数。
16.如权利要求9所述的方法，其中步骤b)包括以下步骤b1)执行目标对象图像的亮度对比度缩放，以便使RIO方框中的目标对象有足够的亮度对比度；b2)在执行了RIO方框中目标对象的二值化之后，从目标对象的二值化区域中去除噪声区域；b3)提取在二值化区域的表面内存在的像素，该二值化区域的表面通过调节RIO方框边界而与RIO方框边界接触；以及b4)通过连接所提取的像素来检测目标对象的轮廓。
全文摘要
本发明涉及通过自动检测目标对象的轮廓而只利用该轮廓内存在的体积数据来形成3D超声图像的3D超声诊断，包括第一单元，用于在2D超声图像上产生感兴趣区域(ROI)方框；第二单元，用于检测该RIO方框中的目标对象的轮廓；以及第三单元，用于通过再现存在于检测到的轮廓中的体积数据来形成3D超声图像。
文档编号G06T15/08GK1676104SQ200510006818
公开日2005年10月5日 申请日期2005年1月28日 优先权日2004年4月1日
发明者沈载允, 郭炳周, 金南哲, 金祥铉 申请人:株式会社美蒂森