CN103222880A - 被检物信息获取装置 - Google Patents
被检物信息获取装置 Download PDFInfo
- Publication number
- CN103222880A CN103222880A CN2013100314019A CN201310031401A CN103222880A CN 103222880 A CN103222880 A CN 103222880A CN 2013100314019 A CN2013100314019 A CN 2013100314019A CN 201310031401 A CN201310031401 A CN 201310031401A CN 103222880 A CN103222880 A CN 103222880A
- Authority
- CN
- China
- Prior art keywords
- light
- checking matter
- light quantity
- sound wave
- photo
- 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
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0059—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
-
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/1702—Systems in which incident light is modified in accordance with the properties of the material investigated with opto-acoustic detection, e.g. for gases or analysing solids
-
- 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
- A61B5/7225—Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Public Health (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Biophysics (AREA)
- Veterinary Medicine (AREA)
- Acoustics & Sound (AREA)
- Analytical Chemistry (AREA)
- Chemical & Material Sciences (AREA)
- Biochemistry (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Signal Processing (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Power Engineering (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physiology (AREA)
- Psychiatry (AREA)
Abstract
本发明公开了被检物信息获取装置。本发明利用如下这样的被检物信息获取装置,该被检物信息获取装置包括:多个光照射器,其将来自光源的光照射到被检物;强度调制器,其对于每个光照射器调制照射光的强度;多个声检测器,其检测从所述被检物产生的声波;和信息处理器,其通过使用所述声波来生成所述被检物的内部特性信息,其中,信息处理器基于每个光照射器中的光的强度来控制生成所述特性信息时每个声检测器的权重。
Description
技术领域
本发明涉及一种被检物信息获取装置。
背景技术
关于光学成像技术的研究已被积极地开展,该光学成像技术通过将来自光源(诸如激光器)的光照射到活体以在体内(in vivo)传播光并检测传播的光等来获得体内信息。已被提出的这些光学成像技术中的一种技术是光声层析成像(PAT)(M.Xu and L.V.Wang″Photoacoustic imaging in biomedicine″,Review of ScientificInstruments,77,041101(2006))。光声层析成像是如下这样一种技术,在该技术中,将从光源产生的脉冲光照射到活体,在多个点处检测从吸收了在体内传播并扩散的光能的体组织产生的声波,并对其信号进行分析和处理以使体内信息可视化。通过该技术,可获取光学特性的值的体内分布,尤其是光能吸收密度的分布。这样的体内成像技术也被称为图像重构。
根据M.Xu and L.V.Wang″Photoacoustic imaging in biomedicine″,Review of ScientificInstruments,77,041101(2006),在光声层析成像中通过光吸收而从体内的吸收体获得的声波(光声波)的声压(P)可用下式(1)表示:
P=Γ·μa·Φ...(1)
其中,Γ是Gruneisen系数,该系数是弹性的特性值,并且是通过将体积膨胀系数(β)和声速(c)的平方的乘积除以比热(Cp)来获得;μa是吸收体的吸收系数,并且Φ是局部区域中的光量(照射到吸收体的光量)。
已知,当组织被识别时Gruneisen系数Γ几乎恒定。因此,通过在多个点处执行对应于声波的程度的声压P的变化的时间解析(time-resolved)测量,可获得吸收系数μa和光量Φ的乘积,即,光能吸收密度的分布。
日本专利申请特开No.2009-018153公开了基于光声层析成像的装置的例子。当对基于光声波的信号进行分析以确定光学特性的值的分布时,该装置利用体内的被照射区域与光吸收体之间的相对位置信息和声波的声压。因此,试图获得具有高分辨率的光学特性的值的体内分布(特别地,光吸收系数或体内平均等效衰减系数的分布)并且使它们准确地可视化。
日本专利申请特开No.2010-088627公开了基于光声层析成像的装置的另一个例子。该装置利用基于活体的形状而确定的光量的体内分布,以便从基于光声波的电信号获得体内信息。
专利文献1:日本专利申请特开No.2009-018153
专利文献2:日本专利申请特开No.2010-088627
非专利文献1:M.Xu and L.V.Wang″Photoacoustic imaging inbiomedicine″,Review of Scientific Instruments,77,041101(2006)
发明内容
在光声层析成像中,优选的是降低运算过程期间的计算的噪声,以便通过图像重构来获得合适的图像。因此,希望的是确保目标中的均匀光量分布。通过确保体内的均匀光量分布,在可照射预定范围内的可能的最强光的同时,可使每单位时间的入射光量限于预定值或更低以减少对活体的伤害。
为了确保体内的均匀光照射,优选的是从光源、从多个方向、在多个照射部位照射光。然而,另一方面,当这些照射部位中的任何一个位于声波检测器附近时,所检测的声波值可能越来越多地包含噪声。这是由声波检测器的电路中的噪声或者来自声波检测器的组件本身的光声效应引起的。作为这样的现象的结果,已知的是重构结果也包含噪声,并且图像中可能包含伪像。
鉴于前述描述,本发明的一个目的是提供一种可降低光声层析成像中的噪声的技术。
本发明提供一种被检物信息获取装置,该被检物信息获取装置包括:
多个光照射器,被配置为将来自光源的光照射到被检物;
强度调制器,被配置为对所述多个光照射器中的每一个的照射光的强度进行调制;
多个声检测器,被配置为检测从被所述光照射的被检物产生的声波;以及
信息处理器,被配置为通过使用所述声波来生成所述被检物的内部特性信息,
其中,信息处理器基于所述多个光照射器中的每一个中的、已被强度调制器调制的光的强度来控制生成所述特性信息时所述多个声检测器中的每一个的权重。
根据本发明,可提供可降低光声层析成像中的噪声的技术。
从以下参照附图对示例性实施例的描述,本发明的进一步的特征将会变得清楚。
附图说明
图1是显示根据第一实施例的生物信息处理装置的配置的平面图;
图2是显示根据第一实施例的生物信息处理装置的配置的侧视图;
图3是根据第一实施例的光量分布的重构的概念图;
图4是根据第一实施例的处理的流程图;和
图5是根据第一实施例的光量强度调制的时序图。
具体实施方式
以下参照附图来描述本发明的优选实施例。然而,以下描述的组件的尺寸、材料和形状及其相对位置等应该根据应用本发明的配置和各种条件来进行适当修改,因此,以下描述并非意图限制本发明的范围。
本文中所描述的生物信息处理装置用于恶性肿瘤或血管疾病的诊断、化疗的跟踪评估等的目的。该生物信息处理装置获取光学特性的值的体内分布和由其获得的特性信息(诸如构成物质的生物组织的浓度分布、被检物中的初始声压分布、从初始声压分布导出的光能吸收密度分布、吸收系数分布等)。物质的浓度分布是指例如氧饱和度的分布、氧合血红蛋白浓度和脱氧血红蛋白浓度的分布等。基于所述特性信息,所述装置可创建被检物内部的图像,即,可执行图像重构。也就是,本发明的生物信息处理装置也可被称为被检物信息获取装置。
本文中所使用的声波典型地指的是超声波,并且包括被称为音波、超声波和声波的弹性波。通过光声效应产生的声波被称为光声波或光超声波。声波检测器(比如,探针)接收在被检物中产生或反射的声波。
<第一实施例>
图1示出了根据本发明的第一实施例的生物信息处理装置的配置。基于图1来描述本发明的第一实施例。
(装置的配置)
生物信息处理装置包括光源102、光耦合器103、声波检测器106、光量计107、重构单元108、显示单元109、控制单元110、强度调制器111和固定构件112。在本实施例中,将被测量的被检物(目标)是活体110。
光源102发射光101,并且在该图中显示了四个光源(102-1至102-4)。
光源产生具有使得被构成活体的成分之中的特定成分吸收的特定波长的光。光源包括可产生几纳秒至几百纳秒的量级的脉冲光的至少一个脉冲光源。光源优选地为激光器,或者可替换地,它可以是发光二极管。可使用各种激光器,诸如固态激光器、气体激光器、染料激光器、半导体激光器等。优选使用的光具有700nm或更大和1100nm或更小的近红外区中的波长,该波长在体内被较少地吸收。然而,当试图确定相当接近于活体表面的生物组织的光学特性的值的分布时,还可使用更宽的波长范围,比如,400nm或更大和1600nm或更小的波长范围。当光源102是半导体激光器时,可通过电流或电压控制来方便地进行强度调制。
光耦合器103是被形成有例如透镜、反射镜、光纤等的光学***。光耦合器103将从光源102发射的、然后被强度调制器111进行强度调制的光作为光101经由固定构件112引入到活体100中。在该图中,与光源的数量一致地示出了四个光耦合器(103-1至103-4)。光耦合器对应于本发明的光照射器。
当在活体100中传播的部分光能被吸收到光吸收体104(诸如血管)中时,光吸收体104产生声波105。声波检测器106检测从光吸收体104产生的声波105,并将该声波信号转换为电信号。声波检测器对应于本发明的声检测器(acoustic detector)。
固定构件112将活体100固定在生物信息处理装置中。固定构件112由传播光(特别是近红外光和声波)的材料形成。具体例子可包括透明树脂和透明玻璃。当在活体100与固定构件112之间存在间隙时,提供填料113。填料113是传播光(特别是近红外光和声波)的材料,并且可以是匹配液(诸如水)。为了确保在活体100与固定构件112之间没有间隙,还可用减压装置(未显示)来降低活体100和固定构件112的压力,或者改变固定构件112的形状。
光量计107是用于测量来自活体100的光量数据的光检测器。光量计107测量如下的光的强度,所述光从光耦合器照射,在体内传播,并在体组织中扩散和吸收之后再次到达活体的表面。在本实施例中,提供多个光量计107(在该图中,四个),并且获得多个光量数据作为测量结果。在控制单元110中,基于从测量得到的光量数据来估计活体100中的光量分布。可通过同时从多个光源102照射光并在每个光量计处测量光量来估计光量分布。还可通过从所述多个光源之一照射光、估计由该光源导致的光量分布并重复这些步骤多次来估计光量分布。
显示单元109是显示由重构单元108重构的图像信息的设备。在本发明的生物信息处理装置中,显示单元109不是必要组件。与生物信息处理装置连接的另一显示设备可被用作显示单元109。
控制单元110控制声波检测器106、光量计107和强度调制器111。可设想控制单元110的多种控制方式。本发明的信息处理器对应于重构单元和控制单元,并且可用包括处理器并遵循程序操作的信息处理器等来实现。
(处理流程)
图4显示了本发明的示意性流程图。在该图中,步骤S100至S103是准备过程,并且在步骤S104至S107中,进行实际的光声测量和体内成像。
首先,将活体100固定在固定构件112上(步骤S100)。
接着,对每个光耦合器103测量光量分布(步骤S101)。以下详细描述具体过程。简要地讲,在这个过程中,对当体内成像时像素中的从每个光耦合器照射的光的近似光量进行估计。
接着,确定照射光量的组合Pm(步骤S102)。也就是说,确定从每个光耦合器照射到活体的实际光量。
然后,确定用于每个Pm的重构权重系数Am(步骤S103)。
其后,启动实际数据收集。首先,进行具有光量Pm的光的照射、从被检物传播的声波的接收和基于该声波的图像重构(步骤S104)。在这种情况下,强度调制器控制每个光耦合器以便照射具有光量Pm的光。在重构时,应用上述权重系数Am。因此,计算重构图像Im。
然后,根据光量分布的估计来校正重构图像Im。具体地讲,通过将重构图像Im除以估计的光量分布来校正光量分布(步骤S105)。
在改变光量Pm的同时,将步骤S104和S105重复预定次数M(步骤S106)。
最后,对校正后的重构图像Im求取平均值以获得最终的重构图像(步骤S107)。
(权重控制)
现在详细地研究以上各个过程。
为了解决本发明的问题(包含由于声波检测器附近的光照射而导致的噪声),可采取简单的方式,在所述简单的方式中,在照射光的光耦合器附近的声波检测器处获得的声学数据(acoustic data)不被用于图像重构。可替换地,可采取另一种方式,在所述另一种方式中,降低来自光耦合器附近的声检测器的声学数据对重构图像的贡献程度。
首先,描述根本不使用来自光耦合器附近的声检测器的声学数据的方式。在图1中,当从光耦合器103-1照射光时,来自声波检测器106-1和106-8的、可能包含噪声的声学数据不被用于图像重构(在表1中,T=1)。如此获得的重构图像1具有降低的噪声。
另一方面,当从光耦合器103-3照射光时,来自声波检测器106-4和106-5附近的声学数据不被用于图像重构(在表1中,T=3)。类似地,在从光耦合器103-2和103-4照射光的情况下,不使用来自附近的声检测器的声学数据。
表1中总结了以上控制方式。
[表1]
计算重构图像1至重构图像4的算术均值,以生成最终的重构图像。通过如以上那样计算算术均值,对于图像的具有低光量和低SN比的部分,多个重构图像可相互补偿,从而提供优选的图像。
接下来,描述降低来自光耦合器附近的声检测器的声学数据对重构图像的贡献程度的方式。在这种情况下,如表2所示,针对各个重构图像,对在各声波检测器处获得的声学数据进行加权。权重意指用于重构图像的算术加法的权重系数。
在表1中所示的例子中,可认为,用于重构的声波检测器的权重为1,而不用于重构的声波检测器的权重为0。例如,在T=1的情况下,可认为,声波检测器106-1和106-8的权重为0,而声波检测器106-2至106-7的权重为1。
还可以例如在0与1之间加权。
表2显示了这样的加权方式。用于确定权重的标准之一是从照射光的光耦合器到每个声波检测器的距离。也就是说,减小用于离光耦合器的距离短的声波检测器的权重,这是因为这些检测器趋向于包含噪声(在T=1的情况下,106-1和106-8)。另一标准是来自光耦合器的光的行进距离。也就是说,离光耦合器的距离增大的部分接收较少量的光,因此,用于这样的位置的权重减小(在T=1的情况下,106-4和106-5)。因此,当将多个重构图像相加时,这些图像相互补偿。
[表2]
可仅考虑从照射光的光耦合器到每个声波检测器的距离来确定特定的权重系数。权重系数还可被确定为使得它与从光耦合器到每个声波检测器的距离成比例。可替换地,权重系数还可通过以下假设被确定:来自照射光的光耦合器的噪声根据该距离的平方根定律或指数定律而衰减。
还可同时从多个光耦合器照射光。例如,以表3中所示的方式,从相对的光耦合器照射光。考虑从照射光的光耦合器到每个声波检测器的距离来确定权重系数。在表3中,接近照射光的光耦合器的声波检测器的权重系数为0。然而,它可以是不为零的一个小的值。在这种情况下,到达接近照射光的光耦合器的声波检测器的数据可在一定程度上被反应在重构中。
[表3]
(体内光量分布的估计)
现在描述本发明要解决的另一个问题,体内均匀光量的获得。这种情况下的光量的估计在实际光声测量(S101)之前;然而,它不是总是被执行。作为从实际测量值估计光量的可替换方案,可基于体内的标准光传播来确定光量。本发明的降噪效果可仅通过如上所述那样执行加权来获得。
首先,选择多个光耦合器103之一,从该光耦合器照射光,并估计该光耦合器导致的体内光量分布。对于所述多个光耦合器,重复光量估计。
当进行实际光声测量时,可同时从多个光耦合器或者从单独的光耦合器照射光。在前一种情况下,通过组合对于各个光耦合器估计的光量,可确定当同时从多个光耦合器照射光时的体内光量分布。从每个光耦合器照射的光量可由强度调制器111控制。因此,通过基于估计值来控制强度调制器,还可对体内光量进行校正。
从所选的光耦合器照射的光在体内传播,以到达每个光量计。照射光量与到达每个光量计的光成比例。然而,当光量太低时,由于用作光量计的光检测器本身的固有噪声,它不能被检测到。在这种情况下,强度调制器改变光的强度,以便照射强度等于或高于光量计的固有噪声的光。光的衰减取决于各个活体的性质。因此,为了确保体内的均匀光分布,优选的是根据每个活体来调制强度。然而,还可使用标准衰减比或者通过相同类型的活体初步测量的值。
现在描述控制单元110的旨在确保体内的均匀光量的控制方式。可通过将从各个光耦合器照射的光的体内光量分布相加来确定体内光量分布。为了执行这个操作,从光耦合器中被选的一个光耦合器将光照射到活体。然后在每个光量计107,对值进行测量。在本实施例中,光量计设置在四个不同的位置(107-1至107-4)。每个光量计是一维或两维光检测器,因此,每个光量计107在多个点处对光量进行测量。
已知在体内传播的光量如以下式(2)所示那样在体内呈指数地衰减:
Φ=Φ0·exp(-μ·d1)...(2)
其中,Φ0是入射光量,μ是从光入射位置到光量计的体内平均衰减系数,并且d1是从光入射点到光量计的距离。
控制单元基于各个光量计处的测量值来估计体内光量分布。通过反投影到三维存储器、内插和低通滤波来估计光量分布。三维存储器可通过控制单元(未显示)中的存储器形成。
图2是图1的侧视图。包括固定构件112的立方体的体积对应于边长为20cm的立方体。当边长为20cm的立方体被解析为200×200×200个像素时,三维存储器的像素的一边为1mm。照射光的光耦合器的位置和光量计的位置在三维存储器上是已知的。控制单元然后将入射光量Φ0、所测量的光量Φ和从式(2)(反投影)计算的光量分配给连接光量计和光耦合器的线上的像素。对所有光量计107执行反投影。因此,这条线上的像素的光量被估计。
在本实施例中,固定构件的所有外表面并未都被光量计覆盖。因此,即使在对所述线上的像素进行反投影之后,仍剩有还未对其估计光量的像素。对于这样的像素可通过内插来确定光量。在本实施例中,在各种内插算法中,使用线性内插。
参照图3来描述线性内插的过程。在模拟活体的立方体的表面上,提供光耦合器103和光量计107。中心附近的像素X的光量还未被估计。从光量计107反投影到光耦合器103的线被称为光束Rn(n=1…N,N为光量计中所包括的光检测器的数量)。假设当从光耦合器103照射光束时光束的数量与构成光量计107的光检测器的数量相同。
从像素X到光耦合器103的距离被指定为L。选择具有最短的从像素X下降到光束Rn的垂直线的长度Pn的两个光束Ri和Rj。光束Ri和Rj在距离L处的光量被称为Ri(L)和Rj(L)。然后,像素X的光量XL可用式(3)表示:
XL=(Ri(L)+Rj(L))/2...(3)
根据式(3)对未被分配光量的所有像素进行计算。然而,一些像素可能需要通过外插来方便地确定Rn(L)。当希望提高光量估计的精度时,可根据从光耦合器到像素的线与从光耦合器到光检测器的线之间的角度来适当地修改该式。
然后,在内插之后,对三维存储器执行三维低通滤波。式(2)描述了光在衰减的同时在直线上行进的情况。然而,光在体内散射,并且散射光也到达光量计。光散射可被认为是降低对比度的低通滤波。因此,对三维存储器执行三维低通滤波。可用的低通滤波器的例子是相对大的内核(kernel),诸如11×11×11像素。通过执行低通滤波,完成对一个光耦合器的光量分布的估计。对所有光耦合器重复相同过程,将各个光耦合器中的各个像素的光量加在一起,使得当同时从多个光源照射光时的光量分布可被确定。如此估计的光量还可被用于强度调制器处的校正或者图像重构时的校正。
(光量校正)
现在描述用于当同时从多个光源照射光时确保体内的均匀光量分布的方式。近似解足以确保均匀光量分布,而不必确定严格解。这是因为,如从图1和图2中的配置可以看到的,活体的中央部分的光量必然低,因此,难以确保严格的均匀性。即使当光量分布均匀时,也有必要考虑声波检测器从光耦合器接收的噪声。各种公式可被设想为评估函数,以便评估体内光量分布的均匀性。在本实施例中,利用在三维存储器中的某些位置(多于一个的位置)处的光量与平均光量之间具有最小差值的评估函数。
当光耦合器103-1至103-4的光量分别被指定为P1、P2、P3和P4时,三维存储器中的特定坐标(X,Y,Z)处的光量Bxyz用式(4)表示:
Bxyz=Fxyz(P1,P2,P3,P4)...(4)
其中,Fxyz是确定坐标(X,Y,Z)处的光量的函数,并且是P1、P2、P3和P4的线性函数。通过每个光耦合器103上的光量分布的估计来确定该线性函数的系数。
通过使用式(4),活体100中的平均光量Ba用式(5)表示:
Ba=1/N·∑Fi(P1,P2,P3,P4)...(5)
其中,i=1…N,N是活体100内的坐标。
评估函数H用式(6)表示:
H=∑(Ba-Fj(P1,P2,P3,P4))2...(6)
其中,j=1…K,表示活体100内的特定坐标。所述特定坐标被选择为使得该特定坐标均匀地分布在活体100内。K被确定为使得匹配光耦合器的数量。在本实施例中,存在四个光耦合器。因此,通过限定K=4,可确定使评估函数H最小的P1、P2、P3和P4。
式(7)是使P1、P2、P3和P4最小的联立方程:
H/θPj=0,(j=1..K)(7)
对于特定坐标的每一组合(在本实施例中,四个元素),可确定满足式(7)的P1、P2、P3和P4。
选择特定坐标的多个组合。当各个组合被指定为Sm(m=1..M)时,对于每个Sm,确定来自四个光耦合器的输出强度Pm=(Pm1,Pm2,Pm3,Pm4)。
对于每个组合Sm确定光量(P1至P4),并根据所确定的光量来执行光声测量。通过一次光声测量来生成重构图像Im,并对Im(m=1..M)求取平均以确定最终的重构图像。因此,根据每个重构图像Im的SN比和最终重构图像的目标SN比来确定M。在本实施例中,M为100。当M=100时,可获得为每个重构图像Im的大约10倍的SN比。
在确定Pm(m=1…M)之后,通过考虑来自光耦合器的噪声的量来确定声波检测器的权重系数Am。在本实施例中,使用八个声波检测器,因此,声波检测器的权重系数被表示为Am=(Am1..Am8)。在本实施例中,通过以下假设来确定Ami(i=1..8):由于光耦合器的存在而导致的声波检测器的噪声的量根据距离的平方根定律而减小。当从光耦合器103-i到声波检测器106-j的距离被指定为Lij时,Aj用公式(8)表示:
Aj=C·∑Pi·(Lij2)...(8)
其中,j=1…8。C是系数,并被确定为使得∑Aj=1。在确定Pm和Am之后,以时间序列的方式驱动Pm(m=1..M),并生成根据所确定的Am的重构图像Im。
(处理的时序图)
图5显示了时序图,在该时序图中,以光量分布S101和各种光量Pm来照射光来对重构图像Im进行重构。
该时序图的左侧对应于实际光声测量之前的准备阶段,并对应于上述流程图的S101至S103。各个光耦合器103-1至103-4照射具有所确定的强度的光,并且光量计处的测量和控制单元处的光量分布的估计被执行。在每个光耦合器的测量之后,确定权重系数Am和光量的组合Pm。
该时序图的右侧对应于实际光声测量,并对应于流程图的S104及其后续步骤。每个光耦合器根据光量Pm和M次的光照射来进行控制,并且伴随照射的声波的检测和基于这些声波的图像重构被执行。对来自各个周期的所有重构图像求取平均,以生成最终的重构图像。
<第二实施例>
当光源是激光光源时,有必要理解用以确保活体的安全性的照射的光量或输出。描述本发明的第二实施例(用于减轻光对活体的伤害的方法)。
上述第一实施例的生物信息处理装置使得可基于测量的光量来估计体内光量分布。该值对应于光输出,即,施加于活体的能量的量。因此,光量分布的估计使得可计算光对活体的伤害。
控制单元基于从光量Pm预测的光量分布的估计的控制使得可将每单位时间的入射光量限制为预定值或更低。所述预定值包括,例如,符合安全标准的值。通过对光量求平均,可在符合安全标准的同时尽可能地增大光量,从而使得可生成好得多的重构图像。
如上所述,本发明使得可控制多个光源和多个声波检测器来使每个光量和每个声波检测器的用于重构的权重最优化。它还可确保被检物中的大致均匀的光量分布。结果,可降低重构图像中的噪声。
另外,通过合适地控制光量,还可减轻对目标的伤害。
因此,本发明可适用于利用光声效应的技术,特别是人类和动物的体内信息成像的领域。
尽管已参照示例性实施例描述了本发明,但是要理解本发明不限于所公开的示例性实施例。要给予权利要求的范围以最宽泛的解释,以便包含所有这样的修改以及等同的结构和功能。
Claims (6)
1.一种被检物信息获取装置,包括:
多个光照射器,被配置为将来自光源的光照射到被检物;
强度调制器,被配置为对于所述多个光照射器中的每一个的照射光的强度进行调制;
多个声检测器,被配置为检测从被所述光照射的被检物产生的声波;以及
信息处理器,被配置为通过使用所述声波来生成所述被检物的内部特性信息,
其中,所述信息处理器基于所述多个光照射器中的每一个中的、已被所述强度调制器调制的光的强度来控制生成所述特性信息时所述多个声检测器中的每一个的权重。
2.根据权利要求1所述的被检物信息获取装置,其中,在生成所述特性信息时,所述信息处理器减小关于由所述多个声检测器检测的声波之中的、由接近于照射了光的光照射器的声检测器检测的声波的权重。
3.根据权利要求1或2所述的被检物信息获取装置,其中,所述信息处理器基于所述被检物中的光量分布来控制生成所述特性信息时用于所述多个声检测器各自的权重。
4.根据权利要求3所述的被检物信息获取装置,还包括光量计,所述光量计被配置为测量从所述光照射器照射的并在所述被检物中传播的光的强度,
其中,从通过所述光量计获得的测量结果估计所述被检物中的光量分布。
5.根据权利要求3所述的被检物信息获取装置,其中,所述强度调制器通过基于所述被检物中的光量分布对所述多个光照射器中的每个的光的强度进行调制,来减小所述被检物中的光量变化。
6.根据权利要求3所述的被检物信息获取装置,其中,所述强度调制器基于所述被检物中的光量分布来对于所述多个光照射器中的每一个的光的强度进行调制,以将每单位时间照射到所述被检物的光量限制为预定值或更小。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2012-018305 | 2012-01-31 | ||
JP2012018305A JP5988598B2 (ja) | 2012-01-31 | 2012-01-31 | 被検体情報取得装置および被検体情報取得方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103222880A true CN103222880A (zh) | 2013-07-31 |
CN103222880B CN103222880B (zh) | 2015-05-20 |
Family
ID=47713838
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310031401.9A Active CN103222880B (zh) | 2012-01-31 | 2013-01-28 | 被检物信息获取装置 |
Country Status (5)
Country | Link |
---|---|
US (1) | US9259158B2 (zh) |
EP (1) | EP2623022A1 (zh) |
JP (1) | JP5988598B2 (zh) |
KR (1) | KR20130088777A (zh) |
CN (1) | CN103222880B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105640493A (zh) * | 2014-11-28 | 2016-06-08 | 佳能株式会社 | 光声装置和用于光声装置的控制方法 |
CN106175677A (zh) * | 2016-07-08 | 2016-12-07 | 华南师范大学 | 基于光纤分束和柔性探测器的一体化光声乳腺成像检测装置和方法 |
CN106667433A (zh) * | 2015-11-10 | 2017-05-17 | 佳能株式会社 | 被检体信息获取装置 |
CN108882881A (zh) * | 2016-03-25 | 2018-11-23 | 索尼公司 | 图像分析设备和图像分析方法 |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2013128661A (ja) | 2011-12-21 | 2013-07-04 | Canon Inc | ステレオx線撮影装置、ステレオx線撮影方法 |
JP6108671B2 (ja) | 2012-03-13 | 2017-04-05 | キヤノン株式会社 | 放射線撮影装置 |
JP2015109948A (ja) * | 2013-10-31 | 2015-06-18 | キヤノン株式会社 | 被検体情報取得装置 |
JP6545190B2 (ja) * | 2014-05-14 | 2019-07-17 | キヤノン株式会社 | 光音響装置、信号処理装置、信号処理方法、プログラム |
JP6456129B2 (ja) * | 2014-12-15 | 2019-01-23 | キヤノン株式会社 | 被検体情報取得装置およびその制御方法ならびに光量制御方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5491329A (en) * | 1993-12-03 | 1996-02-13 | Hamamatsu Photonics K.K. | Photodetecting apparatus having intensity tuneable light irradiating unit |
US20100070233A1 (en) * | 2008-09-17 | 2010-03-18 | Canon Kabushiki Kaisha | Measurement apparatus |
CN101677765A (zh) * | 2007-05-15 | 2010-03-24 | 佳能株式会社 | 生物信息成像设备和方法、生物信息分析方法 |
US20100087733A1 (en) * | 2008-10-07 | 2010-04-08 | Canon Kabushiki Kaisha | Biological information processing apparatus and biological information processing method |
US20110083509A1 (en) * | 2009-10-09 | 2011-04-14 | Nellcor Puritan Bennett Llc | Photoacoustic Spectroscopy With Focused Light |
WO2011093069A1 (en) * | 2010-01-27 | 2011-08-04 | Canon Kabushiki Kaisha | Subject information acquisition apparatus and subject information acquisition method |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4589783A (en) * | 1984-04-04 | 1986-05-20 | Wayne State University | Thermal wave imaging apparatus |
JP4549093B2 (ja) | 2004-04-12 | 2010-09-22 | キヤノン株式会社 | 画像処理装置及びその方法、プログラム |
JP4469903B2 (ja) | 2007-06-11 | 2010-06-02 | キヤノン株式会社 | 生体情報イメージング装置 |
WO2010024290A1 (ja) | 2008-08-27 | 2010-03-04 | キヤノン株式会社 | 生体の光音響情報処理装置および生体の光音響情報処理方法 |
JP5538855B2 (ja) * | 2009-12-11 | 2014-07-02 | キヤノン株式会社 | 光音響装置及び該装置の制御方法 |
SE1050086A1 (sv) | 2010-01-27 | 2011-06-28 | St Jude Medical Systems Ab | Sensorstyrtrådsanordning och system för intravaskulära mätningar av en fysiologisk variabel |
JP5693043B2 (ja) * | 2010-04-28 | 2015-04-01 | キヤノン株式会社 | 被検体情報取得装置、被検体情報取得方法 |
-
2012
- 2012-01-31 JP JP2012018305A patent/JP5988598B2/ja not_active Expired - Fee Related
- 2012-12-19 US US13/719,609 patent/US9259158B2/en not_active Expired - Fee Related
-
2013
- 2013-01-14 EP EP13151163.6A patent/EP2623022A1/en not_active Withdrawn
- 2013-01-28 KR KR1020130009032A patent/KR20130088777A/ko not_active Application Discontinuation
- 2013-01-28 CN CN201310031401.9A patent/CN103222880B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5491329A (en) * | 1993-12-03 | 1996-02-13 | Hamamatsu Photonics K.K. | Photodetecting apparatus having intensity tuneable light irradiating unit |
CN101677765A (zh) * | 2007-05-15 | 2010-03-24 | 佳能株式会社 | 生物信息成像设备和方法、生物信息分析方法 |
US20100070233A1 (en) * | 2008-09-17 | 2010-03-18 | Canon Kabushiki Kaisha | Measurement apparatus |
US20100087733A1 (en) * | 2008-10-07 | 2010-04-08 | Canon Kabushiki Kaisha | Biological information processing apparatus and biological information processing method |
US20110083509A1 (en) * | 2009-10-09 | 2011-04-14 | Nellcor Puritan Bennett Llc | Photoacoustic Spectroscopy With Focused Light |
WO2011093069A1 (en) * | 2010-01-27 | 2011-08-04 | Canon Kabushiki Kaisha | Subject information acquisition apparatus and subject information acquisition method |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105640493A (zh) * | 2014-11-28 | 2016-06-08 | 佳能株式会社 | 光声装置和用于光声装置的控制方法 |
CN106667433A (zh) * | 2015-11-10 | 2017-05-17 | 佳能株式会社 | 被检体信息获取装置 |
CN108882881A (zh) * | 2016-03-25 | 2018-11-23 | 索尼公司 | 图像分析设备和图像分析方法 |
CN106175677A (zh) * | 2016-07-08 | 2016-12-07 | 华南师范大学 | 基于光纤分束和柔性探测器的一体化光声乳腺成像检测装置和方法 |
CN106175677B (zh) * | 2016-07-08 | 2018-12-07 | 华南师范大学 | 基于光纤分束和柔性探测器的一体化光声乳腺成像检测装置和方法 |
Also Published As
Publication number | Publication date |
---|---|
US20130197342A1 (en) | 2013-08-01 |
JP2013154045A (ja) | 2013-08-15 |
KR20130088777A (ko) | 2013-08-08 |
US9259158B2 (en) | 2016-02-16 |
JP5988598B2 (ja) | 2016-09-07 |
EP2623022A1 (en) | 2013-08-07 |
CN103222880B (zh) | 2015-05-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103222880B (zh) | 被检物信息获取装置 | |
US20190350460A1 (en) | Photoacoustic imaging apparatus, photoacoustic imaging method, and program for executing photoacoustic imaging method | |
US8289502B2 (en) | Measurement apparatus and measurement method | |
US20170215739A1 (en) | Object information acquiring apparatus and processing method | |
CN102131463B (zh) | 生物信息处理设备和生物信息处理方法 | |
EP2826417B1 (en) | Object information acquiring apparatus and method for controlling object information acquiring apparatus | |
US9750413B2 (en) | Massively parallel diffuse optical tomography | |
JP5773578B2 (ja) | 被検体情報取得装置、被検体情報取得装置の制御方法およびプログラム | |
US20120302866A1 (en) | Photoacoustic imaging apparatus and photoacoustic imaging method | |
CN103300880A (zh) | 被检体信息获得装置和被检体信息获得方法 | |
CN103300881A (zh) | 对象信息获得装置和对象信息获得方法 | |
WO2014045564A1 (en) | Object information acquiring apparatus and control method for the object information acquiring apparatus | |
EP3178380A1 (en) | Photoacoustic apparatus, display control method, and program | |
CN105078408A (zh) | 光声装置 | |
US20160058295A1 (en) | Photoacoustic wave measurement apparatus and photoacoustic wave measurement method | |
CN105686800B (zh) | 被检体信息获取装置及其控制方法 | |
CN106455994A (zh) | 光声装置 | |
JP2862016B2 (ja) | 被検体通過光中の散乱光成分抑制方法およびその装置 | |
Pal et al. | Transient radiation modeling of short-pulse laser detection of tumors in animal model | |
Ren et al. | Non-contact scanning diffuse optical tomography for three-dimensional vascular imaging in a murine bone graft model | |
JPH0454439A (ja) | 生体計測方法および装置 | |
Ahbbinaya et al. | Diffuse Optical Tomographic imaging using NIR Camera | |
Nouizi et al. | 3D noncontact time-resolved Fluorescent Diffuse Optical Tomography data processing for improving image’s quality |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |