CN115177217A - 基于球形粒子光脉冲激发效应的光声信号仿真方法、装置 - Google Patents

基于球形粒子光脉冲激发效应的光声信号仿真方法、装置 Download PDF

Info

Publication number
CN115177217A
CN115177217A CN202211100325.8A CN202211100325A CN115177217A CN 115177217 A CN115177217 A CN 115177217A CN 202211100325 A CN202211100325 A CN 202211100325A CN 115177217 A CN115177217 A CN 115177217A
Authority
CN
China
Prior art keywords
photoacoustic
array
photoacoustic signal
spherical particle
detection array
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
Application number
CN202211100325.8A
Other languages
English (en)
Other versions
CN115177217B (zh
Inventor
孟彧仟
施钧辉
陈睿黾
李驰野
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang Lab
Original Assignee
Zhejiang Lab
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Zhejiang Lab filed Critical Zhejiang Lab
Priority to CN202211100325.8A priority Critical patent/CN115177217B/zh
Publication of CN115177217A publication Critical patent/CN115177217A/zh
Application granted granted Critical
Publication of CN115177217B publication Critical patent/CN115177217B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Surgery (AREA)
  • Medical Informatics (AREA)
  • Veterinary Medicine (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Molecular Biology (AREA)
  • Artificial Intelligence (AREA)
  • Psychiatry (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physiology (AREA)
  • Signal Processing (AREA)
  • Acoustics & Sound (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

本发明提出了一种基于球形粒子光脉冲激发效应的光声信号仿真方法、装置，所述方法包括：设定包括采样频率、采样时间和光声信号传输速度；设置光声断层成像探测阵列，并计算各阵元的位置坐标；将待测体模分解为多个球形粒子；设定阵元序号循环数；计算光声断层成像探测阵列单个阵元接收到的光声信号，将每个阵元得到的光声信号进行拼接得到光声信号矩阵。本发明方法将待检测体模分解为多个球形粒子，根据光声成像原理，直接考虑待测体模发出的所有光声信号经一定距离传播后在光声断层成像探测阵列单个阵元上的叠加效果，计算得出光声断层成像探测阵列各阵元接收到的上述粒子构成的体模受激光激发所产生的光声信号。

Description

基于球形粒子光脉冲激发效应的光声信号仿真方法、装置
技术领域
本发明涉及光声断层成像领域，尤其涉及一种基于球形粒子光脉冲激发效应的光声信号仿真方法、装置。
背景技术
光声成像（Photoacoustic Imaging, 简称PAI）是一种非入侵式和非电离式的新型生物医学成像方法。当脉冲激光照射到（热声成像则特指用无线电频率的脉冲激光进行照射）生物组织中时，组织的光吸收域将产生超声信号，我们称这种由光激发产生的超声信号为光声信号。生物组织产生的光声信号携带了组织的光吸收特征信息，通过探测光声信号能重建出组织中的光吸收分布图像。光声成像结合了纯光学组织成像中高选择特性和纯超声组织成像中深穿透特性的优点，可得到高分辨率和高对比度的组织图像，从原理上避开了光散射的影响，突破了高分辨率光学成像深度“软极限”，可实现深层活体内组织成像。
光声断层成像（Photoacoustic computed tomography, 简称PACT）是光声成像的一种成像模式，主要针对于较深层组织进行成像。是一种能够清晰获取物体二维断层图像或者三维立体图像的新兴成像技术。该种模式根据光声效应，利用非弹道激光使血流等特异性吸收光的物质产生声学信号，并通过超声换能器接收信号，反演重建出目标的结构和功能信息。在光声断层成像中，通常采用非聚焦大直径脉冲激光束实现组织表面的全场照明。光吸收体对入射光能的吸收导致被照射组织膨胀，随后快速产生宽带超声波。它们传播到组织表面，由机械扫描的非聚焦超声波接收器或接收器阵列检测到，经过放大处理后通过光声图像重建算法对光声信号的传播问题逆向求解。即通过超声换能器检测得到的原始光声信号，反推生物组织成像区域的光吸收分布，从而得到其光声图像。光声层析成像运用了生物组织结构固有的光学对比度，几乎适用于任何具有光吸收特性的生物组织成像。该技术结合了光学成像和超声成像的优点，能够实现分辨率和对比度都较高的生物组织成像。尤其是其非入侵式和非电离式的成像特点，在生物医学领域具有重要应用价值。
光声断层成像技术和设备的发展在一定程度上得益于光声信号仿真软件的研发和应用。目前，在多种光声信号仿真软件中，k-Wave是开发较为完善、应用较为广泛的一款软件。k-Wave软件是MATLAB和C++的开源声学工具箱，由Bradley Treeby和Ben Cox（伦敦大学学院）和Jiri Jaros（布尔诺理工大学）开发。该软件设计用于组织介质中的时域声学和超声仿真。仿真函数基于k空间伪光谱方法，核心是一个可以解释线性和非线性的波传播、非同质材料参数的任意分布和幂律声吸收的数学模型。最重要的功能在于其光声前向模型仿真，通过数值方法求解光声方程，从初始压力分布产生光声信号来模拟光声效应。但是，该软件采用的数值方法复杂，尤其在建立较大尺寸分析模型或者建立较高空间分辨率的模型时，运算时间会大大增加，甚至超出了计算机的运算能力，严重影响了仿真分析的效率，从而制约了该软件的应用场景。
发明内容
针对现有技术不足，本发明的目的是提供一种基于球形粒子光脉冲激发效应的光声信号仿真方法，通过精简的数值计算方式大幅度提高运算效率，以解决k-Wave软件在建立大尺寸或高空间分辨率分析模型时运算效率过低甚至因超过仿真平台能力不能计算的问题。
为达到上述目的，本发明采用的技术方案具体为：
本发明实施例的第一方面提供了一种基于球形粒子光脉冲激发效应的光声信号仿真方法，所述方法包括以下步骤：
S1，设定采样频率和采样时间，得到光声信号传播的时间矢量；根据待测体模所处的传输介质，设定光声信号的传输速度；
S2，设置光声断层成像探测阵列，并计算光声断层成像探测阵列中各阵元的位置坐标；
S3，根据待进行光声断层成像区域的重建分辨率，将待检测体模分解为若干球形粒子，并设定每个球形粒子的初始声压；
S4，根据步骤S1得到的时间矢量、光声信号的传输速度和步骤S3设定的每个球形粒子的初始声压，计算球形粒子向四周辐射的光声信号的声压分布，并将球形粒子向四周辐射的光声信号的声压分布叠加转化为光声断层成像探测阵列中单个阵元接收到的各球形粒子构成的体模所产生的光声信号；
S5，重复步骤S4，将光声断层成像探测阵列中每个阵元接收到的各球形粒子所产生的光声信号组成光声信号矩阵。
进一步地，所述步骤S1中设定采样频率和采样时间，得到光声信号传播的时间矢量包括：根据待测体模和需重建光声断层图像的尺寸、体模和探测阵元之间的最远距离以及需采样的深度在内的需求，设定采样频率和采样时间，从而得到光声信号传播的时间矢量。
进一步地，所述步骤S2中设置光声断层成像探测阵列包括：如果重建二维光声断层图像，则选择二维结构的光声断层成像探测阵列；如果重建三维光声断层图像，则选择三维结构的光声断层成像探测阵列。
进一步地，所述二维结构的光声断层成像探测阵列包括阵元均匀分布的环形阵列；所述三维结构的光声断层成像探测阵列包括碗状阵列。
进一步地，所述碗状阵列中各阵元通过斐波那契格点方法排布。
进一步地，所述步骤S3中根据待进行光声断层成像区域的重建分辨率，待检测体 模分解为若干球形粒子，并设定每个球形粒子的初始声压包括：待进行光声断层成像区域 的重建分辨率为
Figure 656826DEST_PATH_IMAGE001
，则将待测体模共均匀分解为M个半径均为
Figure 42808DEST_PATH_IMAGE002
的球形粒子,并设定各 个粒子的初始声压
Figure 150441DEST_PATH_IMAGE003
,
Figure 14492DEST_PATH_IMAGE004
。
进一步地，所述步骤S4中计算光声断层成像探测阵列中单个阵元接收到的各球形粒子所产生的光声信号的公式如下：
Figure 937187DEST_PATH_IMAGE005
其中，M为球形粒子的个数，
Figure 228491DEST_PATH_IMAGE002
为球形粒子的半径，
Figure 88999DEST_PATH_IMAGE006
为光声信号传播的时间矢量，
Figure 491162DEST_PATH_IMAGE007
为光声信号的传输速度，
Figure 504248DEST_PATH_IMAGE008
为阶跃函数，
Figure 966454DEST_PATH_IMAGE009
是阵元的位置坐标
Figure 314258DEST_PATH_IMAGE010
和球形粒子中心
Figure 520112DEST_PATH_IMAGE011
的距离，
Figure 886240DEST_PATH_IMAGE012
。
进一步地，所述步骤S5中将光声断层成像探测阵列中每个阵元接收到的各球形粒子所产生的光声信号组成光声信号矩阵，光声信号矩阵的公式如下：
Figure 784926DEST_PATH_IMAGE013
其中，
Figure 620027DEST_PATH_IMAGE014
为单个阵元接收到的各球形粒子所产生的光声信号，N为阵元 个数。
本发明实施例的第二方面提供了一种基于球形粒子光脉冲激发效应的光声信号仿真装置，包括一个或多个处理器，用于实现上述的基于球形粒子光脉冲激发效应的光声信号仿真方法。
本发明实施例的第三方面提供了一种计算机可读存储介质，其上存储有程序，该程序被处理器执行时，用于上述的基于球形粒子光脉冲激发效应的光声信号仿真方法。
与现有技术相比，本发明的有益效果为：本发明提供了一种基于球形粒子光脉冲激发效应的光声信号仿真方法，待检测体模分解为若干球形粒子，根据时间矢量、光声信号的传输速度和每个球形粒子的初始声压，计算光声断层成像探测阵列单个阵元接收到的光声信号，将每个阵元得到的光声信号进行拼接得到光声信号矩阵。本发明方法将待测体模受激光照射后产生的光声信号随时间的演化过程进行简化，直接考虑待测体模发出的所有光声信号经一定距离传播后在光声断层成像探测阵列单个阵元上的叠加效果。在不影响仿真信号的准确性的同时，有效降低了数值计算的复杂度，消除了光声信号仿真受到的来自光声断层成像区域的空间尺寸和分辨率的制约，可以快速得到待测体模发送到光声断层成像探测阵列各阵元的光声信号仿真数据，从而大大提高了仿真效率。
附图说明
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例中的光声信号数值计算仿真流程。
图2是本发明实施例中的光声断层成像探测阵列和体模示意图。
图3是本发明实施例中采用k-Wave方法仿真光声信号的结果。
图4是本发明实施例中采用本发明方法仿真光声信号的结果。
图5是本发明实施例中采用k-Wave方法仿真光声信号的xoy截面重建结果。
图6是本发明实施例中采用k-Wave方法仿真光声信号的xoz截面重建结果。
图7是本发明实施例中采用k-Wave方法仿真光声信号的yoz截面重建结果。
图8是本发明实施例中采用本发明方法仿真光声信号的xoy截面重建结果。
图9是本发明实施例中采用本发明方法仿真光声信号的xoz截面重建结果。
图10是本发明实施例中采用本发明方法仿真光声信号的yoz截面重建结果。
图11是本发明实施例提供的基于球形粒子光脉冲激发效应的光声信号仿真装置的结构框图。
具体实施方式
为使本申请的上述目的、特征和优点能够更加明显易懂，下面结合附图对本申请的具体实施方式做详细的说明。在下面的描述中阐述了很多具体细节以便于充分理解本申请。但是本申请能够以很多不同于在此描述的其它方式来实施，本领域技术人员可以在不违背本申请内涵的情况下做类似改进，因此本申请不受下面公开的具体实施的限制。
下面结合附图来具体描述本发明的优选实施例，其中，附图构成本申请一部分，并与本发明的实施例一起用于阐释本发明的原理，并非用于限定本发明的范围。
本发明提出了一种基于球形粒子光脉冲激发效应的光声信号仿真方法，本发明方法的核心是将待检测体模分解为多个微小的球形粒子，根据光声成像相关理论，直接计算得出光声断层成像探测阵列各阵元接收到的上述粒子产生的光声信号。
如果要产生光声信号，通常由一个短脉冲的激光束照射吸收体，吸收体中的分子吸收光子后，当满足一定的条件时，吸收体分子的电子从低能级跃迁到高能级而处于激发态，而处于激发态的电子极不稳定，当电子从高能级向低能级跃迁时，会以光或热量的形式释放能量。在光声成像应用中通常会选择合适波长的激光作为激发源，使吸收的光子的能量转化为热能的效率最大，通常从光能转化为热能的效率可达到90%以上。释放的热量导致吸收体局部温度升高，温度升高后导致热膨胀进一步转换成压力升高，压力升高以超声波的形式传播，这各超声波就被称为光声信号。因此，光声信号的产生过程就是“光能”—“热能”—“机械能”的转化过程。
接下来，如图1所示，对光声断层成像探测阵列各阵元接收到的光声信号进行仿真计算主要包含如下步骤：
S1，设定采样频率和采样时间，得到光声信号传播的时间矢量；根据待测体模所处的传输介质，设定光声信号的传输速度。
具体地，设定仿真所需的基本参数：根据待测体模所处的传输介质的特性，设定好 在该种传输介质中光声信号的传输速度
Figure 363992DEST_PATH_IMAGE015
。并且根据数值计算的需求，如待测体模和需重 建光声断层图像的尺寸，体模和探测阵元之间的最远距离以及需采样的深度等，设定采样 频率和采样时间，从而得到光声信号传播的时间矢量
Figure 273042DEST_PATH_IMAGE006
。
S2，设置光声断层成像探测阵列，并计算光声断层成像探测阵列中各阵元的位置坐标。
具体地，设计光声断层成像的探测阵列：根据仿真要求，如果重建二维光声断层图像，则选择二维结构的光声断层成像探测阵列，如阵元均匀分布的环形阵列；如果重建三维光声断层图像，则选择三维结构的光声断层成像探测阵列，如碗状阵列。此外，还可根据待测体模和需重建光声断层图像的尺寸，以及成像的分辨率等设置光声断层成像探测阵列，计算出光声断层成像阵列各阵元的位置坐标。
S3，根据待进行光声断层成像区域的重建分辨率，将待检测体模分解为若干球形粒子，并设定每个球形粒子的初始声压。
具体地，进行待测体模粒子分解的步骤包括：为了方便后续的数值计算，本发明依 据波的叠加原理对待测体模受激光激发出的光声信号的物理描述进行简化。波的叠加原理 是指，介质中同时存在几列波时，每列波能保持各自的传播规律而不互相干扰。在波的重叠 区域里各点的振动的物理量等于各列波在该点引起的物理量的矢量和。由于要仿真的是各 阵元接收到的待测体模传播光声信号的声压，该物理量为标量。因此，可以认为单个阵元处 接收到的待测体模的光声信号的声压为构成待测体模的数个粒子传播到单个阵元处的光 声信号的声压之和。如果待进行光声断层成像区域的重建分辨率为
Figure 952416DEST_PATH_IMAGE001
，则可将待测体模 在空间坐标系中所占据的区域按照光声断层成像能够重建的最小尺度分解为M个半径均为
Figure 150179DEST_PATH_IMAGE002
的球形粒子,设定好各个粒子的初始声压
Figure 822469DEST_PATH_IMAGE003
,
Figure 789288DEST_PATH_IMAGE004
。
S4，根据步骤S1得到的时间矢量、光声信号的传输速度和步骤S3设定的每个球形粒子的初始声压，计算球形粒子向四周辐射的光声信号的声压分布，并将球形粒子向四周辐射的光声信号的声压分布叠加转化为光声断层成像探测阵列中单个阵元接收到的各球形粒子构成的体模所产生的光声信号。
本发明将球形粒子向四周辐射的光声信号的声压分布
Figure 409537DEST_PATH_IMAGE016
叠加转化为阵元 接收到的各球形粒子所产生的光声信号
Figure 829017DEST_PATH_IMAGE017
，本发明方法将待测体模受激光照射后 产生的光声信号随时间的演化过程进行精简，直接考虑待测体模发出的所有光声信号经一 定距离传播后在光声断层成像探测阵列单个阵元上的叠加效果。
球形粒子向四周辐射的光声信号的声压分布
Figure 304998DEST_PATH_IMAGE016
的计算过程如下：
当一个半径为
Figure 126323DEST_PATH_IMAGE018
的微小球形粒子被
Figure 413079DEST_PATH_IMAGE019
形式的光脉冲加热，球体内部就会产生初 始声压
Figure 54276DEST_PATH_IMAGE020
。假设球体被无吸收的背景介质所包围，设
Figure 333948DEST_PATH_IMAGE021
表示球形粒子吸收体中心所在的 位置，则球形粒子被
Figure 9780DEST_PATH_IMAGE019
形式的光脉冲照射激发并向粒子四周辐射的光声信号的声压分布可 以用下列公示表示：
Figure 231552DEST_PATH_IMAGE022
其中，
Figure 360045DEST_PATH_IMAGE023
是阵元的位置坐标
Figure 177828DEST_PATH_IMAGE010
和球形粒子中心
Figure 708166DEST_PATH_IMAGE011
的距离，
Figure 336725DEST_PATH_IMAGE024
。
Figure 218093DEST_PATH_IMAGE015
是 光声信号的传播速度，
Figure 839567DEST_PATH_IMAGE006
是光声信号传播的时间。
Figure 489992DEST_PATH_IMAGE008
为阶跃函数，当
Figure 522407DEST_PATH_IMAGE025
时，
Figure 891072DEST_PATH_IMAGE026
；当
Figure 722762DEST_PATH_IMAGE027
时，
Figure 617905DEST_PATH_IMAGE028
。
将球形粒子向四周辐射的光声信号的声压分布叠加转化为光声断层成像探测阵列中单个阵元接收到的各球形粒子构成的体模所产生的光声信号的过程具体包括：
叠加计算得到光声断层成像探测阵列第
Figure 447321DEST_PATH_IMAGE029
个阵元接收到的各球形粒子所产生的光 声信号
Figure 178648DEST_PATH_IMAGE030
，公式如下：
Figure 17291DEST_PATH_IMAGE031
其中，M为球形粒子的个数，
Figure 501362DEST_PATH_IMAGE002
为球形粒子的半径，
Figure 767258DEST_PATH_IMAGE006
为光声信号传播的时间矢量，
Figure 953258DEST_PATH_IMAGE007
为光声信号的传输速度，
Figure 595592DEST_PATH_IMAGE008
为阶跃函数，
Figure 199748DEST_PATH_IMAGE009
是阵元的位置坐标
Figure 636546DEST_PATH_IMAGE010
和球形粒子中心
Figure 76886DEST_PATH_IMAGE011
的距离，
Figure 257331DEST_PATH_IMAGE012
。
S5，重复步骤S4，将光声断层成像探测阵列中每个阵元接收到的各球形粒子所产生的光声信号组成光声信号矩阵。
具体包括以下步骤：
S501，设定阵元序号循环数。即先设定阵元序号循环数i=1。每重复步骤S4一次，则i=i
=1。重复步骤S4，并判断循环次数，若i<N，N为阵元个数，则返回S4直至计算完所有阵元接收到的各球形粒子所产生的光声信号，否则进入S502。
S502，求解得到光声断层成像探测阵列光声信号矩阵：通过上述步骤，将光声断层 成像探测阵列中每个阵元接收到的各球形粒子所产生的光声信号组成光声信号矩阵，最终 得到光声断层成像探测阵列所有阵元接收到的光声信号矩阵
Figure 715994DEST_PATH_IMAGE032
：
Figure 323693DEST_PATH_IMAGE033
。
实施例1：本发明实例描述的是通过该方法仿真计算512元碗型阵列接收的来自阵列内部相互垂直的三根线段组成的体模的光声信号。受限于k-Wave仿真对计算资源的要求，碗型阵列的半径设为20mm，包含阵列在内的区域设为40 mm *40 mm*40 mm，重建区域分辨率设为0.1mm。
使用本发明提供的方法按照如下步骤进行仿真：
S1，设定仿真所需的基本参数。根据数值计算的需求，设定采样频率f s =50MHz和采样时间T=50us，得到光声信号传播的时间矢量t=(0:1/fs:T)。并且，根据待测体模所处的传输介质的特性，设定光声信号的传输速度v=0.15cm/us。
S2，完成用于进行光声断层成像的探测阵列的设计。根据仿真要求，计算出光声断层成像阵列各阵元的位置坐标。该实施例中，如图2所示，仿真的是512元碗型超声探测阵列，要求各阵元均匀分布在底部具有规定半径圆孔的半球面上。该实施例中通过斐波那契格点方法构造出均匀分布在设定的碗状球面上的所有阵元。
S3，进行待测体模粒子分解。该实施例中仿真的体模为位于阵列内部分别平行于 笛卡尔坐标系x、y、z坐标轴并且相互垂直的三个线段。假设待进行光声断层成像区域的重 建分辨率为
Figure 172700DEST_PATH_IMAGE001
，则将三根线段构成的待测体模共均匀分解为M=603个半径均为
Figure 282737DEST_PATH_IMAGE034
的球形 粒子,设定好各个粒子的初始声压
Figure 736852DEST_PATH_IMAGE003
,
Figure 640086DEST_PATH_IMAGE004
。本发明实施例1中，重建区域分 辨率设为0.1mm，则球形粒子的半径
Figure 179652DEST_PATH_IMAGE002
为0.05mm。
S4，设定阵元序号循环数。为了通过循环方式获得各阵元的光声信号，先设定阵元序号循环数i=1。
S5，数值计算光声断层成像探测阵列单个阵元接收到的光声信号。通过数值计算 得到光声断层成像探测阵列第
Figure 577266DEST_PATH_IMAGE029
个阵元接收到的各球形粒子所产生的光声信号
Figure 885888DEST_PATH_IMAGE030
：
Figure 960023DEST_PATH_IMAGE031
其中，
Figure 721306DEST_PATH_IMAGE035
。
S6，进行阵元信号计算循环判断。判断是否全部生成各阵元光声信号。首先i=i+1,若i<512，则返回S5，否则进入S7。
S7，求解得到光声断层成像探测阵列光声信号矩阵。通过上述步骤，最终可以得到 光声断层成像探测阵列所有阵元接收到的光声信号矩阵
Figure 155567DEST_PATH_IMAGE036
：
Figure 318695DEST_PATH_IMAGE037
。
本发明首先对比上述方法和k-Wave仿真的光声信号，如图3和图4所示。可以看出，本发明提供的方法可以得到和k-Wave基本相同的仿真信号。两种方法产生的信号无论在信号波形还是幅度方面均十分接近。再对两种方法产生的仿真信号用同一种方法进行光声断层图像重建，如图5~图10所示，对两种方法产生的光声断层成像阵列接收的光声信号进行处理能够得到几乎相同的重建结果。因此，本发明提供的方法能够实现和k-Wave同等质量的仿真计算效果。
进一步，对比两种方法对同一仿真条件的数值计算时间如表1所示：
表1：两种方法的数值计算时间比较表
Figure 563732DEST_PATH_IMAGE038
对同样规模的分析模型进行数值计算，本发明的仿真用时远远少于k-Wave软件。可见本发明提供的方法的计算效率远高于k-Wave软件。
进一步地，在进行上述仿真的计算机上继续开展数值计算，选择更大的仿真区域,或者选择更高的重建区域分辨率,则k-Wave软件仿真会显示计算机内存不够的错误提示并且终止运算。而本发明提供的方法则可以继续正确进行。
综上所述，本发明所提出的基于球形粒子光脉冲激发效应的光声信号仿真方法，由于其仿真机理不受仿真区域尺寸和重建区域分辨率等因素的限制，能够在更大规模或更精细的仿真条件下具有较高的运算效率和运算能力，相较于k-Wave软件等其他光声信号仿真方法具有明显优势。本发明提出的方法可以应用于较大规模复杂结构光声断层成像阵列的成像能力和成像性能的评估，可用于较大尺寸待测体模的光声断层图像重建结果分析，也可用于高分辨率光声断层图像重建情形下的数值分析。因此，本发明提出的方法在光声断层成像技术研究中具有较为广泛且实际的应用价值。
与前述基于球形粒子光脉冲激发效应的光声信号仿真方法的实施例相对应，本发明还提供了基于球形粒子光脉冲激发效应的光声信号装置的实施例。
参见图11，本发明实施例提供的一种基于球形粒子光脉冲激发效应的光声信号装置，包括一个或多个处理器，用于实现上述实施例中的基于球形粒子光脉冲激发效应的光声信号方法。
本发明基于球形粒子光脉冲激发效应的光声信号装置的实施例可以应用在任意具备数据处理能力的设备上，该任意具备数据处理能力的设备可以为诸如计算机等设备或装置。装置实施例可以通过软件实现，也可以通过硬件或者软硬件结合的方式实现。以软件实现为例，作为一个逻辑意义上的装置，是通过其所在任意具备数据处理能力的设备的处理器将非易失性存储器中对应的计算机程序指令读取到内存中运行形成的。从硬件层面而言，如图11所示，为本发明基于球形粒子光脉冲激发效应的光声信号装置所在任意具备数据处理能力的设备的一种硬件结构图，除了图11所示的处理器、内存、网络接口、以及非易失性存储器之外，实施例中装置所在的任意具备数据处理能力的设备通常根据该任意具备数据处理能力的设备的实际功能，还可以包括其他硬件，对此不再赘述。
上述装置中各个单元的功能和作用的实现过程具体详见上述方法中对应步骤的实现过程，在此不再赘述。
对于装置实施例而言，由于其基本对应于方法实施例，所以相关之处参见方法实施例的部分说明即可。以上所描述的装置实施例仅仅是示意性的，其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本发明方案的目的。本领域普通技术人员在不付出创造性劳动的情况下，即可以理解并实施。
本发明实施例还提供一种计算机可读存储介质，其上存储有程序，该程序被处理器执行时，实现上述实施例中的基于球形粒子光脉冲激发效应的光声信号方法。
所述计算机可读存储介质可以是前述任一实施例所述的任意具备数据处理能力的设备的内部存储单元，例如硬盘或内存。所述计算机可读存储介质也可以是任意具备数据处理能力的设备，例如所述设备上配备的插接式硬盘、智能存储卡（Smart Media Card，SMC）、SD卡、闪存卡（Flash Card）等。进一步的，所述计算机可读存储介质还可以既包括任意具备数据处理能力的设备的内部存储单元也包括外部存储设备。所述计算机可读存储介质用于存储所述计算机程序以及所述任意具备数据处理能力的设备所需的其他程序和数据，还可以用于暂时地存储已经输出或者将要输出的数据。
本领域普通技术人员可以理解，以上所述仅为发明的优选实例而已，并不用于限制发明，尽管参照前述实例对发明进行了详细的说明，对于本领域的技术人员来说，其依然可以对前述各实例记载的技术方案进行修改，或者对其中部分技术特征进行等同替换。凡在发明的精神和原则之内，所做的修改、等同替换等均应包含在发明的保护范围之内。

Claims (10)

1.一种基于球形粒子光脉冲激发效应的光声信号仿真方法，其特征在于，所述方法包括以下步骤：
S1，设定采样频率和采样时间，得到光声信号传播的时间矢量；根据待测体模所处的传输介质，设定光声信号的传输速度；
S2，设置光声断层成像探测阵列，并计算光声断层成像探测阵列中各阵元的位置坐标；
S3，根据待进行光声断层成像区域的重建分辨率，将待检测体模分解为若干球形粒子，并设定每个球形粒子的初始声压；
S4，根据步骤S1得到的时间矢量、光声信号的传输速度和步骤S3设定的每个球形粒子的初始声压，计算球形粒子向四周辐射的光声信号的声压分布，并将球形粒子向四周辐射的光声信号的声压分布叠加转化为光声断层成像探测阵列中单个阵元接收到的各球形粒子构成的体模所产生的光声信号；
S5，重复步骤S4，将光声断层成像探测阵列中每个阵元接收到的各球形粒子所产生的光声信号组成光声信号矩阵。
2.根据权利要求1所述的基于球形粒子光脉冲激发效应的光声信号仿真方法，其特征在于，所述步骤S1中设定采样频率和采样时间，得到光声信号传播的时间矢量包括：根据待测体模和需重建光声断层图像的尺寸、体模和探测阵元之间的最远距离以及需采样的深度在内的需求，设定采样频率和采样时间，从而得到光声信号传播的时间矢量。
3.根据权利要求1所述的基于球形粒子光脉冲激发效应的光声信号仿真方法，其特征在于，所述步骤S2中设置光声断层成像探测阵列包括：如果重建二维光声断层图像，则选择二维结构的光声断层成像探测阵列；如果重建三维光声断层图像，则选择三维结构的光声断层成像探测阵列。
4.根据权利要求3所述的基于球形粒子光脉冲激发效应的光声信号仿真方法，其特征在于，所述二维结构的光声断层成像探测阵列包括阵元均匀分布的环形阵列；所述三维结构的光声断层成像探测阵列包括碗状阵列。
5.根据权利要求4所述的基于球形粒子光脉冲激发效应的光声信号仿真方法，其特征在于，所述碗状阵列中各阵元通过斐波那契格点方法排布。
6.根据权利要求1所述的基于球形粒子光脉冲激发效应的光声信号仿真方法，其特征 在于，所述步骤S3中根据待进行光声断层成像区域的重建分辨率，待检测体模分解为若干 球形粒子，并设定每个球形粒子的初始声压包括：待进行光声断层成像区域的重建分辨率 为
Figure 979226DEST_PATH_IMAGE001
，则将待测体模共均匀分解为M个半径均为
Figure 108856DEST_PATH_IMAGE002
的球形粒子,并设定各个粒子的初始 声压
Figure 719966DEST_PATH_IMAGE003
,
Figure 656610DEST_PATH_IMAGE004
。
7.根据权利要求1或4所述的基于球形粒子光脉冲激发效应的光声信号仿真方法，其特征在于，所述步骤S4中计算光声断层成像探测阵列中单个阵元接收到的各球形粒子所产生的光声信号的公式如下：
Figure 720381DEST_PATH_IMAGE005
其中，M为球形粒子的个数，
Figure 388123DEST_PATH_IMAGE002
为球形粒子的半径，
Figure 135630DEST_PATH_IMAGE006
为光声信号传播的时间矢量，
Figure 456890DEST_PATH_IMAGE007
为 光声信号的传输速度，
Figure 945640DEST_PATH_IMAGE008
为阶跃函数，
Figure 525395DEST_PATH_IMAGE009
是阵元的位置坐标
Figure 517622DEST_PATH_IMAGE010
和球形粒子中心
Figure 540942DEST_PATH_IMAGE011
的距 离，
Figure 985829DEST_PATH_IMAGE012
。
8.根据权利要求1所述的基于球形粒子光脉冲激发效应的光声信号仿真方法，其特征在于，所述步骤S5中将光声断层成像探测阵列中每个阵元接收到的各球形粒子所产生的光声信号组成光声信号矩阵，光声信号矩阵的公式如下：
Figure 870740DEST_PATH_IMAGE013
其中，
Figure 983052DEST_PATH_IMAGE014
为单个阵元接收到的各球形粒子所产生的光声信号，N为阵元个数。
9.一种基于球形粒子光脉冲激发效应的光声信号仿真装置，其特征在于，包括一个或多个处理器，用于实现权利要求1-8中任一项所述的基于球形粒子光脉冲激发效应的光声信号仿真方法。
10.一种计算机可读存储介质，其上存储有程序，其特征在于，该程序被处理器执行时，用于实现权利要求1-8中任一项所述的基于球形粒子光脉冲激发效应的光声信号仿真方法。
CN202211100325.8A 2022-09-09 2022-09-09 基于球形粒子光脉冲激发效应的光声信号仿真方法、装置 Active CN115177217B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211100325.8A CN115177217B (zh) 2022-09-09 2022-09-09 基于球形粒子光脉冲激发效应的光声信号仿真方法、装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211100325.8A CN115177217B (zh) 2022-09-09 2022-09-09 基于球形粒子光脉冲激发效应的光声信号仿真方法、装置

Publications (2)

Publication Number Publication Date
CN115177217A true CN115177217A (zh) 2022-10-14
CN115177217B CN115177217B (zh) 2023-01-03

Family

ID=83524586

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211100325.8A Active CN115177217B (zh) 2022-09-09 2022-09-09 基于球形粒子光脉冲激发效应的光声信号仿真方法、装置

Country Status (1)

Country Link
CN (1) CN115177217B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115453459A (zh) * 2022-11-10 2022-12-09 浙江大学湖州研究院 一种基于多模态声信号修正的应急定位方法及系统
CN115619889A (zh) * 2022-11-09 2023-01-17 哈尔滨工业大学（威海） 一种适用于环形阵列的多特征融合光声图像重建方法
CN116468859A (zh) * 2023-06-19 2023-07-21 之江实验室 适用于声速非均匀分布的三维光声断层成像方法、装置
CN116482035A (zh) * 2023-06-21 2023-07-25 之江实验室 一种基于柔性超声探头的光声断层成像方法及装置
CN116740219A (zh) * 2023-08-14 2023-09-12 之江实验室 一种三维光声断层成像方法、装置、设备及可读存储介质
CN117113793A (zh) * 2023-10-23 2023-11-24 之江实验室 磁约束带电粒子成像系统中粒子输运的模拟方法
CN117828957A (zh) * 2024-03-05 2024-04-05 之江实验室 一种超声弹性仿真方法、装置、存储介质及电子设备

Citations (25)

* Cited by examiner, † Cited by third party
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
US20080123083A1 (en) * 2006-11-29 2008-05-29 The Regents Of The University Of Michigan System and Method for Photoacoustic Guided Diffuse Optical Imaging
US20110282181A1 (en) * 2009-11-12 2011-11-17 Ge Wang Extended interior methods and systems for spectral, optical, and photoacoustic imaging
WO2012138965A2 (en) * 2011-04-08 2012-10-11 University Of Florida Research Foundation, Inc. Enhanced image reconstruction in photoacoustic tomography
US20130102865A1 (en) * 2011-10-25 2013-04-25 Andreas Mandelis Systems and methods for frequency-domain photoacoustic phased array imaging
WO2013076986A1 (ja) * 2011-11-22 2013-05-30 富士フイルム株式会社 光音響信号処理装置及び方法
WO2013134782A1 (en) * 2012-03-09 2013-09-12 The Johns Hopkins University Photoacoustic tracking and registration in interventional ultrasound
CN104055541A (zh) * 2014-06-26 2014-09-24 中国科学院苏州生物医学工程技术研究所 一种用于血管内超声多断层剪切波弹性成像方法
US20150245771A1 (en) * 2012-10-18 2015-09-03 Washington University Transcranial photoacoustic/thermoacoustic tomography brain imaging informed by adjunct image data
US20160022150A1 (en) * 2014-07-24 2016-01-28 Canon Kabushiki Kaisha Photoacoustic apparatus
WO2017210675A1 (en) * 2016-06-03 2017-12-07 The General Hospital Corporation System and method for micro laser particles
US20180014732A1 (en) * 2015-01-07 2018-01-18 Canon Kabushiki Kaisha Photoacoustic apparatus, image display method, and program
CN107713990A (zh) * 2017-10-31 2018-02-23 华南师范大学 一种热声、光声、超声三模态乳腺肿瘤检测装置及方法
JP2018033886A (ja) * 2016-09-02 2018-03-08 国立大学法人京都大学 光音響情報処理方法、光音響情報処理プログラムおよび光音響情報処理装置
CN108272439A (zh) * 2017-01-06 2018-07-13 佳能株式会社 被检体信息处理装置和图像显示方法
CN108896489A (zh) * 2018-07-24 2018-11-27 广东工业大学 一种光声成像方法和装置
CN109864707A (zh) * 2019-01-17 2019-06-11 南京科技职业学院 一种在有限视角情况下提高光声断层成像分辨率的方法
CN110456362A (zh) * 2019-07-17 2019-11-15 北京大学 一种基于脉冲对发射的目标声成像及测速方法与系统
JP2020039841A (ja) * 2017-11-24 2020-03-19 炭　親良 計測イメージング装置
US20200163649A1 (en) * 2017-06-05 2020-05-28 Mayo Foundation For Medical Education And Research Estimating Phase Velocity Dispersion in Ultrasound Elastography Using a Multiple Signal Classification
CN112180329A (zh) * 2020-09-07 2021-01-05 黑龙江工程学院 一种基于阵元随机均匀分布球阵反卷积波束形成的汽车噪声源声成像方法
CN112229799A (zh) * 2020-09-21 2021-01-15 华南师范大学 偏振光声穆勒矩阵成像方法和系统
US20210128238A1 (en) * 2019-11-04 2021-05-06 Phononz Inc Mediator-free universal laser light amplification with coaxial propagating focused ultrasound and system
CN113777045A (zh) * 2020-06-10 2021-12-10 复旦大学 基于单粒子多边定位追踪的超分辨功能性光声成像方法
CN114563479A (zh) * 2022-04-20 2022-05-31 之江实验室 一种实时三维高分辨太赫兹光声成像方法和装置

Patent Citations (25)

* Cited by examiner, † Cited by third party
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
US20080123083A1 (en) * 2006-11-29 2008-05-29 The Regents Of The University Of Michigan System and Method for Photoacoustic Guided Diffuse Optical Imaging
US20110282181A1 (en) * 2009-11-12 2011-11-17 Ge Wang Extended interior methods and systems for spectral, optical, and photoacoustic imaging
WO2012138965A2 (en) * 2011-04-08 2012-10-11 University Of Florida Research Foundation, Inc. Enhanced image reconstruction in photoacoustic tomography
US20130102865A1 (en) * 2011-10-25 2013-04-25 Andreas Mandelis Systems and methods for frequency-domain photoacoustic phased array imaging
WO2013076986A1 (ja) * 2011-11-22 2013-05-30 富士フイルム株式会社 光音響信号処理装置及び方法
WO2013134782A1 (en) * 2012-03-09 2013-09-12 The Johns Hopkins University Photoacoustic tracking and registration in interventional ultrasound
US20150245771A1 (en) * 2012-10-18 2015-09-03 Washington University Transcranial photoacoustic/thermoacoustic tomography brain imaging informed by adjunct image data
CN104055541A (zh) * 2014-06-26 2014-09-24 中国科学院苏州生物医学工程技术研究所 一种用于血管内超声多断层剪切波弹性成像方法
US20160022150A1 (en) * 2014-07-24 2016-01-28 Canon Kabushiki Kaisha Photoacoustic apparatus
US20180014732A1 (en) * 2015-01-07 2018-01-18 Canon Kabushiki Kaisha Photoacoustic apparatus, image display method, and program
WO2017210675A1 (en) * 2016-06-03 2017-12-07 The General Hospital Corporation System and method for micro laser particles
JP2018033886A (ja) * 2016-09-02 2018-03-08 国立大学法人京都大学 光音響情報処理方法、光音響情報処理プログラムおよび光音響情報処理装置
CN108272439A (zh) * 2017-01-06 2018-07-13 佳能株式会社 被检体信息处理装置和图像显示方法
US20200163649A1 (en) * 2017-06-05 2020-05-28 Mayo Foundation For Medical Education And Research Estimating Phase Velocity Dispersion in Ultrasound Elastography Using a Multiple Signal Classification
CN107713990A (zh) * 2017-10-31 2018-02-23 华南师范大学 一种热声、光声、超声三模态乳腺肿瘤检测装置及方法
JP2020039841A (ja) * 2017-11-24 2020-03-19 炭　親良 計測イメージング装置
CN108896489A (zh) * 2018-07-24 2018-11-27 广东工业大学 一种光声成像方法和装置
CN109864707A (zh) * 2019-01-17 2019-06-11 南京科技职业学院 一种在有限视角情况下提高光声断层成像分辨率的方法
CN110456362A (zh) * 2019-07-17 2019-11-15 北京大学 一种基于脉冲对发射的目标声成像及测速方法与系统
US20210128238A1 (en) * 2019-11-04 2021-05-06 Phononz Inc Mediator-free universal laser light amplification with coaxial propagating focused ultrasound and system
CN113777045A (zh) * 2020-06-10 2021-12-10 复旦大学 基于单粒子多边定位追踪的超分辨功能性光声成像方法
CN112180329A (zh) * 2020-09-07 2021-01-05 黑龙江工程学院 一种基于阵元随机均匀分布球阵反卷积波束形成的汽车噪声源声成像方法
CN112229799A (zh) * 2020-09-21 2021-01-15 华南师范大学 偏振光声穆勒矩阵成像方法和系统
CN114563479A (zh) * 2022-04-20 2022-05-31 之江实验室 一种实时三维高分辨太赫兹光声成像方法和装置

Non-Patent Citations (9)

* Cited by examiner, † Cited by third party
Title
BONDITA PAUL: "Parametric Reconstruction of Photoacoustic Tomographic Imaging using Gaussian Mixture Model and Evolutionary Methods", 《2021 INTERNATIONAL CONFERENCE ON INTELLIGENT TECHNOLOGIES (CONIT)》 *
JIN WANG等: "A simulation study of the effect of transducer position on photoacoustic spectrum analysis for stochastic microstructure", 《2015 8TH INTERNATIONAL CONFERENCE ON BIOMEDICAL ENGINEERING AND INFORMATICS (BMEI)》 *
PENGFEI ZHANG: "High-resolution Deep Functional Imaging of the Whole Mouse Brain by Photoacoustic Computed Tomography in Vivo", 《JOURNAL OF BIOPHOTONICS》 *
RAVI PRAKASH: "Photoacoustic Signal Simulation Using Discrete Particle Approach and its Application in Tomography", 《IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS, AND FREQUENCY CONTROL》 *
刘高尚: "基于偏振参数成像方法的声场可视化研究", 《中国优秀硕士学位论文全文数据库（工程科技Ⅱ辑）》 *
卢涛等: "基于阵列探测方式的时域光声成像系统", 《纳米技术与精密工程》 *
徐晓辉等: "生物医学光声成像", 《物理》 *
戴娜: "基于非负约束脉冲分解的光声信号与图像增强方法研究", 《中国优秀硕士学位论文全文数据库（信息科技辑）》 *
沈康等: "基于双域神经网络的稀疏视角光声图像重建", 《中国激光》 *

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115619889A (zh) * 2022-11-09 2023-01-17 哈尔滨工业大学（威海） 一种适用于环形阵列的多特征融合光声图像重建方法
CN115619889B (zh) * 2022-11-09 2023-05-30 哈尔滨工业大学（威海） 一种适用于环形阵列的多特征融合光声图像重建方法
CN115453459A (zh) * 2022-11-10 2022-12-09 浙江大学湖州研究院 一种基于多模态声信号修正的应急定位方法及系统
CN116468859A (zh) * 2023-06-19 2023-07-21 之江实验室 适用于声速非均匀分布的三维光声断层成像方法、装置
CN116468859B (zh) * 2023-06-19 2023-09-08 之江实验室 适用于声速非均匀分布的三维光声断层成像方法、装置
CN116482035A (zh) * 2023-06-21 2023-07-25 之江实验室 一种基于柔性超声探头的光声断层成像方法及装置
CN116482035B (zh) * 2023-06-21 2023-11-17 之江实验室 一种基于柔性超声探头的光声断层成像方法及装置
CN116740219A (zh) * 2023-08-14 2023-09-12 之江实验室 一种三维光声断层成像方法、装置、设备及可读存储介质
CN116740219B (zh) * 2023-08-14 2024-01-09 之江实验室 一种三维光声断层成像方法、装置、设备及可读存储介质
CN117113793A (zh) * 2023-10-23 2023-11-24 之江实验室 磁约束带电粒子成像系统中粒子输运的模拟方法
CN117113793B (zh) * 2023-10-23 2024-01-26 之江实验室 磁约束带电粒子成像系统中粒子输运的模拟方法
CN117828957A (zh) * 2024-03-05 2024-04-05 之江实验室 一种超声弹性仿真方法、装置、存储介质及电子设备

Also Published As

Publication number Publication date
CN115177217B (zh) 2023-01-03

Similar Documents

Publication Publication Date Title
CN115177217B (zh) 基于球形粒子光脉冲激发效应的光声信号仿真方法、装置
CN102843960B (zh) 被检体信息获得设备和方法
CN103260502B (zh) 被检体信息获取设备和被检体信息获取方法
CN102258386B (zh) 显示数据获得设备和显示数据获得方法
CN103025248B (zh) 图像信息获取装置和图像信息获取方法
Jin et al. Fast and high-resolution three-dimensional hybrid-domain photoacoustic imaging incorporating analytical-focused transducer beam amplitude
CN102821677B (zh) 光声成像装置和光声成像方法
US9572531B2 (en) Object information acquiring apparatus and control method thereof
Awasthi et al. Sinogram super-resolution and denoising convolutional neural network (SRCN) for limited data photoacoustic tomography
CN104053402B (zh) 被检体信息获得设备和关于被检体的信息的获得方法
CN105078408A (zh) 光声装置
JP6168802B2 (ja) 処理装置、処理方法、及びプログラム
Akhlaghi et al. Multidomain computational modeling of photoacoustic imaging: verification, validation, and image quality prediction
CN116468859B (zh) 适用于声速非均匀分布的三维光声断层成像方法、装置
Saha Solving time-independent inhomogeneous optoacoustic wave equation numerically with a modified Green's function approach
Sun et al. An iterative gradient convolutional neural network and its application in endoscopic photoacoustic image formation from incomplete acoustic measurement
CN106687028A (zh) 光声装置和信息获取装置
CN103815929A (zh) 被检体信息获取装置
Lou et al. Analysis of the use of unmatched backward operators in iterative image reconstruction with application to three-dimensional optoacoustic tomography
Bychkov et al. Toroidally focused sensor array for real-time laser-ultrasonic imaging: The first experimental study
Sahlström et al. Computationally efficient forward operator for photoacoustic tomography based on coordinate transformations
JP2011120780A (ja) 光音響イメージング装置
Tordera Mora et al. Generalized spatial coherence reconstruction for photoacoustic computed tomography
CN113827277B (zh) 一种声致超声成像方法
Tsysar et al. Using a multi-rod waveguide system to create an ultrasound endoscope for imaging in aggressive liquids

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