CN102596050B - 超声波成像装置及超声波成像方法 - Google Patents
超声波成像装置及超声波成像方法 Download PDFInfo
- Publication number
- CN102596050B CN102596050B CN201080046798.3A CN201080046798A CN102596050B CN 102596050 B CN102596050 B CN 102596050B CN 201080046798 A CN201080046798 A CN 201080046798A CN 102596050 B CN102596050 B CN 102596050B
- Authority
- CN
- China
- Prior art keywords
- region
- similar degree
- ultrasonic imaging
- handling part
- statistic
- 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.)
- Expired - Fee Related
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0833—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving detecting or locating foreign bodies or organic structures
- A61B8/085—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving detecting or locating foreign bodies or organic structures for locating body or organic structures, e.g. tumours, calculi, blood vessels, nodules
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5269—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
- G06T7/0014—Biomedical image inspection using an image reference approach
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/20—Image preprocessing
- G06V10/25—Determination of region of interest [ROI] or a volume of interest [VOI]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10132—Ultrasound image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30096—Tumor; Lesion
Abstract
本发明提供一种超声波成像装置、超声波成像方法及超声波成像用程序。从处理接收信号而得到的2帧以上的图像中选择基准帧和比较帧。在基准帧中设定关心区域,在比较帧中设定比关心区域宽的探索区域,在探索区域内设定作为所述关心区域的移动目标的候补的多个候补区域,按每个候补区域计算关心区域内和候补区域内的图像特征值的类似度,求出整个探索区域内的类似度的分布。由此,可基于类似度的分布来判别关心区域是否为噪声区域。例如,采用在类似度的分布中求出将类似度的最小值和类似度的整体值进行比较的统计量的构成。通过比较统计量和阈值,可判定关心区域的可靠度。从而,提供可判别回波信号微弱的噪声区域的超声波成像装置。
Description
技术领域
本发明是涉及在通过超声波进行活体成像时可清楚地识别组织边界的超声波摄像方法及超声波摄像装置的技术。
背景技术
在用于医疗图像诊断的超声波摄像装置中,例如,如专利文献1的记载所示,公知下述方法:基于诊断动态图像(B模式图像)的小区域的变化量来推定组织的弹性系数分布,将硬度变换成彩色图进行显示。然而,例如在肿瘤的边缘部,有时相对于周围的组织而言声阻抗和弹性模量都没有太大不同,该情况下,无论在诊断动态图像中还是在弹性图像中都无法掌握肿瘤和周围的组织之间的边界。
因此,在专利文献2记载的技术中提出了下述方法:通过根据诊断动态图像的运动向量直接作成标量场图像,由此可识别相对于周围而言声阻抗和弹性模量都没有太大不同的组织的边界。
现有技术文献
专利文献
专利文献1:日本特开2004-135929号公报
专利文献2:日本特开2008-79792号公报
发明内容
(发明要解决的课题)
在现有的专利文献2所记载的技术中,通过对2个诊断图像数据进行块匹配处理来进行运动向量的推定,但是在进行推定处理时因图像数据的噪声的影像会产生错误向量。因而,存在着边界的识别度劣化的问题。尤其是,在回波信号微弱的穿透(penetration)界限区域中向量的推定精度会发生较大劣化。
本发明的目的在于提供一种可判别回波信号微弱的噪声区域的超声波成像装置。
(用于解决课题的手段)
为了达成上述目的,根据本发明的第1方式,提供一种如下的超声波成像装置。即,本发明的超声波成像装置具有:发送部,其向对象发送超声波;接收部,其接收来自对象的超声波;和处理部,其处理接收部的接收信号,生成2帧以上的图像。处理部将生成的2帧以上的图像中的1帧作为基准帧,在预先确定的位置或从操作者接受的位置设定关心区域。将另外的1帧作为比较帧,在预先确定的位置或从操作者接受的位置设定比关心区域宽的探索区域,在探索区域内设定作为所述关心区域的移动目标的候补的多个候补区域。按每个候补区域计算关心区域内和候补区域内的图像特征值的类似度,求出整个探索区域内的类似度的分布。由此,可基于类似度的分布来判别关心区域是否为噪声区域。
例如采用下述构成:处理部在类似度的分布中求出将类似度的最小值和类似度的整体值进行比较的统计量,根据该统计量判定关心区域的可靠度。具体而言,例如处理部利用类似度的最小值、平均值及标准偏差求出上述统计量,通过比较求得的统计量和阈值可判定关心区域的可靠度。
例如,处理部能够生成将比较帧中与关心区域对应的位置和类似度为最小的候补区域的位置进行连结的向量,针对被判定为可靠度低的关心区域,将向量置换为零或者规定的向量。由此,能够去除错误向量等,从而能够提高向量的精度。
例如,处理部能够针对类似度分布而算出类似度的平均值、最小值及标准偏差,将平均值与最小值之间的差分除以标准偏差而得到的分离度作为统计量来求出。另外,例如也可针对类似度分布而算出类似度的平均值及标准偏差,将标准偏差除以平均值而得到的变动系数用作统计量。
与上述的统计量进行比较的阈值可如下求出。例如,设定多个关心区域,并针对每个区域求出统计量,求出表示求得的统计量的值的频度的直方图分布,作为阈值,使用直方图分布中的中位数、平均值、或者直方图分布呈多个山形状时山之间的谷的最小值的统计量。
另外,可以生成将类似度分布进行了平滑化处理的平滑化后类似度分 布,求出从平滑化处理前的所述类似度分布中减去平滑化后类似度分布后的差分类似度分布。由此,能够从类似度分布中去除因噪声引起的类似度的变动。
上述平滑化处理例如采用下述方法:对于类似度分布设定预先规定的大小的滤波器,并在使滤波器移动规定量的同时反复进行将该滤波器内的分布平滑化的处理。滤波器的大小如下进行决定。针对多个关心区域而预先生成将比较帧中与关心区域对应的位置和在平滑化处理前的类似度分布中类似度为最小的候补区域的位置进行连结的向量。将生成的向量之中最大的向量长度作为滤波器的大小。
也可用拉普拉斯滤波器对上述类似度分布进行滤波处理来作成轮廓强调分布,通过提取轮廓强调分布中连续的轮廓线来求出对象的边界。
处理部通过采用如下构成,能够判定肿瘤的浸润度。即,处理部采用下述构成,具有:第1处理单元,其针对在活体肿瘤边界附近设定的关心区域而求出类似度分布,生成以该类似度作为图像特征值的类似度分布图像,在该类似度分布图像上以关心区域所对应的位置为中心在多个不同方向上设定规定长度的一维区域;第2处理单元,其按照设定的每个方向计算一维区域内的类似度的总和;第3处理单元,其计算类似度总和最小的方向上的类似度总和与正交于该方向的方向上的一维区域的类似度总和之间的比率;和第4处理单元,其基于比率来判定肿瘤的浸润度。
在第3处理单元计算出的比率小于预先设定的一定值的情况下,判断出该关注像素是构成边界线的点,由此可以求出边界线。
另外,根据本发明的第2方式,提供一种如下的超声波成像方法。即,在超声波成像方法中,向对象发送超声波,处理接收来自对象的超声波而得到的接收信号,并生成2帧以上的图像,从图像中选择基准帧和比较帧,在基准帧中设定关心区域,在比较帧中设定比所述关心区域宽的探索区域,在探索区域内设定作为所述关心区域的移动目标的候补的多个候补区域,按每个候补区域计算关心区域内和候补区域内的图像特征值的类似度,求出整个探索区域内的类似度的分布。
另外,根据本发明的第3方式,提供一种如下的超声波成像用程序。即,本发明的超声波成像用程序用于让计算机执行下述步骤:第1步骤, 从2帧以上的超声波图像中选择基准帧和比较帧;第2步骤,在基准帧中设定关心区域;第3步骤,在比较帧中设定比关心区域宽的探索区域,在探索区域内设定作为关心区域的移动目标的候补的多个候补区域;和第4步骤,按每个候补区域计算关心区域内和候补区域内的图像特征值的类似度,求出整个探索区域内的类似度的分布。
(发明效果)
根据本发明,能够基于类似度分布来判断关心区域是否为噪声多的区域。由此,可以抑制错误向量的发生,且即便在穿透界限区域中也可进行高精度的向量推定。能够提高根据推定的运动向量场而变换的标量场图像的精度,从而能够进行更适当的边界检测。
附图说明
图1是表示实施方式1的超声波成像装置的系统构成例的框图。
图2是表示实施方式1的超声波成像装置涉及的图像生成的处理顺序的流程图。
图3是表示图2的步骤24的块匹配处理的详细的流程图。
图4是以2层构造的仿真图(phantom)来说明图2的步骤24的块匹配处理的图。
图5(a)是表示由实施方式1的超声波成像装置生成的B模式图像例的图,图5(b)是表示由实施方式1的超声波成像装置生成的运动向量分布图像例的图。
图6(a)是表示在图5(b)的位置(3)处设定ROI而得到的SAD分布图像例的图,图6(b)是表示在图5(b)的位置(5)处设定ROI而得到的SAD分布图像例的图,图6(c)是图6(a)所示的SAD值的直方图,图6(d)是图6(b)所示的SAD值的直方图。
图7是在SAD值的直方图上表示实施方式1中的分离度的定义的说明图。
图8是表示在利用分离度进行图2的步骤25的处理的情况下的处理的详细的流程图。
图9是表示将通过图8的处理求得的分离度的分布进行图像化的例子 的图。
图10是表示在利用变动系数进行图2的步骤25的处理的情况下的处理的详细的流程图。
图11是表示将通过图10的处理求得的变动系数的分布进行图像化的例子的图。
图12(a)是表示通过图2的步骤24生成的运动向量分布图像例的图,图12(b)是表示通过图2的步骤25去除了错误向量后的向量分布图像例的图。
图13是在运算实施方式2的SAD分布时去除噪声的处理的流程图。
图14(a)是表示实施方式2的噪声去除前的SAD分布图像的例子的图,图14(b)是表示将SAD分布图像进行平滑化(LPF)处理而得到的SAD分布图像的例子的图,图14(c)是表示噪声去除后的SAD分布图像的例子的图。
图15是表示求出实施方式3的浸润度的处理的流程图。
图16(a)~(h)是表示在图15的处理中用到的SAD分布的区域选择图案的说明图。
图17(a)是表示在图5(b)的位置(1)处设定ROI而得到的SAD分布图像的图,图17(b)是表示在图5(b)的位置(2)处设定ROI而得到的SAD分布图像的图,图17(c)是表示对图17(a)所示的SAD值分布应用拉普拉斯滤波器后的图像的图,图17(d)是表示对图17(b)所示的SAD值分布应用拉普拉斯滤波器后的图像的图。
具体实施方式
下面,说明本发明的一实施方式的超声波成像装置。
(实施方式1)
图1示出本实施方式的超声波成像装置的系统构成。本装置具备超声波边界检测功能。如图1所示,本装置构成为具备:超声波探头(探测器)1、用户接口2、送波波束成形器3、控制系统4、送受切换开关5、受波波束成形器6、包络线检波部7、扫描变换器8、处理部10、参数设定部11、合成部12及显示部13。
超声波元件被一维排列而成的超声波探头1,向活体发送超声波波束(超声波脉冲),并接收被活体反射的回波信号(受波信号)。在控制系统4的控制下,对焦于送波焦点的具有延迟时间的送波信号被送波波束成形器3输出,经由送受切换开关5而送至超声波探头1。在活体内反射或者散射并返回到超声波探头1的超声波波束,通过超声波探头1被变换成电信号,并经由送受切换开关5而作为受波信号被送至受波波束成形器6。
受波波束成形器6是将相位相差90度的2个受波信号进行混合的复波束成形器,在控制系统4的控制下,进行根据接收定时来调整延迟时间的动态聚焦,并输出实部和虚部的RF信号。该RF信号被包络线检波部7检波之后被变换成影像信号,并将其输入至扫描变换器8而变换成图像数据(B模式图像数据)。以上说明的构成与公知的超声波成像装置的构成相同。
在本发明的装置中,通过处理部10实现超声波边界检测处理。处理部10具有CPU10a和存储器10b,CPU10a执行预先存储在存储器10b中的程序,由此进行以下的处理来检测被检者组织的边界。即,基于从扫描变换器8输出的2帧以上的图像数据,在处理部10中首先作成运动向量场。接着,将作成的运动向量场变换处理为标量场。然后,在通过合成部12对原始图像数据和与该原始图像数据对应的运动向量场或者标量场进行合成处理之后,由显示部13进行显示。
在参数设定部11中,进行处理部10中的信号处理用的参数、合成部12中的显示图像的选择设定等。这些参数是由操作员(装置操作者)从用户接口2输入的。作为信号处理用的参数,例如可以从操作员接受期望帧m上的关心区域的设定、不同于帧m的帧m+Δ上的探索区域的设定。作为显示图像的选择设定,例如可以从操作者接受将原始图像和向量场图像(或标量图像)合成为1个图像而在显示器中进行显示、还是并排显示2个图像以上的动态图像这样的选择设定。
图2示出本发明的处理部10及合成部12中的边界检测处理及图像处理的一例的流程图。处理部10首先从扫描变换器8中取得测量信号,然后对其实施通常的信号处理而作成B模式动态图像(步骤21、22)。接着,从B模式动态图像中提取期望帧和与该期望帧的时间不同的帧这两帧(步骤23)。例如,提取期望帧和下一帧这两帧。然后,根据这两个帧算出运动向量场(步骤24)。运动向量场的算出是基于块匹配方法实施的。然后,对算出的运动向量场进行噪声去除处理(步骤25),将进行了噪声去除的运动向量场变换成标量场(步骤26)。然后,对标量场图像和运动向量场图像或B模式图像进行合成显示,结束1图像份的处理(步骤27)。在步骤23中,作为期望帧,在时间序列上依次选择不同的帧,反复进行上述步骤21~27的处理,并连续显示合成图像,由此可以显示合成图像的运动图像。
图3是表示步骤24的详细处理的流程图,图4是说明步骤24的块匹配处理的图。利用图3及图4,具体说明上述步骤24的运动向量场算出用的块匹配处理。在此,设在步骤23中选择了第m帧和第m+Δ帧。例如,设Δ=1帧。首先,处理部10在帧m中如图4所示那样设定规定的像素数N的关心区域ROI(region of interest:基准块)31(步骤51)。将ROI31中包含的像素的亮度分布表示为Pm(i0,j0)。i0、j0表示ROI31内的该像素的位置。然后,处理部10在第m+Δ帧中,在帧m的ROI31所对应的位置及其附近处设定规定的大小的探索区域32(步骤52)。在这里,关于ROI31的设定,说明了处理部10在帧m的整个图像中依次设定ROI31,并设定以其作为中心的规定的大小的探索区域32的构成,但是处理部10也可以设定预先规定的位置及大小的ROI31和在其附近预先规定的大小的探索区域32,或者处理部10也可以在参数设定部11从操作员接受的关心区域(ROI)及探索区域内设定ROI31及探索区域32。
探索区域32被划分成与ROI31相等大小的多个移动候补区域33。处理部10通过计算而求出相对于ROI31的亮度而言类似度最高的移动候补区域33,选择为移动目标区域。作为表示类似度的指标,可采用差分绝对值和、均方根、互相关值等。在这里,以下作为一例而说明利用差分绝对值和的情况。
将探索区域32内的移动候补区域33中包含的像素的亮度分布表示为Pm+Δ(i,j)。i、j表示移动候补区域33内的该像素的位置。处理部10计算ROI31的像素的亮度分布Pm(i0,j0)和移动候补区域33的亮度分布Pm+Δ(i,j)之间的差分绝对值和SAD(sum of absolute difference)(步骤 53)。在这里,SAD由下式(1)进行定义。
[公式1]
处理部10针对探索区域32的所有移动候补区域33求出与ROI31之间的SAD值,将在求得的SAD分布之中SAD值最小的移动候补区域33判定为移动目标区域,决定将ROI31的位置和最小SAD值的移动候补区域33的位置进行连结的运动向量(步骤54)。
之后,处理部10使ROI31在帧m的整个图像中移动的同时,反复进行上述处理,由此针对帧m的图像整体而决定运动向量(步骤55)。通过生成例如用箭头表示所决定的向量的图像,由此得到运动向量场(运动向量分布图像)。
图5(a)、(b)是通过上述处理得到的B模式图像例和运动向量分布图像例。图5(a)的B模式图像是将凝胶基材仿真图41、42叠加成2层,将超声波探测器固定于上侧的仿真图41,在沿着横向移动的同时进行拍摄而得到的图像。将图5(a)的B模式图像作为帧m,通过与下一帧(帧m+Δ,Δ=1帧)之间的块匹配处理(步骤24)而求得的运动向量为图5(b)。
如图5(b)所示,固定了超声波探测器1的上侧仿真图41所对应的上侧1/3左右的区域相对处于静止状态,下侧的仿真图42成为表示横向移动的向量场(朝向侧面的箭头)。然而,关于图5(b)的下侧仿真图42的下部1/3左右的区域,箭头指向斜上方,朝向不固定,发现运动向量紊乱的现象。该现象是由于距探测器1的距离越远检测灵敏度的SN比(SNR)越下降而引起的,表示了穿透界限。即,在远离探测器1的低SNR区域产生了错误向量。
利用图6,说明上述步骤25的运动向量场的噪声去除处理。
图6(a)、(b)是将在图5(b)的位置(3)及(5)分别设定ROI31并通过步骤24求得的下一帧的探索区域32的各移动候补区域33的SAD值,分别表示为移动候补区域33的浓度的图(SAD分布图)。图6(a)、(b)的整体分别表示探索区域32,探索区域32被划分成21×21个移动 候补区域33。设移动候补区域33的块尺寸为30×30像素,设探索区域32为50×50像素,并使移动候补区域33在探索区域32内按每次1像素地移动。探索区域32被设定成:ROI31的位置所对应的移动候补区域33的位置位于探索区域32的中心。
在图6(a)的SAD分布中,从探索区域32的中心偏向右侧旁边的位置的移动候补区域33的SAD值最小。由此可知,针对图5(b)的位置(3),在步骤24中决定了朝右侧面的向量。由图5(b)可确认,位置(3)是2层仿真图41、42的边界的偏下侧(仿真图42侧)附近,显示了朝右侧面的向量。
在这里,关注图6(a)的SAD值的空间分布可知,在SAD值为最小的移动候补区域33的周围,在横向即沿着2层仿真图41、42的边界的方向上,形成了SAD值小的区域(SAD值的谷)。如图6(b)所示,由该现象启示出:即便不使ROI31在帧m的整个区域移动并作成所有的运动向量场,也能够根据图6(a)的SAD分布直接检测出边界。
另一方面,因为图6(b)的SAD值分布是图5(b)的运动向量紊乱的区域内的位置(5)处的SAD值分布,所以可知SAD值最小的移动候补区域33在噪声值变小的靠近探测器1的上部的较宽区域内均匀展宽,因此运动向量朝上,变得容易紊乱。另外,可知在SAD值最小的移动候补区域33的周围应该形成的SAD值小的区域(SAD值的谷)也被整体的噪声变动淹没,因此无法识别出。
在本发明中,如图6(a)所示,在噪声小的位置(3)处设定了ROI31的SAD分布图像中,由于SAD值最小的移动候补区域33采取明显小于周边区域的SAD值,所以能够清晰地识别出,相反,如图6(b)所示,在噪声多的位置(5)处设定了ROI31的探索区域32的SAD分布图像中,无法清晰地识别出SAD值最小的移动候补区域33,利用上述现象来判别噪声多的ROI31。对噪声多的ROI31进行去除所对应的运动向量的处理。
处理部10为了进行上述判别,例如在图2的步骤25中,根据在规定的位置处设定了ROI31的情况下的图6(a)、(b)那样的探索区域32的SAD分布,如图6(c)及(d)所示那样作成表示SAD值的分布的直方图。在图5(b)的位置(3)处设定了ROI31的情况下的直方图如图6(c) 所示,SAD最小值与直方图分布中的频度大的SAD值充分地分离。即,最小的SAD值和频度大的SAD值的范围被充分地分离。另一方面,在图5(b)的位置(5)处设定了ROI31的情况下的直方图如图6(d)所示,因为相对于SAD值的分布而言频度的差并不大,因此直方图分布以宽幅扩展。因此,SAD最小值包含在频度大的SAD值的范围内,最小的SAD值和频度大的SAD值的范围之间的分离不充分。
根据这种特性,通过比较探索区域32的SAD分布的直方图的SAD最小值和频度大的SAD值,能够判别在该探索区域32中对应的ROI31的信号的可靠度(噪声的多少)、及在该探索区域32中决定的运动向量的可靠度。由此,能够将可靠度低的区域判别为低SNR区域,并且也能够判别所对应的运动向量的可靠度。
在本发明中,在SAD值分布的直方图中,为了判断SAD最小值与频度大的SAD值之间的分离度而采用下述指标。首先,说明以分离度参数作为指标的情况下的处理方法。图7示出分离度的定义的概念。分离度是相当于直方图的分布平均值与最小值之间的距离的值,如下式(2)进行定义。
[公式2]
SAD值平均值
Smin:最小SAD值
σs:SAD值标准偏差
在式(2)中,为了避免因分布的差异而带来的影响,用标准偏差实施了规格化。
图8示出判别以分离度作为指标的情况下的噪声较大的区域,并去除运动向量的处理的流程。该处理具体表示图2的步骤25的处理,针对图3的步骤51中在帧m中设定的所有ROI31进行。处理部10首先决定作为对象的ROI31(步骤81),针对所决定的ROI31,利用在图3的步骤52、53中设定及运算出的探索区域32的所有移动候补区域33的SAD值,根 据统计处理的方式计算SAD值的平均值、最小值、标准偏差(步骤82)。然后,求出用上式(2)定义的分离度(步骤83)。针对所有ROI31反复进行该处理(步骤84)。关于通过步骤84求得的分离度低于规定值的ROI31,因为是噪声大、且运动向量的可靠性低的区域,所以将通过图3的步骤54求得的运动向量置换为0(步骤85)。由此,能够根据运动向量图像来判别可靠性低的区域,并去除向量(错误向量)。
在上述步骤85中,作为用于判定分离度是否低的规定值,能够利用预先规定的阈值、或者在步骤84中针对所有ROI31求得的分离度的分布的中位数。或者,生成求得的分离度的直方图,在形成有多个频度的山的情况下,作为阈值,可使用位于分离度最低的一侧的山与位于分离度比其高的一侧的山之间的谷的位置的分离度。
使在步骤83中针对所有ROI31求得的分离度图像化之后的图像在图9中示出。在图9中,在帧m中设定33×51个ROI31,各ROI31的分离度通过浓度进行表示。如图9可知,分离度在帧m的下部的低SNR区域变低,分离度反映出运动向量推定的可靠性。
在上述图8的处理中,虽然利用了分离度,但是作为用于判断SAD最小值与频度大的SAD值之间的分离度的指标,也可使用其他指标。例如,可以使用变动系数。变动系数如下式进行定义,是用平均值对标准偏差进行规格化而得到的统计量,表示分布的偏差的大小(即,最小值的分离难易度)。
[公式3]
图10示出将以变动系数作为指标的情况下的噪声较大的区域的向量去除的处理的流程。该处理与图8的处理流程同样地,决定作为对象的ROI31(步骤81),针对所决定的ROI31,利用在图3的步骤52、53中设定及运算出的探索区域32的所有移动候补区域33的SAD值,根据统计处理的方式计算SAD值的平均值和标准偏差(步骤101)。然后,求出由上式(3)定义的变动系数(步骤102)。针对所有ROI31反复进行上述处理(步骤84)。关于通过步骤102求得的变动系数大于规定值的ROI31,将通过图3的步骤54求得的运动向量置换为0(步骤85)。由此,能够判别噪声大的ROI31,从而能够根据运动向量图像判别可靠性低的区域,并去除向量(错误向量)。
在上述步骤85中,作为用于判定变动系数是否高的规定值,能够利用预先规定的阈值、或者在步骤84中针对所有ROI31求得的变动系数的分布的中位数。或者,生成求得的变动系数的直方图,在频度呈两个山形的情况下,作为规定值而采用两个山之间的谷的最小值也是有效的。
图11示出使在步骤102中针对所有ROI31求得的变动系数图像化之后的图像。在图11中,各ROI31的变动系数通过浓度进行表示。由图11可知,变动系数在帧m的下部的低SNR区域变大,变动系数反映出运动向量推定的可靠性。
图12(a)、(b)示出错误向量去除前和去除后的运动向量分布图像的例子。图12(a)是与图5(b)相同的错误向量去除前的运动向量场,图12(b)是通过图10的处理,根据SAD分布求出变动系数分布,将变动系数分布的中位数作为阈值,将变动系数比该阈值大的ROI判定为可靠性低,并将运动向量设为0(静止状态)的图。将图12(a)、(b)对比可知,关于运动向量紊乱的下侧的区域,其错误向量被明显地去除了,并置换成静止状态。即,能够将下侧的区域判定为未准确得到超声波回波的穿透区域(即可靠性低的区域)。
在通过上述的图8及图10的处理去除了运动向量分布的错误向量之后,通过图2的步骤26、27将运动向量分布变换成标量分布,对标量分布图像和运动向量分布图像(或B模式图像)进行合成显示。
此外,在图8及图10的处理中,去除了可靠性低的区域的运动向量并使之处于静止状态,但是本发明并不限定于该处理方法。例如,也可取代将运动向量设定为静止状态的方式,而采用仍保持之前针对相同区域得到的运动向量的状态的处理方法。
(实施方式2)
在实施方式1中,在求出SAD分布之后,根据SAD分布判别ROI31是否为低SNR区域,并在为低SNR区域的情况下去除运动向量等,而在本实施方式中,在运算SAD分布时去除噪声,并利用噪声去除后的检测灵敏度高的SAD分布来求出可靠性高的运动向量。
图13是在运算本实施方式2的SAD分布时去除噪声的处理流程。图13的处理,在实施方式1的图3的步骤51~54的SAD运算处理中追加了噪声去除处理(步骤132、133)。图14(a)、(b)、(c)表示图13的各处理阶段的SAD分布的例子。
如图13所示,处理部10首先进行与实施方式1的图3的步骤51~53相同的处理,求出SAD值分布。将求得的SAD分布图像的例子表示在图14(a)中。ROI31是图5(b)的位置(4),虽然仿真图42原本沿着横向相对移动,但因噪声的影响使得上侧的移动候补区域33的SAD值最小,所以若直接这样决定运动向量,则会产生向量的误检测。为了避免这种误检测,处理部10对通过步骤51~53得到的SAD分布图像实施平滑化处理(低通滤波器LPF(low pass filter)处理)(步骤132)。平滑化处理,例如将规定大小的滤波器施加于SAD分布图像,削减滤波器内的SAD分布的高频分量而进行平滑化。在使滤波器按规定量移动的同时反复进行该处理。这样,因为仿真图42的移动而带来的SAD值的变化急剧,因此通过SAD值分布图像的平滑化能够去除,相对于此,因为噪声导致的AD值的变动缓慢,因此通过SAD值分布图像的平滑化能够进行提取。图14(b)示出通过平滑化处理而得到的SAD值分布图像。
接着,求出步骤53的原始SAD值分布和步骤132的平滑化处理后的SAD值分布之间的差分(步骤133)。由此,能够得到噪声导致的SAD值的变动被去除后的、因仿真图的移动而带来的原本的SAD值分布。图14(c)示出所得到的分布。利用所得到的SAD值分布来进行图3的步骤54,将最小的SAD值的移动候补区域33判定为移动目标,决定运动向量。运动向量决定后,通过实施方式1的图3的步骤56的处理来生成运动向量分布图像。另外,通过实施方式1的图2的步骤25还可以进行去除运动向量分布的可靠性低的向量等的处理。
如上述,通过本实施方式2的处理,能够利用因噪声导致的SAD值变动被去除后的SAD值分布来决定运动向量,因而能够提高运动向量的可靠性。
此外,在上述平滑化处理中,虽然利用了LPF,但是并不限定于此, 在因被检者(仿真图)的移动而带来的SAD值的分布的空间频率高(即更为复杂的形状)的情况下,应用带通滤波器也是有效的。
在步骤132的滤波处理中用到的滤波器的一边的大小能够如下决定。即,对未进行平滑化处理的SAD分布预先进行图3的步骤54而作成运动向量场,求出该运动向量场的最大向量长度,将最大向量长度设定为低通滤波器或者带通滤波器的一边的大小。
(实施方式3)
下面,作为实施方式3而说明下述的处理方法:利用通过实施方式1的图3的步骤53求得的SAD值分布,直接求出组织的边界,并且判定活体肿瘤相对于正常组织的浸润度。其中,在本实施方式3中,将探索区域32的移动候补区域33简称为区域33。另外,将ROI31也称作关注像素。
由于沿着被检者的组织的边界的区域其组织的类似度高,所以在B模式图像中表现出相近的亮度。为此,SAD的特征在于沿着被检者的组织的边界的区域33表现出比沿着边界的正交方向的区域33小的值。另一方面,若活体肿瘤的浸润扩展,则边界变得不清楚,所以沿着边界的区域33的SAD值变大。利用该特征来判定浸润度。
图15表示本实施方式3的处理部10的处理流程。另外,图16(a)~(h)表示对象方向和与之对应地在SAD值分布图像上选择的区域33的8个图案。
首先,处理部10在B模式图像的期望帧m中想要调查的组织的边界位置设定关注像素(ROI)31,在帧m+Δ中设定探索区域32,求出探索区域32的SAD值分布(步骤151)。帧的选择方法及SAD值分布的运算方法与实施方式1的图2的步骤21~23及图3的步骤51~53同样地进行。
在该SAD值分布中,如图16(a)所示,选择通过探索区域32的中心且沿着规定的对象方向(水平方向)151配置的区域33(步骤63),求出所选择的区域33的SAD值之和(步骤64)。同样地,选择沿着对象方向151的正交方向(垂直方向)152配置的区域33,也求出所选择的区域33的SAD值之和。
反复执行该步骤63、64的处理,直至对图16(a)~(h)的8个图案的所有图案执行完毕为止(步骤62)。在图16(b)的图案中,求出沿 着规定的对象方向(相对于水平方向以逆时针旋转的方式倾斜了约30°的方向)151配置的区域33的SAD值之和,求出沿着其正交方向152配置的区域33的SAD值之和。
在图16(c)~(h)的图案中,分别求出沿着相对于水平方向以逆时针旋转的方式倾斜约45°、约60°、90°、约120°、约135°、约150°的对象方向151配置的区域33的SAD值之和、以及沿着其正交方向152配置的区域33的SAD值之和。
选择求得的各对象方向151的SAD值之和最小的对象方向151(步骤65)。所选择的对象方向151的方向为组织的边界的方向。由此,能够在不求出运动向量的情况下直接检测边界。
接着,选择正交于所选择的对象方向151的方向152(步骤66)。算出所选择的对象方向151的SAD值之和、和与其正交的方向152的SAD值之和的比率(对象方向的SAD值之和/正交方向的SAD值之和)(步骤67)。
在浸润度低而边界明确的情况下,边界方向(所选择的对象方向151)的SAD和变小、正交方向152的SAD和变大,所以作为比率会得到较小的值。另一方面,随着浸润度变大而边界变得不清楚,边界方向(所选择的对象方向151)的SAD和也会增大,所以比率也逐渐增大。因此,能够将比率作为参数来评价浸润度。具体而言,例如将预先规定的多个基准值与比率进行对比来判定浸润的程度,并显示判定结果。
此外,可以在比率小于预先设定的一定值的情况下,识别为该关注像素(ROI31)是构成边界线的点,并显示边界。
此外,在求出各方向的SAD值之和时,能够利用依赖方向的滤波器。依赖方向的滤波器是具有在处理像素的滤波范围(探索区域32)中判定一维方向的浓度变化最小的方向的功能的滤波器。
在图16(a)~(h)中,为了易于图示,示出了由5×5的区域33构成的探索区域32中对象方向151及正交方向152的区域选择图案,但在实际处理中,对应于探索区域32的区域33的数量而设定区域选择图案。
(实施方式4)
作为实施方式4,说明根据探索区域32的SAD值分布在不求出运动 向量的情况下直接检测边界的其他方法。在这里,应用进行与2次微分对应的强调处理的拉普拉斯滤波器。
在图17(a)及(b)中,在图5(b)的位置(1)和(2)处设定ROI31,并示出探索区域32的SAD分布。位置(1)是相对于探测器1而言相对静止的仿真图41的内部的位置。位置(2)位于仿真图41与相对于该仿真图41相对横向移动的仿真图42之间的边界附近。求出位置(1)、(2)处的SAD分布的处理,与实施方式1的图2的步骤21~23及图3的步骤51~53同样地进行。若对得到的SAD分布图像(图17(a)及(b))应用进行空间2次微分的拉普拉斯滤波器,则根据轮廓强调效果,可分别得到SAD值的变动大的部分被强调后的图17(c)及(d)的图像。
在探索区域32中存在边界的情况下,如图17(d)所示,生成了沿着边界的SAD值的变动大的区域被强调并提取为条状的轮廓强调分布。由此,通过2值化处理对轮廓强调分布中连续的轮廓线的区域(移动候补区域33)进行提取显示,由此能够根据SAD值的拉普拉斯图像检测边界。另一方面,在探索区域32中不存在边界的情况下,如图17(c)所示,SAD值的变动大的区域仅是中心区域(与ROI31的位置对应的移动候补区域33),与周围孤立,不是连续的轮廓线。由此可知边界不位于该位置。
根据上述处理,能够对SAD分布实施拉普拉斯处理而直接提取边界。由此,能够省略实施方式1的用于决定运动向量并使其图像化的图3的步骤54、55、以及进行运动向量的噪声去除或向标量分布的变换来推定边界的图2的步骤25~26,因此能大幅削减处理量。
(产业上的可利用性)
本发明能够应用于医用超声波诊断和治疗装置、以及利用包含超声波的电磁波等来测量失真或偏差的所有装置。
符号说明
1-超声波探头(探测器);2-用户接口;3-送波波束成形器;4-控制系统;5-送受切换开关;6-受波波束成形器;7-包络线检波部;8-扫描变换器;10-处理部;10a-CPU;10b-存储器;11-参数设定部;12-合成部;13-显示部。
Claims (12)
1.一种超声波成像装置,其特征在于,具有:
发送部,其向对象发送超声波;
接收部,其接收来自所述对象的超声波;和
处理部,其处理所述接收部的接收信号,生成2帧以上的图像,
所述处理部将所述生成的2帧以上的图像中的1帧作为基准帧,在预先确定的位置或从操作者接受的位置设定关心区域,将另外的1帧作为比较帧,在预先确定的位置或从操作者接受的位置设定比所述关心区域宽的探索区域,在所述探索区域内设定作为所述关心区域的移动目标的候补的多个候补区域,按每个所述候补区域计算所述关心区域内和所述候补区域内的图像特征值的类似度,求出整个所述探索区域内的所述类似度的分布,
所述处理部在所述类似度的分布中求出将所述类似度的最小值与所述类似度的整体值进行比较的统计量,根据该统计量判定所述关心区域的可靠度。
2.根据权利要求1所述的超声波成像装置,其特征在于,
所述处理部利用所述类似度的最小值、平均值及标准偏差求出所述统计量,通过比较求得的所述统计量和阈值来判定所述关心区域的可靠度。
3.根据权利要求2所述的超声波成像装置,其特征在于,
所述处理部生成将所述比较帧中与所述关心区域对应的位置和所述类似度为最小的所述候补区域的位置进行连结的向量,对于被判定为所述可靠度低的所述关心区域,将所述向量置换为零或规定的向量。
4.根据权利要求1所述的超声波成像装置,其特征在于,
所述处理部针对所述类似度分布而算出类似度的平均值、最小值及标准偏差,并求出将所述平均值与所述最小值之间的差分除以所述标准偏差而得到的分离度来作为所述统计量。
5.根据权利要求1所述的超声波成像装置,其特征在于,
所述处理部针对所述类似度分布而算出类似度的平均值及标准偏差,并求出将所述标准偏差除以所述平均值而得到的变动系数来作为所述统计量。
6.根据权利要求2所述的超声波成像装置,其特征在于,
所述处理部设定多个所述关心区域,对各关心区域求出所述统计量,并求出表示求得的统计量的值的频度的直方图分布,
作为所述阈值,使用所述直方图分布中的中位数、平均值、或者直方图分布呈多个山形状时山之间的谷的最小值的统计量。
7.根据权利要求1所述的超声波成像装置,其特征在于,
所述处理部生成将所述类似度分布进行了平滑化处理的平滑化后类似度分布,并求出从平滑化处理前的所述类似度分布中减去所述平滑化后类似度分布而得到的差分类似度分布。
8.根据权利要求7所述的超声波成像装置,其特征在于,
所述平滑化处理是下述处理,即:对于所述类似度分布,设定预先确定的大小的滤波器,并在使滤波器移动规定量的同时反复进行将该滤波器内的分布平滑化的处理,
所述处理部针对多个所述关心区域,生成将所述比较帧中与所述关心区域对应的位置和在所述平滑化处理前的类似度分布中所述类似度为最小的所述候补区域的位置进行连结的向量,并将生成的向量之中最大的向量长度作为所述滤波器的大小。
9.根据权利要求1所述的超声波成像装置,其特征在于,
所述处理部用拉普拉斯滤波器对所述类似度分布进行滤波处理来作成轮廓强调分布,通过提取所述轮廓强调分布中连续的轮廓线来求出所述对象的边界。
10.根据权利要求1所述的超声波成像装置,其特征在于,
所述处理部具有:
第1处理单元,其针对在活体肿瘤边界附近设定的所述关心区域而求出所述类似度分布,生成以该类似度作为图像特征值的类似度分布图像,在该类似度分布图像上以所述关心区域所对应的位置为中心在多个不同方向上设定规定长度的一维区域;
第2处理单元,其按照设定的每个方向计算所述一维区域内的所述类似度的总和;
第3处理单元,其计算所述类似度总和最小的方向上的类似度总和与正交于该方向的方向上的一维区域的类似度总和之间的比率;和
第4处理单元,其基于所述比率来判定肿瘤的浸润度。
11.根据权利要求10所述的超声波成像装置,其特征在于,
在所述第3处理单元计算出的所述比率小于预先设定的一定值的情况下,所述处理部判断该关心区域是构成边界线的点。
12.一种超声波成像方法,其特征在于,
向对象发送超声波,对接收来自所述对象的超声波而得到的接收信号进行处理,生成2帧以上的图像,
从所述图像中选择基准帧和比较帧,
在所述基准帧中设定关心区域,
在所述比较帧中设定比所述关心区域宽的探索区域,在所述探索区域内设定作为所述关心区域的移动目标的候补的多个候补区域,
按每个所述候补区域计算所述关心区域内和所述候补区域内的图像特征值的类似度,求出整个所述探索区域内的所述类似度的分布,
在所述类似度的分布中求出将所述类似度的最小值与所述类似度的整体值进行比较的统计量,根据该统计量判定所述关心区域的可靠度。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2009-246734 | 2009-10-27 | ||
JP2009246734 | 2009-10-27 | ||
PCT/JP2010/068988 WO2011052602A1 (ja) | 2009-10-27 | 2010-10-26 | 超音波イメージング装置、超音波イメージング方法および超音波イメージング用プログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102596050A CN102596050A (zh) | 2012-07-18 |
CN102596050B true CN102596050B (zh) | 2014-08-13 |
Family
ID=43922027
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201080046798.3A Expired - Fee Related CN102596050B (zh) | 2009-10-27 | 2010-10-26 | 超声波成像装置及超声波成像方法 |
Country Status (5)
Country | Link |
---|---|
US (1) | US8867813B2 (zh) |
EP (1) | EP2494924A1 (zh) |
JP (1) | JP5587332B2 (zh) |
CN (1) | CN102596050B (zh) |
WO (1) | WO2011052602A1 (zh) |
Families Citing this family (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101819028B1 (ko) * | 2011-07-11 | 2018-01-16 | 삼성전자주식회사 | 초음파 영상 처리 방법 및 장치 |
US20160213353A1 (en) * | 2011-10-28 | 2016-07-28 | Hironari Masui | Ultrasound imaging apparatus, ultrasound imaging method and ultrasound imaging program |
KR101881924B1 (ko) * | 2011-12-29 | 2018-08-27 | 삼성전자주식회사 | 초음파 영상 처리 방법 및 장치 |
US8897532B2 (en) * | 2012-07-11 | 2014-11-25 | General Electric Company | Systems and methods for performing image type recognition |
AU2013248207A1 (en) * | 2012-11-15 | 2014-05-29 | Thomson Licensing | Method for superpixel life cycle management |
WO2014080833A1 (ja) | 2012-11-21 | 2014-05-30 | 株式会社東芝 | 超音波診断装置、画像処理装置及び画像処理方法 |
WO2014103501A1 (ja) * | 2012-12-28 | 2014-07-03 | 興和株式会社 | 画像処理装置、画像処理方法、画像処理プログラム及びそのプログラムを格納した記録媒体 |
CN103093418B (zh) * | 2013-02-21 | 2015-08-26 | 深圳市晶日盛科技有限公司 | 一种改进的数字图像缩放方法 |
WO2014155272A1 (en) * | 2013-03-28 | 2014-10-02 | Koninklijke Philips N.V. | Real-time quality control for acquisition of 3d ultrasound images |
JP5651258B1 (ja) * | 2014-02-27 | 2015-01-07 | 日立アロカメディカル株式会社 | 超音波診断装置及びプログラム |
US9686470B2 (en) * | 2014-05-30 | 2017-06-20 | Apple Inc. | Scene stability detection |
US11337673B2 (en) * | 2016-06-10 | 2022-05-24 | Koninklijke Philips N.V. | Using reflected shear waves for monitoring lesion growth in thermal ablations |
JP6587583B2 (ja) * | 2016-06-27 | 2019-10-09 | 株式会社日立製作所 | 超音波撮像装置、超音波撮像方法、および、結合状態評価装置 |
CN108225496B (zh) * | 2016-12-15 | 2020-02-21 | 重庆川仪自动化股份有限公司 | 雷达物位计回波信号自动测试装置、方法及系统 |
US11611773B2 (en) * | 2018-04-06 | 2023-03-21 | Qatar Foundation For Education, Science And Community Development | System of video steganalysis and a method for the detection of covert communications |
US11527053B2 (en) * | 2018-12-12 | 2022-12-13 | Samsung Electronics Co., Ltd. | Method and apparatus of processing image |
CN109875607A (zh) * | 2019-01-29 | 2019-06-14 | 中国科学院苏州生物医学工程技术研究所 | 浸润组织检测方法、装置及系统 |
JP6579727B1 (ja) * | 2019-02-04 | 2019-09-25 | 株式会社Qoncept | 動体検出装置、動体検出方法、動体検出プログラム |
CN113556979A (zh) * | 2019-03-19 | 2021-10-26 | 奥林巴斯株式会社 | 超声波观测装置、超声波观测装置的工作方法以及超声波观测装置的工作程序 |
CN116823829B (zh) * | 2023-08-29 | 2024-01-09 | 深圳微创心算子医疗科技有限公司 | 医疗影像的分析方法、装置、计算机设备和存储介质 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101378700A (zh) * | 2006-03-31 | 2009-03-04 | 株式会社东芝 | 超声波诊断装置以及超声波诊断装置控制方法 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4318862B2 (ja) * | 2001-01-25 | 2009-08-26 | 株式会社日立メディコ | 画像診断装置 |
JP4565796B2 (ja) | 2002-07-25 | 2010-10-20 | 株式会社日立メディコ | 画像診断装置 |
JP4596759B2 (ja) | 2002-09-12 | 2010-12-15 | 株式会社日立メディコ | 画像診断装置、及び画像診断装置の作動方法 |
JP4130114B2 (ja) * | 2002-10-09 | 2008-08-06 | 株式会社日立メディコ | 超音波イメージング装置及び超音波信号処理方法 |
JP3932482B2 (ja) | 2002-10-18 | 2007-06-20 | 株式会社日立メディコ | 超音波診断装置 |
JP4321121B2 (ja) | 2003-05-30 | 2009-08-26 | 株式会社日立メディコ | 診断画像における生体組織の動きの追跡方法及びその方法を用いた画像診断装置 |
JP4751282B2 (ja) | 2006-09-27 | 2011-08-17 | 株式会社日立製作所 | 超音波診断装置 |
CN102202580B (zh) | 2008-11-10 | 2013-11-20 | 株式会社日立医疗器械 | 超声波图像处理方法及装置、超声波图像处理程序 |
-
2010
- 2010-10-26 WO PCT/JP2010/068988 patent/WO2011052602A1/ja active Application Filing
- 2010-10-26 US US13/503,858 patent/US8867813B2/en active Active
- 2010-10-26 JP JP2011538440A patent/JP5587332B2/ja not_active Expired - Fee Related
- 2010-10-26 CN CN201080046798.3A patent/CN102596050B/zh not_active Expired - Fee Related
- 2010-10-26 EP EP10826734A patent/EP2494924A1/en not_active Withdrawn
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101378700A (zh) * | 2006-03-31 | 2009-03-04 | 株式会社东芝 | 超声波诊断装置以及超声波诊断装置控制方法 |
Non-Patent Citations (5)
Title |
---|
JP特开2002-222410A 2002.08.09 |
JP特开2004-121834A 2004.04.22 |
JP特开2004-129773A 2004.04.30 |
JP特开2004-351050A 2004.12.16 |
JP特开2004-57275A 2004.02.26 |
Also Published As
Publication number | Publication date |
---|---|
US8867813B2 (en) | 2014-10-21 |
US20120224759A1 (en) | 2012-09-06 |
WO2011052602A1 (ja) | 2011-05-05 |
CN102596050A (zh) | 2012-07-18 |
JP5587332B2 (ja) | 2014-09-10 |
JPWO2011052602A1 (ja) | 2013-03-21 |
EP2494924A1 (en) | 2012-09-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102596050B (zh) | 超声波成像装置及超声波成像方法 | |
EP1117331B1 (en) | Adaptive cross-sectional area computation using statistical signatures | |
US6312385B1 (en) | Method and apparatus for automatic detection and sizing of cystic objects | |
US9585636B2 (en) | Ultrasonic diagnostic apparatus, medical image processing apparatus, and medical image processing method | |
US9445744B2 (en) | Methods, systems, and devices for spine centrum extraction and intervertebral disk dividing | |
KR101121353B1 (ko) | 2차원 초음파 영상에 대응하는 2차원 ct 영상을 제공하는 시스템 및 방법 | |
KR101121396B1 (ko) | 2차원 초음파 영상에 대응하는 2차원 ct 영상을 제공하는 시스템 및 방법 | |
US20030174890A1 (en) | Image processing device and ultrasonic diagnostic device | |
CN103906473B (zh) | 超声波成像装置、超声波成像方法 | |
KR20150120774A (ko) | 관심영역 검출 시스템 및 방법 | |
CN110786880B (zh) | 超声波诊断装置以及超声波图像处理方法 | |
US20140018682A1 (en) | Ultrasonic diagnostic apparatus and ultrasonic diagnostic image rendering method | |
CN106204570B (zh) | 一种基于非因果分数阶梯度算子的角点检测方法 | |
US10188370B2 (en) | Ultrasound imaging system and method | |
US11944499B2 (en) | Tracking contrast elements in ultrasound imaging | |
JP5756812B2 (ja) | 超音波動画像処理方法、装置、およびプログラム | |
JP3662835B2 (ja) | 超音波診断装置 | |
US20080137928A1 (en) | Method and System for Registering CT Data Sets | |
Liu et al. | Automatic segmentation of prostate boundaries in transrectal ultrasound (TRUS) imaging | |
JP2005160916A (ja) | 石灰化陰影判定方法、石灰化陰影判定装置及びプログラム | |
CN116823861A (zh) | 一种机器人识别天花板直线的方法 | |
CN113040822A (zh) | 子宫内膜蠕动的测量方法、用于测量子宫内膜蠕动的设备 | |
CN116058875A (zh) | 超声成像方法和超声成像系统 | |
KR20170001547A (ko) | 움직임 추적을 위한 관심영역 설정 장치 및 방법과 그 기록매체 |
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 | ||
TR01 | Transfer of patent right |
Effective date of registration: 20170327 Address after: Tokyo, Japan, Japan Patentee after: Hitachi Ltd. Address before: Japan Tokyo Chiyoda District Kanda four chome 14 No. 1 Patentee before: Hitachi Medical Corporation |
|
TR01 | Transfer of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140813 Termination date: 20171026 |
|
CF01 | Termination of patent right due to non-payment of annual fee |