JP4188042B2 - X-ray CT apparatus and dose detection apparatus - Google Patents

Description

【０００５】
ここで、Ｔ：スライス厚
Ｄ（ｚ）：線量プロフィール（dose profile）
次に、１スライス及び１ｍＡｓ当たりの規格化された平均線量ｎＣＴＤＩｗは、ファントムの中央部ｃと周辺部ｐの各ＣＴＤＩｗ値を所定の割合で加算（加重平均）することにより（２）式で与えられる。
【０００６】
【数２】

【０００７】
ここで、Ｃ：Ｘ線放射線量（ｍＡｓ）
ＣＴＤＩ_100,ｃ：ファントム中心部ｃの厚さ１００ｍｍに制限されたＣＴ ＤＩ値
ＣＴＤＩ_100,ｐ：ファントム周辺部ｐの厚さ１００ｍｍに制限された４つ のＣＴＤＩ値の平均値
そして、臨床(診断）で使用される１スライス当たりの平均線量ＣＴＤＩｗは、ｎＣＴＤＩｗにｍＡｓを掛けたものとして、（３）式で与えられる。
【０００８】
【数３】

【０００９】
ここで、Ｃ：Ｘ線放射線量（ｍＡｓ）
更に、スライス数Ｎを考慮した場合の被曝線量ＤＬＰ(Dose-Length Product)は（４）式で表される。
【００１０】
【数４】

【００１１】
ここで、ｉ：各シリアルスキャンシーケンス
Ｔ：各スライス厚（ｃｍ）
Ｎ：スライス数
また、ヘリカルスキャンの場合の被曝線量ＤＬＰは（５）式で表される。
【００１２】
【数５】

【００１３】
ここで、Ｔ：公称スライス厚（ｃｍ）
Ａ：管電流（ｍＡ）
ｔ：トータルの曝射時間
このような規格の下で、従来は、予め身体用及び頭部用の各基準ファントム（ＣＴ dosimetry phantom）を使用して測定したＣＴＤＩｗ値に基づき、被検体のＣＴ撮影による被曝線量を、その被検体サイズ（体型，体格等）によらず一律に表示していた。即ち、被検体の頭部を撮影した場合は、予め頭部用ファントムを使用して求められたｎＣＴＤＩｗを基に、被検体の被曝線量（ＣＴＤＩｗ，ＤＬＰ）を求めて表示し、また被検体の胴体部を撮影した場合は、予め胴体用ファントムを使用して求められたｎＣＴＤＩｗを基に、被検体の被曝線量（ＣＴＤＩｗ，ＤＬＰ）を求めて表示していた。
【００１４】
しかし、上記従来方式によると、特にサイズが基準ファントムよりも小さいような幼児や子供等に対しては、被曝線量が実際よりも低く見積もられる傾向があり、非常に危険であった。
【００１５】
この点、本件出願人は、予め被検体サイズを検出することにより、最適のＸ線照射条件（管電流）を決定するＸ線ＣＴ装置を既に提案している(例えば特許文献１，２）。
【００１６】
【特許文献１】
特開２００１−４３９９３号公報（段落「００３４」〜「００３６」，図６）。
【００１７】
【特許文献２】
特開２００１−２１８７６１号公報（段落「００９３」〜「００９７」，図６，図７）。
【００１８】
【発明が解決しようとする課題】
しかし、上記何れの場合も、予め被検体サイズを検出するために、事前の透視撮影（スカウトスキャン）が必要であり、被検体に余分な被曝を与える不利益があった。
【００１９】
本発明は上記従来技術の問題点に鑑みなされたもので、その目的とする所は、被検体の実際のＣＴ撮影に忠実な被曝線量を提供できるＸ線ＣＴ装置及び線量検出方法並びにプログラムを提供することにある。
【００２０】
【課題を解決するための手段】
上記の課題は例えば図１の構成により解決される。即ち、本発明（１）のＸ線ＣＴ装置は、被検体を挟んで相対向するようにＸ線管及びＸ線検出器を備え、Ｘ線検出器の検出信号に基づき被検体のＣＴ断層像を再構成するＸ線ＣＴ装置において、前記被検体のＣＴ断層像につき、所定の組織の部分を抽出する抽出手段と、前記抽出した組織のＣＴ断層像につき求めた画像雑音から対応する位置の線量を求める線量演算手段と、前記求めた線量の情報を出力する出力手段とを備えるものである。
【００２１】
ところで、Ｘ線ＣＴ装置のように、被検体を透過するＸ線光子数をできるだけ少なくしたいような撮影系のＣＴ断層像では、被検体を透過するＸ線光子（フォトン）数の時間的、空間的な揺らぎ（光子雑音）が画像雑音の主な要因となり、この場合に、ＣＴ値の信号成分が線量（光子数の平均値）に比例するとすると、ＣＴ値の雑音成分（標準偏差）は、ポアソン過程に従い、線量の平方根に比例することになる。
【００２２】
本発明（１）では、線量演算手段３は、被検体１００のＣＴ断層像につき求めた画像雑音から対応する位置（部位）の線量を求めるため、実際に被検体１００をＣＴ撮影した際の各位置の線量、即ち、実際の被検体サイズや撮影状況が加味された各位置の被曝線量をより忠実に提供できる。従って、幼児等のサイズの小さい被検体に対して過被曝となってしまうような状況を有効に回避できる。また、従来の様に、単に１スライス当たりの平均線量ＣＴＤＩｗを一律に提供するのではなく、ＣＴ断層像を構成する各位置（組織）毎の被曝線量を提供できるため、臨床上からも、より多様で有用な被曝情報が得られる。
【００２３】
また本発明（２）では、上記本発明（１）において、線量演算手段３は、ＣＴ断層像の画像雑音をＣＴ値の標準偏差σ₂で表し、所定の係数パラメータをｋとするときに、対応する位置の線量ＣＴＤＩ₂を、ＣＴＤＩ₂＝（σ₂／ｋ）² の関係により求めるものである。
【００２４】
本発明（２）によれば、予めＸ線ＣＴ装置に固有の係数パラメータｋを求めておくことにより、被検体の各部位の被曝線量ＣＴＤＩ₂を、各対応する部位のＣＴ値の標準偏差σ₂から容易に求められる。
【００２５】
また本発明（３）では、上記本発明（１）において、線量演算手段３は、被検体１００のＣＴ断層像につき求めた画像のＳ／Ｎ値から対応する位置の線量を求めるものである。
【００２６】
本発明（３）によれば、ＣＴ断層像における各部位の信号対雑音（Ｓ／Ｎ）比は、対応する部位の線量の平方根に比例することになるため、上記本発明（１）と同様に、実際に被検体１００をＣＴ撮影した際の各位置の線量をより忠実に提供できる。また逆に、ＣＴ断層像につき所要の画像のＳ／Ｎ値を得るための線量を有効に制限できる。
【００２７】
また本発明（４）では、上記本発明（１）において、被検体のＣＴ断層像を画面に表示する表示手段と、表示画面上でＣＴ断層像の前記組織を指定させる指定手段とを備え、線量演算手段３は、前記指定された領域のＣＴ断層像につき求めた画像雑音から該指定領域に対する線量を求める。従って、任意の組織（脳，脂肪，肝臓，肺等）に対する被曝線量を容易に検出できる。
【００２８】
また本発明（５）では、上記本発明（１）において、前記抽出手段は、被検体のＣＴ断層像につき、別途に入力されたＣＴ閾値の範囲に含まれる組織の部分を抽出するものであり、前記抽出した組織のＣＴ断層像を他の組織のＣＴ断層像と区別できる態様で画面に表示する表示手段を備えるものである。
【００２９】
本発明（５）によれば、技師はＣＴ閾値の範囲を指定するだけで対応する組織（脳，眼球，骨，脂肪，肝臓，肺等）の複雑な部位を容易に抽出でき、その被曝線量を容易に検出できる。
【００３０】
また本発明（６）では、上記本発明（１）において、線量演算手段３は、被検体のＣＴ断層像につき各画素位置の対応に求めた画像雑音から各画素位置に対応する線量を求め、出力手段４は、前記求めた各線量値で構成される被曝線量画像を出力する。
【００３１】
本発明（６）によれば、ＣＴ断層像の各画素位置の対応に求めた被曝線量が得られるため、被検体の被曝状況を精密（高分解能）に把握できる。
【００３２】
また本発明（７）では、上記本発明（６）において、出力手段４は、線量値に応じて色分けした線量画像を出力する。従って、精密な被曝状況を容易に認識できる。
【００３３】
また、本発明（８）では、上記本発明（１）において、線量演算手段３は、被検体のＣＴ断層像の中心部及び周辺部を含む複数の位置対応に求めた各画像雑音から各対応する位置の線量を求め、これらの重み付け平均により１スライス分の平均線量値を求める。
【００３４】
本発明（８）によれば、従来より一般に採用されているような、上記（２），（３）式による線量値ＣＴＤＩｗと同様の線量値を提供でき、従来システムとの診断の共通性も容易に確保できる。
【００３５】
また本発明（９）では、上記本発明（８）において、線量演算手段３は、スライス毎に求めた各平均線量値をスライス数分加算することで１撮影分の線量値を求める。
【００３６】
本発明（９）によれば、従来と同様に、１撮影分の被曝線量値を容易に提供できる。また、この場合でも、スライス毎に求めた各平均線量値をスライス数分加算することで、１撮影分の線量値をより精密に提供できる。
【００３７】
また、本発明（１０）では、上記本発明（２）において、基準となるファントム６０のＣＴ断層像につき求めた画像雑音をＣＴ値の標準偏差σ₁で表し、かつ前記ファントム中の対応位置で実測された線量をＣＴＤＩ₁とするときに、前記画像雑音σ₁を前記線量ＣＴＤＩ₁に変換するための係数パラメータｋを、ｋ＝σ₁／√ＣＴＤＩ₁の関係により求める係数演算手段２を備えるものである。
【００３８】
本発明（１０）によれば、基準ファントム６０のＣＴ断層像につき求めた画像雑音（標準偏差）σ₁と、基準ファントム６０中で実測された線量との間を、予め装置に固有の係数パラメータｋの形で関係付けておく構成により、その後は、被検体１００のＣＴ断層像に含まれる画像雑音から、実際のＣＴ撮影に応じた被曝線量を容易に正確に求められる。また、本発明方法によりＣＴ断層像から求められる被曝線量値も、基準ファントム６０で測定された実測値に近いものとなる。
【００３９】
また、本発明（１１）では、上記本発明（１０）において、基準となるファントム６０は、均質な材料からなる円柱様体の中心軸及びその外周付近の各所定位置の前記中心軸と平行に複数の線量センサ６２を備えるものである。従って、材料の均質な円柱様体からはＣＴ断層像の画像雑音を精密に求められる。
【００４０】
また、本発明（１２）では、上記本発明（１１）において、線量センサ６２はペンシル型イオンチェンバよりなる。
【００４１】
また、本発明（１３）では、上記本発明（１２）において、線量演算手段３は、前記ファントム６０のＣＴ断層像につきその中心部及び周辺部を含む複数の位置対応に求めた各画像雑音の加重平均をＣＴ値の標準偏差σで表し、かつ前記ファントム中の各対応位置で実測され、前記と同じ方法で加重平均された１スライス分の平均線量をＣＴＤＩとするときに、前記画像雑音σを前記線量ＣＴＤＩに変換するための係数パラメータｋを、ｋ＝σ／√ＣＴＤＩの関係により求める。
【００４２】
本発明（１３）によれば、基準ファントム６０につき上記（２）式と同様の方法で加重平均された画像雑音σを、該（２）式と同じ方法で加重平均された平均線量値ＣＴＤＩに簡単な係数パラメータｋで関係つけることにより、従来より、一般に採用されていると同様の被曝線量値ＣＴＤＩｗを容易に求められる。
【００４３】
また、本発明（１４）の線量検出装置は、被検体のＣＴ断層像につき、所定の組織の部分を抽出する抽出手段と、前記抽出した組織のＣＴ断層像につき求めた画像雑音から対応する位置の線量を求める線量演算手段と、前記求めた線量の情報を出力する出力手段とを備えるものである。
また、本発明の線量検出方法は、被検体１００を挟んで相対向するようにＸ線管４０及びＸ線検出器９０を備え、Ｘ線検出器の検出信号に基づき被検体のＣＴ断層像を再構成するＸ線ＣＴ装置の線量検出方法であって、予め基準となるファントム６０をスキャンしたＣＴ断層像につき求めた画像雑音と、前記ファントム中で実測された線量値とに基づき、前記画像雑音を前記線量値に変換するための係数パラメータを求めるステップと、被検体１００をスキャンしたＣＴ断層像につき求めた画像雑音から前記係数パラメータを用いて対応する位置の線量を求めるステップとを備えるものである。
【００４４】
好ましくは、基準となるファントム６０のＣＴ断層像につき求めた画像雑音をσ_１、かつ該ファントムの対応位置で実測された線量値をＣＴＤＩ_１とするときに、係数パラメータｋを、ｋ＝σ_１／√ＣＴＤＩ_１の関係により求めるステップと、被検体１００のＣＴ断層像につき求めた画像雑音をσ_２とするときに、対応する位置の線量値ＣＴＤＩ_２を、ＣＴＤＩ_２＝（σ_２／ｋ）^２の関係により求めるステップとを備えるものである。
【００４５】
また本発明のプログラムは、コンピュータに上記の線量検出方法を実行させるためのコンピュ−タ実行可能なプログラムである。このようなプログラムは、ＣＤ−ＲＯＭ等の記録媒体に記録して、又は通信回線を介したオンライン通信により提供可能である。
【００４６】
【発明の実施の形態】
以下、添付図面に従って本発明に好適なる実施の形態を詳細に説明する。なお、全図を通して同一符号は同一又は相当部分を示すものとする。
【００４７】
図２は実施の形態によるＸ線ＣＴ装置の要部構成図で、該装置は、Ｘ線ファンビームＸＬＦＢにより被検体１００のアキシャル／ヘリカルスキャン・読取等を行う走査ガントリ部３０と、被検体１００を載せて体軸ＣＬｂの方向に移動させる撮影テーブル２０と、前記走査ガントリ部３０及び撮影テーブル２０の制御を行うと共に、Ｘ線撮影技師が各種の設定・操作を行う遠隔の操作コンソール部１０とを備える。
【００４８】
走査ガントリ部３０において、４０は回転陽極型のＸ線管、４０ＡはＸ線管の管電圧ｋＶ、管電流ｍＡ等を制御するＸ線管制御部、５０はＸ線のスライス厚を制限するコリメータ、５０Ａはコリメータ制御部、９０はチャネルＣＨ方向に並ぶ多数（ｎ＝１０００程度）のＸ線検出素子が体軸ＣＬｂ方向の例えば２列Ｌ１，Ｌ２に配列されているＸ線検出器（マルチディテクタ）、９１はＸ線検出器９０の検出信号に基づき被検体１００の投影データｇ₁（Ｘ，θ），ｇ₂（Ｘ，θ）を生成し、収集するデータ収集部（ＤＡＳ）、３５はこれらＸ線撮影系に係る各機器を体軸ＣＬｂの回りに回転自在に支持するガントリ、３５Ａはガントリ３５の回転制御部である。
【００４９】
操作コンソール部１０において、１１はＸ線ＣＴ装置の主制御・処理（スキャン制御、ＣＴ断層像の再構成、Ｘ線被曝量表示に関する各種処理）を行う中央処理装置、１１ａはそのＣＰＵ、１１ｂはＣＰＵ１１ａが使用するＲＡＭ，ＲＯＭ等からなる主メモリ（ＭＭ）、１２はキーボードやマウス等を含む指令やデータの入力装置、１３はスキャン計画情報、ＣＴ断層像、Ｘ線被曝量に関する各種情報を表示するための表示装置（ＣＲＴ）、１４はＣＰＵ１１ａと走査ガントリ部３０及び撮影テーブル２０等との間で各種制御信号ＣＳやモニタ信号ＭＳのやり取りを行う制御インタフェース、１５はデータ収集部９１からの投影データを一時的に蓄積するデータ収集バッファ、１６はスキャン（投影）データやＣＴ断層像のデータを最終的に蓄積・格納すると共に、Ｘ線ＣＴ装置の運用に必要な各種アプリケーションプログラムや各種演算／補正用のデータファイル等を格納している２次記憶装置（ハードディスク装置等）である。
【００５０】
被検体ＣＴ撮影の基本的な動作を概説すると、Ｘ線管４０からのＸ線ファンビームＸＬＦＢは被検体１００を透過してＸ線検出器９０の検出列Ｌ１，Ｌ２に一斉に入射する。データ収集部９１はＸ線検出器９０の各検出出力に対応する投影データｇ₁（Ｘ，θ），ｇ₂（Ｘ，θ）を生成し、これらをデータ収集バッファ１５に格納する。ここで、ＸはＸ線検出器９０の検出チャネル１〜ｎ、θは体軸ＣＬｂの周りのビュー角を表す。更に、ガントリ３５が僅かに回転した各ビュー角θで上記同様の投影を行い、こうしてガントリ１回転分の投影データを収集・蓄積する。また同時に、アキシャル／ヘリカルスキャン方式に従って撮影テーブル２０を体軸ＣＬｂ方向に間欠的／連続的に移動させ、こうして被検体１００の所要撮影領域についての全投影データを収集・蓄積する。そして、ＣＰＵ１１ａは、上記全スキャンの終了後、又はスキャン実行に追従（並行）して、得られた投影データに基づき被検体１００のＣＴ断層像を再構成し、これを表示装置１３に表示する。そして、本実施の形態ではこのＣＴ断層像のデータを利用して、被検体の被爆線量を求め、画面に表示可能である。以下、詳細に説明する。
【００５１】
図５は実施の形態による装置パラメータ決定処理のフローチャート、図６，図７は該処理のイメージ図（１），（２）であり、予め基準ファントムを読み取った際のＣＴ断層像に含まれる画像雑音から対応する被曝線量を求めるための装置固有のパラメータｋを決定する処理を示している。
【００５２】
図５において、ステップＳ１１では寝台２０の先端部にＣＴ線量測定用の基準ァントムを装着する。その装着イメージを図６（Ａ）に示す。また図７（Ａ）にＣＴ線量測定用ファントム（基準ファントム）の例を示す。図７（Ａ）において、胴部用ファントム６０ｂのボディー６１は、均質なアクリル樹脂（ＰＰＭＡ：polymethylmethacrylate）の円柱様体からなっており、その直径φは、被検体の典型的な胴部サイズを代表するようなφ＝３２ｃｍ、及びその厚さは１４ｃｍ以上となっている。更に、このファントム６０ｂには、ボディー６１のｚ軸（被検体体軸ＣＬｂに相当）と平行な中央部ｃ及びその外周から各１ｃｍ内側の４箇所ｐに、夫々実行長が１０ｃｍのペンシル型イオンチェンバ（pencil ionisation chamber）６２が挿入されている。一方、頭部用ファントム６０ｈの直径φは、被検体の典型的な頭部サイズを代表するようなφ＝１６ｃｍとなっている。その他の構成については上記ファントム６０ｂと同様である。
【００５３】
図５に戻り、ステップＳ１２では技師が基準ファントム（例えば６０ｂ）に対するスキャンパラメータ（管電圧ｋＶ，管電流ｍＡ，Ｘ線曝射時間ｓ等）を設定する。ステップＳ１３では基準ファントム６０ｂのスキャン（１スライス分）を行い、ステップＳ１４ではその投影データを収集する。ステップＳ１５では基準ファントム６０ｂのＣＴ断層像を再構成する。得られたＣＴ断層像の各画素データは均質なアクリル樹脂のＣＴ値からなり、（６）式で表される。
【００５４】
【数６】

【００５５】
ここで、μ_w：水の線減弱係数
μ_m：物質（この例ではアクリル樹脂）の線減弱係数
上記、均質なアクリル樹脂をスキャンしたので、得られたＣＴ値も均一のはずである。しかし、実際にはファントム６０ｂを透過するＸ線光子数の時間的、空間的な揺らぎ（光子雑音）により、再構成したＣＴ値にも光子雑音に起因する揺らぎが生じ、これがＣＴ断層像の画像雑音となって表れる。
【００５６】
図６（Ｂ）に基準ファントム６０ｂの画像雑音を含むＣＴイメージを示す。Ｘ線ＣＴ装置のように、被検体１００を透過するＸ線光子数（被曝線量）をできるだけ少なくしたいような撮影系のＣＴ断層像では、被検体を透過するＸ線光子（フォトン）数の時間的、空間的な揺らぎ（光子雑音）が画像雑音の主な要因となっており、この場合に、ＣＴ値の信号成分が線量（光子数の平均値）に比例するとすると、ＣＴ値の雑音成分（標準偏差）は、ポアソン過程に従い、線量の平方根に比例する。即ち、ＣＴ断層像の信号対雑音比（Ｓ／Ｎ）は、（７）式で表される。
【００５７】
【数７】

【００５８】
ここで、ａ：定数
Dose：線量
また、画像雑音Ｎと線量Doseとの間には（８）式の関係がある。
【００５９】
【数８】

【００６０】
ここで、ｋ：定数
今、基準ファントム６０ｂの各位置ｃ，ｐにおける線量をＤose_c，Ｄose_pで表し、各位置ｃ，ｐにおける画像雑音Ｎを夫々ＣＴ値についての標準偏差σ_c，σ_pで評価したとすると、画像雑音σ_c，σ_pと線量Ｄose_c，Ｄose_pとの間には（９）式の関係がある。
【００６１】
【数９】

【００６２】
ここで、Ｑ：装置構成によって決まる係数
Ｋ：再構成関数によって決まる係数
基準ファントム６０ｂ内の各位置ｃ，ｐにおける線量ＣＴＤＩ_100,ｃ，ＣＴＤＩ_100,ｐは各位置のイオンチェンバ６２により測定可能であるから、図５のステップＳ１６では各イオンチェンバ６２の測定データを収集し、ステップＳ１７では各位置ｃ，ｐにおける線量値ＣＴＤＩ_100,ｃ，ＣＴＤＩ_100,ｐを求める。
【００６３】
図７（Ｂ）に線量プロフィールＤ（ｚ）の例を示す。線量プロフィールＤ（ｚ）は、１回転のスキャンによりイオンチェンバ６２で検出されたｚ軸に沿った線量を表す。中央のスライス厚Ｔの所で高い線量を示すが、実際のＸ線ファンビームＸＬＦＢにはｚ軸方向に広がりがあるため、線量プロフィールＤ（ｚ）にも図示の様な広がりがある。但し、実際上は、実効長１００ｍｍのペンシル型イオンチェンバ６２を使用することで、有効長１００ｍｍ分の測定を行っている。
【００６４】
図５に戻り、一方、ステップＳ１８ではＣＴ断層像上の各位置ｃ，ｐの近傍画素の各ＣＴ値からその標準偏差σ_c，σ_pを求める。図６（Ｂ）の挿入図（ａ）に標準偏差σを求める処理のイメージを示す。例えば注目画素Ｐ（ｘ，ｙ）を囲む５×５のブロック画素ＣＴ_i（ｉ＝１〜２５）につきＣＴ値の標準偏差σ（ｘ，ｙ）を（１０）式により求め、これを注目画素位置の標準偏差σ_pixelとする。
【００６５】
【数１０】

【００６６】
必要なら、注目画素Ｐ（ｘ，ｙ）を１画素づつ移動させた各ブロック領域における各標準偏差σ（ｘ，ｙ）を求め、これらを各注目画素の標準偏差σ_pixelとする。
【００６７】
好ましくは、このような処理を基準ファントム６０ｂの各位置ｃ．ｐの対応に行い、イオンチェンバ６２の線量検出領域ｃ，ｐにつきある程度の広がりがある場合には、夫々に求めた標準偏差σ_pixelの平均をとって各位置ｃ，ｐの標準偏差σ_c，σ_pとする。そして、これらを（１１）式により、各位置ｃ，ｐの規格化された線量ＣＴＤＩに対応付ける。
【００６８】
【数１１】

【００６９】
従って、上記（１１）式の関係より次の（１２）式の関係が得られ、ステップＳ１９では（１２）式の関係に従い、装置固有の係数パラメータＱ・Ｋ（本発明の係数パラメータｋに相当）を求める。
【００７０】
【数１２】

【００７１】
この（１２）式を満足するような一例の各値（σ，ＣＴＤＩ）は、ファントム中心部ｃにおける画像雑音と線量の実測値を夫々σ_c，ＣＴＤＩ_100,c 、ファントム周辺部ｐにおける４つの画像雑音と線量の実測値を夫々σ_p，ＣＴＤＩ_100,p とするときに、上記（２）式と同様の方法で加重平均された画像雑音σと、該（２）式と同じ方法で加重平均された平均線量値ｎＣＴＤＩｗとを（１２）式の関係で関係付けるような、装置固有の係数パラメータｋ（＝Ｑ・Ｋ）として求められる。
【００７２】
なお、線量ＣＴＤＩｗはＸ線の曝射エネルギーＣ（ｍＡｓ）に比例するため、予めｍＡｓ等の撮影条件を変えて上記同様の測定を行い、いくつかの代表的な各ｍＡｓ値に対応する複数の装置パラメータＱ・Ｋを求め、これらをｍＡｓ値の対応に表にしておいても良いし、又は各ｍＡｓ値に対応する複数の装置パラメータＱ・Ｋの平均を求めて、装置の係数パラメータとしても良い。また、必要なら、再構成関数を変えた場合の各係数パラメータを個々に求め、これらを再構成関数の対応に表にしておいても良い。また、Ｘ線ＣＴ装置の出荷時には、予め基準ファントム６０ｂ（６０ｈも同様）を使用して該装置の規格化された被曝線量ＣＴＤＩｗを求めておく。
【００７３】
図３は実施の形態によるＸ線ＣＴ撮影処理のフローチャート、図４は該処理のイメージ図であり、上記求めた装置固有の係数パラメータｋ（＝Ｑ・Ｋ）を使用して、被検体１００の実際のＣＴ撮影状況に忠実な被曝線量を出力する場合を示している。
【００７４】
ステップＳ３１では、技師が、スキャン計画情報（スキャンパラメータ，ＣＴ断層像の再構成パラメータ等）を設定する。ステップＳ３２では設定確認ボタン「ＣＯＮＦＩＲＭ」の入力を待ち、やがて、入力されると、ステップＳ３３では上記設定されたスキャンパラメータに従って被検体１００のスキャン・読取制御を行う。ステップＳ３４では被検体１００の投影データを収集・蓄積する。ステップＳ３５では所要撮影領域についてのスキャン完了か否かを判別し、完了でない場合はステップＳ３３に戻る。こうして、やがて、スキャン完了すると、ステップＳ３６では再構成パラメータに従って被検体１００のＣＴ断層像を再構成する。ステップＳ３７では被検体のＣＴ断層像を表示装置１３に表示する。図４（Ａ）に頭部Ｘ線ＣＴ断層像の表示イメージを示す。
【００７５】
ステップＳ３８では表示ＣＴ断層上で技師が被爆線量を知りたい部位（領域）を指定する。この指定方法には幾つか考えられる。例えば１スライス分のＣＴ断層像の全体を指定する。又は関心のある特定の臓器の部分のみを指定する。従って、特定の臓器，器官についての被曝線量を知ることが可能となる。
【００７６】
ところで、一般に特定の臓器を構成するような組織のＣＴ値は、その公称ＣＴ値を中心とするような所定の範囲内に収まるものと考えられる。例えば脂肪を抽出するためのＣＴ閾値範囲は−１５０＜ＣＴ数＜−５０とできる．そこで、予め、ＣＴ断層像から所望のＣＴ閾値範囲内に含まれるような組織（眼球，骨、脂肪等）を抽出しておくことが可能である。こうすれば、ある組織（例えば眼球）についてのＣＴ値の標準偏差σを求める場合に、予め不要な組織（骨等）のＣＴ値を演算から除外しておくことが可能となるため、ＣＴ値に関する画像雑音の演算の信頼性が大幅に向上する。
【００７７】
ステップＳ３９では、上記（１０）式に従い、注目画素位置近傍の各画素のＣＴ値から注目画素位置についてのＣＴ値の標準偏差σ_pixel）を求める。ステップＳ４０では（１３）式に従って標準偏差σ_pixelを対応する被曝線量値ＣＴＤＩ_pixelに変換する。
【００７８】
【数１３】

【００７９】
ステップＳ４２では画素値を被曝線量値ＣＴＤＩ_pixelとした被爆線量画像を作成する。そして、ステップＳ４２では被爆線量画像を表示部１３に表示する。図４（Ｂ）に頭部Ｘ線被曝線量画像の表示イメージを示す。このとき、好ましくは、元のＣＴイメージの輪郭線等を重ねて表示することで、臓器や組織に対する被曝線量の把握が容易となる。また、各線量値の大きさを色分けして表示すれば、被曝線量値の相違を認識し易い。また、上記（２），（３）式に従って、従来と同様に１スライス分の平均の被曝線量値を求めて、これを数値表示しても良い。更には、ｎスライス分の各平均線量値を加算することで、１撮影分の被曝線量を提供することも可能である。
【００８０】
なお、上記実施の形態では、ＣＴ線量測定用の基準ファントムとしてボディー部の材料がアクリル樹脂のものを使用したが、これに限らない。ボディー部の材料は均質であれば良く、例えば水ファントム等を使用しても良い。
【００８１】
また、上記実施の形態では、ＣＴ断層像の画像雑音をＣＴ値についての標準偏差σで評価したがこれに限らない。ＣＴ断層像の画像雑音をＣＴ値についての分散σ²で評価しても良い。
【００８２】
また、上記各実施の形態では、走査ガントリ３５上にＸ線管３１及び円弧状のＸ線検出器３３を搭載して被検体１００の回りに回転させるような所謂第３世代のＣＴガントリへの適用例を示したが、これに限らない。本発明は、図示しないが、被検体１００の回りに円形のＸ線検出器を固定すると共に、Ｘ線管３１を搭載した走査ガントリ３５を被検体体軸の周りに回転させて被検体１００の断層撮影を行うような所謂第４世代のＣＴガントリにも適用可能であることは明らかである。
【００８３】
また、上記本発明に好適なる実施の形態を述べたが、本発明思想を逸脱しない範囲内で各部の構成、制御、処理及びこれらの組み合わせの様々な変更が行えることは言うまでも無い。
【００８４】
【発明の効果】
以上述べた如く本発明によれば、被検体のサイズ等によらず、実際のＣＴ撮影に即した被曝線量を提供できるため、Ｘ線ＣＴ診断の安全性改善に寄与するところが極めて大きい。
【図面の簡単な説明】
【図１】本発明の原理を説明する図である。
【図２】実施の形態によるＸ線ＣＴ装置の要部構成図である。
【図３】実施の形態によるＸ線ＣＴ撮影処理のフローチャートである。
【図４】実施の形態によるＸ線ＣＴ撮影処理のイメージ図である。
【図５】実施の形態による装置パラメータ決定処理のフローチャートである。
【図６】実施の形態による装置パラメータ決定処理のイメージ図（１）である。
【図７】実施の形態による装置パラメータ決定処理のイメージ図（２）である。
【符号の説明】
１０ 操作コンソール部
２０ 撮影テーブル
３０ 走査ガントリ部
３５ ガントリ
４０ Ｘ線管
５０ コリメータ
６０ ＣＴ線量測定用ファントム
６１ ボディー
６２ ペンシル型イオンチェンバ
９０ Ｘ線検出器[0001]
BACKGROUND OF THE INVENTION
  The present invention relates to an X-ray CT apparatus andDose detectorMore specifically, an X-ray CT apparatus that includes an X-ray tube and an X-ray detector so as to face each other across the subject and reconstructs a CT tomographic image of the subject based on the detection signal of the X-ray detector as well asDose detectorAbout.
[0002]
[Prior art]
In the X-ray CT apparatus, it is recommended to display a CT exposure dose value (CTDIw: Computed Tomography Dose Index) of a subject according to the IEC standard. The outline of the ICE standard will be described below.
[0003]
First, the dose profile D (z) when the X-ray imaging system (gantry) is rotated once around the reference phantom is obtained by an ion chamber at each predetermined position (phantom central part c and peripheral part p) inserted into the phantom. taking measurement. The dose profile D (z) thus obtained represents the dose along the phantom center (z) axis detected by the ion chamber by one rotation scan, and shows a high dose at the slice thickness by X-ray. Based on the dose profile D (z) obtained in this way, a standardized CT dose value per unit length (CTDI: Computed Tomography Dose Index) is given by equation (1).
[0004]
[Expression 1]

[0005]
Where T: slice thickness
D (z): dose profile
Next, the standardized average dose nCTDIw per slice and 1 mAs is given by equation (2) by adding the CTDIw values of the central part c and the peripheral part p of the phantom at a predetermined ratio (weighted average). It is done.
[0006]
[Expression 2]

[0007]
Here, C: X-ray radiation dose (mAs)
CTDI_100,c: CT DI value limited to a thickness of 100 mm at the center c of the phantom
CTDI_100,p: Average of four CTDI values limited to a thickness of 100 mm at the phantom peripheral part p
The average dose CTDIw per slice used in the clinic (diagnosis) is given by equation (3) assuming that nCTDIw is multiplied by mAs.
[0008]
[Equation 3]

[0009]
Here, C: X-ray radiation dose (mAs)
Furthermore, the exposure dose DLP (Dose-Length Product) when the number of slices N is taken into consideration is expressed by equation (4).
[0010]
[Expression 4]

[0011]
Where i: each serial scan sequence
T: Thickness of each slice (cm)
N: Number of slices
Further, the exposure dose DLP in the case of the helical scan is expressed by equation (5).
[0012]
[Equation 5]

[0013]
Where T: nominal slice thickness (cm)
A: Tube current (mA)
t: Total exposure time
Under such a standard, conventionally, based on the CTDIw value measured in advance using each of the reference phantoms for the body and the head (CT dosimetry phantom), the exposure dose by the CT imaging of the subject is measured. Regardless of specimen size (body type, physique, etc.), it was displayed uniformly. That is, when the head of the subject is imaged, the exposure dose (CTDIw, DLP) of the subject is obtained and displayed based on the nCTDIw obtained beforehand using the head phantom, and the subject's exposure is displayed. When the body part was imaged, the exposure dose (CTDIw, DLP) of the subject was obtained and displayed based on the nCTDIw previously obtained using the body phantom.
[0014]
However, according to the conventional method, especially for infants and children whose size is smaller than the reference phantom, the exposure dose tends to be estimated lower than the actual dose, which is very dangerous.
[0015]
In this regard, the present applicant has already proposed an X-ray CT apparatus that determines the optimum X-ray irradiation condition (tube current) by detecting the subject size in advance (for example, Patent Documents 1 and 2).
[0016]
[Patent Document 1]
JP 2001-43993 A (paragraphs “0034” to “0036”, FIG. 6).
[0017]
[Patent Document 2]
Japanese Patent Application Laid-Open No. 2001-218761 (paragraphs “0093” to “0097”, FIGS. 6 and 7).
[0018]
[Problems to be solved by the invention]
However, in any of the above cases, in order to detect the subject size in advance, prior fluoroscopic imaging (scout scan) is necessary, and there has been a disadvantage that extra exposure is given to the subject.
[0019]
The present invention has been made in view of the above-described problems of the prior art, and an object of the present invention is to provide an X-ray CT apparatus, a dose detection method, and a program capable of providing an exposure dose faithful to actual CT imaging of a subject. There is to do.
[0020]
[Means for Solving the Problems]
  The above problem is solved by the configuration of FIG. That is, the X-ray CT apparatus of the present invention (1) includes an X-ray tube and an X-ray detector so as to face each other across the subject, and a CT tomogram of the subject based on the detection signal of the X-ray detector. In the X-ray CT apparatus for reconfiguringExtracting means for extracting a predetermined tissue portion from the CT tomographic image of the subject, and a CT tomographic image of the extracted tissueA dose calculating means for obtaining a dose at a corresponding position from the obtained image noise, and an output means for outputting information on the obtained dose.
[0021]
By the way, in a CT tomographic image of an imaging system in which the number of X-ray photons that pass through the subject is desired to be reduced as much as the X-ray CT apparatus, the time and space of the number of X-ray photons (photons) that pass through the subject. Fluctuation (photon noise) is the main cause of image noise. In this case, if the signal component of CT value is proportional to the dose (average value of the number of photons), the noise component (standard deviation) of CT value is It follows the Poisson process and is proportional to the square root of the dose.
[0022]
In the present invention (1), the dose calculation means 3 obtains the dose at the corresponding position (part) from the image noise obtained for the CT tomogram of the subject 100. It is possible to more faithfully provide the dose at the position, that is, the exposure dose at each position in consideration of the actual subject size and imaging conditions. Therefore, it is possible to effectively avoid a situation in which a small-sized subject such as an infant is overexposed. Moreover, since the average dose CTDIw per slice is not provided uniformly as in the past, the exposure dose for each position (tissue) constituting the CT tomogram can be provided. Diverse and useful exposure information can be obtained.
[0023]
Further, in the present invention (2), in the above-mentioned present invention (1), the dose calculating means 3 converts the image noise of the CT tomogram into the standard deviation σ of the CT value.₂When the predetermined coefficient parameter is k, the dose CTDI at the corresponding position₂CTDI₂= (Σ₂/ K)²  It is determined by the relationship.
[0024]
According to the present invention (2), the dose parameter CTDI of each part of the subject is obtained by obtaining the coefficient parameter k specific to the X-ray CT apparatus in advance.₂Is the standard deviation σ of the CT value of each corresponding part₂Can be easily obtained.
[0025]
In the present invention (3), in the present invention (1), the dose calculation means 3 obtains the dose at the corresponding position from the S / N value of the image obtained for the CT tomographic image of the subject 100.
[0026]
According to the present invention (3), the signal-to-noise (S / N) ratio of each part in the CT tomogram is proportional to the square root of the dose of the corresponding part. In addition, the dose at each position when the subject 100 is actually CT-scanned can be provided more faithfully. Conversely, the dose for obtaining the S / N value of a required image per CT tomogram can be effectively limited.
[0027]
  Further, in the present invention (4), in the present invention (1), the display means for displaying the CT tomogram of the subject on the screen, and the CT tomogram on the display screen.The organizationThe dose calculating means 3 calculates a dose for the designated area from the image noise obtained for the CT tomographic image of the designated area. Therefore, it is possible to easily detect the exposure dose to any tissue (brain, fat, liver, lung, etc.).
[0028]
  In the present invention (5), in the present invention (1),The extraction means is for extracting a portion of a tissue included in a CT threshold range separately input for a CT tomogram of a subject,Display means for displaying the extracted CT tomogram of the tissue on the screen in a manner distinguishable from the CT tomograms of other tissues is provided.
[0029]
According to the present invention (5), the engineer can easily extract a complex part of the corresponding tissue (brain, eyeball, bone, fat, liver, lung, etc.) simply by specifying the CT threshold range, and the exposure dose thereof Can be easily detected.
[0030]
In the present invention (6), in the present invention (1), the dose calculation means 3 obtains a dose corresponding to each pixel position from the image noise obtained for the correspondence of each pixel position with respect to the CT tomographic image of the subject, The output means 4 outputs an exposure dose image composed of the obtained dose values.
[0031]
According to the present invention (6), since the exposure dose obtained for the correspondence of each pixel position of the CT tomographic image can be obtained, the exposure situation of the subject can be accurately grasped (high resolution).
[0032]
In the present invention (7), in the present invention (6), the output means 4 outputs a dose image color-coded according to the dose value. Therefore, it is possible to easily recognize a precise exposure situation.
[0033]
Further, in the present invention (8), in the present invention (1), the dose calculation means 3 corresponds to each of the image noises obtained from a plurality of position correspondences including the central portion and the peripheral portion of the CT tomogram of the subject. The dose at the position to be obtained is obtained, and the average dose value for one slice is obtained by the weighted average.
[0034]
According to the present invention (8), it is possible to provide a dose value similar to the dose value CTDIw according to the above formulas (2) and (3), which is generally adopted from the past, and the commonality of diagnosis with the conventional system Easy to secure.
[0035]
In the present invention (9), in the present invention (8), the dose calculation means 3 obtains a dose value for one imaging by adding each average dose value obtained for each slice by the number of slices.
[0036]
According to the present invention (9), the exposure dose value for one image can be easily provided as in the prior art. Even in this case, the dose values for one radiographing can be provided more precisely by adding the average dose values obtained for each slice by the number of slices.
[0037]
Further, in the present invention (10), the image noise obtained for the CT tomographic image of the reference phantom 60 in the present invention (2) is converted into the standard deviation σ of the CT value.₁And the dose actually measured at the corresponding position in the phantom.₁The image noise σ₁The dose CTDI₁The coefficient parameter k for conversion into k = σ₁/ √CTDI₁The coefficient calculation means 2 calculated | required by these relationships is provided.
[0038]
According to the present invention (10), the image noise (standard deviation) σ obtained for the CT tomogram of the reference phantom 60₁And the dose actually measured in the reference phantom 60 are related in advance in the form of a coefficient parameter k unique to the apparatus, and thereafter, from image noise included in the CT tomographic image of the subject 100 The exposure dose according to actual CT imaging can be easily and accurately obtained. The exposure dose value obtained from the CT tomogram by the method of the present invention is also close to the actual measurement value measured by the reference phantom 60.
[0039]
In the present invention (11), in the present invention (10), the reference phantom 60 is parallel to the central axis of the cylindrical body made of a homogeneous material and the central axis at each predetermined position near the outer periphery thereof. A plurality of dose sensors 62 are provided. Therefore, the image noise of the CT tomographic image can be accurately obtained from the cylindrical body having a uniform material.
[0040]
In the present invention (12), in the present invention (11), the dose sensor 62 is formed of a pencil type ion chamber.
[0041]
Further, in the present invention (13), in the present invention (12), the dose calculation means 3 calculates each of the image noises obtained in correspondence with a plurality of positions including the central portion and the peripheral portion of the CT tomographic image of the phantom 60. The image noise σ is expressed when a weighted average is expressed as a standard deviation σ of CT values, and an average dose for one slice that is actually measured at each corresponding position in the phantom and is weighted averaged by the same method as described above is CTDI. Is obtained from the relationship k = σ / √CTDI.
[0042]
According to the present invention (13), the image noise σ weighted and averaged for the reference phantom 60 by the same method as the above expression (2) is converted into the average dose value CTDI weighted and averaged by the same method as the expression (2). By associating with a simple coefficient parameter k, the same exposure dose value CTDIw as conventionally employed can be easily obtained.
[0043]
  Further, the present invention (14)Dose detectorIsExtraction means for extracting a portion of a predetermined tissue for a CT tomogram of a subject, dose calculation means for obtaining a dose at a corresponding position from image noise obtained for the CT tomogram of the extracted tissue, and the obtained dose Output means for outputting the information.
  The dose detection method of the present invention includes an X-ray tube 40 and an X-ray detector 90 so as to face each other with the subject 100 interposed therebetween, and a CT tomographic image of the subject is obtained based on the detection signal of the X-ray detector. A method for detecting a dose of an X-ray CT apparatus to be reconstructed, wherein the image noise is determined based on an image noise obtained in advance for a CT tomographic image obtained by scanning a reference phantom 60 and a dose value measured in the phantom. And a step of obtaining a coefficient at a corresponding position from the image noise obtained with respect to a CT tomographic image obtained by scanning the object 100 using the coefficient parameter. is there.
[0044]
  Preferably,The image noise obtained for the CT tomogram of the reference phantom 60 is σ.₁, And dose values actually measured at corresponding positions of the phantom₁The coefficient parameter k is k = σ₁/ √CTDI₁And the image noise obtained for the CT tomographic image of the subject 100 as σ.₂The dose value CTDI at the corresponding position₂CTDI₂= (Σ₂/ K)²And a step of obtaining by the relationship.
[0045]
  The program of the present invention is stored in a computer.the aboveThis is a computer-executable program for executing the dose detection method. Such a program can be recorded on a recording medium such as a CD-ROM or provided by online communication via a communication line.
[0046]
DETAILED DESCRIPTION OF THE INVENTION
DESCRIPTION OF EXEMPLARY EMBODIMENTS Hereinafter, preferred embodiments of the invention will be described in detail with reference to the accompanying drawings. Note that the same reference numerals denote the same or corresponding parts throughout the drawings.
[0047]
FIG. 2 is a configuration diagram of a main part of an X-ray CT apparatus according to the embodiment. The apparatus includes an X-ray fan beam XLFB that performs axial / helical scanning / reading of the subject 100, and the subject 100. An imaging table 20 that moves in the direction of the body axis CLb, and a remote operation console unit 10 that controls the scanning gantry unit 30 and the imaging table 20 and performs various settings and operations by an X-ray imaging engineer. Is provided.
[0048]
In the scanning gantry 30, 40 is a rotary anode type X-ray tube, 40 A is an X-ray tube controller for controlling the tube voltage kV, tube current mA, etc. of the X-ray tube, and 50 is a collimator for limiting the slice thickness of the X-ray. , 50A is a collimator control unit, 90 is an X-ray detector (multi-detector) in which a large number (n = 1000) of X-ray detection elements arranged in the channel CH direction are arranged in, for example, two rows L1 and L2 in the body axis CLb direction. ) 91 is the projection data g of the subject 100 based on the detection signal of the X-ray detector 90.₁(X, θ), g₂A data collection unit (DAS) 35 that generates and collects (X, θ), 35 is a gantry that rotatably supports each device related to the X-ray imaging system around the body axis CLb, and 35A is a rotation control of the gantry 35. Part.
[0049]
In the operation console unit 10, 11 is a central processing unit that performs main control and processing (scan control, reconstruction of CT tomographic images, various processes related to X-ray exposure display) of the X-ray CT apparatus, 11a is its CPU, and 11b is Main memory (MM) composed of RAM, ROM, etc. used by the CPU 11a, 12 is a command and data input device including a keyboard, mouse, etc. 13 is a display of scan plan information, CT tomogram, and various information regarding X-ray exposure A display device (CRT) 14 for performing the control interface for exchanging various control signals CS and monitor signals MS between the CPU 11a, the scanning gantry unit 30, the imaging table 20, and the like, and 15 a projection from the data collecting unit 91. A data collection buffer for temporarily accumulating data, 16 is finally used for scan (projection) data and CT tomographic data. While product and storing a to have secondary storage device stores the data file or the like for various application programs and various arithmetic / correction required for the operation of the X-ray CT apparatus (hard disk drive, etc.).
[0050]
An outline of the basic operation of subject CT imaging is as follows. The X-ray fan beam XLFB from the X-ray tube 40 passes through the subject 100 and enters the detection rows L1 and L2 of the X-ray detector 90 all at once. The data collection unit 91 outputs projection data g corresponding to each detection output of the X-ray detector 90.₁(X, θ), g₂(X, θ) are generated and stored in the data collection buffer 15. Here, X represents the detection channels 1 to n of the X-ray detector 90, and θ represents the view angle around the body axis CLb. Further, projection similar to the above is performed at each view angle θ where the gantry 35 is slightly rotated, and thus projection data for one rotation of the gantry is collected and accumulated. At the same time, the imaging table 20 is moved intermittently / continuously in the direction of the body axis CLb in accordance with the axial / helical scan method, and thus all projection data for the required imaging region of the subject 100 is collected and accumulated. Then, the CPU 11a reconstructs a CT tomogram of the subject 100 based on the obtained projection data after displaying all the scans or following (in parallel with) the scan execution, and displays this on the display device 13. . In this embodiment, the CT tomographic image data can be used to determine the exposure dose of the subject and display it on the screen. Details will be described below.
[0051]
FIG. 5 is a flowchart of the apparatus parameter determination process according to the embodiment, and FIGS. 6 and 7 are image diagrams (1) and (2) of the process. The image noise included in the CT tomogram when the reference phantom is read in advance. The process which determines the parameter k intrinsic | native to the apparatus for calculating | requiring the corresponding exposure dose from is shown.
[0052]
In FIG. 5, in step S <b> 11, a reference dose for CT dose measurement is attached to the tip of the bed 20. The mounting image is shown in FIG. FIG. 7A shows an example of a CT dose measurement phantom (reference phantom). 7A, the body 61 of the torso phantom 60b is made of a cylindrical body of a homogeneous acrylic resin (PPMA: polymethylmethacrylate), and its diameter φ is a typical torso size of the subject. As representative, φ = 32 cm, and the thickness is 14 cm or more. Further, the phantom 60b includes a pencil-type ion having an execution length of 10 cm at a central portion c parallel to the z-axis of the body 61 (corresponding to the subject body axis CLb) and at four locations p 1 cm inside from the outer periphery. A chamber (pencil ionisation chamber) 62 is inserted. On the other hand, the diameter φ of the head phantom 60h is φ = 16 cm, which represents a typical head size of the subject. Other configurations are the same as those of the phantom 60b.
[0053]
Returning to FIG. 5, in step S12, the engineer sets scan parameters (tube voltage kV, tube current mA, X-ray exposure time s, etc.) for a reference phantom (for example, 60b). In step S13, the reference phantom 60b is scanned (for one slice), and in step S14, the projection data is collected. In step S15, a CT tomographic image of the reference phantom 60b is reconstructed. Each pixel data of the obtained CT tomographic image is composed of a CT value of a homogeneous acrylic resin, and is expressed by equation (6).
[0054]
[Formula 6]

[0055]
Where μ_w: Water line attenuation coefficient
μ_m: Linear attenuation coefficient of the substance (acrylic resin in this example)
Since the above-mentioned uniform acrylic resin was scanned, the obtained CT value should be uniform. However, actually, due to temporal and spatial fluctuations (photon noise) of the number of X-ray photons transmitted through the phantom 60b, fluctuations due to the photon noise also occur in the reconstructed CT value, and this is an image of a CT tomogram. Appears as noise.
[0056]
FIG. 6B shows a CT image including image noise of the reference phantom 60b. As in an X-ray CT apparatus, in a CT tomographic image of an imaging system in which the number of X-ray photons that pass through the subject 100 (exposure dose) is desired to be reduced as much as possible, the time is the number of X-ray photons (photons) that pass through the subject. And spatial fluctuations (photon noise) are the main causes of image noise. In this case, if the signal component of the CT value is proportional to the dose (average value of the number of photons), the noise component of the CT value (Standard deviation) is proportional to the square root of the dose according to the Poisson process. That is, the signal-to-noise ratio (S / N) of the CT tomogram is expressed by equation (7).
[0057]
[Expression 7]

[0058]
Where a: constant
Dose: Dose
Further, there is a relationship of the formula (8) between the image noise N and the dose Dose.
[0059]
[Equation 8]

[0060]
Where k is a constant
Now, the dose at each position c, p of the reference phantom 60b is Dose._c, Dose_pAnd the image noise N at each of the positions c and p is a standard deviation σ with respect to the CT value._c, Σ_pImage noise σ_c, Σ_pAnd dose Dose_c, Dose_pThere is a relationship of the formula (9).
[0061]
[Equation 9]

[0062]
Where Q: coefficient determined by device configuration
K: Coefficient determined by the reconstruction function
Dose CTDI at each position c and p in the reference phantom 60b_100,c, CTDI_100,Since p can be measured by the ion chamber 62 at each position, measurement data of each ion chamber 62 is collected at step S16 in FIG. 5, and dose values CTDI at the respective positions c and p are collected at step S17._100,c, CTDI_100,Find p.
[0063]
FIG. 7B shows an example of the dose profile D (z). The dose profile D (z) represents the dose along the z-axis detected by the ion chamber 62 by one rotation scan. Although a high dose is shown at the center slice thickness T, since the actual X-ray fan beam XLFB has a spread in the z-axis direction, the dose profile D (z) also has a spread as shown in the figure. However, in practice, the pencil-type ion chamber 62 having an effective length of 100 mm is used to measure the effective length of 100 mm.
[0064]
Returning to FIG. 5, on the other hand, in step S18, the standard deviation σ is calculated from the CT values of the neighboring pixels at the positions c and p on the CT tomographic image._c, Σ_pAsk for. Inset (a) in FIG. 6B shows an image of processing for obtaining the standard deviation σ. For example, a 5 × 5 block pixel CT surrounding the pixel of interest P (x, y)_iThe standard deviation σ (x, y) of the CT value for (i = 1 to 25) is obtained by the equation (10), and this is calculated as the standard deviation σ of the target pixel position._pixelAnd
[0065]
[Expression 10]

[0066]
If necessary, each standard deviation σ (x, y) in each block area where the pixel of interest P (x, y) is moved one pixel at a time is obtained, and these are obtained as the standard deviation σ of each pixel of interest._pixelAnd
[0067]
Preferably, such processing is performed at each position c. Of the reference phantom 60b. When there is a certain extent in the dose detection areas c and p of the ion chamber 62 in correspondence with p, the standard deviation σ obtained respectively is_pixelThe standard deviation σ of each position c and p is averaged_c, Σ_pAnd And these are matched with the standardized dose CTDI of each position c and p by (11) Formula.
[0068]
[Expression 11]

[0069]
Therefore, the relationship of the following equation (12) is obtained from the relationship of the above equation (11), and in step S19, the coefficient parameter Q · K (corresponding to the coefficient parameter k of the present invention) according to the relationship of the equation (12). )
[0070]
[Expression 12]

[0071]
Each example value (σ, CTDI) satisfying the equation (12) is the measured value of image noise and dose at the phantom center c._c, CTDI_{100, c}, The measured values of the four image noises and doses in the phantom peripheral part p are respectively σ_p, CTDI_{100, p}When the image noise σ weighted and averaged by the same method as the above equation (2) and the average dose value nCTDIw weighted and averaged by the same method as the above equation (2) are expressed by the relationship of the equation (12): It is obtained as a coefficient parameter k (= Q · K) specific to the device that is related.
[0072]
Since the dose CTDIw is proportional to the X-ray exposure energy C (mAs), the same measurement is performed in advance by changing imaging conditions such as mAs, and a plurality of representative mAs values corresponding to several mAs values are obtained. The apparatus parameters Q · K may be obtained and tabulated for the correspondence of mAs values, or the average of a plurality of apparatus parameters Q · K corresponding to each mAs value may be obtained and used as the coefficient parameter of the apparatus. good. If necessary, each coefficient parameter when the reconstruction function is changed may be obtained individually and tabulated in correspondence with the reconstruction function. In addition, when the X-ray CT apparatus is shipped, the standardized dose phantom CTDIw of the apparatus is obtained in advance using the reference phantom 60b (same for 60h).
[0073]
FIG. 3 is a flowchart of the X-ray CT imaging process according to the embodiment, and FIG. 4 is an image diagram of the process. Using the coefficient parameter k (= Q · K) unique to the apparatus, the actual measurement of the subject 100 is performed. The case where the exposure dose faithful to the CT imaging situation is output is shown.
[0074]
In step S31, the engineer sets scan plan information (scan parameters, CT tomographic reconstruction parameters, etc.). In step S32, input of the setting confirmation button “CONFIRM” is awaited, and when it is input, scan / reading control of the subject 100 is performed in step S33 according to the set scan parameter. In step S34, projection data of the subject 100 is collected and accumulated. In step S35, it is determined whether or not the scan for the required imaging area is complete. If not, the process returns to step S33. Thus, when the scanning is completed, the CT tomographic image of the subject 100 is reconstructed according to the reconstruction parameter in step S36. In step S37, a CT tomogram of the subject is displayed on the display device 13. FIG. 4A shows a display image of a head X-ray CT tomogram.
[0075]
In step S38, the engineer designates a part (region) on the displayed CT tomography to know the exposure dose. There are several possible ways to specify this. For example, the whole CT tomographic image for one slice is designated. Or specify only the specific organ part of interest. Therefore, it becomes possible to know the exposure dose for a specific organ or organ.
[0076]
By the way, it is generally considered that the CT value of a tissue constituting a specific organ falls within a predetermined range centered on the nominal CT value. For example, the CT threshold range for extracting fat can be −150 <CT number <−50. Therefore, it is possible to extract in advance a tissue (eyeball, bone, fat, etc.) included in a desired CT threshold range from the CT tomogram. In this way, when obtaining the standard deviation σ of the CT value for a certain tissue (for example, the eyeball), the CT value of an unnecessary tissue (such as a bone) can be excluded from the calculation in advance. The reliability of the calculation of image noise is greatly improved.
[0077]
In step S39, the standard deviation σ of the CT value for the target pixel position from the CT value of each pixel in the vicinity of the target pixel position according to the above equation (10)._pixel) In step S40, the standard deviation σ according to the equation (13)._pixelDose value CTDI corresponding to_pixelConvert to
[0078]
[Formula 13]

[0079]
In step S42, the pixel value is converted into the exposure dose value CTDI._pixelCreate an exposure dose image. In step S42, the exposure dose image is displayed on the display unit 13. FIG. 4B shows a display image of the head X-ray exposure image. At this time, preferably, the outline of the original CT image is displayed in an overlapping manner, so that it becomes easy to grasp the exposure dose to the organ or tissue. Moreover, if the magnitude of each dose value is displayed in different colors, it is easy to recognize the difference in the dose value. Further, according to the above formulas (2) and (3), an average exposure dose value for one slice may be obtained in the same manner as in the past, and this may be displayed numerically. Furthermore, it is also possible to provide an exposure dose for one imaging by adding the average dose values for n slices.
[0080]
In the above embodiment, the material of the body part is made of acrylic resin as the reference phantom for CT dose measurement, but the invention is not limited to this. The material of the body part may be homogeneous, and for example, a water phantom or the like may be used.
[0081]
Moreover, in the said embodiment, although the image noise of CT tomogram was evaluated by the standard deviation (sigma) about CT value, it is not restricted to this. The image noise of CT tomogram is expressed as variance σ with respect to CT value.²You may evaluate with.
[0082]
In each of the above-described embodiments, an X-ray tube 31 and an arc-shaped X-ray detector 33 are mounted on the scanning gantry 35 and are rotated around the subject 100 so-called third generation CT gantry. Although an application example was shown, it is not limited to this. Although not shown in the drawings, the present invention fixes a circular X-ray detector around the subject 100 and rotates the scanning gantry 35 on which the X-ray tube 31 is mounted around the subject body axis. It is obvious that the present invention can also be applied to a so-called fourth generation CT gantry that performs tomography.
[0083]
Moreover, although the preferred embodiment of the present invention has been described, it goes without saying that various changes in the configuration, control, processing, and combination of each part can be made without departing from the spirit of the present invention.
[0084]
【The invention's effect】
As described above, according to the present invention, it is possible to provide an exposure dose suitable for the actual CT imaging regardless of the size of the subject, so that it greatly contributes to improving the safety of X-ray CT diagnosis.
[Brief description of the drawings]
FIG. 1 is a diagram illustrating the principle of the present invention.
FIG. 2 is a configuration diagram of a main part of an X-ray CT apparatus according to an embodiment.
FIG. 3 is a flowchart of X-ray CT imaging processing according to the embodiment.
FIG. 4 is an image diagram of X-ray CT imaging processing according to the embodiment.
FIG. 5 is a flowchart of apparatus parameter determination processing according to the embodiment.
FIG. 6 is an image diagram (1) of apparatus parameter determination processing according to the embodiment.
FIG. 7 is an image diagram (2) of the apparatus parameter determination process according to the embodiment.
[Explanation of symbols]
10 Operation console section
20 Shooting table
30 Scanning gantry section
35 Gantry
40 X-ray tube
50 collimator
60 CT phantom for dosimetry
61 body
62 Pencil type ion chamber
90 X-ray detector

Claims (14)

In an X-ray CT apparatus that includes an X-ray tube and an X-ray detector so as to face each other with a subject interposed therebetween, and reconstructs a CT tomogram of the subject based on a detection signal of the X-ray detector,
Extraction means for extracting a predetermined tissue portion for the CT tomogram of the subject;
Dose calculating means for obtaining a dose at a corresponding position from image noise obtained for the CT tomogram of the extracted tissue ;
An X-ray CT apparatus comprising: output means for outputting information on the obtained dose. The dose calculation means represents the image noise of the CT tomogram by the standard deviation σ ₂ of the CT value, and when the predetermined coefficient parameter is k, the dose CTDI ₂ at the corresponding position is
CTDI ₂ = (σ ₂ / k) ²
The X-ray CT apparatus according to claim 1, wherein the X-ray CT apparatus is obtained according to the relationship: The X-ray CT apparatus according to claim 1, wherein the dose calculation means obtains a dose at a corresponding position from an S / N value of an image obtained for a CT tomogram of the subject. Display means for displaying a CT tomogram of the subject on the screen;
Designating means for designating the tissue of the CT tomogram on the display screen,
The X-ray CT apparatus according to claim 1, wherein the dose calculation means obtains a dose for the designated region from image noise obtained for the CT tomographic image of the designated region. The extraction means is for extracting a portion of a tissue included in a CT threshold range separately input for a CT tomogram of a subject,
The X-ray CT apparatus according to claim 1, further comprising display means for displaying the extracted CT tomogram of the tissue on a screen in a manner distinguishable from CT tomograms of other tissues. The dose calculation means obtains the dose corresponding to each pixel position from the image noise obtained for the correspondence of each pixel position for the CT tomographic image of the subject,
The X-ray CT apparatus according to claim 1, wherein the output unit outputs a dose image including the obtained dose values. The X-ray CT apparatus according to claim 6, wherein the output unit outputs a dose image color-coded according to a dose value. The dose calculation means calculates a dose at each corresponding position from each image noise determined for a plurality of positions including the central portion and the peripheral portion of the CT tomographic image of the subject, and averages the dose for one slice by weighted average of these The X-ray CT apparatus according to claim 1, wherein a value is obtained. The X-ray CT apparatus according to claim 8, wherein the dose calculation means obtains a dose value for one imaging by adding each average dose value obtained for each slice by the number of slices. When the image noise obtained for the CT tomographic image of the reference phantom is represented by the standard deviation σ ₁ of the CT value and the dose actually measured at the corresponding position in the phantom is CTDI ₁ , the image noise σ ₁ is The coefficient parameter k for converting to the dose CTDI ₁ is
k = σ ₁ / √CTDI ₁
The X-ray CT apparatus according to claim 2, further comprising a coefficient calculation unit that is obtained based on the relationship of: The X-ray phantom according to claim 10, wherein the reference phantom includes a plurality of dose sensors in parallel with the central axis of a cylindrical body made of a homogeneous material and the central axis at each predetermined position near the outer periphery thereof. Line CT device. The X-ray CT apparatus according to claim 11, wherein the dose sensor is a pencil-type ion chamber. The dose calculation means represents a weighted average of each image noise obtained for a plurality of positions including a central portion and a peripheral portion of the CT tomographic image of the phantom by a standard deviation σ of CT values, and each correspondence in the phantom When an average dose of one slice measured at a position and weighted averaged in the same manner as described above is defined as CTDI, a coefficient parameter k for converting the image noise σ into the dose CTDI is set to k = σ / √CTDI.
The X-ray CT apparatus according to claim 12, wherein the X-ray CT apparatus is obtained according to the relationship: Extraction means for extracting a predetermined tissue portion from a CT tomogram of a subject;
Dose calculating means for obtaining a dose at a corresponding position from image noise obtained for the CT tomogram of the extracted tissue;
An output means for outputting information on the determined dose.

Images