CN107833218A - 一种生物组织剪切波弹性图像质量评估显示方法 - Google Patents
一种生物组织剪切波弹性图像质量评估显示方法 Download PDFInfo
- Publication number
- CN107833218A CN107833218A CN201711175313.0A CN201711175313A CN107833218A CN 107833218 A CN107833218 A CN 107833218A CN 201711175313 A CN201711175313 A CN 201711175313A CN 107833218 A CN107833218 A CN 107833218A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- elastic image
- region
- shearing wave
- 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.)
- Granted
Links
- 238000010008 shearing Methods 0.000 title claims abstract description 41
- 238000000034 method Methods 0.000 title claims abstract description 28
- 238000013441 quality evaluation Methods 0.000 title claims abstract description 16
- 238000006073 displacement reaction Methods 0.000 claims abstract description 52
- 238000005259 measurement Methods 0.000 claims abstract description 21
- 230000007547 defect Effects 0.000 claims abstract description 20
- 238000012360 testing method Methods 0.000 claims description 8
- 238000001914 filtration Methods 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000003384 imaging method Methods 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 3
- 238000002604 ultrasonography Methods 0.000 claims description 2
- 239000003086 colorant Substances 0.000 abstract description 2
- 230000000694 effects Effects 0.000 abstract description 2
- 230000000007 visual effect Effects 0.000 abstract 1
- 230000005855 radiation Effects 0.000 description 4
- 238000013461 design Methods 0.000 description 3
- 238000001514 detection method Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000003908 quality control method Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 206010028980 Neoplasm Diseases 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 210000004204 blood vessel Anatomy 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 239000003989 dielectric material Substances 0.000 description 1
- 238000002091 elastography Methods 0.000 description 1
- 230000003116 impacting effect Effects 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 230000035807 sensation Effects 0.000 description 1
- 238000009331 sowing Methods 0.000 description 1
Classifications
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
-
- 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/30168—Image quality inspection
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Medical Informatics (AREA)
- Quality & Reliability (AREA)
- Radiology & Medical Imaging (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Geometry (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本发明公开一种生物组织剪切波弹性图像质量评估显示方法,该方法通过计算剪切波位移信噪比分布情况,确定弹性图像缺损情况,再结合弹性图像均匀度,综合判断本次弹性测量图像的质量和有效性,最后以彩色条直观表示本次弹性测量效果质量,为医生提供直观的测量可靠性量化指标。本发明通过综合考虑剪切波位移估计信噪比,结合最终弹性图像均匀度,最终给出本次弹性测量有效性的判断并通过质量系数和彩色进度条方式显示,而质量的优中差分别对应绿黄红三种颜色,与人类生活习惯对应,使用者能够通过进度条结果快速判断此次弹性测量结果的有效性。
Description
技术领域
本发明涉及剪切波弹性成像领域,特别涉及一种生物组织剪切波弹性图像质量评估显示方法。
背景技术
剪切波是传播方向与介质质点的振动方向垂直的波,又称横波,S波。
剪切波弹性成像技术可实现生物组织实时硬度定量检测,为临床判断组织病变情况提供依据。它的基本原理如下:声辐射力聚焦冲击能在组织内部产生剪切波,由于剪切波在不同硬度的组织传播速度有差异,通过检测剪切波在不同位置的传播速度可以间接反映该位置的软硬情况。根据剪切波波速大小进行伪彩色映射即为剪切波弹性成像。
剪切波由声辐射力冲击产生,冲击时间非常短,与静压弹性成像相比大大降低了临床操作手法的难度;而冲击深度由超声聚焦深度决定,比静压弹性成像使用范围更加广。但剪切波在复杂的人体组织中传播时,衰减非常快;而且不同介质剪切波质量差异很大,所以弹性图像的质量与ROI(Region Of Interest感兴趣区域,也称为AOI area-of-interest)大小和选取位置密切相关。目前临床一般通过弹性图像的充盈程度或者质控图情况判断弹性测量的有效性,但充盈程度的判断和质控图的好坏最终还是依赖于医生赁经验主观判断,这种方式不能直观反映人体组织中传播剪切波时的弹性图像,这样在诊断过程中可能会带来较大的误差。
发明内容
为了能够更直观的表示弹性图像的质量情况,本发明提出一种生物组织剪切波弹性图像质量评估显示的方法。
本发明为实现其目的所采用的技术方案是:一种生物组织剪切波弹性图像质量评估显示方法,包括以下步骤:
步骤1、对待测量弹性区域中感兴趣区域发射聚焦冲击声束;
步骤2、对待测量弹性区域中感兴趣区域发射采集检测声束;
步骤3、对同一采集位置多次采集回来的信号进行位移估计运算;
步骤4、对位移估计结果进行低通滤波;
步骤5、将滤波后的位移估计结果与滤波前位移估计结果相减得到位移估计噪声结果;
步骤6、计算滤波后的位移估计结果的平方和与位移估计噪声的平方和结果的比值,得到位移估计结果信噪比结果;
步骤7、搜索并记录位移估计结果信噪比小于预设阈值的位置;确定为弹性缺损位置;
步骤8、计算缺损位置面积和弹性测量区域总面积之比值;
步骤9、对待测量弹性区域中感兴趣区域中进行剪切波波速估计;
步骤10、计算待测量弹性区域内所有剪切波波速值的标准差值与均值之比值,得到待测量弹性区域均匀度;
步骤11、根据缺损面积比值和待测量弹性均匀度确定质量系数;
步骤12、根据质量系数量化显示到界面,显示方式为进度条模式。
本发明通过弹性图像质量评估显示功能,为医生快速判断测量结果的有效性提供了直观的依据。
进一步的,上述的生物组织剪切波弹性图像质量评估显示方法中:所述的步骤1中,聚焦冲击声束波形选择与传统B成像模式发射波形,聚焦位置选择待测量弹性区域中感兴趣区域外的一个位置。
进一步的,上述的生物组织剪切波弹性图像质量评估显示方法中:所述的步骤2中,采集检测声束位置选择待测量弹性区域中感兴趣区域内两线或多线所在位置,线间隔1~3mm可选,每线重复采集次数50~200次可选,每线数据采集深度范围选择与待测量弹性区域中感兴趣区域上下边界深度范围一致。
进一步的,上述的生物组织剪切波弹性图像质量评估显示方法中:所述的步骤3中,位移估计方法采用归一化互相关时延估计法;按照下式计算得到被测区域内各点在深度方向上的位移uτ
式中:s为超声射频回波信号,R为互相关系数,τ为对应位移值的超声回波时间延迟,c为超声传播速度。
进一步的,上述的生物组织剪切波弹性图像质量评估显示方法中:所述步骤6中,按照下列公式得到位移估计结果信噪比结果:
式中:Vel_Filter为滤波后位移估计结果矩阵,Vel_Noise为位移估计噪声矩阵,VelSNR是位移估计结果信噪比结果。
进一步的,上述的生物组织剪切波弹性图像质量评估显示方法中:所述的步骤8中,计算缺损位置面积和弹性测量区域总面积之比值就是缺损位置点数个数与待测量弹性区域中感兴趣区域总点数个数之比值;计算公式如下:
式中:Num(Defect)为缺损位置点数个数,Num(ROI)为待测量弹性区域中感兴趣区域总点数个数,Qs是比值。
下面结合具体实施例对本发明作较为详细的描述。
附图说明
图1是本发明实施例1流程图。
图2剪切波位移估计结果。
图3剪切波位移信噪比。
图4各位置“形变-时间”曲线。
图5拟合“时间-距离”直线。
图6弹性质量评估显示。
具体实施方式
本发明实施例通过计算剪切波位移信噪比分布情况,确定弹性图像缺损情况,再结合弹性图像均匀度,综合判断本次弹性测量图像的质量和有效性,最后以彩色条直观表示本次弹性测量效果质量,为医生提供直观的测量可靠性量化指标。剪切波位移信噪比的计算可以反映本次测量中,剪切波在组织介质中传播情况。一般在较为均匀的介质中传播时,剪切波位移信噪比有较大值;但当剪切波经过组织介质复杂如血管,囊肿等组织时,剪切波位移信噪比比较差,此时波速估计结果有效性很低。而弹性图像均匀度能反映组织硬度的部分变化情况,一般组织硬度变化存在一定的过渡性。当弹性图像均匀度较差时,如很多局部小区域硬度变化或变化差异很大时,可能是由于剪切波在复杂介质边界传播时产生分裂等情况导致,此时弹性图像有效性低。综合考虑以上因素,结合这两个特征分析该次弹性测量有效性。具体包括以下步骤:
步骤1:对待测量弹性感兴趣区域(ROI)区域发射聚焦冲击声束。发射波形根据探头中心频率确定,可选中心频率3~5MHz,聚焦位置可选ROI区域外的一个位置,例如:ROI左边界深度中点位置。冲击时长可选例如:100us。
步骤2:待聚焦冲击声束发射完毕后,发射采集检测声束,此方式与传统B成像模式一致,采集声束位置可选ROI内两线或多线所在位置,线间隔可选1~3mm,每线重复采集次数可选50~200次。每线数据采集深度范围可选与ROI上下边界深度范围一致。确定待观察深度的位置,选取两个或多个横向检测位置,根据形变估计结果矩阵作出所选位置对应的“形变-时间”曲线,如图4所示。图中显示了ROI区中五个不同的横向检测位置的“形变-时间”曲线,虽然位置不同但波形基本相似,频率一致,只是波峰到来的时刻不同,另外,可以看出,声辐射力冲击是在50ms。具体聚焦冲击声束和采集检测声束发射模式参考文献:IEEETransactions on Ultrasonics,Ferroelectrics,and Frequency Control,vol.60,no.4,April 2013《Acoustic Radiation Force Elasticity Imaging in DiagnosticUltrasound》。
步骤3:对同一采集位置多次采集回来的信号进行位移估计运算。发射完采集声束后探头接收信号就是采集的过程,这是传统B成像原理,本领域技术人员通过上面的技术可以自动完成。位移估计方法可选例如:归一化互相关时延估计法。通过归一化互相关时延估计法计算得到被测区域内各点在深度方向上的位移uτ(其中s为超声射频回波信号,R为互相关系数,τ为对应位移值的超声回波时间延迟,c为超声传播速度),结果如附图2:图中标记1指示处颜色浅的地方表示波峰。
步骤4:对位移估计结果进行低通滤波,滤波器可选如:FIR类型滤波器,截止频率可选如:1000hz。
步骤5:将滤波后的位移估计结果与滤波前位移估计结果相减得到位移估计噪声结果。
步骤6:计算滤波后的位移估计结果的平方和与位移估计噪声的平方和结果的比值,得到位移估计结果信噪比结果。公式如下:(其中Vel_Filter为滤波后位移估计结果矩阵,Vel_Noise为位移估计噪声矩阵),结果如附图3所示,图中2指示处颜色浅的地方表示信噪比小的地方。
步骤7:搜索并记录位移估计信噪比小于预设阈值的位置。确定为弹性缺损位置。
步骤8:计算缺损位置面积和弹性测量区域总面积之比值(即缺损位置点数个数与ROI总点数个数之比值)。计算公式可如下:
步骤9:剪切波波速估计,估计方法可选用剪切波波速估计经典方法TOF(Time OfFlight)法。TOF原理是选取两个或多个不同检测线位置,分别作出同一深度,不同检测线求取直线斜率即为剪切波波速估计结果,如附图5所示。
步骤10:计算弹性区域内所有剪切波波速值的标准差值与均值之比值,得到弹性区域均匀度,计算公式可如下:
其中Emean表示ROI区域剪切波波速值均值,xi为ROI区域中各位置剪切波波速值,n为剪切波数量。
步骤11:根据缺损面积比值和弹性均匀度确定质量系数,两个参考指标权重可设置为0.5,如最终质量系数计算方法可如下:
QE=1-0.5QS-0.5QI
根据上述方式计算的质量系数,质量系数越大表示弹性图像质量越高,可信度越大。反之代表弹性图像质量低,可信度小。权重值具体可根据不同组织结构特点作出调整。
根据12:根据质量系数量化显示到界面,显示方式为进度条模式。例如:
1)质量优:QE=76~100,进度条充满程度76%~100%,显示颜色绿色。
2)质量良:QE=51~75,进度条充满程度51%~75%,显示颜色橙色。
3)质量中:QE=26~50,进度条充满程度26%~50%,显示颜色黄色。
4)质量优:QE=0~25,进度条充满程度0%~25%,显示颜色红色。如附图6所示。
剪切波弹性成像的帧频一般比较低,而且测量的结果有效性存在不稳定性,在临床应用中,需要靠医生对图像的情况作出主观判断来确定该次测量的有效性,使得应用起来不方便,而且不同使用者之间也存在判断差异。但通过本实施例中,弹性图像质量评估显示功能,为医生快速判断测量结果的有效性提供了直观的依据。本发明通过综合考虑剪切波位移估计信噪比,结合最终弹性图像均匀度,最终给出本次弹性测量有效性的判断并通过质量系数和彩色进度条方式显示,而质量的优中差分别对应绿黄红三种颜色,与人类生活习惯对应,使用者能够通过进度条结果快速判断此次弹性测量结果的有效性。
Claims (6)
1.一种生物组织剪切波弹性图像质量评估显示方法,其特征在于包括以下步骤:
步骤1、对待测量弹性区域中感兴趣区域发射聚焦冲击声束;
步骤2、对待测量弹性区域中感兴趣区域发射采集检测声束;
步骤3、对同一采集位置多次采集回来的信号进行位移估计运算;
步骤4、对位移估计结果进行低通滤波;
步骤5、将滤波后的位移估计结果与滤波前位移估计结果相减得到位移估计噪声结果;
步骤6、计算滤波后的位移估计结果的平方和与位移估计噪声的平方和结果的比值,得到位移估计结果信噪比结果;
步骤7、搜索并记录位移估计结果信噪比小于预设阈值的位置;确定为弹性缺损位置;
步骤8、计算缺损位置面积和弹性测量区域总面积之比值;
步骤9、对待测量弹性区域中感兴趣区域中进行剪切波波速估计;
步骤10、计算待测量弹性区域内所有剪切波波速值的标准差值与均值之比值,得到待测量弹性区域均匀度;
步骤11、根据缺损面积比值和待测量弹性均匀度确定质量系数;
步骤12、根据质量系数量化显示到界面,显示方式为进度条模式。
2.根据权利要求1所述的生物组织剪切波弹性图像质量评估显示方法,其特征在于:所述的步骤1中,聚焦冲击声束波形选择与传统B成像模式发射波形,聚焦位置选择待测量弹性区域中感兴趣区域外的一个位置。
3.根据权利要求1所述的生物组织剪切波弹性图像质量评估显示方法,其特征在于:所述的步骤2中,采集检测声束位置选择待测量弹性区域中感兴趣区域内两线或多线所在位置,线间隔1~3mm可选,每线重复采集次数50~200次可选,每线数据采集深度范围选择与待测量弹性区域中感兴趣区域上下边界深度范围一致。
4.根据权利要求1所述的生物组织剪切波弹性图像质量评估显示方法,其特征在于:所述的步骤3中,位移估计方法采用归一化互相关时延估计法;按照下式计算得到被测区域内各点在深度方向上的位移uτ
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>n</mi>
<mi>c</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msubsup>
<mo>&Integral;</mo>
<mrow>
<mo>-</mo>
<mi>T</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
<mrow>
<mi>T</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mi>r</mi>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>s</mi>
<mi>d</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>+</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
<mo>)</mo>
<mi>d</mi>
<mi>t</mi>
</mrow>
<msqrt>
<mrow>
<msubsup>
<mo>&Integral;</mo>
<mrow>
<mo>-</mo>
<mi>T</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
<mrow>
<mi>T</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mi>r</mi>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mi>d</mi>
<mi>t</mi>
<msubsup>
<mo>&Integral;</mo>
<mrow>
<mo>-</mo>
<mi>T</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
<mrow>
<mi>T</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mi>d</mi>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>+</mo>
<mi>&tau;</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mi>d</mi>
<mi>t</mi>
</mrow>
</msqrt>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>u</mi>
<mi>&tau;</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>c</mi>
<mi>&tau;</mi>
</mrow>
<mn>2</mn>
</mfrac>
</mrow>
式中:Sr和Sd为进行互相关超声射频回波信号,R为互相关系数,τ为对应位移值的超声回波时间延迟,c为超声传播速度,T为进行互相关运算的信号的时间长度。
5.根据权利要求1所述的生物组织剪切波弹性图像质量评估显示方法,其特征在于:所述步骤6中,按照下列公式得到位移估计结果信噪比结果:
<mrow>
<msub>
<mi>Vel</mi>
<mrow>
<mi>S</mi>
<mi>N</mi>
<mi>R</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msup>
<msub>
<mi>&Sigma;Vel</mi>
<mrow>
<mi>F</mi>
<mi>i</mi>
<mi>l</mi>
<mi>t</mi>
<mi>e</mi>
<mi>r</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<msup>
<msub>
<mi>&Sigma;Vel</mi>
<mrow>
<mi>N</mi>
<mi>o</mi>
<mi>i</mi>
<mi>s</mi>
<mi>e</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
</mrow>
式中:Vel_Filter为滤波后位移估计结果矩阵,Vel_Noise为位移估计噪声矩阵,VelSNR是位移估计结果信噪比结果。
6.根据权利要求1所述的生物组织剪切波弹性图像质量评估显示方法,其特征在于:所述的步骤8中,计算缺损位置面积和弹性测量区域总面积之比值就是缺损位置点数个数与待测量弹性区域中感兴趣区域总点数个数之比值;计算公式如下:
<mrow>
<msub>
<mi>Q</mi>
<mi>S</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>N</mi>
<mi>u</mi>
<mi>m</mi>
<mrow>
<mo>(</mo>
<mi>D</mi>
<mi>e</mi>
<mi>f</mi>
<mi>e</mi>
<mi>c</mi>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>N</mi>
<mi>u</mi>
<mi>m</mi>
<mrow>
<mo>(</mo>
<mi>R</mi>
<mi>O</mi>
<mi>I</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
式中:Num(Defect)为缺损位置点数个数,Num(ROI)为待测量弹性区域中感兴趣区域总点数个数,Qs是比值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711175313.0A CN107833218B (zh) | 2017-11-22 | 2017-11-22 | 一种生物组织剪切波弹性图像质量评估显示方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711175313.0A CN107833218B (zh) | 2017-11-22 | 2017-11-22 | 一种生物组织剪切波弹性图像质量评估显示方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107833218A true CN107833218A (zh) | 2018-03-23 |
CN107833218B CN107833218B (zh) | 2021-01-26 |
Family
ID=61653340
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711175313.0A Active CN107833218B (zh) | 2017-11-22 | 2017-11-22 | 一种生物组织剪切波弹性图像质量评估显示方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107833218B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109978823A (zh) * | 2019-02-15 | 2019-07-05 | 无锡海斯凯尔医学技术有限公司 | 弹性成像图像分析及组织粘弹性检测方法及装置 |
WO2020037673A1 (zh) * | 2018-08-24 | 2020-02-27 | 深圳迈瑞生物医疗电子股份有限公司 | 一种超声弹性成像装置及对弹性图像进行处理的方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101708123A (zh) * | 2009-10-28 | 2010-05-19 | 上海理工大学 | 肝纤维化分级研究的磁共振弹性成像检测系统及其方法 |
CN101912278A (zh) * | 2010-08-12 | 2010-12-15 | 陈庆武 | 超声动态弹性成像探头及方法 |
CN103079473A (zh) * | 2010-08-31 | 2013-05-01 | 株式会社日立医疗器械 | 超声波诊断装置以及评价算出方法 |
CN103462643A (zh) * | 2013-09-29 | 2013-12-25 | 深圳市开立科技有限公司 | 一种剪切波速度测量方法、装置和系统 |
CN103845074A (zh) * | 2012-11-28 | 2014-06-11 | 深圳迈瑞生物医疗电子股份有限公司 | 一种超声弹性成像系统和方法 |
CN104055541A (zh) * | 2014-06-26 | 2014-09-24 | 中国科学院苏州生物医学工程技术研究所 | 一种用于血管内超声多断层剪切波弹性成像方法 |
CN104939869A (zh) * | 2014-03-31 | 2015-09-30 | 美国西门子医疗解决公司 | 用于弹性超声成像的获取控制的系统和方法 |
CN107205724A (zh) * | 2014-10-28 | 2017-09-26 | 超声成像 | 成像方法以及进行剪切波弹性成像的装置 |
-
2017
- 2017-11-22 CN CN201711175313.0A patent/CN107833218B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101708123A (zh) * | 2009-10-28 | 2010-05-19 | 上海理工大学 | 肝纤维化分级研究的磁共振弹性成像检测系统及其方法 |
CN101912278A (zh) * | 2010-08-12 | 2010-12-15 | 陈庆武 | 超声动态弹性成像探头及方法 |
CN103079473A (zh) * | 2010-08-31 | 2013-05-01 | 株式会社日立医疗器械 | 超声波诊断装置以及评价算出方法 |
CN103845074A (zh) * | 2012-11-28 | 2014-06-11 | 深圳迈瑞生物医疗电子股份有限公司 | 一种超声弹性成像系统和方法 |
CN103462643A (zh) * | 2013-09-29 | 2013-12-25 | 深圳市开立科技有限公司 | 一种剪切波速度测量方法、装置和系统 |
CN104939869A (zh) * | 2014-03-31 | 2015-09-30 | 美国西门子医疗解决公司 | 用于弹性超声成像的获取控制的系统和方法 |
CN104055541A (zh) * | 2014-06-26 | 2014-09-24 | 中国科学院苏州生物医学工程技术研究所 | 一种用于血管内超声多断层剪切波弹性成像方法 |
CN107205724A (zh) * | 2014-10-28 | 2017-09-26 | 超声成像 | 成像方法以及进行剪切波弹性成像的装置 |
Non-Patent Citations (1)
Title |
---|
曾博等: "基于CUDA的声辐射力弹性成像算法研究", 《计算机工程与应用》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2020037673A1 (zh) * | 2018-08-24 | 2020-02-27 | 深圳迈瑞生物医疗电子股份有限公司 | 一种超声弹性成像装置及对弹性图像进行处理的方法 |
CN112534468A (zh) * | 2018-08-24 | 2021-03-19 | 深圳迈瑞生物医疗电子股份有限公司 | 一种超声弹性成像装置及对弹性图像进行处理的方法 |
CN109978823A (zh) * | 2019-02-15 | 2019-07-05 | 无锡海斯凯尔医学技术有限公司 | 弹性成像图像分析及组织粘弹性检测方法及装置 |
CN109978823B (zh) * | 2019-02-15 | 2021-08-06 | 无锡海斯凯尔医学技术有限公司 | 弹性成像图像分析及组织粘弹性检测方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN107833218B (zh) | 2021-01-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Hindi et al. | Artifacts in diagnostic ultrasound | |
CN104013419B (zh) | 医疗超声波中的自适应声压估计 | |
Needles et al. | Nonlinear contrast imaging with an array-based micro-ultrasound system | |
Mattoon et al. | Fundamentals of diagnostic ultrasound | |
Kim et al. | Attenuation estimation using spectral cross-correlation | |
US8758248B2 (en) | Systems and methods for acoustic radiation force imaging with enhanced performance | |
JP2012508053A (ja) | 時間的パラメトリック・コントラスト画像を使った超音波病変識別 | |
US20140046606A1 (en) | Methods and apparatus for multibeam doppler ultrasound display | |
CN105407809A (zh) | 超声波诊断装置以及弹性评价方法 | |
EP3769690B1 (en) | Method for tissue characterization by ultrasound wave attenuation measurements and ultrasound system for tissue characterization | |
Mast et al. | The effect of abdominal wall morphology on ultrasonic pulse distortion. Part II. Simulations | |
Rivera et al. | Quantification of tricuspid regurgitation by means of the proximal flow convergence method: a clinical study | |
PRICE et al. | Basic concepts of ultrasonic tissue characterization | |
CN107833218A (zh) | 一种生物组织剪切波弹性图像质量评估显示方法 | |
Ahmed et al. | Plane-wave imaging improves single-track location shear wave elasticity imaging | |
Morris et al. | Multiparametric ultrasound for targeting prostate cancer: combining ARFI, SWEI, QUS and B-mode | |
WO2021101635A1 (en) | Additional diagnostic data in parametric ultrasound medical imaging | |
Kutay et al. | On modeling biomedical ultrasound RF echoes using a power-law shot-noise model | |
Rivera et al. | Visual assessment of valvular regurgitation: comparison with quantitative Doppler measurements | |
EP3694418B1 (en) | Ultrasonic shearwave imaging with patient-adaptive shearwave generation | |
CN108392751B (zh) | 一种实时监测高强聚焦超声治疗声空化的方法 | |
CN107049252A (zh) | 一种生物磁光声联合内窥成像方法 | |
Bigelow | Estimating the total ultrasound attenuation along the propagation path by applying multiple filters to backscattered echoes from a single spherically focused source | |
Urban et al. | Harmonic motion detection in a vibrating scattering medium | |
Klimonda et al. | Tissue attenuation estimation by mean frequency downshift and bandwidth limitation |
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 | ||
CB02 | Change of applicant information | ||
CB02 | Change of applicant information |
Address after: 405, building 9, baiwangxin Industrial Park, No. 1002, Songbai Road, Yangguang community, Xili street, Nanshan District, Shenzhen City, Guangdong Province Applicant after: Shenzhen Zhongke Medical Technology Co., Ltd. Address before: Baoan District Shiyan street Shenzhen city Guangdong province 518000 White Pine Road in Yuntai technology industrial factory building 9 building seven floor East Applicant before: Shenzhen Zhongke Medical Technology Co., Ltd. |
|
GR01 | Patent grant | ||
GR01 | Patent grant |