(57)【要約】【課題】肺の低吸収領域（低吸収体積領域）を三次元的に評価する。【解決手段】肺の三次元的ＣＴデータを読み込み（Ｓ１）、肺の三次元的ＣＴデータからボリュームデータを作成し（Ｓ２）、ボリュームデータから肺野体積領域を抽出し（Ｓ３）、肺野体積領域中の低吸収体積領域を抽出し（Ｓ４）、各低吸収体積領域の壁を除去し（Ｓ５）、壁を除去した各低吸収体積領域の体積を基に評価値を算出する（Ｓ６）。また、従来互換の評価値も算出する（Ｓ７）。【効果】肺の低吸収領域（低吸収体積領域）を三次元的に評価できる。【選択図】図２
2004049498.tif 000002
【特許請求の範囲】
【請求項１】
肺の三次元的ＣＴデータからボリュームデータを作成するボリュームデータ作成ステップと、ボリュームデータから肺野体積領域を抽出する肺野体積領域抽出ステップと、肺野体積領域中の低吸収体積領域を抽出する低吸収体積領域抽出ステップと、各低吸収体積領域の体積を基に評価値を算出する評価値算出ステップとをこの順に実行することを特徴とする肺の低吸収領域評価方法。
【請求項２】
請求項１に記載の肺の低吸収領域評価方法において、前記低吸収体積領域抽出ステップと前記評価値算出ステップの間に、各低吸収体積領域の壁部分を各低吸収体積領域から除去する壁除去ステップを実行することを特徴とする肺の低吸収領域評価方法。
【請求項３】
請求項１または請求項２に記載の肺の低吸収領域評価方法において、前記評価値として、各低吸収体積領域の体積の最大値、最小値、平均値、合計値、合計値の肺野体積領域に対する体積比、分布、分散、低吸収体積領域の数の少なくとも一つを算出することを特徴とする肺の低吸収領域評価方法。
【請求項４】
請求項３に記載の肺の低吸収領域評価方法において、前記評価値として、上肺、中肺、下肺における各低吸収領域の面積の最大値、最小値、平均値、合計値、合計値の肺野領域に対する面積比の少なくとも一つを算出することを特徴とする肺の低吸収領域評価方法。
【請求項５】
肺の三次元的ＣＴデータからボリュームデータを作成するボリュームデータ作成手段と、ボリュームデータから肺野体積領域を抽出する肺野体積領域抽出手段と、肺野体積領域中の低吸収体積領域を抽出する低吸収体積領域抽出手段と、各低吸収体積領域の体積を基に評価値を算出する評価値算出手段とを具備することを特徴とする肺の低吸収領域評価装置。
【請求項６】
請求項５に記載の肺の低吸収領域評価装置において、各低吸収体積領域の壁部分を各低吸収体積領域から除去する壁除去手段をさらに具備することを特徴とする肺の低吸収領域評価装置。
【請求項７】
請求項５または請求項６に記載の肺の低吸収領域評価装置において、前記評価値として、各低吸収体積領域の体積の最大値、最小値、平均値、合計値、合計値の肺野体積領域に対する体積比、分布、低吸収体積領域の数の少なくとも一つを算出することを特徴とする肺の低吸収領域評価装置。
【請求項８】
請求項７に記載の肺の低吸収領域評価装置において、前記評価値として、上肺、中肺、下肺における各低吸収領域の面積の最大値、最小値、平均値、合計値、合計値の肺野領域に対する面積比の少なくとも一つを算出することを特徴とする肺の低吸収領域評価装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、肺の低吸収領域評価方法および装置に関し、さらに詳しくは、肺の低吸収領域（Ｌｏｗ　Ａｔｔｅｎｕａｔｉｏｎ　Ａｒｅａ：ＬＡＡ）を三次元的に評価できる肺の低吸収領域評価方法および装置に関する。
【０００２】
【従来の技術】
図４は、従来行われている肺の低吸収領域計測処理を示すフロー図である。
ステップＳ５１では、図５に示すように、Ｘ線ＣＴ装置より上肺、中肺、下肺の各スライス位置でアキシャルスキャンを行って収集した上肺のＣＴデータＤｔ、中肺のＣＴデータＤｍ、下肺のＣＴデータＤｂを読み込む。
ステップＳ５２では、上肺のＣＴデータＤｔ、中肺のＣＴデータＤｍ、下肺のＣＴデータＤｂに対して、ＣＴ値が−８００以下の領域を肺野領域Ｌｔ、Ｌｍ、Ｌｂとして抽出する。
ステップＳ５３では、上肺の肺野領域Ｌｔ、中肺の肺野領域Ｌｍ、下肺の肺野領域Ｌｂに対して、ＣＴ値が−９６０以下の領域を低吸収領域として抽出する。低吸収領域は、上肺、中肺、下肺の各スライスでの低吸収体積領域Ｃの断面である。
ステップＳ５４では、低吸収領域の評価値として、上肺、中肺、下肺における各低吸収領域の面積の最大値、最小値、平均値、合計値、合計値の肺野領域に対する面積比を算出し、出力する。
【０００３】
【発明が解決しようとする課題】
従来行われている肺の低吸収領域計測処理では、上肺、中肺、下肺における各低吸収領域の面積の最大値、最小値、平均値、合計値、合計値の肺野領域に対する面積比を計測している。そして、これらの評価値により、低吸収領域の評価を行っている。
しかし、ここで評価される低吸収領域は、例えば肺気腫のような三次元構造の一断面であるため、スライス位置により結果が左右されてしまい、常に正確な評価を行えない問題点がある。
そこで、本発明の目的は、肺の低吸収領域（低吸収体積領域）を三次元的に評価できる肺の低吸収領域評価方法および装置を提供することにある。
【０００４】
【課題を解決するための手段】
第１の観点では、本発明は、肺の三次元的ＣＴデータからボリュームデータを作成するボリュームデータ作成ステップと、ボリュームデータから肺野体積領域を抽出する肺野体積領域抽出ステップと、肺野体積領域中の低吸収体積領域を抽出する低吸収体積領域抽出ステップと、各低吸収体積領域の体積を基に評価値を算出する評価値算出ステップとをこの順に実行することを特徴とする肺の低吸収領域評価方法を提供する。
上記第１の観点による肺の低吸収領域評価方法では、肺のボリュームデータから肺野体積領域中の低吸収体積領域を三次元的に抽出し、各低吸収体積領域の体積を基に評価値を算出する。従って、スライス位置により結果が左右されず、低吸収領域（低吸収体積領域）を常に正確に評価することが出来る。
【０００５】
第２の観点では、本発明は、上記構成の肺の低吸収領域評価方法において、前記低吸収体積領域抽出ステップと前記評価値算出ステップの間に、各低吸収体積領域の壁部分を各低吸収体積領域から除去する壁除去ステップを実行することを特徴とする肺の低吸収領域評価方法を提供する。
例えば肺気腫の場合、多数の大きさの異なる低吸収体積領域が現れ、各低吸収体積領域は膜のような壁をそれぞれ持っている。この壁は、低吸収体積領域の実体ではないため、これを低吸収体積領域に算入すると、正確な評価に支障を来すおそれがある。
そこで、上記第２の観点による肺の低吸収領域評価方法では、各低吸収体積領域の壁部分を除去してから評価値を算出する。これにより、低吸収体積領域の実体を正確に評価することが出来る。
【０００６】
第３の観点では、本発明は、上記構成の肺の低吸収領域評価方法において、前記評価値として、各低吸収体積領域の体積の最大値、最小値、平均値、合計値、合計値の肺野体積領域に対する体積比、分布、分散、低吸収体積領域の数の少なくとも一つを算出することを特徴とする肺の低吸収領域評価方法を提供する。
上記第３の観点による肺の低吸収領域評価方法では、各低吸収体積領域の体積の最大値、最小値、平均値、合計値、合計値の肺野体積領域に対する体積比、分布、分散、低吸収体積領域の数の少なくとも一つを評価値とすることが出来る。
【０００７】
第４の観点では、本発明は、上記構成の肺の低吸収領域評価方法において、前記評価値として、上肺、中肺、下肺における各低吸収領域の面積の最大値、最小値、平均値、合計値、合計値の肺野領域に対する面積比の少なくとも一つを算出することを特徴とする肺の低吸収領域評価方法を提供する。
上記第４の観点による肺の低吸収領域評価方法では、従来と互換性のある低吸収領域の評価値をも提供することが出来る。
【０００８】
第５の観点では、本発明は、肺の三次元的ＣＴデータからボリュームデータを作成するボリュームデータ作成手段と、ボリュームデータから肺野体積領域を抽出する肺野体積領域抽出手段と、肺野体積領域中の低吸収体積領域を抽出する低吸収体積領域抽出手段と、各低吸収体積領域の体積を基に評価値を算出する評価値算出手段とを具備することを特徴とする肺の低吸収領域評価装置を提供する。上記第５の観点による低吸収領域評価装置では、上記第１の観点による肺の低吸収領域評価方法を好適に実施できる。
【０００９】
第６の観点では、本発明は、上記構成の肺の低吸収領域評価装置において、各低吸収体積領域の壁部分を各低吸収体積領域から除去する壁除去手段をさらに具備することを特徴とする肺の低吸収領域評価装置を提供する。
上記第６の観点による低吸収領域評価装置では、上記第２の観点による肺の低吸収領域評価方法を好適に実施できる。
【００１０】
第７の観点では、本発明は、上記構成の肺の低吸収領域評価装置において、前記評価値として、各低吸収体積領域の体積の最大値、最小値、平均値、合計値、合計値の肺野体積領域に対する体積比、分布、低吸収体積領域の数の少なくとも一つを算出することを特徴とする肺の低吸収領域評価装置を提供する。
上記第７の観点による低吸収領域評価装置では、上記第３の観点による肺の低吸収領域評価方法を好適に実施できる。
【００１１】
第８の観点では、本発明は、上記構成の肺の低吸収領域評価装置において、前記評価値として、上肺、中肺、下肺における各低吸収領域の面積の最大値、最小値、平均値、合計値、合計値の肺野領域に対する面積比の少なくとも一つを算出することを特徴とする肺の低吸収領域評価装置を提供する。
上記第８の観点による低吸収領域評価装置では、上記第４の観点による肺の低吸収領域評価方法を好適に実施できる。
【００１２】
【発明の実施の形態】
以下、図に示す本発明の実施形態により本発明をさらに詳細に説明する。なお、これにより本発明が限定されるものではない。
【００１３】
図１は、本発明の一実施形態にかかるＸ線ＣＴ装置を示す構成図である。
このＸ線ＣＴ装置１００は、スキャナ装置１と，処理装置２と，大容量記憶装置３と，表示モニタ４と，入力装置５とを具備して構成されている。
【００１４】
前記スキャナ装置１は、ヘリカルスキャンまたは連続的にスライス位置をずらせたアキシャルスキャンにより、肺の三次元的ＣＴデータを収集する。
【００１５】
前記処理装置２は、ＣＴデータ取得部２ａと，ボリュームデータ作成部２ｂと，肺野体積領域抽出部２ｃと，低吸収体積領域抽出部２ｄと、壁除去部２ｅと，評価値算出部２ｆとを具備している。
【００１６】
図２は、肺の低吸収領域計測処理を示すフロー図である。
ステップＳ１では、ＣＴデータ取得部２ａは、肺の三次元的ＣＴデータを読み込む。
ステップＳ２では、ボリュームデータ作成部２ｂは、肺の三次元的ＣＴデータからボリュームデータを作成する。
【００１７】
ステップＳ３では、肺野体積領域抽出部２ｃは、ボリュームデータに対して、ＣＴ値が−８００以下の領域を、図３に示すごとき肺野体積領域Ｌとして抽出する。
ステップＳ４では、低吸収体積領域抽出部２ｄは、肺野体積領域Ｌに対して、ＣＴ値が−９６０以下の領域を、図３に示す如き低吸収体積領域Ｃとして抽出する。
ステップＳ５では、壁除去部２ｅは、各低吸収体積領域の境界を三次元的なエッジ検出法（例えば３×３×３のＳｏｂｅｌオペレータまたはラプラシアンオペレータ）を用いて検出する。次いで、境界近傍で隣接するボクセルのＣＴ値の変化量が少ない部分（例えば３０以内）を壁Ｗとみなし、当該低吸収体積領域Ｃから除去する。
ステップＳ６では、評価値算出部２ｆは、各低吸収体積領域Ｃ−Ｗ（壁Ｗを除去した低吸収体積領域ＣをＣ−Ｗと表記する。）の体積の最大値、最小値、平均値、合計値、合計値の肺野体積領域Ｌに対する体積比、各低吸収体積領域Ｃ−Ｗの体積の分布（ヒストグラム）、各低吸収体積領域Ｃ−Ｗの体積の分散、低吸収体積領域Ｃの数を算出し、出力する。
【００１８】
ステップＳ７では、評価値算出部２ｆは、上肺、中肺、下肺のスライス位置に相当する断面での各低吸収体積領域Ｃの面積（低吸収領域の面積）を算出する。そして、上肺、中肺、下肺における各低吸収領域の面積の最大値、最小値、平均値、合計値、合計値の肺野領域に対する面積比を算出し、出力する。
【００１９】
上記Ｘ線ＣＴ装置１００によれば、次の効果が得られる。
（１）肺の低吸収領域（低吸収体積領域）を三次元的に評価することが出来る。
（２）壁Ｗを除去した実体の低吸収体積領域Ｃ−Ｗを評価するので、正確な評価が出来る。
（３）従来と互換性のある評価値を出力することも出来る。
【００２０】
【発明の効果】
本発明の肺の低吸収領域評価方法および装置によれば、肺の低吸収領域（低吸収体積領域）を三次元的に抽出し、それらの体積を基に評価値を算出するから、低吸収領域（低吸収体積領域）を正確に評価することが出来る。
【図面の簡単な説明】
【図１】本発明の一実施形態に係るＸ線ＣＴ装置を示す構成図である。
【図２】本発明の一実施形態に係る低吸収領域計測処理を示すフロー図である。
【図３】肺野体積領域と低吸収体積領域の説明図である。
【図４】従来の低吸収領域計測処理を示すフロー図である。
【図５】肺野領域とＣＴデータの説明図である。
【符号の説明】
１００　　　　　　Ｘ線ＣＴ装置
１　　　　　　　　スキャナ装置
２　　　　　　　　処理装置
２ａ　　　　　　　ＣＴデータ取得部
２ｂ　　　　　　　ボリュームデータ作成部
２ｃ　　　　　　　肺野体積領域抽出部
２ｄ　　　　　　　低吸収体積領域抽出部
２ｅ　　　　　　　壁除去部
２ｆ　　　　　　　評価値算出部
３　　　　　　　　大容量記憶装置
４　　　　　　　　表示モニタ
５　　　　　　　　入力装置
【図１】
2004049498.tif 000003
【図２】
2004049498.tif 000004
【図３】
2004049498.tif 000005
【図４】
2004049498.tif 000006
【図５】
2004049498.tif 000007