CN103135098B - High-accuracy interference specific absorption rate (SAR) system performance analysis method - Google Patents
High-accuracy interference specific absorption rate (SAR) system performance analysis method Download PDFInfo
- Publication number
- CN103135098B CN103135098B CN201310048793.XA CN201310048793A CN103135098B CN 103135098 B CN103135098 B CN 103135098B CN 201310048793 A CN201310048793 A CN 201310048793A CN 103135098 B CN103135098 B CN 103135098B
- Authority
- CN
- China
- Prior art keywords
- sar
- interference sar
- interference
- noise
- ratio
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000004458 analytical method Methods 0.000 title abstract description 3
- 238000010521 absorption reaction Methods 0.000 title abstract 2
- 238000000034 method Methods 0.000 claims abstract description 17
- 238000003384 imaging method Methods 0.000 claims abstract description 15
- 238000001914 filtration Methods 0.000 claims abstract description 6
- 238000004088 simulation Methods 0.000 claims abstract description 6
- 238000013461 design Methods 0.000 claims description 16
- 238000012545 processing Methods 0.000 claims description 16
- 241001310793 Podium Species 0.000 claims description 4
- 230000001427 coherent effect Effects 0.000 claims description 4
- 230000006835 compression Effects 0.000 claims description 4
- 238000007906 compression Methods 0.000 claims description 4
- 230000010354 integration Effects 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 238000012795 verification Methods 0.000 claims description 3
- 230000000007 visual effect Effects 0.000 claims description 3
- 238000004804 winding Methods 0.000 claims description 3
- 238000011156 evaluation Methods 0.000 claims description 2
- VMXUWOKSQNHOCA-UKTHLTGXSA-N ranitidine Chemical compound [O-][N+](=O)\C=C(/NC)NCCSCC1=CC=C(CN(C)C)O1 VMXUWOKSQNHOCA-UKTHLTGXSA-N 0.000 claims description 2
- 230000008569 process Effects 0.000 abstract description 7
- 238000013459 approach Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000007812 deficiency Effects 0.000 description 2
- 238000002592 echocardiography Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000003786 synthesis reaction Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000013501 data transformation Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses a high-accuracy interference specific absorption rate (SAR) system performance analysis method. Interference system performance is overall verified through simulation and treating process of a whole closed loop. The method comprises: step 1, interference SAR system parameters are calculated by adopting conventional theoretical formulas according to indicator requirements of users; step 2, inference SAR echo signal imitation is conducted according to the obtained interference SAR system parameters in step 1; step 3, imaging process and kinematic error compensation are conducted on the emulated echo data in step 2 so that relevant duplication images are obtained; step 4, interference process is conducted on the obtained duplication images in step 3 so as to achieve the purpose that the same feature corresponds to the same pixel on the images, and phase noise is filtered by phase filtering; and step 5, system performance indexes are analyzed according to the imaging results obtained in the step 3 and height inversion results obtained in the step 4, and whether the system can reach the designed requirements or not is verified.
Description
Technical field
The invention belongs to signal processing technology field, be specifically related to a kind of high precision interference SAR system method for analyzing performance.
Background technology
Along with modern radar system is day by day complicated, function is constantly perfect, simulation of Radar System has caused increasing attention.The today particularly developing rapidly in capacity and the speed of computing machine, for meticulous simulation of Radar System provides material base.By Computer Simulation, verify rationality, the realizability of radar system overall design philosophy.SAR needs to carry out population parameter design and demonstration according to technical indicator before development and design, and current traditional radar population parameter design is based on theoretical formula model, and this technical method exists deficiency below:
(1) model is inaccurate: theoretical formula model, obtain, yet actual system behavior index is subject to the impact of multiple error, can not represent by desirable analytic formula through approximate processing.
(2) do not consider various systematic error impacts: in interference SAR system, exist the various error sources such as channel noise, thermonoise, kinematic error, the existence of these errors finally affects interference performance.
(3) do not carry out overall synthetic design: current carrying out population parameter when design, for some or several parameters, design separately respectively, do not carry out overall synthetic design, the error that whole interference SAR echo data is obtained, produced in imaging processing and interference handle link process does not have system demonstration.
Summary of the invention
The present invention overcomes the defect of prior art, for overcoming the deficiency of current interference SAR system population parameter design and demonstration, propose a kind of high precision interference SAR system performance synthesis analytical approach, by emulation and the processing procedure of whole closed loop, totally verify interference system performance.
A method for high precision interference SAR system performance evaluation, comprises the steps:
Step 1, according to user's index demand, adopt conventional theoretical formula to calculate interference SAR system parameter, parameter comprises podium level, beam center visual angle, beam angle, pulse repetition rate, data sampling frequency, antenna parameter, transmitter power, noise figure, system loss, fuzzy ratio, quantizes signal to noise ratio (S/N ratio), thereby obtains interference SAR initial parameter value;
Step 2, the interference SAR system parameter obtaining according to step 1, carry out the simulation of interference SAR echoed signal; Set up interference SAR hyperchannel original echoed signals realistic model, and error source is incorporated in original echo data; Interference SAR echoed signal is based upon on model space geometric basis, and according to SAR data acquisition principle, in orientation, each PRF place obtains original echo data;
Step 3, the echo data of emulation in step 2 is carried out to the compensation of imaging processing and kinematic error, obtain coherence's complex pattern; Imaging processing is carried out azimuth-range compression to original echo data, makes focusing on a synthetic aperture self-energy of each pixel;
Step 4, the complex pattern obtaining in step 3 is interfered to processing, realize same atural object corresponding same pixel on image, phase filtering filtering phase noise, phase unwrapping reverts to original phase by the phase place of winding, according to interference SAR ultimate principle, by launching phase place, resolve landform value, and carry out inverting elevation by the original phase after launching;
Step 5, according to the elevation inversion result obtaining in the imaging results of obtaining in step 3 and step 4, system performance index is analyzed, whether verification system can reach the demand of design.
Since then, completed a kind of high precision interference SAR system performance synthesis analytical approach.
Error source described in step 2 comprises, channel noise, hyperchannel imbalance, thermonoise, coherent speckle noise, fuzzy noise, quantizing noise;
In step 5, design requirement index comprises: vertical accuracy, resolution, mapping bandwidth, operating distance, signal to noise ratio (S/N ratio), peak sidelobe ratio, integration secondary lobe ratio, fuzzy ratio.
The present invention compared with prior art advantage is: the present invention has introduced interference SAR real system error source, and comparing with current method is a kind of analytical approach that more approaches real system; From original echo, getting elevation inverting, the whole active link of interference SAR system is simulated, is a kind of interference SAR system analytical approach.
Accompanying drawing explanation
Fig. 1 interference SAR is comprehensively analyzed process flow diagram;
Fig. 2 interference SAR model space geometric;
Fig. 3 correlativity binary channels analogue echoes process flow diagram schematic diagram;
Fig. 4 imaging processing and kinematic error compensation schematic flow sheet;
Fig. 5 interferes treatment scheme schematic diagram.
Embodiment
As shown in Figure 1, a kind of high precision interference SAR system method for analyzing performance, comprises step below:
Step 1, according to user's index demand, adopt conventional theoretical formula to calculate interference SAR system parameter.The parameter of design comprises podium level, beam center visual angle, beam angle, pulse repetition rate, data sampling frequency, antenna parameter, transmitter power, noise figure, system loss, fuzzy ratio, quantize signal to noise ratio (S/N ratio), thereby obtain interference SAR initial parameter value.
Step 2, the interference SAR parameter obtaining according to step 1, carry out the simulation of interference SAR echoed signal.Set up interference SAR hyperchannel original echoed signals realistic model, and error is incorporated in original echo data.Error source comprises, channel noise, hyperchannel imbalance, thermonoise, coherent speckle noise, fuzzy noise, quantizing noise.Interference SAR echoed signal is based upon on model space geometric basis, and according to SAR data acquisition principle, in orientation, each PRF place obtains original echo data.The flow process of binary channels interference SAR analogue echoes as shown in Figure 3, model model space geometric, interference SAR model space geometric as shown in Figure 2.X is platform speed direction, y be distance to, z is vertical ground direction, H is podium level, h is impact point height value, A1, A2 is two slave antenna positions, B is base length, α is baseline inclination angle.Radar every one pulse-recurrence time emission wavelength be λ, pulsewidth is T
r, the linear FM signal that frequency modulation rate is K:
ω wherein
0be carrier frequency, rect () is rectangular window function, and after coherent demodulation, radar any point (t, τ) in t receives its irradiation area is constantly located base band echo 2D signal and can be expressed as:
T wherein, τ is respectively the fast time of slow time of orientation and distance, and λ is the wavelength that transmits, and c is the light velocity, and σ is target backscattering coefficient, w
a(t) w
r(τ) be respectively orientation to distance to antenna radiation pattern, R is the instantaneous oblique distance of point target and radar.
Step 3, artificial echo data are carried out to the compensation of imaging processing and kinematic error, obtain coherence's complex pattern.Imaging processing is carried out azimuth-range compression to original echo data, makes focusing on a synthetic aperture self-energy of each pixel.
As shown in Figure 4, motion compensation is embedded in imaging processing, comprises step:
1), for airborne interference SAR system, before data are processed, can first utilize phase function to carry out elementary motion compensation.
2) by orientation FFT, data are transformed to range-Dopler domain from signal space.
3) be multiplied by first-phase bit function and carry out linear perturbation for a certain reference distance.
4) by data being transformed to two-dimentional frequency field (being wave number space) from range-Dopler domain apart from FFT
5) being multiplied by second-phase bit function completes band secondary range and compresses and complete range migration correction.
6) by data transformation being returned to range-Dopler domain apart from IFFT.
7) be multiplied by the slow variable-azimuth residual phase that the compensation of third phase bit function is introduced by linear perturbation function.
8) for airborne interference SAR system, by orientation IFFT, data are converted back apart from bearing signal territory from range-Dopler domain, utilize phase function to complete collateral motion compensation, then by orientation FFT layback Doppler domain.
9) be multiplied by the 4th phase function and complete Azimuth Compression, finally by orientation IFFT, obtain SAR complex pattern.
Step 4, complex pattern is interfered to processing.As shown in Figure 5, binary channels complex pattern is interfered and processed acquisition elevation information.First Channel Image is carried out to sub-pixel accuracy registration, make the same plot of a pair of pixel corresponding ground in two sub-pictures.Then interferometric phase is carried out to filtering, owing to there is the noises such as thermonoise, registration error, quantization error in the actual interferometric phase obtaining, these noises bring error to the measurement of higher degree, also to interferometric phase solution tape wrap, come difficult.Then carry out phase unwrapping processing, what obtain due to conjugate multiplication after two sub-picture registrations is interferometric phase winding value, need to carry out two-dimensional phase unwrapping and recover true phase place corresponding to original wave path-difference.
Step 5, system performance analysis and demonstration.According to imaging results and elevation inversion result, system performance index is analyzed, whether verification system can reach the demand of design.These indexs comprise: vertical accuracy, resolution, mapping bandwidth, operating distance, signal to noise ratio (S/N ratio), peak sidelobe ratio, integration secondary lobe ratio, fuzzy ratio.
In sum, these are only preferred embodiment of the present invention, be not intended to limit protection scope of the present invention.Within the spirit and principles in the present invention all, any modification of doing, be equal to replacement, improvement etc., within all should being included in protection scope of the present invention.
Claims (1)
1. a method for high precision interference SAR system performance evaluation, is characterized in that, comprises the steps:
Step 1, according to user's index demand, adopt conventional theoretical formula to calculate interference SAR system parameter, parameter comprises podium level, beam center visual angle, beam angle, pulse repetition rate, data sampling frequency, antenna parameter, transmitter power, noise figure, system loss, fuzzy ratio, quantize signal to noise ratio (S/N ratio), thereby obtain interference SAR initial parameter value;
Step 2, according to interference SAR system parameter, carry out the simulation of interference SAR echoed signal; Set up interference SAR hyperchannel original echoed signals realistic model, and error source is incorporated in original echo data;
Described error source comprises, channel noise, hyperchannel imbalance, thermonoise, coherent speckle noise, fuzzy noise, quantizing noise;
Interference SAR echoed signal is based upon on model space geometric basis, and according to SAR data acquisition principle, in orientation, each PRF place obtains original echo data;
Step 3, the echo data of emulation in step 2 is carried out to the compensation of imaging processing and kinematic error, obtain coherence's complex pattern; Imaging processing is carried out azimuth-range compression to original echo data, makes focusing on a synthetic aperture self-energy of each pixel;
Judge the imaging results requirement that whether touches the mark, if fall short of the target requirement, revise the design parameter of interference SAR and return to step 2; If touched the mark requirement, entered step 4;
Step 4, the complex pattern obtaining in step 3 is interfered to processing, realize same atural object corresponding same pixel on image, phase filtering filtering phase noise, phase unwrapping reverts to original phase by the phase place of winding, according to interference SAR ultimate principle, by launching phase place, resolve landform value, and carry out inverting elevation by the original phase after launching;
Step 5, according to the elevation inversion result obtaining in the imaging results of obtaining in step 3 and step 4, system performance index is analyzed, whether verification system can reach the demand of design;
Whether the precision that judges elevation inversion result reaches measuring accuracy, if do not reach measuring accuracy, revises interference SAR system parameter and returns to step 2; If reached measuring accuracy, completed a kind of high precision interference SAR system method for analyzing performance;
Wherein design requirement index comprises: vertical accuracy, resolution, mapping bandwidth, operating distance, signal to noise ratio (S/N ratio), peak sidelobe ratio, integration secondary lobe ratio, fuzzy ratio.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310048793.XA CN103135098B (en) | 2013-02-06 | 2013-02-06 | High-accuracy interference specific absorption rate (SAR) system performance analysis method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310048793.XA CN103135098B (en) | 2013-02-06 | 2013-02-06 | High-accuracy interference specific absorption rate (SAR) system performance analysis method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103135098A CN103135098A (en) | 2013-06-05 |
| CN103135098B true CN103135098B (en) | 2014-03-26 |
Family
ID=48495180
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310048793.XA Active CN103135098B (en) | 2013-02-06 | 2013-02-06 | High-accuracy interference specific absorption rate (SAR) system performance analysis method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103135098B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104535975A (en) * | 2014-12-17 | 2015-04-22 | 北京无线电计量测试研究所 | Image calibration device and method for synthetic aperture radar echo simulator |
| CN106123894B (en) * | 2016-07-25 | 2019-07-23 | 中国科学院电子学研究所 | Based on the matched InSAR/INS Combinated navigation method of interference fringe |
| CN109490881A (en) * | 2018-11-21 | 2019-03-19 | 中国科学院电子学研究所 | Interference SAR measurement of higher degree system and measurement method based on vortex electromagnetic wave |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102323571A (en) * | 2011-05-31 | 2012-01-18 | 北京航空航天大学 | Distribution method of satellite-borne dual-antenna SAR (Synthetic Aperture Radar) interferometric calibrator with comprehensive overall parameter |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6911931B2 (en) * | 2002-10-24 | 2005-06-28 | The Regents Of The University Of California | Using dynamic interferometric synthetic aperature radar (InSAR) to image fast-moving surface waves |
-
2013
- 2013-02-06 CN CN201310048793.XA patent/CN103135098B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102323571A (en) * | 2011-05-31 | 2012-01-18 | 北京航空航天大学 | Distribution method of satellite-borne dual-antenna SAR (Synthetic Aperture Radar) interferometric calibrator with comprehensive overall parameter |
Non-Patent Citations (4)
| Title |
|---|
| 分布式卫星InSAR系统性能分析与仿真技术研究;孟智勇;《中国优秀硕士学位论文全文数据库》;20120229(第2期);第34页第3段-第39页第3段、第42页第8-13段，第48页第9段、第73页第11段，图3.4-3.12，图5.3 * |
| 孟智勇.分布式卫星InSAR系统性能分析与仿真技术研究.《中国优秀硕士学位论文全文数据库》.2012,(第2期), |
| 干涉合成孔径雷达成像技术;张跃,万昳;《电子科技》;20120331;第25卷(第3期);107-109 * |
| 张跃,万昳.干涉合成孔径雷达成像技术.《电子科技》.2012,第25卷(第3期), |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103135098A (en) | 2013-06-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Franceschetti et al. | SAR sensor trajectory deviations: Fourier domain formulation and extended scene simulation of raw signal | |
| Pu et al. | Joint sparsity-based imaging and motion error estimation for BFSAR | |
| CN108872985B (en) | Near-field circumference SAR rapid three-dimensional imaging method | |
| Gilman et al. | Transionospheric synthetic aperture imaging | |
| Cimmino et al. | Efficient spotlight SAR raw signal simulation of extended scenes | |
| CN102749621B (en) | Bistatic synthetic aperture radar (BSAR) frequency domain imaging method | |
| Pu et al. | Fast factorized backprojection imaging algorithm integrated with motion trajectory estimation for bistatic forward-looking SAR | |
| CN103176170A (en) | SAR (synthetic aperture radar) echo simulating method based on GPU (graphics processing unit) parallel computing | |
| CN105677942A (en) | Rapid simulation method of repeat-pass spaceborne natural scene SAR complex image data | |
| CN103869311A (en) | Real beam scanning radar super-resolution imaging method | |
| CN103018740B (en) | InSAR (Interferometric Synthetic Aperture Radar) imaging method based on curved surface projection | |
| CN105467369B (en) | A kind of target echo simulation method and apparatus | |
| Magnard et al. | Processing of MEMPHIS Ka-band multibaseline interferometric SAR data: From raw data to digital surface models | |
| CN101685159A (en) | Method for constructing spaceborne SAR signal high precision phase-keeping imaging processing platform | |
| CN103323840A (en) | Method for time alignment between interference SAR echo data and platform motion and gesture data | |
| CN103135098B (en) | High-accuracy interference specific absorption rate (SAR) system performance analysis method | |
| CN107607945A (en) | A kind of scanning radar forword-looking imaging method based on spatial embedding mapping | |
| Qiu et al. | A Bistatic SAR Raw Data Simulator Based on Inverse $\omega {-} k $ Algorithm | |
| Ng et al. | Estimation of the total rotational velocity of a non‐cooperative target with a high cross‐range resolution three‐dimensional interferometric inverse synthetic aperture radar system | |
| Kou et al. | Circular SAR processing using an improved omega-k type algorithm | |
| CN102890270B (en) | Echo simulation method of bi-static synthetic aperture radar of fixed station | |
| CN102890277B (en) | Range migration imaging method of shift invariant bi-static synthetic aperture radar | |
| Zhu et al. | Raw signal simulation of synthetic aperture radar altimeter over complex terrain surfaces | |
| Grebner et al. | Simultaneous localization and mapping (SLAM) for synthetic aperture radar (SAR) processing in the field of autonomous driving | |
| Kim et al. | SAR raw data simulation for multiple‐input multiple‐output video synthetic aperture radar using beat frequency division frequency modulated continuous wave |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20201231 Address after: 250101 No.9, Kuangyuan Road, Gongye North Road, Wangsheren street, Licheng District, Jinan City, Shandong Province Patentee after: Jigang Defense Technology Co.,Ltd. Address before: 100190 No. 19 West North Fourth Ring Road, Haidian District, Beijing Patentee before: Aerospace Information Research Institute,Chinese Academy of Sciences Effective date of registration: 20201231 Address after: 100190 No. 19 West North Fourth Ring Road, Haidian District, Beijing Patentee after: Aerospace Information Research Institute,Chinese Academy of Sciences Address before: 100180 No.19, Beisihuan West Road, Haidian District, Beijing Patentee before: Institute of Electronics, Chinese Academy of Sciences |
|
| TR01 | Transfer of patent right |