CN109975804B - Multi-platform constellation SAR fusion coherent imaging method - Google Patents

Multi-platform constellation SAR fusion coherent imaging method Download PDF

Info

Publication number
CN109975804B
CN109975804B CN201910161386.7A CN201910161386A CN109975804B CN 109975804 B CN109975804 B CN 109975804B CN 201910161386 A CN201910161386 A CN 201910161386A CN 109975804 B CN109975804 B CN 109975804B
Authority
CN
China
Prior art keywords
fuzzy
azimuth
platform
channel
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.)
Active
Application number
CN201910161386.7A
Other languages
Chinese (zh)
Other versions
CN109975804A (en
Inventor
赵曜
黄永伟
陈如辉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guangdong University of Technology
Original Assignee
Guangdong University of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Guangdong University of Technology filed Critical Guangdong University of Technology
Priority to CN201910161386.7A priority Critical patent/CN109975804B/en
Publication of CN109975804A publication Critical patent/CN109975804A/en
Application granted granted Critical
Publication of CN109975804B publication Critical patent/CN109975804B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/904SAR modes
    • G01S13/9058Bistatic or multistatic SAR
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/9094Theoretical aspects

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 invention provides a multi-platform constellation SAR fusion coherent imaging method, which comprises the following steps: determining a noise covariance matrix between each platform according to echo data received by a plurality of platforms; constructing a multi-platform azimuth fuzzy suppression filter, and converting multi-channel fuzzy data into single-channel azimuth non-fuzzy data; and imaging the single-channel non-fuzzy data by adopting a compressed sensing method to obtain a multi-platform constellation SAR azimuth fuzzy suppression imaging result. According to the multi-platform constellation SAR fusion coherent imaging method provided by the invention, a mode of weighting a plurality of platform data is adopted, azimuth ambiguity is effectively inhibited, the problem of azimuth ambiguity in the process of realizing wide swath imaging of a single platform is solved, a compressed sensing technology is adopted, single-channel azimuth unambiguous data is imaged, side lobes can be effectively inhibited, and the imaging quality is ensured.

Description

Multi-platform constellation SAR fusion coherent imaging method
Technical Field
The invention relates to the technical field of radar imaging, in particular to a multi-platform constellation SAR fusion coherent imaging method.
Background
The azimuth geometric resolution and the range mapping bandwidth are two important technical indexes for measuring the SAR system. Due to the nyquist sampling theorem, it is difficult to simultaneously realize high resolution in the azimuth direction and wide mapping band in the distance direction, and with the rapid development of microsatellite technology in recent years, the multi-platform constellation SAR imaging method is also an important technical approach for realizing the high resolution and wide amplitude capability of SAR in the future. However, in a single platform, in order to realize a wide swath, a lower PRF is adopted, so that the azimuth direction has a certain azimuth ambiguity.
Disclosure of Invention
The invention provides a multi-platform constellation SAR fusion coherent imaging method, aiming at overcoming the technical defect that azimuth direction has certain azimuth ambiguity in the process of realizing wide swath imaging by using the existing single platform.
In order to solve the technical problems, the technical scheme of the invention is as follows:
a multi-platform constellation SAR fusion coherent imaging method comprises the following steps:
s1: determining a noise covariance matrix between each platform according to echo data received by a plurality of platforms;
s2: constructing a multi-platform azimuth fuzzy suppression filter, and converting multi-channel fuzzy data into single-channel azimuth non-fuzzy data;
s3: and imaging the single-channel non-fuzzy data by adopting a compressed sensing method to obtain a multi-platform constellation SAR azimuth fuzzy suppression imaging result.
Wherein, the step S1 specifically comprises:
for two-dimensional SAR echo data, x (f) is determined in the range-Doppler domain a ) The autocorrelation matrix of the statistical sample of (c), let x (f) a ) The statistical sample covariance matrix of (a) is:
Figure GDA0002072325730000011
wherein, x (f) a ) Signal x (f) indicating the presence of ambiguity for N platforms a )=[x 1 (f a ),x 2 (f a ),…,x N (f a )] H ,[] H Transpose of the representation vector, f a Indicating the azimuthal Doppler frequency, x (n) indicates the pair x (f) a ) Setting the Doppler fuzzy signal of N channels as x in the nth sampling under different range gates 1 (f a ),x 2 (f a ),…,x N (f a )。
Wherein, the step S2 specifically includes the following steps:
s21: constructing a multi-platform azimuth fuzzy suppression filter, wherein the specific form of the filter is as follows:
Figure GDA0002072325730000021
wherein, z (phi) i )=[z 1i ),...,z Ni )] T Representing the corresponding spatial steering vector of the signal,
Figure GDA0002072325730000022
φ i indicating the angle of orientation, Δ x, of the different platforms k Representing the distance difference between the kth platform and the phase origin, and mu represents a regularization parameter;
s22: converting multi-channel fuzzy data into single-channel azimuth non-fuzzy data, and setting weighted output single channelThe channel unambiguous signal is y (f) a ) Then, there are:
y(f a )=w H ·x(f a );
therefore, multi-channel fuzzy data is converted into single-channel azimuth unambiguous data.
Wherein, the step S3 specifically includes the following steps:
s31: constructing a mapping relation between the backscattering coefficient of the target and single-channel azimuth unambiguous data:
y=Φx+N;
wherein x represents the backscattering coefficient of the target; y represents single-channel azimuthal unambiguous data; phi represents a radar observation matrix; n represents system noise;
s32: optimizing the backscattering coefficient x of the target, wherein the specific calculation formula is as follows:
Figure GDA0002072325730000023
wherein λ is a regularization parameter;
s33: and obtaining a multi-platform constellation SAR azimuth fuzzy suppression imaging result according to the optimized backscattering coefficient x.
Compared with the prior art, the technical scheme of the invention has the beneficial effects that:
according to the multi-platform constellation SAR fusion coherent imaging method provided by the invention, a mode of weighting a plurality of platform data is adopted, azimuth ambiguity is effectively inhibited, the problem of azimuth ambiguity in the process of realizing wide swath imaging of a single platform is solved, a compressed sensing technology is adopted, single-channel azimuth unambiguous data is imaged, side lobes can be effectively inhibited, and the imaging quality is ensured.
Drawings
Fig. 1 is a flow diagram of a multi-platform constellation SAR fusion coherent imaging method.
Detailed Description
The drawings are for illustrative purposes only and are not to be construed as limiting the patent;
for the purpose of better illustrating the present embodiments, certain elements of the drawings may be omitted, enlarged or reduced, and do not represent the size of an actual product;
it will be understood by those skilled in the art that certain well-known structures in the drawings and descriptions thereof may be omitted.
The technical solution of the present invention is further described below with reference to the accompanying drawings and examples.
Example 1
As shown in fig. 1, a multi-platform constellation SAR fusion coherent imaging method includes the following steps:
s1: determining a noise covariance matrix between each platform according to echo data received by a plurality of platforms;
s2: constructing a multi-platform azimuth fuzzy suppression filter, and converting multi-channel fuzzy data into single-channel azimuth non-fuzzy data;
s3: and imaging the single-channel non-fuzzy data by adopting a compressed sensing method to obtain a multi-platform constellation SAR azimuth fuzzy suppression imaging result.
More specifically, the step S1 specifically includes:
for two-dimensional SAR echo data, x (f) is determined in the range-Doppler domain a ) The statistical sample autocorrelation matrix of (c), let x (f) a ) The statistical sample covariance matrix of (a) is:
Figure GDA0002072325730000031
wherein, x (f) a ) Signal x (f) indicating the presence of ambiguity for N platforms a )=[x 1 (f a ),x 2 (f a ),…,x N (f a )] H ,[] H Representing the transpose of the vector, f a Indicating the azimuthal Doppler frequency, x (n) indicates the pair x (f) a ) Setting the Doppler fuzzy signal of N channels as x in the nth sampling under different range gates 1 (f a ),x 2 (f a ),…,x N (f a )。
More specifically, the step S2 specifically includes the following steps:
s21: constructing a multi-platform azimuth fuzzy suppression filter, wherein the specific form of the filter is as follows:
Figure GDA0002072325730000032
wherein, z (phi) i )=[z 1i ),...,z Ni )] T Representing the corresponding spatial steering vector of the signal,
Figure GDA0002072325730000041
φ i indicating the angle of orientation, Δ x, of the different platforms k Representing the distance difference between the kth platform and the phase origin, and mu represents a regularization parameter;
s22: converting multi-channel fuzzy data into single-channel azimuth non-fuzzy data, and setting the weighted output single-channel non-fuzzy signal as y (f) a ) Then, there are:
y(f a )=w H ·x(f a );
therefore, multi-channel fuzzy data is converted into single-channel azimuth non-fuzzy data.
More specifically, the step S3 specifically includes the following steps:
s31: constructing a mapping relation between the backscattering coefficient of the target and single-channel azimuth unambiguous data:
y=Φx+N;
wherein x represents the backscattering coefficient of the target; y represents single-channel azimuth unambiguous data; phi represents a radar observation matrix; n represents system noise;
s32: optimizing the backscattering coefficient x of the target, wherein the specific calculation formula is as follows:
Figure GDA0002072325730000042
wherein λ is a regularization parameter;
s33: and obtaining a multi-platform constellation SAR azimuth fuzzy suppression imaging result according to the optimized backscattering coefficient x.
In a specific implementation process, the multi-platform constellation SAR fusion coherent imaging method provided by the invention adopts a mode of weighting a plurality of platform data, effectively inhibits azimuth ambiguity, avoids the problem of azimuth ambiguity in the process of realizing wide swath imaging by a single platform, adopts a compressed sensing technology to image single-channel azimuth unambiguous data, and can effectively inhibit side lobes to ensure imaging quality.
It should be understood that the above-described embodiments of the present invention are merely examples for clearly illustrating the present invention, and are not intended to limit the embodiments of the present invention. Other variations and modifications will be apparent to persons skilled in the art in light of the above description. And are neither required nor exhaustive of all embodiments. Any modification, equivalent replacement, and improvement made within the spirit and principle of the present invention should be included in the protection scope of the claims of the present invention.

Claims (2)

1. A multi-platform constellation SAR fusion coherent imaging method is characterized by comprising the following steps:
s1: determining a noise covariance matrix among the platforms according to echo data received by the platforms;
the step S1 specifically comprises the following steps:
for two-dimensional SAR echo data, x (f) is determined in the range-Doppler domain a ) The statistical sample autocorrelation matrix of (c), let x (f) a ) The statistical sample covariance matrix of (a) is:
Figure FDA0004053536060000011
wherein, x (f) a ) Signal x (f) indicating the presence of ambiguity for N platforms a )=[x 1 (f a ),x 2 (f a ),…,x N (f a )] H ,[] H Representing the transpose of the vector, f a Indicating the azimuthal Doppler frequency, x (n) indicates the pair x (f) a ) Setting the Doppler fuzzy signal of N channels as x in the nth sampling under different range gates 1 (f a ),x 2 (f a ),…,x N (f a );
S2: constructing a multi-platform azimuth fuzzy suppression filter, and converting multi-channel fuzzy data into single-channel azimuth non-fuzzy data;
the step S2 specifically includes the following steps:
s21: constructing a multi-platform azimuth fuzzy suppression filter, wherein the specific form of the filter is as follows:
Figure FDA0004053536060000012
wherein, z (phi) i )=[z 1i ),...,z Ni )] T Representing the corresponding spatial steering vector of the signal,
Figure FDA0004053536060000013
φ i indicating the angle of orientation, Δ x, of the different platforms k Representing the distance difference between the kth platform and the phase origin, and mu represents a regularization parameter;
s22: converting multi-channel fuzzy data into single-channel azimuth direction non-fuzzy data, and setting the weighted output single-channel non-fuzzy signal as y (f) a ) Then, there are:
y(f a )=w H ·x(f a );
so as to convert the multi-channel fuzzy data into single-channel azimuth non-fuzzy data;
s3: and imaging the single-channel non-fuzzy data by adopting a compressed sensing method to obtain a multi-platform constellation SAR azimuth fuzzy suppression imaging result.
2. The multi-platform constellation SAR fusion coherent imaging method according to claim 1, wherein said step S3 specifically comprises the steps of:
s31: constructing a mapping relation between the backscattering coefficient of the target and single-channel azimuth unambiguous data:
y=Φx+N;
wherein x represents the backscattering coefficient of the target; y represents single-channel azimuth unambiguous data; phi represents a radar observation matrix; n represents system noise;
s32: optimizing the backscattering coefficient x of the target, wherein the specific calculation formula is as follows:
Figure FDA0004053536060000021
wherein λ is a regularization parameter;
s33: and obtaining a multi-platform constellation SAR azimuth fuzzy suppression imaging result according to the optimized backscattering coefficient x.
CN201910161386.7A 2019-03-04 2019-03-04 Multi-platform constellation SAR fusion coherent imaging method Active CN109975804B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910161386.7A CN109975804B (en) 2019-03-04 2019-03-04 Multi-platform constellation SAR fusion coherent imaging method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910161386.7A CN109975804B (en) 2019-03-04 2019-03-04 Multi-platform constellation SAR fusion coherent imaging method

Publications (2)

Publication Number Publication Date
CN109975804A CN109975804A (en) 2019-07-05
CN109975804B true CN109975804B (en) 2023-04-07

Family

ID=67077824

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910161386.7A Active CN109975804B (en) 2019-03-04 2019-03-04 Multi-platform constellation SAR fusion coherent imaging method

Country Status (1)

Country Link
CN (1) CN109975804B (en)

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8493262B2 (en) * 2011-02-11 2013-07-23 Mitsubishi Electric Research Laboratories, Inc. Synthetic aperture radar image formation system and method
CN103630897A (en) * 2012-08-28 2014-03-12 中国科学院电子学研究所 Multichannel synthetic aperture radar imaging method
JP6000074B2 (en) * 2012-11-08 2016-09-28 三菱電機株式会社 Image radar processing apparatus and image radar processing method
CN103576151B (en) * 2013-10-15 2015-07-15 中国科学院电子学研究所 Azimuth multi-channel SAR imaging method and system based on compressed sensing
CN103698755B (en) * 2014-01-16 2016-03-09 中国科学院电子学研究所 A kind of method obtaining sparse microwave imaging radar signal to noise ratio (S/N ratio) and sampling ratio
CN103760559A (en) * 2014-01-24 2014-04-30 中国科学院电子学研究所 Sparse microwave imaging method based on online observation
CN104345310A (en) * 2014-10-21 2015-02-11 中国科学院电子学研究所 Method and device for realizing imaging of synthetic aperture radar
CN104865571A (en) * 2015-04-02 2015-08-26 西安电子科技大学 Multi-channel multi-sub-band sliding-spotlight-mode SAR imaging method
CN107229048B (en) * 2017-06-06 2020-07-31 电子科技大学 High-resolution wide-range SAR moving target speed estimation and imaging method
CN107918126B (en) * 2017-11-17 2021-05-18 西安电子科技大学 Multi-feature automatic segmentation-based multi-channel near-coast fuzzy clutter suppression method
CN108845318B (en) * 2018-05-31 2022-04-15 南京航空航天大学 Satellite-borne high-resolution wide-range imaging method based on Relax algorithm
CN109001733B (en) * 2018-06-15 2020-04-28 电子科技大学 GEO satellite-machine bistatic synthetic aperture radar spectrum fuzzy suppression method

Also Published As

Publication number Publication date
CN109975804A (en) 2019-07-05

Similar Documents

Publication Publication Date Title
CN109581352B (en) Super-resolution angle measurement system based on millimeter wave radar
Zhang et al. Adaptive two-step calibration for high-resolution and wide-swath SAR imaging
CN105842693B (en) A kind of method of the Dual-Channel SAR moving-target detection based on compressed sensing
CN111142105A (en) ISAR imaging method for complex moving target
CN110488293B (en) Distributed SAR system with non-uniform spatial configuration
CN112444811A (en) Target detection and imaging method integrating MIMO radar and ISAR
JP2019519784A (en) SAR imaging method for interference analysis
CN109212489B (en) FDA-MIMO radar fuzzy clutter suppression method based on auxiliary pulse
CN111220955B (en) Airborne weather radar ground clutter suppression method based on vertical array covariance matrix characteristic decomposition
CN111175749B (en) Satellite-borne SAR imaging processing method
CN105301589B (en) High-resolution Wide swath SAR Ground moving target imaging method
CN110865346B (en) Satellite-borne SAR time parameter calibration method based on direct positioning algorithm
CN109901162B (en) Long-baseline imaging STAP method suitable for distributed geosynchronous orbit SAR
CN114545411B (en) Polar coordinate format multimode high-resolution SAR imaging method based on engineering realization
JP5634238B2 (en) Radar equipment
CN110488283B (en) Error correction method for multi-channel HRWS-SAR channel
CN116400310B (en) Two-dimensional frequency domain azimuth multi-channel SAR error correction method
CN113126087B (en) Space-borne interference imaging altimeter antenna
CN111722225B (en) Bistatic SAR two-dimensional self-focusing method based on prior phase structure information
CN110133646B (en) NLCS imaging-based multi-channel two-pulse clutter cancellation method for bistatic forward-looking SAR
CN113933804B (en) Image domain azimuth multi-channel SAR error correction method
CN113466797A (en) Bistatic SAR space-time clutter suppression method based on clutter ridge matching sparse recovery
CN109975804B (en) Multi-platform constellation SAR fusion coherent imaging method
CN110658502A (en) Amplitude-phase error correction method
Chen et al. Performance Analysis of Wavenumber Domain Algorithms for Highly Squinted 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