(57)【要約】  局所的に存在する微小な濃度値変化に基づいて対象領域を抽出するとともに、抽出された対象領域を明確に表示することが可能な画像処理装置、領域抽出方法、及び画像処理方法を提供するために、画像処理装置100は、画像から血管領域Aを抽出し、血管領域Aの平均濃度値よりCT値が小さい領域をソフトプラーク領域Bとして抽出する。また、未抽出のソフトプラークについては、血管領域Aとソフトプラーク領域Bとの差分領域に画素対を設定して、画素対間の各画素Pjについて、画素対のCT値よりわずかに小さい値よりさらに小さい画素値であるか否かを判定する。これにより、画素値が局所的にわずかに変化する部分をソフトプラークとして抽出する。
2014132829.tif 000002
【特許請求の範囲】
【請求項１】
画像データを入力する入力部と、
前記入力部により入力された画像データから第1領域を抽出する第1領域抽出部と、
前記第1領域内の濃度値に基づく閾値により前記第1領域を閾値処理し、前記第1領域内から第2領域を抽出する第2領域抽出部と、
前記第1領域と前記第2領域との差分領域内にある2画素の組み合わせである画素対を設定し、各画素対の間にある画素を注目画素とし、前記注目画素の画素値と前記画素対の少なくともいずれか一方の画素値との差が所定の値より大きい場合に、前記注目画素を第3領域として抽出する第3領域抽出部と、
前記第2領域と前記第3領域とを合わせた領域を対象領域として得る対象領域設定部と、
を備えることを特徴とする画像処理装置。
【請求項２】
前記画像データは3次元原画像であり、
前記3次元原画像に基づいて前記対象領域とその周囲を含む範囲の2次元画像を生成する画像生成部と、
前記画像生成部により生成された2次元画像を表示する表示部と、
を更に備えることを特徴とする請求項1に記載の画像処理装置。
【請求項３】
前記第1領域は血管領域であり、前記対象領域はソフトプラーク領域であることを特徴とする請求項1に記載の画像処理装置。
【請求項４】
前記画像生成部は、
前記対象領域の周囲について参照用2次元画像を生成するとともに前記対象領域の陰影付け画像を生成し、前記参照用2次元画像上に前記対象領域の陰影付け画像を重畳または合成した合成画像を生成することを特徴とする請求項2に記載の画像処理装置。
【請求項５】
前記第1領域に対する前記対象領域の占める割合を表示する割合表示部を更に備えることを特徴とする請求項1に記載の画像処理装置。
【請求項６】
前記合成画像の任意の位置を指定する位置指定部と、
前記位置指定部により指定された位置での前記第1領域に直交する断面画像に前記対象領域を重畳表示する断面画像表示部と、
を備えることを特徴とする請求項4に記載の画像処理装置。
【請求項７】
前記対象領域の形状を評価する形状評価部を更に備え、
前記画像生成部は、前記形状評価部による評価の結果、前記対象領域が前記第1領域の内周部と接していない場合に前記対象領域が優先的に表示されるように表示優先度を設定し、前記3次元原画像に基づく2次元画像を生成することを特徴とする請求項2に記載の画像処理装置。
【請求項８】
前記画像生成部は、更に、視点と前記対象領域との位置関係に基づいて前記表示優先度を設定することを特徴とする請求項7に記載の画像処理装置。
【請求項９】
画像から血管領域を抽出する血管領域抽出部と、
前記血管領域内のソフトプラーク領域を抽出するソフトプラーク領域抽出部と、
前記ソフトプラーク領域の形状を評価する形状評価部と、
前記ソフトプラーク領域の形状評価結果に基づいてソフトプラーク領域の表示優先度を設定する表示優先度設定部と、
前記表示優先度に基づいて前記血管領域の2次元画像を生成する画像生成部と、
を備えることを特徴とする画像処理装置。
【請求項１０】
前記形状評価部は、前記血管領域の断面画像に基づいて前記ソフトプラーク領域が血管壁に接しているか否かを判定することを特徴とする請求項9に記載の画像処理装置。
【請求項１１】
前記形状評価部は、更に、前記血管領域の3次元原画像に基づいて前記ソフトプラーク領域の周囲の各方向を3次元に探索し、前記ソフトプラーク領域が前記血管壁に延伸し、到達するか否かを判定することを特徴とする請求項10に記載の画像処理装置。
【請求項１２】
前記表示優先度設定部は、
前記ソフトプラーク領域が前記血管壁と接していない場合は、接している場合より高い示優先度を設定することを特徴とする請求項910記載の画像処理装置。
【請求項１３】
コンピュータを用いて画像データから対象領域を抽出する領域抽出方法であって、
画像データを入力するステップと、
入力された画像データを閾値処理して第1領域を抽出するステップと、
前記第1領域内の濃度値に基づく閾値により前記第1領域を閾値処理し、前記第1領域内から第2領域を抽出するステップと、
前記第1領域と前記第2領域との差分領域内にある2画素の組み合わせである画素対を設定し、各画素対の間にある画素を注目画素とし、前記注目画素の画素値と前記画素対の少なくともいずれか一方の画素値との差が所定の値より大きい場合に、前記注目画素を第3領域として抽出するステップと、
前記第2領域と前記第3領域とを合わせた領域を対象領域として得るステップと、
を含むことを特徴とする領域抽出方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像処理装置、領域抽出方法、及び画像処理方法に係り、詳細には、CT画像またはMR画像における血管ソフトプラークの抽出及び描画に関する。
【背景技術】
【０００２】
血管中に存在するソフトプラークは不安定で血管壁から離脱しやすく、血栓化を引き起こすリスクが大きいと言われている。そのため、ソフトプラークが存在することやその分布を視認しやすく表示することが望まれる。【０００３】
例えば特許文献1には、血管に付随するプラークを解析する画像解析装置について記載されている。この画像解析装置では、画像上の所望の領域を指定すると、指定された領域内の画像の画素値分布に応じてプラークの不安定さを表す指標を表示する。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開2011-115481号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、血管からソフトプラークを正確に抽出することは困難であった。特にCT画像の場合、ソフトプラークのCT値と周囲の血液のCT値は近い。ソフトプラークは局所的に存在する微小な濃度変化として表れるため、周囲の血管から分離して抽出することが困難であった。また、臓器の平均CT値は撮影装置の特性やスライス位置にも依存して変化する。そのため、外部から閾値を設定しにくく、単純な閾値処理で血管とソフトプラークとを分離するのは困難であった。【０００６】
本発明は、以上の問題点に鑑みてなされたものであり、局所的に存在する微小な濃度値変化に基づいて対象領域を抽出するとともに、抽出された対象領域を明確に表示することが可能な画像処理装置、及び領域抽出方法を提供することを目的とする。
【課題を解決するための手段】
【０００７】
前述した目的を達成するために本発明は、画像データを入力する入力部と、前記入力部により入力された画像データから第1領域を抽出する第1領域抽出部と、前記第1領域内の濃度値に基づく閾値により前記第1領域を閾値処理し、前記第1領域内から第2領域を抽出する第2領域抽出部と、前記第1領域と前記第2領域との差分領域内にある2画素の組み合わせである画素対を設定し、各画素対の間にある画素を注目画素とし、前記注目画素の画素値と前記画素対の少なくともいずれか一方の画素値との差が所定の値より大きい場合に、前記注目画素を第3領域として抽出する第3領域抽出部と、前記第2領域と前記第3領域とを合わせた領域を対象領域として得る対象領域設定部と、を備えることを特徴とする画像処理装置である。【０００８】
また、画像から血管領域を抽出する血管領域抽出部と、前記血管領域内のソフトプラーク領域を抽出するソフトプラーク領域抽出部と、前記ソフトプラーク領域の形状を評価する形状評価部と、前記ソフトプラーク領域の形状評価結果に基づいてソフトプラーク領域の表示優先度を設定する表示優先度設定部と、前記表示優先度に基づいて前記血管領域の2次元画像を生成する画像生成部と、を備えることを特徴とする画像処理装置である。【０００９】
また、コンピュータを用いて画像データから対象領域を抽出する領域抽出方法であって、画像データを入力するステップと、入力された画像データを閾値処理して第1領域を抽出するステップと、前記第1領域内の濃度値に基づく閾値により前記第1領域を閾値処理し、前記第1領域内から第2領域を抽出するステップと、前記第1領域と前記第2領域との差分領域内にある2画素の組み合わせである画素対を設定し、各画素対の間にある画素を注目画素とし、前記注目画素の画素値と前記画素対の少なくともいずれか一方の画素値との差が所定の値より大きい場合に、前記注目画素を第3領域として抽出するステップと、前記第2領域と前記第3領域とを合わせた領域を対象領域として得るステップと、を含むことを特徴とする領域抽出方法である。【００１０】
また、コンピュータを用いた画像処理方法であって、被検体の画像から血管領域を抽出するステップと、前記血管領域内のソフトプラーク領域を抽出するステップと、前記ソフトプラーク領域の形状を評価するステップと、前記ソフトプラーク領域の形状評価結果に基づいてソフトプラーク領域の表示優先度を設定するステップと、前記表示優先度に基づいて前記血管領域の2次元画像を生成するステップと、を含むことを特徴とする画像処理方法である。
【発明の効果】
【００１１】
本発明により、局所的に存在する微小な濃度値変化に基づいて対象領域を抽出するとともに、抽出された対象領域を明確に表示することが可能な画像処理装置、及び領域抽出方法を提供することができる。
【図面の簡単な説明】
【００１２】
【図１】画像処理装置100の全体構成を示す図
【図２】画像処理装置100の機能構成を示すブロック図
【図３】(a)血管領域A、(b)ソフトプラーク領域Bの抽出について説明する図
【図４】(a)画素対(Pi、Pk)の設定例、(b)局所的な画素値の微小な変化を表すグラフ
【図５】操作画面3の一例
【図６】画素対を用いた領域抽出条件判定について説明する図
【図７】対象領域(ソフトプラーク領域B、C)の分布例
【図８】抽出した対象領域41pを断層像411上に重畳表示した画像41
【図９】抽出した対象領域42pを骨抜きMIP画像421に重畳または合成表示した2次元画像42及びソフトプラークの面積割合表示42aの例
【図１０】抽出した対象領域43pをMIP画像431に重畳または合成表示した2次元画像43、断面表示43c、及びソフトプラークの面積割合表示43aの例
【図１１】本発明に係る領域抽出処理の手順を説明するフローチャート
【図１２】血管5内のソフトプラークの形状と種類について説明する図
【図１３】画像生成処理の手順を説明するフローチャート
【図１４】形状評価処理(1)の手順を説明するフローチャート
【図１５】表示優先度設定処理の手順を説明するフローチャート
【図１６】ソフトプラークの形状に基づいて設定された表示優先度に応じた描画について説明する図(非接触プラーク53と接触プラーク52とが存在する場合)
【図１７】ソフトプラークの形状に基づいて設定された表示優先度に応じた描画について説明する図(複数の接触プラーク52a、52bが存在する場合)
【図１８】ソフトプラークの形状に基づいて設定された表示優先度に応じた描画について説明する図(非接触プラーク53と接触プラーク52とが存在する場合
【図１９】ソフトプラークの形状に基づいて設定された表示優先度に応じた描画について説明する図(複数の接触プラーク52a、52bが存在する場合)
【図２０】視点60bを血管内に設定した描画例
【図２１】ソフトプラークの形状に応じた色分け表示例。(a)MIP画像451に合成表示(b)骨抜きMIP画像461合成表示
【図２２】片側固定非接触プラーク53Aについて説明する図
【図２３】両側固定非接触プラーク53Bについて説明する図
【図２４】両側固定非接触プラーク53Bの別の例について説明する図
【図２５】形状評価処理(2)の手順を説明するフローチャート
【発明を実施するための形態】
【００１３】
以下図面に基づいて、本発明の実施形態を詳細に説明する。
【００１４】
［第1の実施の形態］
まず、図1を参照して、本発明の画像処理装置100を適用した画像処理システム1の構成について説明する。
【００１５】
図1に示すように、画像処理システム1は、表示装置107、入力装置109を有する画像処理装置100と、画像処理装置100にネットワーク110を介して接続される画像データベース111と、医用画像撮影装置112とを備える。
【００１６】
画像処理装置100は、画像生成、画像解析等の処理を行うコンピュータである。画像処理装置100は、図1に示すように、CPU(Central Processing Unit)101、主メモリ102、記憶装置103、通信インタフェース(通信I/F)104、表示メモリ105、マウス108等の外部機器とのインタフェース(I/F)106を備え、各部はバス113を介して接続されている。
【００１７】
CPU101は、主メモリ102または記憶装置103等に格納されるプログラムを主メモリ102のRAM上のワークメモリ領域に呼び出して実行し、バス113を介して接続された各部を駆動制御し、画像処理装置100が行う各種処理を実現する。
【００１８】
CPU101は、画像から血管のソフトプラーク領域を抽出する領域抽出処理(図3参照)を実行する。領域抽出処理の詳細については後述する。
【００１９】
主メモリ102は、ROM(Read Only Memory)、RAM(Random Access Memory)等により構成される。ROMはコンピュータのブートプログラムやBIOS等のプログラム、データ等を恒久的に保持している。また、RAMは、ROM、記憶装置103等からロードしたプログラム、データ等を一時的に保持するとともに、CPU101が各種処理を行う為に使用するワークメモリ領域を備える。
【００２０】
記憶装置103は、HDD(ハードディスクドライブ)や他の記録媒体へのデータの読み書きを行う記憶装置であり、CPU101が実行するプログラム、プログラム実行に必要なデータ、OS(オペレーティングシステム)等が格納される。プログラムに関しては、OSに相当する制御プログラムや、アプリケーションプログラムが格納されている。これらの各プログラムコードは、CPU101により必要に応じて読み出されて主メモリ102のRAMに移され、各種の手段として実行される。
【００２１】
通信I/F104は、通信制御装置、通信ポート等を有し、画像処理装置100とネットワーク110との通信を媒介する。また通信I/F104は、ネットワーク110を介して、画像データベース111や、他のコンピュータ、或いは、X線CT装置、MRI装置等の医用画像撮影装置112との通信制御を行う。
【００２２】
I/F106は、周辺機器を接続させるためのポートであり、周辺機器とのデータの送受信を行う。例えば、マウス108やスタイラスペン等のポインティングデバイスをI/F106を介して接続させるようにしてもよい。
【００２３】
表示メモリ105は、CPU101から入力される表示データを一時的に蓄積するバッファである。蓄積された表示データは所定のタイミングで表示装置107に出力される。
【００２４】
表示装置107は、液晶パネル、CRTモニタ等のディスプレイ装置と、ディスプレイ装置と連携して表示処理を実行するための論理回路で構成され、表示メモリ105を介してCPU101に接続される。表示装置107はCPU101の制御により表示メモリ105に蓄積された表示データを表示する。
【００２５】
入力装置109は、例えば、キーボード等の入力装置であり、操作者によって入力される各種の指示や情報をCPU101に出力する。操作者は、表示装置107、入力装置109、及びマウス108等の外部機器を使用して対話的に画像処理装置100を操作する。
【００２６】
ネットワーク110は、LAN(Local Area Network)、WAN(Wide Area Network)、イントラネット、インターネット等の各種通信網を含み、画像データベース111やサーバ、他の情報機器等と画像処理装置100との通信接続を媒介する。
【００２７】
画像データベース111は、医用画像撮影装置112によって撮影された画像データを蓄積して記憶するものである。図1に示す画像処理システム1では、画像データベース111はネットワーク110を介して画像処理装置100に接続される構成であるが、画像処理装置100内の例えば記憶装置103に画像データベース111を設けるようにしてもよい。
【００２８】
次に、図2を参照して画像処理装置100の機能構成について説明する。
【００２９】
図2に示すように、画像処理装置100は、入力部21、第1領域抽出部22、閾値算出部23、第2領域抽出部24、第3領域抽出部25、対象領域設定部26、画像生成部27、及び表示部28を有する。
【００３０】
入力部21は、画像データを入力する。画像データは、被検体をX線CT装置やMR装置等を用いて撮影した複数の断層像である。これらの複数の断層像を積み上げたものを3次元原画像データと呼ぶ。以下、入力画像がCT画像である場合を例として説明する。
【００３１】
第1領域抽出部22は、入力部21により入力された画像データから第1領域を抽出する。本実施形態では、第1領域を血管領域とする。第1領域抽出部22は血管領域を抽出するために適切な閾値による閾値処理を行い、複数の血管を含む二値画像を得る。更に領域の円形度等や所在位置等を指定して特定の血管を抽出する。なお、この血管領域の抽出方法は一例であり、その他の方法を用いて血管領域を抽出してもよい。
【００３２】
例えば、操作者のマウス指示により特定の血管を抽出してもよい。図3(a)は第1領域抽出部22により抽出された血管領域(第1領域)Aの一例である。血管領域(第1領域)Aには、血液の領域とソフトプラーク等の領域とが混在することがあるため、濃度値は必ずしも一様ではない。
【００３３】
閾値算出部23は、第1領域抽出部22により抽出した血管領域A内の濃度値に基づいて第2領域を抽出するための閾値を算出する。本実施の形態では、ソフトプラーク(第2領域)を抽出するための閾値として、血管領域Aの平均濃度値を算出する。なお、算出する閾値は、血管領域Aの平均濃度値に限定されず、血管領域A内の濃度値から求められるその他の統計値としてもよい。
【００３４】
第2領域算出部24は、閾値算出部23によって算出された閾値(平均濃度値)により第1領域(血管領域A)を閾値処理し、第2領域を抽出する。例えば、閾値(平均濃度値)より画素値が小さい画素を第2領域として抽出する。図3(b)に示すようなソフトプラーク領域Bが第2領域として抽出される。なお、ソフトプラークの濃度値はムラがあるため、上述の閾値処理ではソフトプラークであっても抽出されない画素が残る。
【００３５】
第3領域抽出部25は、第2領域抽出部24により抽出しそこなったソフトプラーク領域(第3領域)を抽出する。第3領域抽出部25は、図2に示すように、差分領域算出部251、画素対設定部252、比較部253を含む。
【００３６】
差分領域算出部251は、第1領域(血管領域A)と第2領域(ソフトプラーク領域B)との差分領域を求める。画素対設定部252は、差分領域内にある2画素の組み合わせである画素対(Pi、Pk)を設定する。比較部253は、画素対設定部252により設定された各画素対(Pi、Pk)の間にある画素Pjを注目画素とし、注目画素Pjの画素値Gjと、画素対の画素値Gi、Gkよりわずかに小さい値Gi−ΔG、Gk−ΔGとを比較する。注目画素Pjの画素値Gjの方が小さい場合に、この注目画素をソフトプラーク領域として抽出する。
【００３７】
すなわち、第3領域抽出部25は、以下の式(1)、(2)を判定して、少なくともいずれか一方の式を満たす場合に、当該注目画素Pjを第3領域(ソフトプラーク領域C)として抽出する。
【００３８】
Gj＜Gi−ΔG      ・・・(1)
Gj＜Gk−ΔG      ・・・(2)
【００３９】
図4(a)に、画素対(Pi、Pk)と注目画素Pjの例を示す。
図4(a)に示すように画素対(Pi、Pk)を設定した場合、図4(b)に示すように、画素対(Pi、Pk)の間にある注目画素Pjの画素値Gjが、周囲(画素対)の画素値よりわずかに小さい値Gi−ΔG、Gk−ΔGより小さい場合に、注目画素Pjを第3領域(ソフトプラーク領域)Cとして抽出する。
【００４０】
ΔGの値は、一般的な造影CT画像の場合は約20〜30程度とすると好適な結果が得られる。なお、この値は一例である。画像の種類や撮影の仕方、或いは表示装置107への出力レベルの設定に応じて最適な値を設定することが好ましい。
【００４１】
ΔGの値は、操作者により設定可能としてもよい。表示部28は、例えば、図5に示す操作画面3を表示装置107に表示する。図5の操作画面3には、処理対象とする画像(断層像)を表示する画像表示エリア31、微小量ΔGの値を調整するための入力欄32、領域抽出処理の開始を指示する開始ボタン33、領域抽出処理の終了を指示する終了ボタン34等が設けられている。
【００４２】
第3領域抽出部は、差分領域内の各画素対について、上述の条件式(1)、(2)による画素値の比較を繰り返す。
【００４３】
図6(a)は注目画素Pjが第3領域(ソフトプラーク領域)Cとして抽出されない例を示す。差分領域Rに設定した画素対間の注目画素Pjの画素値GjがGj＞Gi−ΔGまたはGj＞Gk−ΔGとなる場合には、注目画素Pjはソフトプラーク(第3領域C)として抽出されない。
【００４４】
図6(b)、(c)は注目画素Pjが第3領域(ソフトプラーク領域)Cとして抽出される例である。差分領域Rに設定した画素対間の注目画素Pjの画素値Gjが、
Gj＜Gi−ΔGまたはGj＜Gk−ΔGとなる場合、及び
Gj＜Gi−ΔGかつGj＜Gk−ΔGとなる場合
は、注目画素Pjはソフトプラーク(第3領域C)として抽出される。
【００４５】
注目画素Pjの画素値を周囲(画素対)の画素値と比較することにより、局所的な微小な画素値の変化を相対的にとらえることができる。その結果、ソフトプラーク領域Cと周囲の血液とを細密に峻別して抽出できる。また、差分領域内に画素対を設定して走査することで、対象を漏れなく探索できる。
【００４６】
なお、上述の条件式(1)、(2)は、造影CT画像を入力画像とした場合の例である。造影CT画像では、ソフトプラークの画素値が血液の領域の画素値より小さく表れるため、上述の条件式(1)、(2)を適用する。
【００４７】
これに対し、MR画像の場合は、撮影の仕方に応じて血液の領域の画素値とソフトプラークの画素値との大きさの関係が反転することがある。ソフトプラークの画素値が、血液の領域の画素値より大きく表れる画像では、上述の条件式(1)、(2)に代えて、以下の条件式(3)、(4)を適用する。
【００４８】
Gj＞Gi＋ΔG      ・・・(3)
Gj＞Gk＋ΔG      ・・・(4)
【００４９】
また、図4(b)のグラフは、縦軸(画素値軸)のプラス方向が逆向きになる。また、第2領域抽出部23では、閾値(平均濃度値)より大きい画素値を持つ領域を第2領域Bとして抽出する。
【００５０】
図2の対象領域設定部26は、第2領域Bと第3領域Cとを合わせた領域を対象領域として得る。対象領域はソフトプラーク領域である。図7は、血管領域A内の対象領域(B、C)の分布を模式的に示した図である。
【００５１】
なお、上述の第1領域抽出部22、第2領域抽出部24、第3領域抽出部25は抽出した各画素について、各領域A、B、Cを区別する符号をメモリ内に書き込むが、各領域抽出部22、24、25がそれぞれ別のメモリを持つものとしてもよいし、ひとつのメモリに順次上書きしてもよい。上書きする場合は、例えば、第2領域Bとして抽出されていた画素が、第3領域抽出部25の処理により第3領域Cとして上書きされることがある。最終的には第2領域B及び第3領域Cを合わせた領域が対象領域として設定される。
【００５２】
画像生成部27は、抽出された対象領域をその周囲の2次元画像上に描画して表示部28に出力する。表示部28は、対象領域が描画された2次元画像を表示装置107に表示する。図8の2次元画像41は、被検体の頸部の断層像411に対象領域41pを色づけ表示した例である。
【００５３】
また、画像生成部27は合成画像生成部271を備えるようにしてよい。合成画像生成部271は、まず3次元原画像データに基づいて、対象領域の周囲の参照用2次元画像421を生成する。図9に示す参照用2次元画像421は、被検体の頸部を左右方向に投影した骨抜きMIP画像である。更に、合成画像生成部271は、対象領域(ソフトプラーク領域)の陰影付け画像42pを生成する。そして、参照用2次元画像421上に対象領域の陰影付け画像42pを重畳または合成した合成画像42を生成する。
【００５４】
図10の参照用2次元画像431は、被検体の頸部を前後方向にMIP(Maximum Intensity Projection)処理したMIP画像である。図10の合成画像43は、参照用2次元画像431上に対象領域の陰影付け画像43pが重畳または合成されている。
【００５５】
なお、参照用2次元画像421、431は、MIP画像や骨抜きMIP画像に限定されない。MPR(Multi-Planar Reconstruction)画像、CPR(Curved Planar Reconstruction)画像、3D画像、ray sum画像、MinIP画像等、3次元原画像から生成可能なその他の種類の2次元画像としてもよい。
【００５６】
また画像生成部27は、割合表示部272を備えるようにしてもよい。割合表示部272は、第1領域(血管領域A)に対して対象領域(ソフトプラーク領域B、C)が占める面積の割合を断面毎に求め、グラフ表示する。図9の割合表示42aや図10の割合表示43aは縦軸を断面位置、横軸を％とする。割合表示42a、43aの断面位置は合成画像42、43の断面位置と一致する。
【００５７】
また画像生成部27は、断面画像表示部273を備えるようにしてもよい。断面画像表示部273は、図10に示すように、合成画像43の任意の位置がマウスポインタ43ｍ等で指定されると、指定された位置における第1領域(血管領域)の芯線に直交する断面画像43cを生成し、表示する。更に、断面画像表示部273は、断面画像43c内に抽出した対象領域を色づけ表示する。
【００５８】
また、画像生成部27はソフトプラーク形状評価部274及び表示優先度設定部275を備えてもよい。ソフトプラーク形状評価部274及び表示優先度設定部275については、第2の実施の形態で説明する。
【００５９】
次に、第1の実施の形態の画像処理装置100が実行する領域抽出処理の流れを図11のフローチャートを参照して説明する。
【００６０】
CPU101(第1領域抽出部22)は、記憶装置103または通信I/F104を介して接続される画像データベース111から処理対象とするCT画像等の画像データを入力する(ステップS101)。そして、入力画像から血管領域A(第1領域)を抽出する(ステップS102)。血管領域Aの抽出は公知の手法により行えばよい。血管を抽出するための閾値は入力装置109等により外部から入力するものとしてもよいし、予め設定されている値としてもよい。閾値処理により、例えば、図3(a)に示す領域Aが抽出される。
【００６１】
なお、閾値処理で得られる二値画像には複数の血管が含まれることがある。その場合は、CPU101は複数の血管の中から特定の血管を抽出する。例えば、領域の円形度や所在位置等を指定して複数の血管の中から特定の血管を抽出したり、マウス等により指示した領域を特定血管として抽出する。
【００６２】
次に、CPU101(閾値算出部23)はステップS102で抽出した血管領域A内の画素値に基づいて、ソフトプラーク領域を抽出するための閾値を求める。例えば、血管領域A内の平均CT値を求める(ステップS103)。CPU101(第2領域抽出部24)は、平均CT値より小さい画素値となる領域を第2領域(ソフトプラーク領域)Bとして抽出する(ステップS104)。例えば、図3(b)の斜線で示す第2領域Bが抽出される。
【００６３】
ステップS104の閾値処理ですべてのソフトプラークを抽出することは困難である。
そこで、次のステップでは、抽出し損なったソフトプラークを局所的な濃度値変化を見つけることにより抽出する。
【００６４】
CPU101(第3領域抽出部25)は、血管領域Aと抽出済みのソフトプラーク領域Bとの差分領域Rから周囲のCT値よりもわずかに低いCT値を示す画素を探索する。CPU101は、まず血管領域AとステップS104で抽出したソフトプラーク領域Bとの差分領域Rに画素対を設定する(ステップS105)。そしてCPU101は、画素対(Pi、Pk)の間にある画素Pjに注目し、注目画素PjのCT値Gjが、以下の条件(1)、(2)を満たすか否かを判定する(CT画像の場合)。以下の条件(1)、(2)のうち少なくともいずれか一方の条件を満たす場合は、注目画素Pjをソフトプラーク領域Cとして抽出する(ステップS106、ステップS107)。
【００６５】
Gj＜Gi−ΔG    ・・・(1)
Gj＜Gk−ΔG    ・・・(2)
CPU101は、差分領域内に設定する画素対の位置をずらしながら、ステップS106〜ステップS107の処理を繰り返し行う(ステップS108；Yes→ステップS106)。
【００６６】
ステップS106、ステップS107によるソフトプラーク領域Cの抽出結果を図7に示す。
【００６７】
ステップS104で抽出されたソフトプラーク領域B(血管の平均CT値よりCT値が小さい領域)の一部は、領域Cとして抽出される。図7において、領域B及び領域Cを含む領域がソフトプラーク領域(対象領域)である(ステップS109)。
【００６８】
CPU101は、断層像や血管断面画像、或いは抽出した対象領域の周囲の参照用画像を生成する。更に、抽出した対象領域を色づけしたり陰影付けしたりして参照用画像上に明確に表示する(ステップS110、ステップS111)。図8は、断層像41上に、抽出したソフトプラーク領域41pを色づけ表示した例である。
【００６９】
また例えば図9、図10に示すように、3次元原画像に基づいて骨抜きMIP画像421やMIP画像431等の参照用2次元画像を生成し、抽出した対象領域(ソフトプラーク領域)については陰影付け2次元画像42pを生成し、これらを重畳または合成して表示してもよい。
【００７０】
また、図10に示すように、マウスポインタ43mで指定した位置での血管断面画像43cを同一表示画面内に表示してもよい。血管断面画像43cにはソフトプラーク領域が明確となるように識別表示(例えば、色づけ表示)される。
【００７１】
また、図9、図10に示すように、CPU101は、血管断面におけるソフトプラークの面積の割合を表示してもよい。各断面位置での血管中ソフトプラークの面積の割合が、横軸を％、縦軸を断面位置とする座標領域内にプロットされる。これにより、操作者はソフトプラークの分布を画像上で観察したり、定量的に観察したりできる。
【００７２】
以上説明したように、第1の実施の形態の画像処理装置100によれば、CPU101は、画像から抽出した血管領域Aから、血管領域Aの平均濃度値よりCT値が小さい領域をソフトプラーク領域Bとして抽出する。また、未抽出のソフトプラークについては、血管領域Aとソフトプラーク領域Bとの差分領域Rに画素対を設定して、画素対間の各画素Pjについて、画素対のCT値よりわずかに小さい値(微小値)よりさらに小さい画素値であるか否かを判定する。これにより、CPU101は画素値の局所的な微小変化を抽出することが可能となる。よって、血管領域からの分離抽出が困難であったソフトプラーク領域を細密に抽出できる。
【００７３】
［第2の実施の形態］
次に、抽出したソフトプラークの形状評価と表示態様について説明する。
【００７４】
まず図12を参照して、状態の異なる2種類のソフトプラークについて説明する。
【００７５】
血管5内のソフトプラークには、断面α内で血管壁51に接していないソフトプラークと、接しているソフトプラークとがある。以下、血管壁51に接しているソフトプラークを接触プラーク52、血管壁51に接していないソフトプラークを非接触プラーク53と呼ぶ。
【００７６】
非接触プラーク53は、着目断面αにおいて周囲がすべて血液領域となっているソフトプラーク領域である。この種のソフトプラークは剥がれたり折れたりしやすく、病理上のリスクが大きい。
【００７７】
接触プラーク52は、血管壁51に接着しているので比較的安定であり、剥がれるリスクが小さい。
【００７８】
第2の実施の形態では、画像処理装置100の画像生成部27は、ソフトプラークの形状を評価して、病理上のリスクが大きい非接触プラーク53を優先的に描画するように画像生成を行う。
【００７９】
第2の実施の形態の画像処理装置100の画像生成部27は、図2に示すソフトプラーク形状評価部274と表示優先度設定部275の機能を用いる。なお、その他の機能構成(入力部21、第1領域抽出部22、閾値算出部23、第2領域抽出部24、第3領域抽出部25、対象領域設定部26)は第1の実施の形態の画像処理装置100と同様であるため、重複する説明を省略する。また、画像処理装置100のハードウエア構成も第1の実施の形態と同様である。
【００８０】
ソフトプラーク形状評価部274は、第1の実施の形態の領域抽出処理により抽出した対象領域が第1領域Aの内周部と接しているか否かを判定する。すなわち、血管内のソフトプラーク領域が血管壁51に接しているか否かを判定する。
【００８１】
表示優先度設定部275は、ソフトプラーク形状評価部274の評価結果に基づいて表示優先度を設定する。対象領域(ソフトプラーク領域)が第1領域Aの内周部(血管壁51)と接していない場合は、対象領域(ソフトプラーク領域)が優先的に描画されるように表示優先度が設定される。例えば、通常のMIP画像等では、投影線上で最大画素値が投影されるが、第2の実施の形態では、非接触プラーク53の表示を優先する。更に、表示優先度設定部275は、視点と対象領域(ソフトプラーク領域)との位置関係を考慮して表示優先度を設定してもよい。例えば、投影線上に非接触プラーク53が2か所ある場合は、視点に近い方の非接触プラークの表示を優先する。
【００８２】
次に、図13〜図20を参照して、第2の実施の形態の画像処理装置100における画像生成処理について説明する。
【００８３】
なお、図13の画像生成処理の対象とする画像は、例えば、第1の実施の形態の領域抽出処理(図11のステップS101〜ステップS109)と同様の手順で血管領域やソフトプラーク領域が抽出されているものとする。
【００８４】
CPU101は、ソフトプラーク領域の形状評価処理を実行する(ステップS201)。形状評価処理は、例えば図14に示す手順(形状評価処理(1))に従って行うものとする。
【００８５】
図14に示す形状評価処理(1)において、CPU101(ソフトプラーク形状評価部274)は、まず血管5に直交する断面αを設定する(ステップS301)。CPU101は、断面α内のソフトプラーク領域が血管壁51に接するか否かを判定する(ステップS302)。血管壁51に接しているか否かは、ソフトプラーク領域の周囲の画素値に基づいて判定する。周囲に血液の画素値でない画素がある場合は、当該ソフトプラーク領域が接触プラーク52であるという結果を返す(ステップS302；Yes→ステップS303)。ソフトプラーク領域の周囲の画素がすべて血液の画素値である場合は、非接触プラーク53であるという結果を返す(ステップS302；No→ステップS304)。
【００８６】
図13の画像生成処理のステップS202に戻る。ステップS202において、CPU101はソフトプラークの形状評価結果に基づいて表示優先度設定処理を行う。表示優先度設定処理は、例えば、図15に示す手順で行うものとする。
【００８７】
なお、図15の表示優先度設定処理の開始に際し、投影面と視点位置とが設定されているものとする。
【００８８】
CPU101は、投影線69上に非接触プラーク53が存在するか否かを判定する(ステップS401)。投影線69上に非接触プラーク53が存在する場合は(ステップS401；Yes)、更に、投影線69上に非接触プラーク53と接触プラーク52とが混在するか否かを判定する(ステップS402)。投影線69上に非接触プラーク53と接触プラーク52とが混在する場合は(ステップS403；Yes)、非接触プラーク53の表示を優先する(ステップS403)。また、非接触プラーク53のみの場合は(ステップS402；非接触プラークのみ)、視点に近い非接触プラーク53の表示を優先する(ステップS404)。
【００８９】
投影線69上に非接触プラーク53がなく、接触プラーク52のみが存在する場合は(ステップS401；接触プラークのみ)、視点に近いソフトプラークの表示を優先する(ステップS404)。
【００９０】
このような手順で各血管直交断面αについて表示優先度設定処理を行うと、図13のステップS203の処理へ戻る。CPU101は、設定された表示優先度に基づいて2次元画像を生成する(ステップS203)。
【００９１】
図16〜図19は、図12に示す血管5の直交断面α(α1〜α4)におけるソフトプラークの配置と投影像との関係を示す。投影面60は、血管直交断面α(α1〜α4)に対して、例えば垂直な面とする。視点は、投影面60に対して血管5を介して対向する位置に配置する。血管5を投影面60に平行投影する例を示す。投影方法はMIPでもよいし、3Dでもよいし、ray-sumでもよいし、その他の方法でもよい。
【００９２】
血管壁51に接触していない非接触プラーク53の表示優先度は他の領域よりも高く設定されている。CPU101は、投影面60に対してソフトプラークを所定の投影方向から投影して投影画像を作成する。
【００９３】
図16に示す断面α1には非接触プラーク53と接触プラーク52とが混在する。この場合、非接触プラーク53を接触プラーク52より優先的に投影面60に投影する。また、非接触プラーク53の投影像63は、病理上リスクの高さを示すため赤系の表示色を使用することが好ましい。具体的には、非接触プラーク53に対応する投影像63を赤系、接触プラーク52に対応する投影像62を緑系で色分け表示する。
【００９４】
通常のMIP処理では投影線上の最大画素値が投影されるため、ソフトプラークより造影血液が描画され、ソフトプラークの存在を投影像上は確認できなかったが、本実施の形態の画像生成処理により、非接触プラーク53が描画されるようになる。操作者は非接触プラーク53の存在をMIP画像上で確認できる。また色分け表示により、病理上のリスクを直感的に察知しやすくなる。
【００９５】
図17に示す断面α2には非接触プラーク53が存在せず、複数の接触プラーク52a、52bが存在する。この場合は、視点に近い方の接触プラーク52aの表示優先度が高く設定される。具体的には、CPU101は、視点に近い接触プラーク52aの投影像62aをもう一つの接触プラーク52bに対応する投影像62bより前面となるように投影する。また、各投影像62a、62bの色は、同じ緑系でも彩度や明度が異なる配色とすれば、複数の領域が重なって投影されていることを表現できる。
【００９６】
図18に示す断面α3には非接触プラーク53と接触プラーク52とが存在する。また投影線69上で接触プラーク52が非接触プラーク53よりも視点寄りの位置に存在する。この場合は、接触プラーク52の投影像62内に非接触プラーク53の投影像63を半透明カラーや破線等で表示する。これにより、視点位置からは隠れているが、投影線上に非接触プラーク53が存在していることを示すことができる。
【００９７】
図19に示す断面α4には非接触プラーク53が存在しない。接触プラーク52a、52bが2か所に存在する。この場合は、図17の例と同様に、視点に近い方の接触プラーク52aの表示優先度を、他の接触プラーク52bより高くする。また、視点に近い接触プラーク52aの投影像62aを他方の接触プラーク52bに対応する投影像62bより前面となるように投影する。各投影像62a、62bの色は、同じ緑系でも彩度や明度が異なる配色とすれば、複数の領域が存在し、重なっていることを表現できる。
【００９８】
図20は中心投影の例を示す。中心投影では、CPU101は血管内部に視点60bを設定し、血管の周囲に血管形状に沿った曲円筒状の投影面60aを設定する。更に、曲円筒状投影面60aの画像を平面状投影面60に展開する。
【００９９】
中心投影の場合も図16と同様に、非接触プラーク53に対応する投影像63を赤系で描画し、接触プラーク52に対応する投影像62を緑系で描画する。これにより、非接触プラーク53の存在とその危険度が視認しやすくなる。
【０１００】
投影線69上に接触プラーク52と非接触プラーク53とが混在する場合は、非接触プラーク53の表示優先度が他より高く設定される。
【０１０１】
図21に第2の実施の形態において生成された画像の例を示す。
図21(a)の画像45は、原CT画像(3次元原画像)を用いて作成したMIP画像451に、ある特定の血管45b内のソフトプラーク領域45cを陰影付けして合成表示している。ソフトプラーク領域45cのうち一部の領域45dは他領域とは異なる表示色または表示形式で表示される。この領域45dが非接触プラーク53を示す。実際の画像では、例えば赤系の表示色で描画することが好ましい。他のソフトプラーク領域45cは、接触プラーク52であり、実際の画像では、緑系で描画することが好ましい。
【０１０２】
図21(b)の画像46は、原CT画像(3次元原画像)を用いて作成した骨抜きMIP画像461に、ある特定の血管46b内のソフトプラーク領域46cを陰影付けして合成表示している。ソフトプラーク領域46cのうち一部の領域46dは他領域とは異なる表示色または表示形式で表示される。この領域46cが非接触プラーク53を示す。実際の画像では、例えば赤系で描画することが好ましい。他のソフトプラーク領域46cは、接触プラーク52であり、実際の画像では、緑系で描画することが好ましい。骨抜きMIP画像461では動脈血管の石灰化領域46e等も描画される。画像46により、石灰化領域46eとソフトプラーク46c、46dとを同時に観察できる。
【０１０３】
以上説明したように、第2の実施の形態の画像処理装置100は、ソフトプラーク領域の形状を判定し、非接触プラーク53については表示優先度を高くして2次元画像を生成する。また、非接触プラーク53を赤系、接触プラーク52を緑系の配色で描画するなど、表示色を変える。これにより、リスクが大きい非接触プラーク53が存在することを明確かつ直感的に示すことができる。
【０１０４】
［第3の実施の形態］
第3の実施の形態では、第2の実施の形態のソフトプラークの形状評価処理(2)(図14参照)に代えて、図25に示すソフトプラーク形状評価処理(2)を行う。
【０１０５】
第3の実施の形態のソフトプラークの形状評価処理(2)では、非接触プラーク53の3次元形状を分析する。具体的には、ある断面で非接触プラーク53と判定されたソフトプラークが別の断面で、血管壁や別のソフトプラークに接着しているか否かを判定する。
【０１０６】
図22〜図24を参照して、ソフトプラークの形状について説明する。
図22に示す非接触プラーク53Aは、一端が接触プラーク52に到達し付着するが、他端側では非接触である。以下、このような形状のソフトプラークを片側固定非接触プラーク53Aと呼ぶ。
【０１０７】
図23に示す非接触プラーク53Bは、両端が接触プラーク52に到達し付着する。以下、このような形状のソフトプラークを両側固定非接触プラーク53Bと呼ぶ。図24に示すソフトプラークのように、接触プラーク52a、52bが2か所に分離して存在し、非接触プラークは一端が接触プラーク52aに付着し、他端側が接触プラーク52bに付着する場合もある。このような形状も両側固定非接触プラーク53Bと呼ぶ。
【０１０８】
ソフトプラークの3次元的な形状を解析するために、第3の実施の形態では図25に示す形状評価処理(2)を実行する。
【０１０９】
以下、形状評価処理(2)について説明する。
まず、CPU101は血管直交断面α内でソフトプラークが血管壁51に接触か否かを判定する(ステップS501)。接触する場合は(ステップS501；Yes)、接触プラーク52と判定する(ステップS502)。血管直交断面α内でソフトプラークが血管壁51に接触しない場合は(ステップS501；No)、ステップS503へ移行する。
【０１１０】
CPU101は、着目しているソフトプラーク内の任意の画素を始点7として、異方性3次元リージョングローイング処理を行う。異方性3次元リージョングローイング処理により領域拡張し、全ての方向(少なくとも2方向以上)で血管壁または接触プラーク52に接触するか否かを判定する(ステップS503)。
【０１１１】
例えば、図22の非接触ソフトプラーク53Aのように、始点7からある方向に領域拡張処理すると非接触のまま拡張処理が終了するが、始点7から別の方向に領域拡張処理すると接触プラークに到達する。このように、ある断面では非接触ソフトプラークであるが3次元でみると片方固定(接触)で片方非接触となる場合は、片側固定非接触プラーク53Aであるという結果を得る(ステップS505)。
【０１１２】
また、図23、図24に示すように、始点7から反対2方向に領域拡張処理した場合に、2方向とも接触プラーク52に到達する場合は(ステップS503；Yes)、両側固定非接触プラークであるという結果を得る(ステップS504)。図23と図24の違いは、図23では非接触プラーク53が一塊の接触プラーク52に固定されるが、図24では、非接触プラーク53が別の塊の接触プラーク52a、52bにそれぞれ固定されている。
【０１１３】
図23、図24に示すように、複数の方向に異方性3次元リージョングローイング処理を行った場合に、少なくとも2方向以上で接触プラーク52に到達すれば、非接触プラークであっても剥がれ落ちるリスクが小さい。そのため、両側固定非接触プラーク53Bについては、CPU101は、例えば、黄系の表示色で描画する。
【０１１４】
一方、複数の方向に異方性3次元リージョングローイングした場合に、いずれかの方向で接触プラーク52に到達しないまま拡張処理が終わった場合(図23の片側固定非接触プラーク53A)、剥がれ落ちるリスクが最も大きい。そのため片側固定非接触プラーク53Aについては、例えば赤系の表示色で描画する。
【０１１５】
以上説明したように、第3の実施の形態では、非接触プラーク53の3次元形状を解析し、端部が接触プラークに固定されているか否かを判定する。そして、片側で固定されている場合と両側で固定されている場合とで、異なる表示カラー等で描画する。ソフトプラークが剥がれ落ちるリスクが大きい形状であるか否かを評価して、その形状に応じて色分け表示できる。
【０１１６】
以上、添付図面を参照しながら、本発明に係る画像処理装置等の好適な実施形態について説明したが、本発明はかかる例に限定されない。例えば、上述の実施の形態では、血管からソフトプラークを抽出する処理について説明したが、その他の臓器の画像を対象としてもよい。当業者であれば、本願で開示した技術的思想の範疇内において、各種の変更例又は修正例に想到し得ることは明らかであり、それらについても当然に本発明の技術的範囲に属するものと了解される。
【符号の説明】
【０１１７】
1  画像処理システム、100  画像処理装置、101  CPU101、102  主メモリ、103  記憶装置、104  通信I/F、105  表示メモリ、106  I/F、107  表示装置、108  マウス、109  入力装置、110  ネットワーク、111  画像データベース、112  医用画像撮影装置、113  バス、21  入力部、22  第1領域抽出部、23  閾値算出部、24  第2領域抽出部、25  第3領域抽出部、251  差分領域算出部、252  画素対設定部、253  比較部、26  対象領域設定部、27  画像生成部、271  合成画像生成部、272  割合表示部、273  断面画像表示部、274  ソフトプラーク形状評価部、275  表示優先度設定部、28  表示部、A  血管領域(第1領域)、B  ソフトプラーク領域(第2領域)、C  ソフトプラーク領域(第3領域)、52  接触プラーク、53  非接触プラーク、53A  片側固定非接触プラーク、53B  両側固定非接触プラーク、60  投影面
【図１】
2014132829.tif 000003
【図２】
2014132829.tif 000004
【図３】
2014132829.tif 000005
【図４】
2014132829.tif 000006
【図５】
2014132829.tif 000007
【図６】
2014132829.tif 000008
【図７】
2014132829.tif 000009
【図８】
2014132829.tif 000010
【図９】
2014132829.tif 000011
【図１０】
2014132829.tif 000012
【図１１】
2014132829.tif 000013
【図１２】
2014132829.tif 000014
【図１３】
2014132829.tif 000015
【図１４】
2014132829.tif 000016
【図１５】
2014132829.tif 000017
【図１６】
2014132829.tif 000018
【図１７】
2014132829.tif 000019
【図１８】
2014132829.tif 000020
【図１９】
2014132829.tif 000021
【図２０】
2014132829.tif 000022
【図２１】
2014132829.tif 000023
【図２２】
2014132829.tif 000024
【図２３】
2014132829.tif 000025
【図２４】
2014132829.tif 000026
【図２５】
2014132829.tif 000027
【手続補正書】
【提出日】平成27年12月21日(2015.12.21)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
画像データを入力する入力部と、
前記入力部により入力された画像データから第1領域を抽出する第1領域抽出部と、
前記第1領域内の濃度値に基づく閾値により前記第1領域を閾値処理し、前記第1領域内から第2領域を抽出する第2領域抽出部と、
前記第1領域と前記第2領域との差分領域内にある2画素の組み合わせである画素対を設定し、各画素対の間にある画素を注目画素とし、前記注目画素の画素値と前記画素対の少なくともいずれか一方の画素値との差が所定の値より大きい場合に、前記注目画素を第3領域として抽出する第3領域抽出部と、
前記第2領域と前記第3領域とを合わせた領域を対象領域として得る対象領域設定部と、
を備えることを特徴とする画像処理装置。
【請求項２】
前記画像データは3次元原画像であり、
前記3次元原画像に基づいて前記対象領域とその周囲を含む範囲の2次元画像を生成する画像生成部と、
前記画像生成部により生成された2次元画像を表示する表示部と、
を更に備えることを特徴とする請求項1に記載の画像処理装置。
【請求項３】
前記第1領域は血管領域であり、前記対象領域はソフトプラーク領域であることを特徴とする請求項1に記載の画像処理装置。
【請求項４】
前記画像生成部は、
前記対象領域の周囲について参照用2次元画像を生成するとともに前記対象領域の陰影付け画像を生成し、前記参照用2次元画像上に前記対象領域の陰影付け画像を重畳または合成した合成画像を生成することを特徴とする請求項2に記載の画像処理装置。
【請求項５】
前記第1領域に対する前記対象領域の占める割合を表示する割合表示部を更に備えることを特徴とする請求項1に記載の画像処理装置。
【請求項６】
前記合成画像の任意の位置を指定する位置指定部と、
前記位置指定部により指定された位置での前記第1領域に直交する断面画像に前記対象領域を重畳表示する断面画像表示部と、
を備えることを特徴とする請求項4に記載の画像処理装置。
【請求項７】
前記対象領域の形状を評価する形状評価部を更に備え、
前記画像生成部は、前記形状評価部による評価の結果、前記対象領域が前記第1領域の内周部と接していない場合に前記対象領域が優先的に表示されるように表示優先度を設定し、前記3次元原画像に基づく2次元画像を生成することを特徴とする請求項2に記載の画像処理装置。
【請求項８】
前記画像生成部は、更に、視点と前記対象領域との位置関係に基づいて前記表示優先度を設定することを特徴とする請求項7に記載の画像処理装置。
【請求項９】
画像から血管領域を抽出する血管領域抽出部と、
前記血管領域内のソフトプラーク領域を抽出するソフトプラーク領域抽出部と、
前記ソフトプラーク領域の形状を評価する形状評価部と、
前記ソフトプラーク領域の形状評価結果に基づいてソフトプラーク領域の表示優先度を設定する表示優先度設定部と、
前記表示優先度に基づいて前記血管領域の2次元画像を生成する画像生成部と、
を備えることを特徴とする画像処理装置。
【請求項１０】
前記形状評価部は、前記血管領域の断面画像に基づいて前記ソフトプラーク領域が血管壁に接しているか否かを判定することを特徴とする請求項9に記載の画像処理装置。
【請求項１１】
前記形状評価部は、更に、前記血管領域の3次元原画像に基づいて前記ソフトプラーク領域の周囲の各方向を3次元に探索し、前記ソフトプラーク領域が前記血管壁に延伸し、到達するか否かを判定することを特徴とする請求項10に記載の画像処理装置。
【請求項１２】
前記表示優先度設定部は、
前記ソフトプラーク領域が前記血管壁と接していない場合は、接している場合より高い示優先度を設定することを特徴とする請求項9記載の画像処理装置。
【請求項１３】
コンピュータを用いて画像データから対象領域を抽出する領域抽出方法であって、
画像データを入力するステップと、
入力された画像データを閾値処理して第1領域を抽出するステップと、
前記第1領域内の濃度値に基づく閾値により前記第1領域を閾値処理し、前記第1領域内から第2領域を抽出するステップと、
前記第1領域と前記第2領域との差分領域内にある2画素の組み合わせである画素対を設定し、各画素対の間にある画素を注目画素とし、前記注目画素の画素値と前記画素対の少なくともいずれか一方の画素値との差が所定の値より大きい場合に、前記注目画素を第3領域として抽出するステップと、
前記第2領域と前記第3領域とを合わせた領域を対象領域として得るステップと、
を含むことを特徴とする領域抽出方法。
【手続補正２】
【補正対象書類名】明細書
【補正対象項目名】０００１
【補正方法】変更
【補正の内容】
【０００１】
本発明は、画像処理装置、及び領域抽出方法に係り、詳細には、CT画像またはMR画像における血管ソフトプラークの抽出及び描画に関する。
【手続補正３】
【補正対象書類名】明細書
【補正対象項目名】００１２
【補正方法】変更
【補正の内容】
【００１２】
【図１】画像処理装置100の全体構成を示す図
【図２】画像処理装置100の機能構成を示すブロック図
【図３】(a)血管領域A、(b)ソフトプラーク領域Bの抽出について説明する図
【図４】(a)画素対(Pi、Pk)の設定例、(b)局所的な画素値の微小な変化を表すグラフ
【図５】操作画面3の一例
【図６】画素対を用いた領域抽出条件判定について説明する図
【図７】対象領域(ソフトプラーク領域B、C)の分布例
【図８】抽出した対象領域41pを断層像411上に重畳表示した画像41
【図９】抽出した対象領域42pを骨抜きMIP画像421に重畳または合成表示した2次元画像42及びソフトプラークの面積割合表示42aの例
【図１０】抽出した対象領域43pをMIP画像431に重畳または合成表示した2次元画像43、断面表示43c、及びソフトプラークの面積割合表示43aの例
【図１１】本発明に係る領域抽出処理の手順を説明するフローチャート
【図１２】血管5内のソフトプラークの形状と種類について説明する図
【図１３】画像生成処理の手順を説明するフローチャート
【図１４】形状評価処理(1)の手順を説明するフローチャート
【図１５】表示優先度設定処理の手順を説明するフローチャート
【図１６】ソフトプラークの形状に基づいて設定された表示優先度に応じた描画について説明する図(非接触プラーク53と接触プラーク52とが存在する場合)
【図１７】ソフトプラークの形状に基づいて設定された表示優先度に応じた描画について説明する図(複数の接触プラーク52a、52bが存在する場合)
【図１８】ソフトプラークの形状に基づいて設定された表示優先度に応じた描画について説明する図(非接触プラーク53と接触プラーク52とが存在する場合
【図１９】ソフトプラークの形状に基づいて設定された表示優先度に応じた描画について説明する図(複数の接触プラーク52a、52bが存在する場合)
【図２０】視点70を血管内に設定した描画例
【図２１】ソフトプラークの形状に応じた色分け表示例。(a)MIP画像451に合成表示(b)骨抜きMIP画像461合成表示
【図２２】片側固定非接触プラーク53Aについて説明する図
【図２３】両側固定非接触プラーク53Bについて説明する図
【図２４】両側固定非接触プラーク53Bの別の例について説明する図
【図２５】形状評価処理(2)の手順を説明するフローチャート
【手続補正４】
【補正対象書類名】明細書
【補正対象項目名】００３０
【補正方法】変更
【補正の内容】
【００３０】
入力部21は、画像データを入力する。画像データは、被検体をX線CT装置やMRI装置等を用いて撮影した複数の断層像である。これらの複数の断層像を積み上げたものを3次元原画像データと呼ぶ。以下、入力画像がCT画像である場合を例として説明する。
【手続補正５】
【補正対象書類名】明細書
【補正対象項目名】００３６
【補正方法】変更
【補正の内容】
【００３６】
差分領域算出部251は、第1領域(血管領域A)と第2領域(ソフトプラーク領域B)との差分領域を求める。画素対設定部252は、差分領域内にある2画素の組み合わせである画素対(Pi、Pk)を設定する。比較部253は、画素対設定部252により設定された各画素対(Pi、Pk)の間にある画素Pjを注目画素とし、注目画素Pjの画素値Gjと、画素対(Gi、Gk)の画素値よりわずかに小さい値Gi−ΔG、Gk−ΔGとを比較する。注目画素Pjの画素値Gjの方が小さい場合に、この注目画素をソフトプラーク領域として抽出する。
【手続補正６】
【補正対象書類名】明細書
【補正対象項目名】００４９
【補正方法】変更
【補正の内容】
【００４９】
また、図4(b)のグラフは、縦軸(画素値軸)のプラス方向が逆向きになる。また、第2領域抽出部24では、閾値(平均濃度値)より大きい画素値を持つ領域を第2領域Bとして抽出する。
【手続補正７】
【補正対象書類名】明細書
【補正対象項目名】００９８
【補正方法】変更
【補正の内容】
【００９８】
図20は中心投影の例を示す。中心投影では、CPU101は血管内部に視点70を設定し、血管の周囲に血管形状に沿った曲円筒状投影面60aを設定する。更に、曲円筒状投影面60aの画像を平面状投影面60bに展開する。
【手続補正８】
【補正対象書類名】明細書
【補正対象項目名】０１０６
【補正方法】変更
【補正の内容】
【０１０６】
図22〜図24を参照して、ソフトプラークの形状について説明する。
図22に示す非接触プラークは、一端が接触プラーク52に到達し付着するが、他端側では非接触である。以下、このような形状のソフトプラークを片側固定非接触プラーク53Aと呼ぶ。
【手続補正９】
【補正対象書類名】明細書
【補正対象項目名】０１０７
【補正方法】変更
【補正の内容】
【０１０７】
図23に示す非接触プラークは、両端が接触プラーク52に到達し付着する。以下、このような形状のソフトプラークを両側固定非接触プラーク53Bと呼ぶ。図24に示すソフトプラークのように、接触プラーク52a、52bが2か所に分離して存在し、非接触プラークは一端が接触プラーク52aに付着し、他端側が接触プラーク52bに付着する場合もある。このような形状も両側固定非接触プラーク53Bと呼ぶ。
【手続補正１０】
【補正対象書類名】明細書
【補正対象項目名】０１１１
【補正方法】変更
【補正の内容】
【０１１１】
例えば、図22のように、始点7からある方向に領域拡張処理すると非接触のまま拡張処理が終了するが、始点7から別の方向に領域拡張処理すると接触プラークに到達する。このように、ある断面では非接触ソフトプラークであるが3次元でみると片方固定(接触)で片方非接触となる場合は、片側固定非接触プラーク53Aであるという結果を得る(ステップS505)。
【手続補正１１】
【補正対象書類名】明細書
【補正対象項目名】０１１４
【補正方法】変更
【補正の内容】
【０１１４】
一方、複数の方向に異方性3次元リージョングローイング処理を行った場合に、いずれかの方向で接触プラーク52に到達しないまま拡張処理が終わった場合(図23の片側固定非接触プラーク53A)、剥がれ落ちるリスクが最も大きい。そのため片側固定非接触プラーク53Aについては、例えば赤系の表示色で描画する。
【手続補正１２】
【補正対象書類名】明細書
【補正対象項目名】０１１７
【補正方法】変更
【補正の内容】
【０１１７】
1  画像処理システム、100  画像処理装置、101  CPU、102  主メモリ、103  記憶装置、104  通信I/F、105  表示メモリ、106  I/F、107  表示装置、108  マウス、109  入力装置、110  ネットワーク、111  画像データベース、112  医用画像撮影装置、113  バス、21  入力部、22  第1領域抽出部、23  閾値算出部、24  第2領域抽出部、25  第3領域抽出部、251  差分領域算出部、252  画素対設定部、253  比較部、26  対象領域設定部、27  画像生成部、271  合成画像生成部、272  割合表示部、273  断面画像表示部、274  ソフトプラーク形状評価部、275  表示優先度設定部、28  表示部、A  血管領域(第1領域)、B  ソフトプラーク領域(第2領域)、C  ソフトプラーク領域(第3領域)、52  接触プラーク、53  非接触プラーク、53A  片側固定非接触プラーク、53B  両側固定非接触プラーク、60  投影面
【手続補正１３】
【補正対象書類名】図面
【補正対象項目名】図２０
【補正方法】変更
【補正の内容】
【図２０】
2014132829.tif 000028
【国際調査報告】
2014132829.tif 000029
2014132829.tif 000030
2014132829.tif 000031
2014132829.tif 000032
2014132829.tif 000033
2014132829.tif 000034
2014132829.tif 000035
2014132829.tif 000036