JP2008058264A - 実河川を対象とした流速観測装置、流速観測方法,および流速観測プログラム - Google Patents
実河川を対象とした流速観測装置、流速観測方法,および流速観測プログラム Download PDFInfo
- Publication number
- JP2008058264A JP2008058264A JP2006238738A JP2006238738A JP2008058264A JP 2008058264 A JP2008058264 A JP 2008058264A JP 2006238738 A JP2006238738 A JP 2006238738A JP 2006238738 A JP2006238738 A JP 2006238738A JP 2008058264 A JP2008058264 A JP 2008058264A
- Authority
- JP
- Japan
- Prior art keywords
- angle
- flow velocity
- measuring
- camera
- river
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Abstract
【課題】実座標が既知である標定点を複数利用する変わりに、少なくとも2点の特徴点の角度と水面から高低差を計測することにより、幾何変換パラメータを決定して実河川の流速計測を行う。
【解決手段】実河川の表面水流を撮影するカメラ手段と、カメラ手段のレンズの投影中心座標からの角度として、少なくとも2点の特徴点の水平角および鉛直角を計測する角度計測手段と、少なくとも水面高さからカメラ手段のレンズの投影中心までの高低差を計測する距離計測手段と、撮影画像データと特徴点の水平角と鉛直角と水面高さからの高低差とを入力する解析手段とを備え、解析手段は、カメラ手段により撮影された実河川の表面水流の斜め画像を、角度計測手段および距離計測手段から決定される幾何変換パラメータを用いて垂直平面画像に幾何補正し、得られた垂直平面画像に対してPIV解析を実行し、実河川の水面各部の流速データを算定する。
【選択図】図2
【解決手段】実河川の表面水流を撮影するカメラ手段と、カメラ手段のレンズの投影中心座標からの角度として、少なくとも2点の特徴点の水平角および鉛直角を計測する角度計測手段と、少なくとも水面高さからカメラ手段のレンズの投影中心までの高低差を計測する距離計測手段と、撮影画像データと特徴点の水平角と鉛直角と水面高さからの高低差とを入力する解析手段とを備え、解析手段は、カメラ手段により撮影された実河川の表面水流の斜め画像を、角度計測手段および距離計測手段から決定される幾何変換パラメータを用いて垂直平面画像に幾何補正し、得られた垂直平面画像に対してPIV解析を実行し、実河川の水面各部の流速データを算定する。
【選択図】図2
Description
本発明は、実河川の平面流速分布や水量を計測する流速観測装置、流速観測方法,および流速観測プログラムに関するものである。
従来から、PIV(Particle Image Velocimetry)技術を応用して、実河川の表面流速を解析するシステムが知られている(特許文献1)。このPIV技術は、時間の経過とともに変化する流体の動画像を解析し、その移動量から流速計測する技術である。実河川の表面流速を解析するために、河岸からビデオカメラで実河川の表面水流を斜めに撮影し、この斜め画像を用い、パターンマッチングによる河川表面の濃淡分布の面的な移動ベクトルの追跡から、表面流速分布を求める。また、斜め画像のためにビデオカメラから近距離の地点と遠方の地点とでは俯角の差により画像の歪み度合いに差が出るので、幾何変換補正により俯角による影響を打ち消し、実河川の川面を垂直に見た平面画像(垂直平面画像)に幾何変換した後に、川幅の仮想的な横断線上の各区分についてPIV技術を用いて解析し、実河川の各区分の表面流速を算出していた。
実河川の表面水流の斜め画像から、川面を垂直に見た平面画像(垂直平面画像)に幾何変換するために必要となる幾何変換パラメータを決定するために、実空間座標上の三次元位置(例えば、x,y,z (m) )が既知である複数の点(標定点)を利用していた。
すなわち、幾何補正を行う際には、画像平面上の位置と三次元物理空間上の位置の対応付けを行う必要があり、その際には図1に示されるようなカメラモデルの座標系を用いた下記式の幾何変換式が用いられる。
ここで、aijは3軸まわりの回転角によって作られる3×3の回転行列の要素であり、下記式によって与えることができる。
ここで、ω,φ,κはそれぞれx軸,y軸,z軸方向の回転角である。また、式(1)中のΔx,Δyはレンズ歪みの補正係数であり、例えば、Δx=dr2x,Δy=dr2y(ここでr=(x2+y2)1/2)を用いることができる。
この幾何変換式のパラメータ(幾何変換パラメータ)を決定するために、特徴点のカメラ画像平面上の位置座標(x,y)(単位は pixel)と、その特徴点の実空間上の位置座標(X,Y,Z)のペアを複数与えるのである。
この幾何変換パラメータを決定することにより、画像情報を実際のスケールに変換するために、画像がどのような位置・角度により撮影されたかを把握することが可能となる。
現在のところ、この幾何変換パラメータを決定するために必要な手順である複数の標定点の測定を行う際には、測量用トータルステーション(例えば、特許文献2を参照)あるいはGPS(Global Positioning System)が利用されている。測量用トータルステーションの場合、標定点の測定を行う際にはプリズムを、GPSの場合、標定点の測定を行う際には計測器を測定点に配置する必要がある。
従来、標定点の測定を行う際には、トータルステーション(TS)あるいはGPSが利用し、プリズムあるいは計測器を測定点に配置する必要があるため、現地の状況によっては、撮影アングル中の適切な位置に標定点を配置することが、物理的な制約や安全上の問題から困難である状況がみられるといった問題があった。
上記問題点に鑑み、この実座標が既知である標定点を複数利用する変わりに、少なくとも2点の特徴点の角度情報と水面から高低差もしくは水面上の特徴点との距離を計測することにより、幾何変換パラメータを決定して実河川の流速計測を行う流速観測装置、流速観測方法,および流速観測プログラムを提供することを目的とする。
本発明は、この実座標が既知である標定点を複数利用する変わりに、画像の撮影位置にトータルステーションを据え付けて、少なくとも2点の特徴点の角度情報と水面から高低差もしくは水面上の特徴点との距離を計測して、この情報に基づき幾何変換パラメータを決定するものである。
上記目的を達成するため、本発明の請求項1の流速観測装置は、実河川の表面水流を撮影するカメラ手段と、前記カメラ手段のレンズの投影中心座標からの角度として、少なくとも2点の特徴点の水平角および鉛直角を計測する角度計測手段と、少なくとも水面高さから前記カメラ手段のレンズの投影中心までの高低差を計測する距離計測手段と、前記カメラ手段の撮影画像データと前記角度計測手段で得られた特徴点の水平角と鉛直角と前記距離計測手段で得られた水面高さからの高低差とを入力する解析手段とを備え、前記解析手段は、前記カメラ手段により撮影された実河川の表面水流の斜め画像を、前記角度計測手段および距離計測手段から決定される幾何変換パラメータを用いて垂直平面画像に幾何補正し、得られた垂直平面画像に対してPIV(Particle Image Velocimetry)解析を実行し、実河川の水面各部の流速データを算定する構成とされたことを特徴としたものである。
上記目的を達成するため、本発明の請求項1の流速観測装置は、実河川の表面水流を撮影するカメラ手段と、前記カメラ手段のレンズの投影中心座標からの角度として、少なくとも2点の特徴点の水平角および鉛直角を計測する角度計測手段と、少なくとも水面高さから前記カメラ手段のレンズの投影中心までの高低差を計測する距離計測手段と、前記カメラ手段の撮影画像データと前記角度計測手段で得られた特徴点の水平角と鉛直角と前記距離計測手段で得られた水面高さからの高低差とを入力する解析手段とを備え、前記解析手段は、前記カメラ手段により撮影された実河川の表面水流の斜め画像を、前記角度計測手段および距離計測手段から決定される幾何変換パラメータを用いて垂直平面画像に幾何補正し、得られた垂直平面画像に対してPIV(Particle Image Velocimetry)解析を実行し、実河川の水面各部の流速データを算定する構成とされたことを特徴としたものである。
上記の構成とすることにより、幾何変換パラメータを決定するために、複数の特徴点のカメラ画像平面上の位置座標(x,y)とその特徴点の実空間上の位置座標(X,Y,Z)の与えることが不要となる。すなわち、カメラ手段のレンズの投影中心座標からの角度として、少なくとも2点の特徴点の水平角および鉛直角を計測し、かつ、少なくとも水面高さからカメラ手段のレンズの投影中心までの高低差を計測することにより、幾何変換パラメータを決定することができるのである。
ここで、本発明のカメラ手段のレンズの投影中心座標からの角度として、少なくとも2点の特徴点の水平角および鉛直角を計測し、かつ、少なくとも水面高さからカメラ手段のレンズの投影中心までの高低差を計測することにより、幾何変換パラメータを決定することができる原理について、数式で示しながら以下に説明を行う。
幾何補正を行う際には、画像平面上の位置と三次元物理空間上の位置の対応付けを行う必要があり、その際には図1に示されるようなカメラモデルの座標系を用いた上記数式1の幾何変換式が用いられることは既に説明した。
本発明は、カメラ手段のレンズの投影中心を実座標の原点に設定することがポイントである。実座標の原点をカメラ手段のレンズの投影中心X0,Y0,Z0と設定することにより、X0,Y0,Z0は0となるため、上記数式1の幾何変換式は下記式のように簡略化することができる。また、下記式では説明を単純化するため、レンズ歪みの補正を割愛している。
本発明は、カメラ手段のレンズの投影中心を実座標の原点に設定することがポイントである。実座標の原点をカメラ手段のレンズの投影中心X0,Y0,Z0と設定することにより、X0,Y0,Z0は0となるため、上記数式1の幾何変換式は下記式のように簡略化することができる。また、下記式では説明を単純化するため、レンズ歪みの補正を割愛している。
上記数式3に示される変換式に含まれるパラメータは、ω,φ,k,cの4つのパラメータである。
これらの4つパラメータを決定するために、通常の幾何変換式のパラメータ決定手順と同様に、標定点を利用することを考える。例えば、標定点1の座標を画面座標(x1,y1) および実座標(X1,Y1,Z1)とする。この実座標は、単位ベクトル(n1x,n1y,n1z)および距離l1を用いて、(l1n1x,l1n1y,l1n1z)と表現することができる。これを上記数式3に代入すると、下記数式となる。
これらの4つパラメータを決定するために、通常の幾何変換式のパラメータ決定手順と同様に、標定点を利用することを考える。例えば、標定点1の座標を画面座標(x1,y1) および実座標(X1,Y1,Z1)とする。この実座標は、単位ベクトル(n1x,n1y,n1z)および距離l1を用いて、(l1n1x,l1n1y,l1n1z)と表現することができる。これを上記数式3に代入すると、下記数式となる。
ここで、上記数式4の中のl1は全て打ち消しあうことから、結局のところ下記数式の変換式になる。
つまり、座標系の原点をカメラ手段のレンズの投影中心とおいた場合には、標定点の実座標は、その標定点の方向を示す単位ベクトルのみが必要となり、距離は全く関係なくなるのである。ここで、実際問題として標定点の計測を行うことを考えると、レンズの投影中心に角度計などの角度計測手段の回転中心を一致させた上で、標定点の水平角φhおよ
び鉛直角φvを計測することになる。
標定点の単位ベクトルは、下記数式で特定されることになるのである。
び鉛直角φvを計測することになる。
標定点の単位ベクトルは、下記数式で特定されることになるのである。
本発明では、最低二点の標定点の角度を計測するだけで、幾何変換パラメータである4つのパラメータω, φ, κ, c を決定することが可能となるのである。
また、少なくとも水面高さからカメラ手段のレンズの投影中心までの高低差を計測するとは、水面に対して幾何補正を行う際には、レンズの投影中心に対する水面の位置を特定する必要があるので、最低一点は、角度だけでなく距離も含めて空間上の実座標(X,Y,Z)を特定する必要がある。水面を水平と仮定すればZのみ、すなわち、水面高さからの高低差のみでもよい。なお、この実座標を計測する点は標定点である必要はないため、カメラ手段の撮影アングル中で確認できない点でも構わない。
以上説明したように、本発明においては従来の標定点の計測が角度のみで済むため、計測点にミラーを据え付けることが物理的や安全面の理由で困難である点も標定点として利用できることとなる。また、距離が遠いため、レーザー距離計により距離計測ができない点も標定点として利用できることとなる。すなわち、カメラ手段の撮影アングル中に含まれる固定点とみなせる特徴点は、容易に標定点として用いることが可能となるのである。
次に、本発明の請求項2の流速観測装置は、実河川の表面水流を撮影するカメラ手段と、前記カメラ手段のレンズの投影中心座標からの角度として、少なくとも2点の特徴点の水平角および鉛直角を計測する角度計測手段と、河川の表面上の特徴点のうち少なくとも1点の距離を計測する距離計測手段と、前記カメラ手段の撮影画像データと前記角度計測手段で得られた特徴点の水平角と鉛直角と前記距離計測手段で得られたデータとを入力する解析手段とを備え、前記解析手段は、前記カメラ手段により撮影された実河川の表面水流の斜め画像を、前記角度計測手段および距離計測手段から決定される幾何変換パラメータを用いて垂直平面画像に幾何補正し、得られた垂直平面画像に対してPIV(Particle Image Velocimetry)解析を実行し、実河川の水面各部の流速データを算定する構成とされたことを特徴としたものである。
これは、請求項1の流速観測装置が、少なくとも水面高さから前記カメラ手段のレンズの投影中心までの高低差を計測するのに対して、請求項2の流速観測装置は、河川の表面上の特徴点のうち少なくとも1点の距離を計測するものである。河川の表面上の特徴点のうち少なくとも1点の水平角と鉛直角と距離とを計測することにより、河川の表面上の特徴点の実空間上の三次元位置座標が求まり、この三次元位置座標を用いて、カメラ手段のレンズの投影中心までの高低差を算出するのである。
次に、本発明の請求項3の流速観測装置は、上記請求項1又は2の流速観測装置において、角度計測手段と距離計測手段として、測距部と測角部とを一体としてなる測量用トータルステーションを用い、カメラ手段とトータルステーションとを一体化された構成とし、カメラ手段のレンズの投影中心座標からの角度として特徴点の水平角および鉛直角を計測することを特徴とするものである。
トータルステーションとは、距離と角度を同時に測定することにより三次元位置を特定する計測機であり、角度の測定は角度計により、距離の測定はレーザー距離計により計測する。レーザー距離計は、計測機から発射されたレーザーが計測点で反射されて計測機に戻ってくるまでの時間から算出している。一般的に、計測点でレーザーを反射するためにプリズムが利用されている。プリズムを利用しない計測も可能であるが、計測点の状況などに左右されるとともに計測距離の制約も大きい。
本発明の流速観測装置では、このトータルステーション、すなわち、光波測距儀と経緯儀とを一体としてなる測量用トータルステーションを用い、カメラ手段とトータルステーションとを一体化された構成とし、カメラ手段のレンズの投影中心座標からの角度として特徴点の水平角および鉛直角を計測することにしたものである。
これにより、本発明の流速観測装置が、容易に構築することが可能となる。
本発明の流速観測装置では、このトータルステーション、すなわち、光波測距儀と経緯儀とを一体としてなる測量用トータルステーションを用い、カメラ手段とトータルステーションとを一体化された構成とし、カメラ手段のレンズの投影中心座標からの角度として特徴点の水平角および鉛直角を計測することにしたものである。
これにより、本発明の流速観測装置が、容易に構築することが可能となる。
また、本発明の請求項4の流速観測装置は、上記請求項1〜3の流速観測装置において、有線もしくは無線の通信ネットワークを介して外部装置と周期的に及び/又はイベント的にデータ通信を行う手段を更に備える構成とされ、実河川の水面各部の流速データ及び/又は撮影画像データを周期的に及び/又はイベント的にネットワーク上に配信することを特徴とするものである。
上記の構成とされることにより、本流速観測装置がスタンドアロンの装置として機能するのではなく、実河川の水面各部の流速データや撮影画像データを所定の周期で外部装置(例えばデータ収集サーバなどのコンピュータ)に送信したり、外部装置からデータ送信要求があった場合に、ネットワーク上に配信することが可能となる。
また、本発明の請求項5の流速観測装置は、上記請求項1〜4の流速観測装置において、水平角と鉛直角の計測対象となる前記特徴点が撮影アングル中に存在しない場合において、河川の表面の見かけを部分的に変化させ得るレーザーポインタ手段を備える構成とされる。
撮影アングル中に特徴的な点がみられない場合は、レーザーポインタ等で、計測対象表面の見かけを部分的に変化させることにより、その点を画像上で確認することができ、角度計測を行う際のターゲットとすることが可能となるのである。
次に、本発明の請求項6の流速観測方法は、下記1)〜6)を備えた構成とされたことを特徴とする。
1)実河川の表面水流をカメラ手段によって撮影するステップ
2)カメラ手段のレンズの投影中心座標からの角度として、少なくとも2点の特徴点の水平角および鉛直角を計測するステップ
3)少なくとも水面高さからカメラ手段のレンズの投影中心までの高低差を計測するステップ
4)カメラ手段の撮影画像データと特徴点の水平角と鉛直角と水面高さからの高低差とを解析手段に入力するステップ
5)解析手段において、カメラ手段により撮影された実河川の表面水流の斜め画像を、前記水平角と鉛直角および高低差から決定される幾何変換パラメータを用いて垂直平面画像に幾何補正するステップ
6)解析手段において、得られた垂直平面画像に対してPIV(Particle Image Velocimetry)解析を実行し、実河川の水面各部の流速データを算定するステップ
1)実河川の表面水流をカメラ手段によって撮影するステップ
2)カメラ手段のレンズの投影中心座標からの角度として、少なくとも2点の特徴点の水平角および鉛直角を計測するステップ
3)少なくとも水面高さからカメラ手段のレンズの投影中心までの高低差を計測するステップ
4)カメラ手段の撮影画像データと特徴点の水平角と鉛直角と水面高さからの高低差とを解析手段に入力するステップ
5)解析手段において、カメラ手段により撮影された実河川の表面水流の斜め画像を、前記水平角と鉛直角および高低差から決定される幾何変換パラメータを用いて垂直平面画像に幾何補正するステップ
6)解析手段において、得られた垂直平面画像に対してPIV(Particle Image Velocimetry)解析を実行し、実河川の水面各部の流速データを算定するステップ
次に、本発明の請求項7の流速観測方法は、下記1)〜6)を備えた構成とされたことを特徴とする。
1)実河川の表面水流をカメラ手段によって撮影するステップ
2)カメラ手段のレンズの投影中心座標からの角度として、少なくとも2点の特徴点の水平角および鉛直角を計測するステップ
3)河川の表面上の特徴点のうち少なくとも1点の距離を計測するステップ
4)カメラ手段の撮影画像データと特徴点の水平角と鉛直角と距離とを解析手段に入力するステップ
5)解析手段において、カメラ手段により撮影された実河川の表面水流の斜め画像を、前記水平角と鉛直角および距離から決定される幾何変換パラメータを用いて垂直平面画像に幾何補正するステップ
6)解析手段において、得られた垂直平面画像に対してPIV(Particle Image Velocimetry)解析を実行し、実河川の水面各部の流速データを算定するステップ
1)実河川の表面水流をカメラ手段によって撮影するステップ
2)カメラ手段のレンズの投影中心座標からの角度として、少なくとも2点の特徴点の水平角および鉛直角を計測するステップ
3)河川の表面上の特徴点のうち少なくとも1点の距離を計測するステップ
4)カメラ手段の撮影画像データと特徴点の水平角と鉛直角と距離とを解析手段に入力するステップ
5)解析手段において、カメラ手段により撮影された実河川の表面水流の斜め画像を、前記水平角と鉛直角および距離から決定される幾何変換パラメータを用いて垂直平面画像に幾何補正するステップ
6)解析手段において、得られた垂直平面画像に対してPIV(Particle Image Velocimetry)解析を実行し、実河川の水面各部の流速データを算定するステップ
また、本発明の請求項8の流速観測方法は、上記請求項6又は7の流速観測方法において、測距部と測角部とを一体としてなる測量用トータルステーションとカメラ手段とが一体化された装置を用いて、実河川の表面水流を前記カメラ手段によって撮影するステップと、前記カメラ手段のレンズの投影中心座標からの角度として、少なくとも2点の特徴点の水平角および鉛直角を計測するステップと、少なくとも水面高さから前記カメラ手段のレンズの投影中心までの高低差を計測するステップを行うものである。
次に、本発明の請求項9の流速観測プログラムは、実河川の表面水流を撮影するカメラ手段と、前記カメラ手段のレンズの投影中心座標からの角度として、少なくとも2点の特徴点の水平角および鉛直角を計測する角度計測手段と、少なくとも水面高さから前記カメラ手段のレンズの投影中心までの高低差もしくは前記特徴点との距離を計測する距離計測手段と、実河川の表面各部の流速分布を算定する解析コンピュータとを備えた流速観測システムにおいて、解析コンピュータに、前記カメラ手段により撮影された実河川の表面水流の斜め画像を前記角度計測手段および距離計測手段から決定される幾何変換パラメータを用いて垂直平面画像に幾何補正させる手順と、得られた垂直平面画像に対してPIV(Particle Image Velocimetry)解析を実行させ、実河川の水面各部の流速データを算定させる手順とを実行させるものである。
また、本発明の請求項10の流速観測プログラムは、上記請求項9の流速観測プログラムにおいて、有線もしくは無線の通信ネットワークを介して外部装置と周期的に及び/又はイベント的にデータ通信を行う手順を更に実行させるものである。
また、本発明の請求項11のコンピュータ読み取り可能な記録媒体は、上記請求項9又は10の流速観測プログラムを記録したものである。
本発明は、トータルステーション等による角度計測では計測点にプリズムを配置する必要がないため、計測点に近づくことが困難である地点を幾何変換パラメータの決定に利用することが可能となる。これにより解析のための画像撮影をより迅速かつ安全に行うことが可能となり、計測アングルの自由度も高まるといった効果を有する。
以下、本発明の実施例について、図面を参照しながら詳細に説明していく。
図2は、本発明の流速観測装置の概念図を示している。また、図3は、本発明の流速観測装置の機能ブロック図を示している。また、図4は、本発明の流速観測装置の処理フローを示している。
本発明の流速観測装置は、図2に示されるように、実河川の水流表面を広範囲に撮影するために、川岸の土手4に設置される。本発明の流速観測装置1は、カメラ手段のレンズの投影中心座標からの角度として、少なくとも2点の特徴点(例えば、図2中のA,Bの2点)の水平角および鉛直角を計測する角度計などの角度計測手段を備える。また、少なくとも水面高さからカメラ手段のレンズの投影中心までの高低差(図2中のhで示す)を計測するレーザー測定器などの距離計測手段を備える。
図2は、本発明の流速観測装置の概念図を示している。また、図3は、本発明の流速観測装置の機能ブロック図を示している。また、図4は、本発明の流速観測装置の処理フローを示している。
本発明の流速観測装置は、図2に示されるように、実河川の水流表面を広範囲に撮影するために、川岸の土手4に設置される。本発明の流速観測装置1は、カメラ手段のレンズの投影中心座標からの角度として、少なくとも2点の特徴点(例えば、図2中のA,Bの2点)の水平角および鉛直角を計測する角度計などの角度計測手段を備える。また、少なくとも水面高さからカメラ手段のレンズの投影中心までの高低差(図2中のhで示す)を計測するレーザー測定器などの距離計測手段を備える。
図3の機能ブロック図に示されるように、カメラ手段21の撮影画像データ(斜め画像データ)と角度計測手段22で得られた特徴点の水平角と鉛直角と距離計測手段23で得られた水面高さからの高低差とを解析手段24へ入力する。解析手段24はコンピュータで実現される。解析手段24では、角度計測手段22および距離計測手段23から幾何変換パラメータを決定し(25)、カメラ手段21により撮影された実河川の表面水流の斜め画像を、決定された幾何変換パラメータを用いて垂直平面画像に幾何補正し(26)、得られた垂直平面画像に対してPIV解析を実行し(27)、実河川の水面各部の流速データを算定する(28)。
また、本発明の流速観測装置は、特徴点の角度計測と水面からの高低差を計測することで、幾何変換パラメータを決定するため、従来の流速観測装置の計測手順によるような、標定パネルの設置・測定や見通しを十分に確保するための草木の伐採を不要とする。
図4の処理フローに示されるように、計測地点の決定(S10)、撮影アングルの決定(S20)、特徴点の測量(S30)、幾何変換パラメータの決定(S40)を行い、カメラ手段による撮影(S50)、撮影画像の幾何補正(S60)、PIV解析(S70)、流速データ算定(S80)、河川表面の流速分布の出力(S90)を繰り返すのである。
図4の処理フローに示されるように、計測地点の決定(S10)、撮影アングルの決定(S20)、特徴点の測量(S30)、幾何変換パラメータの決定(S40)を行い、カメラ手段による撮影(S50)、撮影画像の幾何補正(S60)、PIV解析(S70)、流速データ算定(S80)、河川表面の流速分布の出力(S90)を繰り返すのである。
図5は、トータルステーションの機能ブロック図を示している。本発明の流速観測装置における角度計測手段と距離計測手段として、図5の機能ブロックに示されるような測距部と測角部とを一体としてなる測量用トータルステーションを用い、カメラ手段とトータルステーションとを一体化された構成とし、カメラ手段のレンズの投影中心座標からの角度として特徴点の水平角および鉛直角を計測するのである。
測角部53、測距部54、カメラ画像の各計測の入力情報は、光路変換部51の選択によって1個の望遠レンズ50により計測される。コントローラ56が測角部53,測距部54よびデジタルビデオカメラを制御して測角と測距および画像取得を行う。
操作・表示パネル59を用いて、コントローラ56の操作指示や取得されたデータ表示を行う。通信手段58は外部装置とネットワーク通信を行うものである。また、記録手段57は、測角と測距および画像をハードディスクなどのストレージに記録保存するものである。
測角部53、測距部54、カメラ画像の各計測の入力情報は、光路変換部51の選択によって1個の望遠レンズ50により計測される。コントローラ56が測角部53,測距部54よびデジタルビデオカメラを制御して測角と測距および画像取得を行う。
操作・表示パネル59を用いて、コントローラ56の操作指示や取得されたデータ表示を行う。通信手段58は外部装置とネットワーク通信を行うものである。また、記録手段57は、測角と測距および画像をハードディスクなどのストレージに記録保存するものである。
実施例1では、本発明の流速観測方法の有効性を確認する実験を示し、実際に計測が可能であること、従来法と比較して本発明の流速観測方法の精度を検証する。
本発明の流速観測方法は、屋外計測を対象としており、実河川の表面の如く、ある程度大きなスケールの空間を適当な精度で把握することを目的としている。本実験においては、大学敷地内の水質モニター小屋を計測対象とした。この小屋は矩形の建物であり、形状が明確であるため、誤差評価に適切であるからである。
本発明の流速観測方法は、屋外計測を対象としており、実河川の表面の如く、ある程度大きなスケールの空間を適当な精度で把握することを目的としている。本実験においては、大学敷地内の水質モニター小屋を計測対象とした。この小屋は矩形の建物であり、形状が明確であるため、誤差評価に適切であるからである。
画像の撮影に用いたカメラ手段は、ソニー製HDVカメラ(型番:HDR−FX1)を用いた。このカメラ手段では、1440×1080ピクセル(ピクセルのアスペクト比:1.333)の画像をインターレースで30FPSの動画として記録可能である。動画はMPEG2形式により圧縮されテープに記録される。
また、角度や距離の位置計測手段として、トータルステーション(Leica製,型番:TPS700)を用いた。このトータルステーションの測角の公称精度5秒(0.00027度),光波測距精度は5mm+3ppmである。
評価実験における7点の計測点(1〜7)を図6に示す。また、下記表1に計測点(1〜7)の水平角と鉛直角および斜距離の計測結果を示す。
また、角度や距離の位置計測手段として、トータルステーション(Leica製,型番:TPS700)を用いた。このトータルステーションの測角の公称精度5秒(0.00027度),光波測距精度は5mm+3ppmである。
評価実験における7点の計測点(1〜7)を図6に示す。また、下記表1に計測点(1〜7)の水平角と鉛直角および斜距離の計測結果を示す。
ここで、幾何変換の対象平面は建物屋上としている。建物屋上に相当する計測点1〜4の高さの平均は、−9.38(m)であり、幾何変換の対象平面は、z=−9.38(m)とした。図7は、幾何変換画像((a)幾何変換前、(b)幾何変換後)を示している。図7(a)にみられる建物屋上の歪みが、図7(b)において取り除かれており、建物の本来の平面形である矩形が再現されていることが確認できる。
次に、本発明の流速観測方法と従来の流速観測方法との計測精度を調査するために、幾何変換係数決定において利用する計測点の情報の組み合わせを下記表2に示す7通りに設定した。Case1〜3は本発明の流速観測方法に相当する情報であり、幾何変数パラメータの決定に利用する計測点の数を7点,4点,2点と変えている。Case4は与える情報が一番多く、今回の誤差検討ではこのCase4での結果からの偏差により誤差を評価することとする。Case5〜7は従来の流速観測方法に相当するもので、幾何変数パラメータの決定に利用する計測点の数を7点,4点,3点と変えている。最後にCase8は、従来の流速観測方法において計測点を2点とした場合に相当する。但し、2点の計測点でカメラ位置を逆算することは困難なことから、Case8ではカメラ位置を与えている。
さらに、計測誤差に対する感度を調べる目的で、トータルステーションによって計測された水平角および鉛直角に人為的に誤差を与えた組み合わせを3通り実施した。誤差は、計測点の7点の水平角および鉛直角のすべてに乱数で与えることとした。与えた誤差の標準偏差は0.1度である。
以上により、計測値をそのまま用いる1系列と、人為的な誤差を与える3系列があるため、総ケース数としては4系列×8ケース=32である。本実験の全32ケースのマトリックスを図8に示す。
以上により、計測値をそのまま用いる1系列と、人為的な誤差を与える3系列があるため、総ケース数としては4系列×8ケース=32である。本実験の全32ケースのマトリックスを図8に示す。
本実験の全32ケースの結果をグラフにまとめたものを図9に示す。図9のグラフの縦軸は誤差を示している。誤差の算出方法は、図7(a)中の点1,2,3,4の画像平面上の位置を、それぞれのケースで算出した幾何変換パラメータを用いて物理座標に変換したものを、計測値を用いたCase4で得られた結果と比較し、4点それぞれの位置のずれの最大距離を誤差としている。従って、ここでの誤差が計測値を用いたCase4を基準とした相対誤差であり、計測値自体の誤差やレンズ歪み等の影響は考慮されていない。図9のグラフで確認できるように、誤差を加えたケースでは計測値を用いるケースより誤差が大きくなっていることがわかる。また、用いる計測点の数が少ないCase3,7,8では誤差が大きくなっていることが確認できる。次に、ケースを系統的に整理して精度検討を行う。
本発明の流速観測方法の誤差評価とするために、計測値を用いたCase1〜3の結果をプロットしたグラフを図10に示す。図10のグラフから、計測値に大きな誤差が含まれない場合には、利用する計測点を7点から2点に減らしたとしても、幾何変換への影響はほとんどみられず、数cm程度の高い精度で幾何変換が可能であることが確認できる(計測点数が7点,4点,2点の場合の誤差は、それぞれ0.5cm,2cm,1cmである。)。
次に、図10のグラフに人為的な誤差を加えたケースを追加してプロットしたものを図11のグラフに示す。図11のグラフから、計測点の位置計測における誤差が大きい場合には、利用する計測点数が少なくなると幾何変換への誤差の影響が大きく、標準偏差0.1度程度の誤差が含まれる場合に2点の計測点を用いると最大50cm程度の位置のずれが見られることが確認できる。
次に、本発明の流速観測方法(Case1〜3)と従来の流速観測方法(Case5〜7)との比較のために、横軸を計測点数とし縦軸を誤差としたグラフに両方をプロットしたものを図12に示す。図12のグラフから、計測値に大きな誤差が含まれない場合、人為的な誤差を加えた場合の双方で、概ね、従来の流速観測方法に対して本発明の流速観測方法は同程度か誤差が小さくなっていることが確認できる。
本実験により、角度計測を利用する本発明の流速観測方法により幾何変換が可能であることが示されたことになる。幾何変換の精度については、本発明の流速観測方法を用いる場合、計測点の位置を高精度に観測可能な場合には、用いる計測点が2点であっても7点用いる場合と同程度の精度(数cm)で変換が可能であることが示された。
また、従来の流速観測方法と比較すると本発明の流速観測方法の精度は、同程度か本発明の流速観測方法の方がやや高精度という結果が示された。これは幾何変換パラメータの決定においては角度情報を用いる場合と、位置情報を用いる場合とで本質的な相違はなく、本発明の流速観測方法では角度情報に加えてカメラ位置も暗黙的に考慮されるため、従来の流速観測方法を利用する場合に比べて制約条件が増えるためと考えられる。
また、従来の流速観測方法と比較すると本発明の流速観測方法の精度は、同程度か本発明の流速観測方法の方がやや高精度という結果が示された。これは幾何変換パラメータの決定においては角度情報を用いる場合と、位置情報を用いる場合とで本質的な相違はなく、本発明の流速観測方法では角度情報に加えてカメラ位置も暗黙的に考慮されるため、従来の流速観測方法を利用する場合に比べて制約条件が増えるためと考えられる。
以上により、角度計測を利用する本発明の流速観測方法は、画像の幾何変換を行う際に有効な方法であるとともに、従来の流速観測方法に比べ精度は同程度かやや高精度であることが確認できた。本発明の流速観測方法では従来の流速観測方法で必要となる距離計測が不要となるため、この点で本発明の流速観測方法は測定作業における効率性および精度の両面で従来の流速観測方法に対する優位性を持っていることが理解されよう。
本発明の流速観測装置、流速観測方法および流速観測プログラムは、実河川の流速計測システムに利用することができる。
1 本発明の流速観測装置
2 河川
3 対岸の陸地
4 川岸の土手
5 河川表面の観測範囲
2 河川
3 対岸の陸地
4 川岸の土手
5 河川表面の観測範囲
Claims (11)
- 実河川の表面水流を撮影するカメラ手段と、前記カメラ手段のレンズの投影中心座標からの角度として、少なくとも2点の特徴点の水平角および鉛直角を計測する角度計測手段と、少なくとも水面高さから前記カメラ手段のレンズの投影中心までの高低差を計測する距離計測手段と、前記カメラ手段の撮影画像データと前記角度計測手段で得られた特徴点の水平角と鉛直角と前記距離計測手段で得られた水面高さからの高低差とを入力する解析手段とを備え、前記解析手段は、前記カメラ手段により撮影された実河川の表面水流の斜め画像を、前記角度計測手段および距離計測手段から決定される幾何変換パラメータを用いて垂直平面画像に幾何補正し、得られた垂直平面画像に対してPIV解析を実行し、実河川の水面各部の流速データを算定することを特徴とする流速観測装置。
- 実河川の表面水流を撮影するカメラ手段と、前記カメラ手段のレンズの投影中心座標からの角度として、少なくとも2点の特徴点の水平角および鉛直角を計測する角度計測手段と、河川の表面上の特徴点のうち少なくとも1点の距離を計測する距離計測手段と、前記カメラ手段の撮影画像データと前記角度計測手段で得られた特徴点の水平角と鉛直角と前記距離計測手段で得られたデータとを入力する解析手段とを備え、前記解析手段は、前記カメラ手段により撮影された実河川の表面水流の斜め画像を、前記角度計測手段および距離計測手段から決定される幾何変換パラメータを用いて垂直平面画像に幾何補正し、得られた垂直平面画像に対してPIV解析を実行し、実河川の水面各部の流速データを算定することを特徴とする流速観測装置。
- 前記角度計測手段と距離計測手段として、測距部と測角部を一体としてなる測量用トータルステーションを用い、前記カメラ手段と前記トータルステーションとを一体化して、前記カメラ手段のレンズの投影中心座標からの角度として特徴点の水平角および鉛直角を計測することを特徴とする請求項1又は2に記載の流速観測装置。
- 有線もしくは無線の通信ネットワークを介して外部装置と周期的に及び/又はイベント的にデータ通信を行う手段を更に備え、前記実河川の水面各部の流速データ及び/又は撮影画像データを周期的に及び/又はイベント的にネットワーク上に配信することを特徴とする請求項1乃至3のいずれか1項に記載の流速観測装置。
- 水平角と鉛直角の計測対象となる前記特徴点が撮影アングル中に存在しない場合において、河川の表面の見かけを部分的に変化させ得るレーザーポインタ手段を備えることを特徴とする請求項1乃至4のいずれか1項に記載の流速観測装置。
- 実河川の表面水流をカメラ手段によって撮影するステップと、
前記カメラ手段のレンズの投影中心座標からの角度として、少なくとも2点の特徴点の水平角および鉛直角を計測するステップと、
少なくとも水面高さから前記カメラ手段のレンズの投影中心までの高低差を計測するステップと、
前記カメラ手段の撮影画像データと特徴点の水平角と鉛直角と水面高さからの高低差とを解析手段に入力するステップと、
前記解析手段において、
前記カメラ手段により撮影された実河川の表面水流の斜め画像を、前記水平角と鉛直角および高低差から決定される幾何変換パラメータを用いて垂直平面画像に幾何補正するステップと、
得られた垂直平面画像に対してPIV解析を実行し、実河川の水面各部の流速データを算定するステップとを備えたことを特徴とする流速観測方法。 - 実河川の表面水流をカメラ手段によって撮影するステップと、
前記カメラ手段のレンズの投影中心座標からの角度として、少なくとも2点の特徴点の水平角および鉛直角を計測するステップと、
河川の表面上の特徴点のうち少なくとも1点の距離を計測するステップと、
前記カメラ手段の撮影画像データと特徴点の水平角と鉛直角と距離とを解析手段に入力するステップと、
前記解析手段において、
前記カメラ手段により撮影された実河川の表面水流の斜め画像を、前記水平角と鉛直角および距離から決定される幾何変換パラメータを用いて垂直平面画像に幾何補正するステップと、
得られた垂直平面画像に対してPIV解析を実行し、実河川の水面各部の流速データを算定するステップとを備えたことを特徴とする流速観測方法。 - 測距部と測角部とを一体としてなる測量用トータルステーションとカメラ手段とが一体化された装置を用いて、角度計測と高低差計測もしくは距離計測を行うことを特徴とする請求項6又は7に記載の流速観測方法。
- 実河川の表面水流を撮影するカメラ手段と、前記カメラ手段のレンズの投影中心座標からの角度として、少なくとも2点の特徴点の水平角および鉛直角を計測する角度計測手段と、少なくとも水面高さから前記カメラ手段のレンズの投影中心までの高低差もしくは前記特徴点との距離を計測する距離計測手段と、実河川の表面各部の流速分布を算定する解析コンピュータとを備えた流速観測システムにおいて、
前記解析コンピュータに、前記カメラ手段により撮影された実河川の表面水流の斜め画像を前記角度計測手段および距離計測手段から決定される幾何変換パラメータを用いて垂直平面画像に幾何補正させる手順と、得られた垂直平面画像に対してPIV解析を実行させ、実河川の水面各部の流速データを算定させる手順とを実行させるための流速観測プログラム。 - 有線もしくは無線の通信ネットワークを介して外部装置と周期的に及び/又はイベント的にデータ通信を行う手順を更に実行させることを特徴とする請求項9に記載の流速観測プログラム。
- 請求項9又は10に記載の流速観測プログラムを記録したコンピュータ読み取り可能な記録媒体。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2006238738A JP2008058264A (ja) | 2006-09-04 | 2006-09-04 | 実河川を対象とした流速観測装置、流速観測方法,および流速観測プログラム |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2006238738A JP2008058264A (ja) | 2006-09-04 | 2006-09-04 | 実河川を対象とした流速観測装置、流速観測方法,および流速観測プログラム |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2008058264A true JP2008058264A (ja) | 2008-03-13 |
Family
ID=39241154
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2006238738A Pending JP2008058264A (ja) | 2006-09-04 | 2006-09-04 | 実河川を対象とした流速観測装置、流速観測方法,および流速観測プログラム |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2008058264A (ja) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010123121A (ja) * | 2008-11-17 | 2010-06-03 | Honeywell Internatl Inc | シースルー・ディスプレイに現実世界の対象物の位置をマークする方法及び装置 |
CN102564508A (zh) * | 2011-12-14 | 2012-07-11 | 河海大学 | 基于视频图像的河流流量在线测验实现方法 |
JP2015004660A (ja) * | 2013-06-18 | 2015-01-08 | 財團法人國家實驗研究院National Applied Research Laboratories | モバイル式画像流速測定の方法及びその装置 |
CN107462397A (zh) * | 2017-08-14 | 2017-12-12 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种湖区超大范围表面流场测量方法 |
CN109085385A (zh) * | 2018-10-24 | 2018-12-25 | 内蒙古工业大学 | 一种基于piv靶盘标定多方向调节辅助盘及其测试方法 |
KR101996992B1 (ko) * | 2018-11-13 | 2019-07-08 | 주식회사 하이드로셈 | 옵티컬 플로우 영상 처리를 이용하는 하천 유속 측정 장치 및 방법 |
WO2020101104A1 (ko) * | 2018-11-13 | 2020-05-22 | 주식회사 하이드로셈 | Cctv 영상 기반의 실시간 자동 유량계측 시스템 및 방법 |
CN112067839A (zh) * | 2020-09-18 | 2020-12-11 | 重庆交通大学 | 无控制点条件下的无人机超宽河流表面流场快速测量方法 |
CN112362900A (zh) * | 2020-10-30 | 2021-02-12 | 天地伟业技术有限公司 | 一种无人值守的水位流速监测摄像机及监测方法 |
JP2021189029A (ja) * | 2020-05-29 | 2021-12-13 | 株式会社ハイドロ総合技術研究所 | 河川水流測定装置、該方法および該プログラムならびに記録媒体 |
CN113804916A (zh) * | 2021-09-17 | 2021-12-17 | 河海大学 | 一种基于最大流速先验信息的频域时空图像测速法 |
CN113819974A (zh) * | 2021-09-17 | 2021-12-21 | 河海大学 | 一种免水尺的河流水位视觉测量方法 |
CN116233370A (zh) * | 2023-04-27 | 2023-06-06 | 安徽哈斯特自动化科技有限公司 | 基于水质监测的智能化视频监控方法 |
-
2006
- 2006-09-04 JP JP2006238738A patent/JP2008058264A/ja active Pending
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010123121A (ja) * | 2008-11-17 | 2010-06-03 | Honeywell Internatl Inc | シースルー・ディスプレイに現実世界の対象物の位置をマークする方法及び装置 |
CN102564508A (zh) * | 2011-12-14 | 2012-07-11 | 河海大学 | 基于视频图像的河流流量在线测验实现方法 |
JP2015004660A (ja) * | 2013-06-18 | 2015-01-08 | 財團法人國家實驗研究院National Applied Research Laboratories | モバイル式画像流速測定の方法及びその装置 |
CN104280567A (zh) * | 2013-06-18 | 2015-01-14 | 林圣峰 | 行动式影像流速辨识的方法及其装置 |
CN107462397B (zh) * | 2017-08-14 | 2019-05-31 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种湖区超大范围表面流场测量方法 |
CN107462397A (zh) * | 2017-08-14 | 2017-12-12 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种湖区超大范围表面流场测量方法 |
CN109085385A (zh) * | 2018-10-24 | 2018-12-25 | 内蒙古工业大学 | 一种基于piv靶盘标定多方向调节辅助盘及其测试方法 |
CN109085385B (zh) * | 2018-10-24 | 2023-11-07 | 内蒙古工业大学 | 一种基于piv靶盘标定多方向调节辅助盘及其测试方法 |
US11353354B2 (en) | 2018-11-13 | 2022-06-07 | Hydrosem | System for automatically measuring discharge in real-time based on cctv image and method thereof |
KR101996992B1 (ko) * | 2018-11-13 | 2019-07-08 | 주식회사 하이드로셈 | 옵티컬 플로우 영상 처리를 이용하는 하천 유속 측정 장치 및 방법 |
WO2020101104A1 (ko) * | 2018-11-13 | 2020-05-22 | 주식회사 하이드로셈 | Cctv 영상 기반의 실시간 자동 유량계측 시스템 및 방법 |
WO2020101103A1 (ko) * | 2018-11-13 | 2020-05-22 | 주식회사 하이드로셈 | 옵티컬 플로우 영상 처리를 이용하는 하천 유속 측정 장치 및 방법 |
US11200684B2 (en) | 2018-11-13 | 2021-12-14 | Hydrosem | Apparatus and method for measuring flow velocity of stream using optical flow image processing |
JP2021189029A (ja) * | 2020-05-29 | 2021-12-13 | 株式会社ハイドロ総合技術研究所 | 河川水流測定装置、該方法および該プログラムならびに記録媒体 |
CN112067839A (zh) * | 2020-09-18 | 2020-12-11 | 重庆交通大学 | 无控制点条件下的无人机超宽河流表面流场快速测量方法 |
CN112067839B (zh) * | 2020-09-18 | 2022-04-12 | 重庆交通大学 | 无人机超宽河流表面流场快速测量方法 |
CN112362900A (zh) * | 2020-10-30 | 2021-02-12 | 天地伟业技术有限公司 | 一种无人值守的水位流速监测摄像机及监测方法 |
CN113819974A (zh) * | 2021-09-17 | 2021-12-21 | 河海大学 | 一种免水尺的河流水位视觉测量方法 |
CN113804916A (zh) * | 2021-09-17 | 2021-12-17 | 河海大学 | 一种基于最大流速先验信息的频域时空图像测速法 |
CN113819974B (zh) * | 2021-09-17 | 2022-07-12 | 河海大学 | 一种免水尺的河流水位视觉测量方法 |
CN116233370A (zh) * | 2023-04-27 | 2023-06-06 | 安徽哈斯特自动化科技有限公司 | 基于水质监测的智能化视频监控方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP2008058264A (ja) | 実河川を対象とした流速観測装置、流速観測方法,および流速観測プログラム | |
US10825198B2 (en) | 3 dimensional coordinates calculating apparatus, 3 dimensional coordinates calculating method, 3 dimensional distance measuring apparatus and 3 dimensional distance measuring method using images | |
US11477374B2 (en) | Three dimensional image capture system for imaging building facades using a digital camera, a near-infrared camera, and laser range finder | |
JP4607095B2 (ja) | 測量機器において画像処理を行う方法と装置 | |
US8699005B2 (en) | Indoor surveying apparatus | |
JP5150310B2 (ja) | 地理データ収集装置 | |
JP5027747B2 (ja) | 位置測定方法、位置測定装置、およびプログラム | |
JP6251142B2 (ja) | 測定対象物の非接触検知方法及びその装置 | |
GB2498177A (en) | Apparatus for determining a floor plan of a building | |
US20080267453A1 (en) | Method for estimating the pose of a ptz camera | |
JP2010014443A (ja) | 位置測定方法、位置測定装置、およびプログラム | |
JP5079547B2 (ja) | カメラキャリブレーション装置およびカメラキャリブレーション方法 | |
JP4052382B2 (ja) | 非接触画像計測装置 | |
JPWO2017130700A1 (ja) | 撮影支援装置及び撮影支援方法 | |
US20210264666A1 (en) | Method for obtaining photogrammetric data using a layered approach | |
JP4649192B2 (ja) | ステレオ画像作成方法及び3次元データ作成装置 | |
JP2013015429A (ja) | 移動基準点写真測量装置および方法 | |
JP6725675B2 (ja) | 自己位置推定装置、自己位置推定方法、プログラム、および画像処理装置 | |
JP4359084B2 (ja) | 測量システム | |
JPH11337322A (ja) | 2次元画像比較による外観の測定方法および装置 | |
JP4167509B2 (ja) | 測量システム | |
JP5409451B2 (ja) | 3次元変化検出装置 | |
JP4217083B2 (ja) | 測量システム | |
JP2009069135A (ja) | 撮像系校正装置と記録媒体 | |
Wackrow | Spatial measurement with consumer grade digital cameras |