JP5980197B2 - 時系列流体速度計測システム - Google Patents

Description

本発明は、時系列流体速度計測システムに関する。

流体速度計測法として、流体内に混入した微小な粒子にレーザ光などの光を照射し、それらの散乱光などの画像を取得し、微小時間異なる粒子画像から粒子あるいは粒子群の移動距離を求め二次元平面内の流体速度（流速）を計測する粒子画像流速計（以下、ＰＩＶとよぶ）が知られている。

また、ＰＩＶによる速度場には非常に大きな誤差を含んだベクトル（以下、過誤ベクトルとよぶ）が存在する。過誤ベクトルはレーザライトシートに垂直方向の速度成分、渦度成分によって粒子がレーザライトシートから流出したり、レーザライトシートに流入してきたりすることによって相関が取れなくなることにより発生する。そのため、その発生は避けることができず、非常に大きな誤差を含むことになる。

これらの問題を解決するため、時間的に連続する２枚の画像から取得された各局所流速ベクトル分布に対して各速度定義点における時系列速度をフーリエ変換し、周波数空間においてカットオフ・フィルタ、又は低域フィルタを施し、逆フーリエ変換することにより、流体の最終的な時系列速度を得る技術が提案されている（特許文献１参照）。

特開２００４−２０３８５号公報

しかしながら、特許文献１の技術は、各速度定義点における時系列速度を一個の連続関数としてフーリエ変換等の処理を施しているため、誤差が大きい過誤ベクトル部分の誤差は小さくなるものの、誤差の小さい局所流速ベクトル部分の誤差はむしろ大きくなってしまうなど、流体の時系列速度取得の精度向上の観点から改良の余地があった。

そこで、本発明は、ＰＩＶにおいて、流体の時系列速度取得の精度向上を図り得るシステムを提供することを目的とする。

本発明の時系列流体速度計測システムは、
レーザライトシートを照射するレーザライトシート照射手段と、一又は複数の撮像素子を備える撮像装置と、前記撮像素子によって撮影された画像を解析する画像解析手段とを備える時系列流体速度計測システムであって、
前記画像解析手段は、前記レーザライトシート照射手段が測定対象の流体を連続的又は断続的に照射している間に前記撮像装置によって複数回にわたって撮影された前記流体の画像から、相互相関法により、各検査領域における前記流体の局所流速ベクトルの時系列変化態様を取得し、
１つの検査領域における前記局所流速ベクトルの時系列における基準時刻前及び基準時刻後の少なくとも一方の指標局所流速ベクトルのそれぞれに対し、他の指標局所流速ベクトルとの乖離度が閾値を超えているか否かを判定し、前記乖離度が前記閾値を超えている指標局所流速ベクトルが指標局所流速ベクトル全体に対して所定割合以上存在する場合、前記乖離度が前記閾値を超えている指標局所流速ベクトルが指標局所流速ベクトル全体に対して前記所定割合未満になるまで、指標局所流速ベクトルの数を増やし、
１つの検査領域における前記局所流速ベクトルの時系列において、基準時刻の対象局所流速ベクトルと前記指標局所流速ベクトルとの乖離度が前記閾値を超えた場合、当該対象局所流速ベクトルを前記指標局所流速ベクトルによって修正することにより、前記流体の各検査領域における局所流速ベクトルの時系列変化態様を得るように構成されていることを特徴とする。

本発明の時系列流体速度計測システムによれば、局所的な時系列変化態様から乖離した対象局所流速ベクトルが修正される。この結果、局所的な時系列変化態様と乖離した、すなわち過誤ベクトルである蓋然性が高い局所流速ベクトルが修正されやすくなる一方、局所的な時系列変化態様と乖離していない、すなわち正確である蓋然性が高い局所流速ベクトルは修正されにくくなる。
また、一の指標局所流速ベクトルと他の指標局所流速ベクトルとの前記乖離度が前記閾値を超えている場合、当該一の指標局所流速ベクトルは、過誤ベクトルである蓋然性が高いところ、本発明の時系列流体速度計測システムによれば、過誤ベクトルである蓋然性が高い指標局所流速ベクトルが所定割合以上存在する場合、当該指標局所流速ベクトルが所定割合未満になるまで指標局所流速ベクトルの数が増加される。これにより、正確である蓋然性が高い指標局所流速ベクトルの割合が増加するので、真値に近づける観点から対象局所流速ベクトルがより適切に修正される。

従って、当該構成の時系列流体速度計測システムにより、流体の時系列速度取得の精度向上が図られる。

本発明の時系列流体速度計測システムにおいて、前記画像解析手段は、前記各検査領域それぞれの局所流速ベクトルを算出する前に、前記画像の所定の範囲から、当該所定の範囲から定まる輝度を減算するように構成されていることが好ましい。

当該構成の時系列流体速度計測システムによれば、当該所定の範囲に含まれる輝度情報から定まる輝度が減算されることにより、各粒子間の境界及び背景光と粒子による散乱光との差が明確となり、当該各粒子を含む領域の相互相関関係が精度良く取得される。この結果、流体の時系列速度取得のさらなる精度向上が図られる。

本発明の時系列流体速度計測システムにおいて、前記画像解析手段は、前記各検査領域それぞれの局所流速ベクトルを算出するに当り、画像から分割された複数の重なり合った領域においてファストフーリエ変換処理を行う場合に、重なり合った部分のファストフーリエ変換処理の結果を再利用するように構成されていることが好ましい。

当該構成の時系列流体速度計測システムによれば、各領域の重なった部分のファストフーリエ変換処理（ＦＦＴ処理）の計算結果が再利用されることにより、計算量が大幅に削減される。この結果、複雑な処理の負荷が軽減されるため、より多くのデータが短時間で処理されうる。この結果、多くのデータ解析が必要な流体の時系列速度取得のさらなる精度向上が図られる。

本発明の時系列流体速度計測システムにおいて、前記閾値は、前記指標局所流速ベクトルに応じて定められる値であることが好ましい。

当該構成の時系列流体速度計測システムによれば、前記指標局所流速ベクトルに応じて閾値が定められることにより、あらかじめ閾値が定められている場合に比べ、対象局所流速ベクトルの乖離度が指標局所流速ベクトルの態様に鑑みて柔軟に判断される。この結果、流体の時系列速度取得のさらなる精度向上が図られる。

システム全体の概略を示す図。 カメラの露光期間とレーザ照射のタイミングとの関係を示す図。 流体の時系列的な速度を取得するためのフローチャート。 粒子の像径と相互相関の誤差との関係を示す図。 粒子画像補正処理を説明する図。 （ａ）は粒子画像補正処理が行われなかった場合の局所流速ベクトルの取得可能範囲を示す図、（ｂ）は粒子画像補正処理が行われた場合の局所流速ベクトルの取得可能範囲を示す図、。 第１の画像における各検査領域の説明図。（ａ）は左上隅角の第１検査領域についての説明図、（ｂ）はｘ軸方向へ１つずれた第１検査領域についての説明図、（ｃ）はｙ軸方向へ１つずれた第１検査領域についての説明図、（ｄ）は一般的な第１検査領域についての説明図。 第１の画像及び第２の画像における各検査領域の説明図。（ａ）は第２検査領域の説明図、（ｂ）は相関関数算出の説明図、（ｃ）は局所流速ベクトル算出の説明図。 各検査領域における輝度分布の説明及び、当該輝度分布に対するＦＦＴ処理の説明図。 各検査領域における輝度分布に対するＦＦＴ処理を省略できる部分の説明図。 時系列的誤差修正処理のフローチャート。 （ａ）は時系列的誤差修正処理における各流速ベクトルの修正内容の説明図、（ｂ）は時系列的誤差修正処理における窓長の変更についての説明図。 時系列的誤差修正処理における窓設定についての説明図。 （ａ）は、オリジナルデータ及び時系列的誤差修正処理前の流体の局所流速ベクトルの時系列的変化態様を示す図、（ｂ）は、オリジナルデータ及び従来の時系列的誤差修正処理後の流体の局所流速ベクトルの時系列的変化態様を示す図、（ｃ）は、オリジナルデータ及び本実施形態の時系列的誤差修正処理後の流体の局所流速ベクトルの時系列的変化態様を示す図。

（時系列流体速度計測システムの構成）
以下に本発明に係る平面及び空間の時系列流体速度計測システムを、図面を参照して説明する。図１は、本発明の平面及び空間の時系列流体速度計測システムの実施例の概略を示したものである。

本実施形態の時系列流体速度計測システムは、例えば、車両等の移動体の設計開発におけるシミュレーションの検証のため、移動中における当該移動体の周辺の流体の移動態様の計測を行う場合に使用され得る。

時系列流体速度計測システムは、テストセクション５０内を流れる被測定流体の移動態様を計測するためのシステムである。このため、時系列流体速度計測システムは、テストセクション５０内を流れる被測定流体を照射するためのレーザライトシート照射手段１０と、被測定流体に含まれるトレーサ粒子を撮像する撮像装置２０と、撮像装置２０によって撮像されたトレーサ粒子の粒子画像を解析することにより、被測定流体の移動態様を推定する画像解析手段３０と、レーザライトシート照射手段１０、撮像装置２０及び画像解析手段３０のそれぞれの動作を制御する制御装置４０とを備える。

レーザライトシート照射手段１０は、レーザ光源１１と、レーザライトシート生成光学系１２とを備える。

レーザ光源１１から発振したレーザ光は、レーザライトシート生成光学系１２によりシート状の光に変換されテストセクション５０内を流れる被測定流体に照射される。レーザライトシート生成光学系１２は、第１〜第２円筒面レンズ１２ａ〜１２ｂ及びアウトプットウィンドウ１２ｃから構成される。

撮像装置２０は、撮像装置２０によって撮像される面（撮像面）がテストセクション５０へ照射されたレーザライトシート照射面と平行になるように配置される。

テストセクション５０に照射されたレーザ光の散乱光の画像が撮像装置２０により断続的に撮影される。画像解析手段３０は、撮像装置２０により撮影された画像データを受信し、当該画像に対し、フーリエ変換、相互相関法、及びフィルタリングなどの演算処理を行う。なお、画像解析手段３０は、制御装置４０のハードウェアが所定のプログラムを読み込むことにより構成され得る。

制御装置４０は、レーザ光源１１の種類に応じた種々の同期方法により、レーザライトシート照射手段１０の照射タイミングと撮像装置２０の撮像タイミングとを同期させる。制御装置４０は、コンピュータにより構成されうる。制御装置４０は、図２に示されるように、撮像装置２０のカメラによる奇数枚目（２ｑ＋１枚目、２ｑ＋３枚目、ｑは整数）の露光期間の後半と偶数枚目（２ｑ＋２枚目、２ｑ＋４枚目）の露光期間の前半のそれぞれにおいてレーザーを断続的に照射させる。当該断続期間は、撮像画像における対象物の標準的な移動量が所定ピクセル数（例えば５ピクセル）程度になるように調節される。

テストセクション５０は、撮像装置２０のレンズからレーザライトシート照射面までの距離が所定の距離Ｚ_０となるように透明なアクリルで製作された断面矩形状のパイプを備える。予めトレーサ粒子が混入された流体がパイプ内を流れている。前記トレーサ粒子によるレーザライトシートの散乱光が撮像装置２０によって撮影されることにより、粒子画像が取得される。

本実施形態では、テストセクション５０内に存在するトレーサ粒子又はその粒子群を「粒子群」といい、粒子画像における個々の粒子の像を「粒子像」といい、粒子画像における個々の粒子群の像を「粒子群像」といい、テストセクション５０内における粒子群の位置を表す３次元の座標系（レーザライトシートが照射される平面を構成するＸ軸及びＹ軸並びに撮像装置２０のレンズからの距離Ｚ軸から成る座標系）を「実空間座標系」といい、粒子画像における粒子像等の位置を表す２次元の座標系（撮像の平面を構成するｘ軸及びｙ軸から成る座標系）を「画像座標系」という。

ここで、実空間座標系における座標（Ｘ，Ｙ，Ｚ）と画像座標系における座標（ｘ，ｙ）との関係は、撮像倍率Ｍを用いて次の数１〜３で表わされる。

また、撮影された画像の各ピクセルの位置を各ピクセルに１対１で対応したピクセル座標（ｉ，ｊ）で表わす。各ピクセルの大きさを基にｉ，ｊを上記ｘ，ｙに変換し、その後数１〜数２で変換することにより、各ピクセル（ｉ，ｊ）に対応する実空間上の位置を求めることが出来る。また、各ピクセルの輝度をＩ（ｉ，ｊ）で表わす。

（時系列流体速度計測システムの機能概略）
図３を参照しながら、時系列流体速度計測システムの機能の概略を説明する。

時系列流体速度計測システムは、撮像装置２０によって、テストセクション５０に照射されたレーザ光の粒子群による散乱光の画像を所定の間隔で複数回撮影する（図３／ＳＴＥＰ０１０）。ここで、所定の間隔とは、テストセクション５０における各粒子群の移動態様を観察できる間隔をいう。

時系列流体速度計測システムは、画像解析手段３０によって、撮影された全ての画像について、周辺の環境光の影響の除去処理（図３／ＳＴＥＰ０２１）及び粒子画像補正処理（図３／ＳＴＥＰ０２２）により、画像からノイズを除去する処理を行う（図３／ＳＴＥＰ０２０）。

周辺の環境光の影響の除去処理（図３／ＳＴＥＰ０２１）は、粒子を流す前にレーザー光を照射し、撮影しておいた画像を背景画像として各時刻の画像との差分をとる等の種々の公知の手法により周囲の環境光（例えば車体の反射光等）の影響を除去する処理である。

粒子画像補正処理（図３／ＳＴＥＰ０２２）は、各粒子群像の像径を修正するため、撮像画像のうちの各所定の大きさの範囲から最小の輝度を減算する処理である。

図４のグラフは、乱数により人工的に作った粒子画像を用いてＰＩＶ処理をし、その結果と正解値とを比較した場合の、粒子像の大きさ及び検査領域間の相関値の誤差の関係を示したグラフである。

図４のグラフは、ピクセルを単位として粒子像径の大きさを横軸とし、検査領域間の相関値の誤差を縦軸として、後述する局所流速ベクトル取得処理（図３／ＳＴＥＰ０３０）における検査窓の窓長ＷＬ（ｓ）（ＷＬ（ｓ１）＜ＷＬ（ｓ２）＜ＷＬ（ｓ３））ごとの３種類のグラフを示している。なお、一般に、局所流速ベクトル取得処理における検査窓の窓長は、窓長が大きいほど、検査領域間の相関値の誤差が小さくなる。

図４に示されるグラフからもわかるように、窓長ＷＬ（ｓ）がどのような大きさでも、粒子像の像径がおおよそＡピクセル（２画像による相互相関法に対しては、２ピクセル強）であるときに検査領域間の相関値の誤差が最小となる。この原因は、粒子像の像径が大きすぎると、一の粒子像のすそ野と他の粒子像のすそ野とが重なり合って、粒子像とすそ野との区別がつきにくくなるためと推察される。

粒子画像補正処理（図３／ＳＴＥＰ０２２）では、時系列流体速度計測システムは、図５に示されるように、各（Ａ＋１）×（Ａ＋１）ピクセル（例えば３×３ピクセル）の大きさの範囲に含まれる各ピクセルの輝度Ｉ（ｉ，ｊ）から、当該（Ａ＋１）×（Ａ＋１）ピクセルにおける最小の輝度Ｉ_ｍｉｎ（ｉ_１，ｊ_１）を減算する処理を行う。この結果、粒子像のすそ野の部分の明るさがゼロ輝度値（輝度基準値）とされるので、図５に示されるように、破線で示される当初の粒子像の大きさを、実線で示される粒子像の大きさ（Ａ×Ａピクセル程度）に補正することができる。この結果、検査領域間の相関値の誤差を最小とすることが出来、局所流速ベクトルの取得が向上する。

より具体的には、粒子画像補正処理が行われないと、図６（ａ）に示されるように、各粒子の輪郭がぼやけてしまうため、後述する局所流速ベクトル取得処理（図３／ＳＴＥＰ０３０）において局所流速ベクトルが取得される範囲が比較的狭い。これに対し、粒子画像補正処理が行われることにより、撮像された各粒子像の輪郭が明確になる。このため、図６（ｂ）に示されるように、局所流速ベクトル取得処理（図３／ＳＴＥＰ０３０）において局所流速ベクトルを取得できる範囲が粒子画像補正処理を行わなれない場合と比較して広くなる。

なお、最小輝度の減算に代えて、当該（Ａ＋１）×（Ａ＋１）ピクセルの大きさに含まれる輝度情報から定まる輝度（例えば、平均以下の輝度）が減算されてもよい。

ノイズ除去処理（図３／ＳＴＥＰ２０）の後、時系列流体速度計測システムは、画像解析手段３０によって、撮影された複数の画像から、第１時刻Ｔにおける第１の画像に含まれる第１検査領域と第２時刻Ｔ＋ΔＴにおける第２の画像に含まれる第２検査領域との相関関係に基づいて局所流速ベクトルｖ（ａ，ｂ，Ｔ）を取得する処理を行う（図３／ＳＴＥＰ０３０）。

時間ΔＴは微小な時間であり、第１の画像と第２の画像とは、例えば時刻的に連続する二枚の画像である。

ここで、図７（ａ）に示されるように、第１の画像に含まれる各ピクセルの輝度を、Ｉと上付き数字１とそのピクセルのｘ軸方向のピクセル座標ｉとｙ軸方向のピクセル座標ｊとを用いて、Ｉ^１（ｉ，ｊ）と表す。また、図７（ａ）に示されるように、第２の画像に含まれる各ピクセルの輝度を、Ｉと上付き数字２とそのピクセルのピクセル座標（ｉ，ｊ）とを用いて、Ｉ^２（ｉ，ｊ）と表す。また、第１の画像及び第２の画像の画像種別を区別しないときは、単にＩ（ｉ，ｊ）と表す。

第１の画像は、複数の第１検査領域Ｒ^１（ａ，ｂ）を含む。ここで、ａは、第１の画像左上の隅角を占める検査領域を基準としてｘ軸方向へどれだけずれているかを示す指標値であり、ｂは、第１の画像左上の隅角を占める検査領域を基準としてｙ軸方向へどれだけずれているかを示す指標値である。いずれの第１検査領域Ｒ^１（ａ，ｂ）も、図７（ａ）に太線で囲まれた範囲で示されるように、横（ｘ軸方向）ｍピクセル、縦（ｙ軸方向）ｎピクセルの長方形又は正方形である。

それぞれの第１検査領域Ｒ^１（ａ，ｂ）は、その内部または境界上に一つの第１速度定義点ｐ^１（ａ，ｂ）を有する。それぞれの第１速度定義点ｐ^１（ａ，ｂ）は、対応する第１検査領域Ｒ^１（ａ，ｂ）から見て相対的に同一の位置（例えば中心点）に存在している。当該第１速度定義点ｐ^１（ａ，ｂ）は、局所流速ベクトルｖ（ａ，ｂ，Ｔ）を算出するときの基準点の一つとなる。

一の速度定義点と当該一の速度定義点に隣接する速度定義点とは、ｘ軸方向（またはｙ軸方向）にｄ_ｘ（またはｄ_ｙ）ピクセルだけ離れている。

次に、各第１検査領域の輝度分布について図７（ａ）〜（ｄ）を参照しながら説明する。第１の画像の左上の隅角を占める第１検査領域Ｒ^１（１，１）の輝度分布は、当該第１検査領域Ｒ^１（１，１）に含まれるピクセルの輝度を用いて、以下の数４で表現される。

当該第１検査領域は、当該第１検査領域に対応する一つの第１速度定義点ｐ^１（１，１）を有する。

例えば、図７（ｂ）に示されるように、第１速度定義点ｐ^１（１，１）にｘ軸正方向で隣接する第１速度定義点ｐ^１（２，１）は、第１速度定義点ｐ^１（１，１）とｘ軸正方向にｄ_ｘピクセルだけ離れている。このため、第１速度定義点ｐ^１（２，１）に対応する第１検査領域の輝度分布は、次の数５で表現される。

また、図７（ｃ）に示されるように、第１速度定義点ｐ^１（１，１）にｙ軸正方向で隣接する第１速度定義点ｐ^１（１，２）は、第１速度定義点ｐ^１（１，１）とｙ軸正方向にｄ_ｙピクセルだけ離れている。このため、第１速度定義点ｐ^１（１，２）に対応する第１検査領域の輝度分布は、次の数６で表現される。

図７（ｄ）に示されるように、第１速度定義点ｐ^１（ａ_１，ｂ_１）は、第１速度定義点ｐ^１（１，１）とｘ軸正方向にｄ_ｙ（ａ_１−１）ピクセル、ｙ軸正方向にｄ_ｙ（ｂ_１−１）ピクセル離れている。このため、第１速度定義点ｐ^１（ａ_１，ｂ_１）に対応する第１検査領域の輝度分布は、次の数７で表現される。

第２の画像は、各第１検査領域Ｒ^１（ａ，ｂ）のそれぞれに対し、複数の第２検査領域Ｒ^２（ａ，ｂ，α，β）を含む。ここで、ａ及びｂは、当該第１検査領域Ｒ^１（ａ，ｂ）のａ及びｂと同一の指標値である。αは、図８（ａ）に示すように、対象の第１検査領域Ｒ^１（ａ，ｂ）と第２検査領域Ｒ^２（ａ，ｂ，α，β）とのｘ軸方向との距離（より正確には、第１速度定義点ｐ^１（ａ，ｂ）と、後述する第２速度定義点ｐ^２（ａ，ｂ，α，β）とのｘ軸方向の距離）をピクセル単位で表す指標値である。βは、図８（ａ）に示すように、対象の第１検査領域Ｒ^１（ａ，ｂ）と第２検査領域Ｒ^２（ａ，ｂ，α，β）とのｙ軸方向の距離（より正確には、第１速度定義点ｐ^１（ａ，ｂ）と、後述する第２速度定義点ｐ^２（ａ，ｂ，α，β）とのｙ軸方向の距離）をピクセル単位で表す指標値である。

なお、第１検査領域、第２検査領域の区別をつけないときは、Ｒと当該ａ，ｂと上付き数字の０とを用いて、検査領域Ｒ^０（ａ，ｂ）と表す。

いずれの第２検査領域Ｒ^２（ａ，ｂ，α，β）も、第１検査領域Ｒ^１（ａ，ｂ）と同一の形状及び大きさ、すなわち、横（ｘ軸方向）ｍピクセル、縦（ｙ軸方向）ｎピクセルの長方形又は正方形である。

それぞれの第２検査領域Ｒ^２（ａ，ｂ，α，β）は、その内部または境界上に一つの第２速度定義点ｐ^２（ａ，ｂ，α，β）を有する。それぞれの第２速度定義点ｐ^２（ａ，ｂ，α，β）は、対応する第２検査領域Ｒ^２（ａ，ｂ）から見て相対的に同一の位置（例えば中心点）に存在している。第１検査領域Ｒ^１（ａ，ｂ）と第１速度定義点ｐ^１（ａ，ｂ）との位置関係と、第２検査領域Ｒ^２（ａ，ｂ，α，β）と第２速度定義点ｐ^２（ａ，ｂ，α，β）との位置関係は同一である。当該第２速度定義点ｐ^２（ａ，ｂ，α，β）は、局所流速ベクトルｖ（ａ，ｂ，Ｔ）を算出するときの基準点の一つとなる。

第２検査領域ｐ^２（ａ，ｂ，α，β）は、図８（ａ）に太線で囲まれた範囲で示されるように、第２の画像内の横（ｘ軸方向）ｍピクセル、縦（ｙ軸方向）ｎピクセルの大きさの領域である。

第２検査領域Ｒ^２（ａ，ｂ，α，β）の輝度分布は以下の数８で表現される。

ＳＴＥＰ０３０の処理をより詳細に説明すると、時系列流体速度計測システムは、画像解析手段３０によって、検査領域の大きさｍ，ｎを決定する（図３／ＳＴＥＰ０３１）。

時系列流体速度計測システムは、処理開始時においては、ｍ、ｎを大きくし（例えば、ｍ＝６４、ｎ＝６４）、幅広い範囲でＳＴＥＰ０３２の相関処理を行う。一方、処理を進めた後には、時系列流体速度計測システムは、粒子のｘ方向の最大移動量ｘ＿ｍａｘ＝ｍａｘ（｜ｖ（ａ，ｂ，Ｔ）_ｘ｜）、ｙ方向の最大移動量ｙ＿ｍａｘ＝ｍａｘ（｜ｖ（ａ，ｂ，Ｔ）_ｙ｜）を算出し、以下の数９及び数１０を満たすｍ、ｎのうち、最小のｍ，ｎを設定する。換言すると、時系列流体速度計測システムは、粒子等の最大移動量を考慮して、移動後の粒子等が後述する探索領域ＳＲ（ａ，ｂ）（数１３参照）に含まれるとともに、ＳＲ（ａ，ｂ）の大きさが最小となるように、ｍ、ｎを設定する。

なお、ｘ，ｙ等の下付き文字は、以下の数１１及び数１２に示されるように、各成分ごとの値を表わす。

この結果、処理が進むほど、検査領域の大きさが粒子等の移動量に鑑みて最適な大きさになるので、不必要な計算が削減されうる。

時系列流体速度計測システムは、画像解析手段３０によって、第１の画像の第１検査領域の輝度分布Ｉ^１（ｉ＋ｄ（ａ−１），ｊ＋ｄ（ｂ−１））（ただし１≦ｉ≦ｍ、１≦ｊ≦ｎ）及び第２の画像の第２検査領域の輝度分布Ｉ^２（ｉ＋ｄ（ａ−１）＋α，ｊ＋ｄ（ｂ−１）＋β）（ただし１≦ｉ≦ｍ、１≦ｊ≦ｎ）との相関関数Ｃ（α，β）を、−ｍ／２＜α＜ｍ／２、−ｎ／２＜β＜ｎ／２の範囲で第１の画像に含まれる全ての第１検査領域Ｒ^１（ａ，ｂ）に対し算出する（図３／ＳＴＥＰ０３２）。

より具体的には、時系列流体速度計測システムは、図８（ｂ）に示されるように、ある第１検査領域Ｒ^１（ａ_２，ｂ_２）と、当該第１検査領域Ｒ^１（ａ_２，ｂ_２）に対応して定まる探索範囲ＳＲ（ａ_２，ｂ_２）に含まれる各第２検査領域Ｒ^２（ａ_２，ｂ_２，α_１，β_１）、Ｒ^２（ａ_２，ｂ_２，α_２，β_２）、Ｒ^２（ａ_２，ｂ_２，α_３，β_３）‥との相関値ｃ（α_１，β_１），ｃ（α_２，β_２），ｃ（α_３，β_３）‥を出力とする相関関数Ｃ（α，β）を求める。

ここで、探索範囲ＳＲ（ａ_２，ｂ_２）は、例えば、以下の数１３で表わされる範囲である。

ここで、相関関数Ｃ（α、β）は、例えば、以下の数１４で表わされる。

相関関数Ｃ（α，β）の算出に当たっては、第１検査領域Ｒ^１（ａ_２，ｂ_２）と第２検査領域Ｒ^２（ａ_２，ｂ_２，α，β）とに対し、それぞれ図９に示されるように、縦方向（ｙ方向）、横方向（ｘ方向）それぞれにＦＦＴ処理が行われる。そして、時系列流体速度計測システムは、ＦＦＴ処理によって取得された第１検査領域Ｒ^１（ａ_２，ｂ_２）の輝度分布特性と第２検査領域Ｒ^２（ａ_２，ｂ_２，α，β）の輝度分布特性を比較する。これにより、相関関数Ｃ（α，β）が算出される。

なお、隣接する検査領域の一部が重複しているので、ＦＦＴ処理に際し、重複している部分のＦＦＴ処理結果が再利用されることが好ましい。例えば、時系列流体速度計測システムは、図１０（ａ）に示されるように、検査領域Ｒ^０（ａ_３，ｂ_３）またはＲ^０（ａ_４，ｂ_４）と隣接する検査領域Ｒ^０（ａ_３＋１，ｂ_３）またはＲ^０（ａ_４，ｂ_４＋１））のＦＦＴ処理において、検査領域Ｒ^０（ａ_３，ｂ_３）またはＲ^０（ａ_４，ｂ_４）のＦＦＴ処理結果を再利用する。図１０（ｂ）に示されるように、検査領域Ｒ^０（ａ_３，ｂ_３）と検査領域Ｒ^０（ａ_３＋１，ｂ_３）とが重なっている部分については、検査領域Ｒ^０（ａ_３，ｂ_３）のＦＦＴ処理結果と検査領域Ｒ^０（ａ_３＋１，ｂ_３）のＦＦＴ処理結果とが同一となる。また、図１０（ｃ）に示されるように、検査領域Ｒ^０（ａ_４，ｂ_４）の処理結果のうちｘ方向のＦＦＴ処理結果の所定割合の部分も、検査領域Ｒ^０（ａ_４，ｂ_４＋１）において再利用され得る。このようなＦＦＴ処理結果の再利用により、ＦＦＴ処理結果の再利用がなされない場合に比べ、ほぼ以下の数１５であらわされる処理の演算量が削減され得る。

時系列流体速度計測システムは、以上のようにして求めた相関関数Ｃ（α、β）が最大値ｃ_ｍａｘ（α_０、β_０）をとるようなα_０、β_０の組を求める。第１検査領域Ｒ^１（ａ_２，ｂ_２）と相関関数Ｃ（α、β）が最大となる第２検査領域Ｒ^２（ａ_２，ｂ_２，α_０、β_０）は、当該第１検査領域Ｒ^１（ａ_２，ｂ_２）に含まれる粒子等像と同一又はほぼ同一の配置の粒子等像を含んでいると考えられる。

α_０＝α_３、β_０＝β_３である場合、図８（ｃ）に移動ベクトルＤ（ａ，ｂ）で示されるように、時刻Ｔにおいて第１検査領域Ｒ^１（ａ_２，ｂ_２）に対応する実空間座標系上に存在していた粒子等が、時刻Ｔ＋ΔＴまでに第２検査領域Ｒ^２（ａ_２，ｂ_２，α_３，α_３）に対応する実空間座標系上の位置まで移動したものと推察される。

このため、以下の数１６により当該第１検査領域Ｒ^１（ａ_２，ｂ_２）の時刻Ｔにおける局所流速ベクトルｖ（ａ_２，ｂ_２，Ｔ）が求められる。

ＳＴＥＰ０３２の処理が繰り返されることにより、各第１検査領域Ｒ^１（ａ，ｂ）の時刻Ｔにおける局所流速ベクトルｖ（ａ，ｂ，Ｔ）が取得される。

時系列流体速度計測システムは、画像解析手段３０によって、時刻Ｔの第１検査領域Ｒ^１（ａ，ｂ）における局所流速ベクトルｖ（ａ，ｂ，Ｔ）に基づき、当該第１検査領域Ｒ^１（ａ，ｂ）の周囲の第１検査領域Ｒ^１（ａ_±，ｂ_±）の局所流速ベクトルｖ（ａ_±，ｂ_±，Ｔ）に鑑みて、測定誤差と考えられる局所流速ベクトルｖ（ａ，ｂ，Ｔ）を補正する処理を行う（図３／ＳＴＥＰ０４０）。

より具体的には、時刻Ｔの第１検査領域Ｒ^１（ａ，ｂ）における局所流速ベクトルｖ（ａ，ｂ，Ｔ）の各成分ｖ（ａ，ｂ，Ｔ）_ｘｙが、周囲の第１検査領域Ｒ^１（ａ_±，ｂ_±）の局所流速ベクトルｖ（ａ_±，ｂ_±，Ｔ）の各成分ｖ（ａ，ｂ，Ｔ）_ｘｙから定まる条件（空間的誤差修正条件）を充足するか否かを判定する。空間的誤差修正条件としては、局所流速ベクトルｖ（ａ，ｂ，Ｔ）の各成分ｖ（ａ，ｂ，Ｔ）_ｘｙが、周囲の第１検査領域Ｒ^１（ａ_±，ｂ_±）の局所流速ベクトルｖ（ａ_±，ｂ_±，Ｔ）の各成分ｖ（ａ_±，ｂ_±，Ｔ）_ｘｙの中央値又は平均値から定まる所定範囲（例えば分散を用いて定められる範囲）に含まれないことという条件が採用され得る。

局所流速ベクトルｖ（ａ，ｂ，Ｔ）が空間的誤差修正条件を充足する場合、周囲の局所流速ベクトルの各成分ｖ（ａ_±，ｂ_±，Ｔ）_ｘｙにより、第１検査領域Ｒ^１（ａ_±，ｂ_±）の局所流速ベクトルの各成分ｖ（ａ，ｂ，Ｔ）_ｘｙを置き換える。例えば、当該周囲の局所流速ベクトルの各成分ｖ（ａ_±，ｂ_±，Ｔ）_ｘｙの重み付け平均が補正後の第１検査領域Ｒ^１（ａ_±，ｂ_±）の局所流速ベクトルの各成分ｖ（ａ，ｂ，Ｔ）_ｘｙとされる。

時系列流体速度計測システムは、画像解析手段３０によって、全ての画像に対してＳＴＥＰ０３０〜ＳＴＥＰ０４０の処理を行ったか否か（又は第１時刻Ｔ及び第２時刻Ｔ＋ΔＴの次の時刻（第３時刻及び第４時刻）の画像があるか否か）を判定する（図３／ＳＴＥＰ０５０）。当該判定結果が否定的である場合（図３／ＳＴＥＰ０５０‥ＮＯ）、時系列流体速度計測システムは、第３時刻を新たな第１時刻Ｔとして、ＳＴＥＰ０３０〜ＳＴＥＰ０４０の処理を続行する。

図３／ＳＴＥＰ０３０〜ＳＴＥＰ０５０までの処理によって、図１２（ａ）の実線に示されるように、時刻Ｔの第１検査領域Ｒ^１（ａ，ｂ）における局所流速ベクトルｖ（ａ，ｂ，Ｔ）が得られる。

図３／ＳＴＥＰ０５０の判定結果が肯定的である場合（図３／ＳＴＥＰ０５０‥ＹＥＳ）、時系列流体速度計測システムは、画像解析手段３０によって、時系列的誤差修正処理を行う（図３／ＳＴＥＰ０６０）。

時系列的誤差修正処理は、図１２（ａ）に実線で示される時刻Ｔの第１検査領域Ｒ^１（ａ，ｂ）における対象局所流速ベクトルｖ（ａ，ｂ，Ｔ）を、当該時刻Ｔの前後の時刻Ｔ_±の第１検査領域Ｒ^１（ａ，ｂ）における指標局所流速ベクトルｖ（ａ，ｂ，Ｔ_±）に鑑みて、図１２（ａ）に破線で示される時刻Ｔの第１検査領域Ｒ^１（ａ，ｂ）における対象局所流速ベクトルｖ_＋（ａ，ｂ，Ｔ）へ補正する処理である。

なお、Ｔ−ＷＬ（ｓ）／２≦Ｔ_±≦Ｔ＋ＷＬ（ｓ）／２、ＷＬ（ｓ）は窓長である。

以下においては、時系列流体速度計測システムは、局所流速ベクトルｖ（ａ，ｂ，Ｔ）の各軸ごと（例えば粒子群の局所流速ベクトルのｘ軸）の成分ｖ（ａ，ｂ，Ｔ）_ｘｙに対して処理を行う。

時系列流体速度計測システムは、基準時刻Ｔを初期化する（図１１／ＳＴＥＰ０６１）。

時系列流体速度計測システムは、基準時刻Ｔの窓を設定する（図１１／ＳＴＥＰ０６２）。

窓長ＷＬ（ｓ）が５であれば、図１３に示されるように、時刻Ｔの前後５つ（Ｔ＝ｔ_３の場合、時刻ｔ_１、ｔ_２、ｔ_３、ｔ_４、ｔ_５）の局所流速ベクトル（Ｔ＝ｔ_３の場合、ｖ（ａ，ｂ，ｔ_１）〜ｖ（ａ，ｂ，ｔ_５））が以下の処理で使用される。

次に、時系列流体速度計測システムは、条件判定のための各種のパラメータを定義する（図１１／ＳＴＥＰ０６３）。具体的には、次の数１７及び数１８によりパラメータが定義される。

次に、時系列流体速度計測システムは、時刻Ｔ、第１検査領域Ｒ^１（ａ，ｂ）における対象局所流速ベクトルｖ（ａ，ｂ，Ｔ）が時系列的誤差修正条件を満たすか否かを判定する（図１１／ＳＴＥＰ０６４）。より具体的には、以下の数１９の式が満たされているかが判定される。

なお、数１９の式の左辺が本発明の「乖離度」に該当し、数１９の右辺が本発明の「閾値」に該当する。このように、窓内の各局所流速ベクトルの成分値に応じて閾値が定められることにより、あらかじめ閾値が定められている場合に比べ、局所流速ベクトルの乖離度が時刻ｔの前後の時刻の局所流速ベクトルの態様に鑑みて柔軟に判断される。

乖離度及び閾値は、上記式に限られるものではなく、例えばｖ（ａ，ｂ，Ｔ）_ｘｙと、窓内の他の局所流速ベクトルの平均値等との差分を乖離度とし、あらかじめ実験的に定められた値を閾値としてもよい。

図１１／ＳＴＥＰ０６４の判定結果が否定的であった場合には、時系列流体速度計測システムは、時刻Ｔの対象局所流速ベクトルｖ（ａ，ｂ，Ｔ）（速度データ）を修正せずに次の処理へ進む。

例えば、図１２（ａ）及び図１３に示されるように、時刻ｔ_３における対象局所流速ベクトルｖ（ａ，ｂ，ｔ_３）が、当該時刻ｔ_３の前後ｔ_１，ｔ_２，ｔ_４，ｔ_５における指標局所流速ベクトルｖ（ａ，ｂ，ｔ_１）、ｖ（ａ，ｂ，ｔ_２）、ｖ（ａ，ｂ，ｔ_４）、ｖ（ａ，ｂ，ｔ_５）に鑑みて、数１９を充足していないと判断される場合、時系列流体速度計測システムは、局所流速ベクトルｖ（ａ，ｂ，ｔ_３）の値を変更せずに、次の時刻ｔ_４の処理へ進む。

図１１／ＳＴＥＰ０６４の判定結果が肯定的だった場合、時系列流体速度計測システムは、窓内の他の時刻の指標局所流速ベクトルｖ（ａ，ｂ、Ｔ_±）の各成分の重み付け平均（例えば、時刻Ｔに近い時刻Ｔ−１、Ｔ＋１の指標局所流速ベクトルの重み付け係数を高くする一方、時刻Ｔから離れた時刻Ｔ−２、Ｔ＋２の指標局所流速ベクトルの重み付け係数を低くする等）により、対象局所流速ベクトルｖ（ａ，ｂ，Ｔ）を補正する（図１１／ＳＴＥＰ０６５）。より具体的には、例えば、以下の数２０の式により、対象局所流速ベクトルｖ（ａ，ｂ，Ｔ）が修正される。

例えば、時刻ｔ_４における対象局所流速ベクトルｖ（ａ，ｂ，ｔ_４）が、当該時刻ｔ_４の前後ｔ_２，ｔ_３，ｔ_５，ｔ_６における指標局所流速ベクトルｖ（ａ，ｂ，ｔ_２）、ｖ（ａ，ｂ，ｔ_３）、ｖ（ａ，ｂ，ｔ_５）、ｖ（ａ，ｂ，ｔ_６）に鑑みて、数１９を充足していると判断される場合、時系列流体速度計測システムは、図１２（ａ）に矢印Ｑ１で示されるように、局所流速ベクトルｖ（ａ，ｂ，ｔ_４）を数２０によってｖ_＋（ａ，ｂ，ｔ_４）へ修正する。

また、時系列流体速度計測システムは、時刻ｔ_４以降の処理において、補正後の局所流速ベクトルｖ_＋（ａ，ｂ，ｔ_４）を用いる。例えば、時刻ｔ_５における局所流速ベクトルｖ（ａ，ｂ，ｔ_５）に対しては、局所流速ベクトルｖ（ａ，ｂ，ｔ_５）が、当該時刻ｔ_５の前後ｔ_３，ｔ_４，ｔ_６，ｔ_７における局所流速ベクトルｖ（ａ，ｂ，ｔ_３）、ｖ_＋（ａ，ｂ，ｔ_４）（補正後）、ｖ（ａ，ｂ，ｔ_６）、ｖ（ａ，ｂ，ｔ_７）に鑑みて、数１９を充足していると判断される場合、時系列流体速度計測システムは、図１２（ａ）に矢印Ｑ２で示されるように、局所流速ベクトルｖ（ａ，ｂ，ｔ_５）を数２０によってｖ_＋（ａ，ｂ，ｔ_５）へ修正する。

時系列流体速度計測システムは、図１１／ＳＴＥＰ０６５の処理の前に、指標局所流速ベクトルｖ（ａ，ｂ、Ｔ_±）の全体に対する指標局所流速ベクトルｖ（ａ，ｂ、Ｔ_±）が数１９で表わされる条件を充足する割合（以下、「数１９の充足率」という。）が所定割合（例えば４０％）以上か否かを判定してもよい。当該数１９の充足率が所定割合以上だった場合、時系列流体速度計測システムは、当該数１９の充足率が所定割合未満になるまで窓長ＷＬ（ｓ）を増加させた上で、図１１／ＳＴＥＰ０６３〜ＳＴＥＰ０６５の処理を行ってもよい。このようにすることで、図１２（ｂ）に示されるように、指標局所流速ベクトルの数１９の充足率が高い場合、すなわち、過誤ベクトルである蓋然性が高い指標局所流速ベクトルが多く存在する場合、窓長ＷＬ（ｓ）が増加されることにより、正確である蓋然性が高い指標局所流速ベクトルの割合を多くすることが出来る。この結果、真値に近づける観点から対象局所流速ベクトルがより適切に修正される。

時系列流体速度計測システムは、全ての時間Ｔについて処理を行ったか否かを判定する（図１１／ＳＴＥＰ０６６）。当該判定結果が否定的であった場合、時系列流体速度計測システムは、Ｔを１増加させ（図１１／ＳＴＥＰ０６７）、ＳＴＥＰ０６２〜ＳＴＥＰ０６５の処理を繰り返す。当該判定結果が肯定的だった場合、時系列流体速度計測システムは時系列的誤差修正処理を終了する。

なお、時刻Ｔにおける図１１／ＳＴＥＰ０６６の処理の前または後において、当該時刻Ｔを含む過去の一部又は全部の期間（例えば時刻Ｔ−γ〜時刻Ｔ（γ＞０））における局所流速ベクトルの修正率が所定値以上であるか否か（例えば当該期間の局所流速ベクトルの修正率が５％以上か否か）を判定し、当該判定結果が肯定的な結果であった場合、当該期間に対してＳＴＥＰ０６１〜ＳＴＥＰ０６７の処理を繰り返し実行するように時系列流体速度計測システムが構成されていてもよい。局所流速ベクトルの修正率が高い場合、元々のデータに測定誤差である蓋然性が高い局所流速ベクトル（過誤ベクトル）が多く含まれており、修正後にも過誤ベクトルが多く残存している蓋然性が高い。このような場合に、ＳＴＥＰ０６１〜ＳＴＥＰ０６７の処理を繰り返し実行することにより、過誤ベクトルの残存率を低下させることが出来る。

また、ある期間において局所流速ベクトルの修正率が所定値以下になっていても、当該期間より短い期間において局所流速ベクトルの修正率が高い、すなわち、局所的に過誤ベクトルが集中しているような場合も考えられる。

このことに鑑みて、当該所定値を当該期間の長さが短いほど高くなるように（例えば、当該所定値が、過去の第１期間（時刻Ｔ−γ１〜時刻Ｔ（γ１＞０））の場合には第１所定値（例えば５％）、第１期間よりも短い第２期間（時刻Ｔ−γ２〜時刻Ｔ（γ１＞γ２＞０））の場合には第１所定値よりも高い第２所定値（例えば１０％）となるように）時系列流体速度計測システムが構成されてもよい。当該第１期間は、過去の全部の期間であってもよい。

この結果、特に局所流速ベクトルの修正率が特に高い期間においてのみＳＴＥＰ０６１〜ＳＴＥＰ０６７の処理が繰り返し実行されるので、当該処理についての処理負荷の増加を抑えながら、局所的に過誤ベクトルが集中している蓋然性を低下させることが出来る。

時系列的誤差修正処理（図３／ＳＴＥＰ０６０）の後、時系列流体速度計測システムは、一連の処理を終了する。

（時系列流体速度計測システムの作用効果）
当該構成の時系列流体速度計測システムによれば、局所的な時系列変化態様から乖離した対象局所流速ベクトルｖ（ａ，ｂ，Ｔ）が修正される。この結果、局所的な時系列変化態様と乖離した、すなわち測定誤差である蓋然性が高い対象局所流速ベクトルｖ（ａ，ｂ，Ｔ）（過誤ベクトル）が修正されやすくなる一方、局所的な時系列変化態様と乖離していない、すなわち正確である蓋然性が高い対象局所流速ベクトルｖ（ａ，ｂ，Ｔ）は修正されにくくなる。

従って、当該構成の時系列流体速度計測システムにより、全体として流体の時系列速度取得の精度向上が図られる。

より具体的には、図１４（ａ）にはオリジナルデータ（点線）と時系列誤差修正処理前のデータ（実線）とが示され、図１４（ｂ）にはオリジナルデータ（点線）と特許文献１記載の技術による時系列誤差修正処理後のデータ（実線）とが示され、図１４（ｃ）にはオリジナルデータ（点線）と本実施形態による時系列誤差修正処理後のデータ（実線）とが示されている。

これらからもわかるように、特許文献１記載の技術による時系列誤差修正処理後のデータ（図１４（ｂ）参照）は、オリジナルデータから極端に外れたデータが減少しているものの、全体としてオリジナルデータからはむしろ乖離した結果となっている。一方、本実施形態による時系列誤差修正処理後のデータ（実線）は、誤差データが減少すると共に、全体としてオリジナルデータとの再現性も高まっている。これから、本実施形態の時系列流体速度計測システムにより、流体の時系列速度取得の精度向上が図られているといえる。

本実施形態の流体速度計測システムによれば、最小輝度が減算されることにより、各粒子像間の境界が明確となり、各局所流速ベクトルｖ（ａ，ｂ，Ｔ）が精度良く取得される。この結果、流体の時系列速度取得のさらなる精度向上が図られる。

本実施形態の流体速度計測システムによれば、各検査領域の重なった部分のＦＦＴの計算結果が再利用されることにより、計算量が大幅に削減される。この結果、複雑な処理しても画像解析手段３０等の演算処理装置に大きな負荷がからないため、より多くの処理を短時間で処理することが出来る。この結果、多くのデータ解析が必要な流体の時系列速度取得のさらなる精度向上が図られる。

本実施形態の時系列流体速度計測システムによれば、指標局所流速ベクトルｖ（ａ，ｂ，Ｔ_±）に応じて閾値ｖ＿ｒｍｓが定められることにより、あらかじめ閾値が定められている場合に比べ、対象局所流速ベクトルｖ（ａ，ｂ，Ｔ）の乖離度が指標局所流速ベクトルｖ（ａ，ｂ，Ｔ_±）の態様に鑑みて柔軟に判断される。この結果、流体の時系列速度取得のさらなる精度向上が図られる。

一の指標局所流速ベクトルｖ（ａ，ｂ，Ｔ_±）と他の指標局所流速ベクトルｖ（ａ，ｂ，Ｔ_±）との乖離度が閾値を超えている場合、当該一の指標局所流速ベクトルｖ（ａ，ｂ，Ｔ_±）は、過誤ベクトルである蓋然性が高い。

この点に着目して構成された本実施形態の時系列流体速度計測システムによれば、過誤ベクトルである蓋然性が高い指標局所流速ベクトルｖ（ａ，ｂ，Ｔ_±）が所定割合以上存在する場合、当該指標局所流速ベクトルｖ（ａ，ｂ，Ｔ_±）が所定割合以下になるまで指標局所流速ベクトルｖ（ａ，ｂ，Ｔ_±）の数が増加される。これにより、正確である蓋然性が高い指標局所流速ベクトルｖ（ａ，ｂ，Ｔ_±）の割合が増加するので、真値に近づける観点から対象局所流速ベクトルｖ（ａ，ｂ，Ｔ）がより適切に修正され、流体の時系列速度取得の精度向上が図られる。

１０ レーザライトシート照射手段
１１ レーザ光源
１２ａ 第１円筒面レンズ
１２ｂ 第２円筒面レンズ
１２ｃ アウトプットウィンドウ
１２ レーザライトシート生成光学系
２０ 撮像装置
３０ 画像解析手段
４０ 制御装置
５０ テストセクション
Ｒ^０（ａ，ｂ） 検査領域
Ｒ^１（ａ，ｂ） 第１検査領域
Ｒ^２（ａ，ｂ，α，β） 第２検査領域
ｖ（ａ，ｂ，Ｔ） 局所流速ベクトル

Claims (4)

1. レーザライトシートを照射するレーザライトシート照射手段と、一又は複数の撮像素子を備える撮像装置と、前記撮像素子によって撮影された画像を解析する画像解析手段とを備える時系列流体速度計測システムであって、
   前記画像解析手段は、前記レーザライトシート照射手段が測定対象の流体を連続的又は断続的に照射している間に前記撮像装置によって複数回にわたって撮影された前記流体の画像から、相互相関法により、各検査領域における前記流体の局所流速ベクトルの時系列変化態様を取得し、
   １つの検査領域における前記局所流速ベクトルの時系列における基準時刻前及び基準時刻後の少なくとも一方の指標局所流速ベクトルのそれぞれに対し、他の指標局所流速ベクトルとの乖離度が閾値を超えているか否かを判定し、前記乖離度が前記閾値を超えている指標局所流速ベクトルが指標局所流速ベクトル全体に対して所定割合以上存在する場合、前記乖離度が前記閾値を超えている指標局所流速ベクトルが指標局所流速ベクトル全体に対して前記所定割合未満になるまで、指標局所流速ベクトルの数を増やし、
   １つの検査領域における前記局所流速ベクトルの時系列において、基準時刻の対象局所流速ベクトルと前記指標局所流速ベクトルとの乖離度が前記閾値を超えた場合、当該対象局所流速ベクトルを前記指標局所流速ベクトルによって修正することにより、前記流体の各検査領域における局所流速ベクトルの時系列変化態様を得るように構成されていることを特徴とする時系列流体速度計測システム。
2. 請求項１記載の時系列流体速度計測システムにおいて、
   前記画像解析手段は、前記各検査領域それぞれの局所流速ベクトルを算出する前に、前記画像の所定の範囲から、当該所定の範囲から定まる輝度を減算するように構成されていることを特徴とする時系列流体速度計測システム。
3. 請求項１又は２記載の時系列流体速度計測システムにおいて、
   前記画像解析手段は、前記各検査領域それぞれの局所流速ベクトルを算出するに当り、画像から分割された複数の重なり合った領域においてファストフーリエ変換処理を行う場合に、重なり合った部分のファストフーリエ変換処理の結果を再利用するように構成されていることを特徴とする時系列流体速度計測システム。
4. 請求項１〜３のうちいずれか１項に記載の時系列流体速度計測システムにおいて、
   前記閾値は、前記指標局所流速ベクトルに応じて定められる値であることを特徴とする時系列流体速度計測システム。