CN112147587A - Radar beam azimuth center offshore calibration method - Google Patents
Radar beam azimuth center offshore calibration method Download PDFInfo
- Publication number
- CN112147587A CN112147587A CN202011038017.8A CN202011038017A CN112147587A CN 112147587 A CN112147587 A CN 112147587A CN 202011038017 A CN202011038017 A CN 202011038017A CN 112147587 A CN112147587 A CN 112147587A
- Authority
- CN
- China
- Prior art keywords
- phi
- radar
- azimuth angle
- calculating
- azimuth
- 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
-
- 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
- G01S7/4004—Means for monitoring or calibrating of parts of a radar system
Abstract
The invention discloses a radar beam azimuth center offshore calibration method, which comprises the following steps: step 1, calculating a route; step 2, radar signal acquisition; step 3, calculating the azimuth corresponding to the maximum value of the signal; step 4, error judgment; and 5: adjusting a route; and 6, calculating photoelectric deviation. The radar beam calibration method disclosed by the invention is suitable for radars which are erected on the coast and are difficult to implement by the traditional calibration mode, and can enrich radar calibration means in actual environment.
Description
Technical Field
The invention belongs to the field of radar calibration, and particularly relates to a radar beam azimuth center offshore calibration method in the field.
Background
The measurement accuracy of the radar detection target determines the working capacity of the radar to some extent, and is a very important technical index. In practical applications, the radar measurement errors mainly include random errors and systematic errors. Random errors are usually suppressed by means of data statistics, while systematic errors need to be eliminated by some technical means.
Among the many systematic errors, beam pointing error is the most common one. The radar antenna is installed in place and is caused by the fact that the optical axis and the electric axis are not in the same direction, wherein the optical axis is a visual axis of optical equipment which is installed on the antenna and used as a shaft system reference, and the electric axis is the beam center direction of the radar antenna. Currently, the primary way to calibrate both deviations is based on a reference tower. The erection of the reference tower has strict requirements on distance and height (see Zhang Zhongsheng, Wangzhui, radar calibration technology [ M ], national defense industry publishers, 2017 and GJB 3153-. In far field conditions, the calibration tower height must be higher than the radar antenna to reduce ground multipath effects. However, for radars in certain specific environments, it is difficult to erect calibration towers around the radar site. For example, a sea clutter observation radar erected at the top of a sea island is difficult to realize in the traditional calibration mode, and a new method needs to be found.
Disclosure of Invention
The invention provides a radar beam azimuth center offshore calibration method which is used for solving the problem that a sea clutter observation radar erected at the top of a sea island cannot be calibrated by a calibration tower.
The invention adopts the following technical scheme:
the improvement of the method for calibrating the azimuth center of the radar beam at sea comprises the following steps:
step 1, calculating a route: according to the longitude and latitude LatR, LonR, height H and optical axis azimuth phi of the radar0A downward depression angle alpha of a radar beam and a half-power beam width theta of the radar3dBDesigning a flight path ab vertical to the direction of the optical axis, wherein the distance corresponds to the width of the irradiation area of the radar half-power beam, and the longitude and latitude of the end point of the flight path can be obtained through a Vincenty forward solution formula;
step 3, calculating the azimuth corresponding to the maximum signal value: according to the time t corresponding to the maximum value of the received radar signal, finding corresponding ship longitude and latitude information LatS and LonS in the track information, and according to radar longitude and latitude LatR, LonR and Vincenty inverse solution formula, calculating real-time azimuth angle phi1Drawing a curve of the received signal changing along with the azimuth angle;
step 4, error judgment: calculating the azimuth angle phi corresponding to the maximum value of the received signal1And the azimuth angle phi in step 10If Δ φ is | φ1-φ0|>θ3dB10, then the azimuth angle phi in step 1 is determined0Is changed to phi1Recalculating the course and re-executing the step 2 and the step 3 to obtain the azimuth angle phi2And calculate phi2And the last azimuth angle phi1So as to update a difference of | phi ═ phi1-φ0|>θ3dBThe value of delta phi in/10,
and 5: adjusting a route: after n +1 times of adjustment, the requirement that delta phi is equal to | phi is metn-φn-1|≤θ3dBAt/10 th, phinNamely the azimuth angle of the electric axis;
step 6, calculating photoelectric deviation: azimuth angle phi of optical axis0Azimuth angle phi with respect to the electric axisnThe difference of (a) is the photoelectric deviation delta phi', phi ═ phi0-φn。
The invention has the beneficial effects that:
the radar beam calibration method disclosed by the invention is suitable for radars which are erected on the coast and are difficult to implement by the traditional calibration mode, and can enrich radar calibration means in actual environment.
Drawings
Fig. 1 is a schematic flow chart of a radar beam calibration method disclosed in embodiment 1 of the present invention;
fig. 2 is a diagram of the effect of the amplitude mean estimation method.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention will be described in further detail below with reference to the accompanying drawings and examples. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
Embodiment 1, as shown in fig. 1, this embodiment discloses a radar beam azimuth center maritime calibration method, including the following steps:
step 1, calculating a route: according to the radar longitude and latitude (LatR, LonR), height H and optical axis azimuth angle phi0A downward depression angle alpha of a radar beam and a half-power beam width theta of the radar3dBDesigning a route ab perpendicular to the optical axis direction, as shown in FIG. 2, wherein the distance corresponds to the irradiation area of the radar half-power beamThe width of the domain and the longitude and latitude of the end point of the air route can be obtained through a Vincenty forward solution formula;
step 3, calculating the azimuth corresponding to the maximum signal value: according to the time t corresponding to the maximum value of the received radar signal, finding corresponding ship longitude and latitude information (LatS, LonS) in the track information, and according to radar longitude and latitude (LatR, LonR) and Vincenty inverse solution formula, calculating real-time azimuth angle phi1Drawing a curve of the received signal changing along with the azimuth angle;
step 4, error judgment: calculating the azimuth angle phi corresponding to the maximum value of the received signal1And the azimuth angle phi in step 10If the error is equal to | phi ═ phi1-φ0|>θ3dB10, then the azimuth angle phi in step 1 is determined0Is changed to phi1Recalculating the course and re-executing the step 2 and the step 3 to obtain the azimuth angle phi2And calculate phi2And the last azimuth angle phi1So as to update a difference of | phi ═ phi1-φ0|>θ3dBThe value of delta phi in/10,
and 5: adjusting a route: after n +1 times of adjustment, the error delta phi is equal to | phin-φn-1|≤θ3dBAt/10 th, phinNamely the azimuth angle of the electric axis;
step 6, calculating photoelectric deviation: azimuth angle phi of optical axis0Azimuth angle phi with respect to the electric axisnThe difference of (a) is the photoelectric deviation delta phi', phi ═ phi0-φn。
Claims (1)
1. A radar beam azimuth center offshore calibration method is characterized by comprising the following steps:
step 1, calculating a route: according to the longitude and latitude LatR, LonR, height H and optical axis azimuth phi of the radar0Downward depression angle alpha of radar beam and radar half-power waveBeam width θ3dBDesigning a flight path ab vertical to the direction of the optical axis, wherein the distance corresponds to the width of the irradiation area of the radar half-power beam, and the longitude and latitude of the end point of the flight path can be obtained through a Vincenty forward solution formula;
step 2, radar signal acquisition: guiding a ship to move along a course line at a constant speed through navigation equipment, carrying an electromagnetic signal receiver and a receiving antenna on the ship, ensuring that the beam center of the receiving antenna points to a radar by utilizing a ship-mounted stable platform, receiving and recording radar signals in real time, and synchronously recording track information;
step 3, calculating the azimuth corresponding to the maximum signal value: according to the time t corresponding to the maximum value of the received radar signal, finding corresponding ship longitude and latitude information LatS and LonS in the track information, and according to radar longitude and latitude LatR, LonR and Vincenty inverse solution formula, calculating real-time azimuth angle phi1Drawing a curve of the received signal changing along with the azimuth angle;
step 4, error judgment: calculating the azimuth angle phi corresponding to the maximum value of the received signal1And the azimuth angle phi in step 10If Δ φ is | φ1-φ0|>θ3dB10, then the azimuth angle phi in step 1 is determined0Is changed to phi1Recalculating the course and re-executing the step 2 and the step 3 to obtain the azimuth angle phi2And calculate phi2And the last azimuth angle phi1So as to update a difference of | phi ═ phi1-φ0|>θ3dBThe value of delta phi in/10,
and 5: adjusting a route: after n +1 times of adjustment, the requirement that delta phi is equal to | phi is metn-φn-1|≤θ3dBAt/10 th, phinNamely the azimuth angle of the electric axis;
step 6, calculating photoelectric deviation: azimuth angle phi of optical axis0Azimuth angle phi with respect to the electric axisnThe difference of (a) is the photoelectric deviation delta phi', phi ═ phi0-φn。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011038017.8A CN112147587B (en) | 2020-09-28 | 2020-09-28 | Radar beam azimuth center offshore calibration method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011038017.8A CN112147587B (en) | 2020-09-28 | 2020-09-28 | Radar beam azimuth center offshore calibration method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112147587A true CN112147587A (en) | 2020-12-29 |
CN112147587B CN112147587B (en) | 2022-02-25 |
Family
ID=73894943
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011038017.8A Active CN112147587B (en) | 2020-09-28 | 2020-09-28 | Radar beam azimuth center offshore calibration method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112147587B (en) |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040227658A1 (en) * | 2003-05-14 | 2004-11-18 | Vandenberg Norman | Self-calibrating interferometric synthetic aperture radar altimeter |
CN102955155A (en) * | 2011-08-26 | 2013-03-06 | 中国科学院空间科学与应用研究中心 | Distributed active phased array radar and beam forming method thereof |
CN103091666A (en) * | 2011-11-07 | 2013-05-08 | 中国科学院电子学研究所 | Onboard P-waveband polarization synthetic aperture radar (SAR) calibration method under non-ideal calibrator condition |
CN103116158A (en) * | 2013-01-18 | 2013-05-22 | 中国人民解放军海军航空工程学院 | Pulse passive bistatic radar target locating method |
CN104730503A (en) * | 2015-03-18 | 2015-06-24 | 中国科学院电子学研究所 | Method for determining influence on scaling by high-resolution SAR reference target RCS and compensation method |
US20160334504A1 (en) * | 2015-05-12 | 2016-11-17 | Maxlinear, Inc. | Scalable Architecture for an Automotive Radar System |
CN106249766A (en) * | 2016-09-09 | 2016-12-21 | 中国科学院电子学研究所 | Remotely monitoring corner reflector system |
CN107076835A (en) * | 2014-06-26 | 2017-08-18 | 罗伯特·博世有限公司 | Radar surveying method with different sighting distances |
CN110133609A (en) * | 2019-04-04 | 2019-08-16 | 杭州电子科技大学 | It is mobile under a kind of airborne external sort algorithm Status unknown to receive source path optimization method |
CN110907907A (en) * | 2019-10-19 | 2020-03-24 | 中国电波传播研究所(中国电子科技集团公司第二十二研究所) | Sea clutter Doppler spectrum characteristic analysis and comparison method |
CN111122985A (en) * | 2019-12-19 | 2020-05-08 | 中国电波传播研究所(中国电子科技集团公司第二十二研究所) | Autonomous underwater electromagnetic signal measuring device and measuring method |
CN111566505A (en) * | 2019-06-28 | 2020-08-21 | 深圳市大疆创新科技有限公司 | Calibration method and system for installation angle of millimeter wave radar component and movable platform |
-
2020
- 2020-09-28 CN CN202011038017.8A patent/CN112147587B/en active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040227658A1 (en) * | 2003-05-14 | 2004-11-18 | Vandenberg Norman | Self-calibrating interferometric synthetic aperture radar altimeter |
CN102955155A (en) * | 2011-08-26 | 2013-03-06 | 中国科学院空间科学与应用研究中心 | Distributed active phased array radar and beam forming method thereof |
CN103091666A (en) * | 2011-11-07 | 2013-05-08 | 中国科学院电子学研究所 | Onboard P-waveband polarization synthetic aperture radar (SAR) calibration method under non-ideal calibrator condition |
CN103116158A (en) * | 2013-01-18 | 2013-05-22 | 中国人民解放军海军航空工程学院 | Pulse passive bistatic radar target locating method |
CN107076835A (en) * | 2014-06-26 | 2017-08-18 | 罗伯特·博世有限公司 | Radar surveying method with different sighting distances |
CN104730503A (en) * | 2015-03-18 | 2015-06-24 | 中国科学院电子学研究所 | Method for determining influence on scaling by high-resolution SAR reference target RCS and compensation method |
US20160334504A1 (en) * | 2015-05-12 | 2016-11-17 | Maxlinear, Inc. | Scalable Architecture for an Automotive Radar System |
CN106249766A (en) * | 2016-09-09 | 2016-12-21 | 中国科学院电子学研究所 | Remotely monitoring corner reflector system |
CN110133609A (en) * | 2019-04-04 | 2019-08-16 | 杭州电子科技大学 | It is mobile under a kind of airborne external sort algorithm Status unknown to receive source path optimization method |
CN111566505A (en) * | 2019-06-28 | 2020-08-21 | 深圳市大疆创新科技有限公司 | Calibration method and system for installation angle of millimeter wave radar component and movable platform |
CN110907907A (en) * | 2019-10-19 | 2020-03-24 | 中国电波传播研究所(中国电子科技集团公司第二十二研究所) | Sea clutter Doppler spectrum characteristic analysis and comparison method |
CN111122985A (en) * | 2019-12-19 | 2020-05-08 | 中国电波传播研究所(中国电子科技集团公司第二十二研究所) | Autonomous underwater electromagnetic signal measuring device and measuring method |
Non-Patent Citations (12)
Title |
---|
BJÖRN LUND: "On Shipboard Marine X-Band Radar Near-Surface Current "Calibration"", 《2015 AMERICAN METEOROLOGICAL SOCIETY》 * |
BO LU: "Analysis and Calibration of Crossed-Loop Antenna for Vessel DOA Estimation in HF Radar", 《 IEEE ANTENNAS AND WIRELESS PROPAGATION LETTERS》 * |
XIN LI: "Calibration measurement of shore-based radar with active reflector", 《IEEE XPLORE》 * |
刘冬利: "利用无人机标定雷达精度的新方法", 《电讯技术》 * |
张浙东等: "船载外定标雷达波束照射区域定位方法研究", 《电子设计工程》 * |
张玉石等: "海杂波测量定标的姿态修正数据处理方法", 《电子与信息学报》 * |
张过等: "高分辨率光学/SAR卫星几何辐射定标研究进展", 《测绘学报》 * |
房新兵: "船载测控雷达无塔角度标校新方法", 《电讯技术》 * |
管斌等: "卫星高度计直接绝对定标方法及其最新应用", 《地球物理学进展》 * |
袁兴鹏等: "海面雷达信号环境多路径效应建模与仿真技术研究", 《舰船电子对抗》 * |
赵建农等: "一种星载合成孔径雷达微波辐射场定标接收机", 《火控雷达技术》 * |
陈振林等: "三星时差定位系统四站标校的布站方式研究", 《飞行器测控学报》 * |
Also Published As
Publication number | Publication date |
---|---|
CN112147587B (en) | 2022-02-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103389493B (en) | Multi-beam single-pulse angle measuring method based on beam selection method | |
CN104391281A (en) | Method for improving sky-wave radar sea surface ship target tracking and positioning precision | |
CN111708038A (en) | Unmanned ship laser radar point cloud data correction method based on attitude sensor and GNSS | |
CN111060943B (en) | Weight determining method conforming to practical GNSS observation value | |
CN109738902B (en) | High-precision autonomous acoustic navigation method for underwater high-speed target based on synchronous beacon mode | |
CN105759263A (en) | High resolution satellite-borne squint SAR imaging method in large-scale scene | |
CN103985952B (en) | Boat-carrying A-E-C three-axis satellite communication antenna polarization angle of deviation real-time correcting method | |
CN113960620B (en) | High-precision beam tracking system | |
CN102707269A (en) | Range walk correction method for airborne radar | |
CN103529451B (en) | Method for calibrating coordinate position of seabed transponder of water-surface mother ship | |
CN116148886A (en) | Buoy laser radar-based high-precision wind field inversion algorithm | |
CN111766599A (en) | Laser radar device and method for measuring sea wave height in real time | |
CN112147587B (en) | Radar beam azimuth center offshore calibration method | |
CN105929377A (en) | High-frequency radar ship azimuth estimation method based on monopole intersected loop antenna | |
CN111638514B (en) | Unmanned aerial vehicle height measurement method and unmanned aerial vehicle navigation filter | |
CN116772838B (en) | Inertial navigation error compensation method for mechanical phased array antenna | |
CN100585429C (en) | Passive channel correcting method based on non-linear antenna array | |
CN111190178A (en) | Multi-beam sequence amplitude comparison height measurement method of phase frequency scanning three-coordinate radar | |
CN113534130B (en) | Multi-station radar multi-target data association method based on sight angle | |
CN114355396A (en) | USB shafting parameter calibration method based on Beidou system | |
CN109143219A (en) | A method of improving TWS radar pitch angle measurement precision | |
CN114578310B (en) | Method for judging and processing dispersion effect of satellite-borne SAR ionized layer | |
CN112485769B (en) | Wave beam self-adaptive correction method based on phased array secondary radar | |
RU2791981C2 (en) | Method for dynamic adjustment of antenna array with electronic scanning of ship radar station | |
CN116148840B (en) | Dynamic high-precision positioning environment symbiotic stealth radar system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |