JP6128691B2 - 医用画像計測装置および方法並びにプログラム - Google Patents

Description

本発明は、医用画像内に存在する注目対象の形状の特徴を計測する医用画像計測装置および方法並びにプログラムに関するものである。

従来、放射線画像、ＭＲＩ（magnetic resonance imaging）画像または超音波画像などの医用画像内に含まれる病変または骨や臓器などの注目対象の所定の特徴量を計測し、その計測結果に基づいて、病変の成長率や骨や臓器の異常の有無を確認することが行われている。

たとえば特許文献１においては骨を撮影した放射線画像を用いて骨塩量を計測することが提案されており、特許文献２においては肺を撮影した放射線画像を用いて、肺の小領域の面積の変化を求めることによって異常の有無を判定することが提案されており、特許文献３においては超音波画像を用いて腫瘍の２点間の距離や面積を計測することが提案されている。

特開２０１４−５４４７６号公報 特開２０１３−３９４２７号公報 特開２０１３−１１１４３４号公報

ここで、たとえば上述したような医用画像内に含まれる病変の成長率を確認する場合、過去に撮影した医用画像に含まれる病変の最大径を計測するとともに、今回撮影した医用画像に含まれる病変の最大径を計測し、これらの最大径の長さの変化を算出することによって病変の成長率を確認することが考えられる。

しかしながら、腫瘍などの病変は、円形に近い状態から楕円形へと形を変化させながら成長する場合がある。具体的には、図２に示すように、第１の医用画像内の病変のように円形に近い状態から第２の医用画像内の病変のように所定の方向に延びた楕円に近い状態に変化する場合がある。

このような場合にも上述したように第１の医用画像内の病変の最大径（たとえば点Ｐａ１と点Ｐａ２とによって特定される径）と、第２の医用画像内の病変の最大径（たとえば点Ｐｂ１と点Ｐｂ２によって特定される径）とを計測し、これらの変化量を算出したのでは、腫瘍の正確な成長率を計測できない場合がある。すなわち、腫瘍などの経時変化を考慮しなければ正確な成長率を計測できない場合がある。なお、ここでは病変が拡大する例で説明しているが、病変が縮小する場合についても同様である。

本発明は、上記事情に鑑み、病変などの注目対象の経時変化も考慮した医用画像計測を行うことができる医用画像計測装置および方法並びにプログラムを提供することを目的とするものである。

本発明の医用画像計測装置は、同一患者における同一の注目対象を異なる時点で撮影した第１の医用画像および第２の医用画像を取得する医用画像取得部と、第１の医用画像における注目対象の形状の特徴を計測するために設定された第１の計測パラメータおよび第２の医用画像における注目対象の形状の特徴を計測するために設定された第２の計測パラメータを取得する計測パラメータ取得部と、第１の計測パラメータと第２の計測パラメータとの間の変化を示す評価値を取得する評価値取得部と、その評価値に基づいて、上記変化が予め設定された変化量以上であるか否かを判定する判定部とを備えたことを特徴とする。

また、上記本発明の医用画像計測装置においては、判定部によって上記変化が予め設定された変化量以上であると判定された場合、第１の計測パラメータおよび第２の計測パラメータのうちのいずれか一方の計測パラメータを、他方の計測パラメータを用いて再設定することによって第３の計測パラメータを取得する計測パラメータ再設定部と、第３の計測パラメータを用いて、上記一方の計測パラメータが設定された第１の医用画像または第２の医用画像の注目対象の形状の特徴を計測する計測部とを備えることができる。

また、計測パラメータ再設定部は、第１の医用画像における注目対象の画像と第２の医用画像における注目対象の画像との非剛体レジストレーションの結果および上記他方の計測パラメータに基づいて、第３の計測パラメータを取得することができる。

また、第１の計測パラメータ、第２の計測パラメータおよび第３の計測パラメータとして、計測される注目対象の径を特定する情報を用いることができる。

また、第１の計測パラメータ、第２の計測パラメータおよび第３の計測パラメータとして、注目対象における２点間の距離を計測するために用いられるその２点の座標を用いることができる。

また、計測パラメータ再設定部は、上記他方の計測パラメータである２点の座標を結ぶ直線と、上記一方の計測パラメータが設定された第１の医用画像または第２の医用画像における注目対象の輪郭との交点の座標を第３の計測パラメータとして取得することができる。

また、評価値取得部は、上記他方の計測パラメータである２点の座標を結ぶ直線と、上記一方の計測パラメータである２点の座標のそれぞれとの距離に基づいて評価値を取得することができる。

また、評価値取得部は、第１の計測パラメータである２点の座標を結ぶ直線と、第２の計測パラメータである２点の座標を結ぶ直線とがなす角度の情報を評価値として取得することができる。

また、第１の医用画像と第２の医用画像として３次元画像を取得し、第１の計測パラメータ、第２の計測パラメータおよび第３の計測パラメータとして、注目対象の異なる断面の径を特定する情報を用いることができる。

また、注目対象は、時間の経過によって形状が変化する病変とすることができる。

また、判定部によって上記変化が予め設定された変化量以上であると判定された場合、警告を行う警告部を設けることができる。

本発明の医用画像計測方法は、同一患者における同一の注目対象を異なる時点で撮影した第１の医用画像および第２の医用画像を取得し、第１の医用画像における注目対象の形状の特徴を計測するために設定された第１の計測パラメータおよび第２の医用画像における注目対象の形状の特徴を計測するために設定された第２の計測パラメータを取得し、第１の計測パラメータと第２の計測パラメータとの間の変化を示す評価値を取得し、その評価値に基づいて、上記変化が予め設定された変化量以上であるか否かを判定することを特徴とする。

本発明の医用画像計測プログラムは、コンピュータを、同一患者における同一の注目対象を異なる時点で撮影した第１の医用画像および第２の医用画像を取得する医用画像取得部と、第１の医用画像における注目対象の形状の特徴を計測するために設定された第１の計測パラメータおよび第２の医用画像における注目対象の形状の特徴を計測するために設定された第２の計測パラメータを取得する計測パラメータ取得部と、第１の計測パラメータと第２の計測パラメータとの間の変化を示す評価値を取得する評価値取得部と、その評価値に基づいて、上記変化が予め設定された変化量以上であるか否かを判定する判定部として機能させることを特徴とする。

本発明の医用画像計測装置および方法並びにプログラムによれば、同一患者における同一の注目対象を異なる時点で撮影した第１の医用画像および第２の医用画像を取得し、第１の医用画像における注目対象の形状の特徴を計測するために設定された第１の計測パラメータおよび第２の医用画像における注目対象の形状の特徴を計測するために設定された第２の計測パラメータを取得し、第１の計測パラメータと第２の計測パラメータとの間の変化を示す評価値を取得し、その評価値に基づいて、上記変化が予め設定された変化量以上であるか否かを判定するようにしたので、たとえば上記変化が予め設定された変化量以上である場合には、計測パラメータを自動的に再設定したり、またはその旨を警告することによってユーザが計測パラメータを再設定したりすることによって、より適切な計測を行うことができる。

本発明の医用画像計測装置および方法並びにプログラムの一実施形態を用いた医用画像診断支援システムの概略構成を示すブロック図 過去に撮影した第１の医用画像と最も直近に撮影された第２の医用画像の一例を示す図 第１の計測パラメータと第１の計測パラメータとの間の変化を示す評価値の取得方法の一例を説明するための図 第３の計測パラメータの取得方法の一例を説明するための図 本発明の医用画像計測装置および方法並びにプログラムの一実施形態を用いた医用画像診断支援システムの作用を説明するためのフローチャート ３次元画像内の病変領域に対する第１および第２の計測パラメータの設定方法の一例を説明するための図 ポリープの高さを計測する例を説明するための図 本発明の医用画像計測装置および方法並びにプログラムのその他の実施形態を用いた医用画像診断支援システムの概略構成を示すブロック図

以下、本発明の医用画像計測装置および方法並びにプログラムの一実施形態を用いた医用画像診断支援システムについて、図面を参照しながら詳細に説明する。図１は、本実施形態の医用画像診断支援システムの概略構成を示すブロック図である。

本実施形態の医用画像診断支援システムは、図１に示すように、医用画像計測装置１と、医用画像保管サーバ２と、ディスプレイ３と、入力装置４とを備えている。

医用画像計測装置１は、コンピュータに本実施形態の医用画像計測プログラムがインストールされたものである。

医用画像計測装置１は、中央処理装置（ＣＰＵ（central processing unit））、半導体メモリおよびハードディスクやＳＳＤ(Solid State Drive)等のストレージデバイスを備えている。ストレージデバイスには、本実施形態の医用画像計測プログラムがインストールされており、この医用画像計測プログラムが中央処理装置によって実行されることによって、図１に示すような医用画像取得部１０、計測部１１、計測パラメータ取得部１２、評価値取得部１３、判定部１４、計測パラメータ再設定部１５および表示制御部１６が動作する。

医用画像取得部１０は、予め撮影された患者の医用画像５を取得するものである。医用画像５としては、たとえばＣＴ（Computed Tomography）装置やＭＲＩ（Magnetic Resonance Imaging）装置などによって撮影された断層画像や、複数の断層画像からなる３次元画像などがある。断層画像としては、アキシャル断面画像、サジタル断面画像、コロナル断面画像およびその他のＭＰＲ（multi-planar reconstruction）断面画像などがある。

医用画像５は、医用画像保管サーバ２に患者の識別情報とともに予め保管されており、医用画像取得部１０は、入力装置４を用いてユーザによって入力された患者の識別情報に基づいて、その識別情報を有する医用画像５を医用画像保管サーバ２から読み出すものである。特に、本実施形態の医用画像取得部１０は、同一患者における同一の注目対象を異なる時点で撮影した第１の医用画像と第２の医用画像とを取得するものである。注目対象としては、たとえば時間の経過によって形状が変化する腫瘤やポリープや頭部の脳動脈瘤などがあるが、これらに限らずその他の病変でもよい。また、病変に限らず、特定の骨や関節などでもよい。また、本実施形態の医用画像取得部１０は、第１の医用画像として過去に撮影された医用画像を取得し、第２の医用画像として最も直近に撮影された医用画像を取得するものとする。

第１の医用画像と第２の医用画像は、たとえばユーザが入力装置４を用いてファイル名を入力したり、ディスプレイ３に表示されたサムネイル画像などを選択することによって指定するようにすればよい。

また、たとえば医用画像５に対して頭部、胸部および腹部などといった撮影部位情報を付加しておき、医用画像取得部１０が、入力装置４を用いてユーザによって入力された撮影部位情報に基づいて、同一の撮影部位情報を有する第１の医用画像と第２の医用画像とを医用画像保管サーバ２から自動的に読み出すようにしてもよいし、医用画像５に対して腫瘤やポリープや脳動脈瘤などといった病変情報を付加しておき、入力装置４を用いてユーザによって入力された病変情報に基づいて、同一の病変情報を有する第１の医用画像と第２の医用画像とを医用画像保管サーバ２から読み出すようにしてもよい。また、このような方法に限らず、その他の公知な手法を用いて第１の医用画像と第２の医用画像とを自動的に取得するようにしてもよい。

計測部１１は、医用画像取得部１０によって取得された第１の医用画像および第２の医用画像が入力され、その入力された第１の医用画像および第２の医用画像に含まれる注目対象の形状の特徴を計測するものである。

具体的には、本実施形態の計測部１１は、第１の医用画像および第２の医用画像として断層画像を取得し、各断層画像内の上述した腫瘤などの病変領域の大きさを示す情報として、その病変領域の径を計測するものである。

過去に撮影された第１の医用画像に含まれる病変領域の径と、最も直近に撮影された第１の医用画像に含まれる病変領域の径とを比較することによって病変の成長率または放射線治療などによる病変の縮小効果を知ることができる。

計測パラメータ取得部１２は、計測部１１において第１の医用画像の病変領域の径を計測するために設定された第１の計測パラメータと、第２の医用画像の病変領域の径を計測するために設定された第２の計測パラメータとを取得するものである。具体的には、本実施形態の計測パラメータ取得部１２は、第１の医用画像の病変領域の径の２つの端点の座標を第１の計測パラメータとして取得し、第２の医用画像の病変領域の径の２つの端点の座標を第２の計測パラメータとして取得するものである。なお、本実施形態においては、第１および第２の計測パラメータとして径の２つの端点の座標を取得するようにしたが、端点の座標に限らず、たとえば径を表すベクトル情報を第１および第２の計測パラメータとして取得するようにしてもよい。

以下、第１の計測パラメータと第２の計測パラメータについて、図２を参照しながら具体的に説明する。図２は、同一の病変領域であるが大きさが異なる病変領域を含む第１の医用画像および第２の医用画像の一例を示すものである。

上述したように、まず、第１の医用画像に含まれる病変領域の径を計測するため、その径の端点（Ｐａ１，Ｐａ２）が設定される。径の端点の設定方法としては、たとえばディスプレイ３に表示された第１の医用画像上においてユーザが入力装置４を用いて指定するようにしてもよいし、または、自動抽出された病変領域に基づいて自動的に最大径の端点を特定するようにしてもよい。自動的に最大径の端点を特定する方法としては、既に公知な方法を用いることができる。

また、第２の医用画像に含まれる病変領域の径を計測するため、その径の端点（Ｐｂ１，Ｐｂ２）が設定される。径の端点の設定方法は、第１の医用画像の場合と同様である。

ここで、上述したように病変領域の径を計測して病変領域の成長率などを確認する場合、一般的には病変領域の最大径と思われる径がユーザによって指定されたり、自動的に設定されたりするが、上述したように病変の形の経時変化の仕方によっては、最大径同士を比較したのでは適切な成長率や縮小効果を確認できない場合がある。

そこで、本実施形態においては、計測部１１において第１の医用画像と第２の医用画像の病変領域の径を計測するために設定された第１の計測パラメータと第２の計測パラメータとの間の変化を示す評価値を取得し、その評価値に基づいて、設定された第１の計測パラメータと第２の計測パラメータとが病変の成長率や縮小効果を確認するために適切なものであるか否かを判定する。

評価値取得部１３は、上述した第１の計測パラメータと第２の計測パラメータとの間の変化を示す評価値を取得するものである。以下、図３を参照しながら、本実施形態の評価値取得部１３における評価値の取得方法を説明する。

本実施形態の評価値取得部１３は、具体的には、図３に示すように、第１の医用画像上で設定された端点Ｐａ１と端点Ｐａ２の座標を第２の医用画像上にプロットする。そして、第２の医用画像上に設定された端点Ｐｂ１と端点Ｐｂ２とを通る直線と端点Ｐａ１との距離値ｖ１を下式（１）、上記直線と端点Ｐａ２との距離値ｖ２を下式（２）に基づいて算出する。なお、式（１）および（２）における「×」は、外積を示す記号である。次に、式（１）および（２）によって算出された距離値ｖ１および距離ｖ２に基づいて、下式（３）および（４）によって評価値Ｖを算出する。すなわち、距離値ｖ１と距離値ｖ２を加算した値がゼロでない場合には、下式（３）によって評価値Ｖを算出し、距離値ｖ１と距離値ｖ２を加算した値がゼロである場合には、下式（４）によって評価値Ｖを算出する。なお、下式（４）におけるＶｍｉｎは、予め設定された閾値である。

判定部１４は、上述したようにして評価値取得部１３において算出された評価値Ｖを取得し、この評価値Ｖに基づいて、第１の計測パラメータと第２の計測パラメータとの間の変化が予め設定された変化量以上であるか否かを判定するものである。

具体的には、本実施形態の判定部１４は、評価値Ｖと予め設定された閾値Ｖｍｉｎとを比較し、評価値Ｖが閾値Ｖｍｉｎ以下である場合には、すなわち第２の計測パラメータである端点Ｐｂ１と端点Ｐｂ２を通る直線付近に第１の計測パラメータである端点Ｐａ１と端点Ｐａ２とが存在しない場合には、第１の計測パラメータと第２の計測パラメータとの間の変化が予め設定された変化量以上であると判定する。

一方、評価値Ｖと閾値Ｖｍｉｎとを比較し、評価値Ｖが閾値Ｖｍｉｎより大きい場合には、すなわち第２の計測パラメータである端点Ｐｂ１と端点Ｐｂ２を通る直線付近に、第１の計測パラメータである端点Ｐａ１と端点Ｐａ２とが存在する場合には、第１の計測パラメータと第２の計測パラメータとの間の変化が予め設定された変化量より小さいと判定する。

そして、判定部１４において、第１の計測パラメータと第２の計測パラメータとの間の変化が予め設定された変化量より小さいと判定された場合には、第１の計測パラメータと第２の計測パラメータとを用いてそれぞれ計測された病変領域の径は、病変領域の成長率や縮小効果を確認するものとして適切なものであるとし、これらの径が表示制御部１６に出力され、表示制御部１６によってディスプレイ３に表示される。

一方、判定部１４において、第１の計測パラメータと第２の計測パラメータとの間の変化が予め設定された変化量以上であると判定された場合には、第１の計測パラメータと第２の計測パラメータとを用いてそれぞれ計測された病変領域の径は、病変領域の成長率や縮小効果を確認するものとして不適切なものであるとし、その判定結果は、計測パラメータ再設定部１５に出力される。

計測パラメータ再設定部１５は、入力された判定結果に応じて、第１の医用画像に設定された第１の計測パラメータ（端点Ｐａ１，Ｐａ２）を、第２の計測パラメータ（端点Ｐｂ１，Ｐｂ２）を用いて再設定することによって第３の計測パラメータを取得するものである。

具体的には、本実施形態の計測パラメータ再設定部１５は、まず、図４に示すように第１の医用画像に含まれる病変領域を自動抽出する。なお、病変領域の抽出方法については、既に公知な方法を用いることができる。次に、第２の医用画像上において設定された病変領域の径の端点Ｐｂ１と端点Ｐｂ２を第１の医用画像上にプロットし、その端点Ｐｂ１と端点Ｐｂ２とを通る直線上の点のうち、第１の医用画像上における病変領域の輪郭上の点の座標を第３の計測パラメータとして取得する。すなわち、図４に示す点ＰＡ１と点ＰＡ２の座標を第３の計測パラメータとして取得する。

計測パラメータ再設定部１５において再設定された第３の計測パラメータは、計測部１１に出力される。計測部１１は、入力された点ＰＡ１と点ＰＡ２とを結ぶ径を第１の医用画像の病変領域の大きさを示す情報として計測し、その計測結果を表示制御部１６に出力する。表示制御部１６は、計測部１１において再計測された測定結果を、第１の医用画像の病変領域の大きさを示す情報としてディスプレイに表示させる。

表示制御部１６は、医用画像取得部１０によって取得された第１および第２の医用画像や第１および第２の医用画像に含まれる病変領域の径の計測結果などをディスプレイ３に表示させるものである。

入力装置４は、ユーザによる種々の設定入力を受け付けるものであり、キーボードやマウスなどから構成されるものである。

次に、本実施形態の医用画像計測装置の作用について、図５に示すフローチャートを参照しながら説明する。

まず、ユーザによる患者の識別情報などの入力に基づいて、医用画像取得部１０によって医用画像保管サーバ２から第１の医用画像および第２の医用画像が読み出されて取得される（Ｓ１０）。

医用画像取得部１０によって取得された第１の医用画像および第２の医用画像は、表示制御部１６に出力されてディスプレイ３に表示されるとともに計測部１１に出力される。そして、第１の医用画像に含まれる病変領域について手動または自動で第１の計測パラメータが設定され、第２の医用画像に含まれる病変領域について手動または自動で第２の計測パラメータが設定され、計測パラメータ取得部１２によってその設定された第１および第２の計測パラメータが取得される（Ｓ１２）。

計測パラメータ取得部１２によって取得された第１および第２の計測パラメータは、評価値取得部１３に出力され、評価値取得部１３は、入力された第１および第２の計測パラメータに基づいて上述した評価値Ｖを算出する（Ｓ１４）。

次いで、評価値取得部１３によって算出された評価値Ｖは、判定部１４に出力され、判定部Ｖは、上述したように評価値Ｖが閾値Ｖｍｉｎ以下であるか否かを判定することによって、第１の計測パラメータと第２の計測パラメータとの間の変化が予め設定された変化量以上であるか否かを判定する（Ｓ１６）。

そして、判定部１４において、第１の計測パラメータと第２の計測パラメータとの間の変化が予め設定された変化量より小さいと判定された場合には（Ｓ１８，ＮＯ）、Ｓ１２において設定された第１の計測パラメータを用いて第１の医用画像の病変領域の径を計測した計測結果と、第２の計測パラメータを用いて第２の医用画像の病変領域の径を計測した計測結果とが表示制御部１６に出力され、これらの計測結果がディスプレイ３に表示される（Ｓ２４）。

一方、判定部１４において、第１の計測パラメータと第２の計測パラメータとの間の変化が予め設定された変化量以上と判定された場合には（Ｓ１８，ＹＥＳ）、その判定結果が計測パラメータ再設定部１５に出力され、計測パラメータ再設定部１５は、入力された判定結果に応じて、上述したようにして第３の計測パラメータを取得する（Ｓ２０）。

そして、計測パラメータ再設定部１５によって取得された第３の計測パラメータは計測部１１に出力され、計測部１１は、第３の計測パラメータを用いて第１の医用画像の病変領域の径を計測し、その計測結果を表示制御部１６に出力する（Ｓ２２）。

表示制御部１６は、第１の医用画像の病変領域については、第３の計測パラメータを用いて測定された測定結果を表示させ、第２の医用画像の病変領域については、第２の計測パラメータを用いて測定された測定結果を表示させる（Ｓ２４）。

上記実施形態の医用画像診断支援システムによれば、第１の計測パラメータと第２の計測パラメータとの間の変化を示す評価値を取得し、その評価値に基づいて、上記変化が予め設定された変化量以上であるか否かを判定し、上記変化が予め設定された変化量以上である場合には、注目対象の経時変化を考慮した第３の計測パラメータを自動的に再設定するようにしたので、より適切な計測を行うことができる。

なお、第１の計測パラメータの設定およびその第１の計測パラメータを用いた計測と、第２の計測パラメータの設定およびその第２の計測パラメータを用いた計測とは、第１の医用画像と第２の医用画像との両方が撮影された時点で行うようにしてもよいし、第１の医用画像が撮影された時点で第１の計測パラメータの設定およびその第１の計測パラメータを用いた計測を行い、その後、第２の医用画像が撮影されて時点で第２の計測パラメータの設定およびその第２の計測パラメータを用いた計測を行うようにしてもよい。

また、上記実施形態の医用画像計測装置１においては、図４に示すように、第１の医用画像上に第２の医用画像における病変領域の径の端点Ｐｂ１と端点Ｐｂ２をプロットし、その端点Ｐｂ１と端点Ｐｂ２とを通る直線上の点のうち、第１の医用画像における病変領域の輪郭上の点の座標を第３の計測パラメータとして取得するようにしたが、第３の計測パラメータを取得する方法としては、この方法に限らず、その他の方法を用いるようにしてもよい。

具体的には、たとえば第１の医用画像における病変領域の画像Ｉａと第２の医用画像における病変領域の画像Ｉｂとを抽出し、画像Ｉａと画像Ｉｂとで非剛体レジストレーションを行い、その結果得られた画像変形量Ｗ（Ｐ）（Ｐは画像上の位置）を用いて、下式（５）によって第３の計測パラメータである端点ＰＡ１，ＰＡ２を取得するようにしてもよい。

また、上記実施形態においては、第２の医用画像上に設定された第２の計測パラメータである端点Ｐｂ１と端点Ｐｂ２とを通る直線と、第１の計測パラメータである端点Ｐａ１および端点Ｐａ２とのそれぞれの距離に基づいて評価値を算出するようにしたが、評価値としてはこれに限らず、たとえば第２の計測パラメータである端点Ｐｂ１と端点Ｐｂ２とを通る直線と、第１の計測パラメータである端点Ｐａ１と端点Ｐａ２とを通る直線とがなす角度の情報を評価値として算出するようにしてもよい。この場合、上記角度の情報が大きいほど第１の計測パラメータと第２の計測パラメータとの間の変化量が大きいことになる。

また、上記実施形態の医用画像計測装置１によって取得される第１の医用画像および第２の医用画像は、たとえばこれらの画像の両方がアキシャル断面画像であるような場合には、これらの画像の断面位置は同じ位置であることが望ましいが、すなわち第１の計測パラメータ（端点Ｐａ１，Ｐａ２）と第２の計測パラメータ（端点Ｐｂ１，Ｐｂ２）は、病変領域における同じ断面位置の径を特定するものであることが望ましいが、第１の医用画像および第２の医用画像として３次元画像を取得するような場合には、図６に示すように、第１の医用画像において設定される第１の計測パラメータ（端点Ｐａ１，Ｐａ２）と第２の医用画像において設定される第２の計測パラメータ（端点Ｐｂ１，Ｐｂ２）は、病変領域の異なる断面の径を特定するものでもよい。なお、この場合には、上述したように第１の医用画像における病変領域の画像Ｉａと第２の医用画像における病変領域の画像Ｉｂとの非剛体レジストレーションの結果を用いて第３の計測パラメータを取得するようにすればよい。なお、非剛体レジストレーションは、たとえばDavid Mattes, David R. Haynor, Hubert Vesselle, Thomas K. Lewellen, William Eubank, “Nonrigid multimodality image registration”, Proceedings of the SPIE, volume 4322, p.1609-1620, 2001や、特開２０１１-４１６５６等の公知の手法を用いればよい。

また、上記実施形態においては、病変の径を計測するようにしたが、計測の対象は必ずしも径でなくてもよく、たとえば病変が大腸や胃に形成されたポリープである場合には、そのポリープの大腸表面からの高さを計測するようにしてもよい。

この場合、たとえば、図７に示すように、過去に撮影された第１の医用画像内に存在するポリープの高さＨｃを計測するために設定された計測位置Ｐｃを第１の計測パラメータとして取得するとともに、今回撮影された第２の医用画像内のポリープの高さＨｄを計測するために設定された計測位置Ｐｄを第２の計測パラメータとして取得し、計測位置Ｐｃと計測位置Ｐｄとの距離Ｄを評価値として取得するようにすればよい。

そして、その評価値が予め設定された閾値以上である場合には、計測パラメータの変化が予め設定された変化量以上であると判定し、その場合には、第１の医用画像に設定された第１の計測パラメータである計測位置Ｐｃを第２の医用画像に設定された第２の計測パラメータである計測位置Ｐｄと同じに位置に再設定して第３の計測パラメータとし、計測位置Ｐｄにおける第１の医用画像のポリープの高さＨＤを再計測するようにすればよい。なお、上記説明においてポリープの高さとは、大腸表面から最も高い位置の高さのことをいう。

また、上記実施形態においては、評価値に基づいて、第１の計測パラメータと第２の計測パラメータとの間の変化が、予め設定された変化量以上であるか否かを判定し、予め設定された変化量以上である場合には、第３の計測パラメータを再設定するようにしたが、図８に示すように警告部１７をさらに設け、予め設定された変化量以上である場合に、表示制御部１６によってユーザに対してメッセージを表示したり、警告音を鳴らしたりすることによって警告を行うようにしてもよい。また、予め設定された変化量以上である場合、第３の計測パラメータを再設定することなく、警告のみ行うようにしてもよい。ユーザは警告を受けることによって計測パラメータを再設定することができる。

また、上記実施形態においては、第１の医用画像上に第２の計測パラメータ（Ｐｂ１，Ｐｂ２）をプロットし、その第２の計測パラメータ（Ｐｂ１，Ｐｂ２）を用いて第１の計測パラメータ（Ｐａ１，Ｐａ２）を再設定することによって第３の計測パラメータ（ＰＡ１，ＰＡ２）を取得するようにしたが、第１の計測パラメータと第２の計測パラメータとの関係を逆の関係としてもよい。すなわち、第２の医用画像上に第１の計測パラメータ（Ｐａ１，Ｐａ２）をプロットし、その第１の計測パラメータ（Ｐａ１，Ｐａ２）を用いて第２の計測パラメータ（Ｐｂ１，Ｐｂ２）を再設定することによって第３の計測パラメータ（ＰＢ１，ＰＢ２）を取得するようにしてもよい。第３の計測パラメータの取得方法としては、上記実施形態の取得方法を用いることができ、たとえば第２の医用画像上にプロットされた端点Ｐａ１と端点Ｐａ２とを通る直線上の点のうち、第２の医用画像における病変領域の輪郭上の点の座標を第３の計測パラメータ（ＰＢ１，ＰＢ２）として取得することができる。また、上述したように非剛体レジストレーションの結果を用いるようにしてもよい。

また、上記実施形態では、医用画像に含まれる病変の径を計測するようにしたが、上述したように計測対象は病変に限らず、医用画像に含まれる骨や関節でもよく、第１の医用画像内の骨や関節における２つの解剖学的な特徴点を第１の計測パラメータとして設定するとともに、第２の医用画像内の上記２つの解剖学的な特徴点を第２の計測パラメータとして設定し、上記実施形態と同様に、第１の計測パラメータと第２の計測パラメータとの間の変化が予め設定された変化量以上である場合に、第３の計測パラメータを設定したり、警告を行うようにしてもよい。

１ 医用画像計測装置
２ 医用画像保管サーバ
３ ディスプレイ
４ 入力装置
５ 医用画像
１０ 医用画像取得部
１１ 計測部
１２ 計測パラメータ取得部
１３ 評価値取得部
１４ 判定部
１５ 計測パラメータ再設定部
１６ 表示制御部
１７ 警告部

Claims (10)

1. 同一患者における同一の注目対象を異なる時点で撮影した第１の医用画像および第２の医用画像を取得する医用画像取得部と、
   前記第１の医用画像における前記注目対象の形状の特徴を計測するために設定された第１の計測パラメータおよび前記第２の医用画像における前記注目対象の形状の特徴を計測するために設定された第２の計測パラメータを取得する計測パラメータ取得部と、
   前記第１の計測パラメータと前記第２の計測パラメータとの間の変化を示す評価値を取得する評価値取得部と、
   前記評価値に基づいて、前記変化が予め設定された変化量以上であるか否かを判定する判定部と、
   前記判定部によって前記変化が予め設定された変化量以上であると判定された場合、前記第１の計測パラメータおよび前記第２の計測パラメータのうちのいずれか一方の計測パラメータを、他方の計測パラメータを用いて再設定することによって第３の計測パラメータを取得する計測パラメータ再設定部と、
   前記第３の計測パラメータを用いて、前記一方の計測パラメータが設定された前記第１の医用画像または前記第２の医用画像の前記注目対象の形状の特徴を計測する計測部とを備え、
   前記計測パラメータ再設定部が、前記第１の医用画像における前記注目対象の画像と前記第２の医用画像における前記注目対象の画像との非剛体レジストレーションの結果および前記他方の計測パラメータに基づいて、前記第３の計測パラメータを取得することを特徴とする医用画像計測装置。
2. 前記第１の計測パラメータ、前記第２の計測パラメータおよび前記第３の計測パラメータが、計測される前記注目対象の径を特定する情報である請求項１記載の医用画像計測装置。
3. 前記第１の計測パラメータ、前記第２の計測パラメータおよび前記第３の計測パラメータが、前記注目対象における２点間の距離を計測するために用いられる前記２点の座標である請求項１または２記載の医用画像計測装置。
4. 前記評価値取得部が、前記他方の計測パラメータである２点の座標を結ぶ直線と、前記一方の計測パラメータである２点の座標のそれぞれとの距離に基づいて前記評価値を取得する請求項１から３いずれか１項記載の医用画像計測装置。
5. 前記評価値取得部が、前記第１の計測パラメータである２点の座標を結ぶ直線と、前記第２の計測パラメータである２点の座標を結ぶ直線とがなす角度の情報を前記評価値として取得する請求項１から３いずれか１項記載の医用画像計測装置。
6. 前記第１の医用画像と前記第２の医用画像とが３次元画像であり、前記第１の計測パラメータ、前記第２の計測パラメータおよび第３の計測パラメータが、前記注目対象の異なる断面の径を特定する情報である請求項１から５いずれか１項記載の医用画像計測装置。
7. 前記注目対象が、時間の経過によって形状が変化する病変である請求項１から６いずれか１項記載の医用画像計測装置。
8. 前記判定部によって前記変化が予め設定された変化量以上であると判定された場合、警告を行う警告部を備えた請求項１から７いずれか１項記載の医用画像計測装置。
9. 同一患者における同一の注目対象を異なる時点で撮影した第１の医用画像および第２の医用画像を取得し、
   前記第１の医用画像における前記注目対象の形状の特徴を計測するために設定された第１の計測パラメータおよび前記第２の医用画像における前記注目対象の形状の特徴を計測するために設定された第２の計測パラメータを取得し、
   前記第１の計測パラメータと前記第２の計測パラメータとの間の変化を示す評価値を取得し、
   前記評価値に基づいて、前記変化が予め設定された変化量以上であるか否かを判定し、
   前記変化が予め設定された変化量以上であると判定された場合、前記第１の計測パラメータおよび前記第２の計測パラメータのうちのいずれか一方の計測パラメータを、他方の計測パラメータを用いて再設定することによって第３の計測パラメータを取得し、
   前記第３の計測パラメータを用いて、前記一方の計測パラメータが設定された前記第１の医用画像または前記第２の医用画像の前記注目対象の形状の特徴を計測する場合において、
   前記第１の医用画像における前記注目対象の画像と前記第２の医用画像における前記注目対象の画像との非剛体レジストレーションの結果および前記他方の計測パラメータに基づいて、前記第３の計測パラメータを取得することを特徴とする医用画像計測方法。
10. コンピュータを、同一患者における同一の注目対象を異なる時点で撮影した第１の医用画像および第２の医用画像を取得する医用画像取得部と、
    前記第１の医用画像における前記注目対象の形状の特徴を計測するために設定された第１の計測パラメータおよび前記第２の医用画像における前記注目対象の形状の特徴を計測するために設定された第２の計測パラメータを取得する計測パラメータ取得部と、
    前記第１の計測パラメータと前記第２の計測パラメータとの間の変化を示す評価値を取得する評価値取得部と、
    前記評価値に基づいて、前記変化が予め設定された変化量以上であるか否かを判定する判定部と、
    前記判定部によって前記変化が予め設定された変化量以上であると判定された場合、前記第１の計測パラメータおよび前記第２の計測パラメータのうちのいずれか一方の計測パラメータを、他方の計測パラメータを用いて再設定することによって第３の計測パラメータを取得する計測パラメータ再設定部と、
    前記第３の計測パラメータを用いて、前記一方の計測パラメータが設定された前記第１の医用画像または前記第２の医用画像の前記注目対象の形状の特徴を計測する計測部として機能させる医用画像計測プログラムであって、
    前記計測パラメータ再設定部が、前記第１の医用画像における前記注目対象の画像と前記第２の医用画像における前記注目対象の画像との非剛体レジストレーションの結果および前記他方の計測パラメータに基づいて、前記第３の計測パラメータを取得することを特徴とする医用画像計測プログラム。