(57)【要約】【課題】  本発明は、複数回の光照射により発生する光音響波を利用する光音響装置において、単位時間当たりに得られる光音響波の受信信号のＳＮ比を向上させることを目的とする。【解決手段】  本発明に係る光音響装置は、光照射手段、受信手段、及び処理手段を有し、光照射手段は、複数の符号系列に対応する複数の強度変調光を被検体に照射し、複数の強度変調光のそれぞれは、符号系列の基準タイミングにおいて互いに異なる光量変化量を有し、受信手段は、複数の強度変調光が照射されることにより発生する光音響波を受信することにより受信信号を出力し、処理手段は、複数の符号系列に関する情報に基づいて、受信信号に対して復号処理を行うことにより、複数の強度変調光に対応する光量変化量の少なくとも１つに対応する復号信号を生成する。【選択図】  図１０
2019068914.tif 000002
【特許請求の範囲】
【請求項１】
光照射手段、受信手段、及び処理手段を有し、
前記光照射手段は、第１の符号系列に対応する第１の強度変調光と、前記第１の符号系列とは異なる第２の符号系列に対応する第２の強度変調光とを被検体に照射し、
前記第１の強度変調光は、前記第１の符号系列の基準タイミングにおいて第１の光量変化量を有し、
前記第２の強度変調光は、前記第２の符号系列の基準タイミングにおいて前記第１の光量変化量とは異なる第２の光量変化量を有し、
前記受信手段は、前記第１及び第２の強度変調光が前記被検体に照射されることにより発生する光音響波を受信することにより受信信号を出力し、
前記処理手段は、前記第１及び第２の符号系列に関する情報に基づいて、前記受信信号に対して復号処理を行うことにより、前記第１の光量変化量に対応する第１の復号信号及び前記第２の光量変化量に対応する第２の復号信号の少なくとも一方を生成する
ことを特徴とする光音響装置。
【請求項２】
前記第１の強度変調光の照射期間と前記第２の強度変調光の照射期間とは、少なくとも一部が重なる
ことを特徴とする請求項１に記載の光音響装置。
【請求項３】
前記処理手段は、
前記第１の復号信号を用いて第１の光音響画像を生成し、
前記第２の復号信号を用いて第２の光音響画像を生成する
ことを特徴とする請求項１または２に記載の光音響装置。
【請求項４】
前記第１の光音響画像及び前記第２の光音響画像を表示手段に並べて表示させる、または、切り替えて表示させる
ことを特徴とする請求項３に記載の光音響装置。
【請求項５】
前記受信手段は、第１の受信帯域を有する第１のトランスデューサと、前記第１の受信帯域とは異なる第２の受信帯域を有する第２のトランスデューサとを有し、
前記第１のトランスデューサは、前記光音響波を受信することにより第１の受信信号を出力し、
前記第２のトランスデューサは、前記光音響波を受信することにより第２の受信信号を出力し、
前記処理手段は、
前記第１のトランスデューサから出力された前記第１の受信信号に対して復号処理を行うことにより生成された前記第１の復号信号を用いて、前記第１の光音響画像を生成し、
前記第２のトランスデューサから出力された前記第２の受信信号に対して復号処理を行うことにより生成された前記第２の復号信号を用いて、前記第２の光音響画像を生成する
ことを特徴とする請求項３または４に記載の光音響装置。
【請求項６】
前記第１の受信帯域の中心周波数は、前記第２の受信帯域の中心周波数よりも高く、
前記第１の光量変化量は、前記第２の光量変化量よりも大きい
ことを特徴とする請求項５に記載の光音響装置。
【請求項７】
前記処理手段は、
前記第１の復号信号と前記第２の復号信号との両方を用いて、第３の光音響画像を生成する
ことを特徴とする請求項１から６のいずれか１項に記載の光音響装置。
【請求項８】
前記処理手段は、前記第１の復号信号を第１の重みで重みづけし、前記第２の復号信号を前記第１の重みとは異なる第２の重みで重みづけすることにより、前記第３の光音響画像を生成する
ことを特徴とする請求項７に記載の光音響装置。
【請求項９】
前記第１または第２の光量変化量は、前記基準タイミングにおける前記第１または第２の強度変調光の光量の時間微分値の絶対値であることを特徴とする請求項１から８のいずれか１項に記載の光音響装置。
【請求項１０】
前記第１または第２の光量変化量は、前記第１または第２の強度変調光の光量の時間変化を示すグラフにおける前記基準タイミングでの傾きの絶対値である
ことを特徴とする請求項１から８のいずれか１項に記載の光音響装置。
【請求項１１】
前記第１または第２の光量変化量は、前記基準タイミングを含む所定の期間中の前記第１または第２の強度変調光の光量変化量である
ことを特徴とする請求項１から８のいずれか１項に記載の光音響装置。
【請求項１２】
前記第１または第２の光量変化量は、前記所定の期間中の前記第１または第２の強度変調光の光量の最大値と最小値との差である
ことを特徴とする請求項１１に記載の光音響装置。
【請求項１３】
前記所定の期間は、前記強度変調光において、前記基準タイミングの前に光量が極値となったタイミングと、前記基準タイミングの後に光量が極値となったタイミングとの間の期間である
ことを特徴とする請求項１１または１２に記載の光音響装置。
【請求項１４】
前記第１または第２の光量変化量は、前記基準タイミングにおける前記第１または第２の強度変調光の立ち上がりまたは立ち下がりの時間である
ことを特徴とする請求項１から８のいずれか１項に記載の光音響装置。
【請求項１５】
前記第１及び第２の符号系列のそれぞれは、正の符号要素及び負の符号要素を含む
ことを特徴とする請求項１から１４のいずれか１項に記載の光音響装置。
【請求項１６】
前記光照射手段は、
半導体レーザまたは発光ダイオードで構成された光源と、
前記光源へ駆動電流を投入する駆動手段と、
を有し、
前記駆動手段は、前記正の符号要素に対応する基準タイミングに合わせて、前記正の符号要素に対応する正の強度変調光を発生させるための第１の駆動電流を前記光源に投入し、
前記駆動手段は、前記負の符号要素に対応する基準タイミングに合わせて、前記負の符号要素に対応する負の強度変調光を発生させるための第２の駆動電流を前記光源に投入する
ことを特徴とする請求項１５に記載の光音響装置。
【請求項１７】
前記光照射手段は、前記第１の強度変調光を発する半導体レーザまたは発光ダイオードで構成された第１の光源と、前記第２の強度変調光を発する半導体レーザまたは発光ダイオードで構成された第２の光源とを含む
ことを特徴とする請求項１から１６のいずれか１項に記載の光音響装置。
【請求項１８】
光照射手段、受信手段、及び処理手段を有し、
前記光照射手段は、複数の符号系列に対応する複数の強度変調光を被検体に照射し、
前記複数の強度変調光のそれぞれは、符号系列の基準タイミングにおいて互いに異なる光量変化量を有し、
前記受信手段は、前記複数の強度変調光が照射されることにより発生する光音響波を受信することにより受信信号を出力し、
前記処理手段は、前記複数の符号系列に関する情報に基づいて、前記受信信号に対して復号処理を行うことにより、前記複数の強度変調光に対応する前記光量変化量の少なくとも１つに対応する復号信号を生成する
ことを特徴とする光音響装置。
【請求項１９】
複数の符号系列に対応する複数の強度変調光が被検体に照射されることにより発生する光音響波の受信信号に対して、復号処理を行うことにより復号信号を生成する情報取得方法であって、
前記複数の強度変調光のそれぞれは、符号系列の基準タイミングにおいて互いに異なる光量変化量を有し、
前記複数の符号系列に関する情報に基づいて、前記受信信号に対して復号処理を行うことにより、複数の強度変調光に対応する前記光量変化量の少なくとも１つに対応する前記復号信号を生成することを特徴とする情報取得方法。
【請求項２０】
請求項１９に記載の情報取得方法をコンピュータに実行させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、光音響効果により発生する光音響波を利用する光音響装置に関する。
【背景技術】
【０００２】
被検体に光を照射し、照射光に基づいて得られる被検体内の情報を画像化する光イメージング装置の研究が医療分野で積極的に進められている。この光イメージング技術の一つとして、Ｐｈｏｔｏ  Ａｃｏｕｓｔｉｃ  Ｔｏｍｏｇｒａｐｈｙ（ＰＡＴ）がある。ＰＡＴでは、光源から発生した光を被検体に照射し、被検体内で伝搬、拡散した光のエネルギーを吸収した組織から発生した音響波を検出する。この音響波発生の現象を光音響効果と呼び、発生した音響波を光音響波と呼ぶ。音響波は一般には超音波である。
【０００３】
特許文献１は、パルス光が複数並んだパルス列を被検体に照射することを開示する。また、特許文献１は、光照射により被検体内で発生した光音響波に基づいて、画像データを再構成することを開示する。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１４−３９８０１号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
複数回の光照射を利用した光音響トモグラフィの場合、典型的には光照射により発生した光音響波の受信が完了した後に、次の光照射を行い、次に発生する光音響波の受信を開始する。しかしながら、この方法では、単位時間当たりの光照射回数が制限されてしまうため、単位時間当たりに得られる光音響波の受信信号のＳＮ比の向上が困難である。
【０００６】
そこで、本発明は、複数回の光照射により発生する光音響波を利用する光音響装置において、単位時間当たりに得られる光音響波の受信信号のＳＮ比を向上させることを目的とする。
【課題を解決するための手段】
【０００７】
本発明に係る光音響装置は、本発明に係る光音響装置は、光照射手段、受信手段、及び処理手段を有し、光照射手段は、複数の符号系列に対応する複数の強度変調光を被検体に照射し、複数の強度変調光のそれぞれは、符号系列の基準タイミングにおいて互いに異なる光量変化量を有し、受信手段は、複数の強度変調光が照射されることにより発生する光音響波を受信することにより受信信号を出力し、処理手段は、複数の符号系列に関する情報に基づいて、受信信号に対して復号処理を行うことにより、複数の強度変調光に対応する光量変化量の少なくとも１つに対応する復号信号を生成する。
【発明の効果】
【０００８】
本発明に係る光音響装置によれば、単位時間当たりに得られる光音響波の受信信号のＳＮ比を向上させることができる。
【図面の簡単な説明】
【０００９】
【図１】正負の符号要素に対応する照射光の強度、及び、光音響波の受信信号を模式的に示す図
【図２】符号系列に対応する照射光の強度、及び、光音響波の受信信号を模式的に示す図
【図３】本実施形態に係る光音響装置を示すブロック図
【図４】本実施形態に係るコンピュータとその周辺の構成を示すブロック図
【図５】半導体レーザの特性を示す図
【図６】正の符号要素に対応する受信信号を示す図
【図７Ａ】立ち上がり時間を５０ｎｓから２５０ｎｓまで変化させたときの駆動電流の時間波形と受信信号の時間波形とを示す図
【図７Ｂ】立ち上がり時間を３００ｎｓから５００ｎｓまで変化させたときの駆動電流の時間波形と受信信号の時間波形とを示す図
【図８Ａ】立ち下がり時間を５０ｎｓから２５０ｎｓまで変化させたときの駆動電流の時間波形と受信信号の時間波形とを示す図
【図８Ｂ】立ち下がり時間を３００ｎｓから５００ｎｓまで変化させたときの駆動電流の時間波形と受信信号の時間波形とを示す図
【図９】照射光の立ち上がり時間と受信信号との関係を示す図
【図１０】実施例１に係る符号化のシーケンスを示す図
【図１１】実施例１に係る光の照射と光吸収体の配置を示す図
【図１２】実施例１に係る駆動電流、受信信号を示す図
【図１３】実施例１に係る受信信号を示す図
【図１４】実施例１に係る別の駆動電流、受信信号を示す図
【図１５】実施例１に係る別の受信信号を示す図
【図１６】実施例１に係る復号信号を示す図
【図１７】実施例１に係る別の復号信号を示す図
【図１８】実施例１に係る駆動部の構成を示す図
【図１９】実施例２に係る光音響装置を示すブロック図
【図２０】実施例２に係る符号化のシーケンスを示す図
【発明を実施するための形態】
【００１０】
物質に光を照射すると光音響効果により音響波（光音響波とも呼ぶ）が発生することが知られている。
【００１１】
一般に、小さいあるいは細い光吸収体からは高周波成分を多く含む光音響波が発生する傾向がある。一方、大きいあるいは太い光吸収体からは低周波成分を多く含む光音響波が発生する傾向がある。
【００１２】
また、照射光の単位時間当たりの光量変化量も発生する光音響波の周波数帯域に影響を与える。照射光の単位時間当たりの光量変化量が大きいと、高周波成分を多く含む光音響波が発生する傾向がある。すなわち、光量が急峻に変化する場合、高周波成分を多く含む光音響波が発生する傾向がある。一方、照射光の単位時間当たりの光量変化量が小さいと、低周波成分を多く含む光音響波が発生する傾向がある。照射光の光量は［Ｊ］や［Ｊ／ｍ３］の単位で表されるパラメータであり、照射光の強度や出力とも呼ぶ。
【００１３】
したがって、小さいあるいは細い光吸収体を観察したい場合、単位時間当たりの光量変化量の大きい光を光吸収体に照射することが好ましい。逆に、大きいあるいは太い光吸収体を観察したい場合、単位時間当たりの光量変化量の小さい光を光吸収体に照射することが好ましい。
【００１４】
このように、観察したい光吸収体の大きさに合わせて適切な光量変化量を設定することが好ましい。しかし、物質中に含まれる光吸収体の大きさが未知である場合、適切な光量変化量を設定することが困難となる。その場合、互いに光量変化量の異なる複数の光を照射する方法が考えられる。
【００１５】
ところが、光量変化量の異なる複数の光を光吸収体に照射する場合、ある光量変化量を有する光の照射期間と別の光量変化量を有する光の照射期間とが重なった場合、それぞれの光照射によって発生した光音響波の受信信号を分離することができない場合がある。そこで、ある光量変化量を有する光の照射による光音響波を受信完了した後に、別の光量変化量を有する光の照射による光音響波を受信する、というように時間的に受信信号を分離する手法が知られている。しかしながら、このように光照射ごとに光音響波の受信に十分な時間を確保する手法では、単位時間当たりに得られる光音響波の受信信号のＳＮ比が低下してしまう。
【００１６】
そこで、本発明者は、複数の符号系列に対応する複数の強度変調光を被検体に照射して符号化を行うことにより、複数の光量変化量のそれぞれに対応する信号を分離することを見出した。すなわち、本発明者は、複数の強度変調光のそれぞれが、符号系列の基準タイミングにおいて互いに異なる光量変化量を有するように光照射を行うことを見出した。ここで、複数の強度変調光のそれぞれは、符号化された光として機能する。このように符号化された光照射によって発生する光音響波の受信信号を、符号化に利用した符号系列に関する情報を用いて復号することにより、複数の光量変化量のそれぞれに対応する復号信号を得ることができる。
【００１７】
なお、複数の強度変調光の照射期間は、少なくとも一部が重なっていてもよい。その場合も、後述する復号処理によって、複数の光量変化量のそれぞれに対応する復号信号に分離することができる。また、このような符号化・復号化により、単位時間当たりに得られる光音響波の受信信号のＳＮ比を向上させることができる。
【００１８】
なお、「基準タイミングにおける光量変化量」は、符号化の基準タイミングにおける強度変調光の光量の時間微分値の絶対値であってもよい。また、「基準タイミングにおける光量変化量」は、強度変調光の光量の時間変化を示すグラフにおける基準タイミングでの傾きの絶対値であってもよい。また、「基準タイミングにおける光量変化量」は、基準タイミングを含む所定の期間中の強度変調光の光量変化量（例えば、最大値と最小値との差）であってもよい。例えば、所定の期間とは、強度変調光において、基準タイミングの前に光量が極値となったタイミングと、基準タイミングの後に光量が極値となったタイミングとの間の期間であってもよい。また、ノイズ等による光量のばらつきがある場合、ノイズ等による光量のばらつきを無視したときに極値となるタイミング間の期間を所定の期間としてもよい。また、「基準タイミングにおける光量変化量」は、符号系列の基準タイミングにおける照射光の立ち上がり時間または立ち下がり時間であってもよい。典型的に、ピーク光量が同じである場合、立ち上がり時間または立ち下がり時間が長くなるほど単位時間当たりの光量変化量が小さくなり、立ち上がり時間または立ち下がり時間が短くなるほど単位時間当たりの光量変化量は大きくなる傾向がある。すなわち、基準タイミングにおける立ち上がり時間または立ち下がり時間を変化させることは、基準タイミングにおける光量変化量を変化させることに相当する。
【００１９】
光音響波の受信信号を処理する光音響装置において、照射光の制御により正負それぞれの符号要素を含む符号系列に基づいた符号化を行う方法を説明する。図１は、照射光の強度、及び、その照射光により発生する光音響波の受信信号のレベルの時間変化を模式的に示した図である。典型的に、図１（ａ）に示すように、照射光の強度の時間変化を正とすると、正のレベルの受信信号を得ることができる。一方、図１（ｂ）に示すように照射光の強度の時間変化を負とすると、負のレベルの受信信号を得ることができる。さらに、単位時間当たりの照射光の強度変化が大きいほど、受信信号のレベルも大きくなる傾向がある。なお、図１においては、音源から受信手段までの光音響波の伝搬時間については無視している。
【００２０】
図１に示すように、照射光の強度の時間変化の正負を制御することにより、受信信号のレベルの正負を制御している。すなわち、照射光の強度の時間変化の正負を制御することにより、符号化における符号系列を構成する符号要素の正負を制御している。例えば、図１（ａ）に示す光照射のタイミングにおける符号要素を｛１｝とし、図１（ｂ）に示す光照射のタイミングにおける符号要素を｛−１｝とし、これらの光照射を組み合わせることにより、正負それぞれの符号要素を含む符号系列を定義することができる。本明細書において、正の符号要素に対応する光音響波を発生させるための光を「正の強度変調光」とする。また、負の符号要素に対応する光音響波を発生させるための光を「負の強度変調光」とする。
【００２１】
以下、いくつかのパターンの符号系列に対応する光照射のシーケンスの例を、図２を用いて説明する。図２中の点線は、各符号要素の基準タイミングを示す。
【００２２】
図２（ａ）は、符号系列｛１、１｝に対応する照射光の強度、及び、光音響波の受信信号のレベルの時間変化を模式的に示した図である。図２（ａ）に示す照射光のシーケンスは、強度が短時間に急上昇した後に、時間と共になだらかに下がっていく光（正の強度変調光）を２つ連続させたものである。また、強度が短時間に急上昇するタイミングが、正の符号要素に対応する基準タイミングに対応するように調整されている。例えば、強度が短時間に急上昇する期間の中心のタイミングを基準タイミングと一致させることができる。この場合、基準タイミングに、正の大きな受信信号が得られる。この正の大きな受信信号が、正の符号要素｛１｝に対応する信号となる。
【００２３】
図２（ｂ）は、符号系列｛−１、―１｝に対応する照射光の強度、及び、光音響波の受信信号のレベルの時間変化を模式的に示した図である。図２（ｂ）に示す照射光のシーケンスは、強度が時間と共になだらかに上がっていた後に、短時間に急降下する光（負の強度変調光）を２つ連続させたものである。また、強度が短時間に急降下するタイミングが、負の符号要素に対応する基準タイミングに合うように調整されている。具体的には、強度が短時間に急降下する期間の中心のタイミングを基準タイミングと一致させている。この場合、基準タイミングに、負の大きな受信信号が得られる。この負の大きな受信信号が、負の符号要素｛−１｝に対応する信号となる。
【００２４】
図２（ｃ）は、符号系列｛１、−１｝に対応する照射光の強度、及び、光音響波の受信信号のレベルの時間変化を模式的に示した図である。図２（ｃ）に示す照射光のシーケンスは、図２（ａ）に示す正の強度変調光を照射した後に、図２（ｂ）に示す負の強度変調光を照射するものである。正の強度変調光の強度が急上昇するタイミングが正の符号要素｛１｝の基準タイミングに対応し、負の強度変調光の強度が急降下するタイミングが負の符号要素｛−１｝の基準タイミングと対応するように照射タイミングは制御されている。この場合、各基準タイミングに、正の大きな受信信号、および、負の大きな受信信号が得られる。
【００２５】
なお、正の強度変調光の時間と共になだらかに下がっていく部分と、負の強度変調光の時間と共になだらかに上がっていく部分とが時間的に重なるので、重畳部分は結果的に矩形波となる。このように、符号系列の正負の符号要素が隣り合う場合に、これらの基準タイミング間の光強度をほぼ一定にすることにより、その期間に不要な光音響波が発生しない。これにより、精度良く光照射による符号化を実現することができる。図２（ｃ）では、符号系列｛１、−１｝の例を説明したが、符号系列｛−１、１｝の場合にも同様に基準タイミング間の光強度をほぼ一定としてもよい。なお、基準タイミング間の光強度の時間変化は、光音響波を受信するトランスデューサの受信帯域に基づいて、この受信帯域から外れる周波数を有する光音響波を発生させる程度の所定の範囲内であれば、ほぼ一定とみなしてもよい。図２においては、音源から受信手段までの光音響波の伝搬時間については無視している。
【００２６】
このようにして、正負それぞれの符号要素を含む符号系列に対応する光照射を行い、正負それぞれの符号要素を含む符号化を行うことにより、正負それぞれの符号要素を含む符号系列に基づいて復号化を行う際の復号精度を向上させることができる。特に、固体レーザ等の高出力光源と比較して出力される光強度が小さい半導体レーザやＬＥＤ等の場合、単位時間当たりの照射回数を増やして、受信信号のＳＮ比を向上させたい。このような場合に、正負それぞれの符号要素を含む符号系列に基づいて、先に発生した光音響波の受信が完了する前に、次の光照射を実行して符号化を行うことにより、ＳＮ比の高い復号信号を精度良く取得することができる。
【００２７】
本発明に係る、光音響効果により発生する音響波とは、典型的には超音波であり、音波、音響波と呼ばれるものを含む。本発明は、光音響効果により発生した光音響波に基づいて画像データを取得する光音響装置に適用することができる。
【００２８】
本発明に係る光音響装置により得られる光音響画像は、光照射により発生した光音響波に由来するあらゆる画像を含む概念である。光音響画像は、光音響波の発生音圧（初期音圧）、光吸収エネルギー密度、及び光吸収係数、被検体を構成する物質の濃度（酸素飽和度など）などの少なくとも１つの被検体情報の空間分布を表す画像データである。
【００２９】
以下、図３を用いて本実施形態に係る光音響装置の構成を説明する。図３は、光音響装置全体の概略ブロック図である。本実施形態に係る光音響装置は、光照射部１１０、受信部１２０、信号収集部１４０、コンピュータ１５０、表示部１６０、及び入力部１７０を有する。
【００３０】
光照射部１１０が光を被検体１００に照射し、被検体１００から音響波が発生する。光に起因して光音響効果により発生する音響波を光音響波とも呼ぶ。受信部１２０は、光音響波を受信することによりアナログ信号としての電気信号（光音響信号）を出力する。
【００３１】
信号収集部１４０は、受信部１２０から出力されたアナログ信号をデジタル信号に変換し、コンピュータ１５０に出力する。コンピュータ１５０は、信号収集部１４０から出力されたデジタル信号を、光音響波に由来する信号データとして記憶する。
【００３２】
処理手段としてのコンピュータ１５０は、記憶されたデジタル信号に対して後述する処理を行うことにより、光音響画像を表す画像データを生成する。また、コンピュータ１５０は、得られた画像データに対して表示のための画像処理を施した後に、画像データを表示部１６０に出力する。表示部１６０は、光音響画像を表示する。ユーザーとしての医師や技師等は、表示部１６０に表示された光音響画像を確認することにより、診断を行うことができる。表示画像は、ユーザーやコンピュータ１５０からの保存指示に基づいて、コンピュータ１５０内のメモリや、モダリティとネットワークで接続されたデータ管理システムなどに保存される。
【００３３】
また、コンピュータ１５０は、光音響装置に含まれる構成の駆動制御も行う。また、表示部１６０は、コンピュータ１５０で生成された画像の他にＧＵＩなどを表示してもよい。入力部１７０は、ユーザーが情報を入力できるように構成されている。ユーザーは、入力部１７０を用いて測定開始や終了、作成画像の保存指示などの操作を行うことができる。
【００３４】
以下、本実施形態に係る光音響装置の各構成の詳細を説明する。
【００３５】
（光照射部１１０）
光照射部１１０は、２つの光源として第１の光源１１１ａおよび第２の光源１１１ｂを含む。また、光照射部１１０は、第１の光源１１１ａおよび第２の光源１１１ｂから射出された光を被検体１００へ導く第１の光学系１１２ａおよび第２の光学系１１２ｂを含む。また、光照射部１１０は、第１の光源１１１ａおよび第２の光源１１１ｂの駆動を制御する第１の駆動部１１３ａおよび第２の駆動部１１３ｂを含む。
【００３６】
光の波長として４００ｎｍから１６００ｎｍ程度の範囲の波長であってもよい。血管を高解像度でイメージングする場合は、血管での吸収が大きい波長（４００ｎｍ以上、７００ｎｍ以下）を用いてもよい。生体の深部をイメージングする場合には、生体の背景組織（水や脂肪など）において典型的に吸収が少ない波長（７００ｎｍ以上、１１００ｎｍ以下）の光を用いてもよい。
【００３７】
第１、第２の光源１１１ａ、１１１ｂとしては、レーザや発光ダイオード（ＬＥＤ）を用いることができる。また、波長の変更が可能な光源であってもよい。
【００３８】
例えば、第１、第２の光源１１１ａ、１１１ｂは、周波数１ＭＨｚ以上の鋸波状の駆動波形（駆動電流）に追随して光を発することのできる半導体レーザもしくはＬＥＤを採用してもよい。
【００３９】
第１、第２の光学系１１２ａ、１１２ｂには、レンズ、ミラー、光ファイバ等の光学素子を用いることができる。乳房等を被検体１００とする場合、パルス光のビーム径を広げて照射するために、光学系の光出射部は光を拡散させる拡散板等で構成されていてもよい。一方、光音響顕微鏡においては、解像度を上げるために、第１、第２の光学系１１２ａ、１１２ｂの光出射部はレンズ等で構成し、ビームをフォーカスして照射してもよい。なお、光照射部１１０が光学系１１２ａ、１１２ｂを備えずに、第１、第２の光源１１１ａ、１１１ｂから直接被検体１００に光を照射してもよい。
【００４０】
第１、第２の駆動部１１３ａ、１１３ｂは、第１、第２の光源１１１ａ、１１１ｂを駆動する駆動電流（第１、第２の光源１１１ａ、１１１ｂに投入する電流）を生成するものである。第１、第２の駆動部１１３ａ、１１３ｂには、時間的に第１、第２の光源１１１ａ、１１１ｂへの投入電流を変化させることのできる電源を用いることができる。第１、第２の駆動部１１３ａ、１１３ｂが、第１、第２の光源１１１ａ、１１１ｂの出力を制御することにより、図１に示すような光を発生させて符号化を実現する。なお、第１、第２の駆動部１１３ａ、１１３ｂは、後述するコンピュータ１５０内の制御部１５３により制御されてもよい。また、第１、第２の駆動部１１３ａ、１１３ｂが電流値を制御する制御部を有し、この制御部が投入電流を制御してもよい。駆動電流と照射光の強度との関係については後述する。
【００４１】
（受信部１２０）
受信部１２０は、音響波を受信することにより電気信号を出力するトランスデューサと、トランスデューサを支持する支持体とを含む。
【００４２】
トランスデューサを構成する部材としては、ＰＺＴ（チタン酸ジルコン酸鉛）に代表される圧電セラミック材料や、ＰＶＤＦ（ポリフッ化ビニリデン）に代表される高分子圧電膜材料などを用いることができる。また、圧電素子以外の素子を用いてもよい。例えば、静電容量型トランスデューサ（ＣＭＵＴ：Ｃａｐａｃｉｔｉｖｅ  Ｍｉｃｒｏ−ｍａｃｈｉｎｅｄ  Ｕｌｔｒａｓｏｎｉｃ  Ｔｒａｎｓｄｕｃｅｒｓ）、ファブリペロー干渉計を用いたトランスデューサなどを用いることができる。なお、音響波を受信することにより電気信号を出力できる限り、いかなるトランスデューサを採用してもよい。また、トランスデューサにより得られる信号は時間分解信号である。つまり、トランスデューサにより得られる信号の振幅は、各時刻にトランスデューサで受信される音圧に基づく値（例えば、音圧に比例した値）を表したものである。
【００４３】
光音響波を構成する周波数成分は、典型的には１００ＫＨｚから１００ＭＨｚであり、トランスデューサとして、これらの周波数を検出することのできるものを採用することができる。
【００４４】
支持体は、１Ｄアレイ、１．５Ｄアレイ、１．７５Ｄアレイ、２Ｄアレイと呼ばれるような平面又は曲面内に、複数のトランスデューサを並べて配置してもよい。曲面内に複数のトランスデューサが並べる場合、３次元に配置されたトランスデューサアレイともいえる。
【００４５】
また、受信部１２０が、トランスデューサから出力される時系列のアナログ信号を増幅する増幅器を備えてもよい。また、受信部１２０が、トランスデューサから出力される時系列のアナログ信号を時系列のデジタル信号に変換するＡ／Ｄ変換器を備えてもよい。すなわち、受信部１２０が後述する信号収集部１４０を備えてもよい。
【００４６】
なお、音響波を様々な角度で検出できるようにするために、理想的には被検体１００を全周囲から囲むようにトランスデューサを配置してもよい。ただし、被検体１００が大きく全周囲を囲むようにトランスデューサを配置できない場合は、半球状の支持体上にトランスデューサを配置して全周囲を囲む状態に近づけてもよい。なお、トランスデューサの配置や数及び支持体の形状は被検体に応じて最適化すればよく、本発明に関してはあらゆる受信部１２０を採用することができる。
【００４７】
受信部１２０と被検体１００との間の空間は、光音響波が伝搬することができる媒質で満たしてもよい。この媒質には、音響波が伝搬でき、被検体１００やトランスデューサとの界面において音響特性が整合し、できるだけ光音響波の透過率が高い材料を採用する。例えば、この媒質には、水、超音波ジェルなどを採用することができる。
【００４８】
また、本実施形態に係る装置が、光音響画像に加えて、音響波の送受信により超音波画像も生成する場合、トランスデューサは、音響波を送信する送信手段として機能してもよい。受信手段としてのトランスデューサと送信手段としてのトランスデューサとは、単一（共通）のトランスデューサでもよいし、別々の構成であってもよい。
【００４９】
（信号収集部１４０）
信号収集部１４０は、受信部１２０から出力されたアナログ信号である電気信号を増幅するアンプと、アンプから出力されたアナログ信号をデジタル信号に変換するＡ／Ｄ変換器とを含む。信号収集部１４０は、ＦＰＧＡ（Ｆｉｅｌｄ  Ｐｒｏｇｒａｍｍａｂｌｅ  Ｇａｔｅ  Ａｒｒａｙ）チップなどで構成されてもよい。信号収集部１４０から出力されるデジタル信号は、コンピュータ１５０内の記憶部１５２に記憶される。信号収集部１４０は、Ｄａｔａ  Ａｃｑｕｉｓｉｔｉｏｎ  Ｓｙｓｔｅｍ（ＤＡＳ）とも呼ばれる。本明細書において電気信号は、アナログ信号もデジタル信号も含む概念である。なお、信号収集部１４０は、光照射部１１０の光射出部に取り付けられた光検出センサと接続されており、光が光照射部１１０から射出されたことをトリガーに、同期して処理を開始してもよい。また、信号収集部１４０は、フリーズボタンなどを用いてなされる指示をトリガーとして同期して、当該処理を開始してもよい。
【００５０】
（コンピュータ１５０）
情報処理装置としてのコンピュータ１５０は、演算部１５１、記憶部１５２、制御部１５３を含む。各構成の機能については処理フローの説明の際に説明する。
【００５１】
演算部１５１としての演算機能を担うユニットは、ＣＰＵやＧＰＵ（Ｇｒａｐｈｉｃｓ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）等のプロセッサ、ＦＰＧＡ（Ｆｉｅｌｄ  Ｐｒｏｇｒａｍｍａｂｌｅ  Ｇａｔｅ  Ａｒｒａｙ）チップ等の演算回路で構成されることができる。これらのユニットは、単一のプロセッサや演算回路から構成されるだけでなく、複数のプロセッサや演算回路から構成されていてもよい。演算部１５１は、入力部１７０から、被検体音速や保持部の構成などの各種パラメータを受けて、受信信号を処理してもよい。
【００５２】
記憶部１５２は、ＲＯＭ（Ｒｅａｄ  ｏｎｌｙ  ｍｅｍｏｒｙ）、磁気ディスクやフラッシュメモリなどの非一時記憶媒体で構成することができる。また、記憶部１５２は、ＲＡＭ（Ｒａｎｄｏｍ  Ａｃｃｅｓｓ  Ｍｅｍｏｒｙ）などの揮発性の媒体であってもよい。なお、プログラムが格納される記憶媒体は、非一時記憶媒体である。なお、記憶部１５２は、１つの記憶媒体から構成されるだけでなく、複数の記憶媒体から構成されていてもよい。
【００５３】
記憶部１５２は、後述する方法で演算部１５１により生成される光音響画像を示す画像データを保存することができる。
【００５４】
制御部１５３は、ＣＰＵなどの演算素子で構成される。制御部１５３は、光音響装置の各構成の動作を制御する。制御部１５３は、入力部１７０からの測定開始などの各種操作による指示信号を受けて、光音響装置の各構成を制御してもよい。また、制御部１５３は、記憶部１５２に格納されたプログラムコードを読み出し、光音響装置の各構成の作動を制御する。
【００５５】
コンピュータ１５０は専用に設計されたワークステーションであってもよい。また、コンピュータ１５０の各構成は異なるハードウェアによって構成されてもよい。また、コンピュータ１５０の少なくとも一部の構成は単一のハードウェアで構成されてもよい。
【００５６】
図４は、本実施形態に係るコンピュータ１５０の具体的な構成例を示す。本実施形態に係るコンピュータ１５０は、ＣＰＵ１５４、ＧＰＵ１５５、ＲＡＭ１５６、ＲＯＭ１５７、外部記憶装置１５８から構成される。また、コンピュータ１５０には、表示部１６０としての液晶ディスプレイ１６１、入力部１７０としてのマウス１７１、キーボード１７２が接続されている。
【００５７】
また、コンピュータ１５０および受信部１２０は、共通の筺体に収められた構成で提供されてもよい。また、筺体に収められたコンピュータで一部の信号処理を行い、残りの信号処理を筺体の外部に設けられたコンピュータで行ってもよい。この場合、筺体の内部および外部に設けられたコンピュータを総称して、本実施形態に係るコンピュータとすることができる。すなわち、コンピュータを構成するハードウェアが一つの筺体に収められていなくてもよい。
【００５８】
（表示部１６０）
表示部１６０は、液晶ディスプレイや有機ＥＬ（Ｅｌｅｃｔｒｏ  Ｌｕｍｉｎｅｓｃｅｎｃｅ）などのディスプレイである。コンピュータ１５０により得られた被検体情報等に基づく画像や特定位置の数値等を表示する装置である。表示部１６０は、画像や装置を操作するためのＧＵＩを表示してもよい。なお、被検体情報の表示にあたっては、表示部１６０またはコンピュータ１５０において画像処理（輝度値の調整等）を行った上で表示することもできる。
【００５９】
（入力部１７０）
入力部１７０としては、ユーザーが操作可能な、マウスやキーボードなどで構成される操作コンソールを採用することができる。また、表示部１６０をタッチパネルで構成し、表示部１６０を入力部１７０として利用してもよい。
【００６０】
なお、光音響装置の各構成はそれぞれ別の装置として構成されてもよいし、一体となった１つの装置として構成されてもよい。また、光音響装置の少なくとも一部の構成が一体となった１つの装置として構成されてもよい。
【００６１】
（被検体１００）
被検体１００は光音響装置を構成するものではないが、以下に説明する。本実施形態に係る光音響装置は、人や動物の悪性腫瘍や血管疾患などの診断や化学治療の経過観察などを目的として使用できる。よって、被検体１００としては、生体、具体的には人体や動物の乳房や各臓器、血管網、頭部、頸部、腹部、手指および足指を含む四肢などの診断の対象部位が想定される。例えば、人体が測定対象であれば、オキシヘモグロビンあるいはデオキシヘモグロビンやそれらを含む多く含む血管あるいは腫瘍の近傍に形成される新生血管などを光吸収体の対象としてもよい。また、頸動脈壁のプラークなどを光吸収体の対象としてもよい。また、メチレンブルー（ＭＢ）、インドシニアングリーン（ＩＣＧ）などの色素、金微粒子、またはそれらを集積あるいは化学的に修飾した外部から導入した物質を光吸収体としてもよい。また、刺針や刺針に付された光吸収体を観察対象としてもよい。
【００６２】
＜符号要素｛１｝に相当する光音響波＞
本実施形態に係る光音響装置を用いたときの各符号要素に相当する照射光、および、光音響波の受信信号について考察する。まず、図５または図６を用いて、符号要素｛１｝に相当する照射光、及び、光音響波の受信信号について説明する。図５及び６に示すデータは、シミュレーションにより得られたデータである。
【００６３】
図５（ａ）は、第１の光源１１１ａあるいは第２の光源１１１ｂとして波長８０８ｎｍの半導体レーザを用いた場合の、半導体レーザの電流−光出力特性を示す図である。当該半導体レーザのしきい値電流は２．５Ａ、投入電流が３０Ａのときの光出力は５０Ｗとなる。半導体レーザの場合、典型的に閾値電流以上の電流の領域では、電流−光出力特性は略線形の関係となる。すなわち、半導体レーザの場合、投入電流の時間波形が光出力（照射光の強度）の時間波形となる。
【００６４】
図５（ｂ）は、正の符号要素に対応する光を発生させるための駆動電流（第１の駆動電流）を示しており、５０ｎｓの時間で電流値を０から３０Ａまで上昇させ、１９５０ｎｓの時間で電流値を３０Ａから０Ａまで低下させている。すなわち、第１の駆動電流の中で、正の符号要素に対応するタイミングでの電流の時間変化に比べて、その他のタイミングの方が電流の時間変化が小さい。この結果、正の符号要素に対応する正の強度変調光は、正の符号要素に対応する基準タイミングでの光強度の時間変化が、その他のタイミングの光強度の時間変化よりも大きい。
【００６５】
図５（ｃ）は、図５（ｂ）の駆動電流で半導体レーザを駆動した場合の光出力を示している。前述したように、駆動電流に対して光が略線形に出力されていることが理解される。
【００６６】
図６（ａ）は、点状の光吸収体に光を照射した場合に発生する光音響波を受信帯域が無限大のトランスデューサで受信した場合の受信信号を示している。これは、図５（ｃ）の光出力曲線を時間微分したものに等しい。このように、短時間に光出力が急上昇するときに正の大きな受信信号が得られる。
【００６７】
なお、実際には、トランスデューサの受信帯域が無限大であることはあり得なく、何らかの周波数特性を有している。図６（ｂ）は、中心周波数４ＭＨｚ，６ｄＢ帯域１．６−６．４ＭＨｚの周波数特性を持つトランスデューサの受信特性を示したものである。図６（ｃ）は、図６（ｂ）に示す受信特性を有するトランスデューサで、図５（ｃ）の光照射により発生した光音響波を受信した場合の受信信号を示している。このように、トランスデューサの受信特性を考慮しても、短時間に光出力が急上昇するときに正の大きな受信信号が得られることになる。この正の大きな受信信号が、正の符号要素（例えば、符号要素｛１｝）に対応する受信信号となる。
【００６８】
図６においては、音源から受信手段までの光音響波の伝搬時間については無視している。
【００６９】
＜符号要素｛−１｝に相当する光音響波＞
詳細の説明は割愛するが、図５（ｂ）の駆動電流を時間軸で反転させた駆動電流（第２の駆動電流）で半導体レーザを駆動した場合に得られる光音響波の受信信号は、図６（ｃ）を時間軸で反転させ、かつ、信号レベルの正負を反転させた波形となる。すなわち、第２の駆動電流の中で、負の符号要素に対応するタイミングでの電流の時間変化に比べて、その他のタイミングの方が電流の時間変化が小さい。この結果、負の符号要素に対応する負の強度変調光は、負の符号要素に対応する基準タイミングでの光強度の時間変化が、その他のタイミングの光強度の時間変化よりも大きい。なお、このようにして得られる負の大きな受信信号が、負の符号要素（例えば、符号要素｛−１｝）に対応する受信信号となる。
【００７０】
＜基準タイミングでの駆動電流の時間変化＞
図５（ｂ）において、５０ｎｓの時間で電流値を０から３０Ａまで上昇させ、１９５０ｎｓの時間で電流値を３０Ａから０Ａまで低下させる例を示した。ここで、電流値を０から３０Ａまで上昇させる時間を立ち上がり時間、３０Ａから０Ａまで低下させる時間を立ち下がり時間とする。
【００７１】
図７Ａ及び図７Ｂは、符号要素｛１｝に相当する光音響波の受信信号の時間波形である。図７Ａでは、立ち上がり時間と立ち下がり時間の和を２０００ｎｓと固定し、立ち上がり時間（ｔｒ）を５０ｎｓから２５０ｎｓまで５０ｎｓステップで変えた場合の駆動電流と光音響波の受信信号の時間波形を上から順に示す。図７Ｂでは、立ち上がり時間と立ち下がり時間の和を２０００ｎｓと固定し、立ち上がり時間（ｔｒ）を３００ｎｓから５００ｎｓまで５０ｎｓステップで変えた場合の駆動電流と光音響波の受信信号の時間波形を上から順に示す。基準タイミングは、電流値を０から３０Ａまで上昇させるまでの間に位置する。ここで、立ち上がり時間は、基準タイミングの前に光量が極値となるタイミングと、基準タイミングの後に光量が極値となるタイミングとの間の期間である。
【００７２】
図８Ａ及び図８Ｂは、符号要素｛−１｝に相当する光音響波の受信信号の時間波形である。図８Ａでは、立ち上がり時間と立ち下がり時間の和を２０００ｎｓと固定し、立ち下がり時間（ｔｆ）を５０ｎｓから２５０ｎｓまで５０ｎｓステップで変えた場合の駆動電流と光音響波の受信信号の時間波形を上から順に示す。図８Ｂでは、立ち上がり時間と立ち下がり時間の和を２０００ｎｓと固定し、立ち下がり時間（ｔｆ）を３００ｎｓから５００ｎｓまで５０ｎｓステップで変えた場合の駆動電流と光音響波の受信信号の時間波形を上から順に示す。基準タイミングは、電流値を３０から０Ａまで低下させるまでの間に位置する。ここで、立ち下り時間は、基準タイミングの前に光量が極値となるタイミングと、基準タイミングの後に光量が極値となるタイミングとの間の期間である。
【００７３】
符号要素｛１｝における電流の立ち上がり時間（ｔｒ）と、その時得られる光音響波の受信信号の半値幅（ＦＷＨＭ）の関係を図９（ａ）に示す。これから、立ち上がり時間が１５０ｎｓから３００ｎｓの間ではほぼ線形に半値幅が増加するが、３００ｎｓを超えると図７Ｂからわかるように受信信号がふた山に分離するため線形から外れる。線形から外れた受信信号から画像化を行うと画像に歪が生じる。よって、中心周波数４ＭＨｚ、６ｄＢ帯域１．６−６．４ＭＨｚの周波数特性を持つトランスデューサを用いた場合、立ち上がり時間は３００ｎｓ以下であることが好ましい。ところで、立ち上がり時間の上限値はトランスデューサの受信帯域の中心周波数に反比例する傾向がある。すなわち、トランスデューサの受信帯域の中心周波数をｆＨｚとすると、立ち上がり時間は６／（５ｆ）秒以下とすることが好ましい。ここでは、基準タイミングにおける立ち上り時間の上限について説明したが、基準タイミングにおける立ち下り時間の上限についても同様である。
【００７４】
ある符号系列の符号を発生させる場合、符号要素｛１｝に相当する受信信号と符号要素｛−１｝に相当する受信信号とが正負反転の関係にあることが理想である。よって、ある符号系列において、電流の立ち上がり時間（ｔｒ）と符号要素｛−１｝における電流の立ち下がり時間（ｔｆ）を一致させることが好ましい。立ち上がり時間（ｔｒ）に対する、符号要素｛１｝に相当する受信信号と符号要素｛−１｝に相当する受信信号との相関値の変化を図９（ｂ）に示す。受信信号の波形が完全に反転している場合、相関値が−１になり、相関値が−１に近いほど理想状態に近いことを意味する。図９（ｂ）より、立ち上がり時間（ｔｒ）が大きくなると、相関値の−１からのずれが大きくなることが理解される。相関値の絶対値が０．９以上である範囲を許容すると、立ち上がり時間は４５０ｎｓ以下であることが好ましい。ところで、立ち上がり時間の上限値はトランスデューサの受信帯域の中心周波数に反比例する傾向がある。すなわち、トランスデューサの受信帯域の中心周波数をｆＨｚとすると、立ち上がり時間は９／（５ｆ）秒以下とすることが好ましい。ここでは、基準タイミングにおける立ち上り時間の上限について説明したが、基準タイミングにおける立ち下り時間の上限についても同様である。
【００７５】
＜光音響画像の生成フロー＞
次に、本実施形態に係る光音響装置を用いた符号化・復号化により光音響画像を生成する方法（情報取得方法）を説明する。
【００７６】
Ｓ１．光照射部１１０が、第１の符号系列に対応する、基準タイミングにおいて第１の光量変化量を有する第１の強度変調光を被検体１００に照射する。また、光照射部１１０が、第１の符号系列とは異なる第２の符号系列に対応する、基準タイミングにおいて第２の光量変化量を有する第２の強度変調光を被検体１００に照射する。第１の強度変調光及び第２の強度変調光は、符号化された光として機能する。なお、第１の強度変調光の照射期間と第２の強度変調光の照射期間とは、少なくとも一部が重なっていてもよい。その場合も後述する復号処理によって、第１の光量変化量に対応する復号信号と、第２の光量変化量に対応する復号信号とに分離することができる。
【００７７】
Ｓ２．符号化された光によって発生する光音響波を、受信部１２０に含まれる複数のトランスデューサが受信し、第１の受信信号を出力する。
【００７８】
Ｓ３．光照射部１１０が、第３の符号系列に対応する、基準タイミングにおいて第１の光量変化量を有する第３の強度変調光と、第３の符号系列とは異なる第４の符号系列に対応する、基準タイミングにおいて第２の光量変化量を有する第４の強度変調光とを被検体１００に照射する。なお、第１の強度変調光の照射期間と第２の強度変調光の照射期間とは、少なくとも一部が重なっていてもよい。その場合も後述する復号処理によって、第１の光量変化量に対応する復号信号と、第２の光量変化量に対応する復号信号とに分離することができる。
【００７９】
Ｓ４．符号化された光によって発生する光音響波を、受信部１２０に含まれる複数のトランスデューサが受信し、第２の受信信号を出力する。
【００８０】
Ｓ５．演算部１５１は、複数のトランスデューサから出力された第１、第２の受信信号に対して復号処理を行い、トランスデューサ毎に復号された受信信号（復号信号）を生成する。
【００８１】
Ｓ６．演算部１５１は、複数のトランスデューサに対応する複数の復号信号を用いて、光音響画像を生成する。
【００８２】
具体的な符号化・復号化の方法に関しては、実施例で後述する。
【００８３】
＜再構成手法＞
演算部１５１は、複数の復号信号を計算空間に逆投影（単純逆投影）して、画像データを生成することができる。すなわち、演算部１５１は、時間信号である復号信号を空間分布データに変換してもよい。例えば、演算部１５１は、複数の復号信号を整相加算して深さ方向に直線状の画像データ（１ライン分の画像データ）を取得してもよい。また、演算部１５１は、この処理を複数のラインに対して実行することにより、二次元または三次元の画像データを生成してもよい。なお、整相加算により得られた空間分布データに対して包絡線処理を行って、画像データを生成してもよい。
【００８４】
ところで、ＰＡＴの画像再構成手法として、Ｕｎｉｖｅｒｓａｌ  Ｂａｃｋ  Ｐｒｏｊｅｃｔｉｏｎ（ＵＢＰ）法が知られている。これは、受信部で得られた受信信号を時間微分し、かつ正負を反転したデータを逆投影することで光音響画像を得る手法である。これは、インパルス状のパルス光を照射したときに発生する光音響波がＮ−ｓｈａｐｅと呼ばれるアルファベットのＮ字のような形状を持つ場合に採用できる手法である。
【００８５】
ところで、本実施形態で発生する光音響波は、Ｎ−ｓｈａｐｅの前半部と後半部を分離し、前半部を符号要素｛１｝に相当する光音響波、後半部を符号要素｛−１｝に相当する光音響波としたものと便宜上理解できる。そのため、本実施形態のように符号化し、復号された受信信号では、ＵＢＰ法を適用しても正しい結果は得られない。そのため、本実施形態においては、演算部１５１は、復号された受信信号に対してＵＢＰ法で行われる前処理（時間微分処理等）を行うことなく、位相を整合した後に加算する整相加算処理を行い、逆投影することが好ましい。本明細書では、復号された受信信号に対してＵＢＰ法で行われる前処理を実行することなく逆投影する再構成手法を、単純逆投影と呼ぶ。その他、信号データを３次元のボリュームデータに変換する再構成アルゴリズムとしては、タイムドメインでの逆投影法、フーリエドメインでの逆投影法、モデルベース法（繰り返し演算法）などのあらゆる手法を採用してもよい。
【実施例１】
【００８６】
＜相補符号を適用した符号化・復号化＞
相補符号を適用した、基準タイミングにおいて異なる光量変化量を有する複数の強度変調光の光照射による符号化・復号化について説明する。
【００８７】
符号長Ｎの２つの符号系列｛ａｉ｝、｛ｃｉ｝、（ｉ＝１からＮ、Ｎは符号長、各符号要素は１もしくは−１）において、それぞれの自己相関関数の和が、ピークにおいて２Ｎとなり、ピーク以外のすべての点でＯとなる符号系列の組を相補符号と呼ぶ。
【００８８】
ここで、自己相関関数を、
【００８９】
【数１】
2019068914.tif 000003
のように表す。
【００９０】
例えば、｛ａｉ｝＝｛１，１｝、｛ｃｉ｝＝｛１，−１｝の組は相補符号である。
【００９１】
すなわち、（ａ＊ａ）＝｛１，２，１｝、（ｃ＊ｃ）＝｛−１，２，−１｝、（ａ＊ａ）＋（ｃ＊ｃ）＝｛０，４，０｝となる。
【００９２】
符号長が２のｎ乗、あるいは５×２のｎ乗の場合、相補符号が存在することが知られている（ｎは自然数）。
【００９３】
相補符号である、ある符号系列の組に対し、それぞれの相互相関関数の和が０となる相補符号が存在する。このように、２組の相補符号において、それぞれの相互相関関数の和が０となる関係のことを便宜上「完全直交の関係」と呼ぶ。
【００９４】
ここで、相互相関関数を、
【００９５】
【数２】
2019068914.tif 000004
のように表す。
【００９６】
例えば、相補符号である第１の符号系列｛ａｉ｝＝｛１，１｝、第３の符号系列｛ｃｉ｝＝｛１，−１｝の組に対し、第２の符号系列｛ｂｉ｝＝｛１，−１｝、第４の符号系列｛ｄｉ｝＝｛１，１｝の組は「完全直交の関係」を満たす。すなわち、（ａ＊ｂ）＝｛−１，０，１｝、（ｃ＊ｄ）＝｛１，０，−１｝、（ａ＊ｂ）＋（ｃ＊ｄ）＝０となる。このとき、（ｂ＊ａ）＋（ｄ＊ｃ）＝０も成立する。
【００９７】
符号長８の場合、
第１の符号系列｛ａｉ｝＝｛１，１，−１，１，−１，−１，−１，１｝
第２の符号系列｛ｂｉ｝＝｛１，１，−１，１，１，１，１，−１｝
第３の符号系列｛ｃｉ｝＝｛１，−１，−１，−１，−１，１，−１，−１｝
第４の符号系列｛ｄｉ｝＝｛１，−１，−１，−１，１，−１，１，１｝
とする。このとき、
（ａ＊ａ）＝｛１，０，−３，０，−１，０，−１，８，−１，０，−１，０，−３，０，１｝
（ｃ＊ｃ）＝｛−１，０，３，０，１，０，１，８，１，０，１，０，３，０，−１｝
（ａ＊ａ）＋（ｃ＊ｃ）＝｛０，０，０，０，０，０，０，１６，０，０，０，０，０，０，０｝
となり、第１の符号系列｛ａｉ｝と第３の符号系列｛ｃｉ｝の組は相補符号である。同様に、
（ｂ＊ｂ）＝｛−１，０，３，０，１，０，１，８，１，０，１，０，３，０，−１｝
（ｄ＊ｄ）＝｛１，０，−３，０，−１，０，−１，８，−１，０，−１，０，−３，０，１｝
（ｂ＊ｂ）＋（ｄ＊ｄ）＝｛０，０，０，０，０，０，０，１６，０，０，０，０，０，０，０｝
となり、第２の符号系列｛ｂｉ｝と第４の符号系列｛ｄｉ｝の組も相補符号である。また、
（ａ＊ｂ）＝｛−１，０，３，０，３，０，−１，０，−３，０，１，０，−３，０，１｝
（ｃ＊ｄ）＝｛１，０，−３，０，−３，０，１，０，３，０，−１，０，３，０，−１｝
（ａ＊ｂ）＋（ｃ＊ｄ）＝０
（ｂ＊ａ）＝｛１，０，−３，０，１，０，−３，０，−１，０，３，０，３，０，−１｝
（ｄ＊ｃ）＝｛−１，０，３，０，−１，０，３，０，１，０，−３，０，−３，０，１｝
（ｂ＊ａ）＋（ｄ＊ｃ）＝０
となり、第１の相補符号である符号系列｛ａｉ｝及び｛ｃｉ｝の組と、第２の相補符号である符号系列｛ｂｉ｝及び｛ｄｉ｝の組は「完全直交の関係」を満たす。
【００９８】
このように「完全直交の関係」を満たす符号系列の組を用いると、以下のようなことが実現できる。すなわち、
・第１の符号系列｛ａｉ｝で符号化した信号Ａと第３の符号系列｛ｃｉ｝で符号化した信号Ｃがあるとき、信号Ａを第１の符号系列｛ａｉ｝で復号化した信号と信号Ｃを第３の符号系列｛ｃｉ｝で復号化した信号との和はデルタ関数となる。
・信号Ａを第２の符号系列｛ｂｉ｝で復号化した信号と信号Ｃを第４の符号系列｛ｄｉ｝で復号化した信号との和は０となりキャンセルできる。
・第２の符号系列｛ｂｉ｝で符号化した信号Ｂと第４の符号系列｛ｄｉ｝で符号化した信号Ｄがあるとき、信号Ｂを第２の符号系列｛ｂｉ｝で復号化した信号と信号Ｄを第４の符号系列｛ｄｉ｝で復号化した信号との和はデルタ関数となる。
・信号Ｂを第１の符号系列｛ａｉ｝で復号化した信号と信号Ｄを第３の符号系列｛ｃｉ｝で復号化した信号との和は０となりキャンセルできる。
【００９９】
このような符号系列を、基準タイミングにおいて異なる光量変化量を有する複数の強度変調光の光照射を行う光音響装置に適用する。すなわち、本実施例に係る光音響装置は、第１の符号系列｛ａｉ｝および第３の符号系列｛ｃｉ｝で符号化した、基準タイミングにおいて第１の光量変化量を有する強度変調光を照射する。さらに、本実施例に係る光音響装置は、第２の符号系列｛ｂｉ｝および第４の符号系列｛ｄｉ｝で符号化した、基準タイミングにおいて第２の光量変化量を有する強度変調光を照射する。
【０１００】
この場合、光の照射時間が重複していても、復号処理により、第１の光量変化量に対応する信号と、第２の光量変化量に対応する信号とを分離して取得することができる。
【０１０１】
図３に示す光音響装置を用いて、基準タイミングにおける光量変化量の異なる複数の強度変調光のそれぞれを照射する場合を説明する。本実施例では、第１の光源１１１ａ、第２の光源１１１ｂともに、波長８０８ｎｍ、光出力最大５０Ｗの半導体レーザを用いる。
【０１０２】
受信部１２０には、中心周波数４ＭＨｚ、６ｄＢ帯域１．６−６．４ＭＨｚの周波数特性を持つ圧電素子からなるリニアアレイを用いる。受信部１２０と被検体１００との間は音響マッチングのために超音波ジェルで埋められている。
【０１０３】
本実施例では、符号長８の相補符号を用いる。具体的には、
第１の符号系列｛ａｉ｝＝｛１，１，−１，１，−１，−１，−１，１｝
第２の符号系列｛ｂｉ｝＝｛１，１，−１，１，１，１，１，−１｝
第３の符号系列｛ｃｉ｝＝｛１，−１，−１，−１，−１，１，−１，−１｝
第４の符号系列｛ｄｉ｝＝｛１，−１，−１，−１，１，−１，１，１｝
である。
【０１０４】
本実施例では、図１０に示すシーケンスで、基準タイミングにおいて第１の光量変化量を有する強度変調光、および、基準タイミングにおいて第２の光量変化量を有する強度変調光の照射により発生する光音響波の受信信号の取得を行い、符号化が行われる。
【０１０５】
まず、光照射部１１０が、第１の符号系列｛ａｉ｝に対応する、基準タイミングにおいて第１の光量変化量を有する強度変調光と、第２の符号系列｛ｂｉ｝に対応する、基準タイミングにおいて第２の光量変化量を有する強度変調光とを被検体１００に照射する。これら２つの強度変調光の照射は同期されている。そして、受信部１２０は、それらの強度変調光の照射により発生した光音響波を受信し、受信信号Ｓ１を出力する。
【０１０６】
次に、光照射部１１０が、第３の符号系列｛ｃｉ｝に対応する、基準タイミングにおいて第１の光量変化量を有する強度変調光と、第４の符号系列｛ｄｉ｝に対応する、基準タイミングにおいて第２の光量変化量を有する強度変調光とを被検体１００に照射する。この２つの強度変調光の照射は同期されている。そして、受信部１２０は、それらの強度変調光の照射により発生した光音響波を受信し、受信信号Ｓ２を出力する。なお、受信信号Ｓ１を取得する期間と、受信信号Ｓ２を取得する期間とは、時間的に分離されていることが好ましい。
【０１０７】
以下、本実施例に係る符号化及び復号化について詳述する。
【０１０８】
図１１のように、光照射部１１０から発せられた、基準タイミングにおいて第１の光量変化量を有する強度変調光１１５ａ、および、基準タイミングにおいて第２の光量変化量を有する強度変調光１１５ｂを被検体１００に照射する場合を考える。
【０１０９】
被検体１００中、受信部１２０に近い側に直径０．１ｍｍの球状の光吸収体１９１、受信部１２０から遠い側に直径１ｍｍの球状の光吸収体１９２があると仮定する。
【０１１０】
制御部１５３は、第１の符号系列｛ａｉ｝に関する情報を第１の駆動部１１３ａに送信し、第２の符号系列｛ｂｉ｝に関する情報を第２の駆動部１１３ｂに送信する。
【０１１１】
図１２（ａ）は、第１の符号系列｛ａｉ｝に関する情報に基づいて、第１の駆動部１１３ａが生成した駆動電流である。ここで、基準タイミングの時間間隔（符号要素の周期に相当する）は２０００ｎｓである。正の符号要素に対応する第１の駆動電流の立ち上がり時間は５０ｎｓである。負の符号要素に対応する第２の駆動電流の立ち下がり時間は５０ｎｓである。すなわち、本実施例において、第１の光量変化量は、５０ｎｓの立ち上がり時間または立ち下がり時間内の光量変化量に相当する。
【０１１２】
図１２（ａ）に示す駆動電流によって生成された変調光を光吸収体１９１、１９２に照射した際に発生する光音響波を受信部１２０で受信して得られる受信信号は、図１２（ｂ）のような波形となる。実際には、光吸収体１９１から受信部１２０まで光音響波が伝搬する時間の分だけ時間がずれるが、その時間は無視して表記している。ここで、雑音の抑圧効果を説明するため、平均値０、標準偏差０．１の雑音を付加している。
【０１１３】
また、図１２（ｃ）は、第２の符号系列｛ｂｉ｝に関する情報に基づいて、第２の駆動部１１３ｂが生成した駆動電流である。ここで、図１２（ａ）と同様、基準タイミングの時間間隔（符号要素の周期に相当する）は２０００ｎｓである。正の符号要素に対応する第１の駆動電流の立ち上がり時間は３００ｎｓである。負の符号要素に対応する第２の駆動電流の立ち下がり時間は３００ｎｓである。すなわち、本実施例において、第２の光量変化量は、３００ｎｓの立ち上がり時間または立ち下がり時間内の光量変化量に相当する。
【０１１４】
図１２（ｃ）に示す駆動電流によって生成された変調光を光吸収体１９１、１９２に照射した際に発生する光音響波を受信部１２０で受信して得られる受信信号は、図１２（ｄ）のような波形となる。実際には、光吸収体１９１から受信部１２０まで光音響波が伝搬する時間の分だけ時間がずれるが、その時間は無視して表記している。
【０１１５】
第１の光源１１１ａの光出力と第２の光源１１１ｂの光出力とを同期させた場合（略同じタイミング光照射した場合）、受信部１２０が受信して得られる受信信号は、図１２（ｂ）と図１２（ｄ）との和となる。得られた受信信号をＳ１（ｔ）とし、図１３に示す。
【０１１６】
続いて、制御部１５３は、第３の符号系列｛ｃｉ｝に関する情報を第１の駆動部１１３ａに送信し、第４の符号系列｛ｄｉ｝に関する情報を第２の駆動部１１３ｂに送信する。
【０１１７】
図１４（ａ）は、第３の符号系列｛ｃｉ｝に関する情報に基づいて、第１の駆動部１１３ａが生成した駆動電流である。ここで、基準タイミングの時間間隔（符号要素の周期に相当する）は２０００ｎｓである。正の符号要素に対応する第１の駆動電流の立ち上がり時間は５０ｎｓである。負の符号要素に対応する第２の駆動電流の立ち下がり時間は５０ｎｓである。前述したように、本実施例において、第１の光量変化量は、５０ｎｓの立ち上がり時間または立ち下がり時間内の光量変化量に相当する。
【０１１８】
図１４（ａ）に示す駆動電流によって生成された変調光を光吸収体１９１、１９２に照射した際に発生する光音響波を受信部１２０で受信して得られる受信信号は、図１４（ｂ）のような波形となる。実際には、光吸収体１９１から受信部１２０まで光音響波が伝搬する時間の分だけ時間がずれるが、その時間は無視して表記している。
【０１１９】
また、図１４（ｃ）は、第４の符号系列｛ｄｉ｝に関する情報に基づいて、第２の駆動部１１３ｂが生成した駆動電流である。ここで、図１４（ａ）と同様、基準タイミングの時間間隔（符号要素の周期に相当する）は２０００ｎｓである。正の符号要素に対応する第１の駆動電流の立ち上がり時間は３００ｎｓである。負の符号要素に対応する第２の駆動電流の立ち下がり時間は３００ｎｓである。前述したように、本実施例において、第２の光量変化量は、３００ｎｓの立ち上がり時間または立ち下がり時間内の光量変化量に相当する。
【０１２０】
図１４（ｃ）に示す駆動電流によって生成された変調光を光吸収体１９１、１９２に照射した際に発生する光音響波を受信部１２０で受信して得られる受信信号は、図１４（ｄ）のような波形となる。実際には、光吸収体１９１から受信部１２０まで光音響波が伝搬する時間の分だけ時間がずれるが、その時間は無視して表記している。
【０１２１】
第１の光源１１１ａの光出力と第２の光源１１１ｂの光出力とを同期させた場合（略同じタイミング光照射した場合）、受信部１２０が受信して得られる受信信号は、図１４（ｂ）と図１４（ｄ）との和となる。得られた受信信号をＳ２（ｔ）とし、図１５に示す。
【０１２２】
コンピュータ１５０内の演算部１５１が行う、符号化された受信信号の復号方法について説明する。
【０１２３】
基準タイミングの時間間隔をΔｔとしたとき、演算部１５１は、受信信号Ｓ１及びＳ２に対して（式３）にしたがった復号処理を行うことにより、第１の光量変化量に対応する復号信号ＤＳ１（ｔ）を得る。また、演算部１５１は、受信信号Ｓ１及びＳ２に対して（式４）にしたがった復号処理を行うことにより、第２の光量変化量に対応する復号信号ＤＳ２（ｔ）を得る。
【０１２４】
【数３】
2019068914.tif 000005
【０１２５】
復号信号ＤＳ１は、第１の光量変化量を有する光を用いて光音響波を発生させた場合に得られる受信信号に相当する。また、復号信号ＤＳ２は、第２の光量変化量を有する光を用いて光音響波を発生させた場合に得られる受信信号に相当する。
【０１２６】
図１３の受信波形に対して（式３）の右辺第一項の復号処理を行った結果、図１６（ａ）のような復号された受信信号を得ることができる。また、図１５の受信波形に対して（式３）の右辺第二項の復号処理を行った結果、図１６（ｂ）のような復号された受信信号を得ることができる。そして、図１６（ａ）と図１６（ｂ）の和は、図１６（ｃ）の波形（復号信号ＤＳ１）となる。図１６（ｄ）は、図１６（ｃ）の時間０の近傍を拡大したものである。図１６（ｃ）、図１６（ｄ）は基準タイミングでの立ち上がり時間（立下り時間）５０ｎｓに対応する復号信号ＤＳ１である。復号信号ＤＳ１は、第１の光量変化量を有する光（立ち上がり時間が５０ｎｓの光）を用いて光音響波を発生させた場合に得られる受信信号に相当する。図１６（ｃ）、図１６（ｄ）から、符号化及び復号化によりＳＮ比が向上していることがわかる。
【０１２７】
また、図１３の受信波形に対して（式４）の右辺第一項の復号処理を行った結果、図１７（ａ）のような復号された受信信号を得ることができる。また、図１５の受信波形に対して（式４）の右辺第二項の復号処理を行った結果、図１７（ｂ）のような復号された受信信号を得ることができる。そして、図１７（ａ）と図１７（ｂ）の和は、図１７（ｃ）の波形（復号信号ＤＳ２）となる。図１７（ｄ）は、図１７（ｃ）の時間０の近傍を拡大したものである。図１７（ｃ）、図１７（ｄ）は基準タイミングでの立ち上がり時間（立下り時間）３００ｎｓに対応する復号信号ＤＳ２である。復号信号ＤＳ２は、第２の光量変化量を有する光（立ち上がり時間が３００ｎｓの光）を用いて光音響波を発生させた場合に得られる受信信号に相当する。図１７（ｃ）、図１７（ｄ）から、符号化及び復号化によりＳＮ比が向上していることがわかる。
【０１２８】
また、図１６（ｄ）と図１７（ｄ）を比較すると、図１６（ｄ）では、小さい光吸収体の信号が強調され、図１７（ｄ）では、大きい光吸収体の信号が強調されていることがわかる。このように、符号化に用いる光の立ち上がり時間（立ち下がり時間）を制御することにより、強調させる周波数成分を制御させることができる。
【０１２９】
本実施例により、互いに異なる光量変化量を有する複数の光を同じ期間に符号化して照射して、それを復号化することで異なる光量変化量に対応する光音響信号を分離して受信することが可能となった。
【０１３０】
演算部１５１は、このようにして得られた復号信号を用いて光音響画像を生成することでＳＮ比を向上させた光音響画像を得ることが可能となる。受信部１２０が複数のトランスデューサを備える場合、各トランスデューサから出力された受信信号に対して復号処理を行い、トランスデューサ毎に復号信号が生成される。演算部１５１は、複数のトランスデューサに対応する複数の復号信号を用いて、前述した再構成手法等により光音響画像を生成することができる。
【０１３１】
なお、演算部１５１は、小さな（細い）光吸収体が強調された第１の光音響画像を生成し、表示部１６０に表示させる第１のモードを実行可能であってもよい。また、演算部１５１は、大きな（太い）光吸収体が強調された第２の光音響画像を生成し、表示部１６０に表示させる第２のモードを実行可能であってもよい。また、演算部１５１は、大小両方の光吸収体を表現した第３の光音響画像を生成し、表示部１６０に表示させる第３のモードを実行可能であってもよい。ユーザーが入力部１７０を用いて、これらのモードの中から所望のモードを選ぶことができてもよい。演算部１５１は、ユーザーからの指示に基づいて、ユーザーの所望するモードを実行してもよい。また、演算部１５１が、これらのモードで得られた光音響画像を所定の順番で切り替えて表示部１６０に表示させてもよい。また、演算部１５１が、それぞれのモードで生成された光音響画像を並べて表示部１６０に表示させてもよい。ユーザーが入力部１７０を用いて、並列表示させる画像を生成するモードを選択してもよい。
【０１３２】
演算部１５１は、第１のモードにおいて、復号信号のうち、大きな光量変化量（第１の光量変化量）に対応する復号信号（ＤＳ１）を選択的に用いて、第１の光音響画像を生成することができる。第１のモードでは、高周波成分が多く含まれる復号信号（ＤＳ１）を用いて画像化するため、小さな（細い）光吸収体が強調されやすい画像が得られる。なお、演算部１５１は、第１のモードにおいて、小さな光量変化量（第２の光量変化量）に対応する復号信号（ＤＳ２）についても用いて画像化を行ってもよい。ただし、この場合、大きな光量変化量に対応する復号信号（ＤＳ１）に対する重みを、小さな光量変化量に対応する復号信号（ＤＳ２）に対する重みよりも大きくすることにより、小さな（細い）光吸収体が強調された第１の光音響画像を生成することができる。
【０１３３】
また、演算部１５１は、第２のモードにおいて、復号信号のうち、小さな光量変化量（第２の光量変化量）に対応する復号信号（ＤＳ２）を選択的に用いて、第２の光音響画像を生成することができる。第２のモードでは、低周波成分が多く含まれる復号信号（ＤＳ２）を用いて画像化するため、大きな（太い）光吸収体が強調されやすい画像が得られる。なお、演算部１５１は、第２のモードにおいて、大きな光量変化量（第１の光量変化量）に対応する復号信号（ＤＳ１）についても用いて画像化を行ってもよい。ただし、この場合、小さな光量変化量に対応する復号信号（ＤＳ２）に対する重みを、大きな光量変化量に対応する復号信号（ＤＳ１）に対する重みよりも大きくすることにより、大きな（太い）光吸収体が強調された第２の光音響画像を生成することができる。
【０１３４】
また、演算部１５１は、第３のモードにおいて、大きな光量変化量（第１の光量変化量）に対応する復号信号（ＤＳ１）、及び、小さな光量変化量（第２の光量変化量）に対応する復号信号（ＤＳ２）を用いて、第３の光音響画像を生成することができる。演算部１５１は、両復号信号（ＤＳ１及びＤＳ２）を用いた再構成により第３の光音響画像を生成してもよい。また、演算部１５１は、大きな光量変化量に対応する復号信号（ＤＳ１）を用いた再構成によって得られた光音響画像と、小さな光量変化量に対応する復号信号（ＤＳ２）を用いた再構成によって得られた光音響画像とを合成することにより、第３の光音響画像を生成してもよい。演算部１５１は、入力部１７０を介したユーザーの指示に基づいて、大きな光量変化量に対応する復号信号（ＤＳ１）に対する第１の重みと、小さな光量変化量に対応する復号信号（ＤＳ２）に対する第２の重みとを変更可能であってもよい。これにより、ユーザーが強調させて観察したい光吸収体に適した重みを、ユーザーの指示により決定することができる。このとき、ユーザーの指示に応じた重みを適用した第３の光音響画像を表示装置１６０に表示させることにより、ユーザーが指示の結果を確認しながら、重みを調整することができる。
【０１３５】
なお、重みづけは、上述したように復号信号に対して行われてもよいし、復号信号から生成された光音響画像に対して行われてもよい。本明細書において、「復号信号に対する重みづけ」には、復号信号に重みづけする場合も、復号信号から生成された光音響画像に重みづけすることも含まれる。なお、復号信号や光音響画像をそのままその後の処理に利用することは、重み１が割り当てられたことに相当する。
【０１３６】
また、受信部１２０が互いに異なる受信帯域を有する異なるトランスデューサを備えていてもよい。例えば、受信部１２０が、中心周波数が１０ＭＨｚの受信帯域（第１の受信帯域）を有する第１のトランスデューサと、中心周波数２ＭＨｚの受信帯域（第２の受信帯域）を有する第２のトランスデューサとの両方を備える場合を考える。この場合、演算部１５１は、中心周波数の高い第１のトランスデューサから出力された受信信号（符号信号）に対して、大きな光量変化量（第１の光量変化量）に対応する復号信号（ＤＳ１）を復号する処理を選択的に実行してもよい。一方、演算部１５１は、中心周波数の低い第２のトランスデューサから出力された受信信号（符号信号）に対して、小さな光量変化量（第２の光量変化量）に対応する復号信号（ＤＳ２）を復号する処理を選択的に実行してもよい。このようにトランスデューサの受信帯域に応じて適用する復号処理を変更することにより、トランスデューサの受信帯域に適した周波数成分が多く含まれる信号を復号することができる。
【０１３７】
なお、演算部１５１は、いずれのトランスデューサから出力された受信信号からも、各光量変化量に対応する復号信号を生成してもよい。この場合、演算部１５１は、第１のトランスデューサに対応する復号信号のうち、大きな光量変化量（第１の光量変化量）に対応する復号信号（ＤＳ１）を選択的に用いて、光音響画像を生成してもよい。また、第２のトランスデューサに対応する復号信号のうち、小さな光量変化量（第１の光量変化量）に対応する復号信号（ＤＳ２）を選択的に用いて、光音響画像を生成してもよい。このように、生成された復号信号のうち、トランスデューサの受信帯域に対応する復号信号を選択的に用いて光音響画像を生成してもよい。
【０１３８】
演算部１５１は、受信部１２０からトランスデューサの受信帯域に関する情報を取得し、各トランスデューサに適用する復号処理を決定してもよい。また、演算部１５１は、ＨＩＳやＲＩＳから撮影オーダー情報を取得し、撮影オーダー情報から撮影に使用されるトランスデューサを推定し、トランスデューサの受信帯域に対応する復号処理を決定してもよい。
【０１３９】
なお、各トランスデューサに対応する復号信号（ＤＳ１及びＤＳ２）のそれぞれに重みをつけることにより、トランスデューサの受信帯域と撮影対象とに適した光音響画像を生成することができる。重みは予め設定されていてもよいし、ユーザーが入力部を用いて重みを指示してもよい。
【０１４０】
なお、重みづけは、上述したように復号信号に対して行われてもよいし、復号信号から生成された光音響画像に対して行われてもよい。本明細書において、「復号信号に対する重みづけ」には、復号信号に重みづけする場合も、復号信号から生成された光音響画像に重みづけすることも含まれる。なお、復号信号や光音響画像をそのままその後の処理に利用することは、重み１が割り当てられたことに相当する。
【０１４１】
本実施例では、光照射による符号化の際に、負の符号要素に対応する光照射（負の強度変調光）を行っているため、負の符号要素を含む符号系列に基づいて精度良く符号化することができる。そのため、本実施例によれば、このように符号化された信号を、負の符号要素を含む符号系列に基づいた復号化（例えば、（式３）、（式４）に示す復号処理）によって精度良く復号することができる。このように、負の符号要素に対応する光照射を行うことにより、負の符号要素を０として光照射しない場合と比べて、精度良く復号することができる。
【０１４２】
また、本実施例では、第１の光源１１１ａからの照射光と第２の光源１１１ｂからの照射光とをほぼ同じタイミングに被検体に照射することで、受信時間を共通化している。その結果、２つの光量変化量を有する光による光音響波を時間的に分離して受信する場合と比べ、より少ない時間でＳＮ比を向上させることが可能となる。
【０１４３】
本実施例において、２つの光量変化量を有する光を同期して被検体に照射する場合を示したが、同じタイミングに光照射することは必須ではない。ただし、測定時間を短くするためには、複数の強度変調光による光音響波の受信期間の少なくとも一部を重ねることが好ましい。
【０１４４】
また、複数の強度変調光による光音響波を時間的に分離して受信する場合、その時間の間に被検体が動いたりすると信号に時間的なずれが発生する。一方、本実施例の手法によれば、複数の強度変調光を照射する期間を重ねることで、被検体が動いた場合の信号の時間的なずれを小さくできる。
【０１４５】
本実施例に係る符号要素の基準タイミングの時間間隔の上限について説明する。
【０１４６】
１回の受信信号の取得に要する時間は、受信部から見て、被検体中の観察領域（再構成される領域）の最も遠い部分で発生した光音響波が受信部に到達するのに要する時間と等しい。この時間をＴｔｏｆとする。
【０１４７】
本実施例では、一つの光量変化量に対して符号長８の符号系列を２つ用いており、復号された受信信号においては、信号レベルは１６倍、雑音レベルは４倍となっている。そのため、ＳＮ比は４倍向上することになる。
【０１４８】
一般的に知られるインパルス状のパルス光を照射して光音響波を発生させる手法で同じＳＮ比の向上を得るには、単純に受信信号を１６回取得しそれを平均する必要がある。被検体中の光の伝搬の時間は短いので無視すると、一般的な手法において、受信信号を１６回取得するのに要する時間は１６Ｔｔｏｆとなる。さらに、２つの光量変化量に対応する受信信号を得る必要があるので、測定に要する時間は３２Ｔｔｏｆとなる。
【０１４９】
第１の符号系列｛ａｉ｝に相当する受信信号の取得に要する時間は、第１の符号系列｛ａｉ｝に相当する光を照射する時間に、最後の符号要素に対応する光によって発生した光音響波が受信部に到達するまでの時間を加えたものとなる。すなわち、その時間は７Δｔ＋Ｔｔｏｆとなる。第２〜第４の符号系列｛ｂｉ｝、｛ｃｉ｝、及び｛ｄｉ｝に相当する受信信号の取得に要する時間も同一である。そのため、第１〜第４の符号系列｛ａｉ｝、｛ｂｉ｝、｛ｃｉ｝、及び｛ｄｉ｝に対応する受信信号の取得を順番に（直列に）実行する場合、第１の光量変化量を有する光及び第２の光量変化量を有する光に起因する受信信号の取得に要する時間は２８Δｔ＋４Ｔｔｏｆとなる。
【０１５０】
本実施例では、第１の光量変化量を有する光と第２の光量変化量を有する光とを同時に照射し、受信信号の取得を同時に行っている。これにより、本実施例による受信信号の取得に要する時間は１４Δｔ＋２Ｔｔｏｆとなり、各光量変化量を有する光に対応する信号を時間的に分離する手法と比べて、受信信号の取得に要する時間は削減される。
【０１５１】
本実施例による手法のほうが受信信号の取得に要する時間が短ければ、一般的な手法と比較してＳＮ比の向上の効果が高いといえる。この条件は、１４Δｔ＜３０Ｔｔｏｆである。符号長Ｎを用いて一般化すると、この条件は、
【０１５２】
【数４】
2019068914.tif 000006
【０１５３】
Ｎがある程度大きいとき、Δｔ＜２Ｔｔｏｆであることが望ましい。すなわち、基準タイミングの時間間隔は、受信部から見て被検体の観察領域の最も遠い部分で発生した光音響波が受信部に到達するのに要する時間の２倍よりも小さいことが好ましい。例えば、受信部と被検体の観察領域の最も遠い部分との距離が５ｃｍ、被検体内の音速が１５００ｍ／ｓの場合、被検体の観察領域の最も遠い部分で発生した光音響波が受信部に到達するまでの時間は３３μｓとなる。この場合、基準タイミングの時間間隔は６６μｓよりも短くすることが好ましい。符号長８の場合は、（式５）より、７１μｓよりも短くすることが好ましい。制御部１５３は、受信部と観察領域の最も遠い部分との距離をｄ、被検体内の音速をｃとしたときに、式５中のＴｔｏｆにｄ／ｃを代入して、Δｔを計算してもよい。
【０１５４】
なお、制御部１５３は、ユーザーが入力部１７０を用いて指示した関心領域を表す情報に基づいて、関心領域の最も遠い部分で発生した光音響波が受信部に到達するまでの時間よりも短い基準タイミングの時間間隔となるように変更してもよい。また、制御部１５３は、ユーザーの指示や計算により決定された被検体の音速に応じて、最も遠い部分で発生した光音響波が受信部に到達するまでの時間よりも短い基準タイミングの時間間隔となるように変更してもよい。
【０１５５】
＜駆動部の構成＞
正の強度変調光を発生させる駆動電流を「第１の駆動電流」と呼び、負の強度変調光を発生させる駆動電流を「第２の駆動電流」と呼ぶ。
【０１５６】
第１の駆動部１１３ａあるいは第２の駆動部１１３ｂは、第１の駆動電流と第２の駆動電流の双方を発生可能な１つの電源から構成されてもよい。また、第１の駆動部１１３ａあるいは第２の駆動部１１３ｂは、は、第１の駆動電流を発生可能な第１の電源と、第２の駆動電流を発生可能な第２の電源とで構成されてもよい。ここでは、各駆動電流を別々の電源で生成する例を、図１８を用いて説明する。
【０１５７】
図１８に示す第１の駆動部１１３ａは、第１の駆動電流を発生可能な第１の電源２１０ａ、及び、第２の駆動電流を発生可能な第２の電源２２０ａを含む。制御部１５３は、１と０からなる第１の制御信号２３０と−１と０からなる第２の制御信号２４０を第１の駆動部１１３ａに送信する機能を有する。
【０１５８】
例えば、上記した第１の符号系列｛ａｉ｝＝｛１，１，−１，１，−１，−１，−１，１｝に対応する光照射を行う場合、第１の制御信号｛１，１，０，１，０，０，０，１｝と、第２の制御信号｛０，０，−１，０，−１，−１，−１，０｝に分離し、それぞれを第１の駆動部１１３ａに送信する。すなわち、制御部１５３は、第１の制御信号２３０を第１の電源２１０ａに送信し、第２の制御信号２４０を第２の電源２２０ａに送信する。
【０１５９】
第１の電源２１０ａは、第１の制御信号の符号要素｛１｝のタイミングに合わせて第１の駆動電流を生成し、第１の制御信号の符号要素｛０｝のタイミングでは電流を０とする、もしくは光音響波の発生が抑制された電流を生成する。第２の電源２２０ａは、第２の制御信号の符号要素｛−１｝のタイミングに合わせて第２の駆動電流を生成し、第２の制御信号の符号要素｛０｝のタイミングでは電流を０とする、もしくは光音響波の発生が抑制された電流を生成する。その結果、第１の光源１１１ａには、第１の符号系列｛ａｉ｝に対応する駆動電流（図９（ａ））と同様の電流が投入される。
【０１６０】
図１８に示す第２の駆動部１１３ｂは、第１の駆動電流を発生可能な第３の電源２１０ｂ、及び、第２の駆動電流を発生可能な第４の電源２２０ｂを含む。制御部１５３は、１と０からなる第３の制御信号２５０と−１と０からなる第４の制御信号２６０を第２の駆動部１１３ｂに送信する機能を有する。
【０１６１】
例えば、上記した第２の符号系列｛ｂｉ｝＝｛１，１，−１，１，１，１，１，−１｝に対応する光照射を行う場合、第３の制御信号｛１，１，０，１，１，１，１，０｝と、第４の制御信号｛０，０，−１，０，０，０，０，−１｝に分離し、それぞれを第２の駆動部１１３ｂに送信する。すなわち、制御部１５３は、第３の制御信号２５０を第３の電源２１０ｂに送信し、第４の制御信号２６０を第４の電源２２０ｂに送信する。
【０１６２】
第３の電源２１０ｂは、第３の制御信号の符号要素｛１｝のタイミングに合わせて第１の駆動電流を生成し、第３の制御信号の符号要素｛０｝のタイミングでは電流を０とする、もしくは光音響波の発生が抑制された電流を生成する。第４の電源２２０は、第４の制御信号の符号要素｛−１｝のタイミングに合わせて第２の駆動電流を生成し、第４の制御信号の符号要素｛０｝のタイミングでは電流を０とする、もしくは光音響波の発生が抑制された電流を生成する。その結果、第２の光源１１１ｂには、第３の符号系列｛ｂｉ｝に対応する駆動電流（図９（ｃ））と同様の電流が投入される。
【０１６３】
異なる駆動電流をひとつの電源で発生することに比べ、駆動電流毎に別々の電源を用いる装置の方が、第１の駆動部１１３ａあるいは第２の駆動部１１３ｂの設計を簡略化できる。また、駆動電流毎に別々の電源を用いる場合、異なる駆動電流を高速に切り替えるときの応答性が高い。その結果、異なる符号要素の光を時間的に重複して被検体に照射することができる。これにより、光の照射効率を高めることができ、短時間で高いＳＮ比の復号信号を取得することができる。
【０１６４】
本実施例において、第１の光源１１１ａと第２の光源１１１ｂのピーク光出力の最大強度を同一としたがそれに限ったものではない。本実施例において、第１の光源１１１ａにおける｛１｝のレベルと｛−１｝のレベルは実質的に揃っていることが好ましい。また、第２の光源１１１ｂにおける｛１｝のレベルと｛−１｝のレベルも実質的に揃っていることが好ましい。ここで、実質的に揃っているとは、平均化することでそのばらつきは無視できる程度に揃っている、という意味である。しかしながら、第１の光源１１１ａにおける｛１｝のレベルと第２の光源１１１ｂにおける｛１｝のレベルは一致していなくてよい。例えば、光源の個体差により、同じ投入電流でもその時の光出力が異なる場合がある。その場合、光出力の最大強度が一致するように投入電流を光源ごとにかえてもよい。また、ピーク光出力の最大強度が異なる場合でも、復号された受信信号をそれぞれの光出力の最大ピーク強度で規格化することで補正できる。具体的には、式３で得られたＤＳ１を第１の光源の最大ピーク強度で割り、式４で得られたＤＳ２を第２の光源の最大ピーク強度で割ることで規格化できる。あるいは、受信信号をそれぞれの光出力の最大強度で規格化したのちに復号処理を施してもよい。具体的には、受信信号Ｓ１及びＳ２を第１の光源の最大ピーク強度で割ったものを用いて式３の復号処理を行い、受信信号Ｓ１及びＳ２を第２の光源の最大ピーク強度で割ったものを用いて式４の復号処理を行えばよい。
【０１６５】
また、本実施例において、符号長、基準タイミングの時間間隔はここに示したものに限らず、被検体中の観察領域の深さや光源駆動部の性能などに合わせて、ＳＮ比の向上が図れるよう適切なものを用いればよい。
【実施例２】
【０１６６】
＜直交符号を適用した符号化・復号化＞
「完全直交の関係」を満たす２組の相補符号は存在するが、互いに「完全直交の関係」を満たす３組以上の相補符号は存在しない。よって、実施例１で述べた手法を３つ以上の光量変化量に適用することはできない。
【０１６７】
３つ以上の光量変化量を用いる場合の符号化・復号化の手法について以下に説明する。
【０１６８】
以下に示す、符号長４の互いに直交する４つの符号系列｛ａｉｋ｝（ｋ＝１から４、ｉ＝１から４）を考える。
第１の符号系列｛ａｉ１｝＝｛１，−１，−１，１｝
第２の符号系列｛ａｉ２｝＝｛１，−１，１，−１｝
第３の符号系列｛ａｉ３｝＝｛１，１，−１，−１｝
第４の符号系列｛ａｉ４｝＝｛１，１，１，１｝
【０１６９】
さらに、以下に示す、この４つの符号系列を重複なく順次並べるための順番を表す互いに異なる４つの順列｛ｇｍ｝（ｍ＝１から４、各要素は１から４のいずれか）を考える。
第１の順列｛ｇ１｝＝｛｛ａｉ１｝，｛ａｉ２｝，｛ａｉ３｝，｛ａｉ４｝｝
第２の順列｛ｇ２｝＝｛｛ａｉ２｝，｛ａｉ１｝，｛ａｉ４｝，｛ａｉ３｝｝
第３の順列｛ｇ３｝＝｛｛ａｉ３｝，｛ａｉ４｝，｛ａｉ１｝，｛ａｉ２｝｝
第４の順列｛ｇ４｝＝｛｛ａｉ４｝，｛ａｉ３｝，｛ａｉ２｝，｛ａｉ１｝｝
【０１７０】
ここで、順列ｇｐおよびｇｑに従って４つの符号系列を並べた場合の、それぞれの相互相関関数の総和を考える。式で表すと以下のように表現できる。
【０１７１】
【数５】
2019068914.tif 000007
【０１７２】
（ｐ，ｑは１から４の整数）
ここで、ｐ＝ｑのとき、（式６）は自己相関関数の総和となり、ピークにおいて１６となり、ピーク以外のすべての点で０となる。すなわち、
（ａ１＊  ａ１）＝｛１，−２，−１，４，−１，−２，１｝
（ａ２＊  ａ２）＝｛−１，２，−３，４，−３，２，−１｝
（ａ３＊  ａ３）＝｛−１，−２，１，４，１，−２，−１｝
（ａ４＊  ａ４）＝｛１，２，３，４，３，２，１｝
であり、任意のｐにおいて、
【０１７３】
【数６】
2019068914.tif 000008
が成立する。
【０１７４】
また、ｐ≠ｑのとき、（式６）は相互相関関数の総和となり、ｐ≠ｑであるすべての組み合わせにおいて、すべての点で０となる。例えば、ｐ＝１、ｑ＝２の場合、
（ａ１＊  ａ２）＝｛−１，２，−１，０，１，−２，１｝
（ａ２＊  ａ１）＝｛１，−２，１，０，−１，２，−１｝
（ａ３＊  ａ４）＝｛１，２，１，０，−１，−２，−１｝
（ａ４＊  ａ３）＝｛−１，−２，−１，０，１，２，１｝
であり、
【０１７５】
【数７】
2019068914.tif 000009
が成立する。
【０１７６】
このような特性を有する符号系列の組を用いると、以下のようなことが実現できる。すなわち、
・ある順列ｇｐで決まる順番の４つの符号系列に従って、符号化した信号を順次取得し、取得した信号を同じ順列ｇｐで決まる順番の４つの符号系列に従って順次復号化した場合、復号化した信号の総和はデルタ関数となる。
・ある順列ｇｐで決まる順番の４つの符号系列に従って、符号化した信号を順次取得し、取得した信号をそれと異なる順列ｇｑで決まる順番の４つの符号系列に従って順次復号化した場合、復号化した信号の総和は０となりキャンセルできる。
【０１７７】
このような符号系列とその順番を示す順列を互いに異なる複数の光量変化量を有する光を照射する光音響装置に適用した場合、光の照射時間が重複していても、ある光量変化量に対応する信号と、その他の光量変化量に対応する信号とを分離して取得することができる。
【０１７８】
符号化に用いる光量変化量の数、符号長、及び互いに直交する符号系列の数の好適な関係について説明する。符号化に用いる光量変化量の数が２から４の場合、符号長は４、直交する符号の数は４とするのが好ましい。符号化に用いる光量変化量の数が５から８の場合、符号長は８、直交する符号の数は８とするのが好ましい。直交する符号系列の数は、符号化に用いる光量変化量の数以上の２のべき乗数とするのが好ましい。
【０１７９】
本実施例では、図１９に示す光音響装置を用いて、互いに異なる３つの光量変化量を利用した符号化を行う場合を説明する。図１９は、光音響装置全体の概略ブロック図である。図３と共通する部材には同じ番号を付し、説明は省略する。本実施形態に係る光音響装置は、光照射部３１０、受信部１２０、信号収集部１４０、コンピュータ１５０、表示部１６０、及び入力部１７０を有する。
【０１８０】
光照射部３１０は、第１の光源３１１ａおよび第１の光源３１１ａから射出された光を被検体１００（図中不支持）に導く第１の光学系３１２ａ、および、第１の光源３１１ａの駆動を制御する第１の駆動部３１３ａを含む。また、光照射部３１０は、第２の光源３１１ｂおよび第２の光源３１１ｂから射出された光を被検体１００（図中不支持）に導く第２の光学系３１２ｂ、および、第２の光源３１１ｂの駆動を制御する第２の駆動部３１３ｂを含む。また、光照射部３１０は、第３の光源３１１ｃおよび第３の光源３１１ｃから射出された光を被検体１００（図中不支持）に導く第３の光学系３１２ｃ、および、第３の光源３１１ｃの駆動を制御する第３の駆動部３１３ｃを含む。
【０１８１】
本実施例では、第１の光源３１１ａ、第２の光源３１１ｂ、第３の光源３１１ｃとして、波長８０８ｎｍ、光出力最大５０Ｗの半導体レーザを用いる。また、第１の光源３１１ａは第１の駆動部３１３ａにより、基準タイミングにおける光量変化量が１００ｎｓ（第１の光量変化量）となるように駆動される。また、第２の光源３１１ｂは第２の駆動部３１３ｂにより、基準タイミングにおける光量変化量が２００ｎｓ（第２の光量変化量）となるように駆動される。また、第３の光源３１１ｃは第３の駆動部３１３ｃにより、基準タイミングにおける光量変化量が３００ｎｓ（第３の光量変化量）となるように駆動される。
【０１８２】
受信部１２０には、中心周波数４ＭＨｚ，６ｄＢ帯域１．６−６．４ＭＨｚの周波数特性を持つ圧電素子からなるリニアアレイを用いる。受信部１２０と被検体１００との間は音響マッチングのために超音波ジェルで埋められている。
【０１８３】
本実施例では、符号長４の互いに直交する４つの符号系列を用いる。具体的には、
第１の符号系列｛ａｉ１｝＝｛１，−１，−１，１｝
第２の符号系列｛ａｉ２｝＝｛１，−１，１，−１｝
第３の符号系列｛ａｉ３｝＝｛１，１，−１，−１｝
第４の符号系列｛ａｉ４｝＝｛１，１，１，１｝
である。また、第１、第２、第３の光源からの光それぞれに割り当てた順列を
第１の順列｛ｇ１｝＝｛｛ａｉ１｝，｛ａｉ２｝，｛ａｉ３｝，｛ａｉ４｝｝
第２の順列｛ｇ２｝＝｛｛ａｉ２｝，｛ａｉ１｝，｛ａｉ４｝，｛ａｉ３｝｝
第３の順列｛ｇ３｝＝｛｛ａｉ３｝，｛ａｉ４｝，｛ａｉ１｝，｛ａｉ２｝｝
とする。
【０１８４】
すなわち、第１の順列で決まる順番で４つの符号系列を第１の光源からの光に割り当てる。また、第２の順列で決まる順番で４つの符号系列を第１の光源からの光に割り当てる。また、第３の順列で決まる順番で４つの符号系列を第３の第１の光源からの光に割り当てる。
【０１８５】
本実施例に係る光音響装置の符号化・復号化のフローを説明する。
【０１８６】
本実施例では、図２０に示すシーケンスで、各順列の同じ番号の要素の強度変調光を同期して被検体１００に照射し、符号化が行われる。
【０１８７】
まず、１番目の順列要素については、光照射部１１０が、第１の符号系列｛ａｉ１｝に対応する第１の光量変化量を有する強度変調光、第２の符号系列｛ａｉ２｝に対応する第２の光量変化量を有する強度変調光、及び第３の符号系列｛ａｉ３｝に対応する第３の光量変化量を有する強度変調光を同期して照射する。そして、受信部１２０は、その光照射により発生した光音響波を受信し、受信信号Ｓ１を出力する。
【０１８８】
次に、２番目の順列要素については、光照射部１１０が、第２の符号系列｛ａｉ２｝に対応する第１の光量変化量を有する強度変調光、第１の符号系列｛ａｉ１｝に対応する第２の光量変化量を有する強度変調光、及び第４の符号系列｛ａｉ４｝に対応する第３の光量変化量を有する強度変調光を同期して照射する。そして、受信部１２０は、その光照射により発生した光音響波を受信し、受信信号Ｓ２を出力する。
【０１８９】
次に、３番目の順列要素については、光照射部１１０が、第３の符号系列｛ａｉ３｝に対応する第１の光量変化量を有する強度変調光、第４の符号系列｛ａｉ４｝に対応する第２の光量変化量を有する強度変調光、及び第１の符号系列｛ａｉ１｝に対応する第３の光量変化量を有する強度変調光を同期して照射する。そして、受信部１２０は、その光照射により発生した光音響波を受信し、受信信号Ｓ３を出力する。
【０１９０】
次に、４番目の順列要素については、光照射部１１０が、第４の符号系列｛ａｉ１｝に対応する第１の光量変化量を有する強度変調光、第３の符号系列｛ａｉ３｝に対応する第２の光量変化量を有する強度変調光、及び第２の符号系列｛ａｉ２｝に対応する第３の光量変化量を有する強度変調光を同期して照射する。そして、受信部１２０は、その光照射により発生した光音響波を受信し、受信信号Ｓ４を出力する。
【０１９１】
なお、各順列要素の強度変調光の光照射を完全に同期させずに行ってもよい。ただし、単位時間当たりの取得信号のＳＮ比向上のために、複数の光量変化量を有する強度変調光の照射期間が互いに少なくとも一部重なっていることが好ましい。
【０１９２】
＜Ｓ１．１番目の順列要素に対応する強度変調光照射＞
制御部１５３は、割り当てられた順列に従い、第１の符号系列｛ａｉ１｝に関する情報を第１の駆動部３１３ａに送信する。また、制御部１５３は、割り当てられた順列に従い、第２の符号系列｛ａｉ２｝に関する情報を第２の駆動部３１３ｂに送信する。また、制御部１５３は、割り当てられた順列に従い、第３の符号系列｛ａｉ３｝に関する情報を第３の駆動部３１３ｃに送信する。
【０１９３】
第１の光源３１１ａは、第１の符号系列｛ａｉ１｝に関する情報に基づいて第１の駆動部３１３ａが生成した駆動電流によって駆動される。発生した光は、第１の光学系３１２ａを通して、被検体１００に照射される。
【０１９４】
第２の光源３１１ｂは、第２の符号系列｛ａｉ２｝に関する情報に基づいて第２の駆動部３１３ｂが生成した駆動電流によって駆動される。発生した光は、第２の光学系３１２ｂを通して、被検体１００に照射される。
【０１９５】
第３の光源３１１ｃは、第３の符号系列｛ａｉ３｝に関する情報に基づいて第３の駆動部３１３ｃが生成した駆動電流によって駆動される。発生した光は、第３の光学系３１２ｃを通して、被検体１００に照射される。
【０１９６】
＜Ｓ２．１番目の順列要素に対応する受信信号の取得＞
第１の光源３１１ａ、第２の光源３１１ｂ、第３の光源３１１ｃからの光を被検体１００に照射することにより発生した光音響波を受信部１２０が受信することで受信信号Ｓ１（ｔ）が得られる。すなわち、第１の符号系列｛ａｉ１｝に対応する第１の光量変化量を有する強度変調光、第２の符号系列｛ａｉ２｝に対応する第２の光量変化量を有する強度変調光、及び第３の符号系列｛ａｉ３｝に対応する第３の光量変化量を有する強度変調光を略同じタイミングで被検体１００に照射することによりＳ１（ｔ）が得られる。
【０１９７】
＜Ｓ３．２番目の順列要素に対応する強度変調光照射＞
制御部１５３は、割り当てられた順列に従い、第２の符号系列｛ａｉ２｝に関する情報を第１の駆動部３１３ａに送信する。また、制御部１５３は、割り当てられた順列に従い、第１の符号系列｛ａｉ１｝に関する情報を第２の駆動部３１３ｂに送信する。また、制御部１５３は、割り当てられた順列に従い、第４の符号系列｛ａｉ４｝に関する情報を第３の駆動部３１３ｃに送信する。
【０１９８】
第１の光源３１１ａは、第２の符号系列｛ａｉ２｝に関する情報に基づいて第１の駆動部３１３ａが生成した駆動電流によって駆動される。発生した光は、第１の光学系３１２ａを通して、被検体１００に照射される。
【０１９９】
第２の光源３１１ｂは、第１の符号系列｛ａｉ１｝に関する情報に基づいて第２の駆動部３１３ｂが生成した駆動電流によって駆動される。発生した光は、第２の光学系３１２ｂを通して、被検体１００に照射される。
【０２００】
第３の光源３１１ｃは、第４の符号系列｛ａｉ４｝に関する情報に基づいて第３の駆動部３１３ｃが生成した駆動電流によって駆動される。発生した光は、第３の光学系３１２ｃを通して、被検体１００に照射される。
【０２０１】
＜Ｓ４．２番目の順列要素に対応する受信信号の取得＞
第１の光源３１１ａ、第２の光源３１１ｂ、第３の光源３１１ｃからの光を被検体１００に照射することにより発生した光音響波を受信部１２０が受信することで受信信号Ｓ２（ｔ）が得られる。すなわち、第２の符号系列｛ａｉ２｝に対応する第１の光量変化量を有する強度変調光、第１の符号系列｛ａｉ１｝に対応する第２の光量変化量を有する強度変調光、及び第４の符号系列｛ａｉ４｝に対応する第３の光量変化量を有する強度変調光を略同じタイミングで被検体１００に照射することによりＳ２（ｔ）が得られる。
【０２０２】
＜Ｓ５．３番目の順列要素に対応する強度変調光照射＞
制御部１５３は、割り当てられた順列に従い、第３の符号系列｛ａｉ３｝に関する情報を第１の駆動部３１３ａに送信する。また、制御部１５３は、割り当てられた順列に従い、第４の符号系列｛ａｉ４｝に関する情報を第２の駆動部３１３ｂに送信する。また、制御部１５３は、割り当てられた順列に従い、第１の符号系列｛ａｉ１｝に関する情報を第３の駆動部３１３ｃに送信する。
【０２０３】
第１の光源３１１ａは、第３の符号系列｛ａｉ３｝に関する情報に基づいて第１の駆動部３１３ａが生成した駆動電流によって駆動される。発生した光は、第１の光学系３１２ａを通して、被検体１００に照射される。
【０２０４】
第２の光源３１１ｂは、第４の符号系列｛ａｉ４｝に関する情報に基づいて第２の駆動部３１３ｂが生成した駆動電流によって駆動される。発生した光は、第２の光学系３１２ｂを通して、被検体１００に照射される。
【０２０５】
第３の光源３１１ｃは、第１の符号系列｛ａｉ１｝に関する情報に基づいて第３の駆動部３１３ｃが生成した駆動電流によって駆動される。発生した光は、第３の光学系３１２ｃを通して、被検体１００に照射される。
【０２０６】
＜Ｓ６．３番目の順列要素に対応する受信信号の取得＞
第１の光源３１１ａ、第２の光源３１１ｂ、第３の光源３１１ｃからの光を被検体１００に照射することにより発生した光音響波を受信部１２０が受信することで受信信号Ｓ３（ｔ）が得られる。すなわち、第３の符号系列｛ａｉ３｝に対応する第１の光量変化量を有する強度変調光、第４の符号系列｛ａｉ４｝に対応する第２の光量変化量を有する強度変調光、及び第１の符号系列｛ａｉ１｝に対応する第３の光量変化量を有する強度変調光を略同じタイミングに照射することによりＳ３（ｔ）が得られる。
【０２０７】
＜Ｓ７．４番目の順列要素に対応する強度変調光照射＞
制御部１５３は、割り当てられた順列に従い、第４の符号系列｛ａｉ４｝に関する情報を第１の駆動部３１３ａに送信する。また、制御部１５３は、割り当てられた順列に従い、第３の符号系列｛ａｉ３｝に関する情報を第２の駆動部３１３ｂに送信する。また、制御部１５３は、割り当てられた順列に従い、第２の符号系列｛ａｉ２｝に関する情報を第３の駆動部３１３ｃに送信する。
【０２０８】
第１の光源３１１ａは、第４の符号系列｛ａｉ４｝に関する情報に基づいて第１の駆動部３１３ａが生成した駆動電流によって駆動される。発生した光は、第１の光学系３１２ａを通して、被検体１００に照射される。
【０２０９】
第２の光源３１１ｂは、第３の符号系列｛ａｉ３｝に関する情報に基づいて第２の駆動部３１３ｂが生成した駆動電流によって駆動される。発生した光は、第２の光学系３１２ｂを通して、被検体１００に照射される。
【０２１０】
第３の光源３１１ｃは、第２の符号系列｛ａｉ２｝に関する情報に基づいて第３の駆動部３１３ｃが生成した駆動電流によって駆動される。発生した光は、第３の光学系３１２ｃを通して、被検体１００に照射される。
【０２１１】
＜Ｓ８．４番目の順列要素に対応する受信信号の取得＞
第１の光源３１１ａ、第２の光源３１１ｂ、第３の光源３１１ｃからの光を被検体１００に照射することにより発生した光音響波を受信部１２０が受信することで受信信号Ｓ４（ｔ）が得られる。すなわち、第４の符号系列｛ａｉ１｝に対応する第１の光量変化量を有する強度変調光、第３の符号系列｛ａｉ３｝に対応する第２の光量変化量を有する強度変調光、及び第２の符号系列｛ａｉ２｝に対応する第３の光量変化量を有する強度変調光を略同じタイミングに照射することによりＳ４（ｔ）が得られる。
【０２１２】
＜Ｓ９．復号化＞
コンピュータ１５０内の演算部１５１が行う、符号化された受信信号の復号方法について説明する。
【０２１３】
基準タイミングの時間間隔をΔｔとしたとき、演算部１５１は、（式９）から（式１１）にしたがった復号処理を行うことにより、第１、第２、第３の光量変化量のそれぞれに対応する復号信号ＤＳ１（ｔ）、ＤＳ２（ｔ）、ＤＳ３（ｔ）を得る。ここで、第１、第２、第３の光源それぞれに割り当てたのと同じ順列で決まる順番で、４つの符号系列を復号処理に用いている。
【０２１４】
【数８】
2019068914.tif 000010
【０２１５】
符号長をＮ、互いに直交する符号系列の数（すなわち順列の要素数）をＫとして（式９）から（式１１）を一般化すると、（式１２）のようになる。
【０２１６】
【数９】
2019068914.tif 000011
【０２１７】
ここで、ｉは１以上の自然数、｛ｇｍ（ｊ）｝は複数の光源それぞれに割り当てられた順列、Ｓｊは順列要素に対応する受信信号、ｊは１以上の自然数、ＫはＫ≧Ｍを満たす２のべき乗数である。また、ｍは１以上Ｍ以下の自然数、Ｍは複数の光源の数、ｔは時間、Δｔは符号系列の符号要素の基準タイミングの時間間隔、である。
【０２１８】
（式９）に示された復号処理を行うことにより得られた復号信号ＤＳ１は、第１の光量変化量を有する光（立ち上がりまたは立ち下がりの時間が１００ｎｓの光）を用いて光音響波を発生させた場合に得られる受信信号に相当する。
【０２１９】
（式１０）に示された復号処理を行うことにより得られた復号信号ＤＳ２は、第２の光量変化量を有する光（立ち上がりまたは立ち下がりの時間が２００ｎｓの光）を用いて光音響波を発生させた場合に得られる受信信号に相当する。
【０２２０】
（式１１）に示された復号処理を行うことにより得られた復号信号ＤＳ３は、第３の光量変化量を有する光（立ち上がりまたは立ち下がりの時間が３００ｎｓの光）を用いて光音響波を発生させた場合に得られる受信信号に相当する。
【０２２１】
復号信号ＤＳ１、ＤＳ２、及びＤＳ３を解析することで、光量変化量ごとに光音響画像を複数生成することができる。それをディスプレイなどに表示する際は、それぞれ独立に表示したり、重畳して表示したりすることができる。さらには、暗い部分にはより強い重みを付けて重畳する、などを行ってもよい。演算部１５１は、実施例１と同様に、各光量変化量に対応する復号信号に基づいて、各種の光音響画像を生成することができる。
【０２２２】
なお、実施例１と同様に、制御部１５３は、観察領域（関心領域）や被検体の音速などに応じて、短時間で高ＳＮ比の信号を取得することができるように、基準タイミングの時間間隔を設定してもよい。また、再構成手法に関しても、実施例１と同様の手法を採用してもよい。また、本実施例においても、実施例１で説明した駆動部の構成を採用してもよい。
【実施例３】
【０２２３】
＜少ない符号系列数での符号化・復号化＞
実施例２で説明した手法において、例えば、符号長を４ではなく８にすることで、よりＳＮ比を向上させることが可能である。しかしながら、その場合、互いに直交する符号系列の数も８となり、符号化するための光照射と光音響波の受信についても８回繰り返す必要があり、測定時間が増加する。
【０２２４】
本実施例では、符号系列の数を少なくすることにより、測定時間の増大を抑えつつ、ＳＮ比を向上させる例を説明する。
【０２２５】
以下に示す、符号長８の互いに直交する４つの符号系列｛ａｉｋ｝（ｋ＝１から４、ｉ＝１から８）を考える。
第１の符号系列｛ａｉ１｝＝｛１，１，−１，−１，−１，−１，１，１｝
第２の符号系列｛ａｉ２｝＝｛１，１，−１−１，１，１，−１，−１｝
第３の符号系列｛ａｉ３｝＝｛１，１，１，１，−１，−１，−１，−１｝
第４の符号系列｛ａｉ４｝＝｛１，１，１，１，１，１，１，１｝
【０２２６】
これは、実施例２で説明した各符号系列の｛１｝に｛１，１｝を割り当て、｛−１｝に｛−１，−１｝を割り当てたものである。すなわち、互いに直交する４つの符号系列の各符号要素を２回ずつ繰り返すことにより、符号要素の数を２倍としている。なお、各符号要素の繰り返し回数は、２以上の自然数であればいかなる値を採用してもよい。
【０２２７】
さらに、以下に示す、この４つの符号系列を重複なく順次並べるための順列を表す互いに異なる４つの順列を考える。
第１の順列｛ｇ１｝＝｛｛ａｉ１｝，｛ａｉ２｝，｛ａｉ３｝，｛ａｉ４｝｝
第２の順列｛ｇ２｝＝｛｛ａｉ２｝，｛ａｉ１｝，｛ａｉ４｝，｛ａｉ３｝｝
第３の順列｛ｇ３｝＝｛｛ａｉ３｝，｛ａｉ４｝，｛ａｉ１｝，｛ａｉ２｝｝
第４の順列｛ｇ４｝＝｛｛ａｉ４｝，｛ａｉ３｝，｛ａｉ２｝，｛ａｉ１｝｝
【０２２８】
ここで、（式６）で示したように、順列ｇｐおよびｇｑに従って４つの符号系列を並べた場合の、それぞれの相互相関関数の総和を考える。
【０２２９】
ｐ＝ｑのとき、（式６）は自己相関関数の総和となり、ピーク近傍において１６、３２、１６、それ以外のすべての点でＯとなる。すなわち、
（ａ１＊  ａ１）＝｛１，−２，−１，−４，−３，−２，３，８，３，−２，−３，−４，−１，２，１｝
（ａ２＊  ａ２）＝｛−１，−２，１，４，−１，−６，１，８，１，−６，−１，４，１，−２，−１｝
（ａ３＊  ａ３）＝｛−１，−２，−３，−４，−１，２，５，８，５，２，−１，−４，−３，−２，−１｝
（ａ４＊  ａ４）＝｛１，２，３，４，５，６，７，８，７，６，５，４，３，２，１｝
であり、任意のｐにおいて、
【０２３０】
【数１０】
2019068914.tif 000012
が成立する。
【０２３１】
また、ｐ≠ｑのとき、（式６）は相互相関関数の総和となり、ｐ≠ｑであるすべての組み合わせにおいて、すべての点でＯとなる。例えば、ｐ＝１、ｑ＝２の場合、
（ａ１＊  ａ２）＝｛−１，−２，１，４，１，−２，−１，０，１，２，−１，−４，−１，２，１｝
（ａ２＊  ａ１）＝｛１，２，−１，−４，−１，２，１，０，−１，−２，１，４，１，−２，−１｝
（ａ３＊  ａ４）＝｛１，２，３，４，３，２，１，０，−１，−２，−３，−４，−３，−２，−１｝
（ａ４＊  ａ３）＝｛−１，−２，−３，−４，−３，−２，−１，０，１，２，３，４，３，２，１｝
であり、
【０２３２】
【数１１】
2019068914.tif 000013
が成立する。
【０２３３】
このような符号系列とその順番を示す順列を互いに異なる複数の光量変化量を有する複数の光を照射する光音響装置に適用した場合、光の照射時間が重複していても、ある光量変化量に対応する信号と、その他の光量変化量に対応する信号とを分離して取得することができる。
【０２３４】
実施例２と比べると、Δｔ＜＜Ｔｔｏｆの場合、測定時間をほとんど増加させずに、各光量変化量に対応する復号信号のピーク強度を増加させることが可能となる。その際、サイドローブも増大するが、そのパターンは既知であるので必要に応じ、復号信号に所望の補正をかければよい。例えば、復号信号のうち、あるしきい値以上の信号のみを有効としたりしてもよい。また、復号した受信信号はΔｔの間隔で１：２：１の信号比を有することがわかっているので、これを０：１：０の信号に補正するデコンボリューションフィルタを重畳してもよい。演算部１５１は、実施例１と同様に、各光量変化量に対応する復号信号に基づいて、種々の光音響画像を生成することができる。
【０２３５】
（その他の実施例）
また、本発明は、以下の処理を実行することによっても実現される。即ち、上述した実施形態の機能を実現するソフトウェア（プログラム）を、ネットワーク又は各種記憶媒体を介してシステム或いは装置に供給し、そのシステム或いは装置のコンピュータ（またはＣＰＵやＭＰＵ等）がプログラムを読み出して実行する処理である。
【符号の説明】
【０２３６】
１１０  光照射部
１２０  受信部
１５０  コンピュータ
【図１】
2019068914.tif 000014
【図２】
2019068914.tif 000015
【図３】
2019068914.tif 000016
【図４】
2019068914.tif 000017
【図５】
2019068914.tif 000018
【図６】
2019068914.tif 000019
【図７Ａ】
2019068914.tif 000020
【図７Ｂ】
2019068914.tif 000021
【図８Ａ】
2019068914.tif 000022
【図８Ｂ】
2019068914.tif 000023
【図９】
2019068914.tif 000024
【図１０】
2019068914.tif 000025
【図１１】
2019068914.tif 000026
【図１２】
2019068914.tif 000027
【図１３】
2019068914.tif 000028
【図１４】
2019068914.tif 000029
【図１５】
2019068914.tif 000030
【図１６】
2019068914.tif 000031
【図１７】
2019068914.tif 000032
【図１８】
2019068914.tif 000033
【図１９】
2019068914.tif 000034
【図２０】
2019068914.tif 000035