JP3895323B2 - 多視角x線立体結像の方法及びシステム - Google Patents
多視角x線立体結像の方法及びシステム Download PDFInfo
- Publication number
- JP3895323B2 JP3895323B2 JP2003432314A JP2003432314A JP3895323B2 JP 3895323 B2 JP3895323 B2 JP 3895323B2 JP 2003432314 A JP2003432314 A JP 2003432314A JP 2003432314 A JP2003432314 A JP 2003432314A JP 3895323 B2 JP3895323 B2 JP 3895323B2
- Authority
- JP
- Japan
- Prior art keywords
- image
- ray
- stereoscopic
- angle
- line
- 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.)
- Expired - Fee Related
Links
- 238000003384 imaging method Methods 0.000 title claims description 91
- 238000000034 method Methods 0.000 title claims description 43
- 230000005855 radiation Effects 0.000 claims description 48
- 238000012545 processing Methods 0.000 claims description 38
- 238000005070 sampling Methods 0.000 claims description 32
- 238000004364 calculation method Methods 0.000 claims description 28
- 230000015572 biosynthetic process Effects 0.000 claims description 26
- 238000003786 synthesis reaction Methods 0.000 claims description 23
- 238000010191 image analysis Methods 0.000 claims description 22
- 230000002452 interceptive effect Effects 0.000 claims description 15
- 238000006243 chemical reaction Methods 0.000 claims description 7
- 239000011521 glass Substances 0.000 claims description 6
- 238000004891 communication Methods 0.000 claims description 2
- 230000004044 response Effects 0.000 claims description 2
- 230000002194 synthesizing effect Effects 0.000 claims description 2
- 230000002708 enhancing effect Effects 0.000 claims 2
- 238000010187 selection method Methods 0.000 claims 2
- 239000007787 solid Substances 0.000 claims 2
- 230000003190 augmentative effect Effects 0.000 claims 1
- 238000001914 filtration Methods 0.000 claims 1
- 238000012952 Resampling Methods 0.000 description 12
- 230000000694 effects Effects 0.000 description 9
- 238000005259 measurement Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 238000007689 inspection Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 239000013078 crystal Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 210000000887 face Anatomy 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000009659 non-destructive testing Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000010845 search algorithm Methods 0.000 description 1
- 238000001308 synthesis method Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- 238000007794 visualization technique Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
- G03C—PHOTOSENSITIVE MATERIALS FOR PHOTOGRAPHIC PURPOSES; PHOTOGRAPHIC PROCESSES, e.g. CINE, X-RAY, COLOUR, STEREO-PHOTOGRAPHIC PROCESSES; AUXILIARY PROCESSES IN PHOTOGRAPHY
- G03C9/00—Stereo-photographic or similar processes
-
- G—PHYSICS
- G21—NUCLEAR PHYSICS; NUCLEAR ENGINEERING
- G21K—TECHNIQUES FOR HANDLING PARTICLES OR IONISING RADIATION NOT OTHERWISE PROVIDED FOR; IRRADIATION DEVICES; GAMMA RAY OR X-RAY MICROSCOPES
- G21K4/00—Conversion screens for the conversion of the spatial distribution of X-rays or particle radiation into visible images, e.g. fluoroscopic screens
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/022—Stereoscopic imaging
-
- G—PHYSICS
- G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
- G03C—PHOTOSENSITIVE MATERIALS FOR PHOTOGRAPHIC PURPOSES; PHOTOGRAPHIC PROCESSES, e.g. CINE, X-RAY, COLOUR, STEREO-PHOTOGRAPHIC PROCESSES; AUXILIARY PROCESSES IN PHOTOGRAPHY
- G03C5/00—Photographic processes or agents therefor; Regeneration of such processing agents
- G03C5/16—X-ray, infrared, or ultraviolet ray processes
- G03C2005/168—X-ray material or process
-
- G—PHYSICS
- G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
- G03C—PHOTOSENSITIVE MATERIALS FOR PHOTOGRAPHIC PURPOSES; PHOTOGRAPHIC PROCESSES, e.g. CINE, X-RAY, COLOUR, STEREO-PHOTOGRAPHIC PROCESSES; AUXILIARY PROCESSES IN PHOTOGRAPHY
- G03C5/00—Photographic processes or agents therefor; Regeneration of such processing agents
- G03C5/16—X-ray, infrared, or ultraviolet ray processes
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- High Energy & Nuclear Physics (AREA)
- Image Processing (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Measurement Of Radiation (AREA)
- Nuclear Medicine (AREA)
- Processing Or Creating Images (AREA)
- Image Analysis (AREA)
Description
1)測定或いは較正の方式によって結像システムパラメータを取得するステップと、
2)デジタル結像装置またはその中に置いた物体のいずれか一方を回転させて、前記デジタル結像装置と前記物体が相対的に円周または螺旋状の移動を行わせ、各θ度ごとに一枚の投影画像(Gk(θ)、θは任意数)をサンプリングする投影画像サンプリングステップと、
3)全ての画像画素に対して結像システムパラメータに従って2レベルに分類した索引を作成する画像索引形成ステップと、
4)ユーザがインタラクティブなインタフェースを通じて観察の必要に応じて、異なる角度から観察する効果のある立体画像対を取得することができるように立体画像対の視点パラメータを設定する視点パラメータ設定ステップと、
5)前記視点パラメータより現在の立体画像対を確定し、画像上における各画素に対応する視線パラメータを計算する視線パラメータ計算ステップと、
6)ステップ5において算出された前記視線パラメータによって、ステップ3に形成された画像索引表から前記視線パラメータに近い放射ビームを探し出す画像索引検索ステップと、
7)コンピュータの演算性能およびユーザから画像精度に対する要求に応じて、各種フィルタ補間方式を採用して、近隣放射ビームに対して補間合成の計算を行い、視線(L’ij)に対応する画像画素(p'(i,j))を合成し、ステップ3から5までを繰り返して、立体画像対上における全ての画素点について算出する画素合成ステップと、
8)ユーザの需要に応じてインタラクティブなインタフェースを通じて画像に対して増強処理を行う画像処理ステップと、
9)立体表示装置によって立体画像対の表示を行い、ユーザが左眼で左眼視角に対応する画像を見、右眼で右眼視角に対応する画像を見ることによって、立体画像を形成する立体表示ステップとを備える。
X線源(1)と平板X線検知器(3)からなるX線デジタル結像装置と、
多自由度回転可能載置台(2)と、
走査制御及びデータ獲得手段(4)と、
多自由度制御手段(5)と、
立体表示図形カード(7)と、
画像解析処理手段(6)と、
ディスプレー8と、
立体メガネ9とを備える多視角X線立体結像システムであって、
前記X線デジタル結像装置は、円周軌跡或いは螺旋軌跡に沿って走査することを実現でき、前記走査制御及びデータ獲得手段は、測定或いは較正の方式によって結像システムパラメータを取得し、デジタル結像装置またはその中に置いた物体のいずれか一方を回転させて、前記デジタル結像装置と前記物体を相対的に円周または螺旋状の移動を行わせ、各θ度ごとに一枚の投影画像(Gk(θ)、θが任意数値)をサンプリングし、前記画像解析処理手段は、サンプリングされた画像画素データを需要に応じてデータの圧縮を行ってコンピュータの内部メモリに記憶し、全ての画像画素が結像システムパラメータによって2レベルに分類した索引を作成する画像索引形成と、観察の必要に応じてユーザがインタラクティブなインタフェースを通じて、異なる角度から観察する効果のある立体画像対を取得することができるように立体画像対の視点パラメータを設定する視点パラメータ設定と、視点パラメータより現在の立体画像対を確定し、画像上における各画素に対応する視線パラメータを計算する視線パラメータ計算と、前記ステップにおいて前記算出された視線パラメータによって、前記画像索引表から前記視線パラメータに近い放射ビームを探し出す画像索引検索と、コンピュータの演算性能およびユーザから画像精度に対する要求に応じて、各種フィルタ補間方式を採用して、近隣放射ビームに対して補間合成の計算を行い、視線(L’ij)に対応する画像画素(p'(i,j))を合成し、立体画像対上における全ての画素点について計算を行う画素合成と、ユーザの要求に応じてインタラクティブなインタフェースを通じて画像に対して増強処理を行う画像処理と、立体表示装置によって立体画像対の表示を行い、ユーザが左眼で左眼視角に対応する画像を見させ、右眼で右眼視角に対応する画像を見させることによって、立体画像を形成する立体表示とを行うものである。
本発明の多視角X線立体結像方法は、まず、X線結像システムパラメータとデータ収集パラメータに基づいて円周軌跡及び螺旋軌跡に沿って収集された一連の投影画像データの画素列に対して、角度と距離の索引を確立し、その両索引を素早く検索することができるように順序を付けて記憶し、選択された観察視点や視線方向、視差効果によって、合成すべき立体画像対における画素の角度と距離のパラメータを計算し、記憶されている両索引とその対応する画像データを取り出して画像の合成を行う。本発明における多視角X線立体画像システムは、円周軌跡または螺旋軌跡の走査を実現させることができるX線デジタル結像装置と、コンピュータ装置と、立体デススプレイ装置などを備える。X線デジタル結像装置からX線画像をコンピュータ装置に入力させ、このコンピュータ装置が本発明の多視角立体結像方法によって立体画像を合成し、立体表示装置によって立体画像を表示する。周知のコンピュータ装置は、走査制御及びデータ取得手段と、画像解析処理手段とを含む。
di=D×sin(βi) (3)
ただしDはX線源と物体回転軸間の距離であり、すなわち、サンプリング視点と物体回転軸間の距離であるが、直角三角形の斜辺でもある。βiと同じように、diの符号の正負は射線Lijの回転中心軸における左側或いは右側を示している。本例において、di<0の時、Lijはサンプリング中心射線の左側にあり、di>0の時、Lijはサンプリング中心射線の右側にある。
α'i=ω−(φ+β'i) (5)
ただし、ωは、視点回転方位角で、計算結果からα'iは[0,2π]範囲に正規化される。
d'i=r×sin(φ+β'i) (6)
Pij=r1P11+r2P12+r3P21+r4P22 (7)
ただし、重みファクタr1、r2、r3、r4は、
1)システムパラメータ測定と較正:測定或いは較正方法によって結像システムにおけるパラメータ、すなわちL、D、P0(i0,j0),d0を取得する。
2)円周又は螺旋軌跡投影画像サンプリング:デジタル結像装置とその中に置いた物体を相対的に円周または螺旋状に移動させるようにデジタル結像装置或いは物体のいずれか一方を回転させ、θ度毎に一枚の投影Gk(θ)をサンプリングする。θは任意数値である。
3)2レベルに分類した索引の計算:サンプリングされた画像データをコンピュータの内部メモリに記憶して、必要に応じてデータの圧縮を行う。コンピュータに、全ての画像画素に対して結像システムのパラメータによって角度と距離の索引を作成する。具体的方法として、第k番目の投影画像のサンプリング点は、Pkに位置し、その回転角θは、サンプリング中心射線(すなわち、物体回転中心の射線)と回転平面内にある基準座標軸との平面角であり、Dは、サンプリング円周の半径である。投影画像Gk(θ)上における画素p(i,j)に対応する放射ビームはLijで、結像システムパラメータによってLijとサンプリング中心線との角βi、Lijと基準座標軸との平面角αi、及びLijから物体回転中心までの距離diなどを算出することができる。各画像の各列に対応するαiとdiは2レベルに分類した索引中に記憶して置く。
4)視点パラメータのセット:ユーザは、インタラクティブなインタフェースを通して観察の必要に応じて立体画像対の視点パラメータを設定して、異なる角度から観察する効果の立体画像対を取得する。インタラクティブ操作によってユーザより設定できるパラメータとしては、視点位置、観察視線方向、視景体(body of visual scene)、立体表示の視差パラメータなどがある。
5)視線パラメータの計算:視線パラメータによって現在の立体画像対を決定して、画像上における各画素毎に対応する視線パラメータを計算する。画素点p'(i,j)に対応する視線は、L’ijであり、視線L’ijと基準座標軸との平面角α'i、L’ijから物体回転中心までの距離d'iを計算し、又物体回転中心から源の距離よりそのL’ij、高さに対応する画像行jを計算する。視点パラメータによって左右の眼に対応する立体画像対を計算し、いずれも同じ方法を採用するが、両視点が垂直視線中心方向に距離tcだけずれている。
6)索引表シーク:前記ステップで計算から得られたL’ijと基準座標軸との平面角α'i、L’ijから物体回転中心までの距離d'iによって形成された画像索引表からL’ijパラメータと近隣している放射ビームを探す。隣接放射ビームの数量及び選択方式はステップ7で選ぶ画像画素合成フィルタ方法と関連する。
7)画素補間合成:コンピュータの演算性能およびユーザから画像精度に対する要求に応じて、各種フィルタ補間方式を採用して、近隣放射ビームに対して補間合成の計算を行って、L’ijに対応する画像画素p'(i,j)を合成する。通常のフィルタ補間方法には、領域補間或いは双一次補間、二次補間などがある。ステップ5から7までを繰り返して、立体画像対上における全ての画素点の計算が完了する。
8)画像処理:ユーザの需要に応じて、立体表示効果の強調のために、インタラクティブなインタフェースによって、輝度変換、擬似カラー、エッジ強調、画像増強などの処理の機能が提供される。
9)立体表示:立体表示装置によって立体画像対の表示を行い、ユーザが左眼で左眼視角に対応する画像を見、右眼で右眼視角に対応する画像を見ることによって、立体画像を形成する。ステップ3から9までを繰り返すとユーザ観察視点位置、観察視線方向、視差を絶えずに変えて、多角度立体表示の検査効果を達することができる。
Claims (21)
- 1)X線源を調整して面アレー検知器に直角な放射が被検体(人を除く)の回転軸を通過し、この放射の面アレー検知器上における位置を画素位置P0(i0,j0)として、X線源と面アレー検知器との距離Lと、X線源と被検体(人を除く)回転中心との距離Dと、面アレー検知器における画素間の距離をd0(以下、これらのパラメータP0(i0,j0)、L、D、d0を「結像システムパラメータ」という)を取得するステップと、
2)デジタル結像装置または被検体(人を除く)のいずれか一方を回転させて、前記デジタル結像装置と前記被検体に相対的に円周または螺旋状の移動を行わせ、各θ度毎に各一枚投影画像(Gk(θ)、θは任意数値)を取得する投影画像サンプリングステップと、
3)結像システムパラメータに従って全ての投影画像画素に対して第1レベルの索引として放射ビームから回転中心軸までの距離diと第2レベルの索引として投影角度αiを含む索引(以下、「2レベルに分類した画像索引」という)を作成する画像索引形成ステップと、
4)ユーザがインタラクティブなインタフェースを通じて、観察しながら、異なる角度から観察する効果のある立体画像対を取得することができるように立体画像対の視点パラメータを設定する視点パラメータ設定ステップと、
5)前記視点パラメータにより現在の立体画像対を確定し、画像上における各画素に対応する視線パラメータを計算する視線パラメータ計算ステップと、
6)ステップ5において算出された前記視線パラメータに基づいて、ステップ3で形成された2レベルに分類した画像索引から前記視線パラメータに近い順に複数の放射ビームを探し出す画像索引検索ステップと、
7)フィルタ補間方式を採用して、ステップ6で探し出された複数の放射ビームに対して補間合成の計算を行い、視線(L’ij)に対応する画像画素(p'(i,j))を合成し、立体画像対上における全ての画像画素についてステップ3から5までを繰り返して画素を合成する画素合成ステップと、
8)ユーザの要求に応じてインタラクティブなインタフェースを通じて画像に対して増強処理を行う画像処理ステップと、
9)立体表示装置によって立体画像対の表示を行い、ユーザに左眼で左眼視角に対応する画像を見させ、右眼で右眼視角に対応する画像を見させることによって、立体画像を形成する立体表示ステップと、
を備えることを特徴とする多視角X線立体結像方法。 - 前記ステップ3における前記2レベルに分類した画像索引は、各画像の各列に対応する投影角度αiと放射ビームから回転中心軸までの距離diの索引であることを特徴とする請求項1に記載の多視角X線立体結像方法。
- 前記ステップ3が2レベルに分類した画像索引を計算するステップをさらに含むことを特徴とする請求項1に記載の多視角X線立体結像方法。
- 前記ステップ3は、更に、その回転角度(θ)が回転中心射線に対応したサンプリング中心射線と回転平面内に位置する基準座標軸との平面角であり、投影画像Gk(θ)上における画素(P(i,j))に対応する放射ビームをLijとして、第k番目の投影画像サンプリング点を半径がDのサンプリング円周のサンプリング視点(Pk)に位置せしめ、
結像システムパラメータに基づいて、前記放射ビーム(Lij)とサンプリング中心線との平面角(βi)と、前記放射ビーム(Lij)と基準座標軸との平面角(αi)並びに前記放射ビーム(Lij)から回転中心までの距離(di)を算出し、
各画像の各列に対応する前記平面角(αi)と前記距離(di)を前記の2レベルに分類した画像索引に記憶して索引を作成するステップを含むことを特徴とする請求項3に記載の多視角X線立体結像方法。 - 前記視点パラメータ設定ステップ4において、インタラクティブ操作によってユーザにより設定可能なパラメータが観察視点位置、観察視線方向、視景体と立体表示の視差パラメータであることを特徴とする請求項1に記載の多視角X線立体結像方法。
- 前記視線パラメータ計算ステップ5が、
前記画像画素(p'(i,j))に対応する視線(L’ij)と基準座標軸中心線との平面角(α'i)と、前記視線(L’ij) から回転中心までの距離(d'i)とを計算し、
回転中心からX線源までの距離を用いて前記の視線(L’ij) 高さに対応する画像の行(j)を計算し、
左右眼に対応する画像対を計算するのに、同じ方法を採用して、該両視点を垂直視線中心方向における距離(tc)だけずらして、視点パラメータに基づいて左右眼に対応する画像対を計算するステップを含むことを特徴とする請求項1に記載の多視角X線立体結像方法。 - 前記ステップ6において、隣接放射ビームの数と選択方式が前記のステップ7で選ばれる画像画素合成のフィルタ方法に関連することを特徴とする請求項1に記載の多視角X線立体結像方法。
- 前記ステップ7における前記フィルタ補間方法は近接補間、双一次補間、二次補間を含むことを特徴とする請求項1に記載の多視角X線立体結像方法。
- 前記ステップ8における前記画像処理は、輝度変換、擬似カラー、エッジ強調の立体表示を増強する方法を含むことを特徴とする請求項1に記載の多視角X線立体結像方法。
- 前記ステップ2からステップ7までを繰り返し実行し、多角度立体表示を達成するように、ユーザが視点位置、観察視線方向、視差を絶えず観察しながら変更することを特徴とする請求項1に記載の多視角X線立体結像方法。
- X線源(1)と平板X線検知器(3)を有するX線デジタル結像装置と、
多自由度で回転可能な載置台(2)と、
走査制御及びデータ獲得手段(4)と、
多自由度制御手段(5)と、
立体表示図形カード(7)と、
画像解析処理手段(6)と、
ディスプレー(8)と、
立体メガネ(9)とを備える多視角X線立体結像システムであって、
前記X線デジタル結像装置は円周軌跡或いは螺旋軌跡に沿って走査することが実現でき、
前記走査制御及びデータ獲得手段は、X線源を調整して面アレー検知器に直角な放射が被検体の回転軸を通過し、この放射の面アレー検知器上における位置を画素位置P0(i0,j0)として、X線源と面アレー検知器との距離Lと、X線源と被検体回転中心との距離Dと、面アレー検知器における画素間の距離をd0(以下、これらのパラメータP0(i0,j0)、L、D、d0を「結像システムパラメータ」という)を取得し、デジタル結像装置または被検体のいずれか一方を回転させて、前記デジタル結像装置と前記被検体が相対的に円周または螺旋状の移動を行い、各θ度毎に各一枚の投影画像(Gk(θ)、θは任意数値)を取得し、
前記画像解析処理手段は、画像索引形成と、視点パラメータ設定と、視線パラメータ計算と、画像索引検索と、画素合成と、画像処理と、立体表示のために用いられ、
画像索引形成では、結像システムパラメータに基づいて、全ての画像画素に対して第1レベルの索引として放射ビームから回転中心軸までの距離diと第2レベルの索引として投影角度αiを含む索引(以下、「2レベルに分類した画像索引」という)を作成し、
視点パラメータ設定では、ユーザがインタラクティブなインタフェースを通じて、観察しながら、異なる角度から観察する効果のある立体画像対を取得することができるように立体画像対の視点パラメータを設定し、
視線パラメータ計算では、前記視点パラメータにより現在の立体画像対を確定し、画像上における各画素に対応する視線パラメータを計算し、
画像索引検索では、視線パラメータを計算するステップで算出された前記視線パラメータに基づいて、2レベルに分類した画像索引から前記視線パラメータに近い順に複数の放射ビームを探し出し、
画素合成では、フィルタ補間方式を採用して、前記画像索引検索で探し出された複数の放射ビームに対して補間合成の計算を行い、視線(L’ij)に対応する画像画素(p'(i,j))を合成し、それにより立体画像対上における全ての画像画素について画素を合成し、
画像処理では、ユーザの要求に応じてインタラクティブなインタフェースを通じて画像に対して増強処理を行い、
立体表示では、立体表示装置によって立体画像対の表示を行い、ユーザに左眼で左眼視角に対応する画像を見させ、右眼で右眼視角に対応する画像を見させることによって、立体画像を形成することを特徴とする多視角X線立体結像システム。 - 前記載置台(2)は、前記X線源(1)と平板X線検知器(3)の間に設置され、
前記走査制御及びデータ獲得手段(4)は、データサンプリングカードを介して平板X線検知器(3)から送られた画像データを取得し、かつ、通信ポートによって多自由度制御手段(5)から送られた載置台(2)の走査位置情報を取得し、
前記画像解析処理手段(6)は、走査制御及びデータ獲得手段(4)から送られたデータによって画像処理及び立体画像の合成を完成し、かつ、立体表示図形カード(7)によって立体画像対を反転させてディスプレー(8)の上に表示し、また、立体メガネ(9)を駆動することを特徴とする請求項11に記載の多視角X線立体結像システム。 - 結像システムパラメータに基づき、前記の画像解析処理手段(6)によって形成される前記2レベルに分類した画像索引は、各画像の各列に対応する投影角度αiと放射ビームから回転中心軸までの距離diの索引であることを特徴とする請求項11に記載の多視角X線立体結像システム。
- 前記の画像解析処理手段(6)は、さらに、前記2レベルに分類した画像索引を計算することを特徴とする請求項12に記載の多視角X線立体結像システム。
- 前記の画像解析処理手段(6)は、
第k番目の投影画像サンプリング点が半径がDのサンプリング円周のサンプリング視点(Pk)に位置し、
その回転角度(θ)が回転中心射線に対応したサンプリング中心射線と回転平面内に位置する基準座標軸との平面角であり、
投影画像Gk(θ)上における画素(P(i,j))に対応する放射ビームをLijとして、
結像システムパラメータに基づいて、前記放射ビーム(Lij)とサンプリング中心線との平面角(βi)と、前記放射ビーム(Lij)と基準座標軸との平面角(αi)並びに前記放射ビーム(Lij)から回転中心までの距離(di)を算出し、
各画像の各列に対応する前記平面角(αi)と前記距離(di)を前記の2レベルに分類した画像索引に記憶して索引を作成することを特徴とする請求項13に記載の多視角X線立体結像システム。 - 前記の画像解析処理手段(6)においてインタラクティブ操作によってユーザにより設定可能なパラメータが観察視点位置、観察視線方向、視景体と立体表示の視差パラメータであることを特徴とする請求項12に記載の多視角X線立体結像システム。
- 前記の画像解析処理手段(6)は、
前記画像画素(p'(i,j))に対応する視線(L’ij)と、基準座標軸中心線との平面角(α'i)と、前記視線(L’ij) から回転中心までの距離(d'i)とを計算し、
回転中心からX線源までの距離を用いて前記の視線(L’ij) 高さに対応する画像の行(j)を計算し、
左右眼に対応する画像対を計算するのに、同じ方法を採用して、該両視点を垂直視線中心方向における距離(tc)だけずらして、視点パラメータに基づいて左右眼に対応する画像対を計算する仕方で視線パラメータを設定することを特徴とする請求項11に記載の多視角X線立体結像システム。 - 前記の画像解析処理手段(6)で、隣接放射ビームの数と選択方式は、選ばれた前記画像画素合成のフィルタ方法に関連することを特徴とする請求項11に記載の多視角X線立体結像システム。
- 前記の画像解析処理手段(6)によって行う前記フィルタ補間方法は、近接補間、双一次補間、二次補間を含むことを特徴とする請求項11に記載の多視角X線立体結像システム。
- 前記の画像解析処理手段(6)によって実行する前記画像処理は、輝度変換、擬似カラー、エッジ強調の立体表示を増強する処理を含むことを特徴とする請求項11に記載の多視角X線立体結像システム。
- 前記の画像解析処理手段(6)は、投影画像の再サンプリングと、画像索引作成と、視点パラメータ設定と、画像索引検索と、画素合成とを繰り返し実行し、ユーザが観察視点位置、観察視線方向、視差を絶えず観察しながら変更し、多角度立体表示を達成することを特徴とする請求項11に記載の多視角X線立体結像システム。
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN02159995 | 2002-12-31 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2004333473A JP2004333473A (ja) | 2004-11-25 |
JP3895323B2 true JP3895323B2 (ja) | 2007-03-22 |
Family
ID=33136777
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2003432314A Expired - Fee Related JP3895323B2 (ja) | 2002-12-31 | 2003-12-26 | 多視角x線立体結像の方法及びシステム |
Country Status (5)
Country | Link |
---|---|
US (1) | US6944263B2 (ja) |
JP (1) | JP3895323B2 (ja) |
KR (1) | KR100576228B1 (ja) |
DE (1) | DE10361663A1 (ja) |
RU (1) | RU2253952C1 (ja) |
Families Citing this family (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7456881B2 (en) * | 2006-01-12 | 2008-11-25 | Aptina Imaging Corporation | Method and apparatus for producing Bayer color mosaic interpolation for imagers |
US8137976B2 (en) | 2006-07-12 | 2012-03-20 | Varian Medical Systems, Inc. | Dual angle radiation scanning of objects |
US7692156B1 (en) | 2006-08-23 | 2010-04-06 | Radiation Monitoring Devices, Inc. | Beam-oriented pixellated scintillators for radiation imaging |
JP5498788B2 (ja) * | 2006-08-25 | 2014-05-21 | コーニンクレッカ フィリップス エヌ ヴェ | マルチ管x線検出 |
GB0902138D0 (en) * | 2009-02-10 | 2009-03-25 | Durham Scient Crystals Ltd | Apparatus and method for viewing an object |
JP5677738B2 (ja) * | 2009-12-24 | 2015-02-25 | 株式会社東芝 | X線コンピュータ断層撮影装置 |
GB201008281D0 (en) | 2010-05-19 | 2010-06-30 | Nikonovas Arkadijus | Indirect analysis and manipulation of objects |
DE102010056042A1 (de) * | 2010-12-23 | 2012-06-28 | Yxlon International Gmbh | Verfahren und Vorrichtung zur Sichtprüfung eines mittels Röntgenstrahlung zu überprüfenden Prüfobjekts |
WO2013046709A1 (ja) * | 2011-09-30 | 2013-04-04 | 富士フイルム株式会社 | 放射線画像表示方法および装置 |
CN103197773A (zh) * | 2012-01-09 | 2013-07-10 | 西安智意能电子科技有限公司 | 一种用于检测输入设备的三维位置信息的方法与系统 |
CN103764033B (zh) * | 2012-06-28 | 2017-03-08 | 东芝医疗系统株式会社 | X射线诊断装置 |
US8864763B2 (en) | 2013-03-13 | 2014-10-21 | DePuy Synthes Products, LLC | External bone fixation device |
EP3281592B1 (en) | 2013-03-13 | 2021-04-21 | DePuy Synthes Products, Inc. | External bone fixation device |
US9039706B2 (en) | 2013-03-13 | 2015-05-26 | DePuy Synthes Products, Inc. | External bone fixation device |
JP5751556B2 (ja) * | 2013-03-19 | 2015-07-22 | 株式会社ビームセンス | ステレオ透視装置およびそれを用いたステレオ観察方法 |
CN105784741A (zh) * | 2015-07-10 | 2016-07-20 | 公安部第研究所 | 一种多视角x射线安全检查装置 |
US10835318B2 (en) | 2016-08-25 | 2020-11-17 | DePuy Synthes Products, Inc. | Orthopedic fixation control and manipulation |
CN106725559A (zh) * | 2017-03-17 | 2017-05-31 | 青岛大学 | 一种实时双源x线3d数字成像装置 |
US11439436B2 (en) | 2019-03-18 | 2022-09-13 | Synthes Gmbh | Orthopedic fixation strut swapping |
US11304757B2 (en) | 2019-03-28 | 2022-04-19 | Synthes Gmbh | Orthopedic fixation control and visualization |
US11334997B2 (en) | 2020-04-03 | 2022-05-17 | Synthes Gmbh | Hinge detection for orthopedic fixation |
CN114199905B (zh) * | 2021-12-13 | 2024-02-20 | 中国航发南方工业有限公司 | 一种机匣内部缺陷的空间定位方法及系统 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4672649A (en) * | 1984-05-29 | 1987-06-09 | Imatron, Inc. | Three dimensional scanned projection radiography using high speed computed tomographic scanning system |
-
2003
- 2003-12-22 US US10/743,277 patent/US6944263B2/en not_active Expired - Lifetime
- 2003-12-26 RU RU2003137405/28A patent/RU2253952C1/ru active
- 2003-12-26 JP JP2003432314A patent/JP3895323B2/ja not_active Expired - Fee Related
- 2003-12-30 DE DE10361663A patent/DE10361663A1/de not_active Ceased
- 2003-12-31 KR KR1020030102036A patent/KR100576228B1/ko active IP Right Grant
Also Published As
Publication number | Publication date |
---|---|
JP2004333473A (ja) | 2004-11-25 |
US6944263B2 (en) | 2005-09-13 |
KR100576228B1 (ko) | 2006-05-03 |
RU2253952C1 (ru) | 2005-06-10 |
KR20040062414A (ko) | 2004-07-07 |
DE10361663A1 (de) | 2004-12-30 |
US20040208279A1 (en) | 2004-10-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP3895323B2 (ja) | 多視角x線立体結像の方法及びシステム | |
EP1522051B1 (en) | Discrete linear space sampling method and apparatus for generating digital 3d models | |
US6445807B1 (en) | Image processing method and apparatus | |
AU2011206927B2 (en) | A computed tomography imaging process and system | |
CN1264062C (zh) | 一种多视角x射线立体成像的方法与系统 | |
JP2744490B2 (ja) | 物体内部構造表面の2次元像を表示する装置と方法 | |
AU2011312140B2 (en) | Rapid 3D modeling | |
JP2885842B2 (ja) | 固体内部領域の切断平面を表示する装置と方法 | |
US9117278B2 (en) | System and method for geometric modeling using multiple data acquisition means | |
CN101371786B (zh) | 一种x射线图像三维重构的方法及系统 | |
US20130249901A1 (en) | Systems and methods for geometrically mapping two-dimensional images to three-dimensional surfaces | |
JP4914560B2 (ja) | 環境内から生じる放射線についての情報を得るための方法 | |
KR20040072466A (ko) | 3차원 백 프로젝션 방법 및 x-레이 ct 장치 | |
US6115449A (en) | Apparatus for quantitative stereoscopic radiography | |
CN104374374B (zh) | 基于主动全景视觉的3d环境复制系统及3d全景显示绘制方法 | |
JP4242527B2 (ja) | 立体測定データにおける表面を表示する方法及びシステム | |
Jeschke | Digital close-range photogrammetry for surface measurement | |
CN116804537A (zh) | 一种双目测距系统及方法 | |
WO2013142819A1 (en) | Systems and methods for geometrically mapping two-dimensional images to three-dimensional surfaces | |
JP3637416B2 (ja) | 3次元計測方法、3次元計測システム、画像処理装置、及びコンピュータプログラム | |
Meierhold et al. | Line-based referencing between images and laser scanner data for image-based point cloud interpretation in a CAD-environment | |
Jiang et al. | 3D Surface Reconstruction of Spinal Ultrasound Images Based on 2D Contour Line | |
CN1068670A (zh) | 用定位x线片重建活体颅骨立体形态方法 | |
Xu et al. | A camera network for the voxel data acquiring of the Three-Dimensional Swept Volume Display | |
Yoshimoto et al. | Evaluation of endoscopic entire 3D image acquisition of the digestive tract using a stereo endoscope |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20050621 |
|
A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20050921 |
|
A602 | Written permission of extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A602 Effective date: 20050927 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20051221 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20060228 |
|
A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20060529 |
|
A602 | Written permission of extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A602 Effective date: 20060607 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20060828 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20060926 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20061012 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20061107 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20061128 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20061213 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 Ref document number: 3895323 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101222 Year of fee payment: 4 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111222 Year of fee payment: 5 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121222 Year of fee payment: 6 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20131222 Year of fee payment: 7 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
LAPS | Cancellation because of no payment of annual fees |