(57)【要約】コンピュータ制御Ｘ線撮影画像のモザイクイメージの露光不足領域の強度を補正する方法および装置を提供する。本発明では、ＣＲイメージのセットを取得し（ステップ５０）、ＣＲイメージを処理し（ステップ５１）、イメージ中の露光不足領域を検出する（ステップ５２）。露光不足領域が検出されるとそのＣＲイメージは露光不足部分と正常露光部分とに分離される（ステップ５３）。ついで露光不足部分の強度が正常露光部分の強度にほぼ等しくなるように、またはできる限り近似するように調整される（ステップ５４）。これは非線形強調および／または線形変換の技術を用いて行われる。そののち強度補正された各ＣＲイメージはモザイクイメージとして結合される（ステップ５６）。
2005505960.tif 000002
【特許請求の範囲】
【請求項１】
ＣＲイメージを入力するステップと、
ＣＲイメージを正常露光部分と露光不足部分とに分離するステップと、
露光不足部分を強調するステップと、
強調部分および正常露光部分を用いて強度補正パラメータを求めるステップと、
強度補正パラメータを用いて強調部分のイメージ強度を調整するステップとを有する
ことを特徴とするＣＲイメージの強度の調整方法。
【請求項２】
露光不足部分を強調するステップにおいて、非線形変換を露光不足部分に適用して強度変化のダイナミックレンジを増大する、請求項１記載の方法。
【請求項３】
イメージ強度を調整するステップにおいて、求められた強度補正パラメータを用いて線形変換を強調部分に適用する、請求項１記載の方法。
【請求項４】
強度補正パラメータを求めるステップにおいて、強調部分および正常露光部分の同じ大きさのサンプルに対して線形回帰変換を行う、請求項１記載の方法。
【請求項５】
露光不足部分の強度を調整したＣＲイメージのセットをモザイクイメージとして結合するステップをさらに有する、請求項１記載の方法。
【請求項６】
ＣＲイメージを入力するステップと、
ＣＲイメージを正常露光部分と露光不足部分とに分離するステップと、
露光不足部分のイメージ強度を強調するステップと、
強調部分および正常露光部分を用いて第１の強度補正パラメータのセットを求めるステップと、
露光不足部分および正常露光部分を用いて第２の強度補正パラメータのセットを求めるステップと、
第１の強度補正パラメータのセットおよび第２の強度補正パラメータのセットを評価して最適な強度補正を行える強度補正パラメータのセットを選択するステップと、
選択された強度補正パラメータを用いて露光不足部分のイメージ強度を調整するステップとを有する
ことを特徴とするＣＲイメージの強度の調整方法。
【請求項７】
露光不足部分を強調するステップにおいて、非線形変換を露光不足部分に適用して強度変化のダイナミックレンジを増大する、請求項６記載の方法。
【請求項８】
第１の強度補正パラメータのセットを求めるステップにおいて、強調部分および正常露光部分の同じサイズのサンプルに対して線形回帰変換を行う、請求項６記載の方法。
【請求項９】
第２の強度補正パラメータのセットを求めるステップにおいて、露光不足部分および正常露光部分の同じサイズのサンプルに対して線形回帰変換を行う、請求項６記載の方法。
【請求項１０】
第１の強度補正パラメータのセットおよび第２の強度補正パラメータのセットを評価するステップにおいて、どちらのパラメータセットが最小残留誤差を生じるかを求め、最小残留誤差を生じるほうのセットを選択する、請求項９記載の方法。
【請求項１１】
イメージ強度を調整するステップにおいて、第１の強度補正パラメータのセットが選択された場合、第１の強度補正パラメータのセットを用いて強調部分に線形変換を適用する、請求項１０記載の方法。
【請求項１２】
イメージ強度を調整するステップにおいて、第２の強度補正パラメータのセットが選択された場合、第２の強度補正パラメータのセットを用いて露光不足部分に線形変換を適用する、請求項１０記載の方法。
【請求項１３】
露光不足部分のイメージ強度を調整したＣＲイメージのセットをモザイクイメージとして結合するステップをさらに有する、請求項６記載の方法。
【請求項１４】
ＣＲイメージの強度の調整方法を実行する命令プログラムを有形に備えた機械で読み出し可能なプログラム記憶装置において、
ＣＲイメージを入力するステップの実行命令と、
ＣＲイメージを正常露光部分と露光不足部分とに分離するステップの実行命令と、
露光不足部分を強調するステップの実行命令と、
強調部分および正常露光部分を用いて強度補正パラメータを求めるステップの実行命令と、
強度補正パラメータを用いて強調部分のイメージ強度を調整するステップの実行命令とを含む
ことを特徴とする機械で読み出し可能なプログラム記憶装置。
【請求項１５】
露光不足部分を強調するステップの実行命令は、非線形変換を露光不足部分に適用して強度のバリエーションのダイナミックレンジを増大する命令を含む、請求項１４記載の装置。
【請求項１６】
イメージ強度を調整するステップの実行命令は、求められた強度補正パラメータを用いて線形変換を強調部分に適用する命令を含む、請求項１４記載の装置。
【請求項１７】
強度補正パラメータを求めるステップの実行命令は、強調部分および正常露光部分の同じサイズのサンプルに対して線形回帰変換を行う命令を含む、請求項１４記載の装置。
【請求項１８】
露光不足部分の強度を調整したＣＲイメージのセットをモザイクイメージへ結合するステップの実行命令をさらに含む、請求項１４記載の装置。
【請求項１９】
ＣＲイメージの強度の調整方法を実行する命令プログラムを有形に備えた機械で読み出し可能なプログラム記憶装置において、
ＣＲイメージを入力するステップの実行命令と、
ＣＲイメージを正常露光部分と露光不足部分とに分離するステップの実行命令と、
露光不足部分のイメージ強度を強調するステップの実行命令と、
強調部分および正常露光部分を用いて第１の強度補正パラメータのセットを求めるステップの実行命令と、
露光不足部分および正常露光部分を用いて第２の強度補正パラメータのセットを求めるステップの実行命令と、
第１の強度補正パラメータのセットおよび第２の強度補正パラメータのセットを評価して最適な強度補正を行える強度補正パラメータのセットを選択するステップの実行命令と、
選択された強度補正パラメータを用いて露光不足部分のイメージ強度を調整するステップの実行命令とを有する
ことを特徴とする機械で読み出し可能なプログラム記憶装置。
【請求項２０】
露光不足部分を強調するステップの実行命令は、非線形変換を露光不足部分に適用して強度変化のダイナミックレンジを増大する命令を含む、請求項１９記載の装置。
【請求項２１】
第１の強度補正パラメータのセットを求めるステップの実行命令は、強調部分および正常露光部分の同じサイズのサンプルに対して線形回帰変換を行う命令を含む、請求項１９記載の装置。
【請求項２２】
第２の強度補正パラメータのセットを求めるステップの実行命令は、露光不足部分および正常露光部分の同じサイズのサンプルに対して線形回帰変換を行う命令を含む、請求項２１記載の装置。
【請求項２３】
第１の強度補正パラメータのセットおよび第２の強度補正パラメータのセットを評価するステップの実行命令は、どちらのパラメータセットが最小残留誤差を生じるかを求め、最小残留誤差を生じるほうのセットを選択する命令を含む、請求項２２記載の装置。
【請求項２４】
イメージ強度を調整するステップの実行命令は、第１の強度補正パラメータのセットが選択された場合、第１の強度補正パラメータのセットを用いて強調部分に線形変換を適用する命令を含む、請求項２３記載の装置。
【請求項２５】
イメージ強度を調整するステップの実行命令は、第２の強度補正パラメータのセットが選択された場合、第２の強度補正パラメータのセットを用いて露光不足部分に線形変換を適用する命令を含む、請求項２３記載の装置。
【請求項２６】
露光不足部分のイメージ強度を調整したＣＲイメージのセットをモザイクイメージとして結合するステップの実行命令をさらに含む、請求項１９記載の装置。
【発明の詳細な説明】
【技術分野】
【０００１】
関連出願へのクロスリファレンス
この出願は２００１年１月１０日付提出の米国出願第６０／３２６２３３号明細書の内容に関連している。
【０００２】
技術分野
本発明はコンピュータ制御放射線撮影法（ＣＲイメージング法）に関する。本発明は特に非線形強調技術および／または線形変換技術を用いてＣＲイメージ内の露光不足領域の強度を補正する方法および装置に関する。
【０００３】
発明の背景
ＣＲイメージング法では、リンのストレージプレートを備えた複数のカセットが隣接するものとそれぞれ重なり合うように配置され、１つのプレートでイメージング可能な領域よりも大きな領域を分析できるようになっている。図１には複数のカセットを用いて複数のＣＲイメージを取得する方法が示されている。イメージングすべきオブジェクト１０（例えば人体脚部）はｘ線源１１と複数のイメージ記録カセット１２ａ，１２ｂ，１２ｃのあいだの観察野内に配置される。イメージ記録カセットはそれぞれオブジェクト１０の部分ごとの高分解能画像を取得するように相互に重なり合っているので、オブジェクト１０の完全な記録が保証される。のちに各カセットのＣＲイメージはモザイクイメージとして結合される。こうしたイメージ取得の利点として例えば、人体の動きによるアーティファクトが発生せず、連続する各イメージのアライメントを完璧に行えることが挙げられる。
【０００４】
モザイクイメージのコンポーズ法、すなわち個々のＣＲイメージをシームレスに結合する手段は周知である。例えば２０００年８月に登録された米国特許第６１０１２３８号明細書“System for generating a compound x-ray image for diagnosis”には相互相関プロセスに基づいた手段が記載されている。また１９９９年１１月１６日に登録された米国特許第５９８６２７９号明細書“Method of recording and reading a radiation image of an elongate body”には基準グリッド線を用いる手段が記載されている。ただしこれらの従来技術ではＣＲイメージのモザイクコンポーズの際にアーティファクトが発生してしまう。
【０００５】
より詳細には、従来の手段を用いたモザイクＣＲイメージにアーティファクトが発生する様子が図２に示されている。図２には重なり合う２つのカセット１２ａ，１２ｂが示されており、カセット１２ｂの捉えたイメージ内に露光不足領域が生じている。隣接するカセット１２ａ，１２ｂが重なり合っており、ｘ線源から見るとカセット１２ｂはカセット１２ａの後方に存在しているため、カセット１２ｂで捉えられたイメージに露光不足部分（ホワイトバンド部分）が生じている。後方に隠れたカセット１２ｂの領域は露光のための線量が少なく、カセット１２ｂのイメージには露光不足部分１５が生じるのである。ホワイトバンド部分すなわち露光不足部分１５と正常露光部分１７とが区別される。
【０００６】
これに加えて、カセット１２ａに厚い金属エッジ１６が存在するために、カセット１２ｂの露光不足領域１５の一部がカセット１２ｂの金属エッジに相応にモザイクイメージ１３中にアーティファクト１４として現れている。これをホワイトストラップと称する。
【０００７】
ＣＲイメージ内に露光不足領域が存在するとこれがイメージを検討する医師らに視覚的な障害となるばかりでなく、重要な診断情報を隠してしまうおそれも生じる。したがって露光不足領域のイメージ強度を制御する手段を設け、この領域がイメージの他の部分と同様に明確に表現されるようにすることが望ましい。こうした手段があれば露光不足領域の情報の分析も保証され、診断の際に良好な視認性が得られる。
【０００８】
発明の概要
本発明の課題は、ＣＲイメージ内の露光不足領域または“ホワイトバンド”領域のイメージ強度を自動補正する方法を提供し、診断情報に歪みを生じることなくＣＲイメージからアーティファクトを除去ないし低減することである。
【０００９】
本発明の第１の特徴はモザイクＣＲイメージを形成する方法である。ここではＣＲイメージのセットが取得され、各ＣＲイメージが処理されてＣＲイメージ内の露光不足領域が検出される。露光不足領域が検出された場合、当該のＣＲイメージは露光不足部分と正常露光部分とに分離される。ついで露光不足部分のイメージ強度が線形変換技術および／または非線形強調技術を用いて、正常露光部分の強度とほぼ等しくなるように、または正常露光部分の強度にできるだけ近似するように調整される。強度補正された各ＣＲイメージはモザイクイメージとして結合される。
【００１０】
本発明の第２の特徴はＣＲイメージの強度を調整する方法である。ここではＣＲイメージが入力され、露光不足部分と正常露光部分とが分離され、ついで露光不足部分が強調される。有利には、露光不足部分の強調プロセスにおいて露光不足部分に非線形変換が適用され、この領域の強度変化のダイナミックレンジが増大される。さらに強調部分および正常露光部分を用いて強度補正パラメータのセットが求められる。有利には強度補正パラメータは強調部分および正常露光部分の同じサイズのサンプルに線形回帰変換を行うことにより求められる。強調部分のイメージ強度はこの強度補正パラメータを用いて調整される。有利には、イメージ強度の調整プロセスにおいて、求められた強度補正パラメータを用いて強調部分に線形変換が適用される。
【００１１】
本発明の第３の特徴はＣＲイメージの強度を調整する方法のバリエーションである。ここでは露光不足領域の強調を行うか否かが自動的に判別され、過度の強調が防止される。まずＣＲイメージが入力され、露光不足部分と正常露光部分とが分離され、露光不足部分が強調される。露光不足部分の強調プロセスにおいては露光不足部分に非線形変換が適用されてこの領域の強度変化のダイナミックレンジが増大される。ついで強調部分および正常露光部分を用いて第１の強度補正パラメータのセットが求められる。これは有利には強調部分および正常露光部分の同じサイズのサンプルに線形回帰変換を行うことにより求められる。さらに露光不足部分および正常露光部分を用いて第２の強度補正パラメータのセットが求められる。これは有利には露光不足部分および正常露光部分の同じサイズのサンプルに線形回帰変換を行うことにより求められる。続いて当該の第１のセットと第２のセットとが評価され、最適な強度補正値の得られるほうのセットが選択される。有利には、評価のプロセスにおいてどちらのパラメータセットから最小残留誤差が生じるかが求められ、最小残留誤差を生じるほうのセットが選択される。第１の強度補正パラメータのセットが選択された場合には、露光不足部分のイメージ強度は有利にはこの第１のセットを用いて強調部分に線形変換を適用することにより調整される。第２の強度補正パラメータのセットが選択された場合には、露光不足部分のイメージ強度は有利にはこの第２のセットを用いて露光不足部分に線形変換を適用することにより調整される。
【００１２】
本発明のその他の態様、対象、構成および利点は添付図に則して以下に詳細に説明する有利な実施例から得られる。
【００１３】
図面の簡単な説明
図１には従来のＣＲイメージング法が示されている。図２にはＣＲイメージの露光不足領域に起因してモザイクＣＲイメージ内に発生するアーティファクトの例が示されている。図３にはＣＲイメージ内の露光不足領域の強度を補正する本発明の装置の第１の実施例が示されている。図４にはＣＲイメージ内の露光不足領域の強度を補正する本発明の装置の第２の実施例が示されている。図５には強度補正を行うモザイクＣＲイメージング法のフローチャートが示されている。図６にはＣＲイメージ内の露光不足領域の強度補正方法の第１の実施例のフローチャートが示されている。図７にはＣＲイメージ内の露光不足領域の強度補正方法の第２の実施例のフローチャートが示されている。図８にはＣＲイメージ内の露光不足領域の強調に用いられる非線形変換を表すグラフが示されている。図９のａには露光不足領域の強度補正を行う前の強度値のプロフィルが水平走査線に沿って示されている。図９のｂには非線形強調および線形回帰変換によって補正された後の図９のａと同じ露光不足領域の強度値のプロフィルが水平走査線に沿って示されている。
【００１４】
有利な実施例の説明
本発明はＣＲイメージング法、特に非線形強調技術および／または線形変換技術を用いてＣＲイメージ内の露光不足領域の強度を補正する方法および装置に関する。
【００１５】
本発明に関連して以下に説明する方法および装置は、ハードウェア、ソフトウェア、ファームウェア、専用または汎用のプロセッサ、またはこれらの組み合わせなど、種々の形態でインプリメンテーションできることを理解されたい。有利には本発明は有形のプログラム命令を含むアプリケーションとしてソフトウェアにインプリメンテーションされる。これは１つ以上のプログラム記憶装置、例えば磁気フロッピーディスク、ＲＡＭ、ＣＤＲＯＭ、ＲＯＭおよびフラッシュメモリなどに格納され、適切なアーキテクチャを含むデバイスまたはマシンによって実行可能となる。
【００１６】
また以下に説明する方法および装置は有利には図ではソフトウェアとしてインプリメンテーションされているが、各モジュールの実際の組み合わせかたやプロセスステップの流れなどは、本発明をプログラミングする際の様態に依存して修正可能であることを理解されたい。本明細書を読む当該の分野の技術者は、本発明の範囲内の考察によって種々のインプリメンテーションまたはコンフィグレーションを行うことができるはずである。
【００１７】
図３には本発明のＣＲイメージ内の露光不足領域のイメージ強度を補正する装置の第１の実施例のブロック図が示されている。装置２０はＣＲイメージ３０のセットを入力する入力手段と、ＣＲイメージ３０内の露光不足領域を検出する検出器３１とを有している。ホワイトバンドが検出された場合、それぞれのＣＲイメージごとに検出器３１はイメージを露光不足部分３２と正常露光部分３３とに分離する。
【００１８】
ＣＲイメージ３０の露光不足部分３２内のオブジェクトは正常露光部分３３の他の解剖学的構造から区別しにくい。つまり露光不足部分３２は正常露光部分３３とは異なる強度レンジを有している。したがって露光不足部分３２に単純な線形スケーリングおよび／または輝度シフトを行うのみでは良好な視認性を得るのに充分でないことがある。本発明では露光不足部分３２を非線形強調モジュール３４によって処理する。このモジュールは露光不足部分３２の強度変化のダイナミックレンジを拡張し、この領域内のオブジェクトイメージを強調する。
【００１９】
正常露光部分３３と非線形強調モジュール３４から出力された強調部分とが線形回帰変換モジュール３５で処理される。ここでは線形回帰変換が行われ、線形の強度変換パラメータまたは強度補正パラメータが求められる。
【００２０】
線形回帰変換モジュール３５は最小化関数を用いて強度補正パラメータを求め、これによりホワイトバンド３２のイメージ強度が正常露光部分３３のイメージ強度にできるだけ近似するように調整する。線形回帰変換モジュール３５で求められた線形変換のパラメータは強度補正モジュール３６で用いられ、ここで入力されたＣＲイメージ３０の露光不足部分３２のイメージ強度を補正するプロシージャが実行される。このプロセスは各ＣＲイメージについて行われ、これらがモザイクイメージとして結合される。強度補正された各ＣＲイメージはモザイクイメージモジュール３７において相互に結合される。
【００２１】
ＣＲイメージ中の露光不足部分３２またはホワイトバンドを検出するために種々の適切なプロセスまたは技術を使用可能であることを理解されたい。有利な例として、２００１年５月１５日付提出の米国出願第０９８５５９５６号明細書“Method and System for Automatic Computed Radiography (CR) Image Composition by White Band Detection and Consistency Rechecking”に記載されているホワイトバンド検出器が挙げられる。この明細書の内容は本発明に関連している。
【００２２】
図４には、本発明のＣＲイメージ内の露光不足領域のイメージ強度を補正する装置の第２の実施例のブロック図が示されている。装置２０’は図３の装置２０と類似しているが、イメージ補正のためにホワイトバンド部分３２に非線形強調を適用するか否かを自動的に判別する手段が設けられている点が異なる。実際には露光不足部分３２を強調すると、オブジェクトイメージが過度に強調されてしまうことがある。この場合には線形の強度変換プロセスのみ用いて露光不足部分３２の強度を正常露光部分３３の強度とほぼ等しくする（またはこれにできるだけ近似させる）調整を行えば充分である。このため強調のプロシージャがバイパスされる。換言すれば、図４の実施例では、非線形強調が不要であると判別された場合、強度補正パラメータは入力されたＣＲイメージ３０の露光不足部分３２と正常露光部分３３とのあいだを直接に線形回帰変換することによって求められ、これがホワイトバンド３２のイメージ強度の調整に用いられる。
【００２３】
より詳細には、図４のシステム２０’では、露光不足部分３２および正常露光部分３３またはこれらの一部が直接に線形回帰変換モジュール３５で処理され、一方の線形変換パラメータセットが求められる。さらに正常露光部分３３と非線形強調モジュール３４から出力された強調部分とが線形回帰変換モジュール３５ａで処理され、他方の線形回帰変換パラメータセットが求められる。ここで選択モジュール３８が線形回帰変換モジュール３５，３５ａで求められた２つの強度補正パラメータセットのうち最小残留誤差を有するほうを判別して選択する。選択された強度補正パラメータが強度補正モジュール３６で使用され、露光不足部分３２または強調部分に対して関連するプロシージャにより強度補正が行われる。全てのＣＲイメージについて強度補正が行われたのち、各イメージはモザイクイメージモジュール３７において相互に結合される。
【００２４】
線形回帰変換モジュール３５，３５ａは、図４に示されているように、同じ機能を実行するものであると理解されたい。図４の実施例では並列処理のために複数のモジュール３５，３５ａが用いられている。装置２０’は上述の処理を並列にではなくシーケンシャルに行う１つの線形回帰変換モジュールを有するものであってもよい。
【００２５】
図５には本発明の強度補正を行うモザイクＣＲイメージング法の実施例のフローチャートが示されている。図５に示されているのは図３，図４の装置の１つの動作モードである。図５では、まず最初にステップ５０で相応のＣＲイメージのセットが取得され、装置に入力される。次にステップ５１で処理すべき１つのＣＲイメージが選択される。続いて当該のＣＲイメージからステップ５２で露光不足領域の検出が行われる。露光不足領域が存在する場合、すなわちステップ５２の結果がＹＥＳであった場合には、ステップ５３で、イメージは２つの部分、すなわち相応のホワイトバンド部分とそれ以外の正常露光部分とに分離される。次にステップ５４では、露光不足部分のイメージ強度が正常露光部分のイメージ強度にほぼ等しくなるように、またはこれにできるだけ近似するように調整される。ステップ５４では図６，図７に則して説明するようなプロセスを用いることができる。ステップ５５で強度補正値が各ＣＲイメージの全ての露光不足部分に適用されたことが確認されると、ステップ５６で適切な周知の手法により各ＣＲイメージは相互に結合され、モザイクイメージが形成される。
【００２６】
図６には本発明にしたがってＣＲイメージ内のホワイトバンド部分のイメージ強度を補正する手法の第１の実施例のフローチャートが示されている。より詳細には、図６には、非線形および線形の変換プロセスを用いた強度補正法が示されている。これらの手法は図５の方法のステップ５４で用いられる。
【００２７】
図６によれば非線形のイメージ強度変換法がＣＲイメージのホワイトバンド部分に適用される。このときオブジェクトイメージ（例えば隠れた骨部）を強調する特性を有する種々の非線形変換を用いることができることを理解されたい。有利な実施例では、次のような非線形強調プロセスが用いられる。
【００２８】
【数１】
2005505960.tif 000003
【００２９】
本発明では変換プロセスは式（１）で示されるものに限定されないことを理解されたい。ＣＲイメージ内のオブジェクト（例えば骨部など）を強調する性質を有するものであれば式（１）のプロセスに代えて他の非線形変換をステップ６０で使用することもできる。
【００３０】
非線形変換により露光不足部分の強度は正常露光部分のレンジとは異なるレンジへ拡張されるので、次に線形変換が使用され、強調部分の強度レンジが調整される。第１のステップ６１では全強調部分と同じ大きさの正常露光部分のサンプルとが取り出される。続いてこれらの部分が線形回帰変換により処理され、ステップ６２で強度補正パラメータが求められる。ここで正常露光部分がホワイトバンド部分よりも小さいときには全正常露光部分と同じ大きさの露光不足部分が取り出され、有利には線形回帰変換によって処理されることを理解されたい。
【００３１】
また、露光不足部分（ホワイトバンド部分）の強度レンジを調整して正常露光部分の強度にほぼ等しくさせるかできるだけ近似させることのできる性質を有するものであれば、種々の線形変換を使用できることを理解されたい。本発明の有利な実施例では、次のような線形回帰変換プロセスがステップ６２で行われる。
【００３２】
Ｉ１（ｘ，ｙ）が強調部分すなわち非線形変換後のホワイトバンド部分の強度関数であると仮定する。強度の調整は、有利には、ステップ６３で
Ｉ１’（ｘ，ｙ）＝ａＩ１（ｘ，ｙ）＋ｂ                          （３）
として求められる。ここでＩ１’（ｘ，ｙ）は調整された強調部分の強度関数であり、ａ，ｂはステップ６２の線形回帰変換プロシージャで求められた線形変換パラメータである。
【００３３】
さらにＩ２（ｘ，ｙ）が同じＣＲイメージ内の強調部分Ｉ１（ｘ，ｙ）と同じ大きさの正常露光部分すなわちサブイメージの強度関数であると仮定する。サブイメージＩ２（ｘ，ｙ）は有利にはホワイトバンド部分Ｉ１（ｘ，ｙ）に隣接する箇所で選択される。線形変換パラメータａ，ｂは、有利には
【００３４】
【数２】
2005505960.tif 000004
の最小化にしたがって求められる。この最小化の物理的な意味は、Ｉ１（ｘ，ｙ）の明るさがＩ２（ｘ，ｙ）の明るさにできるだけ近似するように調整されるということである。ここからａ，ｂについての２つの線形の式
【００３５】
【数３】
2005505960.tif 000005
が得られる。ここでＫはＩ１（ｘ，ｙ）のピクセルの全数である。この式がパラメータａ，ｂについて解かれる。
【００３６】
強度補正パラメータａ，ｂが求められると、ステップ６３で強調部分は例えば式（３）のパラメータを用いて調整される。
【００３７】
図９には図６の手法を用いて露光不足部分の強度を補正したときの強度のグラフが示されている。より詳細には、図９のａには強度補正を行う前の露光不足部分の強度プロフィルが露光不足部分の水平走査線に沿って示されている。図９のｂには図９のａと同じ領域に非線形強調および線形回帰変換を行ったのちの強度プロフィルが露光不足部分の水平走査線に沿って示されている。
【００３８】
図７には本発明にしたがってＣＲイメージ内のホワイトバンド部分のイメージ強度を補正する手法の第２の実施例のフローチャートが示されている。図７の方法は図５のステップ５４で用いられ、ここでは露光不足部分に非線形強調を適用すべきか否かが自動的に判別される。図７に示されているのは図４の装置２０’の動作モードである。
【００３９】
より詳細には、図７の方法では、ステップ７０〜７２のプロセスとステップ７３〜７４のプロセスとにより強度補正パラメータの２つのセットが形成される。ステップ７０〜７２のプロセスは図６の方法のステップ６０〜６２と同様である。このステップ７０〜７２のプロセスでは、非線形強調された強調部分および正常露光部分に対して線形回帰変換を行うことによって第１のパラメータセットが形成される。さらにステップ７３〜７４のプロセスでは、オリジナルの露光不足部分および正常露光部分に対して線形回帰変換を行うことによって第２のパラメータセットが形成される。これらのプロセス（ステップ７０〜７２のプロセスおよびステップ７３〜７４のプロセス）は並列に行ってもよいしシーケンシャルに行ってもよいことを理解されたい。
【００４０】
強度補正パラメータの各セットはステップ７５で評価され、どちらのセットが最小残留誤差を生じるかが求められる。そののち最小残留誤差を生じるほうのセットがステップ７６で選択され、ＣＲイメージ内のホワイトバンド部分のイメージ強度補正が行われる。有利にはこの強度補正は前述の式（３）を用いて行われ、ここでＩ１（ｘ，ｙ）は強調部分またはオリジナルの露光不足部分のいずれかを表す。
【００４１】
このように、本発明の強度補正技術によれば、各ＣＲイメージについて例えば図６，図７の方法を用いて露光不足部分の強度が補正される。強度補正された各ＣＲイメージは周知のモザイク技術を用いて相互に結合される。このようにして有利にはモザイクイメージからアーティファクトが低減されるか除去される。
【００４２】
本発明を添付図の実施例に則して説明したが、本発明は前述した方法および装置のみに限定されるものではなく、当該の分野の技術者であれば本発明の範囲内で種々に修正および変更を加えることができることを理解されたい。こうした全ての修正および変更も特許請求の範囲によって定義される本発明の対象となりうる。
【図面の簡単な説明】
【００４３】
【図１】従来のＣＲイメージング方法を示す図である。
【図２】発生するアーティファクトの例を示す図である。
【図３】本発明の装置の第１の実施例を示す図である。
【図４】本発明の装置の第２の実施例を示す図である。
【図５】強度補正を行うモザイクＣＲイメージング法を示すフローチャートである。
【図６】強度補正方法の第１の実施例を示すフローチャートである。
【図７】強度補正方法の第２の実施例を示すフローチャートである。
【図８】非線形変換を表すグラフである。
【図９】強度補正を行う前と行った後の強度値のプロフィルを示すグラフである。
【図１】
2005505960.tif 000006
【図２】
2005505960.tif 000007
【図３】
2005505960.tif 000008
【図４】
2005505960.tif 000009
【図５】
2005505960.tif 000010
【図６】
2005505960.tif 000011
【図７】
2005505960.tif 000012
【図８】
2005505960.tif 000013
【国際公開パンフレット】
2005505960.tif 000014
2005505960.tif 000015
2005505960.tif 000016
2005505960.tif 000017
2005505960.tif 000018
2005505960.tif 000019
2005505960.tif 000020
2005505960.tif 000021
2005505960.tif 000022
2005505960.tif 000023
2005505960.tif 000024
2005505960.tif 000025
2005505960.tif 000026
2005505960.tif 000027
2005505960.tif 000028
2005505960.tif 000029
2005505960.tif 000030
2005505960.tif 000031
2005505960.tif 000032
2005505960.tif 000033
2005505960.tif 000034
2005505960.tif 000035
2005505960.tif 000036
2005505960.tif 000037
2005505960.tif 000038
2005505960.tif 000039
2005505960.tif 000040
2005505960.tif 000041
2005505960.tif 000042
2005505960.tif 000043
2005505960.tif 000044
2005505960.tif 000045
2005505960.tif 000046
【国際調査報告】
2005505960.tif 000047
2005505960.tif 000048
2005505960.tif 000049