JP2016144754A - Apparatus and image generation method - Google Patents
Apparatus and image generation method Download PDFInfo
- Publication number
- JP2016144754A JP2016144754A JP2016100478A JP2016100478A JP2016144754A JP 2016144754 A JP2016144754 A JP 2016144754A JP 2016100478 A JP2016100478 A JP 2016100478A JP 2016100478 A JP2016100478 A JP 2016100478A JP 2016144754 A JP2016144754 A JP 2016144754A
- Authority
- JP
- Japan
- Prior art keywords
- speed
- image
- subject
- detection
- image generation
- 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
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000001514 detection method Methods 0.000 claims abstract description 150
- 239000000523 sample Substances 0.000 claims description 18
- 230000001902 propagating effect Effects 0.000 claims description 10
- 230000007423 decrease Effects 0.000 claims description 5
- 238000012937 correction Methods 0.000 claims description 2
- 238000012545 processing Methods 0.000 abstract description 24
- 230000006866 deterioration Effects 0.000 abstract description 6
- 230000005540 biological transmission Effects 0.000 abstract description 4
- 238000003325 tomography Methods 0.000 description 12
- 230000003287 optical effect Effects 0.000 description 9
- 239000006096 absorbing agent Substances 0.000 description 8
- 238000010586 diagram Methods 0.000 description 6
- 230000035945 sensitivity Effects 0.000 description 6
- 238000003384 imaging method Methods 0.000 description 5
- 210000004204 blood vessel Anatomy 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 241001465754 Metazoa Species 0.000 description 3
- 238000003745 diagnosis Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 102000001554 Hemoglobins Human genes 0.000 description 2
- 108010054147 Hemoglobins Proteins 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 2
- 201000011510 cancer Diseases 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- HFGPZNIAWCZYJU-UHFFFAOYSA-N lead zirconate titanate Chemical compound [O-2].[O-2].[O-2].[O-2].[O-2].[Ti+4].[Zr+4].[Pb+2] HFGPZNIAWCZYJU-UHFFFAOYSA-N 0.000 description 2
- 229910052451 lead zirconate titanate Inorganic materials 0.000 description 2
- 239000004973 liquid crystal related substance Substances 0.000 description 2
- 229920001817 Agar Polymers 0.000 description 1
- 206010006187 Breast cancer Diseases 0.000 description 1
- 208000026310 Breast neoplasm Diseases 0.000 description 1
- 239000008272 agar Substances 0.000 description 1
- 210000000481 breast Anatomy 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000012631 diagnostic technique Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 229940028435 intralipid Drugs 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 230000001172 regenerating effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
- 208000019553 vascular disease Diseases 0.000 description 1
- 238000007794 visualization technique Methods 0.000 description 1
Images
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
Description
本発明は、被検体内部の情報を表す画像を再構成する画像生成装置、画像生成方法、及び、プログラムに関する。 The present invention relates to an image generation apparatus, an image generation method, and a program for reconstructing an image representing information inside a subject.
光源(レーザーなど)からの光を生体などの被検体に照射し、該被検体の内部の情報を可視化する画像化装置の研究が医療分野で積極的に進められている。このような光による可視化技術の一つとして、Photoacoustic Tomography(PAT:光音響トモグラフィー)が
ある。光音響トモグラフィー装置では、被検体内を伝播し、拡散した光のエネルギーを吸収した生体組織から発生される音響波(典型的に超音波である)を、被検体を取り囲む複数の個所で検出する。そして、得られた信号を数学的に解析処理し、被検体内部の光学特性値、特に吸収係数分布に関連した情報を可視化する。近年、この光音響トモグラフィー装置を用いて、小動物の血管像をイメージングする前臨床研究や、この原理を乳がんなどの診断に応用する臨床研究が積極的に進められている。
Research on an imaging apparatus that irradiates a subject such as a living body with light from a light source (such as a laser) and visualizes information inside the subject is being actively promoted in the medical field. One of such light visualization techniques is photoacoustic tomography (PAT). In a photoacoustic tomography apparatus, acoustic waves (typically ultrasonic waves) generated from a living tissue that has propagated through a subject and absorbed diffused light energy are detected at a plurality of locations surrounding the subject. . Then, the obtained signal is mathematically analyzed to visualize the optical characteristic value inside the subject, particularly information related to the absorption coefficient distribution. In recent years, preclinical research for imaging blood vessels of small animals using this photoacoustic tomography apparatus and clinical research for applying this principle to diagnosis of breast cancer and the like have been actively promoted.
このような光音響トモグラフィー装置や、従来から医療分野で使われている超音波診断装置(生体内で反射した音響波を検出し、画像を再構成する装置)では、通常、被検体の平均的な音速(被検体内部における音響波の伝播速度)を用いて、画像が再構成される。一般的に、伝播速度は経験的な値や文献値などにより決定される。しかしながら、伝播速度には個体差があり、また、伝播速度は被検体の保持方法などにより変化する。そのため、画像再構成時に用いた速度と実際の伝播速度が離れていた場合、画像の解像度が大きく劣化する課題があった。 In such a photoacoustic tomography apparatus and an ultrasonic diagnostic apparatus conventionally used in the medical field (an apparatus that detects an acoustic wave reflected in a living body and reconstructs an image), the average of the subject is usually used. An image is reconstructed using a sound velocity (acoustic wave propagation velocity inside the subject). In general, the propagation speed is determined by an empirical value or a literature value. However, there are individual differences in the propagation speed, and the propagation speed changes depending on the method of holding the subject. Therefore, when the speed used at the time of image reconstruction is different from the actual propagation speed, there is a problem that the resolution of the image is greatly deteriorated.
このような課題の解決方法は、例えば特許文献1に開示されている。特許文献1に開示の技術では、各画素(ボクセルまたはピクセル)の輝度値又はコントラストが最大になるように伝播速度が決定される。それにより、画像再構成時に用いた速度と実際の伝播速度との不一致による画質劣化を抑制している。
A solution to such a problem is disclosed in
しかしながら、特許文献1に開示の技術では、画素ごとの輝度値又はコントラストを最大にするため、バックグラウンドノイズの輝度値やコントラストも上昇させてしまう。さらに、検出信号にノイズが含まれている場合には、ノイズ成分と正常な信号成分の合計値が最大となるように伝播速度が決定されるため、正確な伝播速度を得ることができず、画像にボケが表れてしまう。
However, in the technique disclosed in
そこで、本発明は、ノイズおよび画像再構成時に用いた速度と実際の伝播速度との不一致による画質劣化が抑制された、被検体内部の情報を表す画像を再構成することのできる画像生成装置及び画像生成方法を提供することを目的とする。 Accordingly, the present invention provides an image generation apparatus capable of reconstructing an image representing information inside a subject, in which image quality deterioration due to mismatch between noise and the speed used at the time of image reconstruction and the actual propagation speed is suppressed, and An object is to provide an image generation method.
本発明の画像生成装置は、被検体内部から伝播してくる音響波を検出して、検出信号を出力する複数の検出素子を有する探触子と、被検体内部における音響波の伝播速度を決定
する決定手段と、決定手段で決定された伝播速度と、複数の検出素子から得られる複数の検出信号とを用いて被検体内部の情報を表す画像を再構成する画像生成手段と、を有し、決定手段は、被検体内の一部分に注目領域を設定し、各検出素子から注目領域までの距離と、仮定の速度とに基づいて各検出信号の位相を調整し、該位相が調整された複数の検出信号の強度のばらつきを算出する処理を、複数の仮定の速度に対して実行し、複数の仮定の速度のうち、強度のばらつきが最小となる速度を伝播速度に決定することを特徴とする。
The image generation apparatus of the present invention detects an acoustic wave propagating from the inside of a subject, determines a probe having a plurality of detection elements that output detection signals, and the propagation speed of the acoustic wave inside the subject. Determination means, and an image generation means for reconstructing an image representing information inside the subject using the propagation speed determined by the determination means and a plurality of detection signals obtained from a plurality of detection elements. The determination unit sets a region of interest in a part of the subject, adjusts the phase of each detection signal based on the distance from each detection element to the region of interest and the assumed speed, and the phase is adjusted. The process of calculating the intensity variation of a plurality of detection signals is executed for a plurality of hypothetical speeds, and the speed at which the intensity variation is minimum among the plurality of hypothetical speeds is determined as a propagation speed. And
本発明の画像生成方法は、複数の検出素子のそれぞれを用いて被検体内部から伝播してくる音響波を検出して、検出信号を生成するステップと、被検体内部における音響波の伝播速度を決定する決定ステップと、決定ステップで決定された伝播速度と、複数の検出素子から得られる複数の検出信号とを用いて被検体内部の情報を表す画像を再構成する画像生成ステップと、を有し、決定ステップでは、被検体内の一部分に注目領域を設定し、各検出素子から注目領域までの距離と、仮定の速度とに基づいて各検出信号の位相を調整し、該位相が調整された複数の検出信号の強度のばらつきを算出する処理を、複数の仮定の速度に対して実行し、複数の仮定の速度のうち、強度のばらつきが最小となる速度を伝播速度に決定することを特徴とする。 The image generation method of the present invention detects the acoustic wave propagating from the inside of the subject using each of the plurality of detection elements, generates a detection signal, and determines the propagation speed of the acoustic wave inside the subject. A determination step for determining, and an image generation step for reconstructing an image representing information inside the subject using the propagation velocity determined in the determination step and a plurality of detection signals obtained from a plurality of detection elements. In the determination step, a region of interest is set in a part of the subject, and the phase of each detection signal is adjusted based on the distance from each detection element to the region of interest and the assumed speed, and the phase is adjusted. The process of calculating the intensity variation of the plurality of detection signals is executed for a plurality of hypothetical speeds, and the speed at which the intensity variation is minimum is determined as the propagation speed among the plurality of hypothetical speeds. Features .
本発明のプログラムは、コンピュータに、被検体内部における音響波の伝播速度を決定する決定ステップと、複数の検出素子が被検体内部から伝播してくる音響波を検出し、生成した複数の検出信号と、決定ステップで決定された伝播速度とを用いて被検体内部の情報を表す画像を再構成する画像生成ステップと、を実行させるためのプログラムであって、決定ステップでは、被検体内の一部分に注目領域を設定し、各検出素子から注目領域までの距離と、仮定の速度とに基づいて各検出信号の位相を調整し、該位相が調整された複数の検出信号の強度のばらつきを算出する処理を、複数の仮定の速度に対して実行し、複数の仮定の速度のうち、強度のばらつきが最小となる速度を伝播速度に決定することを特徴とする。 The program of the present invention has a determination step for determining a propagation speed of an acoustic wave inside a subject, and a plurality of detection signals generated by detecting acoustic waves propagating from the inside of the subject by a plurality of detection elements. And an image generation step for reconstructing an image representing information inside the subject using the propagation velocity determined in the determination step, wherein the determination step includes a part of the object in the subject. Set the region of interest, adjust the phase of each detection signal based on the distance from each detection element to the region of interest and the assumed speed, and calculate the intensity variation of multiple detection signals that have been adjusted This processing is executed for a plurality of hypothetical speeds, and the speed at which the intensity variation is minimum is determined as the propagation speed among the plurality of hypothetical speeds.
本発明によれば、ノイズおよび画像再構成時に用いた速度と実際の伝播速度との不一致による画質劣化が抑制された、被検体内部の情報を表す画像を再構成することができる。 According to the present invention, it is possible to reconstruct an image representing information inside the subject in which image quality deterioration due to mismatch between the noise and the speed used at the time of image reconstruction and the actual propagation speed is suppressed.
<原理>
以下、本発明の原理について説明する。
図1は本発明の原理を示す模式図である。図1において、符号1は(被検体内部の)音源、符号2はボクセル(あるいはピクセル)、符号3は検出素子を示す。各検出素子は、時間分割測定により被検体内部から伝播してくる音響波を検出して、検出信号を出力する。なお、検出素子3の数はN(Nは2以上の整数)個とする。
一般的に、音響波(超音波)を利用した画像生成装置では、複数の検出素子iから得られる複数の検出信号S(i,t)を用いて被検体内部の情報を表す画像を再構成する。こ
のような画像は再構成画像と呼ばれ、再構成画像が3次元画像(ボクセルの集合体)の場合にはボリュームデータとも呼ばれる。再構成画像が2次元画像の場合にはピクセルデータとも呼ばれる。再構成画像の各画素(ボクセルあるいはピクセル)は、一般的には各検出素子からその画素に対応する位置までの距離と伝播速度(被検体内部における音響波の伝播速度)に基づいて位相が調整された検出信号から算出される。ただし、フーリエドメイン法などでは、周波数空間での演算により、画像が再構成される。なお、図1においてiは検出素子の番号(0からN−1の範囲の整数)、tは時間を示す。
以下、一般的な画像再構成法であるタイムドメイン法について具体的に説明する。まず、検出素子毎に、i番目の検出素子から再構成画像内の或る画素に対応する位置までの距離を伝播速度で除算する。それにより、その画素位置で音響波が発生した場合において、該音響波が検出素子iに到達するまでの時間τi(遅延時間)を算出する(発生時刻をt=0とする)。そして、検出素子毎に、時間τiのときの検出信号の強度S(i,τi)を算出し、それらを足し合わせることにより、再構成画像の画素を生成する(同様に複数の位置について画素を生成することにより、再構成画像が生成される)。このような再構成画像の生成方法において、各ボクセル(またはピクセル)のデータが最大になるように伝播速度、すなわち時間τiを決定するのが特許文献1に開示の技術である。なお、強度S(i,τi)は時間τiでの検出信号S(i,t)の強度(位相が調整された検出信号の強度)を示す。
<Principle>
Hereinafter, the principle of the present invention will be described.
FIG. 1 is a schematic diagram showing the principle of the present invention. In FIG. 1,
In general, in an image generation apparatus using acoustic waves (ultrasound), an image representing information inside a subject is reconstructed using a plurality of detection signals S (i, t) obtained from a plurality of detection elements i. To do. Such an image is called a reconstructed image. When the reconstructed image is a three-dimensional image (an aggregate of voxels), it is also called volume data. When the reconstructed image is a two-dimensional image, it is also called pixel data. The phase of each pixel (voxel or pixel) of the reconstructed image is generally adjusted based on the distance from each detection element to the position corresponding to that pixel and the propagation speed (acoustic wave propagation speed inside the subject). Calculated from the detected signal. However, in the Fourier domain method or the like, an image is reconstructed by calculation in a frequency space. In FIG. 1, i represents the number of the detection element (an integer in the range of 0 to N-1), and t represents time.
Hereinafter, a time domain method that is a general image reconstruction method will be described in detail. First, for each detection element, the distance from the i-th detection element to a position corresponding to a certain pixel in the reconstructed image is divided by the propagation speed. Thereby, when an acoustic wave is generated at the pixel position, a time τi (delay time) until the acoustic wave reaches the detection element i is calculated (the generation time is set to t = 0). Then, for each detection element, the intensity S (i, τi) of the detection signal at time τi is calculated and added together to generate a pixel of the reconstructed image (similarly, the pixels at a plurality of positions are selected). By generating, a reconstructed image is generated). In such a reconstructed image generation method, the technique disclosed in
ここで、伝播速度として予想した速度と、実際の伝播速度とが大きく異なっていた場合には、強度S(i,τi)は大きくばらつく。そこで、本発明では、被検体内の一部分に注目領域5(注目すべき領域)を設定する。その後、各検出素子iから注目領域5までの距離と、仮定の速度とに基づいて各検出信号S(i,t)の位相を調整し、該位相が調整された複数の検出信号の強度S(i,τi)のばらつきを算出する処理を、複数の仮定の速度に対して実行する。そして、複数の仮定の速度のうち、上記強度のばらつきS(i,τi)が最小となる速度を伝播速度に決定し、該伝播速度に基づいて再構成画像を生成する。
Here, when the speed estimated as the propagation speed is significantly different from the actual propagation speed, the strength S (i, τi) varies greatly. Therefore, in the present invention, a region of interest 5 (region to be noted) is set in a part of the subject. Thereafter, the phase of each detection signal S (i, t) is adjusted based on the distance from each detection element i to the region of
例えば、遅延信号値のばらつきの指標として、以下の式(101)で表されるコヒーレントファクター(CF)を用いればよい。
ここで、S(i)は、検出素子iから出力され、位相が調整された検出信号の強度、即ち、ある速度に基づいて算出された時間τiでの検出信号の強度S(i,τi)である。S(i)が検出素子によらず同じ値となる場合、指標CF=1となり、S(i)の値のばらつきが大きいほど、指標CFは0に近づく。つまり、強度S(i,τi)のばらつきが大きいほど、指標CFは0に近い値となり、ばらつきが小さいほど1に近い値となる。本発明では、強度S(i,τi)のばらつきが最小、すなわち指標CFが最大となる速度を伝播速度に決定する。なお、指標はCFに限らない。分散値や標準偏差値など統計分野で使われる尺度を用いてばらつきを評価しても良い。
For example, a coherent factor (CF) represented by the following equation (101) may be used as an index of variation in the delay signal value.
Here, S (i) is the intensity of the detection signal output from the detection element i and whose phase is adjusted, that is, the intensity S (i, τi) of the detection signal at time τi calculated based on a certain speed. It is. When S (i) has the same value regardless of the detection element, the index CF = 1, and the greater the variation in the value of S (i), the closer the index CF is to 0. That is, the larger the variation of the intensity S (i, τi), the closer the index CF is to 0, and the smaller the variation is, the closer to 1 is. In the present invention, the speed at which the variation of the intensity S (i, τi) is minimum, that is, the index CF is maximum is determined as the propagation speed. The index is not limited to CF. The variation may be evaluated using a scale used in the statistical field such as a variance value or a standard deviation value.
このように、本発明では、強度S(i,τi)のばらつきが最小となる速度を伝播速度とすることにより、実際の伝播速度に近い速度に基づいて再構成画像を生成することが可
能となる。また、本発明では、強度S(i,τi)のばらつきを最小とする領域を制限することにより、バックグラウンドノイズの増加を抑制することができる。また、本発明では、強度S(i,τi)のばらつきが最小となる速度を伝播速度とするため、検出信号にノイズがランダム含まれていたとしても、そのノイズを平均化することができる(即ち、ノイズの影響を小さくすることができる)。その結果、検出信号にノイズが含まれていても、実際の伝播速度を精度良く推測することができる。
そのため、本発明によれば、ノイズおよび画像生成時に用いた速度と実際の伝播速度との不一致による画質劣化が抑制された、被検体内部の情報を表す画像を再構成することができる。
As described above, in the present invention, the speed at which the variation in the intensity S (i, τi) is minimized is set as the propagation speed, so that a reconstructed image can be generated based on a speed close to the actual propagation speed. Become. Further, in the present invention, an increase in background noise can be suppressed by limiting a region that minimizes variations in intensity S (i, τi). Further, in the present invention, since the speed at which the variation of the intensity S (i, τi) is minimized is the propagation speed, even if noise is randomly included in the detection signal, the noise can be averaged ( That is, the influence of noise can be reduced). As a result, even if noise is included in the detection signal, the actual propagation speed can be estimated with high accuracy.
Therefore, according to the present invention, it is possible to reconstruct an image representing information inside the subject in which image quality deterioration due to mismatch between the noise and the speed used when generating the image and the actual propagation speed is suppressed.
<実施形態>
次に、本実施形態に係る画像生成装置及び画像生成方法について説明する。図2は、本実施形態に係る画像生成装置の構成の一例を示す。ここでは、被検体内部から伝播してくる音響波を検出し、生体情報を可視化する光音響トモグラフィー技術を用いた画像生成装置(光音響トモグラフィー装置)の場合を例にして説明する。なお、本発明は音響波探触子から音響波(超音波)を送信し、被検体内で反射した反射音響波(超音波)を検出することで被検体内部を表す情報を画像化する超音波診断技術を用いた画像生成装置(超音波診断装置)にも適用できる。
<Embodiment>
Next, an image generation apparatus and an image generation method according to the present embodiment will be described. FIG. 2 shows an example of the configuration of the image generation apparatus according to this embodiment. Here, a case of an image generation apparatus (photoacoustic tomography apparatus) using a photoacoustic tomography technique that detects acoustic waves propagating from the inside of a subject and visualizes biological information will be described as an example. In the present invention, an acoustic wave (ultrasonic wave) is transmitted from an acoustic wave probe, and a reflected acoustic wave (ultrasonic wave) reflected in the subject is detected to image information representing the inside of the subject. The present invention can also be applied to an image generation apparatus (ultrasonic diagnostic apparatus) using an ultrasonic diagnostic technique.
本実施形態に係る画像生成装置は、大まかに、光源11、光学機器13、音響波探触子17、信号収集器18、画像再構成処理部19、信号処理器20、表示装置21から構成される。光源11から発せられた光12(強度変調光)は、光学機器13を介して生体などの被検体15に照射される。被検体15の内部を伝播した光のエネルギーの一部が血管などの光吸収体(結果的に音源となる)14に吸収されると、その光吸収体14の熱膨張により音響波(典型的には超音波)16が発生する。その発生した音響波16は音響波探触子17により検出され、その後の処理により被検体の生体情報が画像化される。
The image generating apparatus according to the present embodiment is roughly composed of a
光源11は、生体を構成する成分のうち特定の成分に吸収される特定の波長の光を発する。なお、光源11は、本実施形態に係る画像生成装置と一体であってもよいし、別体であってもよい。光源11は、数ナノから数百ナノ秒オーダーのパルス光を発生可能な1つ以上のパルス光源で構成される。光源としては大きな出力が得られるレーザーが好ましいが、レーザーのかわりに発光ダイオードなどを用いてもよい。レーザーとしては、固体レーザー、ガスレーザー、色素レーザー、半導体レーザーなど様々なレーザーを使用することができる。照射のタイミング、光の波形、強度などは不図示の制御部によって制御される。
The
光源11から放たれた光12は、光学機器13を介して被検体に導かれる(光導波路などを用いて導かれてもよい)。光学機器13は、例えば、光を反射するミラーや、光を集光したり拡大したりするレンズなどである。このような光学機器は、光源から発せられた光12が被検体15に所望の形状で照射されれば、どのようなものを用いてもかまわない。なお、一般的には、光12をレンズで集光させるより、ある程度の面積まで広げる方が、生体への安全性ならびに診断領域を広げられるという観点で好ましい。また、被検体の光が照射される領域(照射領域)は移動可能(変更可能)であることが好ましい。言い換えると、本実施形態に係る画像生成装置は、光源から発生した光が被検体上を移動可能となるように構成されていることが好ましい。それにより、より広範囲の再構成画像を生成することができる。また、照射領域(被検体に照射される光)は、音響波探触子17と同期して移動可能であることが好ましい。照射領域を移動させる方法としては、可動式ミラー等を用いる方法、光源自体を機械的に移動させる方法などがある。
The light 12 emitted from the
本実施形態に係る画像生成装置は、例えば、人や動物の悪性腫瘍や血管疾患などの診断
に利用される。そのため、被検体15としては、人体や動物の乳房や指・手足などの診断の対象部位が想定される。光音響トモグラフィー装置では、光吸収体(あるいは音源)14は、被検体内の吸収係数が高い部分である。例えば、人体が測定対象である場合には、酸化ヘモグロビン、還元ヘモグロビン、それらを多く含む血管、及び、新生血管を多く含む悪性腫瘍などが光吸収体となる。また、超音波診断装置では、音源14は、周囲との音響インピーダンスが異なる組織界面である。
The image generation apparatus according to the present embodiment is used, for example, for diagnosis of a human or animal malignant tumor or vascular disease. Therefore, the subject 15 is assumed to be a target site for diagnosis such as a breast of a human body or an animal, a finger, a limb or the like. In the photoacoustic tomography apparatus, the light absorber (or sound source) 14 is a portion having a high absorption coefficient in the subject. For example, when the human body is a measurement target, oxygenated hemoglobin, reduced hemoglobin, blood vessels containing many of them, malignant tumors containing many new blood vessels, and the like are light absorbers. In the ultrasonic diagnostic apparatus, the
音響波探触子17は、複数の検出素子を有する。検出素子は、圧電現象を利用したトランスデューサー、光の共振を利用したトランスデューサー、または、容量の変化を利用したトランスデューサーなどで構成される。但し、音響波を検知することができれば、このような構成に限らない。音響波探触子17において、典型的には、複数の検出素子は1次元的あるいは2次元的に配置される。このように配置された複数の検出素子を用いることにより、一度に複数の場所で音響波を検出することができる。そのため、検出時間を短縮できると共に、被検体の振動などの影響を低減できる。なお、1つの検出素子の位置を移動させて、複数箇所で音響波を検出すれば、複数の検出素子を用いた場合と同様の信号(検出信号)を得ることができる。
The
信号収集器18は、音響波探触子17(各検出素子)より得られた電気信号(アナログ信号)を増幅し、デジタル信号(検出信号)に変換する。信号収集器18は、典型的には、増幅器、A/D変換器、FPGA(Field Programmable Gate
Array)チップなどで構成される。なお、複数の検出信号を同時に処理できることが望ましい。それにより、画像を生成するまでの時間を短縮できる。
The
Array) chip. It is desirable that a plurality of detection signals can be processed simultaneously. Thereby, the time until the image is generated can be shortened.
画像再構成処理部19は、複数の検出素子から得られる複数の検出信号(信号収集器18から出力される複数のデジタル信号)を用いて被検体内部の情報を表す画像(再構成画像)を再構成する(画像生成手段)。再構成画像の各画素は、一般的なタイムドメイン法では、各検出素子からその画素に対応する位置までの距離と後述する信号処理器20で決定された伝播速度に基づいて位相が調整された検出信号から算出される。
画像再構成処理部19に接続される信号処理器20は、被検体内部における音響波の伝播速度を決定する(決定手段)。本実施形態では、信号収集器18から出力される複数のデジタル信号と画像再構成処理部19で得られた画像を用いて伝播速度を決定する。
The image
The
表示装置21は画像再構成処理部19で生成された再構成画像を表示する装置である。表示装置21としては、液晶ディスプレイ、プラズマディスプレイ、有機ELディスプレイ、電子放出素子を有するディスプレイなどを利用することができる。
The
次に、画像再構成処理部19と信号処理器20による処理について、図3のフローチャートを用いて説明する。
まず、画像再構成処理部19が、被検体の予想される伝播速度(所定の速度)に基づいて仮の再構成画像(ボリュームデータ)を生成する(ステップ301)。仮の再構成画像の生成方法としては、光音響トモグラフィー装置や超音波診断装置で従来から用いられている手法を用いることができる。例えば、タイムドメインあるいはフーリエドメインでの逆投影法などである。
なお、ステップ301では、超音波診断装置や光音響トモグラフィー装置などのような、音響波を検出して再構成画像を生成する画像生成装置を用いることが好ましい。但し、これに限らず、X線CTやMRIなど、別原理の画像生成装置を用いて、生体情報にかかわるボリュームデータ(再構成画像)を生成してもよい。
Next, processing by the image
First, the image
In
次に、信号処理器20が、ステップ301で得られたボリュームデータ(再構成画像)から生体情報を強く反映した画素(注目画素)を選択し、該注目画素に対応する位置を含
む領域を注目領域に決定する(ステップ302)。注目画素は、光音響トモグラフィー装置の場合、吸収体を構成する画素であり、超音波再構成画像の場合、超音波の反射が大きい領域の画素である。
Next, the
注目画素は、ユーザが再構成画像を確認することにより手動で選択されてもよいし、自動で選択されてもよい。注目画素を自動で選択する場合には、例えば、画像内において、輝度またはコントラストが最も高い画素を注目画素として選択すればよい。 The target pixel may be manually selected by the user confirming the reconstructed image, or may be automatically selected. When the target pixel is automatically selected, for example, a pixel having the highest luminance or contrast in the image may be selected as the target pixel.
注目領域をどの程度の大きさとするかは、後述する仮定の速度の選択範囲に基づいて決定される。例えば、1500m/秒の速度(予想される伝播速度c0;所定の速度)を用いて再構成画像を得た後に、速度を1400〜1600m/秒の範囲内の他の速度(仮定の速度)に切り替え、画像を生成しなおす場合を考える。なお、画素間隔dを0.25mm、被検体の厚さr(注目画素に対応する被検体内の位置から音響波探触子までの距離)を40mmとする。その場合、用いる速度を変化させることにより、注目画素に対応する位置が最大で前後上下左右に約10画素分ずれる可能性がある。そのため、そのような場合には、注目画素の前後上下左右それぞれ10画素を含む範囲(例えば、注目画素を中心として21×21×21個の画素を含む範囲)に対応する領域を注目領域とすればよい。 The size of the attention area is determined based on a hypothetical speed selection range described later. For example, after obtaining a reconstructed image using a velocity of 1500 m / sec (expected propagation velocity c 0 ; predetermined velocity), the velocity is set to another velocity (assumed velocity) within the range of 1400 to 1600 m / sec. Consider the case of switching to, and regenerating the image. Note that the pixel interval d is 0.25 mm, and the thickness r of the subject (the distance from the position in the subject corresponding to the target pixel to the acoustic probe) is 40 mm. In that case, by changing the speed to be used, there is a possibility that the position corresponding to the target pixel is shifted by about 10 pixels in the front, rear, up, down, left, and right directions. Therefore, in such a case, a region corresponding to a range including 10 pixels before, after, above, below, left, and right of the target pixel (for example, a range including 21 × 21 × 21 pixels centered on the target pixel) is set as the target region. That's fine.
具体的には、予想される伝播速度c0(所定の速度)を用いた場合の注目画素に対応する位置は、用いる速度を変化させた場合に、該位置を中心として約幅M画素分だけ変化する。Mの値は以下の式(102)により得られる。なお、式(102)において、cminは仮定の速度のうちの最小値、cmaxは仮定の速度のうちの最大値である。
なお、この式の被検体の厚さrは、具体的には、注目画素に対応する位置と該位置から最も遠い検出素子までの距離である。但し、これに限らず、注目画素に対応する位置と該位置から最も近い検出素子までの距離を被検体の厚さとしても良い。注目画素に対応する位置と各検出素子との間の平均距離を被検体の厚さとしても良い。 The thickness r of the subject in this equation is specifically the distance from the position corresponding to the target pixel and the detection element farthest from the position. However, the present invention is not limited thereto, and the distance from the position corresponding to the target pixel and the closest detection element to the position may be the thickness of the subject. The average distance between the position corresponding to the target pixel and each detection element may be the thickness of the subject.
次に、信号処理器20が、各検出素子から注目領域(例えば、再構成画像の画素に対応する注目領域内の位置)までの距離と、仮定の速度とに基づいて各検出信号の位相を調整する。そして、該位相が調整された複数の検出信号の強度のばらつきを算出する(ステップ303)。なお、ステップ303の処理は、複数の仮定の速度に対して実行される。
そして、信号処理器20が、複数の仮定の速度のうち、強度のばらつきが最小となる速度を伝播速度に決定する(ステップ304)。また、画像再構成処理部19が、該決定された伝播速度に基づいて再構成画像を生成する。
以上の処理により、ノイズおよび画像再構成時に用いた速度と実際の伝播速度との不一致による画質劣化が抑制された、被検体内部の情報を表す画像を再構成することができる。
Next, the
Then, the
Through the above processing, it is possible to reconstruct an image representing information inside the subject in which image quality deterioration due to mismatch between the noise and the speed used at the time of image reconstruction and the actual propagation speed is suppressed.
なお、検出信号は検出素子の面積が大きい場合には指向性の影響を強く受ける。指向性について図4を用いて説明する。指向性とは検出素子と被検体内部の位置(画素に対応する位置;処理対象位置)との相対的な位置関係により、検出素子の検出感度が変化する特性のことである。例えば、図4(A)のように検出素子が処理対象位置の方向を向いている場合には、該位置からの音響波の検出感度が最大となる(検出信号の強度は最大となる)。一方、検出素子の向きと処理対象位置とのずれが大きいほど、検出感度は低下する(検出信号の強度が低下する)。このような指向性の影響は、素子幅が大きい場合や、検出する音響波の周波数が大きい場合に顕著に現れる。指向性R(検出感度)は式(103)で表される。
ここで、検出素子の形状は正方形とする。式(103)において、dは検出素子の幅、θは処理対象位置と検出素子を結ぶ直線と、検出素子の向きとの間の角度、kは2π/λ、λは音響波の波長である。
The detection signal is strongly influenced by directivity when the area of the detection element is large. The directivity will be described with reference to FIG. The directivity is a characteristic in which the detection sensitivity of the detection element changes depending on the relative positional relationship between the detection element and a position inside the subject (a position corresponding to a pixel; a processing target position). For example, when the detection element faces the direction of the processing target position as shown in FIG. 4A, the detection sensitivity of the acoustic wave from the position is maximized (the intensity of the detection signal is maximized). On the other hand, the greater the deviation between the direction of the detection element and the position to be processed, the lower the detection sensitivity (the detection signal strength decreases). Such influence of directivity appears prominently when the element width is large or when the frequency of the acoustic wave to be detected is large. The directivity R (detection sensitivity) is expressed by Expression (103).
Here, the detection element has a square shape. In Expression (103), d is the width of the detection element, θ is the angle between the straight line connecting the processing target position and the detection element, and the direction of the detection element, k is 2π / λ, and λ is the wavelength of the acoustic wave. .
図4(B)に素子幅2mmの検出素子の指向性の一例を示す。ここでは、検出素子のインパルス応答から得られる周波数特性を用いて、各周波数成分の大きさを算出し、上記の式(103)から指向性(最大値が1で規格化された検出感度)を計算している。この図から分かるように、角度θが大きくなると急激に感度(検出信号の強度)が低下することがわかる。このような指向性があると、検出信号間のばらつきが大きくなる。 FIG. 4B shows an example of directivity of a detection element having an element width of 2 mm. Here, the magnitude of each frequency component is calculated using the frequency characteristics obtained from the impulse response of the detection element, and the directivity (detection sensitivity normalized with a maximum value of 1) is calculated from the above equation (103). I'm calculating. As can be seen from this figure, the sensitivity (the intensity of the detection signal) suddenly decreases as the angle θ increases. Such directivity increases the variation between detection signals.
そのため、本実施形態に係る画像生成装置は、画素に対応する位置と、検出素子との相対的な位置関係に応じて、該検出素子の指向性による検出信号の強度の低下を補正する機能をさらに有することが好ましい。そして、信号処理器20が、該強度が補正された検出信号を用いて伝播速度を決定し、画像再構成処理部19が、該強度が補正された検出信号を用いて画像を再構成することが好ましい。それにより、検出素子の指向性による影響を抑制することができる。例えば、指向性による検出信号のばらつきが低減されるため(速度によるばらつきが支配的となるため)、より精度良く伝播速度を決定することができる。また、指向性による検出信号の強度のばらつきが低減されるため、より画質劣化の無い画像を生成することができる。
For this reason, the image generation apparatus according to the present embodiment has a function of correcting a decrease in the intensity of the detection signal due to the directivity of the detection element according to the relative positional relationship between the position corresponding to the pixel and the detection element. Furthermore, it is preferable to have. Then, the
例えば、処理対象位置と検出素子の相対的な位置関係から、角度θを計算し、その角度に応じた指向性の値(例えば、図4(B)の縦軸の値)に基づいて、該検出素子から出力された検出信号の強度を大きくする補正を行えばよい。具体的には、指向性の値の逆数を検出信号に乗算すればよい。それにより、指向性による検出信号の強度のばらつきを抑えることができる。 For example, the angle θ is calculated from the relative positional relationship between the processing target position and the detection element, and based on the directivity value corresponding to the angle (for example, the value on the vertical axis in FIG. 4B), Correction for increasing the intensity of the detection signal output from the detection element may be performed. Specifically, the detection signal may be multiplied by the reciprocal of the directivity value. Thereby, variation in the intensity of the detection signal due to directivity can be suppressed.
なお、多くの場合、画像再構成処理部19、信号処理器20の処理は、コンピュータ(ワークステーションなど)の演算処理装置がソフトウェア(プログラム)を実行することで実現される。そのため、画像再構成処理部19と信号処理器20を区別することができない場合がある。ただし、これらの機能要素はそれぞれ専用のチップで構成されていてもよい。
なお、画像再構成処理部19、信号処理器20の処理は、ソフトウェアではなく、ハードウェアで実現されていてもよい。
In many cases, the processing of the image
Note that the processing of the image
なお、本実施形態では、注目画素の位置と式(102)とを用いて注目領域を決定するものとしたが、この構成に限らない。予め被写体が決まっている場合など、注目領域の大まかな場所が予想できる場合には、注目領域は予め装置内に記憶されていてもよい。被写体の種類(測定部位)毎に予め決められた注目領域が記憶されていてもよい。そのような場合には、本実施形態のステップ301のように仮の再構成画像を生成する必要はない。また、注目領域はどのような大きさであってもよい(1画素分の大きさであってもよい)。例えば、ステップ302で選択された注目画素に対応する位置を含むように注目領域を設定すればよい。音源はある程度の大きさを有することが想定されるため、そのように注目領域を設定すれば、速度を変更することにより注目画素に対応する位置が注目領域の外側に位置することになったとしても、注目領域で音源を捉えることができる。
In this embodiment, the region of interest is determined using the position of the pixel of interest and Expression (102), but the present invention is not limited to this configuration. When a rough place of the attention area can be predicted, such as when a subject is determined in advance, the attention area may be stored in the apparatus in advance. A region of interest determined in advance for each type of subject (measurement site) may be stored. In such a case, it is not necessary to generate a temporary reconstructed image as in
<実施例1>
シミュレーションにより本発明を実施し、その効果を検証した。以下、図5,6を用いて説明する。
被検体として、図5(A)に示すように、X方向:4cm、Y方向:4cm、Z方向:4cmの立方体ファントム43を用いた。立方体ファントム43の中心には、直径が0.05cm、高さが2cmの円柱形状の音源41を、その軸方向がX方向と平行になるように設置した。ここで、音源41は、光音響トモグラフィーでは光吸収体、超音波診断装置では超音波反射体のことを示す。なお、被検体内部における音響波の伝播速度を1500m/秒、音源41で発生する初期発生音圧を3000Pa(パスカル)とした。
<Example 1>
The present invention was implemented by simulation and the effect was verified. This will be described below with reference to FIGS.
As a subject, as shown in FIG. 5A, a
音響波探触子42は2次元アレイ型で、20×20個の正方形型の検出素子で構成されているものとした。また、各素子の一辺の長さを2mm、ピッチ間隔を2mmとした。音響波は、サンプリング周波数20MHz、測定点数1280点で検出されるものとした。
The
以上の条件の元で、物理シミュレーションにより音響波探触子42の各素子で検出された音響波の検出信号を生成した。探触子の中心の検出素子から得られる検出信号の例を図5(B)に示す。図中の大きな信号が音源41からの信号を示している。なお、検出信号は、標準偏差の3倍が20Paとなるようなホワイトノイズを含むものとした。
Under the above conditions, acoustic wave detection signals detected by each element of the
このように得られた検出信号と、速度(所定の速度)1580m/秒を用いて再構成画像を生成した。生成方法としては、公知技術のタイムドメイン法であるユニバーサルバックプロジェクションを用いた。なお、画像化範囲はX方向:3.8cm、Y方向:3.8cm、Z方向:4cmとし、画素間隔(ボクセル間隔)は0.025cmとした。ボクセル総数は画像化範囲をボクセル間隔で除算したものであり、152×152×160となる。 A reconstructed image was generated using the detection signal thus obtained and a speed (predetermined speed) of 1580 m / sec. As a generation method, a universal back projection which is a known time domain method was used. The imaging range was X direction: 3.8 cm, Y direction: 3.8 cm, Z direction: 4 cm, and the pixel interval (voxel interval) was 0.025 cm. The total number of voxels is the imaging range divided by the voxel interval and is 152 × 152 × 160.
図6に再構成された画像を示す。図6に示す画像は、いずれも、音響波探触子42の方向に最大強度のボクセル値を投影したMIP(Maximum Intensity Projection)画像であり、白色領域が高い音圧を示している。また、各々の画像の上側には、その画像のY=2.0cmにおける、X方向の音圧変化を模式的に示している。
FIG. 6 shows the reconstructed image. Each of the images shown in FIG. 6 is a MIP (Maximum Intensity Projection) image obtained by projecting the maximum intensity voxel value in the direction of the
図6において、画像Aは被検体内の音源のMIP像を示している。本実施例では画像Aに近い画像ほどを解像度が高いものとする。
画像Bは、伝播速度として1580m/秒の速度を用いたときに得られる再構成画像である。このように実際の伝播速度1500m/秒と異なる速度を用いると、画像にボケが生じ、解像度が低下することが分かる。
In FIG. 6, an image A shows a MIP image of a sound source in the subject. In this embodiment, an image closer to the image A has a higher resolution.
Image B is a reconstructed image obtained when a velocity of 1580 m / sec is used as the propagation velocity. Thus, it can be seen that when a speed different from the actual propagation speed of 1500 m / sec is used, the image is blurred and the resolution is lowered.
次に、この画像(ボリュームデータ)において輝度が最大となるボクセルを探索した。
その結果、中心のボクセル(152×152×160の全ボクセルにおいて、76×76×80番目)が検出された。そして、このボクセルを中心として21×21×21ボクセルに対応する領域を注目領域として設定した。
Next, the voxel having the maximum luminance in this image (volume data) was searched.
As a result, the center voxel (76 × 76 × 80th in all 152 × 152 × 160 voxels) was detected. Then, an area corresponding to 21 × 21 × 21 voxels with this voxel as a center is set as a region of interest.
そして、従来技術との比較のために、用いる速度を1400m/秒から1600秒/秒まで2m/秒間隔で変化させて、注目領域内の位置(ボクセル)毎に輝度値を算出した。即ち、21×21×21×101個の輝度値を計算した。その結果、輝度値が最大となる速度は1456m/秒であった。そして、この速度1456m/秒を用いて、再び画像再構成を行った。画像Cはそのときに得られた画像である。多少の画像ボケはあるものの、画像Bと比較して、若干であるが解像度が向上していることが分かる。また、最大輝度の値も増加していることが分かる。 For comparison with the prior art, the speed used was changed from 1400 m / sec to 1600 sec / sec at 2 m / sec intervals, and the luminance value was calculated for each position (voxel) in the region of interest. That is, 21 × 21 × 21 × 101 luminance values were calculated. As a result, the speed at which the luminance value was maximized was 1456 m / sec. Then, image reconstruction was performed again using this speed of 1456 m / sec. Image C is an image obtained at that time. Although there are some image blurs, it can be seen that the resolution is slightly improved compared to the image B. It can also be seen that the value of the maximum luminance is increased.
次に、本発明の効果を実証するため、用いる速度を1400m/秒から1600秒/秒まで2m/秒間隔で変化させて、注目領域内の位置(ボクセル)毎に検出信号の強度のばらつき(指標CF)を算出した。即ち、21×21×21×101個の指標CFを計算した。その結果、指標CFが最大となる速度は1492m/秒であった。そして、この速度1492m/秒を用いて、再び画像再構成を行った。画像C’はそのときに得られた画像である。従来技術である画像Cと比較して、画像のボケが少なくなり、解像度が向上していることが分かる。なお、画像Cと比べると最大輝度値は低下している。
以上の結果を表1に示す。
Next, in order to demonstrate the effect of the present invention, the speed used is changed from 1400 m / sec to 1600 sec / sec at intervals of 2 m / sec, and the variation in the intensity of the detection signal for each position (voxel) in the region of interest ( The index CF) was calculated. That is, 21 × 21 × 21 × 101 index CFs were calculated. As a result, the speed at which the index CF was maximum was 1492 m / sec. Then, image reconstruction was performed again using this speed of 1492 m / sec. Image C ′ is an image obtained at that time. It can be seen that the image blur is reduced and the resolution is improved as compared with the image C which is the prior art. Note that the maximum luminance value is lower than that of the image C.
The results are shown in Table 1.
このように、本発明では、実際の伝播速度に近い速度(1492m/秒)を推測することができる。特に、ノイズが含まれた検出信号においては、従来技術で得られる速度(1456m/秒)よりもより実際の伝播速度に近い速度を推測することができる。その結果、解像度が低下することのない画像を得ることができる。
なお、実際の伝播速度が1500m/秒であるのに対し、本実施例では伝播速度が1492m/秒と推定されたのは、ノイズと検出素子の指向性の影響による。
なお、本実施例では、注目領域内の位置毎に検出信号の強度のばらつきを算出するものとしたが、注目領域内の一部の位置についてばらつきを算出する構成であってもよい。例えば、注目領域の中心位置のみについてのばらつきを算出する構成であってもよい。
Thus, in the present invention, a speed (1492 m / sec) close to the actual propagation speed can be estimated. In particular, in the detection signal including noise, it is possible to estimate a speed closer to the actual propagation speed than the speed (1456 m / sec) obtained by the conventional technique. As a result, it is possible to obtain an image whose resolution does not decrease.
The actual propagation speed is 1500 m / second, whereas in this embodiment, the propagation speed is estimated to be 1492 m / second because of the influence of noise and directivity of the detection element.
In this embodiment, the variation in the intensity of the detection signal is calculated for each position in the attention area. However, a configuration may be used in which the dispersion is calculated for some positions in the attention area. For example, it may be configured to calculate the variation only about the center position of the region of interest.
<実施例2>
次に、光音響トモグラフィー装置に本発明を適用した例について、図2を用いて説明する。本実施例では、光源11として波長1064nmで約10ナノ秒のパルス光を発生するQスイッチYAGレーザーを用いた。パルスレーザー光(光12)のエネルギーは0.6Jであり、そのパルス光をミラーとビームエキスパンダーなどの光学機器13を用いて、半径約2cm程度まで広げた。被検体15としては生体を模擬したファントムを用いた。ファントムは1%のイントラリピッドを寒天で固めたものを利用した。ファントムの大きさは幅:12cm、高さ:8cm、奥行き:4cmとした。また、このファントム内には直径0.03cmのゴム製の黒色ワイヤが中心付近に光吸収体14として埋め込まれている。このようなファントムに、半径2cmまで広げた光12を照射した。音響波探触子17としてはPZT(ジルコン酸チタン酸鉛)で作られた超音波トランスデューサーを用
いた。このトランスデューサーは2次元アレイ型で、素子数は18×18個、各素子の形状は正方形形状、素子ピッチは2mmである。また、素子の幅は約2mmである。なお、この音響波探触子は光照射領域と同期して、XY方向に移動可能になっており、大きな領域を画像化できるように構成されている。
<Example 2>
Next, an example in which the present invention is applied to a photoacoustic tomography apparatus will be described with reference to FIG. In this embodiment, a Q-switched YAG laser that generates pulsed light of about 10 nanoseconds at a wavelength of 1064 nm is used as the
このファントムの1つの面にパルス光を照射すると、ファントム内に拡散した光がゴム製ワイヤに吸収され、光音響波が発生する。その光音響波を324チャンネルの各超音波トランスデューサーで検出し、アンプ、A/Dコンバーター、FPGAからなる信号収集器18を用いて、チャンネル毎の光音響波のデジタルデータ(検出信号)を取得する。なお、信号のS/N比を向上させるために、36回レーザーを照射し、得られたすべての検出信号を平均化した。その後、得られたデジタルデータを画像再構成処理器19及び信号処理器20であるワークステーション(WS)へ転送し、WS内に保存する。この保存したデータに対して、離散ウェーブレット変換によるノイズ低減処理を行った後に、生体内部における音響波の平均的な伝播速度である1540m/秒の速度を用いて、画像再構成を行った。ここではタイムドメイン方式であるユニバーサルバックプロジェクション法を用いてボリュームデータを生成した。このとき使用したボクセル間隔は0.05cmとした。画像化範囲は11.8cm×11.8cm×4.0cmとした。得られた画像の一例を図7(A)に示す。
When one surface of the phantom is irradiated with pulsed light, the light diffused in the phantom is absorbed by the rubber wire and a photoacoustic wave is generated. The photoacoustic wave is detected by each of the ultrasonic transducers of 324 channels, and the digital data (detection signal) of the photoacoustic wave for each channel is obtained using the
次に、用いる速度を1400m/秒から1600秒/秒まで2m/秒間隔で変化させて、注目領域内の位置(ボクセル)毎に検出信号の強度のばらつき(指標CF)を算出した。なお、本実施例では、得られたボリュームデータにおいて、ワイヤで発生した音響波(初期音圧)を含むボクセル群(21×21×21個)に対応する領域を注目領域とした。また、本実施例では、検出素子のインパルス応答から算出した指向性を考慮して検出信号を補正し、補正した検出信号を用いて指標CFを計算した。最大となる指標CFは、指向性を考慮することで、0.056から0.078まで向上した。本実施例では、指標CFが最大となる速度は1454m/秒であった。そして、この速度1454m/秒と、PCに保存した検出信号とを用いて再び画像再構成を行った。得られた画像を表示装置21である液晶ディスプレイに表示した。得られた画像の一例を図7(B)に示す。
Next, the used speed was changed from 1400 m / sec to 1600 sec / sec at an interval of 2 m / sec, and the variation (index CF) in the intensity of the detection signal for each position (voxel) in the region of interest was calculated. In the present embodiment, in the obtained volume data, an area corresponding to a voxel group (21 × 21 × 21) including an acoustic wave (initial sound pressure) generated by a wire is set as an attention area. In this embodiment, the detection signal is corrected in consideration of the directivity calculated from the impulse response of the detection element, and the index CF is calculated using the corrected detection signal. The maximum index CF is improved from 0.056 to 0.078 by considering directivity. In this example, the speed at which the index CF becomes maximum was 1454 m / sec. Then, image reconstruction was performed again using this speed 1454 m / sec and the detection signal stored in the PC. The obtained image was displayed on a liquid crystal display as the
以上の結果から、1540m/秒の速度を用いて生成した画像(図7(A))よりも、1454m/秒の速度を用いて生成した画像(図7(B))のほうが、ゴム製ワイヤから発生した音響波(初期音圧)の幅が明らかに狭くなっていることが分かる。また、画像のボケも少ない、つまり、解像度が向上していることがわかる。
このように、本実施例では、被検体内部における音響波の伝播速度を正確に推定することができ、ひいては、再構成画像の解像度を向上させることができる。
From the above results, the image generated using the speed 1454 m / second (FIG. 7B) is more likely to be the rubber wire than the image generated using the speed 1540 m / second (FIG. 7A). It can be seen that the width of the acoustic wave (initial sound pressure) generated from is clearly narrowed. It can also be seen that the image is less blurred, that is, the resolution is improved.
As described above, in this embodiment, it is possible to accurately estimate the propagation speed of the acoustic wave inside the subject, and thus the resolution of the reconstructed image can be improved.
3 検出素子
5 注目領域
15 被検体
16 音響波
17 音響波探触子
19 画像再構成処理部
20 信号処理器
DESCRIPTION OF
本発明に係る画像生成装置は、被検体への光照射に起因して発生する音響波に基づいて光音響画像を生成する装置であって、前記音響波の伝播速度を決定する決定手段と、前記伝播速度と、前記音響波を検出することにより得られる信号とに基づいて、前記光音響画像を生成する画像生成手段と、を有し、前記決定手段は、超音波の送受信により得られる超音波画像に基づいて、注目位置を設定し、前記信号のうち、前記注目位置に対応する信号に基づいて、前記伝播速度を決定することを特徴とする。 Image generating apparatus according to the present invention is an apparatus for generating a photoacoustic image on the basis of acoustic waves caused by the light irradiation to a subject, determining means for determining a propagation speed before Kion HibikiNami If, before Kiden播速degree, on the basis of the signal obtained by detecting the acoustic waves, have a, and an image generating means for generating the photoacoustic image, said determining means, ultrasonic wave transmission and reception of A target position is set based on the ultrasonic image obtained by the above, and the propagation speed is determined based on a signal corresponding to the target position among the signals .
本発明に係る画像生成方法は、被検体への光照射に起因して発生する音響波に基づいて光音響画像を生成する画像生成方法であって、超音波の送受信により得られる超音波画像に基づいて、注目位置を設定する設定ステップと、前記音響波を検出することにより得られる信号のうち、前記注目位置に対応する信号に基づいて、前記音響波の伝播速度を決定する決定ステップと、前記信号及び前記伝播速度に基づいて、前記光音響画像を生成する生成ステップと、を含むことを特徴とする。 An image generation method according to the present invention is an image generation method for generating a photoacoustic image based on an acoustic wave generated due to light irradiation on a subject , and an ultrasonic image obtained by transmitting and receiving ultrasonic waves. A setting step for setting a position of interest, and a determination step for determining a propagation speed of the acoustic wave based on a signal corresponding to the position of interest among signals obtained by detecting the acoustic wave; Generating the photoacoustic image based on the signal and the propagation velocity .
本発明に係るプログラムは、被検体への光照射に起因して発生する音響波に基づいて光音響画像を生成する画像生成方法をコンピュータに実行させるプログラムであって、前記画像生成方法は、超音波の送受信により得られる超音波画像に基づいて、注目位置を設定する設定ステップと、前記音響波を検出することにより得られる信号のうち、前記注目位置に対応する信号に基づいて、前記音響波の伝播速度を決定する決定ステップと、前記信号及び前記伝播速度に基づいて、前記光音響画像を生成する生成ステップと、を含むことを特徴とする。 A program according to the present invention is a program for causing a computer to execute an image generation method for generating a photoacoustic image based on an acoustic wave generated due to light irradiation on a subject. A setting step for setting a position of interest based on an ultrasonic image obtained by transmission and reception of sound waves, and the acoustic wave based on a signal corresponding to the position of interest among signals obtained by detecting the acoustic wave. And a generation step of generating the photoacoustic image based on the signal and the propagation velocity.
Claims (9)
前記被検体内部における音響波の伝播速度を決定する決定手段と、
前記決定手段で決定された伝播速度と、複数の検出素子から得られる複数の検出信号とを用いて前記被検体内部の情報を表す画像を再構成する画像生成手段と、を有し、
前記決定手段は、
被検体内の一部分に注目領域を設定し、
各検出素子から前記注目領域までの距離と、仮定の速度とに基づいて各検出信号の位相を調整し、該位相が調整された複数の検出信号の強度のばらつきを算出する処理を、複数の仮定の速度に対して実行し、
前記複数の仮定の速度のうち、前記強度のばらつきが最小となる速度を前記伝播速度に決定することを特徴とする画像生成装置。 A probe having a plurality of detection elements for detecting acoustic waves propagating from within the subject and outputting detection signals;
Determining means for determining a propagation speed of the acoustic wave inside the object;
Image generation means for reconstructing an image representing information inside the subject using the propagation velocity determined by the determination means and a plurality of detection signals obtained from a plurality of detection elements;
The determining means includes
Set a region of interest in a part of the subject,
A process of adjusting the phase of each detection signal based on the distance from each detection element to the region of interest and an assumed speed, and calculating variations in the intensity of the plurality of detection signals whose phases are adjusted, Run for the assumed speed,
An image generating apparatus characterized in that, among the plurality of assumed speeds, a speed at which the intensity variation is minimum is determined as the propagation speed.
前記決定手段は、前記強度が補正された検出信号を用いて前記伝播速度を決定することを特徴とする請求項1に記載の画像生成装置。 According to the relative positional relationship between the position corresponding to the pixel and the detection element, further includes correction means for correcting a decrease in the intensity of the detection signal due to the directivity of the detection element;
The image generation apparatus according to claim 1, wherein the determination unit determines the propagation speed using a detection signal in which the intensity is corrected.
画素間隔をd、前記探触子から前記注目画素に対応する位置までの距離をr、前記所定の速度をc0とした場合に、前記決定手段は、前記注目画素の周りの、
When the pixel interval is d, the distance from the probe to the position corresponding to the target pixel is r, and the predetermined speed is c 0 , the determining means
前記被検体内部における音響波の伝播速度を決定する決定ステップと、
前記決定ステップで決定された伝播速度と、複数の検出素子から得られる複数の検出信号とを用いて前記被検体内部の情報を表す画像を再構成する画像生成ステップと、を有し、
前記決定ステップでは、
被検体内の一部分に注目領域を設定し、
各検出素子から前記注目領域までの距離と、仮定の速度とに基づいて各検出信号の位相を調整し、該位相が調整された複数の検出信号の強度のばらつきを算出する処理を、複数の仮定の速度に対して実行し、
前記複数の仮定の速度のうち、前記強度のばらつきが最小となる速度を前記伝播速度に決定することを特徴とする画像生成方法。 Detecting acoustic waves propagating from within the subject using each of the plurality of detection elements, and generating a detection signal;
A determination step for determining a propagation speed of an acoustic wave inside the subject;
An image generation step of reconstructing an image representing information inside the subject using the propagation velocity determined in the determination step and a plurality of detection signals obtained from a plurality of detection elements;
In the determination step,
Set a region of interest in a part of the subject,
A process of adjusting the phase of each detection signal based on the distance from each detection element to the region of interest and an assumed speed, and calculating variations in the intensity of the plurality of detection signals whose phases are adjusted, Run for the assumed speed,
An image generation method comprising: determining, as the propagation speed, a speed at which the intensity variation is minimum among the plurality of assumed speeds.
前記被検体内部における音響波の伝播速度を決定する決定ステップと、
複数の検出素子が被検体内部から伝播してくる音響波を検出し、生成した複数の検出信号と、前記決定手段で決定された伝播速度とを用いて前記被検体内部の情報を表す画像を再構成する画像生成ステップと、を実行させるためのプログラムであって、
前記決定ステップでは、
被検体内の一部分に注目領域を設定し、
各検出素子から前記注目領域までの距離と、仮定の速度とに基づいて各検出信号の位相を調整し、該位相が調整された複数の検出信号の強度のばらつきを算出する処理を、複数の仮定の速度に対して実行し、
前記複数の仮定の速度のうち、前記強度のばらつきが最小となる速度を前記伝播速度に決定することを特徴とするプログラム。 On the computer,
A determination step for determining a propagation speed of an acoustic wave inside the subject;
A plurality of detection elements detect acoustic waves propagating from the inside of the subject, and an image representing the information inside the subject is generated using the plurality of detection signals generated and the propagation speed determined by the determining unit. A program for executing an image generation step to be reconstructed,
In the determination step,
Set a region of interest in a part of the subject,
A process of adjusting the phase of each detection signal based on the distance from each detection element to the region of interest and an assumed speed, and calculating variations in the intensity of the plurality of detection signals whose phases are adjusted, Run for the assumed speed,
A program for determining, as the propagation speed, a speed at which the intensity variation is minimum among the plurality of assumed speeds.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2016100478A JP6113330B2 (en) | 2016-05-19 | 2016-05-19 | Apparatus and image generation method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2016100478A JP6113330B2 (en) | 2016-05-19 | 2016-05-19 | Apparatus and image generation method |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2014083739A Division JP5940109B2 (en) | 2014-04-15 | 2014-04-15 | Image generation apparatus, propagation speed determination method, and program |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2017050990A Division JP6381718B2 (en) | 2017-03-16 | 2017-03-16 | Apparatus and image generation method |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2016144754A true JP2016144754A (en) | 2016-08-12 |
JP6113330B2 JP6113330B2 (en) | 2017-04-12 |
Family
ID=56685174
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2016100478A Expired - Fee Related JP6113330B2 (en) | 2016-05-19 | 2016-05-19 | Apparatus and image generation method |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP6113330B2 (en) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007007045A (en) * | 2005-06-29 | 2007-01-18 | Toshiba Corp | Ultrasonograph |
JP2010512940A (en) * | 2006-12-21 | 2010-04-30 | ユニヴァシティート トウェンテ | Imaging apparatus and method |
-
2016
- 2016-05-19 JP JP2016100478A patent/JP6113330B2/en not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007007045A (en) * | 2005-06-29 | 2007-01-18 | Toshiba Corp | Ultrasonograph |
JP2010512940A (en) * | 2006-12-21 | 2010-04-30 | ユニヴァシティート トウェンテ | Imaging apparatus and method |
Also Published As
Publication number | Publication date |
---|---|
JP6113330B2 (en) | 2017-04-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5528083B2 (en) | Image generating apparatus, image generating method, and program | |
US10143381B2 (en) | Object information acquiring apparatus and control method therefor | |
JP6525565B2 (en) | Object information acquisition apparatus and object information acquisition method | |
JP6544910B2 (en) | INFORMATION PROCESSING APPARATUS, OBJECT INFORMATION ACQUIRING APPARATUS, AND METHOD OF DETERMINING SOUND SPEED | |
JP6478572B2 (en) | SUBJECT INFORMATION ACQUISITION DEVICE AND ACOUSTIC WAVE DEVICE CONTROL METHOD | |
EP3143391B1 (en) | Photoacoustic apparatus | |
JP2017119094A (en) | Information acquisition apparatus, information acquisition method, and program | |
JP6742734B2 (en) | Object information acquisition apparatus and signal processing method | |
JP6742745B2 (en) | Information acquisition device and display method | |
JPWO2018043193A1 (en) | Information acquisition apparatus and signal processing method | |
JP6656229B2 (en) | Photoacoustic device | |
JP5645637B2 (en) | Subject information acquisition apparatus and subject information acquisition method | |
US20180228377A1 (en) | Object information acquiring apparatus and display method | |
JP5940109B2 (en) | Image generation apparatus, propagation speed determination method, and program | |
JP6562800B2 (en) | Processing apparatus and processing method | |
JP6469133B2 (en) | Processing apparatus, photoacoustic apparatus, processing method, and program | |
JP6381718B2 (en) | Apparatus and image generation method | |
JP6113330B2 (en) | Apparatus and image generation method | |
US20200085345A1 (en) | Object information acquisition apparatus and method of controlling the same | |
JP6942847B2 (en) | Subject information acquisition device and signal processing method | |
JP2016013189A (en) | Photoacoustic apparatus | |
JP2017086173A (en) | Subject information acquisition device and control method thereof |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20160620 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20160620 |
|
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: 20170214 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20170215 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20170314 |
|
R151 | Written notification of patent or utility model registration |
Ref document number: 6113330 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |
|
S802 | Written request for registration of partial abandonment of right |
Free format text: JAPANESE INTERMEDIATE CODE: R311802 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
LAPS | Cancellation because of no payment of annual fees |