JP3715580B2 - 医用運動解析装置及びその方法 - Google Patents
医用運動解析装置及びその方法 Download PDFInfo
- Publication number
- JP3715580B2 JP3715580B2 JP2002077315A JP2002077315A JP3715580B2 JP 3715580 B2 JP3715580 B2 JP 3715580B2 JP 2002077315 A JP2002077315 A JP 2002077315A JP 2002077315 A JP2002077315 A JP 2002077315A JP 3715580 B2 JP3715580 B2 JP 3715580B2
- Authority
- JP
- Japan
- Prior art keywords
- motion
- analysis
- medical
- component
- heart
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Description
【０００１】
【発明の属する技術分野】
本発明は、臓器などの運動を解析する医用運動解析装置に関する。
【０００２】
【従来の技術】
心臓病等の診断において、臓器運動を正確に測定することは重要である。このため、これまでに臓器部分の２次元または３次元運動情報を取得する手法が提案されている（西浦正英 他、“部分形状拘束輪郭モデルによる超音波心壁動的輪郭抽出法”，電子情報通信学会論文誌（Ｄ−II）, ｖｏｌ．Ｊ８３−Ｄ−II，ｎｏ．1 ，ｐｐ．１８３−１９０，Ｊａｎ．２０００）、（梅原幹雄 他、電子情報通信学会論文誌（Ｄ−II），ｖｏｌ．Ｊ７９−Ｄ−II，ｎｏ．２，ｐｐ．２８６−２９４，Ｆｅｂ．１９９６）。
【０００３】
これらの手法では直接得られる運動情報とは、直交するｘ軸、ｙ軸、ｚ軸方向の位置変化や速度である。
【０００４】
一方、心臓の運動を例にとってみると、心室壁の典型的な運動は、収縮／拡張方向の運動と、ねじれや回転等による収縮／拡張方向とは異なる方向の運動に分けて考えることができる。心臓病の診断等においては、運動の主要成分である収縮／拡張方向を主軸にとった運動情報の方が意味理解に適している。
【０００５】
こうした観点から、従来は事前に仮定された収縮中心へ向かう方向の運動成分を得る方法や、別途得られる輪郭線に直交した方向の運動成分を得る方法がある（特開平１０−９９３３４等）。
【０００６】
【発明が解決しようとする課題】
しかし、これらの方法はあくまでも運動方向は仮定されたものであり、実際の運動方向を反映するものではない。
【０００７】
そこで、本発明は、臓器の運動方向についての仮定を行うこと無しに、診断上有用な運動情報を得ることができる医用運動解析装置及びその方法を提供する。
【０００８】
【課題を解決するための手段】
請求項１の発明は、心臓を含む臓器を解析対象物として、その解析対象物の運動を解析する医用運動解析装置において、前記解析対象物の各部分の時系列運動情報を取得する運動情報取得手段と、 前記時系列運動情報を主成分分析して主成分の方向を求め、この主成分の方向が、前記解析対象物の拡張収縮方向を示す解析座標系の主軸であると意味付けする運動方向推定手段と、前記推定された解析座標系における前記解析対象物の各部分の運動成分を算出する運動成分算出手段と、を具備することを特徴とする医用運動解析装置である。
【０００９】
請求項２の発明は、前記時系列運動情報が、前記解析対象物の各部分の時系列速度データであることを特徴とする請求項１記載の医用運動解析装置である。
【００１０】
請求項３の発明は、前記時系列運動情報が、前記解析対象物の各部分の時系列位置データであることを特徴とする請求項１記載の医用運動解析装置である。
【００１１】
請求項４の発明は、前記時系列運動情報が、前記解析対象物の各部分の時系列組織歪データを含み、前記解析対象物の各部分において前記解析座標系での歪データを算出する歪データ変換手段を具備することを特徴とする請求項１から３のうちいずれか一項に記載の医用運動解析装置である。
【００１２】
請求項５の発明は、前記時系列運動情報が、前記解析対象物の各部分の時系列組織歪速度データを含むことを特徴とする請求項１から３のうちいずれか一項に記載の医用運動解析装置である。
【００１３】
請求項６の発明は、前記解析対象物が心臓であり、前記解析座標系における前記心臓の複数の部分の運動成分から前記心臓の収縮中心を推定する収縮中心推定手段を具備することを特徴とする請求項１から５のうちいずれか一項に記載の医用運動解析装置である。
【００１４】
請求項７の発明は、前記解析対象物が心臓であり、前記時系列運動情報が、前記心臓の心拍から算出された期間内の時系列運動情報であることを特徴とする請求項１から６のうちいずれか一項に記載の医用運動解析装置である。
【００１５】
請求項８の発明は、心臓を含む臓器を解析対象物として、その解析対象物の運動をコンピュータによって解析する医用運動解析方法において、前記コンピュータによって取得された前記解析対象物の各部分の時系列運動情報を主成分分析して主成分の方向を求め、この主成分の方向が、前記解析対象物の拡張収縮方向を示す解析座標系の主軸であると前記コンピュータによって意味付けする運動方向推定ステップと、前前記推定された解析座標系における前記解析対象物の各部分の運動成分を前記コンピュータによって算出する運動成分算出ステップと、を具備することを特徴とする医用運動解析方法である。
【００１６】
請求項９の発明は、心臓を含む臓器を解析対象物として、その解析対象物の運動を解析するための医用運動解析プログラムであって、コンピュータに、前記解析対象物の各部分の時系列運動情報を取得する運動情報取得機能と、前記時系列運動情報を主成分分析して主成分の方向を求め、この主成分の方向が、前記解析対象物の拡張収縮方向を示す解析座標系の主軸であると意味付けする運動方向推定機能と、前記推定された解析座標系における前記解析対象物の各部分の運動成分を算出する運動成分算出機能と、を実現させる医用運動解析プログラムである。
【００１７】
本発明によれは、心臓や胃等の臓器の運動方向についての仮定や、仮定に基づく手動設定作業無しに、診断上有用な運動情報を提供することが可能となる。
【００１８】
【発明の実施の形態】
（第１の実施例）
以下、本発明の第１の実施例を図１から図６に基づいて説明する。
【００１９】
図１は、本実施例の医用運動解析装置１０の構成図である。この医用運動解析装置１０では、臓器の一つである心臓を解析対象物とする。
【００２０】
図１に示すように、医用運動解析装置１０は、解析対象物の運動情報を得る運動情報取得部１２と、得られた運動情報の時系列データを記憶する運動情報記憶部１４と、運動情報の時系列データから運動方向の主要成分方向を推定する運動方向推定部１６と、推定された運動方向を主軸とした運動成分値を算出する運動成分算出部１８と、算出された運動成分値を出力または表示する出力部２０とで構成される。
【００２１】
これら各部１２から２０の機能は、コンピュータに記憶されたプログラムによって実現できる。
【００２２】
以下に医用運動解析装置１０の処理の流れに沿って説明する。
【００２３】
まず、超音波診断装置、ＭＲＩ、ＣＴ等の別の手段によって得られた心臓の動きを撮影した時系列画像を取得する。
【００２４】
そして、この時系列画像を解析して、心臓の各部分の運動情報が運動情報取得部１２に入力される。ここで心臓の各部分の運動情報とは、各部分における位置座標や速度データ等の時系列データである。
【００２５】
時系列の画像から、解析対象物の各点を追跡する手法としては、例えば輝度勾配からオプティカルフローを算出する手法や、パターンマッチングによる追跡法がある。
【００２６】
パターンマッチングでは、ある時刻における画像において、追跡点付近の画像をテンプレート画像として取得し、他の時刻における画像について、テンプレート画像に最も似たパターンを持つ位置を追跡点として検出する。似たパターンの評価は、例えば相関や差分等を用いる方法がある。
【００２７】
このような手法を用いて、解析対象物の各部分の位置座標を得ることができ、位置座標を時間微分することで速度を得ることができる。
【００２８】
以下、本実施例では、各部分の運動情報が２次元の速度データである場合について説明する。速度データは、心臓の各部位各点におけるｘ方向成分値、ｙ方向成分値である。
【００２９】
入力された速度データｖ＝（ｖｘ ，ｖｙ ）は、運動情報記憶部１４に順次記憶される。ここで、ｖｘ ，ｖｙ はそれぞれ、ｘ軸方向、ｙ軸方向の速度成分である。
【００３０】
次に、運動方向推定部１６では、速度データの時系列データから主要な運動方向の推定を行う。運動方向の推定は、例えば主成分分析を用いる。今、時刻ｔにおける速度をｖｔ ＝（ｖｘ ｔ ，ｖｙ ｔ ) と表す。
【００３１】
【数１】
となる。速度データの共分散行列Ｃは、
【数２】
と計算される。主成分分析は、共分散行列Ｃの固有ベクトルを求めることで計算される。つまり、
【数３】
となる固有ベクトルｐｋ を求める。主成分はｐ１ である。
【００３２】
運動方向推定を概念的に図示したものが図２である。
【００３３】
図２において、各時刻における速度がプロットで表されている。主成分分析により、速度のばらつきの最も大きい方向が主成分ｐ１ として算出される。つまり、心臓の収縮／拡張運動のように主要な運動成分が収縮／拡張方向である場合には、主成分ｐ１ としてこの収縮／拡張方向が算出されることになる。
【００３４】
そして、この主成分ｐ１ の方向が、心臓の収縮／拡張方向を示す解析座標系の主軸となる。
【００３５】
次に、主成分分析で求まった解析座標に基づいて、心臓の各部位各点におけるｘ方向成分値、ｙ方向成分値の速度データを変換する。すなわち、心臓の各部位各点における速度ベクトル方向の運動成分Ｖｋ が、下記の式を用いて運動成分算出部１８で計算される。
【００３６】
【数４】
これにより、元はｘ軸方向成分、ｙ軸運動成分で表されていた速度を、診断上の意味を理解しやすい心臓の収縮／拡張方向成分とそれに直交する方向成分の解析座標系で表すことができる。
【００３７】
また、心臓の収縮／拡張方向に注目することで、診断上意味のある正負の符号の付いた速度の大きさを算出することが可能となる。従来のｘ軸方向成分、ｙ軸推定成分で表された速度では、診断上意味のある動き方向が不明なため、意味ある符号付けはできず、速度の大きさは絶対値で表すか（図３参照）、または、仮定された収縮中心へ向かう方向の成分を算出するしかなかった。しかし本実施例によると、診断上意味のある動き方向が速度データ自身から算出されるため、この軸方向に関して正負の符号が付いた速度の大きさを得ることが可能となる（図４参照）。つまり、収縮方向の速度であるのか、拡張方向の速度であるのかを区別することが可能となる。
【００３８】
また、主軸方向以外の方向の速度成分に対する意味付けも可能となる。例えば、主軸方向を収縮／拡張方向とした場合、他方の軸方向は回転方向と意味付けることも可能である。
【００３９】
最後に出力部２０では、主成分方向及びその他方向の方向ベクトルと、主成分方向及びその他方向の速度成分の時系列データが出力される。図示されないディスプレイーやプリンター等の別手段において、速度データのグラフ表示や画像への重畳表示等を行うことが可能である。また、画像上に図５に示すように主成分方向を明示すると、運動状態及び出力される速度データの理解の助けとなる。
【００４０】
なお、図６に示すように、付加的な機能として、心臓壁の複数部分の速度データについての複数の主成分方向から、最小二乗法等により収縮中心の推定を行うことも可能である。
【００４１】
（第２の実施例）
次に本発明の第２の実施例について図７から図１２に基づいて説明する。
【００４２】
本実施例の医用運動解析装置１００は、運動情報として速度データと心筋の歪データが取得される場合であり、主要な運動方向を主軸とする座標系において、歪（伸び縮みの度合い）を求める例である。また本実施例では、運動方向算出に用いる時系列速度データの時間区間を、心拍情報から算出する。
【００４３】
本実施例の医用運動解析装置１００の構成図を図７に示す。
【００４４】
図７に示すように、医用運動解析装置１００は、心臓の各部分における速度データと歪データを得る運動情報取得部１１２と、得られた時系列速度データを記憶する速度データ記憶部１１４と、歪データを記憶する歪データ記憶部１１６と、時系列速度データから運動方向の主要成分方向を推定する運動方向推定部１１８と、推定された運動方向を主軸とした速度成分値を算出する速度成分算出部１２０と、運動方向推定部１１８で推定された主要運動方向を主軸とする座標系に歪データを変換する歪データ変換部１２２と、変換後の速度成分値及び歪データを出力または表示する出力部１２４とで構成される。
【００４５】
運動情報取得部１１２では、速度データと歪データが順次入力される。速度データは、心臓の各部位各点におけるｘ方向成分値、ｙ方向成分値である。歪データは、ｘ方向の伸縮を表すｅｘｘ、ｙ方向の伸縮を表すｅｙｙ、ずれを表すｅｘｙといった、歪テンソルの値が心臓各部位各点について入力される。これらの歪データにより、図９に示すような伸縮とずれを表すことができる。速度データは速度データ記憶部１１４に記憶され、歪データは歪データ記憶部１１６に記憶される。歪テンソルは、次式である。
【００４６】
【数５】
次に、第１の実施例と同様に、時系列の速度データから心臓各部位における運動の主軸が主成分分析により算出される。第２の実施例では、例えば図８に示すように、主成分分析に用いる時系列データの時間幅を心拍情報から算出する。
【００４７】
具体例で説明すると、心拍数が６０であった場合には、１心拍にに要する時間は１秒であるので、過去１秒分の時系列データからの運動の主軸を算出すると良い。あるいは、安定性向上のために過去数心拍分の時系列データを用いても良いし、逆に運動変化への追従性向上のために、半心拍分の時系列データを用いても良い。いずれにしても、心拍情報を基に、使用する時間区間を決定することで、心臓の運動に対する合理的な時間区間を適応的に算出することが可能となる。
【００４８】
心臓の運動では、心筋の運動方向と心筋厚の変化方向は近いと考えることができる。そこで、運動方向推定部１１８で推定された主要運動方向を主軸に持つ座標系に、歪データを変換する。推定された主軸と、元のｘ軸との角度をφとすると、変換後の歪テンソルは次式で計算される。
【００４９】
【数６】
ここで、Ｒは角度をφの回転変換である。このような変換により、運動方向を解析座標系のＸ軸とする歪データを表すことで、心臓の動きをより適切に説明することのできる量、例えば心筋厚の変化量を得ることが可能となる（図１０、１１）。そして、変換後の歪テンソルＥｐ の要素であるｅｐｘｘ は心筋厚の歪を表し、ｅｐｙｙ は心筋に沿った方向の歪を表すことになる。
【００５０】
最後に出力部１２２では、変換された速度データと歪データが出力される。また、図１２のように心筋の部分ごとに心筋厚の歪分布をグラフ化しても良い。このようにすることで、心筋のどの部分が異常であるかを把握しやすい。
【００５１】
（変更例）
２つの実施例においては、２次元の速度データを扱う場合について例示したが、３次元の速度データに対しても同様の処理をすることができる。また速度データや歪データだけでなく、歪速度等、他の量に対しても同様の処理をすることができる。
【００５２】
【発明の効果】
本発明によれば、臓器の運動方向についての仮定や仮定に基づく手動設定作業無しに、診断上有用な運動情報を提供することが可能となり、診断の高精度化及び効率化に大きく寄与する。
【図面の簡単な説明】
【図１】第１の実施例のの医用運動解析装置の構成図である。
【図２】時系列速度ベクトルの主成分方向の例である。
【図３】絶対値で表示された速度の例である。
【図４】主成分方向成分の正負区別された速度の例である。
【図５】主成分方向の表示例である。
【図６】収縮中心推定の説明図である。
【図７】第２の実施例のの医用運動解析装置の構成図である。
【図８】運動方向推定に用いる区間を説明する図である。
【図９】歪の様子を説明する図である。
【図１０】変換前の座標系における歪の説明図である。
【図１１】変換後の座標系における歪の説明図である。
【図１２】心筋厚歪分布の表示例である。
【符号の説明】
１０ 医用運動解析装置
１２ 運動情報取得部
１４ 運動情報記憶部
１６ 運動方向推定部
１８ 運動成分算出部
２０ 出力部
【発明の属する技術分野】
本発明は、臓器などの運動を解析する医用運動解析装置に関する。
【０００２】
【従来の技術】
心臓病等の診断において、臓器運動を正確に測定することは重要である。このため、これまでに臓器部分の２次元または３次元運動情報を取得する手法が提案されている（西浦正英 他、“部分形状拘束輪郭モデルによる超音波心壁動的輪郭抽出法”，電子情報通信学会論文誌（Ｄ−II）, ｖｏｌ．Ｊ８３−Ｄ−II，ｎｏ．1 ，ｐｐ．１８３−１９０，Ｊａｎ．２０００）、（梅原幹雄 他、電子情報通信学会論文誌（Ｄ−II），ｖｏｌ．Ｊ７９−Ｄ−II，ｎｏ．２，ｐｐ．２８６−２９４，Ｆｅｂ．１９９６）。
【０００３】
これらの手法では直接得られる運動情報とは、直交するｘ軸、ｙ軸、ｚ軸方向の位置変化や速度である。
【０００４】
一方、心臓の運動を例にとってみると、心室壁の典型的な運動は、収縮／拡張方向の運動と、ねじれや回転等による収縮／拡張方向とは異なる方向の運動に分けて考えることができる。心臓病の診断等においては、運動の主要成分である収縮／拡張方向を主軸にとった運動情報の方が意味理解に適している。
【０００５】
こうした観点から、従来は事前に仮定された収縮中心へ向かう方向の運動成分を得る方法や、別途得られる輪郭線に直交した方向の運動成分を得る方法がある（特開平１０−９９３３４等）。
【０００６】
【発明が解決しようとする課題】
しかし、これらの方法はあくまでも運動方向は仮定されたものであり、実際の運動方向を反映するものではない。
【０００７】
そこで、本発明は、臓器の運動方向についての仮定を行うこと無しに、診断上有用な運動情報を得ることができる医用運動解析装置及びその方法を提供する。
【０００８】
【課題を解決するための手段】
請求項１の発明は、心臓を含む臓器を解析対象物として、その解析対象物の運動を解析する医用運動解析装置において、前記解析対象物の各部分の時系列運動情報を取得する運動情報取得手段と、 前記時系列運動情報を主成分分析して主成分の方向を求め、この主成分の方向が、前記解析対象物の拡張収縮方向を示す解析座標系の主軸であると意味付けする運動方向推定手段と、前記推定された解析座標系における前記解析対象物の各部分の運動成分を算出する運動成分算出手段と、を具備することを特徴とする医用運動解析装置である。
【０００９】
請求項２の発明は、前記時系列運動情報が、前記解析対象物の各部分の時系列速度データであることを特徴とする請求項１記載の医用運動解析装置である。
【００１０】
請求項３の発明は、前記時系列運動情報が、前記解析対象物の各部分の時系列位置データであることを特徴とする請求項１記載の医用運動解析装置である。
【００１１】
請求項４の発明は、前記時系列運動情報が、前記解析対象物の各部分の時系列組織歪データを含み、前記解析対象物の各部分において前記解析座標系での歪データを算出する歪データ変換手段を具備することを特徴とする請求項１から３のうちいずれか一項に記載の医用運動解析装置である。
【００１２】
請求項５の発明は、前記時系列運動情報が、前記解析対象物の各部分の時系列組織歪速度データを含むことを特徴とする請求項１から３のうちいずれか一項に記載の医用運動解析装置である。
【００１３】
請求項６の発明は、前記解析対象物が心臓であり、前記解析座標系における前記心臓の複数の部分の運動成分から前記心臓の収縮中心を推定する収縮中心推定手段を具備することを特徴とする請求項１から５のうちいずれか一項に記載の医用運動解析装置である。
【００１４】
請求項７の発明は、前記解析対象物が心臓であり、前記時系列運動情報が、前記心臓の心拍から算出された期間内の時系列運動情報であることを特徴とする請求項１から６のうちいずれか一項に記載の医用運動解析装置である。
【００１５】
請求項８の発明は、心臓を含む臓器を解析対象物として、その解析対象物の運動をコンピュータによって解析する医用運動解析方法において、前記コンピュータによって取得された前記解析対象物の各部分の時系列運動情報を主成分分析して主成分の方向を求め、この主成分の方向が、前記解析対象物の拡張収縮方向を示す解析座標系の主軸であると前記コンピュータによって意味付けする運動方向推定ステップと、前前記推定された解析座標系における前記解析対象物の各部分の運動成分を前記コンピュータによって算出する運動成分算出ステップと、を具備することを特徴とする医用運動解析方法である。
【００１６】
請求項９の発明は、心臓を含む臓器を解析対象物として、その解析対象物の運動を解析するための医用運動解析プログラムであって、コンピュータに、前記解析対象物の各部分の時系列運動情報を取得する運動情報取得機能と、前記時系列運動情報を主成分分析して主成分の方向を求め、この主成分の方向が、前記解析対象物の拡張収縮方向を示す解析座標系の主軸であると意味付けする運動方向推定機能と、前記推定された解析座標系における前記解析対象物の各部分の運動成分を算出する運動成分算出機能と、を実現させる医用運動解析プログラムである。
【００１７】
本発明によれは、心臓や胃等の臓器の運動方向についての仮定や、仮定に基づく手動設定作業無しに、診断上有用な運動情報を提供することが可能となる。
【００１８】
【発明の実施の形態】
（第１の実施例）
以下、本発明の第１の実施例を図１から図６に基づいて説明する。
【００１９】
図１は、本実施例の医用運動解析装置１０の構成図である。この医用運動解析装置１０では、臓器の一つである心臓を解析対象物とする。
【００２０】
図１に示すように、医用運動解析装置１０は、解析対象物の運動情報を得る運動情報取得部１２と、得られた運動情報の時系列データを記憶する運動情報記憶部１４と、運動情報の時系列データから運動方向の主要成分方向を推定する運動方向推定部１６と、推定された運動方向を主軸とした運動成分値を算出する運動成分算出部１８と、算出された運動成分値を出力または表示する出力部２０とで構成される。
【００２１】
これら各部１２から２０の機能は、コンピュータに記憶されたプログラムによって実現できる。
【００２２】
以下に医用運動解析装置１０の処理の流れに沿って説明する。
【００２３】
まず、超音波診断装置、ＭＲＩ、ＣＴ等の別の手段によって得られた心臓の動きを撮影した時系列画像を取得する。
【００２４】
そして、この時系列画像を解析して、心臓の各部分の運動情報が運動情報取得部１２に入力される。ここで心臓の各部分の運動情報とは、各部分における位置座標や速度データ等の時系列データである。
【００２５】
時系列の画像から、解析対象物の各点を追跡する手法としては、例えば輝度勾配からオプティカルフローを算出する手法や、パターンマッチングによる追跡法がある。
【００２６】
パターンマッチングでは、ある時刻における画像において、追跡点付近の画像をテンプレート画像として取得し、他の時刻における画像について、テンプレート画像に最も似たパターンを持つ位置を追跡点として検出する。似たパターンの評価は、例えば相関や差分等を用いる方法がある。
【００２７】
このような手法を用いて、解析対象物の各部分の位置座標を得ることができ、位置座標を時間微分することで速度を得ることができる。
【００２８】
以下、本実施例では、各部分の運動情報が２次元の速度データである場合について説明する。速度データは、心臓の各部位各点におけるｘ方向成分値、ｙ方向成分値である。
【００２９】
入力された速度データｖ＝（ｖｘ ，ｖｙ ）は、運動情報記憶部１４に順次記憶される。ここで、ｖｘ ，ｖｙ はそれぞれ、ｘ軸方向、ｙ軸方向の速度成分である。
【００３０】
次に、運動方向推定部１６では、速度データの時系列データから主要な運動方向の推定を行う。運動方向の推定は、例えば主成分分析を用いる。今、時刻ｔにおける速度をｖｔ ＝（ｖｘ ｔ ，ｖｙ ｔ ) と表す。
【００３１】
【数１】
【数２】
【数３】
【００３２】
運動方向推定を概念的に図示したものが図２である。
【００３３】
図２において、各時刻における速度がプロットで表されている。主成分分析により、速度のばらつきの最も大きい方向が主成分ｐ１ として算出される。つまり、心臓の収縮／拡張運動のように主要な運動成分が収縮／拡張方向である場合には、主成分ｐ１ としてこの収縮／拡張方向が算出されることになる。
【００３４】
そして、この主成分ｐ１ の方向が、心臓の収縮／拡張方向を示す解析座標系の主軸となる。
【００３５】
次に、主成分分析で求まった解析座標に基づいて、心臓の各部位各点におけるｘ方向成分値、ｙ方向成分値の速度データを変換する。すなわち、心臓の各部位各点における速度ベクトル方向の運動成分Ｖｋ が、下記の式を用いて運動成分算出部１８で計算される。
【００３６】
【数４】
【００３７】
また、心臓の収縮／拡張方向に注目することで、診断上意味のある正負の符号の付いた速度の大きさを算出することが可能となる。従来のｘ軸方向成分、ｙ軸推定成分で表された速度では、診断上意味のある動き方向が不明なため、意味ある符号付けはできず、速度の大きさは絶対値で表すか（図３参照）、または、仮定された収縮中心へ向かう方向の成分を算出するしかなかった。しかし本実施例によると、診断上意味のある動き方向が速度データ自身から算出されるため、この軸方向に関して正負の符号が付いた速度の大きさを得ることが可能となる（図４参照）。つまり、収縮方向の速度であるのか、拡張方向の速度であるのかを区別することが可能となる。
【００３８】
また、主軸方向以外の方向の速度成分に対する意味付けも可能となる。例えば、主軸方向を収縮／拡張方向とした場合、他方の軸方向は回転方向と意味付けることも可能である。
【００３９】
最後に出力部２０では、主成分方向及びその他方向の方向ベクトルと、主成分方向及びその他方向の速度成分の時系列データが出力される。図示されないディスプレイーやプリンター等の別手段において、速度データのグラフ表示や画像への重畳表示等を行うことが可能である。また、画像上に図５に示すように主成分方向を明示すると、運動状態及び出力される速度データの理解の助けとなる。
【００４０】
なお、図６に示すように、付加的な機能として、心臓壁の複数部分の速度データについての複数の主成分方向から、最小二乗法等により収縮中心の推定を行うことも可能である。
【００４１】
（第２の実施例）
次に本発明の第２の実施例について図７から図１２に基づいて説明する。
【００４２】
本実施例の医用運動解析装置１００は、運動情報として速度データと心筋の歪データが取得される場合であり、主要な運動方向を主軸とする座標系において、歪（伸び縮みの度合い）を求める例である。また本実施例では、運動方向算出に用いる時系列速度データの時間区間を、心拍情報から算出する。
【００４３】
本実施例の医用運動解析装置１００の構成図を図７に示す。
【００４４】
図７に示すように、医用運動解析装置１００は、心臓の各部分における速度データと歪データを得る運動情報取得部１１２と、得られた時系列速度データを記憶する速度データ記憶部１１４と、歪データを記憶する歪データ記憶部１１６と、時系列速度データから運動方向の主要成分方向を推定する運動方向推定部１１８と、推定された運動方向を主軸とした速度成分値を算出する速度成分算出部１２０と、運動方向推定部１１８で推定された主要運動方向を主軸とする座標系に歪データを変換する歪データ変換部１２２と、変換後の速度成分値及び歪データを出力または表示する出力部１２４とで構成される。
【００４５】
運動情報取得部１１２では、速度データと歪データが順次入力される。速度データは、心臓の各部位各点におけるｘ方向成分値、ｙ方向成分値である。歪データは、ｘ方向の伸縮を表すｅｘｘ、ｙ方向の伸縮を表すｅｙｙ、ずれを表すｅｘｙといった、歪テンソルの値が心臓各部位各点について入力される。これらの歪データにより、図９に示すような伸縮とずれを表すことができる。速度データは速度データ記憶部１１４に記憶され、歪データは歪データ記憶部１１６に記憶される。歪テンソルは、次式である。
【００４６】
【数５】
【００４７】
具体例で説明すると、心拍数が６０であった場合には、１心拍にに要する時間は１秒であるので、過去１秒分の時系列データからの運動の主軸を算出すると良い。あるいは、安定性向上のために過去数心拍分の時系列データを用いても良いし、逆に運動変化への追従性向上のために、半心拍分の時系列データを用いても良い。いずれにしても、心拍情報を基に、使用する時間区間を決定することで、心臓の運動に対する合理的な時間区間を適応的に算出することが可能となる。
【００４８】
心臓の運動では、心筋の運動方向と心筋厚の変化方向は近いと考えることができる。そこで、運動方向推定部１１８で推定された主要運動方向を主軸に持つ座標系に、歪データを変換する。推定された主軸と、元のｘ軸との角度をφとすると、変換後の歪テンソルは次式で計算される。
【００４９】
【数６】
【００５０】
最後に出力部１２２では、変換された速度データと歪データが出力される。また、図１２のように心筋の部分ごとに心筋厚の歪分布をグラフ化しても良い。このようにすることで、心筋のどの部分が異常であるかを把握しやすい。
【００５１】
（変更例）
２つの実施例においては、２次元の速度データを扱う場合について例示したが、３次元の速度データに対しても同様の処理をすることができる。また速度データや歪データだけでなく、歪速度等、他の量に対しても同様の処理をすることができる。
【００５２】
【発明の効果】
本発明によれば、臓器の運動方向についての仮定や仮定に基づく手動設定作業無しに、診断上有用な運動情報を提供することが可能となり、診断の高精度化及び効率化に大きく寄与する。
【図面の簡単な説明】
【図１】第１の実施例のの医用運動解析装置の構成図である。
【図２】時系列速度ベクトルの主成分方向の例である。
【図３】絶対値で表示された速度の例である。
【図４】主成分方向成分の正負区別された速度の例である。
【図５】主成分方向の表示例である。
【図６】収縮中心推定の説明図である。
【図７】第２の実施例のの医用運動解析装置の構成図である。
【図８】運動方向推定に用いる区間を説明する図である。
【図９】歪の様子を説明する図である。
【図１０】変換前の座標系における歪の説明図である。
【図１１】変換後の座標系における歪の説明図である。
【図１２】心筋厚歪分布の表示例である。
【符号の説明】
１０ 医用運動解析装置
１２ 運動情報取得部
１４ 運動情報記憶部
１６ 運動方向推定部
１８ 運動成分算出部
２０ 出力部
Claims (9)
- 心臓を含む臓器を解析対象物として、その解析対象物の運動を解析する医用運動解析装置において、
前記解析対象物の各部分の時系列運動情報を取得する運動情報取得手段と、
前記時系列運動情報を主成分分析して主成分の方向を求め、この主成分の方向が、前記解析対象物の拡張収縮方向を示す解析座標系の主軸であると意味付けする運動方向推定手段と、
前記推定された解析座標系における前記解析対象物の各部分の運動成分を算出する運動成分算出手段と、
を具備する
ことを特徴とする医用運動解析装置。 - 前記時系列運動情報が、前記解析対象物の各部分の時系列速度データである
ことを特徴とする請求項１記載の医用運動解析装置。 - 前記時系列運動情報が、前記解析対象物の各部分の時系列位置データである
ことを特徴とする請求項１記載の医用運動解析装置。 - 前記時系列運動情報が、前記解析対象物の各部分の時系列組織歪データを含み、
前記解析対象物の各部分において前記解析座標系での歪データを算出する歪データ変換手段を具備する
ことを特徴とする請求項１から３のうちいずれか一項に記載の医用運動解析装置。 - 前記時系列運動情報が、前記解析対象物の各部分の時系列組織歪速度データを含む
ことを特徴とする請求項１から３のうちいずれか一項に記載の医用運動解析装置。 - 前記解析対象物が心臓であり、
前記解析座標系における前記心臓の複数の部分の運動成分から前記心臓の収縮中心を推定する収縮中心推定手段を具備する
ことを特徴とする請求項１から５のうちいずれか一項に記載の医用運動解析装置。 - 前記解析対象物が心臓であり、
前記時系列運動情報が、前記心臓の心拍から算出された期間内の時系列運動情報である
ことを特徴とする請求項１から６のうちいずれか一項に記載の医用運動解析装置。 - 心臓を含む臓器を解析対象物として、その解析対象物の運動をコンピュータによって解析する医用運動解析方法において、
前記コンピュータによって取得された前記解析対象物の各部分の時系列運動情報を主成分分析して主成分の方向を求め、この主成分の方向が、前記解析対象物の拡張収縮方向を示す解析座標系の主軸であると前記コンピュータによって意味付けする運動方向推定ステップと、
前記推定された解析座標系における前記解析対象物の各部分の運動成分を前記コンピュータによって算出する運動成分算出ステップと、
を具備する
ことを特徴とする医用運動解析方法。 - 心臓を含む臓器を解析対象物として、その解析対象物の運動を解析するための医用運動解析プログラムであって、コンピュータに、
前記解析対象物の各部分の時系列運動情報を取得する運動情報取得機能と、
前記時系列運動情報を主成分分析して主成分の方向を求め、この主成分の方向が、前記解析対象物の拡張収縮方向を示す解析座標系の主軸であると意味付けする運動方向推定機能と、
前記推定された解析座標系における前記解析対象物の各部分の運動成分を算出する運動成分算出機能と、
を実現させる医用運動解析プログラム。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2002077315A JP3715580B2 (ja) | 2002-03-19 | 2002-03-19 | 医用運動解析装置及びその方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2002077315A JP3715580B2 (ja) | 2002-03-19 | 2002-03-19 | 医用運動解析装置及びその方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2003265480A JP2003265480A (ja) | 2003-09-24 |
| JP3715580B2 true JP3715580B2 (ja) | 2005-11-09 |
Family
ID=29205668
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2002077315A Expired - Fee Related JP3715580B2 (ja) | 2002-03-19 | 2002-03-19 | 医用運動解析装置及びその方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3715580B2 (ja) |
Families Citing this family (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006061581A (ja) | 2004-08-30 | 2006-03-09 | Toshiba Corp | 医用運動解析装置、医用運動解析方法および医用運動解析プログラム |
| JP5063216B2 (ja) * | 2007-06-27 | 2012-10-31 | 日立アロカメディカル株式会社 | 超音波診断装置及び処理プログラム |
| US8315354B2 (en) | 2007-12-07 | 2012-11-20 | Konica Minolta Medical & Graphic, Inc. | Dynamic radiographing system |
| CN102056547B (zh) * | 2008-06-03 | 2014-05-14 | 株式会社日立医疗器械 | 医用图像处理装置及医用图像处理方法 |
| US8208703B2 (en) * | 2008-11-05 | 2012-06-26 | Toshiba Medical Systems Corporation | Medical image analysis apparatus and image analysis control program |
| JP5526401B2 (ja) * | 2009-03-31 | 2014-06-18 | 国立大学法人東京農工大学 | 心室壁情報抽出装置 |
| CN104156975B (zh) * | 2013-05-13 | 2018-04-24 | 东芝医疗系统株式会社 | 医学图像分析装置和方法以及医学成像设备 |
| JP7136588B2 (ja) * | 2018-05-14 | 2022-09-13 | キヤノンメディカルシステムズ株式会社 | 超音波診断装置、医用画像診断装置、医用画像処理装置及び医用画像処理プログラム |
-
2002
- 2002-03-19 JP JP2002077315A patent/JP3715580B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JP2003265480A (ja) | 2003-09-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111432733B (zh) | 用于确定超声探头的运动的设备和方法 | |
| US8391546B2 (en) | Method and corresponding apparatus for quantitative measurements on sequences of images, particularly ultrasonic images | |
| US6295464B1 (en) | Apparatus and method for dynamic modeling of an object | |
| JP4348383B2 (ja) | 生物学的構造体の、湾曲した解剖学的ｍモード表示を行う装置 | |
| JP3725442B2 (ja) | 医用画像診断装置及びその方法 | |
| JP2006061581A (ja) | 医用運動解析装置、医用運動解析方法および医用運動解析プログラム | |
| CN106340015B (zh) | 一种关键点的定位方法和装置 | |
| US20030171668A1 (en) | Image processing apparatus and ultrasonic diagnosis apparatus | |
| Hartkens et al. | Using points and surfaces to improve voxel-based non-rigid registration | |
| JP5491747B2 (ja) | 医用画像解析装置及び画像解析用制御プログラム | |
| JP2003508139A (ja) | 非剛体運動画像解析 | |
| US9508140B2 (en) | Quantifying curvature of biological structures from imaging data | |
| JP2007252904A (ja) | ３ｄボリュームデータセット内での仮想の断層位置決め方法および医用撮像システム | |
| McEachen et al. | Multiframe temporal estimation of cardiac nonrigid motion | |
| US20050148844A1 (en) | Biomagnetic measurement apparatus | |
| JP3715580B2 (ja) | 医用運動解析装置及びその方法 | |
| JPWO2020138136A1 (ja) | 画像処理装置、画像処理方法及びプログラム | |
| CN114581599A (zh) | 一种不依赖外部定位设备的超声三维结构快速获取系统的方法 | |
| Wick et al. | Detection of cardiac quiescence from B-mode echocardiography using a correlation-based frame-to-frame deviation measure | |
| JP2012061019A (ja) | 画像処理装置、方法およびプログラム | |
| CN114255200B (zh) | 基于实时电影成像的心肌应力分析方法及装置 | |
| JP2017080157A (ja) | 画像処理装置 | |
| JP6690053B2 (ja) | 血流解析装置および方法並びにプログラム | |
| Foroughi et al. | Elastic registration of 3D ultrasound images | |
| JP4332669B2 (ja) | 関節中心計測装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20050322 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20050519 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20050607 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20050708 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20050725 |
|
| 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: 20050816 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20050825 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080902 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090902 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090902 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100902 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100902 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110902 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120902 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120902 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130902 Year of fee payment: 8 |
|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313117 Free format text: JAPANESE INTERMEDIATE CODE: R313114 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| LAPS | Cancellation because of no payment of annual fees |