JP3258158B2 - 画像処理方法及びその装置 - Google Patents
画像処理方法及びその装置Info
- Publication number
- JP3258158B2 JP3258158B2 JP33286193A JP33286193A JP3258158B2 JP 3258158 B2 JP3258158 B2 JP 3258158B2 JP 33286193 A JP33286193 A JP 33286193A JP 33286193 A JP33286193 A JP 33286193A JP 3258158 B2 JP3258158 B2 JP 3258158B2
- Authority
- JP
- Japan
- Prior art keywords
- image
- magnetic
- energy
- discrete points
- intersection
- 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
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Description
【０００１】
【産業上の利用分野】本発明は、２次元または３次元の
磁気標識付ＭＲ画像における磁気標識部の検出する画像
処理方法及びその装置に関する。
磁気標識付ＭＲ画像における磁気標識部の検出する画像
処理方法及びその装置に関する。
【０００２】
【従来の技術】近年、虚血性心疾患による死亡者の増加
等により、それらの疾患の診断のために、画像による心
壁の動きの解析が行われるようになってきている。
等により、それらの疾患の診断のために、画像による心
壁の動きの解析が行われるようになってきている。
【０００３】ＭＲＩ装置は、無侵襲に生体内の情報を画
像として得ることができ、Ｘ線ＣＴと違って任意の断面
の画像や３次元画像を直接得ることができる。その上、
近年撮像技術の向上により、現実的な時間で動画像を得
ることが可能となってきている。
像として得ることができ、Ｘ線ＣＴと違って任意の断面
の画像や３次元画像を直接得ることができる。その上、
近年撮像技術の向上により、現実的な時間で動画像を得
ることが可能となってきている。
【０００４】動画像を用いて生体内の組織の動き情報を
得る際には画像間の位置の対応をいかにつけるかが問題
であった。それを解決する方法として、ＭＲＩの場合に
は画像中に磁気標識を入れる方法がある。
得る際には画像間の位置の対応をいかにつけるかが問題
であった。それを解決する方法として、ＭＲＩの場合に
は画像中に磁気標識を入れる方法がある。
【０００５】磁気標識とは、ＭＲＩの撮像技術の一つで
あり、ある組織断面に対して選択的に出力信号の位相を
変化させることにより、画像中に線状または格子状の標
識を入れたものである。磁気標識は組織の移動とともに
移動するため、異なる時相の画像間で対応をつけること
が可能なので、心臓など動態部の動き解析に有効であ
る。図２は、心臓部の２次元磁気標識付ＭＲ画像の例で
ある。６は左心室の心壁、７は放射状に入れた磁気標識
である。
あり、ある組織断面に対して選択的に出力信号の位相を
変化させることにより、画像中に線状または格子状の標
識を入れたものである。磁気標識は組織の移動とともに
移動するため、異なる時相の画像間で対応をつけること
が可能なので、心臓など動態部の動き解析に有効であ
る。図２は、心臓部の２次元磁気標識付ＭＲ画像の例で
ある。６は左心室の心壁、７は放射状に入れた磁気標識
である。
【０００６】従来、磁気標識付ＭＲ画像を用いた動きの
解析は、例えば心臓の形状モデルを利用する方法(A.A.Y
oung,et.al: “Three-dimensional Motion and Deforma
tionof the Heart Wall: Estimation with Spatial Mod
ulation of Magnetization− A Model-based Approach
”, Radiology, vol.185, pp.241-247, 1992) がある
が、磁気標識の交点及び心壁の輪郭の検出は手動で行っ
ている。このため、解析に時間と手間がかかるという問
題点がある。
解析は、例えば心臓の形状モデルを利用する方法(A.A.Y
oung,et.al: “Three-dimensional Motion and Deforma
tionof the Heart Wall: Estimation with Spatial Mod
ulation of Magnetization− A Model-based Approach
”, Radiology, vol.185, pp.241-247, 1992) がある
が、磁気標識の交点及び心壁の輪郭の検出は手動で行っ
ている。このため、解析に時間と手間がかかるという問
題点がある。
【０００７】
【発明が解決しようとする課題】上記のように従来は、
磁気標識の検出は手動で行われていたため、時間と手間
がかかる上、誤差が大きいという問題点があった。
磁気標識の検出は手動で行われていたため、時間と手間
がかかる上、誤差が大きいという問題点があった。
【０００８】そこで、本発明は、磁気標識付ＭＲ画像中
の磁気標識と画像中の対象物の輪郭との交点、または、
前記画像中の複数の磁気標識同士の交点の検出を自動化
することを目的とする画像処理方法及びその装置を提供
する。
の磁気標識と画像中の対象物の輪郭との交点、または、
前記画像中の複数の磁気標識同士の交点の検出を自動化
することを目的とする画像処理方法及びその装置を提供
する。
【０００９】
【課題を解決するための手段】第１の発明の画像処理方
法は、磁気標識付ＭＲ画像中の磁気標識と対象物の輪郭
との交点、または、前記画像中の複数の磁気標識同士の
交点の検出方法であって、磁気標識付ＭＲ画像を入力
し、前記画像中に複数個の離散点を配し、前記複数個の
離散点が有する弾性エネルギーを求め、前記複数個の離
散点が有する磁気標識部分の特徴と磁気標識以外の部分
の特徴によって決定される画像エネルギーを求め、少な
くとも前記弾性エネルギーと前記画像エネルギーとの総
和を最小にするように前記複数個の離散点を収束するま
で移動させ、収束した離散点の位置により前記交点の位
置を検出するものである。
法は、磁気標識付ＭＲ画像中の磁気標識と対象物の輪郭
との交点、または、前記画像中の複数の磁気標識同士の
交点の検出方法であって、磁気標識付ＭＲ画像を入力
し、前記画像中に複数個の離散点を配し、前記複数個の
離散点が有する弾性エネルギーを求め、前記複数個の離
散点が有する磁気標識部分の特徴と磁気標識以外の部分
の特徴によって決定される画像エネルギーを求め、少な
くとも前記弾性エネルギーと前記画像エネルギーとの総
和を最小にするように前記複数個の離散点を収束するま
で移動させ、収束した離散点の位置により前記交点の位
置を検出するものである。
【００１０】第２の発明の画像処理方法は、磁気標識付
ＭＲ画像中の磁気標識と対象物の輪郭との交点、また
は、前記画像中の複数の磁気標識同士の交点の検出方法
であって、磁気標識付ＭＲ画像を入力し、前記画像中に
複数個の離散点を配し、前記複数個の離散点が有する弾
性エネルギーを求め、前記複数個の離散点の内、一部の
離散点が有する磁気標識部分の特徴と磁気標識以外の部
分の特徴によって決定されるエネルギー及びその他の離
散点が有する磁気標識以外の部分の特徴によって決定さ
れるエネルギーとよりなる画像エネルギーを求め、少な
くとも前記弾性エネルギーと前記画像エネルギーとの総
和を最小にするように前記複数個の離散点を収束するま
で移動させ、収束した離散点の位置により前記交点の位
置を検出するものである。
ＭＲ画像中の磁気標識と対象物の輪郭との交点、また
は、前記画像中の複数の磁気標識同士の交点の検出方法
であって、磁気標識付ＭＲ画像を入力し、前記画像中に
複数個の離散点を配し、前記複数個の離散点が有する弾
性エネルギーを求め、前記複数個の離散点の内、一部の
離散点が有する磁気標識部分の特徴と磁気標識以外の部
分の特徴によって決定されるエネルギー及びその他の離
散点が有する磁気標識以外の部分の特徴によって決定さ
れるエネルギーとよりなる画像エネルギーを求め、少な
くとも前記弾性エネルギーと前記画像エネルギーとの総
和を最小にするように前記複数個の離散点を収束するま
で移動させ、収束した離散点の位置により前記交点の位
置を検出するものである。
【００１１】
【作 用】上記第１の発明は、画像処理方法における磁
気標識付ＭＲ画像中の磁気標識と対象物の輪郭との交
点、または、前記画像中の複数の磁気標識同士の交点の
検出方法について説明する。
気標識付ＭＲ画像中の磁気標識と対象物の輪郭との交
点、または、前記画像中の複数の磁気標識同士の交点の
検出方法について説明する。
【００１２】磁気標識付ＭＲ画像を入力する。
【００１３】前記入力した画像中に複数個の離散点を配
する。
する。
【００１４】前記複数個の離散点が有する弾性エネルギ
ーを求める。
ーを求める。
【００１５】前記全ての離散点が有する磁気標識部分の
特徴と磁気標識以外の部分の特徴によって決定される画
像エネルギーを求める。
特徴と磁気標識以外の部分の特徴によって決定される画
像エネルギーを求める。
【００１６】少くとも前記弾性エネルギーと前記画像エ
ネルギーとの総和を最小にするように前記複数個の離散
点を収束するまで移動させる。
ネルギーとの総和を最小にするように前記複数個の離散
点を収束するまで移動させる。
【００１７】収束した離散点の位置により前記交点の位
置を検出するなお、使用する離散点は、画像中に配した
複数個の離散点を全て使用してもよいが、その一部分の
みの離散点を使用してもよい。
置を検出するなお、使用する離散点は、画像中に配した
複数個の離散点を全て使用してもよいが、その一部分の
みの離散点を使用してもよい。
【００１８】第２の発明の画像処理方法は、第１の発明
における画像エネルギーに代えて、前記複数個の離散点
の内、一部の離散点が有する磁気標識部分の特徴と磁気
標識以外の部分の特徴によって決定されるエネルギー及
びその他の離散点が有する磁気標識以外の部分の特徴に
よって決定されるエネルギーとよりなる画像エネルギー
を用いる。
における画像エネルギーに代えて、前記複数個の離散点
の内、一部の離散点が有する磁気標識部分の特徴と磁気
標識以外の部分の特徴によって決定されるエネルギー及
びその他の離散点が有する磁気標識以外の部分の特徴に
よって決定されるエネルギーとよりなる画像エネルギー
を用いる。
【００１９】
【実施例】第 １ の 実 施 例 以下、本発明の第１の実施例を図１〜８を用いて説明す
る。
る。
【００２０】［磁気標識検出の対象画像が、２次元デー
タの場合］図１は、本実施例のブロック図とフローチャ
ートを兼ねた図である。この図を用いて図２のような放
射状の磁気標識の入った心臓部の２次元ＭＲ画像から、
磁気標識と心壁の輪郭との交点を求める場合について説
明する。
タの場合］図１は、本実施例のブロック図とフローチャ
ートを兼ねた図である。この図を用いて図２のような放
射状の磁気標識の入った心臓部の２次元ＭＲ画像から、
磁気標識と心壁の輪郭との交点を求める場合について説
明する。
【００２１】（１）画像データ入力部１ 画像データ入力部１において、対象とする磁気標識付Ｍ
Ｒ画像のデータを入力する。
Ｒ画像のデータを入力する。
【００２２】（２）初期値設定部２ 初期値設定部２において、複数個の離散点の位置座標の
初期値を設定する。初期値の例を図３に示す。
初期値を設定する。初期値の例を図３に示す。
【００２３】（３）離散点移動部３ 離散点移動部３において、［数１］においてＥall で表
されるような全エネルギーが最小になるように前記複数
個の離散点を移動させる。
されるような全エネルギーが最小になるように前記複数
個の離散点を移動させる。
【００２４】
【数１】 ｉは、図５で示す複数個の離散点の番号である。ｗtag
（ｉ），ｗext （ｉ）は、重み係数である。
（ｉ），ｗext （ｉ）は、重み係数である。
【００２５】Ｅint （ｉ）は、弾性エネルギーであり、
本実施例では、次式を用いる。
本実施例では、次式を用いる。
【００２６】
【数２】 但し、ｖi （ｘｉ，ｙｉ）は離散点の位置座標である。
【００２７】Ｅtag （ｉ）は、磁気標識の特徴を表す画
像エネルギーであって、画像濃度と、画像のエッジ強度
を表すエネルギーの重みの和を用いる。すなわち、
像エネルギーであって、画像濃度と、画像のエッジ強度
を表すエネルギーの重みの和を用いる。すなわち、
【数３】 となる。
【００２８】但し、Ｉimage （ｖi ）は、画像濃度であ
り、Ｅedge（ｉ）は、エッジ強度を表すエネルギーであ
る。
り、Ｅedge（ｉ）は、エッジ強度を表すエネルギーであ
る。
【００２９】ここで、磁気標識の特徴を表す画像エネル
ギーとして上式を用いたのは、磁気標識部分は画像中の
他の部分と比べて画像濃度が急激に変化して、磁気標識
と輪郭の交点ではエッジ強度が大きくなっているからで
ある。
ギーとして上式を用いたのは、磁気標識部分は画像中の
他の部分と比べて画像濃度が急激に変化して、磁気標識
と輪郭の交点ではエッジ強度が大きくなっているからで
ある。
【００３０】Ｅedge（ｉ）としては、本実施例では、
【数４】
【００３１】を用いる。
【００３２】Ｅext （ｉ）は必要に応じて付加する外部
エネルギーである。本実施例では［数５］は離散点をあ
る点ｖ0 から遠ざけたり、ｖ0 に近付けたりするエネル
ギーである。遠ざけるか近付けるかは係数ｋの正負によ
って決まる。
エネルギーである。本実施例では［数５］は離散点をあ
る点ｖ0 から遠ざけたり、ｖ0 に近付けたりするエネル
ギーである。遠ざけるか近付けるかは係数ｋの正負によ
って決まる。
【００３３】
【数５】
【００３４】エネルギーの総和が最小となるような離散
点の位置を求める解法として、例えば変分法による方法
を用いる。変分法によれば［数１］が最小となるために
は［数６］のような連立方程式が成り立つことが必要で
ある。
点の位置を求める解法として、例えば変分法による方法
を用いる。変分法によれば［数１］が最小となるために
は［数６］のような連立方程式が成り立つことが必要で
ある。
【００３５】
【数６】 但し、各項の重みは定数とし、
【数７】 とする。
【００３６】今、ｘ座標のみを考えるとすると、［数
６］の連立方程式はベクトル表示で次式のように表すこ
とができる。
６］の連立方程式はベクトル表示で次式のように表すこ
とができる。
【００３７】
【数８】 但し、
【数９】 となる。但し、ａ＝β，ｂ＝−α−４β，ｃ＝２α＋６
β，ｆi ＝−ｆ（ｖi ）である。
β，ｆi ＝−ｆ（ｖi ）である。
【００３８】この連立方程式をヤコビ法を用いて逐次近
似法で解く。Ａを次式のように対角行列Ｄとそれ以外の
部分Ｆに分ける。
似法で解く。Ａを次式のように対角行列Ｄとそれ以外の
部分Ｆに分ける。
【００３９】
【数１０】 ｎ回目の近似解をｘｎとすると、ｎ＋１回目の近似解ｘ
n+1 は、
n+1 は、
【数１１】 となる。但し、γは収束速度パラメータである。
【００４０】離散点移動部３の構成例を図４に示す。
【００４１】画像エネルギー算出部９、外部エネルギー
算出部１０において、それぞれ［数１］におけるＥtag
（ｉ）、Ｅext （ｉ）を算出する。
算出部１０において、それぞれ［数１］におけるＥtag
（ｉ）、Ｅext （ｉ）を算出する。
【００４２】パラメータ記憶部１２において、各エネル
ギーの重み係数、［数１１］における行列Ｄ、Ｆの要素
の値、γ等が記憶されている。
ギーの重み係数、［数１１］における行列Ｄ、Ｆの要素
の値、γ等が記憶されている。
【００４３】次座標算出部１１において、［数１１］を
用いて各離散点の次座標を計算して、離散点を移動させ
る。
用いて各離散点の次座標を計算して、離散点を移動させ
る。
【００４４】（４）収束判定部４ 収束判定部４において、離散点を移動させるごとに収束
条件を満たすかどうかを判定し、収束条件を満たさない
場合には再び離散点移動部３で離散点の移動を行ない、
収束条件を満たすまで繰り返す。
条件を満たすかどうかを判定し、収束条件を満たさない
場合には再び離散点移動部３で離散点の移動を行ない、
収束条件を満たすまで繰り返す。
【００４５】収束条件としては、本実施例では［数１
２］で表されるようなｎ回目の繰り返しでのエネルギー
の変化量ΔＥall n が、ある一定量より小さくなったと
きとする。
２］で表されるようなｎ回目の繰り返しでのエネルギー
の変化量ΔＥall n が、ある一定量より小さくなったと
きとする。
【００４６】
【数１２】 収束判定部４の一構成例を図６に示す。
【００４７】エネルギー算出部１３において、［数１］
におけるＥall を算出する。それをｎ回目の値Ｅall n
として、記憶部１４に記憶しておく。
におけるＥall を算出する。それをｎ回目の値Ｅall n
として、記憶部１４に記憶しておく。
【００４８】そして、エネルギー比較部１５で、記憶部
１４に予め記憶しておいたｎ−１回目の値Ｅall n-1 と
の差の絶対値ΔＥall n を求めて、それが一定量εより
小さければ収束したと判定し、大きければ収束していな
いと判定する。ここで収束条件を満たした収束状態は、
各離散点が、全て磁気標識と心壁の輪郭との交点に集ま
る。したがって、この交点の座標は、各離散点の移動後
の座標から容易に求められる。
１４に予め記憶しておいたｎ−１回目の値Ｅall n-1 と
の差の絶対値ΔＥall n を求めて、それが一定量εより
小さければ収束したと判定し、大きければ収束していな
いと判定する。ここで収束条件を満たした収束状態は、
各離散点が、全て磁気標識と心壁の輪郭との交点に集ま
る。したがって、この交点の座標は、各離散点の移動後
の座標から容易に求められる。
【００４９】（５）出力部５ 上記結果を出力部５から出力する。結果の出力方法とし
ては、データとしてファイルに記録したり、ディスプレ
イに表示するなどの方法があるが、本実施例ではディス
プレイに表示する。
ては、データとしてファイルに記録したり、ディスプレ
イに表示するなどの方法があるが、本実施例ではディス
プレイに表示する。
【００５０】［磁気標識検出の対象画像が、３次元デー
タの場合］次に、磁気標識検出の対象画像がボクセルデ
ータまたは複数枚のスライスデータからなる３次元デー
タである場合について説明する。
タの場合］次に、磁気標識検出の対象画像がボクセルデ
ータまたは複数枚のスライスデータからなる３次元デー
タである場合について説明する。
【００５１】３次元データを対象とする場合には、図５
のような離散点の代わりに例えば図７のような格子の上
に離散点を配置して、それを図８のように閉輪郭化した
ものを用いる。この場合、格子間隔は一定である必要は
ない。
のような離散点の代わりに例えば図７のような格子の上
に離散点を配置して、それを図８のように閉輪郭化した
ものを用いる。この場合、格子間隔は一定である必要は
ない。
【００５２】また、全エネルギーは［数１］と同様の次
式のようなものを用いる。
式のようなものを用いる。
【００５３】
【数１３】 ｗtag （ｉ，ｊ）、ｗext （ｉ，ｊ）は重み係数、ｉ，
ｊは、各離散点における図７のような格子上での座標を
表す。
ｊは、各離散点における図７のような格子上での座標を
表す。
【００５４】また、弾性エネルギーＥint （ｉ，ｊ）と
しては次式を用いる。
しては次式を用いる。
【００５５】
【数１４】 但し、ｖi,j ＝（ｘi,j ，ｙi,j ，ｚi,j ）である。磁
気標識の特徴を表すエネルギーＥtag （ｉ，ｊ）として
は［数３］と同様の次式を用いる。
気標識の特徴を表すエネルギーＥtag （ｉ，ｊ）として
は［数３］と同様の次式を用いる。
【００５６】
【数１５】 Ｉimage （ｖi,j ）は、画像濃度であり、Ｅ
edge（ｉ，ｊ）は、エッジ強度を表すエネルギーであ
る。
edge（ｉ，ｊ）は、エッジ強度を表すエネルギーであ
る。
【００５７】Ｅedge（ｉ，ｊ）としては、次式を用い
る。
る。
【００５８】
【数１６】 Ｅext （ｉ，ｊ）としても［数５］と同様に次式を用い
る。
る。
【００５９】
【数１７】 全エネルギーとして［数１］の代わりに［数１３］を用
いる以外は２次元の場合と同様に、エネルギーの総和を
最小になるように離散点を移動させる処理を行うことに
より、磁気標識を検出する。
いる以外は２次元の場合と同様に、エネルギーの総和を
最小になるように離散点を移動させる処理を行うことに
より、磁気標識を検出する。
【００６０】第 ２ の 実 施 例 次に、本発明の第２の実施例について図９〜１３を用い
て説明する。
て説明する。
【００６１】［磁気標識検出の対象画像が、２次元デー
タの場合］図９は、磁気標識検出の対象とする画像また
は対象画像と類似の画像から、画像中の対象物の輪郭を
予め抽出し、その結果を磁気標識の検出に利用する場合
のブロック図とフローチャートを兼ねた図である。この
図により図２のような心臓部の２次元の磁気標識付ＭＲ
画像から、磁気標識と心壁の輪郭との交点を求める場合
について説明する。
タの場合］図９は、磁気標識検出の対象とする画像また
は対象画像と類似の画像から、画像中の対象物の輪郭を
予め抽出し、その結果を磁気標識の検出に利用する場合
のブロック図とフローチャートを兼ねた図である。この
図により図２のような心臓部の２次元の磁気標識付ＭＲ
画像から、磁気標識と心壁の輪郭との交点を求める場合
について説明する。
【００６２】磁気標識検出の対象画像と類似の画像とし
ては、動画像では対象画像と同時相で磁気標識の入って
いない画像、対象画像と近い時相で磁気標識の入ってい
ない画像のいずれかを用いる。
ては、動画像では対象画像と同時相で磁気標識の入って
いない画像、対象画像と近い時相で磁気標識の入ってい
ない画像のいずれかを用いる。
【００６３】（１）輪郭抽出部１６ はじめに輪郭抽出部１６で対象物の輪郭を抽出する。第
１の実施例における磁気標識検出と同様にエネルギー最
小化によって抽出する場合の構成例を図１０に示す。
１の実施例における磁気標識検出と同様にエネルギー最
小化によって抽出する場合の構成例を図１０に示す。
【００６４】まず，画像データ入力部１で輪郭抽出を行
う画像を入力し、次の初期値設定部２で複数個の離散点
の初期座標の値を設定する。
う画像を入力し、次の初期値設定部２で複数個の離散点
の初期座標の値を設定する。
【００６５】そして、離散点移動部３で第１の実施例に
おける磁気標識検出の場合と同様に、エネルギーの総和
が最小になるように離散点を移動させる処理を行う。
おける磁気標識検出の場合と同様に、エネルギーの総和
が最小になるように離散点を移動させる処理を行う。
【００６６】但し、全エネルギーとして次式のＥall を
用いる。
用いる。
【００６７】
【数１８】
【００６８】但し、ｗedge（ｉ）、ｗext （ｉ）は重み
係数である。Ｅint （ｉ）は弾性エネルギーであり、本
実施例では［数２］を用いる。Ｅext （ｉ）は、必要に
応じて付加する外部エネルギーで、本実施例では［数
５］を用いる。Ｅedge（ｉ）は画像エネルギーであり、
本実施例では［数４］のエッジ強度を表すエネルギーを
用いる。
係数である。Ｅint （ｉ）は弾性エネルギーであり、本
実施例では［数２］を用いる。Ｅext （ｉ）は、必要に
応じて付加する外部エネルギーで、本実施例では［数
５］を用いる。Ｅedge（ｉ）は画像エネルギーであり、
本実施例では［数４］のエッジ強度を表すエネルギーを
用いる。
【００６９】（２）離散点配置部１７ 輪郭を抽出した後、離散点配置部１７で、抽出した輪郭
上に複数個の離散点を配置する。但し、前述したような
エネルギー最小化で輪郭抽出を行なった場合には、すで
に離散点が配置されているため、新たに離散点を配置す
る必要はないが、必要に応じて離散点の並べ変えを行
う。具体的には、隣接する離散点が重なってしまった場
合には一方を取り除く、隣接する離散点間の距離が離れ
過ぎた場合には補間するといった処理を行う。
上に複数個の離散点を配置する。但し、前述したような
エネルギー最小化で輪郭抽出を行なった場合には、すで
に離散点が配置されているため、新たに離散点を配置す
る必要はないが、必要に応じて離散点の並べ変えを行
う。具体的には、隣接する離散点が重なってしまった場
合には一方を取り除く、隣接する離散点間の距離が離れ
過ぎた場合には補間するといった処理を行う。
【００７０】（３）磁気標識対応点決定部１８ 次に、それらの配置された離散点の中から、磁気標識対
応点決定部１８で磁気標識対応点を決定する。配置され
た離散点は輪郭線上にあるが、それら離散点のうち一部
の点を磁気標識に対応させて、それ以外の点は輪郭線に
対応するようにすれば、それら複数の離散点は輪郭から
外れることなく、磁気標識部を検出することができる。
このような、磁気標識に対応させる点のことを磁気標識
対応点と呼ぶことにする。
応点決定部１８で磁気標識対応点を決定する。配置され
た離散点は輪郭線上にあるが、それら離散点のうち一部
の点を磁気標識に対応させて、それ以外の点は輪郭線に
対応するようにすれば、それら複数の離散点は輪郭から
外れることなく、磁気標識部を検出することができる。
このような、磁気標識に対応させる点のことを磁気標識
対応点と呼ぶことにする。
【００７１】磁気標識対応点の決定方法としては、図２
に模式的に示す画像が主であるため、磁気標識部分７の
画像濃度が低いので、各点の近傍を輪郭線上に沿って探
索し、画像濃度が低いところが多い場合にはそれを磁気
標識対応点とするという方法を採用する。なお、画像濃
度が低いとは、画像上明るさが暗い箇所をいう。本手法
により、輪郭抽出後の離散点が図１１のようになった場
合には、図１２の×に示す磁気標識対応点２０を決定す
ることができる。この磁気標識対応点２０は、求めたい
磁気標識と心壁の輪郭との交点の近傍には位置している
が、正確な位置ではないため、さらに磁気標識検出部１
９で正確な位置を決定する。
に模式的に示す画像が主であるため、磁気標識部分７の
画像濃度が低いので、各点の近傍を輪郭線上に沿って探
索し、画像濃度が低いところが多い場合にはそれを磁気
標識対応点とするという方法を採用する。なお、画像濃
度が低いとは、画像上明るさが暗い箇所をいう。本手法
により、輪郭抽出後の離散点が図１１のようになった場
合には、図１２の×に示す磁気標識対応点２０を決定す
ることができる。この磁気標識対応点２０は、求めたい
磁気標識と心壁の輪郭との交点の近傍には位置している
が、正確な位置ではないため、さらに磁気標識検出部１
９で正確な位置を決定する。
【００７２】（４）磁気標識検出部１９ 磁気標識対応点が決定された後は、磁気標識検出部１９
で磁気標識と心壁の輪郭との交点を正確に検出する。
で磁気標識と心壁の輪郭との交点を正確に検出する。
【００７３】磁気標識検出部の一構成例を図１３に示
す。
す。
【００７４】これは、図１０における離散点移動部３以
降の処理と同様である。
降の処理と同様である。
【００７５】まず、磁気標識検出の対象画像が、輪郭抽
出した画像と同一でない場合には画像データ入力部１で
同一の画像データを入力する。同一の画像データの場合
は、輪郭抽出部１６で入力した画像データを用いる。
出した画像と同一でない場合には画像データ入力部１で
同一の画像データを入力する。同一の画像データの場合
は、輪郭抽出部１６で入力した画像データを用いる。
【００７６】離散点移動部３で離散点を移動させる処理
を行う。但し、画像エネルギーとしては、磁気標識対応
点に対しては［数１］を特別に与え、その他の点に対し
ては、［数１８］を用いる。すなわち全エネルギーは次
式のように表わされる。
を行う。但し、画像エネルギーとしては、磁気標識対応
点に対しては［数１］を特別に与え、その他の点に対し
ては、［数１８］を用いる。すなわち全エネルギーは次
式のように表わされる。
【００７７】
【数１９】
【００７８】但し、ｗtag （ｉ）、ｗext （ｉ）、ｗ
（ｉ）、ｗedge（ｉ）、ｗ' ext （ｉ）は重み係数であ
る 収束判定部４で収束したと判定したら、出力部５で結果
を出力する。ここで収束条件を満たした収束状態は、各
磁気標識対応点が、全て磁気標識と心壁の輪郭との交点
に正確に集まる。したがって、この交点の座標は、各磁
気標識対応点の移動後の座標から容易に求められる。
（ｉ）、ｗedge（ｉ）、ｗ' ext （ｉ）は重み係数であ
る 収束判定部４で収束したと判定したら、出力部５で結果
を出力する。ここで収束条件を満たした収束状態は、各
磁気標識対応点が、全て磁気標識と心壁の輪郭との交点
に正確に集まる。したがって、この交点の座標は、各磁
気標識対応点の移動後の座標から容易に求められる。
【００７９】このように輪郭抽出結果を磁気標識検出に
利用することにより、画像中の対象物の輪郭線と磁気標
識の交点を精度よく検出することができる。
利用することにより、画像中の対象物の輪郭線と磁気標
識の交点を精度よく検出することができる。
【００８０】［磁気標識検出の対象画像が、３次元デー
タの場合］磁気標識検出の対象画像がボクセルデータま
たは複数枚のスライスデータからなる３次元データであ
る場合についても図７中の格子点上に離散点の配置を行
ない、輪郭抽出の際の全エネルギーを［数１８］の代わ
りに、
タの場合］磁気標識検出の対象画像がボクセルデータま
たは複数枚のスライスデータからなる３次元データであ
る場合についても図７中の格子点上に離散点の配置を行
ない、輪郭抽出の際の全エネルギーを［数１８］の代わ
りに、
【数２０】 を用い、磁気標識検出の際の全エネルギーを［数１９］
の代わりに、
の代わりに、
【数２１】
【００８１】を用いればよい。但し、Ｅint （ｉ，ｊ）
としては［数１４］、Ｅtag （ｉ，ｊ）としては［数１
５］、Ｅext （ｉ，ｊ）としては［数１７］、Ｅ
edge（ｉ，ｊ）としては［数１６］を用いる。
としては［数１４］、Ｅtag （ｉ，ｊ）としては［数１
５］、Ｅext （ｉ，ｊ）としては［数１７］、Ｅ
edge（ｉ，ｊ）としては［数１６］を用いる。
【００８２】第 ３ の 実 施 例 次に、本発明の第３の実施例について図１４を用いて説
明する。
明する。
【００８３】［磁気標識検出の対象画像が、２次元デー
タの場合］図１４は、複数時相の画像を用いて動きの検
出を行う場合のうち、特に輪郭抽出を予め行う場合のブ
ロック図とフローチャートを兼ねた図である。
タの場合］図１４は、複数時相の画像を用いて動きの検
出を行う場合のうち、特に輪郭抽出を予め行う場合のブ
ロック図とフローチャートを兼ねた図である。
【００８４】初期画像において磁気標識検出部１９で磁
気標識検出を行うまでの処理は、第２の実施例の図９と
同様である。まず、輪郭抽出部１６で輪郭を抽出する。
次に離散点配置部１７で、離散点を配置する。磁気標識
対応点決定部１８で磁気標識対応点を決定し、磁気標識
検出部１９で磁気標識を検出する。
気標識検出を行うまでの処理は、第２の実施例の図９と
同様である。まず、輪郭抽出部１６で輪郭を抽出する。
次に離散点配置部１７で、離散点を配置する。磁気標識
対応点決定部１８で磁気標識対応点を決定し、磁気標識
検出部１９で磁気標識を検出する。
【００８５】（１）動き量測定部２１ 次に、動き量測定部２１において、磁気標識検出後の各
磁気標識対応点の座標とそれに対応する前画像における
磁気標識対応点の座標との差を動き量として求める。
磁気標識対応点の座標とそれに対応する前画像における
磁気標識対応点の座標との差を動き量として求める。
【００８６】ｌ番目の画像における磁気標識対応点の点
列をｖl 、ｌ＋１番目の画像における磁気標識対応点の
点列をｖl+1 とすると、ｌ番目の画像からｌ＋１番目の
画像への各磁気標識対応点の動きベクトルである移動量
Δｖl は、［数２２］で表される。
列をｖl 、ｌ＋１番目の画像における磁気標識対応点の
点列をｖl+1 とすると、ｌ番目の画像からｌ＋１番目の
画像への各磁気標識対応点の動きベクトルである移動量
Δｖl は、［数２２］で表される。
【００８７】
【数２２】 求めた移動量Δｖl は、記憶部２２でメモリー等に記憶
しておく。
しておく。
【００８８】（２）終了判定部２３ 終了判定部２３において、画像がまだあるかどうかを判
定し、画像がまだある場合には磁気標識検出部１９に戻
って次の画像の処理を行う。画像がそれ以上ない場合に
は出力部５で結果を出力して終了する。ここで出力の方
法としては、移動量Δｖl を動きベクトルとして表示す
るか、数値化して表示する方法が考えられる。
定し、画像がまだある場合には磁気標識検出部１９に戻
って次の画像の処理を行う。画像がそれ以上ない場合に
は出力部５で結果を出力して終了する。ここで出力の方
法としては、移動量Δｖl を動きベクトルとして表示す
るか、数値化して表示する方法が考えられる。
【００８９】これらの処理を繰り返すことにより、複数
時相の画像から動きの検出を定量的に行うことができ
る。
時相の画像から動きの検出を定量的に行うことができ
る。
【００９０】［磁気標識検出の対象画像が、３次元デー
タの場合］磁気標識検出の対象画像が、ボクセルデータ
または複数枚のスライスデータからなる３次元データで
ある場合についても、例えば、図７のような格子上に離
散点の配置を行い、エネルギーの式として第１、第２の
実施例で示したように３次元用の式を用いれば同様に磁
気標識の検出及び動き量の測定を行うことができる。
タの場合］磁気標識検出の対象画像が、ボクセルデータ
または複数枚のスライスデータからなる３次元データで
ある場合についても、例えば、図７のような格子上に離
散点の配置を行い、エネルギーの式として第１、第２の
実施例で示したように３次元用の式を用いれば同様に磁
気標識の検出及び動き量の測定を行うことができる。
【００９１】第 ４ の 実 施 例 次に、本発明の第４の実施例について図１５〜１９を用
いて説明する。
いて説明する。
【００９２】本実施例は、２次元の磁気標識付ＭＲ画像
において複数の磁気標識の交点を検出する例である。複
数の磁気標識間の交点とは、図１５の２５に示すもので
ある。
において複数の磁気標識の交点を検出する例である。複
数の磁気標識間の交点とは、図１５の２５に示すもので
ある。
【００９３】図１６は、このような複数の磁気標識間の
交点を検出する場合のブロック図とフローチャートを兼
ねた図である。
交点を検出する場合のブロック図とフローチャートを兼
ねた図である。
【００９４】（１）画像データ入力部１ はじめに、画像データ入力部１で対象とする磁気標識付
ＭＲ画像のデータを入力する。
ＭＲ画像のデータを入力する。
【００９５】（２）初期値設定部２ 次に、初期値設定部２で複数個の離散点の位置座標の初
期値を設定する。初期値の例を図１７に示す。
期値を設定する。初期値の例を図１７に示す。
【００９６】初期値の設定方法としては、第２の実施例
で示したようなエネルギー最小化などで心壁輪郭を抽出
した後、磁気標識対応点を決定してそれらの磁気標識対
応点を初期値として用いる。
で示したようなエネルギー最小化などで心壁輪郭を抽出
した後、磁気標識対応点を決定してそれらの磁気標識対
応点を初期値として用いる。
【００９７】（３）離散点移動部３ 次に、離散点移動部３で全エネルギーが最小になるよう
に離散点を移動させる。全エネルギーとしては、［数
１］で表されるものを使用する。但し、Ｅint （ｉ）と
しては［数２］、Ｅtag （ｉ）としては［数３］を用い
る。図１７において、初期値より外側の交点を検出する
には各離散点が外側に行くほどそのエネルギーが低くな
るような外部エネルギーＥext （ｉ）を用いる。本実施
例では［数５］でｖ0 を［数２３］で表わされるような
全離散点の重心とし、係数ｋを負の値としたものを用い
る。
に離散点を移動させる。全エネルギーとしては、［数
１］で表されるものを使用する。但し、Ｅint （ｉ）と
しては［数２］、Ｅtag （ｉ）としては［数３］を用い
る。図１７において、初期値より外側の交点を検出する
には各離散点が外側に行くほどそのエネルギーが低くな
るような外部エネルギーＥext （ｉ）を用いる。本実施
例では［数５］でｖ0 を［数２３］で表わされるような
全離散点の重心とし、係数ｋを負の値としたものを用い
る。
【００９８】
【数２３】 但し、Ｎは、離散点の数を示す。
【００９９】この場合、ｋの値は離散点を収縮させるＥ
int （ｉ）より膨張させるＥext （ｉ）の方が大きくな
るように調整する。
int （ｉ）より膨張させるＥext （ｉ）の方が大きくな
るように調整する。
【０１００】（４）収束判定部４ そして、収束判定部４で収束を判定する。
【０１０１】磁気標識部分は画像中の他の部分と比べて
画像濃度が急激に変化しているため、磁気標識の境界部
分ではエッジ強度が大きくなっている。また、磁気標識
の境界部分に沿ってはエッジ強度があまり変化しないた
めに、［数３］の変化量が少ない。
画像濃度が急激に変化しているため、磁気標識の境界部
分ではエッジ強度が大きくなっている。また、磁気標識
の境界部分に沿ってはエッジ強度があまり変化しないた
めに、［数３］の変化量が少ない。
【０１０２】したがって、［数１］においてＥ
ext （ｉ）の項が支配的な状態となる。このために離散
点移動部３において各離散点は磁気標識の境界部分に沿
って外側に移動する。
ext （ｉ）の項が支配的な状態となる。このために離散
点移動部３において各離散点は磁気標識の境界部分に沿
って外側に移動する。
【０１０３】ところが、離散点が複数の磁気標識の交点
に達すると、［数３］が一時的に低くなり、［数１］で
表される全エネルギーが極小値をとる。そのため、［数
１２］のΔＥall n で表されるような全エネルギーの変
化量や離散点の移動量がある一定値より小さくなるとい
う条件や、繰り返しの回数によって収束させることがで
きる。図１８に離散点の動きの例を示す。矢印の方向が
全エネルギーの低くなる方向すなわち初期値に比べて外
側を示している。前記の収束条件により、離散点は磁気
標識の境界部分に沿ってａ、ｂと移動するが、交点に達
するとｃで止まる。
に達すると、［数３］が一時的に低くなり、［数１］で
表される全エネルギーが極小値をとる。そのため、［数
１２］のΔＥall n で表されるような全エネルギーの変
化量や離散点の移動量がある一定値より小さくなるとい
う条件や、繰り返しの回数によって収束させることがで
きる。図１８に離散点の動きの例を示す。矢印の方向が
全エネルギーの低くなる方向すなわち初期値に比べて外
側を示している。前記の収束条件により、離散点は磁気
標識の境界部分に沿ってａ、ｂと移動するが、交点に達
するとｃで止まる。
【０１０４】収束条件を満たすまで離散点の移動を繰り
返す。
返す。
【０１０５】（５）終了判定部２６ 収束条件を満たしたら、今度は終了判定部２６で、終了
するかどうかを判定する。
するかどうかを判定する。
【０１０６】終了の判定は、検出する複数の磁気標識間
の交点の数を予め設定しておき、その数に達したら終了
する。あるいは離散点が、画像の外に出たら終了すると
する。
の交点の数を予め設定しておき、その数に達したら終了
する。あるいは離散点が、画像の外に出たら終了すると
する。
【０１０７】このように収束計算を複数回繰り返すこと
で、複数個の交点を検出することが可能となる。終了判
定部２６で終了した場合には出力部５で結果を出力す
る。
で、複数個の交点を検出することが可能となる。終了判
定部２６で終了した場合には出力部５で結果を出力す
る。
【０１０８】なお、本発明はここに挙げた実施例に限定
されるものではない。
されるものではない。
【０１０９】第１の実施例において、収束判定部４にお
ける収束条件としては、エネルギーの変化量による方法
を挙げたが、複数個の離散点の移動量がある一定値より
小さくなったとき、あるいはこのような条件で収束しな
い場合もあることを考慮して、繰り返し回数で収束させ
るようにすることもできる。
ける収束条件としては、エネルギーの変化量による方法
を挙げたが、複数個の離散点の移動量がある一定値より
小さくなったとき、あるいはこのような条件で収束しな
い場合もあることを考慮して、繰り返し回数で収束させ
るようにすることもできる。
【０１１０】第２の実施例では、輪郭抽出部１６で、輪
郭抽出方法としてエネルギー最小化による方法の例を示
したが、手で輪郭抽出した結果を用いても差し支えな
い。これは第３の実施例でも同様である。
郭抽出方法としてエネルギー最小化による方法の例を示
したが、手で輪郭抽出した結果を用いても差し支えな
い。これは第３の実施例でも同様である。
【０１１１】第３の実施例では、磁気標識検出部１９に
おける磁気標識検出の方法としては、磁気標識対応点の
みを用いる方法、すなわち、画像エネルギーとしては磁
気標識対応点が磁気標識部分に引き付けられるようなエ
ネルギーのみを用いる方法もある。この場合には全エネ
ルギーとして［数１］を用いる。
おける磁気標識検出の方法としては、磁気標識対応点の
みを用いる方法、すなわち、画像エネルギーとしては磁
気標識対応点が磁気標識部分に引き付けられるようなエ
ネルギーのみを用いる方法もある。この場合には全エネ
ルギーとして［数１］を用いる。
【０１１２】第４の実施例において、複数の磁気標識間
の交点を検出する場合、初期値は手で設定しても差し支
えない。
の交点を検出する場合、初期値は手で設定しても差し支
えない。
【０１１３】また、収束判定部４における収束条件とし
て、画像のエッジ方向の変化量がある一定値よりも大き
くなることを収束条件とすることもできる。次式は勾配
によって求めたエッジ方向の変化量の例である。
て、画像のエッジ方向の変化量がある一定値よりも大き
くなることを収束条件とすることもできる。次式は勾配
によって求めたエッジ方向の変化量の例である。
【０１１４】
【数２４】
【数２５】 図１９にΔθn （ｉ）の例を示す。
【０１１５】収束条件としては、全ての、もしくは一部
の離散点におけるエッジ方向の変化量Δθn （ｉ）
が、ある一定値を越えた時、すなわち、次式の条件が成
立する時とする方法や、
の離散点におけるエッジ方向の変化量Δθn （ｉ）
が、ある一定値を越えた時、すなわち、次式の条件が成
立する時とする方法や、
【数２６】 あるいは、次式のような平均値がある一定値を越えた時
とする方法が考えられる。
とする方法が考えられる。
【０１１６】
【数２７】 但し、Ｎは、離散点の数を示す。
【０１１７】
【発明の効果】第１及び第２の発明であると、磁気標識
付ＭＲ画像中の磁気標識と画像中の対象物の輪郭との交
点、または、前記画像中の複数の磁気標識同士の交点の
検出が容易に行えるようになるため、その検出における
時間と手間が短縮できると共に、画像中の対象物各部の
動きの定量化が正確かつ容易になる。
付ＭＲ画像中の磁気標識と画像中の対象物の輪郭との交
点、または、前記画像中の複数の磁気標識同士の交点の
検出が容易に行えるようになるため、その検出における
時間と手間が短縮できると共に、画像中の対象物各部の
動きの定量化が正確かつ容易になる。
【図１】本発明の第１の実施例の構成図である。
【図２】磁気標識付ＭＲ画像の例（模式図）である。
【図３】初期値の例である。
【図４】離散点移動部の一構成例である。
【図５】複数離散点の番号を表す図（２次元）である。
【図６】収束判定部の一構成例である。
【図７】複数離散点の番号を表す図（３次元）である。
【図８】複数離散点の配置を表す図（３次元）である。
【図９】本発明の第２の実施例の構成図である。
【図１０】輪郭抽出部の一構成例である。
【図１１】輪郭抽出結果の例である。
【図１２】磁気標識対応点の例である。
【図１３】磁気標識検出部の一構成例である。
【図１４】本発明の第３の実施例の構成図である。
【図１５】磁気標識付ＭＲ画像の例である。
【図１６】複数の磁気標識の交点を検出する場合の処理
の流れ図である。
の流れ図である。
【図１７】初期値の例である。
【図１８】離散点の動きの例である。
【図１９】エッジ方向の変化量の例である。
１ 画像データ入力部 ２ 初期値設定部 ３ 離散点移動部 ４ 収束判定部 ５ 出力部 ６ 左心室の心壁 ７ 磁気標識 ８ 離散点（○で表された部分） ９ 画像エネルギー算出部 １０ 外部エネルギー算出部 １１ 次座標算出部 １２ パラメータ記憶部 １３ エネルギー算出部 １４ 記憶部 １５ エネルギー比較部 １６ 輪郭抽出部 １７ 離散点配置部 １８ 磁気標識対応点決定部 １９ 磁気標識検出部 ２０ 磁気標識対応点（×で表された部分） ２１ 動き量測定部 ２２ 記憶部 ２３ 終了判定部 ２４ 磁気標識 ２５ 複数の磁気標識の交点 ２６ 終了判定部 ２７ 初期離散点
フロントページの続き (56)参考文献 特開 平５−204995（ＪＰ，Ａ) 特開 平５−115458（ＪＰ，Ａ) 特表 平６−500035（ＪＰ，Ａ) (58)調査した分野(Int.Cl.7，ＤＢ名) A61B 5/055
Claims (7)
- 【請求項１】磁気標識付ＭＲ画像中の磁気標識と対象物
の輪郭との交点、または、前記画像中の複数の磁気標識
同士の交点の検出方法であって、 磁気標識付ＭＲ画像を入力し、 前記画像中に複数個の離散点を配し、 前記複数個の離散点が有する弾性エネルギーを求め、 前記複数個の離散点が有する磁気標識部分の特徴と磁気
標識以外の部分の特徴によって決定される画像エネルギ
ーを求め、 少なくとも前記弾性エネルギーと前記画像エネルギーと
の総和を最小にするように前記複数個の離散点を収束す
るまで移動させ、 収束した離散点の位置により前記交点の位置を検出する
ことを特徴とする画像処理方法。 - 【請求項２】磁気標識付ＭＲ画像中の磁気標識と対象物
の輪郭との交点、または、前記画像中の複数の磁気標識
同士の交点の検出方法であって、 磁気標識付ＭＲ画像を入力し、 前記画像中に複数個の離散点を配し、 前記複数個の離散点が有する弾性エネルギーを求め、 前記複数個の離散点の内、一部の離散点が有する磁気標
識部分の特徴と磁気標識以外の部分の特徴によって決定
されるエネルギー及びその他の離散点が有する磁気標識
以外の部分の特徴によって決定されるエネルギーとより
なる画像エネルギーを求め、 少なくとも前記弾性エネルギーと前記画像エネルギーと
の総和を最小にするように前記複数個の離散点を収束す
るまで移動させ、 収束した離散点の位置により前記交点の位置を検出する
ことを特徴とする画像処理方法。 - 【請求項３】前記ＭＲ画像、または、前記ＭＲ画像と類
似の画像から、画像中の対象物の輪郭を抽出し、 抽出した輪郭上に複数個の離散点を配し、 前記複数個の離散点の画像中における位置を初期値と
し、 前記複数個の離散点の一部もしくは全部を用いて前記交
点の検出を行うことを特徴とする請求項１もしくは請求
項２記載の画像処理方法。 - 【請求項４】異なる時相毎の画像中の磁気標識と画像中
の対象物の輪郭との交点、または、前記画像中の複数の
磁気標識同士の交点を検出し、 前記異なる時相の画像における交点の検出結果により、
前記交点の動きベクトルを検出することを特徴とする請
求項１もしくは請求項２記載の画像処理方法。 - 【請求項５】前記離散点が有する磁気標識部分によって
決定されるエネルギー及び磁気標識以外の部分によって
決定されるエネルギーよりなる画像エネルギーとして、
画像濃度を表すエネルギーと、画像のエッジ強度を表す
エネルギーの重みの和を用いることを特徴とする請求項
１もしくは請求項２記載の画像処理方法。 - 【請求項６】磁気標識付ＭＲ画像中の磁気標識と対象物
の輪郭との交点、または、前記画像中の複数の磁気標識
同士の交点を検出する装置であって、 磁気標識付ＭＲ画像を入力する入力手段と、 前記画像中に複数個の離散点を配する離散点配置手段
と、 前記複数個の離散点が有する弾性エネルギーを求める弾
性エネルギー算出手段と、 前記複数個の離散点が有する磁気標識部分の特徴と磁気
標識以外の部分の特徴によって決定される画像エネルギ
ーを求める画像エネルギー算出手段と、 少なくとも前記弾性エネルギーと前記画像エネルギーと
の総和を最小にするように前記複数個の離散点を収束す
るまで移動させると、 この離散点移動収束手段により収束した離散点の位置に
より前記交点の位置を検出する交点検出手段とよりなる
ことを特徴とする画像処理装置。 - 【請求項７】磁気標識付ＭＲ画像中の磁気標識と対象物
の輪郭との交点、または、前記画像中の複数の磁気標識
同士の交点を検出する装置であって、 磁気標識付ＭＲ画像を入力する入力手段と、 前記画像中に複数個の離散点を配する離散点配置手段
と、 前記複数個の離散点が有する弾性エネルギーを求める弾
性エネルギー算出手段と、 前記複数個の離散点の内、一部の離散点が有する磁気標
識部分の特徴と磁気標識以外の部分の特徴によって決定
されるエネルギー及びその他の離散点が有する磁気標識
以外の部分の特徴によって決定されるエネルギーとより
なる画像エネルギーを求める画像エネルギー算出手段
と、 少なくとも前記弾性エネルギーと前記画像エネルギーと
の総和を最小にするように前記複数個の離散点を収束す
るまで移動させる離散点移動収束手段と、 この離散点移動収束手段により収束した離散点の位置に
より前記交点の位置を検出する交点検出手段とよりなる
ことを特徴とする画像処理装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP33286193A JP3258158B2 (ja) | 1993-12-27 | 1993-12-27 | 画像処理方法及びその装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP33286193A JP3258158B2 (ja) | 1993-12-27 | 1993-12-27 | 画像処理方法及びその装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH07184877A JPH07184877A (ja) | 1995-07-25 |
| JP3258158B2 true JP3258158B2 (ja) | 2002-02-18 |
Family
ID=18259628
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP33286193A Expired - Fee Related JP3258158B2 (ja) | 1993-12-27 | 1993-12-27 | 画像処理方法及びその装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3258158B2 (ja) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4060615B2 (ja) | 2002-03-05 | 2008-03-12 | 株式会社東芝 | 画像処理装置及び超音波診断装置 |
-
1993
- 1993-12-27 JP JP33286193A patent/JP3258158B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JPH07184877A (ja) | 1995-07-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Saha et al. | A survey on skeletonization algorithms and their applications | |
| US6842638B1 (en) | Angiography method and apparatus | |
| Liu et al. | Shortest path refinement for motion estimation from tagged MR images | |
| JP2018530401A5 (ja) | ||
| CA2486493A1 (en) | Segmentation of 3d medical structures using robust ray propagation | |
| JP2002306483A (ja) | 医用画像診断装置及びその方法 | |
| JP2004105737A (ja) | 心臓磁気共振潅流データのための統合的画像記録方法 | |
| JP2002329216A (ja) | 対象物に含まれる三次元画像をセグメント化する方法 | |
| Makowski et al. | Two-phase active contour method for semiautomatic segmentation of the heart and blood vessels from MRI images for 3D visualization | |
| CN107516314B (zh) | 医学图像超体素分割方法和装置 | |
| Meyer et al. | Tracking myocardial deformation using spatially-constrained velocities | |
| Shi et al. | A model-based integrated approach to track myocardial deformation using displacement and velocity constraints | |
| CN116703994B (zh) | 医学图像配准的方法、计算设备和计算机可读存储介质 | |
| US20070258643A1 (en) | Method, a system, a computer program product and a user interface for segmenting image sets | |
| JP3258158B2 (ja) | 画像処理方法及びその装置 | |
| Uchiyama et al. | Automated classification of cerebral arteries in MRA images and its application to maximum intensity projection | |
| US6934407B2 (en) | Method and system for tag detection and tracking in MRI tagged images | |
| CN109785340A (zh) | 基于多图谱的心脏磁共振图像右心室分割计算设备及方法 | |
| Fan et al. | 3D warping and registration from lung images | |
| JP2002282237A (ja) | 時間経過につれて変形可能な臓器の画像系列のうちの画像に関して定義される単位長あたりの構造の後続する変形 | |
| Macan et al. | Hybrid optical flow and segmentation technique for LV motion detection | |
| JP3639283B2 (ja) | 画像処理装置及び画像処理方法 | |
| EP4270311A1 (en) | Image identification program, image identification method, image identification device | |
| JP3335814B2 (ja) | 画像処理方法及び装置 | |
| Behbahani et al. | Analysing optical flow based methods |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20071207 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20081207 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091207 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091207 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101207 Year of fee payment: 9 |
|
| LAPS | Cancellation because of no payment of annual fees |