(57)【特許請求の範囲】
【請求項１】
時系列に並ぶ複数のフェーズにおける心臓の領域を含み、ＣＴ（Computed Tomography）装置により得られたボリュームデータを取得する画像取得部と、
各フェーズの前記ボリュームデータにおいて、所定の断面で切断された断面画像を導出する断面画像導出部と、
各フェーズの前記断面画像において、心室心筋の輪郭を設定する輪郭設定部と、
各フェーズの前記輪郭に基づいて、心周期における所定期間又は所定時点における前記輪郭の長さの変化を示す変化情報を導出する変化情報導出部と、
前記変化情報を表示する表示部と、
を備えるパラメータ処理装置。
【請求項２】
請求項１に記載のパラメータ処理装置であって、
前記輪郭は、左室心筋の輪郭、左室の下壁における心基部から心尖部までの左室の外壁に沿った輪郭、及び左室の側壁における心尖部から心基部までの左室の外壁に沿った輪郭、のうち少なくとも１つである、パラメータ処理装置。
【請求項３】
請求項１または２に記載のパラメータ処理装置であって、
前記変化情報導出部は、各フェーズの前記輪郭に基づいて、心周期における収縮期を含む収縮期近傍の所定期間又は所定時点における前記輪郭の長さの変化加速度を導出し、
前記表示部は、前記変化加速度の情報を表示する、を備えるパラメータ処理装置。
【請求項４】
請求項１または２に記載のパラメータ処理装置であって、
前記変化情報導出部は、各フェーズの前記輪郭に基づいて、心周期における拡張期から収縮期までの間の所定期間又は所定時点における前記輪郭の長さの変化速度を導出し、
前記表示部は、前記変化速度の情報を表示する、パラメータ処理装置。
【請求項５】
請求項１ないし４のいずれか１項に記載のパラメータ処理装置であって、
前記所定の断面は、左室短軸断面又は左室長軸断面を含む、パラメータ処理装置。
【請求項６】
請求項１ないし５のいずれか１項に記載のパラメータ処理装置であって、
前記所定の断面は複数存在し、
前記輪郭設定部は、各フェーズの複数の断面画像において、前記輪郭を設定し、
前記変化情報導出部は、各フェーズの複数の前記輪郭に基づいて、複数の前記変化情報を導出し、
前記表示部は、複数の前記変化情報を表示する、パラメータ処理装置。
【請求項７】
請求項１ないし６のいずれか１項に記載のパラメータ処理装置であって、更に、
各フェーズの断面画像における前記輪郭の位置を調整するレジストレーション処理部を備える、パラメータ処理装置。
【請求項８】
請求項１ないし７のいずれか１項に記載のパラメータ処理装置であって、更に、
連続するフェーズの間においてフェーズを補間して前記ボリュームデータを導出するフェーズ補間処理部を備え、
前記断面画像導出部は、前記画像取得部により取得されたフェーズ及び前記フェーズ補間処理部により補間されたフェーズのボリュームデータにおいて、前記断面画像を導出する、パラメータ処理装置。
【請求項９】
請求項１ないし８のいずれか１項に記載のパラメータ処理装置であって、
前記表示部は、前記輪郭設定部により設定された前記輪郭の長さの変化を示すグラフを表示する、パラメータ処理装置。
【請求項１０】
パラメータ処理装置におけるパラメータ処理方法であって、
時系列に並ぶ複数のフェーズにおける心臓の領域を含み、ＣＴ装置により得られたボリュームデータを取得するステップと、
各フェーズの前記ボリュームデータにおいて、所定の断面で切断された断面画像を導出するステップと、
各フェーズの前記断面画像において、心室心筋の輪郭を設定するステップと、
各フェーズの前記輪郭に基づいて、心周期における所定期間又は所定時点における前記輪郭の長さの変化を示す変化情報を導出するステップと、
前記変化情報を表示部に表示するステップと、
を備えるパラメータ処理方法。
【請求項１１】
請求項１０に記載のパラメータ処理方法であって、
前記輪郭は、左室心筋の輪郭、左室の下壁における心基部から心尖部までの左室の外壁に沿った輪郭、及び左室の側壁における心尖部から心基部までの左室の外壁に沿った輪郭、のうち少なくとも１つである、パラメータ処理方法。
【請求項１２】
時系列に並ぶ複数のフェーズにおける心臓の領域を含み、ＣＴ装置により得られたボリュームデータを取得するステップと、
各フェーズの前記ボリュームデータにおいて、所定の断面で切断された断面画像を導出するステップと、
各フェーズの前記断面画像において、心室心筋の輪郭を設定するステップと、
各フェーズの前記輪郭に基づいて、心周期における所定期間又は所定時点における前記輪郭の長さの変化を示す変化情報を導出するステップと、
前記変化情報を表示部に表示するステップと、
をコンピュータに実行させるためのパラメータ処理プログラム。
【請求項１３】
請求項１２に記載のパラメータ処理プログラムであって、
前記輪郭は、左室心筋の輪郭、左室の下壁における心基部から心尖部までの左室の外壁に沿った輪郭、及び左室の側壁における心尖部から心基部までの左室の外壁に沿った輪郭、のうち少なくとも１つである、パラメータ処理プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、パラメータ処理装置、パラメータ処理方法、及びパラメータ処理プログラムに関する。
【背景技術】
【０００２】
心臓疾患の１つとして、心臓弁膜症が知られている。心臓弁膜症として、大動脈弁狭窄症（ＡＳ：Ａｏｒｔｉｃ  ｖａｌｖｅ  ｓｔｅｎｏｓｉｓ）及び大動脈弁閉鎖不全（ＡＲ：Ａｏｒｔｉｃ  ｒｅｇｕｒｇｉｔａｔｉｏｎ）が知られている。図１４（Ａ）は、大動脈弁狭窄症と診断される心臓の左室短軸断面図である。図１４（Ｂ）は、大動脈弁閉鎖不全と診断される心臓の左室短軸断面図である。
【０００３】
大動脈弁狭窄症は、大動脈弁が十分に開かなくなり、左室から血液が駆出され難くなる疾患である。大動脈弁閉鎖不全は、大動脈弁が十分に閉まらなくなり、左室から駆出された血液が逆流し易くなる疾患である。大動脈弁狭窄症、大動脈弁閉鎖不全と診断された場合、例えば大動脈弁を人工弁に置換する処置が行われる。図１５では、矢印先端部の位置が、心臓における大動脈弁の位置に相当する。
【０００４】
心臓弁膜症を診断する方法として、例えば、心エコー検査、心臓カテーテル検査、及び心臓ＣＴ（Ｃｏｍｐｕｔｅｄ  Ｔｏｍｏｇｒａｐｈｙ）検査が知られている。
【０００５】
心エコー検査では、超音波を利用する超音波検査装置により、心臓及び大血管の像が抽出される。検査担当者は、心エコー検査により、血流を確認できる。心エコー検査の１つとして、ＴＴＥ（Ｔｒａｎｓｔｈｏｒａｃｉｃ  ｅｃｈｏｃａｒｄｉｏｇｒａｍ）がある。ＴＴＥでは、超音波接触子が心臓の各部に当てられて、心臓付近の像が抽出される。ＴＴＥは、経胸壁心エコー検査とも称される。
【０００６】
心臓カテーテル検査では、心臓に大腿部等の血管からカテーテルが挿入され、心臓カテーテル検査装置により、例えば、心臓内の各場所の血圧、心拍出量、血中の酸素含有量が測定される。心臓カテーテル検査では、必要に応じて造影剤が体内に注入される。
【０００７】
心臓ＣＴ検査では、心臓付近にＸ線が照射され、心臓ＣＴ検査装置により、体内を通過したＸ線量の差がデータとして収集されて、心臓が画像化される。心臓カテーテル検査では、必要に応じて体内に造影剤が注入される。
【０００８】
心臓ＣＴ検査として、心臓ＣＴ血管造影法（ＣＣＴＡ：Ｃａｒｄｉａｃ  Ｃｏｍｐｕｔｅｄ  Ｔｏｍｏｇｒａｐｈｙ  Ａｎｇｉｇｒａｐｈｙ）がある。ＣＣＴＡを用いることで、検査担当者は、複数のフェーズからなるボリュームデータから作られる動画から、心臓全体の動きを認識し易くなり、心臓付近の弁（大動脈弁等）の動きや形状を直接観察できる。検査担当者は、例えば、最大時の心室の大きさと最小時の心室の大きさとを視認し、左室駆出率（ＬＶＥＦ：Ｌｅｆｔ  Ｖｅｎｔｒｉｃｌｅ  Ｅｊｅｃｔｉｏｎ  Ｆｒａｃｔｉｏｎ）を確認できる。
【０００９】
また、上記ＣＴ画像のボリュームデータを用いて、大動脈弁狭窄症の治療方針を補助する医用画像処理装置が知られている（例えば、特許文献１参照）。この医用画像処理装置は、ボリュームデータの画像における大動脈弁本体の領域と石灰化領域とを別々に補間処理して合成することで、弁口面積の変化に対する石灰化部分の影響を判断する。
【先行技術文献】
【非特許文献】
【００１０】
【非特許文献１】特開２０１３−１０６８７０号公報
【発明の概要】
【発明が解決しようとする課題】
【００１１】
従来の心臓弁膜症の各検査では、検査担当者及び患者の負担が大きいことがあった。また、検査により得られるパラメータの精度が不十分であり、心臓弁膜症の診断精度が低下する可能性があった。また、検査担当者により結果がばらつくものもあった。
【００１２】
本発明は、上記事情に鑑みてなされたものであり、心臓弁膜症の検査において、検査担当者及び患者の負担を低減でき、検査により得られるパラメータの精度が向上できるパラメータ処理装置、パラメータ処理方法、及びパラメータ処理プログラムを提供する。
【課題を解決するための手段】
【００１３】
本発明のパラメータ処理装置は、時系列に並ぶ複数のフェーズにおける心臓の領域を含み、ＣＴ（Computed Tomography）装置により得られたボリュームデータを取得する画像取得部と、各フェーズの前記ボリュームデータにおいて、所定の断面で切断された断面画像を導出する断面画像導出部と、各フェーズの前記断面画像において、心室心筋の輪郭を設定する輪郭設定部と、各フェーズの前記輪郭に基づいて、心周期における所定期間又は所定時点における前記輪郭の長さの変化を示す変化情報を導出する変化情報導出部と、前記変化情報を表示する表示部と、を備える。
【００１４】
本発明のパラメータ処理方法は、パラメータ処理装置におけるパラメータ処理方法であって、時系列に並ぶ複数のフェーズにおける心臓の領域を含み、ＣＴ装置により得られたボリュームデータを取得するステップと、各フェーズの前記ボリュームデータにおいて、所定の断面で切断された断面画像を導出するステップと、各フェーズの前記断面画像において、心室心筋の輪郭を設定するステップと、各フェーズの前記輪郭に基づいて、心周期における所定期間又は所定時点における前記輪郭の長さの変化を示す変化情報を導出するステップと、前記変化情報を表示部に表示するステップと、を備える。
【００１５】
本発明のパラメータ処理方法は、時系列に並ぶ複数のフェーズにおける心臓の領域を含み、ＣＴ装置により得られたボリュームデータを取得するステップと、各フェーズの前記ボリュームデータにおいて、所定の断面で切断された断面画像を導出するステップと、各フェーズの前記断面画像において、心室心筋の輪郭を設定するステップと、各フェーズの前記輪郭に基づいて、心周期における所定期間又は所定時点における前記輪郭の長さの変化を示す変化情報を導出するステップと、前記変化情報を表示部に表示するステップと、をコンピュータに実行させるためのパラメータ処理プログラムである。
【発明の効果】
【００１６】
本発明によれば、心臓弁膜症の検査において、検査担当者及び患者の負担を低減でき、検査により得られるパラメータの精度が向上できる。
【図面の簡単な説明】
【００１７】
【図１】実施形態における画像処理装置の構成例を示すブロック図
【図２】（Ａ）〜（Ｃ）実施形態における心臓の領域を含むボリュームデータの３軸断面像を示す模式図
【図３】実施形態における左室短軸断面像における左室心筋の輪郭の一例を示す模式図
【図４】実施形態における左室垂直長軸断面像における左室心筋の輪郭の一例を示す模式図
【図５】実施形態における左室水平長軸断面像における左室心筋の輪郭の一例を示す模式図
【図６】実施形態における心周期の１周期における時間位置と輪郭長との関係性のイメージを示す模式図
【図７】実施形態における左室短軸断面像における輪郭長の変化例を示す模式図
【図８】（Ａ），（Ｂ）実施形態における左室短軸断面像における輪郭の長さの変化速度及び変化加速度の分布例を示す模式図
【図９】（Ａ），（Ｂ）実施形態における実施形態における左室垂直長軸断面像における輪郭の長さの変化速度及び変化加速度の分布例を示す模式図
【図１０】（Ａ），（Ｂ）実施形態における左室水平長軸断面像における輪郭の長さの変化速度及び変化加速度の分布例を示す模式図
【図１１】実施形態における画像処理装置の動作例を示すフローチャート
【図１２】（Ａ）〜（Ｃ）実施形態における画像処理装置がフェーズ補間しない場合を説明するための模式図
【図１３】（Ａ）〜（Ｃ）実施形態における画像処理装置がフェーズ補間する場合を説明するための模式図
【図１４】（Ａ）大動脈弁閉鎖不全と診断される心臓の左室短軸断面図、（Ｂ）大動脈弁閉鎖不全と診断される心臓の左室短軸断面図
【図１５】心臓における大動脈弁の位置を示す模式図
【発明を実施するための形態】
【００１８】
以下、本発明の実施形態について、図面を用いて説明する。
【００１９】
（本発明の一形態を得るに至った経緯）
【００２０】
従来の心エコー検査では、検査時の患者の負担は小さいが、体外から心臓部分に接触させる超音波接触子の角度によって、超音波検査装置により取得される画像が変化する。従って、心エコー検査の検査担当者によって、画像から得られる検査結果が変わる可能性がある。
【００２１】
従来の心臓カテーテル検査では、カテーテルを患者の体内に挿入し、プローブを心臓の弁付近に配置しないと、各種計測ができない。従って、検査担当者及び患者の負担が大きく、患者への侵襲の可能性もある。
【００２２】
従来の心臓ＣＴ検査では、検査担当者は、左室から駆出された血液の逆流や左室の圧力を、ＣＴ画像から確認することが困難である。また、ＣＴ画像を用いた検査結果の評価に手間を要していた。
【００２３】
特許文献１に記載された医用画像処理装置では、心臓弁膜症の診断の重要な指標となる左室圧力及び左室駆出率の情報を取得しないので、心臓弁膜症の診断精度が低下する可能性がある。
【００２４】
以下、心臓弁膜症の検査において、検査担当者及び患者の負担を低減でき、検査により得られるパラメータの精度が向上できるパラメータ処理装置、パラメータ処理方法、及びパラメータ処理プログラムについて説明する。
【００２５】
（実施形態）
図１は、実施形態における画像処理装置１００の構成例を示すブロック図である。画像処理装置１００は、画像取得部１１０、操作部１２０、表示部１３０、制御部１４０及び記憶部１５０を備える。尚、レジストレーション処理部１４２及びフェーズ補間処理部１４５は、省略されてもよい。画像処理装置１００は、ＣＴ（Ｃｏｍｐｕｔｅｄ  Ｔｏｍｏｇｒａｐｈｙ）装置２００からＣＴ画像を取得し、取得されたＣＴ画像に対して処理を行う。
【００２６】
制御部１４０は、設定部１４１、レジストレーション処理部１４２、輪郭情報処理部１４３、画像導出部１４４、及びフェーズ補間処理部１４５を備える。制御部１４０は、例えば、ＣＰＵ（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）又はＤＳＰ（Ｄｉｇｉｔａｌ  Ｓｉｇｎａｌ  Ｐｒｏｃｅｓｓｉｎｇ）を含む。制御部１４０は、例えば、ＲＯＭ（Ｒｅａｄ  Ｏｎｌｙ  Ｍｅｍｏｒｙ）又はＲＡＭ（Ｒａｎｄｏｍ  Ａｃｃｅｓｓ  Ｍｅｍｒｙ）、を含む。ＣＰＵ又はＤＳＰは、ＲＯＭ又はＲＡＭに記憶された画像処理プログラムを実行することにより、制御部１４０の各機能を実現する。
【００２７】
ＣＴ装置２００は、生体へＸ線を照射し、体内の組織によるＸ線の吸収の違いを利用して、画像（ＣＴ画像）を撮像する。ＣＴ画像は、時系列に複数撮像される。本実施形態では、時系列に並ぶ各時点をフェーズとも称する。ＣＴ装置２００は、各フェーズにおいてＣＴ画像を撮像する。ＣＴ画像は、生体内部の任意の箇所の情報を含むボリュームデータを形成する。ＣＴ画像が撮像されることにより、ＣＴ画像における各画素（ボクセル）の画素値（ＣＴ値）が得られる。
【００２８】
また、ＣＴ装置２００は、ＣＴ画像として、ボリュームレンダリングを用いてボリュームデータから三次元画像を生成してもよい。ＣＴ装置２００は、ＣＴ画像として、三次元画像を時系列に配列して三次元動画像（四次元画像）を生成してもよい。尚、三次元画像及び四次元画像の生成は、画像処理装置１００により行われてもよい。
【００２９】
画像取得部１１０は、ＣＴ装置２００からＣＴ画像を取得する。画像取得部１１０は、例えば、ＣＴ装置２００から、有線回線又は無線回線を介して通信によりＣＴ画像を取得してもよいし、任意の記憶媒体（不図示）を介して取得してもよい。
【００３０】
画像導出部１４４は、画像取得部１１０により取得されたＣＴ画像のボリュームデータから、複数の断面により切断された状態を示す断面画像（２次元画像）を導出する。例えば、画像導出部１４４は、３断面（左室短軸断面、左室垂直長軸断面、及び左室水平長軸断面）の断面画像を導出する。画像導出部１４４は、断面画像導出部の一例である。
【００３１】
図２（Ａ）〜（Ｃ）は、心臓の領域を含むボリュームデータの３軸断面像（断面画像）を示す模式図である。図２（Ａ）は、ボリュームデータＢＤが左室短軸断面Ｄ１により切断され、左室短軸断面像が形成されることを例示する。図２（Ｂ）は、ボリュームデータＢＤが左室垂直長軸断面Ｄ２により切断され、左室垂直長軸断面像が形成されることを例示する。図２（Ｃ）は、ボリュームデータＢＤが左室水平長軸断面Ｄ３により切断され、左室水平長軸断面像が形成されることを例示する。
【００３２】
操作部１２０は、例えば、タッチパネル、ポインティングデバイス（例えばマウス）、キーボードを含む。操作部１２０は、画像処理装置１００のユーザ（例えば、医師、放射線技師）から、任意の入力操作を受け付ける。
【００３３】
操作部１２０は、例えば、ボリュームデータの左室短軸断面像１０において、左室１１を包囲する心筋１２（左室心筋）の輪郭１３（図３の白点線部分）を指定する入力操作を受け付ける（図３参照）。
【００３４】
図３は、左室短軸断面像１０における左室心筋の輪郭１３の一例を示す模式図である。図３では、心筋１２の左側には右室１４が位置し、心筋１２の右上側には肺１５（図３では黒色）が位置する。左室１１には造影剤を含む血液が多く含まれており、ＣＴ値が高くなっている。
【００３５】
操作部１２０は、例えば、ボリュームデータの左室垂直長軸断面像２０において、左室１１の下壁１１ａにおける心基部２１から心尖部２２までの左室１１の外壁（心筋１２の外周端部）に沿った輪郭２３（図４の白点線矢印部分）を指定するための入力操作を受け付ける（図４参照）。図４は、左室垂直長軸断面像２０における左室心筋の輪郭２３の一例を示す模式図である。尚、図４における輪郭２３は、図３における左室心筋の輪郭に相当する。
【００３６】
操作部１２０は、例えば、ボリュームデータの左室水平長軸断面像３０において、左室１１の側壁１１ｃにおける心尖部３１から心基部３２までの左室１１の外壁（心筋１２の外周端部）に沿った輪郭３３（図５の白点線矢印部分）を指定するための入力操作を受け付ける（図５参照）。図５は、左室水平長軸断面像３０における左室心筋の輪郭３３の一例を示す模式図である。図５における輪郭３３は、図３における左室心筋の輪郭に相当する。
【００３７】
設定部１４１は、例えば操作部１２０により入力操作を受け付けた場合、入力された情報に基づいて、各種設定を行う。設定部１４１は、例えば、操作部１２０により受け付けられた輪郭１３，２３，３３を設定する。設定された輪郭１３，２３，３３は、左室短軸断面像、左室垂直長軸断面像、及び左室水平長軸断面像における指定された位置に表示される。
【００３８】
尚、設定部１４１は、操作部１２０による入力操作ではなく、ソフトウェアにより輪郭抽出されてもよい。この場合、例えば、制御部１４０が、各断面画像において左室１１、心筋１２、心基部２１、心尖部２２の位置を抽出し、各部の位置に基づいて、輪郭を導出して設定してもよい。設定部１４１は、輪郭設定部の一例である。
【００３９】
レジストレーション処理部１４２は、複数の画像（２次元画像又は３次元画像）において、レジストレーションを実行する。レジストレーションでは、関連する領域（例えば心臓の領域）の位置合わせ（位置の調整）が行われる。この場合、レジストレーション処理部１４２は、３次元画像を考慮して、立体的な位置合わせを行ってもよいし、２次元画像を考慮して、平面的な位置合わせを行ってもよい。
【００４０】
レジストレーション処理部１４２は、例えば、複数のフェーズを含むボリュームデータにおいて、各フェーズ間でレジストレーションを実行する。
【００４１】
レジストレーションの具体的な方法として、例えば公知の方法（例えば参考特許文献１参照）を用いる。
（参考特許文献１：米国特許第８３１１３００号明細書）
【００４２】
輪郭情報処理部１４３は、設定部１４１により設定された、各断面画像における輪郭１３，２３,３３の長さ（輪郭長）を導出（例えば算出）する。輪郭情報処理部１４３は、例えば、表示部１３０の画面上の輪郭１３，２３，３３の長さと、ボリュームデータに含まれるボクセルの物理的な大きさの情報とに基づいて、輪郭１３，２３，３３の長さ（実際の長さ）を導出する。また、輪郭情報処理部１４３は、心周期における輪郭１３，２３，３３の長さの変化を示す変化情報を導出する。輪郭情報処理部１４３は、変化情報に係る画面を生成する。これにより、ユーザは、輪郭情報処理部１４３により得られた輪郭長の情報を解析し、各種病状を診断できる。輪郭情報処理部１４３は、変化情報導出部の一例である。
【００４３】
例えば、輪郭情報処理部１４３は、各フェーズの断面画像に設定された輪郭１３，２３，３３に基づいて、心臓の収縮時における輪郭１３，２３，３３の長さの時間変化率（輪郭長の変化速度）を導出（例えば算出）する。
【００４４】
例えば、輪郭情報処理部１４３は、各フェーズの断面画像に設定された輪郭１３，２３，３３に基づいて、心臓の収縮期における輪郭１３，２３,３３の長さの時間変化率の時間変化率（輪郭長の変化加速度）を導出（例えば算出）する。
【００４５】
また、輪郭情報処理部１４３は、例えば、各断面画像における輪郭１３，２３，３３の長さを示すグラフを生成する。（図６，図７参照）。
【００４６】
尚、輪郭情報処理部１４３は、左室心筋の輪郭１３，２３，３３の長さに係る変化の情報の代わりに、左室心筋の輪郭１３に包囲された領域の面積に係る変化の情報を導出してもよい。つまり、輪郭情報処理部１４３は、各フェーズの断面画像に設定された輪郭１３に基づいて、心臓の収縮時における上記領域の面積の変化速度、心臓の収縮期における上記領域の面積の変化加速度を導出してもよい。
【００４７】
フェーズ補間処理部１４５は、ボリュームデータが生成される複数のフェーズの間において、フェーズを補間し、補間されたフェーズのボリュームデータを導出する（フェーズ補間）。つまり、フェーズ補間では、ボリュームデータが存在しない時点におけるボリュームデータが仮想的に生成される。フェーズ補間の具体的な方法として、例えば公知の方法（例えば参考特許文献２参照）を用いる。フェーズ補間の詳細については後述する。
（参考特許文献２：特許第５４０８４９３号公報）
【００４８】
表示部１３０は、例えばＬＣＤ（Ｌｉｑｕｉｄ  Ｃｒｙｓｔａｌ  Ｄｉｓｐｌａｙ）を含み、各種情報を表示する。表示部１３０は、各種画像（例えば、ＣＴ画像の２次元画像（断面画像）、３次元画像（静止画）、４次元画像（動画）、輪郭が付与された２次元画像）を表示する。また、表示部１３０は、例えば、変化情報を表示し、時間経過に伴う左室心筋の輪郭１３，２３，３３の変化に基づく変化情報（例えば変化加速度、変化速度）を示すグラフや分布図を表示する。また、表示部１３０は、各種操作画面を表示する。
【００４９】
記憶部１５０は、各種画像（例えば、ＣＴ画像の２次元画像（断面画像）、３次元画像（静止画）、４次元画像（動画））、各種データ、各種プログラム、各種情報（例えば手動又は自動的に設定された設定情報）、を記憶する。
【００５０】
次に、各断面画像における輪郭長の変化例について説明する。
【００５１】
図６は、心周期の１周期における時刻（時間位置）と輪郭長との関係性を示すイメージ図である。１つの心周期は、拡張期と収縮期を含む。心周期が反復されると、拡張期と収縮期とが順に反復される。拡張期には、左室１１の容量が最大となり、つまり各輪郭１３，２３，３３の長さが最長となる。収縮期には、左室１１の容量が最小となり、つまり各輪郭１３，２３，３３の長さが最短となる。尚、心周期の１周期における時間位置は、０〜１００（％）で示され、心周期における位相を示す。
【００５２】
輪郭情報処理部１４３は、心臓の収縮時の各断面における輪郭１３，２３，３３の長さの変化速度を導出する。心臓の収縮時とは、拡張期から収縮期までの間の所定期間又は所定時点である。心臓の収縮時の輪郭１３，２３，３３の長さの変化速度は、拡張期（ｅａｒｒｌｙ  ｓｙｓｔｏｌｅ）に左室１１に導入された血液が収縮期（ｅｎｄ−ｓｙｓｔｏｌｅ）に左室１１から駆出される割合に相当するので、左室駆出率に相当すると言える。
【００５３】
輪郭１３，２３，３３の長さの変化速度は、例えば、拡張期（時刻ｔ１）から収縮期（時刻ｔ２）までの所定の瞬間（時点）の瞬間変化速度の平均値又は最小値（負の値の絶対値が最大）でもよい。図６では、時刻ｔ３において、ＡＳの瞬間変化速度が最小となる。また、この変化速度は、例えば、拡張期（時刻ｔ１）における輪郭１３，２３，３３の長さと収縮期（時刻ｔ２）における輪郭１３，２３，３３の長さの傾き（点線Ｌ１の傾き）により示されてもよい。
【００５４】
輪郭情報処理部１４３は、心臓の収縮期における輪郭１３，２３，３３の長さの変化加速度を導出する。この変化加速度から、収縮期における左室１１が血液を駆出する駆出率の変化率が類推できるので、収縮期の左室の圧力に相当すると言える。
【００５５】
輪郭１３，２３，３３の長さの変化加速度は、例えば、収縮期（時刻ｔ２）における瞬間加速度でもよい。また、この変化加速度は、収縮期（時刻ｔ２）近傍に含まれる所定の瞬間（時点）の瞬間加速度の平均値又は最大値でもよい。この最大値の瞬間は、時刻ｔ２からずれていてもよい。この平均値は、収縮期近傍の所定期間の変化加速度を示すことになる。本実施形態では、主に収縮期における変化加速度を例示する。
【００５６】
尚、図６では、輪郭長の変化のグラフが１つであるが、各断面画像における輪郭１３，２３，３３の長さは３つあるので、輪郭長の変化のグラフは３つ存在する。３つのグラフにおいて、上記変化速度及び上記変化加速度が１つずつパラメータとして存在するので、全部で６つのパラメータが存在する。輪郭情報処理部１４３は、この６つのパラメータを導出する。
【００５７】
また、図６では、心周期における輪郭１３，２３，３３の長さの変化をグラフで示したが、心周期における輪郭１３に包囲された領域の面積の変化をグラフで示してもよい。
【００５８】
次に、各断面画像における実際の被験者の心筋１２の輪郭１３，２３，３３の長さの変化について説明する。ここでは、ＣＴ装置２００が実際に被験者の心臓付近のＣＴ画像を撮像し、画像処理装置１００が、ボリュームデータの各断面画像における輪郭１３，２３，３３の長さの変化を解析した。
【００５９】
図７は、左室短軸断面像１０における輪郭１３の長さの変化例を示す模式図である。図７では、横軸において、１つの心周期を０〜１００（％）の時点に分けて示し、縦軸において、心筋１２の輪郭１３の長さを示している。図７の実験データに対する考察は、図８（Ａ），（Ｂ）に示されている。
【００６０】
図８（Ａ），（Ｂ）、図９（Ａ），（Ｂ）、図１０（Ａ），（Ｂ）は、各断面画像における心臓の収縮時の輪郭長の変化速度、又は、心臓の収縮期における輪郭長の変化加速度、の分布例を示す模式図である。
【００６１】
図８（Ａ），（Ｂ）、図９（Ａ），（Ｂ）、図１０（Ａ），（Ｂ）では、「ＡＳ」は大動脈弁狭窄症の被験者に対する実験結果を示し、「ＡＲ」は大動脈弁閉鎖不全の被験者に対する実験結果を示し、「Ｃｏｎｔｒｏｌ」は健常者の被験者に対する実験結果を示す。
【００６２】
図８（Ａ）は、心臓の収縮時の左室短軸断面像１０における輪郭１３の長さの変化速度の分布例を示す模式図である。図８（Ａ）での実験サンプルの変化速度の分布は、「ＡＳ」が５．７±２．１（ｍｍ／ｍｓｅｃ）であり、「ＡＲ」が９．０±１．７（ｍｍ／ｍｓｅｃ）であり、「ｃｏｎｔｒｏｌ」が８．４±３．５（ｍｍ／ｍｓｅｃ）である。尚、「±」の記号の前の数値は、平均値を示し、「±」の記号の後の数値は、標準偏差を示す。
【００６３】
図８（Ｂ）は、心臓の収縮期の左室短軸断面像１０における輪郭１３の長さの変化加速度の分布例を示す模式図である。図９（Ｂ）での実験サンプルの変化加速度の分布は、「ＡＳ」が０．２２±０．１４（ｍｍ／ｍｓｅｃ）であり、「ＡＲ」が０．４０±０．１３（ｍｍ／ｍｓｅｃ）であり、「ｃｏｎｔｒｏｌ」が０．２９±０．１７（ｍｍ／ｍｓｅｃ）である。
【００６４】
図９（Ａ）は、心臓の収縮時の左室垂直長軸断面像２０における輪郭２３の長さの変化速度の分布例を示す模式図である。図９（Ａ）での実験サンプルの変化速度の分布は、「ＡＳ」が２．４±１．０（ｍｍ／ｍｓｅｃ）であり、「ＡＲ」が４．２±０．９（ｍｍ／ｍｓｅｃ）であり、「ｃｏｎｔｒｏｌ」が３．７±１．０（ｍｍ／ｍｓｅｃ）である。
【００６５】
図９（Ｂ）は、心臓の収縮期の左室垂直長軸断面像２０における輪郭２３の長さの変化加速度の分布例を示す模式図である。図９（Ｂ）での実験サンプルの変化加速度の分布は、「ＡＳ」が０．１１±０．０６（ｍｍ／ｍｓｅｃ）であり、「ＡＲ」が０．１６±０．０７（ｍｍ／ｍｓｅｃ）であり、「ｃｏｎｔｒｏｌ」が０．１５±０．０６（ｍｍ／ｍｓｅｃ）である。
【００６６】
図１０（Ａ）は、心臓の収縮時の左室水平長軸断面像３０における輪郭３３の長さの変化速度の分布例を示す模式図である。図１０（Ａ）での実験サンプルの変化速度の分布は、「ＡＳ」が３．０±１．１（ｍｍ／ｍｓｅｃ）であり、「ＡＲ」が３．７±０．８（ｍｍ／ｍｓｅｃ）であり、「ｃｏｎｔｒｏｌ」が３．４±０．７（ｍｍ／ｍｓｅｃ）である。
【００６７】
図１０（Ｂ）は、心臓の収縮期の左室水平長軸断面像３０における輪郭３３の長さの変化加速度の分布例を示す模式図である。図１０（Ｂ）での実験サンプルの変化加速度の分布は、「ＡＳ」が０．１３±０．０８（ｍｍ／ｍｓｅｃ）であり、「ＡＲ」が０．１６±０．０６（ｍｍ／ｍｓｅｃ）であり、「ｃｏｎｔｒｏｌ」が０．１３±０．０４（ｍｍ／ｍｓｅｃ）である。
【００６８】
このように、図８（Ａ），図９（Ａ），図１０（Ａ）では、「ＡＲ」のサンプル（□印で示す）が「ＡＳ」のサンプル（○印）よりも変化速度の値が大きい方に分布していることが理解できる。これは、ＡＳが、大動脈弁が開き難いので、左室１１から血液が駆出され難く、左室１１の収縮具合が比較的穏やかであり、左室駆出率が小さいことに起因すると考えられる。また、ＡＲが、大動脈弁が開き易いので、左室１１から血液が駆出され易く、左室１１の収縮具合が比較的急峻であり、左室駆出率が高いことに起因すると考えられる。尚、図１０（Ａ）では、ＡＳ、ＡＲ、及びＣｏｎｔｒｏｌ間で、変化速度の分布が顕著な傾向を示さなかった。
【００６９】
また、図８（Ｂ），図９（Ｂ），図１０（Ｂ）では、「ＡＲ」のサンプルが「ＡＳ」のサンプルよりも変化加速度の値が大きい方に分布している。これは、ＡＳが、大動脈弁が開きにくいので、左室１１へ血液が逆流し難く、左室１１の拡張具合が比較的緩やかであり、左室１１の圧力がある程度保たれることに起因すると考えられる。また、ＡＲが、大動脈弁が開き易いので、左室１１へ血液が逆流し易く、左室１１の拡張具合が比較的急峻であり、左室１１の圧力が下がることに起因すると考えられる。尚、図１０（Ｂ）では、ＡＳ、ＡＲ、及びＣｏｎｔｒｏｌ間で、変化速度の分布が顕著な傾向を示さなかった。
【００７０】
上記被験者に対する実験結果により、表示部１３０が、左室短軸断面像１０の輪郭１３、左室垂直長軸断面像２０の輪郭２３を表示することで、ユーザによる大動脈弁狭窄症の診断を有効に補助できることが理解できる。また、各パラメータの回帰分析等を更に検討することで、上記実験結果を大動脈弁狭窄症及び大動脈弁閉鎖不全の診断に対して一層有効に活用可能となることが期待できる。
【００７１】
尚、左室駆出率の計測など従来の手法は、心室内腔の体積などを計測しているが、本実施形態では心筋、特に心筋の外壁を計測している。これは、心室内腔には乳頭筋などの組織が存在することによって画像が不明瞭であり、特に心尖部において心室内腔の計測が安定しなかったことを鑑みてのことである。一方、心筋、特に心筋の外壁は明瞭な組織であり、安定した計測結果が期待できる。また、心筋は造影剤を用いないでも計測が出来る。
【００７２】
次に、画像処理装置１００の動作例について説明する。
図１１は、画像処理装置１００の動作例を示すフローチャートである。
【００７３】
ＣＴ装置２００は、人体の心臓付近を含むＣＴ画像を撮像する。画像処理装置１００では、ＣＴ装置２００から、複数フェーズのボリュームデータを取得する（Ｓ１０１）。
【００７４】
操作部１２０は、複数フェーズのうちの１つのフェーズにおいて、ボリュームデータから３断面（例えば、左室短軸断面、左室垂直長軸断面、左室水平長軸断面）を指定するための入力操作を受け付ける（Ｓ１０２）。設定部１４１は、操作部１２０の入力操作に応じて、上記３断面を指定する。画像導出部１４４は、３断面の断面画像（例えば、左室短軸断面像１０、左室垂直長軸断面像２０、左室水平長軸断面像３０）を導出する。
【００７５】
操作部１２０は、導出された左室短軸断面像１０において輪郭１３を指定するための入力操作を受け付ける（Ｓ１０３）。設定部１４１は、輪郭１３を設定する。表示部１３０は、輪郭１３を視認可能にするため、設定された輪郭１３を左室短軸断面像１０に重ねて表示してもよい。
【００７６】
操作部１２０は、導出された左室垂直長軸断面像２０において輪郭２３を指定するための入力操作を受け付ける（Ｓ１０４）。設定部１４１は、輪郭２３を設定する。表示部１３０は、輪郭２３を視認可能にするため、設定された輪郭２３を左室垂直長軸断面像２０に重ねて表示してもよい。
【００７７】
操作部１２０は、導出された左室水平長軸断面像３０において輪郭３３を指定するための入力操作を受け付ける（Ｓ１０５）。設定部１４１は、輪郭３３を設定する。表示部１３０は、輪郭３３を視認可能にするため、設定された輪郭３３を左室水平長軸断面像３０に重ねて表示してもよい。
【００７８】
レジストレーション処理部１４２は、Ｓ１０１において取得されたボリュームデータの複数のフェーズ間において、レジストレーションを実行する（Ｓ１０６）。つまり、レジストレーション処理部１４２は、ボリュームデータの所定の領域（例えば心臓）の位置を複数のフェーズ間で調整する。
【００７９】
設定部１４１は、レジストレーション処理部１４２によるレジストレーション結果、つまり位置合わせの結果に応じて、複数フェーズのうちの上記１つのフェーズ以外のフェーズにおいて、各断面画像における各輪郭１３，２３，３３を設定する（Ｓ１０７）。
【００８０】
輪郭情報処理部１４３は、各フェーズにおける各輪郭１３，２３，３３の長さを導出する（Ｓ１０８）。
【００８１】
輪郭情報処理部１４３は、各フェーズにおける心臓の収縮時の各輪郭１３，２３，３３の長さの変化速度を導出する（Ｓ１０９）。
【００８２】
輪郭情報処理部１４３は、各フェーズにおける心臓の収縮期における各輪郭１３，２３,３３の長さの変化加速度を導出する（Ｓ１１０）。
【００８３】
輪郭情報処理部１４３は、各フェーズにおける各輪郭１３，２３，３３の長さを示すグラフ（例えば、図６、図７参照）を生成する。表示部１３０は、生成されたグラフを表示する（Ｓ１１１）。
【００８４】
輪郭情報処理部１４３は、各フェーズにおける心臓の収縮時の各輪郭１３，２３，３３の長さの変化速度を示す分布図（例えば、図８（Ａ），図９（Ａ），図１０（Ａ））を生成する。表示部１３０は、この生成された輪郭長の変化速度の分布図を表示する（Ｓ１１２）。尚、表示部１３０は、分布図以外によって、輪郭長の変化速度を表してもよい。
【００８５】
輪郭情報処理部１４３は、各フェーズにおける心臓の収縮期における各輪郭１３，２３，３３の長さの変化加速度を示す分布図（例えば、図８（Ｂ），図９（Ｂ），図１０（Ｂ））を生成する。表示部１３０は、この生成された輪郭長の変化加速度の分布図を表示する（Ｓ１１３）。尚、表示部１３０は、分布図以外によって、輪郭長の変化加速度を表してもよい。
【００８６】
画像処理装置１００による動作例によれば、生体内部を撮像した画像から、大動脈弁狭窄症、大動脈弁閉鎖不全の診断の重要な指標となる左室圧力及び左室駆出率に相当するパラメータ（上記変化加速度、上記変化速度）を取得できる。また、画像処理装置１００は、上記変化加速度、上記変化速度の分布状態を導出して表示することで、ユーザは、分布状態の傾向から、大動脈弁狭窄症及び大動脈弁閉鎖不全を容易に区別して診断できる。
【００８７】
次に、フェーズ補間の詳細について説明する。
【００８８】
フェーズ補間処理部１４５は、フェーズ補間において、例えば、断面画像における各領域の位置、断面画像の各画素の画素値、心周期における３次元画像のピッチ（時間間隔）及び３次元画像の数、に基づいて、補間されたフェーズでのボリュームデータを導出する。
【００８９】
例えば、ＣＴ装置２００がＣＴ画像を撮像する場合、１周期の心周期では、例えば、フェーズ数は１０個〜２０個程度、つまり撮像される３次元画像は１０枚〜２０枚である。フェーズ補間処理部１４５がフェーズ補間を行うことで、フェーズ間の画像の変化を滑らかにすることができる。
【００９０】
図１２（Ａ）は、フェーズ補間しない場合の左室短軸断面像１０における心周期の位相と左室面積（左室の断面面積）との関係性を示す模式図である。図１２（Ｂ）は、フェーズ補間しない場合の各フェーズの左室短軸断面像１０を時系列に並べた例を示す模式図である。図１２（Ｃ）は、フェーズ補間しない場合の左室短軸断面像１０の一例を示す模式図である。尚、左室面積は、左室の心筋の輪郭と同様のパラメータとして扱うことができる。
【００９１】
図１３（Ａ）は、フェーズ補間する場合の左室短軸断面像１０の心周期の位相と左室面積との関係性を示す模式図である。図１３（Ｂ）は、フェーズ補間する場合の各フェーズの左室短軸断面像１０を時系列に並べた例を示す模式図である。図１３（Ｃ）は、フェーズ補間する場合における補間されたフレームでの左室短軸断面像１０の一例を示す模式図である。
【００９２】
図１２（Ａ），（Ｂ）及び図１３（Ａ），（Ｂ）を比較すると、図１３（Ａ）の方が、フレーム数が多く、各フレームにより導出される心周期の位相に対する左室面積のサンプル数が多くなり、グラフが滑らかに描かれていることが理解できる。
【００９３】
フェーズ間の画像の変化の加速度は、例えば、下記の（式１）により示される。
【００９４】
【数１】
0006418881.tif 000002
【００９５】
尚、ｎは微分を示し、Φ（ｔ）は時刻ｔにおける輪郭長を示し、±ｈは前後のフレームを示す。つまり、（式１）では、ボリュームデータから連続する３フレームの画像を抽出し、フレーム間の画像の差分をとり、２次微分を行って、時刻ｔにおける加速度を導出することを示している。
【００９６】
一般的に、心臓は時間経過に伴い、収縮又は拡張による移動以外にも、生体が動くことにより位置が多少移動する。従って、ＣＴ画像における心臓の位置が移動することになる。この場合でも、フェーズ補間処理部１４５がフェーズ補間を行うことで、例えば３断面における形状変化を追跡し易くなる。従って、輪郭長の変化を高精度に追跡でき、パラメータの導出精度を向上できる。特に、心臓はねじれながら収縮するので、フェーズ補間を一旦行うことによってフレーム間の対応関係を設定しやくなり、レジストレーションひいては輪郭長の分析精度の向上が期待できる。
【００９７】
このように、画像処理装置１００によれば、生体内部を撮像した画像から、心臓弁膜症（例えば、大動脈弁狭窄症、大動脈弁閉鎖不全）の診断の重要な指標となる左室圧力及び左室駆出率に相当する指標を取得できる。従って、超音波を用いることなく心臓弁膜症の検査を実施でき、検査担当者による検査結果のばらつきを低減できる。また、カテーテルを使用せずに心臓弁膜症の検査を実施できるので、人体への侵襲を低減でき、検査担当者及び患者の手間を低減できる。
【００９８】
このように、心臓弁膜症の検査において、検査担当者及び患者の負担を低減でき、検査により得られるパラメータの精度が向上できる。よって、ユーザは、得られたパラメータを用いて、心臓弁膜症の診断を高精度に実施できる。
【００９９】
なお、本発明は、上記実施形態の構成に限られるものではなく、特許請求の範囲で示した機能、または本実施形態の構成が持つ機能が達成できる構成であればどのようなものであっても適用可能である。
【０１００】
上記実施形態では、画像処理装置１００が、輪郭１３，２３，３３の長さの変化に基づいて、左室駆出率及び左室圧力に相当するパラメータを導出することを例示したが、これに限られない。例えば、画像処理装置１００は、断面画像における左室面積又はボリュームデータにおける左室体積に基づいて、左室駆出率及び左室の収縮期の圧力に相当するパラメータを導出してもよい。例えば、画像処理装置１００は、左室心筋の輪郭と同様に、心臓の収縮時における左室面積の変化速度を導出し、心臓の収縮期における左室面積の変化加速度を導出してもよい。
【０１０１】
上記実施形態では、画像処理装置１００は、輪郭１３，２３，３３の長さの変化速度又は変化加速度を導出する際に用いられるフィルタを備えてもよい。これにより、輪郭情報処理部１４３は、例えば差分を求めて傾きを導出する場合に、フィルタを用いることで、ノイズの影響を低減でき、変化速度又は変化加速度の導出精度を向上できる。
【０１０２】
上記実施形態では、ＣＴ装置２００により画像を撮像し、生体内部の情報を含むボリュームデータを生成することを例示したが、他の装置（例えばＭＲＩ（Ｍａｇｎｅｔｉｃ  Ｒｅｓｏｎａｎｃｅ  Ｉｍａｇｉｎｇ）装置）により画像を撮像し、ボリュームデータを生成してもよい。尚、ＭＲＩ装置を用いる場合、人口弁が体内に挿入された患者に対しては、ボリュームデータの精度が低下する可能性がある。従って、心臓弁膜症の患者に対して検査、人口弁の挿入、人口弁の挿入後の検査を継続的に実施する場合、ＣＴ装置２００を用いることで、検査精度を向上できる。
【０１０３】
上記実施形態では、ボリュームデータを切断する断面として、左室短軸断面、左室垂直長軸断面、左室水平長軸断面の３断面を用いることを例示したが、他の断面（例えば、斜めとなる断面）を用いてもよい。また、それぞれの断面方向について、複数の断面（例えば、複数の左室短軸断面）を用いても良い。また、断面は、曲面（例えば、左室短軸断面の重心を連結した曲線に沿った、いわゆるＣｕｒｖｅｄ  Ｐｌａｎａｒ  Ｒｅｃｏｎｓｔｒｕｃｔｉｏｎに用いる曲面）であっても良い。
【０１０４】
上記実施形態では、生体内部を撮像した画像から、大動脈弁狭窄症及び大動脈弁閉鎖不全を診断するためのパラメータを導出することを主に例示したが、大動脈弁狭窄症及び大動脈弁閉鎖不全以外の心臓弁膜症を診断するためのパラメータを導出してもよい。例えば、画像処理装置１００が、右室の心筋の輪郭の長さ又は輪郭に包囲された領域の面積の変化情報（例えば変化加速度、変化速度）を導出することで、右室の圧力及び駆出率（流出率）を導出できる。これにより、心臓の各弁（例えば僧帽弁、三尖弁、大動脈弁、肺動脈弁）に関する狭窄や閉鎖不全を診断するためのパラメータを、ＣＴ画像等から得ることができる。このとき、例えば、右室短軸断面のなかで左室に接していない部分の輪郭長が優良な指標となり得る。このように、肺動脈弁狭窄症及び肺動脈弁閉鎖不全を診断するためのパラメータを導出してもよい。また、心機能解析に応用することも考えられる。
【０１０５】
（本発明の一態様の概要）
本発明の一態様のパラメータ処理装置は、時系列に並ぶ複数のフェーズにおける心臓の領域を含むボリュームデータを取得する画像取得部と、各フェーズの前記ボリュームデータにおいて、所定の断面で切断された断面画像を導出する断面画像導出部と、各フェーズの前記断面画像において、心室心筋の輪郭を設定する輪郭設定部と、各フェーズの前記輪郭に基づいて、心周期における所定期間又は所定時点における前記輪郭の長さ又は前記輪郭に包囲された領域の面積の変化を示す変化情報を導出する変化情報導出部と、前記変化情報を表示する表示部と、を備える。
【０１０６】
この構成によれば、生体内部を撮像した画像から、心臓弁膜症の診断の指標に相当するパラメータを取得できる。従って、超音波を用いることなく心臓弁膜症の検査を実施でき、検査担当者による検査結果のばらつきを低減できる。また、カテーテルを使用せずに心臓弁膜症の検査を実施できるので、人体への侵襲を低減でき、検査担当者及び患者の手間を低減できる。このように、心臓弁膜症の検査において、検査担当者及び患者の負担を低減でき、検査により得られるパラメータの精度が向上できる。よって、ユーザは、得られたパラメータを用いて、心臓弁膜症の診断を高精度に実施できる。
【０１０７】
本発明の一態様のパラメータ処理装置は、前記変化情報導出部が、各フェーズの前記輪郭に基づいて、心周期における収縮期を含む収縮期近傍の所定期間又は所定時点における前記輪郭の長さ又は前記領域の面積の変化加速度を導出し、前記表示部が、前記変化加速度の情報を表示する、を備える。
【０１０８】
この構成によれば、生体内部を撮像した画像から、心臓弁膜症の診断の重要な指標となる心室（例えば左室又は右室）の圧力を類推可能なパラメータ（上記変化加速度）を取得できる。
【０１０９】
本発明の一態様のパラメータ処理装置は、前記変化情報導出部が、各フェーズの前記輪郭に基づいて、心周期における拡張期から収縮期までの間の所定期間又は所定時点における前記輪郭の長さ又は前記領域の面積の変化速度を導出し、前記表示部が、前記変化速度の情報を表示する。
【０１１０】
この構成によれば、生体内部を撮像した画像から、心臓弁膜症の診断の重要な指標となる心室（例えば左室又は右室）の駆出率を類推可能なパラメータ（上記変化速度）を取得できる。
【０１１１】
本発明の一態様のパラメータ処理装置は、前記所定の断面が、左室短軸断面又は左室長軸断面を含む。
【０１１２】
この構成によれば、一般的な左室の断面を用いて変化情報を容易に導出できる。
【０１１３】
本発明の一態様のパラメータ処理装置は、前記所定の断面が複数存在し、前記輪郭設定部が、各フェーズの複数の断面画像において、前記輪郭を設定し、前記変化情報導出部が、各フェーズの複数の前記輪郭に基づいて、複数の前記変化情報を導出し、前記表示部が、複数の前記変化情報を表示する。
【０１１４】
この構成によれば、複数の断面画像において心室心筋に係る複数の変化情報を導出し、表示するので、ユーザは、複数の変化情報（例えば異なる断面画像を用いた変化情報）を参照して、心臓弁膜症の診断を容易に実施できる。
【０１１５】
本発明の一態様のパラメータ処理装置は、各フェーズの断面画像における前記輪郭の位置を調整するレジストレーション処理部を備える。
【０１１６】
この構成によれば、各フェーズの断面画像において心臓の位置が移動する場合でも、心臓の動きの変化に容易に追従できる。
【０１１７】
本発明の一態様のパラメータ処理装置は、連続するフェーズの間においてフェーズを補間して前記ボリュームデータを導出するフェーズ補間処理部を備え、前記断面画像導出部が、前記画像取得部により取得されたフェーズ及び前記フェーズ補間処理部により補間されたフェーズのボリュームデータにおいて、前記断面画像を導出する。
【０１１８】
この構成によれば、例えば１つの心周期において導出される上記輪郭又は上記領域の面積の数を増大できるので、各フェーズの上記輪郭の長さ又は上記領域の面積に基づいて導出される変化情報の精度を向上できる。従って、心臓の位置が移動する場合でも、心臓の動きの変化に容易に追従できる。
【０１１９】
本発明の一態様のパラメータ処理装置は、前記表示部が、前記輪郭設定部により設定された前記輪郭の長さ又は前記領域の面積の変化を示すグラフを表示する。
【０１２０】
この構成によれば、ユーザは、上記輪郭又は上記領域の面積の変化を容易に視認でき、上記輪郭又は上記領域の面積の変化を直感的に把握できる。
【０１２１】
本発明の一態様のパラメータ処理方法は、時系列に並ぶ複数のフェーズにおける心臓の領域を含むボリュームデータを取得するステップと、各フェーズの前記ボリュームデータにおいて、所定の断面で切断された断面画像を導出するステップと、各フェーズの前記断面画像において、心室心筋の輪郭を設定するステップと、各フェーズの前記輪郭に基づいて、心周期における所定期間又は所定時点における前記輪郭の長さ又は前記輪郭に包囲された領域の面積の変化を示す変化情報を導出するステップと、前記変化情報を表示部に表示するステップと、を備える。
【０１２２】
この方法によれば、生体内部を撮像した画像から、心臓弁膜症の診断の指標に相当するパラメータを取得できる。従って、超音波を用いることなく心臓弁膜症の検査を実施でき、検査担当者による検査結果のばらつきを低減できる。また、カテーテルを使用せずに心臓弁膜症の検査を実施できるので、人体への侵襲を低減でき、検査担当者及び患者の手間を低減できる。このように、心臓弁膜症の検査において、検査担当者及び患者の負担を低減でき、検査により得られるパラメータの精度が向上できる。よって、ユーザは、得られたパラメータを用いて、心臓弁膜症の診断を高精度に実施できる。
【０１２３】
本発明の一態様のパラメータ処理プログラムは、上記パラメータ処理方法の各ステップをコンピュータに実行させるためのプログラムである。
【０１２４】
このプログラムによれば、生体内部を撮像した画像から、心臓弁膜症の診断の指標に相当するパラメータを取得できる。従って、超音波を用いることなく心臓弁膜症の検査を実施でき、検査担当者による検査結果のばらつきを低減できる。また、カテーテルを使用せずに心臓弁膜症の検査を実施できるので、人体への侵襲を低減でき、検査担当者及び患者の手間を低減できる。このように、心臓弁膜症の検査において、検査担当者及び患者の負担を低減でき、検査により得られるパラメータの精度が向上できる。よって、ユーザは、得られたパラメータを用いて、心臓弁膜症の診断を高精度に実施できる。
【産業上の利用可能性】
【０１２５】
本発明は、心臓弁膜症の検査において、検査担当者及び患者の負担を低減でき、検査により得られるパラメータの精度が向上できるパラメータ処理装置、パラメータ処理方法、及びパラメータ処理プログラム等に有用である。
【符号の説明】
【０１２６】
１０  左室短軸断面像
１１  左室
１２  心筋
１３，２３，３３  輪郭
１４  右室
１５  肺
２０  左室垂直長軸断面像
２１，３２  心基部
２２，３１  心尖部
３０  左室水平長軸断面像
１００  画像処理装置
１１０  画像取得部
１２０  操作部
１３０  表示部
１４０  制御部
１４１  設定部
１４２  レジストレーション処理部
１４３  輪郭情報処理部
１４４  画像導出部
１４５  フェーズ補間処理部
１５０  記憶部
２００  ＣＴ装置
【図１】
0006418881.tif 000003
【図２】
0006418881.tif 000004
【図３】
0006418881.tif 000005
【図４】
0006418881.tif 000006
【図５】
0006418881.tif 000007
【図６】
0006418881.tif 000008
【図７】
0006418881.tif 000009
【図８】
0006418881.tif 000010
【図９】
0006418881.tif 000011
【図１０】
0006418881.tif 000012
【図１１】
0006418881.tif 000013
【図１２】
0006418881.tif 000014
【図１３】
0006418881.tif 000015
【図１４】
0006418881.tif 000016
【図１５】
0006418881.tif 000017