JP4719373B2 - Wave height calculation device, wave height calculation method, and computer-readable medium - Google Patents

Wave height calculation device, wave height calculation method, and computer-readable medium Download PDF

Info

Publication number
JP4719373B2
JP4719373B2 JP2001124160A JP2001124160A JP4719373B2 JP 4719373 B2 JP4719373 B2 JP 4719373B2 JP 2001124160 A JP2001124160 A JP 2001124160A JP 2001124160 A JP2001124160 A JP 2001124160A JP 4719373 B2 JP4719373 B2 JP 4719373B2
Authority
JP
Japan
Prior art keywords
wave height
calculating
calculated
water surface
average period
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.)
Expired - Lifetime
Application number
JP2001124160A
Other languages
Japanese (ja)
Other versions
JP2002318115A (en
Inventor
次清 平山
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
JAPAN SHIP TECHNOLOGY RESEARCH ASSOCIATION
IHI Corp
Mitsubishi Heavy Industries Ltd
Mitsui Engineering and Shipbuilding Co Ltd
Nippon Yusen KK
Sumitomo Heavy Industries Ltd
Japan Marine Science Inc
Shin Kurushima Dockyard Co Ltd
Universal Shipbuilding Corp
Mitsui E&S Holdings Co Ltd
Kawasaki Motors Ltd
Original Assignee
JAPAN SHIP TECHNOLOGY RESEARCH ASSOCIATION
IHI Corp
Mitsubishi Heavy Industries Ltd
Mitsui Engineering and Shipbuilding Co Ltd
Nippon Yusen KK
Sumitomo Heavy Industries Ltd
Japan Marine Science Inc
Kawasaki Jukogyo KK
Shin Kurushima Dockyard Co Ltd
Universal Shipbuilding Corp
Mitsui E&S Holdings Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by JAPAN SHIP TECHNOLOGY RESEARCH ASSOCIATION, IHI Corp, Mitsubishi Heavy Industries Ltd, Mitsui Engineering and Shipbuilding Co Ltd, Nippon Yusen KK, Sumitomo Heavy Industries Ltd, Japan Marine Science Inc, Kawasaki Jukogyo KK, Shin Kurushima Dockyard Co Ltd, Universal Shipbuilding Corp, Mitsui E&S Holdings Co Ltd filed Critical JAPAN SHIP TECHNOLOGY RESEARCH ASSOCIATION
Priority to JP2001124160A priority Critical patent/JP4719373B2/en
Publication of JP2002318115A publication Critical patent/JP2002318115A/en
Application granted granted Critical
Publication of JP4719373B2 publication Critical patent/JP4719373B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Images

Classifications

    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30—Assessment of water resources

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Description

【０００１】
【発明の属する技術分野】
本発明は、水面に生じている波の高さを、船舶用レーダーを利用して算出するための波高算出装置及び波高算出方法と、コンピュータを、波高算出装置として動作させるためのコンピュータ可読媒体とに、関する。
【０００２】
【従来の技術】
水面（主として海面）に生じている波の高さ（以下、波高と表記する）に関する統計的及び／或いはリアルタイムな情報は、船舶及び浮動構造の設計に役立つだけではなく、船舶等の保守と船舶のウェザールーティングにも役立ち、当該情報に基づき、操船の安全性及び効率を、向上させることが出来る。
【０００３】
このため、波高を測定するためのさまざまな方法が開発されており、その一つとして、船舶用レーダーを利用して波高を測定（算出）する方法が知られている。
【０００４】
従来の、船舶用レーダーを利用した波高測定方法では、船舶レーダーにより得られた情報に基づきシグナル雑音比（スペクトルピークをバックグラウンド雑音で割った値：ＳＮＲ）が算出され、算出されたＳＮＲから波高が推定されている。
【０００５】
具体的には、（１）式により、有義波高Ｈ1/3（一群の波の中の高い方から数えて全体の１／３の間にある波についての平均の波高を持つ波の高さ）が求められている。
【０００６】
【数１】

Figure 0004719373
なお、α、βは、（１）式により有義波高が得られるようにするために、波浪ブイを利用した方法等の別の方法による測定結果に基づき、それらの値が決定される較正用パラメタである。
【０００７】
【発明が解決しようとする課題】
上記した従来の波高算出手順を用いた場合、比較的に大きな波が生じているときには波高が正確に推定できるが、小さな波が生じているときに、波高が正確に推定できない（測定精度が低くなる）といった問題があった。
【０００８】
そこで、本発明の課題は、より小さな波の波高も精度良く算出（推定）できる船舶用レーダーを利用する波高算出装置及び波高算出方法と、コンピュータを、そのような波高算出装置として動作させるためのコンピュータ可読媒体とを、提供することにある。
【０００９】
【課題を解決するための手段】
上記課題を解決するために、本発明の波高算出装置は、レーダー画像信号に基づき、波高を算出すべき水面領域に相当するレーダー画像上の、信号レベルが予め定められている範囲内のレベルとなっている陰影部の面積を当該水面領域の面積で割った値である陰影比を算出する陰影比算出手段と、この陰影比算出手段によって算出された陰影比を、所定の演算式に代入することにより波高を算出する波高算出手段と、レーダー画像信号に基づき、波高を算出すべき水面領域に生じている波の平均周期を算出する平均周期算出手段とを備え、前記所定の演算式は、波高を算出すべき水面領域に生じている波の平均周期の関数として定められている所定数のパラメータを含み、前記波高算出手段は、前記平均周期算出手段によって算出された平均周期に基づき前記所定数のパラメータの値を求め、求めた値を各パラメータの値とした前記所定の演算式に前記陰影比を代入することにより波高を算出する構成を有する。
【００１０】
また、本発明の波高算出方法では、レーダー画像信号に基づき、波高を算出すべき水面領域に相当するレーダー画像上の、信号レベルが予め定められている範囲内のレベルとなっている陰影部の面積を当該水面領域の面積で割った値である陰影比を算出する陰影比算出ステップと、この陰影比算出ステップによって算出された陰影比を、波高を算出すべき水面領域に生じている波の平均周期の関数として定められている所定数のパラメータを含む所定の演算式に代入することにより波高を算出する波高算出ステップであって、レーダー画像信号に基づき、波高を算出すべき水面領域に生じている波の平均周期を算出し、算出した平均周期に基づき前記所定数のパラメータの値を求め、求めた値を各パラメータの値とした前記所定の演算式に前記陰影比を代入することにより波高を算出する波高算出ステップとにより、波高が算出される。
【００１１】
そして、本発明のコンピュータ可読媒体は、レーダー装置と接続可能なコンピュータを、レーダー画像信号に基づき、波高を算出すべき水面領域に相当するレーダー画像上の、信号レベルが予め定められている範囲内のレベルとなっている陰影部の面積を当該水面領域の面積で割った値である陰影比を算出する陰影比算出手段と、この陰影比算出手段によって算出された陰影比を、波高を算出すべき水面領域に生じている波の平均周期の関数として定められている所定数のパラメータを含む所定の演算式に代入することにより波高を算出する波高算出手段であって、レーダー画像信号に基づき、波高を算出すべき水面領域に生じている波の平均周期を算出し、算出した平均周期に基づき前記所定数のパラメータの値を求め、求めた値を各パラメータの値とした前記所定の演算式に前記陰影比を代入することにより波高を算出する波高算出手段とを、備える装置として動作させるためのプログラムを記録して成る。
【００１２】
すなわち、上記課題を解決するために、鋭意、研究を行なったところ、発明者らは、波高の指標として陰影比を用いると、小さな波の波高をも正確に算出できることを見出した。従って、この原理に基づき波高が算出される本発明の波高算出装置、波高算出方法によれば、従来よりも正確な波高の算出が行なえることになる。また、本発明のコンピュータ可読媒体によれば、レーダー装置と接続可能なコンピュータを、本発明の波高算出装置として機能させることが出来ることになる。
【００１５】
また、本発明の波高算出装置を実現するに際しては、所定の演算式として、陰影比ＳＲと、波高を算出すべき水面領域に生じている波の平均周期の関数として定められている２つのパラメータＡ及びＢとから、波高Ｈを算出する「Ｈ＝Ａ・ＳＲ B 」という式を採用することが出来る。
【００１６】
【発明の実施の形態】
以下、本発明の実施の形態を、図面を参照して詳細に説明する。
【００１７】
図１に、本発明の一実施形態に係る波高算出装置１０の使用形態を示す。
【００１８】
図示したように、本実施形態に係る波高算出装置１０は、インターフェースモジュール５１を介してレーダーユニット５２に接続された状態で使用される。レーダーユニット５２は、送信機、レーダーアンテナ、受信機等からなるシステムである。インターフェースモジュール５１は、レーダーユニット５２の受信機が出力する信号（レーダー画像信号）を、レーダー指示器５３と波高算出装置１０とに、供給するモジュールである。
【００１９】
レーダー指示器５３は、入力されたレーダー画像信号に応じたレーダー画像を表示する機器であり、レーダーユニット５２と直結することも出来る機器となっている。換言すれば、レーダーユニット５２とレーダー指示器５３とからなる部分が、一般に船舶用レーダーと呼ばれているものとなっており、インターフェースモジュール５１は、船舶用レーダー（レーダーユニット５２）から、レーダー画像信号を取り出すためのモジュールとなっている。
【００２０】
波高算出装置１０は、波高算出プログラム２５がＣＤ−ＲＯＭ２０等からインストールされているコンピュータであり、インターフェースモジュール５１からレーダー画像信号を取得するためのインタフェースと、ＧＰＳ信号，船速信号及び進路信号を、図示せぬ操舵装置から取得するためのインタフェースとを、備えている。また、波高算出装置１０は、波高表示装置１１に接続するためのインタフェースを備える。波高表示装置１１は、波高等が表示される液晶ディスプレイを備えた機器であり、波高表示装置１１を接続しなくとも、波高算出装置１０は機能する。なお、波高算出装置１０に接続された波高表示装置１１は、通常、レーダー指示器５３の近傍に、設置される。
【００２１】
以下、波高算出装置１０の動作を説明する。
【００２２】
波高算出装置１０は、電源が投入されると、レーダー画像が表示された領域を含む画面をディスプレイに表示して、波高算出条件がオペレータによって設定されるのを待機する状態を取る。ここで、波高算出条件とは、どの領域のどの時点における波高を算出するかを、オペレータが波高算出装置１０に指示するための情報のことであり、波高算出条件は、波高を複数回算出することを指定できる（波高の算出を行なう時刻及び領域を、予め複数組、指定しておくことができる）ものとなっている。
【００２３】
なお、本波高算出装置１０において、波高の算出対象となる領域（以下、観察対象領域と表記する）は、図２に示してあるように、一方の平行な２辺の中心線が、レーダーアンテナ上（レーダー画像の中心）を通る正方形領域となっており、波高算出条件設定時には、オペレータによって、この観察対象領域の指定を含む作業が行なわれる。
【００２４】
図３に、波高算出条件の設定が完了した後に所定の指示（波高算出処理の開始指示）が入力されたときに、波高算出装置１０が実行する波高算出処理の流れ図を示す。
【００２５】
図示したように、波高算出処理開始時、波高算出装置１０は、設定されている波高算出条件にて、波高を算出すべき時刻として指定されている算出時刻となるまで待機する（ステップS１０１）。そして、波高算出装置１０は、算出時刻となった時にステップS１０１を終了し、レーダーユニット５２が出力するレーダー画像信号を、インターフェースモジュール５１を介して、波高算出条件にて規定されている量、収集して、内蔵しているＨＤＤに記憶する（ステップＳ１０２）。このステップＳ１０２において、波高算出装置１０は、例えば、３２(＝２5)個のレーダー画像分のレーダー画像信号を収集する。その後、波高算出装置１０は、収集したレーダー画像信号に含まれる、波高算出条件にて規定されている観察対象領域に関する信号に対して三次元ＦＦＴ(Fast Fourier Transform）等を行なうことにより、ノイズ（波浪に関係しない情報）を除去した波数・周波数スペクトルを生成する（ステップＳ１０３）。
【００２６】
次いで、波高算出装置１０は、生成した波数・周波数スペクトルから平均波周期を求める（ステップＳ１０４）。さらに、波高算出装置１０は、波数・周波数スペクトルの逆ＦＦＴを行なうことにより、ノイズが除去されたレーダー画像データを生成し、当該レーダー画像データに基づき、陰影比ＳＲを算出した後、平均波周期の関数として用意されているパラメータＡ、Ｂが用いられた関係式Ｈ1/3＝Ａ・ＳＲBに、算出したＳＲを代入することにより有義波高Ｈ1/3の算出を行なう（ステップＳ１０５）。
【００２７】
具体的には、このステップＳ１０５において、波高算出装置１０は、レーダー画像の、低いレベルでその強度がほぼ一定値となっている部分（以下、陰影部と表記する）の面積を算出し、算出した陰影部面積を観察対象領域の全面積で割った値を、陰影比ＳＲとして算出している。また、陰影部の判定条件としては、波浪に関する各種の値（含む有義波高Ｈ1/3）を仮定してシミュレーションにより各種のレーダー画像データを作成し、作成したレーダー画像データから算出される陰影比ＳＲが、仮定した有義波高Ｈ1/3と一致するように定めたものが用いられている。なお、レーダー画像データの作成手順は、“J. of Society of Naval Architects of Japan, vol.187, pp.85-92, 2000”に掲載されているTakaseらの論文に記載のものと同じものであるので、その詳細の説明は省略する。
【００２８】
そして、陰影比ＳＲの算出後、波高算出装置１０は、例えば、平均波周期が12.08秒であったときには、Ａ＝2.8968、Ｂ＝1.1632として、有義波高Ｈ1/3の算出し、ステップＳ１０５を終了する。
【００２９】
ステップＳ１０５の終了後、波高算出装置１０は、算出した有義波高Ｈ1/3を、ディスプレイに表示するとともに、ＨＤＤに記憶する（ステップＳ１０６）。波高表示装置１１が接続されている場合、波高算出装置１０は、このステップＳ１０６において、接続されている波高表示装置１１に、算出した有義波高Ｈ1/3を表示させるための処理も行なう。
【００３０】
その後、波高算出装置１０は、波高算出条件にて指定されている全ての波高の算出が完了しているか否かを判断する（ステップＳ１０７）。波高算出装置１０は、指定されている全ての波高の算出が完了していなかった場合（ステップＳ１０７；ＮＯ）、ステップＳ１０１に戻り、算出時刻がくるのを再び待機する。そして、波高算出装置１０は、指定されている全ての波高の算出が完了したとき（ステップＳ１０７；ＹＥＳ）に、図示した処理を終了する。
【００３１】
以下、波高算出装置１０の性能について説明する。
【００３２】
図４及び図５に、波高算出装置１０により算出された有義波高Ｈ1/3と他の波浪測定デバイスで測定された有義波高Ｈ1/3との関係を示す。なお、図４及び図５に示した算出結果は、図６、７に示してあるように、アンテナ角速度が２２ｒｐｍ、距離分解能が７．５ｍ、水平ビーム幅が０．８°といったレーダー測定条件下、波高算出装置１０が取り付けられた船舶を、図８に示してあるような各条件で運転したときに、得られたものである。
【００３３】
図４に示してあるように、波高算出装置１０により算出された有義波高Ｈ1/3は、波浪ブイにより得られた有義波高Ｈ1/3、及び、視覚観測により得られた有義波高Ｈ1/3と良く一致している。また、波高算出装置１０により算出された有義波高Ｈ1/3は、図５に示してあるように、波浪追算により得られた有義波高Ｈ1/3、及び、船舶の挙動解析により得られた有義波高Ｈ1/3とも良く一致している。
【００３４】
図９に、波高算出装置１０により算出された有義波高Ｈ1/3と波浪追算により得られた有義波高Ｈ1/3との関係を、従来の有義波高算出手順で算出された有義波高Ｈ1/3と波浪追算により得られた有義波高Ｈ1/3との関係と共に示す。
【００３５】
図から明らかなように、本波高算出装置１０で用いられている算出手順（ＳＲ法）により算出された有義波高Ｈ1/3は、従来の有義波高算出手順（ＳＮＲ法）で算出された有義波高Ｈ1/3よりも、波浪追算により得られた有義波高Ｈ1/3との間に強い相関があるものとなっている。
【００３６】
以上、説明したように、本実施形態に係る波高算出装置１０によれば、従来法（ＳＮＲ法）よりも、比較的に小さな波についての有義波高Ｈ1/3を精度良く算出できる。従って、波高算出装置１０を用いれば、船舶及び浮動構造の設計や船舶等の保守と船舶のウェザールーティング等により役立つ情報を得ることが出来ることになり、その結果として、操船の安全性及び効率を向上させることが出来ることになる。
【００３７】
【発明の効果】
本発明の波高算出装置、波高算出方法によれば、従来よりも、比較的に小さな波についての波高を精度良く算出できる。また、本発明のコンピュータ可読媒体によれば、コンピュータを本発明の波高算出装置として機能させることが出来る。従って、本発明によれば、名船舶及び浮動構造の設計や船舶等の保守と船舶のウェザールーティング等により役立つ情報を得ることが出来ることになり、その結果として、操船の安全性及び効率を向上させることが出来ることになる。
【図面の簡単な説明】
【図１】本発明の一実施形態に係る波高算出装置の使用形態の説明図である。
【図２】実施形態に係る波高算出装置が波高の算出を行なう領域の説明図である。
【図３】実施形態に係る波高算出装置が実行する波高算出処理の流れ図である。
【図４】実施形態に係る波高算出装置による波高算出結果を示した図である。
【図５】実施形態に係る波高算出装置による波高算出結果を示した図である。
【図６】図４及び図５に示した結果が得られた測定環境の説明図である。
【図７】図４及び図５に示した結果が得られた測定環境の説明図である。
【図８】図４及び図５に示した結果が得られた測定環境の説明図である。
【図９】実施形態に係る波高算出装置による波高算出結果と従来法による波高算出結果を対比した図である。
【符号の説明】
１０ 波高算出装置
１１ 波高表示装置
２０ ＣＤ−ＲＯＭ
２５ 波高算出プログラム
５１ インターフェースモジュール
５２ レーダーユニット
５３ レーダー指示器[0001]
BACKGROUND OF THE INVENTION
The present invention relates to a wave height calculation device and a wave height calculation method for calculating the height of a wave generated on a water surface using a marine radar, and a computer-readable medium for operating a computer as a wave height calculation device. Related to.
[0002]
[Prior art]
Statistical and / or real-time information on the wave height (hereinafter referred to as wave height) generated on the water surface (mainly the sea surface) is not only useful for the design of ships and floating structures, but also for maintenance of ships and ships. It is also useful for weather routing, and based on this information, the safety and efficiency of maneuvering can be improved.
[0003]
For this reason, various methods for measuring the wave height have been developed, and as one of them, a method for measuring (calculating) the wave height using a marine radar is known.
[0004]
In a conventional wave height measurement method using a ship radar, a signal-to-noise ratio (a value obtained by dividing a spectrum peak by a background noise: SNR) is calculated based on information obtained by the ship radar, and the wave height is calculated from the calculated SNR. Is estimated.
[0005]
Specifically, according to the equation (1), the significant wave height H 1/3 (the wave height having the average wave height for the waves that are between 1/3 of the whole wave counted from the higher one in the group of waves) Sa) is demanded.
[0006]
[Expression 1]
Figure 0004719373
Note that α and β are used for calibration in which their values are determined based on measurement results obtained by another method such as a method using a wave buoy so that a significant wave height can be obtained by the equation (1). It is a parameter.
[0007]
[Problems to be solved by the invention]
When the conventional wave height calculation procedure described above is used, the wave height can be accurately estimated when a relatively large wave is generated, but the wave height cannot be accurately estimated when a small wave is generated (measurement accuracy is low). There was a problem.
[0008]
Accordingly, an object of the present invention is to provide a wave height calculation device and wave height calculation method using a marine radar that can accurately calculate (estimate) the wave height of a smaller wave, and to operate a computer as such a wave height calculation device. A computer-readable medium is provided.
[0009]
[Means for Solving the Problems]
In order to solve the above-described problem, the wave height calculation device according to the present invention is based on a radar image signal and has a signal level on a radar image corresponding to a water surface area where the wave height is to be calculated. The shadow ratio calculation means for calculating the shadow ratio, which is a value obtained by dividing the area of the shaded area by the area of the water surface area, and the shadow ratio calculated by the shadow ratio calculation means are substituted into a predetermined arithmetic expression. A wave height calculating means for calculating the wave height, and an average period calculating means for calculating an average period of the waves generated in the water surface region where the wave height is to be calculated based on the radar image signal. Including a predetermined number of parameters defined as a function of the average period of waves occurring in the water surface region where the wave height is to be calculated, the wave height calculating means is the average calculated by the average period calculating means Obtains a value of the predetermined number of parameters based on the period, it has a configuration to calculate the height by substituting the shading ratio value obtained with the predetermined arithmetic expression the value of each parameter.
[0010]
Further, in the wave height calculation method of the present invention, on the radar image corresponding to the water surface area where the wave height is to be calculated based on the radar image signal, the signal level of the shaded portion is a level within a predetermined range. The shadow ratio calculation step for calculating the shadow ratio, which is a value obtained by dividing the area by the area of the water surface area, and the shadow ratio calculated by this shadow ratio calculation step are used to calculate the wave height generated in the water surface area where the wave height is to be calculated. A wave height calculation step for calculating a wave height by substituting it into a predetermined arithmetic expression including a predetermined number of parameters determined as a function of an average period, and is generated in a water surface region where the wave height is to be calculated based on a radar image signal. The average period of the wave is calculated, the value of the predetermined number of parameters is obtained based on the calculated average period, and the predetermined arithmetic expression with the calculated value as the value of each parameter is calculated in advance. By a wave height calculating a height by substituting the shade ratio, wave height is calculated.
[0011]
The computer-readable medium of the present invention is a computer readable medium that connects a radar device to a radar image corresponding to a water surface area whose wave height is to be calculated based on a radar image signal. The shadow ratio calculation means for calculating the shadow ratio, which is a value obtained by dividing the area of the shadow portion at the level of the water surface area by the area, and the shadow ratio calculated by the shadow ratio calculation means to calculate the wave height. A wave height calculating means for calculating a wave height by substituting into a predetermined arithmetic expression including a predetermined number of parameters defined as a function of an average period of waves generated in a power surface area , based on a radar image signal, The average period of waves generated in the water surface area where the wave height is to be calculated is calculated, the values of the predetermined number of parameters are determined based on the calculated average period, and The height calculation means for calculating the pulse height by substituting the shading ratio to the predetermined arithmetic expression is a value of over data, comprising recording a program for operating the provided device.
[0012]
That is, as a result of diligent research to solve the above-mentioned problems, the inventors have found that the wave height of a small wave can be accurately calculated by using the shadow ratio as a wave height index. Therefore, according to the wave height calculation apparatus and wave height calculation method of the present invention in which the wave height is calculated based on this principle, the wave height can be calculated more accurately than in the past. Further, according to the computer readable medium of the present invention, a computer connectable with the radar device can function as the wave height calculating device of the present invention.
[0015]
Further, when realizing the wave height calculation apparatus of the present invention, two parameters defined as a function of the shadow ratio SR and the average period of the waves generated in the water surface region where the wave height is to be calculated are given as predetermined arithmetic expressions. From A and B, the equation “H = A · SR B ” for calculating the wave height H can be adopted.
[0016]
DETAILED DESCRIPTION OF THE INVENTION
Hereinafter, embodiments of the present invention will be described in detail with reference to the drawings.
[0017]
FIG. 1 shows a usage pattern of a wave height calculation apparatus 10 according to an embodiment of the present invention.
[0018]
As illustrated, the wave height calculation apparatus 10 according to the present embodiment is used in a state of being connected to the radar unit 52 via the interface module 51. The radar unit 52 is a system including a transmitter, a radar antenna, a receiver, and the like. The interface module 51 is a module that supplies a signal (radar image signal) output from the receiver of the radar unit 52 to the radar indicator 53 and the wave height calculation device 10.
[0019]
The radar indicator 53 is a device that displays a radar image corresponding to the input radar image signal, and can be directly connected to the radar unit 52. In other words, a portion composed of the radar unit 52 and the radar indicator 53 is generally called a marine radar, and the interface module 51 receives a radar image from the marine radar (radar unit 52). It is a module for extracting signals.
[0020]
The wave height calculation device 10 is a computer in which the wave height calculation program 25 is installed from the CD-ROM 20 or the like, and includes an interface for acquiring a radar image signal from the interface module 51, a GPS signal, a ship speed signal, and a course signal. And an interface for obtaining from a steering device (not shown). Further, the wave height calculation device 10 includes an interface for connecting to the wave height display device 11. The wave height display device 11 is a device including a liquid crystal display on which the wave height and the like are displayed, and the wave height calculation device 10 functions even if the wave height display device 11 is not connected. The wave height display device 11 connected to the wave height calculation device 10 is usually installed in the vicinity of the radar indicator 53.
[0021]
Hereinafter, the operation of the wave height calculation apparatus 10 will be described.
[0022]
When the power is turned on, the wave height calculation device 10 displays a screen including a region where the radar image is displayed on the display, and waits for the wave height calculation condition to be set by the operator. Here, the wave height calculation condition is information for the operator to instruct the wave height calculation device 10 at which point of time in which region the wave height is to be calculated, and the wave height calculation condition is to calculate the wave height a plurality of times. Can be specified (multiple sets of time and region for calculating the wave height can be specified in advance).
[0023]
Note that, in the present wave height calculation apparatus 10, an area for which the wave height is to be calculated (hereinafter referred to as an observation target area), as shown in FIG. It is a square area that passes above (the center of the radar image), and when setting the wave height calculation conditions, the operator performs work including the designation of the observation target area.
[0024]
FIG. 3 shows a flowchart of the wave height calculation process executed by the wave height calculation device 10 when a predetermined instruction (instruction to start the wave height calculation process) is input after the setting of the wave height calculation conditions is completed.
[0025]
As shown in the figure, at the start of the wave height calculation process, the wave height calculation device 10 stands by until the calculation time designated as the time at which the wave height is to be calculated under the set wave height calculation conditions (step S101). Then, the wave height calculation device 10 ends step S101 when the calculation time comes, and collects the radar image signal output from the radar unit 52 in an amount specified by the wave height calculation condition via the interface module 51. Then, it is stored in the built-in HDD (step S102). In step S102, the wave height calculation apparatus 10 collects radar image signals for 32 (= 2 5 ) radar images, for example. After that, the wave height calculation device 10 performs a three-dimensional FFT (Fast Fourier Transform) or the like on the signal related to the observation target area defined by the wave height calculation condition included in the collected radar image signal, thereby generating noise ( A wave number / frequency spectrum from which information not related to waves is removed is generated (step S103).
[0026]
Next, the wave height calculating apparatus 10 obtains an average wave period from the generated wave number / frequency spectrum (step S104). Further, the wave height calculation device 10 generates radar image data from which noise has been removed by performing inverse FFT of the wave number and frequency spectrum, calculates the shadow ratio SR based on the radar image data, and then calculates the average wave period. The significant wave height H 1/3 is calculated by substituting the calculated SR into the relational expression H 1/3 = A · SR B using the parameters A and B prepared as a function of (Step S105). ).
[0027]
Specifically, in this step S105, the wave height calculating apparatus 10 calculates and calculates the area of a portion (hereinafter referred to as a shaded portion) where the intensity of the radar image has a substantially constant value at a low level. A value obtained by dividing the shaded area by the total area of the observation target area is calculated as the shadow ratio SR. In addition, as a judgment condition of the shadow portion, various radar image data are created by simulation assuming various values (including significant wave height H 1/3 ) relating to waves, and the shadow calculated from the created radar image data. The ratio SR is determined so as to coincide with the assumed significant wave height H 1/3 . The radar image data creation procedure is the same as that described in the paper by Takase et al. Published in “J. of Society of Naval Architects of Japan, vol.187, pp.85-92, 2000”. Since there is, the detailed description is abbreviate | omitted.
[0028]
After calculating the shadow ratio SR, for example, when the average wave period is 12.08 seconds, the wave height calculating apparatus 10 calculates the significant wave height H 1/3 as A = 2.8968, B = 1.1632, and step S105. Exit.
[0029]
After step S105, the wave height calculating apparatus 10 displays the calculated significant wave height H 1/3 on the display and stores it in the HDD (step S106). When the wave height display device 11 is connected, the wave height calculation device 10 also performs processing for displaying the calculated significant wave height H 1/3 on the connected wave height display device 11 in step S106.
[0030]
Thereafter, the wave height calculation device 10 determines whether or not calculation of all the wave heights specified in the wave height calculation conditions has been completed (step S107). When the calculation of all the specified wave heights has not been completed (step S107; NO), the wave height calculation apparatus 10 returns to step S101 and waits again for the calculation time to come. Then, the wave height calculation apparatus 10 ends the illustrated process when the calculation of all the specified wave heights is completed (step S107; YES).
[0031]
Hereinafter, the performance of the wave height calculation apparatus 10 will be described.
[0032]
4 and 5 show the relationship between the significant wave height H 1/3 calculated by the wave height calculation device 10 and the significant wave height H 1/3 measured by another wave measuring device. The calculation results shown in FIGS. 4 and 5 are as shown in FIGS. 6 and 7, under the radar measurement conditions where the antenna angular velocity is 22 rpm, the distance resolution is 7.5 m, and the horizontal beam width is 0.8 °. This is obtained when the ship to which the wave height calculating device 10 is attached is operated under the conditions as shown in FIG.
[0033]
As is shown in FIG. 4, the significant wave height H 1/3 calculated by the pulse height calculating unit 10, the significant wave height H 1/3 obtained by wave buoy, and the significant obtained by visual observation It is in good agreement with the wave height H 1/3 . Moreover, the significant wave height H 1/3 calculated by pulse height calculation device 10, as are shown in Figure 5, the significant wave height obtained by wave additionally calculated H 1/3, and by Analysis of Vessel It also agrees well with the significant wave height H 1/3 obtained.
[0034]
In FIG. 9, the relationship between the significant wave height H 1/3 calculated by the wave height calculation device 10 and the significant wave height H 1/3 obtained by wave addition is calculated by the conventional significant wave height calculation procedure. together with the relationship between the significant wave height H 1/3 obtained by the significant wave height H 1/3 and waves add calculated.
[0035]
As is apparent from the figure, the significant wave height H 1/3 calculated by the calculation procedure (SR method) used in the present wave height calculation device 10 is calculated by the conventional significant wave height calculation procedure (SNR method). than the significant wave height H 1/3 was, thereby providing a strong correlation between the significant wave height H 1/3 obtained by waves add calculation.
[0036]
As described above, according to the wave height calculation apparatus 10 according to the present embodiment, the significant wave height H 1/3 for a relatively small wave can be calculated with higher accuracy than in the conventional method (SNR method). Therefore, if the wave height calculation device 10 is used, useful information can be obtained by designing the ship and the floating structure, maintaining the ship, etc., weather routing of the ship, etc. As a result, the safety and efficiency of the ship maneuvering can be improved. It can be improved.
[0037]
【The invention's effect】
According to the wave height calculation apparatus and the wave height calculation method of the present invention, it is possible to calculate the wave height of a relatively small wave with higher accuracy than in the past. In addition, according to the computer-readable medium of the present invention, the computer can function as the wave height calculation apparatus of the present invention. Therefore, according to the present invention, useful information can be obtained through design of famous ships and floating structures, maintenance of ships, etc., weather routing of ships, etc., and as a result, the safety and efficiency of ship operation are improved. Will be able to.
[Brief description of the drawings]
FIG. 1 is an explanatory diagram of a usage pattern of a wave height calculation apparatus according to an embodiment of the present invention.
FIG. 2 is an explanatory diagram of a region where a wave height calculation apparatus according to an embodiment calculates a wave height.
FIG. 3 is a flowchart of a wave height calculation process executed by the wave height calculation apparatus according to the embodiment.
FIG. 4 is a diagram showing a wave height calculation result by the wave height calculating apparatus according to the embodiment.
FIG. 5 is a diagram showing a wave height calculation result by the wave height calculating apparatus according to the embodiment.
6 is an explanatory diagram of a measurement environment in which the results shown in FIG. 4 and FIG. 5 are obtained.
7 is an explanatory diagram of a measurement environment in which the results shown in FIG. 4 and FIG. 5 are obtained.
FIG. 8 is an explanatory diagram of a measurement environment in which the results shown in FIGS. 4 and 5 are obtained.
FIG. 9 is a diagram comparing a wave height calculation result by a wave height calculation apparatus according to an embodiment with a wave height calculation result by a conventional method.
[Explanation of symbols]
10 Wave height calculation device 11 Wave height display device 20 CD-ROM
25 Wave height calculation program 51 Interface module 52 Radar unit 53 Radar indicator

Claims (4)

レーダー画像信号に基づき、波高を算出すべき水面領域に相当するレーダー画像上の、信号レベルが予め定められている範囲内のレベルとなっている陰影部の面積を当該水面領域の面積で割った値である陰影比を算出する陰影比算出手段と、
この陰影比算出手段によって算出された陰影比を、所定の演算式に代入することにより波高を算出する波高算出手段と、
レーダー画像信号に基づき、波高を算出すべき水面領域に生じている波の平均周期を算出する平均周期算出手段とを、備え、
前記所定の演算式は、波高を算出すべき水面領域に生じている波の平均周期の関数として定められている所定数のパラメータを含み、
前記波高算出手段は、前記平均周期算出手段によって算出された平均周期に基づき前記所定数のパラメータの値を求め、求めた値を各パラメータの値とした前記所定の演算式に前記陰影比を代入することにより波高を算出する
ことを特徴とする波高算出装置。Based on the radar image signal, the area of the shaded area where the signal level is within a predetermined range on the radar image corresponding to the water surface area whose wave height is to be calculated is divided by the area of the water surface area. A shadow ratio calculating means for calculating a shadow ratio which is a value;
A wave height calculating means for calculating the wave height by substituting the shadow ratio calculated by the shadow ratio calculating means into a predetermined arithmetic expression ;
An average period calculating means for calculating an average period of waves generated in the water surface region where the wave height is to be calculated based on the radar image signal;
The predetermined arithmetic expression includes a predetermined number of parameters defined as a function of an average period of waves generated in the water surface region where the wave height is to be calculated,
The wave height calculating means obtains the value of the predetermined number of parameters based on the average period calculated by the average period calculating means, and substitutes the shadow ratio into the predetermined arithmetic expression using the obtained value as the value of each parameter. To calculate the wave height
Wave height calculating device, characterized in that. 前記所定の演算式が、
陰影比ＳＲと、波高を算出すべき水面領域に生じている波の平均周期の関数として定められている２つのパラメータＡ及びＢとから、波高Ｈを算出する「Ｈ＝Ａ・ＳＲB」という式である
ことを特徴とする請求項１に記載の波高算出装置。The predetermined arithmetic expression is
“H = A · SR B ” for calculating the wave height H from the shadow ratio SR and the two parameters A and B defined as a function of the average period of the waves generated in the water surface region where the wave height is to be calculated. The wave height calculation apparatus according to claim 1 , wherein the wave height calculation apparatus is an expression. レーダー画像信号に基づき、波高を算出すべき水面領域に相当するレーダー画像上の、信号レベルが予め定められている範囲内のレベルとなっている陰影部の面積を当該水面領域の面積で割った値である陰影比を算出する陰影比算出ステップと、
この陰影比算出ステップによって算出された陰影比を、波高を算出すべき水面領域に生じている波の平均周期の関数として定められている所定数のパラメータを含む所定の演算式に代入することにより波高を算出する波高算出ステップであって、レーダー画像信号に基づき、波高を算出すべき水面領域に生じている波の平均周期を算出し、算出した平均周期に基づき前記所定数のパラメータの値を求め、求めた値を各パラメータの値とした前記所定の演算式に前記陰影比を代入することにより波高を算出する波高算出ステップ
とを、含むことを特徴とする波高算出方法。Based on the radar image signal, the area of the shaded area where the signal level is within a predetermined range on the radar image corresponding to the water surface area whose wave height is to be calculated is divided by the area of the water surface area. A shadow ratio calculating step for calculating a shadow ratio as a value;
By substituting the shadow ratio calculated by this shadow ratio calculation step into a predetermined arithmetic expression including a predetermined number of parameters defined as a function of the average period of waves generated in the water surface region where the wave height is to be calculated. A wave height calculating step for calculating a wave height, calculating an average period of waves generated in a water surface region where the wave height is to be calculated based on a radar image signal, and calculating the value of the predetermined number of parameters based on the calculated average period. A wave height calculating step of calculating a wave height by substituting the shading ratio into the predetermined arithmetic expression using the calculated values as values of the respective parameters . レーダー装置と接続可能なコンピュータを、
レーダー画像信号に基づき、波高を算出すべき水面領域に相当するレーダー画像上の、信号レベルが予め定められている範囲内のレベルとなっている陰影部の面積を当該水面領域の面積で割った値である陰影比を算出する陰影比算出手段と、
この陰影比算出手段によって算出された陰影比を、波高を算出すべき水面領域に生じている波の平均周期の関数として定められている所定数のパラメータを含む所定の演算式に代入することにより波高を算出する波高算出手段であって、レーダー画像信号に基づき、波高を算出すべき水面領域に生じている波の平均周期を算出し、算出した平均周期に基づき前記所定数のパラメータの値を求め、求めた値を各パラメータの値とした前記所定の演算式に前記陰影比を代入することにより波高を算出する波高算出手段
とを、備える装置として動作させるためのプログラムを記録したことを特徴とするコンピュータ可読媒体。A computer that can be connected to a radar device
Based on the radar image signal, the area of the shaded area where the signal level is within a predetermined range on the radar image corresponding to the water surface area whose wave height is to be calculated is divided by the area of the water surface area. A shadow ratio calculating means for calculating a shadow ratio which is a value;
By substituting the shadow ratio calculated by the shadow ratio calculating means into a predetermined arithmetic expression including a predetermined number of parameters determined as a function of the average period of the waves generated in the water surface region where the wave height is to be calculated. Wave height calculating means for calculating a wave height, calculating an average period of waves generated in a water surface region where the wave height is to be calculated based on a radar image signal, and calculating the value of the predetermined number of parameters based on the calculated average period. A program for operating as a device comprising a wave height calculating means for calculating a wave height by substituting the shading ratio into the predetermined arithmetic expression using the calculated value as a value of each parameter. A computer-readable medium characterized by that.
JP2001124160A 2001-04-23 2001-04-23 Wave height calculation device, wave height calculation method, and computer-readable medium Expired - Lifetime JP4719373B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2001124160A JP4719373B2 (en) 2001-04-23 2001-04-23 Wave height calculation device, wave height calculation method, and computer-readable medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2001124160A JP4719373B2 (en) 2001-04-23 2001-04-23 Wave height calculation device, wave height calculation method, and computer-readable medium

Publications (2)

Publication Number Publication Date
JP2002318115A JP2002318115A (en) 2002-10-31
JP4719373B2 true JP4719373B2 (en) 2011-07-06

Family

ID=18973590

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2001124160A Expired - Lifetime JP4719373B2 (en) 2001-04-23 2001-04-23 Wave height calculation device, wave height calculation method, and computer-readable medium

Country Status (1)

Country Link
JP (1) JP4719373B2 (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103578106A (en) * 2013-11-06 2014-02-12 北京航空航天大学 Out-of-order water surface compiling method based on truthful data
KR101790482B1 (en) * 2015-12-30 2017-10-25 경남대학교 산학협력단 Wave measurement system and method using marine rader
KR101790481B1 (en) 2015-12-30 2017-10-25 경남대학교 산학협력단 Sea surface wind measurement system and method using marine rader

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3992101B2 (en) * 2003-05-16 2007-10-17 よこはまティーエルオー株式会社 Individual wave prediction / warning system
KR100795497B1 (en) * 2006-07-21 2008-01-17 삼성중공업 주식회사 Wave measure method and system using radar
JP6154218B2 (en) 2013-07-03 2017-06-28 古野電気株式会社 Echo signal processing device, wave radar device, echo signal processing method, and echo signal processing program
JP2022091049A (en) 2020-12-08 2022-06-20 ヤマハ発動機株式会社 Boat
CN113030894B (en) * 2021-03-02 2022-06-28 南京信息工程大学 Method for extracting sea wave parameters by using rapidly scanned coherent radar image

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6022681A (en) * 1983-07-19 1985-02-05 Japan Radio Co Ltd Method and device for wave height radar observation
JPS6093317A (en) * 1983-10-27 1985-05-25 Japan Radio Co Ltd Radar crest measuring device
JP2000310678A (en) * 1999-04-28 2000-11-07 Oki Electric Ind Co Ltd Wave observation system with radar
JP2001059869A (en) * 1999-08-24 2001-03-06 Oki Electric Ind Co Ltd Physical quantity observing system by proportional scale of radar reception signal

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6022681A (en) * 1983-07-19 1985-02-05 Japan Radio Co Ltd Method and device for wave height radar observation
JPS6093317A (en) * 1983-10-27 1985-05-25 Japan Radio Co Ltd Radar crest measuring device
JP2000310678A (en) * 1999-04-28 2000-11-07 Oki Electric Ind Co Ltd Wave observation system with radar
JP2001059869A (en) * 1999-08-24 2001-03-06 Oki Electric Ind Co Ltd Physical quantity observing system by proportional scale of radar reception signal

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103578106A (en) * 2013-11-06 2014-02-12 北京航空航天大学 Out-of-order water surface compiling method based on truthful data
CN103578106B (en) * 2013-11-06 2016-08-17 北京航空航天大学 A kind of out of order water surface edit methods based on truthful data
KR101790482B1 (en) * 2015-12-30 2017-10-25 경남대학교 산학협력단 Wave measurement system and method using marine rader
KR101790481B1 (en) 2015-12-30 2017-10-25 경남대학교 산학협력단 Sea surface wind measurement system and method using marine rader

Also Published As

Publication number Publication date
JP2002318115A (en) 2002-10-31

Similar Documents

Publication Publication Date Title
US7532747B2 (en) Method and apparatus for ultrasonic imaging in mmode
US6775617B2 (en) Method for determining hydrographic parameters which describe a sea swell field in situ using a radar device
EP1726262B1 (en) Apparatus and method for processing an ultrasound spectrum image
JP2008128900A (en) Echo image display apparatus
JP2008136860A (en) Ultrasonic diagnostic apparatus and image processing program for it
JP4719373B2 (en) Wave height calculation device, wave height calculation method, and computer-readable medium
US8300909B2 (en) Ultrasonographic device and ultrasonographic method
US10359508B2 (en) System and method for calibration of echo sounding systems and improved seafloor imaging using such systems
CN105241428A (en) Water depth retrieval method with hyper-spectrum
US5485844A (en) Doppler-type ultrasonic diagnostic apparatus
WO1994023652A1 (en) Tissue characterisation using intravascular echoscopy
JP2006288974A (en) Ultrasonic diagnostic equipment
CN113640808B (en) Shallow water submarine cable buried depth detection method and device
JP2003275210A (en) Method for describing spectrum strain and apparatus therefor
CN112882037B (en) Side-scan sonar sea bottom line detection method and device
JP2007040734A (en) Transmission mode recommend system, sonar system, transmission mode recommend method, and its program and storage medium
JP3888671B2 (en) Wave height calculation device, wave height calculation method, recording medium, and ship
JP2003004848A (en) System for conjecturing and visualizing sonar searching range, method using the same and program for the system
US7371219B2 (en) Ultrasound diagnosis apparatus operable in doppler mode
JP4750508B2 (en) Nautical Chart Information Processing Method, Nautical Chart Information System, Nautical Chart Information Program, and Recording Medium
CN116519056A (en) River bottom siltation dynamic monitoring method based on unmanned ship
JP4827330B2 (en) Radar wave measuring method and apparatus
JP7219603B2 (en) Sea level measurement device and sea level measurement program
CN115932884B (en) Wave direction spectrum measurement method and system based on three-dimensional laser radar
JP3589186B2 (en) Onboard marine forecasting device

Legal Events

Date Code Title Description
A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A712

Effective date: 20040419

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20040419

A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A711

Effective date: 20050921

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20050921

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20051116

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20071116

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20100610

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20100615

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20100805

A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A712

Effective date: 20101028

RD03 Notification of appointment of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7423

Effective date: 20101105

RD04 Notification of resignation of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7424

Effective date: 20101105

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20101028

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20101108

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20101214

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20110210

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20110308

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20110404

R150 Certificate of patent or registration of utility model

Ref document number: 4719373

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20140408

Year of fee payment: 3

S531 Written request for registration of change of domicile

Free format text: JAPANESE INTERMEDIATE CODE: R313531

S533 Written request for registration of change of name

Free format text: JAPANESE INTERMEDIATE CODE: R313533

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20140408

Year of fee payment: 3

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313117

R360 Written notification for declining of transfer of rights

Free format text: JAPANESE INTERMEDIATE CODE: R360

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R370 Written measure of declining of transfer procedure

Free format text: JAPANESE INTERMEDIATE CODE: R370

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313117

S533 Written request for registration of change of name

Free format text: JAPANESE INTERMEDIATE CODE: R313533

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313117

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313115

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

S531 Written request for registration of change of domicile

Free format text: JAPANESE INTERMEDIATE CODE: R313531

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313115

S533 Written request for registration of change of name

Free format text: JAPANESE INTERMEDIATE CODE: R313533

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250