(57)【要約】 【課題】トモグラフィックＰＩＶ計測法により可視化流体の３次元速度ベクトル又は３速度成分を取得するための演算処理装置の負荷及びデータ解析時間を軽減又は短縮する。 【解決手段】撮像装置(1-4)によって撮像された２次元画像データに基づいて３次元空間(α)内の各時刻(t,t＋Δt)の粒子分布が再構築され(S3)、連続２時刻の３次元輝度情報が得られる。３次元空間を分割してなる所定厚のスライス領域(β)が設定され、スライス領域内の粒子分布が３次元輝度情報より抽出される(S4)。スライス領域内の粒子分布は、複数の仮想投影面に投影され、連続２時刻の複数且つ仮想の２次元粒子画像が生成される(S6)。複数の２次元粒子画像に２次元ＰＩＶを実行して得られる複数の画像面内の粒子の２次元移動ベクトルに基づいて、３次元空間における可視化流体の３次元速度ベクトル又は３速度成分が取得される(S7,S8,S11)。【選択図】図７
2015161639.tif 000002
【特許請求の範囲】
【請求項１】
微小粒子群により可視化された可視化流体を撮像する複数の撮像装置によって、微小時間間隔を隔てた連続２時刻の粒子像を撮像し、撮像により得られた粒子像の画像データに基づいて３次元空間内の粒子分布を再構築し、前記３次元空間における可視化流体の３次元速度ベクトル又は３速度成分を測定する可視化流体の流速計測方法において、
前記撮像装置によって撮像された２次元画像データに基づいて前記３次元空間内の各時刻の粒子分布を再構築し、連続２時刻の３次元輝度情報を得る工程と、
前記３次元空間を分割してなる所定厚のスライス領域を設定し、該スライス領域の内部の粒子分布を前記３次元輝度情報より抽出して、該粒子分布を複数の仮想投影面に投影することによって、連続２時刻の複数且つ仮想の２次元粒子画像を生成する工程と、
複数の前記２次元粒子画像に２次元ＰＩＶを実行して得られる複数の画像面内の粒子の２次元移動ベクトルに基づいて、前記３次元空間における可視化流体の３次元速度ベクトル又は３速度成分を取得する工程とを有することを特徴とする流速計測方法。 
【請求項２】
微小粒子群により可視化された可視化流体を撮像する複数の撮像装置によって、微小時間間隔を隔てた連続２時刻の粒子像を撮像し、撮像により得られた粒子像の画像データに基づいて３次元空間内の粒子分布を再構築し、前記３次元空間における可視化流体の３次元速度ベクトル又は３速度成分を測定する可視化流体の流速計測方法において、
前記３次元空間を分割してなる所定厚のスライス領域を設定し、前記撮像装置によって撮像された２次元画像データに基づいて前記スライス領域内の各時刻の粒子分布を再構築し、前記スライス領域内における連続２時刻の３次元輝度情報を得る工程と、
前記３次元輝度情報に基づいて前記スライス領域の内部の粒子分布を複数の仮想投影面に投影することによって、連続２時刻の複数且つ仮想の２次元粒子画像を生成する工程と、
複数の前記２次元粒子画像に２次元ＰＩＶを実行して得られる複数の画像面内の粒子の２次元移動ベクトルに基づいて、前記スライス領域における可視化流体の３次元速度ベクトル又は３速度成分を取得する工程と、
各スライス領域における前記３次元速度ベクトル又は３速度成分の情報を集合、統合又は合成して、前記３次元空間における可視化流体の３次元速度ベクトル又は３速度成分を取得する工程とを有することを特徴とする流速計測方法。
【請求項３】
前記スライス領域は、前記３次元空間を等分割した空間であることを特徴とする請求項１又は２に記載の流速計測方法。
【請求項４】
前記スライス領域は、前記可視化流体に照射されるレーザー光の光軸と平行に前記可視化流体の流路を横断する空間であることを特徴とする請求項１乃至３のいずれか１項に記載の流速計測方法。
【請求項５】
前記仮想投影面は、前記撮像装置の光軸と直交する平面であることを特徴とする請求項１乃至４のいずれか１項に記載の流速計測方法。
【請求項６】
前記３次元輝度情報を得る前記工程において、前記２次元画像の画素に基づいて前記スライス領域内の前記３次元輝度情報を重み付け累算で特定するための重み係数を算出し且つ記憶し、
前記２次元粒子画像を生成する前記工程において、前記スライス領域内の粒子分布を前記重み係数に従って前記仮想投影面に投影することを特徴とする請求項２に記載の流速計測方法。
【請求項７】
微小粒子群により可視化された可視化流体を撮像するように配置され、微小時間間隔を隔てた連続２時刻の粒子像を撮像する複数の撮像装置と、該撮像装置によって得られた画像データに基づいて３次元空間内の粒子分布を再構築し、前記３次元空間における可視化流体の３次元速度ベクトル又は３速度成分を演算する演算処理装置とを有する可視化流体の流速計測システムにおいて、
前記演算処理装置は、前記撮像装置によって撮像された２次元画像データに基づいて前記３次元空間内の各時刻の粒子分布を再構築し、連続２時刻の３次元輝度情報を得るように構成された３次元輝度情報取得手段と、
前記３次元空間を分割してなる所定厚のスライス領域を設定し、該スライス領域の内部の粒子分布を前記３次元輝度情報より抽出して、該粒子分布を複数の仮想投影面に投影することによって、連続２時刻の複数且つ仮想の２次元粒子画像を生成する２次元粒子画像生成手段と、
複数の前記２次元粒子画像に２次元ＰＩＶを実行し、これにより得られる複数の画像面内の粒子の２次元移動ベクトルに基づいて、前記３次元空間における可視化流体の３次元速度ベクトル又は３速度成分を演算する演算手段とを有することを特徴とする流速計測システム。 
【請求項８】
微小粒子群により可視化された可視化流体を撮像するように配置され、微小時間間隔を隔てた連続２時刻の粒子像を撮像する複数の撮像装置と、該撮像装置によって得られた画像データに基づいて３次元空間内の粒子分布を再構築し、前記３次元空間における可視化流体の３次元速度ベクトル又は３速度成分を演算する演算処理装置とを有する可視化流体の流速計測システムにおいて、
前記演算処理装置は、前記３次元空間を分割してなる所定厚のスライス領域を設定し、前記撮像装置によって撮像された２次元画像データに基づいて前記スライス領域内の各時刻の粒子分布を再構築し、連続２時刻の３次元輝度情報を得るように構成された３次元輝度情報取得手段と、
前記スライス領域の内部の粒子分布を複数の仮想投影面に投影することによって、該スライス領域内における連続２時刻の複数且つ仮想の２次元粒子画像を生成する２次元粒子画像生成手段と、
複数の前記２次元粒子画像に２次元ＰＩＶを実行し、これにより得られる複数の画像面内の粒子の２次元移動ベクトルに基づいて、前記スライス領域における可視化流体の３次元速度ベクトル又は３速度成分を演算する第１演算手段と、
各スライス領域における前記３次元速度ベクトル又は３速度成分の情報を集合、統合又は合成して、前記３次元空間における可視化流体の３次元速度ベクトル又は３速度成分を取得する第２演算手段とを有することを特徴とする流速計測システム。 
【請求項９】
前記スライス領域は、前記３次元空間を等分割した空間であることを特徴とする請求項７又は８に記載の流速計測システム。
【請求項１０】
前記スライス領域は、前記可視化流体に照射されるレーザー光の光軸と平行に前記可視化流体の流路を横断する空間であることを特徴とする請求項７乃至９のいずれか１項に記載の流速計測システム。
【請求項１１】
前記仮想投影面は、前記撮像装置の光軸と直交する平面であることを特徴とする請求項７乃至１０のいずれか１項に記載の流速計測システム。
【請求項１２】
前記３次元輝度情報取得手段は、前記２次元画像の画素に基づいて前記スライス領域内の前記３次元輝度情報を重み付け累算で特定するための重み係数を算出し且つ記憶し、
前記２次元粒子画像生成手段は、前記スライス領域内の粒子分布を前記重み係数に従って前記仮想投影面に投影することを特徴とする請求項８に記載の流速計測システム。
【請求項１３】
前記可視化流体の流路は、円形の流路断面を有し、画像座標と物理座標との関係を求めるために校正画像を取得すべく前記流路内に配置される円形校正板が用いられ、該校正板は、中心円形領域に配設された直交格子配列の指標と、外周領域に周方向に配設された環状配列の指標とを有することを特徴とする請求項１乃至６のいずれか１項に記載の流速計測方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、可視化流体の流速計測方法及び流速計測システムに関するものであり、より詳細には、微小粒子群により流体を可視化してなる可視化流体の流動場を撮像装置によって撮影し、可視化流体の３次元速度ベクトル又は３速度成分を計測する流速計測方法及び流速計測システムに関するものである。
【背景技術】
【０００２】
流体を可視化するためのトレーサ粒子をマーカとして流体の流れに混入し、流動場の流体速度等を計測する可視化計測技術が知られている。可視化流体計測技術として知られる直接撮影法においては、二重露光撮影又は高速度カメラ撮影や、パルスレーザ及びデジタルＣＣＤカメラを用いたフレームまたぎ撮影等の方法によって、トレーサ粒子を微小時間間隔（例えば、１ms間隔)で連続撮影することにより、粒子の移動距離が測定される。粒子の速度は、移動距離を時間間隔で除すことによって演算される。
【０００３】
可視化流体計測技術の一種として知られたＰＩＶ（Particle Image Velocimetry、粒子像流速計測）は、このような可視化流体の速度分布を調べる直接撮影法として普及した計測技法である。ＰＩＶ計測法は、流体にレーザシート光を照射して、流体に含まれる寸法１０μm程度の微小粒子を微小時間間隔(時刻t及び時刻t＋Δt)の連続撮影により撮像し、粒子又は粒子群の速度を画像解析により求め、これにより、流速分布等を測定する流速計測方法である。このようなＰＩＶ計測法によれば、レーザシート光内における面内２成分の速度成分を計測することができる。
【０００４】
また、レーザシート光（２次元）内の速度３成分を得る計測技法として、ステレオＰＩＶ計測法が知られている。ステレオＰＩＶ計測法は、異なる角度方向から可視化流体を撮影する少なくとも２台のＣＣＤカメラによって流体中の微粒子を微小時間間隔で連続撮影し、複数のカメラ映像の視差に基づいて、粒子又は粒子群の三次元移動量を画像解析により求め、これにより、可視化流体の速度３成分（ｘ軸、ｙ軸及びｚ軸方向）を計測する流速計測方法である。
【０００５】
このようなＰＩＶ計測法又はステレオＰＩＶ計測法は、例えば、特開２０１１−２４７６０１号公報、特開２００４−２８６７３３号公報（特許文献１及び２）等に記載されている。
【０００６】
可視化流体の速度３成分を計測するための他の流速計測方法として、方形断面又は矩形断面を有する所定厚の帯状且つ立体的なレーザ光を可視化流体に照射し、異なる角度方向から可視化流体を撮影する少なくとも３台のＣＣＤカメラによって流体中の微粒子を撮像し、撮像により取得した粒子像の画像データに基づいて３次元空間内の粒子分布を再構築し、３次元空間における可視化流体の３次元速度ベクトル又は３速度成分を測定するトモグラフィックＰＩＶ計測法が知られている（特開２００５−９１３６４号公報、特開２０１３−２１７９０２号公報等）。トモグラフィックＰＩＶ計測法は、ＣＴ（Computed・Tomography)技術とＰＩＶ技術との融合により流体の３次元流速分布を効率的且つ高精度に測定することができる計測技法である。典型的には、４台のＣＣＤカメラが撮影に使用され、各カメラは、連続２時刻（ｔ及びｔ+Δｔ）の各々において流体中の微粒子を同時に撮影する。ステレオＰＩＶ技術が２次元計測技術であるのに対し、トモグラフィックＰＩＶ技術は、３次元計測技術である点で優位性がある。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特開２０１１−２４７６０１号公報
【特許文献２】特開２００４−２８６７３３号公報
【特許文献３】特開２００５−９１３６４号公報
【特許文献４】特開２０１３−２１７９０２号公報
【発明の開示】
【発明が解決しようとする課題】
【０００８】
しかしながら、トモグラフィックＰＩＶ計測法においては、計測システムの画像解析装置を構成する記憶装置は、３次元粒子分布のデータや、３次元voxel（ボクセル）群の輝度データ等の３次元データを解析データとして記憶し且つ保存すべく、過大な記憶容量の記憶部又は記憶手段を保有しなければならない。また、計測システムを構成する演算処理装置には、多量の３次元データのデータ解析のために多大な動作負荷が課せられるので、演算処理装置は、演算処理のために比較的長時間のデータ解析時間を要する。このため、トモグラフィックＰＩＶ計測法においては、高性能且つ大容量の電子機器を記憶装置及び演算処理装置として採用せざるを得ず、この結果、計測システムの初期設備費が高額化するとともに、システム構成が複雑化する傾向がある。
【０００９】
また、３次元粒子分布のデータや、３次元voxel群の輝度データ等の３次元データは、カメラによって撮影された実際の粒子画像とは相違するので、撮影された粒子画像の妥当性を計測過程において視覚的又は感覚的に評価し又は認識し難く、粒子画像に基づく３次元データの再構築の妥当性を視覚的又は感覚的に確認することもできない。
【００１０】
本発明は、このような課題に鑑みてなされたものであり、その目的とするところは、トモグラフィックＰＩＶ計測法により可視化流体の３次元速度ベクトル又は３速度成分を取得するための演算処理装置の負荷及びデータ解析時間を軽減又は短縮することができる流速計測方法及び流速計測システムを提供することにある。
【００１１】
本発明は更に、このような流速計測方法及び流速計測システムにおいて、記憶装置の記憶容量を低減することができる流速計測方法及び流速計測システムを提供することを目的とする。
【課題を解決するための手段】
【００１２】
本発明は、上記目的を達成すべく、微小粒子群により可視化された可視化流体を撮像する複数の撮像装置によって、微小時間間隔を隔てた連続２時刻の粒子像を撮像し、撮像により得られた粒子像の画像データに基づいて３次元空間内の粒子分布を再構築し、前記３次元空間における可視化流体の３次元速度ベクトル又は３速度成分を測定する可視化流体の流速計測方法において、
前記撮像装置によって撮像された２次元画像データに基づいて前記３次元空間内の各時刻の粒子分布を再構築し、連続２時刻の３次元輝度情報を得る工程(S3)と、
前記３次元空間を分割してなる所定厚のスライス領域を設定し、該スライス領域の内部の粒子分布を前記３次元輝度情報より抽出して、該粒子分布を複数の仮想投影面に投影することによって、連続２時刻の複数且つ仮想の２次元粒子画像を生成する工程(S4,S6)と、
複数の前記２次元粒子画像に２次元ＰＩＶを実行して得られる複数の画像面内の粒子の２次元移動ベクトルに基づいて、前記３次元空間における可視化流体の３次元速度ベクトル又は３速度成分を取得する工程(S7,S8,S11)とを有することを特徴とする流速計測方法を提供する。
【００１３】
本発明は又、微小粒子群により可視化された可視化流体を撮像する複数の撮像装置によって、微小時間間隔を隔てた連続２時刻の粒子像を撮像し、撮像により得られた粒子像の画像データに基づいて３次元空間内の粒子分布を再構築し、前記３次元空間における可視化流体の３次元速度ベクトル又は３速度成分を測定する可視化流体の流速計測方法において、
前記３次元空間を分割してなる所定厚のスライス領域を設定し、前記撮像装置によって撮像された２次元画像データに基づいて前記スライス領域内の各時刻の粒子分布を再構築し、前記スライス領域内における連続２時刻の３次元輝度情報を得る工程(S3',S5')と、
前記３次元輝度情報に基づいて前記スライス領域の内部の粒子分布を複数の仮想投影面に投影することによって、連続２時刻の複数且つ仮想の２次元粒子画像を生成する工程(S6)と、
複数の前記２次元粒子画像に２次元ＰＩＶを実行して得られる複数の画像面内の粒子の２次元移動ベクトルに基づいて、前記スライス領域における可視化流体の３次元速度ベクトル又は３速度成分を取得する工程(S7,S8)と、
各スライス領域における前記３次元速度ベクトル又は３速度成分の情報を集合、統合又は合成して、前記３次元空間における可視化流体の３次元速度ベクトル又は３速度成分を取得する工程(S11)とを有することを特徴とする流速計測方法を提供する。
好ましくは、前記３次元輝度情報を得る前記工程においては、前記２次元画像の画素に基づいて前記スライス領域内の前記３次元輝度情報を重み付け累算で特定するための重み係数が算出され且つ記憶され、前記２次元粒子画像を生成する前記工程においては、前記スライス領域内の粒子分布は、前記重み係数に従って前記仮想投影面に投影される。
【００１４】
本発明の上記構成によれば、粒子分布を再構築した３次元空間がスライス領域に分割され、スライス領域の粒子分布が３次元輝度情報より抽出される。スライス領域の粒子分布は、複数の仮想投影面に投影され、連続２時刻の複数且つ仮想の２次元粒子画像が生成される。複数の２次元粒子画像に対して２次元ＰＩＶが実行され、これにより得られる複数の画像面内の粒子の２次元移動ベクトルが求められ、かくして求められた２次元移動ベクトルに基づいて、３次元空間における可視化流体の３次元速度ベクトル又は３速度成分が得られる。即ち、演算処理装置は、３次元のvoxel群輝度データを演算処理するのではなく、２次元のpixel群輝度データを演算処理するにすぎないので、演算処理装置のデータ解析に要する負荷及びデータ処理時間は軽減又は短縮する。従って、上記構成の流速計測方法又は流速計測システムによれば、演算処理装置の解析時間を短縮し、計測を高速化することができる。
【００１５】
また、本発明の上記構成によれば、各スライス領域の粒子群を投影した２次元粒子画像は、シート光照明で撮影される通常のＰＩＶ粒子画像と直接に対比可能な粒子画像であるので、上記撮像装置による粒子撮影の妥当性を評価し得るとともに、３次元輝度情報の再構築の妥当性を視覚的に確認することができる。
【００１６】
更に、本発明の上記構成によれば、各スライス領域の粒子分布を仮想投影面に投影する過程において、ステレオＰＩＶ解析を最も精度良く行うことのできる視線方向を任意に選択することが可能であるので、測定精度の高精度化を達成することが可能となる。
【００１７】
他の観点より、本発明は、微小粒子群により可視化された可視化流体を撮像するように配置され、微小時間間隔を隔てた連続２時刻の粒子像を撮像する複数の撮像装置と、該撮像装置によって得られた画像データに基づいて３次元空間内の粒子分布を再構築し、前記３次元空間における可視化流体の３次元速度ベクトル又は３速度成分を演算する演算処理装置とを有する可視化流体の流速計測システムにおいて、
前記演算処理装置は、前記撮像装置によって撮像された２次元画像データに基づいて前記３次元空間内の各時刻の粒子分布を再構築し、連続２時刻の３次元輝度情報を得るように構成された３次元輝度情報取得手段(S3)と、
前記３次元空間を分割してなる所定厚のスライス領域を設定し、該スライス領域の内部の粒子分布を前記３次元輝度情報より抽出して、該粒子分布を複数の仮想投影面に投影することによって、連続２時刻の複数且つ仮想の２次元粒子画像を生成する２次元粒子画像生成手段(S4,S6)と、
複数の前記２次元粒子画像に２次元ＰＩＶを実行し、これにより得られる複数の画像面内の粒子の２次元移動ベクトルに基づいて、前記３次元空間における可視化流体の３次元速度ベクトル又は３速度成分を演算する演算手段(S7,S8,S11)とを有することを特徴とする流速計測システムを提供する。
【００１８】
本発明は又、微小粒子群により可視化された可視化流体を撮像するように配置され、微小時間間隔を隔てた連続２時刻の粒子像を撮像する複数の撮像装置と、該撮像装置によって得られた画像データに基づいて３次元空間内の粒子分布を再構築し、前記３次元空間における可視化流体の３次元速度ベクトル又は３速度成分を演算する演算処理装置とを有する可視化流体の流速計測システムにおいて、
前記演算処理装置は、前記３次元空間を分割してなる所定厚のスライス領域を設定し、前記撮像装置によって撮像された２次元画像データに基づいて前記スライス領域内の各時刻の粒子分布を再構築し、連続２時刻の３次元輝度情報を得るように構成された３次元輝度情報取得手段(S3',S5')と、
前記スライス領域の内部の粒子分布を複数の仮想投影面に投影することによって、該スライス領域内における連続２時刻の複数且つ仮想の２次元粒子画像を生成する２次元粒子画像生成手段(S6)と、
複数の前記２次元粒子画像に２次元ＰＩＶを実行し、これにより得られる複数の画像面内の粒子の２次元移動ベクトルに基づいて、前記スライス領域における可視化流体の３次元速度ベクトル又は３速度成分を演算する第１演算手段(S7,S8)と、
各スライス領域における前記３次元速度ベクトル又は３速度成分の情報を集合、統合又は合成して、前記３次元空間における可視化流体の３次元速度ベクトル又は３速度成分を取得する第２演算手段(S11)とを有することを特徴とする流速計測システムを提供する。
好ましくは、前記３次元輝度情報取得手段は、前記２次元画像の画素に基づいて前記スライス領域内の前記３次元輝度情報を重み付け累算で特定するための重み係数を算出し且つ記憶し、前記２次元粒子画像生成手段は、前記スライス領域内の粒子分布を前記重み係数に従って前記仮想投影面に投影する。
【００１９】
本発明の上記構成においては、３次元輝度情報の再構築、仮想の２次元粒子画像の生成、２次元ＰＩＶ、３次元速度ベクトル又は３速度成分の算出は、スライス領域単位に実行される。各スライス領域における３次元速度ベクトル又は３速度成分の情報は、集合、統合又は合成され、３次元空間における可視化流体の３次元速度ベクトル又は３速度成分が得られる。このような構成によれば、計測システムを構成する記憶装置は、各スライス領域の３次元輝度情報を記憶する記憶容量を保有すれば良く、従って、記憶装置の記憶容量を低減することができる。即ち、本発明の上記構成によれば、前述した解析時間の短縮、再構築の妥当性の視覚的確認および測定精度の高精度化という優位性のみならず、計測システムを構成する記憶装置の記憶容量を低減することが更に可能となるので、計測システムの初期設備費を更に低廉化することともに、システム構成を更に簡素化することができる。
【発明の効果】
【００２０】
本発明の流速計測方法及び流速計測システムによれば、トモグラフィックＰＩＶ計測法により可視化流体の３次元速度ベクトル又は３速度成分を取得するための演算処理装置の負荷及びデータ解析時間を軽減又は短縮することができる。
【００２１】
また、３次元輝度情報の再構築、２次元粒子画像の生成、２次元ＰＩＶ、３次元速度ベクトル又は３速度成分の算出をスライス領域単位に実行する本発明の流速計測方法及び流速計測システムによれば、演算処理装置の負荷及びデータ解析時間を軽減又は短縮し得るだけではなく、記憶装置の記憶容量を低減することができる。
【図面の簡単な説明】
【００２２】
【図１】図１は、本発明の好適な実施形態に係る流速計測システムの計測部の構成を概略的に示す概略斜視図である。
【図２】図２（Ａ）は、計測システムの初期設定において計測流路内に配置される校正板及び支柱の構造を示す斜視図であり、図２（Ｂ）は、校正板の基準点のパターンを示す平面図である。
【図３】図３（Ａ）は、図２に示す校正板の平面図であり、図３（Ｂ）は、校正板の底面図である。
【図４】図４（Ａ）は、ＣＣＤカメラによって連続２時刻に撮影された粒子画像であり、図４（Ｂ）は、ＣＣＤカメラによって撮影された粒子画像を拡大した画像である。
【図５】図５は、ＣＣＤカメラの撮影対象領域を拡大して示す斜視図である。
【図６】図６（Ａ）は、スライス領域内の粒子分布を仮想投影面に投影することによって得られた２次元粒子画像を例示する画像データであり、図６（Ｂ）は、複数の２次元粒子画像に基づいて得られた可視化流体の３次元速度ベクトルを例示する概略斜視図である。
【図７】図７は、本発明の好適な実施例に係る流速計測方法を示すフローチャートである。
【図８】図８は、本発明の好適な実施例に係る他の流速計測方法を示すフローチャートである。
【発明を実施するための形態】
【００２３】
本発明の好適な実施形態によれば、上記スライス領域は、上記３次元空間を等分割した空間である。好ましくは、スライス領域は、可視化流体に照射されるレーザー光の光軸と平行に可視化流体の流路を横断する空間である。
【００２４】
本発明の或る実施形態においては、上記仮想投影面は、撮像装置の光軸と直交する平面である。
【００２５】
好ましくは、可視化流体の流路は、円形の流路断面を有し、画像座標と物理座標との関係を求めるために校正画像を取得すべく流路内に配置される円形校正板が用いられる。校正板は、中心円形領域に配設された直交格子配列の指標と、外周領域に周方向に配設された環状配列の指標とを有する。このように外周領域の指標を流路の外周面の形態に相応するパターンに配置することにより、流路壁近傍における画像の歪みを適切に補正して画像座標と物理座標との関係を求め、計測領域の形状に最適な座標情報を効率良く取得することができる。
【００２６】
以下、添付図面を参照して、本発明の好適な実施形態について詳細に説明する。
【００２７】
図１は、本発明の好適な実施形態に係る流速計測システムの計測部の構成を概略的に示す概略斜視図である。
【００２８】
本発明に係る流速計測システムは、ＣＴ（Computed・Tomography)技術とステレオＰＩＶ技術との融合により流体の３次元流速分布を効率的且つ高精度に測定するトモグラフィック・ステレオＰＩＶシステム（以下、「ＴＳＰＩＶシステム」という。）である。ＴＳＰＩＶシステムは、微粒子（トレーサ粒子）を注入した可視化流体（気体）Ｐの可視化流路を形成する真円形断面の円筒体１０を有する。微粒子の寸法は、例えば、１０μｍ程度である。円筒体１０は、透明樹脂又は透明ガラス等の管壁１１を有し、円筒体１０の上部及び下部は、流体Ｐの供給路及び排出路（図示せず）に接続される。
【００２９】
ＴＳＰＩＶシステムは、レーザーヘッド及び光学系装置等を含む光源（図示せず）と、記憶装置、演算処理装置、画像表示（ディスプレイ）装置及び各種入・出力装置等を備えた画像処理装置（図示せず）とを有する。光源が照射した帯状且つ立体的なレーザー光Ｌが、円筒体１０の中心軸線と直交する方向に照射される。本例において、円筒体１０の中心軸線は、鉛直方向に配向され、レーザー光Ｌの光軸は、水平方向に配向される。レーザー光Ｌは、幅Ｗ及び高さＨの方形断面又は矩形断面を有し、レーザー光Ｌの幅Ｗは、円筒体１０の内径Ｄよりも大きく、円筒体１０内の垂直流路は、レーザー光Ｌの光軸と直交し、高さＨの範囲内においてレーザー光Ｌに包含される。レーザー光Ｌの高さＨは、数mm〜数cm程度の所定寸法、例えば、５〜１０mmの範囲内の所定値に設定される。可視化流体は、下向き又は上向き（本例では下向き）に円筒体１０内の流路を流動する。
【００３０】
ＣＣＤカメラ２、３がレーザー光Ｌの上側に配置され、ＣＣＤカメラ１、４がレーザー光Ｌの下側に配置される。ＣＣＤカメラ１〜４の光軸は、円筒体１０の中心軸線に対して所定角度をなして傾斜しており、ＣＣＤカメラ１〜４は、角度が異なる４方向から流体Ｐ中の微粒子を撮影するように位置決めされる。
【００３１】
図２（Ａ）は、円筒体１０内に配置される校正板２０及び支柱２１の構造を示す斜視図であり、図２（Ｂ）は、校正板２０の指標パターンを示す平面図である。また、図３（Ａ）は、校正板２０の平面図であり、図３（Ｂ）は、校正板２０の底面図である。
【００３２】
校正板２０は、ＣＣＤカメラ１〜４の画像座標と、流路内空間の物理座標との関係を予め特定するＣＣＤカメラ１〜４の校正を実施するためのものである。校正板２０は、ＣＣＤカメラ１〜４の位置固定後に円筒体１０内に挿入され、カメラ校正作業の完了後に撤去される。図２（Ａ）には、校正板２０及び支柱２１を組み付けた状態が示されている。校正板２０及び支柱２１は、黒色の表面を有する金属部材からなり、支柱２１の頂部は、校正板２０の切欠き部２２に嵌合し、校正板２０を水平に支持する。
【００３３】
図２（Ｂ）に示すように、校正板２０は、黒色の背景面に白色の指標（図３に示す基準点２３）を円形上面及び円形下面の所定位置に配置した構成を有する。背景の色と基準点の色とは輝度が大きく相違すれば良く、従って、白色の背景面に黒の円形基準点を配置することも可能である。図３に示すように、基準点２３は、校正板２０の各面（上面及び下面）に上下同一位置に配置される。基準点２３は、中心円形領域に配設された直交格子配列の基準点２３ａと、外周領域に周方向に配設された環状配列の基準点２３ｂとを有する。基準点２３ａ、２３ｂは、基本的に同一の直径を有する円形輪郭の指標であるが、ＸＹ方向を指示するために中央領域に配列された特定の基準点２３ｃは、基準点２３ａ、２３ｂよりも大きい直径を有する。
【００３４】
このように外周領域の基準点２３ｂを管壁１１の内周面の形態に相応するパターンに配置することにより、管壁１１を透過して得られる画像に歪みが生じないように画像座標及び物理座標の関係（投影関数）を適切に求めることができる。
【００３５】
校正板２０を用いて画像座標と物理座標との関係を求め、ＴＳＰＩＶシステムの初期設定が完了すると、微粒子を含む可視化流体Ｐが供給路（図示せず）から円筒体１０内に供給される。可視化流体Ｐは、円筒体１０内の流路を流動して排出路（図示せず）に流出する。
【００３６】
ＣＣＤカメラ１〜４の撮影は、微小時間間隔Δｔを隔てた２時刻（時刻ｔ、ｔ＋Δｔ）の各々において同時に実行される。図４（Ａ）は、ＣＣＤカメラ１〜４によって時刻ｔ及び時刻ｔ＋Δｔに撮影された粒子画像である。図４（Ｂ）は、ＣＣＤカメラ１によって時刻ｔに撮影された粒子画像を拡大した画像である。また、図５は、ＣＣＤカメラ１〜４の撮影対象領域αを拡大して示す斜視図である。
【００３７】
円筒体１０内の流路を移動する微粒子群は、撮影対象領域αにおいてレーザー光Ｌを反射し、図４に示すような粒子画像としてＣＣＤカメラ１〜４によって撮像される。粒子画像にMART（Multiplicative Algebraic Reconstruction Technique）を適用することにより、異なる時刻（ｔ及びｔ＋Δｔ）における３次元空間内の粒子分布を再構築することができる。再構築された３次元輝度情報より、voxel群（３次元）を要素データとした３次元配列が得られる。この再構築過程において、粒子画像（２次元）の各pixelの輝度は、その輝度の源泉となるvoxel群の輝度の重み付け加算により表され、重み付けの係数が算出される。
【００３８】
再構築された時刻ｔ及び時刻ｔ＋Δｔの３次元輝度情報に対し、任意の位置に定義された所定厚のスライス領域βの粒子分布が抽出される。図５には、厚さｈのスライス領域βが例示されている。スライス領域βは、レーザー光Ｌが照射された円筒体１０内の３次元空間（即ち、撮像対象領域α）をレーザー光Ｌの光軸と平行に等分割した円形平板状の空間であり、レーザー光Ｌの光軸と平行に可視化流体の流路を横断する。スライス領域βの厚さｈは、例えば、レーザー光Ｌの高さＨ×１／５〜高さＨ×１／１０の範囲内の所定値に設定され、各スライス領域βは、撮影対象領域αの全高さＨに亘って均等に分割することより設定される。
【００３９】
スライス領域βにおける粒子分布は、重み付け係数を用いた前述の計算処理の適用により、任意の複数の視線方向に定義された複数の仮想投影面に投影され、この結果、ステレオＰＩＶに最適な連続２時刻（ｔ及びｔ＋Δｔ）の複数の仮想粒子画像が生成される。本例において、仮想投影面は、ＣＣＤカメラ１〜４の光軸と直交する平面である。図６（Ａ）は、仮想投影面に投影された仮想粒子画像を例示する画像データである。
【００４０】
各スライス領域βに位置する粒子群の粒子画像として生成した仮想の投影画像に対し、２次元ＰＩＶが実行され、これにより得られた各画像面内の粒子の２次元移動ベクトルと、スライス領域βにおける３次元速度ベクトルと、別途に生成したカメラパラメータ及び２時刻の時間間隔Δtとの間に成立する関係式に関し、最小二乗法が適用され、スライス領域βにおける３次元速度ベクトルが算出される。このような過程を全てのスライス領域βに対して実行することにより、対象物理空間（撮影対象領域α）における可視化流体Ｐの３次元速度ベクトルを取得することができる。図６（Ｂ）は、上記過程により得られた撮影対象領域α内の３次元速度ベクトルを例示する概略斜視図である。
【００４１】
以上説明したとおり、このようなＴＳＰＩＶシステムによれば、voxel群（３次元）を要素データとした３次元のvoxel群データは、２次元ＰＩＶを適用可能な２次元の複数の粒子画像に変換される。即ち、ＴＳＰＩＶシステムにおいては、ＣＣＤカメラ１〜４によって撮影された画像データは、３次元のvoxel群データとして演算処理されるのではなく、２次元ＰＩＶを適用可能な２次元の複数の仮想粒子画像に変換され、連続２時刻（ｔ及びｔ＋Δｔ）の複数の仮想粒子画像（投影画像）に対して２成分ＰＩＶを実行することにより、粒子画像の移動ベクトル（２次元）が求められる。前述したように、この移動ベクトル、別途生成したカメラパラメータ、２時刻の時間間隔Δtおよび粒子群の３次元の速度ベクトルとの間に成立する関係式に関する最小二乗法の適用により、粒子群の３次元速度ベクトルが求まる。これを流速分布と見做すことによって、３次元流速分布が得られる。
【００４２】
このようなＴＳＰＩＶシステムによれば、演算処理装置が演算すべきデータは、３次元のvoxel群輝度データではなく、２次元のpixel群輝度データであるので、演算処理装置のデータ解析に要する負荷及びデータ処理時間は大幅に軽減又は短縮する（一桁〜二桁程度の軽減又は短縮である）。このため、演算処理装置の解析時間を短縮し、計測を高速化することができる。
【００４３】
また、各スライス領域βの粒子群を投影した粒子画像は、シート光照明で撮影される通常のＰＩＶ粒子画像と直接に対比可能な粒子画像であるので、粒子像の撮影の妥当性を評価し得るとともに、再構築の妥当性を視覚的に確認することができる。
【００４４】
更に、各スライス領域βの粒子群の投影において、ステレオＰＩＶ解析を最も精度良く行うことのできる視線方向を選択することができるので、測定精度の高精度化を達成することが可能となる。
【００４５】
次に、ＴＳＰＩＶシステムの画像解析過程を示すフローチャートに基づいて本発明の流速計測方法について説明する。
【００４６】
図７は、本発明の好適な実施例に係る流速計測方法を示すフローチャートである。
【００４７】
図７に示す如く、ＣＣＤカメラ１〜４は、異なる２時刻（時刻ｔ、ｔ＋Δｔ）において撮影対象領域α内の粒子を同時に撮影する（Ｓ（ステップ）１）。この撮影は、３台以上の撮像装置による流体中粒子の３方向且つ連続２時刻の撮影である。
【００４８】
各撮影時刻（時刻ｔ、ｔ＋Δｔ）においてＣＣＤカメラ１〜４により撮像された各粒子画像に関し、各pixelの輝度値が、その源泉となる３次元空間における粒子群の輝度の重み付け累算で表される下式の重み付け係数ｗｉｊを計算することにより求められる（Ｓ２）。
【数１】
2015161639.tif 000003
【００４９】
上式において、Ｉ（ｘｉ，ｙｉ）はＣＣＤカメラ１〜４によって撮像された２次元画像の２次元座標（ｘｉ，ｙｉ）において、各座標（ｘｉ，ｙｉ）のpixel（ピクセル）の輝度を表しており、ｉは、各２次元画像におけるpixelの番号である。Ｅ（Ｘｊ，Ｘｊ，Ｘｊ）は、３次元座標（Ｘｊ，Ｘｊ，Ｘｊ）を有するvoxel（ボクセル）の輝度を表しており、ｊはvoxel番号である。なお、voxelは、体積の要素であり、pixelが２次元画像データを表す正規格子単位であるのに対し、voxelは、３次元空間における正規格子単位である。輝度Ｅ（Ｘｊ，Ｘｊ，Ｘｊ）が相対的に高い値を示すことは、その座標（Ｘｊ，Ｘｊ，Ｘｊ）を有するj番目のvoxel jにレーザー光を反射する粒子群が比較的多量に存在することを意味する。
【００５０】
各撮影時刻（時刻ｔ、ｔ＋Δｔ）において、全ＣＣＤカメラ１〜４の粒子画像に関し、重み付け係数ｗijを用いてMART（Multiplicative Algebraic Reconstruction Technique）が適用され、各時刻における粒子の３次元輝度分布情報が再構築される（Ｓ３）。
【００５１】
撮影対象領域αは、任意の軸線方向（本例では、流路の中心軸線）に垂直な所定厚のスライス領域βに等分割され、再構築された２時刻（時刻ｔ、ｔ＋Δｔ）の３次元輝度分布は、スライス領域βの３次元輝度分布に分割される（Ｓ４）。
【００５２】
先ず、特定のスライス領域βが着目され（Ｓ５）、このスライス領域βに含まれる粒子の３次元輝度分布情報に基づき、pixelの輝度Ｉ（ｘｉ，ｙｉ）とvoxelの輝度Ｅ（Ｘｊ，Ｘｊ，Ｘｊ）との関係を示す上式を利用して、各ＣＣＤカメラ１〜４の画像を仮想投影面に投影した仮想的な２次元画像が生成される（Ｓ６）。
【００５３】
かくして生成された２時刻（時刻ｔ、ｔ＋Δｔ）の仮想的な２次元画像に関し、２次元ＰＩＶが適用され、各ＣＣＤカメラ１〜４の仮想投影面における粒子の２次元移動ベクトルが算出される（Ｓ７）。
【００５４】
次いで、この２次元移動ベクトルと、粒子の３次元速度ベクトルと、別途生成したカメラパラメータ及び２時刻の時間間隔Δtとの間に成立する関係式に関し、最小二乗法が適用され、粒子の３次元速度ベクトルが算出される（Ｓ８）。３次元速度ベクトルの算出は、撮影対象領域αの全スライス領域βに関して実行される（Ｓ９：Ｓ１０）。
【００５５】
全スライス領域βに関し、このような画像解析及び演算処理が完了した後、各スライス領域βに対して算出した３次元速度ベクトルが統合され、撮影対象空間αにおける粒子速度ベクトルが得られる（Ｓ１１）。所望により、ＣＣＤカメラ１〜４の撮影（Ｓ１）〜粒子速度ベクトル取得（Ｓ１１）の一連の流速計測過程が反復実行され（Ｓ１２、Ｓ１３）、これにより、撮影対象空間αの粒子速度ベクトルが時系列的に取得される。規定数の２時刻粒子速度データが得られた段階で一連の流速計測過程が終了する（Ｓ１２）。
【００５６】
図８は、本発明の好適な実施形態に係る他の流速計測方法を示すフローチャートである。
【００５７】
図７に示す流速計測方法は、撮影対象領域αの粒子分布を再構築して領域α全体の３次元輝度情報を得て、スライス領域β内の粒子分布を３次元輝度情報より抽出した後、仮想２次元粒子画像を生成して２次元ＰＩＶを実行するように構成されているが、図８に示す流速計測方法の如く、各スライス領域βの粒子分布だけを再構築して各スライス領域βの３次元輝度情報を得た後、仮想２次元粒子画像を生成して２次元ＰＩＶを実行するように流速計測方法を構成しても良い。
【００５８】
図８に示す如く、ＣＣＤカメラ１〜４による撮影対象領域α内の粒子の撮影（Ｓ１）により得られた各撮影時刻（時刻ｔ、ｔ＋Δｔ）の粒子画像に基づき、スライス領域βが設定される（Ｓ２'）。特定のスライス領域βが着目され（Ｓ３'）、前述の重み付け係数ｗij（上記数式１）が特定のスライス領域βの粒子群に関して計算され且つ記憶されるとともに（Ｓ４’）、各撮影時刻（時刻ｔ、ｔ＋Δｔ）において、全ＣＣＤカメラ１〜４の粒子画像より、重み付け係数ｗijを用いてMART（Multiplicative Algebraic Reconstruction Technique）が適用され、各時刻におけるスライス領域β内の粒子の３次元輝度分布情報が再構築され、記憶される（Ｓ５'）。
【００５９】
スライス領域βに含まれる粒子の３次元輝度分布情報に基づき、前述した仮想の２次元画像が重み付け係数ｗijに従って生成され（Ｓ６）、かくして得られた２時刻（時刻ｔ、ｔ＋Δｔ）の仮想２次元画像に関し、２次元ＰＩＶが適用され、各ＣＣＤカメラ１〜４の仮想投影面における粒子の２次元移動ベクトルが算出される（Ｓ７）。
【００６０】
次いで、この２次元移動ベクトルと、粒子の３次元速度ベクトルと、別途生成したカメラパラメータ及び２時刻の時間間隔Δtとの間に成立する関係式に関し、最小二乗法が適用され、粒子の３次元速度ベクトルが算出される（Ｓ８）。スライス領域β内の粒子群に関する重み付け係数ｗijの計算（Ｓ４’）〜３次元速度ベクトルの算出（Ｓ８）の過程は、撮影対象領域αの全スライス領域βについて実行される（Ｓ９：Ｓ１０）。
【００６１】
全スライス領域βに関し、このような画像解析及び演算処理が完了した後、各スライス領域βに対して算出した３次元速度ベクトルが統合され、撮影対象空間α全域の粒子速度ベクトルが得られる（Ｓ１１）。所望により、ＣＣＤカメラ１〜４の撮影（Ｓ１）〜粒子速度ベクトル取得（Ｓ１１）の一連の流速計測過程が反復実行され（Ｓ１２、Ｓ１３）、これにより、撮影対象空間αの粒子速度ベクトルが時系列的に取得される。規定数の２時刻粒子速度データが得られた段階で一連の流速計測過程が終了する（Ｓ１２）。
【００６２】
このような流速計測方法によれば、前述の実施形態と同じく、演算処理装置の解析時間の短縮、計測の高速化等の利点が得られる。しかも、上記流速計測方法においては、３次元輝度情報の再構築がスライス領域単位に実行されることから、画像解析過程において記憶すべき３次元輝度情報の情報量又はデータ数が低減するので、記憶装置の記憶容量を低減することができる。
【００６３】
以上、本発明の好適な実施形態について詳細に説明したが、本発明は上記実施形態に限定されるものではなく、特許請求の範囲に記載された本発明の範囲内で種々の変形又は変更が可能である。
【００６４】
例えば、上記実施形態においては、撮影対象空間は、真円形断面の円筒形流路であるが、多角形断面、楕円形断面、３角形断面等の各種断面形態の流路に設計しても良い。
【００６５】
また、上記実施形態では、粒子撮影のために４台の撮像装置（ＣＣＤカメラ）が使用されているが、粒子撮影のために３台又は５台以上の撮像装置を使用しても良い。
【産業上の利用可能性】
【００６６】
本発明は、微小粒子群により流体を可視化してなる可視化流体の流動場を撮像装置によって撮影し、可視化流体の３次元速度ベクトル又は３速度成分を計測する流速計測方法及び流速計測システムに適用される。殊に、本発明は、微小粒子群により可視化された可視化流体を撮像する複数の撮像装置によって、微小時間間隔を隔てた連続２時刻の粒子像を撮像し、撮像により得られた粒子像の画像データに基づいて３次元空間内の粒子分布を再構築し、３次元空間における可視化流体の３次元速度ベクトル又は３速度成分を測定する可視化流体の流速計測方法及び流速計測システムに適用される。本発明の流速計測方法及び流速計測システムは、トモグラフィックＰＩＶ計測法を改良したトモグラフィック・ステレオＰＩＶ計測法である。本発明のトモグラフィック・ステレオＰＩＶ計測法によれば、可視化流体の３次元速度ベクトル又は３速度成分を取得するための演算処理装置の負荷及びデータ解析時間を軽減又は短縮することができるので、本発明の実用的効果は、顕著である。
【符号の説明】
【００６７】
１〜４  ＣＣＤカメラ
１０  円筒体
１１  管壁
２０  校正板
Ｐ  可視化流体
Ｌ  レーザー光
α  撮影対象領域
β  スライス領域
【図１】
2015161639.tif 000004
【図２】
2015161639.tif 000005
【図３】
2015161639.tif 000006
【図４】
2015161639.tif 000007
【図５】
2015161639.tif 000008
【図６】
2015161639.tif 000009
【図７】
2015161639.tif 000010
【図８】
2015161639.tif 000011
【手続補正書】
【提出日】平成26年3月11日(2014.3.11)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】請求項６
【補正方法】変更
【補正の内容】
【請求項６】
前記３次元輝度情報を得る前記工程において、前記２次元画像の画素に基づいて前記スライス領域内の前記３次元輝度情報を特定するための重み係数を算出し且つ記憶し、
前記２次元粒子画像を生成する前記工程において、前記スライス領域内の粒子分布を前記重み係数に従って前記仮想投影面に投影することを特徴とする請求項２に記載の流速計測方法。
【手続補正２】
【補正対象書類名】特許請求の範囲
【補正対象項目名】請求項１２
【補正方法】変更
【補正の内容】
【請求項１２】
前記３次元輝度情報取得手段は、前記２次元画像の画素に基づいて前記スライス領域内の前記３次元輝度情報を特定するための重み係数を算出し且つ記憶し、
前記２次元粒子画像生成手段は、前記スライス領域内の粒子分布を前記重み係数に従って前記仮想投影面に投影することを特徴とする請求項８に記載の流速計測システム。
【手続補正３】
【補正対象書類名】明細書
【補正対象項目名】００１３
【補正方法】変更
【補正の内容】
【００１３】
本発明は又、微小粒子群により可視化された可視化流体を撮像する複数の撮像装置によって、微小時間間隔を隔てた連続２時刻の粒子像を撮像し、撮像により得られた粒子像の画像データに基づいて３次元空間内の粒子分布を再構築し、前記３次元空間における可視化流体の３次元速度ベクトル又は３速度成分を測定する可視化流体の流速計測方法において、
前記３次元空間を分割してなる所定厚のスライス領域を設定し、前記撮像装置によって撮像された２次元画像データに基づいて前記スライス領域内の各時刻の粒子分布を再構築し、前記スライス領域内における連続２時刻の３次元輝度情報を得る工程(S3',S5')と、
前記３次元輝度情報に基づいて前記スライス領域の内部の粒子分布を複数の仮想投影面に投影することによって、連続２時刻の複数且つ仮想の２次元粒子画像を生成する工程(S6)と、
複数の前記２次元粒子画像に２次元ＰＩＶを実行して得られる複数の画像面内の粒子の２次元移動ベクトルに基づいて、前記スライス領域における可視化流体の３次元速度ベクトル又は３速度成分を取得する工程(S7,S8)と、
各スライス領域における前記３次元速度ベクトル又は３速度成分の情報を集合、統合又は合成して、前記３次元空間における可視化流体の３次元速度ベクトル又は３速度成分を取得する工程(S11)とを有することを特徴とする流速計測方法を提供する。
好ましくは、前記３次元輝度情報を得る前記工程においては、前記２次元画像の画素に基づいて前記スライス領域内の前記３次元輝度情報を特定するための重み係数が算出され且つ記憶され、前記２次元粒子画像を生成する前記工程においては、前記スライス領域内の粒子分布は、前記重み係数に従って前記仮想投影面に投影される。
【手続補正４】
【補正対象書類名】明細書
【補正対象項目名】００１８
【補正方法】変更
【補正の内容】
【００１８】
本発明は又、微小粒子群により可視化された可視化流体を撮像するように配置され、微小時間間隔を隔てた連続２時刻の粒子像を撮像する複数の撮像装置と、該撮像装置によって得られた画像データに基づいて３次元空間内の粒子分布を再構築し、前記３次元空間における可視化流体の３次元速度ベクトル又は３速度成分を演算する演算処理装置とを有する可視化流体の流速計測システムにおいて、
前記演算処理装置は、前記３次元空間を分割してなる所定厚のスライス領域を設定し、前記撮像装置によって撮像された２次元画像データに基づいて前記スライス領域内の各時刻の粒子分布を再構築し、連続２時刻の３次元輝度情報を得るように構成された３次元輝度情報取得手段(S3',S5')と、
前記スライス領域の内部の粒子分布を複数の仮想投影面に投影することによって、該スライス領域内における連続２時刻の複数且つ仮想の２次元粒子画像を生成する２次元粒子画像生成手段(S6)と、
複数の前記２次元粒子画像に２次元ＰＩＶを実行し、これにより得られる複数の画像面内の粒子の２次元移動ベクトルに基づいて、前記スライス領域における可視化流体の３次元速度ベクトル又は３速度成分を演算する第１演算手段(S7,S8)と、
各スライス領域における前記３次元速度ベクトル又は３速度成分の情報を集合、統合又は合成して、前記３次元空間における可視化流体の３次元速度ベクトル又は３速度成分を取得する第２演算手段(S11)とを有することを特徴とする流速計測システムを提供する。
好ましくは、前記３次元輝度情報取得手段は、前記２次元画像の画素に基づいて前記スライス領域内の前記３次元輝度情報を特定するための重み係数を算出し且つ記憶し、前記２次元粒子画像生成手段は、前記スライス領域内の粒子分布を前記重み係数に従って前記仮想投影面に投影する。