CN107255817A - ISAR image orientations based on Parameter Estimation Method are to calibrating method - Google Patents
ISAR image orientations based on Parameter Estimation Method are to calibrating method Download PDFInfo
- Publication number
- CN107255817A CN107255817A CN201710529595.3A CN201710529595A CN107255817A CN 107255817 A CN107255817 A CN 107255817A CN 201710529595 A CN201710529595 A CN 201710529595A CN 107255817 A CN107255817 A CN 107255817A
- Authority
- CN
- China
- Prior art keywords
- isar
- point
- mrow
- msub
- special aobvious
- 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
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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/904—SAR modes
-
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/904—SAR modes
- G01S13/9064—Inverse SAR [ISAR]
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The present invention relates to radar image field of locating technology, a kind of ISAR image orientations based on Parameter Estimation Method are specifically provided to calibrating method, ISAR raw radar datas are obtained first and carried out motion compensation and carried out orientation FFT imagings to obtain original I SAR image, then the maximum aobvious point of spy of the top n picture contrast in all special aobvious points of the original I SAR image is obtained by special aobvious point method and be recorded in set φ, the instantaneous Doppler frequency of the strong scattering point of range cell where obtaining each special aobvious point in φ, try to achieve again and the corresponding frequency modulation rate of each special aobvious point in set φ, chirp rate k is obtained after fitting, effective rotational angular velocity ω e estimate is calculated finally by chirp rate k, and then obtain azimuth resolution ρ a estimate, orientation calibration is carried out to ISAR images by way of orientation interpolation, obtain the ISAR images after the ship classification demarcation mark.
Description
Technical field
The present invention relates to radar image field of locating technology, ISAR image orientations more particularly to based on Parameter Estimation Method to
Calibrating method.
Background technology
ISAR (Inverse Synthetic Aperture Radar, abbreviation ISAR) is that one kind can be complete
Weather, round-the-clock, the non-cooperative moving targets such as vehicle, naval vessel, aircraft, guided missile are imaged at a distance, obtain target high
The radar of image in different resolution.ISAR imagings often assume that target rotation angle is smaller in a short time and approximate easy motion, using away from
From-Doppler (R-D) principle imaging.But RD imaging results only reflect projection information of the target in distance-Doppler plane, and
It can not provide target real two-dimensional size information, in order to the further identification of target, the research of ISAR orientation calibration technologies
It is most important.
ISAR image orientations of the prior art are to obtain target by the calculating of the parameters such as distance, angle to exist to calibrating method
Effective corner during imaging, and then laterally calibration, but the robustness of this method is not strong is realized, method performance is by imaging resolution
Etc. the restriction of factor.
The content of the invention
To overcome at least one defect that above-mentioned prior art is present, the invention provides a kind of based on Parameter Estimation Method
ISAR image orientations comprise the following steps to calibrating method:
Step one, ISAR is obtained to the raw radar data of offshore vessel class target imaging, and ISAR is ISAR,
Motion compensation is carried out to the raw radar data, orientation FFT imagings are carried out to the data after compensation, original I SAR is obtained
Image, FFT is fast Fourier transform;
Step 2, the top n image comparison in all special aobvious points of the original I SAR image is obtained by special aobvious point method
The maximum spy of degree shows point, and N number of special aobvious point recorded in set φ;
Step 3, the instantaneous Doppler frequency of the strong scattering point of range cell where obtaining each special aobvious point in the set φ
Rate;
Step 4, tries to achieve a frequency modulation rate corresponding with the aobvious point of each spy in set φ, by institute by the instantaneous Doppler frequency
There is the value fitting of frequency modulation rate in alignment, obtain the chirp rate k of the straight line;
Step 5, when calculating the selected imaging moment in step 3 by the chirp rate k, the ship classification target
Effective rotational angular velocity ω e estimate, calculation formula is:
Wherein c is velocity of electromagnetic wave, and f0 is radar carrier frequency;
The estimate of the ship class target azimuth resolution ρ a, calculation formula are obtained by effective rotational angular velocity ω e
For:
Wherein λ is wavelength, and T is that ISAR is imaged integration time;
Orientation calibration is carried out to ISAR images by way of orientation interpolation, obtained after the ship classification demarcation mark
ISAR images.
It is preferred that, the N in step 2 is not less than 10.
It is preferred that, pass through a distance list where each special aobvious point in Matched Fourier Transform method acquisition set φ in step 3
The instantaneous Doppler frequency of the strong scattering point of member.
The ISAR image orientations based on Parameter Estimation Method that the present invention is provided have the advantages that to calibrating method:
1st, range cell where isolated strong scattering point being chosen from image area using special aobvious point method carries out effective angle of rotation speed
Degree estimation, data volume and operand are smaller, it is easy to Project Realization;
2nd, the instantaneous Doppler frequency of special aobvious point is estimated using the method for matching Fourier's change, frequency modulation is then tried to achieve
Rate, efficiently reduces evaluated error, and precision is higher.
Brief description of the drawings
Fig. 1 is flow chart of the ISAR image orientations based on Parameter Estimation Method to calibrating method.
Embodiment
To make the purpose, technical scheme and advantage of the invention implemented clearer, below in conjunction with the embodiment of the present invention
Accompanying drawing, the technical scheme in the embodiment of the present invention is further described in more detail.
It should be noted that:The embodiments described below with reference to the accompanying drawings are exemplary, it is intended to for explaining this hair
It is bright, and be not considered as limiting the invention.In the accompanying drawings, same or similar label represents same or like from beginning to end
Element or element with same or like function.Described embodiment is a part of embodiment of the invention, rather than entirely
The embodiment in portion, in the case where not conflicting, the feature in embodiment and embodiment in the application can be mutually combined.It is based on
Embodiment in the present invention, it is every other that those of ordinary skill in the art are obtained under the premise of creative work is not made
Embodiment, belongs to the scope of protection of the invention.
As shown in figure 1, a kind of ISAR image orientations based on Parameter Estimation Method that the present invention is provided are to calibrating method, in spy
On the basis of aobvious point method, the instantaneous Doppler frequency of special aobvious point is estimated by Matched Fourier Transform, frequency modulation rate is then tried to achieve,
Effective rotational angular velocity and azimuth resolution are estimated, realizes that orientation is calibrated.The ISAR image orientations are to calibrating method bag
Include following steps:
Step one, ISAR is obtained to the raw radar data of offshore vessel class target imaging, and ISAR is ISAR,
Motion compensation is carried out to the raw radar data, it is therefore an objective to error of the naval vessel relative to radar motion is eliminated, after compensation
Data carry out orientation FFT imagings, obtain original I SAR image, FFT is fast Fourier transform;
Step 2, the top n image comparison in all special aobvious points of the original I SAR image is obtained by special aobvious point method
The maximum spy of degree shows point, N number of special aobvious point recorded in set φ, N is not less than 10 in the present embodiment preferably;
Step 3, selectes imaging moment, the strong scattering point of range cell where obtaining each special aobvious point in the set φ
Instantaneous Doppler frequency, pass through Matched Fourier Transform method and obtain each spy in set φ in the present embodiment preferably
The instantaneous Doppler frequency of the strong scattering point of range cell where aobvious point;
Step 4, tries to achieve a frequency modulation rate corresponding with the aobvious point of each spy in set φ, by institute by the instantaneous Doppler frequency
There is the value fitting of frequency modulation rate in alignment, obtain the chirp rate k of the straight line;
Step 5, when calculating the selected imaging moment in step 3 by the chirp rate k, the ship classification target
Effective rotational angular velocity ωeEstimate, calculation formula is:
Wherein c is velocity of electromagnetic wave, f0For radar carrier frequency;
Pass through effective rotational angular velocity ωeObtain the ship class target azimuth resolution ρaEstimate, calculation formula
For:
Wherein λ is wavelength, and T is that ISAR is imaged integration time;
Orientation calibration is carried out to ISAR images by way of orientation interpolation, i.e., according to distance to pixel size, the side of making
Position is consistent with its to pixel size, obtains the ISAR images after the ship classification demarcation mark.
The foregoing is only a specific embodiment of the invention, but protection scope of the present invention is not limited thereto, any
Those familiar with the art the invention discloses technical scope in, the change or replacement that can be readily occurred in, all should
It is included within the scope of the present invention.Therefore, protection scope of the present invention should using the scope of the claims as
It is accurate.
Claims (3)
1. a kind of ISAR image orientations based on Parameter Estimation Method are to calibrating method, it is characterised in that comprise the following steps:
Step one, ISAR is obtained to the raw radar data of offshore vessel class target imaging, and ISAR is ISAR, to institute
State raw radar data and carry out motion compensation, orientation FFT imagings are carried out to the data after compensation, obtain original I SAR figures
Picture, FFT is fast Fourier transform;
Step 2, by the top n picture contrast in all special aobvious points of the special aobvious point method acquisition original I SAR image most
Big spy shows point, and N number of special aobvious point recorded in set φ;
Step 3, the instantaneous Doppler frequency of the strong scattering point of range cell where obtaining each special aobvious point in the set φ;
Step 4, tries to achieve a frequency modulation rate corresponding with the aobvious point of each spy in set φ by the instantaneous Doppler frequency, is adjusted all
The value fitting of frequency is in alignment, obtains the chirp rate k of the straight line;
Step 5, when calculating the selected imaging moment in step 3 by the chirp rate k, the ship classification target is effective
Rotational angular velocity ωeEstimate, calculation formula is:
<mrow>
<mi>k</mi>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>f</mi>
<mn>0</mn>
</msub>
</mrow>
<mi>c</mi>
</mfrac>
<msubsup>
<mi>&omega;</mi>
<mi>e</mi>
<mn>2</mn>
</msubsup>
<mo>;</mo>
</mrow>
Wherein c is velocity of electromagnetic wave, f0For radar carrier frequency;
Pass through effective rotational angular velocity ωeObtain the ship class target azimuth resolution ρaEstimate, calculation formula is:
<mrow>
<msub>
<mi>&rho;</mi>
<mi>a</mi>
</msub>
<mo>=</mo>
<mfrac>
<mi>&lambda;</mi>
<mrow>
<mn>2</mn>
<msub>
<mi>&omega;</mi>
<mi>e</mi>
</msub>
<mi>T</mi>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
Wherein λ is wavelength, and T is that ISAR is imaged integration time;
Orientation calibration is carried out to ISAR images by way of orientation interpolation, the ISAR figures after the ship classification demarcation mark are obtained
Picture.
2. ISAR image orientations according to claim 1 are to calibrating method, it is characterised in that the N in step 2 is not less than
10。
3. ISAR image orientations according to claim 1 are to calibrating method, it is characterised in that by matching Fu in step 3
In leaf transformation method obtain each special aobvious point place range cell in set φ strong scattering point an instantaneous Doppler frequency.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710529595.3A CN107255817B (en) | 2017-07-02 | 2017-07-02 | ISAR image azimuth direction calibration method based on parameter estimation method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710529595.3A CN107255817B (en) | 2017-07-02 | 2017-07-02 | ISAR image azimuth direction calibration method based on parameter estimation method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107255817A true CN107255817A (en) | 2017-10-17 |
CN107255817B CN107255817B (en) | 2021-04-20 |
Family
ID=60025927
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710529595.3A Active CN107255817B (en) | 2017-07-02 | 2017-07-02 | ISAR image azimuth direction calibration method based on parameter estimation method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107255817B (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109884631A (en) * | 2019-02-21 | 2019-06-14 | 中国科学院电子学研究所 | The method and system of satellite-borne synthetic aperture radar internal calibration data processing |
CN110346797A (en) * | 2019-07-15 | 2019-10-18 | 中国人民解放军32181部队 | A kind of two dimension inverse synthetic aperture radar imaging method and system |
CN112965057A (en) * | 2021-02-08 | 2021-06-15 | 哈尔滨工业大学 | Method for estimating and calibrating rotating speed of near-field turntable target |
CN114966693A (en) * | 2022-07-20 | 2022-08-30 | 南京信息工程大学 | Airborne ship target ISAR refined imaging method based on deep learning |
CN115902895A (en) * | 2022-12-27 | 2023-04-04 | 深圳大学 | Single-bit high-precision signal parameter estimation method for inverse synthetic aperture radar |
WO2024087036A1 (en) * | 2022-10-25 | 2024-05-02 | 华为技术有限公司 | Imaging method and system based on radio frequency signal transmission, and communication apparatus |
CN115902895B (en) * | 2022-12-27 | 2024-07-30 | 深圳大学 | Single-ratio ultra-high precision signal parameter estimation method for inverse synthetic aperture radar |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101000374A (en) * | 2006-12-27 | 2007-07-18 | 西安电子科技大学 | Interference reverse synthetic aperture radarimaging method based on multiple-spot |
JP2011174875A (en) * | 2010-02-25 | 2011-09-08 | Mitsubishi Electric Corp | Passive radar device |
CN103698765A (en) * | 2013-12-27 | 2014-04-02 | 西安电子科技大学 | ISAR (inverse synthetic aperture radar) imaging orientation calibration method |
CN104391286A (en) * | 2014-12-10 | 2015-03-04 | 中国电子科技集团公司第二十研究所 | Retrosynthetic aperture radar azimuth calibration method |
-
2017
- 2017-07-02 CN CN201710529595.3A patent/CN107255817B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101000374A (en) * | 2006-12-27 | 2007-07-18 | 西安电子科技大学 | Interference reverse synthetic aperture radarimaging method based on multiple-spot |
JP2011174875A (en) * | 2010-02-25 | 2011-09-08 | Mitsubishi Electric Corp | Passive radar device |
CN103698765A (en) * | 2013-12-27 | 2014-04-02 | 西安电子科技大学 | ISAR (inverse synthetic aperture radar) imaging orientation calibration method |
CN104391286A (en) * | 2014-12-10 | 2015-03-04 | 中国电子科技集团公司第二十研究所 | Retrosynthetic aperture radar azimuth calibration method |
Non-Patent Citations (2)
Title |
---|
汪玲: "逆合成孔径雷达成像关键技术研究", 《中国博士学位论文全文数据库信息科技辑》 * |
荆腾: "ISAR运动补偿及舰船成像技术研究", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109884631A (en) * | 2019-02-21 | 2019-06-14 | 中国科学院电子学研究所 | The method and system of satellite-borne synthetic aperture radar internal calibration data processing |
CN109884631B (en) * | 2019-02-21 | 2020-11-10 | 中国科学院电子学研究所 | Method and system for processing calibration data in satellite-borne synthetic aperture radar |
CN110346797A (en) * | 2019-07-15 | 2019-10-18 | 中国人民解放军32181部队 | A kind of two dimension inverse synthetic aperture radar imaging method and system |
CN110346797B (en) * | 2019-07-15 | 2021-05-25 | 中国人民解放军32181部队 | Two-dimensional inverse synthetic aperture radar imaging method and system |
CN112965057A (en) * | 2021-02-08 | 2021-06-15 | 哈尔滨工业大学 | Method for estimating and calibrating rotating speed of near-field turntable target |
CN112965057B (en) * | 2021-02-08 | 2023-12-08 | 哈尔滨工业大学 | Method for estimating and calibrating rotation speed of near-field turntable target |
CN114966693A (en) * | 2022-07-20 | 2022-08-30 | 南京信息工程大学 | Airborne ship target ISAR refined imaging method based on deep learning |
CN114966693B (en) * | 2022-07-20 | 2022-11-04 | 南京信息工程大学 | Airborne ship target ISAR refined imaging method based on deep learning |
WO2024087036A1 (en) * | 2022-10-25 | 2024-05-02 | 华为技术有限公司 | Imaging method and system based on radio frequency signal transmission, and communication apparatus |
CN115902895A (en) * | 2022-12-27 | 2023-04-04 | 深圳大学 | Single-bit high-precision signal parameter estimation method for inverse synthetic aperture radar |
CN115902895B (en) * | 2022-12-27 | 2024-07-30 | 深圳大学 | Single-ratio ultra-high precision signal parameter estimation method for inverse synthetic aperture radar |
Also Published As
Publication number | Publication date |
---|---|
CN107255817B (en) | 2021-04-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107255817A (en) | ISAR image orientations based on Parameter Estimation Method are to calibrating method | |
CN106324597B (en) | The translational compensation and imaging method of big corner ISAR radar based on PFA | |
Martorella et al. | ISAR based technique for refocusing non-cooperative targets in SAR images | |
US8994584B2 (en) | Autofocus-based compensation (ABC) system and method for a hovering ground moving target indication (GMTI) sensor | |
CN106597437B (en) | The airborne CW with frequency modulation SAR motion compensation of high frequency and image processing method | |
CN109856635B (en) | CSAR ground moving target refocusing imaging method | |
CN101498788B (en) | Target rotation angle estimating and transverse locating method for inverse synthetic aperture radar | |
CN102590812A (en) | SAR (synthetic aperture radar) real-time imaging method based on frequency modulated continuous wave | |
Li et al. | Autofocus correction of residual RCM for VHR SAR sensors with light-small aircraft | |
CN107390181B (en) | Radar high-resolution imaging method based on multi-beam scanning | |
CN101620272B (en) | Target rotate speed estimation method of inverse synthetic aperture radar (ISAR) | |
CN105974410A (en) | Multi-ship target SAR and ISAR hybrid imaging method of airborne radar | |
US7586434B1 (en) | Method and system for mapping a target scene using scanning radar | |
Kang et al. | Bistatic-ISAR distortion correction and range and cross-range scaling | |
CN105954751A (en) | Ka FMCW SAR moving object detection method and apparatus | |
CN102565772B (en) | Marine dynamic information extraction method on basis of SAR (Synthetic Aperture Radar) sub-aperture sequence images | |
Sjanic | Navigation and SAR Auto-focusing in a Sensor Fusion Framework | |
Kang et al. | ISAR cross-range scaling via joint estimation of rotation center and velocity | |
CN109799502B (en) | Two-dimensional self-focusing method suitable for filtering back projection algorithm | |
CN115166714A (en) | Single-channel SAR moving ship two-dimensional speed estimation and repositioning method and device | |
CN101846741B (en) | Inverse synthetic aperture radar imaging data segment selecting method | |
Felguera-Martin et al. | Interferometric ISAR imaging on maritime target applications: Simulation of realistic targets and dynamics | |
Nie et al. | A beam segmentation based 2-step polar format algorithm for highly squinted SAR | |
Bączyk et al. | An analysis of ISAR image distortions in multistatic passive radars for air surveillance | |
Mao et al. | Signal processing for multi-rotors UAV SAR |
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 |