(57)【要約】【課題】肺を正確に抽出することである。【解決手段】実施形態の画像処理装置は、認識部と、作成部と、表示制御部とを備える。認識部は、肺を含む３次元画像データにおいて、グラディエントとキャビティ構造に対する感度を上げることより、胸壁、血管、結節、及び腫瘍に由来するデータを認識する。作成部は、前記認識部による認識結果に基づいて、前記胸壁に由来するデータを除去する除去処理と、前記血管、前記結節、及び前記腫瘍に由来するデータを残留させる残留処理とを行って、前記肺の領域に対応する肺画像データを作成する。表示制御部は、前記肺画像データを表示させる。【選択図】図１
2018027302.tif 000002
【特許請求の範囲】
【請求項１】
肺を含む３次元画像データにおいて、グラディエントとキャビティ構造に対する感度を上げることより、胸壁、血管、結節、及び腫瘍に由来するデータを認識する認識部と、
前記認識部による認識結果に基づいて、前記胸壁に由来するデータを除去する除去処理と、前記血管、前記結節、及び前記腫瘍に由来するデータを残留させる残留処理とを行って、前記肺の領域に対応する肺画像データを作成する作成部と、
前記肺画像データを表示させる表示制御部と、
を備える、画像処理装置。
【請求項２】
前記３次元画像データから、前記除去処理及び前記残留処理とは異なる前処理により前記肺の領域を抽出した前処理肺画像データを作成する前処理部を更に備え、
前記作成部は、前記前処理部により作成された前処理肺画像データに対して、前記除去処理及び前記残留処理を行って、前記肺画像データを作成する、
請求項１に記載の画像処理装置。
【請求項３】
前記認識部は、グレー値を指数に含む関数を前記グラディエントに対する感度として用いることで、前記胸壁に由来するデータを認識する、
請求項１又は２に記載の画像処理装置。
【請求項４】
前記認識部は、閾値との比較による曲率値の二値化を前記キャビティ構造に対する感度として用いることで、前記血管、前記結節、及び前記腫瘍に由来するデータを認識する、
請求項１〜３のいずれか一つに記載の画像処理装置。
【請求項５】
前記認識部は、等位集合関数モデルにおいて、前記グラディエントに対する感度から前記キャビティ構造に対する感度を減算した値を等位集合正則化項に乗算することで、前記肺画像データに前記キャビティ構造を含むようにする、
請求項１〜４のいずれか一つに記載の画像処理装置。
【請求項６】
前記作成部は、更に、前記血管、前記結節、及び前記腫瘍の少なくとも一つを選択する指示を受け付け、前記指示により選択された前記血管、前記結節、及び前記腫瘍の少なくとも一つを選択的に前記肺画像データから除去する、
請求項１〜５のいずれか一つに記載の画像処理装置。
【請求項７】
前記作成部は、前記指示により選択された前記血管、前記結節、及び前記腫瘍の少なくとも一つを、更に変更する指示を受け付け、前記変更する指示により選択された前記血管、前記結節、及び前記腫瘍の少なくとも一つを前記肺画像データから除去するか否かを変更する、
請求項６に記載の画像処理装置。
【請求項８】
前記作成部は、等位集合関数モデルを用いて、前記肺画像データを作成する、
請求項１〜７のいずれか一つに記載の画像処理装置。
【請求項９】
前記３次元画像データは、磁気共鳴イメージング装置により撮像されたＭＲＩ画像データである、
請求項１〜８のいずれか一つに記載の画像処理装置。
【請求項１０】
肺を含む３次元画像データにおいて、グラディエントとキャビティ構造に対する感度を上げることより、胸壁、血管、結節、及び腫瘍に由来するデータを認識する認識部と、
前記認識部による認識結果に基づいて、前記胸壁に由来するデータを除去する除去処理と、前記血管、前記結節、及び前記腫瘍に由来するデータを残留させる残留処理とを行って、前記肺の領域に対応する肺画像データを作成する作成部と、
前記肺画像データを表示させる表示制御部と、
を備える、磁気共鳴イメージング装置。
【請求項１１】
肺を含む３次元画像データにおいて、グラディエントとキャビティ構造に対する感度を上げることより、胸壁、血管、結節、及び腫瘍に由来するデータを認識し、
前記認識する処理による認識結果に基づいて、前記胸壁に由来するデータを除去する除去処理と、前記血管、前記結節、及び前記腫瘍に由来するデータを残留させる残留処理とを行って、前記肺の領域に対応する肺画像データを作成し、
前記肺画像データを表示させる、
各処理をコンピュータに実行させる、画像処理プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本実施形態の実施形態は、画像処理装置、磁気共鳴イメージング装置及び画像処理プログラムに関する。
【背景技術】
【０００２】
近年、画像認識技術の発展に伴い、医用画像から臓器をセグメンテーションするセグメンテーション技術が注目されている。例えば、肺のセグメンテーションが、肺実質の可視化及び定量分析には極めて重要である。且つ、３Ｄ肺セグメンテーション技術及び肺血管除去技術は、診断過程の加速にも重要な意味を持っている。
【０００３】
従来の技術では、ＭＲＩ（Magnetic  Resonance  Imaging：磁気共鳴イメージング）装置は肺のレンダリングに適せず、頭部のレンダリングに一層適する。従って、通常の場合、ＣＴ（Computerized  Tomography：Ｘ線コンピュータ断層撮影術）装置で採られた画像を用いて、肺のセグメンテーションを行う。現在、ＣＴ画像に基づいた肺のセグメンテーション技術及び肺血管抽出技術がたくさんある。
【０００４】
例えば、特許文献１には、複数器官のセグメンテーション方法及び３次元ＣＴ画像の肺のセグメンテーションについて開示されている。そのうち、マージナルスペース学習の方法が等位集合用の初期化メッシュの生成に用いられ、該等位集合が距離マップに基づいた精確なセグメンテーションの取得に用いられる。
【０００５】
また、特許文献２には、幾何学的変形モデル（geometric  deformation  model）に基づいた３次元肺血管画像のセグメンテーション方法について開示されている。そのうち、高解像度且つ高コントラストのＣＴ画像に対し、等位集合関数が肺血管のセグメンテーションに用いられる。
【０００６】
しかしながら、特許文献１及び特許文献２におけるセグメンテーション方法は、ＭＲＩ装置で採られた画像に適しない。ＭＲＩ装置の場合、従来の技術では、３Ｄ肺セグメンテーション技術及び肺血管除去を精確に施すのが未だ困難である。
【０００７】
例えば、非特許文献１には、肺のセグメンテーションに適用できる方法について開示されている。
【先行技術文献】
【特許文献】
【０００８】
【特許文献１】米国特許第９０４２６２０号明細書
【特許文献２】中国特許第１０２２４３７５９号明細書
【非特許文献】
【０００９】
【非特許文献１】“Distance  Regularized  Level  Set  Evolution  and  Its  Application  to  image  Segmentation[J]”,  Chunming  Li,  Chenyang  Xu  etc.  ,  IEEE  Transactions  on  image  processing.  19(20),2010.
【発明の概要】
【発明が解決しようとする課題】
【００１０】
医用画像に対し肺のセグメンテーションを行う場合、先ず、医用画像収集装置で採られた画像データを入力し、該画像データに基づいて、閾値法又は２値化によるセグメンテーション方法などのラフなセグメンテーション方法を用いてラフな肺のセグメンテーションを施し、肺画像を抽出した後、ラフなセグメンテーション結果を微細修正し、より精確なセグメンテーション結果を得る、ことは一般的に用いられる方法である。
【００１１】
しかしながら、従来のセグメンテーション技術においては、肺実質と胸壁との間の境界における低コントラスト領域に対し、感度が高くない。且つ、キャビティ（くぼみ）構造（例えば、血管及び結節／腫瘍）に対する処理が有効ではない。
【００１２】
図１３は従来の技術における、ＭＲＩ装置により撮像された胸部画像に対し肺のセグメンテーション処理を行った結果を示す例示図である。図１３に示すように、微細修正処理が施された後、抽出された肺画像における矢印Ａが指している領域には依然として、胸壁組織のデータが残っている。一方、微細な血管が含まれているが、例えば、矢印Ｂが指している領域など、比較的粗大な血管は、正確に含まれることができていない。同時、矢印Ｃが指している結節領域なども含まれることができない。
【００１３】
本実施形態は、上記事情に鑑み完成した発明であって、肺を正確に抽出することができる画像処理装置及び画像処理方法を提供することを目的とする。
【課題を解決するための手段】
【００１４】
実施形態に係る画像処理装置は、認識部と、作成部と、表示制御部とを備える。認識部は、肺を含む３次元画像データにおいて、グラディエントとキャビティ構造に対する感度を上げることより、胸壁、血管、結節、及び腫瘍に由来するデータを認識する。作成部は、前記認識部による認識結果に基づいて、前記胸壁に由来するデータを除去する除去処理と、前記血管、前記結節、及び前記腫瘍に由来するデータを残留させる残留処理とを行って、前記肺の領域に対応する肺画像データを作成する。表示制御部は、前記肺画像データを表示させる。
【００１５】
また、本実施形態の別の技術案は、画像収集装置により撮像された画像における肺領域を抽出する画像処理方法であって、肺を含む画像データから、肺領域を抽出し、ラフな肺画像を作成するラフな肺画像作成ステップと、前記画像データに基づいて、前記ラフな肺画像を処理し、精確な肺画像を作成する精確な肺画像作成ステップと、を含み、前記精確な肺画像作成ステップでは、グラディエントに対する感度と、キャビティ構造に対する感度とを用いて、前記精確な肺画像には胸壁を含まず、キャビティ構造を含むように処理する画像処理方法である。
【００１６】
本実施形態によれば、肺のセグメンテーションを精確に施すことができ、肺領域を医用画像からより正確にセグメンテーションすることが可能になる。特に、グラディエントに対する感度と、キャビティ構造に対する感度とを組み合わせて用いることにより、画像にける胸壁とキャビティ構造が位置する領域をより感度高く区分することで、該区分に基づいて胸壁の除去及びキャビティ構造の保留をより良く行うことができる。従って、より精確な肺のセグメンテーション処理を行うことが可能になる。
【００１７】
また、本実施形態に係る画像処理装置では、精確な肺画像作成手段が、精確な肺画像からキャビティ構造を選択的に除去するように、選択処理を行う選択手段を更に備えたものであってもよい。
【００１８】
本実施形態によれば、例えば、血管の画素又は結節の画素がキャビティ構造に対する感度に反映しないようにすることで、セグメンテーション結果には血管を含まないが結節を含むように選択することにより、血管及び結節をより良く区分して別々に処理することが可能になる。血管の除去を選択的に且つ信頼性のより高く行うことが可能であるため、ユーザは血管を除去するか否かを選択し、且つ、結節を保留しながら血管を除去することができ、読影の効率が向上する。
【００１９】
更に、本実施形態は、ＭＲＩ装置による肺画像処理の場合、効果が特に突出する。肺に明確なＭＲ信号を得ることが困難であるため、ＭＲＩ装置が肺の撮像に選択し用いられることが少ない。しかし、ＭＲＩ装置による結果は、肺がんの拡散診断という面で積極的な作用を奏する。本実施形態をＭＲＩ装置に適用することにより、ＭＲＩ装置全体の肺セグメンテーションの正確性を向上させることが可能になり、且つ、心臓及び大動脈に対する分析において、肺マスクの生成に用いることもできる。
【図面の簡単な説明】
【００２０】
【図１】図１は、本実施形態の第１の実施形態に係る画像処理装置を示す構造ブロック図。
【図２】図２は、ＭＲＩ装置により撮像された胸部画像を示す例示図。
【図３】図３は、ラフな肺画像を示す例示図。
【図４】図４は、本実施形態の第１の実施形態に係る肺セグメンテーション処理の結果を示す例示図。
【図５】図５は、本実施形態の第１の実施形態に係る肺セグメンテーション処理を示すフローチャート。
【図６】図６は、本実施形態の第２の実施形態に係る画像処理装置を示す構造ブロック図。
【図７】図７は、本実施形態の第２の実施形態に係る肺セグメンテーション処理を示すフローチャート。
【図８】図８は、本実施形態の第２の実施形態に係る肺セグメンテーション処理の結果を示す例示図。
【図９】図９は、ＭＲＩ装置の画像データを処理する場合、構造セレクタの異なる選択によって生じた異なる処理結果を示す例示図。
【図１０】図１０は、ＭＲＩ装置の画像データを処理する場合、構造セレクタの異なる選択によって生じた異なる処理結果を示す別の例示図。
【図１１】図１１は、ＣＴ装置の画像データを処理する場合、構造セレクタの異なる選択によって生じた異なる処理結果を示す例示図。
【図１２】図１２は、ＣＴ装置の画像データを処理する場合、構造セレクタの異なる選択によって生じた異なる処理結果を示す別の例示図。
【図１３】図１３は、従来の技術における肺セグメンテーション処理の結果を示す例示図。
【図１４】図１４は、本実施形態に係る画像処理装置の構成例を示すブロック図である。
【発明を実施するための形態】
【００２１】
本実施形態は、画像を処理する画像処理装置に関する。この画像処理装置は、Ｘ線装置などの画像収集装置に接続され独立したコンピュータなどのＣＰＵ（central  process  unit）を有する設備が、画像処理装置の各機能を有するソフトウェアを実行することにより実現されてもいいし、また、画像処理装置の各機能を実行可能な回路として、ハードウェアの形態により実現されてもいい。且つ、本実施形態に係る画像処理装置は、ＣＴ（Computerized  Tomography）装置または磁気共鳴イメージング装置などの医用画像収集装置における一部として、上述した医用画像収集装置に予め組み込まれてもよい。
【００２２】
以下、本実施形態を実施するための好適な形態について図面を参照して説明する。各実施形態では、画像収集装置としてＭＲＩ装置を用いることが好ましく、本実施形態は、ＭＲＩ装置で採られた医用画像データに対する肺セグメンテーション処理を用いられることが好ましい。但し、ＣＴ装置、超音波装置などの他の医用画像収集装置の肺セグメンテーション処理に適用されることもできる。
【００２３】
また、下記の各実施形態では、採られた画像データが３次元画像データであることを例として説明するが、本実施形態は二次元などの他の次元の画像データに対する処理に適用されることもできる。且つ、本実施形態に用いられている例示図では、より明確に説明するために、３次元画像におけるある断面図が例示として用いられている。
【００２４】
また、異なる実施形態では、同一の部材について同一の符号（記号）が用いられており、重複の説明を適当に省略する。
【００２５】
（第１の実施形態）
図１は、本実施形態の第１の実施形態に係る画像処理装置を示す構造ブロック図である。図１に示すように、画像処理装置１００は少なくとも、第１画像作成部１０と、第２画像作成部２０とを備える。
【００２６】
第１画像作成部１０は、ＭＲＩ装置が胸部に対し撮像した３次元画像データを受信し、閾値法又は２値化によるセグメンテーション方法などのラフなセグメンテーション方法を用いて、受信された画像データから肺領域を抽出し、ラフな肺画像を作成する。第１画像作成部１０は、上述した機能を実現できる回路又はソフトウェアモジュールであってもよい。言い換えると、３次元画像データは、磁気共鳴イメージング装置により撮像されたＭＲＩ画像データである。また、前処理部としての第１画像作成部１０は、３次元画像データから、除去処理及び残留処理とは異なる前処理により肺の領域を抽出した前処理肺画像データを作成する。
【００２７】
図２は、ＭＲＩ装置により撮像された胸部画像を示す例示図である。図３は、ラフな肺画像を示す例示図である。第１画像作成部１０は、図２に示すような３次元画像データから肺領域を抽出し、図３に示すようなラフな肺画像を作成する。ラフな肺画像において、一般的に、肺実質における共通の特徴を有する画素のみが抽出されるため、血管などの肺に含まれている物質の存在を示すことが困難であり、且つ、マージンに対する処理も雑である。
【００２８】
ここで、血管、結節、腫瘍などの肺実質とは異なる器官又は組織をまとめて、「キャビティ（つぼみ）構造」という。
【００２９】
第２画像作成部２０は、ラフな肺画像を更に処理し、ＭＲＩ装置の３次元画像データに基づいて、ラフな肺画像を修正し、精確な肺画像を作成する。第２画像作成部２０は、上述した機能を実現できる回路又はソフトウェアモジュールであってもよい。
【００３０】
例えば、ラフな肺画像の画素で構築された等位集合を初期の等位集合とし、ＭＲＩ装置の３次元画像データのような元画像データ情報と合わせて、等位集合の進化（evolution）とセグメンテーションを行い、このような反復アルゴリズムにより、精確な肺画像を取得する。このような等位集合の進化を表す等位集合関数を含む偏微分方程式を、等位集合関数モデルとする。
【００３１】
ここで、第２画像作成部２０は、等位集合関数モデルを肺セグメンテーションの精確処理モデルとして、ラフな肺画像を更に修正し、マージン領域がより微細化処理され、且つ血管以外のキャビティ構造を含むことが望まれた精確な肺画像を形成する。
【００３２】
上述の目的を実現するために、第１の実施形態では、第２画像作成部２０は、グラディエントに対する感度と、キャビティ構造に対する感度とを用いて、第２画像作成部２０が生成した精確な肺画像には胸壁を含まず、キャビティ構造を含むように処理する特徴感度算出部２１を備える。言い換えると、認識部としての特徴感度算出部２１は、肺を含む３次元画像データにおいて、グラディエントとキャビティ構造に対する感度を上げることより、胸壁、血管、結節、及び腫瘍に由来するデータを認識する。作成部としての第２画像作成部２０は、認識部による認識結果に基づいて、胸壁に由来するデータを除去する除去処理と、血管、結節、及び腫瘍に由来するデータを残留させる残留処理とを行って、肺の領域に対応する肺画像データを作成する。また、作成部としての第２画像作成部２０は、前処理部により作成された前処理肺画像データに対して、除去処理及び残留処理を行って、肺画像データを作成する。
【００３３】
具体的に、特徴感度算出部２１は、元画像データにおけるグラディエントに敏感なグラディエント感度指標ｇと、元画像データにおけるキャビティ構造に敏感なキャビティ構造感度指標ｓｉｇｎ（ｋ）ｋとを含む感度指標を算出する。そのうち、グラディエント感度指標ｇは下記の式（１）で表してもよい。
【００３４】
【数１】
2018027302.tif 000003
【００３５】
式（１）の中、Ｉは元画像データの３Ｄ画像におけるグレー値のマトリックスであり、▽Ｇはカーネル導関数であり、ｍは要求に応じて任意的に設定された係数である。
【００３６】
グラディエント感度指標ｇをｅインデックスの形態に設定することにより、画像グレースケールとの関係が、グレー値Ｉの増加に伴って、ｇの減少する速度も増加するようになる。特に、従来のグラディエント感度指標に比べて、同一のグレー値Ｉの場合、第１の実施形態におけるｇの値がより小さくなるため、グレー値Ｉの増加に伴って、ｇの値が小さくなる速度がより速い。つまり、グラディエント感度指標ｇは、グレー値Ｉを指数に含む関数（指数関数）により表される。言い換えると、認識部としての特徴感度算出部２１は、グレー値を指数に含む関数（グラディエント感度指標ｇ）をグラディエントに対する感度として用いることで、胸壁に由来するデータを認識する。
【００３７】
このように、式（１）のようなインデックス関数をグラディエント感度指標ｇとして用いると、グラディエントの下がりを加速でき、グラディエントに対しより敏感にする。このような感度指標を等位集合関数φに基づいた肺セグメンテーションの精確処理モデルなどに代入すると、生成された精確な肺画像には胸壁を含まないようにすることが可能になる。
【００３８】
また、キャビティ構造感度指標ｓｉｇｎ（ｋ）ｋにおいて、式（２−２）は平均曲率値であり、ｓｉｇｎ（ｋ）は式（２−１）を満たす。
【００３９】
【数２−１】
2018027302.tif 000004
【数２−２】
2018027302.tif 000005
【００４０】
λは要求に応じて任意的に設定された係数（閾値）である。
【００４１】
キャビティ構造感度指標ｓｉｇｎ（ｋ）ｋはキャビティ構造に対し敏感であるため、精確処理モデルにキャビティ構造感度指標ｓｉｇｎ（ｋ）ｋを加えることにより、血管、結節、腫瘍などのキャビティ構造を監視することができる。つまり、キャビティ構造感度指標ｓｉｇｎ（ｋ）ｋは、曲率値（平均曲率値）と閾値λとの比較結果に応じて、「０」又は「１」を出力する関数である。言い換えると、認識部としての特徴感度算出部２１は、閾値との比較による曲率値の二値化をキャビティ構造に対する感度として用いることで、血管、結節、及び腫瘍に由来するデータを認識する。
【００４２】
第１の実施形態では、特徴感度算出部２１は、グラディエント感度指標ｇとキャビティ構造感度指標ｓｉｇｎ（ｋ）ｋとの差分を取ることにより、キャビティ構造の精確処理モデルにおける拡散エネルギーを強化し、精確処理モデルにより処理された精確な肺画像には血管及び結節を含むようにする。言い換えると、認識部としての特徴感度算出部２１は、等位集合関数モデルにおいて、グラディエントに対する感度からキャビティ構造に対する感度を減算した値を等位集合正則化項に乗算することで、肺画像データにキャビティ構造を含むようにする。
【００４３】
例えば、従来の技術において用いられている、下記の式（３）で表されるような、等位集合正則化項（levelset  distance  regularization）に基づいた精確処理モデルを例に説明すると、特徴感度算出部２１は、グラディエント感度指標ｇ及びキャビティ構造感度指標ｓｉｇｎ（ｋ）ｋと、式（３）とを組み合わせて、特に、式（３）にある元のグラディエント関数の代わりに、本実施形態における式（１）で表されるグラディエント感度指標ｇを用いることで、第２画像作成部２０は、式（４）で表される等位集合正則化項に基づいた肺セグメンテーションの精確処理モデルを用いて、ラフな肺画像を処理する。
【００４４】
【数３】
2018027302.tif 000006
【数４】
2018027302.tif 000007
【００４５】
式（３）の中、周知のように、Φは等位集合関数、∂は偏導関数、▽はグラディエント、ｄｉｖは発散度、｜・・・｜は絶対値、ベクトルの長さであり、μ、λ、α、εは何れも対応する係数である。また、関数（下記の式（５））については、上記の非特許文献１などの既存資料の記載を参照できる。
【００４６】
【数５】
2018027302.tif 000008
【００４７】
図４は本実施形態の第１の実施形態に係る肺セグメンテーション処理の結果を示す例示図である。式（４）で表される精確処理モデルを用いて図３に示すようなラフな肺画像を修正することにより、図４に示すような精確な肺画像が取得される。
【００４８】
図４に示すように、グラディエント感度指標ｇの加速作用のため、矢印Ａが指しているマージン領域には胸壁組織のデータが残っていない。且つ、キャビティ構造感度指標ｓｉｇｎ（ｋ）ｋの共同作用により、矢印Ｂ及び矢印Ｃで示されている血管組織、結節などのキャビティ構造が完全に保留されている。
【００４９】
第１の実施形態では、第１画像作成部１０は「ラフな肺画像作成手段」に対応し、第２画像作成部２０は「精確な肺画像作成手段」に対応し、特徴感度算出部２１は「特徴感度指示手段」に対応する。
【００５０】
以下、図５を参照して画像処理装置１００による肺セグメンテーション処理の流れについて説明する。
【００５１】
図５は本実施形態の第１の実施形態に係る肺セグメンテーション処理を示すフローチャートである。図５に示すように、先ず、ステップＳ５０１では、画像処理装置１００がＭＲＩ装置に生じされた３Ｄ画像データ（図２参照）を受信した後、第１画像作成部１０が３Ｄ画像データについて、ラフな肺領域の抽出処理を行い、肺領域をラフに抽出し、ラフな肺画像を作成する（図３参照）（ステップＳ５０２）。
【００５２】
次に、ステップＳ５０３へ進むと、第２画像作成部２０により、式（４）で示されている精確処理モデルを用いてラフな肺画像の精確修正を更に行い、特に、第２画像作成部２０に含まれた感度指標を利用して、胸壁を含まないが血管及び結節を含む精確な肺画像を生成する（図４参照）。
【００５３】
また、第２画像作成部２０が上記の処理を完了した後、使用者が読影できるように、生成された精確な肺画像を出力することができる。
【００５４】
図５に示すような肺セグメンテーション処理によれば、ＭＲＩ装置の画像に対し肺全体のセグメンテーションを行う時の精確性が向上する。出力される肺画像には胸壁が含まれていないがキャビティ構造が完全に保留されるようにすることができ、特に、結節などの肺診断に比較的重要なデータを保留することができ、使用者の読影に寄与する。
【００５５】
（変形例１）
第１の実施形態では、インデックス関数をグラディエント感度指標として用いており、このようなグラディエント感度指標が感度のより高い指示性能を有する。
【００５６】
しかしながら、グレー値Ｉとｇとの増減関係を満足できれば、他のグラディエント感度指標を用いてもよい。例えば、式（６）を用いてグラディエント感度指標としてもよい。
【００５７】
【数６】
2018027302.tif 000009
【００５８】
また、ガウス畳み込みの乗を増加させてもよい。例えば、式（７）で表される関数をグラディエント感度指標として用いても、何れもグラディエントに対する感度を実現できる。
【００５９】
【数７】
2018027302.tif 000010
【００６０】
（変形例２）
第１の実施形態では、第２画像作成部２０は等位集合正則化項（levelset  distance  regularization）に基づいた精確処理モデルを用いているが、他の精確処理モデルを用いてもよい。グラディエント感度指標と、キャビティ構造感度指標とを含まれれば本実施形態を実施することができる。
【００６１】
例えば、式（８）で表されるような一般的な等位集合関数モデルを精確処理モデルとして用いて、生成される精確な肺画像には胸壁を含まず、キャビティ構造を含むようにすることができる。
【００６２】
【数８】
2018027302.tif 000011
【００６３】
また、式（９）で表されるような、単一井戸ポテンシャル方程式（single-well  potential  equation）のみを等位集合正則化項として、式（４）におけるｐが二重井戸ポテンシャル方程式（double-well  potential  equation）で置換されたモデルを精確処理モデルとして用いて、生成される精確な肺画像には胸壁を含まないようにすることができる。
【００６４】
【数９】
2018027302.tif 000012
【００６５】
それらの変形例は同様に、より精確な肺セグメンテーション処理を実現することができる。
【００６６】
（第２の実施形態）
第２の実施形態は第１の実施形態に基づいたものである。第１の実施形態に対する第２の実施形態の相違点は、画像処理装置２００は表示制御部３０を更に備え、且つ、画像処理装置２００における第２画像作成部２０’が構造セレクタ２２を更に備えることにある。以下、第２の実施形態と第１の実施形態との相違点について主に説明し、重複の説明を適当に省略する。
【００６７】
図６は、本実施形態の第２の実施形態に係る画像処理装置を示す構造ブロック図である。図６に示すように、画像処理装置２００は、第１画像作成部１０と、第２画像作成部２０’と、表示制御部３０とを備える。
【００６８】
第１画像作成部１０は、ＭＲＩ装置が胸部に対し撮像した３次元画像データを受信し、閾値法又は２値化セグメンテーション方法などのラフなセグメンテーション方法を用いて、受信された画像データから肺領域を抽出し、ラフな肺画像を作成する。第１画像作成部１０は、上述した機能を実現できる回路又はソフトウェアモジュールであってもよい。
【００６９】
第２画像作成部２０’は、ラフな肺画像を更に処理し、ＭＲＩ装置の３次元画像データに基づいて、ラフな肺画像を修正し、精確な肺画像を作成する。第２画像作成部２０’は、上述した機能を実現できる回路又はソフトウェアモジュールであってもよい。
【００７０】
例えば、ラフな肺画像の画素で構築された等位集合を初期の等位集合とし、ＭＲＩ装置の３次元画像データのような元画像データ情報と合わせて、等位集合の進化とセグメンテーションを行い、このような反復アルゴリズムにより、精確な肺画像を取得する。このような等位集合の進化を表す等位集合関数を含む偏微分方程式を、等位集合関数モデルとする。
【００７１】
表示制御部３０は、第２画像作成部２０’の処理結果を表示するためのものであって、上述した機能を実現できる回路又はソフトウェアモジュールであってもよい。そして、画像処理装置２００の処理結果を直接出力し、表示制御部３０を省略してもよい。言い換えると、表示制御部３０は、肺画像データを表示させる。
【００７２】
以下、第２画像作成部２０’の構成について詳しく説明する。
【００７３】
第２の実施形態では、第２画像作成部２０’は、グラディエントに対する感度と、キャビティ構造に対する感度とを用いて、第２画像作成部２０’が生成した精確な肺画像には胸壁を含まず、キャビティ構造を含むように処理する特徴感度算出部２１と、構造セレクタ２２とを備える。特徴感度算出部２１は、第１の実施形態で説明されたような、グラディエントに敏感なグラディエント感度指標ｇと、キャビティ構造に敏感なキャビティ構造感度指標ｓｉｇｎ（ｋ）ｋとを含む。
【００７４】
構造セレクタ２２は、第２画像作成部２０’が生成した精確な肺画像からキャビティ構造を選択的に除去するように、選択処理を行う。
【００７５】
具体的に、キャビティ構造における血管と結節とを区分して、構造セレクタ２２は、元画像における血管の画素又は結節の画素がキャビティ構造感度指標ｓｉｇｎ（ｋ）ｋに反映しないようにすることにより、精確な肺画像には、前記血管又は前記結節を含まないように選択する。言い換えると、作成部としての構造セレクタ２２は、血管、結節、及び腫瘍の少なくとも一つを選択する指示を受け付け、指示により選択された血管、結節、及び腫瘍の少なくとも一つを選択的に肺画像データから除去する。
【００７６】
第２の実施形態で用いられる構造セレクタ２２は、下記の式（１０）である。ただし、βは０以上の係数である。
【００７７】
【数１０】
2018027302.tif 000013
【００７８】
上述の構造セレクタ２２は、より選択し易くするように異なるキャビティ構造の特徴を区分するためのキャビティ構造モデルＤ（ＴＨ）と、マトリックスにおける異なるエレメントの設定により、除去すべきキャビティ構造を指定するように選択することが可能なキャビティ構造選択マトリックスＷと、を含む。
【００７９】
構造セレクタ２２は、式（１０）を精確処理モデルに代入することで、選択されたキャビティ構造がキャビティ構造感度指標ｓｉｇｎ（ｋ）ｋに反映しないようにすることにより、セグメンテーション結果から選択されたキャビティ構造を除去する。下記の式（１１）では、第２の実施形態で用いられる等位集合正則化項に基づいた精確処理モデルが示される。
【００８０】
【数１１】
2018027302.tif 000014
【００８１】
また、キャビティ構造モデルＤ（ＴＨ）について、Ｈは画像の本質の特徴を表している、元画像データのヘッセ行列であって、偏導関数∂が用いられる。ＴＨはヘッセ行列の３つの特徴値λ１、λ２、λ３であって、且つλ１＜λ２＜λ３である。３つの固有ベクトルはそれぞれ当該領域に対し相互に直交する３つの方向を表し、固有ベクトルに対応する特徴値の逆数が当該方向の延伸度を表すものである。特徴値の逆数が大きければ、当該方向の延伸度が大きくなる。下記の式（１２−１）及び式（１２−２）で表される関係がある。
【００８２】
【数１２−１】
2018027302.tif 000015
【数１２−２】
2018027302.tif 000016
【００８３】
ただし、ｎは要求に応じて任意的に設定された係数である。λ´１、λ´２、λ´３の異なる組合せにより、キャビティ構造を区分することが可能になる。異なるキャビティ構造に対応する組合せは、下記の表１で示される。
【００８４】
【表１】
2018027302.tif 000017
【００８５】
表１の中、Ｈ＝ｈｉｇｈ，Ｌ＝ｌｏｗ（約０）。
【００８６】
このように、元画像のヘッセ行列を用いることにより、異なるキャビティ構造を区分することが可能になる。
【００８７】
また、キャビティ構造選択マトリックスＷについて、それはキャビティ構造モデルＤ（ＴＨ）と乗算することにより、キャビティ構造モデルＤ（ＴＨ）で表される異なるキャビティ構造を選択するための値である。マトリックスＷを異なる値に定義することにより、異なるキャビティ構造を選択することができる。
【００８８】
例えば、「０」、「１」を用いてマトリックスＷの値を定義する場合、Ｗが取れる値は、下記の式（１３−１），（１３−２），（１３−３），（１３−４）のうちいずれかである。
【００８９】
【数１３−１】
2018027302.tif 000018
【数１３−２】
2018027302.tif 000019
【数１３−３】
2018027302.tif 000020
【数１３−４】
2018027302.tif 000021
【００９０】
如何なる除去を行わず肺全体の画像を用いる場合、Ｗの値を式（１３−１）とする。血管を除去し結節を保留するように選択する場合、Ｗの値を式（１３−２）とする。結節を除去し血管を保留するように選択する場合、Ｗの値を式（１３−３）とする。結節を除去し血管も除去するように選択する場合、Ｗの値を式（１３−４）とする。例えば、管状（λ１＝Ｌ、λ２＝Ｈ、λ３＝Ｌ）の時、血管を除去することが期待される場合、Ｗの値を式（１３−２）とするだけで、下記の式（１３−５）になって、下記の式（１３−６）となる。
【００９１】
【数１３−５】
2018027302.tif 000022
【数１３−６】
2018027302.tif 000023
【００９２】
図８は、血管を除去することが期待される場合の、肺セグメンテーション処理の結果を示す例示図である。図８に示すように、矢印Ａが指しているマージン領域には胸壁組織のデータが残されていない。且つ、キャビティ構造感度指標及び構造セレクタ２２の共同作用により、矢印Ｂで示されている血管組織が完全に除去され、矢印Ｃで示されている結節が保留されている。
【００９３】
第２の実施形態では、第１画像作成部１０は「ラフな肺画像作成手段」に対応し、第２画像作成部２０’は「精確な肺画像作成手段」に対応し、特徴感度算出部２１は「特徴感度指示手段」に対応し、構造セレクタ２２は「選択手段」に対応し、表示制御部３０は「表示手段」に対応する。
【００９４】
以下、図８に示すような精確な肺画像を取得することを例に、図７を参照して画像処理装置２００による肺セグメンテーション処理の流れについて説明する。
【００９５】
図７は本実施形態の第２の実施形態に係る肺セグメンテーション処理を示すフローチャートである。図７に示すように、先ず、ステップＳ７０１では、画像処理装置２００がＭＲＩ装置に生じられた３Ｄ画像データ（図２参照）を受信した後、第１画像作成部１０が３Ｄ画像データについて、ラフな肺領域の抽出処理を行い、肺領域をラフに抽出し、ラフな肺画像を作成する（図３参照）（ステップＳ７０２）。
【００９６】
次に、ステップＳ７０３へ進むと、第２画像作成部２０’により、式（１１）で示されている精確処理モデルを用いて精確修正を更に行い、感度指標及び構造セレクタを利用し、特に、キャビティ構造感度指標及び構造セレクタを同時に用いて、血管を除去すべきキャビティ構造として選択することで、胸壁と血管を含まず、結節を含む精確な肺画像を生成する。
【００９７】
そして、ステップＳ７０４では、表示制御部３０が生成された図８に示すような精確な肺画像を表示する。
【００９８】
図７に示すような肺セグメンテーション処理により、ＭＲＩ装置の画像に対し肺全体のセグメンテーションを行う時の精確性が向上する。血管の画素又は結節の画素がキャビティ構造に対する感度に反映しないようにすることで、セグメンテーション結果に血管が含まれていないが結節が含まれるように選択することにより、血管及び結節をより良く区分して別々に処理することが可能になる。血管の除去を選択的に且つ信頼性のより高く行うことが可能であるため、ユーザは血管を除去するか否かを選択し、且つ、結節を保留しながら血管を除去することができ、読影の効率が向上する。
【００９９】
図８に示すような選択タイプ以外に、他のタイプもたくさんある。例えば、図９は、ＭＲＩ装置の画像データを処理する場合、構造セレクタの異なる選択によって生じた異なる処理結果を示す例示図である。そのうち、図９の（１）は肺全体をセグメンテーションすること、即ち、全てのキャビティ構造を含む場合を示し、図９の（２）は血管を除去するが結節を保留する場合を示し、図９の（３）は結節を除去するが血管を保留する場合を示す。
【０１００】
また、図１０は、ＭＲＩ装置の画像データを処理する場合、構造セレクタの異なる選択によって生じた異なる処理結果を示す別の例示図である。そのうち、図１０の（１）は肺全体をセグメンテーションすること、即ち、全てのキャビティ構造を含む場合を示し、図１０の（２）は血管を除去するが腫瘍を保留する場合を示し、図１０の（３）は腫瘍を除去するが血管を保留する場合を示す。
【０１０１】
（変形例１）
第２の実施形態では、第２画像作成部２０’が、構造セレクタ２２による選択の変更を受け付ける選択変更受付手段を更に備えたものであってもよい。
【０１０２】
つまり、選択変更受付手段はＷの値を変更することにより、キャビティ構造に対する感度に反映するキャビティ構造のタイプを変更し、表示制御部３０に表示された精確な肺画像のタイプを変更することができる。言い換えると、作成部としての第２画像作成部２０’は、指示により選択された血管、結節、及び腫瘍の少なくとも一つを、更に変更する指示を受け付け、変更する指示により選択された血管、結節、及び腫瘍の少なくとも一つを肺画像データから除去するか否かを変更する。
【０１０３】
よって、ユーザは表示しようとする精確な肺セグメンテーション結果のタイプを、いつでも変更することができる。
【０１０４】
（変形例２）
第１の実施形態の変形例も同様に、第２の実施形態に適用できる。例えば、式（１４）で表されるような一般的な等位集合関数モデルを精確処理モデルとして用いてもよい。言い換えると、作成部としての第２画像作成部２０は、等位集合関数モデルを用いて、肺画像データを作成する。
【０１０５】
【数１４】
2018027302.tif 000024
【０１０６】
また、式（１５）で表されるような、単一井戸ポテンシャル方程式（single-well  potential  equation）のみを等位集合正則化項として、式（１１）におけるｐが二重井戸ポテンシャル方程式（double-well  potential  equation）で置換されたモデルを精確処理モデルとして用いてもよい。
【０１０７】
【数１５】
2018027302.tif 000025
【０１０８】
それらの変形例は同様に、より精確な肺セグメンテーション処理を実現することができる。
【０１０９】
（変形例３）
第１の実施形態でも第２の実施形態でも、ＭＲＩ装置の３次元画像を対象として処理したが、本実施形態はＣＴ装置により撮像された画像を処理することもできる。
【０１１０】
図１１は、ＣＴ装置の画像データを処理する場合、構造セレクタの異なる選択によって生じた異なる処理結果を示す例示図である。そのうち、図１１の（１）は肺全体をセグメンテーションすること、即ち、全てのキャビティ構造を含む場合を示し、図１１の（２）は血管を除去するが結節を保留する場合を示し、図１１の（３）は結節を除去するが血管を保留する場合を示す。
【０１１１】
図１２は、ＣＴ装置の画像データを処理する場合、構造セレクタの異なる選択によって生じた異なる処理結果を示す別の例示図である。そのうち、図１２の（１）は肺全体をセグメンテーションすること、即ち、全てのキャビティ構造を含む場合を示し、図１２の（２）は血管を除去するが結節を保留する場合を示し、図１２の（３）は結節を除去するが血管を保留する場合を示す。
【０１１２】
（変形例４）
本実施形態では、等位集合の進化過程における拡張又は収縮を影響するレベルを変更することにより、重み付け面積項の大きさを変更するものであって、キャビティ構造に対する感度指標を加え、構造セレクタにより結果としての構造を選択できるという目的を達成した。本実施形態の発明構想さえ実現できれば、用いられる具体的なモデルを任意的に変更しても構わない。
【０１１３】
また、上記の各数式における係数を画像の要求及び環境に応じて選択することができる。例えば、画像範囲が０〜２５５に圧縮された場合、下記のパラメータデータを用いることが好ましい：μ＝０．００１、λ＝５．０、ε＝１．５５、α＝−６．０、β＝１．０、γ＝−０．５、ｍ＝０．５、ｎ＝１．０。
【０１１４】
そして、進化過程におけるイテレーション数が８０であり、ステップ・サイズｔが１０．０であることが望ましい。各パラメータは対応する実際情況に応じて調整されることができ、範囲が特に明確に規定されていない。実際の効果に応じて調整してもよいが、μ、λ、ε、β、ｍ、ｎ＞０，α、γ＜０という条件を満足する必要がある。
【０１１５】
（変形例５）
上記の各実施形態では、３次元画像を例に説明したが、本実施形態は同様に他の次元の画像に適用できる。例えば、２Ｄ画像に適用される場合、２Ｄ画像について管状及び円形状の２種の構造のみを考慮すればよいので、構造セレクタが異なり、ヘッセ行列及び他の対応する式は下記の式（１６−１），（１６−２）で示されるようになる。
【０１１６】
【数１６−１】
2018027302.tif 000026
【数１６−２】
2018027302.tif 000027
【０１１７】
それに応じて、Ｗが取れる値の範囲が、下記の式（１７−１），（１７−２），（１７−３），（１７−４）のうちいずれかになる。
【０１１８】
【数１７−１】
2018027302.tif 000028
【数１７−２】
2018027302.tif 000029
【数１７−３】
2018027302.tif 000030
【数１７−４】
2018027302.tif 000031
【０１１９】
本実施形態の画像処理装置は、プロセッサとメモリとを含んでなる処理装置であってもよい。メモリには、上記の各実施形態において説明された動作に対応する指令情報が記憶されており、プロセッサはメモリに記憶された指令情報を読み取り、実行することにより、本実施形態を具体的に実施する。
【０１２０】
図１４は、本実施形態に係る画像処理装置の構成例を示すブロック図である。画像処理装置３００は、例えば、パーソナルコンピュータやタブレット端末等の汎用装置である。画像処理装置３００は、入力装置３０１、ディスプレイ３０２、記憶回路３１０、及び処理回路３２０を備える。入力装置３０１、ディスプレイ３０２、記憶回路３１０、及び処理回路３２０は、相互に通信可能に接続される。
【０１２１】
入力装置３０１は、マウス、キーボード、タッチパネル等、操作者からの各種の指示や設定要求を受け付けるための入力装置である。ディスプレイ３０２は、医用画像を表示したり、操作者が入力装置３０１を用いて各種設定要求を入力するためのＧＵＩを表示したりする表示装置である。
【０１２２】
記憶回路３１０は、例えば、ＮＡＮＤ（Not  AND）型フラッシュメモリやＨＤＤ（Hard  Disk  Drive）であり、医用画像データやＧＵＩを表示するための各種のプログラムや、当該プログラムによって用いられる情報を記憶する。
【０１２３】
処理回路３２０は、画像処理装置３００における処理全体を制御する電子機器（プロセッサ）である。処理回路３２０は、特徴感度算出部２１及び第２画像作成部２０に対応する処理を実行する。また、処理回路３２０は、第１画像作成部１０、構造セレクタ２２、又は表示制御部３０に対応する処理を実行してもよい。
【０１２４】
本実施形態の画像処理装置は、各実施形態に説明された機能を実現できる回路として、医用装置に組み込まれてもよいし、コンピュータが実行可能なプログラムとして、磁気ディスク（フロッピーディスク（登録商標：floppy）、ハードディスクなど）、コンパクトディスク（ＣＤ−ＲＯＭ、ＤＶＤなど）、光ディスク（ＭＯ）、半導体メモリなどの記憶媒体に記憶され発行されてもよい。
【０１２５】
そして、記憶媒体からコンピュータにインストールされたプログラムによる指示に基づいてコンピュータ上に実行されるＯＳ（オペレイティングシステム）、データベース管理ソフトウェア、ネットワークソフトウェアなどのＭＷ（ミドルワーク）なども、上述した実施形態を実現するための各処理の一部を実行することができる。
【０１２６】
本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。
【符号の説明】
【０１２７】
１００，２００            画像処理装置
１０                      第１画像作成部
２０、２０’              第２画像作成部
２１                      特徴感度差出部
２２                      構造セレクタ
３０                      表示制御部
【図１】
2018027302.tif 000032
【図２】
2018027302.tif 000033
【図３】
2018027302.tif 000034
【図４】
2018027302.tif 000035
【図５】
2018027302.tif 000036
【図６】
2018027302.tif 000037
【図７】
2018027302.tif 000038
【図８】
2018027302.tif 000039
【図９】
2018027302.tif 000040
【図１０】
2018027302.tif 000041
【図１１】
2018027302.tif 000042
【図１２】
2018027302.tif 000043
【図１３】
2018027302.tif 000044
【図１４】
2018027302.tif 000045