臓器輪郭抽出方法
【課題】臓器輪郭抽出方法において、位置的に連続した複数枚の各断面画像における臓器の輪郭を抽出する際の処理効率を低下させることなく臓器輪郭をより正確に抽出する。
【解決手段】１枚目の断層画像における臓器の輪郭を任意の方法で抽出した後、Ｎ枚目の断層画像Ｇｎにおけるエッジ検出により臓器１Ａの輪郭の一部を成す輪郭部分の候補となる多数の候補輪郭部分Ｒｎ（１）〜Ｒｎ（ｖ）を検出し、Ｎ−１枚目の断層画像Ｇｎ−１において既に抽出された臓器１Ａの輪郭Ｒｎ−１を形状を保持しつつ変形や移動をさせ、多数の候補輪郭部分Ｒｎ（１）〜Ｒｎ（ｖ）のうち、前記形状を保持しつつ変形や移動をさせた輪郭と略一致する複数の候補輪郭部分Ｒｎ（３）、Ｒｎ（５）、Ｒｎ（６）・・・を含む輪郭をＮ枚目における臓器１Ａの輪郭Ｒｎとして抽出する。
【解決手段】１枚目の断層画像における臓器の輪郭を任意の方法で抽出した後、Ｎ枚目の断層画像Ｇｎにおけるエッジ検出により臓器１Ａの輪郭の一部を成す輪郭部分の候補となる多数の候補輪郭部分Ｒｎ（１）〜Ｒｎ（ｖ）を検出し、Ｎ−１枚目の断層画像Ｇｎ−１において既に抽出された臓器１Ａの輪郭Ｒｎ−１を形状を保持しつつ変形や移動をさせ、多数の候補輪郭部分Ｒｎ（１）〜Ｒｎ（ｖ）のうち、前記形状を保持しつつ変形や移動をさせた輪郭と略一致する複数の候補輪郭部分Ｒｎ（３）、Ｒｎ（５）、Ｒｎ（６）・・・を含む輪郭をＮ枚目における臓器１Ａの輪郭Ｒｎとして抽出する。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は臓器輪郭抽出方法に関し、詳しくは、臓器中の互に異なる断面を表す画像を表す各断層画像それぞれについて臓器の輪郭を抽出する臓器輪郭抽出方法に関するものである。
【背景技術】
【０００２】
従来より、診断対象となる臓器中の同一断面を表す断層画像を時系列的に取得し、各断層画像について臓器の輪郭を効率良く抽出する手法が知られている。この手法は、１枚目の断層画像に示される臓器の輪郭を利用して、２枚目以降の断層画像における臓器の輪郭を抽出するものである（特許文献１参照）。より具体的には、この手法は、１枚目の断層画像における臓器の輪郭（以後、省略して臓器輪郭ともいう）の中心点の位置を、これに続く２枚目の断層画像における対応位置に転写してこの２枚目の断層画像における臓器輪郭の中心点の予想位置を先に仮設定しておき、その後、この仮設定された中心点の位置を基準にして２枚目の断層画像における臓器輪郭を抽出するものである。
【０００３】
また、その手法を、診断対象となる臓器の位置的に連続した互に異なる断面を表す各断層画像（スライス画像ともいう）について臓器輪郭を抽出する際に適用することが検討されている。例えば、１枚目の断層画像における臓器輪郭を任意の方法で抽出した後、この抽出された輪郭の中心点を２枚目の断層画像中の対応位置に転写してこの２枚目の断層画像における臓器輪郭の中心点の予想位置を先に仮設定しておく。そして、この仮設定された中心点の位置を基準にして２枚目の断層画像における臓器輪郭を抽出する方式が検討されている。
【０００４】
なお、診断対象となる臓器の断層画像は、被験者の身長方向に対して直交する断面を表す画像を採用するのが一般的である。
【特許文献１】特開平０７−２４９１１５号公報
【発明の開示】
【発明が解決しようとする課題】
【０００５】
ところで、位置的に連続した各断層画像間において臓器輪郭の大きさや位置は大きく変化する。そのため、１枚目の断層画像中の臓器輪郭の中心点を２枚目の断層画像中の対応位置に転写してこの２枚目の断層画像における臓器輪郭の中心点の予想位置を仮設定する際に大きな位置ずれが生じることがある。
【０００６】
すなわち、１枚目の臓器輪郭の中心点を２枚目の断層画像中の対応位置に転写したときに、この２枚目の断層画像に転写され仮設定された中心点の位置が、２枚目の断層画像における臓器輪郭の正しい中心点の位置から大きくずれたり、あるいは臓器輪郭内の領域から外れたりすることがある。このように、仮設定された中心点に大きな位置ずれが生じると、この位置を基準に抽出される臓器輪郭にも誤差が生じたり、さらに臓器輪郭を抽出できないこともある。
【０００７】
また、２枚目の断層画像への中心点の転写を利用することなく臓器輪郭を抽出しようとすると、画像処理の計算量が増大する等のことにより臓器輪郭を抽出する処理の効率が低下してしまうという問題が生じる。
【０００８】
本発明は、上記事情に鑑みてなされたものであり、臓器輪郭を抽出する際の処理効率を低下させることなくより正確に臓器輪郭を抽出することができる臓器輪郭抽出方法を提供することを目的とするものである。
【課題を解決するための手段】
【０００９】
本発明の臓器輪郭抽出方法は、診断対象となる臓器中の位置的に連続した互に異なる断面を表す複数枚の各断層画像について臓器の輪郭を抽出する際、１枚目の断層画像における臓器の輪郭を任意の方法で抽出した後、２枚目以降の断層画像における臓器の輪郭を、この２枚目以降のＮ枚目の断層画像に隣接する既に輪郭が抽出されたＮ−１枚目の断層画像における輪郭を利用して抽出する臓器輪郭抽出方法において、Ｎ枚目の断層画像におけるエッジ検出により臓器の輪郭の一部を成す輪郭部分の候補となる多数の候補輪郭部分を検出し、Ｎ−１枚目の断層画像において抽出された臓器の輪郭を、形状を保持しつつ変形および／または移動させて、Ｎ枚目の断層画像で抽出された多数の候補輪郭部分のうち、前記変形および／または移動させた臓器の輪郭と略一致する複数の候補輪郭部分を含む輪郭をＮ枚目における臓器の輪郭として抽出することを特徴とするものである。
【００１０】
前記Ｎ−１枚目の断層画像から抽出された臓器の輪郭を含むこの輪郭近傍の領域に対応するＮ枚目の断層画像中の領域を、Ｎ−１枚目の断層画像から抽出された臓器の輪郭を形状を保持しつつ変形および／または移動させる領域とすることが望ましい。すなわち、Ｎ−１枚目の断層画像から抽出された臓器の輪郭を含むこの輪郭の近傍の領域をＮ枚目の断層画像へ転写して形成した領域を、Ｎ−１枚目の断層画像から抽出された臓器の輪郭を形状を保持しつつ変形および／または移動させる領域とすることが望ましい。
【００１１】
なお、Ｎ−１枚目の断層画像中の領域とＮ枚目の断層画像中の領域との対応付けは、例えば、Ｎ−１枚目とＮ枚目の各断層画像中に示される背骨、あるいは体内と外部との境界等を基準にして行うことができる。
【００１２】
前記変形は相似拡大または相似縮小とすることができ、前記移動は平行移動および／または回転とすることができる。
【００１３】
前記診断対象となる臓器中の位置的に連続した互に異なる断面としては、例えば被験者の身長方向や被験者の背骨の延びる方向に対して直交する断面を採用することができる。
【発明の効果】
【００１４】
本発明の臓器輪郭抽出方法によれば、Ｎ枚目の断層画像におけるエッジ検出により臓器の輪郭の一部を成す輪郭部分の候補となる多数の候補輪郭部分を検出し、Ｎ−１枚目の断層画像において抽出された臓器の輪郭を、形状を保持しつつ変形および／または移動させて、Ｎ枚目の断層画像で抽出された多数の候補輪郭部分のうち、前記変形および／または移動させた臓器の輪郭と略一致する複数の候補輪郭部分を含む輪郭をＮ枚目における臓器の輪郭として抽出するようにしたので、臓器輪郭を抽出する際の処理効率を低下させることなくより正確に臓器輪郭を抽出することができる。
【００１５】
すなわち、従来の方式では、各断層画像間における臓器輪郭の大きさや位置の変化等を考慮することなく一方の断層画像の臓器輪郭の中心点を他方の断層画像に転写し、この転写された中心点の予想位置を基準にして他方の断層画像における臓器輪郭の抽出を行なっていた。そのため、中心点の予想位置が正しい位置から大きくずれたときには抽出された臓器輪郭に大きな誤差が生じることがあった。
【００１６】
これに対して、本願発明は、各断層画像間において臓器輪郭の大きさや位置は大きく変化するが、その輪郭の大きさや位置の変化に比して形状の変化は小さく、臓器輪郭の形状には臓器固有の特徴があることに着目して臓器輪郭を抽出するものである。
【００１７】
すなわち、一方の断層画像中から既に抽出された臓器輪郭を、形状を保持しつつ変形および／または移動させてこの臓器輪郭と形状が略一致する複数の候補輪郭部分を含む輪郭を他方の断層画像における臓器輪郭として抽出するようにしたので、各断層画像間における臓器輪郭の大きさや位置の変化の影響を受けることなく他方の断層画像における臓器輪郭の予想位置を抽出することができる。この他方の断層画像における臓器輪郭の予想位置を利用して、例えば他方の断層画像中に臓器輪郭の中心点をより正確に仮設定することができ、この仮設定された中心点を基準にして従来より知られている手法等を適用することにより臓器輪郭を抽出することができる。これにより、臓器輪郭を抽出する際の処理効率を低下させることなくより正確に臓器輪郭を抽出することができる。
【００１８】
また、Ｎ−１枚目の断層画像から抽出された臓器輪郭を含むこの輪郭近傍の領域に対応するＮ枚目の断層画像の領域を、Ｎ−１枚目の断層画像から抽出された臓器輪郭を形状を保持しつつ変形および移動させる領域として定めるようにすれば、不要な領域にまで臓器輪郭を変形させたり移動させたりすることなく臓器輪郭の抽出を実施できるので、臓器輪郭を抽出する際の処理効率を高めることができる。
【００１９】
前記変形を相似拡大または相似縮小としたり、前記移動を平行移動や回転とすれば、より確実に臓器輪郭の抽出を実施することができる。
【発明を実施するための最良の形態】
【００２０】
以下、本発明の実施の形態について、図面を用いて説明する。図１は診断対象となる被験者とこの被験者の臓器中の位置的に連続した互に異なる断面を表す複数枚の各断層画像を示す概念図、図２から図５は断層画像から臓器の輪郭を抽出する様子を示す概念図、図６は臓器輪郭を変形させたり移動させたりする領域を示す図である。
【００２１】
ここで、図２（ａ１）は１枚目の断層画像Ｇ１を示す図、図２（ａ２）は１枚目の断層画像Ｇ１中の臓器１Ａの輪郭Ｒ１を表す輪郭画像Ｈ１を示す図である。図３（ｂ１）はＮ−１枚目の断層画像Ｇｎ−１を示す図、図３（ｂ２）はＮ−１枚目の断層画像Ｇｎ−１中の臓器１Ａの輪郭Ｒｎ−１を表す輪郭画像Ｈｎ−１を示す図である。
【００２２】
また、図４（ｃ１）はＮ枚目の断層画像Ｇｎを示す図、図４（ｃ２）はＮ枚目の断層画像Ｇｎをエッジ検出（輪郭検出）して得たエッジ検出画像Ｊｎを示す図、図４（ｃ３）はＮ枚目の断層画像Ｇｎ中の臓器１Ａの輪郭Ｒｎを表す輪郭画像Ｈｎを示す図である。なお、図４（ｂ２）は図３（ｂ２）と同じ図である。
【００２３】
さらに、図５（ｄ１）はＮ＋１枚目の断層画像Ｇｎ＋１を示す図、図５（ｄ２）はＮ＋１枚目の断層画像Ｇｎ＋１をエッジ検出（輪郭検出）して得たエッジ検出画像Ｊｎ＋１を示す図、図５（ｄ３）はＮ＋１枚目の断層画像Ｇｎ＋１中の臓器１Ａの輪郭Ｒｎ＋１を表す輪郭画像Ｈｎ＋１を示す図である。なお、図５（ｃ３）は図４（ｃ３）と同じ図である。
【００２４】
本発明の実施の形態による臓器輪郭抽出方法は、診断対象となる被験者１の臓器、例えば脾臓１Ａについて位置的に連続した互に異なる断面Ｄ１・・・Ｄｎ-1、Ｄｎ・・・Ｄeを表す複数枚の各断層画像Ｇ１・・・Ｇｎ-1、Ｇｎ・・・Ｇeそれぞれにおける脾臓１Ａの輪郭を抽出するものである。
【００２５】
ここで、診断対象となる臓器中の位置的に連続した互に異なる断面は、各断面間における臓器輪郭の大きさや位置の変化に比してこれらの断面間における臓器輪郭の形状の変化が小さい。例えば、被験者の身長方向や被験者の背骨の延びる方向に対して直交する臓器中の位置的に連続した互に異なる各断面は、各断面間における臓器輪郭の大きさや位置の変化に比してこれらの断面間における臓器輪郭の形状の変化が極めて小さいので、ここでは、被験者の身長方向（背骨の延びる方向）に対して直交する各断面を表す断層画像を採用する。
【００２６】
この臓器輪郭抽出方法は、脾臓１Ａの輪郭を抽出する際、はじめに、１枚目の断層画像Ｇ１（図２(ａ１)参照）における脾臓１Ａの正しい輪郭Ｒ１を任意の方法で抽出する（図２(ａ２)参照）。
【００２７】
その後、２枚目以降の断層画像における脾臓１Ａの輪郭を、２枚目以降のＮ枚目の断層画像Ｇｎに対して位置的に隣接するＮ−１枚目の断層画像Ｇｎ−１（図３(ｂ１)参照）中の脾臓１Ａの臓器輪郭Ｒｎ−１（図３(ｂ２)参照）を利用して抽出するものである。
【００２８】
ここで、この臓器輪郭抽出方法では、Ｎ枚目の断層画像Ｇｎ（図４(ｃ１)参照）におけるエッジ検出（輪郭検出）により、断層画像Ｇｎ中の脾臓１Ａにおける臓器輪郭Ｒｎの一部を示す輪郭部分の候補となる多数の候補輪郭部分Ｒｎ（１）〜Ｒｎ（ｖ）を検出する（図４(ｃ２)参照）。
【００２９】
そして、Ｎ−１枚目の断層画像Ｇｎ−１において既に抽出された脾臓１Ａの臓器輪郭Ｒｎ−１（図３(ｂ２)、図４(ｂ２)参照）を、形状を保持しつつ変形させたり移動させたりして、Ｎ枚目の断層画像Ｇｎから抽出された多数の候補輪郭部分Ｒｎ（１）〜Ｒｎ（ｖ）（図４(ｃ２)参照）のうち上記変形させたり移動させたりした臓器輪郭Ｒｎ−１と略一致する複数の候補輪郭部分Ｒｎ（３）、Ｒｎ（５）、Ｒｎ（６）、Ｒｎ（８）等を含む輪郭をＮ枚目における脾臓１Ａの輪郭Ｒｎとして抽出する（図４(ｃ３)参照）。
【００３０】
つづいて、Ｎ枚目の断層画像Ｇｎにおいて抽出された輪郭画像Ｈｎ中の臓器輪郭Ｒｎを、形状を保持しつつ変形させたり移動させたりして、Ｎ＋１枚目の断層画像Ｇｎ+1（図５（ｄ１）参照）で抽出された多数の候補輪郭部分Ｒｎ+１（１）〜Ｒｎ+１（ｗ）（図５（ｄ２）参照）のうち上記変形させたり移動させたりした臓器輪郭Ｒｎ（図４（ｃ３）、図５（ｃ３）参照）と略一致する複数の候補輪郭部分を含む輪郭をＮ+１枚目における脾臓１Ａの輪郭Ｒｎ+１（図５（ｄ３）参照）として抽出する。
【００３１】
このように、１枚目の断層画像Ｇ１から最後の断層画像Ｇeまで、上記と同様の臓器輪郭を抽出する作用を繰り返すことにより、断層画像Ｇ１から断層画像Ｇeまでの各断層画像において予想されるより正確な脾臓１Ａの輪郭を抽出することができる。
【００３２】
各断層画像について抽出した脾臓１Ａの輪郭は、１枚目の画像で抽出した脾臓１Ａの輪郭の形状と相似な形状を有するものであり、各断層画像における脾臓１Ａの正しい輪郭とは多少のずれがある。
【００３３】
そのため、各断層画像について正確な輪郭が必要な場合には、予想される脾臓１Ａの各輪郭について仮の中心点の位置を求め、その仮の中心点の位置を基準にして脾臓１Ａの輪郭を従来より知られている手法（例えば、特開平０７−２４９１１５号公報参照）を用いて正確に定めることができる。
【００３４】
なお、Ｎ−１枚目の断層画像Ｇｎ−１から抽出された臓器輪郭Ｒｎ−１を、形状を保持しつつ変形させたり移動させたりする領域の範囲を以下のように定めるようにしてもよい。図６(a)は断層画像Ｇｎ−１から抽出された輪郭画像Ｈｎ−１上に臓器輪郭Ｒｎ−１を含むこの臓器輪郭Ｒｎ−１の近傍の領域Ｓｎ−１を合成して表す領域合成画像ＫＨｎ−１示す図、図６(ｂ)は断層画像Ｇｎにエッジ検出を施して得られたエッジ検出画像Ｊｎ上に上記断層画像Ｇｎ−１における領域Ｓｎ−１に対応する領域Ｓｎを合成して表す領域合成画像ＫＪｎを示す図である。
【００３５】
上記臓器輪郭を抽出する際に、Ｎ−１枚目の断層画像Ｇｎ−１から抽出された臓器輪郭Ｒｎ−１を含んだこの臓器輪郭Ｒｎ−１の近傍の領域Ｓｎ−１に対応するＮ枚目の断層画像中の領域Ｓｎを、臓器輪郭Ｒｎ−１を形状を保持しつつ変形させたり移動させたりする領域として定めるようにすることができる。これにより、不要な領域にまで臓器輪郭を変形させたり移動させたりすることなく上記臓器輪郭の抽出を実施することができる。
【００３６】
ここで、Ｎ−１枚目の断層画像中の領域とＮ枚目の断層画像中の領域との対応付けは、例えば、Ｎ−１枚目とＮ枚目の各断層画像中に示される背骨、あるいは体内と外部との境界等を基準にして実施することができる。
【００３７】
なお、本発明における臓器輪郭の変形は、相似拡大および相似縮小とすることができる。また、臓器輪郭の移動は、平行移動または回転とすることができる。
【００３８】
なお、断層画像におけるエッジ検出（輪郭検出）は、この断層画像を差分フィルタ（例えば、Ｓｏｂｅｌフィルタ、ラプラシアンフィルタ）に通すことにより実施することができる。
【００３９】
また、臓器輪郭Ｒｎ−１の形状を保持しつつこの臓器輪郭Ｒｎ−１を変形および移動させる領域Ｓｎを表す基準画像ＪＫｎは、既に抽出された臓器輪郭Ｒｎ−１を表す断層画像である輪郭画像Ｈｎ−１に対してエッジをぼかすスムージング処理を施したりこの輪郭画像Ｈｎ−１に対して距離マップ処理を施すことによってその範囲を定めることができる。
【００４０】
なお、上記スムージング処理には、Gaussian smoothingを採用することができる。また、距離マップ処理には、ユークリッド距離変換を採用することができる。このユークリッド距離変換については以下の文献を参考にすることができる。
【００４１】
「斉藤豊文、鳥脇潤一郎著「３次元デジタル画像に対するユークリッド距離変換」 電子情報通信学会論文誌 J76-D-II, No3 pp.445-453 1993年出版」
ここで、断層画像において抽出された臓器輪郭を、形状を保持しつつ変形および／または移動させて、多数の候補輪郭部分のうち上記変形および／または移動させる臓器輪郭と略一致させる処理においては以下のような手法を採用することができる。すなわち、上記形状を保持しつつ変形および／または移動させるパラメータの最適化について、拡大縮小、平行移動、回転移動等の各パラメータを変化させながら最適なパラメータを求める手法、あるいはPowell法（方向集合を用いた多変数関数の極小値の探索方法）等を採用することができる。このPowell法については以下の文献を参照することができる。
【００４２】
「Brent, Richard P. 1973, Algorithms for Minimization without Derivatives (Englewood Cliffs, N.J.: Prentice-Hall), Chapter 7」
【図面の簡単な説明】
【００４３】
【図１】診断対象となる被験者と被験者の臓器中の位置的に連続した互に異なる断面を表す複数枚の各断層画像を示す概念図
【図２】図２（ａ１）は１枚目の断層画像を示す図、図２（ａ２）は１枚目の断層画像中の臓器の輪郭を表す輪郭画像を示す図
【図３】図３（ｂ１）はＮ−１枚目の断層画像を示す図、図３（ｂ２）はＮ−１枚目の断層画像中の臓器の輪郭を表す輪郭画像を示す図
【図４】図４（ｃ１）はＮ枚目の断層画像を示す図、図４（ｃ２）はＮ枚目の断層画像をエッジ検出して得たエッジ検出画像を示す図、図４（ｃ３）はＮ枚目の断層画像中の臓器の輪郭を表す輪郭画像を示す図、図４（ｂ２）は図３（ｂ２）を転記した図
【図５】図５（ｄ１）はＮ＋１枚目の断層画像を示す図、図５（ｄ２）はＮ＋１枚目の断層画像をエッジ検出して得たエッジ検出画像を示す図、図５（ｄ３）はＮ＋１枚目の断層画像中の臓器の輪郭を表す輪郭画像を示す図、図５（ｃ３）は図４（ｃ３）を転記した図
【図６】図６(a)は断層画像から抽出された輪郭画像上に臓器輪郭を含むこの臓器輪郭近傍の領域を合成して表す領域合成画像示す図、図６(ｂ)は断層画像にエッジ検出を施して得られたエッジ検出画像上に前記領域に対応する領域を合成して表す領域合成画像を示す図
【符号の説明】
【００４４】
１ 被験者
１Ａ 臓器
Ｇｎ−１ Ｎ−１枚目の断層画像
Ｇｎ Ｎ枚目の断層画像
Ｒｎ−１ Ｎ−１枚目における臓器の輪郭
Ｒｎ Ｎ枚目における臓器の輪郭
Ｒｎ（１）〜Ｒｎ（ｖ） 候補輪郭部分
【技術分野】
【０００１】
本発明は臓器輪郭抽出方法に関し、詳しくは、臓器中の互に異なる断面を表す画像を表す各断層画像それぞれについて臓器の輪郭を抽出する臓器輪郭抽出方法に関するものである。
【背景技術】
【０００２】
従来より、診断対象となる臓器中の同一断面を表す断層画像を時系列的に取得し、各断層画像について臓器の輪郭を効率良く抽出する手法が知られている。この手法は、１枚目の断層画像に示される臓器の輪郭を利用して、２枚目以降の断層画像における臓器の輪郭を抽出するものである（特許文献１参照）。より具体的には、この手法は、１枚目の断層画像における臓器の輪郭（以後、省略して臓器輪郭ともいう）の中心点の位置を、これに続く２枚目の断層画像における対応位置に転写してこの２枚目の断層画像における臓器輪郭の中心点の予想位置を先に仮設定しておき、その後、この仮設定された中心点の位置を基準にして２枚目の断層画像における臓器輪郭を抽出するものである。
【０００３】
また、その手法を、診断対象となる臓器の位置的に連続した互に異なる断面を表す各断層画像（スライス画像ともいう）について臓器輪郭を抽出する際に適用することが検討されている。例えば、１枚目の断層画像における臓器輪郭を任意の方法で抽出した後、この抽出された輪郭の中心点を２枚目の断層画像中の対応位置に転写してこの２枚目の断層画像における臓器輪郭の中心点の予想位置を先に仮設定しておく。そして、この仮設定された中心点の位置を基準にして２枚目の断層画像における臓器輪郭を抽出する方式が検討されている。
【０００４】
なお、診断対象となる臓器の断層画像は、被験者の身長方向に対して直交する断面を表す画像を採用するのが一般的である。
【特許文献１】特開平０７−２４９１１５号公報
【発明の開示】
【発明が解決しようとする課題】
【０００５】
ところで、位置的に連続した各断層画像間において臓器輪郭の大きさや位置は大きく変化する。そのため、１枚目の断層画像中の臓器輪郭の中心点を２枚目の断層画像中の対応位置に転写してこの２枚目の断層画像における臓器輪郭の中心点の予想位置を仮設定する際に大きな位置ずれが生じることがある。
【０００６】
すなわち、１枚目の臓器輪郭の中心点を２枚目の断層画像中の対応位置に転写したときに、この２枚目の断層画像に転写され仮設定された中心点の位置が、２枚目の断層画像における臓器輪郭の正しい中心点の位置から大きくずれたり、あるいは臓器輪郭内の領域から外れたりすることがある。このように、仮設定された中心点に大きな位置ずれが生じると、この位置を基準に抽出される臓器輪郭にも誤差が生じたり、さらに臓器輪郭を抽出できないこともある。
【０００７】
また、２枚目の断層画像への中心点の転写を利用することなく臓器輪郭を抽出しようとすると、画像処理の計算量が増大する等のことにより臓器輪郭を抽出する処理の効率が低下してしまうという問題が生じる。
【０００８】
本発明は、上記事情に鑑みてなされたものであり、臓器輪郭を抽出する際の処理効率を低下させることなくより正確に臓器輪郭を抽出することができる臓器輪郭抽出方法を提供することを目的とするものである。
【課題を解決するための手段】
【０００９】
本発明の臓器輪郭抽出方法は、診断対象となる臓器中の位置的に連続した互に異なる断面を表す複数枚の各断層画像について臓器の輪郭を抽出する際、１枚目の断層画像における臓器の輪郭を任意の方法で抽出した後、２枚目以降の断層画像における臓器の輪郭を、この２枚目以降のＮ枚目の断層画像に隣接する既に輪郭が抽出されたＮ−１枚目の断層画像における輪郭を利用して抽出する臓器輪郭抽出方法において、Ｎ枚目の断層画像におけるエッジ検出により臓器の輪郭の一部を成す輪郭部分の候補となる多数の候補輪郭部分を検出し、Ｎ−１枚目の断層画像において抽出された臓器の輪郭を、形状を保持しつつ変形および／または移動させて、Ｎ枚目の断層画像で抽出された多数の候補輪郭部分のうち、前記変形および／または移動させた臓器の輪郭と略一致する複数の候補輪郭部分を含む輪郭をＮ枚目における臓器の輪郭として抽出することを特徴とするものである。
【００１０】
前記Ｎ−１枚目の断層画像から抽出された臓器の輪郭を含むこの輪郭近傍の領域に対応するＮ枚目の断層画像中の領域を、Ｎ−１枚目の断層画像から抽出された臓器の輪郭を形状を保持しつつ変形および／または移動させる領域とすることが望ましい。すなわち、Ｎ−１枚目の断層画像から抽出された臓器の輪郭を含むこの輪郭の近傍の領域をＮ枚目の断層画像へ転写して形成した領域を、Ｎ−１枚目の断層画像から抽出された臓器の輪郭を形状を保持しつつ変形および／または移動させる領域とすることが望ましい。
【００１１】
なお、Ｎ−１枚目の断層画像中の領域とＮ枚目の断層画像中の領域との対応付けは、例えば、Ｎ−１枚目とＮ枚目の各断層画像中に示される背骨、あるいは体内と外部との境界等を基準にして行うことができる。
【００１２】
前記変形は相似拡大または相似縮小とすることができ、前記移動は平行移動および／または回転とすることができる。
【００１３】
前記診断対象となる臓器中の位置的に連続した互に異なる断面としては、例えば被験者の身長方向や被験者の背骨の延びる方向に対して直交する断面を採用することができる。
【発明の効果】
【００１４】
本発明の臓器輪郭抽出方法によれば、Ｎ枚目の断層画像におけるエッジ検出により臓器の輪郭の一部を成す輪郭部分の候補となる多数の候補輪郭部分を検出し、Ｎ−１枚目の断層画像において抽出された臓器の輪郭を、形状を保持しつつ変形および／または移動させて、Ｎ枚目の断層画像で抽出された多数の候補輪郭部分のうち、前記変形および／または移動させた臓器の輪郭と略一致する複数の候補輪郭部分を含む輪郭をＮ枚目における臓器の輪郭として抽出するようにしたので、臓器輪郭を抽出する際の処理効率を低下させることなくより正確に臓器輪郭を抽出することができる。
【００１５】
すなわち、従来の方式では、各断層画像間における臓器輪郭の大きさや位置の変化等を考慮することなく一方の断層画像の臓器輪郭の中心点を他方の断層画像に転写し、この転写された中心点の予想位置を基準にして他方の断層画像における臓器輪郭の抽出を行なっていた。そのため、中心点の予想位置が正しい位置から大きくずれたときには抽出された臓器輪郭に大きな誤差が生じることがあった。
【００１６】
これに対して、本願発明は、各断層画像間において臓器輪郭の大きさや位置は大きく変化するが、その輪郭の大きさや位置の変化に比して形状の変化は小さく、臓器輪郭の形状には臓器固有の特徴があることに着目して臓器輪郭を抽出するものである。
【００１７】
すなわち、一方の断層画像中から既に抽出された臓器輪郭を、形状を保持しつつ変形および／または移動させてこの臓器輪郭と形状が略一致する複数の候補輪郭部分を含む輪郭を他方の断層画像における臓器輪郭として抽出するようにしたので、各断層画像間における臓器輪郭の大きさや位置の変化の影響を受けることなく他方の断層画像における臓器輪郭の予想位置を抽出することができる。この他方の断層画像における臓器輪郭の予想位置を利用して、例えば他方の断層画像中に臓器輪郭の中心点をより正確に仮設定することができ、この仮設定された中心点を基準にして従来より知られている手法等を適用することにより臓器輪郭を抽出することができる。これにより、臓器輪郭を抽出する際の処理効率を低下させることなくより正確に臓器輪郭を抽出することができる。
【００１８】
また、Ｎ−１枚目の断層画像から抽出された臓器輪郭を含むこの輪郭近傍の領域に対応するＮ枚目の断層画像の領域を、Ｎ−１枚目の断層画像から抽出された臓器輪郭を形状を保持しつつ変形および移動させる領域として定めるようにすれば、不要な領域にまで臓器輪郭を変形させたり移動させたりすることなく臓器輪郭の抽出を実施できるので、臓器輪郭を抽出する際の処理効率を高めることができる。
【００１９】
前記変形を相似拡大または相似縮小としたり、前記移動を平行移動や回転とすれば、より確実に臓器輪郭の抽出を実施することができる。
【発明を実施するための最良の形態】
【００２０】
以下、本発明の実施の形態について、図面を用いて説明する。図１は診断対象となる被験者とこの被験者の臓器中の位置的に連続した互に異なる断面を表す複数枚の各断層画像を示す概念図、図２から図５は断層画像から臓器の輪郭を抽出する様子を示す概念図、図６は臓器輪郭を変形させたり移動させたりする領域を示す図である。
【００２１】
ここで、図２（ａ１）は１枚目の断層画像Ｇ１を示す図、図２（ａ２）は１枚目の断層画像Ｇ１中の臓器１Ａの輪郭Ｒ１を表す輪郭画像Ｈ１を示す図である。図３（ｂ１）はＮ−１枚目の断層画像Ｇｎ−１を示す図、図３（ｂ２）はＮ−１枚目の断層画像Ｇｎ−１中の臓器１Ａの輪郭Ｒｎ−１を表す輪郭画像Ｈｎ−１を示す図である。
【００２２】
また、図４（ｃ１）はＮ枚目の断層画像Ｇｎを示す図、図４（ｃ２）はＮ枚目の断層画像Ｇｎをエッジ検出（輪郭検出）して得たエッジ検出画像Ｊｎを示す図、図４（ｃ３）はＮ枚目の断層画像Ｇｎ中の臓器１Ａの輪郭Ｒｎを表す輪郭画像Ｈｎを示す図である。なお、図４（ｂ２）は図３（ｂ２）と同じ図である。
【００２３】
さらに、図５（ｄ１）はＮ＋１枚目の断層画像Ｇｎ＋１を示す図、図５（ｄ２）はＮ＋１枚目の断層画像Ｇｎ＋１をエッジ検出（輪郭検出）して得たエッジ検出画像Ｊｎ＋１を示す図、図５（ｄ３）はＮ＋１枚目の断層画像Ｇｎ＋１中の臓器１Ａの輪郭Ｒｎ＋１を表す輪郭画像Ｈｎ＋１を示す図である。なお、図５（ｃ３）は図４（ｃ３）と同じ図である。
【００２４】
本発明の実施の形態による臓器輪郭抽出方法は、診断対象となる被験者１の臓器、例えば脾臓１Ａについて位置的に連続した互に異なる断面Ｄ１・・・Ｄｎ-1、Ｄｎ・・・Ｄeを表す複数枚の各断層画像Ｇ１・・・Ｇｎ-1、Ｇｎ・・・Ｇeそれぞれにおける脾臓１Ａの輪郭を抽出するものである。
【００２５】
ここで、診断対象となる臓器中の位置的に連続した互に異なる断面は、各断面間における臓器輪郭の大きさや位置の変化に比してこれらの断面間における臓器輪郭の形状の変化が小さい。例えば、被験者の身長方向や被験者の背骨の延びる方向に対して直交する臓器中の位置的に連続した互に異なる各断面は、各断面間における臓器輪郭の大きさや位置の変化に比してこれらの断面間における臓器輪郭の形状の変化が極めて小さいので、ここでは、被験者の身長方向（背骨の延びる方向）に対して直交する各断面を表す断層画像を採用する。
【００２６】
この臓器輪郭抽出方法は、脾臓１Ａの輪郭を抽出する際、はじめに、１枚目の断層画像Ｇ１（図２(ａ１)参照）における脾臓１Ａの正しい輪郭Ｒ１を任意の方法で抽出する（図２(ａ２)参照）。
【００２７】
その後、２枚目以降の断層画像における脾臓１Ａの輪郭を、２枚目以降のＮ枚目の断層画像Ｇｎに対して位置的に隣接するＮ−１枚目の断層画像Ｇｎ−１（図３(ｂ１)参照）中の脾臓１Ａの臓器輪郭Ｒｎ−１（図３(ｂ２)参照）を利用して抽出するものである。
【００２８】
ここで、この臓器輪郭抽出方法では、Ｎ枚目の断層画像Ｇｎ（図４(ｃ１)参照）におけるエッジ検出（輪郭検出）により、断層画像Ｇｎ中の脾臓１Ａにおける臓器輪郭Ｒｎの一部を示す輪郭部分の候補となる多数の候補輪郭部分Ｒｎ（１）〜Ｒｎ（ｖ）を検出する（図４(ｃ２)参照）。
【００２９】
そして、Ｎ−１枚目の断層画像Ｇｎ−１において既に抽出された脾臓１Ａの臓器輪郭Ｒｎ−１（図３(ｂ２)、図４(ｂ２)参照）を、形状を保持しつつ変形させたり移動させたりして、Ｎ枚目の断層画像Ｇｎから抽出された多数の候補輪郭部分Ｒｎ（１）〜Ｒｎ（ｖ）（図４(ｃ２)参照）のうち上記変形させたり移動させたりした臓器輪郭Ｒｎ−１と略一致する複数の候補輪郭部分Ｒｎ（３）、Ｒｎ（５）、Ｒｎ（６）、Ｒｎ（８）等を含む輪郭をＮ枚目における脾臓１Ａの輪郭Ｒｎとして抽出する（図４(ｃ３)参照）。
【００３０】
つづいて、Ｎ枚目の断層画像Ｇｎにおいて抽出された輪郭画像Ｈｎ中の臓器輪郭Ｒｎを、形状を保持しつつ変形させたり移動させたりして、Ｎ＋１枚目の断層画像Ｇｎ+1（図５（ｄ１）参照）で抽出された多数の候補輪郭部分Ｒｎ+１（１）〜Ｒｎ+１（ｗ）（図５（ｄ２）参照）のうち上記変形させたり移動させたりした臓器輪郭Ｒｎ（図４（ｃ３）、図５（ｃ３）参照）と略一致する複数の候補輪郭部分を含む輪郭をＮ+１枚目における脾臓１Ａの輪郭Ｒｎ+１（図５（ｄ３）参照）として抽出する。
【００３１】
このように、１枚目の断層画像Ｇ１から最後の断層画像Ｇeまで、上記と同様の臓器輪郭を抽出する作用を繰り返すことにより、断層画像Ｇ１から断層画像Ｇeまでの各断層画像において予想されるより正確な脾臓１Ａの輪郭を抽出することができる。
【００３２】
各断層画像について抽出した脾臓１Ａの輪郭は、１枚目の画像で抽出した脾臓１Ａの輪郭の形状と相似な形状を有するものであり、各断層画像における脾臓１Ａの正しい輪郭とは多少のずれがある。
【００３３】
そのため、各断層画像について正確な輪郭が必要な場合には、予想される脾臓１Ａの各輪郭について仮の中心点の位置を求め、その仮の中心点の位置を基準にして脾臓１Ａの輪郭を従来より知られている手法（例えば、特開平０７−２４９１１５号公報参照）を用いて正確に定めることができる。
【００３４】
なお、Ｎ−１枚目の断層画像Ｇｎ−１から抽出された臓器輪郭Ｒｎ−１を、形状を保持しつつ変形させたり移動させたりする領域の範囲を以下のように定めるようにしてもよい。図６(a)は断層画像Ｇｎ−１から抽出された輪郭画像Ｈｎ−１上に臓器輪郭Ｒｎ−１を含むこの臓器輪郭Ｒｎ−１の近傍の領域Ｓｎ−１を合成して表す領域合成画像ＫＨｎ−１示す図、図６(ｂ)は断層画像Ｇｎにエッジ検出を施して得られたエッジ検出画像Ｊｎ上に上記断層画像Ｇｎ−１における領域Ｓｎ−１に対応する領域Ｓｎを合成して表す領域合成画像ＫＪｎを示す図である。
【００３５】
上記臓器輪郭を抽出する際に、Ｎ−１枚目の断層画像Ｇｎ−１から抽出された臓器輪郭Ｒｎ−１を含んだこの臓器輪郭Ｒｎ−１の近傍の領域Ｓｎ−１に対応するＮ枚目の断層画像中の領域Ｓｎを、臓器輪郭Ｒｎ−１を形状を保持しつつ変形させたり移動させたりする領域として定めるようにすることができる。これにより、不要な領域にまで臓器輪郭を変形させたり移動させたりすることなく上記臓器輪郭の抽出を実施することができる。
【００３６】
ここで、Ｎ−１枚目の断層画像中の領域とＮ枚目の断層画像中の領域との対応付けは、例えば、Ｎ−１枚目とＮ枚目の各断層画像中に示される背骨、あるいは体内と外部との境界等を基準にして実施することができる。
【００３７】
なお、本発明における臓器輪郭の変形は、相似拡大および相似縮小とすることができる。また、臓器輪郭の移動は、平行移動または回転とすることができる。
【００３８】
なお、断層画像におけるエッジ検出（輪郭検出）は、この断層画像を差分フィルタ（例えば、Ｓｏｂｅｌフィルタ、ラプラシアンフィルタ）に通すことにより実施することができる。
【００３９】
また、臓器輪郭Ｒｎ−１の形状を保持しつつこの臓器輪郭Ｒｎ−１を変形および移動させる領域Ｓｎを表す基準画像ＪＫｎは、既に抽出された臓器輪郭Ｒｎ−１を表す断層画像である輪郭画像Ｈｎ−１に対してエッジをぼかすスムージング処理を施したりこの輪郭画像Ｈｎ−１に対して距離マップ処理を施すことによってその範囲を定めることができる。
【００４０】
なお、上記スムージング処理には、Gaussian smoothingを採用することができる。また、距離マップ処理には、ユークリッド距離変換を採用することができる。このユークリッド距離変換については以下の文献を参考にすることができる。
【００４１】
「斉藤豊文、鳥脇潤一郎著「３次元デジタル画像に対するユークリッド距離変換」 電子情報通信学会論文誌 J76-D-II, No3 pp.445-453 1993年出版」
ここで、断層画像において抽出された臓器輪郭を、形状を保持しつつ変形および／または移動させて、多数の候補輪郭部分のうち上記変形および／または移動させる臓器輪郭と略一致させる処理においては以下のような手法を採用することができる。すなわち、上記形状を保持しつつ変形および／または移動させるパラメータの最適化について、拡大縮小、平行移動、回転移動等の各パラメータを変化させながら最適なパラメータを求める手法、あるいはPowell法（方向集合を用いた多変数関数の極小値の探索方法）等を採用することができる。このPowell法については以下の文献を参照することができる。
【００４２】
「Brent, Richard P. 1973, Algorithms for Minimization without Derivatives (Englewood Cliffs, N.J.: Prentice-Hall), Chapter 7」
【図面の簡単な説明】
【００４３】
【図１】診断対象となる被験者と被験者の臓器中の位置的に連続した互に異なる断面を表す複数枚の各断層画像を示す概念図
【図２】図２（ａ１）は１枚目の断層画像を示す図、図２（ａ２）は１枚目の断層画像中の臓器の輪郭を表す輪郭画像を示す図
【図３】図３（ｂ１）はＮ−１枚目の断層画像を示す図、図３（ｂ２）はＮ−１枚目の断層画像中の臓器の輪郭を表す輪郭画像を示す図
【図４】図４（ｃ１）はＮ枚目の断層画像を示す図、図４（ｃ２）はＮ枚目の断層画像をエッジ検出して得たエッジ検出画像を示す図、図４（ｃ３）はＮ枚目の断層画像中の臓器の輪郭を表す輪郭画像を示す図、図４（ｂ２）は図３（ｂ２）を転記した図
【図５】図５（ｄ１）はＮ＋１枚目の断層画像を示す図、図５（ｄ２）はＮ＋１枚目の断層画像をエッジ検出して得たエッジ検出画像を示す図、図５（ｄ３）はＮ＋１枚目の断層画像中の臓器の輪郭を表す輪郭画像を示す図、図５（ｃ３）は図４（ｃ３）を転記した図
【図６】図６(a)は断層画像から抽出された輪郭画像上に臓器輪郭を含むこの臓器輪郭近傍の領域を合成して表す領域合成画像示す図、図６(ｂ)は断層画像にエッジ検出を施して得られたエッジ検出画像上に前記領域に対応する領域を合成して表す領域合成画像を示す図
【符号の説明】
【００４４】
１ 被験者
１Ａ 臓器
Ｇｎ−１ Ｎ−１枚目の断層画像
Ｇｎ Ｎ枚目の断層画像
Ｒｎ−１ Ｎ−１枚目における臓器の輪郭
Ｒｎ Ｎ枚目における臓器の輪郭
Ｒｎ（１）〜Ｒｎ（ｖ） 候補輪郭部分
【特許請求の範囲】
【請求項１】
診断対象となる臓器中の位置的に連続した互に異なる断面を表す複数枚の各断層画像について前記臓器の輪郭を抽出する際、１枚目の断層画像における前記臓器の輪郭を任意の方法で抽出した後、
２枚目以降の断層画像における前記臓器の輪郭を、該２枚目以降のＮ枚目の断層画像に隣接するＮ−１枚目の断層画像において抽出された前記臓器の輪郭を利用して抽出する臓器輪郭抽出方法において、
前記Ｎ枚目の断層画像におけるエッジ検出により前記臓器の輪郭の一部を成す輪郭部分の候補となる多数の候補輪郭部分を検出し、
前記Ｎ−１枚目の断層画像において抽出された前記臓器の輪郭を、形状を保持しつつ変形および／または移動させて、前記Ｎ枚目の断層画像で抽出された多数の候補輪郭部分のうち、前記変形および／または移動させた前記臓器の輪郭と略一致する複数の候補輪郭部分を含む輪郭をＮ枚目における前記臓器の輪郭として抽出することを特徴とする臓器輪郭抽出方法。
【請求項２】
前記Ｎ−１枚目の断層画像から抽出された前記臓器の輪郭を含む該輪郭近傍の領域に対応する前記Ｎ枚目の断層画像中の領域を、前記形状を保持しつつ変形および／または移動させる領域とすることを特徴とする請求項１記載の臓器輪郭抽出方法。
【請求項３】
前記変形が、相似拡大または相似縮小であることを特徴とする請求項１または２記載の臓器輪郭抽出方法。
【請求項４】
前記移動が、平行移動および／または回転であることを特徴とする請求項１から３のいずれか１項記載の臓器輪郭抽出方法。
【請求項１】
診断対象となる臓器中の位置的に連続した互に異なる断面を表す複数枚の各断層画像について前記臓器の輪郭を抽出する際、１枚目の断層画像における前記臓器の輪郭を任意の方法で抽出した後、
２枚目以降の断層画像における前記臓器の輪郭を、該２枚目以降のＮ枚目の断層画像に隣接するＮ−１枚目の断層画像において抽出された前記臓器の輪郭を利用して抽出する臓器輪郭抽出方法において、
前記Ｎ枚目の断層画像におけるエッジ検出により前記臓器の輪郭の一部を成す輪郭部分の候補となる多数の候補輪郭部分を検出し、
前記Ｎ−１枚目の断層画像において抽出された前記臓器の輪郭を、形状を保持しつつ変形および／または移動させて、前記Ｎ枚目の断層画像で抽出された多数の候補輪郭部分のうち、前記変形および／または移動させた前記臓器の輪郭と略一致する複数の候補輪郭部分を含む輪郭をＮ枚目における前記臓器の輪郭として抽出することを特徴とする臓器輪郭抽出方法。
【請求項２】
前記Ｎ−１枚目の断層画像から抽出された前記臓器の輪郭を含む該輪郭近傍の領域に対応する前記Ｎ枚目の断層画像中の領域を、前記形状を保持しつつ変形および／または移動させる領域とすることを特徴とする請求項１記載の臓器輪郭抽出方法。
【請求項３】
前記変形が、相似拡大または相似縮小であることを特徴とする請求項１または２記載の臓器輪郭抽出方法。
【請求項４】
前記移動が、平行移動および／または回転であることを特徴とする請求項１から３のいずれか１項記載の臓器輪郭抽出方法。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図２】
【図３】
【図４】
【図５】
【図６】
【公開番号】特開２００８−２５９６７７（Ｐ２００８−２５９６７７Ａ）
【公開日】平成２０年１０月３０日（２００８．１０．３０）
【国際特許分類】
【出願番号】特願２００７−１０４８１０（Ｐ２００７−１０４８１０）
【出願日】平成１９年４月１２日（２００７．４．１２）
【出願人】（３０６０３７３１１）富士フイルム株式会社 (25,513)
【Ｆターム（参考）】
【公開日】平成２０年１０月３０日（２００８．１０．３０）
【国際特許分類】
【出願日】平成１９年４月１２日（２００７．４．１２）
【出願人】（３０６０３７３１１）富士フイルム株式会社 (25,513)
【Ｆターム（参考）】
[ Back to top ]