(57)【要約】  画像から特定の領域を認識する処理において、画像内に複数存在し濃度値にゆらぎのある領域を、簡易な操作で精度よく認識することが可能な画像処理装置等を提供するために、CPU101は、3次元原画像30に含まれる注目画素31とその周囲の複数の画素(判定範囲33)に対して所定の閾値条件を適用して閾値判定し、閾値条件を満たす場合に注目画素を認識画素とする。閾値条件は注目画素31と同一平面にある画素と、別の平面にある画素とで異なる閾値を適用することが望ましい。また、注目画素31(判定範囲33)を順次移動して上述の閾値判定を繰り返し実行することにより3次元原画像30全体に対して上述の閾値判定を行うものとする。これにより、操作者が開始点を設定する必要がなく、自動で軟骨のように濃度値にゆらぎがある組織をより精度よく認識することができる。
2014038428.tif 000002
【特許請求の範囲】
【請求項１】
画像から特定の領域を認識する処理を実行する画像処理装置であって、
前記画像に含まれる画素のうち、注目画素と周囲の複数の画素を含む判定範囲に対して所定の閾値条件を適用して閾値判定し、前記判定範囲内のいずれかの画素で前記閾値条件を満たす場合に前記注目画素を認識画素とする閾値判定部と、
前記判定範囲を順次移動して前記閾値判定を繰り返し実行する処理実行部と、
を備えることを特徴とする画像処理装置。
【請求項２】
前記閾値判定部は、前記注目画素と同一平面にある画素と、別の平面にある画素とで異なる閾値条件を適用し、前記判定範囲内のいずれかの画素で前記閾値条件を満たすと判定した場合に、前記注目画素を前記認識画素とすることを特徴とする請求項1に記載の画像処理装置。
【請求項３】
前記閾値条件として、前記注目画素からの距離に応じて異なる閾値を設定することを特徴とする請求項2に記載の画像処理装置。
【請求項４】
前記閾値判定部は、前記判定範囲内のすべての画素に対し、同一の閾値条件を適用し、いずれかの画素で前記閾値条件を満たすと判定した場合に、前記注目画素を前記認識画素とすることを特徴とする請求項1に記載の画像処理装置。
【請求項５】
前記画像の濃度値分布のヒストグラムに基づいて前記閾値条件を算出する閾値自動算出部を更に備え、
前記閾値判定部は、前記閾値自動算出部により算出された閾値条件を適用して前記閾値判定を行うことを特徴とする請求項1に記載の画像処理装置。
【請求項６】
前記閾値自動算出部は、
前記画像の濃度値のヒストグラムの最大のピーク位置を求め、前記ピーク位置より濃度値の低い部分のデータを使用して前記ピーク位置を頂点として左右対称な形状の分布曲線を求め、前記ヒストグラムの分布曲線から前記左右対称な形状の分布曲線を減算することにより減算後分布曲線を求め、減算後分布曲線の特徴に基づいて前記閾値条件を求めることを特徴とする請求項5に記載の画像処理装置。
【請求項７】
前記閾値条件を操作者に設定させるための閾値手動設定部を更に備え、
前記閾値判定部は、前記閾値手動設定部により設定された閾値条件を適用して前記閾値判定を行うことを特徴とする請求項1に記載の画像処理装置。
【請求項８】
前記閾値判定における処理モードとして、前記注目画素単体について閾値判定を行うか、前記判定範囲に含まれる全ての画素を同一の閾値条件を用いて閾値判定するか、前記注目画素と同一平面にある画素と、別の平面にある画素とで異なる閾値条件を用いて閾値判定するか、を選択する処理モード選択部を更に備えることを特徴とする請求項1に記載の画像処理装置。
【請求項９】
前記閾値条件を操作者に手動設定させるか、或いは前記閾値条件を前記画像の濃度値分布のヒストグラムに基づいて算出するかを選択する閾値設定方法選択部を更に備え、
前記閾値判定部は、前記閾値設定方法選択部により選択された方法で前記閾値条件を設定することを特徴とする請求項1に記載の画像処理装置。
【請求項１０】
前記閾値判定部により認識された認識画素を前記画像から除去または抽出した診断用画像を生成する診断用画像生成部を更に備えることを特徴とする請求項1に記載の画像処理装置。
【請求項１１】
コンピュータを用いて、画像から特定の領域を認識する画像処理方法であって、
前記画像に含まれる画素のうち、注目画素と周囲の複数の画素を含む判定範囲に対して所定の閾値条件を適用して閾値判定し、前記判定範囲内のいずれかの画素で前記閾値条件を満たす場合に前記注目画素を認識画素とする閾値判定処理を、前記判定範囲を順次移動して繰り返し実行することを特徴とする画像処理方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像から特定の領域を認識する画像処理装置等に関する。
【背景技術】
【０００２】
従来から、例えばX線CT(Computed Tomography)装置やMRI(Magnetic Resonance Imaging)装置等によって撮影される断層像群から3次元原画像(3次元ボリュームデータ)を生成し、3次元原画像から特定の組織を抽出または除去し、診断に適した画像(以下、診断用画像という)を生成して表示する手法が知られている。診断用画像としては、例えば、3次元ボリュームレンダリング画像やMIP(Maximum Intensity Projection)画像等が作成される。具体的には、頭部の血管走行を観察する場合には、上述の断層像群から頭蓋骨や軟骨等を含む骨領域を除去した骨抜きMIP画像を生成することがある。【０００３】
ところで、上述したように特定組織を抽出または除去した診断用画像を生成するために、コンピュータ等を用いて元となる画像から特定の臓器領域を自動抽出する方法が提案されている。その抽出手法としては、例えば、領域拡張法(リージョングローイング法)等がある。領域拡張法とは、画像上で開始点となる画素を操作者が指定すると、コンピュータがその開始点の画素またはその画素を含む周囲の画素について所定の閾値条件を満たすか否かを判定し、条件を満たす場合に、その画素を抽出領域として拡張していくものである。【０００４】
特許文献1には、領域拡張法により領域抽出を行う際に、わずかな画素連結による過剰抽出を防止する方法が開示されている。この方法では、既抽出領域に連なり、複数の画素からなるセンス領域内のうち、所定の割合の画素が条件を満たしたときにだけ領域を拡張する。また、特許文献2には、骨の内部組織である海綿骨と骨髄等のように分離しにくい組織を区別するための方法が開示されている。具体的には、海綿骨と骨髄を含むMR画像の画素値ヒストグラムを求め、これを3つの正規分布曲線(海綿骨、骨髄、両者の中間値)をもってフィッティングすることにより、両者を分離してより正確に海綿骨体積率を算出することが記載されている。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特許4538260号公報
【特許文献２】特許4487080号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかしながら、上述したような領域拡張法や組織の分離手法を用いても、軟骨のように濃度値のゆらぎ(変動)が大きい組織を画像中から正確に区別して認識してすることは困難であった。また、頭部の画像において、軟骨領域はすべてが連なっているわけでなく点在するため、領域拡張法を用いる場合には、操作者がひとつずつ処理の開始点を指定する必要があった。そのため、操作が非常に煩雑であった。【０００７】
本発明は、以上の問題点に鑑みてなされたものであり、画像から特定の領域を認識する処理において、画像内に複数存在し、濃度値のゆらぎの大きい領域を、簡易な操作で精度よく認識することが可能な画像処理装置及び画像処理方法を提供することを目的とする。
【課題を解決するための手段】
【０００８】
前述した目的を達成するために、第1の発明は、画像から特定の領域を認識する処理を実行する画像処理装置であって、前記画像に含まれる画素のうち、注目画素と周囲の複数の画素を含む判定範囲に対して所定の閾値条件を適用して閾値判定し、前記判定範囲内のいずれかの画素で前記閾値条件を満たす場合に前記注目画素を認識画素とする閾値判定部と、前記判定範囲を順次移動して前記閾値判定を繰り返し実行する処理実行部と、を備えることを特徴とする画像処理装置である。【０００９】
なお、上述の「認識」とは、該当する画素を3次元原画像内の他の画素と区別して保持することである。その具体的方法としては、例えば、(a)元の3次元原画像中に目印となる特定値を直接書き込む。または、(b)3次元原画像とは別のメモリに二値画像として認識画素を記録することである。【００１０】
第2の発明は、コンピュータを用いて、画像から特定の領域を認識する画像処理方法であって、前記画像に含まれる画素のうち、注目画素と周囲の複数の画素を含む判定範囲に対して所定の閾値条件を適用して閾値判定し、前記判定範囲内のいずれかの画素で前記閾値条件を満たす場合に前記注目画素を認識画素とする閾値判定処理を、前記判定範囲を順次移動して繰り返し実行することを特徴とする画像処理方法である。
【発明の効果】
【００１１】
本発明の画像処理装置及び画像処理方法により、画像から特定の領域を認識する処理において、画像内に複数存在し、濃度値のゆらぎの大きい領域を、簡易な操作で精度よく認識することが可能となる。例えば、軟骨領域を除去する画像処理を、より簡易かつ精度よく行えるようになる。
【図面の簡単な説明】
【００１２】
【図１】画像処理装置100の全体構成を示す図
【図２】処理対象とする3次元原画像30内の注目画素31と、注目画素31の周囲の画素であって閾値判定対象とする範囲(判定範囲33)の一例を示す図
【図３】本発明に係る画像処理装置100が実行する診断用画像生成処理の流れを説明するフローチャート
【図４】操作画面20の一例
【図５】各処理モードにおける閾値条件の適用の仕方を説明する図であり、(a)は単体モード、(b)は等方的近隣判定モード、(c)は異方的近隣判定モードについて説明する図
【図６】異方的近隣判定モードによる閾値判定処理の流れを説明するフローチャート
【図７】閾値を自動設定する場合に用いる濃度値ヒストグラムについて説明する図
【図８】(a)濃度値ヒストグラムのピーク位置から左側(CT値の低い側)のデータを使った左右対称な分布曲線72について説明する図。(b)元の濃度値ヒストグラムから左右対称な分布曲線を減算した減算後分布曲線73。
【図９】異方的近隣判定モード及び閾値自動設定を用いて閾値判定し、認識画素を除去して生成した頭部の骨抜きMIP画像
【図１０】従来の方法(画素単体モード及び閾値手動設定)を用いて骨を除去した頭部のMIP画像
【図１１】異方的近隣判定モードにおける別の閾値条件設定例
【図１２】異方的近隣判定モードにおいて、注目画素からの距離に応じて異なる閾値条件を設定する例
【発明を実施するための形態】
【００１３】
以下図面に基づいて、本発明の実施の形態を詳細に説明する。
【００１４】
［第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参照)において、処理対象とする画像から特定の領域を認識する処理を行う。本実施の形態では、処理対象とする画像は、複数の断層像を積み上げてなる3次元原画像とする。領域を認識する処理において、CPU101は、3次元原画像に含まれる画素のうち注目画素の周囲の複数の画素を判定範囲とし、判定範囲の各画素に対して所定の閾値条件を適用して閾値判定する。そして、閾値条件を満たす場合に注目画素を認識画素とする。また、CPU101は、3次元原画像全体を走査して上述の閾値判定を行う。すなわち、判定範囲を順次移動して上述の閾値判定を繰り返し実行することにより、3次元原画像全体に対して閾値判定を行う。
【００２０】
ここで、閾値判定の対象となる「注目画素の周囲の複数画素(判定範囲)」とは、図2に示すように、3次元原画像30に含まれる注目画素31を中心とするX方向、Y方向、スライス(SLS)方向の所定距離範囲内にある画素である。すなわち、注目画素31の近傍に位置する画素である。本発明では、閾値判定の対象となる注目画素31とその周囲の複数画素を判定範囲33と呼ぶ。例えば、判定範囲33は、注目画素31を中心とするスライス(SLS)方向、左右(X)方向、奥行き(Y)方向の3×3×3の画素群とする。なお、判定範囲33の広さはこれに限定されず、5×3×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を設けるようにしてもよい。
【００３１】
次に、図3〜図8を参照して、画像処理装置100の動作について説明する。
【００３２】
画像処理装置100のCPU101は、主メモリ102または記憶装置103から図3の診断用画像生成に関するプログラム及びデータを読み出し、このプログラム及びデータに基づいて処理を実行する。
【００３３】
なお、診断用画像生成処理の実行開始に際して、処理対象とする画像データは画像データベース111等からネットワーク110及び通信I/F104を介して取り込まれ、画像処理装置100の記憶装置103に記憶されているものとする。
【００３４】
処理対象とする画像データは、例えば、対象領域を含む複数の断層像を積み上げてなる3次元原画像とする。例えば、図2に示すように複数スライスSLSn-2、SLSn-1、SLSn、SLSn+1、・・・の断層像からなる3次元原画像30が処理対象の画像データとして読み込まれるものとする。また、画像データの好適な例として、CT画像またはMR画像等が挙げられる。本実施の形態では、血管の石灰化領域を残しつつ、骨や軟骨を適切に除去したMIP画像(診断用画像)を最終的に生成することを目的とする。そのため、処理対象とする3次元原画像30は頭部の画像としているが、本発明の適用範囲は頭部画像に限定されない。例えば、骨の断面画像について骨の内部組織である海綿骨を骨髄と分けて認識したり、胸部画像について濃度値にゆらぎのある肋骨を認識したりする場合等にも適用可能である。
【００３５】
CPU101は、読み込んだ画像データのうち、処理対象とする断層像群(3次元原画像)を読み込む(ステップS101)。
【００３６】
次に、CPU101は、閾値条件の決定の仕方、閾値判定処理の処理モードの設定を行う(ステップS102)。ステップS102において、CPU101は、例えば図4に示すような操作画面20を表示装置107に表示し、操作者による各種設定、入力操作を受け付ける。
【００３７】
図4に示すように、操作画面20には、画像表示エリア21、処理モード選択欄22、閾値決定欄23、演算ボタン24、終了ボタン25等の各種ボタンが設けられる。
【００３８】
画像表示エリア21には、処理対象とする3次元原画像30の全部または一部や、処理結果として生成された診断用画像等が表示される。
【００３９】
処理モード選択欄22は、閾値判定の処理方法(処理モード)を選択する欄であり、単体モード、等方的近隣判定モード、異方的近隣判定モードの各ラジオボタン221、222、223が設けられる。
【００４０】
単体モードとは、注目画素単体について閾値判定を行う処理モードである。従来の一般的な閾値処理と同様に、注目画素だけを対象として所定の閾値で閾値処理するものである。
【００４１】
等方的近隣判定モードとは、注目画素とその周囲の画素とを同一の閾値条件を用いて閾値判定するモードである。
【００４２】
異方的近隣判定モードとは、注目画素を含む平面にある画素と、別の平面にある画素とで異なる閾値条件を用いて閾値判定するモードである。
【００４３】
図5は各モードにおける判定範囲と、適用する閾値条件の違いを説明する図である。
【００４４】
図5(a)に示すように、単体モードでは、注目画素31のみを判定範囲とする。また閾値は各注目画素で同じ閾値t1を用いる。CPU101は注目画素31を閾値t1で判定し、条件を満たす場合に注目画素31を認識画素とする。
【００４５】
また、図5(b)に示すように、等方的近隣判定モードでは、注目画素31を中心とする例えば3×3×3の画素範囲を判定範囲33とし、判定範囲33に含まれる各画素を全て閾値t1で判定する。そして、判定範囲33内のいずれかの画素が上述の閾値条件(t1)を満たす場合に注目画素31を認識画素とする。周囲の画素を判定範囲33に含むことで、軟骨のように濃度値にゆらぎのある組織を抽出しやすくなる。
【００４６】
また、図5(c)に示すように、異方的近隣判定モードでは、注目画素31を中心とする例えば3×3×3の画素範囲を判定範囲33とし、注目画素31を含む平面(スライスSLSn)には閾値t1を適用し、別の平面(スライスSLSn-1、SLSn+1)には異なる閾値t2を適用して、閾値判定する。そして、いずれかの画素において上述の閾値条件を満たす場合に注目画素31を認識画素とする。周囲の画素を判定範囲に含み、かつ注目画素を含む平面と別の平面とで異なる閾値を設定することで、閾値条件を広くとることができ、濃度値にゆらぎのある組織を更に抽出しやすくなる。これにより、更に軟骨を正確に認識できるようになる。
【００４７】
図4の説明に戻る。
【００４８】
操作画面20の閾値決定欄23には、閾値を自動設定する場合に操作される閾値自動設定ラジオボタン231、閾値を手動設定する場合に操作される閾値手動設定ラジオボタン232、閾値を手動設定する際に閾値を変更するためのスライダ233が設けられる。スライダ233はマウスポインタ28により移動させられる。閾値の設定については後述する。
【００４９】
演算ボタン24は、閾値判定処理を開始する際に操作されるボタンである。演算ボタン24が押下されると、CPU101は、処理モード選択欄22で設定された処理モード及び閾値決定欄23で決定された閾値条件にて、閾値判定処理を開始する。
【００５０】
終了ボタン25は、操作画面20を閉じ、診断用画像生成処理を終了する際に操作されるボタンである。
【００５１】
操作画面20において、閾値の決定方法及び処理モードが選択される(図23のステップS102)と、CPU101は処理モードに応じた判定範囲33を設定する。そして、CPU101は、選択された閾値決定方法及び処理モードに応じた閾値判定処理を実行する(ステップS103)。閾値判定処理の詳細については、後述する。
【００５２】
CPU101は、ステップS103の閾値判定処理により認識した画素を除去または抽出し、診断用画像を生成する(ステップS104)。診断用画像は、例えば、血管観察用であれば、3次元原画像から骨を除去したMIP画像や3次元ボリュームレンダリング画像等が一般的である。CPU101は、生成した診断用画像を記憶部103に保存するとともに表示装置107に表示し(ステップS105)、処理を終了する。
【００５３】
次に、図6を参照して、ステップS103の閾値判定処理について説明する。
【００５４】
閾値判定処理は、上述したように、単体モード、等方的近隣判定モード、及び異方的近隣判定モードの3つのモードのうちいずれかを用いて行われる。
【００５５】
本実施の形態では、頭部の血管の画像を描画するために骨や軟骨を除去した診断用画像(MIP画像)を生成するための最も好適な設定例として、異方的近隣判定モードが選択され、閾値は自動設定されることが望ましい。
【００５６】
図6は、異方的近隣判定モードによる閾値判定処理の流れを説明するフローチャートである。また閾値判定に用いる第1の閾値t1、第2の閾値t2は、予め手動設定されるか、自動算出されているものとする。閾値自動設定の方法については後述する。
【００５７】
図6に示すように、まずCPU101は、初めに処理する注目画素31と判定範囲33を指定する(ステップS201)。初めに処理する注目画素31等は、領域拡張法等で行われるように操作者が指定するものでなく、CPU101が決定する。例えば、CPU101は、3次元原画像30の端点から判定範囲33を順に移しながら走査し、最終的に3次元原画像30全体を閾値判定する。以下の説明では、一例として、注目画素31を中心としてスライスSLS方向、Y方向、X方向に1点ずつ広げた画素群(3×3×3＝27画素)を判定範囲33とする。
【００５８】
異方的近隣判定モードにおいて、CPU101は、まず、注目画素31を含むスライスSLSnの画素9点のどれかが、第1の閾値t1を満たすか否かを判定する(ステップS202)。注目画素31を含むスライスSLSnの画素9点のうち、いずれかの画素が第1の閾値t1を満たす場合は(ステップS202；Yes)、スライスSLSnの中心点、すなわち注目画素31を認識画素とする(ステップS203)。注目画素31を含む平面(スライスSLSn)の画素9点のうち、いずれの画素も第1の閾値t1を満たさない場合は(ステップS202；No)、別の平面(スライスSLSn-1)の画素9点のうち、いずれかの画素が第2の閾値t2を満たすか否かを判定する(ステップS204)。
【００５９】
スライスSLSn-1の画素9点のうち、いずれかの画素が第2の閾値t2を満たす場合は(ステップS204；Yes)、スライスSLSnの中心点、すなわち注目画素31を認識画素とする(ステップS203)。スライスSLSn-1の画素9点のうち、いずれの画素も第2の閾値t2を満たさない場合は(ステップS204；No)、別の平面(スライスSLSn+1)の画素9点のうち、いずれかの画素が第2の閾値t2を満たすか否かを判定する(ステップS205)。スライスSLSn+1の画素9点のうち、いずれかの画素が第2の閾値t2を満たす場合は(ステップS205；Yes)、スライスSLSnの中心点、すなわち注目画素31を認識画素とする(ステップS203)。スライスSLSn+1の画素9点のうち、いずれかの画素も第2の閾値t2を満たさない場合は(ステップS205；No)、注目画素31を認識画素としないで、次の判定範囲33の処理に移る。
【００６０】
すなわち、CPU101は、3次元原画像30のすべての画素について閾値判定が終了したか否かを判定し(ステップS206)、終了していない場合は(ステップS206；No)、判定範囲30を1画素分または所定画素数分隣に移動する(ステップS207)。
【００６１】
CPU101は、次の判定範囲33を対象に、第1の閾値t1、第2の閾値t2による閾値判定を繰り返し行う(ステップS2012〜ステップS205)。すなわち、次の判定範囲33のうち、その中心画素(注目画素31)を含む平面に対して第1の閾値t1を適用し、別の平面に対して第2の閾値t2を適用し、いずれかの画素が閾値条件を満たす場合に、注目画素31を認識画素とする。
【００６２】
3次元原画像30のすべての画素について閾値判定が終了すると(ステップS206；Yes)、閾値判定処理を終了する。
【００６３】
なお、上述の処理において、認識画素を保持する方法としては例えば、以下の(a)、(b)に示すような方法があり、どの方法を用いてもよい。
(a)3次元原画像中に目印となる特定値を直接書き込む。
(b)3次元原画像とは別のメモリに二値画像として記録する。
【００６４】
また、診断用画像(MIP画像等)を演算中に、診断用画像の生成に用いる画素について上述の閾値判定処理を実行して各画素を認識または非認識するようにしてもよい。この場合は、診断用画像の生成において不要な画素については閾値判定処理を行わないものとしてもよい。
【００６５】
以上説明した閾値判定処理は、異方的近隣判定モードの場合であるが、等方的近隣判定モードの場合は、判定範囲33(例えば、注目画素を中心とする3×3×3の画素)内のすべての画素に対して同一の閾値t1を適用し、いずれかの画素が閾値t1を満たす場合に、注目画素31を認識画素とする。この閾値判定を、判定範囲33を移動させながら画像全体にわたって繰り返し実行する。
【００６６】
次に、閾値を自動設定する方法について、図7、図8を参照して説明する。
図7(a)は、頭部のCT画像の濃度値(CT値)のヒストグラム60を模式的に示す図であり、図7(b)は、胸部のCT画像の濃度値(CT値)のヒストグラム63を模式的に示す図である。
【００６７】
図7(a)に示すように、頭部の濃度値ヒストグラム60には、頭部軟部組織を示すピーク61と、頭蓋骨や軟骨に対応する小さなピーク62が現れるという特徴がある。CPU101は、この頭蓋骨や軟骨に対応する小さなピーク62の位置情報(CT値範囲)を基に骨領域の閾値を自動的に決定することができる。
【００６８】
一方、図7(b)に示すように、胸部の濃度値ヒストグラム63には、胸部や頚部の軟部組織を示すピーク64は現れるが、肋骨や軟骨のCT値(図7(b)の破線部65)を示すはっきりとしたピークが現れない。このように、区別したい部位(肋骨や軟骨)の特徴が現れない場合は、図8に示す方法で閾値t1(閾値t2)を決定する。
【００６９】
すなわち、CPU101は、まず、図8(a)に示すようにCT画像のヒストグラム7を求める。求めたヒストグラム7のピーク(頂点70)とピーク位置(頂点70)におけるCT値を求める。このヒストグラム7の分布曲線のうち、頂点70よりCT値が低い分布曲線を71a、頂点70よりCT値が高い分布曲線を71bと呼ぶ。
【００７０】
なお、図8(a)は、分布のピークが一つの例であるが、ピークが複数現れる場合もある。この場合は、例えば、骨の閾値を求める場合はピーク位置のCT値が最も高いピークを使うことができる。
【００７１】
次に、CPU101は、頂点70を境界位置として、CT値が低い側の分布曲線71aに対称な分布曲線72を、CT値が高い側に求め、対称分布曲線71a、72とする。そして、元のヒストグラム7の分布曲線71a、71bから、対称分布曲線71a、72を減算し、減算後の分布曲線73を求める。図8(b)に減算後分布曲線73の例を示す。
図8に示す分布曲線の横軸はCT値、縦軸は頻度の差分値である。
CPU101は、減算後分布曲線73に基づいて閾値を算出する。
【００７２】
例えば、図8(b)に示す減算後分布曲線73の頂点位置CT1と頂点位置CT1における差分値A1を求め、差分値がA1/NとなるCT値のうち値が大きい方のCT値(CT2)を、第1の閾値t1と設定する。ここで、Nは所定の定数である。また、閾値t1は、減算後分布曲線73の期待値expCTを使って、f(expCT)として求めてもよい。なお、f()は任意の関数を示し、期待値expCTは、例えば次式を用いて算出すればよい。
【００７３】
【数１】
2014038428.tif 000003
【００７４】
ここで、iはCT値の階級を表すインデックス、CTiはインデックスiに対応するCT値(図8(b)の横軸の各値)、Fiはインデックスiに対応する頻度(図8(b)の縦軸の各値)である。
【００７５】
また、第2の閾値t2は、以下の式(2)のように求めればよい。
【００７６】
t2＝t1−B    ・・・(2)
ここで、Bは、軟骨の濃度値の変動幅(30〜50程度)に含まれる値とすればよい。
【００７７】
図9は、図8の方法で自動決定した閾値と、異方的近隣判定モードで閾値判定処理により骨を除去したMIP画像を示す。閾値t2と閾値t1との差は50とした(t2＝t1−50)。
【００７８】
比較例として、単体モード、閾値手動設定により閾値判定処理を行った場合に得た骨除去後MIP画像を図10に示す。
【００７９】
図9に示す画像は、図10に示す画像と比較して、特に画像中央部の鼻の周囲の軟骨が除去され、血管走行が観察しやすくなっている。また、血管の石灰化領域(画像の濃度値が高い(白い)部分)は残されている。
【００８０】
なお、閾値を手動設定する場合は、操作者は、図4に示す操作画面20の閾値決定欄23で閾値手動設定ラジオボタン232を選択し、閾値設定用のスライダ233を操作して、閾値を適宜設定する。
【００８１】
以上説明したように、本発明に係る画像処理装置100によれば、CPU101は、複数の断層像を積み上げてなる3次元原画像30に含まれる注目画素31の周囲の複数の画素(判定範囲33)に対して所定の閾値条件を適用して閾値判定し、閾値条件を満たす場合に注目画素を認識画素とする。また、CPU101は、従来の領域拡張法のように操作者が指定した開始点から領域を拡張するのではなく、注目画素を順次別の画素に移動して上述の閾値判定を繰り返し実行することにより、3次元原画像30全体に対して閾値判定を行うものとする。
【００８２】
これにより、軟骨のように、濃度値にゆらぎ(変動幅)がある組織をより精度よく認識することができる。更に、画像全体を走査して閾値判定処理を行うため、開始点の指定操作が不要となる。従来の領域拡張法を用いる場合は、点在する領域を抽出(認識)する際には、それぞれの領域に操作者が開始点を指定する必要があったが、本発明を適用すれば、開始点の指定操作は不要となり操作が簡単である。
【００８３】
また、閾値判定処理では、判定範囲33のうち、注目画素31と同一平面にある画素と、別の平面にある画素とで異なる閾値を適用し、判定範囲33内のいずれかの画素で閾値条件を満たすと判定した場合に、注目画素31を認識画素とする(異方的近隣判定モード)。これにより、閾値条件を広くとることができ、特に、頭部の軟骨や骨の内部組織、肋骨等のように、濃度値にゆらぎのある組織を精度よく認識できるようになる。
【００８４】
また、閾値判定処理では、注目画素と周囲の複数の画素(判定範囲33内の各画素)に対し、同一の閾値を適用し、いずれかの画素で閾値条件を満たすと判定した場合に、注目画素31を認識画素とするようにしてもよい(等方的近隣判定モード)。これにより、注目画素31自体が閾値条件を満たしていなくても、近傍の画素が閾値条件を満たしていれば注目画素を認識画素と判定できるため、濃度値にゆらぎのある組織の認識に好適である。
【００８５】
また、閾値条件は、3次元原画像30の濃度値のヒストグラムに基づいて算出されることが望ましい。これにより、元となる画像の濃度値分布の特徴に適合した閾値を設定することができるため、より精度のよい結果を得ることができる。
【００８６】
また、閾値条件を算出する際は、3次元原画像30の濃度値のヒストグラム(71a、71b)のピーク位置(頂点70)を求め、ピーク位置(頂点70)より濃度値の低い部分のデータ(分布曲線71a)を使用してピーク位置を頂点として左右対称な形状の分布曲線72を求め、元のヒストグラムの分布曲線71a、71bから左右対称な形状の分布曲線71a、72を減算することにより減算後分布曲線73を求め、減算後分布曲線73の特徴に基づいて閾値条件を求めるようにしてもよい。
【００８７】
これにより、例えば胸部画像のように、認識したい組織を示す特徴的な分布形状を持たないヒストグラムからでも好適な閾値を自動設定できる。
【００８８】
また、閾値条件を操作者に設定させるために、例えば操作画面20に閾値を変更するスライダ233等のGUI(Graphical User Interface)を更に備えるようにすることが望ましい。
【００８９】
また、操作画面20には、閾値判定の処理モードとして、注目画素単体について閾値判定を行う「単体モード」か、注目画素と周囲の画素とを同一の閾値条件を用いて閾値判定する「等方的近隣判定モード」か、注目画素を含む平面と別の平面とで異なる閾値条件を用いて閾値判定する「異方的近隣判定モード」か、を選択する処理モード選択欄22を設けることが望ましい。これにより、容易に処理モードを切り替えて閾値判定処理を実行することができる。
【００９０】
また、操作画面20には、閾値条件を操作者に手動設定させるか、或いは自動設定するかを選択する閾値決定欄23を設けることが望ましい。これにより、容易に閾値の設定方法を切り替えて閾値判定処理を実行することができる。
【００９１】
また、上述の閾値判定処理により認識された認識画素を3次元原画像30から除去または抽出して、3次元画像またはMIP画像等の診断用画像を生成することが望ましい。これにより、診断用画像の生成のために本発明を利用できる。例えば、血管走行観察用には、頭蓋骨や軟骨を適切に除去した好適なMIP画像を生成できるようになる。
【００９２】
なお、上述の第1の実施の形態では、処理モードとして異方的近隣判定モードを設定し、閾値は自動設定とする例について説明したが、この組み合わせに限定されるものではない。例えば、異方的近隣判定モードで閾値は手動設定としてもよい。また、等方的近隣判定モードで閾値を自動設定としてもよいし、等方的近隣判定モードで閾値は手動設定としてもよい。
【００９３】
また、上述の説明では、異方的近隣判定モードでは、スライス方向に異なる閾値を適用する例を示したが、これに限定されない。例えば、図11に示すように、X方向に異なる閾値(第1の閾値t1、第2の閾値t2)を適用するようにしてもよい。すなわち、注目画素を含むY-SLS平面(X＝Xの面)には閾値t1を適用し、その隣のY-SLS平面(X＝X−1の面及びX＝X＋1の面)には閾値t2を適用するようにしてもよい。同様に、Y方向に異なる閾値を適用するようにしてもよい。
【００９４】
また、上述の説明では、閾値判定により領域を認識(抽出/除去)した後に診断用画像(MIP画像)を生成すると説明したが、閾値判定のタイミングと診断用画像生成のタイミングとの時間的前後関係は必ずしもこの順番に限定されない。診断用画像の演算を行うと同時に上述の閾値判定処理を行いつつ、認識画素を活用/無視して診断用画像を生成するようにしてもよい。
【００９５】
［第2の実施の形態］
次に、図12を参照して、第2の実施の形態について説明する。
第2の実施の形態に係る画像処理装置100では、異方的近隣判定モードで判定範囲を広くとる場合に、注目画素からの距離に応じて異なる閾値条件を設定する。
【００９６】
なお、第2の実施の形態の画像処理装置100のハードウエア構成は、第1の実施の形態と同様であるため、重複する説明を省略し、同一の各部は同一の符号を付して説明する。
【００９７】
すなわち、第2の実施の形態では、閾値判定処理において、CPU101は、図12に示す判定範囲33aのように、注目画素の周囲に比較的広い判定範囲を設定する。そして、注目画素を含む平面SLSnには第1の閾値t1を適用し、1画素隣の平面SLSn-1、SLSn+1には第2の閾値t2を適用し、2画素隣の平面SLSn-2、SLSn+2には第3の閾値t3を適用する。そして、CPU101は、判定範囲33a内のいずれかの画素で閾値条件を満たすと判定した場合に、注目画素を認識画素とする。
【００９８】
閾値を自動設定する場合は、第1の実施の形態と同様に、3次元原画像30の濃度値ヒストグラムに基づいて認識しようとする組織のCT値範囲を決定し、閾値t1を設定する。また、別の閾値t2は、上述の式(2)(t2＝t1-B)を用いて算出し、閾値t3は、以下の式(3)を用いて算出すればよい。
【００９９】
t3＝t1−C    ・・・(3)
ここで、Cは、Bとは異なる値とする。例えば、軟骨を認識するためには、軟骨の濃度値の変動幅(30〜50程度)の範囲内の数値とする。
【０１００】
以上説明したように、第2の実施の形態の画像処理装置100は、異方的近隣判定モードにおいて、注目画素からの距離に応じて異なる閾値条件を設定し、判定範囲33aのいずれかの画素で閾値条件を満たす場合には、注目画素を認識画素とする。
【０１０１】
これにより、濃度値の変動幅が大きく、また分布が複雑な組織を更に精度よく認識できるようになる。
【０１０２】
なお、図12の例では、3つの閾値t1、t2、t3を設定するものとしたが、判定範囲をより広くとる場合には、閾値条件を更に多く設定し、注目画素からの距離に応じて徐々に閾値条件を変更するようにすればよい。また閾値の適用範囲も1画素ずつに限定されず、ある程度の幅(画素数)毎に異なる閾値を適用するようにしてもよい。
【０１０３】
以上、添付図面を参照しながら、本発明に係る画像処理装置の好適な実施形態について説明したが、本発明はかかる例に限定されない。当業者であれば、本願で開示した技術的思想の範疇内において、各種の変更例又は修正例に想到し得ることは明らかであり、それらについても当然に本発明の技術的範囲に属するものと了解される。
【符号の説明】
【０１０４】
1  画像処理システム、100  画像処理装置、101  CPU、102  主メモリ、103  記憶装置、104  通信I/F、105  表示メモリ、106  I/F、107  表示装置、108  マウス、109  入力装置、110  ネットワーク、111  画像データベース、112  医用画像撮影装置、113  バス、20  操作画面、21  画像表示領域、22  処理モード選択欄、23  閾値決定欄、30  3次元原画像、31  注目画素、33、33a  判定範囲(注目画素の周囲の画素)、7  濃度値ヒストグラム、70  濃度値ヒストグラムの頂点、71a、71b  分布曲線、72  対称分布曲線、73  減算後分布曲線、t1、t2、t3  閾値
【図１】
2014038428.tif 000004
【図２】
2014038428.tif 000005
【図３】
2014038428.tif 000006
【図４】
2014038428.tif 000007
【図５】
2014038428.tif 000008
【図６】
2014038428.tif 000009
【図７】
2014038428.tif 000010
【図８】
2014038428.tif 000011
【図９】
2014038428.tif 000012
【図１０】
2014038428.tif 000013
【図１１】
2014038428.tif 000014
【図１２】
2014038428.tif 000015
【手続補正書】
【提出日】平成27年4月21日(2015.4.21)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
画像から特定の領域を認識する処理を実行する画像処理装置であって、
前記画像に含まれる画素のうち、注目画素と周囲の複数の画素を含む判定範囲に対して所定の閾値条件を適用して閾値判定し、前記判定範囲内のいずれかの画素で前記閾値条件を満たす場合に前記注目画素を認識画素とする閾値判定部と、
前記判定範囲を順次移動して前記閾値判定を繰り返し実行する処理実行部と、
を備えることを特徴とする画像処理装置。
【請求項２】
前記閾値判定部は、前記注目画素と同一平面にある画素と、別の平面にある画素とで異なる閾値条件を適用し、前記判定範囲内のいずれかの画素で前記閾値条件を満たすと判定した場合に、前記注目画素を前記認識画素とすることを特徴とする請求項1に記載の画像処理装置。
【請求項３】
前記閾値条件として、前記注目画素からの距離に応じて異なる閾値を設定することを特徴とする請求項2に記載の画像処理装置。
【請求項４】
前記閾値判定部は、前記判定範囲内のすべての画素に対し、同一の閾値条件を適用し、いずれかの画素で前記閾値条件を満たすと判定した場合に、前記注目画素を前記認識画素とすることを特徴とする請求項1に記載の画像処理装置。
【請求項５】
前記画像の濃度値分布のヒストグラムに基づいて前記閾値条件を算出する閾値自動算出部を更に備え、
前記閾値判定部は、前記閾値自動算出部により算出された閾値条件を適用して前記閾値判定を行うことを特徴とする請求項1に記載の画像処理装置。
【請求項６】
前記閾値自動算出部は、
前記画像の濃度値のヒストグラムの最大のピーク位置を求め、前記ピーク位置より濃度値の低い部分のデータを使用して前記ピーク位置を頂点として左右対称な形状の分布曲線を求め、前記ヒストグラムの分布曲線から前記左右対称な形状の分布曲線を減算することにより減算後分布曲線を求め、減算後分布曲線の特徴に基づいて前記閾値条件を求めることを特徴とする請求項5に記載の画像処理装置。
【請求項７】
前記閾値条件を操作者に設定させるための閾値手動設定部を更に備え、
前記閾値判定部は、前記閾値手動設定部により設定された閾値条件を適用して前記閾値判定を行うことを特徴とする請求項1に記載の画像処理装置。
【請求項８】
前記閾値判定における処理モードとして、前記注目画素の単体について閾値判定を行うか、前記判定範囲に含まれる全ての画素を同一の閾値条件を用いて閾値判定するか、前記注目画素と同一平面にある画素と、別の平面にある画素とで異なる閾値条件を用いて閾値判定するか、を選択する処理モード選択部を更に備えることを特徴とする請求項1に記載の画像処理装置。
【請求項９】
前記閾値条件を操作者に手動設定させるか、或いは前記閾値条件を前記画像の濃度値分布のヒストグラムに基づいて算出するかを選択する閾値設定方法選択部を更に備え、
前記閾値判定部は、前記閾値設定方法選択部により選択された方法で前記閾値条件を設定することを特徴とする請求項1に記載の画像処理装置。
【請求項１０】
前記閾値判定部により認識された認識画素を前記画像から除去または抽出した診断用画像を生成する診断用画像生成部を更に備えることを特徴とする請求項1に記載の画像処理装置。
【請求項１１】
コンピュータを用いて、画像から特定の領域を認識する画像処理方法であって、
前記画像に含まれる画素のうち、注目画素と周囲の複数の画素を含む判定範囲に対して所定の閾値条件を適用して閾値判定し、前記判定範囲内のいずれかの画素で前記閾値条件を満たす場合に前記注目画素を認識画素とする閾値判定処理を、前記判定範囲を順次移動して繰り返し実行することを特徴とする画像処理方法。
【手続補正２】
【補正対象書類名】明細書
【補正対象項目名】０００６
【補正方法】変更
【補正の内容】
【０００６】
しかしながら、上述したような領域拡張法や組織の分離手法を用いても、軟骨のように濃度値のゆらぎ(変動)が大きい組織を画像中から正確に区別して認識することは困難であった。また、頭部の画像において、軟骨領域はすべてが連なっているわけでなく点在するため、領域拡張法を用いる場合には、操作者がひとつずつ処理の開始点を指定する必要があった。そのため、操作が非常に煩雑であった。
【手続補正３】
【補正対象書類名】明細書
【補正対象項目名】００３７
【補正方法】変更
【補正の内容】
【００３７】
図4に示すように、操作画面20には、画像表示エリア21、判定処理モード選択欄22、閾値決定欄23、演算ボタン24、終了ボタン25等の各種ボタンが設けられる。
【手続補正４】
【補正対象書類名】明細書
【補正対象項目名】００３９
【補正方法】変更
【補正の内容】
【００３９】
判定処理モード選択欄22は、閾値判定の処理方法(処理モード)を選択する欄であり、単体モード、等方的近隣判定モード、異方的近隣判定モードの各ラジオボタン221、222、223が設けられる。
【手続補正５】
【補正対象書類名】明細書
【補正対象項目名】００４９
【補正方法】変更
【補正の内容】
【００４９】
演算ボタン24は、閾値判定処理を開始する際に操作されるボタンである。演算ボタン24が押下されると、CPU101は、判定処理モード選択欄22で設定された処理モード及び閾値決定欄23で決定された閾値条件にて、閾値判定処理を開始する。
【手続補正６】
【補正対象書類名】明細書
【補正対象項目名】００５１
【補正方法】変更
【補正の内容】
【００５１】
操作画面20において、閾値の決定方法及び処理モードが選択される(図3のステップS102)と、CPU101は処理モードに応じた判定範囲33を設定する。そして、CPU101は、選択された閾値決定方法及び処理モードに応じた閾値判定処理を実行する(ステップS103)。閾値判定処理の詳細については、後述する。
【手続補正７】
【補正対象書類名】明細書
【補正対象項目名】００６０
【補正方法】変更
【補正の内容】
【００６０】
すなわち、CPU101は、3次元原画像30のすべての画素について閾値判定が終了したか否かを判定し(ステップS206)、終了していない場合は(ステップS206；No)、判定範囲33を1画素分または所定画素数分隣に移動する(ステップS207)。
【手続補正８】
【補正対象書類名】明細書
【補正対象項目名】００６１
【補正方法】変更
【補正の内容】
【００６１】
CPU101は、次の判定範囲33を対象に、第1の閾値t1、第2の閾値t2による閾値判定を繰り返し行う(ステップS202〜ステップS205)。すなわち、次の判定範囲33のうち、その中心画素(注目画素31)を含む平面に対して第1の閾値t1を適用し、別の平面に対して第2の閾値t2を適用し、いずれかの画素が閾値条件を満たす場合に、注目画素31を認識画素とする。
【手続補正９】
【補正対象書類名】明細書
【補正対象項目名】００７１
【補正方法】変更
【補正の内容】
【００７１】
次に、CPU101は、頂点70を境界位置として、CT値が低い側の分布曲線71aに対称な分布曲線72を、CT値が高い側に求め、対称分布曲線71a、72とする。そして、元のヒストグラム7の分布曲線71a、71bから、対称分布曲線71a、72を減算し、減算後の分布曲線73を求める。図8(b)に減算後分布曲線73の例を示す。
図8(b)に示す分布曲線の横軸はCT値、縦軸は頻度の差分値である。
CPU101は、減算後分布曲線73に基づいて閾値を算出する。
【手続補正１０】
【補正対象書類名】明細書
【補正対象項目名】００８９
【補正方法】変更
【補正の内容】
【００８９】
また、操作画面20には、閾値判定の処理モードとして、注目画素単体について閾値判定を行う「単体モード」か、注目画素と周囲の画素とを同一の閾値条件を用いて閾値判定する「等方的近隣判定モード」か、注目画素を含む平面と別の平面とで異なる閾値条件を用いて閾値判定する「異方的近隣判定モード」か、を選択する判定処理モード選択欄22を設けることが望ましい。これにより、容易に処理モードを切り替えて閾値判定処理を実行することができる。
【手続補正１１】
【補正対象書類名】明細書
【補正対象項目名】０１０４
【補正方法】変更
【補正の内容】
【０１０４】
1  画像処理システム、100  画像処理装置、101  CPU、102  主メモリ、103  記憶装置、104  通信I/F、105  表示メモリ、106  I/F、107  表示装置、108  マウス、109  入力装置、110  ネットワーク、111  画像データベース、112  医用画像撮影装置、113  バス、20  操作画面、21  画像表示領域、22  判定処理モード選択欄、23  閾値決定欄、30  3次元原画像、31  注目画素、33、33a  判定範囲(注目画素の周囲の画素)、7  濃度値ヒストグラム、70  濃度値ヒストグラムの頂点、71a、71b  分布曲線、72  対称分布曲線、73  減算後分布曲線、t1、t2、t3  閾値
【手続補正１２】
【補正対象書類名】図面
【補正対象項目名】図３
【補正方法】変更
【補正の内容】
【図３】
2014038428.tif 000016
【手続補正１３】
【補正対象書類名】図面
【補正対象項目名】図４
【補正方法】変更
【補正の内容】
【図４】
2014038428.tif 000017
【手続補正１４】
【補正対象書類名】図面
【補正対象項目名】図８
【補正方法】変更
【補正の内容】
【図８】
2014038428.tif 000018
【国際調査報告】
2014038428.tif 000019
2014038428.tif 000020
2014038428.tif 000021
2014038428.tif 000022