(57)【特許請求の範囲】
【請求項１】
相互に異なる２以上の波長を含む所定の波長系列で複数のパルスレーザ光を順次に出射するレーザ光源ユニットであって、
レーザロッドと、
前記レーザロッドに励起光を照射する励起光源と、
前記レーザロッドを挟んで対向する一対のミラーを含む光共振器と、
前記光共振器内に挿入されたＱスイッチと、
前記光共振器の内部に挿入され、回転変位に伴って前記光共振器の発振波長を変化させる複屈折フィルタと、
前記波長系列に含まれる波長の数に応じた所定の回転速度で前記複屈折フィルタを回転させる回転制御部と、
前記励起光源から前記レーザロッドに励起光を照射した後、前記複屈折フィルタの回転変位位置が、出射すべきパルスレーザ光の波長に対応した位置となるタイミングで前記Ｑスイッチをオンにしてパルスレーザ光を出射させる発光制御部とを備えたことを特徴とするレーザ光源ユニット。
【請求項２】
前記所定の回転速度は、前記複屈折フィルタにおける回転変位位置に対する発振波長の変化特性と、前記波長系列に含まれる波長の数、及び単位時間当たりのパルスレーザ光の発光回数とに基づいて決定される請求項１に記載のレーザ光源ユニット。
【請求項３】
前記複屈折フィルタが１回転の間に繰り返す自由スペクトルレンジの回数をｋ［回／回転］とし、前記波長系列に含まれる波長の数をｎ［個］とし、単位時間当たりのパルスレーザ光の発光回数をｍ［回／秒］とする場合、前記所定の回転速度が、ｖ＝ｍ／（ｋ×ｎ）［回転／秒］で計算された値に決定される請求項１又は２に記載のレーザ光源ユニット。
【請求項４】
前記回転制御部は、前記複屈折フィルタを所定の方向に前記所定の回転速度で連続的に回転させる請求項１から３何れかに記載のレーザ光源ユニット。
【請求項５】
前記複屈折フィルタを回転させる駆動手段と、前記複屈折フィルタの回転変位を検出する回転変位検出手段とを更に備え、
前記回転制御部は、所定時間の間に前記回転変位検出手段が検出した回転変位の量が前記所定の回転速度に応じた量になるように前記駆動手段を制御する請求項１から４何れかに記載のレーザ光源ユニット。
【請求項６】
前記複屈折フィルタはステッピングモータにより回転されており、前記回転制御部が、前記複屈折フィルタが前記所定の回転速度で回転するように、前記ステッピングモータに供給するパルス信号のパルス間隔を制御する請求項１から４何れかに記載のレーザ光源ユニット。
【請求項７】
相互に異なる２以上の波長を含む所定の波長系列で複数のパルスレーザ光を順次に出射するレーザ光源ユニットと、
前記所定の波長系列に含まれる各波長のパルスレーザ光が被検体に照射されたときに被検体内で生じた光音響信号を検出し、各波長に対応した光音響データを生成する検出手段と、
前記各波長に対応した光音響データ間の相対的な信号強度の大小関係を抽出する強度比抽出手段と、
前記抽出された大小関係に基づいて光音響画像を生成する光音響画像構築手段とを備え、
前記レーザ光源ユニットが、
レーザロッドと、
前記レーザロッドに励起光を照射する励起光源と、
前記レーザロッドを挟んで対向する一対のミラーを含む光共振器と、
前記光共振器内に挿入されたＱスイッチと、
前記光共振器の内部に挿入され、回転変位に伴って前記光共振器の発振波長を変化させる複屈折フィルタと、
前記波長系列に含まれる波長の数に応じた所定の回転速度で前記複屈折フィルタを回転させる回転制御部と、
前記励起光源から前記レーザロッドに励起光を照射した後、前記複屈折フィルタの回転変位位置が、出射すべきパルスレーザ光の波長に対応した位置となるタイミングで前記Ｑスイッチをオンにしてパルスレーザ光を出射させる発光制御部とを有することを特徴とする光音響画像生成装置。
【請求項８】
前記各波長に対応した光音響データに基づいて信号強度を示す強度情報を生成する強度情報抽出手段を更に備え、
前記光音響画像構築手段が、前記光音響画像の各画素の階調値を前記強度情報に基づいて決定すると共に、各画素の表示色を前記抽出された大小関係に基づいて決定する請求項７に記載の光音響画像生成装置。
【請求項９】
前記所定の波長系列が第１の波長と第２の波長を含み、
前記第１の波長のパルスレーザ光が照射されたときに検出された光音響信号に対応する第１の光音響データと、前記第２の波長のパルスレーザ光が照射されたときに検出された光音響信号に対応する第２の光音響データとのうちの何れか一方を実部、他方を虚部とした複素数データを生成する複素数化手段と、
前記複素数データからフーリエ変換法により再構成画像を生成する光音響画像再構成手段とを更に備え、
前記強度比抽出手段が前記再構成画像から前記大小関係を抽出し、前記強度情報抽出手段が前記再構成画像から前記強度情報を抽出する請求項８に記載の光音響画像生成装置。
【請求項１０】
前記検出手段が、更に、被検体に送信された音響波に対する反射音響波を検出して反射音響波データを生成し、
前記反射音響波データに基づいて音響波画像を生成する音響波画像生成手段を更に備える請求項７から９何れかに記載の光音響画像生成装置。
【請求項１１】
相互に異なる２以上の波長を含む所定の波長系列で複数のパルスレーザ光を順次に出射するレーザ光源ユニットの制御方法であって、
レーザロッドを挟んで対向する一対のミラーを含む光共振器内に挿入され、回転変位に伴って前記光共振器の発振波長を変化させる複屈折フィルタを所定の回転速度で回転させるステップと、
前記レーザロッドに励起光を照射し、励起光の照射後、前記複屈折フィルタの回転変位位置が、出射すべきパルスレーザ光の波長に対応した位置となるタイミングで前記光共振器内に挿入されたＱスイッチをオンにしてパルスレーザ光を出射させるステップとを有することを特徴とするレーザ光源ユニットの制御方法。
【請求項１２】
相互に異なる２以上の波長を含む所定の波長系列で複数のパルスレーザ光を順次に出射させるステップであって、レーザロッドを挟んで対向する一対のミラーを含む光共振器内に挿入され、回転変位に伴って前記光共振器の発振波長を変化させる複屈折フィルタを所定の回転速度で回転させつつ、前記レーザロッドに励起光を照射し、励起光の照射後、前記複屈折フィルタの回転変位位置が、出射すべきパルスレーザ光の波長に対応した位置となるタイミングで前記光共振器内に挿入されたＱスイッチをオンにしてパルスレーザ光を出射させるステップと、
前記所定の波長系列に含まれる各波長のパルスレーザ光が被検体に照射されたときに被検体内で生じた光音響信号を検出し、各波長に対応した光音響データを生成するステップと、
前記各波長に対応した光音響データ間の相対的な信号強度の大小関係を抽出するステップと、
前記抽出された大小関係に基づいて光音響画像を生成するステップとを有することを特徴とする光音響画像生成方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、レーザ光源ユニット及びその制御方法に関し、更に詳しくは、複数の波長のレーザ光を切り替えて出射可能なレーザ光源ユニット、及びその制御方法に関する。
【０００２】
また、本発明は、光音響画像生成装置及び方法に関し、更に詳しくは、被検体に複数の波長のレーザ光を照射して光音響信号を検出し、検出された光音響信号に基づいて光音響画像を生成する光音響画像生成装置及び方法に関する。
【背景技術】
【０００３】
従来、例えば特許文献１や非特許文献１に示されているように、光音響効果を利用して生体の内部を画像化する光音響画像化装置が知られている。この光音響画像化装置においては、例えばパルスレーザ光等のパルス光が生体に照射される。このパルス光の照射を受けた生体内部では、パルス光のエネルギーを吸収した生体組織が熱によって体積膨張し、音響波が発生する。この音響波を超音波プローブなどで検出し、検出された信号（光音響信号）に基づいて生体内部を可視像化することが可能となっている。光音響画像化方法では、特定の光吸収体において音響波が発生するため、生体における特定の組織、例えば血管等を画像化することができる。
【０００４】
ところで、生体組織の多くは光吸収特性が光の波長に応じて変わり、また一般に、その光吸収特性も組織ごとに特有のものとなっている。例えば図１１に、ヒトの動脈に多く含まれる酸素化ヘモグロビン(酸素と結合したヘモグロビン:oxy-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  等に記載されている従来公知の方法を適用することができる。光音響画像再構成手段２５は、再構成画像を示すフーリエ変換のデータを位相情報抽出手段２６と強度情報抽出手段２７とに入力する。
【００３３】
位相情報抽出手段２６は、強度比抽出手段であり、各波長に対応した光音響データ間の相対的な信号強度の大小関係を抽出する。本実施形態では、位相情報抽出手段２６は、光音響画像再構成手段２５で再構成された再構成画像を入力データとし、複素数データである入力データから、実部と虚部とを比較したときに、相対的に、どちらがどれくらい大きいかを示す位相情報を生成する。位相情報抽出手段２６は、例えば複素数データがＸ＋ｉＹで表わされるとき、θ＝ｔａｎ−１（Ｙ／Ｘ）を位相情報として生成する。なお、Ｘ＝０の場合はθ＝９０°とする。実部を構成する第１の光音響データ（Ｘ）と虚部を構成する第２の光音響データ（Ｙ）とが等しいとき、位相情報はθ＝４５°となる。位相情報は、相対的に第１の光音響データが大きいほどθ＝０°に近づいていき、第２の光音響データが大きいほどθ＝９０°に近づいていく。
【００３４】
強度情報抽出手段２７は、各波長に対応した光音響データに基づいて信号強度を示す強度情報を生成する。本実施形態では、強度情報抽出手段２７は、光音響画像再構成手段２５で再構成された再構成画像を入力データとし、複素数データである入力データから、強度情報を生成する。強度情報抽出手段２７は、例えば複素数データがＸ＋ｉＹで表わされるとき、（Ｘ２＋Ｙ２）１／２を、強度情報として抽出する。検波・対数変換手段２８は、強度情報抽出手段２７で抽出された強度情報を示すデータの包絡線を生成し、次いでその包絡線を対数変換してダイナミックレンジを広げる。
【００３５】
光音響画像構築手段２９は、位相情報抽出手段２６から位相情報を入力し、検波・対数変換手段２８から検波・対数変換処理後の強度情報を入力する。光音響画像構築手段２９は、入力された位相情報と強度情報とに基づいて、光吸収体の分布画像である光音響画像を生成する。光音響画像構築手段２９は、例えば入力された強度情報に基づいて、光吸収体の分布画像における各画素の輝度（階調値）を決定する。また、光音響画像構築手段２９は、例えば位相情報に基づいて、光吸収体の分布画像における各画素の色（表示色）を決定する。光音響画像構築手段２９は、例えば例えば位相０°から９０°の範囲を所定の色に対応させたカラーマップに用いて、入力された位相情報に基づいて各画素の色を決定する。
【００３６】
ここで、位相０°から４５°の範囲は、第１の光音響データが第２の光音響データよりも大きい範囲であるため、光音響信号の発生源は、波長７９８ｎｍに対する吸収よりも波長７５６ｎｍに対する吸収の方が大きい脱酸素化ヘモグロビンを主に含む血液が流れている静脈であると考えられる。一方、位相４５°から９０°の範囲は、第２の光音響データが第１の光音響データよりも小さい範囲であるため、光音響信号の発生源は、波長７９８ｎｍに対する吸収よりも波長７５６ｎｍに対する吸収の方が小さい酸素化ヘモグロビンを主に含む血液が流れている動脈であると考えられる。
【００３７】
そこで、カラーマップとして、例えば位相が０°が青色で、位相が４５°に近づくに連れて白色になるように色が徐々に変化すると共に、位相９０°が赤色で、位相が４５°に近づくに連れて白色になるように色が徐々に変化するようなカラーマップを用いる。この場合、光音響画像上で、動脈に対応した部分を赤色で表わし、静脈に対応した部分を青色で表わすことができる。強度情報を用いずに、位相情報に従って動脈に対応した部分と静脈に対応した部分との色分けを行うだけでもよい。画像表示手段１４は、光音響画像構築手段２９が生成した光音響画像を、表示画面上に表示する。
【００３８】
次いで、レーザ光源ユニット１３の構成を詳細に説明する。図２は、レーザ光源ユニット１３の構成を示す。レーザ光源ユニット１３は、レーザロッド５１、フラッシュランプ５２、ミラー５３、５４、Ｑスイッチ５５、複屈折フィルタ（ＢＲＦ：birefringent filter）５６、駆動手段５７、回転変位検出手段５８、及び制御部５９を有する。レーザロッド５１は、レーザ媒質である。レーザロッド５１には、例えばアレキサンドライト結晶やＣｒ：ＬｉＳＡＦ（Ｃｒ：ＬｉＳｒＡｌＦ６），Ｃｒ：ＬｉＣＡＦ（Ｃｒ：ＬｉＣａＡｌＦ６）結晶，Ｔｉ：Ｓａｐｐｈｉｒｅ結晶を用いることができる。フラッシュランプ５２は、励起光源であり、レーザロッド５１に励起光を照射する。フラッシュランプ５２以外の光源、例えば半導体レーザや固体レーザを、励起光源として用いてもよい。
【００３９】
ミラー５３、５４は、レーザロッド５１を挟んで対向しており、ミラー５３、５４により光共振器が構成される。ミラー５４が出力側であるものとする。光共振器内には、Ｑスイッチ５５及び複屈折フィルタ５６が挿入される。Ｑスイッチ５５により、光共振器内の挿入損失を損失大（低Ｑ）から損失小（高Ｑ）へと急速に変化させることで、パルスレーザ光を得ることができる。複屈折フィルタ５６は、回転変位に伴って透過波長を変化させ、光共振器の発振波長を変化させる。駆動手段５７は、複屈折フィルタ５６を回転させる。回転変位検出手段５８は、複屈折フィルタ５６の回転変位を検出する。
【００４０】
制御部５９は、回転制御部６０と発光制御部６１とを含む。回転制御部６０は、出射すべきパルスレーザ光の波長系列に含まれる波長の数に応じた所定の回転速度で複屈折フィルタを連続的に回転させる。複屈折フィルタの回転速度は、例えば複屈折フィルタ５６における回転変位位置に対する発振波長の変化特性と、波長系列に含まれる波長の数、及び単位時間当たりのパルスレーザ光の発光回数（パルスレーザ光間の時間間隔）とに基づいて決定できる。
【００４１】
発光制御部６１は、フラッシュランプトリガ信号をフラッシュランプ５２に出力し、フラッシュランプ５２からレーザロッド５１に励起光を照射させる。励起光の照射後、発光制御部６１は、複屈折フィルタ５６の回転変位位置が、出射すべきパルスレーザ光の波長に対応した位置となるタイミングでＱスイッチ５５にＱスイッチトリガ信号を出力する。また、発光制御部６１は、Ｑスイッチトリガ信号の出力タイミングに合わせて、Ｑスイッチ同期信号を超音波ユニット１２に出力する。Ｑスイッチ５５が、Ｑスイッチトリガ信号に応答して光共振器内の挿入損失を損失大から損失小に急激に変化させることで（Ｑスイッチがオンすることで）、出力側のミラー５４からパルスレーザ光が出射する。
【００４２】
図１に戻り、制御手段３１は、超音波ユニット１２内の各部の制御を行う。トリガ制御回路３０は、レーザ光源ユニット１３に対して、複屈折フィルタ５６（図２）の回転速度を制御するためのＢＲＦ制御信号を出力する。また、トリガ制御回路３０は、レーザ光源ユニット１３に、フラッシュランプ５２の発光を制御するためのフラッシュランプスタンバイ信号を出力する。トリガ制御回路３０は、例えばレーザ光源ユニット１３の回転制御部６０から複屈折フィルタ５６の現在の回転変位位置を受け取っており、受け取った回転変位位置に基づくタイミングでフラッシュランプスタンバイ信号を出力する。
【００４３】
トリガ制御回路３０は、レーザ光源ユニット１３から、Ｑスイッチがオンになるタイミング、すなわちレーザ発光タイミングを示すＱスイッチ同期信号を入力する。トリガ制御回路３０は、Ｑスイッチ同期信号を受け取ると、ＡＤ変換手段２２にサンプリングトリガ信号（ＡＤトリガ信号）を出力する。ＡＤ変換手段２２は、サンプリングトリガ信号にと基づいて光音響信号のサンプリングを開始する。
【００４４】
図３は、複屈折フィルタ５６、駆動手段５７、及び回転変位検出手段５８の構成例を示す。この例では、駆動手段５７はサーボモータであり、回転変位検出手段５８はロータリーエンコーダーである。複屈折フィルタ５６は、サーボモータの出力軸の回転に従って回転する。ロータリーエンコーダーは、サーボモータの出力軸に取り付けられたスリット入りの回転板と、透過型フォトインタラプタとで、複屈折フィルタ５６の回転変位を検出する。回転制御部６０（図２）は、例えば所定時間の間にロータリーエンコーダーで検出されるサーボモータの回転軸の回転変位量が所定の量に保たれるように、サーボモータに供給する電圧などを制御することで、複屈折フィルタ５６を所定の速度で回転させる。
【００４５】
図４は、複屈折フィルタ５６の回転変位に対する波長透過特性（発振波長特性）の一例を示す。複屈折フィルタ５６は、例えば７００ｎｍから８４０ｎｍまで間で光共振器の発振波長を変化させる。複屈折フィルタ５６は、例えば回転変位位置０°から９０°までの間に（１／４回転で）自由スペクトルレンジ（ＦＳＲ:Free Spectral Range）を３回繰り返し、１回転につきＦＳＲを１２回繰り返す。
【００４６】
図５は、上記の複屈折フィルタ５６を１秒間に１回転の速度で回転したときの発振波長特性を示す。図４に示す波長透過特性を持つ複屈折フィルタ５６を１秒あたり１回転の回転速度で回転させると、複屈折フィルタ５６は、１／４秒の間にＦＳＲを３回繰り返し、１秒間に１２回（１２Ｈｚ）でＦＳＲを繰り返す。複屈折フィルタ５６の回転速度を速めるほど、１秒間におけるＦＳＲの繰り返し回数は多くなり、回転速度を遅くすると、１秒間におけるＦＳＲの繰り返し回数は少なくなる。
【００４７】
図６は、各種トリガと発光タイミングとを示すタイミングチャートである。（ａ）は、時間変化に対する光共振器の発振波長特性（複屈折フィルタ５６の透過波長特性）を示す。（ｂ）はフラッシュランプトリガを示し、（ｃ）はＱスイッチトリガを示す。（ｄ）は、フラッシュランプの発光タイミング及びパルスレーザ光の出射タイミングを示す。なお、図５では、説明簡略化のために、フラッシュランプ５２及びＱスイッチ５５がトリガに対して瞬時に応答するものとしているが、実際には遅延時間が存在する。だだし、その遅延は数μ秒から１００μ秒程度であるので、遅延は無視できる。
【００４８】
超音波ユニット１２のトリガ制御回路３０は、まず、レーザ光源ユニット１３から波長７５０ｎｍのパルスレーザ光を出射させるために、フラッシュランプスタンバイ信号を出力する。レーザ光源ユニット１３の発光制御部６１は、フラッシュランプスタンバイ信号を受け取ると、時刻ｔ１で、フラッシュランプ５２にフラッシュランプトリガ信号を出力し（ｂ）、フラッシュランプ５２を点灯させる（ｄ）。その後、発光制御部６１は、複屈折フィルタ５６の回転変位位置が波長７５０ｎｍに対応した位置となる時刻ｔ２でＱスイッチトリガ信号を出力し（ｃ）、Ｑスイッチ５５をオンにすることで光共振器から波長７５０ｎｍのパルスレーザ光を出射させる。
【００４９】
次いで、超音波ユニット１２のトリガ制御回路３０は、レーザ光源ユニット１３から波長８００ｎｍのパルスレーザ光を出射させるために、フラッシュランプスタンバイ信号を出力する。発光制御部６１は、フラッシュランプスタンバイ信号を受け取ると、時刻ｔ３で、フラッシュランプ５２にフラッシュランプトリガ信号を出力し（ｂ）、フラッシュランプ５２を点灯させる（ｄ）。その後、発光制御部６０は、複屈折フィルタ５６の回転変位位置が波長８００ｎｍに対応した位置となる時刻ｔ４でＱスイッチトリガ信号を出力し（ｃ）、Ｑスイッチ５５をオンにすることで光共振器から８００ｎｍのパルスレーザ光を出射させる。
【００５０】
以降同様に、発光制御部６１は、時刻ｔ５、ｔ７、ｔ９、ｔ１１で、フラッシュランプ５２にフラッシュランプトリガ信号を出力する。また、時刻ｔ６、ｔ８、ｔ１０、ｔ１２で、Ｑスイッチ５５にＱスイッチトリガ信号を出力し、各時刻における複屈折フィルタ５６の透過波長に応じた波長のパルスレーザ光を出射させる。時刻ｔ６及びｔ１０における複屈折フィルタ５６の透過波長は７５０ｎｍであり、時刻ｔ８及びｔ１２における複屈折フィルタの透過波長は８００ｎｍであるので、レーザ光源ユニット１３は、７５０ｎｍと８００ｎｍのパルスレーザ光を順次にこの順で繰り返し出射する。
【００５１】
図６の例では、レーザ光源ユニット１３は、１／１２秒の間に、波長７５０ｎｍのパルスレーザ光と波長８００ｎｍのパルスレーザ光との２つを交互に出射する。レーザ光源ユニット１３は、２つの波長を切り替えつつパルスレーザ光を１秒間に２４回出射する（２４Ｈｚ動作）。言い換えると、１セット２波長のパルスレーザ光を、１秒間に１２セット出射する。
【００５２】
図７は、光音響画像生成装置１０の動作手順を示す。ここでは、被検体のレーザ光が照射される領域が複数の部分領域に分割されているものとして説明する。超音波ユニット１２のトリガ制御回路３０は、被検体に対するパルスレーザ光照射に先立って、所定の回転速度でレーザ光源ユニット１３内の複屈折フィルタ５６を回転させる旨のＢＲＦ制御信号をレーザ光源ユニット１３に出力する（ステップＡ１）。例えば複屈折フィルタ５６が１回転の間にＦＳＲを１２回繰り返すものであり、１／１２秒の間に波長７５０ｎｍのパルスレーザ光と波長８００ｎｍのパルスレーザ光とを順次に出射させる場合（２４Ｈｚ動作の場合）、トリガ制御回路３０は、複屈折フィルタ５６を１秒間で１回転させる旨のＢＲＦ制御信号を出力する。
【００５３】
トリガ制御回路３０は、光音響信号の受信準備が整うと、波長系列を構成する１つ目の波長（７５０ｎｍ）のパルスレーザ光を出射させるべく、所定のタイミングでレーザ光源ユニット１３にフラッシュランプスタンバイ信号を出力する（ステップＡ２）。レーザ光源ユニット１３の発光制御部６１は、フラッシュランプスタンバイ信号が入力されると、フラッシュランプ５２を点灯させ、レーザロッド５１の励起を開始させる（ステップＡ３）。発光制御部６１は、例えば複屈折フィルタ５６の回転変位位置が波長７５０ｎｍに対応した位置になるタイミングから逆算されたタイミングで、フラッシュランプ５２を点灯させる。
【００５４】
発光制御部６１は、フラッシュランプ５２の点灯後、複屈折フィルタ５６の回転変位位置が波長系列を構成する１つ目の波長（７５０ｎｍ）に対応した位置になるタイミングでＱスイッチ５５をオンにする（ステップＡ４）。Ｑスイッチ５５がオンになることで、レーザ光源ユニット１３は、波長７５０ｎｍのパルスレーザ光を出射する。発光制御部６１は、Ｑスイッチ５５をオンにするタイミングに同期して、Ｑスイッチ同期信号を超音波ユニット１２に出力する（ステップＡ５）。
【００５５】
レーザ光源ユニット１３から出射した波長７５０ｎｍのパルスレーザ光は、例えばプローブ１１まで導光され、プローブ１１から被検体の１つ目の部分領域に照射される。被検体内では、光吸収体が照射されたパルスレーザ光のエネルギーを吸収することで、光音響信号が発生する。プローブ１１は、被検体内で発生した光音響信号を検出する。プローブ１１で検出された光音響信号は、受信回路２１にて受信される。
【００５６】
トリガ制御回路３０は、フラッシュランプトリガ信号の出力後、Ｑスイッチ同期信号を受け取るまで待機する。トリガ制御回路３０は、Ｑスイッチ同期信号を受け取ると、ＡＤ変換手段２２にサンプリングトリガ信号を出力する。ＡＤ変換手段２２は、受信回路２１で受信された光音響信号を、所定のサンプリング周期でサンプリングする（ステップＡ６）。ＡＤ変換手段２２でサンプリングされた光音響信号は、受信メモリ２３に第１の光音響データとして格納される。
【００５７】
制御手段３１は、残りの波長があるか、つまり波長系列を構成する全ての波長のパルスレーザ光を出射させたか否かを判断する（ステップＡ７）。残りの波長がある場合、次の波長のパルスレーザ光を出射させるために、ステップＡ２に戻り、トリガ制御回路３０からレーザ光源ユニット１３にフラッシュランプスタンバイ信号を出力する。発光制御部６１は、ステップＡ３でフラッシュランプ５２を点灯し、ステップＡ４で、複屈折フィルタ５６が波長系列を構成する２つ目の波長（８００ｎｍ）に対応した回転変位位置になるタイミングでＱスイッチ５５をオンにし、パルスレーザ光を出射させる。
【００５８】
レーザ光源ユニット１３から出射した波長８００ｎｍのパルスレーザ光は、例えばプローブ１１まで導光され、プローブ１１から被検体の１つ目の部分領域に照射される。プローブ１１は、被検体内の光吸収体が波長８００ｎｍのパルスレーザ光を吸収することで発生した光音響信号を検出する。トリガ制御回路３０は、Ｑスイッチ同期信号を受け取ると、ＡＤ変換手段２２にサンプリングトリガ信号を出力し、ＡＤ変換手段２２は、ステップＡ６で光音響信号のサンプリングを行う。ＡＤ変換手段２２でサンプリングされた光音響信号は、受信メモリ２３に第２の光音響データとして格納される。光音響画像生成装置１０は、波長系列を構成する各波長に対してステップＡ１からＡ６を実行し、波長系列を構成する各波長のパルスレーザ光を被検体に照射して、被検体から光音響信号を検出する。
【００５９】
制御手段３１は、ステップＡ７で残りの波長がないと判断すると、部分領域を全て選択したか否かを判断する（ステップＡ８）。選択すべき部分領域が残っておりときは、ステップＡ２に戻る。光音響画像生成装置１０は、各部分領域に対してステップＡ２からＡ７を実行し、各部分領域に波長系列を構成する各波長（７５０ｎｍ、８００ｎｍ）のパルスレーザ光を順次照射し、各部分領域に対応した第１の光音響データと第２の光音響データとを受信メモリ２３に格納する。全ての部分領域に対してパルスレーザ光の照射及び光音響信号の検出を行うと、１フレームの光音響画像を生成するために必要な光音響データが揃う。
【００６０】
制御手段３１は、ステップＡ８で全ての部分領域を選択したと判断すると、光音響画像の生成に処理を移す。複素数化手段２４は、受信メモリ２３から第１の光音響データと第２の光音響データとを読み出し、第１の光音響画像データを実部とし、第２の光音響画像データを虚部とした複素数データを生成する（ステップＡ９）。光音響画像再構成手段２５は、ステップＡ９で複素数化された複素数データから、フーリエ変換法（ＦＴＡ法）により画像再構成を行う（ステップＡ１０）。
【００６１】
位相情報抽出手段２６は、再構成された複素数データ（再構成画像）から位相情報を抽出する（ステップＡ１１）。位相情報抽出手段２６は、例えば再構成された複素数データがＸ＋ｉＹで表わされるとき、θ＝ｔａｎ−１（Ｙ／Ｘ）を位相情報として抽出する（ただし、Ｘ＝０の場合はθ＝９０°）。強度情報抽出手段２７は、再構成された複素数データから強度情報を抽出する（ステップＡ１２）。強度情報抽出手段２７は、例えば再構成された複素数データがＸ＋ｉＹで表わされるとき、（Ｘ２＋Ｙ２）１／２を強度情報として抽出する。
【００６２】
検波・対数変換手段２８は、ステップＡ１２で抽出された強度情報に対して検波・対数変換処理を施す。光音響画像構築手段２９は、ステップＡ１１で抽出された位相情報と、ステップＡ１２で抽出された強度情報に対して検波・対数変換処理を施したものとに基づいて、光音響画像を生成する（ステップＡ１３）。光音響画像構築手段２９は、例えば強度情報に基づいて光吸収体の分布画像における各画素の輝度（階調値）を決定し、位相情報に基づいて各画素の色を決定することで、光音響画像を生成する。生成された光音響画像は、画像表示手段１４に表示される。
【００６３】
ここで、複屈折フィルタ５６の回転速度は、出射すべくパルスレーザ光の波長系列に含まれる波長の数に応じて適宜決定すればよい。以下、波長系列が６つの波長（７２０ｎｍ、７４０ｎｍ、７６０ｎｍ、７８０ｎｍ、８００ｎｍ、８２０ｎｍ）を含む場合について説明する。図８は、波長系列が６つの波長を含む場合の各種トリガと発光タイミングとを示すタイミングチャートである。（ａ）は、時間変化に対する光共振器の発振波長特性（複屈折フィルタ５６の透過波長特性）を示す。（ｂ）はフラッシュランプトリガを示し、（ｃ）はＱスイッチトリガを示す。（ｄ）は、フラッシュランプの発光タイミング及びパルスレーザ光の出射タイミングを示す。
【００６４】
図４に示す、１回転でＦＳＲを１２回繰り返す複屈折フィルタ５６を４秒で１回転させると、複屈折フィルタ５６は１秒間に１／４回転し、１秒間にＦＳＲを３回繰り返す（ａ）。トリガ制御回路３０は、まず、レーザ光源ユニット１３から波長７２０ｎｍのパルスレーザ光を出射させるために、フラッシュランプスタンバイ信号を出力する。発光制御部６１は、フラッシュランプスタンバイ信号を受け取ると、時刻ｔ２１で、フラッシュランプ５２にフラッシュランプトリガ信号を出力し（ｂ）、フラッシュランプ５２を点灯させる（ｄ）。その後、発光制御部６１は、複屈折フィルタ５６の回転変位位置が波長７２０ｎｍに対応した位置となる時刻ｔ２２で、Ｑスイッチトリガ信号を出力し（ｃ）、Ｑスイッチ５５をオンにすることで光共振器から波長７２０ｎｍのパルスレーザ光を出射させる。
【００６５】
次いで、トリガ制御回路３０は、レーザ光源ユニット１３から波長７４０ｎｍのパルスレーザ光を出射させるために、フラッシュランプスタンバイ信号を出力する。発光制御部６１は、フラッシュランプスタンバイ信号を受け取ると、時刻ｔ２３で、フラッシュランプ５２にフラッシュランプトリガ信号を出力し（ｂ）、フラッシュランプ５２を点灯させる（ｄ）。その後、発光制御部６１は、複屈折フィルタ５６の回転変位位置が波長７４０ｎｍに対応した位置となる時刻ｔ２４でＱスイッチトリガ信号を出力し（ｃ）、Ｑスイッチ５５をオンにすることで光共振器から７４０ｎｍのパルスレーザ光を出射させる。
【００６６】
以降同様に、発光制御部６１は、時刻ｔ２５、ｔ２７、ｔ２９、ｔ３１で、フラッシュランプ５２にフラッシュランプトリガ信号を出力する。また、時刻ｔ２６、ｔ２８、ｔ３０、ｔ３２で、Ｑスイッチ５５にＱスイッチトリガ信号を出力し、各時刻における複屈折フィルタ５６の透過波長に応じた波長のパルスレーザ光を出射させる。時刻ｔ２６、ｔ２８、ｔ３０、及びｔ３２おける複屈折フィルタ５６の透過波長はそれぞれ７６０ｎｍ、７８０ｎｍ、８００ｎｍ、及び８２０ｎｍであり、レーザ光源ユニット１３は、１／３秒間に、７２０ｎｍから８２０ｎｍまでの範囲で波長が２０ｎｍずつ長くなる６つのパルスレーザ光を出射する。
【００６７】
図８の例では、レーザ光源ユニット１３は、１／３秒の間に、波長７２０ｎｍから８２０ｎｍまでの６つの波長のパルスレーザ光を出射する。レーザ光源ユニット１３は、６つの波長を切り替えつつパルスレーザ光を１秒間に１８回出射する（１８Ｈｚ動作）。言い換えると、１セット６波長のパルスレーザ光を、１秒間に３セット出射する
【００６８】
複屈折フィルタ５６の回転速度は、１回のＦＳＲで、波長系列を構成する各波長のパルスレーザ光が出射できるように設定されていることが好ましい。例えば、複屈折フィルタ５６が１回転の間に繰り返すＦＳＲの回数をｋ［回／回転］とし、波長系列に含まれる波長の数をｎ［個］とし、単位時間当たりのパルスレーザ光の発行回数をｍ［回／秒］としたとき、複屈折フィルタ５６の回転速度を、ｖ＝ｍ／（ｋ×ｎ）［回転／秒］で決定した値とすることができる。この場合、１ＦＳＲあたりｎ個の波長を切り替えつつ、１秒間にｍ個のパルスレーザを出射できる（ｍＨｚ動作）。
【００６９】
本実施形態では、複屈折フィルタ５６を所定の回転速度で回転させつつ、フラッシュランプ５２を点灯してレーザロッド５１を励起し、その励起後、複屈折フィルタ５６の回転変位位置が、出射すべきパルスレーザ光の波長に対応した位置となるタイミングでＱスイッチ５５をオンにする。複屈折フィルタ５６の回転速度を遅くするほど、例えば複屈折フィルタ５６の１回分のＦＳＲの中で選択可能な発振波長の数を増やすことができる。一方、波長系列に含まれる波長が２つであるような場合は、複屈折フィルタ５６の回転速度を速めることで、２つの波長の切り替えを高速化することができる。このように、本実施形態では、複屈折フィルタ５６の回転速度を制御することで、レーザ光源ユニット１３から所望の波長系列でパルスレーザ光を出射することができる。
【００７０】
本実施形態では、２つの波長で得られた第１の光音響データと、第２の光音響データとの何れか一方を実部、他方を虚部とした複素数データを生成し、その複素数データからフーリエ変換法により再構成画像を生成している。このようにする場合、第１の光音響データと第２の光音響データとを別々に再構成する場合に比べて、再構成を効率的に行うことができる。複数の波長のパルスレーザ光を照射し、各波長のパルスレーザ光を照射したときの光音響信号（光音響データ）を用いることで、各光吸収体の光吸収特性が波長に応じて異なることを利用した機能イメージングを行うことができる。
【００７１】
また、本実施形態では、例えば光照射領域が３つの部分領域に分かれているときには、第１の部分領域に対して第１の波長のパルスレーザ光、第２の波長のパルスレーザ光を順次に照射し、次いで、第２の部分領域に対して第１の波長のパルスレーザ光、第２の波長のパルスレーザ光を順次に照射し、その後、第３の部分領域に対して第１の波長のパルスレーザ光、第２の波長のパルスレーザ光を順次に照射する。本実施形態では、ある部分領域に対して第１の波長のパルスレーザ光、第２の波長のパルスレーザ光を連続的に照射した後に、次の部分領域に移っている。この場合、第１の波長のパルスレーザ光を３つの部分領域に照射した後に、第２の波長のパルスレーザ光を３つの部分領域に照射する場合に比べて、同じ位置において第１の波長のパルスレーザ光を照射してから第２のパルスレーザ光が照射されるまでの間の時間を短くすることができる。第１の波長のパルスレーザ光が照射されてから第２の波長のパルスレーザ光が照射されるまでの間の時間を短縮することで、第１の光音響データと第２の光音響データとの不整合を抑制することができる。
【００７２】
引き続き、本発明の第２実施形態を説明する。図９は、本発明の第２実施形態の光音響画像生成装置を示す。本実施形態の光音響画像生成装置１０ａは、超音波ユニット１２ａが、図１に示す第１実施形態の光音響画像生成装置１０における超音波ユニット１２の構成に加えて、データ分離手段３２、超音波画像再構成手段３３、検波・対数変換手段３４、超音波画像構築手段３５、画像合成手段３６、及び送信制御回路３７を有する。本実施形態の光音響画像生成装置１０は、光音響画像に加えて、超音波画像を生成する点で第１実施形態と相違する。その他の部分は、第１実施形態と同様でよい。
【００７３】
本実施形態では、プローブ１１は、光音響信号の検出に加えて、被検体に対する超音波の出力（送信）、及び送信した超音波に対する被検体からの反射超音波の検出（受信）を行う。トリガ制御回路３０は、超音波画像の生成時は、送信制御回路３７に超音波送信を指示する旨の超音波送信トリガ信号を送る。送信制御回路３７は、トリガ信号を受けると、プローブ１１から超音波を送信させる。プローブ１１は、超音波の送信後、被検体からの反射超音波を検出する。
【００７４】
プローブ１１が検出した反射超音波は、受信回路２１を介してＡＤ変換手段２２に入力される。トリガ制御回路３０は、超音波送信のタイミングに合わせてＡＤ変換手段２２にサンプリグトリガ信号を送り、反射超音波のサンプリングを開始させる。ＡＤ変換手段２２は、反射超音波のサンプリングデータ（反射超音波データ）を受信メモリ２３に格納する。
【００７５】
データ分離手段３２は、受信メモリ２３に格納された反射超音波データと、第１及び第２の光音響データとを分離する。データ分離手段３２は、反射超音波データを超音波画像再構成手段３３に渡し、第１及び第２の光音響データを複素数化手段２４に渡す。第１及び第２の光音響データに基づく光音響画像の生成は第１実施形態と同様である。データ分離手段３２は、分離した反射超音波のサンプリングデータを、超音波画像再構成手段３３に入力する。
【００７６】
超音波画像再構成手段３３は、プローブ１１の複数の超音波振動子で検出された反射超音波（そのサンプリングデータ）に基づいて、超音波画像の各ラインのデータを生成する。超音波画像再構成手段３３は、例えばプローブ１１の６４個の超音波振動子からのデータを、超音波振動子の位置に応じた遅延時間で加算し、１ライン分のデータを生成する（遅延加算法）。
【００７７】
検波・対数変換手段３４は、超音波画像再構成手段３３が出力する各ラインのデータの包絡線を求め、求めた包絡線を対数変換する。超音波画像構築手段３５は、対数変換が施された各ラインのデータに基づいて、超音波画像を生成する。超音波画像再構成手段３３、検波・対数変換手段３４、及び超音波画像構築手段３５は、反射超音波に基づいて超音波画像を生成する超音波画像生成手段を構成する。
【００７８】
画像合成手段３６は、光音響画像と超音波画像とを合成する。画像合成手段３６は、例えば光音響画像と超音波画像とを重畳することで画像合成を行う。その際、画像合成手段３６は、光音響画像と超音波画像とで、対応点が同一の位置となるように位置合わせをすることが好ましい。合成された画像は、画像表示手段１４に表示される。画像合成を行わずに、画像表示手段１４に、光音響画像と超音波画像とを並べて表示し、或いは光音響画像と超音波画像とを切り替えてすることも可能である。
【００７９】
図１０は、光音響画像生成装置１０ａの動作手順を示している。以下では、被検体のパルスレーザ光が照射される領域が複数の部分領域に分割されているものとして説明する。超音波ユニット１２のトリガ制御回路３０は、所定の回転速度でレーザ光源ユニット１３内の複屈折フィルタ５６を回転させる旨のＢＲＦ制御信号をレーザ光源ユニット１３に出力する（ステップＢ１）。
【００８０】
トリガ制御回路３０は、光音響信号の受信準備が整うと、波長系列を構成する１つ目の波長（７５０ｎｍ）のパルスレーザ光を出射させるべく、フラッシュランプスタンバイ信号を出力する（ステップＢ２）。レーザ光源ユニット１３の発光制御部６１は、フラッシュランプスタンバイ信号が入力されると、フラッシュランプ５２を点灯させ、レーザロッド５１の励起を開始させる（ステップＢ３）。
【００８１】
発光制御部６１は、フラッシュランプ５２の点灯後、複屈折フィルタ５６の回転変位位置が波長系列を構成する１つ目の波長（７５０ｎｍ）に対応した位置になるタイミングでＱスイッチ５５をオンにする（ステップＢ４）。Ｑスイッチ５５がオンになることで、レーザ光源ユニット１３は、波長７５０ｎｍのパルスレーザ光を出射する。発光制御部６１は、Ｑスイッチ５５をオンにするタイミングに同期して、Ｑスイッチ同期信号を超音波ユニット１２に出力する（ステップＢ５）。
【００８２】
レーザ光源ユニット１３から出射した波長７５０ｎｍのパルスレーザ光は、例えばプローブ１１まで導光され、プローブ１１から被検体の１つ目の領域に照射される。被検体内では、光吸収体が照射されたパルスレーザ光のエネルギーを吸収することで、光音響信号が発生する。プローブ１１は、被検体内で発生した光音響信号を検出する。トリガ制御回路３０は、Ｑスイッチ同期信号を受け取ると、ＡＤ変換手段２２にサンプリングトリガ信号を出力する。ＡＤ変換手段２２は、受信回路２１を介してプローブ１１で検出された光音響信号を受信し、所定のサンプリング周期で光音響信号をサンプリングする（ステップＢ６）。ＡＤ変換手段２２でサンプリングされた光音響信号は、受信メモリ２３に第１の光音響データとして格納される。
【００８３】
制御手段３１は、残りの波長があるか、つまり波長系列を構成する全ての波長のパルスレーザ光を出射させたか否かを判断する（ステップＢ７）。残りの波長がある場合、次の波長のパルスレーザ光を出射させるために、ステップＢ２に戻り、トリガ制御回路３０からレーザ光源ユニット１３にフラッシュランプスタンバイ信号を出力する。発光制御部６１は、ステップＢ３でフラッシュランプ５２を点灯し、ステップＢ４で、複屈折フィルタ５６が波長系列を構成する２つ目の波長（８００ｎｍ）に対応した回転変位位置になるタイミングでＱスイッチ５５をオンにし、パルスレーザ光を出射させる。
【００８４】
レーザ光源ユニット１３から出射した波長８００ｎｍのパルスレーザ光は、例えばプローブ１１まで導光され、プローブ１１から被検体の１つ目の部分領域に照射される。プローブ１１は、被検体内の光吸収体が波長８００ｎｍのパルスレーザ光を吸収することで発生した光音響信号を検出する。トリガ制御回路３０は、Ｑスイッチ同期信号を受け取ると、ＡＤ変換手段２２にサンプリングトリガ信号を出力し、ＡＤ変換手段２２は、ステップＢ６で光音響信号のサンプリングを行う。ＡＤ変換手段２２でサンプリングされた光音響信号は、受信メモリ２３に第２の光音響データとして格納される。光音響画像生成装置１０は、波長系列を構成する各波長に対してステップＢ１からＢ６を実行し、波長系列を構成する各波長のパルスレーザ光を被検体に照射して、被検体から光音響信号を検出する。ステップＢ１からＢ６は、図７のステップＡ１からＡ６と同様でよい。
【００８５】
制御手段３１は、ステップＢ７で残りの波長がないと判断すると、超音波の送受信に処理を移す。トリガ制御回路３０は、送信制御回路３７を介してプローブ１１から被検体に超音波を送信させる（ステップＢ８）。ステップＢ８では、被検体のパルスレーザ光が照射された部分領域と同じ領域に対して超音波を送信する。プローブ１１は、送信した超音波に対する反射超音波を検出する（ステップＢ９）。検出された反射超音波は、受信回路２１を経てＡＤ変換手段２２でサンプリングされ、受信メモリ２３に反射超音波データとして格納される。
【００８６】
制御手段３１は、部分領域を全て選択したか否かを判断する（ステップＢ１０）。選択すべき部分領域が残っておりときは、ステップＢ２に戻る。光音響画像生成装置１０は、各部分領域に対してステップＢ２からＢ７を実行し、各部分領域に波長系列を構成する各波長（７５０ｎｍ、８００ｎｍ）のパルスレーザ光を順次照射して、第１の光音響データと第２の光音響データとを受信メモリ２３に格納する。また、ステップＢ８及びＢ９を実行して、反射超音波データを受信メモリ２３に格納する。全ての部分領域に対してパルスレーザ光の照射と光音響信号の検出、及び超音波の送受信を行うと、１フレームの光音響画像及び超音波画像を生成するために必要なデータが揃う。
【００８７】
制御手段３１は、ステップＢ１０で全ての部分領域を選択したと判断すると、光音響画像及び超音波画像の生成に処理を移す。データ分離手段３２は、第１及び第２の光音響データと反射超音波データとを分離する。データ分離手段３２は、分離した第１及び第２の光音響データを複素数化手段２４に渡し、反射超音波データを超音波画像再構成手段３３に渡す。複素数化手段２４は、第１の光音響画像データを実部とし、第２の光音響画像データを虚部とした複素数データを生成する（ステップＢ１１）。光音響画像再構成手段２５は、ステップＢ１１で複素数化された複素数データから、フーリエ変換法（ＦＴＡ法）により画像再構成を行う（ステップＢ１２）。
【００８８】
位相情報抽出手段２６は、再構成された複素数データから位相情報を抽出する（ステップＢ１３）。強度情報抽出手段２７は、再構成された複素数データから強度情報を抽出する（ステップＢ１４）。検波・対数変換手段２８は、ステップＢ１４で抽出された強度情報に対して検波・対数変換処理を施す。光音響画像構築手段２９は、ステップＢ１３で抽出された位相情報と、ステップＢ１４で抽出された強度情報に対して検波・対数変換処理を施したものとに基づいて、光音響画像を生成する（ステップＢ１５）。ステップＢ１１からＢ１５は、図７のステップＡ９からＡ１３と同様でよい。
【００８９】
超音波画像再構成手段３３は、例えば遅延加算法により、超音波画像の各ラインのデータを生成する。検波・対数変換手段３４は、超音波画像再構成手段３３が出力する各ラインのデータの包絡線を求め、求めた包絡線を対数変換する。超音波画像構築手段３５は、対数変換が施された各ラインのデータに基づいて、超音波画像を生成する（ステップＢ１６）。画像合成手段３６は、光音響画像と超音波画像とを合成し、合成した画像を画像表示手段１４に表示する（ステップＢ１７）。
【００９０】
本実施形態では、光音響画像生成装置は、光音響画像に加えて超音波画像を生成する。超音波画像を参照することで、光音響画像では画像化することができない部分を観察することができる。その他の効果は、第１実施形態と同様である。
【００９１】
なお、上記各実施形態では、第１の光音響データと第２の光音響データとを複素数化する例について説明したが、複素数化せずに、第１の光音響データと第２の光音響データとを別々に再構成してもよい。さらに、ここでは、複素数化して位相情報を用いて第１の光音響データと第２の光音響データとの比を計算しているが、両者の強度情報から比を計算しても同様の効果が得られる。また、強度情報は、第１の再構成画像における信号強度と、第２の再構成画像における信号強度とに基づいて生成できる。
【００９２】
光音響画像の生成に際して、被検体に照射されるパルスレーザ光の波長の数は２つには限られず、３以上のパルスレーザ光を被検体に照射し、各波長に対応する光音響データに基づいて光音響画像を生成してもよい。その場合、例えば位相情報抽出手段２６は、各波長に対応する光音響データ間での相対的な信号強度の大小関係を位相情報として生成すればよい。また、強度情報抽出手段２７は、例えば各波長に対応する光音響データにおける信号強度を１つにまとめたものを強度情報として生成すればよい。
【００９３】
図３では、駆動手段５７としてサーボモータを用い、サーボモータを回転させることで複屈折フィルタ５６を回転させる例について説明したが、駆動手段５７は複屈折フィルタ５６を回転させるものであればよく、サーボモータには限定されない。例えば駆動手段５７に、ステッピングモータを用いることもできる。その場合、回転制御部６０は、複屈折フィルタ５６が所定の回転速度で回転するように、ステッピングモータに供給するパルス信号のパルス間隔を制御すればよい。ステッピングモータを用いる場合、回転変位検出手段５８は省略してもよい。
【００９４】
以上、本発明をその好適な実施形態に基づいて説明したが、本発明のレーザ光源ユニット及びその制御方法、並びに光音響画像生成装置及び方法は、上記実施形態にのみ限定されるものではなく、上記実施形態の構成から種々の修正及び変更を施したものも、本発明の範囲に含まれる。
【符号の説明】
【００９５】
１０：光音響画像生成装置
１１：プローブ
１２：超音波ユニット
１３：レーザ光源ユニット
１４：画像表示手段
２１：受信回路
２２：ＡＤ変換手段
２３：受信メモリ
２４：複素数化手段
２５：光音響画像再構成手段
２６：位相情報抽出手段
２７：強度情報抽出手段
２８：検波・対数変換手段
２９：光音響画像構築手段
３０：トリガ制御回路
３１：制御手段
３２：データ分離手段
３３：超音波画像再構成手段
３４：検波・対数変換手段
３５：超音波画像構築手段
３６：画像合成手段
３７：送信制御回路
５１：レーザロッド
５２：フラッシュランプ
５３、５４：ミラー
５５：Ｑスイッチ
５６：複屈折フィルタ
５７：駆動手段
５８：回転変位検出手段
５９：制御部
６０：回転制御部
６１：発光制御部
【図１】
0005662973.tif 000002
【図２】
0005662973.tif 000003
【図３】
0005662973.tif 000004
【図４】
0005662973.tif 000005
【図５】
0005662973.tif 000006
【図６】
0005662973.tif 000007
【図７】
0005662973.tif 000008
【図８】
0005662973.tif 000009
【図９】
0005662973.tif 000010
【図１０】
0005662973.tif 000011
【図１１】
0005662973.tif 000012