CN113768533A - 超声显影装置和超声显影方法 - Google Patents
超声显影装置和超声显影方法 Download PDFInfo
- Publication number
- CN113768533A CN113768533A CN202010524116.0A CN202010524116A CN113768533A CN 113768533 A CN113768533 A CN 113768533A CN 202010524116 A CN202010524116 A CN 202010524116A CN 113768533 A CN113768533 A CN 113768533A
- Authority
- CN
- China
- Prior art keywords
- ultrasonic
- beams
- module
- speckle
- imaging
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 19
- 238000003384 imaging method Methods 0.000 claims abstract description 74
- 230000008859 change Effects 0.000 claims abstract description 52
- 238000002592 echocardiography Methods 0.000 claims abstract description 36
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 25
- 238000003786 synthesis reaction Methods 0.000 claims abstract description 25
- 238000012545 processing Methods 0.000 claims abstract description 16
- 238000002604 ultrasonography Methods 0.000 claims description 30
- 230000006870 function Effects 0.000 claims description 27
- 238000001914 filtration Methods 0.000 claims description 22
- 230000005540 biological transmission Effects 0.000 claims description 10
- 238000012285 ultrasound imaging Methods 0.000 claims description 8
- 230000002194 synthesizing effect Effects 0.000 claims description 6
- 108091026890 Coding region Proteins 0.000 claims description 5
- 230000009131 signaling function Effects 0.000 claims description 4
- 230000003068 static effect Effects 0.000 claims description 4
- 230000017531 blood circulation Effects 0.000 abstract description 13
- 210000004204 blood vessel Anatomy 0.000 abstract description 9
- 238000006243 chemical reaction Methods 0.000 abstract description 2
- 201000010099 disease Diseases 0.000 abstract description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 abstract description 2
- 210000005075 mammary gland Anatomy 0.000 abstract 1
- 230000001427 coherent effect Effects 0.000 description 12
- 238000010586 diagram Methods 0.000 description 10
- 239000011159 matrix material Substances 0.000 description 7
- 239000002131 composite material Substances 0.000 description 6
- 230000005284 excitation Effects 0.000 description 6
- 239000000523 sample Substances 0.000 description 5
- 239000013598 vector Substances 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 239000012530 fluid Substances 0.000 description 4
- 230000004044 response Effects 0.000 description 4
- 238000013329 compounding Methods 0.000 description 3
- 239000012216 imaging agent Substances 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 206010020751 Hypersensitivity Diseases 0.000 description 1
- 230000036770 blood supply Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 210000001685 thyroid gland Anatomy 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/44—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
- A61B8/4483—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/06—Measuring blood flow
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/44—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
- A61B8/4483—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer
- A61B8/4494—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer characterised by the arrangement of the transducer elements
-
- 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/5207—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
-
- 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/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
Abstract
本发明涉及超声成像技术领域，涉及超声显影装置和超声显影方法，可以用于乳腺等疾病的诊断中，超声显影装置包括：编码发射模块用于激励换能器晶片阵列向待测区域连续发射N组超声平面波束，N为大于1的正整数；接收解码模块用于接收每组超声平面波束的回波；波束合成模块用于对接收到的各组超声平面波束的回波进行波束合成，对应获得连续的N帧显影图像；成像处理模块用于依次对N帧B模式图像进行斑点追踪，确定斑点在N帧B模式图像中的斑点变化信息；图像合成模块，用于将斑点变化信息的成像数据进行数据转换，获得包含斑点变化信息的B模式超声图像序列。超声显影装置和超声显影方法提高超声图像的分辨率，能够清楚显示微小血管血流状态。
Description
技术领域
本发明涉及超声成像技术领域，具体涉及超声显影装置和超声显影方法。
背景技术
在医用超声成像系统中，利用超声波在扫描区域的组织界面处发生反射和散射，通过接收和处理包含有组织特征信息的回波信号，获得扫描区域中的超声图像。
由于血流成像存在分辨力和灵敏度的限制，以及人体血管组织结构复杂，尤其对于血流量小、血流速低的微小血管，存在显像困难的问题；然而，对微小血管进行稳定可靠的血流显影，能够辅助医生达到对疾病进程进行鉴别判断的目的，例如根据颈动脉斑块的内部供血情况判断斑块稳定程度。
目前相关技术采用的超声显影技术，是通过显影剂来使得微小血管血流状态清楚显示。然而，由于显影剂直接作用于血循环中，部分人群对显影剂会产生不同程度过敏反应，风险难以把控。
发明内容
为了解决上述技术问题，本发明提供一种超声显影装置和超声显影方法，能够清楚显示微小血管血流状态。
作为本发明的第一方面，提供一种超声显影装置，包括：
编码发射模块，所述编码发射模块用于激励换能器晶片阵列向待测区域连续发射N组超声平面波束，N为大于1的正整数；
接收解码模块，所述接收解码模块用于接收每组所述超声平面波束的回波；
波束合成模块，所述波束合成模块用于对接收到的各组所述超声平面波束的回波进行波束合成，对应获得连续的N帧显影图像；
成像处理模块，所述成像处理模块包括B模式成像模块和斑点追踪模块；所述B模式成像模块，用于依次对所述N帧显影图像进行B模式成像，获得N帧B模式图像；所述斑点追踪模块用于依次对N帧B模式图像进行斑点追踪，确定斑点在N帧B模式图像中的斑点变化信息；
图像合成模块，所述图像合成模块用于将所述斑点变化信息与N帧显影图像进行合成，获得包含所述斑点变化信息的B模式超声图像序列。
可选的，所述编码发射模块用于通过脉冲编码序列，激励换能器晶片阵列向待测区域连续发射N组超声平面波束；
所述接收解码模块为接收解码模块，用于对接收到的每组所述超声平面波束的回波进行匹配解码。
可选的，所述编码发射模块，用于通过脉冲编码序列，激励换能器晶片阵列沿多个不同发射角度，向待测区域发射第n组超声平面波束中的多个超声平面波束，n为大于1且小于N的正整数。
可选的，所述波束合成模块，用于将在第n组超声平面波束中，多个超声平面波束的回波进行波束合成，获得第n帧显影图像，所述第n帧显影图像表达式为：
其中，x为第n帧显影图像像素点在与换能器晶片阵列平行方向上的位置；
z为第n帧显影图像像素点在待测区域深度方向上的位置；
x1为在所述第n帧显影图像像素点所对应介质点产生的回波，被换能器晶片阵列接收的位置；
α为所述换能器晶片阵列发射超声平面波束的发射角度；
m为大于1的正整数；
c为超声波在待测区域中的传播速度；
τ函数表示：换能器晶片阵列向第n帧显影图像像素点所对应介质点发出超声平面波束后，至接收到回波时的时间间隔函数；
RF函数表示：在所述第n帧显影图像像素点位置处产生的回波，被换能器晶片阵列接收形成的回波信号函数；
s函数表示：在特定发射角度下，所述待测区域的显影图像函数；
I函数表示：叠加在多个不同发射角度下所述待测区域的显影图像函数，获得的第n帧显影图像函数。
可选的，所述斑点追踪模块，用于计算相邻帧B模式图像中对应像素点成像数据的成像方差；根据所述成像方差确定静态像素点和动态像素点；追踪所述动态像素点成像数据的变化信息，形成斑点变化信息。
可选的，所述成像处理模块还包括滤波模块，所述滤波模块用于依次对所述N帧B模式图像进行滤波处理，抑制斑点噪声；
所述滤波模块将滤波处理后的N帧B模式图像输出给斑点追踪模块；
所述斑点追踪模块用于依次对抑制斑点噪声后的N帧B模式图像进行斑点追踪，确定在所述N帧显影图像中的斑点变化信息；
所述图像合成模块用于将显影图像斑点的斑点变化信息重迭显示在对应帧的所述B模式图像上。
可选的，所述滤波模块采用均值滤波器、中值滤波器、Lee滤波器、Sigma滤波器、Frost滤波器中的任意一种或多种对所述斑点变化信息进行滤波处理，抑制斑点噪声。
作为本发明的第二方面，提供一种超声显影方法，包括以下步骤：
激励换能器晶片阵列向待测区域连续发射N组超声平面波束，N为大于1的正整数；
接收每组所述超声平面波束的回波；
对接收到的各组所述超声平面波束的回波进行波束合成，对应获得连续的N帧显影图像；
依次对所述N帧显影图像进行B模式成像，获得N帧B模式图像；
依次对N帧B模式图像进行斑点追踪，确定斑点在N帧B模式图像中的斑点变化信息；
将所述斑点变化信息与N帧显影图像进行合成，获得包含所述斑点变化信息的B模式超声图像序列。
可选的，所述激励换能器晶片阵列向待测区域连续发射N组超声平面波束，包括：
通过脉冲编码序列，激励换能器晶片阵列向待测区域连续发射N组超声平面波束；
所述接收每组所述超声平面波束的回波步骤，还包括：对接收到的每组所述超声平面波束的回波进行匹配解码。
可选的，所述激励换能器晶片阵列向待测区域连续发射N组超声平面波束，包括：
通过脉冲编码序列，激励换能器晶片阵列沿多个不同发射角度，向待测区域发射第n组超声平面波束中的多个超声平面波束，n为大于1且小于N的正整数。
本发明的优点：本发明采用N组超声平面波束获得以较快速度获得N帧显影图像，缩短相邻帧的时间间隔，同时结合波束合成和斑点追踪从而提高超声图像的分辨率，能够清楚显示微小血管血流状态。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案，下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明第一方面实施例的结构框图。
图2为本发明第二方面实施例的流程图。
图3(左)为探头接收到传统聚焦波束回波示意图。
图3(右)为探头接收到超声平面波束回波示意图。
图4(左)为换能器晶片阵列收发传统聚焦波束示意图。
图4(右)为换能器晶片阵列收发超声平面波束回波示意图。
图5为发射角度为0度时，超声平面波束在待测区域中的传播示意图。
图6为发射角度为α度(α不为0)时，超声平面波束在待测区域中的传播示意图。
图7为传统单脉冲激励成像方式示意图。
图8为脉冲编码激励成像方式示意图。
图9为甲状腺超声图像中斑点的分部图。
图10为图9中A区域的放大图。
图11为图9的A区域中一个斑点中心的超声示意图。
具体实施方式
下面将结合附图，对本发明中的技术方案进行清楚、完整的描述，显然，所描述的实施例是本发明的一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在不做出创造性劳动的前提下所获得的所有其它实施例，都属于本发明保护的范围。
在本发明的描述中，需要说明的是，术语“中心”、“上”、“下”、“左”、“右”、“竖直”、“水平”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。此外，术语“第一”、“第二”、“第三”仅用于描述目的，而不能理解为指示或暗示相对重要性。
在本发明的描述中，需要说明的是，除非另有明确的规定和限定，术语“安装”、“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电气连接；可以是直接相连，也可以通过中间媒介间接相连，还可以是两个元件内部的连通，可以是无线连接，也可以是有线连接。对于本领域的普通技术人员而言，可以具体情况理解上述术语在本发明中的具体含义。
此外，下面所描述的本发明不同实施方式中所涉及的技术特征只要彼此之间未构成冲突就可以相互结合。
作为本发明的第一方面，提供一种超声显影装置。
实施例1：
参照图1，本实施例中的超声显影装置包括：
编码发射模块，所述编码发射模块用于激励换能器晶片阵列向待测区域连续发射N组超声平面波束，N为大于1的正整数。
接收模块，所述接收模块用于接收每组所述超声平面波束的回波。
波束合成模块，所述波束合成模块用于对接收到的各组所述超声平面波束的回波进行波束合成，对应获得连续的N帧显影图像，并通过伪像去除算法对去除所述N帧显影图像中的伪像；通过对各组所述超声平面波束的回波进行波束合成，能够提升所获得的各帧显影图像的对比度和信噪比。
成像处理模块，成像处理模块包括B模式成像模块和斑点追踪模块；所述B模式成像模块，用于依次对所述N帧显影图像进行B模式成像，获得N帧B模式图像；所述斑点追踪模块用于依次对所述N帧B模式图像进行斑点追踪，确定在所述N帧B模式成像中的斑点变化信息。
图像合成模块，所述图像合成模块用于将所述斑点变化信息与N帧B模式图像进行合成，获得包含所述斑点变化信息的B模式超声图像序列。
所述图像合成模块能够将合成后得到的包含所述斑点变化信息的B模式超声图像序列存储在图像存储器中。显示模块通过扫描图像存储器以读出其中的数据进行显示。
传统的超声成像是通过发射多次聚焦波束照射组织来成像(参照图3左和图4)，以单波束接收为例，聚焦波每获得一根接收线需要付出一次发射的代价，如果一帧图像有M根扫描线，则要形成一帧图像共需发射M次；而在平面波成像模式(参照图3右和图4右)下，同时激励换能器所有阵元，产生一个平行于换能器的超声平面波，此平面波产生的声场能够有效覆盖所有感兴趣区域，经组织散射后，所有阵元也同时接收回波信号，通过这些回波信号就可以得到整个区域的超声图像。因此平面波只需要进行一次发射即可得到一帧图像，和聚焦波相比，发射次数降为了1/M。即使聚焦波采用双波束、四波束等多波束并行处理技术，平面波在减少发射次数方面的优势仍然是非常可观的。
利用超声平面波束能够快速地对组织进行成像(优选每秒1000帧的速度)，通过在相邻帧之间的较短时间里，明显移动会显示为半波长量级的斑点，对这些斑块的分布进行高斯拟合，并定位其中心所在位置，通过重叠N帧B模式图像斑点变化信息，以获得N帧B模式图像中斑点变化的高分辨率合成图像，实现斑点追踪。通过斑点变化的高分辨率合成图像可以推断出待测区域中流体的流动速度及方向。
实施例2：
参照图7和图8，对于血流量小、血流速低的微小血管，存在显像困难的问题，为了提高此类组织的回波强度的同时不降低分辨力，本实施例基于实施例1，所提供超声显影装置中的编码发射模块用于通过脉冲编码序列，激励换能器晶片阵列向待测区域连续发射N组超声平面波束。所述接收解码模块用于对接收到的每组所述超声平面波束的回波进行与编码发射模块匹配的解码操作。
传统的单脉冲激励超声成像系统(参照图7)中，探头的换能器晶片收到脉冲激励发射超声波，图像的分辨率取决于换能器的脉冲响应。在编码激励成像系统(参照图8)中，采用持续时间较长的编码信号激励换能器，脉冲编码信号可以是一串长脉冲序列(如巴克码)，也可以是一定持续时间的连续信号(如线性调频信号)，由于脉冲编码信号的持续时间远长于换能器的脉冲响应时间，因此可以增加信号携带的能量，提高超声信号的平均功率。
需要解释的是，在通过相位调制或者频率调制的方式获得脉冲编码序列，以激励晶片阵列发射超声平面波束，同样地，在接收超声平面波束的回波时，需要使用对应的解调方式对所接收到的信号进行解调。
实施例3：
参照图5和图6，由于单一角度平面波没有聚焦效果，成像对比度和信噪比低，因此，本实施例在实施例1的基础上，提供超声显影装置中的编码发射模块，用于通过脉冲编码序列，激励换能器晶片阵列沿多个不同发射角度，向待测区域发射第n组超声平面波束中的多个超声平面波束，n为大于1且小于N的正整数。
需要解释的是，在一组超声平面波束中，包括多个超声平面波束，每个所述超声平面波束分别沿各自的发射角度发射，即，一组超声平面波束的回波包括多个超声平面波束的回波，每个所述超声平面波束的回波沿着其各自的反射角度或者散射角度进行回波。欲获得特定一帧的显影图像，需要将超声平面波束的回波所有超声平面波束的回波进行波束合成。即通过发射一系列不同偏转角度的平面波，将这些不同角度发射得到的回波数据进行波束合成，然后再相干叠加来提升最终一帧图像的对比度和信噪比。
具体地，将线阵探头置于待测区域上方(直接与待测区域介质表面接触)，超声显影装置中的波束合成模块，用于将在第n组超声平面波束中，多个超声平面波束的回波进行波束合成，获得第n帧显影图像，所述第n帧显影图像表达式为：
其中，x为第n帧显影图像像素点在与换能器晶片阵列平行方向上的位置；
z为第n帧显影图像像素点在待测区域深度方向上的位置；
x1为在所述第n帧显影图像像素点所对应介质点产生的回波，被换能器晶片阵列接收的位置。
需要解释的是，在所发射的超声平面波束覆盖空间范围内的组织能够以显影图像的形式进行显示，其中，像素点为显影图像的最小单位，介质点为超声平面波束覆盖空间范围内组织的最小单位，通过介质点的回声特性的不同，从而使得对应像素点的亮暗等级不同。
α为所述换能器晶片阵列发射超声平面波束的发射角度；
c为超声波在待测区域中的传播速度；
τ函数表示：换能器晶片阵列向第n帧显影图像像素点所对应介质点发出超声平面波束后，至接收到回波时的时间间隔函数；
RF函数表示：在所述第n帧显影图像像素点位置处产生的回波，被换能器晶片阵列接收形成的回波信号函数；
s函数表示：在特定发射角度下，所述待测区域的显影图像函数。
参照图5和图6，对于无偏转情况，及α＝0度，代入式为(3)得到
对于式(2)表示在单次超声平面波发射后，将换能器晶片阵列沿横向各个位置阵元获得的回波信号经过延迟τ(α，x1，x，z)后叠加起来。
将不同发射角度下得到的显影图像结果进行相干叠加，从而能够得到该帧最终显影图像，提升该帧最终显影图像的成像质量。其中，采用一系列不同的发射角度αi(i＝1，…m)，对每一个特定的发射角度，通过式(2)可得到该特定发射角度对应的显影图像，将在m个不同发射角度下输出的显影图像进行相干叠加，即可获得一帧最终的显影图像，即第n帧显影图像。
与传统聚焦波逐线扫描成像相比，超声平面波成像技术发射一次就可以获得一帧图像。将多个不同发射角度下的平面波发射获得的信号相干复合可以有效提升平面波图像的横向分辨率和信噪比。实验结果表明，9个角度的相干叠加即可达到与聚焦波相当的信噪比；17个角度的相干叠加可达到与聚焦波相当的横向分辨率，而与聚焦波要形成一帧图像所需的83次发射相比，平面波的发射次数显著减少。
多角度相干复合能够在保持帧率优势的前提下有效提升平面波成像的图像质量，使其能够应用于对时间分辨率要求较高的临床应用场景，基于多角度相干复合的超声平面波技术将进一步促进超声成像在各临床应用领域的发展。
实施例4：
本实施例在实施例1的基础上，提供超声显影装置中的所述斑点追踪模块，用于计算相邻帧B模式图像中对应像素点成像数据的成像方差；根据所述成像方差确定静态像素点和动态像素点；追踪所述动态像素点成像数据的变化信息，形成斑点变化信息。
具体地，对于本实施例中所述的“追踪所述动态像素点成像数据的变化信息，形成斑点变化信息”步骤，具体包括：
通过所述动态像素点及其对应的成像数据，确定动态斑点。
通过二维高斯曲面拟合确定位于所述动态斑点中心的中心动态像素点的坐标。
确定所述中心动态像素点成像数据的变化信息为所述动态斑点的斑点变化信息。
若需要获得第n帧B模式图像中的动态斑点的中心动态像素点坐标如图9至图11，其二维高斯曲面拟合函数为：
对式(4)两边取对数，并展开平方项，整理后得到：
若所述动态斑点参与拟合的动态像素点的有Z个，则将这Z个动态像素点写成矩阵的形式：A＝B C，其中：
A为Z*1的向量，其元素为：
ai＝fi×lnfi (其中i＝0、1、2……Z)；
B为Z*5的矩阵：
C为一个由高斯参数组成的向量：
对于式(4)和式(5)：
G为高斯分布的幅值；
x为第n帧B模式图像像素点在与换能器晶片阵列平行方向；
z为第n帧B模式图像像素点在待测区域深度方向上；
对于x0，其为所述中心动态像素点在x方向上的坐标；
对于z0，为所述中心动态像素点在z方向上的坐标；
对于σx，其为第n帧B模式图像各像素点在x方向上的标准差；
对于σz，其为第n帧B模式图像各个像素点在z方向上的标准差；
上述Z个参与拟合动态像素点误差的列向量为：E＝A-BC，用最小二乘法拟合，使其Z个动态像素点的均方差最小，即：
在图像数据处理时，像素点多数据量比较大，为减小计算量，将矩阵B进行QR分解，即：B＝QR，分解后Q为一个Z*Z的正交矩阵，R为一个Z*5的上三角矩阵，对E＝A-BC进行如下推导：
QTE＝QTA-QTBC＝QTA-RC (8)；
由于Q为正交矩阵，可以得到：
令：
对于以上式(7)至式(9)：其中，S为一个5维列向量；T为一个Z-5维列向量；R1为一个5*5的上三角方阵，则均方差：
对于上述式(10)，当S＝R1C时取得最小值，因此只需解出
即可求出CT中的x0,z0,σx,σz。
利用超声平面波束能够快速地对组织进行成像，通过在相邻帧之间较短时间里，明显移动会显示为半波长量级的斑点，对这些斑块的分布进行高斯拟合，并定位其中心所在位置，通过重叠N帧B模式图像斑点变化信息，以获得N帧B模式图像中斑点变化的高分辨率合成图像，实现斑点追踪。通过斑点变化的高分辨率合成图像可以推断出待测区域中流体的流动速度及方向。
实施例5：
本实施例在实施例1的基础上，提供超声显影装置中的所述成像处理模块还包括滤波模块，所述滤波模块用于依次对N帧B模式图像进行滤波处理，抑制斑点噪声；所述滤波模块将滤波处理后的N帧B模式图像输出给斑点追踪模块。
所述斑点追踪模块用于依次对抑制斑点噪声后的N帧B模式图像进行斑点追踪，确定在所述N帧B模式图像中的斑点变化信息。
所述图像合成模块用于将B模式图像斑点的斑点变化信息重迭显示在对应帧的B模式图像上。通过将斑点变化信息与B模式图像进行重迭合成，使得图像合成模块存储在图像存储器中的B模式超声图像序列包含斑点变化信息，以使得最终由显示模块显示的图像中能够显示出B模式图像的斑点变化信息。
所述滤波模块采用均值滤波器、中值滤波器、Lee滤波器、Sigma滤波器、Frost滤波器中的任意一种或多种对所述斑点变化信息进行滤波处理，抑制斑点噪声。
具体地，当所述成像处理模块通过采用Lee滤波器对所述斑点变化信息进行滤波处理，抑制斑点噪声时，使用5x5或者3x3或其他的滤波窗口内的局域均值和局域方差来估计斑点变化信息的先验均值和方差，从而实现对斑点变化信息进行滤波处理，抑制斑点噪声。
作为本发明的第二方面，提供一种超声显影方法。
实施方式1：
本实施方式提供的超声显影方法包括以下步骤：
S1：激励换能器晶片阵列向待测区域连续发射N组超声平面波束，N为大于1的正整数；
S2：接收每组所述超声平面波束的回波；
S3：对接收到的各组所述超声平面波束的回波进行波束合成，对应获得连续的N帧显影图像，并通过伪像去除算法对去除所述N帧显影图像中的伪像；通过对各组所述超声平面波束的回波进行波束合成，能够提升所获得的各帧显影图像的对比度和信噪比。
S4：依次对所述N帧显影图像进行B模式成像，获得N帧B模式图像；
S5：依次对N帧B模式图像进行斑点追踪，确定斑点在N帧B模式图像中的斑点变化信息；
S6：将所述斑点变化信息与N帧B模式图像进行合成，获得包含所述斑点变化信息的B模式超声图像序列。并将数据变换后得到的包含所述斑点变化信息的B模式超声图像序列存储在图像存储器中。显示模块通过扫描图像存储器以读出其中的数据进行显示。
传统的超声成像是通过发射多次聚焦波束照射组织来成像，以单波束接收为例，聚焦波每获得一根接收线需要付出一次发射的代价，如果一帧图像有M根扫描线，则要形成一帧图像共需发射M次；而在平面波成像模式下，同时激励换能器所有阵元，产生一个平行于换能器的超声平面波，此平面波产生的声场能够有效覆盖所有感兴趣区域，经组织散射后，所有阵元也同时接收回波信号，通过这些回波信号就可以得到整个区域的超声图像。因此平面波只需要进行一次发射即可得到一帧图像，和聚焦波相比，发射次数降为了1/M。即使聚焦波采用双波束、四波束等多波束并行处理技术，平面波在减少发射次数方面的优势仍然是非常可观的。
利用超声平面波束能够快速地对组织进行成像(优选每秒1000帧的速度)，通过在相邻帧之间的较短时间里，明显移动会显示为半波长量级的斑点，对这些斑块的分布进行高斯拟合，并定位其中心所在位置，通过重叠N帧B模式图像斑点变化信息，以获得N帧B模式图像中斑点变化的高分辨率合成图像，实现斑点追踪。通过斑点变化的高分辨率合成图像可以推断出待测区域中流体的流动速度及方向。
实施方式2：
本实施方式是在实施方式1的基础上，对于实施方式1中的步骤S1：激励换能器晶片阵列向待测区域连续发射N组超声平面波束，N为大于1的正整数，具体包括以下步骤：
通过脉冲编码序列，激励换能器晶片阵列向待测区域连续发射N组超声平面波束。在接收每组所述超声平面波束的回波时，对接收到的每组所述超声平面波束的回波进行匹配解码，其解码的算法与编码发射时的编码算法匹配。
由于对于血流量小、血流速低的微小血管，存在显像困难的问题，为了提高此类组织的回波强度的同时不降低分辨力。传统的单脉冲激励超声成像系统中，探头的换能器晶片收到脉冲激励发射超声波，图像的分辨率取决于换能器的脉冲响应。在编码激励成像系统中，采用持续时间较长的编码信号激励换能器，脉冲编码信号可以是一串长脉冲序列(如巴克码)，也可以是一定持续时间的连续信号(如线性调频信号)，由于脉冲编码信号的持续时间远长于换能器的脉冲响应时间，因此可以增加信号携带的能量，提高超声信号的平均功率。
需要解释的是，在通过相位调制或者频率调制的方式获得脉冲编码序列，以激励晶片阵列发射超声平面波束，同样地，在接收超声平面波束的回波时，需要使用对应的解调方式对所接收到的信号进行解调。
实施方式3：
本实施方式是在实施方式1的基础上，对于实施方式2中的步骤S1：激励换能器晶片阵列向待测区域连续发射N组超声平面波束，N为大于1的正整数，具体包括以下步骤：
通过脉冲编码序列，激励换能器晶片阵列沿多个不同发射角度，向待测区域发射第n组超声平面波束中的多个超声平面波束，n为大于1且小于N的正整数。
由于单一角度平面波没有聚焦效果，成像对比度和信噪比低，通过发射一系列不同偏转角度的平面波，将这些不同角度发射得到的回波数据进行波束合成，然后再相干叠加来提升最终一帧图像的对比度和信噪比。
实施方式4：
本实施方式是在实施方式3的基础上，对于实施方式3中的步骤S3：对接收到的各组所述超声平面波束的回波进行波束合成，对应获得连续的N帧显影图像，包括：将在第n组超声平面波束中，多个超声平面波束的回波进行波束合成，获得第n帧显影图像，所述第n帧显影图像表达式为：
其中，x为第n帧显影图像像素点在与换能器晶片阵列平行方向上的位置；
z为第n帧显影图像像素点在待测区域深度方向上的位置；
x1为在所述第n帧显影图像像素点所对应介质点产生的回波，被换能器晶片阵列接收的位置；
α为所述换能器晶片阵列发射超声平面波束的发射角度；
c为超声波在待测区域中的传播速度；
τ函数表示：换能器晶片阵列向第n帧显影图像像素点所对应介质点发出超声平面波束后，至接收到回波时的时间间隔函数；
RF函数表示：在所述第n帧显影图像像素点位置处产生的回波，被换能器晶片阵列接收形成的回波信号函数；
s函数表示：在特定发射角度下，所述待测区域的显影图像函数。
对接收到的各组所述超声平面波束的回波进行波束合成，对应获得连续的N帧显影图像；具体地：将线阵探头置于待测区域上方(直接与待测区域介质表面接触)，超声显影装置中的波束合成模块，用于将在第n组超声平面波束中，多个超声平面波束的回波进行波束合成，获得所述第n帧显影图像。
对于无偏转情况，及α＝0度，代入式为(3)得到
对于式(2)表示在单次超声平面波发射后，将换能器晶片阵列沿横向各个位置阵元获得的回波信号经过延迟τ(α，x1，x，z)后叠加起来。
将不同发射角度下得到的显影图像结果进行相干叠加，从而能够得到该帧最终显影图像，提升该帧最终显影图像的成像质量。其中，采用一系列不同的发射角度αi(i＝1，。。。m)，对每一个特定的发射角度，通过式(2)可得到该特定发射角度对应的显影图像，将在m个不同发射角度下输出的显影图像进行相干叠加，即可获得一帧最终的显影图像，即第n帧显影图像。
实施方式5：
本实施方式在实施方式1的基础上，对于实施方式1中的步骤S5：依次对N帧B模式图像进行斑点追踪，确定斑点在N帧B模式图像中的斑点变化信息，包括：
S51：计算相邻帧B模式图像中对应像素点成像数据的成像方差；
S52：根据所述成像方差确定静态像素点和动态像素点；
S53：追踪所述动态像素点成像数据的变化信息，形成斑点变化信息。
利用超声平面波束能够快速地对组织进行成像，通过在相邻帧之间较短时间里，明显移动会显示为半波长量级的斑点，对这些斑块的分布进行高斯拟合，并定位其中心所在位置，通过重叠N帧B模式图像斑点变化信息，以获得N帧B模式图像中斑点变化的高分辨率合成图像，实现斑点追踪。通过斑点变化的高分辨率合成图像可以推断出待测区域中流体的流动速度及方向。
显然，上述实施例仅仅是为清楚地说明所作的举例，而并非对实施方式的限定。对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。而由此所引伸出的显而易见的变化或变动仍处于本发明创造的保护范围之中。
Claims (10)
1.一种超声显影装置，其特征在于，包括：
编码发射模块，所述编码发射模块用于激励换能器晶片阵列向待测区域连续发射N组超声平面波束，N为大于1的正整数；
接收解码模块，所述接收解码模块用于接收每组所述超声平面波束的回波；
波束合成模块，所述波束合成模块用于对接收到的各组所述超声平面波束的回波进行波束合成，对应获得连续的N帧显影图像；
成像处理模块，所述成像处理模块包括B模式成像模块和斑点追踪模块；所述B模式成像模块，用于依次对所述N帧显影图像进行B模式成像，获得N帧B模式图像；所述斑点追踪模块用于依次对所述N帧B模式图像进行斑点追踪，确定在所述N帧B模式图像中的斑点变化信息；
图像合成模块，所述图像合成模块用于将所述斑点变化信息与N帧显影图像进行合成，获得包含所述斑点变化信息的B模式超声图像序列。
2.如权利要求1所述的超声显影装置，其特征在于，所述编码发射模用于通过脉冲编码序列，激励换能器晶片阵列向待测区域连续发射N组超声平面波束；
所述接收解码模块用于对接收到的每组所述超声平面波束的回波进行匹配解码。
3.如权利要求1所述的超声显影装置，其特征在于，所述编码发射模块，用于通过脉冲编码序列，激励换能器晶片阵列沿多个不同发射角度，向待测区域发射第n组超声平面波束中的多个超声平面波束，n为大于1且小于N的正整数。
4.如权利要求3所述的超声显影装置，其特征在于，所述波束合成模块，用于将在第n组超声平面波束中，多个超声平面波束的回波进行波束合成，获得第n帧显影图像，所述第n帧显影图像表达式为：
其中，x为第n帧显影图像像素点在与换能器晶片阵列平行方向上的位置；
z为第n帧显影图像像素点在待测区域深度方向上的位置；
x1为在所述第n帧显影图像像素点所对应介质点产生的回波，被换能器晶片阵列接收的位置；
α为所述换能器晶片阵列发射超声平面波束的发射角度；
m为大于1的正整数；
c为超声波在待测区域中的传播速度；
τ函数表示：换能器晶片阵列向第n帧显影图像像素点所对应介质点发出超声平面波后，至接收到回波时的时间间隔函数；
RF函数表示：在所述第n帧显影图像像素点位置处产生的回波，被换能器晶片阵列接收后形成的回波信号函数；
s函数表示：在特定发射角度下，所述待测区域的显影图像函数；
I函数表示：叠加在多个不同发射角度下所述待测区域的显影图像函数，获得的第n帧显影图像函数。
5.如权利要求1所述的超声显影装置，其特征在于，所述斑点追踪模块，用于计算相邻帧B模式图像中对应像素点成像数据的成像方差；根据所述成像方差确定静态像素点和动态像素点；追踪所述动态像素点成像数据的变化信息，形成斑点变化信息。
6.如权利要求1所述的超声显影装置，其特征在于，所述成像处理模块还包括滤波模块，所述滤波模块用于依次对N帧B模式图像进行滤波处理，抑制斑点噪声；
所述滤波模块将滤波处理后的N帧B模式图像输出给斑点追踪模块；
所述斑点追踪模块用于依次对抑制斑点噪声后的N帧B模式图像进行斑点追踪，确定在所述N帧显影图像中的斑点变化信息；
所述图像合成模块用于将显影图像斑点的斑点变化信息重迭显示在对应帧的B模式图像上。
7.如权利要求6所述的超声显影装置，其特征在于，所述滤波模块采用均值滤波器、中值滤波器、Lee滤波器、Sigma滤波器、Frost滤波器中的任意一种或多种对所述斑点变化信息进行滤波处理，抑制斑点噪声。
8.一种超声显影方法，其特征在于，包括以下步骤：
激励换能器晶片阵列向待测区域连续发射N组超声平面波束，N为大于1的正整数；
接收每组所述超声平面波束的回波；
对接收到的各组所述超声平面波束的回波进行波束合成，对应获得连续的N帧显影图像；
依次对所述N帧显影图像进行B模式成像，获得N帧B模式图像；
依次对所述N帧B模式图像进行斑点追踪，确定在所述N帧B模式图像中的斑点变化信息；
将所述斑点变化信息与N帧显影图像进行合成，获得包含所述斑点变化信息的B模式超声图像序列。
9.如权利要求8所述的超声显影方法，其特征在于，所述激励换能器晶片阵列向待测区域连续发射N组超声平面波束，包括：
通过脉冲编码序列，激励换能器晶片阵列向待测区域连续发射N组超声平面波束；
所述接收每组所述超声平面波束的回波步骤，还包括：对接收到的每组所述超声平面波束的回波进行匹配解码。
10.如权利要求9所述的超声显影方法，其特征在于，所述激励换能器晶片阵列向待测区域连续发射N组超声平面波束，包括：
通过脉冲编码序列，激励换能器晶片阵列沿多个不同发射角度，向待测区域发射第n组超声平面波束中的多个超声平面波束，n为大于1且小于N的正整数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010524116.0A CN113768533A (zh) | 2020-06-10 | 2020-06-10 | 超声显影装置和超声显影方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010524116.0A CN113768533A (zh) | 2020-06-10 | 2020-06-10 | 超声显影装置和超声显影方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113768533A true CN113768533A (zh) | 2021-12-10 |
Family
ID=78834693
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010524116.0A Pending CN113768533A (zh) | 2020-06-10 | 2020-06-10 | 超声显影装置和超声显影方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113768533A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116458925A (zh) * | 2023-06-15 | 2023-07-21 | 山东百多安医疗器械股份有限公司 | 一种便携式无盲区多模态超声心电系统 |
Citations (22)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6123670A (en) * | 1998-12-15 | 2000-09-26 | General Electric Company | Ultrasound imaging with optimal image quality in region of interest |
| US6181810B1 (en) * | 1998-07-30 | 2001-01-30 | Scimed Life Systems, Inc. | Method and apparatus for spatial and temporal filtering of intravascular ultrasonic image data |
| CN1413559A (zh) * | 2001-10-16 | 2003-04-30 | 株式会社东芝 | 计算关于局部血液流量循环指数的方法和装置 |
| JP2003250799A (ja) * | 2002-03-06 | 2003-09-09 | Shimadzu Corp | 超音波診断装置 |
| CN1892696A (zh) * | 2005-07-08 | 2007-01-10 | 深圳迈瑞生物医疗电子股份有限公司 | 超声图像边缘锐化与斑点抑制方法 |
| US20070049828A1 (en) * | 2005-08-29 | 2007-03-01 | Shenzhen Mindray Bio-Medical Electronics Co., Ltd. | Method and apparatus for ultrasonic imaging in mmode |
| WO2007069156A1 (en) * | 2005-12-14 | 2007-06-21 | Koninklijke Philips Electronics, N.V. | Transducer cuff for guidance and application of high intensity focused ultrasound for control of bleeding due to severed limbs |
| CN101330874A (zh) * | 2005-12-14 | 2008-12-24 | 皇家飞利浦电子股份有限公司 | 用于引导和施加高强度聚焦超声波以控制因截肢而出血的方法和设备 |
| JP2010082198A (ja) * | 2008-09-30 | 2010-04-15 | Fujifilm Corp | 超音波信号処理装置および方法 |
| CN101854537A (zh) * | 2009-03-31 | 2010-10-06 | 深圳迈瑞生物医疗电子股份有限公司 | 一种超声图像数据优化及其造影定量分析方法和系统 |
| US20120050074A1 (en) * | 2010-02-26 | 2012-03-01 | Bechtel Jon H | Automatic vehicle equipment monitoring, warning, and control system |
| CN102423265A (zh) * | 2011-09-28 | 2012-04-25 | 无锡祥生医学影像有限责任公司 | 超声诊断仪的复合成像方法 |
| CN103747742A (zh) * | 2011-04-14 | 2014-04-23 | 明尼苏达大学评议会 | 使用超声成像的脉管表征 |
| WO2016049681A1 (en) * | 2014-09-29 | 2016-04-07 | Signostics Limited | Ultrasound image processing system and method |
| CN106102588A (zh) * | 2015-09-06 | 2016-11-09 | 深圳迈瑞生物医疗电子股份有限公司 | 超声灰阶成像系统及方法 |
| US20160331335A1 (en) * | 2002-03-19 | 2016-11-17 | Medtronic Navigation, Inc. | Systems And Methods For Imaging Large Field-Of-View Objects |
| CN107403438A (zh) * | 2017-08-07 | 2017-11-28 | 河海大学常州校区 | 改进模糊聚类算法的超声图像病灶区分割方法 |
| US20180286199A1 (en) * | 2017-03-31 | 2018-10-04 | Qualcomm Incorporated | Methods and systems for shape adaptation for merged objects in video analytics |
| US20190008485A1 (en) * | 2017-07-09 | 2019-01-10 | The Board Of Trustees Of The Leland Stanford Junior University | Ultrasound imaging with spectral compounding for speckle reduction |
| CN110197713A (zh) * | 2019-05-10 | 2019-09-03 | 上海依智医疗技术有限公司 | 一种医疗影像的处理方法、装置、设备和介质 |
| CN110811688A (zh) * | 2019-12-02 | 2020-02-21 | 云南大学 | 多角度平面波重复复合的超快超声多普勒血流估计方法 |
| CN111012377A (zh) * | 2019-12-06 | 2020-04-17 | 北京安德医智科技有限公司 | 超声心动图心脏参数计算以及心肌应变测量方法、装置 |
-
2020
- 2020-06-10 CN CN202010524116.0A patent/CN113768533A/zh active Pending
Patent Citations (22)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6181810B1 (en) * | 1998-07-30 | 2001-01-30 | Scimed Life Systems, Inc. | Method and apparatus for spatial and temporal filtering of intravascular ultrasonic image data |
| US6123670A (en) * | 1998-12-15 | 2000-09-26 | General Electric Company | Ultrasound imaging with optimal image quality in region of interest |
| CN1413559A (zh) * | 2001-10-16 | 2003-04-30 | 株式会社东芝 | 计算关于局部血液流量循环指数的方法和装置 |
| JP2003250799A (ja) * | 2002-03-06 | 2003-09-09 | Shimadzu Corp | 超音波診断装置 |
| US20160331335A1 (en) * | 2002-03-19 | 2016-11-17 | Medtronic Navigation, Inc. | Systems And Methods For Imaging Large Field-Of-View Objects |
| CN1892696A (zh) * | 2005-07-08 | 2007-01-10 | 深圳迈瑞生物医疗电子股份有限公司 | 超声图像边缘锐化与斑点抑制方法 |
| US20070049828A1 (en) * | 2005-08-29 | 2007-03-01 | Shenzhen Mindray Bio-Medical Electronics Co., Ltd. | Method and apparatus for ultrasonic imaging in mmode |
| WO2007069156A1 (en) * | 2005-12-14 | 2007-06-21 | Koninklijke Philips Electronics, N.V. | Transducer cuff for guidance and application of high intensity focused ultrasound for control of bleeding due to severed limbs |
| CN101330874A (zh) * | 2005-12-14 | 2008-12-24 | 皇家飞利浦电子股份有限公司 | 用于引导和施加高强度聚焦超声波以控制因截肢而出血的方法和设备 |
| JP2010082198A (ja) * | 2008-09-30 | 2010-04-15 | Fujifilm Corp | 超音波信号処理装置および方法 |
| CN101854537A (zh) * | 2009-03-31 | 2010-10-06 | 深圳迈瑞生物医疗电子股份有限公司 | 一种超声图像数据优化及其造影定量分析方法和系统 |
| US20120050074A1 (en) * | 2010-02-26 | 2012-03-01 | Bechtel Jon H | Automatic vehicle equipment monitoring, warning, and control system |
| CN103747742A (zh) * | 2011-04-14 | 2014-04-23 | 明尼苏达大学评议会 | 使用超声成像的脉管表征 |
| CN102423265A (zh) * | 2011-09-28 | 2012-04-25 | 无锡祥生医学影像有限责任公司 | 超声诊断仪的复合成像方法 |
| WO2016049681A1 (en) * | 2014-09-29 | 2016-04-07 | Signostics Limited | Ultrasound image processing system and method |
| CN106102588A (zh) * | 2015-09-06 | 2016-11-09 | 深圳迈瑞生物医疗电子股份有限公司 | 超声灰阶成像系统及方法 |
| US20180286199A1 (en) * | 2017-03-31 | 2018-10-04 | Qualcomm Incorporated | Methods and systems for shape adaptation for merged objects in video analytics |
| US20190008485A1 (en) * | 2017-07-09 | 2019-01-10 | The Board Of Trustees Of The Leland Stanford Junior University | Ultrasound imaging with spectral compounding for speckle reduction |
| CN107403438A (zh) * | 2017-08-07 | 2017-11-28 | 河海大学常州校区 | 改进模糊聚类算法的超声图像病灶区分割方法 |
| CN110197713A (zh) * | 2019-05-10 | 2019-09-03 | 上海依智医疗技术有限公司 | 一种医疗影像的处理方法、装置、设备和介质 |
| CN110811688A (zh) * | 2019-12-02 | 2020-02-21 | 云南大学 | 多角度平面波重复复合的超快超声多普勒血流估计方法 |
| CN111012377A (zh) * | 2019-12-06 | 2020-04-17 | 北京安德医智科技有限公司 | 超声心动图心脏参数计算以及心肌应变测量方法、装置 |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116458925A (zh) * | 2023-06-15 | 2023-07-21 | 山东百多安医疗器械股份有限公司 | 一种便携式无盲区多模态超声心电系统 |
| CN116458925B (zh) * | 2023-06-15 | 2023-09-01 | 山东百多安医疗器械股份有限公司 | 一种便携式无盲区多模态超声心电系统 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Jensen et al. | Ultrasound vector flow imaging—Part II: Parallel systems | |
| US8111810B2 (en) | Method for producing highly constrained ultrasound images | |
| KR101820422B1 (ko) | 포커싱되지 않은 송신 빔들을 이용한 높은 프레임 레이트의 정량적 도플러 흐름 촬영 | |
| US7347820B2 (en) | Phased array acoustic system for 3D imaging of moving parts | |
| JP4297555B2 (ja) | 超音波カラー・ドップラー速度／方向撮像 | |
| RU2480147C2 (ru) | Комбинированная система фотоакустического и ультразвукового формирования изображений | |
| JP4730125B2 (ja) | 血流画像表示装置 | |
| Forsberg | Ultrasonic biomedical technology; marketing versus clinical reality | |
| JP5341352B2 (ja) | 超音波画像化システムのクラッター抑制本出願は、２００４年１２月３０日に出願された米国仮特許出願Ｎｏ．６０／６４０，３６８に基づき優先権の利益を主張する。本出願は、２００４年１月７日に出願した米国仮特許出願Ｎｏ．６０／５３４，３９０に関連し、明細はここに参照のため組み込まれている。 | |
| US20110144495A1 (en) | Perfusion Imaging of a Volume in Medical Diagnostic Ultrasound | |
| JP2009505771A (ja) | 実時間空間合成によるフローイメージングのための超音波画像形成システム及び方法 | |
| JP2014512217A (ja) | 超音波イメージング方法及び超音波イメージング装置 | |
| JP2010523253A (ja) | 高速超音波による厚いスライスのイメージング | |
| JP2005177494A (ja) | フローパラメータイメージングのための方法及び装置 | |
| US20220292637A1 (en) | Methods for High Spatial and Temporal Resolution Ultrasound Imaging of Microvessels | |
| EP2134262A2 (en) | Quantified perfusion studies with ultrasonic thick slice imaging | |
| KR20190087041A (ko) | 초음파 영상장치 및 그 제어방법 | |
| WO2019179758A1 (en) | Ultrasound system for shear wave imaging in three dimensions | |
| JP2021522011A (ja) | 超音波スクリーニングのためのシステム及び方法 | |
| US8235906B2 (en) | System and method for accelerated focused ultrasound imaging | |
| CN113768533A (zh) | 超声显影装置和超声显影方法 | |
| Jakovljevic et al. | Short-lag spatial coherence imaging on matrix arrays, Part II: Phantom and in vivo experiments | |
| Hyun et al. | A GPU-based real-time spatial coherence imaging system | |
| WO2008127927A1 (en) | Tissue border detection in ultrasonic thick slice imaging | |
| JP7345678B2 (ja) | ３ｄベクトルフロー場を取得するための方法及びシステム |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |