JP5879098B2 - 体内臓器の動画像生成装置および体内臓器の動画像生成方法 - Google Patents

体内臓器の動画像生成装置および体内臓器の動画像生成方法 Download PDF

Info

Publication number
JP5879098B2
JP5879098B2 JP2011239897A JP2011239897A JP5879098B2 JP 5879098 B2 JP5879098 B2 JP 5879098B2 JP 2011239897 A JP2011239897 A JP 2011239897A JP 2011239897 A JP2011239897 A JP 2011239897A JP 5879098 B2 JP5879098 B2 JP 5879098B2
Authority
JP
Japan
Prior art keywords
cross
image
mark
photographing
sectional
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
Application number
JP2011239897A
Other languages
English (en)
Other versions
JP2013097586A (ja
Inventor
嶐一 岡
嶐一 岡
覚 前場
覚 前場
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Aizu
Original Assignee
University of Aizu
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by University of Aizu filed Critical University of Aizu
Priority to JP2011239897A priority Critical patent/JP5879098B2/ja
Publication of JP2013097586A publication Critical patent/JP2013097586A/ja
Application granted granted Critical
Publication of JP5879098B2 publication Critical patent/JP5879098B2/ja
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Magnetic Resonance Imaging Apparatus (AREA)
  • Processing Or Creating Images (AREA)

Description

本発明は、体内臓器の動画像生成装置および体内臓器の動画像生成方法に関し、より詳細には、体内臓器を一定時間間隔で連続的に撮影した複数の撮影画像に基づいて、臓器の所定部位の状態を３次元的な動画像として提供することが可能な体内臓器の動画像生成装置および体内臓器の動画像生成方法に関する。
従来より、体内臓器の状態を医師が判断するために、超音波ＣＴ（Computed Tomography：コンピュータ断層撮影）、弱Ｘ線ＣＴ、ＭＲＩ（Magnetic Resonance Imaging:磁気共鳴画像法）などの一般的な臓器撮影機器による臓器撮影方法を用いて、体内臓器の状態を２次元的な断面画像として撮影する（画像解析により断面画像を求める）方法が知られている。医師は、撮影タイミングあるいは撮影位置の異なる複数の２次元画像を比較・検討することにより、臓器の状態を判断している。
ここで、上述した臓器撮影機器は、身体の一断面の画像、あるいは一方向からの画像を撮影するものであるため、画像全体から診断対象となる体内臓器の特定を行うことが容易ではないという問題があった。
また、上述した臓器撮影機器では、断面画像において明確（鮮明）に映し出す（撮影する）ことが困難な体内臓器が存在する。例えば、弱Ｘ線ＣＴでは、臓器が骨部の内側にあると、その骨によって臓器部分が遮断されてしまうという問題があった。また超音波ＣＴでは、画像の解像度が十分でなく不鮮明な画像になってしまうという問題があった。さらにＭＲＩでは、骨による臓器部分の遮断はないが、臓器の筋肉（例えば心臓の筋肉）や腱（例えば心臓の三尖弁など）の部位が画像において黒く映ってしまい、画像から直接筋肉や腱を特定することが困難であるという問題があった。
ＭＲＩを用いる場合、心臓の心房のように血液で満たされている器官に関しては、その臓器の存在を比較的明確に断面画像に映し出す（撮影する）ことができる。しかしながら、血液で満たされていない腱や筋肉などは、ＭＲＩ画像に鮮明に映し出す（撮影する）ことが困難であった。従って、医師は、明暗レベル調整などを行うことによりＭＲＩ画像の改善を図ったうえで、不鮮明な部位の予測・特定を行い、診断対象となる部位を診断しなければならかった。このため、腱や筋肉などの部位の状態を判断することが困難であり、医師の経験や技量によって判断が異なってしまうおそれがあるという問題があった。
さらに、上述した臓器撮影機器は、臓器の２次元的な画像を撮影するものであるため、本来３次元的な動きを行う臓器の活動状況を３次元的な動きとして映像化することが困難であった。
例えば、心臓の三尖弁の動きを３次元的に観察することができれば、弁が正常に働いているか否かの判断を容易に行うことが可能となる。しかしながら、上述したように三尖弁は腱であって、ＣＴ画像やＭＲＩ画像において鮮明に表示することが困難な部位であるため、３次元的な動きはもちろんのこと、２次元画像においてもその部位を特定することが容易ではなかった。
このため、複数のＣＴ画像から三尖弁の位置を推測して手作業で三尖弁リングの形状を推測する研究がなされている（例えば、非特許文献１参照）。
しかしながら、上述した非特許文献１の手法では、複数のＣＴ画像のそれぞれにおいて三尖弁リングの位置を推測し、位置の推測が行われたそれぞれの複数のＣＴ画像から３次元的な位置を算出して、１つ１つ手作業で三尖弁リングの３次元的な位置を求める必要がある。このため、三尖弁リングの３次元画像を求めるだけで、多大な労力が必要となってしまう。さらに、このようにして求められる三尖弁リングの３次元画像を動画像にするためには、さらなる負担が必要とされ、３次元的な動画像を得ることが大変困難なものであった。
本発明は上記課題に鑑みてなされたものであり、ＭＲＩ画像等の断面画像において必ずしも鮮明に表示されない体内臓器の所定部位を、複数の断面画像において自動的に求めて、該当する所定部位における状態変化を３次元的な動画像として示すことが可能な体内臓器の動画像生成装置および体内臓器の動画像生成方法を提供することを課題とする。
上記課題を解決するために、本発明に係る体内臓器の動画像生成装置は、体内臓器における断面画像の撮影位置を当該断面画像の積層方向へ移動させることができ、かつ、前記断面画像を所定時間間隔で連続的に撮影することが可能な断面画像撮影手段と、該断面画像撮影手段により撮影された前記断面画像を、同一撮影位置において連続して撮影された前記断面画像の撮影順番と、当該断面画像の撮影位置とに関連付けて記録する記録手段と、前記断面画像撮影手段における前記撮影位置および撮影タイミングを制御する撮影制御手段と、前記記録手段に記録された前記断面画像の画像解析を行う画像解析手段と、前記断面画像に対してユーザがマーク位置を設定することにより、前記体内臓器における所定部位の位置をマークによって当該断面画像に記録させることが可能なマーク付加手段とを有し、前記撮影制御手段は、前記断面画像撮影手段に、前記体内臓器の前記断面画像を一の撮影位置において連続的に撮影させ、かつ、当該連続的な前記断面画像の撮影を、前記体内臓器が収まる積層方向範囲において少しずつ撮影位置をずらして当該撮影位置毎に行わせ、前記画像解析手段は、前記記録手段に記録された断面画像を読み出して、一の撮影位置において連続的に撮影された複数の前記断面画像のうち当該撮影位置において最初に撮影された断面画像に対して、前記マーク付加手段により設定されたマーク位置に前記マークを付加し、２次元連続ＤＰを用いて、前記マークが付加された前記断面画像のマーク位置に基づいて、前記マークが付加された前記断面画像の次に撮影された断面画像にマークを付加することにより、連続的に撮影された全ての撮影位置における全ての断面画像に自動的にマークを付加し、前記断面画像に付加された前記マークの重心位置を、当該マークのピクセル群に基づいて全ての断面画像について算出し、算出された重心位置に基づいて、前記断面画像の画像平面および前記積層方向からなる３次元空間における前記体内臓器の前記部位の座標情報を、異なる撮影位置であって前記連続的に撮影された断面画像のうち同じ撮影順番で撮影された断面画像毎に求め、同じ撮影順番で撮影された前記座標情報に基づいて３次元的なスプライン曲線を求めて、撮影順番毎に前記部位の３次元的な形状を求め、前記撮影順番毎に求められた前記部位の３次元的な形状を、前記撮影順番に従って変化させることにより、前記部位の３次元的な形状の時間的な変化を求めて、前記部位における３次元的な動画像を生成することを特徴とする。
また、本発明に係る体内臓器の動画像生成方法は、体内臓器における断面画像の撮影位置を当該断面画像の積層方向へ移動させることができ、かつ、前記断面画像を所定時間間隔で連続的に撮影することが可能な断面画像撮影手段と、該断面画像撮影手段により撮影された前記断面画像を、同一撮影位置において連続して撮影された前記断面画像の撮影順番と、当該断面画像の撮影位置とに関連付けて記録する記録手段と、前記断面画像撮影手段における前記撮影位置および撮影タイミングを制御する撮影制御手段と、前記記録手段に記録された前記断面画像の画像解析を行う画像解析手段と、前記断面画像に対してユーザがマーク位置を設定することにより、前記体内臓器における所定部位の位置をマークによって当該断面画像に記録させることが可能なマーク付加手段とを有する体内臓器の動画像生成装置において、前記撮影制御手段が、前記断面画像撮影手段に、前記体内臓器の前記断面画像を一の撮影位置において連続的に撮影させ、かつ、当該連続的な前記断面画像の撮影を、前記体内臓器が収まる積層方向範囲において少しずつ撮影位置をずらして当該撮影位置毎に行わせ、前記画像解析手段が、前記記録手段に記録された断面画像を読み出して、一の撮影位置において連続的に撮影された複数の前記断面画像のうち当該撮影位置において最初に撮影された断面画像に対して、前記マーク付加手段により設定されたマーク位置に前記マークを付加し、２次元連続ＤＰを用いて、前記マークが付加された前記断面画像のマーク位置に基づいて、前記マークが付加された前記断面画像の次に撮影された断面画像にマークを付加することにより、連続的に撮影された全ての撮影位置における全ての断面画像に自動的にマークを付加し、前記断面画像に付加された前記マークの重心位置を、当該マークのピクセル群に基づいて全ての断面画像について算出し、算出された重心位置に基づいて、前記断面画像の画像平面および前記積層方向からなる３次元空間における前記体内臓器の前記部位の座標情報を、異なる撮影位置であって前記連続的に撮影された断面画像のうち同じ撮影順番で撮影された断面画像毎に求め、同じ撮影順番で撮影された前記座標情報に基づいて３次元的なスプライン曲線を求めて、撮影順番毎に前記部位の３次元的な形状を求め、前記撮影順番毎に求められた前記部位の３次元的な形状を、前記撮影順番に従って変化させることにより、前記部位の３次元的な形状の時間的な変化を求めて、前記部位における３次元的な動画像を生成することを特徴とする。
本発明に係る体内臓器の動画像生成装置および体内臓器の動画像生成方法では、一の撮影位置において連続的に撮影された複数の断面画像のうち、当該撮影位置において最初に（一番初めに）撮影された断面画像に対して、マーク付加手段により設定されたマーク位置にマークが付加され、マークが付加された断面画像の次に撮影された断面画像におけるマークが、当該マークが付加された断面画像のマーク位置に基づいて、２次元連続ＤＰにより付加される。このように、連続的に撮影された全ての撮影位置における全ての断面画像に対して、２次元連続ＤＰを用いて順番にマークを自動的に付加することが可能となる。従って、例えば、断面画像において鮮明（明確）に映し出されていない体内臓器の部位であっても、その部位に該当する部位を連続的に特定してマークを付加することが可能となる。
このようにして、断面画像において鮮明（明確）に映し出されていない体内臓器の部位であっても、部位の位置および状態変化を、連続した断面画像より簡単に特定することができるので、体内臓器の部位の診断結果が医師の経験や技量によって異なってしまうことを低減させることができる。
また、本発明に係る体内臓器の動画像生成装置および体内臓器の動画像生成方法では、断面画像に付加されたマークの重心位置に基づいて体内臓器の部位の位置を求めて、その部位の位置を断面画像の画像平面および積層方向からなる３次元空間における座標位置として求めて、さらに、その座標位置から３次元的なスプライン曲線を求めて、体内臓器の３次元的な形状を求めることができる。
このため、従来のように、医師が、体内臓器の部位に関する２次元的な断面画像に基づいて３次元的な形状を推測することにより部位の形状を判断していた場合に比べて、部位の３次元的な形状の信頼性を向上させることが可能になる。また、形状判断負担の軽減を図ることが可能となり、診断精度の向上を図ることが可能となる。
さらに、本発明に係る体内臓器の動画像生成装置および体内臓器の動画像生成方法では、撮影順番毎に求められた部位の３次元的な形状を、撮影順番に従って変化させることにより、部位の３次元的な形状の時間的な変化を求めて、部位における３次元的な動画像を生成することができる。このため、医師が部位の動きを容易に把握することが可能となり、３次元的な静止形状だけでなく、３次元的な動きを考慮した上で、体内臓器の診断を行うことができ、診断精度の向上を図ることが可能となる。また、３次元的な部位の動画像から、部位の可動体積を容易に算出することができるため、可動体積に基づいて統計的に正常あるいは異常を判断するための判断基準（検定指針など）を構築し、診断に利用することが可能となる。
また、３次元的な動画像を求めることによって診断精度・判断手法に飛躍的な効果を果たすことが期待される体内臓器は心臓である。心臓は、常に鼓動に伴って動いているが、その動きが不自然である部位には何らかの異常が生じている可能性が高い。例えば、心臓の肺動脈弁、三尖弁、僧帽弁、大動脈弁などの弁機能部は、超音波ＣＴ、弱Ｘ線ＣＴ、ＭＲＩなどで鮮明（明確）に撮影することが難しく、断面画像において部位を特定するだけでも困難であることから、その状態を３次元的に求めることが困難であった。このため、本発明に係る体内臓器の動画像生成装置および体内臓器の動画像生成方法により、弁機能部の形状を３次元的に求めるだけでなくその動画像をも生成することにより、心臓弁膜症などの診断精度および診断手法を飛躍的に向上させることが可能である。
なお、断面画像撮影手段とは、例えば、超音波ＣＴ、弱Ｘ線ＣＴ、ＭＲＩなどの一般的な臓器撮影機器を意味する。また、断面画像とは、これらの臓器撮影機器により撮影された解析画像であって、被検者の体内臓器の断面状態をスライス状に示した画像を意味する。
また、２次元連続ＤＰ（Dynamic Programming）とは、対応する２つの画像において共通するピクセルの対応を自動的に求めることが可能な画像解析技術であって、『岡隆一，外２名「連続ＤＰの一般スキームについて−画像スポッティングのための全画素最適マッチング−」，信学技報，電子情報通信学会，IEICE Technical Report, PRMU2010-87,IBISML2010-59 (2010-09)』において公開され、広く知られた技術である。また、２次元連続ＤＰに関しては、特開２０１０−１６５１０４号公報において説明がなされており、既に知られた技術に該当する。
また、上述した体内臓器の動画像生成装置または体内臓器の動画像生成方法において、前記画像解析手段が、前記マークの重心位置を全ての断面画像について算出する前に、当該マークが付加された前記断面画像に対して、ダイレーションとエロージョンとによる画像処理を施すものであってもよい。
このように、画像解析手段が、マークの重心位置を全ての断面画像について算出する前に、当該マークが付加された断面画像に対して、ダイレーションとエロージョンとによる画像処理を施す。この画像処理により、断面画像に付加されたマークを構成するピクセルの孤立点（孤立したピクセル）の除去、不連続なピクセルの接続と穴埋めを行うことが可能となる。このため、マークを示すピクセルをピクセル群として一塊にまとめて明確化することが可能となる。
なお、ダイレーション（dilation）とエロージョン（erosion）とは、モルフォロジ処理と呼ばれる画像処理の代表的な処理である。上述したように、孤立点（孤立したピクセル）の除去、不連続なピクセルの接続と穴埋めのために利用される処理である。
さらに、上述した体内臓器の動画像生成装置または体内臓器の動画像生成方法において、前記撮影制御手段が、前記断面画像撮影手段に前記断面画像を連続的に撮影させる場合に、前記体内臓器を備える被検者の一の鼓動の開始から次の鼓動が開始されるまでの断面画像を連続的に撮影させるものであってもよい。
３次元的な動画像を求めることによって体内臓器の診断結果・判断手法に飛躍的な効果を果たすことが期待される心臓は、常に鼓動に伴って動くものである。このため、心臓の部位における３次元の動画像を生成する場合に、被検者の心臓における一の鼓動の開始から次の鼓動が開始されるまでの断面画像を連続的に撮影して３次元的な動画像を生成することにより、心臓の鼓動に伴う部位の一連の動きを動画像として提供することが可能となる。このため、心臓に関連する部位の診断における精度の向上を図ることが可能となり、従来判断し難かった心臓の疾患などの早期発見・診断を実現することが可能となる。
本発明に係る体内臓器の動画像生成装置および体内臓器の動画像生成方法では、一の撮影位置において連続的に撮影された複数の断面画像のうち、当該撮影位置において最初に（一番初めに）撮影された断面画像に対して、マーク付加手段により設定されたマーク位置にマークが付加され、マークが付加された断面画像の次に撮影された断面画像におけるマークが、当該マークが付加された断面画像のマーク位置に基づいて、２次元連続ＤＰにより付加される。このように、連続的に撮影された全ての撮影位置における全ての断面画像に対して、２次元連続ＤＰを用いて順番にマークを自動的に付加することが可能となる。従って、例えば、断面画像において鮮明（明確）に映し出されていない体内臓器の部位であっても、その部位に該当する部位を連続的に特定してマークを付加することが可能となる。
このようにして、断面画像において鮮明（明確）に映し出されていない体内臓器の部位であっても、部位の位置および状態変化を、連続した断面画像より簡単に特定することができるので、体内臓器の部位の診断結果が医師の経験や技量によって異なってしまうことを低減させることができる。
また、本発明に係る体内臓器の動画像生成装置および体内臓器の動画像生成方法では、断面画像に付加されたマークの重心位置に基づいて体内臓器の部位の位置を求めて、その部位の位置を断面画像の画像平面および積層方向からなる３次元空間における座標位置として求めて、さらに、その座標位置から３次元的なスプライン曲線を求めて、体内臓器の３次元的な形状を求めることができる。
このため、従来のように、医師が、体内臓器の部位に関する２次元的な断面画像に基づいて３次元的な形状を推測することにより部位の形状を判断していた場合に比べて、部位の３次元的な形状の信頼性を向上させることが可能になる。また、形状判断負担の軽減を図ることが可能となり、診断精度の向上を図ることが可能となる。
さらに、本発明に係る体内臓器の動画像生成装置および体内臓器の動画像生成方法では、撮影順番毎に求められた部位の３次元的な形状を、撮影順番に従って変化させることにより、部位の３次元的な形状の時間的な変化を求めて、部位における３次元的な動画像を生成することができる。このため、医師が部位の動きを容易に把握することが可能となり、３次元的な静止形状だけでなく、３次元的な動きを考慮した上で、体内臓器の診断を行うことができ、診断精度の向上を図ることが可能となる。また、３次元的な部位の動画像から、部位の可動体積を容易に算出することができるため、可動体積に基づいて定量的かつ統計的に正常あるいは異常を判断するための判断基準（検定指針など）を構築し、診断に利用することが可能となる。
また、３次元的な動画像を求めることによって診断精度・判断手法に飛躍的な効果を果たすことが期待される体内臓器は心臓である。心臓は、常に鼓動に伴って動いているが、その動きが不自然である部位には何らかの異常が生じている可能性が高い。例えば、心臓の肺動脈弁、三尖弁、僧帽弁、大動脈弁などの弁機能部は、超音波ＣＴ、弱Ｘ線ＣＴ、ＭＲＩなどで鮮明（明確）に撮影することが難しく、断面画像において部位を特定するだけでも困難であることから、その状態を３次元的に求めることが困難であった。このため、本発明に係る体内臓器の動画像生成装置および体内臓器の動画像生成方法により、弁機能部の形状を３次元的に求めるだけでなくその動画像をも生成することにより、心臓弁膜症などの診断精度および診断手法を飛躍的に向上させることが可能である。
実施の形態に係るＭＲＩ装置の概略構成を示したブロック図である。 実施の形態に係るＭＲＩ装置のコンピュータの概略構成を示したブロック図である。 心臓の構造を模式的に示した図である。 実施の形態に係るＣＰＵの処理を示したフローチャートである。 （ａ）は、ＭＲＩ装置により撮影された心臓のスライス画像の一例を示した図であり、（ｂ）は心臓におけるＭＲＩ画像の撮影位置をＬ＝１〜６までずらして撮影する場合の位置関係を示した心臓の模式図である。 （ａ）は、時間ｔ＝１のＭＲＩ画像においてマーク処理を施す前の画像を示し、（ｂ）は、マーク処理を施した後の画像を示している。 （ａ）は時間ｔ＝１のＭＲＩ画像を示し、（ｂ）は、時間ｔ＝２におけるＭＲＩ画像であって、２次元連続ＤＰによりマーク位置が求められた画像を示し、（ｃ）は、（ｂ）のＭＲＩ画像に対してダイレーションおよびエロージョンを実行した後の画像を示している。 （ａ）は、三尖弁の形状を３次元的に示した図であり、（ｂ）は、三尖弁の可動体積を示した図である。
以下、本発明に係る体内臓器の動画像生成装置について、一例であるＭＲＩ（Magnetic Resonance Imaging）装置（磁気共鳴画像法装置）を提示し、図面を用いて詳細に説明する。
図１は、ＭＲＩ装置の概略構成を示したブロック図である。
ＭＲＩ装置１は、静磁場発生磁石２と、高周波送信用コイル３と、高周波受信用コイル４と、傾斜磁場コイル５，５と、コンピュータ６と、シーケンサ７と、高周波発生装置８と、変調器９と、増幅器１０と、傾斜磁場発生装置１１と、増幅器１２と、位相検波器１３と、ＡＤ変換器１４と、表示部１５と、操作入力部（マーク付加手段）１６とを備えている。
静磁場発生磁石２は、被検者Ｘが置かれる空間に均一な静磁場を発生する役割を有している。傾斜磁場コイル５，５は、被検者Ｘが置かれる空間に直交３軸方向の傾斜磁場を発生する役割を有しており、傾斜磁場コイル５，５には、傾斜磁場発生装置１１によって駆動電力が給電されている。
高周波発生装置８は、高周波信号を出力する装置であり、高周波発生装置８から出力される高周波信号は、変調器９において適切な信号に変調され、増幅器１０において増幅されて高周波送信用コイル３に供給される。高周波送信用コイル３に対して高周波信号が供給されることにより、被検者Ｘの生体組織を構成する静磁場中の原子核（プロトン）に核磁気共鳴を起こせる高周波パルスを照射することが可能となる。
高周波送信用コイル３により高周波パルスが被検者Ｘに照射されると、プロトンに核磁気共鳴が起こり、エコー信号などの磁気共鳴信号が出力される。この磁気共鳴信号は、ＭＲＩ装置１の高周波受信用コイル４によって受波され、増幅器１２において増幅された後、位相検波器１３においてアナログ信号として検出される。検出されたアナログ信号は、ＡＤ変換器１４によってデジタル信号の計測データに変換され、コンピュータ６へ出力される。
シーケンサ７は、コンピュータ６の制御指示に応じて、磁気共鳴撮影を行うのに必要な命令を、ＭＲＩ装置１の磁場発生系統（傾斜磁場コイル５、５や傾斜磁場発生装置１１など）や高周波送受信系統（高周波発生装置８、変調器９、増幅器１０、１２、高周波送信用コイル３、高周波受信用コイル４、位相検波器１３、ＡＤ変換器１４など）に送り、これらの動作を制御する役割を有している。
コンピュータ６の制御指示に応じて、ＭＲＩ装置１の磁場発生系統および高周波送受信系統が動作させることにより、被検者Ｘの体内臓器の断面画像を連続的に撮影することが可能となる。また、コンピュータ６の制御指示に従って、断面画像の撮影位置を被検者Ｘの身長方向（撮影される断面画像の積層方向）へ少しずつ移動させて、断面画像を撮影することが可能となっている。
なお、このようにして断面画像を撮影するための処理を行う静磁場発生磁石２、高周波送信用コイル３、高周波受信用コイル４、傾斜磁場コイル５，５、シーケンサ７、高周波発生装置８、変調器９、増幅器１０、傾斜磁場発生装置１１、増幅器１２、位相検波器１３およびＡＤ変換器１４は、本発明に係る断面画像撮影手段に該当する。
表示部１５は、コンピュータ６により画像処理された体内臓器の部位画像および部位動画を表示する役割を有しており、コンピュータ６に接続されている。また、操作入力部１６は、被検者ＸにおけるＭＲＩ画像の撮影位置設定・調整操作および撮影されたＭＲＩ画像に対して所定の画像処理を行うために必要な情報の入力操作など（後述するマーク位置の設定など）を行う役割を有している。
コンピュータ６は、図２に示すように、ＣＰＵ（撮影制御手段、画像解析手段）２０と、ＲＯＭ２１と、ＲＡＭ２２、記憶部（記録手段）２３を有している。ＲＯＭ２１は、ＣＰＵ２０における処理内容を記録したプログラム（例えば、図４に示すフローチャートからなるプログラム）や、ＭＲＩ装置１における初期設定情報などを記録する役割を有している。ＲＡＭ２２は、ＣＰＵ２０の処理において利用されるワークエリアとしての役割を有している。
ＣＰＵ２０は、ＲＯＭ２１に記録されるプログラムに従って、臓器における所定部位のＭＲＩ画像を、所定時間毎に連続的に撮影する処理を行う。また、ＣＰＵ２０は、医師等が操作入力部１６を操作して一枚目のＭＲＩ画像にマークされた所定部位を基準として、その後に撮影された複数のＭＲＩ画像から所定部位の特定を行う処理（該当するマーク位置を設定してマークを追加（記録）する処理）を行う。さらに、ＣＰＵ２０は、所定部位におけるＭＲＩ画像の撮影位置を少しずつずらして撮影を行うことにより、所定部位における３次元的な動画像の生成を行う。
また、ＣＰＵ２０は、処理内容に応じてシーケンサ７に制御指示を出力する。シーケンサ７は、ＣＰＵ２０による制御指示応じて、ＭＲＩ装置１における磁場発生系統（傾斜磁場コイル５、５や傾斜磁場発生装置１１など）や高周波送受信系統（高周波発生装置８、変調器９、増幅器１０、１２、高周波送信用コイル３、高周波受信用コイル４、位相検波器１３、ＡＤ変換器１４など）などの動作を行う。
次に、ＣＰＵ２０の処理について説明する。本実施の形態においては、体内臓器の所定部位の一例として、心臓の三尖弁の動きを動画像として生成する場合について説明をする。ただし、臓器の所定部位として三尖弁の動きを検出することは、単なる一例であって、臓器の所定部位は必ずしも心臓の三尖弁に限定されるものではない。
心臓の外科手術において弁膜症が関与する割合は大変高く、７０％以上の手術において弁膜が起因しているというデータも存在する。心臓には、図３に示すように、肺動脈弁、三尖弁、僧帽弁、大動脈弁からなる４種類の弁が存在し、それらの弁の中で、弁の動きの低下による疾患の大きな原因に該当し得るものが、三尖弁である。三尖弁は、一種の腱であり、ＣＴやＭＲＩにより撮影された画像では鮮明にその状態を捕らえることが困難であった。そのため、ＭＲＩ画像などから三尖弁の疾患の様相を視覚的に捕らえることが難しく、疾患が生じているか否かの判断基準となる三尖弁の動きを観察することはさらに困難なものとなっている。
本実施の形態に係るＭＲＩ装置１では、ＭＲＩ画像において鮮明に捕らえることが困難な三尖弁の３次元的な動画像を生成する処理を行う。
図４は、ＭＲＩ装置１のＣＰＵ２０における処理内容を示したフローチャートである。ＣＰＵ２０は、心臓における特定の位置において複数のＭＲＩ画像を撮影し（例えば、図５（ａ）のような画像）、各ＭＲＩ画像において、三尖弁の位置を自動的に特定する処理を行う（図４のステップＳ．１１〜ステップＳ．２１）。ＣＰＵ２０は、三尖弁の位置検出処理を、図５（ｂ）に示すように、少しずつ位置（図５（ｂ）における撮影位置Ｌ）をずらして数回測定を行う（撮影位置Ｌを図５（ｂ）に示すようにＬ＝１からＬ＝６までずらして測定する、図４のステップＳ２２〜Ｓ２４）。本実施の形態に係るＭＲＩ装置１では、心臓の上下方向に１２ｍｍずつ検出位置をずらして合計６箇所の位置においてＭＲＩ画像の撮影を行う。
そして、ＣＰＵ２０は、測定された全てのＭＲＩ画像より三尖弁の位置を検出し、検出された三尖弁の位置に基づいて、三尖弁の状態を３次元的に示した動画像を生成する（図４のステップＳ．２５〜２８）。これらの詳細な処理内容を、以下に説明する。
まず、ＣＰＵ２０は、記憶部２３に記録された初期情報を読み出して初期設定を行う（ステップＳ．１１）。具体的には、ＭＲＩ画像の撮影タイミング（撮影時間）を示すｔの値を１に設定し（ｔ＝１）、ＭＲＩ画像を撮影する位置を示すＬの値を１に設定する（撮影位置Ｌ＝１）処理を行う。この初期設定値は、ＲＡＭ２２に一時的に記録され、ＣＰＵ２０における処理に応じて適宜更新されることになる。
次に、ＣＰＵ２０が、シーケンサ７に対して制御指示を出力することにより、ＭＲＩ装置１の磁場発生系統や高周波送受信系統を動作させて、被検者Ｘの心臓の２次元画像（心臓が撮影範囲に含まれる２次元画像）をＭＲＩ装置１によって撮影する（ステップＳ．１２）。２次元画像をＭＲＩ装置１により撮影する場合には、心臓の同じ位置（同じ撮影位置）における複数のＭＲＩ画像を一定の時間間隔で複数枚撮影する。
撮影された画像は、記憶部２３に一時的に記録される。撮影されたＭＲＩ画像は、最初に撮影されたＭＲＩ画像を、時間ｔ＝１において撮影されたＭＲＩ画像とし、その後に撮影されたＭＲＩ画像を、時間ｔ＝２のＭＲＩ画像、時間ｔ＝３のＭＲＩ画像、・・・時間ｔ＝ＴにおけるＭＲＩ画像として区別するものとする。また、本実施の形態においては、１つの鼓動（ビート）が開始された時をｔ＝１とし、次の鼓動が開始される直前がｔ＝Ｔとなるようにする。このようにして、１つの鼓動間隔で、ｔ＝１からｔ＝ＴまでのＴ枚のＭＲＩ画像を撮影することにより、１つの鼓動の開始から終了までの心臓の状態を撮影することができる。
次に、ＣＰＵ２０は、記憶部２３に記録された２次元のＭＲＩ画像のそれぞれに対してラプラシアンと呼ばれる画像マスク処理を行う（ステップＳ．１３）。ＭＲＩ画像に画像マスク処理を施すことによって、ＭＲＩ画像のエッジを強調することができ、画像を鮮明にすることが可能となる。なお、ラプラシアンによる画像マスク処理は、既に画像処理技術として一般的に知られた方法であるため、詳細なマスク処理の内容については説明を省略する。
そして、ＣＰＵ２０は、時間ｔ＝１のＭＲＩ画像に対して、操作入力部１６を介して医師等が三尖弁の位置をマークすることを可能な状態とする（医師等によるマークの許可処理）（ステップＳ．１４）。ＭＲＩ画像には、腱である三尖弁が鮮明には表示されにくい傾向がある。しかしながら、医師等は、ＭＲＩ画像において鮮明に表示される部位（例えば、血液が多い部位）などの位置関係から、三尖弁の位置を推測することが可能である。このため、ＣＰＵ２０は、最初に撮影されたＭＲＩ画像（ｔ＝１において撮影されたＭＲＩ画像）に対して、医師等が三尖弁の位置をマークすることができるようにする（マーク処理を医師等に行わせる）ことにより、ＭＲＩ画像における三尖弁の位置（マーク位置）を高い精度で特定することが可能となる。
図６（ａ）は、ｔ＝１のＭＲＩ画像に対して、マーク処理が行われる前の画像を示し、図６（ｂ）は、マーク処理が行われた後の画像を示している。三尖弁は円環状（リング状）の腱であるため、円環状の端部であれば、ＭＲＩ画像においてマークされる箇所が１箇所となり、円環状の中央部（少なくとも端部以外）であれば、２箇所となる。図６（ｂ）においては、２箇所のマークＭが示されている。
医師等は、操作入力部１６の操作により三尖弁の位置をＭＲＩ画像にマークし終わった場合には、マークを終了した旨の操作（例えば、操作入力部１６において予め設けられるマーク終了操作ボタンの押下など）を行う。マークを終了した旨の操作が行われていない場合（ステップＳ．１５においてＮｏの場合）、ＣＰＵ２０は、三尖弁の位置のマーク処理が行われている最中であると判断して、処理をステップＳ．１４の処理へ移行する。一方で、マークを終了した旨の操作が行われた場合（ステップＳ．１５においてＹｅｓの場合）、ＣＰＵ２０は、三尖弁の位置のマーク処理が終了したものと判断する。
三尖弁の位置のマーク処理が終了した場合（ステップＳ．１５においてＹｅｓの場合）、ＣＰＵ２０は、２次元連続ＤＰ(weak spotting version)と呼ばれるアルゴリズムに基づいて、三尖弁の位置がマークされた時間ｔのＭＲＩ画像を参照画像とし、時間ｔ＋１のＭＲＩ画像を入力画像とすることにより、参照画像および入力画像における全ピクセルの対応関係を求める（ステップＳ．１６）。ここで、２次元連続ＤＰ（Dynamic Programming）は、既に知られた技術であり、『岡隆一，外２名「連続ＤＰの一般スキームについて−画像スポッティングのための全画素最適マッチング−」，信学技報，電子情報通信学会，IEICE Technical Report, PRMU2010-87,IBISML2010-59 (2010-09)』および特開２０１０−１６５１０４号公報おいて詳細な説明がなされているため、具体的な説明を省略する。

そして、ＣＰＵ２０は、２次元連続ＤＰに基づく全ピクセルの対応関係から、時間ｔのＭＲＩ画像においてマークした三尖弁の全ピクセルに対応するピクセルを、時間ｔ＋１のＭＲＩ画像において求める（ステップＳ．１７）。
次に、ＣＰＵ２０は、２次元連続ＤＰによりマークの対応関係が求められた時間ｔ＋１のＭＲＩ画像のピクセル群に対して、ダイレーション（dilation）というピクセルの増大処理とエロージョン（erosion）というピクセル削減処理とを施す（ステップＳ．１８）。時間ｔ＋１のＭＲＩ画像に対して、このダイレーションおよびエロージョンを施すことにより、時間ｔ＋１において三尖弁に対応するピクセルとして求められたピクセル群の孤立点（孤立したピクセル）の除去、不連続なピクセルの接続と穴埋めが行われ、時間ｔ＋１のＭＲＩ画像のマーク（マーク位置におけるピクセル群）が明確になる。
本実施の形態に係るＭＲＩ装置１では、一例として、ダイレーション（dilation：ピクセルの増大処理）を２回行い、エロージョン（erosion：ピクセル削減処理）を４回行う。ただし、ダイレーション（dilation）とエロージョン（erosion）とのそれぞれの処理回数は、特に限定されるものではない。
なお、ダイレーション（dilation）とエロージョン（erosion）とによるピクセルの増大・削減処理は、モルフォロジ処理と呼ばれる画像処理の代表的な処理であり、上述したように、孤立点（孤立したピクセル）の除去、不連続なピクセルの接続と穴埋めのために利用される処理である。この画像処理は既に知られた処理方法であるため、ここでの詳細な説明は省略する。
図７は、ｔ＝１のＭＲＩ画像においてマークされた三尖弁のピクセル群から２次元連続ＤＰに基づいて求められたｔ＝２のＭＲＩ画像のマークを、ダイレーション（dilation）とエロージョン（erosion）とをそれぞれ複数回繰り返すことにより明確化する一例を示したＭＲＩ画像である。
まず、図７（ａ）は、既に説明したように、ｔ＝１のＭＲＩ画像であって、既に三尖弁の位置にマークＭ１がなされた画像である。一方で図７（ｂ）に示すＭＲＩ画像は、２次元連続ＤＰに基づいて対応するマークが求められたｔ＝２のＭＲＩ画像を示している。図７（ｂ）の拡大画像に示すように、２次元連続ＤＰに基づいて求められたｔ＝２のＭＲＩ画像のマークＭ２は、（ａ）のマークＭ１と異なり、細かい点の集まりとして示されている。
図７（ｂ）に示すようにして求められたＭＲＩ画像に対して、ＣＰＵ２０がダイレーション（dilation）とエロージョン（erosion）とを複数回実行する。ダイレーションとエロージョンとを実行することにより、図７（ｃ）に示すように、マークＭ３が一塊のピクセル群となり、マークＭ３の明確化を実現することが可能となる。
その後、ＣＰＵ２０は、マークの明確化が図られたＭＲＩ画像を、記憶部２３に記録する（ステップＳ．１９）。なお、このＭＲＩ画像は、撮影位置Ｌの値と、撮影時間ｔの値とを関連付けて記録される。
次に、ＣＰＵ２０は、時間ｔ＋１のＭＲＩ画像が時間ＴのＭＲＩ画像であるか否か（つまり、ｔ＋１＝Ｔであるか否か、ｔ＋１の画像が各撮影位置における最後のＭＲＩ画像であるか否か）を判断する（ステップＳ．２０）。時間ｔ＋１のＭＲＩ画像が時間ＴのＭＲＩ画像でない場合（ステップＳ．２０においてＮｏの場合）、ＣＰＵ２０は、ｔ＝ｔ＋１の処理によって、ｔの値を１つ増加させて（ステップＳ．２１）、処理を上述したステップＳ．１６へと移行させて、上述した処理を繰り返し実行する。
一方で、時間ｔ＋１のＭＲＩ画像が時間ＴのＭＲＩ画像である場合（ステップＳ．２０においてＹｅｓの場合）、ＣＰＵ２０は、ＭＲＩ画像の撮影位置Ｌが６箇所目の撮影位置であるか否か、つまりＬ＝６であるか否かの判断を行う（ステップＳ．２２）。
ＭＲＩ画像の撮影位置Ｌが６箇所目でない場合（ステップＳ．２２においてＮｏの場合）、ＣＰＵ２０は、シーケンサ７に対して制御指示を出力することにより、ＭＲＩ装置１におけるＭＲＩ画像の撮影位置を１２ｍｍだけ移動させる（ステップＳ．２３）。そして、ＣＰＵ２０は、撮影位置Ｌの値を１つ増加させて（ステップＳ．２４）、処理をステップＳ．１２へと移行させる。
ここで、撮影位置が移動する方向をＺ軸方向とし、ＭＲＩ画像の撮影が行われる心臓の断面（スライス面）をＸＹ平面とする。ＭＲＩ画像は、一般的に５１２ピクセル×５１２ピクセルの解像度で撮影され、この平面は、物理的に３５ｃｍ×３５ｃｍの大きさに対応するものである。従って、１ピクセルの区間は、０．６８４ｍｍに該当することになる。
５１２ピクセル×５１２ピクセルのＭＲＩ画像の中で、心臓が映っている領域を切り出すと、２５０ピクセル×２５０ピクセルの大きさとなり、三尖弁の位置検出に利用されるＭＲＩ画像の大きさも２５０ピクセル×２５０ピクセルの大きさに含まれることになる。従って、切り出された領域の画像の一辺の実際の長さは、２５０ピクセル×０．６８４ｍｍ＝１７０．９ｍｍとなる。
このことから、三尖弁の位置検出において求められる３次元空間は、Ｘ軸方向およびＹ軸方向の長さが１７０．９ｍｍ、Ｚ軸報告の長さが、１２ｍｍ×５（撮影位置Ｌ＝１からＬ＝６までの距離）＝６０ｍｍとなる。従って、三尖弁の位置検出が行われる物理的な空間は、１７１ｍｍ×１７１ｍｍ×６０ｍｍの空間と相似したものとなる。なお、Ｚ軸方向の長さ６０ｍｍは、ピクセル換算すると、６０ｍｍ÷０．６８４ｍｍ（１ピクセルの長さ）＝８８となり、８８ピクセルに対応する。
ＭＲＩ画像の撮影位置Ｌが６箇所目である場合（ステップＳ．２２においてＹｅｓの場合）、つまり、最後の撮影位置（Ｌ＝６）におけるＴ枚目のＭＲＩ画像のマーク処理が終わった場合、ＣＰＵ２０は、撮影された全てのＭＲＩ画像のマークから、マーク領域を構成するピクセル群の重心位置を算出する（ステップＳ．２５）。このように各ＭＲＩ画像におけるマークの重心位置を求めることにより、心臓の鼓動に伴って動く三尖弁の位置を時系列的に抽出することが可能となる。
また、ＭＲＩ画像による三尖弁のスライス画像は、上述した処理から明らかなように、一度に６箇所同時に撮影することができない。しかしながら、各撮影位置におけるＭＲＩ画像の撮影は、心臓の鼓動が開始されたタイミングで時間ｔ＝１のＭＲＩ画像が撮影され、次の鼓動が開始されるまでの間に時間ｔ＝ＴのＭＲＩ画像の撮影が行われる。この撮影タイミングは、撮影位置Ｌ＝１からＬ＝６まで同じタイミングであるため、鼓動の動きに対応してそれぞれの撮影位置で同じタイミングのＭＲＩ画像を撮影したものと判断することができる。従って、鼓動が開始したタイミングである時間ｔ＝１のＭＲＩ画像を初期時刻として同一視すれば、複数の撮影位置において撮影されたＭＲＩ画像を、同時期に撮影されたＭＲＩ画像として処理することが可能となり、三尖弁の３次元的な動きを求めることが可能となる。
なお、三尖弁はリング状を呈するため、三尖弁のＭＲＩ画像（スライス画像）には、三尖弁のリング状断面が１断面あるいは２断面撮影されることになる。このリング状断面は、撮影位置によって異なり、三尖弁の端部側を撮影する場合には、三尖弁の断面が１断面となる。従って、撮影位置の異なる６箇所の撮影箇所のＭＲＩ画像から、時間ｔに応じて１０〜１２個の三尖弁のリング位置（マークの重心位置）を３次元的に求めることができる。３次元的に示される１０〜１２個の三尖弁のリング位置を、時間ｔに応じて、ｔ＝１からｔ＝Ｔまで変化させることにより、３次元的に求められる三尖弁の状態を時系列に求めることが可能となり、三尖弁の動きを求めることができる。
しかしながら、上述のようにして求められる三尖弁のリング位置は、時間ｔにおいて１０〜１２個の断片的な位置情報であって、３次元における全体的な三尖弁の位置情報の一部でしかない。一方で、多くの位置情報を得るためには多数の撮影位置においてＭＲＩ装置１で三尖弁のスライス撮影を行うことが好ましいが、現実的には困難である。
このため、求められた時間ｔ毎の三尖弁のリング位置を三尖弁の３次元的なサンプル点とし、３次元空間における三尖弁のリング形状が、サンプル点を繋ぐことにより形成される滑らかな３次元上の曲線によるものと想定することにより、合理性を備えた三尖弁のリング形状を求めることができる。
従って、ＣＰＵ２０は、時間ｔ毎に、それぞれのＭＲＩ画像から求められた１０〜１２個の三尖弁のリング位置（マークの重心位置）を通り滑らかな閉曲線を形成するスプライン曲線を算出する（ステップＳ．２６）。いくつかのサンプル点に基づいてスプライン曲線を算出する手法は既に知られているため、ここでの詳細な説明は省略する。このようにして３次元的なスプライン曲線を求めることによって、３次元的な三尖弁のリング形状を、例えば１０００点の空間点のつながりによって表現することが可能となる。
図８（ａ）は、サンプル点（白い球体）を通るスプライン曲線状の各点に、小さな球体を配置することにより三尖弁のリング形状を表現した図である。心臓の１つの鼓動（ビート）の開始から次の鼓動の開始直前までの全体の時間を、前述した既出の時間Ｔとして考えると、Ｔ個のスプライン曲線を得ることが可能となる。従って、時間ｔ＝１からＴまでのＴ個のスプライン曲線に基づいて、１つの鼓動の開始から終了までの時間における三尖弁の３次元的な動きを求めることができる。
ＣＰＵ２０では、三尖弁の３次元的な座標情報を、時間ｔ＝１からＴまでのＴ個のスプライン曲線に基づいて求めて（ステップＳ．２７）、求められた三尖弁の３時次元的な座標情報の変化に基づいて、１つ鼓動に対応する三尖弁の３次元動画像を生成する（ステップＳ．２８）。
このようにして求められた三尖弁の３次元動画を繰り返し表示部１５に表示させる（再生表示させる）ことにより、ＭＲＩ装置１において鮮明に撮影できない腱などの部位の状態変化、本実施の形態では三尖弁の状態変化を明確に示した映像として抽出し、医師などに提供することが可能となる。このため、医師は、三尖弁の鮮明な３次元画像および３次元動画像に基づいて、リング形状に異常がないか否か、また、動きに異変がないか否かを、視覚的に判断することが可能となる。
さらに、三尖弁の動画像は、三尖弁のリングの全て位置情報とその位置変化（動き）が、３次元的な電子データとして記録されている。このため、３次元空間を小さなサイコロ状に分割し、三尖弁の動画像に基づいて、１つの鼓動において三尖弁の移動の際に接触するサイコロ状の分割部分の個数を加算することにより、三尖弁の可動体積を算出することができる。サイコロ状の分割間隔は、既に説明したようにミリメートル（ｍｍ）単位で可能であるため、立法ミリメートル単位での可動体積を求めることが可能となる。
図８（ｂ）は、三尖弁の可動体積を示した３次元的な図である。正常な三尖弁の可動体積と異常がある三尖弁の可動体積とを、信頼できるサンプルに基づいて多数収集することにより、統計的な検定によって疾患の有無を判定することも可能である。
以上、本発明に係る体内臓器の動画像生成装置および体内臓器の動画像生成方法について、ＭＲＩ装置１を一例として用い、ＭＲＩ装置１で三尖弁を撮影する場合について説明を行ったが、本発明に係る体内臓器の動画像生成装置および体内臓器の動画像生成方法は、上述した実施の形態に示したＭＲＩ装置１の構成およびその処理方法だけに限定されず、また、三尖弁を撮影する場合だけに限定されるものではない。
例えば、実施の形態に係るＭＲＩ装置１では、心臓の三尖弁のリング形状の動画像を生成する場合について説明したが、他の部位を時間ｔ＝１のＭＲＩ画像にマークすることにより、マーク位置に該当する部位の動画像を生成することが可能である。従って、例えば、心臓の心房などの部位にマークを設定することにより心房における動画像を生成することができ、心房の動作状況の可視化や、可動体積の算出などを同様にして行うことが可能である。
特に、本発明に係る体内臓器の動画像生成装置では、ＭＲＩ画像において必ずしも鮮明に示されていない部位であっても、時間ｔ＝１のＭＲＩ画像においてマークを設定することによって、２次元連続ＤＰにより対応関係が求められたピクセル群を自動的に求めることができる。このため、時間ｔ＝１のＭＲＩ画像において部位のマーク処理を確実に行うことにより、該当するマーク位置の部位の状態を自動的に求めることができる。従って、従来では視覚的に確認することができなかった部位の形状（様子）およびその動作状況を可視化することができ、臨床診断における診断精度の向上とその診断負担の軽減を図ることが可能となる。
また、本実施の形態に係るＭＲＩ装置１では、断面画像の画像マスク処理や、２次元連続ＤＰに基づくマークの特定処理や、ダイレーションおよびエロージョンによる画像処理や、３次元的なスプライン曲線を求める処理などを、ＣＰＵ２０だけで行う場合について説明した。しかしながら、これらの全ての処理を、１つのＣＰＵだけに行わせる構成には限定されない。各処理に応じてそれぞれ異なるＣＰＵを設けて、複数のＣＰＵが順次処理を行うことによって、結果として体内臓器の３次元的な動画像を生成する構成であってもよい。
１ …ＭＲＩ装置
２ …静磁場発生磁石（断面画像撮影手段）
３ …高周波送信用コイル（断面画像撮影手段）
４ …高周波受信用コイル（断面画像撮影手段）
５ …傾斜磁場コイル（断面画像撮影手段）
６ …コンピュータ
７ …シーケンサ（断面画像撮影手段）
８ …高周波発生装置（断面画像撮影手段）
９ …変調器（断面画像撮影手段）
１０ …増幅器（断面画像撮影手段）
１１ …傾斜磁場発生装置（断面画像撮影手段）
１２ …増幅器（断面画像撮影手段）
１３ …位相検波器（断面画像撮影手段）
１４ …ＡＤ変換器（断面画像撮影手段）
１５ …表示部
１６ …操作入力部（マーク付加手段）
２０ …ＣＰＵ（撮影制御手段、画像解析手段）
２１ …ＲＯＭ
２２ …ＲＡＭ
２３ …記憶部（記録手段）
Ｘ …被検者

Claims (6)

  1. 体内臓器における断面画像の撮影位置を当該断面画像の積層方向へ移動させることができ、かつ、前記断面画像を所定時間間隔で連続的に撮影することが可能な断面画像撮影手段と、
    該断面画像撮影手段により撮影された前記断面画像を、同一撮影位置において連続して撮影された前記断面画像の撮影順番と、当該断面画像の撮影位置とに関連付けて記録する記録手段と、
    前記断面画像撮影手段における前記撮影位置および撮影タイミングを制御する撮影制御手段と、
    前記記録手段に記録された前記断面画像の画像解析を行う画像解析手段と、
    前記断面画像に対してユーザがマーク位置を設定することにより、前記体内臓器における所定部位の位置をマークによって当該断面画像に記録させることが可能なマーク付加手段と
    を有し、
    前記撮影制御手段は、前記断面画像撮影手段に、前記体内臓器の前記断面画像を一の撮影位置において連続的に撮影させ、かつ、当該連続的な前記断面画像の撮影を、前記体内臓器が収まる積層方向範囲において少しずつ撮影位置をずらして当該撮影位置毎に行わせ、
    前記画像解析手段は、
    前記記録手段に記録された断面画像を読み出して、
    一の撮影位置において連続的に撮影された複数の前記断面画像のうち当該撮影位置において最初に撮影された断面画像に対して、前記マーク付加手段により設定されたマーク位置に前記マークを付加し、
    ２次元連続ＤＰを用いて、前記マークが付加された前記断面画像のマーク位置に基づいて、前記マークが付加された前記断面画像の次に撮影された断面画像にマークを付加することにより、連続的に撮影された全ての撮影位置における全ての断面画像に自動的にマークを付加し、
    前記断面画像に付加された前記マークの重心位置を、当該マークのピクセル群に基づいて全ての断面画像について算出し、
    算出された重心位置に基づいて、前記断面画像の画像平面および前記積層方向からなる３次元空間における前記体内臓器の前記部位の座標情報を、異なる撮影位置であって前記連続的に撮影された断面画像のうち同じ撮影順番で撮影された断面画像毎に求め、
    同じ撮影順番で撮影された前記座標情報に基づいて３次元的なスプライン曲線を求めて、撮影順番毎に前記部位の３次元的な形状を求め、
    前記撮影順番毎に求められた前記部位の３次元的な形状を、前記撮影順番に従って変化させることにより、前記部位の３次元的な形状の時間的な変化を求めて、前記部位における３次元的な動画像を生成する
    ことを特徴とする体内臓器の動画像生成装置。
  2. 前記画像解析手段は、前記マークの重心位置を全ての断面画像について算出する前に、当該マークが付加された前記断面画像に対して、ダイレーションとエロージョンとによる画像処理を施すこと
    を特徴とする請求項１に記載の体内臓器の動画像生成装置。
  3. 前記撮影制御手段は、前記断面画像撮影手段に前記断面画像を連続的に撮影させる場合に、前記体内臓器を備える被検者の一の鼓動の開始から次の鼓動が開始されるまでの断面画像を連続的に撮影させること
    を特徴とする請求項１または請求項２に記載の体内臓器の動画像生成装置。
  4. 体内臓器における断面画像の撮影位置を当該断面画像の積層方向へ移動させることができ、かつ、前記断面画像を所定時間間隔で連続的に撮影することが可能な断面画像撮影手段と、
    該断面画像撮影手段により撮影された前記断面画像を、同一撮影位置において連続して撮影された前記断面画像の撮影順番と、当該断面画像の撮影位置とに関連付けて記録する記録手段と、
    前記断面画像撮影手段における前記撮影位置および撮影タイミングを制御する撮影制御手段と、
    前記記録手段に記録された前記断面画像の画像解析を行う画像解析手段と、
    前記断面画像に対してユーザがマーク位置を設定することにより、前記体内臓器における所定部位の位置をマークによって当該断面画像に記録させることが可能なマーク付加手段と
    を有する体内臓器の動画像生成装置において、
    前記撮影制御手段が、
    前記断面画像撮影手段に、前記体内臓器の前記断面画像を一の撮影位置において連続的に撮影させ、かつ、当該連続的な前記断面画像の撮影を、前記体内臓器が収まる積層方向範囲において少しずつ撮影位置をずらして当該撮影位置毎に行わせ、
    前記画像解析手段が、
    前記記録手段に記録された断面画像を読み出して、
    一の撮影位置において連続的に撮影された複数の前記断面画像のうち当該撮影位置において最初に撮影された断面画像に対して、前記マーク付加手段により設定されたマーク位置に前記マークを付加し、
    ２次元連続ＤＰを用いて、前記マークが付加された前記断面画像のマーク位置に基づいて、前記マークが付加された前記断面画像の次に撮影された断面画像にマークを付加することにより、連続的に撮影された全ての撮影位置における全ての断面画像に自動的にマークを付加し、
    前記断面画像に付加された前記マークの重心位置を、当該マークのピクセル群に基づいて全ての断面画像について算出し、
    算出された重心位置に基づいて、前記断面画像の画像平面および前記積層方向からなる３次元空間における前記体内臓器の前記部位の座標情報を、異なる撮影位置であって前記連続的に撮影された断面画像のうち同じ撮影順番で撮影された断面画像毎に求め、
    同じ撮影順番で撮影された前記座標情報に基づいて３次元的なスプライン曲線を求めて、撮影順番毎に前記部位の３次元的な形状を求め、
    前記撮影順番毎に求められた前記部位の３次元的な形状を、前記撮影順番に従って変化させることにより、前記部位の３次元的な形状の時間的な変化を求めて、前記部位における３次元的な動画像を生成する
    ことを特徴とする体内臓器の動画像生成方法。
  5. 前記画像解析手段が、
    前記マークの重心位置を全ての断面画像について算出する前に、当該マークが付加された前記断面画像に対して、ダイレーションとエロージョンとによる画像処理を施すこと
    を特徴とする請求項４に記載の体内臓器の動画像生成方法。
  6. 前記撮影制御手段が、
    前記断面画像撮影手段に前記断面画像を連続的に撮影させる場合に、前記体内臓器を備える被検者の一の鼓動の開始から次の鼓動が開始されるまでの断面画像を連続的に撮影させること
    を特徴とする請求項４または請求項５に記載の体内臓器の動画像生成方法。
JP2011239897A 2011-11-01 2011-11-01 体内臓器の動画像生成装置および体内臓器の動画像生成方法 Expired - Fee Related JP5879098B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2011239897A JP5879098B2 (ja) 2011-11-01 2011-11-01 体内臓器の動画像生成装置および体内臓器の動画像生成方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2011239897A JP5879098B2 (ja) 2011-11-01 2011-11-01 体内臓器の動画像生成装置および体内臓器の動画像生成方法

Publications (2)

Publication Number Publication Date
JP2013097586A JP2013097586A (ja) 2013-05-20
JP5879098B2 true JP5879098B2 (ja) 2016-03-08

Family

ID=48619464

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2011239897A Expired - Fee Related JP5879098B2 (ja) 2011-11-01 2011-11-01 体内臓器の動画像生成装置および体内臓器の動画像生成方法

Country Status (1)

Country Link
JP (1) JP5879098B2 (ja)

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001155019A (ja) * 1999-11-25 2001-06-08 Olympus Optical Co Ltd 類似画像検索装置
JP4630564B2 (ja) * 2004-03-30 2011-02-09 国立大学法人浜松医科大学 手術支援装置、方法及びプログラム
JP5526401B2 (ja) * 2009-03-31 2014-06-18 国立大学法人東京農工大学 心室壁情報抽出装置
JP5068336B2 (ja) * 2010-03-19 2012-11-07 富士フイルム株式会社 医用画像変換装置および方法並びにプログラム

Also Published As

Publication number Publication date
JP2013097586A (ja) 2013-05-20

Similar Documents

Publication Publication Date Title
JP6238651B2 (ja) 超音波診断装置及び画像処理方法
JP4884528B2 (ja) 腔部の画像を評定する方法、装置ならびにコンピュータプログラム製品
JP5031758B2 (ja) 医用画像診断装置
US8213699B2 (en) Medical image processing apparatus and medical image diagnosis apparatus
Li et al. Quantification and MRI validation of regional contractile dysfunction in mice post myocardial infarction using high resolution ultrasound
JP5105981B2 (ja) 医用画像処理表示装置、およびその処理プログラム
CN102665560B (zh) X射线ct装置以及基于x射线ct装置的图像显示方法
JP5580030B2 (ja) 画像処理装置、および画像位置合せ方法
US20200327979A1 (en) Medical information processing apparatus and medical information processing method
CN102232845B (zh) 用ct系统自动识别患者血管中造影剂流动的方法及ct系统
CN101243971B (zh) 用于确定心脏兴奋传导的方法以及对应的医学装置
JP2012205899A (ja) ３次元的モデルを利用した身体臓器の映像生成方法及び装置並びにコンピュータ読み取り可能な記録媒体
JP2008086767A (ja) ３次元及び４次元コントラスト撮像のためのシステム及び方法
CN104080408A (zh) 超声波诊断装置、图像处理装置以及图像处理方法
JP6739318B2 (ja) 超音波診断装置
JP2007167656A (ja) 対象の運動の解析方法および断層撮影装置
WO2013125276A1 (ja) Ｘ線ｃｔ装置、画像表示装置、画像表示方法
JP2010094181A (ja) 超音波診断装置および超音波診断装置のデータ処理プログラム
JP2004121834A (ja) 生体組織の動き追跡方法、その追跡方法を用いた画像診断装置及び生体組織の動き追跡プログラム
JP2004313551A (ja) 心臓の運動表示方法及び画像診断装置
JP2010154982A (ja) Ｘ線コンピュータ断層撮影装置および画像処理装置
JP6863774B2 (ja) 超音波診断装置、画像処理装置及び画像処理プログラム
CN111317508B (zh) 超声波诊断装置、医用信息处理装置、计算机程序产品
JP5885234B2 (ja) 疾患判定装置および超音波画像形成装置における画像解析方法
JP2004313545A (ja) 心臓の運動表示方法及び画像診断装置

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20141001

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20141001

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20150803

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20150818

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20151005

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: 20160119

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20160201

R150 Certificate of patent or registration of utility model

Ref document number: 5879098

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

LAPS Cancellation because of no payment of annual fees