画像処理装置、方法およびプログラム
【課題】心臓を撮影して得られた３次元画像データから心筋梗塞の疑いがある部分をより正確に特定する。
【解決手段】心臓を拡縮状態が異なる２つの時相で撮影して得られた３次元画像Ｖ１、Ｖ２を記憶し、３次元画像Ｖ１において心臓の心筋を表す心筋領域を特定し、特定された心筋領域の各画素について、その画素の３次元画像Ｖ１上の位置、および、その画素に解剖学的に対応する３次元画像Ｖ２上の点の位置をそれぞれ取得し、取得された３次元画像Ｖ１、Ｖ２上の位置を用いて、各画素により表される心筋部分の２つの時相間での移動を表す移動ベクトルを取得し、各画素にその画素により表される心筋部分の移動ベクトルを与えてなるベクトル場を各画素において空間的に微分した微分値を取得し、出力する。
【解決手段】心臓を拡縮状態が異なる２つの時相で撮影して得られた３次元画像Ｖ１、Ｖ２を記憶し、３次元画像Ｖ１において心臓の心筋を表す心筋領域を特定し、特定された心筋領域の各画素について、その画素の３次元画像Ｖ１上の位置、および、その画素に解剖学的に対応する３次元画像Ｖ２上の点の位置をそれぞれ取得し、取得された３次元画像Ｖ１、Ｖ２上の位置を用いて、各画素により表される心筋部分の２つの時相間での移動を表す移動ベクトルを取得し、各画素にその画素により表される心筋部分の移動ベクトルを与えてなるベクトル場を各画素において空間的に微分した微分値を取得し、出力する。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、心臓を撮影して得られた３次元画像データから心筋梗塞の疑いがある部分を特定するための画像処理装置、方法、およびプログラムに関する。
【背景技術】
【０００２】
心臓は心筋の周期的な収縮・拡張運動によって血液を体内に循環させる。心筋の運動が健常である部分は拡張期末期に壁厚が最小になり、その後収縮期にかけて壁厚が増加する一方、心筋梗塞部分においては収縮期における壁の増加が健常部分と比べ低下することから、従来より、心臓を時系列的に撮影して得られた複数の３次元画像データから心筋の壁厚の変化量を測定し、その測定結果を用いて心筋梗塞の疑いがある部分を特定する技術が知られている（特許文献１参照）。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特表２０１０−５０２２５５号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかし、上記従来技術では、心筋壁の各位置における壁厚を用いてその位置における内壁から外壁に亘る壁全体を対象として心筋梗塞の疑いがある部分であるか否かを評価しているため、内壁から外壁にかけて心筋梗塞部分と健常部分が混在している場合には、心筋梗塞の疑いがある部分を正しく特定できないという問題がある。
【０００５】
本発明は、上記事情に鑑み、心筋梗塞の疑いがある部分をより正確に特定できる画像処理装置、方法およびプログラムを提供することを目的とするものである。
【課題を解決するための手段】
【０００６】
本発明の画像処理装置は、心臓を拡縮状態が異なる２つの時相で撮影して得られた第１および第２の３次元画像を記憶する画像記憶手段と、第１の３次元画像において心臓の心筋を表す心筋領域を特定する心筋領域特定手段と、特定された心筋領域の各画素について、該画素の第１の３次元画像上の位置、および、該画素に解剖学的に対応する第２の３次元画像上の点の位置をそれぞれ取得する位置取得手段と、前記各画素について取得された第１および第２の３次元画像上の位置を用いて、該画素により表される心筋部分の２つの時相間での移動を表す移動ベクトルを取得するベクトル取得手段と、前記各画素に該画素により表される心筋部分の前記取得された移動ベクトルを与えてなるベクトル場を前記各画素において空間的に微分した微分値を取得する微分値取得手段と、取得された微分値を出力する出力手段とを備えたことを特徴とするものである。
【０００７】
ここで、ベクトル場を空間的に微分するとは、ベクトル場の空間が２以上の座標軸により定義されているとき、その少なくとも１つの座標軸に関して微分することを意味する。
【０００８】
また、ベクトル場をある画素において微分した微分値というのは、その画素のベクトルと近傍画素のベクトルとの差分を意味する。
【０００９】
また、出力とは、データ記録装置に対する出力（メディアへの記録）、表示装置に対する出力（モニタ表示）、プリンタ装置への出力（プリントアウト）等を広く意味する。
【００１０】
上記装置において、出力手段は、前記取得された微分値が所定の閾値より小さい画素領域を他の領域と区別可能に表示するものであってもよい。
【００１１】
この表示は、３次元画像を表示する場合に限らず、３次元画像の一断層面を表示する場合にも適用され得る。この表示の一例として、３次元画像の一断層面を、微分値が所定の閾値より小さい画素領域に該当する部分に他の領域とは異なる色調を付与して表示することが考えられる。
【００１２】
また、出力手段は、前記取得された微分値を、画素毎の微分値の大きさの相違を視覚的に識別可能な態様で表示するものであってもよい。
【００１３】
また、前記各画素の位置が、座標系(Ｘ,Ｙ,Ｚ)により定義されるものであり、座標（ｘ、ｙ、ｚ）に位置する画素における微分値Ｄが、下記式（１）により求められるものであってもよい。
【数１】
ここで、Ｖ（ｘ、ｙ、ｚ）は、座標（ｘ、ｙ、ｚ）に位置する画素における移動ベクトルである。
【００１４】
また、本発明の画像処理方法は、上記画像処理装置の各手段が行う処理を、少なくとも１台のコンピュータにより実行する方法である。
【００１５】
また、本発明の画像処理プログラムは、上記画像処理装置の各手段が行う処理を、少なくとも１台のコンピュータに実行させるプログラムである。このプログラムは、ＣＤ−ＲＯＭ，ＤＶＤなどの記録メディアに記録され、またはサーバコンピュータに付属するストレージやネットワークストレージにダウンロード可能な状態で記録されて、ユーザに提供される。
【発明の効果】
【００１６】
本発明の画像処理装置、方法およびプログラムによれば、心臓を拡縮状態が異なる２つの時相で撮影して得られた第１および第２の３次元画像を記憶し、第１の３次元画像において心臓の心筋を表す心筋領域を特定し、特定された心筋領域の各画素について、該画素の第１の３次元画像上の位置、および、該画素に解剖学的に対応する第２の３次元画像上の点の位置をそれぞれ取得し、前記各画素について取得された第１および第２の３次元画像上の位置を用いて、該画素により表される心筋部分の２つの時相間での移動を表す移動ベクトルを取得し、前記各画素に該画素により表される心筋部分の前記取得された移動ベクトルを与えてなるベクトル場を前記各画素において空間的に微分した微分値を取得し、取得された微分値を出力することができる。この微分値は、心筋の各部分における心筋の運動がその近傍部分における心筋の運動に対してどの程度の独立性を有しているかを表すものであることから、医師等のユーザは、この微分値を用いて、心筋の運動を３次元画像の画素単位で評価することができ、心筋梗塞の疑いがある部分をより正確に特定することができる。
【００１７】
上記装置、方法およびプログラムにおいて、微分値の出力が、取得された微分値が所定の閾値より小さい画素領域を他の領域と区別可能に表示するものである場合には、その表示によって、医師等のユーザは、心臓の心筋梗塞の疑いがある部分を容易に把握することができる。
【００１８】
また、微分値の出力が、取得された微分値を、画素毎の微分値の大きさの相違を視覚的に識別可能な態様で表示するものである場合にはその表示によって、医師等のユーザは、心筋の各部分における心筋の運動がその近傍部分における心筋の運動に対してどの程度の独立性を有しているかを容易に把握することができる。
【図面の簡単な説明】
【００１９】
【図１】本発明の画像処理装置の概略構成を示す図
【図２】画像処理装置が実行する処理の概要を示す図
【図３】収縮期の心臓が撮影された３次元画像を示す図
【図４】拡大期の心臓が撮影された３次元画像を示す図
【図５】ベクトル取得処理によって取得された移動ベクトルの一例を示す図
【図６】微分値取得処理によって取得された微分値の一例を示す図
【図７】表示制御処理を説明するための図
【発明を実施するための形態】
【００２０】
以下、図面を参照して、本発明の画像処理装置、方法、およびプログラムの実施形態について説明する。
【００２１】
本発明の画像処理装置は、一台のコンピュータに、画像処理プログラムをインストールしたものである。コンピュータは、診断を行う医師が直接操作するワークステーションやパソコンでもよいし、もしくは、それらとネットワークを介して接続されたサーバコンピュータでもよい。画像処理プログラムは、ＤＶＤ，ＣＤ−ＲＯＭ等の記録メディアに格納されて配布され、その記録媒体からコンピュータにインストールされる。もしくは、ネットワークに接続されたサーバコンピュータの記憶装置、あるいはネットワークストレージに、外部からアクセス可能な状態で記憶され、要求に応じて医師が使用するコンピュータにダウンロードされ、インストールされる。
【００２２】
図１は、ワークステーションに画像処理プログラムをインストールすることにより実現された画像処理装置の概略構成を示す図である。同図が示すように、画像処理装置１は、標準的なワークステーションの構成として、ＣＰＵ２、メモリ３およびストレージ４を備えている。また、画像処理装置１には、ディスプレイ５と、マウス６等の入力装置が接続されている。
【００２３】
ストレージ４には、心臓Ｈを一心拍周期内の複数の時相において撮影して得られた一連の３次元画像が記憶されている。この各３次元画像は、ＣＴ装置、ＭＲＩ装置、超音波診断装置などで撮影されたマルチスライス画像をボクセル（画素）に分割し、３次元座標空間内に配列してなるものであり、各ボクセルの位置は、被写体の左右方向をＸ軸、前後方向をＹ軸、上下方向をＺ軸とする３次元座標系で定義され、各ボクセルの画素値は、そのボクセルの位置の座標と関連づけられている。
【００２４】
また、メモリ３には、画像処理プログラムと画像処理プログラムが参照するデータ（処理パラメータ等）が記憶されている。画像処理プログラムは、ＣＰＵ２に実行させる処理として、心筋領域特定処理、位置取得処理、ベクトル取得処理、微分値取得処理および表示制御処理を規定している。そして、ＣＰＵ２がプログラムに従いこれらの処理を実行することで、汎用のワークステーションは、心筋領域特定手段、位置取得手段、ベクトル取得手段、微分値取得手段および表示制御手段として機能することになる。
【００２５】
以下、画像処理装置１が実行する各処理について詳細に説明する。図２は、画像処理装置１が実行する処理の流れを示す図である。画像処理装置１は、まず、ストレージ４から心臓の収縮した状態を撮影して得られた３次元画像Ｖ１と拡張した状態を撮影して得られた３次元画像Ｖ２をそれぞれメモリ３にロードする。図３は３次元画像Ｖ１を、図４は３次元画像Ｖ２を示す図である。そして、３次元画像データＶ１において心臓Ｈの心筋を表す心筋領域Ｍを特定する心筋領域特定処理１１を実行する。
【００２６】
心筋領域特定処理１１において、心筋領域Ｍは、その領域の輪郭を決定することにより特定される。三次元形態画像７を構成する各ボクセルの画素値について、心筋領域の輪郭らしさを表す特徴量を算出し、算出した特徴量を機械学習により予め取得された評価関数に基づいて評価することで、そのボクセルが心筋領域の輪郭を表すものであるか否かを判断する。この判断を繰り返すことにより、心筋領域全体の輪郭を表すボクセル群が抽出される。本実施形態では、評価関数の取得にアダブーストアルゴリズムを用いている。なお、心臓領域の抽出は、他のマシンラーニング法や統計解析法、例えば線形判別法やニューラルネットワーク、サポートベクターマシン等を用いて行ってもよい。
【００２７】
続いて、画像処理装置１は、位置取得処理１２を実行する。位置取得処理１２では、３次元画像Ｖ１とＶ２を位置合わせすることによって、心筋領域Ｍ上のボクセルＰそれぞれに対して、そのボクセルＰの３次元画像Ｖ１上の位置（座標）、および、ボクセルＰに解剖学的に対応する３次元画像Ｖ２上の点Ｐ´の位置（座標）をそれぞれ取得する。
【００２８】
位置合わせは、例えば、３次元空間内に一定の間隔で３次元の格子状にＢ−Ｓｐｌｉｎｅの制御点を配置し、この制御点を移動させることで、３次元画像Ｖ１およびＶ２における心臓の形状同士が一致するように３次元空間を変形させることにより行うことができる。なお、３次元空間の変形により３次元画像同士を位置合わせする方法の詳細は、例えば、特開２００８−２８９７９９号公報に開示されている。
【００２９】
また、位置合わせは、たとえば特表２００５―５２８９７４号公報に記載されているように、生理学的拘束条件を用いて数学的にモデル化された非剛体心臓運動モデルを、各画像Ａ，Ｂに適用することによって行うようにしてもよい。
【００３０】
続いて、画像処理装置１は、ベクトル取得処理１３を実行する。ベクトル取得処理１３では、心筋領域Ｍ上のボクセルＰそれぞれについて、そのボクセルにより表される心筋部分の移動ベクトルＶを取得する。各ボクセルＰについての移動ベクトルＶは、ボクセルＰに解剖学的に対応する３次元画像Ｖ２上の点Ｐ´の位置ベクトルからボクセルＰの位置ベクトルを減算することで求められる。ここで、位置ベクトルとは、それぞれの点またはボクセルの位置をベクトルとして表現したものをいう。
【００３１】
図５は、ベクトル取得処理１３によって取得された移動ベクトルＶを示す概念図である。なお、図５においては、理解を容易にするため、一部の移動ベクトルのみをその各ボクセルを基点とする矢印で示している。このように、心筋領域Ｍ上のすべてのボクセルについて移動ベクトルが取得されると、各ボクセルにそのボクセルにより表される心筋部分の移動ベクトルを与えてなるベクトル場が形成される。
【００３２】
続いて、画像処理装置１は、微分値取得処理１４を実行する。微分値取得処理１４では、心筋領域Ｍ上のボクセルそれぞれについて、下記式（１）により、微分値Ｄを取得する。
【数２】
ここで、Ｖ（ｘ、ｙ、ｚ）は、座標（ｘ、ｙ、ｚ）に位置するボクセルにおける移動ベクトルである。
【００３３】
これにより、各ボクセルについて取得された微分値Ｄは、そのボクセルの位置の座標と関連づけられて、メモリ３に記憶される。図６は、微分値取得処理１４によって取得された微分値を示す概念図である。なお、図６においては、理解を容易にするため、一部の微分値のみを、そのボクセルを基点とする矢印で示している。
【００３４】
なお、ここでは、微分値取得処理１４が、式（１）により微分値を取得するものである場合について説明したが、微分値取得処理１４は、上述したようなベクトル場をその各ボクセル位置において空間的に微分して微分値を取得するものであればよく、その具体的な微分方法は上記方法に限定されない。
【００３５】
また、ここでは、Ｘ、ＹまたはＺ軸方向にそれぞれ隣接する３個の隣接ボクセルの移動ベクトルを用いて微分値を取得する場合について説明したが、Ｘ、Ｙ、Ｚ軸のうち１または２つの座標軸方向に隣接するボクセルを用いて微分値を取得するようにしてもよいし、隣接ボクセルではなく１ボクセル以上離れた位置に存在する近傍ボクセルの移動ベクトルを用いて微分値を取得するようにしてもよい。
【００３６】
続いて、画像処理装置１は、表示制御処理１５を実行する。表示制御処理１５では、
メモリ３から３次元画像Ｖ１の各ボクセルに関連付けられている微分値Ｄを取得し、取得された微分値Ｄを所定の閾値と比較することにより、微分値が所定の閾値より小さいボクセル領域ＭＩを特定する。そして、例えば図７に示すように、３次元画像Ｖ１を、領域ＭＩに他の領域とは異なる色調を付与してディスプレイ５に表示させる。ここで、閾値としては、ほぼゼロに近い値を与える。
【００３７】
なお、ここでは、表示制御処理１５が、３次元画像Ｖ１を、微分値が所定の閾値より小さい画素領域ＭＩを他の領域と区別可能に表示するものである場合について説明したが、これに限らず、３次元画像Ｖ２、同じ心臓を撮影して得られた他の３次元画像、またはそのいずれの３次元画像の一断層面において、微分値が所定の閾値より小さい画素領域に該当する部分を他の部分と区別可能に表示した画像を表示するものであってもよい。
【００３８】
また、取得された微分値を、ボクセル毎の微分値の大きさの相違を視覚的に識別可能な態様で表示するものであってもよい。たとえば、心臓の３次元形状を表す画像の心筋を表す各ボクセルに、そのボクセルにおける微分値の大きさに予め対応付けられている色を配した画像を生成、表示することができる。
【００３９】
以上に説明したとおり、本実施形態の画像処理装置１によれば、心筋領域Ｍの各ボクセルにより表される心筋部分の運動がその近傍部分における心筋の運動に対してどの程度の独立性を有しているかを表す、移動ベクトルの微分値を出力しているので、医師等のユーザは、この微分値を用いて心筋の運動を３次元画像のボクセル単位で評価することができる。これにより、心筋梗塞の疑いがある部分をより正確に特定することができる。
【００４０】
また、画像処理装置１によれば、微分値が所定の閾値より小さい画素領域ＭＩを他の領域と区別可能に表示しているので、その表示によって、医師等のユーザは、心臓の心筋梗塞の疑いがある部分を容易に把握することができる。
【００４１】
なお、上記実施の形態では、心臓の収縮した状態を撮影して得られた３次元画像Ｖ１と拡張した状態を撮影して得られた３次元画像Ｖ２を用いて本発明による各種処理を行う場合について説明したが、これに限らず、ストレージ４に記憶されている一連の３次元画像のうち心臓を拡縮状態（拡大または収縮している程度）が異なるいずれか２つの３次元画像を用いるようにしてもよい。
【符号の説明】
【００４２】
１ 画像処理装置
２ ＣＰＵ
３ メモリ
４ ストレージ
５ ディスプレイ
６ マウス
Ｍ 心筋領域
Ｄ 微分値
Ｖ 移動ベクトル
Ｖ１、Ｖ２ ３次元画像
【技術分野】
【０００１】
本発明は、心臓を撮影して得られた３次元画像データから心筋梗塞の疑いがある部分を特定するための画像処理装置、方法、およびプログラムに関する。
【背景技術】
【０００２】
心臓は心筋の周期的な収縮・拡張運動によって血液を体内に循環させる。心筋の運動が健常である部分は拡張期末期に壁厚が最小になり、その後収縮期にかけて壁厚が増加する一方、心筋梗塞部分においては収縮期における壁の増加が健常部分と比べ低下することから、従来より、心臓を時系列的に撮影して得られた複数の３次元画像データから心筋の壁厚の変化量を測定し、その測定結果を用いて心筋梗塞の疑いがある部分を特定する技術が知られている（特許文献１参照）。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特表２０１０−５０２２５５号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかし、上記従来技術では、心筋壁の各位置における壁厚を用いてその位置における内壁から外壁に亘る壁全体を対象として心筋梗塞の疑いがある部分であるか否かを評価しているため、内壁から外壁にかけて心筋梗塞部分と健常部分が混在している場合には、心筋梗塞の疑いがある部分を正しく特定できないという問題がある。
【０００５】
本発明は、上記事情に鑑み、心筋梗塞の疑いがある部分をより正確に特定できる画像処理装置、方法およびプログラムを提供することを目的とするものである。
【課題を解決するための手段】
【０００６】
本発明の画像処理装置は、心臓を拡縮状態が異なる２つの時相で撮影して得られた第１および第２の３次元画像を記憶する画像記憶手段と、第１の３次元画像において心臓の心筋を表す心筋領域を特定する心筋領域特定手段と、特定された心筋領域の各画素について、該画素の第１の３次元画像上の位置、および、該画素に解剖学的に対応する第２の３次元画像上の点の位置をそれぞれ取得する位置取得手段と、前記各画素について取得された第１および第２の３次元画像上の位置を用いて、該画素により表される心筋部分の２つの時相間での移動を表す移動ベクトルを取得するベクトル取得手段と、前記各画素に該画素により表される心筋部分の前記取得された移動ベクトルを与えてなるベクトル場を前記各画素において空間的に微分した微分値を取得する微分値取得手段と、取得された微分値を出力する出力手段とを備えたことを特徴とするものである。
【０００７】
ここで、ベクトル場を空間的に微分するとは、ベクトル場の空間が２以上の座標軸により定義されているとき、その少なくとも１つの座標軸に関して微分することを意味する。
【０００８】
また、ベクトル場をある画素において微分した微分値というのは、その画素のベクトルと近傍画素のベクトルとの差分を意味する。
【０００９】
また、出力とは、データ記録装置に対する出力（メディアへの記録）、表示装置に対する出力（モニタ表示）、プリンタ装置への出力（プリントアウト）等を広く意味する。
【００１０】
上記装置において、出力手段は、前記取得された微分値が所定の閾値より小さい画素領域を他の領域と区別可能に表示するものであってもよい。
【００１１】
この表示は、３次元画像を表示する場合に限らず、３次元画像の一断層面を表示する場合にも適用され得る。この表示の一例として、３次元画像の一断層面を、微分値が所定の閾値より小さい画素領域に該当する部分に他の領域とは異なる色調を付与して表示することが考えられる。
【００１２】
また、出力手段は、前記取得された微分値を、画素毎の微分値の大きさの相違を視覚的に識別可能な態様で表示するものであってもよい。
【００１３】
また、前記各画素の位置が、座標系(Ｘ,Ｙ,Ｚ)により定義されるものであり、座標（ｘ、ｙ、ｚ）に位置する画素における微分値Ｄが、下記式（１）により求められるものであってもよい。
【数１】
ここで、Ｖ（ｘ、ｙ、ｚ）は、座標（ｘ、ｙ、ｚ）に位置する画素における移動ベクトルである。
【００１４】
また、本発明の画像処理方法は、上記画像処理装置の各手段が行う処理を、少なくとも１台のコンピュータにより実行する方法である。
【００１５】
また、本発明の画像処理プログラムは、上記画像処理装置の各手段が行う処理を、少なくとも１台のコンピュータに実行させるプログラムである。このプログラムは、ＣＤ−ＲＯＭ，ＤＶＤなどの記録メディアに記録され、またはサーバコンピュータに付属するストレージやネットワークストレージにダウンロード可能な状態で記録されて、ユーザに提供される。
【発明の効果】
【００１６】
本発明の画像処理装置、方法およびプログラムによれば、心臓を拡縮状態が異なる２つの時相で撮影して得られた第１および第２の３次元画像を記憶し、第１の３次元画像において心臓の心筋を表す心筋領域を特定し、特定された心筋領域の各画素について、該画素の第１の３次元画像上の位置、および、該画素に解剖学的に対応する第２の３次元画像上の点の位置をそれぞれ取得し、前記各画素について取得された第１および第２の３次元画像上の位置を用いて、該画素により表される心筋部分の２つの時相間での移動を表す移動ベクトルを取得し、前記各画素に該画素により表される心筋部分の前記取得された移動ベクトルを与えてなるベクトル場を前記各画素において空間的に微分した微分値を取得し、取得された微分値を出力することができる。この微分値は、心筋の各部分における心筋の運動がその近傍部分における心筋の運動に対してどの程度の独立性を有しているかを表すものであることから、医師等のユーザは、この微分値を用いて、心筋の運動を３次元画像の画素単位で評価することができ、心筋梗塞の疑いがある部分をより正確に特定することができる。
【００１７】
上記装置、方法およびプログラムにおいて、微分値の出力が、取得された微分値が所定の閾値より小さい画素領域を他の領域と区別可能に表示するものである場合には、その表示によって、医師等のユーザは、心臓の心筋梗塞の疑いがある部分を容易に把握することができる。
【００１８】
また、微分値の出力が、取得された微分値を、画素毎の微分値の大きさの相違を視覚的に識別可能な態様で表示するものである場合にはその表示によって、医師等のユーザは、心筋の各部分における心筋の運動がその近傍部分における心筋の運動に対してどの程度の独立性を有しているかを容易に把握することができる。
【図面の簡単な説明】
【００１９】
【図１】本発明の画像処理装置の概略構成を示す図
【図２】画像処理装置が実行する処理の概要を示す図
【図３】収縮期の心臓が撮影された３次元画像を示す図
【図４】拡大期の心臓が撮影された３次元画像を示す図
【図５】ベクトル取得処理によって取得された移動ベクトルの一例を示す図
【図６】微分値取得処理によって取得された微分値の一例を示す図
【図７】表示制御処理を説明するための図
【発明を実施するための形態】
【００２０】
以下、図面を参照して、本発明の画像処理装置、方法、およびプログラムの実施形態について説明する。
【００２１】
本発明の画像処理装置は、一台のコンピュータに、画像処理プログラムをインストールしたものである。コンピュータは、診断を行う医師が直接操作するワークステーションやパソコンでもよいし、もしくは、それらとネットワークを介して接続されたサーバコンピュータでもよい。画像処理プログラムは、ＤＶＤ，ＣＤ−ＲＯＭ等の記録メディアに格納されて配布され、その記録媒体からコンピュータにインストールされる。もしくは、ネットワークに接続されたサーバコンピュータの記憶装置、あるいはネットワークストレージに、外部からアクセス可能な状態で記憶され、要求に応じて医師が使用するコンピュータにダウンロードされ、インストールされる。
【００２２】
図１は、ワークステーションに画像処理プログラムをインストールすることにより実現された画像処理装置の概略構成を示す図である。同図が示すように、画像処理装置１は、標準的なワークステーションの構成として、ＣＰＵ２、メモリ３およびストレージ４を備えている。また、画像処理装置１には、ディスプレイ５と、マウス６等の入力装置が接続されている。
【００２３】
ストレージ４には、心臓Ｈを一心拍周期内の複数の時相において撮影して得られた一連の３次元画像が記憶されている。この各３次元画像は、ＣＴ装置、ＭＲＩ装置、超音波診断装置などで撮影されたマルチスライス画像をボクセル（画素）に分割し、３次元座標空間内に配列してなるものであり、各ボクセルの位置は、被写体の左右方向をＸ軸、前後方向をＹ軸、上下方向をＺ軸とする３次元座標系で定義され、各ボクセルの画素値は、そのボクセルの位置の座標と関連づけられている。
【００２４】
また、メモリ３には、画像処理プログラムと画像処理プログラムが参照するデータ（処理パラメータ等）が記憶されている。画像処理プログラムは、ＣＰＵ２に実行させる処理として、心筋領域特定処理、位置取得処理、ベクトル取得処理、微分値取得処理および表示制御処理を規定している。そして、ＣＰＵ２がプログラムに従いこれらの処理を実行することで、汎用のワークステーションは、心筋領域特定手段、位置取得手段、ベクトル取得手段、微分値取得手段および表示制御手段として機能することになる。
【００２５】
以下、画像処理装置１が実行する各処理について詳細に説明する。図２は、画像処理装置１が実行する処理の流れを示す図である。画像処理装置１は、まず、ストレージ４から心臓の収縮した状態を撮影して得られた３次元画像Ｖ１と拡張した状態を撮影して得られた３次元画像Ｖ２をそれぞれメモリ３にロードする。図３は３次元画像Ｖ１を、図４は３次元画像Ｖ２を示す図である。そして、３次元画像データＶ１において心臓Ｈの心筋を表す心筋領域Ｍを特定する心筋領域特定処理１１を実行する。
【００２６】
心筋領域特定処理１１において、心筋領域Ｍは、その領域の輪郭を決定することにより特定される。三次元形態画像７を構成する各ボクセルの画素値について、心筋領域の輪郭らしさを表す特徴量を算出し、算出した特徴量を機械学習により予め取得された評価関数に基づいて評価することで、そのボクセルが心筋領域の輪郭を表すものであるか否かを判断する。この判断を繰り返すことにより、心筋領域全体の輪郭を表すボクセル群が抽出される。本実施形態では、評価関数の取得にアダブーストアルゴリズムを用いている。なお、心臓領域の抽出は、他のマシンラーニング法や統計解析法、例えば線形判別法やニューラルネットワーク、サポートベクターマシン等を用いて行ってもよい。
【００２７】
続いて、画像処理装置１は、位置取得処理１２を実行する。位置取得処理１２では、３次元画像Ｖ１とＶ２を位置合わせすることによって、心筋領域Ｍ上のボクセルＰそれぞれに対して、そのボクセルＰの３次元画像Ｖ１上の位置（座標）、および、ボクセルＰに解剖学的に対応する３次元画像Ｖ２上の点Ｐ´の位置（座標）をそれぞれ取得する。
【００２８】
位置合わせは、例えば、３次元空間内に一定の間隔で３次元の格子状にＢ−Ｓｐｌｉｎｅの制御点を配置し、この制御点を移動させることで、３次元画像Ｖ１およびＶ２における心臓の形状同士が一致するように３次元空間を変形させることにより行うことができる。なお、３次元空間の変形により３次元画像同士を位置合わせする方法の詳細は、例えば、特開２００８−２８９７９９号公報に開示されている。
【００２９】
また、位置合わせは、たとえば特表２００５―５２８９７４号公報に記載されているように、生理学的拘束条件を用いて数学的にモデル化された非剛体心臓運動モデルを、各画像Ａ，Ｂに適用することによって行うようにしてもよい。
【００３０】
続いて、画像処理装置１は、ベクトル取得処理１３を実行する。ベクトル取得処理１３では、心筋領域Ｍ上のボクセルＰそれぞれについて、そのボクセルにより表される心筋部分の移動ベクトルＶを取得する。各ボクセルＰについての移動ベクトルＶは、ボクセルＰに解剖学的に対応する３次元画像Ｖ２上の点Ｐ´の位置ベクトルからボクセルＰの位置ベクトルを減算することで求められる。ここで、位置ベクトルとは、それぞれの点またはボクセルの位置をベクトルとして表現したものをいう。
【００３１】
図５は、ベクトル取得処理１３によって取得された移動ベクトルＶを示す概念図である。なお、図５においては、理解を容易にするため、一部の移動ベクトルのみをその各ボクセルを基点とする矢印で示している。このように、心筋領域Ｍ上のすべてのボクセルについて移動ベクトルが取得されると、各ボクセルにそのボクセルにより表される心筋部分の移動ベクトルを与えてなるベクトル場が形成される。
【００３２】
続いて、画像処理装置１は、微分値取得処理１４を実行する。微分値取得処理１４では、心筋領域Ｍ上のボクセルそれぞれについて、下記式（１）により、微分値Ｄを取得する。
【数２】
ここで、Ｖ（ｘ、ｙ、ｚ）は、座標（ｘ、ｙ、ｚ）に位置するボクセルにおける移動ベクトルである。
【００３３】
これにより、各ボクセルについて取得された微分値Ｄは、そのボクセルの位置の座標と関連づけられて、メモリ３に記憶される。図６は、微分値取得処理１４によって取得された微分値を示す概念図である。なお、図６においては、理解を容易にするため、一部の微分値のみを、そのボクセルを基点とする矢印で示している。
【００３４】
なお、ここでは、微分値取得処理１４が、式（１）により微分値を取得するものである場合について説明したが、微分値取得処理１４は、上述したようなベクトル場をその各ボクセル位置において空間的に微分して微分値を取得するものであればよく、その具体的な微分方法は上記方法に限定されない。
【００３５】
また、ここでは、Ｘ、ＹまたはＺ軸方向にそれぞれ隣接する３個の隣接ボクセルの移動ベクトルを用いて微分値を取得する場合について説明したが、Ｘ、Ｙ、Ｚ軸のうち１または２つの座標軸方向に隣接するボクセルを用いて微分値を取得するようにしてもよいし、隣接ボクセルではなく１ボクセル以上離れた位置に存在する近傍ボクセルの移動ベクトルを用いて微分値を取得するようにしてもよい。
【００３６】
続いて、画像処理装置１は、表示制御処理１５を実行する。表示制御処理１５では、
メモリ３から３次元画像Ｖ１の各ボクセルに関連付けられている微分値Ｄを取得し、取得された微分値Ｄを所定の閾値と比較することにより、微分値が所定の閾値より小さいボクセル領域ＭＩを特定する。そして、例えば図７に示すように、３次元画像Ｖ１を、領域ＭＩに他の領域とは異なる色調を付与してディスプレイ５に表示させる。ここで、閾値としては、ほぼゼロに近い値を与える。
【００３７】
なお、ここでは、表示制御処理１５が、３次元画像Ｖ１を、微分値が所定の閾値より小さい画素領域ＭＩを他の領域と区別可能に表示するものである場合について説明したが、これに限らず、３次元画像Ｖ２、同じ心臓を撮影して得られた他の３次元画像、またはそのいずれの３次元画像の一断層面において、微分値が所定の閾値より小さい画素領域に該当する部分を他の部分と区別可能に表示した画像を表示するものであってもよい。
【００３８】
また、取得された微分値を、ボクセル毎の微分値の大きさの相違を視覚的に識別可能な態様で表示するものであってもよい。たとえば、心臓の３次元形状を表す画像の心筋を表す各ボクセルに、そのボクセルにおける微分値の大きさに予め対応付けられている色を配した画像を生成、表示することができる。
【００３９】
以上に説明したとおり、本実施形態の画像処理装置１によれば、心筋領域Ｍの各ボクセルにより表される心筋部分の運動がその近傍部分における心筋の運動に対してどの程度の独立性を有しているかを表す、移動ベクトルの微分値を出力しているので、医師等のユーザは、この微分値を用いて心筋の運動を３次元画像のボクセル単位で評価することができる。これにより、心筋梗塞の疑いがある部分をより正確に特定することができる。
【００４０】
また、画像処理装置１によれば、微分値が所定の閾値より小さい画素領域ＭＩを他の領域と区別可能に表示しているので、その表示によって、医師等のユーザは、心臓の心筋梗塞の疑いがある部分を容易に把握することができる。
【００４１】
なお、上記実施の形態では、心臓の収縮した状態を撮影して得られた３次元画像Ｖ１と拡張した状態を撮影して得られた３次元画像Ｖ２を用いて本発明による各種処理を行う場合について説明したが、これに限らず、ストレージ４に記憶されている一連の３次元画像のうち心臓を拡縮状態（拡大または収縮している程度）が異なるいずれか２つの３次元画像を用いるようにしてもよい。
【符号の説明】
【００４２】
１ 画像処理装置
２ ＣＰＵ
３ メモリ
４ ストレージ
５ ディスプレイ
６ マウス
Ｍ 心筋領域
Ｄ 微分値
Ｖ 移動ベクトル
Ｖ１、Ｖ２ ３次元画像
【特許請求の範囲】
【請求項１】
心臓を拡縮状態が異なる２つの時相で撮影して得られた第１および第２の３次元画像を記憶する画像記憶手段と、
前記第１の３次元画像において前記心臓の心筋を表す心筋領域を特定する心筋領域特定手段と、
前記特定された心筋領域の各画素について、該画素の前記第１の３次元画像上の位置、および、該画素に解剖学的に対応する前記第２の３次元画像上の点の位置をそれぞれ取得する位置取得手段と、
前記各画素について取得された第１および第２の３次元画像上の位置を用いて、該画素により表される心筋部分の前記２つの時相間での移動を表す移動ベクトルを取得するベクトル取得手段と、
前記各画素に該画素により表される心筋部分の前記取得された移動ベクトルを与えてなるベクトル場を前記各画素において空間的に微分した微分値を取得する微分値取得手段と、
前記取得された微分値を出力する出力手段とを備えたことを特徴とする画像処理装置。
【請求項２】
前記出力手段が、前記取得された微分値が所定の閾値より小さい画素領域を他の領域と区別可能に表示するものであることを特徴とする請求項１記載の画像処理装置。
【請求項３】
前記出力手段が、前記取得された微分値を、前記画素毎の微分値の大きさの相違を視覚的に識別可能な態様で表示するものであることを特徴とする請求項１記載の画像処理装置。
【請求項４】
前記各画素の位置が、座標系(Ｘ,Ｙ,Ｚ)により定義されるものであり、
座標（ｘ、ｙ、ｚ）に位置する画素における前記微分値Ｄが、下記式（１）により求められるものであることを特徴とする請求項１から３のいずれか１項記載の画像処理装置。
【数１】
前記Ｖ（ｘ、ｙ、ｚ）は、座標（ｘ、ｙ、ｚ）に位置する画素における前記移動ベクトル
【請求項５】
心臓を拡縮状態が異なる２つの時相で撮影して得られた第１および第２の３次元画像が記憶されたた記憶媒体から、前記第１および第２の３次元画像を読み出し、
前記第１の３次元画像において前記心臓の心筋を表す心筋領域を特定し、
前記特定された心筋領域の各画素について、該画素の前記第１の３次元画像上の位置、および、該画素に解剖学的に対応する前記第２の３次元画像上の点の位置をそれぞれ取得し、
前記各画素について取得された第１および第２の３次元画像上の位置を用いて、該画素により表される心筋部分の前記２つの時相間での移動を表す移動ベクトルを取得し、
前記各画素に該画素により表される心筋部分の前記取得された移動ベクトルを与えてなるベクトル場を前記各画素において空間的に微分した微分値を取得し、
前記取得された微分値を出力することを特徴とする画像処理方法。
【請求項６】
前記微分値の出力が、前記取得された微分値が所定の閾値より小さい画素領域を他の領域と区別可能に表示するものであることを特徴とする請求項５記載の画像処理方法。
【請求項７】
前記微分値の出力が、前記取得された微分値を、前記画素毎の微分値の大きさの相違を視覚的に識別可能な態様で表示するものであることを特徴とする請求項５記載の画像処理方法。
【請求項８】
コンピュータを、
心臓を拡縮状態が異なる２つの時相で撮影して得られた第１および第２の３次元画像を記憶する画像記憶手段、
前記第１の３次元画像において前記心臓の心筋を表す心筋領域を特定する心筋領域特定手段、
前記特定された心筋領域の各画素について、該画素の前記第１の３次元画像上の位置、および、該画素に解剖学的に対応する前記第２の３次元画像上の点の位置をそれぞれ取得する位置取得手段、
前記各画素について取得された第１および第２の３次元画像上の位置を用いて、該画素により表される心筋部分の前記２つの時相間での移動を表す移動ベクトルを取得するベクトル取得手段、
前記各画素に該画素により表される心筋部分の前記取得された移動ベクトルを与えてなるベクトル場を前記各画素において空間的に微分した微分値を取得する微分値取得手段、および
前記取得された微分値を出力する出力手段として機能させることを特徴とする画像処理プログラム。
【請求項９】
前記出力手段が、前記取得された微分値が所定の閾値より小さい画素領域を他の領域と区別可能に表示するものであることを特徴とする請求項８記載の画像処理プログラム。
【請求項１０】
前記出力手段が、前記取得された微分値を、前記画素毎の微分値の大きさの相違を視覚的に識別可能な態様で表示するものであることを特徴とする請求項８記載の画像処理プログラム。
【請求項１】
心臓を拡縮状態が異なる２つの時相で撮影して得られた第１および第２の３次元画像を記憶する画像記憶手段と、
前記第１の３次元画像において前記心臓の心筋を表す心筋領域を特定する心筋領域特定手段と、
前記特定された心筋領域の各画素について、該画素の前記第１の３次元画像上の位置、および、該画素に解剖学的に対応する前記第２の３次元画像上の点の位置をそれぞれ取得する位置取得手段と、
前記各画素について取得された第１および第２の３次元画像上の位置を用いて、該画素により表される心筋部分の前記２つの時相間での移動を表す移動ベクトルを取得するベクトル取得手段と、
前記各画素に該画素により表される心筋部分の前記取得された移動ベクトルを与えてなるベクトル場を前記各画素において空間的に微分した微分値を取得する微分値取得手段と、
前記取得された微分値を出力する出力手段とを備えたことを特徴とする画像処理装置。
【請求項２】
前記出力手段が、前記取得された微分値が所定の閾値より小さい画素領域を他の領域と区別可能に表示するものであることを特徴とする請求項１記載の画像処理装置。
【請求項３】
前記出力手段が、前記取得された微分値を、前記画素毎の微分値の大きさの相違を視覚的に識別可能な態様で表示するものであることを特徴とする請求項１記載の画像処理装置。
【請求項４】
前記各画素の位置が、座標系(Ｘ,Ｙ,Ｚ)により定義されるものであり、
座標（ｘ、ｙ、ｚ）に位置する画素における前記微分値Ｄが、下記式（１）により求められるものであることを特徴とする請求項１から３のいずれか１項記載の画像処理装置。
【数１】
前記Ｖ（ｘ、ｙ、ｚ）は、座標（ｘ、ｙ、ｚ）に位置する画素における前記移動ベクトル
【請求項５】
心臓を拡縮状態が異なる２つの時相で撮影して得られた第１および第２の３次元画像が記憶されたた記憶媒体から、前記第１および第２の３次元画像を読み出し、
前記第１の３次元画像において前記心臓の心筋を表す心筋領域を特定し、
前記特定された心筋領域の各画素について、該画素の前記第１の３次元画像上の位置、および、該画素に解剖学的に対応する前記第２の３次元画像上の点の位置をそれぞれ取得し、
前記各画素について取得された第１および第２の３次元画像上の位置を用いて、該画素により表される心筋部分の前記２つの時相間での移動を表す移動ベクトルを取得し、
前記各画素に該画素により表される心筋部分の前記取得された移動ベクトルを与えてなるベクトル場を前記各画素において空間的に微分した微分値を取得し、
前記取得された微分値を出力することを特徴とする画像処理方法。
【請求項６】
前記微分値の出力が、前記取得された微分値が所定の閾値より小さい画素領域を他の領域と区別可能に表示するものであることを特徴とする請求項５記載の画像処理方法。
【請求項７】
前記微分値の出力が、前記取得された微分値を、前記画素毎の微分値の大きさの相違を視覚的に識別可能な態様で表示するものであることを特徴とする請求項５記載の画像処理方法。
【請求項８】
コンピュータを、
心臓を拡縮状態が異なる２つの時相で撮影して得られた第１および第２の３次元画像を記憶する画像記憶手段、
前記第１の３次元画像において前記心臓の心筋を表す心筋領域を特定する心筋領域特定手段、
前記特定された心筋領域の各画素について、該画素の前記第１の３次元画像上の位置、および、該画素に解剖学的に対応する前記第２の３次元画像上の点の位置をそれぞれ取得する位置取得手段、
前記各画素について取得された第１および第２の３次元画像上の位置を用いて、該画素により表される心筋部分の前記２つの時相間での移動を表す移動ベクトルを取得するベクトル取得手段、
前記各画素に該画素により表される心筋部分の前記取得された移動ベクトルを与えてなるベクトル場を前記各画素において空間的に微分した微分値を取得する微分値取得手段、および
前記取得された微分値を出力する出力手段として機能させることを特徴とする画像処理プログラム。
【請求項９】
前記出力手段が、前記取得された微分値が所定の閾値より小さい画素領域を他の領域と区別可能に表示するものであることを特徴とする請求項８記載の画像処理プログラム。
【請求項１０】
前記出力手段が、前記取得された微分値を、前記画素毎の微分値の大きさの相違を視覚的に識別可能な態様で表示するものであることを特徴とする請求項８記載の画像処理プログラム。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【公開番号】特開２０１２−６１０１９（Ｐ２０１２−６１０１９Ａ）
【公開日】平成２４年３月２９日（２０１２．３．２９）
【国際特許分類】
【出願番号】特願２０１０−２０５１９４（Ｐ２０１０−２０５１９４）
【出願日】平成２２年９月１４日（２０１０．９．１４）
【出願人】（３０６０３７３１１）富士フイルム株式会社 (25,513)
【Ｆターム（参考）】
【公開日】平成２４年３月２９日（２０１２．３．２９）
【国際特許分類】
【出願日】平成２２年９月１４日（２０１０．９．１４）
【出願人】（３０６０３７３１１）富士フイルム株式会社 (25,513)
【Ｆターム（参考）】
[ Back to top ]