CN115856886A - Cross-domain distributed MISO-ISAR radar forward-looking imaging detection method - Google Patents

Cross-domain distributed MISO-ISAR radar forward-looking imaging detection method Download PDF

Info

Publication number
CN115856886A
CN115856886A CN202211528970.XA CN202211528970A CN115856886A CN 115856886 A CN115856886 A CN 115856886A CN 202211528970 A CN202211528970 A CN 202211528970A CN 115856886 A CN115856886 A CN 115856886A
Authority
CN
China
Prior art keywords
aperture
fusion
sub
echoes
isar
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.)
Pending
Application number
CN202211528970.XA
Other languages
Chinese (zh)
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.)
CASIC Microelectronic System Research Institute Co Ltd
Original Assignee
CASIC Microelectronic System Research Institute Co Ltd
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 CASIC Microelectronic System Research Institute Co Ltd filed Critical CASIC Microelectronic System Research Institute Co Ltd
Priority to CN202211528970.XA priority Critical patent/CN115856886A/en
Publication of CN115856886A publication Critical patent/CN115856886A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention relates to a cross-domain distributed MISO-ISAR radar forward-looking imaging detection method, which comprises the following steps: a plurality of movable transmitting nodes are arranged, and each transmitting node transmits orthogonal waveforms to form a transmitting array; the transmitting array moves according to the requirement, and a plurality of transmitting nodes transmit echoes; and setting a receiving station, receiving echoes by the receiving station, recovering a plurality of groups of two-dimensional imaging data sets, representing an observation sub-aperture by each group of echoes, and then carrying out multi-sub-aperture echo coherent fusion processing to obtain a fused optimal aperture image. The invention innovatively provides a multi-sub-aperture echo coherent fusion processing method, which is used for preprocessing on the basis of the MISO-ISAR technology, obtains the optimal features through fusion processing of image features, thereby obtaining the optimal aperture image, achieving the effect of high-definition two-dimensional imaging of the target, and realizing high-definition imaging resolution of the low-slow small target and the dense target in a detection airspace.

Description

Cross-domain distributed MISO-ISAR radar forward-looking imaging detection method
Technical Field
The invention belongs to the technical field of radar detection imaging, and particularly relates to a cross-domain distributed MISO-ISAR radar forward-looking imaging detection method.
Background
The cross-domain distributed MISO-ISAR radar forward-looking imaging detection technology is a radar detection technology which collaboratively constructs an air-ground combined MISO-ISAR (multiple-transmitting single-receiving-inverse synthetic aperture) imaging configuration through an air distributed type irradiation platform and a ground silent receiving imaging platform in time and space, performs air multi-platform distributed type transmission and ground platform silent receiving, and performs ISAR two-dimensional high-resolution imaging on a target.
A time-varying distributed MISO-ISAR imaging system is adopted, a distributed array is assumed to be motionless, the radar is equivalent to a fixed station-distributed MISO-ISAR radar, transmitting nodes simultaneously transmit orthogonal radar waveforms, a B-ISAR system is formed between each transmitting node and each receiving node, a plurality of B-ISAR systems form a plurality of equivalent phase centers in space, each phase center simultaneously observes a target from a plurality of angles, virtual apertures formed by rotation between each equivalent phase center and the target can be spliced in a synthetic aperture time range, and then the distributed MISO-ISAR can obtain a synthetic aperture larger than the B-ISAR, so that the higher direction resolution capability of the target can be realized.
The existing research of the MISO-ISAR synthetic aperture, such as the research content in the literature, "moving target distributed radar imaging technology research", only relates to the synthesis between two apertures, and the research on the synthesis method between a plurality of apertures is still lacked.
Disclosure of Invention
The invention provides a forward-looking imaging detection method for a cross-domain distributed MISO-ISAR radar, which realizes multi-aperture feature fusion on the basis of distributed MISO-ISAR.
The technical scheme adopted by the invention is as follows: a cross-domain distributed MISO-ISAR radar forward-looking imaging detection method comprises the following steps:
a plurality of movable transmitting nodes are arranged, and each transmitting node transmits orthogonal waveforms to form a transmitting array;
the transmitting array moves according to requirements, and a plurality of transmitting nodes transmit echoes;
and setting a receiving station, receiving echoes by the receiving station, recovering a plurality of groups of two-dimensional imaging data sets, wherein each group of echoes represents an observation sub-aperture, and then carrying out multi-sub-aperture echo coherent fusion processing to obtain a fused optimal aperture image (namely the fused optimal two-dimensional imaging data set).
Further, the multi-sub-aperture echo coherent fusion processing method comprises:
respectively generating basic characteristic functions K according to the imaged characteristics m (x),
K m (x)=f(x 1 ,x 2 ...x n ),
Wherein x is 1 ,x 2 ...x n Respectively corresponding to different characteristics of the imaging process;
obtaining a comprehensive characteristic function according to the basic characteristic function:
Figure BDA0003973885330000021
with the constraint of
Figure BDA0003973885330000022
Wherein, K m (x) Is a basic feature function, M represents the mth basic feature function, M is the number of the basic feature functionsNumber, d m A weight representing the corresponding underlying feature function; according to d m Solving a comprehensive characteristic function; and performing feature fusion according to the comprehensive feature function, and solving each optimal fusion feature to realize optimal aperture image fusion.
Further, the method for solving each optimal fusion feature includes:
the synthetic feature function K (x) is expressed as follows:
Figure BDA0003973885330000031
the optimal fusion algorithm is equivalently found to be an optimization problem, to
Figure BDA0003973885330000032
Finding the minimum value, and d corresponding to the minimum value m Substituting the comprehensive characteristic function to obtain an optimal comprehensive characteristic function K y (x);
Obtaining
Figure BDA0003973885330000033
The minimum constraints are as follows:
Figure BDA0003973885330000034
wherein ξ i Is a slack variable, C is a penalty factor, y i Are category labels. Xi i Is a relaxation variable, C is a penalty factor, y i The class target is given by performing multiple optimization assignments through a large number of samples and examples and combining engineering experience.
Further, the different features of the imaging process include: pixels, resolution, single synthetic aperture size, separation distance of apertures for the imaging process.
Further, the method also comprises the following step of fusing the sub-aperture:
and calculating a virtual aperture delta theta formed after multi-aperture fusion, and selecting a fusion method according to the virtual aperture delta theta to fuse the sub-apertures.
Further, the fusing the sub-apertures comprises: a total of P groups of echoes when
Figure BDA0003973885330000041
When the temperature of the water is higher than the set temperature,
the fusion is performed using a conventional sub-aperture fusion method, where θ p The aperture representing the p-th set of echoes,
Figure BDA0003973885330000042
β 0 representing the double base angle of the radar at time Tp = 0.
Further, the fusing the sub-apertures comprises: a total of P groups of echoes when
Figure BDA0003973885330000043
Time of flight
The following fusion model was used:
s=Aσ+e
wherein the content of the first and second substances,
Figure BDA0003973885330000044
s is a total echo signal formed by combining a plurality of echo signals S, A is an observation matrix corresponding to S, and the observation matrix A is determined by the rotating speed of the target relative to the P radars and the double base angle of the P-1 radars when Tp = 0; sigma is a single radar echo signal; e is the noise vector.
Compared with the prior art, the invention has the beneficial effects that:
the invention innovatively provides a multi-sub-aperture echo coherent fusion processing method, which is used for preprocessing on the basis of the MISO-ISAR technology, obtains the optimal features through fusion processing of image features, thereby obtaining the optimal aperture image, achieving the effect of high-definition two-dimensional imaging of the target, and realizing high-definition imaging resolution of the low-slow small target and the dense target in a detection airspace.
Drawings
Fig. 1 is a schematic diagram of multi-subaperture synthetic aperture fusion.
Detailed Description
The invention is further illustrated below with reference to the figures and examples.
The invention discloses a cross-domain distributed MISO-ISAR radar forward-looking imaging detection method, which comprises the following steps:
arranging a plurality of unmanned aerial vehicles in the air, wherein each unmanned aerial vehicle is a transmitting node, each transmitting node transmits orthogonal waveforms, and the unmanned aerial vehicles in the air form a transmitting array;
controlling the transmitting array to move according to requirements, and transmitting echoes by a plurality of transmitting nodes; and arranging a receiving station on the ground, receiving echoes by the receiving station, recovering a plurality of groups of two-dimensional imaging data sets, wherein each group of echoes represents an observation sub-aperture, and then carrying out multi-sub-aperture echo coherent fusion processing to obtain a fused optimal aperture image.
It is well known that the doppler resolution of moving object imaging is inversely proportional to the size of the synthetic aperture, with larger synthetic apertures giving higher resolution. In order to realize high-resolution imaging of a 'cluster low-slow small' target, an ultra-large synthetic aperture array must be synthesized, in ISAR imaging, the moving speed of the target, particularly the rotating angle of a relative radar in a synthetic aperture process, directly influences the size of a synthetic aperture, for the slow target, only single-base ISAR imaging is adopted, the slow movement of the target within the synthetic aperture time forms a very small rotating angle, and a large synthetic aperture cannot be formed for the small target.
The imaging performance and the realizability and the usability of the system are comprehensively considered, a time-varying distributed MISO-ISAR imaging system is provided, a plurality of unmanned aerial vehicles are arranged in the air, each unmanned aerial vehicle has the capability of transmitting and receiving echoes, each transmitting node transmits an orthogonal waveform, the aerial unmanned aerial vehicles form a transmitting array, and the transmitting array can be flexibly arranged and moved based on the system detection requirement; the ground vehicle is a receiving station, only receives and does not transmit, keeps the radio silent, and has high anti-detective and anti-interference capability.
A time-varying distributed MISO-ISAR imaging system is adopted, a distributed array is assumed to be motionless, the radar is equivalent to a fixed station-distributed MISO-ISAR radar, transmitting nodes simultaneously transmit orthogonal radar waveforms, a B-ISAR system is formed between each transmitting node and each receiving node, a plurality of B-ISAR systems form a plurality of equivalent phase centers in space, each phase center simultaneously observes a target from a plurality of angles, virtual apertures formed by rotation between each equivalent phase center and the target can be spliced in a synthetic aperture time range, and then the distributed MISO-ISAR can obtain a synthetic aperture larger than the B-ISAR, so that the higher direction resolution capability of the target can be realized.
Due to the small and slow characteristics of the target, particularly for small and slow unmanned aerial vehicles with the flight speed below subsonic speed and the target size not exceeding one meter, the virtual synthetic aperture formed by the distributed MISO-ISAR is not enough to form a high-resolution image for subsequent classification and identification on the target, so that on the basis of distributed MISO-ISAR imaging, the distributed transmitting nodes are enabled to move as required, and the accumulation rotation angle of azimuth high-resolution imaging is provided.
In the synthetic aperture time range, a plurality of nodes transmit echoes, in a ground receiving station, after orthogonal demodulation, a plurality of groups of two-dimensional echo data sets can be recovered, each group of echoes represents an observation sub-aperture, and the echo signals are obtained from equivalent phase centers of different observation angles in space, so the sub-aperture echoes are positioned on different azimuth angles. If only these data are imaged separately, the imaging resolution is low because the synthetic aperture is short.
In order to achieve high-resolution imaging of clustered small targets, the coherent fusion method of multiple echo arrays obtained by MISO-ISAR must be studied, while existing studies only discuss synthesis between two apertures, and synthesis methods between multiple apertures are still lacking. Therefore, the invention provides a multi-sub-aperture echo coherent fusion processing method, which is used for fusing imaging images of multiple echoes to obtain an optimal image.
The multi-sub-aperture echo coherent fusion processing method comprises the following steps:
respectively generating basic characteristic functions K according to the imaged characteristics m (x),
K m (x)=f(x 1 ,x 2 ...x n ),
Wherein x is 1 ,x 2 ...x n Each corresponding to a different characteristic of the imaging process including, but not limited to, the pixel, resolution, single synthetic aperture size, separation distance of apertures, etc. of the imaging process.
Obtaining a comprehensive characteristic function according to the basic characteristic function:
Figure BDA0003973885330000071
with the constraint of
Figure BDA0003973885330000072
Wherein, K m (x) As a basis feature function, M denotes the mth basis feature function, M is the number of basis feature functions, d m A weight representing the corresponding underlying feature function; according to d m Solving a comprehensive characteristic function; performing feature fusion according to the comprehensive feature function, and solving each optimal fusion feature to realize optimal aperture image fusion;
the synthetic feature function K (x) may be expressed as follows:
Figure BDA0003973885330000073
the optimal fusion algorithm is equivalently found to be an optimization problem, pair
Figure BDA0003973885330000074
Finding the minimum value, and d corresponding to the minimum value m Substituting the comprehensive characteristic function to obtain an optimal comprehensive characteristic function K y (x);
Obtaining
Figure BDA0003973885330000075
The minimum constraints are as follows:
Figure BDA0003973885330000081
wherein ξ i Is a relaxation variable, C is a penalty factor, y i Are category labels. Xi i 、C、y i And performing multiple optimization assignment giving through a large number of samples and examples and combining engineering experience.
The method has the advantages that the multiple images are fused through the algorithm, so that the fused optimal characteristics are obtained, the optimal aperture image is obtained, the effect of high-definition two-dimensional imaging of the target is achieved, and high-definition imaging resolution of the low-slow small target and the dense target in the detection space is realized.
Fusing the sub-aperture:
calculating the virtual aperture delta theta formed after multi-aperture fusion, and combining the sub-apertures delta theta as shown in FIG. 1 1 To delta theta P Obtaining delta theta after fusion, and selecting a fusion method according to the virtual aperture delta theta.
Fusing the sub-aperture specifically includes: assuming a common P set of echoes (i.e. P radars),
when the temperature is higher than the set temperature
Figure BDA0003973885330000082
When the temperature of the water is higher than the set temperature,
the fusion is performed using a conventional sub-aperture fusion method, where p represents the pth group of echoes, θ p The aperture representing the p-th set of echoes,
Figure BDA0003973885330000083
β 0 representing the double base angle of the radar at time Tp = 0. The conventional sub-aperture fusion method has been reported in the existing research, and reference can be specifically made to the third chapter of the article "research on distributed radar imaging technology for moving targets".
When in use
Figure BDA0003973885330000091
When the temperature of the water is higher than the set temperature,
the following fusion model was used: s = A σ + e
Wherein the content of the first and second substances,
Figure BDA0003973885330000092
s is the total echo signal obtained by combining a plurality of echo signals S, A is an observation matrix corresponding to S, and the observation matrix A is formed by an unknown number omega 12 ···ω P And beta 0203 ···β 0P It is determined that ω is the rotational speed of the target relative to the radar (P radars correspond to P ω), β 0203 ···β 0P Is the bistatic angle of the target to P-1 radars at the time Tp = 0; sigma is a single radar echo signal; e is the noise vector.

Claims (7)

1. A cross-domain distributed MISO-ISAR radar forward-looking imaging detection method is characterized by comprising the following steps:
a plurality of movable transmitting nodes are arranged, and each transmitting node transmits orthogonal waveforms to form a transmitting array;
the transmitting array moves according to requirements, and a plurality of transmitting nodes transmit echoes;
and setting a receiving station, receiving echoes by the receiving station, recovering a plurality of groups of two-dimensional imaging data sets, representing an observation sub-aperture by each group of echoes, and then carrying out multi-sub-aperture echo coherent fusion processing to obtain a fused optimal aperture image.
2. The cross-domain distributed MISO-ISAR radar forward-looking imaging detection method of claim 1, wherein the multi-sub-aperture echo coherent fusion processing method is as follows:
respectively generating basic characteristic functions K according to the imaged characteristics m (x),
K m (x)=f(x 1 ,x 2 ...x n ),
Wherein x is 1 ,x 2 ...x n Respectively corresponding to different characteristics of the imaging process;
obtaining a comprehensive characteristic function according to the basic characteristic function:
Figure FDA0003973885320000011
with the constraint of
Figure FDA0003973885320000012
Wherein, K m (x) As a basis feature function, M denotes the mth basis feature function, M is the number of basis feature functions, d m A weight representing the corresponding underlying feature function; according to d m Solving a comprehensive characteristic function; and performing feature fusion according to the comprehensive feature function, and solving each optimal fusion feature to realize optimal aperture image fusion.
3. The cross-domain distributed MISO-ISAR radar forward-looking imaging detection method of claim 2, wherein said method for finding each optimal fusion feature is:
the synthetic feature function K (x) is expressed as follows:
Figure FDA0003973885320000021
the optimal fusion algorithm is equivalently found to be an optimization problem, pair
Figure FDA0003973885320000022
Finding the minimum value, and d corresponding to the minimum value m Substituting the comprehensive characteristic function to obtain an optimal comprehensive characteristic function K y (x);
Obtaining
Figure FDA0003973885320000023
The minimum constraints are as follows: />
Figure FDA0003973885320000024
Wherein ξ i For the relaxation variable, C is a penalty factor and yi is a category label.
4. The cross-domain distributed MISO-ISAR radar forward-looking imaging detection method of claim 2, wherein the different features of the imaging process include: pixels, resolution, single synthetic aperture size, separation distance of apertures for the imaging process.
5. The cross-domain distributed MISO-ISAR radar forward-looking imaging detection method of claim 1, further comprising fusing sub-apertures:
and firstly, calculating a virtual aperture delta theta formed after the multi-aperture fusion, and selecting a fusion method delta theta according to the virtual aperture to fuse the sub-apertures.
6. The cross-domain distributed MISO-ISAR radar forward-looking imaging detection method of claim 5, wherein the fusing sub-apertures comprises: a total of P groups of echoes when
Figure FDA0003973885320000031
When the utility model is used, the water is discharged,
the fusion is performed using a conventional sub-aperture fusion method, where θ p An aperture representing the pth group of echoes;
Figure FDA0003973885320000032
β 0 representing the double base angle of the radar at Tp = 0.
7. The cross-domain distributed MISO-ISAR radar forward-looking imaging detection method of claim 5, wherein said fusing sub-apertures comprises: a total of P groups of echoes when
Figure FDA0003973885320000033
When the utility model is used, the water is discharged,
the following fusion model was used:
s=Aσ+e
wherein the content of the first and second substances,
Figure FDA0003973885320000034
s is a total echo signal formed by combining a plurality of echo signals, A is an observation matrix corresponding to S, the observation matrix A is determined by the rotating speed of a target relative to P radars and the bistatic angle of P-1 radars when Tp =0, and sigma is a single radar echo signal; e is the noise vector.
CN202211528970.XA 2022-11-30 2022-11-30 Cross-domain distributed MISO-ISAR radar forward-looking imaging detection method Pending CN115856886A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211528970.XA CN115856886A (en) 2022-11-30 2022-11-30 Cross-domain distributed MISO-ISAR radar forward-looking imaging detection method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211528970.XA CN115856886A (en) 2022-11-30 2022-11-30 Cross-domain distributed MISO-ISAR radar forward-looking imaging detection method

Publications (1)

Publication Number Publication Date
CN115856886A true CN115856886A (en) 2023-03-28

Family

ID=85668806

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211528970.XA Pending CN115856886A (en) 2022-11-30 2022-11-30 Cross-domain distributed MISO-ISAR radar forward-looking imaging detection method

Country Status (1)

Country Link
CN (1) CN115856886A (en)

Similar Documents

Publication Publication Date Title
Zhang et al. Multi-channel back-projection algorithm for mmwave automotive MIMO SAR imaging with Doppler-division multiplexing
US11378654B2 (en) Recurrent super-resolution radar for autonomous vehicles
JP6437925B2 (en) Multi-swath strip map SAR imaging
CN108828593B (en) Random radiation radar correlation imaging method
CN111352107B (en) Single pulse tracking and imaging method based on multi-channel digital sum and difference
KR101652395B1 (en) ISAR Cross-range scaling for a maneuvering target
CN103954964B (en) The method that multi-angle data of synthetic aperture radar obtains
CN112835009B (en) Data set expanding method based on RAI (radar-ranging) of FMCW (frequency modulated continuous wave)
CN109814100A (en) SAR Ground moving target imaging method based on sub-aperture parameter Estimation
US20220044359A1 (en) Super-resolution radar for autonomous vehicles
Stolz et al. A new antenna array and signal processing concept for an automotive 4D radar
CN110441772B (en) Satellite-borne sparse re-navigation SAR three-dimensional imaging method under large orbit distribution range
CN110261852B (en) Frequency diversity inverse synthetic aperture radar imaging system and method
CN115856886A (en) Cross-domain distributed MISO-ISAR radar forward-looking imaging detection method
CN116359921A (en) Quick time domain imaging method based on acceleration track double-base forward looking synthetic aperture radar
CN114185047B (en) Double-base SAR moving target refocusing method based on optimal polar coordinate transformation
CN114422060B (en) Method and system for constructing unmanned aerial vehicle communication channel model
WO2023070267A1 (en) Electronic device and related device
CN116148891A (en) Satellite-ground double-base sea surface two-dimensional flow field measurement method based on along-track interference
CN114966684A (en) Fast decomposition projection imaging method for group multi-base synthetic aperture radar
CN116184343A (en) Three-dimensional space swarm target detection and information estimation method based on phased array radar
Schlichenmaier et al. Clustering and subsequent contour and motion estimation of automotive objects using a network of cooperative radar sensors
Heidenreich et al. Introduction to the issue on recent advances in automotive radar signal processing
Kim SAR video generation of MIMO video SAR with beat frequency division FMCW
CN111522010A (en) Method and system for processing automobile anti-collision radar signals

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