CN104614722A - 一种基于信噪比识别雷达遮挡的方法 - Google Patents
一种基于信噪比识别雷达遮挡的方法 Download PDFInfo
- Publication number
- CN104614722A CN104614722A CN201510036478.4A CN201510036478A CN104614722A CN 104614722 A CN104614722 A CN 104614722A CN 201510036478 A CN201510036478 A CN 201510036478A CN 104614722 A CN104614722 A CN 104614722A
- Authority
- CN
- China
- Prior art keywords
- radar
- signal
- data
- noise ratio
- threshold value
- 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
- 238000000034 method Methods 0.000 title claims abstract description 23
- 230000000052 comparative effect Effects 0.000 claims description 2
- 238000013461 design Methods 0.000 claims description 2
- 238000010408 sweeping Methods 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 abstract description 2
- 238000012545 processing Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 abstract 2
- 230000005540 biological transmission Effects 0.000 abstract 1
- 230000005855 radiation Effects 0.000 description 9
- 230000005457 Black-body radiation Effects 0.000 description 8
- VMXUWOKSQNHOCA-UKTHLTGXSA-N ranitidine Chemical compound [O-][N+](=O)\C=C(/NC)NCCSCC1=CC=C(CN(C)C)O1 VMXUWOKSQNHOCA-UKTHLTGXSA-N 0.000 description 6
- 241001598107 Imperata Species 0.000 description 5
- 230000000903 blocking effect Effects 0.000 description 5
- 230000000694 effects Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000003595 spectral effect Effects 0.000 description 2
- 241000132092 Aster Species 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000005670 electromagnetic radiation Effects 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/40—Means for monitoring or calibrating
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于信噪比识别雷达遮挡的方法,包含如下步骤:关闭雷达发射机,设定雷达发射信号的参数;将雷达天线仰角对准90度,发射电磁波遮挡;让雷达采用体扫模式进行扫描;保存雷达接收机收到的每个仰角层面的I/Q数据;将I/Q数据进行处理得到雷达天线处于不同仰角、不同方位上,雷达接收机输出的噪声功率图;利用已确定的门限值,将任意方向上接收到的信号与门限值比较,超过该门限,就可识别出该方向上有遮挡;根据比较结果绘制雷达周围不同仰角不同方位的遮挡图。本发明所述方法相比现有的技术,提高了计算的精度,解决了较近距离的遮挡计算问题。
Description
技术领域
本发明涉及一种识别雷达周围遮挡物分布的方法,尤其涉及一种利用信噪比(SNR,Signalto Noise Ratio)来计算和识别雷达周围遮挡物分布的方法,属于雷达信号处理研究领域。
背景技术
在固定站雷达选址和建站时,首要考虑的就是周围地形遮挡和建筑物遮挡问题。过多的遮挡会大大降低雷达的探测能力,使雷达发挥不出原有的作用,造成人力和财力的浪费。而对于移动式雷达来说,选择一个不受遮挡的观测场地是一项非常重要的工作。因此,雷达遮挡的计算合识别显得非常重要。
识别地物遮挡的方法有很多,基于数字高程模型(DEM)计算是当前的主流方法。万玉发等(2000)以人工读取高分辨率地图上四周地形的经度、纬度和海拔高度,辅以经纬仪从不同高度和不同方位测量四周新建高大建筑物相对于雷达站的方位和遮挡仰角,并转换回经纬度和海拔高度,通过余弦定理计算雷达单站的遮蔽角图,该方法虽然可行,但是效率较低。Kucera et al.(2004)分别采用10米、30米和100米分辨率的DEM数据,研究关岛地区的气象雷达的地形遮挡,发现地形资料分辨率越高,地形遮挡模拟效果越好。
但是,使用DEM数据计算雷达遮挡依然存在一些问题,首先是对数据的准确性要求很高,否则很容易产生较大偏差,并且对于距离雷达较近的建筑物等物体,由于DEM的数据精度问题,往往无法计算出正确的遮挡位置。此外,雷达波束在空间传输时,并不是射线,并且还有旁瓣的影响。这样,单纯从DEM数据分析雷达的遮挡,虽然具有很高的指导价值,但无法得到准确的结果。
任何物体都具有不断辐射、吸收、发射波束的本领。辐射出去的波束在各个波段是不同的,也就是具有一定的谱分布。这种谱分布与物体本身的特性及其温度有关,因而被称之为热辐射。为了研究不依赖于物质具体物性的热辐射规律,物理学定义了一种理想物体,称为黑体,以此作为热辐射研究的标准物体,并引入黑体辐射的概念。黑体辐射是指由理想放射物放射出来的辐射,在特定温度及特定波长放射最大量之辐射。理论上黑体会放射频谱上所有波长之波束。维恩位移定律是描述黑体电磁辐射能流密度的峰值波长与自身温度关系的定律。
在雷达工作时,雷达接收机的噪声会随方位而变化,规律为:有物体遮挡雷达波束的地方,茅草值为-66dBm左右;没有物体遮挡雷达波束的地方,即能看到回波的地方,茅草值为-67.4dBm左右,参见图1。这个现象是很容易解释,对雷达波束有遮挡的建筑物有温度,有黑体辐射。辐射被雷达的接收机收到,自然就增加了茅草值;而当雷达波束发射的方位上没有遮挡时,即天线对准太空,并且雷达天线仰角超过10度时,则太空的黑体辐射是很小的,因此茅草值就可以低1.5dB;当天线仰角比较低时,由于受地面的黑体辐射噪声,茅草值只能低0.5dB左右。
在雷达发射机不工作的情况下,在某个仰角做扫描时,雷达接收机的信噪比会随着周围地物的分布而变化,这主要是由于建筑物等遮挡物的黑体辐射作用。在有遮挡的时候,信噪比会变大,反之,则信噪比变小。因此,可以通过雷达接收机噪声随天线方位和仰角的变化,得出雷达四周近距离各个建筑物对雷达波束的遮挡分布图。
发明内容
本发明针对目前技术的不足,直接利用雷达接收机给出的信噪比来对遮挡进行计算和识别。该方法不但可以很好地计算雷达远处大面积的地形遮挡,对于雷达近处的建筑物等地物遮挡,也可以很好地计算与识别。
为实现以上的技术目的,本发明将采取以下的技术方案:
①关闭雷达发射机,设定雷达发射信号的参数。如重复周期,脉冲宽度,重复周期取500Hz至2000Hz,脉冲宽度取0.2μs至1μs。
②将雷达天线仰角对准90度,发射电磁波遮挡。该方法的目的是为了确定门限。
③让雷达采用体扫模式进行扫描;
④保存雷达接收机收到的每个仰角层面的I/Q数据;
⑤将I/Q数据进行处理得到雷达天线处于不同仰角、不同方位上,雷达接收机输出的噪声功率图。可以利用Fun_Read_NJU_CPol_IQ_RadarCharacter.m程序对I/Q数据进行处理。
⑥利用已确定的门限值,将任意方向上接收到的信号与门限值比较,超过了该门限,就可识别出该方向上有遮挡;
⑦根据比较结果绘制雷达周围不同仰角不同方位的遮挡图。
根据据以上的技术方案,可以实现以下的有益效果:
本发明所述方法,自动化程度高,方便易行,既可以作为雷达站建站前的评估,也可以用于建站后的效果测试。本发明所述方法相比现有的技术,提高了计算的精度,解决了较近距离的遮挡计算问题。
附图说明
图1为南京仙林雷达30km量程下的DEM数据图;
图2为南京仙林雷达四周30km范围内,各个地物相对于雷达的仰角图;
图3为南京仙林雷达在不同方位上,各个距离点的地物相对于雷达的仰角图;
图4为南京仙林雷达的遮挡角直角坐标显示图;
图5为南京仙林雷达的遮挡角极坐标显示图。
具体实施方式
下面结合附图对本发明的技术方案进行详细说明。
图1~5是以中国江苏南京仙林雷达(架设高度为60m)的遮挡情况为例,利用30m的DEM数据计算的结果。采用方案如下:
关闭雷达发射机,重复周期设为1000Hz,脉冲宽度设为0.5μs。
雷达天线仰角转至90°,即对准天顶方向。天线对准天顶,是因为天顶方向肯定没有遮挡,以便确定一个门限值。
设定一个雷达体扫模式,包括雷达扫描的层数和雷达天线的转速,让雷达按照设定好的体扫模式运转。
保存接收机输出的IQ数据,每个仰角层面产生的数据保存为一个文件,采用30m分辨率的DEM数据。
当雷达天线运转到某个方位/仰角时,如果两个接收通道的噪声功率出现变化,说明该方位/仰角处存在遮挡。原因为功率明显增加的区域是由于周围遮挡物的黑体辐射造成,利用确定的门限值进行比较,超过这个门限,就是有地物。
当雷达完成一个体扫后,就得到了各方位上不同仰角的噪声功率分布,利用黑体辐射的信噪比,就能计算出各仰角的遮挡情况。
绘制雷达周围的遮挡图。
上述技术方案中DEM数据的提取方法如下:
DEM数据是按照1度或5度,分成了一个个的文件。
因此需要先从庞大的DEM数据库中,根据用户设定的中心经纬度和量程、分辨率,选择需要的DEM数据。
程序见Fun_Select_DEM_Data.m。
注意:该程序包括了将不同的DEM文件的数据进行拼接的功能
程序流程如下:
先按照0、90、180、270度,计算四个方向上,经纬度的边界
如果使用30m分辨率的ASTER GDEM的DEM数据,则采用DEM_File_Directory_30m目录下的数据
如果使用90m分辨率的SRTM3的DEM数据,则采用DEM_File_Directory_90m目录下的数据,90m分辨率的DEM文件名的命令规则比较奇怪,见网上的介绍。
再根据四个边界,选择哪几个DEM文件
然后再将这些文件中,需要的数据拷贝到以雷达经纬度为中心的数组之中去。
上述方案中计算遮挡角的方法如下:
程序见Generate_Radar_Shield_Sector.m。
主要计算过程如下:
根据用户设定的经纬度和量程,选择需要的DEM数据
得到对应DEM数据中各个点的经纬度
调用LonToPolar,计算各个点相对于雷达的仰角
用imagesc绘制各个点相对于雷达的仰角
以横轴是方位(以0.1度为分辨率),竖轴是距离,用imagesc进行显示
计算各个方位上的最大的仰角,然后绘制极坐标下的遮挡图
注意:在设定雷达高度的时候,不能直接利用GPS得到的结果,而是要从DEM数据中,找出雷达所在经纬度的那个数据点,看看这个数据点所表示的高度是多少。和GPS对比,看两者相差多少。如果相差很多,则说明这个点的GPS值是不是搞错了。此时要利用Google Earth来进行确认。
Claims (3)
1.一种基于信噪比识别雷达遮挡的方法,其特征在于包含如下步骤:
关闭雷达发射机,设定雷达发射信号的参数;
将雷达天线仰角对准90度,发射电磁波遮挡;
让雷达采用体扫模式进行扫描;
保存雷达接收机收到的每个仰角层面的I/Q数据;
将I/Q数据进行处理得到雷达天线处于不同仰角、不同方位上,雷达接收机输出的噪声功率图;
利用已确定的门限值,将任意方向上接收到的信号与门限值比较,超过该门限,就可识别出该方向上有遮挡;
根据比较结果绘制雷达周围不同仰角不同方位的遮挡图。
2.一种基于信噪比识别雷达遮挡的方法,其特征在于所述设定雷达发射信号的参数为重复周期,脉冲宽度,重复周期取500Hz至2000Hz,脉冲宽度取0.2μs至1μs。
3.一种基于信噪比识别雷达遮挡的方法,其特征在于所述I/Q数据利用Fun_Read_NJU_CPol_IQ_RadarCharacter.m程序处理。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510036478.4A CN104614722B (zh) | 2015-01-20 | 2015-01-20 | 一种基于信噪比识别雷达遮挡的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510036478.4A CN104614722B (zh) | 2015-01-20 | 2015-01-20 | 一种基于信噪比识别雷达遮挡的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104614722A true CN104614722A (zh) | 2015-05-13 |
CN104614722B CN104614722B (zh) | 2018-05-08 |
Family
ID=53149256
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510036478.4A Active CN104614722B (zh) | 2015-01-20 | 2015-01-20 | 一种基于信噪比识别雷达遮挡的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104614722B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107796483A (zh) * | 2016-09-07 | 2018-03-13 | 桓达科技股份有限公司 | 远距自动强化信噪比的液位感测装置 |
WO2019024775A1 (zh) * | 2017-08-04 | 2019-02-07 | 维沃移动通信有限公司 | 天线波束切换感测系统、方法及移动终端 |
CN111650564A (zh) * | 2020-02-26 | 2020-09-11 | 惠州市德赛西威智能交通技术研究院有限公司 | 一种cfar后峰值点统计的车载雷达失效检测方法 |
CN111965607A (zh) * | 2020-07-09 | 2020-11-20 | 惠州市德赛西威智能交通技术研究院有限公司 | 一种车载雷达功能失效检测方法、装置及车辆 |
CN114690138A (zh) * | 2022-05-30 | 2022-07-01 | 南京隼眼电子科技有限公司 | 雷达性能监测方法、装置及存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6801153B2 (en) * | 2001-11-02 | 2004-10-05 | Spectrum Target Detection, Inc. | Spread spectrum radar with leak compensation at baseband |
CN101135729A (zh) * | 2007-09-04 | 2008-03-05 | 西安电子科技大学 | 基于支撑矢量机的合成孔径雷达遮挡目标识别方法 |
-
2015
- 2015-01-20 CN CN201510036478.4A patent/CN104614722B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6801153B2 (en) * | 2001-11-02 | 2004-10-05 | Spectrum Target Detection, Inc. | Spread spectrum radar with leak compensation at baseband |
CN101135729A (zh) * | 2007-09-04 | 2008-03-05 | 西安电子科技大学 | 基于支撑矢量机的合成孔径雷达遮挡目标识别方法 |
Non-Patent Citations (3)
Title |
---|
张记龙 等: ""太阳辐射对相干激光告警接收机信噪比影响的理论分析"", 《应用基础与工程科学学报》 * |
胡志群 等: ""C波段偏振雷达几种系统误差标定方法对比分析"", 《高原气象》 * |
荀平 等: ""山地雷达站营区防电磁辐射污染规划"", 《四川兵工学报》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107796483A (zh) * | 2016-09-07 | 2018-03-13 | 桓达科技股份有限公司 | 远距自动强化信噪比的液位感测装置 |
WO2019024775A1 (zh) * | 2017-08-04 | 2019-02-07 | 维沃移动通信有限公司 | 天线波束切换感测系统、方法及移动终端 |
CN111650564A (zh) * | 2020-02-26 | 2020-09-11 | 惠州市德赛西威智能交通技术研究院有限公司 | 一种cfar后峰值点统计的车载雷达失效检测方法 |
CN111965607A (zh) * | 2020-07-09 | 2020-11-20 | 惠州市德赛西威智能交通技术研究院有限公司 | 一种车载雷达功能失效检测方法、装置及车辆 |
CN111965607B (zh) * | 2020-07-09 | 2023-11-17 | 惠州市德赛西威智能交通技术研究院有限公司 | 一种车载雷达功能失效检测方法、装置及车辆 |
CN114690138A (zh) * | 2022-05-30 | 2022-07-01 | 南京隼眼电子科技有限公司 | 雷达性能监测方法、装置及存储介质 |
CN114690138B (zh) * | 2022-05-30 | 2022-08-02 | 南京隼眼电子科技有限公司 | 雷达性能监测方法、装置及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN104614722B (zh) | 2018-05-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104614722A (zh) | 一种基于信噪比识别雷达遮挡的方法 | |
US10871457B2 (en) | Determining material category based on the polarization of received signals | |
CN104133216A (zh) | 一种获取低空风廓线的雷达探测方法及装置 | |
CN103837867A (zh) | 一种利用ais信息进行高频雷达天线通道校正的方法 | |
JP6466480B2 (ja) | 探知装置、および、探知方法 | |
CN1633050A (zh) | 基于公众移动通信网络实现对空中目标监视的系统及方法 | |
CN110109072A (zh) | 一种多基地海上小目标雷达探测方法 | |
CN102798858B (zh) | 全息型主动式微波成像方法 | |
Yue et al. | Preliminary experimental results by the prototype of Sanya Incoherent Scatter Radar | |
Belmonte et al. | Analysis of wind turbines blockage on Doppler weather radar beams | |
CN109655828B (zh) | 一种多频一维综合孔径微波辐射计sst深度学习反演方法 | |
JP2011214975A (ja) | 物体変位監視装置、及び物体変位監視方法 | |
US20180209920A1 (en) | Water vapor observing apparatus | |
CN104898127A (zh) | 一种多发多收合成孔径激光雷达系统及其体制设计方法 | |
US10809369B2 (en) | Radar signal processing apparatus and radar signal processing method | |
CN112068087B (zh) | 一种岸基多通道雷达仿机载海杂波测量方法 | |
Kashiwayanagi et al. | Rapid 3D scanning high resolution X-band weather radar with active phased array antenna | |
Zhao et al. | A new inversion algorithm for backscatter ionogram and its experimental validation | |
US3553691A (en) | Long range position determination system | |
Kurdzo | Pulse compression waveforms and applications for weather radar | |
Tran et al. | Combining color and texture for a robust interactive segmentation algorithm | |
CN109031280A (zh) | 一种基于相控阵天线的时间式信标回收机及回收方法 | |
CN209387871U (zh) | 一种基于Vega的低空监测预警系统 | |
Yang et al. | The vertical ionosphere parameters inversion for high frequency surface wave radar | |
Belmonte et al. | Impact analysis of wind turbines blockage on Doppler weather radar |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |