JP2016153728A - レーダ信号処理装置及びレーダ装置 - Google Patents
レーダ信号処理装置及びレーダ装置 Download PDFInfo
- Publication number
- JP2016153728A JP2016153728A JP2015031334A JP2015031334A JP2016153728A JP 2016153728 A JP2016153728 A JP 2016153728A JP 2015031334 A JP2015031334 A JP 2015031334A JP 2015031334 A JP2015031334 A JP 2015031334A JP 2016153728 A JP2016153728 A JP 2016153728A
- Authority
- JP
- Japan
- Prior art keywords
- radar
- precipitation
- radio wave
- received power
- plot
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
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/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
【課題】 探査地域における降水の強度を算出し、探査地域における目標の確認を容易にするレーダ信号処理装置およびレーダ装置を得る。
【解決手段】探査地域を分割するメッシュ毎のレーダ反射電波の受信電力から閾値超過メッシュを検出する検出部と、隣接する閾値超過メッシュ同士をまとめた検出プロットを作成する検出プロット作成部と、検出プロットの面積を予め決められた下限値及び上限値と比較し、目標プロットであるか降水プロットであるかを判別するクラッタ判別部と、探査地域に設定された区画において、降水プロットの合計面積が占める割合を誤検出頻度として算出し、降水反射電波の受信電力の確率分布関数から、受信電力が電力閾値を超過する確率が誤検出頻度になる降水反射電波の推定平均受信電力を求める反射強度推定部と、推定平均受信電力とレーダ方程式により算出したレーダ反射因子から降水強度を推定する降水強度推定部とを備える。
【選択図】 図2
【解決手段】探査地域を分割するメッシュ毎のレーダ反射電波の受信電力から閾値超過メッシュを検出する検出部と、隣接する閾値超過メッシュ同士をまとめた検出プロットを作成する検出プロット作成部と、検出プロットの面積を予め決められた下限値及び上限値と比較し、目標プロットであるか降水プロットであるかを判別するクラッタ判別部と、探査地域に設定された区画において、降水プロットの合計面積が占める割合を誤検出頻度として算出し、降水反射電波の受信電力の確率分布関数から、受信電力が電力閾値を超過する確率が誤検出頻度になる降水反射電波の推定平均受信電力を求める反射強度推定部と、推定平均受信電力とレーダ方程式により算出したレーダ反射因子から降水強度を推定する降水強度推定部とを備える。
【選択図】 図2
Description
この発明は、レーダ電波の反射電波から、探査地域に存在する目標物を検出するレーダ信号処理装置及びレーダ装置に関する。
探査地域内の目標検出を行うレーダ装置は、レーダ電波を探査地域に照射し、探査地域から反射された反射電波を受信し、受信した反射電波から車両、航空機、船舶等の目標物を検出する。探査地域では、レーダ装置が検出対象とする目標物の他に、雪や雨などの降水もレーダ電波を反射する。レーダ装置は、降水による反射についても、目標物からの反射と同様に検出してしまい、検出した結果は、クラッタと呼ばれる誤目標となる。このため、雪や雨によるレーダ電波の反射が甚だしいときは、レーダ装置の検出対象とする目標は、クラッタに埋もれてしまい、レーダ表示画面上で目標を確認することが困難になる。
レーダ装置において、誤目標の発生を抑えるためには、例えば、受信部の増幅器の入出力特性を制御して反射信号の増幅率を抑える方法が有効であるが、この方法は反面、目標の検出も抑圧してしまう。このため、目標検出の抑圧を必要最小限にするよう、降水の強度を検出し、検出した降水の強度にあわせてレーダ装置の各種特性を制御することが必要となる。
降水の強度にあわせて誤目標の抑圧を行なう従来のレーダ装置としては、検出した目標を追尾して誤目標を判定し、誤目標と判定された数をもとに、レーダ装置の増幅器の入出力特性を、調整する構成としたものがある(例えば、特許文献1参照)。この構成のレーダ装置は、降水によりランダムに反射される反射波が誤目標として検出される頻度をもとに、土砂降り、小降りなどの雨の降り方の違いを検出し、降水の状態にあわせてレーダ装置の制御を行い、誤目標の発生を抑圧する。
また、レーダ装置により降水の状態を求める従来の技術としては、降水により反射されるレーダ電波の受信電力を測定し、受信電力からレーダ方程式により算出される反射因子と降水の強度との関係から降水の強度を算出するものがある(例えば、非特許文献1参照)。
吉田孝 監修、「改訂 レーダ技術」、社団法人電子情報通信学会、平成8年10月1日、P238−253
特許文献1の様に構成されたレーダ信号処理装置やレーダ装置では、降水による誤目標の発生頻度は、降水により反射されるレーダ反射電波の受信電力に影響されるため、降水の状態の他に、レーダアンテナからの距離にも左右される。例えば、降水の状態が同様であっても、レーダアンテナからの距離の近い場所とレーダアンテナからの距離の遠い場所では、誤目標の検出数は異なるものとなる。このため、探査地域がある程度の広さを持ち、探査地域中のそれぞれの地点でレーダアンテナからの距離が大きく異なる場合、反射波が誤目標として検出される頻度のみからは、降水の状態を正しく判定できないという課題があった。また、非特許文献1の様に構成されたレーダ信号処理装置やレーダ装置では、降雨によるレーダ電波の反射波の受信電力を正確に測定するために、広いダイナミックレンジが要求される他、測定値の校正を行なう手段等も必要になる。さらに、降雨によるレーダ電波の反射波の受信電力はレイリー散乱等により分布するが、分布する反射波の受信電力の平均値を正確に測定するため、長時間の観測と、観測データの積分も要する。このため、主として移動する目標物の検出を目的とした、短い周期で目標を観測するレーダ信号処理装置やレーダ装置に、非特許文献1と同様な機能を加えることは、レーダ信号処理装置やレーダ装置の規模の増大や性能の劣化をもたらすため、困難であるという課題があった。
この発明は、上記のような課題を解決するためのものであって、探査地域における降水の強度を算出し、探査地域における目標の確認を容易にするレーダ信号処理装置およびレーダ装置を得るものである。
この発明に係るレーダ信号処理装置は、レーダアンテナが受信したレーダ電波の探査地域からの反射であるレーダ反射電波の、前記探査地域を分割するメッシュ毎の受信電力を、予め決められた電力閾値と比較し、前記レーダ反射電波の受信電力が前記電力閾値を超える前記メッシュである閾値超過メッシュを検出する検出部と、隣接する前記閾値超過メッシュ同士をまとめた前記メッシュのまとまりである検出プロットを作成する検出プロット作成部と、前記検出プロットの面積を予め決められた下限値及び上限値と比較し、前記検出プロットの面積が前記下限値以上かつ前記上限値以下である場合は、前記検出プロットは前記目標物に反射された目標反射電波が検出された目標プロットであり、前記検出プロットの面積が前記下限値未満かまたは前記上限値を越える場合は、前記検出プロットは降水に反射された降水反射電波が検出された降水プロットであると判別するクラッタ判別部と、前記探査地域に設定された、前記レーダアンテナから決められた距離の範囲および決められた方位の範囲の複数の前記メッシュをまとめた区画において、この区画に含まれる前記降水プロットの合計面積がこの区画の面積に対して占める割合を誤検出頻度として算出し、前記降水反射電波についての平均受信電力により決まる受信電力の確率分布関数から、前記降水反射電波の受信電力が前記電力閾値を超過する確率が前記誤検出頻度になる平均受信電力を前記降水反射電波の推定平均受信電力とする反射強度推定部と、前記区画についての前記レーダアンテナとの間の距離と前記推定平均受信電力とレーダ方程式により、降水が前記レーダ電波を反射する強度を表すレーダ反射因子を算出し、前記レーダ反射因子と降水強度との間の予め求められた関係により降水強度を推定する降水強度推定部とを備えるものである。
この発明に係るレーダ装置は、探査地域に向けてレーダ電波を送信し、反射されたレーダ反射電波を受信するレーダアンテナと、前記反射電波の前記探査地域を分割するメッシュ毎の受信電力を、予め決められた電力閾値と比較し、前記レーダ反射電波の受信電力が前記電力閾値を超える前記メッシュである閾値超過メッシュを検出する検出部、隣接する前記閾値超過メッシュ同士をまとめた前記メッシュのまとまりである検出プロットを作成する検出プロット作成部、前記検出プロットの面積を予め決められた下限値及び上限値と比較し、前記検出プロットの面積が前記下限値以上かつ前記上限値以下である場合は、前記検出プロットは前記目標物に反射された目標反射電波が検出された目標プロットであり、前記検出プロットの面積が前記下限値未満かまたは前記上限値を越える場合は、前記検出プロットは降水に反射された降水反射電波が検出された降水プロットであると判別するクラッタ判別部、前記探査地域に設定された、前記レーダアンテナから決められた距離の範囲および決められた方位の範囲の複数の前記メッシュをまとめた区画において、この区画に含まれる前記降水プロットの合計面積がこの区画の面積に対して占める割合を誤検出頻度として算出し、前記降水反射電波についての平均受信電力により決まる受信電力の確率分布関数から、前記降水反射電波の受信電力が前記電力閾値を超過する確率が前記誤検出頻度になる平均受信電力を前記降水反射電波の推定平均受信電力とする反射強度推定部、前記区画についての前記レーダアンテナとの間の距離と前記推定平均受信電力とレーダ方程式により、降水が前記レーダ電波を反射する強度を表すレーダ反射因子を算出し、前記レーダ反射因子と降水強度との間の予め求められた関係により降水強度を推定する降水強度推定部を有するレーダ信号処理装置とを備えるものである。
この発明によれば、探査地域における降水の強度を算出し、探査地域における目標の確認を容易にするレーダ信号処理装置およびレーダ装置を得ることができる。
実施の形態1.
図1は、本発明の実施の形態1に係るレーダ信号処理装置及びレーダ装置の信号処理の概念図である。図1において、探査地域1に存在する車両、航空機、船舶等の目標物を検出するため、レーダ装置は、Aと記した位置に配置したレーダアンテナからレーダ電波を送信し、その反射であるレーダ反射電波を受信し、目標物を検出する。目標物の検出は、探査地域1を距離、方位等で分割するメッシュ2と呼ばれる小区間毎に、レーダアンテナが受信した、それぞれのメッシュ2に該当する場所から反射したレーダ反射電波の受信電力を、予め定めた電力閾値と比較することにより行なう。メッシュ2の大きさは、主として、レーダ装置の空間分解能をもとに定められる。レーダ装置の空間分解能は、レーダアンテナのビーム幅や、レーダ電波の送信パルス長などで決まる、尚、図1では、探査地域1を扇状の形状をした地域として図示したが、Aの地点を中心とした全周の地域としても良い。
図1は、本発明の実施の形態1に係るレーダ信号処理装置及びレーダ装置の信号処理の概念図である。図1において、探査地域1に存在する車両、航空機、船舶等の目標物を検出するため、レーダ装置は、Aと記した位置に配置したレーダアンテナからレーダ電波を送信し、その反射であるレーダ反射電波を受信し、目標物を検出する。目標物の検出は、探査地域1を距離、方位等で分割するメッシュ2と呼ばれる小区間毎に、レーダアンテナが受信した、それぞれのメッシュ2に該当する場所から反射したレーダ反射電波の受信電力を、予め定めた電力閾値と比較することにより行なう。メッシュ2の大きさは、主として、レーダ装置の空間分解能をもとに定められる。レーダ装置の空間分解能は、レーダアンテナのビーム幅や、レーダ電波の送信パルス長などで決まる、尚、図1では、探査地域1を扇状の形状をした地域として図示したが、Aの地点を中心とした全周の地域としても良い。
図1における3は、それぞれ、目標物が検出された、すなわちレーダ反射電波の受信電力が電力閾値を超過したメッシュ2(以後、閾値超過メッシュと呼ぶ)のうちの、隣接するものをまとめたもので、検出プロットと呼ばれるものである。検出された検出プロット3の位置は、レーダ反射電波5Bを反射する反射物の位置を表し、検出プロット3の大きさは、レーダ反射電波5Bを反射する反射物の大きさまたは、反射物の分布の広がりを表す。特に、検出プロット3を、レーダ反射電波5Bを反射した反射物により区別する場合は、検出プロット3の内の、車両、航空機、船舶等の目標物が検出された検出プロット3を目標プロット3A、降水が検出された検出プロットを降水プロット3Bとする。目標プロット3Aは、目標物の大きさに対応する面積を有する。これに対し、降水プロット3Bの面積は、目標プロット3Aに比べて面積が広かったり、逆に狭かったりする。これは、降水などの現象は、目標物とは異なり広範囲に分布し、分布する範囲内でレーダ反射電波5Bを反射する反射の強度も上述した電力閾値より小さい値から大きい値まで変動するためである。
尚、車両、航空機、船舶等の目標物を検出するためのレーダ信号処理装置やレーダ装置においては、検出プロット3Bの様に、本来検出するべき対象でない降水などの現象を検出してしまった場合を、誤検出と呼ぶ。誤検出された降水プロット3Bは、クラッタとも呼ばれる。4は、降水強度を推定するために探査地域内に設定した、レーダアンテナから決められた距離範囲および決められた方位範囲にあるメッシュ2をまとめた区域であり、詳細については、後述する。
図2は、本発明の実施の形態1に係るレーダ装置の構成を示す図である。本発明の実施の形態1に係るレーダ装置は、レーダアンテナ6と、レーダ受信機7と、レーダ信号処理装置8と、レーダ表示器9とを備えている。レーダアンテナ6は、探査地域1にレーダ電波5Aを送信し、探査地域1からのレーダ電波5Aの反射波であるレーダ反射電波5Bを受信する。なお、レーダ反射電波5Bについて、反射物を区別する場合は、以後、目標物からの反射電波を目標反射電波と呼び、降水からの反射電波を降水反射電波と呼ぶ。レーダ受信機7は、レーダアンテナ6が受信したレーダ反射電波5Bを増幅する。レーダ信号処理装置8は、レーダ受信機7により増幅されたレーダ反射波5Bから目標を検出し、検出した目標の位置及び大きさを表す検出プロット3を作成する。レーダ表示器9は、信号処理装置8が作成した検出プロット3を表示画面上に表示する。
また、信号処理装置8は、検出部81と、プロット作成部82と、クラッタ判別部83と、反射強度推定部84と、降水強度推定部85と、パラメータ設定部86とを備えている。検出部81および、プロット作成部82は、レーダ表示器9に表示する検出プロット3を作成する部分であり、探査地域1における目標物を検出するレーダ装置が一般に備える。検出部81は、レーダアンテナ6が受信し、受信機7により増幅されたレーダ反射電波5Bのそれぞれのメッシュ2ごとの受信電力を、予め決められた電力閾値と比較し、レーダ反射電波5Bの受信電力が電力閾値を超過する閾値超過メッシュを検出する。プロット作成部82は、検出部81が検出した閾値超過メッシュの隣接するもの同士のまとまりである検出プロット3を作成する。
クラッタ判別部83、反射強度推定部84、および降水強度推定部85は、プロット作成部82により作成された検出プロット3から降水の強度を推定する部分である。クラッタ判別部83は、それぞれの検出プロット3について、検出プロットの面積(以下、領域サイズと言う)をもとに目標反射電波による目標プロット3Aであるか、降水反射電波による降水プロット3Bであるかの判定を行なう。目標物からの電波反射であれば、検出プロット3の面積は、目標物の面積に対応した一定の下限値から一定の上限値値の間になるため、検出プロット3の面積を、一定の下限値及び一定の上限値と比較し、検出プロットの領域サイズがこの下限値以上かつ上限値以下である場合は、目標反射電波、下限値未満かまたは上限値を超過する場合は降水反射電波と判定するものである。
反射強度推定部84は、クラッタ判別部83で降水からの反射電波と判定された降水プロット3Bの面積の合計である合計面積を基に、降水反射電波の平均受信電力の推定値を算出する。降水強度推定部85は、反射強度推定部84が算出した降水反射電波の平均受信電力の推定値をもとに、降水の強度を推定する。
パラメータ設定部86は、降水強度推定部85が推定する降水強度をもとに、受信機7、検出部81、レーダ表示器9のそれぞれのパラメータの設定を変更する。変更するパラメータは、検出部81が降水反射電波を検出する誤検出確率を調整するためのパラメータとして、検出部81の電力閾値およびレーダ受信機7の増幅率、レーダ表示器9上で降水反射電波を検出した検出プロット3Bの表示のされ方を調整し、視認性を改善するためのパラメータとして、レーダ表示器9の輝度、ブラインドレベルがある。
次に、図2に示した構成のレーダ信号処理装置およびレーダ装置において、降水の強度を推定する原理を説明する。レーダアンテナ6から探査地域1にレーダ電波5Aを送信すると、探査地域1の目標物や、降水により探査地域1の上空の空気中に含まれる水滴により反射されたレーダ反射電波5Bがレーダアンテナ6に受信される。レーダ反射電波5Bが、降水により反射されたものである場合、受信される降水反射電波の受信電力の平均値(Er)は、気象レーダのレーダ方程式により、式(1)のように表される。
ここで、rは探査地域1中のレーダ電波を反射した各部分(メッシュ2)のレーダアンテナ6からの距離、Cはレーダ装置により決まる定数、BRβは、レーダ反射因子と呼ばれる、降水がレーダ電波3Aを反射する強度を示す値であり、Rが降水の強度を表す。このため、受信された降水反射電波の平均受信電力(Er)が求まれば、式(1)を、式(2)のように変形して、降水がレーダ電波3Aを反射する強度を示すレーダ反射因子(BRβ)を算出し、降水の強度を算出することが可能である。
また、降水反射電波は、空気中に分布する多数の水滴からのランダムな反射が重畳されたものである。降水反射電波の受信電力の分布は、レイリー分布で近似が可能である。具体的には、降水反射電波の受信電力が、決められた電力閾値(χ2)を超過する確率分布関数(Pd)は、式(3)で表される。式(3)で表される降水反射電波の受信電力の分布は、降水反射電波の受信電力の平均値である平均受信電力(Er)に対して決まるものである。なお、目標物を検出するレーダ信号処理装置やレーダ装置においては、このように、本来は検出対象でない降雨反射電波が検出されることを、誤検出と呼び、確率(Pd)を誤検出確率と呼ぶ。
上記の(3)式を逆に解いて式(4)のようにすると、誤検出確率(Pd)と電力閾値(χ2)とをもとに、誤検出確率(降雨反射電波の受信電力が決められた電力閾値(χ2)を超過する確率(Pd))が、ある与えられた値になる降水反射電波の平均受信電力(Er)を算出することができる。
このため、降水反射電波の誤検出確率が求まれば、式(4)および式(2)により、降水の強さに対応したレーダ反射因子(BRβ)を算出することができる。なお、レーダ反射因子(BRβ)と降水強度Rとの関係は、非特許文献1等にも記載されている通り、これまでに多くの観測結果等から明らかにされている。具体的には、種々の観測結果よりBおよびβの値は求められているため、レーダ反射因子(BRβ)が求まれば、降水強度Rも算出することができる。
ここで、降雨反射電波の誤検出確率(Pd)は、図1に示した区画4に含まれるメッシュ2で受信される降水反射電波が誤検出された(降水反射電波の受信電力が電力閾値を超過した)誤検出頻度により以下のように推定することができる。
区画4は、降水反射電波の誤検出頻度を算出するために設定されたものであり、レーダアンテナ6から決められた距離の範囲および決められた方位の範囲にあるメッシュ2をまとめたものである。区画4における降水反射電波の誤検出頻度は、後述のように、誤検出確率(Pd)の推定値として使用される。区画4は、探査地域1全体を一定の距離範囲と一定の方位範囲により区切って設定しても良いし、探査地域1の一部の決められた部分のみに、決められた距離の範囲および決められた方位の範囲にあるメッシュ2をまとめて設定しても良い。また、区画4は、探査地域1内で山や建造物などによるグランドクラッタの発生が多い場所を避けて、グランドクラッタの少ないメッシュ2をまとめて設定しても良い。また、区画4は、レーダアンテナ6からの距離が決められた距離範囲内にある、決められた方位範囲のメッシュ2を集めても良いし、レーダアンテナ6からの距離が決められた距離範囲内にある、全ての方位範囲のメッシュを集めても良い。さらに、区画4は複数設定され、複数の区画4で推定された結果を集めて降水の状況が推定されることが望ましいが、区画4が探査地域1内に1つのみ設定され、1つの区画4のみによって推定された結果によって降水の状況が推定されても良い。
区画4内のそれぞれのメッシュ2では、レーダアンテナ6からの距離がほぼ同じであるため、(1)式より、受信電力は、同じ平均受信電力の周りに分布している。このため、区画4内では、同じ電力閾値(χ2)を超過する確率(Pd)もそれぞれのメッシュ2で等しい。この様な場合、区画4内のメッシュ2の数が十分大きければ、誤検出確率(Pd)の推定値は、区画4内の全てのメッシュ2の合計面積に占める、区画4内の降雨反射電波が検出されたメッシュ2の合計面積の割合である誤検出頻度とほぼ一致する。
この様に、本発明の実施の形態1に係るレーダ信号処理装置及びレーダ装置は、降雨反射電波の平均受信電力を直接測定するのではなく、降雨反射電波の誤検出頻度をもとに降雨反射電波の平均受信電力を算出する。このようにすることで、受信電力を精密に測定するための特別な受信機、信号処理装置や、校正のための特別な設備を備える必要が無く、目標物を検出するレーダ装置の通常の受信機、信号処理装置の構成により、降水反射電波の平均受信電力を簡単に求めることができる。また、本発明の実施の形態1に係るレーダ信号処理装置及びレーダ装置は、誤検出確率(Pd)の推定値である誤検出頻度を、区画4内の全てのメッシュ2の合計面積に占める、区画4内の降雨反射電波が検出されたメッシュ2の合計面積の割合として算出する。この様にすることで、多数のメッシュ2で降水反射電波の受信電力と電力閾値との比較を行なった結果を短時間で集めることが可能であり、長時間の観測も不要になる。このため、例えば、レーダアンテナ6を1回スキャンさせた分などの短時間のデータにより降水の強度を推定できる。
図3は、本発明の実施の形態1に係るレーダ信号処理装置における、降水強度の算出に関する処理の流れを表すフローチャートである。図3は、上記の原理に基づき、クラッタ判別部83、反射強度推定部84、および降水強度推定部85が、検出された目標の検出プロットから降水の強度を推定する処理の流れを示す。図3において、クラッタ判定部83に、プロット作成部82が作成した検出プロットが入力される(ST−1)と、クラッタ判定部83は、それぞれの検出プロットに含まれるメッシュ2について、検出プロットの面積(領域サイズ)をもとに、検出されたレーダ反射電波5Bが、目標反射電波であるか、降水反射電波であるかの判定を行なう。判定方法は、検出プロット3の領域サイズが目標プロット3Aについて決められた下限値から上限値の範囲内であれば、目標プロット3A、そうで無ければ降水プロット3Bと判定する(ST−2)。
次に、反射強度推定部84が、それぞれの区画4について、区画4内の降水反射電波が検出された閾値超過メッシュの面積を合計し、合計面積を算出する(ST−3)。更に、反射強度推定部84は、それぞれの区画4について、ST−3で算出した閾値超過メッシュの面積の合計面積がそれぞれの区画4に占める割合を算出する。この算出した割合が誤検出頻度であり、降雨反射電波の受信電力が電力閾値(χ2)を超過する確率(Pd)の推定値である(ST−4)。降水強度推定部84は、算出した誤検出頻度を、式(4)の誤検出確率(Pd)に当てはめて、それぞれの区画4について、降水反射電波の受信電力が電力閾値を超過する確率(Pd)が算出した誤検出頻度になる、平均受信電力(Er)の推定値を算出する(ST−5)。
降雨強度推定部85は、それぞれの区画4について、反射強度推定部84が算出した、の平均受信電力(Er)の推定値をもとに、式(2)により、降水がレーダ電波5Aを反射する強度を表すレーダ反射因子(BRβ)を算出する(ST−6)。次に、降雨強度推定部85は、複数の区画4について、レーダ反射因子(BRβ)の特異な値を除いた上で、平均値を算出する(ST−7)。レーダ反射因子(BRβ)は、レーダアンテナ6からの距離によらない値であるため、降雨強度推定部85は、レーダアンテナ6からの距離が同じか異なるかに関わらず、それぞれの区画4について算出したレーダ反射因子(BRβ)を平均することにより、降水以外の原因によるレーダ反射電波の影響や、観測のばらつきにより生じる誤差を低減する。更に、Bおよびβについての一般に知られている値や、実測をもとに解析した値をもとに、レーダ反射因子(BRβ)の平均値より、降水強度(R)の値を算出する(ST−8)。
図4は、本発明の実施の形態1に係るレーダ信号処理装置8における、降水強度の推定に基づくパラメータ調整に関する処理の流れを表すフローチャートである。このフローチャートにおける降雨強度の算出(ST−14)は、図3に示した、クラッタ判別部83、反射強度推定部84、および降水強度推定部85が、検出された目標のプロットから降水の強度を推定する処理の流れ(ST−1〜ST−8)を表す。パラメータ設定部86は、降水強度推定部85が推定する降水強度をもとに、受信機7、検出部81、レーダ表示器9のそれぞれのパラメータの設定を変更する。これらのパラメータの内の、特に受信機7や検出部81のパラメータを変更する場合、検出部81において閾値超過メッシュを検出する際の検出条件を変更することになる。その結果、降雨強度の推定結果が不安定になる可能性がある。このため、図4に示すように、降雨強度の推定は一定時間毎に、受信機7、検出部81のパラメータを同じ条件に初期化して実施する。
図4において、レーダ信号処理装置8は、起動されると、受信機7、検出部81、レーダ表示器9のパラメータを初期化する(ST−11)。初期化されるパラメータは、検出部81の電力閾値およびレーダ受信機7の増幅率、レーダ表示器9の輝度、ブラインドレベル等である。パラメータが初期化されると、検出部81は、受信機7から入力されるレーダ反射電波5Bの受信電力を、ぞれぞれのメッシュ2毎に予め決められた電力閾値と比較し、レーダ反射電波の受信電力が電力閾値を超過する閾値超過メッシュを検出する(ST−12)。プロット作成部82は、検出部81が検出した閾値超過メッシュの隣接するもの同士のまとまりである検出プロット3A、3Bを作成する(ST−13)。
プロット作成部により検出プロットが作成されると、クラッタ判定部83、反射強度推定部84、および降水強度推定部85は、降水強度の推定値を算出する(ST−14)。降水強度の推定値を算出する処理の詳細は、図3のST−1〜ST−8の通りである。
降水強度の推定値が算出されると、パラメータ設定部86は、降水強度の推定値に従い、誤目標の発生を抑えるよう、受信機7、検出部81、レーダ表示器9のパラメータを設定する(ST−15)。パラメータ設定部86は、例えば、降水強度の値に対して適した受信機7、検出部81、レーダ表示器9のパラメータを実験的に求めて記憶しておき、記憶されているパラメータから、降水強度の推定値の値に該当するパラメータを選択し、設定する。パラメータの設定を行なうと、パラメータ設定部86は、次に降水強度の推定値の算出、パラメータの設定を行なうまでの時間を設定したタイマを開始させる(ST−16)。以後、タイマに設定した時間まで経過するまでの間(ST−17でNoである間)は、レーダ信号処理装置8は、通常の目標検出を運用しており、パラメータ設定部86が設定したパラメータにより、レーダ反射電波から目標の検出を行なう(ST−18)。ST−18で実施する目標の検出は、検出部の閾値超過メッシュの検出(ST−12)、プロット作成部82による検出プロットの作成(ST−13)と、処理に使用されるパラメータの設定が異なる以外は、全く同様である。タイマに設定した時間が経過(ST−17でYes)すると、レーダ信号処理装置8は、受信機7、検出部81、レーダ表示器9のパラメータを初期化して降水強度の推定を実施する(ST−11〜ST−15)。
図4において、ST−11〜ST−15は、降水強度を推定し、降水強度の推定値に従い、誤目標の発生を抑えるよう、受信機7、検出部81、レーダ表示器9のパラメータを設定する処理である。また、(ST−18)は、レーダ信号処理装置8が、目標物の検出を行なう処理である。図4の様に、レーダ信号処理装置8は、通常の目標検出処理(ST−18)の間に、定期的に降水強度の推定とパラメータの変更(ST−11〜ST−15)を実施するようにすることで、常に、気象の変化に合わせて、誤目標の発生を抑えることができる。
また、降水強度を推定するにあたり、降水反射電波の受信電力を精密に測定するのではなく、降水反射電波の受信電力が検出された検出プロットの面積により降水の強度を推定するため、長時間の測定や積分も不要であり、また、精密な測定のための設備を要することなく、従来のレーダ装置の有する構成そのままで実施することができる。
1 探査地域
2 メッシュ
3 検出プロット
3A 目標プロット
3B 降水プロット
4 区域
5A レーダ電波
5B レーダ反射電波
6 レーダアンテナ
7 レーダ受信機
8 レーダ信号処理装置
9 レーダ表示器
81 検出部
82 プロット作成部
83 クラッタ判別部
84 反射強度推定部
85 降水強度推定部
86 パラメータ設定部
2 メッシュ
3 検出プロット
3A 目標プロット
3B 降水プロット
4 区域
5A レーダ電波
5B レーダ反射電波
6 レーダアンテナ
7 レーダ受信機
8 レーダ信号処理装置
9 レーダ表示器
81 検出部
82 プロット作成部
83 クラッタ判別部
84 反射強度推定部
85 降水強度推定部
86 パラメータ設定部
Claims (4)
- レーダアンテナが受信したレーダ電波の探査地域からの反射であるレーダ反射電波の、前記探査地域を分割するメッシュ毎の受信電力を、予め決められた電力閾値と比較し、前記レーダ反射電波の受信電力が前記電力閾値を超える前記メッシュである閾値超過メッシュを検出する検出部と、
隣接する前記閾値超過メッシュ同士をまとめた前記メッシュのまとまりである検出プロットを作成する検出プロット作成部と、
前記検出プロットの面積を予め決められた下限値及び上限値と比較し、前記検出プロットの面積が前記下限値以上かつ前記上限値以下である場合は、前記検出プロットは前記目標物に反射された目標反射電波が検出された目標プロットであり、前記検出プロットの面積が前記下限値未満かまたは前記上限値を越える場合は、前記検出プロットは降水に反射された降水反射電波が検出された降水プロットであると判別するクラッタ判別部と、
前記探査地域に設定された、前記レーダアンテナから決められた距離の範囲および決められた方位の範囲の複数の前記メッシュをまとめた区画において、この区画に含まれる前記降水プロットの合計面積がこの区画の面積に対して占める割合を誤検出頻度として算出し、前記降水反射電波についての平均受信電力により決まる受信電力の確率分布関数から、前記降水反射電波の受信電力が前記電力閾値を超過する確率が前記誤検出頻度になる平均受信電力を前記降水反射電波の推定平均受信電力とする反射強度推定部と、
前記区画についての前記レーダアンテナとの間の距離と前記推定平均受信電力とレーダ方程式により、降水が前記レーダ電波を反射する強度を表すレーダ反射因子を算出し、前記レーダ反射因子と降水強度との間の予め求められた関係により降水強度を推定する降水強度推定部と
を備えたレーダ信号処理装置。 - 前記反射強度推定部は、複数の前記区画においてそれぞれ推定平均受信電力を算出し、
前記降水強度推定部は、それぞれの前記区画についての前記レーダアンテナとの間の距離と前記推定平均受信電力とレーダ方程式により算出した値を平均化した前記レーダ反射因子を算出し、前記レーダ反射因子と降水強度との間の予め求められた関係により降水強度を推定する
請求項1に記載のレーダ信号処理装置。 - 前記降水強度推定部が推定する降水強度に基づき、前記電力閾値を変更するパラメータ設定部を更に備えた
請求項1または2に記載のレーダ信号処理装置。 - 探査地域に向けてレーダ電波を送信し、反射されたレーダ反射電波を受信するレーダアンテナと、
前記反射電波の前記探査地域を分割するメッシュ毎の受信電力を、予め決められた電力閾値と比較し、前記レーダ反射電波の受信電力が前記電力閾値を超える前記メッシュである閾値超過メッシュを検出する検出部、
隣接する前記閾値超過メッシュ同士をまとめた前記メッシュのまとまりである検出プロットを作成する検出プロット作成部、
前記検出プロットの面積を予め決められた下限値及び上限値と比較し、前記検出プロットの面積が前記下限値以上かつ前記上限値以下である場合は、前記検出プロットは前記目標物に反射された目標反射電波が検出された目標プロットであり、前記検出プロットの面積が前記下限値未満かまたは前記上限値を越える場合は、前記検出プロットは降水に反射された降水反射電波が検出された降水プロットであると判別するクラッタ判別部、
前記探査地域に設定された、前記レーダアンテナから決められた距離の範囲および決められた方位の範囲の複数の前記メッシュをまとめた区画において、この区画に含まれる前記降水プロットの合計面積がこの区画の面積に対して占める割合を誤検出頻度として算出し、前記降水反射電波についての平均受信電力により決まる受信電力の確率分布関数から、前記降水反射電波の受信電力が前記電力閾値を超過する確率が前記誤検出頻度になる平均受信電力を前記降水反射電波の推定平均受信電力とする反射強度推定部、
前記区画についての前記レーダアンテナとの間の距離と前記推定平均受信電力とレーダ方程式により、降水が前記レーダ電波を反射する強度を表すレーダ反射因子を算出し、前記レーダ反射因子と降水強度との間の予め求められた関係により降水強度を推定する降水強度推定部を有するレーダ信号処理装置と
を備えたレーダ装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2015031334A JP2016153728A (ja) | 2015-02-20 | 2015-02-20 | レーダ信号処理装置及びレーダ装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2015031334A JP2016153728A (ja) | 2015-02-20 | 2015-02-20 | レーダ信号処理装置及びレーダ装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2016153728A true JP2016153728A (ja) | 2016-08-25 |
Family
ID=56761176
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2015031334A Pending JP2016153728A (ja) | 2015-02-20 | 2015-02-20 | レーダ信号処理装置及びレーダ装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2016153728A (ja) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101770742B1 (ko) | 2016-12-28 | 2017-08-23 | 엘아이지넥스원 주식회사 | 클러터를 억제하는 표적 탐지 장치 및 그 방법 |
CN111830481A (zh) * | 2020-07-09 | 2020-10-27 | 中国人民解放军海军航空大学 | 雷达回波单分量幅度分布模型参数估计方法及装置 |
CN112255614A (zh) * | 2020-12-23 | 2021-01-22 | 北京海兰信数据科技股份有限公司 | 一种手动抑制导航雷达海杂波的方法及系统 |
CN113923594A (zh) * | 2021-10-12 | 2022-01-11 | 中国人民解放军国防科技大学 | 一种基于时分长期演进网络的晴雨区分方法 |
CN115267787A (zh) * | 2022-08-17 | 2022-11-01 | 北京城市气象研究院 | 一种天气雷达定量降雪的估测方法 |
-
2015
- 2015-02-20 JP JP2015031334A patent/JP2016153728A/ja active Pending
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101770742B1 (ko) | 2016-12-28 | 2017-08-23 | 엘아이지넥스원 주식회사 | 클러터를 억제하는 표적 탐지 장치 및 그 방법 |
CN111830481A (zh) * | 2020-07-09 | 2020-10-27 | 中国人民解放军海军航空大学 | 雷达回波单分量幅度分布模型参数估计方法及装置 |
CN112255614A (zh) * | 2020-12-23 | 2021-01-22 | 北京海兰信数据科技股份有限公司 | 一种手动抑制导航雷达海杂波的方法及系统 |
CN113923594A (zh) * | 2021-10-12 | 2022-01-11 | 中国人民解放军国防科技大学 | 一种基于时分长期演进网络的晴雨区分方法 |
CN113923594B (zh) * | 2021-10-12 | 2023-10-03 | 中国人民解放军国防科技大学 | 一种基于时分长期演进网络的晴雨区分方法 |
CN115267787A (zh) * | 2022-08-17 | 2022-11-01 | 北京城市气象研究院 | 一种天气雷达定量降雪的估测方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP2016153728A (ja) | レーダ信号処理装置及びレーダ装置 | |
US9348015B2 (en) | Integrated rainfall estimation method using X-band dual-polarimetric radar measurement data | |
Upton et al. | Microwave links: The future for urban rainfall measurement? | |
Nuimura et al. | Elevation changes of glaciers revealed by multitemporal digital elevation models calibrated by GPS survey in the Khumbu region, Nepal Himalaya, 1992-2008 | |
CN105393135B (zh) | 机动车的雷达传感器的俯仰失调角的确定 | |
US7132974B1 (en) | Methods and systems for estimating three dimensional distribution of turbulence intensity using radar measurements | |
US7515087B1 (en) | Weather radar system and method using data from a lightning sensor | |
US9459348B2 (en) | Millimeter wave radar system for and method of weather detection | |
US9097805B2 (en) | Systems and methods for calibrating dual polarization radar systems | |
JP4667426B2 (ja) | 気象レーダ装置 | |
KR100922130B1 (ko) | 도플러 기상레이더의 이차 에코 제거 방법 | |
Goudenhoofdt et al. | Generation and verification of rainfall estimates from 10-yr volumetric weather radar measurements | |
JP5582732B2 (ja) | 降雪検知装置および降雪検知方法 | |
Vulpiani et al. | Characterization of Mediterranean hail-bearing storms using an operational polarimetric X-band radar | |
KR101538368B1 (ko) | 엑스밴드 이중편파 레이더 원시자료를 이용한 강수체 구분 방법 | |
US20160161608A1 (en) | Micro climate corrections for radar interferometry measurements | |
Bellon et al. | Real-time comparisons of VPR-corrected daily rainfall estimates with a gauge mesonet | |
RU2549192C1 (ru) | Способ распознавания цели (варианты) | |
Einfalt et al. | Quality control of precipitation data | |
JP7110004B2 (ja) | 処理装置、処理方法、およびプログラム | |
JP5429493B2 (ja) | 降雨観測設備の投資計画評価方法 | |
JP2020076665A (ja) | 積雪特性を測定する方法及びその装置及びこの積雪特性を測定する方法を利用した融雪災害の予測監視方法及びその装置 | |
KR101840651B1 (ko) | 구름레이더를 이용한 강우 강도 산출 시스템 및 이를 이용한 강우 강도 산출 방법 | |
Cyr | Estimation of ZR relationship and comparative analysis of precipitation data from colocated rain-gauge, vertical radar and disdrometer | |
McRoberts | Minimizing Biases in Radar Precipitation Estimates |