(57)【要約】【課題】本発明は、画像処理装置及び画像処理方法を提供する。【解決手段】該画像処理装置は、所定の順に並び替えられた複数のスライス画像を取得する取得部と、該複数のスライス画像に対して検出を順次行い、参考スライス画像及び該参考スライス画像における参考気管領域を決定する選択部と、気管と気管支との分岐点の所在する分岐点スライス画像を決定するまで、該参考気管領域をシード領域とし、領域成長法を用いて該参考スライス画像の後のスライス画像の気管領域を順次決定し、該参考スライス画像の後のスライス画像の気管領域の連結性を順次決定する分岐点決定部と、を含み、該分岐点スライス画像の気管領域は、連結していない２つの領域を含む。本発明の画像処理装置及び画像処理方法によれば、人間の介在を低減でき、分岐点の位置をより正確に決定できる。【選択図】図１
2018134389.tif 000002
【特許請求の範囲】
【請求項１】
所定の順に並び替えられた複数のスライス画像を取得する取得手段と、
前記複数のスライス画像に対して検出を順次行い、参考スライス画像及び前記参考スライス画像における参考気管領域を決定する選択手段と、
気管と気管支との分岐点の所在する分岐点スライス画像を決定するまで、前記参考気管領域をシード領域とし、領域成長法を用いて前記参考スライス画像の後のスライス画像の気管領域を順次決定し、前記参考スライス画像の後のスライス画像の気管領域の連結性を順次決定する分岐点決定手段と、を含み、
前記分岐点スライス画像の気管領域は、連結していない２つの領域を含む、画像処理装置。
【請求項２】
前記選択手段は、
前記参考気管領域を決定するまで、前記複数のスライス画像の各スライス画像について、前記スライス画像の胸郭領域における１つ又は複数の穴領域を順次決定し、前記１つ又は複数の穴領域のうち第１の所定条件を満たす穴領域を前記参考気管領域として決定する参考気管領域決定手段と、
前記参考気管領域の所在するスライス画像を前記参考スライス画像とする参考スライス画像決定手段と、を含む、請求項１に記載の画像処理装置。
【請求項３】
前記第１の所定条件は、前記穴領域の円形度、空気密度及びサイズに関する制約を含む、請求項２に記載の画像処理装置。
【請求項４】
前記分岐点決定手段は、前記分岐点の前記分岐点スライス画像における位置を決定する、請求項１に記載の画像処理装置。
【請求項５】
前記分岐点決定手段は、前記分岐点の前記分岐点スライス画像における位置を決定する際に、
前記分岐点スライス画像の気管領域に含まれる連結していない２つの領域の質量中心をそれぞれ決定し、
２つの質量中心の連結線の中点を前記分岐点の前記分岐点スライス画像における位置として決定する、請求項４に記載の画像処理装置。
【請求項６】
前記分岐点スライス画像の気管領域に含まれる連結していない２つの領域をシード領域とし、領域成長法を用いて前記分岐点スライス画像を初めとして前記分岐点スライス画像、その前のスライス画像及びその後のスライス画像の肺領域を決定する肺決定手段、をさらに含む、請求項１に記載の画像処理装置。
【請求項７】
前記分岐点決定手段は、領域成長法を用いて前記気管領域を決定する際に、第１の閾値を画素差分値の閾値として用い、
前記肺決定手段は、領域成長法を用いて前記肺領域を決定する際に、第２の閾値を画素差分値の閾値として用い、
前記第２の閾値は、前記第１の閾値よりも大きい、請求項６に記載の画像処理装置。
【請求項８】
前記肺決定手段は、
領域成長法を用いて決定された肺領域に基づいて肺の三次元形状を決定し、
前記肺の三次元形状及び前記分岐点の位置に基づいて、左肺及び右肺の三次元形状を決定し、
前記左肺の三次元形状に基づいて前記左肺の上縁及び下縁を決定し、前記右肺の三次元形状に基づいて前記右肺の上縁及び下縁を決定する、請求項６に記載の画像処理装置。
【請求項９】
前記肺決定手段は、前記左肺の上縁の所在するスライス画像、前記左肺の下縁の所在するスライス画像、前記右肺の上縁の所在するスライス画像及び前記右肺の下縁の所在するスライス画像を決定する、請求項８に記載の画像処理装置。
【請求項１０】
所定の順に並び替えられた複数のスライス画像を取得するステップと、
前記複数のスライス画像に対して検出を順次行い、参考スライス画像及び前記参考スライス画像における参考気管領域を決定するステップと、
気管と気管支との分岐点の所在する分岐点スライス画像を決定するまで、前記参考気管領域をシード領域とし、領域成長法を用いて前記参考スライス画像の後のスライス画像の気管領域を順次決定し、前記参考スライス画像の後のスライス画像の気管領域の連結性を順次決定するステップと、を含み、
前記分岐点スライス画像の気管領域は、連結していない２つの領域を含む、画像処理方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明の実施例は、画像処理の分野に関し、具体的には、気管と気管支との分岐点の位置を決定できる画像処理装置及び画像処理方法、並びに肺領域を決定できる画像処理装置及び画像処理方法に関する。
【背景技術】
【０００２】
この部分は、本発明に関連する背景情報を提供するが、必ずしも従来技術ではない。
【０００３】
肺の三次元形状を決定する際に、気管と気管支との分岐点、並びに左肺及び右肺の上縁及び下縁を含む肺の三次元マークポイントを用いる必要がある。従来の気管と気管支との分岐点の検出方法では、肺全体の三次元画像から気管領域を抽出し、細線化アルゴリズムにより気管の骨格構造を解析し、気管と主気管支との接合点を分岐点とする。このような方法では、気管領域を抽出するために、領域成長のためのシード（ｓｅｅｄ）点を手動で選択する必要があり、大量の画像データの処理に適用されない。また、用いられる細線化アルゴリズムは十分なロバスト性を有しないため、気管の骨格構造のノイズに繋がり、分岐点の誤検出の原因となる場合がある。
【０００４】
従来の肺の上縁及び下縁の検出方法では、肺画像から所定範囲内のグレースケール値を有する画素を肺領域として抽出することにより、肺領域の上縁及び下縁を決定する。このような方法では、画素のグレースケール情報のみを考慮し、類似のグレースケール情報を有する他の部分、例えば腸領域も肺領域として認識してしまう可能性があり、このような誤検出された領域は肺の上縁及び下縁の検出に影響を及ぼす。
【０００５】
上記の技術的問題を鑑み、本発明は、分岐点の位置及び肺領域を正確に検出できるものを提供することを目的とする。
【発明の概要】
【発明が解決しようとする課題】
【０００６】
この部分は、本発明の一般的な概要を提供し、その全範囲又はその全ての特徴を完全に開示するものではない。
【０００７】
本発明は、分岐点の位置及び肺領域を正確に検出できる画像処理装置及び画像処理方法を提供することを目的とする。
【課題を解決するための手段】
【０００８】
本発明の１つの態様では、所定の順に並び替えられた複数のスライス画像を取得する取得手段と、前記複数のスライス画像に対して検出を順次行い、参考スライス画像及び前記参考スライス画像における参考気管領域を決定する選択手段と、気管と気管支との分岐点の所在する分岐点スライス画像を決定するまで、前記参考気管領域をシード領域とし、領域成長法を用いて前記参考スライス画像の後のスライス画像の気管領域を順次決定し、前記参考スライス画像の後のスライス画像の気管領域の連結性を順次決定する分岐点決定手段と、を含み、前記分岐点スライス画像の気管領域は、連結していない２つの領域を含む、画像処理装置を提供する。
【０００９】
本発明のもう１つの態様では、所定の順に並び替えられた複数のスライス画像を取得するステップと、前記複数のスライス画像に対して検出を順次行い、参考スライス画像及び前記参考スライス画像における参考気管領域を決定するステップと、気管と気管支との分岐点の所在する分岐点スライス画像を決定するまで、前記参考気管領域をシード領域とし、領域成長法を用いて前記参考スライス画像の後のスライス画像の気管領域を順次決定し、前記参考スライス画像の後のスライス画像の気管領域の連結性を順次決定するステップと、を含み、前記分岐点スライス画像の気管領域は、連結していない２つの領域を含む、画像処理方法を提供する。
【００１０】
本発明のもう１つの態様では、機器が読み取り可能な命令コードを含むプログラムプロダクトであって、前記命令コードがコンピュータにより読み取られ、実行される際に、前記コンピュータに、本発明の画像処理方法を実行させることができる、プログラムプロダクトを提供する。
【００１１】
本発明のもう１つの態様では、機器が読み取り可能な命令コードを含むプログラムプロダクトを記憶する機器読み取り可能な記憶媒体であって、前記命令コードがコンピュータにより読み取られ、実行される際に、前記コンピュータに、本発明の画像処理方法を実行させることができる、記憶媒体を提供する。
【００１２】
本発明の画像処理装置及び画像処理方法によれば、複数のスライス画像から参考スライス画像及びシード領域として用いられる参考気管領域を決定でき、領域成長法を用いて参考スライス画像の後のスライス画像における気管領域の連結性を順次決定でき、分岐点スライス画像を決定できる。これによって、シード領域として用いられる参考気管領域を自動的に決定でき、人間の介在を低減でき、シード領域の正確性を向上できる。さらに、気管領域の連結性に基づいて分岐点を決定することで、より正確な分岐点を決定でき、ノイズの導入を回避でき、後続の肺領域の検出もより正確に行うことができる。
【００１３】
この部分における説明及び特定の例は、単なる例示するためのものであり、本発明の範囲を限定するものではない。
【図面の簡単な説明】
【００１４】
ここで説明される図面は、好ましい実施例を例示するためのものであり、全ての可能な実施例ではなく、本発明の範囲を限定するものではない。
【図１】本発明の実施例に係る画像処理装置の構成を示すブロック図である。
【図２】本発明の実施例に係る所定の順に並び替えられた複数のスライス画像を示す図である。
【図３】本発明の実施例に係る画像処理装置における選択部の構成を示すブロック図である。
【図４】本発明の実施例に係るスライス画像の二値画像を示す図である。
【図５】本発明の実施例に係るスライス画像の胸郭領域における穴領域を示す図である。
【図６】本発明の実施例に係る決定された参考気管領域を示す図である。
【図７】図２におけるＮ１スライス画像を示す図である。
【図８】図２におけるＮ２スライス画像を示す図である。
【図９】図２におけるＮ３スライス画像を示す図である。
【図１０】本発明のもう１つの実施例に係る画像処理装置の構成を示すブロック図である。
【図１１】本発明の実施例に係る決定された気管及び肺の三次元形状を示す図である。
【図１２】本発明の実施例に係る決定された肺の三次元形状を示す図である。
【図１３】本発明の実施例に係る画像処理方法を示すフローチャートである。
【図１４】本発明に係る画像処理方法を実施するための汎用パーソナルコンピュータの例示的な構成を示すブロック図である。  本発明の上記及び他の利点及び特徴を説明するために、以下は、図面を参照しながら本発明の具体的な実施形態をさらに詳細に説明する。該図面及び下記の詳細な説明は本明細書に含まれ、本明細書の一部を形成するものである。同一の機能及び構成を有するユニットは、同一の符号で示されている。なお、これらの図面は、本発明の典型的な例を説明するためのものであり、本発明の範囲を限定するものではない。
【発明を実施するための形態】
【００１５】
以下は、図面を参照しながら本発明の例示的な実施例を詳細に説明する。以下の説明は単なる例示的なものであり、本発明、応用及び用途を限定するものではない。
【００１６】
説明の便宜上、明細書には実際の実施形態の全ての特徴が示されていない。なお、実際に実施する際に、開発者の具体的な目標を実現するために、特定の実施形態を変更してもよく、例えばシステム及び業務に関する制限条件に応じて実施形態を変更してもよい。また、開発作業が非常に複雑であり、且つ時間がかかるが、本公開の当業者にとって、この開発作業は単なる例の作業である。
【００１７】
以下は、図１を参照しながら本発明に係る画像処理装置１００を説明する。
【００１８】
図１に示すように、画像処理装置１００は、取得部１１０、選択部１２０及び分岐点決定部１３０を含む。
【００１９】
本発明の実施例では、取得部１１０は、所定の順に並び替えられた複数のスライス画像を取得してもよい。ここで、取得部１１０は、取得された複数のスライス画像を選択部１２０に送信してもよい。
【００２０】
本発明の実施例では、選択部１２０は、複数のスライス画像に対して検出を順次行い、参考スライス画像及び前記参考スライス画像における参考気管領域を決定してもよい。ここで、選択部１２０は、取得部１１０から複数のスライス画像を取得し、決定された参考スライス画像及び参考気管領域を分岐点決定部１３０に送信してもよい。
【００２１】
本発明の実施例では、分岐点決定部１３０は、気管と気管支との分岐点の所在する分岐点スライス画像を決定するまで、参考気管領域をシード領域とし、領域成長法（ｒｅｇｉｏｎ  ｇｒｏｗｉｎｇ）を用いて参考スライス画像の後のスライス画像の気管領域を順次決定し、参考スライス画像の後のスライス画像の気管領域の連結性を順次決定してもよい。ここで、分岐点スライス画像の気管領域は、連結していない２つの領域を含む。本発明の実施例では、分岐点決定部１３０は、選択部１２０から参考スライス画像及び参考気管領域を取得し、決定された分岐点スライス画像を出力してもよい。
【００２２】
このように、本発明の画像処理装置１００によれば、複数のスライス画像から参考スライス画像及びシード領域として用いられる参考気管領域を決定でき、領域成長法を用いて参考スライス画像の後のスライス画像における気管領域の連結性を順次決定でき、分岐点スライス画像を決定できる。これによって、シード領域として用いられる参考気管領域を自動的に決定でき、人間の介在を低減でき、シード領域の正確性を向上できる。さらに、気管領域の連結性に基づいて分岐点を決定することで、より正確な分岐点を決定でき、ノイズの導入を回避でき、後続の肺領域の検出もより正確に行うことができる。
【００２３】
本発明の実施例では、取得部１１０は、人体の肺に対してコンピュータ断層撮影（ＣＴ：Ｃｏｍｐｕｔｅｒ  Ｔｏｍｏｇｒａｐｈｙ）を行うことで複数のスライス画像を取得し、複数のスライス画像を人体の頭から足までの順に並び替えてもよい。本発明は気管と気管支との分岐点及び肺領域を決定することを目的とするため、複数のスライス画像は、気管、気管支及び肺領域を含む複数のスライス画像、例えば人体の頸から腹部までの部分、より具体的に人体の肺の複数のスライス画像であってもよい。さらに、この複数のスライス画像は、人体の頭から足までの順に並び替えられてもよい。即ち、選択部１２０及び分岐点決定部１３０は、スライス画像に対して検出を順次行う際に、人体の頭から足までの順にスライス画像に対して検出を順次行う。１つの態様では、人体の肺に対してコンピュータ断層撮影を行い、複数のスライス画像を取得し、複数のスライス画像を頭から足までの順に並び替える。
【００２４】
図２は本発明の実施例に係る所定の順に並び替えられた複数のスライス画像を示す図である。図２に示すように、気管領域を含む３つのスライス画像の番号はＮ１、Ｎ２及びＮ３である。ここで、説明の便宜上、３つのスライス画像のみを示しているが、複数のスライス画像の数は３よりも多くてもよい。また、図２に示すスライス画像は何れも気管領域を含むが、気管領域の下方の気管支領域及び肺領域を含むスライス画像が示されていない。
【００２５】
本発明の実施例では、取得部１１０が複数のスライス画像を取得した後に、選択部１２０は、複数のスライス画像に対して検出を順次行い、参考スライス画像及び前記参考スライス画像における参考気管領域を決定してもよい。ここで、選択部１２０は、参考スライス画像及び参考気管領域を決定するまで、スライス画像の並び替えの順に複数のスライス画像に対して検出を順次行ってもよい。即ち、選択部１２０は、参考スライス画像を決定した場合、参考スライス画像の後のスライス画像の検出を行わない。
【００２６】
図３は本発明の実施例に係る画像処理装置１００における選択部１２０の構成を示すブロック図である。図３に示すように、選択部１２０は、参考気管領域決定部１２１及び参考スライス画像決定部１２２を含んでもよい。
【００２７】
本発明の実施例では、参考気管領域決定部１２１は、参考気管領域を決定するまで、複数のスライス画像の各スライス画像について、スライス画像の胸郭領域における１つ又は複数の穴領域を順次決定し、１つ又は複数の穴領域のうち第１の所定条件を満たす穴領域を参考気管領域として決定してもよい。
【００２８】
即ち、参考気管領域決定部１２１は、１枚目のスライス画像を初めとして、スライス画像の胸郭領域における穴領域を順次決定し、第１の所定条件を満たす穴領域が存在するか否かを判定する。第１の所定条件を満たす穴領域があると判定された場合、該穴領域を参考気管領域としてもよい。
【００２９】
本発明の実施例では、参考気管領域決定部１２１は、スライス画像の胸郭領域を決定し、スライス画像の胸郭領域を二値画像に変換してもよい。
【００３０】
図４は本発明の実施例に係るスライス画像の二値画像を示す図である。図４に示すように、白色領域は胸郭領域である。
【００３１】
本発明の実施例では、参考気管領域決定部１２１は、スライス画像の胸郭領域における１つ又は複数の穴領域を決定してもよい。例えば、参考気管領域決定部１２１は、二値画像の胸郭領域における画素値が０に近い領域又は点を穴領域として決定してもよい。例えば、二値画像の胸郭領域において画素値の０との差分値が所定閾値よりも小さい画素を検索し、これらの画素により構成された領域又は点を穴領域とする。ここで、穴領域は１つ又は複数の画素を含んでもよく、穴領域が１つの画素を含む場合、穴領域は実際に点で表され、穴領域が複数の画素を含む場合、穴領域は１つの領域で表される。
【００３２】
図５は本発明の実施例に係るスライス画像の胸郭領域における穴領域を示す図である。図５における白色部分に示すように、図４における胸郭領域は３つの穴領域を含む。
【００３３】
本発明の実施例では、参考気管領域決定部１２１は、１つ又は複数の穴領域に第１の所定条件を満たす穴領域が存在するか否かを決定してもよい。
【００３４】
本発明の実施例では、第１の所定条件は、穴領域の円形度、空気密度及びサイズに関する制約を含む。即ち、参考気管領域決定部１２１は、円形度、空気密度及びサイズが何れも所定条件を満たす穴領域を、第１の所定条件を満たす穴領域として決定し、参考気管領域とする。非制限的な一例として、第１の所定条件を満たす穴領域は、円形度が円形度閾値よりも大きく、空気密度が空気密度閾値よりも大きく、且つサイズが第１サイズ閾値より大きく、第２サイズ閾値よりも小さい。
【００３５】
本発明の実施例では、参考気管領域決定部１２１は、以下の式に従って穴領域の円形度を算出してもよい。
【００３６】
Ｃｉｒｃｕｌａｒｉｔｙ＝４πＳ２／Ｃ２
ここで、Ｃｉｒｃｕｌａｒｉｔｙは穴領域の円形度を表し、Ｓは穴領域の面積を表し、Ｃは穴領域の周囲を表す。
【００３７】
ここで、参考気管領域決定部１２１は、円形度が以下の条件を満たす穴領域を、第１の所定条件を満たす穴領域として決定してもよい。
【００３８】
Ｃｉｒｃｕｌａｒｉｔｙ＞ＴＨＣ
ここで、ＴＨＣは、円形度閾値を表し、実際の要求に応じて設定されてもよい。
【００３９】
本発明の実施例では、参考気管領域決定部１２１は、以下の式に従って穴領域の空気密度を算出してもよい。
【００４０】
Ｄｅｎｓｉｔｙ＝Ｎａｉｒ／Ｎｔｏｔａｌ
ここで、Ｄｅｎｓｉｔｙは穴領域の空気密度を表し、Ｎａｉｒは穴領域における画素値が空気密度範囲内の画素の数を表し、Ｎｔｏｔａｌは穴領域における画素の総数を表す。ここで、実際の要求又は経験値に基づいて空気密度範囲内の画素値を設定し、穴領域における画素値が空気密度範囲内の全ての画素の数を算出し、Ｎａｉｒの値を取得してもよい。さらに、穴領域における全ての画素の数を算出し、Ｎｔｏｔａｌの値を取得してもよい。
【００４１】
ここで、参考気管領域決定部１２１は、空気密度が以下の条件を満たす穴領域を、第１の所定条件を満たす穴領域として決定してもよい。
【００４２】
Ｄｅｎｓｉｔｙ＞ＴＨＤ
ここで、ＴＨＤは、空気密度閾値を表し、実際の要求に応じて設定されてもよい。
【００４３】
本発明の実施例では、参考気管領域決定部１２１は、サイズＳｉｚｅが以下の条件を満たす穴領域を、第１の所定条件を満たす穴領域として決定してもよい。
【００４４】
ＴＨＳ１＞Ｓｉｚｅ＞ＴＨＳ２
ここで、ＴＨＳ１は第２サイズ閾値を表し、ＴＨＳ２は第１サイズ閾値を表し、実際の要求に応じて設定されてもよい。ここで、各種のパラメータ、例えば穴領域の面積、外接円形の半径、外接矩形の長さ及び幅などを用いて穴領域のサイズを表してもよいが、本発明はこれに限定されない。
【００４５】
以上は、参考気管領域決定部１２１による参考気管領域の決定方法の一例を説明しているが、本発明はこの態様に限定されず、参考気管領域決定部１２１は、他の方法を用いて参考気管領域を決定してもよい。上述したように、参考気管領域決定部１２１は、１枚目のスライス画像を初めとして、その後のスライス画像に第１の所定条件を満たす穴領域が存在するか否かを順次決定してもよい。第１の所定条件を満たす穴領域があると決定された場合、該穴領域を参考気管領域として決定する。
【００４６】
本発明の実施例では、参考スライス画像決定部１２２は、参考気管領域の所在するスライス画像を参考スライス画像としてもよい。
【００４７】
図６は本発明の実施例に係る決定された参考気管領域を示す図である。図６に示すように、図５における比較的に大きい穴が参考気管領域として決定された場合は、図５の所在するスライス画像は参考スライス画像である。
【００４８】
本発明の実施例では、選択部１２０が参考気管領域及び参考スライス画像を決定した後に、分岐点決定部１３０は、気管と気管支との分岐点の所在する分岐点スライス画像を決定するまで、参考スライス画像を初めとして参考スライス画像の後のスライス画像について、直前のスライス画像の気管領域をシード領域とし、領域成長法を用いて現在のスライス画像の気管領域を決定し、現在のスライス画像の気管領域の連結性を決定してもよい。
【００４９】
例えば、分岐点決定部１３０は、参考スライス画像における参考気管領域をシード領域とし、領域成長法を用いて参考スライス画像の直後のスライス画像の気管領域を決定し、該気管領域の連結性を決定してもよい。参考スライス画像の直後のスライス画像の気管領域が連結していない２つの領域を含まない場合は、参考スライス画像の直後のスライス画像の気管領域をシード領域とし、領域成長法を用いて次のスライス画像の気管領域を決定し、該気管領域が連結していない２つの領域を含むか否かを決定する。このように、連結していない２つの領域を含む気管領域が見つかるまで上記の処理を行う。そして、分岐点決定部１３０は、取得された連結していない２つの領域を含む気管領域の所在するスライス画像を分岐点スライス画像としてもよい。
【００５０】
以下は、図７〜図９を一例にして分岐点決定部１３０の処理を詳細に説明する。図７、図８及び図９は、図２におけるＮ１スライス画像、Ｎ２スライス画像及びＮ３スライス画像をそれぞれ示す図である。
【００５１】
ここで、Ｎ１スライス画像が参考スライス画像であると仮定し、図７はその参考気管領域を示し、ここの参考気管領域は胸郭に位置する。本発明の実施例では、分岐点決定部１３０は、図８に示すように、参考気管領域をシード領域とし、領域成長法を用いてＮ２スライス画像における気管領域を決定してもよい。図８において、Ｎ２スライス画像における気管領域は左肺と右肺との間に位置する。そして、分岐点決定部１３０は、Ｎ２スライス画像における気管領域の連結性を決定し、該気管領域が連結していない２つの領域を含まないと決定する。そして、分岐点決定部１３０は、図９に示すように、Ｎ２スライス画像における気管領域をシード領域とし、領域成長法を用いてＮ３スライス画像における気管領域を決定する。そして、分岐点決定部１３０は、Ｎ３スライス画像における気管領域の連結性を決定し、該気管領域が連結していない２つの領域を含むと決定する。本発明の実施例では、分岐点決定部１３０は、Ｎ３スライス画像が分岐点スライス画像であると決定する。さらに、分岐点スライス画像を決定した後に、分岐点決定部１３０は、分岐点スライス画像の後のスライス画像に対する検出を停止する。
【００５２】
上述したように、分岐点決定部１３０は、分岐点の所在する分岐点スライス画像を決定してもよい。ここで、例えばスライス画像の番号を用いて分岐点スライス画像を表してもよい。さらに、分岐点決定部１３０は、分岐点の分岐点スライス画像における位置を決定してもよい。
【００５３】
本発明の実施例では、分岐点決定部１３０は、分岐点の分岐点スライス画像における位置を決定する際に、分岐点スライス画像の気管領域に含まれる連結していない２つの領域の質量中心を決定し、２つの質量中心の連結線の中点を分岐点の分岐点スライス画像における位置として決定してもよい。例えば、分岐点決定部１３０は、図９に示す連結していない２つの領域の質量中心を決定し、質量中心の連結線の中点が分岐点の分岐点スライス画像における位置であると決定してもよい。
【００５４】
本発明の実施例では、分岐点決定部１３０は、分岐点の分岐点スライス画像における位置を決定する際に、分岐点スライス画像の気管領域に含まれる連結していない２つの領域における距離が最も近い点を決定し、距離が最も近い点の連結線の中点を分岐点の分岐点スライス画像における位置として決定してもよい。
【００５５】
以上は、分岐点決定部１３０による分岐点の分岐点スライス画像における位置の決定方法の２つの実施例を説明した。なお、この２つの実施例は単なる例示的なものであり、本発明を限定するものではない。分岐点決定部１３０は、他の方法を用いて分岐点の分岐点スライス画像における位置を決定してもよい。ここで、例えば二次元座標を用いて分岐点の分岐点スライス画像における位置を表してもよい。
【００５６】
本発明の画像処理装置１００によれば、分岐点の所在する分岐点スライス画像及び分岐点の分岐点スライス画像における位置を含む分岐点の位置を決定できる。これによって、本発明において決定されたシード領域は選択部１２０により自動的に取得され、人間の指定によるものではないため、人間の介在を低減でき、シード領域の正確性を向上できる。さらに、分岐点決定部１３０は、細線化アルゴリズムを用いることなく、気管及び肺領域全体の骨格構造を決定することなく、分岐点を検出した場合に他のスライス画像の検出を停止することで、より正確な分岐点を決定でき、ノイズの導入を回避でき、後続の肺領域の検出もより正確に行うことができる。
【００５７】
以上は、本発明の実施例に係る気管と気管支との分岐点の位置を決定する画像処理装置１００を詳細に説明した。本発明の実施例では、画像処理装置１００は、肺領域をさらに決定してもよい。
【００５８】
図１０は本発明のもう１つの実施例に係る画像処理装置の構成を示すブロック図である。図１０に示すように、画像処理装置１００は、取得部１１０、選択部１２０、分岐点決定部１３０及び肺決定部１４０を含んでもよい。取得部１１０、選択部１２０及び分岐点決定部１３０は以上で詳細に説明しているため、ここでその説明を省略する。
【００５９】
本発明の実施例では、肺決定部１４０は、分岐点スライス画像の気管領域に含まれる連結していない２つの領域をシード領域とし、領域成長法を用いて分岐点スライス画像を初めとして、その前及びその後のスライス画像（分岐点スライス画像を含む）の肺領域を決定してもよい。
【００６０】
上述したように、分岐点決定部１３０は、分岐点スライス画像の気管領域に含まれる連結していない２つの領域を決定してもよい。本発明の実施例では、肺決定部１４０は、分岐点決定部１３０から、分岐点スライス画像（例えば該スライス画像の番号）及び分岐点スライス画像の気管領域に含まれる連結していない２つの領域を取得してもよい。そして、肺決定部１４０は、この連結していない２つの領域をシード領域とし、領域成長法を用いて分岐点スライス画像を初めとして、スライス画像の並び替えの順及びスライス画像の並び替えの順とは反対する順という２つの方向に沿って、肺領域を決定してもよい。
【００６１】
ここで、肺決定部１４０が領域成長法を用いて肺領域を決定する処理は、分岐点決定部１３０が気管領域を決定する処理と類似する。即ち、直前のスライス画像の肺領域をシード領域として用い、現在のスライス画像の肺領域を検索する。相違点としては、分岐点決定部１３０は、気管領域を決定する際に、参考スライス画像を初めとしてスライス画像の並び替えの順に気管領域を順次決定し、即ちスライス画像の並び替えの順という１つの方法に沿って気管領域を決定し、一方、肺決定部１４０は、肺領域を決定する際に、分岐点スライス画像を初めとして、スライス画像の並び替えの順及びスライス画像の並び替えの順とは反対する順という２つの方向に沿って肺領域を決定する。
【００６２】
また、肺決定部１４０が領域成長法を用いて肺領域を決定する処理と分岐点決定部１３０が気管領域を決定する処理とは、もう１つの相違点があり、即ち領域成長法に用いられる閾値は異なる。
【００６３】
領域成長法の主な原理としては、シード領域の位置及びシード領域に含まれる画素の画素値が既知であり、シード領域周辺の連結領域から画素値がシード領域における画素値と類似する画素を検索し、これらの画素により構成された領域を次のシード領域として検索を行い続ける。本発明に用いられる領域成長法は３Ｄ分野に適用されてもよく、即ち、シード領域周辺の各方向において条件を満たす画素を検索してもよい。通常、閾値を設定することで画素値がシード領域における画素値と類似する画素を判定する。即ち、画素の画素値とシード領域における画素値との差分値が所定閾値よりも小さい場合は、このような画素により構成された領域を次のシード領域としてもよい。従って、閾値は、領域成長法における重要なパラメータである。
【００６４】
本分野では、領域成長法を用いて気管領域を決定する際に、肺組織の密度が気管領域の密度に近いため、肺領域へ「漏れ」やすく、即ち肺領域を気管領域として誤って決定しやすい。この効果は領域成長法の「漏れ」特性と称される。本発明は、領域成長法の「漏れ」特性を利用し、気管領域を決定した後に、分岐点の連結領域をシード領域として、領域成長法を用いて肺領域を決定し続ける。しかし、気管に気体が満ちるが、肺領域に気体の以外に肺胞及び組織など他の成分があるため、気管領域と肺領域とは密度が類似するが、両者は完全に同一ではない。このため、本発明では、領域成長法を用いて気管領域及び肺領域を決定する際に、用いられる閾値を異ならせる必要がある。
【００６５】
本発明の実施例では、分岐点決定部１３０は、領域成長法を用いて気管領域を決定する際に、第１の閾値を画素差分値の閾値として用い、肺決定部１４０は、領域成長法を用いて肺領域を決定する際に、第２の閾値を画素差分値の閾値として用い、第２の閾値は第１の閾値よりも大きい。
【００６６】
即ち、領域成長法を用いて肺領域を決定する場合の判定条件は、気管領域を決定する場合の判定条件よりも緩い。
【００６７】
ここで、分岐点決定部１３０が気管領域を正確に決定し、且つ肺決定部１４０が肺領域を正確に決定するように、第１閾値及び第２閾値を合理的に設定してもよい。好ましくは、第２閾値と第１閾値との差分値は２００以上である。より好ましくは、第１閾値は１００であり、第２閾値は３００である。
【００６８】
上述したように、肺決定部１４０は、領域成長法を用いて、分岐点スライス画像、その前及びその後のスライス画像における肺領域を決定できる。
【００６９】
本発明の実施例では、分岐点決定部１３０は、参考スライス画像と分岐点スライス画像との間の各スライス画像（参考スライス画像及び分岐点スライス画像を含む）における気管領域に基づいて気管の三次元形状を決定してもよく、肺決定部１４０は、分岐点スライス画像、その前及びその後のスライス画像における肺領域に基づいて肺の三次元形状を決定してもよい。さらに、本発明の実施例では、肺決定部１４０は、肺の三次元形状を画像処理装置１００から出力してもよい。
【００７０】
図１１は本発明の実施例に係る決定された気管及び肺の三次元形状を示す図である。図１２は本発明の実施例に係る決定された肺の三次元形状を示す図である。図１２において、気管領域を示しておらず、肺の三次元形状のみを示している。
【００７１】
本発明の実施例では、肺決定部１４０は、肺の三次元形状及び分岐点の位置に基づいて、左肺及び右肺の三次元形状をさらに決定してもよい。即ち、肺決定部１４０は、分岐点の左側に位置する領域を左肺として決定し、分岐点の右側に位置する領域を右肺として決定してもよい。
【００７２】
本発明の実施例では、肺決定部１４０は、左肺の三次元形状に基づいて左肺の上縁及び下縁を決定し、右肺の三次元形状に基づいて右肺の上縁及び下縁を決定してもよい。具体的には、肺決定部１４０は、左肺の上縁の所在するスライス画像、左肺の下縁の所在するスライス画像、右肺の上縁の所在するスライス画像及び右肺の下縁の所在するスライス画像を決定してもよい。ここで、例えばスライス画像の番号を用いて、左肺及び右肺の上縁及び下縁の所在するスライス画像を表してもよい。
【００７３】
上述したように、本発明の画像処理装置１００によれば、肺の三次元形状を決定でき、左肺の上縁の所在するスライス画像、左肺の下縁の所在するスライス画像、右肺の上縁の所在するスライス画像及び右肺の下縁の所在するスライス画像を含む肺の位置を決定できる。これによって、本発明では、領域成長法の「漏れ」特性を用いて肺領域を抽出し、分岐点スライス画像の気管領域をシード領域として肺領域を徐々に検索することで、肺の位置情報を十分に考慮し、他の部分の誤検出を回避し、肺領域をより正確に検出できる。
【００７４】
以上は本発明の実施例に係る画像処理装置１００を詳細に説明した。以下は、本発明の実施例に係る画像処理方法を詳細に説明する。
【００７５】
図１３は本発明の実施例に係る画像処理方法を示すフローチャートである。
【００７６】
図１３に示すように、ステップＳ１３１０において、所定の順に並び替えられた複数のスライス画像を取得する。
【００７７】
そして、ステップＳ１３２０において、複数のスライス画像に対して検出を順次行い、参考スライス画像及び参考スライス画像における参考気管領域を決定する。
【００７８】
そして、ステップＳ１３３０において、気管と気管支との分岐点の所在する分岐点スライス画像を決定するまで、参考気管領域をシード領域とし、領域成長法を用いて参考スライス画像の後のスライス画像の気管領域を順次決定し、参考スライス画像の後のスライス画像の気管領域の連結性を順次決定する。
【００７９】
ここで、分岐点スライス画像の気管領域は、連結していない２つの領域を含む。
【００８０】
好ましくは、参考スライス画像及び前記参考気管領域を決定するステップは、参考気管領域を決定するまで、複数のスライス画像の各スライス画像について、スライス画像の胸郭領域における１つ又は複数の穴領域を順次決定し、１つ又は複数の穴領域のうち第１の所定条件を満たす穴領域を参考気管領域として決定するステップと、参考気管領域の所在するスライス画像を前記参考スライス画像とするステップと、を含む。
【００８１】
好ましくは、第１の所定条件は、穴領域の円形度、空気密度及びサイズに関する制約を含む。
【００８２】
好ましくは、該方法は、分岐点の分岐点スライス画像における位置を決定するステップ、をさらに含む。
【００８３】
好ましくは、分岐点の分岐点スライス画像における位置を決定するステップは、分岐点スライス画像の気管領域に含まれる連結していない２つの領域の質量中心をそれぞれ決定するステップと、２つの質量中心の連結線の中点を分岐点の分岐点スライス画像における位置として決定するステップと、を含む。
【００８４】
好ましくは、該方法は、分岐点スライス画像の気管領域に含まれる連結していない２つの領域をシード領域とし、領域成長法を用いて分岐点スライス画像を初めとして分岐点スライス画像、その前のスライス画像及びその後のスライス画像の肺領域を決定するステップ、をさらに含む。
【００８５】
好ましくは、領域成長法を用いて気管領域を決定する際に、第１の閾値を画素差分値の閾値として用い、領域成長法を用いて肺領域を決定する際に、第２の閾値を画素差分値の閾値として用い、第２の閾値は第１の閾値よりも大きい。
【００８６】
好ましくは、該方法は、領域成長法を用いて決定された肺領域に基づいて肺の三次元形状を決定するステップと、肺の三次元形状及び分岐点の位置に基づいて、左肺及び右肺の三次元形状を決定するステップと、左肺の三次元形状に基づいて左肺の上縁及び下縁を決定し、右肺の三次元形状に基づいて右肺の上縁及び下縁を決定するステップと、をさらに含む。
【００８７】
好ましくは、左肺の上縁及び下縁を決定するステップは、左肺の上縁の所在するスライス画像及び左肺の下縁の所在するスライス画像を決定するステップ、を含み、右肺の上縁及び下縁を決定するステップは、右肺の上縁の所在するスライス画像及び右肺の下縁の所在するスライス画像を決定するステップ、を含む。
【００８８】
上述した画像処理方法は本発明の実施例に係る画像処理装置１００により実現されてもよいため、上述した画像処理装置１００の各態様は全て適用でき、ここでその説明を省略する。
【００８９】
以上から分かるように、本発明の画像処理装置及び画像処理方法によれば、シード領域として用いられる参考気管領域を自動的に決定でき、人間の介在を低減でき、シード領域の正確性を向上できる。さらに、気管領域の連結性に基づいて分岐点を決定することで、より正確な分岐点を決定でき、ノイズの導入を回避でき、後続の肺領域の検出もより正確に行うことができる。また、領域成長法の「漏れ」特性を用いて肺領域を抽出し、分岐点スライス画像の気管領域をシード領域として肺領域を徐々に検索することで、肺の位置情報を十分に考慮し、他の部分の誤検出を回避できる。よって、本発明の画像処理装置及び画像処理方法は、分岐点及び肺領域をより正確に検出できる。
【００９０】
なお、本発明の画像処理方法の各処理は、各種の機器読み取り可能な記憶媒体に記憶されたコンピュータが実行可能なプログラムにより実現されてもよい。
【００９１】
また、本発明の目的は、上記実行可能なプログラムコードを記憶した記憶媒体をシステム又は装置に直接的又は間接的に提供し、該システム又は装置におけるコンピュータ又は中央処理装置（ＣＰＵ）が該プログラムコードを読み出して実行することによって実現されてもよい。この場合は、該システム又は装置はプログラムを実行可能な機能を有すればよく、本発明の実施形態はプログラムに限定されない。また、該プログラムは任意の形であってもよく、例えばオブジェクトプログラム、インタプリタによって実行されるプログラム、又はオペレーティングシステムに提供されるスクリプトプログラム等であってもよい。
【００９２】
上記の機器が読み取り可能な記憶媒体は、各種のメモリ、記憶部、半導体装置、光ディスク、磁気ディスク及び光磁気ディスクのようなディスク、及び他の情報を記憶可能な媒体等を含むが、これらに限定されない。
【００９３】
また、コンピュータがインターネット上の対応するウェブサイトに接続し、本発明のコンピュータプログラムコードをコンピュータにダウンロードしてインストールして実行することによって、本発明の実施形態を実現することもできる。
【００９４】
図１４は本発明に係る画像処理方法を実施するための汎用パーソナルコンピュータの例示的な構成を示すブロック図である。
【００９５】
図１４において、ＣＰＵ１４０１は、読み出し専用メモリ（ＲＯＭ）１４０２に記憶されているプログラム、又は記憶部１４０８からランダムアクセスメモリ（ＲＡＭ）１４０３にロードされたプログラムにより各種の処理を実行する。ＲＡＭ１４０３には、必要に応じて、ＣＰＵ１４０１が各種の処理を実行するに必要なデータが記憶されている。ＣＰＵ１４０１、ＲＯＭ１４０２、及びＲＡＭ１４０３は、バス１４０４を介して互いに接続されている。入力／出力インターフェース１４０５もバス１４０４に接続されている。
【００９６】
入力部１４０６（キーボード、マウスなどを含む）、出力部１４０７（ディスプレイ、例えばブラウン管（ＣＲＴ）、液晶ディスプレイ（ＬＣＤ）など、及びスピーカなどを含む）、記憶部１４０８（例えばハードディスクなどを含む）、通信部１４０９（例えばネットワークのインタフェースカード、例えばＬＡＮカード、モデムなどを含む）は、入力／出力インターフェース１４０５に接続されている。通信部１４０９は、ネットワーク、例えばインターネットを介して通信処理を実行する。必要に応じて、ドライブ部１４１０は、入力／出力インターフェース１４０５に接続されてもよい。取り外し可能な媒体１４１１は、例えば磁気ディスク、光ディスク、光磁気ディスク、半導体メモリなどであり、必要に応じてドライブ部１４１０にセットアップされて、その中から読みだされたコンピュータプログラムは必要に応じて記憶部１４０８にインストールされている。
【００９７】
ソフトウェアにより上記処理を実施する場合、ネットワーク、例えばインターネット、又は記憶媒体、例えば取り外し可能な媒体１４１１を介してソフトウェアを構成するプログラムをインストールする。
【００９８】
なお、これらの記憶媒体は、図１４に示されている、プログラムを記憶し、機器と分離してユーザへプログラムを提供する取り外し可能な媒体１４１１に限定されない。取り外し可能な媒体１４１１は、例えば磁気ディスク（フロッピーディスク（登録商標）を含む）、光ディスク（光ディスク−読み出し専用メモリ（ＣＤ−ＲＯＭ）、及びデジタル多目的ディスク（ＤＶＤ）を含む）、光磁気ディスク（ミニディスク（ＭＤ）（登録商標））及び半導体メモリを含む。或いは、記憶媒体は、ＲＯＭ１４０２、記憶部１４０８に含まれるハードディスクなどであってもよく、プログラムを記憶し、それらを含む機器と共にユーザへ提供される。
【００９９】
なお、本発明のシステム及び方法では、各ユニット又は各ステップを分解且つ、或いは再組み合わせてもよい。これらの分解及び／又は再組み合わせは、本発明と同等であると見なされる。また、本発明の方法は、明細書に説明された時間的順序で実行するものに限定されず、他の時間的順序で順次、並行、又は独立して実行されてもよい。このため、本明細書に説明された方法の実行順序は、本発明の技術的な範囲を限定するものではない。
【０１００】
以上は本発明の具体的な実施例の説明を通じて本発明を開示するが、上記の全ての実施例及び例は例示的なものであり、制限的なものではない。当業者は、特許請求の範囲の主旨及び範囲内で本発明に対して各種の修正、改良、均等的なものに変更してもよい。これらの修正、改良又は均等的なものに変更することは本発明の保護範囲に含まれるものである。
【０１０１】
また、上述の各実施例を含む実施形態に関し、更に以下の付記を開示する。
（付記１）
所定の順に並び替えられた複数のスライス画像を取得する取得手段と、
前記複数のスライス画像に対して検出を順次行い、参考スライス画像及び前記参考スライス画像における参考気管領域を決定する選択手段と、
気管と気管支との分岐点の所在する分岐点スライス画像を決定するまで、前記参考気管領域をシード領域とし、領域成長法を用いて前記参考スライス画像の後のスライス画像の気管領域を順次決定し、前記参考スライス画像の後のスライス画像の気管領域の連結性を順次決定する分岐点決定手段と、を含み、
前記分岐点スライス画像の気管領域は、連結していない２つの領域を含む、画像処理装置。
（付記２）
前記選択手段は、
前記参考気管領域を決定するまで、前記複数のスライス画像の各スライス画像について、前記スライス画像の胸郭領域における１つ又は複数の穴領域を順次決定し、前記１つ又は複数の穴領域のうち第１の所定条件を満たす穴領域を前記参考気管領域として決定する参考気管領域決定手段と、
前記参考気管領域の所在するスライス画像を前記参考スライス画像とする参考スライス画像決定手段と、を含む、付記１に記載の画像処理装置。
（付記３）
前記第１の所定条件は、前記穴領域の円形度、空気密度及びサイズに関する制約を含む、付記２に記載の画像処理装置。
（付記４）
前記分岐点決定手段は、前記分岐点の前記分岐点スライス画像における位置を決定する、付記１に記載の画像処理装置。
（付記５）
前記分岐点決定手段は、前記分岐点の前記分岐点スライス画像における位置を決定する際に、
前記分岐点スライス画像の気管領域に含まれる連結していない２つの領域の質量中心をそれぞれ決定し、
２つの質量中心の連結線の中点を前記分岐点の前記分岐点スライス画像における位置として決定する、付記４に記載の画像処理装置。
（付記６）
前記分岐点スライス画像の気管領域に含まれる連結していない２つの領域をシード領域とし、領域成長法を用いて前記分岐点スライス画像を初めとして前記分岐点スライス画像、その前のスライス画像及びその後のスライス画像の肺領域を決定する肺決定手段、をさらに含む、付記１に記載の画像処理装置。
（付記７）
前記分岐点決定手段は、領域成長法を用いて前記気管領域を決定する際に、第１の閾値を画素差分値の閾値として用い、
前記肺決定手段は、領域成長法を用いて前記肺領域を決定する際に、第２の閾値を画素差分値の閾値として用い、
前記第２の閾値は、前記第１の閾値よりも大きい、付記６に記載の画像処理装置。
（付記８）
前記肺決定手段は、
領域成長法を用いて決定された肺領域に基づいて肺の三次元形状を決定し、
前記肺の三次元形状及び前記分岐点の位置に基づいて、左肺及び右肺の三次元形状を決定し、
前記左肺の三次元形状に基づいて前記左肺の上縁及び下縁を決定し、前記右肺の三次元形状に基づいて前記右肺の上縁及び下縁を決定する、付記６に記載の画像処理装置。
（付記９）
前記肺決定手段は、前記左肺の上縁の所在するスライス画像、前記左肺の下縁の所在するスライス画像、前記右肺の上縁の所在するスライス画像及び前記右肺の下縁の所在するスライス画像を決定する、付記８に記載の画像処理装置。
（付記１０）
前記取得手段は、人体の肺に対してコンピュータ断層撮影を行うことで前記複数のスライス画像を取得する、付記１に記載の画像処理装置。
（付記１１）
所定の順に並び替えられた複数のスライス画像を取得するステップと、
前記複数のスライス画像に対して検出を順次行い、参考スライス画像及び前記参考スライス画像における参考気管領域を決定するステップと、
気管と気管支との分岐点の所在する分岐点スライス画像を決定するまで、前記参考気管領域をシード領域とし、領域成長法を用いて前記参考スライス画像の後のスライス画像の気管領域を順次決定し、前記参考スライス画像の後のスライス画像の気管領域の連結性を順次決定するステップと、を含み、
前記分岐点スライス画像の気管領域は、連結していない２つの領域を含む、画像処理方法。
（付記１２）
前記参考スライス画像及び前記参考気管領域を決定するステップは、
前記参考気管領域を決定するまで、前記複数のスライス画像の各スライス画像について、前記スライス画像の胸郭領域における１つ又は複数の穴領域を順次決定し、前記１つ又は複数の穴領域のうち第１の所定条件を満たす穴領域を前記参考気管領域として決定するステップと、
前記参考気管領域の所在するスライス画像を前記参考スライス画像とするステップと、を含む、付記１１に記載の画像処理方法。
（付記１３）
前記第１の所定条件は、前記穴領域の円形度、空気密度及びサイズに関する制約を含む、付記１２に記載の画像処理方法。
（付記１４）
前記分岐点の前記分岐点スライス画像における位置を決定するステップ、をさらに含む、付記１１に記載の画像処理方法。
（付記１５）
前記分岐点の前記分岐点スライス画像における位置を決定するステップは、
前記分岐点スライス画像の気管領域に含まれる連結していない２つの領域の質量中心をそれぞれ決定するステップと、
２つの質量中心の連結線の中点を前記分岐点の前記分岐点スライス画像における位置として決定するステップと、を含む、付記１４に記載の画像処理方法。
（付記１６）
前記分岐点スライス画像の気管領域に含まれる連結していない２つの領域をシード領域とし、領域成長法を用いて前記分岐点スライス画像を初めとして前記分岐点スライス画像、その前のスライス画像及びその後のスライス画像の肺領域を決定するステップ、をさらに含む、付記１１に記載の画像処理方法。
（付記１７）
領域成長法を用いて前記気管領域を決定する際に、第１の閾値を画素差分値の閾値として用い、
領域成長法を用いて前記肺領域を決定する際に、第２の閾値を画素差分値の閾値として用い、
前記第２の閾値は、前記第１の閾値よりも大きい、付記１６に記載の画像処理方法。
（付記１８）
領域成長法を用いて決定された肺領域に基づいて肺の三次元形状を決定するステップと、
前記肺の三次元形状及び前記分岐点の位置に基づいて、左肺及び右肺の三次元形状を決定するステップと、
前記左肺の三次元形状に基づいて前記左肺の上縁及び下縁を決定し、前記右肺の三次元形状に基づいて前記右肺の上縁及び下縁を決定するステップと、をさらに含む、付記１６に記載の画像処理方法。
（付記１９）
前記左肺の上縁及び下縁を決定するステップは、前記左肺の上縁の所在するスライス画像及び前記左肺の下縁の所在するスライス画像を決定するステップ、を含み、
前記右肺の上縁及び下縁を決定するステップは、前記右肺の上縁の所在するスライス画像及び前記右肺の下縁の所在するスライス画像を決定するステップ、を含む、付記１８に記載の画像処理方法。
（付記２０）
機器が読み取り可能な命令コードを含むプログラムプロダクトを記憶する機器読み取り可能な記憶媒体であって、前記命令コードがコンピュータにより読み取られ、実行される際に、前記コンピュータに、付記１１乃至１９の何れかに記載の画像処理方法を実行させることができる、記憶媒体。
【図１】
2018134389.tif 000003
【図２】
2018134389.tif 000004
【図３】
2018134389.tif 000005
【図４】
2018134389.tif 000006
【図５】
2018134389.tif 000007
【図６】
2018134389.tif 000008
【図７】
2018134389.tif 000009
【図８】
2018134389.tif 000010
【図９】
2018134389.tif 000011
【図１０】
2018134389.tif 000012
【図１１】
2018134389.tif 000013
【図１２】
2018134389.tif 000014
【図１３】
2018134389.tif 000015
【図１４】
2018134389.tif 000016