JP2016077352A - 光音響装置および光音響装置の処理方法 - Google Patents
光音響装置および光音響装置の処理方法 Download PDFInfo
- Publication number
- JP2016077352A JP2016077352A JP2014209193A JP2014209193A JP2016077352A JP 2016077352 A JP2016077352 A JP 2016077352A JP 2014209193 A JP2014209193 A JP 2014209193A JP 2014209193 A JP2014209193 A JP 2014209193A JP 2016077352 A JP2016077352 A JP 2016077352A
- Authority
- JP
- Japan
- Prior art keywords
- absorption coefficient
- subject
- calculation method
- distribution
- photoacoustic apparatus
- 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
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0093—Detecting, measuring or recording by applying one single type of energy and measuring its conversion into another type of energy
- A61B5/0095—Detecting, measuring or recording by applying one single type of energy and measuring its conversion into another type of energy by applying light and detecting acoustic waves, i.e. photoacoustic measurements
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/145—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
- A61B5/14542—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue for measuring blood gases
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/43—Detecting, measuring or recording for evaluating the reproductive systems
- A61B5/4306—Detecting, measuring or recording for evaluating the reproductive systems for evaluating the female reproductive systems, e.g. gynaecological evaluations
- A61B5/4312—Breast evaluation or disorder diagnosis
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/70—Means for positioning the patient in relation to the detecting, measuring or recording means
- A61B5/708—Breast positioning means
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/1702—Systems in which incident light is modified in accordance with the properties of the material investigated with opto-acoustic detection, e.g. for gases or analysing solids
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/06—Visualisation of the interior, e.g. acoustic microscopy
- G01N29/0654—Imaging
- G01N29/0672—Imaging by acoustic tomography
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/22—Details, e.g. general constructional or apparatus details
- G01N29/24—Probes
- G01N29/2418—Probes using optoacoustic interaction with the material, e.g. laser radiation, photoacoustics
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
- G01N29/46—Processing the detected response signal, e.g. electronic circuits specially adapted therefor by spectral analysis, e.g. Fourier analysis or wavelet analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/1702—Systems in which incident light is modified in accordance with the properties of the material investigated with opto-acoustic detection, e.g. for gases or analysing solids
- G01N2021/1706—Systems in which incident light is modified in accordance with the properties of the material investigated with opto-acoustic detection, e.g. for gases or analysing solids in solids
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/01—Indexing codes associated with the measuring variable
- G01N2291/015—Attenuation, scattering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/024—Mixtures
- G01N2291/02466—Biological material, e.g. blood
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Pathology (AREA)
- General Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Immunology (AREA)
- General Physics & Mathematics (AREA)
- Biochemistry (AREA)
- Analytical Chemistry (AREA)
- Molecular Biology (AREA)
- Biomedical Technology (AREA)
- Animal Behavior & Ethology (AREA)
- Surgery (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- Medical Informatics (AREA)
- Heart & Thoracic Surgery (AREA)
- Acoustics & Sound (AREA)
- Signal Processing (AREA)
- Optics & Photonics (AREA)
- Mathematical Physics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Gynecology & Obstetrics (AREA)
- Reproductive Health (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
Description
また、これらを併用することで、被検体内の形態情報(例えば血管構造)と、機能情報(例えば酸素飽和度)の双方を取得することができる。
このような技術として、例えば、特許文献1には、画像上に発生したアーティファクトを削除する処理が記載されている。また、特許文献2には、モデルに基づいて信号を復元することで、真の信号部分を強調する為の手法が記載されている。また、特許文献3には、ブラインドデコンボリューションによって画像の視認性を向上させる手法が記載されている。
被検体に対して複数の波長のパルス光を照射可能な光源と、前記パルス光を照射された被検体から発生する音響波を受信し、電気信号に変換する音響波受信部と、前記電気信号に基づいて、第一の計算手法を用いて、前記被検体内における第一の吸収係数分布を取得する第一の情報取得手段と、前記電気信号に基づいて、第二の計算手法を用いて、前記被検体内における第二の吸収係数分布を取得する第二の情報取得手段と、複数の波長のパルス光をそれぞれ照射して得られた、複数の前記第一の吸収係数分布に基づいて、前記被検体内の機能情報の分布を算出する第三の情報取得手段と、前記機能情報の分布を、前記第二の吸収係数分布に基づいてマスクした画像を生成する画像生成手段と、を有し、前記第二の計算手法は、前記第一の計算手法と比べて、吸収係数分布を画像化した際により高い視認性が得られる手法であることを特徴とする。
被検体に対して複数の波長のパルス光を照射可能な光源と、前記パルス光を照射された被検体から発生する音響波を受信し、電気信号に変換する音響波受信部と、を有する光音響装置の処理方法であって、前記電気信号に基づいて、第一の計算手法を用いて、前記被検体内における第一の吸収係数分布を取得する第一の情報取得ステップと、前記電気信号に基づいて、第二の計算手法を用いて、前記被検体内における第二の吸収係数分布を取得する第二の情報取得ステップと、複数の波長のパルス光をそれぞれ照射して得られた、複数の前記第一の吸収係数分布に基づいて、前記被検体内の機能情報の分布を算出する第三の情報取得ステップと、前記機能情報の分布を、前記第二の吸収係数分布に基づいてマスクした画像を生成する画像生成ステップと、を含み、前記第二の計算手法は、前記第一の計算手法と比べて、吸収係数分布を画像化した際により高い視認性が得られる手法であることを特徴とする。
適宜変更されるべきものであり、発明の範囲を限定するものではない。
まず、第1および第2の実施形態にて、装置の構成と、大まかな処理の概要を説明したうえで、実施例1〜5にて、具体的な処理の内容を説明する。
第1の実施形態に係る光音響装置は、パルス光を被検体に照射し、当該パルス光に起因して被検体内で発生した光音響波を受信および解析することで、被検体内の構造情報と機能情報を可視化、すなわち画像化する装置である。構造情報とは、初期音圧分布や光吸収エネルギー密度分布、あるいは、そこから導かれる吸収係数分布に関連した被検体情報を指し、主に被検体内の光吸収体構造情報、特に、血管の構造情報である。また、機能情報とは、複数の波長で取得された光音響信号とスペクトル情報を用いて算出された分光情報である。また、機能情報とは、主に被検体内の物質濃度、特に、血管内の血液中に含まれる酸素濃度や、脂肪、コラーゲン、ヘモグロビンの濃度などの、生体機能に関する情報である。
図1を参照しながら、本実施形態に係る光音響測定装置の構成を説明する。本実施形態に係る光音響測定装置は、光照射部1、保持板21および22、音響波受信部4、信号処理部5、計算処理部6、表示部7を有している。また、計算処理部6は、第一光学特性取得部61、第二光学特性取得部62、酸素飽和度算出部63を有している。
以下、本実施形態に係る光音響測定装置を構成する各手段を説明しながら、測定の方法について概要を説明する。
光照射部1は、パルス光を発生させ、被検体に照射する手段であり、光源と、照射光学系からなる。
光源は、大出力を得るためレーザ光源であることが望ましいが、レーザの代わりに発光ダイオードやフラッシュランプ等を用いることもできる。光源としてレーザを用いる場合、固体レーザ、ガスレーザ、色素レーザ、半導体レーザなど様々なものが使用できる。
理想的には、出力が強く連続的に波長を変えられる、Nd:YAG励起のOPOレーザ、色素レーザ、Ti:saレーザや、アレキサンドライトレーザを用いるとよい。また、
異なる波長の単波長レーザを複数有していてもよい。
パルス光の波長は、被検体を構成する成分のうち特定の成分に吸収される特定の波長であって、被検体内部まで光が伝搬する波長であることが望ましい。具体的には、被検体が生体である場合、500nm以上1200nm以下であることが望ましい。
また、光音響波を効果的に発生させるためには、被検体の熱特性に応じて十分短い時間に光を照射させなければならない。被検体が生体である場合、光源から発生するパルス光のパルス幅は5ナノから50ナノ秒程度が好適である。光源から発生するパルス光を以下、照射光と称する。
なお、光源は、必ずしも本実施形態に係る光音響測定装置の一部である必要はなく、外部に接続されていてもよい。
なお、被検体表面において照射光学系を移動させるための走査機構を設けてもよい。ま
た、光源から直接、所望の形状のパルス光を照射可能である場合、照射光学系は必ずしも用いる必要はない。
保持板21および22は、被検体を保持する手段である。具体的には、二枚の平板状の保持部材のうち、いずれか片方、または両方が、X軸方向に移動することで被検体を圧迫保持する。
光照射部1から出射した計測光は、保持板21を介して被検体表面に照射されるため、保持板21は、計測光に対する透過率が高く、減衰率が低い素材であることが好ましい。典型的にはガラス、ポリメチルペンテン、ポリカーボネート、アクリルなどが好適である。
また、被検体内で発生した音響波は、保持板22を介して音響波受信部4に入射する。そのため、保持板22は、被検体との界面において音響波を反射させず、かつ、音響波を伝達しやすい素材であることが望ましい。具体的には、被検体との音響インピーダンスの差が小さい材質であることが好ましい。このような材質には、ポリメチルペンテンなどの樹脂材料がある。
被検体31および光吸収体32は、装置を構成するものではないが、ここで説明する。本実施形態に係る光音響測定装置は、血管の造影、人や動物の悪性腫瘍や血管疾患などの診断や化学治療の経過観察などを主な目的とする。よって、被検体31としては生体、具体的には人体や動物の乳房や指、手足などの診断の対象部位が想定される。なお、被検体が小動物である場合は、特定の部位だけではなく、全体を対象としてもよい。
被検体31の内部に存在する光吸収体32は、被検体内で相対的に光に対する吸収係数が高い部位である。例えば、測定の対象が人体である場合、オキシヘモグロビンまたはデオキシヘモグロビン、赤血球を含む血管、新生血管を含む悪性腫瘍などが光吸収体32に該当する。また、被検体31の表面にあるメラニンなども光吸収体32に該当する。なお、光吸収体32は、メチレンブルー(MB)、インドシニアングリーン(ICG)などの色素や金微粒子、および、それらを集積あるいは化学的に修飾した物質であってもよい。
音響波受信部4は、被検体内部で発生した音響波を受信し、電気信号に変換する手段である。音響波受信部は、音響波検出器、トランスデューサとも呼ばれる。なお、本発明における音響波とは、典型的には超音波であり、音波、超音波、光音響波、光超音波と呼ばれる弾性波を含む。
生体から発生する音響波は、100KHzから100MHzの超音波であるため、音響波受信部4には、上記の周波数帯を受信できる音響素子を用いる。具体的には、圧電現象を用いたトランスデューサ、光の共振を用いたトランスデューサ、容量の変化を用いたトランスデューサなどを用いることができる。また、音響波受信部4は、感度が高く、周波数帯域が広いものが望ましい。
また、音響波受信部4は、複数の音響波受信素子が1次元、または2次元状に配列されたものであってもよい。複数の位置で同時に音響波を受信することで、測定時間を短縮することができ、被検体の振動などの影響を低減することができる。また、音響波受信部4は、走査機構によって機械的に走査可能に構成されていてもよい。また、音響波受信部4は音響レンズを有していてもよい。
信号処理部5は、取得した電気信号(以下、光音響信号)を増幅してデジタル信号に変換する手段である。信号処理部5は、典型的には、増幅器、A/D変換器、FPGA(Field Programmable Gate Array)チップなどで構成される。音響波受信部4から得られる
信号が複数である場合は、同時に複数の信号を処理できることが望ましい。なお、本明細書における光音響信号とは、音響波受信部4で得られたアナログの電気信号と、信号処理機構によって変換されたデジタル信号の双方を含む概念である。
また、被検体に対して同じ位置で受信した光音響信号を積算し、一つの信号にしてもよい。積算の方法は、信号同士を足し合わせるものであってもよいし、平均を取るものであってもよい。また、信号にそれぞれ重みを付けて足しあわせるものであってもよい。
計算処理部6は、被検体に照射する光の強度や、光の照射タイミング、音響波の受信タイミング、受信した音響波に基づいて被検体内の情報を取得する処理等を制御する手段である。また、取得した光音響信号に基づいて、光量分布の算出や、画像再構成などを行うことにより、被検体内部の光学特性に関連した情報を取得する手段である。計算処理部6が、本発明における第一の情報取得手段、第二の情報取得手段、および、画像生成手段である。
計算処理部6は、典型的にはワークステーションであり、あらかじめ記憶されたソフトウェアによって前述した処理が行われる。
本実施形態では、計算処理部6が、第一光学特性取得部61、第二光学特性取得部62、酸素飽和度算出部63を有しており、これらの手段はソフトウェアによって実装される。
なお、本実施形態では、計算処理部6をワークステーションとしたが、計算処理部6は、複数のハードウェアの集合であってもよい。この場合、それぞれのハードウェアを総じて計算処理部6と呼ぶ。
表示部7は、計算処理部6が生成した情報を画像として表示する装置であり、典型的には液晶ディスプレイなどが利用されるが、プラズマディスプレイ、有機ELディスプレイ、FEDなど、他の方式のディスプレイであってもよい。
なお、表示部7は、必ずしも本実施形態に係る光音響測定装置の一部である必要はなく、外部に接続されていてもよい。
次に、本実施形態に係る光音響測定装置によって、被検体である生体を測定する方法について、処理フローチャートである図2を参照しながら説明する。
まず、光照射部1から発せられた特定波長のパルス光を被検体に照射する。パルス光は、例えばレンズ、ミラー、光ファイバ、拡散板などからなる照射光学系によって所望の形状に加工されながら被検体の表面に導かれ、照射される。被検体31の内部を伝播した光のエネルギーの一部が血管などの光吸収体32に吸収されると、当該光吸収体32の熱膨張により光音響波(典型的には超音波)が発生する。
発生した音響波は、被検体内を伝播し、保持板22を介して音響波受信部4で受信され、光音響信号に変換される(S1)。
本実施形態に係る光音響測定装置は、これらの測定を、パルス光の波長を変えて複数回行う。各波長に対応する光音響信号は、計算処理部6によって一時的に記憶される。
P0=Γ・μa・Φ・・・(式1)
ここで、Γはグルナイゼン係数であり、体積膨張係数βと音速cの二乗の積を定圧比熱CPで割ったものである。Γは被検体が決まれば、ほぼ一定の値をとることが知られてい
る。また、μaは吸収体の光吸収係数であり、Φは局所的な領域での光量(吸収体に照射された光量で、光フルエンスとも言う)である。
被検体内における酸素飽和度StOは、式2によって表すことができる。
StO=[(E_HbR(λ2)−{A(λ2)/A(λ1)}×E_HbR(λ1))]/[{E_HbR(λ2)−E_HbО(λ2)}−{A(λ2)/A(λ1)}×{E_HbR(λ1)−E_HbO(λ1)}]×100 ・・・(式2)
ここでStOは酸素飽和度、E_HbR(λ1)とE_HbR(λ2)は、それぞれ波長λ1と波長λ2に対する還元ヘモグロビンの1モル/リットルあたりの吸収係数である。また、E_HbO(λ1)とE_HbO(λ2)は、同様に酸化ヘモグロビンの1モル/リットルあたりの吸収係数である。また、A(λ1)は、波長λ1に対応する被検体内の吸収係数であり、A(λ2)は、波長λ2に対応する被検体内の吸収係数である。
そして、酸素飽和度算出部63が、当該第一の吸収係数分布に基づいて、酸素飽和度分布を算出する(S3)。
また、第二光学特性取得部62が、記憶された光音響信号に基づいて、第二の計算手法によって吸収係数分布を算出する(S4)。第二の計算手法とは、算出した吸収係数分布を画像化した際に、視認性が高くなるような手法である。すなわち、血管構造の提示に適した算出手法である。このようにして求めた吸収係数分布を、第二の吸収係数分布と称する。
第一光学特性取得部61は、初期音圧P(λ1)およびP(λ2)に対して、第一の計算手法を適用した処理を行い、第一の吸収係数分布B(λ1)、B(λ2)を取得する。
第二光学特性取得部62は、初期音圧P(λ1)およびP(λ2)に対して、第二の計算手法を適用した処理を行い、第二の吸収係数分布C(λ1)、C(λ2)を取得する。
第2の実施形態は、半球状の音響波受信部を用い、人の乳房に対して測定を行う光音響測定装置である。図3は、第2の実施形態に係る光音響測定装置のシステム構成図である。
光源11から出射されたパルス光は、ミラーやビームエキスパンダーなどの光学部材からなる光学系12を用いて、ある程度半径を広げた後に、被検体に照射される。
127mmの半球であり、複数の音響波受信素子42は、全て当該半球の曲率中心を向くように配置される。音響波受信素子42のサイズは2mmで、帯域は2MHz100%の
cMUT素子である。
また、光学系12から出射したパルス光は、支持体41の底部から、被検体に向けてY軸正方向に照射される。
また、支持体41は、移動部43によって、図中点線で示したY軸を中心として、X−Z平面状で回転可能な構成となっている。すなわち、支持体41を回転させることで、支持体41上に配置された複数の音響波受信素子が、被検体に対して移動可能となっている。このようにすることで、被検体に対して複数の位置で音響波を受信することが可能になり、測定の精度を向上させることができる。
計算処理部6は、第1の実施形態と同様に、第一光学特性取得部61、第二光学特性取得部62、酸素飽和度算出部63を有している。これらの手段は、第一の実施形態と同様に、ワークステーションで動作するソフトウェアとして実装される。
第2の実施形態では、信号処理部5が取得した光音響信号が、計算処理部6内のメモリに記憶され、ソフトウェアによって処理される。以降の処理は、第1の実施形態と同様である。
実施例1では、第一光学特性取得部61が、再構成視野角を10度として、UBP法によって初期音圧分布を算出したうえで、当該初期音圧分布に基づいて吸収係数分布を算出する。当該手法が、実施例1における第一の計算手法である。
また、第二光学特性取得部62が、再構成視野角を15度として、UBP法によって初期音圧分布を算出したうえで、当該初期音圧分布に基づいて吸収係数分布を算出する。当該手法が、実施例1における第二の計算手法である。
なお、再構成視野角とは、ある素子で受信された信号を再構成領域にバックプロジェクションする際に処理対象となる、素子の感度が最大となる方向に対する最大角度である。
図4(A)は、波長756nmのパルス光を被検体に照射して得られた吸収係数分布であり、図4(B)は、波長797nmのパルス光を被検体に照射して得られた吸収係数分
布である。被検体内に吸収体は2つあり、左にある吸収体401が、動脈に相当する、酸素飽和度96%の吸収体である。当該吸収体の吸収係数は、波長756nmにおいて0.138/mm、波長797nmにおいて0.189/mmである。
また、右にある吸収体402が、静脈に相当する、酸素飽和度76%の吸収体である。当該吸収体の吸収係数は、波長756nmにおいて0.185/mm、波長797nmにおいて0.189/mmである。
一方、右側の吸収体は、支持体の曲率中心点から約22mmの半径内に入っていないため、共に再構成アーティファクトが発生し、形状が精度よく再現できていないうえ、コントラストも悪くなっている。
すなわち、第二の計算手法では、より大きい再構成視野角を採用することが好ましいことがわかる。このようにすることで、吸収係数分布を画像化した際に、より高い視認性を得ることができる。
図5(A)は、図4(C)および図4(D)に示した吸収係数分布に基づいて算出された酸素飽和度分布であり、図5(B)は、図4(E)および図4(F)に示した吸収係数分布に基づいて算出された酸素飽和度分布である。
また、図5(C)は、図5(A)の左側にある吸収体に対応する、酸素飽和度のヒストグラムであり、図5(E)は、図5(A)の右側にある吸収体に対応する、酸素飽和度のヒストグラムである。
同様に、図5(D)は、図5(B)の左側にある吸収体に対応する、酸素飽和度のヒストグラムであり、図5(F)は、図5(B)の右側にある吸収体に対応する、酸素飽和度のヒストグラムである。
図5(C)を参照すると、平均値が96.21%、分散が0.00926%という結果が得られており、図5(E)を参照すると、平均値が77.63%、分散が0.0285%という結果が得られている。
この結果を、図5(D)および図5(F)と比較すると、図5(C)および図5(E)、すなわち、再構成視野角が10度である場合のほうが、酸素飽和度が真値に近く、分散も小さいことがわかる。これは、再構成視野角を小さくすることにより、他の吸収体のストリークアーティファクトの影響を受けにくくなり、酸素飽和度の誤差が小さくなるためである。
すなわち、第一の計算手法では、より小さい再構成視野角を採用することが好ましいことがわかる。このようにすることで、酸素飽和度分布を画像化した際に、より正確な結果を得ることができる。
実施例2は、第一の計算手法に、ホワイトノイズを軽減するための平均化処理を含ませた実施例である。実施例2に係る光音響測定装置の構成は、以下に説明する点を除き、実施例1と同様である。
また、実施例2において、移動部43は、支持体41をX−Z方向に平行移動させるための手段である。
また、第二光学特性取得部62が、第一の吸収係数分布算出機構と同様にUBP法を用いて初期音圧分布を算出し、光量分布を除することによって吸収係数分布を算出する。こちらは、前述したような加算平均を行わない。
このように、実施例2では、周囲のボクセル同士で加算平均を行う処理を、第一の計算手法にのみ含め、第二の計算手法には含めない。
図6(A)は、波長756nmのパルス光を照射して得られた吸収係数分布であり、図6(B)は、波長797nmのパルス光を照射して得られた吸収係数分布である。
各吸収係数分布は、吸収係数が0.138/mmおよび0.189/mm、酸素飽和度が96%である、Φ2mmの動脈血管相当の吸収係数分布である。
また、当該吸収係数分布は、分散値が20である、正規分布をもったホワイトノイズを仮想的な吸収係数として与えたものである。図6に示した各図の大きさは、縦5mm×横5mmであり、ボクセルサイズは0.1mmである。
また、図8(C)が、図7(C)および図7(D)に示した吸収係数分布を用いて算出した酸素飽和度分布であり、図8(D)が、対応するヒストグラムである。
実施例3は、第二の計算手法に、音響波の多重反射によって被検体の内部に結像されるアーティファクトを取り除くための信号処理を含ませた実施例である。
実施例3に係る光音響測定装置は、第1の実施形態に対応する実施例である。すなわち、平板状の保持板を用いて被検体を圧迫保持し、当該保持板越しにパルス光を照射し、二次元探触子を用いて音響波を取得することで測定を行う装置である。
そこで、実施例3では、第二光学特性取得部62が、音響波の多重反射に起因して被検体内部に発生するアーティファクトを取り除くための信号処理を行う。一方、第一光学特性取得部61は、このような信号処理を行わない。
このように、実施例3では、アーティファクトを除去する処理を、第二の計算手法にのみ含め、第一の計算手法には含めない。
そのため、XY方向に、DC成分を含む低周波成分を消すためのフィルタ(ハイパスフィルタ)をかけることで、これらの平面波アーティファクトを消すことができ、血管イメージング精度を上げることができる。このような手法は、特許文献1に記載されている。
間の吸収係数の比が崩れるおそれがある。
そこで、実施例3では、第一光学特性取得部61が、このような平面波アーティファクト削除処理を行わずに吸収係数分布を算出し、かつ、第二光学特性取得部62が、平面波アーティファクト削除処理を行ったうえで吸収係数分布を算出する。
実施例4は、第二光学特性取得部が、信号のインパルス応答補正処理、ブラインドデコンボリューション処理、空間インパルス応答補正処理等の最適化処理、または、モデルベース再構成法等の再構成処理を行うことによって、吸収係数分布を生成する実施例である。
なお、最適化処理は、信号のインパルス応答補正、ブラインドデコンボリューション、空間インパルス応答補正などの追加処理であってもよく、モデルベース再構成など、再構成の手法そのものであってもよい。
実施例5は、第二光学特性取得部が、線状の物体の形状を強調する画像処理(以下、線強調処理)を行う実施例である。線強調処理により、血管部分を強調することができる。
そこで、実施例5では、第二光学特性取得部62のみが、算出された吸収係数分布に対して線強調処理を行い、第一光学特性取得部61は、通常の方法によって吸収係数分布を算出する。
なお、各実施形態の説明や、実施例の説明は本発明を説明する上での例示であり、本発明は、発明の趣旨を逸脱しない範囲で適宜変更または組み合わせて実施することができる。
例えば、本発明は、上記処理の少なくとも一部を含む光音響装置として実施することもできる。また、上記処理の少なくとも一部を含む光音響装置が実行する処理方法として実施することもできる。上記処理や手段は、技術的な矛盾が生じない限りにおいて、自由に組み合わせて実施することができる。
Claims (13)
- 被検体に対して複数の波長のパルス光を照射可能な光源と、
前記パルス光を照射された被検体から発生する音響波を受信し、電気信号に変換する音響波受信部と、
前記電気信号に基づいて、第一の計算手法を用いて、前記被検体内における第一の吸収係数分布を取得する第一の情報取得手段と、
前記電気信号に基づいて、第二の計算手法を用いて、前記被検体内における第二の吸収係数分布を取得する第二の情報取得手段と、
複数の波長のパルス光をそれぞれ照射して得られた、複数の前記第一の吸収係数分布に基づいて、前記被検体内の機能情報の分布を算出する第三の情報取得手段と、
前記機能情報の分布を、前記第二の吸収係数分布に基づいてマスクした画像を生成する画像生成手段と、を有し、
前記第二の計算手法は、前記第一の計算手法と比べて、吸収係数分布を画像化した際により高い視認性が得られる手法である
ことを特徴とする、光音響装置。 - 前記第二の計算手法は、アーティファクトを低減する処理、または、解像度を向上させる処理を含む
ことを特徴とする、請求項1に記載の光音響装置。 - 前記第二の計算手法は、前記電気信号に対してインパルス応答補正を行う処理を含む
ことを特徴とする、請求項2に記載の光音響装置。 - 前記第二の計算手法は、空間インパルス応答補正を行う処理を含む
ことを特徴とする、請求項2に記載の光音響装置。 - 前記第二の計算手法は、ブラインドデコンボリューション処理を含む
ことを特徴とする、請求項2に記載の光音響装置。 - 前記第二の計算手法は、吸収係数の算出における再構成視野角を、前記第一の計算手法に比べて大きくしたものである
ことを特徴とする、請求項2に記載の光音響装置。 - 前記第二の計算手法は、モデルベース再構成法によって吸収係数を算出するものである
ことを特徴とする、請求項2に記載の光音響装置。 - 前記第一の計算手法は、前記第二の計算手法と比べて、波長間における吸収係数の比をより精度よく算出できる手法である
ことを特徴とする、請求項1から7のいずれか1項に記載の光音響装置。 - 前記第一の計算手法は、吸収係数の算出における再構成視野角を、前記第二の計算手法に比べて小さくしたものである
ことを特徴とする、請求項8に記載の光音響装置。 - 前記第一の計算手法は、周囲のピクセルまたはボクセル同士で加算平均を行う処理を含む
ことを特徴とする、請求項8に記載の光音響装置。 - 前記画像生成手段は、前記機能情報の分布が表す値を色相として、前記第二の吸収係数
分布が表す値を明度としてそれぞれ表した画像を生成する
ことを特徴とする、請求項1から10のいずれか1項に記載の光音響装置。 - 前記第三の情報取得手段は、前記機能情報の分布として、酸素飽和度分布を取得する
ことを特徴とする、請求項1から11のいずれか1項に記載の光音響装置。 - 被検体に対して複数の波長のパルス光を照射可能な光源と、前記パルス光を照射された被検体から発生する音響波を受信し、電気信号に変換する音響波受信部と、を有する光音響装置の処理方法であって、
前記電気信号に基づいて、第一の計算手法を用いて、前記被検体内における第一の吸収係数分布を取得する第一の情報取得ステップと、
前記電気信号に基づいて、第二の計算手法を用いて、前記被検体内における第二の吸収係数分布を取得する第二の情報取得ステップと、
複数の波長のパルス光をそれぞれ照射して得られた、複数の前記第一の吸収係数分布に基づいて、前記被検体内の機能情報の分布を算出する第三の情報取得ステップと、
前記機能情報の分布を、前記第二の吸収係数分布に基づいてマスクした画像を生成する画像生成ステップと、を含み、
前記第二の計算手法は、前記第一の計算手法と比べて、吸収係数分布を画像化した際により高い視認性が得られる手法である
ことを特徴とする、光音響装置の処理方法。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2014209193A JP6486056B2 (ja) | 2014-10-10 | 2014-10-10 | 光音響装置および光音響装置の処理方法 |
PCT/JP2015/005098 WO2016056237A1 (en) | 2014-10-10 | 2015-10-07 | Photoacoustic apparatus and processing method for photoacoustic apparatus |
US15/510,743 US20170273568A1 (en) | 2014-10-10 | 2015-10-07 | Photoacoustic apparatus and processing method for photoacoustic apparatus |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2014209193A JP6486056B2 (ja) | 2014-10-10 | 2014-10-10 | 光音響装置および光音響装置の処理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2016077352A true JP2016077352A (ja) | 2016-05-16 |
JP6486056B2 JP6486056B2 (ja) | 2019-03-20 |
Family
ID=54361133
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2014209193A Expired - Fee Related JP6486056B2 (ja) | 2014-10-10 | 2014-10-10 | 光音響装置および光音響装置の処理方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US20170273568A1 (ja) |
JP (1) | JP6486056B2 (ja) |
WO (1) | WO2016056237A1 (ja) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2018126389A (ja) * | 2017-02-09 | 2018-08-16 | キヤノン株式会社 | 情報処理装置、情報処理方法、およびプログラム |
US10544811B2 (en) * | 2017-02-21 | 2020-01-28 | University Of Electronic Science And Technology Of China | Photoacoustic layer disposed on a substrate generating directional ultrasound waves |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012217554A (ja) * | 2011-04-06 | 2012-11-12 | Canon Inc | 光音響装置およびその制御方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6053265B2 (ja) * | 2011-09-01 | 2016-12-27 | キヤノン株式会社 | 情報処理装置および情報処理方法 |
-
2014
- 2014-10-10 JP JP2014209193A patent/JP6486056B2/ja not_active Expired - Fee Related
-
2015
- 2015-10-07 US US15/510,743 patent/US20170273568A1/en not_active Abandoned
- 2015-10-07 WO PCT/JP2015/005098 patent/WO2016056237A1/en active Application Filing
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012217554A (ja) * | 2011-04-06 | 2012-11-12 | Canon Inc | 光音響装置およびその制御方法 |
Also Published As
Publication number | Publication date |
---|---|
US20170273568A1 (en) | 2017-09-28 |
JP6486056B2 (ja) | 2019-03-20 |
WO2016056237A1 (en) | 2016-04-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10709419B2 (en) | Dual modality imaging system for coregistered functional and anatomical mapping | |
AU2012332233B2 (en) | Dual modality imaging system for coregistered functional and anatomical mapping | |
US9757092B2 (en) | Method for dual modality optoacoustic imaging | |
US20100087733A1 (en) | Biological information processing apparatus and biological information processing method | |
JP6132466B2 (ja) | 被検体情報取得装置及び被検体情報取得方法 | |
JP2012024460A (ja) | 画像情報取得装置及びその制御方法 | |
JP2013255697A (ja) | 被検体情報取得装置およびその制御方法 | |
JP2011217914A (ja) | 光音響イメージング装置、光音響イメージング方法およびプログラム | |
JP2017042612A (ja) | 被検体情報を取得する装置および方法 | |
JP6222936B2 (ja) | 装置および画像生成方法 | |
JP2017029610A (ja) | 光音響装置、信頼度取得方法、プログラム | |
JP6486056B2 (ja) | 光音響装置および光音響装置の処理方法 | |
JP6664176B2 (ja) | 光音響装置、情報処理方法、およびプログラム | |
JP6598528B2 (ja) | 被検体情報取得装置および被検体情報取得方法 | |
JP2018126454A (ja) | 被検体情報取得装置および表示方法 | |
JP6521761B2 (ja) | 情報処理装置および表示制御方法 | |
JP2013188489A (ja) | 被検体情報処理装置およびその作動方法 | |
JP2019083887A (ja) | 情報処理装置および情報処理方法 | |
JP6513121B2 (ja) | 処理装置、被検体情報取得装置、光音響画像の表示方法、及びプログラム | |
JP2017042603A (ja) | 被検体情報取得装置 | |
JP6643108B2 (ja) | 被検体情報取得装置および被検体情報取得方法 | |
JP2018161467A (ja) | 画像処理装置および画像処理方法 | |
EP2773267B1 (en) | Dual modality imaging system for coregistered functional and anatomical mapping | |
JP2019136520A (ja) | 処理装置、光音響画像の表示方法、及びプログラム | |
JP2019155004A (ja) | 光音響装置および被検体情報取得方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20171006 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20180911 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20180918 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20181115 |
|
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: 20190122 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20190219 |
|
R151 | Written notification of patent or utility model registration |
Ref document number: 6486056 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |
|
LAPS | Cancellation because of no payment of annual fees |