WO2005039416A1

WO2005039416A1 - 治療支援用画像処理装置

Info

Publication number: WO2005039416A1
Application number: PCT/JP2004/015692
Authority: WO
Inventors: Osamu Miyazaki; Hiromitsu Hayashi
Original assignee: Hitachi Medical Corporation
Priority date: 2003-10-23
Filing date: 2004-10-22
Publication date: 2005-05-06
Also published as: US20080027314A1; JPWO2005039416A1; US7817828B2

Abstract

　治療の目標となる治療範囲を予め設定しておいて、治療中に既に治療が完了したと考えられる既治療領域の画像を取得し、この既治療領域に基づいて目的治療領域全体の治療が終了するまでの時間を推定する。また、治療中の部位の変化を色で表現させて表示する。また、複数の治療時間について治療の完了したイメージをシミュレートすることでよりわかりやすく治療が終了するまでの時間を決定する。以上の特徴により、治療精度を向上させ、治療に要する総時間や被曝の低減を図ることができる。

Description

明細書

治療支援用画像処理装置

技術分野

[0001] 本発明は、被検体の対象部位を治療するに際し、その対象部位の画像をモニタリングして、治療時の計画性を向上させることのできる治療支援用画像処理装置に関する。

背景技術

[0002] X線 CT装置や MRI装置などは周知のように有用な画像診断装置であるが、最近では単なる画像を表示するための診断装置としての利用以外に外科手術の治療モニタリングや生検のガイドなどに広く応用されるようになってきている。生検の場合は、穿刺針と病巣の位置関係の把握が目的であり、再構成処理を工夫し逐次画像を更新することでリアルタイムに位置関係が把握できるような CTZMR透視等が用いられている。治療の場合は、レーザーや RF(Radio Frequency)を用いた加温療法やガスを用いた冷凍治療など様々な手技の有効性が検討されている。また、モニタリングに用いるモダリティも CT装置、 MRI装置、 US装置など各種画像診断機器を用いて試みられており、それぞれ長所短所がある。このような治療モニタリングに関するものとして、特許文献 1一 3に示すようなものがある。

特許文献 1：特開 2003 - 10228号公報

特許文献 2：特開 2000 - 300591号公報

特許文献 3：特開 2003 - 88508号公報

発明の開示

発明が解決しょうとする課題

[0003] 例えば、レーザーによる椎間板ヘルニアの治療では椎間板にレーザーを照射することで椎間板内の水分を蒸散させ、神経の圧迫を弱める。 CT装置ではその蒸散ガス発生の程度などが確認できれば蒸散の様子を推定することができる。また、例えば、腫瘍などの治療では病巣を加熱あるいは冷凍してその組織を壊死させることが行われている。その一つである冷凍治療 (Cryosurgery)は、凍った領域 (アイスボール (Ice Ball))が病巣を十分に覆っているかを確認しながら治療を進めるという手法を採っている。

[0004] 超音波診断装置 (US)ではこのアイスボールの陰が観察できないため、組織コントラストの高い MRI装置がモニタリングに用いられることが多い。また、 CT装置をモニタリングに用いることも可能である。モニタリングの際の一つの問題点は、治療の合間に逐次画像を取得しなければならないという煩雑さであり、さらに X線 CT装置の場合は画像取得回数に応じて患者の被曝量が増加するという問題を有する。また、治療領域は、病巣よりもある程度大きくしなければならないが、正常な組織はできるだけ残存させた、と、う要求があるため、現時点では治療領域を目視などで決定するしかなかつた。しかし、目視の場合、治療が進むにつれ凍結領域と病巣境界との正確な区別がっきにくくなり、その把握が困難になってくるなどの問題があった。さらに、病理学的に組織温度がどれぐらいになれば壊死するかと、う基礎データはある力定量的にそれらを測定し、し力も病巣と温度分布の関係を確認する手段は無力たため、治療手技を確実に進める上での問題点となっていた。

[0005] この発明は、上述の点に鑑みなされたものであり、治療のモニタリング時に治療が終了するまでに要する時間を推定することのできる治療支援用画像処理装置を提供することにある。

[0006] また、この発明は、上述の点に鑑みなされたものであり、治療のモニタリング時に病巣と治療効果を期待できる温度の領域を容易に把握することのできる治療支援用画像処理装置を提供することにある。

課題を解決するための手段

[0007] 本発明によれば、被検体の断層像を取得し治療経過をモニタリングする治療支援用画像処理装置において、前記断層像上力モニタリング画像を選択し、その上に目的とする治療領域を設定する目的治療領域設定手段と、前記モニタリング画像上に既に治療したと考えられる既治療領域を設定する既治療領域設定手段と、前記目的治療領域設定手段により設定した目的治療領域の治療が終了する時間を前記既治療領域に基づ!、て推定する総治療時間推定手段と、を備えた治療支援用画像処理装置を提供する。 [0008] また、前記総治療時間推定手段は、複数のモニタリング画像内の既治療領域について、治療の開始地点である基準点力前記既治療領域の境界までの距離をそれぞれ求め、前記距離の変化に基づいて前記総治療時間を推定することを特徴とする

[0009] また、前記総治療時間推定手段は、前記目的治療領域上、前記基準点から等角度間隔で放射状に伸びる直線を基準線とし、前記基準線上における治療進度をもとに前記総治療時間を推定することを特徴とする。

[0010] また、前記総治療時間推定手段は、前記基準点から前記目的治療領域を二等分する少なくとも一本の基準線を指定し、前記基準線を所定間隔に分割する垂線が前記目的治療領域と交差する点を求め、これら点と前記基準点を結ぶ線をさらに基準線とし、前記基準線上における治療進度をもとに前記総治療時間を推定することを特徴とする。

[0011] また、前記総治療時間推定手段は、前記目的治療領域の輪郭を表示する前記表示手段上の画素中の所定間隔の画素と前記基準点を結ぶ線を基準線とし、前記基準線上における治療進度をもとに前記総治療時間を推定することを特徴とする。

[0012] また、前記目的治療領域設定手段は、複数の閉曲線を合成してから前記閉曲線の交差部分をスムージングして設定し、前記総治療時間推定手段は、前記閉曲線毎に基準線を設定し、前記複数の基準線上における治療進度をもとに前記総治療時間を推定することを特徴とする。

[0013] 本発明によれば、被検体の断層像を取得し治療経過をモニタリングする治療支援用画像処理装置において、前記断層像上力モニタリング画像を選択し、その上に目的とする治療領域を設定する目的治療領域設定手段と、前記モニタリング画像に既に治療したと考えられる既治療領域を設定する既治療領域設定手段と、上記既治療領域に基づいて、所定時間経過後の治療領域を予測し、この予想される治療領域を少なくとも 1つの色情報で表示する表示手段と、を備えた治療支援用画像処理装置を提供する。

[0014] また、前記表示手段は、前記目的治療領域と前記既治療領域と前記予測治療領域とにはそれぞれ異なる色情報を割り当てて表示する。 [0015] また、前記表示手段は、上記目的治療範囲に色情報を割り当てて表示することを特徴とする。

[0016] また、前記表示手段は、上記目的治療領域に色情報を割り当てて表示するとともに

、治療中の部位の変化を色で表示することを特徴とする。

[0017] また、前記表示手段は、複数時点における上記既治療領域にそれぞれ異なる色情報を割り当てて表示することを特徴とする。

[0018] また、前記表示手段は、上記予測された治療領域に色情報を段階的に割り当てて表示することを特徴とする。

[0019] また、前記表示手段は、複数時点における上記予測された治療領域をそれぞれ所望の色階調で表示することを特徴とする。

[0020] また、前記表示手段は、前記予測された治療領域と、前記既治療領域とのうち少なくとも一つを、前記目的治療領域に重ね合わせて表示することを特徴とする。

[0021] 本発明によれば、被検体の断層像を取得し治療経過をモニタリングする治療支援用画像処理装置におヽて、前記断層像を作成するためのデータを記憶する画像メモリと、前記断層像作成用のデータから計測値を求める計測値算出手段と、前記計測値を温度に読み替えする計測値 -温度変換手段と、前記読み替えされた温度に対して色情報を割り当て温度色表示をするカラーテーブルと、前記画像メモリに記憶された断層像作成用のデータにグレースケール階調を割り当てるルックアップテープルと、前記温度色表示とグレースケール階調が割り振られた断層像作成用のデータとの合成を指示するモニタリングモードスィッチと、前記合成された温度色表示と断層像作成用のデータを画像として表示する表示手段と、を備えた治療支援用画像処理装置を提供する。

[0022] また、前記計測値温度変換手段は、 CT値を計測値とするものであり、 25度付近での CT値が 0である組織が 0度以下にお!、て— 20力 50の CT値を有するとして、 CT値を温度に読み替えることを特徴とする。

[0023] また、治療支援のための表示モードを選択する治療支援表示モード選択手段と、表示したい治療進行イメージの治療開始後の時間を入力する時間選択入力手段と、前記選択された時間における治療進行イメージを作成する治療進行イメージ作成手段と、前記目的治療領域と前記既治療領域のうち少なくとも一方と前記治療進行ィメ一ジとを表示する表示手段と、を備えたことを特徴とする。

[0024] また、前記目的治療領域と前記既治療領域のうち少なくとも一方と前記治療進行ィメージとの間の差分画像を作成する差分画像作成手段と、を備え、前記表示手段は、前記目的治療領域と前記既治療領域のうち少なくとも一方と前記差分画像とを表示することを特徴とする。

[0025] また、前記表示手段は、前記目的治療領域と、前記既治療領域と、前記治療進行イメージと、前記差分画像とに、それぞれ異なる色情報を割り当てて表示することを特徴とする。

発明の効果

[0026] この発明によれば、治療のモニタリング時に治療が終了するまでに要する時間を推定することで、残りの治療計画の精度向上、さらには画像取得回数を低減することができる。また、治療のモニタリング時に病巣と治療効果を期待できる温度の領域を容易に把握することができると、う効果がある。

図面の簡単な説明

[0027] [図 1]本発明による治療支援用画像処理装置の実施例を示すブロック図である。

[図 2]本発明の治療支援用画像処理装置の実施例に係る画像処理装置の構成を示す図である。

[図 3]肝臓の冷凍治療時の治療支援表示の一例である。

[図 4]治療支援用画像処理装置の動作を示すフローチャート図である。

[図 5]モニタリング画像における治療時間と、治療領域を示す距離との関係を示すグラフである。

[図 6]基準線ゃ閉領域の設定方法の別の例を示す図である。

[図 7]複数のプローブを用いた場合の閉領域の設定方法を示す図である。

[図 8]本発明の治療支援用画像処理装置の別の実施例に係る画像処理装置の構成を示す図である。

[図 9]図 8の計測値 -温度変換手段の一例を示す図である。

[図 10]実施例 3に関わる治療支援用画像処理装置の動作を示すフローチャート図である。

[図 11]図 10の治療支援用画像処理装置の画像処理部の構成を示す図である。発明を実施するための最良の形態

[0028] [実施例 1]

以下、本発明の実施例 1を添付図面に基づいて詳細に説明する。

図 1は、本発明による治療支援用画像処理装置の実施例を示すブロック図である。この治療支援用画像処理装置は、システム全体を統括制御するホストコンピュータ 1 と、 X線源 2、 X線制御部 3、 X線検出器 4、計測回路 5などの計測部を搭載した回転盤 14と、この回転盤 14の回転走査を制御するスキャナ制御部 6と、被検体の位置決め時やらせん走査時の搬送用被検体テーブル 7及びテーブル制御部 8と、前処理，再構成処理をはじめとする各種画像処理を実施する画像処理装置 9などからなる。計測制御部 10は、ホストコンピュータ 1からの指示に従って、スキャナの回転盤 14に搭載された X線制御部 3及び計測回路 5の動作を制御する。 X線制御部 3及び計測回路 5は、計測制御部 10力の指示に従って X線の曝射及びデータ計測の計測動作を開始する。外部記憶装置 11は、装置全体の制御プログラムを格納したり、計測回路 5から出力される計測データやそれを処理して得られた断層像データや各種プログラム等を格納した磁気ディスク、フロッピディスクドライブ、ハードディスクドライブ、 CD— ROM ドライブ、光磁気ディスク (MO)ドライブ、 ZIPドライブ、 PDドライブ、 DVDドライブなどである。表示装置 12は、画像データを一時記憶する表示メモリと、この表示メモリからの画像データに基づ、て画像を表示する CRTディスプレイなどである。入力デバイス 13 は、画面上のソフトスィッチを操作するマウスやそのコントローラ、各種パラメータ設定用のキーやスィッチを備えたキーボードなどで構成され、各種の指令や情報などをホストコンピュータ 1に入力するものである。ホストコンピュータ 1は、このような構成の治療支援用画像処理装置を通信インターフェイスを介して LAN (ローカルエリアネットヮーク)やインターネット、電話回線などの種々の通信ネットワーク上に接続可能とし、他のコンピュータや CT装置などとの間で画像データのやりとりを行えるようにしてもよい

[0029] 図 2は、この発明の治療支援用画像処理装置の一実施の形態に係る画像処理装置の構成を示す図である。この画像処理装置は、目的治療領域設定部 21、基準点設定部 22、既治療領域設定部 23、総治療時間推定部 24及び治療時間提示部 25からなる。目的治療領域設定部 21は、治療時のモニタリング断面で病巣を十分に覆うことのできる領域を指定するものであり、領域指定は円形、矩形、自由曲線などの様々な形状で指定することができるようになって、る。冷凍治療の場合のように楕円形の治療領域 (アイスボール)が成長する場合などその治療方法に応じて使い分けるようになっている。基準点設定部 22は、一般的には治療の原点を指定するものであり、冷凍治療では最も温度の低、ところなどで、例えば治療用プローブの先端部などを基準点として設定する。既治療領域設定部 23は、既に組織が壊死するであろう温度に達したと考えられる領域を設定するものであり、画像上で操作者自身が指定したり、閾値処理などで自動的に抽出するようにしても良い。総治療時間推定部 24は、目的治療領域設定部 23で設定された領域を治療するのに必要な所要治療時間を算出するものである。治療時間提示部 25は、総治療時間推定部 24で算出された推定終了時間を表示装置 12などに表示して観察者に提示するものである。

[0030] 図 3は、肝臓の冷凍治療時の表示部での表示イメージである。図 3には、点線で示した目的治療領域 31、時刻 tlにおける既治療領域 Bl、時刻 t2における既治療領域 B2、治療の原点を示す基準点 P、及び治療終了時間を推定するための基準線 L1, L2のそれぞれの関係が示してある。すなわち、プローブ 30の先端から放出されたガスが放射口の近傍である基準点 Pから順次組織を凍らせ、次第に病巣部を覆って行くようすが図 3には示されている。図 4は、この治療支援用画像処理装置の動作を示すフローチャート図である。以下、この治療支援用画像処理装置の動作について説明する。

[0031] ステップ S40では、この実施の形態に係る治療支援用画像処理装置は、治療に先立って穿刺する断面、すなわちモニタリングする断面を決定し、それをモニタリング画像 IMG0として取得する。

[0032] ステップ S41では、目的治療領域設定部 21が GUI(Graphical User Interface)を用いてモニタリング画像 IMG0に目的治療領域 31を設定する。

[0033] ステップ S42では、治療が開始されるので、その治療経過時の画像を観察するために適宜治療経過画像が取得される。治療が開始してから時刻 tl経過時の治療経過画像としてモニタリング画像 IMG1と、時刻 t2経過時の治療経過画像としてモニタリング画像 IMG2とがそれぞれ取得される。

[0034] ステップ S43では、治療経過画像の中力治療による凍結領域力例えば MRIでは黒い低信号として描出されるので、既治療領域設定部 23は、モニタリング画像 IMG1

, IMG2から凍結領域を抽出し、その境界部座標を各モニタリング画像 IMG1, IMG2における既治療領域 Bl, B2として設定する。

[0035] ステップ S44では、基準点設定部 22が、モニタリング画像 IMG1に基づ、て、基準点

Pを定義する。

[0036] ステップ S45—ステップ S49までの処理は、総治療時間推定部 24が行う処理である。

[0037] ステップ S45では、治療領域の拡大方向に基準点 Pから時間推定に用いる直線 Ll( 治療領域を算出するための基準線)を定義する。最初に定義する基準線 L1は図 3に示すようにガスの放射方向に沿った基準線 L1である。基準線は 1本でも良いが、複数定義することで時間推定精度を向上することができる。従って、基準線を複数定義する場合は基準線 L1となす角度 Θをパラメータとするか、あるいは境界を定義した個々の座標を用いて求めるようにしてもよい。図 3では、基準線 L1にたいして角度 θ 1の基準線 L2を定義している。

[0038] ステップ S46では、基準点 Pから定義した基準線 L1と、各々のモニタリング画像 IMG1

, IMG2の既治療領域 Bl, B2との交点 CI, C2を求める。

[0039] ステップ S47では、基準点 Pから各交点 CI, C2までの距離 dl, d2を求める。

[0040] ステップ S48では、対象画像であるモニタリング画像 IMG1及び IMG2の少なくとも 2枚の画像に対してステップ S46, 47の処理が終了したか否かの判定を行い、終了した場合にはステップ S49に進み、終了していない場合にはステップ S46にリターンし、同様の処理を繰り返し実行する。

[0041] ステップ S49では、治療経過を示すモニタリング画像 IMG1, IMG2の凍結領域である既治療領域 Bl, B2に達するまでに、それぞれの治療時間 tl, t2を要したことは既知なので、あら力じめ設定した目的治療領域 31と基準線 L1の交点 Cまでの距離 Dを求め、基準点 Pから各交点 CI, C2までの距離 dl, d2に到達するのに要した治療時間 tl , t2との関係から、治療に要する時間 Tを算出する。図 5は、上述のモニタリング画像 IMG1, IMG2における治療時間 tl, t2と、治療領域を示す距離 dl, d2との関係を示すグラフである。図 5からも明らかなように、これらの比力も (あるいは回帰によって)目的治療領域 31の治療に要する時間 Tを算出することができる。なお、複数のプローブを用いた場合には、その基準点をプローブの数だけ設定し、各々について基準線を設定して求めれば良い。例えば、基準線を複数設定した場合には、それぞれの基準線において前述の時間 Tを算出し、最も大きい時間 Tを治療に要する時間としたり、各基準線において算出された時間 Tの平均を治療に要する時間とすることができる。

[0042] ステップ S4Aでは、治療時間提示部 25が推定すべく算出された総治療時間を数値で観察者に表示する。なお、この治療時間提示部 25は、既治療領域と 1分後、 2分後の治療領域を予測し、その予測治療領域を断層像上に段階的に色情報で示すようにしてもよい。また、表示する色やグレーには透過度を設定できるようにし、モニタリング画像の濃淡情報と重畳表示可能にしてもよい。ここで色情報とは、白、黒及びダレ一を含めたカラーの色合いやその濃度および模様を含むものとする。また、モニタリング画像では無ぐ治療前の画像に重畳表示する方が病巣との位置関係の把握が容易な場合は、治療前の画像を用いても良い。

[0043] 上述の実施の形態では、 2枚のモニタリング画像 IMG1, IMG2に基づいて総治療時間を算出する場合について説明したが、複数 n枚のモニタリング画像に基づいて判断することによって、より高精度に総治療時間を推定することが可能となる。

[0044] なお、基準線ゃ閉領域の設定方法は種々考えられるが、図 6(A)に示すように閉領域の重心位置 GOから等角度間隔 (図では約 15° 間隔)で放射状に交点を求める。交点と基準点を結ぶ線を基準線とし、基準線上における治療進度をもとに前記総治療時間を推定することになる。

[0045] あるいは、図 6(B)のように基準線 L1に垂直な直線を等間隔 (不当間隔でも当然良いが)で設定して交点を求めることも考えられる。基準点から前記目的治療領域を二等分する少なくとも一本の基準線 L1を指定し、基準線を所定間隔 (等間隔でも不等間隔でもよい)に分割する垂線が目的治療領域と交差する点を求め、これら点と基準点を結ぶ線をさらに基準線とし、基準線上における治療進度をもとに前記総治療時間を推定することになる。

[0046] あるいは、目的治療領域の輪郭を表示する表示部上の画素中の所定間隔の画素と基準点を結ぶ線を基準線とし、基準線上における治療進度をもとに総治療時間を推定してもよい。これはすなわち、閉曲線の境界画素の画素数を等しくサンプルしてこれらサンプル点に向力つて追カ卩の基準線を引いて推定に用いるものである。

このようにして、求めた交点と基準点 Pを結ぶことによって基準線を複数形成することができるので、前述のように所定時間経過後の治療領域を予測し、その予測治療領域を高精度に描画することができるようになる。

[0047] 図 7は、複数のプローブを用いた場合の閉領域の設定方法を示す概念図である。

図 7に示すように、複数のプローブ 30a, 30bを用いた場合には、複数の基準線 La, Lb と、両方の閉曲線 70a, 70bを含むような閉領域を設定することが好ましいが、両方の閉領域 70a, 70bの論理和をとつた領域が実際の治療領域であり、なおかつ図 7に示すように両方の閉領域 70a, 70bの論理和をとつた閉領域の場合はその境界付近がくぼんだ領域となる力実際にはこのくぼんだ領域にはさらに広く治療領域が及んでいると考えられるので、図示のようにくぼんだ領域はスムースな閉領域として点線のようなスムージンク処理を施した合成領域 71(点線で示す領域)を定義するのが望ましい。

[0048] スムージングを施した合成領域 71は、境界線の連続性、角度、曲率などの，情報を用いて決定すればよい。具体的なスムージングの例としては、両閉領域の法線のなす角度など定量値を基にスムージングの度合いを決定しても良いが、重心力もの距離に対して移動平均を取り、距離が大きくなる領域のみスムージング結果を採用するのが最も簡易である。総治療時間の予測に関しては複数の基準線上での治療進度の平均値、最大値、及び最小値等任意のものを使い予測可能である。

[0049] また、既治療領域を複数時点で取得して関数近似することで総治療時間を推測してもよい。たとえば、図 3の既治療領域 Bl, B2の外形を座標上の時間の関数として近似式を作成し、時間変数に所定の時間を代入して治療時間と治療範囲の関係を求めてもょ、。ここで近似する手法としてはカーブフィッティングなど公知のものが使用可能である。

[0050] [実施例 2] 本実施例は基本的に実施例 1と同様の構成や処理を用いており、特に断らない限り実施例 1におけると同じ処理をする。また、同じ参照番号は実施例 1におけると同一部分を示す。

[0051] 図 8は、この発明の治療支援用画像処理装置の別の実施例に係る画像処理装置の構成を示す図である。この画像処理部 9は、図 8に示すように、計測値算出部 82、計測値温度変換部 83、カラーテーブル 84、入力スィッチ 85、ルックアップテーブル (LUT)86、画像合成部 87、及び画像メモリ 81から構成される。画像メモリ 81は、計測回路 5から出力される計測データ、計測値、又はこれを再構成して得られた断層像を一時的に記憶するものである。ルックアップテーブル (LUT)86は、 CT画像を表示するためのグレースケールを与えるルックアップテーブルを常に保存しており、例えば、表示階調の中心が CT値で 40、幅が 100であれば、—10— 90までの CT値を 0— 255のグレースケール階調に割り振るようになつている。例えば、 CT値 40の画素は 127レベルの階調値になり、 RGB値はすべてが等しくなる。ルックアップテーブル (LUT)86は、計測データ又は断層像から計測値を算出後、上記のようにグレースケール階調に割り振るようになつていてもよい。

[0052] 計測値算出手段 82は、画像メモリ 81からの計測データに基づいて CT値を算出する。計測値温度変換手段 83は、たとえば図 9に示すような CT値温度変換テーブルを有しており、計測値算出手段 82力出力される CT値を対象部位の温度に変換する。この CT値-温度変換テーブルは、あら力じめ動物実験を実施し、治療中の温度を測定しながら治療シミュレーションを行って求めておくのが最も信頼性が高いと考えられる。別の手法としては、水の温度を変化させながらそれを撮影して代用することも可能である。いずれにしても、図 9に示すような横軸に温度、縦軸に CT値とした対応曲線を求め、それをテーブルィ匕したものである。カラーテーブル 84は、温度に対して色調、濃度、透明度等の色情報を割り当てるもので、例えば温度が上昇すると赤に近づき、下降すると青に近づくなどの設定がなされているカラーテーブルを備え、変換後のカラー情報をモニタリングモードスィッチ 85を介して画像合成手段 87に出力する

[0053] X線 CT装置の場合、再構成画像は CT値と呼ばれる絶対値をもとに作成されるので、再現性良く定量ィ匕が可能である。実際には、冷凍治療中の凍結状態から常温までの間に 20— 50程度の CT値の変化があると言われている。モニタリングモードスィッチ 85は、モニタリングモードを切り替えするものである。図 8は、画像合成部 87は、モニタリングモードスィッチ 85がオン/オフで制御されるものを例示しており、モニタリングモ一ドスイッチ 85がオフの場合は、通常通り、表示部上や画像メモリ上の画像を LUT86 によって指定された階調で表示し、オンの場合は、適宜、温度情報等を元にして決定された色情報と合成して表示する。

[0054] 図 8の画像処理装置の動作例を説明する。まず、治療中にモニタリングモードを選択すると、モニタリングモードスィッチ 85がオンとなる。そして、画像メモリ 81に記憶されている画像の関心領域が入力デバイス 13によって設定されると、計測値算出部 82 はその領域の各画素の CT値を読み出し、 CT値温度変換部 83に出力する。 CT値温度変換部 83は、入力される CT値を温度情報に変換してカラーテーブル 84に出力する。ノイズが多い場合は断面内、あるいは複数のスライスが取得されていればスライス方向にフィルタ処理を施すことで、近傍画素との平均的な温度として算出する場合もある。

[0055] 冷凍治療 (Cryosurgery)時には、治療前を青 (ある!/、は透明)、治療効果が期待できる温度で赤に変化するように、ある、は連続的に近づ、て、くように表示することが好ましい。このような色の設定においては種々の変更が可能で、これに限定されるものではないが、治療効果の期待できる温度に達したことが容易に視認できることが望ましい。

[0056] 画像合成部 87は、温度を元にカラーテーブル 84で変換されたカラー情報と、元の CT画像の表示色とを合成して表示装置 12に表示する。例えば、透過度パラメータ q を用い、各 RGB値の加重和を新たな色として表示する。例えば、元の階調を Cgray、温度対応色が Cred, Cblue, Cgreenとすれば、新しい表示色 R、 G、 Bは、

R=(l. 0-q) X Cred + q X Cgray

G = (i . 0— q) X Cgreen + q Xし gray

B = (l. 0-q) X Cb lue + q X Cgray

となる。 [0057] 本発明はこのような色の合成方法に限定されるものではない。また、表示部にはァキシャルだけでなぐサジタル、コロナルなどの断面を同時に表示しておくことが望ましい。特に、 X線ビーム幅の広いコーンビーム CTにおいては一回の撮影でアイスボールの成長の様子が把握できるため特に好適である。

[0058] なお、時系列画像を比較し、アイスボールの成長速度を評価する手段を設け、病巣部を十分アイスボールが覆うまでに要する時間を推定するような付加機能を追加してちょい。

[0059] また、そのときの成長速度の評価手段は、ガスの照射中心から治療領域の端までの距離と時間との関係で定義することができる。表示手段の好適な例としては、現状の治療領域と 1分後、 2分後の治療領域で段階的に色又は透過度を変更するように表示することが好ましい。表示部には、設定された色と治療前や治療中の画像の色とを合成したものであって、画像は断面変換した画像 (MPR)を含んでもょ、。

[0060] [実施例 3]

上記実施例 1は治療完了までの時間を予測して表示するものであった。時間表示に換えて、あるいは時間表示に加えて、治療の進行していくイメージを作成して表示してもよい。こうすることで、たとえばアイスボールが十分病巣部を覆うまでに要する時間や条件を順次に一目で確認できる。

[0061] 図 10と 11を用いて実施例 3を説明する。なお、本実施例は基本的に実施例 1と同様の構成や処理を用いており、特に断らない限り実施例 1におけると同じ処理をする。また、実施例 1と同じ参照番号は同一部分を示す。

[0062] 上記の治療進行イメージの表示は、たとえば実施例 1に記載した総治療時間推定部 94の機能を改変して、たとえば 1分後や 2分後など分刻みでは治療はどこまで進行するかを推測し、これら複数時点で治療が完了したイメージを元の断層像に重ね合わせて表示する。

[0063] 図 11に示す通り、本実施例の画像処理部 9には、画像再構成処理部 9Z、画像メモリ 81、実施例 1において説明したと同様の画像解析部 96に加えて、支援画像作成部 9N が含まれる。支援画像作成部 9Nには、治療支援モード選択部 9G、表示治療進行ィメージの時間選択入力部 9H、治療進行イメージ作成部 91、患部モニタリング画像との差分画像作成部 9J、および画像作成部 9Kが含まれる。

[0064] 図 10のステップ 40力ステップ 49までにお!/、て、たとえば図 3に示すような目的治療領域 31の設定に基づいて、既治療領域 Bl, Β2の基準点 Ρから各交点 CI, C2までの距離 dl, d2と、既治療領域 Bl, B2が凍結するまでの治療時間 tl,t2との関係から、治療の進度を把握する。ここまでのステップは実施例 1と同様であり、実施例 1に記載の通り基準線 L2も併用したり、図 6のようにさらに多くの基準線を使用したり、図 7のように複数のプローブを使用してもよい。これらの場合でも、それぞれの基準線中、最も速く患部を治療完了する基準線上の進度や、最も遅く患部を治療完了する基準線上の進度や、平均化した進度を把握可能であり、それらをもとにして計算した時間を総治療時間として使用する。

[0065] ステップ 4Bで図 11の治療支援表示モード選択部 9Gは治療支援表示のモードを選択する。選択可能なモードとしては、実施例 1と同じ総治療時間 τ、後述する任意の指定時間における治療進行イメージ、任意の指定時間における治療進行イメージと目的治療領域との差分画像、任意の指定時間における治療進行イメージと既治療領域との差分画像のいずれか、後述するグラフ、治療誤差画像、あるいはこれらのうち任意のものとの組み合わせなどがある。

[0066] 指定時間における治療進行イメージと指定時間における治療進行イメージとの差分画像のうち!ヽずれか一方を選択した場合は、表示すべき治療進行イメージの時間を図 11での表示治療進行イメージの時間選択入力部 9Ηに入力して指定する。この際、ステップ 4Αにて総治療時間 Τを事前に提示しておいて、その時間を参照して治療進行イメージを得たヽ時間を指定してもよヽし、何も指定しな、場合はデフォルトで総治療時間 Τが指定されるようにしておいてもよい。また、この時間は複数指定できるようにしてちょい。

[0067] ステップ 4Dでは、図 11の治療進行イメージ作成部 91にて、上記ステップ 4Cで入力された時間における治療進行イメージを作成する。

[0068] ステップ 4Εでは、図 11の患部モニタリング画像との差分画像作成部 9Jにて、上記ステツプ 4Dで作成された治療進行イメージと治療前の画像 IMG0との差分画像を作成する。なお、ここで差分する画像は、治療前の画像 IMG0に限らず、さらに一つ前の治療進行イメージであってよ、。

[0069] ステップ 4Fでは、ステップ 4Bで選択された治療支援表示のモードに応じて、総治療時間 T、指定時間における治療進行イメージ、指定時間における治療進行イメージとの差分画像の、ずれか、あるいはこれらのうち任意のものの組み合わせの、ずれかの画像表示を合成する。上記合成は図 11の画像作成部 9Κにて行う。なお、画像作成に際しては、異なった時間の治療イメージを別の色で表示したり、ある時間における治療部位周辺の温度分布も実施例 4に記載した手法にてさらに重ね合わせて表示してもよい。こうすることで、たとえば、現在と 1分後と 2分後の治療領域をそれぞれ別の色や透過度で表示することができる。また、現在と 1分後と 2分後の治療領域をそれぞれ別の色で表示するとともに、その温度分布をそれぞれの色の透過度で表示してちょい。

[0070] そのほか、色座標や領域を描画する線の幅や線種を駆使して、時間に応じた治療進行イメージや温度分布を様々に表示してもよい。

[0071] 以上のように表示することで、治療部位全体とその周辺の正常組織にとって最適な治療時間を決定するのに有用となる。

[0072] ステップ 4Εにおける治療進行イメージと治療前の画像 IMG0との差分画像を作成すると、治療中の進拔状況を把握しやすくすることができる。差分画像は微分成分であるので、治療進度の加速度が評価しやすい。たとえば、差分画像ではプローブ 30やプローブにより生じうるアーチファクト成分も都合よく除去される。また、データを差分した結果、治療が進行した領域にのみ画像が現れるので好都合である。

[0073] 差分画像作成時に差分に用いた、たとえば IMG0の画像上に重畳表示することもできる。あるいは、差分の面積や体積 (3次元像の場合)を求めて、その面積変化や体積変化を時間軸に対してグラフ表示することで治療進度の変化が把握しやすくなる。

[0074] 差分画像の作成にあたっては、治療中にプローブがずれることがあるので、閾値処理ゃノイズ除去処理が必要となる。ノイズ除去処理はラベリング後に面積を求め、特定の面積以上のもののみを対象とする公知の手法を使用できる。

[0075] なお、本実施例で示した差分画像作成部 9Jは、 X線 CT装置によって治療経過をモユタする場合、治療経過を予測表示するために用いるだけでなぐ更新して取得した既治療領域と治療進行イメージとの治療誤差画像を作成し、誤差から血管等の存在によって温度伝達が遅れている領域の抽出に用いることが可能である。なお、目的治療領域、既治療領域、治療進行イメージ、及び差分画像には、それぞれ異なる色情報を割り当ててもよい。また、割り当てる色は、上述した計測値温度変換手段によつて求めた温度に応じて変化させてもよい。

[0076] この場合、曲率や不連続点などを画像上で認識して、差異の領域を抽出していけばよ、が、これに限らず差異を認識して抽出できれば、かなる方法であってもよ、。こうすることで、治療の遅れた部分に血管の存在を疑う機会が治療をする者に与えられる。治療者はこれを判断し治療計画を変更可能であり、治療対象ではない血管が閉塞してしまうことなどを避けることができる。

[0077] なお、以上、ずれの実施例にお!、ても、使用する画像は断面変換した画像 (MPR) でもよい。また、表示部にはアキシャルだけでなぐサジタル、コロナルなどの断面や三次元像を同時に表示しておくことが望ましい。

[0078] 特に、 X線ビーム幅の広いコーンビーム CTを使用した場合には、一回の撮影で三次元像およびアキシャル、サジタル、およびコロナルなどの視点からのアイスボールの成長の様子を把握できるため治療が一層容易となる。

産業上の利用可能性

[0079] 以上、本発明によれば、治療時にモニタリングをしながら、治療が終了するまでに要する時間を推定可能とすることで、残りの治療計画の精度を向上可能である。さらに計画に従って治療を進行するために、画像取得回数を低減して、迅速な熱治療ができる。特に X線を使用した画像を使用する場合は、必要回数だけ撮影するので被曝量を低減可能である。また、治療のモニタリング時に病巣と治療対象領域の位置合わせ、治療対象領域内の血管の存在、温度などの治療条件が実現を容易に把握することができるので、治療時に非常に有用な支援情報を提供できる。

Claims

請求の範囲

[1] 被検体の断層像を取得し治療経過をモニタリングする治療支援用画像処理装置において、

前記断層像上力モニタリング画像を選択し、その上に目的とする治療領域を設定する目的治療領域設定手段と、

前記モニタリング画像上に既に治療したと考えられる既治療領域を設定する既治療領域設定手段と、

前記目的治療領域設定手段により設定した目的治療領域の治療が終了する時間を前記既治療領域に基づいて推定する総治療時間推定手段と、

を備えたことを特徴とする治療支援用画像処理装置。

[2] 前記総治療時間推定手段は、複数のモニタリング画像内の既治療領域にっ、て、治療の開始地点である基準点力前記既治療領域の境界までの距離をそれぞれ求め、前記距離の変化に基づいて前記総治療時間を推定することを特徴とする請求項 1に記載の治療支援用画像処理装置。

[3] 前記総治療時間推定手段は、前記目的治療領域上、前記基準点から等角度間隔で放射状に伸びる直線を基準線とし、前記基準線上における治療進度をもとに前記総治療時間を推定することを特徴とする請求項 2に記載の治療支援用画像処理装置

[4] 前記総治療時間推定手段は、前記基準点から前記目的治療領域を二等分する少なくとも一本の基準線を指定し、前記基準線を所定間隔に分割する垂線が前記目的治療領域と交差する点を求め、これら点と前記基準点を結ぶ線をさらに基準線とし、前記基準線上における治療進度をもとに前記総治療時間を推定することを特徴とする請求項 2に記載の治療支援用画像処理装置。

[5] 前記総治療時間推定手段は、前記目的治療領域の輪郭を表示する前記表示手段上の画素中の所定間隔の画素と前記基準点を結ぶ線を基準線とし、前記基準線上における治療進度をもとに前記総治療時間を推定することを特徴とする請求項 2に記載の治療支援用画像処理装置。

[6] 前記目的治療領域設定手段は、複数の閉曲線を合成してから前記閉曲線の交差部分をスムージングして設定し、

前記総治療時間推定手段は、前記閉曲線毎に基準線を設定し、前記複数の基準線上における治療進度をもとに前記総治療時間を推定することを特徴とする請求項 3 から 5のいずれかに記載の治療支援用画像処理装置。

[7] 被検体の断層像を取得し治療経過をモニタリングする治療支援用画像処理装置において、

前記モニタリング画像に既に治療したと考えられる既治療領域を設定する既治療領域設定手段と、

上記既治療領域に基づいて、所定時間経過後の治療領域を予測し、この予想される治療領域を少なくとも 1つの色情報で表示する表示手段と、

を備えたことを特徴とする治療支援用画像処理装置。

[8] 前記表示手段は、前記目的治療領域と前記既治療領域と前記予測治療領域とにはそれぞれ異なる色情報を割り当てて表示することを特徴とする請求項 7に記載の治療支援用画像処理装置。

[9] 前記表示手段は、上記目的治療範囲に色情報を割り当てて表示することを特徴とする請求項 7に記載の治療支援用画像処理装置。

[10] 前記表示手段は、上記目的治療領域に色情報を割り当てて表示するとともに、治療中の部位の変化を色で表示することを特徴とする請求項 7又は 9に記載の治療支援用画像処理装置。

[11] 前記表示手段は、複数時点における上記既治療領域にそれぞれ異なる色情報を割り当てて表示することを特徴とする請求項 7又は 9に記載の治療支援用画像処理装

[12] 前記表示手段は、上記予測された治療領域に色情報を段階的に割り当てて表示することを特徴とする請求項 7、 9、 10又は 11のいずれかに記載の治療支援用画像処

[13] 前記表示手段は、複数時点における上記予測された治療領域をそれぞれ所望の色階調で表示することを特徴とする請求項 7、 9、 10又は 11のいずれかに記載の治療支援用画像処理装置。

[14] 前記表示手段は、前記予測された治療領域と、前記既治療領域とのうち少なくとも一つを、前記目的治療領域に重ね合わせて表示することを特徴とする請求項 7から 13のいずれかに記載の治療支援用画像処理装置。

[15] 被検体の断層像を取得し治療経過をモニタリングする治療支援用画像処理装置において、

前記断層像を作成するためのデータを記憶する画像メモリと、

前記断層像作成用のデータから計測値を求める計測値算出手段と、

前記計測値を温度に読み替えする計測値 -温度変換手段と、

前記読み替えされた温度に対して色情報を割り当て温度色表示をするカラーテーブルと、

前記画像メモリに記憶された断層像作成用のデータにグレースケール階調を割り当てるルックアップテーブルと、

前記温度色表示とグレースケール階調が割り振られた断層像作成用のデータとの合成を指示するモニタリングモードスィッチと、

前記合成された温度色表示と断層像作成用のデータを画像として表示する表示手段と、を備えたことを特徴とする治療支援用画像処理装置。

[16] 前記計測値温度変換手段は、 CT値を計測値とするものであり、 25度付近での CT 値が 0である組織力^度以下にぉ、て 20から 50の CT値を有するとして、 CT値を温度に読み替えることを特徴とする請求項 15に記載の治療支援用画像処理装置。

[17] 治療支援のための表示モードを選択する治療支援表示モード選択手段と、

表示したい治療進行イメージの治療開始後の時間を入力する時間選択入力手段と前記選択された時間における治療進行イメージを作成する治療進行イメージ作成手段と、

前記目的治療領域と前記既治療領域のうち少なくとも一方と前記治療進行イメージとを表示する表示手段と、を備えた請求項 1又は 2に記載の治療支援用画像処理装置。

[18] 前記目的治療領域と前記既治療領域のうち少なくとも一方と前記治療進行イメージとの間の差分画像を作成する差分画像作成手段と、を備え、

前記表示手段は、前記目的治療領域と前記既治療領域のうち少なくとも一方と前記差分画像とを表示する請求項 17に記載の治療支援用画像処理装置。

[19] 前記表示手段は、前記目的治療領域と、前記既治療領域と、前記治療進行ィメージと、前記差分画像とに、それぞれ異なる色情報を割り当てて表示することを特徴とする請求項 17又は 18に記載の治療支援用画像処理装置。