JP2004305486A - 画像処理装置及び方法及びシステム - Google Patents

Abstract

【課題】動態画像に基づいて診断に有効な画像を生成すること。
【解決手段】時系列的に連続して取得された複数のＸ線画像の少なくともいくつかの画像を表示し、それらの各々に着目点を示すマークを指示可能なインターフェースと呈示する（ステップＳ２０１）。このインターフェースを介して、ユーザが複数のフレーム画像の所望の位置にマークを設定すると、複数のＸ線画像が時間軸方向に並べられた３次元空間において上記設定された着目点を含む曲面が生成される（ステップＳ２０２）。そして、時間軸に沿って並べられた複数のＸ線画像の各々より、上記曲面との交線に沿った画素値群を取得し、この画素値群に基づいて新たな画像を生成する（ステップＳ２０３、Ｓ２０４）。この新たな画像は例えば平面上に投影された状態でディスプレイに表示される（ステップＳ２０５）。
【選択図】 図１

Description

【０００１】
【発明の属する技術分野】
本発明は画像処理技術に関し、特に医療分野における放射線撮影画像の処理に好適な画像処理技術に関する。
【０００２】
【従来の技術】
Ｘ線に代表される放射線の物質透過能力を用いて、その透過強度分布を画像化する技術は、近代医療技術発展の基本となるものである。Ｘ線の発見以来、その強度分布の画像化については、Ｘ線強度分布を蛍光体により可視光に変換した後、銀塩フィルムで潜像を作り現像するという方法が採られてきた。近年、Ｘ線画像をデジタル化する際に輝尽性蛍光体を用い、Ｘ線照射による輝尽性蛍光体上の蓄積エネルギ分布としての潜像をレーザ光で励起して読み出し、デジタル画像化する、いわゆるイメージングプレートを用いる方法も一般化してきた。さらに、半導体技術の進歩により人体の大きさをカバーできる大判の固体撮像素子、いわゆるフラットパネルディテクタも開発され、潜像をつくることなく直接にＸ線画像をデジタル化し、効率のよい診断が行えるようになって来た。
【０００３】
また一方、光電子増倍管（イメージインテンシファイア）に代表される高感度の撮像素子により微弱なＸ線による蛍光を画像化し、人体内部の動態を観察することも可能であり、一般に用いられてきている。最新のフラットパネルディテクタは、そのようなイメージインテンシファイアにも匹敵する感度を持ち、人体の広範囲における動態を撮影することも可能になってきている。
【０００４】
一般に、医療用のＸ線撮影で最も有効であるのは人体の胸部撮影である。腹部を含む胸部の広範囲を撮影すれば、肺疾患を含む多くの疾病の発見に役立つため、通常の健康診断では胸部Ｘ線撮影は不可欠なものになっている。また、近年、健康診断のために撮影された膨大な量の胸部Ｘ線画像を効率よく診断するため、胸部デジタルＸ線画像に対して計算機を用いて画像解析を行い、医師の初期診断を補助するいわゆる計算機支援診断（ＣＡＤ：Ｃｏｍｐｕｔｅｒ Ａｉｄｅｄ Ｄｉａｇｎｏｓｉｓ）も実用化されつつある。
なお、尚、放射線画像に現われる異常陰影を検出する異常陰影検出装置が特許文献１及び特許文献２において知られている。
【０００５】
【特許文献１】
特許第２５８２６６５号公報
【特許文献２】
特許第２５８２６６６号公報
【０００６】
【発明が解決しようとする課題】
さて、健康診断で胸部Ｘ線撮影等により何らかの疾病の疑いがあるとの所見が得られた場合は、いわゆる精密診断により確定診断が行われる。このような精密診断には、多くの場合ＣＴ、ＭＲスキャンが行われる。このＣＴスキャンやＭＲスキャンは、通常のＸ線撮影と比して数倍の診断コストがかかるものである。このような精密診断により初期の診断が誤診であり何の疾病も無いという場合も多く見られるが、このような誤診は無駄な医療費を費やすものであり、医療費高騰の一因ともなっている。これを防ぐには、初期診断である健康診断の正確度を向上させることが重要である。
【０００７】
コスト上昇を抑えて診断の正確度を向上させる方法としては、前述の大判のフラットパネルディテクタを用いて、呼吸などにより動態を示す胸部動画像を取得し、動態観察を行うことが有効である。動態観察では従来の静止画を観察する場合とは異なり、画像を時系列で切り換えて観察することになる。しかしながら、このような動態観察では、静止画観察に比べて観察者の主観に依存するところが多く、あいまいな点が残りやすい。また、多大な観察時間も費やされることになるため、診断効率も悪い。
【０００８】
本発明は上記の課題に鑑みてなされたものであり、動態画像に基づいて診断に有効な画像を生成することを１つの目的とする。
【０００９】
【課題を解決するための手段】
上記の目的を達成するための本発明による画像処理装置は以下の構成を備える。すなわち、
対象物の動態画像を構成する複数の画像を時間軸に沿って並べた時空間３次元空間に設定された曲面と前記複数の画像の各々との交線に沿って並ぶ画素列を取得する取得手段と、
前記取得手段により前記複数の画像から取得された複数の画素列に基づいて画像を生成する生成手段とを備える。
また、本発明の他の態様による画像処理装置は以下の構成を備える。すなわち、
対象物の動態画像を構成する複数の画像の各々から、設定された矩形領域の一辺の方向に並ぶ複数画素を統合することにより画素列を生成する統合手段と、
前記統合手段により前記複数の画像から生成された複数の画素列に基づいて画像を生成する生成手段とを備える。
【００１０】
また、本発明によれば、上記画像処理装置と撮像装置を備えたシステムが提供される。
また、本発明によれば、上記画像処理装置によって実行される画像処理方法が提供される。
【００１１】
【発明の実施の形態】
以下、添付の図面を参照して本発明の好適な実施形態について説明する。
【００１２】
［第１実施形態］
上述した課題を解決するために、本実施形態では、時系列的に連続して取得された動態画像（本実施形態では呼吸動態画像）を、時間方向を奥行きに取った「時空間３次元画像」として扱い、その横断像、矢状断像（ｓａｇｉｔｔａｌ ｐｌａｎｅ）を作成し、観察可能に表示する。これにより、動態画像の、ひいては医療用動画像の静的な観察および解析を可能にする。
【００１３】
図１は第１実施形態による呼吸動態撮影装置の構成を模式的に示すブロック図である。図１において、１はＸ線強度を画像化するＸ線画像センサであり、連続的な画像取得が可能である。Ｘ線画像センサ１としてはフラットパネルディテクタを用いることができる。２はＸ線発生装置でありＸ線を発生し、被検体に照射する。３は制御装置であり、Ｘ線画像センサ１およびＸ線発生装置２を制御し、Ｘ線発生と画像取得の適切な同期をとる。４は被写体（被検体）としての人体を示している。
【００１４】
５は連続画像取得ユニットであり、Ｘ線画像センサ１から出力される画像データを連続的に取得する。６は動画表示ユニットであり、連続画像取得ユニット５によって連続的に取得される人体の胸部画像を実時間でディスプレイ上に表示する。７はストレージ装置であり、連続画像取得ユニット５で取得された画像を格納する。なお、動画像表示ユニット６はストレージ装置７から随時データを引き出し、順次参照することが可能である。すなわち、動画像表示ユニット６は、実時間での動画的表示や、録画再生のように連続的に取得された画像を後に動画的表示させることができる。
【００１５】
８は時空間３次元画像作成ユニットであり、連続画像取得ユニット５により連続的に取得された時系列画像から時間方向を奥行きとして３次元画像を作成する。断面作成指示ＵＩユニット９は、時空間３次元画像作成ユニット８によって作成された３次元画像において、表示すべき断面を指示するためのユーザインターフェース（ＵＩ）を提供する。１０は断面作成ユニットであり、時空間３次元画像作成ユニット８によって生成された時空間３次元画像の、時間方向（奥行き方向）を含む任意の断面（この断面は断面作成指示ＵＩユニット９により指定される）を作成する。１１は断面表示ユニットであり、断面作成ユニット１０によって生成された断面をディスプレイに表示する。なお、動画像表示ユニット６及び断面表示ユニット１１が表示に用いるディスプレイは共通のものであってもよい。
【００１６】
本システムで得られる呼吸動態画像は複数の人体胸部Ｘ線画像である。以下、実施形態ではこのような動態画像を構成する複数の画像の各々をフレーム画像という。図２では仮に５つのフレーム画像が存在するものとして、各タイミングＴ１〜Ｔ５で得られたフレーム画像を並べて示している。図２のように動画像のフレーム画像を並べて示すことは各フレーム画像を精査する場合は有効である。しかしながら、一般に多数のフレーム画像を並べなければならないため、診断には不適切である。したがって、通常はフレームサイズの画面にフレーム画像を連続して切り換え表示し、動画像として観察する。このように動画像として観察する方法が、観察者に直感的に動態を把握させるのに有利であることは明らかであるが、動画観察には時間を要する。また動画像の場合には客観的な診断記録を残すのが困難である。すなわち、動画像は多くの情報をもっているが、動的に表示して観察しなければならないため、書類などに客観的診断記録を残すことが困難である。これに対し、静止画であれば、従来どおり誰でも、いつでも観察できるし、マーキングなども容易なため、客観的な記録をのこすことが可能である。
【００１７】
そこで、本実施形態では、動的な様子を客観的に示すために、たとえば図２のＡの横ライン上の各フレーム画像の画素値を取り出して画素列を得て、取り出された画素列を取り出し元のフレーム画像の時系列順で並べ直し、図３に示すように一つの画像を形成する。図３では下からＴ１、Ｔ２、…、Ｔ５というように並べ直して画像を形成している。このような図３の画像によれば、図２のＳ１で示される陰影の時間的な動き・変動が１つの画像として得られることになる。
【００１８】
上記のような画像によれば、たとえば、図２のＳ１が腫瘍陰影であるとして、Ｓ１が向かって左側の胸壁に癒着しているのかどうかという解剖学的判断が容易に行える（図２、図３の例では左胸壁に癒着していないことがわかる）。このような判断は１枚の静止画像からでは、撮影の方向・タイミングなどさまざまな要因があって、通常は判断し難いものである。図３の画像のように、あるタイミングでの撮影で、癒着がないことが確認できることは医学的に重要であり、図３に示すような画像表示の診断への貢献度は高い。また、さらに、このような１枚の画像を形成すれば、診断した様子の記録を容易に残すことができる。
【００１９】
同様に、図２のＢで示される縦ライン上の各フレーム画像の画素列を時系列順に並べ直し、図４に示すような１枚の画像に表し、観察、記録することで、陰影Ｓ１とＳ２が１つの陰影ではなく、離れた２つの陰影であることが確認できる。なお、図４では、左から右へＴ１からＴ５の順に各フレーム画像から取得された画素列を並べている。
【００２０】
以上説明した図３、図４に示される画像は、時系列的に連続した複数の画像を時間方向を奥行きに取った「時空間３次元画像（図５、図６により詳述する）」として扱うことにより得られた３次元画像の横断面像及び縦断面像に相当する。以下、図３、図４のような画像を時系列断面画像と称する。
【００２１】
通常、上述したような癒着等の診断はＣＴスキャン等による人体の３次元的な情報からのみ得られるものであったが、本方式のような時系列断面画像を観察することで、それに相当するような診断も可能になる。
【００２２】
なお、上記では図示及び説明の都合で５つのフレーム画像を示したが、このような少ないフレーム数ではなく、少なくとも３０フレーム以上ある動画像を対象とするものとする。また、指定されたライン上の１画素列を配置する際に、画素列間に１つ又は複数の補間画素列を挿入して図３或いは図４のような時系列断面画像を形成してもよい。また、ライン上の画素列を含む複数画素列を各フレームより抽出し、これらを用いて、図３、図４に示すような時系列断面画像を並べて生成、表示するようにしてもよい。
【００２３】
断面作成指示ＵＩユニット９では、例えば、図２に示すように複数枚のフレーム画像を並べて表示し、ラインＡの位置をポインティングデバイス或いはキーボードからの指示によって移動させ、所望の断面位置を指定するようにすればよい。また、ラインＢを指定する場合には、図２に示すように各フレーム画像にラインＢを示す表示を行なえばよい。この場合例えば１つのフレーム画像中のラインＢ（例えばＴ１のラインＢ）を移動すると、他のフレーム画像中のラインＢも同じ量だけ移動する。或いは、ラインＢの指定時には各フレーム画像を縦方向に並べて表示し、ラインＢを設定させるようにしてもよい。なお、このようなＵＩでは、特に時空間３次元画像作成ユニット８によって作成された時空間３次元画像は用いなくてよいので、時空間３次元画像作成ユニット８は省略することも可能である。
【００２４】
さて、以上のようにして断面を指定するＵＩを構成することも可能であるが、本実施形態では、さらに、このような断面を医師もしくは機器操作者が画像を観察しながら容易に形成できるＵＩを提案する。このＵＩでは、時空間３次元画像作成ユニット８によって作成された３次元画像を表示させ、所望の断面をユーザに指定させる。
【００２５】
図５は、時空間３次元画像作成ユニット８により作成された時空間３次元画像を表し、通常の２次元画像のＸ，Ｙ方向に加えて、奥行きに時間軸であるＴ方向が加えられて３次元空間を構成している。したがって、図５に示される空間の１つの点は（ｘ，ｙ，Ｔ）の３つの値で表現できるものになる。３次元画像の表示方法としては、レイトレーシングによる半透明な物体としての表現、後ろの画像を前面の画像で覆い隠すような表現方法などがある。
【００２６】
図１の断面作成指示ＵＩユニット９では、図５のような３次元画像が表示され、ユーザはこの表示を見ながら自分の所望する断面像を得るため、断面の指示を行う。たとえば、断面の方程式として
Ｔ＝Ｃ ； Ｃは定数 （数式１）
という平面を選択し、Ｃを可変にすれば、通常の各フレーム画像が観察可能になる。
【００２７】
また、断面の方程式として
Ｙ＝Ａ ； Ａは定数 （数式２）
という平面を選択すれば、図６のＡで示すようなＸ−Ｔ平面に平行な平面の断面像が得られる。すなわち、図３で示したような断面像が得られて、観察・診断・記録が可能になる。
【００２８】
また、断面の方程式として
Ｘ＝Ｂ ； Ｂは定数 （数式３）
という図６のＢで示すような平面を選択すれば、Ｙ−Ｔ平面に平行な平面の断面像が得られる。すなわち、図４で説明したような断面画像が得られて観察・診断・記録が可能になる。
【００２９】
一般にユーザは、任意の平面が選択でき、その平面は下記のような一般的な方程式、
Ｔ＝Ｋｘ・Ｘ＋Ｋｙ・Ｙ＋Ｋ１；Ｋｘ，Ｋｙ，Ｋ１は定数 （数式４）
もしくは、同値でありながら、別の表現で表記できる２つの方程式、
Ｘ＝Ｋｙ・Ｙ＋Ｋｔ・Ｔ＋Ｋ２；Ｋｙ，Ｋｔ，Ｋ２は定数 （数式５）
Ｙ＝Ｋｘ・Ｘ＋Ｋｔ・Ｔ＋Ｋ３；Ｋｘ，Ｋｔ，Ｋ３は定数 （数式６）
で表される。
【００３０】
操作者は、３次元画像を観察しながら、たとえばマウス操作により画面に表示された平面（例えば、図６の平面ＡやＢ）を平行、回転移動させることにより、所望の断面を選択できる。
【００３１】
図７は本実施形態の呼吸動態撮影装置における時系列断面画像の生成処理を示すフローチャートである。
【００３２】
まず、ステップＳ１１において、時空間３次元画像作成ユニット８はストレージ装置７に格納された時間的に連続する複数のＸ線画像群を用いて図５に示すような時空間３次元画像を生成する。断面作成指示ＵＩユニット９はこれを不図示のディスプレイに表示する。なお、時空間３次元画像の生成においては、各フレーム画像間に補間フレーム画像を生成して挿入してもよい。
【００３３】
ステップＳ１２において、ユーザはこの３次元表示を参照しながら、当該３次元画像の所望の断面を指定するべく切り出し用平面を設定する。切り出し用平面は、図６に示したＡ及びＢの平面の他に、図８に示すようなＸ、Ｙ、Ｔの各軸に対して傾いた平面を設定することも可能である。切り出し用平面の移動、回転等の操作はマウスにより行なうものとするが、このような３次元画像と平面の表示処理は当業者には明らかであろう。
【００３４】
ステップＳ１２による切り出し用平面の設定により、数式４〜６のいずれかにおける各定数が決定され、処理はステップＳ１３へ進む。ステップＳ１３では、各フレーム画像から、フレーム画像と切り出し用平面との交線を取得し、当該交線に沿った画素列を取得する。この様子を図９に示す。図９では、切り出し用平面９１と各フレーム画像（Ｔ１〜Ｔ４）との交線として、Ｌ１、Ｌ２、Ｌ３、Ｌ４が示されている。なお、これらの交線はステップＳ１２で決定された切り出し用平面の式を用いることで取得できる。
【００３５】
続いてステップＳ１４において、切り出し元のフレーム画像の時系列順に従ってステップＳ１３で取得した画素列を配置して時系列断面画像（図９の９２）を生成する。なお画素列の配置としては、時系列順に下→上、上→下、左→右、右→左等をユーザが設定可能とすればよい。なお、上述したように時系列断面画像の生成においては、取得された画素列間に補間画素列を生成して挿入してもよい。
【００３６】
以上のようにして生成された時系列断面画像９２は断面表示ユニット１１によってディスプレイ上に表示される。
【００３７】
［第２実施形態］
断面の表現方法は、第１実施形態で示した図３や図４の形態に限られるものではない。例えば、図１０に示すような形態を用いることも可能である。図１０では、時系列断面として１０３が指定された場合の第２実施形態による時系列断面画像の表示方法を説明する図である。フレームＴ１の指定された断面より上側の部分画像１０１とフレームＴＮの指定された断面より下側の部分画像１０２を離間させて表示し、部分画像１０１と部分画像１０２をつなぐように断面画像１０３を配置、表示する。
【００３８】
図１１の（ａ）は図３に相当する横断面の画像を図１０に示した方法に従って表現したものである。また、図１１の（ｂ）は図４に相当するような縦断面を図１０に示した方法に従って表現したものである。このような３次元的な表示により、観察者は実際にどの断面の時間的変化を観察しているのかをより明確に把握しながら観察することができ、この画像を記録することにより、より客観的な診断記録の保存が可能になる。
【００３９】
［第３実施形態］
以上第１及び第２実施形態では、動画像観察に頼らず、静止画的に動態を診断する方法として、時空間３次元像を作成し、その平面断面を観察する方法を説明した。しかしながら、呼吸による人体部位の動きにより、平面断面で診断することが困難な場合もある。そこで、第３実施形態では、取得された動態画像（例えば、呼吸動態画像、胸部動態画像ともいう）を、時間方向を奥行きに取った時空間３次元画像として扱い、注目部位を含む任意の曲面を作成し、その曲面断面画像を表示する。このような曲面断面画像の観察によれば、例えば呼吸による人体の動きに良好に対応でき、より効果的に診断を補助することができる。
【００４０】
なお、第３実施形態においても、図１に示した構成の呼吸動態撮影装置を用いるものとする。以下、第３実施形態による画像処理手順について図１２乃至図１５を参照して説明する。
【００４１】
図１２は第３実施形態の呼吸動態撮影装置における時系列断面画像の生成処理を示すフローチャートである。まずステップＳ２０１において、複数のＸ線画像の所望の位置マークを指示するためのインターフェースを呈示する。例えば、図１３に示すように各タイミングＴ１〜Ｔ５で得られたフレーム画像を並べて表示する。なお表示するフレーム画像は、全てを表示してもよいし、それらのうちから選択されたフレーム画像群（例えば所定間隔で抽出されたフレーム画像）であってもよい。図１３は、呼吸動態の各フレーム画像を模式的に示したものであるが、ここで符号Ｓ３，Ｓ４で示す陰影の時間的変動を追跡したいという要求が出たとする。操作者は、図１４に示すように各フレーム画像において、目的とする陰影にマーク（×）をつける。この操作は、腫瘍陰影を自動的に抽出するプログラム（上述の特許文献１及び２を参照）を用いて行うことも可能である。また、目的とする陰影を囲むように指示してもよい。この場合、例えば指示された領域の重心位置を上記のマーク位置として扱う。なお、ステップＳ２０１及びＳ２０２では、全フレームを表示してもよいし、適当な間隔で間引いて表示してもよい。或いは、複数の画像を頁切り換えのように１枚ずつ切換表示するようにしてもよい。
【００４２】
このように複数のフレーム画像においてマーク（×）を付けられ、その総数をｎ個とすると、時系列３次元空間（ｘ，ｙ，Ｔ）中の座標点の集合｛（ｘｉ，ｙｉ，Ｔｉ）；ｉ＝１．．．ｎ｝が得られることになる。そして、ステップＳ２０２において、ｎ点を補間する曲面を作成する。曲面を作成する方法としては、多項式スプラインによるもの、最小２乗による近似多項式平面によるもの、ラグランジュ多項式によるものが考えられ、どれを用いても構わない。ここでは、簡単な一例としてラグランジュ多項式平面を用いる。ラグランジュ多項式平面は、下式の（数式７）に示されるラグランジュ係数ｌ_ｊ（ｘ，Ｔ）に対して、曲面方程式ｙ（ｘ，Ｔ）を（数式８）の方程式で表される関数で補間し、曲面を得る。
【００４３】

【００４４】
次に、ステップＳ２０３において、ステップＳ２０２で求めた曲面上の画素値を求め、ステップＳ２０４で、図１５の符号４１で示すような濃度を持つ曲断面画像を形成する。より具体的には、各画像と曲面との交線を求め、当該交線に沿った画素列を各フレーム画像より取得する。そして、取得した画素列を取得元のフレーム画像の時系列に従って配置し、曲断面画像を得る。なお、第１実施形態で述べたように、各画素列間に補間画素列を挿入してもよい。或いは、時間軸を奥行き方向とする時空間３次元画像を生成しておき、この時空間３次元画像の上記曲面による曲断面を取得することで曲断面画像を構成してもよい。
【００４５】
次に、ステップＳ２０５において、ステップＳ２０４で生成した曲断面画像を表示する。曲断面を表示する場合は、２次元画像にする必要があるので、図１５の（ａ）において符号４２で示すように、平面上へ投影した画像、もしくは図１５の（ｂ）に示すように平面上に広げた（展開した）画像を形成する。なお、曲面を平面的に表示する方法は周知の各種方法を適用することができる。或いは、３次元的にそのまま画像表示する（例えば図１５の４１のような形態で表示する）ことも可能である。このような画像を観察・解析することにより所望部位の動態の様子を客観的に判断できるようになる。
【００４６】
以上のような第３実施形態によれば、以下のような効果が得られる。すなわち、胸部動態画像は呼吸による部位の移動を基本としているが、その部位の移動は上記第１実施形態で示したような平面によって表されるような単純な動きではない場合も多く、そのような場合は第１、第２実施形態のような平面断面では把握し難い。これに対して、第３実施形態では、任意の曲面の断面画像を形成し、その断面を曲面的または平面的に観察できるので、移動する部位を的確に追跡して診断を行うことが可能になる。
なお、曲面の観点からいえば、第１及び第２実施形態では平面によって時空間３次元画像を切り取ったが、平面との交線が２次曲線や３次曲線になるような曲平面とすることも可能であることは明らかであろう。
【００４７】
［第４実施形態］
時系列に取得された複数の画像から、腫瘍等の動きを呈示するための静止画像として、第１及び第３実施形態では所望の平面或いは曲面によって得られる時系列３次元画像の断面を表示した。第４実施形態では時系列に取得された複数の画像から、腫瘍等の動きを呈示するための静止画像として、以下に説明する統合処理を採用する。
【００４８】
図１６は第４実施形態の呼吸動態撮影装置における時系列断面画像の生成処理を示すフローチャートである。まずステップＳ３０１において、複数のＸ線画像に所望の範囲を指定させるためのインターフェースを呈示する。例えば、図１７に示すようなインターフェース画面を表示する。図１７は図１３と同じフレーム画像であるが、第３の実施形態では部位の移動を追跡して曲面を形成したのに対し、第４実施形態では、移動範囲の画素値を統合した画像を作成する。そのために、第４実施形態のインターフェースでは、４５に示すような範囲をユーザに設定させる。
【００４９】
ステップＳ３０２では、範囲の指定方向（図１７ではＹ方向）に沿って並ぶ画素の画素値を統合する。そして、ステップＳ３０３では、得られた統合値をＹ方向と直交するＸ方向に並べて、各画像毎に統合値列を形成する。画素値の統合の方法としては、最大値投影、最小値投影、平均、メディアンフィルタなどさまざまな方法があり、画像の特性、目的により使い分ける必要がある。図１７では符号４５で示す範囲の画素値を統合して、Ｘ方向に並ぶ統合値列を生成する。例えば、画像において腫瘍部分の濃度が高ければ、最大値投影を用いて統合することにより腫瘍の存在部分の濃度が高い統合値列が生成されることになる。
【００５０】
ステップＳ３０４では、ステップＳ３０３で得られた画素列（統合値列）を、取得元のフレーム画像の時系列順に配置して、統合画像を形成する。そして、ステップＳ３０５において、生成した統合画像を表示する。
【００５１】
以上の操作は、図１８の４６で示す時空間３次元画像に対して、４７で示す部分の画素値を上下方向に統合することに相当する。したがって、図１７に示すようなインターフェースの代りに、断面作成指示ＵＩユニット９は、図１８に示す時系列３次元画像と切り出し部分を示す直方体（平行六面体）４７を表示し、直方体４７の操作（直方体の位置、回転角、厚み等を変更すること）によって所望の部位を指定するようにしてもよい。なお、図１８の時空間３次元画像は時空間３次元画像作成ユニット８で作成される。４８で示すように、部位の移動に関わり無く統合された画像が得られ、この画像を観察・解析することにより所望部位の動態の様子を客観的に判断できるようになる。
【００５２】
以上のように第４実施形態によれば、時空間３次元像中の注目する部位の移動範囲を最大値、最小値、平均値、中間値、合計値などの統計量（統計値）で代表させることにより、統合画像を作成し、その画像を観察・解析することで、呼吸動態画像、ひいては医療用動画像の静的な観察および解析が可能になる。
【００５３】
なお、本発明の目的は、前述した実施形態の機能を実現するソフトウェアのプログラムコードを記録した記憶媒体を、システムあるいは装置に供給し、そのシステムあるいは装置のコンピュータ（またはＣＰＵやＭＰＵ）が記憶媒体に格納されたプログラムコードを読出し実行することによっても、達成されることは言うまでもない。
【００５４】
この場合、記憶媒体から読出されたプログラムコード自体が前述した実施形態の機能を実現することになり、当該プログラムコード、及び当該プログラムコードを記憶した記憶媒体は本発明を構成することになる。
【００５５】
プログラムコードを供給するための記憶媒体としては、例えば、フレキシブルディスク，ハードディスク，光ディスク，光磁気ディスク，ＣＤ−ＲＯＭ，ＣＤ−Ｒ，磁気テープ，不揮発性のメモリカード，ＲＯＭなどを用いることができる。
【００５６】
また、コンピュータが読出したプログラムコードを実行することにより、前述した実施形態の機能が実現されるだけでなく、そのプログラムコードの指示に基づき、コンピュータ上で稼働しているＯＳ（オペレーティングシステム）などが実際の処理の一部または全部を行い、その処理によって前述した実施形態の機能が実現される場合も、本発明の実施態様に含まれることは言うまでもない。
【００５７】
さらに、記憶媒体から読出されたプログラムコードが、コンピュータに挿入された機能拡張ボードやコンピュータに接続された機能拡張ユニットに備わるメモリに書込まれた後、そのプログラムコードの指示に基づき、その機能拡張ボードや機能拡張ユニットに備わるＣＰＵなどが実際の処理の一部または全部を行い、その処理によって前述した実施形態の機能が実現される場合も、本発明の実施態様に含まれることは言うまでもない。
【００５８】
【発明の効果】
以上説明したように本発明によれば、動態画像に基づいて診断に有効な画像を生成することができる。
【図面の簡単な説明】
【図１】本発明を実施した第１の実施形態のブロック図である。
【図２】胸部動態画像の各フレームを模式的に表した図である。
【図３】特定横ラインの時間変動を画像化した図である。
【図４】特定縦ラインの時間変動を画像化した図である。
【図５】時空間３次元画像の模式図である。
【図６】時空間３次元画像の平面断面を示す図である。
【図７】第１実施形態の呼吸動態撮影装置における時系列断面画像の生成処理を示すフローチャートである。
【図８】切り出し用平面による時空間３次元画像の断面の設定を説明する図である。
【図９】各フレームから時系列断面上の画素列を取得し、時系列断面画像を生成する様子を説明する図である。
【図１０】第２実施形態による時系列断面画像の表示形態を説明する図である。
【図１１】第２実施形態による時系列断面画像の表示例を示す図である。
【図１２】第３実施形態の呼吸動態撮影装置における時系列断面画像の生成処理を示すフローチャートである。
【図１３】胸部動態画像の各フレームを模式的に表した図である。
【図１４】各フレームから注目する画素位置を指定する様子を示した図である。
【図１５】指定された画素位置から、時空間上の任意の曲面を形成し、２次元画像化することを説明する図である。
【図１６】第４実施形態の呼吸動態撮影装置における時系列断面画像の生成処理を示すフローチャートである。
【図１７】胸部動態画像の各フレームを模式的に表した図および部位の移動範囲を示す図である。
【図１８】時空間３次元像で図１１の移動範囲を表した図およびその範囲の画素値を統合し画像を形成した図である。

Claims (21)

1. 対象物の動態画像を構成する複数の画像を時間軸に沿って並べた時空間３次元空間に設定された曲面と前記複数の画像の各々との交線に沿って並ぶ画素列を取得する取得手段と、
   前記取得手段により前記複数の画像から取得された複数の画素列に基づいて画像を生成する生成手段とを備えることを特徴とする画像処理装置。
2. 前記生成手段は、前記複数の画素列を、取得元の画像の時系列順に基づいて並べることを特徴とする請求項１に記載の画像処理装置。
3. 前記曲面を操作者による入力に基づいて設定する設定手段を備えることを特徴とする請求項１に記載の画像処理装置。
4. 前記生成手段は、前記曲面上の前記複数の画素列を平面に投影して得られる２次元画像を生成することを特徴とする請求項１に記載の画像処理装置。
5. 前記生成手段は、前記曲面上の前記複数の画素列を平面に展開して得られる２次元画像を生成することを特徴とする請求項１に記載の画像処理装置。
6. 前記生成手段は、前記時空間３次元空間における前記曲面上の前記複数の画素列の２次元表現としての２次元画像を生成することを特徴とする請求項１に記載の画像処理装置。
7. 前記生成手段により生成された画像を表示する表示手段を備えることを特徴とする項１乃至６のいずれかに記載の画像処理装置。
8. 前記動態画像は、前記対象物を透過した放射線の強度分布に基づいて生成されたものであることを特徴とする請求項１乃至７のいずれかに記載の画像処理装置。
9. 請求項１乃至８のいずれかに記載の画像処理装置と、対象物の動態を撮影する撮影装置とを備えることを特徴とするシステム。
10. 対象物の動態画像を構成する複数の画像の各々から、設定された矩形領域の一辺の方向に並ぶ複数画素を統合することにより画素列を生成する統合手段と、
    前記統合手段により前記複数の画像から生成された複数の画素列に基づいて画像を生成する生成手段とを備えることを特徴とする画像処理装置。
11. 前記生成手段は、前記複数の画素列を、取得元の画像の時系列順に基づいて並べることを特徴とする請求項１０に記載の画像処理装置。
12. 前記統合手段は、前記複数の画像を時間軸に沿って並べた時空間３次元空間に設定された平行六面体と前記複数の画像の各々との交面を前記矩形領域として前記画素列を生成することを特徴とする請求項１０に記載の画像処理装置。
13. 前記平行六面体を操作者による入力に基づいて設定する設定手段を備えることを特徴とする請求項１２に記載の画像処理装置。
14. 前記設定手段は、前記時空間３次元空間における前記動態画像及び前記平行六面体の２次元的表現を与える表示手段を含むことを特徴とする請求項１３に記載の画像処理装置。
15. 前記統合手段は、前記複数画素の画素値の統計量を画素値とする１つの画素で前記複数画素を置換することを特徴とする請求項１０に記載の画像処理装置。
16. 前記動態画像は、前記対象物を透過した放射線の強度分布に基づいて生成されたものであることを特徴とする請求項１０乃至１５のいずれかに記載の画像処理装置。
17. 請求項１０乃至１６のいずれかに記載の画像処理装置と、対象物の動態を撮影する撮影装置とを備えることを特徴とするシステム。
18. 対象物の動態画像を構成する複数の画像を時間軸に沿って並べた時空間３次元空間に設定された曲面と前記複数の画像の各々との交線に沿って並ぶ画素列を取得する取得工程と、
    前記取得工程において前記複数の画像から取得された複数の画素列に基づいて画像を生成する生成工程とを備えることを特徴とする画像処理方法。
19. 対象物の動態画像を構成する複数の画像の各々から、設定された矩形領域の一辺の方向に並ぶ複数画素を統合することにより画素列を生成する統合工程と、
    前記統合手段において前記複数の画像から生成された複数の画素列に基づいて画像を生成する生成工程とを備えることを特徴とする画像処理方法。
20. 請求項１８又は１９に記載の方法をコンピュータに実行させることを特徴とするプログラム。
21. 請求項２０に記載のプログラムを格納したことを特徴とするコンピュータ可読記憶媒体。

Images