JP4532005B2 - Ｘ線ｃｔ装置及びその画像表示方法 - Google Patents
Ｘ線ｃｔ装置及びその画像表示方法 Download PDFInfo
- Publication number
- JP4532005B2 JP4532005B2 JP2001067150A JP2001067150A JP4532005B2 JP 4532005 B2 JP4532005 B2 JP 4532005B2 JP 2001067150 A JP2001067150 A JP 2001067150A JP 2001067150 A JP2001067150 A JP 2001067150A JP 4532005 B2 JP4532005 B2 JP 4532005B2
- Authority
- JP
- Japan
- Prior art keywords
- subject
- ray
- image
- ray source
- current
- 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
Links
- 238000000034 method Methods 0.000 title claims description 49
- 230000005540 biological transmission Effects 0.000 claims description 95
- 230000008859 change Effects 0.000 claims description 76
- 238000009826 distribution Methods 0.000 claims description 62
- 238000003384 imaging method Methods 0.000 claims description 43
- 238000001514 detection method Methods 0.000 claims description 12
- 238000012545 processing Methods 0.000 claims description 10
- 230000001678 irradiating effect Effects 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 description 30
- 238000010586 diagram Methods 0.000 description 15
- 230000008569 process Effects 0.000 description 12
- 238000013170 computed tomography imaging Methods 0.000 description 11
- 238000004458 analytical method Methods 0.000 description 9
- 230000006870 function Effects 0.000 description 7
- 210000000056 organ Anatomy 0.000 description 7
- 230000005855 radiation Effects 0.000 description 7
- 210000001015 abdomen Anatomy 0.000 description 6
- 210000004072 lung Anatomy 0.000 description 5
- 238000001228 spectrum Methods 0.000 description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 210000000988 bone and bone Anatomy 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000002591 computed tomography Methods 0.000 description 3
- 210000004185 liver Anatomy 0.000 description 3
- 238000002360 preparation method Methods 0.000 description 3
- 210000001519 tissue Anatomy 0.000 description 3
- 238000010521 absorption reaction Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000012886 linear function Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000001959 radiotherapy Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 238000002083 X-ray spectrum Methods 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 210000001185 bone marrow Anatomy 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 210000001035 gastrointestinal tract Anatomy 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Images
Classifications
-
- 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/46—Arrangements for interfacing with the operator or the patient
- A61B6/461—Displaying means of special interest
- A61B6/463—Displaying means of special interest characterised by displaying multiple images or images and diagnostic data on one display
-
- 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/032—Transmission computed tomography [CT]
-
- 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/46—Arrangements for interfacing with the operator or the patient
- A61B6/461—Displaying means of special interest
- A61B6/466—Displaying means of special interest adapted to display 3D data
-
- 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/46—Arrangements for interfacing with the operator or the patient
- A61B6/467—Arrangements for interfacing with the operator or the patient characterised by special input means
- A61B6/469—Arrangements for interfacing with the operator or the patient characterised by special input means for selecting a region of interest [ROI]
-
- 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/48—Diagnostic techniques
- A61B6/488—Diagnostic techniques involving pre-scan acquisition
-
- 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/54—Control of apparatus or devices for radiation diagnosis
- A61B6/542—Control of apparatus or devices for radiation diagnosis involving control of exposure
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Pathology (AREA)
- Heart & Thoracic Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Physics & Mathematics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Veterinary Medicine (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Human Computer Interaction (AREA)
- Pulmonology (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Description
【０００１】
【発明の属する技術分野】
本発明は、撮影中のＸ線管電流（以下、管電流と略称する）を制御して被検体の被曝線量を抑制するＸ線ＣＴ装置に係り、特に被検体の被曝線量と画質とを考慮して、撮影中の管電流の変化曲線を再設定できるＸ線ＣＴ装置に関する。
【０００２】
【従来の技術】
従来のＸ線ＣＴ装置は、同一断層面では、同一のＣＴスキャン（以下、スキャンと略称する）条件（Ｘ線管電圧（以下、管電圧と略称する）や管電流等）で、撮影を行うようにしている。また、近年被検体を螺旋状にスキャンして撮影するヘリカルスキャンが広く用いられるようになってきているが、体軸方向のスキャン条件もスキャン中一定である。
【０００３】
従って、例えば被検体の断面が、ＣＴスキャナ（以下、スキャナと略称する）の回転軸に対して同心円ではなく楕円である場合、Ｘ線源の回転角度位置によって、被検体におけるＸ線の透過長が大きく変化するため、同一断層面内で、透過するＸ線量の過不足が発生する問題点を有していた。
【０００４】
また、肺等の胸部の低密度の臓器と、肝臓等の腹部の高密度の臓器とでは、Ｘ線の吸収係数が大きく異なるため、胸部から上腹部へ連続的にスキャンを行う場合、肺に適するようなＸ線量を設定すると、肝臓では不足し、肝臓に適するようなＸ線量を設定すると、肺では過剰となる事態が生じていた。
【０００５】
透過するＸ線量が不足する場合、Ｘ線検出器（以下、検出器と略称する）によって検出されるＸ線光子の量の減少によりＳ／Ｎ（ＳＮ比）が悪化し、結果として画像再構成によって得られる断層画像全体のＳ／Ｎが悪化する。逆に、透過するＸ線量が多すぎる場合には、被検体に対して、無効なＸ線被曝がなされていることになる。
【０００６】
これらの問題点を解決する方法として、特開昭53−110495号公報で開示された管電圧を制御する方法や、特開平9−108209号公報、特開平10−309271号公報で開示された管電流を制御する方法が提案されている。
【０００７】
【発明が解決しようとする課題】
しかしながら、特開昭53−110495号公報の管電圧を制御する方法は、スキャン中に管電圧を変化させるためにＸ線のスペクトルが変化し、ＣＴ値が決定できないという問題点がある。このため、現在では、管電流を制御する方法が主流となっている。
【０００８】
管電流を被検体に応じて最適に制御しようとする方法として、特開平10−309271号公報のように、スキャナ回転の半周期前の透過Ｘ線量データを用いて管電流を制御する方法と、特開平9−108209号公報のように、異なる2方向から撮影したスキャノグラムをもとにして、予め管電流を制御するパターンを被検体の位置に応じて決定しておく方法がある。
【０００９】
しかし、特開平10−309271号のスキャナ回転の半周期前の透過Ｘ線量データを用いる方法は、特にヘリカルスキャンにて、スキャンピッチを大きくする場合に、透過Ｘ線量データのずれが大きくなる問題点がある。また、横隔膜前後のように、被検体のＸ線吸収特性が大きく変わる領域では、対応することができない。
【００１０】
特開平9−108209号の異なる2方向からのスキャノグラムを取得する方法は、スキャノグラム撮影を2回行うことにより、被検体への無用なＸ線被曝を増加させることになり、管電流制御による被曝線量の低減の目的に相反するものである。
【００１１】
また、本発明の発明者達も、管電流を被検体に応じて最適に制御する一方法を特願2000−100501号にて提案している。特願2000−100501号の発明は、被検体への無用なＸ線被曝を抑えて、被検体の低被曝化を実現したＸ線ＣＴ装置に関するもので、被検体に関するスキャナの回転角度とＸ線透過長との関係を示すモデルをメモリに格納しておき、被検体に対してのスキャン計測時に、このモデルから定まるスキャナの回転角度毎の設定管電流によるＸ線を、同じ被検体に照射してスキャン計測を行い、断層画像を再構成するものである。
【００１２】
特願2000−100501号の発明では、被検体のＸ線透過長モデルの生成と、このＸ線透過長モデルのＸ線透過長に基づいて管電流を設定することに重点が置かれ、設定された管電流によって被検体に照射されるＸ線量や、被検体内の臓器へのＸ線被曝などについては考慮されていなかった。
【００１３】
上記の問題点を考慮し、本発明では、被検体のＸ線透過長モデルから自動的に設定した管電流の変化パターンによる被検体の内外における被曝線量を算定し、この被検体のＸ線被曝を考慮に入れて操作者が管電流の変化パターンを再設定することができるＸ線ＣＴ装置を提供することを目的とする。
【００１４】
【課題を解決するための手段】
上記目的を達成するため、本発明のX線CT装置は、被検体にX線を照射するX線源と、X線源に高電圧及び電流を供給する高電圧発生装置と、X線源と対向配置され被検体の透過X線量データを検出するX線検出器と、被検体の断層像を得るためのスキャン条件を設定するスキャン条件設定手段と、被検体のスキャノグラム像を表示する表示装置と、表示されたスキャノグラム像上に被検体の断層像の撮影位置を設定する撮影位置設定手段と、設定された撮影位置を含む被検体の周囲に設定されたスキャン条件によりX線源を回転しながらそのX線源からX線を照射し、X線検出に検出された被検体の透過X線データから断層像を再構成する断層像再構成手段と、再構成された断層像を表示装置に表示させる制御手段とを備えたX線CT装置において、スキャノグラム像に基づき被検体のX線透過長を計算する透過長計算手段と、計算されたX線透過長とスキャン条件とから被検体の撮影部位に応じたX線をX線源に発生させるための電流値を高電圧発生装置に設定する電流設定手段とを備えたことを特徴とする。
【００１５】
また本発明のX線CT装置は、被検体にX線を照射するX線源と、X線源に高電圧及び電流を供給する高電圧発生装置と、X線源と対向配置され被検体の透過X線量データを検出するX線検出器と、被検体の断層像を得るためのスキャン条件を設定するスキャン条件設定手段と、被検体のスキャノグラム像を表示する表示装置と、表示されたスキャノグラム像上に被検体の断層像の撮影位置を設定する撮影位置設定手段と、CTスキャナのスキャン条件を設定する条件設定手段と、設定された撮影位置を含む被検体の周囲に設定されたスキャン条件によりX線源を回転しながらそのX線源からX線を照射し、X線検出に検出された被検体の透過X線データから断層像を再構成する断層像再構成手段と、再構成された断層像を表示装置に表示させる制御手段とを備えたX線CT装置において、スキャノグラム像に基づき被検体のX線透過長を計算する透過長計算手段と、計算されたX線透過長とスキャン条件とから被検体の撮影部位に応じたX線をX線源に発生させるための電流値に応じた線量を計算する線量計算手段とを備え、制御手段は、計算した線量をスキャノグラム像と並置して、又は重ねて表示装置に表示させることを特徴とする。
【００１６】
また本発明のX線CT装置は、被検体にX線を照射するX線源と、X線源に高電圧及び電流を供給する高電圧発生装置と、X線源と対向配置され被検体の透過X線量データを検出するX線検出器と、被検体の断層像を得るためのスキャン条件を設定するスキャン条件設定手段と、被検体のスキャノグラム像を表示する表示装置と、表示されたスキャノグラム像上に被検体の断層像の撮影位置を設定する撮影位置設定手段と、CTスキャナのスキャン条件を設定する条件設定手段と、設定された撮影位置を含む被検体の周囲に設定されたスキャン条件によりX線源を回転しながらそのX線源からX線を照射し、X線検出に検出された被検体の透過X線データから断層像を再構成する断層像再構成手段と、再構成された断層像を表示装置に表示させる制御手段とを備えたX線CT装置において、スキャノグラム像に基づき被検体のX線透過長を計算する透過長計算手段と、計算されたX線透過長とスキャン条件とから被検体の撮影部位に応じたX線をX線源に発生させるための電流値に応じた被検体内の線量分布を計算する線量分布計算手段とを備え、制御手段は、計算した線量分布をスキャノグラム像と並置して、又は重ねて表示装置に表示させることを特徴とする。
【００１７】
また本発明のX線CT装置は、被検体にX線を照射するX線源に高電圧と電流を供給する高電圧発生装置と、X線源と対向配置され被検体の透過X線を検出するX線検出器と、X線源とX線検出器を被検体の周囲を回転させるスキャナを回転駆動するスキャナ駆動手段と、被検体の断層像を得るためのスキャン条件を設定する操作手段と、被検体のスキャノグラム像を得て、そのスキャノグラム像を解析し、被検体のX線透過長を求めるスキャノグラム解析手段と、スキャン条件とX線透過長とから電流値を高電圧発生装置に設定する電流設定手段と、スキャナ駆動手段にスキャナを回転させながら、設定された電流値に基づきX線源にX線を照射させて被検体の断層像を得るためのデータを収集する制御装置と、収集されたデータを再構成して断層像を再構成する画像処理装置と、再構成された断層像を表示する表示装置とを備えたことを特徴とする。
【００１８】
また本発明のX線CT装置の画像表示方法は、被検体の断層像を得るためのスキャノグラム像を得て、そのスキャノグラム像を表示装置に表示するステップと、表示されたスキャノグラム像と対応づけてX線源への第1の供給電流の変化のパターンを表示するステップと、そのスキャノグラム像を参照しながら第1の供給電流の供給電流の変化パターンから第2の供給電流の変化パターンへ操作手段によって編集するステップと、編集された第2の供給電流の変化パターンに基づいて被検体の断層像を得て、その断層像を表示装置に表示するステップとを備える。
【００１９】
また本発明のX線CT装置の画像表示方法は、X線源への供給電流を含む被検体の断層像を得るためのスキャン条件を設定するステップと、被検体のスキャノグラム像を撮影して、そのスキャノグラム像を表示装置に表示するステップと、表示されたスキャノグラム像を解析して被検体のX線透過長データを求めるステップと、X線源への供給電流と被検体のX線透過長データとから被検体の撮影部位に応じた電流パターンを計算するステップと、計算された電流パターンに基づく電流値を高電圧発生装置に設定し、その電流をX線源に供給し、その電流とスキャン条件によって被検体の断層像を得て、その断層像を表示装置に表示するステップとを備える。
【００２０】
また本発明のX線CT装置の画像表示方法は、X線源への供給電流を含む被検体の断層像を得るためのスキャン条件を設定するステップと、被検体のスキャノグラム像を撮影して、そのスキャノグラム像を表示装置に表示するステップと、表示されたスキャノグラム像を解析して被検体のX線透過長データを求めるステップと、X線源への供給電流と被検体のX線透過長データとから被検体の撮影部位に応じた電流パターンと被検体内の線量分布を計算するステップと、計算された電流パターンと被検体内の線量分布に基づく電流値を高電圧発生装置に設定し、その電流をX線源に供給し、その電流とスキャン条件によって被検体の断層像を得て、その断層像を表示装置に表示するステップとを備える。
【００２７】
【発明の実施の形態】
以下、添付図面を用いて、本発明の実施例について説明する。
図1は、本発明に係るＸ線ＣＴ装置の全体構成を示すブロック図である。図1に示すように、このＸ線ＣＴ装置は、主としてＸ線源12と検出器13等を搭載し、被検体15に対してその周囲を連続回転可能なスキャナ11を内蔵するガントリ10と、装置全体を総括するホストコンピュータ20と、Ｘ線源12に高電圧を供給する高電圧発生装置22と、画像データの前処理や画像再構成処理、あるいは各種の解析処理を行う画像処理装置24と、画像を表示する表示装置25と、被検体15を載せるテーブル装置18と、操作者がスキャン条件等を入力する操作手段21等から成る。なお、スキャナ11と被検体15とは、相対的に回転することができればよいので、被検体15が静止してスキャナ11が回転してもよいし、スキャナ11が静止して被検体15の方が回転するとしてもよい。
【００２８】
図2は、本発明に係るＸ線ＣＴ装置の要部構成要素を示す。図2を用いて、先ず、スキャナ11の詳細について説明する。図2において、スキャナ11には、Ｘ線源12と検出器13とが180度対向した位置関係で配置されている。このＸ線源12から放射されたＸ線ビーム14は、コリメータ19によってビームの幅及び厚さが制限されたファン状のＸ線ビーム14となって、被検体15に照射される。Ｘ線源12は、高電圧発生装置22を介してホストコンピュータ20によって制御される。スキャナ11全体は、スキャナ角度検出手段17によって回転角度を検出し、検出した回転角度に基づいてホストコンピュータ20が、スキャナ駆動手段16を制御し、スキャナ11を駆動する。検出器13は、被検体15を透過したＸ線14を検出し、検出データは被検体15によるＸ線の減衰量を示す投影データとして取り込まれる。投影データは、画像処理装置24において、ホストコンピュータ20の持つスキャナ角度等のデータと照合され、画像再構成等の処理をされた後に、表示装置25にて断層画像として表示される。
【００２９】
次に、本発明に係るスキャナ11以外の要部構成要素について説明する。図2において、装置全体を総括するホストコンピュータ20には、スキャナ11はスキャナ駆動手段16を介して、Ｘ線源12は高電圧発生装置22を介して、検出器13は画像処理装置24を介して、それぞれ間接的に接続されており、操作手段21と、管電流設定手段22と、画像処理装置24と、スキャノグラム解析手段26は直接接続されている。ホストコンピュータ20とスキャナ11、Ｘ線源12、検出器13との接続により、ホストコンピュータ20はＸ線源12による被検体15へのＸ線照射と、検出器13による投影データ（検出データ）の取り込みを制御する。画像処理装置24はホストコンピュータ20の指令に応じて、取り込まれた投影データに基づき、断層画像を順次再構成する。
【００３０】
本発明に係るＸ線ＣＴ装置では、被検体の断層画像を取得する本スキャンの前に、スキャン条件を設定するために種々の準備操作を行う。この準備操作としては、被検体の位置決めのためのスキャノグラム画像の撮影、管電流設定のためのスキャノグラム画像データの解析、スキャン条件としての管電流の変化パターンの決定などが、ホストコンピュータ20の介在のもとで行われる。
【００３１】
これらの準備操作に関与する主な構成要素としては、図2において、ホストコンピュータ20と、操作手段21と、スキャノグラム解析手段26と、管電流設定手段22と、Ｘ線源12と、検出器13などである。この準備操作において、先ず、操作手段21は主として管電流の設定値(最大値、最小値)などのスキャン条件をシステムに入力する。Ｘ線源12と検出器13はスキャナ11を回転させずに、スキャノグラム画像の撮影を行い、画像データをホストコンピュータ20に保存する。スキャノグラム解析手段26はスキャノグラム画像データを解析し、被検体のＸ線透過長を体軸方向のスライス位置毎及びスキャナの回転角度毎に算出可能な、3次元形状データとしてモデル化し、このモデル（以下、被検体の3次元的Ｘ線透過長モデルという）のデータをホストコンピュータ20に保存する。管電流設定手段22は操作手段21から入力された管電流設定値と被検体の3次元的Ｘ線透過長モデルのデータを基にして、スキャン中に被検体の撮影部位のＸ線透過長の変化に応じて経時的に変化する一連の管電流値すなわち、管電流の変化パターンを自動的に決定する。このように決定された管電流の変化パターンは、ホストコンピュータ20に保存され、本スキャン時に被検体の撮影部位に応じて順次呼び出されて、Ｘ線源12の管電流を変化させる。
【００３２】
本発明では、更に上記で決定された管電流の変化パターンを基に、被検体15に照射されるＸ線量を本スキャン前に予め計算する。Ｘ線装置では、通常照射したＸ線量に対応する量としてmAs値が使用されているが、この計算でもmAs値を採用している。このmAs値は、管電流（mA）と照射時間（ｓ）の積であり、管電圧が一定である場合（Ｘ線ＣＴ装置では管電圧が一定で使用される場合が多い）には、Ｘ線源12より照射されるＸ線量の総和に比例するために、Ｘ線量の基準として用いられる。ここで計算された被検体15に照射されるＸ線量については被検体15への予測被曝線量として操作者によって評価されることになる。
【００３３】
図3に、本発明に係るＸ線ＣＴ装置を使用した第1のスキャン操作の一連の動作のフローチャートを示す。このスキャン操作では、ステップ103、106、108の3次元データ生成、ステップ110の管電流パターン生成、ステップ111のmAs計算、ステップ114のmAs表示に特徴がある。以下、図2を参照しながら、図3の第1のスキャン操作の詳細について説明する。
【００３４】
図3において、先ず、ステップ101のスキャノグラム撮影の工程では、被検体15のスキャノグラム画像を撮影する。被検体15のスキャノグラム画像を撮影する構成と断層画像を撮影する構成とは基本的には同じである。本ステップでは、スキャノグラム画像データは、スキャナ11を回転させずに、被検体15に正面方向からＸ線14を照射して、検出器13によって検出データを取り込むことによって得られる。このとき得られるスキャノグラム画像は正面方向のものである。このスキャノグラム画像データは、検出器13からホストコンピュータ20に送られる。このスキャノグラム画像データは、本スキャン時の被検体15の位置決めのために利用される他、本発明では特に管電流制御のための管電流の変化パターンの決定のために利用される。
【００３５】
次に、ステップ102のスキャノグラムデータ解析の工程及びステップ103の第1の3次元データ生成の工程では、スキャノグラム画像データがホストコンピュータ20に接続されたスキャノグラム解析手段26によって解析され、被検体15の3次元的Ｘ線透過長モデルが生成される。この3次元的Ｘ線透過長モデルは、被検体15をＣＴ撮影する場合の、被検体15の位置とＸ線透過長との関係を示すモデルである。被検体15の3次元的Ｘ線透過長モデルの作成方法については特願2000−100501号にも開示されている。
【００３６】
以下、被検体15の3次元的Ｘ線透過長モデルの作成方法の一例について説明する。図4はスキャノグラム画像とスライス位置及び投影データ例との対応を示す図、図5は3次元的Ｘ線透過長モデルの1スライス位置でのモデルを示す図である。図4（a）は、ステップ101で撮影された被検体のスキャノグラム画像29を示す。このスキャノグラム画像29は胸部から腹部の中間位置までの領域を撮影領域としている。このようなスキャノグラム画像の撮影領域の中からスライス位置が選ばれる。図示の場合n個のスライス位置が選ばれている。図中、P1、…、Pi、…、Pj、…、Pnがスライス位置である。
【００３７】
図4（b），（c），（d）は3次元的Ｘ線透過長モデルのモデル決定の説明図である。任意の2つのスライス位置Pi、PjでのＣＴ断層像が図4（b）の如くなっているものと仮定すると、その縦方向（図示の上下方向）のＸ線減衰量の投影データは図4（c）の如くなるはずである。人体の体幹部の横断面は通常楕円形に近いものであるので、任意のスライス位置Pi、PjのＣＴ断層像を図4（b）で仮定するのは、大きな誤差はないものと判断される。そこで、図4（c）の投影データについては、Ｘ線透過長データに換算し、その後横軸に沿って積分して面積を求める。このとき、Ｘ線減衰量の投影データのＸ線透過長への換算にあたっては、簡単のため人体が水と等価であるとみなしてデータの変換を行う。Ｘ線減衰量をc、Ｘ線透過長をb、水の線減弱係数をμwとしたとき、両者の関係はb＝log c／μwで表される。また、横軸については、Ｘ線減衰量データの存在する領域全体の幅が人体の幅寸法に一致するように変換する。図4（d）は、図4（c）の投影データから変換したスライス位置Pi、Pjにおける被検体15のＸ線透過長データの分布図である。図4（d）においてスライス位置Pi、Pjにおける最大Ｘ線透過長はbi、bj、面積はSi、Sjとなる。図4（ｄ）のＸ線透過長データについて、最大Ｘ線透過長bi、bj及び面積Si、Sjに注目すると、bi、Siはスライス位置Piでの断層画像のＸ線透過状況を反映し、bj、Sjはスライス位置Pjでの断層画像のＸ線透過状況を反映した値と見なすことができる。
【００３８】
そこで、被検体15の3次元的Ｘ線透過長モデルとして、各スライス位置でのスライス断面を図5に示すような楕円形30でモデル化することにした。このモデル化では、スライス位置Pi、Pjにおける楕円形モデル30i、30jの面積をSi、Sj、短軸をbi、bjとしている。この結果、楕円形モデル30i、30jの長軸をai、ajとすると、楕円形モデル30i、30jの面積が［数1］で表されることから、
【数１】
長軸ai、ajは［数2］によって求められる。
【数２】
【００３９】
上記によって、各スライス位置での断層画像に対応するＸ線透過長モデルとしての楕円形モデル30が求められたので、これらの楕円形モデル30を体軸方向に配列することによって3次元的Ｘ線透過長モデル30を作成することができる。体軸方向のスライス位置のピッチが粗いときには、例えば隣り合う楕円形モデル間で最小自乗法により、その途中の1つ又は2つ以上の楕円形モデルを補間で求める。以上の如き手順により被検体15の3次元的Ｘ線透過長モデル30のデータとして、3次元座標（X、Y、Z）系での被検体15のＸ線透過長データT＝T（X、Y、Z）が生成される。
【００４０】
次に、ステップ104からステップ110の工程では、上記の3次元的Ｘ線透過長モデル30を使用してＸ線源12に付与する管電流の変化パターンを設定することになるが、その前に3次元的Ｘ線透過長モデル30を使用した管電流の求め方について説明する。図4、図5で求めた3次元的Ｘ線透過長モデル30は、被検体の各スライス位置での断層画像のＸ線透過長を反映したものである。3次元的Ｘ線透過長モデル30のデータは一旦ホストコンピュータ20のレジスタを含むメモリに格納されているので、撮影範囲及びテーブルピッチなどのスキャン条件が決定するとその範囲のモデルのデータがメモリから取り出され、第2、第3の3次元データの生成及び管電流の変化パターンの決定に使用される。
【００４１】
管電流は、管電流設定手段23によって各スライス位置でスキャナ回転角度毎に3次元的Ｘ線透過長モデル30から得られるＸ線透過長に基づいて決定される。図6には、スライス位置（体軸方向の位置）ZにおけるＸ線透過長モデル30を示す。また、あるスキャナ回転角度における管電流は、通常そのスキャナ回転角度における3次元的Ｘ線透過長モデルのＸ線透過長のうちの最大値に対応付けて決定される。この最大値を示すＸ線透過長は図6の楕円形モデル30の中心0を通過するパスで得られるので、管電流を設定するにあたっては、スキャナ回転角度毎にこの楕円形モデル30の中心0を通過するパスのＸ線透過長のみを考慮すればよい。従って、図6においてスライス位置をZ、スキャナ回転角度をθ（θの始点は楕円形モデル30の短軸方向とする）としたとき、その位置における最大Ｘ線透過長Tは、Zとθの関数としてT＝T（Z、θ）と表すことができる。
【００４２】
この最大Ｘ線透過長T（Z、θ）は、楕円形モデル30の中心位置0を通るパスの長さであるので、長軸をa、短軸をb、スキャナ回転角度をθとした場合、［数3］の如く表すことができる。
【数３】
ここで、a、bは［数1］、［数2］のai、aj、及びbi、bjと対応する。
【００４３】
次に、管電流の設定方法の一例について説明する。先ず、被検体をスキャンする全範囲におけるパスの最大値（全スライス位置P1〜Pnの中でのパスの最大値）をTmax、最小値（同じくパスの最小値）をTminとする。これらの値は3次元的Ｘ線透過長モデル30を作るとき既知である。管電流を最大値Imax（mA）と最小値Imin（mA）の範囲で変化させる場合、本例では管電流の最大値、最小値とパスの最大値、最小値をそれぞれ対応させて、管電流とパスとの間に直線関係を持たせるものである。管電流IとパスTとの関係は［数4］の如く表される。
【数４】
ここで、パスTはT(Z、θ)に対応するので、管電流Iは(Z、θ)の1次関数となり、スライス位置Z及びスキャナ回転角度θごとに求められる。
【００４４】
次に、図3のフローチャートに戻って説明する。ステップ104及びステップ105の工程では、操作者がスキャノグラム画像を参照して操作手段21からスキャン条件としてのテーブルピッチ及びスキャン開始位置を入力する。これらのデータにより、被検体のＣＴ撮影範囲とスライス位置とスキャン回転角度が決定される。このときの座標系としては、上記の如く(Z、θ)座標系がよく、スキャン条件のデータも、(Z、θ)座標系のデータで入力するのがよい。
【００４５】
次に、ステップ106の工程では、第2の3次元的Ｘ線透過長モデルのデータを生成する。この工程で生成するデータは、各スライス位置Z、スキャン回転角度θ毎の最大Ｘ線透過長であり、第1の3次元的Ｘ線透過長モデルのデータから［数3］によって求めることができるので、ホストコンピュータ20のメモリから第1の3次元的Ｘ線透過長モデルのデータを呼び出して演算する。この演算結果は、T＝T(Z、θ)で表される。
【００４６】
次に、ステップ107の工程では、スキャン条件としてのスキャン時間を操作手段21より入力する。スキャン開始位置と、テーブルピッチと、スキャン時間が決定すると、スキャン中のＸ線源12の位置(Z、θ)はスキャン開始後の経過時間tの関数として表すことができるので、各スキャン位置での被検体15の第2の3次元的Ｘ線透過長モデル、すなわち最大Ｘ線透過長Tも時間tの関数T＝T(t)として表すことができる。このため、ステップ108の第3の3次元的Ｘ線透過長モデル生成の工程では、最大Ｘ線透過長Tの関数を、T＝T(Z、θ)からT＝T(t)に変換する。
【００４７】
次に、ステップ109の管電流設定値入力の工程では、操作者が操作手段21より、管電流の設定値、例えばスキャン中の管電流の最大値Imaxと最小値Iminを入力する。ステップ110の管電流パターン生成の工程では、管電流設定手段23がホストコンピュータ20から3次元的Ｘ線透過長モデルのデータT(t)を呼び出し、上記の管電流設定値に基づき、被検体15の撮影部位に応じた管電流の変化パターンを自動的に決定する。このとき、Ｘ線透過長T(t)に対応付けてスキャン中の管電流の値を設定することになるが、Ｘ線透過長T(t)が最小の時には最小の管電流を設定し、Ｘ線透過長T(t)が最大の時には最大の管電流を設定するように管電流の変化パターンを決定する。また、Ｘ線透過長T(t)と管電流の値との関係として［数4］に示した1次関数の他、種々のものがある(但し、［数4］のTはT(t)に対応する)。
【００４８】
上記の如くして被検体15の3次元的Ｘ線透過長モデルに合わせて、管電流が時間tの関数として決定される。従って、管電流の変化パターンは、I＝I（t）と表すことができる。このように決定された管電流の変化パターンI＝I（t）はホストコンピュータ20に保存され、本スキャン時に被検体15の撮影部位に応じて順次呼び出されて、高電圧発生装置22を介してスキャン中の管電流を制御する。
【００４９】
図7に、管電流の変化パターンの表示例を示す(本例も含めて以下の管電流の変化パターンの表示例では、簡単のためスキャナ回転角度θの変化に伴う大略周期的な管電流の変化については省略し、スライス位置による変化のみ示している)。これは、表示装置25の画面上にスキャノグラム画像29と対比して表示したものである。管電流の変化パターン31では、縦軸に管電流値（mA）、横軸にスキャン開始後の経過時間tをとっている。表示例の場合、管電流はスキャンの初期（腹部）には中程度の値で、中期（腹部と胸部との間）には最小値となり、終期（胸部）には最大値となっている。本表示例の如く、管電流の変化パターン31とスキャノグラム画像29を同一画面上に並置することにより、管電流値と撮影部位との対比を一目で行うことができるので、管電流値の妥当性の判断に有効である。
【００５０】
次に、ステップ111のmAs計算の工程では、ステップ110で決定された管電流の変化パターンに基づき、スキャン中に被検体15に照射されるＸ線量を計算する。ここで、被検体15に照射されるＸ線量の基準としては上述の如く管電流（mA）と照射時間（s）との積であるmAsが用いられる。従って、本ステップでは、管電流の変化パターンI＝I（t）を時間で積分して、被検体15に照射されるＸ線量mAsを求める。この積分は［数5］によって行われる。
【数５】
ここで求めたmAs値は、あくまでも被検体15に照射されるＸ線量に相当する量であるので、実験等によってＸ線量とmAs値との正確な対応をとり、両者の間の換算ができるようにしておく必要がある。
【００５１】
図7に示した管電流の変化パターン31を例に上げて、被検体15への照射Ｘ線量mAsを計算することを考える。この場合、管電流の変化パターン31は、I＝I（t）で、時間の関数であるから、これを積分してmAsを計算することは、管電流の変化パターン31の図の面積を計算することになり、管電流の変化パターン31の図の面積SがmAsに相当する。
【００５２】
次に、ステップ112のmAs計算値表示の工程では、ステップ111にて求めたmAsの計算値を表示装置25の画面に表示する。ステップ111では、管電流の変化パターンが生成された段階で、被検体15のスキャン範囲の全領域にわたって被検体に照射されるＸ線量に対応するmAsが計算されるので、このステップ112でこのmAs値がスキャンを開始してもよいか否かを操作者が判断するための資料として操作者に提示される。
【００５３】
ステップ113のmAs判断の工程では、操作者が全体としてのmAsの妥当性について判断する。すなわち、操作者はこれから行われるＣＴ撮影による利益と被検体15へのＸ線被曝による不利益とを比較考量して、全体としてのmAsが大きすぎないか否かを判断し、全体としてのmAsが大きすぎると判断した場合には管電流の設定値を下げることになる。この場合には、ステップ109に戻り、管電流設定値を再度入力し、管電流の変化パターンを再設定する。
【００５４】
図8には、管電流の変化パターンの例を示す。図8（a）に示した管電流の変化パターン31aは、通常の管電流一定（I0）の場合で、管電流の最小値（Imin）と最大値（Imax）が同じI0である場合の例である。図8（a）の管電流の変化パターン31aのグラフでは、mAsは面積Saである。次に、図8（b）示した管電流の変化パターン31bでは、管電流は初期には図8（a）と同じI0で、中期にI0より低い最小値Iminとなり、終期にI0より高い最大値Imaxとなっている。図8（b）の管電流の変化パターン31aのグラフでのmAsは面積Sbであるが、この面積Sbは図8（a）の面積Saより小さくなっており、被検体15の被曝線量は低減している。
【００５５】
図8（c）に示した管電流の変化パターン31cでは、図8（b）のグラフと比べて、終期の最大値Imaxを小さくし、I0とほぼ同じとしたものである。この場合、スキャン全体として管電流が低目に抑えられている。図8（c）の管電流の変化パターン31cのグラフでのmAsは面積Scであるが、この面積Scは図8（b）の面積Sbより更に小さくなっており、被検体15の被曝線量は更に低減している。
【００５６】
図8（d）に示した管電流の変化パターン31dでは、図8（b）に比べて中期の最小値Iminを更に小さくし、IminがI0に対し大幅に低下するようにしたものである。図8（d）の管電流の変化パターン31dでのmAsは面積Sdであるが、図8（c）の面積Scより更に小さくなっている。
【００５７】
図8（c）の場合の如く管電流設定値の最大値を下げた場合には、被検体15の厚さが厚い部分、すなわちＸ線透過長が大きい部分での管電流を減らすことになる。従って、図8（c）の管電流の変化パターン31cは、特に肺のように低密度の部位での画質を重視し、腹部のように高密度の領域での被曝線量を低減したい場合に適する。
【００５８】
図8（d）の場合の如く管電流設定値の最小値を下げた場合には、被検体15の厚さが薄い部分、すなわちＸ線透過長が小さい部分での管電流を減らすことになる。従って、図8（d）の管電流の変化パターン31dは、特に骨周辺や実質部のように高密度の部位での画質を重視し、低密度の領域での被曝線量を低減したい場合に適する。
【００５９】
図8（c），図8（d）では、管電流設定値の最大値Imax又は最小値Iminのいずれか一方のみを小さくしているが、両方を小さくして、被検体15の全領域について平均的に被曝線量を低減することも可能である。上記の如く新たに設定した管電流の変化パターンについては、再度mAsを計算し、問題が無い場合には、そのまま採用する。
【００６０】
次に、ステップ114のスキャンの工程では、操作者は上記で決定した管電流の変化パターンを含めたスキャン条件にてスキャンを実行する。
【００６１】
次に、ステップ115のmAs積算値表示の工程では、上記の［数５］に従って、スキャン中のmAsの積算値を逐次計算し、表示装置25の画面にリアルタイムで表示し、操作者に提示する。mAs積算値表示方法としては、全体としてのmAs値に対する比率としての相対値を表示する方法と、mAs積算値の絶対値を表示する方法のいずれかを選択することができる。もちろん、その両方を同時に表示することも可能である。
【００６２】
次に、本発明に係るＸ線ＣＴ装置の第2のスキャン操作例について説明する。図9は、第2のスキャン操作例の一連の動作のフローチャートを示したものである。このスキャン操作例では、図3の第1のスキャン操作例のフローチャートに対し、ステップ111のmAs計算及びステップ112のmAs表示の工程の後に、被検体15内の線量分布の計算及び表示の工程を追加し、操作者が被検体15内の線量分布をも見てスキャン実行の判断を下せるようにしたものである。このため、本操作例の説明では、ステップ201〜205の被検体15内の線量分布の計算及び表示の工程を重点に説明する。
【００６３】
以下、図9の第2のスキャン操作例の中のステップ201からステップ205について説明するが、ステップの内容説明に入る前に、被検体15内の線量分布の計算結果の表示例を図10に示す。図10（a）は被検体の横断面35での線量分布の表示例、図10（b）は被検体の体軸方向の側面36での線量分布の表示例である。両図とも被検体15内の分布線量の等しい等線量線38a〜38c、39a〜39cが示されており、体表に近いほど高線量になっている。本操作例では、被検体15の線量分布が操作者に提示されるため、操作者は被検体のＸ線被曝についてより詳細な評価を行うことができるということが大きな特徴となる。
【００６４】
本操作例では、ステップ203で被検体内の線量分布の計算を行う前に、準備として、ステップ201で被検体のＣＴ値モデルの生成、ステップ202で被検体のμモデルの生成を行い、その後で被検体のμモデルデータと被検体に照射されたＸ線量データを基にして被検体内の線量分布を計算する。
【００６５】
先ず、ステップ201の被検体ＣＴ値モデル生成の工程では、予め標準的な人体のＣＴ値分布モデル（以下、標準人体ＣＴ値モデルという）データを取得して、ホストコンピュータ20の記憶手段に保存しておき、この標準人体ＣＴ値モデルデータを、ステップ101で取得した被検体15のスキャノグラム画像のデータに基づいて補正を行うことによって、被検体15のＣＴ値モデル（以下、被検体ＣＴ値モデルという）を作成する。上記の標準人体ＣＴ値モデルデータとしては、例えば標準的な人体ファントムなどをＣＴ撮影した断層画像から得られる3次元的ＣＴ値分布データが用いられる。断層画像はＣＴ値の分布を表し、実効エネルギー（通常60keV）のＸ線に対する線減弱係数の分布を表しているので、この断層画像を3次元に再構成をした3次元的ＣＴ値分布データは線減弱係数の3次元的な空間分布のデータであり、被検体に照射されたＸ線の減弱量の計算に利用することができる。
【００６６】
次に、被検体ＣＴ値モデルの生成方法の一例について、図11を用いて説明する。この例では、ステップ101で被検体15を撮影して取得して実測したスキャノグラム画像データと、上記の標準人体ＣＴ値モデルデータとから、被検体15の3次元的ＣＴ値分布を表す被検体ＣＴ値モデルデータを生成する。この被検体ＣＴ値モデルの生成にあたっては、被検体15のスキャノグラム画像データと標準的人体のスキャノグラム画像データが媒体として利用される。
【００６７】
図11は、ステップ201の被検体ＣＴ値モデルの作成手順を説明するための図である。図11において、図11(a)は標準人体ＣＴ値モデルデータ41の例を、図11(b)は図11(a)の標準人体ＣＴ値モデルデータ41から計算で求めた標準人体のスキャノグラム画像データ42の例を、図11(c)は被検体の実測のスキャノグラム画像データ43の例を、図11(d)は計算で求めた被検体ＣＴ値モデルデータ44の例を示す。図11(a)の標準人体ＣＴ値モデルデータ41は、人体ファントムなどの標準的な人体の体幹部のＣＴ値分布モデルで、肩から腹部までのスライス位置ごとの断面のＣＴ値分布モデルを示している。
【００６８】
スキャノグラム画像は、その3次元的ＣＴ値分布モデルから計算によって生成することができるので、図11(b)の標準人体スキャノグラム画像データ42は図11(a)の標準人体ＣＴ値モデルデータ41について正面方向から投影したデータを求めることによって得られる。図11(c)の被検体の実測スキャノグラム画像データ43は、被検体15の体幹部について標準人体スキャノグラム画像データ42と同じ領域を正面方向から撮影したスキャノグラム画像データである。この画像については、以下被検体スキャノグラム画像と呼ぶことにする。
【００６９】
図11においては、体幹部の標準人体スキャノグラム画像データ42と被検体スキャノグラム画像データ43とを並置して対比できるように示してあるが、両者は普通寸法及びＣＴ値とも異なるものである。このため、標準人体スキャノグラム画像データ42と被検体スキャノグラム画像データ43とを対比しながら、両者の差異に基づいて、一致する部分はそのままとし、異なる部分については変形させて、被検体15に合うように標準人体ＣＴ値モデルデータ41を補正して、被検体ＣＴ値モデルデータ44を生成する。
【００７０】
図11の体幹部の例では、先ず体軸方向に関して、標準人体スキャノグラム画像データ42と被検体スキャノグラム画像データ43の肩から横隔膜までの長さAと、横隔膜から腸管までの長さBとに分けて、それぞれの差異を基に、標準人体ＣＴ値モデルデータ41を補間、伸長したり、あるいは間引き、短縮したりすることで、標準人体ＣＴ値モデルデータ41の体軸方向のＣＴ値分布を被検体15の実状に近似させる。左右方向に関しても、同様に、体軸を基準にして左と右に分けて、それぞれの差異を基に、左右の広がりを補正し、被検体15の実状に近似させる。前後方向に関しては、被検体スキャノグラム画像データ43から推定される前後方向のＸ線透過長を基に、標準人体ＣＴ値モデルデータ41の前後方向のデータを線形に補間する。このようにして、2つのスキャノグラム画像データ42、43を基に、標準人体ＣＴ値モデルデータ41を実際の被検体15に合わせ込むことにより、被検体ＣＴ値モデルデータ44を生成する。
【００７１】
次に、ステップ202の被検体μモデル生成の工程では、ステップ201で生成した被検体ＣＴ値モデルデータ44のＣＴ値を線減弱係数μに変換して、被検体15の3次元的μ値分布モデルを生成する。ＣＴ値から線減弱係数μへの変換は下記の如く行われる。
【００７２】
ＣＴ値は実効エネルギー(通常60keVを使用)のＸ線に対する線減弱係数によって決定され、水＝0、空気＝1000、平均的な骨＝1000として定義されている。今、被検体ＣＴ値モデルの位置XにおけるＣＴ値をＣＴXとすると、その位置Ｘにおける実効エネルギー(60keV)での線減弱係数μXは［数6］、［数7］で表される。
【数６】
【数７】
ここで、μwは水の線源弱係数(＝0.206cm-1)、μairは空気の線減弱係数(＝0.00025cm-1)、μboneは骨の線減弱係数(＝0.567cm-1、ただし密度1.8g／cm3の場合)である。
【００７３】
次に、ステップ203の線量分布計算の工程では、図12に示す如く、ステップ111で求めた被検体に照射されたＸ線量のデータ(図12(a))とステップ202で求めた被検体μモデルデータ45(図12(b))を用いて、被検体15内の線量分布(図12(c))を計算する。このステップ203の計算では、被検体に照射されるＸ線のエネルーギースペクトルを考慮して、被検体15内でのＸ線の減衰を計算し、被検体15内の線量の空間的分布を計算する。被検体15に任意の方向からＸ線を照射した場合のＸ線の減弱量は、被検体の3次元的線減弱係数(μ)のモデルである被検体μモデルデータを用いることにより解析的に計算可能であり、このような計算手法は既に他の分野、例えば放射線治療計画装置などの分野でも行われている(参考文献1、稲邑清也、放射線治療計画システム、P．90〜92、P．113〜115、篠原出版、平成4月20日発行)。
【００７４】
Ｘ線の減衰の計算にあたっては、先ず、被検体15内の注目する位置XからＸ線源12に向かってＸ線が透過する実効距離δを計算する。実効距離は、Ｘ線が透過する媒質によって1／ｅに減衰する距離を1と定義される。Ｘ線のエネルギースペクトルを考慮しない場合、注目位置ＸからＸ線源12までに、組成ｉの媒質が実距離でdi(cm)あり、かつ組成ｉの線減弱係数がμi(cm-1)であるとき、実効距離δは［数8］で表される。
【数８】
【００７５】
しかし、Ｘ線のエネルギースペクトルを考慮すると実効距離はＸ線のエネルギーに応じて異なる値となる。組成iのＸ線のエネルギーjに対する線減弱係数をμij(cm-1)とすると、Ｘ線のエネルギーjに対する実効距離δjは［数9］で表される。
【数９】
ここで、μijは組成iのエネルギーjのＸ線に対する線減弱系数(cm-1)、diは組成i中のＸ線の透過距離(cm)である。μijについては、Ｘ線のエネルギーjに応じて被検体15のμ値モデルから求める必要がある。
【００７６】
次に、被検体15内の注目位置Xでの線量を計算する。Ｘ線源12から注目位置Xまでの距離をrx(m)、距離1mにおける線量をI0(C／kg：Cはクーロン)とする。I0については、例えば実験的に求める。Ｘ線のエネルギースペクトル、つまりエネルギーjの成分比をSjとすると、注目位置Ｘでの線量Ix(C／kg)は［数10］で表される。
【数１０】
ここで、I0は線量(Ｘ線源12から単位距離における空中の線量、単位C／kg)、rxはＸ線源12から注目位置Ｘまでの距離、単位m)、SjはＸ線のエネルギースペクトルである。
【００７７】
［数10］による計算の結果、Ｘ線源12がある位置Q(Z、θ)にある場合の被検体15内の任意に位置Ｘにおける線量が求められる。Ｘ線源12を1回転した場合の任意の位置Ｘにおける線量は、Ｘ線源12の位置Q(Z、θ)を被検体15の周囲で回転させ、上記の線量を1回転分(θ=0〜2π)積算することによって得られる。上記の如き手順で、被検体15内の設定された各位置での計算を行うことにより、スライス位置Zでの被検体15内の線量分布を求めることができる。また、他のスライス位置における線量分布も同様に計算によって求めることができるので、被検体15の体軸方向について撮影部位の全領域の計算を進めることにより、被検体15内の3次元的な線量分布が得られる。
【００７８】
1スライスの断層画像を撮影する場合(2次元の場合)には上記の計算によって被検体15内の線量の分布が精度良く求められるが、1回のスキャン中に複数スライスのＣＴ撮影を行う場合(3次元の場合)には散乱Ｘ線の撮影を考慮しないと計算精度が低下する恐れがある。Ｘ線ＣＴ装置では、Ｘ線のエネルギーが100keV以下のオーダーであるため、散乱線としてはコンプトン(Compton)散乱のみ考慮すればよい(参考文献１参照)。このコンプトン散乱を考慮することにより、計算精度をより高くすることができる。
【００７９】
上記の如き手順で、被検体15内の各スライス位置での各設定点における線量を計算することにより、計算された被検体15の3次元的線量分布が得られる。この被検体15の線量分布の計算結果はホストコンピュータ20に一時的に保存され、操作者などにとって見やすい図象、例えば図10に示したような図象で表示される。
【００８０】
次に、ステップ204の線量分布表示の工程ではステップ203の計算結果が表示装置25に表示される。本実施例での表示例としては、図10に示したような被検体15の断層面35内の線量分布(図10(a))又は被検体15の側面36の線量分布(図10(b))が上げられる。これらの図では、被検体15の臓器と線量分布を示す等線量線を重ねて表示しているので、一目で各臓器への被曝線量を認識することができるので、被検体15へのＸ線被曝を評価する上では有効である。
【００８１】
次に、ステップ205の被曝線量判断の工程では、ステップ204で表示された被検体15内の線量分布の計算結果を操作者が見て、被検体15内の臓器へのＸ線被曝が過剰になる恐れがないかどうかを判断し、Yesと判断した場合にはステップ114のスキャン実行の工程に進み、スキャンを開始することになり、Noと判断した場合にはステップ109の管電流設定値入力の工程に戻り、管電流設定値の再入力、管電流パターンの再検討を行うことになる。
【００８２】
また、ステップ201での被検体15のＣＴ値モデルデータの生成方法としては、上記の標準人体ＣＴ値モデルデータを使う方法以外に、過去に撮影した同じ被検体15のＣＴ撮影データを使用する方法も実施可能である。この場合には、実際に同一の被検体15のＣＴ値分布データを使用するために、標準人体ＣＴ値モデルデータ41の形状を補正する手順が必要なくなるという利点がある。しかし初回のＣＴ撮影には適さないため、過去にＣＴ撮影を行った被検体に関して、2回目以降のＣＴ撮影を行う場合が対象となる。
【００８３】
上記したように、ＣＴスキャン前に被検体15の体内の線量分布をシミュレーション計算して、計算結果を、例えば図10に示す如く表示することにより、操作者は事前に撮影手技に応じた被検体の体内の線量分布を近似的に知ることが可能となる。
【００８４】
この結果、例えば、単純に被検体の全ての組織に関して一律に被曝線量を減らすのではなく、骨髄や肺等の放射線感受性の高い組織に関しては特に被曝線量を低減し、逆に脂肪や筋肉等の放射線感受性の比較的低い組織には、画質が満足できる程度に被曝線量のレベルを維持するというような詳細な設定が可能となる。
【００８５】
次に、図13を用いて、被検体をＣＴ撮影する管電流の変化パターンの編集例を説明する。図13は、被検体のスキャノグラム画像上に管電流の変化パターンを重畳して表示したものである。図12において、スキャノグラム画像データ29aは体幹部のもの、管電流の変化パターンは、編集前の初期の管電流の変化パターン46aと編集後の修正された管電流の変化パターン46bである。
【００８６】
この管電流の変化パターンの編集工程では、表示装置25の画面においてスキャノグラム画像29a上に表示された初期に設定された管電流の変化パターン46aに対し、スキャノグラム画像データ29aを参照しながら、また場合によっては被検体15の内部の照射線量分布を参照して、操作手段21によって修正を加えて、新しい管電流の変化パターン46bを編集する。この編集操作によって任意の部位の管電流の変化パターンを再設定する。
【００８７】
この編集操作において、例えば、自動的な管電流の変化パターンの設定では、横隔膜付近のように密度が大きく変化する領域では、管電流を平均的な値に設定するが、被曝線量が増えても画質を向上させる必要がある領域などでは、管電流を部分的に高く設定する。管電流の変化パターンは、上記の如くスキャン条件が設定されていれば、時間tのみの関数になるので、任意時刻の管電流の値を変化させることができる。図13の例では、初期の管電流の変化パターン46aに対し、肺の領域の管電流を少し低下させ、横隔膜の領域の管電流を少し増加させることで、修正後の管電流の変化パターン46bに編集している。
【００８８】
図14は、管電流と被検体の厚さ(Ｘ線透過長に相当)との関係を示した図である。上記のステップ110の工程の説明では、管電流の最大値、最小値と被検体の厚さの最大値、最小値とを一致させて、両者間で線形の関係を持つものとして説明したが、両者の関係については操作者の設定により非線形の関係を持つようにすることが可能である。図14に示した管電流Iと被検体の厚さTとの関係は、［数11］で表される。
【数１１】
ここで、Imax、Iminは管電流の最大値と最小値、Tmax、Tminは被検体の厚さの最大値と最小値、γは定数である。γについては以下ガンマと呼ぶことにする。
【００８９】
図14において、グラフ50はガンマ=1の場合で、管電流と被検体の厚さとの関係は線形であり、グラフ51はガンマ＜1の場合、グラフ52はガンマ＞1の場合で、共に管電流と被検体の厚さとの関係は非線形である。図14の場合、ガンマの値を決めることによって管電流と被検体の厚さとの関係が一義的に決まるので、［数11］のような関係式を装置に組み込んでおくことにより、操作者は操作手段21からガンマの値を入力することによって、図14の如く、管電流Iと被検体の厚さTとの関係を変化させることができる。また、実際の操作にあたっては、例えば操作者が特別な設定をしない初期設定では管電流の最大値、最小値と被検体の厚さの最大値、最小値を一致させて、線形の関係を持たせることにし、操作者の設定によりガンマを入力することで、非線形の関係を持たせることができる。図14において、ガンマ＝1を基準にした場合、ガンマ＞1のときは被検体の被曝低減を重視する場合とみられ、ガンマ＜1のときは画質を重視する場合とみられる。
【００９０】
【発明の効果】
以上説明した如く、本発明のＸ線ＣＴ装置では、装置のスキャン条件を設定する操作手段や被検体のスキャノグラム画像データから被検体の3次元的Ｘ線透過長モデルを生成するスキャノグラム解析手段や、スキャン条件と被検体の3次元的Ｘ線透過長モデルとから被検体の撮影部位に応じた管電流の変化パターンを自動的に設定する管電流設定手段や、管電流の変化パターンに基づき被検体に照射される線量を計算し、表示する線量計算手段などを備えているので、スキャン条件として管電流の最大値、最小値を入力することにより、スキャン中の管電流の変化パターンを自動的に設定することができ、かつ被検体へのＸ線被曝も評価することができる。更に、被検体へのＸ線被曝が過剰になる恐れがある場合には、管電流の変化パターンを再設定することも可能である。
【００９１】
また、本発明のＸ線ＣＴ装置では、スキャン中の管電流の変化パターンと予め生成した被検体の3次元ＣＴ値モデルデータに基づいて、被検体の体内の線量分布を計算し、表示する線量分布計算手段を備えているので、ＣＴ撮影による被検体の体内での線量分布を予め知ることができ、注目する臓器のＸ線被曝の程度を考慮して、スキャン実行の可否を判断することができる。
【００９２】
また、本発明のＸ線ＣＴ装置では、上記の被検体の3次元ＣＴ値モデルデータを人体ファントムなどをＣＴ撮影して取得した標準人体ＣＴ値モデルデータと被検体のスキャノグラム画像データに基づいて生成する被検体ＣＴ値モデル生成手段を備えているので、以前にＣＴ撮影を実施したことのない被検体についても、予備撮影としての1回のスキャノグラム画像データの取得のみによって、被検体ＣＴ値モデルデータの生成が可能である。
【００９３】
また、本発明のＸ線ＣＴ装置では、表示手段の同一画面に被検体のスキャノグラム画像と管電流の変化パターンを並置又は重ねて表示しているので、操作者などは被検体の撮影部位を見ながら管電流の変化パターンの編集を行うことが可能となり、撮影部位に適した管電流の設定を容易に行うことができる。
【図面の簡単な説明】
【図１】本発明に係るＸ線ＣＴ装置の全体構成を示すブロック図。
【図２】本発明に係るＸ線ＣＴ装置の要部構成要素。
【図３】本発明に係るＸ線ＣＴ装置を使用した第1のスキャン操作例の一連の動作のフローチャート。
【図４】スキャノグラム画像とスライス位置及び投影データ例との対応を示す図。
【図５】3次元的Ｘ線透過長モデルの1スライス位置でのモデルを示す図。
【図６】スライス位置ZにおけるＸ線透過長モデル。
【図７】管電流の変化パターンの表示例。
【図８】管電流の変化パターンの例。
【図９】本発明に係るＸ線ＣＴ装置を使用した第2のスキャン操作例の一連の動作のフローチャート。
【図１０】被検体内の線量分布の計算結果の表示例。
【図１１】被検体ＣＴ値モデルの作成手順を説明するための図。
【図１２】照射Ｘ線量分布計算の手順を説明するための図。
【図１３】被検体のスキャノグラム画像上に管電流の変化パターンを重畳して表示したもの。
【図１４】管電流と被検体の厚さとの関係を示した図。
【符号の説明】
10…ガントリ
11…ＣＴスキャナ(スキャナ)
12…Ｘ線源
13…Ｘ線検出器(検出器)
14…Ｘ線ビーム(Ｘ線)
15…被検体
16…スキャナ駆動手段
17…スキャナ角度検出手段
18…テーブル
19…コリメータ
20…ホストコンピュータ
21…操作手段
22…高電圧発生装置
23…管電流設定手段
24…画像処理装置
25…表示装置
26…スキャノグラム解析手段
29、29a…スキャノグラム画像
30…Ｘ線透過長モデル(楕円形モデル)
31、31ａ、31ｂ、31ｃ、31ｄ、47ａ、47ｂ…管電流の変化パターン
35…被検体横断面
36…被検体側面
38ａ、38ｂ、38ｃ、39ａ、39ｂ、39ｃ、39ｄ、46ａ、46ｂ、46ｃ…等線量線
41…標準人体ＣＴ値モデルデータ
42…計算スキャノグラム画像データ
43…実測スキャノグラム画像データ
44…被検体ＣＴ値モデルデータ
45…被検体μモデルデータ
51、52、53…グラフ
【発明の属する技術分野】
本発明は、撮影中のＸ線管電流（以下、管電流と略称する）を制御して被検体の被曝線量を抑制するＸ線ＣＴ装置に係り、特に被検体の被曝線量と画質とを考慮して、撮影中の管電流の変化曲線を再設定できるＸ線ＣＴ装置に関する。
【０００２】
【従来の技術】
従来のＸ線ＣＴ装置は、同一断層面では、同一のＣＴスキャン（以下、スキャンと略称する）条件（Ｘ線管電圧（以下、管電圧と略称する）や管電流等）で、撮影を行うようにしている。また、近年被検体を螺旋状にスキャンして撮影するヘリカルスキャンが広く用いられるようになってきているが、体軸方向のスキャン条件もスキャン中一定である。
【０００３】
従って、例えば被検体の断面が、ＣＴスキャナ（以下、スキャナと略称する）の回転軸に対して同心円ではなく楕円である場合、Ｘ線源の回転角度位置によって、被検体におけるＸ線の透過長が大きく変化するため、同一断層面内で、透過するＸ線量の過不足が発生する問題点を有していた。
【０００４】
また、肺等の胸部の低密度の臓器と、肝臓等の腹部の高密度の臓器とでは、Ｘ線の吸収係数が大きく異なるため、胸部から上腹部へ連続的にスキャンを行う場合、肺に適するようなＸ線量を設定すると、肝臓では不足し、肝臓に適するようなＸ線量を設定すると、肺では過剰となる事態が生じていた。
【０００５】
透過するＸ線量が不足する場合、Ｘ線検出器（以下、検出器と略称する）によって検出されるＸ線光子の量の減少によりＳ／Ｎ（ＳＮ比）が悪化し、結果として画像再構成によって得られる断層画像全体のＳ／Ｎが悪化する。逆に、透過するＸ線量が多すぎる場合には、被検体に対して、無効なＸ線被曝がなされていることになる。
【０００６】
これらの問題点を解決する方法として、特開昭53−110495号公報で開示された管電圧を制御する方法や、特開平9−108209号公報、特開平10−309271号公報で開示された管電流を制御する方法が提案されている。
【０００７】
【発明が解決しようとする課題】
しかしながら、特開昭53−110495号公報の管電圧を制御する方法は、スキャン中に管電圧を変化させるためにＸ線のスペクトルが変化し、ＣＴ値が決定できないという問題点がある。このため、現在では、管電流を制御する方法が主流となっている。
【０００８】
管電流を被検体に応じて最適に制御しようとする方法として、特開平10−309271号公報のように、スキャナ回転の半周期前の透過Ｘ線量データを用いて管電流を制御する方法と、特開平9−108209号公報のように、異なる2方向から撮影したスキャノグラムをもとにして、予め管電流を制御するパターンを被検体の位置に応じて決定しておく方法がある。
【０００９】
しかし、特開平10−309271号のスキャナ回転の半周期前の透過Ｘ線量データを用いる方法は、特にヘリカルスキャンにて、スキャンピッチを大きくする場合に、透過Ｘ線量データのずれが大きくなる問題点がある。また、横隔膜前後のように、被検体のＸ線吸収特性が大きく変わる領域では、対応することができない。
【００１０】
特開平9−108209号の異なる2方向からのスキャノグラムを取得する方法は、スキャノグラム撮影を2回行うことにより、被検体への無用なＸ線被曝を増加させることになり、管電流制御による被曝線量の低減の目的に相反するものである。
【００１１】
また、本発明の発明者達も、管電流を被検体に応じて最適に制御する一方法を特願2000−100501号にて提案している。特願2000−100501号の発明は、被検体への無用なＸ線被曝を抑えて、被検体の低被曝化を実現したＸ線ＣＴ装置に関するもので、被検体に関するスキャナの回転角度とＸ線透過長との関係を示すモデルをメモリに格納しておき、被検体に対してのスキャン計測時に、このモデルから定まるスキャナの回転角度毎の設定管電流によるＸ線を、同じ被検体に照射してスキャン計測を行い、断層画像を再構成するものである。
【００１２】
特願2000−100501号の発明では、被検体のＸ線透過長モデルの生成と、このＸ線透過長モデルのＸ線透過長に基づいて管電流を設定することに重点が置かれ、設定された管電流によって被検体に照射されるＸ線量や、被検体内の臓器へのＸ線被曝などについては考慮されていなかった。
【００１３】
上記の問題点を考慮し、本発明では、被検体のＸ線透過長モデルから自動的に設定した管電流の変化パターンによる被検体の内外における被曝線量を算定し、この被検体のＸ線被曝を考慮に入れて操作者が管電流の変化パターンを再設定することができるＸ線ＣＴ装置を提供することを目的とする。
【００１４】
【課題を解決するための手段】
上記目的を達成するため、本発明のX線CT装置は、被検体にX線を照射するX線源と、X線源に高電圧及び電流を供給する高電圧発生装置と、X線源と対向配置され被検体の透過X線量データを検出するX線検出器と、被検体の断層像を得るためのスキャン条件を設定するスキャン条件設定手段と、被検体のスキャノグラム像を表示する表示装置と、表示されたスキャノグラム像上に被検体の断層像の撮影位置を設定する撮影位置設定手段と、設定された撮影位置を含む被検体の周囲に設定されたスキャン条件によりX線源を回転しながらそのX線源からX線を照射し、X線検出に検出された被検体の透過X線データから断層像を再構成する断層像再構成手段と、再構成された断層像を表示装置に表示させる制御手段とを備えたX線CT装置において、スキャノグラム像に基づき被検体のX線透過長を計算する透過長計算手段と、計算されたX線透過長とスキャン条件とから被検体の撮影部位に応じたX線をX線源に発生させるための電流値を高電圧発生装置に設定する電流設定手段とを備えたことを特徴とする。
【００１５】
また本発明のX線CT装置は、被検体にX線を照射するX線源と、X線源に高電圧及び電流を供給する高電圧発生装置と、X線源と対向配置され被検体の透過X線量データを検出するX線検出器と、被検体の断層像を得るためのスキャン条件を設定するスキャン条件設定手段と、被検体のスキャノグラム像を表示する表示装置と、表示されたスキャノグラム像上に被検体の断層像の撮影位置を設定する撮影位置設定手段と、CTスキャナのスキャン条件を設定する条件設定手段と、設定された撮影位置を含む被検体の周囲に設定されたスキャン条件によりX線源を回転しながらそのX線源からX線を照射し、X線検出に検出された被検体の透過X線データから断層像を再構成する断層像再構成手段と、再構成された断層像を表示装置に表示させる制御手段とを備えたX線CT装置において、スキャノグラム像に基づき被検体のX線透過長を計算する透過長計算手段と、計算されたX線透過長とスキャン条件とから被検体の撮影部位に応じたX線をX線源に発生させるための電流値に応じた線量を計算する線量計算手段とを備え、制御手段は、計算した線量をスキャノグラム像と並置して、又は重ねて表示装置に表示させることを特徴とする。
【００１６】
また本発明のX線CT装置は、被検体にX線を照射するX線源と、X線源に高電圧及び電流を供給する高電圧発生装置と、X線源と対向配置され被検体の透過X線量データを検出するX線検出器と、被検体の断層像を得るためのスキャン条件を設定するスキャン条件設定手段と、被検体のスキャノグラム像を表示する表示装置と、表示されたスキャノグラム像上に被検体の断層像の撮影位置を設定する撮影位置設定手段と、CTスキャナのスキャン条件を設定する条件設定手段と、設定された撮影位置を含む被検体の周囲に設定されたスキャン条件によりX線源を回転しながらそのX線源からX線を照射し、X線検出に検出された被検体の透過X線データから断層像を再構成する断層像再構成手段と、再構成された断層像を表示装置に表示させる制御手段とを備えたX線CT装置において、スキャノグラム像に基づき被検体のX線透過長を計算する透過長計算手段と、計算されたX線透過長とスキャン条件とから被検体の撮影部位に応じたX線をX線源に発生させるための電流値に応じた被検体内の線量分布を計算する線量分布計算手段とを備え、制御手段は、計算した線量分布をスキャノグラム像と並置して、又は重ねて表示装置に表示させることを特徴とする。
【００１７】
また本発明のX線CT装置は、被検体にX線を照射するX線源に高電圧と電流を供給する高電圧発生装置と、X線源と対向配置され被検体の透過X線を検出するX線検出器と、X線源とX線検出器を被検体の周囲を回転させるスキャナを回転駆動するスキャナ駆動手段と、被検体の断層像を得るためのスキャン条件を設定する操作手段と、被検体のスキャノグラム像を得て、そのスキャノグラム像を解析し、被検体のX線透過長を求めるスキャノグラム解析手段と、スキャン条件とX線透過長とから電流値を高電圧発生装置に設定する電流設定手段と、スキャナ駆動手段にスキャナを回転させながら、設定された電流値に基づきX線源にX線を照射させて被検体の断層像を得るためのデータを収集する制御装置と、収集されたデータを再構成して断層像を再構成する画像処理装置と、再構成された断層像を表示する表示装置とを備えたことを特徴とする。
【００１８】
また本発明のX線CT装置の画像表示方法は、被検体の断層像を得るためのスキャノグラム像を得て、そのスキャノグラム像を表示装置に表示するステップと、表示されたスキャノグラム像と対応づけてX線源への第1の供給電流の変化のパターンを表示するステップと、そのスキャノグラム像を参照しながら第1の供給電流の供給電流の変化パターンから第2の供給電流の変化パターンへ操作手段によって編集するステップと、編集された第2の供給電流の変化パターンに基づいて被検体の断層像を得て、その断層像を表示装置に表示するステップとを備える。
【００１９】
また本発明のX線CT装置の画像表示方法は、X線源への供給電流を含む被検体の断層像を得るためのスキャン条件を設定するステップと、被検体のスキャノグラム像を撮影して、そのスキャノグラム像を表示装置に表示するステップと、表示されたスキャノグラム像を解析して被検体のX線透過長データを求めるステップと、X線源への供給電流と被検体のX線透過長データとから被検体の撮影部位に応じた電流パターンを計算するステップと、計算された電流パターンに基づく電流値を高電圧発生装置に設定し、その電流をX線源に供給し、その電流とスキャン条件によって被検体の断層像を得て、その断層像を表示装置に表示するステップとを備える。
【００２０】
また本発明のX線CT装置の画像表示方法は、X線源への供給電流を含む被検体の断層像を得るためのスキャン条件を設定するステップと、被検体のスキャノグラム像を撮影して、そのスキャノグラム像を表示装置に表示するステップと、表示されたスキャノグラム像を解析して被検体のX線透過長データを求めるステップと、X線源への供給電流と被検体のX線透過長データとから被検体の撮影部位に応じた電流パターンと被検体内の線量分布を計算するステップと、計算された電流パターンと被検体内の線量分布に基づく電流値を高電圧発生装置に設定し、その電流をX線源に供給し、その電流とスキャン条件によって被検体の断層像を得て、その断層像を表示装置に表示するステップとを備える。
【００２７】
【発明の実施の形態】
以下、添付図面を用いて、本発明の実施例について説明する。
図1は、本発明に係るＸ線ＣＴ装置の全体構成を示すブロック図である。図1に示すように、このＸ線ＣＴ装置は、主としてＸ線源12と検出器13等を搭載し、被検体15に対してその周囲を連続回転可能なスキャナ11を内蔵するガントリ10と、装置全体を総括するホストコンピュータ20と、Ｘ線源12に高電圧を供給する高電圧発生装置22と、画像データの前処理や画像再構成処理、あるいは各種の解析処理を行う画像処理装置24と、画像を表示する表示装置25と、被検体15を載せるテーブル装置18と、操作者がスキャン条件等を入力する操作手段21等から成る。なお、スキャナ11と被検体15とは、相対的に回転することができればよいので、被検体15が静止してスキャナ11が回転してもよいし、スキャナ11が静止して被検体15の方が回転するとしてもよい。
【００２８】
図2は、本発明に係るＸ線ＣＴ装置の要部構成要素を示す。図2を用いて、先ず、スキャナ11の詳細について説明する。図2において、スキャナ11には、Ｘ線源12と検出器13とが180度対向した位置関係で配置されている。このＸ線源12から放射されたＸ線ビーム14は、コリメータ19によってビームの幅及び厚さが制限されたファン状のＸ線ビーム14となって、被検体15に照射される。Ｘ線源12は、高電圧発生装置22を介してホストコンピュータ20によって制御される。スキャナ11全体は、スキャナ角度検出手段17によって回転角度を検出し、検出した回転角度に基づいてホストコンピュータ20が、スキャナ駆動手段16を制御し、スキャナ11を駆動する。検出器13は、被検体15を透過したＸ線14を検出し、検出データは被検体15によるＸ線の減衰量を示す投影データとして取り込まれる。投影データは、画像処理装置24において、ホストコンピュータ20の持つスキャナ角度等のデータと照合され、画像再構成等の処理をされた後に、表示装置25にて断層画像として表示される。
【００２９】
次に、本発明に係るスキャナ11以外の要部構成要素について説明する。図2において、装置全体を総括するホストコンピュータ20には、スキャナ11はスキャナ駆動手段16を介して、Ｘ線源12は高電圧発生装置22を介して、検出器13は画像処理装置24を介して、それぞれ間接的に接続されており、操作手段21と、管電流設定手段22と、画像処理装置24と、スキャノグラム解析手段26は直接接続されている。ホストコンピュータ20とスキャナ11、Ｘ線源12、検出器13との接続により、ホストコンピュータ20はＸ線源12による被検体15へのＸ線照射と、検出器13による投影データ（検出データ）の取り込みを制御する。画像処理装置24はホストコンピュータ20の指令に応じて、取り込まれた投影データに基づき、断層画像を順次再構成する。
【００３０】
本発明に係るＸ線ＣＴ装置では、被検体の断層画像を取得する本スキャンの前に、スキャン条件を設定するために種々の準備操作を行う。この準備操作としては、被検体の位置決めのためのスキャノグラム画像の撮影、管電流設定のためのスキャノグラム画像データの解析、スキャン条件としての管電流の変化パターンの決定などが、ホストコンピュータ20の介在のもとで行われる。
【００３１】
これらの準備操作に関与する主な構成要素としては、図2において、ホストコンピュータ20と、操作手段21と、スキャノグラム解析手段26と、管電流設定手段22と、Ｘ線源12と、検出器13などである。この準備操作において、先ず、操作手段21は主として管電流の設定値(最大値、最小値)などのスキャン条件をシステムに入力する。Ｘ線源12と検出器13はスキャナ11を回転させずに、スキャノグラム画像の撮影を行い、画像データをホストコンピュータ20に保存する。スキャノグラム解析手段26はスキャノグラム画像データを解析し、被検体のＸ線透過長を体軸方向のスライス位置毎及びスキャナの回転角度毎に算出可能な、3次元形状データとしてモデル化し、このモデル（以下、被検体の3次元的Ｘ線透過長モデルという）のデータをホストコンピュータ20に保存する。管電流設定手段22は操作手段21から入力された管電流設定値と被検体の3次元的Ｘ線透過長モデルのデータを基にして、スキャン中に被検体の撮影部位のＸ線透過長の変化に応じて経時的に変化する一連の管電流値すなわち、管電流の変化パターンを自動的に決定する。このように決定された管電流の変化パターンは、ホストコンピュータ20に保存され、本スキャン時に被検体の撮影部位に応じて順次呼び出されて、Ｘ線源12の管電流を変化させる。
【００３２】
本発明では、更に上記で決定された管電流の変化パターンを基に、被検体15に照射されるＸ線量を本スキャン前に予め計算する。Ｘ線装置では、通常照射したＸ線量に対応する量としてmAs値が使用されているが、この計算でもmAs値を採用している。このmAs値は、管電流（mA）と照射時間（ｓ）の積であり、管電圧が一定である場合（Ｘ線ＣＴ装置では管電圧が一定で使用される場合が多い）には、Ｘ線源12より照射されるＸ線量の総和に比例するために、Ｘ線量の基準として用いられる。ここで計算された被検体15に照射されるＸ線量については被検体15への予測被曝線量として操作者によって評価されることになる。
【００３３】
図3に、本発明に係るＸ線ＣＴ装置を使用した第1のスキャン操作の一連の動作のフローチャートを示す。このスキャン操作では、ステップ103、106、108の3次元データ生成、ステップ110の管電流パターン生成、ステップ111のmAs計算、ステップ114のmAs表示に特徴がある。以下、図2を参照しながら、図3の第1のスキャン操作の詳細について説明する。
【００３４】
図3において、先ず、ステップ101のスキャノグラム撮影の工程では、被検体15のスキャノグラム画像を撮影する。被検体15のスキャノグラム画像を撮影する構成と断層画像を撮影する構成とは基本的には同じである。本ステップでは、スキャノグラム画像データは、スキャナ11を回転させずに、被検体15に正面方向からＸ線14を照射して、検出器13によって検出データを取り込むことによって得られる。このとき得られるスキャノグラム画像は正面方向のものである。このスキャノグラム画像データは、検出器13からホストコンピュータ20に送られる。このスキャノグラム画像データは、本スキャン時の被検体15の位置決めのために利用される他、本発明では特に管電流制御のための管電流の変化パターンの決定のために利用される。
【００３５】
次に、ステップ102のスキャノグラムデータ解析の工程及びステップ103の第1の3次元データ生成の工程では、スキャノグラム画像データがホストコンピュータ20に接続されたスキャノグラム解析手段26によって解析され、被検体15の3次元的Ｘ線透過長モデルが生成される。この3次元的Ｘ線透過長モデルは、被検体15をＣＴ撮影する場合の、被検体15の位置とＸ線透過長との関係を示すモデルである。被検体15の3次元的Ｘ線透過長モデルの作成方法については特願2000−100501号にも開示されている。
【００３６】
以下、被検体15の3次元的Ｘ線透過長モデルの作成方法の一例について説明する。図4はスキャノグラム画像とスライス位置及び投影データ例との対応を示す図、図5は3次元的Ｘ線透過長モデルの1スライス位置でのモデルを示す図である。図4（a）は、ステップ101で撮影された被検体のスキャノグラム画像29を示す。このスキャノグラム画像29は胸部から腹部の中間位置までの領域を撮影領域としている。このようなスキャノグラム画像の撮影領域の中からスライス位置が選ばれる。図示の場合n個のスライス位置が選ばれている。図中、P1、…、Pi、…、Pj、…、Pnがスライス位置である。
【００３７】
図4（b），（c），（d）は3次元的Ｘ線透過長モデルのモデル決定の説明図である。任意の2つのスライス位置Pi、PjでのＣＴ断層像が図4（b）の如くなっているものと仮定すると、その縦方向（図示の上下方向）のＸ線減衰量の投影データは図4（c）の如くなるはずである。人体の体幹部の横断面は通常楕円形に近いものであるので、任意のスライス位置Pi、PjのＣＴ断層像を図4（b）で仮定するのは、大きな誤差はないものと判断される。そこで、図4（c）の投影データについては、Ｘ線透過長データに換算し、その後横軸に沿って積分して面積を求める。このとき、Ｘ線減衰量の投影データのＸ線透過長への換算にあたっては、簡単のため人体が水と等価であるとみなしてデータの変換を行う。Ｘ線減衰量をc、Ｘ線透過長をb、水の線減弱係数をμwとしたとき、両者の関係はb＝log c／μwで表される。また、横軸については、Ｘ線減衰量データの存在する領域全体の幅が人体の幅寸法に一致するように変換する。図4（d）は、図4（c）の投影データから変換したスライス位置Pi、Pjにおける被検体15のＸ線透過長データの分布図である。図4（d）においてスライス位置Pi、Pjにおける最大Ｘ線透過長はbi、bj、面積はSi、Sjとなる。図4（ｄ）のＸ線透過長データについて、最大Ｘ線透過長bi、bj及び面積Si、Sjに注目すると、bi、Siはスライス位置Piでの断層画像のＸ線透過状況を反映し、bj、Sjはスライス位置Pjでの断層画像のＸ線透過状況を反映した値と見なすことができる。
【００３８】
そこで、被検体15の3次元的Ｘ線透過長モデルとして、各スライス位置でのスライス断面を図5に示すような楕円形30でモデル化することにした。このモデル化では、スライス位置Pi、Pjにおける楕円形モデル30i、30jの面積をSi、Sj、短軸をbi、bjとしている。この結果、楕円形モデル30i、30jの長軸をai、ajとすると、楕円形モデル30i、30jの面積が［数1］で表されることから、
【数１】
【数２】
上記によって、各スライス位置での断層画像に対応するＸ線透過長モデルとしての楕円形モデル30が求められたので、これらの楕円形モデル30を体軸方向に配列することによって3次元的Ｘ線透過長モデル30を作成することができる。体軸方向のスライス位置のピッチが粗いときには、例えば隣り合う楕円形モデル間で最小自乗法により、その途中の1つ又は2つ以上の楕円形モデルを補間で求める。以上の如き手順により被検体15の3次元的Ｘ線透過長モデル30のデータとして、3次元座標（X、Y、Z）系での被検体15のＸ線透過長データT＝T（X、Y、Z）が生成される。
【００４０】
次に、ステップ104からステップ110の工程では、上記の3次元的Ｘ線透過長モデル30を使用してＸ線源12に付与する管電流の変化パターンを設定することになるが、その前に3次元的Ｘ線透過長モデル30を使用した管電流の求め方について説明する。図4、図5で求めた3次元的Ｘ線透過長モデル30は、被検体の各スライス位置での断層画像のＸ線透過長を反映したものである。3次元的Ｘ線透過長モデル30のデータは一旦ホストコンピュータ20のレジスタを含むメモリに格納されているので、撮影範囲及びテーブルピッチなどのスキャン条件が決定するとその範囲のモデルのデータがメモリから取り出され、第2、第3の3次元データの生成及び管電流の変化パターンの決定に使用される。
【００４１】
管電流は、管電流設定手段23によって各スライス位置でスキャナ回転角度毎に3次元的Ｘ線透過長モデル30から得られるＸ線透過長に基づいて決定される。図6には、スライス位置（体軸方向の位置）ZにおけるＸ線透過長モデル30を示す。また、あるスキャナ回転角度における管電流は、通常そのスキャナ回転角度における3次元的Ｘ線透過長モデルのＸ線透過長のうちの最大値に対応付けて決定される。この最大値を示すＸ線透過長は図6の楕円形モデル30の中心0を通過するパスで得られるので、管電流を設定するにあたっては、スキャナ回転角度毎にこの楕円形モデル30の中心0を通過するパスのＸ線透過長のみを考慮すればよい。従って、図6においてスライス位置をZ、スキャナ回転角度をθ（θの始点は楕円形モデル30の短軸方向とする）としたとき、その位置における最大Ｘ線透過長Tは、Zとθの関数としてT＝T（Z、θ）と表すことができる。
【００４２】
この最大Ｘ線透過長T（Z、θ）は、楕円形モデル30の中心位置0を通るパスの長さであるので、長軸をa、短軸をb、スキャナ回転角度をθとした場合、［数3］の如く表すことができる。
【数３】
【００４３】
次に、管電流の設定方法の一例について説明する。先ず、被検体をスキャンする全範囲におけるパスの最大値（全スライス位置P1〜Pnの中でのパスの最大値）をTmax、最小値（同じくパスの最小値）をTminとする。これらの値は3次元的Ｘ線透過長モデル30を作るとき既知である。管電流を最大値Imax（mA）と最小値Imin（mA）の範囲で変化させる場合、本例では管電流の最大値、最小値とパスの最大値、最小値をそれぞれ対応させて、管電流とパスとの間に直線関係を持たせるものである。管電流IとパスTとの関係は［数4］の如く表される。
【数４】
【００４４】
次に、図3のフローチャートに戻って説明する。ステップ104及びステップ105の工程では、操作者がスキャノグラム画像を参照して操作手段21からスキャン条件としてのテーブルピッチ及びスキャン開始位置を入力する。これらのデータにより、被検体のＣＴ撮影範囲とスライス位置とスキャン回転角度が決定される。このときの座標系としては、上記の如く(Z、θ)座標系がよく、スキャン条件のデータも、(Z、θ)座標系のデータで入力するのがよい。
【００４５】
次に、ステップ106の工程では、第2の3次元的Ｘ線透過長モデルのデータを生成する。この工程で生成するデータは、各スライス位置Z、スキャン回転角度θ毎の最大Ｘ線透過長であり、第1の3次元的Ｘ線透過長モデルのデータから［数3］によって求めることができるので、ホストコンピュータ20のメモリから第1の3次元的Ｘ線透過長モデルのデータを呼び出して演算する。この演算結果は、T＝T(Z、θ)で表される。
【００４６】
次に、ステップ107の工程では、スキャン条件としてのスキャン時間を操作手段21より入力する。スキャン開始位置と、テーブルピッチと、スキャン時間が決定すると、スキャン中のＸ線源12の位置(Z、θ)はスキャン開始後の経過時間tの関数として表すことができるので、各スキャン位置での被検体15の第2の3次元的Ｘ線透過長モデル、すなわち最大Ｘ線透過長Tも時間tの関数T＝T(t)として表すことができる。このため、ステップ108の第3の3次元的Ｘ線透過長モデル生成の工程では、最大Ｘ線透過長Tの関数を、T＝T(Z、θ)からT＝T(t)に変換する。
【００４７】
次に、ステップ109の管電流設定値入力の工程では、操作者が操作手段21より、管電流の設定値、例えばスキャン中の管電流の最大値Imaxと最小値Iminを入力する。ステップ110の管電流パターン生成の工程では、管電流設定手段23がホストコンピュータ20から3次元的Ｘ線透過長モデルのデータT(t)を呼び出し、上記の管電流設定値に基づき、被検体15の撮影部位に応じた管電流の変化パターンを自動的に決定する。このとき、Ｘ線透過長T(t)に対応付けてスキャン中の管電流の値を設定することになるが、Ｘ線透過長T(t)が最小の時には最小の管電流を設定し、Ｘ線透過長T(t)が最大の時には最大の管電流を設定するように管電流の変化パターンを決定する。また、Ｘ線透過長T(t)と管電流の値との関係として［数4］に示した1次関数の他、種々のものがある(但し、［数4］のTはT(t)に対応する)。
【００４８】
上記の如くして被検体15の3次元的Ｘ線透過長モデルに合わせて、管電流が時間tの関数として決定される。従って、管電流の変化パターンは、I＝I（t）と表すことができる。このように決定された管電流の変化パターンI＝I（t）はホストコンピュータ20に保存され、本スキャン時に被検体15の撮影部位に応じて順次呼び出されて、高電圧発生装置22を介してスキャン中の管電流を制御する。
【００４９】
図7に、管電流の変化パターンの表示例を示す(本例も含めて以下の管電流の変化パターンの表示例では、簡単のためスキャナ回転角度θの変化に伴う大略周期的な管電流の変化については省略し、スライス位置による変化のみ示している)。これは、表示装置25の画面上にスキャノグラム画像29と対比して表示したものである。管電流の変化パターン31では、縦軸に管電流値（mA）、横軸にスキャン開始後の経過時間tをとっている。表示例の場合、管電流はスキャンの初期（腹部）には中程度の値で、中期（腹部と胸部との間）には最小値となり、終期（胸部）には最大値となっている。本表示例の如く、管電流の変化パターン31とスキャノグラム画像29を同一画面上に並置することにより、管電流値と撮影部位との対比を一目で行うことができるので、管電流値の妥当性の判断に有効である。
【００５０】
次に、ステップ111のmAs計算の工程では、ステップ110で決定された管電流の変化パターンに基づき、スキャン中に被検体15に照射されるＸ線量を計算する。ここで、被検体15に照射されるＸ線量の基準としては上述の如く管電流（mA）と照射時間（s）との積であるmAsが用いられる。従って、本ステップでは、管電流の変化パターンI＝I（t）を時間で積分して、被検体15に照射されるＸ線量mAsを求める。この積分は［数5］によって行われる。
【数５】
【００５１】
図7に示した管電流の変化パターン31を例に上げて、被検体15への照射Ｘ線量mAsを計算することを考える。この場合、管電流の変化パターン31は、I＝I（t）で、時間の関数であるから、これを積分してmAsを計算することは、管電流の変化パターン31の図の面積を計算することになり、管電流の変化パターン31の図の面積SがmAsに相当する。
【００５２】
次に、ステップ112のmAs計算値表示の工程では、ステップ111にて求めたmAsの計算値を表示装置25の画面に表示する。ステップ111では、管電流の変化パターンが生成された段階で、被検体15のスキャン範囲の全領域にわたって被検体に照射されるＸ線量に対応するmAsが計算されるので、このステップ112でこのmAs値がスキャンを開始してもよいか否かを操作者が判断するための資料として操作者に提示される。
【００５３】
ステップ113のmAs判断の工程では、操作者が全体としてのmAsの妥当性について判断する。すなわち、操作者はこれから行われるＣＴ撮影による利益と被検体15へのＸ線被曝による不利益とを比較考量して、全体としてのmAsが大きすぎないか否かを判断し、全体としてのmAsが大きすぎると判断した場合には管電流の設定値を下げることになる。この場合には、ステップ109に戻り、管電流設定値を再度入力し、管電流の変化パターンを再設定する。
【００５４】
図8には、管電流の変化パターンの例を示す。図8（a）に示した管電流の変化パターン31aは、通常の管電流一定（I0）の場合で、管電流の最小値（Imin）と最大値（Imax）が同じI0である場合の例である。図8（a）の管電流の変化パターン31aのグラフでは、mAsは面積Saである。次に、図8（b）示した管電流の変化パターン31bでは、管電流は初期には図8（a）と同じI0で、中期にI0より低い最小値Iminとなり、終期にI0より高い最大値Imaxとなっている。図8（b）の管電流の変化パターン31aのグラフでのmAsは面積Sbであるが、この面積Sbは図8（a）の面積Saより小さくなっており、被検体15の被曝線量は低減している。
【００５５】
図8（c）に示した管電流の変化パターン31cでは、図8（b）のグラフと比べて、終期の最大値Imaxを小さくし、I0とほぼ同じとしたものである。この場合、スキャン全体として管電流が低目に抑えられている。図8（c）の管電流の変化パターン31cのグラフでのmAsは面積Scであるが、この面積Scは図8（b）の面積Sbより更に小さくなっており、被検体15の被曝線量は更に低減している。
【００５６】
図8（d）に示した管電流の変化パターン31dでは、図8（b）に比べて中期の最小値Iminを更に小さくし、IminがI0に対し大幅に低下するようにしたものである。図8（d）の管電流の変化パターン31dでのmAsは面積Sdであるが、図8（c）の面積Scより更に小さくなっている。
【００５７】
図8（c）の場合の如く管電流設定値の最大値を下げた場合には、被検体15の厚さが厚い部分、すなわちＸ線透過長が大きい部分での管電流を減らすことになる。従って、図8（c）の管電流の変化パターン31cは、特に肺のように低密度の部位での画質を重視し、腹部のように高密度の領域での被曝線量を低減したい場合に適する。
【００５８】
図8（d）の場合の如く管電流設定値の最小値を下げた場合には、被検体15の厚さが薄い部分、すなわちＸ線透過長が小さい部分での管電流を減らすことになる。従って、図8（d）の管電流の変化パターン31dは、特に骨周辺や実質部のように高密度の部位での画質を重視し、低密度の領域での被曝線量を低減したい場合に適する。
【００５９】
図8（c），図8（d）では、管電流設定値の最大値Imax又は最小値Iminのいずれか一方のみを小さくしているが、両方を小さくして、被検体15の全領域について平均的に被曝線量を低減することも可能である。上記の如く新たに設定した管電流の変化パターンについては、再度mAsを計算し、問題が無い場合には、そのまま採用する。
【００６０】
次に、ステップ114のスキャンの工程では、操作者は上記で決定した管電流の変化パターンを含めたスキャン条件にてスキャンを実行する。
【００６１】
次に、ステップ115のmAs積算値表示の工程では、上記の［数５］に従って、スキャン中のmAsの積算値を逐次計算し、表示装置25の画面にリアルタイムで表示し、操作者に提示する。mAs積算値表示方法としては、全体としてのmAs値に対する比率としての相対値を表示する方法と、mAs積算値の絶対値を表示する方法のいずれかを選択することができる。もちろん、その両方を同時に表示することも可能である。
【００６２】
次に、本発明に係るＸ線ＣＴ装置の第2のスキャン操作例について説明する。図9は、第2のスキャン操作例の一連の動作のフローチャートを示したものである。このスキャン操作例では、図3の第1のスキャン操作例のフローチャートに対し、ステップ111のmAs計算及びステップ112のmAs表示の工程の後に、被検体15内の線量分布の計算及び表示の工程を追加し、操作者が被検体15内の線量分布をも見てスキャン実行の判断を下せるようにしたものである。このため、本操作例の説明では、ステップ201〜205の被検体15内の線量分布の計算及び表示の工程を重点に説明する。
【００６３】
以下、図9の第2のスキャン操作例の中のステップ201からステップ205について説明するが、ステップの内容説明に入る前に、被検体15内の線量分布の計算結果の表示例を図10に示す。図10（a）は被検体の横断面35での線量分布の表示例、図10（b）は被検体の体軸方向の側面36での線量分布の表示例である。両図とも被検体15内の分布線量の等しい等線量線38a〜38c、39a〜39cが示されており、体表に近いほど高線量になっている。本操作例では、被検体15の線量分布が操作者に提示されるため、操作者は被検体のＸ線被曝についてより詳細な評価を行うことができるということが大きな特徴となる。
【００６４】
本操作例では、ステップ203で被検体内の線量分布の計算を行う前に、準備として、ステップ201で被検体のＣＴ値モデルの生成、ステップ202で被検体のμモデルの生成を行い、その後で被検体のμモデルデータと被検体に照射されたＸ線量データを基にして被検体内の線量分布を計算する。
【００６５】
先ず、ステップ201の被検体ＣＴ値モデル生成の工程では、予め標準的な人体のＣＴ値分布モデル（以下、標準人体ＣＴ値モデルという）データを取得して、ホストコンピュータ20の記憶手段に保存しておき、この標準人体ＣＴ値モデルデータを、ステップ101で取得した被検体15のスキャノグラム画像のデータに基づいて補正を行うことによって、被検体15のＣＴ値モデル（以下、被検体ＣＴ値モデルという）を作成する。上記の標準人体ＣＴ値モデルデータとしては、例えば標準的な人体ファントムなどをＣＴ撮影した断層画像から得られる3次元的ＣＴ値分布データが用いられる。断層画像はＣＴ値の分布を表し、実効エネルギー（通常60keV）のＸ線に対する線減弱係数の分布を表しているので、この断層画像を3次元に再構成をした3次元的ＣＴ値分布データは線減弱係数の3次元的な空間分布のデータであり、被検体に照射されたＸ線の減弱量の計算に利用することができる。
【００６６】
次に、被検体ＣＴ値モデルの生成方法の一例について、図11を用いて説明する。この例では、ステップ101で被検体15を撮影して取得して実測したスキャノグラム画像データと、上記の標準人体ＣＴ値モデルデータとから、被検体15の3次元的ＣＴ値分布を表す被検体ＣＴ値モデルデータを生成する。この被検体ＣＴ値モデルの生成にあたっては、被検体15のスキャノグラム画像データと標準的人体のスキャノグラム画像データが媒体として利用される。
【００６７】
図11は、ステップ201の被検体ＣＴ値モデルの作成手順を説明するための図である。図11において、図11(a)は標準人体ＣＴ値モデルデータ41の例を、図11(b)は図11(a)の標準人体ＣＴ値モデルデータ41から計算で求めた標準人体のスキャノグラム画像データ42の例を、図11(c)は被検体の実測のスキャノグラム画像データ43の例を、図11(d)は計算で求めた被検体ＣＴ値モデルデータ44の例を示す。図11(a)の標準人体ＣＴ値モデルデータ41は、人体ファントムなどの標準的な人体の体幹部のＣＴ値分布モデルで、肩から腹部までのスライス位置ごとの断面のＣＴ値分布モデルを示している。
【００６８】
スキャノグラム画像は、その3次元的ＣＴ値分布モデルから計算によって生成することができるので、図11(b)の標準人体スキャノグラム画像データ42は図11(a)の標準人体ＣＴ値モデルデータ41について正面方向から投影したデータを求めることによって得られる。図11(c)の被検体の実測スキャノグラム画像データ43は、被検体15の体幹部について標準人体スキャノグラム画像データ42と同じ領域を正面方向から撮影したスキャノグラム画像データである。この画像については、以下被検体スキャノグラム画像と呼ぶことにする。
【００６９】
図11においては、体幹部の標準人体スキャノグラム画像データ42と被検体スキャノグラム画像データ43とを並置して対比できるように示してあるが、両者は普通寸法及びＣＴ値とも異なるものである。このため、標準人体スキャノグラム画像データ42と被検体スキャノグラム画像データ43とを対比しながら、両者の差異に基づいて、一致する部分はそのままとし、異なる部分については変形させて、被検体15に合うように標準人体ＣＴ値モデルデータ41を補正して、被検体ＣＴ値モデルデータ44を生成する。
【００７０】
図11の体幹部の例では、先ず体軸方向に関して、標準人体スキャノグラム画像データ42と被検体スキャノグラム画像データ43の肩から横隔膜までの長さAと、横隔膜から腸管までの長さBとに分けて、それぞれの差異を基に、標準人体ＣＴ値モデルデータ41を補間、伸長したり、あるいは間引き、短縮したりすることで、標準人体ＣＴ値モデルデータ41の体軸方向のＣＴ値分布を被検体15の実状に近似させる。左右方向に関しても、同様に、体軸を基準にして左と右に分けて、それぞれの差異を基に、左右の広がりを補正し、被検体15の実状に近似させる。前後方向に関しては、被検体スキャノグラム画像データ43から推定される前後方向のＸ線透過長を基に、標準人体ＣＴ値モデルデータ41の前後方向のデータを線形に補間する。このようにして、2つのスキャノグラム画像データ42、43を基に、標準人体ＣＴ値モデルデータ41を実際の被検体15に合わせ込むことにより、被検体ＣＴ値モデルデータ44を生成する。
【００７１】
次に、ステップ202の被検体μモデル生成の工程では、ステップ201で生成した被検体ＣＴ値モデルデータ44のＣＴ値を線減弱係数μに変換して、被検体15の3次元的μ値分布モデルを生成する。ＣＴ値から線減弱係数μへの変換は下記の如く行われる。
【００７２】
ＣＴ値は実効エネルギー(通常60keVを使用)のＸ線に対する線減弱係数によって決定され、水＝0、空気＝1000、平均的な骨＝1000として定義されている。今、被検体ＣＴ値モデルの位置XにおけるＣＴ値をＣＴXとすると、その位置Ｘにおける実効エネルギー(60keV)での線減弱係数μXは［数6］、［数7］で表される。
【数６】
【００７３】
次に、ステップ203の線量分布計算の工程では、図12に示す如く、ステップ111で求めた被検体に照射されたＸ線量のデータ(図12(a))とステップ202で求めた被検体μモデルデータ45(図12(b))を用いて、被検体15内の線量分布(図12(c))を計算する。このステップ203の計算では、被検体に照射されるＸ線のエネルーギースペクトルを考慮して、被検体15内でのＸ線の減衰を計算し、被検体15内の線量の空間的分布を計算する。被検体15に任意の方向からＸ線を照射した場合のＸ線の減弱量は、被検体の3次元的線減弱係数(μ)のモデルである被検体μモデルデータを用いることにより解析的に計算可能であり、このような計算手法は既に他の分野、例えば放射線治療計画装置などの分野でも行われている(参考文献1、稲邑清也、放射線治療計画システム、P．90〜92、P．113〜115、篠原出版、平成4月20日発行)。
【００７４】
Ｘ線の減衰の計算にあたっては、先ず、被検体15内の注目する位置XからＸ線源12に向かってＸ線が透過する実効距離δを計算する。実効距離は、Ｘ線が透過する媒質によって1／ｅに減衰する距離を1と定義される。Ｘ線のエネルギースペクトルを考慮しない場合、注目位置ＸからＸ線源12までに、組成ｉの媒質が実距離でdi(cm)あり、かつ組成ｉの線減弱係数がμi(cm-1)であるとき、実効距離δは［数8］で表される。
【数８】
しかし、Ｘ線のエネルギースペクトルを考慮すると実効距離はＸ線のエネルギーに応じて異なる値となる。組成iのＸ線のエネルギーjに対する線減弱係数をμij(cm-1)とすると、Ｘ線のエネルギーjに対する実効距離δjは［数9］で表される。
【数９】
【００７６】
次に、被検体15内の注目位置Xでの線量を計算する。Ｘ線源12から注目位置Xまでの距離をrx(m)、距離1mにおける線量をI0(C／kg：Cはクーロン)とする。I0については、例えば実験的に求める。Ｘ線のエネルギースペクトル、つまりエネルギーjの成分比をSjとすると、注目位置Ｘでの線量Ix(C／kg)は［数10］で表される。
【数１０】
【００７７】
［数10］による計算の結果、Ｘ線源12がある位置Q(Z、θ)にある場合の被検体15内の任意に位置Ｘにおける線量が求められる。Ｘ線源12を1回転した場合の任意の位置Ｘにおける線量は、Ｘ線源12の位置Q(Z、θ)を被検体15の周囲で回転させ、上記の線量を1回転分(θ=0〜2π)積算することによって得られる。上記の如き手順で、被検体15内の設定された各位置での計算を行うことにより、スライス位置Zでの被検体15内の線量分布を求めることができる。また、他のスライス位置における線量分布も同様に計算によって求めることができるので、被検体15の体軸方向について撮影部位の全領域の計算を進めることにより、被検体15内の3次元的な線量分布が得られる。
【００７８】
1スライスの断層画像を撮影する場合(2次元の場合)には上記の計算によって被検体15内の線量の分布が精度良く求められるが、1回のスキャン中に複数スライスのＣＴ撮影を行う場合(3次元の場合)には散乱Ｘ線の撮影を考慮しないと計算精度が低下する恐れがある。Ｘ線ＣＴ装置では、Ｘ線のエネルギーが100keV以下のオーダーであるため、散乱線としてはコンプトン(Compton)散乱のみ考慮すればよい(参考文献１参照)。このコンプトン散乱を考慮することにより、計算精度をより高くすることができる。
【００７９】
上記の如き手順で、被検体15内の各スライス位置での各設定点における線量を計算することにより、計算された被検体15の3次元的線量分布が得られる。この被検体15の線量分布の計算結果はホストコンピュータ20に一時的に保存され、操作者などにとって見やすい図象、例えば図10に示したような図象で表示される。
【００８０】
次に、ステップ204の線量分布表示の工程ではステップ203の計算結果が表示装置25に表示される。本実施例での表示例としては、図10に示したような被検体15の断層面35内の線量分布(図10(a))又は被検体15の側面36の線量分布(図10(b))が上げられる。これらの図では、被検体15の臓器と線量分布を示す等線量線を重ねて表示しているので、一目で各臓器への被曝線量を認識することができるので、被検体15へのＸ線被曝を評価する上では有効である。
【００８１】
次に、ステップ205の被曝線量判断の工程では、ステップ204で表示された被検体15内の線量分布の計算結果を操作者が見て、被検体15内の臓器へのＸ線被曝が過剰になる恐れがないかどうかを判断し、Yesと判断した場合にはステップ114のスキャン実行の工程に進み、スキャンを開始することになり、Noと判断した場合にはステップ109の管電流設定値入力の工程に戻り、管電流設定値の再入力、管電流パターンの再検討を行うことになる。
【００８２】
また、ステップ201での被検体15のＣＴ値モデルデータの生成方法としては、上記の標準人体ＣＴ値モデルデータを使う方法以外に、過去に撮影した同じ被検体15のＣＴ撮影データを使用する方法も実施可能である。この場合には、実際に同一の被検体15のＣＴ値分布データを使用するために、標準人体ＣＴ値モデルデータ41の形状を補正する手順が必要なくなるという利点がある。しかし初回のＣＴ撮影には適さないため、過去にＣＴ撮影を行った被検体に関して、2回目以降のＣＴ撮影を行う場合が対象となる。
【００８３】
上記したように、ＣＴスキャン前に被検体15の体内の線量分布をシミュレーション計算して、計算結果を、例えば図10に示す如く表示することにより、操作者は事前に撮影手技に応じた被検体の体内の線量分布を近似的に知ることが可能となる。
【００８４】
この結果、例えば、単純に被検体の全ての組織に関して一律に被曝線量を減らすのではなく、骨髄や肺等の放射線感受性の高い組織に関しては特に被曝線量を低減し、逆に脂肪や筋肉等の放射線感受性の比較的低い組織には、画質が満足できる程度に被曝線量のレベルを維持するというような詳細な設定が可能となる。
【００８５】
次に、図13を用いて、被検体をＣＴ撮影する管電流の変化パターンの編集例を説明する。図13は、被検体のスキャノグラム画像上に管電流の変化パターンを重畳して表示したものである。図12において、スキャノグラム画像データ29aは体幹部のもの、管電流の変化パターンは、編集前の初期の管電流の変化パターン46aと編集後の修正された管電流の変化パターン46bである。
【００８６】
この管電流の変化パターンの編集工程では、表示装置25の画面においてスキャノグラム画像29a上に表示された初期に設定された管電流の変化パターン46aに対し、スキャノグラム画像データ29aを参照しながら、また場合によっては被検体15の内部の照射線量分布を参照して、操作手段21によって修正を加えて、新しい管電流の変化パターン46bを編集する。この編集操作によって任意の部位の管電流の変化パターンを再設定する。
【００８７】
この編集操作において、例えば、自動的な管電流の変化パターンの設定では、横隔膜付近のように密度が大きく変化する領域では、管電流を平均的な値に設定するが、被曝線量が増えても画質を向上させる必要がある領域などでは、管電流を部分的に高く設定する。管電流の変化パターンは、上記の如くスキャン条件が設定されていれば、時間tのみの関数になるので、任意時刻の管電流の値を変化させることができる。図13の例では、初期の管電流の変化パターン46aに対し、肺の領域の管電流を少し低下させ、横隔膜の領域の管電流を少し増加させることで、修正後の管電流の変化パターン46bに編集している。
【００８８】
図14は、管電流と被検体の厚さ(Ｘ線透過長に相当)との関係を示した図である。上記のステップ110の工程の説明では、管電流の最大値、最小値と被検体の厚さの最大値、最小値とを一致させて、両者間で線形の関係を持つものとして説明したが、両者の関係については操作者の設定により非線形の関係を持つようにすることが可能である。図14に示した管電流Iと被検体の厚さTとの関係は、［数11］で表される。
【数１１】
【００８９】
図14において、グラフ50はガンマ=1の場合で、管電流と被検体の厚さとの関係は線形であり、グラフ51はガンマ＜1の場合、グラフ52はガンマ＞1の場合で、共に管電流と被検体の厚さとの関係は非線形である。図14の場合、ガンマの値を決めることによって管電流と被検体の厚さとの関係が一義的に決まるので、［数11］のような関係式を装置に組み込んでおくことにより、操作者は操作手段21からガンマの値を入力することによって、図14の如く、管電流Iと被検体の厚さTとの関係を変化させることができる。また、実際の操作にあたっては、例えば操作者が特別な設定をしない初期設定では管電流の最大値、最小値と被検体の厚さの最大値、最小値を一致させて、線形の関係を持たせることにし、操作者の設定によりガンマを入力することで、非線形の関係を持たせることができる。図14において、ガンマ＝1を基準にした場合、ガンマ＞1のときは被検体の被曝低減を重視する場合とみられ、ガンマ＜1のときは画質を重視する場合とみられる。
【００９０】
【発明の効果】
以上説明した如く、本発明のＸ線ＣＴ装置では、装置のスキャン条件を設定する操作手段や被検体のスキャノグラム画像データから被検体の3次元的Ｘ線透過長モデルを生成するスキャノグラム解析手段や、スキャン条件と被検体の3次元的Ｘ線透過長モデルとから被検体の撮影部位に応じた管電流の変化パターンを自動的に設定する管電流設定手段や、管電流の変化パターンに基づき被検体に照射される線量を計算し、表示する線量計算手段などを備えているので、スキャン条件として管電流の最大値、最小値を入力することにより、スキャン中の管電流の変化パターンを自動的に設定することができ、かつ被検体へのＸ線被曝も評価することができる。更に、被検体へのＸ線被曝が過剰になる恐れがある場合には、管電流の変化パターンを再設定することも可能である。
【００９１】
また、本発明のＸ線ＣＴ装置では、スキャン中の管電流の変化パターンと予め生成した被検体の3次元ＣＴ値モデルデータに基づいて、被検体の体内の線量分布を計算し、表示する線量分布計算手段を備えているので、ＣＴ撮影による被検体の体内での線量分布を予め知ることができ、注目する臓器のＸ線被曝の程度を考慮して、スキャン実行の可否を判断することができる。
【００９２】
また、本発明のＸ線ＣＴ装置では、上記の被検体の3次元ＣＴ値モデルデータを人体ファントムなどをＣＴ撮影して取得した標準人体ＣＴ値モデルデータと被検体のスキャノグラム画像データに基づいて生成する被検体ＣＴ値モデル生成手段を備えているので、以前にＣＴ撮影を実施したことのない被検体についても、予備撮影としての1回のスキャノグラム画像データの取得のみによって、被検体ＣＴ値モデルデータの生成が可能である。
【００９３】
また、本発明のＸ線ＣＴ装置では、表示手段の同一画面に被検体のスキャノグラム画像と管電流の変化パターンを並置又は重ねて表示しているので、操作者などは被検体の撮影部位を見ながら管電流の変化パターンの編集を行うことが可能となり、撮影部位に適した管電流の設定を容易に行うことができる。
【図面の簡単な説明】
【図１】本発明に係るＸ線ＣＴ装置の全体構成を示すブロック図。
【図２】本発明に係るＸ線ＣＴ装置の要部構成要素。
【図３】本発明に係るＸ線ＣＴ装置を使用した第1のスキャン操作例の一連の動作のフローチャート。
【図４】スキャノグラム画像とスライス位置及び投影データ例との対応を示す図。
【図５】3次元的Ｘ線透過長モデルの1スライス位置でのモデルを示す図。
【図６】スライス位置ZにおけるＸ線透過長モデル。
【図７】管電流の変化パターンの表示例。
【図８】管電流の変化パターンの例。
【図９】本発明に係るＸ線ＣＴ装置を使用した第2のスキャン操作例の一連の動作のフローチャート。
【図１０】被検体内の線量分布の計算結果の表示例。
【図１１】被検体ＣＴ値モデルの作成手順を説明するための図。
【図１２】照射Ｘ線量分布計算の手順を説明するための図。
【図１３】被検体のスキャノグラム画像上に管電流の変化パターンを重畳して表示したもの。
【図１４】管電流と被検体の厚さとの関係を示した図。
【符号の説明】
10…ガントリ
11…ＣＴスキャナ(スキャナ)
12…Ｘ線源
13…Ｘ線検出器(検出器)
14…Ｘ線ビーム(Ｘ線)
15…被検体
16…スキャナ駆動手段
17…スキャナ角度検出手段
18…テーブル
19…コリメータ
20…ホストコンピュータ
21…操作手段
22…高電圧発生装置
23…管電流設定手段
24…画像処理装置
25…表示装置
26…スキャノグラム解析手段
29、29a…スキャノグラム画像
30…Ｘ線透過長モデル(楕円形モデル)
31、31ａ、31ｂ、31ｃ、31ｄ、47ａ、47ｂ…管電流の変化パターン
35…被検体横断面
36…被検体側面
38ａ、38ｂ、38ｃ、39ａ、39ｂ、39ｃ、39ｄ、46ａ、46ｂ、46ｃ…等線量線
41…標準人体ＣＴ値モデルデータ
42…計算スキャノグラム画像データ
43…実測スキャノグラム画像データ
44…被検体ＣＴ値モデルデータ
45…被検体μモデルデータ
51、52、53…グラフ
Claims (8)
- 被検体にＸ線を照射するＸ線源と、Ｘ線源に高電圧及び電流を供給する高電圧発生装置と、Ｘ線源と対向配置され被検体の透過Ｘ線量データを検出するＸ線検出器と、被検体の断層像を得るためのスキャン条件を設定するスキャン条件設定手段と、被検体のスキャノグラム像を表示する表示装置と、表示されたスキャノグラム像上に被検体の断層像の撮影位置を設定する撮影位置設定手段と、設定された撮影位置を含む被検体の周囲に設定されたスキャン条件によりＸ線源を回転しながらそのＸ線源からＸ線を照射し、Ｘ線検出器に検出された被検体の透過Ｘ線データから断層像を再構成する断層像再構成手段と、再構成された断層像を表示装置に表示させる制御手段とを備えたＸ線ＣＴ装置において、
スキャノグラム像上に設定された撮影位置に対応する投影データを変換したＸ線透過長の分布データの最大値と面積から楕円形モデルを作成する楕円形モデル作成手段と、前記楕円モデルとスキャン条件とから被検体の撮影部位に応じたＸ線をＸ線源に発生させるための電流値を高電圧発生装置に設定する電流設定手段とを備えたことを特徴とするＸ線ＣＴ装置。 - 上記電流設定手段は、Ｘ線透過長とＸ線源に供給する電流値の特性曲線のガンマを用いて、被検体の部位に応じた電流値を高電圧発生装置に設定することを特徴とする請求項１に記載のＸ線ＣＴ装置。
- 被検体にＸ線を照射するＸ線源と、Ｘ線源に高電圧及び電流を供給する高電圧発生装置と、Ｘ線源と対向配置され被検体の透過Ｘ線量データを検出するＸ線検出器と、被検体の断層像を得るためのスキャン条件を設定するスキャン条件設定手段と、被検体のスキャノグラム像を表示する表示装置と、表示されたスキャノグラム像上に被検体の断層像の撮影位置を設定する撮影位置設定手段と、ＣＴスキャナのスキャン条件を設定する条件設定手段と、設定された撮影位置を含む被検体の周囲に設定されたスキャン条件によりＸ線源を回転しながらそのＸ線源からＸ線を照射し、Ｘ線検出器に検出された被検体の透過Ｘ線データから断層像を再構成する断層像再構成手段と、再構成された断層像を表示装置に表示させる制御手段とを備えたＸ線ＣＴ装置において、
スキャノグラム像上に設定された撮影位置に対応する投影データを変換したＸ線透過長の分布データの最大値と面積から楕円形モデルを作成する楕円形モデル作成手段と、前記楕円モデルとスキャン条件とから被検体の撮影部位に応じたＸ線をＸ線源に発生させるための電流値に応じた線量を計算する線量計算手段とを備え、制御手段は、計算した線量をスキャノグラム像と並置して、又は重ねて表示装置に表示させることを特徴とするＸ線ＣＴ装置。 - 被検体にＸ線を照射するＸ線源と、Ｘ線源に高電圧及び電流を供給する高電圧発生装置と、Ｘ線源と対向配置され被検体の透過Ｘ線量データを検出するＸ線検出器と、被検体の断層像を得るためのスキャン条件を設定するスキャン条件設定手段と、被検体のスキャノグラム像を表示する表示装置と、表示されたスキャノグラム像上に被検体の断層像の撮影位置を設定する撮影位置設定手段と、ＣＴスキャナのスキャン条件を設定する条件設定手段と、設定された撮影位置を含む被検体の周囲に設定されたスキャン条件によりＸ線源を回転しながらそのＸ線源からＸ線を照射し、Ｘ線検出器に検出された被検体の透過Ｘ線データから断層像を再構成する断層像再構成手段と、再構成された断層像を表示装置に表示させる制御手段とを備えたＸ線ＣＴ装置において、
スキャノグラム像上に設定された撮影位置に対応する投影データを変換したＸ線透過長の分布データの最大値と面積から楕円形モデルを作成する楕円形モデル作成手段と、前記楕円モデルとスキャン条件とから被検体の撮影部位に応じたＸ線をＸ線源に発生させるための電流値に応じた被検体内の線量分布を計算する線量分布計算手段とを備え、制御手段は、計算した線量分布をスキャノグラム像と並置して、又は重ねて表示装置に表示させることを特徴とするＸ線ＣＴ装置。 - 被検体にＸ線を照射するＸ線源に高電圧と電流を供給する高電圧発生装置と、Ｘ線源と対向配置され被検体の透過Ｘ線を検出するＸ線検出器と、Ｘ線源とＸ線検出器を被検体の周囲を回転させるスキャナを回転駆動するスキャナ駆動手段と、被検体の断層像を得るためのスキャン条件を設定する操作手段と、被検体のスキャノグラム像を得て、そのスキャノグラム像上に設定された撮影位置に対応する投影データを変換したＸ線透過長の分布データの最大値と面積から楕円形モデルを作成する楕円形モデル作成手段と、前記楕円モデルとスキャン条件とから電流値を高電圧発生装置に設定する電流設定手段と、スキャナ駆動手段にスキャナを回転させながら、設定された電流値に基づきＸ線源にＸ線を照射させて被検体の断層像を得るためのデータを収集する制御装置と、収集されたデータを再構成して断層像を再構成する画像処理装置と、再構成された断層像を表示する表示装置とを備えたことを特徴とするＸ線ＣＴ装置。
- Ｘ線源への供給電流を含むスキャン条件の設定を受付けるステップと、被検体のスキャノグラム像を撮影して、そのスキャノグラム像を表示装置に表示するステップと、表示されたスキャノグラム像上に設定された撮影位置に対応する投影データを変換したＸ線透過長の分布データの最大値と面積から楕円形モデルを作成するステップと、Ｘ線源への供給電流と前記楕円モデルとから被検体の撮影部位に応じた電流パターンを計算するステップと、計算された電流パターンに基づく電流値を高電圧発生装置に設定し、その電流をＸ線源に供給し、その電流とスキャン条件によってX線源がX線を放射し、Ｘ線検出器により検出された被検体の透過Ｘ線データから断層像が再構成されて、その断層像を表示装置に表示するステップとを備えたＸ線ＣＴ装置の画像表示方法。
- Ｘ線源への供給電流を含むスキャン条件の設定を受付けるステップと、被検体のスキャノグラム像を撮影して、そのスキャノグラム像を表示装置に表示するステップと、表示されたスキャノグラム像上に設定された撮影位置に対応する投影データを変換したＸ線透過長の分布データの最大値と面積から楕円形モデルを作成するステップと、Ｘ線源への供給電流と前記楕円モデルとから被検体の撮影部位に応じた電流パターンと被検体内の線量分布を計算するステップと、計算された電流パターンと被検体内の線量分布に基づく電流値を高電圧発生装置に設定し、その電流をＸ線源に供給し、その電流とスキャン条件によってX線源がX線を放射し、Ｘ線検出器により検出された被検体の透過Ｘ線データから断層像が再構成されて、その断層像を表示装置に表示するステップとを備えたＸ線ＣＴ装置の画像表示方法。
- 上記表示ステップは、被検体の部位に応じて設定される電流値の変化パターンと被検体内の線量分布と被検体のスキャノグラム画像とを並置して、又は重ねて表示することを特徴とする請求項６に記載のＸ線ＣＴ装置の画像表示方法
Priority Applications (5)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001067150A JP4532005B2 (ja) | 2001-03-09 | 2001-03-09 | Ｘ線ｃｔ装置及びその画像表示方法 |
| EP08018586A EP2047799B1 (en) | 2001-03-09 | 2002-03-08 | X-ray CT image display method |
| EP02702836A EP1393681B1 (en) | 2001-03-09 | 2002-03-08 | X-ray ct device and image displaying method therefor |
| US10/468,711 US7103139B2 (en) | 2001-03-09 | 2002-03-08 | X-ray CT device and image displaying method therefor |
| PCT/JP2002/002190 WO2002071946A1 (fr) | 2001-03-09 | 2002-03-08 | Tomodensitometre et son procede d'affichage |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001067150A JP4532005B2 (ja) | 2001-03-09 | 2001-03-09 | Ｘ線ｃｔ装置及びその画像表示方法 |
Related Child Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2009099614A Division JP4644292B2 (ja) | 2009-04-16 | 2009-04-16 | Ｘ線ｃｔ装置とその画像表示方法 |
Publications (3)
| Publication Number | Publication Date |
|---|---|
| JP2002263097A JP2002263097A (ja) | 2002-09-17 |
| JP2002263097A5 JP2002263097A5 (ja) | 2008-04-24 |
| JP4532005B2 true JP4532005B2 (ja) | 2010-08-25 |
Family
ID=18925543
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2001067150A Expired - Lifetime JP4532005B2 (ja) | 2001-03-09 | 2001-03-09 | Ｘ線ｃｔ装置及びその画像表示方法 |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US7103139B2 (ja) |
| EP (2) | EP1393681B1 (ja) |
| JP (1) | JP4532005B2 (ja) |
| WO (1) | WO2002071946A1 (ja) |
Families Citing this family (86)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4154990B2 (ja) * | 2002-10-17 | 2008-09-24 | 株式会社島津製作所 | Ｘ線ｃｔ装置 |
| JP2004173924A (ja) * | 2002-11-27 | 2004-06-24 | Ge Medical Systems Global Technology Co Llc | Ｘ線制御方法およびｘ線画像撮影装置 |
| JP4316231B2 (ja) * | 2002-12-26 | 2009-08-19 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | Ｘ線ｃｔシステムならびに画像処理装置および方法 |
| WO2004071301A1 (ja) * | 2003-02-14 | 2004-08-26 | Hitachi Medical Corporation | X線ct装置 |
| JP4087274B2 (ja) * | 2003-04-01 | 2008-05-21 | 株式会社東芝 | Ｘ線厚さ測定装置 |
| JP4490645B2 (ja) * | 2003-04-09 | 2010-06-30 | 株式会社東芝 | Ｘ線コンピュータ断層撮影装置 |
| JP3999176B2 (ja) * | 2003-08-20 | 2007-10-31 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | Ｘ線ｃｔ装置、情報処理方法ならびに記憶媒体、プログラム |
| JP2005074000A (ja) * | 2003-09-01 | 2005-03-24 | Ge Medical Systems Global Technology Co Llc | Ｘ線ｃｔ装置および被曝線量分布算出方法 |
| JP4429694B2 (ja) * | 2003-11-13 | 2010-03-10 | 株式会社日立メディコ | Ｘ線ｃｔ装置 |
| JP4554185B2 (ja) * | 2003-11-18 | 2010-09-29 | 株式会社日立メディコ | Ｘ線ｃｔ装置 |
| JP2005168870A (ja) * | 2003-12-12 | 2005-06-30 | Shimadzu Corp | Ｘ線ｃｔ装置 |
| JP2005185718A (ja) * | 2003-12-26 | 2005-07-14 | Ge Medical Systems Global Technology Co Llc | 放射線断層撮像装置および撮像方法 |
| US6990172B2 (en) * | 2004-02-19 | 2006-01-24 | General Electric Company | Method and apparatus to determine tube current modulation profile for radiographic imaging |
| WO2005099577A1 (en) * | 2004-04-13 | 2005-10-27 | Koninklijke Philips Electronics, N.V. | Dynamic dose control for computed tomography |
| US7082183B2 (en) * | 2004-07-21 | 2006-07-25 | General Electric Company | Computed tomography dose indexing phantom selection for dose reporting |
| JP4731151B2 (ja) * | 2004-10-22 | 2011-07-20 | 株式会社日立メディコ | Ｘ線管電流決定方法及びｘ線ｃｔ装置 |
| JP4739738B2 (ja) * | 2004-12-01 | 2011-08-03 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 線量評価方法およびｘ線ｃｔ装置 |
| DE102005005285B4 (de) * | 2005-02-04 | 2014-02-13 | Siemens Aktiengesellschaft | Verfahren für einen ein Aufnahmesystem aufweisenden Computertomographen zur Abtastung eines Objektes |
| CN101128153B (zh) * | 2005-02-25 | 2010-09-29 | 株式会社日立医药 | X射线ct设备 |
| JP4954062B2 (ja) * | 2005-04-04 | 2012-06-13 | 株式会社日立メディコ | Ｘ線ｃｔ装置 |
| JP2006320523A (ja) * | 2005-05-19 | 2006-11-30 | Ge Medical Systems Global Technology Co Llc | シャトルモードヘリカルスキャンのスキャンパラメータ設定方法およびｘ線ｃｔ装置 |
| JP2007054372A (ja) * | 2005-08-25 | 2007-03-08 | Ge Medical Systems Global Technology Co Llc | Ｘ線ｃｔ装置 |
| US7532702B2 (en) * | 2005-11-23 | 2009-05-12 | General Electric Company | Method and system for performing CT image reconstruction with motion artifact correction |
| JP5317389B2 (ja) * | 2006-02-03 | 2013-10-16 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 放射線断層撮影装置 |
| WO2007138979A1 (ja) * | 2006-05-25 | 2007-12-06 | Hitachi Medical Corporation | Ｘ線ｃｔ装置 |
| US20080019476A1 (en) * | 2006-07-24 | 2008-01-24 | Armen Mirzayan, Dds, Inc. | Apparatus and Method for Producing Medical X-ray Images |
| JP4783256B2 (ja) * | 2006-10-13 | 2011-09-28 | 株式会社リガク | Ｃｔデータ処理装置およびｃｔデータ処理プログラム |
| JP2008113960A (ja) * | 2006-11-07 | 2008-05-22 | Ge Medical Systems Global Technology Co Llc | 放射線撮影装置 |
| WO2008056814A1 (en) | 2006-11-09 | 2008-05-15 | Canon Kabushiki Kaisha | Radiographic imaging control apparatus using multi radiation generating apparatus |
| WO2008096813A1 (ja) * | 2007-02-08 | 2008-08-14 | Hitachi Medical Corporation | Ｘ線ｃt装置 |
| JP5039064B2 (ja) * | 2007-02-14 | 2012-10-03 | 株式会社日立メディコ | Ｘ線ｃｔ装置 |
| JP4729519B2 (ja) * | 2007-03-09 | 2011-07-20 | ゼネラル・エレクトリック・カンパニイ | 器官に基づく放射線プロファイル設定を設けた放射線撮像の方法及びシステム |
| JP5184834B2 (ja) * | 2007-07-18 | 2013-04-17 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | Ｘ線ｃｔ装置 |
| JP2009089739A (ja) * | 2007-10-03 | 2009-04-30 | Toshiba Corp | 乳房撮影検査用ｘ線診断装置 |
| JP5389345B2 (ja) * | 2007-10-04 | 2014-01-15 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | Ｘ線ｃｔ装置 |
| KR20100133950A (ko) | 2008-01-30 | 2010-12-22 | 더 리전츠 오브 더 유니버시티 오브 캘리포니아 | 동적인 제약들에 따른 오브젝트 주변의 사용을 통한 단층 촬영에 있어서의 양 감소 및 이미지 강화 |
| JP4542591B2 (ja) * | 2008-03-05 | 2010-09-15 | 富士フイルム株式会社 | 画像処理システム、画像保管装置、および、画像処理装置、並びに、画像処理方法、プログラム |
| DE102008014738A1 (de) * | 2008-03-18 | 2009-09-24 | Siemens Aktiengesellschaft | Verfahren zur medizinischen Bildgebung sowie medizinische Bildgebungsvorrichtung |
| JP5523726B2 (ja) * | 2008-04-04 | 2014-06-18 | 株式会社東芝 | Ｘ線ｃｔ装置 |
| JP5455903B2 (ja) * | 2008-07-04 | 2014-03-26 | 株式会社日立メディコ | Ｘ線ｃｔ装置 |
| US8300765B2 (en) * | 2008-08-04 | 2012-10-30 | Koninklijke Philips Electronics N.V. | Interventional imaging and data processing |
| JP5404767B2 (ja) * | 2009-03-04 | 2014-02-05 | 株式会社日立メディコ | Ｘ線ｃｔ装置及び画像修正プログラム |
| JP2010213798A (ja) * | 2009-03-13 | 2010-09-30 | Toshiba Corp | 循環器用ｘ線診断システム |
| JP4644292B2 (ja) * | 2009-04-16 | 2011-03-02 | 株式会社日立メディコ | Ｘ線ｃｔ装置とその画像表示方法 |
| JP5279637B2 (ja) * | 2009-07-02 | 2013-09-04 | 株式会社日立製作所 | ベッド位置決めシステム及びベッド位置決め方法 |
| JP5442363B2 (ja) * | 2009-08-28 | 2014-03-12 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | Ｘ線ｃｔ装置 |
| US8494814B2 (en) | 2009-11-02 | 2013-07-23 | Siemens Corporation | Methods for inlet and outlet modelling of vent as large as possible for hearing aid devices |
| JP5481185B2 (ja) * | 2009-12-28 | 2014-04-23 | 株式会社東芝 | Ｘ線診断装置 |
| JP2011172847A (ja) * | 2010-02-25 | 2011-09-08 | Fujifilm Corp | 放射線画像生成装置及び放射線画像生成方法 |
| US8693634B2 (en) | 2010-03-19 | 2014-04-08 | Hologic Inc | System and method for generating enhanced density distribution in a three dimensional model of a structure for use in skeletal assessment using a limited number of two-dimensional views |
| DE102010025512B4 (de) * | 2010-06-29 | 2016-02-25 | Siemens Aktiengesellschaft | Verfahren zum Darstellen einer Patientendosis und Röntgengerät |
| JP5642439B2 (ja) * | 2010-07-06 | 2014-12-17 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 放射線断層撮像装置 |
| US8768029B2 (en) * | 2010-10-20 | 2014-07-01 | Medtronic Navigation, Inc. | Selected image acquisition technique to optimize patient model construction |
| JP5718014B2 (ja) * | 2010-10-22 | 2015-05-13 | 株式会社日立メディコ | Ｘ線ｃｔ装置 |
| US8503613B2 (en) * | 2010-11-24 | 2013-08-06 | General Electric Company | Dose level indication |
| DK3524159T3 (da) * | 2010-12-08 | 2021-04-06 | Bayer Healthcare Llc | Frembringelse af en egnet model til estimering af strålingsdosis modtaget af patient under scanninger til medicinsk billeddannelse |
| EP2465435B1 (en) * | 2010-12-14 | 2019-12-04 | General Electric Company | Selection of optimal viewing angle to optimize anatomy visibility and patient skin dose |
| US9089265B2 (en) * | 2011-01-14 | 2015-07-28 | The Cleveland Clinic Foundation | Automated parameter selection for a tomographic imaging device |
| JP5985836B2 (ja) * | 2011-03-03 | 2016-09-06 | ゼネラル・エレクトリック・カンパニイ | イメージング・システムによって放出される放射線量を減少させる方法 |
| US9326738B2 (en) * | 2011-06-30 | 2016-05-03 | General Electric Company | Method and system for reduced dose X-ray imaging |
| JP6125206B2 (ja) * | 2011-12-02 | 2017-05-10 | 東芝メディカルシステムズ株式会社 | Ｘ線画像診断装置 |
| JP6129474B2 (ja) * | 2012-02-09 | 2017-05-17 | 東芝メディカルシステムズ株式会社 | Ｘ線診断装置 |
| EP2854644A4 (en) * | 2012-06-05 | 2016-03-16 | Rapiscan Systems Inc | SOURCE TRIGGER PATTERN OPTIMIZATION FOR X-RAY SCANNING SYSTEMS |
| JP6257948B2 (ja) * | 2012-08-07 | 2018-01-10 | 東芝メディカルシステムズ株式会社 | Ｘ線撮影システム |
| KR20140044174A (ko) * | 2012-10-04 | 2014-04-14 | 주식회사바텍 | 엑스선 촬영장치 및 촬영방법 |
| JP6117524B2 (ja) | 2012-11-22 | 2017-04-19 | 富士フイルム株式会社 | 放射線画像撮影システム、放射線画像撮影方法及びプログラム |
| JP6242683B2 (ja) * | 2012-12-27 | 2017-12-06 | 東芝メディカルシステムズ株式会社 | Ｘ線ｃｔ装置及び制御方法 |
| JP5805689B2 (ja) | 2013-03-08 | 2015-11-04 | 株式会社モリタ製作所 | Ｘ線ｃｔ撮影装置及びｘ線ｃｔ撮影方法 |
| JP6165511B2 (ja) * | 2013-06-12 | 2017-07-19 | 東芝メディカルシステムズ株式会社 | Ｘ線コンピュータ断層撮影装置 |
| KR101534098B1 (ko) * | 2013-09-13 | 2015-07-07 | 삼성전자주식회사 | Ｃｔ 장치 및 이를 이용한 엑스선 제어 방법 |
| DE102013219249A1 (de) * | 2013-09-25 | 2015-03-26 | Siemens Aktiengesellschaft | Verfahren und System zur automatischen Auswahl eines Scanprotokolls |
| WO2015065843A1 (en) * | 2013-10-31 | 2015-05-07 | Hexagon Metrology, Inc. | Parametric control of object scanning |
| CN104287768A (zh) * | 2014-09-30 | 2015-01-21 | 沈阳东软医疗系统有限公司 | 一种ct扫描剂量控制方法及系统 |
| GB2533632B (en) * | 2014-12-24 | 2018-01-03 | Gen Electric | Method and system for obtaining low dose tomosynthesis and material decomposition images |
| WO2017103238A1 (en) * | 2015-12-17 | 2017-06-22 | Koninklijke Philips N.V. | Method for estimating the radiation dose received by an organ during a computed tomography scan |
| US10085698B2 (en) * | 2016-01-26 | 2018-10-02 | Genereal Electric Company | Methods and systems for automated tube current modulation |
| JP6342437B2 (ja) * | 2016-02-22 | 2018-06-13 | ゼネラル・エレクトリック・カンパニイ | 放射線断層撮影システム及びその制御プログラム |
| JP7055599B2 (ja) * | 2016-06-06 | 2022-04-18 | キヤノンメディカルシステムズ株式会社 | Ｘ線ｃｔ装置 |
| US10845615B2 (en) * | 2016-09-29 | 2020-11-24 | Marel Iceland Ehf | Method of generating a three dimensional surface profile of a food object |
| CN106725570B (zh) * | 2016-12-30 | 2019-12-20 | 上海联影医疗科技有限公司 | 成像方法及系统 |
| JP6125126B1 (ja) * | 2017-01-18 | 2017-05-10 | 典人 畠山 | Ｘ線撮影条件補正方法、プログラム、及び、ｘ線システム |
| JP6815209B2 (ja) * | 2017-01-24 | 2021-01-20 | キヤノンメディカルシステムズ株式会社 | Ｘ線ｃｔ装置 |
| US10806418B2 (en) | 2017-05-31 | 2020-10-20 | Canon Medical Systems Corporation | X-ray CT apparatus and imaging condition calculating method |
| US10973489B2 (en) * | 2017-09-29 | 2021-04-13 | General Electric Company | CT imaging system and method using a task-based image quality metric to achieve a desired image quality |
| JP7282487B2 (ja) * | 2018-06-07 | 2023-05-29 | キヤノンメディカルシステムズ株式会社 | 医用画像診断装置 |
| GB2604553B (en) * | 2019-11-14 | 2023-08-02 | Eos Imaging | Radiological imaging method |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2001043993A (ja) * | 1999-07-30 | 2001-02-16 | Ge Yokogawa Medical Systems Ltd | 管電流調節方法および装置並びにｘ線ｃｔ装置 |
Family Cites Families (19)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS53110495A (en) | 1977-03-09 | 1978-09-27 | Toshiba Corp | Radiation tomograph ic device |
| JPS5894833A (ja) * | 1981-12-01 | 1983-06-06 | 株式会社東芝 | Ｘ線ｃｔ装置 |
| JPS59155238A (ja) * | 1983-02-24 | 1984-09-04 | 株式会社東芝 | Ｘ線ｃｔ装置 |
| JPS6152860A (ja) * | 1984-08-24 | 1986-03-15 | 株式会社 日立メデイコ | Ｘ線ｃｔ装置 |
| JPH0636793B2 (ja) * | 1988-05-20 | 1994-05-18 | 株式会社東芝 | Ｘ線ｃｔスキヤナ |
| JP2933107B2 (ja) * | 1992-02-14 | 1999-08-09 | 株式会社東芝 | Ｃｔ装置 |
| JP3449561B2 (ja) * | 1993-04-19 | 2003-09-22 | 東芝医用システムエンジニアリング株式会社 | Ｘ線ｃｔ装置 |
| JP3466678B2 (ja) * | 1993-11-01 | 2003-11-17 | 株式会社東芝 | Ｘ線ｃｔスキャナ |
| US5450462A (en) * | 1993-11-19 | 1995-09-12 | General Electric Company | Modulation of x-ray tube current during CT scanning with modulation limit |
| JP3483157B2 (ja) * | 1994-01-24 | 2004-01-06 | 株式会社日立メディコ | Ｘ線ｃｔ装置 |
| US5485494A (en) * | 1994-08-03 | 1996-01-16 | General Electric Company | Modulation of X-ray tube current during CT scanning |
| JP2914891B2 (ja) * | 1995-07-05 | 1999-07-05 | 株式会社東芝 | Ｘ線コンピュータ断層撮影装置 |
| US5625662A (en) * | 1995-11-20 | 1997-04-29 | General Electric Company | Modulating x-ray tube current in a CT system |
| JP2768932B2 (ja) * | 1996-10-14 | 1998-06-25 | 株式会社東芝 | Ｘ線ｃｔ装置 |
| US5867555A (en) | 1997-03-04 | 1999-02-02 | Siemens Aktiengesellschaft | Adaptive dose modulation during CT scanning |
| JP4497570B2 (ja) * | 1998-01-22 | 2010-07-07 | 株式会社東芝 | 画像診断装置 |
| JPH11206755A (ja) * | 1998-01-30 | 1999-08-03 | Shimadzu Corp | Ｘ線ｃｔ装置 |
| US6307912B1 (en) * | 1999-11-29 | 2001-10-23 | General Electric Company | Methods and apparatus for optimizing CT image quality with optimized data acquisition |
| JP4519254B2 (ja) | 2000-04-03 | 2010-08-04 | 株式会社日立メディコ | Ｘ線ｃｔ装置 |
-
2001
- 2001-03-09 JP JP2001067150A patent/JP4532005B2/ja not_active Expired - Lifetime
-
2002
- 2002-03-08 US US10/468,711 patent/US7103139B2/en not_active Expired - Lifetime
- 2002-03-08 WO PCT/JP2002/002190 patent/WO2002071946A1/ja active Application Filing
- 2002-03-08 EP EP02702836A patent/EP1393681B1/en not_active Expired - Lifetime
- 2002-03-08 EP EP08018586A patent/EP2047799B1/en not_active Expired - Lifetime
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2001043993A (ja) * | 1999-07-30 | 2001-02-16 | Ge Yokogawa Medical Systems Ltd | 管電流調節方法および装置並びにｘ線ｃｔ装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| EP2047799A2 (en) | 2009-04-15 |
| WO2002071946A1 (fr) | 2002-09-19 |
| JP2002263097A (ja) | 2002-09-17 |
| US20040086076A1 (en) | 2004-05-06 |
| EP1393681A1 (en) | 2004-03-03 |
| EP2047799A3 (en) | 2009-06-24 |
| EP1393681A4 (en) | 2006-05-31 |
| EP2047799B1 (en) | 2012-10-10 |
| EP1393681B1 (en) | 2011-08-17 |
| US7103139B2 (en) | 2006-09-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4532005B2 (ja) | Ｘ線ｃｔ装置及びその画像表示方法 | |
| JP4519254B2 (ja) | Ｘ線ｃｔ装置 | |
| JP2024100832A (ja) | イメージングにおける散乱評価および散乱補正を向上させるための方法および機器 | |
| JP3124254B2 (ja) | 放射線断層撮影装置 | |
| US8483363B2 (en) | Movable wedge for improved image quality in 3D X-ray imaging | |
| WO2007074772A1 (ja) | Ｘ線ｃｔ装置 | |
| US7344306B2 (en) | Systems and methods for compensating for table sag | |
| JP5039064B2 (ja) | Ｘ線ｃｔ装置 | |
| CN107252318A (zh) | 采集协议评价设备 | |
| JP2006068538A (ja) | 放射線療法のための送出修正システム | |
| JP4468352B2 (ja) | コンピュータトモグラフィにおける局所的患者線量の再構成 | |
| JP4729519B2 (ja) | 器官に基づく放射線プロファイル設定を設けた放射線撮像の方法及びシステム | |
| JP5329256B2 (ja) | ベッド位置決めシステム、放射線治療システム及びベッド位置決め方法 | |
| JP4429694B2 (ja) | Ｘ線ｃｔ装置 | |
| JP2003070780A (ja) | コンピュータトモグラフによる測定データの取得方法およびコンピュータトモグラフ | |
| JP4554185B2 (ja) | Ｘ線ｃｔ装置 | |
| JP2001218767A (ja) | スキャノグラムを形成するコンピュータ断層撮影方法 | |
| WO2006090321A1 (en) | Determination of the coverage of a ct scan | |
| JP5027909B2 (ja) | Ｘ線ｃｔ装置 | |
| JP4644292B2 (ja) | Ｘ線ｃｔ装置とその画像表示方法 | |
| JP3789728B2 (ja) | プロジェクションデータ補正方法および装置並びに放射線断層撮像装置 | |
| CN107809954B (zh) | 计算机断层切片图像相对于要被成像的对象的深度位置的显示 | |
| WO2023176264A1 (ja) | 医用画像処理装置、治療システム、医用画像処理方法、およびプログラム | |
| JP5384293B2 (ja) | Ｘ線ｃｔ装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20080307 |
|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20080307 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20090615 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20090812 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20091028 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20091224 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20100308 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20100427 |
|
| 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: 20100607 |
|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20100610 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 4532005 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130618 Year of fee payment: 3 |
|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313111 |
|
| S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| EXPY | Cancellation because of completion of term |