JP3763159B2

JP3763159B2 - SPECT absorption correction method

Info

Publication number: JP3763159B2
Application number: JP8743396A
Authority: JP
Inventors: 敏夫春日; 良彦熊澤; 廣巳田島
Original assignee: Shimadzu Corp
Current assignee: Shimadzu Corp
Priority date: 1996-03-15
Filing date: 1996-03-15
Publication date: 2006-04-05
Anticipated expiration: 2016-03-15
Also published as: JPH09251077A

Description

【０００１】
【発明の属する技術分野】
この発明は、核医学診断分野あるいは理工学分野で利用されるＳＰＥＣＴ装置（シングルフォトンエミッションＣＴ装置）に関し、とくに、その吸収補正方法に関する。
【０００２】
【従来の技術】
ＳＰＥＣＴ装置は、シングルフォトン放出性の放射性核種の放射性同位元素（ＲＩ）を用い、それから放出されるガンマ線を検出して核種の分布像を撮影するものである。たとえば人体にシングルフォトン放出性の放射性核種で標識された薬剤を投与すると、特定の臓器に集積する。そのとき人体の外部に放出されてくるガンマ線を、人体外に配置した検出器で検出する。検出器に入射するガンマ線の方向をコリメータで規制しながらその入射位置を検出し、その位置ごとに放射線事象をカウントすれば、投影データを収集することができる。このような投影データを多方向から得れば、逆投影法などのアルゴリズムを用いて人体内のＲＩの濃度分布像を再構成することができる。
【０００３】
このＳＰＥＣＴ装置としては、ガンマカメラ回転型や、リング型などが知られている。被検体外部でガンマ線を検出する検出器としてガンマカメラを用い、これを被検体の周囲に回転させるようにしたものが、ガンマカメラ回転型のＳＰＥＣＴ装置である。リング型ＳＰＥＣＴ装置では、多数の検出器をリング型に配列し、そのリングの中央部に被検体を置く。
【０００４】
ところで、これらＳＰＥＣＴ装置では、被検体の内部の核種からの放射線を外部において検出するため、その放射線が被検体の内部で吸収されてしまうことの影響を受けることが避けられない。そこで、再構成画像では被検体の中央部の濃度が異常に低いものとなったり、定量的な測定ができず精度が低いなどの問題が生じるので、その吸収の影響を補正する必要が生じる。
【０００５】
この吸収補正方法としては、均一吸収補正法と不均一吸収補正法の２種類に分類される。均一吸収補正法は、均一な吸収体（つまり被検体内の吸収係数が均一）を仮定して補正する方法であり、Ｓｏｒｅｎｓｏｎ法( Sorenson J.A. : Methods for quantitative measurement of radioactivity in whole body counting. Instrumentation in Nuclear Medicine, eds. Hine G.J. and Sorenson J.A., Vol.2, pp.311-348, Academic Press, New York, 1974 )、Ｃｈａｎｇ法( Chang L.T. : A method for attenuation correction in radionuclide computed tomography. IEEE Trans. Nucl. Sci. NS-25 : 635-643, 1978 )、ＷＢＰ( weighted back projection )法、ＲＰＣ( radial post-correction method )法( Tanaka E., Toyama H. and Murayama H. : Convolutional image reconstruction for quantitative single photon emission computed tomography. Phys. Med. Biol. 29 : 1489-1500, 1984 )などが知られている。これらによれば、頭部ＳＰＥＣＴの場合には比較的良好な結果が得られるためそのような場合に使用されている。
【０００６】
一方、不均一吸収補正法は、トランスミッションＣＴによって実際の被検体ごとに求めた吸収係数分布データを用いて厳密に吸収補正する方法である。
【０００７】
【発明が解決しようとする課題】
しかしながら、従来の吸収補正方法はいずれも問題を持っている。均一補正法は、頭部ＳＰＥＣＴなどでは比較的良好な結果が得られるが、心筋ＳＰＥＣＴのように周囲に肺、筋肉、骨などが存在して不均一な吸収体となっている場合、補正精度が悪く、再構成画像上にアーティファクト（偽像）を生じるため、不適である。
【０００８】
また、不均一吸収補正法は補正精度は高いが、トランスミッションＣＴデータ収集を別個に行なうため、特別なハードウェアが必要であり、さらに収集時間が長くかかったり、トランスミッションＣＴ用の外部線源による放射線照射により患者の被曝も増え、また処理時間が長くかかるなどの欠点がある。
【０００９】
この発明は、上記に鑑み、周囲に肺、筋肉、骨などが存在して不均一な吸収体となっている場合でも比較的高い精度で吸収補正を行なうことができ、かつ再構成画像上のアーティファクトも小さく抑えることができる、新規なＳＰＥＣＴの吸収補正方法を提供することを目的とする。
【００１０】
【課題を解決するための手段】
上記の目的を達成するため、この発明によるＳＰＥＣＴの吸収補正方法においては、被検体内に分布している核種から放出される複数のエネルギーピークの放射線を被検体外部において検出し、検出位置ごとのカウントデータを各エネルギーピークごとに収集し、各検出位置ごとのカウントデータについて、各エネルギーに関するカウントデータの比に、被検体の高エネルギー側ガンマ線に関する吸収係数をｍＨ、低エネルギー側特性Ｘ線に関する吸収係数をｍＬとしたときに
ａ＝ｍＬ／（ｍＬ−ｍＨ）
で求められる定数の因子ａを、べき乗演算することによって各検出位置ごとの吸収補正係数を求め、該吸収補正係数を、対応する検出位置の一方のエネルギーに関するカウントデータに乗算して該カウントデータの吸収補正を行なうことが特徴となっている。
【００１１】
各検出位置ごとのカウントデータの各エネルギーについての比に対する定数の因子のべき乗を考えると、その対応する検出位置の一方のエネルギーに関するカウントデータの補正係数を導くことができる。そして、上記の定数の因子というのは、各エネルギーについての吸収係数の差に対する一方のエネルギーについての吸収係数の比として考えられる。つまり、この定数の因子ａは、被検体の高エネルギー側ガンマ線に関する吸収係数をｍＨ、低エネルギー側特性Ｘ線に関する吸収係数をｍＬとしたときに
ａ＝ｍＬ／（ｍＬ−ｍＨ）
で求められる。そのため、この因子ａは、水、筋肉、肺などであまり異ならない値となる。そこで、同じ値の定数の因子を用いて吸収補正係数を求め、これにより吸収補正しても、水、筋肉、肺などの吸収係数の異なる吸収体が周囲にある場合でも良好に吸収補正することが可能となる。
【００１２】
【発明の実施の形態】
つぎに、この発明の実施の形態について図面を参照しながら詳細に説明する。図１はこの発明をガンマカメラ回転型ＳＰＥＣＴ装置に適用した実施形態を示すものである。この図において、被検体１０の周囲にガンマカメラ２０が矢印のように回転させられるようになっている（回転機構は図示しない）。このガンマカメラ２０は、被検体１０内のＲＩの濃度分布を検出面（シンチレータ２１の面）に投影した、２次元の投影像を撮影するものである。ガンマカメラ２０の回転角度は回転角度検出器２７によって検出される。
【００１３】
ガンマカメラ２０は、平板状のシンチレータ２１と、その背面にライトガイド２２を介して光学的に結合された多数のフォトマルチプライア（ＰＭＴ）２３とを備える。シンチレータ２１の前面には放射線の入射方向を規制するためのコリメータ２５が取り付けられる。たとえば、コリメータ２５によってシンチレータ２１に直角な方向の放射線のみが入射するようにさせられているとする。
【００１４】
入射した放射線は、シンチレータ２１において発光し、その光がライトガイド２２を経てＰＭＴ２３に入射する。ＰＭＴ２３の出力は位置およびエネルギー演算回路２４に送られる。各ＰＭＴ２３の出力はそれらへの入射光量に対応しているため、発光位置に近いＰＭＴ２３ほど大きな出力となる。そこで、位置およびエネルギー演算回路２４は、この関係を利用してＰＭＴ２３の出力をその位置に応じて重み付けして加算するなどにより、シンチレータ２１の平面方向における発光位置を表わす位置信号Ｘ，Ｙを得る。また、全ＰＭＴ２３の出力を加算すれば入射放射線のエネルギーに対応した波高値を持つパルスが得られるので、位置およびエネルギー演算回路２４は、全ＰＭＴ２３の出力を加算して得たエネルギー信号Ｚを出力する。
【００１５】
ここで、被検体１０にはたとえば²⁰¹Ｔｌを核種として含む放射性薬剤が投与されているものとし、心筋ＳＰＥＣＴや腫瘍ＳＰＥＣＴを行なうものとする。この核種の場合、７０ｋｅＶ（および８０ｋｅＶ）をピークとする特性Ｘ線と、１６７ｋｅＶをピークとするガンマ線とが放出され、そのエネルギースペクトルは図２のようになる。そこで、たとえばエネルギー分析器２６におけるエネルギーウインドウを、これらに合わせて中心７０ｋｅＶ、幅２０％と、中心１６７ｋｅＶ、幅２０％の２つに設定する。そして、エネルギー信号Ｚのパルス波高値がこれらのウインドウのそれぞれに入ると、エネルギー分析器２６は、低い方のウインドウに入ったときはデータ収集メモリ２８に、高い方のウインドウに入ったときはデータ収集メモリ２９に、それぞれ信号を送る。
【００１６】
データ収集メモリ２８、２９では、エネルギー分析器２６からの信号ごとに、位置およびエネルギー演算回路２４からの位置信号Ｘ，Ｙおよび回転角度検出器２７からの回転角度信号で指定されるアドレスに「＋１」を格納する（該アドレスでカウントする）。こうして、２つのデータ収集メモリ２８、２９には、ガンマカメラ２０の各回転角度ごとの２次元投影データが、エネルギーウインドウごとに収集されることになる。
【００１７】
補正回路３０は、この２つのデータ収集メモリ２８、２９に収集されたデータを用いてメモリ２８、２９の一方のデータの吸収補正を行なう。ここでは、低エネルギー側放出確率は高エネルギー側に比較して１０倍ほどであるから、この低エネルギー側のデータを画像再構成に用いることとし、これの吸収補正を行なうことにする。
【００１８】
７０ｋｅＶについての２次元投影データの各画素のデータ（測定カウント値）をＮ_Lとすると、つぎの式（１）により補正後のデータ（放出カウント値）Ｎ_L0を得る。
Ｎ_L0＝Ｎ_L・Ｒ^a （１）
ここで、Ｒは各画素ごとにつぎの式（２）で求められる。
Ｒ＝（Ｎ_H／ｋ）／Ｎ_L （２）
Ｎ_Hは、高エネルギー側つまり１６７ｋｅＶについての２次元投影データの各画素のデータ（測定カウント値）である。ｋは、高エネルギー側（１６７ｋｅＶ）のガンマ線の放出確率と、低エネルギー側（７０ｋｅＶ）の特性Ｘ線の放出確率との比であり、上記の通り０．１ほどである。
【００１９】
Ｒをべき乗する因子ａは、高エネルギー側ガンマ線に関する吸収係数ｍＨと低エネルギー側特性Ｘ線に関する吸収係数ｍＬとから、次式（３）により定数として求められる。
ａ＝ｍＬ／（ｍＬ−ｍＨ）（３）
ｍＬ，ｍＨは人体各組織によって異なり、ａもそれに応じて異なるがあまり変わらない値となっている。これを表に示すと、つぎのようになる。

【００２０】
そこで、ａを２≦ａ≦4.2の適当な定数として入力し、上記（２）式のように放出確率の比ｋで校正した高エネルギー側測定カウント値Ｎ_Hと低エネルギー側測定カウント値Ｎ_Lとの比Ｒにａをべき乗して求めた補正係数Ｒ^a を、上記（１）式のように低エネルギー側測定カウント値Ｎ_Lに乗じて補正後のデータを得る。この補正後のデータは画像再構成装置３１に送られ、再構成しようとするスライス面のデータのみを用いて逆投影などの処理を行なうことによりそのスライス面での核種の分布像が再構成される。
【００２１】
つぎに、この吸収補正方法の原理について説明する。図３に示すように、線源４１における低エネルギー側（７０ｋｅＶ）の特性Ｘ線の放出カウントをＮ_L0、高エネルギー側（１６７ｋｅＶ）のガンマ線の放出カウントをＮ_H0とし、これらが長さ（厚さ）ｘの吸収体４２を通って吸収されたとする。吸収後のカウント値として、低エネルギー側につきＮ_L、高エネルギー側につきＮ_Hがそれぞれ測定されたとすると、つぎの関係が成り立つ。
Ｎ_L＝Ｎ_L0・ｅ^-mL・x （４）
Ｎ_H＝Ｎ_H0・ｅ^-mH・x （５）
そのため、Ｎ_Lを補正してＮ_L0を得るには、
Ｎ_L0＝Ｎ_L・ｅ^mL・x （６）
の演算を行なえばよい。
【００２２】
一方、Ｎ_L0とＮ_H0との比がｋであるということは、つぎの式のように表わされる。
Ｎ_H0＝ｋ・Ｎ_L0 （７）
そこで、（４）、（５）式と（７）式とを用いれば、（２）式はつぎのようになる。
Ｒ＝（Ｎ_H／ｋ）／Ｎ_L＝ｅ^(mL-mH)・x （８）
この（８）式の両辺をａ乗したことを考えてみると、
Ｒ^a＝ｅ^{a・(mL-mH)・x} （９）
であり、ａは（３）式の通りであるから、結局
Ｒ^a＝ｅ^mL・x （１０）
となる。つまり、（６）式における吸収補正係数ｅ^mL・xは、Ｒ^aとなる。
【００２３】
このようにして吸収補正することの効果をコンピュータシュミレーションで確認したところ、非常に良好な結果が得られることが分かった。これは、上記の表に示すように、水や筋肉と比べて肺や骨は吸収係数がかなり異なっているが、因子ａについては肺も水や筋肉とほとんど同じ値になるため、同じ値の因子ａを用いても、水、筋肉、肺について高い精度で吸収補正ができることによる。骨についてはａの値が小さいため水のａ値で吸収補正すると骨の部分がやや過補正となるが、全体として、従来の均一補正法のどれよりも格段に補正精度が向上することが確認できた。
【００２４】
なお、上記では、核種として²⁰¹Ｔｌを用いたが、¹¹¹Ｉｎや⁶⁷Ｇａなどの複数エネルギーピークを有する放射線を放出するものであれば、他の核種を用いることもできる。また、上記ではガンマカメラ回転型ＳＰＥＣＴ装置について説明したが、リング型ＳＰＥＣＴ装置にも同様に適用できることはもちろんである。
【００２５】
【発明の効果】
以上説明したように、この発明のＳＰＥＣＴの吸収補正方法によれば、心筋ＳＰＥＣＴのような不均一な吸収体が周囲に存在する場合でも、比較的高い精度で吸収補正することができ、再構成画像上のアーティファクトも抑えることができて、定量性が改善され、診断精度が向上する。また、不均一吸収補正法のようにトランスミッションＣＴデータを収集する必要がないので、特別なハードウェアを用いないで容易に吸収補正を行なうことができる。また、収集時間が余分にかかることや、患者の放射線被曝の問題、あるいは処理時間がかかる問題にも、無縁である。
【図面の簡単な説明】
【図１】この発明の実施の形態を示すブロック図。
【図２】放射線エネルギースペクトルを示すグラフ。
【図３】原理説明のための模式図。
【符号の説明】
１０被検体
２０ガンマカメラ
２１シンチレータ
２２ライトガイド
２３フォトマルチプライア
２４位置およびエネルギー演算回路
２５コリメータ
２６エネルギー分析器
２７回転角度検出器
２８低エネルギー側データ収集メモリ
２９高エネルギー側データ収集メモリ
３０補正回路
３１画像再構成装置
４１線源
４２吸収体[0001]
BACKGROUND OF THE INVENTION
The present invention relates to a SPECT apparatus (single photon emission CT apparatus) used in the field of nuclear medicine diagnosis or science and engineering, and more particularly to an absorption correction method thereof.
[0002]
[Prior art]
The SPECT apparatus uses a radioisotope (RI) of a single photon-emitting radionuclide, detects gamma rays emitted from the radioisotope (RI), and takes a radionuclide distribution image. For example, when a drug labeled with a single photon-emitting radionuclide is administered to the human body, it accumulates in a specific organ. At that time, the gamma rays emitted outside the human body are detected by a detector arranged outside the human body. Projection data can be collected by detecting the incident position while regulating the direction of gamma rays incident on the detector with a collimator and counting radiation events for each position. If such projection data is obtained from multiple directions, an RI density distribution image in the human body can be reconstructed using an algorithm such as a back projection method.
[0003]
As this SPECT apparatus, a gamma camera rotation type, a ring type, and the like are known. A gamma camera rotation type SPECT apparatus uses a gamma camera as a detector for detecting gamma rays outside the object and rotates it around the object. In the ring type SPECT apparatus, a large number of detectors are arranged in a ring shape, and a subject is placed at the center of the ring.
[0004]
By the way, in these SPECT apparatuses, since the radiation from the nuclide inside the subject is detected outside, it is inevitable that the radiation is influenced by being absorbed inside the subject. Therefore, in the reconstructed image, problems such as abnormally low concentration in the center of the subject or low accuracy due to inability to quantitatively measure occur, and it is necessary to correct the influence of the absorption.
[0005]
This absorption correction method is classified into two types, a uniform absorption correction method and a non-uniform absorption correction method. The uniform absorption correction method is a correction method that assumes a uniform absorber (that is, the absorption coefficient in the subject is uniform). The Sorenson method (Sorenson JA: Methods for quantitative measurement of radioactivity in whole body counting. Instrumentation in Nuclear Medicine, eds.Hine GJ and Sorenson JA, Vol.2, pp.311-348, Academic Press, New York, 1974), Chang method (Chang LT: A method for attenuation correction in radionuclide computed tomography.IEEE Trans. Nucl NS-25: 635-643, 1978), WBP (weighted back projection) method, RPC (radial post-correction method) method (Tanaka E., Toyama H. and Murayama H.: Convolutional image reconstruction for quantitative single Photon emission computed tomography. Phys. Med. Biol. 29: 1489-1500, 1984). According to these, in the case of the head SPECT, since a relatively good result is obtained, it is used in such a case.
[0006]
On the other hand, the non-uniform absorption correction method is a method of strictly correcting absorption using the absorption coefficient distribution data obtained for each actual subject by transmission CT.
[0007]
[Problems to be solved by the invention]
However, all the conventional absorption correction methods have problems. In the uniform correction method, relatively good results can be obtained with the head SPECT. However, the correction accuracy is improved when the lungs, muscles, bones, etc. are present in the surrounding area as in the myocardial SPECT, resulting in a non-uniform absorber. Is unsuitable because it produces artifacts on the reconstructed image.
[0008]
In addition, the non-uniform absorption correction method has high correction accuracy, but transmission CT data collection is performed separately, so special hardware is required, and it takes a long time to collect radiation from an external source for transmission CT. Irradiation increases the patient's exposure and takes a long time.
[0009]
In view of the above, the present invention makes it possible to perform absorption correction with relatively high accuracy even when lungs, muscles, bones, and the like are present in the surrounding area, and is a non-uniform absorber. It is an object of the present invention to provide a novel SPECT absorption correction method capable of suppressing artifacts.
[0010]
[Means for Solving the Problems]
In order to achieve the above object, in the SPECT absorption correction method according to the present invention, radiation of a plurality of energy peaks emitted from nuclides distributed in a subject is detected outside the subject, and is detected at each detection position. Count data is collected for each energy peak, the count data for each detection position is the ratio of the count data for each energy, the absorption coefficient for the high energy side gamma rays of the subject is mH, and the absorption for the low energy side characteristic X-rays. When the coefficient is mL, a = mL / (mL-mH)
An absorption correction coefficient for each detection position is obtained by performing a power operation on the constant factor a obtained in step (a), and the absorption correction coefficient is multiplied by count data relating to one energy of the corresponding detection position to It is characterized by performing absorption correction.
[0011]
Considering the power of a constant factor with respect to the ratio of the count data for each detection position for each energy, a correction coefficient for the count data regarding one energy of the corresponding detection position can be derived. The constant factor is considered as the ratio of the absorption coefficient for one energy to the difference in absorption coefficient for each energy. That is, the constant factor a is a = mL / (mL-mH) where mH is the absorption coefficient for the high energy side gamma ray of the subject and mL is the absorption coefficient for the low energy side characteristic X-ray.
Is required. Therefore, the factor a is a value that does not differ so much between water, muscle, lungs, and the like. Therefore, the absorption correction coefficient is calculated using constant factors of the same value, and even if absorption correction is performed by this, it is possible to correct absorption well even when there are absorbers with different absorption coefficients such as water, muscle, lung, etc. Is possible.
[0012]
DETAILED DESCRIPTION OF THE INVENTION
Next, embodiments of the present invention will be described in detail with reference to the drawings. FIG. 1 shows an embodiment in which the present invention is applied to a gamma camera rotating SPECT apparatus. In this figure, the gamma camera 20 is rotated around the subject 10 as indicated by an arrow (the rotation mechanism is not shown). The gamma camera 20 captures a two-dimensional projection image obtained by projecting the RI concentration distribution in the subject 10 onto the detection surface (the surface of the scintillator 21). The rotation angle of the gamma camera 20 is detected by a rotation angle detector 27.
[0013]
The gamma camera 20 includes a flat scintillator 21 and a large number of photomultipliers (PMTs) 23 optically coupled to the back surface of the gamma camera 20 via light guides 22. A collimator 25 for restricting the incident direction of radiation is attached to the front surface of the scintillator 21. For example, it is assumed that only radiation in a direction perpendicular to the scintillator 21 is incident by the collimator 25.
[0014]
The incident radiation is emitted in the scintillator 21, and the light enters the PMT 23 through the light guide 22. The output of the PMT 23 is sent to the position and energy calculation circuit 24. Since the output of each PMT 23 corresponds to the amount of light incident thereon, the PMT 23 closer to the light emission position has a larger output. Accordingly, the position and energy calculation circuit 24 obtains position signals X and Y representing the light emission position in the plane direction of the scintillator 21 by using the relationship and weighting and adding the output of the PMT 23 according to the position. . Further, if the outputs of all the PMTs 23 are added, a pulse having a peak value corresponding to the energy of the incident radiation can be obtained. Therefore, the position and energy calculation circuit 24 outputs an energy signal Z obtained by adding the outputs of all the PMTs 23. To do.
[0015]
Here, for example, it is assumed that a radiopharmaceutical containing ²⁰¹ Tl as a nuclide is administered to the subject 10, and myocardial SPECT or tumor SPECT is performed. In the case of this nuclide, a characteristic X-ray peaking at 70 keV (and 80 keV) and a gamma ray peaking at 167 keV are emitted, and the energy spectrum is as shown in FIG. Therefore, for example, the energy window in the energy analyzer 26 is set to two of the center 70 keV and the width 20%, the center 167 keV and the width 20% in accordance with these. When the pulse peak value of the energy signal Z enters each of these windows, the energy analyzer 26 stores data in the data acquisition memory 28 when entering the lower window and data when entering the higher window. Each signal is sent to the collection memory 29.
[0016]
In the

data collection memories

28 and 29, for each signal from the energy analyzer 26, “+1” is added to the address specified by the position signals X and Y from the position and energy calculation circuit 24 and the rotation angle signal from the rotation angle detector 27. Is stored (counts at the address). Thus, the two

data collection memories

28 and 29 collect two-dimensional projection data for each rotation angle of the gamma camera 20 for each energy window.
[0017]
The correction circuit 30 performs absorption correction of one data in the

memories

28 and 29 using the data collected in the two

data collection memories

28 and 29. Here, since the low energy side emission probability is about 10 times that of the high energy side, the data on the low energy side is used for image reconstruction, and the absorption correction is performed.
[0018]
When the data (measurement count value) of each pixel of the two-dimensional projection data for 70 keV is N _L , corrected data (release count value) N _L0 is obtained by the following equation (1).
N _L0 = N _L · R ^a (1)
Here, R is obtained by the following equation (2) for each pixel.
R = (N _H / k) / N _L (2)
N _H is data (measurement count value) of each pixel of the two-dimensional projection data on the high energy side, that is, 167 keV. k is a ratio between the emission probability of gamma rays on the high energy side (167 keV) and the emission probability of characteristic X rays on the low energy side (70 keV), and is about 0.1 as described above.
[0019]
The factor a that raises R to the right is obtained as a constant from the absorption coefficient mH related to the high energy side gamma ray and the absorption coefficient mL related to the low energy side characteristic X-ray by the following equation (3).
a = mL / (mL-mH) (3)
mL and mH are different depending on each human body tissue, and a is different depending on the tissue but is not so different. This is shown in the table below.

[0020]
Therefore, a is input as an appropriate constant of 2 ≦ a ≦ 4.2, and the high energy side measurement count value _NH and the low energy side measurement count value N _L calibrated with the ratio k of the emission probability as in the above equation (2). correction coefficient R ^a determined by powers of a the ratio R between the obtain data corrected by multiplying the low energy side measurement count value N _L as described above (1). The corrected data is sent to the image reconstruction device 31, and the nuclide distribution image on the slice plane is reconstructed by performing processing such as back projection using only the slice plane data to be reconstructed. The
[0021]
Next, the principle of this absorption correction method will be described. As shown in FIG. 3, the low energy side (70 keV) characteristic X-ray emission count of the radiation source 41 is N _L0 , and the high energy side (167 keV) gamma ray emission count is N _H0 , which are length (thickness). Suppose x is absorbed through the absorber 42 of x. Assuming that N _L is measured on the low energy side and N _H is measured on the high energy side as the count values after absorption, the following relationship is established.
N _L = N _L0 · e ^{-mL · x} (4)
N _H = N _H0 · e ^{-mH · x} (5)
Therefore, in order to obtain N _L0 by correcting the N _L is
N _L0 = N _L · e ^{mL · x} (6)
It is sufficient to perform the operation.
[0022]
On the other hand, the fact that the ratio of N _L0 to N _H0 is k is expressed as the following equation.
N _H0 = k · N _L0 (7)
Therefore, using the equations (4), (5) and (7), the equation (2) becomes as follows.
R = (N _H / k) / N _L = e ^{(mL-mH) · x} (8)
Considering that both sides of this equation (8) are raised to the power a,
R ^a = e ^{a · (mL-mH) · x} (9)
And a is as shown in the equation (3), so that R ^a = e ^{mL · x} (10)
It becomes. That is, the absorption correction coefficient e ^{mL · x} in the equation (6), the R ^a.
[0023]
Thus, when the effect of carrying out the absorption correction was confirmed by computer simulation, it was found that very good results were obtained. As shown in the table above, the absorption coefficient of lungs and bones is considerably different from that of water and muscle. However, the factor a is almost the same as that of water and muscles. Even if the factor a is used, the absorption correction can be performed with high accuracy for water, muscle, and lung. For bone, because the value of a is small, the bone part is slightly overcorrected when it is absorbed and corrected with the a value of water, but overall, it is confirmed that the correction accuracy is significantly improved over any of the conventional uniform correction methods. did it.
[0024]
In the above, ²⁰¹ Tl is used as the nuclide, but other nuclides may be used as long as they emit radiation having multiple energy peaks such as ¹¹¹ In and ⁶⁷ Ga. In the above description, the gamma camera rotation type SPECT apparatus has been described. However, the present invention can also be applied to a ring type SPECT apparatus.
[0025]
【The invention's effect】
As described above, according to the SPECT absorption correction method of the present invention, even when a non-uniform absorber such as myocardial SPECT is present in the vicinity, absorption correction can be performed with relatively high accuracy, and reconstruction Artifacts on the image can also be suppressed, quantitativeness is improved, and diagnostic accuracy is improved. Further, since there is no need to collect transmission CT data unlike the non-uniform absorption correction method, absorption correction can be easily performed without using special hardware. In addition, there is no need for extra collection time, patient radiation exposure, or processing time.
[Brief description of the drawings]
FIG. 1 is a block diagram showing an embodiment of the present invention.
FIG. 2 is a graph showing a radiation energy spectrum.
FIG. 3 is a schematic diagram for explaining the principle.
[Explanation of symbols]
DESCRIPTION OF SYMBOLS 10 Subject 20 Gamma camera 21 Scintillator 22 Light guide 23 Photomultiplier 24 Position and energy calculation circuit 25 Collimator 26 Energy analyzer 27 Rotation angle detector 28 Low energy side data collection memory 29 High energy side data collection memory 30 Correction circuit 31 Image reconstruction device 41 line source 42 absorber

Claims

Radiation of multiple energy peaks emitted from nuclides distributed within the subject is detected outside the subject, and count data for each detection position is collected for each energy peak. Count data for each detection position When the absorption coefficient for the high energy side gamma ray of the subject is mH and the absorption coefficient for the low energy side characteristic X-ray is mL, a = mL / (mL-mH)
An absorption correction coefficient for each detection position is obtained by performing a power operation on the constant factor a obtained in step (a), and the absorption correction coefficient is multiplied by count data relating to one energy of the corresponding detection position to An absorption correction method of SPECT, wherein absorption correction is performed.