特開 2023-110463 温度分布情報取得装置およびプログラム

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】公開特許公報(A)

(11)【公開番号】P2023110463

(43)【公開日】2023-08-09

(54)【発明の名称】温度分布情報取得装置およびプログラム

(51)【国際特許分類】

   G01K 11/30 20060101AFI20230802BHJP

   A61B 5/00 20060101ALI20230802BHJP

   A61B 6/03 20060101ALI20230802BHJP

【ＦＩ】

G01K11/30

A61B5/00 G

A61B6/03 370A

A61B6/03 360Z

A61B6/03 377

【審査請求】未請求

【請求項の数】6

【出願形態】ＯＬ

(21)【出願番号】P 2022011936

(22)【出願日】2022-01-28

(71)【出願人】

【識別番号】598041566

【氏名又は名称】学校法人北里研究所

(74)【代理人】

【識別番号】100106909

【弁理士】

【氏名又は名称】棚井 澄雄

(74)【代理人】

【識別番号】100188558

【弁理士】

【氏名又は名称】飯田 雅人

(72)【発明者】

【氏名】原 秀剛

(72)【発明者】

【氏名】水上 慎也

【テーマコード（参考）】

4C093

4C117

【Ｆターム（参考）】

4C093AA22

4C093AA25

4C093AA26

4C093CA37

4C093FD11

4C093FF03

4C093FF08

4C093FF09

4C093FF28

4C093FF35

4C093FG01

4C117XA07

4C117XB01

4C117XD27

4C117XE23

4C117XE44

4C117XG22

4C117XG24

4C117XJ01

4C117XK05

4C117XR07

(57)【要約】

【課題】画像ノイズの影響を低減した温度分布情報を取得することができる温度分布情報取得装置を提供する。
【解決手段】温度分布情報取得装置は、ＣＴ画像から検査対象のＣＴ値のばらつきを算出する算出部と、前記検査対象のＣＴ値と温度との関係を表す変換テーブルを用いて、前記ＣＴ値のばらつきを温度のばらつきに変換し、当該温度のばらつきに基づいて前記検査対象についての温度分解能を設定する設定部と、前記変換テーブルを用いて前記ＣＴ画像における前記検査対象のＣＴ値を温度に変換する変換部と、設定した前記温度分解能に基づく温度範囲毎に異なる表示態様で、前記検査対象の前記温度の分布を画像化した温度分布情報を生成する生成部と、を備える。
【選択図】図２

【特許請求の範囲】

【請求項1】

ＣＴ画像から検査対象のＣＴ値のばらつきを算出する算出部と、
前記検査対象のＣＴ値と温度との関係を表す変換テーブルを用いて、前記ＣＴ値のばらつきを温度のばらつきに変換し、当該温度のばらつきに基づいて前記検査対象についての温度分解能を設定する設定部と、
前記変換テーブルを用いて前記ＣＴ画像における前記検査対象のＣＴ値を温度に変換する変換部と、
設定した前記温度分解能に基づく温度範囲毎に異なる表示態様で、前記検査対象の前記温度の分布を画像化した温度分布情報を生成する生成部と、
を備える温度分布情報取得装置。

【請求項2】

前記算出部は、前記検査対象内の関心領域の指定を受け付けて、当該関心領域内における前記ＣＴ値のばらつきを算出する、
請求項１に記載の温度分布情報取得装置。

【請求項3】

前記生成部は、前記温度範囲毎に異なる色を設定したルックアップテーブルを作成し、当該ルックアップテーブルを用いて前記温度分布情報を生成する、
請求項１または２に記載の温度分布情報取得装置。

【請求項4】

前記ＣＴ画像に画像ノイズを低減するための画像処理を行う画像処理部をさらに備え、
前記算出部は、前記画像処理後の前記ＣＴ画像から前記検査対象の前記ＣＴ値のばらつきを算出し、
前記変換部は、前記画像処理後の前記ＣＴ画像における前記検査対象のＣＴ値を前記温度に変換する、
請求項１から３の何れか一項に記載の温度分布情報取得装置。

【請求項5】

前記温度分布情報を前記ＣＴ画像に重畳して表示する表示部をさらに備える、
請求項１から４の何れか一項に記載の温度分布情報取得装置。

【請求項6】

ＣＴ画像から検査対象のＣＴ値のばらつきを算出するステップと、
前記検査対象のＣＴ値と温度との関係を表す変換テーブルを用いて、前記ＣＴ値のばらつきを温度のばらつきに変換し、当該温度のばらつきに基づいて前記検査対象についての温度分解能を設定するステップと、
前記変換テーブルを用いて前記ＣＴ画像における前記検査対象のＣＴ値を温度に変換するステップと、
設定した前記温度分解能に基づく温度範囲毎に異なる表示態様で、前記検査対象の前記温度の分布を画像化した温度分布情報を生成するステップと、
を温度分布情報取得装置に実行させるプログラム。

【発明の詳細な説明】

【技術分野】

【0001】

本開示は、温度分布情報取得装置およびプログラムに関する。

【背景技術】

【0002】

特許文献１には、Ｘ線ＣＴ装置が測定した対象物の各部（例えば筋肉および脂肪）のＣＴ値を温度に換算して、対象物の温度分布情報を取得する装置が開示されている。

【先行技術文献】

【特許文献】

【0003】

【特許文献1】特許第６４８３４７３号公報

【発明の概要】

【発明が解決しようとする課題】

【0004】

Ｘ線ＣＴ装置が撮影した画像にはノイズが多く含まれるため、その影響に起因するＣＴ値のバラツキ（標準偏差）が発生し、ＣＴ値から換算した温度値にもバラツキが反映される。特許文献１に記載の装置で取得した温度分布情報において、温度値に色付け（mapping）をした画像情報を生成した場合、画像ノイズの影響による実際の温度と異なる表示色が示される可能性がある。すなわち、各部位の温度値を正確に認識することが困難となる可能性がある。

【0005】

本開示は、このような課題に鑑みてなされたものであって、画像ノイズの影響を低減した温度分布情報を取得することができる温度分布情報取得装置およびプログラムを提供する。

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

【0006】

（１）本開示の第１の態様によれば、温度分布情報取得装置は、ＣＴ画像から検査対象のＣＴ値のばらつきを算出する算出部と、前記検査対象のＣＴ値と温度との関係を表す変換テーブルを用いて、前記ＣＴ値のばらつきを温度のばらつきに変換し、当該温度のばらつきに基づいて前記検査対象についての温度分解能を設定する設定部と、前記変換テーブルを用いて前記ＣＴ画像における前記検査対象のＣＴ値を温度に変換する変換部と、設定した前記温度分解能に基づく温度範囲毎に異なる表示態様で、前記検査対象の前記温度の分布を画像化した温度分布情報を生成する生成部と、を備える。

【0007】

（２）本開示の第２の態様によれば、第１の態様に係る温度分布情報取得装置において、前記算出部は、前記検査対象内の関心領域の指定を受け付けて、当該関心領域内における前記ＣＴ値のばらつきを算出する。

【0008】

（３）本開示の第３の態様によれば、第１または第２の態様に係る温度分布情報取得装置において、前記生成部は、前記温度範囲毎に異なる色を設定したルックアップテーブルを作成し、当該ルックアップテーブルを用いて前記温度分布情報を生成する。

【0009】

（４）本開示の第４の態様によれば、第１から第３の何れか一の態様に係る温度分布情報取得装置は、前記ＣＴ画像に画像ノイズを低減するための画像処理を行う画像処理部をさらに備え、前記算出部は、前記画像処理後の前記ＣＴ画像から前記検査対象の前記ＣＴ値のばらつきを算出し、前記変換部は、前記画像処理後の前記ＣＴ画像における前記検査対象のＣＴ値を前記温度に変換する。

【0010】

（５）本開示の第５の態様によれば、第１から第４の何れか一の態様に係る温度分布情報取得装置は、前記温度分布情報を前記ＣＴ画像に重畳して表示する表示部をさらに備える。

【0011】

（６）本開示の第６の態様によれば、プログラムは、ＣＴ画像から検査対象のＣＴ値のばらつきを算出するステップと、前記検査対象のＣＴ値と温度との関係を表す変換テーブルを用いて、前記ＣＴ値のばらつきを温度のばらつきに変換し、当該温度のばらつきに基づいて前記検査対象についての温度分解能を設定するステップと、前記変換テーブルを用いて前記ＣＴ画像における前記検査対象のＣＴ値を温度に変換するステップと、設定した前記温度分解能に基づく温度範囲毎に異なる表示態様で、前記検査対象の前記温度の分布を画像化した温度分布情報を生成するステップと、を温度分布情報取得装置に実行させる。

【発明の効果】

【0012】

本開示に係る温度分布情報取得装置およびプログラムによれば、画像ノイズの影響を低減した温度分布情報を取得することができる。

【図面の簡単な説明】

【0013】

【図1】本発明の温度分布情報取得システムの全体構成を示す図である。

【図2】本発明の温度分布情報取得装置のプロセッサの機能構成を説明するためのブロック図である。

【図3】本発明の温度分布情報取得装置の処理の一例を示す第１のフローチャートである。

【図4】本発明の温度分布情報取得装置の機能を説明するための第１の図である。

【図5】本発明の変換テーブルの一例を示す図である。

【図6】本発明の器官毎の温度分解能の一例を示す図である。

【図7】本発明の温度分布情報取得装置の処理の一例を示す第２のフローチャートである。

【図8】本発明の温度分布情報取得装置の機能を説明するための第２の図である。

【図9】本発明の変形例における温度分布情報取得装置のプロセッサの機能構成を説明するためのブロック図である。

【発明を実施するための形態】

【0014】

＜第１の実施形態＞
以下、本開示の第１の実施形態に係る温度分布情報取得システム１および温度分布情報取得装置２０について、図１～図８を参照しながら説明する。

【0015】

（温度分布情報取得システムの全体構成）
図１は、本発明の温度分布情報取得システムの全体構成を示す図である。
本発明に係る温度分布情報取得システム１は、検査対象の内部の温度分布を画像表示するシステムである。検査対象は、例えば、温熱療法、焼灼術、凍結療法などの治療を受ける患者などの生体である。なお、検査対象は人間に限られることはなく、他の実施形態では動物を検査対象としてもよい。また、さらに他の実施形態は、食品（例えば、発酵食品）などの非生体を検査対象としてもよい。なお、本実施形態では、検査対象が生体である態様を例として説明する。図１に示すように、温度分布情報取得システム１は、Ｘ線ＣＴ（Computed Tomography）装置１０と、温度分布情報取得装置２０とを備えている。

【0016】

Ｘ線ＣＴ装置１０は、生体にＸ線を照射し、Ｘ線透過量（生体を透過したＸ線の量）を測定する。Ｘ線ＣＴ装置１０は、複数の方向から生体にＸ線を照射して方向毎にＸ線透過量を測定し、測定結果を解析して生体の各部におけるＣＴ値を算出する。そして、Ｘ線ＣＴ装置１０は、算出したＣＴ値を検出可能なＣＴ画像を温度分布情報取得装置２０へ送信する。Ｘ線ＣＴ装置１０として、既存のＸ線ＣＴ装置を用いることができる。

【0017】

ＣＴ値は、Ｘ線を吸収する度合いを、空気の場合を－１０００とし、水の場合を０とした相対値で示す値であり、単位としてＨＵ（Hounsfield Unit）が用いられる。ＣＴ値を算出する方法は、例えば特許文献１に記載の方法を用いることができる。

【0018】

温度分布情報取得装置２０は、Ｘ線ＣＴ装置１０から受信したＣＴ値に基づいて、生体の温度分布情報を算出し、対象物の温度分布を画像にて表示する。図１に示すように、温度分布情報取得装置２０は、例えばプロセッサ２００、メモリ２０１、ストレージ２０２、インタフェース２０３、表示部２０４、および操作受付部２０５からなるコンピュータを用いて構成される。

【0019】

プロセッサ２００は、温度分布情報取得装置２０の動作全体を司る。

【0020】

メモリ２０１は、プロセッサ２００がプログラムに基づいて動作するための命令及びデータが展開される。

【0021】

ストレージ２０２は、いわゆる補助記憶装置であって、例えば、ＨＤＤ（Hard Disk Drive）、ＳＳＤ（Solid State Drive）などであってよい。

【0022】

インタフェース２０３は、Ｘ線ＣＴ装置１０と通信可能に接続するためのインタフェース（通信インタフェース）である。

【0023】

表示部２０４は、温度分布情報を表示するディスプレイである。

【0024】

操作受付部２０５は、温度分布情報取得装置２０の操作者による入力操作を受け付ける入力デバイスであり、例えば、一般的なマウス、キーボード、タッチセンサなどであってよい。

【0025】

（温度分布情報取得装置の機能構成）
図２は、本発明の温度分布情報取得装置のプロセッサの機能構成を説明するためのブロック図である。
図２に示すように、温度分布情報取得装置２０のプロセッサ２００は、所定のプログラムにしたがって動作することにより、算出部２００１、設定部２００２、変換部２００３、および生成部２００４としての機能を発揮する。

【0026】

算出部２００１は、ＣＴ画像から生体のＣＴ値のばらつきを算出する。

【0027】

設定部２００２は、生体のＣＴ値と温度との関係を表す変換テーブルを用いて、ＣＴ値のばらつきを温度のばらつきに変換し、当該温度のばらつきに基づいて生体についての温度分解能を設定する。

【0028】

変換部２００３は、変換テーブルを用いてＣＴ画像における生体のＣＴ値を温度に変換する。

【0029】

生成部２００４は、設定した温度分解能に基づく温度範囲毎に異なる表示態様で、生体の温度の分布を画像化した温度分布情報を生成する。生成部２００４は、例えば、温度範囲毎に異なる色を設定したルックアップテーブルを作成し、当該ルックアップテーブルを用いて温度分布情報を生成する。

【0030】

（温度分布情報取得装置の処理フロー）
図３は、本発明の温度分布情報取得装置の処理の一例を示す第１のフローチャートである。
まず、図３を参照しながら、温度分布情報取得装置２０による温度分解能設定処理の流れについて詳細に説明する。ここでは、温度分布情報取得装置２０が、温熱療法、焼灼術、凍結療法などの治療を受ける患者の体内（特に、治療を施す部位）の温度分布情報を取得する例について説明する。

【0031】

温度分布情報取得装置２０は、Ｘ線ＣＴ装置１０からＣＴ画像を取得する（ステップＳ０１）。このＣＴ画像は、温度分布情報取得装置２０の表示部２０４に表示される。

【0032】

次に、温度分布情報取得装置２０の操作受付部２０５は、操作者によるＣＴ画像上の対象領域Ｒ１および関心領域Ｒ２を指定する操作を受け付ける（ステップＳ０２）。

【0033】

図４は、本発明の温度分布情報取得装置の機能を説明するための第１の図である。
図４に示すように、表示部２０４には、Ｘ線ＣＴ装置１０から取得した患者のＣＴ画像が表示されている。操作者は、表示部２０４に表示されたＣＴ画像を見ながら、患者の温度分布情報を生成する対象領域Ｒ１、および関心領域Ｒ２を指定する操作を行う。対象領域Ｒ１は、例えば、治療を施す器官を示す領域である。関心領域Ｒ２は、対象領域Ｒ１内においてＣＴ値がより均一に近い領域であり、例えば、操作者がＣＴ画像を目視確認して、関心領域Ｒ２を設定する位置を決める。

【0034】

このとき、対象領域Ｒ１および関心領域Ｒ２の範囲は、操作者が操作受付部２０５（マウスなど）を使って手動で指定してよい。

【0035】

また、操作者がＣＴ画像上のポイントを指定する操作を行うと、温度分布情報取得装置２０の算出部２００１がそのポイントに基づいて対象領域Ｒ１および関心領域Ｒ２の範囲を自動的に設定するようにしてもよい。例えば、算出部２００１は、既知の画像処理を使ってＣＴ画像上の指定されたポイントの周囲に存在する器官（例えば、筋肉、脂肪、肝臓、腎臓など）の輪郭を抽出し、この輪郭内の領域を対象領域Ｒ１として設定する。また、算出部２００１は、ＣＴ画像上の指定されたポイントから所定距離の円形領域を関心領域Ｒ２として設定する。所定距離は、予め器官毎に固定値が設定されていてもよいし、操作者が任意の値を指定してもよい。また、関心領域は円形に限られることはなく、矩形など他の形状であってもよい。

【0036】

なお、対象領域Ｒ１は複数設定されてもよい。また対象領域Ｒ１が複数設定された場合、対象領域Ｒ１それぞれについて１箇所ずつ関心領域Ｒ２が設定される。

【0037】

次に、算出部２００１は、対象領域Ｒ１の関心領域Ｒ２内のＣＴ値のばらつきを算出する（ステップＳ０３）。ここでは、算出部２００１は、ＣＴ値の標準偏差ＳＤ（σ_ＣＴ）を算出する。

【0038】

なお、算出部２００１は、複数の対象領域Ｒ１（例えば、器官Ｔ１～Ｔ４の４つ）が設定されている場合、各対象領域Ｒ１について、ＣＴ値の標準偏差ＳＤ（σ_ＣＴ１～σ_ＣＴ４）を算出する。

【0039】

また、設定部２００２は、ＣＴ値のばらつきを、温度のばらつきに変換する（ステップＳ０４）。具体的には、設定部２００２は、以下の式（１）により、ＣＴ値の標準偏差ＳＤ（σ_ＣＴ）を温度の標準偏差ＳＤ（σ_Ｔ）に変換する。

【0040】

【数1】

【0041】

なお、式（１）において、ｆ´（Ｔ）はＣＴ値と温度の関係を表す関数ｆ（Ｔ）（以下、「変換テーブル」とも記載する。）の傾きである。

【0042】

図５は、本発明の変換テーブルの一例を示す図である。
図５は、生体の複数の器官Ｔ１～Ｔ４それぞれのＣＴ値と温度との関係を線形近似したグラフである。例えば、器官Ｔ１～Ｔ４は、それぞれ肝臓、筋肉、腎臓、脂肪である。これらのグラフは、事前に各器官Ｔ１～Ｔ４（例えば、肝臓、筋肉、腎臓、脂肪）のＣＴ値および温度の計測値に基づき求めておいたものであり、ＣＴ値から温度への変換に用いる変換テーブルとして、ストレージ２０２に記録されている。変換テーブルは、全ての患者に対し同じものを使用してよい。

【0043】

また、例えば性別、年齢、身長、体重などで分類したグループ毎に変換テーブルを用意しておき、患者が属するグループに対応する変換テーブルを使用してもよい。

【0044】

例えば、ステップＳ０２において、器官Ｔ１が対象領域Ｒ１として設定されたとする。この場合、設定部２００２は、器官Ｔ１の変換テーブルの傾きｆ´（Ｔ）を用いて、器官Ｔ１のＣＴ値の標準偏差ＳＤ（σ_ＣＴ）を、温度の標準偏差ＳＤ（σ_Ｔ）に変換する。

【0045】

なお、設定部２００２は、複数の対象領域Ｒ１（例えば、器官Ｔ１～Ｔ４の４つ）が設定されている場合、各対象領域Ｒ１について、ＣＴ値の標準偏差ＳＤ（σ_ＣＴ１～σ_ＣＴ４）を温度の標準偏差ＳＤ（σ_Ｔ１～σ_Ｔ４）に変換する。

【0046】

次に、設定部２００２は、温度の標準偏差ＳＤ（σ_Ｔ）に基づいて、対象領域Ｒ１における温度分解能を設定する（ステップＳ０５）。例えば、設定部２００２は、温度の標準偏差ＳＤ（σ_Ｔ）の２倍の値である、２σ_Ｔ（すなわち、±σ_Ｔの範囲）を対象領域Ｒ１の温度分解能として設定する。

【0047】

図６は、本発明の器官毎の温度分解能の一例を示す図である。
図６の例のように、複数の対象領域Ｒ１（例えば、器官Ｔ１～Ｔ４の４つ）が設定されている場合、設定部２００２は、各対象領域の温度分解能を設定する。

【0048】

図６の例では、設定部２００２は、算出部２００１が算出した器官Ｔ１のＣＴ値の標準偏差ＳＤ（σ_ＣＴ１＝１．３）と、器官Ｔ１の変換テーブルの傾き（－０．５）とを上記した式（１）に代入し、温度の標準偏差ＳＤ（σ_Ｔ１＝２．６）に変換する。また、設定部２００２は、器官Ｔ１の温度分解能を２σ_Ｔである「５．２」に設定する。同様に、設定部２００２は、他の対象領域Ｒ１（器官Ｔ２、Ｔ３、Ｔ４）の温度分解能をそれぞれ設定する。

【0049】

なお、１つのＣＴ画像に複数の対象領域Ｒ１が含まれる場合、設定部２００２は、これら複数の対象領域Ｒ１のうち最も温度分解能が大きい（悪い）ものに合わせて、ＣＴ画像全体に適用する１つの温度分解能を設定してもよい。図６の例では、設定部２００２は、器官Ｔ３の温度分解能（９．４）が最大となるため、この器官Ｔ３の温度分解能を切り上げた値（１０．０）をＣＴ画像全体の温度分解能として設定する。

【0050】

温度分布情報取得装置２０は、図３に示す一連の処理を、例えば患者の治療を行う直前に１回実施する。また、温度分布情報取得装置２０は、他の患者の治療を行う際には、改めて図３の一連の処理を実施して、患者に応じた温度分解能を設定する。

【0051】

また、患者の治療中にＸ線ＣＴ装置１０の撮影条件などを変更した場合、ＣＴ値のばらつきが変化する可能性がある。この場合、温度分布情報取得装置２０は、図３の一連の処理を再度実施して、この患者についての温度分解能を更新してもよい。

【0052】

図７は、本発明の温度分布情報取得装置の処理の一例を示す第２のフローチャートである。
次に、図７を参照しながら、温度分布情報取得装置２０による温度分布情報生成処理の流れについて詳細に説明する。図３の一連の処理によって温度分解能の設定が完了すると、温度分布情報取得装置２０は、治療時における対象領域の温度分布情報生成処理を実行する。

【0053】

まず、図７に示すように、温度分布情報取得装置２０は、Ｘ線ＣＴ装置１０からＣＴ画像を取得する（ステップＳ１１）。このＣＴ画像は、温度分布情報取得装置２０の表示部２０４に表示される。

【0054】

次に、変換部２００３は、対象領域Ｒ１の変換テーブルｆ（Ｔ）を用いて、対象領域Ｒ１のＣＴ値を温度に変換する（ステップＳ１２）。

【0055】

次に、生成部２００４は、設定部２００２が設定した温度分解能に基づく温度範囲毎に異なる色を割り当てたルックアップテーブルＤ１（図８）を作成し、このルックアップテーブルＤ１を用いて対象領域Ｒ１の温度の分布を色分けした温度分布情報Ｄ２（図８）を生成する（ステップＳ１３）。

【0056】

図８は、本発明の温度分布情報取得装置の機能を説明するための第２の図である。
図８に示すように、ルックアップテーブルＤ１は、所定の温度範囲を温度分解能２σ_Ｔ毎の階調に分割し、階調毎に異なる色を設定したテーブルである。図８の例では、温度範囲は－２０℃～６０℃であり、温度分解能は１０．０℃である。したがって、ルックアップテーブルＤ１は８つの階調に分割され、それぞれに異なる色が設定される。図８に示すように、ルックアップテーブルＤ１は、操作者が色と温度の対応を理解できるように、表示部２０４に表示されてもよい。また、ルックアップテーブルＤ１には、ＣＴ画像を撮影したときの患者の体温を示す情報が含まれていてもよい。

【0057】

生成部２００４は、ルックアップテーブルＤ１を参照し、対象領域Ｒ１の各部（各ピクセル）を、温度に対応する色で塗った画像である温度分布情報Ｄ２を生成する。

【0058】

また、表示部２０４は、ＣＴ画像上に、生成部２００４が生成した温度分布情報Ｄ２を重畳表示する（ステップＳ１４）。これにより、操作者は、ＣＴ画像上において、対象領域Ｒ１（患者の治療部位）の温度分布を容易に把握することができる。

【0059】

なお、温度分布情報取得装置２０は、例えば患者に温熱療法を実施する前、実施中、および実施後の各時点でＣＴ画像を撮影する度に、図７に示す一連の処理を実施するようにしてもよい。これにより、操作者は、治療部位が適切な温度となっているかを確認しながら、治療を行うことができる。

【0060】

（温度分布情報取得装置の作用効果）
以上のように、本発明に係る温度分布情報取得装置２０は、ＣＴ画像における生体（検査対象）のＣＴ値のばらつきを温度のばらつきを求め、この温度のばらつきに基づいて生体についての温度分解能を設定する。また、温度分布情報取得装置２０は、温度分解能に基づく温度範囲毎に異なる表示態様（色）で、生体のＣＴ値から変換した温度の分布を画像化した温度分布情報を生成する。

【0061】

このようにすることで、温度分布情報取得装置２０は、画像ノイズに起因する温度のばらつきを低減した温度分布情報を取得することができる。

【0062】

また、温度分布情報取得装置２０は、この温度分布情報をＣＴ画像に重畳表示する。

【0063】

このようにすることで、温度分布情報取得装置２０は、温熱療法などの治療により、生体内の温度分布が適切な状態になっているかを、操作者に容易に把握させることができる。

【0064】

＜変形例＞
図９は、本発明の変形例における温度分布情報取得装置のプロセッサの機能構成を説明するためのブロック図である。
図９に示すように、変形例に係る温度分布情報取得装置２０のプロセッサ２００は、画像処理部２００５としての機能をさらに発揮してもよい。

【0065】

画像処理部２００５は、Ｘ線ＣＴ装置１０から受信したＣＴ画像に画像ノイズを低減させるための画像処理を行う。画像処理は、例えば平滑化処理である。

【0066】

本発明の変形例において、画像処理部２００５は、図３のステップＳ０１において、Ｘ線ＣＴ装置１０から取得したＣＴ画像に対し、画像ノイズを低減させるための画像処理（平滑化処理など）を行う。この場合、図３のステップＳ０３において、算出部２００１は、画像処理後のＣＴ画像に基づいてＣＴ値のばらつきを算出する。

【0067】

同様に、画像処理部２００５は、図７のステップＳ１１において、Ｘ線ＣＴ装置１０から取得したＣＴ画像に対し、画像ノイズを低減させるための画像処理を行う。この場合、図７のステップＳ１２において、変換部２００３は、画像処理後のＣＴ画像に基づいてＣＴ値から温度への変換を行う。

【0068】

このようにすることで、温度分布情報取得装置２０は、より画像ノイズの影響を低減した温度分布情報を取得することができる。

【0069】

以上のとおり、本開示に係るいくつかの実施形態を説明したが、これら全ての実施形態は、例として提示したものであり、発明の範囲を限定することを意図していない。これらの実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で種々の省略、置き換え、変更を行うことができる。これらの実施形態及びその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれる。

【符号の説明】

【0070】

１ 温度分布情報取得システム
１０ Ｘ線ＣＴ装置
２０ 温度分布情報取得装置
２００ プロセッサ
２００１ 算出部
２００２ 設定部
２００３ 変換部
２００４ 生成部
２００５ 画像処理部
２０１ メモリ
２０２ ストレージ
２０３ インタフェース
２０４ 表示部
２０５ 操作受付部

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】