JP2000205930A

JP2000205930A - 水位の周期的変動を計測する方法、および、同計測装置

Info

Publication number: JP2000205930A
Application number: JP11004537A
Authority: JP
Inventors: Yoichi Takagi; 陽市高木; Akio Tsujikawa; 秋雄辻川
Original assignee: Hitachi Ltd; Hitachi Process Computer Engineering Inc
Current assignee: Hitachi Ltd; Hitachi Information and Control Systems Inc
Priority date: 1999-01-11
Filing date: 1999-01-11
Publication date: 2000-07-28

Abstract

(57)【要約】【課題】河川４などの水面に立てられた量水板３（標
尺）をＴＶカメラ２で撮像して得られた画像情報信号を
画像処理装置で演算処理して、前記水面の水位を計測す
る技術を改良して、大型で高度なコンピュータを用いる
必要なく、水面の波浪（水位の周期的変動）の周期と振
幅とをリアルタイムで算出できるようにする。【解決手段】画像処理装置１の中に本発明に係る演算
回路および／または記憶回路を設けて、水位の周期的変
動を表す波動関数をフーリエ変換する。ただし、理論的
に完全なフーリエ解析するには高度な大型コンピュータ
と長時間の演算時間とを必要とするので、本発明におい
ては近似計算を行なうことによって、パソコン程度のコ
ンピュータによって迅速な算出を行なう。

Description

【発明の詳細な説明】

【０００１】

【発明の属する技術分野】本発明は、例えば河川，湖
沼，海洋などのように波浪を生じる水面の水位を計測す
る技術に関するものである。

【０００２】

【従来の技術】河川，湖沼，海洋などのように自然環境
の中に在る水面の水位を自動的に計測して監視し、また
は記録して研究資料とするための装置を大別すると、浮
子や圧力計などのように、水と直接的に接触するセンサ
を有する接触方式と、ＴＶカメラなどを用いて、水から
離れた地点から遠隔的に計測する非接触式とが有る。両
者にはそれぞれ長短が有るが、河川や湖沼は増水する場
合が有り、海洋には津波や高潮を発生する危険性が有る
ので、精密な機器を岸壁から離間させて設置し得る非接
触式が安全である。これらの水位を監視したり記録した
りすることは、平時におけるよりも、災害ないし準災害
時においてより重要であるから、増水や津波に関して安
全性が高いことの意義は大きい。

【０００３】図１２は、公知の撮像−画像処理方式の水
位計の構成と作用とを説明するために示した模式的な系
統図である。河川４の水面に、画像処理対象物である量
水板３が立てられており、該標尺の水面付近をＴＶカメ
ラ２が撮像して画像情報信号を出力する。

【０００４】上記画像情報信号は、画像処理装置１に入
力されて演算，解析され、水位信号を出力する。この水
位信号は、データ処理結果５としてプリントすることも
でき、ディスプレイ６で映像化することもできる。図示
を省略するが、上記水位信号を有線もしくは無線で遠隔
地へ送信することもできる。

【０００５】図１３は、前掲の図１２に示した画像処理
装置１の詳細な構成を描いた模式的な系統図である。主
メモリ１−２にはプログラムが格納されており、このプ
ログラムはＣＰＵ１−１や画像処理プロセッサ１−４で
実行される。画像処理プロセッサ１−４で算出された水
位計測値は、外部インターフェース１−６を介してプリ
ンタ２０に出力され、データ処理結果５がアウトプット
される。次に示す表１は、２秒間隔で３００秒間（５分
間）計測した水位計測データの１例である。

【０００６】

【表１】

【０００７】図１４は、時間を横軸にとって、表１に示
されている水位計測データをプロットしたグラフであ
る。

【０００８】

【発明が解決しようとする課題】河川，湖沼，海洋な
ど、自然環境の中に在る水面は一般に波浪を生じ、か
つ、該波浪は天候，地震などの影響を受けて著しく変動
する。波浪の有る水面の水位を計測したデータには、一
般に複数の周期と複数の振幅とを有する波が含まれてい
る。しかし、従来の技術によっては、波浪の周期や振幅
を情報として採取することが出来なかった。すなわち、
従来技術においては前掲の図１４に示したように周期的
に変動する水位をこの図のように表すか、または１分間
ないし１０分間毎の平均値（移動平均）を求めるという
程度のデータ加工しか行なわれていなかった。従来技術
における水位の平均化処理（移動平均）の１例を次の数
式１に示す。水位の平均値をＷａとすると、

【０００９】

【数１】

【００１０】ここに、Ｗ（ｔ）：水位の計測瞬時値ｔ：水位計測時刻ｔ−３０：平均処理開始時刻（直前３０秒前から平均処
理開始）。

【００１１】その上、例えばフロート式の水位計によれ
ば、浮子部材の慣性が作用するので波浪が有るときは２
〜７センチメートルの誤差を生じると言われている。ま
た、水圧計によって水位を測定しようとしても、水面に
生じた波の影響は、深い所に位置している水圧計に対し
ては平均化して伝えられるので、波浪の振幅を即時的に
計測することは原理的に無理である（非常に大きい誤差
の混入を避け難い）。本発明は上述の事情に鑑みて為さ
れたものであって、その目的とするところは、水面の波
の周期と振幅とを自動的に計測し得る方法と、上記の方
法を実施するに好適な装置とを提供するにある。

【００１２】

【課題を解決するための手段】上記の目的を達成するた
め、本発明においては水面位置（水位）を、一定の時間
間隔で繰り返して、画像処理によって計測する。上記一
定の時間間隔は、予想される波浪の周期よりも短いこと
が必要である。その理由は、波浪の周期と同程度もしく
はそれ以上の時間間隔で計測した結果を解析しても、波
浪の周期と正確に算出することは理論的に不可能だから
である。通常の計測条件においては２秒間隔程度で充分
である。現在の画像処理技術によれば０．２秒間隔の計
測も可能である。そして、本発明では、画像処理によっ
て得られた水位の時系列データを解析すること、例えば
フーリエ変換することにより波浪の周期や振幅を算出す
る。ただし、解析、例えばフーリエ変換を理論式どおり
に遂行するには高度のコンピュータと長時間とを必要と
するので、本発明は、パソコン程度の演算回路によって
リアルタイムに演算処理できるよう、演算処理量を節減
した。

【００１３】以下に、本発明の原理について述べる。水
位は、量水板（水面に立てた標尺）に接している水面の
瞬時値を計測する。時刻ｔ（秒）における水位計測結果
をｆ（ｔ）とする。このｆ（ｔ）は、例えば２秒毎にサ
ンプリングされた時系列的なデータである。上記の関数
ｆ（ｔ）をフーリエ級数によって近似する。ここで、ｆ
（ｘ）が、変数ｘに対して、２πを周期とする周期関数
であると仮定すると、次に揚げる数式２となる。

【００１４】

【数２】

【００１５】この関数式の右辺をフーリエ級数で近似す
るのであるが、この関数ｆ（ｘ）が、一定時間間隔毎に
採取されるデータであることに着目して、最小２乗法を
用いて次のように係数を定める。右辺のフーリエ級数の
式をＧ（ｘ）で表すと、

【００１６】

【数３】

【００１７】上掲の数式３を最小ならしめる係数を選択
すれば良い。

【００１８】係数ａ０を定めるには、次の数式４が成立
するようにａ０の値を求める。

【００１９】

【数４】

【００２０】以下に詳細計算する。

【００２１】

【数５】

【００２２】ここで、次の数式６に着目すると、

【００２３】

【数６】

【００２４】

【数７】

【００２５】上掲の数式７が得られる。

【００２６】前掲の数式２に示したｆ（ｘ）の関数式の
右辺における余弦の項は、

【００２７】

【数８】

【００２８】となるようにａｋを求めれば良い。以下、
詳細に計算する。

【００２９】

【数９】

【００３０】ここで、

【００３１】

【数１０】

【００３２】ただし、ｐ，ｑは任意の整数である。上掲
の数式１０に着目すると、ａｋは、

【００３３】

【数１１】

【００３４】上掲の数式１１によって得られる。

【００３５】次に、正弦の項は、

【００３６】

【数１２】

【００３７】上掲の数式１２が成立するようなｂｋの値
を求めることによって得られる。以下に、詳細な計算を
示す。

【００３８】

【数１３】

【００３９】ここで、次の数式１４に着目すると、

【００４０】

【数１４】

【００４１】ただし、ｐ，ｑは任意の整数である。

【００４２】

【数１５】

【００４３】上掲の数式１５が得られる。

【００４４】前掲の数式２，数式７，数式１１，および
数式１５により、ｆ（ｘ）の近似式を、次の数式１６の
如く立てる。

【００４５】

【数１６】

【００４６】その係数は、次式に掲げる数式１７によっ
て計算できることを確認した。

【００４７】

【数１７】

【００４８】時間Ｔを周期として水位が変動している場
合、上掲の近似式を水位に適用することにする。すなわ
ち、水位が周期Ｔ（秒）で変動するものとする。基準時
刻からの経過時間をｔ、水位をｗ（ｔ）で表すものとし
て、以下に、水位のフーリエ変換を求める。

【００４９】次掲の数式１８のような近似式を定義す
る。

【００５０】

【数１８】

【００５１】上掲の数式１８に、ｘ＝（２π／Ｔ）ｔ、
または、ｔ＝Ｔ＊ｘ／２πの関係を代入する。

【００５２】

【数１９】

【００５３】上掲の数式１９の係数は、

【００５４】

【数２０】

【００５５】この数式２０によって計算することができ
る。

【００５６】前掲の図１４と表１とに示した水位変動に
当てはめると、

【００５７】

【数２１】

【００５８】この場合、計測間隔時間Δｔ＝２秒であ
る。Ｔ＝３００（秒）を代入する。

【００５９】

【数２２】

【００６０】ここで、ｔ１：処理開始時刻ｔ２：処理終了時刻とする（Ｔ＝ｔ２−ｔ１）０〜３００秒の間を２秒毎に１５０等分してある。ｋは
１〜９９について計算し、高調波の有無を検知すること
にする。

【００６１】

【数２３】

【００６２】計算結果を次の表２の左半部に示す。

【００６３】

【表２】

【００６４】図１０は、横軸に時間をとって、表２の左
半部に結算結果として示された係数の値を表した棒グラ
フである。このグラフにおいて卓絶した大きい値を有す
る係数は、

【００６５】

【数２４】

【００６６】であり、他の係数は無視できる程度に小さ
い値である。

【００６７】

【数２５】

【００６８】計算結果は前掲の表２の右半部に示した。
図１１は、横軸に時間をとって、表２の右半部に結算結
果として示された係数の値を表した棒グラフである。こ
のグラフにおいて卓絶した大きい値を有する係数は、ｂ
１２＝２０．１６１（周期＝２５秒）のみであって、他
の係数は無視できる程度に小さい値である。

【００６９】以上を整理すると、水位は次の数式２６で
表される。

【００７０】

【数２６】

【００７１】以上に述べたようにして、画像による水位
計測に伴って波浪の周期や振幅を算出することができ
る。

【００７２】

【発明の実施の形態】本発明は、前掲の図１２および図
１３を参照して説明した、撮像→画像処理方式の水位計
測に併用して実施される。本発明が図１３（従来例）に
比して異なるところは、この図１３のような概要的な系
統図には表されない。すなわち、この図１３に示した主
メモリ１−２の中に、以下に述べるようなソフト処理機
能を与えれば、本発明方法を実施することのできる本発
明装置に改良される。図１は、本発明の１実施形態にお
けるソフト処理部の概要的な構成を描いたブロック図で
ある。本図に示したソフト処理部７はコンピュータプロ
グラムで実現される部分であり、前記の主メモリ１−２
に格納されていて、ＣＰＵ１−１（図１３参照）や画像
処理プロセッサ１−４で実行される。

【００７３】ソフト処理部７（図１参照）は、大きく分
けて、水位計測部７−１と、波浪計測部７−５、および
テーブル類８，９，１０から構成されている。水位計測
部７−１は、カメラ画像入力手段７−２，水面検知手段
７−３，および水位決定手段７−４から構成されてい
る。波浪計測部７−５は、データ範囲の更新手段７−
７，余弦０次常数項の更新手段７−７，要素データ計算
手段７−８，余弦０次以外の常数項および、正弦の常数
項更新手段７−９、並びに、データ変動周期の計算手段
７−１０から構成されている。図１に示した符号７−１
〜７−１０の手段は、それぞれ演算回路および／または
記憶回路によって構成されている。

【００７４】図２は、波浪計測部の主要な機能を示すブ
ロック図である。入力は計測データである。そして、こ
の波浪計測部７−５は、データ範囲の更新機能７−６，
０次の余弦常数項の更新機能７−７、要素データの計算
機能７−８，余弦０次以外の常数項および正弦常数項の
更新機能７−９，データ変動周期の計算機能７−１０に
分割されている。図３は、水位計測部の処理の１例を示
すフロー図である。カメラ画像を取り込み（ボックス
Ａ）、水面位置を検知し（ボックスＢ）、水位を決定す
る（ボックスＣ）。計測結果は計測データテーブル８に
記憶する。計測データ構成は、水位計測時刻８−１，水
位計測データ８−２，および、前回計測との計測時間間
隔８−３である。図３について以上に説明した計測デー
タテーブル８は、前掲の図１においても同一符号，名称
で表されている。

【００７５】図４は、水位計測処理の処理手順の１例を
示す系統図である。本図４に示す処理は水位計測毎に起
動され、処理途中のデータはテーブル類８，９，１０
に、一時的に記憶される。まず、最新の水位計測値を取
り込む（ボックスＤ）。データ範囲は、周期Ｔに相当す
る範囲である。この周期Ｔは、予想される波浪の周期に
比して充分に大きく取ることが望ましく、最少限、波浪
周期よりも大きく取る。この値は更新可能（ボックス
Ｅ）である。次に、余弦０次の常数項の更新を行なう
（ボックスＦ）。さらに、要素データの計算を行ない
（ボックスＧ）。要素データを使って余弦、正弦の常数
項の更新を行なう（ボックスＨ）。最後にデータ周期と
振幅とを算出する（ボックスＪ）。

【００７６】図５は、前掲の図４にステップＥとして示
したデータ範囲の更新処理の１例を表した模式図であ
る。この処理により、データ採取の間隔は、任意でも対
応することができる。また、例えば２秒毎に計測し、計
測範囲を固定することも可能であって、本実施例によれ
ば、任意範囲演算も固定範囲演算もできる。

【００７７】図６は、前掲の図４にステップＦとして示
したＯ次の余弦常数項を更新する補正処理の１例を示す
模式図である。この処理は、処理範囲外となる古いデー
タを減算するとともに、新しいデータを加算して行な
う。この処理によってデータが更新された１例を次の数
式２７に示す。

【００７８】

【数２７】

【００７９】上掲の数式２７で、変更前のデータ（第１
項）から古いデータ（第２項）を減算し、新しいデータ
（第３項）を加算することによって新しい係数が得られ
る。これにより、計算処理量が大幅に減少する。

【００８０】図７は、前掲の図４にステップＧとして示
した要素データの計算処理の１例を表した模式図であ
る。要素データは、余弦常数，正弦常数を更新するため
の要素データであって、最新の水位計測値に対する次の
数式２８である。

【００８１】

【数２８】

【００８２】この計測結果は、不要になるまで保存する
ものとする。図８は、前掲の図４においてステップＨと
して示したＯ次以外の余弦常数項と正弦常数項との更新
処理を表す模式図である。更新処理は、古いデータ部分
を減算し、新しいデータ部分を加算することによって行
なわれる（ボックスＨ１０〜Ｈ４０）。

【００８３】更新前の余弦の常数項は、

【００８４】

【数２９】

【００８５】更新後の余弦の常数項は、

【００８６】

【数３０】

【００８７】更新後の正弦の常数項は、

【００８８】

【数３１】

【００８９】更新後の正弦の常数項は、

【００９０】

【数３２】

【００９１】である。上掲の更新式に示したΔＡｋ，Δ
Ｂｋ、要素データとして、水位計測毎に補正されて一時
的に保存される。Ｔ時間後には不要となるので廃棄して
も良い。このようにして、要素データの加算，減算によ
る常数項の計算を行なうので、計算回数が大幅に減少す
る。その結果、高度のコンピュータを用いなくても、パ
ソコン程度のコンピュータによって、実用上リアルタイ
ムと見做し得るような迅速な演算が可能である。

【００９２】図９は、データの変動周期と振幅とを算出
するに際して、余弦常数項および正弦常数項のそれぞれ
について、比較的大きい値のみを取り出す処理を表した
模式図である。本例では、大きい順に５個を取り出すと
ともに、取り出された中で所定値よりも小さいものは捨
てる。これにより、演算の精度を余り下げずに演算処理
を迅速ならしめることができる。

【００９３】

【発明の効果】以上に本発明の実施形態を挙げて、その
構成と作用とを明らかならしめたように、本発明方法に
よれば、ＴＶカメラで撮像された画像情報信号を画像処
理して水位を計測する際、比較的簡単なコンピュータを
用いて、水面の周期的変動（波浪）の周期と振幅とを迅
速に算出することができる。また、本発明装置によれば
上記の発明方法を容易に実施して、その効果を充分に発
揮せしめることができる。

【図面の簡単な説明】

【図１】本発明の１実施形態におけるソフト処理部の概
要的な構成を描いたブロック図である。

【図２】波浪計測部の１実施形態における主要な機能を
示すブロック図である。

【図３】１実施形態における水位計測部の処理の１例を
示すフロー図である。

【図４】１実施形態ににおける水位計測処理の手順を示
す系統図である。

【図５】前掲の図４にステップＥとして示したデータ範
囲の更新処理の１例を表した模式図である。

【図６】前掲の図４にステップＦとして示した０次の余
弦常数項を更新する補正処理の１例を示す模式図であ
る。

【図７】前掲の図４にステップＧとして示した正弦常数
の計算処理の１例を表した模式図である。

【図８】前掲の図４にステップＨとして示した０次以外
の余弦常数項の正弦常数項との更新処理を表す模式図で
ある。

【図９】データの変動周期と振幅とを算出する演算にお
いて、余弦常数項および正弦常数項のそれぞれについ
て、比較的大きい値のみを取り出す処理を表した模式図
である。

【図１０】横軸に時間をとって、表２の左半部に計算結
果として示した係数の値を表した棒グラフである。

【図１１】横軸に時間をとって、表２の右半部に計算結
果として示した係数の値を表した棒グラフである。

【図１２】公知の撮像式水位計の構成と作用とを説明す
るために示した模式的な系統図である。

【図１３】前掲の図１２に示した画像処理装置の詳細な
構成を描いた模式的な系統図である。

【図１４】時間を横軸にとって、表１に示されている水
位計測データをプロットしたグラフである。

【符号の説明】

１…画像処理装置、１−１…ＣＰＵ、１−２…主メモ
リ、１−３…バス、１−４…画像処理プロセッサ、１−
５…画像メモリ、１−６…外部インターフェース、２…
ＴＶカメラ、３…画像処理対象物である量水板、４…河
川、５…データ処理結果、６…ディスプレイ、７…ソフ
ト処理部、７−１…水位計測部、７−２…カメラ画像入
力手段、７−３…水面検知手段、７−４…水位決定手
段、７−５…波浪計測部、７−６…データ範囲の更新手
段、７−７…余弦０次常数項の更新手段、７−８…要素
データ計算手段、７−９…０次以外の余弦常数項および
正弦常数項の更新手段、７−１０…大きい順に常数項を
選出する機能、８…計測データテーブル、９…時系列デ
ータテーブル、１０…時系列データ累積テーブル。

───────────────────────────────────────────────────── フロントページの続き (72)発明者辻川秋雄茨城県日立市大みか町五丁目２番１号株式会社日立製作所大みか工場内Ｆターム(参考） 2F014 FA04 GA01 2F065 AA22 AA24 BB01 BB05 BB27 CC00 DD06 FF04 FF26 FF31 FF61 JJ03 JJ19 QQ01 QQ16 QQ17 QQ23

Claims

【特許請求の範囲】

【請求項１】水面に立てられた量水板をＴＶカメラで
撮像して得られた画像情報を、画像処理プロセッサによ
り画像処理して前記水面の水位を計測する方法におい
て、所定の時間間隔で、繰り返して水位を計測して、水位を
時系列的に把握し、上記水位の時系列データを解析して、水面に発生してい
る波浪の周期および振幅を算出することを特徴とする、
水位の周期的変動を計測する方法。
【請求項２】前記の解析をフーリエ変換によって行な
い、かつ、フーリエ級数のうちで絶対の小さい項を削除
して近似値を求めることを特徴とする、請求項１に記載
した水位の周期的変動を計測する方法。
【請求項３】前記の水位を計測する時間間隔を、「発
生すると予想される波浪の周期」よりも短く設定し、複数の波浪が発生する可能性が有る場合は、これらの波
浪の中で最も短周期の予想波浪の周期よりも短く設定す
ることを特徴とする、請求項１または２に記載した水位
の周期的変動を計測する方法。
【請求項４】解析された近似フーリエ級数の各項を、
「波浪の周期よりも長い時間Ｔ」を周期として累積記憶
しておき、前記所定の時間間隔で行なわれた水位計測の計測値を得
る度に、最新の水位計測値を取り込み（ステップＤ）、前記近似フーリエ級数中の、余弦の０次常数項要素デー
タを、最新の水位計測値に基づいて補正して更新し（ス
テップＦ）、余弦の０次以外の常数項と、正弦の常数項との要素デー
タを補正して更新（ステップＨ）することにより、水面
水位の変動周期と振幅とを算出する（ステップＪ）こと
を特徴とする、請求項３に記載した水位の周期的変動を
計測する方法。
【請求項５】前記の余弦０次常数項要素データの更新
（ステップＦ）は、累積されている周期Ｔの期間内の余
弦０次常数項の中から、最も古いデータを減算し（ステ
ップＦ１０）、余弦０次常数項の累積値に対して最新の余弦０次常数項
データを加算（ステップＦ１０）して行なうことを特徴
とする、請求項４に記載した水位の周期的変動を計測す
る方法。
【請求項６】前記余弦の０次以外の常数項と正弦の常
数項との要素データの更新（ステップＨ）は、累積されている周期Ｔの期間内の余弦の０次以外の常数
項の中から最も古いデータを減算して、最新の余弦０次
以外の常数項を加算するとともに、累積されている周期Ｔの期間内の正弦の常数項の中から
最も古いデータを減算して、最新の正弦の常数項を加算
することによって行なう演算であることを特徴とする、
請求項４もしくは請求項５に記載した水位の周期的変動
を計測する方法。
【請求項７】前記の、水面水位の変動周期と振幅との
算出（ステップＪ）は、余弦の常数項の中から、大きい
順に所定個数の余弦常数項を取り出すとともに、取り出
した余弦常数項の中で所定値以下のものを削除し、か
つ、正弦の常数項の中から、大きい順に所定個数の正弦常数
項を取り出すとともに、取り出した正弦常数項の中で所
定値以下のものを削除する演算であることを特徴とす
る、請求項４ないし請求項６の何れかに記載した水位の
周期的変動を計測する方法。
【請求項８】水面に立てられた標尺を撮像して水面付
近の画像情報信号を出力するＴＶカメラと、上記画像情
報信号を演算解析して水位信号を出力する画像処理装置
とを具備している水位計測設備において、前記の画像処理装置は、所定の時間間隔で、水位の瞬時
値を繰り返し出力する機能を備えており、かつ、出力された多数の水位瞬時値の経時的変動をフー
リエ級数に解析する演算手段と、解析されたフーリエ級数の多数の項の中で主要な項を選
択して微小な項を削除する演算手段と、選択され、削除されなかった項に基づいて一定の時間間
隔Ｔで前記フーリエ級数の近似値を算出して、前記水面
に発生している波浪の周期や振幅を算出する演算手段
と、を具備していることを特徴とする、水位の周期的変
動を計測する装置。
【請求項９】前記の画像処理装置が水位の瞬時値を出
力する時間間隔は、発生を予想される波浪の周期よりも
短く設定されており、かつ、前記フーリエ級数の近似値を算出する時間間隔Ｔ
は、発生を予想される波浪の周期よりも長く設定されて
いることを特徴とする、請求項８に記載した水位の周期
的変動を計測する装置。
【請求項１０】前記の画像処理装置が出力した水位に
関するデータ、および／または上記水位に基づく演算途
中のデータを、前記の時間間隔Ｔの期間中、一時的に記
憶するデータテーブル、並びに、上記の時間間隔Ｔを周期として前記データテー
ブルに記憶されている数値を、新たに算出された水位デ
ータに基づいて補正，更新する演算回路を具備している
ことを特徴とする、請求項８もしくは請求項９に記載し
た水位の周期的変動を計測する装置。
【請求項１１】前記のデータテーブルに記憶されてい
る数値を補正，更新する演算回路は、フーリエ解析された余弦の常数項データおよび正弦の常
数項データを、新たな水位データが算出される度に、該新たな算出デー
タに基づいて更新する機能を有するものであることを特
徴とする、請求項１０に記載した水位の周期的変動を計
測する装置。
【請求項１２】前記のデータテーブルに記憶されてい
るデータに基づいて波浪の周期や振幅を算出する演算手
段が設けられていて、この演算手段は、フーリエ解析された余弦の常数項および正弦の常数項の
それぞれについて、大きい順から所定個数の常数項を取り出すとともに、取
り出された常数項の中で、所定の値よりも小さい常数項
を削除する機能を有するものであることを特徴とする、
請求項１０もしくは請求項１１に記載した水位の周期的
変動を計測する装置。