エコー粒子画像流速測定法

要約

光学的に不透明な流体中または光学的に不透明なジオメトリを通じ、速度の2次元フィールドを取得することが可能なエコー粒子画像流速（EPIV）システムが記載されており、パイプフローの検証測定が報告されています。

要約

流体の流れの質量、運動量、エネルギーの輸送、最終的に流体速度場の時空間分布によって決定されます^。1したがって、理解のための前提条件、予測、および流体の流れを制御することにより、十分な空間と速度場を測定するための機能であり、時間分解能²速度測定のために光学的に不透明な流体中または光学的に不透明なジオメトリを通じ、エコー粒子画像流速（EPIV）は速度の"瞬間"二次元のフィールドを生成するための魅力的な診断手法です。これで^3,4,5,6紙、ハーゲン·ポアズイユのオペレーティングソフトウェア^8が記載されている商業用粒子画像流速測定法（PIV）を実行しているPCと、商業、医療用超音波マシン^7を統合^することによって構築EPIVシステムのプロトコル、および検証の測定（ すなわち 、層流管）流れが報告されている。

EPIVメジャーのメンツ、医療用超音波マシンに接続されているフェーズドアレイプローブは、異なる時間に圧電プローブ要素にパルスを与えることにより、二次元超音波画像を生成するために使用されます。各プローブ要素には、流体中に超音波パルスを送信し、流体（天然またはシードのいずれか）のトレーサー粒子はバック彼らが記録されているプローブに超音波エコーを反映しています。伝送に対する反射超音波とその時間遅延の振幅は、Bモード（輝度モード）は、2次元超音波画像と呼ばれるものを作成するために使用されます。具体的には、時間遅延は、流体中の散乱体の位置を決定するために使用されており、振幅は散乱体に強度を割り当てるために使用されます。シングルBモード画像tを得るために必要な時間は、それはフェイズドアレイ探触子のすべての要素をパルス化するのに要する時間によって決定される。複数のBモード画像を取得するため、秒あたりのフレーム数でシステムのフレームレート（fps）= 1 /＆ドLTA;トン。 （超音波イメージングのレビューのために9を参照）。

典型的なEPIV実験では、フレームレートがフローの状態に応じて、20〜60 fpsの間にあり、フロー内のトレーサー粒子の空間分布の100から1000 Bモード画像が取得される。一旦獲得されると、Bモード超音波画像は、PIV商用ソフトウェアを実行しているPCへのイーサネット接続を介して転送されます。 PIVソフトウェア、トレーサ粒子の変位場は、D（x、y）を [ピクセル]を使用することで、（ここで、xとyはそれぞれ、超音波画像の水平方向と垂直方向の空間的な位置を示す）連続超音波に相互相関アルゴリズムを適用することによって取得され、B-モード画像¹⁰速度場は、u（x、y）は [m / s]と、画像のペアの間の時間ステップを知っているので、変位フィールドから決定され、ΔT[s]とし、画像の倍率は、M [メートル/ピクセル]、 すなわち、U（x、y）は = MDの（x、y）/ΔT。タイムステップbetween画像ΔT= 1/fps + D（x、y）はA / B、Bは [ピクセル/ s]は、超音波探触子がイメージの幅を横切って掃引するのにかかる時間である。本研究では、M = 77〔μm/ pixel]で、FPS = 49.5 [1 / s]とし、B = 25047 [画素/秒]。一旦獲得されると、速度場は、関心のある流れの数量を計算するために分析することができます。

プロトコル

1。測定可能なフローの作成

1. EPIV検証測定はグリセリン水溶液（ - 50％水50％グリセリン）のパイプフローで実証されます。実験装置の概略を図1に示します。
2. 10μmの呼び径を有する中空のガラス球は百万分の約17重量部の濃度で流体に追加されます。中空ガラス球は、超音波造影剤として機能し、その大きさと密度が、彼 ​​らは受動的に流体の流れをたどるように選択される^。10
3. 固定電圧は、既知の流量を導入するポンプに供給されています。流量はUはパイプ内の平均速度であるU <<ΔX/ΔTは、Δxは EPIV測定体積の直線長さであり、ΔTは、画像間の時間ステップ、 すなわちであるように選択され、流れがある必要が超音波syのFPSに比べて"遅い"ステム^。3

2。超音波のキャリブレーションを

1. 外側パイプ壁に超音波探触子をマウントします。水ベースの局所ゲルは、プローブ面と管壁の間に超音波ビームの伝送損失を最小限にするために、超音波プローブに適用されます。
2. 超音波マシンをパワーオンします。超音波画像のライブストリームは、すべてのシステムの負荷と自動的に開始されます。
3. 超音波マシンのコントロールパネル上の奥行きコントロールノブを使って、イメージの深さを設定します。
4. 超音波マシンのコントロールパネルの2Dゲイン·ノブを使用して、画像全体のゲインを設定します。
5. 管壁からの散乱を減衰させるために、超音波信号の深関連の減衰を補償するために、 時間ゲイン補正（TGC）スライダを調整します。
6. 画像の幅、フォーカス、プローブ動作周波数 、 フレームレートが割り当て可能なコントロールノブを使って調整されます。これらのコントロールパネルの左上に位置する4つのノブは、システムが実行されているモードに応じて異なります。 2Dモード（現在使用されているように）で、左から右のノブは、それぞれ幅、フォーカス、周波数、フレームレートに対応しています。超音波イメージング⁹の基本的な原理に起因しが、これら4つのパラメータは、本質的に結合されていることに注意してください。したがって、与えられた超音波画像スキャン（ すなわち 、EPIV実験）のための空間と時間分解能はトレードオフの関係にある。
7. 10μmの中空ガラス球を播種したパイプ流の代表的な超音波像については図2を参照してください。限られた横方向の解像度が原因で、ガラス球が左右方向にまみれて、画像中の楕円のように表示されていることに注意してください。

3。データ収集

1. 新しい実験を開始するために超音波のコントロールパネルの[ 新規試験ボタンを押します。
2. 作る姓と名と患者IDのテスト番号で今日の日付でパイプ流を入力して、新しい"患者"。
3. "患者"の作成に続いて1000-1500の画像間の予め設定された最大に達するまで、超音波スキャンが開始され、新しいスキャンループが開始された後。超音波制御パネル上のフリーズボタンを 2回押すと、画像の最大プリセット数に到達する前にいつでもスキャンを再開します。
4. 超音波画像の良いセット（ すなわち 、シャープシード粒子画像と十分なシード粒子密度）を取得した後、画像の集録を停止するには、超音波のコントロールパネルのFREEZEボタンを押してください。
5. 超音波のコントロールパネルのシネループボタンを押します。セット内の最初の画像を選択するために超音波のコントロールパネル上の最初のサイクル·ノブを使用して分析される超音波画像のセット、および選択するための最終サイクルのノブを選択セット内の最後の画像。
6. 超音波画像の選択されたセットを保存するために超音波のコントロールパネルのイメージStoreボタンを押します。
7. 超音波のコントロールパネルの[アーカイブ]ボタンを押して、 終了試験を選択するには、マウスカーソルを使用します。これにより、ユーザーはローカルのハードドライブに保存する画像やcineloopsを選択するようにプロンプ​​トが表示されます。興味のあるシネ（複数可）を選択し、次に試験を終了します。
8. 超音波のコントロールパネルの[アーカイブ]ボタンを押して、最初の[詳細]を選択し 、[ ディスクの管理を選択するには、マウスカーソルを使用します。ディスク管理では、PIVソフトウェアを実行しているPCに保存されてシネ（s）を転送します。

4。ファイルタイプを変換する

1. 超音波画像は、超音波マシン上薬（DICOM）のファイル·タイプのデジタル·イメージング·コミュニケーションとして格納されます。開かおよびPIVソフトウェアによって読み取られるようにするために、DICOMファイルを画像ファイルに変換する必要があります。現在、Matlabのスクリプト実行DICOM2JPG.mは DICOMファイルジェイペグ（JPEG）ファイルタイプに変換するために使用されています。
2. JPEGの超音波画像は、その後LaVisionからDAVISソフトウェアを用いて分析する。

5。 D（x、y）はデービスを使用して、変位場を計算する

1. ダブルマウスは、PC上でデービスのアイコンをクリックしてください。 [新しいプロジェクト]を選択します。 PIVを選択します。
2. ツールバーのイメージのインポートを選択し、 番号付きファイルを経由してインポート]を選択します。プルダウンメニューで、JPEGの超音波画像が保存されているフォルダを探して、ダブルセットの最初の画像をクリックします。これは、このセット内のすべての番号が付けられた超音波画像をインポートします。
3. 典型的には、画像マスクを処理する超音波画像における関心領域（ROI）を単離するために定義されます。パイプ流の場合、マスクはパイプ壁（ すなわち 、流体）の間のROIを定義するために使用されます。
4. デイビスのメインコントロールパネルに移動し、インポートされた画像を含む現在のプロジェクトの下に配置]タブを選択し、タブのラベルの付いたバッチ処理を選択します 。これはインポートされた超音波画像のバッチ処理のためデイビスのベクトル処理ウィンドウを有効にします。
5. 操作リストから、PIV-時系列ツリーを使用して、 ベクトル計算のパラメータを選択して 、ベクトル処理に使用するパラメータを選択します。マスクを使用する場合は、ボックスのデータ範囲を確認する=ベクトル計算パラメータメニューでマスクされた領域を使用します 。 ベクトル演算パラメータの最適な選択が流動形状、流動特性、画像解像度、トレーサー粒子の密度、および所望の定量的なフロー分析に依存していることに注意してください^。10
   パイプ流測定では、通常、最良の結果が得られているパラメータは32×32ピクセルの²から8にxの減少尋問サイズでマルチアール50％のオーバーラップで8画素^2、。相対ベクトルの範囲の制限はすべて±（ウィンドウサイズ/ 2）と絶対値ベクトル範囲の制限が±5ピクセルに設定されていたに設定されていました。最後に、3×3pixel ²メディアンフィルタは、ノイズを抑制し、ベクトル場を滑らかにするために使用されていました。
6. バッチ処理の画面の左側に、処理すべき画像の合計金額を選択して、 起動処理を選択します 。これは、相互相関アルゴリズムを使用して連続した超音波画像との間で、A、D（x、y）を変位場を計算します。

6。ベクトル場の分析

1. 後処理とデータ解析のために、EPIVベクトルフィールドがtxtファイルとしてデイビスからエクスポートされます。これは、プロジェクト画面でJPEG画像の分岐の下ベクトル変位の分岐を選択することによって達成される。ツールバーで、[エクスポート]タブを選択し、ファイルタイプがASCIIを選択します。txtを選択して、/ exportフォルダを作成し、dは、[エクスポート]を選択します 。
2. エクスポートされたベクトル場はBxxxxx.txt、Bはバッファを表すと00001≤≤99999 xxxxxは、名前が付けられています。各ファイルには、4つのデータ列が含まれます。画像内のベクトルの（1）のx位置、画像内のベクトルの（2）のy位置を、（3）変位のx成分（ すなわち流れの変位）、（4）変位のy成分（ すなわち 、壁垂直変位）。 Bxxxxx.txtファイルは、最初の画像のペアの間の時間ステップを知ることは、ΔT[s]とし、画像の倍率は、M [メートル/ピクセル]、 すなわちU（X、yで、速度場を計算するために、MATLABで開かれ、処理され）= MDの（x、y）/ΔTは 、。画像間の時間ステップΔT= 1/fps + D（x、y）はA / B、Bは [ピクセル/ s]は、超音波探触子がイメージの幅を横切って掃引するのにかかる時間である。本研究では、M = 77〔μm/ pixel]で、FPS = 49.5 [1 / s]とし、B = 25047 [画素/秒]。次に、ensem関心のある他のフロー量のうちBLE平均速度ベクトル場は、平均速度の壁垂直プロファイルは、計算されます。 （セクション描写結果を参照してください。）

結果

瞬時エコー粒子画像流速（EPIV）ベクトル場を図3に示します。ベクトルプロットは、すべての4番目の列の速度ベクトルを示しており、背景色コンターマップは速度の大きさに対応しています。 1000年瞬時EPIVベクトルプロットにわたって平均アンサンブル平均ベクトルプロットを図4に示します。パイプ流と一致して、速度ベクトルは、流れ方向に主に、最大速度はパ?...

ディスカッション

光学的に不透明な流体中または光学的に不透明なジオメトリを通じ、速度の2次元フィールドを取得することが可能なエコー粒子画像流速（EPIV）システムの動作プロトコルが記述されていた。 EPIVの実用化は不透明な流体の流れは非常に多くのアプリケーションで発生した産業やバイオフローシステムの研究に適しています。ここに提示され、特定のシステムを意図的にリグノセルロース系エ...

資料

Name	Company	Catalog Number	Comments
試薬の名称	会社	カタログ番号	コメント（オプション）
超音波装置	GE	鮮やかな7 Proの
リニア超音波アレイ	GE	10リットル
DCの水ポンプ	KNF	NF 10 KPDC
ベクトル処理ソフトウェア	Lavision	デイビス7.2
処理ソフトウェアを投稿	Mathworks社	MATLAB 7.12
アクリルチューブ	マクマスター·カーR	8486K531
超音波ゲル	パーカー	Aquasonic 100