r>[0055] 进一步，如图3所示，所述步骤S104将所述求取的所有灰度值差与阈值T相比较； 将灰度值差大于阈值T的点设置为造影成像区域，包括如下步骤：
[0056] S1041，将所述求取的灰度值差与所述阈值T相比较，若大于灰度值差阈值T，将该 点作为有效点记录。
[0057] 用后一卷造影成像数据，减去前一卷造影成像数据对应位置点的灰度值，若 大于阈值T，则将该点的位置记录下来。记当前卷造影成像数据被记录下的点集为 ValidRegioni?
[0058]ValidRegioni={pointcurVol (i，j，k)GV〇lcur | (p〇intcurVol (i，j，k)-pointpreVol (i， j，k)) >T}
[0059]point^^Q，j，k)为两卷数据中的后一卷坐标为（i，j，k)的点；pointpreVcil(i，j， k)为两卷数据中前一卷坐标为（i，j，k)的点，ValidRegion$根据上文方法判断出的造影 剂灌注区域点集；
[0060] S1042,顺序遍历所述采集的N卷造影成像数据，每相邻两卷求出一个点集 ValidRegionp直到最后一卷的有效点集ValidRegionn丨被求取出来。
[0061] S1043,将所述有效点并集的补集设为非造影区域。
[0062] 设第2卷被记录的点集为ValidRegior^，第3卷为ValidRegion2,......，第n卷 被记录的点集为ValidRegionni。将背景点区域设为：
[0063]
[0064] point(i，j，k)point(i，j，k)为所采集的超声图像坐标为（i，j，k)的点； ValidRegionValidRegion为根据上文方法判断出的造影剂灌注区域点集；BKRegion BKRegion为背景区域点集。
[0065] 即非造影区域为有效区域并集的补集；
[0066] 进一步，所述设置灰度值差阈值T的方法包括，统计造影剂灌注过程中相同位置 各点相邻两卷超声造影数据之间的平均灰度值差，将该灰度值差作为阈值T，所述平均灰度 差求取方法可以采用如下两个公式求取：
[0067]
[0068]
[0069] 其中，Vdi，j，k)为相邻两卷数据中的前一卷坐标为（i，j，k)点的灰度值；V2(i， j，k)为相邻两卷数据中的后一卷数据坐标为（i，j，k)点的灰度值。frms为一卷数据中的 总帧数，rows为一卷数据中单帧数据的行数；cols为一卷数据中单帧数据的列数。T为所 求阈值，M为相邻两卷数据的平均灰度值。
[0070] 实施例二、
[0071] 如图4所示，本发明具体实施例还提供一种超声造影成像非造影区域自动去除的 装置300。所述装置包括：
[0072] 采集单元301，用于采集灌注造影剂过程中的N卷超声造影成像数据。
[0073] 第一计算单元302,用于相邻两卷的超声造影成像数据的灰度值差的阈值T;
[0074] 第二计算单元303,用于求所述N卷超声造影成像数据中相邻两卷所有相同位置 的超声造影成像数据的灰度值差；
[0075] 第三计算单元304,用于将所述求的灰度值差与所述阈值T相比较，将所述灰度值 差大于所述阈值T的点设置为造影成像区域，则剩余部分为非造影成像区域；
[0076] 控制单元305,用于去除非造影成像区域。
[0077] 进一步，如图5所示，所述装置400还包括第四计算406单元，用于对采集的超声 数据进行降噪处理。
[0078] 进一步，如图6所示，所述第三计算单元304包括：第一计算装置3041、第二计算 装置3042、第三计算装置3043。
[0079] 所述第一计算装置，用于将所述求的灰度值差与所述阈值T相比较，将所述灰度 值差大于所述阈值T的点作为造影区域点记录；
[0080] 所述第二计算装置，用于顺序遍历所述采集的N卷造影成像数据，每相邻两卷求 出一个点集，直到最后一卷的有效点集被求取出来；
[0081] 所述第三计算装置，用于将所述有效点并集的补集设为背景点区域。
[0082] 所述各装置的工作过程参见实施例一，在此不再赘述。
[0083] 实施例三、
[0084] 为本发明实施例所述的一种设备，所述设备为超声设备，所述超声设备包括如任 意实施例二所述的装置，所述装置参见具体实施例二，在此不再赘述。
[0085] 以上对本发明所提供的一种超声造影成像非造影区域自动去除方法、装置及设备 进行了详细介绍，对于本领域的一般技术人员，依据本发明实施例的思想，在【具体实施方式】 及应用范围上均会有改变之处，因此，本说明书内容不应理解为对本发明的限制。
【主权项】
1. 一种超声造影成像非造影区域自动去除方法，其特征在于，所述方法包括： 采集灌注造影剂过程中的N卷超声造影成像数据，N多2 ; 设置相邻两卷的超声造影成像数据的灰度值差的阈值T ; 求所述N卷超声造影成像数据中相邻两卷所有相同位置的超声造影成像数据的灰度 值差； 将所述求的灰度值差与所述阈值T相比较，将所述灰度值差大于所述阈值T的点设置 为造影成像区域，则剩余部分为非造影成像区域； 去除非造影成像区域。2. 根据权利要求1所述的方法，其特征在于，所述设置相邻两卷的超声造影成像数据 的灰度值差阈值T之前还包括，对所述采集的超声造影成像数据进行降噪处理。3. 根据权利要求1或2所述的方法，其特征在于，所述将所述求取的灰度值差与所述阈 值T相比较，将所述灰度值差大于所述阈值T的点设置为造影成像区域，则剩余部分为非造 影成像区域包括： 将所述求取的灰度值差与所述阈值T相比较，将所述灰度值差大于所述阈值T的点作 为造影区域点记录； 顺序遍历所述采集的N卷造影成像数据，每相邻两卷求出一个点集，直到最后一卷的 有效点集被求取出来； 将所述有效点并集的补集设为非造影成像区域。4. 根据权利要求1或2所述的方法，其特征在于，所述设置灰度值差阈值T的方法包 括，统计造影剂灌注过程中相同采图位置的相邻两卷超声造影成像数据之间的平均灰度值 差，将所述平均灰度值差作为当前阈值T。5. -种超声造影成像非造影区域自动去除的装置，所述装置包括： 采集单元，用于采集灌注造影剂过程中的N卷超声造影成像数据。 第一计算单元，用于相邻两卷的超声造影成像数据的灰度值差的阈值T ; 第二计算单元，用于求所述N卷超声造影成像数据中相邻两卷所有相同位置的超声造 影成像数据的灰度值差； 第三计算单元，用于将所述求的灰度值差与所述阈值T相比较，将所述灰度值差大于 所述阈值T的点设置为造影成像区域，则剩余部分为非造影成像区域； 控制单元，用于去除非造影成像区域。6. 根据权利要求5所述的装置，其特征在于，所述装置还包括第四计算单元，用于对所 述采集的超声造影成像数据进行降噪处理。7. 根据权利要求5或6所述的装置，其特征在于，所述第三计算单元包括：第一计算装 置、第二计算装置；第三计算装置 所述第一计算装置，将所述求取的灰度值差与所述阈值T相比较，将所述灰度值差大 于所述阈值T的点作为造影区域点记录； 所述第二计算装置，用于顺序遍历所述采集的N卷造影成像数据，每相邻两卷求出一 个点集，直到最后一卷的有效点集被求取出来； 所述第三计算装置，用于将所述有效点并集的补集设为背景点区域。8. -种设备，其特征在于，所述设备包括如权利要求5-7任意一项所述的装置。
【专利摘要】本发明提出一种超声造影成像非造影区域自动去除方法。所述方法包括：采集灌注造影剂过程中的N卷超声造影成像数据，N≥2；设置相邻两卷的超声造影成像数据的灰度值差阈值T；求所述N卷超声造影成像数据中相邻两卷所有相同位置的超声造影成像数据的灰度值差；将所述求的灰度值差与所述阈值T相比较，将所述灰度值差大于所述阈值T的点设置为造影成像区域，通过求取N卷造影成像区域并集的补集来估计背景区域并将其去除。本发明还提供相应的装置及设备。采用本发明根据灌注造影剂过程中相同位置的相邻两卷数据的灰度值差估计背景区域，对超声图像每一卷数据去除背景区域，从而免去超声造影成像中，医生手动去除背景区域操作的麻烦。
【IPC分类】A61B8/00
【公开号】CN105105785
【申请号】CN201510439790
【发明人】黄柳倩, 孙慧, 刘旭江
【申请人】深圳开立生物医疗科技股份有限公司
【公开日】2015年12月2日
【申请日】2015年7月23日