JP2008194239A - 画像処理装置及びその方法 - Google Patents
画像処理装置及びその方法 Download PDFInfo
- Publication number
- JP2008194239A JP2008194239A JP2007032755A JP2007032755A JP2008194239A JP 2008194239 A JP2008194239 A JP 2008194239A JP 2007032755 A JP2007032755 A JP 2007032755A JP 2007032755 A JP2007032755 A JP 2007032755A JP 2008194239 A JP2008194239 A JP 2008194239A
- Authority
- JP
- Japan
- Prior art keywords
- image
- center
- cavity
- filter
- boundary
- 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.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/12—Edge-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/40—Extraction of image or video features
- G06V10/46—Descriptors for shape, contour or point-related descriptors, e.g. scale invariant feature transform [SIFT] or bags of words [BoW]; Salient regional features
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30048—Heart; Cardiac
Landscapes
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Multimedia (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
【課題】臓器にある空洞部の境界推定を可能とする画像処理装置を提供する。
【解決手段】空洞部を有する臓器を撮像した断面画像を入力する画像入力部110と、断面画像に対して2つ以上のスケールを変えた空間フィルタを用いて2つ以上の出力画像を得るフィルタ処理部120と、出力画像を用いて空洞部の中心位置を推定する物体中心推定部130と、出力画像と推定した中心位置とを利用して空洞部の壁面位置の初期境界を推定する初期境界推定部140と、初期境界を初期値として境界位置を決定する境界決定部150とを備える。
【選択図】 図1
【解決手段】空洞部を有する臓器を撮像した断面画像を入力する画像入力部110と、断面画像に対して2つ以上のスケールを変えた空間フィルタを用いて2つ以上の出力画像を得るフィルタ処理部120と、出力画像を用いて空洞部の中心位置を推定する物体中心推定部130と、出力画像と推定した中心位置とを利用して空洞部の壁面位置の初期境界を推定する初期境界推定部140と、初期境界を初期値として境界位置を決定する境界決定部150とを備える。
【選択図】 図1
Description
本発明は、内部に空洞を有する臓器を撮像した画像からその空洞の輪郭を自動的に推定する画像処理装置及びその方法に関する。
従来の画像処理装置には、特許文献1に開示されているように、関心物体の輪郭を得るために、関心物体を含む関心領域(ROI)を指定した上で、画像中の関心領域に対して2値化処理を施して輪郭を得る方法があった。しかし、関心領域は手動で指定しなければならないという問題点があった。
また、特許文献2に開示されているように、心筋の内壁境界に基づいて、動的輪郭モデルによって心筋の外壁境界を推定するものがあった。しかし、内壁境界については何らかの方法で求めなくてはならないという問題点があった。
さらに、特許文献3に開示されているように、診断画像の中心点を求め、この中心点に楕円アークモデルをあてはめ、曲線境界を求める方法があった。しかし、診断画像を直接用いて検出するため、その中心点などが検出し難いという問題点があった。
特開平8‐336503公報
特開平10‐229979公報
特許第3194741号公報
上述したように、従来技術には、心筋の内外壁面境界を推定するために、何らかの初期値が必要であったり、診断装置の作業従事者になんからの操作を求める必要があったり、診断画像を直接用いるため検出し難いという問題点があった。
そこで本発明は、上記問題点を解決するためになされたものであって、内部に空洞を有する臓器を撮像した画像からその空洞の輪郭を自動的に推定する画像処理装置及びその方法を提供することを目的とする。
本発明は、空洞部を有する臓器を撮像した画像を入力する画像入力部と、前記空洞部に対応する閉領域の画素情報とその閉領域の周囲領域の画素情報を比較して、前記閉領域の中心位置の画素情報を強調して出力する空間フィルタを、前記入力画像に施して出力画像を取得するフィルタ処理部と、前記出力画像を用いて前記空洞部の中心位置を推定する物体中心推定部と、前記出力画像と前記推定した中心位置とを用いて、前記空洞部の壁面位置である境界位置を決定する境界決定部と、を備える画像処理装置である。
本発明によれば、手動入力による初期値が必要ない臓器の空洞部の内外壁面推定ができる。
以下、本発明の実施形態に係わる画像処理装置について図面に基づいて説明する。
(第1の実施形態)
本発明の第1の実施形態に係わる画像処理装置について図1から図5に基づいて説明する。本実施形態では、対象臓器が心臓である場合を例にし、空洞部である左心室を関心物体とした左心室心筋内外の境界を推定する場合について述べる。
本発明の第1の実施形態に係わる画像処理装置について図1から図5に基づいて説明する。本実施形態では、対象臓器が心臓である場合を例にし、空洞部である左心室を関心物体とした左心室心筋内外の境界を推定する場合について述べる。
(1)画像処理装置の構成
図1は、本実施形態に係わる画像処理装置を示すブロック図である。
図1は、本実施形態に係わる画像処理装置を示すブロック図である。
画像処理装置は、心臓の断面画像を取得する画像入力部110と、断面画像に対して空間フィルタを施して出力画像を得るフィルタ処理部120と、出力画像から物体中心を推定する物体中心推定部130と、推定された物体中心と出力画像を用いて空洞部の初期境界を推定する初期境界推定部140と、得られた初期境界を初期値として最終的な境界を定める境界決定部150とを備えている。
(2)画像処理装置の動作
次に図1及び図2を用いて、本実施形態に係わる画像処理装置の動作について説明する。なお、図2は、本実施形態に係わる画像処理装置の動作を示すフローチャートである。
次に図1及び図2を用いて、本実施形態に係わる画像処理装置の動作について説明する。なお、図2は、本実施形態に係わる画像処理装置の動作を示すフローチャートである。
(2−1)画像入力部110
画像入力部110では、空洞部を含む断面画像を取得する(図2のステップA1参照)。
画像入力部110では、空洞部を含む断面画像を取得する(図2のステップA1参照)。
例えばここでは、超音波診断装置を用いて心臓の2次元断面像を取得する。プローブの位置や角度により得られる断面像が異なるが、ここでは、非特許文献4(日本医師会編、"心エコーのABC"、pp.6-7、中山書店、1995)にあるような、被験者を半左側臥位に寝かせ、胸骨左縁第3、4肋間からアプローチした際に得られる乳頭筋レベルでの胸骨左縁短軸断面像を例として説明する。
乳頭筋レベルでの左心室断面像の模式図を図3に示す。断面画像内には空洞部である左心室510以外に、右心室520、心筋530が撮像されている。
(2−2)フィルタ処理部120
次に、フィルタ処理部120では、断面画像に対して空間フィルタを施す。
次に、フィルタ処理部120では、断面画像に対して空間フィルタを施す。
図3に示すように、短軸像では左心室領域内部が比較的輝度が低くほぼ円形状として撮像され、心筋部分は比較的輝度が高くなる。そこでフィルタ処理部120では、2つの領域の輝度を比較することができる空間フィルタとして、非特許文献1(Tony Lindeberg, "Feature Detection with Automatic Scale Selection", International Journal of Computer Vision, vol.30, no.2, pp.79-116, 1998 )に記載のラプラシアンガウシアンフィルタ(Laplacian Of Gaussian : LOG)を用いる。ラプラシアンガウシアンフィルタの式を式(1)に示す。
ここでI(x,y)は入力断面画像、G(σ)は2次元ガウシアンフィルタ、Lは2次元ラプラシアンフィルタ、*は畳み込み積分を示す記号、F(x,y)は出力画像であり、σはガウシアンフィルタのボケ量を表すパラメータである。
2次元のラプラシアンガウシアンフィルタは図4に示すような形状のフィルタであり、対象画素を中心とする領域と、さらにその周辺部の2つの領域の重みつき輝度値の差で計算される。
パラメータσは空間フィルタのスケールを調整するスケールパラメータであり、σによって比較する領域の大きさを調節することができる。そして、ラプラシアンガウシアンフィルタの出力の絶対値が大きくなるのは、2つの領域の差が大きい場合である。すなわち、適切なスケールパラメータσが設定された場合は、左心室中心部と周縁心筋部が比較されて出力が大きくなる。事前知識により、心臓の大きさや心筋の厚さなどが推定可能な場合は、左心室中心を推定するのに最適な大きさであるスケールパラメータσを予め決定する(ステップA2参照)。しかし、スケールパラメータσを一意に決めることができない場合は、複数のσを用意してもよい。
定められたスケールパラメータσで入力断面画像に対してラプラシアンガウシアンフィルタにより処理された出力画像を得る(ステップA3参照)。複数のスケールパラメータσを設定している場合は、各スケールパラメータによって複数の出力画像をそれぞれ得る。
なお、空間フィルタは2つの領域の輝度の比較結果を出力できるようなものであれば良く、例えば、ラプラシアンガウシアンフィルタの代わりに、非特許文献2(David G. Lowe, "Distinctive Image Features from Scale-Invariant Keypoints", International Journal of Computer Vision, vol.60, no.2, pp.91-110, 2004)に記載の差分ガウシアンフィルタ(Difference Of Gaussian : DOG)や、特許第3279913号公報に記載の分離度フィルタなどを用いても、同様の効果が得られる。
(2−3)物体中心推定部130
次に、物体中心推定部130では、得られた出力画像に基づいて関心物体である空洞部の中心位置を推定する。
次に、物体中心推定部130では、得られた出力画像に基づいて関心物体である空洞部の中心位置を推定する。
左心室の中心付近は輝度が低く、周辺部の心筋部分は輝度が高いため、適切なスケールパラメータσを与えた際のラプラシアンガウシアンフィルタによる出力は左心室領域の中心付近で大きくなる。よってここでは、出力画像の各画素の中から、8近傍を比較して画素値が最大となる位置を空洞部の中心候補として取得する(ステップA4参照)。複数の出力画像がある場合は、各出力画像で空洞部の中心候補をそれぞれ取得する。
得られた空洞部の中心候補の中から空洞部の中心を決定する(ステップA4参照)。ここでは、得られた空洞部の中心候補における出力画像の値と、出力画像の中心位置からの距離の2つの値の加重和からなる値が最大となる候補点を中心位置とする。加重和の重み係数は事前の実験から適当な値を求めておく。なお、出力画像の中心位置を加重和に用いるのは、医師が心臓の左心室を撮影する場合には、通常、画像の中心と左心室の中心位置を一致させるか。または、一致させようとするからである。これにより画像の端部近傍の中心候補を除外することができる。
(2−4)初期境界推定部140
次に、初期境界推定部140では、決定された中心位置とフィルタ処理された出力画像を用いて、左心室心筋の内外壁面の境界の初期値を推定する。式(2)に内壁決定時のエネルギー関数を、式(3)に外壁決定時のエネルギー関数を示す。
次に、初期境界推定部140では、決定された中心位置とフィルタ処理された出力画像を用いて、左心室心筋の内外壁面の境界の初期値を推定する。式(2)に内壁決定時のエネルギー関数を、式(3)に外壁決定時のエネルギー関数を示す。
ここで、cは推定された物体中心、rはcを中心とした円の半径、Fminは出力画像の最小値を示す。
図5に示すように、出力画像を利用して、推定された物体中心を中心とする半径rの円形部の線積分から計算されるエネルギーを定義し、定義された内壁、外壁のエネルギーが最小となるrを決定する(ステップA6参照)。
ここで、出力画像が複数ある場合は、全ての出力画像の加重和をとった平均出力画像を用いてエネルギーを計算してもよいし、中心位置が選択された出力画像を代表として用いてもよい。
(2−5)境界決定部150
最後に、境界決定部150では、初期境界位置を用いて最終的な境界位置を決定する(ステップA7参照)。
最後に、境界決定部150では、初期境界位置を用いて最終的な境界位置を決定する(ステップA7参照)。
ここでは、非特許文献3(M. Kass, A. Witkin and D. Terzopoulos, "Snakes : Active Contour Models", International Journal of Computer Vision, 1, pp.321-331, 1988)に記載の動的輪郭モデルを用いる。
動的輪郭モデルによる輪郭抽出結果は初期値に大きく左右されるが、本実施形態で得られる輪郭位置を初期境界として用いることにより、安定した境界抽出が可能となる。また、ここで示した動的輪郭モデル以外にも既存の手法を利用することができ、例えば、特許第3468869号公報に記載の輪郭抽出方法を適用することもできる。
(3)効果
このように第1の実施形態に係わる画像処理装置によれば、入力画像をフィルタ処理して得られる出力画像から空洞部の中心位置を推定し、中心推定に利用した出力画像を用いて初期輪郭推定に必要なエネルギー関数を定義し、得られた中心位置を中心とした円形状を変形させて初期境界を取得し、得られた初期境界を初期値とした動的輪郭モデルを適用することによって、最終的な境界抽出を自動的に行うことができる。
このように第1の実施形態に係わる画像処理装置によれば、入力画像をフィルタ処理して得られる出力画像から空洞部の中心位置を推定し、中心推定に利用した出力画像を用いて初期輪郭推定に必要なエネルギー関数を定義し、得られた中心位置を中心とした円形状を変形させて初期境界を取得し、得られた初期境界を初期値とした動的輪郭モデルを適用することによって、最終的な境界抽出を自動的に行うことができる。
(第2の実施形態)
次に、本発明の第2の実施形態に係わる画像処理装置について図6と図7に基づいて説明する。
次に、本発明の第2の実施形態に係わる画像処理装置について図6と図7に基づいて説明する。
(1)本実施形態の特徴
第1の実施形態では、入力画像をフィルタ処理した1つ以上の出力画像から物体中心候補を取得した。この方法では、フィルタ処理におけるスケールパラメータσが複数設定された場合、出力画像が複数得られ、各出力画像から物体中心候補を取得するため、候補点数が増加する。候補点数が多いほど、正しい中心位置を選択することが難しくなる。また初期境界推定には、適切なスケールパラメータσによる出力画像を用いる方が境界推定を安定に行える。よって、中心位置を決定する前に適切なスケールパラメータσを判断することによって、物体中心推定の失敗を減らし、初期境界推定の精度を向上させることができる。
第1の実施形態では、入力画像をフィルタ処理した1つ以上の出力画像から物体中心候補を取得した。この方法では、フィルタ処理におけるスケールパラメータσが複数設定された場合、出力画像が複数得られ、各出力画像から物体中心候補を取得するため、候補点数が増加する。候補点数が多いほど、正しい中心位置を選択することが難しくなる。また初期境界推定には、適切なスケールパラメータσによる出力画像を用いる方が境界推定を安定に行える。よって、中心位置を決定する前に適切なスケールパラメータσを判断することによって、物体中心推定の失敗を減らし、初期境界推定の精度を向上させることができる。
そこで、図6のブロック図に示すように、本実施形態の画像処理装置は、第1の実施形態における物体中心推定部130に代わり、フィルタ処理された出力画像から物体中心候補を取得する物体中心候補取得部131と、出力画像と中心候補から中心推定に最適な出力画像を選択するスケール評価部132と、スケール評価部132で最適とされた出力画像から得られて中心候補から中心を選択する物体中心決定部133とを備えることを特徴とする。
(2)画像処理装置の動作
次に、図6及び図7を用いて、本実施形態に係わる画像処理装置の動作について説明する。なお、図7は、本実施形態に係わる画像処理装置の動作を示すフローチャートである。
次に、図6及び図7を用いて、本実施形態に係わる画像処理装置の動作について説明する。なお、図7は、本実施形態に係わる画像処理装置の動作を示すフローチャートである。
画像入力部110では、空洞部を含む画像を取得する。ここでは第1の実施形態と同様に、乳頭筋レベルでの胸骨左縁短軸断面像を例として説明する(図7のステップA1参照)。
次に、フィルタ処理部120では、入力画像に対して空間フィルタを施す。空間フィルタとしてラプラシアンガウシアンフィルタを用いる。なお、スケールパラメータσを事前の実験により適当な初期値に定める(ステップA2参照)。
次に、初期値、または、変化したスケールパラメータσで入力画像に対してラプラシアンガウシアンフィルタにより処理された出力画像を得る(ステップA3参照)。
次に、物体中心候補取得部131では、得られた出力画像に基づいて物体の中心位置を推定する。出力画像の各画素の中から、8近傍を比較して画素値が最大となる位置を中心候補として取得する(ステップA4参照)。
次に、スケール評価部132では、物体中心候補取得部131で得られた中心候補の数に基づきスケールパラメータσが適切であるかを判断する。
これは、撮影する空洞部が左心室であって、画像中心付近には低輝度の画素の大きな塊が描写され、左心室の外側に同様の低輝度の塊がごく小数(例えば、左心房や撮影範囲外の画像の端の部分)であることを仮定している。
このような画像の大局的な構造を見るためには、ある程度大きなスケールパラメータσを与えると良いため、ステップA2で与えたスケールパラメータσが条件を満たすまで大きくしていく。
具体的には、中心候補の点数が一定数以上であれば、スケールパラメータσが小さすぎると判断し、ステップB2に移る。中心候補の点数が一定数以下であれば、大局的な構造が得られたと判断し、ステップA5に移る。ここで中心候補数に対する閾値は事前の実験から適当に値に定める(ステップB1参照)。
次に、ステップB2に移った場合は、スケールパラメータσを一定倍してステップA3に戻る。スケール評価部132では、スケールパラメータσが適切であると判断されるまで、スケールパラメータの変更と中心候補の抽出を繰り返す(ステップB2参照)。
次に、物体中心決定部133では、スケール評価部132で適切であると判断されたスケールパラメータσによって得られた出力画像より取得された中心候補の中から中心位置を決定する。その方法は、得られた中心候補における出力画像の値と、出力画像の中心からの距離の2つの値の加重和からなる値が最大となる候補点を物体中心とする。加重和の重み係数は事前の実験から適当な値を求めておく(ステップA5参照)。
次に、初期境界推定部140では、スケール評価部132で最適と判断された出力画像と物体中心決定部133で得た物体中心から、第1の実施形態と同様に初期境界位置を推定する(ステップA6参照)。
最後に、境界決定部150では、第1の実施形態と同様に境界位置を決定する(ステップA7参照)。
(3)効果
このように第2の実施形態に係わる画像処理装置によれば、入力画像に対するフィルタ処理のスケールパラメータを適切な値に定めることができる。
このように第2の実施形態に係わる画像処理装置によれば、入力画像に対するフィルタ処理のスケールパラメータを適切な値に定めることができる。
また、定めたスケールパラメータの空間フィルタより出力画像を取得し、取得した出力画像から空洞部の中心位置を推定し、中心推定に利用した出力画像を用いて初期輪郭推定に必要なエネルギー関数を定義して初期境界を取得し、得られた初期境界を初期値とした動的輪郭モデルを適用することによって、最終的な境界抽出を自動的に行うことができる。
(変更例)
なお、本発明は上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。また、上記実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。さらに、異なる実施形態にわたる構成要素を適宜組み合わせてもよい。
なお、本発明は上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。また、上記実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。さらに、異なる実施形態にわたる構成要素を適宜組み合わせてもよい。
(1)変更例1
第1の実施形態では、心臓の左心室短軸断面像を入力とする場合について述べたが、例えば心尖部四腔断面像において左心室を対象とする場合は、フィルタ処理部120や物体中心推定部130のパラメータを適宜変更し、初期境界推定部140で当てはめたモデル形状を円形から楕円形や、任意の曲線に変更することで本発明を適用できる。
第1の実施形態では、心臓の左心室短軸断面像を入力とする場合について述べたが、例えば心尖部四腔断面像において左心室を対象とする場合は、フィルタ処理部120や物体中心推定部130のパラメータを適宜変更し、初期境界推定部140で当てはめたモデル形状を円形から楕円形や、任意の曲線に変更することで本発明を適用できる。
(2)変更例2
第1の実施形態では、2次元画像である断面像を入力画像とする方法について述べたが、入力が3次元画像の場合は、フィルタ処理部120で3次元の空間フィルタを用い、物体中心推定部130のパラメータを適宜変更し、初期境界推定部140で当てはめたモデル形状を3次元曲面とすることで本実施形態に適用できる。
第1の実施形態では、2次元画像である断面像を入力画像とする方法について述べたが、入力が3次元画像の場合は、フィルタ処理部120で3次元の空間フィルタを用い、物体中心推定部130のパラメータを適宜変更し、初期境界推定部140で当てはめたモデル形状を3次元曲面とすることで本実施形態に適用できる。
(3)変更例3
上記各実施形態では、臓器として心臓で説明したが、これに限らず空洞部を有する臓器であればよく、例えば、血管、胃、子宮などでもよい。
上記各実施形態では、臓器として心臓で説明したが、これに限らず空洞部を有する臓器であればよく、例えば、血管、胃、子宮などでもよい。
110 画像入力部
120 フィルタ処理部
130 物体中心推定部
140 初期境界推定部
150 境界決定部150
120 フィルタ処理部
130 物体中心推定部
140 初期境界推定部
150 境界決定部150
Claims (13)
- 空洞部を有する臓器を撮像した画像を入力する画像入力部と、
前記空洞部に対応する閉領域の画素情報とその閉領域の周囲領域の画素情報を比較して、前記閉領域の中心位置の画素情報を強調して出力する空間フィルタを、前記入力画像に施して出力画像を取得するフィルタ処理部と、
前記出力画像を用いて前記空洞部の中心位置を推定する物体中心推定部と、
前記出力画像と前記推定した中心位置とを用いて、前記空洞部の壁面位置である境界位置を決定する境界決定部と、
を備える画像処理装置。 - 前記フィルタ処理部は、前記空間フィルタとして、前記閉領域と前記周囲領域のそれぞれの明度の加重和の違いを求める空間フィルタを用いる、
請求項1記載の画像処理装置。 - 前記フィルタ処理部は、前記空間フィルタとして、ラプラシアンガウシアンフィルタを用いる、
請求項1記載の画像処理装置。 - 前記フィルタ処理部は、前記空間フィルタとして、差分ガウシアンフィルタを用いる、
請求項1記載の画像処理装置。 - 前記フィルタ処理部は、前記空間フィルタとして、分離度フィルタを用いる、
請求項1記載の画像処理装置。 - 前記フィルタ処理部は、前記空間フィルタにおいて前記閉領域の大きさを表す1つのスケールパラメータを用いて、1つの前記出力画像を取得する、
請求項1記載の画像処理装置。 - 前記フィルタ処理部は、前記空間フィルタにおいて前記閉領域の大きさを表すスケールパラメータを2つ以上用いて、前記出力画像をそれぞれ取得し、
前記物体中心推定部は、それぞれの前記出力画像について中心位置を推定し、
前記境界決定部は、それぞれの前記出力画像とそれぞれの前記中心位置とを用いて、前記境界位置を決定する、
請求項1記載の画像処理装置。 - 前記フィルタ処理部は、前記空間フィルタにおいて前記閉領域の大きさを表すスケールパラメータを2つ以上用いて、前記出力画像をそれぞれ取得し、
前記物体中心推定部は、
前記複数の出力画像のそれぞれから前記空洞部の中心候補をそれぞれ取得する物体中心候補取得部と、
前記各中心候補の数が閾値より小さい出力画像を選択するスケール評価部と、
前記選択した出力画像の中心候補を用いて前記空洞部の中心位置を選択する物体中心決定部と、
を備える請求項1記載の画像処理装置。 - 前記初期境界推定部は、前記選択した出力画像と前記中心位置とを用いて、前記境界位置を決定する、
請求項8記載の画像処理装置。 - 空洞部を有する臓器を撮像した画像を入力する画像入力ステップと、
前記空洞部に対応する閉領域の画素情報とその閉領域の周囲領域の画素情報を比較して、前記閉領域の中心位置の画素情報を強調して出力する空間フィルタを、前記入力画像に施して出力画像を取得するフィルタ処理ステップと、
前記出力画像を用いて前記空洞部の中心位置を推定する物体中心推定ステップと、
前記出力画像と前記推定した中心位置とを用いて、前記空洞部の壁面位置である境界位置を決定する境界決定ステップと、
を備える画像処理方法。 - 前記フィルタ処理ステップにおいて、前記空間フィルタにおいて前記閉領域の大きさを表す1つのスケールパラメータを用いて、1つの前記出力画像を取得する、
請求項10記載の画像処理方法。 - 前記フィルタ処理ステップにおいて、前記空間フィルタにおいて前記閉領域の大きさを表すスケールパラメータを2つ以上用いて、前記出力画像をそれぞれ取得し、
前記物体中心推定ステップにおいて、それぞれの前記出力画像について中心位置を推定し、
前記境界決定ステップにおいて、それぞれの前記出力画像とそれぞれの前記中心位置とを用いて、前記境界位置を決定する、
請求項10記載の画像処理方法。 - 前記フィルタ処理ステップにおいて、前記空間フィルタにおいて前記閉領域の大きさを表すスケールパラメータを2つ以上用いて、前記出力画像をそれぞれ取得し、
前記物体中心推定ステップは、
前記複数の出力画像のそれぞれから前記空洞部の中心候補をそれぞれ取得する物体中心候補取得ステップと、
前記各中心候補の数が閾値より小さい出力画像を選択するスケール評価ステップと、
前記選択した出力画像の中心候補を用いて前記空洞部の中心位置を選択する物体中心決定ステップと、
を備える請求項10記載の画像処理方法。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2007032755A JP2008194239A (ja) | 2007-02-13 | 2007-02-13 | 画像処理装置及びその方法 |
US12/028,954 US20080192998A1 (en) | 2007-02-13 | 2008-02-11 | Image processing apparatus and image processing method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2007032755A JP2008194239A (ja) | 2007-02-13 | 2007-02-13 | 画像処理装置及びその方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2008194239A true JP2008194239A (ja) | 2008-08-28 |
Family
ID=39685851
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2007032755A Pending JP2008194239A (ja) | 2007-02-13 | 2007-02-13 | 画像処理装置及びその方法 |
Country Status (2)
Country | Link |
---|---|
US (1) | US20080192998A1 (ja) |
JP (1) | JP2008194239A (ja) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2014090958A (ja) * | 2012-11-06 | 2014-05-19 | Toshiba Corp | 医用画像処理装置 |
JP2016511043A (ja) * | 2013-03-01 | 2016-04-14 | ボストン サイエンティフィック サイムド,インコーポレイテッドBoston Scientific Scimed,Inc. | 血管内超音波シーケンスにおける管腔境界検出のためのシステム及び方法 |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7460696B2 (en) | 2004-06-01 | 2008-12-02 | Lumidigm, Inc. | Multispectral imaging biometrics |
US20090105578A1 (en) * | 2007-10-19 | 2009-04-23 | Siemens Medical Solutions Usa, Inc. | Interactive Medical Imaging Processing and User Interface System |
US20110105931A1 (en) * | 2007-11-20 | 2011-05-05 | Siemens Medical Solutions Usa, Inc. | System for Determining Patient Heart related Parameters for use in Heart Imaging |
US8553952B2 (en) * | 2008-11-29 | 2013-10-08 | Supratik Kumar Moulik | Method and system for automated interpretation of computer tomography scan data |
US8625902B2 (en) | 2010-07-30 | 2014-01-07 | Qualcomm Incorporated | Object recognition using incremental feature extraction |
WO2012107050A1 (en) | 2011-02-08 | 2012-08-16 | Region Nordjylland, Aalborg Sygehus | A system for determining flow properties of a blood vessel |
US9235890B2 (en) | 2013-03-15 | 2016-01-12 | Yahoo! Inc. | Identifying regions characterized by labeled measurements |
JP6260215B2 (ja) * | 2013-11-13 | 2018-01-17 | 富士通株式会社 | 特徴点追跡装置、及び、特徴点追跡方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH1185996A (ja) * | 1997-09-12 | 1999-03-30 | Toshiba Corp | 分離度フィルタリング処理装置 |
JP2000126182A (ja) * | 1998-10-27 | 2000-05-09 | Mitani Sangyo Co Ltd | 腫瘍診断方法 |
JP2002224116A (ja) * | 2001-01-31 | 2002-08-13 | Matsushita Electric Ind Co Ltd | 超音波診断装置及び画像処理装置 |
JP2007241579A (ja) * | 2006-03-07 | 2007-09-20 | Toshiba Corp | 特徴点検出装置及びその方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5360006A (en) * | 1990-06-12 | 1994-11-01 | University Of Florida Research Foundation, Inc. | Automated method for digital image quantitation |
US5574764A (en) * | 1995-06-06 | 1996-11-12 | General Electric Company | Digital brightness detector |
US5669382A (en) * | 1996-11-19 | 1997-09-23 | General Electric Company | System for measuring myocardium in cardiac images |
EP1037166A1 (de) * | 1999-03-16 | 2000-09-20 | Philips Corporate Intellectual Property GmbH | Verfahren zur Detektion von Konturen in einem Röntgenbild |
US7542622B1 (en) * | 2003-06-02 | 2009-06-02 | The Trustees Of Columbia University In The City Of New York | Spatio-temporal treatment of noisy images using brushlets |
-
2007
- 2007-02-13 JP JP2007032755A patent/JP2008194239A/ja active Pending
-
2008
- 2008-02-11 US US12/028,954 patent/US20080192998A1/en not_active Abandoned
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH1185996A (ja) * | 1997-09-12 | 1999-03-30 | Toshiba Corp | 分離度フィルタリング処理装置 |
JP2000126182A (ja) * | 1998-10-27 | 2000-05-09 | Mitani Sangyo Co Ltd | 腫瘍診断方法 |
JP2002224116A (ja) * | 2001-01-31 | 2002-08-13 | Matsushita Electric Ind Co Ltd | 超音波診断装置及び画像処理装置 |
JP2007241579A (ja) * | 2006-03-07 | 2007-09-20 | Toshiba Corp | 特徴点検出装置及びその方法 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2014090958A (ja) * | 2012-11-06 | 2014-05-19 | Toshiba Corp | 医用画像処理装置 |
JP2016511043A (ja) * | 2013-03-01 | 2016-04-14 | ボストン サイエンティフィック サイムド,インコーポレイテッドBoston Scientific Scimed,Inc. | 血管内超音波シーケンスにおける管腔境界検出のためのシステム及び方法 |
Also Published As
Publication number | Publication date |
---|---|
US20080192998A1 (en) | 2008-08-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP2008194239A (ja) | 画像処理装置及びその方法 | |
US9230331B2 (en) | Systems and methods for registration of ultrasound and CT images | |
US8798344B2 (en) | Image processing apparatus, image processing method and computer-readable recording device | |
KR101017611B1 (ko) | 초음파 - ct 영상 정합을 위한 3차원 초음파 간 영상의 해부학적 특징을 추출하는 시스템 및 방법 | |
Loizou | A review of ultrasound common carotid artery image and video segmentation techniques | |
JP5498299B2 (ja) | 2次元超音波映像に対応する2次元ct映像を提供するシステムおよび方法 | |
JP5579527B2 (ja) | 2次元超音波映像に対応する2次元ct映像を提供するシステムおよび方法 | |
CN104665872B (zh) | 一种基于超声图像的颈动脉内中膜厚度测量方法和装置 | |
CN108090909B (zh) | 一种基于统计偏微分模型的超声造影图像分割方法 | |
JP6361776B2 (ja) | 診断支援装置、及び診断支援装置における画像処理方法、並びにプログラム | |
JP5722414B1 (ja) | 骨粗鬆症診断支援装置 | |
JP6570145B2 (ja) | 画像を処理する方法、プログラム、代替的な投影を構築する方法および装置 | |
WO2009038376A1 (en) | Extraction method of tongue region using graph-based approach and geometric properties | |
US20110058718A1 (en) | Extracting method and apparatus of blood vessel crossing/branching portion | |
CN108601509B (zh) | 图像处理装置、图像处理方法以及记录有程序的介质 | |
JPWO2017086433A1 (ja) | 医用画像処理方法及び装置及びシステム及びプログラム | |
WO2016208016A1 (ja) | 画像処理装置、画像処理方法、及び画像処理プログラム | |
EP3449464B1 (en) | Retinal image processing | |
Alattar et al. | Myocardial segmentation using constrained multi-seeded region growing | |
JP6603709B2 (ja) | 画像処理装置、画像処理方法および画像処理プログラム | |
KR101366341B1 (ko) | 비모수적 확률 모델과 스무딩 함수를 이용한 혈관 내 초음파 영상에서 내강 경계면 자동 분할 방법 | |
JP6082337B2 (ja) | 心胸郭比算出装置 | |
CN112826535A (zh) | 一种超声成像中自动定位血管的方法和装置及设备 | |
JP6099436B2 (ja) | 画像処理装置 | |
US20230245442A1 (en) | Image processing apparatus, image processing method, and non-transitory computer readable medium |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20090326 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20110525 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20110607 |
|
A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20120110 |