Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JPH06105295B2 - Scattered ray correction method of ECT device - Google Patents
[go: Go Back, main page]

JPH06105295B2 - Scattered ray correction method of ECT device - Google Patents

Scattered ray correction method of ECT device

Info

Publication number
JPH06105295B2
JPH06105295B2 JP61180624A JP18062486A JPH06105295B2 JP H06105295 B2 JPH06105295 B2 JP H06105295B2 JP 61180624 A JP61180624 A JP 61180624A JP 18062486 A JP18062486 A JP 18062486A JP H06105295 B2 JPH06105295 B2 JP H06105295B2
Authority
JP
Japan
Prior art keywords
data
correction
scattered
radiation
filter function
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Lifetime
Application number
JP61180624A
Other languages
Japanese (ja)
Other versions
JPS6311886A (en
Inventor
昌治 天野
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shimadzu Corp
Original Assignee
Shimadzu Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shimadzu Corp filed Critical Shimadzu Corp
Publication of JPS6311886A publication Critical patent/JPS6311886A/en
Publication of JPH06105295B2 publication Critical patent/JPH06105295B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Landscapes

  • Apparatus For Radiation Diagnosis (AREA)
  • Nuclear Medicine (AREA)

Description

【発明の詳細な説明】 産業上の利用分野 この発明は、放射性同位体(RI)で標識された化合物を
患者に投与して、その2次元的な分布を断層像として得
るECT装置(エミッション型コンピュータ断層撮影装
置)に関し、とくにその散乱線によって誤差を生じたデ
ータを補正する方法に関する。
TECHNICAL FIELD The present invention relates to an ECT device (emission type) in which a compound labeled with a radioisotope (RI) is administered to a patient and a two-dimensional distribution thereof is obtained as a tomographic image. The present invention relates to a computer tomography apparatus), and more particularly to a method for correcting data having an error due to its scattered rays.

従来の技術 患者に投与された標識化合物から放出される放射線は、
患者身体内を通過して外部に出るため、身体内で吸収さ
れて減衰したり、散乱されて誤った情報を提供したりす
ることにもなる。これらを補正するための方法が種々模
索されているが、なかでも散乱線の補正法は確立されて
いず、いろいろと試みられている。
Radiation emitted from labeled compounds administered to patients is
As it passes through the patient's body to the outside, it may be absorbed and attenuated in the body, or scattered to provide erroneous information. Various methods for correcting these have been sought, but among them, a method for correcting scattered rays has not been established, and various methods have been tried.

たとえば、PET(ポジトロンECT)の場合、第10図Aに示
すようにリング状検出器配列1の中に、RIを含む部分7
と含まない部分8とを有するファントム6を配置したと
きの中心を通るX軸上のプロファイルデータをとってみ
ると第10図Bのようになり、本来RIが存在しない部分に
も散乱線のために誤って得られたデータが存在している
ことが分る。そこで、ファントム6外の部分9のデータ
が散乱線による誤ったデータのみであるとして、この部
分9の平均値を求め、この平均値をプロファイルデータ
の全体から差し引き演算して散乱線の補正をしようとい
うのである。
For example, in the case of PET (positron ECT), a portion 7 containing RI in the ring-shaped detector array 1 as shown in FIG. 10A.
Taking the profile data on the X-axis passing through the center when the phantom 6 having the part 8 and the part 8 not containing is arranged, it becomes as shown in FIG. It can be seen that there is data obtained by mistake. Therefore, assuming that the data of the portion 9 outside the phantom 6 is only erroneous data due to scattered rays, obtain the average value of this portion 9 and subtract this average value from the entire profile data to correct the scattered rays. That is why.

発明が解決しようとする問題点 しかしながら、上記の従来の一例として示した補正法で
は、単に一定値を一律に差し引くというだけであるから
精度が悪いという問題がある。
Problems to be Solved by the Invention However, the correction method shown as the above conventional example has a problem that the accuracy is poor because it simply subtracts a constant value uniformly.

他の従来の種々の補正法についても、補正用データの収
集に時間や手間がかかったり、特別な装置を必要として
コストが嵩んだり、あるいは補正精度が悪いとか、補正
用データを用いて実際の被写体のデータを補正するデー
タ処理に時間がかかる等の問題がある。
Regarding other various conventional correction methods, it takes time and effort to collect the correction data, the cost increases due to the need for a special device, or the correction accuracy is poor. There is a problem that the data processing for correcting the subject data takes time.

この発明は、補正用データ収集に時間や手間がかから
ず、そのために特別な装置を必要とすることもなく安価
で、補正精度も高くて、補正用データを用いて実際の被
写体のデータを補正するデータ処理を短時間で行なうこ
とができる、ECT装置の散乱線補正法を提供することを
目的とする。
The present invention does not require any time or labor to collect the correction data, does not require a special device for that, is inexpensive, has a high correction accuracy, and uses the correction data to obtain the actual subject data. An object of the present invention is to provide a scattered radiation correction method for an ECT device, which can perform data processing for correction in a short time.

問題点を解決するための手段 この発明によるECT装置の散乱線補正法では、断層像を
得ようとする面内で1点を占めるように線源を配置し、
この線源のデータを、周囲に散乱体がある状態とない状
態とについて収集し、これらのデータの間の関係から散
乱線補正用フィルタ関数を求め、実際の被写体について
得られたデータに対してこのフィルタ関数を作用させて
散乱線の補正を行なうことが特徴となっている。
Means for Solving the Problems In the scattered radiation correction method of the ECT apparatus according to the present invention, the radiation source is arranged so as to occupy one point in the plane to obtain a tomographic image,
The data of this radiation source was collected for the condition with and without scatterers in the surroundings, and the filter function for scattered radiation correction was obtained from the relationship between these data, and for the data obtained for the actual subject. The feature is that the scattered radiation is corrected by applying this filter function.

作用 点線源のデータを、周囲に散乱体がある状態とない状態
とについて収集すると、散乱体がある場合は散乱線によ
ってぼけた像が得られるため、プロファイルデータも散
乱線がない場合に比べてなまったものとなる。これら両
データの比をとってみると、それを散乱線によってぼけ
たデータを本来のデータに戻すためのデータとして使え
ることが分る。
Action Collecting point source data with and without scatterers in the surroundings gives a blurred image due to scatterers in the presence of scatterers, so the profile data is also compared to the case without scatterers. It becomes blunt. If you take the ratio of these two data, you can see that it can be used as data to restore the original data to the data blurred by scattered radiation.

そこで、この両データの比をある関数で近似し、この関
数をフィルタ関数として実際の被写体のデータに作用さ
せれば、散乱線の補正をすることができる。
Therefore, the scattered radiation can be corrected by approximating the ratio of these two data by a certain function and applying this function as a filter function to the actual data of the object.

そして、通常、プロファイルデータを逆投影しに画像を
再構成する際に、各プロファイルデータに対してフィル
タリングを行なうので、このとき同時に上記の散乱線補
正用フィルタ関数を作用させることができるため、補正
のために必要な時間が少なくて済む。
In general, when the image is reconstructed by backprojecting the profile data, each profile data is filtered. Therefore, at the same time, the above-mentioned scattered radiation correction filter function can be actuated. Less time is needed for

実施例 PETに適用した一実施例について説明する。まず、第1
図Aのようにリング型検出器配列1の中にライン状線源
2を断層面に垂直に配置する。すると、実質的に、断層
面に点状線源を置いたのと同じになる。そしてこの状態
でデータ収集を行なう。X方向のプロファイルデータは
第1図Bのようになる。
Example An example applied to PET will be described. First, the first
As shown in FIG. A, the linear radiation source 2 is arranged in the ring-type detector array 1 perpendicularly to the tomographic plane. Then, it is substantially the same as placing a point source on the fault plane. Then, data collection is performed in this state. The profile data in the X direction is as shown in FIG. 1B.

つぎに第2図Aのようにこのライン線源2の周囲に散乱
体3を置く。この散乱体3は、人体に似た散乱特性を有
する水やアクリルにより構成される。この状態でデータ
収集を行なうと、線源2より放出された放射線の一部が
散乱体3で散乱するためデータにぼけが加わり、同じX
方向のプロファイルデータは第2図Bのように散乱体3
がないときのデータ(第1図B)とほとんど同じである
がそれに比較して多少なまって広がったものとなる。
Next, as shown in FIG. 2A, the scatterer 3 is placed around the line source 2. The scatterer 3 is made of water or acrylic having a scattering characteristic similar to that of a human body. If data is collected in this state, a part of the radiation emitted from the radiation source 2 is scattered by the scatterer 3, and the data is blurred, and the same X
The profile data of the direction is scatterer 3 as shown in FIG. 2B.
It is almost the same as the data when there is no (Fig. 1B), but it becomes a little wider than that and spreads.

これら第1図Bおよび第2図Bのプロファイルデータを
フーリエ変換したものも、それぞれ第3図および第4図
のようにほとんど同じになり、適当なスケーリングを行
なった後重ね合せればほとんど重なり、低周波部分で散
乱線成分による差異が表われているに過ぎない。そこ
で、これらフーリエ変換後のデータの比をとってみる
と、第5図Aのようなデータが得られる。
The Fourier transforms of the profile data shown in FIGS. 1B and 2B are almost the same as shown in FIGS. 3 and 4, respectively, and if they are appropriately scaled and then overlapped, they almost overlap each other. Only the difference due to the scattered ray component appears in the low frequency part. Therefore, by taking the ratio of the data after the Fourier transform, the data as shown in FIG. 5A is obtained.

ところで、この第5図Aのデータは、散乱線によって誤
差を含むようになったデータを元の誤差のない状態に戻
すための補正データとして考えることができる。そこ
で、第5図Aのデータをある関数で近似すれば、第5図
Bのような関数が得られ、この関数は散乱線補正用フィ
ルタ関数として使用できる。なお、第1図Bで示した散
乱体のない場合のデータは、第2図Bで示した散乱体3
のある場合のデータから曲線近似によって作成すること
もできる。
By the way, the data shown in FIG. 5A can be considered as correction data for returning the data including an error due to scattered radiation to the original error-free state. Therefore, if the data of FIG. 5A is approximated by a certain function, a function as shown in FIG. 5B is obtained, and this function can be used as a scattered radiation correction filter function. The data in the case without the scatterer shown in FIG. 1B is the scatterer 3 shown in FIG. 2B.
It can also be created by curve fitting from the data in the case of.

このようにして散乱線補正用フィルタ関数が得られた
ら、実際の被写体に関して収集されたプロファイルデー
タにこの関数を作用させてフィルタリングすれば、散乱
線の補正ができる。画像再構成は、通常、プロファイル
データに対してフィルタリングした後プロファイルデー
タを逆投影して行なわれるので、そのフィルタリングに
際して同時に上記の散乱線補正用フィルタ関数を作用さ
せればよく、こうすることによって補正のために必要な
データ処理に要する時間を短縮することができる。
When the scattered radiation correction filter function is obtained in this way, the scattered radiation can be corrected by applying the function to the profile data collected for the actual subject to perform filtering. Since image reconstruction is usually performed by back-projecting profile data after filtering the profile data, it is sufficient to apply the above scattered radiation correction filter function at the same time when performing the filtering. Therefore, it is possible to shorten the time required for the data processing required.

なお、この散乱線補正用フィルタ関数は、中央の1つの
位置に点状線源を置いたときに得られるデータから求め
ているが、この散乱線補正用フィルタ関数は他の位置に
関して収集したデータから求めてもあまり変らないもの
と思われる。しかしより正確を期するならば、いくつか
の位置について求めることも望ましい。
The scattered radiation correction filter function is obtained from the data obtained when the point radiation source is placed at one central position. However, this scattered radiation correction filter function is data collected at other positions. It seems that it doesn't change much even if asked from. However, for more accuracy, it is also desirable to determine for some positions.

ところで、上記の散乱線補正方法では、断層面以外から
の散乱線の影響が無視できない場合には精度が悪くな
る。このような場合には、断層面外からの散乱線成分を
見積って上記の補正用フィルタ関数を修正すればよい。
これには、たとえば、従来の散乱線補正方法で用いた第
10図Aのようなファントム6を用いる。前述のように、
このファントム6は、RIを含む部分7と含まない部分8
とを有し、そのX軸上のプロファイルデータは第10図B
のようになるので、本来RIの存在しない部分9のデータ
が散乱線成分のデータということになる。そこで、この
部分9のデータによって断層面外からの散乱線成分を見
積り、上記の散乱線補正用フィルタ関数を修正する。た
とえば、部分9のデータ量に応じて修正することにより
第5図Cのような補正用フィルタ関数を得る。
By the way, in the above scattered ray correction method, the accuracy is deteriorated when the influence of scattered rays from other than the tomographic plane cannot be ignored. In such a case, the above-mentioned correction filter function may be corrected by estimating the scattered radiation component from outside the tomographic plane.
This includes, for example, the
10 Use a phantom 6 as shown in FIG. As aforementioned,
This phantom 6 has a portion 7 including RI and a portion 8 not including RI.
, And the profile data on the X axis is shown in FIG. 10B.
Therefore, the data of the portion 9 where RI does not originally exist is the data of the scattered radiation component. Therefore, the scattered ray component from outside the tomographic plane is estimated from the data of this portion 9 and the above scattered ray correcting filter function is corrected. For example, the correction filter function as shown in FIG. 5C is obtained by correcting the data amount of the portion 9.

上記では、実際の被写体のデータにつき散乱線補正を行
なう場合を説明したが、吸収補正用データを収集する散
乱線の影響が現われるので、吸収補正用データから散乱
線による誤差を除くことにも利用できる。すなわち、吸
収補正用データは、まず、第6図に示すようにリング状
線源4をリング型検出器配列1の中に置いてデータ収集
し、つぎに第7図に示すようにリング状線源4中に実際
の被写体5を入れてデータ収集し、それらの比を求めて
得る。すなわち、第6図の場合のデータには被写体5に
よる吸収の影響が現われず、第7図の場合のデータには
その吸収の影響が現われているので、それらの比を求め
ればそれが吸収補正用データとなるからである。ところ
が、この第7図の場合のデータには、リング状線源4か
ら放出された放射線が被写体5で散乱することの影響が
ある。そこで第8図、第9図のようにライン状線源2
の、被写体5がない場合とある場合とのデータを収集
し、上記と同様に、これらのデータをフーリエ変換しそ
の比から散乱線補正用フィルタ関数を作成する。このフ
ィルタ関数を第7図で得たデータに作用させてフィルタ
リングすれば、散乱線の補正ができる。こうして、散乱
線の影響を除いた吸収補正用データを得ることができ
る。
In the above, the case of performing scattered ray correction on the actual subject data was explained, but since the influence of scattered rays that collects absorption correction data appears, it is also used to remove the error due to scattered rays from the absorption correction data. it can. That is, the absorption correction data is first collected by placing the ring-shaped radiation source 4 in the ring-shaped detector array 1 as shown in FIG. 6, and then collecting the ring-shaped radiation as shown in FIG. The actual subject 5 is placed in the source 4, data is collected, and their ratio is obtained. That is, since the influence of absorption by the subject 5 does not appear in the data in the case of FIG. 6 and the influence of absorption by the subject 5 appears in the data in the case of FIG. This is because it becomes the data for use. However, the data in the case of FIG. 7 is affected by the radiation emitted from the ring-shaped radiation source 4 being scattered by the subject 5. Therefore, as shown in FIG. 8 and FIG.
The data for the case where there is no subject 5 and the case where there is the subject 5 are collected, and similarly to the above, Fourier transform is performed on these data and a scattered radiation correction filter function is created from the ratio. Scattered rays can be corrected by applying this filter function to the data obtained in FIG. 7 for filtering. In this way, absorption correction data that excludes the influence of scattered radiation can be obtained.

発明の効果 この発明にかかるECT装置の散乱線補正法によれば、補
正用データは1度だけ収集しておけばよく、各被写体毎
に、その都度補正用データ収集するという時間や手間が
かかることがない。また、そのために特別な装置を必要
とすることもなく安価で、簡便に補正用データを得るこ
とができる。さらに、単に一定値を一律に差し引くこと
に比較して補正の精度を向上させることができる。ま
た、補正用データを用いて実際の被写体のデータを補正
するデータ処理を短時間で行なうことができる。
EFFECTS OF THE INVENTION According to the scattered radiation correction method of the ECT apparatus according to the present invention, the correction data only needs to be collected once, and it takes time and effort to collect the correction data for each subject each time. Never. Further, for that reason, a special device is not required, and the correction data can be easily obtained at a low cost. Further, the accuracy of correction can be improved as compared with the case where a constant value is simply subtracted. Further, it is possible to perform the data processing for correcting the data of the actual subject using the correction data in a short time.

【図面の簡単な説明】[Brief description of drawings]

第1図乃至第5図はこの発明の一実施例を説明するため
のもので、第1図Aはライン状線源の散乱体のない状態
でのデータを収集することを示す模式図、第1図Bは第
1図Aで得られた1つのプロファイルデータを示すグラ
フ、第2図Aはライン状線源の散乱体のある状態でのデ
ータを収集することを示す模式図、第2図Bは第2図A
で得られた1つのプロファイルデータを示すグラフ、第
3図及び第4図は第1図B及び第2図Bのデータをフー
リエ変換したデータを示すグラフ、第5図Aは第3図お
よび第4図のデータの比を示すグラフ、第5図B、Cは
補正用フィルタ関数をそれぞれ表わすグラフ、第6図、
第7図、第8図および第9図はこの発明を吸収補正用デ
ータの散乱線補正に適用することを説明するための、線
源や被写体の配置をそれぞれ示す模式図、第10図は従来
の散乱線補正法を説明するためのもので、第10図Aはフ
ァントムの配置を示す模式図、第10図Bは第10図Aで得
られた1つのプロファイルデータを示すグラフである。 1……リング型検出器配列、2……ライン状線源 3……散乱体、4……リング状線源 5……被写体、6……ファントム
FIGS. 1 to 5 are for explaining one embodiment of the present invention, and FIG. 1A is a schematic diagram showing that data is collected in a state where there is no scatterer of a linear source, FIG. 1B is a graph showing one profile data obtained in FIG. 1A, FIG. 2A is a schematic view showing collecting data in a state where a scatterer of a linear source is collected, and FIG. B is A in FIG.
A graph showing one profile data obtained in FIG. 3, FIGS. 3 and 4 are graphs showing data obtained by Fourier transforming the data of FIGS. 1B and 2B, and FIG. 5A is FIGS. 4 is a graph showing the ratio of the data in FIG. 4, FIGS. 5B and 5C are graphs showing the correction filter function, FIG. 6,
FIGS. 7, 8 and 9 are schematic diagrams showing the arrangement of a radiation source and a subject, respectively, for explaining the application of the present invention to scattered radiation correction of absorption correction data, and FIG. 10A is a schematic diagram showing the arrangement of phantoms, and FIG. 10B is a graph showing one profile data obtained in FIG. 10A. 1 ... Ring type detector array, 2 ... Line source 3 ... Scatterer, 4 ... Ring source 5 ... Subject, 6 ... Phantom

Claims (2)

【特許請求の範囲】[Claims] 【請求項1】断層像を得ようとする面内で1点を占める
ように線源を配置し、この線源のデータを、周囲に散乱
体がある状態とない状態とについて収集し、これらのデ
ータの間の関係から散乱線補正用フィルタ関数を求め、
実際の被写体について得られたデータに対してこのフィ
ルタ関数を作用させて散乱線の補正を行なうことを特徴
とするECT装置の散乱線補正法。
1. A radiation source is arranged so as to occupy one point in a plane where a tomographic image is to be obtained, and data of this radiation source are collected with and without a scatterer around them. The filter function for scattered radiation correction is obtained from the relationship between the data of
A scattered radiation correction method for an ECT device, characterized in that this filter function is applied to the data obtained for an actual subject to correct scattered radiation.
【請求項2】上記フィルタ関数は、上記断層像を得よう
とする面内に線源を置いたときの、該線源が存在しない
領域で得られるデータにより修正されることを特徴とす
る特許請求の範囲第1項記載のECT装置の散乱線補正
法。
2. The filter function is corrected by data obtained in a region where the radiation source does not exist when the radiation source is placed in a plane to obtain the tomographic image. A scattered radiation correction method for an ECT device according to claim 1.
JP61180624A 1986-03-27 1986-07-31 Scattered ray correction method of ECT device Expired - Lifetime JPH06105295B2 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP6941886 1986-03-27
JP61-69418 1986-03-27

Publications (2)

Publication Number Publication Date
JPS6311886A JPS6311886A (en) 1988-01-19
JPH06105295B2 true JPH06105295B2 (en) 1994-12-21

Family

ID=13402044

Family Applications (1)

Application Number Title Priority Date Filing Date
JP61180624A Expired - Lifetime JPH06105295B2 (en) 1986-03-27 1986-07-31 Scattered ray correction method of ECT device

Country Status (1)

Country Link
JP (1) JPH06105295B2 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003057348A (en) * 2001-08-10 2003-02-26 Shimadzu Corp Positron ECT device

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111257377A (en) * 2020-02-17 2020-06-09 青岛理工大学 ECT sensor calibration method

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003057348A (en) * 2001-08-10 2003-02-26 Shimadzu Corp Positron ECT device

Also Published As

Publication number Publication date
JPS6311886A (en) 1988-01-19

Similar Documents

Publication Publication Date Title
EP1522045B1 (en) Motion artifact correction of tomographical images
US7623691B2 (en) Method for helical windmill artifact reduction with noise restoration for helical multislice CT
Lewitt Processing of incomplete measurement data in computed tomography
EP3067864B1 (en) Iterative reconstruction with enhanced noise control filtering
EP1308898A2 (en) System and method for image reconstruction
US7057178B1 (en) Method and system for imaging using a filter for Time-of-Flight PET
AU2019271915A1 (en) Method and system for motion correction in CT imaging
US8126248B2 (en) Procedures for the presentation of medical images
JPH07508347A (en) How to map slices of the human body to regions with constant linear attenuation coefficients to correct images obtained during nuclear medicine examinations for attenuation artifacts
JPH02125371A (en) Method and apparatus for estimating contour of object with sector beam computer type tomograph
JPH0375045A (en) Ct data compensation for object outside visual field of regenerated image
US7623615B2 (en) X-ray CT image reconstruction method and X-ray CT system
JPH03225292A (en) Artifact/sensitivity correcting device for spect device
JPH06105295B2 (en) Scattered ray correction method of ECT device
US6763081B2 (en) Cone beam computed tomography imaging system and method providing efficient utilization of detector area
JP3381384B2 (en) Correction method for nuclear medicine imaging device
CN117671074B (en) Image reconstruction method, device, electronic equipment and CT imaging system
JPH06105296B2 (en) Scattered ray correction method of ECT device
JP2006133204A (en) SPECT equipment
JP4439022B2 (en) Nuclear medicine diagnostic apparatus and image processing method of nuclear medicine diagnostic apparatus
JPH06308243A (en) Scintillation camera and radiation diagnostic method using it
Lodge et al. A clinical evaluation of the quantitative accuracy of simultaneous emission/transmission scanning in whole-body positron emission tomography
JP4387758B2 (en) SPECT apparatus and SPECT image reconstruction method
JPS6346136A (en) X-ray ct apparatus
Xia et al. Dose reduction achieved by dynamically collimating the redundant rays in fan-beam and cone-beam CT
<