JP6415650B2 - 機能情報取得装置、機能情報取得方法、及びプログラム - Google Patents
機能情報取得装置、機能情報取得方法、及びプログラム Download PDFInfo
- Publication number
- JP6415650B2 JP6415650B2 JP2017133072A JP2017133072A JP6415650B2 JP 6415650 B2 JP6415650 B2 JP 6415650B2 JP 2017133072 A JP2017133072 A JP 2017133072A JP 2017133072 A JP2017133072 A JP 2017133072A JP 6415650 B2 JP6415650 B2 JP 6415650B2
- Authority
- JP
- Japan
- Prior art keywords
- image data
- image
- light
- processing apparatus
- resolution
- 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.)
- Expired - Fee Related
Links
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
Description
酸素飽和度とは、血液中の全ヘモグロビンに対する酸素と結合したヘモグロビンの含有率である。酸素飽和度を検出することによって心肺機能が正常に作動しているか否かを測定することができる。また、酸素飽和度は腫瘍の良悪性を見分ける指標になることから、悪性腫瘍の効率的発見手段として期待されている。
ただし、この方法によって機能情報を取得するには、同一箇所について行った複数回の測定結果を比較演算する必要があり、体の動きなどにより測定位置が合致しなかった場合、演算の結果が誤ったものとなるという可能性がある。
図1は、本発明の光音響イメージングの第一の実施形態について示したものである。ここでは、図1に基づいて、本発明を実施するための最良の形態について説明する。
まず、波長Aの光を用いた測定のとき、電気信号処理回路9から送られてきたデジタル信号を吸収係数算出部101において再構成することにより吸収係数分布Aを示すデータ(第1のデータ)が算出され(S3)、算出された吸収係数分布Aを示すデータはメモリA102に蓄えられる。また、波長Bの光を用いた測定についても同様に吸収係数分布Bを示すデータ(第3のデータ)が算出され(S6)、メモリB103に蓄えられる。次に、波長Aの光を用いた測定における光吸収体5の位置と波長Bの光を用いた測定における光吸収体5の位置との間の位置ずれ量をずれ量入力装置11に入力し、その値に基づき分解能変化量決定部108において分解能変化量を決定する(S7)。分解能変化処理部104は、メモリに蓄えられた吸収係数分布を示すデータのうち少なくとも一つの吸収係数分布を示すデータについて、決定された分解能変化量だけ、画像データ上の画像空間分解能を低減させ、低減後の吸収係数分布を得る(S8)。本発明は、この、分解能が低減された吸収係数分布のデータ(第2のデータ)を用いて酸素飽和度などの被検体情報を求めることを特徴とする。なお、本発明は2次元の画像データ(ピクセルデータ)でも3次元の画像データ(ボクセルデータ)にも適用できる。なお、画像空間分解能とは、音響検出器7の素子のサイズで決まる分解能ではなく、画像空間での分解能のことである。本明細書中において3次元画像データ上の空間的な分解能をボクセル空間分解能と呼び、2次元画像データ上の空間的な分解能をピクセル空間分解能と呼ぶこととする。また、ボクセル空間分解能とピクセル空間分解能とをあわせて画像空間分解能と定義する。図2ではメモリB103に蓄えられた吸収係数分布B(波長Bの光を照射した際に算出され得る吸収係数分布)を示すデータの画像空間分解能を低減しているが、複数の吸収係数分布のデータのうちいずれの画像空間分解能を低減してもよい。また、複数の吸収係数分布のデータの画像空間分解能を低減してもよい。
なお、画像空間分解能を変化させる吸収係数分布の数に関わらず、分解能を低減させた結果の低減後の吸収係数分布における光吸収体の像の範囲内に、分解能を低減していない吸収係数分布(全ての吸収係数分布の分解能を低減した場合は、分解能を低減する前のいずれか1つの吸収係数分布)における光吸収体の像が少なくとも含まれるように、位置ずれ量に対する画像空間分解能の変化量を決定する。このとき、画像空間分解能を低減させるそれぞれの吸収係数分布について独立して変化量を決定してもよいし、画像空間分解能を低減させるすべての吸収係数分布について一律に変化量を定めてもよい。位置ずれ量を得る方法は特に限定されず、いかなる公知の方法によっても得ることができる。また、位置ずれ量を機械的な測定や画像からの測定から得てもよいし、入力は手動、自動のいずれでもよい。なお、位置ずれ量に対して変化させるべき画像空間分解能の量は、分解能を変化させる手法ごとに異なるので、あらかじめ分解能を変化させる手法ごとに位置ずれ量と分解能を変化させる量との関係をテーブル化もしくは関係式化しておき、このあらかじめ用意されたテーブルもしくは関係式を用いて分解能変化量を決定してもよい。
酸素飽和度の算出方法は後述する。
この結果を、表示装置12で表示する(S11)。
得られる吸収係数分布の画像空間分解能を低減させて第2のデータを得るために、[実施形態1]で述べた空間フィルターに代わる手段として、音響波検出器で得られた信号に帯域制限を加える手法について図2および図5を用いて述べる。
複数の光を用いた測定の間に光吸収体の位置ずれがない場合、位置ずれがあるが位置ずれに関する処理を行わない場合、位置ずれがあるが[実施形態1]を実施した場合のそれぞれについて酸素飽和度を算出するシミュレーションをした。
実施例1と同様のシミュレーションを行い、吸収係数分布のボクセル空間分解能の低減手法として、複数の位置で得られた音響信号を足し合わせて一箇所の信号として扱う方法を用いた例について説明する。
このときの計算時間は従来手法に比べて無視できるほどの増加しかなかった。
2 光
3 被検体
4 光学部品
5 光吸収体
6 音響波
7 音響検出器
8 制御装置
9 電気信号処理回路
10 データ処理装置
11 ずれ量入力装置
12 表示装置
Claims (24)
- 第1の光が被検体に照射されることにより発生する第1の音響波に由来する第1の画像データを取得し、第2の光が前記被検体に照射されることにより発生する第2の音響波に由来する第2の画像データを取得する画像データ取得手段と、
前記第1の画像データの空間分解能を低減することにより第3の画像データを生成し、前記第2の画像データの空間分解能を低減することにより第4の画像データを生成する分解能低減手段と、
前記第3の画像データおよび前記第4の画像データを用いて第5の画像データを生成する画像データ生成手段と、
前記第5の画像データに基づいた第1の画像を表示手段に表示させる表示制御手段と、
を有し、
前記表示制御手段は、前記第1の画像データまたは前記第2の画像データに基づいて、前記第1の画像の表示態様を決定する
ことを特徴とする処理装置。 - 第1の光が被検体に照射されることにより発生する第1の音響波に由来する第1の画像データを取得し、第2の光が前記被検体に照射されることにより発生する第2の音響波に由来する第2の画像データを取得する画像データ取得手段と、
前記第1の画像データの空間分解能を低減することにより第3の画像データを生成する分解能低減手段と、
前記第2の画像データおよび前記第3の画像データを用いて第5の画像データを生成する画像データ生成手段と、
前記第5の画像データに基づいた第1の画像を表示手段に表示させる表示制御手段と、
を有し、
前記表示制御手段は、前記第1の画像データまたは前記第2の画像データに基づいて、前記第1の画像の表示態様を決定する
ことを特徴とする処理装置。 - 前記画像データ生成手段は、前記第3の画像データと前記第4の画像データとを除算することにより前記第5の画像データを生成する
ことを特徴とする請求項1に記載の処理装置。 - 前記画像データ生成手段は、前記第2の画像データと前記第3の画像データとを除算することにより前記第5の画像データを生成する
ことを特徴とする請求項2に記載の処理装置。 - 前記分解能低減手段は、前記第1の画像データと前記第2の画像データとの位置ずれ量を取得し、前記位置ずれ量に基づいて前記第1の画像データに対する空間分解能の変化量を決定する
ことを特徴とする請求項1から4のいずれか1項に記載の処理装置。 - 前記分解能低減手段は、前記第1の画像データにフィルターを畳み込むことにより前記第3の画像データを生成する
ことを特徴とする請求項1から5のいずれか1項に記載の処理装置。 - 前記フィルターは移動平均フィルターまたはガウスフィルターである
ことを特徴とする請求項6に記載の処理装置。 - 前記第1の画像データは、前記第1の音響波に由来する光吸収係数分布を示すデータであり、
前記第2の画像データは、前記第2の音響波に由来する光吸収係数分布を示すデータである
ことを特徴とする請求項1から7のいずれか1項に記載の処理装置。 - 前記表示制御手段は、前記第1の画像と、前記第1の画像データまたは前記第2の画像データに基づいた第2の画像とを前記表示手段に重畳して表示させる
ことを特徴とする請求項1から8のいずれか1項に記載の処理装置。 - 前記表示制御手段は、前記第1の画像データまたは前記第2の画像データに基づいて、前記第1の画像の色属性を決定することを特徴とする請求項1から9のいずれか1項に記載の処理装置。
- 前記色属性は、明度であることを特徴とする請求項10に記載の処理装置。
- 前記表示制御手段は、前記第5の画像データの値と、前記第1の画像データまたは前記第2の画像データの値とを色相、明度、及び彩度のうち互いに異なる少なくとも1つに対応させた前記第1の画像を前記表示手段に表示させる
ことを特徴とする請求項1から9のいずれか1項に記載の処理装置。 - 前記表示制御手段は、前記第1の画像データまたは前記第2の画像データの値を明度に、前記第5の画像データの値を色相または彩度に対応させた前記第1の画像を前記表示手段に表示させる
ことを特徴とする請求項1から9のいずれか1項に記載の処理装置。 - 前記表示制御手段は、前記第1の画像データまたは前記第2の画像データにおいて閾値より値が低いピクセルまたはボクセルに対応する前記第5の画像データが表示されないように前記第1の画像の表示態様を決定する
ことを特徴とする請求項1から13のいずれか1項に記載の処理装置。 - 前記表示制御手段は、前記第1の画像データまたは前記第2の画像データにおいて閾値より値が高いピクセルまたはボクセルに対応する前記第5の画像データが選択的に表示されるように前記第1の画像の表示態様を決定する
ことを特徴とする請求項1から13のいずれか1項に記載の処理装置。 - 前記閾値は、前記第1の画像データまたは前記第2の画像データの最大値に基づいた値である
ことを特徴とする請求項14または15に記載の処理装置。 - 前記表示制御手段は、前記第1の画像データまたは前記第2の画像データに基づいて光吸収体の位置を決定し、当該光吸収体の位置に対応しない前記第5の画像データが表示されないように前記第1の画像の表示態様を決定する
ことを特徴とする請求項1から16のいずれか1項に記載の処理装置。 - 前記表示制御手段は、前記第1の画像データまたは前記第2の画像データに基づいて光吸収体の位置を決定し、当該光吸収体の位置に対応する前記第5の画像データが選択的に表示されるように前記第1の画像の表示態様を決定する
ことを特徴とする請求項1から16のいずれか1項に記載の処理装置。 - 前記光吸収体は、ヘモグロビンを含む
ことを特徴とする請求項17または18に記載の処理装置。 - 前記第5の画像データは、酸素飽和度分布を示すデータである
ことを特徴とする請求項1から19のいずれか1項に記載の処理装置。 - 前記第2の光の波長は、前記第1の光の波長とは異なる
ことを特徴とする請求項1から20のいずれか1項に記載の処理装置。 - 光照射により発生する音響波に由来する画像データを処理する処理装置が行う処理方法であって、
第1の光が被検体に照射されることにより発生する第1の音響波に由来する第1の画像データを取得し、
第2の光が前記被検体に照射されることにより発生する第2の音響波に由来する第2の画像データを取得し、
前記第1の画像データの空間分解能を低減することにより第3の画像データを生成し、
前記第2の画像データの空間分解能を低減することにより第4の画像データを生成し、
前記第3の画像データおよび前記第4の画像データを用いて第5の画像データを生成し、
前記第5の画像データに基づいた第1の画像を表示手段に表示させ、
前記第1の画像データまたは前記第2の画像データに基づいて、前記第1の画像の表示態様を決定する
ことを特徴とする処理方法。 - 光照射により発生する音響波に由来する画像データを処理する処理装置が行う処理方法であって、
第1の光が被検体に照射されることにより発生する第1の音響波に由来する第1の画像データを取得し、
第2の光が前記被検体に照射されることにより発生する第2の音響波に由来する第2の画像データを取得し、
前記第1の画像データの空間分解能を低減することにより第3の画像データを生成し、
前記第2の画像データおよび前記第3の画像データを用いて第5の画像データを生成し、
前記第5の画像データに基づいた第1の画像を表示手段に表示させ、
前記第1の画像データまたは前記第2の画像データに基づいて、前記第1の画像の表示態様を決定する
ことを特徴とする処理方法。 - 請求項22または23に記載の処理方法の各工程をコンピュータに実行させるためのプログラム。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010022892 | 2010-02-04 | ||
JP2010022892 | 2010-02-04 |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2015195955A Division JP6173402B2 (ja) | 2010-02-04 | 2015-10-01 | 機能情報取得装置、機能情報取得方法、及びプログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2017170243A JP2017170243A (ja) | 2017-09-28 |
JP6415650B2 true JP6415650B2 (ja) | 2018-10-31 |
Family
ID=55265081
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2015195955A Expired - Fee Related JP6173402B2 (ja) | 2010-02-04 | 2015-10-01 | 機能情報取得装置、機能情報取得方法、及びプログラム |
JP2017133072A Expired - Fee Related JP6415650B2 (ja) | 2010-02-04 | 2017-07-06 | 機能情報取得装置、機能情報取得方法、及びプログラム |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2015195955A Expired - Fee Related JP6173402B2 (ja) | 2010-02-04 | 2015-10-01 | 機能情報取得装置、機能情報取得方法、及びプログラム |
Country Status (1)
Country | Link |
---|---|
JP (2) | JP6173402B2 (ja) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2018050775A (ja) | 2016-09-27 | 2018-04-05 | キヤノン株式会社 | 光音響装置、情報処理方法、及びプログラム |
JP2018050776A (ja) | 2016-09-27 | 2018-04-05 | キヤノン株式会社 | 光音響装置、情報処理方法、及びプログラム |
JP6759032B2 (ja) | 2016-09-27 | 2020-09-23 | キヤノン株式会社 | 光音響装置、情報処理方法、及びプログラム |
CN110403576B (zh) * | 2019-08-01 | 2024-01-30 | 中国医学科学院北京协和医院 | 三维光声成像在乳腺肿瘤评分系统中的应用 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5982915A (en) * | 1997-07-25 | 1999-11-09 | Arch Development Corporation | Method of detecting interval changes in chest radiographs utilizing temporal subtraction combined with automated initial matching of blurred low resolution images |
JP2004350962A (ja) * | 2003-05-29 | 2004-12-16 | Olympus Corp | 超音波診断装置のデータ処理方法 |
JP4583068B2 (ja) * | 2004-05-11 | 2010-11-17 | 株式会社東芝 | 超音波診断装置 |
WO2007088709A1 (ja) * | 2006-01-31 | 2007-08-09 | Kansai Technology Licensing Organization Co., Ltd. | 光音響断層撮影装置及び光音響断層撮影方法 |
WO2008103982A2 (en) * | 2007-02-23 | 2008-08-28 | The Regents Of The University Of Michigan | System and method for monitoring photodynamic therapy |
JP5284129B2 (ja) * | 2008-02-06 | 2013-09-11 | キヤノン株式会社 | イメージング装置、及び解析方法 |
JP5235477B2 (ja) * | 2008-04-14 | 2013-07-10 | キヤノン株式会社 | 超音波による画像形成方法及び超音波装置 |
-
2015
- 2015-10-01 JP JP2015195955A patent/JP6173402B2/ja not_active Expired - Fee Related
-
2017
- 2017-07-06 JP JP2017133072A patent/JP6415650B2/ja not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JP6173402B2 (ja) | 2017-08-02 |
JP2016019847A (ja) | 2016-02-04 |
JP2017170243A (ja) | 2017-09-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5818444B2 (ja) | 機能情報取得装置、機能情報取得方法、及びプログラム | |
JP6320594B2 (ja) | 被検体情報取得装置及び被検体情報取得方法 | |
JP6669521B2 (ja) | マルチスペクトル光音響イメージングのためのデバイスおよび方法 | |
JP5704998B2 (ja) | 光音響装置およびその制御方法 | |
JP6415650B2 (ja) | 機能情報取得装置、機能情報取得方法、及びプログラム | |
US9339254B2 (en) | Object information acquiring apparatus | |
JP6071260B2 (ja) | 被検体情報取得装置および情報処理方法 | |
EP2638850B1 (en) | Subject information obtaining device, subject information obtaining method, and program | |
JP6362301B2 (ja) | 被検体情報取得装置、および、被検体情報取得装置の作動方法 | |
US20190216331A1 (en) | Object information acquiring apparatus and display method | |
KR20170068388A (ko) | 광음향 장치, 표시 제어 방법 및 기억매체 | |
JP2015142740A (ja) | 光音響装置、情報処理装置および表示方法 | |
US10695006B2 (en) | Apparatus and display control method | |
JP2017000897A (ja) | 方法 | |
JP6425438B2 (ja) | 被検体情報取得装置および画像処理方法 | |
JP6776115B2 (ja) | 処理装置および処理方法 | |
D’Alessandro et al. | Depth-dependent hemoglobin analysis from multispectral transillumination images | |
EP3655972B1 (en) | Imaging method, controller and imaging system, for monitoring a patient post evar | |
JP6490171B2 (ja) | データ処理装置およびデータ処理方法 | |
JP6419281B2 (ja) | 情報処理装置および方法 | |
JP7142865B2 (ja) | 癌発生疑い部位特定装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20170804 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20170804 |
|
RD05 | Notification of revocation of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7425 Effective date: 20171214 |
|
RD04 | Notification of resignation of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7424 Effective date: 20180126 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20180516 |
|
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: 20180822 |
|
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: 20180904 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20181002 |
|
R151 | Written notification of patent or utility model registration |
Ref document number: 6415650 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |
|
LAPS | Cancellation because of no payment of annual fees |