JPH08304309A - 物体の放射線量分布図を該物体による減衰に対する補正の後に製作する方法 - Google Patents
物体の放射線量分布図を該物体による減衰に対する補正の後に製作する方法Info
- Publication number
- JPH08304309A JPH08304309A JP8111724A JP11172496A JPH08304309A JP H08304309 A JPH08304309 A JP H08304309A JP 8111724 A JP8111724 A JP 8111724A JP 11172496 A JP11172496 A JP 11172496A JP H08304309 A JPH08304309 A JP H08304309A
- Authority
- JP
- Japan
- Prior art keywords
- radiation
- attenuation
- measurement value
- radiation dose
- propagation
- 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.)
- Pending
Links
- 230000005855 radiation Effects 0.000 title claims abstract description 153
- 238000012937 correction Methods 0.000 title claims abstract description 40
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 9
- 238000005259 measurement Methods 0.000 claims abstract description 74
- 238000000034 method Methods 0.000 claims abstract description 63
- 238000010586 diagram Methods 0.000 claims abstract description 21
- 238000001514 detection method Methods 0.000 claims abstract description 12
- 230000000644 propagated effect Effects 0.000 claims abstract description 3
- 238000013016 damping Methods 0.000 claims description 6
- 238000003325 tomography Methods 0.000 description 11
- 230000002238 attenuated effect Effects 0.000 description 9
- 230000002457 bidirectional effect Effects 0.000 description 5
- 238000002591 computed tomography Methods 0.000 description 5
- 238000002059 diagnostic imaging Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 238000003384 imaging method Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 3
- 230000002285 radioactive effect Effects 0.000 description 3
- 210000001519 tissue Anatomy 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 208000019622 heart disease Diseases 0.000 description 2
- 239000002901 radioactive waste Substances 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 235000007849 Lepidium sativum Nutrition 0.000 description 1
- 244000211187 Lepidium sativum Species 0.000 description 1
- 238000007476 Maximum Likelihood Methods 0.000 description 1
- GKLVYJBZJHMRIY-OUBTZVSYSA-N Technetium-99 Chemical compound [99Tc] GKLVYJBZJHMRIY-OUBTZVSYSA-N 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 210000001367 artery Anatomy 0.000 description 1
- 230000008033 biological extinction Effects 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 230000000747 cardiac effect Effects 0.000 description 1
- 230000002490 cerebral effect Effects 0.000 description 1
- 238000004980 dosimetry Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 230000005489 elastic deformation Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 210000004072 lung Anatomy 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000009659 non-destructive testing Methods 0.000 description 1
- 238000009206 nuclear medicine Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000005258 radioactive decay Effects 0.000 description 1
- 239000000700 radioactive tracer Substances 0.000 description 1
- 229910052704 radon Inorganic materials 0.000 description 1
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical group [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 description 1
- 238000002603 single-photon emission computed tomography Methods 0.000 description 1
- 210000004872 soft tissue Anatomy 0.000 description 1
- 238000010561 standard procedure Methods 0.000 description 1
- BKVIYDNLLOSFOA-OIOBTWANSA-N thallium-201 Chemical compound [201Tl] BKVIYDNLLOSFOA-OIOBTWANSA-N 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/16—Measuring radiation intensity
- G01T1/161—Applications in the field of nuclear medicine, e.g. in vivo counting
- G01T1/1615—Applications in the field of nuclear medicine, e.g. in vivo counting using both transmission and emission sources simultaneously
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/42—Arrangements for detecting radiation specially adapted for radiation diagnosis
- A61B6/4208—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector
- A61B6/4258—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector for detecting non x-ray radiation, e.g. gamma radiation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/16—Measuring radiation intensity
- G01T1/161—Applications in the field of nuclear medicine, e.g. in vivo counting
- G01T1/164—Scintigraphy
- G01T1/1641—Static instruments for imaging the distribution of radioactivity in one or two dimensions using one or several scintillating elements; Radio-isotope cameras
- G01T1/1642—Static instruments for imaging the distribution of radioactivity in one or two dimensions using one or several scintillating elements; Radio-isotope cameras using a scintillation crystal and position sensing photodetector arrays, e.g. ANGER cameras
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/29—Measurement performed on radiation beams, e.g. position or section of the beam; Measurement of spatial distribution of radiation
- G01T1/2914—Measurement of spatial distribution of radiation
- G01T1/2985—In depth localisation, e.g. using positron emitters; Tomographic imaging (longitudinal and transverse section imaging; apparatus for radiation diagnosis sequentially in different planes, steroscopic radiation diagnosis)
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/037—Emission tomography
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Biomedical Technology (AREA)
- General Physics & Mathematics (AREA)
- Optics & Photonics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- Biophysics (AREA)
- Radiology & Medical Imaging (AREA)
- Heart & Thoracic Surgery (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Nuclear Medicine (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
(57)【要約】
【課題】 従来の減衰分布図の生成を要しない、減衰補
正された放射線量分布図製作を実施するための方法を提
案すること。 【解決手段】 物体(2)による放射線の放射量の分布
図製作を、物体による放射線の減衰に対して補正して、
実施するための方法であって、物体に対して幾つかの位
置を採り得る放射線伝播源(3)から物体にフォトンを
送出する方法において、放射源の各位置に対して、検出
手段(7)を用いて、放射線源から放射され且つ物体に
より伝播されたフォトンの伝播計測値(N)の取得と、
物体から放射されたフォトンの放射量計測値の投影図の
取得(E1)とを実施し、伝播計測値の取得と放射線計
測値の投影図の取得とが同じ配置で実施され、更に、放
射線源の各位置に対して、放射線量計測値の投影図(E
6)を補正するために、物体による放射線の減衰補正係
数(E4)を決定し、更に、減衰補正された放射量計測
値の投影図に基づき放射線量分布図(E2)を構成す
る。
正された放射線量分布図製作を実施するための方法を提
案すること。 【解決手段】 物体(2)による放射線の放射量の分布
図製作を、物体による放射線の減衰に対して補正して、
実施するための方法であって、物体に対して幾つかの位
置を採り得る放射線伝播源(3)から物体にフォトンを
送出する方法において、放射源の各位置に対して、検出
手段(7)を用いて、放射線源から放射され且つ物体に
より伝播されたフォトンの伝播計測値(N)の取得と、
物体から放射されたフォトンの放射量計測値の投影図の
取得(E1)とを実施し、伝播計測値の取得と放射線計
測値の投影図の取得とが同じ配置で実施され、更に、放
射線源の各位置に対して、放射線量計測値の投影図(E
6)を補正するために、物体による放射線の減衰補正係
数(E4)を決定し、更に、減衰補正された放射量計測
値の投影図に基づき放射線量分布図(E2)を構成す
る。
Description
【0001】
【発明の属する技術分野】本発明は、物体の放射線量分
布図の製作を、該物体による減衰に対して補正して、実
施するための方法に関する。本発明は、対象の放射線量
の真の分布図製作を要する全ての分野、特に様々な心臓
疾患を研究するための医学画像の分野、及び、特に放射
線廃棄物を監視するための非破壊検査の分野で、使用可
能である。
布図の製作を、該物体による減衰に対して補正して、実
施するための方法に関する。本発明は、対象の放射線量
の真の分布図製作を要する全ての分野、特に様々な心臓
疾患を研究するための医学画像の分野、及び、特に放射
線廃棄物を監視するための非破壊検査の分野で、使用可
能である。
【0002】
【従来の技術】対象即ち物体の放射線量分布図製作を実
施するために、三次元(3D)断層撮影法としても知ら
れる三次元(3D)断層撮影装置を使用することが公知
である。この撮影法は、対象の周りの一連の範囲及び任
意の構成に応じて、対象それ自体から放射された放射線
量の計測値、及び外部放射線源から該対象を通って伝播
された放射線量の計測値の、二次元(2D)取得を得る
ことを可能とする。
施するために、三次元(3D)断層撮影法としても知ら
れる三次元(3D)断層撮影装置を使用することが公知
である。この撮影法は、対象の周りの一連の範囲及び任
意の構成に応じて、対象それ自体から放射された放射線
量の計測値、及び外部放射線源から該対象を通って伝播
された放射線量の計測値の、二次元(2D)取得を得る
ことを可能とする。
【0003】放射線量分布図製作は、各点における放射
線の放射能壊変速度の局部値を得ることを特徴とする。
放射期間は、放射素子の化学的性質によって決まる。放
射性廃棄物の場合、放射性崩壊時間は、一般に長い。医
学画像の場合、放射性トレーサを含む分子が、患者の体
内に注入され、血液により検査すべき組織まで搬送され
る。該組織は、それ自体、フォトンを放出する。この場
合、放射性同位元素は、数時間という短い放射期間を有
する。
線の放射能壊変速度の局部値を得ることを特徴とする。
放射期間は、放射素子の化学的性質によって決まる。放
射性廃棄物の場合、放射性崩壊時間は、一般に長い。医
学画像の場合、放射性トレーサを含む分子が、患者の体
内に注入され、血液により検査すべき組織まで搬送され
る。該組織は、それ自体、フォトンを放出する。この場
合、放射性同位元素は、数時間という短い放射期間を有
する。
【0004】図1は、対象の放射線量分布図製作及び減
衰分布図製作に従来用いられてきた三次元断層撮影装置
を示す。この断層撮影装置1は、検査対象2にフォトン
を供給可能な放射線源3と、放射線源3に整列して対象
2の反対側に位置するコリメータ5と、該コリメータ5
に当接して配設された検出手段7と、を備えている。検
出手段は、幾つかの検出器を備えたガンマカメラであっ
てもよい。
衰分布図製作に従来用いられてきた三次元断層撮影装置
を示す。この断層撮影装置1は、検査対象2にフォトン
を供給可能な放射線源3と、放射線源3に整列して対象
2の反対側に位置するコリメータ5と、該コリメータ5
に当接して配設された検出手段7と、を備えている。検
出手段は、幾つかの検出器を備えたガンマカメラであっ
てもよい。
【0005】放射線源3は、円形軌道Tc上を移動し、
対象2に対して幾つかの位置を取り得る。軌道Tcは、
回転軸線Arに対して垂直な平面である、横軸平面内に
位置している。コリメータ/検出手段組立体は、軌道T
c上で放射線源3と同じ回転動作を行う。
対象2に対して幾つかの位置を取り得る。軌道Tcは、
回転軸線Arに対して垂直な平面である、横軸平面内に
位置している。コリメータ/検出手段組立体は、軌道T
c上で放射線源3と同じ回転動作を行う。
【0006】図1は、対象2の中心を通過する回転軸線
を示しているが、対象2に対して中心からずらしてもよ
い。より詳細には、コリメータ5は、一般に、複数の孔
を有する鉛板から構成されており、各孔は、検出手段7
の一点で計測された放射量の取得コーンを形成し、積分
投影支持体即ち積分線を画定する軸線を有する。断層撮
影装置1において、コリメータ5は、計測値の取得構造
を形成する。この構造は、平行構造、円錐形構造、扇形
構造又は他の公知の構造であってよい。
を示しているが、対象2に対して中心からずらしてもよ
い。より詳細には、コリメータ5は、一般に、複数の孔
を有する鉛板から構成されており、各孔は、検出手段7
の一点で計測された放射量の取得コーンを形成し、積分
投影支持体即ち積分線を画定する軸線を有する。断層撮
影装置1において、コリメータ5は、計測値の取得構造
を形成する。この構造は、平行構造、円錐形構造、扇形
構造又は他の公知の構造であってよい。
【0007】伝播源3は、選定された取得構造に応じ
て、異なる形状をとることができる。伝播源3は、円錐
型コリメータの焦点に位置する点源でもよいし、平面源
或いは線源でもよい。図1に示した装置では、取得構造
が平行型であるので、放射線源3は平面源である。かか
る平行構造の場合、放射線源3は、伝播面をシミュレー
トするように平面に沿って変移させた、線源であっても
よい。
て、異なる形状をとることができる。伝播源3は、円錐
型コリメータの焦点に位置する点源でもよいし、平面源
或いは線源でもよい。図1に示した装置では、取得構造
が平行型であるので、放射線源3は平面源である。かか
る平行構造の場合、放射線源3は、伝播面をシミュレー
トするように平面に沿って変移させた、線源であっても
よい。
【0008】放射線源として線源を備えた図1と同様の
断層撮影装置は、タンP.による文献D1の「シングル
フォトン・放射形コンピュータ断層撮影法(SPEC
T)における放射及び伝播同時計測用の走査線源」(1
993年、JNM、第34巻(第10号)、1752乃
至1760頁)に記載されている。
断層撮影装置は、タンP.による文献D1の「シングル
フォトン・放射形コンピュータ断層撮影法(SPEC
T)における放射及び伝播同時計測用の走査線源」(1
993年、JNM、第34巻(第10号)、1752乃
至1760頁)に記載されている。
【0009】このような三次元断層撮影装置により、幾
つかの放射線量分布図製作方法が実現される。最も簡単
で標準的な方法は、対象から放射された放射線の線量計
測値の投影図の取得と、それに続く対応放射線量分布図
の再構成とから成る。これらの分布図の三次元的再構成
は、選定された取得構造に応じて、解析学的方法又は双
方向の代数学的方法により、行うことができる。
つかの放射線量分布図製作方法が実現される。最も簡単
で標準的な方法は、対象から放射された放射線の線量計
測値の投影図の取得と、それに続く対応放射線量分布図
の再構成とから成る。これらの分布図の三次元的再構成
は、選定された取得構造に応じて、解析学的方法又は双
方向の代数学的方法により、行うことができる。
【0010】平行構造又は扇形構造の場合、一般に、
F.ナテラーによる文献D2の「コンピュータ断層撮影
の数学」(1986年、ジョン・ウイリー&サンズ社
刊)に記載されたフィルタリング/背面投影のアルゴリ
ズムを用いるなどした解析学的方法が利用される。
F.ナテラーによる文献D2の「コンピュータ断層撮影
の数学」(1986年、ジョン・ウイリー&サンズ社
刊)に記載されたフィルタリング/背面投影のアルゴリ
ズムを用いるなどした解析学的方法が利用される。
【0011】円錐形構造の場合、L.フェルドカンプ、
K.C.デヴィッド、及びJ.W.クレスによる文献D
3の「実用的円錐ビームアルゴリズム」(1984年、
米国光学学会誌、第1巻(第6号)、612乃至619
頁)に記載されたような、フィルタリング及び背面投影
による直接解析法を用いることができる。
K.C.デヴィッド、及びJ.W.クレスによる文献D
3の「実用的円錐ビームアルゴリズム」(1984年、
米国光学学会誌、第1巻(第6号)、612乃至619
頁)に記載されたような、フィルタリング及び背面投影
による直接解析法を用いることができる。
【0012】かかる構造の場合、B.D.スミスによる
文献D4の「円錐ビーム投影からの画像の再構成:必要
十分条件と再構成法」(医学画像に関する米国電気電子
学会紀要、MI−第4巻(第1号)、14乃至25頁)
並びにP.グランジートの欧州特許出願公開第0292
402号公報に記載されたラドン領域内の再構成による
直接解析法を用いることも可能である。
文献D4の「円錐ビーム投影からの画像の再構成:必要
十分条件と再構成法」(医学画像に関する米国電気電子
学会紀要、MI−第4巻(第1号)、14乃至25頁)
並びにP.グランジートの欧州特許出願公開第0292
402号公報に記載されたラドン領域内の再構成による
直接解析法を用いることも可能である。
【0013】双方向代数法は、選定した構造に関わりな
く、使用することができる。これらは、構成すべき画像
と線形等価装置により計測されたデータとの結合関係、
及び反復による該装置の解像度についての記載に根拠を
置いている。例えば、かかる方法は、GT.ガルバー
グ、G.L.ゼング、B.M.W.ツイ、及びJ.T.
ハジウスによる文献D5の「円錐ビーム構造を用いたシ
ングルフォトン・放射形コンピュータ断層撮影法(SP
ECT)用の双方向再構成アルゴリズム」(1989
年、画像装置及び技術学会国際大会講演集、第1巻、1
69乃至186頁)に記載されている。
く、使用することができる。これらは、構成すべき画像
と線形等価装置により計測されたデータとの結合関係、
及び反復による該装置の解像度についての記載に根拠を
置いている。例えば、かかる方法は、GT.ガルバー
グ、G.L.ゼング、B.M.W.ツイ、及びJ.T.
ハジウスによる文献D5の「円錐ビーム構造を用いたシ
ングルフォトン・放射形コンピュータ断層撮影法(SP
ECT)用の双方向再構成アルゴリズム」(1989
年、画像装置及び技術学会国際大会講演集、第1巻、1
69乃至186頁)に記載されている。
【0014】しかしながら、かかる双方向代数法は、大
量の計算を要し、工程が長く面倒になる。なお、対象が
放射線を放射する場合、放射線は対象により減衰される
ので、得られる分布図は、減衰された放射線量分布図で
ある、ことに留意されたい。
量の計算を要し、工程が長く面倒になる。なお、対象が
放射線を放射する場合、放射線は対象により減衰される
ので、得られる分布図は、減衰された放射線量分布図で
ある、ことに留意されたい。
【0015】上述した断層撮影装置により使用される方
法は、対象による減衰に対して放射線量分布図を補正す
ることを目指している。これらの方法の一つは、補正さ
れた放射線量分布図の三次元再構成作業のために、均一
減衰分布図製作モデルを用いる。かかる方法は、S.ベ
リーニ、M.ピアセンチニ、C.カフォリオ及びF.ロ
ッカによる文献D6「放射線断層撮影法における組織吸
収の補償」(1979年、米国電気電子学会紀要、音
響、音声及び信号処理、第27巻(第3号)、213乃
至218頁)に記載されている。しかしながら、この方
法は、減衰分布図があまり複雑でない場合に、補正分布
図を提供するだけである。例えば、医学画像分野では、
肺、骨及び柔組織を有する胸部の減衰分布図は余りに複
雑で、補正された放射線量分布図の良好な再構成を得る
ことができない。従って、本方法は、心臓医学画像分野
で使用することができない。
法は、対象による減衰に対して放射線量分布図を補正す
ることを目指している。これらの方法の一つは、補正さ
れた放射線量分布図の三次元再構成作業のために、均一
減衰分布図製作モデルを用いる。かかる方法は、S.ベ
リーニ、M.ピアセンチニ、C.カフォリオ及びF.ロ
ッカによる文献D6「放射線断層撮影法における組織吸
収の補償」(1979年、米国電気電子学会紀要、音
響、音声及び信号処理、第27巻(第3号)、213乃
至218頁)に記載されている。しかしながら、この方
法は、減衰分布図があまり複雑でない場合に、補正分布
図を提供するだけである。例えば、医学画像分野では、
肺、骨及び柔組織を有する胸部の減衰分布図は余りに複
雑で、補正された放射線量分布図の良好な再構成を得る
ことができない。従って、本方法は、心臓医学画像分野
で使用することができない。
【0016】他の方法は、補正された放射線量分布図の
三次元再構成のために、真の減衰分布図を用いる。放射
線量計測値の投影図の取得に加えて、当該方法は、三次
元放射線源から対象を通って放射された放射線の伝播計
測値の取得を実施することを含む。次に、これらの取得
は、そこから減衰投影計測値を導出するように、予め処
理される。これらの方法は、更に、減衰分布図及び減衰
された放射線量分布図に基づいて、減衰補正された放射
線量分布図を作製するように、前記減衰分布図を再構成
することを含む。
三次元再構成のために、真の減衰分布図を用いる。放射
線量計測値の投影図の取得に加えて、当該方法は、三次
元放射線源から対象を通って放射された放射線の伝播計
測値の取得を実施することを含む。次に、これらの取得
は、そこから減衰投影計測値を導出するように、予め処
理される。これらの方法は、更に、減衰分布図及び減衰
された放射線量分布図に基づいて、減衰補正された放射
線量分布図を作製するように、前記減衰分布図を再構成
することを含む。
【0017】かくして、これらの方法は、減衰計測値の
真の分布図を使用する。しかしながら、かかる減衰分布
図製作法は、伝播計測用に使用されたエネルギに対して
定義された減衰係数を提供する。従って、該減衰係数
は、放射線量計測に使用される放射性素子、即ち対象
(例えば放射性廃棄物)内に自然に生じたり或いは(特
に医学画像において)対象又は人体内に導入される放射
性素子、に対応するエネルギに適合するように、補正さ
れなければならない。
真の分布図を使用する。しかしながら、かかる減衰分布
図製作法は、伝播計測用に使用されたエネルギに対して
定義された減衰係数を提供する。従って、該減衰係数
は、放射線量計測に使用される放射性素子、即ち対象
(例えば放射性廃棄物)内に自然に生じたり或いは(特
に医学画像において)対象又は人体内に導入される放射
性素子、に対応するエネルギに適合するように、補正さ
れなければならない。
【0018】かかる方法では、放射線量分布図の減衰補
正は、以下のいずれかにより行われる。 a) 減衰投影図に基づいて再構成された放射線量の後
補正による。この場合、放射線量分布図の各ヴォクセル
(基本量)は、ヴォクセルの周りの全取得方向に応じた
平均減衰量に連係した係数により、補正される。補正
は、再構成−背面投影により一度だけ、或いは反復的に
行うことができる。この方法は、大きな減衰系内に配設
された小さい源にも適用される。しかしながら、逆の場
合に使用すれば、対象の中心で過補正が生じることにな
る。かかる補正方法は、L.T.チャンによる文献D7
の「放射性核種コンピュータ断層撮影法(RCT)にお
ける減衰方法」(1978年、米国電気電子学会紀要、
核科学、第25巻(第1号)、198乃至123頁)に
記載されている。
正は、以下のいずれかにより行われる。 a) 減衰投影図に基づいて再構成された放射線量の後
補正による。この場合、放射線量分布図の各ヴォクセル
(基本量)は、ヴォクセルの周りの全取得方向に応じた
平均減衰量に連係した係数により、補正される。補正
は、再構成−背面投影により一度だけ、或いは反復的に
行うことができる。この方法は、大きな減衰系内に配設
された小さい源にも適用される。しかしながら、逆の場
合に使用すれば、対象の中心で過補正が生じることにな
る。かかる補正方法は、L.T.チャンによる文献D7
の「放射性核種コンピュータ断層撮影法(RCT)にお
ける減衰方法」(1978年、米国電気電子学会紀要、
核科学、第25巻(第1号)、198乃至123頁)に
記載されている。
【0019】b) 或いは、同じ投影ビーム上で計算さ
れた平均減衰量による放射線量の投影図の前補正によ
る。この方法は、反復的に使用することもでき、対応す
る減衰系と同じ形状を有して均質な放射線源にも適用さ
れる。この方法は、マゼット他による文献D8の「SP
ECTにおける非均一減衰分布用の反復的再構成法」
(1993年、JNM、第34巻(第7号)、1204
乃至1209頁に記載されている。
れた平均減衰量による放射線量の投影図の前補正によ
る。この方法は、反復的に使用することもでき、対応す
る減衰系と同じ形状を有して均質な放射線源にも適用さ
れる。この方法は、マゼット他による文献D8の「SP
ECTにおける非均一減衰分布用の反復的再構成法」
(1993年、JNM、第34巻(第7号)、1204
乃至1209頁に記載されている。
【0020】c) 或いは、放射線量投影図の減衰補正
係数の反復的計算による。収束時、検出手段の画素用の
この係数は、放射量に基づいて計算された背面投影図と
放射線量及び減衰量に基づいて計算された減衰背面投影
図との比である。かかる方法は、T.モロズミ、M.ナ
カジマ、K.オガワ、及びS.ユタによる文献D9の
「シングルフォトン・放射形コンピュータ断層撮影法
(SPECT)のための、減衰分布の情報を用いた減衰
補正法」(1984年、医学画像技術、第2号、20乃
至28頁)に記載されている。
係数の反復的計算による。収束時、検出手段の画素用の
この係数は、放射量に基づいて計算された背面投影図と
放射線量及び減衰量に基づいて計算された減衰背面投影
図との比である。かかる方法は、T.モロズミ、M.ナ
カジマ、K.オガワ、及びS.ユタによる文献D9の
「シングルフォトン・放射形コンピュータ断層撮影法
(SPECT)のための、減衰分布の情報を用いた減衰
補正法」(1984年、医学画像技術、第2号、20乃
至28頁)に記載されている。
【0021】d) 或いは、評価及び背面投影段階を考
慮に入れて減衰を補正する最大尤度反復再構成法によ
る。計算時間で特に有利なこの方法は、ガルバーグによ
る文献D5に記載されている。
慮に入れて減衰を補正する最大尤度反復再構成法によ
る。計算時間で特に有利なこの方法は、ガルバーグによ
る文献D5に記載されている。
【0022】図2は、補正された放射線量分布図製作を
実施するための二種類の公知の方法の機能図である。こ
れらの方法の一方は、破線形で示した経路に対応する。
この方法は、減衰投影図E4を用いて構成された真の
(即ち計測された)減衰分布図E5から供給された情報
に基づいて、補正放射量投影図E6内で、減衰された放
射量投影図E1を補正することを含む。この方法は、マ
ゼットによる文献D8に記載されたものと実質的に同一
である。
実施するための二種類の公知の方法の機能図である。こ
れらの方法の一方は、破線形で示した経路に対応する。
この方法は、減衰投影図E4を用いて構成された真の
(即ち計測された)減衰分布図E5から供給された情報
に基づいて、補正放射量投影図E6内で、減衰された放
射量投影図E1を補正することを含む。この方法は、マ
ゼットによる文献D8に記載されたものと実質的に同一
である。
【0023】図2に示した他方の方法は、実線で示した
経路に対応する。この方法は、一方で、減衰された放射
線量分布図E3を構成し、他方で、減衰分布図E5を構
成し、これらの二つの分布図から、補正された放射線量
分布図E2を得ることを含む。この方法は、チャンによ
る文献D7に記載されたものと実質的に同一である。こ
れらの二つの方法を含め、減衰分布図に基づいて放射線
量分布図を補正する全て方法は、減衰分布図を得ること
が難しいため、実施が困難である。
経路に対応する。この方法は、一方で、減衰された放射
線量分布図E3を構成し、他方で、減衰分布図E5を構
成し、これらの二つの分布図から、補正された放射線量
分布図E2を得ることを含む。この方法は、チャンによ
る文献D7に記載されたものと実質的に同一である。こ
れらの二つの方法を含め、減衰分布図に基づいて放射線
量分布図を補正する全て方法は、減衰分布図を得ること
が難しいため、実施が困難である。
【0024】かくして、減衰分布図は、放射用と同じガ
ンマカメラを用いて得られる伝播計測値を取得して計算
することができる。この特定の場合では、ガンマカメラ
の検出器は、放射線の計測値を取得するために十分な寸
法を有するが、伝播計測値の取得を省略することも多
い。この省略のため、再構成された減衰分布図は、補正
の質に影響する瑕疵を含むことになる。
ンマカメラを用いて得られる伝播計測値を取得して計算
することができる。この特定の場合では、ガンマカメラ
の検出器は、放射線の計測値を取得するために十分な寸
法を有するが、伝播計測値の取得を省略することも多
い。この省略のため、再構成された減衰分布図は、補正
の質に影響する瑕疵を含むことになる。
【0025】減衰分布図は、伝播計測値の取得を省略し
ないX線スキャナに基づいて計算することもできる。し
かしながら、この場合、放射線量計測値の取得と伝播計
測値の取得との間の対象の(医学画像の場合は患者の)
不可避の弾性変形により、幾何配置の基準位置における
伝播及び放射データ構成の再設定に困難が伴う。更に、
減衰分布図に基づいたこれらの方法は、一般に双方向の
減衰補正による減衰分布図の再構成と放射線量分布図の
再構成とを要するので、計算量が増加する。
ないX線スキャナに基づいて計算することもできる。し
かしながら、この場合、放射線量計測値の取得と伝播計
測値の取得との間の対象の(医学画像の場合は患者の)
不可避の弾性変形により、幾何配置の基準位置における
伝播及び放射データ構成の再設定に困難が伴う。更に、
減衰分布図に基づいたこれらの方法は、一般に双方向の
減衰補正による減衰分布図の再構成と放射線量分布図の
再構成とを要するので、計算量が増加する。
【0026】
【発明が解決しようとする課題】本発明は、上述した方
法の欠点を除去することを目的とする。このため、本発
明は、従来の減衰分布図の生成を要しない、減衰補正さ
れた放射線量分布図製作を実施するための方法を提案す
る。
法の欠点を除去することを目的とする。このため、本発
明は、従来の減衰分布図の生成を要しない、減衰補正さ
れた放射線量分布図製作を実施するための方法を提案す
る。
【0027】
【課題を解決するための手段】本発明は、物体による放
射線の放射量の分布図製作を、前記物体による放射線の
減衰に対して補正して、実施するための方法であって、
前記物体に対して幾つかの位置を採り得る放射線伝播源
から前記物体にフォトンを送出する方法において、前記
放射源の各位置に対して、検出手段を用いて、前記放射
線源から放射され且つ前記物体により伝播されたフォト
ンの伝播計測値の取得と、前記物体から放射されたフォ
トンの放射量計測値の投影図の取得とを実施し、前記伝
播計測値の取得と前記放射線計測値の投影図の取得とが
同じ配置で実施され、更に、放射線源の各位置に対し
て、放射線量計測値の投影図を補正するために、物体に
よる放射線の減衰補正係数を決定し、更に、減衰補正さ
れた放射量計測値の投影図に基づき放射線量分布図を構
成する、ことを特徴とする方法に関する。
射線の放射量の分布図製作を、前記物体による放射線の
減衰に対して補正して、実施するための方法であって、
前記物体に対して幾つかの位置を採り得る放射線伝播源
から前記物体にフォトンを送出する方法において、前記
放射源の各位置に対して、検出手段を用いて、前記放射
線源から放射され且つ前記物体により伝播されたフォト
ンの伝播計測値の取得と、前記物体から放射されたフォ
トンの放射量計測値の投影図の取得とを実施し、前記伝
播計測値の取得と前記放射線計測値の投影図の取得とが
同じ配置で実施され、更に、放射線源の各位置に対し
て、放射線量計測値の投影図を補正するために、物体に
よる放射線の減衰補正係数を決定し、更に、減衰補正さ
れた放射量計測値の投影図に基づき放射線量分布図を構
成する、ことを特徴とする方法に関する。
【0028】本発明によれば、減衰量の決定は、物体が
無い場合に得られる基準伝播計測値Noの取得を実行
し、Emを計測された放射線量計測値の投影図、Ecを
減衰補正された放射線量計測値の投影図としてEc=E
m×Cが成立するように補正係数Cを取得値N及びNo
の関数として決定することを含む。
無い場合に得られる基準伝播計測値Noの取得を実行
し、Emを計測された放射線量計測値の投影図、Ecを
減衰補正された放射線量計測値の投影図としてEc=E
m×Cが成立するように補正係数Cを取得値N及びNo
の関数として決定することを含む。
【0029】好ましくは、前記伝播計測値の取得と前記
放射線量計測値の投影図の取得とは、同時に行われ、前
記放射線伝播源から放射されるフォトンと前記物体から
放射されるフォトンとは、異なるエネルギを有する。
放射線量計測値の投影図の取得とは、同時に行われ、前
記放射線伝播源から放射されるフォトンと前記物体から
放射されるフォトンとは、異なるエネルギを有する。
【0030】物体が同心の減衰領域と放射領域とを有す
る本発明の一実施形態において、補正係数CはNを伝播
計測値の取得値、Noを基準伝播計測値の取得値とし
て、
る本発明の一実施形態において、補正係数CはNを伝播
計測値の取得値、Noを基準伝播計測値の取得値とし
て、
【0031】
【数3】
【0032】により近似される。
【0033】物体が一定幅の減衰領域とその減衰領域に
対して中心がずれた放射領域とを有する本発明の別の実
施形態において、補正係数Cは、Lcを前記検出手段に
最も近い減衰領域の縁部と放射領域の中心との間の距
離、Nを伝播計測値の取得値、Noを基準伝播計測値の
取得値として、
対して中心がずれた放射領域とを有する本発明の別の実
施形態において、補正係数Cは、Lcを前記検出手段に
最も近い減衰領域の縁部と放射領域の中心との間の距
離、Nを伝播計測値の取得値、Noを基準伝播計測値の
取得値として、
【0034】
【数4】
【0035】により近似される。
【0036】
【発明の実施形態】本発明に係る方法は、対象からの放
射線に対して減衰補正された、放射線量分布図製作法を
実施すること即ち対象の分布図製作から成る。方法を使
用する場所に応じて、対象は、射出性物体でも、或い
は、放射性素子を導入した患者の身体でもよく、自然に
又は非自然に放射線を放射し得る他の型の物体でもよ
い。本明細書の以下の部分では、任意に対象即ち物体と
呼ぶ。
射線に対して減衰補正された、放射線量分布図製作法を
実施すること即ち対象の分布図製作から成る。方法を使
用する場所に応じて、対象は、射出性物体でも、或い
は、放射性素子を導入した患者の身体でもよく、自然に
又は非自然に放射線を放射し得る他の型の物体でもよ
い。本明細書の以下の部分では、任意に対象即ち物体と
呼ぶ。
【0037】図3は、本方法の主要な段階を示す機能図
である。本発明の方法は、先ず、減衰された放射線量計
測値の投影図の取得(段階E1)及び平行伝播計測値の
取得を実施することから成る、ことが分かる。E4で得
られる伝播計測値は、真の放射線量計測値の投影図を得
る(段階E6)ように、E1の減衰された放射線量計測
値の投影図を補正するために使用される。次に、段階E
6の真の、即ち補正された放射線量計測値の投影図に基
づいて公知の方法により、補正された放射線量分布図が
製作される(段階E2)。
である。本発明の方法は、先ず、減衰された放射線量計
測値の投影図の取得(段階E1)及び平行伝播計測値の
取得を実施することから成る、ことが分かる。E4で得
られる伝播計測値は、真の放射線量計測値の投影図を得
る(段階E6)ように、E1の減衰された放射線量計測
値の投影図を補正するために使用される。次に、段階E
6の真の、即ち補正された放射線量計測値の投影図に基
づいて公知の方法により、補正された放射線量分布図が
製作される(段階E2)。
【0038】本発明の方法によれば、放射線源3は、物
体の周りの幾つかの位置に順次配設される。コリメータ
/検出手段組立体は、前記放射線源3の位置に関わりな
く、前記放射線源3の直径方向反対側に位置する。本方
法は、放射線源3の各位置で、物体から放射された放射
線の放射線量計測値の投影図の取得、及び物体2を介し
て放射線源3から放射された放射線の伝播計測値の取得
を実施する工程を含む。これら二種類の取得は、同一構
成の断層撮影装置に対しても実施される。
体の周りの幾つかの位置に順次配設される。コリメータ
/検出手段組立体は、前記放射線源3の位置に関わりな
く、前記放射線源3の直径方向反対側に位置する。本方
法は、放射線源3の各位置で、物体から放射された放射
線の放射線量計測値の投影図の取得、及び物体2を介し
て放射線源3から放射された放射線の伝播計測値の取得
を実施する工程を含む。これら二種類の取得は、同一構
成の断層撮影装置に対しても実施される。
【0039】放射線量の計測値の取得及び伝播の計測値
の取得は、順次実施し得る、ことは明らかである。しか
しながら、二つの取得型の間で物体が変移する危険を避
けるために、これらは同時に実施することが好ましい。
このような実施形態では、物体により放射されるフォト
ンと比較して、異なるエネルギフォトンを放射する源
を、放射線源3として選定することができる。次に、放
射線源3より放射されるフォトンから物体2より放射さ
れるフォトンを分離し得るように、検出手段7を選定す
る。この分離は、エネルギ差により、或いは放射線源3
が移動する場合は検出手段に対する関与領域の摺動によ
り、或いは他の公知の手段により、行うことができる。
の取得は、順次実施し得る、ことは明らかである。しか
しながら、二つの取得型の間で物体が変移する危険を避
けるために、これらは同時に実施することが好ましい。
このような実施形態では、物体により放射されるフォト
ンと比較して、異なるエネルギフォトンを放射する源
を、放射線源3として選定することができる。次に、放
射線源3より放射されるフォトンから物体2より放射さ
れるフォトンを分離し得るように、検出手段7を選定す
る。この分離は、エネルギ差により、或いは放射線源3
が移動する場合は検出手段に対する関与領域の摺動によ
り、或いは他の公知の手段により、行うことができる。
【0040】例えば、前記フォトン分離は、二つのエネ
ルギ窓、即ち、放射フォトン(例えば、物体2から放射
された)のエネルギ上に中心決めされた一方の窓と、伝
送フォトン(例えば、放射線源3から放射された)のエ
ネルギ上に中心決めされた他方の窓とを用いて、行うこ
とができる。
ルギ窓、即ち、放射フォトン(例えば、物体2から放射
された)のエネルギ上に中心決めされた一方の窓と、伝
送フォトン(例えば、放射線源3から放射された)のエ
ネルギ上に中心決めされた他方の窓とを用いて、行うこ
とができる。
【0041】伝送計測値の取得として、物体による減衰
の生じた後に検出手段7上で検出されたフォトン数の計
測となることは言うまでもない。この減衰は、物体を構
成する材料及び使用した放射線源に係る物理的パラメー
タである。減衰計測値の投影図即ちより簡単に言えば減
衰量の投影図は、選定線に沿った物体による減衰量を合
計したものである。
の生じた後に検出手段7上で検出されたフォトン数の計
測となることは言うまでもない。この減衰は、物体を構
成する材料及び使用した放射線源に係る物理的パラメー
タである。減衰計測値の投影図即ちより簡単に言えば減
衰量の投影図は、選定線に沿った物体による減衰量を合
計したものである。
【0042】この減衰量は、物体による二つの伝播計測
値の取得(即ち簡単に言えば、伝播の取得)により決定
される。 − 第一の伝播の取得は、物体が移動するときに実施さ
れる基準取得であり、実際には、放射線源3により放射
される放射線量の計測値の取得である。この取得値をN
oで表す。 − 第二の伝播の取得は、上述したような通常の伝送状
態で実施されるものであり、この第二の計測値の取得値
をNで表す。 次に、以下の式に基づいて、減衰量μを決定する。
値の取得(即ち簡単に言えば、伝播の取得)により決定
される。 − 第一の伝播の取得は、物体が移動するときに実施さ
れる基準取得であり、実際には、放射線源3により放射
される放射線量の計測値の取得である。この取得値をN
oで表す。 − 第二の伝播の取得は、上述したような通常の伝送状
態で実施されるものであり、この第二の計測値の取得値
をNで表す。 次に、以下の式に基づいて、減衰量μを決定する。
【0043】
【数5】
【0044】ここで、Noは、基準伝播の取得値であ
り、Nは、計測された伝播の取得値であり、Xμは、選
定線に沿った減衰量μの積分値である。
り、Nは、計測された伝播の取得値であり、Xμは、選
定線に沿った減衰量μの積分値である。
【0045】この減衰量μに基づいて、減衰量μの放射
線量計測値の投影図の補正を可能とする補正係数Cを決
定することができる。かくして、補正された放射線量計
測値の投影図Ecは、以下の等式により、計測された放
射線量計測値の投影図Emと結合される。 Ec=Em × C (1) ここで、 Ec=Xf (2) Em=Xμf Xfは、ある線に沿った放射線量計測値の積分値であ
り、Xufは、ある線に沿った減衰放射線量計測値の積
分値である。
線量計測値の投影図の補正を可能とする補正係数Cを決
定することができる。かくして、補正された放射線量計
測値の投影図Ecは、以下の等式により、計測された放
射線量計測値の投影図Emと結合される。 Ec=Em × C (1) ここで、 Ec=Xf (2) Em=Xμf Xfは、ある線に沿った放射線量計測値の積分値であ
り、Xufは、ある線に沿った減衰放射線量計測値の積
分値である。
【0046】かくして、この補正係数Cは、研究対象で
ある物体の形状及び位置に従属する。Cの等価物は、放
射線量ファントム又はゴーストと呼ばれる物体の一般形
状の関数として、決定された。本明細書の以下におい
て、この係数Cを、研究すべき対象即ち物体に対応する
ファントム型の関数として、どのように決定するかにつ
いて説明する。
ある物体の形状及び位置に従属する。Cの等価物は、放
射線量ファントム又はゴーストと呼ばれる物体の一般形
状の関数として、決定された。本明細書の以下におい
て、この係数Cを、研究すべき対象即ち物体に対応する
ファントム型の関数として、どのように決定するかにつ
いて説明する。
【0047】最初に考察するファントムを図4の(A)
に示す。ファントム9は、直径Laを有する第一の円筒
即ち減衰円筒8と、該減衰円筒8と同じ軸線Arを有し
て直径Laより小さい直径Leを有する第二の円筒即ち
放射円筒10と、により構成されている。かかるファン
トム9の場合、円筒軸線Arに垂直な投影ビームRに対
し、その中心で補正係数Cを分析的に計算することがで
きる。係数Cは、次の式で表すことができる。
に示す。ファントム9は、直径Laを有する第一の円筒
即ち減衰円筒8と、該減衰円筒8と同じ軸線Arを有し
て直径Laより小さい直径Leを有する第二の円筒即ち
放射円筒10と、により構成されている。かかるファン
トム9の場合、円筒軸線Arに垂直な投影ビームRに対
し、その中心で補正係数Cを分析的に計算することがで
きる。係数Cは、次の式で表すことができる。
【0048】
【数6】
【0049】ここで、この第一のファントムの特定の場
合において、∫E・dlはXfに対応し、
合において、∫E・dlはXfに対応し、
【0050】
【数7】
【0051】はXufに対応する。
【0052】指数関数
【0053】
【数8】
【0054】の差を有限展開して等式(3)で置換する
ことにより、以下の次数の補正係数Cが得られる。
ことにより、以下の次数の補正係数Cが得られる。
【0055】
【数9】
【0056】核医学画像診断に応用する場合、大脳取得
構造を研究するために、直径Leを直径Laから極僅か
に異ならせたこの種のファントムを、使用することがで
きる。このように直径Leが直径Laより極めて小さい
ファントムを使用することにより、患者の動脈と比較し
て放射線量が小さい心臓を用いて、心臓取得構造を研究
することが可能となる。
構造を研究するために、直径Leを直径Laから極僅か
に異ならせたこの種のファントムを、使用することがで
きる。このように直径Leが直径Laより極めて小さい
ファントムを使用することにより、患者の動脈と比較し
て放射線量が小さい心臓を用いて、心臓取得構造を研究
することが可能となる。
【0057】後者の場合補正係数Cの式は前述した場合
と同様に近似することができる。
と同様に近似することができる。
【0058】
【数10】
【0059】図4の(B)は、ファントムの別の例を示
す。前記ファントム9(図4の(A))と同様、図4の
(B)のファントム11は、直径Laの減衰円筒8と、
Laより小さい直径Leを有する放射円筒10と、から
構成されている。このファントムの場合も、投影ビーム
Rは、円筒の中心を通過する。しかしながら、この投影
ビームRは、円筒の軸線Arに対して垂直ではない。こ
の場合、更にLe<<Laであるとき、補正係数Cは、
やはり、
す。前記ファントム9(図4の(A))と同様、図4の
(B)のファントム11は、直径Laの減衰円筒8と、
Laより小さい直径Leを有する放射円筒10と、から
構成されている。このファントムの場合も、投影ビーム
Rは、円筒の中心を通過する。しかしながら、この投影
ビームRは、円筒の軸線Arに対して垂直ではない。こ
の場合、更にLe<<Laであるとき、補正係数Cは、
やはり、
【0060】
【数11】
【0061】として表すことができる。
【0062】同様に、ファントム9と同じ形状を有する
が投影ビームRが円筒軸線Arを通過しないファントム
の場合も、補正係数Cを等式(6)と同様に表すことが
できる。
が投影ビームRが円筒軸線Arを通過しないファントム
の場合も、補正係数Cを等式(6)と同様に表すことが
できる。
【0063】従って、以上に定義したような補正係数C
は、検出手段7の任意の点で適用可能である。かくし
て、放射線量計測値の投影図Emは、以下のように、補
正することができる。
は、検出手段7の任意の点で適用可能である。かくし
て、放射線量計測値の投影図Emは、以下のように、補
正することができる。
【0064】
【数12】
【0065】又は
【0066】
【数13】
【0067】ここで、Ecは、減衰補正された放射線量
計測値の投影図であり、Paは、Pa=μ・Laとし
て、検出手段の同一画素に対する減衰量の投影図であ
る。
計測値の投影図であり、Paは、Pa=μ・Laとし
て、検出手段の同一画素に対する減衰量の投影図であ
る。
【0068】図4の(C)は、減衰領域12と該減衰領
域12内に位置決めされた放射領域14とを有する、フ
ァントム13を示す。放射領域14は、減衰領域12に
対して中心がずれている。減衰領域12は長さLaを有
し、放射領域14は長さLeを有する。検出手段7に近
い方の減衰領域12の縁部と放射領域14との間の、投
影ビームRに沿った長さはLcである。
域12内に位置決めされた放射領域14とを有する、フ
ァントム13を示す。放射領域14は、減衰領域12に
対して中心がずれている。減衰領域12は長さLaを有
し、放射領域14は長さLeを有する。検出手段7に近
い方の減衰領域12の縁部と放射領域14との間の、投
影ビームRに沿った長さはLcである。
【0069】この場合、LeがLaよりずっと小さいと
すると、補正係数Cは、
すると、補正係数Cは、
【0070】
【数14】
【0071】として表すことができる。
【0072】補正された放射線量計測値の投影図Ec
は、放射線量計測値の投影図Emの関数として、以下の
ように表される。 Ec = Em(ePa)Lc/La (7) 又は Ec = Em(No/N)Lc/La ここでPa=μ・Laである。
は、放射線量計測値の投影図Emの関数として、以下の
ように表される。 Ec = Em(ePa)Lc/La (7) 又は Ec = Em(No/N)Lc/La ここでPa=μ・Laである。
【0073】かくして、放射線量の投影図は、減衰分布
図の再構成を要することなく、伝播計測値に基づいて直
接、減衰補正される。
図の再構成を要することなく、伝播計測値に基づいて直
接、減衰補正される。
【0074】この方法は、従来技術におけるように、非
省略式の放射線量投影図を取得し得る検出手段7を備え
た装置によっても、実施可能である。しかしながら、こ
の方法によれば、補正に有用な伝播の取得情報は直接導
出され、省略情報は補正に有用ではないので、補正量に
影響を及ぼすことなく、伝播の取得を省略することがで
きる。
省略式の放射線量投影図を取得し得る検出手段7を備え
た装置によっても、実施可能である。しかしながら、こ
の方法によれば、補正に有用な伝播の取得情報は直接導
出され、省略情報は補正に有用ではないので、補正量に
影響を及ぼすことなく、伝播の取得を省略することがで
きる。
【0075】本発明の一実施形態によれば、対象又は患
者の物理的輪郭の計測を可能とする装置、即ち自動輪郭
描写装置を備えたガンマカメラを使用している。この輪
郭を用いて、ガンマカメラは、対象又は患者に接触する
ことなくできるだけ接近して検出手段7を回転させる。
者の物理的輪郭の計測を可能とする装置、即ち自動輪郭
描写装置を備えたガンマカメラを使用している。この輪
郭を用いて、ガンマカメラは、対象又は患者に接触する
ことなくできるだけ接近して検出手段7を回転させる。
【0076】例えば、心臓画像診断では、かかる自動輪
郭描写装置を有するソファ・メジカル(SOPHA-MEDICAL
、登録商標)ガンマカメラを使用することができる。
この例では、テクネチウム99の均質平面放射線源をガ
ンマカメラと整列させて、且つ物体の反対側に位置決め
している。放射線源と物体との距離は、例えば、60c
mとすることができる。
郭描写装置を有するソファ・メジカル(SOPHA-MEDICAL
、登録商標)ガンマカメラを使用することができる。
この例では、テクネチウム99の均質平面放射線源をガ
ンマカメラと整列させて、且つ物体の反対側に位置決め
している。放射線源と物体との距離は、例えば、60c
mとすることができる。
【0077】この装置は、種々の心臓疾患を診断するた
めに用いることができる。使用するプロトコルは、タリ
ウム201で標識化した薬品放射性素子の注入としても
よい。断層撮影装置は、トレーサの分布の三次元分布図
を製作する。この場合、エネルギ補正係数の値は、K=
1.210である。
めに用いることができる。使用するプロトコルは、タリ
ウム201で標識化した薬品放射性素子の注入としても
よい。断層撮影装置は、トレーサの分布の三次元分布図
を製作する。この場合、エネルギ補正係数の値は、K=
1.210である。
【0078】同様に、減衰分布図は(既に説明したよう
に)エネルギ補正されなければならないので、前記補正
係数Cを増倍率Kにより補正する必要がある。増倍率K
は、放射線量計測に用いるエネルギに対する水の減衰係
数の比、及び伝播計測に用いるエネルギに対する水の減
衰係数の比であり、前記水の減衰係数は、既知の物理量
である。
に)エネルギ補正されなければならないので、前記補正
係数Cを増倍率Kにより補正する必要がある。増倍率K
は、放射線量計測に用いるエネルギに対する水の減衰係
数の比、及び伝播計測に用いるエネルギに対する水の減
衰係数の比であり、前記水の減衰係数は、既知の物理量
である。
【0079】
【発明の効果】以上述べたように、本発明によれば、減
衰に起因する瑕疵及び変形に対して補正がなされた心臓
の放射線量分布図を調べることができるので、より良い
診断が可能となる。更に、補正は、迅速に行われ、検出
器が小さく省略式の伝播のみを検出し得る場合でも、可
能である。
衰に起因する瑕疵及び変形に対して補正がなされた心臓
の放射線量分布図を調べることができるので、より良い
診断が可能となる。更に、補正は、迅速に行われ、検出
器が小さく省略式の伝播のみを検出し得る場合でも、可
能である。
【図1】減衰補正された又は非補正の放射線量分布図製
作の実施を可能とする三次元断層撮影装置を示す概略斜
視図。
作の実施を可能とする三次元断層撮影装置を示す概略斜
視図。
【図2】減衰補正された放射線量分布図製作の実施する
ための二つの公知の方法を示した機能図。
ための二つの公知の方法を示した機能図。
【図3】減衰分布図を用いることなく、減衰補正された
放射線量分布図製作法を実施するための本発明の方法を
示した機能図。
放射線量分布図製作法を実施するための本発明の方法を
示した機能図。
【図4】(A)、(B)、及び(C)は、それぞれ、本
発明の方法で用いるファントム即ちゴーストの例を示し
た図。
発明の方法で用いるファントム即ちゴーストの例を示し
た図。
Claims (5)
- 【請求項1】 物体(2)による放射線の放射量の分布
図を該物体による減衰に対する補正の後に製作する方法
であって、 前記物体に対して幾つかの位置を採り得る放射線伝播源
(3)から前記物体にフォトンを送出し、 前記放射源の各位置に対して、検出手段(7)を用い
て、前記放射線源から放射され且つ前記物体により伝播
されたフォトンの伝播計測値Nの取得と、前記物体から
放射されたフォトンの放射量計測値の投影図の取得(E
1)とを実施し、 前記伝播計測値の取得と前記放射線計測値の投影図の取
得とが同じ配置で実施され、 放射線源の各位置に対して、放射線量計測値の投影図
(E6)を補正するために、物体による放射線の減衰補
正係数(E4)を決定し、 減衰補正された放射量計測値の投影図に基づき放射線量
分布図(E2)を構成する、 ことを特徴とする方法。 - 【請求項2】 減衰量の決定において、物体が無い場合
に得られる基準伝播計測値Noの取得を実行し、Emを
計測された放射線量計測値の投影図、Ecを減衰補正さ
れた放射線量計測値の投影図としてEc=Em×Cが成
立するように補正係数Cを取得値N及びNoの関数とし
て決定することを含む、減衰補正係数の決定動作から成
る、ことを特徴とする請求項1記載の方法。 - 【請求項3】 前記伝播計測値の取得と前記放射線量計
測値の投影図の取得とが、同時に行われ、前記放射線伝
播源から放射されるフォトンと前記物体から放射される
フォトンとが、異なるエネルギを有することを特徴とす
る請求項1記載の方法。 - 【請求項4】 前記物体が、互いの中心が一致する減衰
領域(8)と放射領域(10)とによりモデル化され、
補正係数Cが、Nを伝播計測値の取得値、Noを基準伝
播計測値の取得値として 【数1】 により近似されることを特徴とする請求項2記載の方
法。 - 【請求項5】 放射線量の計測の一方の方向において、
前記物体が、幅(La)の減衰領域(12)と該減衰領
域に対して中心がずれた放射領域(14)とによりモデ
ル化され、 補正係数Cが、前記方向で、Lcを前記検出手段に最も
近い減衰領域(12)の縁部と放射領域(14)の中心
との間の距離、Nを伝播計測値の取得値、Noを基準伝
播計測値の取得値として、 【数2】 により近似されることを特徴とする請求項2記載の方
法。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FR9505342A FR2733841B1 (fr) | 1995-05-04 | 1995-05-04 | Procede de realisation de la cartographie d'emission d'un corps corrigee de l'attenuation par ce corps |
| FR95-05342 | 1995-05-04 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPH08304309A true JPH08304309A (ja) | 1996-11-22 |
Family
ID=9478702
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP8111724A Pending JPH08304309A (ja) | 1995-05-04 | 1996-05-02 | 物体の放射線量分布図を該物体による減衰に対する補正の後に製作する方法 |
Country Status (5)
| Country | Link |
|---|---|
| US (1) | US5814817A (ja) |
| EP (1) | EP0741308A1 (ja) |
| JP (1) | JPH08304309A (ja) |
| FR (1) | FR2733841B1 (ja) |
| IL (1) | IL118104A0 (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005193018A (ja) * | 2003-12-16 | 2005-07-21 | Ge Medical Systems Global Technology Co Llc | 画像処理のシステム及び方法 |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5682036A (en) * | 1995-12-29 | 1997-10-28 | Adac Laboratories | Method and apparatus for accurately calibrating an attenuation map for emission computed tomography |
| JP4170449B2 (ja) * | 1998-07-07 | 2008-10-22 | 株式会社東芝 | トランスミッションctのトランケーション補正装置、核医学診断装置及びトランケーション補正方法 |
| US6310968B1 (en) * | 1998-11-24 | 2001-10-30 | Picker International, Inc. | Source-assisted attenuation correction for emission computed tomography |
| GB0230110D0 (en) * | 2002-12-24 | 2003-01-29 | British Nuclear Fuels Plc | Improvements in and relating to apparatus and methods for materials investigation |
| WO2012010959A2 (en) * | 2010-07-20 | 2012-01-26 | Szabolcs Osvath | Improved single photon emission computed tomography system and method |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4633398A (en) * | 1983-12-07 | 1986-12-30 | General Electric Company | Attenuation compensated emission reconstruction with simultaneous attenuation factor evaluation |
| JPS61204582A (ja) * | 1985-03-08 | 1986-09-10 | Hitachi Ltd | 放射能分布測定方法及び装置 |
| FR2650397B1 (fr) * | 1989-07-28 | 1991-10-04 | Informatek Ste Nle | Dispositif a scintillation utilisable pour mesurer l'attenuation par tomographie de transmission |
| US5155365A (en) * | 1990-07-09 | 1992-10-13 | Cann Christopher E | Emission-transmission imaging system using single energy and dual energy transmission and radionuclide emission data |
| US5338936A (en) * | 1991-06-10 | 1994-08-16 | Thomas E. Kocovsky, Jr. | Simultaneous transmission and emission converging tomography |
| US5210421A (en) * | 1991-06-10 | 1993-05-11 | Picker International, Inc. | Simultaneous transmission and emission converging tomography |
| US5390225A (en) * | 1992-06-29 | 1995-02-14 | Siemens Medical Systems, Inc. | Mapping slices of the human body into regions having a constant linear attenuation coefficient for correcting images acquired during a nuclear medicine study for attenuation artifacts |
| US5461232A (en) * | 1993-01-08 | 1995-10-24 | Mccandless; Brian K. | Pulse transmission scintigraphic imaging |
-
1995
- 1995-05-04 FR FR9505342A patent/FR2733841B1/fr not_active Expired - Fee Related
-
1996
- 1996-05-01 IL IL11810496A patent/IL118104A0/xx unknown
- 1996-05-02 EP EP96400934A patent/EP0741308A1/fr not_active Withdrawn
- 1996-05-02 JP JP8111724A patent/JPH08304309A/ja active Pending
- 1996-05-03 US US08/642,432 patent/US5814817A/en not_active Expired - Fee Related
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005193018A (ja) * | 2003-12-16 | 2005-07-21 | Ge Medical Systems Global Technology Co Llc | 画像処理のシステム及び方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| US5814817A (en) | 1998-09-29 |
| EP0741308A1 (fr) | 1996-11-06 |
| FR2733841A1 (fr) | 1996-11-08 |
| FR2733841B1 (fr) | 1997-05-30 |
| IL118104A0 (en) | 1996-09-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US5376795A (en) | Emission-transmission imaging system using single energy and dual energy transmission and radionuclide emission data | |
| EP1006370B1 (en) | Image reconstruction in tomography | |
| US8472683B2 (en) | Motion correction in tomographic images | |
| RU2606561C2 (ru) | Компенсация усечения для итерационной реконструкции в компьютерной томографии (кт) с коническим пучком в комбинированных системах офэкт/кт | |
| US8098916B2 (en) | System and method for image-based attenuation correction of PET/SPECT images | |
| US5818050A (en) | Collimator-free photon tomography | |
| JP3128634B2 (ja) | 同時透過・放出型集束断層撮影法 | |
| US5338936A (en) | Simultaneous transmission and emission converging tomography | |
| CN101849247B (zh) | 用于生成衰减图的设备和方法 | |
| CN101528131B (zh) | 对带有运动伪影的图像进行的伪影校正 | |
| Ter-Pogossian | Basic principles of computed axial tomography | |
| US7737406B2 (en) | Compensating for truncated CT images for use as attenuation maps in emission tomography | |
| US20080135769A1 (en) | Attenuation correction of pet image using image data acquired with an mri system | |
| US5461232A (en) | Pulse transmission scintigraphic imaging | |
| CN104093360A (zh) | 核成像系统 | |
| EP0648341B1 (en) | Spect artifact correction for compensating attenuation coefficient variations through reconstruction of a scatter image | |
| CN101080651B (zh) | 成像系统及成像方法 | |
| US6539103B1 (en) | Method and apparatus for image reconstruction using a knowledge set | |
| US6921902B2 (en) | Scatter correction device for radiative tomographic scanner | |
| JP6021347B2 (ja) | 医用画像撮像装置及び医用画像撮像方法 | |
| US6410919B1 (en) | Positron gun for pet attenuation measurements | |
| US20040210132A1 (en) | Simulation of nuclear medical imaging | |
| US7262416B2 (en) | Gamma camera with dynamic threshold | |
| JPH08304309A (ja) | 物体の放射線量分布図を該物体による減衰に対する補正の後に製作する方法 | |
| US7323689B2 (en) | Attenuation correction in nuclear medicine studies by simultaneous transmission and emission data measurement and estimation of emission-to-transmission crosstalk |