JP2011214972A - 気象レーダ信号処理装置及びそのグランドクラッタ除去方法 - Google Patents
気象レーダ信号処理装置及びそのグランドクラッタ除去方法 Download PDFInfo
- Publication number
- JP2011214972A JP2011214972A JP2010082879A JP2010082879A JP2011214972A JP 2011214972 A JP2011214972 A JP 2011214972A JP 2010082879 A JP2010082879 A JP 2010082879A JP 2010082879 A JP2010082879 A JP 2010082879A JP 2011214972 A JP2011214972 A JP 2011214972A
- Authority
- JP
- Japan
- Prior art keywords
- weather echo
- echo
- index
- weather
- data
- 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.)
- Granted
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
【解決手段】MTI処理A1により、目標反射波受信信号のIQデータからグランドクラッタ成分を除去することで気象エコー成分を抽出する。続いて、スペクトル評価処理A2により、気象エコー成分から順次正規化したパルスペア合成ベクトルを求め、そのパルスペア合成ベクトルの総和の絶対値からスペクトルの平坦さを評価する評価指標を算出し、この評価指標に基づいてグランドクラッタの除去と気象エコーの保持を行う。最終的に、ノイズ除去処理A3により、処理A2の評価指標に対する振幅成分の電力値に対する閾値処理を行って孤立点を除去することでノイズを除去する。このようにしてスペクトルの平坦さを評価し、この評価指標に対する閾値処理によりグランドクラッタ除去を行う。
【選択図】 図2
Description
図3は、上記スペクトル評価処理(ステップA2)の第1の実施例を示すブロック図である。すなわち、この実施例のスペクトル評価処理では、位相変換部S1、パルスペア合成ベクトル算出部(和)S2、グランドクラッタ除去部S3で構成される。
IQA(n)=r(n)ejθ(n) (n=0,1,…,N-1) …(1)
とすると、位相変換後のIQデータは、IQデータから算出した平均速度をv、ナイキスト速度をVnyqとすると
IQB(n)=r(n)ej(θ(n)-πv(n-1)/Vnyq) …(2)
となる。
次に、スペクトルの平坦さを評価する指標Jを次式より算出する(ステップS23)。
J=|ΣIQC(n)| (nは0〜N-2) …(4)
続いて、グランドクラッタ除去判別処理を行い、上記指標Jが予め決められた閾値未満か否かを判断する(ステップS24)。上記指標Jが閾値未満の場合には(yes)、スペクトルが平坦であってそのメッシュには気象エコーが存在しないと判断してその時系列データを除去する(ステップS25)。ステップS24において、指標Jが閾値以上の場合には(no)、そのメッシュには気象エコーが存在すると判断してそのまま時系列データを出力する。
図10は、図2に示したスペクトル評価処理S2において、第1の実施例により簡易な手法で処理する場合の第2の実施例の構成を示すブロック図である。すなわち、この実施例のスペクトル評価処理では、パルスペア合成ベクトル算出部(積)S4及びグランドクラッタ除去部S5で構成される。
z(n)=IQA *(n)・IQA(n+1)
=|IQA(n)||IQA(n+1)|ej[arg{IQA(n+1)}-arg{IQA(n)}] …(5)
*:共役複素数
次に、スペクトルの平坦さを評価する指標Zを次式から求める(ステップS43)。
Z=|Σz(n)| (nは0〜N-2) …(6)
続いて、グランドクラッタ除去判別処理を行い、上記指標Zが予め決められた閾値未満か否かを判断する(ステップ44)。上記指標Zが閾値未満の場合には(yes)、スペクトルが平坦であってそのメッシュには気象エコーが存在しないと判断してその時系列データを除去する(ステップS45)。ステップS44において、指標Zが閾値以上の場合には(no)、そのメッシュには気象エコーが存在すると判断してそのまま時系列データを出力する。
図16は、従来の速度幅フィルタによる処理例と上記第1及び第2の実施例それぞれの処理例に関するシミュレーション結果を示す特性図である。図において、実線部分はMTI処理後のグランドクラッタから算出した特性、点線部分はMTI処理後の速度幅5m/sの気象エコーから算出した特性を示しており、図16(a)は従来の速度幅フィルタを用いた場合、図16(b)は第1の実施例の場合、図16(c)は第2の実施例の場合を示している。
Claims (10)
- 一定周期で繰り返し得られるレーダパルス反射波の受信信号のIQデータから地形エコーによるグランドクラッタ成分を演算し除去することで気象エコー成分のIQデータを抽出するMTI処理手段と、
前記MTI処理手段で繰り返し得られる気象エコー成分の時系列IQデータから相互間のベクトル成分を求め、そのベクトル成分の総和の絶対値からスペクトルの平坦さを評価する評価指標を求めるスペクトル評価処理手段と、
前記スペクトル評価処理手段で得られた評価指標に基づいて前記MTI処理手段で得られた気象エコー成分のIQデータから前記グランドクラッタ成分の消え残りを判別し削除するグランドクラッタ除去手段と
を具備することを特徴とする気象レーダ信号処理装置。 - 前記スペクトル評価処理手段は、
前記MTI処理手段で繰り返し得られる気象エコー成分IQデータの速度を0に位相変換して正規化する正規化手段と、
前記正規化された気象エコー成分IQデータを順次入力して時系列データの和で処理する場合の正規化パルスペア合成ベクトルを求め、この正規化パルスペア合成ベクトルからスペクトルの平坦さを評価する指標を算出する評価指標算出手段と、
前記指標を予め決められた閾値と比較して閾値未満であるときは、スペクトルが平坦であってそのメッシュには気象エコーが存在しないと判断してその時系列データを除去し、閾値以上であるときは、そのメッシュには気象エコーが存在すると判断してそのまま時系列データを出力することを特徴とする請求項1記載の気象レーダ信号処理装置。 - 前記スペクトル評価処理手段は、
前記MTI処理手段で得られた気象エコー成分IQデータを順次入力して時系列データの積で処理する場合のパルスペア自己相関ベクトルを求める自己相関ベクトル算出手段と、
前記パルスペア自己相関ベクトルからスペクトルの平坦さを評価する指標を算出する評価指標算出手段と、
前記指標を予め決められた閾値と比較して閾値未満であるときは、スペクトルが平坦であってそのメッシュには気象エコーが存在しないと判断してその時系列データを除去し、閾値以上であるときは、そのメッシュには気象エコーが存在すると判断してそのまま時系列データを出力することを特徴とする請求項1記載の気象レーダ信号処理装置。 - 前記グランドクラッタ除去手段は、速度幅一定、降水量一定を条件とする模擬気象エコーを生成し、この模擬気象エコーをフーリエ変換して周波数領域に変換し、この周波数領域の模擬気象エコーに対して、平均値と標準偏差が速度、速度幅に対応するガウス関数を乗じて、逆フーリエ変換を施して時間領域の模擬気象エコーを得る模擬気象エコー生成手段と、
前記模擬気象エコー生成手段で生成された模擬気象エコーにMTI処理を行い、正規化パルスペア合成ベクトルの算出を行い、スペクトルの平坦さを評価する指標を算出する模擬気象エコー評価指標算出手段と、
前記模擬気象エコー生成手段で生成される複数回の模擬気象エコーについて前記模擬気象エコー評価指標算出手段によりそれぞれの評価指標を算出し、保持したい気象エコーの指標の分布を求め、それを保持できるように分布の下限値に閾値を設定する閾値設定手段とを備え、
前記閾値設定手段で設定される閾値に基づいて前記MTI処理手段で得られた気象エコー成分IQデータからグランドクラッタ成分を除去することを特徴とする請求項2記載の気象レーダ信号処理装置。 - 前記グランドクラッタ除去手段は、速度幅一定、降水量一定を条件とする模擬気象エコーを生成し、この模擬気象エコーをフーリエ変換して周波数領域に変換し、この周波数領域の模擬気象エコーに対して、平均値と標準偏差が速度、速度幅に対応するガウス関数を乗じて、逆フーリエ変換を施して時間領域の模擬気象エコーを得る模擬気象エコー生成手段と、
前記模擬気象エコー生成手段で生成された模擬気象エコーにMTI処理を行い、パルスペア自己相関ベクトルの算出を行い、スペクトルの平坦さを評価する指標を算出する模擬気象エコー評価指標算出手段と、
前記模擬気象エコー生成手段で生成される複数回の模擬気象エコーについて前記模擬気象エコー評価指標算出手段によりそれぞれの評価指標を算出し、保持したい気象エコーの指標の分布を求め、それを保持できるように分布の下限値に閾値を設定する閾値設定手段とを備え、
前記閾値設定手段で設定される閾値に基づいて前記MTI処理手段で得られた気象エコー成分IQデータからグランドクラッタ成分を除去することを特徴とする請求項3記載の気象レーダ信号処理装置。 - 一定周期で繰り返し得られるレーダパルス反射波の受信信号のIQデータから地形エコーによるグランドクラッタ成分を演算し除去することで気象エコー成分のIQデータを抽出し、
前記繰り返し得られる気象エコー成分の時系列IQデータから相互間のベクトル成分を求め、そのベクトル成分の総和の絶対値からスペクトルの平坦さを評価する評価指標を求め、
前記評価指標に基づいて前記気象エコー成分のIQデータから前記グランドクラッタ成分の消え残りを判別し削除することを特徴とする気象レーダ信号処理装置のグランドクラッタ除去方法。 - 前記スペクトルの平坦さの評価処理は、
前記繰り返し得られる気象エコー成分IQデータの速度を0に位相変換して正規化し、
前記正規化された気象エコー成分IQデータを順次入力して時系列データの和で処理する場合の正規化パルスペア合成ベクトルを求め、
この正規化パルスペア合成ベクトルからスペクトルの平坦さを評価する指標を算出し、
前記指標を予め決められた閾値と比較して閾値未満であるときは、スペクトルが平坦であってそのメッシュには気象エコーが存在しないと判断してその時系列データを除去し、閾値以上であるときは、そのメッシュには気象エコーが存在すると判断してそのまま時系列データを出力することを特徴とする請求項6記載の気象レーダ信号処理装置のグランドクラッタ除去方法。 - 前記スペクトルの平坦さの評価処理は、
前記気象エコー成分IQデータを順次入力して時系列データの積で処理する場合のパルスペア自己相関ベクトルを求め、
前記パルスペア自己相関ベクトルからスペクトルの平坦さを評価する指標を算出し、
前記指標を予め決められた閾値と比較して閾値未満であるときは、スペクトルが平坦であってそのメッシュには気象エコーが存在しないと判断してその時系列データを除去し、閾値以上であるときは、そのメッシュには気象エコーが存在すると判断してそのまま時系列データを出力することを特徴とする請求項6記載の気象レーダ信号処理装置のグランドクラッタ除去方法。 - 前記グランドクラッタの除去処理は、
速度幅一定、降水量一定を条件とする模擬気象エコーを生成し、この模擬気象エコーをフーリエ変換して周波数領域に変換し、この周波数領域の模擬気象エコーに対して、平均値と標準偏差が速度、速度幅に対応するガウス関数を乗じて、逆フーリエ変換を施して時間領域の模擬気象エコーを複数回取得し、
前記模擬気象エコーにMTI処理を行い、正規化パルスペア合成ベクトルの算出を行い、スペクトルの平坦さを評価する指標を算出し、
前記複数回の模擬気象エコーについてそれぞれの評価指標を算出し、保持したい気象エコーの指標の分布を求め、それを保持できるように分布の下限値に閾値を設定し、
前記閾値に基づいて前記気象エコー成分IQデータからグランドクラッタ成分を除去することを特徴とする請求項7記載の気象レーダ信号処理装置のグランドクラッタ除去方法。 - 前記グランドクラッタの除去処理は、
速度幅一定、降水量一定を条件とする模擬気象エコーを生成し、この模擬気象エコーをフーリエ変換して周波数領域に変換し、この周波数領域の模擬気象エコーに対して、平均値と標準偏差が速度、速度幅に対応するガウス関数を乗じて、逆フーリエ変換を施して時間領域の模擬気象エコーを複数回取得し、
前記模擬気象エコーにMTI処理を行い、パルスペア自己相関ベクトルの算出を行い、スペクトルの平坦さを評価する指標を算出し、
前記複数回の模擬気象エコーについてそれぞれの評価指標を算出し、保持したい気象エコーの指標の分布を求め、それを保持できるように分布の下限値に閾値を設定し、
前記閾値に基づいて前記気象エコー成分IQデータからグランドクラッタ成分を除去することを特徴とする請求項8記載の気象レーダ信号処理装置のグランドクラッタ除去方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010082879A JP5701511B2 (ja) | 2010-03-31 | 2010-03-31 | 気象レーダ信号処理装置及びそのグランドクラッタ除去方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010082879A JP5701511B2 (ja) | 2010-03-31 | 2010-03-31 | 気象レーダ信号処理装置及びそのグランドクラッタ除去方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2011214972A true JP2011214972A (ja) | 2011-10-27 |
JP5701511B2 JP5701511B2 (ja) | 2015-04-15 |
Family
ID=44944871
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2010082879A Active JP5701511B2 (ja) | 2010-03-31 | 2010-03-31 | 気象レーダ信号処理装置及びそのグランドクラッタ除去方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP5701511B2 (ja) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2015068700A (ja) * | 2013-09-27 | 2015-04-13 | パナソニック株式会社 | レーダ装置及び物体検出方法 |
CN109298424A (zh) * | 2018-11-28 | 2019-02-01 | 广东电网有限责任公司 | 一种基于雷达拼图的杆塔电力线路天气监测方法和装置 |
CN111624556A (zh) * | 2020-06-08 | 2020-09-04 | 河海大学 | 基于形态成分分析的气象雷达wtc抑制方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS6150283U (ja) * | 1984-09-07 | 1986-04-04 | ||
JPS62129774A (ja) * | 1985-11-30 | 1987-06-12 | Nec Corp | ウエザクラツタ抑圧mti回路 |
JPS62231185A (ja) * | 1986-03-31 | 1987-10-09 | Nec Corp | 擬似移動目標発生装置 |
JPH032433B2 (ja) * | 1984-02-14 | 1991-01-16 | Nippon Electric Co | |
JPH06214016A (ja) * | 1993-01-18 | 1994-08-05 | Mitsubishi Electric Corp | ドップラレーダの信号処理装置 |
JPH09257910A (ja) * | 1996-03-21 | 1997-10-03 | Mitsubishi Electric Corp | レーダ信号処理装置 |
JP2001091633A (ja) * | 1999-09-20 | 2001-04-06 | Toshiba Corp | 模擬クラッタ発生装置 |
JP2006226954A (ja) * | 2005-02-21 | 2006-08-31 | Toshiba Corp | レーダ装置とその信号処理方法 |
JP2006234580A (ja) * | 2005-02-24 | 2006-09-07 | Toshiba Corp | 気象レーダシステム及びそれに用いられる信号処理方法 |
JP2007298406A (ja) * | 2006-04-28 | 2007-11-15 | Toshiba Corp | レーダ装置とその信号処理方法 |
-
2010
- 2010-03-31 JP JP2010082879A patent/JP5701511B2/ja active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH032433B2 (ja) * | 1984-02-14 | 1991-01-16 | Nippon Electric Co | |
JPS6150283U (ja) * | 1984-09-07 | 1986-04-04 | ||
JPS62129774A (ja) * | 1985-11-30 | 1987-06-12 | Nec Corp | ウエザクラツタ抑圧mti回路 |
JPS62231185A (ja) * | 1986-03-31 | 1987-10-09 | Nec Corp | 擬似移動目標発生装置 |
JPH06214016A (ja) * | 1993-01-18 | 1994-08-05 | Mitsubishi Electric Corp | ドップラレーダの信号処理装置 |
JPH09257910A (ja) * | 1996-03-21 | 1997-10-03 | Mitsubishi Electric Corp | レーダ信号処理装置 |
JP2001091633A (ja) * | 1999-09-20 | 2001-04-06 | Toshiba Corp | 模擬クラッタ発生装置 |
JP2006226954A (ja) * | 2005-02-21 | 2006-08-31 | Toshiba Corp | レーダ装置とその信号処理方法 |
JP2006234580A (ja) * | 2005-02-24 | 2006-09-07 | Toshiba Corp | 気象レーダシステム及びそれに用いられる信号処理方法 |
JP2007298406A (ja) * | 2006-04-28 | 2007-11-15 | Toshiba Corp | レーダ装置とその信号処理方法 |
Non-Patent Citations (2)
Title |
---|
CSNC201100306010; 吉田孝: '第9章 気象レーダ' 改訂 レーダ技術 , 19961001, 238-253頁 * |
JPN6013052063; 吉田孝: '第9章 気象レーダ' 改訂 レーダ技術 , 19961001, 238-253頁 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2015068700A (ja) * | 2013-09-27 | 2015-04-13 | パナソニック株式会社 | レーダ装置及び物体検出方法 |
US10126417B2 (en) | 2013-09-27 | 2018-11-13 | Panasonic Corporation | Radar apparatus and object detecting method |
CN109298424A (zh) * | 2018-11-28 | 2019-02-01 | 广东电网有限责任公司 | 一种基于雷达拼图的杆塔电力线路天气监测方法和装置 |
CN111624556A (zh) * | 2020-06-08 | 2020-09-04 | 河海大学 | 基于形态成分分析的气象雷达wtc抑制方法 |
Also Published As
Publication number | Publication date |
---|---|
JP5701511B2 (ja) | 2015-04-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20120179047A1 (en) | Method and apparatus for automatic optimization of doppler imaging parameters | |
JP6415288B2 (ja) | レーダ装置 | |
KR101953697B1 (ko) | 균일 선형 배열 안테나 기반의 가상 수신신호 생성을 이용한 도래각 추정 방법 및 장치 | |
WO2022061597A1 (zh) | 一种信号噪声滤除方法、装置、存储介质及激光雷达 | |
JP5699404B2 (ja) | レーダ受信信号処理装置とその方法 | |
US10061014B2 (en) | Radar signal processing method and apparatus for compensating for in-phase/quadrature-phase (I/Q) imbalance | |
JP2008014837A (ja) | レーダ装置とその信号処理方法 | |
JP5752481B2 (ja) | パッシブレーダ装置、誘導装置及び電波検出方法 | |
JP6279193B2 (ja) | 物体検出装置及びセンサ装置 | |
JP5701511B2 (ja) | 気象レーダ信号処理装置及びそのグランドクラッタ除去方法 | |
EP2990820B1 (en) | Clutter suppressing device and radar apparatus provided with the same | |
JP5597411B2 (ja) | 気象レーダ信号処理装置及びそのグランドクラッタ除去方法 | |
JP2011237338A (ja) | レーダ装置 | |
KR101041990B1 (ko) | 레이더 모의표적장치의 도플러 주파수 구현방법 | |
JP5766896B1 (ja) | ピーク周波数検出装置、方法およびプログラム | |
Vizinho et al. | Evaluation of the use of the modified-covariance method in HF radar ocean measurement | |
JP3881078B2 (ja) | 周波数推定方法、周波数推定装置、ドップラソナーおよび潮流計 | |
JP6398689B2 (ja) | 信号処理装置、信号処理方法及びプログラム | |
EP3835822B1 (en) | Signal processing device and signal processing method | |
JP2012149951A (ja) | レーダ装置 | |
KR101282692B1 (ko) | 충격음의 음장 표시 방법 | |
JP2005017143A (ja) | 気象レーダ信号処理装置 | |
Lu et al. | Impulsive noise excision and performance analysis | |
Sinitsyn | Copula based detection algorithm for MIMO ultrawideband noise radars | |
KR101390679B1 (ko) | 창 함수 비선형 에너지 연산자를 이용한 투과 초음파 단층 촬영기에서 최초 도달 펄스를 검출하기 위한 방법 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20121018 |
|
RD02 | Notification of acceptance of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7422 Effective date: 20130731 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20131009 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20131022 |
|
RD02 | Notification of acceptance of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7422 Effective date: 20131024 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20131224 |
|
RD07 | Notification of extinguishment of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7427 Effective date: 20140319 |
|
RD07 | Notification of extinguishment of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7427 Effective date: 20140508 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20140715 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20140807 |
|
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: 20150120 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20150218 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 5701511 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
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 |