JP2004160208A - 結合フル画像再現用の部分放射線画像の自動配列決定方法 - Google Patents
結合フル画像再現用の部分放射線画像の自動配列決定方法 Download PDFInfo
- Publication number
- JP2004160208A JP2004160208A JP2003348652A JP2003348652A JP2004160208A JP 2004160208 A JP2004160208 A JP 2004160208A JP 2003348652 A JP2003348652 A JP 2003348652A JP 2003348652 A JP2003348652 A JP 2003348652A JP 2004160208 A JP2004160208 A JP 2004160208A
- Authority
- JP
- Japan
- Prior art keywords
- sub
- images
- image
- virtual
- screen
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Withdrawn
Links
- 238000000034 method Methods 0.000 title claims abstract description 45
- 230000005855 radiation Effects 0.000 title abstract description 4
- 238000005314 correlation function Methods 0.000 claims abstract description 20
- 238000005259 measurement Methods 0.000 claims abstract description 12
- 230000035945 sensitivity Effects 0.000 claims abstract description 9
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 claims description 66
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 238000012163 sequencing technique Methods 0.000 claims description 2
- 238000001514 detection method Methods 0.000 description 14
- 230000008569 process Effects 0.000 description 14
- 230000007704 transition Effects 0.000 description 14
- 230000000875 corresponding effect Effects 0.000 description 11
- 238000010586 diagram Methods 0.000 description 7
- 238000003384 imaging method Methods 0.000 description 7
- 238000006073 displacement reaction Methods 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 4
- 238000003491 array Methods 0.000 description 4
- 239000002131 composite material Substances 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000013459 approach Methods 0.000 description 3
- 230000002596 correlated effect Effects 0.000 description 3
- 238000002601 radiography Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 238000012886 linear function Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 230000000737 periodic effect Effects 0.000 description 2
- 210000001519 tissue Anatomy 0.000 description 2
- 210000000988 bone and bone Anatomy 0.000 description 1
- 239000006227 byproduct Substances 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 239000011888 foil Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000001454 recorded image Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 206010039722 scoliosis Diseases 0.000 description 1
Images
Classifications
-
- G06T3/147—
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
- G06T7/32—Determination of transform parameters for the alignment of images, i.e. image registration using correlation-based methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10116—X-ray image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30008—Bone
- G06T2207/30012—Spine; Backbone
Abstract
【課題】結合された完全な画像を自動的に構成する。
【解決手段】長い対象物の少なくとも2つの放射サブ画像を、重なり合う記録メディアにより取得してデジタル画像へ変換する。サブ画像の順番、向き、重なりの考えられる全ての配列を含む仮想リストを構築し、各仮想におけるサブ画像の対に対して重なり領域を検出し、対毎に複数の異なる計測を実行し、計測結果が範囲外の場合に前記リストから当該仮想を消去する。残存する各仮想における重なり領域の相関関数を算出し、関数最大値及び対のサブ画像間の水平距離を求め、連続したサブ画像の各対の水平距離の大きさを検出し、該大きさが範囲外の場合に前記仮想リストから当該仮想を消去する。相関関数の最大値の合算値に基づき前記仮想リストに残存する各仮想に対して全体の最小感度を確立し、最大の最小感度の仮想を最も良好な候補に選択し、選択された仮想に基づいてサブ画像を結合処理する。
【選択図】図5
【解決手段】長い対象物の少なくとも2つの放射サブ画像を、重なり合う記録メディアにより取得してデジタル画像へ変換する。サブ画像の順番、向き、重なりの考えられる全ての配列を含む仮想リストを構築し、各仮想におけるサブ画像の対に対して重なり領域を検出し、対毎に複数の異なる計測を実行し、計測結果が範囲外の場合に前記リストから当該仮想を消去する。残存する各仮想における重なり領域の相関関数を算出し、関数最大値及び対のサブ画像間の水平距離を求め、連続したサブ画像の各対の水平距離の大きさを検出し、該大きさが範囲外の場合に前記仮想リストから当該仮想を消去する。相関関数の最大値の合算値に基づき前記仮想リストに残存する各仮想に対して全体の最小感度を確立し、最大の最小感度の仮想を最も良好な候補に選択し、選択された仮想に基づいてサブ画像を結合処理する。
【選択図】図5
Description
本発明は、デジタル放射線透過写真術に関し、特に、記憶蛍光体ベースのコンピュータX線撮影法システムを用いて、脊椎や足のような、人体の長い部分の撮像に関する。
人体の長い部分が従来のスクリーンフィルム技術を用いて撮像されたとき、14インチ×36インチ及び14インチ×51インチのような、拡大された長さの特別なカセット及びフィルムが使用される。医療機関は、アナログスクリーンフィルムからコンピュータX線写真(CR)のようなデジタルモダリティに乗り換えているので、これらの種の試験は、多大な挑戦となる。これは、デジタル検出器の大きさが制限されるからである。例えば、幾つかの主要なCR供給業者からの最も大きな記憶蛍光体カセットは、14インチ×17インチに制限され、これは、一回に長い体の部分の一部しか撮像できない。
幾つかの方法が、既存の標準CRカセットを幾つかスタックすることによってCR撮像範囲を拡大すべく提案されている。例えば、特許文献1,2は、幾つかの記憶蛍光体カセットを近傍に重ねる方法を開示する。カセットは、交互(図1A)、階段状(図1B)若しくは傾斜状(図1C)の配列でありうる。X線を当てている間、全ての部分的に重なるカセットは、同時に晒され、それ故に、対応するカセット内にある各記憶蛍光体スクリーン102が、長い体部の画像の一部を記録する。特許文献3,4には、同様のアプローチとして、撮像範囲を拡大するための、複数の重なる記憶蛍光体スクリーンを基本とした方法が開示されている。スクリーン102は、同様に、単一の細長いカセット103内で、交互(図1D)、階段状(図1E)若しくは傾斜状(図1F)で重なる配列で構成されうる。同様のアプローチが特許文献5,6にも提案されている。更に、カセット及びスクリーンは、図1Gに示すように、交互配列で共に使用されてよく、カセット101及び記憶蛍光体スクリーン102は、スクリーン102が常にカセット101の前側に位置する状態で、部分的に重なる交互配列で配置される。この方法は、取得した画像からカセットの陰を無くし、カセットから取り出し元に戻す必要がある記憶蛍光体スクリーンの数を低減する。
蛍光体スクリーンは、基本的な撮像記録装置であるので、スクリーンが個々のカセット内にパッケージされているか否かに関係なく、用語“記憶蛍光体スクリーン”、“蛍光体スクリーン”若しくは“スクリーン”が、蛍光体スクリーン自体若しくはカセット内に移送された蛍光体スクリーンを表わすために、以下で使用される。それ故に、図1における異なるシナリオは、スクリーン面間の距離が、スクリーンがカセット内にあるかないかに依存して各シナリオに対して異なることを除いて、重なる蛍光体スクリーンの交互、階段状及び傾斜状の配列と考えられる。
患者のX線写真画像の取得方法の概略が図2に示される。患者(要素203)は、X線源(要素201)と複数のスクリーン(要素205)との間に配置される。図1Aないし図1Gに示すスクリーン配列方法の何れが撮像に使用されてよい。選択的な散乱防止用グリッドが、患者とスクリーンとの間に配置されてよい。グリッドは、固定タイプ若しくは往復動タイプであってよい。X線露光中、X線は平行化されてよく、患者の組織への診断に無関係な放射を最小化する。X線生成器が点火され、カセットが露光された後、患者の画像は、潜在的なX線写真信号として複数のスクリーンにより記録される。各スクリーンは、患者の画像の一部のみ捕捉する。スクリーンは、CRリーダーに供給され、潜在X線写真信号が、電子画像に変換される。個々のスクリーンから取得される電子画像を、サブ画像と称する。
個々の記憶蛍光体スクリーンにより取得されるサブ画像は、複合の完全な画像を生成するために結合されなければならない。X線露光中に使用された蛍光体スクリーンの空間的な順番、向き、重なり配列に関する情報が、複合画像を結合するために必要とされる。X線露光後、蛍光体スクリーンは、CRリーダで任意の順番で走査されてよい。従って、走査したサブ画像を、画像取得で使用された物理的設定に対応する順番に再構成することが必要である。また、連続したスクリーン間の重なり配列が正確に知られる必要がある。例えば、スクリーンは、階段状及び傾斜状スクリーン配列の場合上部若しくは下部に重なりうり、交互スクリーン配列の場合、X線源により近いスクリーンが、遠いものとは異なる態様で重なりうる。スクリーンは、2つの長手寸法の何れかを上として、(横向き)水平な向きで画像捕捉のために等価的に配設されてよい。結果として、走査されるサブ画像は、近傍のサブ画像から180°回転されてよい。サブ画像の回転の検出がそれ故に必要である。
特許文献7,8は、それぞれ、サブ画像のセットから複合のX線写真画像を再構成するための方法を開示する。前者の方法は、全てのサブ画像が、所定の空間的な位置及び向きで連続的に取得及び走査されることを前提とする。後者の方法は、ハードウェア位置センサに依存して、各サブ画像の取得中の相対位置を決定する。これらの何れの方法であっても、走査されたサブ画像のシーケンスが、取得中に使用されたシーケンスと一致しないときの状況には適用できない。特許文献9,10は、サブ画像のアライメントのために参照マーカーのパターンを利用する結合方法を提案する。しかし、この特許は、どのようにしてサブ画像を順列化するかについて教示が無く、サブ画像が結合処理前に正確に方向付けられていることを必要とする。階段状配列の蛍光体シートから読み出されたサブ画像を自動的に結合するための、その他の結合方法は、Wei他により提案されている(例えば、非特許文献1参照)。しかし、この方法は、サブ画像が結合前に既に対状に連続配列されている。
特許文献11は、全ての蛍光体スクリーンに、対応するスクリーンの向き及び位置を決定するための補助情報として、付した2つの異なる識別ラベルに依存する方法を開示する。特許文献11は、更に、どのようにして補助情報を結合処理で使用しうるかについて教示する。しかし、補助情報を全く使用しないことが望ましく、というのは、これは、標準のカセットが修正される必要があり、若しくは、標準の蛍光体スクリーンIDが置換される必要があることを意味し、何れの場合も、カセット及びスクリーンを他の一般用途に互換性のないものとするからである。
欧州特許出願0919856号A1
米国特許第6,273,606 B1
欧州特許出願0866342 A1
米国特許第5,986,279号 B1
日本特許出願2000-267210 A
日本特許出願2000-241920 A
米国特許第4,613,983号
米国特許第5,833,607号
欧州特許出願0919858 A1
米国特許第6,269,177号 B1
日本特許出願2000-258861 A
日本特許出願2000-232976 A
「A new fully automatic method forCR image composition by white band detection and consistency rechecking」Guo−Qing Wei他、Proceedings of SPIE Medical Imaging, 2001, vol.4322, pp 1570-1573
従って、取得されたサブ画像に直接的に基づいてX線露光で使用される蛍光体スクリーンの空間的な順番、向き、重なり配列を決定する自動方法を発展させることは望ましい。
本発明によれば、上述の問題点に対する解決策が提供される。
本発明の特徴によれば、結合された完全な画像を構成するための、部分的な放射線画像の自動配列決定方法であって、
長い対象物の少なくとも2つの放射サブ画像を、重なり合う記録メディアにより、取得する取得ステップと、
前記記録されたサブ画像をデジタル画像の変換し、保存する変換ステップと、
サブ画像の順番、向き及び重なりの考えられる全ての固有の配列を含む仮想リストを構築するステップと、
各仮想における各連続したサブ画像の対に対して重なり領域を検出するステップと、
各仮想の連続したサブ画像の対毎に複数の異なる計測を実行し、計測結果が範囲外の場合に前記リストから当該仮想を消去するステップと、
前記仮想リストに残存する各仮想における各連続したサブ画像に対して重なり領域の相関関数を算出し、関数最大値及び対のサブ画像間の水平距離を求めるステップと、
連続したサブ画像の各対の水平距離の大きさをチェックし、該大きさが範囲外の場合に前記仮想リストから当該仮想を消去するステップと、
相関関数の最大値の合算値に基づいて前記仮想リストに残存する各仮想に対して全体の最小感度を確立するステップと、
最も良好な候補に対して最大の最小感度の仮想を選択するステップと、
結合されたサブ画像からなる出力画像を生成するため、選択された仮想に基づいてサブ画像を結合処理するステップとを含む、方法が提供される。
長い対象物の少なくとも2つの放射サブ画像を、重なり合う記録メディアにより、取得する取得ステップと、
前記記録されたサブ画像をデジタル画像の変換し、保存する変換ステップと、
サブ画像の順番、向き及び重なりの考えられる全ての固有の配列を含む仮想リストを構築するステップと、
各仮想における各連続したサブ画像の対に対して重なり領域を検出するステップと、
各仮想の連続したサブ画像の対毎に複数の異なる計測を実行し、計測結果が範囲外の場合に前記リストから当該仮想を消去するステップと、
前記仮想リストに残存する各仮想における各連続したサブ画像に対して重なり領域の相関関数を算出し、関数最大値及び対のサブ画像間の水平距離を求めるステップと、
連続したサブ画像の各対の水平距離の大きさをチェックし、該大きさが範囲外の場合に前記仮想リストから当該仮想を消去するステップと、
相関関数の最大値の合算値に基づいて前記仮想リストに残存する各仮想に対して全体の最小感度を確立するステップと、
最も良好な候補に対して最大の最小感度の仮想を選択するステップと、
結合されたサブ画像からなる出力画像を生成するため、選択された仮想に基づいてサブ画像を結合処理するステップとを含む、方法が提供される。
本発明は、次の効果を有する。1.本発明は、取得したサブ画像に直接的に基づいて、X線露光で使用される蛍光体スクリーンの空間的順番、向き、重なり配列の自動決定を可能とする。2.本方法は、蛍光体スクリーンが任意に重なり、180°の多様性で配置され、あらゆる順番でCRレコーダーで走査されることを許容する。3.本発明が自動画像結合方法と組み合わせられる時、複合画像をサブ画像のセットから生成する全プロセスを完全に自動化できる。
一般に、本発明は、例えば脊柱側弯症の診断のための完全な背骨若しくは人体の完全な足のような細長い物体のX線放射撮像に関する。
図3Aは、3つの記憶蛍光体スクリーンが露光される場合に対する一例を示す。画像は、図2の放射露光方法及び図1Gの蛍光体スクリーン構成を用いて取得される。第1の画像(要素301)及び第3の画像(要素303)を捕捉する2つの蛍光体スクリーンは、第2の画像(要素302)を捕捉するスクリーンの背後に配置される。第1のスクリーン301及び第2のスクリーン302は、部分的に重なり、第2のスクリーン302は、入射X線に全体的に傾斜していないので、第1のスクリーン301は、依然としてスクリーン重なり領域(307)の患者の画像を捕捉する。しかし、重なり領域で第1のスクリーン301上に捕捉される画像の信号対雑音比は、第2のスクリーン302によるX線減衰により比較的低い。第2のスクリーン302の上縁は、要素305によって示すように、第1の画像内に明確な影を付与する。同様に、第2のスクリーン302の下縁は、スクリーン303に要素306により示すように、第3の画像内に影を付与する。全ての3つの画像において、平行化露光領域と非平行化露光領域との間の境界は、要素304により示される。
上述の如く、露光された記憶蛍光体スクリーンは、CRリーダーで読み出され、デジタル化される電子画像を生成し、画像は、ピクセルのマトリックスからなる。本発明によれば、記憶蛍光体スクリーンは、記憶蛍光体スクリーンの露光画像領域の上縁及び下縁を超えてオーバースキャンされる。
如何なる2つの近傍のスクリーン間の関係は、図1Aないし図1Gで示すスクリーン構成に対して等価である。スクリーンは、部分的に重なり合い、一のスクリーンがX線源により近い位置に配置される。用語“フロントスクリーン”が、X線源に近くにあるスクリーンを指すために使用され、用語“フロント画像”は、フロントスクリーンで捕捉された画像を指すために使用される。同様に、用語“バックスクリーン”は、X線源から遠い側のスクリーンを指すために使用され、用語“バック画像”は、“バックスクリーン”で捕捉された画像を指すために使用される。N個の近傍の重なるスクリーンの集合が、X線露光のための使用されるとき(N≧2)、スクリーンは、N−1個のスクリーン対へと分解でき、それぞれの対は、2つの隣接するスクリーン、即ち1つのフロントスクリーンと1つのバックスクリーンからなる。この定義により、一のスクリーンは、一の“スクリーン対”のフロントスクリーンであると共に、次のスクリーン対のバックスクリーン若しくはフロントスクリーンでありえる。
図4A及び図4Bは、本発明を説明する際に有用な追加の用語を示し、スクリーン303,402のバック重なり縁部及びフロント重なり縁部が、要素406,407によりそれぞれ指示されている。装置401は、X線源である。
本発明の次の説明において、用語“鉛直”は、蛍光体スクリーンがスタックされる方向を意味するために使用され、用語“水平”は、“鉛直”に垂直な方向を意味するために使用される。スクリーンがX線露光中にあらゆる方向に向けられてもよいことは理解されるべきである。
図5を参照するに、取得されたサブ画像からの直接的な、X線露光の蛍光体スクリーンの空間的な順番、向き及び重なり配列の自動決定のための本発明の一実施例が示される。本方法は、細長い対象のサブ画像を取得し(501);、各サブ画像を読み出し及び記憶し(502);、サブ画像の順番、向き及び重なりの考えられる固有の配列を全て含む仮想リストを構築し(503);、各仮想における連続したサブ画像の各対に対して領域を検出し(504);、各仮想の連続したサブ画像の対毎に複数の計測を実行し、例えば、重なり領域のサイズ、重なり領域におけるフロント画像及びバック画像間での平均ピクセル強度差、バック画像におけるフロントスクリーンの重なりエッジの影の全体の遷移大きさ、若しくは、フロント画像における支配的な水平ラインの全体の遷移大きさのような、計測結果が範囲外の場合、前記リストから当該仮想を消去し(505);、前記仮想リストに残存する各仮想における各連続したサブ画像に対して重なり領域の相関関数を算出し、関数最大値及び対のサブ画像間の水平距離を求め(506);、各連続したサブ画像の対の間の水平距離の大きさをチェックし、該大きさが範囲外の場合に前記仮想リストから当該仮想を消去し(507);、相関関数(複数でも可)の最大値の合算値に基づいて前記仮想リストに残存する各仮想に対して全体の最小感度を確立し(508);、最大の最小感度の仮想を最も良好な候補に選択し(509);、結合されたサブ画像からなる出力画像を生成するため、選択された仮想に基づいてサブ画像を結合処理する(510)を含む。505及び507のステップは、連続するサブ画像の対の局所レベルに基づいて動作し、これらの誤り仮想を拒絶するために大まかに定義された閾値で測定値を用いる。508のステップは、仮想の全体としてのレベルに基づいて動作し、509のステップは、最小感度に基づいて最も可能性のある候補を特定するため、全体のレベルに基づいて動作する。アルゴリズムを異なる解析レベルで動作させることによって、コンピュータ効率を大幅に向上しつつ、検出のロバストさに対する要求を満足することができる。
この本発明の詳細な実施例が、図1Gのスクリーン配列に基づいて説明される。この作業の精神は、交互、階段状、傾斜状若しくは他の任意の重なり配列に適用可能であることは理解されるべきである。
本発明の方法は、スクリーン画像をデジタルサブ画像に変換するCRリーダー、モニタ上で表示できる若しくは横長としてプリント出力できる出力画像を生成するために前記記憶から処理まで実行する、デジタルコンピュータやマイクロプロセッサのような画像プロセッサにより実行される。
本発明の初期ステップは、サブ画像の順番、向き及び重なりの考えられる固有の配列を全て含む仮想リストを構築することである。例えば、図8は、行A乃至H、列1−6で、交互に重なるスクリーンの3つのスクリーンに対するサブ画像の考えられる全ての配列を示す。各仮想図では、正方形1,2,3は、個々のサブ画像を表わし、中の数字は、蛍光体スクリーンがCRリーダーにより走査されたときに固有に割り当てられる。白抜きの正方形は、対応するサブ画像が、X線源に近い側の蛍光体スクリーンから取得されたことを意味し、影のついた正方形は、X線源に遠い側の蛍光体スクリーンから取得されたことを意味する。サブ画像が回転されたとき、対応する数は、記号“’”が付される。例えば、“2’”は、サブ画像2が180°回転されたことを意味する。
より具体的には、行A,B,C,Dは、上及び下のサブ画像がX線源に近く、中央のサブ画像がX線源から遠い、種々のサブ画像配列を示している。行E,F,G,Hは、上及び下のサブ画像がX線源から遠く、中央のサブ画像がX線源に近い、種々のサブ画像配列を示している。A行1列、乃至、A行6列は、それぞれ、上から下までの画像配列を示す。即ち、1,2,3;1,3,2;2,1,3;2,3,1;3,1,2;及び3,2,1である。行B及びCでは、サブ画像の3つの配列が、各ボックス内に示されている。B行1列は、サブ画像配列11,2,3;1,21,3;1,2,31を示し、B行2列は、サブ画像配列11,3,2;1,31,2;1,3,21を示す。同様に、C行1列は、サブ画像配列11,21,3;11,2,31;1,21,31を示し、C行2列は、サブ画像配列11,31,2;11,3,21;1,31,21を示す。更に進み、単一のサブ画像配列が、行D,E,Hのボックス内に示され、3つのサブ画像配列が、行F,Hのボックス内に示されている。行B列3,4,5;行C列3,4,5;行D列3,4,5;行F列3,4,5;行G列3,4,5;行H列3,4,5のボックスに対しては、示されていないが、示したサブ画像から容易に導出可能である。
次のステップは、各仮想における連続したサブ画像の各対に対して重なり領域を検出することである(図5の504参照)。これは、クリティカルなステップである、というのは、多くの解析が、重なり領域内で実行されるからである。バック画像における重なり領域は、比較的低い信号大きさを有し、フロントスクリーン重なり縁の影により画像の残りから簡易に分離される。フロントスクリーン重なり縁が、フロント画像において知られている場合、それをバック画像における影と対応させることができ、フロント画像における重なり領域は、バック画像の重なり領域のサイズに基づいて決定できる。フロント画像のフロントスクリーン重なり縁を確定論的に位置特定する一方法は、CRリーダーに重なり縁を越えてフロント蛍光体スクリーンをオーバースキャンさせることであり、これは、通常的には利用可能でなく若しくは従来のCRリーダーで保証されていない。例えば図3Bに示す修正画像(要素310)を使用した場合、蛍光体スクリーンの上下の縁部が、オーバースキャンされ、それ故に、両縁部が、サブ画像ピクセルマトリックス中で完全に見ることができる(要素311,312)。
特定のサブ画像の対の重なる領域を検出するための次の説明は、CRリーダーが全ての蛍光体スクリーンの上下の縁部の双方をオーバースキャンできることを想定する。多くの変形が本記載の精神に基づいて考え出されることは理解されるべきである。
スクリーン重なり縁部は、先ず、フロント及びバック画像から同様に位置特定され、次いで、バック画像内のフロントスクリーン重なり縁部の影が検出される。これらのステップは、図6に示されている。スクリーン重なり縁部を越える画像領域のピクセル値は、CRリーダーの基線ノイズレベルを反映する。これは、蛍光体スクリーンから信号寄与が無いからである。結果的に、これらの領域でのピクセル値は、露光された画像領域に比して低く、それ故に、画像のスクリーン重なり縁部を横切る急なピクセル値増分/不連続がある。このピクセル値不連続は、スクリーン重なり縁部の方向及び位置を検出するために使用され、これは、多くの方法により実現できる。本発明の好ましい実施例では、検出は、(1)スクリーン重なり縁部近傍の大きなエッジ遷移ピクセルの全てを計算し、(2)スクリーン重なり縁部の候補遷移ピクセルのライン表現を実行することにより、実現される。
例としてフロント画像を使用して、図6は、検出プロセスの好ましい実施例を説明する。第1に、狭いバンド602が、フロント画像600の端部から抽出される。どのようにして蛍光体スクリーンがCRリーダーで走査されているかに依存して、スクリーン端部縁部601の向きは、各走査間で取得画像において幾つかの度合いの変動を有することができる。それ故に、狭いバンドのサイズは、全体のスクリーン端縁が高い信頼性で抽出できるほど十分に大きくなければならない。2048ピクセルの幅を有する画像に対しては、狭いバンドのサイズは、約200×2048ピクセルであるべきである。
第2に、画像の一次導関数が、鉛直方向で計算される。[−1,0,1]、[−1,0,0,0,1]若しくは[−1,0,0,0,0,0]等のような一次導関数演算子が好ましい、というのは、ピクセル値不連続性は、常に略水平であるエッジ方向を横切って発生するだけであるからであり、また、計算効率的な効果による。好ましい閾値が、より大きな大きさであり且つ下降する傾きである候補縁部遷移ピクセルのみを選択するために使用される。要素603は、このステップからの結果を示す。狭いバンドに2次元のローパスフィルタ処理を施すことは、検出精度を向上することを助けることがわかる。この本発明の好ましい実施例では、狭いバンドはローパスフィルタ処理される。
第3に、線形関数が、候補縁部ピクセルに適合され、最も適合するパラメータが、最小自乗誤りに到達したときに得られる。要素604は、縁部遷移ピクセルの上部に重畳された適合線形関数を示す。このプロセスは、バック画像の開始時に狭いバンド614内に上昇縁部遷移ピクセルが代りに探索される以外、同様に、バック画像610に対して実行される。
スクリーン重なり縁部がフロント画像でうまく見出されると、それは、画像位置合わせのため、バック画像のその影と比較される。バック画像におけるフロントスクリーン重なり縁部の影を位置特定するため、スクリーン重なり縁部を検出する際に使用された方法と類似するアプローチが使用される。これは、バック画像のピクセル値が、フロントスクリーンによる入射X線の高減衰に起因してスクリーン重なり領域内で強い信号強度減少を受けることから、可能である。フロント画像重なり縁部の影を検出するための狭いバンドの位置は、プライオリであるべきである、蛍光体スクリーン間の適切な重なりサイズに基づいて算出できる。狭いバンドのサイズは、探索範囲を広げるために拡大されてもよい。
本発明の全体のプロセスの最初に、スクリーン重なり縁部及びフロントスクリーンの候補影を検出することは、計算効率がよい。各サブ画像は、上部若しくは下部に重なりを備える仮想リスト内のバック画像として処理される機会を有するので、スクリーン重なり縁部及びフロントスクリーンの候補影の検出は、各画像の上部及び下部にて一回限りで実現できる。結果は、保存され後に使用される。用語“フロントスクリーンの候補影”は、連続するサブ画像の対のフロント画像に対して影が実際には存在せず、この場合、特定された線が検出領域内の単なる支配的な水平線であることから、使用される。それ故に各サブ画像に対してラインの2セットが存在し、一は、画像の上部に関連し、他の一は画像の下部に関連し、各セットは、スクリーン重なり縁部ライン及びフロントスクリーン重なり縁部の候補影ラインからなる。これらの4つのラインは、対応して、続く解析で使用される。ライン検出プロセスの副産物として、各ラインを横切る遷移大きさが、計算され、これは、ある範囲内のラインの2つの側での平均ピクセル値の差である。
理論的には、フロント画像におけるフロントスクリーン重なり及びバック画像におけるその影に対する適合関数の検出傾斜は、同一であるべきである。しかし、それらは、実際には、X線露光中の2つの蛍光体スクリーンの間のミスアラインメントや、読み出しプロセス中のCRリーダーでのスクリーン位置変動のような、幾つかの理由により数度だけ異なりうる。偏差は、フロントとバックの画像間の向きミスアラインメントを表わす。このミスアラインメントは、その大きさが全体の検出アルゴリズムの信頼性が影響を受けるようなレベルを超えた場合に、特に連続するサブ画像の対の間の相互関数の算出において、補正されなければならない。
ミスアラインメント補正は、フロント若しくはバック画像を回転させることによって実現される。フロント画像が回転される場合、回転角度は、θ=atan(k)−atan(kf)であり、kf及びkは、それぞれ、フロント画像におけるフロントスクリーン重なり縁部、及び、バック画像におけるその影に対する適合関数の傾斜である。バック画像が回転される場合、角度は、θ=atan(kf)−atan(k)であり、両画像が回転される場合、角度は、フロント及びバック画像に対してそれぞれ−atan(kf)及びatan(k)である。
バック画像におけるスクリーン重なり領域は、バックスクリーン重なり縁部とフロントスクリーン重なり縁部の影の間に位置し、これに基づいて、更なる解析のためにこのバック画像のスクリーン重なり領域が抽出される。フロント画像の重なり領域は、フロントスクリーン重なり縁部から始まり、バック画像の重なり領域の大きさと同一の大きさであり、これに基づいて、このフロント画像の重なり領域が抽出される。
補償ファクターが、CRリーダーがスクリーン縁部を越えて読み出さないときに、重なり領域のサイズに適用できる。例えば、図3Aの要素320は、CRリーダーがリン光縁部から画像ピクセルのmラインを読み損ねたことを示す。重なり領域は、それ故に、mラインだけより小さくなり、これに応じて補償がなされうる。これは、CR走査プロセス中においてスクリーンの向きの多様性が非常に小さい場合にのみ機能する。
次に、複数の測定が、各仮想の連続サブ画像の対毎に実行され、仮想が、測定結果が範囲外の場合にリストから外される。これらの測定は、重なり領域のサイズ、フロント画像及びバック画像内の重なり領域間の平均ピクセル強度差、バック画像におけるフロントスクリーンの重なり縁部の影の全体の遷移大きさ、若しくは、フロント画像における支配的な水平ラインの全体の遷移大きさ等を含む。
重なり領域のサイズは、設計によるカセット若しくはスクリーン保持具により定まる固定の最小及び最大範囲内であるべきである。例えば、100から300ピクセルラインの大まかに定義される範囲が、200ピクセルラインの典型的な重なりサイズに対して使用される。
バック画像の重なり領域は、フロントスクリーンによる減衰に見舞われる。図9は、種々のX線ビームエネルギ(kVp)での平均ピクセル値差を示す。このデータセットは、散乱防止用に各蛍光体の背後にリードホイルを有するコダック社の蛍光体スクリーンに基づく。この図は、フロントスクリーンの重なり領域でのピクセル値が、バックスクリーンのものよりも高いが、その差は、X線ビームエネルギの関数である、ある一定範囲内であるべきである、ということを示す。図9は、最大の差に特に関心が向けられている、というのは、平均入射X線ビームエネルギは、ビーム硬化効果に起因して撮像物体によりより高い端部に向けて押される傾向にあり、これがピクセル値の差を低減するからである。
各連続サブ画像の対間の重なり領域でのピクセル値差は、計算効率を改善するために、全体的及び局所的の双方で検査されてよい。全体のチェックを行うため、重なり領域における全体の平均ピクセル値が、フロント画像及びバック画像に対して別個に計算される。典型的に70−90kVpのX線が、完全な背骨及び完全な足の検査に使用され、対応したピクセル値差は、図9に基づいて590−750周辺のコード値である。最大値1000は、それ故に、平均ピクセル値差に対して選択される。更に、最小値40kVpが、使用され、フロント画像の重なり領域での平均ピクセル値がバック画像よりも高いことが保証される。最小値及び最大値の双方は大まかに定義され、ビームエネルギの多様性を考慮し、また、有効な仮想を早まって排除しないようにする。
重なり領域での平均ピクセル値は、また、検査の特定性を改善するために局所的に算出される。連続サブ画像の各対に対して、フロント画像及びバック画像の重なり領域は、8のような、同一の数のより小さい領域へと仕切られ、各小さな区画内の平均ピクセル値が、計算され、対応して比較される。最大値1000は、全体の平均ピクセル値差の検査に対して同様の理由により選択される。重なり領域において記録画像情報間に幾らかの水平変位がありうるので、正確な変位の大きさを知らないと、局所平均ピクセル値の比較は、異なる画像コンテンツを比較するという危険に陥る。これを克服するために、局所領域のサイズは、最大水平変位よりも十分大きくあるべきであり(例えば、4倍)、これにより、対応する局所区画間の画像コンテンツ不一致が比較的小さくなることが保証される。更に、局所平均ピクセル値差に対して選択される許容範囲は、−100のような、低い側の端部でより広くあるべきである。
仮想リストが、重なり領域の大きさのチェック及び重なり領域での平均ピクセル値差のチェックを受けた後、多数の誤り仮想が、リストから削除され、残りが、続く解析に通されることになる。
バック画像におけるフロントスクリーン重なり縁部の影の全体の遷移大きさのチェックは、リストに残存する各仮想の連続サブ画像の各対に対して実行される。この検査は、フロント蛍光体スクリーンから取得された幾つかのサブ画像がバック画像として扱われる仮想を除去することである。影ラインは、フロントスクリーン重なり縁部により導かれると予測されるので、図9に示すように、最小の遷移大きさは、最小閾値よりも大きくあるべきである。
フロント画像における支配的な水平ラインの全体の遷移大きさのチェックも、リストに残存する各仮想の連続サブ画像の各対に対して実行される。フロント画像はバック画像が重なり領域で受けるのと同一の減衰を受けず、また、フロント画像の画像特徴は、一般的に連続的であるべきであるので、水平ラインの全体の遷移大きさは、ある最大値を超えるべきでない。
尚、全体の検出プロセスは、上述のチェックプロセスの何れかの終了時にリストに唯一の仮想が残存する場合、完了できる。この場合、最後の仮想は、結合のためのサブ画像の最も可能性のある配列を含む。
連続のサブ画像間で相互関数を算出することは、全体の検出プロセスにおいて最も計算的に進歩性のあるステップである。それは、それ故に、最後の2,3の処理ステップ内に配設される。連続サブ画像の各対間の重なり領域で記録される画像コンテンツは、対応するピクセル間の、幾つかの水平変位、y_offsetを除いて、同一である。この特性は、相関関数により最も良く特徴化される、というのは、それは、2つの重なり領域が最も一致する場所にて、即ちy_offset点で、最大となるからである。相関関数は、水平方向でのみ算出される必要がある、というのは、2つの重なり領域は、この方向でしか変位しないからである。関数は、次式を用いて算出される。
本発明の好ましい実施例では、正規化相関関数が、次式を用いて算出される。
本発明の好ましい実施例では、相関関数はフーリエ法を用いて算出される。
図7は、この動作の好ましい実現を示す。第1に、重なり領域702,703が、それぞれ、フロント画像及びバック画像から抽出される。第2に、要素704は、702の一部を抽出することにより得られ、次いで、703と相関が取られ、相関関数c(Δ)、706が生成される。同様の結果は、703の一部を702に対して相関を取ることによって、達成できる。第3に、関数の最大値cmaxが、検出され、対応するΔの値が、y_offset、707として特定される。
スキンライン、組織境界、骨のエッジ、平行化境界(図3Aの要素304)及びハードウェアラベル等を含む702及び703でのエッジ情報は、相関に対して最も有用な情報をもたらすので、本発明の好ましい実施例では、低周波数のコンテンツが、相関のロバスト性を改善するために、702及び703から取り除かれる。平行化の影(要素304)の存在が、y_offsetを見出す際のアルゴリズムロバスト性を改善する補助をすることがわかった。それ故に、X線露光中の平行化を使用することが推奨される。
通常的に相関関数は、円滑であるが、固定のグリッドがX線露光中に使用される場合、それは、取得画像に周期的ラインパターンのアーチファクトを付与する。このアーチファクトは、グリッドが鉛直方向に向くときに特に支配的となり、自身で相関を取ることができ、背景の相関関数の上部に周期的な小さなスパイクを導く原因となる。このアーチファクトは、真の関数最大値の位置を決定する際の精度にマイナスの影響を及ぼす。この問題に対処するため、相関関数のローパスフィルタ処理が、最大値探索前に使用される。
y_offsetは、ある一定の範囲内にあるべきである、というのは、それは、2つのサブ画像間の水平変位の大きさがどのくらいあるかを指示するからである。この範囲は、カセット/スクリーンを保持するハードウェア固定具により定まる。y_offsetの値のチェックは、それ故に、実行され、y_offsetが、+/−70ピクセルのような、最大の許容範囲を超えた場合に、対応する仮想が拒絶される。
特別な状況は、c(Δ)が一定であることであり、これは、2つの連続サブ画像の重なり領域が相関しないことを意味する。対応した仮想は拒絶される。
最後に、各仮想に対して最小感度(フィガ・オブ・メリット)関数が確立され、最も高いスコアの仮想が、全体の検出プロセスに対する最も可能性のある候補として選択される。上述のプロセスの全ては個々の仮想の局所レベルで動作し、この段階のみが、全ての考えられる仮想を比較しようとする全体レベルで動作する。本発明の一実施例では、仮想の最小感度関数は、全ての連続サブ画像の対の正規化相関関数の最大値の合算値に基づいて構築される。特に以下の式で表わされる。
検出された仮想に基づいて、サブ画像は、必要な場合180°回転され、順に再配列され、次いで、共に結合される。
本発明は、デジタル化されている従来のX線写真画像を重ねることで得られるデジタル画像にも適用可能であることは、理解されるべきである。
101 カセット
102 記憶蛍光体スクリーン
600 フロント画像
601 スクリーン端部縁部
610 バック画像
102 記憶蛍光体スクリーン
600 フロント画像
601 スクリーン端部縁部
610 バック画像
Claims (3)
- 結合された完全な画像を構成するための、部分的な放射線画像の自動配列決定方法であって、
長い対象物の少なくとも2つの放射サブ画像を、重なり合う記録メディアにより、取得する取得ステップと、
前記記録されたサブ画像をデジタル画像の変換し、保存する変換ステップと、
サブ画像の順番、向き及び重なりの考えられる全ての固有の配列を含む仮想リストを構築するステップと、
各仮想における各連続したサブ画像の対に対して重なり領域を検出するステップと、
各仮想の連続したサブ画像の対毎に複数の異なる計測を実行し、計測結果が範囲外の場合に前記リストから当該仮想を消去するステップと、
前記仮想リストに残存する各仮想における各連続したサブ画像に対して重なり領域の相関関数を算出し、関数最大値及び対のサブ画像間の水平距離を求めるステップと、
連続したサブ画像の各対の水平距離の大きさを検出し、該大きさが範囲外の場合に前記仮想リストから当該仮想を消去するステップと、
相関関数の最大値の合算値に基づいて前記仮想リストに残存する各仮想に対して全体の最小感度を確立するステップと、
最大の最小感度の仮想を最も良好な候補に選択するステップと、
結合されたサブ画像からなる出力画像を生成するため、選択された仮想に基づいてサブ画像を結合処理するステップとを含む、方法。 - 前記取得ステップにおいて、前記長い対象物は、個人の足若しくは背骨を含む人体部である、請求項1記載の方法。
- 前記取得ステップは、重なる記憶蛍光体スクリーンにより行われ、前記変換ステップは、デジタルメモリに記憶されるデジタルサブ画像を生成する記憶蛍光体CRリーダにより実行される、請求項1記載の方法。
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US10/268,452 US6793390B2 (en) | 2002-10-10 | 2002-10-10 | Method for automatic arrangement determination of partial radiation images for reconstructing a stitched full image |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2004160208A true JP2004160208A (ja) | 2004-06-10 |
Family
ID=32030367
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2003348652A Withdrawn JP2004160208A (ja) | 2002-10-10 | 2003-10-07 | 結合フル画像再現用の部分放射線画像の自動配列決定方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US6793390B2 (ja) |
EP (1) | EP1408450A3 (ja) |
JP (1) | JP2004160208A (ja) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012045172A (ja) * | 2010-08-26 | 2012-03-08 | Fujifilm Corp | 放射線画像撮影システム、放射線画像撮影方法、及びプログラム |
JP2016032590A (ja) * | 2014-07-31 | 2016-03-10 | 富士フイルム株式会社 | 放射線画像撮影システム |
Families Citing this family (33)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7010152B2 (en) * | 2002-01-22 | 2006-03-07 | Canon Kabushiki Kaisha | Radiographic image composition and use |
US8194946B2 (en) * | 2005-07-28 | 2012-06-05 | Fujifilm Corporation | Aligning apparatus, aligning method, and the program |
CN101237818A (zh) * | 2005-08-08 | 2008-08-06 | 皇家飞利浦电子股份有限公司 | 用于固定焦点长格式数字放射线照相术的系统和方法 |
US7801351B2 (en) * | 2005-11-22 | 2010-09-21 | General Electric Company | Method and system to manage digital medical images |
DE102005058394A1 (de) * | 2005-12-07 | 2007-06-14 | Forschungsgesellschaft für Angewandte Naturwissenschaften e.V.(FGAN) | Verfahren und Einrichtung zur Stabilisierung von digitalen Bildsequenzen nach einer bildgestützten Transformationsvorschrift |
CN101479767A (zh) * | 2006-06-30 | 2009-07-08 | Nxp股份有限公司 | 视频缝合的方法和设备 |
US8098956B2 (en) | 2007-03-23 | 2012-01-17 | Vantana Medical Systems, Inc. | Digital microscope slide scanning system and methods |
US20110090327A1 (en) * | 2009-10-15 | 2011-04-21 | General Electric Company | System and method for imaging with enhanced depth of field |
US8314837B2 (en) * | 2009-10-15 | 2012-11-20 | General Electric Company | System and method for imaging with enhanced depth of field |
US8614484B2 (en) * | 2009-12-24 | 2013-12-24 | Taiwan Semiconductor Manufacturing Company, Ltd. | High voltage device with partial silicon germanium epi source/drain |
US9004757B2 (en) * | 2010-03-24 | 2015-04-14 | PaloDEx Grou Oy | Systems, assemblies, computer readable media and methods for medical imaging |
CN102855613B (zh) * | 2011-07-01 | 2016-03-02 | 株式会社东芝 | 图像处理设备和方法 |
US8774549B2 (en) * | 2011-09-30 | 2014-07-08 | Stmicroelectronics, Inc. | Compression error handling for temporal noise reduction |
JP6413534B2 (ja) * | 2014-09-17 | 2018-10-31 | コニカミノルタ株式会社 | 放射線画像撮影システム |
EP3023057B1 (en) * | 2014-11-19 | 2021-03-31 | Fujifilm Corp | RADIATION IMAGING DEVICE |
JP6072096B2 (ja) * | 2015-01-30 | 2017-02-01 | キヤノン株式会社 | 放射線撮影システム、制御方法、制御方法、及びプログラム |
JP6072102B2 (ja) * | 2015-01-30 | 2017-02-01 | キヤノン株式会社 | 放射線撮影システム及び放射線撮影方法 |
JP6072097B2 (ja) * | 2015-01-30 | 2017-02-01 | キヤノン株式会社 | 放射線撮影装置、制御装置、長尺撮影システム、制御方法、及びプログラム |
CN105832353B (zh) * | 2015-01-30 | 2020-11-06 | 佳能株式会社 | 放射线摄像系统 |
US10368823B2 (en) * | 2015-01-30 | 2019-08-06 | Canon Kabushiki Kaisha | Radiographing apparatus, control apparatus, control method, and storage medium |
JP6412815B2 (ja) * | 2015-02-26 | 2018-10-24 | 富士フイルム株式会社 | 放射線画像撮影システム、撮影台、及び撮影方法 |
GB2566748B (en) * | 2017-09-26 | 2022-08-17 | Focal Point Positioning Ltd | A method and system for calibrating a system parameter |
US10748319B1 (en) | 2016-09-19 | 2020-08-18 | Radlink, Inc. | Composite radiographic image that corrects effects of parallax distortion |
US10650561B2 (en) * | 2016-09-19 | 2020-05-12 | Radlink, Inc. | Composite radiographic image that corrects effects of parallax distortion |
WO2019076436A1 (en) | 2017-10-16 | 2019-04-25 | Hp Indigo B.V. | APPARATUS FOR PROCESSING IMAGES |
JP7050518B2 (ja) * | 2018-02-15 | 2022-04-08 | キヤノン株式会社 | 放射線撮影装置、放射線撮影システム、及び放射線撮影システムの制御方法 |
US10888294B2 (en) | 2018-12-27 | 2021-01-12 | Medtronic Navigation, Inc. | System and method for imaging a subject |
US10881371B2 (en) | 2018-12-27 | 2021-01-05 | Medtronic Navigation, Inc. | System and method for imaging a subject |
US11071507B2 (en) | 2018-12-27 | 2021-07-27 | Medtronic Navigation, Inc. | System and method for imaging a subject |
US11253212B2 (en) * | 2020-01-07 | 2022-02-22 | General Electric Company | Tileable X-ray detector cassettes |
WO2021168690A1 (en) * | 2020-02-26 | 2021-09-02 | Shenzhen Xpectvision Technology Co., Ltd. | Image sensors and methods of operating the same |
EP3996033A1 (en) * | 2020-11-06 | 2022-05-11 | Koninklijke Philips N.V. | Apparatus for scan direction and stiching sequence determination of a plurality of x-ray images |
CN113305017B (zh) * | 2021-05-28 | 2022-11-08 | 柳州源创电喷技术有限公司 | 全自动阀芯综合智能检测分选方法 |
Family Cites Families (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2565053B1 (fr) | 1984-05-25 | 1986-08-22 | Thomson Cgr | Procede de traitement d'images radiologiques |
US5528043A (en) * | 1995-04-21 | 1996-06-18 | Thermotrex Corporation | X-ray image sensor |
US5986279A (en) | 1997-03-21 | 1999-11-16 | Agfa-Gevaert | Method of recording and reading a radiation image of an elongate body |
DE69732829T2 (de) | 1997-03-21 | 2006-04-13 | Agfa-Gevaert | Verfahren zur Aufzeichnung und Wiedergabe eines Strahlungsbildes von einem länglichen Körper |
US6078701A (en) * | 1997-08-01 | 2000-06-20 | Sarnoff Corporation | Method and apparatus for performing local to global multiframe alignment to construct mosaic images |
EP0919856B1 (en) | 1997-12-01 | 2005-07-06 | Agfa-Gevaert | Method and assembly for recording a radiation image of an elongate body |
US6273606B1 (en) | 1997-12-01 | 2001-08-14 | Agfa-Gevaert | Method and assembly for recording a radiation image of an elongate body |
EP0919858B1 (en) | 1997-12-01 | 2004-08-25 | Agfa-Gevaert | Method for reconstructing a radiation image of a body from partial radiation images |
JP4076684B2 (ja) | 1998-12-15 | 2008-04-16 | 富士フイルム株式会社 | 画像合成システム |
JP2000241925A (ja) | 1999-02-18 | 2000-09-08 | Konica Corp | ハロゲン化銀写真乳剤の製造方法および該乳剤を用いたハロゲン化銀写真感光材料 |
JP2000250153A (ja) | 1999-03-04 | 2000-09-14 | Fuji Photo Film Co Ltd | 蓄積性蛍光体シート用カセッテ |
JP2000258861A (ja) | 1999-03-08 | 2000-09-22 | Fuji Photo Film Co Ltd | 蓄積性蛍光体シート及びこれを収容するカセッテ |
JP2000267210A (ja) | 1999-03-15 | 2000-09-29 | Fuji Photo Film Co Ltd | 蓄積性蛍光体シート用カセッテ |
US6459094B1 (en) * | 2000-12-20 | 2002-10-01 | Eastman Kodak Company | Method for stitching partial radiation images to reconstruct a full image |
US6895106B2 (en) * | 2001-09-11 | 2005-05-17 | Eastman Kodak Company | Method for stitching partial radiation images to reconstruct a full image |
-
2002
- 2002-10-10 US US10/268,452 patent/US6793390B2/en not_active Expired - Lifetime
-
2003
- 2003-09-29 EP EP03078092A patent/EP1408450A3/en not_active Withdrawn
- 2003-10-07 JP JP2003348652A patent/JP2004160208A/ja not_active Withdrawn
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012045172A (ja) * | 2010-08-26 | 2012-03-08 | Fujifilm Corp | 放射線画像撮影システム、放射線画像撮影方法、及びプログラム |
JP2016032590A (ja) * | 2014-07-31 | 2016-03-10 | 富士フイルム株式会社 | 放射線画像撮影システム |
US9697923B2 (en) | 2014-07-31 | 2017-07-04 | Fujifilm Corporation | Radiation image capturing system |
Also Published As
Publication number | Publication date |
---|---|
US20040071269A1 (en) | 2004-04-15 |
US6793390B2 (en) | 2004-09-21 |
EP1408450A3 (en) | 2009-07-01 |
EP1408450A2 (en) | 2004-04-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP2004160208A (ja) | 結合フル画像再現用の部分放射線画像の自動配列決定方法 | |
US6895106B2 (en) | Method for stitching partial radiation images to reconstruct a full image | |
US6775399B1 (en) | ROI segmentation image processing system | |
US6101238A (en) | System for generating a compound x-ray image for diagnosis | |
JP2008018047A (ja) | 欠陥領域補正装置および方法並びにプログラム、放射線検出装置 | |
EP1768069A1 (en) | Method of processing radiographic image and processing apparatus | |
JP4110074B2 (ja) | 放射線画像処理装置、放射線画像処理方法、プログラム及びコンピュータ可読媒体 | |
JP2013512037A (ja) | 一連のct射影画像においてx線マーカを追跡するための方法 | |
JPH07323024A (ja) | 画像診断支援装置 | |
JP5329077B2 (ja) | ハイブリッド医用撮像ユニットの作動方法およびハイブリッド医用撮像ユニット | |
JP4010375B2 (ja) | 放射線画像処理方法 | |
JP3888046B2 (ja) | 放射線画像処理方法および放射線画像処理装置 | |
JPH10268451A (ja) | 細長対象物の放射画像の記録及び読取方法 | |
JP4408863B2 (ja) | 医用画像処理装置及び方法 | |
JP6677256B2 (ja) | 読影支援プログラム、断面画像生成装置及び断面画像生成方法 | |
JP6642048B2 (ja) | 医療画像表示システム、医療画像表示プログラム及び医療画像表示方法 | |
JP3901188B2 (ja) | 放射線画像処理方法 | |
JP3731400B2 (ja) | 画像処理方法および画像処理装置 | |
Wang et al. | Fully automatic and reference-marker-free image stitching method for full-spine and full-leg imaging with computed radiography | |
JP2014054478A (ja) | 骨塩量計測装置および方法 | |
JP2006239270A (ja) | 異常陰影検出装置、異常陰影検出方法およびそのプログラム | |
WO2012111458A1 (ja) | X線画像診断装置及び画像表示方法 | |
JP2002143141A (ja) | 画像診断支援装置 | |
JP2004305486A (ja) | 画像処理装置及び方法及びシステム | |
JP4756756B2 (ja) | 画像処理方法、およびプログラム |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20061004 |
|
A711 | Notification of change in applicant |
Free format text: JAPANESE INTERMEDIATE CODE: A711 Effective date: 20080310 |
|
A761 | Written withdrawal of application |
Free format text: JAPANESE INTERMEDIATE CODE: A761 Effective date: 20081204 |