知財求人 - 知財ポータルサイト「IP Force」
▶ ジーイー・プレシジョン・ヘルスケア・エルエルシーの特許一覧
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】特許公報(B1)
(11)【特許番号】
(24)【登録日】2025-09-26
(45)【発行日】2025-10-06
(54)【発明の名称】医用システムおよび記憶媒体
(51)【国際特許分類】
A61B 6/58 20240101AFI20250929BHJP
A61B 6/03 20060101ALI20250929BHJP
【ＦＩ】
A61B6/58 500Z
A61B6/03 573
【請求項の数】
10
(21)【出願番号】P 2024063670
(22)【出願日】2024-04-10
【審査請求日】2024-04-15
【前置審査】
(73)【特許権者】
【識別番号】319011672
【氏名又は名称】ジーイー・プレシジョン・ヘルスケア・エルエルシー
(74)【代理人】
【識別番号】100105588
【弁理士】
【氏名又は名称】小倉 博
(74)【代理人】
【識別番号】100151286
【弁理士】
【氏名又は名称】澤木 亮一
(72)【発明者】
【氏名】星野 耕一
(72)【発明者】
【氏名】今井 靖浩
【審査官】蔵田 真彦
(56)【参考文献】
【文献】特開２０２３－１１９３７５（ＪＰ，Ａ）
【文献】特開２０１４－１２８４５６（ＪＰ，Ａ）
【文献】特開２０２４－０１３１３８（ＪＰ，Ａ）
【文献】特開２００１－０７０２９７（ＪＰ，Ａ）
【文献】特開昭５４－０２２１９２（ＪＰ，Ａ）
【文献】特開２０２２－１０３６１５（ＪＰ，Ａ）
【文献】国際公開第２０１５／０６４４４６（ＷＯ，Ａ１）
【文献】米国特許出願公開第２０１７／０３１１９１７（ＵＳ，Ａ１）
(58)【調査した分野】(Int.Cl.，ＤＢ名)
Ａ６１Ｂ ６／００－６／５８
(57)【特許請求の範囲】
【請求項1】
Ｘ線発生装置と、前記Ｘ線発生装置から照射されたＸ線を検出する検出器と、
１つ又は複数のプロセッサであって、
所定のキャリブレーションスキャンが実行されるたびに、前記検出器で検出されたＸ線のデータに基づいて、前記Ｘ線発生装置から照射されたＸ線の特性値を計算すること、
前記特性値の経時的な変化に基づいて、アラートを出力するかどうか判断すること
を実行する１つ又は複数のプロセッサと
を含む、医用システムであって、
前記医用システムは、管電圧を第１の管電圧と第２の管電圧との間で交互に切り替えながら前記所定のキャリブレーションスキャンを定期的に実行し、
前記特性値は、第１の基準物質の吸収係数と、第２の基準物質の吸収係数とを含む、医用システム。
【請求項2】
前記１つ又は複数のプロセッサは、前記所定のキャリブレーションスキャンが実行されるたびに、第１の基準物質の吸収係数と、第２の基準物質の吸収係数を計算する、請求項１に記載の医用システム。
【請求項3】
前記１つ又は複数のプロセッサは、アラートを出力するかどうかを判断するための基準となる第１の基準吸収係数を計算すること、
前記第１の基準吸収係数と前記第１の基準物質の吸収係数との第１の差を計算すること、および
前記第１の差と第１の閾値とを比較し、その比較結果に基づいて、アラートを出力するかどうか判断すること
を実行する、請求項２に記載の医用システム。
【請求項4】
前記１つ又は複数のプロセッサは、アラートを出力するかどうかを判断するための基準となる第２の基準吸収係数を計算すること、
前記第２の基準吸収係数と前記第２の基準物質の吸収係数との第２の差を計算すること、および
前記第２の差と第２の閾値とを比較し、その比較結果に基づいて、アラートを出力するかどうか判断すること
を実行する、請求項３に記載の医用システム。
【請求項5】
前記１つ又は複数のプロセッサは、前記所定のキャリブレーションスキャンを実行するたびに、過去に実行された前記所定のキャリブレーションスキャンにより得られた第１の基準物質の吸収係数に基づいて、前記第１の基準吸収係数を計算すること、又は
過去に実行された前記所定のキャリブレーションスキャンにより得られた第２の基準物質の吸収係数に基づいて、前記第２の基準吸収係数を計算すること
を実行する、請求項４に記載の医用システム。
【請求項6】
前記第１の基準物質は水であり、前記第２の基準物質はヨードである、又は前記第１の基準物質は水であり、前記第２の基準物質はカルシウムである、請求項５に記載の医用システム。
【請求項7】
前記アラートが出力された場合、前記医用システムは、吸収係数の時間変化を表すデータをバックオフィスに送る、請求項１に記載の医用システム。
【請求項8】
前記第１の基準物質の吸収係数および前記第２の基準物質の吸収係数は、前記検出器の各列のチャネルごとに計算される、請求項１に記載の医用システム。
【請求項9】
前記医用システムは、複数のキャリブレーションスキャンを実行し、前記複数のキャリブレーションスキャンのうちの１つのキャリブレーションスキャンが、前記所定のキャリブレーションスキャンとして実行される、請求項１に記載の医用システム。
【請求項10】
Ｘ線発生装置と、前記Ｘ線発生装置から照射されたＸ線を検出する検出器と、
１つ又は複数のプロセッサと
を含む医用システムに含まれる又は前記医用システムと通信可能である非一時的なコンピュータ可読記憶媒体であって、
前記記憶媒体に記憶された命令が、前記１つ又は複数のプロセッサによって実行されると、前記１つ又は複数のプロセッサに、
所定のキャリブレーションスキャンが実行されるたびに、前記検出器で検出されたＸ線のデータに基づいて、前記Ｘ線発生装置から照射されたＸ線の特性値を計算すること、
前記特性値の経時的な変化に基づいて、アラートを出力するかどうか判断すること
を実行させ、
前記医用システムは、管電圧を第１の管電圧と第２の管電圧との間で交互に切り替えながら前記所定のキャリブレーションスキャンを定期的に実行し、
前記特性値は、第１の基準物質の吸収係数と、第２の基準物質の吸収係数とを含む、非一時的なコンピュータ可読記憶媒体。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、Ｘ線を照射する医用システム、および当該医用システムを制御するための命令が記録された記憶媒体に関する。
【背景技術】
【0002】
被検体を非侵襲的に撮影する医用システムとしてＣＴシステムが知られている。ＣＴシステムは、短いスキャン時間で被検体の断層画像を取得することができるので、病院などの医療施設に普及している。
【0003】
ＣＴシステムは、Ｘ線管の陰極－陽極管に所定の電圧を印加し、Ｘ線を発生させる。発生したＸ線は被検体を透過して検出器で検出される。ＣＴシステムは、検出器で検出されたデータに基づいて、被検体のＣＴ画像を再構成している。
【先行技術文献】
【特許文献】
【0004】
【文献】特許第６０３１６１８号公報
【発明の概要】
【発明が解決しようとする課題】
【0005】
ＣＴシステムの撮影技術として、ＳＥＣＴ（Single Energy CT）が知られている。ＳＥＣＴは、Ｘ線管の陰極－陽極管に所定の電圧（例えば、１２０ｋＶ）を印加し、Ｘ線を発生させ、被検体のＣＴ画像を得る方法である。しかし、ＳＥＣＴでは、異なる物質であってもＣＴ値が近い値になることがあり、異なる物質の同定が難しい場合がある。
【0006】
そこで、ＤＥＣＴ（Dual Energy CT）の技術が研究、開発されている。ＤＥＣＴは、異なるエネルギー領域のＸ線を利用して物質の弁別を行うことができる技術であり、臨床現場での診断に有益な画像を取得することができ、広く普及し始めている。ＤＥＣＴの技術では、Ｘ線管の管電圧を、低い管電圧と高い管電圧との間で切り替えるｋＶスイッチング技術が知られている。
【0007】
一方、ＤＥＣＴのｋＶスイッチング法による物質の弁別精度は、Ｘ線発生装置（例えば、Ｘ線管、ジェネレータなど）の経年劣化によって低下する。したがって、物質の弁別精度を低下させないようにするために、ＣＴシステムのユーザは、Ｘ線発生装置などの経年劣化をできるだけ早急に知ることが重要となる。しかし、Ｘ線発生装置の劣化を検出することは行われておらず、ＣＴシステムのユーザは、Ｘ線発生装置の劣化がＣＴ画像にノイズやアーチファクトとして現れることによって、Ｘ線発生装置の劣化に気が付く。したがって、ＣＴ画像の画質が低下する前に（ＣＴ画像にノイズやアーチファクトが現れる前に）、Ｘ線発生装置に経年劣化が発生している恐れがあることをユーザに知らせる技術が要求されている。
【課題を解決するための手段】
【0008】
本発明の第１の観点は、Ｘ線発生装置と、
前記Ｘ線発生装置から照射されたＸ線を検出する検出器と、
１つ又は複数のプロセッサであって、
所定のキャリブレーションスキャンが実行されるたびに、前記検出器で検出されたＸ線のデータに基づいて、前記Ｘ線発生装置から照射されたＸ線の特性値を計算すること、
前記特性値の経時的な変化に基づいて、アラートを出力するかどうか判断すること
を実行する１つ又は複数のプロセッサと
を含む、医用システムである。
前記Ｘ線発生装置から照射されたＸ線を検出する検出器と、
１つ又は複数のプロセッサであって、
所定のキャリブレーションスキャンが実行されるたびに、前記検出器で検出されたＸ線のデータに基づいて、前記Ｘ線発生装置から照射されたＸ線の特性値を計算すること、
前記特性値の経時的な変化に基づいて、アラートを出力するかどうか判断すること
を実行する１つ又は複数のプロセッサと
を含む、医用システムである。
【0009】
また、本発明の第２の観点は、Ｘ線発生装置と、
前記Ｘ線発生装置から照射されたＸ線を検出する検出器と、
１つ又は複数のプロセッサと
を含む医用システムに含まれる又は前記医用システムと通信可能である非一時的なコンピュータ可読記憶媒体であって、
前記記憶媒体に記憶された命令が前記１つ又は複数のプロセッサによって実行されると、前記１つ又は複数のプロセッサに、
所定のキャリブレーションスキャンが実行されるたびに、前記検出器で検出されたＸ線のデータに基づいて、前記Ｘ線発生装置から照射されたＸ線の特性値を計算すること、
前記特性値の経時的な変化に基づいて、アラートを出力するかどうか判断すること
を実行させる、非一時的なコンピュータ可読記憶媒体である。
前記Ｘ線発生装置から照射されたＸ線を検出する検出器と、
１つ又は複数のプロセッサと
を含む医用システムに含まれる又は前記医用システムと通信可能である非一時的なコンピュータ可読記憶媒体であって、
前記記憶媒体に記憶された命令が前記１つ又は複数のプロセッサによって実行されると、前記１つ又は複数のプロセッサに、
所定のキャリブレーションスキャンが実行されるたびに、前記検出器で検出されたＸ線のデータに基づいて、前記Ｘ線発生装置から照射されたＸ線の特性値を計算すること、
前記特性値の経時的な変化に基づいて、アラートを出力するかどうか判断すること
を実行させる、非一時的なコンピュータ可読記憶媒体である。
【発明の効果】
【0010】
本発明では、前記検出器で検出されたＸ線の情報を含むデータに基づいて、前記Ｘ線発生装置から照射されたＸ線の特性値を計算し、前記特性値の経時的な変化に基づいてアラートを出力するかどうか判断される。したがって、特性値の経時的な変化が大きい場合にアラートが出力されるようにすることで、ＣＴ画像の画質が低下する前に、Ｘ線発生装置に劣化の徴候があることをユーザに知らせることができる。
【図面の簡単な説明】
【0011】
【図1】第１の実施形態におけるＣＴシステム１０のブロック図である。
【図2】各診療日に実行されるフローを示す図である。
【図3】キャリブレーションの説明図である。
【図4】診療日１において取得されたキャリブレーションデータの説明図である。
【図5】吸収係数データｇ１１の拡大図である。
【図6】吸収係数データｈ１１の拡大図である。
【図7】診療日２に得られた第１の基準物質の吸収係数データｇ２１を示す図である。
【図8】診療日２に得られた第２の基準物質の吸収係数データｈ２１を示す図である。
【図9】ステップＳＴ２０のフローの一例を示す図である。
【図11】チャネル２に対してステップＳＴ２０１およびＳＴ２０２を実行するときの説明図である。
【図12】チャネルｊに対してステップＳＴ２０１およびＳＴ２０２を実行するときの説明図である。
【図13】チャネルｎに対してステップＳＴ２０１およびＳＴ２０２を実行するときの説明図である。
【図14】第２の基準物質の吸収係数データに対してステップＳＴ２０を実行するときの説明図である。
【図15】診療日１～３に得られた第１の基準物質の吸収係数データｇ１１～ｇ３１を示す図である。
【図16】診療日１～３に得られた第２の基準物質の吸収係数データｈ１１～ｈ３１を示す図である。
【図17】ステップＳＴ２０の説明図である。
【図18】第２の基準物質の吸収係数データに対してステップＳＴ２０１～ＳＴ２０４を実行するときの説明図である。
【図19】診療日ｉに得られた第１の基準物質の吸収係数データｇｉ１を示す図である。
【図20】診療日ｉに得られた第２の基準物質の吸収係数データｈｉ１を示す図である。
【図21】ステップＳＴ２０の説明図である。
【図22】ＣＴシステムの稼働再開後のフローの説明図である。
【図23】診療日ｐに得られた第１の基準物質の吸収係数データｇｐ１を示す図である。
【図24】診療日ｐに得られた第２の基準物質の吸収係数データｈｐ１を示す図である。
【図25】診療日ｑに得られた第１の基準物質の吸収係数データｇｑ１を示す図である。
【図26】診療日ｑに得られた第２の基準物質の吸収係数データｈｑ１を示す図である。
【図27】ステップＳＴ２０の説明図である。
【図28】第２の基準物質の吸収係数データに対してステップＳＴ２０１～ＳＴ２０４を実行するときの説明図である。
【図29】複数列の検出器が使用された場合のステップＳＴ２０のフローの一例を示す図である。
【図30】診療日ごとに基準吸収係数を更新する場合のステップＳＴ２０のフローの説明図である。
【図31】診療日１～診療日ｉで取得された第１の基準物質の吸収係数データｇ１１～ｇｉ１を示す図である。
【図32】第２の実施形態で実施されるフローの説明図である。
【図33】Ｘ線スペクトルの説明図である。
【図34】診療日２で生成された低ｋＶに対応するＸ線スペクトルＳＬ２と、診療日２で生成された高ｋＶに対応するＸ線スペクトルＳＨ２を概略的に示す図である。
【図35】ステップＳＴ６０のフローの一例を示す図である。
【図37】ステップＳＴ６０４～ＳＴ６０６の説明図である。
【図38】診療日ｉで生成された低ｋＶに対応するＸ線スペクトルＳＬｉと、診療日ｉで生成された高ｋＶに対応するＸ線スペクトルＳＨｉを概略的に示す図である。
【図39】ステップＳＴ６０の説明図である。
【図40】ＣＴシステムの稼働再開後のフローの説明図である。
【図41】診療日ごとにエネルギー基準値を計算する場合のステップＳＴ６０のフローの説明図である。
【図42】診療日１～診療日ｉで取得された低ｋＶに対応するＸ線スペクトルＳＬ１～ＳＬｉを示す図である。
【発明を実施するための形態】
【0012】
以下、発明を実施するための形態について説明するが、本発明は、以下の形態に限定されることはない。
【0013】
【0014】
ガントリ１０２はボア１０７を有しており、そのボア１０７に被検体１１２が搬送され、被検体１１２がスキャンされる。
【0015】
ガントリ１０２には、Ｘ線発生装置１０４、フィルタ部１０３、前置コリメータ１０５、および検出器１０８などが取り付けられている。
【0016】
Ｘ線発生装置１０４は、Ｘ線管１０４Ａと、ジェネレータ１０４Ｂとを含んでいる。ジェネレータ１０４Ｂは、Ｘ線管１０４Ａに電力を供給する。Ｘ線管１０４Ａは、陰極－陽極管に所定の電圧が印加されることにより、Ｘ線を出力する。Ｘ線管１０４は、ＸＹ面内において、回転軸２０６を中心とした経路上を回転することができるように構成されている。ここで、Ｚ方向は体軸方向を表し、Ｙ方向は鉛直方向（テーブル１１６の高さ方向）を表し、Ｘ方向は、Ｚ方向およびＹ方向に対して垂直の方向を表している。本実施形態では、Ｘ線管１０４Ａは、Ｘ線管に印加される管電圧が第１の管電圧と第２の管電圧との間で交互に切り替えることができるｋＶスイッチング方式に対応したＸ線管である。尚、本実施形態では、ＣＴシステム１０は１つのＸ線管１０４Ａを備えているが、２つのＸ線管１０４Ａを備えてもよい。
【0017】
フィルタ部１０３は、例えば、平板フィルタおよび／又はボウタイフィルタを含んでいる。
前置コリメータ１０５は、不要な領域にＸ線が照射されないようにＸ線の照射範囲を絞り込むための部材である。
前置コリメータ１０５は、不要な領域にＸ線が照射されないようにＸ線の照射範囲を絞り込むための部材である。
【0018】
検出器１０８は複数の検出器素子２０２を含んでいる。複数の検出器素子２０２は、Ｘ線管１０４Ａから照射され、患者などの被検体１１２を通過するＸ線１０６を検出する。したがって、Ｘ線検出器１０８は、ビューごとに投影データを取得することができる。
【0019】
検出器１０８により検出された投影データは、ＤＡＳ２１４で収集される。ＤＡＳ２１４は、収集した投影データに対して、サンプリング、デジタル変換などを含む所定の処理を実行する。処理された投影データは、コンピュータ２１６に送信される。コンピュータ２１６は、ＤＡＳ２１４からのデータを記憶装置２１８に記憶する。記憶装置２１８は、プログラムや、プロセッサで実行される命令などを記録する1つ以上の記憶媒体を含むものである。記憶媒体は、例えば、1つ以上の非一時的なコンピュータ可読記憶媒体とすることができる。記憶装置２１８は、例えば、ハードディスクドライブ、フロッピーディスクドライブ、コンパクトディスク読み出し／書き込み（ＣＤ－Ｒ／Ｗ）ドライブ、デジタル多用途ディスク（ＤＶＤ）ドライブ、フラッシュドライブ、および／またはソリッドステート記録ドライブを含むことができる。
【0020】
コンピュータ２１６は１つ又は複数のプロセッサ２１７を含んでいる。コンピュータ２１６は、１つ又は複数のプロセッサを使用して、ＤＡＳ２１４、Ｘ線コントローラ２１０、および／又はガントリモータコントローラ２１２に、コマンドおよびパラメータを出力し、データ取得および／または処理などのシステム動作を制御する。また、コンピュータ２１６は、１つ又は複数のプロセッサを使用して、後述するフローの各ステップにおいて、信号処理、データ処理、画像処理など、様々な処理を実行する。尚、図１では、１つ又は複数のプロセッサ２１７はコンピュータ２１６に含まれているが、１つ又は複数のプロセッサ２１７は、コンピュータ２１６と他の構成要素（例えば、Ｘ線コントローラ２１０、ガントリモータコントローラ２１２、テーブルコントローラ１１８など）に分散させるように設けられていてもよい。
【0021】
コンピュータ２１６には、オペレータコンソール２２０が結合されている。オペレータは、オペレータコンソール２２０を操作することにより、ＣＴシステム１０の動作に関連する所定のオペレータ入力をコンピュータ２１６に入力することができる。コンピュータ２１６は、オペレータコンソール２２０を介して、コマンドおよび／またはスキャンパラメータを含むオペレータ入力を受信し、そのオペレータ入力に基づいてシステム動作を制御する。オペレータコンソール２２０は、オペレータがコマンドおよび／またはスキャンパラメータを指定するためのキーボード（図示せず）またはタッチスクリーンを含むことができる。
【0022】
Ｘ線コントローラ２１０は、コンピュータ２１６からの命令に基づいてＸ線発生装置１０４を制御する。また、ガントリモータコントローラ２１２は、コンピュータ２１６からの命令に基づいて、Ｘ線管１０４Ａおよび検出器１０８などの構成要素が回転するように、ガントリモータを制御する。
【0023】
図１は、１つのオペレータコンソール２２０のみを示しているが、２つ以上のオペレータコンソールをコンピュータ２１６に結合してもよい。
【0024】
また、ＣＴシステム１０は、例えば、有線ネットワークおよび／又は無線ネットワークを介して、遠隔に位置する複数のディスプレイ、プリンタ、ワークステーション、および／もしくは同様のデバイスが結合されるようにしてもよい。
【0025】
一実施形態では、例えば、ＣＴシステム１０は、画像保管通信システム（ＰＡＣＳ）２２４を含んでいてもよいし、ＰＡＣＳ２２４に結合されていてもよい。例示的な実施態様では、ＰＡＣＳ２２４は、放射線科情報システム、病院情報システム、および／または内部もしくは外部ネットワーク（図示せず）などの遠隔システムに結合されていてもよい。
【0026】
コンピュータ２１６は、テーブルモータコントローラ１１８に、テーブル１１６を制御するための命令を供給する。テーブルモータコントローラ１１８は、受け取った命令に基づいて、テーブル１１６が移動するように、テーブルモータを制御することができる。例えば、テーブルモータコントローラ１１８は、被検体１１２が撮影に適した位置に位置決めされるように、テーブル１１６を移動させることができる。
【0027】
前述のように、ＤＡＳ２１４は、検出器素子２０２によって取得された投影データをサンプリングしてデジタル変換する。その後、画像再構成器２３０が、サンプリングされデジタル変換されたデータを使用して画像を再構成する。画像再構成器２３０は１つ又は複数のプロセッサを含んでおり、このプロセッサが画像再構成の処理を実行することができる。図１では、画像再構成器２３０は、コンピュータ２１６とは別個の構成要素として示されているが、画像再構成器２３０は、コンピュータ２１６の一部を形成するものであってもよい。また、コンピュータ２１６が、画像再構成器２３０の１つまたは複数の機能を実施してもよい。さらに、画像再構成器２３０は、ＣＴシステム１０から離れた位置に設けられ、有線ネットワークまたは無線ネットワークを使用してＣＴシステム１０に動作可能に接続されるようにしてもよい。
【0028】
画像再構成器２３０は、再構成された画像を記憶装置２１８に記憶することができる。また、画像再構成器２３０は、再構成された画像をコンピュータ２１６に送信してもよい。コンピュータ２１６は、再構成された画像および／または患者情報を、コンピュータ２１６および／または画像再構成器２３０に通信可能に結合された表示装置２３２に送信することができる。
【0029】
本明細書で説明される様々な方法およびプロセスは、ＣＴシステム１０に含まれる又はＣＴシステム１０と通信可能である非一時的なコンピュータ可読記憶媒体に実行可能命令として記録することができる。この実行可能命令は、１つの記憶媒体に記録されていてもよいし、複数の記憶媒体に分散させて記録されるようにしてもよい。ＣＴシステム１０に備えられる１つ以上のプロセッサは、記憶媒体に記録された命令に従って、本明細書で説明される様々な方法、ステップ、およびプロセスを実行する。
【0030】
ＣＴシステム１０は上記のように構成されている。本実施形態のＣＴシステムは、ＤＥＣＴ（Dual Energy CT）のｋＶスイッチング法」に対応しており、物質の弁別を行うことができる。しかし、ＤＥＣＴのｋＶスイッチング法による物質の弁別精度は、Ｘ線発生装置の経年劣化によって低下する。したがって、物質の弁別精度を低下させないようにするために、ＣＴシステムのユーザは、Ｘ線発生装置の経年劣化をできるだけ早急に知ることが重要となる。しかし、現状はＸ線発生装置の劣化を検出することは行われていない。したがって、ＣＴシステムのユーザは、ＣＴ画像の画質が低下していることを視覚的に認識することによって、Ｘ線発生装置の劣化に気が付く。このため、ＣＴ画像の画質が低下する前に、Ｘ線発生装置に経年劣化が発生している恐れがあることをユーザに知らせる技術が要求されている。
【0031】
そこで、発明者等は、鋭意研究し、ＣＴシステムのユーザがＣＴ画像の画質低下を視覚的に認識する前に、ユーザに、Ｘ線発生装置に経年劣化が発生している恐れがあることを報知する方法を考え出した。以下に、本方法について説明する。
【0032】
本実施形態について具体的に説明する前に、本実施形態の特徴を要約すると、以下の通りである。本実施形態では、定期的にキャリブレーションスキャンを実行する。プロセッサは、キャリブレーションスキャンが実行されることにより検出器１０８で検出されたＸ線のデータに基づいて、Ｘ線発生装置１０４から照射されたＸ線の特性値を計算する。そして、プロセッサは、特性値の経時的な変化に基づいて、アラートを出力するかどうかを判断する。したがって、特性値の経時的な変化が大きい場合にアラートが出力されるようにすることで、ＣＴ画像の画質が低下する前に、Ｘ線発生装置に劣化の徴候があることをユーザに知らせることができる。
【0033】
また、キャリブレーションスキャンは、安定したＣＴ画像を得るために、定期的に（例えば、毎朝）実行することが推奨されているスキャンである。したがって、これまで通りにキャリブレーションスキャンを定期的に実行していけば、必要に応じてアラートが出力され、ユーザに知らせることができる。このため、キャリブレーションスキャンとは別に、追加のスキャンを実行しなくても、ＣＴ画像の画質が低下する前に、Ｘ線発生装置に劣化の徴候があることをユーザに知らせることができる。
【0034】
以下に、本実施形態のＣＴシステムにおいて、Ｘ線の特性値の経時的な変化に基づいてアラートを出力するかどうかを判断する方法について説明する。尚、以下では、Ｘ線の特性値として、２つの特性値、即ち、第１の基準物質の吸収係数および第２の基準物質の吸収係数を考えることにする。
【0035】
図２は、各診療日に実行されるフローを示す図である。以下では、説明の便宜上、診療日にキャリブレーションスキャンが一回実行する例を取り上げて、本方法を説明するが、キャリブレーションスキャンが実行される頻度は、一日一回に限定されることはない。例えば、キャリブレーションスキャンは、一日に複数回（例えば、朝と昼）実行されてもよいし、キャリブレーションスキャンが実行されない日を定期的に設定してもよい。
【0036】
先ず、診療日１について説明する。診療日１では、フロー１００が実行される。以下、フロー１００について説明する。
【0037】
ステップＳＴ１１では、Ｘ線管１０４Ａに印加される管電圧が第１の管電圧（低ｋＶ）と第２の管電圧（高ｋＶ）との間で交互に切り替えられるｋＶスイッチングを使用したキャリブレーションスキャンが実行される。
【0038】
図３は、キャリブレーションの説明図である。
ｋＶスイッチング法のキャリブレーションでは、例えば、３つのパラメータ（回転速度，コーン角，管電流）の値の組合せとしてｚ通りの組合せＰ１～Ｐｚが用意されている。そして、ｚ通りの組合せＰ１～ＰｚをプリセットしてＣＴシステムに保存しておき、各組合せに対してキャリブレーションスキャンＳ１～Ｓｚが実行される。キャリブレーションスキャンを実行することにより、検出器１０８でＸ線が検出される。プロセッサは、検出器１０８で検出されたＸ線のデータに基づいて、キャリブレーションデータを生成する。図２では、ｚ個のプリセットＰ１～Ｐｚに対して、ｚ個のキャリブレーションデータＤ１～Ｄｚが生成された例が示されている。各キャリブレーションデータには、例えば、吸収係数、ビームハードニング用の補正係数など、様々な係数などが含まれているが、ここでは、キャリブレーションデータとして、本実施形態の説明に関係する吸収係数を表すキャリブレーションデータのみを考えることにする。
ｋＶスイッチング法のキャリブレーションでは、例えば、３つのパラメータ（回転速度，コーン角，管電流）の値の組合せとしてｚ通りの組合せＰ１～Ｐｚが用意されている。そして、ｚ通りの組合せＰ１～ＰｚをプリセットしてＣＴシステムに保存しておき、各組合せに対してキャリブレーションスキャンＳ１～Ｓｚが実行される。キャリブレーションスキャンを実行することにより、検出器１０８でＸ線が検出される。プロセッサは、検出器１０８で検出されたＸ線のデータに基づいて、キャリブレーションデータを生成する。図２では、ｚ個のプリセットＰ１～Ｐｚに対して、ｚ個のキャリブレーションデータＤ１～Ｄｚが生成された例が示されている。各キャリブレーションデータには、例えば、吸収係数、ビームハードニング用の補正係数など、様々な係数などが含まれているが、ここでは、キャリブレーションデータとして、本実施形態の説明に関係する吸収係数を表すキャリブレーションデータのみを考えることにする。
【0039】
尚、吸収係数を解析する場合、全てのキャリブレーションデータＤ１～Ｄｚの吸収係数を解析する必要は無く、キャリブレーションスキャンＳ１～Ｓｚのうちの特定のキャリブレーションスキャンで得られたキャリブレーションデータの吸収係数の経時的な変化を解析すればよい。そこで、以下では、キャリブレーションスキャンＳ１により得られたキャリブレーションデータＤ１にのみ着目し、吸収係数を説明することにする。
【0040】
吸収係数には、第１の基準物質の吸収係数と、第２の基準物質の吸収係数が含まれている。したがって、キャリブレーションスキャンを実行することにより、第１の基準物質の吸収係数と、第２の基準物質の吸収係数が計算される。２つの基準物質は、例えば、水とヨードの組合せや、水とカルシウムの組合せが考えられる。図４は、診療日１において取得されたキャリブレーションデータの説明図である。
【0041】
図４には、診療日１において取得されたキャリブレーションデータの第１の基準物質の吸収係数を表す吸収係数データｇ１１～ｇ１ｍと、第２の基準物質の吸収係数を表す吸収係数データｈ１１～ｈ１ｍの概略が示されている。
【0042】
吸収係数データｇ１１は、検出器の第１列の各チャネルにおける第１の基準物質の吸収係数を表しており、吸収係数データｇ１２～ｇ１ｍは、それぞれ、検出器の第２列～第ｍ列の各チャネルにおける第１の基準物質の吸収係数を表している。また、吸収係数データｈ１１は、検出器の第１列の各チャネルにおける第２の基準物質の吸収係数を表しており、吸収係数データｈ１２～ｈ１ｍは、それぞれ、検出器の第２列～第ｍ列の各チャネルにおける第２の基準物質の吸収係数を表している。例えば、ｍ＝６４、即ち、検出器が６４列の多列構造を有している場合、第１の基準物質については、６４個の吸収係数データｇ１１～ｇ１，６４が生成され、第２の基準物質については６４個の吸収係数データｈ１１～ｈ１，６４が生成される。尚、以下では、吸収係数データの説明を簡略化するため、ｍ＝１、即ち、検出器が１列構造を有している場合を考えることにする。したがって、第１の基準物質については、吸収係数データｇ１１のみを考え、第２の基準物質については吸収係数データｈ１１のみを考えることにする。
【0043】
図５は、吸収係数データｇ１１の拡大図であり、図６は、吸収係数データｈ１１の拡大図である。
尚、以下に示す吸収係数データの波形は、実施形態を説明するために示されたものであり、実際の波形とは異なっていることに留意されたい。
尚、以下に示す吸収係数データの波形は、実施形態を説明するために示されたものであり、実際の波形とは異なっていることに留意されたい。
【0044】
図５では、検出器のチャネル１～ｎにおける第１の基準物質の吸収係数のうち、代表して、チャネル１、チャネル２、チャネルｊ、および、チャネルｎの吸収係数のみを示している。チャネル１では、吸収係数は「ａ１１」で示されており、チャネル２では、吸収係数は「ａ１２」で示されており、チャネルｊでは、吸収係数は「ａ１ｊ」で示されており、チャネルｎでは、吸収係数は「ａ１ｎ」で示されている。
【0045】
図６では、検出器のチャネル１～ｎにおける第２の基準物質の吸収係数のうち、代表して、チャネル１、チャネル２、チャネルｊ、および、チャネルｎの吸収係数のみを示している。チャネル１では、吸収係数は「ｂ１１」で示されており、チャネル２では、吸収係数は「ｂ１２」で示されており、チャネルｊでは、吸収係数は「ｂ１ｊ」で示されており、チャネルｎでは、吸収係数は「ｂ１ｎ」で示されている。
【0046】
本実施形態では、診療日１に計算された第１の基準物質の吸収係数ａ１１～ａ１ｎは、アラートを出力するかどうかを判断するための基準となる第１の基準吸収係数として保存される。また、診療日１に計算された第２の基準物質の吸収係数ｂ１１～ｂ１ｎは、アラートを出力するかどうかを判断するための基準となる第２の基準吸収係数として保存される。
【0047】
ステップＳＴ１１が実行された後は、ステップＳＴ１２に進み、診療日１の検査スケジュールに従って被検体の検査が実行される。
【0048】
次に、診療日２について説明する。診療日２では、フロー２００が実行される。
先ず、ステップＳＴ１１において、キャリブレーションスキャンが実行される。プロセッサは、キャリブレーションスキャンで得られたデータに基づいて、キャリブレーションデータ（吸収係数データなど）を生成する。
先ず、ステップＳＴ１１において、キャリブレーションスキャンが実行される。プロセッサは、キャリブレーションスキャンで得られたデータに基づいて、キャリブレーションデータ（吸収係数データなど）を生成する。
【0049】
【0050】
尚、図７には、診療日２で得られた吸収係数データｇ２１の他に、診療日１で得られた吸収係数データｇ１１も示されている。また、図８には、診療日２で得られた吸収係数データｈ２１の他に、診療日１で得られた吸収係数データｈ１１も示されている。
ステップＳＴ１１を実行した後、ステップＳＴ２０に進む。
ステップＳＴ１１を実行した後、ステップＳＴ２０に進む。
【0051】
ステップＳＴ２０では、吸収係数の経時的な変化に基づいて、アラートを出力するかどうか判断する。
図９は、ステップＳＴ２０のフローの一例を示す図、図１０は、図９のフローの説明図である。図１０には、図７に示す第１の基準物質の吸収係数データｇ１１およびｇ２１が示されている。
図９は、ステップＳＴ２０のフローの一例を示す図、図１０は、図９のフローの説明図である。図１０には、図７に示す第１の基準物質の吸収係数データｇ１１およびｇ２１が示されている。
【0052】
ステップＳＴ２０１において、プロセッサは、基準吸収係数（診療日１の吸収係数）と診療日２の吸収係数との差を計算する。具体的には、プロセッサは、チャネル１における診療日１の基準吸収係数ａ１１と診療日２の吸収係数ａ２１との差Δａ２１を計算する。吸収係数の差Δａ２１を計算したら、ステップＳＴ２０２に進む。
【0053】
ステップＳＴ２０２では、プロセッサは、Δａ２１と閾値ＴＨ１とを比較し、その比較結果に基づいて、アラートを出力するかどうか判断する。例えば、Δａ２１が閾値ＴＨ１を超えた場合（Δａ２１＞ＴＨ１）、プロセッサは、吸収係数の経時的な変化が大きいので、アラートを出力すると判断する。この場合、ステップＳＴ２１に進み、プロセッサは、表示部にアラートを出力させる。
【0054】
一方、吸収係数の差Δａ２１が閾値ＴＨ１を超えていない場合（Δａ２１≦ＴＨ１）、吸収係数の経時的な変化は小さいので、プロセッサは、アラートは出力しないと判断する。この場合、ステップＳＴ２０３に進む。
【0055】
ステップＳＴ２０３では、プロセッサは、全チャネルに対してステップＳＴ２０１およびＳＴ２０２が実行されたかどうかを判断する。ここでは、チャネル１に対してステップＳＴ２０１およびＳＴ２０２が実行されたが、他のチャネル２～チャネルｎに対しては、まだ実行されていない。したがって、ステップＳＴ２０４に進み、プロセッサは、チャネルをチャネル１からチャネル２にインクリメントする。チャネルがインクリメントされると、ステップＳＴ２０１に戻る。
【0056】
図１１は、チャネル２に対してステップＳＴ２０１およびＳＴ２０２を実行するときの説明図である。
ステップＳＴ２０１において、プロセッサは、チャネル２における診療日１の吸収係数ａ１２と診療日２の吸収係数ａ２２との差Δａ２２を計算する。吸収係数の差Δａ２２を計算したら、ステップＳＴ２０２に進む。
ステップＳＴ２０１において、プロセッサは、チャネル２における診療日１の吸収係数ａ１２と診療日２の吸収係数ａ２２との差Δａ２２を計算する。吸収係数の差Δａ２２を計算したら、ステップＳＴ２０２に進む。
【0057】
ステップＳＴ２０２では、プロセッサは、Δａ２２と閾値ＴＨ１とを比較し、その比較結果に基づいて、アラートを出力するかどうか判断する。例えば、Δａ２２が閾値ＴＨ１を超えた場合（Δａ２２＞ＴＨ１）、プロセッサは、吸収係数の経時的な変化が大きいので、アラートを出力すると判断する。この場合、ステップＳＴ２１に進み、プロセッサは、表示部にアラートを出力させる。
【0058】
一方、吸収係数の差Δａ２２が閾値ＴＨ１を超えていない場合（Δａ２２≦ＴＨ１）、吸収係数の経時的な変化は小さいので、プロセッサは、アラートは出力しないと判断する。この場合、ステップＳＴ２０３に進む。
【0059】
ステップＳＴ２０３では、プロセッサは、全チャネルに対してステップＳＴ２０１およびＳＴ２０２が実行されたかどうか判断する。ここでは、チャネル１および２に対してステップＳＴ２０１およびＳＴ２０２が実行されたが、他のチャネル３～チャネルｎに対しては、まだ実行されていない。したがって、ステップＳＴ２０４に進み、プロセッサは、チャネルをチャネル２からチャネル３にインクリメントする。そして、ステップＳＴ２０１に戻る。
【0060】
以下同様に、プロセッサは、ステップＳＴ２０４でチャネルがインクリメントされるたびに、ステップＳＴ２０１に戻り、ステップＳＴ２０１～ＳＴ２０４のループを繰り返し実行する（図１２参照）。
【0061】
図１２は、チャネルｊに対してステップＳＴ２０１およびＳＴ２０２を実行するときの説明図である。
ステップＳＴ２０１において、プロセッサは、チャネルｊにおける診療日１の吸収係数ａ１ｊと診療日２の吸収係数ａ２ｊとの差Δａ２ｊを計算する。吸収係数の差Δａ２ｊを計算したら、ステップＳＴ２０２に進む。
ステップＳＴ２０１において、プロセッサは、チャネルｊにおける診療日１の吸収係数ａ１ｊと診療日２の吸収係数ａ２ｊとの差Δａ２ｊを計算する。吸収係数の差Δａ２ｊを計算したら、ステップＳＴ２０２に進む。
【0062】
ステップＳＴ２０２では、プロセッサは、Δａ２ｊと閾値ＴＨ１とを比較し、その比較結果に基づいて、アラートを出力するかどうか判断する。例えば、Δａ２ｊが閾値ＴＨ１を超えた場合（Δａ２ｊ＞ＴＨ１）、プロセッサは、吸収係数の経時的な変化が大きいので、アラートを出力すると判断する。この場合、ステップＳＴ２１に進み、プロセッサは、表示部にアラートを出力させる。
【0063】
一方、吸収係数の差Δａ２ｊが閾値ＴＨ１を超えていない場合（Δａ２ｊ≦ＴＨ１）、吸収係数の経時的な変化が小さいので、プロセッサは、アラートは出力しないと判断する。この場合、ステップＳＴ２０３に進む。
【0064】
以下同様に、プロセッサは、ステップＳＴ２０４でチャネルがインクリメントされるたびに、ステップＳＴ２０１に戻り、ステップＳＴ２０１～ＳＴ２０４のループを繰り返し実行する（図１３参照）。
【0065】
図１３は、チャネルｎに対してステップＳＴ２０１およびＳＴ２０２を実行するときの説明図である。
ステップＳＴ２０１において、プロセッサは、チャネルｎにおける診療日１の吸収係数ａ１ｎと診療日２の吸収係数ａ２ｎとの差Δａ２ｎを計算する。吸収係数の差Δａ２ｎを計算したら、ステップＳＴ２０２に進む。
ステップＳＴ２０１において、プロセッサは、チャネルｎにおける診療日１の吸収係数ａ１ｎと診療日２の吸収係数ａ２ｎとの差Δａ２ｎを計算する。吸収係数の差Δａ２ｎを計算したら、ステップＳＴ２０２に進む。
【0066】
ステップＳＴ２０２では、プロセッサは、Δａ２ｎと閾値ＴＨ１とを比較し、その比較結果に基づいて、アラートを出力するかどうか判断する。例えば、Δａ２ｎが閾値ＴＨ１を超えた場合（Δａ２ｎ＞ＴＨ１）、プロセッサは、吸収係数の経時的な変化が大きいので、アラートを出力すると判断する。この場合、ステップＳＴ２１に進み、プロセッサは、表示部にアラートを出力させる。
【0067】
一方、吸収係数の差Δａ２ｎが閾値ＴＨ１を超えていない場合（Δａ２ｎ≦ＴＨ１）、吸収係数の経時的な変化が小さいので、プロセッサは、アラートは出力しないと判断する。この場合、ステップＳＴ２０３に進む。
【0068】
ステップＳＴ２０３では、プロセッサは、全チャネルに対してステップＳＴ２０１～ＳＴ２０４が実行されたかどうか判断する。ここでは、全チャネル１～ｎに対してステップＳＴ２０１～ＳＴ２０４が実行されている。したがって、ステップＳＴ２０５に進む。
【0069】
ステップＳＴ２０５では、プロセッサは、第１および第２の基準物質の両方の吸収係数の解析が終了したかどうか判断する。ここでは、第１の基準物質の吸収係数は解析したが、第２の基準物質の吸収係数は解析されていない。したがって、ステップＳＴ２０１に戻る。
【0070】
ステップＳＴ２０１に戻ると、プロセッサは、第１の基準物質に対して説明したように、第２の基準物質の吸収係数データｈ１１およびｈ２１（図８参照）に対して、ステップＳＴ２０１～ＳＴ２０４の処理を実行する。
【0071】
【0072】
ステップＳＴ２０１において、プロセッサは、チャネル１における診療日１の基準吸収係数ｂ１１と診療日２の吸収係数ｂ２１との差Δｂ２１を計算する。吸収係数の差Δｂ２１を計算したら、ステップＳＴ２０２に進む。
【0073】
ステップＳＴ２０２では、プロセッサは、Δｂ２１と閾値ＴＨ２とを比較し、その比較結果に基づいて、アラートを出力するかどうか判断する。例えば、Δｂ２１が閾値ＴＨ２を超えた場合（Δｂ２１＞ＴＨ２）、プロセッサは、吸収係数の経時的な変化が大きいので、アラートを出力すると判断する。この場合、ステップＳＴ２１に進み、プロセッサは、表示部にアラートを出力させる。
【0074】
一方、吸収係数の差Δｂ２１が閾値ＴＨ２を超えていない場合（Δｂ２１≦ＴＨ２）、吸収係数の経時的な変化が小さいので、プロセッサは、アラートは出力しないと判断する。この場合、ステップＳＴ２０３に進む。
【0075】
ステップＳＴ２０３では、プロセッサは、全チャネルに対してステップＳＴ２０１およびＳＴ２０２が実行されたかどうかを判断する。ここでは、チャネル１に対してステップＳＴ２０１およびＳＴ２０２が実行されたが、他のチャネル２～チャネルｎに対しては、まだ実行されていない。したがって、ステップＳＴ２０４に進み、プロセッサは、チャネルをチャネル１からチャネル２にインクリメントする。チャネルがインクリメントされると、ステップＳＴ２０１に戻る。
【0076】
以下同様に、ステップＳＴ２０１～ＳＴ２０４のループが実行される。そして、吸収係数データｈ１１およびｈ２１のチャネル１～ｎに対してスステップＳＴ２０１～ＳＴ２０４の処理を実行している間に、あるチャネルにおいて、吸収係数の差が閾値ＴＨ２を超えた場合、アラートを出力すると判断する。
【0077】
一方、プロセッサは、第２の基準物質の吸収係数データｈ１１およびｈ１２の全チャネルに対してステップＳＴ２０１～ＳＴ２０４を実行したが、どのチャネルにおいても吸収係数の差が閾値ＴＨ２以下である場合、ステップＳＴ２０５に進む。
【0078】
ステップＳＴ２０５では、プロセッサは、第１および第２の基準物質の両方の吸収係数の解析が終了したかどうか判断する。ここでは、第１および第２の基準物質の両方の吸収係数の解析が終了している。したがって、プロセッサは、アラートの出力は不要であると判断して、ステップＳＴ２０を終了する。この場合、ステップＳＴ２２に進み、診療日２の検査スケジュールに従って、被検体の検査が実行される。
【0079】
本実施形態では、診療日２では、アラートは出力されなかったとする。したがって、ステップＳＴ２０の終了後、被検体の検査が実行される。
図２に戻って説明を続ける。
図２に戻って説明を続ける。
【0080】
次に、診療日３について説明する。診療日３では、診療日２と同様に、フロー２００が実行される。
先ず、ステップＳＴ１１において、キャリブレーションスキャンが実行される。プロセッサは、キャリブレーションスキャンで得られたデータに基づいて、キャリブレーションデータ（吸収係数データなど）を生成する。
先ず、ステップＳＴ１１において、キャリブレーションスキャンが実行される。プロセッサは、キャリブレーションスキャンで得られたデータに基づいて、キャリブレーションデータ（吸収係数データなど）を生成する。
【0081】
【0082】
ステップＳＴ１１を実行した後、ステップＳＴ２０に進む。
ステップＳＴ２０では、吸収係数の経時的な変化に基づいて、アラートを出力するかどうか判断する。
ステップＳＴ２０では、吸収係数の経時的な変化に基づいて、アラートを出力するかどうか判断する。
【0083】
図１７は、ステップＳＴ２０の説明図である。
ステップＳＴ２０１において、プロセッサは、チャネル１における診療日１の基準吸収係数ａ１１と診療日２の吸収係数ａ３１との差Δａ３１を計算する。吸収係数の差Δａ３１を計算したら、ステップＳＴ２０２に進む。
ステップＳＴ２０１において、プロセッサは、チャネル１における診療日１の基準吸収係数ａ１１と診療日２の吸収係数ａ３１との差Δａ３１を計算する。吸収係数の差Δａ３１を計算したら、ステップＳＴ２０２に進む。
【0084】
ステップＳＴ２０２では、プロセッサは、Δａ３１と閾値ＴＨ１とを比較し、その比較結果に基づいて、アラートを出力するかどうか判断する。例えば、Δａ３１が閾値ＴＨ１を超えた場合（Δａ３１＞ＴＨ１）、プロセッサは、吸収係数の経時的な変化が大きいので、アラートを出力すると判断する。この場合、ステップＳＴ２１に進み、プロセッサは、表示部にアラートを出力させる。
【0085】
一方、吸収係数の差Δａ３１が閾値ＴＨ１を超えていない場合（Δａ３１≦ＴＨ１）、吸収係数の経時的な変化が小さいので、プロセッサは、アラートは出力しないと判断する。この場合、ステップＳＴ２０３に進む。
【0086】
ステップＳＴ２０３では、プロセッサは、全チャネルに対してステップＳＴ２０１およびＳＴ２０２が実行されたかどうかを判断する。ここでは、チャネル１に対してステップＳＴ２０１およびＳＴ２０２が実行されたが、他のチャネル２～チャネルｎに対しては、まだ実行されていない。したがって、ステップＳＴ２０４に進み、プロセッサは、チャネルをチャネル１からチャネル２にインクリメントする。チャネルがインクリメントされると、ステップＳＴ２０１に戻る。
【0087】
以下同様に、ステップＳＴ２０１～ＳＴ２０４のループが実行される。そして、吸収係数データｇ１１およびｇ３１のチャネル１～ｎに対してステップＳＴ２０１～ＳＴ２０４の処理を実行している間に、あるチャネルにおいて、吸収係数の差が閾値ＴＨ１を超えた場合、アラートを出力すると判断する。
【0088】
一方、プロセッサは、第１の基準物質の吸収係数データｇ１１およびｇ３１の全チャネルに対してステップＳＴ２０１～ＳＴ２０４を実行したが、どのチャネルにおいても吸収係数の差が閾値ＴＨ１以下である場合、ステップＳＴ２０５に進む。ステップＳＴ２０５では、プロセッサは、第１および第２の基準物質の両方の吸収係数の解析が終了したかどうか判断する。ここでは、第１の基準物質の吸収係数の解析は終了したが、第２の基準物質の吸収係数は解析されていない。したがって、ステップＳＴ２０１に戻る。
【0089】
ステップＳＴ２０１に戻ると、プロセッサは、第２の基準物質の吸収係数に対して、ステップＳＴ２０１～ＳＴ２０４の処理を実行する。
【0090】
図１８は、第２の基準物質の吸収係数データに対してステップＳＴ２０１～ＳＴ２０４を実行するときの説明図である。
【0091】
ステップＳＴ２０１において、プロセッサは、チャネル１における診療日１の吸収係数ｂ１１と診療日２の吸収係数ｂ２１との差Δｂ２１を計算する。吸収係数の差Δｂ３１を計算したら、ステップＳＴ２０２に進む。
【0092】
ステップＳＴ２０２では、プロセッサは、Δｂ３１と閾値ＴＨ２とを比較し、その比較結果に基づいて、アラートを出力するかどうか判断する。例えば、Δｂ３１が閾値ＴＨ２を超えた場合（Δｂ３１＞ＴＨ２）、プロセッサは、吸収係数の経時的な変化が大きいので、アラートを出力すると判断する。この場合、ステップＳＴ２１に進み、プロセッサは、表示部にアラートを出力させる。
【0093】
一方、吸収係数の差Δｂ３１が閾値ＴＨ２を超えていない場合（Δｂ３１≦ＴＨ２）、吸収係数の経時的な変化が小さいので、プロセッサは、アラートは出力しないと判断する。この場合、ステップＳＴ２０３に進む。
【0094】
ステップＳＴ２０３では、プロセッサは、全チャネルに対してステップＳＴ２０１およびＳＴ２０２が実行されたかどうかを判断する。ここでは、チャネル１に対してステップＳＴ２０１およびＳＴ２０２が実行されたが、他のチャネル２～チャネルｎに対しては、まだ実行されていない。したがって、ステップＳＴ２０４に進み、プロセッサは、チャネルをチャネル１からチャネル２にインクリメントする。チャネルがインクリメントされると、ステップＳＴ２０１に戻る。
【0095】
以下同様に、ステップＳＴ２０１～ＳＴ２０４のループが実行される。そして、吸収係数データｈ１１およびｈ３１のチャネル１～ｎに対してステップＳＴ２０１～ＳＴ２０４の処理を実行している間に、あるチャネルにおいて、吸収係数の差が閾値ＴＨ２を超えた場合、アラートを出力すると判断する。
【0096】
一方、プロセッサは、第２の基準物質の吸収係数データｈ１１およびｈ１３の全チャネルに対してステップＳＴ２０１～ＳＴ２０４を実行したが、どのチャネルにおいても吸収係数の差が閾値ＴＨ２以下である場合、ステップＳＴ２０５に進む。
【0097】
ステップＳＴ２０５では、プロセッサは、第１および第２の基準物質の両方の吸収係数の解析が終了したかどうか判断する。ここでは、第１の基準物質および第２の基準物質の両方の吸収係数の解析が終了している。したがって、プロセッサは、アラートの出力は不要であると判断して、ステップＳＴ２０を終了する。この場合、ステップＳＴ２２に進み、診療日３の検査スケジュールに従って、被検体の検査が実行される。
【0098】
本実施形態では、診療日３では、アラートは出力されなかったとする。したがって、ステップＳＴ２０の終了後、被検体の検査が実行される。
【0099】
次に、診療日ｉについて説明する。診療日ｉでは、診療日２と同様に、フロー２００が実行される。
ステップＳＴ１１では、キャリブレーションスキャンが実行される。プロセッサは、キャリブレーションスキャンで得られたデータに基づいて、キャリブレーションデータ（吸収係数データなど）を生成する。図１９には、診療日ｉに得られた第１の基準物質の吸収係数データｇｉ１が示されており、図２０には、診療日ｉに得られた第２の基準物質の吸収係数データｈｉ１が示されている。
ステップＳＴ１１では、キャリブレーションスキャンが実行される。プロセッサは、キャリブレーションスキャンで得られたデータに基づいて、キャリブレーションデータ（吸収係数データなど）を生成する。図１９には、診療日ｉに得られた第１の基準物質の吸収係数データｇｉ１が示されており、図２０には、診療日ｉに得られた第２の基準物質の吸収係数データｈｉ１が示されている。
【0100】
ステップＳＴ１１を実行した後、ステップＳＴ２０に進む。
【0101】
ステップＳＴ２０では、吸収係数の経時的な変化に基づいて、アラートを出力するかどうか判断する。
【0102】
【0103】
ステップＳＴ２０１において、プロセッサは、チャネル１における診療日１の吸収係数ａ１１と診療日ｉの吸収係数ａｉ１との差Δａｉ１を計算する。吸収係数の差Δａｉ１を計算したら、ステップＳＴ２０２に進む。
【0104】
ステップＳＴ２０２では、プロセッサは、Δａｉ１と閾値ＴＨ１とを比較し、その比較結果に基づいて、アラートを出力するかどうか判断する。例えば、Δａｉ１が閾値ＴＨ１を超えた場合（Δａｉ１＞ＴＨ１）、プロセッサは、吸収係数の経時的な変化が大きいので、アラートを出力すると判断する。この場合、ステップＳＴ２１に進み、プロセッサは、表示部にアラートを出力させる。
【0105】
一方、吸収係数の差Δａｉ１が閾値ＴＨ１を超えていない場合（Δａｉ１≦ＴＨ１）、吸収係数の経時的な変化が小さいので、プロセッサは、アラートは出力しないと判断する。この場合、ステップＳＴ２０３に進む。
【0106】
ここでは、チャネル１において計算された吸収係数の差Δａｉ１が、Δａｉ１＞ＴＨ１を満たすとする。したがって、アラートが出力される。
【0107】
アラートが出力されると、ユーザは、例えば、サービスに点検１１（図２参照）を依頼することができる。また、アラートが出力された場合、ＣＴシステムは、吸収係数の時間変化を表すデータをバックオフィスに送り、バックオフィスで、吸収係数の時間変化を表すデータを解析するようにしてもよい。
【0108】
サービスによる点検１１の結果、Ｘ線発生装置に障害が発生している又はＸ線発生装置が故障していると認められた場合、Ｘ線発生装置が修理又は交換される。ＣＴシステムは、Ｘ線発生装置の修理又は交換の後、動作チェックが行われ、正常に動作することが確認できたら、ＣＴシステムの稼働を再開する。図２２は、ＣＴシステムの稼働再開後のフローの説明図である。
ここでは、診療日ｉに点検等が完了し、次の診療日ｐにＣＴシステムの稼働が再開できたとする。
ここでは、診療日ｉに点検等が完了し、次の診療日ｐにＣＴシステムの稼働が再開できたとする。
【0109】
以下、診療日ｐについて説明する。
診療日ｐでは、診療日１と同じフロー１００が実行される。したがって、ステップＳＴ１１において、キャリブレーションスキャンが実行される。プロセッサは、キャリブレーションスキャンで得られたデータに基づいて、キャリブレーションデータ（吸収係数データなど）を生成する。
診療日ｐでは、診療日１と同じフロー１００が実行される。したがって、ステップＳＴ１１において、キャリブレーションスキャンが実行される。プロセッサは、キャリブレーションスキャンで得られたデータに基づいて、キャリブレーションデータ（吸収係数データなど）を生成する。
【0110】
【0111】
吸収係数データｇｐ１は、検出器の第１列の各チャネルにおける第１の基準物質の吸収係数を表している。また、吸収係数データｈｐ１は、検出器の第１列の各チャネルにおける第２の基準物質の吸収係数を表している。図２３および図２４では、検出器のチャネル１～ｎにおける第１の基準物質の吸収係数のうち、代表して、チャネル１、チャネル２、チャネルｊ、および、チャネルｎの吸収係数のみを示している。
【0112】
診療日１～診療日ｉでは、診療日１に取得された第１の基準物質の吸収係数ａ１１～ａ１ｎが、診療日２～診療日ｉで取得された第１の基準物質の吸収係数（例えば、吸収係数ａ２１、ａ３１、ａｉ１）の経時的な変化量を求めるための基準値として使用されていた。しかし、Ｘ線発生装置は診療日ｉに点検され、次の診療日ｐにＣＴシステムの稼働が再開されている。したがって、Ｘ線発生装置の点検前に取得された診療日１の吸収係数ａ１１～ａ１ｎは、第１の基準物質の吸収係数の経時的な変化量を求めるための基準値として使用することはできない。そこで、プロセッサは、診療日ｐに取得された吸収係数ａｐ１～ａｐｎを、第１の基準物質の吸収係数の経時的な変化量を求めるための新たな第１の基準値として保存する。
【0113】
同様に、プロセッサは、診療日ｐに取得された吸収係数ｂｐ１～ｂｐｎを、第２の基準物質の吸収係数の経時的な変化量を求めるための新たな第２の基準値として保存する。
【0114】
ステップＳＴ１１が実行された後は、ステップＳＴ１２に進み、被検体の検査が実行される。
【0115】
次に、診療日ｑについて説明する。
診療日ｑでは、診療日２と同様に、フロー２００が実行される。したがって、ステップＳＴ１１において、キャリブレーションスキャンが実行される。プロセッサは、キャリブレーションスキャンで得られたデータに基づいて、キャリブレーションデータ（吸収係数データなど）を生成する。図２５には、診療日ｑに得られた第１の基準物質の吸収係数データｇｑ１が示されており、図２６には、診療日ｑに得られた第２の基準物質の吸収係数データｈｑ１が示されている。
診療日ｑでは、診療日２と同様に、フロー２００が実行される。したがって、ステップＳＴ１１において、キャリブレーションスキャンが実行される。プロセッサは、キャリブレーションスキャンで得られたデータに基づいて、キャリブレーションデータ（吸収係数データなど）を生成する。図２５には、診療日ｑに得られた第１の基準物質の吸収係数データｇｑ１が示されており、図２６には、診療日ｑに得られた第２の基準物質の吸収係数データｈｑ１が示されている。
【0116】
ステップＳＴ１１を実行した後、ステップＳＴ２０に進む。
ステップＳＴ２０では、吸収係数の経時的な変化に基づいて、アラートを出力するかどうか判断する。
ステップＳＴ２０では、吸収係数の経時的な変化に基づいて、アラートを出力するかどうか判断する。
【0117】
図２７は、ステップＳＴ２０の説明図である。
ステップＳＴ２０１において、プロセッサは、チャネル１における診療日ｐの吸収係数ａｐ１と診療日ｑの吸収係数ａｑ１との差Δａｑ１を計算する。吸収係数の差Δａｑ１を計算したら、ステップＳＴ２０２に進む。
ステップＳＴ２０１において、プロセッサは、チャネル１における診療日ｐの吸収係数ａｐ１と診療日ｑの吸収係数ａｑ１との差Δａｑ１を計算する。吸収係数の差Δａｑ１を計算したら、ステップＳＴ２０２に進む。
【0118】
ステップＳＴ２０２では、プロセッサは、Δａｑ１と閾値ＴＨ１とを比較し、その比較結果に基づいて、アラートを出力するかどうか判断する。例えば、Δａｑ１が閾値ＴＨ１を超えた場合（Δａｑ１＞ＴＨ１）、プロセッサは、吸収係数の経時的な変化が大きいので、アラートを出力すると判断する。この場合、ステップＳＴ２１に進み、プロセッサは、表示部にアラートを出力させる。
【0119】
一方、吸収係数の差Δａｑ１が閾値ＴＨ１を超えていない場合（Δａｑ１≦ＴＨ１）、吸収係数の経時的な変化が小さいので、プロセッサは、アラートは出力しないと判断する。この場合、ステップＳＴ２０３に進む。
【0120】
ステップＳＴ２０３では、プロセッサは、全チャネルに対してステップＳＴ２０１およびＳＴ２０２が実行されたかどうかを判断する。ここでは、チャネル１に対してステップＳＴ２０１およびＳＴ２０２が実行されたが、他のチャネル２～チャネルｎに対しては、まだ実行されていない。したがって、ステップＳＴ２０４に進み、プロセッサは、チャネルをチャネル１からチャネル２にインクリメントする。チャネルがインクリメントされると、ステップＳＴ２０１に戻る。
【0121】
以下同様に、ステップＳＴ２０１～ＳＴ２０４のループが実行される。そして、吸収係数データｇｐ１およびｇｑ１のチャネル１～ｎに対してステップＳＴ２０１～ＳＴ２０４の処理を実行している間に、あるチャネルにおいて、吸収係数の差が閾値ＴＨ１を超えた場合、アラートを出力すると判断する。
【0122】
一方、プロセッサは、第１の基準物質の吸収係数データｇ１１およびｇ３１の全チャネルに対してステップＳＴ２０１～ＳＴ２０４を実行したが、どのチャネルにおいても吸収係数の差が閾値ＴＨ１以下である場合、ステップＳＴ２０５に進む。ステップＳＴ２０５では、プロセッサは、第１および第２の基準物質の両方の吸収係数の解析が終了したかどうか判断する。ここでは、第１の基準物質の吸収係数は解析されたが、第２の基準物質の吸収係数は解析されていない。したがって、ステップＳＴ２０１に戻る。
【0123】
ステップＳＴ２０１に戻ると、プロセッサは、第２の基準物質の吸収係数に対して、ステップＳＴ２０１～ＳＴ２０４の処理を実行する。
【0124】
図２８は、第２の基準物質の吸収係数データに対してステップＳＴ２０１～ＳＴ２０４を実行するときの説明図である。
【0125】
ステップＳＴ２０１において、プロセッサは、チャネル１における診療日ｐの吸収係数ｂｐ１と診療日ｑの吸収係数ｂｑ１との差Δｂｑ１を計算する。吸収係数の差Δｂｑ１を計算したら、ステップＳＴ２０２に進む。
【0126】
ステップＳＴ２０２では、プロセッサは、Δｂｑ１と閾値ＴＨ２とを比較し、その比較結果に基づいて、アラートを出力するかどうか判断する。例えば、Δｂｑ１が閾値ＴＨ２を超えた場合（Δｂｑ１＞ＴＨ２）、プロセッサは、吸収係数の経時的な変化が大きいので、アラートを出力すると判断する。この場合、ステップＳＴ２１に進み、プロセッサは、表示部にアラートを出力させる。
【0127】
一方、吸収係数の差Δｂｑ１が閾値ＴＨ２を超えていない場合（Δｂｑ１≦ＴＨ２）、吸収係数の経時的な変化が小さいので、プロセッサは、アラートは出力しないと判断する。この場合、ステップＳＴ２０３に進む。
【0128】
ステップＳＴ２０３では、プロセッサは、全チャネルに対してステップＳＴ２０１およびＳＴ２０２が実行されたかどうかを判断する。ここでは、チャネル１に対してステップＳＴ２０１およびＳＴ２０２が実行されたが、他のチャネル２～チャネルｎに対しては、まだ実行されていない。したがって、ステップＳＴ２０４に進み、プロセッサは、チャネルをチャネル１からチャネル２にインクリメントする。チャネルがインクリメントされると、ステップＳＴ２０１に戻る。
【0129】
以下同様に、ステップＳＴ２０１～ＳＴ２０４のループが実行される。そして、吸収係数データｈ１１およびｈｑ１のチャネル１～ｎに対してステップＳＴ２０１～ＳＴ２０４の処理を実行している間に、あるチャネルにおいて、吸収係数の差が閾値ＴＨ２を超えた場合、アラートを出力すると判断する。
【0130】
一方、プロセッサは、第２の基準物質の吸収係数データｈ１１およびｈ１３の全チャネルに対してステップＳＴ２０１～ＳＴ２０４を実行したが、どのチャネルにおいても吸収係数の差が閾値ＴＨ２以下である場合、ステップＳＴ２０５に進む。
【0131】
ステップＳＴ２０５では、プロセッサは、第１および第２の基準物質の両方の吸収係数の解析が終了したかどうか判断する。ここでは、第１および第２の基準物質の両方の吸収係数の解析が終了している。したがって、プロセッサは、アラートの出力は不要であると判断して、ステップＳＴ２０を終了する。この場合、ステップＳＴ２２に進み、診療日ｑの検査スケジュールに従って、被検体の検査が実行される。
【0132】
本実施形態では、診療日ｑでは、アラートは出力されなかったとする。したがって、ステップＳＴ２０の終了後、被検体の検査が実行される。
【0133】
以下同様に、診療日ｑ以降もフロー２００が実行される。そして、吸収係数の差が閾値を超えた場合、アラートが出力され、Ｘ線発生装置の点検が実行され、必要に応じて、Ｘ線管の修理、Ｘ線管の交換が行われる。そして、ＣＴシステムが再稼働されると、再稼働の初日に得られたキャリブレーションデータの吸収係数が新たな基準吸収係数として決定される。したがって、ＣＴシステムが再稼働した後は、再稼働の初日に得られた吸収係数を基準値として、吸収係数の経時的な変化が解析される。
【0134】
以上説明したように、本実施形態によれば、キャリブレーションスキャンを実行するたびに、吸収係数データが取得される。そして、吸収係数の差を計算することにより、吸収係数の経時的な変化量を求め、吸収係数の差が閾値（ＴＨ１又はＴＨ２）よりも大きいかどうかが判断される。閾値ＴＨ１又はＴＨ２は、例えば、Ｘ線発生装置に経年劣化が発生したかどうかを判断するのに適した値に設定することができる。したがって、Ｘ線発生装置に経年劣化が発生した場合、アラートが出力されるので、ユーザは、Ｘ線発生装置の劣化がＣＴ画像にノイズやアーチファクトとして現れる前に、Ｘ線発生装置に障害又は故障が発生している恐れがあることを事前に把握することができる。このため、ユーザは、Ｘ線発生装置の劣化がＣＴ画像にノイズやアーチファクトとして現れる前に、サービスに点検依頼するなど、必要な処置をとることができる。
【0135】
尚、第１の実施形態では、説明の便宜上、単列の検出器を取り上げて、アラートを出力するかどうかを判断する方法を説明したが、本発明は、複数列の検出器にも適用することができる。
【0136】
図２９は、複数列の検出器が使用された場合のステップＳＴ２０のフローの一例を示す図である。
図２９に示すステップＳＴ２０は、図９に示すステップＳＴ２０と比較すると、ステップＳＴ２０５とステップＳＴ２２との間に、ステップＳＴ２０５１が設けられている点が異なっている。
図２９に示すステップＳＴ２０は、図９に示すステップＳＴ２０と比較すると、ステップＳＴ２０５とステップＳＴ２２との間に、ステップＳＴ２０５１が設けられている点が異なっている。
【0137】
ステップＳＴ２０５１は、プロセッサが、検出器の全ての列に対して、ステップＳＴ２０１～ＳＴ２０５の処理が完了したかどうかを判断するステップである。ステップＳＴ２０１～ＳＴ２０５の処理が完了していない列がある場合は、ステップＳＴ１に戻る。一方、検出器の全ての列に対してステップＳＴ２０１～ＳＴ２０５の処理が完了した場合、ステップＳＴ２２に進む。
【0138】
このステップＳＴ２０５１を設けておくことにより、複数列の検出器の場合にも対応することができる。
【0139】
尚、本実施形態では、チャネル１～ｎの順序でステップＳＴ２０１～ＳＴ２０４の処理が実行されているが、必ずしもチャネル１～ｎの順序でステップＳＴ２０１～ＳＴ２０４の処理を実行する必要は無く、任意のチャネルの順序でステップＳＴ２０１～ＳＴ２０４の処理を実行することができる。また、複数列の検出器が使用されている場合、第１列～第ｍ列の順序でステップＳＴ２０１～ＳＴ２０５の処理を実行する必要は無く、任意の列の順序でステップＳＴ２０１～ＳＴ２０５を実行することができる。また、検出器の或る列の１つ以上のチャネルに対してステップＳＴ２０１およびＳＴ２０２の処理を実行した後、検出器の別の列の１つ以上のチャネルに対してステップＳＴ２０１およびＳＴ２０２の処理を実行してもよい。
【0140】
また、本実施形態では、ステップＳＴ２０において、第１の基準物質の吸収係数データを解析する場合、診療日１に取得された吸収係数ａ１１～ａ１ｎを第１の基準吸収係数として使用して、アラートを出力するかどうかを判断している。また、ステップＳＴ２０において、第２の基準物質の吸収係数データを解析する場合、診療日１に取得された吸収係数ｂ１１～ｂ１ｎを第２の基準吸収係数として使用して、アラートを出力するかどうかを判断している。しかし、診療日ごとに基準吸収係数を更新してもよい（図３０参照）。
【0141】
図３０は、診療日ごとに基準吸収係数を更新する場合のステップＳＴ２０のフローの説明図である。
図３０では、ステップＳＴ２０は、ステップＳＴ１１とステップＳＴ２０１との間に、基準吸収係数を計算するステップＳＴ２００を含んでいる。ステップＳＴ２００を設けておくことにより、診療日ごとに基準吸収係数を更新してステップＳＴ２０を実行することができる。以下に、図３０に示すステップＳＴ２０の処理の流れについて、図３１を参照しながら説明する。
図３０では、ステップＳＴ２０は、ステップＳＴ１１とステップＳＴ２０１との間に、基準吸収係数を計算するステップＳＴ２００を含んでいる。ステップＳＴ２００を設けておくことにより、診療日ごとに基準吸収係数を更新してステップＳＴ２０を実行することができる。以下に、図３０に示すステップＳＴ２０の処理の流れについて、図３１を参照しながら説明する。
【0142】
図３１は、診療日１～診療日ｉで取得された第１の基準物質の吸収係数データｇ１１～ｇｉ１を示す図である。
【0143】
ステップＳＴ１１において、診療日ｉにおける第１の基準物質の吸収係数データｇｉ１を取得したら、ステップＳＴ２０に進む。
【0144】
ステップＳＴ２０では、吸収係数の経時的な変化に基づいて、アラートを出力するかどうか判断する。
先ず、ステップＳＴ２００において、プロセッサは、診療日１～診療日ｉ－１のチャネル１の吸収係数ａ１１～ａｉ－１，１（過去に実行されたキャリブレーションスキャンにより得られた吸収係数）に基づいて、基準吸収係数ｒｉ１を計算する。この基準吸収係数ｒｉ１は、チャネル１における第１の基準吸収係数として使用される。基準吸収係数ｒｉ１は、例えば、吸収係数ａ１１～ａｉ－１，１の平均値でもよいし、吸収係数ａ１１～ａｉ－１，１の加重平均であってもよい。基準吸収係数ｒｉ１を計算したら、ステップＳＴ２０１に進む。
先ず、ステップＳＴ２００において、プロセッサは、診療日１～診療日ｉ－１のチャネル１の吸収係数ａ１１～ａｉ－１，１（過去に実行されたキャリブレーションスキャンにより得られた吸収係数）に基づいて、基準吸収係数ｒｉ１を計算する。この基準吸収係数ｒｉ１は、チャネル１における第１の基準吸収係数として使用される。基準吸収係数ｒｉ１は、例えば、吸収係数ａ１１～ａｉ－１，１の平均値でもよいし、吸収係数ａ１１～ａｉ－１，１の加重平均であってもよい。基準吸収係数ｒｉ１を計算したら、ステップＳＴ２０１に進む。
【0145】
ステップＳＴ２０１では、プロセッサは、基準吸収係数ｒ１１と診療日ｉの吸収係数ａｉ１との差Δｄｉ１を計算する。吸収係数の差Δｄｉ１を計算したら、ステップＳＴ２０２に進む。
【0146】
ステップＳＴ２０２では、プロセッサは、Δｄｉ１と閾値ＴＨ３とを比較し、その比較結果に基づいて、アラートを出力するかどうか判断する。例えば、Δｄｉ１が閾値ＴＨ３を超えた場合（Δｄｉ１＞ＴＨ３）、プロセッサは、吸収係数の経時的な変化が大きいので、アラートを出力すると判断する。この場合、ステップＳＴ２１に進み、プロセッサは、表示部にアラートを出力させる。
【0147】
一方、吸収係数の差Δｄｉ１が閾値ＴＨ３を超えていない場合（Δｄｉ１≦ＴＨ３）、吸収係数の経時的な変化が小さいので、プロセッサは、アラートは出力しないと判断する。この場合、ステップＳＴ２０３に進む。
【0148】
ステップＳＴ２０３では、プロセッサは、全チャネルに対してステップＳＴ２００～ＳＴ２０２が実行されたかどうかを判断する。ここでは、チャネル１に対してステップＳＴ２００～ＳＴ２０２が実行されたが、他のチャネル２～チャネルｎに対しては、まだ実行されていない。したがって、ステップＳＴ２０４に進み、プロセッサは、チャネルをチャネル１からチャネル２にインクリメントする。チャネルがインクリメントされると、ステップＳＴ２００に戻る。
【0149】
以下同様に、ステップＳＴ２００～ＳＴ２０４のループが繰り返し実行されるので、チャネルごとに、基準吸収係数を計算し、アラートを出力するかどうかを判断することができる。
【0150】
尚、図３１では、第１の基準物質の吸収係数の経時的な変化を解析するときに使用される第１の基準吸収係数を計算することが説明されているが、同様の方法で、第２の基準物質の吸収係数の経時的な変化を解析するときに使用される第２の基準吸収係数を計算することもできる。
【0151】
また、図３１では、診療日１～診療日ｉ－１の吸収係数ａ１１～ａｉ－１，１の全てを使用して基準吸収係数ｒｉ１を計算している。しかし、必ずしも、全ての吸収係数ａ１１～ａｉ－１，１を使用して基準吸収係数ｒｉ１を計算する必要は無く、吸収係数ａ１１～ａｉ－１，１のうちの少なくとも１つの吸収係数を使用して基準吸収係数ｒｉ１を計算してもよい。
【0152】
尚、図３１では、診療日ｉにおける基準吸収係数ｒｉ１を計算する方法が示されているが、他の診療日も、診療日ｉと同様に、過去の吸収係数に基づいて基準吸収係数を計算することができる。
【0153】
（第２の実施形態）
第１の実施形態では、Ｘ線の特性値の例として吸収係数を計算する例を説明したが、第２の実施形態では、Ｘ線の特性値の例としてＸ線スペクトルの代表値を計算する例を説明する。
第１の実施形態では、Ｘ線の特性値の例として吸収係数を計算する例を説明したが、第２の実施形態では、Ｘ線の特性値の例としてＸ線スペクトルの代表値を計算する例を説明する。
【0154】
尚、第２の実施形態について具体的に説明する前に、第２の実施形態の特徴を要約すると、以下の通りである。つまり、第２の実施形態では、キャリブレーションスキャンを実行し、キャリブレーションスキャンにより得られたデータに基づいて、高ｋＶに対応するＸ線スペクトルと、低ｋＶに対応するＸ線スペクトルを求める。そして、プロセッサは、高ｋＶに対応するＸ線スペクトルの代表値と、低ｋＶに対応するＸ線スペクトルの代表値を計算する。Ｘ線スペクトルの代表値は、例えば、エネルギー平均値、エネルギー積分値である。プロセッサは、Ｘ線スペクトルの代表値の経時的な変化に基づいて、アラートを出力するかどうかを判断する。したがって、Ｘ線スペクトルの代表値の経時的な変化が大きい場合にアラートが出力されるようにすることで、ＣＴ画像の画質が低下する前に、Ｘ線発生装置に劣化の徴候があることをユーザに知らせることができる。
【0155】
図３２は、第２の実施形態で実施されるフローの説明図である。
先ず、診療日１について説明する。診療日１では、フロー３００が実行される。
先ず、診療日１について説明する。診療日１では、フロー３００が実行される。
【0156】
ステップＳＴ５１では、Ｘ線管１０４Ａに印加される管電圧が第１の管電圧（低ｋＶ）と第２の管電圧（高ｋＶ）との間で交互に切り替えられるｋＶスイッチングを使用したキャリブレーションスキャンが実行される。キャリブレーションスキャンを実行した後、プロセッサは、キャリブレーションスキャンで得られたデータに基づいて、Ｘ線スペクトルを生成する（図３３参照）。
【0157】
図３３は、Ｘ線スペクトルの説明図である。
Ｘ線スペクトルＳＬ１は、低ｋＶに対応するＸ線スペクトルであり、Ｘ線スペクトルＳＨ１は、高ｋＶに対応するＸ線スペクトルである。Ｘ線スペクトルＳＬ１およびＳＨ１を生成した後、ステップＳＴ５２に進む。
Ｘ線スペクトルＳＬ１は、低ｋＶに対応するＸ線スペクトルであり、Ｘ線スペクトルＳＨ１は、高ｋＶに対応するＸ線スペクトルである。Ｘ線スペクトルＳＬ１およびＳＨ１を生成した後、ステップＳＴ５２に進む。
【0158】
ステップＳＴ５２では、プロセッサは、Ｘ線スペクトルＳＬ１の代表値を計算する。第２の実施形態では、プロセッサは、Ｘ線スペクトルＳＬ１の代表値として、Ｘ線スペクトルＳＬ１のエネルギー平均値ＲＬ１を計算する。また、プロセッサは、Ｘ線スペクトルＳＨ１の代表値を計算する。第２の実施形態では、プロセッサは、Ｘ線スペクトルＳＨ１の代表値として、Ｘ線スペクトルＳＨ１のエネルギー平均値ＲＨ１を計算する。しかし、Ｘ線スペクトルＳＬ１の代表値は、エネルギー平均値ＲＬ１に限定されることはなく、エネルギー積分値など、Ｘ線スペクトルＳＬ１の特徴を表す任意の値を、Ｘ線スペクトルＳＬ１の代表値として使用することができる。同様に、Ｘ線スペクトルＳＨ１の代表値は、エネルギー平均値ＲＨ１に限定されることはなく、エネルギーの積分値など、Ｘ線スペクトルＳＨ１の特徴を表す任意の値を、Ｘ線スペクトルＳＨ１の代表値として使用することができる。
【0159】
尚、診療日１に計算されたエネルギー平均値ＲＬ１は、アラートを出力するかどうかを判断するための基準となる第１のエネルギー基準値として保存される。また、診療日１に計算されたエネルギー平均値ＲＨ１は、アラートを出力するかどうかを判断するための基準となる第２のエネルギー基準値として保存される。
【0160】
ステップＳＴ５２が実行された後は、ステップＳＴ５３に進み、被検体の検査が実行される。
【0161】
次に、診療日２について説明する。診療日２では、フロー４００が実行される。
ステップＳＴ５１では、キャリブレーションスキャンが実行される。キャリブレーションスキャンを実行した後、プロセッサは、キャリブレーションスキャンで得られたデータに基づいて、Ｘ線スペクトルを生成する。図３４に、診療日２で生成された低ｋＶに対応するＸ線スペクトルＳＬ２と、診療日２で生成された高ｋＶに対応するＸ線スペクトルＳＨ２を概略的に示す。
ステップＳＴ５１では、キャリブレーションスキャンが実行される。キャリブレーションスキャンを実行した後、プロセッサは、キャリブレーションスキャンで得られたデータに基づいて、Ｘ線スペクトルを生成する。図３４に、診療日２で生成された低ｋＶに対応するＸ線スペクトルＳＬ２と、診療日２で生成された高ｋＶに対応するＸ線スペクトルＳＨ２を概略的に示す。
【0162】
ステップＳＴ５１を実行した後、ステップＳＴ６０に進む。
【0163】
ステップＳＴ６０では、Ｘ線スペクトルのエネルギー平均値の経時的な変化に基づいてアラートを出力するかどうかを判断する。
【0164】
図３５は、ステップＳＴ６０のフローの一例を示す図、図３６は、図３５のフローの説明図である。
ステップＳＴ６０１では、プロセッサは、診療日２のＸ線スペクトルＳＬ２のエネルギー平均値ＲＬ２を計算する。エネルギー平均値ＲＬ２を計算したら、ステップＳＴ６０２に進む。
ステップＳＴ６０１では、プロセッサは、診療日２のＸ線スペクトルＳＬ２のエネルギー平均値ＲＬ２を計算する。エネルギー平均値ＲＬ２を計算したら、ステップＳＴ６０２に進む。
【0165】
ステップＳＴ６０２では、プロセッサは、エネルギー基準値ＲＬ１とエネルギー平均値ＲＬ２とのエネルギー差ΔＤＬ１を計算する。エネルギー差ΔＤＬ１を計算した後、ステップＳＴ６０３に進む。
【0166】
ステップＳＴ６０３では、エネルギー差ΔＤＬ１と閾値ＴＨ５とを比較し、その比較結果に基づいて、アラートを出力するかどうか判断する。例えば、エネルギー差ΔＤＬ１が閾値ＴＨ５を超えた場合（ΔＤＬ１＞ＴＨ５）、プロセッサは、エネルギー差ΔＤＬ１が大きいので、アラートを出力すると判断する。この場合、ステップＳＴ６１に進み、プロセッサは、表示部にアラートを出力させる。
【0167】
一方、エネルギー差ΔＤＬ１が閾値ＴＨ５を超えていない場合（ΔＤＬ１≦ＴＨ５）、エネルギー差ΔＤＬ１が小さいので、プロセッサは、アラートは出力しないと判断する。この場合、ステップＳＴ６０４に進む（図３７参照）。
【0168】
図３７は、ステップＳＴ６０４～ＳＴ６０６の説明図である。
ステップＳＴ６０４において、プロセッサは、診療日２のＸ線スペクトルＳＨ２のエネルギー平均値ＲＨ２を計算する。エネルギー平均値ＲＨ２を計算したら、ステップＳＴ６０５に進む。
ステップＳＴ６０４において、プロセッサは、診療日２のＸ線スペクトルＳＨ２のエネルギー平均値ＲＨ２を計算する。エネルギー平均値ＲＨ２を計算したら、ステップＳＴ６０５に進む。
【0169】
ステップＳＴ６０５では、プロセッサは、エネルギー基準値ＲＨ１とエネルギー平均値ＲＨ２とのエネルギー差ΔＤＨ１を計算する。エネルギー差ΔＤＨ１を計算した後、ステップＳＴ６０６に進む。
【0170】
ステップＳＴ６０６では、エネルギー差ΔＤＨ１と閾値ＴＨ６とを比較し、その比較結果に基づいて、アラートを出力するかどうか判断する。例えば、エネルギー差ΔＤＨ１が閾値ＴＨ６を超えた場合（ΔＤＨ１＞ＴＨ６）、プロセッサは、エネルギー差ΔＤＨ１が大きいので、アラートを出力すると判断する。この場合、ステップＳＴ６１に進み、プロセッサは、表示部にアラートを出力させる。
【0171】
一方、エネルギー差ΔＤＨ１が閾値ＴＨ６を超えていない場合（ΔＤＨ１≦ＴＨ６）、エネルギー差ΔＤＬ１が小さいので、アラートは出力しないと判断する。この場合、低ｋＶのエネルギー平均値の経時的な変化と、高ｋＶのエネルギー平均値の経時的な変化の両方を解析したが、アラートは出力しないと判断されたので、ステップＳＴ６２に進み、診療日２の検査スケジュールに従って、被検体の検査が実行される。
【0172】
以下同様に、診療日３～診療日ｉも、診療日２と同様に、フロー４００が実行される。したがって、診療日３～診療日ｉも、診療日２と同様の手順で、Ｘ線スペクトルのエネルギー平均値の経時的な変化に基づいて、アラートを出力するかどうかが判断される。
【0173】
次に、診療日ｉについて説明する。
ステップＳＴ５１では、キャリブレーションスキャンが実行される。キャリブレーションスキャンを実行した後、プロセッサは、キャリブレーションスキャンで得られたデータに基づいて、Ｘ線スペクトルを生成する。図３８に、診療日ｉで生成された低ｋＶに対応するＸ線スペクトルＳＬｉと、診療日ｉで生成された高ｋＶに対応するＸ線スペクトルＳＨｉを概略的に示す。
ステップＳＴ５１では、キャリブレーションスキャンが実行される。キャリブレーションスキャンを実行した後、プロセッサは、キャリブレーションスキャンで得られたデータに基づいて、Ｘ線スペクトルを生成する。図３８に、診療日ｉで生成された低ｋＶに対応するＸ線スペクトルＳＬｉと、診療日ｉで生成された高ｋＶに対応するＸ線スペクトルＳＨｉを概略的に示す。
【0174】
ステップＳＴ５１を実行した後、ステップＳＴ６０に進む。
ステップＳＴ６０では、プロセッサは、Ｘ線スペクトルのエネルギー平均値の経時的な変化に基づいて、アラートを出力するかどうかを判断する。
ステップＳＴ６０では、プロセッサは、Ｘ線スペクトルのエネルギー平均値の経時的な変化に基づいて、アラートを出力するかどうかを判断する。
【0175】
図３９は、ステップＳＴ６０の説明図である。
ステップＳＴ６０１において、プロセッサは、診療日ｉのＸ線スペクトルＳＬｉのエネルギー平均値ＲＬｉを計算する。エネルギー平均値ＲＬｉを計算したら、ステップＳＴ６０２に進む。
ステップＳＴ６０１において、プロセッサは、診療日ｉのＸ線スペクトルＳＬｉのエネルギー平均値ＲＬｉを計算する。エネルギー平均値ＲＬｉを計算したら、ステップＳＴ６０２に進む。
【0176】
ステップＳＴ６０２では、プロセッサは、エネルギー基準値ＲＬｉとエネルギー平均値ＲＬｉとのエネルギー差ΔＤＬｉを計算する。エネルギー差ΔＤＬｉを計算した後、ステップＳＴ６０３に進む。
【0177】
ステップＳＴ６０３では、エネルギー差ΔＤＬｉと閾値ＴＨ５とを比較し、その比較結果に基づいて、アラートを出力するかどうか判断する。例えば、エネルギー差ΔＤＬｉが閾値ＴＨ５を超えた場合（ΔＤＬｉ＞ＴＨ５）、プロセッサは、エネルギー差ΔＤＬｉが大きいので、アラートを出力すると判断する。
【0178】
ここでは、エネルギー差ΔＤＬｉがΔＤＬｉ＞ＴＨ５を満たすとする。したがって、アラートが出力される。
【0179】
アラートが出力されると、ユーザは、例えば、サービスに点検１１（図２参照）を依頼することができる。サービスによる点検１１の結果、Ｘ線発生装置に障害が発生している又はＸ線発生装置が故障していると認められた場合、Ｘ線発生装置が修理又は交換される。ＣＴシステムは、Ｘ線発生装置の修理又は交換の後、動作チェックが行われ、正常に動作することが確認できたら、ＣＴシステムの稼働を再開する。図４０は、ＣＴシステムの稼働再開後のフローの説明図である。
ここでは、診療日ｉに点検等が完了し、次の診療日ｐにＣＴシステムの稼働が再開できたとする。
ここでは、診療日ｉに点検等が完了し、次の診療日ｐにＣＴシステムの稼働が再開できたとする。
【0180】
診療日ｐでは、診療日１と同じフロー３００が実行され、Ｘ線スペクトルのエネルギー平均値が計算される。そして、診療日ｑ以降は、診療日２と同じフロー４００が実行される。したがって、診療日ｐで計算されたエネルギー平均値を基準にして、エネルギー差が計算され、エネルギー差が閾値を超えた場合は、アラートが出力される。したがって、診療日ｑ以降も、Ｘ線スペクトルのエネルギー平均値の経時的な変化に基づいて、アラートを出力するかどうか判断することができる。
【0181】
以上説明したように、第２の実施形態によれば、キャリブレーションスキャンを実行するたびに、Ｘ線スペクトルが生成される。そして、エネルギー差を計算することにより、エネルギー平均値の経時的な変化量を求め、エネルギー差が閾値よりも大きいかどうかが判断される。したがって、エネルギー差が大きい場合、アラートが出力されるので、ユーザは、ＣＴ画像の画質が低下する前に、Ｘ線発生装置に障害又は故障が発生している恐れがあることを事前に把握することができる。このため、ユーザは、サービスにＸ線発生装置の点検を依頼するなど、必要な処置をとることができる。
【0182】
尚、第２の実施形態では、診療日１で計算されたＸ線スペクトルのエネルギー平均値をエネルギー基準値として決定し、アラートを出力するかどうかを判断している。しかし、過去に取得したＸ線スペクトルのエネルギー平均値に基づいて、診療日ごとにエネルギー基準値を計算するようにしてもよい（図４１参照）。
【0183】
図４１は、診療日ごとにエネルギー基準値を計算する場合のステップＳＴ６０のフローの説明図である。
【0184】
図４１では、ステップＳＴ６０は、ステップＳＴ６０１とステップＳＴ６０２との間に、低ｋＶのＸ線スペクトルのエネルギー基準値を計算するステップＳＴ６０１１を含んでいる。更に、ステップＳＴ６０は、ステップＳＴ６０４とステップＳＴ６０５との間に、高ｋＶのＸ線スペクトルのエネルギー基準値を計算するステップＳＴ６０４１を含んでいる。ステップＳＴ６０１１およびステップＳＴ６０４１を設けておくことにより、診療日ごとにエネルギー基準値を計算してステップＳＴ６０を実行することができる。以下に、図４１に示すステップＳＴ６０の処理の流れについて、図４２を参照しながら説明する。
【0185】
図４２は、診療日１～診療日ｉで取得された低ｋＶに対応するＸ線スペクトルＳＬ１～ＳＬｉを示す図である。
【0186】
ステップＳＴ５１において、診療日ｉの低ｋＶに対応するＸ線スペクトルＳＬｉを取得したら、ステップＳＴ６０１に進む。
【0187】
ステップＳＴ６０１では、プロセッサは、診療日ｉのＸ線スペクトルＳＬｉのエネルギー平均値ＲＬｉを計算する。エネルギー平均値ＲＬｉを計算したら、ステップＳＴ６０１１に進む。
【0188】
ステップＳＴ６０１１において、プロセッサは、診療日１～診療日ｉ－１のエネルギー平均値ＲＬ１～ＲＬｉ－１に基づいて、第１のエネルギー基準値ｖｉを計算する。このエネルギー基準値ｖｉは、診療日ｉにおいて、アラートを出力するかどうかを判断するための基準値として使用される。エネルギー基準値ｖｉは、例えば、エネルギー平均値ＲＬ１～ＲＬｉ－１の平均値でもよいし、エネルギー平均値ＲＬ１～ＲＬｉ－１の加重平均であってもよい。エネルギー基準値ｖｉを計算したら、ステップＳＴ６０２に進む。
【0189】
ステップＳＴ６０２では、プロセッサは、エネルギー基準値ｖｉとエネルギー平均値ＲＬｉとのエネルギー差ΔＲＬｉを計算する。エネルギー差ΔＲＬｉを計算した後、ステップＳＴ６０３に進む。
【0190】
ステップＳＴ６０３では、エネルギー差ΔＲＬｉと閾値ＴＨ５とを比較し、その比較結果に基づいて、アラートを出力するかどうか判断する。例えば、エネルギー差ΔＲＬｉが閾値ＴＨ５を超えた場合（ΔＲＬｉ＞ＴＨ５）、プロセッサは、エネルギー差ΔＲＬｉが大きいので、アラートを出力すると判断する。この場合、ステップＳＴ６１に進み、プロセッサは、表示部にアラートを出力させる。
【0191】
一方、エネルギー差ΔＲＬｉが閾値ＴＨ５を超えていない場合（ΔＲＬｉ≦ＴＨ５）、エネルギー差ΔＲＬｉが小さいので、プロセッサは、アラートは出力しないと判断する。この場合、ステップＳＴ６０４に進む。
したがって、アラートを出力するかどうかを判断することができる。
したがって、アラートを出力するかどうかを判断することができる。
【0192】
尚、図４２では、低ｋＶに対応するＸ線スペクトルのエネルギー平均値の経時的な変化を解析するときに使用される第１のエネルギー基準値を計算することが説明されているが、同様の方法で、高ｋＶに対応するＸ線スペクトルのエネルギー平均値の経時的な変化を解析するときに使用される第２のエネルギー基準値を計算することもできる。
【0193】
また、図４２では、診療日１～診療日ｉ－１のエネルギー平均値ＲＬ１～ＲＬｉ－１の全てを使用してエネルギー基準値ｖｉを計算している。しかし、診療日１～診療日ｉ－１のエネルギー平均値ＲＬ１～ＲＬｉ－１のうちの少なくとも１つのエネルギー平均値を使用してエネルギー基準値ｖｉを計算してもよい。
【0194】
尚、図４２では、診療日ｉにおけるエネルギー基準値ｖｉを計算する方法が示されているが、他の診療日も、診療日ｉと同様に、過去のＸ線スペクトルから得られたエネルギー平均値に基づいて、エネルギー基準値を計算することができる。
【0195】
尚、本実施形態では、医用システムとしてＣＴシステムを使用した例が示されている。しかし、本発明は、ＣＴシステムに限定されることはなく、Ｘ線源を被検体１１２に照射する医用システムであれば、ＣＴシステム以外の他のシステム（例えば、ＰＥＴ－ＣＴシステム）に適用することができる。
【符号の説明】
【0196】
１０ ＣＴシステム
１００ フロー
１０２ ガントリ
１０３ フィルタ部
１０４ Ｘ線発生装置
１０４Ａ Ｘ線管
１０４Ｂ ジェネレータ
１０６ Ｘ線
１０７ ボア
１０８ 検出器
１１６ テーブル
１１８ テーブルモータコントローラ
２００ フロー
２０２ 検出器素子
２１０ Ｘ線コントロ-ラ
２１２ ガントリモータコントローラ
２１４ ＤＡＳ
２１６ コンピュータ
２１８ 記憶装置
２２０ オペレータコンソール
２２４ ＰＡＣＳ
２３０ 画像再構成器
２３２ 表示装置
３００ フロー
４００ フロー
１００ フロー
１０２ ガントリ
１０３ フィルタ部
１０４ Ｘ線発生装置
１０４Ａ Ｘ線管
１０４Ｂ ジェネレータ
１０６ Ｘ線
１０７ ボア
１０８ 検出器
１１６ テーブル
１１８ テーブルモータコントローラ
２００ フロー
２０２ 検出器素子
２１０ Ｘ線コントロ-ラ
２１２ ガントリモータコントローラ
２１４ ＤＡＳ
２１６ コンピュータ
２１８ 記憶装置
２２０ オペレータコンソール
２２４ ＰＡＣＳ
２３０ 画像再構成器
２３２ 表示装置
３００ フロー
４００ フロー
【要約】 （修正有）
【課題】ＣＴ画像の画質が低下する前に、Ｘ線発生装置に経年劣化が発生している恐れがあることをユーザに知らせる。
【解決手段】Ｘ線発生装置と、前記Ｘ線発生装置から照射されたＸ線を検出する検出器と、１つ又は複数のプロセッサであって所定のキャリブレーションスキャンが実行されるたびに、前記検出器で検出されたＸ線のデータに基づいて、前記Ｘ線発生装置から照射されたＸ線の特性値を計算すること、前記特性値の経時的な変化に基づいて、アラートを出力するかどうか判断することを実行する１つ又は複数のプロセッサとを含む、医用システム。
【選択図】図１０
【解決手段】Ｘ線発生装置と、前記Ｘ線発生装置から照射されたＸ線を検出する検出器と、１つ又は複数のプロセッサであって所定のキャリブレーションスキャンが実行されるたびに、前記検出器で検出されたＸ線のデータに基づいて、前記Ｘ線発生装置から照射されたＸ線の特性値を計算すること、前記特性値の経時的な変化に基づいて、アラートを出力するかどうか判断することを実行する１つ又は複数のプロセッサとを含む、医用システム。
【選択図】図１０
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【図14】
【図15】
【図16】
【図17】
【図18】
【図19】
【図20】
【図21】
【図22】
【図23】
【図24】
【図25】
【図26】
【図27】
【図28】
【図29】
【図30】
【図31】
【図32】
【図33】
【図34】
【図35】
【図36】
【図37】
【図38】
【図39】
【図40】
【図41】
【図42】