CN109674490A - 一种超声引导的低反射伪像光声显微镜成像方法 - Google Patents
一种超声引导的低反射伪像光声显微镜成像方法 Download PDFInfo
- Publication number
- CN109674490A CN109674490A CN201910043220.5A CN201910043220A CN109674490A CN 109674490 A CN109674490 A CN 109674490A CN 201910043220 A CN201910043220 A CN 201910043220A CN 109674490 A CN109674490 A CN 109674490A
- Authority
- CN
- China
- Prior art keywords
- signal
- photoacoustic
- ultrasonic
- ultrasonic transducer
- sample
- 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
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/4416—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device related to combined acquisition of different diagnostic modalities, e.g. combination of ultrasound and X-ray acquisitions
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0033—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
- A61B5/0035—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for acquisition of images from more than one imaging mode, e.g. combining MRI and optical tomography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0033—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
- A61B5/004—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for image acquisition of a particular organ or body part
- A61B5/0042—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for image acquisition of a particular organ or body part for the brain
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0093—Detecting, measuring or recording by applying one single type of energy and measuring its conversion into another type of energy
- A61B5/0095—Detecting, measuring or recording by applying one single type of energy and measuring its conversion into another type of energy by applying light and detecting acoustic waves, i.e. photoacoustic measurements
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Animal Behavior & Ethology (AREA)
- Radiology & Medical Imaging (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Pathology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Acoustics & Sound (AREA)
- Neurology (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本发明公开了一种超声引导的低反射伪像光声显微镜成像方法,利用超声波作引导,实现低伪像的穿颅光声显微镜成像方法。在做小动物脑部光声成像时,光声信号在穿过头盖骨时会产生多次的反射,从而产生严重的伪像。本发明利用光声‑超声双模显微镜***,采集每个位置的光声信号和超声信号。利用超声反射信号和光声信号穿过颅骨时类似的传递特性,进而利用超声反射信号估计光声传递函数,并通过反卷积算法,去除光声图像中的反射伪像。本发明提出的成像方法,可以透过颅骨非侵入地获得高质量的三维光声图像。
Description
技术领域
本发明涉及一种利用超声回波作引导,实现低伪像的穿颅光声显微镜成像方法。具体来说是在一种在结合超声显微镜技术和光声显微镜技术的双模显微镜的基础上,利用超声回波信号与光声信号具有类似的穿颅传递函数,采用超声回波估计光声传递函数,并修正光声成像结果,从而抑制光声穿颅成像的反射伪像,改善成像质量。
背景技术
光声成像是一种快速发展的生物医学成像技术,它既具有丰富的光学吸收对比度,又具有较大的成像深度和良好的空间分辨率。声学分辨率光声显微镜是光声成像的一种重要的实现方式。在脉冲激光的照射下,由于光声效应,组织内部的光学吸收体吸收电磁能而产生宽带超声波。光声信号由组织外的点聚焦超声换能器接收。根据接收到的光声信号的幅度和时延,可以重建沿换能器声轴的一维光学吸收信息。通过二维平面扫描,可以获得三维图像。受益于生物组织的低声学散射,这种成像方式可以打破光学散射的限制,提供深层组织声学分辨率的光吸收图像。这些优点使得声学分辨率的光声显微镜在生物医学成像方面具有很大潜力,例如脑监测、癌症检测、微血管可视化、细胞成像等。
脑成像一直是医学成像较为关注的问题,光声显微镜可以获得一定深度上足够分辨率的小动物的脑部三维图像。但是头盖骨会导致光声信号的多次反射,从而严重影响成像的质量。由于头盖骨较高的声阻抗,光声信号在穿过头盖骨时会在骨中多次反射,从而在真实的光声信号后面会有很多延迟发射信号,这在光声图像上会表现为伪像。这种伪像会和真实的图像混叠,导致真实的信息难以被辨别。这一问题会严重限制此技术的活体脑成像的效果。因此开发一种在非侵入成像情况下,抑制反射伪像,改善光声成像质量的方法就显得意义重大。
发明内容
发明目的:针对现有光声显微镜透颅成像存在反射伪像的问题,本发明提供一种利用超声回波作引导消除反射伪像的新型的光声显微镜成像方法,显著地提升光声成像质量;该方法结合了光声显微镜和超声显微镜两种技术,能够有效地去除光声信号中存在的反射伪像,为准确识别小动物脑部信息提供有利的依据。
技术方案:为实现上述目的,本发明采用的技术方案为:
一种超声引导的低反射伪像光声显微镜成像方法,采用一个三通同时连通超声换能器、信号发生器和数据采集卡,其中超声换能器为球聚焦水浸式超声换能器;该方法为:脉冲激光器发出的脉冲激光通过耦合光纤照射在样品上,样品产生的光声信号经放大后被超声换能器接收;脉冲激光同时也被光电二极管接收,光电二极管产生的光电信号触发数据采集卡进行采集信号,采集的信号延迟一段时间后触发信号发生器;信号发生器每触发一次就发出一个窄脉冲电信号,窄脉冲信号经超声换能器转换为超声脉冲发射至样品上,超声换能器同时接收来自样品的超声回波信号;超声回波信号经放大后被数据采集卡采集;该过程称为双模态成像过程;双模态成像过程同步获取了样品每个位置点的光声信号和超声回波信号,根据信号的幅度和时延,即可重建沿着超声换能器声轴的光吸收和声反射的一维信息;最后通过电机的逐点扫描便可得到样品的光声-超声三维图像。
该方法的具体计算过程包括:
(1)估算头盖骨上(x,y)位置处的声传递函数
记(x,y)位置点产生的超声回波信号为u(t),超声换能器发出的超声脉冲为u0(t),能够得到(x,y)位置处的冲激响应函数hus(t):
hus(t)=F-1[F(u)/F(u0)] (1)
其中:F[·]和F-1[·]表示对·进行傅里叶变换和傅里叶逆变换,t表示时间;
(2)消除光声信号中的反射波成分
设(x,y)位置处接收到的光声信号为p(t)、脉冲激光器发出的脉冲激光激发光吸收体而产生的光声信号为p0(t)、光声信号的传递函数为hap(t),则有p(t)=p0(t)*hap(t),其中*代表卷积;考虑到超声回波信号有数量关系光声信号有数量关系其中:R和T分别表示头骨的声反射和声透射系数,τ表示信号在头骨中反射导致的时间延迟,n表示反射次数,获得hap(t)≈hus(t);
因此根据超声回波信号与光声信号的传递函数,并利用hap(t)≈hus(t),同时对p(t)做反卷积,以消除p(t)中的反射波成分:
P0(t)=F-1[F(P)/F(hpa)]≈F-1[F(P)/F(hus)] (2)
其中:P(t)是p(t)的希尔伯特变换;
(3)计算(x,y)位置处的去除伪像的三维图像A(x,y;z)=||P0(z)||,z=tz/c,c为头骨中的声速,tz为光声信号的传播时间,零点时刻为脉冲激光激发时刻;
(4)对于扫描平面上的每个位置(x,y),重复步骤(1)~(3),得到去除伪像的三维图像A(x,y,z)。
利用本发明的成像***,不仅可以在一次扫描的情况下同时获得光声和超声两幅三维图像,而且采集到的超声回波信号还可以用来修正光声信号,去除光声图像中因头盖骨等引起的反射伪像。利用这样的方法,我们可以在不破坏头盖骨的情况下获得高质量的光声成像结果。
有益效果:本发明提供的超声回波引导的光声显微镜成像方法,利用超声回波信号引导光声图像重构,可以去除颅骨等引起的反射伪像,有效提高光声成像质量;本发明还可以在无损、非侵入的情况下获得较高质量的小动物脑部的光声图像;另外,本发明的***还同时提供了三维超声图像和光声图像。
附图说明
图1为本发明使用的光声-超声双模显微镜***示意图;
图2为颅骨等产生反射伪像的原理、以及光声信号和超声回波信号激发接收的示意图;
图3为超声显微镜、声学分辨率光声显微镜、以及本方法的获得的一维信号和B平面(x-y平面)图像对比;
图4为传统光声显微镜和本方法获得的三层不同深度的C平面(x-y平面)图像对比。
具体实施方式
下面结合附图对本发明作更进一步的说明。
一种超声引导的光声显微镜成像方法,采用一个三通同时连通超声换能器、信号发生器和数据采集卡,其中超声换能器为球聚焦水浸式超声换能器;该方法为:脉冲激光器发出的脉冲激光通过耦合光纤照射在样品上,样品产生的光声信号经放大后被超声换能器接收;脉冲激光同时也被光电二极管接收,光电二极管产生的光电信号触发数据采集卡进行采集信号,采集的信号延迟一段时间后触发信号发生器;信号发生器每触发一次就发出一个窄脉冲电信号,窄脉冲信号经超声换能器转换为超声脉冲发射至样品上,超声换能器同时接收来自样品的超声回波信号;超声回波信号经放大后被数据采集卡采集;该过程称为双模态成像过程;双模态成像过程同步获取了样品每个位置点的光声信号和超声回波信号,根据信号的幅度和时延,即可重建沿着超声换能器声轴的光吸收和声反射的一维信息;最后通过电机的逐点扫描便可得到样品的光声-超声三维图像。该方法的具体计算过程包括:
(1)估算头盖骨上(x,y)位置点的传递函数
记(x,y)位置点产生的超声回波信号为u(t),超声换能器发出的超声脉冲为u0(t),能够得到(x,y)位置处的冲激响应函数hus(t):
hus(t)=F-1[F(u)/F(u0)] (1)
其中:F[·]和F-1[·]表示对·进行傅里叶变换和傅里叶逆变换,t表示时间;
(2)消除光声信号中的反射波成分
设(x,y)位置处接收到的光声信号为p(t)、脉冲激光器发出的脉冲激光为p0(t)、光声信号的传递函数为hap(t),则有p(t)=p0(t)*hap(t),其中*代表卷积;考虑到超声回波信号有数量关系光声信号有数量关系其中:R和T分别表示头骨的声反射和声透射系数,τ表示信号在头骨中反射导致的时间延迟,n表示反射次数,获得hap(t)≈hus(t);
因此根据超声回波信号与光声信号的传递函数,并利用hap(t)≈hus(t),同时对p(t)做反卷积,以消除p(t)中的反射波成分:
P0(t)=F-1[F(P)/F(hpa)]≈F-1[F(P)/F(hus)] (2)
其中:P(t)是p(t)的希尔伯特变换;
(3)计算(x,y)位置处的去除伪像的三维图像A(x,y;z)=||P0(z)||,z=tz/c,c为头骨中的声速,tz为光声信号的传播时间,零点时刻为脉冲激光激发时刻;
(4)对于扫描平面上的每个位置(x,y),重复步骤(1)~(3),得到去除伪像的三维图像A(x,y,z)。
图3和图4给出了一种基于本发明的具体实施例。选用若干钨丝(模拟脑部血管)随机排布,并用一张亚克力薄膜(模拟骨骼)压在钨丝之上。用本发明的***对该样品进行扫描,得到其三维图像。并且利用我们的方法进行处理,从图3可明显看出在一维,二维图像上伪像都得到较好去除。为了更好地体现此方法的效果,图4展示了该三维图像的三层C平面图像以及处理之后的结果。
以上所述仅是本发明的优选实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (2)
1.一种超声引导的低反射伪像光声显微镜成像方法,其特征在于:采用一个三通同时连通超声换能器、信号发生器和数据采集卡,其中超声换能器为球聚焦水浸式超声换能器;该方法为:脉冲激光器发出的脉冲激光通过耦合光纤照射在样品上,样品产生的光声信号经放大后被超声换能器接收;脉冲激光同时也被光电二极管接收,光电二极管产生的光电信号触发数据采集卡进行采集信号,采集的信号延迟一段时间后触发信号发生器;信号发生器每触发一次就发出一个窄脉冲电信号,窄脉冲信号经超声换能器转换为超声脉冲发射至样品上,超声换能器同时接收来自样品的超声回波信号;超声回波信号经放大后被数据采集卡采集;该过程称为双模态成像过程;双模态成像过程同步获取了样品每个位置点的光声信号和超声回波信号,根据信号的幅度和时延,即可重建沿着超声换能器声轴的光吸收和声反射的一维信息;最后通过电机的逐点扫描便可得到样品的光声-超声三维图像。
2.根据权利要求1所述的超声引导的低反射伪像光声显微镜成像方法,其特征在于:该方法的具体计算过程包括:
(1)估算头盖骨上(x,y)位置处的声传递函数
记(x,y)位置点产生的超声回波信号为u(t),超声换能器发出的超声脉冲为u0(t),能够得到(x,y)位置处的冲激响应函数hus(t):
hus(t)=F-1[F(u)/F(u0)] (1)
其中:F[·]和F-1[·]表示对·进行傅里叶变换和傅里叶逆变换,t表示时间;
(2)消除光声信号中的反射波成分
设(x,y)位置处接收到的光声信号为p(t)、脉冲激光器发出的脉冲激光激发光吸收体而产生的光声信号为p0(t)、光声信号的传递函数为hap(t),则有p(t)=p0(t)*hap(t),其中*代表卷积;考虑到超声回波信号有数量关系光声信号有数量关系其中:R和T分别表示头骨的声反射和声透射系数,τ表示信号在头骨中反射导致的时间延迟,n表示反射次数,获得hap(t)≈hus(t);
因此根据超声回波信号与光声信号的传递函数,并利用hap(t)≈hus(t),同时对p(t)做反卷积,以消除p(t)中的反射波成分:
P0(t)=F-1[F(P)/F(hpa)]≈F-1[F(P)/F(hus)] (2)
其中:P(t)是p(t)的希尔伯特变换;
(3)计算(x,y)位置处的去除伪像的三维图像A(x,y;z)=||P0(z)||,z=tz/c,c为头骨中的声速,tz为光声信号的传播时间,零点时刻为脉冲激光激发时刻;
(4)对于扫描平面上的每个位置(x,y),重复步骤(1)~(3),得到去除伪像的三维图像A(x,y,z)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910043220.5A CN109674490B (zh) | 2019-01-17 | 2019-01-17 | 一种超声引导的低反射伪像光声显微镜成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910043220.5A CN109674490B (zh) | 2019-01-17 | 2019-01-17 | 一种超声引导的低反射伪像光声显微镜成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109674490A true CN109674490A (zh) | 2019-04-26 |
CN109674490B CN109674490B (zh) | 2021-09-10 |
Family
ID=66193326
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910043220.5A Active CN109674490B (zh) | 2019-01-17 | 2019-01-17 | 一种超声引导的低反射伪像光声显微镜成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109674490B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111398175A (zh) * | 2020-03-05 | 2020-07-10 | 南京大学 | 一种声速自适应的光声-超声双模态显微镜成像方法 |
CN111759276A (zh) * | 2020-06-11 | 2020-10-13 | 南昌大学 | 一种基于光纤延时的多通道光声显微成像方法 |
CN111938579A (zh) * | 2020-06-29 | 2020-11-17 | 南昌大学 | 一种基于光纤延时的提升光声图像质量的多脉冲光声显微成像方法及其装置 |
CN115919362A (zh) * | 2023-03-15 | 2023-04-07 | 深圳英美达医疗技术有限公司 | 超声成像***的伪像去除方法、装置、设备及存储介质 |
Citations (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060184042A1 (en) * | 2005-01-22 | 2006-08-17 | The Texas A&M University System | Method, system and apparatus for dark-field reflection-mode photoacoustic tomography |
CN103300880A (zh) * | 2012-03-13 | 2013-09-18 | 佳能株式会社 | 被检体信息获得装置和被检体信息获得方法 |
US8875582B2 (en) * | 2010-03-25 | 2014-11-04 | Canon Kabushiki Kaisha | Apparatus and method for photoacoustic imaging |
CN104490365A (zh) * | 2015-01-12 | 2015-04-08 | 南京大学 | 一种使用窄带传感器重建宽带光声图像的方法 |
WO2015095721A1 (en) * | 2013-12-20 | 2015-06-25 | Fujifilm Sonosite, Inc. | High frequency ultrasound transducers |
US20150182122A1 (en) * | 2012-07-26 | 2015-07-02 | The Institute Of Cancer Research: Royal Cancer Hospital | Ultrasonic imaging |
CN104771136A (zh) * | 2014-01-14 | 2015-07-15 | 三星麦迪森株式会社 | 光声设备和操作所述光声设备的方法 |
CN105395170A (zh) * | 2015-12-15 | 2016-03-16 | 同济大学 | 一种光声超声双模态同步成像***及方法 |
WO2016182435A1 (en) * | 2015-05-08 | 2016-11-17 | Universiteit Twente | Artifact reduction in photoacoustic and thermoacoustic imaging |
CN106901693A (zh) * | 2017-01-19 | 2017-06-30 | 南京大学 | 基于环形超声换能器阵列的大景深光声‑声速双模成像法 |
CN107157448A (zh) * | 2017-05-25 | 2017-09-15 | 睿芯生命科技(深圳)有限公司 | 用于浅表部位成像的光声与超声同步成像***和方法 |
CN107223035A (zh) * | 2017-01-23 | 2017-09-29 | 深圳迈瑞生物医疗电子股份有限公司 | 一种成像***、方法及超声成像*** |
WO2018056184A1 (ja) * | 2016-09-21 | 2018-03-29 | 富士フイルム株式会社 | 光音響画像生成装置 |
CN108024794A (zh) * | 2015-09-29 | 2018-05-11 | 富士胶片株式会社 | 光声测量装置及光声测量装置的信号处理方法 |
CN105548022B (zh) * | 2015-12-23 | 2018-08-10 | 深圳先进技术研究院 | 一种暗场照明声学分辨率光声显微镜 |
CN108888238A (zh) * | 2018-05-11 | 2018-11-27 | 南京大学深圳研究院 | 一种基于带孔超声反射板的光声显微镜及成像法 |
CN109091109A (zh) * | 2018-07-02 | 2018-12-28 | 南京大学 | 基于全矩阵滤波和时间反转算子的优化型光声断层成像的图像重构方法 |
CN109199333A (zh) * | 2018-09-28 | 2019-01-15 | 同济大学 | 基于光纤传导式的光声和超声双模态内窥探头装置及方法 |
-
2019
- 2019-01-17 CN CN201910043220.5A patent/CN109674490B/zh active Active
Patent Citations (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060184042A1 (en) * | 2005-01-22 | 2006-08-17 | The Texas A&M University System | Method, system and apparatus for dark-field reflection-mode photoacoustic tomography |
US8875582B2 (en) * | 2010-03-25 | 2014-11-04 | Canon Kabushiki Kaisha | Apparatus and method for photoacoustic imaging |
CN103300880A (zh) * | 2012-03-13 | 2013-09-18 | 佳能株式会社 | 被检体信息获得装置和被检体信息获得方法 |
US20150182122A1 (en) * | 2012-07-26 | 2015-07-02 | The Institute Of Cancer Research: Royal Cancer Hospital | Ultrasonic imaging |
WO2015095721A1 (en) * | 2013-12-20 | 2015-06-25 | Fujifilm Sonosite, Inc. | High frequency ultrasound transducers |
CN104771136A (zh) * | 2014-01-14 | 2015-07-15 | 三星麦迪森株式会社 | 光声设备和操作所述光声设备的方法 |
CN104490365A (zh) * | 2015-01-12 | 2015-04-08 | 南京大学 | 一种使用窄带传感器重建宽带光声图像的方法 |
WO2016182435A1 (en) * | 2015-05-08 | 2016-11-17 | Universiteit Twente | Artifact reduction in photoacoustic and thermoacoustic imaging |
CN108024794A (zh) * | 2015-09-29 | 2018-05-11 | 富士胶片株式会社 | 光声测量装置及光声测量装置的信号处理方法 |
CN105395170A (zh) * | 2015-12-15 | 2016-03-16 | 同济大学 | 一种光声超声双模态同步成像***及方法 |
CN105548022B (zh) * | 2015-12-23 | 2018-08-10 | 深圳先进技术研究院 | 一种暗场照明声学分辨率光声显微镜 |
WO2018056184A1 (ja) * | 2016-09-21 | 2018-03-29 | 富士フイルム株式会社 | 光音響画像生成装置 |
CN106901693A (zh) * | 2017-01-19 | 2017-06-30 | 南京大学 | 基于环形超声换能器阵列的大景深光声‑声速双模成像法 |
CN107223035A (zh) * | 2017-01-23 | 2017-09-29 | 深圳迈瑞生物医疗电子股份有限公司 | 一种成像***、方法及超声成像*** |
CN107157448A (zh) * | 2017-05-25 | 2017-09-15 | 睿芯生命科技(深圳)有限公司 | 用于浅表部位成像的光声与超声同步成像***和方法 |
CN108888238A (zh) * | 2018-05-11 | 2018-11-27 | 南京大学深圳研究院 | 一种基于带孔超声反射板的光声显微镜及成像法 |
CN109091109A (zh) * | 2018-07-02 | 2018-12-28 | 南京大学 | 基于全矩阵滤波和时间反转算子的优化型光声断层成像的图像重构方法 |
CN109199333A (zh) * | 2018-09-28 | 2019-01-15 | 同济大学 | 基于光纤传导式的光声和超声双模态内窥探头装置及方法 |
Non-Patent Citations (9)
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111398175A (zh) * | 2020-03-05 | 2020-07-10 | 南京大学 | 一种声速自适应的光声-超声双模态显微镜成像方法 |
CN111759276A (zh) * | 2020-06-11 | 2020-10-13 | 南昌大学 | 一种基于光纤延时的多通道光声显微成像方法 |
CN111938579A (zh) * | 2020-06-29 | 2020-11-17 | 南昌大学 | 一种基于光纤延时的提升光声图像质量的多脉冲光声显微成像方法及其装置 |
CN115919362A (zh) * | 2023-03-15 | 2023-04-07 | 深圳英美达医疗技术有限公司 | 超声成像***的伪像去除方法、装置、设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN109674490B (zh) | 2021-09-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109674490A (zh) | 一种超声引导的低反射伪像光声显微镜成像方法 | |
EP2328480B1 (en) | Low-cost device for c-scan photoacoustic imaging | |
US9492139B2 (en) | Non-imaging low frequency ultrasonic testing and diagnostic evaluation system | |
CN211511858U (zh) | 一种光学超声双导管内窥成像*** | |
US6926672B2 (en) | Electret acoustic transducer array for computerized ultrasound risk evaluation system | |
US7285092B2 (en) | Computerized ultrasound risk evaluation system | |
US7025725B2 (en) | Three-dimensional ultrasound computed tomography imaging system | |
CN100446730C (zh) | 基于声透镜的光声成像和层析成像方法及其装置 | |
CN102944521B (zh) | 非接触式光声和光学相干断层双成像装置及其检测方法 | |
EP2514364B1 (en) | Measurement system, and image forming method and program | |
CN102697526B (zh) | 超声浅表组织与器官容积扫描断层成像设备 | |
CN100434042C (zh) | 生物组织光学和超声的采集和层析成像的方法及其装置 | |
CN101529241B (zh) | 使用超声波换能器对生物组织粘弹特性测量的方法 | |
CN102440806A (zh) | 用于电磁噪声消除的超声方法和探头 | |
CN107174202A (zh) | 一种基于主动探测的磁声成像方法及*** | |
CA3127900A1 (en) | Method and system for ultrasonic characterization of a medium | |
Hemmsen et al. | Tissue harmonic synthetic aperture ultrasound imaging | |
RU2378989C2 (ru) | Способ диагностики с помощью ультразвуковых, звуковых и электромагнитных волн | |
CA3127822A1 (en) | Method and system for ultrasonic characterization of a medium | |
CN102688071A (zh) | 超声浅表组织与器官容积扫描断层成像方法 | |
CN108601583A (zh) | 声波图像生成装置及声波图像生成方法 | |
CN109864707A (zh) | 一种在有限视角情况下提高光声断层成像分辨率的方法 | |
Thompson et al. | Laser-induced synthetic aperture ultrasound imaging | |
CA3127873A1 (en) | Method and system for ultrasonic characterization of a medium | |
Jin et al. | Attenuation compensation for high-frequency acoustic-resolution photoacoustic imaging |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |