JP4190194B2 - 断層撮影法によって得た像から３次元モデルを再構成するための較正方法 - Google Patents

Description

【０００１】
【発明の属する技術分野】
本発明は、Ｘ線断層撮影法によって得た像から３次元モデルを画像再構成するための較正方法に関する。
【０００２】
【発明の背景】
従来、断層撮影システムは、図１に示されているように、固定してもよい（例えば、Ｘ線に感応するスクリーンのような）検出手段１と、移動可能であって、検出手段１に対して且つ放射線検査を受けるべき人体又は物体に対して幾つかの位置を取り、また該人体又は物体の固定の点を中心にして回転する（例えば、Ｘ線源のような）放射線供給手段２とを含んでいる。
【０００３】
断層撮影２次元像から３次元モデルを再構成するための多数の手法が知られている。例えば、文献「Image reconstruction from projections: the fundamentales of computerized tomography 」, T. Herman, Academic Press, New York (1980)に記載の再構成アルゴリズムを参照されたい。このような再構成法は一般的に撮像システムの精密な「幾何学的較正」を必要とし、この較正は３次元空間と様々な２次元投影によって得られた２次元情報とを結びつける。もしこの較正が精確でないと、再構成された３次元モデルの品質に欠陥を有し、特に、小さな構造が明瞭にならない。場合によっては、この較正は、線源と検出手段との間の距離、線源の角度位置などのような、撮像システム自身によって供給される情報から直接的に実行される。この形式の較正の代替法では、（患者を含めない）撮像用幾何学的配置の事前較正を行い、また線源を予め較正した相次ぐ位置に置くことを含んでいる。これらの形式の較正では、機械的な理由のため、満足のいく精確さ又は結果が得られない。
【０００４】
他の較正手法では、空間内で基準として作用するマーカーを使用しており、これらのマーカーは撮像の際にＸ線の場の中にあり、従って投影像上に現れる。３次元空間内でのこれらのマーカーの位置は既知であると想定されており、各投影についての撮像用幾何学的配置は、投影像上のマーカーの位置から導き出された連立方程式の置き換えによって推測される。理論的には、これらのマーカーを用いる手法は、マーカーを用いない手法よりもより良い精確さが得られるはずである。実際には、空間内のマーカーの位置を精確に決定することは困難な場合が多い。特に、乳房撮影検査の際、マーカーは圧迫板に対して固定されているが、検出器に対する圧迫板の位置が高精度で分かっているわけではなく、また圧迫板は患者の胸に押し付けられたとき僅かに動き易い。
【０００５】
【発明の概要】
本発明の一実施形態は、マーカーを使用して、空間内でのそれらの位置を実質的に精確に決定する方法である。
【０００６】
本発明の一実施形態は、検出手段と該検出手段に対して移動可能である放射線供給手段とを含んでいる断層撮影装置によって得た像から３次元モデルを再構成するための較正方法であり、該方法では、１つの同じ支持体上に担持されている一組のマーカーを放射線の場の中に配置し、撮像した像へのマーカーの投影の位置を処理して、撮像時における放射線供給手段の位置を推測する。また、マーカーを支持体に対して固定し、撮像した像へのマーカーの投影の位置を処理することにより、空間内でのマーカーの位置を推測する。
【０００７】
特に、本発明の一実施形態では、像へのマーカーの投影の位置を処理して、これから、マーカーの想定した初期位置と変位した後にマーカーが取る位置との間のマーカーの剛性変位を特徴付ける変換パラメータを推測する。
【０００８】
特に、一実施形態によれば、様々な撮像の際の線源の位置を決定すると共に、全てのマーカー及び撮像について、一方では撮像から生じる像への各マーカーの投影の理論的位置と他方では像内でのその実際の位置との間の距離の関数である総合誤差を最小にする変換パラメータの値を決定する。
【０００９】
本方法の実施形態により、最初に空間内におけるマーカーの精確な位置を知ることを必要とせずに、実質的に精確な較正を得ることができることが理解されよう。
【００１０】
【発明の詳しい説明】
本発明の実施形態は、単に例示し且つそれに限定するものではない以下の説明を、添付の図面を参照して読むことにより、更に明らかになろう。
【００１１】
図２及び図３に示されるように、本発明を実施するためのマーカー４ａ〜４ｄに対する可能な１つの支持体は、平坦な板（例えば、圧迫板）３である。マーカー４ａ〜４ｄは例えば金属の球であり、板３はＸ線透過性の材料で作られる。撮像の際、撮像対象物（例えば、乳房）はスクリーン１とマーカー４ａ〜４ｄを支持している板３との間に圧迫される。線源２の各位置において、線源から放出された放射線が、検出手段を構成するスクリーン１上へ、放射線検査を受ける対象物の領域の像と同時にマーカー４ａ〜４ｄの像を投影する。板の正確な位置は厳密には分からない。特に、スクリーン１と板３とを乳房に対して配置するときに加える圧迫により板の変位が生じ、この変位は患者毎に大幅に変わり得る。既知であるデータは、様々なマーカー４ａ〜４ｄの互いに対する相対位置だけである。一旦圧迫板３が所定位置に置かれると、それは実行される一連の撮像全体を通じて固定された状態に留まると考えられる。
【００１２】
較正は、各々の撮像毎に、線源２の位置と、板３の初期位置に対する変位とを実質的に精確に決定することを含む。これは、Ｎを撮像数として、３Ｎ＋６個の撮像用幾何学的配置パラメータ、すなわち、各々の撮像についての線源２の３つの座標（これらの座標は、デカルト座標（図におけるｘ，ｙ，ｚ）又は極座標とすることができる）と、板の初期位置から実際の位置までの変位に対応する剛性変換を記述する６個のパラメータ（３つの平行移動パラメータ及び３つの回転パラメータ）とを決定することに対応する。上記のパラメータの決定を行うために、マーカーが互いに対して固定されていると云う制約を用いる。
【００１３】
例えば、全てのマーカー及び全ての撮像について、撮像から生じる像への各マーカーの投影の理論的位置Ｐｍと像内でのその実際の位置Ｐｄとの間の距離の和を最小にすることのできるような撮像パラメータが決定される。理論的位置Ｐｍは次のように定義される。
【００１４】
（１）線源の位置と仮定した位置（ｘ，ｙ，ｚ）から、線源の場の中にある点の理論的投影の法則を定義する。このような法則は、例えば、ビューイング・ソフトウエアにおける従来の場合と同様に、３×４投影マトリクスと要約される；
（２）板の変位に対応すると想定される剛体変換を全てのマーカーに適用することによって、各マーカーの位置を決定する；
（３）計算された投影をマーカーの計算された位置に適用し、これから、投影された像上のマーカーの理論的な位置Ｐｍを推測する。従来の非線形最小化手法、例えば、パウエル法によって、最小化を達成し得る。
【００１５】
多数の未知のパラメータがある場合、計算時間が望ましくないほど長くなることがある。その解決のために、上記の最小化が、好ましくは、２つの異なるステップにおいて、一方では剛体変換のパラメータを決定し、且つ他方では線源の取る様々な位置（ｘｉ，ｙｉ，ｚｉ）に対応する座標を決定することによって、実行される。第１に、剛体変換のパラメータに初期化値を与え、且つ総合誤差を最小にする全ての撮像についての線源位置を決定する。第２に、このように第１のステップで決定された線源位置を用いて、総合誤差を最小にする６つの剛性変換パラメータを決定する。
【００１６】
この反復処理において、様々な撮像の際に線源の様々な位置を決定することを含むステップがｎ個の独立の最小化に分解されるようにすることができる。この較正方法では、空間内のマーカーの位置を決定する。マーカーが圧迫板（乳房撮影検査の場合）に固定されているとき、上述の方法は胸部を圧迫する際の板の動きに対してかなり不感応である。
【００１７】
当業者にとって、特許請求の範囲に記載したような本発明の趣旨及び範囲から離れることなく構造や処理ステップや機能に様々な変更をなすことができよう。
【図面の簡単な説明】
【図１】断層撮影装置の一般的な原理を例示する略図である。
【図２】本発明の可能な一実施形態を示す略図である。
【図３】本発明の可能な一実施形態を示す略図である。
【符号の説明】
１ スクリーン
２ 放射線源
３ 圧迫板
４ａ〜４ｄ マーカー

Claims (10)

1. 検出手段と、該検出手段に対して移動可能である放射線源を構成する手段と、撮影対象物を前記検出手段に向かって圧迫する支持体（３）とを含んでいる断層撮影システムによって撮像した像から３次元モデルを構成するための較正方法であって、
   一組の複数マーカー（４ａ、４ｂ、４ｃ、４ｄ）を前記線源の場の中に配置するステップであって、前記複数のマーカー（４ａ、４ｂ、４ｃ、４ｄ）が同じ前記支持体（３）上に担持されている、当該ステップと、
   様々な撮像の際の前記線源の位置を決定すると共に、全ての複数マーカー（４ａ、４ｂ、４ｃ、４ｄ）及び撮像について、撮像から生じる像への各マーカーの投影の理論的位置（Ｐｍ）と像内でのその実際の位置（Ｐｄ）との間の距離の関数である総合誤差を最小にする変換パラメータの値を決定するステップと、
   撮像した像への前記複数マーカー（４ａ、４ｂ、４ｃ、４ｄ）の投影の位置を処理して、空間内での前記複数マーカー（４ａ、４ｂ、４ｃ、４ｄ）の位置を推測するステップと、
   を有し、
   撮像した像への前記複数のマーカー（４ａ、４ｂ、４ｃ、４ｄ）の投影の位置を処理して、前記複数のマーカー（４ａ、４ｂ、４ｃ、４ｄ）の想定した初期位置と変位した後に前記複数のマーカー（４ａ、４ｂ、４ｃ、４ｄ）が取る位置との間の前記複数のマーカーの剛性変位を特徴付ける変換パラメータを推測する、較正方法。
2. 前記最小にする処理が、全ての撮像及び前記複数のマーカー（４ａ、４ｂ、４ｃ、４ｄ）について変換パラメータ及び線源の位置に所与の値を割り当て、そして、これらの変換パラメータ値について大域誤差の最小化を決定し、このように決定された線源位置について総合誤差を最小にする変換パラメータ値を決定すること、を含んでいる、請求項１記載の較正方法。
3. 線源の位置の決定のため、１つを除く全ての線源に所与の値が割り当てられ、その線源については、全ての複数マーカー（４ａ、４ｂ、４ｃ、４ｄ）について、また変換パラメータに割り当てられた値について、また他の線源に割り当てられた位置について、大域誤差を最小にする位置が決定され、この決定が各々の線源について順に実行される、請求項２記載の較正方法。
4. 前記支持体が圧迫板である、請求項１記載の較正方法。
5. 前記剛性変位を特徴付ける変換パラメータが、３つの平行移動パラメータ及び３つの回転パラメータを含む請求項１記載の較正方法。
6. 検出手段（１）と、
   前記検出手段に対して移動可能である放射線源を構成する手段（２）と、
   撮影対象物を前記検出手段に向かって圧迫する支持体（３）と、
   前記放射線源の場の中に配置された一組の複数マーカー（４ａ、４ｂ、４ｃ、４ｄ）と、
   前記複数のマーカーを１つの同じ前記支持体上に担持する手段と、
   処理手段と、
   を有し、
   前記処理手段は、撮像した像への前記複数のマーカー（４ａ、４ｂ、４ｃ、４ｄ）の投影の位置を提供して、前記複数のマーカーの想定した初期位置と変位した後に前記複数のマーカーが取る位置との間の前記複数のマーカーの剛性変位を特徴付ける変換パラメータを推測し、
   前記処理手段は、様々な撮像の際の前記線源の位置を決定すると共に、全ての複数マーカー及び撮像について、撮像から生じる像への各マーカーの投影の理論的位置（Ｐｍ）と像内でのその実際の位置（Ｐｄ）との間の距離の関数である総合誤差を最小にする変換パラメータの値を決定する、
   撮像した像の再構成のための装置。
7. 前記処理手段は、総合誤差を最小にするため、全ての撮像及び前記複数のマーカー（４ａ、４ｂ、４ｃ、４ｄ）について変換パラメータ及び線源の位置に所与の値を割り当て、そして、これらの変換パラメータ値について大域誤差の最小化を決定し、このように決定された線源位置について総合誤差を最小にする変換パラメータ値を決定すること、を含んでいる、請求項６記載の装置。
8. 前記処理手段は、線源の位置の決定のため、１つを除く全ての線源に所与の値が割り当てられ、その線源については、全ての複数マーカー（４ａ、４ｂ、４ｃ、４ｄ）について、また変換パラメータに割り当てられた値について、また他の線源に割り当てられた位置について、大域誤差を最小にする位置が決定されるようにし、この決定が各々の線源について順に実行される、請求項７記載の装置。
9. 前記担持させる手段が圧迫板である、請求項６記載の装置。
10. 前記剛性変位を特徴付ける変換パラメータが、３つの平行移動パラメータ及び３つの回転パラメータを含む請求項６記載の装置。

Images