(57)【特許請求の範囲】
【請求項１】
頭部の医用画像において脳出血部位をセグメンテーションする装置の作動方法であって、
前記頭部の医用画像において、頭骨の内側領域を区分する工程と、
該頭骨の内側領域を区分する工程において区分された領域の中から、脳出血部位が含まれる可能性のある領域を区分する工程と、 
該脳出血部位が含まれる可能性のある領域を区分する工程により区分された領域の中から、脳出血部位を判定する工程と
を備え、
前記脳出血部位が含まれる可能性のある領域を区分する工程は、
前記頭骨の内側領域を区分する工程において区分された領域のうち、ＣＴ値が緩やかに変化する領域を区分する工程と、
該ＣＴ値が緩やかに変化する領域を区分する工程において区分された領域のうち、画素数が所定の数よりも大きい領域を区分する工程と、 
該画素数が所定の数よりも大きい領域を区分する工程において区分された領域から、比較的ＣＴ値の大きい領域を区分する工程と、 
該比較的ＣＴ値の大きい領域を区分する工程において抽出された領域について、ＣＴ値が緩やかに変化する領域を区分する工程と、 
を含んでいる脳出血部位をセグメンテーションする装置の作動方法。
【請求項２】
前記脳出血部位を判定する工程は、
前記脳出血部位が含まれる可能性のある領域を区分する工程において区分された領域について、該領域のＣＴ値及び大きさに基づいて、該領域が脳出血部位か否かを判定する第１の工程と、 
該第１の工程において脳出血部位と判定された領域以外の領域のうち、脳出血部位の可能性がある領域について、該領域とその周辺の領域とのＣＴ値の差又は該領域のＣＴ値に基づいて、該領域が脳出血部位か否かを判定する第２の工程と、 
を含んでいる請求項１に記載の脳出血部位をセグメンテーションする装置の作動方法。
【請求項３】
前記脳出血部位を判定する工程において脳出血部位であると判定された領域についてバーシャルボリューム効果の影響を補正する工程をさらに含んでいる請求項１に記載の脳出血部位をセグメンテーションする装置の作動方法。 
【請求項４】
頭部のＸ線ＣＴ画像について脳出血部位をセグメンテーションする方法であって、
前処理工程と、１次探索工程と、分析・調整工程と、フィルタリング工程と、２次探索工程と、分析・判定工程と、後処理工程とを有し、 
前記前処理工程では、頭部のＸ線ＣＴ画像について、
第１の設定値より大きいＣＴ値を持つ画素と第２の設定値より小さいＣＴ値を持つ画素を除外し、
第３の設定値に基づいて頭骨の境界を特定し、
前記特定した境界に基づいて頭骨の外側領域を除外し、
前記１次探索工程では、前記前処理が行われた画像について、
ＣＴ値変化が緩やかな領域を探索し、
前記探索により判明した領域ごとにラベルを付与し、
前記分析・調整工程では、前記１次探索が行われた画像について、
前記領域ごとに画素数を求め、
画素数が第４の設定値より少ない領域を除外し、
前記フィルタリング工程では、前記分析・調整が行われた画像について、
前記領域ごとに、
２次元座標ｉ，ｊが増加する方向にそれぞれ隣接する画素間のＣＴ値の差の絶対値の和を画素ごとに求め、
前記和が第５の設定値以上となる画素位置を特定し、
前記特定した全ての画素位置におけるＣＴ値の平均値を求め、
前記平均値以下のＣＴ値を持つ画素を除外し、
前記２次探索工程では、前記フィルタリングが行われた画像について、
ＣＴ値変化が緩やかな領域を探索し、
前記探索により判明した領域ごとにあらためてラベルを付与し、
前記分析・判定工程では、前記２次探索においてラベルが付与された画像について、
前記領域の画素のＣＴ値をCTregion、第６の設定値をCTmax、第７の設定値をCTmin、第８の設定値をCTbloodとしたとき、 
第１のインデックスIndexCTを、 
CTmin＜=CTregion＜=CTbloodならば
【数１】
0005383486.tif 000030
とし、
CTblood＜ CTregion＜= CTmaxならば
【数２】
0005383486.tif 000031
とし、
それ以外は
【数３】
0005383486.tif 000032
とし、
前記領域の面積および周長をそれぞれAreaおよびPerimeterで表し、前記領域の特性値を 
【数４】
0005383486.tif 000033
としたとき、
第２のインデックスindexradianを、 
Radianが第９の設定値Radianmaxより大きければindexradian =1とし、
Radianが第１０の設定値Radianminより小さければindexradian =0とし、
Radianが第９の設定値Radianmax以下で第１０の設定値Radianmin以上ならば 
【数５】
0005383486.tif 000034
とし、
第３のインデックスIndex0を、
【数６】
0005383486.tif 000035
とし、
Index0＞=20%ならば、その領域は脳出血部位であると判定し、 
Index0＜=3%ならば、その領域は脳出血部位でないと判定し、 
3%＜Index0＜20%ならば、そのような領域の周囲の領域の画素のＣＴ値をCTAroundRegionとし、そのような全ての領域の平均ＣＴ値をCTAverage -All-Regionとしたとき、
第４のインデックスIndexSubを、 
【数７】
0005383486.tif 000036
とし、
第５のインデックスIndexOderを、 
【数８】
0005383486.tif 000037
とし、
IndexSub＞=8またはIndexOrder＞ =10ならば、その領域は脳出血部位であると判定し、
IndexSub＜=0またはIndexOrder＜ =-5ならば、その領域は脳出血部位でないと判定し、
0＜IndexSub＜ 8かつ-5＜ IndexOrder＜10ならば、 
第６のインデックスIndexFinalを、 
【数９】
0005383486.tif 000038
とし、
IndexFinal＞=50%ならば、その領域は脳出血部位であると判定し 
IndexFinal＜50%ならば、その領域は脳出血部位でないと判定し、 
前記後処理工程では、前記分析・判定により脳出血部位であると判定された領域についてパーシャルボリューム効果の影響を補正する、 
ことを特徴とする脳出血部位セグメンテーションする装置の作動方法。
【請求項５】
前記前処理における頭骨境界の特定は、前記第３の設定値より小さい値から大きい値へのＣＴ値変化点または前記第３の設定値より大きい値から小さい値へのＣＴ値変化点を検出することにより行う 
ことを特徴とする請求項４に記載の脳出血部位セグメンテーションする装置の作動方法。
【請求項６】
前記１次探索および前記２次探索における領域探索は、隣接する画素間のＣＴ値の相違が５以下である領域を探索することによって行う 
ことを特徴とする請求項４に記載の脳出血部位セグメンテーションする装置の作動方法。
【請求項７】
前記後処理における前記補正は、前記領域に関するダイレーション演算によって行う
ことを特徴とする請求項４に記載の脳出血部位セグメンテーションする装置の作動方法。
【請求項８】
前記第１の設定値は２４５であり、前記第２の設定値は３０であり、前記第３の設定値は１９０であり、前記第４の設定値は３００であり、前記第５の設定値は４であり、前記第６の設定値は１００であり、前記第７の設定値は４０であり、前記第８の設定値は７０であり、前記第９の設定値は０．０１５であり、前記第１０の設定値は０．００３であることを特徴とする請求項４に記載の脳出血部位セグメンテーションする装置の作動方法。 
【請求項９】
頭部のＸ線ＣＴ画像において脳出血部位をセグメンテーションする装置であって、
前記頭部のＸ線ＣＴ画像において、頭骨の内側領域を区分する手段と、
該頭骨の内側領域を区分する手段により区分された領域の中から、脳出血部位が含まれる可能性のある領域を区分する手段と、 
該脳出血部位が含まれる可能性のある領域を区分する手段により区分された領域の中から、脳出血部位を判定する手段と
を備え、
前記脳出血部位が含まれる可能性のある領域を区分する手段は、
前記頭骨の内側領域を区分する手段により区分された領域のうち、ＣＴ値が緩やかに変化する領域を区分する手段と、
該ＣＴ値が緩やかに変化する領域を区分する手段により区分された領域のうち、画素数が所定の数よりも大きい領域を区分する手段と、 
該画素数が所定の数よりも大きい領域を区分する手段により区分された領域から、比較的ＣＴ値の大きい領域を区分する手段と、 
該比較的ＣＴ値の大きい領域を区分する手段により区分された領域について、ＣＴ値が緩やかに変化する領域を区分する手段と、 
を含んでいる脳出血部位をセグメンテーションする装置。
【請求項１０】
前記脳出血部位を判定する手段は、
前記脳出血部位が含まれる可能性のある領域を区分する手段により区分された領域について、該領域のＣＴ値及び大きさに基づいて、該領域が脳出血部位か否かを判定する第１の手段と、 
該第１の手段により脳出血部位と判定された領域以外の領域のうち、脳出血部位の可能性がある領域について、該領域とその周辺の領域とのＣＴ値の差又は該領域のＣＴ値に基づいて、該領域が脳出血部位か否かを判定する第２の手段と、 
を含んでいる請求項９に記載の脳出血部位をセグメンテーションする装置。
【請求項１１】
前記脳出血部位を判定する手段により脳出血部位であると判定された領域についてバーシャルボリューム効果の影響を補正する手段をさらに含んでいる請求項９に記載の脳出血部位をセグメンテーションする装置。 
頭部のＸ線ＣＴ画像について脳出血部位をセグメンテーションする装置であって、
前処理手段と、１次探索手段と、分析・調整手段と、フィルタリング手段と、２次探索手段と、分析・判定手段と、後処理手段とを有し、 
前記前処理手段は、頭部のＸ線ＣＴ画像について、
第１の設定値より大きいＣＴ値を持つ画素と第２の設定値より小さいＣＴ値を持つ画素を除外し、
第３の設定値に基づいて頭骨の境界を特定し、
前記特定した境界に基づいて頭骨の外側領域を除外し、
前記１次探索手段は、前記前処理が行われた画像について、
ＣＴ値変化が緩やかな領域を探索し、
前記探索により判明した領域ごとにラベルを付与し、
前記分析・調整手段は、前記１次探索が行われた画像について、
前記領域ごとに画素数を求め、
画素数が第４の設定値より少ない領域を除外し、
前記フィルタリング手段は、前記分析・調整が行われた画像について、
前記領域ごとに、
２次元座標ｉ，ｊが増加する方向にそれぞれ隣接する画素間のＣＴ値の差の絶対値の和を画素ごとに求め、
前記和が第５の設定値以上となる画素位置を特定し、
前記特定した全ての画素位置におけるＣＴ値の平均値を求め、
前記平均値以下のＣＴ値を持つ画素を除外し、
前記２次探索手段は、前記フィルタリングが行われた画像について、
ＣＴ値変化が緩やかな領域を探索し、
前記探索により判明した領域ごとにあらためてラベルを付与し、
前記分析・判定手段は、
前記２次探索においてラベルが付与された画像について、
前記領域の画素のＣＴ値をCTregion、第６の設定値をCTmax、第７の設定値をCTmin、第８の設定値をCTbloodとしたとき、 
第１のインデックスIndexCTを、 
CTmin＜= CTregion＜= CTbloodならば
【数１０】
0005383486.tif 000039
とし、
CTblood＜ CTregion＜= CTmaxならば
【数１１】
0005383486.tif 000040
とし、
それ以外は
【数１２】
0005383486.tif 000041
とし、
前記領域の面積および周長をそれぞれAreaおよびPerimeterで表し、前記領域の特性値を 
【数１３】
0005383486.tif 000042
としたとき、
第２のインデックスindexradianを、 
Radianが第９の設定値Radianmaxより大きければindexradian =1とし、
Radianが第１０の設定値Radianminより小さければindexradian =0とし、
Radianが第９の設定値Radianmax以下で第１０の設定値Radianmin以上ならば 
【数１４】
0005383486.tif 000043
とし、
第３のインデックスIndex0を、
【数１５】
0005383486.tif 000044
とし、
Index0＞=20%ならば、その領域は脳出血部位であると判定し、 
Index0＜=3%ならば、その領域は脳出血部位でないと判定し、 
3%＜Index0＜20%ならば、そのような領域の周囲の領域の画素のＣＴ値をCTAroundRegionとし、そのような全ての領域の平均ＣＴ値をCTAverage -All-Regionとしたとき、
第４のインデックスIndexSubを、 
【数１６】
0005383486.tif 000045
とし、
第５のインデックスIndexOderを、 
【数１７】
0005383486.tif 000046
とし、
IndexSub＞=8またはIndexOrder＞ =10ならば、その領域は脳出血部位であると判定し、
IndexSub＜=0またはIndexOrder＜ =-5ならば、その領域は脳出血部位でないと判定し、
0＜IndexSub＜ 8かつ-5＜ IndexOrder＜10ならば、 
第６のインデックスIndexFinalを、 
【数１８】
0005383486.tif 000047
とし、
IndexFinal＞=50%ならば、その領域は脳出血部位であると判定し 
IndexFinal＜50%ならば、その領域は脳出血部位でないと判定し、 
前記後処理手段は、前記分析・判定により脳出血部位であると判定された領域についてパーシャルボリューム効果の影響を補正する、 
ことを特徴とする脳出血部位セグメンテーション装置。
【請求項１２】
前記前処理における頭骨境界の特定は、前記第３の設定値より小さい値から大きい値へのＣＴ値変化点または前記第３の設定値より大きい値から小さい値へのＣＴ値変化点を検出することにより行う 
ことを特徴とする請求項１１に記載の脳出血部位セグメンテーション装置。
【請求項１３】
前記１次探索および前記２次探索における領域探索は、隣接する画素間のＣＴ値の相違が５以下である領域を探索することによって行う 
ことを特徴とする請求項１１に記載の脳出血部位セグメンテーション装置。
【請求項１４】
前記後処理における前記補正は、前記領域に関するダイレーション演算によって行う
ことを特徴とする請求項１１に記載の脳出血部位セグメンテーション装置。
【請求項１５】
前記第１の設定値は２４５であり、前記第２の設定値は３０であり、前記第３の設定値は１９０であり、前記第４の設定値は３００であり、前記第５の設定値は４であり、前記第６の設定値は１００であり、前記第７の設定値は４０であり、前記第８の設定値は７０であり、前記第９の設定値は０．０１５であり、前記第１０の設定値は０．００３であることを特徴とする請求項１１に記載の脳出血部位セグメンテーション装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、脳出血部位セグメンテーション(cerebral hemorrhage segmentation)方法および装置に関し、特に、脳出血を発症している頭部のＸ線ＣＴ(Computed Tomography)画像について、脳出血部位を特定する方法および装置に関する。
【背景技術】
【０００２】
脳出血(cerebral hemorrhage)の診断および治療のために、頭部のＸ線ＣＴ画像について、脳出血部位のセグメンテーションが行われる。セグメンテーションは専門医が介在して人為的に行われる（例えば、特許文献１参照）。
【特許文献１】特開２００５−１１８５１０号公報
【発明の開示】
【発明が解決しようとする課題】
【０００３】
専門医が介在する人為的なセグメンテーションは、労力と時間がかかり、しかもセグメンテーションの結果は担当者の技量等に左右される。また、自動化するにしても、脳出血部位のＣＴ値は症状によって様々な値を取り、正常な他の部位のＣＴ値と重なることもあるので、単一の閾値によるセグメンテーションでは目的が達せられない。
【０００４】
そこで、本発明の目的は、脳出血部位のセグメンテーションを適切に行う方法および装置を実現することである。
【課題を解決するための手段】
【０００５】
上記の課題を解決するための第１の観点での発明は、
頭部の医用画像において脳出血部位をセグメンテーションする方法であって、
前記頭部の医用画像において、頭骨の内側領域を区分する工程と、
該頭骨の内側領域を区分する工程において区分された領域の中から、脳出血部位が含まれる可能性のある領域を区分する工程と、
該脳出血部位が含まれる可能性のある領域を区分する工程により区分された領域の中から、脳出血部位を判定する工程と
を備えた脳出血部位をセグメンテーションする方法である。
【０００６】
さらに好ましくは、前記脳出血部位が含まれる可能性のある領域を区分する工程は、
前記頭骨の内側領域を区分する工程において区分された領域のうち、ＣＴ値が緩やかに変化する領域を区分する工程と、
該ＣＴ値が緩やかに変化する領域を区分する工程において区分された領域のうち、画素数が所定の数よりも大きい領域を区分する工程と、
該画素数が所定の数よりも大きい領域を区分する工程において区分された領域から、比較的ＣＴ値の大きい領域を区分する工程と、
該比較的ＣＴ値の大きい領域を区分する工程において抽出された領域について、ＣＴ値が緩やかに変化する領域を区分する工程と、
を含んでいる。
【０００７】
さらに好ましくは、
前記脳出血部位を判定する工程は、
前記脳出血部位が含まれる可能性のある領域を区分する工程において区分された領域について、該領域のＣＴ値及び大きさに基づいて、該領域が脳出血部位か否かを判定する第１の工程と、
該第１の工程において脳出血部位と判定された領域以外の領域のうち、脳出血部位の可能性がある領域について、該領域とその周辺の領域とのＣＴ値の差又は該領域のＣＴ値に基づいて、該領域が脳出血部位か否かを判定する第２の工程と、
を含んでいる。
【０００８】
さらに好ましくは、
脳出血部位をセグメンテーションする方法は、
前記脳出血部位を判定する工程において脳出血部位であると判定された領域についてバーシャルボリューム効果の影響を補正する工程をさらに含んでいる。
【０００９】
さらに特定的には、本発明は、頭部のＸ線ＣＴ画像について脳出血部位をセグメンテーションする方法であって、前処理工程と、１次探索工程と、分析・調整工程と、フィルタリング工程と、２次探索工程と、分析・判定工程と、後処理工程とを有し、
前記前処理工程では、頭部のＸ線ＣＴ画像について、第１の設定値より大きいＣＴ値を持つ画素と第２の設定値より小さいＣＴ値を持つ画素を除外し、第３の設定値に基づいて頭骨の境界を特定し、前記特定した境界に基づいて頭骨の外側領域を除外し、
前記１次探索工程では、前記前処理が行われた画像について、ＣＴ値変化が緩やかな領域を探索し、前記探索により判明した領域ごとにラベルを付与し、
前記分析・調整工程では、前記１次探索が行われた画像について、前記領域ごとに画素数を求め、画素数が第４の設定値より少ない領域を除外し、
前記フィルタリング工程では、前記分析・調整が行われた画像について、前記領域ごとに、２次元座標ｉ，ｊが増加する方向にそれぞれ隣接する画素間のＣＴ値の差の絶対値の和を画素ごとに求め、前記和が第５の設定値以上となる画素位置を特定し、前記特定した全ての画素位置におけるＣＴ値の平均値を求め、前記平均値以下のＣＴ値を持つ画素を除外し、
前記２次探索工程では、前記フィルタリングが行われた画像について、ＣＴ値変化が緩やかな領域を探索し、前記探索により判明した領域ごとにあらためてラベルを付与し、
前記分析・判定工程では、前記２次探索においてラベルが付与された画像について、前記領域の画素のＣＴ値をCTregion、第６の設定値をCTmax、第７の設定値をCTmin、第８の設定値をCTbloodとしたとき、
第１のインデックスIndexCTを、
CTmin＜=CTregion＜=CTbloodならば
【００１０】
【数１９】
0005383486.tif 000002
【００１１】
とし、
CTblood＜CTregion＜=CTmaxならば
【００１２】
【数２０】
0005383486.tif 000003
【００１３】
とし、
それ以外は
【００１４】
【数２１】
0005383486.tif 000004
【００１５】
とし、
前記領域の面積および周長をそれぞれAreaおよびPerimeterで表し、前記領域の特性値を
【００１６】
【数２２】
0005383486.tif 000005
【００１７】
としたとき、
第２のインデックスindexradianを、
Radianが第９の設定値Radianmaxより大きければindexradian=1とし、
Radianが第１０の設定値Radianminより小さければindexradian=0とし、
Radianが第９の設定値Radianmax以下で第１０の設定値Radianmin以上ならば
【００１８】
【数２３】
0005383486.tif 000006
【００１９】
とし、
第３のインデックスIndex0を、
【００２０】
【数２４】
0005383486.tif 000007
【００２１】
とし、
Index0＞=20%ならば、その領域は脳出血部位であると判定し、
Index0＜=3%ならば、その領域は脳出血部位でないと判定し、
3%＜Index0＜20%ならば、そのような領域の周囲の領域の画素のＣＴ値をCTAroundRegionとし、そのような全ての領域の平均ＣＴ値をCTAverage-All-Regionとしたとき、
第４のインデックスIndexSubを、
【００２２】
【数２５】
0005383486.tif 000008
【００２３】
とし、
第５のインデックスIndexOderを、
【００２４】
【数２６】
0005383486.tif 000009
【００２５】
とし、
IndexSub＞=8またはIndexOrder＞=10ならば、その領域は脳出血部位であると判定し、
IndexSub＜=0またはIndexOrder＜=-5ならば、その領域は脳出血部位でないと判定し、
0＜IndexSub＜8かつ-5＜IndexOrder＜10ならば、
第６のインデックスIndexFinalを、
【００２６】
【数２７】
0005383486.tif 000010
【００２７】
とし、
IndexFinal＞=50%ならば、その領域は脳出血部位であると判定し
IndexFinal＜50%ならば、その領域は脳出血部位でないと判定し、
前記後処理工程では、前記分析・判定により脳出血部位であると判定された領域についてパーシャルボリューム効果の影響を補正する、
ことを特徴とする脳出血部位セグメンテーション方法である。
【００２８】
上記の課題を解決するための第２の観点での発明は、
頭部のＸ線ＣＴ画像において脳出血部位をセグメンテーションする装置であって、
前記頭部のＸ線ＣＴ画像において、頭骨の内側領域を区分する手段と、
該頭骨の内側領域を区分する手段により区分された領域の中から、脳出血部位が含まれる可能性のある領域を区分する手段と、
該脳出血部位が含まれる可能性のある領域を区分する手段により区分された領域の中から、脳出血部位を判定する手段と
を備えた脳出血部位をセグメンテーションする装置である。
【００２９】
さらに好ましくは、前記脳出血部位が含まれる可能性のある領域を区分する手段は、
前記頭骨の内側領域を区分する手段により区分された領域のうち、ＣＴ値が緩やかに変化する領域を区分する手段と、
該ＣＴ値が緩やかに変化する領域を区分する手段により区分された領域のうち、画素数が所定の数よりも大きい領域を区分する手段と、
該画素数が所定の数よりも大きい領域を区分する手段により区分された領域から、比較的ＣＴ値の大きい領域を区分する手段と、
該比較的ＣＴ値の大きい領域を区分する手段により区分された領域について、ＣＴ値が緩やかに変化する領域を区分する手段と、
を含んでいる。
【００３０】
さらに好ましくは、前記脳出血部位を判定する手段は、
前記脳出血部位が含まれる可能性のある領域を区分する手段により区分された領域について、該領域のＣＴ値及び大きさに基づいて、該領域が脳出血部位か否かを判定する第１の手段と、
該第１の手段により脳出血部位と判定された領域以外の領域のうち、脳出血部位の可能性がある領域について、該領域とその周辺の領域とのＣＴ値の差又は該領域のＣＴ値に基づいて、該領域が脳出血部位か否かを判定する第２の手段と、
を含んでいる。
【００３１】
さらに好ましくは、脳出血部位をセグメンテーションする装置は、
前記脳出血部位を判定する手段により脳出血部位であると判定された領域についてバーシャルボリューム効果の影響を補正する手段をさらに含んでいる。
【００３２】
さらに特定的には、本発明は、頭部のＸ線ＣＴ画像について脳出血部位をセグメンテーションする装置であって、前処理手段と、１次探索手段と、分析・調整手段と、フィルタリング手段と、２次探索手段と、分析・判定手段と、後処理手段とを有し、
前記前処理手段は、頭部のＸ線ＣＴ画像について、第１の設定値より大きいＣＴ値を持つ画素と第２の設定値より小さいＣＴ値を持つ画素を除外し、第３の設定値に基づいて頭骨の境界を特定し、前記特定した境界に基づいて頭骨の外側領域を除外し、
前記１次探索手段は、前記前処理が行われた画像について、ＣＴ値変化が緩やかな領域を探索し、前記探索により判明した領域ごとにラベルを付与し、
前記分析・調整手段は、前記１次探索が行われた画像について、前記領域ごとに画素数
を求め、画素数が第４の設定値より少ない領域を除外し、
前記フィルタリング手段は、前記分析・調整が行われた画像について、前記領域ごとに、２次元座標ｉ，ｊが増加する方向にそれぞれ隣接する画素間のＣＴ値の差の絶対値の和を画素ごとに求め、前記和が第５の設定値以上となる画素位置を特定し、前記特定した全ての画素位置におけるＣＴ値の平均値を求め、前記平均値以下のＣＴ値を持つ画素を除外し、
前記２次探索手段は、前記フィルタリングが行われた画像について、ＣＴ値変化が緩やかな領域を探索し、前記探索により判明した領域ごとにあらためてラベルを付与し、
前記分析・判定手段は、前記２次探索においてラベルが付与された画像について、前記領域の画素のＣＴ値をCTregion、第６の設定値をCTmax、第７の設定値をCTmin、第８の設定値をCTbloodとしたとき、
第１のインデックスIndexCTを、
CTmin＜=CTregion＜=CTbloodならば
【００３３】
【数２８】
0005383486.tif 000011
【００３４】
とし、
CTblood＜CTregion＜=CTmaxならば
【００３５】
【数２９】
0005383486.tif 000012
【００３６】
とし、
それ以外は
【００３７】
【数３０】
0005383486.tif 000013
【００３８】
とし、
前記領域の面積および周長をそれぞれAreaおよびPerimeterで表し、前記領域の特性値を
【００３９】
【数３１】
0005383486.tif 000014
【００４０】
としたとき、
第２のインデックスindexradianを、
Radianが第９の設定値Radianmaxより大きければindexradian=1とし、
Radianが第１０の設定値Radianminより小さければindexradian=0とし、
Radianが第９の設定値Radianmax以下で第１０の設定値Radianmin以上ならば
【００４１】
【数３２】
0005383486.tif 000015
【００４２】
とし、
第３のインデックスIndex0を、
【００４３】
【数３３】
0005383486.tif 000016
【００４４】
とし、
Index0＞=20%ならば、その領域は脳出血部位であると判定し、
Index0＜=3%ならば、その領域は脳出血部位でないと判定し、
3%＜Index0＜20%ならば、そのような領域の周囲の領域の画素のＣＴ値をCTAroundRegionとし、そのような全ての領域の平均ＣＴ値をCTAverage-All-Regionとしたとき、
第４のインデックスIndexSubを、
【００４５】
【数３４】
0005383486.tif 000017
【００４６】
とし、
第５のインデックスIndexOderを、
【００４７】
【数３５】
0005383486.tif 000018
【００４８】
とし、
IndexSub＞=8またはIndexOrder＞=10ならば、その領域は脳出血部位であると判定し、
IndexSub＜=0またはIndexOrder＜=-5ならば、その領域は脳出血部位でないと判定し、
0＜IndexSub＜8かつ-5＜IndexOrder＜10ならば、
第６のインデックスIndexFinalを、
【００４９】
【数３６】
0005383486.tif 000019
【００５０】
とし、
IndexFinal＞=50%ならば、その領域は脳出血部位であると判定し
IndexFinal＜50%ならば、その領域は脳出血部位でないと判定し、
前記後処理手段は、前記分析・判定により脳出血部位であると判定された領域についてパーシャルボリューム効果の影響を補正する、
ことを特徴とする脳出血部位セグメンテーション装置である。
【００５１】
前記前処理における頭骨境界の特定は、前記第３の設定値より小さい値から大きい値へのＣＴ値変化点または前記第３の設定値より大きい値から小さい値へのＣＴ値変化点を検出することにより行うことが、境界特定を適切に行う点で好ましい。
【００５２】
前記１次探索および前記２次探索における領域探索は、隣接する画素間のＣＴ値の相違が５以下である領域を探索することによって行うことが、領域探索を適切に行う点で好ましい。
【００５３】
前記後処理における前記補正は、前記領域に関するダイレーション演算によって行うことが、領域補正を適切に行う点で好ましい。
前記第１の設定値は２４５であり、前記第２の設定値は３０であり、前記第３の設定値は１９０であり、前記第４の設定値は３００であり、前記第５の設定値は４であり、前記第６の設定値は１００であり、前記第７の設定値は４０であり、前記第８の設定値は７０であり、前記第９の設定値は０．０１５であり、前記第１０の設定値は０．００３であることが、セグメンテーションを高精度に行う点で好ましい。
【発明の効果】
【００５４】
上記各観点での本発明によれば、候補領域の画素のＣＴ値に関するインデックスIndexCTと、候補領域の特性に関するインデックスindexradianを用いてインデックスIndex0を生成し、Index0＞=20%ならばその領域は脳出血部位であると判定し、Index0＜=3%ならばその領域は脳出血部位でないと判定する。3%＜Index0＜20%のときは、そのような領域についてインデックスIndexSubとインデックスIndexOderを生成し、IndexSub＞=8またはIndexOrder＞=10ならばその領域は脳出血部位であると判定し、IndexSub＜=0またはIndexOrder＜=-5ならばその領域は脳出血部位でないと判定する。0＜IndexSub＜8かつ-5＜IndexOrder＜10のときは、インデックスIndexFinalを生成し、IndexFinal＞=50%ならばその領域は脳出血部位であると判定しIndexFinal＜50%ならばその領域は脳出血部位でないと判定する。これにより、脳出血部位のセグメンテーションを適切に行う方法および装置を実現することができる。
【発明を実施するための最良の形態】
【００５５】
以下、図面を参照して発明を実施するための最良の形態を説明する。なお、本発明は、発明を実施するための最良の形態に限定されるものではない。図１に画像処理装置の構成をブロック(block)図によって示す。
【００５６】
本装置は発明を実施するための最良の形態の一例である。本装置の構成によって、脳出血部位セグメンテーション装置に関する発明を実施するための最良の形態の一例が示される。本装置の動作によって、脳出血部位セグメンテーション方法に関する発明を実施するための最良の形態の一例が示される。
【００５７】
図１に示すように、本装置は、データ(data)処理部１０、表示部２０、操作部３０、記憶部４０および入出力部５０を有する。データ処理部１０は、表示部２０および操作部３０を通じてのユーザー(user)によるインタラクティブ(interactive)な操作に基づいて、記憶部４０に記憶されたデータについて所定のデータ処理を行う。
【００５８】
データ処理部１０は、また、入出力部５０を通じて外部装置に対するデータ入出力を行う。脳出血部位セグメンテーションを行うべきＸ線ＣＴ画像も、入出力部５０を通じて外部装置から入力される。
【００５９】
外部装置は例えばＸ線ＣＴ装置や医用画像サーバー(server)等である。なお、本装置は、Ｘ線ＣＴ装置や医用画像サーバーの一部であってもよい。その場合は、入出力部５０は必ずしも必要でない。図２に、脳出血部位セグメンテーション対象となるＸ線ＣＴ画像の一例を示す。
【００６０】
本装置によって行われる脳出血部位セグメンテーションについて説明する。図３に、脳出血部位セグメンテーションの工程を示す。図３に示すように、脳出血部位セグメンテーションは７つの工程Ｐ１−Ｐ７によって行われる。
【００６１】
工程Ｐ１は前処理工程である。工程Ｐ２は１次探索工程である。工程Ｐ３は分析・調整工程である。工程Ｐ４はフィルタリング工程である。工程Ｐ５は２次探索工程である。工程Ｐ６は分析・判定工程である。工程Ｐ７は後処理工程である。
【００６２】
これらの工程は、データ処理部１０によって遂行される。データ処理部１０は、本発明における前処理手段の一例であり、１次探索手段の一例であり、分析・調整手段の一例であり、フィルタリング手段の一例であり、２次探索手段の一例であり、分析・判定手段の一例であり、後処理手段の一例である。以下、各工程について詳細に説明する。
【００６３】
図４に、前処理工程Ｐ１の詳細なステップ(step)を示す。図４示すように、前処理工程Ｐ１では、頭部のＸ線ＣＴ画像について、ステップ１０１で、画素選別を行う。画素選別は、第１の設定値より大きいＣＴ値を持つ画素と、第２の設定値より小さいＣＴ値を持つ画素を、画像全体にわたって除外することによって行われる。第１の設定値は例えば２４５であり、第２の設定値は例えば３０である。これによって、２４５より大きいＣＴ値を持つ画素と、３０より小さいＣＴ値を持つ画素が除外される。
【００６４】
ステップ１０２で、頭骨境界特定が行われる。頭骨境界特定は、第３の設定値に基づいて行われる。頭骨境界の特定は、第３の設定値より小さい値から大きい値へのＣＴ値変化点または第３の設定値より大きい値から小さい値へのＣＴ値変化点を、画像全体にわたって検出することにより行われる。第３の設定値は例えば１９０である。特定された頭骨境界に基づいて、ステップ１０３で、頭骨の外側領域を除外することにより、頭骨の内側領域を区分する。
【００６５】
ステップ１０１−１０３の処理によって、２４５より大きいＣＴ値を持つ画素と３０より小さいＣＴ値を持つ画素、および、頭骨の外側領域が除外される。このような前処理によって、例えば図５に示すような画像が得られる。図５において、黒塗り部分は、除外されずに残った部分を示す。この部分が、次工程の処理の対象となる。
【００６６】
図６に、１次探索工程Ｐ２の詳細なステップを示す。図６示すように、１次探索工程Ｐ２では、前処理が行われた画像について、ステップ２０１で、ＣＴ値変化が緩やかな領域を探索する。ＣＴ値変化が緩やかな領域の探索は、隣接する画素間のＣＴ値の相違が例えば５以下である領域を探索することによって行われる。そして、ステップ２０２で、探索により判明した領域ごとにラベルを付与することにより、ＣＴ値が緩やかに変化する領域を区分する。
【００６７】
図７に、分析・調整工程Ｐ３の詳細なステップを示す。図７示すように、分析・調整工程Ｐ３では、１次探索が行われた画像について、ステップ３０１で、領域ごとに画素数を求める。そして、ステップ３０２で、画素数が第４の設定値より少ない領域を除外することにより、画素数が第４の設定値以上の領域を区分する。
【００６８】
第４の設定値は、画像のデファインドＦＯＶ(defined Field of View)が２５ｃｍのとき例えば３００である。なお、この設定値はデファインドＦＯＶが２５ｃｍ以外のときは３００以外の適宜の値としてよい。
【００６９】
このような１次探索および分析・調整によって、例えば図８に示すような画像が得られる。図８において、黒塗り部分は、ＣＴ値変化が緩やかで画素数が３００以上の部分を示す。この部分が、次工程の処理の対象となる。
【００７０】
図９に、フィルタリング工程Ｐ４の詳細なステップを示す。図９示すように、フィルタリング工程Ｐ４では、分析・調整が行われた画像について、ステップ４０１で、領域ごとに、２次元座標ｉ，ｊが増加する方向にそれぞれ隣接する画素間のＣＴ値の差の絶対値の和を画素ごとに求める。すなわち、次式による演算を行う。
【００７１】
【数３７】
0005383486.tif 000020
【００７２】
ここで、F[i,j]は２次元座標ｉ，ｊにおける画素のＣＴ値であり、F[i+1,j]は座標ｉが増加する方向に隣接する画素のＣＴ値であり、F[i,j+1]は座標jが増加する方向に隣接する画素のＣＴ値である。ＣＴ値F[i,j]、F[i+1,j]、 F[i,j+1]およびF[i+1,j+1]の空間的位置関係は図１０に示すとおりである。
【００７３】
ステップ４０２で、G[i,j]が第５の設定値以上となる画素位置を特定する。第５の設定値は例えば４である。ステップ４０３で、特定された全ての画素位置におけるＣＴ値の平均値を求め、ステップ４０４で、平均値以下のＣＴ値を持つ画素を除外することにより、残りの画素を区分する。
【００７４】
ＣＴ値の平均値は、元画像について、領域ごとに求められる。また、平均値以下のＣＴ値を持つ画素の除外は、元画像について、領域ごとに行われる。これによって、バックグラウンド(back ground)の除去が領域ごとに行われる。
【００７５】
このようなフィルタリングによって、例えば図１１に示すような画像が得られる。図１１において、黒塗り部分は、除去されなかった部分を示す。この部分が、次工程の処理の対象となる。
【００７６】
図１２に、２次探索工程Ｐ５の詳細なステップを示す。図１２示すように、２次探索工程Ｐ５では、フィルタリングが行われた画像について、ステップ５０１で、ＣＴ値変化が緩やかな領域を探索する。ＣＴ値変化が緩やかな領域の探索は、隣接する画素間のＣＴ値の相違が例えば５以下である領域を探索することによって行われる。そして、ステップ５０２で、探索により判明した領域ごとにあらためてラベルを付与することによって、ＣＴ値が緩やかに変化する領域を区分する。このような２次探索によって、例えば図１３に示すような画像が得られる。
【００７７】
図１４に、分析・判定工程Ｐ６の詳細なステップを示す。図１４示すように、分析・判定工程Ｐ６では、２次探索によりラベルが付与された画像について、ステップ６０１で、第１のインデックスIndexCTを求める。
【００７８】
インデックスIndexCTは、領域の画素のＣＴ値をCTregion、第６の設定値をCTmax、第７の設定値をCTmin、第８の設定値をCTbloodとしたとき、
CTmin＜=CTregion＜=CTbloodならば
【００７９】
【数３８】
0005383486.tif 000021
【００８０】
とし、
CTblood＜CTregion＜=CTmaxならば
【００８１】
【数３９】
0005383486.tif 000022
【００８２】
とし、
それ以外は
【００８３】
【数４０】
0005383486.tif 000023
【００８４】
とする。なお、第６の設定値CTmaxは例えば１００であり、第７の設定値CTminは例えば４０であり、第８の設定値CTbloodは７０である。
ステップ６０２で、第２のインデックスindexradianを求める。第２のインデックスindexradianは、領域の面積および周長をそれぞれAreaおよびPerimeterで表し、領域の特性値を
【００８５】
【数４１】
0005383486.tif 000024
【００８６】
としたとき、
Radianが第９の設定値Radianmaxより大きければindexradian=1とし、Radianが第１０の設定値Radianminより小さければindexradian=0とし、Radianが第９の設定値Radianmax以下で第１０の設定値Radianmin以上ならば
【００８７】
【数４２】
0005383486.tif 000025
【００８８】
とする。なお、第９の設定値Radianmaxは例えば０．０１５であり、第１０の設定値Radianminはたとえばあ０．００３である。
ステップ６０３で、第３のインデックスIndex0を、
【００８９】
【数４３】
0005383486.tif 000026
【００９０】
とする。
ステップ６０４で、第３のインデックスIndex0の値に応じて、その領域が脳出血部位であるか否かを判定する。すなわち、Index0＞=20%ならば、その領域は脳出血部位である(CH)と判定し、Index0＜=3%ならば、その領域は脳出血部位でない(not CH)と判定し、3%＜Index0＜20%ならば、次のステップ６０５に移行する。
【００９１】
ステップ６０５では、第４のインデックスIndexSubと第５のインデックスIndexOderを求める。すなわち、3%＜Index0＜20%である領域の周囲の領域の画素のＣＴ値をCTAroundRegionとし、3%＜Index0＜20%である全ての領域の平均ＣＴ値をCTAverage-All-Regionとしたとき、
第４のインデックスIndexSubを、
【００９２】
【数４４】
0005383486.tif 000027
【００９３】
とし、
第５のインデックスIndexOderを、
【００９４】
【数４５】
0005383486.tif 000028
【００９５】
とする。
ステップ６０６で、第４のインデックスIndexSubと第５のインデックスIndexOderの値に応じて、その領域が脳出血部位であるか否かを判定する。すなわち、IndexSub＞=8またはIndexOrder＞=10ならば、その領域は脳出血部位である(CH)と判定し、IndexSub＜=0またはIndexOrder＜=-5ならば、その領域は脳出血部位でない(not CH)と判定し、0＜IndexSub＜8かつ-5＜IndexOrder＜10ならば、次のステップ６０７移行する。
【００９６】
ステップ６０７では、第６のインデックスIndexFinalを、
【００９７】
【数４６】
0005383486.tif 000029
【００９８】
とする。
ステップ６０８で、第６のインデックスIndexFinalの値に応じて、その領域が脳出血部位であるか否かを判定する。すなわち、IndexFinal＞=50%ならば、その領域は脳出血部位
である(CH)と判定し、IndexFinal＜50%ならば、その領域は脳出血部位でない(not CH)と判定する。
【００９９】
以上のような３段階の分析・判定によって、脳出血部位の正確なセグメンテーションが行われ、例えば図１５に示すような画像が得られる。図１５において、黒塗り部分は脳出血部位を示す。黒塗り部分は、実際の脳出血部位の位置と範囲を精度良く示すものとなる。
【０１００】
脳出血部位が頭骨に近接しているときは、パーシャルボリューム効果(partial volume effect)のために脳出血部位のＣＴ値が変化して、脳出血部位のセグメンテーションが不正確になることがあり得る。
【０１０１】
図１６に、そのような例を示す。図１６において、黒塗り部分がセグメンテーションされた領域であるが、この領域の右側の輪郭は頭骨境界に達していない。これは、パーシャルボリューム効果でＣＴ値が変化した部分が除外されたためである。このようなセグメンテーション結果については、後処理工程Ｐ７によって、パーシャルボリューム効果の影響を補正する。
【０１０２】
図１７に、後処理工程Ｐ７の詳細なステップを示す。図１７示すように、後処理工程Ｐ７では、分析・判定により脳出血部位であると判定された領域すなわちＣＨ領域について、ステップ７０１で、周囲の画素のＣＴ値がＣＨ領域の平均ＣＴ値より大きいか否かを判定し、大きいときは、ＣＨ領域の外縁が頭骨境界に到達するまでダイレーション(dilation)演算を行う。このような後処理によって、例えば図１８に示すように、ＣＨ領域についてパーシャルボリューム効果の影響が補正された画像が得られる。
【図面の簡単な説明】
【０１０３】
【図１】本発明を実施するための最良の形態の一例の画像処理装置の構成を示す図である。
【図２】Ｘ線ＣＴ画像の一例を中間調の写真で示す図である。
【図３】脳出血部位セグメンテーションを行うための工程を示す図である。
【図４】前処理工程の詳細を示す図である。
【図５】前処理後のＸ線ＣＴ画像の一例を中間調の写真で示す図である。
【図６】１次探索工程の詳細を示す図である。
【図７】分析・調査工程の詳細を示す図である。
【図８】１次探索および分析・調査後のＸ線ＣＴ画像の一例を中間調の写真で示す図である。
【図９】フィルタリング工程の詳細を示す図である。
【図１０】画素値の空間的位置関係を示す図である。
【図１１】フィルタリング後のＸ線ＣＴ画像の一例を中間調の写真で示す図である。
【図１２】２次探索工程の詳細を示す図である。
【図１３】２次探索後のＸ線ＣＴ画像の一例を中間調の写真で示す図である。
【図１４】分析・判定工程の詳細を示す図である。
【図１５】分析・判定後のＸ線ＣＴ画像の一例を中間調の写真で示す図である。
【図１６】分析・判定後のＸ線ＣＴ画像の一例を中間調の写真で示す図である。
【図１７】後処理工程の詳細を示す図である。
【図１８】後処理後のＸ線ＣＴ画像の一例を中間調の写真で示す図である。
【符号の説明】
【０１０４】
１０  データ処理部
２０  表示部
３０  操作部
４０  記憶部
５０  入出力部
0005383486.tif 000048
0005383486.tif 000049
0005383486.tif 000050
0005383486.tif 000051
0005383486.tif 000052
0005383486.tif 000053
0005383486.tif 000054
0005383486.tif 000055
0005383486.tif 000056
0005383486.tif 000057
0005383486.tif 000058
0005383486.tif 000059
0005383486.tif 000060
0005383486.tif 000061
0005383486.tif 000062
0005383486.tif 000063
0005383486.tif 000064
【図１】
0005383486.tif 000065
【図３】
0005383486.tif 000066
【図４】
0005383486.tif 000067
【図６】
0005383486.tif 000068
【図７】
0005383486.tif 000069
【図９】
0005383486.tif 000070
【図１０】
0005383486.tif 000071
【図１２】
0005383486.tif 000072
【図１４】
0005383486.tif 000073
【図１７】
0005383486.tif 000074
【図２】
0005383486.tif 000075
【図５】
0005383486.tif 000076
【図８】
0005383486.tif 000077
【図１１】
0005383486.tif 000078
【図１３】
0005383486.tif 000079
【図１５】
0005383486.tif 000080
【図１６】
0005383486.tif 000081
【図１８】
0005383486.tif 000082