JP4196591B2 - レンダリングされたボリュームの画像に対する変換関数を非管理下で生成するための方法及びシステム - Google Patents
レンダリングされたボリュームの画像に対する変換関数を非管理下で生成するための方法及びシステム Download PDFInfo
- Publication number
- JP4196591B2 JP4196591B2 JP2002149914A JP2002149914A JP4196591B2 JP 4196591 B2 JP4196591 B2 JP 4196591B2 JP 2002149914 A JP2002149914 A JP 2002149914A JP 2002149914 A JP2002149914 A JP 2002149914A JP 4196591 B2 JP4196591 B2 JP 4196591B2
- Authority
- JP
- Japan
- Prior art keywords
- volume
- image processing
- dimensional data
- dimensional
- data
- 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
- 238000000034 method Methods 0.000 title claims description 56
- 230000006870 function Effects 0.000 title claims description 51
- 230000009466 transformation Effects 0.000 title claims description 28
- 238000009877 rendering Methods 0.000 claims description 38
- 238000006243 chemical reaction Methods 0.000 claims description 21
- 238000012545 processing Methods 0.000 claims description 20
- 238000003384 imaging method Methods 0.000 claims description 16
- 238000009607 mammography Methods 0.000 claims description 14
- 238000002591 computed tomography Methods 0.000 claims description 10
- 238000002604 ultrasonography Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 2
- 238000005070 sampling Methods 0.000 claims description 2
- 230000005855 radiation Effects 0.000 description 43
- 210000000481 breast Anatomy 0.000 description 34
- 230000033001 locomotion Effects 0.000 description 15
- 210000001519 tissue Anatomy 0.000 description 10
- 208000004434 Calcinosis Diseases 0.000 description 9
- 230000001133 acceleration Effects 0.000 description 7
- 206010006187 Breast cancer Diseases 0.000 description 6
- 208000026310 Breast neoplasm Diseases 0.000 description 6
- 238000009826 distribution Methods 0.000 description 6
- 230000005484 gravity Effects 0.000 description 6
- 210000004872 soft tissue Anatomy 0.000 description 5
- 238000005266 casting Methods 0.000 description 4
- 230000006835 compression Effects 0.000 description 4
- 238000007906 compression Methods 0.000 description 4
- 238000001514 detection method Methods 0.000 description 4
- 230000014509 gene expression Effects 0.000 description 4
- 238000013507 mapping Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000002059 diagnostic imaging Methods 0.000 description 3
- 239000012530 fluid Substances 0.000 description 3
- 230000003902 lesion Effects 0.000 description 3
- 206010028980 Neoplasm Diseases 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 201000011510 cancer Diseases 0.000 description 2
- 239000002131 composite material Substances 0.000 description 2
- 208000031513 cyst Diseases 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 230000002452 interceptive effect Effects 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 238000012552 review Methods 0.000 description 2
- 238000000844 transformation Methods 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 241000699670 Mus sp. Species 0.000 description 1
- DERRCVPMYGOTDB-UHFFFAOYSA-M P.[I-].I.I.I.[Cs+] Chemical compound P.[I-].I.I.I.[Cs+] DERRCVPMYGOTDB-UHFFFAOYSA-M 0.000 description 1
- 239000006096 absorbing agent Substances 0.000 description 1
- 238000007792 addition Methods 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- OFHCOWSQAMBJIW-AVJTYSNKSA-N alfacalcidol Chemical compound C1(/[C@@H]2CC[C@@H]([C@]2(CCC1)C)[C@H](C)CCCC(C)C)=C\C=C1\C[C@@H](O)C[C@H](O)C1=C OFHCOWSQAMBJIW-AVJTYSNKSA-N 0.000 description 1
- 229910021417 amorphous silicon Inorganic materials 0.000 description 1
- 230000003042 antagnostic effect Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000007599 discharging Methods 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 230000001678 irradiating effect Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000002156 mixing Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 238000007794 visualization technique Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/90—Dynamic range modification of images or parts thereof
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/008—Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T15/00—3D [Three Dimensional] image rendering
- G06T15/08—Volume rendering
-
- 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/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Graphics (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Generation (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Description
【発明の背景】
本発明は、全般的には、医用イメージングにおけるボリューム・レンダリングに関し、さらに詳細には、レンダリングしたボリュームの描出のための変換関数を容易に生成させる方法及びシステムに関する。
【0002】
マンモグラフィは、がんを初期段階で発見するのに望ましいような患者の乳房の一つまたは複数の画像を作成する低線量X線手技である。図1は従来技術によるマンモグラフィ装置10の一例を表している。マンモグラフィ装置10は一般に、支柱16に回動自在に装着されたアーム14に装着されているX線管12と、同じく支柱16に回動自在に装着されたアーム20に装着されているフィルムプレートまたはディジタル検出器18と、を含んでいる。X線管12及びアーム14と、ディジタル検出器18及びアーム20とは、拮抗平衡状態(counterbalanced)とし、X線管12とディジタル検出器18が手動により容易に回動でき、上下移動でき、さらに様々な角度方向で適所にロックできるようにしている。
【0003】
典型的なマンモグラフィ手技には概ね30分を要する。この手技には一般に、患者の各側の乳房からそれぞれ2枚の画像(上方向から1枚と横方向から1枚)を取得することを含む。例えば、図1に示すようにX線管12とディジタル検出器18を軸A(すなわち、頭尾方向(cranio−caudal))に沿った垂直向き配置に位置させた状態で患者の乳房のそれぞれの側に対して別々に画像を取得する。さらに、例えば、患者の一方の乳房では軸B1(すなわち、中外側斜位方向(medio−lateral oblique))に沿い、患者のもう一方の乳房では軸B2に沿うなど、X線管12とディジタル検出器18をある角度をつけた状態として、患者の乳房のそれぞれの側に対して別々に画像を取得している。
【0004】
この手技の中で、患者の乳房は、圧迫用パドル22(例えば、放射線に対して透明な一片のプラスチック)と、ディジタル検出器カバー18との間で圧迫させて乳房を平坦にし、これにより厚さを減少させると共に乳房組織を押し広げて、乳房が容易に画像化できるようにしている。上方向からと側面方向からのいずれかの方向で画像を取得する際に、X線管12はフィルムプレートまたはディジタル検出器18に対して概して直交、すなわち垂直に整列させている。次いで、担当医または放射線医はこの乳房画像(すなわち、マンモグラム)を精査して乳がんを発見している。
【0005】
上述の手技は乳がんの初期形態を発見する最善の方法の一つであるが、それでも、マンモグラムを精査する担当医や放射線医が乳がんを見落とすおそれがある。例えば、マンモグラム内の関心対象構造に重なっている放射線に対してデンス(高密度)で線維腺性の(fibroglandular)乳房組織によって不鮮明となることにより乳がんが見落とされることがある。
【0006】
乳がんの早期発見を改善させようとの試みから、静止した乳房及び静止したディジタル検出器に対してX線源を弧を描くように移動させながら複数の画像(すなわち、投影放射線像)を収集しているトモシンセシス乳房イメージング(Tomosynthesis breast imaging)が研究されてきた。複数の投影放射線像に対するシフト、スケール調整及び加算により、撮影対象乳房の検出器と平行な任意の面を再構成し、これにより重複している組織を関心対象構造から「除去(remove)」することができる。
【0007】
微小石灰化やマス(mass)、嚢胞、及びその他乳房の診断上重要な構造をトモシンセシス乳房イメージングから収集した一連の2次元(2D)面で描画することにより、重要な診断情報が提供される。しかし、データのボリュームは概して大きく、かつかなりの範囲にわたるデータ容量(data content)を含んでいる。したがって、放射線医は本質的に、例えば、概ね60〜80枚の画像(通常、圧迫した乳房厚の各1cmあたり10枚の画像であり、かつそれぞれが典型的には2304×1800のマトリックス・サイズを有する)を観察することにより、これらの面を概念的に3D構造に再構成させようと試みるが、こうした作業は困難でありかつ時間がかかる。
【0008】
ボリューム・レンダリングによれば、対象ボリュームの3次元(3D)描画を提供できる。ボリューム・レンダリングは、3D断層撮影データ(例えばコンピュータ断層(CT)、磁気共鳴(MR)、またごく最近では、3DディジタルX線マンモグラフィ(3DDM))を用いた患者のスクリーニングや診断のため医用イメージングにおいて特に有用であることが分かっている。ボリューム・レンダリング技法の使用を妨げる問題点の一つは、診断において有用となるボリューム・データの容量を示す初期変換関数の定義である。
【0009】
変換関数とは、個々のボクセルのデータ値を色相や不透明度にマッピングさせるためのメカニズムである。典型的には、この変換関数は3D断層撮影データに関する事前の知見に基づいて生成される。例えば、CTでは、ハウンスフィールドCT値を単位として計測される強度値の範囲は、様々な組織の種類ごとに概ね既知でありかつ再現性(reproducible)がある。こうした事前の知見を使用することにより、ユーザはレンダリングされたボリュームを描画させるのに必要な変換関数を決定しかつ操作することが可能となる。
【0010】
CTデータと異なり、MRにおける組織の強度値では、ハウンスフィールド値に対応する考え方がないため一般に再現性がない。一般に、MR画像における軟部組織は、特にデータを高超伝導マグネットにより収集する際にコントラストが高く、また信号対雑音比(SNR)も高くなる。したがって、変換関数は一般に、最小強度値と最大強度値の間の単純な直線傾斜を用いて3Dデータ組から設定することによりレンダリングされたボリュームを十分に描画することができる。しかし、コントラストが低いデータやSNRが低いデータ(例えば3DDM)では、適当な方式でデータをボリューム・レンダリングするのにこうした直線傾斜では容認できない場合がある。変換関数を手作業により操作しようとすることは、時間がかかり臨床上実現不可能であるため、典型的には実用的でない。
【0011】
したがって、低SNRでありかつ低コントラストである3Dデータをボリューム・レンダリングするのに適した変換関数を容易に決定するための方法及びシステムが必要である。
【0012】
【課題を解決するための手段】
対象ボリュームの3次元データのボリューム・レンダリングで使用するための変換関数を生成させるための方法及びシステムを提供する。
【0013】
第1の態様では、本方法は、対象ボリュームの3次元データを取得すること、この3次元データの複数のサンプルに関して選択した特性を評価すること、並びに、選択した特性に基づいて3次元データをボリューム・レンダリングするための変換関数範囲を計算すること、を含む。
【0014】
第2の態様では、本システムはイメージング装置及び画像処理装置を備える。この画像処理装置は、取得した画像データに基づいて対象ボリュームの少なくとも一つの3次元画像データ組を作成するように適合されると共に、さらにボリューム・レンダリングで使用するための変換関数範囲を計算するように適合されている。
【0015】
【発明の実施の形態】
本発明の特徴及び利点は、添付の図面と共に以下の本発明の詳細な説明を読めば明らかとなろう。
【0016】
がんの早期発見のため患者の乳房などの被検体102に対するトモシンセシス・データを作成するために使用されるX線投影放射線像を得るための装置100を図2及び3に図示している。この例示的な実施形態では、本装置は、放射線検出器及び患者の乳房に対する複数の位置を通過するように放射線源を滑らかに移動させ患者の乳房に関する複数の投影放射線像を取得するように重力と協同して動作しており、この放射線像は担当医または放射線医による検討を受けるように処理することができる。
【0017】
アクチュエータは、放射線源が複数の位置において実質的に停止するように放射線源の動きを制御するように動作可能であり、かつこれら複数の位置の間では放射線源が滑らかに移行するように動作可能である。この放射線源の加速度は、時間の経過と共に概して滑らかに(例えば、正弦波的に)変化することが望ましい。
【0018】
別法として、放射線源には第1の端部部分を有するアームを装着し、またその第2の端部部分は支柱に回動自在に装着して、放射線源が一つの弧に沿って移動可能にしている。このアームと動作可能に接続した油圧システムにより放射線源は重力の影響下で放射線検出器に対して少なくとも一つの高い位置から少なくとも一つのより低い位置まで制御可能に移動することができ、この際、少なくとも一つの高い位置と少なくとも一つのより低い位置の間にある複数の位置において放射線源から放出されて被検体を通過する放射線は放射線検出器により検出可能である。本明細書で使用する場合、「重力の影響下で(under the influence of gravity)」などの表現は、所望のトモシンセシス・データを得るためにX線源を移動させる際に、重力がX線源をある位置から順序が次の位置まで移行させるための起動力となることを意味している。さらに、本明細書で使用する場合、「動作可能に接続される(operably connected)」などの表現は、構成要素が互いに結合されてこれらの構成要素により所望の動作または動きが提供されることを意味しており、例えば、構成要素の部分間に着脱可能な結合を提供するような締結具(例えば、コッターピン結合、バネ解錠式ラッチ、その他)による機械的結合を含むが、これに限るものではない。
【0019】
被検体のトモシンセシス・データを取得するための方法は、放射線源と放射線検出器の間に配置した被検体に対して放射線源を高い位置に位置決めするステップと、放射線源が重力の影響下で高い位置からより低い位置まで移動することを可能にしながら被検体に複数の照射位置で放射線を照射するステップと、この複数の位置で被検体を通過した放射線を放射線検出器によって検出するステップと、を含む。本方法はさらに、被検体を通過した放射線を複数の位置の各々において検出することにより得た投影放射線像から、間隔をあけた複数の被検体透過平面画像を作成するステップを含んでいる。
【0020】
こうした装置は、例えば、ディジタル検出器を基準としてX線源を物理的に移動させ位置決めするために一つまたは複数のサーボモータを使用しているようなマンモグラフィ装置と比べてより簡単かつより堅牢である。こうしたモータ式システムでは、装置に対して推力性の(thrusting)動きや振動性の(jolting)動きが加えられ、手技の全体にわたって検出器及び患者の乳房を静止の(すなわち、固定した)位置に保つ機能に悪影響が及んでいる。本装置では装置全体にわたる急峻な推力性や振動性の動きを制限でき、手技の全体にわたって検出器及び患者の乳房をより適正に静止位置に保つことができるので有利である。
【0021】
例示的な装置100は、X線源110と、検出器120と、重力の影響下で検出器120を基準としてX線源110を制御可能に移動できるようにするためのアクチュエータ130と、を含んでいる。アクチュエータ130は、支柱131と、支柱131とX線源110とに装着した第1のアーム132(X線源アームともいう)と、支柱131と検出器120とに装着した第2のアーム134(X線検出器アームともいう)と、作動システム140と、を含んでいる。
【0022】
X線源110は典型的にはX線管であり、また検出器120は典型的にはディジタルX線検出器である。例えば、検出器120は、例えば、アモルファス・シリコン製トランジスタ−フォトダイオード・アレイ上にヨウ化セシウム蛍光体(シンチレータ)を有するような半導体放射線イメージャとすることがある。適当なその他の検出器としては、一つまたは複数の電荷結合素子(CCD)やX線を直接ディジタル信号に変換している直接型ディジタル検出器を含むことができる。検出器120は平坦でありかつ平坦な画像面を規定しているように図示しているが、例えば、湾曲した画像面を有する湾曲形状のディジタルX線検出器などディジタルX線検出器のその他の構成も適正に利用することができる。
【0023】
X線源110と検出器120とは互いに別々に移動可能である。例えば、アーム132はX線源110に装着された第1の端部211または部分と、支柱131に回動自在に装着された第2の端部212または部分とを有し、X線源110が一つの点または軸Rの周りで湾曲した両側矢印Cで示す方向に回転可能となるようにしている。同様に、アーム134は、検出器120に装着された第1の端部221または部分と、支柱131に回動自在に装着された第2の端部222または部分とを有し、検出器120が軸Rの周りで両側矢印Dで示す方向に回転可能となるようにしている。
【0024】
X線源110とアーム132は、アーム132の上側部分211(すなわち、軸Rの上方の部分)に対するモーメントが、アーム132の下側部分212(すなわち、軸Rの下方の部分)のモーメントと比べてより大きくなるような適当な重さとする。X線源110が第1の位置(すなわち、高い上側の位置201(図2及び3に示す))に位置決めされているとき、アーム132及びX線源110は、X線源110が第2の位置(すなわち、より低い位置202(図2及び3において想像線で示す))まで弧Cに沿って移動できるように軸Rの周りで回転しようとする傾向をもつことが望ましい。
【0025】
この例示的な実施形態では、作動システム140は、第1及び第2の油圧シリンダ(それぞれ、142及び144)を有する油圧システムと、リザーバ160と、ポンプ170と、バルブ172及び174と、を備えている。油圧シリンダ142は、支柱131に回動自在に装着されたシリンダ本体141と、X線源110を検出器120に対して位置決めするために上側端部の位置でアーム132の一方の側面に回動自在に装着可能としたピストン143と、を備えている。油圧シリンダ144は、支柱131に回動自在に装着されたシリンダ本体145と、X線源110を検出器120に対して位置決めするために上側端部の位置でアーム支柱132のもう一方の側面に回動自在に装着可能としたピストン147と、を備えている。
【0026】
リザーバ160、ポンプ170、並びにバルブ172及び174は、油圧シリンダ142及び144と動作可能に接続させて、流体をシリンダ本体141及び145内に導入させたり該本体内から排除したりしてアーム132の位置決め、したがって、検出器120に対するX線源110の高い位置(例えば、位置201)での位置決めをすることができる。油圧シリンダ142及び144がアーム132のそれぞれの側面に装着されていないとき、あるいはアーム132及びアーム134を支柱131に対して固定の位置に適当にロックさせることができるような場合である装置100が標準の頭尾方向または中外側斜位方向のマンモグラムを得るように動作しているときは、油圧シリンダ142及び144は概して垂直の配置位置に適当に安定化、すなわち保持されている。
【0027】
限定ではなく一例として、X線源、支柱及びアームは、典型的には、General Electric Medical SystemsのModel Senographe 2000D Digital Mammography Systemなどのタイプのものであり、これらは上述のように適当に修正されている。当業者であれば、油圧シリンダの代わりに空気圧式シリンダも適正に利用できることを理解するであろう。空気圧式シリンダでは、リザーバは不要となる。画像を取得するための上述した動きをもたらすために、装置100において別の様々な機械的構成要素(例えば、ボールねじドライブ、非油圧式システム、その他)を使用することもできることも理解されよう。
【0028】
図2を参照しながら、トモシンセシス・データを取得するための装置100に対する初期設定を以下に示す。X線検出器アーム134(したがって、検出器120)を、垂直方向(図面では基準線Yで示す)から固定の角度Bの位置に配置させる。中外側方向斜位像では、角度Bは典型的には、垂直方向から約50度から約70度までの範囲にあることが望ましい。角度Bを選択し終えた後、アーム134を動かないように適所に固定する。油圧シリンダ142はX線源アーム132に装着させる。ポンプ170に通電しシリンダ本体143を満たし、X線源110を図2に示すような高い位置201(例えば、検出器120に対して垂線方向に延びる軸L1から周知の技法に従った最大角度の位置)に上昇させて位置決めする。
【0029】
撮影対象ボリューム102(例えば、患者の一方の乳房)は、圧迫用パドル122と検出器120の間に位置決めし、複数の画像を取得している期間である高い位置からより低い位置202(想像線で示す)までX線源110を移動させる間は静止した固定位置に保持させる。この手技は、患者のもう一方の乳房の画像取得に関しても同様であるが、この場合には、油圧シリンダ144をアーム132に接続させ、アームの動きが図3に示すような対応するより上方の位置201とより低い位置202との間となるようにする。
【0030】
複数の画像の取得に関する一態様では、流体をシリンダ本体141からリザーバ160まで排出させるようにバルブ172を開き、X線源110が重力の影響下で高い位置201からより低い位置202まで滑らかに下降するようにすることができる。本明細書で使用する場合、「滑らかに下降する(smooth descent)」その他の表現は、X線源が、途中で停止せず、かつ加速度の急峻な変化を伴わない方式によりそれぞれの照射位置間を移動することを意味している。X線源110のこの下降の間で、X線源が撮影しようとする被検体と検出器に向けて放射線を放出することにより複数の画像(望ましくは8枚から15枚までの画像)が取得される。例えば、±24度の最大角度では、6度の間隔(図2に示す位置「a」から位置「i」まで)を使用することにより9枚の画像が得られる。線源が撮影しようとする被検体に向けて放射線を放出する各位置がそれぞれの照射位置となる。各位置において検出器により得る画像は、低い放射線量により収集することが望ましい。画像のすべてに対する全放射線量は、標準の単一像のマンモグラムの放射線量と比べて、同等であるか、あるいは医学的に問題とならない程度で高いだけとすることが望ましい。
【0031】
装置100は、アクチュエータ・システムを制御し、かつ検出器120に対するX線源110の複数の位置で取得した複数の画像を処理するためのプロセッサ150を含むことが望ましい(これについては以下でさらに説明する)。
【0032】
例えば、概ね±0.1度の感度を有する適当な傾斜計180をアーム132に装着しプロセッサ150と結合させることができる。さらに、プロセッサ150は典型的には、X線源110とインタフェースさせ、X線源を高い位置201からより低い位置202まで下降させながら適正な位置(例えば、位置「a〜i」)においてX線源110を付勢させるための信号を供給している。装置100は、アーム132を位置202で穏やかに停止させるためにバネや加圧したリザーバなどの適当な緩衝器を含むと有利である。また別法として、プロセッサ150は、X線源110がスキャンの終了時に位置202に接近するのに伴ってX線源110を減速させるようにバルブ172をゆっくりと閉じるようにプログラムされている。
【0033】
複数の画像の取得に関する別の態様では、X線源110は適正な位置の各々(例えば、位置「a〜i」)で静止させる、すなわち停止させるように移動させ、この時点で(もはや移動していない状態において)X線源を付勢させて照射を発生させることができる。例えば、プロセッサ150は典型的には、X線源110の下降を制御してX線源を所望の位置の各々において静止させるように、並びにこれら所望位置の各々にX線源を停止させている間(すなわち静止を維持している間)にX線源110を付勢させるための信号を提供するように傾斜計180及びバルブ172とインタフェースさせている。
【0034】
本発明から、当業者であれば、傾斜計ではなく直線エンコーダその他同様の装置を油圧シリンダに装着して角度計測値が得られるように較正することもできることを理解するであろう。同様に、このエンコーダは、コンピュータとインタフェースさせると共に、X線源の位置を制御するために使用することができる。
【0035】
図4〜6は、上側位置すなわち高い位置201から下側位置すなわち最終位置202までの滑らかな移行が得られることが望ましいようなX線源110の例示的な動きをグラフで表したものである。図4は、アームの停止位置(例えば、位置201と202)の間での距離または角度の変化を時間の経過に対して表したグラフである。図5及び6は、停止位置間での速度及び加速度のそれぞれの変化を時間の経過に対して表したグラフである。高い位置とより低い位置の間での時間経過に対するX線源の加速度及び減速度は、アームの動きに所望の滑らかさを提供するように本質的に概して滑らか、例えば正弦波的、とすることが望ましい。
【0036】
図7に示すように、複数の位置に関して時間の経過に対してそれぞれの位置「a〜i」の間の距離または角度を変化させることにより、X線源110が固定の位置、すなわち静止した位置に保持され、この時点でX線源を付勢できることが望ましい。照射間の各時間期間によりプロセッサはディジタル検出器からデータを読み出すことができる。X線源は、検出器の読み取りに要する最小時間の間で迅速かつ滑らかに移動させると有利である。
【0037】
アクチュエータの別の実施形態では、油圧シリンダのそれぞれからの排出ラインの中間に可変サイズのリザーバ対を使用し、このリザーバ対の充満と排出を交互に繰り返させることができる。各リザーバの体積は、X線源110が一つの位置から次の位置まで移動できるように油圧シリンダから排出させる、すなわち排除する流体の量を収容できるように構成させることが望ましい。例えば、一方のリザーバが満たされた時点で、X線源を付勢させて照射を発生させることができる。次いで、スキャンを続けるように油圧シリンダをもう一方のリザーバに接続し、この間に最初の方からこの下側リザーバ内に排出させることになる。
【0038】
この油圧システムとX線源の動きにより、その装置を臨床的に適用可能とするような十分な迅速さである概ね3〜5秒でトモシンセシス画像を収集することができる。典型的には、イメージング手技の間に得られた別々の投影放射線像は、オペレータが使用できるように検出器と結合させたプロセッサ150やその他の読み出し様式により処理している。
【0039】
装置100が(X線源110を位置「a〜i」に配置させて)作成した別々の複数の投影放射線像に対しては、前処理の後、検出器120の面と平行な被検体の一つまたは複数の面(すなわち、スライス)Sを再構成させるように、シフト、スケール調整、並びに加算(例えば、変換)を実施する。当業者であれば、所与のイメージング用途に適した様々な再構成技法から選択することになることを理解するであろう。本明細書で使用する場合、この複数の面またはスライスのことをトモシンセシス・データと呼ぶ。トモシンセシス・データは再構成させたデータを記述しており、一方、投影放射線像はX線曝射(x−ray shot)を示している。トモシンセシスにより、撮影された被検体内の任意のレベル(例えば、約1.0mm間隔)で、「焦点の合った(in−focus)」断層撮影面(すなわち、スライス)Sの再構成が可能となる。典型的には、スライス間隔が約1mmの等間隔のスライスが利用される。しかし、別の間隔及び/または不均等な間隔も代替的な実施形態となる。トモシンセシス・データ組の値は、撮影した被検体内の対応する構造の線形減衰値と明白な関数関係を有しているのが理想的である。本発明によれば、担当医または放射線医は正常な乳房組織からの「構造雑音」を通した観察により、乳がんの検出及び特徴付けを改善させることができるので有利である。
【0040】
限定ではなく一例として、断層撮影面(すなわち、スライス)Sの作成は、別法として、Niklasonらに対する米国特許第5,872,828号に開示された処理法、あるいは当技術分野で周知のその他適当な再構成アルゴリズムによっても取得することができる。
【0041】
上述した断層イメージング法を利用した後において、得られた断層撮影データのボリュームは一般に大きく、またそのデータの容量も一般にボクセル強度値のかなりのダイナミックレンジを包含している。本発明の実施形態では、かなりのデータ・ボリューム及びデータ容量を有するトモシンセシス・データの(例えば、診断を目的とした)描画を可能にする方法を提供する。データ・ボリュームに関しては、複数の間隔をあけた画像が、100マイクロメートルの平面画像分解能あるいはこれより良好な分解能(乳房撮影の場合)、並びに1000マイクロメートルの画像間隔分解能(image spacing resolution)あるいはこれより良好な分解能を有することが望ましい。典型的には、各画像は概ね2304×1800ボクセルのマトリックス・サイズを表している。本明細書で使用する場合、ボクセル(voxel)とは、X、Y及びZ軸のそれぞれに沿ったNx×Ny×Nzのデータサンプルからなる一つのデータ・ボリューム範囲内の単一のデータ・ポイントのことをいう。典型的には、圧迫した乳房の厚さを典型的には約4cm〜約7cmの範囲にした状態で圧迫厚さ1cmあたりに10画像が存在する。したがって、本明細書に記載した断層イメージングにより大きなデータ・ボリュームを提供できる。データ容量に関しては、上述した技法により取得されるデータのダイナミックレンジは、撮影した構造の減衰係数のレンジ差を2バイトの情報(ボクセル強度)を利用して記憶していることが望ましい。ボクセル強度値と撮影した構造の減衰係数の間にはある関係が存在する。本再構成技法の出力はこの関係を比例関係または反比例関係のいずれにもすることができる。反比例関係の場合、乳房内の構造の減衰係数が大きいほど、その強度値は小さくなる。例えば、微小石灰化は、繊維組織に関するボクセル値と比べてより小さいボクセル値を有する。ボクセル強度値と撮影した構造の減衰との間の関係は別の関数関係とすることも可能である。
【0042】
対象ボリューム(例えば、患者の乳房)の3次元画像を作成する一方法は、対象ボリュームのトモシンセシス・データを上述のようにして取得すること、並びにこのトモシンセシス・データに3D描画技法を適用することにより3次元画像を作成すること、を含む。本明細書で使用する場合、3D描画技法(3D visualization technique)とは、周知の技法に従って任意の観察角度で得られる、最大、最小、合成(composite)、表面、あるいは平均化(averaged)の投影など様々な周知のボリューム・レンダリング技法のことをいう。
【0043】
ボリューム・レンダリングとは、撮影対象ボリュームの3次元レンダリングを生成させるために様々な画像作成面上で投影画像を作成するための様々な技法のことをいう。上述の方法により取得したトモシンセシス・データ(すなわち、複数の間隔をあけた平面画像)は協同して、対象ボリュームを表す一つの3次元データ・ボリュームを形成する。図8を参照すると、0から(n−1)までの番号を付けたn枚のスライスを含む例示的なデータ・ボリューム40を表している。投影画像は番号36で示している。このデータ・ボリュームから3次元画像を作成するために、図2のプロセッサ150は、データ・ボリュームに対してボリューム・レンダリング技法を用いた変換を行うように適合させている。例示的な実施の一形態では、ボリューム・レンダリングは最大不透明度レイ・キャスティング(ray casting)により実施する。本明細書で使用する場合、不透明度(opacity)とは、透明度の一尺度と定義する。この実施形態では、レイ(ray)38がデータ・ボリューム40を通ってキャストされる。レイ38がデータ・ボリューム40の各面を通過するのに伴って、データ・ボリュームの当該スライスの各々からの画素とレイ38との交差部分から強度を計算し累積することによりレンダリングされた画像が形成される。プロセッサ150は周知の方法に従ったレイ・キャスティング式ボリューム・レンダリングを繰り返し実行する。一般に、2次元スライスの各方向での画素数(例えば、x方向の画素数及びy方向のスライス数)は分解能要件に基づいて事前に決定されている。典型的には、最終画像の分解能は事前に決定されている。しかし、処理時間を短縮するため対話式作業の間ではより低い分解能を利用することが多い。3D画像を作成するための方法に関する実施の一形態では、ユーザによる画像操作(例えば、様々な分解能での方向操作や観察角度操作)を可能とするようにユーザ入力に応答して選択した3D描画技法を適用している。別の実施形態では、プロセッサ150は、より高分解能の画像に対する対話式作業を可能にするように適合させている。代替的な実施形態では、トモシンセシス・データを用いて乳房の3次元画像を作成するために、2次元または3次元のテクスチャ・マッピング、最大強度値投影、最小強度値投影、合成レイ・キャスティングなど別のボリューム・レンダリング技法を選択している。また別の実施形態では、サーフェス・レンダリング技法を、単独で、あるいは様々なボリューム・レンダリング技法と組み合わせて使用することができる。サーフェス・レンダリング技法により、トモシンセシス・データを検討する際に望ましい機能である形状が規則的及び不規則な関心対象ボリュームを容易に抽出する機能が可能となることを理解されたい。換言すると、サーフェス・レンダリングは、何らかの幾何学物体(例えば、球形、楕円形、円筒形など)に対して使用されることになる。3Dデータがこれら幾何学物体と交差することにより、撮影した3Dボリュームによる関心対象ボリュームの再構成が可能となる。
【0044】
3次元レンダリングされた画像を診断目的で表示させかつ操作するために、レイ・キャスティング技法における各ボクセルに対して、色相と不透明度という二つの主要なパラメータを設定することが望ましい。典型的には、これらのパラメータの設定に使用できる二つの区分線形変換関数が存在する。変換関数の一つは色相(例えば、赤、緑及び青、すなわちRGB)に関する設定に使用され、もう一つの変換関数は不透明度に関する設定に使用される。一般に、レンダリングされた3D画像はほとんどの放射線医ではグレー表示としているため、赤、緑及び青の各値はすべて同じであるのが普通である。
【0045】
描画の分野では、RGBAカラー(赤、緑、青及びアルファ)をもつカラーモニタに言及するのが典型的である。これらの色相はモニタ内の電子銃に当てはまる。赤の電子銃だけをオンにすると画面上の画素は赤色で表示され、赤と緑の電子銃をオンにするとその画素は黄色で表示され、また、すべての電子銃をオンにするとその画素は白で表示される。これらの電子銃の強度の大きさにより様々なレベルの濃淡(shading)が生成される。すべての電子銃の強度を同じ量だけ変化させた場合(したがって、同じ値とした場合)、その画素は様々なグレイスケール値となる。典型的には、放射線医はグレイスケール値を有するフィルムを観察するように訓練されており、放射線医はディジタル画像であってもグレイスケール値で観察することをより好むが、放射線医がカラー画像で訓練を受けていれば状況は違ったものとなり得る。アルファ値(不透明度)は、二つの面を融合させること(blending)に関係する。一方の面は画面から生成される色相のことを指し、もう一方は不透明度面(すなわち、アルファ面)である。したがって、各画素は、その画素に関連する一つのアルファ値を有する。一般に、変換関数が単調増加の場合では、値1は画面上の画素の不透明度が極めて高いことを意味し、また値0はその画素が極めて透明であること意味する。
【0046】
各変換関数は、利用可能なディスプレイの深度(depth)に応じて、原データの全ダイナミックレンジを8ビットまたはこれを超えるビット数のいずれかにマッピングしている。このデータの全ダイナミックレンジは典型的には15ビットであるが、情報の深度を失わずに関心対象ボリューム(VOI)をレンダリングするのに要するダイナミックレンジでは、必要となるビット数は一般にこれより少なくて済む。例えば、微小石灰化と、マス(mass)、嚢胞その他の病変部や軟部組織構造との両方を全深度分解能で観察する場合と比べて、微小石灰化のみを全深度分解能で観察するのに必要なビット数はより少なくて済む。したがって、本発明の実施形態を利用して乳房構造を高分解能で描画することができる。
【0047】
微小石灰化の強度値分布は軟部組織構造の強度値分布と比べてより分散されていることに留意されたい。したがって、15ビット深度の分解能(すべての強度値が微小石灰化と軟部組織構造の両方に対する分布を含んでいる場合)を8ビット深度の分解能にマッピングまたはスケール変換しようとすると、情報の深度が失なわれることになる。しかし、一つの強度値分布(例えば、軟部組織構造)のカバー範囲は15ビット未満であり、したがって、一つの強度値分布を8ビット表示にスケール変換することは、その特定の分布に対しては情報の深度の損失を少なくすることにつながる。
【0048】
図9は二つの区分線形変換関数を表している。その一方は、乳房内の微小石灰化を伴う関心対象ボリュームからのボクセル値のマッピングに使用し、もう一方はスピキュラ状病変を伴うVOIのマッピングに使用する。この二つの変換関数に対応してボリューム・レンダリングした画像により、各構造(微小石灰化やスピキュラ状病変)の描画が可能となる。
【0049】
本発明の実施の一形態では、変換関数を生成するための方法は、対象ボリュームに対応する3次元(3D)データ組を取得すること、この3Dデータ組をサンプリングすること、そのサンプルの強度値に対する平均値と標準偏差を計測すること、並びに直線傾斜の変換関数の上限及び下限を計算すること、を含む。
【0050】
例示的な実施の一形態では、その3次元データ組は、トモシンセシス及びボリューム・レンダリングに関して上述した方法により収集している。したがって、この3Dデータ組は、サンプル数nでランダムにサンプリングし、次いでこれらのサンプルの強度値に対する平均値と標準偏差を計測している。その後、この平均値と標準偏差に基づいて直線傾斜の変換関数の下限及び上限を次式により決定している。
【0051】
下限=LL=μ−N×σ
上限=UL=μ+N×σ (式1)
上式において、μ、σ、Nはそれぞれ、平均値、標準偏差、及び標準偏差の数である。対応する不透明度の変換関数は次式となる。
【0052】
TF(x)=A (x≦LLの場合)
TF(x)=(B−A)x/(UL−LL)+C (LL<x<ULの場合)
TF(x)=B (x≧ULの場合) (式2)
TF(x)は変換関数、xはボクセル強度値、A及びBは不透明度値、またこれらの値は0と1の間の値をもつ。Cは不透明度y軸に対するy切片である。この線形変換関数は組織の強度値に応じて単調増加または単調減少とすることができる。例えば、3DDMデータでは、ボリューム・レンダリングに最大強度値投影(MIP)を使用する場合、下限に対しては1という値を割り当て、また上限に対しては0という値を割り当てている。別のボリューム・レンダリング技法(例えば合成ボリューム・レンダリング)を用いる場合には、1未満の値(例えば、0.5)を使用する。その理由は、密度がより大きい組織は密度がより小さい組織と比べてより小さい強度値を有することに基づく(例えば、微小石灰化の強度値は小さいが、放射線医にとって診断上極めて重要である)。
【0053】
MIPの3D描画では、色相変換関数は不透明度変換関数と同じ形状を有する場合や、単に1という一定の値を有する場合もあり得る。別のボリューム・レンダリング技法(例えば、合成ボリューム・レンダリング)では、色相変換関数が同じ直線傾斜特性を有することもあるが、最大可能な色相範囲を利用するためには、その最小値をゼロとしかつその最大値を1とする必要がある。
【0054】
そのボリュームから採取するサンプル数が多いほど変換関数の再現性はより高くなる。しかし、犠牲となるのは処理時間の増加である。したがって、再現性対速度のトレードオフが生ずる。
【0055】
3DDMデータでは、そのSNRが低い。したがって、(式1)において六つ分の標準偏差を用いて変換関数を適正に評価できるようにすることが望ましい。
【0056】
プロセッサ150は、典型的には、X線手技の間に検出器120から受け取ったイメージング信号を処理するようにプログラムされたディジタル計算機である。本発明の実施の一形態では、プロセッサ150はさらに、トモシンセシス・データに基づいて、上でより詳細に説明したようにして3D画像を作成するように適合されている。本明細書で使用する場合、「ように適合され(adaptedto)」、「構成された(configured)」などの表現は、与えられた入力信号に応答して出力を提供させる一連操作を実行するようにプログラムされているような、アナログ式やディジタル式の計算機、特定用途向け装置(例えば、特定用途向け集積回路(ASIC))などの電気的構成要素の動作機能のことを示している。例示的なプロセッサ150は、ハードドライブ、キーボードなどの入力デバイス、マウス、磁気記憶媒体(例えば、テープ・カートリッジやディスク)、光学記憶媒体(例えば、CD−ROM)、及びディスプレイやプリンタなどの出力装置を備えた適当なコンピュータ(例えば、Sun MicrosystemsのPENTIUM(商標)ベースのコンピュータからなるSPARC20ワークステーションなど)を含むことがある。上でより詳細に説明したように、プロセッサ150には、変換を実行するための適当なプログラムがインストール(すなわち、移植)されている。
【0057】
別の実施形態では、対象ボリュームの3次元画像を作成するためのシステムを提供する。本システムは、イメージング装置及び画像処理装置を備える。このイメージング装置は、上述したように、当該ボリュームに対して複数の照射位置において放射線を照射し、さらに、この複数の照射位置において放射線を検出することにより導き出される対象ボリュームを透過させた複数の間隔をあけた平面画像を作成している。この画像処理装置(例えば、図2のプロセッサ150)は、この間隔をあけた平面画像に基づいて対象ボリュームの3次元画像を作成している。画像処理装置のさらに別の実施形態では、その画像処理装置は上述のようにして変換関数を計算するように適合されている。この画像処理装置は、ユーザの介入に頼らずに取得データに基づいて変換関数を計算し、最適な変換関数を「オンザフライで」計算できるように適合させることが望ましい。本明細書で使用する場合、「オンザフライ(on the fly)」での計算とは、データ収集と実質的に同時事象であることを示す意図である。したがって、本発明の実施形態を利用したシステムにより、放射線医その他のユーザは、変換関数に関する追加の操作を実質的に行うことなくレンダリングされたボリュームを描画することができる。
【0058】
3D画像は、システムのユーザに対して表示させるために利用可能である(表示装置は図示していない)。別の実施形態では、その画像処理装置は、所望の表示形式に関するシステムのユーザからの入力に対して対話式に応答するように適合させている。例えば、放射線医などのユーザは、レンダリングされた3Dボリュームの特定の領域または観察角度を選択することができ、また別法として、強度などのパラメータに基づいた3D画像の再検討を希望することもできる。プロセッサ150は、適当なプログラムを用いて、レンダリングされた3D画像の所望の表示に関するユーザ入力に応答するように適合させていることを理解されたい。さらに別の実施形態では、プロセッサ150はレンダリングされた3D画像の動画映像表示を提供するように適合させている。
【0059】
この説明から、本装置は、乳房撮影に限定することなく様々な被検体を画像化するように構成できることを理解されたい。さらに、当業者であれば、この変換関数の計算は、ボリューム・レンダリングが利用されるような任意の3D用途にも適用可能であること、並びに様々な医用イメージング様式のうちの任意の様式(例えば、超音波)にも適用可能であることを理解するであろう。
【0060】
本発明の好ましい実施形態について図示すると共に本明細書に記載してきたが、こうした実施形態が例としてのみ提供されていることは明らかであろう。当業者であれば、ここに記載した本発明を逸脱することなく、多くの変形形態、変更形態及び置換形態を生じさせるであろう。したがって、本発明は本特許請求の範囲の精神及び趣旨によってのみ限定させようとする意図である。
【図面の簡単な説明】
【図1】従来技術によるX線マンモグラフィ装置の斜視図である。
【図2】本発明の実施形態が適用可能な、患者の乳房のトモシンセシス・データを取得するための装置の正面図である。
【図3】図2に示す装置を、患者のもう一方の乳房のトモシンセシス・データを取得するように構成させている正面図である。
【図4】断層撮影データを生成するために使用する投影放射線像を取得するための手技の間で、図2及び3に示すX線源の動きに対する距離または角度、速度、並びに加速度のそれぞれを時間の経過に対して図示した部分的グラフである。
【図5】断層撮影データを生成するために使用する投影放射線像を取得するための手技の間で、図2及び3に示すX線源の動きに対する距離または角度、速度、並びに加速度のそれぞれを時間の経過に対して図示した部分的グラフである。
【図6】断層撮影データを生成するために使用する投影放射線像を取得するための手技の間で、図2及び3に示すX線源の動きに対する距離または角度、速度、並びに加速度のそれぞれを時間の経過に対して図示した部分的グラフである。
【図7】図4と同様のグラフを複数の画像を取得するための時間期間に及ぶように図示したグラフである。
【図8】本発明の実施形態に従ってデータ・ボリュームを通過させてキャストしたレイを表した模式図である。
【図9】本発明の実施形態において有用な二つの区分線形変換関数を図示したグラフである。
【符号の説明】
10 従来技術によるマンモグラフィ装置
12 X線管
14 アーム
16 支柱
18 フィルムプレート、ディジタル検出器
20 アーム
22 圧迫用パドル
100 トモシンセシス・データを取得する装置
102 撮影対象ボリューム、患者の乳房
110 X線源
120 検出器
122 圧迫用パドル
130 アクチュエータ
131 支柱
132 X線源アーム、第1のアーム
134 X線検出器アーム、第2のアーム
140 作動システム
141 シリンダ本体
142 第1の油圧シリンダ
143 シリンダ本体
144 第2の油圧シリンダ
145 シリンダ本体
147 ピストン
150 プロセッサ
160 リザーバ
170 ポンプ
172 バルブ
174 バルブ
180 傾斜計
201 上側の位置
202 より低い位置
211 アーム132の上側部分
212 アーム132の下側部分
221 アーム134の第1の端部
222 アーム134の第2の端部
Claims (11)
- 対象ボリュームの3次元データのボリューム・レンダリングで使用するための変換関数を画像処理装置が生成させる方法であって、
対象ボリュームから取得した3次元データの複数のサンプルに関して選択した特性を前記画像処理装置が評価するステップと、
選択した特性に基づいて3次元データをボリューム・レンダリングするための変換関数範囲を前記画像処理装置が計算するステップと、を含み、
変換関数範囲を計算する前記ステップが、TF(x)を変換関数、xをボクセル強度値、A及びBを不透明度値、Cを不透明度y軸でのy切片とし、さらにLLを変換関数範囲の下限、ULを変換関数範囲の上限として、
TF(x)=A (x≦LLの場合)
TF(x)=(B−A)x/(UL−LL)+C (LL<x<ULの場合)
TF(x)=B (x≧ULの場合)
に従って実行されており、
前記下限及び上限が、μ、σ、Nをそれぞれ平均値、標準偏差、標準偏差の数として、
下限=LL=μ−N×σ
上限=UL=μ+N×σ
に従って計算されている、方法。 - 前記3次元データがノイズ及び低コントラストの影響を受けている、請求項1に記載の方法。
- 前記3次元データが、コンピュータ断層(CT)、磁気共鳴(MR)、超音波、及び3次元ディジタルX線マンモグラフィ(3DDM)のうちの少なくとも一つによって収集されている、請求項1に記載の方法。
- 前記選択した特性が平均値と標準偏差である、請求項1に記載の方法。
- 対象ボリュームの3次元データ組のボリューム・レンダリングで使用するための変換関数を画像処理装置が生成させる方法であって、
前記画像処理装置が対象ボリュームから取得した3次元データ組をサンプリングして対象ボリュームに対する複数のサンプルを生成させるステップと、
前記複数のサンプルに対する平均値と標準偏差を前記画像処理装置が計測するステップと、
サンプルに対する平均値及び標準偏差の前記計測に基づいて、3次元データをボリューム・レンダリングするために使用する変換関数の範囲を前記画像処理装置が生成させるステップと、を含み、
前記変換関数範囲が上限と下限を備えており、
前記上限をULとしかつ前記下限をLLとして、UL及びLLが、μ、σ、Nをそれぞれ平均値、標準偏差、標準偏差の数として、
下限=LL=μ−N×σ
上限=UL=μ+N×σ
に従って計算されている、方法。 - 前記3次元データ組がノイズ及び低コントラストの影響を受けている、請求項5に記載の方法。
- 前記3次元データ組がコンピュータ断層(CT)、磁気共鳴(MR)、及び3次元ディジタルX線マンモグラフィ(3DDM)のうちの少なくとも一つによって収集されている、請求項5に記載の方法。
- 対象ボリュームの3次元データのボリューム・レンダリングのためのシステムであって、対象ボリュームの画像データを取得するためのイメージング装置と、取得した画像データに基づいて対象ボリュームの少なくとも一つの3次元画像データ組を作成するように前記イメージング装置と結合されており、ボリューム・レンダリングで使用するための変換関数範囲を計算するように適合されている画像処理装置と、を備え、
前記変換関数範囲が上限と下限を備えており、
前記上限をULとしかつ前記下限をLLとして、UL及びLLが、μ、σ、Nをそれぞれ平均値、標準偏差、標準偏差の数として、
下限=LL=μ−N×σ
上限=UL=μ+N×σ
に従って計算されている、システム。 - 前記イメージング装置が、コンピュータ断層(CT)装置、磁気共鳴(MR)イメージング装置、超音波、及び3次元ディジタルX線マンモグラフィ(3DDM)装置のうちの少なくとも一つを含む、請求項8に記載のシステム。
- 前記画像処理装置がさらに、対象ボリュームの画像データの取得と実質的に同時に前記変換関数範囲を計算するように適合されている、請求項8に記載のシステム。
- 前記3次元データがノイズ及び低コントラストの影響を受けている、請求項8に記載のシステム。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US09/682157 | 2001-07-27 | ||
US09/682,157 US7085406B2 (en) | 2001-07-27 | 2001-07-27 | Method and system for unsupervised transfer function generation for images of rendered volumes |
Publications (3)
Publication Number | Publication Date |
---|---|
JP2003116825A JP2003116825A (ja) | 2003-04-22 |
JP2003116825A5 JP2003116825A5 (ja) | 2008-07-31 |
JP4196591B2 true JP4196591B2 (ja) | 2008-12-17 |
Family
ID=24738475
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2002149914A Expired - Fee Related JP4196591B2 (ja) | 2001-07-27 | 2002-05-24 | レンダリングされたボリュームの画像に対する変換関数を非管理下で生成するための方法及びシステム |
Country Status (3)
Country | Link |
---|---|
US (1) | US7085406B2 (ja) |
JP (1) | JP4196591B2 (ja) |
DE (1) | DE10223185A1 (ja) |
Families Citing this family (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE10301941B4 (de) * | 2003-01-20 | 2005-11-17 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. | Kamera und Verfahren zur optischen Aufnahme eines Schirms |
JP4709604B2 (ja) * | 2004-08-03 | 2011-06-22 | 株式会社東芝 | 画像表示装置、画像表示方法、記憶媒体及びプログラム |
US9024949B2 (en) * | 2004-10-13 | 2015-05-05 | Sony Corporation | Object representation using distance functions |
GB0523084D0 (en) * | 2005-11-11 | 2005-12-21 | Cancer Res Inst Royal | Imaging method and apparatus |
DE102006026490B4 (de) * | 2006-06-07 | 2010-03-18 | Siemens Ag | Radiotherapievorrichtung mit Angiographie-CT-Vorrichtung |
US7903113B2 (en) * | 2006-07-14 | 2011-03-08 | Siemens Medical Solutions Usa, Inc. | Systems and methods of image rendering from datasets |
JP2008068032A (ja) * | 2006-09-15 | 2008-03-27 | Toshiba Corp | 画像表示装置 |
JP4857070B2 (ja) * | 2006-10-11 | 2012-01-18 | キヤノン株式会社 | 乳房撮影用x線ct装置 |
US7990378B2 (en) * | 2007-05-07 | 2011-08-02 | General Electric Company | Methods and apparatus for volume rendering |
JP5288737B2 (ja) * | 2007-07-20 | 2013-09-11 | キヤノン株式会社 | X線画像撮影装置 |
US7920669B2 (en) * | 2007-07-25 | 2011-04-05 | Siemens Aktiengesellschaft | Methods, apparatuses and computer readable mediums for generating images based on multi-energy computed tomography data |
DE102011076929A1 (de) * | 2011-06-03 | 2012-12-06 | Siemens Ag | Verfahren und Vorrichtung zur Darstellung von Volumendaten für eine Untersuchung von Dichteeigenschaften |
DE102011076930A1 (de) * | 2011-06-03 | 2012-12-06 | Siemens Aktiengesellschaft | Verfahren und Vorrichtung zur Anpassung der Darstellung von Volumendaten eines Objektes |
JP5439453B2 (ja) * | 2011-10-20 | 2014-03-12 | 株式会社東芝 | 画像表示装置 |
KR20150016032A (ko) * | 2013-08-02 | 2015-02-11 | 삼성전자주식회사 | 영상 복원 모드 선택이 가능한 영상 복원 방법 및 그 장치 |
JP6438408B2 (ja) * | 2013-10-24 | 2018-12-12 | キヤノン株式会社 | 情報処理装置、情報処理方法、制御装置、制御システム、制御方法、トモシンセシス撮影装置、x線撮影装置、画像処理装置、画像処理システム、画像処理方法、及び、コンピュータプログラム |
US10565774B2 (en) * | 2015-09-03 | 2020-02-18 | Siemens Healthcare Gmbh | Visualization of surface-volume hybrid models in medical imaging |
US10726608B2 (en) * | 2016-11-23 | 2020-07-28 | 3D Systems, Inc. | System and method for real-time rendering of complex data |
JP6759172B2 (ja) * | 2017-09-19 | 2020-09-23 | キヤノン株式会社 | 乳房撮影装置、乳房画像処理装置及び乳房画像処理方法 |
US10593099B2 (en) | 2017-11-14 | 2020-03-17 | Siemens Healthcare Gmbh | Transfer function determination in medical imaging |
US11030742B2 (en) * | 2019-03-29 | 2021-06-08 | GE Precision Healthcare LLC | Systems and methods to facilitate review of liver tumor cases |
Family Cites Families (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS6058504B2 (ja) | 1980-11-17 | 1985-12-20 | 株式会社東芝 | 再構成装置 |
US5146414A (en) * | 1990-04-18 | 1992-09-08 | Interflo Medical, Inc. | Method and apparatus for continuously measuring volumetric flow |
US5920319A (en) * | 1994-10-27 | 1999-07-06 | Wake Forest University | Automatic analysis in virtual endoscopy |
JP3704652B2 (ja) * | 1995-09-08 | 2005-10-12 | 株式会社日立メディコ | 3次元画像処理方法 |
EP1006880A4 (en) * | 1996-06-03 | 2005-04-06 | Mallinckrodt Inc | METHOD FOR CALIBRATING ULTRASONIC EQUIPMENT AND CONTRASTING EQUIPMENT AND SYSTEM THEREFOR |
DE69739995D1 (de) * | 1996-07-23 | 2010-10-28 | Gen Hospital Corp | Anordnung zur mammographie mittels tomosynthese |
JP3522488B2 (ja) * | 1997-04-10 | 2004-04-26 | アロカ株式会社 | 超音波診断装置 |
US6008813A (en) | 1997-08-01 | 1999-12-28 | Mitsubishi Electric Information Technology Center America, Inc. (Ita) | Real-time PC based volume rendering system |
US5934288A (en) | 1998-04-23 | 1999-08-10 | General Electric Company | Method and apparatus for displaying 3D ultrasound data using three modes of operation |
US6142942A (en) * | 1999-03-22 | 2000-11-07 | Agilent Technologies, Inc. | Ultrasound imaging system and method employing an adaptive filter |
US6102861A (en) | 1999-04-23 | 2000-08-15 | General Electric Company | Method and apparatus for three-dimensional ultrasound imaging using surface-enhanced volume rendering |
FI991036A (fi) | 1999-05-06 | 2000-11-07 | Planmeca Oy | Tomosynteesimenetelmä ja laite tomosynteesikuvantamisessa |
US6842502B2 (en) * | 2000-02-18 | 2005-01-11 | Dilliam Beaumont Hospital | Cone beam computed tomography with a flat panel imager |
US20020172409A1 (en) * | 2001-05-18 | 2002-11-21 | Motoaki Saito | Displaying three-dimensional medical images |
JP3730143B2 (ja) * | 2001-07-17 | 2005-12-21 | アロカ株式会社 | 超音波画像処理装置 |
-
2001
- 2001-07-27 US US09/682,157 patent/US7085406B2/en not_active Expired - Lifetime
-
2002
- 2002-05-24 DE DE10223185A patent/DE10223185A1/de not_active Ceased
- 2002-05-24 JP JP2002149914A patent/JP4196591B2/ja not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
DE10223185A1 (de) | 2003-02-06 |
US20050074155A1 (en) | 2005-04-07 |
US7085406B2 (en) | 2006-08-01 |
JP2003116825A (ja) | 2003-04-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4424579B2 (ja) | マンモグラフィ画像の高分解能3d描画のための方法及びシステム | |
JP4196591B2 (ja) | レンダリングされたボリュームの画像に対する変換関数を非管理下で生成するための方法及びシステム | |
Wu et al. | Tomographic mammography using a limited number of low‐dose cone‐beam projection images | |
JP5400326B2 (ja) | トモシンセシス画像を表示するための方法 | |
US5872828A (en) | Tomosynthesis system for breast imaging | |
AU767518B2 (en) | Apparatus and method for cone beam volume computed tomography breast imaging | |
JP3984822B2 (ja) | Ct潅流法のためのカラー・パラメトリック合成マップ | |
US6487432B2 (en) | Method and system for selecting and displaying medical image data | |
US7142633B2 (en) | Enhanced X-ray imaging system and method | |
JP4965433B2 (ja) | トランケートされた投影と事前に取得した3次元ct画像を用いる円錐ビームct装置 | |
CN104939850B (zh) | 成像断层合成系统、特别是乳腺造影系统 | |
EP2205157B1 (en) | System for quantification of neovasculature in ct volumes | |
US8041087B2 (en) | Radiographic imaging display apparatus and method | |
JP2015518392A (ja) | 光子計数撮像装置で投影データをビニングして画像を生成する方法とシステム{projectiondatabinningandimagegenerationinphotoncountingimagingmodality} | |
CN1487478A (zh) | Ct中调节清晰度和噪声的图像回溯滤波及窗口控制滤波 | |
KR101971625B1 (ko) | Ct 영상을 처리하는 장치 및 방법 | |
CN1513419A (zh) | 利用多个探测器面板的容积测定ct系统和方法 | |
KR101798083B1 (ko) | 단층 영상 처리 장치, 방법, 및 그 방법에 관련된 기록매체 | |
KR20180041862A (ko) | 관심 영역의 동적 설정이 가능한 2차원 엑스선 디텍터와 이를 포함하는 콘빔 ct 장치 및 그 동작 방법 | |
US20220015710A1 (en) | Systems and methods for patient positioning for imaging acquisition | |
CN100583161C (zh) | 用于显示以立体数据组成像的对象的方法 | |
Fahrig et al. | Fast tomosynthesis for lung cancer detection using the SBDX geometry | |
US11162909B2 (en) | System and method for colorizing a radiograph from cabinet X-ray systems | |
US10937161B2 (en) | System and method for colorizing a radiograph from cabinet X-ray systems |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20050520 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20080616 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20080701 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20080723 |
|
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: 20080826 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20080922 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111010 Year of fee payment: 3 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 4196591 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121010 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: 20131010 Year of fee payment: 5 |
|
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 |