(57)【要約】【課題】本発明は、肝臓の線維化の状態の変化を確認可能な画像処理装置および画像処理装置の作動方法並びに画像処理プログラムを提供する。【解決手段】被検者の肝臓を形態的変化が生じた２つの状態で撮影を行った第１のボリュームデータと第２のボリュームデータの肝臓領域間で非剛体位置合わせ処理を行うことにより、２つの肝臓領域内の対応する位置が変形した方向および距離を表す変形ベクトルを解析して取得した肝臓の硬さを表す指標を出力する。【選択図】図８
2019088672.tif 000002
【特許請求の範囲】
【請求項１】
被検者の肝臓を形態的変化が生じた２つの状態で撮影を行った第１のボリュームデータと第２のボリュームデータを取得する取得部と、
前記第１のボリュームデータの肝臓領域および前記第２のボリュームデータの肝臓領域間で非剛体位置合わせ処理を行うことにより、２つの肝臓領域内の対応する位置が変形した方向および距離を表す変形ベクトルを解析して取得した肝臓の硬さを表す指標を出力する指標出力部と、
を備えた画像処理装置。
【請求項２】
前記指標出力部が、前記指標を表す指標値として、前記変形ベクトルの大きさに基づいて取得した変形量を出力する請求項１記載の画像処理装置。
【請求項３】
前記変形量が、前記変形ベクトルの大きさの平均値である請求項２項記載の画像処理装置。
【請求項４】
前記第１のボリュームデータは、前記被検者が仰臥位の状態で撮影されたデータであり、
前記第２のボリュームデータは、前記被検者が伏臥位の状態で撮影されたデータである請求項１〜３のいずれか１項記載の画像処理装置。
【請求項５】
前記第１のボリュームデータは、前記被検者の肺野が拡大した状態で撮影されデータであり、
前記第２のボリュームデータは、前記被検者の肺野が前記第１のボリュームデータの肺野より縮小した状態で撮影されたデータである請求項１〜３のいずれか１項記載の画像処理装置。
【請求項６】
前記指標出力部が、前記第１のボリュームデータの肝臓領域の上端位置と前記第２のボリュームデータの肝臓領域の上端位置との体軸方向の差分値を用いて、前記指標が表す硬さのレベルを前記差分値が大きいほど前記指標が表す硬さのレベルを高くする調整、または前記差分値が小さくなるにしたがって前記指標が表す硬さのレベルを低くする調整を行う請求項５記載の画像処理装置。
【請求項７】
前記指標出力部が、前記指標を表す指標値として、前記変形ベクトルの大きさに基づいて取得した変形量を前記差分値で割った値を出力する請求項６記載の画像処理装置。
【請求項８】
前記指標出力部が、前記指標を予め用意したテーブルを用いて調整する請求項６または請求項７記載の画像処理装置。
【請求項９】
前記指標出力部が、前記第１のボリュームデータの肝臓領域および前記第２のボリュームデータの肝臓領域を構成する複数の肝臓の区域ごとに前記指標を取得する請求項１〜８のいずれか１項記載の画像処理装置。
【請求項１０】
前記指標出力部が、前記第１のボリュームデータの肝臓領域および前記第２のボリュームデータの肝臓領域の肝臓右葉上部の変形ベクトルを解析して前記指標を取得する請求項１〜８のいずれか１項記載の画像処理装置。
【請求項１１】
前記指標が肝線維化レベルであり、
前記指標出力部が、前記肝線維化レベルおよび患者の肝臓を形態的変化が生じた２つの状態で撮影を行った第１のボリュームデータと第２のボリュームデータとを有する教師データセットを用いて前記肝線維化レベルと前記変形ベクトルとの関係を学習させた判別器に、前記被検者の前記第１のボリュームデータと前記第２のボリュームデータを入力して前記肝線維化レベルを取得する請求項１記載の画像処理装置。
【請求項１２】
前記指標出力部が、前記第１のボリュームデータの肝臓領域および前記第２のボリュームデータの肝臓領域の各位置の前記変形ベクトルの大きさを前記指標として取得し、前記第１のボリュームデータの肝臓領域および前記第２のボリュームデータの肝臓領域の各位置の前記指標に応じた色で肝臓領域の各位置に対応する肝臓領域の外形を表す外形形状上の各位置を色分けして表示する請求項１記載の画像処理装置。
【請求項１３】
前記指標出力部が、前記第１のボリュームデータの肝臓領域および前記第２のボリュームデータの肝臓領域の各位置の前記変形ベクトルを前記指標として取得し、前記第１のボリュームデータの肝臓領域および前記第２のボリュームデータの肝臓領域の各位置の前記変形ベクトルを肝臓領域の各位置に対応する肝臓領域の外形を表す外形形状上の各位置に表示する請求項１記載の画像処理装置。
【請求項１４】
前記指標出力部は、前記非剛体位置合せの前に剛体位置合せを行う請求項１〜１３のいずれか１項記載の画像処理装置。
【請求項１５】
前記剛体位置合せは、前記第１のボリュームデータの肝臓領域の重心の位置と前記第２のボリュームデータの肝臓領域の重心の位置を一致させる請求項１４記載の画像処理装置。
【請求項１６】
前記剛体位置合せは、前記第１のボリュームデータの肝臓領域のランドマークと前記第２のボリュームデータの肝臓領域のランドマークを用いて行う請求項１４記載の画像処理装置。
【請求項１７】
前記第１のボリュームデータおよび前記第２のボリュームデータは、ＣＴ装置、ＭＲＩ装置、および超音波装置のいずれかの撮影により得られたデータである請求項１〜１６のいずれか１項記載の画像処理装置。
【請求項１８】
取得部と指標出力部とを備えた画像処理装置の作動方法であって、
前記取得部が、被検者の肝臓を形態的変化が生じた２つの状態で撮影を行った第１のボリュームデータと第２のボリュームデータを取得し、
前記指標出力部が、前記第１のボリュームデータの肝臓領域および前記第２のボリュームデータの肝臓領域間で非剛体位置合わせ処理を行うことにより、２つの肝臓領域内の対応する位置が変形した方向および距離を表す変形ベクトルを解析して取得した肝臓の硬さを表す指標を出力する画像処理装置の作動方法。
【請求項１９】
コンピュータを、
被検者の肝臓を形態的変化が生じた２つの状態で撮影を行った第１のボリュームデータと第２のボリュームデータを取得する取得部と、
前記第１のボリュームデータの肝臓領域および前記第２のボリュームデータの肝臓領域間で非剛体位置合わせ処理を行うことにより、２つの肝臓領域内の対応する位置が変形した方向および距離を表す変形ベクトルを解析して取得した肝臓の硬さを表す指標を出力する指標出力部として機能させる画像処理プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、医用画像の画像処理に関し、特に、肝臓を撮影した画像の解析に関する画像処理装置および画像処理装置の作動方法並びに画像処理装置プログラムに関する。
【背景技術】
【０００２】
線維化疾患は、全身のあらゆる臓器におこる疾患であり、様々な臓器において臓器不全を起こす。線維化とは、炎症などで受けた損傷を修復するときにできる「線維（コラーゲン）」というタンパク質が増加してひろがった状態をさすが、線維化が過剰になると、組織の伸縮ができなくなり、臓器が正常な機能を果たせなくなる。肝硬変は線維化疾患の１つである。
【０００３】
従来から、肝硬変の状態を確認するためには肝生検が行われていた。あるいは、肝硬変が進行すると肝表面の凹凸になるため、腹部ＣＴ（Ｃｏｍｐｕｔｅｄ  Ｔｏｍｏｇｒａｐｈｙ）で得られた画像を使って肝表面の凹凸の程度および／または肝臓の右葉と左葉の大きさのバランスなどを調べることによって肝硬変の診断が行われている。また、特許文献１では、ＭＲＩ（Ｍａｇｎｅｔｉｃ  Ｒｅｓｏｎａｎｃｅ  Ｉｍａｇｉｎｇ）装置を用いて肝臓が形態的に変化する前後（吸気と呼気など）の画像を撮影して、それぞれの画像から２枚の断面画像を取得して、２枚の断面画像中に表れた肝臓の辺縁を抽出し、辺縁の変化を表す指標から肝臓の硬さに関する情報を取得するようにしている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開平６−３１１９７８号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
実際に肝硬変の程度を知るには、肝生検を行うのが最も正確だが、肝生検では同じ位置のサンプリングを再現するはむずかしいため経時変化の観察には適さない。一方、特許文献１の手法では、断面画像の２次元画像の解析を行っているだけであり、断面の選び方によって結果がかわってくると考えられる。また、辺縁の情報だけ利用しているため、臓器内部の情報は活用されておらず正確な結果が得られていない可能性がある。このように、肝臓の硬さ、あるいは、肝臓の線維化の状態の変化を簡便かつ客観的に確認する有効な手段が確立されていない。
【０００６】
そこで、本発明では、上述のような問題を解決するために、肝臓の線維化の状態を簡便かつ客観的に確認することが可能な画像処理装置および画像処理装置の作動方法並びに画像処理装置プログラムを提供することを目的とする。
【課題を解決するための手段】
【０００７】
本発明の画像処理装置は、被検者の肝臓を形態的変化が生じた２つの状態で撮影を行った第１のボリュームデータと第２のボリュームデータを取得する取得部と、第１のボリュームデータの肝臓領域および第２のボリュームデータの肝臓領域間で非剛体位置合わせ処理を行うことにより、２つの肝臓領域内の対応する位置が変形した方向および距離を表す変形ベクトルを解析して取得した肝臓の硬さを表す指標を出力する指標出力部と、を備える。
【０００８】
本発明の画像処理装置の作動方法は、取得部と指標出力部とを備えた画像処理装置の作動方法であって、取得部が、被検者の肝臓を形態的変化が生じた２つの状態で撮影を行った第１のボリュームデータと第２のボリュームデータを取得し、指標出力部が、第１のボリュームデータの肝臓領域および第２のボリュームデータの肝臓領域間で非剛体位置合わせ処理を行うことにより、２つの肝臓領域内の対応する位置が変形した方向および距離を表す変形ベクトルを解析して取得した肝臓の硬さを表す指標を出力する。
【０００９】
本発明の画像処理プログラムは、コンピュータを、被検者の肝臓を形態的変化が生じた２つの状態で撮影を行った第１のボリュームデータと第２のボリュームデータを取得する取得部と、第１のボリュームデータの肝臓領域および第２のボリュームデータの肝臓領域間で非剛体位置合わせ処理を行うことにより、２つの肝臓領域内の対応する位置が変形した方向および距離を表す変形ベクトルを解析して取得した肝臓の硬さを表す指標を出力する指標出力部として機能させる。
【００１０】
また、指標出力部が、指標を表す指標値として、変形ベクトルの大きさに基づいて取得した変形量を出力することが望ましい。
【００１１】
また、変形量が、変形ベクトルの大きさの平均値であってもよい。
【００１２】
また、第１のボリュームデータは、被検者が仰臥位の状態で撮影されたデータであり、第２のボリュームデータは、被検者が伏臥位の状態で撮影されたデータであってもよい。
【００１３】
また、第１のボリュームデータは、被検者の肺野が拡大した状態で撮影されたデータであり、第２のボリュームデータは、被検者の肺野が第１のボリュームデータの肺野より縮小した状態で撮影されたデータであってもよい。
【００１４】
また、指標出力部が、第１のボリュームデータの肝臓領域の上端と第２のボリュームデータの肝臓領域の上端位置との体軸方向の差分値を用いて、指標が表す硬さのレベルを差分値が大きいほど指標が表す硬さのレベルを高くする調整、または差分値が小さくなるにしたがって指標が表す硬さのレベルを低くする調整を行うことが望ましい。
【００１５】
「硬さのレベルが高い」とは、「硬さのレベルが低い」場合より肝臓が硬い状態を指す。「指標が表す硬さのレベルを高くする調整」は、指標が現在の指標より肝臓がより硬い状態を指すように調整を行い、「指標が表す硬さのレベルを低くする調整」は、指標が現在の指標より肝臓がより柔らかい状態を指すように調整を行うことをいう。
【００１６】
また、指標出力部が、指標を表す指標値として、変形ベクトルの大きさに基づいて取得した変形量を差分値で割った値を出力してもよい。
【００１７】
また、指標出力部が、指標を予め用意したテーブルを用いて調整してもよい。
【００１８】
また、指標出力部が、第１のボリュームデータの肝臓領域および第２のボリュームデータの肝臓領域を構成する複数の肝臓の区域ごとに指標を取得してもよい。
【００１９】
また、指標出力部が、第１のボリュームデータの肝臓領域および第２のボリュームデータの肝臓領域の肝臓右葉上部の変形ベクトルを解析して指標を取得してもよい。
【００２０】
また、指標が肝線維化レベルであり、指標出力部が、肝線維化レベルおよび患者の肝臓を形態的変化が生じた２つの状態で撮影を行った第１のボリュームデータと第２のボリュームデータとを有する教師データセットを用いて肝線維化レベルと変形ベクトルとの関係を学習させた判別器に、被検者の第１のボリュームデータと第２のボリュームデータを入力して肝線維化レベルを取得してもよい。
【００２１】
また、指標出力部が、第１のボリュームデータの肝臓領域および第２のボリュームデータの肝臓領域の各位置の変形ベクトルの大きさを指標として取得し、第１のボリュームデータの肝臓領域および第２のボリュームデータの肝臓領域の各位置の指標に応じた色で肝臓領域の各位置に対応する肝臓領域の外形を表す外形形状上の各位置を色分けして表示してもよい。
【００２２】
また、指標出力部が、第１のボリュームデータの肝臓領域および第２のボリュームデータの肝臓領域の各位置の変形ベクトルを指標として取得し、第１のボリュームデータの肝臓領域および第２のボリュームデータの肝臓領域の各位置の変形ベクトルを肝臓領域の各位置に対応する肝臓領域の外形を表す外形形状上の各位置に表示してもよい。
【００２３】
また、指標出力部は、非剛体位置合せの前に剛体位置合せを行うのが望ましい。
【００２４】
また、剛体位置合せは、第１のボリュームデータの肝臓領域の重心の位置と第２のボリュームデータの肝臓領域の重心の位置を一致させることが望ましい。
【００２５】
また、剛体位置合せは、第１のボリュームデータの肝臓領域のランドマークと第２のボリュームデータの肝臓領域のランドマークを用いて行うようにしてもよい。
【００２６】
また、第１のボリュームデータおよび第２のボリュームデータは、ＣＴ装置、ＭＲＩ装置、および超音波装置のいずれかの撮影により得られたデータが望ましい。
【００２７】
本発明の他の画像処理装置は、コンピュータに実行させるための命令を記憶するメモリと、記憶された命令を実行するよう構成されたプロセッサとを備え、プロセッサは、被検者の肝臓を形態的変化が生じた２つの状態で撮影を行った第１のボリュームデータと第２のボリュームデータを取得し、第１のボリュームデータの肝臓領域および第２のボリュームデータの肝臓領域間で非剛体位置合わせ処理を行うことにより、２つの肝臓領域内の対応する位置が変形した方向および距離を表す変形ベクトルを取得し、変形ベクトルを解析して肝臓の硬さを表す指標を取得し、指標を出力する処理を実行する。
【発明の効果】
【００２８】
本発明によれば、被検者の肝臓を形態的変化が生じた２つの状態で撮影を行った第１のボリュームデータと第２のボリュームデータの肝臓領域間で非剛体位置合わせ処理を行うことにより、２つの肝臓領域内の対応する位置が変形した方向および距離を表す変形ベクトルを取得して肝臓の硬さを表す指標を出力するようにしたので、肝臓の線維化の状態の変化を簡便かつ客観的に確認することが可能になる。
【図面の簡単な説明】
【００２９】
【図１】医療情報システムの概略構成を表す図
【図２】本発明の第１の実施の形態の画像処理装置の概略構成を表す図
【図３Ａ】被検者が仰臥位の状態と伏臥位の状態で肝臓を撮影した正常例の画像
【図３Ｂ】被検者が仰臥位の状態と伏臥位の状態で肝臓を撮影した非正常例の画像
【図４Ａ】被検者が呼気の状態と吸気の状態で肝臓を撮影した正常例の画像
【図４Ｂ】被検者が呼気の状態と吸気の状態で肝臓を撮影した非正常例の画像
【図５】非剛体位置合せを説明するための図
【図６】変形ベクトルの大きさの分布の一例
【図７】肝臓領域を構成する複数の区域を説明するための図
【図８】肝臓領域内部の各位置の変形ベクトルを矢印で表示した例
【図９】肝臓領域内の複数の位置の変形ベクトルの大きさを色分け表示した例
【図１０】第１の実施の形態の画像処理装置の処理の流れを示すフローチャート
【図１１】本発明の第２の実施の形態の画像処理装置の概略構成を表す図
【図１２】第２の実施の形態の画像処理装置の処理の流れを示すフローチャート
【発明を実施するための形態】
【００３０】
以下、図面を参照して本発明の第１の実施の形態の画像処理装置を備えた医療情報システムについて説明する。図１は、本実施の形態の医療情報システムの概略構成を示すブロック図である。
【００３１】
本実施の形態の医療情報システムは、図１に示すように、画像処理装置１、医用画像保管サーバ２、および、撮影装置３（以下、モダリティという）がネットワーク４を介して互いに通信可能な状態で接続されて構成されている。
【００３２】
モダリティ３は、たとえばＣＴ装置、ＭＲＩ装置および超音波装置などであり、撮影された３次元のボリュームデータは、ＤＩＣＯＭ（Ｄｉｇｉｔａｌ  Ｉｍａｇｉｎｇ  ａｎｄ  ＣＯｍｍｕｎｉｃａｔｉｏｎｉｎ  Ｍｅｄｉｃｉｎｅ）規格に準拠した格納フォーマットおよび通信規格に従って、ネットワーク４を介して医用画像保管サーバ２に送信されて格納される。
【００３３】
画像処理装置１は、汎用のコンピュータであり、ＣＰＵ（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）、メモリ（主記憶装置）、ストレージ（補助記憶装置）、入出力インターフェース、通信インターフェース、入力装置、表示装置、および、データバスなどの周知のハードウェア構成を備え、周知のオペレーションシステムなどがインストールされている。また、表示装置として液晶ディスプレイなどを有し、入力装置としてキーボードおよび／またはマウスなどのポインティングデバイスを有している。ストレージは、ハードディスクまたはＳＳＤ(Ｓｏｌｉｄ  Ｓｔａｔｅ  Ｄｒｉｖｅ)などで構成される。なお、必要に応じてコンピュータにＧＰＵ（Ｇｒａｐｈｉｃｓ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）を設けるようにしてもよい。このコンピュータに、本実施の形態の画像処理プログラムをインストールすることにより画像処理装置１として機能する。また、画像処理装置１は、医用画像保管サーバ２に対する画像の送信要求、および、医用画像保管サーバ２から画像の受信を行う機能を備え、各機能のためのソフトウェアプログラムを実行することにより行われる。
【００３４】
画像処理プログラムは、ＤＶＤ（Ｄｉｇｉｔａｌ  Ｖｅｒｓａｔｉｌｅ  Ｄｉｓｃ）またはＣＤ−ＲＯＭ（Ｃｏｍｐａｃｔ  Ｄｉｓｃ  Ｒｅａｄ  Ｏｎｌｙ  Ｍｅｍｏｒｙ）などの記録媒体に記録されて配布され、その記録媒体からコンピュータにインストールされる。または、画像処理プログラムは、ネットワークに接続されたサーバコンピュータの記憶装置もしくはネットワークストレージに対して、外部からアクセス可能な状態で記憶され、外部からの要求に応じてコンピュータにダウンロードされた後に、インストールされるようにしてもよい。
【００３５】
第１の実施の形態の画像処理装置１は、図２に示すように、取得部１０、肝臓領域抽出部１１、グローバル位置合せ部１２、ローカル位置合せ部１３、指標取得部１４、および指標出力部１５を備える。
【００３６】
取得部１０は、被検者の肝臓を形態的変化が生じた２つの状態で撮影を行った第１のボリュームデータＶ１と第２のボリュームデータＶ２を医用画像保管サーバ２から取得する。ボリュームデータＶ１，Ｖ２は、本実施の形態においては、ＣＴ装置、ＭＲＩ装置または超音波撮影装置などによって撮影された３次元の画像データである。ボリュームデータは、３次元空間を細かく区切ったボクセル位置の物理量を表わすデータで構成され、データは各ボクセルの位置に存在する臓器または組織などを放射線または磁気が透過した透過量に応じた値、あるいは、臓器または組織などで反射された超音波に応じた値で表される。
【００３７】
第１のボリュームデータＶ１と第２のボリュームデータＶ２は、肝臓の形態的変化が生じる前後の２つの状態で撮影を行った画像である。具体的には、肝臓の形態的変化は、体位の変化または呼吸の影響などによって生じる。
【００３８】
例えば、被検者が、仰向け（仰臥位）になっているか、うつ伏せ（伏臥位）になっているかによって肝臓に形態的変化が生じる。そこで、被検者が仰臥位の状態で撮影された画像を第１のボリュームデータＶ１とし、被検者が伏臥位の状態で撮影された画像を第２のボリュームデータＶ２とする。図３Ａは、正常な肝臓を被検者が仰臥位の状態と伏臥位の状態で撮影した画像である。左側のａの下にある画像が仰臥位の状態で撮影したアキシャル画像、サジタル画像、およびコロナル画像である。右側のｂの下にある画像が伏臥位の状態で撮影したアキシャル画像、サジタル画像、およびコロナル画像である。図３Ｂは、肝硬変が進んだ肝臓を被検者が仰臥位の状態（左側のａ）と伏臥位の状態（右側のｂ）で撮影した画像である。図３Ａの正常例の画像の肝臓の上部の形状と下部の形状（黒い線）が大きく変化しているが、図３Ｂの非正常例の画像の肝臓の上部の形状（黒い線）の変化が少なく、２つを比較すると、非正常例は正常例に比べて形状変化が少ない。
【００３９】
あるいは、呼吸をすることで肺野の大きさが変わることによっても肝臓に形態的変化が生じる。そこで、被検者の肺野が拡大した状態で撮影した画像を第１のボリュームデータＶ１とし、被検者の肺野が第１のボリュームデータＶ１の肺野より縮小した状態で撮影された画像を第２のボリュームデータＶ２としてもよい。第１のボリュームデータＶ１と第２のボリュームデータＶ２は肺野の大きさの変化が大きい方が好ましく、肺野の大きさが最大になる最大吸気の状態と、肺野の大きさが最小になる最小呼気の状態で撮影するのが好ましい。図４Ａは、正常な肝臓を被検者が吸気の状態と呼気の状態で撮影した画像である。左側のａの下にある画像が吸気の状態で撮影したアキシャル画像、サジタル画像、およびコロナル画像である。右側のｂの下にある画像が呼気の状態で撮影したアキシャル画像、サジタル画像、およびコロナル画像である。図４Ｂは、肝硬変が進んだ肝臓を被検者が吸気の状態（左側のａ）と呼気の状態（右側のｂ）で撮影した画像である。図４Ａの正常例の画像の肝臓の上部の形状（白い矢印）と、図４Ｂの非正常例の画像の肝臓の上部の形状（白い矢印）を比較すると、非正常例は正常例に比べて形状変化が少ない。
【００４０】
第１のボリュームデータＶ１と第２のボリュームデータＶ２は、肺野の大きさが大きく変わる最大吸気と最小呼気の状態で撮影するのが好ましいが、肺野の大きさが異なっていればよい。例えば、造影剤を投与して複数のフェーズで撮影を行った時の画像を利用することができる。通常、造影剤を用いてＣＴ撮影を行なう際には造影剤が拡散する状態に応じて４相のフェーズで撮影が行われる。まず、息を吸った状態で呼吸停止を行い２相（早期動脈相、後期動脈相）の撮影を行い、その後息を吐いた状態で残りの２相の撮影を行うようにすることで、肺野の大きさが変化した状態でボリュームデータを撮影することができる。この４相のボリュームデータ中から肺野の大きさが大きく変わったボリュームデータを選んで、第１のボリュームデータＶ１と第２のボリュームデータＶ２としてもよい。
【００４１】
ボリュームデータＶ１，Ｖ２は、医用画像保管サーバ２に患者の識別情報とともに予め保管されている。取得部１０は、キーボードなどの入力装置を用いてユーザによって入力された患者の識別情報に基づいて、診断対象の患者の第１のボリュームデータＶ１と第２のボリュームデータＶ２を医用画像保管サーバ２から読み出してストレージ（不図示）に記憶する。
【００４２】
肝臓領域抽出部１１は、第１のボリュームデータＶ１と第２のボリュームデータＶ２のそれぞれから肝臓領域を抽出する。肝臓領域は、撮影されたモダリティに応じた画素値などに応じて抽出され、既知の様々な手法を用いることができる。
【００４３】
グローバル位置合せ部１２は、剛体位置合せによって肝臓の大域的な位置合せを行う。剛体位置合せでは、肝臓の形を変えないように平行移動と回転移動だけを用いた位置合わせが行われる。具体的には、第１のボリュームデータＶ１の肝臓領域の重心の位置と第２のボリュームデータＶ２の肝臓領域の重心の位置を一致させるように、Ｘ，Ｙ，およびＺのそれぞれの軸方向に平行移動を行う。さらに、剛体変換モデルのＸ，Ｙ，およびＺのそれぞれの軸方向の平行移動と回転移動の６つのパラメータを最適化して、２つの肝臓領域の重なる体積が最も大きくなるように位置合わせを行うようにする。あるいは、第１のボリュームデータＶ１の肝臓領域のランドマークに対応する第２のボリュームデータＶ２の肝臓領域のランドマークが近づくように、平行移動および／または回転移動を行う。ランドマークは、肝臓を代表する位置を用いればよく、肝臓領域の上端、右葉と左葉の境界、門脈、または肝動脈などの位置を用いることができ、これらのランドマークのうち１つ以上のランドマークを用いて、２つの肝臓のランドマーク間の距離の合計が最も小さくなるように位置合せを行う。
【００４４】
次に、ローカル位置合せ部１３は、グローバル位置合せ部１２で剛体位置合わせが行われた肝臓領域に対して、非剛体位置合せによって肝臓領域内の詳細な位置合せを行う。非剛体位置合せでは、位置合わせする２つの肝臓のうち少なくとも一方の形を変えて位置合せが行われる。２つの異なるボリュームデータ間で位置合わせをするための非剛体位置合せとして、アフィン変換などの線形位置合せと、Ｂ−ｓｐｌｉｎｅを用いた非線形位置合せが知られている。このような線形位置合わせと非線形位置合わせを用いて医用画像の位置合せが行われる。
【００４５】
医用画像の位置合せは、第１のボリュームデータＶ１の画像ＩＦ内の各点ｘＦは固定したままにして、第２のボリュームデータＶ２の画像ＩＭ内の各点ｘＭを移動させることで、画像ＩＭの各点ｘＭの位置を解剖学的に対応する画像ＩＦの各点ｘＦの位置に一致させるための幾何学的変換の関数Ｔを求める問題である。この問題は、図５に示すように、一方の画像ＩＦと、他方の画像ＩＭをＴによって変形した画像との一致度を評価関数として、これを最大化させる最適化問題として解くことができる。大域的な最適解を得る方法はないので、ある初期値を設定した後に、勾配法などを使って反復的によりよい解を探索することによって解を求める。一致度には、相互情報量を用いることができる。相互情報量はＣＴ画像とＭＲＩ画像のように異なるモダリティで撮影された画像の位置合わせに用いられるが、同じモダリティで撮影した画像間の場合には、相互相関関数を用いるようにしてもよい。
【００４６】
具体的には、まず、幾何学的変換のモデルとして、Ｂ−ｓｐｌｉｎｅ変換モデルを用いてボリュームデータの中に格子点を設け、格子点を移動させて一致する点を捜すことで粗い位置合わせを行う。その後、Ｄｉｆｆｅｃｏｍｏｒｐｈｉｃ（微分同相）変換モデルで画素単位で詳細な位置合わせ行う。つまり、まずＢ−ｓｐｌｉｎｅ変換モデルの最適解を求めた後、その結果を初期値としてＤｉｆｆｅｃｏｍｏｒｐｈｉｃ変換モデルの最適解を得る。これにより、剛体位置合わせ後の第１のボリュームデータと第２のボリュームデータの２つの肝臓領域間において、２つの肝臓領域内の対応する位置を見つけることができる。
【００４７】
指標取得部１４は、ローカル位置合せ部１３の非剛体位置合せで得られた２つの肝臓領域内の対応する位置から、２つの肝臓領域が変形した方向および距離を表す変形ベクトルを取得する。つまり、非剛体位置合わせによって見つけられた画像ＩＦの肝臓領域内の各点ｘＦに対応する画像ＩＭの肝臓領域内の各点ｘＭを結ぶベクトルを変形ベクトルとして取得する（なお、画像ＩＦと画像ＩＭは肝臓領域が剛体位置合わせされた後の画像である）。さらに、変形ベクトル場から肝臓の硬さを表す指標を取得する。例えば、この変形ベクトルの変形量を指標値として算出する。
【００４８】
図６に、変形ベクトルの大きさの分布の例を示す。図６は横軸が変形ベクトルの大きさ（変形ベクトル長）を示し、縦軸が相対頻度を示す。また、実線は、移植が必要なくらい線維化が進んだ状態の肝臓の変形ベクトルの分布の例を示し、破線は、線維化はしているが重篤ではない状態の肝臓の変形ベクトルの分布の例を示し、一点鎖線は、健康な状態の肝臓の変形ベクトルの分布の例を示す。このような変形ベクトルの大きさの平均値を求めて指標値とすることができる。変形ベクトルの大きさの平均値は、値が小さいほど線維化が進んだ状態を示す。また、変形ベクトルの大きさの平均値と標準偏差を、肝臓の硬さを表す指標値として算出してもよい。あるいは、変形ベクトルの大きさの分布と肝線維化レベルとの関係を多数のデータから統計的に解析して、線維化ステージを表す肝線維化レベルを指標としてもよい。肝線維化レベルは線維化ステージに対応させるのが好ましい。
【００４９】
また、指標は、数値でなくてもよく、肝臓の硬さのレベルを文字で表してもよい。例えば、変形ベクトルの大きさの平均値または肝線維化レベルに応じて、肝硬変の所見がない健康な肝臓の場合には「Ａ」、軽度の肝硬変の場合には「Ｂ」、中程度の肝硬変の場合には「Ｃ」、重度の肝硬変の場合には「Ｄ」というように、肝臓の硬さのレベルを文字で表わした指標を用いてもよい。さらに、指標は、文字に限らず、数字であってもよいし、文字と数字の組み合わせであってもよい。さらに、指標が「中程度の肝硬変」というように文章であってもよい。
【００５０】
あるいは、指標取得部１４が、肝臓領域を構成する複数の区域ごとに指標を取得するようにしてもよい。肝臓は、肝鎌状間膜を境に左葉と右葉に分かれ、さらに、外側区域、内側区域、前区域、および後区域の４区域に分かれる。これらの区域をさらに２区域ずつに分けて、Ｓ１の尾状葉区域、Ｓ２の後外側区域（外側上区域）、Ｓ３の前外側区域（外側下区域）、Ｓ４の内側上下区域（方形葉）、Ｓ５の前下区域、Ｓ６の後下区域、Ｓ７の後上区域、およびＳ８の前上区域の８区域に分けられる（図７参照）。これらのそれぞれの区域ごとに指標を求めるようにしてもよい。
【００５１】
また、肝臓の右葉上部は肝臓の中でも薄く平べったい形状をしているため、被検者の体位や呼吸の状態によって形態的変化が大きく表れる。そこで、肝臓領域の肝臓右葉上部の変形ベクトルを解析して指標を取得することが好ましい。
【００５２】
さらに、指標取得部１４は、判別器を用いて肝線維化レベルを指標値として取得するようにしてもよい。例えば、患者の肝臓を形態的変化が生じた２つの状態で撮影を行った第１のボリュームデータＶ１と第２のボリュームデータＶ２に対して肝線維化レベルが対応付けられたデータを教師データセットとして多数用意して、肝線維化レベルと変形ベクトルとの関係を機械学習させた判別器をあらかじめ生成しておき、判別器を用いて、被検者の第１のボリュームデータＶ１と第２のボリュームデータＶ２を入力して肝線維化レベルを取得するようにしてもよい。
【００５３】
指標出力部１５は、液晶ディスプレイなどの表示装置に指標値または指標を表す文字を出力して表示する。具体的には、図３Ａ，図３Ｂ，図４Ａ，図４Ｂのようにボリュームデータを表示した画面上に硬さのレベルを表す数値や文字を表示する。あるいは、指標をプリンターに出力して、紙などに印刷するようにしてもよい。
【００５４】
あるいは、指標取得部１４が、指標として肝臓領域内の複数の位置の変形ベクトル自体を取得してもよい。この場合には、指標出力部１５は、肝臓領域の各位置の変形ベクトルを肝臓領域の各位置に対応する肝臓領域の外形を表す外形形状上の各位置に表示するようにしてもよい。図８は、肝臓のあるアキシャル断面における肝臓領域内部の各位置の変形ベクトルを表示した例である。また、図８では、Ｚ方向（体軸方向）の変位を黒矢印が＋方向、白抜きの矢印が−方向を表している。また、図８は、便宜上、Ｚ方向の変位を黒の矢印と白抜きの矢印で表しているが、赤と青など異なる色でディスプレイ上に表示するようにしてもよい。
【００５５】
さらに、指標取得部１４が、指標として肝臓領域内の複数の位置の変形ベクトルから得られた大きさを取得してもよい。この場合には、指標出力部１５は、肝臓領域の各位置の指標に応じた色で、肝臓領域の外形を表す外形形状上の各位置を色分けして表示するようにしてもよい。図９は、便宜上、色分けをハッチングで表した例である。図９の異なるハッチングがそれぞれ赤と青の色分けに対応し、赤の場所は変位ベクトルの大きさが所定の値以上であり、青の場所は変位ベクトルの大きさが所定の値以下であることを示すようにしてもよい。
【００５６】
次に本実施の形態の画像処理装置１の処理の流れについて、図１０に示すフローチャートを参照しながら説明する。
【００５７】
まず、ユーザによって患者の識別情報が入力されると、その患者の造影ＣＴ画像から２つのボリュームデータＶ１，Ｖ２を取得部１０で取得する（ＳＴ１０）。
【００５８】
次に、取得部１０で取得された２つのボリュームデータＶ１，Ｖ２から肝臓領域抽出部１１を用いて肝臓領域を抽出する（ＳＴ１１）。グローバル位置合せ部１２で、肝臓の重心が一致するように剛体変換のＸ，Ｙ，およびＺのそれぞれの軸方向に平行移動させる。さらに、肝臓の重心を一致させる平行移動のパラメータを初期値として、Ｘ，Ｙ，およびＺのそれぞれの軸方向の平行移動と回転移動の６つのパラメータを最適化して剛体位置合わせを行う（ＳＴ１２）。
【００５９】
さらに、ローカル位置合せ部１３は、グローバル位置合せ部１２で剛体位置合わせが行われた肝臓領域に対して非剛体位置合せを行うことで肝臓領域内の詳細な位置合せを行ない、第１のボリュームデータＶ１と第２のボリュームデータＶ２の肝臓領域内の対応する位置を見つける（ＳＴ１３）。
【００６０】
指標取得部１４は、非剛体位置合せの結果に基づいて、第１のボリュームデータＶ１と第２のボリュームデータＶ２の２つの肝臓領域が変形した方向および距離を表す変形ベクトルを取得する（ＳＴ１４）。さらに、変形ベクトルから指標を取得する（ＳＴ１５）。指標出力部１５は、指標をディスプレイ上に表示する（ＳＴ１６）。
【００６１】
以上詳細に説明したように、非剛体位置合わせにより肝臓の変形ベクトルを取得して、肝臓の硬さのレベルを取得するようにしたので、肝臓の線維化の状態を簡便かつ客観的に確認することが可能になる。
【００６２】
上述の第１の実施の形態では、非剛体位置合わせの前に剛体位置合わせを行う場合について説明したが、例えば、仰臥位の状態と伏臥位の状態で撮影する場合のように、肝臓の体軸方向の位置が大きく変わらないで２つのボリュームデータの撮影が行われる場合には、剛体位置合わせを行うことなく、非剛体位置合わせのみ行うようにしてもよい。
【００６３】
次に、第２の実施の形態について説明する。本実施の形態では、肺野の大きさに応じて肝臓の形態的変化が生じる場合に、より正確に硬さのレベルを取得する手法について説明する。第１の実施の形態と同じ構成については、第１の実施の形態と同一符号を付して、詳細な説明は省略し、異なる構成についてのみ詳細な説明を行う。
【００６４】
第２の実施の形態の医療情報システムは、第１の実施の形態と同様であるので詳細な説明は省略する。第２の実施の形態の画像処理装置１ａは、図１１に示すように、取得部１０、肝臓領域抽出部１１、グローバル位置合せ部１２、ローカル位置合せ部１３、指標取得部１４ａ、および指標出力部１５を備える。また、本実施の形態の指標取得部１４ａは、指標の調整を行う調整部１６を備える点で第１の実施の形態とは異なる。
【００６５】
指標取得部１４ａは、第１の実施の形態と同様に、ローカル位置合せ部１３によって行われた非剛体位置合せの結果に基づいて、第１のボリュームデータＶ１と第２のボリュームデータＶ２の２つの肝臓領域が変形した方向および距離を表す変形ベクトルを取得し、変形ベクトルに基づいて指標を求める。
【００６６】
しかし、同じ被検者でも、呼気や吸気の状態が異なると結果も異なってしまう。例えば、２つのボリュームデータＶ１，Ｖ２を撮影した時の呼吸の状態が最大吸気と最小呼気であれば、肝硬変が進行していないときは肺野の大きさが大きく変わるため肝臓領域の変形量は大きな値になるが、２つのボリュームデータＶ１，Ｖ２を撮影した時の肺野の大きさの変化が小さいときには、肝硬変が進行していなくても肝臓領域の変形量は比較的小さい値になる。一方、肝硬変が進んでいる場合には、肺野の大きさが大きく変わっていても、肺野の大きさがあまり変わらなくても変形量は小さい値となる。
【００６７】
呼気と吸気では、横隔膜の位置が大きく変わり、それに伴って肝臓の位置も上下に移動する。呼吸の状態を変えた２つのボリュームデータＶ１，Ｖ２を撮影するときは、撮影装置に対する被検者の体位の移動はないため、肺野の大きさの変化を、肝臓領域の体軸方向の位置の変化から推定することが可能である。
【００６８】
そこで、調整部１６は、グローバル位置合せ部１２による位置合せが行われる前の第１のボリュームデータＶ１の肝臓領域および第２のボリュームデータＶ２の肝臓領域の上端位置の体軸方向の差分値を取得し、指標が表す硬さのレベルを差分値が大きいほど指標が表す硬さのレベルを高くする調整を行い、差分値が小さくなるにしたがって指標が表す硬さのレベルを低くする調整を行う。つまり、肺野の大きさの大きく変わる場合には、肺野の大きさがあまり変わらない場合よりも、変異ベクトルの大きさが大きくなり指標は肝臓が柔らかいこと（硬さのレベルが低い）を指すので、硬さのレベルを高くする調整を行う。一方、肺野の大きさがあまり変わらない場合には、肺野の大きさが大きく変わる場合よりも、変異ベクトルの大きさが小さくなり指標は実際より肝臓が硬いこと（硬さのレベルが高い）を指すので、硬さのレベルを低くする調整を行う。これによって、呼気や吸気の状態によって変わる肺野の大きさの差異を吸収するような補正を加えることができる。
【００６９】
具体的には、指標取得部１４ａが、指標を表す指標値として求めた変形ベクトルの大きさに基づいて取得した変形量を差分値で割った値を指標値とする。あるいは、変形量を差分値で割った値に応じて、硬さのレベルを「Ａ」〜「Ｄ」文字で表した指標で表してもよい。
【００７０】
指標取得部１４ａの調整部１６が、指標を予め変形ベクトルの変形量と差分値に応じた指標（指標値、または、硬さのレベルを文字で表した指標）を表すテーブルを用意して、指標を調整するようにしてもよい。あるいは、変形ベクトルから得られた指標と差分値に応じた指標値を表すテーブルを用意して、指標を調整するようにしてもよい。
【００７１】
次に本実施の形態の画像処理装置１の処理の流れについて、図１２に示すフローチャートを参照しながら説明する。
【００７２】
まず、ユーザが患者の識別情報を入力すると、その患者の造影ＣＴ画像から２つのボリュームデータＶ１，Ｖ２を取得部１０で取得する（ＳＴ２０）。
【００７３】
次に、取得部１０で取得された２つのボリュームデータＶ１，Ｖ２から肝臓領域抽出部１１を用いて肝臓領域を抽出する（ＳＴ２１）。ここで、調整部１６は、２つのボリュームデータＶ１，Ｖ２の肝臓領域の上端位置の体軸方向の差分値を取得する（ＳＴ２２）。
【００７４】
次に、グローバル位置合せ部１２で、肝臓の重心が一致するように剛体変換のＸ，Ｙ，およびＺのそれぞれの軸方向に平行移動させる。さらに、肝臓の重心を一致させる平行移動のパラメータを初期値として、Ｘ，Ｙ，およびＺのそれぞれの軸方向の平行移動と回転移動の６つのパラメータを最適化して剛体位置合わせを行う（ＳＴ２３）。
【００７５】
さらに、ローカル位置合せ部１３は、グローバル位置合せ部１２で剛体位置合わせが行われた肝臓領域に対して非剛体位置合せを行うことで肝臓領域内の詳細な位置合せを行ない、第１のボリュームデータＶ１と第２のボリュームデータＶ２の肝臓領域内の対応する位置を見つける（ＳＴ２４）。
【００７６】
指標取得部１４ａは、非剛体位置合せの結果に基づいて、第１のボリュームデータＶ１と第２のボリュームデータＶ２の２つの肝臓領域が変形した方向および距離を表す変形ベクトルを取得する（ＳＴ２５）。さらに、変形ベクトルから指標を取得する（ＳＴ２６）。ここで、調整部１６で指標が表す硬さのレベルを差分値が大きいほど指標が表す硬さのレベルを高くし、差分値が小さくなるにしたがって指標が表す硬さのレベルを低くするように調整を行う（ＳＴ２７）。指標出力部１５は、調整された指標をディスプレイ上に表示する（ＳＴ２８）。
【００７７】
以上詳細に説明したように、非剛体位置合わせにより肝臓の変形ベクトルを取得して、肝臓の硬さのレベルを取得し、さらに、呼吸の状態によって肝臓の変形の程度が異なってくることを考慮することで、肝臓の線維化の状態をより正確に取得することが可能になる。
【００７８】
上記では、１つのコンピュータを画像処理装置として機能させる場合について説明したが、複数のコンピュータに機能を分散させるようにしてもよい。
【００７９】
また、上記では、汎用コンピュータが画像処理装置として機能する場合について説明したが、専用コンピュータによって実施されてもよい。専用コンピュータは、内蔵されたＲＯＭ(Ｒｅａｄ  Ｏｎｌｙ  Ｍｅｍｏｒｙ)やフラッシュメモリなど、不揮発メモリに記録されたプログラムを実行するファームウェアであってもよい。さらに、この画像処理装置の少なくとも一部の機能を実行するためのプログラムを永久的に記憶するＡＳＩＣ（Ａｐｐｌｉｃａｔｉｏｎ  Ｓｐｅｃｉｆｉｃ  Ｉｎｔｅｇｒａｔｅｄ  Ｃｉｒｃｕｉｔ：特定用途向け集積回路）やＦＰＧＡ（ｆｉｅｌｄ  ｐｒｏｇｒａｍｍａｂｌｅ  ｇａｔｅ  ａｒｒａｙｓ）などの専用回路を設けるようにしてもよい。あるいは、専用回路に記憶されたプログラム命令と、専用回路のプログラムを利用するようにプログラムされた汎用のＣＰＵによって実行されるプログラム命令と組み合わせるようにしてもよい。以上のように、コンピュータのハードウェア構成をどのように組み合わせてプログラム命令を実行してもよい。
【符号の説明】
【００８０】
１、１ａ      画像処理装置
２     医用画像保管サーバ
３     撮影装置
４     ネットワーク
１０   取得部
１１   肝臓領域抽出部
１２   グローバル位置合せ部
１３   ローカル位置合せ部
１４、１４ａ  指標取得部
１５   指標出力部
１６   調整部
Ｖ１   第１のボリュームデータ
Ｖ２   第２のボリュームデータ
【図１】
2019088672.tif 000003
【図２】
2019088672.tif 000004
【図３Ａ】
2019088672.tif 000005
【図３Ｂ】
2019088672.tif 000006
【図４Ａ】
2019088672.tif 000007
【図４Ｂ】
2019088672.tif 000008
【図５】
2019088672.tif 000009
【図６】
2019088672.tif 000010
【図７】
2019088672.tif 000011
【図８】
2019088672.tif 000012
【図９】
2019088672.tif 000013
【図１０】
2019088672.tif 000014
【図１１】
2019088672.tif 000015
【図１２】
2019088672.tif 000016