CN116930969A - Metamaterial tag target positioning and imaging method based on synthetic aperture radar - Google Patents
Metamaterial tag target positioning and imaging method based on synthetic aperture radar Download PDFInfo
- Publication number
- CN116930969A CN116930969A CN202311175498.0A CN202311175498A CN116930969A CN 116930969 A CN116930969 A CN 116930969A CN 202311175498 A CN202311175498 A CN 202311175498A CN 116930969 A CN116930969 A CN 116930969A
- Authority
- CN
- China
- Prior art keywords
- frequency
- modulation
- signal
- tag
- metamaterial
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 34
- 230000006870 function Effects 0.000 claims abstract description 16
- 238000001514 detection method Methods 0.000 claims abstract description 10
- 238000000034 method Methods 0.000 claims abstract description 9
- 238000001228 spectrum Methods 0.000 claims description 13
- 230000009466 transformation Effects 0.000 claims description 8
- 230000006835 compression Effects 0.000 claims description 5
- 238000007906 compression Methods 0.000 claims description 5
- 230000033001 locomotion Effects 0.000 claims description 5
- 238000012937 correction Methods 0.000 claims description 4
- 230000005012 migration Effects 0.000 claims description 4
- 238000013508 migration Methods 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims 3
- 230000008685 targeting Effects 0.000 claims 3
- 230000001427 coherent effect Effects 0.000 claims 1
- 238000001914 filtration Methods 0.000 claims 1
- 230000010355 oscillation Effects 0.000 claims 1
- 230000001131 transforming effect Effects 0.000 claims 1
- 238000013461 design Methods 0.000 abstract description 4
- 230000008447 perception Effects 0.000 abstract description 4
- 238000005516 engineering process Methods 0.000 description 10
- 230000000694 effects Effects 0.000 description 5
- 238000004891 communication Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 206010047513 Vision blurred Diseases 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/904—SAR modes
- G01S13/9052—Spotlight mode
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/904—SAR modes
- G01S13/9047—Doppler beam sharpening mode
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
本发明公开了一种基于合成孔径雷达的超材料标签目标定位与成像方法,属于目标成像与探测领域,采用调频连续波合成孔径雷达(SAR)系统作为探测主机,发射点频与线性调频信号的联合波形,将基于超材料调制编码方法实现的标签设计等效为具有对雷达探测信号调制功能的点目标。利用点频段的回波信号解算超材料点目标的调制频率,并构建补偿因子,对波数域成像算法进行距离向与方位向二维调制补偿,实现调制目标的聚焦成像。本发明可在标签识别与定位的同时,实现标签周围环境的成像与感知。
The invention discloses a metamaterial tag target positioning and imaging method based on synthetic aperture radar, which belongs to the field of target imaging and detection. It uses a frequency modulated continuous wave synthetic aperture radar (SAR) system as the detection host to emit point frequency and linear frequency modulation signals. With the combined waveform, the tag design based on the metamaterial modulation coding method is equivalent to a point target with the function of modulating the radar detection signal. The echo signal in the point frequency band is used to calculate the modulation frequency of the metamaterial point target, and a compensation factor is constructed to perform two-dimensional modulation compensation in the range and azimuth directions for the wavenumber domain imaging algorithm to achieve focused imaging of the modulated target. The invention can realize imaging and perception of the surrounding environment of the tag while identifying and positioning the tag.
Description
技术领域Technical field
本发明属于目标成像与探测领域,具体涉及一种基于合成孔径雷达的超材料标签目标定位与成像方法。The invention belongs to the field of target imaging and detection, and specifically relates to a metamaterial tag target positioning and imaging method based on synthetic aperture radar.
背景技术Background technique
标签技术通过标签端的信息编码与主机端的信号检测完成信息交互。一方面通过信息编码实现标签的唯一性识别,另一方面通过信号检测与通信时间延迟,实现标签的发现与定位。通过对携带自定义标签,可实现人、物体或设备的识别与定位,进而为人类的生活、生产、管理等提供更好的保障。Tag technology completes information interaction through information encoding on the tag side and signal detection on the host side. On the one hand, the unique identification of the tag is achieved through information encoding, and on the other hand, the discovery and positioning of the tag is achieved through signal detection and communication time delay. By carrying custom tags, people, objects or equipment can be identified and positioned, thereby providing better protection for human life, production, management, etc.
传统RFID标签技术,主机通过读取标签的编码信息可实现对个体的准确识别,这种方法技术成熟,标签具有廉价、低功耗、重量轻的优势,但存在作用距离不足、定位精度差等问题。另一代表性技术为超宽带无线通信技术(UWB),UWB是通过主机与标签之间通信完成标签的识别与定位,其具有低功耗、低成本的优势,但在定位精度、刷新率方面也存在不足。With traditional RFID tag technology, the host can accurately identify individuals by reading the coded information of the tag. This method is technically mature and the tags have the advantages of cheapness, low power consumption, and light weight. However, they have insufficient range and poor positioning accuracy. question. Another representative technology is ultra-wideband wireless communication technology (UWB). UWB completes the identification and positioning of tags through communication between the host and the tag. It has the advantages of low power consumption and low cost, but it suffers from disadvantages in terms of positioning accuracy and refresh rate. There are also shortcomings.
合成孔径雷达技术可实现远距离、高分辨二维成像,将合成孔径雷达技术与标签技术相结合,有望在实现标签的远距离、高精度、高刷新率定位的同时,实现对异步调制标签的聚焦成像。相比传统标签技术,可极大提升标签周围环境的成像与感知能力,对复杂场景的目标高精度定位与环境感知具有重要意义。Synthetic aperture radar technology can achieve long-distance, high-resolution two-dimensional imaging. Combining synthetic aperture radar technology with tag technology is expected to achieve long-distance, high-precision, and high-refresh rate positioning of tags while also enabling asynchronous modulation of tags. Focused imaging. Compared with traditional tag technology, it can greatly improve the imaging and perception capabilities of the environment around the tag, which is of great significance for high-precision target positioning and environmental perception in complex scenes.
发明内容Contents of the invention
为解决上述技术问题,本发明提供一种基于合成孔径雷达的超材料标签目标定位与成像方法,采用调频连续波合成孔径雷达(SAR)系统作为探测主机,发射点频与线性调频信号的联合波形,将基于超材料调制编码方法实现的标签设计等效为具有对雷达探测信号调制功能的点目标。利用点频段的回波信号解算超材料点目标的调制频率,并构建补偿因子,对波数域成像算法进行距离向与方位向二维调制补偿,实现调制目标的聚焦成像。In order to solve the above technical problems, the present invention provides a metamaterial tag target positioning and imaging method based on synthetic aperture radar, which uses a frequency modulated continuous wave synthetic aperture radar (SAR) system as the detection host to emit a combined waveform of point frequency and linear frequency modulation signals. , the tag design based on the metamaterial modulation coding method is equivalent to a point target with the function of modulating the radar detection signal. The echo signal in the point frequency band is used to calculate the modulation frequency of the metamaterial point target, and a compensation factor is constructed to perform two-dimensional modulation compensation in the range and azimuth directions for the wavenumber domain imaging algorithm to achieve focused imaging of the modulated target.
为达到上述目的,本发明采用如下技术方案:In order to achieve the above objects, the present invention adopts the following technical solutions:
一种基于合成孔径雷达的超材料标签目标定位与成像方法,包括如下步骤:A metamaterial tag target positioning and imaging method based on synthetic aperture radar, including the following steps:
步骤1、构建点频信号与线性调频信号的联合波形;Step 1. Construct the joint waveform of the point frequency signal and the linear frequency modulation signal;
步骤2、线性调频段回波信号相干接收后,根据去调频接收原理构建超材料标签调制的回波信号模型;Step 2. After the linear frequency modulation band echo signal is coherently received, a metamaterial tag modulated echo signal model is constructed based on the de-FM reception principle;
步骤3、将点频段回波信号变换到二维频域,根据超材料标签调制的回波信号模型解算目标方位调制频率;Step 3. Transform the point frequency band echo signal into the two-dimensional frequency domain, and calculate the target azimuth modulation frequency based on the echo signal model modulated by the metamaterial tag;
步骤4、利用超材料标签的调制参数和步骤3中解算出的目标方位调制频率分别对线性调频段回波信号进行距离向和方位向的调制补偿,再通过构建一致性压缩的参考函数与Stolt插值实现全局徙动校正与方位压缩,最后通过二维逆傅里叶变化得到聚焦图像结果。Step 4: Use the modulation parameters of the metamaterial tag and the target azimuth modulation frequency calculated in step 3 to perform modulation compensation in the range and azimuth direction of the chirp band echo signal, and then construct a consistent compression reference function and Stolt Interpolation realizes global migration correction and azimuth compression, and finally the focused image result is obtained through two-dimensional inverse Fourier transformation.
有益效果:Beneficial effects:
本发明提出的方法结合了合成孔径雷达技术和超材料标签调制的特点,可实现对异步调制标签的定位与聚焦成像,与传统标签技术相比,可在标签识别与定位的同时,实现标签周围环境的成像与感知。The method proposed by the present invention combines the characteristics of synthetic aperture radar technology and metamaterial tag modulation, and can realize the positioning and focused imaging of asynchronously modulated tags. Compared with traditional tag technology, it can realize tag identification and positioning at the same time as tag surroundings. Imaging and perception of the environment.
附图说明Description of the drawings
图1为点频信号与线性调频信号联合波形图;Figure 1 is a joint waveform diagram of point frequency signal and linear frequency modulation signal;
图2 为本发明的基于合成孔径雷达的超材料调制标签目标定位与成像方法的流程图。Figure 2 is a flow chart of the synthetic aperture radar-based metamaterial modulation tag target positioning and imaging method of the present invention.
具体实施方式Detailed ways
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。In order to make the purpose, technical solutions and advantages of the present invention more clear, the present invention will be further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described here are only used to explain the present invention and are not intended to limit the present invention. In addition, the technical features involved in the various embodiments of the present invention described below can be combined with each other as long as they do not conflict with each other.
如图2所示,本发明的基于合成孔径雷达的超材料标签目标定位与成像方法包括如下步骤:As shown in Figure 2, the metamaterial tag target positioning and imaging method based on synthetic aperture radar of the present invention includes the following steps:
步骤1、构建点频信号与线性调频信号的联合波形;Step 1. Construct the joint waveform of the point frequency signal and the linear frequency modulation signal;
步骤2、线性调频段回波信号相干接收后,根据去调频接收原理构建超材料标签调制的回波信号模型;Step 2. After the linear frequency modulation band echo signal is coherently received, a metamaterial tag modulated echo signal model is constructed based on the de-FM reception principle;
步骤3、将点频段回波信号变换到二维频域,根据超材料标签调制的回波信号模型解算目标方位调制频率;Step 3. Transform the point frequency band echo signal into the two-dimensional frequency domain, and calculate the target azimuth modulation frequency based on the echo signal model modulated by the metamaterial tag;
步骤4、利用超材料标签的调制参数和步骤3中解算出的目标方位调制频率分别对线性调频段回波信号进行距离向和方位向的调制补偿,再通过构建一致性压缩的参考函数与Stolt插值实现全局徙动校正与方位压缩,最后通过二维逆傅里叶变化得到聚焦图像结果。Step 4: Use the modulation parameters of the metamaterial tag and the target azimuth modulation frequency calculated in step 3 to perform modulation compensation in the range and azimuth direction of the chirp band echo signal, and then construct a consistent compression reference function and Stolt Interpolation realizes global migration correction and azimuth compression, and finally the focused image result is obtained through two-dimensional inverse Fourier transformation.
具体地,所述步骤1 包括:Specifically, step 1 includes:
超材料调制标签在SAR回波的距离向和方位向都将产生调制作用,在聚焦成像前需要消除调制效应的影响。距离向的调制效应可以根据超材料调制标签的调制参数予以消除,方位向的调制效应需要专门的发射波形设计以提取方位调制参数。为此,本发明设计了一种由点频和线性调频组成的联合波形,如图1所示,上图为幅度随时间变化关系图,下图为频率随时间变化关系图,波形重复周期为,在一个周期内,点频信号时间长度为/>,线性调频信号时间长度为/>,/>为/>时间内点频信号频率,同时也表示线性调频段的起始频率,/>为线性调频段的截止频率,/>为线性调频率利用点频信号回波解算方位调制频率,利用线性调频信号实现目标成像。The metamaterial modulation tag will produce modulation effects in both the range and azimuth directions of the SAR echo, and the influence of the modulation effect needs to be eliminated before focusing imaging. The modulation effect in the range direction can be eliminated according to the modulation parameters of the metamaterial modulation tag, while the modulation effect in the azimuth direction requires special transmission waveform design to extract the azimuth modulation parameters. To this end, the present invention designs a joint waveform composed of point frequency and linear frequency modulation, as shown in Figure 1. The upper picture shows the relationship between amplitude and time, and the lower picture shows the relationship between frequency and time. The waveform repetition period is , in one cycle, the time length of the point frequency signal is/> , the time length of the chirp signal is/> ,/> for/> The frequency of the dot frequency signal within the time period also represents the starting frequency of the chirp segment,/> is the cutoff frequency of the chirp band,/> For the linear frequency modulation frequency, the point frequency signal echo is used to solve the azimuth modulation frequency, and the linear frequency modulation signal is used to achieve target imaging.
(1) (1)
其中,为发射信号、/>、/>分别为点频信号和线性调频信号幅度、/>为载频、/>为发射信号传播时间、/>为平台运动时间、/>为线性调频率、/>为自然数、exp[]代表e的指数幂运算,j表示虚数、/>为点频信号距离向包络函数、为线性调频信号距离向包络函数、rect()为矩形函数(对任意变量/>,若/>,/>,若/>,/>)、为方位向包络函数、/>为斜距平面内测得的目标视角、/>为波束宽度、/>为脉冲函数,/>为不等于零的任意变量。in, To transmit signals,/> ,/> are the amplitude of the point frequency signal and linear frequency modulation signal respectively,/> is the carrier frequency,/> is the propagation time of the transmitted signal,/> For the platform movement time,/> is the chirp frequency,/> is a natural number, exp[] represents the exponential power operation of e, j represents an imaginary number,/> is the range envelope function of the point frequency signal, is the range envelope function of the chirp signal, and rect() is the rectangular function (for any variable/> , if/> ,/> , if/> ,/> ), is the azimuth envelope function,/> is the target angle of view measured in the slant distance plane,/> is the beam width,/> is the pulse function,/> is any variable not equal to zero.
具体地,所述步骤2包括:Specifically, the step 2 includes:
本发明采用调制编码方法实现标签,其可等效为具有对雷达探测信号调制功能的点目标,设调制频率为、第n个调制信号的初始相位为/>,n取正整数,调制方式为相位调制。根据本步骤1中的联合波形,线性调频段回波信号/>表示为:The present invention uses a modulation coding method to implement tags, which can be equivalent to a point target with the function of modulating radar detection signals. Suppose the modulation frequency is , the initial phase of the nth modulated signal is/> , n takes a positive integer, and the modulation method is phase modulation. According to the joint waveform in step 1 of this step, the chirp band echo signal/> Expressed as:
(2) (2)
其中,为回波幅度参数、/>表示回波的延时,in, is the echo amplitude parameter,/> Represents the delay of the echo,
为雷达与目标的实时斜距,目标在场景中坐标为/>,/>表示最小斜距,/>为目标方位位置,/>为平台运动速度,全时间,/>为斜视角,/>电磁波的在空气中的传播速度。 is the real-time slant distance between the radar and the target. The coordinates of the target in the scene are/> ,/> Indicates the minimum slope distance,/> is the target azimuth position,/> is the platform movement speed, full time ,/> For oblique angle,/> The propagation speed of electromagnetic waves in the air.
采用发射信号作为去斜接收的本振信号,去斜后的中频信号为,其表达式为:The transmitted signal is used as the local oscillator signal for deskewing reception. The deskewed intermediate frequency signal is , its expression is:
(3) (3)
其中,为去斜接收的中频信号幅度参数。in, It is the amplitude parameter of the intermediate frequency signal received by deskewing.
对公式(3)进行距离向与方位向的二维傅里叶变换得到接收中频回波的二维频域信号,其表达式为:Perform two-dimensional Fourier transform on the range and azimuth directions of equation (3) to obtain the two-dimensional frequency domain signal of the received intermediate frequency echo. , its expression is:
(4) (4)
其中,为表征目标距离的频率值、/>表示点目标的在零多普勒时刻的距离,/>为距离向频率轴变量,/>为方位向频率轴变量,/>为二维频域信号的幅度参数。in, is the frequency value that represents the target distance,/> Represents the distance of a point target at zero Doppler time,/> is the distance-to-frequency axis variable,/> is the azimuth frequency axis variable,/> is the amplitude parameter of the two-dimensional frequency domain signal.
为方位频谱包络,/>为最小斜距,/>为平台运动速度,/>为点目标调制导致的方位向频谱偏移量,mod为取余函数。公式(4)表明,标签调制信号使得回波信号在距离频域产生/>的偏移,在方位频域产生/>的偏移,将导致成像处理无法聚焦,在进行成像处理前需要解算出偏移量并进行补偿。 is the azimuth spectrum envelope,/> is the minimum slope distance,/> is the movement speed of the platform,/> is the azimuth spectrum offset caused by point target modulation, and mod is the remainder function. Equation (4) shows that the tag modulation signal causes the echo signal to be generated in the range frequency domain/> offset, generated in the azimuth frequency domain/> The offset will cause the imaging process to be unable to focus, and the offset needs to be calculated and compensated before imaging processing.
具体地,所述步骤3包括:Specifically, step 3 includes:
根据步骤1的联合波形,采用去斜接收后得到的点频信号段回波信号表示为:According to the joint waveform in step 1, the point frequency signal segment echo signal obtained after deskewing reception is expressed as:
(5) (5)
其中,为去斜接收的点频信号的时域形式,/>为去斜接收的点频信号幅度参数;in, is the time domain form of the deskewed received point frequency signal,/> It is the amplitude parameter of the deskewed received point frequency signal;
对去斜接收的点频信号段回波信号做二维傅里叶变换,获得去斜接收后点频信号的二维频谱信号理论模型,其表达式为:Perform a two-dimensional Fourier transform on the echo signal of the deskewed received point frequency signal segment to obtain a two-dimensional spectrum signal theoretical model of the deskewed received point frequency signal. , its expression is:
(6) (6)
其中,为距离频谱包络,/>为去斜接收后点频信号二维频谱的幅度参数。in, is the distance spectrum envelope,/> is the amplitude parameter of the two-dimensional spectrum of the point frequency signal after deskewing reception.
由公式(6)可知,对实际接收的点频信号段回波信号进行二维傅里叶变换,通过实际接收的点频信号段回波信号的二维频谱与去斜接收后得到的点频信号段回波信号的二维频谱信号理论模型的偏差,解算出超材料标签调制导致的方位向频谱偏移量。具体地,所述步骤4包括:It can be seen from formula (6) that the two-dimensional Fourier transform is performed on the echo signal of the actually received point frequency signal segment. The two-dimensional spectrum of the echo signal of the actually received point frequency signal segment and the point frequency obtained after deslope reception are Deviation of the two-dimensional spectrum signal theoretical model of the signal segment echo signal, and the azimuth spectrum offset caused by metamaterial tag modulation is solved . Specifically, step 4 includes:
基于已知的和解算得到的/>,对线性调频段进行补偿成像。based on what is known and solved/> , perform compensation imaging on the linear frequency modulation band.
首先构建距离向调制补偿因子,对公式(2)的线性调频段回波信号进行距离调制补偿,得到距离向补偿后的回波信号/>,其表达式为:First, construct the distance modulation compensation factor , perform distance modulation compensation on the linear frequency modulation band echo signal of formula (2), and obtain the range-compensated echo signal/> , its expression is:
(7) (7)
对公式(7)补偿残余视频相位因子,得到去除残余视频相位后的回波信号/>,其表达式为:Compensating the residual video phase factor for equation (7) , get the echo signal after removing the residual video phase/> , its expression is:
(8) (8)
其中,表示幅度因子,由于在聚焦成像中信号相位的处理最为关键,因此在后续成像算法分析中不再考虑该项。in, Represents the amplitude factor. Since signal phase processing is the most critical in focused imaging, this item will no longer be considered in subsequent imaging algorithm analysis.
采用波束域算法进行聚焦成像,对公式(8)进行时频变换,得到时频变换后的信号,其表达式为:The beam domain algorithm is used for focused imaging, and formula (8) is subjected to time-frequency transformation to obtain the time-frequency transformed signal. , its expression is:
(9) (9)
其中表示距离向波数参量,/>为方位向位置参量。in Represents the distance wavenumber parameter,/> is the azimuth position parameter.
对公式(9)在方位向进行傅里叶变换,获取二维空间频域信号,其表达式为:Perform Fourier transform on equation (9) in the azimuth direction to obtain the two-dimensional spatial frequency domain signal. , its expression is:
(10) (10)
其中,表示方位波数,该式中利用/>对方位波数进行了补偿,消除方位向的调制影响。in, Represents the azimuthal wave number, which is used in this formula/> The azimuth wavenumber is compensated to eliminate the modulation effect in the azimuth direction.
公式(10)中最后一项为脉内运动引入的多普勒频移,可通过复乘补偿因子/>抵消,得到补偿多普勒频移后的二维频域信号/>:The last term in formula (10) The Doppler shift introduced for intrapulse motion can be multiplied by a compensation factor/> Cancel, and obtain the two-dimensional frequency domain signal after the Doppler frequency shift is compensated/> :
(11) (11)
构建参考函数,其中参考距离/>取雷达到景中心的距离,参考函数与公式(12)相乘得到在参考距离处聚焦的二维频域信号:Build reference function , where the reference distance/> Taking the distance from the radar to the center of the scene, the reference function is multiplied by formula (12) to obtain the two-dimensional frequency domain signal focused at the reference distance. :
(12) (12)
令对公式(12)进行Stolt插值变换,得到变换后信号/>,其表达式为:make Perform Stolt interpolation transformation on formula (12) to obtain the transformed signal/> , its expression is:
(13) (13)
其中,表示Stolt插值变换以后的纵坐标。in, Represents the ordinate after Stolt interpolation transformation.
最后对上式进行二维快速傅里叶逆变换(2D-IFFT),即可得到完全聚焦的二维图像。Finally, a two-dimensional inverse fast Fourier transform (2D-IFFT) is performed on the above equation to obtain a fully focused two-dimensional image.
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。It is easy for those skilled in the art to understand that the above descriptions are only preferred embodiments of the present invention and are not intended to limit the present invention. Any modifications, equivalent substitutions and improvements, etc., made within the spirit and principles of the present invention, All should be included in the protection scope of the present invention.
Claims (5)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311175498.0A CN116930969B (en) | 2023-09-13 | 2023-09-13 | Metamaterial tag target positioning and imaging method based on synthetic aperture radar |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311175498.0A CN116930969B (en) | 2023-09-13 | 2023-09-13 | Metamaterial tag target positioning and imaging method based on synthetic aperture radar |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116930969A true CN116930969A (en) | 2023-10-24 |
CN116930969B CN116930969B (en) | 2023-12-26 |
Family
ID=88394361
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311175498.0A Active CN116930969B (en) | 2023-09-13 | 2023-09-13 | Metamaterial tag target positioning and imaging method based on synthetic aperture radar |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116930969B (en) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101111852A (en) * | 2004-06-22 | 2008-01-23 | 伍比克公司 | Rfid system utilizing parametric reflective technology |
JP2009236545A (en) * | 2008-03-26 | 2009-10-15 | Mitsubishi Space Software Kk | Image processing apparatus, image processing system, image processing method, and image processing program |
US20120256730A1 (en) * | 2011-04-11 | 2012-10-11 | Lone Star Ip Holdings, Lp | Interrogator and System Employing the Same |
CN110187333A (en) * | 2019-05-23 | 2019-08-30 | 天津大学 | A RFID tag positioning method based on synthetic aperture radar technology |
CN110187344A (en) * | 2019-05-31 | 2019-08-30 | 天津大学 | Recognition and positioning method of marked items on conveyor belt based on inverse synthetic aperture radar |
CN110850401A (en) * | 2019-08-27 | 2020-02-28 | 天津大学 | RFID label positioning method based on motion model and synthetic aperture |
-
2023
- 2023-09-13 CN CN202311175498.0A patent/CN116930969B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101111852A (en) * | 2004-06-22 | 2008-01-23 | 伍比克公司 | Rfid system utilizing parametric reflective technology |
JP2009236545A (en) * | 2008-03-26 | 2009-10-15 | Mitsubishi Space Software Kk | Image processing apparatus, image processing system, image processing method, and image processing program |
US20120256730A1 (en) * | 2011-04-11 | 2012-10-11 | Lone Star Ip Holdings, Lp | Interrogator and System Employing the Same |
CN110187333A (en) * | 2019-05-23 | 2019-08-30 | 天津大学 | A RFID tag positioning method based on synthetic aperture radar technology |
CN110187344A (en) * | 2019-05-31 | 2019-08-30 | 天津大学 | Recognition and positioning method of marked items on conveyor belt based on inverse synthetic aperture radar |
CN110850401A (en) * | 2019-08-27 | 2020-02-28 | 天津大学 | RFID label positioning method based on motion model and synthetic aperture |
Also Published As
Publication number | Publication date |
---|---|
CN116930969B (en) | 2023-12-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106842210B (en) | A kind of new multiple submatrixes synthetic aperture sonar fast imaging algorithm | |
CN107229048B (en) | A high-resolution wide-range SAR moving target velocity estimation and imaging method | |
CN103091674B9 (en) | High-resolution imaging method of space target based on HRRP sequence | |
CN103760558B (en) | Terahertz radar ISAR imaging method | |
CN103760548B (en) | A kind of signal processing method based on laser radar waveform of coherent system | |
CN103454632B (en) | Imaging method for one-station fixed FM continuous wave bistatic SAR | |
CN103969628B (en) | A kind of synthetic aperture radar PFA imaging algorithm based on compressed sensing | |
CN103454638B (en) | Circular synthetic aperture radar three-dimension layer tomographic imaging method | |
CN102176011A (en) | Method for realizing three-dimensional coherent imaging by ground penetrating radar under near field condition | |
CN107607952B (en) | Three-dimensional synthetic aperture radar imaging method based on electromagnetic vortex wave | |
JP6179940B2 (en) | Doppler imaging signal transmitter, Doppler imaging signal receiver, Doppler imaging system and method | |
CN107102327A (en) | SAR imaging methods based on LFM PC multiplex modulated signals and polar format algorithm | |
CN109444901A (en) | A kind of multiple submatrixes SAS sub-aperture image method under isomerous environment | |
CN102121990B (en) | Estimation Method of Inverse Synthetic Aperture Radar's Target Rotational Speed Based on Space-Time Analysis | |
CN109164421A (en) | A kind of object detection method based on Two-Dimensional Reconstruction algorithm | |
CN109358330B (en) | A Linear Time Synchronization Error Estimation Method for Bistatic Forward Looking SAR | |
CN109597076B (en) | Data processing method and device for ground-based synthetic aperture radar | |
CN103901425A (en) | Radar and method for measuring transverse velocity of moving object | |
CN108008387B (en) | An airborne array antenna down-view three-dimensional imaging method | |
CN110879391B (en) | Radar image dataset creation method based on electromagnetic simulation and missile echo simulation | |
CN104237886A (en) | High-precision synthetic aperture radar imaging method | |
CN101793960A (en) | High-accuracy dynamic threshold target monitoring method with chirping characteristics | |
AU2018313118A1 (en) | Sensor array imaging device | |
CN116930969A (en) | Metamaterial tag target positioning and imaging method based on synthetic aperture radar | |
CN109738893A (en) | Method and device for generating echo data of dual-base synthetic aperture radar |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |