WO2022085044A1 - Radar signal processing device and radar signal processing method - Google Patents
Radar signal processing device and radar signal processing method Download PDFInfo
- Publication number
- WO2022085044A1 WO2022085044A1 PCT/JP2020/039224 JP2020039224W WO2022085044A1 WO 2022085044 A1 WO2022085044 A1 WO 2022085044A1 JP 2020039224 W JP2020039224 W JP 2020039224W WO 2022085044 A1 WO2022085044 A1 WO 2022085044A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- unit
- channel signals
- imbalance correction
- channel
- band
- Prior art date
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
- 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
Definitions
- This disclosure relates to a radar signal processing device.
- the multi-channel synthetic aperture radar (hereinafter referred to as multi-channel SAR) is provided with a plurality of receiving antennas, so that a plurality of signals (hereinafter referred to as a plurality of channel signals) are received via a plurality of channels.
- the multi-channel synthetic aperture radar generates a SAR image by performing a restoration process based on a plurality of received channel signals.
- the multi-channel synthetic aperture radar corrects the amount of deviation between channel signals (hereinafter referred to as imbalance) caused by the difference in channels (for example, Non-Patent Document 1).
- the azimuth ambiguity is multiplexed, and the above-mentioned is performed with sufficient accuracy. There is a problem that it becomes difficult to correct the imbalance between the channel signals of.
- the present disclosure has been made to solve the above-mentioned problems, and an object of the present disclosure is to provide a technique for improving the accuracy of imbalance correction between channel signals.
- the hit region band limiting unit and the hit region band limiting unit that limit the bandwidth in the hit region hit a plurality of channel signals obtained by observation by the multi-channel SAR, respectively. It includes an imbalance correction unit that performs imbalance correction between channel signals based on a plurality of channel signals for which band limitation has been performed in a region.
- the accuracy of imbalance correction between channel signals can be improved.
- FIG. 4A is a block diagram showing a hardware configuration that realizes the functions of the radar signal processing device according to the first embodiment.
- FIG. 4B is a block diagram showing a hardware configuration for executing software that realizes the functions of the radar signal processing device according to the first embodiment.
- FIG. 1 is a block diagram showing a configuration of a radar signal processing device 100 according to the first embodiment.
- the radar signal processing device 100 includes a power imbalance correction unit 1, a fast Fourier transform unit 2, a Doppler frequency domain band limiting unit 3, an inverse fast Fourier transform unit 4, a hit region band limiting unit 5, and an in. It includes a balance correction unit 6, a sub-aperture division unit 9, a restoration algorithm unit 10, an image reproduction unit 11, and a sub-aperture coupling unit 12.
- the power imbalance correction unit 1 corrects the power imbalance for a plurality of channel signals (raw data) obtained by observation by the multi-channel SAR.
- the power imbalance correction unit 1 outputs a plurality of channel signals that have undergone power imbalance correction to the fast Fourier transform unit 2.
- the "plurality of channel signals” means signals obtained from the corresponding channels among the plurality of channels in the multi-channel SAR, respectively.
- power imbalance correction means to correct the amount of power deviation between channel signals caused by a channel difference.
- the observation by the multi-channel SAR is a sliding spot observation.
- the configuration in which the observation by the multi-channel SAR is the sliding spot observation will be described, but the observation by the multi-channel SAR may be, for example, a full spotlight observation or the like.
- the fast Fourier transform unit 2 performs Doppler frequency conversion for each of a plurality of channel signals obtained by observation by the multi-channel SAR.
- the fast Fourier transform unit 2 converts a plurality of channel signals obtained by observation by the multi-channel SAR into signals in the Doppler frequency domain by performing a fast Fourier transform. do.
- the fast Fourier transform unit 2 performs Doppler frequency conversion for each of the plurality of channel signals for which the power imbalance correction unit 1 has performed the power imbalance correction.
- the fast Fourier transform unit 2 outputs a plurality of channel signals that have undergone Doppler frequency conversion to the Doppler frequency domain band limiting unit 3.
- the Doppler frequency domain band limiting unit 3 limits the band in the Doppler frequency domain for each of the plurality of channel signals for which the fast Fourier transform unit 2 has undergone the Doppler frequency conversion. More specifically, in the first embodiment, the Doppler frequency domain band limiting unit 3 extracts a channel signal in a band centered on the Doppler center frequency by band limiting in the Doppler frequency domain. The Doppler frequency domain band limiting unit 3 outputs a plurality of channel signals band-limited in the Doppler frequency domain to the inverse fast Fourier transform unit 4.
- the inverse fast Fourier transform unit 4 (inverse Doppler frequency converter) performs inverse Doppler frequency conversion for each of a plurality of channel signals for which the Doppler frequency domain band limiting unit 3 has band-limited in the Doppler frequency domain.
- the inverse fast Fourier transform unit 4 performs an inverse fast Fourier transform on each of a plurality of channel signals for which the Doppler frequency domain band limiting unit 3 has band-limited in the Doppler frequency domain, thereby performing a fast Fourier transform.
- the transforming unit 2 converts the signal into the signal in the original time domain before performing the fast Fourier transform.
- the inverse fast Fourier transform unit 4 outputs a plurality of channel signals subjected to the inverse Doppler frequency conversion to the hit region band limiting unit 5.
- the hit area band limiting unit 5 limits the band in the hit area for each of the plurality of channel signals obtained by the observation by the multi-channel SAR. More specifically, in the first embodiment, the hit region band limiting unit 5 limits the bandwidth in the hit region for each of the plurality of channel signals subjected to the inverse Doppler frequency conversion by the inverse fast Fourier transform unit 4. .. More specifically, in the first embodiment, the hit region band limiting unit 5 extracts a channel signal in a band centered on the hit center by band limiting in the hit region. The hit area band limiting unit 5 outputs a plurality of channel signals for which the band is limited in the hit area to the imbalance correction unit 6.
- the imbalance correction unit 6 performs imbalance correction between channel signals based on a plurality of channel signals for which the hit region band limiting unit 5 has band-limited in the hit region. More specifically, in the first embodiment, the imbalance correction unit 6 corrects the phase imbalance as the imbalance correction between the channel signals. The imbalance correction unit 6 outputs a plurality of channel signals for which imbalance correction has been performed to the sub-aperture division unit 9.
- phase imbalance means the amount of phase shift between channel signals caused by the difference in channels.
- the imbalance correction unit 6 corrects the phase imbalance as the imbalance correction will be described, but the imbalance correction unit 6 corrects the range sample deviation as the imbalance correction between the channel signals. May be done.
- the imbalance correction unit 6 includes a Doppler center frequency estimation unit 7 and a phase imbalance correction unit 8.
- the Doppler center frequency estimation unit 7 estimates the Doppler center frequencies of a plurality of channel signals for which the hit region band limiting unit 5 has band-limited in the hit region.
- the Doppler center frequency estimation unit 7 outputs the estimated Doppler center frequency to the phase imbalance correction unit 8.
- the phase imbalance correction unit 8 is based on the Doppler center frequency estimated by the Doppler center frequency estimation unit 7, and is based on a plurality of channel signals for which the hit region band limiting unit 5 has band-limited in the hit region. The phase imbalance of is corrected.
- the phase imbalance correction unit 8 outputs a plurality of channel signals for which the phase imbalance has been corrected to the sub-aperture division unit 9.
- the sub-aperture division unit 9 performs sub-aperture division for a plurality of channel signals for which the imbalance correction unit 6 has performed imbalance correction. More specifically, in the first embodiment, the sub-aperture division unit 9 divides the sub-aperture into a plurality of channel signals for which the phase imbalance correction unit 8 of the imbalance correction unit 6 has corrected the phase imbalance. I do. The sub-aperture division unit 9 outputs a plurality of channel signals obtained by sub-aperture division to the restoration algorithm unit 10.
- the restoration algorithm unit 10 performs a restoration algorithm on a plurality of channel signals for which the sub-aperture division unit 9 has performed sub-aperture division.
- the restoration algorithm unit 10 outputs a plurality of channel signals subjected to the restoration algorithm to the image reproduction unit 11.
- the image reproduction unit 11 reproduces an image on a plurality of channel signals for which the restoration algorithm unit 10 has performed the restoration algorithm.
- the image reproduction unit 11 outputs a plurality of channel signals for which image reproduction has been performed to the sub-aperture coupling unit 12.
- the sub-aperture coupling unit 12 generates a SAR image by performing sub-aperture coupling on a plurality of channel signals for which the image reproduction unit 11 has reproduced an image. For example, the SAR image generated by the sub-aperture coupling portion 12 is displayed by a display (not shown).
- FIG. 2 is a flowchart showing a radar signal processing method by the radar signal processing device 100 according to the first embodiment. Before each step described below, it is assumed that the radar signal processing device 100 has acquired a plurality of channel signals obtained by observation by the multi-channel SAR.
- the power imbalance correction unit 1 performs power imbalance correction on a plurality of channel signals obtained by observation by the multi-channel SAR (step ST1).
- the power imbalance correction unit 1 outputs a plurality of channel signals that have undergone power imbalance correction to the fast Fourier transform unit 2.
- the fast Fourier transform unit 2 performs Doppler frequency conversion for each of the plurality of channel signals for which the power imbalance correction unit 1 has performed the power imbalance correction (step ST2).
- the fast Fourier transform unit 2 outputs a plurality of channel signals that have undergone Doppler frequency conversion to the Doppler frequency domain band limiting unit 3.
- the Doppler frequency domain band limiting unit 3 limits the band in the Doppler frequency domain for each of the plurality of channel signals subjected to the Doppler frequency conversion by the fast Fourier transform unit 2 (step ST3).
- the Doppler frequency domain band limiting unit 3 outputs a plurality of channel signals band-limited in the Doppler frequency domain to the inverse fast Fourier transform unit 4.
- the inverse fast Fourier transform unit 4 performs inverse Doppler frequency conversion for each of the plurality of channel signals for which the Doppler frequency domain band limiting unit 3 has band-limited in the Doppler frequency domain (step ST4).
- the inverse fast Fourier transform unit 4 outputs a plurality of channel signals subjected to the inverse Doppler frequency conversion to the hit region band limiting unit 5.
- the hit region band limiting unit 5 limits the bandwidth in the hit region for each of the plurality of channel signals for which the inverse fast Fourier transform unit 4 has undergone the inverse Doppler frequency conversion (step ST5).
- the hit area band limiting unit 5 outputs a plurality of channel signals for which the band is limited in the hit area to the imbalance correction unit 6.
- the Doppler center frequency estimation unit 7 of the imbalance correction unit 6 estimates the Doppler center frequency of a plurality of channel signals for which the hit region band limiting unit 5 has band-limited in the hit region (step ST6).
- the Doppler center frequency estimation unit 7 outputs the estimated Doppler center frequency to the phase imbalance correction unit 8.
- the phase imbalance correction unit 8 of the imbalance correction unit 6 corrects the phase imbalance between the channel signals based on the Doppler center frequency estimated by the Doppler center frequency estimation unit 7 (step ST7).
- the phase imbalance correction unit 8 outputs a plurality of channel signals for which the phase imbalance has been corrected to the sub-aperture division unit 9.
- the sub-aperture division unit 9 performs sub-aperture division for a plurality of channel signals for which the imbalance correction unit 6 has corrected the phase imbalance (step ST8).
- the sub-aperture division unit 9 outputs a plurality of channel signals obtained by sub-aperture division to the restoration algorithm unit 10.
- the restoration algorithm unit 10 performs a restoration algorithm on a plurality of channel signals for which the sub-aperture division unit 9 has undergone sub-aperture division (step ST9).
- the restoration algorithm unit 10 outputs a plurality of channel signals subjected to the restoration algorithm to the image reproduction unit 11.
- the image reproduction unit 11 reproduces an image on a plurality of channel signals for which the restoration algorithm unit 10 has performed the restoration algorithm (step ST10).
- the image reproduction unit 11 outputs a plurality of channel signals for which image reproduction has been performed to the sub-aperture coupling unit 12.
- the sub-aperture coupling unit 12 generates a SAR image by performing sub-aperture coupling on a plurality of channel signals for which the image reproduction unit 11 has reproduced an image (step ST11). For example, the SAR image generated by the sub-aperture coupling unit 12 in step ST11 is displayed by a display (not shown).
- FIG. 3 is a diagram showing the progress of a channel signal processed by the radar signal processing method by the radar signal processing device 100 according to the first embodiment.
- the plurality of channel signals obtained by the observation by the multi-channel SAR are the signals obtained by the two-channel sliding spotlight observation.
- the first channel signal (channel 0 signal) before being processed by the radar signal processing device 100 is s 0 ( ⁇ , ⁇ ), and the second channel signal (channel 1 signal) is s 1 ( ⁇ , ⁇ ). ).
- ⁇ is the range time and ⁇ is a hit.
- the first channel signal will be used as a reference channel.
- the power imbalance correction unit 1 has the following with respect to the first channel signal s 0 ( ⁇ , ⁇ ) and the second channel signal s 1 ( ⁇ , ⁇ ). Power imbalance correction is performed according to the method described in (For example, Non-Patent Document J. H. Kim, M. Younis, P. P. Iraola, M. Gebele, and G. Krieger, “First spaceborne demonstration of digital beamforming for azimuth ambiguity suppression”, IEEE. See Trans. Geosci. Remote Sens., Vol.51, no.1, pp.579-590, Jan., 2013).
- the power imbalance correction unit 1 uses the first channel signal s 0 ( ⁇ , ⁇ ) and the second channel signal s 1 ( ⁇ , ⁇ ) to correct the power imbalance in the range frequency region. On the other hand, range frequency conversion is performed. As a result, the power imbalance correction unit 1 obtains S 0 ( fr , ⁇ ) and S 1 ( fr , ⁇ ) with the range frequency as fr .
- the power imbalance correction unit 1 sets the noise power generated by the first channel signal as ⁇ 0 2 , and the noise power generated by the second channel signal as ⁇ 1 2 , according to the following equation (1).
- the correction coefficient ⁇ ( fr ) of the power (amplitude) imbalance is calculated.
- 2 ) indicates the average of
- 2 ) indicates the average of
- the power imbalance correction unit 1 corrects the amplitude of S 1 ( fr , ⁇ ) according to the following equation (2) based on the calculated correction coefficient ⁇ (fr).
- the power imbalance correction unit 1 corrects the power imbalance by performing range time conversion on S 1 '( fr , ⁇ ) obtained by the equation (2).
- the power imbalance correction unit 1 may use a method other than the method of the specific example as a method of correcting the power imbalance.
- FIG. 3A shows a graph of the first channel signal s 0 ( ⁇ , ⁇ ) after the processing of step ST1.
- the vertical axis in A indicates a hit
- the horizontal axis indicates a range.
- the raw data of the seven point images are multiplexed in the hit-range space.
- the fast Fourier transform unit 2 has a Doppler frequency with respect to the first channel signal s 0 ( ⁇ , ⁇ ) and the second channel signal s 1 '( ⁇ , ⁇ ), respectively.
- FIG. 3B shows a graph of the first channel signal after the processing of step ST2.
- the vertical axis in B indicates the Doppler frequency
- the horizontal axis indicates the range.
- the first channel signal s 0 ( ⁇ , ⁇ ) is multiplexed with azimuth ambiguity.
- the azimuth ambiguity is similarly multiplexed for the second channel signal s 1 '( ⁇ , ⁇ ).
- FIG. 3C shows a graph of the first channel signal after the processing of step ST3.
- the vertical axis in C in FIG. 3 indicates the Doppler frequency, and the horizontal axis indicates the range.
- a signal near the Doppler center frequency is extracted by the process of step ST3.
- a signal near the Doppler center frequency is similarly extracted by the process of step ST3.
- step ST4 the inverse fast Fourier transform unit 4 performs inverse Doppler frequency conversion on the first channel signal and the second channel signal obtained by the processing of step ST3.
- FIG. 3D shows a graph of the first channel signal after the processing of step ST4.
- the vertical axis in D in FIG. 3 indicates a hit, and the horizontal axis indicates a range.
- D in FIG. 3 in the first channel signal s 0 ( ⁇ , ⁇ ) after the processing of step ST4, the signals are discretely dispersed.
- the signals of the second channel signal s 1 '( ⁇ , ⁇ ) after the processing of step ST4 are also discretely dispersed.
- the hit region band limiting unit 5 has a hit region for each hit-range component of the first channel signal and the second channel signal obtained by the processing of step ST4.
- FIG. 3E shows a graph of the first channel signal after the processing of step ST5.
- the vertical axis in E in FIG. 3 indicates a hit, and the horizontal axis indicates a range.
- E in FIG. 3 by suppressing the azimuth ambiguity, a signal from which the azimuth ambiguity is separated is obtained.
- a signal from which the azimuth ambiguity is separated can be obtained similarly for the second channel signal s 1 '( ⁇ , ⁇ ).
- the first channel signal obtained by the process of step ST5 is s 0 ( ⁇ , ⁇ )
- the second channel signal is s 1 '( ⁇ , ⁇ ).
- the imbalance correction unit 6 performs the first channel signal s 0 ( ⁇ , ⁇ ) and the second channel signal s 1 '( ⁇ ) obtained by the processing of step ST5. , ⁇ ), and corrects the phase imbalance between the channel signals according to the method described below (eg, H. Fan, Z. Zhang, and R. Wang, “phase mismatch calibration for multichannel sliding spotlight SAR). imaging with extended azimuth cross correlation ”, see IGARSS 2019).
- step ST6 the Doppler center frequency estimation unit 7 of the imbalance correction unit 6 is based on the first channel signal s 0 ( ⁇ , ⁇ ) and the second channel signal s 1 '( ⁇ , ⁇ ). Then, two interference data C 1,0 ( ⁇ ) and C 0,1 ( ⁇ ) are calculated according to the following equations (3) and (4).
- the PRF is the pulse repetition frequency and corresponds to one sample shift in the azimuth direction in the channel signal.
- E ⁇ () indicates the average of the equations in parentheses in the hit direction, as described above.
- the Doppler center frequency estimation unit 7 calculates the product C 1,0 ( ⁇ ) C 0,1 ( ⁇ ) of both the interference data C 1,0 ( ⁇ ) and C 0,1 ( ⁇ ). Therefore, the phase imbalance contained in s 1 '( ⁇ , ⁇ ) is offset. Assuming that the Doppler center frequency is f dc ( ⁇ ), the relationship of the following equation (5) is obtained.
- the Doppler center frequency estimation unit 7 calculates the Doppler center frequency f dc ( ⁇ ) based on the calculated C 1,0 ( ⁇ ) C 0,1 ( ⁇ ) according to the following equation (6).
- the phase imbalance correction unit 8 of the imbalance correction unit 6 has the following equation (7) based on the Doppler center frequency f dc ( ⁇ ) calculated by the Doppler center frequency estimation unit 7.
- the phase imbalance ⁇ ( ⁇ ) is calculated according to the above.
- vs represents the platform speed
- d represents the distance between the reception phase centers of the two channels.
- the phase imbalance correction unit 8 generates a weight for the obtained phase imbalance ⁇ ( ⁇ ) and performs a weighted averaging process so as to take an expected value in the hit direction. This process reduces the error in the estimated value. More specifically, first, the phase imbalance correction unit 8 calculates the average of
- the phase imbalance correction unit 8 takes a weighted average according to the following equation (9) and calculates the phase imbalance ⁇ .
- E ⁇ ⁇ indicates the average of the equations in parentheses in the range time direction.
- the phase imbalance correction unit 8 corrects the phase imbalance between the channel signals according to the following equation (10) based on the calculated phase imbalance ⁇ . More specifically, in the specific example, the phase imbalance correction unit 8 has the second channel signal s after the power imbalance correction described above according to the following equation (10) based on the calculated phase imbalance ⁇ . Correct the phase imbalance for 1 '( ⁇ , ⁇ ).
- the imbalance correction unit 6 may use a method other than the method of the specific example as a method of correcting the phase imbalance.
- the sub-aperture dividing unit 9 has the second channel signal s 1 '' ( ⁇ , ⁇ ) and the first channel signal s 0 ( ⁇ , ⁇ ) after the processing of step ST7.
- the method of sub-aperture division for example, the non-patent document J. Mittermayer, R. Lord, and E. B ⁇ oner, “Sliding spotlight SAR processing for TerraSAR-X using a new” fomulation of the extended chirp scaling algorithm, ”IGARSS2003, vol.3, p.1462-1464, 2003.).
- the restoration algorithm unit 10 performs a restoration algorithm on the first channel signal and the second channel signal to which the sub-aperture division unit 9 has performed the sub-aperture division (reconstruction algorithm).
- a restoration algorithm for the method, see, for example, Non-Patent Document 1 described above).
- the image reproduction unit 11 reproduces an image with respect to the first channel signal and the second channel signal to which the restoration algorithm unit 10 has performed the restoration algorithm.
- the image reproduction method include extended chirp scaling (for example, the above-mentioned non-patent documents J. Mittermayer, R. Lord, and E. B ⁇ oner, “Sliding spotlight SAR processing for TerraSAR-X using a”. See new fomulation of the extended chirp scaling algorithm, ”IGARSS2003, vol.3, p.1462-1464, 2003.).
- the sub-aperture coupling unit 12 obtains a SAR image by performing sub-aperture coupling to the first channel signal and the second channel signal for which the image reproduction unit 11 has reproduced the image. Generate.
- the imbalance correction unit 6 only corrects the phase imbalance in the subsequent stage of the band limitation by the bandpass filter twice, but the imbalance correction unit 6 corrects the range sample deviation and the like. Other imbalance corrections may be made.
- the band path filters of the Doppler frequency domain band limiting unit 3 and the hit region band limiting unit 5 performed the band limitation twice, the Doppler frequency domain band limiting unit 3 did not limit the band of the Doppler frequency region.
- the hit area band limiting unit 5 may limit the band only in the hit area.
- Each function of the frequency estimation unit 7, the phase imbalance correction unit 8, the sub-aperture division unit 9, the restoration algorithm unit 10, the image reproduction unit 11, and the sub-aperture coupling unit 12 is realized by a processing circuit. That is, the radar signal processing device 100 includes a processing circuit for executing the processing of each step shown in FIG. This processing circuit may be dedicated hardware, or may be a CPU (Central Processing Unit) that executes a program stored in the memory.
- CPU Central Processing Unit
- FIG. 4A is a block diagram showing a hardware configuration that realizes the functions of the radar signal processing device 100.
- FIG. 4B is a block diagram showing a hardware configuration for executing software that realizes the functions of the radar signal processing device 100.
- the processing circuit 20 may be, for example, a single circuit, a composite circuit, a programmed processor, a parallel programmed processor, or an ASIC (Application Specific Integrated Circuitd). Circuit), FPGA (Field-Programmable Gate Array) or a combination thereof is applicable.
- ASIC Application Specific Integrated Circuitd
- FPGA Field-Programmable Gate Array
- Doppler center of power imbalance correction unit 1 Doppler center of power imbalance correction unit 1, fast Fourier transform unit 2, Doppler frequency domain band limiting unit 3, inverse fast Fourier transform unit 4, hit region band limiting unit 5, and imbalance correction unit 6 in the radar signal processing device 100.
- the functions of the frequency estimation unit 7, the phase imbalance correction unit 8, the sub-aperture division unit 9, the restoration algorithm unit 10, the image reproduction unit 11, and the sub-aperture coupling unit 12 may be realized by separate processing circuits, or these may be realized. The functions may be collectively realized by one processing circuit.
- the power imbalance correction unit 1 When the processing circuit is the processor 21 shown in FIG. 4B, the power imbalance correction unit 1, the fast Fourier transform unit 2, the Doppler frequency domain band limiting unit 3, the inverse fast Fourier transform unit 4, and the hit in the radar signal processing device 100.
- the software or firmware is described as a program and stored in the memory 22.
- the processor 21 By reading and executing the program stored in the memory 22, the processor 21 reads and executes the power imbalance correction unit 1, the fast Fourier transform unit 2, the Doppler frequency domain band limiting unit 3, and the inverse high-speed Fourier in the radar signal processing device 100. Conversion unit 4, hit region band limiting unit 5, Doppler center frequency estimation unit 7 and phase imbalance correction unit 8 of imbalance correction unit 6, sub-aperture division unit 9, restoration algorithm unit 10, image reproduction unit 11, and sub-aperture coupling. Each function of the part 12 is realized. That is, the radar signal processing device 100 includes a memory 22 for storing a program in which the processing of each step shown in FIG. 2 is executed as a result when each of these functions is executed by the processor 21.
- These programs include a power imbalance correction unit 1, a high-speed Fourier transform unit 2, a Doppler frequency domain band limiting unit 3, an inverse fast Fourier transform unit 4, a hit region band limiting unit 5, and an imbalance correction in the radar signal processing device 100.
- the computer is made to execute each procedure or method of the Doppler center frequency estimation unit 7 of the unit 6, the phase imbalance correction unit 8, the sub-aperture division unit 9, the restoration algorithm unit 10, the image reproduction unit 11, and the sub-aperture coupling unit 12.
- the memory 22 uses the computer as a power imbalance correction unit 1, a fast Fourier transform unit 2, a Doppler frequency domain band limiting unit 3, an inverse fast Fourier transform unit 4, a hit region band limiting unit 5, and an algorithm in the radar signal processing device 100.
- a program for functioning as the Doppler center frequency estimation unit 7 of the balance correction unit 6, the phase imbalance correction unit 8, the sub-aperture division unit 9, the restoration algorithm unit 10, the image reproduction unit 11, and the sub-aperture coupling unit 12 is stored. It may be a computer-readable storage medium.
- the processor 21 corresponds to, for example, a CPU (Central Processing Unit), a processing device, a computing device, a processor, a microprocessor, a microcomputer, a DSP (Digital Signal Processor), or the like.
- the memory 22 may include, for example, a RAM (Random Access Memory), a ROM (Read Only Memory), a flash memory, an EPROM (Erasable Programmable Read Only Memory), an EPROM (Electrically-volatile) semiconductor, or an EPROM (Electrically-EPROM).
- a RAM Random Access Memory
- ROM Read Only Memory
- flash memory an EPROM (Erasable Programmable Read Only Memory)
- EPROM Electrically-volatile semiconductor
- EPROM Electrically-EPROM
- Doppler center of power imbalance correction unit 1 Doppler center of power imbalance correction unit 1, fast Fourier transform unit 2, Doppler frequency domain band limiting unit 3, inverse fast Fourier transform unit 4, hit region band limiting unit 5, and imbalance correction unit 6 in the radar signal processing device 100.
- Some of the functions of the frequency estimation unit 7, the phase imbalance correction unit 8, the sub-aperture division unit 9, the restoration algorithm unit 10, the image reproduction unit 11, and the sub-aperture coupling unit 12 are realized by dedicated hardware, and some of them are realized. May be realized by software or firmware.
- each function of the power imbalance correction unit 1, the fast Fourier transform unit 2, the Doppler frequency domain band limiting unit 3, the inverse fast Fourier transform unit 4, and the hit region band limiting unit 5 is a processing circuit as dedicated hardware. Realize the function.
- the processor 21 is stored in the memory 22. The function may be realized by reading and executing the stored program.
- the processing circuit can realize each of the above functions by hardware, software, firmware or a combination thereof.
- the radar signal processing device 100 has the hit area band limiting unit 5 that limits the band in the hit area for each of the plurality of channel signals obtained by the observation by the multi-channel SAR. And an imbalance correction unit 6 that performs imbalance correction between channel signals based on a plurality of channel signals for which the hit region band limiting unit 5 has band-limited in the hit region.
- a signal in which the azimuth ambiguity is separated can be obtained due to the band limitation in the hit region.
- the accuracy of imbalance correction between channel signals can be improved.
- the hit region band limiting unit 5 in the radar signal processing device 100 extracts a channel signal in a band centered on the hit center by band limiting in the hit region. According to the above configuration, a signal in which the azimuth ambiguity is separated can be obtained by band limitation for extracting the channel signal in the band centered on the hit center. By performing imbalance correction on the signal obtained thereby, the accuracy of imbalance correction between channel signals can be improved.
- the radar signal processing device 100 has a high-speed Fourier conversion unit 2 that performs Doppler frequency conversion on a plurality of channel signals, and a plurality of channels that the high-speed Fourier conversion unit 2 performs Doppler frequency conversion.
- the Doppler frequency domain band limiting unit 3 that limits the band in the Doppler frequency domain and the plurality of channel signals that the Doppler frequency domain band limiting section 3 band limits in the Doppler frequency domain, respectively.
- Inverse high-speed Fourier conversion unit 4 that performs inverse Doppler frequency conversion, and the hit region band limiting unit 5 for each of a plurality of channel signals that the inverse high-speed Fourier conversion unit 4 has performed inverse Doppler frequency conversion.
- Band limitation in the hit area For each of the signals, the Doppler frequency domain band limiting unit 3 that limits the band in the Doppler frequency domain and the plurality of channel signals that the Doppler frequency domain band limiting section 3 band limits in the Doppler frequency domain.
- a signal in which the azimuth ambiguity is separated can be obtained due to the band limitation in the hit region and the band limitation in the Doppler frequency region.
- the Doppler frequency domain band limiting unit 3 in the radar signal processing device 100 extracts a channel signal in a band centered on the Doppler center frequency by band limiting in the Doppler frequency domain.
- a signal in which the azimuth ambiguity is separated can be obtained by band limitation for extracting the channel signal in the band centered on the Doppler center frequency.
- the imbalance correction unit 6 in the radar signal processing device 100 corrects the phase imbalance as the imbalance correction. According to the above configuration, the accuracy of the phase imbalance correction between the channel signals can be improved by performing the phase imbalance correction on the signal whose azimuth ambiguity is separated by the band limitation.
- the imbalance correction unit 6 in the radar signal processing device 100 may correct the range sample deviation as the imbalance correction. According to the above configuration, it is possible to improve the correction accuracy of the range sample deviation between the channel signals by correcting the range sample deviation for the signal whose azimuth ambiguity is separated by the band limitation.
- the observation by the multi-channel SAR in the radar signal processing apparatus 100 according to the first embodiment is a sliding spot observation or a full spotlight observation. According to the above configuration, when the azimuth ambiguity is multiplexed by sliding spot observation or full spotlight observation, the accuracy of imbalance correction between channel signals can be improved.
- a hit area band limiting step and a hit area band limiting step for band limiting in the hit region are performed for a plurality of channel signals obtained by observation by the multi-channel SAR, respectively.
- the radar processing device can be used in the multi-channel SAR technique because it can improve the accuracy of imbalance correction between channel signals.
- 1 power imbalance correction unit 1 power imbalance correction unit, 2 high-speed Fourier transform unit, 3 Doppler frequency domain band limiting unit, 4 inverse fast Fourier transform unit, 5 hit region band limiting unit, 6 imbalance correction unit, 7 Doppler center frequency estimation unit, 8 phase Imbalance correction unit, 9 sub-aperture division unit, 10 restoration algorithm unit, 11 image reproduction unit, 12 sub-aperture coupling unit, 20 processing circuit, 21 processor, 22 memory, 100 radar signal processing device.
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
A radar signal processing device (100) comprises: a hit-region band limiting unit (5) for limiting a band in a hit region with respect to each of a plurality of channel signals obtained by observation using a multichannel SAR; and an imbalance correction unit (6) for correcting imbalance between channel signals on the basis of the plurality of channel signals that have been subjected to the band limitation in the hit region by the hit-region band limiting unit (5).
Description
本開示は、レーダ信号処理装置に関する。
This disclosure relates to a radar signal processing device.
マルチチャネル合成開口レーダ(以下、マルチチャネルSARと称する)は、複数の受信アンテナを備えていることにより、複数のチャネルを介して複数の信号(以下、複数のチャネル信号と称する)を受信する。マルチチャネル合成開口レーダは、受信した複数のチャネル信号に基づいて、復元処理を行うことによりSAR画像を生成する。その際、マルチチャネル合成開口レーダは、チャネルの相違に伴って生じるチャネル信号間のずれ量(以下、インバランスと称する)を補正する(例えば、非特許文献1)。
The multi-channel synthetic aperture radar (hereinafter referred to as multi-channel SAR) is provided with a plurality of receiving antennas, so that a plurality of signals (hereinafter referred to as a plurality of channel signals) are received via a plurality of channels. The multi-channel synthetic aperture radar generates a SAR image by performing a restoration process based on a plurality of received channel signals. At that time, the multi-channel synthetic aperture radar corrects the amount of deviation between channel signals (hereinafter referred to as imbalance) caused by the difference in channels (for example, Non-Patent Document 1).
上述のマルチチャネル合成開口レーダにおいて、例えば、さらなるアジマス高分解能化を実現するためにスライディングスポットライト観測又はフルスポットライト観測等を行う場合、アジマスアンビギュイティが多重化してしまい、十分な精度で上述のチャネル信号間のインバランスを補正することが困難になるという問題がある。
In the above-mentioned multi-channel synthetic aperture radar, for example, when performing sliding spotlight observation or full spotlight observation in order to realize further high resolution of azimuth, the azimuth ambiguity is multiplexed, and the above-mentioned is performed with sufficient accuracy. There is a problem that it becomes difficult to correct the imbalance between the channel signals of.
本開示は、上記のような問題点を解決するためになされたものであり、チャネル信号間のインバランス補正の精度を向上させる技術を提供することを目的とする。
The present disclosure has been made to solve the above-mentioned problems, and an object of the present disclosure is to provide a technique for improving the accuracy of imbalance correction between channel signals.
本開示に係るレーダ信号処理装置は、マルチチャネルSARによる観測によって得られた複数のチャネル信号に対して、それぞれ、ヒット領域における帯域制限を行うヒット領域帯域制限部と、ヒット領域帯域制限部がヒット領域における帯域制限を行った複数のチャネル信号に基づいて、チャネル信号間のインバランス補正を行うインバランス補正部と、を備えている。
In the radar signal processing apparatus according to the present disclosure, the hit region band limiting unit and the hit region band limiting unit that limit the bandwidth in the hit region hit a plurality of channel signals obtained by observation by the multi-channel SAR, respectively. It includes an imbalance correction unit that performs imbalance correction between channel signals based on a plurality of channel signals for which band limitation has been performed in a region.
本開示によれば、チャネル信号間のインバランス補正の精度を向上させることができる。
According to the present disclosure, the accuracy of imbalance correction between channel signals can be improved.
以下、本開示をより詳細に説明するため、本開示を実施するための形態について、添付の図面に従って説明する。
実施の形態1.
図1は、実施の形態1に係るレーダ信号処理装置100の構成を示すブロック図である。図1が示すように、レーダ信号処理装置100は、電力インバランス補正部1、高速フーリエ変換部2、ドップラー周波数領域帯域制限部3、逆高速フーリエ変換部4、ヒット領域帯域制限部5、インバランス補正部6、サブアパーチャー分割部9、復元アルゴリズム部10、画像再生部11、及びサブアパーチャー結合部12を備えている。 Hereinafter, in order to explain the present disclosure in more detail, a mode for carrying out the present disclosure will be described with reference to the accompanying drawings.
Embodiment 1.
FIG. 1 is a block diagram showing a configuration of a radarsignal processing device 100 according to the first embodiment. As shown in FIG. 1, the radar signal processing device 100 includes a power imbalance correction unit 1, a fast Fourier transform unit 2, a Doppler frequency domain band limiting unit 3, an inverse fast Fourier transform unit 4, a hit region band limiting unit 5, and an in. It includes a balance correction unit 6, a sub-aperture division unit 9, a restoration algorithm unit 10, an image reproduction unit 11, and a sub-aperture coupling unit 12.
実施の形態1.
図1は、実施の形態1に係るレーダ信号処理装置100の構成を示すブロック図である。図1が示すように、レーダ信号処理装置100は、電力インバランス補正部1、高速フーリエ変換部2、ドップラー周波数領域帯域制限部3、逆高速フーリエ変換部4、ヒット領域帯域制限部5、インバランス補正部6、サブアパーチャー分割部9、復元アルゴリズム部10、画像再生部11、及びサブアパーチャー結合部12を備えている。 Hereinafter, in order to explain the present disclosure in more detail, a mode for carrying out the present disclosure will be described with reference to the accompanying drawings.
FIG. 1 is a block diagram showing a configuration of a radar
電力インバランス補正部1は、マルチチャネルSARによる観測によって得られた複数のチャネル信号(生データ)に対して、電力インバランス補正を行う。電力インバランス補正部1は、電力インバランス補正を行った複数のチャネル信号を高速フーリエ変換部2に出力する。
The power imbalance correction unit 1 corrects the power imbalance for a plurality of channel signals (raw data) obtained by observation by the multi-channel SAR. The power imbalance correction unit 1 outputs a plurality of channel signals that have undergone power imbalance correction to the fast Fourier transform unit 2.
なお、本明細書において、「複数のチャネル信号」は、それぞれ、マルチチャネルSARにおける複数のチャネルのうちの対応するチャネルから得られた信号を意味する。また、本明細書において、「電力インバランス補正」は、チャネルの相違に伴って生じるチャネル信号間の電力のずれ量を補正することを意味する。実施の形態1では、マルチチャネルSARによる観測は、スライディングスポット観測である。なお、実施の形態1では、マルチチャネルSARによる観測がスライディングスポット観測である構成について説明するが、マルチチャネルSARによる観測は、例えば、フルスポットライト観測等であってもよい。
In the present specification, the "plurality of channel signals" means signals obtained from the corresponding channels among the plurality of channels in the multi-channel SAR, respectively. Further, in the present specification, "power imbalance correction" means to correct the amount of power deviation between channel signals caused by a channel difference. In the first embodiment, the observation by the multi-channel SAR is a sliding spot observation. In the first embodiment, the configuration in which the observation by the multi-channel SAR is the sliding spot observation will be described, but the observation by the multi-channel SAR may be, for example, a full spotlight observation or the like.
高速フーリエ変換部2(ドップラー周波数変換部)は、マルチチャネルSARによる観測によって得られた複数のチャネル信号に対して、それぞれ、ドップラー周波数変換を行う。換言すれば、高速フーリエ変換部2(ドップラー周波数変換部)は、マルチチャネルSARによる観測によって得られた複数のチャネル信号に対して、それぞれ、高速フーリエ変換を行うことによりドップラー周波数領域の信号に変換する。より詳細には、実施の形態1では、高速フーリエ変換部2は、電力インバランス補正部1が電力インバランス補正を行った複数のチャネル信号に対して、それぞれ、ドップラー周波数変換を行う。高速フーリエ変換部2は、ドップラー周波数変換を行った複数のチャネル信号をドップラー周波数領域帯域制限部3に出力する。
The fast Fourier transform unit 2 (Doppler frequency conversion unit) performs Doppler frequency conversion for each of a plurality of channel signals obtained by observation by the multi-channel SAR. In other words, the fast Fourier transform unit 2 (Doppler frequency transform unit) converts a plurality of channel signals obtained by observation by the multi-channel SAR into signals in the Doppler frequency domain by performing a fast Fourier transform. do. More specifically, in the first embodiment, the fast Fourier transform unit 2 performs Doppler frequency conversion for each of the plurality of channel signals for which the power imbalance correction unit 1 has performed the power imbalance correction. The fast Fourier transform unit 2 outputs a plurality of channel signals that have undergone Doppler frequency conversion to the Doppler frequency domain band limiting unit 3.
ドップラー周波数領域帯域制限部3(バンドパスフィルタ)は、高速フーリエ変換部2がドップラー周波数変換を行った複数のチャネル信号に対して、それぞれ、ドップラー周波数領域における帯域制限を行う。より詳細には、実施の形態1では、ドップラー周波数領域帯域制限部3は、ドップラー周波数領域における帯域制限によって、ドップラー中心周波数を中心とした帯域のチャネル信号を抽出する。ドップラー周波数領域帯域制限部3は、ドップラー周波数領域における帯域制限を行った複数のチャネル信号を逆高速フーリエ変換部4に出力する。
The Doppler frequency domain band limiting unit 3 (bandpass filter) limits the band in the Doppler frequency domain for each of the plurality of channel signals for which the fast Fourier transform unit 2 has undergone the Doppler frequency conversion. More specifically, in the first embodiment, the Doppler frequency domain band limiting unit 3 extracts a channel signal in a band centered on the Doppler center frequency by band limiting in the Doppler frequency domain. The Doppler frequency domain band limiting unit 3 outputs a plurality of channel signals band-limited in the Doppler frequency domain to the inverse fast Fourier transform unit 4.
逆高速フーリエ変換部4(逆ドップラー周波数変換部)は、ドップラー周波数領域帯域制限部3がドップラー周波数領域における帯域制限を行った複数のチャネル信号に対して、それぞれ、逆ドップラー周波数変換を行う。換言すれば、逆高速フーリエ変換部4は、ドップラー周波数領域帯域制限部3がドップラー周波数領域における帯域制限を行った複数のチャネル信号に対して、それぞれ、逆高速フーリエ変換を行うことにより、高速フーリエ変換部2が高速フーリエ変換を行う前の元の時間領域の信号に変換する。逆高速フーリエ変換部4は、逆ドップラー周波数変換を行った複数のチャネル信号をヒット領域帯域制限部5に出力する。
The inverse fast Fourier transform unit 4 (inverse Doppler frequency converter) performs inverse Doppler frequency conversion for each of a plurality of channel signals for which the Doppler frequency domain band limiting unit 3 has band-limited in the Doppler frequency domain. In other words, the inverse fast Fourier transform unit 4 performs an inverse fast Fourier transform on each of a plurality of channel signals for which the Doppler frequency domain band limiting unit 3 has band-limited in the Doppler frequency domain, thereby performing a fast Fourier transform. The transforming unit 2 converts the signal into the signal in the original time domain before performing the fast Fourier transform. The inverse fast Fourier transform unit 4 outputs a plurality of channel signals subjected to the inverse Doppler frequency conversion to the hit region band limiting unit 5.
ヒット領域帯域制限部5(バンドパスフィルタ)は、マルチチャネルSARによる観測によって得られた複数のチャネル信号に対して、それぞれ、ヒット領域における帯域制限を行う。より詳細には、実施の形態1では、ヒット領域帯域制限部5は、逆高速フーリエ変換部4が逆ドップラー周波数変換を行った複数のチャネル信号に対して、それぞれ、ヒット領域における帯域制限を行う。さらに詳細には、実施の形態1では、ヒット領域帯域制限部5は、ヒット領域における帯域制限によって、ヒット中心を中心とした帯域のチャネル信号を抽出する。ヒット領域帯域制限部5は、ヒット領域における帯域制限を行った複数のチャネル信号をインバランス補正部6に出力する。
The hit area band limiting unit 5 (bandpass filter) limits the band in the hit area for each of the plurality of channel signals obtained by the observation by the multi-channel SAR. More specifically, in the first embodiment, the hit region band limiting unit 5 limits the bandwidth in the hit region for each of the plurality of channel signals subjected to the inverse Doppler frequency conversion by the inverse fast Fourier transform unit 4. .. More specifically, in the first embodiment, the hit region band limiting unit 5 extracts a channel signal in a band centered on the hit center by band limiting in the hit region. The hit area band limiting unit 5 outputs a plurality of channel signals for which the band is limited in the hit area to the imbalance correction unit 6.
インバランス補正部6は、ヒット領域帯域制限部5がヒット領域における帯域制限を行った複数のチャネル信号に基づいて、チャネル信号間のインバランス補正を行う。より詳細には、実施の形態1では、インバランス補正部6は、チャネル信号間のインバランス補正として、位相インバランスの補正を行う。インバランス補正部6は、インバランス補正を行った複数のチャネル信号をサブアパーチャー分割部9に出力する。
The imbalance correction unit 6 performs imbalance correction between channel signals based on a plurality of channel signals for which the hit region band limiting unit 5 has band-limited in the hit region. More specifically, in the first embodiment, the imbalance correction unit 6 corrects the phase imbalance as the imbalance correction between the channel signals. The imbalance correction unit 6 outputs a plurality of channel signals for which imbalance correction has been performed to the sub-aperture division unit 9.
なお、本明細書において、「位相インバランス」は、チャネルの相違に伴って生じるチャネル信号間の位相のずれ量を意味する。実施の形態1では、インバランス補正部6がインバランス補正として位相インバランスの補正を行う構成について説明するが、インバランス補正部6は、チャネル信号間のインバランス補正として、レンジサンプルずれの補正を行ってもよい。
In the present specification, "phase imbalance" means the amount of phase shift between channel signals caused by the difference in channels. In the first embodiment, the configuration in which the imbalance correction unit 6 corrects the phase imbalance as the imbalance correction will be described, but the imbalance correction unit 6 corrects the range sample deviation as the imbalance correction between the channel signals. May be done.
さらに詳細には、実施の形態1では、インバランス補正部6は、ドップラー中心周波数推定部7、及び位相インバランス補正部8を備えている。
ドップラー中心周波数推定部7は、ヒット領域帯域制限部5がヒット領域における帯域制限を行った複数のチャネル信号のドップラー中心周波数を推定する。ドップラー中心周波数推定部7は、推定したドップラー中心周波数を位相インバランス補正部8に出力する。 More specifically, in the first embodiment, theimbalance correction unit 6 includes a Doppler center frequency estimation unit 7 and a phase imbalance correction unit 8.
The Doppler center frequency estimation unit 7 estimates the Doppler center frequencies of a plurality of channel signals for which the hit regionband limiting unit 5 has band-limited in the hit region. The Doppler center frequency estimation unit 7 outputs the estimated Doppler center frequency to the phase imbalance correction unit 8.
ドップラー中心周波数推定部7は、ヒット領域帯域制限部5がヒット領域における帯域制限を行った複数のチャネル信号のドップラー中心周波数を推定する。ドップラー中心周波数推定部7は、推定したドップラー中心周波数を位相インバランス補正部8に出力する。 More specifically, in the first embodiment, the
The Doppler center frequency estimation unit 7 estimates the Doppler center frequencies of a plurality of channel signals for which the hit region
位相インバランス補正部8は、ドップラー中心周波数推定部7が推定したドップラー中心周波数に基づいて、ヒット領域帯域制限部5がヒット領域における帯域制限を行った複数のチャネル信号に基づいて、チャネル信号間の位相インバランスの補正を行う。位相インバランス補正部8は、位相インバランスの補正を行った複数のチャネル信号をサブアパーチャー分割部9に出力する。
The phase imbalance correction unit 8 is based on the Doppler center frequency estimated by the Doppler center frequency estimation unit 7, and is based on a plurality of channel signals for which the hit region band limiting unit 5 has band-limited in the hit region. The phase imbalance of is corrected. The phase imbalance correction unit 8 outputs a plurality of channel signals for which the phase imbalance has been corrected to the sub-aperture division unit 9.
サブアパーチャー分割部9は、インバランス補正部6がインバランス補正を行った複数のチャネル信号に対して、サブアパーチャー分割を行う。より詳細には、実施の形態1では、サブアパーチャー分割部9は、インバランス補正部6の位相インバランス補正部8が位相インバランスの補正を行った複数のチャネル信号に対して、サブアパーチャー分割を行う。サブアパーチャー分割部9は、サブアパーチャー分割を行った複数のチャネル信号を復元アルゴリズム部10に出力する。
The sub-aperture division unit 9 performs sub-aperture division for a plurality of channel signals for which the imbalance correction unit 6 has performed imbalance correction. More specifically, in the first embodiment, the sub-aperture division unit 9 divides the sub-aperture into a plurality of channel signals for which the phase imbalance correction unit 8 of the imbalance correction unit 6 has corrected the phase imbalance. I do. The sub-aperture division unit 9 outputs a plurality of channel signals obtained by sub-aperture division to the restoration algorithm unit 10.
復元アルゴリズム部10は、サブアパーチャー分割部9がサブアパーチャー分割を行った複数のチャネル信号に対して、復元アルゴリズムを行う。復元アルゴリズム部10は、復元アルゴリズムを行った複数のチャネル信号を画像再生部11に出力する。
The restoration algorithm unit 10 performs a restoration algorithm on a plurality of channel signals for which the sub-aperture division unit 9 has performed sub-aperture division. The restoration algorithm unit 10 outputs a plurality of channel signals subjected to the restoration algorithm to the image reproduction unit 11.
画像再生部11は、復元アルゴリズム部10が復元アルゴリズムを行った複数のチャネル信号に対して、画像再生を行う。画像再生部11は、画像再生を行った複数のチャネル信号をサブアパーチャー結合部12に出力する。
The image reproduction unit 11 reproduces an image on a plurality of channel signals for which the restoration algorithm unit 10 has performed the restoration algorithm. The image reproduction unit 11 outputs a plurality of channel signals for which image reproduction has been performed to the sub-aperture coupling unit 12.
サブアパーチャー結合部12は、画像再生部11が画像再生を行った複数のチャネル信号に対してサブアパーチャー結合を行うことによりSAR画像を生成する。例えば、サブアパーチャー結合部12が生成したSAR画像は、図示しないディスプレイによって表示される。
The sub-aperture coupling unit 12 generates a SAR image by performing sub-aperture coupling on a plurality of channel signals for which the image reproduction unit 11 has reproduced an image. For example, the SAR image generated by the sub-aperture coupling portion 12 is displayed by a display (not shown).
以下で、実施の形態1に係るレーダ信号処理装置100の動作について図面を参照して説明する。図2は、実施の形態1に係るレーダ信号処理装置100によるレーダ信号処理方法を示すフローチャートである。なお、以下で説明する各ステップの前に、レーダ信号処理装置100は、マルチチャネルSARによる観測によって得られた複数のチャネル信号を取得したものとする。
Hereinafter, the operation of the radar signal processing device 100 according to the first embodiment will be described with reference to the drawings. FIG. 2 is a flowchart showing a radar signal processing method by the radar signal processing device 100 according to the first embodiment. Before each step described below, it is assumed that the radar signal processing device 100 has acquired a plurality of channel signals obtained by observation by the multi-channel SAR.
図2が示すように、電力インバランス補正部1は、マルチチャネルSARによる観測によって得られた複数のチャネル信号に対して、電力インバランス補正を行う(ステップST1)。電力インバランス補正部1は、電力インバランス補正を行った複数のチャネル信号を高速フーリエ変換部2に出力する。
As shown in FIG. 2, the power imbalance correction unit 1 performs power imbalance correction on a plurality of channel signals obtained by observation by the multi-channel SAR (step ST1). The power imbalance correction unit 1 outputs a plurality of channel signals that have undergone power imbalance correction to the fast Fourier transform unit 2.
次に、高速フーリエ変換部2は、電力インバランス補正部1が電力インバランス補正を行った複数のチャネル信号に対して、それぞれ、ドップラー周波数変換を行う(ステップST2)。高速フーリエ変換部2は、ドップラー周波数変換を行った複数のチャネル信号をドップラー周波数領域帯域制限部3に出力する。
Next, the fast Fourier transform unit 2 performs Doppler frequency conversion for each of the plurality of channel signals for which the power imbalance correction unit 1 has performed the power imbalance correction (step ST2). The fast Fourier transform unit 2 outputs a plurality of channel signals that have undergone Doppler frequency conversion to the Doppler frequency domain band limiting unit 3.
次に、ドップラー周波数領域帯域制限部3は、高速フーリエ変換部2がドップラー周波数変換を行った複数のチャネル信号に対して、それぞれ、ドップラー周波数領域における帯域制限を行う(ステップST3)。ドップラー周波数領域帯域制限部3は、ドップラー周波数領域における帯域制限を行った複数のチャネル信号を逆高速フーリエ変換部4に出力する。
Next, the Doppler frequency domain band limiting unit 3 limits the band in the Doppler frequency domain for each of the plurality of channel signals subjected to the Doppler frequency conversion by the fast Fourier transform unit 2 (step ST3). The Doppler frequency domain band limiting unit 3 outputs a plurality of channel signals band-limited in the Doppler frequency domain to the inverse fast Fourier transform unit 4.
次に、逆高速フーリエ変換部4は、ドップラー周波数領域帯域制限部3がドップラー周波数領域における帯域制限を行った複数のチャネル信号に対して、それぞれ、逆ドップラー周波数変換を行う(ステップST4)。逆高速フーリエ変換部4は、逆ドップラー周波数変換を行った複数のチャネル信号をヒット領域帯域制限部5に出力する。
Next, the inverse fast Fourier transform unit 4 performs inverse Doppler frequency conversion for each of the plurality of channel signals for which the Doppler frequency domain band limiting unit 3 has band-limited in the Doppler frequency domain (step ST4). The inverse fast Fourier transform unit 4 outputs a plurality of channel signals subjected to the inverse Doppler frequency conversion to the hit region band limiting unit 5.
次に、ヒット領域帯域制限部5は、逆高速フーリエ変換部4が逆ドップラー周波数変換を行った複数のチャネル信号に対して、それぞれ、ヒット領域における帯域制限を行う(ステップST5)。ヒット領域帯域制限部5は、ヒット領域における帯域制限を行った複数のチャネル信号をインバランス補正部6に出力する。
Next, the hit region band limiting unit 5 limits the bandwidth in the hit region for each of the plurality of channel signals for which the inverse fast Fourier transform unit 4 has undergone the inverse Doppler frequency conversion (step ST5). The hit area band limiting unit 5 outputs a plurality of channel signals for which the band is limited in the hit area to the imbalance correction unit 6.
次に、インバランス補正部6のドップラー中心周波数推定部7は、ヒット領域帯域制限部5がヒット領域における帯域制限を行った複数のチャネル信号のドップラー中心周波数を推定する(ステップST6)。ドップラー中心周波数推定部7は、推定したドップラー中心周波数を位相インバランス補正部8に出力する。
Next, the Doppler center frequency estimation unit 7 of the imbalance correction unit 6 estimates the Doppler center frequency of a plurality of channel signals for which the hit region band limiting unit 5 has band-limited in the hit region (step ST6). The Doppler center frequency estimation unit 7 outputs the estimated Doppler center frequency to the phase imbalance correction unit 8.
次に、インバランス補正部6の位相インバランス補正部8は、ドップラー中心周波数推定部7が推定したドップラー中心周波数に基づいて、チャネル信号間の位相インバランスの補正を行う(ステップST7)。位相インバランス補正部8は、位相インバランスの補正を行った複数のチャネル信号をサブアパーチャー分割部9に出力する。
Next, the phase imbalance correction unit 8 of the imbalance correction unit 6 corrects the phase imbalance between the channel signals based on the Doppler center frequency estimated by the Doppler center frequency estimation unit 7 (step ST7). The phase imbalance correction unit 8 outputs a plurality of channel signals for which the phase imbalance has been corrected to the sub-aperture division unit 9.
次に、サブアパーチャー分割部9は、インバランス補正部6が位相インバランスの補正を行った複数のチャネル信号に対して、サブアパーチャー分割を行う(ステップST8)。サブアパーチャー分割部9は、サブアパーチャー分割を行った複数のチャネル信号を復元アルゴリズム部10に出力する。
Next, the sub-aperture division unit 9 performs sub-aperture division for a plurality of channel signals for which the imbalance correction unit 6 has corrected the phase imbalance (step ST8). The sub-aperture division unit 9 outputs a plurality of channel signals obtained by sub-aperture division to the restoration algorithm unit 10.
次に、復元アルゴリズム部10は、サブアパーチャー分割部9がサブアパーチャー分割を行った複数のチャネル信号に対して、復元アルゴリズムを行う(ステップST9)。復元アルゴリズム部10は、復元アルゴリズムを行った複数のチャネル信号を画像再生部11に出力する。
Next, the restoration algorithm unit 10 performs a restoration algorithm on a plurality of channel signals for which the sub-aperture division unit 9 has undergone sub-aperture division (step ST9). The restoration algorithm unit 10 outputs a plurality of channel signals subjected to the restoration algorithm to the image reproduction unit 11.
次に、画像再生部11は、復元アルゴリズム部10が復元アルゴリズムを行った複数のチャネル信号に対して、画像再生を行う(ステップST10)。画像再生部11は、画像再生を行った複数のチャネル信号をサブアパーチャー結合部12に出力する。
Next, the image reproduction unit 11 reproduces an image on a plurality of channel signals for which the restoration algorithm unit 10 has performed the restoration algorithm (step ST10). The image reproduction unit 11 outputs a plurality of channel signals for which image reproduction has been performed to the sub-aperture coupling unit 12.
サブアパーチャー結合部12は、画像再生部11が画像再生を行った複数のチャネル信号に対してサブアパーチャー結合を行うことによりSAR画像を生成する(ステップST11)。例えば、サブアパーチャー結合部12がステップST11で生成したSAR画像は、図示しないディスプレイによって表示される。
The sub-aperture coupling unit 12 generates a SAR image by performing sub-aperture coupling on a plurality of channel signals for which the image reproduction unit 11 has reproduced an image (step ST11). For example, the SAR image generated by the sub-aperture coupling unit 12 in step ST11 is displayed by a display (not shown).
以下で、実施の形態1に係るレーダ信号処理装置100によるレーダ信号処理方法の具体例について図面を参照して説明する。図3は、実施の形態1に係るレーダ信号処理装置100によるレーダ信号処理方法によって処理されたチャネル信号の途中経過を示す図である。
Hereinafter, a specific example of the radar signal processing method by the radar signal processing device 100 according to the first embodiment will be described with reference to the drawings. FIG. 3 is a diagram showing the progress of a channel signal processed by the radar signal processing method by the radar signal processing device 100 according to the first embodiment.
当該具体例では、マルチチャネルSARによる観測によって得られた複数のチャネル信号は、2チャネルのスライディングスポットライト観測で得られた信号である。レーダ信号処理装置100によって処理される前の第1のチャネル信号(チャネル0の信号)をs0(τ,η)とし、第2のチャネル信号(チャネル1の信号)をs1(τ,η)とする。τは、レンジ時間であり、ηは、ヒットである。なお、以下では、第1のチャネル信号を基準チャネルとする。
In the specific example, the plurality of channel signals obtained by the observation by the multi-channel SAR are the signals obtained by the two-channel sliding spotlight observation. The first channel signal (channel 0 signal) before being processed by the radar signal processing device 100 is s 0 (τ, η), and the second channel signal (channel 1 signal) is s 1 (τ, η). ). τ is the range time and η is a hit. In the following, the first channel signal will be used as a reference channel.
当該具体例では、上述のステップST1において、電力インバランス補正部1は、第1のチャネル信号s0(τ,η)、及び第2のチャネル信号s1(τ,η)に対して、以下で説明する方法に従って電力インバランス補正を行う(例えば、非特許文献J. H. Kim, M. Younis, P. P. Iraola, M. Gebele, and G. Krieger, “First spaceborne demonstration of digital beamforming for azimuth ambiguity suppression”, IEEE Trans. Geosci. Remote Sens., vol.51, no.1, pp.579-590, Jan., 2013を参照)。
In the specific example, in the above-mentioned step ST1, the power imbalance correction unit 1 has the following with respect to the first channel signal s 0 (τ, η) and the second channel signal s 1 (τ, η). Power imbalance correction is performed according to the method described in (For example, Non-Patent Document J. H. Kim, M. Younis, P. P. Iraola, M. Gebele, and G. Krieger, “First spaceborne demonstration of digital beamforming for azimuth ambiguity suppression”, IEEE. See Trans. Geosci. Remote Sens., Vol.51, no.1, pp.579-590, Jan., 2013).
まず、電力インバランス補正部1は、レンジ周波数領域で電力インバランスの補正を行うために、第1のチャネル信号s0(τ,η)及び第2のチャネル信号s1(τ,η)に対してレンジ周波数変換を行う。これにより、電力インバランス補正部1は、レンジ周波数をfrとしてS0(fr,η)、及びS1(fr,η)を得る。
First, the power imbalance correction unit 1 uses the first channel signal s 0 (τ, η) and the second channel signal s 1 (τ, η) to correct the power imbalance in the range frequency region. On the other hand, range frequency conversion is performed. As a result, the power imbalance correction unit 1 obtains S 0 ( fr , η) and S 1 ( fr , η) with the range frequency as fr .
次に、電力インバランス補正部1は、第1のチャネル信号で発生する雑音電力をσ0
2とし、第2のチャネル信号で発生する雑音電力をσ1
2として、以下の式(1)に従って、電力(振幅)インバランスの補正係数α(fr)を算出する。
Next, the powerimbalance correction unit 1 sets the noise power generated by the first channel signal as σ 0 2 , and the noise power generated by the second channel signal as σ 1 2 , according to the following equation (1). , The correction coefficient α ( fr ) of the power (amplitude) imbalance is calculated.
Next, the power
式(1)において、Eη(|S0(fr,η)|2)は、|S0(fr,η)|2のヒット方向の平均を示し、Eη(|S1(fr,η)|2)は、|S1(fr,η)|2のヒット方向の平均を示す。ただし、雑音電力が未知である場合には、σ0
2=σ1
2=0としてもよい。
In the equation (1), E η (| S 0 ( fr , η) | 2 ) indicates the average of | S 0 ( fr , η) | 2 in the hit direction, and E η (| S 1 (f)). r , η) | 2 ) indicates the average of | S 1 ( fr , η) | 2 in the hit direction. However, when the noise power is unknown, σ 0 2 = σ 1 2 = 0 may be set.
次に、電力インバランス補正部1は、算出した補正係数α(fr)に基づいて、以下の式(2)に従って、S1(fr,η)の振幅を補正する。
Next, the powerimbalance correction unit 1 corrects the amplitude of S 1 ( fr , η ) according to the following equation (2) based on the calculated correction coefficient α (fr).
Next, the power
次に、電力インバランス補正部1は、式(2)により得られたS1’(fr,η)に対してレンジ時間変換を行うことにより、電力インバランスの補正が行われた第2のチャネル信号s1’(τ,η)を得る。
なお、電力インバランス補正部1は、電力インバランスを補正する方法として、当該具体例の方法以外の方法を用いてもよい。 Next, the powerimbalance correction unit 1 corrects the power imbalance by performing range time conversion on S 1 '( fr , η) obtained by the equation (2). Channel signal s 1 '(τ, η) of.
The powerimbalance correction unit 1 may use a method other than the method of the specific example as a method of correcting the power imbalance.
なお、電力インバランス補正部1は、電力インバランスを補正する方法として、当該具体例の方法以外の方法を用いてもよい。 Next, the power
The power
図3のAは、ステップST1の処理後の第1のチャネル信号s0(τ,η)のグラフを示す。図3のAにおける縦軸は、ヒットを示し、横軸は、レンジを示す。図3のAが示すように、7点の点像の生データがヒット-レンジ空間に多重化されている。
FIG. 3A shows a graph of the first channel signal s 0 (τ, η) after the processing of step ST1. In FIG. 3, the vertical axis in A indicates a hit, and the horizontal axis indicates a range. As shown by A in FIG. 3, the raw data of the seven point images are multiplexed in the hit-range space.
次に、上述のステップST2において、高速フーリエ変換部2は、第1のチャネル信号s0(τ,η)及び第2のチャネル信号s1’(τ,η)に対して、それぞれ、ドップラー周波数変換を行う。図3のBは、ステップST2の処理後の第1のチャネル信号のグラフを示す。図3のBにおける縦軸は、ドップラー周波数を示し、横軸は、レンジを示す。図3のBが示すように、第1のチャネル信号s0(τ,η)は、アジマスアンビギュイティが多重化している。なお、図示しないが、第2のチャネル信号s1’(τ,η)についても同様に、アジマスアンビギュイティが多重化している。
Next, in step ST2 described above, the fast Fourier transform unit 2 has a Doppler frequency with respect to the first channel signal s 0 (τ, η) and the second channel signal s 1 '(τ, η), respectively. Perform the conversion. FIG. 3B shows a graph of the first channel signal after the processing of step ST2. In FIG. 3, the vertical axis in B indicates the Doppler frequency, and the horizontal axis indicates the range. As shown by B in FIG. 3, the first channel signal s 0 (τ, η) is multiplexed with azimuth ambiguity. Although not shown, the azimuth ambiguity is similarly multiplexed for the second channel signal s 1 '(τ, η).
次に、上述のステップST3において、ドップラー周波数領域帯域制限部3は、ステップST2の処理によって得られた第1のチャネル信号及び第2のチャネル信号の各ドップラー周波数成分に対して、ドップラー周波数領域における帯域制限を行う。図3のCは、ステップST3の処理後の第1のチャネル信号のグラフを示す。図3のCにおける縦軸は、ドップラー周波数を示し、横軸は、レンジを示す。図3のCが示すように、ステップST3の処理によりドップラー中心周波数付近の信号が抽出される。なお、図示しないが、第2のチャネル信号s1’(τ,η)についても同様に、ステップST3の処理によりドップラー中心周波数付近の信号が抽出される。
Next, in step ST3 described above, the Doppler frequency domain band limiting unit 3 in the Doppler frequency domain for each Doppler frequency component of the first channel signal and the second channel signal obtained by the processing of step ST2. Bandwidth limitation. FIG. 3C shows a graph of the first channel signal after the processing of step ST3. The vertical axis in C in FIG. 3 indicates the Doppler frequency, and the horizontal axis indicates the range. As shown by C in FIG. 3, a signal near the Doppler center frequency is extracted by the process of step ST3. Although not shown, for the second channel signal s 1 '(τ, η), a signal near the Doppler center frequency is similarly extracted by the process of step ST3.
次に、上述のステップST4において、逆高速フーリエ変換部4は、ステップST3の処理によって得られた第1のチャネル信号及び第2のチャネル信号に対して逆ドップラー周波数変換を行う。図3のDは、ステップST4の処理後の第1のチャネル信号のグラフを示す。図3のDにおける縦軸は、ヒットを示し、横軸は、レンジを示す。図3のDが示すように、ステップST4の処理後の第1のチャネル信号s0(τ,η)は、信号が離散的に分散されている。なお、図示しないが、ステップST4の処理後の第2のチャネル信号s1’(τ,η)についても同様に、信号が離散的に分散されている。
Next, in step ST4 described above, the inverse fast Fourier transform unit 4 performs inverse Doppler frequency conversion on the first channel signal and the second channel signal obtained by the processing of step ST3. FIG. 3D shows a graph of the first channel signal after the processing of step ST4. The vertical axis in D in FIG. 3 indicates a hit, and the horizontal axis indicates a range. As shown by D in FIG. 3, in the first channel signal s 0 (τ, η) after the processing of step ST4, the signals are discretely dispersed. Although not shown, the signals of the second channel signal s 1 '(τ, η) after the processing of step ST4 are also discretely dispersed.
次に、上述のステップST5において、ヒット領域帯域制限部5は、ステップST4の処理によって得られた第1のチャネル信号及び第2のチャネル信号の各ヒット-レンジ成分に対して、ヒット領域での帯域制限を行う。図3のEは、ステップST5の処理後の第1のチャネル信号のグラフを示す。図3のEにおける縦軸は、ヒットを示し、横軸は、レンジを示す。図3のEが示すように、アジマスアンビギュイティを抑圧することによりアジマスアンビギュイティが分離された信号が得られる。なお、図示しないが、第2のチャネル信号s1’(τ,η)についても同様に、アジマスアンビギュイティが分離された信号が得られる。以下では、ステップST5の処理によって得られた第1のチャネル信号を、s
0(τ,η)とし、第2のチャネル信号をs
1’(τ,η)とする。
Next, in step ST5 described above, the hit region band limiting unit 5 has a hit region for each hit-range component of the first channel signal and the second channel signal obtained by the processing of step ST4. Bandwidth limitation. FIG. 3E shows a graph of the first channel signal after the processing of step ST5. The vertical axis in E in FIG. 3 indicates a hit, and the horizontal axis indicates a range. As shown by E in FIG. 3, by suppressing the azimuth ambiguity, a signal from which the azimuth ambiguity is separated is obtained. Although not shown, a signal from which the azimuth ambiguity is separated can be obtained similarly for the second channel signal s 1 '(τ, η). In the following, the first channel signal obtained by the process of step ST5 is s 0 (τ, η), and the second channel signal is s 1 '(τ, η).
次に、上述のステップST6及びステップST7において、インバランス補正部6は、ステップST5の処理によって得られた第1のチャネル信号s
0(τ,η)及び第2のチャネル信号s
1’(τ,η)に基づいて、以下で説明する方法に従って、チャネル信号間の位相インバランスの補正を行う(例えば、H. Fan, Z. Zhang, and R. Wang, “phase mismatch calibration for multichannel sliding spotlight SAR imaging with extended azimuth cross correlation”, IGARSS 2019を参照)。
Next, in step ST6 and step ST7 described above, the imbalance correction unit 6 performs the first channel signal s 0 (τ, η) and the second channel signal s 1 '(τ) obtained by the processing of step ST5. , Η), and corrects the phase imbalance between the channel signals according to the method described below (eg, H. Fan, Z. Zhang, and R. Wang, “phase mismatch calibration for multichannel sliding spotlight SAR). imaging with extended azimuth cross correlation ”, see IGARSS 2019).
まず、上述のステップST6において、インバランス補正部6のドップラー中心周波数推定部7は、第1のチャネル信号s
0(τ,η)及び第2のチャネル信号s
1’(τ,η)に基づいて、以下の式(3)及び式(4)に従って、2つの干渉データC1,0(τ)及びC0,1(τ)を算出する。
First, in step ST6 described above, the Doppler center frequency estimation unit 7 of theimbalance correction unit 6 is based on the first channel signal s 0 (τ, η) and the second channel signal s 1 '(τ, η). Then, two interference data C 1,0 (τ) and C 0,1 (τ) are calculated according to the following equations (3) and (4).
First, in step ST6 described above, the Doppler center frequency estimation unit 7 of the
式(4)において、PRFは、パルス繰り返し周波数であり、チャネル信号内のアジマス方向の1サンプルずれに相当する。式(3)及び式(4)において、Eη( )は、上述の通り、括弧内の数式のヒット方向の平均を示す。
In equation (4), the PRF is the pulse repetition frequency and corresponds to one sample shift in the azimuth direction in the channel signal. In equations (3) and (4), E η () indicates the average of the equations in parentheses in the hit direction, as described above.
次に、ドップラー中心周波数推定部7は、干渉データC1,0(τ)及びC0,1(τ)の両者の積C1,0(τ)C0,1(τ)を算出することにより、s
1’(τ,η)に含まれる位相インバランスを相殺する。ドップラー中心周波数をfdc(τ)とすると,以下の式(5)の関係が得られる。
Next, the Doppler center frequency estimation unit 7 calculates the product C 1,0 (τ) C 0,1 (τ) of both the interference data C 1,0 (τ) and C 0,1 (τ). Therefore, the phase imbalance contained in s 1 '(τ, η) is offset. Assuming that the Doppler center frequency is f dc (τ), the relationship of the following equation (5) is obtained.
Next, the Doppler center frequency estimation unit 7 calculates the product C 1,0 (τ) C 0,1 (τ) of both the interference data C 1,0 (τ) and C 0,1 (τ). Therefore, the phase imbalance contained in s 1 '(τ, η) is offset. Assuming that the Doppler center frequency is f dc (τ), the relationship of the following equation (5) is obtained.
ドップラー中心周波数推定部7は、算出したC1,0(τ)C0,1(τ)に基づいて、以下の式(6)に従って、ドップラー中心周波数fdc(τ)を算出する。
The Doppler center frequency estimation unit 7 calculates the Doppler center frequency f dc (τ) based on the calculated C 1,0 (τ) C 0,1 (τ) according to the following equation (6).
The Doppler center frequency estimation unit 7 calculates the Doppler center frequency f dc (τ) based on the calculated C 1,0 (τ) C 0,1 (τ) according to the following equation (6).
次に、上述のステップST7において、インバランス補正部6の位相インバランス補正部8は、ドップラー中心周波数推定部7が算出したドップラー中心周波数fdc(τ)に基づいて、以下の式(7)に従って、位相インバランスφ(τ)を算出する。
式(7)において、vsは、プラットフォーム速度を示し、dは、2チャネルの受信位相中心間距離を示す。 Next, in step ST7 described above, the phase imbalance correction unit 8 of theimbalance correction unit 6 has the following equation (7) based on the Doppler center frequency f dc (τ) calculated by the Doppler center frequency estimation unit 7. The phase imbalance φ (τ) is calculated according to the above.
In equation (7), vs represents the platform speed and d represents the distance between the reception phase centers of the two channels.
式(7)において、vsは、プラットフォーム速度を示し、dは、2チャネルの受信位相中心間距離を示す。 Next, in step ST7 described above, the phase imbalance correction unit 8 of the
In equation (7), vs represents the platform speed and d represents the distance between the reception phase centers of the two channels.
次に、位相インバランス補正部8は、得られた位相インバランスφ(τ)に対して,重みを生成し,ヒット方向の期待値をとるように加重平均処理を行う。この処理により,推定値の誤差を低減する。より詳細には、まず、位相インバランス補正部8は、以下の式(8)に従って、|s
0(τ,η)|のヒット方向の平均を算出する。なお、| |は、括弧内の数式の絶対値を示す。
Next, the phase imbalance correction unit 8 generates a weight for the obtained phase imbalance φ (τ) and performs a weighted averaging process so as to take an expected value in the hit direction. This process reduces the error in the estimated value. More specifically, first, the phase imbalance correction unit 8 calculates the average of | s 0 (τ, η) | in the hit direction according to the following equation (8). Note that || indicates the absolute value of the formula in parentheses.
Next, the phase imbalance correction unit 8 generates a weight for the obtained phase imbalance φ (τ) and performs a weighted averaging process so as to take an expected value in the hit direction. This process reduces the error in the estimated value. More specifically, first, the phase imbalance correction unit 8 calculates the average of | s 0 (τ, η) | in the hit direction according to the following equation (8). Note that || indicates the absolute value of the formula in parentheses.
次に、位相インバランス補正部8は、以下の式(9)に従って、加重平均をとり、位相インバランスφを算出する。
式(9)において、Eτ{ }は、括弧内の数式のレンジ時間方向の平均を示す。 Next, the phase imbalance correction unit 8 takes a weighted average according to the following equation (9) and calculates the phase imbalance φ .
In equation (9), E τ {} indicates the average of the equations in parentheses in the range time direction.
式(9)において、Eτ{ }は、括弧内の数式のレンジ時間方向の平均を示す。 Next, the phase imbalance correction unit 8 takes a weighted average according to the following equation (9) and calculates the phase imbalance φ .
In equation (9), E τ {} indicates the average of the equations in parentheses in the range time direction.
次に、位相インバランス補正部8は、算出した位相インバランスφに基づいて、以下の式(10)に従って、チャネル信号間の位相インバランスの補正を行う。より詳細には、当該具体例では、位相インバランス補正部8は、算出した位相インバランスφに基づいて、以下の式(10)に従って、上述の電力インバランス補正後の第2のチャネル信号s1’(τ,η)に対して位相インバランスの補正を行う。
なお、インバランス補正部6は、位相インバランスを補正する方法として、当該具体例の方法以外の方法を用いてもよい。 Next, the phase imbalance correction unit 8 corrects the phase imbalance between the channel signals according to the following equation (10) based on the calculated phase imbalance φ . More specifically, in the specific example, the phase imbalance correction unit 8 has the second channel signal s after the power imbalance correction described above according to the following equation (10) based on the calculated phase imbalance φ . Correct the phase imbalance for 1 '(τ, η).
Theimbalance correction unit 6 may use a method other than the method of the specific example as a method of correcting the phase imbalance.
なお、インバランス補正部6は、位相インバランスを補正する方法として、当該具体例の方法以外の方法を用いてもよい。 Next, the phase imbalance correction unit 8 corrects the phase imbalance between the channel signals according to the following equation (10) based on the calculated phase imbalance φ . More specifically, in the specific example, the phase imbalance correction unit 8 has the second channel signal s after the power imbalance correction described above according to the following equation (10) based on the calculated phase imbalance φ . Correct the phase imbalance for 1 '(τ, η).
The
次に、上述のステップST8において、サブアパーチャー分割部9は、ステップST7の処理後の第2のチャネル信号s1’’(τ,η)、及び第1のチャネル信号s0(τ,η)に対して、サブアパーチャー分割を行う(サブアパーチャー分割の方法については、例えば、非特許文献J. Mittermayer, R. Lord, and E. B¨oner,“Sliding spotlight SAR processing for TerraSAR-X using a new fomulation of the extended chirp scaling algorithm,” IGARSS2003, vol.3, p.1462-1464, 2003.を参照)。
Next, in step ST8 described above, the sub-aperture dividing unit 9 has the second channel signal s 1 '' (τ, η) and the first channel signal s 0 (τ, η) after the processing of step ST7. (For the method of sub-aperture division, for example, the non-patent document J. Mittermayer, R. Lord, and E. B ¨oner, “Sliding spotlight SAR processing for TerraSAR-X using a new” fomulation of the extended chirp scaling algorithm, ”IGARSS2003, vol.3, p.1462-1464, 2003.).
次に、上述のステップST9において、復元アルゴリズム部10は、サブアパーチャー分割部9がサブアパーチャー分割を行った第1のチャネル信号及び第2のチャネル信号に対して、復元アルゴリズムを行う(復元アルゴリズムの方法については、例えば、上述の非特許文献1を参照)。
Next, in step ST9 described above, the restoration algorithm unit 10 performs a restoration algorithm on the first channel signal and the second channel signal to which the sub-aperture division unit 9 has performed the sub-aperture division (reconstruction algorithm). For the method, see, for example, Non-Patent Document 1 described above).
次に、上述のステップST10において、画像再生部11は、復元アルゴリズム部10が復元アルゴリズムを行った第1のチャネル信号及び第2のチャネル信号に対して、画像再生を行う。なお、画像再生方法の例として、拡張チャープスケーリング等が挙げられる(例えば、上述の非特許文献J. Mittermayer, R. Lord, and E. B¨oner,“Sliding spotlight SAR processing for TerraSAR-X using a new fomulation of the extended chirp scaling algorithm,” IGARSS2003, vol.3, p.1462-1464, 2003.を参照)。
Next, in step ST10 described above, the image reproduction unit 11 reproduces an image with respect to the first channel signal and the second channel signal to which the restoration algorithm unit 10 has performed the restoration algorithm. Examples of the image reproduction method include extended chirp scaling (for example, the above-mentioned non-patent documents J. Mittermayer, R. Lord, and E. B¨oner, “Sliding spotlight SAR processing for TerraSAR-X using a”. See new fomulation of the extended chirp scaling algorithm, ”IGARSS2003, vol.3, p.1462-1464, 2003.).
次に、上述のステップST11において、サブアパーチャー結合部12は、画像再生部11が画像再生を行った第1のチャネル信号及び第2のチャネル信号に対してサブアパーチャー結合を行うことによりSAR画像を生成する。
Next, in step ST11 described above, the sub-aperture coupling unit 12 obtains a SAR image by performing sub-aperture coupling to the first channel signal and the second channel signal for which the image reproduction unit 11 has reproduced the image. Generate.
当該具体例の構成によれば、スライディングスポットライト観測によってアジマスアンビギュイティが多重化する場合にも,マルチチャネルSARにおけるチャネル信号間のインバランスを補正することが可能である。当該具体例では,マルチチャネルSARによる観測をスライディングスポットライト観測としたが、マルチチャネルSARがフルスポットライト観測を行う場合にも適用可能である。なお、本具体例では、2回のバンドパスフィルタによる帯域制限の後段で、インバランス補正部6が位相インバランスの補正のみを行ったが、インバランス補正部6は、レンジサンプルずれの補正等のその他のインバランスの補正を行ってもよい。また、ドップラー周波数領域帯域制限部3及びヒット領域帯域制限部5の各バンドパスフィルタによる2回の帯域制限を行ったが、ドップラー周波数領域帯域制限部3によるドップラー周波数領域の帯域制限は行わず、ヒット領域帯域制限部5によってヒット領域のみ帯域制限を行ってもよい。
According to the configuration of the specific example, it is possible to correct the imbalance between the channel signals in the multi-channel SAR even when the azimuth ambiguity is multiplexed by the sliding spotlight observation. In this specific example, the observation by the multi-channel SAR is defined as the sliding spotlight observation, but it can also be applied to the case where the multi-channel SAR performs the full spotlight observation. In this specific example, the imbalance correction unit 6 only corrects the phase imbalance in the subsequent stage of the band limitation by the bandpass filter twice, but the imbalance correction unit 6 corrects the range sample deviation and the like. Other imbalance corrections may be made. Further, although the band path filters of the Doppler frequency domain band limiting unit 3 and the hit region band limiting unit 5 performed the band limitation twice, the Doppler frequency domain band limiting unit 3 did not limit the band of the Doppler frequency region. The hit area band limiting unit 5 may limit the band only in the hit area.
レーダ信号処理装置100における、電力インバランス補正部1、高速フーリエ変換部2、ドップラー周波数領域帯域制限部3、逆高速フーリエ変換部4、ヒット領域帯域制限部5、インバランス補正部6のドップラー中心周波数推定部7及び位相インバランス補正部8、サブアパーチャー分割部9、復元アルゴリズム部10、画像再生部11並びにサブアパーチャー結合部12の各機能は、処理回路により実現される。すなわち、レーダ信号処理装置100は、図2に示した各ステップの処理を実行するための処理回路を備える。この処理回路は、専用のハードウェアであってもよいが、メモリに記憶されたプログラムを実行するCPU(Central Processing Unit)であってもよい。
Doppler center of power imbalance correction unit 1, fast Fourier transform unit 2, Doppler frequency domain band limiting unit 3, inverse fast Fourier transform unit 4, hit region band limiting unit 5, and imbalance correction unit 6 in the radar signal processing device 100. Each function of the frequency estimation unit 7, the phase imbalance correction unit 8, the sub-aperture division unit 9, the restoration algorithm unit 10, the image reproduction unit 11, and the sub-aperture coupling unit 12 is realized by a processing circuit. That is, the radar signal processing device 100 includes a processing circuit for executing the processing of each step shown in FIG. This processing circuit may be dedicated hardware, or may be a CPU (Central Processing Unit) that executes a program stored in the memory.
図4Aは、レーダ信号処理装置100の機能を実現するハードウェア構成を示すブロック図である。図4Bは、レーダ信号処理装置100の機能を実現するソフトウェアを実行するハードウェア構成を示すブロック図である。
FIG. 4A is a block diagram showing a hardware configuration that realizes the functions of the radar signal processing device 100. FIG. 4B is a block diagram showing a hardware configuration for executing software that realizes the functions of the radar signal processing device 100.
上記処理回路が図4Aに示す専用のハードウェアの処理回路20である場合、処理回路20は、例えば、単一回路、複合回路、プログラム化したプロセッサ、並列プログラム化したプロセッサ、ASIC(Application Specific Integrated Circuit)、FPGA(Field-Programmable Gate Array)又はこれらを組み合わせたものが該当する。
When the processing circuit is the processing circuit 20 of the dedicated hardware shown in FIG. 4A, the processing circuit 20 may be, for example, a single circuit, a composite circuit, a programmed processor, a parallel programmed processor, or an ASIC (Application Specific Integrated Circuitd). Circuit), FPGA (Field-Programmable Gate Array) or a combination thereof is applicable.
レーダ信号処理装置100における、電力インバランス補正部1、高速フーリエ変換部2、ドップラー周波数領域帯域制限部3、逆高速フーリエ変換部4、ヒット領域帯域制限部5、インバランス補正部6のドップラー中心周波数推定部7及び位相インバランス補正部8、サブアパーチャー分割部9、復元アルゴリズム部10、画像再生部11並びにサブアパーチャー結合部12の各機能を別々の処理回路で実現してもよいし、これらの機能をまとめて1つの処理回路で実現してもよい。
Doppler center of power imbalance correction unit 1, fast Fourier transform unit 2, Doppler frequency domain band limiting unit 3, inverse fast Fourier transform unit 4, hit region band limiting unit 5, and imbalance correction unit 6 in the radar signal processing device 100. The functions of the frequency estimation unit 7, the phase imbalance correction unit 8, the sub-aperture division unit 9, the restoration algorithm unit 10, the image reproduction unit 11, and the sub-aperture coupling unit 12 may be realized by separate processing circuits, or these may be realized. The functions may be collectively realized by one processing circuit.
上記処理回路が図4Bに示すプロセッサ21である場合、レーダ信号処理装置100における、電力インバランス補正部1、高速フーリエ変換部2、ドップラー周波数領域帯域制限部3、逆高速フーリエ変換部4、ヒット領域帯域制限部5、インバランス補正部6のドップラー中心周波数推定部7及び位相インバランス補正部8、サブアパーチャー分割部9、復元アルゴリズム部10、画像再生部11並びにサブアパーチャー結合部12の各機能は、ソフトウェア、ファームウェア又はソフトウェアとファームウェアとの組み合わせによって実現される。
なお、ソフトウェア又はファームウェアは、プログラムとして記述されてメモリ22に記憶される。 When the processing circuit is theprocessor 21 shown in FIG. 4B, the power imbalance correction unit 1, the fast Fourier transform unit 2, the Doppler frequency domain band limiting unit 3, the inverse fast Fourier transform unit 4, and the hit in the radar signal processing device 100. Functions of the region band limiting unit 5, the Doppler center frequency estimation unit 7 of the imbalance correction unit 6, the phase imbalance correction unit 8, the sub-aperture division unit 9, the restoration algorithm unit 10, the image reproduction unit 11, and the sub-aperture coupling unit 12. Is realized by software, firmware or a combination of software and firmware.
The software or firmware is described as a program and stored in thememory 22.
なお、ソフトウェア又はファームウェアは、プログラムとして記述されてメモリ22に記憶される。 When the processing circuit is the
The software or firmware is described as a program and stored in the
プロセッサ21は、メモリ22に記憶されたプログラムを読み出して実行することにより、レーダ信号処理装置100における、電力インバランス補正部1、高速フーリエ変換部2、ドップラー周波数領域帯域制限部3、逆高速フーリエ変換部4、ヒット領域帯域制限部5、インバランス補正部6のドップラー中心周波数推定部7及び位相インバランス補正部8、サブアパーチャー分割部9、復元アルゴリズム部10、画像再生部11並びにサブアパーチャー結合部12の各機能を実現する。すなわち、レーダ信号処理装置100は、これらの各機能がプロセッサ21によって実行されるときに、図2に示した各ステップの処理が結果的に実行されるプログラムを記憶するためのメモリ22を備える。
By reading and executing the program stored in the memory 22, the processor 21 reads and executes the power imbalance correction unit 1, the fast Fourier transform unit 2, the Doppler frequency domain band limiting unit 3, and the inverse high-speed Fourier in the radar signal processing device 100. Conversion unit 4, hit region band limiting unit 5, Doppler center frequency estimation unit 7 and phase imbalance correction unit 8 of imbalance correction unit 6, sub-aperture division unit 9, restoration algorithm unit 10, image reproduction unit 11, and sub-aperture coupling. Each function of the part 12 is realized. That is, the radar signal processing device 100 includes a memory 22 for storing a program in which the processing of each step shown in FIG. 2 is executed as a result when each of these functions is executed by the processor 21.
これらのプログラムは、レーダ信号処理装置100における、電力インバランス補正部1、高速フーリエ変換部2、ドップラー周波数領域帯域制限部3、逆高速フーリエ変換部4、ヒット領域帯域制限部5、インバランス補正部6のドップラー中心周波数推定部7及び位相インバランス補正部8、サブアパーチャー分割部9、復元アルゴリズム部10、画像再生部11並びにサブアパーチャー結合部12の各手順又は方法をコンピュータに実行させる。メモリ22は、コンピュータを、レーダ信号処理装置100における、電力インバランス補正部1、高速フーリエ変換部2、ドップラー周波数領域帯域制限部3、逆高速フーリエ変換部4、ヒット領域帯域制限部5、インバランス補正部6のドップラー中心周波数推定部7及び位相インバランス補正部8、サブアパーチャー分割部9、復元アルゴリズム部10、画像再生部11並びにサブアパーチャー結合部12として機能させるためのプログラムが記憶されたコンピュータ可読記憶媒体であってもよい。
These programs include a power imbalance correction unit 1, a high-speed Fourier transform unit 2, a Doppler frequency domain band limiting unit 3, an inverse fast Fourier transform unit 4, a hit region band limiting unit 5, and an imbalance correction in the radar signal processing device 100. The computer is made to execute each procedure or method of the Doppler center frequency estimation unit 7 of the unit 6, the phase imbalance correction unit 8, the sub-aperture division unit 9, the restoration algorithm unit 10, the image reproduction unit 11, and the sub-aperture coupling unit 12. The memory 22 uses the computer as a power imbalance correction unit 1, a fast Fourier transform unit 2, a Doppler frequency domain band limiting unit 3, an inverse fast Fourier transform unit 4, a hit region band limiting unit 5, and an algorithm in the radar signal processing device 100. A program for functioning as the Doppler center frequency estimation unit 7 of the balance correction unit 6, the phase imbalance correction unit 8, the sub-aperture division unit 9, the restoration algorithm unit 10, the image reproduction unit 11, and the sub-aperture coupling unit 12 is stored. It may be a computer-readable storage medium.
プロセッサ21には、例えば、CPU(Central Processing Unit)、処理装置、演算装置、プロセッサ、マイクロプロセッサ、マイクロコンピュータ、またはDSP(Digital Signal Processor)などが該当する。
The processor 21 corresponds to, for example, a CPU (Central Processing Unit), a processing device, a computing device, a processor, a microprocessor, a microcomputer, a DSP (Digital Signal Processor), or the like.
メモリ22には、例えば、RAM(Random Access Memory)、ROM(Read Only Memory)、フラッシュメモリ、EPROM(Erasable Programmable Read Only Memory)、EEPROM(Electrically-EPROM)などの不揮発性又は揮発性の半導体メモリ、ハードディスク、フレキシブルディスク等の磁気ディスク、フレキシブルディスク、光ディスク、コンパクトディスク、ミニディスク、CD(Compact Disc)、DVD(Digital Versatile Disc)などが該当する。
The memory 22 may include, for example, a RAM (Random Access Memory), a ROM (Read Only Memory), a flash memory, an EPROM (Erasable Programmable Read Only Memory), an EPROM (Electrically-volatile) semiconductor, or an EPROM (Electrically-EPROM). This includes hard disks, magnetic disks such as flexible disks, flexible disks, optical discs, compact disks, mini disks, CDs (Compact Disc), DVDs (Digital Versailles Disc), and the like.
レーダ信号処理装置100における、電力インバランス補正部1、高速フーリエ変換部2、ドップラー周波数領域帯域制限部3、逆高速フーリエ変換部4、ヒット領域帯域制限部5、インバランス補正部6のドップラー中心周波数推定部7及び位相インバランス補正部8、サブアパーチャー分割部9、復元アルゴリズム部10、画像再生部11並びにサブアパーチャー結合部12の各機能について一部を専用のハードウェアで実現し、一部をソフトウェア又はファームウェアで実現してもよい。
Doppler center of power imbalance correction unit 1, fast Fourier transform unit 2, Doppler frequency domain band limiting unit 3, inverse fast Fourier transform unit 4, hit region band limiting unit 5, and imbalance correction unit 6 in the radar signal processing device 100. Some of the functions of the frequency estimation unit 7, the phase imbalance correction unit 8, the sub-aperture division unit 9, the restoration algorithm unit 10, the image reproduction unit 11, and the sub-aperture coupling unit 12 are realized by dedicated hardware, and some of them are realized. May be realized by software or firmware.
例えば、電力インバランス補正部1、高速フーリエ変換部2、ドップラー周波数領域帯域制限部3、逆高速フーリエ変換部4及びヒット領域帯域制限部5の各機能は、専用のハードウェアとしての処理回路で機能を実現する。インバランス補正部6のドップラー中心周波数推定部7及び位相インバランス補正部8、サブアパーチャー分割部9、復元アルゴリズム部10、画像再生部11並びにサブアパーチャー結合部12については、プロセッサ21がメモリ22に記憶されたプログラムを読み出して実行することにより機能を実現してもよい。
このように、処理回路は、ハードウェア、ソフトウェア、ファームウェア又はこれらの組み合わせにより上記機能のそれぞれを実現することができる。 For example, each function of the powerimbalance correction unit 1, the fast Fourier transform unit 2, the Doppler frequency domain band limiting unit 3, the inverse fast Fourier transform unit 4, and the hit region band limiting unit 5 is a processing circuit as dedicated hardware. Realize the function. For the Doppler center frequency estimation unit 7 and phase imbalance correction unit 8, sub-aperture division unit 9, restoration algorithm unit 10, image reproduction unit 11, and sub-aperture coupling unit 12 of the imbalance correction unit 6, the processor 21 is stored in the memory 22. The function may be realized by reading and executing the stored program.
As described above, the processing circuit can realize each of the above functions by hardware, software, firmware or a combination thereof.
このように、処理回路は、ハードウェア、ソフトウェア、ファームウェア又はこれらの組み合わせにより上記機能のそれぞれを実現することができる。 For example, each function of the power
As described above, the processing circuit can realize each of the above functions by hardware, software, firmware or a combination thereof.
以上のように、実施の形態1に係るレーダ信号処理装置100は、マルチチャネルSARによる観測によって得られた複数のチャネル信号に対して、それぞれ、ヒット領域における帯域制限を行うヒット領域帯域制限部5と、ヒット領域帯域制限部5がヒット領域における帯域制限を行った複数のチャネル信号に基づいて、チャネル信号間のインバランス補正を行うインバランス補正部6と、を備えている。
As described above, the radar signal processing device 100 according to the first embodiment has the hit area band limiting unit 5 that limits the band in the hit area for each of the plurality of channel signals obtained by the observation by the multi-channel SAR. And an imbalance correction unit 6 that performs imbalance correction between channel signals based on a plurality of channel signals for which the hit region band limiting unit 5 has band-limited in the hit region.
上記の構成によれば、ヒット領域における帯域制限によりアジマスアンビギュイティが分離された信号が得られる。これにより得られた信号に対してインバランス補正を行うことによって、チャネル信号間のインバランス補正の精度を向上させることができる。
According to the above configuration, a signal in which the azimuth ambiguity is separated can be obtained due to the band limitation in the hit region. By performing imbalance correction on the signal obtained thereby, the accuracy of imbalance correction between channel signals can be improved.
実施の形態1に係るレーダ信号処理装置100におけるヒット領域帯域制限部5は、ヒット領域における帯域制限によって、ヒット中心を中心とした帯域のチャネル信号を抽出する。
上記の構成によれば、ヒット中心を中心とした帯域のチャネル信号を抽出する帯域制限によりアジマスアンビギュイティが分離された信号が得られる。これにより得られた信号に対してインバランス補正を行うことによって、チャネル信号間のインバランス補正の精度を向上させることができる。 The hit regionband limiting unit 5 in the radar signal processing device 100 according to the first embodiment extracts a channel signal in a band centered on the hit center by band limiting in the hit region.
According to the above configuration, a signal in which the azimuth ambiguity is separated can be obtained by band limitation for extracting the channel signal in the band centered on the hit center. By performing imbalance correction on the signal obtained thereby, the accuracy of imbalance correction between channel signals can be improved.
上記の構成によれば、ヒット中心を中心とした帯域のチャネル信号を抽出する帯域制限によりアジマスアンビギュイティが分離された信号が得られる。これにより得られた信号に対してインバランス補正を行うことによって、チャネル信号間のインバランス補正の精度を向上させることができる。 The hit region
According to the above configuration, a signal in which the azimuth ambiguity is separated can be obtained by band limitation for extracting the channel signal in the band centered on the hit center. By performing imbalance correction on the signal obtained thereby, the accuracy of imbalance correction between channel signals can be improved.
実施の形態1に係るレーダ信号処理装置100は、複数のチャネル信号に対して、それぞれ、ドップラー周波数変換を行う高速フーリエ変換部2と、高速フーリエ変換部2がドップラー周波数変換を行った複数のチャネル信号に対して、それぞれ、ドップラー周波数領域における帯域制限を行うドップラー周波数領域帯域制限部3と、ドップラー周波数領域帯域制限部3がドップラー周波数領域における帯域制限を行った複数のチャネル信号に対して、それぞれ、逆ドップラー周波数変換を行う逆高速フーリエ変換部4と、をさらに備え、ヒット領域帯域制限部5は、逆高速フーリエ変換部4が逆ドップラー周波数変換を行った複数のチャネル信号に対して、それぞれ、ヒット領域における帯域制限を行う。
The radar signal processing device 100 according to the first embodiment has a high-speed Fourier conversion unit 2 that performs Doppler frequency conversion on a plurality of channel signals, and a plurality of channels that the high-speed Fourier conversion unit 2 performs Doppler frequency conversion. For each of the signals, the Doppler frequency domain band limiting unit 3 that limits the band in the Doppler frequency domain and the plurality of channel signals that the Doppler frequency domain band limiting section 3 band limits in the Doppler frequency domain, respectively. , Inverse high-speed Fourier conversion unit 4 that performs inverse Doppler frequency conversion, and the hit region band limiting unit 5 for each of a plurality of channel signals that the inverse high-speed Fourier conversion unit 4 has performed inverse Doppler frequency conversion. , Band limitation in the hit area.
上記の構成によれば、ヒット領域における帯域制限、及びドップラー周波数領域における帯域制限によりアジマスアンビギュイティが分離された信号が得られる。これにより得られた信号に対してインバランス補正を行うことによって、チャネル信号間のインバランス補正の精度を向上させることができる。
According to the above configuration, a signal in which the azimuth ambiguity is separated can be obtained due to the band limitation in the hit region and the band limitation in the Doppler frequency region. By performing imbalance correction on the signal obtained thereby, the accuracy of imbalance correction between channel signals can be improved.
実施の形態1に係るレーダ信号処理装置100におけるドップラー周波数領域帯域制限部3は、ドップラー周波数領域における帯域制限によって、ドップラー中心周波数を中心とした帯域のチャネル信号を抽出する。
The Doppler frequency domain band limiting unit 3 in the radar signal processing device 100 according to the first embodiment extracts a channel signal in a band centered on the Doppler center frequency by band limiting in the Doppler frequency domain.
上記の構成によれば、ドップラー中心周波数を中心とした帯域のチャネル信号を抽出する帯域制限によりアジマスアンビギュイティが分離された信号が得られる。これにより得られた信号に対してインバランス補正を行うことによって、チャネル信号間のインバランス補正の精度を向上させることができる。
According to the above configuration, a signal in which the azimuth ambiguity is separated can be obtained by band limitation for extracting the channel signal in the band centered on the Doppler center frequency. By performing imbalance correction on the signal obtained thereby, the accuracy of imbalance correction between channel signals can be improved.
実施の形態1に係るレーダ信号処理装置100におけるインバランス補正部6は、インバランス補正として、位相インバランスの補正を行う。
上記の構成によれば、帯域制限によりアジマスアンビギュイティが分離された信号に対して位相インバランス補正を行うことによって、チャネル信号間の位相インバランス補正の精度を向上させることができる。 Theimbalance correction unit 6 in the radar signal processing device 100 according to the first embodiment corrects the phase imbalance as the imbalance correction.
According to the above configuration, the accuracy of the phase imbalance correction between the channel signals can be improved by performing the phase imbalance correction on the signal whose azimuth ambiguity is separated by the band limitation.
上記の構成によれば、帯域制限によりアジマスアンビギュイティが分離された信号に対して位相インバランス補正を行うことによって、チャネル信号間の位相インバランス補正の精度を向上させることができる。 The
According to the above configuration, the accuracy of the phase imbalance correction between the channel signals can be improved by performing the phase imbalance correction on the signal whose azimuth ambiguity is separated by the band limitation.
実施の形態1に係るレーダ信号処理装置100におけるインバランス補正部6は、インバランス補正として、レンジサンプルずれの補正を行ってもよい。
上記の構成によれば、帯域制限によりアジマスアンビギュイティが分離された信号に対してレンジサンプルずれの補正を行うことによって、チャネル信号間のレンジサンプルずれの補正精度を向上させることができる。 Theimbalance correction unit 6 in the radar signal processing device 100 according to the first embodiment may correct the range sample deviation as the imbalance correction.
According to the above configuration, it is possible to improve the correction accuracy of the range sample deviation between the channel signals by correcting the range sample deviation for the signal whose azimuth ambiguity is separated by the band limitation.
上記の構成によれば、帯域制限によりアジマスアンビギュイティが分離された信号に対してレンジサンプルずれの補正を行うことによって、チャネル信号間のレンジサンプルずれの補正精度を向上させることができる。 The
According to the above configuration, it is possible to improve the correction accuracy of the range sample deviation between the channel signals by correcting the range sample deviation for the signal whose azimuth ambiguity is separated by the band limitation.
実施の形態1に係るレーダ信号処理装置100におけるマルチチャネルSARによる観測は、スライディングスポット観測、又はフルスポットライト観測である。
上記の構成によれば、スライディングスポット観測、又はフルスポットライト観測によってアジマスアンビギュイティが多重化した場合に、チャネル信号間のインバランス補正の精度を向上させることができる。 The observation by the multi-channel SAR in the radarsignal processing apparatus 100 according to the first embodiment is a sliding spot observation or a full spotlight observation.
According to the above configuration, when the azimuth ambiguity is multiplexed by sliding spot observation or full spotlight observation, the accuracy of imbalance correction between channel signals can be improved.
上記の構成によれば、スライディングスポット観測、又はフルスポットライト観測によってアジマスアンビギュイティが多重化した場合に、チャネル信号間のインバランス補正の精度を向上させることができる。 The observation by the multi-channel SAR in the radar
According to the above configuration, when the azimuth ambiguity is multiplexed by sliding spot observation or full spotlight observation, the accuracy of imbalance correction between channel signals can be improved.
実施の形態1に係るレーダ信号処理方法は、マルチチャネルSARによる観測によって得られた複数のチャネル信号に対して、それぞれ、ヒット領域における帯域制限を行うヒット領域帯域制限ステップと、ヒット領域帯域制限ステップでヒット領域における帯域制限を行った複数のチャネル信号に基づいて、チャネル信号間のインバランス補正を行うインバランス補正ステップと、を含む。
上記の構成によれば、実施の形態1に係るレーダ信号処理装置100が奏する効果と同様の効果を奏する。
なお、実施の形態の任意の構成要素の変形、もしくは実施の形態の任意の構成要素の省略が可能である。 In the radar signal processing method according to the first embodiment, a hit area band limiting step and a hit area band limiting step for band limiting in the hit region are performed for a plurality of channel signals obtained by observation by the multi-channel SAR, respectively. Includes an imbalance correction step of performing imbalance correction between channel signals based on a plurality of channel signals for which band limitation has been performed in the hit region.
According to the above configuration, the same effect as that of the radarsignal processing device 100 according to the first embodiment is obtained.
It is possible to modify any component of the embodiment or omit any component of the embodiment.
上記の構成によれば、実施の形態1に係るレーダ信号処理装置100が奏する効果と同様の効果を奏する。
なお、実施の形態の任意の構成要素の変形、もしくは実施の形態の任意の構成要素の省略が可能である。 In the radar signal processing method according to the first embodiment, a hit area band limiting step and a hit area band limiting step for band limiting in the hit region are performed for a plurality of channel signals obtained by observation by the multi-channel SAR, respectively. Includes an imbalance correction step of performing imbalance correction between channel signals based on a plurality of channel signals for which band limitation has been performed in the hit region.
According to the above configuration, the same effect as that of the radar
It is possible to modify any component of the embodiment or omit any component of the embodiment.
本開示に係るレーダ処理装置は、チャネル信号間のインバランス補正の精度を向上させることができるため、マルチチャネルSARの技術に利用可能である。
The radar processing device according to the present disclosure can be used in the multi-channel SAR technique because it can improve the accuracy of imbalance correction between channel signals.
1 電力インバランス補正部、2 高速フーリエ変換部、3 ドップラー周波数領域帯域制限部、4 逆高速フーリエ変換部、5 ヒット領域帯域制限部、6 インバランス補正部、7 ドップラー中心周波数推定部、8 位相インバランス補正部、9 サブアパーチャー分割部、10 復元アルゴリズム部、11 画像再生部、12 サブアパーチャー結合部、20 処理回路、21 プロセッサ、22 メモリ、100 レーダ信号処理装置。
1 power imbalance correction unit, 2 high-speed Fourier transform unit, 3 Doppler frequency domain band limiting unit, 4 inverse fast Fourier transform unit, 5 hit region band limiting unit, 6 imbalance correction unit, 7 Doppler center frequency estimation unit, 8 phase Imbalance correction unit, 9 sub-aperture division unit, 10 restoration algorithm unit, 11 image reproduction unit, 12 sub-aperture coupling unit, 20 processing circuit, 21 processor, 22 memory, 100 radar signal processing device.
Claims (8)
- マルチチャネルSARによる観測によって得られた複数のチャネル信号に対して、それぞれ、ヒット領域における帯域制限を行うヒット領域帯域制限部と、
前記ヒット領域帯域制限部が前記ヒット領域における帯域制限を行った複数のチャネル信号に基づいて、チャネル信号間のインバランス補正を行うインバランス補正部と、を備えていることを特徴とする、レーダ信号処理装置。 A hit area band limiting unit that limits the band in the hit area for a plurality of channel signals obtained by observation by the multi-channel SAR, respectively.
The radar is characterized in that the hit region band limiting unit includes an imbalance correction unit that performs imbalance correction between channel signals based on a plurality of channel signals for which band limitation is performed in the hit region. Signal processing device. - 前記ヒット領域帯域制限部は、前記ヒット領域における帯域制限によって、ヒット中心を中心とした帯域のチャネル信号を抽出することを特徴とする、請求項1に記載のレーダ信号処理装置。 The radar signal processing device according to claim 1, wherein the hit region band limiting unit extracts a channel signal in a band centered on the hit center by band limiting in the hit region.
- 前記複数のチャネル信号に対して、それぞれ、ドップラー周波数変換を行う高速フーリエ変換部と、
前記高速フーリエ変換部が前記ドップラー周波数変換を行った複数のチャネル信号に対して、それぞれ、ドップラー周波数領域における帯域制限を行うドップラー周波数領域帯域制限部と、
前記ドップラー周波数領域帯域制限部が前記ドップラー周波数領域における帯域制限を行った複数のチャネル信号に対して、それぞれ、逆ドップラー周波数変換を行う逆高速フーリエ変換部と、をさらに備え、
前記ヒット領域帯域制限部は、前記逆高速フーリエ変換部が前記逆ドップラー周波数変換を行った複数のチャネル信号に対して、それぞれ、前記ヒット領域における帯域制限を行うことを特徴とする、請求項1に記載のレーダ信号処理装置。 A fast Fourier transform unit that performs Doppler frequency conversion for each of the plurality of channel signals,
A Doppler frequency domain band limiting section that limits the band in the Doppler frequency domain for a plurality of channel signals that have undergone the Doppler frequency transform by the fast Fourier transform section, respectively.
The Doppler frequency domain band limiting section further includes an inverse fast Fourier transform section that performs inverse Doppler frequency conversion for each of a plurality of channel signals for which band limiting has been performed in the Doppler frequency domain.
The hit region band limiting unit is characterized in that the inverse fast Fourier transform unit performs band limitation in the hit region for each of a plurality of channel signals subjected to the inverse Doppler frequency conversion. The radar signal processing device described in. - 前記ドップラー周波数領域帯域制限部は、前記ドップラー周波数領域における帯域制限によって、ドップラー中心周波数を中心とした帯域のチャネル信号を抽出することを特徴とする、請求項3に記載のレーダ信号処理装置。 The radar signal processing device according to claim 3, wherein the Doppler frequency domain band limiting unit extracts a channel signal in a band centered on the Doppler center frequency by band limiting in the Doppler frequency domain.
- 前記インバランス補正部は、前記インバランス補正として、位相インバランスの補正を行うことを特徴とする、請求項1に記載のレーダ信号処理装置。 The radar signal processing device according to claim 1, wherein the imbalance correction unit corrects a phase imbalance as the imbalance correction.
- 前記インバランス補正部は、前記インバランス補正として、レンジサンプルずれの補正を行うことを特徴とする、請求項1に記載のレーダ信号処理装置。 The radar signal processing device according to claim 1, wherein the imbalance correction unit corrects a range sample deviation as the imbalance correction.
- 前記観測は、スライディングスポット観測、又はフルスポットライト観測であることを特徴とする、請求項1に記載のレーダ信号処理装置。 The radar signal processing device according to claim 1, wherein the observation is a sliding spot observation or a full spotlight observation.
- マルチチャネルSARによる観測によって得られた複数のチャネル信号に対して、それぞれ、ヒット領域における帯域制限を行うヒット領域帯域制限ステップと、
前記ヒット領域帯域制限ステップで前記ヒット領域における帯域制限を行った複数のチャネル信号に基づいて、チャネル信号間のインバランス補正を行うインバランス補正ステップと、を含むことを特徴とする、レーダ信号処理方法。 A hit region bandwidth limiting step that limits the bandwidth in the hit region for a plurality of channel signals obtained by observation by the multi-channel SAR, respectively.
Radar signal processing comprising an imbalance correction step of performing imbalance correction between channel signals based on a plurality of channel signals for which band limitation is performed in the hit region in the hit area band limiting step. Method.
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2022556832A JP7214061B2 (en) | 2020-10-19 | 2020-10-19 | Radar signal processing device and radar signal processing method |
PCT/JP2020/039224 WO2022085044A1 (en) | 2020-10-19 | 2020-10-19 | Radar signal processing device and radar signal processing method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/JP2020/039224 WO2022085044A1 (en) | 2020-10-19 | 2020-10-19 | Radar signal processing device and radar signal processing method |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2022085044A1 true WO2022085044A1 (en) | 2022-04-28 |
Family
ID=81290261
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2020/039224 WO2022085044A1 (en) | 2020-10-19 | 2020-10-19 | Radar signal processing device and radar signal processing method |
Country Status (2)
Country | Link |
---|---|
JP (1) | JP7214061B2 (en) |
WO (1) | WO2022085044A1 (en) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006343290A (en) * | 2005-06-10 | 2006-12-21 | Mitsubishi Electric Corp | Imaging radar device |
JP2015001428A (en) * | 2013-06-14 | 2015-01-05 | 三菱電機株式会社 | Radar signal processor and radar signal processing method |
WO2016148104A1 (en) * | 2015-03-16 | 2016-09-22 | 三菱電機株式会社 | Radar moving image creation device and method |
-
2020
- 2020-10-19 JP JP2022556832A patent/JP7214061B2/en active Active
- 2020-10-19 WO PCT/JP2020/039224 patent/WO2022085044A1/en active Application Filing
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006343290A (en) * | 2005-06-10 | 2006-12-21 | Mitsubishi Electric Corp | Imaging radar device |
JP2015001428A (en) * | 2013-06-14 | 2015-01-05 | 三菱電機株式会社 | Radar signal processor and radar signal processing method |
WO2016148104A1 (en) * | 2015-03-16 | 2016-09-22 | 三菱電機株式会社 | Radar moving image creation device and method |
Non-Patent Citations (3)
Title |
---|
FAN, HUAITAO ET AL.: "Phase Mismatch Calibration for Multichannel Sliding Spotlight SAR Imaging with Extended Azimuth Cross Correlation", IGARSS 2019-2019 IEEE INTERNATIONAL GEOSCIENCE AND REMOTE SENSING SYMPOSIUM, August 2019 (2019-08-01), pages 2531 - 2534, XP033656916, DOI: 10.1109/IGARSS.2019.8900075 * |
WU, DI ET AL.: "A Channel Calibration Algorithm Based on Isolated Scatterers for Multi-Channel HRWS-SAR", IEEE ACCESS, vol. 7, 2019, pages 135665 - 135677, XP011747900, DOI: 10.1109/ACCESS.2019.2941203 * |
ZHOU, YASHI ET AL.: "A Novel Approach to Doppler Centroid and Channel Errors Estimation in Azimuth Multi-Channel SAR", IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING, vol. 57, no. 11, November 2019 (2019-11-01), pages 8430 - 8444, XP011755607, DOI: 10.1106/TGRS.2019.2921094 * |
Also Published As
Publication number | Publication date |
---|---|
JP7214061B2 (en) | 2023-01-27 |
JPWO2022085044A1 (en) | 2022-04-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109085590B (en) | Broadband direct acquisition data ISAR imaging method based on ARP (Address resolution protocol) segmented coherent accumulation | |
CN108387900B (en) | Vibration error compensation method for helicopter-mounted rotary synthetic aperture radar | |
CN108051812B (en) | Satellite-borne SAR moving target detection method based on two-dimensional speed search | |
CN108279404B (en) | Two-channel SAR phase error correction method based on spatial spectrum estimation | |
KR20200010816A (en) | Apparatus and method for processing radar signal for a multi channel antenna | |
CN113030962B (en) | Airborne terahertz synthetic aperture radar and imaging method | |
Liu et al. | A novel channel phase bias estimation method for spaceborne along-track multi-channel HRWS SAR in time-domain | |
Feng et al. | An improved scheme of digital beam-forming in elevation for spaceborne SAR | |
WO2012006226A1 (en) | Processing radar return signals to detect targets | |
JP5489813B2 (en) | Radar image processing device | |
JP6246338B2 (en) | Angle measuring device and angle measuring method | |
JP2010127771A (en) | Synthetic aperture sonar, and method and program for correcting phase error of synthetic aperture sonar | |
WO2022085044A1 (en) | Radar signal processing device and radar signal processing method | |
JP2015532848A (en) | Image tissue motion estimation | |
CN112379379A (en) | Channel phase deviation estimation method and device of TOPS SAR system | |
CN105974413B (en) | The self-focusing method of more base external illuminators-based radar imaging systems | |
CN109725313B (en) | SAR (synthetic aperture radar) sea wave imaging method and system, electronic equipment and medium | |
WO2014084933A1 (en) | Enhanced radar range resolution | |
JP2009180650A (en) | Radar system | |
WO2021181590A1 (en) | Calibration value calculation device, radar calibration device, and calibration value calculation method | |
JP7123670B2 (en) | Radar system and signal processing method | |
JP7499995B1 (en) | Radar signal processing device, radar signal processing method, and radar device | |
Yang et al. | An auto-focusing algorithm for monopulse imaging technique | |
JP7399368B2 (en) | Radar signal processing device, radar signal processing method, and target observation system | |
JP6735632B2 (en) | Signal processing device, signal processing method, and signal processing program |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 20958600 Country of ref document: EP Kind code of ref document: A1 |
|
ENP | Entry into the national phase |
Ref document number: 2022556832 Country of ref document: JP Kind code of ref document: A |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 20958600 Country of ref document: EP Kind code of ref document: A1 |