CN105640496B - 光声装置和被检体信息取得方法 - Google Patents
光声装置和被检体信息取得方法 Download PDFInfo
- Publication number
- CN105640496B CN105640496B CN201510846561.8A CN201510846561A CN105640496B CN 105640496 B CN105640496 B CN 105640496B CN 201510846561 A CN201510846561 A CN 201510846561A CN 105640496 B CN105640496 B CN 105640496B
- Authority
- CN
- China
- Prior art keywords
- subject
- period
- subject information
- signal
- unit
- 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.)
- Active
Links
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
公开了光声装置和被检体信息取得方法。光声装置包括:被配置为利用脉冲光多次照射被检体的光照射单元;被配置为接收通过利用来自光照射单元的脉冲光多次照射被检体产生的光声波并且输出与所述多次光照射对应的多个信号的接收单元;被配置为取得关于被检体的血液量的信息的血液信息取得单元;及被配置为基于所述多个信号取得在被检体中的目标区域的被检体信息的被检体信息取得单元,其中,被检体信息取得单元基于在血液量的重复的波动周期中的各波动周期内的共同时段中取得的多个信号来取得被检体信息。
Description
技术领域
本发明涉及利用光声效应的光声装置。
背景技术
光声成像技术是一种使用光的成像技术。在光声成像中,首先,利用由光源产生的脉冲光照射被检体。照射光被传播和扩散到被检体中并且在被检体的多个部位处被吸收以产生光声波。换能器将光声波转换成电信号,并且处理装置对电信号执行分析处理来取得关于在被检体中的光学特征值的信息。
从被检体中的光吸收体产生的光声波的产生的声压P0(在下文中也被称为“初始声压”)可以通过以下公式来表达:
P0=Γ·μa·Γ·φ (1),
其中,Γ是Gruneisen系数,它是将体积膨胀系数β与声速c的平方的乘积除以在恒定压力Cp下的比热的商,并且Φ是在某位置处(局部区域)的光量(已到达吸收体的光量,其也被称为“光注量”)。
初始声压P0可以利用从已接收到光声波的探头输出的接收信号(光声(PA)信号)来计算。
已经知道如果组织被确定,则Gruneisen系数的值是基本上恒定的。因此,光吸收系数μa与光量Φ(即,光能量吸收密度)的乘积可以通过在多个部分处测量和分析PA信号中的时间变化来获得。
日本专利申请特许公开号2013-248077讨论了被配置为基于由光诱导的光声波生成血管的光声图像的光声图像生成装置。
同时,在其中测量目标是活生物体的情况下,通过光声测量获得的信号会受到活生物体的搏动的影响。例如,在其中光吸收体是血红蛋白的情况下,如果信号是在血管中的血液量多的时间取得的,那么由于在被测区域中存在的血红蛋白的量多,因此根据公式(1),所产生的光声波的声压高。因此,可以预测,获得的信号的信噪比(S/N)相对高。另一方面,如果信号是在其中血管中的血液量少的时间取得的,那么由于在被测区域中存在的血红蛋白量少,因此即使被测区域是相同的,所产生的光声波的声压也低。因此,可以预测,在其中血液量少的时间产生的光声波的接收信号的S/N比相对低。换句话说,所取得的被检体信息的精度取决于血管中的血液量而不同。
发明内容
根据本发明的一方面,光声装置包括:被配置为利用脉冲光多次照射被检体的光照射单元;被配置为接收通过利用来自光照射单元的脉冲光多次照射被检体产生的光声波并且输出与所述多次光照射对应的多个信号的接收单元;被配置为取得关于被检体的血液量的信息的血液信息取得单元;及被配置为基于所述多个信号取得在被检体中的目标区域的被检体信息的被检体信息取得单元,其中,被检体信息取得单元基于在血液量的重复的波动周期中的各波动周期内的共同(common)时段中取得的多个信号来取得被检体信息。
参考附图,本发明的其它特征将从以下示例性实施例的说明中变得清楚。
附图说明
图1示出了根据示例性实施例的光声装置的构造。
图2是示出根据第一示例性实施例的被检体信息取得方法的流程图。
图3A至3D示出了根据第一示例性实施例的各种类型的序列。
图4是示出根据第二示例性实施例的被检体信息取得方法的流程图。
图5示出了根据第二示例性实施例的显示例子。
具体实施方式
本发明的各种示例性实施例将在以下参考附图进行详细描述。原则上,相同的部件被给予相同的标号,并且其说明被省略。
如上所述,在其中血液量少的区域中产生的光声波的接收信号的S/N比相对低。因此,当光声装置取得目标区域的被检体信息时,在其中血液量少的时段中取得的被检体信息的精度可能低。在第一示例性实施例中,将给出其中基于心电图信号来估计在被检体区域中的血液量并且基于在其中血液量相对多的时段中产生的声波来取得被检体信息的例子的描述。
根据本示例性实施例的光声装置是被配置为基于光声波的接收信号取得被检体信息的装置。根据本示例性实施例的被检体信息指从由光声效应产生的光声波的接收信号取得的关于被检体的信息。具体而言,被检体信息是产生的声压(初始声压)、光能量吸收密度、光吸收系数、形成组织的物质的浓度等等。物质的浓度指氧饱和度、氧合血红蛋白浓度、脱氧血红蛋白浓度、总血红蛋白浓度等等。总血红蛋白浓度指氧合血红蛋白浓度与脱氧血红蛋白浓度的和。此外,被检体信息可以是分布数据,诸如光吸收系数分布、氧饱和度分布等。
(基本构造)
以下参考图1描述根据本示例性实施例的光声装置的基本构造。
图1是示出根据本示例性实施例的光声装置的构造的示意图。根据本示例性实施例的光声装置包括光照射单元110、声波接收单元130、心电图取得单元150、输入单元170、显示单元180以及处理单元190。光照射单元110包括光源111和光学系统113。其详细构造将在下面进行描述。
首先,来自光源111的脉冲光112被光学系统113引导。从光学系统113中发射的脉冲光112被照射到被检体120上并且到达被检体120中的光吸收体121。光吸收体121通常是活体内的血管,尤其是诸如在血管中存在的血红蛋白、肿瘤等物质。光吸收体121吸收光的能量以产生光声波122。产生的光声波122在被检体中传播并且到达声波接收单元130。
声波接收单元130接收光声波122以输出时间序列接收信号。从声波接收单元130输出的接收信号被顺序地输入到处理单元190。上述步骤在多次光照射中的每一次光照射中执行以取得对应于这多次光照射的多个时间序列接收信号。
处理单元190利用这多个输入的时间序列接收信号产生目标区域的被检体信息。然后,处理单元190将产生的被检体信息数据发送到显示单元180,以使显示单元180显示目标区域的被检体信息的图像和/或数值。目标区域可以被预先设置或者由用户经由输入单元170输入。目标区域被设置为使得包括被检体120的至少一部分。被检体信息取得方法的细节将在下面进行描述。
同时,在其中光吸收体是血红蛋白的情况下,由于血红蛋白的量在其中血管中的血液量少的区域中小,因此该区域的光吸收系数相对低。因此,根据公式(1),产生的光声波的声压低。换句话说,在其中血液量少的区域中产生的光声波的接收信号的S/N比相对低。此外,在其中血液量非常少的情况下,光声波的接收信号会被淹没在噪声中。因此,当光声装置取得目标区域的被检体信息时,其中血液量少的区域取得的被检体信息的精度可能低。
鉴于上述问题,根据本示例性实施例的光声装置包括被配置为取得被检体120的心电图信号的心电图取得单元150。基于由心电图取得单元150取得的心电图信号的波形,可以估计被检体120的心脏的状态,并且可以相应地估计被检体120的血流状态。因此,处理单元190在不使用对应于多次光照射的多个接收信号中的当目标区域中的血液量少时产生的光声波的接收信号的情况下,基于被检体120的心电图信号来取得目标区域的被检体信息。换句话说,处理单元190利用对应于多次光照射的多个接收信号中的当目标区域中的血液量多时产生的光声波的接收信号的至少一部分来取得目标区域的被检体信息。在本示例性实施例中,心电图取得单元150对应于血液信息取得单元。
被用于取得被检体信息的信号以以下方式被提取,即,使得被检体信息可以利用更多的当血液量多时(即当作为光吸收体的血红蛋白的量多时)产生的具有高S/N比的光声波的接收信号来取得的方式。此外,根据本示例性实施例,可以在不使用当血液量少时产生的具有低S/N比的光声波的接收信号的情况下取得被检体信息。这允许具有高精度的对目标区域的被检体信息的取得。信号提取定时的细节将在下面进行描述。
以下描述根据本示例性实施例的光声装置的每个部件块。
(光源111)
光源111优选地是能够产生纳秒或微秒级脉冲光的脉冲光源。具体而言,脉冲宽度优选地为大约1至100纳秒。而且,波长优选地在大约400nm至大约1600nm的范围内。尤其在高分辨率成像活生物体的表面附近的血管的情况下,光的波长优选地是在可见光的范围内(400nm至700nm,其中包括400nm和700nm)。另一方面,在成像活生物体的深的部位的情况下,优选地使用具有较不可能被活生物体的背景组织吸收的波长(700nm至1100nm,其中包括700nm和1100nm)的光。但是,也可以使用太赫兹波、微波和无线电波的范围。
具体而言,光源111优选地是激光器。而且,在利用多个波长的光测量的情况下,能够发射可变波长的激光器是更优选的。在其中要用多个波长的光照射被检体120的情况下,可以通过切换发射光的激光器或者通过使激光器交替地发射光,使用发射彼此具有不同波长的光的多个激光器。在利用多个激光器的情况下,这多个激光器被统称为光源。
可以使用各种激光器,诸如固态激光器、气体激光器、染料激光器和半导体激光器。尤其诸如Nd:YAG激光器和变石(alexandrite)激光器的脉冲激光器是优选的。此外,可以使用利用Nd:YAG激光器光作为激发光的Ti:sa激光器,或者光参量振荡器(OPO)激光器。此外,可以使用发光二极管等代替激光器。
(光学系统113)
光学系统113将脉冲光112从光源111发送到被检体120。诸如透镜、反射镜和光纤的光学元件可以被用作光学系统113。此外,根据本示例性实施例的光学系统113包括用于改变脉冲光112的行进方向的光学反射镜114、光调整单元115和扩散板116。
在其被检体是乳房等的生物信息取得装置中,光学系统113的光发射单元优选地发射光束直径通过扩散板116等而变宽的脉冲光。另一方面,在光声显微镜中,光学系统113的光发射单元优选地包括透镜等并且以聚焦的光束直径照射光以提高分辨率。
此外,光学系统113可以包括能够调整从光源111发射的脉冲光112的衰减量的光调整单元115。诸如机械快门和液晶快门的能够调整脉冲光112的衰减量的任何单元都可以被用作光调整单元115。
此外,光学系统113可以相对于被检体120被移动以允许成像被检体120的宽广范围。
此外,光源111可以直接将光照射到被检体120上而无需利用光学系统113。
(被检体120)
以下描述被检体120,但是被检体120并不构成根据本示例性实施例的光声装置的一部分。使用根据本示例性实施例的光声装置的主要目的是人或动物的恶性肿瘤、血管疾病等的诊断,化疗的随访观察,等等。因此,被检体120被假定为活生物体,并且特定的诊断目标区域被假定为人体或动物的乳房、颈部、腹部等。
此外,被检体120中的光吸收体121优选地是在被检体120中具有相对高光吸收系数的光吸收体。例如,在其中测量目标是人体的情况下,光吸收体121可以是氧合血红蛋白、脱氧血红蛋白、含有大量氧合血红蛋白或脱氧血红蛋白的血管、或者在肿瘤附近形成的新生血管。
(声波接收单元130)
声波接收单元130包括一个或多个转换元件和外壳。可以使用能够接收声波并将声波转换为电信号的任何转换元件。这种转换元件的例子包括利用压电现象的压电元件,诸如锆钛酸铅(PZT),利用光的共振的转换元件,以及诸如电容式微加工超声换能器(CMUT)的静电电容转换元件。在其中声波接收单元130包括多个转换元件的情况下,这多个转换元件优选地被排列在平坦或弯曲的表面上,被称为1D阵列、1.5D阵列、1.75D阵列、2D阵列等等。
此外,为了取得大范围的被检体信息,声波接收单元130优选地被配置为通过扫描机构(未示出)相对于被检体120被机械地移动。此外,光学系统113(脉冲光112的照射位置)和声波接收单元130优选地彼此同步进行移动。
此外,在其中声波接收单元130是手持式声波接收单元的情况下,声波接收单元130包括用户利用其把持声波接收单元130的把持单元。此外,可以在声波接收单元130的接收表面上提供声透镜。此外,声波接收单元130可以包括多个转换元件。
此外,声波接收单元130可以包括被配置为放大从转换元件输出的时间序列模拟信号的放大器。
(心电图取得单元150)
心电图取得单元150取得被检体120的心电图信号。通常,心电图取得单元150包括用于提取心电图信号的感应电极、放大器、模数(A/D)转换器等。例如,在日本专利申请特许公开号2014-128455或2014-100244中讨论的装置可以被用作心电图取得单元150。基于由心电图取得单元150取得的心电图信号,可以估计被检体120的心脏的状态。此外,基于根据心电图信号估计的心脏的状态,可以估计血管中的血流。
(输入单元170)
输入单元170从用户(主要是诸如医疗人员的检查者)接收各种类型的输入,并且将输入信息经由系统总线发送到诸如处理单元190的部件。例如,利用输入单元170,用户可以设置与成像有关的参数、输入成像开始指令、设置诸如目标区域的范围和形状的观察参数以及进行与图像有关的其它图像处理操作。
输入单元170包括鼠标、键盘、触摸面板等,并且根据用户操作向诸如在控制单元193上运行的操作系统(OS)的软件执行事件通知。此外,在手持式光声装置的情况下,手持式光声装置优选地包括用于输入光照射单元110的驱动指令的输入单元170。作为这种输入单元170,可以采用被提供给探头的按钮开关、脚踏开关等。
(显示单元180)
显示单元180可以是显示器,诸如液晶显示器(LCD)、阴极射线管(CRT)和有机电致发光(EL)显示器。替代地,显示单元180不是被包括在根据本示例性实施例的光声装置中,而是可以被准备为单独的设备并且连接到光声装置。
(处理单元190)
作为计算机的处理单元190包括计算单元191、存储单元192和控制单元193。
计算单元191收集从声波接收单元130输出的时间序列接收模拟信号,并且执行信号处理,诸如放大接收信号、AD转换接收模拟信号以及存储数字化的接收信号。一般地,被称为“数据取得系统(DAS)”的电路可以被用作配置为执行上述处理的计算单元191。具体而言,计算单元191包括被配置为放大接收信号的放大器、被配置为数字化接收模拟信号的AD转换器等。
此外,计算单元191可以利用接收信号取得在被检体中的每个位置处的产生的声压信息。在被检体中的每个位置处的产生的声压信息也被称为被检体中的初始声压分布。在其中光声装置是光声层析成像(tomography)装置的情况下,计算单元191利用取得的接收信号执行图像重构,以取得对应于在二维或三维空间坐标上的位置的产生的声压数据。计算单元191可以使用公知的图像重构方法作为图像重构方法,诸如通用反投影(universal back projection,UBP)、滤波反投影(filtered back projection,FBP)和模型库方法。此外,计算单元191可以使用延迟与求和(delay-and-sum)处理作为图像重构方法。
此外,计算单元191可以相对于时间变化向取得的接收信号应用包络检测,将在时间轴方向上的各个光脉冲的包络检测信号的幅度值转换为在转换元件的深度方向上的幅度值,并且在空间坐标的指向方向(通常,深度方向)上绘制转换的幅度值。计算单元191为转换元件的每个位置执行上述处理,从而获得初始声压分布数据。使用上述方法是优选的,尤其在其中光声装置是光声显微镜的情况下。
诸如中央处理单元(CPU)和图形处理单元(GPU)的处理器或者诸如现场可编程门阵列(FPGA)芯片的计算电路可以被用作配置为执行处理以取得产生的声压信息的计算单元191。计算单元191可以包括单个处理器或计算电路,或者多个处理器或计算电路。
存储单元192可以存储已经历AD转换的接收信号、各种类型的分布数据、显示图像数据、各种类型的测量参数等等。此外,要在以下描述的被检体信息取得方法中执行的每个处理可以被存储在存储单元192中作为要被处理单元190中的控制单元193执行的程序。其中要存储程序的存储单元192是非暂时性记录介质。存储单元192通常是存储介质,诸如先入先出(FIFO)存储器、只读存储器(ROM)、随机存取存储器(RAM)和硬盘。存储单元192可以包括单个存储介质或多个存储介质。
此外,处理单元190包括用于控制光声装置的每个部件块的操作的控制单元193。控制单元193经由总线向光声装置的每个部件块提供必要的控制信号和数据。具体而言,控制单元193提供用于指示光源111来发射光的光发射控制信号、用于声波接收单元130中的转换元件的接收控制信号等等。控制单元193通常是CPU。
处理单元190的部件可以被集成到单个设备中或者可以是单独的设备。此外,计算单元191和控制单元193可以被包括在单个设备中。换句话说,处理单元190可以包括被配置为执行计算单元191和控制单元193的功能的单个设备。
[被检体信息取得方法]
以下参考图2描述由根据本示例性实施例的光声装置执行的被检体信息取得的流程。控制单元193读取其中描述了被检体信息取得方法的并且存储在存储单元192中的程序,并且使光声装置执行以下被检体信息取得方法。
(步骤S100:取得通过多次光照射产生的光声波的接收信号的步骤)
在步骤S100中,光照射单元110利用脉冲光112照射被检体120。然后,声波接收单元130接收通过脉冲光112的照射产生的光声波122,并且输出时间序列接收模拟信号。计算单元191收集从声波接收单元130输出的时间序列接收模拟信号,并且对接收信号执行放大处理和对接收模拟信号执行AD转换处理。然后,计算单元191将数字化的接收信号存储到存储单元192中。存储在存储单元192中的时间序列接收信号数据也被称为光声数据。在本公开内容中,术语“接收信号”的范围包括模拟信号和数字信号两者。
此外,在步骤S100中,光照射单元110执行多次光照射,使得对应于这多次光照射的多个时间序列接收信号被存储在存储单元192中。
在其中光源111是利用容易产生热的灯激发的固态激光器的情况下,为了实现光源111的稳定驱动,优选地以恒定的重复频率发射光并且在被检体120上执行多次光照射。图3A示出了根据本示例性实施例的光源111的驱动顺序。如在图3A中所示出的,在本示例性实施例中,光源111以预定的重复频率(大约5Hz)发射光。
(步骤S200:取得在多次光照射期间的心电图信号的步骤)
在步骤S200中,心电图取得单元150取得被检体120的心电图信号,并且将心电图信号发送到处理单元190。在心电图取得单元150中包括的电极被适当地布置,使得与心脏有关的肌电图信号(心电图信号)可以被取得。
图3B示出了由心电图取得单元150取得的心电图信号的典型波形。在图3B中示出的心电图信号具有大约1.2秒的周期的波形。通常,心电图信号的波形由P波、Q波、R波、S波和T波的组合形成。一般地,从R波的顶点到T波的结束附近的时段对应于心室收缩阶段,并且血液被泵送到动脉。此外,从T波的结束附近到R波的顶点的时段对应于心室扩张阶段。如本文所使用的,从R波的顶点到T波的结束附近的时间t1将被称为“从R波的产生定时(generation timing)到T波的产生定时的第一时间”。
图3C是示出在目标区域中的血液量的变化的图。如上所述,关于由心电图信号的R波触发的心脏的状态,开始心室收缩,并且相应地血液开始被泵送到动脉中。但是,如可以从图3C中理解的,目标区域中的血液量不是在R波的产生定时增加,而是具有从R波的产生定时到与心室收缩对应的血流到达目标区域的定时的时间t2的时间滞后。然后,在时间t2从R波的产生定时起逝去之后,其中血液量多的时段被考虑为保持达心室收缩时段的时间t1。
如本文所使用的,从R波的产生定时到与心室收缩对应的血流到达目标区域的定时的时间t2也被称为“延迟时间”。
(步骤S300:基于心电图信号提取当目标区域中的血液量多时取得的接收信号的步骤)
在步骤S300中,作为被检体信息取得单元的计算单元191基于在步骤S200中取得的心电图信号从在步骤S100中取得的对应于多次光照射的多个时间序列接收信号中提取要被用于取得被检体信息的信号。
基于由心电图取得单元150取得的心电图信号,计算单元191确定在目标区域中的血液量多的定时。然后,计算单元191从存储单元192中读取在该定时产生的光声波的接收信号。另一方面,计算单元191不从存储单元192中读取在目标区域中的血液量少的定时产生的光声波的接收信号,并且在该定时产生的接收信号不被用于取得被检体信息。
图3D示出了由计算单元191执行的信号提取顺序,并且指示计算单元191提取在“读取”时取得的接收信号。计算单元191从在步骤S100中取得的多个时间序列接收信号中提取在从延迟时间t2自心电图信号的R波的产生定时起逝去的定时到心室收缩的时间t1逝去的定时的时段中产生的光声波的接收信号。换句话说,计算单元191从存储单元192中读取通过在其中心室收缩引起的血流在目标区域中存在的时段中(即其中血液量多的时段中)执行的光照射产生的光声波的接收信号。另一方面,计算单元191不从存储单元192中读取通过在除了其中血液量多的时段之外的时段中执行的光照射产生的光声波的接收信号。
在步骤S300中提取的信号是在血液量由于心室收缩而增加的定时产生的光声波的接收信号。因此,所提取的信号包括许多具有高S/N比的信号。
由于光速比光声波的速度快得多,因此可以认为光声波在利用脉冲光112照射被检体120的定时在目标区域的各个位置处被同时产生。在本说明书中,利用脉冲光112照射被检体的定时被称为通过脉冲光112产生光声波的定时。
此外,已知通常从R波的产生定时到T波的产生定时的时间t1是大于等于0.3秒并且小于等于0.45秒。因此,计算单元191可以从存储单元192中读取在其中血液量多的时段中(即,在时间t2自R波产生定时起逝去之后大于等于0.3秒且小于等于0.45秒的预定时间中)产生的光声波的接收信号。
同时,计算单元191可以从心电图信号中检测产生诸如R波和T波的特定波的定时。例如,计算单元191可以把具有比预定幅度大的幅度的心电图信号的波检测为R波。此外,例如,计算单元191可以利用存储在存储单元192中的R波和T波的模板波形执行心电图信号的模板匹配,并且具有高相似性的波可以被检测为R波或T波。用于检测特定波的方法可以是任何方法,只要诸如R波和T波的特征波形可以被检测到。
由于R波已被产生,因此在对应于将心脏和目标区域之间的血管长度除以血流速度获得的值的时间逝去之后,对应于心室收缩的血流到达目标区域。因此,计算单元191可以基于R波的产生定时、关于心脏和目标区域之间的血管长度的信息和关于血流速度的信息来确定要用于取得被检体信息的信号的提取开始定时。但是,在上述方式中,为了确定要使用的信号的提取开始定时,必须为每个被检体测量心脏与目标区域之间的距离和血流速度。这会导致装置尺寸的增加。
因此,优选地从为对应于被检体的不同部位的不同区域预先确定的那些提取开始定时中选择提取开始定时。更具体而言,存储单元192优选地包括指示出目标区域的类型和延迟时间t2之间的关系的关系表。此外,光声装置优选地包括被配置为允许用户输入目标区域的类型的输入单元170。例如,输入单元170可以被配置为允许用户从在显示单元180上显示的多种类型或被检体的部位中选择目标区域的类型。然后,计算单元191可以从存储在存储单元192中的关系表中读出对应于经由输入单元170输入的类型的延迟时间t2。计算单元191从心电图信号中检测R波的产生定时,并且可以从在自检测到的产生定时起从存储单元192中读取的延迟时间t2逝去之后产生的光声波的接收信号中提取期望的信号。
虽然以上描述了目标区域的类型作为确定延迟时间t2所必需的信息,但是确定延迟时间t2所必需的信息不限于目标区域的类型。例如,可以认为,即使目标区域的类型是相同的,延迟时间t2也可取决于例如被检体的年龄而不同。因此,除目标区域的类型之外,输入单元170优选地被配置为还允许输入诸如被检体的年龄等的信息。换句话说,输入单元170优选地被配置为允许输入至少目标区域的类型。此外,控制单元193优选地从关系表中读取对应于诸如被检体年龄的输入信息的延迟时间t2。
此外,在其中光声装置的目标区域是预定的情况下,存储单元192优选地存储关于预先获得的延迟时间t2的信息。然后,计算单元191可以从心电图信号中检测R波的产生定时,并且从在自检测到的产生定时起存储在存储单元192中的延迟时间t2逝去之后产生的光声波的接收信号中提取期望的信号。
在其中从R波的产生定时到对应于心室收缩的血流到达目标区域的定时的时段可以被忽略的情况下,R波的产生定时可以被用作提取开始定时。换句话说,在这种情况下,延迟时间t2可以为t2=0。
在本示例性实施例中,接收信号提取定时是基于假设血液量在自R波的产生定时起延迟时间t2逝去之后的心室收缩期阶段的时间t1期间增加来设置的。但是,提取定时设置不限于此。例如,定时t1对应于心室收缩阶段的时间,并且取决于储存的血液量,可能存在其中大部分血液的泵送是在时间t1逝去之前完成的情况。换句话说,可能存在其中心室收缩阶段的时间与泵送血液所需的时间不匹配的情况。在这种情况下,血液量只可能在比时间t1短的时段增加。在这种情况下,处理单元190优选地利用在自R波的产生定时起延迟时间t2逝去之后在时间t1逝去之前产生的光声波的接收信号的至少一部分取得被检体信息。更具体而言,处理单元190优选地使用在自R波的产生定时起延迟时间t2逝去之后在时间t1逝去之前取得的接收信号中的、在时间t1的一半逝去之前取得的多个接收信号。
此外,例如,即使在自R波的产生定时起延迟时间t2逝去之后的时间t1期间,在血液量的增加不足够多的定时,要被取得的接收信号也可能不具有足够高的S/N比。因此,计算单元191优选地提取在自R波的产生定时起延迟时间t2逝去之后的时间t1期间具有比预定值大的幅度的接收信号。以这种方式,具有特别高S/N比的接收信号可以被选择性地在其中血液量增加的时段中提取。
在图3D中示出的序列可以与图3A至3C中示出的序列实时地并行执行,或者可以在图3A至3C中示出的序列在所有时段中被完成之后执行。
在本示例性实施例中,期望的信号是从存储在存储单元192中的多个时间序列接收信号中提取的。但是,这不是限制性的例子,而是任何方法都可以被使用,只要被检体信息可以通过选择性地利用期望的信号来取得。例如,在从声波接收单元130中输出的模拟电信号中,对应于当血液量少时产生的光声波的接收信号的模拟电信号可以从要存储在存储单元192中的信号中排除。因此,当血液量多时产生的光声波的接收信号被选择性地存储在存储单元192中。然后,计算单元191可以通过选择性地利用存储在存储单元192中的当血液量多时产生的光声波的接收信号来取得被检体信息。
(步骤S400:基于提取的接收信号取得目标区域的被检体信息的步骤)
在步骤S400中,计算单元191基于在步骤S300中提取的接收信号取得目标区域的被检体信息。在本示例性实施例中,计算单元191计算在目标区域中的每个位置处的光声波的产生的声压信息(即,初始声压分布)作为被检体信息,并且将该被检体信息存储到存储单元192中。
由于在步骤S400中取得的初始声压分布是基于在步骤S300中提取的具有高S/N比的信号计算的,因此精度高。因此,如果计算单元191使显示单元180显示存储在存储单元192中的初始声压分布的图像,则具有高的诸如分辨率和对比度的图像质量的图像可以被提供给用户。
计算单元191可以计算光注量,即,已到达目标区域中的每个位置的脉冲光112的光量分布。在本示例性实施例中,计算单元191可以通过求解在Bin Luo and Sailing He,Optics Express,Vol.15,Issue 10,pp.59055918(2007)中讨论的光扩散方程来取得关于在目标区域中的脉冲光112的光量分布的信息,并且可以将取得的信息存储在存储单元192中。计算单元191可以利用任何方法取得光量分布,只要在目标区域中的光量分布可以被取得。
然后,计算单元191可以利用在存储单元192中存储的目标区域中的初始声压分布和光量分布根据公式(1)取得在目标区域中的光吸收系数分布作为被检体信息。
在步骤S400,计算单元191可以取得在步骤S300中提取的信号中的通过单脉冲光照射取得的时间序列接收信号中的一个帧的被检体信息。此外,计算单元191可以取得在步骤S300中提取的信号中的通过多次光照射取得的多个时间序列接收信号中的一个帧的被检体信息。换句话说,计算单元191只需要利用当血液量多时产生的光声波的接收信号的至少一部分来取得被检体信息。
根据上述被检体信息取得方法,可以抑制血管中血液量对取得的被检体信息的精度的影响。
此外,根据本示例性实施例的光声装置可以利用具有不同波长的光执行上述步骤来类似地取得光吸收系数分布。然后,计算单元191可以利用对应于具有彼此不同波长的光的多个光吸收系数分布取得在被检体120中包括的物质的浓度分布信息作为被检体信息。
但是,在其中使用单个光源来产生多个波长的光的情况下,可能需要时间来切换波长。如果波长是在血液量多时被切换的,那么当血液量多时可以被执行的光照射的次数降低,从而导致被检体信息的精度较低。因此,优选地在血液量少时切换波长。例如,光照射单元110在心电图信号的一个周期期间利用第一波长入1的光照射被检体120,其中一个周期是从R波到下一个R波。然后,在其中血液量在该周期中少的时段中,光源111中的波长改变机构被驱动,使得光源111准备好产生第二波长λ2的光。然后,在下一周期中,光照射单元110利用第二波长λ2的光照射被检体120。
以这种方式,波长可以在其中血液量少并且其接收信号被确定为不被使用的时段中进行切换。这使得当血液量多时利用多个波长的光高效地照射被检体120,而不会降低在其中血液量多并且其接收信号被确定为要使用的时段中可以被执行的光照射的次数。此外,在本示例性实施例中,可以有效地确保可以被提取用于取得被检体信息的信号,使得可以有效地增加被检体信息取得的精度。
如上所述,在本示例性实施例中,被检体信息是基于多个接收信号取得的,这多个接收信号是在与被检体的脉搏相关联地被重复多次的血液量的每个波动周期中的共同时段中取得的。被检体信息是从同一定时的接收信号中取得的。这可以抑制血液量的变化对要取得的被检体信息的影响。尤其,在本示例性实施例中,使用在其中血液量在波动周期中多的时段中取得的接收信号,使得可以取得具有高精度的被检体信息的图像。
以下描述第二示例性实施例。
在第一示例性实施例中,被检体信息只根据在其中血液量多的时段中取得的接收信号产生。但是,用户可能希望在其中血液量少的时段中观察被检体区域。因此,在本示例性实施例中,将给出对以下情况的描述,即,其中取得在不被用于第一示例性实施例的被检体信息取得的其中血液量少的时段中取得的接收信号并且基于取得的接收信号取得被检体信息的情况。
根据本示例性实施例的光声装置具有与第一示例性实施例中所描述的类似的构造。
以下参考图4描述由根据本示例性实施例的光声装置取得被检体信息的流程。步骤S100和S200类似于在第一示例性实施例中的那些步骤。
(步骤S500:基于心电图信号提取当目标区域中的血液量多时取得的接收信号作为第一接收信号和在其他定时取得的接收信号作为第二接收信号的步骤)
在步骤S500中,计算单元191作为被检体信息取得单元基于在步骤S200中取得的心电图信号,从在步骤S100中取得的对应于多次光照射的多个时间序列接收信号中提取要用于取得被检体信息的信号。
计算单元191基于由心电图取得单元150取得的心电图信号估计其中在目标区域中的血液量多的时段。然后,计算单元191从存储单元192中读取在该时段中产生的光声波的接收信号作为第一接收信号。另一方面,计算单元191从存储单元192中读取在目标区域中的血液量少的定时产生的光声波的接收信号作为第二接收信号。换句话说,计算单元191提取在图3D中的“读取”时取得的接收信号作为第一接收信号,并且提取在“读取”之外的定时取得的接收信号作为第二接收信号。以这种方式提取为不同类型的接收信号的每一种被称为接收信号分组。
在步骤S500中提取的信号将在血液量由于心室收缩而增加的定时产生的光声波和在血液量由于心室扩张而降低的定时产生的光声波读作不同的接收信号。因此,在产生信号的各个定时的血液量的波动在每个接收信号分组中小。
(步骤S600:基于提取的第一和第二接收信号来取得和显示目标区域的被检体信息的步骤)
在步骤S600中,计算单元191基于在步骤S500中提取的第一和第二接收信号中的每一个取得目标区域的被检体信息。在本示例性实施例中,计算单元191计算在目标区域中的每个位置处的光声波的产生的声压信息(即,初始声压分布)作为被检体信息,并且将该被检体信息存储在存储单元192中。
如在第一示例性实施例中所描述的,被检体信息不限于初始声压分布,并且可以是在被检体120中包含的物质的光吸收系数分布信息或浓度分布信息。
存在已知的平均从多次激光照射取得的光吸收系数分布或物质浓度分布的数值以增加这些数值的S/N比的方法。同时,光吸收系数分布和物质浓度分布与血红蛋白的量成比例,因此,可以说光吸收系数分布和物质浓度分布与血液量成比例。因此,与从与血管中的血液量无关地取得的多个接收信号中取得的被检体信息的平均值相比,具有较小血液量波动的被检体信息可以根据从第一接收信号中取得的被检体信息的平均值和从第二接收信号中取得的被检体信息的平均值取得。
以下参考附图5描述显示取得的被检体信息的方法。
区域220是其中显示从第一接收信号产生的初始声压分布的图像的区域。区域221是其中显示从第二接收信号产生的初始声压分布的图像的区域。一般地,初始声压分布被取得为三维(3D)图像。因此,至于初始声压分布的显示,可以显示3D图像,或者可以显示3D图像的截面图,某个范围的最大强度投影(MIP)图像,等等。
区域200是这样的区域,其中显示接收信号在其中被提取为第一接收信号的范围和接收信号在其中被提取为第二接收信号的范围。范围204是指示在第一示例性实施例中描述的时间t2的范围。范围205是指示在第一示例性实施例中描述的时间t1的范围。范围201指示从P波的开始到自R波的顶点起时间t2逝去的定时的时段,并且是其中接收信号被提取为第二接收信号的范围。范围202指示时间t1,并且是其中接收信号被提取为第一接收信号的范围。范围203指示从时间t1的结束到P波的开始的范围,并且是其中接收信号被提取为第二接收信号的范围。典型的心电图波形可以被显示在区域200的背景中以使得用户容易理解。此外,可以显示由心电图取得单元150实际取得的最新心电图波形。典型的心电图波形和最新的心电图波形可以被并排(alongside)显示在区域200中。
项210是利用其可以输入/改变时间t1的用户界面(UI)部分。项211是利用其可以输入/改变时间t2的UI部分。范围205和204的宽度可以结合利用项210和211做出的改变进行改变。此外,步骤S500的第一和第二接收信号提取处理基于结合项210和211中的变化而最新设置的时间t1和t2被再次执行。此外,被检体信息被再次取得,并且要被显示在区域220和221上的图像被更新。在图5所示的例子中,如果用户选择被指示为项210的上部或下部三角形,则在文本框中指示的时间t1的长度被更新,并且限定时间t1的虚线之间的距离在区域200上被改变。用户可以通过直接将数字输入到文本框中或者通过利用鼠标等拖动区域200中的虚线等来设置时间t1的范围,从而确定时间t1的开始时间和结束时间。时间t2可以利用类似于用于时间t1的方法的方法来设置。
要被显示的被检体信息不限于初始声压分布并且可以是光吸收系数分布或物质浓度分布。
在本示例性实施例中,已经给出了以下示例情况的描述,即,其中心电图波形的一个周期被划分为两个时段,其中一个是其中被检体区域中的血液量多的时段(图5中的时段(1)),并且其中另一个是其中血液量少的时段(图5中的时段(2))。可替代地,一个周期可以被划分为三个或更多个片段。如果一个周期被划分为更多数量的片段,则在每个片段中血液量的波动变得较小。因此,可以获得更均匀信号电平的接收信号。另一方面,如果心电图波形的划分数量增加,则在每个片段中取得的接收信号的数量减少,因此通过平均接收信号产生的增加S/N比的效果减小。在确定一个周期要被划分的划分数量时考虑这种权衡是所期望的。
虽然基于第一接收信号的图像和基于第二接收信号的图像在本示例性实施例中被显示为彼此并排,但是任何其它显示方法也可以被使用。例如,两个图像可以被重叠并且显示。此外,图像可以被交替地显示。如果心电图波形的划分数量增加,则取得的图像可以像运动图像那样来显示。由于通过多次照射取得的接收信号的平均值的图像被顺序地显示,因此与其中每个初始声压分布都是从通过单次激光照射取得的接收信号中取得的初始声压分布像运动图像那样被顺序地显示的情况相比,基于具有高S/N比的图像的运动图像可以被提供给用户。
虽然在本示例性实施例中,被检体信息在光声波信号和心电图信号的取得之后被取得和显示,但是被检体信息也可以在光声波信号和心电图信号被取得的同时实时地被取得和显示。此外,在其中被检体信息被实时地取得和显示的情况下,不是利用从取得开始时取得的所有光声波信号和心电图信号,而是只有部分光声波信号和心电图信号可被使用。例如,可以考虑使用在从当前时间之前给定时间的时间到当前时间的时段期间取得的光声波信号和心电图信号。以这种方式,即使被检体被移动,也可以取得较准确的被检体信息。
上述第二示例性实施例产生了与第一示例性实施例类似的优点。此外,根据本示例性实施例中,可以向用户提供在血液量的波动周期中的各不同时段的被检体信息。动脉在其中血液量增加的定时搏动,并且认为动脉和静脉可以通过比较多条被检体信息进行彼此区分。
本发明的(一个或多个)实施例也可以由系统或装置的计算机来实现,其中计算机读出并执行记录在存储介质(也可以更完整地被称为“非临时性计算机可读存储介质”)上的计算机可执行指令(例如,一个或多个程序),以执行上述(一个或多个)实施例之中的一个或多个实施例的功能,和/或计算机包括用于执行上述(一个或多个)实施例之中的一个或多个实施例的功能的一个或多个电路(例如,专用集成电路(ASIC)),并且由系统或装置的计算机通过例如从存储介质中读出并执行计算机可执行指令以便执行上述(一个或多个)实施例之中的一个或多个实施例的功能和/或控制一个或多个电路执行上述(一个或多个)实施例之中的一个或多个实施例的功能所执行的方法来实现。计算机可以包括一个或多个处理器(例如,中央处理单元(CPU)、微处理单元(MPU))并且可以包括独立计算机或独立处理器的网络来读出和执行计算机可执行指令。计算机可执行指令可以从例如网络或存储介质中提供给计算机。存储介质可以包括例如硬盘、随机存取存储器(RAM)、只读存储器(ROM)、分布式计算系统的存储器、光盘(诸如压缩盘(CD)、数字多样化盘(DVD)或蓝光盘(BD)TM)、闪存存储器设备、存储卡等等当中的一种或多种。
其它实施例
本发明的实施例还可以通过如下的方法来实现,即,通过网络或者各种存储介质将执行上述实施例的功能的软件(程序)提供给系统或装置,该系统或装置的计算机或是中央处理单元(CPU)、微处理单元(MPU)读出并执行程序的方法。
虽然前面描述了本发明的特定示例性实施例,但是应当理解,上述特定示例性实施例不是要限制本发明的范围,并且在本发明的技术构思内可以进行修改。
例如,虽然前面描述了其中心电图仪被用作配置为确定被测目标中的血液量增加/减少的血液信息取得单元的情况,但是被测目标中的血液量可以利用用于测量被检体脉搏的任何其它方法来确定。例如,可以使用红外脉博计。
虽然本发明已经参考示例性实施例进行了描述,但是应当理解,本发明不限于所公开的示例性实施例。以下权利要求的范围是要赋予最广泛的解释,从而涵盖所有此类修改和等同结构及功能。
Claims (7)
1.一种光声装置,其特征在于,包括:
光照射单元,被配置为利用脉冲光多次照射被检体;
接收单元,被配置为接收通过利用来自光照射单元的脉冲光多次照射被检体而产生的光声波,并且输出与所述多次光照射对应的多个信号;
血液信息取得单元,被配置为取得关于被检体的血液量的信息;及
被检体信息取得单元,被配置为基于所述多个信号取得在被检体中的目标区域的被检体信息,
其中,被检体信息取得单元基于在血液量的重复的波动周期中的各波动周期内的共同时段中取得的多个信号来取得被检体信息,
并且其中,被检体信息取得单元取得基于在所述周期中的各周期内的第一时段中取得的多个信号的第一被检体信息,并且取得基于在所述周期中的各周期内的与第一时段不同的第二时段中取得的多个信号的第二被检体信息。
2.如权利要求1所述的光声装置,其中,被检体信息取得单元使显示单元将第一被检体信息和第二被检体信息彼此并排显示。
3.如权利要求2所述的光声装置,其中,被检体信息取得单元使显示单元显示第一被检体信息和第二被检体信息以使得第一被检体信息和第二被检体信息重叠。
4.如权利要求1所述的光声装置,其中,被检体信息取得单元使显示单元显示被检体信息和指示出在所述周期中的各周期内的第一时段和第二时段的信息。
5.如权利要求1所述的光声装置,其中,血液信息取得单元包括被配置为取得被检体的心电图信号的心电图取得单元,及
其中,被检体信息取得单元使显示单元显示心电图信号。
6.一种基于多个信号来取得被检体信息的被检体信息取得方法,所述多个信号与多次光照射对应并且通过接收利用脉冲光照射被检体所述多次所产生的光声波而被取得,其特征在于,该方法包括:
基于在血液量的重复的波动周期中的各波动周期内的共同时段中取得的多个信号来取得被检体信息,
其中,基于在所述周期中的各周期内的第一时段中取得的多个信号的第一被检体信息被取得,并且基于在所述周期中的各周期内的与第一时段不同的第二时段中取得的多个信号的第二被检体信息被取得。
7.如权利要求6所述的被检体信息取得方法,其中,第一被检体信息是通过从所述多个信号中选择性地利用与如下光声波对应的信号的至少一部分而被取得的,该光声波是在从基于从被检体取得的心电图信号的R波的产生定时确定的开始定时起直到从心电图信号的R波的产生定时到心电图信号的T波的产生定时的第一时间逝去为止的所述第一时段中产生的。
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2014242454 | 2014-11-28 | ||
JP2014-242454 | 2014-11-28 | ||
JP2015204123A JP6614910B2 (ja) | 2014-11-28 | 2015-10-15 | 光音響装置 |
JP2015-204123 | 2015-10-15 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105640496A CN105640496A (zh) | 2016-06-08 |
CN105640496B true CN105640496B (zh) | 2019-09-03 |
Family
ID=56121431
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510846561.8A Active CN105640496B (zh) | 2014-11-28 | 2015-11-27 | 光声装置和被检体信息取得方法 |
Country Status (2)
Country | Link |
---|---|
JP (1) | JP6614910B2 (zh) |
CN (1) | CN105640496B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6594355B2 (ja) | 2017-01-06 | 2019-10-23 | キヤノン株式会社 | 被検体情報処理装置および画像の表示方法 |
JP2019033806A (ja) | 2017-08-10 | 2019-03-07 | キヤノン株式会社 | 光音響装置および被検体情報取得方法 |
CN111387969B (zh) * | 2020-03-26 | 2023-09-19 | 苏州市东方电子仪器厂 | 食管心脏电生理食管导联同步记录技术 |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4549552A (en) * | 1981-03-06 | 1985-10-29 | Siemens Gammasonics, Inc. | Heart sound detector and cardiac cycle data are combined for diagnostic reliability |
JP3713329B2 (ja) * | 1996-06-04 | 2005-11-09 | 株式会社東芝 | 超音波ドプラ診断装置 |
GB0126804D0 (en) * | 2001-11-07 | 2002-01-02 | Univ London | Flow velocity measurement |
JP2005342006A (ja) * | 2004-05-31 | 2005-12-15 | Toshiba Corp | 超音波診断装置、超音波画像処理装置、及び超音波信号処理プログラム |
EP2219514A1 (en) * | 2007-11-14 | 2010-08-25 | Koninklijke Philips Electronics N.V. | Systems and methods for detecting flow and enhancing snr performance in photoacoustic imaging applications |
JP5619347B2 (ja) * | 2007-12-20 | 2014-11-05 | 東芝メディカルシステムズ株式会社 | 超音波診断装置、超音波画像解析装置及び超音波画像解析プログラム |
JP2010172376A (ja) * | 2009-01-27 | 2010-08-12 | Toshiba Corp | 超音波画像診断装置および画像処理プログラム |
JP2010269018A (ja) * | 2009-05-22 | 2010-12-02 | Shimadzu Corp | 超音波診断装置 |
US9241643B2 (en) * | 2012-05-31 | 2016-01-26 | Covidien Lp | Methods and systems for power optimization in a medical device |
EP2754388B1 (en) * | 2013-01-15 | 2020-09-09 | Helmholtz Zentrum München Deutsches Forschungszentrum für Gesundheit und Umwelt GmbH | System and method for quality-enhanced high-rate optoacoustic imaging of an object |
JP2016042922A (ja) * | 2014-08-20 | 2016-04-04 | プレキシオン株式会社 | 光音響画像化装置 |
-
2015
- 2015-10-15 JP JP2015204123A patent/JP6614910B2/ja active Active
- 2015-11-27 CN CN201510846561.8A patent/CN105640496B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN105640496A (zh) | 2016-06-08 |
JP6614910B2 (ja) | 2019-12-04 |
JP2016107069A (ja) | 2016-06-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6335612B2 (ja) | 光音響装置、処理装置、処理方法、及びプログラム | |
JP6498036B2 (ja) | 光音響装置、信号処理方法、及びプログラム | |
EP3266378A1 (en) | Apparatus, method, and program for obtaining information derived from ultrasonic waves and photoacoustic waves | |
WO2018008439A1 (en) | Apparatus, method and program for displaying ultrasound image and photoacoustic image | |
CN105640496B (zh) | 光声装置和被检体信息取得方法 | |
CN105640494B (zh) | 光声装置、被检体信息获取方法和程序 | |
JP6742745B2 (ja) | 情報取得装置および表示方法 | |
JP6656229B2 (ja) | 光音響装置 | |
US20180228377A1 (en) | Object information acquiring apparatus and display method | |
US20160150990A1 (en) | Photoacoustic apparatus, subject information acquisition method, and program | |
JP2023123874A (ja) | 光音響イメージングシステム、光音響イメージングシステムの制御方法、および、プログラム | |
EP3329843B1 (en) | Display control apparatus, display control method, and program | |
CN108403083A (zh) | 信息处理装置、信息处理方法和非暂时性存储介质 | |
JP7108985B2 (ja) | 画像処理装置、画像処理方法、プログラム | |
US20190099083A1 (en) | Information processing apparatus and information processing method | |
US10617319B2 (en) | Photoacoustic apparatus | |
US20210169397A1 (en) | Image processing apparatus, image processing method, and non-transitory computer-readable medium | |
JP2019000387A (ja) | 情報処理装置、情報処理方法、及びプログラム | |
CN105640499B (zh) | 光声装置 | |
JP6929204B2 (ja) | 情報処理装置、情報処理方法、およびプログラム | |
US20190321005A1 (en) | Subject information acquisition apparatus, subject information processing method, and storage medium using probe to receive acoustic wave | |
JP2020162746A (ja) | 画像処理装置、画像処理方法及びプログラム | |
JP2018161467A (ja) | 画像処理装置および画像処理方法 | |
JP2020028661A (ja) | 画像処理装置、画像処理方法、プログラム | |
JP2020028671A (ja) | 画像処理装置、画像処理方法、プログラム |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |