ベッド負荷実験における画像ベースのラグランジュ粒子トラッキング

注：堆積物輸送実験は、ミラノ工科大学のレッコキャンパスにあるMountain Hydraulics Labの水路で行われた。水路は透明なアクリル材料で完全に構成されており、5.2 x 0.3 x 0.45 m 3です。チャンネルは2本のスチールビームで支えられており、ヒンジとスクリュージャックのために異なる傾斜で操作することができます。一連の蓋により、水路が閉じた導管として機能することが可能になる。これは、覆われた流れの構成であり、この作業に用いられるチャネルである。 1. Flume Slopeの測定と設定水路の出口をシールし、まだ水で満たしてください。 注：シーリング方法はチャネルの機能に依存します。ここに示したケースでは、チャンネルの終端部にカップリングフランジが装備されているため、ゴム製のワッシャを介してフランジにプラスチック製のプラグをボルト止めすることによってシールが得られました。 ネジジャックを操作する任意のチャネル勾配を設定する。 水が静かになるように少なくとも30分待ってください。振動の減衰と静水状態の確認を行うには、水面の壁にポイントゲージを置き、自由表面の高さを繰り返し測定します。 チャンネルの壁に複数の場所にポイントゲージを置き、フリーサーフェスの標高を読み取ります。 注：水面が水平である間に水路が傾いているので、測定値は互いに異なります。 ポイントゲージによる測定値の線形補間に基づいてスロープ値を計算します。 チャネルを支えるビームの1つに取り付けられ、床を指すレーザー距離センサーで測定を行います。 注記：ここでは、レーザーセンサーはチャンネルに恒久的に接続されているため、設定にプロトコルを使用する必要はありません。 様々な勾配条件に対して、1.2から1.6までのステップを繰り返します。 フィットする横軸にレーザ測定を行い、縦軸にチャネル勾配測定を用いて、デカルト平面内の実験点に対して線形方程式を適用する。レーザー測定からチャネル傾斜までの伝達関数を決定する（逆もまた同様）。 2.作業設定のセットアップチャネル底部の寸法と等しい総寸法を有する鋼板（2mm厚）のセットを準備する。 注：ここでは、1.3×0.3m 2のサイズを有する4つのプレートを使用した。 プレートの上に沈殿物の粒子を接着させることによって、固定された粗いベッドを作りなさい。これを行うには、プレートを2液ポリエステル樹脂で塗装し、その上に1〜1.5cmの土層を広げます。 注：ここで使用された沈殿物粒子は、白色で製造され、3mmの大きさのポリブチレンテレフタレート（PBT）粒子であった。上記の手順では、プレートに接着された沈降物の厚さは粒径に類似している。 W樹脂を乾燥させるために少なくとも24時間放置し、その後、プレートを傾けて余分な沈降物を除去し、沈降物を逃げさせる。防水スプレー塗料で黒い土砂の表面をペイントします。塗料を乾燥させるために少なくとも10時間待ちます。 土砂被覆プレート（ステップ2.1-2.3を参照）を水路に、ポリ塩化ビニル（PVC）サポートに置き、作業セクションの底を作ります。ベッドの連続性を確保するために連続するプレートの配置に注意してください。 注：PVCサポートの使用はオプションです。このケースでは、ゆるい堆積物の層で堆積物輸送実験を実施するために、本来、水深（幅0.3m、高さ0.45m、断面積比1未満）の断面が設計されていた。この原稿で提示された実験は固定ベッドで実行されていたため、幅と深さの比が大きすぎるという珍しい値になっていました。従って、PVC支持体は、この値のr値。 一連の透明なアクリル製の蓋（チャンネルの全長を覆う）をチャンネルの壁面に置き、カバーされた作業セクションを作成します。 注：ここに記載されているフロートでは、ふたは高さ20cmの内部ボックスで、単に水路の壁に敷設されています。したがって、実験中に蓋側に水が存在するが、水路内の流れを大きく変化させない。以下に示す実験に用いた作業部は0.3 MX 0.105メートル2と大きかったです。 3.定常流れ条件の確立ポンプのスイッチを入れ、流路に水を満たし、調整バルブを使用して流量を設定します。 注：この作業では、流量は、デリバリパイプに沿って配置された電磁流量計によって測定されました。 テール・ウォーター・レギュレーターを使用して、圧力ヘッドの高さをチャンネルの蓋よりわずかに上に設定し、カバーされたフローが存在することを確認します蓋に大きな浮力を与える。 注：ここに示されているチャンネルでは、水路の端にあるスティックのアレイによってテール・ウォーター・レギュレーションが実現されています。 フロー条件の安定性をチェックするために、流量と圧力ヘッドを繰り返し測定します。 4.フロー分布の特徴付け 様々な場所で流れ方向の速度成分の垂直プロファイルを測定する。 適切なプローブホルダーを使用して、水蓋の上に超音波速度プロファイラ（UVP）のプローブを置きます。チャネルの入口に向かってテールが付いた選択された傾斜にプローブを置きます。プローブチップと蓋の間の空間に適切な超音波カップリングジェルを適用し、超音波が空気中を通過するのを避けます。プローブを取得モジュールに接続します。 注：この場合、プローブホルダーはPVC製であり、傾斜したベースガイド付属。この支持は、プローブ角度が選択された後に構築された。 いくつかの瞬時速度プロファイルをサンプルします。 注：この実験で使用した機器では、このステップでは、放射された超音波の周波数、パルス繰返し周波数、分解能、および所望の瞬時プロファイル数を手動で設定する必要がありました。所望の数に達すると、プロファイルを連続的に取得して保存した。 ステップ4.1.1と4.1.2を繰り返します。ただし、チャンネルの出口に向かってプローブを尾部に置きます。 取得された流速の時間発展を目視検査することによって、信号脱スパイク24の必要性を評価する。必要に応じて、データ処理中にデスパイキングを実行します。 UVPの各測定位置（ゲート）の平均速度値を計算して、2つの測定値からのプローブ – ワイズ速度成分の時間平均プロファイルを得る（プローブの尾部をチャネルの入口と出口）。 任意の標高で測定された2つのプローブ方向速度の三角関数の構成を使用して、流れ方向および垂直方向の速度成分を得る。 ステップ4.1.2および4.1.3で測定した時間平均速度としてv upおよびv downを用いて、流れ方向（ u ）および垂直方向速度（ v ）成分を次のように決定する： ここで、αはチャネルに対するプローブの傾きである。 音響ビーム25を走行し、それを通して異なる媒体（ゲル、アクリル及び水）を占めプロファイルにおける測定位置に対する距離値を調整します。 すべての測定場所について、4.1.1〜4.1.7の手順を繰り返します。 <li> 流れ方向の速度成分の垂直プロファイルからせん断速度を測定する。 流れ方向速度成分のプロファイルが半対数プロットの線形傾向を示す隆起の範囲を決定する（ 図2参照）。 次のように対数方程式を当てはめることにより、測定されたプロファイルからせん断速度u sを推定する： ここで、 u （ z ）は床からの一定の高さzでの時間平均流速であり、 κは0.4に等しいカルマン定数であり、 z 0は流体力学的粗さの長さである。 せん断速度26の推定における不確実性を以下のように定量化する： ここでNはcに使用される値の数ですJは1からNの範囲のカウンタである。 5.土砂輸送実験の実施 目的のカメラパラメータ（解像度、周波数）を設定します。カメラの製造元のサポートを使用して、2つのアクションカメラを2つの流れ方向の位置でチャネルの底に面する蓋の側壁に取り付けます。 2つのカメラのフォーカスエリアが重なっていることを確認してください。 試行錯誤してカメラの位置と向きを調整します。 2つのフォーカスエリアが重ならないか、またはカメラビューがフルートとうまく整列していない場合、各カメラから短いビデオをキャプチャし、ビデオを見て、カメラの位置または向きを変更します。 注：現在の作業では、カメラは1920×1080ピクセルの解像度で30fpsで動作しました。 手順3.1〜3.3で説明したように、安定した流れ状態を確立します。白い粒子を供給しなさい帽子をステップ2.2および2.3で黒色に塗りつぶして塗りつぶした）を水路入口の流れに注入した。いくつかの粒子を選び、黒い床の上に白い粒子の濃度が低く維持されたまま（数秒ごとに粒子が放出されます）。実験の全期間中、飼育を続ける。 注：低粒子濃度が高い濃度18との状況と比較して17のトラッキングプロセスを単純化します。実際、その後の画像における同じ粒子と特定の画像内の1個の粒子との間の一致は、粒子19、27の元の位置の周りの探索ウィンドウに基づいています。濃度が高いほど、検索ウインドウ内の複数の粒子を検出する可能性が高くなり、ミスマッチが発生します。 ルームライトが消灯している場合は、カメラを同期させる必要があるため、室内照明をオンにします。カメラをトリガして撮影を開始する適切なコントロールを使用してください。両方のカメラが撮影を開始したことを確認した後、ルームライトを消灯します。 注記：この方法では、ライトがオフにされた時間は、画像内のピクセル強度の履歴からはっきりと分かるので、取得した2つのムービー間の時間シフトを（1フレーム未満の精度で）可能にします決断される。もちろん、視覚化は暗闇の中で行うことはできません。今回の実験では、自然照明のみを使用していました（部屋の側壁は完全にガラス製です）。実験中に照明条件が大きく変化した場合、照明は以下に説明する粒子の同定と追跡に影響するため、実験を繰り返す必要があります。 希望の時間（ここでは15分）撮影を続けてから、カメラを停止してください。 注記：代表的な結果は、以下に100秒間表示されます。 他の所望の水力学的条件について、5.2から5.7までのステップを繰り返す（for例えば、異なる流量）。 注：ここに記載されている実験の間、撮影後にベッドにトラップされた粒子がいくつか存在していました。それらは、別の構成がテストされる前に、蓋を持ち上げ、ブラシを使用して粒子を除去し、蓋を取り替えることによって除去する必要があります。 6.画像の前処理ピクセル座標に放射状の変換を適用し、画像を再マッピングすることによって、レンズの長さが限られていることによる画像の歪みを補正します。変換を適用するのに必要な較正係数をトライアルで調整します。 注：任意のピクセルから画像中心までの距離をrとすると、変換距離は次のように計算できます。 ここで、 kは、画像内に直線として現れるはずの底面の目視検査に基づいて、試験28によって調整する必要がある較正係数であるs。 水深の既知の距離と仰角に配置されたターゲットを使用して、ピクセルから実距離までの線形画像変換を決定する。 注：これは、ターゲットをクリックして、画素内の座標を取得するために、イメージファイルを開くコマンド関数ginputをコマンド関数imreadを使用して、（ 材料の表を参照）プログラミング言語、例えば達成との間の比率を求めることができますそのように決定されたピクセル座標と実際のピクセル座標を含む。 7.粒子の特定と追跡注記：両方のカメラで収集された画像に対して、次の操作はすべて実行する必要があります。粒子の同定および追跡は、Streams 29を使用して行った。このソフトウェアは、開発者に問い合わせると自由に入手できます。 Streamsはベッドロアのためのいくつかの実験で既に著者によって採用されました異なる条件におけるDの土砂運搬3、16、17、18、28、30。 画像→画像シーケンスの作成をクリックしてフレームをインポートします 。 2つのフレーム間の時間ステップと画像較正係数を入力します。シーケンスに含めるイメージファイルを選択します。 [ OK]をクリックします。 注：この場合、時間間隔は1/30秒（ステップ5.1で述べたフレームレートに相当）であり、画像較正係数は0.5mm /ピクセルであった。 パーティクルを特定する 白色粒子に対応するピクセルの典型的な強度値（0〜255のスケールで）を見つけるために、ランダムに選択されたいくつかの画像の強度マップを作成します。 作成された画像シーケンスを右クリックし、 Opを選択しますイメージビュー。 Shiftキーを押しながら 、任意の長方形をイメージ上に描画します。四角形を右クリックし、 ShowIntensityMatrixを選択します。適切なしきい値を選択します。 注：この場合、しきい値は80に設定されています。 バイナリ画像内の白色点の強度値および予想されるサイズに基づいて画像を閾値化する。 作成された画像シーケンスを右クリックし、[ プロセスビューを開く]を選択します。 Newをクリックし、 Filter images pipelineを選択してOKをクリックします。プロセスに名前を付け、シーケンス名を入力して[ OK ]をクリックします。 フィルタパイプラインをダブルクリックし、 新規をクリックし、 背景を削除を選択してOKをクリックします。 Processウィンドウで、 Newをクリックし、 Identify particlesを選択してOKをクリックします。 Single thresholdを選択し、 OKをクリックします。 アルゴリズムの中で </em>ページで、しきい値強度を入力し、スポットの直径を制限して、 OKをクリックします。作成したプロセスを選択し、[ 追加 ]をクリックしてパイプラインし 、次に[ 実行 ]をクリックします。 注：ここでのスポットサイズは0.5~8mmであった。スポットサイズは粒子サイズに関連するが、照明条件にも依存する。実際、スポットは粒子自体ではなく粒子によって反射された光に対応する。 パーティクルを追跡する 作成されたパーティクルレコードを右クリックし、[ 画像ビューを開く]を選択します。 順方向および逆方向の矢印ボタンをクリックして、フレームをスクロールします。表示されたスポット上にカーソルを移動して座標を読み取ることにより、連続する画像間の典型的な粒子の変位を観察する。それに応じて適切な検索ウィンドウを決定します。 作成されたParticleレコードを右クリックし、 Open process viewを選択します。 </em>。 [ 新規 ]をクリックし、 PTV分析パイプラインを選択して[ OK ]をクリックします。プロセスに名前を付け、[ OK ]をクリックします。 PTV分析パイプラインをダブルクリックし、 新規をクリックします。 原価計算ページで距離を選択します。 [ 最適化]ページで、ストリームの横と横の位置と検索ウィンドウのディメンションを入力し、[ OK ]をクリックします。 作成したプロセスを選択し、[ 追加 ]をクリックしてパイプラインし 、次に[ 実行 ]をクリックします。 注：パーティクルトラックは、理想的には、フォーカスエリアの最後または最終的な実験時間に到達する必要があります。同様に、フォーカスエリアの先頭または最初の時刻に開始する必要があります。しかしながら、測定されたトラックは、典型的には、2つの連続するフレーム間での粒子の不一致のために、強度が低いか、またはそれほど頻繁ではないため、一部のフレームで粒子が検出されないため、トラックがtiの境界に達することなく終了した場合再接続の候補を検索することができ、中断前の最後のフレームの終わりに近い位置と直後の位置でフレームを選択します。このようにして、パーティクルの単一トラックをステップ7.3.6で詳述したように再構成することができる。 測定されたトラックに中断がある場合は、その目的に合わせた検索ウィンドウを使用してトラックの再接続によって修正します。 作成されたParticleレコードを右クリックし、 Open process viewを選択します。 Newをクリックし、 Create Lagrangian path fieldを選択し、 OKをクリックします。 [ OK]をもう一度クリックします。プロセスをパイプラインに追加し、「 実行」をクリックします。 作成されたParticleレコード （2番目のもの）を右クリックし、 Open process viewを選択します。 Newをクリックし、 Join Lagrangian pathsを選択し、 OKをクリックします。 「 パラメータ 」ページで、検索ウィンドウの詳細を入力し、cOKを押してください。プロセスをパイプラインに追加し、「 実行」をクリックします。 作成されたパーティクルレコード （2番目のもの）を右クリックし、 テキストファイルへのパスを保存を選択します 。入力パスとファイル名を入力し、[ OK ]をクリックします。 8.異なるカメラからの軌跡に参加する注記：これは、測定領域のサイズを拡大するために複数のカメラを使用することを利用するために必要な操作です。手順は、著者によって開発されたグラフィカルユーザーインターフェイス（ 補足コードファイルを参照）でMatLabコード（ join_cameras.m ）によって実行されます 。 コンピュータのフォルダを参照し、両方のカメラのトラックファイルを検索し 、[ トラックプロパティを検索 ]をクリックします 。 2つのカメラの（ x 、 y ）参照を、下流のカメラからのデータへの座標変換を適用することによって統一します。決定する両方のカメラからの画像で目に見えるターゲットから2つの方向に沿って移動するために使用される定数。必要な値を入力し、 参照システムの統一をクリックします。 注記：この時点まで、画像の左下隅に（0,0）を持つカメラには独立参照が使用され、 x軸は右向き、 y軸は上向きです。この作業で使用された座標の並進は、それぞれ流れ方向および横方向に760.15および-1.5ピクセルであった。 2つのデータ間の重複領域を特定し、適切なテキストボックスに制限値を入力します。 注：この作品におけるオーバーラップの長さは、流れ方向に760.15ピクセルから880.11ピクセルまでで、チャンネル幅全体をカバーしていました。 分析を単純な状況に限定するために、重複領域の長さよりも短いすべての軌跡をサンプルから削除します（ 図1参照）。に関連するテキストボックスにスレッショルドの長さを入力し、[ 削除 ]をクリックします。 トラックに参加する マージテーブルをクリックして、2つのカメラからパーティクルトラックデータベースをマージします。ストリーム方向と横方向のオーバーラップ（ここでは10 mm）の許容値を入力し、[ JOIN ]をクリックします。 注記： JOINをクリックした後、次の操作を実行します。オーバーラップする領域で終わるトラックを見つけるまで、トラックをスクロールします。重複領域から開始して出てくるトラックの中から参加する候補を検索します。候補が見つかった場合は、2つのトラックが重複領域内にある時間を比較します。これらの時間間隔が一致する場合は、可能なすべての点で2つのトラックのパーティクル座標の差を計算します。 xとyの平均二乗差の平方根を取って2つの粒子軌跡は等しくなるようにする。このインジケータがしきい値よりも低い場合は、トラックに参加してください。オーバーラップした部分で、前の2つのトラックの平均に対応する新しいトラックを構築します。新しいジョイントができなくなるまで、これらの操作をすべて繰り返します。この作業で使用された閾値は、両方向で10mmであった。 結果ファイルに名前を付け、結合されたトラックを保存をクリックして結合されたトラックを保存します 。 図1.トラック結合の状況上流のカメラからのトラックは赤色であり、下流のカメラからのトラックは緑色である（簡略化のために1次元表現）。垂直の破線は、オーバーラップの名目上の長さに結びついています。トラックの中断の可能性があるため、様々な結果は、予想される直接的なtr第1のカメラからの重複領域に到達するトラックと、それを離れる第2のカメラからのトラックとを有するack（最初の4つのスケッチされたケースに対応する）。理論的に可能な13の状況の総数が提示される。解析を簡略化するために、重複領域の長さより短いトラックは予備データから除外されます。 この図の拡大版を見るには、ここをクリックしてください。 9.土砂輸送のキネマティクスの分析任意の時点で取られた（ x 、 y ）位置に関して表現された、運搬された各粒子のトラックを含む得られたデータセットから出発して、ベッド荷重粒子運動学3の描写についての様々な統計分析を行い、 4、5、8 </suP>、16、17、18。