JPWO2020075571A1 - 流体解析装置、方法およびプログラム - Google Patents

Abstract

第１解析部が、第１の３次元画像を解析して、血管内の各位置における血液の流れに関する第１の結果を導出する。第２解析部が、血管内の各位置における血液の流れをシミュレーションして、血液の流れに関する第２の結果を導出する。一致度導出部が、第１の結果および第２の結果のいずれか一方を基準とした、第１の結果と第２の結果との対応する位置における一致度を導出する。表示制御部が、一致度に応じて第１の結果および第２の結果の少なくとも一方を表示部に表示する。

Description

本開示は、流体の流れを解析して表示する流体解析装置、方法およびプログラムに関する。

近年、心臓および脳等を撮影した医用画像を用いて、血管内の血流を解析することが行われている。このような医用画像を用いた血流解析方法としては、例えば、実際の血流を４次元的に測定する４Ｄフローの手法が用いられている。４Ｄフローは、例えば３次元シネ位相コントラスト磁気共鳴法によって撮影された３次元のＭＲＩ（Magnetic Resonance Imaging）画像を用いて、ボクセル毎、ピクセル毎または領域毎に流速ベクトルを導出し、これを時間の流れと合わせて動的に表示する手法である。また、数値流体力学を用いた血流解析（ＣＦＤ：Computational Fluid Dynamics）により、血液の流れをシミュレーションにより把握する手法も用いられている。

ところで、４Ｄフローは実際の血流が表現されたＭＲＩ画像に基づいて導出されるが、空間分解能および時間分解能に制約がある。また、乱流等が起こっている部分は血流が過小評価されるため、その部分においてはＭＲＩ画像に速度情報が現れてこない場合がある。一方、ＣＦＤは空間分解能および時間分解能を任意に設定可能であるが、あくまでもシミュレーションであるため、その結果が妥当であるかを十分に見極める必要がある。

このため、４Ｄフローにより導出された結果である流速ベクトルと、ＣＦＤにより導出された結果である流速ベクトルとを組み合わせて、血流を分析する手法が提案されている（特開２０１６−１３５２６５号公報参照）。また、４Ｄフローにより導出された結果と、ＣＦＤにより導出された結果とを表示して、ユーザに確認させるための手法も提案されている（van Ooij P, et al., "3D cine phase-contrast MRI at 3T in intracranial aneurysms compared with patient-specific computational fluid dynamics", AJNR Am J Neuroradiol. 2013 Sep;34(9):1785-91）。特開２０１６−１３５２６５号公報およびOoijPらの文献に記載された手法によれば、４Ｄフローにより導出された結果と、ＣＦＤにより取得された結果との双方を用いて、両者の長所を生かしつつ、診断および治療を行うことが可能となる。

しかしながら、OoijPらの文献に記載された手法は、４Ｄフローのように画像を解析することにより導出された結果と、ＣＦＤのようにシミュレーションにより導出された結果とが別々に表示される。このため、画像を解析することにより導出された結果とシミュレーションにより導出された結果とを同時に確認することができず、両者の相違が分かりにくい。

本開示は上記事情に鑑みなされたものであり、画像を解析することにより導出された結果とシミュレーションにより導出された結果との相違を容易に認識できるようにすることを目的とする。

本開示による流体解析装置は、内部に流体が流れる構造物を含む被写体を撮影することにより取得した画像を解析して、構造物内の各位置における流体の流れに関する第１の結果を導出する第１解析部と、
構造物内の各位置における流体の流れをシミュレーションして、流体の流れに関する第２の結果を導出する第２解析部と、

第１の結果および第２の結果のいずれか一方を基準とした、第１の結果と第２の結果との対応する位置における一致度を導出する一致度導出部と、
一致度に応じて第１の結果および第２の結果の少なくとも一方を表示部に表示する表示制御部とを備える。

「各位置」とは画像のピクセル位置およびボクセル位置のみならず、複数ピクセル（ボクセル）からなる領域の位置をも含む。

「流体の流れに関する結果」としては、流速ベクトル、流体の圧力、壁面せん断応力（ＷＳＳ(Wall Shear Stress)）、渦度およびヘリシティ等が挙げられる。

なお、本開示による流体解析装置においては、表示制御部は、第１の結果および第２の結果のいずれか一方を表示し、一致度に応じて、表示された第１の結果および第２の結果のいずれか一方の表示態様を変更するものであってもよい。

また、本開示による流体解析装置においては、表示制御部は、第１の結果および第２の結果のいずれか一方を表示し、一致度に応じた表示態様により、他方の結果をさらに表示するものであってもよい。

また、本開示による流体解析装置においては、表示制御部は、一致度が小さいほど他方の結果を強調表示するものであってもよい。

また、本開示による流体解析装置においては、表示制御部は、一致度が予め定められたしきい値以上となる位置において、第１の結果および第２の結果のいずれか一方を表示し、一致度がしきい値未満となる位置において、一致度に応じた表示態様により、他方の結果をさらに表示するものであってもよい。

また、本開示による流体解析装置においては、画像は、３次元シネ位相コントラスト磁気共鳴法によって被写体を撮影することにより取得された３次元画像であり、
第１解析部は、３次元画像を解析することにより取得される流体の流速ベクトルを、第１の結果として導出するものであってもよい。

また、本開示による流体解析装置においては、第２解析部は、数値流体力学を用いた解析により流体の流れをシミュレーションすることにより取得される流体の流速ベクトルを、第２の結果として導出するものであってもよい。

また、本開示による流体解析装置においては、構造物が血管であり、流体が血液であってもよい。

また、本開示による流体解析装置においては、表示制御部は、表示部に構造物の形態画像を表示し、形態画像上において、一致度に応じて第１の結果および第２の結果の少なくとも一方を表示するものであってもよい。

本開示による流体解析方法は、内部に流体が流れる構造物を含む被写体を撮影することにより取得した画像を解析して、構造物内の各位置における流体の流れに関する第１の結果を導出し、
構造物内の各位置における流体の流れをシミュレーションして、流体の流れに関する第２の結果を導出し、
第１の結果および第２の結果のいずれか一方を基準とした、第１の結果と第２の結果との対応する位置における一致度を導出し、
一致度に応じて第１の結果および第２の結果の少なくとも一方を表示部に表示する。

なお、本開示による流体解析方法をコンピュータに実行させるためのプログラムとして提供してもよい。

本開示による他の流体解析装置は、コンピュータに実行させるための命令を記憶するメモリと、
記憶された命令を実行するよう構成されたプロセッサとを備え、プロセッサは、
内部に流体が流れる構造物を含む被写体を撮影することにより取得した画像を解析して、構造物内の各位置における流体の流れに関する第１の結果を導出し、
構造物内の各位置における流体の流れをシミュレーションして、流体の流れに関する第２の結果を導出し、
第１の結果および第２の結果のいずれか一方を基準とした、第１の結果と第２の結果との対応する位置における一致度を導出し、
一致度に応じて第１の結果および第２の結果の少なくとも一方を表示部に表示する処理を実行する。

本開示によれば、画像を解析することにより導出された結果とシミュレーションにより導出された結果との相違を容易に認識することができる。

本開示の実施形態による流体解析装置を適用した、診断支援システムの概要を示すハードウェア構成図 本開示の実施形態による流体解析装置の概略構成を示す図 ３次元シネ位相コントラスト磁気共鳴法によって撮影された第１の３次元画像を示す図 同一の血管領域における第１の結果および第２の結果を示す図 補間演算により第２の結果と空間分解能を一致させた第１の結果を示す図 一致度に応じた第１の結果および第２の結果の少なくとも一方の表示画面を示す図 血管の形態画像において、第１の結果および第２の結果の少なくとも一方が表示された状態を示す図 本実施形態において行われる処理を示すフローチャート 一致度に応じた第１の結果および第２の結果の少なくとも一方の表示画面の他の例を示す図 一致度に応じた第１の結果および第２の結果の少なくとも一方の表示画面の他の例を示す図 一致度に応じた第１の結果および第２の結果の少なくとも一方の表示画面の他の例を示す図 一致度に応じた第１の結果および第２の結果の少なくとも一方の表示画面の他の例を示す図

以下、図面を参照して本開示の実施形態について説明する。図１は、本開示の実施形態による流体解析装置を適用した、診断支援システムの概要を示すハードウェア構成図である。図１に示すように、診断支援システムでは、本実施形態による流体解析装置１、３次元画像撮影装置２、および画像保管サーバ３が、ネットワーク４を経由して通信可能な状態で接続されている。

３次元画像撮影装置２は、被検体の診断対象となる部位を撮影することにより、その部位を表す３次元画像を生成する装置であり、具体的には、ＣＴ装置、ＭＲＩ装置、およびＰＥＴ（Positron Emission Tomography）装置等である。３次元画像撮影装置２により生成された３次元画像は画像保管サーバ３に送信され、保存される。なお、本実施形態では、患者の心臓（本開示の被写体に相当する）の３次元画像を取得する場合について説明するが、これに限らず、肺、肝臓および頭部等その他の臓器でもよい。とくに、本実施形態においては、複数種類の３次元画像を取得するものとする。具体的には、ＭＲＩ装置において３次元シネ位相コントラスト磁気共鳴法によって被写体を撮影したＭＲＩ画像を第１の３次元画像Ｇ１として取得し、ＣＴ装置において造影剤を用いて被写体を撮影した造影ＣＴ画像を第２の３次元画像Ｇ２として取得するものとする。しかしながら、取得される３次元画像の種類はこれらに限定されるものではない。また、心臓の血管が本開示の構造物に、血液が本開示の流体にそれぞれ対応する。

画像保管サーバ３は、各種データを保存して管理するコンピュータであり、大容量外部記憶装置およびデータベース管理用ソフトウェアを備えている。画像保管サーバ３は、有線あるいは無線のネットワーク４を介して他の装置と通信を行い、画像データ等を送受信する。具体的には３次元画像撮影装置２で生成された３次元画像の画像データを含む各種データをネットワーク経由で取得し、大容量外部記憶装置等の記録媒体に保存して管理する。なお、画像データの格納形式およびネットワーク４経由での各装置間の通信は、ＤＩＣＯＭ（Digital Imaging and Communication in Medicine）等のプロトコルに基づいている。

流体解析装置１は、１台のコンピュータに、本開示の流体解析プログラムをインストールしたものである。コンピュータは、診断を行う医師が直接操作するワークステーションまたはパーソナルコンピュータでもよいし、それらとネットワークを介して接続されたサーバコンピュータでもよい。流体解析プログラムは、ＤＶＤ（Digital Versatile Disc）あるいはＣＤ−ＲＯＭ（Compact Disc Read Only Memory）等の記録媒体に記録されて配布され、その記録媒体からコンピュータにインストールされる。または、ネットワークに接続されたサーバコンピュータの記憶装置、もしくはネットワークストレージに、外部からアクセス可能な状態で記憶され、要求に応じて医師が使用するコンピュータにダウンロードされ、インストールされる。

図２は、コンピュータに流体解析プログラムをインストールすることにより実現される流体解析装置の概略構成を示す図である。図２に示すように、流体解析装置１は、標準的なワークステーションの構成として、ＣＰＵ（Central Processing Unit）１１、メモリ１２およびストレージ１３を備えている。また、流体解析装置１には、液晶ディスプレイ等の表示部１４、並びにキーボードおよびマウス等の入力部１５が接続されている。

ストレージ１３は、ハードディスクドライブ等からなり、ネットワーク４を経由して画像保管サーバ３から取得した３次元画像、並びに処理に必要な情報を含む各種情報が記憶されている。

また、メモリ１２には、流体解析プログラムが記憶されている。流体解析プログラムは、ＣＰＵ１１に実行させる処理として、被写体の第１および第２の３次元画像Ｇ１，Ｇ２を取得する画像取得処理、第１の３次元画像Ｇ１を解析して、血管内の各位置における血液の流れに関する第１の結果を導出する第１解析処理、第２の３次元画像Ｇ２を用いて血管内の各位置における血液の流れをシミュレーションして、血液の流れに関する第２の結果を導出する第２解析処理、第１の結果および第２の結果のいずれか一方を基準とした、第１の結果と第２の結果との対応する位置における一致度を導出する一致度導出処理、並びに一致度に応じて第１の結果および第２の結果の少なくとも一方を表示部１４に表示する表示制御処理を規定する。

そして、ＣＰＵ１１がプログラムに従いこれらの処理を実行することで、コンピュータは、画像取得部２０、第１解析部２１、第２解析部２２、一致度導出部２３および表示制御部２４として機能する。

画像取得部２０は、第１および第２の３次元画像Ｇ１，Ｇ２を画像保管サーバ３から取得する。なお、第１および第２の３次元画像Ｇ１，Ｇ２が既にストレージ１３に記憶されている場合には、画像取得部２０は、ストレージ１３から第１および第２の３次元画像Ｇ１，Ｇ２を取得するようにしてもよい。

第１解析部２１は、第１の３次元画像Ｇ１を解析して、血管内の各位置における血液の流れに関する第１の結果Ｒ１を導出する。本実施形態においては、第１解析部２１は、まず、第１の３次元画像Ｇ１から血管領域を抽出する。具体的には、第１解析部２１は、第１の３次元画像Ｇ１に対して多重解像度変換を行い、各解像度の画像に対してヘッセ行列の固有値解析を行い、各解像度の画像における解析結果を統合することによって、第１の３次元画像Ｇ１に含まれる心臓領域中の様々なサイズの線構造（血管）の集合体として、冠動脈の領域を血管領域として抽出する（例えばY Sato, et al.、「Three-dimensional multi-scale line filter for segmentation and visualization of curvilinear structures in medical images.」、Medical Image Analysis、1998年6月、Vol.2、No.2、p.p.143-168等参照）。また、さらに最小全域木アルゴリズム等を用いて、抽出された各線構造の中心点を連結することにより、冠動脈を表す木構造のデータを生成し、抽出された冠動脈の中心点を結ぶ芯線上の各点（木構造データの各ノード）において、芯線に直交する断面を求め、各断面において、グラフカット法等の公知のセグメンテーション手法を用いて冠動脈の輪郭を認識し、その輪郭を表す情報を木構造データの各ノードに関連づけることによって、冠動脈の領域を血管領域として抽出するようにしてもよい。

なお、血管領域の抽出方法としては上記の方法に限らず、領域拡張法等のその他の公知の手法を用いるようにしてもよい。

そして、第１解析部２１は、第１の３次元画像Ｇ１から抽出された血管領域内の速度情報を用いて、血管内の各ボクセル位置における流速ベクトルを第１の結果Ｒ１として導出する。図３は、３次元シネ位相コントラスト磁気共鳴法によって撮影された第１の３次元画像Ｇ１を示す図である。図３に示すように、３次元シネ位相コントラスト磁気共鳴法によって撮影された第１の３次元画像Ｇ１の画像データは、マグニチュードデータＭ、Ｘ方向の位相データＰｈｘ、Ｙ軸方向の位相データＰｈｙ、およびＺ軸方向の位相データＰｈｚを、時間ｔに沿って所定の周期（例えば心周期）で得た３次元データを含む。Ｘ方向の位相データＰｈｘ、Ｙ軸方向の位相データＰｈｙ、およびＺ軸方向の位相データＰｈｚは、マグニチュードデータＭをＸ軸方向、Ｙ軸方向、およびＺ軸方向にエンコード（ＶＥＮＣ：velocity encoding）することにより生成される。Ｘ方向の位相データＰｈｘ、Ｙ軸方向の位相データＰｈｙ、およびＺ軸方向の位相データＰｈｚは各軸方向の流速を表すデータである。第１解析部２１は、３つの位相データから第１の３次元画像Ｇ１の各ボクセル位置の３次元流速ベクトル（以下、流速ベクトルとする）を第１の結果Ｒ１として導出する。なお、流速ベクトルは、第１の３次元画像Ｇ１の各ボクセル位置において導出しているが、これに限定されるものではなく、予め定められたボクセル間隔（例えば５ボクセル、１０ボクセル等）にて導出してもよい。また、複数ボクセルからなる領域において１つの流速ベクトルを導出するようにしてもよい。

なお、画像取得部２０において、ドップラー計測によって時系列に撮影された３次元の超音波画像を取得し、その超音波画像に基づいて取得された血管領域内の速度情報を用いて流速ベクトルを取得して、第１の結果Ｒ１を導出するようにしてもよい。

第２解析部２２は、第２の３次元画像Ｇ２を用いて、血管内の各位置における血液の流れをシミュレーションして、血液の流れに関する第２の結果Ｒ２を導出する。このために、第２解析部２２は、まず第２の３次元画像Ｇ２から血管領域を抽出する。血管領域の抽出は、上述した第１解析部２１と同様に行えばよい。そして、第２解析部２２は、抽出された血管領域を用いてＣＦＤ（Computational Fluid Dynamics）による血流解析を行うことにより、第２の３次元画像Ｇ２の各ボクセル位置における流速ベクトルを第２の結果Ｒ２として取得する。なお、第２の結果Ｒ２についても、第２の３次元画像Ｇ２における予め定められたボクセル間隔にて導出してもよく、複数ボクセルからなる領域において１つの流速ベクトルを導出するようにしてもよい。

なお、第１の３次元画像Ｇ１は、ＭＲＩ装置における制約により、第２の３次元画像Ｇ２よりも空間分解能および時間分解能が低い。このため、空間的に見ると、第１の結果Ｒ１は第２の結果Ｒ２よりも大きい間隔にて導出される。図４は同一の血管領域における第１の結果および第２の結果を示す図である。図４に示すように、第１の結果Ｒ１の流速ベクトルは、第２の結果Ｒ２の流速ベクトルと比較すると、その間隔は第１の結果Ｒ１の方が空間的に広いものとなる。

一致度導出部２３は、第１の結果Ｒ１および第２の結果Ｒ２のいずれか一方を基準とした、第１の結果Ｒ１と第２の結果Ｒ２との対応するボクセル位置における一致度Ｃ０を導出する。なお、本実施形態においては、第２の結果Ｒ２を基準とした一致度Ｃ０を導出するものとするが、これに限定されるものではない。ここで、上述したように、第１の結果Ｒ１の流速ベクトルの空間分解能は、第２の結果Ｒ２の流速ベクトルの空間分解能よりも低い。このため、一致度導出部２３は、第１の結果Ｒ１について第２の結果Ｒ２と同一の空間分解能となるように、第１の結果Ｒ１と第２の結果Ｒ２との空間分解能を一致させる。具体的には、第２の結果Ｒ２において流速ベクトルが導出されたすべてのボクセル位置における第１の結果Ｒ１の流速ベクトルを、第１解析部２１が導出した第１の結果Ｒ１の流速ベクトルを補間することにより導出する。図５は補間演算により第２の結果Ｒ２と空間分解能を一致させた第１の結果Ｒ１を示す図である。

そして、一致度導出部２３は、第１の結果Ｒ１と第２の結果Ｒ２との対応するボクセル位置における流速ベクトルの内積を一致度Ｃ０として導出する。なお、一致度Ｃ０は、対応するボクセル位置における流速ベクトル同士で導出してもよいが、各ボクセル位置を基準として予め定められた範囲内の流速ベクトルの平均値を各ボクセル位置の流速ベクトルとして算出し、算出した流速ベクトルを用いて一致度Ｃ０を導出してもよい。なお、平均値は、注目するボクセル位置から離れるほど重み付けを小さくする重み付け平均値としてもよい。

表示制御部２４は、一致度導出部２３が導出した一致度Ｃ０に応じて、第１の結果Ｒ１および第２の結果Ｒ２の少なくとも一方を表示部１４に表示する。本実施形態においては、第２の結果Ｒ２を表示するものとし、表示された第２の結果Ｒ２の表示態様を一致度Ｃ０に応じて変更する。なお、第１の結果Ｒ１を表示し、表示された第１の結果Ｒ１の表示態様を一致度Ｃ０に応じて変更してもよい。図６は一致度Ｃ０に応じた第１の結果および第２の結果の表示画面を示す図である。図６に示すように、表示画面３０においては、第２の結果Ｒ２は、第１の結果Ｒ１との一致度Ｃ０が大きいほど、流速ベクトルの濃度が大きいものとなっている。表示画面３０においては、右側ほど第１の結果Ｒ１と第２の結果Ｒ２との一致度Ｃ０が大きく、左に向かうほど一致度Ｃ０が小さいものとなっている。

なお、実際には、図７に示すように、表示部１４に表示された血管の形態画像３１を表示した表示画面３０Ａにおいて、第２の結果Ｒ２の流速ベクトルが、一致度Ｃ０に応じた表示態様で表示されることとなる。ここで、図７においては、説明のために流速ベクトルの間隔を大きくしているが、実際には第１の結果Ｒ１および第２の結果Ｒ２を取得したボクセル位置毎または複数のボクセルからなる領域毎に、流速ベクトルが表示されることとなる。

また、図６においては、一致度Ｃ０に応じて流速ベクトルの濃度を変えているが、色を変えるようにしてもよく、透明度を変えるようにしてもよく、流速ベクトルの線の種類（破線および実線等）を変えるようにしてもよい。また、これらを組み合わせて表示態様を変更してもよい。

次いで、本実施形態において行われる処理について説明する。図８は本実施形態において行われる処理を示すフローチャートである。まず、画像取得部２０が第１および第２の３次元画像Ｇ１，Ｇ２を画像保管サーバ３から取得する（画像取得；ステップＳＴ１）。そして、第１解析部２１が、第１の３次元画像Ｇ１を解析して、血管内の各位置における血液の流れに関する第１の結果Ｒ１を導出する（ステップＳＴ２）。また、第２解析部２２が、血管内の各位置における血液の流れをシミュレーションして、血液の流れに関する第２の結果Ｒ２を導出する（ステップＳＴ３）。なお、ステップＳＴ３の処理を先に行ってもよく、ステップＳＴ２の処理とステップＳＴ３の処理とを並列に行ってもよい。

次いで、一致度導出部２３が、第１の結果および第２の結果のいずれか一方を基準とした、第１の結果Ｒ１と第２の結果Ｒ２との対応する位置における一致度Ｃ０を導出する（ステップＳＴ４）。そして、表示制御部２４が、一致度Ｃ０に応じて第１の結果Ｒ１および第２の結果Ｒ２の少なくとも一方を表示部１４に表示し（結果表示；ステップＳＴ５）、処理を終了する。

このように、本実施形態によれば、第１の結果Ｒ１および第２の結果Ｒ２のいずれか一方を基準とした、第１の結果Ｒ１と第２の結果Ｒ２との対応する位置における一致度Ｃ０を導出し、一致度Ｃ０に応じて第１の結果Ｒ１および第２の結果Ｒ２の少なくとも一方を表示するようにした。このため、ユーザは、表示された第１の結果Ｒ１および第２の結果Ｒ２の少なくとも一方を見ることにより、血管内の各位置において、第１の結果Ｒ１と第２の結果Ｒ２とがどの程度一致しているかを認識できる。したがって、本実施形態によれば、第１の結果Ｒ１と第２の結果Ｒ２との相違、すなわち、第１の３次元画像Ｇ１画像を解析することにより導出された結果とシミュレーションにより導出された結果との相違を容易に認識することができる。

なお、上記実施形態においては、図６に示すように一致度Ｃ０に応じて表示態様が変更された第２の結果Ｒ２の表示画面３０を表示部１４に表示しているが、これに限定されるものではない。例えば、図９に示すように、表示制御部２４は、第１の結果Ｒ１および第２の結果Ｒ２のいずれを表示するかを切り替えるためのバー４０を含む表示画面３２を表示部１４に表示してもよい。図９に示す表示画面３２においては、ＣＦＤにより導出した第２の結果Ｒ２を表示するようにバー４０が切り替えられている。バー４０を４Ｄフローである第１の結果Ｒ１を表示するように切り替えた場合、図１０の表示画面３３に示すように、第１の結果Ｒ１が表示され、第２の結果Ｒ２との一致度Ｃ０が大きいほど流速ベクトルの濃度が濃いものとなる。

また、表示制御部２４は、第１の結果Ｒ１および第２の結果Ｒ２のいずれか一方を表示し、一致度Ｃ０に応じた表示態様により、他方の結果をさらに表示するようにしてもよい。例えば、図１１の表示画面３４に示すように、第２の結果Ｒ２を表示し、一致度Ｃ０に応じた表示態様により第１の結果Ｒ１をさらに表示してもよい。なお、第１の結果Ｒ１を表示し、一致度Ｃ０に応じた表示態様により第２の結果Ｒ２を表示してもよい。図１１に示す表示画面３４おいては、図６と同様に、第２の結果Ｒ２である流速ベクトルが一致度Ｃ０に応じた表示態様により表示され、さらに、一致度Ｃ０が小さいほど、第１の結果Ｒ１である流速ベクトルが一致度Ｃ０に応じた表示態様にて表示されている。なお、図１１に示す表示画面３４において、第１の結果Ｒ１の流速ベクトルは、一致度Ｃ０が小さいほど太い破線により強調表示されている。したがって、表示画面３４には、徐々に一致度Ｃ０が小さくなることを認識可能に、第２の結果Ｒ２と併せて第１の結果Ｒ１が表示されることとなる。

なお、第１の結果Ｒ１の流速ベクトルの表示態様の変更については、一致度Ｃ０に応じて破線の太さを変えているが、色を変えるようにしてもよく、透明度を変えるようにしてもよい。また、これらを組み合わせて表示態様を変更してもよい。

また、図１２の表示画面３５に示すように、一致度Ｃ０が小さい位置については、第１の結果Ｒ１および第２の結果Ｒ２のいずれを表示するかを選択するためのバー４５を表示するようにしてもよい。図１２においては、バー４５において４Ｄフロー、すなわち第１の結果Ｒ１が選択されているため、例えば一致度Ｃ０が予め定められたしきい値Ｔｈ１以下であるように、一致度Ｃ０が小さい位置については、第２の結果Ｒ２に代えて第１の結果Ｒ１の流速ベクトルが表示される。なお、図１２に示す表示画面３５おいては、第１の結果Ｒ１の流速ベクトルを破線の矢印により示している。なお、図１２に示す表示画面３５おいて、バー４５をＣＦＤの側に切り替えれば、図６と同様に第２の結果Ｒ２のみが表示される。なお、しきい値Ｔｈ１は、入力部１５からの入力により、任意の値に変更できるようにしてもよい。

また、上記実施形態においては、血管内の各位置における流速ベクトルを第１の結果Ｒ１および第２の結果Ｒ２として導出しているが、これに限定されるものではない。流速ベクトルの他、例えば流体の圧力、壁面せん断応力（ＷＳＳ(Wall Shear Stress)）、渦度およびヘリシティ等を第１の結果Ｒ１および第２の結果Ｒ２として用いてもよい。

また、上記実施形態においては、内部に流体が流れる構造物として血管を用いているがこれに限定されるものではない。例えば、脳脊髄液の流れを可視化することを考えた場合、内部に髄液が流れる構造物として、頭蓋内での脳室、とくにくも膜下腔を、脊柱管内では脊髄くも膜下腔を、内部に流体が流れる構造物として用いてもよい。また、内部にリンパ液が流れるリンパ管を用いてもよい。

また、上記実施形態においては、人体を対象とした画像を対象としているが、これに限定されるものではない。例えば、配管内を流れる流体の流れを解析する際にも、本開示を適用できることはもちろんである。

また、上記実施形態において、例えば、画像取得部２０、第１解析部２１、第２解析部２２、一致度導出部２３および表示制御部２４といった各種の処理を実行する処理部（Processing Unit）のハードウェア的な構造としては、次に示す各種のプロセッサ（Processor）を用いることができる。上記各種のプロセッサには、上述したように、ソフトウェア（プログラム）を実行して各種の処理部として機能する汎用的なプロセッサであるＣＰＵに加えて、ＦＰＧＡ（Field Programmable Gate Array）等の製造後に回路構成を変更可能なプロセッサであるプログラマブルロジックデバイス（Programmable Logic Device :PLD）、ＡＳＩＣ（Application Specific Integrated Circuit）等の特定の処理を実行させるために専用に設計された回路構成を有するプロセッサである専用電気回路等が含まれる。

１つの処理部は、これらの各種のプロセッサのうちの１つで構成されてもよいし、同種または異種の２つ以上のプロセッサの組み合わせ（例えば、複数のＦＰＧＡの組み合わせまたはＣＰＵとＦＰＧＡとの組み合わせ）で構成されてもよい。また、複数の処理部を１つのプロセッサで構成してもよい。

複数の処理部を１つのプロセッサで構成する例としては、第１に、クライアントおよびサーバ等のコンピュータに代表されるように、１つ以上のＣＰＵとソフトウェアとの組み合わせで１つのプロセッサを構成し、このプロセッサが複数の処理部として機能する形態がある。第２に、システムオンチップ（System On Chip:SoC）等に代表されるように、複数の処理部を含むシステム全体の機能を１つのＩＣ（Integrated Circuit）チップで実現するプロセッサを使用する形態がある。このように、各種の処理部は、ハードウェア的な構造として、上記各種のプロセッサの１つ以上を用いて構成される。

さらに、これらの各種のプロセッサのハードウェア的な構造としては、より具体的には、半導体素子等の回路素子を組み合わせた電気回路（Circuitry）を用いることができる。

１ 流体解析装置
２ ３次元画像撮影装置
３ 画像保管サーバ
４ ネットワーク
１１ ＣＰＵ
１２ メモリ
１３ ストレージ
１４ 表示部
１５ 入力部
２０ 画像取得部
２１ 第１解析部
２２ 第２解析部
２３ 一致度導出部
２４ 表示制御部
３０，３０Ａ，３２〜３５ 表示画面
３１ 形態画像
４０，４５ バー
Ｍ マグニチュードデータ
Ｐｈｘ 位相データ
Ｐｈｙ 位相データ
Ｐｈｚ 位相データ
Ｒ１ 第１の結果
Ｒ２ 第２の結果

Claims (11)

1. 内部に流体が流れる構造物を含む被写体を撮影することにより取得した画像を解析して、前記構造物内の各位置における前記流体の流れに関する第１の結果を導出する第１解析部と、
   前記構造物内の各位置における前記流体の流れをシミュレーションして、前記流体の流れに関する第２の結果を導出する第２解析部と、
   前記第１の結果および前記第２の結果のいずれか一方を基準とした、前記第１の結果と前記第２の結果との対応する位置における一致度を導出する一致度導出部と、
   前記一致度に応じて前記第１の結果および前記第２の結果の少なくとも一方を表示部に表示する表示制御部とを備えた流体解析装置。
2. 前記表示制御部は、前記第１の結果および前記第２の結果のいずれか一方を表示し、前記一致度に応じて、表示された前記第１の結果および前記第２の結果のいずれか一方の表示態様を変更する請求項１に記載の流体解析装置。
3. 前記表示制御部は、前記第１の結果および前記第２の結果のいずれか一方を表示し、前記一致度に応じた表示態様により、他方の結果をさらに表示する請求項１または２に記載の流体解析装置。
4. 前記表示制御部は、前記一致度が小さいほど前記他方の結果を強調表示する請求項３に記載の流体解析装置。
5. 前記表示制御部は、前記一致度が予め定められたしきい値以上となる位置において、前記第１の結果および前記第２の結果のいずれか一方を表示し、前記一致度が前記しきい値未満となる位置において、前記一致度に応じた表示態様により、他方の結果をさらに表示する請求項１または２に記載の流体解析装置。
6. 前記画像は、３次元シネ位相コントラスト磁気共鳴法によって前記被写体を撮影することにより取得された３次元画像であり、
   前記第１解析部は、前記３次元画像を解析することにより取得される前記流体の流速ベクトルを、前記第１の結果として導出する請求項１から５のいずれか１項に記載の流体解析装置。
7. 前記第２解析部は、数値流体力学を用いた解析により前記流体の流れをシミュレーションすることにより取得される前記流体の流速ベクトルを、前記第２の結果として導出する請求項１から６のいずれか１項に記載の流体解析装置。
8. 前記構造物が血管であり、前記流体が血液である請求項１から７のいずれか１項に記載の流体解析装置。
9. 前記表示制御部は、前記表示部に前記構造物の形態画像を表示し、該形態画像上において、前記一致度に応じて前記第１の結果および前記第２の結果の少なくとも一方を表示する請求項１から８のいずれか１項に記載の流体解析装置。
10. 内部に流体が流れる構造物を含む被写体を撮影することにより取得した画像を解析して、前記構造物内の各位置における前記流体の流れに関する第１の結果を導出し、
    前記構造物内の各位置における前記流体の流れをシミュレーションして、前記流体の流れに関する第２の結果を導出し、
    前記第１の結果および前記第２の結果のいずれか一方を基準とした、前記第１の結果と前記第２の結果との対応する位置における一致度を導出し、
    前記一致度に応じて前記第１の結果および前記第２の結果の少なくとも一方を表示部に表示する流体解析方法。
11. 内部に流体が流れる構造物を含む被写体を撮影することにより取得した画像を解析して、前記構造物内の各位置における前記流体の流れに関する第１の結果を導出する手順と、
    前記構造物内の各位置における前記流体の流れをシミュレーションして、前記流体の流れに関する第２の結果を導出する手順と、
    前記第１の結果および前記第２の結果のいずれか一方を基準とした、前記第１の結果と前記第２の結果との対応する位置における一致度を導出する手順と、
    前記一致度に応じて前記第１の結果および前記第２の結果の少なくとも一方を表示部に表示する手順とをコンピュータに実行させる流体解析プログラム。

Images