(57)【要約】      （修正有）【課題】観察に適した動画像を生成することを目的の一つとする。【解決手段】画像処理装置１３００は、被検体への光照射により発生する光音響波に由来するボリュームデータである光音響画像データを取得する第１の取得手段と、前記光音響画像データを少なくとも２つ以上合成することにより合成画像データを取得する第２の取得手段と、前記合成画像データに基づいて、動画データの生成条件を決定する決定手段と、前記生成条件に基づいて、前記合成画像データから前記動画データを生成する生成手段と、を備えることを特徴とする。【選択図】図１
2020162746.tif 000002
【特許請求の範囲】
【請求項１】
被検体への光照射により発生する光音響波に由来するボリュームデータである光音響画像データを取得する第１の取得手段と、
前記光音響画像データを少なくとも２つ以上合成することにより合成画像データを取得する第２の取得手段と、
前記合成画像データに基づいて、動画データの生成条件を決定する決定手段と、
前記生成条件に基づいて、前記合成画像データから前記動画データを生成する生成手段と、
を備えることを特徴とする画像処理装置。
【請求項２】
前記第２の取得手段は、前記被検体とプローブの相対位置が異なる前記光音響画像データを少なくとも２つ以上合成することにより前記合成画像データを取得することを特徴とする請求項１に記載の画像処理装置。
【請求項３】
前記第２の取得手段は、前記光音響画像データから複数の前記合成画像データを生成することを特徴とする請求項１または２に記載の画像処理装置。
【請求項４】
前記第２の取得手段は、前記相対位置の変位量が互いに小さい前記光音響画像データを少なくとも２つ以上合成することにより前記合成画像データを取得することを特徴とする請求項１乃至３のいずれか１項に記載の画像処理装置。
【請求項５】
前記第２の取得手段は、観察対象となる領域の特徴量に基づいた数の前記光音響画像データを合成することにより前記合成画像データを取得することを特徴とする請求項１乃至４のいずれか１項に記載の画像処理装置。
【請求項６】
前記第２の取得手段は、前記観察対象となる血管が細くなるほど、合成する前記光音響画像データの数を多くすることを特徴とする請求項５に記載の画像処理装置。
【請求項７】
前記決定手段は、前記動画データを生成する画像領域、最大値投影方向、合成する画像データの数または前記画像データの階調のうち少なくともいずれか１つを前記生成条件として決定することを特徴とする請求項１乃至６のいずれか１項に記載の画像処理装置。
【請求項８】
前記決定手段は、閾値以上の画素値を有する画素を含まない領域を、前記動画データを生成する画像領域として決定することを特徴とする請求項１乃至７のいずれか１項に記載の画像処理装置。
【請求項９】
前記生成手段は、前記合成画像データに対するレンダリング処理により得られる２次元画像を用いて前記動画データを生成することを特徴とする請求項１乃至８のいずれか１項に記載の画像処理装置。
【請求項１０】
前記生成手段は、前記決定された最大値投影方向に基づいて前記合成画像データを最大値投影することにより得られる最大値投影画像データを用いて前記動画データを生成することを特徴とする請求項９に記載の画像処理装置。
【請求項１１】
前記生成手段は、観察対象となる血管の太さに基づいた数の前記光音響画像データを用いて前記動画データを生成することを特徴とする請求項１乃至１０のいずれか１項に記載の画像処理装置。
【請求項１２】
前記第２の取得手段は、前記生成手段が前記動画データの生成に用いる前記光音響画像データの数よりも多い前記光音響画像データを合成することにより前記合成画像データを生成することを特徴とする請求項１乃至１１のいずれか１項に記載の画像処理装置。
【請求項１３】
ユーザーからの入力を受け付ける受付手段をさらに備え、
前記決定手段は、前記入力に基づいて前記生成条件を決定することを特徴とする請求項１乃至１２のいずれか１項に記載の画像処理装置。
【請求項１４】
前記光音響画像データ、前記合成画像データまたは前記動画データのうち少なくとも１つを画像または動画として表示部に表示させる表示制御手段をさらに備えることを特徴とする請求項１乃至１３のいずれか１項に記載の画像処理装置。
【請求項１５】
前記表示制御手段は、前記合成画像データを参照画像として前記表示部に表示し、
前記決定手段は、前記参照画像をユーザーに提示して、前記動画データの生成条件について、該ユーザーからの入力を受けることにより前記動画条件を決定することを特徴とする請求項１４に記載の画像処理装置。
【請求項１６】
被検体への光照射により発生する光音響波に由来するボリュームデータである光音響画像データを取得する第１の取得手段と、
前記光音響画像データを少なくとも２つ以上合成することにより合成画像データを取得する第２の取得手段と、
前記合成画像データに基づいて、動画の生成条件を決定する決定手段と、
前記生成条件に基づいて、前記合成された前記光音響画像データとは異なる前記光音響画像データを用いて動画データを生成する生成手段と、
を備えることを特徴とする画像処理装置。
【請求項１７】
前記生成手段は、前記合成された前記光音響画像データとは異なる波長の光照射から得られた前記光音響画像データを用いて動画データを生成することを特徴とする請求項１６に記載の画像処理装置。
【請求項１８】
被検体への光照射により発生する光音響波に由来するボリュームデータである光音響画像データを取得する第１の取得工程と、
前記光音響画像データを少なくとも２つ以上合成することにより合成画像データを取得する第２の取得工程と、
前記合成画像データに基づいて、動画の生成条件を決定する決定工程と、
前記生成条件に基づいて、前記合成画像データから動画データを生成する生成工程と、
を備えることを特徴とする画像処理方法。
【請求項１９】
被検体への光照射により発生する光音響波に由来するボリュームデータである光音響画像データを取得する第１の取得工程と、
前記光音響画像データを少なくとも２つ以上合成することにより合成画像データを取得する第２の取得工程と、
前記合成画像データに基づいて、動画の生成条件を決定する決定工程と、
前記生成条件に基づいて、前記合成された前記光音響画像データとは異なる前記光音響画像データを用いて動画データを生成する生成工程と、
を備えることを特徴とする画像処理方法。
【請求項２０】
請求項１乃至１７のいずれか１項に記載の画像処理装置の各手段の制御をコンピュータに実行させることを特徴とするプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本明細書の開示は、画像処理装置、画像処理方法及びプログラムに関する。
【背景技術】
【０００２】
被検体に光を照射することに起因して生じた光音響波を検知することにより、被検体内の特性情報を取得する光音響装置が知られている。特許文献１には、光音響波の受信信号である光音響信号に基づいて、３次元のボリュームデータを生成することが記載されている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１４−６８７０１号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかしながら、時系列に取得した複数の３次元のボリュームデータを動画化する場合、生成される動画の観察断面方向や画像範囲、時間分解能が必ずしも観察に適したものではないという課題があった。
【０００５】
本明細書の開示は上記課題に鑑みてなされたものであり、観察に適した動画像を生成することを目的の一つとする。
【０００６】
なお、前記目的に限らず、後述する発明を実施するための形態に示す各構成により導かれる作用効果であって、従来の技術によっては得られない作用効果を奏することも本明細書の開示の他の目的の１つとして位置付けることができる。
【課題を解決するための手段】
【０００７】
本明細書に開示の画像処理装置は、被検体への光照射により発生する光音響波に由来するボリュームデータである光音響画像データを取得する第１の取得手段と、前記光音響画像データを少なくとも２つ以上合成することにより合成画像データを取得する第２の取得手段と、前記合成画像データに基づいて、動画データの生成条件を決定する決定手段と、前記生成条件に基づいて、前記合成画像データから前記動画データを生成する生成手段と、を備えることを特徴とする。
【発明の効果】
【０００８】
本明細書の開示によれば、観察に適した動画データを生成することができる。
【図面の簡単な説明】
【０００９】
【図１】本実施形態に係るシステムの構成の一例を示すブロック図。
【図２】本実施形態に係る画像処理装置とその周辺構成の一例を示すブロック図。
【図３】本実施形態に係る光音響装置の詳細な構成の一例を示すブロック図。
【図４】本実施形態に係るプローブの一例を示す模式図。
【図５】第１の実施形態に係る光音響装置が実行する画像処理方法の一例を示すフロー図。
【図６】第１の実施形態に係るデータの流れの一例を示す図。
【図７】第２の実施形態に係る光音響装置が実行する画像処理方法の一例を示すフロー図。
【発明を実施するための形態】
【００１０】
以下に図面を参照しつつ、本明細書に開示の画像処理装置の好適な実施の形態について説明する。ただし、以下に記載されている構成部品の寸法、材質、形状およびそれらの相対配置などは、発明が適用される装置の構成や各種条件により適宜変更されるべきものである。よって、この発明の範囲を以下の記載に限定する趣旨のものではない。
【００１１】
本明細書に開示の画像処理装置に係るシステムにより得られる光音響画像は、光エネルギーの吸収量や吸収率を反映している。光音響画像は、光音響波の発生音圧（初期音圧）、光吸収エネルギー密度、及び光吸収係数などの少なくとも１つの被検体情報の空間分布を表す画像である。光音響画像は、２次元の空間分布を表す画像であってもよいし、３次元の空間分布を表す画像（ボリュームデータ）であってもよい。なお、光音響画像は被検体表面から深さ方向の２次元の空間分布を表す画像または３次元の空間分布を表す画像であってもよい。
【００１２】
また、本明細書に開示の画像処理装置に係るシステムは、複数の波長に対応する複数の光音響画像を用いて被検体の機能画像を生成することができる。機能画像は、グルコース濃度、コラーゲン濃度、メラニン濃度、脂肪や水の体積分率など、被検体を構成する物質の濃度に相当する情報を示す画像であってもよい。機能画像は、第１波長λ１の光照射により発生した光音響波に基づいた吸収係数画像μａλ１（ｒ）および第２波長λ２の光照射により発生した光音響波に基づいた吸収係数画像μａλ２（ｒ）を用いて生成された酸素飽和度画像ＳＯ２（ｒ）であってもよい。例えば、本発明に係るシステムは、式（１）にしたがって、酸素飽和度画像ＳＯ２（ｒ）を機能画像として生成してもよい。
【００１３】
【数１】
2020162746.tif 000003
【００１４】
ここで、εＨｂλ１は第１波長λ１に対応するデオキシヘモグロビンのモラー吸収係数［ｍｍ−１ｍｏｌ−１］であり、εＨｂλ２は第２波長λ２に対応するデオキシヘモグロビンのモラー吸収係数［ｍｍ−１ｍｏｌ−１］である。εＨｂＯλ１は第１波長λ１に対応するオキシヘモグロビンのモラー吸収係数［ｍｍ−１ｍｏｌ−１］であり、εＨｂＯλ２は第２波長λ２に対応するオキシヘモグロビンのモラー吸収係数［ｍｍ−１ｍｏｌ−１］である。ｒは位置である。
【００１５】
また、本明細書に開示の画像処理装置に係るシステムは、第１波長λ１の光照射により発生した光音響波に基づいた第１光音響画像および第２波長λ２の光照射により発生した光音響波に基づいた第２光音響画像の比を示す画像を機能画像としてもよい。すなわち、第１波長λ１の光照射により発生した光音響波に基づいた第１光音響画像および第２波長λ２の光照射により発生した光音響波に基づいた第２光音響画像の比に基づいた画像を機能画像としてよい。なお、式（１）の変形式にしたがって生成される画像も、第１光音響画像および第２光音響画像の比によって表現できるため、第１光音響画像および第２光音響画像の比に基づいた画像（機能画像）といえる。
【００１６】
なお、機能画像は被検体表面から深さ方向の２次元の空間分布を表す画像または３次元の空間分布を表す画像であってもよい。
【００１７】
（第１の実施形態）
以下、本実施形態のシステムの構成及び画像処理方法について説明する。
【００１８】
図１を用いて本実施形態に係るシステムを説明する。図１は、本実施形態に係るシステムの構成の一例を示すブロック図である。本実施形態に係るシステムは、光音響装置１１００、記憶装置１２００、画像処理装置１３００、表示装置１４００、及び入力装置１５００を備える。装置間のデータの送受信は有線で行われてもよいし、無線で行われてもよい。
【００１９】
光音響装置１１００は、被検体１００を撮影することにより光音響画像を生成し、記憶装置１２００もしくは画像処理装置１３００に出力する。光音響装置１１００は、光照射により発生した光音響波を受信することにより得られる受信信号を用いて、被検体１００内の複数位置のそれぞれに対応する特性値の情報を取得する装置である。すなわち、光音響装置１１００は、被検体１００への光照射により、被検体内で発生した光音響波を受信して得られる光音響信号に基づいて、被検体の光学特性に関連した機能情報を可視化した画像データ（光音響画像）を取得する装置である。
【００２０】
記憶装置１２００は、ＲＯＭ（Ｒｅａｄ  ｏｎｌｙ  ｍｅｍｏｒｙ）、磁気ディスクやフラッシュメモリなどの記憶媒体であってもよい。また、記憶装置１２００は、ＰＡＣＳ（Ｐｉｃｔｕｒｅ  Ａｒｃｈｉｖｉｎｇ  ａｎｄ  Ｃｏｍｍｕｎｉｃａｔｉｏｎ  Ｓｙｓｔｅｍ）等のネットワークを介した記憶サーバであってもよい。すなわち、記憶装置１２００は記憶部の一例に相当する。
【００２１】
画像処理装置１３００は、記憶装置１２００に記憶された光音響画像や光音響画像の付帯情報等の情報を処理する装置である。
【００２２】
画像処理装置１３００の演算機能を担うユニットは、ＣＰＵやＧＰＵ（Ｇｒａｐｈｉｃｓ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）等のプロセッサ、ＦＰＧＡ（Ｆｉｅｌｄ  Ｐｒｏｇｒａｍｍａｂｌｅ  Ｇａｔｅ  Ａｒｒａｙ）チップ等の演算回路で構成されることができる。これらのユニットは、単一のプロセッサや演算回路から構成されるだけでなく、複数のプロセッサや演算回路から構成されていてもよい。
【００２３】
画像処理装置１３００の記憶機能を担うユニットは、ＲＯＭ（Ｒｅａｄ  ｏｎｌｙ  ｍｅｍｏｒｙ）、磁気ディスクやフラッシュメモリなどの非一時記憶媒体で構成することができる。また、記憶機能を担うユニットは、ＲＡＭ（Ｒａｎｄｏｍ  Ａｃｃｅｓｓ  Ｍｅｍｏｒｙ）などの揮発性の媒体であってもよい。なお、プログラムが格納される記憶媒体は、非一時記憶媒体である。なお、記憶機能を担うユニットは、１つの記憶媒体から構成されるだけでなく、複数の記憶媒体から構成されていてもよい。すなわち、画像処理装置１３００の記憶機能を担うユニットは記憶部の一例に相当する。
【００２４】
画像処理装置１３００の制御機能を担うユニットは、ＣＰＵなどの演算素子で構成される。制御機能を担うユニットは、システムの各構成の動作を制御する。制御機能を担うユニットは、入力部１７０からの測定開始などの各種操作による指示信号を受けて、システムの各構成を制御してもよい。また、制御機能を担うユニットは、記憶部に格納されたプログラムコードを読み出し、システムの各構成の作動を制御してもよい。
【００２５】
表示装置１４００は、液晶ディスプレイや有機ＥＬ（Ｅｌｅｃｔｒｏ  Ｌｕｍｉｎｅｓｃｅｎｃｅ）などのディスプレイである。また、表示装置１４００は、画像や装置を操作するためのＧＵＩを表示してもよい。
【００２６】
入力装置１５００としては、ユーザーが操作可能な、マウスやキーボードなどで構成される操作コンソールを採用することができる。また、表示装置１４００をタッチパネルで構成し、表示装置１４００を入力装置１５００として利用してもよい。
【００２７】
図２は、図１に示す本実施形態に係るシステムの具体的な構成例を示す。画像処理装置１３００は、ＣＰＵ１３１０、ＧＰＵ１３２０、ＲＡＭ１３３０、ＲＯＭ１３４０、外部記憶装置１３５０から構成される。画像処理装置１３００には、表示装置１４００としての液晶ディスプレイ１４１０が接続されている。また、入力装置１５００としてのマウス１５１０、キーボード１５２０が接続されている。さらに、画像処理装置１３００は、ＰＡＣＳ（Ｐｉｃｔｕｒｅ  Ａｒｃｈｉｖｉｎｇ  ａｎｄ  Ｃｏｍｍｕｎｉｃａｔｉｏｎ  Ｓｙｓｔｅｍ）などの記憶装置１２００としての画像サーバ１２１０と接続されている。これにより、画像データを画像サーバ１２１０上に保存したり、画像サーバ１２１０上の画像データを液晶ディスプレイ１４１０に表示したりすることができる。
【００２８】
次に、本実施形態に係るシステムに含まれる装置の構成例を説明する。
【００２９】
図３は、本実施形態に係るシステムに含まれる装置の概略ブロック図である。
【００３０】
本実施形態に係る光音響装置１１００は、駆動部１３０、信号収集部１４０、コンピュータ１５０、及びプローブ１８０を有する。プローブ１８０は、光照射部１１０、及び受信部１２０を有する。また、図４は、本実施形態に係るプローブ１８０の模式図を示す。測定対象は、被検体１００である。駆動部１３０は、光照射部１１０と受信部１２０を駆動し、機械的な走査を行う。光照射部１１０が光を被検体１００に照射し、被検体１００内で音響波が発生する。光に起因して光音響効果により発生する音響波を光音響波とも呼ぶ。受信部１２０は、光音響波を受信することによりアナログ信号としての電気信号（光音響信号）を出力する。
【００３１】
信号収集部１４０は、受信部１２０から出力されたアナログ信号をデジタル信号に変換し、コンピュータ１５０に出力する。コンピュータ１５０は、信号収集部１４０から出力されたデジタル信号を、光音響波に由来する信号データとして記憶する。
【００３２】
コンピュータ１５０は、記憶されたデジタル信号に対して信号処理を行うことにより、光音響画像を生成する。また、コンピュータ１５０は、得られた光音響画像に対して画像処理を施した後に、光音響画像を表示部１６０に出力する。表示部１６０は、光音響画像に基づいた画像を表示する。表示画像は、ユーザーやコンピュータ１５０からの保存指示に基づいて、コンピュータ１５０内のメモリや、モダリティとネットワークで接続されたデータ管理システムなどの記憶装置１２００に保存される。
【００３３】
また、コンピュータ１５０は、光音響装置に含まれる構成の駆動制御も行う。また、表示部１６０は、コンピュータ１５０で生成された画像の他にＧＵＩなどを表示してもよい。入力部１７０は、ユーザーが情報を入力できるように構成されている。ユーザーは、入力部１７０を用いて測定開始や終了、生成画像の保存指示などの操作を行うことができる。
【００３４】
以下、本実施形態に係る光音響装置１１００の各構成の詳細を説明する。
【００３５】
（光照射部１１０）
光照射部１１０は、光を発する光源１１１と、光源１１１から射出された光を被検体１００へ導く光学系１１２とを含む。なお、光は、いわゆる矩形波、三角波などのパルス光を含む。
【００３６】
光源１１１が発する光のパルス幅としては、１ｎｓ以上、１００ｎｓ以下のパルス幅であってもよい。また、光の波長として４００ｎｍから１６００ｎｍ程度の範囲の波長であってもよい。血管を高解像度でイメージングする場合は、血管での吸収が大きい波長（４００ｎｍ以上、７００ｎｍ以下）を用いてもよい。生体の深部をイメージングする場合には、生体の背景組織（水や脂肪など）において典型的に吸収が少ない波長（７００ｎｍ以上、１１００ｎｍ以下）の光を用いてもよい。
【００３７】
光源１１１としては、レーザーや発光ダイオードを用いることができる。また、複数波長の光を用いて測定する際には、波長の変更が可能な光源であってもよい。なお、複数波長を被検体１００に照射する場合、互いに異なる波長の光を発生する複数台の光源を用意し、それぞれの光源から交互に照射することも可能である。複数台の光源を用いた場合もそれらをまとめて光源として表現する。レーザーとしては、固体レーザー、ガスレーザー、色素レーザー、半導体レーザーなど様々なレーザーを使用することができる。例えば、Ｎｄ：ＹＡＧレーザーやアレキサンドライトレーザーなどのパルスレーザーを光源として用いてもよい。また、Ｎｄ：ＹＡＧレーザー光を励起光とするＴｉ：ｓａレーザーやＯＰＯ（Ｏｐｔｉｃａｌ  Ｐａｒａｍｅｔｒｉｃ  Ｏｓｃｉｌｌａｔｏｒｓ）レーザーを光源として用いてもよい。また、光源１１１としてフラッシュランプや発光ダイオードを用いてもよい。また、光源１１１としてマイクロウェーブ源を用いてもよい。
【００３８】
光学系１１２には、レンズ、ミラー、光ファイバ等の光学素子を用いることができる。乳房等を被検体１００とする場合、パルス光のビーム径を広げて照射するために、光学系の光出射部は光を拡散させる拡散板等で構成されていてもよい。一方、光音響顕微鏡においては、解像度を上げるために、光学系１１２の光出射部はレンズ等で構成し、ビームをフォーカスして照射してもよい。
【００３９】
なお、光照射部１１０が光学系１１２を備えずに、光源１１１から直接被検体１００に光を照射してもよい。
【００４０】
（受信部１２０）
受信部１２０は、音響波を受信することにより電気信号を出力するトランスデューサ１２１と、トランスデューサ１２１を支持する支持体１２２とを含む。また、トランスデューサ１２１は、音響波を送信する送信手段としてもよい。受信手段としてのトランスデューサと送信手段としてのトランスデューサとは、単一（共通）のトランスデューサでもよいし、別々の構成であってもよい。
【００４１】
トランスデューサ１２１を構成する部材としては、ＰＺＴ（チタン酸ジルコン酸鉛）に代表される圧電セラミック材料や、ＰＶＤＦ（ポリフッ化ビニリデン）に代表される高分子圧電膜材料などを用いることができる。また、圧電素子以外の素子を用いてもよい。例えば、静電容量型トランスデューサ（ＣＭＵＴ：Ｃａｐａｃｉｔｉｖｅ  Ｍｉｃｒｏ−ｍａｃｈｉｎｅｄ  Ｕｌｔｒａｓｏｎｉｃ  Ｔｒａｎｓｄｕｃｅｒｓ）、ファブリペロー干渉計を用いたトランスデューサなどを用いることができる。なお、音響波を受信することにより電気信号を出力できる限り、いかなるトランスデューサを採用してもよい。また、トランスデューサにより得られる信号は時間分解信号である。つまり、トランスデューサにより得られる信号の振幅は、各時刻にトランスデューサで受信される音圧に基づく値（例えば、音圧に比例した値）を表したものである。
【００４２】
光音響波を構成する周波数成分は、典型的には１００ＫＨｚから１００ＭＨｚであり、トランスデューサ１２１として、これらの周波数を検出することのできるものを採用してもよい。
【００４３】
支持体１２２は、機械的強度が高い金属材料などから構成されていてもよい。照射光を被検体１００に多く入射させるために、支持体１２２の被検体１００側の表面に、鏡面加工もしくは光散乱させる加工が行われていてもよい。本実施形態において支持体１２２は半球殻形状であり、半球殻上に複数のトランスデューサ１２１を支持できるように構成されている。この場合、支持体１２２に配置されたトランスデューサ１２１の指向軸は半球の曲率中心付近に集まる。そして、複数のトランスデューサ１２１から出力された信号を用いて画像化したときに曲率中心付近の画質が高くなる。なお、支持体１２２はトランスデューサ１２１を支持できる限り、いかなる構成であってもよい。支持体１２２は、１Ｄアレイ、１．５Ｄアレイ、１．７５Ｄアレイ、２Ｄアレイと呼ばれるような平面又は曲面内に、複数のトランスデューサを並べて配置してもよい。複数のトランスデューサ１２１が複数の受信手段に相当する。
【００４４】
また、支持体１２２は音響マッチング材を貯留する容器として機能してもよい。すなわち、支持体１２２をトランスデューサ１２１と被検体１００との間に音響マッチング材を配置するための容器としてもよい。
【００４５】
また、受信部１２０が、トランスデューサ１２１から出力される時系列のアナログ信号を増幅する増幅器を備えてもよい。また、受信部１２０が、トランスデューサ１２１から出力される時系列のアナログ信号を時系列のデジタル信号に変換するＡ／Ｄ変換器を備えてもよい。すなわち、受信部１２０が後述する信号収集部１４０を備えてもよい。
【００４６】
なお、音響波を様々な角度で検出できるようにするために、理想的には被検体１００を全周囲から囲むようにトランスデューサ１２１を配置してもよい。ただし、被検体１００が大きく全周囲を囲むようにトランスデューサを配置できない場合は、半球状の支持体１２２上にトランスデューサを配置して全周囲を囲む状態に近づけてもよい。なお、トランスデューサの配置や数及び支持体１２２の形状は被検体１００に応じて最適化すればよく、本発明に関してはあらゆる受信部１２０を採用することができる。
【００４７】
受信部１２０と被検体１００との間の空間は、光音響波が伝播することができる媒質で満たす。この媒質には、音響波が伝搬でき、被検体１００やトランスデューサ１２１との界面において音響特性が整合し、できるだけ光音響波の透過率が高い材料を採用する。例えば、この媒質には、水、超音波ジェルなどを採用することができる。
【００４８】
図４は、プローブ１８０の側面図を示す。本実施形態に係るプローブ１８０は、開口を有する半球状の支持体１２２に複数のトランスデューサ１２１が３次元に配置された受信部１２０を有する。また、支持体１２２の底部には、光学系１１２の光射出部が配置されている。
【００４９】
本実施形態においては、図４に示すように被検体１００は、保持部２００に接触することにより、その形状が保持される。
【００５０】
受信部１２０と保持部２００の間の空間は、光音響波が伝播することができる媒質（音響マッチング材）で満たされる。この媒質には、光音響波が伝搬でき、被検体１００やトランスデューサ１２１との界面において音響特性が整合し、できるだけ光音響波の透過率が高い材料を採用する。例えば、この媒質には、水、超音波ジェルなどを採用することができる。
【００５１】
保持手段としての保持部２００は被検体１００の形状を測定中に保持するために使用される。保持部２００により被検体１００を保持することによって、被検体１００の動きの抑制および被検体１００の位置を保持部２００内に留めることができる。保持部２００の材料には、ポリカーボネートやポリエチレン、ポリエチレンテレフタレート等、樹脂材料を用いることができる。つまり、保持部２００は、被検体１００を保持できる高度を有する材料であることが好ましい。また、保持部２００は、測定に用いる光を透過する材料であってもよい。さらに、保持部２００は、インピーダンスが被検体１００と同程度の材料で構成されていてもよい。また、乳房等の曲面を有するものを被検体１００とする場合、凹型に成型した保持部２００であってもよい。この場合、保持部２００の凹部分に被検体１００を挿入することができる。
【００５２】
保持部２００は、取り付け部２０１に取り付けられている。取り付け部２０１は、被検体の大きさに合わせて複数種類の保持部２００を交換可能に構成されていてもよい。例えば、取り付け部２０１は、曲率半径や曲率中心などの異なる保持部に交換できるように構成されていてもよい。
【００５３】
また、保持部２００には保持部２００の情報が登録されたタグとタグに登録された情報を読み取る読み取り部が設置されていてもよい。例えば、タグには、保持部２００の曲率半径、曲率中心、音速、識別ＩＤ等の情報を登録することができる。タグに登録された情報は、読み取り部により読み出され、コンピュータ１５０に転送される。保持部２００が取り付け部２０１に取り付けられたときに容易にタグを読み取るために、読み取り部は取り付け部２０１に設置されていてもよい。例えば、タグはバーコードであり、読み取り部はバーコードリーダである。
【００５４】
（駆動部１３０）
駆動部１３０は、被検体１００と受信部１２０との相対位置を変更する部分である。駆動部１３０は、駆動力を発生させるステッピングモータなどのモータと、駆動力を伝達させる駆動機構と、受信部１２０の位置情報を検出する位置センサとを含む。駆動機構としては、リードスクリュー機構、リンク機構、ギア機構、油圧機構、などを用いることができる。また、位置センサとしては、エンコーダー、可変抵抗器、リニアスケール、磁気センサ、赤外線センサ、超音波センサなどを用いたポテンショメータなどを用いることができる。
【００５５】
なお、駆動部１３０は被検体１００と受信部１２０との相対位置をＸＹ方向（２次元）に変更させるものに限らず、１次元または３次元に変更させてもよい。
【００５６】
移動経路は平面的にスパイラル状やライン＆スペースで走査してもよいし、さらに三次元的に体表に沿うように傾けてもよい。また、被検体１００の表面からの距離を一定に保つようにしてプローブ１８０を移動させてもよい。このとき駆動部１３０は、モータの回転数をモニターするなどしてプローブの移動量を計測してもよい。
【００５７】
なお、駆動部１３０は、被検体１００と受信部１２０との相対的な位置を変更できれば、受信部１２０を固定し、被検体１００を移動させてもよい。被検体１００を移動させる場合は、被検体１００を保持する保持部を動かすことで被検体１００を移動させる構成などが考えられる。また、被検体１００と受信部１２０の両方を移動させてもよい。
【００５８】
駆動部１３０は、相対位置を連続的に移動させてもよいし、ステップアンドリピートによって移動させてもよい。駆動部１３０は、プログラムされた軌跡で移動させる電動ステージであってもよいし、手動ステージであってもよい。
【００５９】
また、本実施形態では、駆動部１３０は光照射部１１０と受信部１２０を同時に駆動して走査を行っているが、光照射部１１０だけを駆動したり、受信部１２０だけを駆動したりしてもよい。
【００６０】
なお、プローブ１８０が、把持部が設けられたハンドヘルドタイプである場合、光音響装置１１００は駆動部１３０を有していなくてもよい。
【００６１】
（信号収集部１４０）
信号収集部１４０は、トランスデューサ１２１から出力されたアナログ信号である電気信号を増幅するアンプと、アンプから出力されたアナログ信号をデジタル信号に変換するＡ／Ｄ変換器とを含む。信号収集部１４０から出力されるデジタル信号は、コンピュータ１５０に記憶される。信号収集部１４０は、Ｄａｔａ  Ａｃｑｕｉｓｉｔｉｏｎ  Ｓｙｓｔｅｍ（ＤＡＳ）とも呼ばれる。本明細書において電気信号は、アナログ信号もデジタル信号も含む概念である。なお、フォトダイオードなどの光検出センサが、光照射部１１０から光射出を検出し、信号収集部１４０がこの検出結果をトリガーに同期して上記処理を開始してもよい。
【００６２】
（コンピュータ１５０）
情報処理装置としてのコンピュータ１５０は、画像処理装置１３００と同様のハードウェアで構成されている。すなわち、コンピュータ１５０の演算機能を担うユニットは、ＣＰＵやＧＰＵ（Ｇｒａｐｈｉｃｓ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）等のプロセッサ、ＦＰＧＡ（Ｆｉｅｌｄ  Ｐｒｏｇｒａｍｍａｂｌｅ  Ｇａｔｅ  Ａｒｒａｙ）チップ等の演算回路で構成されることができる。これらのユニットは、単一のプロセッサや演算回路から構成されるだけでなく、複数のプロセッサや演算回路から構成されていてもよい。
【００６３】
コンピュータ１５０の記憶機能を担うユニットは、ＲＡＭ（Ｒａｎｄｏｍ  Ａｃｃｅｓｓ  Ｍｅｍｏｒｙ）などの揮発性の媒体であってもよい。なお、プログラムが格納される記憶媒体は、非一時記憶媒体である。なお、コンピュータ１５０の記憶機能を担うユニットは、１つの記憶媒体から構成されるだけでなく、複数の記憶媒体から構成されていてもよい。
【００６４】
コンピュータ１５０の制御機能を担うユニットは、ＣＰＵなどの演算素子で構成される。コンピュータ１５０の制御機能を担うユニットは、光音響装置の各構成の動作を制御する。コンピュータ１５０の制御機能を担うユニットは、入力部１７０からの測定開始などの各種操作による指示信号を受けて、光音響装置の各構成を制御してもよい。すなわち、コンピュータ１５０は、ユーザーからの入力を受け付ける受付手段の一例に相当する。また、コンピュータ１５０の制御機能を担うユニットは、記憶機能を担うユニットに格納されたプログラムコードを読み出し、光音響装置の各構成の作動を制御する。すなわち、コンピュータ１５０は、本実施形態に係るシステムの制御装置として機能することができる。
【００６５】
本実施形態では、コンピュータ１５０は、その機能的な構成として、ワンショットボリュームデータ生成部１５１と参照画像生成部１５２と動画生成部１５３を備える。
【００６６】
なお、コンピュータ１５０と画像処理装置１３００は同じハードウェアで構成されていてもよい。１つのハードウェアがコンピュータ１５０と画像処理装置１３００の両方の機能を担っていてもよい。すなわち、コンピュータ１５０が、画像処理装置１３００の機能を担ってもよい。また、画像処理装置１３００が、情報処理装置としてのコンピュータ１５０の機能を担ってもよい。
【００６７】
（ワンショットボリュームデータ生成部１５１）
ワンショットボリュームデータ生成部１５１は、光照射ごとにそれぞれ異なる位置で取得した受信信号を用いて画像再構成処理を実行し、複数のボリュームデータ（３次元医用画像データ）を生成する。すなわち、ワンショットボリュームデータ生成部１５１は、被検体１００への光照射により発生する光音響波に由来するボリュームデータである光音響画像データを取得する第１の取得手段に相当する。この画像再構成は光照射ごとに、得られた光音響信号に対して個別に実施するため、再構成される光音響画像データをワンショットボリュームデータと呼ぶ。
【００６８】
（参照画像生成部１５２）
参照画像生成部１５２は、動画条件を決定するために参照する参照画像データを生成する。なお、参照画像データはボリュームデータである。参照画像生成部１５２は、ワンショットボリュームデータ生成部１５１で生成したワンショットボリュームデータを複数合成することにより、ＳＮ比を向上させた参照用の画像データを生成する。すなわち、参照画像生成部１５２は、光音響画像データを少なくとも２つ以上合成することにより合成画像データを取得する第２の取得手段の一例に相当する。
【００６９】
動画条件とは、レンダリング処理により２次元画像データを生成する際の対象領域、最大値投影方向、合成に用いるワンショットボリュームデータの数など、ワンショットボリュームデータ群から動画を生成するための条件である。ユーザーは、決定した動画条件を例えば入力部１７０を介して動画生成部１５３に入力する。動画生成部１５３は、入力された動画条件に基づいて動画を生成する。
【００７０】
（動画生成部１５３）
動画生成部１５３は、指定された動画条件に基づいて、動画（動画データ）を生成する。動画生成部１５３は、指定された数のワンショットボリュームデータを合成して合成画像データを生成する。なお、合成画像データの生成に用いるワンショットボリュームデータの数は少ないほど時間分解能の高い画像が得られる。一方で、合成画像データの生成に用いるワンショットボリュームデータの数が多いほどＳＮ比の高い画像が得られる。また、観察対象の血管に合わせて所望のＳＮ比と時間分解能が得られるワンショットボリュームデータの合成数を指定することで、観察に適した動画を得ることができる。
【００７１】
続いて、動画生成部１５３は、指定された領域を合成画像データから切り出し、指定された方向からの最大値投影画像データを生成する。
【００７２】
合成画像データに、毛や皮膚などが含まれる場合、これらのボクセル値が観察対象の血管よりも大きいと、血管は画像化されずに毛や皮膚などが画像化される。また、観察対象の血管が生体の深部にある場合、ボクセル値の大きな表在血管が画像化され、深部血管が観察できない場合がある。このような場合には、毛や皮膚、表在血管などが画像化領域に含まれないように、合成画像データから画像化領域を指定することで、観察に適した画像を得ることができる。
【００７３】
最後に、動画生成部１５３は、生成した最大輝度投影画像を撮影の時系列順に並べて、動画用画像データを生成する。
【００７４】
（表示部１６０）
表示部１６０は、参照画像生成部１５２で生成された参照画像データに基づいた画像および動画生成部１５３で生成された動画用画像データに基づいた画像を表示する。動画用画像データに基づいた画像を順番に更新表示することで動画として表示する。動画生成部１５３は、表示部１６０に動画用画像データに基づいた画像を表示させる表示制御手段として機能してもよい。
【００７５】
また、表示部１６０は、液晶ディスプレイや有機ＥＬ（Ｅｌｅｃｔｒｏ  Ｌｕｍｉｎｅｓｃｅｎｃｅ）ＦＥＤ、メガネ型ディスプレイ、ヘッドマウントディスプレイなどのディスプレイである。また、表示部１６０は、コンピュータ１５０により得られた被検体情報等に基づく画像や特定位置の数値等を表示する装置である。さらに、表示部１６０は、画像や装置を操作するためのＧＵＩを表示してもよい。
【００７６】
なお、表示部１６０と表示装置１４００は同じディスプレイであってもよい。すなわち、１つのディスプレイが表示部１６０と表示装置１４００の両方の機能を担っていてもよい。
【００７７】
（入力部１７０）
入力部１７０としては、ユーザーが操作可能な、マウスやキーボードなどで構成される操作コンソールを採用することができる。また、表示部１６０をタッチパネルで構成し、表示部１６０を入力部１７０として利用してもよい。
【００７８】
入力部１７０は、観察したい位置や深さの情報などを入力できるように構成されていてもよい。入力方法としては、数値を入力してもよいし、スライダーバーを操作することにより入力してもよい。また、入力された情報に応じて表示部１６０に表示される画像が更新されていってもよい。これにより、ユーザーは自身の操作によって決定されたパラメータにより生成された画像を確認しながら、適切なパラメータに設定できる。
【００７９】
なお、入力部１７０と入力装置１５００は同じ装置であってもよい。すなわち、１つの装置が入力部１７０と入力装置１５００の両方の機能を担っていてもよい。
【００８０】
（被検体１００）
被検体１００はシステムを構成するものではないが、以下に説明する。本実施形態に係るシステムは、人や動物の悪性腫瘍や血管疾患などの診断や化学治療の経過観察などを目的として使用できる。よって、被検体１００としては、生体、具体的には人体や動物の乳房や各臓器、血管網、頭部、頸部、腹部、手指または足指を含む四肢などの診断の対象部位が想定される。例えば、人体が測定対象であれば、オキシヘモグロビンあるいはデオキシヘモグロビンやそれらを多く含む血管あるいは腫瘍の近傍に形成される新生血管などを光吸収体の対象としてもよい。また、頸動脈壁のプラークなどを光吸収体の対象としてもよい。また、皮膚等に含まれるメラニン、コラーゲン、脂質などを光吸収体の対象としてもよい。さらに、生体を模したファントムを被検体１００としてもよい。
【００８１】
なお、光音響装置の各構成はそれぞれ別の装置として構成されてもよいし、一体となった１つの装置として構成されてもよい。また、光音響装置の少なくとも一部の構成が一体となった１つの装置として構成されてもよい。
【００８２】
なお、本実施形態に係るシステムを構成する各装置は、それぞれが別々のハードウェアで構成されていてもよいし、全ての装置が１つのハードウェアで構成されていてもよい。
【００８３】
すなわち、本実施形態に係るシステムの機能は、いかなるハードウェアで構成されていてもよい。
【００８４】
以下、本実施形態に係るシステムが実行する画像処理方法のフローを説明する。図３は、本実施形態に係るシステムが実行する画像処理方法のフロー図である。以下の本実施形態における光音響装置は、人や動物の血管疾患や悪性腫瘍などの診断や化学治療の経過観察などを主な目的とする。よって、被検体１００としては生体の一部が想定される。
【００８５】
（ステップＳ５０１：ワンショットボリュームデータの生成）
ステップＳ５０１において、本実施形態に係る光音響装置は、被検体を撮影することにより複数のワンショットボリュームデータを取得する。
【００８６】
具体的には、コンピュータ１５０は、ステップＳ３０１で指定された制御パラメータに基づいて、駆動部１３０に、プローブ１８０を指定の位置へ移動させる。ステップＳ３１０において複数位置での撮像が指定された場合には、駆動部１３０は、まずプローブ１８０を最初の指定位置へ移動させる。
【００８７】
そして、光照射部１１０は、Ｓ３１０で指定された制御パラメータに基づいて、被検体１００および機能情報マーカー１０１に光を照射する。
【００８８】
光源１１１から発生した光は、光学系１１２を介してパルス光として被検体１００に照射される。そして、被検体１００内部でパルス光が吸収され、光音響効果により光音響波が生じる。光照射部１１０はパルス光の伝送と併せて信号収集部１４０へ同期信号を送信する。
【００８９】
次に、信号収集部１４０は、光照射部１１０から送信された同期信号を受信すると、信号収集の動作を開始する。すなわち、信号収集部１４０は、受信部１２０から出力された、音響波に由来するアナログ電気信号を、増幅・ＡＤ変換することにより、増幅されたデジタル電気信号を生成し、コンピュータ１５０へ出力する。コンピュータ１５０は、信号収集部１４０から送信された信号を記憶部に保存する。なお、コンピュータ１５０は、発光をトリガーとして、発光時の受信部１２０の位置情報を駆動部１３０の位置センサからの出力に基づいて取得し、記憶してもよい。
【００９０】
そして、コンピュータ１５０を構成するワンショットボリュームデータ生成部１５１は、記憶部に記憶された信号データを用いて再構成処理を実行し、ワンショットボリュームデータを取得し、記憶部に保存する。この再構成処理は、光照射ごとに光音響信号に対して個別に実施される。なお、このとき、コンピュータ１５０は、信号データに加え、プローブ１８０の位置などの制御パラメータに基づいて、ワンショットボリュームデータを取得してもよい。
【００９１】
信号データを空間分布としてのボリュームデータに変換する再構成アルゴリズムとしては、タイムドメインでの逆投影法やフーリエドメインでの逆投影法などの解析的な再構成法やモデルベース法（繰り返し演算法）を採用することができる。例えば、タイムドメインでの逆投影法として、Ｕｎｉｖｅｒｓａｌ  ｂａｃｋ−ｐｒｏｊｅｃｔｉｏｎ（ＵＢＰ）、Ｆｉｌｔｅｒｅｄ  ｂａｃｋ−ｐｒｏｊｅｃｔｉｏｎ（ＦＢＰ）、または整相加算（Ｄｅｌａｙ−ａｎｄ−Ｓｕｍ）などが挙げられる。
【００９２】
また、ワンショットボリュームデータ生成部１５１は、被検体１００の内部における被検体１００に照射された光の光フルエンス分布を計算し、初期音圧分布を光フルエンス分布で除算することにより、吸収係数分布情報を取得してもよい。この場合、吸収係数分布情報を光音響画像データとして取得してもよい。ワンショットボリュームデータ生成部１５１は、光を吸収、散乱する媒質における光エネルギーの挙動を示す輸送方程式や拡散方程式を数値的に解く方法により、被検体１００の内部における光フルエンスの空間分布を算出することができる。数値的に解く方法としては、有限要素法、差分法、モンテカルロ法等を採用することができる。例えば、ワンショットボリュームデータ生成部１５１は、式（２）に示す光拡散方程式を解くことにより、被検体１００の内部における光フルエンスの空間分布を算出してもよい。
【００９３】
【数２】
2020162746.tif 000004
【００９４】
ここで、Ｄは拡散係数、μａは吸収係数、Ｓは照射光の入射強度、φは到達する光フルエンス、ｒは位置、ｔは時間を示す。
【００９５】
また、本工程でワンショットボリュームデータ生成部１５１は、複数の波長の光のそれぞれに対応する吸収係数分布情報を取得してもよい。そして、ワンショットボリュームデータ生成部１５１は、複数の波長の光のそれぞれに対応する吸収係数分布情報に基づいて、被検体１００を構成する物質の濃度の空間分布情報を機能情報として取得してもよい。
【００９６】
なお、モダリティとは別の装置である情報処理装置としてのコンピュータ１５０が、本実施形態に係る画像処理方法を実行してもよい。この場合、コンピュータ１５０が、事前にモダリティで生成された画像データを、ＰＡＣＳ（Ｐｉｃｔｕｒｅ  Ａｒｃｈｉｖｉｎｇ  ａｎｄ  Ｃｏｍｍｕｎｉｃａｔｉｏｎ  Ｓｙｓｔｅｍ）等の記憶部から読み出すことにより取得し、この画像データに対して本実施形態に係る画像処理方法を適用する。このように、本実施形態に係る画像処理方法は、事前に生成された画像データに対しても適用することができる。
【００９７】
ワンショットボリュームデータ生成部１５１で生成されたワンショットボリュームデータは参照画像生成部１５２および動画生成部１５３へ送られる。
【００９８】
（ステップＳ５０２：参照画像データの生成）
ステップＳ５０２において、参照画像生成部１５２は、ユーザーが動画条件を決定するための判断材料となる参照画像データを、ステップＳ５０１で生成したワンショットボリュームデータを用いて生成する。
【００９９】
参照画像生成部１５２は、複数のワンショットボリュームデータを加算、加算平均、相加相乗平均するなどして合成してもよい。これらの処理により、アーチファクトが抑制された３次元の参照画像データが生成される。この参照画像データに基づく画像を表示部１６０に表示させ、ユーザーからの動画条件に関する指示を受け付けることができる。なお、ここで生成する参照画像データは複数であってもよく、合成に用いるワンショットボリュームデータの数を変えたときのパターンに応じて参照画像データを複数生成してもよい。
【０１００】
このとき、例えば、参照画像生成部１５２は、複数のワンショットボリュームデータのうち、観察対象となる血管の変動が少ないボリュームを選択的に合成することにより参照画像データを生成してもよい。すなわち、被検体１００とプローブ１８０の相対位置の変位量が互いに小さい光音響画像データを合成することにより合成画像データを取得してもよい。具体的には、ユーザーが、ワンショットボリュームデータから観察対象となる血管が含まれる領域を指定する。参照画像生成部１５２は、複数のボリュームデータに対して、指定された領域のボクセル値（画素値）を比較し、それぞれのボリュームデータの変位量を算出する。参照画像生成部１５２は、変位量の少ないボリュームデータを選択的に合成することにより参照画像を取得する。このような合成を行うことで、撮影中の位置変動による対象血管のぼけを抑制した参照画像データを得ることができる。
【０１０１】
参照画像生成部１５２で生成された参照画像データは表示部で参照画像として表示され、動画条件の決定に用いられる。
【０１０２】
（ステップＳ５０３：動画条件の決定）
ステップＳ５０３において、動画生成部１５３は、ユーザーは、表示部１６０に表示された参照画像を元に動画条件を決定する。ここで決定する動画条件とは、例えば、動画を生成する画像の領域、合成画像データを２次元画像データ化する際の視線方向（最大値投影の場合は、最大値投影方向）、合成に使用する画像データの数、または合成に使用する画像データの階調などを含む。
【０１０３】
ユーザーは、レンダリング後に動画のフレームとして使用されるレンダリング後の２次元画像データを、対象血管の観察に適したものとするために、２次元画像データにレンダリング処理する領域を３次元の参照画像データの中から決定してもよい。領域の指定は参照画像から座標値やボクセル番号などの数値を読み取ってもよいし、表示画面上の領域をマウスなどの操作コンソールによるドラッグ操作などで直接指定してもよい。レンダリング後の２次元画像が対象血管の観察に適したものとなるように、例えば、毛や皮膚を表す画素値が高い画素などが含まれないように画像化領域を決定する。すなわち、例えば、閾値以上の画素値を有する画素を含まない領域を、動画データを生成する画像領域として決定してもよい。
【０１０４】
レンダリング手法として、例えば、最大値投影法（ＭＩＰ：Ｍａｘｉｍｕｍ  Ｉｎｔｅｎｓｉｔｙ  Ｐｒｏｊｅｃｔｉｏｎ）がある。ユーザーは、参照画像を元に、最大値を投影する方向を任意に指定することで自由断面画像を得ることができる。ユーザーは、動画で観察する対象血管の観察が容易となる方向に最大値投影方向（ＭＩＰ化方向）を決定する。なお、レンダリング手法は、上記に限定されず公知の種々の手法により実現される。
【０１０５】
また、合成に用いる画像データの数は、観察対象となる領域の特徴量や合成画像データのＳＮ比と時間分解能を考慮して決定してもよい。一般的に、合成に用いるボリュームデータの数が多いほどアーチファクトの低減効果およびシステムノイズの低減効果が高まり、画質の良い合成ボリューム画像が確認できるが、時間分解能が低下する。そのため、例えば、観察対象となる領域が太い血管である場合は、合成画像データのＳＮ比が低くても観察可能なため、合成にかかる時間を短縮するために合成に用いる画像データの数を少なくする。一方、観察対象となる領域が細い血管である場合は、よりＳＮ比が高い画像を必要とするため、合成に用いる画像データの数を多くする。すなわち、観察対象となる血管が細くなるほど、合成する画像データの数を多くしてもよい。なお、上記は一例であって、必ずしも血管の細さに比例して画像データの数を増加させなくてもよい。
【０１０６】
さらに、合成に用いるボリュームデータの数を変更したときの合成ボリューム画像と共に、その数を表示してもよい。これにより、ユーザーへの利便性を向上させることができる。  また、合成に用いている光音響信号の取得時間幅から算出できる周期（Ｈｚ）を表示してもよい。これにより、得られる画像の時間分解能を認識しやすくなり、ユーザーが合成に用いるボリュームデータの数の決定をしやすくなる。
【０１０７】
ユーザーは決定した動画条件を、入力部１７０を介して動画生成部１５３に入力する。また一連のユーザーが行う指定動作は、表示部１６０に表示されたＧＵＩを操作することにより行われてもよい。
【０１０８】
（ステップＳ５０４：合成画像データの生成）
動画生成部１５３は、ステップＳ５０１で生成した複数のワンショットボリュームデータを合成して合成画像データを生成する。動画生成部１５３は、複数のワンショットボリュームデータを加算、加算平均、相加相乗平均するなどして合成してもよい。これらの処理により、アーチファクトが抑制された合成画像データが生成される。また、ボリュームデータが光音響波に由来するものである場合、動画生成部１５３は、互いに異なる波長の光照射によって生成された複数の波長に対応する複数のワンショットボリュームデータを合成することにより、分光情報を取得してもよい。例えば、２波長の光に対応する吸収係数の空間分布を示すワンショットボリュームデータを用いて、分光情報として酸素飽和度の空間分布を取得する。この場合、動画生成部１５３は、式（１）にしたがって酸素飽和度ＳＯ２を取得する合成処理を行ってもよい。
【０１０９】
以下、図６を用いて、信号収集部１４０、コンピュータ１５０によって行われるデータ処理の関係を説明する。被検体１００とプローブ１８０の相対位置がＰｏｓ１、Ｐｏｓ２，Ｐｏｓ３，・・・ＰｏｓＮのときに取得される光音響信号データをＳｉｇ１、Ｓｉｇ２、Ｓｉｇ３、・・・ＳｉｇＮで示す。各相対位置で光音響信号が受信回路系によって取得されると、取得された光音響信号を用いた画像再構成により、ワンショットボリュームデータ生成部１５１がワンショットボリュームデータＶ１、Ｖ２、・・・、ＶＮを生成する。ここで、参照画像生成部１５２は、複数のワンショットボリュームデータを用いて参照画像用ボリュームデータＶｒｅｆを生成する。被検体１００とプローブ１８０の相対位置がそれぞれ異なる複数のワンショットボリュームデータを合成することで、アーチファクトを抑制した画像を生成することができる。また、複数のワンショットボリュームのうち、観察対象となる血管の変動が少ないボリュームを選択的に合成することで、撮影中の変動による対象血管のぼけを抑制した参照画像データを得ることができる。生成された参照画像データは表示部１６０に表示され、ユーザーにより動画条件が決定される。
【０１１０】
動画生成部１５３は、決定された動画条件に基づき、複数のワンショットボリュームデータを用いて合成画像データＶｉｎｔ１、Ｖｉｎｔ２、・・・ＶｉｎｔＮ−２を生成する。合成画像データの生成に用いるワンショットボリュームデータの数が少ないほど時間分解能の高い画像が得られ、逆に合成画像データの生成に用いるワンショットボリュームデータの数が多いほどＳＮ比の高い画像が得られる。動画生成部１５３は、入力部１７０を介して合成数や合成対象のボリュームデータの組み合わせの情報を取得し、合成画像の生成に用いるワンショットボリュームデータの組み合わせを変更することができる。
【０１１１】
（ステップＳ５０５：動画用画像データの生成）
ステップＳ５０５は、動画生成部１５３は、ステップＳ５０４で生成した３次元の合成画像データをレンダリングすることにより、動画の１フレームに対応する動画用２次元画像データを生成する。複数の合成画像データに対してレンダリングを行い、動画用画像データとして時系列順に並べられた２次元画像データ群を生成する。すなわち、画像処理装置１０２は、合成画像データＶｉｎｔ１、Ｖｉｎｔ２，・・・，ＶｉｎｔＮ−２から、２次元の動画用画像データＩｍ１，Ｉｍ２，・・・，ＩｍＮ−２を生成する。
【０１１２】
レンダリング手法としては、最大値投影法（ＭＩＰ：Ｍａｘｉｍｕｍ  Ｉｎｔｅｎｓｉｔｙ  Ｐｒｏｊｅｃｔｉｏｎ）、ボリュームレンダリング、及びサーフェイスレンダリングなどのあらゆる方法を採用することができる。ここで、３次元の合成画像データを２次元画像にレンダリングする際の領域および視線方向は、ユーザーが指定し動画条件により決定される。
【０１１３】
上記で説明した動作を一般化して説明すると、本実施例においては、被検体１００に対する相対位置の異なるＮ箇所（Ｎは３以上の整数）のそれぞれで光音響信号を取得し、これによってＮ個のボリュームデータを得る。
【０１１４】
そして、Ｎ個のボリュームデータのうち任意の数のボリュームデータを合成して参照画像データを得る。次に、参照画像データに基づいて動画条件を決定する。そして、決定された動画条件に基づいて、第ｉ〜第（ｉ＋ｍ）のワンショットボリュームデータ（ｉ＋ｍ＜Ｎ；ｉ，ｍはともに自然数）のうちの少なくとも２つ以上のボリュームデータを合成することで第１の合成画像データが生成される。さらに、第（ｉ＋ｎ）〜第（ｉ＋ｎ＋ｍ）のワンショットボリュームデータ（ｎは自然数）のうちの少なくとも２つ以上のボリュームデータを合成することで第２の合成画像データが生成される。その後、第１の合成画像データに基づく画像を第２の合成画像データに基づく画像で更新することにより、第１の合成画像と第２の合成画像が順次表示される。
【０１１５】
（ステップＳ５０６：動画の表示）
表示部１６０は、ステップＳ５０５で生成された動画用画像データを順次更新表示する。ユーザーは表示される動画を観察し、所望の動画が得られたかどうかを判断する。ユーザーは、画像化領域、最大値投影方向、階調設定、観察対象のＳＮ比および時間変化の様子を確認し、変更の必要があると判断した場合には、ステップＳ５０３に進み、合成画像データに使用するワンショットボリュームデータの数など、動画生成条件の再設定を行う。
【０１１６】
（ステップＳ５０７：終了？）
ユーザーは、ステップＳ５０７で表示部１６０に表示された動画を観察した結果、動画生成条件の変更が必要であると判断した場合、ステップＳ５０３に進み、動画生成条件を再設定する。
【０１１７】
ユーザーは、動画条件の変更が必要であると判断した場合、入力部１７０を用いて動画条件の変更を指示することができる。動画生成部１５３は、入力部１７０を介して変更された動画条件の情報を受け取った場合、動画条件の再設定を行うことができる。そして、ステップＳ５０４に戻り、再度動画生成を行う。一方、動画条件の変更指示がない場合、設定済みの動画条件での動画表示を継続する。
【０１１８】
例えば、動画生成部１５３は、合成画像データに使用するワンショットボリュームデータの数や動画用画像データのレンダリング条件（画像化領域、最大値投影方向、諧調設定）などを動画条件として変更することができる。
【０１１９】
以上説明したように、本実施形態に係る光音響装置によれば、参照画像を元に動画条件を決定することができ、観察に適した動画を生成することができる。
【０１２０】
（第２の実施形態）
第１の実施形態では、参照画像の生成に用いるワンショットボリュームデータと、動画データの生成に用いるワンショットボリュームデータが同じ場合を説明した。しかしながら、ワンショットボリュームデータは、参照画像生成に用いるデータと動画データの生成に用いるデータとで異なっていてもよい。
【０１２１】
図７は本実施例における動画生成までのフローを表した図である。ここでは第１の実施形態との差分について説明する。
【０１２２】
（ステップＳ７０１：参照画像用ワンショットボリュームデータを生成する工程）
ステップＳ７０１において、光音響装置は、被検体１００が指定の姿勢に設置された後、被検体１００に対して光を照射する。ここでは参照画像を生成することを目的に、動画用の計測とは別に、参照画像用の計測を行う。
【０１２３】
動画データは時間分解能を高めるために、１フレームの画像をできるだけ少ない数のワンショットボリュームで合成する。動画データで可視化する被検体１００の領域を、プローブ１８０の高感度領域と同等のサイズ程度とすることで、少ない合成数でＳＮ比を向上させることが可能となる。動画で可視化される領域がプローブ１８０の高感度領域程度の比較的狭い領域であることで、観察対象となる血管が動画可視化領域に適切に入らない場合がある。
【０１２４】
ここで、参照画像は動画よりも広い領域を可視化する画像とする。光音響装置は、参照画像用に動画データ用計測よりも広い範囲で被検体情報を取得する。参照画像は静止画であるため、時間分解能の制約がなく、多数の位置でのワンショットボリュームデータを使用した広い領域の画像とすることができる。動画よりも広い領域を可視化した参照画像を用いることで、観察対象血管を容易に探索することができる。探索した観察対象血管の位置を指定して動画用データの計測を行うことで、観察対象血管を適切に動画領域に含めることができる。
【０１２５】
また、参照画像用計測は動画用計測と異なる波長で行われてもよい。例えば動画計測よりもメラニンの光吸収が大きい短い波長で計測することで、皮膚をより強く可視化し、動画領域の指定を容易にするなどの方法が考えられる。
【０１２６】
ワンショットボリュームデータ生成部１５１は参照画像用のワンショットボリュームデータを生成し、参照画像生成部１５２に出力する。
【０１２７】
（ステップＳ７０３：動画条件の決定）
ユーザーは、表示部１６０に表示された参照画像を元に動画条件を決定する。ここで動画条件とは、画像化領域、ＭＩＰ化方向、合成画像に使用するワンショットボリュームの数に加えて、動画計測を行う位置情報を含む。本実施例では動画計測を行う位置を決定した後に動画用データ生成のための被検体計測を開始する。
【０１２８】
ユーザーは、表示部１６０に表示された参照画像を参照して、動画条件を決定する。表示部１６０は、動画化される領域の位置とサイズを示す枠線を、参照画像に重畳表示することで、ユーザーが動画計測位置を決定する際の利便性を向上させることができる。決定した動画計測位置情報はシステム制御部に入力され、動画用の計測が行われる。
【０１２９】
（ステップＳ７０４：動画用ワンショットボリュームデータの生成）
光音響装置は、動画計測位置情報を元に動画用計測を行う。ワンショットボリュームデータ生成部１５１は、プローブ１８０が受信した光音響波の受信信号を用いて画像再構成処理を実行し、ワンショットボリュームデータを生成する。
【０１３０】
（ステップＳ７０９：動画条件の再設定）
ユーザーは、表示部１６０に表示された動画を観察した結果、ステップＳ７０８で動画生成条件の変更が必要であると判断した場合、ステップＳ７０９に進み、動画生成条件を再設定する。例えば、動画生成部１５３は、合成画像データに使用するワンショットボリュームデータの数や動画用画像データのレンダリング条件（画像化領域、最大値投影方向、諧調設定）などを動画条件として変更することができる。
【０１３１】
以上説明したように、本実施形態に係る光音響装置によれば、参照画像を元に動画条件を決定することができ、観察に適した動画を生成することができる。
【図１】
2020162746.tif 000005
【図２】
2020162746.tif 000006
【図３】
2020162746.tif 000007
【図４】
2020162746.tif 000008
【図５】
2020162746.tif 000009
【図６】
2020162746.tif 000010
【図７】
2020162746.tif 000011