著录项信息
专利名称 | 用于超声成像的设备和方法 |
申请号 | CN200780011611.4 | 申请日期 | 2007-04-02 |
法律状态 | 暂无 | 申报国家 | 中国 |
公开/公告日 | 2009-06-03 | 公开/公告号 | CN101448460 |
优先权 | 暂无 | 优先权号 | 暂无 |
主分类号 | A61B8/00 | IPC分类号 | A;6;1;B;8;/;0;0查看分类表>
|
申请人 | 阿洛卡株式会社 | 申请人地址 | 日本东京
变更
专利地址、主体等相关变化,请及时变更,防止失效 |
权利人 | 株式会社日立制作所 | 当前权利人 | 株式会社日立制作所 |
发明人 | 田村正 |
代理机构 | 北京泛华伟业知识产权代理有限公司 | 代理人 | 王勇 |
摘要
使用实际的超声组织衰减来计算增益补偿轮廓,该增益补偿轮廓用于产生均匀图像。轴向、侧向、高度增益轮廓用于修正在每个方向上的衰减和超声变化。此外,描述了自动增益补偿的自动激活。
1.一种方法,用于从超声图像数据推导增益补偿数据,包括:接收对应图像扫描的超声图像数据;
从超声图像数据组合图像帧,所述图像帧具有图像数据值的列和行,其中,列表示侧向位置而行表示轴向位置;
将图像帧分成包括有效图像数据值的多个部分;
对于所述多个部分的每一个,基于该部分的有效图像数据值计算代表信号值;
通过如下步骤计算每对相邻部分的代表信号值之间的增益变化:
选择所述多个部分中的一个作为参考部分;
选择分布方向;
对于在分布方向上进行的增益变化计算,从比某一部分距离参考部分更远的相邻部分的代表信号值减去所述某一部分的代表信号值;以及对于不在分布方向上进行的增益变化计算,从与某一部分距离参考部分更近的相邻部分的代表信号值减去所述某一部分的代表信号值;
从基于分布方向计算的增益变化计算增益变化分布序列值;
从增益变化分布序列值计算增益变化曲线;以及
从增益变化曲线计算增益补偿曲线,其中,增益补偿曲线的值与位置对应。
2.权利要求1所述的方法,还包括:
检查所述帧中的每个图像数据值,找到在预定的下限和上限之外的值;以及将在所述预定的上限和下限之外的所述图像数据值标记为无用。
3.权利要求2所述的方法,其中,所述超声图像数据包括RF数据、IF数据、基带信号数据或线数据。
4.权利要求3所述的方法,其中,每个部分的代表信号值使用柱状图、百分值、平均值、中间值或众数来计算。
5.权利要求4所述的方法,其中,每个部分包括相同数目的图像数据值、不同数目的图像数据值或与相邻部分共享图像数据值。
6.权利要求1所述的方法,其中,计算增益变化分布序列值包括:
使用参考部分的代表信号值作为一个增益变化分布序列值;以及
如果在分布方向上,通过将增益变化加到之前的增益变化分布序列值,或者如果不在分布方向上,通过从之前的增益变化分布序列值减去增益变化,来产生新的增益变化分布序列值。
7.权利要求6所述的方法,其中,增益变化曲线是通过曲线拟合或顺序的增益变化分布序列值的插值计算的。
8.权利要求7所述的方法,其中,增益补偿曲线是通过改变增益变化曲线的极性计算的。
9.权利要求1所述的方法,还包括将每个增益变化值的绝对值同预定的限制相比较,其中,如果绝对增益变化值大于预定限制,用0代替该增益变化值。
10.权利要求1所述的方法,还包括将每个增益变化值的绝对值同预定的限制相比较,其中,如果绝对增益变化值大于预定限制,用预定值代替该增益变化值。
11.权利要求1所述的方法,其中,多个部分是轴向部分。
12.权利要求1所述的方法,其中,多个部分是侧向部分。
13.一种方法,用于为超声图像体推导增益补偿数据,包括:
接收对应多个连续图像扫描的超声图像数据;
从超声图像数据组合图像体,所述图像体具有三维图像数据值,其中,列表示侧向位置、行表示轴向位置而帧表示高度位置;
将图像体分成包括有效图像数据值的多个体部分;
对于所述多个体部分的每一个,基于该体部分的有效图像数据值计算代表信号值;
通过如下步骤计算每对相邻体部分的代表信号值之间的增益变化:
选择所述多个体部分中的一个作为参考体部分;
选择分布方向;
对于在分布方向上进行的增益变化计算,从比某一体部分距离参考体部分更远的相邻体部分的代表信号值减去所述某一体部分的代表信号值;以及
对于不在分布方向上进行的增益变化计算,从比某一体部分距离参考体部分更近的相邻体部分的代表信号值减去所述某一体部分的代表信号值;
从基于分布方向计算的增益变化计算增益变化分布序列值;
从增益变化分布序列值计算增益变化曲线;以及
从增益变化曲线计算增益补偿曲线,其中,增益补偿曲线的值与位置对应。
14.权利要求13所述的方法,还包括:
检查所述图像体中的每个图像数据值,找到在预定的下限和上限之外的值;以及将在所述预定的上限和下限之外的所述图像数据值标记为无用。
15.权利要求14所述的方法,其中,超声图像数据包括RF数据、IF数据、基带信号数据或线数据。
16.权利要求15所述的方法,其中,每个体部分的代表信号值使用柱状图、百分值、平均值、中间值或众数来计算。
17.权利要求16所述的方法,其中,每个体部分包括相同数目的图像数据值、不同数目的图像数据值或与相邻体部分共享图像数据值。
18.权利要求13所述的方法,其中,计算增益变化分布序列值包括:使用参考体部分的代表信号值作为一个增益变化分布序列值;以及
如果在分布方向上,通过将增益变化加到之前的增益变化分布序列值,或者如果不在分布方向上,通过从之前的增益变化分布序列值减去增益变化,来产生新的增益变化分布序列值。
19.权利要求18所述的方法,其中,增益变化曲线是通过顺序的增益变化分布序列值的曲线拟合计算的。
20.权利要求19所述的方法,其中,增益补偿曲线是通过改变增益变化曲线的极性计算的。
21.权利要求13所述的方法,还包括将每个增益变化值的绝对值同预定的限制相比较,其中,如果绝对增益变化值大于预定限制,用0代替增益变化值。
22.权利要求13所述的方法,还包括将每个增益变化值的绝对值同预定的限制相比较,其中,如果绝对增益变化值大于预定限制,用预定值代替该增益变化值。
23.权利要求13所述的方法,其中,多个体部分是高度体部分。
24.权利要求13所述的方法,其中,多个体部分是轴向体部分。
25.权利要求13所述的方法,其中,多个体部分是侧向体部分。
26.用于推导用于超声系统的增益补偿数据的自动增益补偿系统,其中,所述超声系统具有配置用于在获取超声图像数据时调整所接收的信号增益的放大器,所述自动增益补偿系统包括:
自动增益处理器,配置用于:
接收对应图像扫描的超声图像数据;
从超声图像数据组合图像帧,图像帧具有图像数据值的列和行,其中,列表示侧向位置而行表示轴向位置;
将图像帧分成包括有效图像数据值的多个部分;
对于多个部分的每一个,基于该部分的有效图像数据值计算代表信号值;
通过如下步骤计算每对相邻部分的代表信号值之间的增益变化:
选择多个部分中的一个作为参考部分;
选择分布方向;
对于在分布方向上进行的增益变化计算,从比某一部分距离参考部分更远的相邻部分的代表信号值减去所述某一部分的代表信号值;以及
对于不在分布方向上进行的增益变化计算,从比某一部分距离参考部分更近的相邻部分的代表信号值减去所述某一部分的代表信号值;
从基于分布方向计算的增益变化计算增益变化分布序列值;
从增益变化分布序列值计算增益变化曲线;和
从增益变化曲线计算增益补偿曲线,其中,增益补偿曲线的值与位置对应。
27.用于推导用于超声系统的增益补偿数据的自动增益补偿系统,其中,所述超声系统具有配置用于在获取超声图像数据时调整所接收的信号增益的放大器,所述自动增益补偿系统包括:
自动增益处理器,配置用于:
接收对应多个连续的图像扫描的超声图像数据;
从超声图像数据组合图像体,其中,所述图像体具有三维图像数据值,其中,列表示侧向位置、行表示轴向位置而帧表示高度位置;
将所述图像体分成包括有效图像数据值的多个体部分;
对于多个体部分的每一个,基于该体部分的有效图像数据值计算代表信号值;
通过如下步骤计算相邻的体部分代表信号值之间的增益变化:
选择多个体部分中的一个作为参考体部分;
选择分布方向;
对于在分布方向上进行的增益变化计算,从比某一体部分距离参考体部分更远的相邻体部分的代表信号值减去所述某一体部分的代表信号值;以及
对于不在分布方向上进行的增益变化计算,从比某一体部分距离参考体部分更近的相邻体部分的代表信号值减去所述某一体部分的代表信号值;
从基于分布方向计算的增益变化计算增益变化分布序列值;
从增益变化分布序列值计算增益变化曲线;和
从增益变化曲线计算增益补偿曲线,其中,增益补偿曲线值与位置对应。
用于超声成像的设备和方法\n[0001] 相关申请的交叉引用\n[0002] 本申请要求提交于2006年3月31日的美国临时申请No.60/744,019和提交于\n2006年3月31日的美国临时申请60/744,020的优先权,这两个申请均通过全文引用的方\n式包括于此。\n技术领域\n[0003] 本发明涉及超声成像。更具体地,本发明涉及当执行超声成像时对所接收到的超\n声信号增益的自动调整。\n背景技术\n[0004] 超声被用于对人体进行成像,以诊断诸如胎儿、心脏、肝脏、肾和其它器官的各种\n医学情况。超声由超声换能器经皮肤传输到在人体内的组织。超声被超声散射体散射并被\n超声换能器接收。所接收的超声被超声换能器转换成电信号并被处理以产生组织的图像。\n[0005] 超声经人类组织以大约0.5dB/MHz/cm的速度衰减。超声波束的强度在穿透组织\n时减弱。因此,在不同深度的两个相同目标将产生不同的回声,由较近目标产生的回声比另\n一个大。这个问题将通过使用时间-增益补偿(TGC)避开,在TGC中所接收信号放大器的\n增益以深度(时间)的函数形式增加以补偿由于衰减造成的能量的损失。已经采用多种形\n式的时间-增益补偿。在现代的扫描器中,TGC形状可以方便地调整以最优化图像或应用。\n[0006] 传统的超声成像系统通常装备有TGC用于对这种类型的衰减进行补偿。然而,大\n部分系统使用固定的TGC或操作者可调的TGC。如果不对所接收的超声信号进行补偿,所得\n的超声图像将在小深度处较亮,在较大深度处较暗,导致不均匀的图像。\n[0007] 固定的TGC对所有的成像使用预定的衰减曲线,这不是最优化的,因为在每个患\n者体内的衰减程度是不同的。可调的TGC是由操作者通过滑动在控制台上的分压计来控制\n的。典型的超声系统包括数个分压计,以便操作者调整来产生提供均匀图像的衰减补偿曲\n线。操作者需要为每个患者和人类身体的每个部位设置TGC分压计来成像。如果TGC的设\n置设定得不准确,诊断质量将受影响。\n[0008] 操作者可调的TGC耗时,并且可能不利地影响对患者的诊断。这就存在了使用来\n自患者的信号对每个超声图像自动调整TGC的需要。\n发明内容\n[0009] 发明人揭示出期望有一种系统和方法,其使用实际的超声组织衰减来计算增益补\n偿轮廓,然后该增益补偿轮廓用于补偿以修正并产生用于显示的均匀图像。轴向、侧向、高\n度(用于体成像)及其结合的增益轮廓用于修正增益并减小跨过图像的变化。\n[0010] 除了提供增益补偿外,还教示了增益补偿的自动激活。在一个实施例中,运动传感\n器用于探测探头运动并激活增益补偿系统。在另一个实施例中,运动由超声图像探测。\n[0011] 本发明的一个方面提供一种用于从超声图像数据推导增益补偿数据的方法。根据\n本发明的该方面的方法包括接收对应图像扫描的超声图像数据;从超声图像数据组合图像\n帧,所述图像帧具有图像数据值的列和行,其中,列表示侧向位置而行表示轴向位置;将图\n像数据帧分隔成多个部分;计算每个部分的代表信号值;计算相邻部分代表信号值之间的\n增益变化;从所述各部分之间的增益变化计算增益变化分布;从增益变化分布计算增益变\n化曲线;以及从增益变化曲线计算增益补偿曲线,其中,补偿增益值对应位置。\n[0012] 在本发明的另一方面中,超声成像数据包括RF数据、IF数据、基带信号数据、所探\n测到的基带信号数据、对数压缩的探测到的基带信号数据、线数据或图像数据。\n[0013] 本发明的另一方面提供一种用于推导用于超声图像体的增益补偿数据的方法。根\n据本发明该方面的方法包括接收对应多个连续图像扫描的超声图像数据;从多个图像扫描\n组合图像体,该图像体具有三维图像数据值,其中,列表示侧向位置、行表示轴向位置而帧\n表示高度位置;将该体分成多个体部分,计算每个体部分的代表信号值;计算相邻体部分\n代表信号值之间的增益变化;从体部分之间的增益变化计算增益变化分布;从增益变化分\n布计算增益变化曲线;以及从增益变化曲线计算增益补偿曲线,其中,增益补偿值对应位\n置。\n[0014] 本发明的另一方面提供一种方法,用于激活用于超声系统的自动增益补偿数据系\n统。根据该方面的方法包括接收对应连续图像扫描的超声图像数据,并且执行对每个所接\n收的扫描的图像数据之间的相关性分析,其中,该相关性分析返回代表探头运动的值,并且\n如果所探测到的改变大于预定值,可以认为探头被移动,从而启动新的增益补偿计算。\n[0015] 另一方面,相关性分析是相关性、绝对差异D的和(SAD)、差异平方的和(SSD)、差\nq\n异立方的和以及差异的幂D 的和。\n[0016] 本发明的另一方面提供一种用于激活用于超声系统的自动增益补偿数据系统的\n方法。根据该方面的方法包括接收对应连续图像扫描的超声图像数据,并观测轴向部分、侧\n向部分或高度部分之间的增益变化的改变,其中,如果所探测到的改变大于预定值,可以认\n为探头被移动,启动新的增益补偿计算。\n[0017] 本发明的另一方面提供一种用于激活用于超声系统的自动增益补偿数据系统的\n方法。根据该方面的方法包括接收对应连续图像扫描的超声图像数据,并通过相关性、SAD\n或SSD、SCD或差异的幂的和探测运动,其中,如果所探测到的改变大于预定值,可以认为探\n头被移动,启动新的增益补偿计算,以对新的增益改变进行补偿。\n[0018] 本发明的另一方面提供一种用于激活用于超声系统的自动增益补偿数据系统的\n方法。根据该方面的方法包括使用附着在超声探头上的运动传感器探测运动,其中,如果运\n动大于预定值,可以认为探头被移动,启动新的增益补偿计算。\n[0019] 本发明的另一方面提供在获取超声图像数据的同时,用于控制TGC放大器和超声\n系统的轴向、侧向以及高度可变增益放大器的增益的方法。根据本方面的方法包括从超声\n图像接收轴向、侧向和高度增益补偿数据;从超声系统接收超声信号;从所接收的超声信\n号计算总的增益值;从之前的轴向、侧向和高度结合增益数据和所接收的轴向、侧向和高度\n增益补偿数据以及总的增益数据计算TGC、轴向、侧向和高度放大器增益数据;并分配TGC、\n轴向、侧向和高度放大器增益数据至各个放大器。\n[0020] 本发明的另一方面是用于为超声系统推导增益补偿数据的自动增益补偿系统,其\n中,该超声系统具有配置为在获取超声图像数据的同时调整所接收的信号增益的放大器。\n根据该方面的系统包括自动增益处理器,该自动增益处理器配置为用于接收对应图像扫描\n的超声图像数据;从超声图像数据组合图像数据帧,图像帧具有图像数据值的列和行,其\n中,列表示侧向位置而行表示轴向位置;将图像数据帧分隔成多个部分;计算每个部分的\n代表信号值;计算相邻部分代表信号值之间的增益变化;从部分之间的增益变化计算增益\n变化分布;从增益变化分布计算增益变化曲线;和从增益变化曲线计算增益补偿曲线,其\n中,补偿增益值与位置对应。\n[0021] 本发明的另一方面是用于推导用于超声系统的增益补偿数据的自动增益补偿系\n统,其中,该超声系统具有配置为在获取超声图像数据的同时调整所接收的信号增益的放\n大器。根据该方面的系统包括自动增益处理器,该自动增益处理器配置为用于接收对应多\n个连续的图像扫描的超声图像数据;从多个图像扫描组合图像体,图像体具有三维图像数\n据值,其中,列表示侧向位置、行表示轴向位置而帧表示高度位置;将体分隔成多个体部分;\n计算每个体部分的代表信号值;计算相邻体部分代表信号值之间的增益变化;从体部分之\n间的增益变化计算增益变化分布;从增益变化分布计算增益变化曲线;和从增益变化曲线\n计算增益补偿曲线,其中,补偿增益值与位置对应。\n[0022] 本发明的另一方面是一种用于激活用于超声系统的自动增益补偿数据系统的系\n统。根据该方面的系统包括用于接收对应连续图像扫描的超声图像数据的装置,和用于执\n行对每个所接收的扫描的图像数据之间的相关性分析的装置,其中,相关性分析返回表示\n探头运动的值,并且如果所探测到的改变大于预定值,可以认为探头被移动,启动新的增益\n补偿计算。\n[0023] 本系统的另一方面,相关性分析是相关性、绝对差异D的和(SAD)、差异平方的和\nq\n(SSD)、差异立方的和以及差异的幂D 的和。\n[0024] 本发明的另一方面是一种用于激活用于超声系统的自动增益补偿数据系统的系\n统。根据该方面的系统包括接收对应连续图像扫描的超声图像数据的装置,以及用于观测\n轴向、侧向或高度部分之间的增益变化的改变的装置,其中,如果所探测到的改变大于预定\n值,可以认为探头被移动,启动新的增益补偿计算。\n[0025] 本发明的另一方面是一种用于激活用于超声系统的自动增益补偿数据系统的系\n统。根据该方面的系统包括接收对应连续图像扫描的超声图像数据的装置,以及通过相关\n性、SAD或SSD、SCD或差异的幂的和来探测运动的装置,其中,如果所探测到的改变大于预\n定值,可以认为探头被移动,启动新的增益补偿计算,以对新的增益改变进行补偿。\n[0026] 本发明的另一方面是一种用于激活用于超声系统的自动增益补偿数据系统的系\n统。根据该方面的系统包括使用附着在超声探头上的运动传感器探测运动的装置,其中,如\n果运动大于预定值,可以认为探头被移动,启动新的增益补偿计算。\n[0027] 该系统的另一方面,其中,运动传感器是速度计、加速计、陀螺仪、运动跟踪设备和位置传感器。\n[0028] 本发明的另一方面是在获取超声图像数据的同时,用于控制TGC放大器和超声系\n统的轴向、侧向以及高度可变增益放大器的增益的系统。根据本方面的系统包括处理器,该\n处理器配置为用于从超声图像接收轴向、侧向和高度增益补偿数据;从超声系统接收超声\n信号;从所接收的超声信号计算总的增益值;从之前的轴向、侧向和高度结合增益数据和\n所接收的轴向、侧向和高度增益补偿数据以及总的增益数据计算TGC、轴向、侧向和高度放\n大器增益数据;并分配TGC、轴向、侧向和高度放大器增益数据至各个放大器。\n[0029] 下面在附图和说明书中描述本发明的一个或多个实施例的细节。本发明的其它特\n征、目的和优势将通过说明书、附图以及权利要求阐明。\n附图说明\n[0030] 图1是具有自动增益补偿系统的典型的超声系统的示例性实施例。\n[0031] 图2是示例性的图像数据帧。\n[0032] 图3是示例性的方法。\n[0033] 图4是在图2中所示的图像数据帧的示例性的轴向分割。\n[0034] 图5是示例性的轴向增益补偿曲线。\n[0035] 图6是在轴向增益补偿之前的具有10个轴向部分柱状图的超声图像。\n[0036] 图7是在轴向增益补偿之后的具有10个轴向部分柱状图的图6中所示的超声图\n像,示出了更加均匀的垂直强度分布。\n[0037] 图8是图2中所示的图像数据帧的示例性的侧向分割。\n[0038] 图9是示例性的侧向增益补偿曲线。\n[0039] 图10是在侧向增益补偿之前的具有四个侧向部分柱状图的超声图像。\n[0040] 图11是在侧向增益补偿之后的具有四个侧向柱状图的在图10中所示的超声图\n像,示出了更加均匀的水平强度分布。\n[0041] 图12是五个图像数据帧体的示例性高度分割。\n[0042] 图13是示例性的高度增益补偿曲线。\n[0043] 图14是5个图像数据帧体的示例性轴向分割。\n[0044] 图15是示例性的轴向体增益补偿曲线。\n[0045] 图16是5个图像数据帧体的示例性的侧向分割。\n[0046] 图17是示例性侧向体增益补偿曲线。\n具体实施方式\n[0047] 下面将参照附图描述本发明的实施例,其中,类似的标号表示类似的元件。在对本\n发明的实施例详细解释之前,应该理解,本发明并未将其应用限制于在下面的描述或附图\n中所示的例子中的细节。本发明可以具有其它实施例并且能够在多种应用和以多种方式实\n践或执行。而且,应该理解,这里使用的措辞和术语是用于描述的目的并且不应该理解为限\n制。使用“包括”、“包含”或“具有”以及这里所使用的这些词的变体是指包括其后所列项目和其等价物,以及附加项目。术语“安装”、“连接”和“耦合”使用其广义并且包括直接和间接的安装、连接和耦合。此外,“连接的”和“耦合的”并不限定是物理或机械连接或耦合。\n[0048] 值得注意的是本发明并不限于所描述或图中所暗示的任何特定的软件语言。本领\n域的普通技术人员将会理解,多种可选的软件语言可以用于实现本发明。还应该理解,一些\n部件和项目是将其当成硬件元件来进行说明和描述的,这是本领域中的常用手法。然而,本\n领域的普通技术人员,在阅读该详细说明的基础上,将能够理解在至少一个实施例中,方法\n和系统中的部件可以以软件或硬件的方式实现。\n[0049] 本发明的实施例提供方法、系统和存储计算机可读指令的计算机可用的介质,该\n机算机可读指令有效地处理超声图像数据以得到轴向、侧向和高度增益补偿曲线,用于实\n时诊断成像应用。本发明从一个超声图像数据帧有效地分析数据,或者对于图像数据体,分\n析多个图像数据帧,并且推导用于每个图像维度的增益补偿轮廓或曲线。本发明是模块化\n结构并且可被配置为驻存在具有机载电源的封装中的硬件,或配置为诸如实在地包含在程\n序存储设备上的用于由处理器执行的应用程序的软件。用于执行的应用代码可以驻存在多\n个不同类型的计算机可读介质上。\n[0050] 作为背景技术,超声波扫描术(超声波检查法)使用探头,该探头包括一个或多个\n声学换能器用于发射声音脉冲到材料中。通常声波由装在探头中的压电换能器阵列通过产\n生短的、强声音脉冲而形成。用于医学成像的频率通常在从1至20MHz的范围内,这属于中\n到高超声频率,并且可以从由换能器发射的所有单独脉冲的和产生单个的、聚焦的声波。较\n高频率具有相应较低的波长并产生较高的空间分辨率图像。超声波检查法可以使用采用基\n频的系统,也可以使用采用基频的谐波的系统。\n[0051] 当声波遇到具有不同声阻抗的材料时,部分声波反射,换能器探测这部分声波作\n为回声。返回的声波振动换能器的元件并将该振动转变成电脉冲,该电脉冲从探头被发送\n至处理器,并在处理器中被处理并转换成图像。回声传输回探头所使用的时间被测量并用\n于计算产生回声的组织分界面的深度。\n[0052] 为了产生二维图像,使用阵列换能器(例如,相控阵列,线性阵列、凸阵列和其它\n阵列)电子扫描超声束。可选地,单个元件换能器或角阵列换能器也可使用并且机械地扫\n描以产生二维图像。所接收的RF数据还进一步被处理并用于构建超声图像。阵列换能器\n可以装备有机械扫描设备,机械扫描设备允许机械地扫描换能器以辅助阵列的电子扫描,\n从而可以得到用于三维成像的体数据。\n[0053] 假定固定的声音速度为大约1540m/s,超声系统可以基于所测得的超声传输时间\n确定器官或目标组织的位置。大回声信号在图像中将表现得亮,而小回声信号将表现得暗。\n[0054] 所接收到的数据指的是RF数据值并且其表现形式类似于矩阵。例如,用i表示轴\n向行并用j表示侧向列,其中,i=1,2,3,...,M并且j=1,2,3,...,N。RF数据值ai,j通常是有两极的(±)多比特值。例如,2048×128RF数据帧可以具有262,144个图像数据值\nai,j。\n[0055] 图1中所示的是典型的超声系统101和自动增益补偿系统103的示例性实施例。\n自动增益补偿系统103可以用于任何超声成像系统101。\n[0056] 超声系统101包括具有换能器阵列的探头105,换能器阵列包括多个元件\n107x-107X,其中,X是可以布置在线中作为线性阵列的矩形元件的数目,是布置在行和列中\n作为二维阵列的方形元件的数目,或是同心地布置为环形阵列的环形元件的数目。探头105\n耦合到通常具有X个同轴电缆的微同轴电缆109。\n[0057] 电缆109可以耦合至高压多路复用器111以减少硬件复杂性,或可以直接耦合至\n发射/接收开关113,其阻止高传输电压脉冲115。发射/接收开关113的输出耦合至多个\n通道Y。每个通道可以包括耦合至时间增益补偿(TGC)放大器119y的低噪声前置放大器\n117y。每个TGC放大器119y包括可变增益放大器部分和可编程增益放大器部分,其可以产\n生用于轴向方向的固定增益。\n[0058] 所有通道Y(119y-119Y)的输出耦合至数字束形成器121,其转向并聚焦所述多个\n单独通道y的输出,使其成为一个束。该束形成器121可以包括多个A/D转换器123y-123Y,\n其每个耦合至可变延迟125y-125Y,例如FIFO和/或相转子。该延迟的通道耦合至数字加\n法器127并输出。\n[0059] 束形成器121的输出耦合至级联放大器部分129,该级联放大器部分129可包括探\n测和对数放大器131用于在下游的增益和图像处理之前压缩动态范围。该探测和对数放大\n器131被耦合至可变增益轴向放大器133,可变增益轴向放大器133耦合至可变增益侧向放\n大器135。该放大器部分129耦合至图像处理器137用于B模式成像处理,耦合至彩色多普\n勒处理器139用于流动模式处理,耦合至光谱多普勒处理器138,和耦合至体处理器142,该\n体处理器142具有可变增益高度放大器143和三维扫描转换器145。处理器137、139和转\n换器145的输出被耦合至显示器141。处理器137、139和转换器145存储图像帧值和组合\n的图像帧。\n[0060] 如上所述,阵列探头105可以包含从32到300个换能器元件107x,可以通过适当\n地延迟去往用于传输的元件的信号并且在接收后延迟(125y)信号来聚焦并转向这些换能\n器元件107x。聚焦的目的是提高侧向图像分辨率。激发探头105的中心元件的脉冲115相\n对于激发探头105的周围的元件的脉冲115延迟时间周期,使得所有换能器脉冲107x同时\n到达在组织中的焦点P。在接收过程中,换能器元件107x可以从浅区域到深区域连续地被\n聚焦。\n[0061] 在接收时,来自多个探头105的换能器元件107x的超声信号被束形成器121聚焦\n并转向为一束。重复这个过程,形成与之前的束相邻的另一束,直到到达末端位置。从32\n至300个换能器元件107x,可以通过聚焦并转向形成50至500个侧向束,每束在轴向方向\n包括从大约100至1000个图像数据值ai,j。在扫描转换以后,多个束和图像数据值形成一\n个二维图像数据帧。\n[0062] 示例性的超声系统101使用数字束形成器121。然而,在其他实施例中,可以使用\n采用模拟延迟部件的模拟束形成器。此外,也可以使用采用声学透镜来聚焦的使用机械扫\n描、单个元件换能器107的超声系统换能器。\n[0063] 自动增益补偿系统103接收从束形成器121输出的图像数据147、从可变增益侧向\n放大器129输出的图像数据148和从可变增益高度放大器143输出的图像数据149作为输\n入。自动增益补偿系统103输出控制信号到可变增益TGC放大器119y、可变增益轴向放大\n器133、可变增益侧向放大器135和可变增益高度放大器143。\n[0064] 自动增益补偿系统103可以包括耦合至帧组合器153的用于压缩并处理从束形成\n器121输出的RF信号的探测器/对数放大器151。从可变增益侧向放大器135和可变增益\n高度放大器143输出的图像数据148、149被耦合至帧组合器153。图像数据可以包括幅值\n2 3 4 q\n(p)、强度(p)、(p)、(p)、p,或由幅值或强度推导出的任何信号,其中,q是实数。在一个\n实施例中,耦合至帧组合器153的是配置用于并行处理的轴向增益变化引擎155、侧向增益\n变化引擎157、以及高度增益变化引擎159。可选地,该轴向、侧向和高度引擎可以功能上结\n合到一个引擎中并串行地处理。这三个引擎155、157、159输出补偿数据耦合至控制器161,\n控制器161可分配补偿数据输出至可变增益TGC放大器119y、可变增益轴向放大器133、可\n变增益侧向放大器135和可变增益高度放大器143,用于为在一个或多个方向上的信号增\n益变化提供补偿。\n[0065] 为了执行增益补偿数据计算,自动增益补偿系统103可以由用户、运动探测启动。\n对于运动探测,处理器165耦合至探测器151和帧组合器153。该运动处理器165可以使用\n从束形成器121输出的RF信号来执行相关,或使用基带、线数据或来自帧组合器153的图\n像数据来感测何时超声探头105从一个位置移动至另一个位置。可以在帧之间或运动之前\n和之后执行相关。\n[0066] \n[0067] \n[0068] 可选地,处理器165可以使用绝对差异的和(SAD)。\n[0069] \n[0070] 或差异平方的和(SSD)。\n[0071] \n[0072] 或差异立方的和(SCD)\n[0073] \n[0074] 此外,差异的幂的和(SPD)\n[0075] \n1 2\n[0076] 也可以采用,Sx,y是在第一帧的x,y处的超声信号,Sx+X,y+Y是在第二帧的x+X,y+Y处的超声信号,q是实数,m和n限定相关窗口尺寸,而 和 为在第一帧和第二帧的窗口中\n的平均值。在一个实施例中,这些等式可以用于找到最小值以追踪超声信号来获取运动(X,\nY)。对于相关性的情况,执行对最大值相关性的寻找以获取运动。如果运动大于预定值,认\n为探头移动,启动自动增益补偿系统103。可选地,这些X=0,Y=0的等式可以用于计算\n两个帧之间的总的相关性以示出静止图像的质量。除了相关性为(1)的情况外,如果这些\n值大于预设值,认为探头移动,启动自动增益补偿系统103。在相关性为(1)的情况下,如果\n相关性小于预设值,认为探头移动,启动自动增益补偿系统103。\n[0077] TGC、轴向、侧向和高度放大器119y,133,135,143的增益补偿的启动可以在增益补偿数据改变时执行。来自轴向引擎155,侧向引擎157或高度引擎159的增益变化数据可以\n由运动处理器165来监测。第一组织区域与另一个相比可以显著地改变超声信号以指示探\n头运动,这里需要新的增益变化数据。如果所探测到的运动大于预定值,可以认为探头105\n运动,启动增益补偿,其中,增益补偿包括使用新的增益补偿曲线设定放大器增益。在可选\n实施例中,只有在探测到运动时,才启动增益变化和补偿计算。新的增益补偿改变施加到放\n大器的任何之前的增益。\n[0078] 还可以使用附着到探头105的运动传感器167来探测运动。运动传感器可以通过\n运动传感器系统166耦合至运动探测处理器165,所述运动传感器例如速度计、加速计、陀\n螺仪或运动跟踪装置或诸如Ascension Technology公司的flock of bird传感器的位置\n传感器。如果运动、速度或加速度大于预定值,认为探头运动,并且执行增益补偿计算。在\n另一个实施例中,可以周期性地激活自动增益补偿系统103,例如以预定的时间间隔。\n[0079] 轴向引擎155、侧向引擎157和高度引擎159可以是DSP、ASIC(特定用途集成电\n路)、FPGA(现场可编程门阵列)、通用处理器、存储器或诸如加法器、乘法器、除法器和其\n它设备的分立设备,或以上部件类型的组合。自动增益补偿系统103可以包括DSP、ASIC、\nFPGA、通用处理器或分立设备,或以上部件类型的组合。自动增益补偿系统103可以接收线\n数据、基带数据或RF信号数据作为图像数据输入。在可选实施例中,系统103还可以使用\n从处理器137输出的超声图像数据,例如规则超声B模式(灰阶)线数据。该线数据可以\n用于计算增益补偿曲线。如果使用束形成的RF数据147,系统103执行探测151以获得信\n号的幅值。探测151可以使用积分探测、混频器(乘法器)、低通滤波、希尔伯特转换或幅值\n探测来执行。\n[0080] 自动增益补偿系统103可以执行图像数据值的数种不同的分析。对于一个图像\n帧,这些分析可以包括轴向增益分析和/或侧向增益分析。可以施加一个或两个增益补偿\n到从随后图像扫描所接收的数据信号。对于三维或体成像,需要预定系列的图像帧。体成\n像可以包括高度体增益分析、轴向体分析和/或侧向体分析,其中,高度体增益分析用于确\n定帧之间或帧的组之间的增益。这三种分析的任意的组合可以应用于下一个所接收的体。\n[0081] 包括本发明的增益补偿方法的软件可以被加载到数据存储163中并在不同的引\n擎155、157、159中执行以实现增益变化和增益补偿曲线的计算。因为增益变化曲线是基于\n实际患者的超声数据计算的,通道TGC放大器119y通常对于补偿之后的随后的图像扫描执\n行得最好,提供最好的动态范围和信号噪声比(SNR)。\n[0082] 为了推导轴向或侧向增益补偿曲线,帧图像数据值ai,j必须根据所期望的分析进\n行分区。\n[0083] 如图2中所示的是示例性图像数据帧201。图3示出了本发明的方法。为了简\n便,示例性帧是从15束组合,N=15,每束具有15个图像值,M=15,总计225个图像值ai,\nj(步骤305、310)。每个图像数据值ai,j可以是多位值,包括例如8位表示255个图像等级\nn\n(2-1,其中,n=8)。\n[0084] 确定下限limitlower和上限limitupper并将其应用到帧中的每个图像数据值ai,j。从血液或其它流体返回的超声信号与从健康软组织返回的信号相比显示出较低的信号幅值。\n因此,可以识别并标记具有比下限小的幅值的信号ai,j≤limitlower。相反,从器官边界或骨骼返回的超声显示出极其高的信号幅值。具有大于上限的幅值的信号ai,j≥limitupper可\n以被识别并标记为无用,因为这些信号不表示软组织。所标记的图像数据值被认为不是有\n效的limitlower
limit1,代入0,或预定的值\n[0098] 如果|ΔAS C-AS B|>limitl,代入0,或预定的值\n[0099] 如果|ΔAS D-AS C|>limit1,代入0,或预定的值\n[0100] 如果|ΔAS E-AS D|>limit1,代入0,或预定的值(6)\n[0101] 如果一个部分的代表信号值和与其相邻的部分的代表信号值之间的绝对差异大\n于预定limit1,将不使用增益变化△值(ΔAS B-AS A,ΔAS C-AS B,ΔAS D-AS C和ΔAS E-AS D,或ΔAS D-AS E ,ΔAS C-AS D,ΔAS B-AS C和ΔAS B-AS A),使用零或一个预定值来代替该增益变化△值。这种替代阻止了在分布中的明显扭曲。在可选的实施例中,可能不发生这种比较和替代过程。\n[0102] 从增益变化△’s计算多个增益变化分布序列值。参考位置部分的代表信号值用\n于分布序列值中的一个。新的分布序列值如果在分布方向是通过将增益变化△加到之前的\n增益变化分布序列值来计算的,或者如果不在分布方向是通过从之前的增益变化分布序列\n值中减去增益变化△来计算的。分布序列值的整个序列是顺序布置的(步骤345)。\n[0103] 使用最浅的轴向部分AS A作为具有从顶部到底部的分布方向的参考位置。\n[0104] \n[0105] IAS A =axial gainvar dist1\n[0106] IAS A+ΔAS B-AS A =axial gainvar dist2\n[0107] IAS A+ΔAS B-AS A+ΔAS C-AS B =axial gainvar dist3 (7)\n[0108] IAS A+ΔAS B-AS A+ΔAS C-AS B+ΔAS D-AS C =axial gainvar dist4[0109] IAS A+ΔAS B-AS A+ΔAS C-AS B+ΔAS B-AS C+ΔAS E-AS D =axial gainvar dist5[0110] 使用最深的轴向部分AS E作为具有从顶部到底部的分布方向的参考位置。\n[0111] \n[0112] IAS E =axial gainvar dist5\n[0113] IAS E+ΔAS D-AS E =axial gainvar dist4\n[0114] IAS E+ΔAS D-AS E+ΔAS C-AS D =axial gainvar dist3 (8)[0115] IAS E+ΔAS D-AS E+ΔAS C-AS D+ΔAS B-AS C =axial gainvar dist2[0116] IAS E+ΔAS D-AS E+ΔAS C-AS D+ΔAS B-AS C+ΔAS A-AS B =axial gainvar dist1[0117] 使用中间的轴向部分AS C作为具有从顶部到底部的分布方向的参考位置。\n[0118] \n[0119] IAS C-ΔAS C-AS B-ΔAS B-AS A =axial gainvar dist1\n[0120] IAS C-ΔAS C-AS B =axial gainvar dist2\n[0121] IAS C =axial gainvar dist3 \n(9)\n[0122] IAS C+ΔAS D-AS C =axial gainvar dist4\n[0123] IAS C+ΔAS D-AS C+ΔAS E-AS D =axial gainvar dist5\n[0124] 以1,2,3,4和5的顺序的分布序列值的整个序列使用诸如线性、三次样条以及其\n它的插值函数组成连续增益变化曲线(步骤350)。轴向增益补偿曲线是通过改变增益变化\n曲线的极性(即,加或减)从增益变化曲线推导出来的。如果该值是在线性标度中,该增益\n补偿曲线是增益变化曲线的倒数(步骤355)。图5中示出了在图像数据帧轴向分成20个\n部分的示例性轴向增益补偿曲线。\n[0125] 轴向增益补偿曲线用于使用由轴向引擎155输出和控制器161分配的增益补偿曲\n线数据经由可变增益TGC放大器119y和可变增益轴向放大器133补偿到达信号。与位置\n(轴向深度)相关的来自轴向补偿曲线的增益的量可以被分配。增益可以在可变增益TGC\n放大器119y和可变增益轴向放大器133之间分配,而不会存在可变增益TGC放大器119y的\n增益潜在的过载。轴向增益补偿可以施加到下一个所接收的图像数据帧(步骤360)。\n[0126] 图6示出了实际的超声图像、分成10个轴向部分。每个部分旁边示出的是其分别\n的柱状图。图7示出了轴向增益补偿之后的同一超声图像。在轴向增益补偿之后,分析同\n一部分,显示所得的柱状图。可以发现,在补偿之后图像具有垂直的均匀强度。在补偿之后\n的柱状图同样示出了部分与部分之间的均匀强度分布。\n[0127] 同样的过程被用于获取侧向增益补偿曲线。为了推导侧向增益补偿曲线,图像数\n据帧被侧向地分区。图8中所示的为在图2中所示的图像帧的示例性的侧向分区。对于本\n例,该帧被分成5个相等侧向部分(LS)F,G,H,I和J,每个部分具有三个独特的列用于二\n维分析(步骤325)。每个侧向部分可以包括多达45个有效图像数据值ai,j。然而,如上所\n述,为了轴向分析,该帧可以被分成任意数目的侧向部分,具有任意数目的有效图像数据值\nai,j,并且可以为各种形状。一些部分可以与其他部分重叠,共享有效图像数据值。\n[0128] 侧向引擎157使用每个部分中的有效图像数据值ai,j为每个限定的侧向部分LS \nF,LS G,LS H,LS I和LS J计算代表信号值(I)。将为每个侧向部分计算的代表信号值与\n相邻侧向部分相比较。该侧向部分LSF,LS G,LS H,LS I和LS J代表信号值之前如果不\n是分贝可以被转换成分贝或者如果不是对数压缩可以转换成对数压缩的信号。该代表信号\n值的计算可以使用任何上面所列的方法来实现。\n[0129] 侧向部分代表信号值之间的增益变化Δ’s基于所期望的分布方向来计算(步骤\n335)。\n[0130] 结合增益变化Δ,将每个侧向部分LS F,LS G,LS H,LS I和LS J的代表信号值\n与相邻侧向部分的代表信号值比较并和预定的limit2比较(步骤340)。\n[0131] \n[0132] 如果|ΔLS F-LS G|>limit2,代入0,或预定的值\n[0133] 如果|ΔLS G-LS H|>limit2,代入0,或预定的值\n[0134] 如果|ΔLS H-LS I|>limit2,代入0,或预定的值\n[0135] 如果|ΔLS I-LS J|>limit2,代入0,或预定的值 (10)\n[0136] 如果一个部分的代表信号值和与其相邻的部分的代表信号值之间的绝对差异大\n于预定limit2,将不使用增益变化△值(ΔLS G-LS F,ΔLS H-LS G,ΔLS I-LS H和ΔLS J-LS I,或ΔLS I-LS J,ΔLS H-LS I,ΔLS G-LS H和ΔLS F-LS G),使用零或一个预定值来代替该增益变化Δ值。这种替代阻止了在分布中的明显扭曲。在可选的实施例中,可能不执行这种比较和替代过程。\n[0137] 侧向部分的增益变化分布是侧向累积序列。该序列一直延续直到所有的侧向增益\n变化Δ值均被考虑(步骤345)。\n[0138] 使用右端的侧向部分LS F作为具有从右到左分布方向的参考位置。\n[0139] \n[0140] ILS F =lat gainvar dist1\n[0141] ILS F+ΔLS G-LS F =lat gainvar dist2\n[0142] ILS F+ΔLS G-LS F+ΔLS H-LS G =lat gainvar dist3 (11)[0143] ILS F+ΔLS G-LS F+ΔLS H-LS G+ΔLS I-LS H =lat gainvar dist4[0144] ILS F+ΔLS G-LS F+ΔLS H-LS G+ΔLS I-LS H+ΔLS J-LS I =lat gainvar dist5[0145] 使用左端的侧向部分LS J作为具有从左到右分布方向的参考位置。\n[0146] \n[0147] ILS J =lat gainvar dist5\n[0148] ILS J+ΔLS I-LS J =lat gainvar dist4\n[0149] ILS J+ΔLS I-LS J+ΔLS H-LS I =lat gainvar dist3(12)[0150] ILS J+ΔLS I-LS J+ΔLS H-LS I+ΔLS G-LS H =lat gainvar dist2[0151] ILS J+ΔLS I-LS J+ΔLS H-LS I+ΔLS G-LS H+ΔLS F-LS G =lat gainvar dist1[0152] 以1,2,3,4和5顺序的分布序列值的整个序列使用插值函数组成连续增益变化曲\n线(步骤350)。侧向增益补偿曲线是通过改变增益变化曲线的极性(即,加或减)从增益\n变化曲线推导出来的。如果该值是在线性标度中,该增益补偿曲线是增益变化曲线的倒数\n(步骤355)。图9中示出了在图像数据帧侧向上分成12个部分的示例性侧向增益补偿曲\n线。\n[0153] 侧向增益补偿曲线用于使用由侧向引擎157和控制器161输出的侧向增益补偿曲\n线数据经由可变增益TGC放大器119y和可变增益侧向放大器135补偿到达信号。来自侧\n向补偿曲线的增益的量可以在放大器之间被分配。使用高速TGC放大器,侧向增益补偿曲\n线可以应用到TGC放大器用于每个所接收到的超声信号,实现侧向增益补偿。通常,该可变\n增益侧向补偿放大器135用于对侧向增益变化进行补偿。增益可以在可变增益TGC放大器\n119y和可变增益侧向放大器135之间分配,而不会过载可变增益TGC放大器119y的增益。\n侧向增益补偿可以施加到一个或多个输入的图像数据帧(步骤360)。\n[0154] 图10示出了分成4个侧向部分的实际的超声图像,每个分别的部分的所得的柱状\n图在旁边显示。图11示出了侧向增益补偿之后的同一超声图像。在侧向增益补偿之后,分\n析同一部分,显示了所得的柱状图。可以发现,在处理之后图像具有均匀的水平强度。柱状\n图同样示出了部分与部分之间的均匀强度分布。\n[0155] 因为高度增益补偿曲线的计算是用于体成像,所以需要一系列帧(步骤320,\n365)。图12中所示是用于一系列五个连续的超声图像帧a,b,c,d和e的示例性高度计\n算。通常,50-100帧或更多可以包括体。下限和上限limitlower,limitupper被应用于每个\n图像数据值limitlowerlimit3,代入0,或预定的值\n[0162] 如果|ΔES2-ES3|>limit3,代入0,或预定的值(13)\n[0163] 如果一个高度部分代表信号值和与其相邻的高度部分代表信号值之间的绝对差\n异大于预定的limite3,将不使用增益变化Δ值(ΔES1-ES2和ΔES2-ES3,或ΔES3-ES2和ΔES2-ES1),使用零或一个预定值来代替该增益变化Δ值。在可选的实施例中,可能不执行这种比较和\n替代过程。\n[0164] 用于高度部分的增益变化分布是高度累积序列。该序列将一直延续直到所有的高\n度部分增益变化Δ值均被考虑(步骤390)。\n[0165] 使用前面的高度部分ES1作为具有从前到后分布方向的参考位置。\n[0166] \n[0167] IES1 =elevation gainvar dist1\n[0168] IES1+ΔES2-ES1 =elevation gainvar dist2 (14)\n[0169] IES1+ΔES2-ES1+ΔES3-ES2 =elevation gainvar dist3\n[0170] 使用后面的高度部分ES3作为具有从后到前分布方向的参考位置。\n[0171] \n[0172] IES3 =elevation gainvar dist3\n[0173] IES3+ΔES2-ES3 =elevation gainvar dist2 (15)\n[0174] IES3+ΔES2-ES3+ΔES1-ES2 =elevation gainvar dist1\n[0175] 以1,2和3顺序的分布序列值的整个序列使用插值函数组成连续增益变化曲线\n(步骤395)。高度增益补偿曲线是通过改变增益变化曲线的极性(即,加或减)从增益变化\n曲线推导出来的。如果该值是在线性标度中,该增益补偿曲线是增益变化曲线的倒数(步\n骤397)。示例性高度补偿曲线如图13中所示。\n[0176] 高度增益补偿曲线用于使用由高度引擎159和控制器161输出的高度增益补偿曲\n线数据使用可变增益TGC放大器119y和可变增益高度放大器143补偿到达信号。高度增\n益补偿可以施加到体中的下一个序列的帧(步骤360)。\n[0177] 图14中所示是5个图像帧的示例性轴向分区。为了轴向体增益补偿,轴向引擎\n155为遍及该体的同一轴向部分计算代表信号值(步骤365,370)。\n[0178] 遍及帧b,c,d和e,轴向体部分(AVS)可以是在帧a中的图像数据值的同一区域。\n该体可以被分成任意数目的轴向部分,这些部分具有相同或不同数目的图像数据值。如前\n所述,以圆柱、椭圆柱或其它形状的形式的取样可以遍及该体作为一个部分使用。此外,一\n些部分可以与其他部分重叠。使用在限定部分内的有效图像数据值。\n[0179] 轴向引擎155为每个限定的轴向体部分AVS A,AVS B,AVSC,AVS D和AVS E计算\n代表信号值(I)(步骤370),并将其与相邻轴向体部分代表信号值比较。该轴向体部分AVS \nA,AVS B,AVS C,AVS D和AVSE的代表信号值如果不是分贝可以转换成分贝,或者转换成对\n数压缩的信号。如上所述,该代表信号值的分析可以使用数个不同方法来执行。\n[0180] 相邻轴向体部分之间的增益变化Δ’s基于所期望的分布方向来计算(步骤375)。\n[0181] 结合增益变化Δ,将每个轴向体部分AVS A,AVS B,AVS C,AVSD和AVS E的代表\n信号值与相邻轴向体部分的代表信号值比较并和预定的limit4比较(步骤385)。\n[0182] \n[0183] 如果|ΔAVS A-AVS B|>limit4,代入0,或预定的值\n[0184] 如果|ΔAVS B-AVS C|>limit4,代入0,或预定的值\n[0185] 如果|ΔAVS C-AVS D|>limit4,代入0,或预定的值\n[0186] 如果|ΔAVS D-AVS E|>limit4,代入0,或预定的值(16)\n[0187] 如果一个轴向体部分的代表信号值和与其相邻的轴向体部分的代表信号值之间\n的绝对差异大于预定limit4,将不使用增益变化Δ值(ΔAVS B-AVS A,ΔAVS C-AVS B,ΔAVS D-AVS C和ΔAVS E-AVS D,或ΔAVS A-AVS B,ΔAVS B-AVS C,ΔAVS C-AVS D,ΔAVS D-AVS E),使用零或一个预定值来代替该增益变化Δ值。这种替代阻止了在分布中的明显扭曲。在可选地实施例中,可能不执行这\n种比较和替代过程。\n[0188] 轴向体部分的增益变化分布是轴向累积序列。该序列一直延续直到所有的轴向体\n增益变化Δ值均被考虑(步骤390)。\n[0189] 使用最浅的轴向体部分AVS A作为具有从顶部到底部分布方向的参考位置。\n[0190] \n[0191] IAVS A =axvol gainvar \ndistl\n[0192] IAVS A+ΔAVS B-AVS A =axvol gaindist var2[0193] IAVS A+ΔAVS B-AVS A+ΔAVS C-AVS B =axvolgaindist var3(17)\n[0194] IAVS A+ΔAVS B-AVS A+ΔAVS C-AVS B+ΔAVS D-AVS C =axvol gaindist var4[0195] IAVS A+ΔAVS B-AVS A+ΔAVS C-AVS B+ΔAVS D-AVS C+ΔAVS E-AVS D =axvol gaindist var5[0196] 使用最深的轴向体部分AVS E作为具有从底部到顶部分布方向的参考位置。\n[0197] \n[0198] IAVS E =ax vol gainvar \ndist5\n[0199] IAVS E+ΔAVS D-AVS E =ax vol gainvar dist4\n[0200] IAVS E+ΔAVS D-AVS E+ΔAVS C-AVS D =ax vol gainvar dist3 (18)\n[0201] IAVS E+ΔAVS D-AVS E+ΔAVS C-AVS D+ΔAVS B-AVS C =ax vol gainvar dist2[0202] IAVS E+ΔAVS D-AVS E+ΔAVS C-AVS D+ΔAVS B-AVS C+ΔAVS A-AVS B =ax vol gainvar dist1[0203] 以1,2,3,4和5顺序的分布序列值的整个序列使用插值函数组成连续轴向体增益\n变化曲线(步骤395)。轴向体增益补偿曲线是通过改变轴向体增益变化曲线的极性(即,\n加或减)从轴向体增益变化曲线计算的(步骤397)。如果该值是在线性标度中,该增益补\n偿曲线是增益变化曲线的倒数。图15中示出了示例性的轴向体补偿曲线。\n[0204] 轴向体增益补偿曲线用于使用由轴向引擎155输出和控制器161分配的轴向体增\n益补偿曲线数据经由可变增益TGC放大器119y和可变增益轴向放大器133补偿到达信号。\n轴向体补偿可以施加到下一个体(步骤360)。\n[0205] 图16中所示为5个图像数据帧体的示例性侧向体分区。为了侧向体增益补偿,侧\n向引擎157为遍及该体的同一侧向体部分计算代表信号值(步骤365)。\n[0206] 遍及帧b,c,d和e,侧向体部分(LVS)可以是在帧a中的图像数据值的同一区域。\n该体可以被分成任意数目的侧向部分,这些部分具有相同或不同数目的图像数据值。如轴\n向体分析所述,以圆柱或其它形状的形式的取样可以遍及该体作为一部分使用。一些部分\n可以与其他部分重叠。使用在限定部分内的有效图像数据值。\n[0207] 侧向引擎157为每个限定的侧向体部分LVS F,LVS G,LVSH,LVS I和LVS J计算\n代表信号值(I),并将其与相邻侧向体部分代表信号值比较。该侧向体部分LVS F,LVS G,\nLVS H,LVS I和LVS J的代表信号值如果不是分贝可以转换成分贝或者转换成对数压缩的\n信号。该代表信号值的计算可以使用上面列出的任何方法来执行。\n[0208] 相邻侧向体部分之间的增益变化Δ’s基于所期望的方向来计算(步骤380)。\n[0209] 结合增益变化Δ,将每个侧向体部分LVS F,LVS G,LVS H,LVSI和LVS J的代表\n信号值与相邻侧向体部分的代表信号值比较并和预定的limit5比较(步骤385)。\n[0210] \n[0211] 如果|ΔLVS F-LVS G|>limit5,代入0,或预定的值\n[0212] 如果|ΔLVS G-LVS H|>limit5,代入0,或预定的值\n[0213] 如果|ΔLVS H-LVS I|>limit5,代入0,或预定的值\n[0214] |ΔLVS I-LVS J|>limit5,\n[0215] 如果,代入0,或预定的值(19)\n[0216] 如果一个侧向体部分的代表信号值和与其相邻的侧向体部分的代表信号值之间\n的绝对差异大于预定limit5,将不使用增益变化Δ值(ΔLVS F-LVS G,ΔLVS G-LVS H,ΔLVS H-LVS I和ΔLVS I-LVS J,或ΔLVS G-LVS F,ΔLVS H-LVS G,ΔLVS I-LVS H和ΔLVS J-LVS I),使用零或一个预定值来代替该增益变化Δ值。这种替代阻止了在分布中的明显扭曲。在可选的实施例中,可能不执行\n这种比较和替代过程。\n[0217] 侧向体部分之间的增益变化分布是侧向累积序列。该序列一直延续直到所有的侧\n向体增益变化Δ值均被考虑(步骤397)。\n[0218] 使用右端的侧向体部分LVS F作为具有从右到左分布方向的参考位置。\n[0219] \n[0220] ILVS F =latvol gainvar \ndist1\n[0221] ILVS F+ΔLVS G-LVS F =latvol gainvar dist2\n[0222] ILVS F+ΔLVS G-LVS F+ΔLVS H-LVS G =latvol gainvar dist3(20)\n[0223] ILVS F+ΔLVS G-LVS F+ΔLVS H-LVS G+ΔLVS I-LVS H =latvol gainvar dist4[0224] ILVS F+ΔLVS G-LVS F+ΔLVS H-LVS G+ΔLVS I-LVS H+ΔLVS J-LVS I =latvol gainvar dist5[0225] 使用左端的侧向体部分LVS J作为具有从左到右分布方向的参考位置。\n[0226] \n[0227] ILVS J =latvol gainvar \ndist5\n[0228] ILVS J+ΔLVS I-LVS J =latvol gainvar dist4\n[0229] ILVS J+ΔLVS I-LVS J+ΔLVS H-LVS I =latvol gainvar dist3\n[0230] (21)\n[0231] ILVS J+ΔLVS I-LVS J+ΔLVS H-LVS I+ΔLVS G-LVS H =latvol gainvar dist2[0232] ILVS J+ΔLVS I-LVS J+ΔLVS H-LVS I+ΔLVS G-LVS H+ΔLVS F-LVS G =latvol gainvar dist1[0233] 以1,2,3,4和5顺序的分布序列值的整个序列使用插值函数组成连续侧向体增益\n变化曲线(步骤395)。侧向体增益补偿曲线是通过改变增益变化曲线的极性(即,加或减)\n从增益变化曲线推导出来的(步骤397)。如果该值是在线性标度中,该增益补偿曲线是增\n益变化曲线的倒数。图17中示出了示例性的侧向体增益变化补偿曲线。\n[0234] 侧向体增益补偿曲线用于使用由侧向引擎157输出和控制器161分配的侧向体增\n益补偿曲线数据经由可变增益TGC放大器119y和可变增益侧向放大器135侧向补偿到达\n信号。侧向补偿可以施加到下一个体(步骤355)。\n[0235] 如在二维轴向和/或侧向增益补偿曲线的情况下,高度和轴向和/或侧向体分区\n的结合可以得到三个体增益补偿曲线。每个曲线可以施加到下一个体(步骤335)。\n[0236] 高度、轴向和侧向增益补偿的任意结合都可以应用到下一个输入的图像帧。该新\n计算的增益补偿曲线示出了与之前放大器的增益曲线的差别,其中,之前放大器的增益曲\n线提供输出以控制每个放大器119y,133,135,143的增益。推导的该增益补偿曲线修改预\n先存在的增益曲线。\n[0237] 回到图1,来自轴向引擎155、侧向引擎157和高度引擎159的输出被输入到控制\n器161,其中,增益补偿曲线在控制器161处被处理。\n[0238] 因为新的增益补偿曲线升级或修改先前计算的增益补偿曲线,每个控制器161的\n轴向、侧向和高度输入耦合至相应的增益加法器169a,b,c。每个加法器169a,b,c的输出\n耦合至相应的延迟171a,b,c,延迟171a,b,c被耦合至存储器170a,b,c。存储器170a,b,c存储在初始化时的之前的增益补偿曲线,及其之后升级的结合增益。增益加法器169a,b,\nc输出增益,该增益是新计算的增益曲线和之前或升级的增益的结合。\n[0239] 轴向延迟171a用于延迟TGC加上轴向放大器增益曲线的总和(即,TGC与轴向放\n大器增益的结合,并且称之为轴向结合增益)一个帧,从而使得之前的帧的轴向结合增益\n曲线存储在存储器170a中,之前的帧的轴向结合增益曲线加到新的轴向补偿曲线以获取\n新的轴向结合增益曲线。\n[0240] 类似的,侧向延迟171b用于延迟TGC加上侧向放大器增益曲线的总和(即,TGC与\n侧向放大器增益的结合,并且称之为侧向结合增益)一个帧,从而使得之前的帧的侧向结\n合增益曲线存储在存储器170b中,之前的帧的侧向结合增益曲线加到新的侧向补偿曲线\n以获取新的侧向结合增益曲线。\n[0241] 高度延迟171c用于延迟TGC加上高度放大器增益曲线的总和(即,TGC与高度放\n大器增益的结合,并且称之为高度结合增益)一个体,从而使得之前的体的高度结合增益\n曲线存储在存储器170c中,之前的体的高度结合增益曲线加到新的高度补偿曲线以获取\n新的高度结合增益曲线。·\n[0242] 增益加法器169a,b,c的输出被耦合至第一权重乘法器173a,b,c和第二权重乘\n法器175a,b,c。该第一权重乘法器173a,b,c被耦合至分配权重α(轴向),β(侧向)和\nγ(高度),而第二权重乘法器175a,b,c被耦合至分配权重1-α177a,β177b和1-γ177c。\n[0243] 权重α,β和γ确定在时间t在轴向可变增益放大器133、侧向可变增益放大器\n135和高度可变增益放大器143,以及TGC放大器119y的可变增益部分之间分配的各自的\n增益补偿值是多少。权重α,β和γ使用优化每个权重α,β和γ的分配处理器179来\n推导。确定权重α,β和γ的优化过程的讨论不在本发明内容的范围内。\n[0244] 来自第二权重乘法器175a,b,c的输出被输出至相应的轴向可变增益放大器133、\n侧向可变增益放大器135和高度可变增益放大器143。\n[0245] 轴向权重乘法器173a的输出被耦合至加法器179和平均值计算器181。平均值计\n算器181获取来自轴向乘法器173a的输出的平均值,该平均值被从轴向权重乘法器173a\n输出中减去并耦合至公用加法器183。加法器179从TGC增益曲线中减去TGC增益曲线平\n均值,使得恒定的增益施加到TGC放大器119y的固定增益部分,并且只有曲线增益施加到\nTGC放大器119y的可变增益部分。本操作保持TGC放大器的动态范围。在可选实施例中,\n可以不使用平均值和减法。\n[0246] 轴向结合增益曲线乘以α以便将增益曲线的一部分施加到TGC放大器119y可变\n增益部分。加法器179的输出被耦合至DAC(数字到模拟转换器)185,如果需要,并输出至\nTGC放大器119y可变增益部分。\n[0247] 每束更新来自存储器170b的侧向结合增益以实现侧向增益曲线的效用,并乘以\nβ以便将部分增益曲线施加到TGC放大器119y的可变增益。侧向延迟187被用于同步所\n有的增益曲线。\n[0248] 每帧更新来自存储器170c的高度结合增益以实现高度增益曲线的效用,并乘以\nγ以便将部分增益曲线施加到TGC放大器119y的可变增益。高度延迟189被用于同步所\n有的增益曲线。\n[0249] 延迟187、189的输出和平均值计算器181的输出被耦合至公用加法器183。加法\n器183结合所有TGC放大器119y的固定增益部分。所结合的增益被输出至另一个加法器\n191并通过增益计算器193与总增益输出结合。\n[0250] 对于三维成像的情况,侧向放大器135或高度放大器143的输出被耦合至总增益\n计算器193以通过使用柱状图、百分数、平均值、众数、中间数来获得总增益。以dB或对数\n标度形式的总增益被改变极性(减或加)并加上一个预设值以设定最终图像亮度到达预设\n值。该操作设定总的图像亮度到大部分用户喜欢的预设亮度值。\n[0251] 从加法器191得到的结合的输出是TGC放大器119y的固定增益部分。该TGC放\n大器119y的固定增益部分耦合至在TGC放大器119y芯片中的固定的或可编程的增益放大\n器。\n[0252] 轴向、侧向和高度增益变化由上面的操作进行补偿,并且总的图像增益被设定至\n用户可能喜欢的预定图像亮度。本发明提供所有图像维数上的均匀性。\n[0253] 上面描述了本发明的一个或多个实施例。然而,应该理解,在不脱离本发明精神和\n范围的情况下可以进行多种修改。因此,其它实施例也在下面权利要求的范围内。法律信息
- 2022-03-15
未缴年费专利权终止
IPC(主分类): A61B 8/00
专利号: ZL 200780011611.4
申请日: 2007.04.02
授权公告日: 2011.09.14
- 2017-02-22
专利权的转移
登记生效日: 2017.02.03
专利权人由日立阿洛卡医疗株式会社变更为株式会社日立制作所
地址由日本东京变更为日本东京
- 2012-03-28
专利权人的姓名或者名称、地址的变更
专利权人由阿洛卡株式会社变更为日立阿洛卡医疗株式会社
地址由日本东京变更为日本东京
- 2011-09-14
- 2009-07-29
- 2009-06-03
引用专利(该专利引用了哪些专利)
序号 | 公开(公告)号 | 公开(公告)日 | 申请日 | 专利名称 | 申请人 |
1
| | 暂无 |
2004-05-31
| | |
2
| |
2005-07-27
|
2003-04-01
| | |
3
| | 暂无 |
1995-03-21
| | |
被引用专利(该专利被哪些专利引用)
序号 | 公开(公告)号 | 公开(公告)日 | 申请日 | 专利名称 | 申请人 | 该专利没有被任何外部专利所引用! |