《_弹性成像技术_超声图像处理与分析_陆敏华.pdf》由会员分享，可在线阅读，更多相关《_弹性成像技术_超声图像处理与分析_陆敏华.pdf（7页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、doi:103969/j issn1674 1242201302 005医学超声关键技术研究和进展(二)超声弹性成像技术、超声图像处理与分析陆敏华1，2，3，陈昕1，2，3，倪东1，2，3#，张灵1，2，3，林浩铭1，2，3，汪天富1，2，3，陈思平1，2，3*1 医学超声关键技术国家地方联合工程实验室2 广东省生物医学信息检测与超声成像重点实验室3 深圳大学医学院生物医学工程系，广东省深圳市，518060作者简介:陆敏华，陈昕，倪东(共同第一作者)通信作者:陈思平，E mail:chensiping szu edu cn【摘要】超声弹性成像技术是近 20 年来在医学超声领域出现的新技术，它利
2、用准静态压力、外加低频振动、声辐射力等形式激励组织振动，然后用超声方法检测组织的形变或者由振动引起的剪切波传播，最终得到组织的生物力学特性。超声弹性成像已经被广泛应用在乳腺、肝脏、前列腺等领域，显示了良好的临床应用前景。此外，超声图像处理与分析是超声成像系统的重要组成部分，主要包括超声图像的三维可视化、预处理、分割、配准及图像理解等，旨在为临床医生诊断提供有效的辅助手段。【关键词】超声弹性成像，粘弹性，超声图像处理与分析【中图分类号】445 1【文献标志码】A文章编号:1674 1242(2013)02 0082 07Progress in the esearch on Key Technol
3、ogies for Medical Ultrasound，Part Two:Ultrasound Elastography，Ultrasound Image Processing and AnalysisLu Minhua1，2，3，Chen Xin1，2，3，Ni Dong1，2，3#，Zhang Ling1，2，3，Lin Haoming1，2，3，Wang Tianfu1，2，3，Chen Siping1，2，3*1 National egional Key Technology Engineering Laboratory for Medical Ultrasound(Shenzhen
4、，518060)2 Guangdong Key Laboratory for Biomedical Measurements and Ultrasound Imaging(Shenzhen，518060)3 Department of Biomedical Engineering，School of Medicine，Shenzhen University(Shenzhen，518060)【Abstract】Ultrasound elastography has become a new technique in the medical ultrasound field during thel
5、ast decades It stimulates human tissues vibration through different methods，such as quasi static compress，lowfrequency vibration and acoustic radiation force，and detects the tissue strain or shear wave propagation by ultra-sound to final obtain biomechanical properties of tissues Ultrasound elastogr
6、aphy has been widely applied in the di-agnosis of breast，liver and prostate，demonstrating its great potential in clinical applications Furthermore，as oneof the important components of ultrasound imaging systems，ultrasound image processing and analysis，including 3Drendering，pre processing，segmentatio
7、n，registration，image understanding and so on，is designed to help the di-agnosis of clinicians as an effective auxiliary method【Key Words】ultrasound elastography，viscoelasticity，ultrasound image processing and analysis医学超声技术的发展使得超声成像在临床诊断领域发挥着巨大的作用。由于超声成像具有安全、无创、便携、易用、价格便宜等优势，20 世纪末超声检查已占据各类医学影像检查方式的
8、四分之一。然而，许多新的医学超声技术仍然在不断涌现，并且从诊断领域跨入了治疗领域。本文将从超声弹性成像技术、超声图像处理与分析技术方面进行简述和探讨。1超声弹性成像进展28生物医学工程学进展 2013 年第 34 卷第 2 期综合述评人体组织病变与其硬度的改变有直接的关系，因此人体组织的硬度信息很早就被临床用来帮助诊断，医生用手指按压可触及的病变区域，感知肿块的大小、形状、可活动性和硬度，以此来诊断肿瘤的良恶性，其中硬度是区分良恶性最具特异性的指标。这就是广泛用于临床的检查手法，称为“触诊”(pal-pation)。“触诊”高度依赖于医生的经验，较为主观，并且无法检测难以触及的深在病灶，也很难
9、发现早期的微小病变。目前认为，粘弹性(viscoelasticity)受病理生理过程影响最大的人体组织的生物力学参数。Fung 等人对物体应力 应变(stress strain)关系进行了大量的研究，发现粘弹性作为人体软组织的基本属性是反映人体组织力学特性的最佳指标 1。目前临床所用的超声、CT、MI 等影像学手段虽然具有较高的空间分辨率，但不能提供组织弹性改变的信息 2。人体不同组织弹性模量的差异可达4 个数量级，如果能用影像手段对人体组织粘弹性进行成像或测量，则有望获得直接反映组织力学特性的高对比度分辨率的图像，从而为临床诊断提供强有力的工具。超声弹性成像是建立在组织弹性力学理论、声波在组
10、织内的非线性传播理论、信号调制与处理等基础上的系统问题，以下将针对其关键技术的研究进展与现状进行综述。1 1现有超声弹性成像方法超声弹性成像的核心是用一定的方法激励组织使其动起来，并用超声影像技术来检测组织对于激励的响应，从中可提取组织的粘弹性系数。基于不同的激励方法，有相应的组织弹性分布的重建算法。因此，在组织弹性成像中，组织激励方法是极为关键的问题。以下根据不同的激励方法，对各种超声弹性成像方法进行介绍。1 1 1组织的自然激励最早的研究利用了组织的自然激励，比如利用心脏运动引起附近肝脏组织的运动来观察肝的癌变情况，心脏运动引起的胎儿肺部组织运动来评估肺部的发育成熟度等，更为常见的是利用血
11、管的收缩和舒张来测量血管壁在一定压力下的应变情况或者弹性模量，并在此基础上发展了实时血管内超声弹性成像技术(Intravascular ultrasound elastography)3 6。1 1 2准静态激励方法研究美国Texas 大学Ophir 教授领导的研究组于1991年首次报道了利用准静态施压(quasi static com-pression)的方法获得的组织内部应变分布图，并首次提出超声弹性成像(Ultrasound Elastography)的概念 7。这篇文献影响深广，突破了超声只用于解剖成像和血流成像的限制，证明超声还可对组织的力学特性成像。准静态弹性成像以美国 Texas
12、 大学 Ophir 教授、Michigan 大学 ODonnell 教授 8、日本 Tsukuba 大学的 Shiina 教授 9、英国剑桥大学的 Treece 教授 10 等各自领导的研究组为代表。这种方法直接利用 B超探头同时完成施压和信号采集，原理简单，操作方便，因此在之后的 20 年里被广泛应用于乳腺、肾脏、前列腺、眼球等静态器官 11 15，并已经获得 FDA 批准用于临床乳腺癌的早期诊断 16 17。图 1肾脏的图片，从左往右分别是超声图像、弹性图像、病理切片11 Fig1From left to right，longitudinal sonogram，elastogram and
13、gross pathological specimen from an ovine kidney in vitro1 1 3低频振动激励方法研究与准静态弹性成像同步进行的，是另一种基于低频振动方法(low frequency vibration)的弹性成像技术，以美国 ochester 大学 Parker 教授18 20、日本东京理工学院(Tokyo Institute of Technology)的 Sato教授21，22 等各自领导的研究组为代表。他们的算法各有不同，但均采用了体外低频(20 1000)Hz振动源的激励方法，即将低频振动器置于体表，或者将振动器和超声探头集成在一起，以便于操
14、作。通过振动源向体内传播低频振动，用超声检测振动引起的剪切波在体内传播时的幅度、相位和速度信息，当体内有硬化组织时，波的传播模式会改变，由此可检测到组织内部的异变。以上两类方法都采用了接触式的机械激励(me-chanical stimulation)方法，作用力施加于人体的体表，以激励皮下组织，直接模拟了医生用手指以一定频率按压人体组织实施检查的过程，方法较为直观简单，但是由于两种方法都是对组织整体施加激励，因此在38生物医学工程学进展 2013 年第 34 卷第 2 期综合述评后续的重构组织弹性图像的过程中遇到了很多难题，也很难对深层的组织或腔内器官进行激励。1 1 4声辐射力激励方法研究为
15、了解决上述难题，以日本东京理工学院的Sugimoto 教授23、美国 Duke 大学的 Nightingale 教授24 25、俄 罗 斯 Moscow 大 学 的 Sarvazyan 教授26 27、法国巴黎大学的 Fink 教授28 29 和美国Mayo Clinic 的 Greenleaf30 32 教授为代表的研究组提出了最新的方法，采用声辐射力(acoustic radiationforce)的方法来局部激励人体内部组织，并针对声辐射力引发的局部组织振动引起的体内剪切波进行检测，剪切波的振幅、传播速度、衰减距离等与组织的力学特性直接相关，因此通过测定剪切波的声学参数可推导组织的力学特
16、性。这种方法借助了声波在组织内传播时与组织相互作用引起的声辐射力来对组织施加激励，通常采用聚焦超声在焦点处引起组织振动，并向外传播剪切波。现有的基于声辐射力的超声弹性成像方法有多种，包括:声辐射力脉冲成像(AFI)25，超音剪切波弹性成像(SSI)33，剪切波频散振动成像(SDUV)34。1 2超声弹性成像方法与现有商用超声系统的集成最早取得进展的是准静态弹性成像方法。日本的 Hitachi 公司与日本 Tsukuba 大学经过多年的合作，于 2005 年推出了世界上第一台实时组织弹性成像超声诊断仪(Hitachi EUB 8500)，并且继超声各种成像模式包括 A(amplitude)超、B
17、(Brightness)超、D(Doppler)超、M(motion)超之后，将弹性成像模式命名为 E(Elastography)超。另一个成功的例子是声辐射力弹性成像方法的商业化。西门子公司在2008 年推出 ACUSON S2000 最新高档彩超，包含了“虚拟触诊”(Virtual Touch)功能，其原理就是利用AFI 方法，通过调控激励脉冲和追踪脉冲的间隔来实现对组织的声辐射力远程激励和信号追踪。法国Echosens 公司根据瞬态弹性成像技术，研发了 Fi-broScan 产品，可提供组织的弹性模量，但该系统与现有彩超成像模式不能兼容，而且该设备只能专用于肝脏测量，其得到的结果是整个肝
18、脏的平均弹性程度，还不能精细到肝脏内部每一点的弹性系数。但鉴于其重要的临床意义，该产品在中国的市场较好，目前已经在多家三甲医院用于肝硬化的检测。1 3超声弹性成像小结综上所述，超声对人体组织硬度测量是一种新的成像方法，在肿瘤早期诊断及组织定征甚至材料力学特性检测中均有广阔的应用前景。在研究者提出的各种超声方法中，准静态弹性成像方法是最早开始研究，也是理论和算法最成熟的一种，但这种方法本身存在固有的缺点，使得它的推广受到一定的限制。声辐射力弹性成像方法是目前研究的热点，也是最有应用前景的一种方法。国外有多个团队正在分别进行这方面的研究，国内陈思平教授带领的团队也开展了利用 SDUV 方法进行大鼠
19、肝脏粘弹性测量的初步研究 35，希望能早日实现临床上的应用。图 2基于声辐射力的超声弹性成像系统框图Fig 2Diagram of ultrasound elastography systemusing acoustic radiation force2超声图像处理与分析过去 50 年间，超声医学图像处理及分析技术得以迅猛发展。作为非侵入式成像设备的代表，超声成像具有安全、无创、实时、廉价、便携等优势，早在上世纪末，超声检查占各类医学影像检查方式的比例已超过四分之一36，相应的超声图像可视化及分析技术也取得了长足进展。起初，为方便后续的可视化及分析，图像可视化及预处理技术是研究的热点。自 19
20、70 年开始，随着计算机视觉和模式识别技术的发展和逐渐成熟，研究者们把注意力转移到了超声医学图像分割、配准和理解等任务上。本节我们将分别对与超声医学图像相关的可视化、预处理、分割、配准和理解等技术方法进行综述。2 1超声图像三维可视化超声图像可视化最具代表性的成果是在产科中有广泛应用的三维超声成像，其解决了人体组织表面形态及解剖结构的三维重建，为医生提供了一种直观地诊断手段。三维超声成像技术由三部分组成:获取、重建和显示。48生物医学工程学进展 2013 年第 34 卷第 2 期综合述评1)图像的获取是精确三维超声成像的先决条件，其有三种方式，分别是机械驱动扫查、自由臂扫查和二维阵列探头扫查。
21、2)三维图像可由两种重建方式获得 37 38:a)基于特征的重建，即通过对每一幅二维图像分割出感兴趣的特征(如脏器边界)，重建出目标的表面信息。b)基于体素的重建，即通过将二维平面图像中的每一个像素都转换到一个三维坐标系中，来构建基于体素的立体图，此方法是目前最常用的重建方法。3)三维图像的显示实际上是一个体数据的可视化问题，它除了能显示组织的立体形态和结构，还可以显示该组织的任意剖面。在实际应用中，大多数图像都采用体绘制技术产生39，该技术最大的优点是可以显示实质性脏器的内部结构。直接体绘制算法大体可以分为三个大类40:以对象空间(ObjectSpace)为序的体绘制算法、以图像空间(Ima
22、geSpace)为序的体绘制算法以及对象和物体空间混合体绘制算法。这三类算法都是属于纯软件的算法，随着普通显卡里的三维加速功能的逐渐强大，基于图形硬件的体绘制算法正逐渐成为主流。现在最活跃的研究领域是彻底在 GPU 上实现光线投射算法41，从而实现基于硬件的体绘制算法可以达到和基于软件的算法同样的绘制效果。2 2超声图像预处理超声医学图像预处理主要包括图像滤波和增强。由于超声图像中的斑点噪声大大降低了图像的对比度和组织内可获取的细节信息，因此，关于超声图像的滤波技术一直是一个热门话题。这里主要介绍超声图像滤波算法。传统滤波方法如均值滤波、中值滤波、Lee 滤波、同质滤波等，实现简单，滤波效果也
23、有限。Donoho42 提出的通过修改变换小波域中的小波系数来滤除噪声的方法和 Yu 和 Acton43 提出的各向异性扩散方程(SAD)，在平滑噪声和保留图像细节上均有良好性能，从而成为这一时期的研究热点。Comaniciu 和 Meer44 提出的均值移位滤波(mean shift)算法也能在滤除噪声的同时保持边界。上述滤波技术都将超声图像中的斑点看作噪声来滤除，实际上，斑点包含了关于成像组织中散射体密度这样很有价值的信息。飞利浦医疗推出了一套实时自适应灰度图像滤波方法 XES45，其基于多分辨率算法，通过对局部梯度的估计来选择是做平滑处理还是增强处理，能提高超声图像的质量。事实上，没有一
24、个滤波算法能保证在所有应用中都是最优的。对于特定图像，以定性和定量的方式比较各种滤波算法的性能以择优选择，是一种较简单可行的方案。2 3超声图像分割图像分割一直都是超声医学图像分析领域研究的热点46，其目标是按照某种规则将图像中的像素划归为不同的类别。早期的技术如阈值化、区域生长和聚类，利用的是图像的灰度信息，然而这类方法难以分开拥有类似灰度值的两类物体;基于边缘检测的方法利用图像的梯度信息，然而这类方法对噪声敏感，且不能正确提取边界模糊的物体。由于超声图像质量较差，对于复杂的待分割目标，通常需要利用先验信息才有能取得较好的分割效果。形状先验信息随着活动轮廓模型47 的提出而引起了广泛关注，2
25、000 年以后提出的标准割48、均值移位44 和图割49 等算法则考虑到了如区域内的连续性和区域间的差异性等信息。基于神经网络的方法50 则更加专注于特定的分割目标，可融入多种图像特征。2005 年出现了一种全新的分割模式 数据库引导的分割51，其通过学习一个判别分类器来区分目标和背景，从而利用隐含在专家标注数据库中的信息实现分割。上述分割算法各有其优缺点，通常说来，对于特定的分割任务，如何在某种分割算法中融入特定目标的图像底层信息如:轮廓、纹理和某些统计特征参数，以及先验信息如:解剖结构形状先验、平滑先验和成像物理等先验才是一个分割算法成功的关键。2 4超声图像配准图像配准是医学图像处理中的
26、基本任务之一52，是指将不同时间或者不同成像设备获取的图像或体数据进行匹配，并求出图像或体数据间的空间变换关系。按照图像来源来分，医学图像配准可分为单模态和多模态配准两种;按照空间变换的种类来分，医学图像配准主要可以分为刚体和非刚体配准两种，刚体配准通常作为非刚体配准的初始条件。超声图像之间的配准主要用于超声宽景成像53 54(如图 3 所示)、空间复合成像55 以及自由臂三维超声成像的自动校准56 等。近年来，由于超声成像具有无辐射、实时及廉价等优点，超声成像被广泛应用于手术引导，但是超声图像固有的斑点噪声、分辨率低及对比度差等问题影响了手术引导的58生物医学工程学进展 2013 年第 34
27、 卷第 2 期综合述评效果，因此术中超声成像与术前其它模态影像的结合正成为当前的研究热点，超声图像与其它成像模式间的多模态配准正是其中的关键技术。由于成像原理的不同，超声图像与其它模态图像(如 CT，MI)之间差异较大，因此配准是非常困难的，通常需要基于定位器来提供自动配准的初始条件。目前常用的超声与其它模态图像间的配准方法主要可以分为两种:(1)基于灰度的方法，如 Correlation atio 57 58，CrossCorrelation59，互信息60 等;(2)基于表面(Surface)匹配的方法61 62 等。超声图像与其它模态图像配准的关键是根据特定的临床应用，寻找两者之间的有效
28、相似测度，如Wein63 提出通过将 CT 模拟成超声图像再与超声图像进行配准的方法可有效提高配准的鲁棒性。虽然现有方法已经取得了一定进展，但是超声图像与其它模态影像间的刚体和非刚体配准仍将是未来研究的重点和难点。图 3基于 3D SIFT 特征匹配的三维超声宽景成像53 Fig 3Volumetric ultrasound panorama based on 3D SIFT2 5超声图像理解机器学习在医学图像理解任务中起到了至关重要的作用，其成功地促成了大量的应用，包括计算机辅助诊断、病案检索、组织器官分割、定量分析等64。对于那些缺乏精确知识和难以分析推导的问题来说，从实例中自动学习的算法
29、就能发挥其重要作用。近年来，由于许多高性价比的超声成像设备正越来越广地应用于临床，从而积累了大量的医学图像数据库，这使得机器学习在医学中的应用变得更加可行。用于超声医学图像理解的机器学习算法中，核学习和概率模型扮演着关键角色:核学习(支持向量机)能给超声医学计算机辅助检测及诊断任务提供最佳分类器;概率模型(朴素贝叶斯、图模型、马尔可夫随机场)为医学图像分析任务如图像重建、分割、配准等提供理论框架。对于学习分类器来说，线性模型(线性判别分析)、人工神经网络(误差反向传播、自组织映射)和集成学习(Boos-ting、AdaBoost)是除核学习之外的选择。降维和特征选择是计算机辅助诊断的关键技术。
30、聚类分析可用于识别超声医学图像中的相似病灶。国内汪天富教授的团队较早开展了超声图像理解方面的研究工作，主要基于人工神经网络方法，分别研究了基于组织多普勒图像的窦房结定量分析 65、基于 B 超图像的脂肪肝与正常肝自动识别 66 和基于彩色多普勒血流图像和 B 超图像相结合的乳腺肿瘤计算机辅助诊断 67 等;最近，陈思平教授的团队提出二维超声“智能扫查”理念，并基于集成学习方法，实现了二维超声影像中早孕囊标准切面的自动提取和径线的自动测量，如图 4 所示 68，从而有望提高产前超声检测的标准化和规范化程度。图 4二维超声“智能扫查”示意图Fig4Illustration of 2D intell
31、igent ultrasound scanning method医生手持超声探头扫过孕妇腹部，产生一系列切面图像(a、b、c)，计算机则自动从中找出标准切面图像(即 b，该切面上所含的孕囊最大)68。在模式分类领域，我们没有理由偏好任何一个分类器69。因此，实际应用时，通常需要对比多种分类器的分类性能来选取一个对该应用、该组数据、该组特征较优的分类器。Forsyth 和 Ponce70 指出，可先尝试核学习，然后线性模型和神经网络等。超声图像处理与分析涉及内容广泛，远不止以上综述的五个部分，随着计算机科学和超声成像设备的发展，超声图像处理和分析也将向智能化、鲁棒性强和实时快速等方向发展。参考文
32、献 1Fung Y C Biomechanics:Mechanical properties of livingtissues M 2nd Edition ed:Springer 1994 2Greenleaf J F，Fatemi M，Insana M Selected methods forimaging elastic properties of biological tissues J Annu68生物医学工程学进展 2013 年第 34 卷第 2 期综合述评ev Biomed Eng，2003，5:57 78 3de Korte C L，Sierevogel M J，Mastik F
33、，et al Identifica-tion of atherosclerotic plaque components with intravascularultrasound elastography in vivo a Yucatan pig studyJ Circulation，2002，105(14):1627 1630 4Lee W N，Ingrassia C M，Fung Kee Fung S D，et alTheoretical quality assessment of myocardial elastographywith in vivo validationJ IEEE T
34、rans Ultrason Ferro-electr Freq Control，2007，54(11):2233 2245 5Schaar J A，de Korte C L，Mastik F，et al Characterizingvulnerable plaque features with intravascular elastography J Circulation，2003，108(21):2636 2641 6Zhang P，Zhang Y The technical principle and clinical ap-plication of intravascular ultr
35、asound elastogramJ J Bi-omed Eng，2005，22(1):163 7Ophir J，Cespedes I，Ponnekanti H，et al Elastography:aquantitative method for imaging the elasticity of biologicaltissues J Ultrasonic imaging，1991，13(2):111 134 8 Lubinski M A，Emelianov S Y，ODonnell M Speckletracking methods for ultrasonic elasticity i
36、maging usingshort time correlationJ IEEE Trans Ultrason Ferro-electr Freq Control，1999，46(1):82 96 9Yamakawa M，Shiina T Strain estimation using the extend-ed combined autocorrelation method J Jpn J Appl Phys，2001，40:3872 10Treece G M，Lindop J E，Gee A H，et al Freehand ul-trasound elastography with a
37、3 D probeJ UltrasoundMed Biol，2008，34(3):463 474 11Ophir J，Alam S K，Garra B S，et al Elastography:ima-ging the elastic properties of soft tissues with ultrasound J J Med，2002，29:156 12罗建文，白净 超声弹性成像的研究进展 J 中国医疗器械信息，2005，11(005):23 31 13Itoh A，Ueno E，Tohno E，et al Breast Disease:ClinicalApplication of
38、US Elastography for Diagnosis J adiol-ogy，2006，239(2):341 350 14关少卿，罗葆明，欧冰等 超声弹性成像诊断乳腺局灶性病变价值的初步探讨J 临床超声医学杂志，2006，8(10):598 601 15王怡，王涌，张希敏等 实时组织弹性成像技术在鉴别诊断乳腺良恶性肿块中的价值评估J 中华超声影像学杂志，2005，14(12):911 913 16Ginat D T，Destounis S V，Barr G，et al US Elastog-raphy of Breast and Prostate LesionsJ adiographic
39、s，2009，29(7):2007 2016 17Salomon G，K?llerman J，Thederan I，et al Evaluationof prostate cancer detection with ultrasound real timeelastography:a comparison with step section pathologicalanalysis after radical prostatectomy J Eur Urol，2008，54(6):1354 1362 18 Huang S，Lerner M，Parker K J Time domainDoppl
40、er estimators of the amplitude of vibrating targets J J Acoust Soc Am，1992，91:965 19ubens D J，Hadley M A，Alam S K，et al Sonoelastici-ty imaging of prostate cancer:in vitro results J adiol-ogy，1995，195(2):379 383 20Parker K J，Taylor L S，Gracewski S，et al A unifiedview of imaging the elastic propertie
41、s of tissueJ JAcoust Soc Am，2005，117:2705 21Yamakoshi Y，Sato J，Sato T Ultrasonic imaging of inter-nal vibration of soft tissue under forced vibrationJ IEEE Trans Ultrason Ferroelectr Freq Control，1990，37(2):45 53 22Yamakoshi Y，Otaki H，Shinozuka N，et al Internal tis-sue displacement measurement based
42、 on ultrasonic waveDoppler signal digital detection and its application to fetalmovement monitoringJ Ultrasonics，1996，34(8):769 775 23Sugimoto T，Ueha S，Itoh K Tissue hardness measure-ment using the radiation force of focused ultrasoundC IEEE Ultrasonics SymposiumHonolulu，HI:IEEE，1990 24Nightingale K
43、，Kornguth P J，Trahey G E The use ofacoustic streaming in breast lesion diagnosis:a clinicalstudy J Ultrasound Med Biol，1999，25(1):75 87 25Nightingale K，Soo M S，Nightingale，et al Acousticradiation force impulse imaging:in vivo demonstration ofclinical feasibilityJ Ultrasound Med Biol，2002，28(2):227 2
44、35 26udenko O，Sarvazyan A，Emelianov S Y Acoustic radi-ation force and streaming induced by focused nonlinear ul-trasound in a dissipative mediumJ J Acoust Soc Am，1996，99:2791 27Sarvazyan A P，udenko O V，Swanson S D，et al Shearwave elasticity imaging:a new ultrasonic technology ofmedical diagnostics J
45、 Ultrasound Med Biol，1998，24(9):1419 1435 28Muller M，Gennisson J L，Deffieux T，et al Quantitativeviscoelasticity mapping of human liver using supersonicshear imaging:preliminary in vivo feasibility study J Ultrasound Med Biol，2009，35(2):219 229 29Bavu E，Gennisson J，Couade M，et al Noninvasive in vi-vo
46、 liver fibrosis evaluation using supersonic shear imaging:a clinical study on 113 hepatitis C virus patients J Ul-trasound Med Biol，2011，37(9):1361 1373 30Fatemi M，Greenleaf J F Ultrasound stimulated vibro acoustic spectrography J Science，1998，280(5360):82 85 31Fatemi M，Greenleaf J F Vibro acoustogr
47、aphy:An ima-ging modality based on ultrasound stimulated acoustic e-missionJ Proc Natl Acad Sci U S A，1999，96(12):6603 6608 32Fatemi M，Wold L E，Alizad A，et al Vibro acoustictissue mammographyJ IEEE Trans Med Imaging，2002，21(1):1 8 33Bercoff J，Tanter M，Fink M Supersonic shear imaging:a new technique
48、for soft tissue elasticity mappingJ IEEE Trans Ultrason Ferroelectr Freq Control，2004，51(4):396 409 34Chen S，Urban M，Pislaru C，et al Shearwave dispersionultrasound vibrometry(SDUV)for measuring tissue elas-ticity and viscosityJ IEEE Trans Ultrason FerroelectrFreq Control，2009，56(1):55 62 35Wang Y，Ch
49、en S，Wang T，et al Development of a ge-neric ultrasound vibro acoustic imaging platform for tis-sue elasticity and viscosityJ J Innovative Opt Health78生物医学工程学进展 2013 年第 34 卷第 2 期综合述评Sci，2012，5(1):1 7 36陈思平 医学超声影像产业现状和发展J 应用声学，2005(4):201 207 37Fenster A，Downey D B Three dimensional ultrasoundimagingJ
50、 Annu ev Biomed Eng，2000，2:457 475 38Solberg O V，Lindseth F，Torp H，et al Freehand 3D ul-trasound reconstruction algorithms a reviewJ Ul-trasound Med Biol，2007，33(7):991 1009 39汪天富，郑昌琼，李德玉等 超声心脏图像多维重建研究进展J 生物医学工程学杂志，2001，18(1):133 137 40Preim B，Bartz D Visualization in medicine:theory，algo-rithms，and