(19)【発行国】日本国特許庁(JP)
(12)【公報種別】公開特許公報(A)
(11)【公開番号】P2022114137
(43)【公開日】2022-08-05
(54)【発明の名称】波形情報算出装置、波形情報計算方法およびプログラム
(51)【国際特許分類】
G01S 17/95 20060101AFI20220729BHJP
G01C 13/00 20060101ALI20220729BHJP
G01D 3/02 20060101ALI20220729BHJP
【ＦＩ】
G01S17/95
G01C13/00 W
G01D3/02 N
【審査請求】未請求
【請求項の数】9
【出願形態】ＯＬ
(21)【出願番号】P 2021010302
(22)【出願日】2021-01-26
(71)【出願人】
【識別番号】599161890
【氏名又は名称】ＮＥＣネットワーク・センサ株式会社
(74)【代理人】
【識別番号】100106909
【弁理士】
【氏名又は名称】棚井 澄雄
(74)【代理人】
【識別番号】100134544
【弁理士】
【氏名又は名称】森 隆一郎
(74)【代理人】
【識別番号】100149548
【弁理士】
【氏名又は名称】松沼 泰史
(74)【代理人】
【識別番号】100162868
【弁理士】
【氏名又は名称】伊藤 英輔
(72)【発明者】
【氏名】小佐野 卓
【テーマコード（参考）】
2F075
5J084
【Ｆターム（参考）】
2F075AA05
2F075BB09
2F075EE15
2F075EE16
2F075EE18
2F075FF02
5J084AB13
5J084AC04
5J084CA22
5J084CA26
5J084CA31
5J084CA49
5J084EA11
(57)【要約】
【課題】波高または波長など波形に関する値を算出する際、波形における局所的な情報にとらわれて波の大きさを過少に評価する可能性を低減させられるようにする。
【解決手段】波形情報算出装置が、波形の測定データを取得するデータ取得手段と、前記測定データにおける極点を検出する極点検出手段と、複数の前記極点の位置関係に基づいて、所定の条件を満たす前記極点を除外する極点除外手段と、除外されなかった前記極点に基づいて、前記波形における波高または波長の少なくとも何れかを算出する波形情報算出手段と、を備える。
【選択図】図１
【解決手段】波形情報算出装置が、波形の測定データを取得するデータ取得手段と、前記測定データにおける極点を検出する極点検出手段と、複数の前記極点の位置関係に基づいて、所定の条件を満たす前記極点を除外する極点除外手段と、除外されなかった前記極点に基づいて、前記波形における波高または波長の少なくとも何れかを算出する波形情報算出手段と、を備える。
【選択図】図１
【特許請求の範囲】
【請求項1】
波形の測定データを取得するデータ取得手段と、
前記測定データにおける極点を検出する極点検出手段と、
複数の前記極点の位置関係に基づいて、所定の条件を満たす前記極点を除外する極点除外手段と、
除外されなかった前記極点に基づいて、前記波形における波高または波長の少なくとも何れかを算出する波形情報算出手段と、
を備える波形情報算出装置。
前記測定データにおける極点を検出する極点検出手段と、
複数の前記極点の位置関係に基づいて、所定の条件を満たす前記極点を除外する極点除外手段と、
除外されなかった前記極点に基づいて、前記波形における波高または波長の少なくとも何れかを算出する波形情報算出手段と、
を備える波形情報算出装置。
【請求項2】
前記極点除外手段は、連続する４つの前記極点のうち内側の２つの前記極点について、これら２つの極点における極値の差、または、少なくともこれら２つの極点に基づいて算出される波長の少なくとも何れかに基づいて、これら２つの極点を除外するか否かを決定する、
請求項１に記載の波形情報算出装置。
請求項１に記載の波形情報算出装置。
【請求項3】
前記極点除外手段は、前記４つの極点のうち内側の２つの前記極点について、これら２つの極点における極値の差の大きさが所定の閾値以下であり、かつ、これら２つの極点の何れも、前記４つの極点における最大値または最小値の何れにも該当しない場合、これら２つの極点を除外する、
請求項２に記載の波形情報算出装置。
請求項２に記載の波形情報算出装置。
【請求項4】
前記極点除外手段は、前記４つの極点のうち内側の２つの前記極点における極値の差の大きさの、前記４つの極点のうち外側の２つの前記極点における極値の差の大きさに対する割合が所定の閾値以下であり、かつ、前記内側の２つの極点の何れも、前記４つの極点における最大値または最小値の何れにも該当しない場合、前記内側の２つの極点を除外する、
請求項２または請求項３に記載の波形情報算出装置。
請求項２または請求項３に記載の波形情報算出装置。
【請求項5】
前記極点除外手段は、前記４つの極点のうち内側の２つの前記極点について、これら２つの極点と前記４つの極点のうち外側のいずれか１つの極点とに基づいて算出される波長が所定の閾値以下であり、かつ、前記内側の２つの極点の何れも、前記４つの極点における最大値または最小値の何れにも該当しない場合、前記内側の２つの極点を除外する、
請求項２から４の何れか一項に記載の波形情報算出装置。
請求項２から４の何れか一項に記載の波形情報算出装置。
【請求項6】
前記極点除外手段は、前記測定データにおける極点のうち外側の前記極点について、前記外側の極点と、その極点に隣接する極点との極値の差の大きさが所定の閾値以下であり、かつ、前記外側の極点が、その極点を含めて外側から３つの極点における最大値または最小値の何れにも該当しない場合、前記外側の極点を除外する、
請求項１から５の何れか一項に記載の波形情報算出装置。
請求項１から５の何れか一項に記載の波形情報算出装置。
【請求項7】
前記極点除外手段は、極大点を中心に連続する５つの極点について、所定の波高算出方法で波高を算出できないか、または、所定の波長算出方法で波長を算出できない場合、その５つの極点の組み合わせを波形情報算出対象から除外し、
前記波形情報算出手段は、除外されなかった前記連続する５つの極点に基づいて、中心の前記極大点における波高および波形を算出する、
請求項１から６の何れか一項に記載の波形情報算出装置。
前記波形情報算出手段は、除外されなかった前記連続する５つの極点に基づいて、中心の前記極大点における波高および波形を算出する、
請求項１から６の何れか一項に記載の波形情報算出装置。
【請求項8】
波形の測定データを取得することと、
前記測定データにおける極点を検出することと、
複数の前記極点の位置関係に基づいて、所定の条件を満たす前記極点を除外することと、
除外されなかった前記極点に基づいて、前記波形における波高または波長の少なくとも何れかを算出することと、
を含む波形情報算出方法。
前記測定データにおける極点を検出することと、
複数の前記極点の位置関係に基づいて、所定の条件を満たす前記極点を除外することと、
除外されなかった前記極点に基づいて、前記波形における波高または波長の少なくとも何れかを算出することと、
を含む波形情報算出方法。
【請求項9】
コンピュータに、
波形の測定データを取得することと、
前記測定データにおける極点を検出することと、
複数の前記極点の位置関係に基づいて、所定の条件を満たす前記極点を除外することと、
除外されなかった前記極点に基づいて、前記波形における波高または波長の少なくとも何れかを算出することと、
を実行させるためのプログラム。
波形の測定データを取得することと、
前記測定データにおける極点を検出することと、
複数の前記極点の位置関係に基づいて、所定の条件を満たす前記極点を除外することと、
除外されなかった前記極点に基づいて、前記波形における波高または波長の少なくとも何れかを算出することと、
を実行させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、波形情報算出装置、波形情報計算方法およびプログラムに関する。
【背景技術】
【0002】
波形に関する処理について、幾つかの技術が提案されている。
例えば、特許文献１には、水面に対して斜めにレーザビームを照射した場合に、レーザビームが波浪の谷部に届いているか否かを判定し、届いていないと判定した場合は、所定の式を用いて波高を計算することが記載されている。
また、特許文献２には、センサからサンプリングされるデータの傾きに基づいて極値のデータを検出し、検出した極値のデータのみを送信することで、通信量を削減することが記載されている。
例えば、特許文献１には、水面に対して斜めにレーザビームを照射した場合に、レーザビームが波浪の谷部に届いているか否かを判定し、届いていないと判定した場合は、所定の式を用いて波高を計算することが記載されている。
また、特許文献２には、センサからサンプリングされるデータの傾きに基づいて極値のデータを検出し、検出した極値のデータのみを送信することで、通信量を削減することが記載されている。
【先行技術文献】
【特許文献】
【0003】
【特許文献1】特開２０１７－１３３９０２号公報
【特許文献2】特開２０１７－０１６１８７号公報
【発明の概要】
【発明が解決しようとする課題】
【0004】
波高または波長など波形に関する値を算出する際、波形における局所的な情報にとらわれて波の大きさを過少に評価することを避けたい場合が考えられる。
【0005】
本発明の目的の一例は、上述の課題を解決することのできる波形情報算出装置、波形情報計算方法およびプログラムを提供することである。
【課題を解決するための手段】
【0006】
本発明の第１の態様によれば、波形情報算出装置は、波形の測定データを取得するデータ取得手段と、前記測定データにおける極点を検出する極点検出手段と、複数の前記極点の位置関係に基づいて、所定の条件を満たす前記極点を除外する極点除外手段と、除外されなかった前記極点に基づいて、前記波形における波高または波長の少なくとも何れかを算出する波形情報算出手段と、を備える。
【0007】
本発明の第２の態様によれば、波形情報算出方法は、波形の測定データを取得することと、前記測定データにおける極点を検出することと、複数の前記極点の位置関係に基づいて、所定の条件を満たす前記極点を除外することと、除外されなかった前記極点に基づいて、前記波形における波高または波長の少なくとも何れかを算出することと、を含む。
【0008】
本発明の第３の態様によれば、プログラムは、コンピュータに、波形の測定データを取得することと、前記測定データにおける極点を検出することと、複数の前記極点の位置関係に基づいて、所定の条件を満たす前記極点を除外することと、除外されなかった前記極点に基づいて、前記波形における波高または波長の少なくとも何れかを算出することと、を実行させるためのプログラムである。
【発明の効果】
【0009】
この発明によれば、波高または波長など波形に関する値を算出する際、波形における局所的な情報にとらわれて波の大きさを過少に評価する可能性を低減させることができる。
【図面の簡単な説明】
【0010】
【図1】実施形態に係る波形情報算出装置の機能構成の例を示す概略ブロック図である。
【図2】波形の例を示す図である。
【図3】波形にローパスフィルタを適用した場合の例を示す図である。
【図4】測定データにおける、連続する４つの極点の第１例を示す図である。
【図5】測定データにおける、連続する４つの極点の第２例を示す図である。
【図6】測定データにおける、連続する４つの極点の第３例を示す図である。
【図7】測定データにおける、連続する４つの極点の第４例を示す図である。
【図8】測定データにおける、連続する４つの極点の第５例を示す図である。
【図9】測定データにおける、連続する４つの極点の第６例を示す図である。
【図10】測定データにおける、連続する４つの極点の第７例を示す図である。
【図11】測定データにおける、連続する４つの極点の第８例を示す図である。
【図12】測定データにおける、連続する４つの極点の第９例を示す図である。
【図13】測定データにおける、連続する４つの極点の第１０例を示す図である。
【図14】測定データにおける、連続する５つの極点の例を示す図である。
【図15】実施形態に係る波形情報算出装置が波形情報を算出する処理の手順の第１例を示すフローチャートである。
【図16】実施形態に係る波形情報算出装置が波形情報を算出する処理の手順の第２例を示すフローチャートである。
【図17】実施形態に係る制御部の構成例を示す概略ブロック図である。
【図18】実施形態に係る極点除外部が除外する極値の例を示す図である。
【図19】実施形態に係る波形情報算出装置の構成例を示す図である。
【図20】実施形態に係る波形情報算出方法における処理の手順の例を示すフローチャートである。
【図21】少なくとも１つの実施形態に係るコンピュータの構成を示す概略ブロック図である。
【発明を実施するための形態】
【0011】
以下、本発明の実施形態を説明するが、以下の実施形態は請求の範囲にかかる発明を限定するものではない。また、実施形態の中で説明されている特徴の組み合わせの全てが発明の解決手段に必須であるとは限らない。
図１は、実施形態に係る波形情報算出装置の機能構成の例を示す概略ブロック図である。図１に示す構成で、波形情報算出装置１００は、通信部１１０と、表示部１２０と、操作入力部１３０と、記憶部１８０と、制御部１９０とを備える。制御部１９０は、データ取得部１９１と、極点検出部１９２と、極点除外部１９３と、波形情報算出部１９４とを備える。
図１は、実施形態に係る波形情報算出装置の機能構成の例を示す概略ブロック図である。図１に示す構成で、波形情報算出装置１００は、通信部１１０と、表示部１２０と、操作入力部１３０と、記憶部１８０と、制御部１９０とを備える。制御部１９０は、データ取得部１９１と、極点検出部１９２と、極点除外部１９３と、波形情報算出部１９４とを備える。
【0012】
波形情報算出装置１００は、波形の測定データを取得して、波高または波長の少なくとも何れかを算出する。波形情報算出装置１００は、例えばパソコン（Personal Computer；ＰＣ）などのコンピュータを用いて構成されていてもよい。
以下では、波形情報算出装置１００が、海面の波を対象とする場合を例に説明する。ただし、波形情報算出装置１００が対象とする波は、特定のものに限定されない。
波高または波長など波形に関する値を示す情報を、波形情報と称する。波形に関する値を算出することを、波形情報を算出すると称する。
以下では、波形情報算出装置１００が、海面の波を対象とする場合を例に説明する。ただし、波形情報算出装置１００が対象とする波は、特定のものに限定されない。
波高または波長など波形に関する値を示す情報を、波形情報と称する。波形に関する値を算出することを、波形情報を算出すると称する。
【0013】
波形情報を算出する際の問題点として、波形における局所的な情報にとらわれて波の大きさを過少に評価することが考えられる。
図２は、波形の例を示す図である。図２のグラフの横軸（ｘ軸）は、水平方向（波長の方向）の位置を示す座標軸である。あるいは、同一地点における波の変化を観測する場合、図２のグラフの横軸は時間軸であってもよい。図２のグラフの縦軸（ｙ軸）は、垂直方向（波高の方向）の位置を示す座標軸である。
図２は、波形の例を示す図である。図２のグラフの横軸（ｘ軸）は、水平方向（波長の方向）の位置を示す座標軸である。あるいは、同一地点における波の変化を観測する場合、図２のグラフの横軸は時間軸であってもよい。図２のグラフの縦軸（ｙ軸）は、垂直方向（波高の方向）の位置を示す座標軸である。
【0014】
図２の波形について、点Ｐ１１を波の谷底とし、点Ｐ１４を波の山頂として波高ｈ１３を算出したい場合について考える。
これに対して、点Ｐ１２を波の山頂として扱うと、波高ｈ１１を算出することが考えられる。これは、波の局所的な極大点Ｐ１２を波の山頂として扱ったことで、波高をｈ１３より小さいｈ１１と算出したものと言える。
また、点Ｐ１３を波の谷底として扱うと、波高ｈ１２を算出することが考えられる。これは、波の局所的な極小点Ｐ１３を波の谷底として扱ったことで、波高をｈ１３より小さいｈ１２と算出したものと言える。
これに対して、点Ｐ１２を波の山頂として扱うと、波高ｈ１１を算出することが考えられる。これは、波の局所的な極大点Ｐ１２を波の山頂として扱ったことで、波高をｈ１３より小さいｈ１１と算出したものと言える。
また、点Ｐ１３を波の谷底として扱うと、波高ｈ１２を算出することが考えられる。これは、波の局所的な極小点Ｐ１３を波の谷底として扱ったことで、波高をｈ１３より小さいｈ１２と算出したものと言える。
【0015】
点Ｐ１２およびＰ１３のような局所的な極点は、測定対象の波に実際に生じる場合だけでなく、波形の測定状況に応じて生じる場合が考えられる。例えば、波形測定のための装置が飛行艇などの移動体に搭載される場合、移動体の動揺によって、測定結果の波形に局所的な極点が生じる可能性がある。また、海面の反射点の変動によって、測定結果の波形に局所的な極点が生じる可能性がある。
【0016】
波形の局所的な変化を無視して大局的な値を算出するための方法の１つとして、ローパスフィルタ（Low Pass Filter）を用いて、波形の局所的な変化の情報を除去する（波形をなまらせる）ことが考えられる。しかしながら、ローパスフィルタの適用によって波形が変化することで、波高などの波形情報に誤差が生じることが考えられる。
【0017】
図３は、波形にローパスフィルタを適用した場合の例を示す図である。図３のグラフの横軸（ｘ軸）および縦軸（ｙ軸）は、図２の場合と同様である。
図３は、図２の波形にローパスフィルタを適用した場合の例を示している。線Ｌ１１は、図２の波形を示しており、ローパスフィルタ適用前の波形に該当する。線Ｌ１２は、ローパスフィルタ適用後の波形を示す。
図３は、図２の波形にローパスフィルタを適用した場合の例を示している。線Ｌ１１は、図２の波形を示しており、ローパスフィルタ適用前の波形に該当する。線Ｌ１２は、ローパスフィルタ適用後の波形を示す。
【0018】
線Ｌ１１の波形が、点Ｐ１１からＰ１４までの４つの極点を有するのに対し、線Ｌ１２の波形は、点Ｐ２１およびＰ２２の２つの極点を有する。このように、線Ｌ１２の波形では、点Ｐ１２およびＰ１３に相当する極点が無くなっている。これにより、図２の波高ｈ１１およびｈ１２に相当する、いわば局所的な波高を算出することを回避できる。
【0019】
一方、ローパスフィルタ適用後の波形（線Ｌ１２）は、ローパスフィルタ適用前の波形（線Ｌ１１）から変化している。これにより、線Ｌ１２の波高ｈ２３は、線Ｌ１３の波高ｈ１３よりも小さくなっている。
ローパスフィルタ適用後のデータに基づいて波高を算出する場合、波高ｈ１３と波高ｈ２３との差のように、ローパスフィルタの適用によって波形が変化することによる誤差が生じる。
ローパスフィルタ適用後のデータに基づいて波高を算出する場合、波高ｈ１３と波高ｈ２３との差のように、ローパスフィルタの適用によって波形が変化することによる誤差が生じる。
【0020】
そこで、波形情報算出装置１００は、波形における極点のうち所定の条件を満たす極点を、波形情報の算出に用いる極点から除外する。これにより、波形情報算出装置１００では、波形における局所的な情報にとらわれて波の大きさを過少に評価する可能性を低減させ、かつ、波形の変化（歪み）による誤差を回避するか、あるいは低減させて、波形情報を算出できると期待される。波形情報算出装置１００によれば、この点で、波形情報を高精度に算出することができる。
図２の点Ｐ１２を波の山頂とする波のように、除外される極点に基づく波を不要波とも称する。
図２の点Ｐ１２を波の山頂とする波のように、除外される極点に基づく波を不要波とも称する。
【0021】
通信部１１０は、他の装置と通信を行う。例えば、通信部１１０は、データ取得部１９１の制御に従って波形測定用のセンサと通信を行い、波形の測定値を示すセンサ信号を受信する。
表示部１２０は、例えば液晶パネルまたはＬＥＤ（Light Emitting Diode、発光ダイオード）パネルなどの表示画面を備え、各種画像を表示する。表示部１２０が、波形の測定データに基づく波形と、波形情報算出装置１００が算出する波形情報とを表示するようにしてもよい。
操作入力部１３０は、例えば表示部１２０の表示画面に設けられてタッチパネルを構成するタッチセンサなどの入力デバイスを備え、ユーザ操作を受け付ける。例えば操作入力部１３０が、波形情報の算出を指示するユーザ操作を受け付けるようにしてもよい。
表示部１２０は、例えば液晶パネルまたはＬＥＤ（Light Emitting Diode、発光ダイオード）パネルなどの表示画面を備え、各種画像を表示する。表示部１２０が、波形の測定データに基づく波形と、波形情報算出装置１００が算出する波形情報とを表示するようにしてもよい。
操作入力部１３０は、例えば表示部１２０の表示画面に設けられてタッチパネルを構成するタッチセンサなどの入力デバイスを備え、ユーザ操作を受け付ける。例えば操作入力部１３０が、波形情報の算出を指示するユーザ操作を受け付けるようにしてもよい。
【0022】
記憶部１８０は、各種情報を記憶する。例えば、記憶部１８０は、波形の測定データを記憶する。記憶部１８０は、波形情報算出装置１００が備える記憶デバイスを用いて構成される。
制御部１９０は、波形情報算出装置１００の各部を制御して各種処理を行う。制御部１９０の機能は、例えば、波形情報算出装置１００が備えるＣＰＵ（Central Processing Unit、中央処理装置）が、記憶部１８０からプログラムを読み出して実行することで実行される。
制御部１９０は、波形情報算出装置１００の各部を制御して各種処理を行う。制御部１９０の機能は、例えば、波形情報算出装置１００が備えるＣＰＵ（Central Processing Unit、中央処理装置）が、記憶部１８０からプログラムを読み出して実行することで実行される。
【0023】
データ取得部１９１は、波形の測定データを取得する。例えばデータ取得部１９１は、通信部１１０が受信するセンサ信号から波形の測定データを読み出す。
データ取得部１９１は、データ取得手段の例に該当する。
データ取得部１９１が制御部１９０の外部に設けられていてもよい。例えば、波形情報算出装置１００が、波形測定用のセンサを備え、そのセンサがデータ取得部１９１を構成していてもよい。
データ取得部１９１は、データ取得手段の例に該当する。
データ取得部１９１が制御部１９０の外部に設けられていてもよい。例えば、波形情報算出装置１００が、波形測定用のセンサを備え、そのセンサがデータ取得部１９１を構成していてもよい。
【0024】
データ取得部１９１が取得する波形の測定データは、波形のサンプリング点ごとに、水平方向（波長の方向）の座標値と、垂直方向（波高の方向）の座標値とを示すデータであってもよい。あるいは、波形測定用のセンサが、同一地点における波の変化を測定する場合、データ取得部１９１が取得する波形の測定データは、波形のサンプリング点ごとに、測定時刻と、垂直方向（波高の方向）の座標値とを示すデータであってもよい。
以下では、水平方向の座標または時刻座標をｘ座標と称し、垂直方向の座標をｙ座標と称する。波形の測定データを、単に測定データとも称する。波形のサンプリング点を単にサンプリング点とも称する。
以下では、水平方向の座標または時刻座標をｘ座標と称し、垂直方向の座標をｙ座標と称する。波形の測定データを、単に測定データとも称する。波形のサンプリング点を単にサンプリング点とも称する。
【0025】
極点検出部１９２は、データ取得部１９１が取得する測定データにおける極点を検出する。
極点検出部１９２は、極点検出手段の例に該当する。
極点検出部１９２が、ｘ座標値が小さい順にサンプリング点のｙ座標値を参照し、ｙ座標値の変曲点を極点として検出するようにしてもよい。ここでいう変曲点は、ｙ座標値の増減が切り替わる点である。極点検出部１９２が、ｙ座標値が増加から減少に転じる点を極大点として検出し、ｙ座標値が減少から増加に転じる点が極小点として検出するようにしてもよい。
極点検出部１９２は、極点検出手段の例に該当する。
極点検出部１９２が、ｘ座標値が小さい順にサンプリング点のｙ座標値を参照し、ｙ座標値の変曲点を極点として検出するようにしてもよい。ここでいう変曲点は、ｙ座標値の増減が切り替わる点である。極点検出部１９２が、ｙ座標値が増加から減少に転じる点を極大点として検出し、ｙ座標値が減少から増加に転じる点が極小点として検出するようにしてもよい。
【0026】
極点除外部１９３は、極点検出部１９２が検出する複数の極点の位置関係に基づいて、所定の条件を満たす極点を除外する。
極点除外部１９３は、極点除外手段の例に該当する。
極点除外部１９３は、極点除外手段の例に該当する。
【0027】
波形情報算出部１９４は、極点除外部１９３によって除外されなかった極点に基づいて、波形における波高または波長の少なくとも何れかを算出する。
波形情報算出部１９４は、波形情報算出手段の例に該当する。
波形情報算出部１９４が、隣り合う極大点と極小点とのｙ座標値の差を波高として算出するようにしてもよい。ここでいう、極大点と極小点が隣り合うことは、ｘ座標方向（ｘ軸方向）において、その極大点と極小点との間に他の極点が存在し無いことである。
波形情報算出部１９４は、波形情報算出手段の例に該当する。
波形情報算出部１９４が、隣り合う極大点と極小点とのｙ座標値の差を波高として算出するようにしてもよい。ここでいう、極大点と極小点が隣り合うことは、ｘ座標方向（ｘ軸方向）において、その極大点と極小点との間に他の極点が存在し無いことである。
【0028】
波形情報算出部１９４が、極大点とｘ座標方向に次の極大点とのｘ座標値の差を波長として算出するようにしてもよい。ここでいう、第１の極大点からｘ座標方向に次の極大点が第２の極大点であることは、ｘ座標方向において、それら２つの極大点の間に１つの極小点のみが存在し、他の極点が存在しないことである。
波形情報算出部１９４が、極小点とｘ座標方向に次の極小点とのｘ座標値の差を波長として算出するようにしてもよい。ここでいう、第１の極小点からｘ座標方向に次の極小点が第２の極小点であることは、ｘ座標方向において、それら２つの極小点の間に１つの極大点のみが存在し、他の極点が存在しないことである。
波形情報算出部１９４が、極小点とｘ座標方向に次の極小点とのｘ座標値の差を波長として算出するようにしてもよい。ここでいう、第１の極小点からｘ座標方向に次の極小点が第２の極小点であることは、ｘ座標方向において、それら２つの極小点の間に１つの極大点のみが存在し、他の極点が存在しないことである。
【0029】
極点除外部１９３が、連続する４つの極点のうち内側の２つの極点について、これら２つの極点における極値の差、または、少なくともこれら２つの極点に基づいて算出される波長の少なくとも何れかに基づいて、これら２つの極点を除外するか否かを決定するようにしてもよい。
ここでいう、複数の極点が連続することは、ｘ座標方向において、これら複数の極点の間に他の極点が存在しないことである。ここでいう、極点を除外することは、その極点を、波形情報算出部１９４が波形情報を計算する際に用いる極点から除外することである。
ここでいう、複数の極点が連続することは、ｘ座標方向において、これら複数の極点の間に他の極点が存在しないことである。ここでいう、極点を除外することは、その極点を、波形情報算出部１９４が波形情報を計算する際に用いる極点から除外することである。
【0030】
例えば、極点除外部１９３が、連続する４つの極点のうち内側の２つの極点について、これら２つの極点における極値の差の大きさが所定の閾値以下であり、かつ、これら２つの極点の何れも、４つの極点における最大値または最小値の何れにも該当しない場合、これら２つの極点を除外するようにしてもよい。
【0031】
また、極点除外部１９３が、連続する４つの極点のうち内側の２つの極点における極値の差の大きさの、４つの極点のうち外側の２つの極点における極値の差の大きさに対する割合が所定の閾値以下であり、かつ、内側の２つの極点の何れも、４つの極点における最大値または最小値の何れにも該当しない場合、内側の２つの極点を除外するようにしてもよい。
【0032】
また、極点除外部１９３が、連続する４つの極点のうち内側の２つの極点について、これら２つの極点と４つの極点のうち外側のいずれか１つの極点とに基づいて算出される波長が所定の閾値以下であり、かつ、内側の２つの極点の何れも、４つの極点における最大値または最小値の何れにも該当しない場合、内側の２つの極点を除外するようにしてもよい。
【0033】
図４は、測定データにおける、連続する４つの極点の第１例を示す図である。図４のグラフの横軸は上述したｘ軸を示し、縦軸はｙ軸を示す。
以下では、連続する４つの極点を、ｘ座標値の小さい順に点ａ、ｂ、ｃ、ｄと表記する。図４の例では、点ａおよびｃが極大点に該当し、点ｂおよびｄが極小点に該当する。
以下では、連続する４つの極点を、ｘ座標値の小さい順に点ａ、ｂ、ｃ、ｄと表記する。図４の例では、点ａおよびｃが極大点に該当し、点ｂおよびｄが極小点に該当する。
【0034】
点ａ、ｂ、ｃ、ｄのｘ座標値を、それぞれｘａ、ｘｂ、ｘｃ、ｘｄと表記する。点ａ、ｂ、ｃ、ｄのｙ座標値を、それぞれｙａ、ｙｂ、ｙｃ、ｙｄと表記する。
点ａおよび点ｄに基づく波高を波高ｈａｄと表記する。点ａが極大点の場合、ｈａｄ＝ｙａ－ｙｄと表される。点ｂおよび点ｃに基づく波高を波高ｈｂｃと表記する。点ａが極大点の場合、点ｂが極小点、点ｃが極大点であり、ｈｂｃ＝ｙｃ－ｙｂと表される。
点ａおよび点ｄに基づく波高を波高ｈａｄと表記する。点ａが極大点の場合、ｈａｄ＝ｙａ－ｙｄと表される。点ｂおよび点ｃに基づく波高を波高ｈｂｃと表記する。点ａが極大点の場合、点ｂが極小点、点ｃが極大点であり、ｈｂｃ＝ｙｃ－ｙｂと表される。
【0035】
点ａから点ｃまでは、波の１周期に該当し、ｘｃ－ｘａは、波長に該当する。この波長を、λａｃと表記する。点ｂから点ｄまでは、波の１周期に該当し、ｘｄ－ｘｂは、波長に該当する。この波長を、λｂｄと表記する。
波形情報算出部１９４が、これらの波高および波長を算出するようにしてもよい。あるいは、極点除外部１９３が、これらの波高および波長を算出するようにしてもよい。
極点除外部１９３が、式（１）が成立する場合に、点ｂおよび点ｃを除外するようにしてもよい。
波形情報算出部１９４が、これらの波高および波長を算出するようにしてもよい。あるいは、極点除外部１９３が、これらの波高および波長を算出するようにしてもよい。
極点除外部１９３が、式（１）が成立する場合に、点ｂおよび点ｃを除外するようにしてもよい。
【0036】
【数1】
【0037】
「∧」は、「かつ」（論理積）を表す。「∨」は、「または」（論理和）を表す。
「ｈｂｃ≦０．３」の「０．３」は、０．３メートル（ｍ）を示す。この「０．３」は、波高の閾値の例に該当する定数である。ただし、波高の閾値は、特定の値に限定されない。
「λａｃ<５」の「５」、および、「λｂｄ<５」の「５」は、何れも５メートルを示す。これらの「５」は、波長の閾値の例に該当する定数である。ただし、波長の閾値は、特定の値に限定されない。
「ｈｂｃ≦０．３」の「０．３」は、０．３メートル（ｍ）を示す。この「０．３」は、波高の閾値の例に該当する定数である。ただし、波高の閾値は、特定の値に限定されない。
「λａｃ<５」の「５」、および、「λｂｄ<５」の「５」は、何れも５メートルを示す。これらの「５」は、波長の閾値の例に該当する定数である。ただし、波長の閾値は、特定の値に限定されない。
【0038】
「ｈｂｃ≦０．３」は、２つの極点における極値の差の大きさが所定の閾値以下であることの例に該当する。この場合、点ｂおよびｃが、２つの極点の例に該当する。
「ｙａ≧ｙｃ∧ｙｂ≧ｙｄ」は、点ａが極大点に該当する場合における、２つの極点の何れも、４つの極点における最大値または最小値の何れにも該当しないことの例に該当する。この場合、点ｂおよびｃが、２つの極点の例に該当する。
「ｙａ≧ｙｃ∧ｙｂ≧ｙｄ」は、点ａが極大点に該当する場合における、２つの極点の何れも、４つの極点における最大値または最小値の何れにも該当しないことの例に該当する。この場合、点ｂおよびｃが、２つの極点の例に該当する。
【0039】
「λａｃ≦５」は、内側の２つの極点と４つの極点のうち外側のいずれか１つの極点とに基づいて算出される波長が所定の閾値以下であることの例に該当する。この場合、点ｂおよびｃが、内側の２つの極点の例に該当する。点ａが、４つの極点のうち外側のいずれか１つの極点の例に該当する。
「λｂｄ≦５」は、内側の２つの極点と４つの極点のうち外側のいずれか１つの極点とに基づいて算出される波長が所定の閾値以下であることの例に該当する。この場合、点ｂおよびｃが、内側の２つの極点の例に該当する。点ｄが、４つの極点のうち外側のいずれか１つの極点の例に該当する。
極点除外部１９３が、式（２）が成立する場合に、点ｂおよび点ｃを除外するようにしてもよい。
「λｂｄ≦５」は、内側の２つの極点と４つの極点のうち外側のいずれか１つの極点とに基づいて算出される波長が所定の閾値以下であることの例に該当する。この場合、点ｂおよびｃが、内側の２つの極点の例に該当する。点ｄが、４つの極点のうち外側のいずれか１つの極点の例に該当する。
極点除外部１９３が、式（２）が成立する場合に、点ｂおよび点ｃを除外するようにしてもよい。
【0040】
【数2】
【0041】
「ｈｂｃ／ｈａｄ≦０．２」の「０．２」は、波高の割合の閾値の例に該当する。ただし、波高の割合の閾値は、特定の値に限定されない。
「ｈｂｃ／ｈａｄ≦０．２」は、４つの極点のうち内側の２つの極点における極値の差の大きさの、４つの極点のうち外側の２つの極点における極値の差の大きさに対する割合が所定の閾値以下であることの例に該当する。この場合、点ｂおよびｃが、４つの極点のうち内側の２つの極点の例に該当する。点ａおよびｄが、４つの極点のうち外側の２つの極点の例に該当する。
「ｈｂｃ／ｈａｄ≦０．２」は、４つの極点のうち内側の２つの極点における極値の差の大きさの、４つの極点のうち外側の２つの極点における極値の差の大きさに対する割合が所定の閾値以下であることの例に該当する。この場合、点ｂおよびｃが、４つの極点のうち内側の２つの極点の例に該当する。点ａおよびｄが、４つの極点のうち外側の２つの極点の例に該当する。
【0042】
図５は、測定データにおける、連続する４つの極点の第２例を示す図である。図５のグラフの横軸は上述したｘ軸を示し、縦軸はｙ軸を示す。
図５の例では、点ａおよびｃが極大点に該当し、点ｂおよびｄが極小点に該当する。
図５の例では、点ｃが、４つの極点における最大値に該当し、ｙａ＜ｙｃである。この場合、「ｙａ≧ｙｃ∧ｙｂ≧ｙｄ」が成立しないことから、極点除外部１９３は、上記の式（１）または式（２）の何れに基づいても、点ｂおよびｃを除外しない。
図５の例では、点ａおよびｃが極大点に該当し、点ｂおよびｄが極小点に該当する。
図５の例では、点ｃが、４つの極点における最大値に該当し、ｙａ＜ｙｃである。この場合、「ｙａ≧ｙｃ∧ｙｂ≧ｙｄ」が成立しないことから、極点除外部１９３は、上記の式（１）または式（２）の何れに基づいても、点ｂおよびｃを除外しない。
【0043】
図６は、測定データにおける、連続する４つの極点の第３例を示す図である。図６のグラフの横軸は上述したｘ軸を示し、縦軸はｙ軸を示す。
図６の例では、点ａおよびｃが極大点に該当し、点ｂおよびｄが極小点に該当する。
図６の例では、点ｂが、４つの極点における最小値に該当し、ｙｂ＜ｙｄである。この場合、「ｙａ≧ｙｃ∧ｙｂ≧ｙｄ」が成立しないことから、極点除外部１９３は、上記の式（１）または式（２）の何れに基づいても、点ｂおよびｃを除外しない。
図６の例では、点ａおよびｃが極大点に該当し、点ｂおよびｄが極小点に該当する。
図６の例では、点ｂが、４つの極点における最小値に該当し、ｙｂ＜ｙｄである。この場合、「ｙａ≧ｙｃ∧ｙｂ≧ｙｄ」が成立しないことから、極点除外部１９３は、上記の式（１）または式（２）の何れに基づいても、点ｂおよびｃを除外しない。
【0044】
図７は、測定データにおける、連続する４つの極点の第４例を示す図である。図７のグラフの横軸は上述したｘ軸を示し、縦軸はｙ軸を示す。
図７の例では、点ａおよびｃが極小点に該当し、点ｂおよびｄが極大点に該当する。
図７の例では、点ａおよびｃが極小点に該当し、点ｂおよびｄが極大点に該当する。
【0045】
点ａが極小点の場合、ｈａｄ＝ｙｄ－ｙａと表される。点ａが極小点の場合、点ｂが極大点、点ｃが極小点であり、ｈｂｃ＝ｙｂ－ｙｃと表される。
極点除外部１９３が、式（３）が成立する場合に、点ｂおよび点ｃを除外するようにしてもよい。
極点除外部１９３が、式（３）が成立する場合に、点ｂおよび点ｃを除外するようにしてもよい。
【0046】
【数3】
【0047】
式（３）の「ｈｂｃ≦０．３∨λａｃ≦５∨λｂｄ≦５」は、式（１）の場合と同様である。
「ｙａ≦ｙｃ∧ｙｂ≦ｙｄ」は、点ａが極小点に該当する場合における、２つの極点の何れも、４つの極点における最大値または最小値の何れにも該当しないことの例に該当する。この場合、点ｂおよびｃが、２つの極点の例に該当する。
極点除外部１９３が、式（４）が成立する場合に、点ｂおよび点ｃを除外するようにしてもよい。
「ｙａ≦ｙｃ∧ｙｂ≦ｙｄ」は、点ａが極小点に該当する場合における、２つの極点の何れも、４つの極点における最大値または最小値の何れにも該当しないことの例に該当する。この場合、点ｂおよびｃが、２つの極点の例に該当する。
極点除外部１９３が、式（４）が成立する場合に、点ｂおよび点ｃを除外するようにしてもよい。
【0048】
【数4】
【0049】
式（４）の「ｈｂｃ／ｈａｄ≦０．２」は、式（２）の場合と同様である。
【0050】
図８は、測定データにおける、連続する４つの極点の第５例を示す図である。図８のグラフの横軸は上述したｘ軸を示し、縦軸はｙ軸を示す。
図８の例では、点ａおよびｃが極小点に該当し、点ｂおよびｄが極大点に該当する。
図８の例では、点ｃが、４つの極点における最小値に該当し、ｙａ＞ｙｃである。この場合、「ｙａ≦ｙｃ∧ｙｂ≦ｙｄ」が成立しないことから、極点除外部１９３は、上記の式（３）または式（４）の何れに基づいても、点ｂおよびｃを除外しない。
図８の例では、点ａおよびｃが極小点に該当し、点ｂおよびｄが極大点に該当する。
図８の例では、点ｃが、４つの極点における最小値に該当し、ｙａ＞ｙｃである。この場合、「ｙａ≦ｙｃ∧ｙｂ≦ｙｄ」が成立しないことから、極点除外部１９３は、上記の式（３）または式（４）の何れに基づいても、点ｂおよびｃを除外しない。
【0051】
図９は、測定データにおける、連続する４つの極点の第６例を示す図である。図９のグラフの横軸は上述したｘ軸を示し、縦軸はｙ軸を示す。
図９の例では、点ａおよびｃが極小点に該当し、点ｂおよびｄが極大点に該当する。
図９の例では、点ｂが、４つの極点における最大値に該当し、ｙｂ＞ｙｄである。この場合、「ｙａ≦ｙｃ∧ｙｂ≦ｙｄ」が成立しないことから、極点除外部１９３は、上記の式（３）または式（４）の何れに基づいても、点ｂおよびｃを除外しない。
図９の例では、点ａおよびｃが極小点に該当し、点ｂおよびｄが極大点に該当する。
図９の例では、点ｂが、４つの極点における最大値に該当し、ｙｂ＞ｙｄである。この場合、「ｙａ≦ｙｃ∧ｙｂ≦ｙｄ」が成立しないことから、極点除外部１９３は、上記の式（３）または式（４）の何れに基づいても、点ｂおよびｃを除外しない。
【0052】
極点除外部１９３が、測定データにおける極点のうち外側の極点について、その極点と、その極点に隣接する極点との極値の差の大きさが所定の閾値以下であり、かつ、外側の極点が、その極点を含めて外側から３つの極点における最大値または最小値の何れにも該当しない場合、外側の極点を除外するようにしてもよい。
【0053】
図１０は、測定データにおける、外側の極点の第１例を示す図である。図１０のグラフの横軸は上述したｘ軸を示し、縦軸はｙ軸を示す。
図１０の例では、測定データ全体における極点のうち先頭から順に３つの点が示されている。これら３つの極点を、先頭から順に点ａ、ｂ、ｃと表記する。ここでいうデータの先頭からの順番は、ｘ座標値が小さい順である。
図１０の例では、測定データ全体における極点のうち先頭から順に３つの点が示されている。これら３つの極点を、先頭から順に点ａ、ｂ、ｃと表記する。ここでいうデータの先頭からの順番は、ｘ座標値が小さい順である。
【0054】
点ａは、測定データにおける極点のうち外側の極点の例に該当する。点ａ、ｂ、ｃは、測定データにおける極点のうち外側から３つの極点の例に該当する。図１０の例では、点ａおよびｃが極大点に該当し、点ｂが極小点に該当する。
連続する４つの極点の場合について上述したのと同様、点ａ、ｂ、ｃのｘ座標値を、それぞれｘａ、ｘｂ、ｘｃと表記する。点ａ、ｂ、ｃのｙ座標値を、それぞれｙａ、ｙｂ、ｙｃと表記する。
連続する４つの極点の場合について上述したのと同様、点ａ、ｂ、ｃのｘ座標値を、それぞれｘａ、ｘｂ、ｘｃと表記する。点ａ、ｂ、ｃのｙ座標値を、それぞれｙａ、ｙｂ、ｙｃと表記する。
【0055】
点ａおよび点ｂに基づく波高を波高ｈａｂと表記する。点ａが極大点の場合、ｈａｂ＝ｙａ－ｙｂと表される。波形情報算出部１９４が波高ｈａｂを算出するようにしてもよい。あるいは、極点除外部１９３が、波高ｈａｂを算出するようにしてもよい。
極点除外部１９３が、式（５）が成立する場合に、点ａを除外するようにしてもよい。
極点除外部１９３が、式（５）が成立する場合に、点ａを除外するようにしてもよい。
【0056】
【数5】
【0057】
「ｈａｂ≦０．３」の「０．３」は、０．３メートルを示す。この「０．３」は、波高の閾値の例に該当する定数である。上述したように、波高の閾値は、特定の値に限定されない。
「ｈａｂ≦０．３」は、外側の極点と、その極点に隣接する極点との極値の差の大きさが所定の閾値以下であることの例に該当する。上述したように、点ａが、外側の極点の例に該当する。また、点ｂが、外側の極点に隣接する極点の例に該当する。
「ｙａ≦ｙｃ」は、点ａが極大点に該当する場合における、点ａが、その極点を含めて外側から３つの極点における最大値または最小値の何れにも該当しないことの例に該当する。
「ｈａｂ≦０．３」は、外側の極点と、その極点に隣接する極点との極値の差の大きさが所定の閾値以下であることの例に該当する。上述したように、点ａが、外側の極点の例に該当する。また、点ｂが、外側の極点に隣接する極点の例に該当する。
「ｙａ≦ｙｃ」は、点ａが極大点に該当する場合における、点ａが、その極点を含めて外側から３つの極点における最大値または最小値の何れにも該当しないことの例に該当する。
【0058】
波形情報計算装置が、波形データにおける末尾の極点に対しても、上述した先頭の極点に対する処理と同様の処理を行うようにしてもよい。
図１１は、測定データにおける、外側の極点の第２例を示す図である。図１１のグラフの横軸は上述したｘ軸を示し、縦軸はｙ軸を示す。
図１１の例では、測定データ全体における極点のうち末尾から順に３つの点が示されている。これら３つの極点を、末尾から順に点ａ、ｂ、ｃと表記する。ここでいうデータの末尾からの順番は、ｘ座標値が大きい順である。
図１１は、測定データにおける、外側の極点の第２例を示す図である。図１１のグラフの横軸は上述したｘ軸を示し、縦軸はｙ軸を示す。
図１１の例では、測定データ全体における極点のうち末尾から順に３つの点が示されている。これら３つの極点を、末尾から順に点ａ、ｂ、ｃと表記する。ここでいうデータの末尾からの順番は、ｘ座標値が大きい順である。
【0059】
先頭の極点に対する処理について上述したのと同様、点ａは、測定データにおける極点のうち外側の極点の例に該当する。点ａ、ｂ、ｃは、測定データにおける極点のうち外側から３つの極点の例に該当する。図１１の例では、点ａおよびｃが極大点に該当し、点ｂが極小点に該当する。
【0060】
先頭の極点に対する処理について上述したのと同様、点ａ、ｂ、ｃのｘ座標値を、それぞれｘａ、ｘｂ、ｘｃと表記する。点ａ、ｂ、ｃのｙ座標値を、それぞれｙａ、ｙｂ、ｙｃと表記する。点ａおよび点ｂに基づく波高を波高ｈａｂと表記する。点ａが極大点の場合、点ｃが極大点、点ｄが極小点であり、ｈａｂ＝ｙａ－ｙｂと表される。波形情報算出部１９４が波高ｈｃｄを算出するようにしてもよい。あるいは、極点除外部１９３が、波高ｈｃｄを算出するようにしてもよい。
極点除外部１９３が、上記の式（５）が成立する場合に、点ａを除外するようにしてもよい。
極点除外部１９３が、上記の式（５）が成立する場合に、点ａを除外するようにしてもよい。
【0061】
図１２は、測定データにおける、外側の極点の第３例を示す図である。図１２のグラフの横軸は上述したｘ軸を示し、縦軸はｙ軸を示す。
図１０の場合と同様、図１２の例では、測定データ全体における極点のうち先頭から順に３つの点が示されている。これら３つの極点を、先頭から順に点ａ、ｂ、ｃと表記する。
図１２の例では、点ａおよびｃが極小点に該当し、点ｂが極大点に該当する。
点ａが極小点の場合、ｈａｂ＝ｙｂ－ｙａと表される。
極点除外部１９３が、式（６）が成立する場合に、点ａを除外するようにしてもよい。
図１０の場合と同様、図１２の例では、測定データ全体における極点のうち先頭から順に３つの点が示されている。これら３つの極点を、先頭から順に点ａ、ｂ、ｃと表記する。
図１２の例では、点ａおよびｃが極小点に該当し、点ｂが極大点に該当する。
点ａが極小点の場合、ｈａｂ＝ｙｂ－ｙａと表される。
極点除外部１９３が、式（６）が成立する場合に、点ａを除外するようにしてもよい。
【0062】
【数6】
【0063】
式（６）の「ｈａｂ≦０．３」は、式（５）の場合と同様である。
「ｙａ≧ｙｃ」は、点ａが極小点に該当する場合における、点ａが、その極点を含めて外側から３つの極点における最大値または最小値の何れにも該当しないことの例に該当する。
「ｙａ≧ｙｃ」は、点ａが極小点に該当する場合における、点ａが、その極点を含めて外側から３つの極点における最大値または最小値の何れにも該当しないことの例に該当する。
【0064】
図１３は、測定データにおける、外側の極点の第４例を示す図である。図１３のグラフの横軸は上述したｘ軸を示し、縦軸はｙ軸を示す。
図１１の場合と同様、図１３の例では、測定データ全体における極点のうち末尾から順に３つの点が示されている。これら３つの極点を、末尾から順に点ａ、ｂ、ｃと表記する。
図１３の例では、点ａおよびｃが極小点に該当し、点ｂが極大点に該当する。
極点除外部１９３が、上記の式（６）が成立する場合に、点ｄを除外するようにしてもよい。
図１１の場合と同様、図１３の例では、測定データ全体における極点のうち末尾から順に３つの点が示されている。これら３つの極点を、末尾から順に点ａ、ｂ、ｃと表記する。
図１３の例では、点ａおよびｃが極小点に該当し、点ｂが極大点に該当する。
極点除外部１９３が、上記の式（６）が成立する場合に、点ｄを除外するようにしてもよい。
【0065】
極点除外部１９３が、極大点を中心に連続する５つの極点について、所定の波高算出方法で波高を算出できないか、または、所定の波長算出方法で波長を算出できない場合、その５つの極点の組み合わせを波形情報算出対象から除外するようにしてもよい。
波形情報算出部１９４が、除外されなかった連続する５つの極点に基づいて、中心の極大点における波高および波形を算出するようにしてもよい。
波形情報算出部１９４が、除外されなかった連続する５つの極点に基づいて、中心の極大点における波高および波形を算出するようにしてもよい。
【0066】
図１４は、測定データにおける、連続する５つの極点の例を示す図である。図１４のグラフの横軸は上述したｘ軸を示し、縦軸はｙ軸を示す。
以下では、連続する５つの極点を、ｘ座標値の小さい順に点ａ、ｂ、ｃ、ｄ，ｅと表記する。ｘａ、ｘｂ、ｘｃ、ｘｄ、ｙａ、ｙｂ、ｙｃ、および、ｙｄは、連続する４つの極点の場合について上述したのと同様である。点ｅのｘ座標値をｘｅと表記する。点ｅのｙ座標値をｙｅと表記する。
以下では、連続する５つの極点を、ｘ座標値の小さい順に点ａ、ｂ、ｃ、ｄ，ｅと表記する。ｘａ、ｘｂ、ｘｃ、ｘｄ、ｙａ、ｙｂ、ｙｃ、および、ｙｄは、連続する４つの極点の場合について上述したのと同様である。点ｅのｘ座標値をｘｅと表記する。点ｅのｙ座標値をｙｅと表記する。
【0067】
図１４の例では、点ａ、ｃ、および、ｅが極大点に該当し、点ｂおよびｄが極小点に該当する。点ｃが、中心の極大点に該当する。
波高ｈｂｃ、波高ｈｃｄ、波長λａｃ、および、波長λｂｄは、連続する４つの極点の場合について上述したのと同様である。
点ｃから点ｅまでは、波の１周期に該当し、ｘｅ－ｘｃは、波長に該当する。この波長をλｃｅと表記する。
波高ｈｂｃ、波高ｈｃｄ、波長λａｃ、および、波長λｂｄは、連続する４つの極点の場合について上述したのと同様である。
点ｃから点ｅまでは、波の１周期に該当し、ｘｅ－ｘｃは、波長に該当する。この波長をλｃｅと表記する。
【0068】
極点除外部１９３が、観測データにおける極大点をｘ座標値の小さい順に１つずつ選択し、選択した極大点を中心に連続する５つの極点の組み合わせを選択可能か否かを判定するようにしてもよい。５つの極点の組み合わせを選択できた場合、極点除外部１９３が、さらに、図１４に例示される波高ｈｂｃおよびｈｃｄと、波長λａｃ、λｂｄ、および、λｃｅとを算出可能か否かを判定するようにしてもよい。
【0069】
例えば、極点除外部１９３が選択した極大点が、観測データにおける極大点のうちｘ座標方向の端に位置する極大点である場合、その極大点を中心に連続する５つの極点の組み合わせを選択できないことが考えられる。また、極点除外部１９３が選択した５つの極点の組み合わせに、測定エラーによってｙ座標値を特定できない極点が含まれる場合、上記の波高または波長を算出できないことが考えられる。
【0070】
極点除外部１９３が５つの極点の組み合わせを選択可能、かつ、上記の波高および波長を算出可能と判定した場合、波形情報算出部１９４が、５つの極点のデータに基づいて、中心の極大点における波高および波長を算出するようにしてもよい。図１４の例では、点ｃが、中心の極大点の例に該当する。
【0071】
例えば、波形情報算出部１９４が、波高ｈｂｃおよびｈｃｄのうち何れか大きい方を、中心の極大点における波高として算出するようにしてもよい。この場合、中心の極大点における波高は、ＭＡＸ（ｈａｃ，ｈｃｄ）と表される。「ＭＡＸ」は、引数のうち最大のものの値を出力する関数である。
【0072】
また、波形情報算出部１９４が、波長λａｃ、λｂｄ、および、λｃｅのうち何れか最消のものを、中心の極大点における波長として算出するようにしてもよい。この場合、中心の極大点における波長は、ＭＩＮ（λａｃ，λｂｄ，λｃｅ）と表される。「ＭＩＮ」は、引数のうち最小のものの値を出力する関数である。
【0073】
図１５は、波形情報算出装置１００が波形情報を算出する処理の手順の第１例を示すフローチャートである。
図１５の処理で、データ取得部１９１は、測定データを取得する（ステップＳ１１１）。
極点検出部１９２は、データ取得部１９１が取得した測定データにおける極点を検出する（ステップＳ１１２）。
図１５の処理で、データ取得部１９１は、測定データを取得する（ステップＳ１１１）。
極点検出部１９２は、データ取得部１９１が取得した測定データにおける極点を検出する（ステップＳ１１２）。
【0074】
極点除外部１９３は、極点検出部１９２が検出した極点のうち連続する４点を決定する（ステップＳ１１３）。そして、波形情報算出部１９４は、極点除外部１９３決定した４点における波長および波高を算出する（ステップＳ１１４）。あるいは、極点除外部１９３が波長および波高を算出するようにしてもよい。
極点除外部１９３は、算出された波長および波高に基づいて、所定の条件を満たす極点を除外する（ステップＳ１１５）。図４から図９までを参照して説明したように、連続する４つの極点が式（５）または式（６）を満たす場合に、極点除外部１９３が、４つの極点のうち内側の２つの極点を除外するようにしてもよい。
極点除外部１９３は、算出された波長および波高に基づいて、所定の条件を満たす極点を除外する（ステップＳ１１５）。図４から図９までを参照して説明したように、連続する４つの極点が式（５）または式（６）を満たす場合に、極点除外部１９３が、４つの極点のうち内側の２つの極点を除外するようにしてもよい。
【0075】
そして、極点除外部１９３は、極点を除外する処理の終了条件が成立しているか否かを判定する（ステップＳ１１６）。極点を除外する処理の終了条件は、特定の条件に限定されない。例えば、極点を除外する処理の終了条件は、連続する４点をどのように選んでも削除される極点が無いことであってもよいが、これに限定されない。
【0076】
極点を除外する処理の終了条件が成立していないと極点除外部１９３が判定した場合（ステップＳ１１６：ＮＯ）、処理がステップＳ１１３へ戻る。
一方、極点を除外する処理の終了条件が成立していると判定した場合（ステップＳ１１６：ＹＥＳ）、極点除外部１９３は、測定データ全体における外側の極点に関して所定の条件が成立する場合に、その極点を除外する（ステップＳ１１７）。図１０から図１３までを参照して説明したように、先頭から連続する３つの極点が式（５）または式（６）を満たす場合に、極点除外部１９３が、先頭の極点を除外するようにしてもよい。また、末尾から連続する３つの極点が式（５）または式（６）を満たす場合に、極点除外部１９３が、末尾の極点を除外するようにしてもよい。
一方、極点を除外する処理の終了条件が成立していると判定した場合（ステップＳ１１６：ＹＥＳ）、極点除外部１９３は、測定データ全体における外側の極点に関して所定の条件が成立する場合に、その極点を除外する（ステップＳ１１７）。図１０から図１３までを参照して説明したように、先頭から連続する３つの極点が式（５）または式（６）を満たす場合に、極点除外部１９３が、先頭の極点を除外するようにしてもよい。また、末尾から連続する３つの極点が式（５）または式（６）を満たす場合に、極点除外部１９３が、末尾の極点を除外するようにしてもよい。
【0077】
次に、波形情報算出部１９４は、除外されなかった極点に基づいて波形情報を算出する（ステップＳ１１８）。
波形情報算出装置１００は、算出した波形情報を出力する（ステップＳ１１９）。波形情報算出装置１００が、波形情報を出力する方法は、特定の方法に限定されない。例えば、表示部１２０が制御部１９０の制御に従って波形情報を表示するようにしてもよいが、これに限定されない。
ステップＳ１１９の後、波形情報算出装置１００は、図１５の処理を終了する。
波形情報算出装置１００は、算出した波形情報を出力する（ステップＳ１１９）。波形情報算出装置１００が、波形情報を出力する方法は、特定の方法に限定されない。例えば、表示部１２０が制御部１９０の制御に従って波形情報を表示するようにしてもよいが、これに限定されない。
ステップＳ１１９の後、波形情報算出装置１００は、図１５の処理を終了する。
【0078】
図１６は、波形情報算出装置１００が波形情報を算出する処理の手順の第２例を示すフローチャートである。波形情報算出装置１００が、図１５のステップＳ１１９で図１６の処理を行うようにしてもよい。あるいは、波形情報算出装置１００が、極点検出部１９２が検出する極点に基づいて、図１５の処理とは別に図１６の処理を行うようにしてもよい。
【0079】
図１６の処理で、極点除外部１９３は、測定データにおける極点のうち極大点を１つ選択する（ステップＳ２１１）。例えば、極点除外部１９３が、観測データにおける極大点をｘ座標値の小さい順に１つずつ選択するようにしてもよいが、これに限定されない。
次に、極点除外部１９３は、選択した極大点を中心に連続する５つの極点を選択可能か否かを判定する（ステップＳ２１２）。
次に、極点除外部１９３は、選択した極大点を中心に連続する５つの極点を選択可能か否かを判定する（ステップＳ２１２）。
【0080】
５つの極点を選択不可能と極点除外部１９３が判定した場合（ステップＳ２１２：ＮＯ）、処理がステップＳ２１１へ戻る。
一方、５つの極点を選択可能と極点除外部１９３が判定した場合（ステップＳ２１２：ＹＥＳ）、波形情報算出部１９４は、図１４に例示される波高および波長の算出を試みる（ステップＳ２１３）。そして、波形情報算出部１９４は、波高および波長の算出に成功したか否かを判定する（ステップＳ２１４）。波形情報算出部１９４に代えて極点除外部１９３が、ステップＳ２１３およびＳ２１４の処理を行うようにしてもよい。
一方、５つの極点を選択可能と極点除外部１９３が判定した場合（ステップＳ２１２：ＹＥＳ）、波形情報算出部１９４は、図１４に例示される波高および波長の算出を試みる（ステップＳ２１３）。そして、波形情報算出部１９４は、波高および波長の算出に成功したか否かを判定する（ステップＳ２１４）。波形情報算出部１９４に代えて極点除外部１９３が、ステップＳ２１３およびＳ２１４の処理を行うようにしてもよい。
【0081】
波高および波長の算出に失敗したと波形情報算出部１９４が判定した場合（ステップＳ２１４：ＮＯ）、処理がステップＳ２１１へ戻る。
一方、波高および波長の算出に成功したと判定した場合（ステップＳ２１４：ＹＥＳ）、波形情報算出部１９４は、ステップＳ２１３で算出した波高および波長に基づいて、中心の極大点における波高および波長を算出する（ステップＳ２１５）。
一方、波高および波長の算出に成功したと判定した場合（ステップＳ２１４：ＹＥＳ）、波形情報算出部１９４は、ステップＳ２１３で算出した波高および波長に基づいて、中心の極大点における波高および波長を算出する（ステップＳ２１５）。
【0082】
次に、極点除外部１９３は、波形情報算出の終了条件が成立しているか否かを判定する（ステップＳ２１６）。波形情報算出の終了条件は、特定の条件に限定されない。例えば、極点除外部１９３が、測定データにおける全ての極大点についてステップＳ２１１からＳ２１５までの処理を行ったか否かを判定するようにしてもよいが、これに限定されない。
【0083】
終了条件が成立していないと極点除外部１９３が判定した場合（ステップＳ２１６：ＮＯ）、処理がステップＳ２１１へ戻る。
一方、終了条件が成立していると極点除外部１９３が判定した場合（ステップＳ２１６：ＹＥＳ）、波形情報算出装置１００は、図１６の処理を終了する。
図１５の処理とは別に図１６の処理を行う場合、波形情報算出装置１００は、図１６の処理の終了後、算出した波形情報を出力する。波形情報算出装置１００が、波形情報を出力する方法は、特定の方法に限定されない。例えば、表示部１２０が制御部１９０の制御に従って波形情報を表示するようにしてもよいが、これに限定されない。
一方、終了条件が成立していると極点除外部１９３が判定した場合（ステップＳ２１６：ＹＥＳ）、波形情報算出装置１００は、図１６の処理を終了する。
図１５の処理とは別に図１６の処理を行う場合、波形情報算出装置１００は、図１６の処理の終了後、算出した波形情報を出力する。波形情報算出装置１００が、波形情報を出力する方法は、特定の方法に限定されない。例えば、表示部１２０が制御部１９０の制御に従って波形情報を表示するようにしてもよいが、これに限定されない。
【0084】
図１７は、制御部１９０の構成例を示す概略ブロック図である。図１７は、波形情報算出装置１００が飛行艇に搭載される場合の例を示している。
図１７に示す構成では、制御部１９０は、図１に示される極点検出部１９２と、極点除外部１９３と、波形情報算出部１９４とに加えて、加速度補正部２０１と、ピッチ・ロール補正部２０２と、ローパスフィルタ２０３と、測定結果画面出力部２０４とを備える。
図１７は、データ取得部１９１が制御部１９０の外部の構成となっている場合の例を示している。あるいは、制御部１９０が、図１７の各部に加えてさらにデータ取得部１９１を備えるようにしてもよい。
図１７に示す構成では、制御部１９０は、図１に示される極点検出部１９２と、極点除外部１９３と、波形情報算出部１９４とに加えて、加速度補正部２０１と、ピッチ・ロール補正部２０２と、ローパスフィルタ２０３と、測定結果画面出力部２０４とを備える。
図１７は、データ取得部１９１が制御部１９０の外部の構成となっている場合の例を示している。あるいは、制御部１９０が、図１７の各部に加えてさらにデータ取得部１９１を備えるようにしてもよい。
【0085】
加速度補正部２０１は、飛行艇の加速度に応じて測定データを補正する。例えば、加速度補正部２０１が、飛行艇の加減速にかかわらずｘ座標方向におけるサンプル点の間隔が一定になるようにサンプル点の間引きおよび補間を行うようにしてもよい。
ピッチ・ロール補正部２０２は、飛行艇の機体の傾斜に応じて測定データを補正する。例えば、ピッチ・ロール補正部２０２が、飛行艇の機体の傾斜に応じて測定データのｙ座標値を補正することで、機体が傾斜していない場合と同様のデータを得られるようにしてもよい。
ピッチ・ロール補正部２０２は、飛行艇の機体の傾斜に応じて測定データを補正する。例えば、ピッチ・ロール補正部２０２が、飛行艇の機体の傾斜に応じて測定データのｙ座標値を補正することで、機体が傾斜していない場合と同様のデータを得られるようにしてもよい。
【0086】
ローパスフィルタ２０３は、測定データにおける不要な高周波成分を除去する。制御部１９０が、ローパスフィルタ２０３のカットオフ周波数（Fc；Cutoff Frequency）を、飛行艇の速度と計測に不要な最小波長から決定するようにしてもよい。
また、制御部１９０が、波形の測定に要求される精度に応じて、ローパスフィルタ２０３の適用による波形の歪みを抑制するように、カットオフ周波数を設定する、あるいは、ローパスフィルタ２０３をオフにするようにしてもよい。
また、制御部１９０が、波形の測定に要求される精度に応じて、ローパスフィルタ２０３の適用による波形の歪みを抑制するように、カットオフ周波数を設定する、あるいは、ローパスフィルタ２０３をオフにするようにしてもよい。
【0087】
測定結果画面出力部２０４は、波形の測定結果を表示部１２０に表示させる。例えば、測定結果画面出力部２０４が、測定データが示す波形、および、波形情報算出部１９４が計算する波形情報を表示部１２０に表示させるようにしてもよい。
【0088】
図１８は、極点除外部１９３が除外する極値の例を示す図である。図１８のグラフの横軸はｘ座標を示し、縦軸はｙ座標を示す。
線Ｌ２１は、ローパスフィルタ２０３を適用する前の波形の例を示す。線Ｌ２２は、ローパスフィルタ２０３適用後の波形の例を示す。
線Ｌ２１は、ローパスフィルタ２０３を適用する前の波形の例を示す。線Ｌ２２は、ローパスフィルタ２０３適用後の波形の例を示す。
【0089】
ローパスフィルタ２０３適用後の波形の極点のうち、極点除外部１９３が除外する極点を四角で示し、残りの極点（除外されない極点）を丸で示している。
ローパスフィルタ２０３の適用、および、極点除外部１９３による極点の除外によって、波形における局所的な情報にとらわれて波の大きさを過少に評価する可能性を低減させることができる。
ローパスフィルタ２０３の適用、および、極点除外部１９３による極点の除外によって、波形における局所的な情報にとらわれて波の大きさを過少に評価する可能性を低減させることができる。
【0090】
以上のように、データ取得部１９１は、波形の測定データを取得する。極点検出部１９２は、測定データにおける極点を検出する。極点除外部１９３は、複数の極点の位置関係に基づいて、所定の条件を満たす極点を除外する。波形情報算出部１９４は、除外されなかった極点に基づいて、波形における波高または波長の少なくとも何れかを算出する。
波形情報算出装置１００では、波形における局所的な情報にとらわれて波の大きさを過少に評価する可能性を低減させ、かつ、波形の変化（歪み）による誤差を回避するか、あるいは低減させて、波形情報を算出できると期待される。波形情報算出装置１００によれば、この点で、波形情報を高精度に算出することができる。
波形情報算出装置１００では、波形における局所的な情報にとらわれて波の大きさを過少に評価する可能性を低減させ、かつ、波形の変化（歪み）による誤差を回避するか、あるいは低減させて、波形情報を算出できると期待される。波形情報算出装置１００によれば、この点で、波形情報を高精度に算出することができる。
【0091】
極点除外部１９３による極値選択処理を連続的に実施することにより、波形情報算出装置１００側の動揺や海面の反射点の変動等の影響を軽減すると共に、高速かつ波高を損なうことなく波高データを計測することができる。また、極点除外部１９３が極値を選択する（選択しない極値を除外する）ことで、ローパスフィルタ適用による波高データのなまりを回避または軽減することができ、波形情報を高精度に算出することができる。
また、極点除外部１９３が行う処理は、フィルター処理のような全データを対象とせず、極値に絞った処理であり、処理対象データが極めて少なくなる為、高速な処理が可能となる。
また、極点除外部１９３が行う処理は、フィルター処理のような全データを対象とせず、極値に絞った処理であり、処理対象データが極めて少なくなる為、高速な処理が可能となる。
【0092】
また、極点除外部１９３は、連続する４つの極点のうち内側の２つの極点について、これら２つの極点における極値の差、または、少なくともこれら２つの極点に基づいて算出される波長の少なくとも何れかに基づいて、これら２つの極点を除外するか否かを決定する。
極点除外部１９３は、極値の差、および、極点のｘ座標の差に基づく波長を用いて極点を除外するか否かを決定することができ、この点で、極点除外部１９３の負荷が軽くて済む。
極点除外部１９３は、極値の差、および、極点のｘ座標の差に基づく波長を用いて極点を除外するか否かを決定することができ、この点で、極点除外部１９３の負荷が軽くて済む。
【0093】
また、極点除外部１９３は、４つの極点のうち内側の２つの極点について、これら２つの極点における極値の差の大きさが所定の閾値以下であり、かつ、これら２つの極点の何れも、４つの極点における最大値または最小値の何れにも該当しない場合、これら２つの極点を除外する。
極点除外部１９３は、極値の差の大きさと閾値との比較、および、極点が４つの極点における最大値または最小値か否かの判定といった比較的簡単な処理で極点を除外するか否かを決定することができ、この点で、極点除外部１９３の負荷が軽くて済む。
極点除外部１９３は、極値の差の大きさと閾値との比較、および、極点が４つの極点における最大値または最小値か否かの判定といった比較的簡単な処理で極点を除外するか否かを決定することができ、この点で、極点除外部１９３の負荷が軽くて済む。
【0094】
また、極点除外部１９３は、４つの極点のうち内側の２つの極点における極値の差の大きさの、４つの極点のうち外側の２つの極点における極値の差の大きさに対する割合が所定の閾値以下であり、かつ、内側の２つの極点の何れも、４つの極点における最大値または最小値の何れにも該当しない場合、内側の２つの極点を除外する。
極点除外部１９３は、極値の差の大きさの割合と閾値との比較、および、極点が４つの極点における最大値または最小値か否かの判定といった比較的簡単な処理で極点を除外するか否かを決定することができ、この点で、極点除外部１９３の負荷が軽くて済む。
極点除外部１９３は、極値の差の大きさの割合と閾値との比較、および、極点が４つの極点における最大値または最小値か否かの判定といった比較的簡単な処理で極点を除外するか否かを決定することができ、この点で、極点除外部１９３の負荷が軽くて済む。
【0095】
また、極点除外部１９３は、４つの極点のうち内側の２つの極点について、これら２つの極点と４つの極点のうち外側のいずれか１つの極点とに基づいて算出される波長が所定の閾値以下であり、かつ、内側の２つの極点の何れも、４つの極点における最大値または最小値の何れにも該当しない場合、内側の２つの極点を除外する。
極点除外部１９３は、極点のｘ座標に基づく波長と閾値との比較、および、極点が４つの極点における最大値または最小値か否かの判定といった比較的簡単な処理で極点を除外するか否かを決定することができ、この点で、極点除外部１９３の負荷が軽くて済む。
極点除外部１９３は、極点のｘ座標に基づく波長と閾値との比較、および、極点が４つの極点における最大値または最小値か否かの判定といった比較的簡単な処理で極点を除外するか否かを決定することができ、この点で、極点除外部１９３の負荷が軽くて済む。
【0096】
また、極点除外部１９３は、連続する４つの極点のうち外側の極点について、外側の極点と、その極点に隣接する極点との極値の差の大きさが所定の閾値以下であり、かつ、外側の極点が、４つの極点における最大値または最小値の何れにも該当しない場合、外側の極点を除外する。
極点除外部１９３は、極値の差の大きさと閾値との比較、および、極点が４つの極点における最大値または最小値か否かの判定といった比較的簡単な処理で極点を除外するか否かを決定することができ、この点で、極点除外部１９３の負荷が軽くて済む。
極点除外部１９３は、極値の差の大きさと閾値との比較、および、極点が４つの極点における最大値または最小値か否かの判定といった比較的簡単な処理で極点を除外するか否かを決定することができ、この点で、極点除外部１９３の負荷が軽くて済む。
【0097】
また、極点除外部１９３は、極大点を中心に連続する５つの極点について、所定の波高算出方法で波高を算出できないか、または、所定の波長算出方法で波長を算出できない場合、その５つの極点の組み合わせを波形情報算出対象から除外する。波形情報算出部１９４は、除外されなかった連続する５つの極点に基づいて、中心の極大点における波高および波形を算出する。
波形情報算出装置１００によれば、波形データの極大点ごとに波高および波長を算出することができる。
波形情報算出装置１００によれば、波形データの極大点ごとに波高および波長を算出することができる。
【0098】
図１９は、実施形態に係る波形情報算出装置の構成例を示す図である。図１９に示す波形情報算出装置６１０は、データ取得部６１１と、極点検出部６１２と、極点除外部６１３と、波形情報算出部６１４とを備える。
かかる構成で、データ取得部６１１は、波形の測定データを取得する。極点検出部６１２は、測定データにおける極点を検出する。極点除外部６１３は、複数の極点の位置関係に基づいて、所定の条件を満たす極点を除外する。波形情報算出部６１４は、除外されなかった極点に基づいて、波形における波高または波長の少なくとも何れかを算出する。
かかる構成で、データ取得部６１１は、波形の測定データを取得する。極点検出部６１２は、測定データにおける極点を検出する。極点除外部６１３は、複数の極点の位置関係に基づいて、所定の条件を満たす極点を除外する。波形情報算出部６１４は、除外されなかった極点に基づいて、波形における波高または波長の少なくとも何れかを算出する。
【0099】
データ取得部６１１は、データ取得手段の例に該当する。極点検出部６１２は、極点検出手段の例に該当する。極点除外部６１３は、極点除外手段の例に該当する。波形情報算出部６１４は、波形情報算出手段の例に該当する。
波形情報算出装置６１０では、波形における局所的な情報にとらわれて波の大きさを過少に評価する可能性を低減させ、かつ、波形の変化（歪み）による誤差を回避するか、あるいは低減させて、波形情報を算出できると期待される。波形情報算出装置６１０によれば、この点で、波形情報を高精度に算出することができる。
波形情報算出装置６１０では、波形における局所的な情報にとらわれて波の大きさを過少に評価する可能性を低減させ、かつ、波形の変化（歪み）による誤差を回避するか、あるいは低減させて、波形情報を算出できると期待される。波形情報算出装置６１０によれば、この点で、波形情報を高精度に算出することができる。
【0100】
図２０は、実施形態に係る波形情報算出方法における処理の手順の例を示すフローチャートである。
図２０に示す処理は、データを取得すること（ステップＳ６１１）と、極点を検出すること（ステップＳ６１２）と、極点を除外すること（ステップＳ６１３）と、波形情報を算出すること（ステップＳ６１４）とを含む。
データを取得すること（ステップＳ６１１）と、極点を検出すること（ステップＳ６１２）と、極点を除外すること（ステップＳ６１３）と、波形情報を算出すること（ステップＳ６１４）とを含む。
図２０に示す処理は、データを取得すること（ステップＳ６１１）と、極点を検出すること（ステップＳ６１２）と、極点を除外すること（ステップＳ６１３）と、波形情報を算出すること（ステップＳ６１４）とを含む。
データを取得すること（ステップＳ６１１）と、極点を検出すること（ステップＳ６１２）と、極点を除外すること（ステップＳ６１３）と、波形情報を算出すること（ステップＳ６１４）とを含む。
【0101】
データを取得すること（ステップＳ６１１）では、波形の測定データを取得する。極点を検出すること（ステップＳ６１２）では、測定データにおける極点を検出する。極点を除外すること（ステップＳ６１３）では、複数の極点の位置関係に基づいて、所定の条件を満たす極点を除外する。波形情報を算出すること（ステップＳ６１４）では、除外されなかった極点に基づいて、波形における波高または波長の少なくとも何れかを算出する。
図２０に示す処理では、波形における局所的な情報にとらわれて波の大きさを過少に評価する可能性を低減させ、かつ、波形の変化（歪み）による誤差を回避するか、あるいは低減させて、波形情報を算出できると期待される。図２０に示す処理によれば、この点で、波形情報を高精度に算出することができる。
図２０に示す処理では、波形における局所的な情報にとらわれて波の大きさを過少に評価する可能性を低減させ、かつ、波形の変化（歪み）による誤差を回避するか、あるいは低減させて、波形情報を算出できると期待される。図２０に示す処理によれば、この点で、波形情報を高精度に算出することができる。
【0102】
図２１は、少なくとも１つの実施形態に係るコンピュータの構成を示す概略ブロック図である。
図７に示す構成で、コンピュータ７００は、ＣＰＵ（Central Processing Unit、中央処理装置）７１０と、主記憶装置７２０と、補助記憶装置７３０と、インタフェース７４０とを備える。
図７に示す構成で、コンピュータ７００は、ＣＰＵ（Central Processing Unit、中央処理装置）７１０と、主記憶装置７２０と、補助記憶装置７３０と、インタフェース７４０とを備える。
【0103】
上記の波形情報算出装置１００および波形情報算出装置６１０のうち何れか１つ以上が、コンピュータ７００に実装されてもよい。その場合、上述した各処理部の動作は、プログラムの形式で補助記憶装置７３０に記憶されている。ＣＰＵ７１０は、プログラムを補助記憶装置７３０から読み出して主記憶装置７２０に展開し、当該プログラムに従って上記処理を実行する。また、ＣＰＵ７１０は、プログラムに従って、上記処理に用いられる記憶領域を主記憶装置７２０に確保する。
【0104】
波形情報算出装置１００がコンピュータ７００に実装される場合、制御部１９０およびその各部の動作は、プログラムの形式で補助記憶装置７３０に記憶されている。ＣＰＵ７１０は、プログラムを補助記憶装置７３０から読み出して主記憶装置７２０に展開し、当該プログラムに従って上記処理を実行する。
また、ＣＰＵ７１０は、プログラムに従って、記憶部１８０に対応する記憶領域を主記憶装置７２０に確保する。
また、ＣＰＵ７１０は、プログラムに従って、記憶部１８０に対応する記憶領域を主記憶装置７２０に確保する。
【0105】
通信部１１０による通信は、インタフェース７４０が通信機能を有し、ＣＰＵ７１０の制御に従って通信を行うことで実行される。
表示部１２０による表示は、インタフェース７４０が表示画面を備え、ＣＰＵ７１０の制御に従って各種画像を表示することで実行される。操作入力部１３０によるユーザ操作の受け付けは、インタフェース７４０が入力デバイスを備えてユーザ操作を受け付け、受け付けたユーザ操作を示す信号をＣＰＵ７１０に出力することで実行される。
表示部１２０による表示は、インタフェース７４０が表示画面を備え、ＣＰＵ７１０の制御に従って各種画像を表示することで実行される。操作入力部１３０によるユーザ操作の受け付けは、インタフェース７４０が入力デバイスを備えてユーザ操作を受け付け、受け付けたユーザ操作を示す信号をＣＰＵ７１０に出力することで実行される。
【0106】
波形情報算出装置６１０がコンピュータ７００に実装される場合、データ取得部６１１と、極点検出部６１２と、極点除外部６１３と、波形情報算出部６１４との動作は、プログラムの形式で補助記憶装置７３０に記憶されている。ＣＰＵ７１０は、プログラムを補助記憶装置７３０から読み出して主記憶装置７２０に展開し、当該プログラムに従って上記処理を実行する。
また、データ取得部６１１が波形の測定データを取得する処理は、例えば、インタフェース７４０が通信機能を有し、ＣＰＵ７１０の制御に従ってセンサなど他の装置と通信を行うことで実行される。
また、データ取得部６１１が波形の測定データを取得する処理は、例えば、インタフェース７４０が通信機能を有し、ＣＰＵ７１０の制御に従ってセンサなど他の装置と通信を行うことで実行される。
【0107】
なお、波形情報算出装置１００、および、波形情報算出装置６１０の全部または一部の機能を実現するためのプログラムをコンピュータ読み取り可能な記録媒体に記録して、この記録媒体に記録されたプログラムをコンピュータシステムに読み込ませ、実行することにより各部の処理を行ってもよい。ここでいう「コンピュータシステム」とは、ＯＳ（オペレーティングシステム）や周辺機器等のハードウェアを含む。
「コンピュータ読み取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ＲＯＭ（Read Only Memory）、ＣＤ－ＲＯＭ（Compact Disc Read Only Memory）等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。また上記プログラムは、前述した機能の一部を実現するためのものであってもよく、さらに前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるものであってもよい。
「コンピュータ読み取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ＲＯＭ（Read Only Memory）、ＣＤ－ＲＯＭ（Compact Disc Read Only Memory）等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。また上記プログラムは、前述した機能の一部を実現するためのものであってもよく、さらに前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるものであってもよい。
【0108】
以上、本発明の実施形態について図面を参照して詳述してきたが、具体的な構成はこの実施形態に限られるものではなく、この発明の要旨を逸脱しない範囲の設計変更等も含まれる。
【符号の説明】
【0109】
１００、６１０ 波形情報算出装置
１１０ 通信部
１２０ 表示部
１３０ 操作入力部
１８０ 記憶部
１９０ 制御部
１９１、６１１ データ取得部
１９２、６１２ 極点検出部
１９３、６１３ 極点除外部
１９４、６１４ 波形情報算出部
２０１ 加速度補正部
２０２ ピッチ・ロール補正部
２０３ ローパスフィルタ
２０４ 測定結果画面出力部
１１０ 通信部
１２０ 表示部
１３０ 操作入力部
１８０ 記憶部
１９０ 制御部
１９１、６１１ データ取得部
１９２、６１２ 極点検出部
１９３、６１３ 極点除外部
１９４、６１４ 波形情報算出部
２０１ 加速度補正部
２０２ ピッチ・ロール補正部
２０３ ローパスフィルタ
２０４ 測定結果画面出力部
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【図14】
【図15】
【図16】
【図17】
【図18】
【図19】
【図20】
【図21】