JP2018001013A - データ処理装置およびデータ処理方法 - Google Patents
データ処理装置およびデータ処理方法 Download PDFInfo
- Publication number
- JP2018001013A JP2018001013A JP2017197818A JP2017197818A JP2018001013A JP 2018001013 A JP2018001013 A JP 2018001013A JP 2017197818 A JP2017197818 A JP 2017197818A JP 2017197818 A JP2017197818 A JP 2017197818A JP 2018001013 A JP2018001013 A JP 2018001013A
- Authority
- JP
- Japan
- Prior art keywords
- subject
- information
- region
- interest
- distribution
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000012545 processing Methods 0.000 title description 20
- 238000003672 processing method Methods 0.000 title description 7
- 238000009826 distribution Methods 0.000 claims abstract description 157
- 238000004364 calculation method Methods 0.000 claims abstract description 47
- 230000031700 light absorption Effects 0.000 claims abstract description 19
- 230000001678 irradiating effect Effects 0.000 claims abstract description 5
- 238000000034 method Methods 0.000 claims description 29
- 238000005259 measurement Methods 0.000 claims description 22
- 230000003287 optical effect Effects 0.000 claims description 9
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 abstract description 110
- 229910052760 oxygen Inorganic materials 0.000 abstract description 110
- 239000001301 oxygen Substances 0.000 abstract description 110
- 239000012491 analyte Substances 0.000 abstract 10
- 238000010521 absorption reaction Methods 0.000 description 63
- 239000006096 absorbing agent Substances 0.000 description 12
- 239000000126 substance Substances 0.000 description 11
- 108010054147 Hemoglobins Proteins 0.000 description 10
- 102000001554 Hemoglobins Human genes 0.000 description 10
- 238000012937 correction Methods 0.000 description 9
- 239000008280 blood Substances 0.000 description 8
- 210000004369 blood Anatomy 0.000 description 8
- 238000003325 tomography Methods 0.000 description 7
- 238000003384 imaging method Methods 0.000 description 6
- 108010064719 Oxyhemoglobins Proteins 0.000 description 5
- XUMBMVFBXHLACL-UHFFFAOYSA-N Melanin Chemical compound O=C1C(=O)C(C2=CNC3=C(C(C(=O)C4=C32)=O)C)=C2C4=CNC2=C1C XUMBMVFBXHLACL-UHFFFAOYSA-N 0.000 description 4
- 239000003921 oil Substances 0.000 description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 230000006870 function Effects 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 238000004611 spectroscopical analysis Methods 0.000 description 3
- 102000008186 Collagen Human genes 0.000 description 2
- 108010035532 Collagen Proteins 0.000 description 2
- 206010028980 Neoplasm Diseases 0.000 description 2
- 239000002033 PVDF binder Substances 0.000 description 2
- 210000004204 blood vessel Anatomy 0.000 description 2
- 229920001436 collagen Polymers 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 150000002632 lipids Chemical class 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 229920002981 polyvinylidene fluoride Polymers 0.000 description 2
- 230000008929 regeneration Effects 0.000 description 2
- 238000011069 regeneration method Methods 0.000 description 2
- 238000009877 rendering Methods 0.000 description 2
- 230000008827 biological function Effects 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 201000011510 cancer Diseases 0.000 description 1
- 239000004359 castor oil Substances 0.000 description 1
- 235000019438 castor oil Nutrition 0.000 description 1
- 239000000919 ceramic Substances 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000002872 contrast media Substances 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000000499 gel Substances 0.000 description 1
- ZEMPKEQAKRGZGQ-XOQCFJPHSA-N glycerol triricinoleate Natural products CCCCCC[C@@H](O)CC=CCCCCCCCC(=O)OC[C@@H](COC(=O)CCCCCCCC=CC[C@@H](O)CCCCCC)OC(=O)CCCCCCCC=CC[C@H](O)CCCCCC ZEMPKEQAKRGZGQ-XOQCFJPHSA-N 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 239000003068 molecular probe Substances 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 238000006213 oxygenation reaction Methods 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 229920000306 polymethylpentene Polymers 0.000 description 1
- 239000011116 polymethylpentene Substances 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 239000011347 resin Substances 0.000 description 1
- 229920005989 resin Polymers 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000012306 spectroscopic technique Methods 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Images
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
Description
近赤外光は、生体の大部分を構成する水を透過しやすく、血液中のヘモグロビンで吸収されやすい性質を持つため、パルス光に近赤外光を用いることで、血管像を取得することができる。さらに、異なる波長のパルス光によって取得された血管像を比較することによって、機能情報である血液中の酸素飽和度を測定することができる。悪性腫瘍周辺の血液は良性腫瘍周辺の血液より酸素飽和度が低くなっていると考えられているので、酸素飽和度を知ることによって腫瘍の良悪鑑別を行えるようになることが期待されている。
非特許文献1には、測定した被検体内の吸収係数に基づいて、信号部分、すなわち血液の存在が推定される部分と、その他の部分であるバックグラウンド部分を区別し、信号部分の酸素飽和度分布のみを作成および表示する光音響トモグラフィ装置が開示されている。
波長λの光を用いた測定によって得られる、被検体内の吸収係数μa(λ)は、酸化ヘ
モグロビンの吸収係数εox(λ)と酸化ヘモグロビンの存在比Coxの積と、還元ヘモグロビンの吸収係数εde(λ)と還元ヘモグロビンの存在比Cdeの積の和である。すなわち、μa(λ)は式1で表すことができる。
されている。式1における未知数は、CoxおよびCdeの二つであるので、異なる波長の光を用いて測定を少なくとも二回行うことで、連立方程式によってCoxおよびCdeを算出することができる。さらに多くの波長で測定を行った場合には、最小自乗法によるフィッティングによってCoxおよびCdeを得ることができる。
そして、酸素飽和度SO2は、全ヘモグロビン中の酸化ヘモグロビンの比なので、式2
によって算出することができる。
素飽和度分布が三次元である場合、三次元空間内の分布を二次元画像としてディスプレイ上に表示するため、一目で分布が捉えにくく、酸素飽和度値を読み取ることがさらに困難になる。
被検体に異なる波長のパルス光をそれぞれ照射する光照射手段と、前記パルス光に起因して前記被検体から発生した音響波を前記波長ごとに受信し、電気信号に変換する音響波受信手段と、前記変換された電気信号に基づいて、前記被検体内の各位置における光吸収情報を前記波長ごとに算出する光吸収情報算出手段と、前記被検体内の各位置における光吸収情報を前記波長ごとに比較することで、前記被検体内の各位置における特性情報を示す情報である特性分布を取得する特性情報算出手段と、前記被検体内における特性分布からヒストグラムを生成する統計量算出手段と、前記被検体内における特性分布と、前記ヒストグラムを、同一の表示装置に表示させるための画像情報を生成し、出力する画像生成手段と、を有することを特徴とする。
被検体に計測光を照射し、前記被検体内で発生した光音響波を受信し、当該光音響波に基づいて、生体内部の光学特性に関する情報を可視化した画像を生成する被検体情報取得装置が行う被検体情報取得方法であって、異なる波長のパルス光をそれぞれ照射された被検体から発生した音響波を前記波長ごとに受信し、電気信号に変換する音響波受信ステップと、前記変換された電気信号に基づいて、被検体内の各位置における光吸収情報を前記波長ごとに算出する光吸収情報算出ステップと、前記被検体内の各位置における光吸収情報を前記波長ごとに比較することで、前記被検体内の各位置における特性情報を示す情報である特性分布を取得する特性情報算出ステップと、前記被検体内における特性分布からヒストグラムを生成する統計量算出ステップと、前記被検体内における特性分布と、前記ヒストグラムを、同一の表示装置に表示させるための画像情報を生成し、出力する画像生成ステップとを含むことを特徴とする。
本発明の第一の実施形態に係る光音響診断装置は、計測光を被検体に照射し、計測光に起因して発生した光音響波を解析することで、被検体である生体内部の酸素飽和度の分布を画像化する装置である。具体的には、図2(A)〜(D)に示すように、酸素飽和度分布を画像表示し、設定された関心領域中の統計量を表示する装置である。まず、構成要素について述べ、その後、処理方法を説明し、最後に効果について述べる。
酸素飽和度や吸収係数などは、被検体内の特性情報を表す被検体情報と言える。すなわち、実施形態に係る光音響診断装置は、被検体情報取得装置とも呼べる。
まず、図1を参照しながら、第一の実施形態に係る光音響診断装置の構成を説明する。第一の実施形態に係る光音響診断装置は、第一光源1、第二光源2、光照射装置3、被検体4、音響検出器5、電気信号処理装置6、データ処理装置7、関心領域設定装置12、表示装置14からなる。また、データ処理装置7には、画像再構成部8、光分布補正部9、比較演算部10、統計量算出部11、画像生成部13が含まれる。
第一光源1、および第二光源2はパルス光を発生させる装置である。各光源は、大出力を得るためレーザ光源であることが望ましいが、発光ダイオードなどでもよい。光音響波を効果的に発生させるためには、被検体の熱特性に応じて十分短い時間に光を照射させなければならない。被検体が生体である場合、光源から発生するパルス光のパルス幅は数十ナノ秒以下にすることが望ましい。また、パルス光の波長は、700nm〜1200nm程度であることが望ましい。この領域の光は比較的生体深部まで到達することができ、深部の情報を得ることができるため、生体の窓と呼ばれている。また、パルス光の波長は、観測対象に対して吸収係数が高いことが望ましい。
また、第一光源1によるパルス光の照射タイミングと、第二光源2によるパルス光の照射タイミングは異なるタイミングであることが好ましい。それぞれの波長によって、吸収係数を個別に測定しなければならないためである。
なお、本実施形態では光源を二つ用いるが、光源が波長可変レーザである場合、第一光源と第二光源を兼ねた単一の光源を使用してもよい。
光照射装置3は、第一光源1および第二光源2で発生したパルス光を被検体4へ導く装置である。具体的には、光ファイバーやレンズ、ミラー、拡散板などからなる光学機器である。これらの光学機器を用いて、パルス光の照射形状、光密度、被検体への照射方向などの照射条件を任意のものに設定することができる。なお、光照射装置3を構成する光学機器は、ここにあげたものだけに限定されず、上記のような機能を満たすものであれば、どのようなものであってもよい。
第一光源1、第二光源2、光照射手段3が、本発明における光照射手段である。
被検体4は、本発明を構成するものではないが、ここで説明する。
被検体4は、測定対象物である。被検体4は、典型的には生体であるが、生体の音響特性と光学特性を模擬したファントムであってもよい。光音響診断装置では、被検体4の内
部に存在する光吸収係数の大きい光吸収体をイメージングできる。被検体が生体である場合、イメージングの対象はヘモグロビン、水、メラニン、コラーゲン、脂質などである。被検体がファントムである場合は、前述した物質の光学特性を模擬した物質を光吸収体の代替として内部に封入する。さらに、造影剤や分子プローブなどを、被検体である生体またはファントム中に注入してもよい。
音響検出器5は、被検体4の内部で発生した音響波をアナログの電気信号に変換する手段である。音響検出器5は、単一の音響検出器からなってもよいし、複数の音響検出器からなってもよい。また、音響検出器5は、音響波の反射、減衰の影響を排除するために、被検体4と音響的に結合されるように設置される必要がある。例えば、音響検出器5と被検体4の間には音響マッチングジェルや水、オイルなどの音響整合材を設けることが望ましい。
また、音響検出器5は、感度が高く、周波数帯域が広いものが望ましい。具体的にはPZT(圧電セラミックス)、PVDF(ポリフッ化ビニリデン樹脂)、cMUT(容量性マイクロマシン超音波トランスデューサ)、ファブリペロー干渉計を用いたものなどが挙げられる。ただし、ここに挙げたものだけに限定されず、機能を満たすものであれば、どのようなものであってもよい。
電気信号処理装置6は、音響検出器5で得られたアナログの電気信号を増幅し、デジタル信号に変換する手段である。電気信号処理装置6は、具体的には電気回路から成るアンプやAnalog-Digital Converter(ADC)などである。効率的にデータを取得するため、音響検出器5の受信素子数と同じだけのアンプ及びADCを設けることが望ましいが、一つずつ用意したアンプ及びADCを順次につなぎ換えて使用してもよい。
音響検出器5、および電気信号処理装置6が、本発明における音響波受信手段である。
データ処理装置7は、デジタル化された信号または数値データを処理し、表示装置14に表示する画像を生成する手段である。データ処理装置7には、画像再構成部8、光分布補正部9、比較演算部10、統計量算出部11、画像生成部13が含まれている。
また、データ処理装置7は、不図示のCPUとDRAM、および不揮発メモリ、制御ポートから構成されるコンピュータである。不揮発メモリに格納されたプログラムがCPUで実行されることにより、各モジュールの制御が行われる。本実施形態では、データ処理装置7はコンピュータであるが、専用に設計されたハードウェアであってもよい。
被検体内の光強度分布は光の波長によって異なるので、吸収係数分布は波長ごとに計算される。本実施形態では、第一光源の波長に対応する吸収係数分布を第一の吸収係数分布、第二光源の波長に対応する吸収係数分布を第二の吸収係数分布と呼ぶ。
算出される吸収係数が、本発明における光吸収情報であり、画像再構成部8、および光分布補正部9が、本発明における光吸収情報算出手段である。
酸素飽和度の分布は、前述したように、式3に、得られた吸収係数分布を代入することで得ることができる。また、3つ以上の光源を用いた場合は、酸素飽和度分布は、最小二乗法によってCoxおよびCdeを求め、式2に代入することによって得ることができる。
比較演算部10が、本発明における特性情報算出手段であり、生成される酸素飽和度分布が、本発明における特性分布である。
統計量の算出は、画素ごとの酸素飽和度値を用いて行うことが望ましいが、ある範囲の画素集合の代表値や、画素を補間することによって得られた補間点での値を計算に用いてもよい。本実施形態では、画素ごとの値を用いて統計量を算出する。
統計量算出部11が、本発明における統計量算出手段である。
なお、生成される画像は、同一の表示装置によって表示されるものであればよく、必ずしも同一の画像内に全てを含める必要はない。例えば、メインウインドウに酸素飽和度分布を表示するための第一の画像と、サブウインドウに統計量を表示するための第二の画像を生成するようにしてもよい。
画像生成部13が、本発明における画像生成手段であり、画像生成部13が生成する画像が、本発明における画像情報である。
関心領域設定装置12は、比較演算部10が生成した酸素飽和度分布に対し、統計量を算出する範囲を設定する手段である。関心領域設定装置12は、具体的にはキーボード、マウス、タッチペン、タッチパネルなどの入力インターフェースや、ディスプレイを備えたコンピュータなどである。
関心領域は、操作者が着目している領域であるため、手動で設定されることが望ましいが、関心領域が前回の測定と同じである場合や、ある法則に基づいて設定できる場合は、所定の場所や法則をあらかじめ記憶しておき、自動的に設定してもよい。また、手動設定で大まかな設定を行い、詳細な設定は自動で行ってもよい。またはその逆であってもよい。手動設定と自動設定は組み合わせることができる。
関心領域設定装置12が、本発明における関心領域設定手段である。
また、表示装置14は、画像生成部13が生成した画像を表示する手段であり、典型的にはディスプレイ装置である。表示装置14は、関心領域を指定する際の出力インターフェースとして使用することもできる。表示装置14は、本発明の必須構成ではない。
次に、統計量算出部11が出力する統計量を例示し、それぞれを算出する方法を説明する。
例えば、ヒストグラムを微分し、微分値の正負が反転する箇所の数をカウントすることで、当該値を算出することができる。酸素飽和度が異なる複数種類の吸収体が関心領域内に存在する場合、ヒストグラムは多峰性になり、関心領域内に単一種類の吸収体のみが存在する場合は、ヒストグラムは単峰性になる。
従って、ヒストグラムの単峰性、多峰性を表す値を統計量として算出することで、関心領域内にある吸収帯が、複数種類であるか単一種類であるかを判断することができる。また、ヒストグラムが多峰性である場合は、それぞれの峰ごとに統計量を算出することが望ましい。
このように、数値で表わされる統計量を要約統計量と称する。
例えば、関心領域が狭く、単一種類の吸収体しか存在しないと予測される領域など、酸素飽和度の変化が小さくなるべき領域において、酸素飽和度がばらついている場合、測定の誤差が大きいことがわかる。すなわち、酸素飽和度のばらつきを求めることで、測定の信頼性を得ることができる。
標準偏差は式7で計算され、分散は標準偏差の二乗で計算される。また、最大値と最小値の差は、関心領域内の全画素の酸素飽和度値のうち最大のものと最小の者の差分を計算することによって得られる。
次に、統計量算出部11が算出した各種統計量を画像化する例について、図2を参照しながら説明する。
図2は、画像生成部13によって生成され、表示装置14に出力される画像を表した例である。図2(A)〜(D)は、算出した酸素飽和度を輝度値で表した画像と、当該画像に重畳された各種統計量の例である。
信号部分とは、吸収体の像であると考えられる部分であり、バックグラウンド部分とは吸収体の像ではないと考えられる部分である。両者の区別は、非特許文献1のように、吸収係数分布に閾値を設けることによって行ってもよいし、吸収係数分布や初期音圧分布の信号の形状から信号らしさを算出し、それに閾値を設けることによって行ってもよい。
また、この際に用いる吸収係数分布や初期音圧分布は、酸化ヘモグロビンと還元ヘモグロビンの吸収係数が近い波長の光源の吸収係数分布や、初期音圧分布を用いることが望ましい。また、信号らしさを算出する場合には、データ処理装置7の中に信号確率算出部を設け、そこで処理するようにしてもよい。また、本実施形態では、信号部分とバックグラウンド部分とを分離して画像を生成しているが、分離は必ずしも行わなくてもよい。
また、図2(A)〜(D)に示した画像は例示である。画像生成部13は、統計量算出部11が算出できる統計量であれば、任意の統計量を組み合わせて画像を生成することができる。
第一の実施形態に係る光音響装置が行う処理について、処理フローチャートである図3を参照しながら説明する。
まず、ステップS1で、第一光源1が被検体にパルス光を照射する。
そして、ステップS2で、被検体から発生した音響波を音響検出器5が受信し、ステップS3で、電気信号処理装置6が当該音響波をデジタル信号に変換する。そして、画像再構成部8が、当該デジタル信号から初期音圧分布を生成し、光分布補正部9が、第一の吸収係数分布を生成する。生成された第一の吸収係数分布は、比較演算部10および画像生成部13へ送信される。
次に、ステップS7で、比較演算部10が、得られた第一及び第二の吸収係数分布を用いて酸素飽和度分布を生成する。生成された酸素飽和度分布は、統計量算出部11および画像生成部13へ送信される。
ステップS8では、画像生成部13が、光分布補正部9から受信した吸収係数分布を用いて、比較演算部10から受信した酸素飽和度に対して、信号部分とバックグラウンド部分の分離を実施する。そして、信号部分のみの酸素飽和度分布からなる画像を生成し、表示装置14へ出力する。この時点では、表示装置14に表示される画像には、統計量は重畳されていない。ステップS8を処理すると、操作者による関心領域の指定が可能になる。
そして、ステップS10で、統計量算出部11が、比較演算部10から受信した酸素飽和度分布と、関心領域設定装置12から受信した関心領域の情報を用いて、関心領域内部の統計量を算出する。算出された統計量は、画像生成部13に送信される。
ステップS11では、画像生成部13が、受信した統計量を用いて、表示画像を更新する。具体的には、受信した統計量を表す図形情報または文字情報を生成し、ステップS8で生成した画像に重畳させ、新しい画像を生成する。生成された画像は、表示装置14へ送信される。これにより、酸素飽和度分布と、関心領域内部の統計量の両方が操作者に提示される。
第二の実施形態は、操作者に対して、酸素飽和度分布ではなく、吸収係数分布を提示する実施形態である。
第二の実施形態に係る光音響診断装置のシステム構成を図4に示す。第二の実施形態に係る光音響診断装置のシステム構成は、第一の実施形態と比較して、比較演算部10から画像生成部13へ酸素飽和度分布が送信されないという点のみが相違する。
第二の実施形態では、ステップS7で、比較演算部10が、得られた第一及び第二の吸収係数分布を用いて酸素飽和度分布を生成する。生成された酸素飽和度分布は、統計量算出部11のみに送信される。
また、ステップS8では、画像生成部13が、光分布補正部9から受信した吸収係数分布を用いて、吸収係数分布を表す画像を生成する。この結果、操作者には、吸収係数分布を表す画像のみが提示される。
すなわち、ステップS9で関心領域設定装置12が取得する関心領域は、吸収係数分布に基づいて設定されたものとなる。吸収係数分布と酸素飽和度分布は同じ座標系であるので、吸収係数分布を見ながら設定した関心領域は、酸素飽和度分布に対する関心領域として使用することができる。
ステップS10およびS11では、第一の実施形態と同様に、統計量算出部11が、関心領域内の酸素飽和度分布の統計量を計算し、画像生成部13が、統計量を含んだ画像を生成し、表示装置14に出力する。統計量は、第一の実施形態と同様に、任意のものを用いることができる。
なお、本実施形態では、吸収係数分布のみを表示したが、酸素飽和度分布と吸収係数分布のどちらかを選択的に表示させるようにしてもよい。また、表示する吸収係数分布は、第一の吸収係数分布であってもよいし、第二の吸収係数分布であってもよい。
本発明の効果を実験にて確認した。本実施例では、被検体として半球状のファントムを使用した。ファントム母材の音響特性及び光学特性は、生体に近いものを使用し、ファントムの内部には、被検体表面から25mmの位置に光吸収体を設置した。光吸収体の光吸収係数は、ファントム母材に対して約5倍のものを使用した。光吸収体は均一であり、酸素飽和度の平均は80.3%である。
また、ファントムの両側を厚さ10mmのポリメチルペンテンからなる保持版二枚で挟み込んで密着させ、さらに片側の保持板の向こう側に1mmの油層を設け、油層を通して音響検出器を密着させた。
また、音響検出器をXYステージに接続し、音響検出器の面と同じ面方向に走査することができるようにした。
また、第一および第二光源にはTiSレーザを用い、第一光源の波長を756nm、第二光源の波長を797nmとした。この光源によるナノ秒オーダーのパルス光を、音響検出器と同じ面と、被検体を挟んで逆の面から同時に被検体に照射した。
そして、バックプロジェクションを用いて画像再構成、光分布補正、酸素飽和度の算出を行い、3次元データの初期音圧分布、吸収係数分布、酸素飽和度分布を得た。
信号部分とバックグラウンド部分の分離には第二の吸収係数分布を用い、信号部分の酸素飽和度分布を表示させて、関心領域を設定した。
認することができ、また、平均値と標準偏差を表示することによって、定量的な情報による診断を行えることが確認できた。
なお、各実施形態の説明は本発明を説明する上での例示であり、本発明は、発明の趣旨を逸脱しない範囲で適宜変更または組み合わせて実施することができる。本発明は、上記処理の少なくとも一部を含む被検体情報取得装置の制御方法として実施することもできるし、これらの方法を被検体情報取得装置に実行させるプログラムとして実施することもできる。上記処理や手段は、技術的な矛盾が生じない限りにおいて、自由に組み合わせて実施することができる。
また、各実施形態では光源をそれぞれ二つずつ使用したが、波長が異なる光源をさらに加え、三つ以上の光源を使用するようにしてもよい。
また、本発明の課題も、算出方法が同じである以上、酸素飽和度特有のものではなく、分光を用いた含有率算出技術に共通のものである。従って、本発明の範囲は酸素飽和度のみに限定されず、分光を用いる他の物質の含有率を含みうる。本発明の説明における酸素飽和度は、他の物質の含有率と読み替えることができる。
被検体を構成する物質の含有量または含有率の空間分布を表す第一の画像データ、および、前記含有量または含有率とは異なる前記被検体に関する情報の空間分布を表す第二の画像データを取得する画像データ取得手段と、前記第二の画像データに基づいて、第一の関心領域を設定する関心領域設定手段と、前記第一の関心領域に対応する前記第一の画像データの統計情報を取得する統計情報取得手段と、前記統計情報に基づいた画像を表示手段に表示させる表示制御手段と、を有することを特徴とする。
また、本発明の第二の形態に係るデータ処理方法は、
被検体を構成する物質の含有量または含有率の空間分布を表す第一の画像データ、および、前記含有量または含有率とは異なる前記被検体に関する情報の空間分布を表す第二の画像データを取得する画像データ取得手段と、前記第二の画像データに基づいた画像を表示手段に表示させる表示制御手段と、前記表示手段に前記第二の画像データに基づいた画像が表示されているときに操作者により指定された第一の関心領域に関する情報に基づいて、前記第一の関心領域を設定する関心領域設定手段と、前記第一の画像データに基づいて、前記第一の関心領域に対応する前記第一の画像データの統計情報を取得する統計情報取得手段と、を有し、前記表示制御手段は、前記統計情報に基づいた画像を前記表示手段に表示させることを特徴とする。
被検体を構成する物質の含有量または含有率の空間分布を表す第一の画像データ、および、前記含有量または含有率とは異なる前記被検体に関する情報の空間分布を表す第二の画像データを取得し、前記第二の画像データに基づいて、第一の関心領域を設定し、前記第一の関心領域に対応する前記第一の画像データの統計情報を取得し、前記統計情報に基づいた画像を表示手段に表示させることを特徴とする。
また、本発明の別形態に係るデータ処理方法は、
被検体を構成する物質の含有量または含有率の空間分布を表す第一の画像データ、および、前記含有量または含有率とは異なる前記被検体に関する情報の空間分布を表す第二の画像データを取得し、前記第二の画像データに基づいた画像を表示手段に表示させ、前記表示手段に前記第二の画像データに基づいた画像が表示されているときに操作者により指定された第一の関心領域に関する情報に基づいて、前記第一の関心領域を設定し、前記第一の画像データに基づいて、前記第一の関心領域に対応する前記第一の画像データの統計情報を取得し、前記統計情報に基づいた画像を前記表示手段に表示させることを特徴とす
る。
Claims (11)
- 被検体に異なる波長のパルス光をそれぞれ照射する光照射手段と、
前記パルス光に起因して前記被検体から発生した音響波を前記波長ごとに受信し、電気信号に変換する音響波受信手段と、
前記変換された電気信号に基づいて、前記被検体内の各位置における光吸収情報を前記波長ごとに算出する光吸収情報算出手段と、
前記被検体内の各位置における光吸収情報を前記波長ごとに比較することで、前記被検体内の各位置における特性情報を示す情報である特性分布を取得する特性情報算出手段と、
前記被検体内における特性分布からヒストグラムを生成する統計量算出手段と、
前記被検体内における特性分布と、前記ヒストグラムを、同一の表示装置に表示させるための画像情報を生成し、出力する画像生成手段と、
を有することを特徴とする、被検体情報取得装置。 - 被検体内における関心領域を設定する関心領域設定手段をさらに有し、
前記統計量算出手段は、前記関心領域内の特性分布に基づいてヒストグラムを生成する
ことを特徴とする、請求項1に記載の被検体情報取得装置。 - 前記統計量算出手段は、前記生成したヒストグラムの単峰性、または多峰性を示す値をさらに算出し、
前記画像生成手段は、当該算出した値を前記画像情報にさらに含ませる
ことを特徴とする、請求項2に記載の被検体情報取得装置。 - 前記統計量算出手段は、前記関心領域内における前記特性分布の統計量を示す値をさらに算出し、
前記画像生成手段は、当該算出した値を前記画像情報にさらに含ませる
ことを特徴とする、請求項2または3のいずれか1項に記載の被検体情報取得装置。 - 前記統計量算出手段は、前記関心領域設定手段で設定された関心領域が変化するごとにヒストグラムの再生成を行い、
前記画像生成手段は、前記ヒストグラムが再生成されるごとに前記画像情報を再生成する
ことを特徴とする、請求項2から4のいずれか1項に記載の被検体情報取得装置。 - 被検体に計測光を照射し、前記被検体内で発生した光音響波を受信し、当該光音響波に基づいて、生体内部の光学特性に関する情報を可視化した画像を生成する被検体情報取得装置が行う被検体情報取得方法であって、
異なる波長のパルス光をそれぞれ照射された被検体から発生した音響波を前記波長ごとに受信し、電気信号に変換する音響波受信ステップと、
前記変換された電気信号に基づいて、被検体内の各位置における光吸収情報を前記波長ごとに算出する光吸収情報算出ステップと、
前記被検体内の各位置における光吸収情報を前記波長ごとに比較することで、前記被検体内の各位置における特性情報を示す情報である特性分布を取得する特性情報算出ステップと、
前記被検体内における特性分布からヒストグラムを生成する統計量算出ステップと、
前記被検体内における特性分布と、前記ヒストグラムを、同一の表示装置に表示させるための画像情報を生成し、出力する画像生成ステップと
を含むことを特徴とする、被検体情報取得方法。 - 被検体内における関心領域を設定する関心領域設定ステップをさらに含み、
前記統計量算出ステップでは、前記関心領域内における特性分布に基づいてヒストグラムを生成する
ことを特徴とする、請求項6に記載の被検体情報取得方法。 - 前記統計量算出ステップでは、前記生成したヒストグラムの単峰性、または多峰性を示す値をさらに算出し、
前記画像生成ステップでは、当該算出した値を前記画像情報にさらに含ませる
ことを特徴とする、請求項7に記載の被検体情報取得方法。 - 前記統計量算出ステップでは、前記関心領域内における前記特性分布の統計量を示す値をさらに算出し、
前記画像生成ステップでは、当該算出した値を前記画像情報にさらに含ませる
ことを特徴とする、請求項7または8のいずれか1項に記載の被検体情報取得方法。 - 前記統計量算出ステップでは、前記関心領域設定ステップで設定された関心領域が変化するごとにヒストグラムの再生成を行い、
前記画像生成ステップでは、前記ヒストグラムが再生成されるごとに前記画像情報を再生成する
ことを特徴とする、請求項7から9のいずれか1項に記載の被検体情報取得方法。 - 請求項6から10のいずれか1項に記載の、被検体情報取得方法の各ステップをコンピュータに実行させるためのプログラム。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2017197818A JP6490171B2 (ja) | 2017-10-11 | 2017-10-11 | データ処理装置およびデータ処理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2017197818A JP6490171B2 (ja) | 2017-10-11 | 2017-10-11 | データ処理装置およびデータ処理方法 |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2012286549A Division JP6226523B2 (ja) | 2012-12-28 | 2012-12-28 | 被検体情報取得装置、表示方法、およびデータ処理装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2018001013A true JP2018001013A (ja) | 2018-01-11 |
JP6490171B2 JP6490171B2 (ja) | 2019-03-27 |
Family
ID=60945410
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2017197818A Expired - Fee Related JP6490171B2 (ja) | 2017-10-11 | 2017-10-11 | データ処理装置およびデータ処理方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP6490171B2 (ja) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2020024741A (ja) * | 2018-04-27 | 2020-02-13 | 株式会社東芝 | 読取システム、読取方法、プログラム、及び記憶媒体 |
US12087154B2 (en) | 2018-04-27 | 2024-09-10 | Kabushiki Kaisha Toshiba | Reading system, reading method, and storage medium |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH06292666A (ja) * | 1993-04-09 | 1994-10-21 | Fujitsu Ltd | 超音波診断装置 |
JP2002263079A (ja) * | 2001-03-06 | 2002-09-17 | Toshiba Corp | 磁気共鳴映像装置 |
JP2004121652A (ja) * | 2002-10-04 | 2004-04-22 | Aloka Co Ltd | 超音波診断装置 |
JP2005095376A (ja) * | 2003-09-25 | 2005-04-14 | Ge Medical Systems Global Technology Co Llc | 超音波診断装置 |
JP2006110028A (ja) * | 2004-10-14 | 2006-04-27 | Ge Medical Systems Global Technology Co Llc | 超音波イメージング装置 |
JP2011526514A (ja) * | 2008-06-30 | 2011-10-13 | ネルコア・ピユーリタン・ベネツト・アイルランド | 基点を用いて勾配を決定するための信号処理システムおよび方法 |
JP2012024227A (ja) * | 2010-07-22 | 2012-02-09 | Canon Inc | 画像情報取得装置、画像情報取得方法および画像情報取得プログラム |
JP2012217554A (ja) * | 2011-04-06 | 2012-11-12 | Canon Inc | 光音響装置およびその制御方法 |
EP2529665A1 (en) * | 2011-05-30 | 2012-12-05 | University of Graz | System and method of determining thickness of body tissue |
-
2017
- 2017-10-11 JP JP2017197818A patent/JP6490171B2/ja not_active Expired - Fee Related
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH06292666A (ja) * | 1993-04-09 | 1994-10-21 | Fujitsu Ltd | 超音波診断装置 |
JP2002263079A (ja) * | 2001-03-06 | 2002-09-17 | Toshiba Corp | 磁気共鳴映像装置 |
JP2004121652A (ja) * | 2002-10-04 | 2004-04-22 | Aloka Co Ltd | 超音波診断装置 |
JP2005095376A (ja) * | 2003-09-25 | 2005-04-14 | Ge Medical Systems Global Technology Co Llc | 超音波診断装置 |
JP2006110028A (ja) * | 2004-10-14 | 2006-04-27 | Ge Medical Systems Global Technology Co Llc | 超音波イメージング装置 |
JP2011526514A (ja) * | 2008-06-30 | 2011-10-13 | ネルコア・ピユーリタン・ベネツト・アイルランド | 基点を用いて勾配を決定するための信号処理システムおよび方法 |
JP2012024227A (ja) * | 2010-07-22 | 2012-02-09 | Canon Inc | 画像情報取得装置、画像情報取得方法および画像情報取得プログラム |
JP2012217554A (ja) * | 2011-04-06 | 2012-11-12 | Canon Inc | 光音響装置およびその制御方法 |
EP2529665A1 (en) * | 2011-05-30 | 2012-12-05 | University of Graz | System and method of determining thickness of body tissue |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2020024741A (ja) * | 2018-04-27 | 2020-02-13 | 株式会社東芝 | 読取システム、読取方法、プログラム、及び記憶媒体 |
US12087154B2 (en) | 2018-04-27 | 2024-09-10 | Kabushiki Kaisha Toshiba | Reading system, reading method, and storage medium |
Also Published As
Publication number | Publication date |
---|---|
JP6490171B2 (ja) | 2019-03-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6226523B2 (ja) | 被検体情報取得装置、表示方法、およびデータ処理装置 | |
JP6320594B2 (ja) | 被検体情報取得装置及び被検体情報取得方法 | |
US20190133451A1 (en) | Apparatus for acquiring biofunctional information, method for acquiring biofunctional information, and program therefor | |
US9032800B2 (en) | Photoacoustic imaging apparatus and photoacoustic imaging method | |
JP6504826B2 (ja) | 情報処理装置および情報処理方法 | |
JP6071260B2 (ja) | 被検体情報取得装置および情報処理方法 | |
US10653322B2 (en) | Photoacoustic apparatus, method of acquiring subject information, and non-transitory computer readable medium | |
US9330462B2 (en) | Object information acquiring apparatus and control method of object information acquiring apparatus | |
JP2013188310A (ja) | 被検体情報取得装置に用いる信号処理装置および被検体情報取得方法 | |
JP2014140718A (ja) | 被検体情報取得装置、表示方法、およびプログラム | |
JP2013188311A (ja) | 被検体情報取得装置に用いる信号処理装置および被検体情報取得方法 | |
JP6362420B2 (ja) | 被検体情報取得装置、被検体情報取得方法、および、プログラム | |
JP2017029610A (ja) | 光音響装置、信頼度取得方法、プログラム | |
US20160058295A1 (en) | Photoacoustic wave measurement apparatus and photoacoustic wave measurement method | |
JP6490171B2 (ja) | データ処理装置およびデータ処理方法 | |
EP3329843B1 (en) | Display control apparatus, display control method, and program | |
JP6425438B2 (ja) | 被検体情報取得装置および画像処理方法 | |
JP6486056B2 (ja) | 光音響装置および光音響装置の処理方法 | |
JP2017164222A (ja) | 処理装置および処理方法 | |
JP6512969B2 (ja) | 処理装置、光音響装置、処理方法、及びプログラム | |
JP6556300B2 (ja) | 被検体情報取得装置、被検体情報取得方法およびプログラム | |
JP6625182B2 (ja) | 情報処理装置、情報処理方法およびプログラム | |
JP6501820B2 (ja) | 処理装置、処理方法、及びプログラム | |
US20180299763A1 (en) | Information processing apparatus, object information acquiring apparatus, and information processing method | |
JP6391249B2 (ja) | 被検体情報取得装置および信号処理方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20171110 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20171110 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20180530 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20180626 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20180827 |
|
RD02 | Notification of acceptance of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7422 Effective date: 20181116 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20190129 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20190226 |
|
R151 | Written notification of patent or utility model registration |
Ref document number: 6490171 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |
|
LAPS | Cancellation because of no payment of annual fees |