(57)【要約】【課題】従来に比して解剖学的構造を好適に視覚化することができる医用画像処理装置及び医用画像処理プログラムを提供すること。【解決手段】実施形態に係る医用画像処理装置は、曲面決定部と、画像生成部とを具備する。曲面決定部は、複数のサブ構造を有する解剖学的構造に関する三次元データに対し、複数のサブ構造を通過する第１の曲面を決定する。画像生成部は、第１の曲面上の複数のサンプル点についての三次元データを用いて、第1の曲面に対応する二次元画像を生成する。【選択図】図２
2017164496.tif 000002
【特許請求の範囲】
【請求項１】
複数のサブ構造を有する解剖学的構造に関する三次元データに対し、前記複数のサブ構造を通過する第１の曲面を決定する曲面決定部と、
前記第１の曲面上の複数のサンプル点についての前記三次元データを用いて、前記第1の曲面に対応する二次元画像又は三次元画像を生成する画像生成部と、
を具備する医用画像処理装置。
【請求項２】
前記画像生成部は、レイキャスティング処理における複数のレイと前記第１の曲面との交点を前記複数のサンプル点として設定する請求項１記載の医用画像処理装置。
【請求項３】
前記画像生成部は、前記複数のレイの少なくとも一つに対して、当該少なくとも一つのレイと前記第１の曲面との交点における前記三次元データの値から、前記二次元画像上の画素値又は前記三次元画像上のボクセル値を決定する請求項２記載の医用画像処理装置。
【請求項４】
前記画像生成部は、前記解剖学的構造の少なくとも一部と交わる曲面を決定する処理、
前記解剖学的構造の内側に少なくとも部分的に落ち込む曲面の決定する処理、
前記解剖学的構造によって明確にされた領域を取り囲む曲面を決定する処理、
複数のデータサンプリングパスに沿った前記三次元データの値の変化に基づいて曲面を決定する処理、
のいずれかを少なくとも実行することで、前記複数のサブ構造を通過する前記第１の曲面を決定する、
請求項１乃至３のうちいずれか一項記載の医用画像処理装置。
【請求項５】
前記画像生成部は、前記解剖学的構造の厚さにわたる前記三次元データの値の予想プロフィールに基づいて、前記複数のサブ構造を通過する前記第１の曲面を決定する請求項１乃至４のうちいずれか一項記載の医用画像処理装置。
【請求項６】
前記画像生成部は、前記三次元データの各位置における値に基づいて前記第１の曲面の少なくとも一部を決定する処理、補間処理または補外処理によって前記第１の曲面の他の部分を決定する処理を実行することで、前記複数のサブ構造を通過する前記第１の曲面を決定する請求項１乃至５のうちいずれか一項記載の医用画像処理装置。
【請求項７】
前記曲面決定部は、前記複数のサブ構造を通過し境界及び位置の少なくとも一方が前記第１の曲面と異なる第２の曲面を決定し、
前記画像生成部は、前記第２の曲面上の複数のサンプル点についての前記三次元データを用いて、前記第２の曲面を基準とする前記二次元画像又は前記三次元画像を生成する請求項１乃至６のうちいずれか一項記載の医用画像処理装置。
【請求項８】
前記第２の曲面の前記境界及び前記位置の少なくとも一方を入力する入力部をさらに具備する請求項７記載の医用画像処理装置。
【請求項９】
前記第２の曲面は、前記第１の曲面を拡大または縮小した曲面である請求項７又は８記載の医用画像処理装置。
【請求項１０】
前記曲面決定部は、前記第１の曲面をスラブ曲面として決定し、
前記画像生成部は、前記スラブ曲面の厚みにわたる投影処理を実行することにより、前記二次元画像又は前記三次元画像を生成する請求項１乃至９のうちいずれか一項記載の医用画像処理装置。
【請求項１１】
前記画像生成部は、軸周辺の複数の面のそれぞれに対し前記レイキャスティング処理を実行する請求項２記載の医用画像処理装置。
【請求項１２】
前記複数の面は、前記軸に対して傾けられて設定される請求項１１記載の医用画像処理装置。
【請求項１３】
前記画像生成部は、前記解剖学的構造の展開図であるカーブドＭＰＲ画像又は前記解剖学的構造の投影画像として、前記二次元画像又は前記三次元画像を生成する請求項１乃至１２のうちいずれか一項記載の医用画像処理装置。
【請求項１４】
前記解剖学的構造は、少なくとも一本の肋骨、胸郭の少なくとも一部、頭がい骨の少なくとも一部、少なくとも一本の歯、少なくとも一つの顎、臓器の少なくとも一部、骨盤の少なくとも一部のうち、少なくともいずれかを含む請求項１乃至１３のうちいずれか一項記載の医用画像処理装置。
【請求項１５】
前記画像生成部は、最大値投影処理、最小値投影処理、ボリュームレンダリング処理を用いて前記二次元画像又は前記三次元画像を生成する請求項１乃至１４のうちいずれか一項記載の医用画像処理装置。
【請求項１６】
コンピュータに、
複数のサブ構造を有する解剖学的構造に関する三次元データに対し、前記複数のサブ構造を通過する第１の曲面を決定させる曲面決定機能と、
前記第１の曲面上の複数のサンプル点についての前記三次元データを用いて、前記第1の曲面に対応する二次元画像を生成させる画像生成機能と、
を実現させる医用画像処理プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本実施形態は、例えば、解剖学的構造を視覚化することができる医用画像処理装置及び医用画像処理プログラムに関する。
【背景技術】
【０００２】
臨床医が胸部検査を行う際に、変位や変位のない、骨折または骨病変を識別し数えられることが必要である。例えば、全ての肋骨が検査されていることを確実にし、その中の肋骨のどこに骨折や損傷が位置するのかを決定する必要がある。
【０００３】
例えばコンピュータ断層撮影法（ＣＴ）、磁気共鳴映像法（ＭＲＩ）や陽電子放射断層撮影法（ＰＥＴ）など、幅広い医用撮像モダリティは、診断またはその他の目的のための患者またはその他の被検体を描出する三次元医用撮像データを取得するための標準的な技術となった。
【０００４】
患者やその他の被検体の肋骨を観察して何か骨折があるかを特定するために、ＣＴデータまたはその他の医用撮像モダリティから生成された、体軸断面や斜断面での多断面再構成（ＭＰＲ）ビューを使うことが知られている。体軸または斜断面ＭＰＲビューは、例えば患者の胴体にわたるスライスなど、三次元画像データのセットを通して、体軸または斜断面スライスを具備していてもよい。臨床医は、順に体軸スライスまたは斜断面スライスのシリーズの段階を踏むことで、患者の肋骨を観察することが出来る。それぞれのスライスは、少なくとも一本またはそれ以上の肋骨を含むかもしれない。異なる肋骨が異なるスライス上に示されるかもしれない。また個々の肋骨が多数のスライス上に現れるかもしれない。
【０００５】
実際に、このような技術を使うことで、全ての肋骨が適切に検査されてきたことを臨床医に対して保証するのは困難になる可能性がある。さらに、一旦骨折が確認されても、一体特定の肋骨のどこにその骨折があるのかを突き止めることは、困難な可能性がある。それぞれのスライスは、特定の肋骨または複数の肋骨しか表さないため、どの肋骨が所定のスライス上に示されているかを見分けるのは、難しい可能性がある。また、臨床医にとって、体軸または斜断面ＭＰＲビューのセットを通して段階を踏むことは、時間のかかることかもしれない。
【０００６】
従来から、それぞれの肋骨についての曲面再構成（Ｃｕｒｖｅｄ  ＭＰＲ：「ＣＰＲ」、「カーブドＭＰＲ」とも言う。）画像を明確にする肋骨中心線追跡アルゴリズム等の技術が提案されている。合成画像はその後、肋骨のＣＰＲ画像集と背骨のＣＰＲ画像とから生成される。効率的に、それぞれの肋骨は画像において真っ直ぐにされた形で示されている。原理上、このような技術は、変位のない骨折をより簡単に可視化する。
【発明の概要】
【発明が解決しようとする課題】
【０００７】
しかし、従来技術による画像における肋骨の直線化は、肋骨の相対的な位置や方位の把握を困難なものとしてしまう可能性がある。肋骨の直線化は、画像の見た目を不自然なものにしてしまうかもしれない。さらに、合成画像が生成されてきたいくつかの方法において、生成された合成画像に、背骨と肋骨との間に断絶のようなものが現れるかもしれない。
【０００８】
骨折を特定するために、展開された、胸郭の最大輝度投影（ＭＩＰ）画像を生成することも提案されてきた。この様な展開されたＭＩＰ画像は、原則として、胸郭の優れた全体ビューを与えることが出来る。しかし、ＭＩＰなど、輝度投影技術を使い生成されたビューは、変位のない骨折、特に変位のない骨折が部分的な骨折の場合、巧みに隠してしまうことがある。
【０００９】
ＭＩＰは、例えばレイに沿った最高輝度を返すことで、最大輝度の点を選択する。骨折した骨、特に骨折が部分的な骨折に対し、輝度の最大値は、骨折してない点で起きるかもしれない。従ってある状況では、ＭＩＰは骨折の線に沿った点を描かないかもしれず、結果画像において骨折の可視を困難にするまたは不可視にしてしまう。
【００１０】
変位のない骨折は、例えば、ＭＩＰやその他の輝度投影技術など、肋骨の表面から取得された画像データ値により描出している、または画像データ値により支配された画像を生成する、撮像技術には見えない可能性があるので、特に特定が困難になることがある。
【００１１】
本実施形態は上記事情に着目してなされたもので、従来に比して解剖学的構造を好適に視覚化することができる医用画像処理装置及び医用画像処理プログラムを提供しようとするものである。
【課題を解決するための手段】
【００１２】
一実施形態に係る医用画像処理装置は、曲面決定部と、画像生成部とを具備する。曲面決定部は、複数のサブ構造を有する解剖学的構造に関する三次元データに対し、複数のサブ構造を通過する第１の曲面を決定する。画像生成部は、第１の曲面上の複数のサンプル点についての三次元データを用いて、第1の曲面に対応する二次元画像を生成する。
【図面の簡単な説明】
【００１３】
【図１】図１は、実施形態に係る医用撮像データ処理システムの概略図である。
【図２】図２は、図１のシステムの作動モード全体を描いているフローチャートである。
【図３】図３は、軸とデータサンプリングパスを描いている軸の周辺の複数の面と複数の線とにオーバーレイした被検体の胸郭の概要図である。
【図４】図４は、複数の面の一部を描いた略図であり、軸に対して垂直方向に対して傾けられた面を示している。
【図５ａ】図５ａは、軸から放射状に延びており、撮像データセットに示された解剖学的構造と交わる、データサンプリングパスを示す略図である。
【図５ｂ】図５ｂは、解剖学的構造を通過するデータサンプリングパスを示す図５ａの一部分を拡大した図であり、データサンプリングパスと解剖学的構造との交点の範囲のインジケーターを含んでいる。
【図６】図６は、多様体と、軸と、軸からレイパスが続き多様体と交差するレイと、を通過するスライスは、多様体が生成された画像データを通過するスライスの描出とオーバーレイされたスライスの概略図。
【図７】図７は、図１のシステムの作動モード概要を描いているフローチャート。
【発明を実施するための形態】
【００１４】
実施形態は、これより付随する図面を参照に例に限定されない方法で、説明される。
【００１５】
実施形態に係る医用画像処理装置１０は、図１に示されている。医用画像処理装置１０は、この場合はパーソナルコンピュータ（ＰＣ）またはワークステーションである、計算装置１２を具備しており、ＣＴスキャナ１４、一つ以上の表示スクリーン１６またはその他の表示デバイス、そしてコンピュータキーボード、マウスまたはトラックボールなど、入力装置または複数の入力装置１８に接続されている。
【００１６】
ＣＴスキャナ１４は、患者またはその他の被検体の領域のボリューメトリックＣＴデータを取得するように構成されている、任意のＣＴスキャナでもよい。患者またはその他の被検体の領域は、関心の少なくとも一つの解剖学的構造を具備している可能性がある。本実施形態において、患者の領域は胴体で、関心の解剖学的構造は肋骨である。その他の実施形態において、解剖学的構造は任意の適切な構造であってもよく、例えば頭がい骨または骨盤などの骨や、肝臓、心臓または胃などの臓器である。患者の領域のボリューメトリックＣＴデータは、画像処理に最適な一つ以上のボリューメトリック撮像データセット中へと再構成される可能性がある。
【００１７】
代替的な実施形態において、ＣＴスキャナ１４は、例えばＭＲＩ（磁気共鳴撮像）スキャナ、Ｘ線スキャナ、ＰＥＴ（陽電子放射断層撮影法）スキャナ、ＳＰＥＣＴ（単光子放出コンピュータ断層撮影）スキャナ、または超音波スキャナなど、三次元医用撮像データを提供するように構成された撮像モダリティであれば、任意のボリューメトリックスキャナと置き換えられてもよいし、または補われてもよい。
【００１８】
本実施形態では、ＣＴスキャナ１４によって取得されたボリューメトリックＣＴデータセットは、メモリ２０に格納され、その後に計算装置１２に提供される。代替的な実施形態において、ボリューメトリックＣＴデータセットは、画像保管通信システム（ＰＡＣＳ）の一部を形成する可能性のある、遠隔データストア（図示せず）から供給される。メモリ２０または遠隔データストアは、メモリストレージの任意の適切な形態を具備していてもよい。
【００１９】
計算装置１２は、撮像データセットを自動的にまたは半自動的に処理するための処理リソースを提供し、中央処理ユニット（ＣＰＵ）を具備する。本開示において、計算装置１２は、多様体決定回路２４と、画像生成回路２８とを含む。
【００２０】
本実施形態において、多様体決定回路２４と画像生成回路２８は、実施形態の方法を実行することが可能なコンピュータ可読指示を有するコンピュータプログラムの手段により、計算装置１２でそれぞれ実行される。例えば、多様体決定回路２４と画像生成回路２８は、例えばＣＰＵ２２によって、計算装置１２で実行可能な個別のコンピュータプログラムまたはアルゴリズムとして、それぞれ実行される可能性がある。しかし、その他の実施形態において、回路は、一つ以上のＡＳＩＣｓ（特定用途向け集積回路）またはＦＰＧＡｓ（フィールド・プログラマブル・ゲート・アレイ）として、実行されてもよい。
【００２１】
計算装置１２は、ハードドライブや、ＲＡＭ、ＲＯＭ、データバス、様々なデバイスドライバを含むオペレーティング・システムや、グラフィックス・カードを含むハードウェアデバイス、も包括するＰＣのその他構成要素も持つ。このような構成要素は、明確にするために図１には図示されていない。
【００２２】
図１のシステムの工程は、図２のフローチャートにおける概略で示されている。プロセスの第一ステージ１００で、計算装置１２は医用撮像データセットを受信する。この実施形態において、医用撮像データセットは、スキャンボリュームにおける対応位置でのＣＴスキャンにより決定された、Ｘ線吸収レベルを表すそれぞれのボクセルの輝度値を伴う、ボクセルのセットを具備する三次元ＣＴデータを具備する。
【００２３】
代替的な実施形態において、医用撮像データセットを取得するために、任意のその他適切なモダリティが使われてもよいが、例えば磁気共鳴（ＭＲ）撮像または陽電子放出断層撮影（ＰＥＴ）撮像がある。
【００２４】
プロセスの次の段階１０２で、多様体決定回路２４は、多様体決定プロセスにおいて使用するための軸を決定する。軸２００は、図３および４に概略的に示されている。本実施形態において、軸２００は、医用撮像データセットの被検体であり、背骨と実質的に位置合わせされた、患者の上下（ｓｕｐｅｒｉｏｒ−ｉｎｆｅｒｉｏｒ）軸に従う。この場合、軸２００は、ＣＴスキャナの長軸および患者の上下軸とは実質的に一直線であるという理解に基づき、医用撮像データセットの座標系から設定される。
【００２５】
その他の実施形態において、任意の適切なプロセスが、軸を決定するために使われてもよい。例えば、セグメンテーションプロセスは、背骨またはその他の解剖学的構造の位置または方位を決定するために医用撮像データセットに実行されてもよいし、セグメンテーションに基づいて、軸２００の位置または方位が決定されてもよい。
【００２６】
プロセスの次の段階１０４で、多様体決定回路２４は、軸２００の周りに拡がっている複数の面２１０を決定する。軸２００と垂直方向に延長するよりも、軸２００に対して面２１０が傾けられるのが、本実施形態の特徴である。一般に、ヒト患者のそれぞれの個々の肋骨は、背骨の方向に対し傾けられた面と位置合わせされており、面２１０を傾けることで、面２１０は肋骨とより一直線に並べられるだろう。その他の実施形態において、面２１０は、任意の解剖学的構造と実質的に位置合わせされて傾けられてもよい。
【００２７】
本実施形態において、面２１０の全ては、図４に角度アルファで示されているように、同じ角度で傾いている。その他の実施形態において、異なる面２１０は、異なるチルト角を有していても良い。実施形態において、面２１０のチルト角が固定されてもよく、所定のチルト角、またはチルト角が自動的か半自動的に決定されてよく、例えば、医用撮像データセットのセグメンテーションに基づく、および／またはユーザ入力に従って傾き角度が決められてもよい。
【００２８】
プロセスのステージ１０６で、多様体決定回路２４は、それぞれの面２１０上の複数のデータサンプリングパス２２０を決定する。データサンプリングパスの例のいくつかが図３に示されている。データサンプリングパス２２０は、軸２００に対して放射状に拡がっている。
【００２９】
本実施形態において、それぞれの面２１０に対する複数のデータサンプリングパス２２０は、軸２００の周りを実質３６０度にわたって分配されている。その他の実施形態において、データサンプリングパス２２０は、３６０度よりも小さい角度にわたって分配されているかもしれない。データサンプリングパス２２０は、面の周辺を実質的に等しい角度間隔で空いている可能性がある。
【００３０】
図５ａは、軸２００が面２１０の中心にある、単一面２１０の概略図である。図５ａは、面２１０に対応する医用撮像データセットを通したスライスの描出とオーバーレイされている。図５ａは、軸２００から放射状に延びる複数のデータサンプリングパス２２０を示している。明確にするために、図５ａにはプロセスのステージ１０６において決定される、データサンプリングパス２２０のサブセットのみが描かれている。さらに面２１０上のデータサンプリングパス２２０は、存在しているものの図５ａには図示されていない。
【００３１】
図５ｂは、図５ａの一部分の拡大図で、二つのデータサンプリングパス２２０ａおよび２２０ｂの一部を示している。それぞれのデータサンプリングパス２２０ａおよび２２０ｂは、医用画像データセットにおいて骨の領域２３０と交差している（この場合、骨は肋骨の一部）。骨の領域２３０は、骨の境界が骨の真ん中よりも高い輝度を有しているため、図５ａおよび図５ｂにおいて輪のような見た目になっている。
【００３２】
ステージ１０８で、多様体決定回路２４は、軸２２０からそれぞれのデータサンプリングパス２２０に沿ってレイをキャスティングすることにより、それぞれのデータサンプリングパス２２０に沿った複数の点のそれぞれで、医用撮像データセットをサンプルする。多様体決定回路２４は、医用撮像データセットの輝度値からそれぞれのデータサンプリングパス２２０に沿って、複数の点のそれぞれに対する輝度値を決定する。複数の点に対する輝度値を決定する際、多様体決定回路２４は、複数のデータサンプリングパス２２０のそれぞれに沿って、医用撮像データの輝度値の変化の度合いを決定する可能性がある。
【００３３】
ステージ１１０から１１４で、多様体決定回路２４は、少なくとも一つの解剖学的構造の形を描出する形状を有する曲面を決定する。ステージ１１０から１１４で、多様体決定回路２４は、少なくとも一つの解剖学的構造の形状を描出する形状を持つ曲面を決定する。少なくとも一つの解剖学的構造は、複数の分離されたサブ構造を有するものであり、今の場合胸郭が複数の肋骨を含んでいる。記載の実施形態では、曲面は多様体４００と呼ばれ、例えば曲面または複合３Ｄ多様体と多様体４００とは、全ての関心解剖学的構造、本実施形態では全ての肋骨と交差する。
【００３４】
問題の解剖学的構造または複数の構造の形状の任意のその他適切な曲面、例えば任意の適切な多様体は、代替的な実施形態に従って決定される可能性があり、また曲面は任意の適切な二次元または三次元形状の可能性もある。曲面は最小、ごくわずかか全く厚みがなく、このような曲面は曲面二次元形状と考えられ得る。代替的に、曲面がゼロではない、またわずかでない、または最小の厚みよりもより大きな厚みを有してもよく、このようにスラブなど三次元形状でもよく、例えば湾曲したスラブである。このようなスラブは、例えば、一ボクセルよりも、より厚みのある可能性がある。対象的に、全く厚みや幅のないラインは、ラインが三次元空間を通して曲面パスに従う場合であっても、一次元形状と考えられ得る。
【００３５】
図１の実施形態において、多様体決定回路２４は曲面を決定するために使用されるが、この場合曲面は、多様体４００である。代替的な実施形態において、任意の適切な回路は、曲面を決定するために使用されてもよい。
【００３６】
ステージ１１０で、少なくとも一本の肋骨で交わるそれぞれのデータサンプリングパス２２０に対し、多様体決定回路２４は、データサンプリングパス２２０に沿ったレイキャストが肋骨と交差する最初の点（つまり、軸に最も近い点）と、データサンプリングパス２２０に沿ったレイキャストが肋骨と交差する最後の点（つまり、軸から最も遠い点）と、を決定する。多様体決定回路２４は、データサンプリングパス２２０に対する推定される多様体位置を決定する。この推定される多様体位置は、最初の点と最後の点との中間点である。この中間点は、レイと肋骨との交点の始まりと終わりの中間である。その他の実施形態において、推定される多様体位置を決定するのに、異なる方法が使用されてもよい。例えば、実施形態の中には、レイと肋骨との間の交差への終点がないものがある。実施形態には、推定される多様体位置は、レイと肋骨との間の最初の交点に関連して、決定されるものもある。
【００３７】
図５ａは、データサンプリングパス２２０ａについて決定された、最初の点２４０、最後の点２４２、そして中間点２４４と、データサンプリングパス２２０ｂについて決定された、最初の点２５０、最後の点２５２、中間点２５４を示している。
【００３８】
中間点は、骨の骨髄の内にある可能性がある。ある状況下では、最初の点と最後の点とのそれぞれが中間点よりも好適に正義されている可能性があるので、中間点を探すために最初の点と最後の点とを使用することは有益な場合がある。最初の点と最後の点とのそれぞれは、骨の外側部分が骨髄よりも大きな輝度を持つだろうから、最初の点と最後の点とは中間点よりもよりコントラストを持つだろう。
【００３９】
本実施形態において、中間点を決定するために最初の点と最後の点とを使用することによる多様体の決定は、肋骨の外側部分は肋骨の内側部分よりも大きな輝度を有することに伴い、肋骨の厚さにわたる輝度の値の予想プロフィールに基づく。その他の実施形態において、多様体の決定は、任意の解剖学的構造にわたる任意の適切なパラメータの予想プロフィールに基づいていてもよい。
【００４０】
本実施形態において、多様体決定回路２４は、骨を描出する輝度を有する点を特定するために輝度しきい値を使用することで、また骨を描出する輝度を有する最初の点と最後の点とを検出することで、それぞれのデータサンプリングパス２２０に対する最初の点と最後の点とを決定する。その他の実施形態において、肋骨はステージ１０８の前に医用撮像データセットからセグメントされ、中間点の特定はセグメントされた肋骨を使用して実行される。最初の点と最後の点とは、セグメンテーションにより、肋骨としてラベル付けられた最初の点と最後の点である。更なる実施形態において、関心解剖学的構造の部分である、点の特定の任意の適切な方法が使用されてもよい。
【００４１】
実施形態の中には、中間点を計算する前に、最初の点と最後の点とがフィルタされるものもある。交点は、レイに沿った距離として格納される。このフィルタリングは、格納された距離の近接をその近接と畳み込むことにより、格納された距離に影響を及ぼすことが出来る。フィルタされた距離の中間点は、多様体を明確にするために使用することが出来る。最初の点と最後の点とをフィルタリングすることは、レイに対して選択された、推定される多様体位置が、連続する曲面とその近辺を形成することを確かなものにするために、一役買う可能性がある。
【００４２】
肋骨と交差する、それぞれのデータサンプリングパス２２０に対して推定される多様体位置は、ステージ１１０で決定された。しかし、本実施形態において、データサンプリングパス２２０の中には、肋骨と交差しないものがある。例えば、何本かのデータサンプリングパス２２０は、肋骨の間の隙間を通過する可能性があるからだ。
【００４３】
ステージ１１２で、多様体決定回路２４は、ステージ１１０で中間点として決定された、推定される多様体位置から補間または補外することにより、肋骨と交差しないデータサンプリングパス２２０のそれぞれに対する推定される多様体位置を決定する。どの骨とも交差しないデータサンプリングパス２２０に対し、多様体決定回路２４は、骨に当らない近接のレイとの間の曲面表面上に位置を補間する。
【００４４】
多様体決定回路２４は、曲面をステージ１１０で決定された中間点へフィッティングすることで、補間を実行する可能性がある。多様体決定回路２４は、近接のレイを処理することで、肋骨領域の中についてもっともらしいＣＰＲ（曲面プラナー再フォーマット）表面（即ち、多様体４００の少なくとも一部のもっともらしい推定）を検出する可能性がある。肋骨と交差しないデータサンプリングパス２２０に対する推定される多様体位置決定の方法は、二次元信号再構成のための規格化畳み込みの使用を具備する可能性がある。（例えば、Ａｎｄｅｒｓ  Ｌａｎｄｓｔｒｏｍ  共著、『重点増分規格化畳み込みによる画像再構成、画像分析』（Ｉｍａｇｅ  Ｒｅｃｏｎｓｔｒｕｃｔｉｏｎ  ｂｙ  Ｐｒｉｏｒｉｔｉｚｅｄ  Ｉｎｃｒｅｍｅｎｔａｌ  Ｎｏｒｍａｌｉｚｅｄ  Ｃｏｎｖｏｌｕｔｉｏｎ，  Ｉｍａｇｅ  Ａｎａｌｙｓｉｓ）：１７ｔｈ  スカンジナビア協議会、ＳＣＩＡ  ２０１１、Ｙｓｔａｄ、スウェーデン、２０１１年５月、Ｐｒｏｃｅｅｄｉｎｇｓ、１７６―１８５頁、Ｈ．Ｋｎｕｔｓｓｏｎ、Ｃ．―Ｆ  Ｗｅｓｔｉｎ.『規格化および微分畳み込み』（Ｎｏｒｍａｌｉｚｅｄ  ａｎｄ  ｄｉｆｆｅｒｅｎｔｉａｌ  ｃｏｎｖｏｌｕｔｉｏｎ）  Ｉｎ  Ｐｒｏｃｅｅｄｉｎｇｓ  ｏｆ  ＩＥＥＥ  『コンピュータビジョンおよびパターン認識における社会協議』（Ｓｏｃｉｅｔｙ  Ｃｏｎｆｅｒｅｎｃｅ  ｏｎ  Ｃｏｍｐｕｔｅｒ  Ｖｉｓｉｏｎ  ａｎｄ  Ｐａｔｔｅｒｎ  Ｒｅｃｏｇｎｉｔｉｏｎ）  ＣＶＰＲ‘９３、５１５―５２３頁、１９９３年６月、を参照。）
ステージ１１４で、多様体決定回路２４は、ステージ１１０から１１２の推定される多様体位置に基づいて、多様体４００を決定する。実施形態の中に、多様体４００は、ステージ１１０と１１２とで決定された、全ての推定される多様体位置を含むものがある。実施形態には、多様体決定回路２４が、推定される多様体位置の少なくともいくつかを調節し、その多様体は調整された位置を含むものもある。例えば、多様体決定回路２４は、推定される多様体位置との整合性に対する推定される多様体位置を検査し、もし整合性がないと判明した場合に、少なくとも推定される多様体位置のいくつかを調整する可能性がある。実施形態の中には、多様体決定回路２４は、偽陽性を除去するために、データサンプリングパス２２０のそれぞれに対する候補肋骨領域のリスト作成し、隣接のレイから候補領域を比較するものがある。
【００４５】
中には、決定された多様体４００がフル多様体であり、連続的な表面を具備するかもしれない実施形態がある。その他の実施形態において、多様体４００は部分的な多様体であり、少なくとも一つの非連続的な表面を具備する可能性がある。例えば、多様体４００は、肋骨上に決められるかもしれないし、また肋骨の間の領域に対して決められないかもしれない。本実施形態において、多様体４００は、全ての肋骨と交差する表面であり、胸郭の少なくとも一部の形を表す。その他の実施形態において、三次元多様体４００は、任意の適切な解剖学的構造と交差してもよい。
【００４６】
多様体決定回路２４は、多様体の要求される基準と一致する、点に対するボリュームを検索するように構成されている。図２のステージ１００から１１４に関する上述のプロセスにおいて、多様体決定回路２４は、円筒投影法を使って検索を開始し、円筒投影法と肋骨とによって決められたレイの交点の開始と終わりとの間の中ほどに点を決定する。その他の実施形態において、多様体上に点を決める任意の適切な方法が使われてよい。例えば、異なる投影システムが使われてもよい。本実施形態において、多様体４００は自動的に導出される。その他の実施形態において、多様体４００はユーザ入力を使って決定されてもよい。
【００４７】
多様体決定回路２４は、多様体４００を通過して画像生成回路２８へ達する。
【００４８】
ステージ１１６から１２２において、多様体４００は、多様体４００上の点の輝度値を使って、肋骨の曲面ＭＰＲ画像を取得するために使用される。その肋骨の曲面ＭＰＲ画像は、肋骨の展開されたビューを表す、二次元画像である。投影システム（この場合、傾けられた円筒投影法を有する）は、生成される予定の局面ＭＰＲ画像の点を、多様体４００の表面上の点へマップするために使用される。ステージ１１６で、画像生成回路２８は、投影システムを使用して複数のレイパス５００を決める。レイパス５００の数は、生成されるピクセルの数と対応する。
【００４９】
レイパス５００決定における第一のステップとして、カメラ軸が患者の上下軸に沿って決められる。本実施形態において、決められたカメラ軸は、図３および図４に示された軸２００と同じである。その他の実施形態において、カメラ軸は、多様体４００の決定において使用される、軸２００とは異なるものでもよい。
【００５０】
点の集まりO1…Nは、カメラ軸に沿って決められる。それぞれの点Oiは、出力画像においてピクセルの横の列に対応するようになる。それぞれの点Oiに対し、垂直面Piは決められ、面Piは角度に基づいてチルトされる。本実施形態において、面P1…Nは図３および図４に示される面２１０と同じであり、上記チルト角は図４に示される角度αである。チルト面を使用することで、もしチルト面がチルトしていない場合よりも、結果画像において肋骨をより水平に見えるようになる。
【００５１】
レイパス５００の収集は、Oiを中心として基点Aiを含む面Pi内を回転することで、面Piにおいて定義される。本実施形態において、レイパス５００は、カメラ軸を中心として、実質３６０度に等しい範囲に分配されている。それぞれのレイパス５００は、列に対する出力画像において、異なる水平ピクセルに対応する。
【００５２】
本実施形態において、レイパス５００は、データサンプリングパス２２０と同じである。その他の実施形態において、レイパス５００はデータサンプリングパス２２０とは異なる。例えば、多様体４００を取得するために使用される投影とは、異なる投影を用いてもよい。データサンプリングパス２２０よりも、多いまたは少ないレイ５００があってもよい。例えば、多様体４００を決定するために使用されるデータサンプリングパス２２０の数は、画像におけるピクセル数より多くなるだろうし、従って画像生成回路２８によって使われるレイ５００の数よりも多くなるだろう。
【００５３】
実施形態の中に、ステージ１０６で決定されたデータサンプリングパス２２０は、レイパス５００として使用され、ステージ１１６が省略されるものがある。
【００５４】
ステージ１１８で、画像生成回路２８は、レイパス５００のそれぞれに沿ったレイをキャスティングする。それぞれのレイがどのように処理されるかが、出力画像に後々影響するようになる。
【００５５】
本実施形態において、画像生成回路２８は、それぞれのレイパス５００に対する点を決定し、その点とは、レイパス５００に沿ったレイキャストが多様体と交差する箇所である。本実施形態で、多様体４００は、レイパス５００に沿った方向に厚みがないため、多様体４００とそれぞれのレイとの交点は、多様体４００上の一点５１０である。画像生成回路２８は、従ってレイパス５００の数に対応する数の点５１０のセットを決定する（また、それにより画像のピクセル数に対応する数を決定する。）
本実施形態において、点５１０は、レイパス５００に沿ったレイキャストが多様体４００と交差する点である。その他の実施形態において、多様体４００上の点５１０は、レイキャスティングを含んでも含まなくてもよい、任意の適切な方法を使用して、決定されてもよい。任意の適した投影方法が使われる可能性がある。複合多様体を二次元撮像平面に展開するという問題については、地図投影を取得する際に直面する問題と同様の問題を持つために、複数の意義において考慮した方が良い。地図投影に対して多くの異なるアプローチが知られているからである。所望により、地図投影法に使用する適切な投影方法が、本願に対して使用される可能性もある。
【００５６】
本願実施形態の一方で、レイ５００はデータサンプリングパス２２０と実質的に接触するが、その他の実施形態において、レイ５００とデータサンプリングパス２２０とは異なる。画像生成回路２８は、点５１０のそれぞれで医用撮像データセットをサンプルする。本実施形態において、画像生成回路２８は、医用撮像データセットから点５１０のそれぞれに対するボクセル輝度値を決定する。その他の実施形態において、任意の適切なパラメータの値が点５１０でサンプリングすることにより、決定されてもよい。
【００５７】
図６は、多様体４００にわたるスライスの概略図である。図６は、軸５２０とレイパス５００を示す。レイは、軸５２０からレイパス５００をたどり、多様体４００と交差する。図６の描出は、画像データから多様体４００が生成され、その画像データを通したスライスとオーバーレイされる。
【００５８】
レイは、取得される肋骨の展開されたビューを可能にする、円筒投影法を提供する。レイは、画像におけるポイントを多様体表面上のポイントにマップするために使用される。それぞれのレイは、出力画像の個別のピクセルに対応する。
【００５９】
ステージ１２０で、画像生成回路２８は、ステージ１１８でサンプルしたボクセル輝度値を使用して、複数のレイパス５００のそれぞれに対する個別のピクセル値を決定する。本実施形態において、画像生成回路２８は、そのレイパス５００に沿ったレイキャストが多様体と交差する点５１０についてのボクセル輝度値に、伝達関数を適用することで、それぞれのレイ５００に対するピクセル値を決定する。それぞれの点５１０に対するボクセル輝度値は、医用撮像データセットをサンプリングすることにより、ステージ１１８で決定されている。
【００６０】
その他の実施形態において、それぞれのレイ５００に対するピクセル値を取得する任意の方法が使用されてよく、曲面はスラブの形状であってもよい。そのスラブは、多様体４００にレイに沿った方向の厚みを与えることで、取得される（この場合、放射状方向）。実施形態の中には、スラブの厚みが自動的に決定されるものがある。またある実施形態は、ユーザによって厚みが選択されるものもある。
【００６１】
ある実施形態で、スラブサンプルは、導出された直交／デカルトボリュームを出力するために使われる。この導出された直交／デカルトボリューム中において、新しいボリュームＸＹは、投影された画像空間に対応し、Ｘはレイ方向にあるスラブオフセットにおいて使用されるサンプルのセットである。この導出されたボリュームはその場合、例えば影付きボリュームレンダリング、ＩＰボリュームレンダリング、全体照明レンダリング、ＭＰＲなど、任意のその他現行のボリューム撮像技術と使用することが出来る。
【００６２】
画像生成回路２８は、レイがオリジナル多様体４００と交差する最も近い交点である、それぞれのレイに沿った多数のボクセルまたはその他の点（例えば、３、４、または５つのボクセル）をサンプルする。画像生成回路２８は、それぞれのレイに対するピクセル値を取得するために、スラブの投影を使用する。
【００６３】
実施形態の中に、画像生成回路２８は、それぞれのレイに対しピクセル値を決定するスラブの最大輝度投影を使用するものがある。例えば、それぞれのレイ５００のピクセル値は、レイ５００について画像生成回路２８によってサンプルされた、３、４、または５つのボクセル値の最大値に基づいてもよい。画像生成回路２８は、最終ピクセル色を決定するために、スラブと交差するそれぞれのレイの一部に沿ったＩＰ集積を実行し得る。実施形態の中には、画像生成回路２８は、ピクセル値がスラブにおけるレイの一部分にわたる信号変化率に依存する、分散投影法を使用する。一般的に、サンプル点の収集、レイ方向におけるそれぞれのオフセットは、スラブ画像を生成するのに使用される可能性があり、ＩＰ投影方法は出力結果を決定するために使われる可能性がある。スラブの厚みにわたる任意の適切な投影が使われてよく、実施形態は最大輝度投影の使用に限られたものではない。
【００６４】
図２のプロセスに再び戻って、ステージ１２２で画像生成回路２８は、ステージ１２０で決定されたピクセル値を使って、表示スクリーン１６上に表示するための二次元出力画像を生成する。二次元画像のそれぞれのピクセルは、個別のレイパス５００に対応する。レイパス５００が傾けられた円筒座標系で定義されるので、二次元画像は、少なくとも何本かの肋骨が実質的に水平に見える可能性のある箇所において、肋骨の展開されたビューを提供する。
【００６５】
本実施形態において、ユーザが関心領域上に焦点を合わせられるように、ユーザは、拡大縮小ツールを使って、出力画像の見た目をコントロールすることができる。画像生成回路２８は、ユーザ入力に従って、二次元出力画像をアップデートしてよい。同じデータの展開されたビューと体軸ＭＰＲビューとの間のナビゲーションを可能とする、追加の三角測量ツールが使用されてもよい。例えば、ユーザは、体軸ＭＰＲビューをクリックしてよいし、展開されたビュー上の対応する点がハイライトされてもよいし、その逆もまた然りである。
【００６６】
二次元画像は、患者の肋骨の全てと交差する、発展型曲面ＭＰＲビュー（advanced curved MPR）（またはＣＰＲビュー、曲面再構成ビュー）として、記載される可能性がある。発展型曲面ＭＰＲビューによって、全ての肋骨は同時に視覚化できるようになるだろう。同時に全ての肋骨を視覚化することで、骨折や病変が発生した肋骨を臨床医が特定しやすくなるだろう。体軸ＭＰＲビューまたは斜断面ＭＰＲビューのセットを通して段階を踏むよりもむしろ、曲面ＭＰＲビューを使うことで同時に全ての肋骨を見る方が、臨床医にとってより速いだろう。
【００６７】
二次元画像は、二段階プロセスを使って、取得される。二段階プロセスの第一段階は、全ての関心構造が交差する多様体決定を具備する。二段階プロセスの第二段階は、最終ピクセル色を決定するためのカラールックアップに使用するために、２Ｄ対象画像におけるポイントを多様体表面の点へのマッピングと、多様体表面のマッピングする点のサンプリングを具備する。
【００６８】
多様体の決定を含む二段階プロセスを使用することによって、肋骨など特定の生体組織がユーザにとって便利な方法によって画像化されることを可能にする。例えば、本実施形態において、二段階プロセスによって、肋骨骨折の効果的な視覚化を可能にするかもしれない。その他の実施形態において、二段階プロセスによって、その他の解剖学的構造の効果的な視覚化を可能にするかもしれない。
【００６９】
肋骨を通過する多様体の決定によって、例えば、肋骨の真ん中を通過する表面の視覚化など、取得される肋骨を通過する表面の視覚化が可能になるだろう。発生したビューは、例えば、多様体の決定なしで最大輝度投影を使う方法によるなど、その他の方法によって発生されるビューにとっては異なるものになる可能性がある。肋骨を通る多様体４００上の点に対する輝度に基づいて画像を示すことは、最大輝度に基づく画像よりも、骨折の改善されたビューを提供し得る。
【００７０】
図２の方法は、変位のない骨折の改善された可視性を提供することがある。薄いＣＰＲビューは、変位のない骨折を可視化するだろう。図２の方法は、もし骨折が背骨近くにあれば、特に便利かもしれない。背骨に近い骨折は、もし何らかの既知の方法を使う場合よりも、より見やすい可能性がある。図２の方法は、骨病変の改善されたビューを提供するだろう。
【００７１】
展開された投影は、図２の処理で使用される。患者の肋骨の全てと交差する複合多様体４００は、自動的に抽出される。画像を作成するために使用される多様体４００は、実質的に連続的であってよい。投影システムは、多様体表面上の点に対応するピクセルの箇所で、二次元画像を表示するのに使用される。
【００７２】
図２の方法は、肋骨と交差する多様体４００の取得を具備する。多様体は、肋骨の中間点に基づくので、肋骨の真ん中を通過または中間を通過するかもしれない。その他の実施形態において、多様体決定回路２４は、多様体の境界線を変化させることで、修正された多様体を取得するように構成されている。多様体の境界変化は、多様体の拡大または縮小を具備し得る。多様体４００を拡大させるまたは縮小させるために、オフセットが多様体４００に適用される可能性がある。
【００７３】
例えば、実施形態の中には、多様体決定回路２４は、拡大された多様体を取得するために、軸に対して多様体外側上の点をオフセットするものがある。図２に関連して上記で述べた通り、画像生成回路２８は、拡大された多様体を使用して、レイキャスティング処理を実行する。画像生成回路２８は、それにより拡大された多様体に基づいて、更なる画像を取得する。拡大された多様体は、それぞれの肋骨の中心を通過するよりも、例えば、それぞれの肋骨の外側表面または外側表面近くなど、それぞれの肋骨の外側部分を通過する可能性がある。
【００７４】
実施形態には、多様体決定回路２４が、縮小された多様体を取得するために、軸に向かって多様体内側上の点をオフセットするものもある。縮小された多様体は、例えば、それぞれの肋骨の内側表面または内側表面近くなど、それぞれの肋骨の内側部分を通過するだろう。画像生成回路２８は、縮小された多様体に基づく画像を取得する。
【００７５】
実施形態の中には、多様体決定回路２４が、ユーザにより選択されてきたオフセット値を具備する、ユーザ入力を受信するものがある。多様体決定回路２４は選択されたオフセットを使用して修正された多様体（この多様体は、拡大されたまたは縮小された多様体であってもよい）を決定し、画像生成回路２８は修正された多様体に対する画像を生成する。
【００７６】
修正された多様体は、薄い多様体（それぞれのレイが一点で多様体と交差してもよい）でもよいし、また上述のようなスラブでもよい。
【００７７】
拡大されたまたは縮小された多様体に基づいて画像を形成することにより、肋骨表面の画像が形成されることを可能にする。肋骨表面を撮像することによって、肋骨全体にわたって拡がらない骨折をユーザが特定することを助けるかもしれない。また部分的な骨折が特定されるかもしれない。
【００７８】
修正された多様体を作成するために使用されたオフセットの選択をユーザに許可することによって、例えば骨折など、関心特徴の最適なビューを提供し得るオフセットを、ユーザが選択することを可能にする。もしユーザがオリジナル画像において骨折を上手く見分けられない場合、更なる画像を作るために、ユーザは多様体を変化させてもよい。
【００７９】
実施形態の中には、多様体決定回路２４が、それぞれが異なるオフセットを有する、複数の修正された多様体を自動的に決定する。画像生成回路２８は、それぞれが異なる修正された多様体を使用する、画像のシリーズを自動的に決定する。例えば、画像のシリーズは、肋骨の全体的な厚みをカバーするだろう。画像のシリーズを通して段階的に踏むことにより、ユーザは肋骨の外側表面から内側表面まで、肋骨の最大拡張を観察することが可能になる。このような画像のシリーズを自動的に提供することによって、ユーザにとっての使い勝手を改善するかもしれない。肋骨の最大厚みを通して、画像の段階を踏む能力は、例えば表面骨折など、部分骨折の検出に改善をもたらす可能性がある。
【００８０】
図２を参考にした上述の実施形態において、医用画像データセットは、多様体４００を決定するために使用される。多様体決定回路２４は、フィルタされたおよび／またはリサンプルされたデータセットを取得するために、医用撮像データセットをフィルタするおよび／またはリサンプルし、多様体４００の決定においてフィルタされたおよび／またはリサンプルされたデータセットを使用する。医用撮像データセットのフィルタリングは、ノイズ低減フィルタおよび／またはスペクル提言フィルタを使用することを具備するだろう。フィルタされたおよび／またはリサンプルされたデータセットを使用することにより、肋骨を検出するために使用されるアルゴリズムを、より強固なものにする可能性がある。
【００８１】
本実施形態において、肋骨は輝度しきい値を使用して検出される。フィルタリングは、高い値のノイズボクセルの数を減らす可能性があるので、もしオリジナル医用撮像データセットよりもフィルタされたデータセットが肋骨ＣＰＲ表面（つまり、多様体４００）を検出するために使用されれば、レイは高い値のノイズボクセルに当たることが少なくなるかもしれない。
【００８２】
高い値のノイズボクセルは、ある状況下において解析フェーズを混乱させる可能性がある。例えば、高い値のノイズボクセルが骨のボクセルにも見えるような、輝度を有するかもしれないのだ。もし高い値のノイズボクセルが存在するならば、肋骨を描出するボクセルを正しく特定するのはより難しくなるだろうし、従って肋骨との交差の中間点を見つけるのもより難しくなるだろう。
【００８３】
実施形態の中に、多様体４００の決定は、フィルタされたおよび／またはリサンプルされたデータセットを使用してレイ５００の交点５１０をサンプリングして実行され、多様体４００は、オリジナル且つフィルタされていない、医用撮像データセットを使って、実行されるものがある。多様体４００の決定に対してフィルタされたおよび／またはリサンプルされたデータセットを使用することと、点のサンプリングに対するフィルタされていないデータセットとを使用することで、もしフィルタされたおよび／またはリサンプルされたデータセットが全てのステージに対して使われた場合、またはもしフィルタされていないデータセットが全てのステージに対して使われた場合、に生み出された画像よりも、ある状況下においてはより良い画像を生み出す可能性がある。
【００８４】
例えば、多様体４００の決定に対してフィルタされたデータセットを使用することは、計算から高い値のノイズピクセルを除去することで、より良い多様体４００を生み出すかもしれない。しかし、多様体を伴う交点が決定されるステージにおいて、サンプリングのための位置は多様体において既に決められているので、高い値のノイズピクセルの存在は、それほど重要ではないかもしれない。例えば、フィルタされていないデータセットからのオリジナル輝度値など、正確な輝度値を有することが、より大事になるかもしれない。
【００８５】
実施形態の中には、骨および／または背骨セグメンテーションは、複合多様体および患者の肋骨の検出に制約を課すために使用され得るものがある。また実施形態の中には、骨セグメンテーションが解析から骨でない領域を除外するために使用されるものもある。例えば、多様体決定回路２４は、多様体位置を推定する際、骨を描出する医用画像データセットの部分のみを使う可能性がある。
【００８６】
中には、背骨セグメンテーションが肋骨検出から背骨を除外するために使用される実施形態がある。背骨セグメンテーションの使用は、背骨に近い肋骨の部分上にある骨折の視覚化を試みる際に、特に関連のある可能性がある。背骨を描出するボクセルを除去することは、例えば肋骨の中間点を見つけるなど、多様体位置を推定することをより簡単にするかもしれない。また背骨を描出するボクセルを除去することは、アルゴリズムが誤って背骨を肋骨と特定する可能性を減らすかもしれず、それにより誤って多様体を配置する可能性を減らすかもしれない。
【００８７】
実施形態の中には、背骨セグメンテーションが使用されなければ、多様体は背骨の近くに決定することが難しくなるだろう。ステージ１０８で骨として特定された最初の点は肋骨である可能性があり、骨として特定された最後の点は背骨である可能性があり、その逆もまた然りである。最初の点および最後の点のうち一つが骨であり、その他の点が背骨であれば、多様体は、正しく肋骨の間を通過するよりも、肋骨と背骨の間の柔組織を通過することで結局終わるであろう。もし、背骨がセグメンテーションによって代わりに除去されれば、多様体は一本の肋骨をたどり、それから肋骨の改善された描出を提供する可能性のある、その他の肋骨に飛び移るかもしれない。
【００８８】
実施形態の中には、データセットのセグメントされたバージョンが多様体４００を取得するために使用され、オリジナルのセグメントされていない医用撮像データセットはレイ５００が多様体４００と交差する、点５１０のサンプリングに対して使用される。
【００８９】
フィルタリング、骨セグメンテーションおよび／または背骨セグメンテーションの使用は、肋骨検出の方法をよりロバストなものにする可能性があり、従って多様体の構成方法もよりロバストなものとなる可能性がある。
【００９０】
図１のシステムの更なる操作モードは、図７のフローチャートに概要で描かれている。図７のプロセスにおいて、フィルタされたデータは多様体を決定するために使用される。
【００９１】
多様体決定回路２４は、入力ボリューム６００を受信する。入力ボリューム６００は、医用撮像データセットである。ステージ６１０で、多様体決定回路２４は、フィルタされたボリューム６２０を取得するために、入力ボリューム６００をフィルタする。フィルタされたボリューム６２０は、入力ボリューム６００と比較した場合、ノイズが低減されている。フィルタされたボリューム６２０は、骨を特定するためにおよび／または背骨と肋骨とを見分けるために、セグメントされてもよい。
【００９２】
ステージ６３０で、多様体決定回路２４は、レイと肋骨との交点を決定するために、フィルタされたボリューム６２０と円筒投影情報６３０のセットとを使用する。例えば、円筒投影情報６３０は、図３および図４を参考に上述したように、傾けられた面２１０のセット上のデータサンプリングパス２２０など、データサンプリングパスのセットを具備する。レイと肋骨の交点は、データサンプリングパス２２０に沿ったレイキャストがフィルタされたボリューム６２０において肋骨と交差する箇所で、点のセットを具備する。
【００９３】
肋骨と交差するそれぞれのデータサンプリングパス２２０に対し、フローチャートはステージ６４０へと進む。ステージ６４０で、多様体決定回路２４は、データサンプリングパス２２０が肋骨と交差にわたる、交点間隔を決定する。多様体決定回路２４は、交点間隔の中間点を決定するだろう。
【００９４】
肋骨と交差しないそれぞれのデータサンプリングパスに対し、フローチャートはステージ６５０へと進む。ステージ６５０で、多様体決定回路２４は、交差しないデータサンプリングパスの近隣から位置を補間する。多様体決定回路２４は、例えばすぐ近くのデータサンプリングパスからの中間点など、すぐ近くのデータサンプリングパスからの交点情報を使って、位置を補間する可能性がある。
【００９５】
ステージ６６０で、多様体決定回路２４は、三次元複合多様体を決定するために、ステージ６４０とステージ６５０の出力を使用する。例えば、多様体決定回路２４は、三次元複合多様体を決定するために、ステージ６４０で決定された中間点と、ステージ６５０で決定された位置とを使用するだろう。
【００９６】
ステージ６７０で、画像生成回路２８は、ステージ６６０で決定された多様体とオリジナル入力ボリューム６００とを受信する。画像生成回路２８は、投影情報６３０に依存してレイをキャスティングすることで、投影と多様体の交点を実行し、これにより多様体上に点のセットを決定する。これは、多様体を得て、２Ｄレンダリングされた画像を生成する、画像生成投影である。ステージ６８０で、画像生成回路２８は、ステージ６７０で決定された、点のセットで画像ボリューム６００をサンプルする。
【００９７】
ステージ６９０で、画像生成回路２８は、例えば変換関数を使うことで、サンプルされた点に対応するピクセル値を取得し、ピクセル値を使用して最終レンダリングされた画像を表示する。最終レンダリングされた画像のビューは、必要に応じスラブされたり、またはオフセットされたりする場合がある。例えば、レンダリングされた画像を観察する上で（on）、ユーザはスラブ厚みおよび／またはオフセットを選択するだろうし、そして画像生成回路２８はその厚みおよび／またはオフセットを使用してさらなる画像をレンダリングするだろう。
【００９８】
図２と図７の実施形態は、複合多様体を決定する第１のレイキャスティングと、レンダリングされた画像を決定する第２のレイキャスティングとを具備する。一般的に、任意の適切な方法が、投影情報の取得（例えば、レイのセット）や、少なくとも一つの解剖学的構造を描出する曲面の取得のために、使用され得る。投影情報と曲面との間の交点を取得する、任意の適切な方法が使用されてよい。ボリューメトリック医用撮像データはその後、最終レンダリングされた画像を取得するために、交差の点でサンプルされる。
【００９９】
実施形態は、肋骨の撮像に関して、説明されている。しかし、上述の方法は、例えば複数の分離したサブ構造を具備するなど、任意の適切な解剖学的構造または関心構造に対して使用されてもよい。その方法は、例えば頭蓋骨や骨盤など、展開されたビューが有益だと思われる、任意の骨または複数の骨を画像化するために使われ得る。異なる投影が、異なる解剖学的構造に対して使用されてもよい。図３や図４に関連して上述された通り、肋骨に対して使用される投影は、傾けられた円筒投影として述べられてもよい。頭蓋骨の場合は、球体投影が使用されてもよい。
【０１００】
中には、図２の方法が歯の撮像、または膿瘻の撮像のために使用される実施形態もある。その方法は、少なくとも一本の歯または歯のセットを画像化するために使用され得る。実施形態の中には、図２の方法が柔組織を画像化するために使用される。例えば、図２の方法は、心臓、肝臓または胃などの臓器を画像化するために使用される可能性がある。図２の方法は、心臓の冠動脈を画像化するために使用されてもよい。
【０１０１】
決定された多様体は、少なくとも一つの関心の解剖学的構造と、少なくとも部分的に交わるだろう。多様体は、少なくとも一つの関心の解剖学的構造の内側において、少なくとも部分的に落ち込む（fall）可能性がある。多様体は、少なくとも一つの関心の解剖学的構造によって決められた領域を、少なくとも部分的に囲む可能性がある。
【０１０２】
実施形態に係る方法は、任意の適切なボリューメトリック撮像データ、つまり関心構造が見える任意のボリューメトリック撮像データ、と共に使用され得る。上記の実施形態は、ＣＴ撮像に対して述べられてきたが、この方法は、例えばＭＲ撮像データセットまたはＰＥＴデータセットの使用など、その他のモダリティにおいても使用されてよい。
【０１０３】
実施形態に係る方法は、例えばヒトの被検体または動物の被検体など、任意の適切な被検体の画像取得のために使用されるだろう。この被検体とは、患者であってもよい。例えば、被検体は、医療処置中でもよい。
【０１０４】
特定の実施形態は、患者のスキャン取得が可能なボリューメトリックスキャナを具備する、医用撮像システムを提供する；ボリュームにおいて関心構造と交わる、複合３Ｄ多様体構造の導出が可能な医用処理システム；出力画像におけるピクセルを、ボリューム空間におけるレイにマップする投影システム；そしてレイの交差と、その後対応する交点でボリュームをサンプルする複合三次元多様体と、の検出が可能なレイ処理システム。サンプルされたデータボリュームはその後、レイに対する最終ピクセル色を決定するために、使用される。
【０１０５】
複合三次元多様体は、患者の肋骨の全てと交差するように設定されている可能性がある。それぞれのレイに対し、サンプル点の収集、レイ方向におけるそれぞれのオフセットは、スラブされた画像を生成するために使用されてよい。ＩＰ投影方法が出力結果を決定するために使用される可能性がある。骨および背骨セグメンテーション情報は、複合多様体レイおよび患者の肋骨に対する検査に、制約を課すために使用され得る。ボリューメトリックフィルタリングまたはリサンプリングは、レイと患者の肋骨との交差のための検査を簡素化するために使用されてもよい。
【０１０６】
上述した実施形態において、最終的に表示される画像は、二次元画像に限定されない。例えば、設定した断面（又は多様体）を基準として、肋骨全体を含む或いは肋骨の所定領域を含む厚さを持った領域を設定し、当該厚さを持った領域に含まれる三次元データを用いて最大値投影処理、最小値投影処理、ボリュームレンダリング処理等を実行し、表示するようにしてもよい。
【０１０７】
本明細書では特定の回路が説明されたが、代替的な実施形態においては、これらの回路の一つ以上の機能性が単一の処理リソースまたはその他構成要素によって提供されることも可能で、あるいは単一の回路によって提供される機能性が二つ以上の処理リソースまたはその他の構成要素の組み合わせで提供されることも可能である。単一回路への言及は、多数の構成要素が互いに離れているか否かにかかわらず、単一の回路の機能性を提供する多数の構成要素を包含し、複数回路への言及は、複数の回路の機能性を提供する単一の構成要素を包含する。
【０１０８】
特定の実施形態が説明されたが、これらの実施形態は単に一例として提示されているだけであり、本発明の範囲を限定することを意図としていない。実際、本明細書に説明された新規の方法およびシステムは、他の様々な形態で実施されてもよい。さらに、本明細書に記載される方法およびシステムの形態において、様々な省略、置換、および変さらが本発明の趣旨から逸脱することなく行われてもよい。添付の特許請求の範囲およびその均等物は、本発明の範囲内に入るそのような形態および修正形態を含むように意図されている。
【符号の説明】
【０１０９】
１０…医用画像処理装置、１２…計算装置、１４…ＣＴスキャナ、１６…表示スクリーン、１８…入力装置、２０…メモリ、２２…ＣＰＵ、２４…多様体決定回路、２８…画像生成回路、２００…軸、２１０…面、２２０ａ、２２０ｂ…データサンプリングパス、２３０…骨の領域、２４４…中間点、２５４…中間点、４００…多様体、５００…レイ、５１０…交点、５２０…軸
【図１】
2017164496.tif 000003
【図２】
2017164496.tif 000004
【図３】
2017164496.tif 000005
【図４】
2017164496.tif 000006
【図５ａ】
2017164496.tif 000007
【図５ｂ】
2017164496.tif 000008
【図６】
2017164496.tif 000009
【図７】
2017164496.tif 000010