JP2018061716A - 情報処理装置、情報処理方法、及びプログラム - Google Patents
情報処理装置、情報処理方法、及びプログラム Download PDFInfo
- Publication number
- JP2018061716A JP2018061716A JP2016202044A JP2016202044A JP2018061716A JP 2018061716 A JP2018061716 A JP 2018061716A JP 2016202044 A JP2016202044 A JP 2016202044A JP 2016202044 A JP2016202044 A JP 2016202044A JP 2018061716 A JP2018061716 A JP 2018061716A
- Authority
- JP
- Japan
- Prior art keywords
- information
- sound speed
- subject
- correction amount
- unit
- 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.)
- Withdrawn
Links
Images
Classifications
-
- 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
-
- 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
-
- 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/47—Scattering, i.e. diffuse reflection
- G01N21/4795—Scattering, i.e. diffuse reflection spatially resolved investigating of object in scattering medium
-
- 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/02—Analysing fluids
- G01N29/024—Analysing fluids by measuring propagation velocity or propagation time of acoustic waves
-
- 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
-
- 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/07—Analysing solids by measuring propagation velocity or propagation time of acoustic waves
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2562/00—Details of sensors; Constructional details of sensor housings or probes; Accessories for sensors
- A61B2562/04—Arrangements of multiple sensors of the same type
- A61B2562/046—Arrangements of multiple sensors of the same type in a matrix array
-
- 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
- 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
- G01N2021/178—Methods for obtaining spatial resolution of the property being measured
- G01N2021/1785—Three dimensional
-
- 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/011—Velocity or travel time
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Pathology (AREA)
- General Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- Chemical & Material Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Acoustics & Sound (AREA)
- Heart & Thoracic Surgery (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Optics & Photonics (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
【課題】本発明は、少ない計算量で音響波の屈折によって生じる影響が抑制された被検体情報を取得する情報処理装置を提供する。【解決手段】情報処理装置は、被検体内の注目位置から伝搬する音響波が複数の受信手段で受信されることにより得られる複数の信号を用いて、注目位置の被検体情報を取得する情報処理装置であって、音響波の伝搬経路上の媒質の音速情報を取得する音速取得手段と、音速情報に対する補正量を取得し、当該補正量を用いて、音速情報を補正する補正手段と、補正された音速情報に基づいて、注目位置から複数の受信手段までの伝搬経路を直線近似して音響波の伝搬時間を決定し、複数の信号及び伝搬時間に基づいて、注目位置の被検体情報を取得する情報取得手段と、を有する。【選択図】図1
Description
本発明は、音響波の受信信号に基づいて被検体情報を取得する情報処理装置、情報処理方法に関する。
被検体から伝搬する音響波の受信信号に基づいてイメージングを行う装置が知られている。特許文献1には、光音響効果により発生した音響波(光音響波とも呼ぶ)の受信信号を用いて、被検体内の情報を取得する光音響装置が記載されている。
また、被検体内で発生した音響波が、被検体とは音速値が異なる媒質を伝搬するときに屈折することが知られている。この屈折によって生じる音響波の伝搬時間への影響を考慮せずに再構成を行うと、得られる画像の画質の低下を招いてしまう。
そこで、特許文献1では、屈折した音響波の伝搬経路(音線とも呼ぶ)をスネルの法則に基づいて算出し、この伝搬経路に沿った音響波の伝搬時間を算出している。また、特許文献1では、屈折を考慮して算出された伝搬時間に基づいて再構成を行い、光音響イメージングによる画像化を行っている。
ところが、特許文献1に記載されたように、屈折した音響波の伝搬経路をスネルの法則に基づいて算出する場合、計算量が膨大となり、画像化に多くの時間を要してしまう。
そこで、本発明は、少ない計算量で音響波の屈折によって生じる影響が抑制された被検体情報を取得する情報処理装置を提供することを目的とする。
本発明に係る情報処理装置は、被検体内の注目位置から伝搬する音響波が複数の受信手段で受信されることにより得られる複数の信号を用いて、注目位置の被検体情報を取得する情報処理装置であって、音響波の伝搬経路上の媒質の音速情報を取得する音速取得手段と、音速情報に対する補正量を取得し、当該補正量を用いて、音速情報を補正する補正手段と、補正された音速情報に基づいて、注目位置から複数の受信手段までの伝搬経路を直線近似して音響波の伝搬時間を決定し、複数の信号及び伝搬時間に基づいて、注目位置の被検体情報を取得する情報取得手段と、を有する。
本発明に係る情報処理装置によれば、少ない計算量で音響波の屈折によって生じる影響が抑制された被検体情報を取得することができる。
以下に図面を参照しつつ、本発明の好適な実施の形態について説明する。ただし、以下に記載されている構成部品の寸法、材質、形状およびそれらの相対配置などは、発明が適用される装置の構成や各種条件により適宜変更されるべきものである。よって、この発明の範囲を以下の記載に限定する趣旨のものではない。
本発明に係る音響波とは、典型的には超音波であり、音波、音響波と呼ばれるものを含み、光音響効果により発生した音響波である光音響波や、被検体内部の超音波を送信して反射された超音波を含む。本発明は、光音響効果により発生した音響波に基づいて画像データを取得する光音響装置や、被検体に超音波を送受信することにより画像データを取得する超音波装置に適用することができる。
本発明に係る光音響装置により得られる光音響画像は、光照射により発生した光音響波に由来するあらゆる画像を含む概念である。光音響画像は、光音響波の発生音圧(初期音圧)、光吸収エネルギー密度、及び光吸収係数、被検体を構成する物質の濃度(酸素飽和度など)などの少なくとも1つの被検体情報の空間分布を表す画像データである。また、本発明に係る超音波装置により得られる超音波画像は、超音波の送信により発生したエコー波に由来するあらゆる画像を含む概念である。超音波画像は、Bモード画像、ドップラー画像、エラスト画像の少なくとも一つの画像を含む。これらの超音波画像は、被検体の音響インピーダンス、血流、弾性などの被検体の特性に基づいた被検体情報の空間分布を表す画像データである。
以下、本実施形態の光音響装置の構成及び情報処理方法について説明する。
本実施形態において、光音響装置を用いた例を説明する。図1を用いて本実施形態に係る光音響装置の構成を説明する。図1は、光音響装置全体の概略ブロック図である。本実施形態に係る光音響装置は、駆動部130、信号収集部140、コンピュータ150、表示部160、入力部170、プローブ180、及び温度センサ190を有する。プローブ180は、光照射部110、及び受信部120を有する。図2は、本実施形態に係るプローブ180の模式図を示す。測定対象は、被検体100である。駆動部130は、光照射部110と受信部120を駆動し、機械的な走査を行う。光照射部110が光を被検体100に照射し、被検体100内で音響波が発生する。光に起因して光音響効果により発生する音響波を光音響波とも呼ぶ。受信部120は、光音響波を受信することによりアナログ信号としての電気信号(光音響信号)を出力する。
信号収集部140は、受信部120から出力されたアナログ信号をデジタル信号に変換し、コンピュータ150に出力する。コンピュータ150は、信号収集部140から出力されたデジタル信号を、超音波または光音響波に由来する信号データとして記憶する。
コンピュータ150は、記憶されたデジタル信号に対して信号処理を行うことにより、光音響画像を表す画像データを生成する。また、コンピュータ150は、得られた画像データに対して画像処理を施した後に、画像データを表示部160に出力する。表示部160は、光音響画像を表示する。ユーザーとしての医師や技師等は、表示部160に表示された光音響画像を確認することにより、診断を行うことができる。表示画像は、ユーザーやコンピュータ150からの保存指示に基づいて、コンピュータ150内のメモリや、モダリティとネットワークで接続されたデータ管理システムなどに保存される。
また、コンピュータ150は、光音響装置に含まれる構成の駆動制御も行う。また、表示部160は、コンピュータ150で生成された画像の他にGUIなどを表示してもよい。入力部170は、ユーザーが情報を入力できるように構成されている。ユーザーは、入力部170を用いて測定開始や終了、作成画像の保存指示などの操作を行うことができる。
本実施形態に係る情報処理方法では、画像データの生成に用いるために取得された音速情報を所定の補正量で補正し、補正された音速情報を用いて、音響波の伝搬経路を直線近似して再構成する。これにより、測定対象である被検体と音響マッチング材との音速差に起因して屈折した音響波の伝搬経路を算出することなく、少ない計算量で音響波の屈折によって生じる画質への影響を抑制することができる。
以下、本実施形態に係る光音響装置の各構成の詳細を説明する。
(光照射部110)
光照射部110は、光を発する光源111と、光源111から射出された光を被検体100へ導く光学系112とを含む。なお、光は、いわゆる矩形波、三角波などのパルス光を含む。
光照射部110は、光を発する光源111と、光源111から射出された光を被検体100へ導く光学系112とを含む。なお、光は、いわゆる矩形波、三角波などのパルス光を含む。
光源111が発する光のパルス幅としては、1ns以上、100ns以下のパルス幅であってもよい。また、光の波長として400nmから1600nm程度の範囲の波長であってもよい。血管を高解像度でイメージングする場合は、血管での吸収が大きい波長(400nm以上、700nm以下)を用いてもよい。生体の深部をイメージングする場合には、生体の背景組織(水や脂肪など)において典型的に吸収が少ない波長(700nm以上、1100nm以下)の光を用いてもよい。
光源111としては、レーザーや発光ダイオードを用いることができる。また、複数波長の光を用いて測定する際には、波長の変更が可能な光源であってもよい。なお、複数波長を被検体に照射する場合、互いに異なる波長の光を発生する複数台の光源を用意し、それぞれの光源から交互に照射することも可能である。複数台の光源を用いた場合もそれらをまとめて光源として表現する。レーザーとしては、固体レーザー、ガスレーザー、色素レーザー、半導体レーザーなど様々なレーザーを使用することができる。例えば、Nd:YAGレーザーやアレキサンドライトレーザーなどのパルスレーザーを光源として用いてもよい。また、Nd:YAGレーザー光を励起光とするTi:saレーザーやOPO(Optical Parametric Oscillators)レーザーを光源として用いてもよい。また、光源111としてフラッシュランプや発光ダイオードを用いてもよい。また、光源111としてマイクロウェーブ源を用いてもよい。
光学系112には、レンズ、ミラー、光ファイバ等の光学素子を用いることができる。乳房等を被検体100とする場合、パルス光のビーム径を広げて照射するために、光学系の光出射部は光を拡散させる拡散板等で構成されていてもよい。一方、光音響顕微鏡においては、解像度を上げるために、光学系112の光出射部はレンズ等で構成し、ビームをフォーカスして照射してもよい。
なお、光照射部110が光学系112を備えずに、光源111から直接被検体100に光を照射してもよい。
(受信部120)
受信部120は、音響波を受信することにより電気信号を出力するトランスデューサ121と、トランスデューサ121を支持する支持体122とを含む。また、トランスデューサ121は、音響波を送信する送信手段としてもよい。受信手段としてのトランスデューサと送信手段としてのトランスデューサとは、単一(共通)のトランスデューサでもよいし、別々の構成であってもよい。
受信部120は、音響波を受信することにより電気信号を出力するトランスデューサ121と、トランスデューサ121を支持する支持体122とを含む。また、トランスデューサ121は、音響波を送信する送信手段としてもよい。受信手段としてのトランスデューサと送信手段としてのトランスデューサとは、単一(共通)のトランスデューサでもよいし、別々の構成であってもよい。
トランスデューサ121を構成する部材としては、PZT(チタン酸ジルコン酸鉛)に代表される圧電セラミック材料や、PVDF(ポリフッ化ビニリデン)に代表される高分子圧電膜材料などを用いることができる。また、圧電素子以外の素子を用いてもよい。例えば、静電容量型トランスデューサ(CMUT:Capacitive Micro−machined Ultrasonic Transducers)、ファブリペロー干渉計を用いたトランスデューサなどを用いることができる。なお、音響波を受信することにより電気信号を出力できる限り、いかなるトランスデューサを採用してもよい。また、トランスデューサにより得られる信号は時間分解信号である。つまり、トランスデューサにより得られる信号の振幅は、各時刻にトランスデューサで受信される音圧に基づく値(例えば、音圧に比例した値)を表したものである。
光音響波を構成する周波数成分は、典型的には100KHzから100MHzであり、トランスデューサ121として、これらの周波数を検出することのできるものを採用することができる。
支持体122は、機械的強度が高い金属材料などから構成されていてもよい。照射光を被検体に多く入射させるために、支持体122の被検体100側の表面に鏡面もしくは光散乱させる加工が行われていてもよい。本実施形態において支持体122は半球殻形状であり、半球殻上に複数のトランスデューサ121を支持できるように構成されている。この場合、支持体122に配置されたトランスデューサ121の指向軸は半球の曲率中心付近に集まる。そして、複数のトランスデューサ121から出力された信号を用いて画像化したときに曲率中心付近の画質が高くなる。なお、支持体122はトランスデューサ121を支持できる限り、いかなる構成であってもよい。支持体122は、1Dアレイ、1.5Dアレイ、1.75Dアレイ、2Dアレイと呼ばれるような平面又は曲面内に、複数のトランスデューサを並べて配置してもよい。複数のトランスデューサ121が複数の受信手段に相当する。
また、支持体122は音響マッチング材を貯留する容器として機能してもよい。すなわち、支持体122をトランスデューサ121と被検体100との間に音響マッチング材を配置するための容器としてもよい。
また、受信部120が、トランスデューサ121から出力される時系列のアナログ信号を増幅する増幅器を備えてもよい。また、受信部120が、トランスデューサ121から出力される時系列のアナログ信号を時系列のデジタル信号に変換するA/D変換器を備えてもよい。すなわち、受信部120が後述する信号収集部140を備えてもよい。
なお、音響波を様々な角度で検出できるようにするために、理想的には被検体100を全周囲から囲むようにトランスデューサ121を配置してもよい。ただし、被検体100が大きく全周囲を囲むようにトランスデューサを配置できない場合は、半球状の支持体122上にトランスデューサを配置して全周囲を囲む状態に近づけてもよい。
なお、トランスデューサの配置や数及び支持体の形状は被検体に応じて最適化すればよく、本発明に関してはあらゆる受信部120を採用することができる。
受信部120と被検体100との間の空間は、光音響波が伝播することができる媒質で満たす。この媒質には、音響波が伝搬でき、被検体100やトランスデューサ121との界面において音響特性が整合し、できるだけ光音響波の透過率が高い材料を採用する。例えば、この媒質には、水、超音波ジェルなどを採用することができる。
図2(a)は、プローブ180の側面図を示し、図2(b)は、プローブ180の上面図(図2(a)の紙面上方向から見た図)を示す。図2に示された本実施形態に係るプローブ180は、開口を有する半球状の支持体122に複数のトランスデューサ121が3次元に配置された受信部120を有する。また、図2に示されたプローブ180は、支持体122の底部に光学系112の光射出部が配置されている。また、図2に示されたプローブ180は、支持体122の内部に貯留された音響マッチング材210の温度を測定する温度センサ190を有する。図2では、温度センサがプローブ180に取り付けられ、音響マッチング材210と接触して温度を測定する例を示したが、プローブ180とは別に構成された温度センサにより非接触に音響マッチング材210の温度を測定してもよい。例えば、温度センサ190としては、サーミスタ式温度センサなどの接触型の温度センサや放射型温度センサなどの非接触型の温度センサを用いることができる。
本実施形態においては、図2に示すように被検体100は、保持部200に接触することにより、その形状が保持される。本実施形態では、被検体100が乳房の場合に、伏臥位の被検者を支持する寝台に乳房を挿入するための開口を設けて、開口から鉛直方向に垂らされた乳房を測定する形態を想定している。
受信部120と保持部200の間の空間は、光音響波が伝播することができる媒質(音響マッチング材210)で満たされる。この媒質には、光音響波が伝搬でき、被検体100やトランスデューサ121との界面において音響特性が整合し、できるだけ光音響波の透過率が高い材料を採用する。例えば、この媒質には、水、超音波ジェルなどを採用することができる。
保持手段としての保持部200は被検体100の形状を測定中に保持するために使用される。保持部200により被検体100を保持することによって、被検体100の動きの抑制および被検体100の位置を保持部200内に留めることができる。保持部200の材料には、ポリカーボネートやポリエチレン、ポリエチレンテレフタレート等、樹脂材料を用いることができる
保持部200は、被検体100を保持できる硬度を有する材料であることが好ましい。保持部200は、測定に用いる光を透過する材料であってもよい。保持部200は、インピーダンスが被検体100と同程度の材料で構成されていてもよい。乳房等の曲面を有するものを被検体100とする場合、凹型に成型した保持部200であってもよい。この場合、保持部200の凹部分に被検体100を挿入することができる。
保持部200は、取り付け部201に取り付けられている。取り付け部201は、被検体の大きさに合わせて複数種類の保持部200を交換可能に構成されていてもよい。例えば、取り付け部201は、曲率半径や曲率中心などの異なる保持部に交換できるように構成されていてもよい。
また、保持部200には保持部200の情報が登録されたタグ202が設置されていてもよい。例えば、タグ202には、保持部200の曲率半径、曲率中心、音速、識別ID等の情報を登録することができる。タグ202に登録された情報は、読み取り部203により読み出され、コンピュータ150に転送される。保持部200が取り付け部201に取り付けられたときに容易にタグ202を読み取るために、読み取り部203は取り付け部201に設置されていてもよい。例えば、タグ202はバーコードであり、読み取り部203はバーコードリーダである。
(駆動部130)
駆動部130は、被検体100と受信部120との相対位置を変更する部分である。本実施形態では、駆動部130は、支持体122をXY方向に移動させる装置であり、ステッピングモータを搭載した電動のXYステージある。駆動部130は、駆動力を発生させるステッピングモータなどのモータと、駆動力を伝達させる駆動機構と、受信部120の位置情報を検出する位置センサとを含む。駆動機構としては、リードスクリュー機構、リンク機構、ギア機構、油圧機構、などを用いることができる。また、位置センサとしては、エンコーダー、可変抵抗器、などを用いたポテンショメータなどを用いることができる。
駆動部130は、被検体100と受信部120との相対位置を変更する部分である。本実施形態では、駆動部130は、支持体122をXY方向に移動させる装置であり、ステッピングモータを搭載した電動のXYステージある。駆動部130は、駆動力を発生させるステッピングモータなどのモータと、駆動力を伝達させる駆動機構と、受信部120の位置情報を検出する位置センサとを含む。駆動機構としては、リードスクリュー機構、リンク機構、ギア機構、油圧機構、などを用いることができる。また、位置センサとしては、エンコーダー、可変抵抗器、などを用いたポテンショメータなどを用いることができる。
なお、駆動部130は被検体100と受信部120との相対位置をXY方向(二次元)に変更させるものに限らず、一次元または三次元に変更させてもよい。
なお、駆動部130は、被検体100と受信部120との相対的な位置を変更できれば、受信部120を固定し、被検体100を移動させてもよい。被検体100を移動させる場合は、被検体100を保持する保持部を動かすことで被検体100を移動させる構成などが考えられる。また、被検体100と受信部120の両方を移動させてもよい。
駆動部130は、相対位置を連続的に移動させてもよいし、ステップアンドリピートによって移動させてもよい。駆動部130は、プログラムされた軌跡で移動させる電動ステージであってもよいし、手動ステージであってもよい。
また、本実施形態では、駆動部130は光照射部110と受信部120を同時に駆動して走査を行っているが、光照射部110だけを駆動したり、受信部120だけを駆動したりしてもよい。
(信号収集部140)
信号収集部140は、トランスデューサ121から出力されたアナログ信号である電気信号を増幅するアンプと、アンプから出力されたアナログ信号をデジタル信号に変換するA/D変換器とを含む。信号収集部140は、FPGA(Field Programmable Gate Array)チップなどで構成されてもよい。信号収集部140から出力されるデジタル信号は、コンピュータ150内の記憶部152に記憶される。信号収集部140は、Data Acquisition System(DAS)とも呼ばれる。本明細書において電気信号は、アナログ信号もデジタル信号も含む概念である。なお、信号収集部140は、光照射部110の光射出部に取り付けられた光検出センサと接続されており、光が光照射部110から射出されたことをトリガーに、同期して処理を開始してもよい。また、信号収集部140は、フリーズボタンなどを用いてなされる指示をトリガーに同期して、当該処理を開始してもよい。
信号収集部140は、トランスデューサ121から出力されたアナログ信号である電気信号を増幅するアンプと、アンプから出力されたアナログ信号をデジタル信号に変換するA/D変換器とを含む。信号収集部140は、FPGA(Field Programmable Gate Array)チップなどで構成されてもよい。信号収集部140から出力されるデジタル信号は、コンピュータ150内の記憶部152に記憶される。信号収集部140は、Data Acquisition System(DAS)とも呼ばれる。本明細書において電気信号は、アナログ信号もデジタル信号も含む概念である。なお、信号収集部140は、光照射部110の光射出部に取り付けられた光検出センサと接続されており、光が光照射部110から射出されたことをトリガーに、同期して処理を開始してもよい。また、信号収集部140は、フリーズボタンなどを用いてなされる指示をトリガーに同期して、当該処理を開始してもよい。
(コンピュータ150)
情報処理装置としてのコンピュータ150は、演算部151、記憶部152、制御部153を含む。各構成の機能については処理フローの説明の際に説明する。
情報処理装置としてのコンピュータ150は、演算部151、記憶部152、制御部153を含む。各構成の機能については処理フローの説明の際に説明する。
演算部151としての演算機能を担うユニットは、CPUやGPU(Graphics Processing Unit)等のプロセッサ、FPGA(Field Programmable Gate Array)チップ等の演算回路で構成されることができる。これらのユニットは、単一のプロセッサや演算回路から構成されるだけでなく、複数のプロセッサや演算回路から構成されていてもよい。演算部151は、入力部170から、被検体音速や保持部の構成などの各種パラメータを受けて、受信信号を処理してもよい。
記憶部152は、ROM(Read only memory)、磁気ディスクやフラッシュメモリなどの非一時記憶媒体で構成することができる。また、記憶部152は、RAM(Random Access Memory)などの揮発性の媒体であってもよい。なお、プログラムが格納される記憶媒体は、非一時記憶媒体である。なお、記憶部152は、1つの記憶媒体から構成されるだけでなく、複数の記憶媒体から構成されていてもよい。
記憶部152は、後述する方法で演算部151により生成される光音響画像を示す画像データを保存することができる。
制御部153は、CPUなどの演算素子で構成される。制御部153は、光音響装置の各構成の動作を制御する。制御部153は、入力部170からの測定開始などの各種操作による指示信号を受けて、光音響装置の各構成を制御してもよい。また、制御部153は、記憶部152に格納されたプログラムコードを読み出し、光音響装置の各構成の作動を制御する。
コンピュータ150は専用に設計されたワークステーションであってもよい。また、コンピュータ150の各構成は異なるハードウェアによって構成されてもよい。また、コンピュータ150の少なくとも一部の構成は単一のハードウェアで構成されてもよい。
図3は、本実施形態に係るコンピュータ150の具体的な構成例を示す。本実施形態に係るコンピュータ150は、CPU154、GPU155、RAM156、ROM157、外部記憶装置158から構成される。また、コンピュータ150には、表示部160としての液晶ディスプレイ161、入力部170としてのマウス171、キーボード172が接続されている。
また、コンピュータ150および複数のトランスデューサ121は、共通の筺体に収められた構成で提供されてもよい。ただし、筺体に収められたコンピュータで一部の信号処理を行い、残りの信号処理を筺体の外部に設けられたコンピュータで行ってもよい。この場合、筺体の内部および外部に設けられたコンピュータを総称して、本実施形態に係るコンピュータとすることができる。すなわち、コンピュータを構成するハードウェアが一つの筺体に収められていなくてもよい。
(表示部160)
表示部160は、液晶ディスプレイや有機EL(Electro Luminescence)などのディスプレイである。コンピュータ150により得られた被検体情報等に基づく画像や特定位置の数値等を表示する装置である。表示部160は、画像や装置を操作するためのGUIを表示してもよい。なお、被検体情報の表示にあたっては、表示部160またはコンピュータ150において画像処理(輝度値の調整等)を行った上で表示することもできる。
表示部160は、液晶ディスプレイや有機EL(Electro Luminescence)などのディスプレイである。コンピュータ150により得られた被検体情報等に基づく画像や特定位置の数値等を表示する装置である。表示部160は、画像や装置を操作するためのGUIを表示してもよい。なお、被検体情報の表示にあたっては、表示部160またはコンピュータ150において画像処理(輝度値の調整等)を行った上で表示することもできる。
(入力部170)
入力部170としては、ユーザーが操作可能な、マウスやキーボードなどで構成される操作コンソールを採用することができる。また、表示部160をタッチパネルで構成し、表示部160を入力部170として利用してもよい。
入力部170としては、ユーザーが操作可能な、マウスやキーボードなどで構成される操作コンソールを採用することができる。また、表示部160をタッチパネルで構成し、表示部160を入力部170として利用してもよい。
入力部170は、観察したい位置や深さなどの補正量を決定するための情報や、補正量そのものなどを入力できるように構成されていてもよい。すなわち、入力部170は、これらの補正量に関する情報を入力できるように構成されていてもよい。入力方法としては、数値を入力してもよいし、スライダーバーを操作することにより入力ができてもよい。また、入力された情報に応じて表示部160に補正された画像が更新されていってもよい。これにより、ユーザーは自身の操作によって決定された補正量によって補正された画像を確認しながら、適切な補正量に設定できる。
なお、光音響装置の各構成はそれぞれ別の装置として構成されてもよいし、一体となった1つの装置として構成されてもよい。また、光音響装置の少なくとも一部の構成が一体となった1つの装置として構成されてもよい。
(被検体100)
被検体100は光音響装置を構成するものではないが、以下に説明する。本実施形態に係る光音響装置は、人や動物の悪性腫瘍や血管疾患などの診断や化学治療の経過観察などを目的として使用できる。よって、被検体100としては、生体、具体的には人体や動物の乳房や各臓器、血管網、頭部、頸部、腹部、手指および足指を含む四肢などの診断の対象部位が想定される。例えば、人体が測定対象であれば、オキシヘモグロビンあるいはデオキシヘモグロビンやそれらを含む多く含む血管あるいは腫瘍の近傍に形成される新生血管などを光吸収体の対象としてもよい。また、頸動脈壁のプラークなどを光吸収体の対象としてもよい。また、メチレンブルー(MB)、インドシニアングリーン(ICG)などの色素、金微粒子、またはそれらを集積あるいは化学的に修飾した外部から導入した物質を光吸収体としてもよい。
被検体100は光音響装置を構成するものではないが、以下に説明する。本実施形態に係る光音響装置は、人や動物の悪性腫瘍や血管疾患などの診断や化学治療の経過観察などを目的として使用できる。よって、被検体100としては、生体、具体的には人体や動物の乳房や各臓器、血管網、頭部、頸部、腹部、手指および足指を含む四肢などの診断の対象部位が想定される。例えば、人体が測定対象であれば、オキシヘモグロビンあるいはデオキシヘモグロビンやそれらを含む多く含む血管あるいは腫瘍の近傍に形成される新生血管などを光吸収体の対象としてもよい。また、頸動脈壁のプラークなどを光吸収体の対象としてもよい。また、メチレンブルー(MB)、インドシニアングリーン(ICG)などの色素、金微粒子、またはそれらを集積あるいは化学的に修飾した外部から導入した物質を光吸収体としてもよい。
(被検体情報を取得するためのフロー)
次に、本実施形態に係る情報処理を含む被検体情報の取得方法の各工程を、図4を参照して説明する。なお、各工程は、コンピュータ150が光音響装置の構成の動作を制御することにより実行される。
次に、本実施形態に係る情報処理を含む被検体情報の取得方法の各工程を、図4を参照して説明する。なお、各工程は、コンピュータ150が光音響装置の構成の動作を制御することにより実行される。
(S310:制御パラメータを指定する工程)
ユーザーが、被検体情報の取得のために必要な光照射部110の照射条件(繰り返し周波数や波長など)やプローブ180の位置などの制御パラメータを、入力部170を用いて指定する。
ユーザーが、被検体情報の取得のために必要な光照射部110の照射条件(繰り返し周波数や波長など)やプローブ180の位置などの制御パラメータを、入力部170を用いて指定する。
(S320:プローブを指定位置に移動させる工程)
制御部153が、ステップS301で指定された制御パラメータに基づいて、駆動部130にプローブ180を指定の位置へ移動させる。ステップS301において複数位置での撮像が指定された場合には、駆動部130は、まずプローブ180を最初の指定位置へ移動させる。なお、駆動部130は、測定の開始指示がなされたときに、あらかじめプログラムされた位置にプローブ180を移動させてもよい。
制御部153が、ステップS301で指定された制御パラメータに基づいて、駆動部130にプローブ180を指定の位置へ移動させる。ステップS301において複数位置での撮像が指定された場合には、駆動部130は、まずプローブ180を最初の指定位置へ移動させる。なお、駆動部130は、測定の開始指示がなされたときに、あらかじめプログラムされた位置にプローブ180を移動させてもよい。
(S330:光を照射する工程)
光照射部110は、S310で指定された制御パラメータに基づいて、被検体100に光を照射する。
光源111から発生した光は、光学系112を介してパルス光として被検体100に照射される。そして、被検体100内部でパルス光が吸収され、光音響効果により光音響波が生じる。光照射部110はパルス光の伝送と併せて信号収集部140へ同期信号を送信する。
光照射部110は、S310で指定された制御パラメータに基づいて、被検体100に光を照射する。
光源111から発生した光は、光学系112を介してパルス光として被検体100に照射される。そして、被検体100内部でパルス光が吸収され、光音響効果により光音響波が生じる。光照射部110はパルス光の伝送と併せて信号収集部140へ同期信号を送信する。
(S340:光音響波を受信する工程)
信号収集部140は、光照射部110から送信された同期信号を受信すると、信号収集の動作を開始する。すなわち、信号収集部140は、受信部120から出力された、音響波に由来するアナログ電気信号を、増幅・AD変換することにより、増幅されたデジタル電気信号を生成し、コンピュータ150へ出力する。コンピュータ150は、信号収集部140から送信された信号を記憶部152に保存する。ステップS301で複数の走査位置での撮像を指定した場合には、指定した走査位置において、S320−S340のステップを繰り返し実行し、パルス光の照射と音響波に由来するデジタル信号の生成を繰り返す。
信号収集部140は、光照射部110から送信された同期信号を受信すると、信号収集の動作を開始する。すなわち、信号収集部140は、受信部120から出力された、音響波に由来するアナログ電気信号を、増幅・AD変換することにより、増幅されたデジタル電気信号を生成し、コンピュータ150へ出力する。コンピュータ150は、信号収集部140から送信された信号を記憶部152に保存する。ステップS301で複数の走査位置での撮像を指定した場合には、指定した走査位置において、S320−S340のステップを繰り返し実行し、パルス光の照射と音響波に由来するデジタル信号の生成を繰り返す。
(S350:音響波の伝搬経路上の媒質の音速情報を取得する工程)
音速取得手段としてのコンピュータ150は、音響波の伝搬経路上の媒質の音速情報を取得する。本実施形態において、コンピュータ150は、被検体100の音速情報及び音響マッチング材210の音速情報を取得する。コンピュータ150は、被検体100の音速値ch及び音響マッチング材210の音速値cmを音速情報として取得する。
音速取得手段としてのコンピュータ150は、音響波の伝搬経路上の媒質の音速情報を取得する。本実施形態において、コンピュータ150は、被検体100の音速情報及び音響マッチング材210の音速情報を取得する。コンピュータ150は、被検体100の音速値ch及び音響マッチング材210の音速値cmを音速情報として取得する。
音響マッチング材210の音速情報について、予め測定された音速情報を記憶部152に保存しておき、本工程において演算部151が記憶部152から読み出すことにより取得してもよい。また、音響マッチング材210の温度に対する音速情報の関係式または関係テーブルが記憶部152に予め保存されていてもよい。そして、本工程において温度センサ190が音響マッチング材210の温度を測定し、コンピュータ150が測定された温度に対応する音速情報を関係式または関係テーブルに従って取得してもよい。また、ユーザーが入力部170により音響マッチング材210の音速情報を入力し、コンピュータ150がその情報を受け取ることにより取得してもよい。
被検体100での音速情報についても、上述した音響マッチング材210の音速情報と同様の方法で取得してもよい。ただし、被検体100の音速情報は被検体100毎に異なるため、被検体100毎に新たなデータを取得することが好ましい。コンピュータ150は、被検体100から発生した音響波に由来する信号を用いて被検体100の音速情報を取得してもよい。この方法によれば、装置規模を大きくすることなく、被検体100の固有の音速情報を取得することができる。また、被検体100から発生した音響波に由来する信号を用いて音響マッチング材210の音速情報を取得してもよい。
また、音速値が文献等で知られている媒質については、演算部151は、記憶部152に保存された既知の音速情報を読み出すことにより、取得してもよい。また、音速情報の取得方法としては、超音波の透過などを利用して実験的に測定してもよいし、超音波CTなど他モダリティを用いて測定してもよい。
なお、本工程では音速値を求めることのできるパラメータであればいかなるパラメータを音速情報として取得してもよい。例えば、密度ρと体積弾性率Kとから音速値は求められるため、本工程において密度ρ及び体積弾性率Kを音速情報として取得し、それらのパラメータから音速値を求めてもよい。また、被検体100や音響マッチング材210等を伝搬する音響波の伝搬経路における音速値の空間分布を示す情報を、音速情報として取得してもよい。
なお、コンピュータ150は、その他公知の手法により被検体100または音響マッチング材210の音速情報を取得してもよい。
(S360:音速情報を補正する工程)
補正手段としてのコンピュータ150は、ステップS350で取得された音速情報に対する補正量を取得する。そして、補正手段としてのコンピュータ150は、取得された補正量を用いて、S350で取得された音速情報を補正する。本実施形態では、コンピュータ150が、音響マッチング材210の音速値cmに対する補正量Δcmを取得し、音響マッチング材210の補正された音速値cm+Δcmを補正された音速情報として取得する。すなわち、コンピュータ150は、S350で取得された音速値に、補正量を加減算することにより、補正された音速値を取得している。
補正手段としてのコンピュータ150は、ステップS350で取得された音速情報に対する補正量を取得する。そして、補正手段としてのコンピュータ150は、取得された補正量を用いて、S350で取得された音速情報を補正する。本実施形態では、コンピュータ150が、音響マッチング材210の音速値cmに対する補正量Δcmを取得し、音響マッチング材210の補正された音速値cm+Δcmを補正された音速情報として取得する。すなわち、コンピュータ150は、S350で取得された音速値に、補正量を加減算することにより、補正された音速値を取得している。
補正量の取得方法としては、演算部151が記憶部152に保存された所定の補正量を読み出すことにより、取得してもよい。また、ユーザーが入力部170を用いて補正量を決定するための指示を行い、演算部151が入力部170から送信された指示情報に基づいて、補正量を取得してもよい。また、演算部151は、記憶部152に保存された、補正量を求めるための情報と補正量との関係を示す関係テーブルまたは関係式を用いて、補正量を取得してもよい。
以下、記憶部152に保存する補正量、関係テーブル、または関係式を、図5に示す計算モデルを用いたシミュレーションにより決定する方法を説明する。まず、計算モデルの詳細を説明する。受信部120は、512個の複数のトランスデューサ121が半球状の支持体122の内周面に配置された部材である。複数のトランスデューサ121はφ1.5mmの振動子とした。また、複数のトランスデューサ121は、図2(b)に示すように半球上で3次元スパイラルを形成するように配置した。また、被検体100の音速値を1600m/sに設定した。また、被検体100の内部、かつ、支持体122の曲率中心点に、対称性を有するφ1mmの球形をした音源を設置した。被検体100と受信部120との間は音響マッチング材210である水で満たされている。音響マッチング材210の音速値は1500m/sに設定した。なお、φ1mmの球形をした音源から、被検体100と音響マッチング材210との境界までの距離は15mmである。
注目位置101(典型的にはピクセルまたはボクセル)に存在する音源から発生した音響波を、複数のトランスデューサ121で受信し、40MHzのサンプリング周波数でサンプリングして得られる受信信号を生成する。図6は光照射を行ったタイミングにおける受信部120の位置(半球状の支持体の曲率中心の位置)をプロットしたものである。図6に示すように、受信部120は、XY平面内でスパイラル状に移動し、受信部120と被検体100との相対位置を変化させながら、1024点で光照射することにより発生した音響波を複数のトランスデューサ121で受信する。
受信信号を生成するために、演算部151は、被検体100と音響マッチング材210の境界面上の点Xを取りうる座標の範囲で変更させて、注目位置101と点Xの間との伝搬時間と、点Xとトランスデューサ121との間の伝搬時間とを算出する。演算部151は、注目位置101と点Xの間との伝搬時間と、点Xとトランスデューサ121との間の伝搬時間の合計が最小となる点Xを算出する。これを計算対象としているすべてのトランスデューサ121及びすべての注目位置101について求める。なお、注目位置101―点X―トランスデューサ121の伝搬時間が最小となるときの伝搬経路が、屈折した場合の伝搬経路103であると仮定している。そして、演算部151は、この伝搬経路103に相当する伝搬時間にトランスデューサ121が音響波を受信したものとして受信信号を生成する。
そして、演算部151は、生成された受信信号を用いて、Universal back−projection(UBP)法によって再構成し、再構成画像を生成する。UBP法の詳細に関して後述する。再構成の際のボクセル数(X×Y×Z)は、512×512×512とし、ボクセルピッチは0.1mmとした。
ただし、注目位置101と複数のトランスデューサ121を結ぶ直線102(図5の点線)と、実際の屈折した伝搬経路103とは異なる。そのため、音響波の伝搬経路を直線近似して再構成した場合(音響波の伝搬経路を直線102として伝搬時間を計算して再構成する場合)には、屈折による画質の低下が生じることが理解される。
そこで、演算部151は、音響マッチング材210の音速値を複数パターン変更して再構成を行い、再構成画像の画質が高くなる音速値を決定する。さらに、演算部151は、このように再構成画像を評価して決定された、画質が高くなる音速値と、計算モデルの音速値の設定値と差分を、音響マッチング材210の音速値に対する補正量として決定する。演算部151は、被検体100の音速値、音響マッチング材210の音速値、及び注目位置101の位置の複数の組み合わせに対してこの処理を行い、各場合の補正量を算出する。そして、被検体100の音速値、音響マッチング材210の音速値、及び注目位置101の位置の組み合わせと、音響マッチング材210に対する補正量との関係を示す関係式または関係テーブルを記憶部152に保存する。あるいは、複数の組み合わせに対して共通で利用できる補正量の値自体を、記憶部152に保存してもよい。そして、後述するように、複数の組み合わせに対して共通の補正量を記憶部152から読み出し、各組み合わせに共通の補正量を適用してもよい。
なお、駆動部による、被検体100と受信部120の相対位置の移動がXY平面のみである場合は、補正量はXY方向よりも、Z方向(深さ方向)の方が敏感度が高い傾向がある。そのため、注目位置を含む深さの平面については共通の補正量を適用してもよい。また、注目位置の含む深さの近傍の深さについても共通の補正量を適用してもよい。すなわち、特定の深さ、または、特定の深さ及びその近傍の深さに存在する複数の注目位置に対して共通の補正量を適用してもよい。記憶部152に保存する補正量、関係テーブル、または関係式のデータ容量を抑制することができる。
なお、再構成画像の画質の評価方法については、ユーザーが再構成画像に関する情報(画像そのものや画像強度のラインプロファイルなど)を確認し、入力部170を用いて補正量の基準となる画質の高い再構成画像を指定してもよい。あるいは、演算部151が、再構成画像に基づいて画質を評価してもよい。例えば、演算部151は、再構成画像の画像強度のラインプロファイルを生成し、画像強度のピークの両端に現れる、ラインプロファイルの2つの谷の極値の画像強度差を評価し、この画像強度差が小さいほど画質が高いと評価してもよい。この手法では、UBP法による再構成により発生する特有の画像特徴を利用してラインプロファイルの対称性を評価し、画質を評価している。この他にも、再構成画像の画質を評価できる方法であればあらゆる方法を採用することができる。
前述したように、ユーザーが入力部170を用いて補正量を決定するための指示を行い、演算部151が入力部170から送信された指示情報に基づいて、補正量を取得してもよい。ユーザーは、補正量の値自体を指示することや、補正量を決定するために必要な情報を、入力部170を用いて指示することができる。例えば、補正量を決定するために必要な情報としては、音響波の伝搬経路上の媒質の音速情報(例えば、被検体100や音響マッチング材210の音速情報)、補正したい注目位置などである。演算部151は、これらの指示情報と、記憶部152に保存されている、これらの情報と補正量との関係を示す関係テーブルまたは関係式とに基づいて、補正量を決定することができる。以下、ユーザーが入力部170を用いて補正量を決定するための指示を行う場合に表示部160に表示されるGUIの例を、図7を参照して説明する。
例えば、図7(a)に示すように、GUI上で補正量の値自体をテキストボックス701に入力部170を用いて入力してもよい。図7(b)に示すように、ユーザーが画像中の深さ(Z方向の位置)を指定可能なスライドバー702を操作して、補正したい注目位置を含む深さを指示してもよい。また、図7(c)に示すように、ユーザーがスライドバー703を操作して、補正したい注目位置を含む深さを複数にしてもよい。さらに、図7(d)に示すように、ユーザーがスライドバー704を操作することにより、図7(b)や図7(c)で指定された深さに対する補正量を指示できるようにしてもよい。
ここで、深さ方向とは、受信部120の音響波に対する指向方向のことを指す。本実施形態において深さ方向とは、半球状の支持体122の開口の面外方向(図2(a)の紙面上下方向)のことを指す。また、複数のトランスデューサが2次元に配置される場合、トランスデューサの配置面の面外方向が深さ方向となる。
補正量は、音響波の伝搬経路上のいかなる媒質に対する補正量であってもよい。1つの媒質に対して補正量を取得してもよいし、複数の媒質に対して補正量を取得してもよい。なお、屈折に伴う伝搬時間に与える影響の大きい、伝搬距離の長い媒質に対する補正量を取得することが好ましい。本実施形態の場合、被検体100よりも音響マッチング材210を伝搬する距離が典型的には長いため、音響マッチング材210の音速情報のみを補正するための補正量を取得してもよい。また、音響マッチング材210を構成する材料の温度と音速値との関係が既知である場合、音響マッチング材210の温度情報に基づいて決定された音速値を、本工程における補正量で補正してもよい。一方、被検体100については、構成される物質に個体差があるため、温度情報を用いる方法とは別の方法で音速情報が取得され、補正されなくてもよい。
(S370:補正された音速情報に基づいて、音響波の伝搬経路を直線近似して再構成する工程)
情報取得手段としてのコンピュータ150は、S340で取得された受信信号及びステップS360にて補正された音速情報に基づいて、音響波の伝搬経路を直線近似して再構成することにより注目位置の被検体情報を取得する。コンピュータ150は、複数の注目位置の被検体情報を取得することにより、被検体情報の空間分布(画像データ)を取得することができる。
情報取得手段としてのコンピュータ150は、S340で取得された受信信号及びステップS360にて補正された音速情報に基づいて、音響波の伝搬経路を直線近似して再構成することにより注目位置の被検体情報を取得する。コンピュータ150は、複数の注目位置の被検体情報を取得することにより、被検体情報の空間分布(画像データ)を取得することができる。
例えば、演算部151は、被検体情報として音響波の発生音圧(初期音圧)を取得する再構成の手法として、式(1)で表されるUniversal back−projection(UBP)法を採用してもよい。
ここで、r0は再構成する位置(再構成位置、注目位置とも呼ぶ)を示す位置ベクトル、p0(r0,t)は再構成する位置での初期音圧、cは伝搬経路の音速を示す。また、ΔΩiは再構成する位置からi番目のトランスデューサ121を見込む立体角、Nは再構成に用いるトランスデューサ121の個数を示す。式(1)は、受信信号p(ri,t)に微分等の処理を行い、それらに立体角の加重をかけて整相加算すること(逆投影)を示している。式(1)のtは、注目位置とトランスデューサ121とを結ぶ音線を光音響波が伝搬する時間(伝搬時間)である。なお、b(ri、t)の計算においては、他にも演算処理を施してもよい。例えば、周波数フィルタリング(ローパス、ハイパス、バンドパス等)、デコンボリューション、包絡線検波、ウェーブレットフィルタリング、等である。また、本発明においては、トランスデューサと注目位置とを結ぶ音線の伝搬時間を求めて再構成する方法であれば、どのような再構成アルゴリズムを用いても良い。例えば、タイムドメインでの逆投影法としてFiltered back−projection(FBP)や整相加算(Delay−and−Sum)などを採用してもよい。
図5において、符号102は、注目位置101からトランスデューサ121までの音響波の伝搬経路を直線近似したときの音線(伝搬経路)を表す。直線近似したときの被検体100内の伝搬距離をD1とし、音響マッチング材210内の伝搬距離をD2とする。この場合、演算部151は、式(1)の伝搬時間tは、式(2)で表すことができる。
ここで、D1は直線近似したときの被検体100内の伝搬距離、D2は直線近似したときの音響マッチング材210内の伝搬距離を示し、音響波を受信したときのトランスデューサ121の位置情報に基づいて決定される。なお、トランスデューサ121の位置情報は光照射部110からの同期信号を受信したときに駆動部130から出力された受信部120の位置情報に基づいて取得されてもよい。また、受信部120に備えられた磁気センサなどの位置取得手段からの出力に基づいて、トランスデューサ121の位置情報を取得してもよい。
また、chは被検体100の音速値を示し、本実施形態ではステップS350で取得される値を用いる。cmは音響マッチング材210の音速値を示し、Δcmは音響マッチング材210の音速値に対する補正量を示す。また、cm+Δcmは補正された音響マッチング材210の音速値を示し、本実施形態ではステップS360で取得される値を用いる。
本実施形態では、演算部151が、式(1)及び式(2)で示したように、音響波の伝搬経路上の媒質の補正された音速情報を用いて、直線近似して再構成することにより、注目位置の被検体情報を取得する。この方法によれば、屈折した音響波の伝搬経路を算出するために必要な計算を省略することができる。さらに、音速情報を補正し、補正された音速情報に基づいて直線近似した再構成を行うことにより、直線近似による被検体情報の取得精度の低下を抑制することができる。
(S380:再構成画像を表示する工程)
表示制御手段としてのコンピュータ150は、S370で取得された画像データを表示部160に送信し、画像データを再構成画像として表示部160に表示させる。なお、コンピュータ150は、被検体情報の画像データを画像として表示させる他に、画像データの特定位置の被検体情報の数値を表示させてもよい。また、補正前の音速情報に基づいて得られる再構成画像と、補正後の音速情報に基づいて得られる再構成画像とを並列表示したり、切り換えて表示したりしてもよい。
表示制御手段としてのコンピュータ150は、S370で取得された画像データを表示部160に送信し、画像データを再構成画像として表示部160に表示させる。なお、コンピュータ150は、被検体情報の画像データを画像として表示させる他に、画像データの特定位置の被検体情報の数値を表示させてもよい。また、補正前の音速情報に基づいて得られる再構成画像と、補正後の音速情報に基づいて得られる再構成画像とを並列表示したり、切り換えて表示したりしてもよい。
ここで、本実施形態に係る情報処理方法の効果を説明するために行ったシミュレーションの計算結果を説明する。本シミュレーションでは、図5に示す計算モデルを用いた。
まず、比較例として、被検体100の音速値を1600m/s、音響マッチング材210の音速値を1500m/sというシミュレーションの設定値を用いて、直線近似で再構成した場合に得られる再構成画像を図8(a)に示す。また、図8(a)は、Y=0におけるZ方向の画像強度のラインプロファイルも示す。図8(a)に示すラインプロファイルを見ると、対象性を有する音源が非対称に画像化されており、+Z方向にアーチファクトが発生していることからも、屈折の影響による画質低下が生じていることが理解される。
続いて、本実施形態に係る手法で得られた再構成画像を図8(b)に示す。本方法では、被検体100の音速値1600m/s、音響マッチング材210の音速値1500m/s、注目位置101に相当する深さ15mmの組み合わせに対応する、音響マッチング材210の音速値に対する補正量+4.5m/sを用いる。この補正量を用いて、被検体100の音速値を1600m/s、音響マッチング材210の補正後の音速値を1504.5m/sとし、伝搬経路を直線近似して再構成した。その結果、図8(b)に示す再構成画像が得られた。図8(b)に示す再構成画像は、図8(a)と比較して、Y=0におけるZ方向の画像強度のラインプロファイルが対称的になり、+Z方向のアーチファクトは減少している。すなわち、本実施形態に係る手法では、屈折による再構成画像の画質低下が抑制されていることが理解される。
なお、本シミュレーションでは、複数のトランスデューサ121が512個、受信部120の走査位置が1024点、再構成領域のボクセル数が512×512×512である。そのため、各注目位置についてスネルの法則に基づいて屈折による音響波の伝搬経路を算出しようとすると、512×1024×512×512通りの計算をしなければならない。そのため、この方法では、その計算量は膨大であり、再構成時間の長時間化に繋がる。
一方、本実施形態に係る手法の場合、補正された音速情報に基づいて直線近似して再構成を行うことにより、計算量が膨大である屈折した音響波の伝搬経路を算出することなく、屈折による画質低下を抑制することができる。
(その他の実施例)
また、本発明は、以下の処理を実行することによっても実現される。即ち、上述した実施形態の機能を実現するソフトウェア(プログラム)を、ネットワーク又は各種記憶媒体を介してシステム或いは装置に供給し、そのシステム或いは装置のコンピュータ(またはCPUやMPU等)がプログラムを読み出して実行する処理である。
また、本発明は、以下の処理を実行することによっても実現される。即ち、上述した実施形態の機能を実現するソフトウェア(プログラム)を、ネットワーク又は各種記憶媒体を介してシステム或いは装置に供給し、そのシステム或いは装置のコンピュータ(またはCPUやMPU等)がプログラムを読み出して実行する処理である。
120 受信部
150 コンピュータ
150 コンピュータ
Claims (13)
- 被検体内の注目位置から伝搬する音響波が複数の受信手段で受信されることにより得られる複数の信号を用いて、前記注目位置の被検体情報を取得する情報処理装置であって、
音響波の伝搬経路上の媒質の音速情報を取得する音速取得手段と、
前記音速情報に対する補正量を取得し、当該補正量に基づいて、前記音速情報を補正する補正手段と、
補正された前記音速情報に基づいて、前記注目位置から前記複数の受信手段までの伝搬経路を直線近似して音響波の伝搬時間を決定し、前記複数の信号及び前記伝搬時間に基づいて、前記注目位置の前記被検体情報を取得する情報取得手段と、
を有する
ことを特徴とする情報処理装置。 - 前記補正量が保存された記憶手段を有し、
前記補正手段は、前記記憶手段から前記補正量を読み出すことにより取得する
ことを特徴とする請求項1に記載の情報処理装置。 - 前記音速情報、前記複数の受信手段の位置情報、前記注目位置の位置情報、及び前記補正量の関係を示す関係テーブルまたは関係式が保存された記憶手段を有し、
前記補正手段は、
前記複数の受信手段の位置情報及び前記注目位置の位置情報を取得し、
前記関係テーブルまたは前記関係式を前記記憶手段から読み出し、
前記音速情報、前記複数の受信手段の前記位置情報、前記注目位置の前記位置情報、及び前記関係テーブルまたは前記関係式に基づいて、前記補正量を取得する
ことを特徴とする請求項1に記載の情報処理装置。 - 前記補正手段は、
ユーザーの指示に基づいて決定される前記注目位置の前記位置情報を取得する
ことを特徴とする請求項3に記載の情報処理装置。 - 前記補正手段は、
ユーザーの指示に基づいて前記補正量を取得する
ことを特徴とする請求項1から4のいずれか1項に記載の情報処理装置。 - 前記音速取得手段は、温度センサにより測定された前記媒質の温度情報に基づいて、前記音速情報を取得し、
前記補正手段は、前記温度情報に基づいて取得された前記音速情報を補正する
ことを特徴とする請求項1から5のいずれか1項に記載の情報処理装置。 - 前記媒質は、前記被検体、及び、前記被検体と前記複数の受信手段との間に配置された音響マッチング材を含み、
前記音速取得手段は、
前記温度センサにより測定された前記音響マッチング材の温度情報に基づいて、前記音響マッチング材の音速情報を取得し、
前記補正手段は、前記温度情報に基づいて取得された前記音響マッチング材の音速情報を補正し、
前記情報取得手段は、前記音響マッチング材の補正された前記音速情報と、前記音速取得手段により取得された前記被検体の前記音速情報とに基づいて、前記注目位置から前記複数の受信手段までの伝搬経路を直線近似して音響波の伝搬時間を決定する
ことを特徴とする請求項6に記載の情報処理装置。 - 前記補正手段は、
特定の深さに存在する複数の前記注目位置に対して、前記補正量に基づいて、前記音速情報を補正し、
前記情報取得手段は、
複数の前記注目位置に対応する補正された前記音速情報に基づいて、複数の前記注目位置から前記複数の受信手段までの伝搬経路を直線近似して、複数の前記注目位置に対応する前記伝搬時間を決定し、
前記複数の信号及び複数の前記注目位置に対応する前記伝搬時間に基づいて、複数の前記注目位置の前記被検体情報を取得する
ことを特徴とする請求項1から7のいずれか1項に記載の情報処理装置。 - 前記補正手段は、
特定の深さまたは前記特定の深さの近傍の深さに存在する複数の前記注目位置に対して、前記補正量に基づいて、前記音速情報を補正し、
前記情報取得手段は、
複数の前記注目位置に対応する補正された前記音速情報に基づいて、複数の前記注目位置から前記複数の受信手段までの伝搬経路を直線近似して、複数の前記注目位置に対応する前記伝搬時間を決定し、
前記複数の信号及び複数の前記注目位置に対応する前記伝搬時間に基づいて、複数の前記注目位置の前記被検体情報を取得する
ことを特徴とする請求項1から7のいずれか1項に記載の情報処理装置。 - 前記被検体に光照射を行う光照射手段を有し、
前記複数の受信手段は、前記光照射に起因して発生する音響波を受信することにより、前記複数の信号を出力する
ことを特徴とする請求項1から9のいずれか1項に記載の情報処理装置。 - 前記被検体に超音波を送信する送信手段を有し、
前記複数の受信手段は、前記超音波が前記被検体の内部で反射することにより発生する音響波を受信することにより、前記複数の信号を出力する
ことを特徴とする請求項1から9のいずれか1項に記載の情報処理装置。 - 被検体内の注目位置から伝搬する音響波が複数の受信手段で受信されることにより得られる複数の信号を用いて、前記注目位置の被検体情報を取得する情報処理方法であって、
音響波の伝搬経路上の媒質の音速情報を取得し、
前記音速情報に対する補正量を取得し、
前記補正量を用いて、前記音速情報を補正し、
補正された前記音速情報に基づいて、前記注目位置から前記複数の受信手段までの伝搬経路を直線近似して音響波の伝搬時間を決定し、
前記複数の信号及び前記伝搬時間に基づいて、前記注目位置の前記被検体情報を取得する
ことを特徴とする情報処理方法。 - 請求項12に記載の情報処理方法をコンピュータに実行させるためのプログラム。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2016202044A JP2018061716A (ja) | 2016-10-13 | 2016-10-13 | 情報処理装置、情報処理方法、及びプログラム |
US15/727,904 US10436706B2 (en) | 2016-10-13 | 2017-10-09 | Information processing apparatus, information processing method, and storage medium |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2016202044A JP2018061716A (ja) | 2016-10-13 | 2016-10-13 | 情報処理装置、情報処理方法、及びプログラム |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2018061716A true JP2018061716A (ja) | 2018-04-19 |
Family
ID=61903826
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2016202044A Withdrawn JP2018061716A (ja) | 2016-10-13 | 2016-10-13 | 情報処理装置、情報処理方法、及びプログラム |
Country Status (2)
Country | Link |
---|---|
US (1) | US10436706B2 (ja) |
JP (1) | JP2018061716A (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2020162745A (ja) * | 2019-03-28 | 2020-10-08 | キヤノン株式会社 | 画像処理装置、画像処理方法及びプログラム |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114264721B (zh) * | 2021-11-26 | 2024-04-19 | 海鹰企业集团有限责任公司 | 新型声速处理系统 |
CN117158911B (zh) * | 2023-10-25 | 2024-01-23 | 杭州励影光电成像有限责任公司 | 一种多声速自适应光声层析图像重建方法 |
Family Cites Families (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6607489B2 (en) * | 2001-04-05 | 2003-08-19 | General Electric Company | Focus correction for ultrasound imaging through mammography compression plate |
JP5692986B2 (ja) * | 2008-12-25 | 2015-04-01 | キヤノン株式会社 | 被検体情報取得装置 |
EP2818117A1 (en) * | 2013-06-26 | 2014-12-31 | Canon Kabushiki Kaisha | Object information obtaining system, signal processing method, and program |
US20160213257A1 (en) * | 2013-09-04 | 2016-07-28 | Canon Kabushiki Kaisha | Photoacoustic apparatus |
JP6223129B2 (ja) * | 2013-10-31 | 2017-11-01 | キヤノン株式会社 | 被検体情報取得装置、表示方法、被検体情報取得方法、及びプログラム |
US9730589B2 (en) * | 2013-10-31 | 2017-08-15 | Canon Kabushiki Kaisha | Examined-portion information acquisition apparatus |
JP6544910B2 (ja) * | 2014-11-07 | 2019-07-17 | キヤノン株式会社 | 情報処理装置、被検体情報取得装置及び音速決定方法 |
JP6478572B2 (ja) * | 2014-11-10 | 2019-03-06 | キヤノン株式会社 | 被検体情報取得装置および音響波装置の制御方法 |
JP6512969B2 (ja) * | 2015-07-06 | 2019-05-15 | キヤノン株式会社 | 処理装置、光音響装置、処理方法、及びプログラム |
US20170086680A1 (en) * | 2015-09-29 | 2017-03-30 | Canon Kabushiki Kaisha | Apparatus for acquiring object information, information processing method, and non-transitory computer-readable storage medium storing program |
JP2017070488A (ja) * | 2015-10-07 | 2017-04-13 | キヤノン株式会社 | 被検体情報取得装置 |
JP6598667B2 (ja) * | 2015-12-17 | 2019-10-30 | キヤノン株式会社 | 被検体情報取得装置およびその制御方法 |
CN106913314B (zh) * | 2015-12-25 | 2020-01-21 | 佳能株式会社 | 信息获取装置和信息获取方法 |
JP2017184972A (ja) * | 2016-04-04 | 2017-10-12 | キヤノン株式会社 | 処理システム、情報取得装置、信号処理方法、プログラム |
JP2018061725A (ja) * | 2016-10-14 | 2018-04-19 | キヤノン株式会社 | 被検体情報取得装置および信号処理方法 |
JP6777511B2 (ja) * | 2016-11-22 | 2020-10-28 | 株式会社日立製作所 | 超音波撮像装置 |
JP6594355B2 (ja) * | 2017-01-06 | 2019-10-23 | キヤノン株式会社 | 被検体情報処理装置および画像の表示方法 |
JP2018126389A (ja) * | 2017-02-09 | 2018-08-16 | キヤノン株式会社 | 情報処理装置、情報処理方法、およびプログラム |
-
2016
- 2016-10-13 JP JP2016202044A patent/JP2018061716A/ja not_active Withdrawn
-
2017
- 2017-10-09 US US15/727,904 patent/US10436706B2/en not_active Expired - Fee Related
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2020162745A (ja) * | 2019-03-28 | 2020-10-08 | キヤノン株式会社 | 画像処理装置、画像処理方法及びプログラム |
JP7277212B2 (ja) | 2019-03-28 | 2023-05-18 | キヤノン株式会社 | 画像処理装置、画像処理方法及びプログラム |
Also Published As
Publication number | Publication date |
---|---|
US20180106716A1 (en) | 2018-04-19 |
US10436706B2 (en) | 2019-10-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10582910B2 (en) | Information acquisition apparatus and information acquisition method | |
US20200205749A1 (en) | Display control apparatus, display control method, and non-transitory computer-readable medium | |
JP6472437B2 (ja) | 光音響装置、及び音響波受信装置 | |
JP6742745B2 (ja) | 情報取得装置および表示方法 | |
US10436706B2 (en) | Information processing apparatus, information processing method, and storage medium | |
KR20170074171A (ko) | 광음향장치, 정보 취득장치, 정보 취득방법, 및 프로그램 | |
JP2017529913A (ja) | 光音響装置 | |
EP3329843B1 (en) | Display control apparatus, display control method, and program | |
US11510630B2 (en) | Display control apparatus, image display method, and non-transitory computer-readable medium | |
KR20200067254A (ko) | 정보 처리장치, 정보 처리방법, 및 비일시적인 기억매체 | |
JP7013215B2 (ja) | 情報処理装置、情報処理方法、プログラム | |
WO2018230409A1 (ja) | 情報処理装置、情報処理方法、及びプログラム | |
JP7277212B2 (ja) | 画像処理装置、画像処理方法及びプログラム | |
JP6929048B2 (ja) | 表示制御装置、表示方法、及びプログラム | |
US20200305727A1 (en) | Image processing device, image processing method, and program | |
JP6929204B2 (ja) | 情報処理装置、情報処理方法、およびプログラム | |
JP6882108B2 (ja) | 画像生成装置、画像生成方法、及びプログラム | |
JP6614968B2 (ja) | 情報処理装置、被検体情報取得装置、情報処理方法、及びプログラム | |
JP2019097804A (ja) | 情報処理装置、情報処理方法、プログラム |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20190924 |
|
A761 | Written withdrawal of application |
Free format text: JAPANESE INTERMEDIATE CODE: A761 Effective date: 20200703 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20200710 |