CN112147587B - Radar beam azimuth center offshore calibration method - Google Patents
Radar beam azimuth center offshore calibration method Download PDFInfo
- Publication number
- CN112147587B CN112147587B CN202011038017.8A CN202011038017A CN112147587B CN 112147587 B CN112147587 B CN 112147587B CN 202011038017 A CN202011038017 A CN 202011038017A CN 112147587 B CN112147587 B CN 112147587B
- 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.)
- Active
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
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
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 flight path ab perpendicular to the optical axis direction, wherein the distance corresponds to the width of a radar half-power beam irradiation area as shown in fig. 2, and the longitude and latitude of a flight path endpoint 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 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 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 timeAzimuth 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 CN112147587A (en) | 2020-12-29 |
CN112147587B true 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 (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
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 |
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 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7218268B2 (en) * | 2003-05-14 | 2007-05-15 | Veridian Systems | Self-calibrating interferometric synthetic aperture radar altimeter |
US10175352B2 (en) * | 2015-05-12 | 2019-01-08 | Maxlinear, Inc. | Scalable architecture for an automotive radar system |
-
2020
- 2020-09-28 CN CN202011038017.8A patent/CN112147587B/en active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
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 |
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 |
---|
Analysis and Calibration of Crossed-Loop Antenna for Vessel DOA Estimation in HF Radar;Bo Lu;《 IEEE ANTENNAS AND WIRELESS PROPAGATION LETTERS》;20180131;第17卷(第1期);42-45 * |
BJÖRN LUND.On Shipboard Marine X-Band Radar Near-Surface Current "Calibration".《2015 American Meteorological Society》.2015,1928-1944. * |
Calibration measurement of shore-based radar with active reflector;Xin Li;《IEEE Xplore》;20181231;1-3 * |
一种星载合成孔径雷达微波辐射场定标接收机;赵建农等;《火控雷达技术》;20070625(第02期);56-59 * |
三星时差定位系统四站标校的布站方式研究;陈振林等;《飞行器测控学报》;20111231(第04期);86-90 * |
利用无人机标定雷达精度的新方法;刘冬利;《电讯技术》;20150228;第55卷(第2期);146-150 * |
卫星高度计直接绝对定标方法及其最新应用;管斌等;《地球物理学进展》;20171215(第06期);52-56 * |
海杂波测量定标的姿态修正数据处理方法;张玉石等;《电子与信息学报》;20150315(第03期);103-108 * |
海面雷达信号环境多路径效应建模与仿真技术研究;袁兴鹏等;《舰船电子对抗》;20131225(第06期);74-78 * |
船载外定标雷达波束照射区域定位方法研究;张浙东等;《电子设计工程》;20180605(第11期);52-57、62 * |
船载测控雷达无塔角度标校新方法;房新兵;《电讯技术》;20101231;第50卷(第12期);116-120 * |
高分辨率光学/SAR卫星几何辐射定标研究进展;张过等;《测绘学报》;20191215(第12期);134-153 * |
Also Published As
Publication number | Publication date |
---|---|
CN112147587A (en) | 2020-12-29 |
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 | |
CN109738902B (en) | High-precision autonomous acoustic navigation method for underwater high-speed target based on synchronous beacon mode | |
CN103985952B (en) | Boat-carrying A-E-C three-axis satellite communication antenna polarization angle of deviation real-time correcting method | |
CN105759263A (en) | High resolution satellite-borne squint SAR imaging method in large-scale scene | |
CN113960620B (en) | High-precision beam tracking system | |
CN102707269A (en) | Range walk correction method for airborne radar | |
CN108152811B (en) | Method for dynamically calibrating shipborne satellite communication earth station tracking phase based on probability | |
CN111766599B (en) | Laser radar device and method for measuring sea wave height in real time | |
CN105929377B (en) | A kind of higher-frequency radar ship azimuth method of estimation for intersecting loop antenna based on monopole | |
CN116148886A (en) | Buoy laser radar-based high-precision wind field inversion algorithm | |
CN112147587B (en) | Radar beam azimuth center offshore calibration method | |
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 | |
CN107643514A (en) | A kind of array calibration method of buoy base/shipborne radar based on direct wave | |
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 | |
CN112068087A (en) | Shore-based multi-channel radar aircraft-simulated sea clutter measurement method | |
CN109143219A (en) | A method of improving TWS radar pitch angle measurement precision | |
CN111025358A (en) | Orientation method based on single short baseline of navigation satellite signal | |
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 |
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 |