JP2006239195A - コンピュータによる画像診断支援方法および画像診断支援装置ならびにプログラム - Google Patents

コンピュータによる画像診断支援方法および画像診断支援装置ならびにプログラム Download PDF

Info

Publication number
JP2006239195A
JP2006239195A JP2005060028A JP2005060028A JP2006239195A JP 2006239195 A JP2006239195 A JP 2006239195A JP 2005060028 A JP2005060028 A JP 2005060028A JP 2005060028 A JP2005060028 A JP 2005060028A JP 2006239195 A JP2006239195 A JP 2006239195A
Authority
JP
Japan
Prior art keywords
image
displacement
respiratory
images
region
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Withdrawn
Application number
JP2005060028A
Other languages
English (en)
Inventor
Satoru Osawa
哲 大沢
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fujifilm Holdings Corp
Original Assignee
Fuji Photo Film Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Fuji Photo Film Co Ltd filed Critical Fuji Photo Film Co Ltd
Priority to JP2005060028A priority Critical patent/JP2006239195A/ja
Priority to US11/367,300 priority patent/US20060239530A1/en
Publication of JP2006239195A publication Critical patent/JP2006239195A/ja
Withdrawn legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/0002—Inspection of images, e.g. flaw detection
    • G06T7/0012—Biomedical image inspection
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing
    • G06T2207/30061—Lung
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V2201/00—Indexing scheme relating to image or video recognition or understanding
    • G06V2201/03—Recognition of patterns in medical or anatomical images

Landscapes

  • Engineering & Computer Science (AREA)
  • Quality & Reliability (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Health & Medical Sciences (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Image Processing (AREA)
  • Image Analysis (AREA)

Abstract

【課題】異なる呼吸相を表す複数の胸部医用画像に基づいて被検者の呼吸機能の異常の有無を判断する際に、異常箇所の見落としを減らすとともに、読影時間の短縮による診断効率の向上を図る。
【解決手段】被検者の呼吸活動に伴う呼吸相の変化の順に並べられた胸部医用画像Ｐ1からＰNについて、変位算出手段２が、胸部医用画像Ｐ1中の局所領域の各々の位置を、画像Ｐ2から順に画像ＰNまでの対応する位置に合わせる位置合わせ処理を順次行って、各画像間の呼吸相の変化における各局所領域の変位を算出し、散布度算出手段３が、各画像間における局所領域の変位毎に変位の差分の散布度Ｖ1からＶN-2を算出し、判定手段４が、散布度Ｖ1からＶN-2のいずれかが閾値Thよりも大きい場合には、被検者の呼吸機能に異常があると判定する。
【選択図】 図４

Description

本発明は、コンピュータによる画像診断支援方法および画像診断支援装置ならびにプログラムに関し、より詳細には、異なる呼吸相を表す複数の胸部医用画像に基づいて被検者の呼吸機能の異常の有無を判定する方法および装置、ならびに、コンピュータにこの判定を行わせるプログラムに関するものである。
近年、吸気から呼気までの複数の呼吸段階（呼吸相）を撮影した胸部単純Ｘ線画像を用いて呼吸機能の診断を行おうとする試みが行われている。
例えば、Ｉ.Ｉ.(Image Intensifier）を用いて吸気から呼気までの複数の呼吸相を撮影した動画像から換気異常を検出することが可能であり、このような換気異常の検出をより高い精度で行うためには複数の画像間の位置合わせが必要であるとの見解が報告されている（例えば、非特許文献１）。
また、動画対応のＦＰＤ(Flat-Panel Detector；平面検出器）を搭載したＸ線画像デジタル撮影装置を用いた検討も行われており、Ｉ.Ｉ.よりも高解像度の胸部動画像を用いることで、より精密に肺野の濃度の変化から呼吸機能情報を取得できることが示唆されており、さらに、撮影された動画像のうち隣接する画像間の位置合わせを行い、差分を計算した画像による差分動画像を作成した解析も試みられている（例えば、非特許文献２）。
このような背景の中で、本出願人は、異なる呼吸相を表す３以上の胸部医用画像をその呼吸相の変化の順に並べた場合に隣接しない２つの注目画像の各々について、注目画像に含まれる被写体の位置を、隣接する画像に含まれる被写体の位置に順次合わせていき、最終的に位置合わせの基準となる基準画像に含まれる被写体の位置に合わせることによって、異なる呼吸相の注目画像間の位置合わせの精度を向上するとともに、複数の呼吸相の画像間での濃度変化や、画像中の局所領域の位置の変化を視覚化した画像を生成することによって、複数の呼吸相を撮影した胸部医用画像の観察を容易にすることを提案している（例えば、特許文献１）。
特開２００５-２０７１９号公報 J.Liang, et.al, "Dynamic chest image analysis: model-based ventilation study with pyramid images",（米国）, Proceedings of SPIE Medical Imaging 1997: Physiology and Function from Multidimensional Images, SPIE, May 1997, Vol.3033, p.81-92 田中，真田，外５名，「胸部動画像を対象とした呼吸性動態の定量化」，医用画像情報学会雑誌，医用画像情報学会，2003年1月，第20巻，第1号, p.13-19
しかしながら、上記文献では、読影者が、複数の呼吸相を表す画像による動画や、それらの画像の差分画像による動画、濃度変化や位置の変化を表す画像を見て呼吸機能の異常の有無を判断しているため、異常箇所を見落とす可能性がある。また、このような見落としを減らすためには、時間をかけて画像を読影する必要があり、そうすると、診断効率が低下してしまうことになる。
本発明は、上記事情に鑑みてなされたものであり、読影者が、異なる呼吸相を表す複数の胸部医用画像に基づいて被検者の呼吸機能の異常の有無を判断する際に、異常箇所の見落としを減らすとともに、読影時間の短縮による診断効率の向上を図る方法および装置ならびにプログラムを提供することを目的とするものである。
本発明によるコンピュータを用いた画像診断支援方法は、被検者の呼吸運動における任意の３つの呼吸相を、その呼吸運動に伴う呼吸相の変化の順に第１、第２、第３の呼吸相とし、第１、第２、第３の呼吸相の各々を表す胸部医用画像を第１、第２、第３の画像とした場合において、第１の画像を構成する複数の局所領域の各々の位置を、第２の画像中の対応する位置に合わせる位置合わせ処理を行って、第１の呼吸相から第２の呼吸相への変化における局所領域の変位を表す第１の変位を局所領域毎に算出し、位置合わせ処理後の局所領域の各々の位置を、第３の画像中の対応する位置に合わせる位置合わせ処理を行って、第２の呼吸相から第３の呼吸相への変化における局所領域の変位を表す第２の変位を局所領域毎に算出し、局所領域毎の第１の変位と第２の変位の差分の散布度を算出し、算出された散布度が所定の閾値よりも大きい場合に被検者の呼吸機能に異常があると判定することを特徴とするものである。
また、本発明による画像診断支援装置は、この方法を実現するものである。すなわち、第１の画像を構成する複数の局所領域の各々の位置を、第２の画像中の対応する位置に合わせる位置合わせ処理を行って、第１の呼吸相から第２の呼吸相への変化における局所領域の変位を表す第１の変位を局所領域毎に算出するとともに、位置合わせ処理後の局所領域の各々の位置を、第３の画像中の対応する位置に合わせる位置合わせ処理を行って、第２の呼吸相から第３の呼吸相への変化における局所領域の変位を表す第２の変位を局所領域毎に算出する変位算出手段と、局所領域毎の第１の変位と第２の変位の差分の散布度を算出する散布度算出手段と、算出された散布度が所定の閾値よりも大きい場合に被検者の呼吸機能に異常があると判定する判定手段とを設けたことを特徴とするものである。
さらに、本発明による画像診断支援プログラムは、コンピュータに上記の方法を行わせる（上記の各手段として機能させる）ものである。すなわち、コンピュータを、第１の画像を構成する複数の局所領域の各々の位置を、第２の画像中の対応する位置に合わせる位置合わせ処理を行って、第１の呼吸相から第２の呼吸相への変化における局所領域の変位を表す第１の変位を局所領域毎に算出するとともに、位置合わせ処理後の局所領域の各々の位置を、第３の画像中の対応する位置に合わせる位置合わせ処理を行って、第２の呼吸相から第３の呼吸相への変化における局所領域の変位を表す第２の変位を局所領域毎に算出する変位算出手段と、局所領域毎の第１の変位と第２の変位の差分の散布度を算出する散布度算出手段と、算出された散布度が所定の閾値よりも大きい場合に被検者の呼吸機能に異常があると判定する判定手段として機能させることを特徴とするものである。
次に、本発明の詳細について説明する。
「第１、第２、第３の画像」は、異なる呼吸相を表す複数の胸部医用画像の各々について、画像中の肺野領域を検出し、検出された肺野領域の面積と、各画像から検出された肺野領域の面積の最大値または最小値との比率を各画像について算出し、算出された比率に基づいて選択されたものであることが好ましい。さらに、第１、第３の画像として、肺野領域の面積が最大となる最大吸気画像と、肺野領域の面積が最小となる最大呼気画像とを選択することが好ましい。
「局所領域」の大きさは、呼吸機能の観察に適した大きさとすることが好ましい。具体的には、２次小葉と呼ばれる肺胞が集まった６〜１０mm程度の組織と同程度の大きさにすることにより、局所的な換気機能の観察が効率的に行えるようになる。なお、各局所領域間で重複する部分が生じていてもよい。また、局所領域を１画素の領域としてもよい。
「位置合わせ処理」は、画像中の局所領域を個別にシフトさせる非線形変換処理であるが、前述の特許文献１記載の経時サブトラクション手法のように、線形の位置変換である概略的な位置合わせ処理（グローバルマッチング）と非線形の位置変換である局所的な位置合わせ処理（ローカルマッチング）とを組み合わせたものであってもよい。
「位置合わせ処理後の局所領域の各々の位置を、第３の画像中の対応する位置に合わせる位置合わせ処理」では、位置合わせ処理後の第１の画像中の局所領域の各々の位置を第３の画像中の対応する位置に合わせるようにしてもよいし、位置合わせ処理後の局所領域を第２の画像中に設定し、第２の画像中の局所領域の各々の位置を第３の画像中の対応する位置に合わせるようにしてもよい。
「第１の変位と第２の変位の差分」とは、第１の変位と第２の変位の間での変位の変化量を意味する。
「変位」と「差分」は、向きおよび／または大きさによって表現される量である。
「差分の散布度」は、差分の散らばりの度合を表すものであり、具体例としては、差分の分散や標準偏差等が考えられる。なお、ここでは、散布度が、その値が大きいほど散らばりが大きいという性質を有するものであることを前提としている。散布度の値が大きいほど差分の散らばりが小さいという関係となる場合には、後続の「判定」処理における大小関係を逆転させればよい。
「呼吸機能に異常があると判定する」のは「散布度が所定の閾値よりも大きい場合」だけでなく、散布度が所定の閾値と等しい場合も含めてもよい。すなわち、散布度が所定の閾値以上の場合に呼吸機能に異常があると判定してもよい。
「所定の閾値」は、その被検者の呼吸機能が正常な時の胸部医用画像に基づいて得られた散布度（以下、正常時散布度という）に基づいて設定された値とすることが好ましい。具体的には、正常時散布度に対してある程度の許容範囲をもたせた値を閾値とすることが考えられる。なお、この正常時散布度は予め求めておいてもよいし、判定対象の画像と正常時の画像の両方を準備しておき、判定の都度、正常時散布度を求めるようにしてもよい。
また、第１、第２、第３の画像中の対応する位置に注目領域を設定し、設定された注目領域内の局所領域について、散布度を算出するようにしてもよい。ここで、この注目領域を肺野領域内に設定することが好ましい。肺野領域は、公知の手法（例えば、特開２００３−００６６６１号公報）により検出することができる。なお、「注目領域」は、複数の局所領域が含まれるように設定する必要がある。また、注目領域は複数設定してもよく、その場合には、注目領域毎に散布度を算出してもよい。さらに、複数の注目領域の各々の間で重複する部分が生じていてもよい。
また、被検者の呼吸機能に異常があると判定された場合において、前記の差分が所定の閾値よりも大きい局所領域における第１および／または第２の変位を他の局所領域における第１および／または第２の変位と識別可能な態様にして表示することが好ましい。さらに、前述の注目領域内の局所領域について散布度を算出する態様では、前記の散布度が所定の閾値よりも大きい注目領域を他の領域と識別可能な態様にして表示するようにしてもよい。ここで、「識別可能な態様にして表示」の具体例としては、他とは異なる色にして表示することが考えられる。
本発明は、被検者の呼吸機能の異常の有無により、呼吸相の変化に伴う胸部医用画像中の局所領域毎の変位の方向や大きさの変化量が異なるという点に着目したものである。図１（ａ）は、呼吸機能が正常な場合において、呼気状態から吸気状態への呼吸相の変化に伴う胸部医用画像中の局所領域毎の変位をベクトルで模式的に表したものであり、図１（ｂ）は、このベクトルの向きと大きさの分布を表したものである。また、図２（ａ）（ｂ）は、呼吸機能に異常がある場合において、図１と同様に表したものである。図１（ａ）、図２（ａ）の上段では、呼吸相ｎからｎ＋２に変化するにつれて肺野内に空気が取り込まれ、肺野領域が拡張していく様子を表しており、図の下段では、特に図の上段の矩形領域に注目し、各呼吸相間における矩形領域内の局所領域の位置の変化をベクトル（以下、変位ベクトル）で表したものである。図に示したように、呼吸機能が正常な場合には、変位ベクトルの大きさと向きは呼吸相間でほぼ一定の分布となるのに対し（図１（ｂ）参照）、呼吸機能に異常がある場合には、呼吸相間で変位ベクトルの大きさや向きにばらつきが生じる（図２（ｂ）参照）。
本発明の画像診断支援方法および画像診断支援装置ならびにプログラムによれば、被検者の呼吸運動における任意の３つの呼吸相を表す胸部医用画像中の複数の局所領域の各々の位置を合わせる位置合わせ処理を行って、各局所領域の画像間の変位を求め、その変位の差分（変化量）の散布度が所定の閾値よりも大きい場合に被検者の呼吸機能に異常があると自動的に判定するので、判定結果に基づいて読影者の注意を喚起することが可能になり、異常箇所の見落としの減少や、読影時間の短縮による診断効率の向上に資する。
また、異なる呼吸相を表す複数の画像中から検出された肺野領域の面積と、各画像中の肺野領域の面積の最大値または最小値との比率に基づいて上記の処理対象の３つの画像を選択するようにした場合には、上記の比率に基づいて選択候補となる各画像の呼吸相が推定できるため、最大呼気画像や最大吸気画像のように、より診断に適した画像を上記の判定処理に用いることが可能になり、判定精度が向上し、読影者の診断精度や診断効率の向上が図られる。
上記の判定処理における所定の閾値を、同じ被検者の正常時の画像を入力として得られた散布度に基づいて設定された値にした場合には、異なる被検者間での個人差によって生じる判定誤差を減少させることが可能になり、さらなる判定精度の向上、診断精度や診断効率の向上に資する。
また、第１、第２、第３の画像中の対応する位置に設定された注目領域内の局所領域について、散布度を算出するようにした場合には、診断上より重要な領域における局所領域毎の呼吸相間での変位の向きや大きさの変化に基づいて、呼吸機能の異常の有無の判定ができるので、さらなる判定精度の向上、診断精度や診断効率の向上に資する。
被検者の呼吸機能に異常があると判定された場合において、第１の変位と第２の変位の差分が所定の閾値よりも大きい局所領域における第１および／または第２の変位を他の局所領域における第１および／または第２の変位と識別可能な態様にして表示したり、散布度が所定の閾値よりも大きい注目領域を他の領域と識別可能な態様にして表示したりすれば、読影者に対して、異常の有無だけでなく、異常箇所を指摘することが可能になるので、異常箇所の視認性が向上するため、さらなる診断効率の向上に資する。
以下、図面を参照しつつ、本発明の実施の形態について説明する。
図３は、本発明の実施形態となる胸部画像診断支援システムの構成の概要を示したものである。図に示すように、画像撮影・読取システム２０と、画像管理システム３０、画像処理システム１０、画像表示システム４０が、ＬＡＮ等のネットワーク５０を経由して通信可能な状態で接続されている。
画像撮影・読取システム２０は、患者の胸部を被写体として、複数の呼吸段階を撮影した放射線画像を取得するものであり、ＣＲ(Computed Radiography)装置２１や、Ｉ.Ｉ.２２、動画対応のＦＰＤを搭載したＸ線画像デジタル撮影装置（以下、ＦＰＤ装置という）２３等が含まれる。
画像処理システム１０は、画像撮影・読取システム２０で撮影された放射線画像に対して画像処理を行い、診断者の読影に適した画像を生成するものであり、画像処理サーバ１１等が含まれる。本発明の画像診断支援方法（装置）は、この画像処理サーバ１１に実装される。
画像管理システム３０は、画像撮影・読取システム２０や画像処理システム１０で生成された画像を保存・管理するものであり、画像管理サーバ３１と大容量外部記憶装置３２、データベース管理ソフトウェア（例えば、ＯＲＤＢ(Object Relational DataBase)管理ソフトウェア）等が含まれる。各画像は、被検者ＩＤ、被検者名、性別、生年月日等の被検者情報や、検査日時、検査部位、検査結果等の検査情報等と関連づけられて、データベース管理ソフトウェアによる管理の下で、大容量外部記憶装置３２に記憶されている。
画像表示システム４０は、画像撮影・読取システム２０や画像処理システム１０で生成された画像を表示するものであり、クライアントＰＣ４１や高精細液晶ディスプレイ４２等が含まれる。また、クライアントＰＣ４１のキーボードやマウス等の入力装置を用いて、読影対象の画像やその画像に対する画像処理方法等の読影条件の設定を行うことも可能になっている。
本発明の第１の実施形態では、被検者の呼吸運動における複数の呼吸相をその変化の順に表す胸部医用画像Ｐ1からＰNに基づいて、呼吸機能の異常を自動的に判定し、異常と判定された場合には画像中の局所領域の変位を表す局所変位ベクトル画像を生成する機能を実現する。なお、以下の説明では、画像Ｐ1が最大呼気状態（呼吸相１）を表す画像、画像ＰNが最大吸気状態（呼吸相Ｎ）を表す画像、また、画像Ｐnは呼吸相ｎを表す画像とする。
図４は、この機能が実装される画像処理サーバ１１の論理的な構成とデータの流れとを表したブロック図である。図に示すように、画像管理システム３０から画像Ｐ1、・・・、Ｐn、・・・、ＰNや後述の正常時散布度Ｖ0を取得する取得手段１と、取得された画像Ｐ1からＰNのうち、画像Ｐ1を複数の局所領域Ｒ1[1,1]、・・・、Ｒ1[i,j]、・・・、Ｒ1[I,J]に分割し、分割された局所領域Ｒ1[1,1]からＲ1[I,J]の各々の位置を、画像Ｐ2から順に画像ＰNまでの対応する位置に合わせる位置合わせ処理を順次行って、各画像間の呼吸相の変化における各局所領域Ｒ1[1,1]からＲ1[I,J]の変位Δ1[1,1]からΔ1[I,J]、・・・、Δn[1,1]からΔn[I,J]、・・・、ΔN-1[1,1]からΔN-1[I,J]を算出する変位算出手段２と、各画像間における局所領域の変位毎に変位の差分の散布度Ｖ1、・・・、Ｖn、・・・、ＶN-2を算出する散布度算出手段３と、算出された散布度Ｖ1からＶN-2の各々が、取得手段１によって取得された正常時散布度Ｖ0に基づいて設定された閾値Thよりも大きい場合に被検者の呼吸機能に異常があると判定する判定手段４と、各画像間における各局所領域の変位を可視化した局所変位ベクトル画像Ｑ1、・・・、Ｑn、・・・、ＱN-1を生成する局所変位ベクトル画像生成手段５と、呼吸機能が正常であると判定された場合には最大吸気画像ＰNをクライアントＰＣ４１に送信し、異常であると判定された場合には最大吸気画像ＰNと局所変位ベクトル画像Ｑ1からＱN-1を送信する送信手段６から構成される。
上記の各手段は、ＣＤ−ＲＯＭ等の記憶媒体から画像処理サーバ１１にインストールされたプログラムが実行されることによって、画像処理サーバ１１のＣＰＵ、主記憶装置、外部記憶装置、入出力インターフェース、オペレーティングシステム等と連携して実現される。また、このプログラムによって、上記の各手段による処理順序が制御される。
次に、本発明の第１の実施形態において行われる処理の流れについて説明する。
まず、画像撮影・読取システム２０において、例えば、動画対応のＦＰＤ装置２３によって、患者の胸部の撮影が行われ、複数の呼吸段階を表す画像Ｐ1（最大呼気）から画像ＰN（最大吸気）のＮ個（Ｎ≧３）の画像が生成される。生成された画像は、各々、被検者情報や検査情報等の付随情報と関連づけられ、画像毎にファイルに格納され、Ｎ個のファイルが出力される。出力されたＮ個のファイルは、画像管理システム３０にネットワーク５０経由で送信される。なお、本実施形態では、撮影時に最大吸気画像と最大呼気画像とが特定され、付随情報の中に最大吸気画像や最大呼気画像を特定する情報も含まれているものとする。
画像管理システム３０では、画像管理サーバ３１が、送信されてきたＮ個のファイルを受信し、データベース管理ソフトウェアで定められたデータ形式、データ構造に基づき、その受信したＮ個の画像ファイル中の画像データを付随情報と関連づけて大容量外部記憶装置３２に記憶する。記憶された画像データは、被検者ＩＤや検査日時等の付随情報の一部または全部を検索キーとして、検索可能な状態となっている。
読影者は、クライアントＰＣ４１の操作により、患者ＩＤや撮影日時等の付随情報と、読影パターン（呼吸機能異常の自動判定）とを指定し、処理の実行を要求する。ここで、上記の付随情報を指定することによって読影対象の画像Ｐ1からＰNが特定され、読影パターンを指定することによってその読影に必要な画像処理の内容が特定される。
クライアントＰＣ４１で指定された付随情報と読影パターンは、画像処理サーバ１１に送信され、画像処理サーバ１１では、受信した読影パターンに基づいて、その読影に必要な画像を生成する画像処理プログラム、すなわち、画像処理サーバ１１を図４の各手段として機能させるプログラムが起動される。
図５は、起動された画像処理プログラムによって実行される処理の流れを表すフローチャートである。また、図６は、変位算出手段２と散布度算出手段３で行われる処理を模式的に表したものである。
まず、取得手段１が、クライアントＰＣ４１から受信した付随情報に基づく検索条件による画像管理システム３０のデータベースに対する検索要求を送信する。画像管理サーバ３１では、受信した検索要求に応じてデータベース検索を行い、検索条件に合致する画像Ｐ1からＰNと、後述の正常時散布度Ｖ0を取得し、画像処理サーバ１１に送信する。画像処理サーバ１１の取得手段１は、画像Ｐ1からＰNと正常時散布度Ｖ0を受信し、画像処理サーバ１１の主記憶装置または外部記憶装置に一時的に記憶する（ステップa1）。
次に、画像を識別する添字ｎの値が１に設定された後（ステップa2）、変位算出手段２が、画像Ｐn（＝Ｐ1）を読み込み（ステップa3）、横Ｉ×縦Ｊ個の局所領域に分割し、局所領域Ｒ1[1,1]、・・・、Ｒ1[i,j]、・・・、Ｒ1[I,J]を設定する（ステップa4）。ここで、個々の局所領域は２次小葉と呼ばれる肺胞が集まった６〜１０mm程度の組織が含まれる程度の大きさの矩形領域であり、局所領域の数はこの大きさを基準として決定される。
そして、画像を識別する添字ｎの値が１加算され、ｎ＝２となった後（ステップa5）、変位算出手段２は、Ｐn（＝Ｐ2）を読み込み（ステップa6）、局所領域Ｒ1[1,1]からＲ1[I,J]の各々の位置を、画像Ｐ2中の対応する位置に合わせる位置合わせ処理を行って、呼吸相１から呼吸相２への変化における各局所領域Ｒ1[1,1]からＲ1[I,J]の変位を表す変位Δ1[1,1]からΔ1[I,J]を算出する（ステップa7）。ここでは、前述の特許文献１記載の経時サブトラクション手法を用いる。以下、詳細に説明する。
まず、画像Ｐ１と画像Ｐ２との概略的な位置合わせ処理（グローバルマッチング）を行う。これは、画像Ｐ2を画像Ｐ1に一致させるように、Ｐ２に対してアフィン変換（回転、平行シフト）を施す処理であり、この処理により画像Ｐ２は、図７に示すようにＰ２′に変換される。
グローバルマッチング処理の完了後、局所的な位置合わせ処理（ローカルマッチング）を行う。
具体的には、図８に示すように、画像Ｐ１中の各局所領域Ｒ1[1,1]からＲ1[I,J]の中心画素をそれぞれｘ−ｙ座標系（ｘ，ｙ）により表す。また、画像Ｐ２′に探索ＲＯＩＲ2'[1,1]、・・・、Ｒ2'[i,j]、・・・、Ｒ2'[I,J]を設定する。この探索ＲＯＩは、画像Ｐ１の各局所領域Ｒ1[1,1]からＲ1[I,J]に対応して設定され、同一の中心座標（ｘ，ｙ）を有し、各局所領域Ｒ1[1,1]からＲ1[I,J]よりも大きい領域である。ここでは、各局所領域Ｒ1[1,1]からＲ1[I,J]の４倍（縦横ともに２倍）の大きさの領域としている。
画像Ｐ２′に設定された各探索ＲＯＩの中で、画像Ｐ１の各局所領域Ｒ1[1,1]からＲ1[I,J]を移動させて、局所領域のマッチング度合いが最も高くなる位置（各局所領域Ｒ1[1,1]からＲ1[I,J]の中心位置（ｘ′，ｙ′））を探索ＲＯＩ毎に求める。マッチング度合いの高低を示す指標値としては、最小自乗法や相互相関による指標値を用いることができる。
そして、各局所領域Ｒ1[1,1]からＲ1[I,J]の中心画素（ｘ，ｙ）毎のシフト値（Δｘ，Δｙ）（ただし、Δｘ＝ｘ′−ｘ，Δｙ＝ｙ′−ｙ）を求め、呼吸相１から呼吸相２への変化における各局所領域（Ｒ1[1,1]からＲ1[I,J]）の変位（Δ1[1,1]からΔ1[I,J]）とする。図９は、画像Ｐ１と画像Ｐ２′の各局所領域の中心画素（ｘ，ｙ）と（ｘ′，ｙ′）、画像間での変位Δ1[1,1]を模式的に表したものである。図に示したとおり、局所領域の変位は、向きと大きさの両方を有するベクトル量である。
さらに、添字ｎの値が１加算されてｎ＝３となり（ステップa5）、変位算出手段２は、上記と同様にして、位置合わせ処理後の各局所領域Ｒ1[1,1]からＲ1[I,J]を画像Ｐ3′に設定された各探索ＲＯＩの中で移動させて、マッチング度合いが最も高くなる位置を探索ＲＯＩ毎に求め、さらに局所領域の中心画素毎のシフト値を求め、呼吸相２から呼吸相３への変化における各局所領域（Ｒ1[1,1]からＲ1[I,J]）の変位（Δ2[1,1]からΔ2[I,J]）とする（ステップa6, a7）。
そして、添字ｎの値がＮになるまでは（ステップa8）、上記と同様にして、添字ｎの値が１加算され（ステップa5）、変位算出手段２が、直前の位置合わせ処理後の局所領域Ｒ1[1,1]からＲ1[I,J]を画像Ｐn′に設定された各探索ＲＯＩの中で移動させて、マッチング度合いが最も高くなる位置を探索ＲＯＩ毎に求め、さらに局所領域の中心画素毎のシフト値を求め、呼吸相ｎ−１から呼吸相ｎへの変化における各局所領域（Ｒ1[1,1]からＲ1[I,J]）の変位（Δn-1[1,1]からΔn-1[I,J]）とする処理（ステップa6, a7）が繰り返し行われる。
添字の値がＮになったら（ステップa8）、すなわち、呼吸相１から順に呼吸相Ｎまでの変化における各局所領域（Ｒ1[1,1]からＲ1[I,J]）の変位（Δ1[1,1]からΔ1[I,J]、・・・、Δn[1,1]からΔn[I,J]、・・・、ΔN-1[1,1]からΔN-1[I,J]）が求まったら、散布度算出手段３が、１からＮ−２までの自然数について、各画像間における局所領域の変位毎に変位の差分（Δn＋1[1,1]−Δn[1,1]）、・・・、（Δn＋1[i,j]−Δn[i,j]）、・・・、（Δn＋1[I,J]−Δn[I,J]）を求め、これらの差分の分散を散布度Ｖnとして算出する（ステップa9）。これにより、散布度Ｖ1、・・・、Ｖn、・・・、ＶN-2の値が求まる。なお、局所領域の変位が向きと大きさを有するベクトル量であることから、その変位の差分もベクトル量となる。したがって、ここでは、この変位の差分を表すベクトルがｘ軸方向となす角度を表す値（rad、tan等）の分散を向きの散布度としたり、この変位の差分を表すベクトルの長さの分散を大きさの散布度としたりすることができる。
次に、判定手段４は、取得手段１が取得した正常時散布度Ｖ0を読み込み、αＶ0＋β（α、βは、臨床データに基づいて得られた許容係数）の値を閾値Thとし、散布度Ｖ1からＶN-2の中にその値が閾値Thよりも大きいものがあれば、被検者の呼吸機能に異常があると判定し、すべての散布度Ｖ1からＶN-2が閾値Th以下であれば、被検者の呼吸機能は正常であると判定する（ステップa10）。ここで、正常時散布度Ｖ0とは、診断対象の被検者が過去において呼吸機能が正常であると診断された時の胸部医用画像に基づいて、上記と同様にして算出された散布度の平均値である。この正常であると診断された時の画像は、今回と同様の方法（撮影と呼吸相のタイミング、撮影画像数等）で撮影されたものであることが好ましい。
判定手段４により正常と判定された場合には、送信手段６が、最大吸気画像ＰNの画像データのみをクライアントＰＣ４１に送信する（ステップa11）。
一方、判定手段４により異常と判定された場合には、局所変位ベクトル画像生成手段５が、１からＮ−１までの自然数ｎについて、画像Ｐn中の各局所領域の中心画素の位置を始点として、変位Δn[1,1]からΔn[I,J]を変位ベクトルとして向きや大きさを視覚化した局所変位ベクトル画像Ｑnを生成する（ステップa12）。図１０は、生成された局所変位ベクトル画像の一例である。図に示すように、各局所領域の変位は、各局所領域の中心画素の位置を始点とする棒によってその向きと大きさが表現されている。さらに、送信手段６が、生成された局所変位ベクトル画像Ｑ1からＱN-1と最大吸気画像ＰNの画像データをクライアントＰＣ４１に送信する（ステップa13）。
クライアントＰＣ４１は、判定手段４により正常と判定された場合には、画像処理サーバ１１から受信した画像データに基づいて最大吸気画像ＰNを高精細液晶ディスプレイ４２に表示させる。一方、判定手段４により異常と判定された場合には、画像処理サーバ１１から受信した画像データに基づいて最大吸気画像ＰNを高精細液晶ディスプレイ４２の画面の１つに表示させるとともに、別の画面には、局所変位ベクトル画像Ｑ1からＱN-1を呼吸相の変化の順に動画的に表示させる。
以上のように、本発明の第１の実施形態では、画像処理サーバ１１において、変位算出手段２が、胸部医用画像Ｐ1中の局所領域（Ｒ1[1,1]からＲ1[I,J]）の各々の位置を、画像Ｐ2から順に画像ＰNまでの対応する位置に合わせる位置合わせ処理を順次行って、各画像間の呼吸相の変化における各局所領域（Ｒ1[1,1]からＲ1[I,J]）の変位（Δ1[1,1]からΔ1[I,J]、・・・、Δn[1,1]からΔn[I,J]、・・・、ΔN-1[1,1]からΔN-1[I,J]）を算出し、散布度算出手段３が、各画像間における局所領域の変位毎に変位の差分の散布度Ｖ1からＶN-2を算出し、判定手段４が、散布度Ｖ1からＶN-2が閾値Thよりも大きい場合には、被検者の呼吸機能に異常があると自動判定する。したがって、判定結果に基づいて読影者の注意を喚起することが可能になり、異常箇所の見落としの減少や、読影時間の短縮による診断効率の向上に資する。
また、判定手段４において、その被検者の正常時の胸部医用画像を入力として算出された正常時散布度Ｖ0に基づいて閾値Thを設定するので、異なる被検者での個人差によって生じる判定誤差を減少させることが可能になり、さらなる判定精度の向上、診断精度や診断効率の向上に資する。
さらに、正常と判定された場合には、最大吸気画像ＰNをクライアントＰＣ４１の高精細液晶ディスプレイ４２に表示し、異常と判定された場合には、局所変位ベクトル画像生成手段５が、画像Ｐ1からＰNと、変位Δ1[1,1]からΔ1[I,J]、・・・、Δn[1,1]からΔn[I,J]、・・・、ΔN-1[1,1]からΔN-1[I,J]に基づいて、呼吸相の変化に伴う画像中の各局所領域の変位を視覚化した局所変位ベクトル画像Ｑ1からＱN-1を生成し、生成された画像がクライアントＰＣ４１の高精細液晶ディスプレイ４２に表示されるので、判定結果に応じて読影に供する画像を変えることにより、診断効率の向上が図られる。
なお、局所変位ベクトル画像生成手段５は、１からＮ−２までの自然数ｎについて、変位の差分（Δn＋1[1,1]−Δn[1,1]）、・・・、（Δn＋1[i,j]−Δn[i,j]）、・・・、（Δn＋1[I,J]−Δn[I,J]）の各々と所定の閾値ThΔとを比較し、変位の差分が所定の閾値ThΔよりも大きい局所領域では、変位ベクトルの視覚化の際に、他の領域とは異なる色にして局所変位ベクトル画像Ｑnを生成してもよい。これにより、読影者に対して、異常の有無だけでなく、異常箇所を指摘することが可能になるとともに、異常箇所の視認性が向上するため、さらなる診断効率の向上に資する。
また、局所変位ベクトル画像生成手段５は、１からＮ−１までの自然数ｎについて、各局所領域の中心画素の変位Δn[1,1]からΔn[I,J]に基づいて、２次元１０次のフィッティング多項式による近似処理を行うことによって、呼吸相ｎから呼吸相ｎ＋１への変化における画像Ｐn中のすべての画素の変位を求め、さらに、画像Ｐnの各画素を呼吸相ｎから呼吸相Ｎまでの各画素の変位の和に基づいてシフトさせる非線形歪変換処理（ワーピング）を行うことによって、画像Ｐ1からＰN-1中の被写体の位置を最大吸気画像ＰNの対応する位置に合わせたワープ画像Ｐ1″からＰN-1″を生成するようにし、１からＮ−１までの自然数ｎについて、ワープ画像Ｐn″中の各局所領域の中心画素の位置を始点として、変位Δn[1,1]からΔn[I,J]を変位ベクトルとして向きや大きさを視覚化した位置合わせ済み局所変位ベクトル画像Ｑn′を生成するようにしてもよい。このようにすれば、位置合わせ済み局所変位ベクトル画像Ｑ1′からＱN-1′は、画像中の各局所領域の位置が最大吸気画像ＰNの対応する位置に合わせられたものとなっているため、各局所領域の変位の視認性が向上し、診断精度や診断効率の向上に資する。
本発明の第２の実施形態では、被検者の呼吸運動における異なる呼吸相を表す３以上の胸部医用画像（以下、診断時画像）から選択された最大呼気画像Ｐ11、中間画像Ｐ12、最大吸気画像Ｐ13と、同じ被検者の正常時の３以上の胸部医用画像（以下、正常時画像）から選択された最大呼気画像Ｐ21、中間画像Ｐ22、最大吸気画像Ｐ23に基づいて、呼吸機能の異常を自動的に判定し、異常と判定された場合には、その異常部分を強調した強調画像を生成する機能を実現する。ここで、「正常時画像」とは、その被検者が過去に「異常なし」と診断された際に読影された胸部医用画像をいう。また、画像を表す符号Ｐtnの添字ｔは画像タイプで、「１」が診断時画像、「２」が正常時画像を意味し、添字ｎは呼吸相で、「１」が最大呼気状態、「３」が最大吸気状態、「２」がその中間状態を意味するものとする。以下、他の符号の添字についても同様の意味づけとなっている。
図１１は、この機能が実装される画像処理サーバ１１の論理的な構成とデータの流れとを表したブロック図である。ここでは、第１の実施形態と同様の機能を実現する手段は、同じ符号、名称としている。図に示すように、画像管理システム３０から複数の診断時画像と複数の正常時画像とを取得する取得手段１と、取得した診断時画像と正常時画像の各々について、画像中の肺野領域を検出し、検出された肺野領域の面積と、各画像から検出された肺野領域の面積の最大値との比率を算出し、算出された比率に基づいて、診断時の最大呼気画像Ｐ11、中間画像Ｐ12、最大吸気画像Ｐ13と、正常時の最大呼気画像Ｐ21、中間画像Ｐ22、最大吸気画像Ｐ23を選択する選択手段７と、取得された画像のうち診断時の最大呼気画像Ｐ11を複数の局所領域Ｒ11[1,1]、・・・、Ｒ11[i,j]、・・・、Ｒ11[I,J]に分割し、分割された局所領域Ｒ11[1,1]からＲ11[I,J]の各々の位置を、画像Ｐ12とＰ13の対応する位置に順次合わせる位置合わせ処理を行って、各画像間の呼吸相の変化における各局所領域（Ｒ11[1,1]からＲ11[I,J]）の変位（Δ11[1,1]からΔ11[I,J]、Δ12[1,1]からΔ12[I,J]）を算出するとともに、正常時の最大呼気画像Ｐ21を複数の局所領域Ｒ21[1,1]、・・・、Ｒ21[i,j]、・・・、Ｒ21[I,J]に分割し、分割された局所領域Ｒ21[1,1]からＲ21[I,J]の各々の位置を、画像Ｐ22とＰ23の対応する位置に順次合わせる位置合わせ処理を行って、各画像間の呼吸相の変化における各局所領域（Ｒ21[1,1]からＲ21[I,J]）の変位（Δ21[1,1]からΔ21[I,J]、Δ22[1,1]からΔ22[I,J]）を算出する変位算出手段２と、選択手段７によって選択された診断時と正常時の各画像中の注目領域Ｓ1n[1]からＳ1n[6]、Ｓ2n[1]からＳ2n[6]（n=1,2,3）を設定する注目領域設定手段８と、診断時と正常時の各々について、各画像間における局所領域の変位毎に変位の差分を算出し、注目領域毎に変位の差分の散布度（診断時散布度：Ｖ1[1]からＶ1[6]、正常時散布度：Ｖ2[1]からＶ2[6]）を算出する散布度算出手段３と、算出された診断時散布度Ｖ1[1]からＶ1[6]が、各々、正常時散布度Ｖ2[1]からＶ2[6]に基づいて設定された閾値Th[1]からTh[6]よりも大きい場合に被検者の呼吸機能に異常があると判定する判定手段４と、診断時の各画像について、呼吸機能に異常があると判定された注目領域を半透明色に着色した強調画像Ｔ11、Ｔ12、Ｔ13を生成する強調画像生成手段９と、呼吸機能が正常であると判定された場合には診断時の最大吸気画像Ｐ13をクライアントＰＣ４１に送信し、異常であると判定された場合には強調画像Ｔ11、Ｔ12、Ｔ13と診断時の各画像Ｐ11、Ｐ12、Ｐ13を送信する送信手段６から構成される。
上記の各手段は、ＣＤ−ＲＯＭ等の記憶媒体から画像処理サーバ１１にインストールされたプログラムが実行されることによって、画像処理サーバ１１のＣＰＵ、主記憶装置、外部記憶装置、入出力インターフェース、オペレーティングシステム等と連携して実現される。また、このプログラムによって、上記の各手段による処理順序が制御される。
次に、本発明の第２の実施形態において行われる処理の流れについて、第１の実施形態と相違する点を中心に説明する。
まず、画像撮影・読取システム２０における撮影や、画像管理システム３０における画像データの記憶については、第１の実施形態と同様である。
読影者が、第１の実施形態と同様に、クライアントＰＣ４１の操作により、患者ＩＤや撮影日時等の付随情報と、読影パターン（呼吸機能異常の自動判定）とを指定し、処理の実行を要求すると、その要求が画像処理サーバ１１に送信され、画像処理サーバ１１では、受信した読影パターンに基づいて、その読影に必要な画像を生成する画像処理プログラムが起動される。
図１２は、起動された画像処理プログラムによって実行される処理の流れを表すフローチャートである。
まず、取得手段１が、クライアントＰＣ４１から送信された付随情報に基づく検索条件による画像管理システム３０のデータベース検索で得られた３以上の診断時画像と３以上の正常時画像を取得し、画像処理サーバ１１の主記憶装置または外部記憶装置に一時的に記憶する（ステップb1）。なお、診断時画像も正常時画像も、最大呼気状態から最大吸気状態までの異なる複数の呼吸相を表したものとなっている。
次に、選択手段７が、診断時画像と正常時画像の各々について、最大呼気画像、中間画像、最大吸気画像を処理対象として選択する（ステップb2）。具体的には、公知の手法（例えば、特開２００３−００６６６１号公報）により各画像中の肺野領域を検出し、検出された肺野領域の面積を求め、診断時画像と正常時画像の各々について、肺野領域の面積が最小のものを最大呼気画像、最大のものを最大吸気画像とする。さらに、診断時画像と正常時画像の各々について、各画像中の肺野領域の面積と肺野領域の面積の最大値との比率を求め、比率が０．５に最も近い画像を中間画像とする。
次に、処理対象の画像が診断時のものか正常時のものかを識別する画像タイプに「１（診断時）」が設定され（ステップb3）、処理対象の画像の呼吸相の値が１（最大呼気）に設定される（ステップb4）。これにより、処理対象の診断時最大呼気画像Ｐ11が特定される。
以下、第１の実施形態と同様にして、変位算出手段２は、診断時最大呼気画像Ｐ11を読み込み（ステップb5）、局所領域Ｒ11[1,1]、・・・、Ｒ11[i,j]、・・・、Ｒ11[I,J]に分割する（ステップb6）。ここで、呼吸相を識別する添字ｎの値が１加算され、ｎ＝２となった後（ステップb7）、変位算出手段２は、診断時中間画像Ｐ12を読み込み（ステップb8）、局所領域Ｒ11[1,1]からＲ11[I,J]の各々の位置を、画像Ｐ12中の対応する位置に合わせる位置合わせ処理を行って、最大呼気状態から中間状態への呼吸相の変化における各局所領域Ｒ11[1,1]からＲ11[I,J]の変位を表す変位Δ11[1,1]からΔ11[I,J]を算出する（ステップb9）。
さらに、添字ｎの値が１加算されてｎ＝３となり（ステップB10,b7）、変位算出手段２は、上記と同様にして、中間状態から最大吸気状態への呼吸相の変化における各局所領域（Ｒ11[1,1]からＲ11[I,J]）の変位（Δ12[1,1]からΔ12[I,J]）を算出する（ステップb8, b9）。
この時点で添字ｎの値が３であるから変位算出手段２による処理が終了し（ステップb10）、注目領域設定手段８が、ｎ＝１，２，３について、診断時画像Ｐ1n中に注目領域を設定する（ステップb11）。具体的には、選択手段７で行われた肺野領域の検出結果を利用し、肺野領域内に右肺上葉を表す注目領域Ｓ1n[1]、右肺中葉を表す注目領域Ｓ1n[2]、右肺下葉を表す注目領域Ｓ1n[3]、左肺上葉を表す注目領域Ｓ1n[4]、左肺中葉を表す注目領域Ｓ1n[5]、左肺下葉を表す注目領域Ｓ1n[6]を設定することができる。なお、注目領域の設定とは、その注目領域に属する局所領域の範囲を決定することを意味する。また、[ ]を付した添字は各注目領域を識別するものである。なお、注目領域の数はいくつであってもよい。
散布度算出手段３は、第１の実施形態と同様にして、各画像間における局所領域の変位毎に変位の差分（Δ2[1,1]−Δ1[1,1]）、・・・、（Δ2[i,j]−Δ1[i,j]）、・・・、（Δ2[I,J]−Δ1[I,J]）を求める。そして、注目領域設定手段８で設定された注目領域毎に、これらの変位の差分の分散を散布度Ｖ1[1]からＶ1[6]として算出する（ステップb12）。
ここで、画像タイプｔの値が１加算されて、ｔ＝２（正常時）となり（ステップb13、b14）、正常時画像Ｐ21からＰ23について、ステップｂ４からｂ１２の処理が行われ、正常時散布度Ｖ2[1]からＶ2[6]の値が求まる。
さらに画像タイプｔの値が１加算されて、ｔ＝３となるので（ステップb13、b14）、次の判定処理に進む。
判定手段４は、１から６までの自然数ｋについて、散布度算出手段３によって求められた正常時散布度Ｖ2[k]を読み込み、αＶ2[k]＋β（α、βは、臨床データに基づいて得られた許容係数）の値を閾値Th[k]とし、診断時散布度Ｖ1[k]が閾値Th[k]よりも大きければ、被検者の注目領域Ｓ1n[k]（n=1,2,3)に異常があると判定し、診断時散布度Ｖ1[k]が閾値Th[k]以下であれば正常であると判定する（ステップb15）。
判定手段４によりすべての注目領域において正常と判定された場合には、送信手段６が、診断時最大吸気画像Ｐ13の画像データのみをクライアントＰＣ４１に送信する（ステップb16）。一方、判定手段４により異常と判定された場合には、強調画像生成手段９が、診断時の各画像について、呼吸機能に異常があると判定された注目領域を半透明色に着色した強調画像Ｔ11、Ｔ12、Ｔ13を生成する（ステップb17）。さらに、送信手段６が、生成された強調画像Ｔ11、Ｔ12、Ｔ13と診断時の各画像Ｐ11、Ｐ12、Ｐ13の画像データをクライアントＰＣ４１に送信する（ステップb18）。
クライアントＰＣ４１は、判定手段４により正常と判定された場合には、画像処理サーバ１１から受信した画像データに基づいて診断時最大吸気画像Ｐ13を高精細液晶ディスプレイ４２に表示させる。一方、判定手段４により異常と判定された場合には、画像処理サーバ１１から受信した画像データに基づいて、強調画像Ｔ11、Ｔ12、Ｔ13と診断時画像Ｐ11、Ｐ12、Ｐ13を高精細液晶ディスプレイ４２の画面に表示させる。
以上のように、本発明の第２の実施形態では、画像処理サーバ１１において、選択手段７によって、診断時画像から選択された最大呼気画像Ｐ11、中間画像Ｐ12、最大吸気画像Ｐ13と、同じ被検者の正常時画像から選択された最大呼気画像Ｐ21、中間画像Ｐ22、最大吸気画像Ｐ23に基づいて、注目領域設定手段８によって設定された注目領域Ｓ1n[1]からＳ1n[6]毎に（n=1,2,3）、正常時画像の対応する領域Ｓ2n[1]からＳ2n[6]との散布度の比較により呼吸機能の異常を自動的に判定する。したがって、判定結果に基づいて読影者の注意を喚起することが可能になり、異常箇所の見落としの減少や、読影時間の短縮による診断効率の向上に資する。
また、判定手段４において、その被検者の正常時の胸部医用画像を入力として算出された正常時散布度Ｖ2[1]からＶ2[6]に基づいて閾値Th[1]からTh[6]を設定するので、異なる被検者での個人差によって生じる判定誤差を減少させることが可能になり、さらなる判定精度の向上、診断精度や診断効率の向上に資する。
さらに、正常と判定された場合には、診断時最大吸気画像Ｐ13をクライアントＰＣ４１の高精細液晶ディスプレイ４２に表示し、異常と判定された場合には、強調画像生成手段９が、診断時の各画像について、呼吸機能に異常があると判定された注目領域を半透明色に着色した強調画像Ｔ11、Ｔ12、Ｔ13を生成し、生成された画像はクライアントＰＣ４１の高精細液晶ディスプレイ４２に表示されるので、判定結果に応じて読影に供する画像を変えることにより、診断効率の向上が図られる。
なお、上記の実施形態に第１の実施形態における局所変位ベクトル画像生成手段５と同様の処理を付加し、変位の差分（Δ2[1,1]−Δ1[1,1]）、・・・、（Δ2[i,j]−Δ1[i,j]）、・・・、（Δ2[I,J]−Δ1[I,J]）のうち、異常と判定された注目領域に属するものについて、変位の差分と所定の閾値ThΔとを比較し、変位の差分が所定の閾値ThΔよりも大きい局所領域では、変位ベクトルの視覚化の際に、他の局所領域とは異なる色にするようにしてもよい。これにより、注目領域レベルでの異常の有無に加えて、局所領域レベルでの異常の有無を読影者に対して指摘することが可能になるとともに、異常箇所の視認性が向上するため、さらなる診断効率の向上に資する。
また、強調画像生成手段９は、ｎ＝１，２について、各局所領域の中心画素の変位Δ1n[1,1]からΔ1n[I,J]に基づいて、２次元１０次のフィッティング多項式による近似処理を行うことによって、最大呼気から中間、中間から最大吸気への変化における画像Ｐ1n中のすべての画素の変位を求め、さらに、最大呼気画像Ｐ11および中間画像Ｐ12の各画素を最大吸気状態までの各画素の変位の和に基づいてシフトさせる非線形歪変換処理（ワーピング）を行うことによって、最大呼気画像Ｐ11および中間画像Ｐ12の被写体の位置を最大吸気画像Ｐ12の対応する位置に合わせたワープ画像Ｐ11″、Ｐ12″を生成するようにし、これらのワープ画像に注目領域を設定し、異常と判定された注目領域をワープ画像上で強調するようにしてもよい。このようにすれば、強調画像中の各画素の位置が最大吸気画像Ｐ12の対応する位置に合わせられたものとなっているため、各局所領域の変位の視認性が向上し、診断精度や診断効率の向上に資する。
上記の実施形態において、呼吸機能の異常の有無の判定が行われた後、診断時画像を画像管理システム３０の大容量外部記憶装置３２に保管する際、正常と判定された画像は高圧縮率でデータ圧縮し、異常と判定された画像は低圧縮率でデータ圧縮するか、または非圧縮としてもよい。
上記の実施形態において、呼吸機能に異常があると判定され、複数の呼吸相の画像を表示する際、時間方向（呼吸相の変化の方向）を奥行き方向と考えて、表示対象の画像を３次元的に表示してもよい。
（ａ）呼吸機能が正常な場合において、呼気状態から吸気状態への呼吸相の変化に伴う胸部医用画像中の局所領域毎の変位をベクトルで模式的に表した図（ｂ）このベクトルの向きと大きさの分布を表した図 （ａ）呼吸機能に異常がある場合において、呼気状態から吸気状態への呼吸相の変化に伴う胸部医用画像中の局所領域毎の変位をベクトルで模式的に表した図（ｂ）このベクトルの向きと大きさの分布を表した図 本発明の実施形態となる胸部画像診断支援システムの構成の概要を示した図 本発明の第１の実施形態による画像処理機能が実装される画像処理サーバの論理的な構成とデータの流れとを表したブロック図 本発明の第１の実施形態による画像処理の流れを表すフローチャート 変位算出手段と散布度算出手段で行われる処理を模式的に表した図 変位算出手段が行うグローバルマッチング処理を説明する図 変位算出手段が行うローカルマッチング処理を説明する図 変位算出手段が行うローカルマッチング処理により求めた各局所領域の中心画素のシフトの様子を模式的に表した図 局所変位ベクトル画像の一例を示す図 本発明の第２の実施形態による画像処理機能が実装される画像処理サーバの論理的な構成とデータの流れとを表したブロック図 本発明の第２の実施形態による画像処理の流れを表すフローチャート
符号の説明
１ 取得手段
２ 変位算出手段
３ 散布度算出手段
４ 判定手段
５ 局所変位ベクトル画像生成手段
６ 送信手段
７ 選択手段
８ 注目領域設定手段
９ 強調画像生成手段
１０ 画像処理システム
１１ 画像処理サーバ
２０ 画像撮影・読取システム
２１ ＣＲ装置
２２ Ｉ.Ｉ.
２３ ＦＰＤを搭載したＸ線画像デジタル撮影装置
３０ 画像管理システム
３１ 画像管理サーバ
３２ 大容量外部記憶装置
４０ 画像表示システム
４１ クライアントＰＣ
４２ 高精細液晶ディスプレイ
５０ ネットワーク

Claims (24)

  1. 被検者の呼吸運動における任意の３つの呼吸相を、該呼吸運動に伴う呼吸相の変化の順に第１、第２、第３の呼吸相とし、該第１、第２、第３の呼吸相の各々を表す胸部医用画像を第１、第２、第３の画像とした場合において、
    前記第１の画像を構成する複数の局所領域の各々の位置を、前記第２の画像中の対応する位置に合わせる位置合わせ処理を行って、前記第１の呼吸相から前記第２の呼吸相への変化における前記局所領域の変位を表す第１の変位を前記局所領域毎に算出し、
    前記位置合わせ処理後の前記局所領域の各々の位置を、前記第３の画像中の対応する位置に合わせる位置合わせ処理を行って、前記第２の呼吸相から前記第３の呼吸相への変化における前記局所領域の変位を表す第２の変位を前記局所領域毎に算出し、
    該局所領域毎の前記第１の変位と前記第２の変位の差分の散布度を算出し、
    該散布度が所定の閾値よりも大きい場合に前記被検者の呼吸機能に異常があると判定することを特徴とするコンピュータによる画像診断支援方法。
  2. 異なる呼吸相を表す複数の胸部医用画像の各々について、該画像中の肺野領域を検出し、検出された肺野領域の面積と、前記各画像から検出された肺野領域の面積の最大値または最小値との比率を前記各画像について算出し、算出された比率に基づいて前記第１、前記第２、前記第３の画像を選択し、選択された各画像に基づいて、前記第１および前記第２の変位の算出を行うものであることを特徴とする請求項１記載のコンピュータによる画像診断支援方法。
  3. 前記第１、前記第３の画像として、肺野領域の面積が最大となる最大吸気画像と、該面積が最小となる最大呼気画像とを選択することを特徴とする請求項２記載のコンピュータによる画像診断支援方法。
  4. 前記所定の閾値が、前記被検者の呼吸機能が正常な時の胸部医用画像に基づいて得られた前記散布度に基づいて設定された値であることを特徴とする第１項から第３項のいずれか１項記載のコンピュータによる画像診断支援方法。
  5. 前記第１、前記第２、前記第３の画像中の対応する位置に注目領域を設定し、
    設定された注目領域内の前記局所領域について、前記散布度を算出することを特徴とする第１項から第４項のいずれか１項記載のコンピュータによる画像診断支援方法。
  6. 前記第１、前記第２、前記第３の画像中の肺野領域を検出し、
    検出された肺野領域内に前記注目領域を設定することを特徴とする請求項５記載のコンピュータによる画像診断支援方法。
  7. 前記被検者の呼吸機能に異常があると判定された場合において、前記散布度が所定の閾値よりも大きい前記注目領域を他の領域と識別可能な態様にして表示することを特徴とする請求項５または６記載のコンピュータによる画像診断支援方法。
  8. 前記被検者の呼吸機能に異常があると判定された場合において、前記差分が所定の閾値よりも大きい前記局所領域における前記第１および／または前記第２の変位を他の局所領域における前記第１および／または前記第２の変位と識別可能な態様にして表示することを特徴とする第１項から第７項のいずれか１項記載のコンピュータによる画像診断支援方法。
  9. 被検者の呼吸運動における任意の３つの呼吸相を、該呼吸運動に伴う呼吸相の変化の順に第１、第２、第３の呼吸相とし、該第１、第２、第３の呼吸相の各々を表す胸部医用画像を第１、第２、第３の画像とした場合において、
    前記第１の画像を構成する複数の局所領域の各々の位置を、前記第２の画像中の対応する位置に合わせる位置合わせ処理を行って、前記第１の呼吸相から前記第２の呼吸相への変化における前記局所領域の変位を表す第１の変位を前記局所領域毎に算出するとともに、前記位置合わせ処理後の前記局所領域の各々の位置を、前記第３の画像中の対応する位置に合わせる位置合わせ処理を行って、前記第２の呼吸相から前記第３の呼吸相への変化における前記局所領域の変位を表す第２の変位を前記局所領域毎に算出する変位算出手段と、
    該局所領域毎の前記第１の変位と前記第２の変位の差分の散布度を算出する散布度算出手段と、
    該散布度が所定の閾値よりも大きい場合に前記被検者の呼吸機能に異常があると判定する判定手段とを備えたことを特徴とする画像診断支援装置。
  10. 異なる呼吸相を表す複数の胸部医用画像の各々について、該画像中の肺野領域を検出し、検出された肺野領域の面積と、前記各画像から検出された肺野領域の面積の最大値または最小値との比率を前記各画像について算出し、算出された比率に基づいて前記第１、前記第２、前記第３の画像を選択する選択手段をさらに備えたことを特徴とする請求項９記載の画像診断支援装置。
  11. 前記選択手段が、前記第１、前記第３の画像として、肺野領域の面積が最大となる最大吸気画像と、該面積が最小となる最大呼気画像とを選択するものであることを特徴とする請求項１０記載の画像診断支援装置。
  12. 前記所定の閾値が、前記被検者の呼吸機能が正常な時の胸部医用画像に基づいて得られた散布度に基づいて設定されたものであることを特徴とする第９項から第１１項のいずれか１項記載の画像診断支援装置。
  13. 前記第１、前記第２、前記第３の画像中の対応する位置に注目領域を設定する設定手段をさらに備え、
    前記散布度算出手段が、設定された注目領域内の前記局所領域について、前記散布度を算出するものであることを特徴とする第９項から第１２項のいずれか１項記載の画像診断支援装置。
  14. 前記第１、前記第２、前記第３の画像中の肺野領域を検出する検出手段をさらに備え、
    前記設定手段が、検出された肺野領域内に前記注目領域を設定するものであることを特徴とする請求項１３記載の画像診断支援装置。
  15. 前記被検者の呼吸機能に異常があると判定された場合において、前記散布度が所定の閾値よりも大きい前記注目領域を他の領域と識別可能な態様にして表示する手段をさらに備えたことを特徴とする請求項１３または１４記載の画像診断支援装置。
  16. 前記被検者の呼吸機能に異常があると判定された場合において、前記差分が所定の閾値よりも大きい前記局所領域における前記第１および／または前記第２の変位を他の局所領域における前記第１および／または前記第２の変位と識別可能な態様にして表示する手段をさらに備えたことを特徴とする第９項から第１５項のいずれか１項記載の画像診断支援装置。
  17. 被検者の呼吸運動における任意の３つの呼吸相を、該呼吸運動に伴う呼吸相の変化の順に第１、第２、第３の呼吸相とし、該第１、第２、第３の呼吸相の各々を表す胸部医用画像を第１、第２、第３の画像とした場合において、
    コンピュータを、
    前記第１の画像を構成する複数の局所領域の各々の位置を、前記第２の画像中の対応する位置に合わせる位置合わせ処理を行って、前記第１の呼吸相から前記第２の呼吸相への変化における前記局所領域の変位を表す第１の変位を前記局所領域毎に算出するとともに、前記位置合わせ処理後の前記局所領域の各々の位置を、前記第３の画像中の対応する位置に合わせる位置合わせ処理を行って、前記第２の呼吸相から前記第３の呼吸相への変化における前記局所領域の変位を表す第２の変位を前記局所領域毎に算出する変位算出手段と、
    該局所領域毎の前記第１の変位と前記第２の変位の差分の散布度を算出する散布度算出手段と、
    該散布度が所定の閾値よりも大きい場合に前記被検者の呼吸機能に異常があると判定する判定手段として機能させることを特徴とする画像診断支援プログラム。
  18. 前記コンピュータを
    異なる呼吸相を表す複数の胸部医用画像の各々について、該画像中の肺野領域を検出し、検出された肺野領域の面積と、前記各画像から検出された肺野領域の面積の最大値または最小値との比率を前記各画像について算出し、算出された比率に基づいて前記第１、前記第２、前記第３の画像を選択する選択手段としてさらに機能させることを特徴とする請求項１７記載の画像診断支援プログラム。
  19. 前記選択手段が、前記第１、前記第３の画像として、肺野領域の面積が最大となる最大吸気画像と、該面積が最小となる最大呼気画像とを選択するように前記コンピュータを機能させることを特徴とする請求項１８記載の画像診断支援プログラム。
  20. 前記所定の閾値が、前記被検者の呼吸機能が正常な時の胸部医用画像に基づいて得られた散布度に基づいて設定されたものであることを特徴とする第１７項から第１９項のいずれか１項記載の画像診断支援プログラム。
  21. 前記コンピュータを、
    前記第１、前記第２、前記第３の画像中の対応する位置に注目領域を設定する設定手段としてさらに機能させるとともに、
    前記散布度算出手段が、設定された注目領域内の前記局所領域について、前記散布度を算出するように機能させることを特徴とする第１７項から第２０項のいずれか１項記載の画像診断支援プログラム。
  22. 前記コンピュータを、
    前記第１、前記第２、前記第３の画像中の肺野領域を検出する検出手段としてさらに機能させるとともに、
    前記設定手段が、検出された肺野領域内に前記注目領域を設定するように機能させることを特徴とする請求項２１記載の画像診断支援プログラム。
  23. 前記コンピュータを、
    前記被検者の呼吸機能に異常があると判定された場合において、前記散布度が所定の閾値よりも大きい前記注目領域を他の領域と識別可能な態様にして表示させる手段としてさらに機能させることを特徴とする請求項２１または２２記載の画像診断支援プログラム。
  24. 前記コンピュータを、
    前記被検者の呼吸機能に異常があると判定された場合において、前記差分が所定の閾値よりも大きい前記局所領域における前記第１および／または前記第２の変位を他の局所領域における前記第１および／または前記第２の変位と識別可能な態様にして表示させる手段としてさらに機能させることを特徴とする第１７項から第２３項のいずれか１項記載の画像診断支援プログラム。
JP2005060028A 2005-03-04 2005-03-04 コンピュータによる画像診断支援方法および画像診断支援装置ならびにプログラム Withdrawn JP2006239195A (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP2005060028A JP2006239195A (ja) 2005-03-04 2005-03-04 コンピュータによる画像診断支援方法および画像診断支援装置ならびにプログラム
US11/367,300 US20060239530A1 (en) 2005-03-04 2006-03-06 Computer-aided image diagnosis

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2005060028A JP2006239195A (ja) 2005-03-04 2005-03-04 コンピュータによる画像診断支援方法および画像診断支援装置ならびにプログラム

Publications (1)

Publication Number Publication Date
JP2006239195A true JP2006239195A (ja) 2006-09-14

Family

ID=37046178

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2005060028A Withdrawn JP2006239195A (ja) 2005-03-04 2005-03-04 コンピュータによる画像診断支援方法および画像診断支援装置ならびにプログラム

Country Status (2)

Country Link
US (1) US20060239530A1 (ja)
JP (1) JP2006239195A (ja)

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2009078297A1 (ja) * 2007-12-19 2009-06-25 Konica Minolta Medical & Graphic, Inc. 動態画像処理システム
JP2009148336A (ja) * 2007-12-19 2009-07-09 Konica Minolta Medical & Graphic Inc 動態画像診断支援システム
JP2009153677A (ja) * 2007-12-26 2009-07-16 Konica Minolta Medical & Graphic Inc 動態画像処理システム
WO2009090894A1 (ja) * 2008-01-15 2009-07-23 Konica Minolta Medical & Graphic, Inc. 動態画像診断支援システム
JP2009273605A (ja) * 2008-05-14 2009-11-26 Konica Minolta Medical & Graphic Inc 動態画像診断支援システム
JP2010057795A (ja) * 2008-09-05 2010-03-18 Konica Minolta Holdings Inc 画像表示装置およびプログラム
JP2010158288A (ja) * 2009-01-06 2010-07-22 Konica Minolta Holdings Inc 画像表示装置およびプログラム
JP2012081037A (ja) * 2010-10-12 2012-04-26 Fujifilm Corp 診断支援装置、診断支援プログラムおよび診断支援方法
JP2012192255A (ja) * 2012-07-13 2012-10-11 Konica Minolta Holdings Inc 画像表示装置およびプログラム
WO2013099036A1 (ja) 2011-12-29 2013-07-04 トヨタ自動車株式会社 自動車用エアバッグシステム
JP2014516637A (ja) * 2011-04-28 2014-07-17 コーニンクレッカ フィリップス エヌ ヴェ 医療画像処理システム及び方法
JP2017196410A (ja) * 2016-04-26 2017-11-02 東芝メディカルシステムズ株式会社 医用画像処理装置及び医用画像診断装置
JP2018130264A (ja) * 2017-02-15 2018-08-23 コニカミノルタ株式会社 動態解析システム
JP2019171102A (ja) * 2019-05-30 2019-10-10 コニカミノルタ株式会社 画像処理装置、画像処理方法、およびプログラム

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8194946B2 (en) * 2005-07-28 2012-06-05 Fujifilm Corporation Aligning apparatus, aligning method, and the program
JP5537262B2 (ja) * 2010-05-28 2014-07-02 株式会社東芝 Ｘ線画像診断装置
JP6264949B2 (ja) * 2014-03-05 2018-01-24 富士ゼロックス株式会社 画像処理装置及びプログラム
JP6318739B2 (ja) * 2014-03-17 2018-05-09 コニカミノルタ株式会社 画像処理装置、およびプログラム
TWI524878B (zh) * 2014-09-16 2016-03-11 國立臺灣大學 疾病檢測方法及穿戴式裝置
JP6482250B2 (ja) * 2014-11-20 2019-03-13 キヤノン株式会社 画像処理装置、画像処理方法、およびプログラム
JP2017176202A (ja) * 2016-03-28 2017-10-05 コニカミノルタ株式会社 動態解析システム
JP6686733B2 (ja) * 2016-06-23 2020-04-22 コニカミノルタ株式会社 動態解析システム

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6980679B2 (en) * 1998-10-23 2005-12-27 Varian Medical System Technologies, Inc. Method and system for monitoring breathing activity of a subject
JP3697233B2 (ja) * 2002-04-03 2005-09-21 キヤノン株式会社 放射線画像処理方法及び放射線画像処理装置
WO2005112757A2 (de) * 2004-05-17 2005-12-01 Thomas Rupprecht Verfahren zur diagnose funktioneller lungenerkrankungen

Cited By (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8483456B2 (en) 2007-12-19 2013-07-09 Konica Minolta Medical & Graphic, Inc. Dynamic image processing system
JP2009148336A (ja) * 2007-12-19 2009-07-09 Konica Minolta Medical & Graphic Inc 動態画像診断支援システム
WO2009078297A1 (ja) * 2007-12-19 2009-06-25 Konica Minolta Medical & Graphic, Inc. 動態画像処理システム
US9125621B2 (en) 2007-12-19 2015-09-08 Konica Minolta, Inc. Dynamic image processing system
JP2009153677A (ja) * 2007-12-26 2009-07-16 Konica Minolta Medical & Graphic Inc 動態画像処理システム
JP2017124325A (ja) * 2008-01-15 2017-07-20 コニカミノルタ株式会社 動態画像解析装置
WO2009090894A1 (ja) * 2008-01-15 2009-07-23 Konica Minolta Medical & Graphic, Inc. 動態画像診断支援システム
JP2013039427A (ja) * 2008-01-15 2013-02-28 Konica Minolta Medical & Graphic Inc 動態画像診断支援システム及び胸部診断支援情報生成方法
JP5136562B2 (ja) * 2008-01-15 2013-02-06 コニカミノルタエムジー株式会社 動態画像診断支援システム
JP2009273605A (ja) * 2008-05-14 2009-11-26 Konica Minolta Medical & Graphic Inc 動態画像診断支援システム
JP2010057795A (ja) * 2008-09-05 2010-03-18 Konica Minolta Holdings Inc 画像表示装置およびプログラム
JP2010158288A (ja) * 2009-01-06 2010-07-22 Konica Minolta Holdings Inc 画像表示装置およびプログラム
JP2012081037A (ja) * 2010-10-12 2012-04-26 Fujifilm Corp 診断支援装置、診断支援プログラムおよび診断支援方法
JP2014516637A (ja) * 2011-04-28 2014-07-17 コーニンクレッカ フィリップス エヌ ヴェ 医療画像処理システム及び方法
WO2013099036A1 (ja) 2011-12-29 2013-07-04 トヨタ自動車株式会社 自動車用エアバッグシステム
JP2012192255A (ja) * 2012-07-13 2012-10-11 Konica Minolta Holdings Inc 画像表示装置およびプログラム
JP2017196410A (ja) * 2016-04-26 2017-11-02 東芝メディカルシステムズ株式会社 医用画像処理装置及び医用画像診断装置
JP2018130264A (ja) * 2017-02-15 2018-08-23 コニカミノルタ株式会社 動態解析システム
US10765391B2 (en) 2017-02-15 2020-09-08 Konica Minolta, Inc. Dynamic analysis system
JP2019171102A (ja) * 2019-05-30 2019-10-10 コニカミノルタ株式会社 画像処理装置、画像処理方法、およびプログラム

Also Published As

Publication number Publication date
US20060239530A1 (en) 2006-10-26

Similar Documents

Publication Publication Date Title
JP2006239195A (ja) コンピュータによる画像診断支援方法および画像診断支援装置ならびにプログラム
JP4493408B2 (ja) 画像読影支援方法及び装置並びにプログラム
US10417517B2 (en) Medical image correlation apparatus, method and storage medium
JP4911029B2 (ja) 異常陰影候補検出方法、異常陰影候補検出装置
US7620229B2 (en) Method and apparatus for aiding image interpretation and computer-readable recording medium storing program therefor
JP6318739B2 (ja) 画像処理装置、およびプログラム
JP2006014928A (ja) 画像表示方法および装置並びにプログラム
JP2011125568A (ja) 画像処理装置、画像処理方法、プログラム及び画像処理システム
KR20130018168A (ko) 의료 이미지 데이터 세트의 정합 품질의 시각화 방법 및 장치
JP2007159643A (ja) 画像処理装置および方法
US9530203B2 (en) Image processing apparatus and image processing method
US20180018772A1 (en) Dynamic analysis apparatus
JP4755863B2 (ja) 読影支援装置、読影支援方法、およびそのプログラム
JP6772123B2 (ja) 画像処理装置、画像処理方法、画像処理システムおよびプログラム
JP4493436B2 (ja) 画像読影支援方法および装置ならびにプログラム
JP2005124617A (ja) 医用画像診断支援システム
RU2568262C2 (ru) Соотнесение полученных изображений с объектами
JP5364334B2 (ja) 医用画像処理装置
JP2010057795A (ja) 画像表示装置およびプログラム
JP5100041B2 (ja) 画像処理装置及び画像処理プログラム
JP7349870B2 (ja) 医用画像処理装置、断層撮影装置、医用画像処理方法及びプログラム
US9918685B2 (en) Medical image processing apparatus and method for medical image processing
JP6852545B2 (ja) 画像表示システム及び画像処理装置
JP2019180899A (ja) 医用画像処理装置
JP2006247293A (ja) 画像処理方法および画像処理装置ならびにプログラム

Legal Events

Date Code Title Description
A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A712

Effective date: 20061208

A300 Application deemed to be withdrawn because no request for examination was validly filed

Free format text: JAPANESE INTERMEDIATE CODE: A300

Effective date: 20080513