JP2009178423A - 画像処理装置及びその方法 - Google Patents

画像処理装置及びその方法 Download PDF

Info

Publication number
JP2009178423A
JP2009178423A JP2008021456A JP2008021456A JP2009178423A JP 2009178423 A JP2009178423 A JP 2009178423A JP 2008021456 A JP2008021456 A JP 2008021456A JP 2008021456 A JP2008021456 A JP 2008021456A JP 2009178423 A JP2009178423 A JP 2009178423A
Authority
JP
Japan
Prior art keywords
teaching
diagnostic
heart
myocardial
image
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.)
Granted
Application number
JP2008021456A
Other languages
English (en)
Other versions
JP5072625B2 (ja
Inventor
Masahide Nishiura
正英 西浦
Satoyuki Takeguchi
智行 武口
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.)
Toshiba Corp
Original Assignee
Toshiba Corp
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 Toshiba Corp filed Critical Toshiba Corp
Priority to JP2008021456A priority Critical patent/JP5072625B2/ja
Priority to EP09250141A priority patent/EP2085930A3/en
Priority to CN2009101267131A priority patent/CN101504727B/zh
Priority to US12/360,277 priority patent/US8170309B2/en
Publication of JP2009178423A publication Critical patent/JP2009178423A/ja
Application granted granted Critical
Publication of JP5072625B2 publication Critical patent/JP5072625B2/ja
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/12Edge-based segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/149Segmentation; Edge detection involving deformable models, e.g. active contour models
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/20Image preprocessing
    • G06V10/26Segmentation of patterns in the image field; Cutting or merging of image elements to establish the pattern region, e.g. clustering-based techniques; Detection of occlusion
    • G06V10/267Segmentation of patterns in the image field; Cutting or merging of image elements to establish the pattern region, e.g. clustering-based techniques; Detection of occlusion by performing operations on regions, e.g. growing, shrinking or watersheds
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/70Arrangements for image or video recognition or understanding using pattern recognition or machine learning
    • G06V10/74Image or video pattern matching; Proximity measures in feature spaces
    • G06V10/75Organisation of the matching processes, e.g. simultaneous or sequential comparisons of image or video features; Coarse-fine approaches, e.g. multi-scale approaches; using context analysis; Selection of dictionaries
    • G06V10/752Contour matching
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20081Training; Learning
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30048Heart; Cardiac
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V2201/00Indexing scheme relating to image or video recognition or understanding
    • G06V2201/03Recognition of patterns in medical or anatomical images
    • G06V2201/031Recognition of patterns in medical or anatomical images of internal organs

Abstract

【課題】解剖学的特徴点が画像上に存在しない場合においても心筋セグメンテーションを可能にする画像処理装置を提供する。
【解決手段】複数の教示用画像について、基準三角形の3個の頂点の位置関係を有するように心筋境界上に指定されている3点の教示用位置に対応した教示用心筋領域境界データを格納し、診断画像中の心筋境界上に、基準三角形の各頂点に対応する3点の診断用位置を指定し、複数の教示用画像のそれぞれについて、3点の教示用位置と3点の診断用位置とを一致させて、診断画像と教示用心筋領域境界データとを比較して類似度を求め、最も高い類似度を有する教示用心筋領域境界データを出力する。
【選択図】 図1

Description

本発明は、医用画像診断装置で得られた心臓画像から、心筋領域をセグメンテーションまたは輪郭抽出する画像処理装置及びその方法に関する。
心臓病の診断を支援する画像処理技術として、心筋の領域を識別するセグメンテーション技術や輪郭抽出技術は重要である。
特許文献1には、心臓断面像における弁輪の位置を解剖学的特徴点として利用して、左心室の輪郭抽出を行う方法が記載されている。
また、特許文献2には、弁輪や乳頭筋といった解剖学的特徴点を用いて輪郭同士の対応付けを行う方法が記載されている。
しかし、例えば左心室の短軸像で心筋のセグメンテーション(または輪郭抽出)を行う場合には、解剖学的特徴点が利用できない場合がある。短軸像の断面位置によっては、解剖学的特徴点となり得る乳頭筋が画像上に含まれない場合があり、このような場合には解剖学的特徴点を利用できない。また、乳頭筋の位置自体も個人差等による変動があり、確実な特徴点とはならない。このような場合においては、解剖学的特徴点を利用した心筋セグメンテーションを行うことはできなかった。
特開2000−217818号公報 特開平10−99334号公報
上記したように、従来技術では解剖学的特徴点を利用した心筋セグメンテーションを行えない場合があるという問題点があった。
そこで本発明は、上記問題点を解決するためになされたものであって、解剖学的特徴点が画像上に存在しない場合においても心筋セグメンテーションを可能にする画像処理装置及びその方法を提供することを目的とする。
本発明は、教示用心臓が写った複数の教示用画像のそれぞれについて、予め定めた基準三角形の3個の頂点の位置関係を有するように前記心臓の心筋境界上に指定されている3点の教示用位置に対応した前記教示用心臓の教示用心筋領域境界データを、学習データとして格納した記憶部と、診断用心臓が写った診断画像を入力する入力部と、前記診断画像中の前記診断用心臓の心筋境界上に、前記基準三角形の相似形と類似するように、前記基準三角形の各頂点に対応する3点の診断用位置を指定する位置指定部と、前記複数の教示用画像のそれぞれについて、前記3点の教示用位置と前記3点の診断用位置とを一致させて、前記診断画像と前記教示用心筋領域境界データとを比較して類似度を求める計算部と、最も高い前記類似度を有する前記教示用心筋領域境界データを、前記診断用心臓の心筋領域境界データとして出力する出力部と、を有する画像処理装置である。
また、本発明は、教示用心臓が写った複数の教示用画像のそれぞれについて、予め定めた基準直線の両端の点の位置関係を有するように前記心臓の心筋境界上に指定されている2個の教示用位置に対応した前記教示用心臓の教示用心筋領域境界データを、学習データとして格納した記憶部と、診断用心臓が写った診断画像を入力する入力部と、前記診断画像中の前記診断用心臓の心筋境界上に、前記基準直線に対応する2点の診断用位置を指定する位置指定部と、前記複数の教示用画像のそれぞれについて、前記2点の教示用位置と前記2点の診断用位置とを一致させて、前記診断画像と前記教示用心筋領域境界データとを比較して類似度を求める計算部と、最も高い前記類似度を有する前記教示用心筋領域境界データを、前記診断用心臓の心筋領域境界データとして出力する出力部と、を有する画像処理装置である。
本発明によれば、明確な解剖学的特徴点がない断面像においても心筋セグメンテーションができる。
(第1の実施形態)
以下、本発明の第1の実施形態の画像処理装置である心筋セグメンテーション装置(以下、本装置という)10について図1〜図5に基づいて説明する。
(1)本装置10の構成
図1は、本実施形態に係わる心解析を行う本装置10を示すブロック図である。
本装置10は、画像が入力する画像入力部22、この画像を記憶するメモリ24、画像を表示する表示部26、画像上の位置を指定する位置指定部12と、指定された位置を用いて心臓の画像上での姿勢を推定処理する姿勢推定部14と、事前に学習した輪郭辞書20を用いて心筋の内側及び外側の境界を抽出するセグメンテーション部16と、セグメンテーションされた境界形状データを出力する出力部18とを備えている。
なお、本装置10は、例えば、汎用のコンピュータ装置を基本ハードウェアとして用いることでも実現することが可能である。すなわち、位置指定部12、姿勢推定部14、セグメンテーション部16及び出力部18は、上記のコンピュータ装置に搭載されたプロセッサにプログラムを実行させることにより実現することができる。このとき、本装置10は、上記のプログラムをコンピュータ装置に予めインストールすることで実現してもよいし、CD−ROMなどの記憶媒体に記憶して、あるいはネットワークを介して上記のプログラムを配布して、このプログラムをコンピュータ装置に適宜インストールすることで実現してもよい。また、メモリは、上記のコンピュータ装置に内蔵あるいは外付けされたメモリ、ハードディスクもしくはCD−R、CD−RW、DVD−RAM、DVD−Rなどの記憶媒体などを適宜利用して実現することができる。
次に、図1及び図2を用いて、本装置10の動作について説明する。図2は、本装置10の動作を示すフローチャートである。
(2)学習処理
まず、心筋セグメンテーションの処理を行う事前の処理として、輪郭辞書20の学習を行う。輪郭辞書20の学習は心筋セグメンテーション処理時には行う必要はない。
以下に輪郭辞書20の学習処理の手順について説明する。
(2−1)基準三角形の設定
まず、心筋境界上に指定する3点の位置関係を表す基準三角形を定める。基準三角形は、抽出対象とする心筋の形状である断面カテゴリーに応じて規定する。
心筋の形状は、観察する断面の位置によって変化する。例えば、心筋形状が円形に近いような短軸断面像では、基準三角形を正三角形とすると、円周上に偏りなく頂点が配置されることとなり、精度向上の効果が得られる。
また、長軸断面像では二等辺三角形あるいは鋭角三角形とすると、心筋の両端をカバーすることができ、精度向上の効果が得られる。
このように、基準三角形を断面カテゴリー毎に切り替えると効率的である。
(2−2)3点の教示用位置の指定
次に、複数枚の教示用画像を準備する。そして、各教示用画像における教示用心筋領域境界上において前記基準三角形の相似形に類似するように3点の教示用位置、すなわち、三角形の頂点の位置を指定する。例として、図3を用いて説明する。図3では、短軸断面像において基準三角形が正三角形の場合について図示している。
この場合、教示用画像における教示用心筋境界上で正三角形に類似するように3点の位置を指定する。この3点の位置は解剖学的特徴点である必要はない。
(2−3)学習データの作成
次に、各教示用画像における、少なくとも3点の教示用位置と教示用心筋領域境界データを用いて輪郭辞書20に格納する学習データを作成する。
まず、教示用画像における3点の教示用位置を、予め定めた正規化三角形の頂点位置に重なるような座標変換式を求める。ここで、「正規化三角形」とは、基準三角形とは異なる概念であり、後から説明する診断用画像と教示用画像との類似度を求めるときに、互いに正規化して比較し易くするために生まれた概念である。この正規化三角形は、仮想の正規化座標系に存在する3点の座標の位置から構成され、その形状は必ずしも基準三角形と相似形である必要はなく、異なる形状でもよい。そして、この座標変換式は、教示用画像の座標系から正規化座標系に変換する式である。
次に、この座標変換式により教示用画像における教示用輪郭データの座標を正規化した上で、輪郭形状の学習を行う。
輪郭形状の学習には、例えば一般的なActive Shape Model(文献1、T.F.Cootes, C.J.Taylor, D.H.Cooper and J.Graham, 「Active shape models - their training and application」, Computer Vision and Image Understanding, 61(1):38-59, Jan. 1995.)などを用いることができる。
Active Shape Modelでは、教示用輪郭形状データとして、輪郭形状の分布と輪郭周辺の画像の輝度プロファイルとを学習する。
輪郭形状の分布の学習では、例えば、輪郭を表現する頂点の座標データを並べたベクトルの分布を主成分分析し、固有値の大きい上位の形状ベクトルなどを輪郭形状データとして輪郭辞書20に保存する。
画像輝度のプロファイルの学習では、例えば、輪郭に直交する直線上の画像輝度データ(輝度プロファイル)の分布を主成分分析した結果を、輝度プロファイルとして輪郭辞書に保存する。
(2−4)まとめ
以上の学習処理によって、複数枚の教示用画像のそれぞれについて、位置指定部12で定める3点を基準とした輪郭形状の学習が行われ学習データが生成される。つまり、心筋の画像中での姿勢の分布までをも含めた輪郭学習を行うことができる。
生成された学習データを保存しておき、後述する心筋セグメンテーション処理に使用する。
(3)心筋セグメンテーション処理
次に、心筋セグメンテーション処理の手順について説明する。
(3−1)診断用画像の入力
まず、診断用心臓が写った診断用画像が画像入力部22によりメモリ24に格納される。診断用画像は超音波診断装置、X線CT装置、MRI装置等の医用画像診断装置から入力しても良いし、画像サーバー等に保存された画像を入力としても良い。
(3−2)3点の指定
メモリ24に格納された診断用画像は表示部26に表示され、ユーザは表示された診断用画像を参照しながら、位置指定部12を用いて、診断用心臓の心筋境界上の3点の診断用位置を指定する。指定する3点は、心筋境界上で、事前に規定された基準三角形の相似形に類似するように配置する。
例えば、図4に示すように診断用画像中の心臓短軸像において解剖学的特徴点が明確でない場合には、心筋境界のうちで最も上部になる点と、その点から例えばおおよそ120度の2点を指定する。
また、解剖学的特徴点は明確ではないが、短軸像のおおよその回転姿勢をユーザが判断できる場合には、その回転角に合わせて入力する3点(基準三角形の頂点)を回転させた状態とする(図5参照)。
なお、3点の診断用位置の入力時に入力位置関係のガイドとなるように、表示部26上にアイコンを表示してもよい。例えば、図4中に示すようにY字型または三角形状のアイコンを表示し3点の角度の目安をガイドすることによって、三角形の位置関係での入力が容易となる。
また、3点の診断用位置の位置関係が、診断画像の画像座標系において心筋境界上で最上部、最下部、最右部、最左部の1点と、心筋境界上の他の2点によって基準三角形の相似形に類似する位置関係となるように定めても良い。例えば、超音波画像を経胸壁で得た場合には、心臓の前壁が画像上部に撮像されるなど、心臓の大まかな姿勢が一定の条件下にあることも多く、このような場合に上記のような位置関係を定めることでより安定した心筋セグメンテーションが可能となる。
(3−3)姿勢推定
次に、姿勢推定部14では、診断用画像の指定された3点の診断用位置が、学習時に使用した正規化座標系の正規化三角形の頂点位置に重なるような座標変換式を求める。
この座標変換式は、例えば3点の診断用位置の座標と正規化三角形の頂点座標の関係からアフィン変換行列を計算することによって求められる。
また、3点の診断用位置の座標と正規化三角形の頂点座標の関係から、位置、スケール、回転のパラメータを計算することによって求められる。
座標変換式を計算することによって、診断用画像上での心筋の位置、スケール、回転角といった姿勢が求められたことになる。
(3−4)境界輪郭の抽出
セグメンテーション部16では、位置指定部12で指定された診断用画像の3点の診断位置と、診断用画像と、学習データから心筋の境界輪郭を抽出する。以下では、Active Shape Modelを利用した例について説明する。
まず、輪郭辞書20の中の一つの学習データの教示用輪郭形状データを、正規化三角形の各頂点の位置と姿勢推定部14で求めた座標変換式の逆変換により、診断用画像上の座標へ変換する。
次に、診断用画像中の輪郭周辺の輝度プロファイルと、教示用輪郭形状データの輝度プロファイルとの類似度を求める。
次に、この処理を輪郭辞書20に格納されている複数枚の教示用画像に関する輪郭形状データの全てについて繰り返し行う、
次に、輪郭辞書20の輪郭形状の分布範囲において最も輝度プロファイルの類似度が、最も高い輪郭形状データの探索を行う。
次に、この最も類似度が高い輪郭形状データが、診断用画像の診断用心臓の心筋の境界輪郭を表すものとして抽出する。
(3−5)出力
最後に、出力部18は、抽出された輪郭形状データを出力する。
(4)効果
以上のように本実施形態であると、3点の教示用位置を用いた輪郭形状の学習と、3点位置指定による輪郭の抽出を用いれば、診断用画像上で解剖学的特徴点が明確でない場合においても、ユーザ入力に応じた適切な心筋セグメンテーション結果を得ることができる。これにより、関心領域の指定などのユーザ操作を簡単化する点で大きな効果を得ることができる。
(第2の実施形態)
第2の実施形態の本装置10について図6に基づいて説明する。
本実施形態は、心筋輪郭形状の学習時に、心筋の部分領域の境界位置を合わせて学習しておくことで、3点位置指定により心筋セグメンテーションと部分領域境界の位置指定を同時に行うことができる(図6参照)。
「心筋の部分領域」とは、解剖学的に使用される領域であって、例えば、心基部前壁、心基部後壁などと呼ばれる領域である。
(1)学習処理
まず、学習処理について説明する。
心筋輪郭形状の学習処理時には、教示用輪郭形状データと、教示用画像と、教示用部分領域境界位置データを用いる。
第1の実施形態と同様に、複数枚の教示用画像のそれぞれについて、規定の基準三角形の相似形に類似するように3点の教示用位置を指定する。
次に、3点の教示用位置が正規化三角形の頂点位置に重なるような座標変換式により教示用輪郭データ及び教示用部分領域境界位置データの座標を変換し、輪郭形状と部分境界位置を合わせて学習する。
具体的には、輪郭形状データも部分境界位置データも頂点の座標データであるので、これら座標データを羅列して、1つのベクトルとした上で学習に用いる。学習により生成される学習データには、輪郭形状データ(輪郭形状の分布と輝度プロファイル)と部分領域境界位置データが含まれることになる。
(2)心筋セグメンテーション処理
次に、心筋セグメンテーション処理について説明する。
心筋セグメンテーション部16では、位置指定部12で指定された3点位置と診断用画像と輪郭辞書20から心筋の境界輪郭及び心筋の部分領域境界位置を抽出する。
まず、第1の実施形態と同様に、指定された3点の診断用位置が学習時に使用した正規化三角形の頂点位置に重なるような座標変換式を求める。
次に、学習データ中の輪郭形状データを座標逆変換式により診断用画像上の座標へ変換する。
次に、診断用画像における輪郭周辺の輝度プロファイルと、輪郭辞書20中の輝度プロファイルとの類似度を評価する。
次に、この処理を複数枚の教示用画像のついてぞれぞれ行う。
次に、輪郭辞書20中の輪郭形状の分布範囲において最も輝度プロファイルの類似度が高い輪郭形状の探索を行うことで心筋の境界輪郭を抽出する。
次に、抽出された輪郭形状と組みで保存されている部分領域境界位置データを読み出すことで、部分領域境界の位置を得ることができる。
(3)効果
以上により第2の実施形態であると、心筋の部分領域の境界位置を合わせて学習しておくことで、同じ3点位置指定によって心筋セグメンテーションと部分領域境界の位置指定をも同時に行うことができる。このことは、関心領域の指定などのユーザ操作を簡単化する点でさらに大きな効果を有する。
(第3の実施形態)
第3の実施形態の本装置10について図8に基づいて説明する。
本実施形態では、基準三角形に代えて基準直線を用いて、輪郭上の2点の位置指定により心筋セグメンテーションを行う。
本装置10について、上記実施形態と異なる点を中心に説明する。
(1)学習処理
事前に行う学習処理では、心筋境界上に指定する2点を通る直線の傾きを表す基準直線を定める。例えば、輪郭上の最上点と最下点を指定する場合には基準直線は縦直線と規定する。
まず、複数枚の教示用画像のそれぞれについて、基準直線を教示用心筋領域の想定重心を通るように配置した場合に心筋領域境界と交差する2点(すなわち、基準直線の両端の点)を2点の教示用位置とする(図8参照)。なお、2点の教示用位置は解剖学的特徴点となる必要はない。
次に、少なくとも前記教示用2点位置と教示用心筋領域境界データを用いて学習データを作成する。2点の教示用位置の座標を用いて、この2点間距離が、正規化直線上の単位長となるような座標変換式を求め正規化処理を行い、輪郭形状を学習し、輪郭辞書20データを作成する。
ここで、「正規化直線」とは、基準三角形と同じ概念であって、基準直線とは異なる概念であり、後から説明する診断用画像と教示用画像との類似度を求めるときに、互いに正規化して比較し易くするために生まれた概念である。この正規化直線は、仮想の正規化座標系に存在する2点の座標の位置から構成される。
(2)心筋セグメンテーション処理
まず、心筋セグメンテーション処理では、心筋領域の想定重心に規定の基準直線を配置した場合に心筋境界と交差する2点の診断用位置を指定する。
次に、この2点間距離が正規化直線上の単位長となるような座標変換式を求める。
次に、学習データ中の輪郭形状データを座標変換式の逆変換により、診断用画像上の座標へ変換する。
次に、診断用画像中の輪郭周辺の輝度プロファイルと、輪郭辞書20中の輝度プロファイルとの類似度を求める。
次に、輪郭辞書20中の輪郭形状の分布範囲において最も輝度プロファイルの類似度が最も高い、教示用画像の輪郭形状の探索を行うことで心筋の境界輪郭を抽出する。
(3)効果
以上により本実施形態によれば、基準図形として直線を用い、2点の入力により心筋セグメンテーションを行うことで、学習時の形状変動幅が大きくなり精度面で不利になる反面、指定する点数を減らすことができる。また、断面像内に解剖学的特徴点のない場合においても、より少ない2点の入力により心筋セグメンテーションを行うことができる。
(変更例)
なお、本発明は上記各実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。また、上記実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。
例えば、上記各実施形態の説明においては心臓短軸断面像に対する使用例を示したが、短軸断面像に限らない。例えば、長軸断面像に対しては図7に示すように3点を指定することで、同様の効果を得ることができる。
この場合に、基準三角形または基準直線を、心臓の断面カテゴリー(長軸断面像、短軸断面像)毎に切り替える切り替え部を有するようにしてもよい。
本発明の第1の実施形態に係わる心筋セグメンテーション装置の構成を示すブロック図である。 同じく動作を示すフローチャートである。 辞書生成時の心臓短軸像における3点指定を説明する図である。 心臓短軸像における3点指定を説明する図である。 心臓短軸像においてユーザが意図的に心臓の姿勢を入力する場合の3点指定を説明する図である。 第2の実施形態における、心筋の輪郭と部分領域の境界位置の同時抽出を説明する図である。 変更例における、心臓長軸断面像における3点指定の例を説明する図である。 第3の実施形態における、心臓短軸像における2点指定を説明する図である。
符号の説明
10 心筋セグメンテーション装置
12 位置指定部
14 姿勢推定部
16 セグメンテーション部
18 出力部
20 輪郭辞書
22 画像入力部
24 メモリ
26 表示部

Claims (13)

  1. 教示用心臓が写った複数の教示用画像のそれぞれについて、予め定めた基準三角形の3個の頂点の位置関係を有するように前記心臓の心筋境界上に指定されている3点の教示用位置に対応した前記教示用心臓の教示用心筋領域境界データを、学習データとして格納した記憶部と、
    診断用心臓が写った診断画像を入力する入力部と、
    前記診断画像中の前記診断用心臓の心筋境界上に、前記基準三角形の相似形と類似するように、前記基準三角形の各頂点に対応する3点の診断用位置を指定する位置指定部と、
    前記複数の教示用画像のそれぞれについて、前記3点の教示用位置と前記3点の診断用位置とを一致させて、前記診断画像と前記教示用心筋領域境界データとを比較して類似度を求める計算部と、
    最も高い前記類似度を有する前記教示用心筋領域境界データを、前記診断用心臓の心筋領域境界データとして出力する出力部と、
    を有する画像処理装置。
  2. 教示用心臓が写った複数の教示用画像のそれぞれについて、予め定めた基準直線の両端の点の位置関係を有するように前記心臓の心筋境界上に指定されている2個の教示用位置に対応した前記教示用心臓の教示用心筋領域境界データを、学習データとして格納した記憶部と、
    診断用心臓が写った診断画像を入力する入力部と、
    前記診断画像中の前記診断用心臓の心筋境界上に、前記基準直線に対応する2点の診断用位置を指定する位置指定部と、
    前記複数の教示用画像のそれぞれについて、前記2点の教示用位置と前記2点の診断用位置とを一致させて、前記診断画像と前記教示用心筋領域境界データとを比較して類似度を求める計算部と、
    最も高い前記類似度を有する前記教示用心筋領域境界データを、前記診断用心臓の心筋領域境界データとして出力する出力部と、
    を有する画像処理装置。
  3. 前記教示用画像について、前記基準三角形の相似形と類似するように、前記教示用心臓の心筋境界上に、3点の前記教示用位置を指定する教示用位置指定部と、
    3点の前記教示用位置に対応した前記教示用心臓の教示用心筋領域境界データを前記教示用画像から抽出して前記記憶部に格納する学習部と、
    を有する請求項1記載の画像処理装置。
  4. 前記教示用画像について、前記基準直線に対応するように、前記教示用心臓の心筋境界上に、2点の前記教示用位置を指定する教示用位置指定部と、
    2点の前記教示用位置に対応した前記教示用心臓の教示用心筋領域境界データを前記教示用画像から抽出して前記記憶部に格納する学習部と、
    を有する請求項2記載の画像処理装置。
  5. 前記記憶部が格納している前記各学習データが、前記教示用心臓の心筋輪郭形状データ及び心筋部分領域の境界データを含み、
    前記出力部は、最も高い前記類似度を有する前記教示用心筋領域境界データに対応した前記心筋輪郭形状データ及び前記心筋部分領域の境界データに基づいて、前記診断用心臓の心筋の境界形状と心筋部分領域の境界位置を出力する、
    請求項1または2記載の画像処理装置。
  6. 前記診断画像が心臓短軸断面像である場合に、
    前記位置指定部における前記3点の診断用位置が、前記診断画像における画像座標系において心筋境界上の最上部、最下部、最右部、または、最左部の中の1点と、前記心筋境界上の他の2点である、
    請求項1記載の画像処理装置。
  7. 前記心筋領域境界データが、心筋の内側及び外側の境界輪郭形状を示している、
    請求項1または2記載の画像処理装置。
  8. 前記基準三角形を、心臓の断面カテゴリー毎に切り替える切り替え部を有する、
    請求項1記載の画像処理装置。
  9. 前記基準直線を、心臓の断面カテゴリー毎に切り替える切り替え部を有する、
    請求項2記載の画像処理装置。
  10. 教示用心臓が写った複数の教示用画像のそれぞれについて、予め定めた基準三角形の3個の頂点の位置関係を有するように前記心臓の心筋境界上に指定されている3点の教示用位置に対応した前記教示用心臓の教示用心筋領域境界データを、学習データとして格納する記憶ステップと、
    診断用心臓が写った診断画像を入力する入力ステップと、
    前記診断画像中の前記診断用心臓の心筋境界上に、前記基準三角形の相似形と類似するように、前記基準三角形の各頂点に対応する3点の診断用位置を指定する位置指定ステップと、
    前記複数の教示用画像のそれぞれについて、前記3点の教示用位置と前記3点の診断用位置とを一致させて、前記診断画像と前記教示用心筋領域境界データとを比較して類似度を求める計算ステップと、
    最も高い前記類似度を有する前記教示用心筋領域境界データを、前記診断用心臓の心筋領域境界データとして出力する出力ステップと、
    を有する画像処理方法。
  11. 教示用心臓が写った複数の教示用画像のそれぞれについて、予め定めた基準直線の両端の点の位置関係を有するように前記心臓の心筋境界上に指定されている2個の教示用位置に対応した前記教示用心臓の教示用心筋領域境界データを、学習データとして格納する記憶ステップと、
    診断用心臓が写った診断画像を入力する入力ステップと、
    前記診断画像中の前記診断用心臓の心筋境界上に、前記基準直線に対応する2点の診断用位置を指定する位置指定ステップと、
    前記複数の教示用画像のそれぞれについて、前記2点の教示用位置と前記2点の診断用位置とを一致させて、前記診断画像と前記教示用心筋領域境界データとを比較して類似度を求める計算ステップと、
    最も高い前記類似度を有する前記教示用心筋領域境界データを、前記診断用心臓の心筋領域境界データとして出力する出力ステップと、
    を有する画像処理方法。
  12. 教示用心臓が写った複数の教示用画像のそれぞれについて、予め定めた基準三角形の3個の頂点の位置関係を有するように前記心臓の心筋境界上に指定されている3点の教示用位置に対応した前記教示用心臓の教示用心筋領域境界データを、学習データとして格納する記憶機能と、
    診断用心臓が写った診断画像を入力する入力機能と、
    前記診断画像中の前記診断用心臓の心筋境界上に、前記基準三角形の相似形と類似するように、前記基準三角形の各頂点に対応する3点の診断用位置を指定する位置指定機能と、
    前記複数の教示用画像のそれぞれについて、前記3点の教示用位置と前記3点の診断用位置とを一致させて、前記診断画像と前記教示用心筋領域境界データとを比較して類似度を求める計算機能と、
    最も高い前記類似度を有する前記教示用心筋領域境界データを、前記診断用心臓の心筋領域境界データとして出力する出力機能と、
    をコンピュータによって実現する画像処理プログラム。
  13. 教示用心臓が写った複数の教示用画像のそれぞれについて、予め定めた基準直線の両端の点の位置関係を有するように前記心臓の心筋境界上に指定されている2個の教示用位置に対応した前記教示用心臓の教示用心筋領域境界データを、学習データとして格納する記憶機能と、
    診断用心臓が写った診断画像を入力する入力機能と、
    前記診断画像中の前記診断用心臓の心筋境界上に、前記基準直線に対応する2点の診断用位置を指定する位置指定機能と、
    前記複数の教示用画像のそれぞれについて、前記2点の教示用位置と前記2点の診断用位置とを一致させて、前記診断画像と前記教示用心筋領域境界データとを比較して類似度を求める計算機能と、
    最も高い前記類似度を有する前記教示用心筋領域境界データを、前記診断用心臓の心筋領域境界データとして出力する出力機能と、
    をコンピュータによって実現する画像処理プログラム。
JP2008021456A 2008-01-31 2008-01-31 画像処理装置及びその方法 Expired - Fee Related JP5072625B2 (ja)

Priority Applications (4)

Application Number Priority Date Filing Date Title
JP2008021456A JP5072625B2 (ja) 2008-01-31 2008-01-31 画像処理装置及びその方法
EP09250141A EP2085930A3 (en) 2008-01-31 2009-01-20 Image processing apparatus and method thereof
CN2009101267131A CN101504727B (zh) 2008-01-31 2009-01-23 图像处理设备及其方法
US12/360,277 US8170309B2 (en) 2008-01-31 2009-01-27 Image processing apparatus and method thereof

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2008021456A JP5072625B2 (ja) 2008-01-31 2008-01-31 画像処理装置及びその方法

Publications (2)

Publication Number Publication Date
JP2009178423A true JP2009178423A (ja) 2009-08-13
JP5072625B2 JP5072625B2 (ja) 2012-11-14

Family

ID=40640263

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2008021456A Expired - Fee Related JP5072625B2 (ja) 2008-01-31 2008-01-31 画像処理装置及びその方法

Country Status (4)

Country Link
US (1) US8170309B2 (ja)
EP (1) EP2085930A3 (ja)
JP (1) JP5072625B2 (ja)
CN (1) CN101504727B (ja)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20160012596A1 (en) * 2013-03-21 2016-01-14 Koninklijke Philips N.V. View classification-based model initialization
JP2017148438A (ja) * 2016-02-26 2017-08-31 東芝メディカルシステムズ株式会社 医用画像処理装置、超音波診断装置、及び医用画像処理プログラム
JP2021085757A (ja) * 2019-11-27 2021-06-03 国立大学法人神戸大学 顕微鏡による焦点画像群を用いた形状計測方法及び装置

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6465569B2 (ja) * 2014-06-11 2019-02-06 キヤノン株式会社 画像処理方法、および画像処理装置
US10682181B2 (en) * 2017-09-06 2020-06-16 Biosense Webster (Israel) Ltd. Methods and systems for modeling and registration of 3-dimensional images of the heart
EP3570249A1 (en) 2018-05-14 2019-11-20 Koninklijke Philips N.V. Preprocessing medical image data for machine learning
CN109947394B (zh) * 2019-03-11 2023-03-10 湖北工业大学 一种基于三角形覆盖的图像二维分形维数计算方法
CN117243637B (zh) * 2023-10-19 2024-04-19 河北港口集团有限公司秦皇岛中西医结合医院 一种超声心动图图像识别方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050238233A1 (en) * 2002-08-20 2005-10-27 Miguel Mulet Parada Computation of contour
US20060171586A1 (en) * 2004-11-08 2006-08-03 Bogdan Georgescu Method of database-guided segmentation of anatomical structures having complex appearances

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5734739A (en) * 1994-05-31 1998-03-31 University Of Washington Method for determining the contour of an in vivo organ using multiple image frames of the organ
JP3406785B2 (ja) 1996-09-26 2003-05-12 株式会社東芝 心機能解析支援装置
US6346124B1 (en) * 1998-08-25 2002-02-12 University Of Florida Autonomous boundary detection system for echocardiographic images
JP3668629B2 (ja) 1999-01-29 2005-07-06 株式会社東芝 画像診断装置及び画像処理方法
US6917826B2 (en) * 2002-02-01 2005-07-12 Siemens Corporate Research, Inc. Automatic and quantitative assessment of image quality in SPECT myocardial perfusion imaging
US7352888B2 (en) * 2004-07-30 2008-04-01 Carestream Health, Inc. Method for computer recognition of projection views and orientation of chest radiographs

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050238233A1 (en) * 2002-08-20 2005-10-27 Miguel Mulet Parada Computation of contour
US20060171586A1 (en) * 2004-11-08 2006-08-03 Bogdan Georgescu Method of database-guided segmentation of anatomical structures having complex appearances

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20160012596A1 (en) * 2013-03-21 2016-01-14 Koninklijke Philips N.V. View classification-based model initialization
JP2016519363A (ja) * 2013-03-21 2016-06-30 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. ビュー分類に基づいたモデル初期化
US10109072B2 (en) * 2013-03-21 2018-10-23 Koninklijke Philips N.V. View classification-based model initialization
JP2017148438A (ja) * 2016-02-26 2017-08-31 東芝メディカルシステムズ株式会社 医用画像処理装置、超音波診断装置、及び医用画像処理プログラム
JP2021085757A (ja) * 2019-11-27 2021-06-03 国立大学法人神戸大学 顕微鏡による焦点画像群を用いた形状計測方法及び装置

Also Published As

Publication number Publication date
US20090196474A1 (en) 2009-08-06
US8170309B2 (en) 2012-05-01
EP2085930A2 (en) 2009-08-05
CN101504727A (zh) 2009-08-12
CN101504727B (zh) 2012-11-07
EP2085930A3 (en) 2011-10-26
JP5072625B2 (ja) 2012-11-14

Similar Documents

Publication Publication Date Title
JP5072625B2 (ja) 画像処理装置及びその方法
Tobon-Gomez et al. Benchmark for algorithms segmenting the left atrium from 3D CT and MRI datasets
JP5610129B2 (ja) 3次元テンプレート変形方法、装置及びプログラム
US9275190B2 (en) Method and system for generating a four-chamber heart model
JP6059261B2 (ja) マルチモーダル画像統合におけるレジストレーション精度を改善するインテリジェントなランドマーク選択
US9155470B2 (en) Method and system for model based fusion on pre-operative computed tomography and intra-operative fluoroscopy using transesophageal echocardiography
CN106462963B (zh) 用于自适应放射治疗中自动勾画轮廓的系统和方法
US9292917B2 (en) Method and system for model-based fusion of computed tomography and non-contrasted C-arm computed tomography
RU2711140C2 (ru) Редактирование медицинских изображений
US8948484B2 (en) Method and system for automatic view planning for cardiac magnetic resonance imaging acquisition
US10891787B2 (en) Apparatus and method for creating biological model
JP5832938B2 (ja) 画像処理装置、方法及びプログラム
EP1851721B1 (en) A method, a system and a computer program for segmenting a surface in a multidimensional dataset
JP2007530088A (ja) オブジェクト特有のセグメンテーション
JP2011120825A (ja) 医用画像表示装置および方法並びにプログラム
WO2017180097A1 (en) Deformable registration of intra and preoperative inputs using generative mixture models and biomechanical deformation
JP6032610B2 (ja) 表示処理プログラム、表示処理方法及び表示処理装置
JP5403431B2 (ja) 断層画像処理方法及び装置
US9367913B2 (en) Choosing anatomical variant model for image segmentation
KR102629331B1 (ko) Cta 영상에서 계층적 변형을 이용한 심혈관 비강체 정합 방법, 이를 수행하기 위한 기록매체 및 장치
EP4270311A1 (en) Image identification program, image identification method, image identification device
Meyer et al. A multi-modality segmentation framework: application to fully automatic heart segmentation
Aase Graph Convolutional Neural Networks for Automated Echocardiography View Recognition: A Holistic Approach
JP2023547525A (ja) 解剖学的構造の3d表現のレンダリング及び表示
Mcleod Reduced-order statistical models of cardiac growth, motion and blood flow: application to the tetralogy of Fallot heart

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20101019

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20120705

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

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20120821

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20150831

Year of fee payment: 3

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