JP4651190B2

JP4651190B2 - デジタルｘ線イメージング検出器のダイナミックレンジ拡張システム

Info

Publication number: JP4651190B2
Application number: JP2000399850A
Authority: JP
Inventors: ポール・リチャード・グランフォース
Original assignee: General Electric Co
Current assignee: General Electric Co
Priority date: 1999-12-30
Filing date: 2000-12-28
Publication date: 2011-03-16
Anticipated expiration: 2020-12-28
Also published as: EP1120745A3; JP2001314396A; EP1120745A2; US6721441B1

Description

【０００１】
【発明の属する技術分野】
本発明は、一般的には、デジタル検出器を使用するＸ線システムなどのデジタル・イメージング・システムの分野に関する。更に特定すれば、本発明は、様々な組織及び患者の内蔵器官の形態の均一且つ有用な画像データを提供するために、Ｘ線イメージング・システムにおけるデジタル検出器構造のダイナミックレンジを拡張する技法に関する。
【０００２】
【従来の技術】
デジタル・イメージング・システムはいくつかの分野で、特に医療診断の分野でその重要性を増しており、従来の技法より好ましい場合が多い。例えば、Ｘ線イメージングの場合、従来のシステムは写真フィルムを使用し、Ｘ線検査中にフィルムを露光する。その結果得られる画像のコントラストと詳細は患者の身体を通過するＸ線放射の様々な吸収レベルにより規定される。例えば、骨やそれに類する組織はＸ線写真では明るく見えるが、肺などの軟組織はそれより暗く見える。通常は、特に関心のある特徴のコントラストを向上させて、最も有用な画像を得るために、露光のパラメータの均衡をはかろうとする試みがなされている。デジタル・イメージング・システムでも、同様の試みが提示されている。
【０００３】
現在のデジタルＸ線イメージング・システムにおいては、Ｘ線源からの放射線を被写体、医療診断に適用する場合には通常は患者に向けて導く。放射線の一部は患者の身体を通過して、検出器に入射する。検出器の面は放射線を光子に変換し、その光子が感知される。検出器は個別の画素から成る行列に分割されており、各画素領域に入射した放射線の量、すなわち、強度に基づいて出力信号を符号化する。放射線の強さは患者の身体を放射線が通過するにつれて変化するので、出力信号に基づいて再構成される画像は、従来の写真フィルム技法により利用できる患者の組織の投射に類似する投射画像を提供する。
【０００４】
デジタルＸ線イメージング・システムは、Ｘ線技師や診断を下すべき医師が要求する画像に再構成できると共に、必要な時期が来るまでデジタル記憶又はアーカイビングしておくことができるデジタル・データを収集する能力があるという意味で、特に有用である。従来のフィルムを利用するＸ線撮影技法においては、フィルムそのものを準備し、露光し、現像した後、Ｘ線技師が使用する場合に備えてフィルムを保管しておく。特に重要である器官の詳細な情報を捕捉できるという理由によりフィルムはすぐれた診断ツールであるが、本来、撮影を行う施設又は部門から様々な医療現場へ搬送するなど、異なる場所の間でフィルムを搬送することは困難である。これに対し、直接デジタルＸ線システムにより生成されるデジタル・データは処理、画質向上、格納及びネットワークを介する送信が可能であり、デジタル・データから再構成した画像を別の場所でモニタに表示したり、その他のソフトコピー表示装置で出力することも可能である。従来のフィルムによるＸ線撮影画像をデジタル・データに変換するデジタル化システムによっても、同様の利点が得られる。
【０００５】
画像データを取得し格納し送信する上で有効であるにもかかわらず、デジタルＸ線システムはいくつかの克服すべき課題を以前としてかかえている。例えば、フィルムを利用するシステムでは、画像を見る場合にライトボックス又はそれに類するバックライトを利用することが多いが、そのような場合、再構成画像のいくつかの領域が明るすぎるか、又は暗すぎるかのいずれかの状態で再生されてしまうために、それらの領域が飽和するという問題が起こる。そのような画像における根本的な問題は、画像データのダイナミックレンジと出力装置のダイナミックレンジとが一致しないことに原因があるといえるであろう。従来のイメージ・インテンシファイア及び撮像管システムなどのアナログの分野では、ダイナミックレンジを拡張する技法が開発され、実施されている。しかし、デジタルＸ線画像を表示する際に使用すべきデジタル方式でそれらの技法を実現するのは容易ではない。
【０００６】
デジタル・イメージング検出器を使用するシステムにおける問題は、システムの様々な構成要素のデジタル・ダイナミックレンジに起因している。すなわち、デジタル検出器のダイナミックレンジは表示装置又は出力装置のダイナミックレンジとは異なる。そのため、検出器により生成される画像データのダイナミックレンジを出力装置のダイナミックレンジに適応させることが必要になると考えられる。例えば、現在のシステムの１つにおいては、データ収集回路により各画素を１２ビットのデータに符号化する。しかし、表示装置はそれとは異なるダイナミックレンジを有し、この例で言えば通常は１２ビットのレンジより小さくなる。現時点で利用できるシステムの１つにおいては、コンピュータ・モニタなどのソフトコピー表示装置はわずか８ビットの有効ダイナミックレンジを提供できるにすぎない。画像の種類、画像中で識別できる器官又は組織及びそれに類する様々な要因によって、ダイナミックレンジを適正に調整しないと、所望のコントラストが損なわれ、その結果、いくつかの領域（例えば、ダイナミックレンジの下端部）で特徴を識別することが不可能になってしまうこともある。また、画像のその他の領域の飽和も起こり、同様に、それらの領域（例えば、ダイナミックレンジの上端部）でも詳細が失われてしまう。そのような問題は特に胸部Ｘ線画像で発生する。この場合、ダイナミックレンジを適正に調整しないと、再構成画像の心臓の領域ではコントラストが低下する一方、肺の領域では飽和が発生する。
【０００７】
【発明が解決しようとする課題】
従って、デジタル・イメージング・システム、特にデジタル検出器を採用するＸ線システムの有効ダイナミックレンジを拡張するための改良された技法が必要とされている。更に、画像の特定の領域におけるコントラストの向上によって有効な画像データを提供する一方で、暗い領域でのコントラストの低下または明るい領域における飽和の発生を回避するように、画像の特定の領域に適用されるダイナミックレンジを適合させるために、プログラム可能なフィルタリングにより実現可能である改良されたシステムも現在求められている。特に、異なるダイナミックレンジを有する出力装置で再構成画像を出力させるときに、各領域のコントラストの損失及び飽和を回避するように、収集された画像データのダイナミックレンジを調整するための簡単な技法が必要とされている。
【０００８】
【課題を解決するための手段】
本発明は、デジタル・イメージング・システム、特にＸ線システムにおいて使用するのに適するダイナミックレンジ拡張技法を提供する。この技法は、デジタル検出器から収集される入力データに基づいて、システムの信号処理要素又はデジタル・フィルタリング要素を適切にプログラミングすることによって実現できる。すなわち、本発明の技法はソフトウェアのレトロフィット（遡及的修復）などにより既存のシステムでも実現できるし、また、ダイナミックレンジを拡張することが望ましい新たなシステムにも適用できる。この技法は直接デジタルＸ線イメージング・システムに特に適しているが、写真フィルムなどの従来の支援手段から符号化された画像の表示などの同様の分野にも適用可能である。更に一般的には、本発明の技法はＸ線イメージング以外の分野及び医療診断用イメージング以外の分野にも適用できるであろう。
【０００９】
本発明の技法は純粋階調(pure gray level)処理と、純粋周波数処理との均衡をとる方法を採用している。すなわち、入力データを出力データに変換するために使用されるパラメータは、純粋階調処理によりコントラストの限界を越えてしまう状況と、純粋周波数処理により「リンギング」アーティファクトが生成されてしまうという事態とを互いに相殺させるように働く。パラメータは、入力信号値の所望の閾値に対し入力信号値の関数として入力範囲にわたり変化するという非直線性を示す。階調処理制御パラメータと周波数処理制御パラメータの双方を入力値の範囲にわたって適用することにより、その結果得られるダイナミックレンジは拡張され、再構成画像の暗い領域ではコントラストが向上する一方、より明るい領域の飽和は回避される。
【００１０】
通常、この技法に従った処理により、再構成画像の明るい部分では個々の画素値、すなわち、強さは減少するが、それらの領域における詳細はある程度のレベルまで維持される。そこで、原画像データは、１２ビットのダイナミックレンジなどの第１のダイナミックレンジから表示装置の８ビット・レンジなどの所望のレンジに調整される。従って、胸部Ｘ線写真の心臓領域のような画像の暗い領域のコントラストは維持される一方、肺の詳細部におけるような画像の明るい部分の飽和が回避される。
【００１１】
【発明の実施の形態】
図１は、離散的画素画像データを収集、処理するイメージング・システム１０を概略的に示す。図示した実施例においては、イメージング・システム１０は原画像データを収集すると共に、データを表示するために本発明の技法に従って画像データを処理するように設計されたデジタルＸ線システムである。以下の説明を通して、デジタルＸ線システムについて基本的な基礎情報を与えるが、本発明による技法の様々な面を他の種類のシステムや、既存の画像のデジタル化によって原画像データ・セットを収集するようなシステムにも適用できることに留意すべきである。同様に、別途収集、格納されていた画像データを受信又はアクセスするように設計されたビューイング（視覚）ステーションにも本発明による技法の様々な面を適用できるであろう。
【００１２】
図１に示す実施例では、イメージング・システム１０は、コリメータ１４に隣接して配置されるＸ線放射線源１２を含む。コリメータ１４はＸ線放射線の流れ１６を患者１８などの被写体が位置決めされている領域へ導く。放射線の一部２０は被写体そのもの又はその周囲を通過し、全体を図中符号２２で示されるデジタルＸ線検出器に入射する。以下に更に詳細に説明するが、Ｘ線検出器２２はその面で受け取ったＸ線光子をよりエネルギーの低い光子に変換し、その後、電気信号に変換する。それらの電気信号は収集された後、被写体内部の各特徴の画像を再構成するために処理される。
【００１３】
Ｘ線放射線源１２は、検査シーケンスごとに電力と制御信号の双方を供給する給電／制御回路２４により制御される。更に、Ｘ線検出器２２は検出器制御部２６に結合している。検出器制御部２６は検出器で発生される信号の収集を指令する。検出器制御部２６は、ダイナミックレンジの初期調整のために、デジタル画像データをインタリーブするなどの様々な信号処理機能やフィルトレーション（フィルタ処理）機能を更に実行しても良い。給電／制御回路２４と検出器制御部２６は共にシステム制御部２８からの信号に応答する。通常、システム制御部２８は、検査プロトコルを実行し且つ収集した画像データを処理するためのイメージング・システムの動作を指令する。本発明に関して言えば、システム制御部２８は、一般的には汎用コンピュータ又は特定用途向けデジタル・コンピュータに基づく信号処理回路と、そのコンピュータにより実行されるプログラム及びルーチン、並びに構成パラメータ及び画像データを格納する関連メモリ回路と、インタフェース回路などを含む。図１に示す実施例では、システム制御部２８は、図中符号３０で示す表示装置又はプリンタなどの少なくとも１つの出力装置にリンクしている。出力装置は標準型コンピュータ・モニタ又は専用コンピュータ・モニタと、関連処理回路とを含んでいても良い。システム・パラメータの出力、検査要求、画像表示などのために、システムに１つ又は複数のオペレータ・ワークステーション３２を更にリンクさせても良い。通常、システム内部に設けられる表示装置、プリンタ、ワークステーション及びそれに類する装置はデータ収集のための構成要素に対してはローカルであるが、研究機関又は病院、あるいは全く別の場所などの、データ収集用構成要素から離れた場所にあって、インターネット、仮想専用ネットワークなどの１つ又は複数の構成可能ネットワークを介して画像収集システムにリンクされていても良い。
【００１４】
図２は、デジタルＸ線検出器２２の機能要素を概略的に示す図である。図２には、通常は検出器制御部２６の内部に構成されるイメージング検出器制御部（ＩＤＣ）３４も示されている。ＩＤＣ３４はＣＰＵ又はデジタル信号処理装置と、メモリ回路とを含み、検出器からの感知信号の収集を指令する。ＩＤＣ３４は２方向光ファイバ導線を介して検出器２２内部の検出器制御回路３６に結合している。これにより、動作中、ＩＤＣ３４は検出器内部で指令信号と画像データを交換する。
【００１５】
検出器制御回路３６は、図中符号３８で示される電源からＤＣ電力を受け取る。検出器制御回路３６は、システムの動作のデータ収集段階で信号を送信するために使用される行ドライバ及び列ドライバに対してタイミング指令及び制御指令を発生するように構成されている。従って、検出器制御回路３６は電力と制御信号を基準／調整器回路４０へ送信すると共に、この回路４０からデジタル画像画素データを受信する。
【００１６】
現時点で好ましい実施例では、Ｘ線検出器２２は、検査中に検出器の面で受け取られたＸ線光子をよりエネルギーの低い（光の弱い）光子に変換するシンチレータから構成されている。そこで、光検出器のアレイは光子を検出器の面の個々の画素領域に入射した光子の数、すなわち、放射線の強さを表す電気信号に変換する。以下に説明するように、電子読み取り回路は得られたアナログ信号を処理、格納することができるデジタル値に変換する。それらのデジタル値は、画像の再構成後、表示装置３０又はワークステーション３２などで表示される。現時点で好ましい実施例においては、光検出器のアレイは１つのアモルファス・シリコン基板の上に形成される。アレイ要素は行列を成し、各要素は１つのフォトダイオードと、薄膜トランジスタとから構成されている。各ダイオードの陰極はトランジスタのソースに接続し、全てのダイオードの陽極は負バイアス電圧に接続している。各行のトランジスタのゲートは互いに接続し、行電極は以下に説明する電子走査回路に接続している。１つの列のトランジスタのドレインは互いに接続し、各列の電極は電子読み取り回路に接続している。
【００１７】
図２に示す特定の実施例では、一例として、行バス４２は、検出器の様々な列からの読み取りを可能にし、且つ必要に応じて行の動作を停止し、選択された行に電荷補正電圧を印加するための複数の導線を含む。列バス４４は、行が順次イネーブル(enable)されている間に列からの読み取りを指令するための追加の導線を含む。行バス４２は、それぞれが検出器の一連の行のイネーブルを指令する一連の行ドライバ４６に結合している。同様に、電子読み取り回路４８は検出器の全ての列の読み取りを指令するために列バス４４に結合している。
【００１８】
図示した実施例では、行ドライバ４６と電子読み取り回路４８は検出器パネル５０に結合している。このパネルは複数の部分５２に分割されていても良い。各部分５２は１つの行ドライバ４６に結合し、複数の行を含む。同様に、各列ドライバ４８は一連の列に結合している。これにより、先に述べたフォトダイオードと、薄膜トランジスタから成る構造は、行５６と列５８を成して配列された一連の画素、すなわち、離散的画素５４を規定する。これらの行と列は、高さ６２と幅６４を有する画像行列６０を規定する。
【００１９】
更に、図２に示すように、各画素５４は通常は１つの行と１つの列との交差部分に規定されており、その交差部分において列電極６８は行電極７０と交差する。先に述べたように、画素ごとに、各々の交差場所に１つの薄膜トランジスタ７２が設けられ、フォトダイオード７４も同様である。行ドライバ４６により各行がイネーブルされると、各々のフォトダイオードからの信号は電子読み取り回路４８によりアクセスされ、その後の処理及び画像再構成に備えてデジタル信号に変換される。
【００２０】
図３は、図２に概略的に示した構成要素の物理的配置の一例を示す。図３に示すように、検出器はガラス基板７６を含み、その上に以下に説明する構成要素が配置されている。列電極６８と行電極７０は基板上に設けられ、先に説明した薄膜トランジスタと、フォトダイオードとを含むアモルファス・シリコン・フラットパネル・アレイ７８が規定されている。このアモルファス・シリコン・アレイの上に、先に説明したように検査シーケンス中に放射線を受け取るシンチレータ８０が設けられている。列電極及び行電極との間で信号を通信するためのコンタクト・フィンガ８２が形成され、コンタクト・フィンガと外部回路との間で信号を通信するためのコンタクト・リード８４が設けられている。
【００２１】
先に説明したようなデジタルＸ線イメージング／表示システムにおいては、表示装置３０又はオペレータ・ワークステーション３２で表示された画像を見るとき、又はその他の方法で画像を出力するときに、出力画像のダイナミックレンジがシステムの画像データ収集要素のダイナミックレンジと異なるという特定の問題が起こる。例えば、このデジタルＸ線イメージング・システムの実施例では、画像データ収集要素は１２ビットのダイナミックレンジを有するが、表示装置のダイナミックレンジはそれより少ない８ビットである。そのため、収集されて格納されている画像データの場合、画像行列の各画素は長さ１２ビットの一連の２進値により表現される。再構成画像においては、表示装置又は出力装置のダイナミックレンジについて適切な修正又は対応を行わないと、再構成画像のある領域が非常に暗く見える一方で、他の領域は非常に明るく見えるという状況が起こる。いずれの場合にも、所望のコントラストや識別可能な詳細部が失われるか、あるいはユーザによる識別が不可能になってしまう。この事態は、肺の領域は非常に明るく映るが、暗い領域、特に心臓の領域のコントラストが低下する胸部Ｘ線写真の場合にとりわけ顕著に現れる。
【００２２】
本発明は、出力装置又は表示装置に対応するために、イメージング・システムにおいて収集又は格納された画像データを調整するメカニズムを提供する。現時点で考えられている実施例においては、この技法は、システム制御部２８、オペレータ・ステーション３２（図１を参照）、又は画像データ収集システム又は画像表示システムの同様の構成要素の内部に格納された適切なプログラミング・コードに基づいて実行される。プログラミング・コードは画像データ収集ルーチン又は画像表示ルーチンの一部であっても良いし、あるいは指令に応じて自動的に実行されても良く、また、レントゲン技師がワークステーション３２で画像を見る場合のようにオペレータの要求に応じて実行されても良い。更に、本発明の技法を実現するプログラミング・コードは、磁気記憶装置、光記憶装置、ランダム・アクセス・メモリ、読み取り専用メモリなどの適切な媒体で支援されていても良い。また、構成可能ネットワーク又はオープン・ネットワークなどの何らかの適切なプロセスを介して、あるいは専有ネットワークを介して、コードを機械で読み取り可能な媒体にロードしても良い。
【００２３】
本発明の技法によれば、原画像データを後続するいくつかの過程を経て処理し、表示又は画像再構成に適する出力データを得る。通常、この技法は画像の明るい部分における画素値を減少させ、画像のより暗い部分のコントラストの低下又は飽和を回避する。現時点で好ましいアルゴリズムは、アンシャープ(unsharp) マスキング及び輝度制御のために２組のパラメータ値を制御する。その結果、輝度は、コントラストが厳密に制限される純粋階調処理と、「リンギング」アーティファクトが生成される可能性のある純粋周波数処理との間で変化することになる。
【００２４】
そこで、好ましい処理技法をより詳細に説明すると、プロセスの第１のステップとして、システム制御部２２のメモリ又はアーカイブ、あるいはその他の記憶装置から原画像データをアクセスする。尚、この原画像データ・セットは生データから構成されていても良いし、一部処理済のデータであっても良い。データは、画像データ収集要素及び処理要素のダイナミックレンジにわたり変化する値を有する、画像行列中の各画素を表すコードから構成される。それらの値は、関係式
Out＝GAMMALUT(In-BOOSTLUT(Smooth)*Smooth) （式１）
に従って出力値に変換される。式中、「Ｉｎ」は入力画素の値、「Ｏｕｔ」は同じ画素の出力値である。式１のその他のパラメータについては、以下に更に詳細に説明する。
【００２５】
式１により示唆される処理は、自動輝度制御ウィンドウ又はそれに類するユーザ調整技法によって平均値が固定されるように画素値がスケーリング（拡縮）されているような入力画像データ・セット、すなわち、初期画像データ・セットに基づいているといえる。この初期データ・セットをボックスカー平滑化(boxcar smoothing)技法により平滑化して、式１では「Ｓｍｏｏｔｈ」と表されている１組の画像値を得る。当業者には理解されるであろうが、この平滑化演算は一般には低域フィルタリングプロセスとして機能し、そのプロセスにおいては、ここの画素の元の値と、隣接する画素の値から取り出された新たな値が画素値に割り当てられる。この技法では、適切な何らかの平滑化手続きを採用すれば良く、ボックスカー・フィルタのサイズも適切に設定されれば良い。
【００２６】
式１では「ＢＯＯＳＴＬＵＴ」及び「ＧＡＭＭＡＬＵＴ」と表されているパラメータ値又はパラメータ・セットは、この技法においては、アンシャープ・マスキング及び輝度制御のための演算を実行するために使用される。すなわち、「ＢＯＯＳＴＬＵＴ」は、アンシャープ・マスキングの量を信号レベルの関数として制御するルックアップ・テーブルから構成される。以下に説明するように、この実施例においては、「ＢＯＯＳＴＬＵＴ」の関数は、閾値画素レベルより低いときはアンシャープ・マスキング値は０であり（すなわち、この閾値より低いときはアンシャープ・マスキングは実行されない）、閾値を越えると単調に増加するような関数であるのが好ましい。パラメータ・セット「ＧＡＭＭＡＬＵＴ」は、画像の明るい部分をロールオフする同様のルックアップ・テーブルである。この実施例では、「ＧＡＭＭＡＬＵＴ」の値は定義済み閾値より低いときは線形関数を定義し、閾値を越えると、ｙ＝（ｘのγ乗）の出力対入力関係を有する。尚、γは０と１の間の値である。
【００２７】
先に述べた通り、「ＢＯＯＳＴＬＵＴ」及び「ＧＡＭＭＡＬＵＴ」の値は、純粋階調処理と純粋周波数処理との均衡をとるように決定するのが好ましい。特定のシステムでは一般にそのような確定を実験により要求するが、それらのパラメータの値の例を図４及び図５に画素値の関数として示す。すなわち、図４は、画素値を水平軸９２に沿って定義し、「ＢＯＯＳＴＬＵＴ」値を垂直軸９４に沿って定義して、一連の「ＢＯＯＳＴＬＵＴ」値又は関係を示す。一連の線９６、９８及び１００はそれらの入力値と出力値との関係を定義する。式１に従えば、「ＢＯＯＳＴＬＵＴ」の選択された値は処理の初期段階で生成された平滑化画素値に基づいているのが好ましいことに留意すべきである。
【００２８】
線９６は、プロトタイプシステムで当初に採用された値を表す。プロトタイプシステムでは、画像データを平滑化するために、２５×２５から１５１×１５１までの平滑化カーネル・サイズを使用した。「ＢＯＯＳＴＬＵＴ」の値は、この値のアンシャープ・マスキングが１００カウントの閾値（平滑画素強さ値）より低いときは０であり、閾値を越えると１．８×１０^-3の傾きで直線的に増加するように規定された。
【００２９】
線９８は、蛍光画像の場合に使用される別の「ＢＯＯＳＴＬＵＴ」値を表す。
この例では、「ＢＯＯＳＴＬＵＴ」の値は１００カウントに満たない画素値については０レベルにとどまり、その後、図４に図中符号１０４で示すように、５×１０^-4の傾きで直線的に増加する。この第２の演算においては、７５×７５の平滑化カーネルを使用した。
【００３０】
図４の線１００は、胸部Ｘ線写真の肺の部分の領域が暗くならないようにするために使用される「ＢＯＯＳＴＬＵＴ」の第３の代替関数を示す。この実施例では、「ＢＯＯＳＴＬＵＴ」の値は、画素値が大きいところではアンシャープ・マスキングの強さを減少させる一方、適度な画素値のところでは強さを増すように修正されている。すなわち、この場合にも、「ＢＯＯＳＴＬＵＴ」関数は１００カウントの閾値を有し、それに満たない場合には「ＢＯＯＳＴＬＵＴ」値は０である。しかし、閾値を越えると、「ＢＯＯＳＴＬＵＴ」の値は初めの増加領域１０６においては相対的に速い速度で増加し、その後、遷移領域１０８に至る。その後、「ＢＯＯＳＴＬＵＴ」の値は、図４に図中符号１１０で示すようにゆるやかな速度で、画素値の関数として増加する。
【００３１】
図５は、本発明の技法において、式１により定義される関係を適用する場合に使用される「ＧＡＭＭＡＬＵＴ」の値の例を示す。図中の値は、図５に図中符号１１２で示されている関係により規定されており、このグラフの場合にも、画素値は水平軸１１４に沿って表され、「ＧＡＭＭＡＬＵＴ」値は垂直軸１１６に沿って示されている。
【００３２】
最初のプロトタイプでは、図４の線９６により定義される「ＢＯＯＳＴＬＵＴ」値に対応して、「ＧＡＭＭＡＬＵＴ」値は下記の関係式により定義される。
【００３３】
ｘ＜ｂ／ｍの場合：ｙ＝ｍｘ
ｂ／ｍ＜ｘ＜２５５の場合：
【００３４】
【数１】

【００３５】
（式２）
式中、ｙは「ＧＡＭＭＡＬＵＴ」の値、ｘは入力画素値、ｍ及びｂはそれぞれ１と１４０の値に設定された実験パラメータ、γは、９０の表示最適カウントの９倍で２５５の出力カウントに到達するように調整された値である。図からわかるように、関数はｂ／ｍの閾値より低いときは傾きｍをもって直線的に増加し、滑らかにロールオフして、γによって決まるｘの値で２５５の最大値に達する。
【００３６】
「ＢＯＯＳＴＬＵＴ」の線９８に対応する次の実施例では、「ＧＡＭＭＡＬＵＴ」の値は図５の線１１８により表されるように規定されている。この実施例においては、初期領域１２２は初期速度で画素値に伴って増加し、この増加速度は後続領域１２４に入ると徐々に低下している。
【００３７】
更に、図５は、図４の線１００により「ＢＯＯＳＴＬＵＴ」の値が示されているような実施例において採用される「ＧＡＭＭＡＬＵＴ」の値の例も示している。この実施例では、「ＧＡＭＭＡＬＵＴ」の値は式２により定義される元来の「ＧＡＭＭＡＬＵＴ」の関係に従っているが、図５で図中符号１２６により示される初めの線形部分の傾きはわずかに小さくなり、表示信号レベルを約２％減少させている（すなわち、式２においてｍは０．９７８に変更されている）。
【００３８】
上述の技法を一連の異なる患者に対して適用した。それぞれの患者は、通常、デジタルＸ線イメージング・システムで撮影された２つのレコード（記録）シーケンスと、２つの蛍光透視画像シーケンスを有していた。各シーケンスには数百枚の画像が含まれる。上述の技法に従って処理した画像をＸ線技師による無作為比較試験の基礎資料として使用した結果、画像の明るい領域と暗い領域において適切な強さレベルを維持すると共に、その一方で詳細部を維持するという課題に対し、すぐれた成果を示した。
【００３９】
本発明の特定の実施例を例として添付の図面に示し且つ詳細に説明したが、本発明は様々な変形及び代替構成をとることが可能である。本発明がここで開示した特定の形態に限定されるものでないということを理解すべきである。本発明は特許請求の範囲により定義される本発明の趣旨に含まれる全ての変形、等価物及び代替構成を包含する。
【図面の簡単な説明】
【図１】本発明による技法を取り入れたデジタルＸ線イメージング・システムの概略図。
【図２】再構成に備えて画像データを生成するために図１のシステムの検出器において画像データを生成する機能回路の一部の略図。
【図３】画像データを生成する検出器構造の一例を示す部分断面図。
【図４】本発明の技法による画像データのアンシャープ・マスキングのための数値の一例を示すグラフ。
【図５】本発明の技法によるアンシャープ・マスキングされた画素値の輝度制御のための数値を示すグラフ。
【符号の説明】
１０イメージング・システム
１２Ｘ線放射線源
２２デジタルＸ線検出器
３０表示装置／プリンタ
３２オペレータ・ワークステーション
５４画素
６０画像行列
９６，９８，１００一連の線
１１８線

Claims

離散的画素イメージング・システムでデジタル画像データを処理する方法において、
画像行列（６０）における対応する複数の画素位置（５４）に対する複数の入力データ値を含む入力データのセットをアクセスする過程と、
関係式 Out＝GAMMALUT(In-BOOSTLUT(Smooth)*Smooth)に従って、
１）前記画像行列（５０）の隣接する入力データ値によって前記入力データ値を平滑化する過程と、
２）前記平滑化データ値、前記入力データ値及びアンシャープ・マスキング・パラメータ値のセット（１００）に基づいて画素ごとにアンシャープ・マスキング値を生成する過程と、
３）前記アンシャープ・マスキング値及び輝度制御パラメータのセット（１１８、１２０）に基づいて画素ごとに出力値を生成する過程と、を有し、
ここで、
Ｉｎは、前記入力データ値、
Ｏｕｔは、前記出力値、
Smoothは、前記平滑化データ値、
BOOSTLUT()は、前記アンシャープ・マスキング・パラメータ値のセット（１００）に対応する信号レベルの関数、
GAMMALUT()は、前記輝度制御パラメータのセット（１１８、１２０）に対応する関数であり、
前記アンシャープ・マスキング・パラメータ（１００）は、
閾値平滑化データ値より小さいときはほぼ０に等しい値（１０２）を有し、
画素値が前記閾値平滑化データ値を越えた後、相対的に速い速度で増加する増加領域（１０６）と、その後の遷移領域（１０８）と、その後、相対的にゆるやかな速度で増加する増加領域（１１０）に対応する値を有し、
前記遷移領域（１０８）が、前記相対的に速い速度で増加する増加領域（１０６）と、前記相対的にゆるやかな速度で増加する増加領域（１１０）とを連続する曲線により接続する、方法。
前記入力データ値は第１のダイナミックレンジを有し、前記出力値は前記第１のダイナミックレンジより小さい第２のダイナミックレンジを有し、且つ前記アンシャープ・マスキング・パラメータ（１００）及び前記輝度制御パラメータ（１１８、１２０）は前記第１のダイナミックレンジを前記第２のダイナミックレンジへマッピングする請求項１記載の方法。
前記第１のダイナミックレンジは１２ビットのレンジであり、前記第２のダイナミックレンジは８ビットのレンジである請求項２記載の方法。
前記入力データは、Ｘ線放射線源（１２）からの放射線を受け取り且つ前記画像行列内の各画素の場所で受け取られた放射線に基づいて画素ごとに前記入力データ値を生成するように構成されたデジタルＸ線検出器（２２）から収集される請求項１乃至３のいずれかに記載の方法。
前記デジタル画像データが、胸部Ｘ線写真である請求項４に記載の方法。
デジタル画像データから再構成された画像を表示するシステムにおいて、
画像（６０）の画素（５４）に対して第１のダイナミックレンジを有する画像データ・セットを複数の入力データ値を含む入力データのセットとして格納するデータ記憶装置と、
前記第１のダイナミックレンジとは異なる第２のダイナミックレンジを有する表示装置（３０）と、
関係式 Out＝GAMMALUT(In-BOOSTLUT(Smooth)*Smooth)に従って、前記データ記憶装置から画像データ・セットをアクセスし、前記画像データ・セットの画素に対し前記入力データ値を平滑化し、前記平滑化データ値からアンシャープ・マスキング値を生成し、且つ前記アンシャープ・マスキング値から前記第２のダイナミックレンジに適合する出力値を生成するように構成されたデータ処理装置（３２）と、を具備し、
ここで、
Ｉｎは、前記入力データ値、
Ｏｕｔは、前記出力値、
Smoothは、前記平滑化データ値、
BOOSTLUT()は、前記アンシャープ・マスキング・パラメータ値に対応する信号レベルの関数、
GAMMALUT()は、前記輝度制御パラメータのセット（１１８、１２０）に対応する関数であり、
前記アンシャープ・マスキング・パラメータ値は、
閾値平滑化データ値より小さいときはほぼ０に等しい値（１０２）を有し、
画素値が前記閾値平滑化データ値を越えた後、相対的に速い速度で増加する増加領域（１０６）と、その後の遷移領域（１０８）と、その後、相対的にゆるやかな速度で増加する増加領域（１１０）に対応する値を有し、
前記遷移領域（１０８）が、前記相対的に速い速度で増加する増加領域（１０６）と、前記相対的にゆるやかな速度で増加する増加領域（１１０）とを連続する曲線により接続する、システム。
前記出力値は前記第１のダイナミックレンジより小さい前記第２のダイナミックレンジを有し、且つ前記アンシャープ・マスキング・パラメータ値及び前記輝度制御パラメータ（１１８、１２０）は前記第１のダイナミックレンジを前記第２のダイナミックレンジへマッピングする請求項６記載のシステム。
前記第１のダイナミックレンジは１２ビットのレンジであり、前記第２のダイナミックレンジは８ビットのレンジである請求項７記載のシステム。
前記画像データ・セットを生成する画像収集システム（１０）と、
Ｘ線源（１２）と、
デジタルＸ線検出器（２２）とを含む請求項６乃至８のいずれかに記載のシステム。
前記デジタル画像データが、胸部Ｘ線写真である請求項９に記載のシステム。