スキャニングソナー装置
【課題】外乱による反射波の強度の揺らぎによる影響を低減し標的の俯角方向に対する移動の向き検出の精度を改善したスキャニングソナー装置を提供する。
【解決手段】スキャニングソナー装置１は、標的像抽出部３１と特徴量抽出部３２と移動方向判定部３３とを具備し、標的像抽出部３１が、音波が標的で反射した反射波を受波して得られた反射波の強さを示す複数のサンプリング値から標的の形状を示す標的像を抽出し、特徴量抽出部３２が、標的像抽出部が抽出した標的像に含まれるサンプリング値の空間的平均値から標的の特徴量を算出する。標的像抽出部３１が、水面に対して複数の異なる俯角に送波された音波それぞれに対して複数の異なる俯角ごとの標的像を抽出し、特徴量抽出部３２が、当該複数の標的像それぞれの特徴量を算出し、移動方向判定部３３が、算出された複数の特徴量を比較することにより標的の俯角方向に対する移動の向きを判定する。
【解決手段】スキャニングソナー装置１は、標的像抽出部３１と特徴量抽出部３２と移動方向判定部３３とを具備し、標的像抽出部３１が、音波が標的で反射した反射波を受波して得られた反射波の強さを示す複数のサンプリング値から標的の形状を示す標的像を抽出し、特徴量抽出部３２が、標的像抽出部が抽出した標的像に含まれるサンプリング値の空間的平均値から標的の特徴量を算出する。標的像抽出部３１が、水面に対して複数の異なる俯角に送波された音波それぞれに対して複数の異なる俯角ごとの標的像を抽出し、特徴量抽出部３２が、当該複数の標的像それぞれの特徴量を算出し、移動方向判定部３３が、算出された複数の特徴量を比較することにより標的の俯角方向に対する移動の向きを判定する。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、スキャニングソナー装置、主に設定した標的を自動追尾するスキャニングソナー装置に関する。
【背景技術】
【０００２】
従来、漁労などに用いられるスキャニングソナー装置では、例えば、図２に示すように船体などに設けられ水中に位置する送受波機５０が、一定の周波数、例えば、２０ｋＨｚから３００ｋＨｚまでの間の周波数の音波を、ユーザによって設定される海面からθの角度（俯角という）に送波して、受波したエコー（反射波）の強度情報（サンプリング値）を用いて標的である魚群などの探査を行い、その位置検出をしている。また、スキャニングソナー装置は、反射波の強度を予め定めた閾値を超える領域を図３（ａ）又は図３（ｂ）に示すように表示装置に画像として表示する。ここで、図３（ａ）は、縦軸方向及び横軸方向共に中心から周囲へ向かういずれの方向に対しても俯角θを有する円錐体の傘状の表面における送受波機５０からの距離を示している。図３（ｂ）は、縦軸方向が俯角θを有する円錐体の傘状の表面における送受波機５０からの距離を示し、横軸方向が船首又は船尾を基準とし海面に対する垂直な軸ｚ（図２）の回転方向の角度φを示している。
【０００３】
ここで、送受波機５０は、図２に示すように狭い範囲からの反射波を受波する指向性を有するペンシル形状の受波ビームＰを、海面に対する垂直方向（ｚ軸方向）に対して時計回り又は反時計回りに回転させ全周３６０度からの反射波を受信して円錐体の傘状の表面に対してスキャンを行う。
このようなスキャニングソナー装置において、標的の移動方向を検出する際には、異なる俯角に対して複数回の送波及び受波を行い、受波した反射波の強度分布や、標的像の面積、及び、寸法（図３（ｂ）におけるΔφ及びΔＲ）を特徴量として抽出して標的の俯角方向に対する移動の向きの推定を行うものもある（特許文献１）。
【０００４】
ところで、反射波を受波する送受波機の移動及び波による揺れ、標的である魚などの移動、及び、潮流や温度により音波の伝搬路の状態が時間的変化することによって生じる外乱により、受波する反射波の強度は、主にスペックルと同様の原理により空間的及び時間的に揺らぐという現象がある。このため、空間的及び時間的に揺らいだ反射波の強度から標的像を抽出すると、抽出した標的像の形状が、実際の標的像に比べ歪むと共に、抽出した標的像の分布が、実際の標的の分布に比べ揺らいで安定しない。上述のようなスキャニングソナー装置においては、複数の異なる俯角で音波を送波し、その反射波を受波して標的を立体的に捉えようとする場合、上述の揺らぎの影響を低減させるために、複数の異なる俯角それぞれにおいて複数回の探査を行って得られた標的像の平均値を算出して、時間的な揺らぎの影響を緩和した標的像の抽出と、特徴量の算出とが行われている。
【特許文献１】特開２００８−２０３２２７号公報
【発明の開示】
【発明が解決しようとする課題】
【０００５】
しかしながら、平均値を算出するために複数回探査を行おうとすると、特徴量算出に要する時間が長くなり、その間に標的が移動するなどして、実際の標的の俯角方向に対する移動の向きと、得られた特徴量から得られる標的の俯角方向に対する移動の向きとの乖離が大きくなり、標的の俯角方向に対する移動の向き検出の精度が悪くなるという問題がある。
【０００６】
本発明は、上記問題を解決すべくなされたもので、その目的は、外乱による反射波の強度の揺らぎによる影響を低減し、標的の俯角方向に対する移動の向き検出の精度を改善したスキャニングソナー装置を提供することにある。
【課題を解決するための手段】
【０００７】
上記問題を解決するために、本発明は、音波を標的に送波し、その反射波を受波して得られる反射波の強度に基づく複数のサンプリング値から前記標的の形状を示す標的像を抽出する標的像抽出部と、前記標的像抽出部が抽出した前記標的像に含まれる前記サンプリング値から前記標的像の特徴量を算出する特徴量算出部と、前記標的像抽出部が算出した複数の俯角ごとの前記標的像から前記特徴量抽出部によって算出された当該複数の俯角それぞれに対応した前記特徴量を比較して前記標的の俯角方向に対する移動の向きを判定する移動方向判定部と、を具備し、前記特徴量は、前記標的像に含まれる前記サンプリング値の平均値と前記標的像の面積とを乗じた数値であることを特徴とするスキャニングソナー装置である。
【０００８】
また、本発明は、音波を標的に送波し、その反射波を受波して得られた反射波の強度に基づく複数のサンプリング値から前記標的の形状を示す標的像を抽出する標的像抽出部と、前記標的像抽出部が抽出した前記標的像に含まれる前記サンプリング値から前記標的像の特徴量を算出する特徴量算出部と、前記標的像抽出部が算出した複数の俯角ごとの前記標的像から前記特徴量抽出部によって算出された当該複数の俯角それぞれに対応した前記特徴量を比較して、前記標的の俯角方向に対する移動の向きを判定する移動方向判定部と、を具備し、前記特徴量は、前記標的像を示す反射波の前記サンプリング値を総和した数値であることを特徴とするスキャニングソナー装置である。
【０００９】
また、本発明は、上記に記載の発明において、前記標的像の面積が予め定めた数値より小さいとき、前記特徴量算出部は、前記標的像抽出部が抽出した前記複数の異なる俯角ごとの複数の前記標的像から前記複数の異なる俯角ごとに複数の前記特徴量を算出し、算出した前記複数の特徴量から前記複数の異なる俯角ごとに前記特徴量の平均値を算出して、前記標的の俯角方向に対する移動の向きを判定することを特徴とする。
【発明の効果】
【００１０】
この発明によれば、スキャニングソナー装置は、水面に対して複数の異なる俯角で送波された音波が標的で反射した反射波（エコー）に基づいて、複数の異なる俯角ごとに標的の形状を示す標的像を示す反射波のサンプリング値を用いて特徴量を算出し、算出した特徴量を俯角間で比較することで標的の俯角方向に対する移動の向きを推定する。このとき、特徴量算出部は、空間平均として、標的像を示す反射波のサンプリング値の平均値と標的像の面積とを乗じた値、あるいは、標的像を示す反射波のサンプリング値を総和した値を特徴量として算出する。外乱による反射波の強度の揺らぎは、時間平均を用いずにエルゴード性を適用し空間平均を用いることで低減させて、複数回の送波及び受波を繰り返して時間平均を算出する代わりに空間平均を算出することで移動方向判定時間の短縮を行うと共に、反射波の強度の揺らぎを低減させることができる。これにより、スキャニングソナー装置は、探査回数を増やすことなく反射波の強度揺らぎを低減させ、標的の俯角方向に対する移動の向き検出の精度を改善することが可能となる。
【発明を実施するための最良の形態】
【００１１】
以下、本発明の実施形態によるスキャニングソナー装置を図面を参照して説明する。
【００１２】
＜第１実施形態＞
図１は、第１実施形態におけるスキャニングソナー装置１の構成を示す概略ブロック図である。以下、スキャニングソナー装置１は、魚群を標的として漁労に用いられる場合を例として説明する。
スキャニングソナー装置１は、入力装置１０、操作制御部２０、標的追尾制御部３０、送受波機制御部４０、送受波機５０、表示情報変換部６０、及び、表示装置７０を具備している。
入力装置１０は、ユーザの操作によりスキャニングソナー装置１への入力が行える、例えば、キーボード、マウス、トラックボールなどである。操作制御部２０は、入力装置１０から入力された情報の種類に応じて標的追尾制御部３０、送受波機制御部４０、表示情報変換部６０に出力する。ここで、入力装置１０から入力される情報としては、スキャニングソナー装置１に対して標的の追尾の開始を示す動作開始情報、音波を送波する俯角を示す俯角情報、標的をスキャンする範囲（レンジ）情報、標的（魚群）の設定情報、表示装置７０に表示するスキャン範囲を示す情報などである。
【００１３】
標的追尾制御部３０は、標的像抽出部３１、特徴量抽出部３２、移動方向判定部３３、及び俯角制御部３４を含んで構成され、操作制御部２０から動作開始情報が入力されると動作を開始し、操作制御部２０から入力される標的の設定情報及び俯角情報を初期値として標的の追尾処理を開始する。また、標的追尾制御部３０は、標的までの距離と標的の位置に応じた俯角とを含むデータを送受波機制御部４０に出力し、送受波機制御部４０から反射波（エコー）を受信することで得られた複数のサンプリング値が入力される。
ここで、サンプリング値は、図３（ｂ）に示される海面の垂直方向に対して時計回り又は反時計回り方向の角度φと距離Ｒとを用いて表されるＲ−φ座標系の各座標における受波した反射波の強さに基づく値（強度情報）である。なお、角度φは、船首や船尾など船体の一部を基準とし０°として定められる。また、角度φに対して０°から３６０°までの一周分の全探査距離に渡るサンプリング値を１フレームという。
【００１４】
標的像抽出部３１は、送受波機制御部４０から入力される１フレーム分のサンプリング値に対して、標的の形状を示す標的像の抽出を行う。ここで、標的像抽出部３１は、操作制御部２０を介して入力されるユーザの指定した閾値を基準とし、当該閾値を超えるサンプリング値により構成される領域を標的像として抽出を行う。例えば、標的像は、図３（ｂ）に図示される標的Ｔの領域に含まれるＲ−φ座標系の情報として出力される。また、ユーザの指定する閾値とは、漁労の対象となる魚種や観測環境ごとに定められる値であり、ユーザの経験などに基づいて定められる。また、標的像抽出部３１は、送受波機５０から抽出した標的像までの距離を送受波機制御部４０に出力する。
特徴量抽出部３２は、標的像抽出部３１が抽出した標的像の領域に含まれるサンプリング値から標的の移動方向判断に用いる特徴量Ｆを算出する。特徴量Ｆは、次式（１）により求める。
【００１５】
【数１】
【００１６】
ここで、式（１）における、Ａｐは、標的像の面積を表し、ｎｐは、標的像に含まれるサンプル数を表し、ｘ（ｐｉ）は、ｉ番目のサンプリング値を表す。すなわち、標的像に含まれる全サンプリング値の平均と標的像の面積とを乗じた数値が特徴量Ｆである。なお、標的像の面積は、Ｒ−φ座標から算出される。
【００１７】
移動方向判定部３３は、標的像抽出部３１から入力される複数の異なる俯角ごとの特徴量Ｆを比較して、魚群の俯角方向における移動の向きを判断する。判断方法としては、俯角θ−Δθ、θ、θ＋Δθ（Δθ＞０）それぞれの特徴量Ｆを比較し、俯角θ−Δθの特徴量Ｆが最も大きいとき、魚群は俯角が小さくなる方向に移動していると判断し、俯角θの特徴量Ｆが最も大きいとき、魚群は移動していないと判断し、俯角θ＋Δθの特徴量Ｆが最も大きいとき、魚群は俯角が大きくなる方向に移動していると判断する。すなわち、移動方向判定部３３は、特徴量Ｆの大小比較により魚群の俯角方向に対する移動の向き判断する。魚群の魚の量は、平均的な反射波の強度と相関が高いのでこの値が大きい領域を追うことで、魚群における魚が密集している領域を追尾することができ、判定の誤りによる魚群を見失うことを防ぐことが可能になる。
なお、移動方向判定部３３が、３つの異なる俯角ごとの特徴量Ｆを用いて魚群の移動の向きを判断する例を示したが、２つ以上の異なる俯角ごとの特徴量Ｆを用いて魚群の移動の向きを判断するようにしてもよい。
【００１８】
俯角制御部３４は、移動方向判定部３３により判断された魚群の移動の向きに応じて、送受波機５０が送波及び受波する俯角情報を更新して、送受波機制御部４０に出力する。
例えば、移動方向判定部３３が上述のように判断を行い、俯角が小さくなる方向に移動しているとき、現在の俯角情報θに対して俯角θ−Δθを次の俯角情報として更新する。
【００１９】
送受波機制御部４０は、操作制御部２０から動作開始情報が入力されると動作を開始し、標的追尾制御部３０から入力される標的までの距離と標的の俯角とに基づいて、入力される標的の俯角を基準として含む複数の俯角で送受波機５０に音波を送波させてその反射波を受波させる制御を行うと共に、送受波機５０から入力される受波した反射波の強度を示すサンプリング値を標的追尾制御部３０と表示情報変換部６０とに出力する。
【００２０】
送受波機５０は、送波部５１と受波部５２とを含んで構成される。送波部５１は、送受波機制御部４０の制御に従い、標的までの距離に応じた送波間隔で、送受波機制御部４０から入力された俯角で周囲３６０度の方向、すなわち、図２に示されるように円錐の表面方向に音波を送波する。受波部５２は、送波された音波が魚群で反射された反射波を俯角ごとに受波し、受波した反射波の強度を示す複数のサンプリング値を送受波機制御部４０に出力する。
また、受波部５２は、図２に示すように、狭い範囲からの反射波を受波する指向性を有するペンシル形状の受波ビームＰを海面に対する垂直方向（ｚ軸方向）に対して時計回り又は反時計回りに回転させて反射波を受波する。
【００２１】
また、送受波機５０は、漁船の船底に設置されるか、ブイなどに繋がれて海中（水中）に設置される。ここで、俯角は、図２に示されるように、海面（水面）に対してθの角度で示される角度である。送波間隔は、水中における音速（約１５００メートル毎秒）に基づいて定められ、例えば、標的まで１５００メートルの場合、音波が送波部５１から魚群までの伝搬に要する時間と、魚群から受波部５２までの伝搬に要する時間と合わせた約２秒、あるいは、魚群を中心とした範囲で探索する場合、魚群までの距離と探索する範囲とに応じた時間となる。
【００２２】
表示情報変換部６０は、送受波機制御部４０から出力されたサンプリング値を１フレームごとに表示装置７０に表示するための画像情報に変換して表示装置７０に出力する。このとき、表示情報変換部６０は、魚群の標的像を認識しやすい色を配置したりなどしてもよい。
表示装置７０は、例えば、ＣＲＴ（Cathode Ray Tube）や液晶表示装置などであり、表示情報変換部６０が出力した画像情報を表示して出力する。
【００２３】
次に、スキャニングソナー装置１の動作について説明する。
まず、スキャニングソナー装置１において、入力装置１０は、ユーザの操作により魚群の追尾開始を示す動作開始情報、俯角を示す俯角情報、及び、標的像抽出部３１が標的を抽出する際に用いる閾値を含む初期閾値情報が入力される。操作制御部２０は、入力装置１０に入力された動作開始情報を標的追尾制御部３０、送受波機制御部４０、及び、表示情報変換部６０に出力し、入力装置１０に入力された俯角情報と閾値情報とを標的追尾制御部３０に出力する。俯角制御部３４は、操作制御部２０から入力された初期閾値情報を閾値情報として記憶して、操作制御部２０から入力された俯角情報を送受波機制御部４０に出力する。
上述のようにして、スキャニングソナー装置１は、魚群追尾処理を開始する。
【００２４】
続いて、魚群追尾処理を説明する。送受波機制御部４０は、俯角制御部３４から入力された俯角情報に含まれる俯角θを基準とした３つの俯角情報（俯角θ−Δθ、θ、θ＋Δθ）を送受波機５０に出力する。送受波機５０は、送受波機制御部４０から入力された俯角情報に応じた俯角で順次音波を送波し、受波した当該音波の反射波に対する複数のサンプリング値を送受波機制御部４０に出力する。送受波機制御部４０は、送受波機５０から入力された３つの異なる俯角それぞれに対する１フレーム分のサンプリング値を標的追尾制御部３０と表示情報変換部６０とに出力する。
表示情報変換部６０は、送受波機制御部４０から入力される３つのフレーム分のサンプリング値を極座標変換やその他の情報を重畳した表示情報に変換し、表示装置７０に出力して表示させる。
【００２５】
標的像抽出部３１は、操作制御部２０から入力された閾値情報に基づいて、送受波機制御部４０から入力された３フレーム（３つの俯角）分のサンプリング値からフレームごとに標的像を抽出して、特徴量抽出部３２に出力する。特徴量抽出部３２は、３つのフレームそれぞれの特徴量Ｆを式（１）により算出して、移動方向判定部３３に出力する。移動方向判定部３３は、入力された３つの特徴量Ｆの大小比較を行い、魚群の移動方向を判定して判定結果を俯角制御部３４に出力する。
俯角制御部３４は、移動方向判定部３３の判定結果により、現在の俯角情報を更新して、送受波機制御部４０に出力する。以後、スキャニングソナー装置１は、魚群追尾処理を繰り返して行う。
【００２６】
上述のように動作することにより、スキャニングソナー装置１は、複数の異なる俯角ごとの特徴量Ｆを算出して、算出した複数の特徴量Ｆの大小比較を行うことで魚群（標的）のなかでサンプリング値が大きい領域、すなわち、魚が密集している領域の俯角方向に対する移動の向きを判断し、標的とした魚群の俯角方向の追尾を行うことができる。
【００２７】
なお、標的像抽出部３１の標的像抽出において、複数の標的像が抽出される場合、ユーザが、入力装置１０を介して選択した複数の標的像のうちいずれか１つに対して移動方向の判定を行うようにしてもよい。
また、標的像が小さい場合、揺らぎをさらに軽減するために、３つの異なる俯角（θ−Δθ、θ、θ＋Δθ）それぞれに対するフレームを１フレームでなく、複数フレーム分スキャンして、３つの俯角ごとに平均値を算出しても良い。その場合、送受波機制御部４０は、複数フレームのサンプリング値を取得するために、送受波機５０に当該回数分の送波と受波とを行わせる制御を行う。これにより、空間的平均と時間的平均とを併用して特徴量Ｆを算出することになり、移動方向の判定制度をさらに改善することができる。なお、各フレームの平均値を算出するための探査回数は、空間的平均を用いない場合に比べ、減らすことができるので、時間的平均値の算出を行ったとしても、従来の算出方法より短時間で特徴量Ｆを算出することが可能である。
【００２８】
＜第２実施形態＞
第２実施形態は、特徴量抽出部３２が行う特徴量Ｆの算出方法が第１実施形態と異なり、次式（２）を用いて算出する。式（２）では、標的像の面積Ａｐを標的像内に含まれるサンプル数ｎｐで代用することによって計算を簡略化した。
【００２９】
【数２】
【００３０】
ここで、式（２）における、ｎｐは、標的像に含まれるサンプル数を表し、ｘ（ｐｉ）は、ｉ番目のサンプリング値を表す。すなわち、標的像に含まれる全サンプリング値の平均と標的像の面積とを乗じた数値を特徴量Ｆとする。すなわち、特徴量Ｆは、サンプリング値を総和した数値となる。
これにより、第２実施形態のスキャニングソナー装置１は、第１実施形態のスキャニングソナー装置１に比べ、特徴量抽出部３２が標的像の面積を算出せずに特徴量Ｆの算出を行うので、算出に要する時間を短縮することが可能となる。
【００３１】
上述のように、第１実施形態又は第２実施形態におけるスキャニングソナー装置１において、観測点である送受波機５０が移動又は浮遊し、観測対象である魚群が移動し、かつ、音波の伝搬状態（海水の温度及び潮流）が時間経過と共に変化する環境では、スペックルによる反射波の強度の揺らぎは、観測される各サンプルにおいて時間的にランダムなノイズと考えられる。そこで、平均数をｎｔとする各サンプルの時間的な平均を算出することで揺らぎを１／√（ｎｔ）倍に低減することができる。
【００３２】
一方、同様の環境において、スペックルによる反射波の強度の揺らぎは、空間的にも発生するため、時間的ノイズと空間的ノイズとの間にエルゴード性（次式（３））を仮定し、時間平均の代わりに空間平均を用いることで、スペックルによる揺らぎを緩和して、特徴量Ｆを算出することが上述の第１実施形態及び第２実施形態のスキャニングソナー装置１の特徴の１つである。
特徴量Ｆの算出において、空間平均を用いることで、特徴量抽出部３２が特徴量Ｆを算出する際に要する時間を短縮すると共に、サンプリング値の揺らぎに対して空間的平均値を用いることで低減させて、標的の俯角方向に対する移動の向き検出の精度を改善することが可能である。
【００３３】
【数３】
【００３４】
上述の標的追尾制御部３０は内部に、コンピュータシステムを有していてもよい。その場合、上述した標的像抽出部３１、特徴量抽出部３２、移動方向判定部３３、俯角制御部３４が行う処理過程は、プログラムの形式でコンピュータ読み取り可能な記録媒体に記憶されており、このプログラムをコンピュータが読み出して実行することによって、上記処理が行われることになる。ここでコンピュータ読み取り可能な記録媒体とは、磁気ディスク、光磁気ディスク、ＣＤ−ＲＯＭ、ＤＶＤ−ＲＯＭ、半導体メモリ等をいう。また、このコンピュータプログラムを通信回線によってコンピュータに配信し、この配信を受けたコンピュータが当該プログラムを実行するようにしても良い。
【図面の簡単な説明】
【００３５】
【図１】本実施形態によるスキャニングソナー装置の構成を示すブロック図である。
【図２】従来例におけるスキャニングソナー装置における反射波の受波の概念図である。
【図３】従来例における標的像の表示例を示す概略図である。
【符号の説明】
【００３６】
１…スキャニングソナー装置
１０…入力装置、２０…操作制御部
３０…標的追尾制御部、３１…標的像抽出部、３２…特徴量抽出部
３３…移動方向判定部、３４…俯角制御部
４０…送受波機制御部、５０…送受波機、５１…送波部、５２…受波部
６０…表示情報変換部、７０…表示装置
１Ａ…スキャニングソナー装置
【技術分野】
【０００１】
本発明は、スキャニングソナー装置、主に設定した標的を自動追尾するスキャニングソナー装置に関する。
【背景技術】
【０００２】
従来、漁労などに用いられるスキャニングソナー装置では、例えば、図２に示すように船体などに設けられ水中に位置する送受波機５０が、一定の周波数、例えば、２０ｋＨｚから３００ｋＨｚまでの間の周波数の音波を、ユーザによって設定される海面からθの角度（俯角という）に送波して、受波したエコー（反射波）の強度情報（サンプリング値）を用いて標的である魚群などの探査を行い、その位置検出をしている。また、スキャニングソナー装置は、反射波の強度を予め定めた閾値を超える領域を図３（ａ）又は図３（ｂ）に示すように表示装置に画像として表示する。ここで、図３（ａ）は、縦軸方向及び横軸方向共に中心から周囲へ向かういずれの方向に対しても俯角θを有する円錐体の傘状の表面における送受波機５０からの距離を示している。図３（ｂ）は、縦軸方向が俯角θを有する円錐体の傘状の表面における送受波機５０からの距離を示し、横軸方向が船首又は船尾を基準とし海面に対する垂直な軸ｚ（図２）の回転方向の角度φを示している。
【０００３】
ここで、送受波機５０は、図２に示すように狭い範囲からの反射波を受波する指向性を有するペンシル形状の受波ビームＰを、海面に対する垂直方向（ｚ軸方向）に対して時計回り又は反時計回りに回転させ全周３６０度からの反射波を受信して円錐体の傘状の表面に対してスキャンを行う。
このようなスキャニングソナー装置において、標的の移動方向を検出する際には、異なる俯角に対して複数回の送波及び受波を行い、受波した反射波の強度分布や、標的像の面積、及び、寸法（図３（ｂ）におけるΔφ及びΔＲ）を特徴量として抽出して標的の俯角方向に対する移動の向きの推定を行うものもある（特許文献１）。
【０００４】
ところで、反射波を受波する送受波機の移動及び波による揺れ、標的である魚などの移動、及び、潮流や温度により音波の伝搬路の状態が時間的変化することによって生じる外乱により、受波する反射波の強度は、主にスペックルと同様の原理により空間的及び時間的に揺らぐという現象がある。このため、空間的及び時間的に揺らいだ反射波の強度から標的像を抽出すると、抽出した標的像の形状が、実際の標的像に比べ歪むと共に、抽出した標的像の分布が、実際の標的の分布に比べ揺らいで安定しない。上述のようなスキャニングソナー装置においては、複数の異なる俯角で音波を送波し、その反射波を受波して標的を立体的に捉えようとする場合、上述の揺らぎの影響を低減させるために、複数の異なる俯角それぞれにおいて複数回の探査を行って得られた標的像の平均値を算出して、時間的な揺らぎの影響を緩和した標的像の抽出と、特徴量の算出とが行われている。
【特許文献１】特開２００８−２０３２２７号公報
【発明の開示】
【発明が解決しようとする課題】
【０００５】
しかしながら、平均値を算出するために複数回探査を行おうとすると、特徴量算出に要する時間が長くなり、その間に標的が移動するなどして、実際の標的の俯角方向に対する移動の向きと、得られた特徴量から得られる標的の俯角方向に対する移動の向きとの乖離が大きくなり、標的の俯角方向に対する移動の向き検出の精度が悪くなるという問題がある。
【０００６】
本発明は、上記問題を解決すべくなされたもので、その目的は、外乱による反射波の強度の揺らぎによる影響を低減し、標的の俯角方向に対する移動の向き検出の精度を改善したスキャニングソナー装置を提供することにある。
【課題を解決するための手段】
【０００７】
上記問題を解決するために、本発明は、音波を標的に送波し、その反射波を受波して得られる反射波の強度に基づく複数のサンプリング値から前記標的の形状を示す標的像を抽出する標的像抽出部と、前記標的像抽出部が抽出した前記標的像に含まれる前記サンプリング値から前記標的像の特徴量を算出する特徴量算出部と、前記標的像抽出部が算出した複数の俯角ごとの前記標的像から前記特徴量抽出部によって算出された当該複数の俯角それぞれに対応した前記特徴量を比較して前記標的の俯角方向に対する移動の向きを判定する移動方向判定部と、を具備し、前記特徴量は、前記標的像に含まれる前記サンプリング値の平均値と前記標的像の面積とを乗じた数値であることを特徴とするスキャニングソナー装置である。
【０００８】
また、本発明は、音波を標的に送波し、その反射波を受波して得られた反射波の強度に基づく複数のサンプリング値から前記標的の形状を示す標的像を抽出する標的像抽出部と、前記標的像抽出部が抽出した前記標的像に含まれる前記サンプリング値から前記標的像の特徴量を算出する特徴量算出部と、前記標的像抽出部が算出した複数の俯角ごとの前記標的像から前記特徴量抽出部によって算出された当該複数の俯角それぞれに対応した前記特徴量を比較して、前記標的の俯角方向に対する移動の向きを判定する移動方向判定部と、を具備し、前記特徴量は、前記標的像を示す反射波の前記サンプリング値を総和した数値であることを特徴とするスキャニングソナー装置である。
【０００９】
また、本発明は、上記に記載の発明において、前記標的像の面積が予め定めた数値より小さいとき、前記特徴量算出部は、前記標的像抽出部が抽出した前記複数の異なる俯角ごとの複数の前記標的像から前記複数の異なる俯角ごとに複数の前記特徴量を算出し、算出した前記複数の特徴量から前記複数の異なる俯角ごとに前記特徴量の平均値を算出して、前記標的の俯角方向に対する移動の向きを判定することを特徴とする。
【発明の効果】
【００１０】
この発明によれば、スキャニングソナー装置は、水面に対して複数の異なる俯角で送波された音波が標的で反射した反射波（エコー）に基づいて、複数の異なる俯角ごとに標的の形状を示す標的像を示す反射波のサンプリング値を用いて特徴量を算出し、算出した特徴量を俯角間で比較することで標的の俯角方向に対する移動の向きを推定する。このとき、特徴量算出部は、空間平均として、標的像を示す反射波のサンプリング値の平均値と標的像の面積とを乗じた値、あるいは、標的像を示す反射波のサンプリング値を総和した値を特徴量として算出する。外乱による反射波の強度の揺らぎは、時間平均を用いずにエルゴード性を適用し空間平均を用いることで低減させて、複数回の送波及び受波を繰り返して時間平均を算出する代わりに空間平均を算出することで移動方向判定時間の短縮を行うと共に、反射波の強度の揺らぎを低減させることができる。これにより、スキャニングソナー装置は、探査回数を増やすことなく反射波の強度揺らぎを低減させ、標的の俯角方向に対する移動の向き検出の精度を改善することが可能となる。
【発明を実施するための最良の形態】
【００１１】
以下、本発明の実施形態によるスキャニングソナー装置を図面を参照して説明する。
【００１２】
＜第１実施形態＞
図１は、第１実施形態におけるスキャニングソナー装置１の構成を示す概略ブロック図である。以下、スキャニングソナー装置１は、魚群を標的として漁労に用いられる場合を例として説明する。
スキャニングソナー装置１は、入力装置１０、操作制御部２０、標的追尾制御部３０、送受波機制御部４０、送受波機５０、表示情報変換部６０、及び、表示装置７０を具備している。
入力装置１０は、ユーザの操作によりスキャニングソナー装置１への入力が行える、例えば、キーボード、マウス、トラックボールなどである。操作制御部２０は、入力装置１０から入力された情報の種類に応じて標的追尾制御部３０、送受波機制御部４０、表示情報変換部６０に出力する。ここで、入力装置１０から入力される情報としては、スキャニングソナー装置１に対して標的の追尾の開始を示す動作開始情報、音波を送波する俯角を示す俯角情報、標的をスキャンする範囲（レンジ）情報、標的（魚群）の設定情報、表示装置７０に表示するスキャン範囲を示す情報などである。
【００１３】
標的追尾制御部３０は、標的像抽出部３１、特徴量抽出部３２、移動方向判定部３３、及び俯角制御部３４を含んで構成され、操作制御部２０から動作開始情報が入力されると動作を開始し、操作制御部２０から入力される標的の設定情報及び俯角情報を初期値として標的の追尾処理を開始する。また、標的追尾制御部３０は、標的までの距離と標的の位置に応じた俯角とを含むデータを送受波機制御部４０に出力し、送受波機制御部４０から反射波（エコー）を受信することで得られた複数のサンプリング値が入力される。
ここで、サンプリング値は、図３（ｂ）に示される海面の垂直方向に対して時計回り又は反時計回り方向の角度φと距離Ｒとを用いて表されるＲ−φ座標系の各座標における受波した反射波の強さに基づく値（強度情報）である。なお、角度φは、船首や船尾など船体の一部を基準とし０°として定められる。また、角度φに対して０°から３６０°までの一周分の全探査距離に渡るサンプリング値を１フレームという。
【００１４】
標的像抽出部３１は、送受波機制御部４０から入力される１フレーム分のサンプリング値に対して、標的の形状を示す標的像の抽出を行う。ここで、標的像抽出部３１は、操作制御部２０を介して入力されるユーザの指定した閾値を基準とし、当該閾値を超えるサンプリング値により構成される領域を標的像として抽出を行う。例えば、標的像は、図３（ｂ）に図示される標的Ｔの領域に含まれるＲ−φ座標系の情報として出力される。また、ユーザの指定する閾値とは、漁労の対象となる魚種や観測環境ごとに定められる値であり、ユーザの経験などに基づいて定められる。また、標的像抽出部３１は、送受波機５０から抽出した標的像までの距離を送受波機制御部４０に出力する。
特徴量抽出部３２は、標的像抽出部３１が抽出した標的像の領域に含まれるサンプリング値から標的の移動方向判断に用いる特徴量Ｆを算出する。特徴量Ｆは、次式（１）により求める。
【００１５】
【数１】
【００１６】
ここで、式（１）における、Ａｐは、標的像の面積を表し、ｎｐは、標的像に含まれるサンプル数を表し、ｘ（ｐｉ）は、ｉ番目のサンプリング値を表す。すなわち、標的像に含まれる全サンプリング値の平均と標的像の面積とを乗じた数値が特徴量Ｆである。なお、標的像の面積は、Ｒ−φ座標から算出される。
【００１７】
移動方向判定部３３は、標的像抽出部３１から入力される複数の異なる俯角ごとの特徴量Ｆを比較して、魚群の俯角方向における移動の向きを判断する。判断方法としては、俯角θ−Δθ、θ、θ＋Δθ（Δθ＞０）それぞれの特徴量Ｆを比較し、俯角θ−Δθの特徴量Ｆが最も大きいとき、魚群は俯角が小さくなる方向に移動していると判断し、俯角θの特徴量Ｆが最も大きいとき、魚群は移動していないと判断し、俯角θ＋Δθの特徴量Ｆが最も大きいとき、魚群は俯角が大きくなる方向に移動していると判断する。すなわち、移動方向判定部３３は、特徴量Ｆの大小比較により魚群の俯角方向に対する移動の向き判断する。魚群の魚の量は、平均的な反射波の強度と相関が高いのでこの値が大きい領域を追うことで、魚群における魚が密集している領域を追尾することができ、判定の誤りによる魚群を見失うことを防ぐことが可能になる。
なお、移動方向判定部３３が、３つの異なる俯角ごとの特徴量Ｆを用いて魚群の移動の向きを判断する例を示したが、２つ以上の異なる俯角ごとの特徴量Ｆを用いて魚群の移動の向きを判断するようにしてもよい。
【００１８】
俯角制御部３４は、移動方向判定部３３により判断された魚群の移動の向きに応じて、送受波機５０が送波及び受波する俯角情報を更新して、送受波機制御部４０に出力する。
例えば、移動方向判定部３３が上述のように判断を行い、俯角が小さくなる方向に移動しているとき、現在の俯角情報θに対して俯角θ−Δθを次の俯角情報として更新する。
【００１９】
送受波機制御部４０は、操作制御部２０から動作開始情報が入力されると動作を開始し、標的追尾制御部３０から入力される標的までの距離と標的の俯角とに基づいて、入力される標的の俯角を基準として含む複数の俯角で送受波機５０に音波を送波させてその反射波を受波させる制御を行うと共に、送受波機５０から入力される受波した反射波の強度を示すサンプリング値を標的追尾制御部３０と表示情報変換部６０とに出力する。
【００２０】
送受波機５０は、送波部５１と受波部５２とを含んで構成される。送波部５１は、送受波機制御部４０の制御に従い、標的までの距離に応じた送波間隔で、送受波機制御部４０から入力された俯角で周囲３６０度の方向、すなわち、図２に示されるように円錐の表面方向に音波を送波する。受波部５２は、送波された音波が魚群で反射された反射波を俯角ごとに受波し、受波した反射波の強度を示す複数のサンプリング値を送受波機制御部４０に出力する。
また、受波部５２は、図２に示すように、狭い範囲からの反射波を受波する指向性を有するペンシル形状の受波ビームＰを海面に対する垂直方向（ｚ軸方向）に対して時計回り又は反時計回りに回転させて反射波を受波する。
【００２１】
また、送受波機５０は、漁船の船底に設置されるか、ブイなどに繋がれて海中（水中）に設置される。ここで、俯角は、図２に示されるように、海面（水面）に対してθの角度で示される角度である。送波間隔は、水中における音速（約１５００メートル毎秒）に基づいて定められ、例えば、標的まで１５００メートルの場合、音波が送波部５１から魚群までの伝搬に要する時間と、魚群から受波部５２までの伝搬に要する時間と合わせた約２秒、あるいは、魚群を中心とした範囲で探索する場合、魚群までの距離と探索する範囲とに応じた時間となる。
【００２２】
表示情報変換部６０は、送受波機制御部４０から出力されたサンプリング値を１フレームごとに表示装置７０に表示するための画像情報に変換して表示装置７０に出力する。このとき、表示情報変換部６０は、魚群の標的像を認識しやすい色を配置したりなどしてもよい。
表示装置７０は、例えば、ＣＲＴ（Cathode Ray Tube）や液晶表示装置などであり、表示情報変換部６０が出力した画像情報を表示して出力する。
【００２３】
次に、スキャニングソナー装置１の動作について説明する。
まず、スキャニングソナー装置１において、入力装置１０は、ユーザの操作により魚群の追尾開始を示す動作開始情報、俯角を示す俯角情報、及び、標的像抽出部３１が標的を抽出する際に用いる閾値を含む初期閾値情報が入力される。操作制御部２０は、入力装置１０に入力された動作開始情報を標的追尾制御部３０、送受波機制御部４０、及び、表示情報変換部６０に出力し、入力装置１０に入力された俯角情報と閾値情報とを標的追尾制御部３０に出力する。俯角制御部３４は、操作制御部２０から入力された初期閾値情報を閾値情報として記憶して、操作制御部２０から入力された俯角情報を送受波機制御部４０に出力する。
上述のようにして、スキャニングソナー装置１は、魚群追尾処理を開始する。
【００２４】
続いて、魚群追尾処理を説明する。送受波機制御部４０は、俯角制御部３４から入力された俯角情報に含まれる俯角θを基準とした３つの俯角情報（俯角θ−Δθ、θ、θ＋Δθ）を送受波機５０に出力する。送受波機５０は、送受波機制御部４０から入力された俯角情報に応じた俯角で順次音波を送波し、受波した当該音波の反射波に対する複数のサンプリング値を送受波機制御部４０に出力する。送受波機制御部４０は、送受波機５０から入力された３つの異なる俯角それぞれに対する１フレーム分のサンプリング値を標的追尾制御部３０と表示情報変換部６０とに出力する。
表示情報変換部６０は、送受波機制御部４０から入力される３つのフレーム分のサンプリング値を極座標変換やその他の情報を重畳した表示情報に変換し、表示装置７０に出力して表示させる。
【００２５】
標的像抽出部３１は、操作制御部２０から入力された閾値情報に基づいて、送受波機制御部４０から入力された３フレーム（３つの俯角）分のサンプリング値からフレームごとに標的像を抽出して、特徴量抽出部３２に出力する。特徴量抽出部３２は、３つのフレームそれぞれの特徴量Ｆを式（１）により算出して、移動方向判定部３３に出力する。移動方向判定部３３は、入力された３つの特徴量Ｆの大小比較を行い、魚群の移動方向を判定して判定結果を俯角制御部３４に出力する。
俯角制御部３４は、移動方向判定部３３の判定結果により、現在の俯角情報を更新して、送受波機制御部４０に出力する。以後、スキャニングソナー装置１は、魚群追尾処理を繰り返して行う。
【００２６】
上述のように動作することにより、スキャニングソナー装置１は、複数の異なる俯角ごとの特徴量Ｆを算出して、算出した複数の特徴量Ｆの大小比較を行うことで魚群（標的）のなかでサンプリング値が大きい領域、すなわち、魚が密集している領域の俯角方向に対する移動の向きを判断し、標的とした魚群の俯角方向の追尾を行うことができる。
【００２７】
なお、標的像抽出部３１の標的像抽出において、複数の標的像が抽出される場合、ユーザが、入力装置１０を介して選択した複数の標的像のうちいずれか１つに対して移動方向の判定を行うようにしてもよい。
また、標的像が小さい場合、揺らぎをさらに軽減するために、３つの異なる俯角（θ−Δθ、θ、θ＋Δθ）それぞれに対するフレームを１フレームでなく、複数フレーム分スキャンして、３つの俯角ごとに平均値を算出しても良い。その場合、送受波機制御部４０は、複数フレームのサンプリング値を取得するために、送受波機５０に当該回数分の送波と受波とを行わせる制御を行う。これにより、空間的平均と時間的平均とを併用して特徴量Ｆを算出することになり、移動方向の判定制度をさらに改善することができる。なお、各フレームの平均値を算出するための探査回数は、空間的平均を用いない場合に比べ、減らすことができるので、時間的平均値の算出を行ったとしても、従来の算出方法より短時間で特徴量Ｆを算出することが可能である。
【００２８】
＜第２実施形態＞
第２実施形態は、特徴量抽出部３２が行う特徴量Ｆの算出方法が第１実施形態と異なり、次式（２）を用いて算出する。式（２）では、標的像の面積Ａｐを標的像内に含まれるサンプル数ｎｐで代用することによって計算を簡略化した。
【００２９】
【数２】
【００３０】
ここで、式（２）における、ｎｐは、標的像に含まれるサンプル数を表し、ｘ（ｐｉ）は、ｉ番目のサンプリング値を表す。すなわち、標的像に含まれる全サンプリング値の平均と標的像の面積とを乗じた数値を特徴量Ｆとする。すなわち、特徴量Ｆは、サンプリング値を総和した数値となる。
これにより、第２実施形態のスキャニングソナー装置１は、第１実施形態のスキャニングソナー装置１に比べ、特徴量抽出部３２が標的像の面積を算出せずに特徴量Ｆの算出を行うので、算出に要する時間を短縮することが可能となる。
【００３１】
上述のように、第１実施形態又は第２実施形態におけるスキャニングソナー装置１において、観測点である送受波機５０が移動又は浮遊し、観測対象である魚群が移動し、かつ、音波の伝搬状態（海水の温度及び潮流）が時間経過と共に変化する環境では、スペックルによる反射波の強度の揺らぎは、観測される各サンプルにおいて時間的にランダムなノイズと考えられる。そこで、平均数をｎｔとする各サンプルの時間的な平均を算出することで揺らぎを１／√（ｎｔ）倍に低減することができる。
【００３２】
一方、同様の環境において、スペックルによる反射波の強度の揺らぎは、空間的にも発生するため、時間的ノイズと空間的ノイズとの間にエルゴード性（次式（３））を仮定し、時間平均の代わりに空間平均を用いることで、スペックルによる揺らぎを緩和して、特徴量Ｆを算出することが上述の第１実施形態及び第２実施形態のスキャニングソナー装置１の特徴の１つである。
特徴量Ｆの算出において、空間平均を用いることで、特徴量抽出部３２が特徴量Ｆを算出する際に要する時間を短縮すると共に、サンプリング値の揺らぎに対して空間的平均値を用いることで低減させて、標的の俯角方向に対する移動の向き検出の精度を改善することが可能である。
【００３３】
【数３】
【００３４】
上述の標的追尾制御部３０は内部に、コンピュータシステムを有していてもよい。その場合、上述した標的像抽出部３１、特徴量抽出部３２、移動方向判定部３３、俯角制御部３４が行う処理過程は、プログラムの形式でコンピュータ読み取り可能な記録媒体に記憶されており、このプログラムをコンピュータが読み出して実行することによって、上記処理が行われることになる。ここでコンピュータ読み取り可能な記録媒体とは、磁気ディスク、光磁気ディスク、ＣＤ−ＲＯＭ、ＤＶＤ−ＲＯＭ、半導体メモリ等をいう。また、このコンピュータプログラムを通信回線によってコンピュータに配信し、この配信を受けたコンピュータが当該プログラムを実行するようにしても良い。
【図面の簡単な説明】
【００３５】
【図１】本実施形態によるスキャニングソナー装置の構成を示すブロック図である。
【図２】従来例におけるスキャニングソナー装置における反射波の受波の概念図である。
【図３】従来例における標的像の表示例を示す概略図である。
【符号の説明】
【００３６】
１…スキャニングソナー装置
１０…入力装置、２０…操作制御部
３０…標的追尾制御部、３１…標的像抽出部、３２…特徴量抽出部
３３…移動方向判定部、３４…俯角制御部
４０…送受波機制御部、５０…送受波機、５１…送波部、５２…受波部
６０…表示情報変換部、７０…表示装置
１Ａ…スキャニングソナー装置
【特許請求の範囲】
【請求項１】
音波を標的に送波し、その反射波を受波して得られる反射波の強度に基づく複数のサンプリング値から前記標的の形状を示す標的像を抽出する標的像抽出部と、
前記標的像抽出部が抽出した前記標的像に含まれる前記サンプリング値から前記標的像の特徴量を算出する特徴量算出部と、
前記標的像抽出部が算出した複数の俯角ごとの前記標的像から前記特徴量抽出部によって算出された当該複数の俯角それぞれに対応した前記特徴量を比較して前記標的の俯角方向に対する移動の向きを判定する移動方向判定部と、
を具備し、
前記特徴量は、
前記標的像に含まれる前記サンプリング値の平均値と前記標的像の面積とを乗じた数値である
ことを特徴とするスキャニングソナー装置。
【請求項２】
音波を標的に送波し、その反射波を受波して得られる反射波の強度に基づく複数のサンプリング値から前記標的の形状を示す標的像を抽出する標的像抽出部と、
前記標的像抽出部が抽出した前記標的像に含まれる前記サンプリング値から前記標的像の特徴量を算出する特徴量算出部と、
前記標的像抽出部が算出した複数の俯角ごとの前記標的像から前記特徴量抽出部によって算出された当該複数の俯角それぞれに対応した前記特徴量を比較して、前記標的の俯角方向に対する移動の向きを判定する移動方向判定部と、
を具備し、
前記特徴量は、
前記標的像を示す反射波の前記サンプリング値を総和した数値である
ことを特徴とするスキャニングソナー装置。
【請求項３】
前記標的像の面積が予め定めた数値より小さいとき、
前記特徴量算出部は、前記標的像抽出部が抽出した前記複数の異なる俯角ごとの複数の前記標的像から前記複数の異なる俯角ごとに複数の前記特徴量を算出し、算出した前記複数の特徴量から前記複数の異なる俯角ごとに前記特徴量の平均値を算出して、前記標的の俯角方向に対する移動の向きを判定する
ことを特徴とする請求項１又は請求項２のスキャニングソナー装置。
【請求項１】
音波を標的に送波し、その反射波を受波して得られる反射波の強度に基づく複数のサンプリング値から前記標的の形状を示す標的像を抽出する標的像抽出部と、
前記標的像抽出部が抽出した前記標的像に含まれる前記サンプリング値から前記標的像の特徴量を算出する特徴量算出部と、
前記標的像抽出部が算出した複数の俯角ごとの前記標的像から前記特徴量抽出部によって算出された当該複数の俯角それぞれに対応した前記特徴量を比較して前記標的の俯角方向に対する移動の向きを判定する移動方向判定部と、
を具備し、
前記特徴量は、
前記標的像に含まれる前記サンプリング値の平均値と前記標的像の面積とを乗じた数値である
ことを特徴とするスキャニングソナー装置。
【請求項２】
音波を標的に送波し、その反射波を受波して得られる反射波の強度に基づく複数のサンプリング値から前記標的の形状を示す標的像を抽出する標的像抽出部と、
前記標的像抽出部が抽出した前記標的像に含まれる前記サンプリング値から前記標的像の特徴量を算出する特徴量算出部と、
前記標的像抽出部が算出した複数の俯角ごとの前記標的像から前記特徴量抽出部によって算出された当該複数の俯角それぞれに対応した前記特徴量を比較して、前記標的の俯角方向に対する移動の向きを判定する移動方向判定部と、
を具備し、
前記特徴量は、
前記標的像を示す反射波の前記サンプリング値を総和した数値である
ことを特徴とするスキャニングソナー装置。
【請求項３】
前記標的像の面積が予め定めた数値より小さいとき、
前記特徴量算出部は、前記標的像抽出部が抽出した前記複数の異なる俯角ごとの複数の前記標的像から前記複数の異なる俯角ごとに複数の前記特徴量を算出し、算出した前記複数の特徴量から前記複数の異なる俯角ごとに前記特徴量の平均値を算出して、前記標的の俯角方向に対する移動の向きを判定する
ことを特徴とする請求項１又は請求項２のスキャニングソナー装置。
【図１】
【図２】
【図３】
【図２】
【図３】
【公開番号】特開２０１０−１４５２２４（Ｐ２０１０−１４５２２４Ａ）
【公開日】平成２２年７月１日（２０１０．７．１）
【国際特許分類】
【出願番号】特願２００８−３２２３５５（Ｐ２００８−３２２３５５）
【出願日】平成２０年１２月１８日（２００８．１２．１８）
【出願人】（０００００４３３０）日本無線株式会社 (1,186)
【Ｆターム（参考）】
【公開日】平成２２年７月１日（２０１０．７．１）
【国際特許分類】
【出願日】平成２０年１２月１８日（２００８．１２．１８）
【出願人】（０００００４３３０）日本無線株式会社 (1,186)
【Ｆターム（参考）】
[ Back to top ]