JP5713968B2 - 光音響画像生成装置及び音響波ユニット - Google Patents
光音響画像生成装置及び音響波ユニット Download PDFInfo
- Publication number
- JP5713968B2 JP5713968B2 JP2012160427A JP2012160427A JP5713968B2 JP 5713968 B2 JP5713968 B2 JP 5713968B2 JP 2012160427 A JP2012160427 A JP 2012160427A JP 2012160427 A JP2012160427 A JP 2012160427A JP 5713968 B2 JP5713968 B2 JP 5713968B2
- Authority
- JP
- Japan
- Prior art keywords
- photoacoustic
- wavelength
- unit
- laser light
- photoacoustic image
- 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
- 230000003287 optical effect Effects 0.000 claims description 100
- 238000006243 chemical reaction Methods 0.000 claims description 46
- 238000001514 detection method Methods 0.000 claims description 32
- 230000005284 excitation Effects 0.000 claims description 29
- 230000005540 biological transmission Effects 0.000 claims description 28
- 238000000605 extraction Methods 0.000 claims description 24
- 238000006073 displacement reaction Methods 0.000 claims description 22
- 238000000034 method Methods 0.000 claims description 21
- 238000010276 construction Methods 0.000 claims description 20
- 239000000284 extract Substances 0.000 claims description 12
- 230000008859 change Effects 0.000 claims description 8
- 230000001678 irradiating effect Effects 0.000 claims description 6
- 229960001716 benzalkonium Drugs 0.000 claims 1
- CYDRXTMLKJDRQH-UHFFFAOYSA-N benzododecinium Chemical compound CCCCCCCCCCCC[N+](C)(C)CC1=CC=CC=C1 CYDRXTMLKJDRQH-UHFFFAOYSA-N 0.000 claims 1
- 239000000523 sample Substances 0.000 description 34
- 238000005070 sampling Methods 0.000 description 19
- 238000002604 ultrasonography Methods 0.000 description 15
- 108010054147 Hemoglobins Proteins 0.000 description 13
- 102000001554 Hemoglobins Human genes 0.000 description 13
- 239000006096 absorbing agent Substances 0.000 description 11
- 238000010521 absorption reaction Methods 0.000 description 8
- 210000001367 artery Anatomy 0.000 description 8
- 210000003462 vein Anatomy 0.000 description 8
- 238000003384 imaging method Methods 0.000 description 7
- 238000012545 processing Methods 0.000 description 7
- 230000031700 light absorption Effects 0.000 description 6
- 238000000926 separation method Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 230000004044 response Effects 0.000 description 5
- 238000013459 approach Methods 0.000 description 4
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 4
- 238000009826 distribution Methods 0.000 description 4
- 239000000203 mixture Substances 0.000 description 4
- 229910052760 oxygen Inorganic materials 0.000 description 4
- 239000001301 oxygen Substances 0.000 description 4
- 239000013078 crystal Substances 0.000 description 3
- 238000003780 insertion Methods 0.000 description 3
- 230000037431 insertion Effects 0.000 description 3
- 230000002194 synthesizing effect Effects 0.000 description 3
- 239000008280 blood Substances 0.000 description 2
- 210000004369 blood Anatomy 0.000 description 2
- 210000004204 blood vessel Anatomy 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000010355 oscillation Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 108010064719 Oxyhemoglobins Proteins 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 230000001151 other effect Effects 0.000 description 1
- 238000010895 photoacoustic effect Methods 0.000 description 1
- 239000002243 precursor Substances 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 229910052594 sapphire Inorganic materials 0.000 description 1
- 239000010980 sapphire Substances 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 230000005236 sound signal Effects 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000003325 tomography 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
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/13—Tomography
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01S—DEVICES USING THE PROCESS OF LIGHT AMPLIFICATION BY STIMULATED EMISSION OF RADIATION [LASER] TO AMPLIFY OR GENERATE LIGHT; DEVICES USING STIMULATED EMISSION OF ELECTROMAGNETIC RADIATION IN WAVE RANGES OTHER THAN OPTICAL
- H01S3/00—Lasers, i.e. devices using stimulated emission of electromagnetic radiation in the infrared, visible or ultraviolet wave range
- H01S3/10—Controlling the intensity, frequency, phase, polarisation or direction of the emitted radiation, e.g. switching, gating, modulating or demodulating
- H01S3/106—Controlling the intensity, frequency, phase, polarisation or direction of the emitted radiation, e.g. switching, gating, modulating or demodulating by controlling devices placed within the cavity
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/145—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
- A61B5/14542—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue for measuring blood gases
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0891—Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of blood vessels
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/46—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
- A61B8/461—Displaying means of special interest
- A61B8/463—Displaying means of special interest characterised by displaying multiple images or images and diagnostic data on one display
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5207—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
- A61B8/5238—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image
- A61B8/5261—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image combining images from different diagnostic modalities, e.g. ultrasound and X-ray
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/54—Control of the diagnostic device
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01S—DEVICES USING THE PROCESS OF LIGHT AMPLIFICATION BY STIMULATED EMISSION OF RADIATION [LASER] TO AMPLIFY OR GENERATE LIGHT; DEVICES USING STIMULATED EMISSION OF ELECTROMAGNETIC RADIATION IN WAVE RANGES OTHER THAN OPTICAL
- H01S3/00—Lasers, i.e. devices using stimulated emission of electromagnetic radiation in the infrared, visible or ultraviolet wave range
- H01S3/05—Construction or shape of optical resonators; Accommodation of active medium therein; Shape of active medium
- H01S3/08—Construction or shape of optical resonators or components thereof
- H01S3/08018—Mode suppression
- H01S3/08022—Longitudinal modes
- H01S3/08027—Longitudinal modes by a filter, e.g. a Fabry-Perot filter is used for wavelength setting
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01S—DEVICES USING THE PROCESS OF LIGHT AMPLIFICATION BY STIMULATED EMISSION OF RADIATION [LASER] TO AMPLIFY OR GENERATE LIGHT; DEVICES USING STIMULATED EMISSION OF ELECTROMAGNETIC RADIATION IN WAVE RANGES OTHER THAN OPTICAL
- H01S3/00—Lasers, i.e. devices using stimulated emission of electromagnetic radiation in the infrared, visible or ultraviolet wave range
- H01S3/09—Processes or apparatus for excitation, e.g. pumping
- H01S3/091—Processes or apparatus for excitation, e.g. pumping using optical pumping
- H01S3/0915—Processes or apparatus for excitation, e.g. pumping using optical pumping by incoherent light
- H01S3/092—Processes or apparatus for excitation, e.g. pumping using optical pumping by incoherent light of flash lamp
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01S—DEVICES USING THE PROCESS OF LIGHT AMPLIFICATION BY STIMULATED EMISSION OF RADIATION [LASER] TO AMPLIFY OR GENERATE LIGHT; DEVICES USING STIMULATED EMISSION OF ELECTROMAGNETIC RADIATION IN WAVE RANGES OTHER THAN OPTICAL
- H01S3/00—Lasers, i.e. devices using stimulated emission of electromagnetic radiation in the infrared, visible or ultraviolet wave range
- H01S3/10—Controlling the intensity, frequency, phase, polarisation or direction of the emitted radiation, e.g. switching, gating, modulating or demodulating
- H01S3/11—Mode locking; Q-switching; Other giant-pulse techniques, e.g. cavity dumping
- H01S3/1123—Q-switching
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- General Health & Medical Sciences (AREA)
- Animal Behavior & Ethology (AREA)
- Surgery (AREA)
- Electromagnetism (AREA)
- Acoustics & Sound (AREA)
- Plasma & Fusion (AREA)
- Optics & Photonics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Lasers (AREA)
Description
本発明は、光音響画像生成装置及び音響波ユニットに関し、更に詳しくは、複数波長のレーザ光を被検体に照射したときに各波長に対して検出された光音響信号に基づいて光音響画像を生成する光音響画像生成装置及び音響波ユニットに関する。
従来、例えば特許文献１や非特許文献１に示されているように、光音響効果を利用して生体の内部を画像化する光音響画像化装置が知られている。この光音響画像化装置においては、例えばパルスレーザ光等のパルス光が生体に照射される。このパルス光の照射を受けた生体内部では、パルス光のエネルギーを吸収した生体組織が熱によって体積膨張し、音響波が発生する。この音響波を超音波プローブなどで検出し、検出された信号（光音響信号）に基づいて生体内部を可視像化することが可能となっている。光音響画像化方法では、特定の光吸収体において音響波が発生するため、生体における特定の組織、例えば血管等を画像化することができる。
ところで、生体組織の多くは光吸収特性が光の波長に応じて変わり、また一般に、その光吸収特性も組織ごとに特有のものとなっている。例えば図８に、ヒトの動脈に多く含まれる酸素化ヘモグロビン(酸素と結合したヘモグロビン:oxy-Hb)と、静脈に多く含まれる
脱酸素化ヘモグロビン（酸素と結合していないヘモグロビンdeoxy-Hb)の光波長ごとの分
子吸収係数を示す。動脈の光吸収特性は、酸素化ヘモグロビンのそれに対応し、静脈の光吸収特性は、脱酸素化ヘモグロビンのそれに対応する。この波長に応じた光吸収率の違いを利用して、互いに異なる２種の波長の光を血管部分に照射し、動脈と静脈とを区別して画像化する光音響画像化方法が知られている（例えば特許文献２参照）。
脱酸素化ヘモグロビン（酸素と結合していないヘモグロビンdeoxy-Hb)の光波長ごとの分
子吸収係数を示す。動脈の光吸収特性は、酸素化ヘモグロビンのそれに対応し、静脈の光吸収特性は、脱酸素化ヘモグロビンのそれに対応する。この波長に応じた光吸収率の違いを利用して、互いに異なる２種の波長の光を血管部分に照射し、動脈と静脈とを区別して画像化する光音響画像化方法が知られている（例えば特許文献２参照）。
ここで、可変波長レーザに関して、特許文献３には、複数種の波長のレーザ光を容易に切り換えて出力することが可能な多色の固体レーザ装置が記載されている。特許文献３では、特定のピーク波長光だけを選択的に透過させるフィルタを、レーザ活性媒体と光共振器ミラーのうちの一方との間の光路上に配置する。フィルタを、選択すべきピーク波長の分だけ用意し、用意されたフィルタのうちの何れかを光路上に配置することで、複数の波長のレーザ光を切り替えて出射することができる。
A High-Speed Photoacoustic Tomography System based on a Commercial Ultrasound and a Custom Transducer Array, Xueding Wang, Jonathan Cannata, Derek DeBusschere, Changhong Hu, J. Brian Fowlkes, and Paul Carson, Proc. SPIE Vol. 7564, 756424 (Feb.23, 2010)
光音響イメージングにおいては、複数フレームにわたって光音響画像を生成する場合などに、例えば１セットあたり２つの波長で構成されるパルスレーザ光を複数セット続けて被検体に照射したい場合がある。特許文献３に記載されているのは、単に、光路中に挿入されるフィルタを切り替えることでレーザ光の波長を切り替えることでしかなく、特許文献３は、複数波長のレーザ光を連続的に切り換えて出射するための手段を提供していない。また、引用文献３は、レーザ光源の外部から、レーザ光源が出射するレーザ光の波長を制御するための手段を提供していない。
本発明は、上記に鑑み、光音響イメージングにおいて、レーザ光源の外部から複数の波長のレーザ光を切り替えて出射させることが可能な音響波ユニット、及びそのような音響波ユニットを含む光音響画像生成装置を提供することを目的とする。
上記目的を達成するために、本発明は、相互に異なる複数の波長のパルスレーザ光を出射するレーザ光源ユニットであって、レーザロッドと、該レーザロッドに励起光を照射する励起光源と、レーザロッドを挟んで対向する一対のミラーを含む光共振器と、光共振器内に挿入されたＱスイッチと、透過波長が相互に異なる複数のバンドパスフィルタを含み、該複数のバンドパスフィルタを光共振器の光路上に選択的に挿入する波長選択手段と、光共振器の光路上に挿入されるバンドパスフィルタが所定の順序で順次に切り替わるように波長選択手段を駆動する駆動手段とを有するレーザ光源ユニットと、光音響画像の生成を行う音響波ユニットであって、複数の波長のパルスレーザ光が被検体に照射されたときに被検体内で生じる光音響信号を検出し、各波長に対応した光音響データを生成する検出手段と、各波長に対応した光音響データ間の相対的な信号強度の大小関係を抽出する強度比抽出手段と、抽出された大小関係に基づいて光音響画像を生成する光音響画像構築手段と、光共振器の光路上に挿入されるバンドパスフィルタが所定の切替え速度で切り替わるように駆動手段を制御しつつ、励起光源からレーザロッドに励起光を照射させ、該励起光の照射後、波長選択手段が、出射すべきパルスレーザ光の波長に対応した透過波長のバンドパスフィルタを光路上に挿入しているタイミングでＱスイッチをオンにしてパルスレーザ光を出射させるトリガ制御回路とを有する音響波ユニットとを備えたことを特徴とする光音響画像生成装置を提供する。
本発明では、波長選択手段が、回転変位に伴って光共振器の光路上に選択的に挿入するバンドパスフィルタを切り替えるフィルタ回転体で構成され、駆動手段がフィルタ回転体を回転駆動するものとすることができる。
トリガ制御回路は、フィルタ回転体が、所定の方向に所定の回転速度で連続的に回転するように駆動手段を制御することが好ましい。
所定の回転速度が、出射すべきパルスレーザ光の波長の数と単位時間当たりのパルスレーザ光の発光回数とに基づいて決定されるものとすることができる。
トリガ制御回路は、フィルタ回転体の回転変位位置を示すフィルタ状態情報に基づいて、励起光を照射させるタイミング及びＱスイッチをオンにするタイミングを決定してもよい。
トリガ制御回路が、フィルタ状態情報が、出射すべきパルスレーザ光の波長に対応したバンドパスフィルタが光路上に挿入されるフィルタ回転体の位置からレーザロッドの励起に要する時間の間にフィルタ回転体が回転変位する量を差し引いた位置を示す情報となると、レーザロッドに励起光を照射させるものとすることができる。
トリガ制御回路が、所定時間の間におけるフィルタ状態情報の変化の量が所定の回転速度に応じた変化量となるように、駆動手段を制御する構成を採用することができる。
レーザ光源ユニットが、所定時間の間にフィルタ回転体の回転変位の量が所定の回転速度に応じた量になるように駆動手段を制御する回転制御部を更に有する構成を採用でき、トリガ制御回路が、回転制御部を通じて駆動手段を制御することとしてもよい。
音響波ユニットが、各波長に対応した光音響データに基づいて信号強度を示す強度情報を生成する強度情報抽出手段を更に備える構成を採用でき、光音響画像構築手段が、光音響画像の各画素の階調値を強度情報に基づいて決定すると共に、各画素の表示色を抽出された大小関係に基づいて決定するものとすることができる。
レーザ光源ユニットが出射すべきパルスレーザ光の複数の波長が第１の波長と第２の波長を含み、音響波ユニットが、第１の波長のパルスレーザ光が照射されたときに検出された光音響信号に対応する第１の光音響データと、第２の波長のパルスレーザ光が照射されたときに検出された光音響信号に対応する第２の光音響データとのうちの何れか一方を実部、他方を虚部とした複素数データを生成する複素数化手段と、複素数データからフーリエ変換法により再構成画像を生成する光音響画像再構成手段とを更に備える構成を採用でき、強度比抽出手段が再構成画像から大小関係としての位相情報を抽出し、強度情報抽出手段が再構成画像から強度情報を抽出するものとすることができる。
検出手段が、更に、被検体に送信された音響波に対する反射音響波を検出して反射音響波データを生成するものであり、音響波ユニットが、反射音響波データに基づいて音響波画像を生成する音響波画像生成手段を更に備える構成としてもよい。
本発明は、また、相互に異なる複数の波長のパルスレーザ光が被検体に照射されたときに被検体内で生じた光音響信号を検出し、各波長に対応した光音響データを生成する検出手段と、各波長に対応した光音響データ間の相対的な信号強度の大小関係を抽出する強度比抽出手段と、抽出された大小関係に基づいて光音響画像を生成する光音響画像構築手段と、透過波長が相互に異なる複数のバンドパスフィルタを含む波長選択手段が、レーザロッドを挟んで対向する一対のミラーを含む光共振器の光路上に複数のバンドパスフィルタを所定の順序で選択的に挿入するように前駆波長選択手段を駆動する駆動手段を、光共振器の光路上に挿入されるバンドパスフィルタが所定の切替え速度で切り替わるように制御しつつ、励起光源からレーザロッドに励起光を照射させ、該励起光の照射後、波長選択手段が、出射すべきパルスレーザ光の波長に対応した透過波長のバンドパスフィルタを光路上に挿入しているタイミングで光共振器内に挿入されたＱスイッチをオンにしてパルスレーザ光を出射させるトリガ制御回路とを備えたことを特徴とする音響波ユニットを提供する。
本発明の光音響画像生成装置及び音響波ユニットは、透過波長が相互に異なる複数のバンドパスフィルタを含む波長選択手段を、光共振器の光路上に選択的に挿入される複数のバンドパスフィルタが所定の速度で切り替わるように駆動しつつ、ラッシュランプを点灯してレーザロッドを励起し、その励起後、光共振器の光路上に、出射すべきパルスレーザ光の波長の光を透過するバンドパスフィルタが挿入されているときにＱスイッチをオンにする。例えば、透過波長が異なる２つのバンドパスフィルタを含む波長選択手段を連続的に駆動することで、２つのバンドパスフィルタを連続的、かつ選択的に光共振器の光路上に挿入することができ、レーザ光源ユニットから複数波長のパルスレーザ光を連続的に切り換えて出射することができる。また、本発明では、Ｑスイッチをオンにするときに、何れのバンドパスフィルタが光共振器の光路上に挿入されているかに応じてパルスレーザ光の波長が決まり、音響波ユニット側から、レーザ光源ユニットが出射するパルスレーザ光の波長を制御することができる。更に、本発明では、音響波ユニットがパルスレーザ光の出射タイミングを決めており、光音響信号のサンプリング開始にあたり、レーザ光源ユニットからレーザ出射を示す同期信号のような信号を取得する必要がない。
以下、図面を参照し、本発明の実施の形態を詳細に説明する。なお、本発明の実施例では、音響波として超音波を用いるが、被検対象や測定条件等に応じて適切な周波数を選択することにより、可聴周波数の音響波であっても良い。
図１は、本発明の第１実施形態の光音響画像生成装置を示す。光音響画像生成装置１０は、超音波探触子（プローブ）１１と、超音波ユニット１２と、レーザ光源ユニット１３とを備える。レーザ光源ユニット１３は、被検体に照射すべきパルスレーザ光を出射する。レーザ光源ユニット１３は、相互に異なる複数の波長のパルスレーザ光を切り替えて出射する。以下の説明においては、主に、レーザ光源ユニット１３が、第１の波長のパルスレーザ光と第２の波長のパルスレーザ光とを順次に出射するものとして説明する。
例えば、第１の波長（中心波長）として約７５０ｎｍを考え、第２の波長として約８００ｎｍを考える。先に説明した図８を参照すると、ヒトの動脈に多く含まれる酸素化ヘモグロビン(酸素と結合したヘモグロビン:oxy-Hb)の波長７５０ｎｍにおける分子吸収係数は、波長８００ｎｍにおける分子吸収係数よりも低い。一方、静脈に多く含まれる脱酸素化ヘモグロビン（酸素と結合していないヘモグロビンdeoxy-Hb)の波長７５０ｎｍにおける分子吸収係数は、波長８００ｎｍにおける分子吸収係数よりも高い。この性質を利用し、波長８００ｎｍで得られた光音響信号に対して、波長７５０ｎｍで得られた光音響信号が相対的に大きいのか小さいのかを調べることで、動脈からの光音響信号と静脈からの光音響信号とを判別することができる。
レーザ光源ユニット１３から出射したパルスレーザ光は、例えば光ファイバなどの導光手段を用いてプローブ１１まで導光され、プローブ１１から被検体に向けて照射される。パルスレーザ光の照射位置は特に限定されず、プローブ１１以外の場所からパルスレーザ光の照射を行ってもよい。被検体内では、光吸収体が照射されたパルスレーザ光のエネルギーを吸収することで超音波（音響波）が生じる。プローブ１１は、超音波検出器を含む。プローブ１１は、例えば一次元的に配列された複数の超音波検出器素子（超音波振動子）を有し、その一次元配列された超音波振動子により、被検体内からの音響波（光音響信号）を検出する。
超音波ユニット１２は、受信回路２１、ＡＤ変換手段２２、受信メモリ２３、複素数化手段２４、光音響画像再構成手段２５、位相情報抽出手段２６、強度情報抽出手段２７、検波・対数変換手段２８、光音響画像構築手段２９、トリガ制御回路３０、及び、制御手段３１を有する。受信回路２１は、プローブ１１が検出した光音響信号を受信する。ＡＤ変換手段２２は検出手段であり、受信回路２１が受信した光音響信号をサンプリングし、デジタルデータである光音響データを生成する。ＡＤ変換手段２２は、ＡＤクロック信号に同期して、所定のサンプリング周期で光音響信号のサンプリングを行う。
ＡＤ変換手段２２は、光音響データを受信メモリ２３に格納する。ＡＤ変換手段２２は、レーザ光源ユニット１３から出射されるパルスレーザ光の各波長に対応した光音響データを受信メモリ２３に格納する。つまり、ＡＤ変換手段２２は、被検体に第１の波長のパルスレーザ光が照射されたときにプローブ１１で検出された光音響信号をサンプリングした第１の光音響データと、第２のパルスレーザ光が照射されたときにプローブ１１で検出された光音響信号をサンプリングした第２の光音響データとを、受信メモリ２３に格納する。
複素数化手段２４は、受信メモリ２３から第１の光音響データと第２の光音響データとを読み出し、何れか一方を実部、他方を虚部とした複素数データを生成する。以下では、複素数化手段２４が、第１の光音響データを実部とし、第２の光音響データを虚部とした複素数データを生成するものとして説明する。
光音響画像再構成手段２５は、複素数化手段２４から複素数データを入力する。光音響画像再構成手段２５は、入力された複素数データから、フーリエ変換法（ＦＴＡ法）により画像再構成を行う。フーリエ変換法による画像再構成には、例えば文献”Photoacoustic Image Reconstruction-A Quantitative Analysis”Jonathan I.Sperl et al. SPIE-OSA
Vol.6631 663103 等に記載されている従来公知の方法を適用することができる。光音響画像再構成手段２５は、再構成画像を示すフーリエ変換のデータを位相情報抽出手段２６と強度情報抽出手段２７とに入力する。
Vol.6631 663103 等に記載されている従来公知の方法を適用することができる。光音響画像再構成手段２５は、再構成画像を示すフーリエ変換のデータを位相情報抽出手段２６と強度情報抽出手段２７とに入力する。
位相情報抽出手段２６は、各波長に対応した光音響データ間の相対的な信号強度の大小関係を抽出する。本実施形態では、位相情報抽出手段２６は、光音響画像再構成手段２５で再構成された再構成画像を入力データとし、複素数データである入力データから、実部と虚部とを比較したときに、相対的に、どちらがどれくらい大きいかを示す位相情報を抽出する。位相情報抽出手段２６は、例えば複素数データがＸ＋ｉＹで表わされるとき、θ＝ｔａｎ−１（Ｙ／Ｘ）を位相情報として生成する。なお、Ｘ＝０の場合はθ＝９０°とする。実部を構成する第１の光音響データ（Ｘ）と虚部を構成する第２の光音響データ（Ｙ）とが等しいとき、位相情報はθ＝４５°となる。位相情報は、相対的に第１の光音響データが大きいほどθ＝０°に近づいていき、第２の光音響データが大きいほどθ＝９０°に近づいていく。
強度情報抽出手段２７は、各波長に対応した光音響データに基づいて信号強度を示す強度情報を生成する。本実施形態では、強度情報抽出手段２７は、光音響画像再構成手段２５で再構成された再構成画像を入力データとし、複素数データである入力データから、強度情報を生成する。強度情報抽出手段２７は、例えば複素数データがＸ＋ｉＹで表わされるとき、（Ｘ２＋Ｙ２）１／２を、強度情報として抽出する。検波・対数変換手段２８は、強度情報抽出手段２７で抽出された強度情報を示すデータの包絡線を生成し、次いでその包絡線を対数変換してダイナミックレンジを広げる。
光音響画像構築手段２９は、位相情報抽出手段２６から位相情報を入力し、検波・対数変換手段２８から検波・対数変換処理後の強度情報を入力する。光音響画像構築手段２９は、入力された位相情報と強度情報とに基づいて、光吸収体の分布画像である光音響画像を生成する。光音響画像構築手段２９は、例えば入力された強度情報に基づいて、光吸収体の分布画像における各画素の輝度（階調値）を決定する。また、光音響画像構築手段２９は、例えば位相情報に基づいて、光吸収体の分布画像における各画素の色（表示色）を決定する。光音響画像構築手段２９は、例えば例えば位相０°から９０°の範囲を所定の色に対応させたカラーマップに用いて、入力された位相情報に基づいて各画素の色を決定する。
ここで、位相０°から４５°の範囲は、第１の光音響データが第２の光音響データよりも大きい範囲であるため、光音響信号の発生源は、波長７９８ｎｍに対する吸収よりも波長７５６ｎｍに対する吸収の方が大きい脱酸素化ヘモグロビンを主に含む血液が流れている静脈であると考えられる。一方、位相４５°から９０°の範囲は、第１の光音響データが第２の光音響データよりも小さい範囲であるため、光音響信号の発生源は、波長７９８ｎｍに対する吸収よりも波長７５６ｎｍに対する吸収の方が小さい酸素化ヘモグロビンを主に含む血液が流れている動脈であると考えられる。
そこで、カラーマップとして、例えば位相が０°が青色で、位相が４５°に近づくに連れて無色（白色）になるように色が徐々に変化すると共に、位相９０°が赤色で、位相が４５°に近づくに連れて白色になるように色が徐々に変化するようなカラーマップを用いる。この場合、光音響画像上で、動脈に対応した部分を赤色で表わし、静脈に対応した部分を青色で表わすことができる。強度情報を用いずに、階調値は一定として、位相情報に従って動脈に対応した部分と静脈に対応した部分との色分けを行うだけでもよい。画像表示手段１４は、光音響画像構築手段２９が生成した光音響画像を表示画面上に表示する。
次いで、レーザ光源ユニット１３の構成を詳細に説明する。図２は、レーザ光源ユニット１３の構成を示す。レーザ光源ユニット１３は、レーザロッド５１、フラッシュランプ５２、ミラー５３、５４、Ｑスイッチ５５、波長選択手段５６、駆動手段５７、及び駆動状態検出手段５８を有する。レーザロッド５１は、レーザ媒質である。レーザロッド５１には、例えばアレキサンドライト結晶やＣｒ：ＬｉＳＡＦ（Ｃｒ：ＬｉＳｒＡｌＦ６），Ｃｒ：ＬｉＣＡＦ（Ｃｒ：ＬｉＣａＡｌＦ６）結晶，Ｔｉ：Ｓａｐｐｈｉｒｅ結晶を用いることができる。フラッシュランプ５２は、励起光源であり、レーザロッド５１に励起光を照射する。フラッシュランプ５２以外の光源、例えば半導体レーザや固体レーザを、励起光源として用いてもよい。
ミラー５３、５４は、レーザロッド５１を挟んで対向しており、ミラー５３、５４により光共振器が構成される。ミラー５４が出力側であるものとする。光共振器内には、Ｑスイッチ５５が挿入される。Ｑスイッチ５５により、光共振器内の挿入損失を損失大（低Ｑ）から損失小（高Ｑ）へと急速に変化させることで、パルスレーザ光を得ることができる。
波長選択手段５６は、透過波長が相互に異なる複数のバンドパスフィルタ（ＢＰＦ：Band Pass Filter）を含む。波長選択手段５６は、複数のバンドパスフィルタを光共振器の光路上に選択的に挿入する。波長選択手段５６は、例えば波長７５０ｎｍ（中心波長）の光を透過させる第１のバンドパスフィルタと、波長８００ｎｍ（中心波長）の光を透過させる第２のバンドパスフィルタとを含む。光共振器の光路上に第１のバンドパスフィルタを挿入することで、光発振器の発振波長を７５０ｎｍとすることができ、光共振器の光路上に第２のバンドパスフィルタを挿入することで、光発振器の発振波長を８００ｎｍとすることができる。
駆動手段５７は、光共振器の光路上に挿入されるバンドパスフィルタが所定の順序で順次に切り替わるように波長選択手段５６を駆動する。例えば波長選択手段５６が、回転変位に伴って光共振器の光路上に選択的に挿入するバンドパスフィルタを切り替えるフィルタ回転体で構成されているとき、駆動手段５７は、波長選択手段５６を構成するフィルタ回転体を回転駆動する。駆動状態検出手段５８は、波長選択手段５６の駆動状態を検出する。駆動状態検出手段５８は、例えばフィルタ回転体である波長選択手段５６の回転変位を検出する。駆動状態検出手段５８は、フィルタ回転体の回転変位位置を示すＢＰＦ状態情報を超音波ユニット１２に出力する。
図１に戻り、制御手段３１は、超音波ユニット１２内の各部の制御を行う。トリガ制御回路３０は、レーザ光源ユニット１３内の波長選択手段５６が光共振器の光路上に挿入するバンドパスフィルタが所定の切替え速度で切り替わるように駆動手段５７を制御する。トリガ制御回路３０は、例えば、波長選択手段５６を構成するフィルタ回転体が、所定の方向に所定の回転速度で連続的に回転するように駆動手段５７を制御する。フィルタ回転体の回転速度は、例えばレーザ光源ユニット１３から出射すべきパルスレーザ光の波長の数（バンドパスフィルタの透過波長の数）と、単位時間当たりのパルスレーザ光の発光回数とに基づいて決定できる。
トリガ制御回路３０は、波長選択手段５６の駆動を制御するためのＢＰＦ制御信号を出力する。レーザ光源ユニット１３の駆動手段５７は、ＢＰＦ制御信号に応じて波長選択手段５６を駆動する。トリガ制御回路３０は、ＢＰＦ制御信号を通じて、例えば所定時間の間におけるＢＰＦ状態情報の変化の量が所定のバンドパスフィルタの切替え速度（フィルタ回転体の回転速度）に応じた変化量となるように、駆動手段５７を制御する。
上記に加えて、トリガ制御回路３０は、レーザ光源ユニット１３に、フラッシュランプ５２の発光を制御するためのフラッシュランプトリガ信号を出力し、フラッシュランプ５２からレーザロッド５１に励起光を照射させる。トリガ制御回路３０は、ＢＰＦ状態信号に基づいて、フラッシュランプトリガ信号を出力する。例えばトリガ制御回路３０は、ＢＰＦ状態情報が、出射すべきパルスレーザ光の波長に対応したバンドパスフィルタが光共振器の光路上に挿入される波長選択手段５６の駆動位置から、レーザロッド５１の励起に要する時間の間に波長選択手段５６が変位する量を差し引いた位置を示す情報になるとフラッシュランプトリガ信号を出力し、レーザロッド５１に対して励起光を照射させる。
トリガ制御回路３０は、励起光の照射後、波長選択手段５６が、出射すべきパルスレーザ光の波長に対応した透過波長のバンドパスフィルタを光共振器の光路上に挿入しているタイミングでＱスイッチ５５にＱスイッチトリガ信号を出力する。例えば波長選択手段５６がフィルタ回転体で構成されるとき、トリガ制御回路３０は、ＢＰＦ状態情報が、出射すべきパルスレーザ光の波長の光を透過するバンドパスフィルタが光共振器の光路上に挿入されていることを示す位置となっているときに、Ｑスイッチトリガ信号を出力する。Ｑスイッチ５５が、Ｑスイッチトリガ信号に応答して光共振器内の挿入損失を損失大から損失小に急激に変化させることで（Ｑスイッチがオンすることで）、出力側のミラー５４からパルスレーザ光が出射する。
トリガ制御回路３０は、Ｑスイッチトリガ信号のタイミング、すなわちパルスレーザ光の出射タイミングに合わせて、ＡＤ変換手段２２にサンプリングトリガ信号（ＡＤトリガ信号）を出力する。ＡＤ変換手段２２は、サンプリングトリガ信号にと基づいて光音響信号のサンプリングを開始する。
図３は、波長選択手段５６、駆動手段５７、及び駆動状態検出手段５８の構成例を示す。この例では、波長選択手段５６は２つのバンドパスフィルタを含むフィルタ回転体であり、駆動手段５７はサーボモータである。また、駆動状態検出手段５８はロータリーエンコーダーである。波長選択手段５６は、サーボモータの出力軸の回転に従って回転する。波長選択手段５６を構成するフィルタ回転体の半分（例えば回転変位位置０°から１８０°）は波長７５０ｎｍの光を透過させる第１のバンドパスフィルタであり、残りの半分（例えば回転変位位置１８０°から３６０°）は波長８００ｎｍの光を透過させる第２のバンドパスフィルタである。このようなフィルタ回転体を回転させることで、光共振器の光路上に、第１のバンドパスフィルタと第２のバンドパスフィルタとを、フィルタ回転体の回転速度に応じた切り替え速度で交互に挿入することができる。
ロータリーエンコーダーは、サーボモータの出力軸に取り付けられたスリット入りの回転板と透過型フォトインタラプタとでフィルタ回転体の回転変位を検出し、ＢＰＦ状態情報を生成する。トリガ制御回路３０は、例えばＢＰＦ状態情報をモニタし、所定時間の間にロータリーエンコーダーで検出されるサーボモータの回転軸の回転変位量が所定の量に保たれるように、ＢＰＦ制御信号を通じてサーボモータに供給する電圧などを制御することで、フィルタ回転体を所定の速度で回転させる。
図４は、光音響画像生成装置１０の動作手順を示す。ここでは、被検体のレーザ光が照射される領域が複数の部分領域に分割されているものとして説明する。トリガ制御回路３０は、被検体に対するパルスレーザ光照射に先立って、所定の回転速度でレーザ光源ユニット１３内の波長選択手段（フィルタ回転体）５６を回転させる旨のＢＰＦ制御信号をレーザ光源ユニット１３に出力する（ステップＡ１）。例えば図３に示すフィルタ回転体を用いる場合で、かつ１秒間に２４回のパルスレーザ光を出射する場合、フィルタ回転体の１回転の間に７５０ｎｍと８００ｎｍの２つのパルスレーザ光が出射可能であるから、１秒あたり２４／２＝１２回転の回転速度でフィルタ回転体を回転させればよい。
トリガ制御回路３０は、光音響信号の受信準備が整うと、１つ目の波長（例えば７５０ｎｍ）のパルスレーザ光を出射させるべく、所定のタイミングでレーザ光源ユニット１３にフラッシュランプトリガ信号を出力する（ステップＡ２）。レーザ光源ユニット１３のフラッシュランプ５２は、フラッシュランプトリガ信号に応答して点灯し、レーザロッド５１の励起が開始される（ステップＡ３）。トリガ制御回路３０は、ＢＰＦ状態情報に基づいて、例えば波長選択手段５６の回転変位位置が波長７５０ｎｍの光を透過させるバンドパスフィルタを光共振器の光路上に挿入している位置となっているタイミングから逆算されたタイミングで、フラッシュランプ５２を点灯させる。
トリガ制御回路３０は、フラッシュランプ５２の点灯後、ＢＰＦ状態情報に基づいて、波長選択手段５６の回転変位位置が波長７５０ｎｍの光を透過するバンドパスフィルタが光共振器の光路上に挿入される位置となっているタイミングでＱスイッチ５５をオンにする（ステップＡ４）。Ｑスイッチ５５がオンになったとき、光共振器の光路上には透過波長７５０ｎｍのバンドパスフィルタが挿入されているため、レーザ光源ユニット１３は、波長７５０ｎｍのパルスレーザ光を出射する。
レーザ光源ユニット１３から出射した波長７５０ｎｍのパルスレーザ光は、例えばプローブ１１まで導光され、プローブ１１から被検体の１つ目の部分領域に照射される。被検体内では、光吸収体が照射されたパルスレーザ光のエネルギーを吸収することで、光音響信号が発生する。プローブ１１は、被検体内で発生した光音響信号を検出する。プローブ１１で検出された光音響信号は、受信回路２１にて受信される。
トリガ制御回路３０は、Ｑスイッチトリガ信号を出力するタイミングに合わせて、ＡＤ変換手段２２にサンプリングトリガ信号を出力する。ＡＤ変換手段２２は、受信回路２１で受信された光音響信号を、所定のサンプリング周期でサンプリングする（ステップＡ５）。ＡＤ変換手段２２でサンプリングされた光音響信号は、受信メモリ２３に第１の光音響データとして格納される。
制御手段３１は、残りの波長があるか、つまり出射すべき複数の波長のパルスレーザ光の全てを出射したか否かを判断する（ステップＡ６）。残りの波長がある場合、次の波長のパルスレーザ光を出射させるために、ステップＡ２に戻り、トリガ制御回路３０からレーザ光源ユニット１３にフラッシュランプトリガ信号を出力する。フラッシュランプ５２は、ステップＡ３でフラッシュランプトリガ信号に応答して点灯し、トリガ制御回路３０は、ステップＡ４で、波長選択手段５６の回転変位位置が、２つ目の波長（８００ｎｍ）の光を透過させるバンドパスフィルタが光共振器の光路上に挿入される位置となっているタイミングでＱスイッチ５５をオンにし、パルスレーザ光を出射させる。
レーザ光源ユニット１３から出射した波長８００ｎｍのパルスレーザ光は、例えばプローブ１１まで導光され、プローブ１１から被検体の１つ目の部分領域に照射される。プローブ１１は、被検体内の光吸収体が波長８００ｎｍのパルスレーザ光を吸収することで発生した光音響信号を検出する。トリガ制御回路３０は、Ｑスイッチトリガ信号の出力に合わせてＡＤ変換手段２２にサンプリングトリガ信号を出力し、ＡＤ変換手段２２は、ステップＡ５で光音響信号のサンプリングを行う。ＡＤ変換手段２２でサンプリングされた光音響信号は、受信メモリ２３に第２の光音響データとして格納される。光音響画像生成装置１０は、被検体に照射すべきパルスレーザ光の各波長に対してステップＡ１からＡ５を実行し、各波長のパルスレーザ光を被検体に照射して、被検体から光音響信号を検出する。
制御手段３１は、ステップＡ６で残りの波長がないと判断すると、部分領域を全て選択したか否かを判断する（ステップＡ７）。選択すべき部分領域が残っておりときは、ステップＡ２に戻る。光音響画像生成装置１０は、各部分領域に対してステップＡ２からＡ６を実行し、各部分領域に各波長（７５０ｎｍ、８００ｎｍ）のパルスレーザ光を順次照射し、各部分領域に対応した第１の光音響データと第２の光音響データとを受信メモリ２３に格納する。全ての部分領域に対してパルスレーザ光の照射及び光音響信号の検出を行うと、１フレームの光音響画像を生成するために必要な光音響データが揃う。
制御手段３１は、ステップＡ７で全ての部分領域を選択したと判断すると、光音響画像の生成に処理を移す。複素数化手段２４は、受信メモリ２３から第１の光音響データと第２の光音響データとを読み出し、第１の光音響画像データを実部とし、第２の光音響画像データを虚部とした複素数データを生成する（ステップＡ８）。光音響画像再構成手段２５は、ステップＡ８で複素数化された複素数データから、フーリエ変換法（ＦＴＡ法）により画像再構成を行う（ステップＡ９）。
位相情報抽出手段２６は、再構成された複素数データ（再構成画像）から位相情報を抽出する（ステップＡ１０）。位相情報抽出手段２６は、例えば再構成された複素数データがＸ＋ｉＹで表わされるとき、θ＝ｔａｎ−１（Ｙ／Ｘ）を位相情報として抽出する（ただし、Ｘ＝０の場合はθ＝９０°）。強度情報抽出手段２７は、再構成された複素数データから強度情報を抽出する（ステップＡ１１）。強度情報抽出手段２７は、例えば再構成された複素数データがＸ＋ｉＹで表わされるとき、（Ｘ２＋Ｙ２）１／２を強度情報として抽出する。
検波・対数変換手段２８は、ステップＡ１１で抽出された強度情報に対して検波・対数変換処理を施す。光音響画像構築手段２９は、ステップＡ１０で抽出された位相情報と、ステップＡ１１で抽出された強度情報に対して検波・対数変換処理を施したものとに基づいて、光音響画像を生成する（ステップＡ１２）。光音響画像構築手段２９は、例えば強度情報に基づいて光吸収体の分布画像における各画素の輝度（階調値）を決定し、位相情報に基づいて各画素の色を決定することで、光音響画像を生成する。生成された光音響画像は、画像表示手段１４に表示される。
本実施形態では、波長選択手段５６を、光共振器の光路上に選択的に挿入される複数のバンドパスフィルタが所定の速度で切り替わるように駆動しつつ、フラッシュランプ５２を点灯してレーザロッド５１を励起し、その励起後、光共振器の光路上に、出射すべきパルスレーザ光の波長の光を透過するバンドパスフィルタが挿入されているときに、Ｑスイッチ５５をオンにする。例えば、透過波長が異なる２つのバンドパスフィルタを含む波長選択手段５６を連続的に駆動することで、２つのバンドパスフィルタを連続的、かつ選択的に光共振器の光路上に挿入することができ、レーザ光源ユニット１３から複数波長のレーザ光を連続的に切り換えて出射することができる。また、トリガ制御回路３０は、波長選択手段５６の駆動状態を把握しており、所望の透過波長のバンドパスフィルタが光共振器の光路上に挿入されているときにＱスイッチ５５をオンすることで、超音波ユニット１２側から、レーザ光源ユニット１３が出射するパルスレーザ光の波長を制御することができる。更に、本実施形態では、超音波ユニット１２がパルスレーザ光の出射タイミングを決めており、光音響信号のサンプリング開始にあたり、レーザ光源ユニット１３からレーザ出射を示す同期信号のような信号を取得する必要がない。
本実施形態では、２つの波長で得られた第１の光音響データと、第２の光音響データとの何れか一方を実部、他方を虚部とした複素数データを生成し、その複素数データからフーリエ変換法により再構成画像を生成している。このようにする場合、第１の光音響データと第２の光音響データとを別々に再構成する場合に比べて、再構成を効率的に行うことができる。複数の波長のパルスレーザ光を照射し、各波長のパルスレーザ光を照射したときの光音響信号（光音響データ）を用いることで、各光吸収体の光吸収特性が波長に応じて異なることを利用した機能イメージングを行うことができる。
また、本実施形態では、例えば光照射領域が３つの部分領域に分かれているときには、第１の部分領域に対して第１の波長のパルスレーザ光、第２の波長のパルスレーザ光を順次に照射し、次いで、第２の部分領域に対して第１の波長のパルスレーザ光、第２の波長のパルスレーザ光を順次に照射し、その後、第３の部分領域に対して第１の波長のパルスレーザ光、第２の波長のパルスレーザ光を順次に照射する。本実施形態では、ある部分領域に対して第１の波長のパルスレーザ光、第２の波長のパルスレーザ光を連続的に照射した後に、次の部分領域に移っている。この場合、第１の波長のパルスレーザ光を３つの部分領域に照射した後に、第２の波長のパルスレーザ光を３つの部分領域に照射する場合に比べて、同じ位置において第１の波長のパルスレーザ光を照射してから第２のパルスレーザ光が照射されるまでの間の時間を短くすることができる。第１の波長のパルスレーザ光が照射されてから第２の波長のパルスレーザ光が照射されるまでの間の時間を短縮することで、第１の光音響データと第２の光音響データとの不整合を抑制することができる。
引き続き、本発明の第２実施形態を説明する。図５は、本発明の第２実施形態の光音響画像生成装置を示す。本実施形態の光音響画像生成装置１０ａは、超音波ユニット１２ａが、図１に示す第１実施形態の光音響画像生成装置１０における超音波ユニット１２の構成に加えて、データ分離手段３２、超音波画像再構成手段３３、検波・対数変換手段３４、超音波画像構築手段３５、画像合成手段３６、及び送信制御回路３７を有する。本実施形態の光音響画像生成装置１０は、光音響画像に加えて、超音波画像を生成する点で第１実施形態と相違する。その他の部分は、第１実施形態と同様でよい。
本実施形態では、プローブ１１は、光音響信号の検出に加えて、被検体に対する超音波の出力（送信）、及び送信した超音波に対する被検体からの反射超音波の検出（受信）を行う。トリガ制御回路３０は、超音波画像の生成時は、送信制御回路３７に超音波送信を指示する旨の超音波送信トリガ信号を送る。送信制御回路３７は、トリガ信号を受けると、プローブ１１から超音波を送信させる。プローブ１１は、超音波の送信後、被検体からの反射超音波を検出する。
プローブ１１が検出した反射超音波は、受信回路２１を介してＡＤ変換手段２２に入力される。トリガ制御回路３０は、超音波送信のタイミングに合わせてＡＤ変換手段２２にサンプリグトリガ信号を送り、反射超音波のサンプリングを開始させる。ＡＤ変換手段２２は、反射超音波のサンプリングデータ（反射超音波データ）を受信メモリ２３に格納する。
データ分離手段３２は、受信メモリ２３に格納された反射超音波データと、第１及び第２の光音響データとを分離する。データ分離手段３２は、反射超音波データを超音波画像再構成手段３３に渡し、第１及び第２の光音響データを複素数化手段２４に渡す。第１及び第２の光音響データに基づく光音響画像の生成は第１実施形態と同様である。データ分離手段３２は、分離した反射超音波のサンプリングデータを、超音波画像再構成手段３３に入力する。
超音波画像再構成手段３３は、プローブ１１の複数の超音波振動子で検出された反射超音波（そのサンプリングデータ）に基づいて、超音波画像の各ラインのデータを生成する。超音波画像再構成手段３３は、例えばプローブ１１の６４個の超音波振動子からのデータを、超音波振動子の位置に応じた遅延時間で加算し、１ライン分のデータを生成する（遅延加算法）。
検波・対数変換手段３４は、超音波画像再構成手段３３が出力する各ラインのデータの包絡線を求め、求めた包絡線を対数変換する。超音波画像構築手段３５は、対数変換が施された各ラインのデータに基づいて、超音波画像を生成する。超音波画像再構成手段３３、検波・対数変換手段３４、及び超音波画像構築手段３５は、反射超音波に基づいて超音波画像を生成する超音波画像生成手段を構成する。
画像合成手段３６は、光音響画像と超音波画像とを合成する。画像合成手段３６は、例えば光音響画像と超音波画像とを重畳することで画像合成を行う。その際、画像合成手段３６は、光音響画像と超音波画像とで、対応点が同一の位置となるように位置合わせをすることが好ましい。合成された画像は、画像表示手段１４に表示される。画像合成を行わずに、画像表示手段１４に、光音響画像と超音波画像とを並べて表示し、或いは光音響画像と超音波画像とを切り替えてすることも可能である。
図６は、光音響画像生成装置１０ａの動作手順を示している。以下では、被検体のレーザ光が照射される領域が複数の部分領域に分割されているものとして説明する。トリガ制御回路３０は、所定の回転速度でレーザ光源ユニット１３内の波長選択手段（フィルタ回転体）５６を回転させる旨のＢＰＦ制御信号をレーザ光源ユニット１３に出力する（ステップＢ１）。
トリガ制御回路３０は、光音響信号の受信準備が整うと、１つ目の波長（例えば７５０ｎｍ）のパルスレーザ光を出射させるべく、フラッシュランプトリガ信号を出力する（ステップＢ２）。フラッシュランプ５２は、フラッシュランプトリガ信号に応答して点灯し、レーザロッド５１の励起が開始される（ステップＢ３）。
トリガ制御回路３０は、フラッシュランプ５２の点灯後、ＢＰＦ状態情報に基づいて、波長選択手段５６の回転変位位置が波長７５０ｎｍの光を透過するバンドパスフィルタが光共振器の光路上に挿入される位置となっているタイミングでＱスイッチ５５をオンにする。Ｑスイッチ５５がオンになることで、レーザ光源ユニット１３は、波長７５０ｎｍのパルスレーザ光を出射する。
レーザ光源ユニット１３から出射した波長７５０ｎｍのパルスレーザ光は、例えばプローブ１１まで導光され、プローブ１１から被検体の１つ目の領域に照射される。被検体内では、光吸収体が照射されたパルスレーザ光のエネルギーを吸収することで、光音響信号が発生する。プローブ１１は、被検体内で発生した光音響信号を検出する。トリガ制御回路３０は、Ｑスイッチトリガ信号の出力に合わせて、ＡＤ変換手段２２にサンプリングトリガ信号を出力する。ＡＤ変換手段２２は、受信回路２１を介してプローブ１１で検出された光音響信号を受信し、所定のサンプリング周期で光音響信号をサンプリングする（ステップＢ５）。ＡＤ変換手段２２でサンプリングされた光音響信号は、受信メモリ２３に第１の光音響データとして格納される。
制御手段３１は、残りの波長があるか、つまり出射すべき複数の波長のパルスレーザ光の全てを出射させたか否かを判断する（ステップＢ６）。残りの波長がある場合、次の波長のパルスレーザ光を出射させるために、ステップＢ２に戻り、トリガ制御回路３０からレーザ光源ユニット１３にフラッシュランプトリガ信号を出力する。フラッシュランプ５２は、ステップＢ３でフラッシュランプトリガ信号に応答して点灯し、トリガ制御回路３０は、ステップＢ４で、波長選択手段５６の回転変位位置が、２つ目の波長（８００ｎｍ）の光を透過させるバンドパスフィルタが光共振器の光路上に挿入される位置となっているタイミングでＱスイッチ５５をオンにし、パルスレーザ光を出射させる。
レーザ光源ユニット１３から出射した波長８００ｎｍのパルスレーザ光は、例えばプローブ１１まで導光され、プローブ１１から被検体の１つ目の部分領域に照射される。プローブ１１は、被検体内の光吸収体が波長８００ｎｍのパルスレーザ光を吸収することで発生した光音響信号を検出する。トリガ制御回路３０は、Ｑスイッチトリガ信号の出力に合わせてＡＤ変換手段２２にサンプリングトリガ信号を出力し、ＡＤ変換手段２２は、ステップＢ５で光音響信号のサンプリングを行う。ＡＤ変換手段２２でサンプリングされた光音響信号は、受信メモリ２３に第２の光音響データとして格納される。光音響画像生成装置１０は、被検体に照射すべきパルスレーザ光の各波長に対してステップＢ１からＢ５を実行し、各波長のパルスレーザ光を被検体に照射して、被検体から光音響信号を検出する。ステップＢ１からＢ５は、図４のステップＡ１からＡ５と同様でよい。
制御手段３１は、ステップＢ６で残りの波長がないと判断すると、超音波の送受信に処理を移す。トリガ制御回路３０は、送信制御回路３７を介してプローブ１１から被検体に超音波を送信させる（ステップＢ７）。ステップＢ７では、被検体のパルスレーザ光が照射された部分領域と同じ領域に対して超音波を送信する。プローブ１１は、送信した超音波に対する反射超音波を検出する（ステップＢ８）。検出された反射超音波は、受信回路２１を経てＡＤ変換手段２２でサンプリングされ、受信メモリ２３に反射超音波データとして格納される。
制御手段３１は、部分領域を全て選択したか否かを判断する（ステップＢ９）。選択すべき部分領域が残っておりときは、ステップＢ２に戻る。光音響画像生成装置１０は、各部分領域に対してステップＢ２からＢ６を実行し、各部分領域に各波長（７５０ｎｍ、８００ｎｍ）のパルスレーザ光を順次照射して、第１の光音響データと第２の光音響データとを受信メモリ２３に格納する。また、ステップＢ７及びＢ８を実行して、反射超音波データを受信メモリ２３に格納する。全ての部分領域に対してパルスレーザ光の照射と光音響信号の検出、及び超音波の送受信を行うと、１フレームの光音響画像及び超音波画像を生成するために必要なデータが揃う。
制御手段３１は、ステップＢ９で全ての部分領域を選択したと判断すると、光音響画像及び超音波画像の生成に処理を移す。データ分離手段３２は、第１及び第２の光音響データと反射超音波データとを分離する。データ分離手段３２は、分離した第１及び第２の光音響データを複素数化手段２４に渡し、反射超音波データを超音波画像再構成手段３３に渡す。複素数化手段２４は、第１の光音響画像データを実部とし、第２の光音響画像データを虚部とした複素数データを生成する（ステップＢ１０）。光音響画像再構成手段２５は、ステップＢ１０で複素数化された複素数データから、フーリエ変換法（ＦＴＡ法）により画像再構成を行う（ステップＢ１１）。
位相情報抽出手段２６は、再構成された複素数データから位相情報を抽出する（ステップＢ１２）。強度情報抽出手段２７は、再構成された複素数データから強度情報を抽出する（ステップＢ１３）。検波・対数変換手段２８は、ステップＢ１３で抽出された強度情報に対して検波・対数変換処理を施す。光音響画像構築手段２９は、ステップＢ１２で抽出された位相情報と、ステップＢ１３で抽出された強度情報に対して検波・対数変換処理を施したものとに基づいて、光音響画像を生成する（ステップＢ１４）。ステップＢ１０からＢ１４は、図４のステップＡ８からＡ１２と同様でよい。
超音波画像再構成手段３３は、例えば遅延加算法により、超音波画像の各ラインのデータを生成する。検波・対数変換手段３４は、超音波画像再構成手段３３が出力する各ラインのデータの包絡線を求め、求めた包絡線を対数変換する。超音波画像構築手段３５は、対数変換が施された各ラインのデータに基づいて、超音波画像を生成する（ステップＢ１５）。画像合成手段３６は、光音響画像と超音波画像とを合成し、合成した画像を画像表示手段１４に表示する（ステップＢ１６）。
本実施形態では、光音響画像生成装置は、光音響画像に加えて超音波画像を生成する。超音波画像を参照することで、光音響画像では画像化することができない部分を観察することができる。その他の効果は、第１実施形態と同様である。
なお、上記各実施形態では、第１の光音響データと第２の光音響データとを複素数化する例について説明したが、複素数化せずに、第１の光音響データと第２の光音響データとを別々に再構成してもよい。さらに、ここでは、複素数化して位相情報を用いて第１の光音響データと第２の光音響データの比を計算しているが、両者の強度情報から比を計算しても同様の効果が得られる。また、強度情報は、第１の再構成画像における信号強度と、第２の再構成画像における信号強度とに基づいて生成できる。
光音響画像の生成に際して、被検体に照射されるパルスレーザ光の波長の数は２つには限られず、３以上のパルスレーザ光を被検体に照射し、各波長に対応する光音響データに基づいて光音響画像を生成してもよい。その場合、例えば位相情報抽出手段２６は、各波長に対応する光音響データ間での相対的な信号強度の大小関係を位相情報として生成すればよい。また、強度情報抽出手段２７は、例えば各波長に対応する光音響データにおける信号強度を１つにまとめたものを強度情報として生成すればよい。
上記各実施形態では、トリガ制御回路３０がＢＰＦ状態情報をモニタし、ＢＰＦ制御信号を通して波長選択手段５６（フィルタ回転体）の回転速度を所定の回転速度に制御するものとして説明したが、これには限定されない。図７は、レーザ光源ユニットの変形例を示している。このレーザ光源ユニット１３ａは、図２に示すレーザ光源ユニット１３の構成に加えて回転制御部５９を有する。回転制御部５９は、所定時間の間に駆動状態検出手段５８が検出した回転変位の量が、フィルタ回転体の所定の回転速度に応じた量になるように駆動手段５７に供給する電圧などを制御する。トリガ制御回路３０は、ＢＰＦ制御信号を通じて、回転制御部５９に、フィルタ回転体の回転速度を指示する。回転制御部５９は、フィルタ回転体の回転速度が指示された回転速度となるように駆動手段５７を駆動する。
上記各実施形態では、主に、波長選択手段５６が、図３に示されるような２つのバンドパスフィルタ領域を含むフィルタ回転体で構成される例を説明したが、波長選択手段５６は、透過波長が異なる複数のバンドパスフィルタを、光共振器の光路上に選択的に挿入歌のであればよく、図３に示したフィルタ回転体の構成には限定されない。例えば、波長選択手段５６を、複数のバンドパスフィルタを円周状に配置した回転体で構成してもよい。波長選択手段５６は回転体である必要はなく、例えば、複数のバンドパスフィルタを一列に並べたものでもよい。その場合、複数のバンドパスフィルタが循環的に光共振器の光路上に挿入されるよう波長選択手段５６を駆動してもよいし、一列に並べられた複数のバンドパスフィルタが光共振器の光路上を横切るように波長選択手段５６を往復駆動させてもよい。
以上、本発明をその好適な実施形態に基づいて説明したが、本発明の超音波ユニット及び光音響画像生成装置は、上記実施形態にのみ限定されるものではなく、上記実施形態の構成から種々の修正及び変更を施したものも、本発明の範囲に含まれる。
１０：光音響画像生成装置
１１：プローブ
１２：超音波ユニット
１３：レーザ光源ユニット
１４：画像表示手段
２１：受信回路
２２：ＡＤ変換手段
２３：受信メモリ
２４：複素数化手段
２５：光音響画像再構成手段
２６：位相情報抽出手段
２７：強度情報抽出手段
２８：検波・対数変換手段
２９：光音響画像構築手段
３０：トリガ制御回路
３１：制御手段
３２：データ分離手段
３３：超音波画像再構成手段
３４：検波・対数変換手段
３５：超音波画像構築手段
３６：画像合成手段
３７：送信制御回路
５１：レーザロッド
５２：フラッシュランプ
５３、５４：ミラー
５５：Ｑスイッチ
５６：波長選択手段
５７：駆動手段
５８：駆動状態検出手段
５９：回転制御部
１１：プローブ
１２：超音波ユニット
１３：レーザ光源ユニット
１４：画像表示手段
２１：受信回路
２２：ＡＤ変換手段
２３：受信メモリ
２４：複素数化手段
２５：光音響画像再構成手段
２６：位相情報抽出手段
２７：強度情報抽出手段
２８：検波・対数変換手段
２９：光音響画像構築手段
３０：トリガ制御回路
３１：制御手段
３２：データ分離手段
３３：超音波画像再構成手段
３４：検波・対数変換手段
３５：超音波画像構築手段
３６：画像合成手段
３７：送信制御回路
５１：レーザロッド
５２：フラッシュランプ
５３、５４：ミラー
５５：Ｑスイッチ
５６：波長選択手段
５７：駆動手段
５８：駆動状態検出手段
５９：回転制御部
Claims (12)
- 相互に異なる複数の波長のパルスレーザ光を出射するレーザ光源ユニットであって、レーザロッドと、該レーザロッドに励起光を照射する励起光源と、前記レーザロッドを挟んで対向する一対のミラーを含む光共振器と、前記光共振器内に挿入されたＱスイッチと、透過波長が相互に異なる複数のバンドパスフィルタを含み、該複数のバンドパスフィルタを前記光共振器の光路上に選択的に挿入する波長選択手段と、前記光共振器の光路上に挿入されるバンドパスフィルタが所定の順序で順次に切り替わるように前記波長選択手段を駆動する駆動手段とを有するレーザ光源ユニットと、
光音響画像の生成を行う音響波ユニットであって、前記複数の波長のパルスレーザ光が被検体に照射されたときに被検体内で生じる光音響信号を検出し、各波長に対応した光音響データを生成する検出手段と、前記各波長に対応した光音響データ間の相対的な信号強度の大小関係を抽出する強度比抽出手段と、前記抽出された大小関係に基づいて光音響画像を生成する光音響画像構築手段と、前記光共振器の光路上に挿入されるバンドパスフィルタが所定の切替え速度で切り替わるように前記駆動手段を制御しつつ、前記励起光源から前記レーザロッドに励起光を照射させ、該励起光の照射後、前記波長選択手段が、出射すべきパルスレーザ光の波長に対応した透過波長のバンドパスフィルタを前記光路上に挿入しているタイミングで前記Ｑスイッチをオンにしてパルスレーザ光を出射させるトリガ制御回路とを有する音響波ユニットとを備えたことを特徴とする光音響画像生成装置。 - 前記波長選択手段が、回転変位に伴って前記光共振器の光路上に選択的に挿入するバンドパスフィルタを切り替えるフィルタ回転体で構成され、前記駆動手段は前記フィルタ回転体を回転駆動する請求項１に記載の光音響画像生成装置。
- 前記トリガ制御回路は、前記フィルタ回転体が、所定の方向に所定の回転速度で連続的に回転するように前記駆動手段を制御する請求項２に記載の光音響画像生成装置。
- 前記所定の回転速度は、出射すべきパルスレーザ光の波長の数と単位時間当たりのパルスレーザ光の発光回数とに基づいて決定される請求項３に記載の光音響画像生成装置。
- 前記トリガ制御回路は、前記フィルタ回転体の回転変位位置を示すフィルタ状態情報に基づいて、前記励起光を照射させるタイミング及び前記Ｑスイッチをオンにするタイミングを決定する請求項２から４何れかに記載の光音響画像生成装置。
- 前記トリガ制御回路は、前記フィルタ状態情報が、前記出射すべきパルスレーザ光の波長に対応したバンドパスフィルタが前記光路上に挿入される前記フィルタ回転体の位置から前記レーザロッドの励起に要する時間の間に前記フィルタ回転体が回転変位する量を差し引いた位置を示す情報になると、前記レーザロッドに励起光を照射させる請求項５に記載の光音響画像生成装置。
- 前記トリガ制御回路は、所定時間の間における前記フィルタ状態情報の変化の量が前記所定の回転速度に応じた変化量となるように、前記駆動手段を制御する請求項５又は６に記載の光音響画像生成装置。
- 前記レーザ光源ユニットは、所定時間の間に前記フィルタ回転体の回転変位の量が前記所定の回転速度に応じた量になるように前記駆動手段を制御する回転制御部を更に有し、前記トリガ制御回路は、前記回転制御部を通じて前記駆動手段を制御する請求項２から６何れかに記載の光音響画像生成装置。
- 前記音響波ユニットは、前記各波長に対応した光音響データに基づいて信号強度を示す強度情報を生成する強度情報抽出手段を更に備え、
前記光音響画像構築手段が、前記光音響画像の各画素の階調値を前記強度情報に基づいて決定すると共に、各画素の表示色を前記抽出された大小関係に基づいて決定する請求項１から８何れかに記載の光音響画像生成装置。 - 前記レーザ光源ユニットが出射すべきパルスレーザ光の複数の波長が第１の波長と第２の波長を含み、
前記音響波ユニットは、前記第１の波長のパルスレーザ光が照射されたときに検出される光音響信号に対応する第１の光音響データと、前記第２の波長のパルスレーザ光が照射されたときに検出される光音響信号に対応する第２の光音響データとのうちの何れか一方を実部、他方を虚部とした複素数データを生成する複素数化手段と、前記複素数データからフーリエ変換法により再構成画像を生成する光音響画像再構成手段とを更に備え、
前記強度比抽出手段が前記再構成画像から前記大小関係としての位相情報を抽出し、前記強度情報抽出手段が前記再構成画像から前記強度情報を抽出する請求項９に記載の光音響画像生成装置。 - 前記検出手段は、更に、被検体に送信された音響波に対する反射音響波を検出して反射音響波データを生成し、
前記音響波ユニットは、前記反射音響波データに基づいて音響波画像を生成する音響波画像生成手段を更に備える請求項１から１０何れかに記載の光音響画像生成装置。 - 相互に異なる複数の波長のパルスレーザ光が被検体に照射されたときに被検体内で生じた光音響信号を検出し、各波長に対応した光音響データを生成する検出手段と、
前記各波長に対応した光音響データ間の相対的な信号強度の大小関係を抽出する強度比抽出手段と、
前記抽出された大小関係に基づいて光音響画像を生成する光音響画像構築手段と、
波長選択手段を駆動する駆動手段を制御し、かつレーザロッドを挟んで対向する一対のミラーを含む光共振器内に挿入されたＱスイッチを制御するトリガ制御回路とを備え、
前記波長選択手段は、透過波長が相互に異なる複数のバンドパスフィルタを含むものであり、かつ前記駆動手段は、前記光共振器の光路上に複数のバンドパスフィルタが所定の順序で選択的に挿入されるように前記波長選択手段を駆動するものであり、
前記トリガ制御回路は、前記光共振器の光路上に挿入されるバンドパスフィルタが所定の切替え速度で切り替わるように前記駆動手段を制御しつつ、励起光源から前記レーザロッドに励起光を照射させ、該励起光の照射後、前記波長選択手段が、出射すべきパルスレーザ光の波長に対応した透過波長のバンドパスフィルタを前記光路上に挿入しているタイミングで前記Ｑスイッチをオンにしてパルスレーザ光を出射させることを特徴とする音響波ユニット。
Priority Applications (5)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2012160427A JP5713968B2 (ja) | 2011-07-29 | 2012-07-19 | 光音響画像生成装置及び音響波ユニット |
| EP12820334.6A EP2737855B1 (en) | 2011-07-29 | 2012-07-26 | Photoacoustic image-generating apparatus and acoustic unit |
| CN201280037835.3A CN103717141B (zh) | 2011-07-29 | 2012-07-26 | 光声图像生成装置及声波单元 |
| PCT/JP2012/068904 WO2013018627A1 (ja) | 2011-07-29 | 2012-07-26 | 光音響画像生成装置及び音響波ユニット |
| US14/167,808 US9486144B2 (en) | 2011-07-29 | 2014-01-29 | Photoacoustic image generation apparatus and acoustic wave unit |
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2011166978 | 2011-07-29 | ||
| JP2011166978 | 2011-07-29 | ||
| JP2012160427A JP5713968B2 (ja) | 2011-07-29 | 2012-07-19 | 光音響画像生成装置及び音響波ユニット |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2013051405A JP2013051405A (ja) | 2013-03-14 |
| JP5713968B2 true JP5713968B2 (ja) | 2015-05-07 |
Family
ID=47629152
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2012160427A Active JP5713968B2 (ja) | 2011-07-29 | 2012-07-19 | 光音響画像生成装置及び音響波ユニット |
Country Status (5)
| Country | Link |
|---|---|
| US (1) | US9486144B2 (ja) |
| EP (1) | EP2737855B1 (ja) |
| JP (1) | JP5713968B2 (ja) |
| CN (1) | CN103717141B (ja) |
| WO (1) | WO2013018627A1 (ja) |
Families Citing this family (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5662973B2 (ja) * | 2011-07-29 | 2015-02-04 | 富士フイルム株式会社 | レーザ光源ユニット、その制御方法、光音響画像生成装置及び方法 |
| JP6223014B2 (ja) * | 2013-06-26 | 2017-11-01 | キヤノン株式会社 | 被検体情報取得装置 |
| JP6041832B2 (ja) * | 2013-08-02 | 2016-12-14 | 富士フイルム株式会社 | 光音響画像生成装置及びその作動方法 |
| EP3167812B1 (en) * | 2014-07-08 | 2023-05-31 | FUJIFILM Corporation | Photoacoustic image generation device and insert |
| CN104505702A (zh) * | 2014-11-26 | 2015-04-08 | 江苏中科四象激光科技有限公司 | 一种实现激光输出功率缓升缓降控制的方法 |
| US20190021639A1 (en) * | 2016-01-21 | 2019-01-24 | Sunnybrook Research Institute | Methemoglobin detection using photoacoustic imaging |
| CN108601583B (zh) * | 2016-02-08 | 2021-05-11 | 富士胶片株式会社 | 声波图像生成装置及声波图像生成方法 |
| TWI648035B (zh) * | 2018-03-22 | 2019-01-21 | 亞洲大學 | 最小割脊椎分節演算法 |
Family Cites Families (20)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4272733A (en) * | 1978-10-20 | 1981-06-09 | Allied Chemical Corporation | Broadly tunable chromium-doped beryllium aluminate lasers and operation thereof |
| US4809283A (en) * | 1988-02-26 | 1989-02-28 | Allied-Signal Inc. | Method of manufacturing chromium-doped beryllium aluminate laser rod and lasers incorporating the rods therein |
| US5175736A (en) * | 1988-10-21 | 1992-12-29 | Spectra-Physics | Tunable dye laser with thin birefringent filter for improved tuning |
| JPH1065260A (ja) * | 1996-08-23 | 1998-03-06 | Mitsubishi Cable Ind Ltd | 固体レーザ装置 |
| WO2003028177A1 (en) * | 2001-09-24 | 2003-04-03 | Giga Tera Ag | Pulse-generating laser |
| US20030072333A1 (en) * | 2001-10-11 | 2003-04-17 | International Business Machines Corporation | NRZ to RZ conversion in mode-locked laser using wavelength locked loops |
| JP4406226B2 (ja) | 2003-07-02 | 2010-01-27 | 株式会社東芝 | 生体情報映像装置 |
| JP4643153B2 (ja) * | 2004-02-06 | 2011-03-02 | 株式会社東芝 | 非侵襲生体情報映像装置 |
| US7852961B2 (en) * | 2004-05-20 | 2010-12-14 | Samsung Electronics Co., Ltd. | Digital broadcasting transmission/reception devices capable of improving a receiving performance and signal processing method thereof |
| JP2007081065A (ja) | 2005-09-13 | 2007-03-29 | Sunx Ltd | レーザ装置、レーザ加工装置及び測定装置 |
| CN2915036Y (zh) * | 2005-12-16 | 2007-06-27 | 华南师范大学 | 基于声透镜的光声成像和层析成像装置 |
| EP2275022B8 (en) * | 2006-01-20 | 2016-10-19 | Olympus Corporation | Apparatus for analyzing characteristic information of object with the use of mutual interaction between ultrasound wave and light |
| US9561017B2 (en) * | 2006-12-19 | 2017-02-07 | Koninklijke Philips N.V. | Combined photoacoustic and ultrasound imaging system |
| JP2009142318A (ja) * | 2007-12-11 | 2009-07-02 | Nippon Telegr & Teleph Corp <Ntt> | 成分濃度測定装置 |
| JP4829934B2 (ja) * | 2008-07-11 | 2011-12-07 | キヤノン株式会社 | 検査装置 |
| JP5460000B2 (ja) | 2008-08-20 | 2014-04-02 | キヤノン株式会社 | イメージング装置およびイメージング方法 |
| JP5692988B2 (ja) * | 2009-10-19 | 2015-04-01 | キヤノン株式会社 | 音響波測定装置 |
| CN102095685B (zh) | 2010-12-02 | 2012-10-03 | 华南师范大学 | 基于光谱编码的光声组分解析成像方法及装置 |
| JP5662973B2 (ja) * | 2011-07-29 | 2015-02-04 | 富士フイルム株式会社 | レーザ光源ユニット、その制御方法、光音響画像生成装置及び方法 |
| JP5681675B2 (ja) * | 2011-07-29 | 2015-03-11 | 富士フイルム株式会社 | 光音響画像生成装置及び音響波ユニット |
-
2012
- 2012-07-19 JP JP2012160427A patent/JP5713968B2/ja active Active
- 2012-07-26 CN CN201280037835.3A patent/CN103717141B/zh active Active
- 2012-07-26 EP EP12820334.6A patent/EP2737855B1/en active Active
- 2012-07-26 WO PCT/JP2012/068904 patent/WO2013018627A1/ja active Application Filing
-
2014
- 2014-01-29 US US14/167,808 patent/US9486144B2/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| EP2737855B1 (en) | 2017-08-23 |
| CN103717141B (zh) | 2015-06-24 |
| CN103717141A (zh) | 2014-04-09 |
| US9486144B2 (en) | 2016-11-08 |
| US20140148682A1 (en) | 2014-05-29 |
| EP2737855A4 (en) | 2015-04-01 |
| WO2013018627A1 (ja) | 2013-02-07 |
| EP2737855A1 (en) | 2014-06-04 |
| JP2013051405A (ja) | 2013-03-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP5713968B2 (ja) | 光音響画像生成装置及び音響波ユニット | |
| JP5796896B2 (ja) | 断層画像生成装置及び方法 | |
| JP6010306B2 (ja) | 光音響計測装置 | |
| JP5730253B2 (ja) | レーザ光源ユニット及び光音響画像生成装置 | |
| JP5662973B2 (ja) | レーザ光源ユニット、その制御方法、光音響画像生成装置及び方法 | |
| JP5681675B2 (ja) | 光音響画像生成装置及び音響波ユニット | |
| WO2012114729A1 (ja) | 光音響画像生成装置、及び方法 | |
| JP5662974B2 (ja) | レーザ光源ユニット、その制御方法、光音響画像生成装置及び方法 | |
| JP5694991B2 (ja) | 光音響画像化方法および装置 | |
| WO2013132976A1 (ja) | レーザ装置及び光音響計測装置 | |
| WO2013129108A1 (ja) | レーザ光源ユニット、その制御方法および光音響画像生成装置 | |
| JP2013214722A (ja) | レーザ装置及び光音響計測装置 | |
| US20140316239A1 (en) | Laser light source unit, method of controlling the same, and photoacoustic image generation device | |
| WO2013129105A1 (ja) | レーザ光源ユニット、その制御方法および光音響画像生成装置 | |
| JP2014064733A (ja) | 光音響計測装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20140116 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20150106 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20150213 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20150303 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20150310 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 5713968 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |