CN102238903A - 光声成像设备和光声成像方法 - Google Patents
光声成像设备和光声成像方法 Download PDFInfo
- Publication number
- CN102238903A CN102238903A CN200980148696.XA CN200980148696A CN102238903A CN 102238903 A CN102238903 A CN 102238903A CN 200980148696 A CN200980148696 A CN 200980148696A CN 102238903 A CN102238903 A CN 102238903A
- Authority
- CN
- China
- Prior art keywords
- light
- absorption coefficient
- image
- sound wave
- absorber
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 61
- 239000006096 absorbing agent Substances 0.000 claims abstract description 50
- 238000012545 processing Methods 0.000 claims abstract description 44
- 230000008859 change Effects 0.000 claims abstract description 12
- 238000010521 absorption reaction Methods 0.000 claims description 128
- 238000001514 detection method Methods 0.000 claims description 21
- 230000002093 peripheral effect Effects 0.000 claims description 7
- 230000004044 response Effects 0.000 claims description 5
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 238000000034 method Methods 0.000 abstract description 17
- 230000003287 optical effect Effects 0.000 abstract description 15
- 239000000463 material Substances 0.000 description 10
- 230000015572 biosynthetic process Effects 0.000 description 9
- 238000005259 measurement Methods 0.000 description 9
- 239000000126 substance Substances 0.000 description 8
- 210000004204 blood vessel Anatomy 0.000 description 3
- 239000003086 colorant Substances 0.000 description 3
- 238000012937 correction Methods 0.000 description 3
- 238000001727 in vivo Methods 0.000 description 3
- 238000003672 processing method Methods 0.000 description 3
- OMOVVBIIQSXZSZ-UHFFFAOYSA-N [6-(4-acetyloxy-5,9a-dimethyl-2,7-dioxo-4,5a,6,9-tetrahydro-3h-pyrano[3,4-b]oxepin-5-yl)-5-formyloxy-3-(furan-3-yl)-3a-methyl-7-methylidene-1a,2,3,4,5,6-hexahydroindeno[1,7a-b]oxiren-4-yl] 2-hydroxy-3-methylpentanoate Chemical compound CC12C(OC(=O)C(O)C(C)CC)C(OC=O)C(C3(C)C(CC(=O)OC4(C)COC(=O)CC43)OC(C)=O)C(=C)C32OC3CC1C=1C=COC=1 OMOVVBIIQSXZSZ-UHFFFAOYSA-N 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000000630 rising effect Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 206010059245 Angiopathy Diseases 0.000 description 1
- 102000008186 Collagen Human genes 0.000 description 1
- 108010035532 Collagen Proteins 0.000 description 1
- WQZGKKKJIJFFOK-GASJEMHNSA-N Glucose Natural products OC[C@H]1OC(O)[C@H](O)[C@@H](O)[C@@H]1O WQZGKKKJIJFFOK-GASJEMHNSA-N 0.000 description 1
- 102000001554 Hemoglobins Human genes 0.000 description 1
- 108010054147 Hemoglobins Proteins 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000002308 calcification Effects 0.000 description 1
- 201000011510 cancer Diseases 0.000 description 1
- 230000000973 chemotherapeutic effect Effects 0.000 description 1
- 239000008103 glucose Substances 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 230000031700 light absorption Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000033116 oxidation-reduction process Effects 0.000 description 1
- NRNCYVBFPDDJNE-UHFFFAOYSA-N pemoline Chemical compound O1C(N)=NC(=O)C1C1=CC=CC=C1 NRNCYVBFPDDJNE-UHFFFAOYSA-N 0.000 description 1
- 238000010895 photoacoustic effect Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000010189 synthetic method Methods 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
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/0059—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
- A61B5/0073—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence by tomography, i.e. reconstruction of 3D images from 2D projections
-
- 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
-
- 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
- G01N2021/1706—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 in solids
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Pathology (AREA)
- General Health & Medical Sciences (AREA)
- Animal Behavior & Ethology (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Biophysics (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Acoustics & Sound (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
一种光声成像设备执行光吸收体的成像。该光声成像设备包括:光源;检测器,所述检测器被配置为检测从已经吸收从所述光源发射的光的能量的光吸收体产生的声波;以及信号处理单元,所述信号处理单元被配置为形成光吸收体的图像。所述信号处理单元在对由检测器检测的声波执行波形处理之前存储指示由检测器检测的声波的压力的变化率是正还是负的信息。
Description
技术领域
本发明涉及使用光声效应的成像设备以及光声成像方法。
背景技术
近年来,已经提出了用于通过使用与光相比在活体内散射更少的超声波的特性来以高分辨率获得在活体内的光学特性值的分布的方法。
根据PTL,用由光源产生的脉冲光照射活体,检测由生物组织通过脉冲光的能量吸收而产生的声波,并且分析与所检测的声波对应的信号,由此获得在活体内的光学特性值的分布。使用通过用光照射活体而获得的声波的成像一般被称为光声成像。
关于光声成像方法,已知以下技术。即,在如图9A中所示出地由声波检测器20检测通过用光照射球形光吸收体10而产生的光声波时,如果光吸收体10的光吸收是均匀的,则获得图9B中示出的N形(N-shape)的声压信息(参见NPL)。
通过将N形的波形的时间宽度乘以音速获得的值是反映了光吸收体10的尺寸(这里为球的直径)的值。此外,检测到N形的波形时的时间反映了光吸收体10的位置信息。此外,在到达光吸收体10的光强度相等的情况下,具有N形波形的信号的幅值与光吸收体10的吸收系数成比例。
如上所述,在光声成像中,通过使用从光声波获得的数据来重建光吸收体10的图像。
在通过检测通过光学吸收产生的超声波而执行的上述光声成像中,对具有比光吸收体周围的介质的光吸收系数高的光吸收系数的组织成像。例如,活体中的血管具有比周围介质的光吸收系数高的光吸收系数。已经研究了血管的成像。
作为用于处理从检测器获得的光声波信号的方法,可以使用诸如包络检波之类的波形处理。通过在波形处理之后执行图像形成处理,可以对活体内的光学特性分布成像。
在使用由设置在各个位置处的检测器测量的光声波信号时,可以通过使用图像重建和诸如包络检波之类的波形处理来对活体内的光学特性分布成像,所述图像重建使用包含诸如延迟求和(delay andsum)之类的孔径合成处理的方法。
在光声成像中使用的约700到1100nm的近红外区中,存在具有比周围组织的光吸收系数高的光吸收系数的组织(诸如血管),并且还存在具有比周围组织的光吸收系数低的光吸收系数的组织(诸如钙化的物质)。
为此,在上述的简单地使用包络检波的图像形成方法中,只可以以绝对值的形式检测对象和周围介质之间的光吸收系数的差。因此,难以确定对象的光吸收系数是比周围介质的光吸收系数高还是低。
引文列表
专利文献
PTL 1:美国专利No.5,713,356
非专利文献
NPL 1:L.V.Wang等人的“Biomedical Optics-Principles andImaging”,Wiley,Ch.12,2007
发明内容
本发明提供了能够通过将具有比周围介质的光吸收系数低的光吸收系数的组织和具有比周围介质的光吸收系数高的光吸收系数的组织彼此区分来执行成像的光声成像设备和光声成像方法。
根据本发明的方面的光声成像设备包括:光源;检测器,所述检测器被配置为检测从已经吸收从所述光源发射的光的能量的光吸收体产生的声波;以及信号处理单元,所述信号处理单元被配置为通过在对由检测器检测的声波执行波形处理之前存储指示所述声波的压力的变化率是正还是负的信息来形成光吸收体的图像。
根据本发明的方面的光声成像方法包括:从光源发射光的步骤;检测从已经吸收从所述光源发射的光的能量的光吸收体产生的声波的步骤;以及在对所检测的声波执行波形处理之前存储指示所述声波的压力的变化率是正还是负的信息的步骤。
根据本发明的方面,能够提供如下的光声成像设备和光声成像方法,所述光声成像设备和光声成像方法能够通过将具有比周围介质的光吸收系数低的光吸收系数的组织和具有比周围介质的光吸收系数高的光吸收系数的组织彼此区分来执行成像。
从结合附图进行的以下描述中本发明的其它特征和优点将明显,在附图中相似的附图标记指示在所有附图中相同或类似的部件。
附图说明
图1示出根据本发明实施例的光声成像设备的示例性配置。
图2A是用于说明具有比背景的光吸收系数高的光吸收系数的光吸收体的成像方法的图。
图2B是用于说明具有比背景的光吸收系数高的光吸收系数的光吸收体的成像方法的图。
图2C是用于说明具有比背景的光吸收系数高的光吸收系数的光吸收体的成像方法的图。
图2D是用于说明具有比背景的光吸收系数高的光吸收系数的光吸收体的成像方法的图。
图3A是用于说明具有比背景的光吸收系数低的光吸收系数的光吸收体的成像方法的图。
图3B是用于说明具有比背景的光吸收系数低的光吸收系数的光吸收体的成像方法的图。
图3C是用于说明具有比背景的光吸收系数低的光吸收系数的光吸收体的成像方法的图。
图3D是用于说明具有比背景的光吸收系数低的光吸收系数的光吸收体的成像方法的图。
图3E是用于说明具有比背景的光吸收系数低的光吸收系数的光吸收体的成像方法的图。
图4是根据第一实施例的流程图。
图5是根据第二实施例的流程图。
图6是根据第三实施例的流程图。
图7A是用于说明根据本发明的示例的模拟模型的图。
图7B是用于说明根据本发明的示例的模拟模型的图。
图8A是用于说明根据本发明的示例的模拟结果的图。
图8B是用于说明根据本发明的示例的模拟结果的图。
图9A是用于说明具有N形波形的光声信号的图。
图9B是用于说明具有N形波形的光声信号的图。
具体实施方式
现在将根据附图详细描述本发明的实施例。
图2A~2D是用于说明在球形组织210的光吸收系数比周围介质220的光吸收系数高的情况下的图像处理方法的图。如图2A中所示出的,在由检测器240检测通过用脉冲光230照射组织210而产生的光声波时,获得图2B中示出的波形数据250。在光被吸收体均匀地吸收时,波形数据是N形的。
其后执行的用于处理的方法具有两个可替代方案。
在它们中的一种方法中,如图2C中所示出的,对图2B中示出的波形数据执行基于包络检波的处理。基于处理后的波形数据执行图像形成处理,并且如图2D中所示出地在显示单元270上显示反映光吸收系数的光吸收体的图像280。
在另一种方法中,在使用孔径合成处理等的算法的图像重建方法中将在各个位置处检测的多段波形数据相加,然后执行包络检波,并且在显示单元270上显示反映光吸收系数的光吸收体的图像280。
另一方面,图3A~3E是用于说明在球形组织310的光吸收系数比周围介质220的光吸收系数低的情况下的图像处理方法的图。
本发明的发明人通过研究已经发现了,在组织310的光吸收系数比周围介质220的光吸收系数低的这种情况下,如图3B中所示出的,检测到具有相对于图2B中示出的波形正/负被颠倒的N形波形(颠倒的N形波形)的数据350。
即,在目标部分的光吸收系数比周围组织的光吸收系数高的情况下,光声波信号在信号开始点处上升(即,开始时压力的变化率是正的)。另一方面,在目标部分的光吸收系数比周围组织的光吸收系数低的情况下,光声波信号在信号开始点处下降(即,开始时压力的变化率是负的)。
已经使用上述发现来进行了本发明,并且本发明的特征在于,在将具有比周围介质的光吸收系数高的光吸收系数的组织与具有比周围介质的光吸收系数低的光吸收系数的组织区分开的情况下执行成像。
第一实施例
图1示出根据本发明第一实施例的光声成像设备100的示例性配置。
根据本实施例的光声成像设备100出于诊断恶性肿瘤、血管疾病等或者观察化学治疗的进展的目的而能够对活体内的光学特性值和构成从关于那些值的信息获得的生物组织的物质的浓度分布成像。特别地,光声成像设备100能够对具有比周围介质的光吸收系数低的光吸收系数的光吸收体成像。
在根据本实施例的光声成像设备100中,从光源130发射的光通过光纤140传播,使得用该光照射活体120。作为周围介质的活体120具有光吸收体110。由光吸收体110吸收的光的能量被转换为声波160。由检测器150检测声波160。即使在光吸收体110的光吸收系数比构成对象的周围介质的光吸收系数低的情况下,光吸收体110也根据光吸收体110与周围介质之间的光吸收系数的差而产生声波。
此外,根据本实施例的光声成像设备100包括信号处理单元170,该信号处理单元170通过使用从检测器150获得的电信号而获得关于光学特性值的分布信息。
光源130产生脉冲光。脉冲光在几纳秒到数百纳秒的量级,并且其波长应该为700nm以上并且1100nm以下。使用激光器作为光源130,但是可以使用发光二极管等来代替激光器。在使用能够转换振荡波长的OPO(光学参数振荡器)或者染料激光器时,能够测量由波长引起的光学特性值的分布的差异。
每个检测器150吸收照射活体120的光的能量,根据光吸收系数的差来检测从活体120中的光吸收体110产生的声波,并且将声波转换成电信号。作为检测器150,可以使用能够只接收从特定区域产生的声波的聚焦换能器(focus transducer)。在使用聚焦换能器时,可以指定光吸收体的位置,使得不必执行诸如延迟求和之类的图像重建的处理。此外,可以使用任何类型的声波检测器(例如,使用压电现象的换能器、使用光谐振的换能器或者使用电容变化的换能器)作为检测器150,只要该检测器能够检测声波信号即可。在本实施例中,在活体120的表面上放置多个检测器150。可替代地,可以使用单个检测器来扫描活体120的表面。
信号处理单元170分析电信号,由此获得关于活体120的光学特性值的分布信息。可以使用任何装置作为信号处理单元170,只要声波的时间变化和强度可以被存储并且通过使用计算单元被转换为光学特性值的分布的数据即可。例如,可以使用示波器和能够分析存储在示波器中的数据的计算机。
在图2A中示出的球形组织210的情况下,如图2B中所示出的,从具有比周围介质220的光吸收系数高的光吸收系数的组织210产生的光声信号在开始点处上升(N形的波形)。
另一方面,在图3A中示出的具有比周围介质220的光吸收系数低的光吸收系数的球形组织310的情况下,如图3B中所示出的,信号在开始点处下降(颠倒的N形的波形)。
将参考图4中的关于基于上述发现对具有比周围介质的光吸收系数低的光吸收系数的组织成像的信号处理方法的流程图来进行描述。
首先,用作检测器150的聚焦换能器检测光声信号(S410)。
随后,信号处理单元170存储指示光声信号的压力的变化率是具有正值还是具有负值的信息(S420)。为了使光声信号的数据适合于光吸收体的吸收特性,执行包络检波处理(例如,在希尔伯特(Hilbert)变换之后获得绝对值)作为波形处理(S430)。在图2C和图3C中示出了通过包络检波处理获得的声压信息的示例。
随后,确定在光声信号的压力变化中光声信号是否在信号开始点处上升(S440)。应该依赖于使用的检测器的特性而在确定之前根据检测器的频率响应特性来执行校正。即,应该根据使用的检测器的频率响应特性来对接收的光声信号执行去卷积。
这里,如果在S440中确定光声信号上升,则以与传统光声成像设备中相同的方式执行图像形成(S460)。另一方面,如果在步骤S440中确定光声信号下降,则包络检波后的声压数据的正/负被颠倒(S450)。利用S450中的该处理,可以获得图3D中示出的波形数据360。随后,以与在信号上升的情况中相同的方式来执行图像形成(S460)。为了执行图像形成,波形数据360的时间数据可以被转换为位置数据,并且随后被绘制(plot)。
最后,从信号处理单元170输出在S460中获得的图像的信号,由此在图像显示单元180和270上显示光吸收体的图像。在该情况下,如图2D和图3E中所示出的,可以以不同的显示形式在图像显示单元270上显示具有比周围组织的光吸收系数高的光吸收系数的组织的图像280和具有比周围组织的光吸收系数低的光吸收系数的组织的图像380。为了表示具有不同的光吸收系数的不同类型的组织,可以使用不同的显示颜色或者不同的色调(例如,改变浓淡度)。可替代地,可以把识别符号(例如,+或-)给予具有比周围组织的光吸收系数高或低的光吸收系数的组织。
根据第一实施例的光声成像设备100的使用使得即使在使用包络检波时也能够不仅对具有比活体内的周围介质的光吸收系数高的光吸收系数的组织成像而且对具有比周围介质的光吸收系数低的光吸收系数的组织成像。
已知的是,物质的光吸收系数根据光的波长而变化。在使用多个波长的光的情况下,计算对于各个波长的活体内的光吸收系数。然后,将那些系数与构成生物组织的物质(葡萄糖、胶原蛋白、氧化-还原血红蛋白等)特有的波长相关性进行比较,使得可以对构成活体的物质的浓度分布成像。
此外,在通过使用根据本实施例的光声成像设备100利用多个波长执行测量的情况下,可以确定在物质的光吸收系数与周围介质的光吸收系数之间的高/低关系,使得可以指定作为测量目标的物质。
现在,将对于确定其光吸收系数与周围介质的光吸收系数不同的物质A、B和C的存在的示例来进行描述。假设在某两个波长lambda-1和lambda-2两者中物质A的光吸收系数比周围介质的光吸收系数高。此外,假设物质B的光吸收系数在波长lambda-1中较高而在波长lambda-2中较低。此外,假设物质C的光吸收系数在波长lambda-1和lambda-2两者中较低。
在该情况下,在通过使用两个波长lambda-1和lambda-2来执行测量时,在两个波长中从物质A输出正信号(在信号开始点处上升的信号)。至于物质B,在使用波长lambda-1时获得正信号,而在使用波长lambda-2时输出负信号(在信号开始点处下降的信号)。至于物质C,在两个波长中输出负信号。通过比较这些结果,可以相互区分物质A、B和C。
如上所述,通过使用包括能够发射多个波长的光的光源和信号处理单元的光声成像设备,可以基于图像重建结果的正/负值来识别物质A、B和C或者可以容易地确定物质A、B和C的分布。
即,在通过用波长lambda-1的光的照射获得的第一声波的形状与通过用波长lambda-2的光的照射获得的第二声波的形状之间的比较使得能够分析物质。可以由上述信号处理单元170或另一单元执行这种处理。
第二实施例
根据第二实施例的光声成像设备与根据第一实施例的光声成像设备100的不同之处在于,在已经执行了通常的包络检波处理之后执行图像形成处理。即,在第一实施例中,在已经使包络检波后的声压数据的正/负颠倒之后执行图像形成。另一方面,在第二实施例中,在不执行这种处理的情况下执行图像形成。根据第二实施例的设备的配置除了信号处理和图像重建之外都与第一实施例的配置相同,并且因此省略了其描述。
图5是根据第二实施例的图像重建和信号处理的流程图。
首先,用作检测器150的聚焦换能器检测光声信号(S510)。
随后,信号处理单元170存储光声信号的压力变化信息(S520)。在已经执行了包络检波处理以便使光声信号的数据适合于光吸收体的吸收特性(S530)之后,执行图像形成处理(S540)。在图像形成处理中,波形数据360的时间数据被转换为位置数据,并且被绘制。
在根据第一实施例的图像形成处理中,在通过包络检波获得的数据的正/负已经在检测到在信号开始点处的下降时被颠倒之后,执行图像重建。然而,根据第二实施例,在没有使通过包络检波获得的数据的正/负颠倒的情况下执行图像形成处理。
随后,基于在S520中存储的光声信号的压力变化信息来确定光声信号是否下降(S550)。应该依赖于使用的检测器的特性在确定之前根据检测器的频率响应特性来执行校正。
如果在S550中确定光声信号上升,则在图像显示单元180和270上显示在S540中获得的图像(S570)。
另一方面,如果在S550中确定光声信号下降,则添加指示组织的光吸收系数比背景的光吸收系数低的信息(S560)。其后,从信号处理单元170输出添加有该信息的图像的信号,并且在图像显示单元180和270上显示该图像(S570)。
因此,可以以不同的显示形式显示具有比周围组织的光吸收系数高的光吸收系数的组织的图像和具有比周围组织的光吸收系数低的光吸收系数的组织的图像。为了表示具有不同的光吸收系数的不同类型的组织,可以使用不同的显示颜色或者不同的色调。可替代地,可以把识别符号给予具有比周围组织的光吸收系数低的光吸收系数的组织。
在根据第二实施例的光声成像设备中,与在根据第一实施例的设备中一样,可以通过使用多个波长的光来指定作为测量目标的物质。
第三实施例
根据第三实施例的光声成像设备与根据第一实施例和第二实施例的光声成像设备的不同之处在于,通过使用基于延迟求和的算法或者基于傅里叶变换的算法来执行图像重建。作为根据本实施例的检测器,应该使用能够检测来自各个区域的信号的检测器。除了上述点之外,根据第三实施例的设备的配置与根据第一实施例的设备的配置相同,并且因此省略了其描述。
图6是根据第三实施例的图像重建和信号处理的流程图。
首先,检测器150检测光声信号(S610)。
随后,信号处理单元170存储光声信号的压力变化信息(S620),对光声信号的数据执行希尔伯特变换(S630),并且执行图像重建(S640)。作为用于图像重建的方法,使用基于延迟求和的算法或者基于傅里叶变换的算法。其后,获得所获得的图像数据的绝对值(S650)。
随后,基于在S620中存储的光声信号的压力变化信息来确定光声信号是否下降(S660)。应该依赖于使用的检测器的特性在确定之前根据检测器的频率响应特性来执行校正。
如果在S660中确定光声信号上升,则从信号处理单元170输出所获得的图像的信号,并且在图像显示单元180和270上显示该图像(S680)。
另一方面,如果在S660中确定光声信号下降,则添加指示组织的光吸收系数比背景的光吸收系数低的信息(S670)。其后,从信号处理单元170输出添加有该信息的图像的信号,并且在图像显示单元180和270上显示该图像(S680)。
因此,可以以不同的显示形式显示具有比周围组织的光吸收系数高的光吸收系数的组织的图像和具有比周围组织的光吸收系数低的光吸收系数的组织的图像。为了表示具有不同的光吸收系数的不同类型的组织,可以使用不同的显示颜色或者不同的色调。可替代地,可以给予识别符号(例如,+或-)。
可替代地,可以通过使光吸收体的光吸收系数的正/负颠倒来在图像显示单元180和270上显示通过包络检波后的图像重建获得的数据。
在根据第三实施例的光声成像设备中,与在根据第一实施例的设备中一样,可以通过使用多个波长的光来指定作为测量目标的物质。
示例1
对于在下面的情况下获得的光声信号执行模拟:其中测量目标的光吸收系数比构成对象的周围组织的光吸收系数高的情况;以及其中测量目标的光吸收系数比构成对象的周围组织的光吸收系数低的情况。
在图7A和图7B中示出了模拟模型。具体来说,设定二维空间,并且对于在执行用来自图的底部的光630的照射时的点640处获得的光声信号执行模拟。作为具有1cm的直径的圆形的作为测量目标的光吸收体610被放置在距离光照射点1cm处。音速是1500m/s。
背景区域620的光吸收系数是0.1cm-1,并且其等效的散射系数是10cm-1。光吸收体610的光吸收系数是1.0cm-1,并且其等效的散射系数是10cm-1。光吸收体710的光吸收系数是0cm-1,并且其等效的散射系数是10cm-1。
在图8A和图8B中示出了关于在点640处的声压的时移(time-lapse)变化的模拟结果。在具有比区域620的光吸收系数高的光吸收系数的光吸收体610的情况下,获得其中来自测量目标的信号在开始点处上升的正的光声信号(图8A)。另一方面,在具有比区域620的光吸收系数低的光吸收系数的光吸收体710的情况下,获得其中来自测量目标的信号在开始点处下降的负的光声信号(图8B)。
因此,应理解,在光吸收系数比周围区域的光吸收系数高时产生在开始点处上升的正的光声信号,而在光吸收系数比周围区域的光吸收系数低时产生在开始点处下降的负的光声信号。
虽然已经参考示例性实施例描述了本发明,但是应当理解,本发明不限于所公开的示例性实施例。以下权利要求的范围将被给予最宽的解释从而包括所有这样的修改、等同的结构与功能。
本申请要求2008年12月11日提交的日本专利申请No.2008-316042的权益,该日本专利申请其整体通过参考被并入于此。
Claims (10)
1.一种光声成像设备,包括:
光源;
检测器,所述检测器被配置为检测从已经吸收从所述光源发射的光的能量的光吸收体产生的声波;以及
信号处理单元,所述信号处理单元被配置为通过在对由检测器检测的声波执行波形处理之前存储指示所述由检测器检测的声波的压力的变化率是正还是负的信息来形成所述光吸收体的图像。
2.根据权利要求1所述的光声成像设备,其中所述波形处理是包络检波。
3.根据权利要求1或2所述的光声成像设备,其中,在声波的压力的变化率在开始时为负的情况下,所述信号处理单元输出用于显示指示所述光吸收体具有比所述光吸收体的周围区域的光吸收系数低的光吸收系数的信息的信号。
4.根据权利要求1到3中的任何一个所述的光声成像设备,其中,在声波的压力的变化率在开始时为负的情况下,所述信号处理单元将通过所述波形处理获得的声压信息从正值转换为负值或者从负值转换为正值,并且基于转换后的声压信息来形成光吸收体的图像。
5.根据权利要求1到4中的任何一个所述的光声成像设备,其中所述检测器是聚焦换能器。
6.根据权利要求3到5中的任何一个所述的光声成像设备,还包括图像显示单元,
其中所述图像显示单元通过使用与周围区域的显示颜色不同的显示颜色或者通过使用与周围区域的色调不同的色调来显示具有比周围区域的光吸收系数低的光吸收系数的光吸收体的图像。
7.根据权利要求3到5中的任何一个所述的光声成像设备,还包括图像显示单元,
其中所述图像显示单元显示具有比周围区域的光吸收系数低的光吸收系数的组织的图像,其中识别符号被添加到所述图像。
8.根据权利要求1到7中的任何一个所述的光声成像设备,其中在存储指示声波的压力的变化率是正还是负的信息之前,所述信号处理单元根据所述检测器的频率响应来对由所述检测器检测的声波执行校正。
9.根据权利要求1到8中的任何一个所述的光声成像设备,
其中所述光源能够发射具有第一波长的光和具有与第一波长不同的第二波长的光,并且
其中所述信号处理单元通过将通过用具有第一波长的光的照射获得的第一声波的形状与通过用具有第二波长的光的照射获得的第二声波的形状进行比较来识别所述光吸收体。
10.一种光声成像方法,包括:
从光源发射光的步骤;
检测从已经吸收从所述光源发射的光的能量的光吸收体产生的声波的步骤;以及
在对检测的声波执行波形处理之前存储指示所述检测的声波的压力的变化率是正还是负的信息的步骤。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2008-316042 | 2008-12-11 | ||
JP2008316042A JP5241465B2 (ja) | 2008-12-11 | 2008-12-11 | 光音響イメージング装置および光音響イメージング方法 |
PCT/JP2009/006758 WO2010067608A1 (en) | 2008-12-11 | 2009-12-10 | Photoacoustic imaging apparatus and photoacoustic imaging method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102238903A true CN102238903A (zh) | 2011-11-09 |
CN102238903B CN102238903B (zh) | 2014-09-24 |
Family
ID=41716386
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200980148696.XA Expired - Fee Related CN102238903B (zh) | 2008-12-11 | 2009-12-10 | 光声成像设备和光声成像方法 |
Country Status (5)
Country | Link |
---|---|
US (2) | US20110239766A1 (zh) |
EP (1) | EP2375966B1 (zh) |
JP (1) | JP5241465B2 (zh) |
CN (1) | CN102238903B (zh) |
WO (1) | WO2010067608A1 (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104780832B (zh) * | 2012-11-12 | 2017-06-13 | 佳能株式会社 | 被检体信息获取装置及被检体信息获取装置的控制方法 |
CN111316087A (zh) * | 2017-10-27 | 2020-06-19 | 柏林洪堡大学 | 借助具有预定的波长范围的测量光的、用于确定不均匀样品的特征的光声方法 |
CN111699376A (zh) * | 2017-12-28 | 2020-09-22 | Avl李斯特有限公司 | 用于求得测量气体的测量参数的测量装置 |
CN111743511A (zh) * | 2019-03-26 | 2020-10-09 | 株式会社爱德万测试 | 光声波测定装置 |
CN111973161A (zh) * | 2020-09-09 | 2020-11-24 | 南京诺源医疗器械有限公司 | 一种肿瘤检测和成像系统 |
CN113812921A (zh) * | 2020-06-18 | 2021-12-21 | 株式会社爱德万测试 | 光超声波测量装置、方法、程序、记录介质 |
Families Citing this family (33)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7750536B2 (en) | 2006-03-02 | 2010-07-06 | Visualsonics Inc. | High frequency ultrasonic transducer and matching layer comprising cyanoacrylate |
US9380967B2 (en) | 2007-04-11 | 2016-07-05 | The Board Of Regents Of The University Of Texas System | Systems and methods for measuring fetal cerebral oxygenation |
US10226206B2 (en) | 2007-04-11 | 2019-03-12 | The Board Of Regents Of The University Of Texas System | Systems and methods for measuring neonatal cerebral oxygenation |
JP5818444B2 (ja) * | 2010-02-04 | 2015-11-18 | キヤノン株式会社 | 機能情報取得装置、機能情報取得方法、及びプログラム |
JP5751769B2 (ja) * | 2010-07-27 | 2015-07-22 | キヤノン株式会社 | 画像情報取得装置及びその制御方法 |
US9289191B2 (en) | 2011-10-12 | 2016-03-22 | Seno Medical Instruments, Inc. | System and method for acquiring optoacoustic data and producing parametric maps thereof |
JP5653882B2 (ja) * | 2010-10-27 | 2015-01-14 | 富士フイルム株式会社 | 光音響撮像装置およびその作動方法 |
JP5677044B2 (ja) * | 2010-11-19 | 2015-02-25 | キヤノン株式会社 | 光音響測定装置及び方法 |
JP5777394B2 (ja) * | 2011-04-28 | 2015-09-09 | 富士フイルム株式会社 | 光音響画像化方法および装置 |
KR101287099B1 (ko) | 2011-10-11 | 2013-07-16 | 서강대학교산학협력단 | 석회화 조직을 검출하기 위한 광음향 영상 생성 방법 및 장치 |
US9733119B2 (en) | 2011-11-02 | 2017-08-15 | Seno Medical Instruments, Inc. | Optoacoustic component utilization tracking |
US20130289381A1 (en) | 2011-11-02 | 2013-10-31 | Seno Medical Instruments, Inc. | Dual modality imaging system for coregistered functional and anatomical mapping |
WO2013067419A1 (en) * | 2011-11-02 | 2013-05-10 | Seno Medical Instruments, Inc. | Dual modality imaging system for coregistered functional and anatomical mapping |
US11287309B2 (en) | 2011-11-02 | 2022-03-29 | Seno Medical Instruments, Inc. | Optoacoustic component utilization tracking |
US11191435B2 (en) | 2013-01-22 | 2021-12-07 | Seno Medical Instruments, Inc. | Probe with optoacoustic isolator |
US10433732B2 (en) | 2011-11-02 | 2019-10-08 | Seno Medical Instruments, Inc. | Optoacoustic imaging system having handheld probe utilizing optically reflective material |
JP2013128759A (ja) * | 2011-11-22 | 2013-07-04 | Fujifilm Corp | 光音響信号処理装置及び方法 |
JP2013128761A (ja) * | 2011-11-22 | 2013-07-04 | Fujifilm Corp | 光音響信号処理装置及び方法 |
WO2013134772A2 (en) | 2012-03-09 | 2013-09-12 | Seno Medical Instruments, Inc. | Statistical mapping in an optoacoustic imaging system |
JP6053339B2 (ja) * | 2012-06-06 | 2016-12-27 | キヤノン株式会社 | 被検体情報取得装置および被検体情報取得方法 |
US9610043B2 (en) | 2012-06-13 | 2017-04-04 | Seno Medical Instruments, Inc. | System and method for producing parametric maps of optoacoustic data |
MX2014015274A (es) * | 2012-06-13 | 2015-02-20 | Seno Medical Instr Inc | Sistema y metodo para producir mapas parametricos de datos optoacusticos. |
US20140182383A1 (en) * | 2012-12-28 | 2014-07-03 | Canon Kabushiki Kaisha | Object information obtaining device, display method, and non-transitory computer-readable storage medium |
JP6112861B2 (ja) * | 2012-12-28 | 2017-04-12 | キヤノン株式会社 | 被検体情報取得装置、信号処理装置および表示装置 |
JP6177530B2 (ja) * | 2013-01-18 | 2017-08-09 | 富士フイルム株式会社 | ドプラ計測装置およびドプラ計測方法 |
JP5936559B2 (ja) * | 2013-01-18 | 2016-06-22 | 富士フイルム株式会社 | 光音響画像生成装置および光音響画像生成方法 |
TWI662540B (zh) | 2014-03-12 | 2019-06-11 | 美商富士膠片索諾聲公司 | 具有超音波透鏡與整合中央匹配層的高頻超音波轉換器 |
CN106999113A (zh) | 2014-07-08 | 2017-08-01 | 德克萨斯大学系统董事会 | 用于测量胎儿脑氧合的系统和方法 |
AU2016233575A1 (en) | 2015-03-14 | 2017-07-13 | The Board Of Regents Of The University Of Texas System | Systems and methods for measuring neonatal cerebral oxygenation |
KR101709272B1 (ko) | 2015-07-08 | 2017-02-23 | 연세대학교 산학협력단 | 유도 광음향 분자진동 이미징 시스템 |
KR102450680B1 (ko) * | 2016-12-06 | 2022-10-06 | 한국전자통신연구원 | 입자 빔의 깊이 프로파일 측정방법 |
TWI617294B (zh) * | 2016-12-15 | 2018-03-11 | 財團法人金屬工業研究發展中心 | 用於光聲顯像裝置測試的仿生模擬系統 |
JP7277212B2 (ja) * | 2019-03-28 | 2023-05-18 | キヤノン株式会社 | 画像処理装置、画像処理方法及びプログラム |
Family Cites Families (36)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4385634A (en) * | 1981-04-24 | 1983-05-31 | University Of Arizona Foundation | Radiation-induced thermoacoustic imaging |
US4607341A (en) * | 1984-03-05 | 1986-08-19 | Canadian Patents And Development Limited | Device for determining properties of materials from a measurement of ultrasonic absorption |
NL8401041A (nl) * | 1984-04-03 | 1985-11-01 | Philips Nv | Elektro-akoestische inrichting. |
JPH0827264B2 (ja) * | 1988-09-21 | 1996-03-21 | 工業技術院長 | マルチ変調周波数による光音響撮像方法 |
JP3076945B2 (ja) * | 1993-06-15 | 2000-08-14 | 松下電器産業株式会社 | 吸音装置 |
US6216540B1 (en) * | 1995-06-06 | 2001-04-17 | Robert S. Nelson | High resolution device and method for imaging concealed objects within an obscuring medium |
US5657754A (en) * | 1995-07-10 | 1997-08-19 | Rosencwaig; Allan | Apparatus for non-invasive analyses of biological compounds |
US5840023A (en) * | 1996-01-31 | 1998-11-24 | Oraevsky; Alexander A. | Optoacoustic imaging for medical diagnosis |
US6309352B1 (en) * | 1996-01-31 | 2001-10-30 | Board Of Regents, The University Of Texas System | Real time optoacoustic monitoring of changes in tissue properties |
US5944687A (en) * | 1996-04-24 | 1999-08-31 | The Regents Of The University Of California | Opto-acoustic transducer for medical applications |
JPH1038856A (ja) * | 1996-07-22 | 1998-02-13 | Nikon Corp | 光吸収率測定装置及び測定方法 |
US5713356A (en) * | 1996-10-04 | 1998-02-03 | Optosonics, Inc. | Photoacoustic breast scanner |
GB9704737D0 (en) * | 1997-03-07 | 1997-04-23 | Optel Instr Limited | Biological measurement system |
GB9712524D0 (en) * | 1997-06-16 | 1997-08-20 | Nycomed Imaging As | Method |
US6139543A (en) * | 1998-07-22 | 2000-10-31 | Endovasix, Inc. | Flow apparatus for the disruption of occlusions |
IT1306874B1 (it) * | 1999-06-25 | 2001-10-11 | Consiglio Nazionale Ricerche | Metodo acustico non invasivo per la rivelazione e la mappatura deidistacchi nei dipinti ad affresco e nei dipinti murali in genere, e |
US6466806B1 (en) * | 2000-05-17 | 2002-10-15 | Card Guard Scientific Survival Ltd. | Photoacoustic material analysis |
IL138073A0 (en) * | 2000-08-24 | 2001-10-31 | Glucon Inc | Photoacoustic assay and imaging system |
AU2002341372A1 (en) * | 2001-10-09 | 2003-04-22 | Glucon Inc. | Method and apparatus for determining absorption of electromagnetic radiation by a material |
EP1593091B1 (en) * | 2003-01-25 | 2013-07-17 | Seno Medical Instruments, Inc. | High contrast optoacoustical imaging using non-spherical nanoparticles |
AU2004264839B2 (en) * | 2003-06-13 | 2009-06-04 | Becton, Dickinson And Company | Improved intra-dermal delivery of biologically active agents |
US20050107694A1 (en) * | 2003-11-17 | 2005-05-19 | Jansen Floribertus H. | Method and system for ultrasonic tagging of fluorescence |
JP4643153B2 (ja) * | 2004-02-06 | 2011-03-02 | 株式会社東芝 | 非侵襲生体情報映像装置 |
CN101305903B (zh) * | 2004-05-06 | 2011-09-28 | 日本电信电话株式会社 | 成分浓度测定装置及成分浓度测定装置控制方法 |
US20090054763A1 (en) * | 2006-01-19 | 2009-02-26 | The Regents Of The University Of Michigan | System and method for spectroscopic photoacoustic tomography |
US9439571B2 (en) * | 2006-01-20 | 2016-09-13 | Washington University | Photoacoustic and thermoacoustic tomography for breast imaging |
US8239006B2 (en) * | 2006-07-06 | 2012-08-07 | The University Of Connecticut | Method and apparatus for medical imaging using near-infrared optical tomography and fluorescence tomography combined with ultrasound |
WO2008011112A2 (en) * | 2006-07-19 | 2008-01-24 | University Of Connecticut | Method and apparatus for medical imaging using combined near-infrared optical tomography, fluorescent tomography and ultrasound |
US20080123083A1 (en) * | 2006-11-29 | 2008-05-29 | The Regents Of The University Of Michigan | System and Method for Photoacoustic Guided Diffuse Optical Imaging |
KR20090088909A (ko) * | 2006-12-19 | 2009-08-20 | 코닌클리케 필립스 일렉트로닉스 엔.브이. | 결합된 광음향 및 초음파 이미징 시스템 |
US20080173093A1 (en) * | 2007-01-18 | 2008-07-24 | The Regents Of The University Of Michigan | System and method for photoacoustic tomography of joints |
WO2008100386A2 (en) * | 2007-02-09 | 2008-08-21 | Board Of Regents, The University Of Texas System | Intravascular photoacoustic and ultrasound echo imaging |
WO2008103982A2 (en) * | 2007-02-23 | 2008-08-28 | The Regents Of The University Of Michigan | System and method for monitoring photodynamic therapy |
JP4739363B2 (ja) * | 2007-05-15 | 2011-08-03 | キヤノン株式会社 | 生体情報イメージング装置、生体情報の解析方法、及び生体情報のイメージング方法 |
EP2002784B1 (en) * | 2007-06-11 | 2018-07-11 | Canon Kabushiki Kaisha | Intravital-information imaging apparatus |
CN102131463B (zh) * | 2008-08-27 | 2013-01-16 | 佳能株式会社 | 生物信息处理设备和生物信息处理方法 |
-
2008
- 2008-12-11 JP JP2008316042A patent/JP5241465B2/ja not_active Expired - Fee Related
-
2009
- 2009-12-10 CN CN200980148696.XA patent/CN102238903B/zh not_active Expired - Fee Related
- 2009-12-10 WO PCT/JP2009/006758 patent/WO2010067608A1/en active Application Filing
- 2009-12-10 US US13/132,998 patent/US20110239766A1/en not_active Abandoned
- 2009-12-10 EP EP09798972.7A patent/EP2375966B1/en active Active
-
2014
- 2014-01-14 US US14/154,504 patent/US9032800B2/en not_active Expired - Fee Related
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104780832B (zh) * | 2012-11-12 | 2017-06-13 | 佳能株式会社 | 被检体信息获取装置及被检体信息获取装置的控制方法 |
US9782080B2 (en) | 2012-11-12 | 2017-10-10 | Canon Kabushiki Kaisha | Object information acquiring apparatus and control method for the object information acquiring apparatus |
CN111316087A (zh) * | 2017-10-27 | 2020-06-19 | 柏林洪堡大学 | 借助具有预定的波长范围的测量光的、用于确定不均匀样品的特征的光声方法 |
CN111316087B (zh) * | 2017-10-27 | 2021-09-28 | 柏林洪堡大学 | 借助具有预定的波长范围的测量光的、用于确定不均匀样品的特征的光声方法 |
CN111699376A (zh) * | 2017-12-28 | 2020-09-22 | Avl李斯特有限公司 | 用于求得测量气体的测量参数的测量装置 |
CN111699376B (zh) * | 2017-12-28 | 2023-09-05 | Avl李斯特有限公司 | 用于求得测量气体的测量参数的测量装置 |
CN111743511A (zh) * | 2019-03-26 | 2020-10-09 | 株式会社爱德万测试 | 光声波测定装置 |
CN111743511B (zh) * | 2019-03-26 | 2023-07-21 | 株式会社爱德万测试 | 光声波测定装置 |
CN113812921A (zh) * | 2020-06-18 | 2021-12-21 | 株式会社爱德万测试 | 光超声波测量装置、方法、程序、记录介质 |
CN111973161A (zh) * | 2020-09-09 | 2020-11-24 | 南京诺源医疗器械有限公司 | 一种肿瘤检测和成像系统 |
Also Published As
Publication number | Publication date |
---|---|
US20140128718A1 (en) | 2014-05-08 |
EP2375966A1 (en) | 2011-10-19 |
US9032800B2 (en) | 2015-05-19 |
JP2010136887A (ja) | 2010-06-24 |
JP5241465B2 (ja) | 2013-07-17 |
WO2010067608A1 (en) | 2010-06-17 |
CN102238903B (zh) | 2014-09-24 |
EP2375966B1 (en) | 2022-02-23 |
US20110239766A1 (en) | 2011-10-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102238903B (zh) | 光声成像设备和光声成像方法 | |
CN103300880B (zh) | 被检体信息获得装置和被检体信息获得方法 | |
EP1313396B1 (en) | Photoacoustic assay and imaging system | |
JP5235586B2 (ja) | 生体情報処理装置及び生体情報処理方法 | |
US10653322B2 (en) | Photoacoustic apparatus, method of acquiring subject information, and non-transitory computer readable medium | |
JP6501474B2 (ja) | 被検体情報取得装置 | |
EP2946725B1 (en) | Photoacoustic apparatus | |
CN106560160A (zh) | 被检体信息获取装置及其控制方法 | |
JP6226523B2 (ja) | 被検体情報取得装置、表示方法、およびデータ処理装置 | |
CN103126727A (zh) | 超声成像设备和显示超声图像的方法 | |
US10582910B2 (en) | Information acquisition apparatus and information acquisition method | |
KR20170077037A (ko) | 정보 취득장치, 신호 처리방법, 및 프로그램 | |
CN105395165A (zh) | 被检体信息获取装置 | |
JP2011092631A (ja) | 生体情報処理装置及び生体情報処理方法 | |
US10436706B2 (en) | Information processing apparatus, information processing method, and storage medium | |
JP2017529913A (ja) | 光音響装置 | |
US20150327769A1 (en) | Photoacoustic apparatus | |
JP2017164198A (ja) | 情報処理システムおよび表示制御方法 | |
CN105640496A (zh) | 光声装置和被检体信息取得方法 | |
CN108403083A (zh) | 信息处理装置、信息处理方法和非暂时性存储介质 | |
JP5575293B2 (ja) | 被検体情報取得装置及び被検体情報取得方法 | |
JP6490171B2 (ja) | データ処理装置およびデータ処理方法 | |
CN106455994A (zh) | 光声装置 | |
JP2015216983A (ja) | 光音響装置 | |
US10172524B2 (en) | Photoacoustic apparatus |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140924 |
|
CF01 | Termination of patent right due to non-payment of annual fee |