CN108120980A - A kind of implementation method of the FPGA of satellite-borne SAR multi-modal imaging signal processing algorithm - Google Patents

A kind of implementation method of the FPGA of satellite-borne SAR multi-modal imaging signal processing algorithm Download PDF

Info

Publication number
CN108120980A
CN108120980A CN201711325281.8A CN201711325281A CN108120980A CN 108120980 A CN108120980 A CN 108120980A CN 201711325281 A CN201711325281 A CN 201711325281A CN 108120980 A CN108120980 A CN 108120980A
Authority
CN
China
Prior art keywords
data
satellite
sub
borne sar
fpga
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201711325281.8A
Other languages
Chinese (zh)
Inventor
朱岱寅
崔爱欣
金微微
贺雪莉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN201711325281.8A priority Critical patent/CN108120980A/en
Publication of CN108120980A publication Critical patent/CN108120980A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/03Details of HF subsystems specially adapted therefor, e.g. common to transmitter and receiver
    • G01S7/032Constructional details for solid-state radar subsystems
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/9004SAR image acquisition techniques
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The present invention discloses a kind of implementation method of the FPGA of satellite-borne SAR multi-modal imaging signal processing algorithm, includes the following steps:Step 1, satellite-borne SAR raw radar data and parameter profit are sent to fpga chip by host computer and gigabit Ethernet;Step 2, echo data is subjected to sub-aperture division, and calculates required parameter in sub-aperture PFA processing procedures;Step 3, Estimation of Doppler central frequency is carried out according to echo data;Step 4, sub- pore size data is carried out based on the distance of dimensional variation to processing;Step 5, orientation high-precision Sinc interpolation is carried out;Step 6, bidimensional FFT imagings are carried out and data are sent back into host computer by Ethernet;Step 7, radiant correction processing and weighting splicing are done to sub- subaperture image in host computer, is shown after doing histogram equalization in host computer.Such method can realize the unitized imaging to the wide data of satellite-borne SAR multi-mode wide swath, have the characteristics that high-precision, it is efficient.

Description

A kind of implementation method of the FPGA of satellite-borne SAR multi-modal imaging signal processing algorithm
Technical field
The invention belongs to radar imaging technology field and digital signal processing technique field, more particularly to a kind of satellite-borne SAR The implementation method of the FPGA of multi-modal imaging signal processing algorithm realizes that the satellite-borne SAR based on sub-aperture stitching is more using FPGA Mode imaging processing.
Background technology
Satellite-borne synthetic aperture radar (SAR) is a kind of high-resolution imaging radar, it can be worked with round-the-clock round-the-clock, tool Have the ability for penetrating soil and vegetation, have been widely used for battle reconnaissance, Hitting Effect Evaluation, target identification, resource exploration, The key areas of many national defence such as environmental protection, the condition of a disaster detection and national economy.High resolution wide swath is satellite-borne SAR development One of important directions.The remote sensing such as one side earthquake and disaster scenarios it assessment, target detection identification, military attack recruitment evaluation should Very high requirement is proposed with the spatial resolution index to SAR image;Exploration that another aspect targets in ocean monitors, disaster area is met an urgent need The ability that remote sensing applications requirement radar system is waited to possess wide cut imaging.The satellite-borne SAR of big mapping scene domain and more imaging patterns Real Time Image System is increasingly becoming current research hotspot.
SAR Real-time processings system can carry out accident under conditions of round-the-clock, round-the-clock real-time on star Wide area is observed and quick acquisition of information, so as to make correct decisions within the shorter time.Therefore research satellite-borne SAR is real-time Imaging processing techniques have realistic meaning.Current Space-borne SAR Imaging treatment technology, typical imaging algorithm have range Doppler (RD) algorithm, Chirp Scaling (CS) algorithms and wave-number domain (Omega-K) algorithm etc..And based on the sub-aperture stitching of PFA into As algorithm, imaging precision is high, and processing speed is fast, and Irnaging procedures can be to band, slip a variety of satellite-borne SAR patterns such as pack and TOPS Echo data carries out unitized processing, without being changed over according to operating mode as processing framework.
Field programmable gate array (Filed Programmable Gate Array, FPGA) has abundant look-up table And register resources, using parallel processing manner, arithmetic speed is fast, low in energy consumption.SAR echo datas have that data volume is big, data The characteristics of rate is high compared to DSP, using FPGA processing, can be rapidly completed FFT, complex multiplication addition, cordic algorithm Etc. Digital Signal Processing basic operation.
The content of the invention
To realize high-precision, the processing of efficient satellite-borne SAR multi-modal imaging, the present invention provides a kind of satellite-borne SAR multi-mode The implementation method of the FPGA of imaging signal Processing Algorithm is realized and unitized to the wide data of satellite-borne SAR multi-mode wide swath Imaging.
In order to achieve the above objectives, solution of the invention is:
A kind of implementation method of the FPGA of satellite-borne SAR multi-modal imaging signal processing algorithm, includes the following steps:
Step 1, satellite-borne SAR raw radar data and parameter profit are sent to FPGA cores by host computer and gigabit Ethernet Piece;
Step 2, echo data is subjected to sub-aperture division, and calculates required parameter in sub-aperture PFA processing procedures;
Step 3, Estimation of Doppler central frequency is carried out according to echo data;
Step 4, sub- pore size data is carried out based on the distance of dimensional variation to processing;
Step 5, orientation high-precision Sinc interpolation is carried out;
Step 6, bidimensional FFT imagings are carried out and data are sent back into host computer by Ethernet;
Step 7, radiant correction processing and weighting splicing are done to sub- subaperture image in host computer, do after histogram equalization Host computer is shown.
In above-mentioned steps 2, when echo data is carried out sub-aperture division, if the satellite-borne SAR data of big bandwidth are divided into It is dry to have the sub-aperture of lap and complete imaging.
The detailed content of above-mentioned steps 3 is:It realizes that dislocation is multiplied using shift register first and obtains auto-correlation function, then Data deposit FIFO after multiplication is waited into addition corresponding with next line result, the data line finally obtained is tired out using floating-point Device is added to carry out full flowing water floating-point to add up, the horizontal echo-signal for obtaining multiple range cells is averaged, then does by CORDIC cores ATAN operates to obtain accurate doppler centroid.
The detailed content of above-mentioned steps 4 is:Satellite-borne SAR distance is read one by one first from DDR3SDRAM to carry out to pulse Motion compensation, then using PCS principles multiplied again with FFT computings, often handled waited in the write-in DDR of a pulse transposition into The processing of row orientation.
The detailed content of above-mentioned steps 5 is:By spaceborne data, Pulse by Pulse transposition is read from DDR first, double by SAR64 Accuracy floating-point echo data turns fixed point IP kernel by floating-point and is converted to fixed-point data, then the integer part for pinpointing echo data is made For address ram, piecemeal storage echo data and fractional part, the weighted sum that echo data is completed by addressing operation obtain SINC interpolation results.
The detailed content of above-mentioned steps 7 is:Orientation antenna direction is multiplied by orientation data field first in sub-aperture The inverse of figure eliminates influence of the antenna radiation pattern to the point target echo amplitude of different azimuth position in sub-aperture, using weighting Stitching algorithm be averaged to realize seamlessly transitting for overlapping region, obtains full aperture image, obtained full aperture image is converted into Gray-scale map shows last imaging results figure in host computer after doing histogram equalization.
After using the above program, the present invention realizes Estimation of Doppler central frequency, sub-aperture using based on FPGA technology A series of imaging signal process flows such as PFA imaging algorithms, radiant correction, image mosaic.First, it is wide to wide swath spaceborne SAR echo datas carry out sub-aperture PFA imagings, and PFA algorithm performances influence from angle of squint, pass through orientation Keystone is converted, and can be solved the more resolution elements walk problem away from imaging area centers scatter point, can be met spaceborne SAR is comprehensive, high-resolution requirement;The scallop effect of sub-aperture image, and splicing energy are eliminated by radiant correction Enough obtain a wide range of scene imaging result of satellite-borne SAR.The Irnaging procedures can be a variety of spaceborne to band, slip pack and TOPS etc. SAR mode echoes data carry out unitized processing, without being changed over according to operating mode as processing framework.
Description of the drawings
Fig. 1 is the principle of the present invention Organization Chart;
Fig. 2 is Estimation of Doppler central frequency module frame chart;
Fig. 3 is sub-aperture PFA imaging flow charts;
Fig. 4 is orientation high-precision Sinc interpolating module block diagrams;
Fig. 5 is image mosaic display module figure;
Fig. 6 is MATLAB handling result figures;
Fig. 7 is FPGA handling result figures
Fig. 8 is the partial enlarged view of FPGA handling result figures;
Wherein, (a), (b), (c) are the partial enlarged view on airport in imaging results, harbour and mountain range respectively.
Specific embodiment
Below with reference to attached drawing, technical scheme and advantageous effect are described in detail.
As shown in Figure 1, it is the FPGA hardware realization system block diagram of the present invention.Space-borne SAR Imaging system is mainly by four cores Core module is formed:The data transmission transposition module (caching of the transmission of satellite-borne SAR echo data to fpga chip and data in itself And transposition), Estimation of Doppler central frequency module (accurate estimating Doppler centre frequency to improve picture quality), sub-aperture PFA algorithm processing modules (cascade based on PCS and SINC interpolation is realized is compressed into picture to the bidimensional of sub- pore size data) and image Tiled display module (multiple sub-aperture images do registration and the splicing of radiant correction and its image).
Data transmission transposition module is illustrated first.
PC ends and the IP address of fpga chip hardware end in host computer are set, satellite-borne SAR echo data and track are joined Number is sent to fpga chip by data transmission module.
Ethernet module make use of tri-state Ethernet (Tri-Mode Ethernet MAC) IP kernel that Xilinx companies provide And User Datagram Protocol (UDP) and receive buffering area (RX_FIFO) and send buffering area (TX_FIFO) come complete data receipts Send out function.
DDR modules have used Xilinx companies to provide DDR3 memory interfaces (MIG) solution.The solution provides Memory Controller Hub and physical layer design, user, which is based on the program, directly to be operated DDR memories, and orbit parameter is deposited It stores up in parameter register, in satellite-borne SAR echo data deposit DDR3SDRAM.
Satellite-borne SAR data volume is big, and algorithm processing module is related to the ranks operation of data, to lead to after processing in distance Cross the transposition that DDR realizes matrix.Using the direct-type fragmented storage method of high transposition efficiency, stored in same a line of SDRAM The line feed time for reading data is greatly reduced in the row data and column data of equal length while being slightly increased and writing data line feed number Number.
Next Estimation of Doppler central frequency module is illustrated.
According to wiener-Xin Xin theorems, the power spectrum function of random process is the Fourier transformation of its auto-correlation function.
Doppler center f is contained in the phase of auto-correlation functiondcInformation, therefore by estimating rτPhase, so that it may To obtain fdc, it is shown below:
In specific implementation, horizontal echo is sampled with pulse recurrence frequency PRF, sampling interval T=1/PRF, related Function is r (KT), is had at this time
Radar return data is are done shift LD conjugate multiplication by specific FPGA realization methods first, by the knot after multiplication Fruit deposit FIFO cachings wait addition corresponding with next line result, and obtained last column data are carried out full flowing water floating number tires out Add, constantly enter data into a ends of floating-point adder, b ends are the feedback at adder result end, when the delay of adder is 12 In the clock cycle, b ends result is set to 0 when not feeding back initially, acquired results are input in ring register simultaneously, left final 12 as a result, gained is last accumulation result.
Finally by CORDIC (Coordinate Rotation Digital Computer Coordinate Rotation Digitals calculating sides Method) IP (Intelligence Property intellectual properties) core seeks arg { r (τ) }, f is obtained further according to formuladc
Table 1 estimates the comparison of test result and theoretical value for Doppler center.
Table 1
Next sub-aperture PFA algorithm processing modules are illustrated.
The distance of sub-aperture image processing is multiplied by QP function first to the PFA algorithms based on change of scale are used:
In formula, c is the light velocity, and τ is the fast time, and k is chirp rate, δrIt is distance to the change of scale factor, RaFor antenna phase Center is to the instantaneous distance of scene center.
Carry out Fourier transformation is done, multiplied by with filter function:
It is multiplied by QP function:
Wherein, fcFor carrier frequency.
Fourier transformation is carried out, is multiplied by motion compensation factor t=fc/(fc+fτ)·t:
Realized in FPGA it is above-mentioned it is a series of multiply again and Fourier transformation, and data transposition is stored in DDR, wherein When doing FFT processing, FFT is handled using an IP kernel using time-multiplexed mode three times, saves great amount of hardware resources.
Orientation processing converts the correction of further remaining linear range walk using Keystone, i.e., by high-precision Sinc interpolation operations realize conversion:
T=fc/(fc+fτ)·t' (7)
Wherein, t' is the time variable after orientation conversion.
It is converted by the Keystone of orientation, the more resolution elements away from imaging area centers scatter point can be solved and walked Dynamic problem, can meet that satellite-borne SAR is comprehensive, high-resolution requirement.
Fig. 4 is specific implementation block diagram of the orientation interpolation based on FPGA.Orientation processing is using based on piecemeal in FPGA The SINC interpolation methods of RAM groups, transposition reads distance to the floating data after processing from DDR, and floating data is converted to Fixed-point data pinpoints the integer part of echo data as address ram, piecemeal storage echo data and fractional part, by seeking The weighted sum that echo data is completed in location operation obtains SINC interpolation results.
Above-mentioned distance to after being disposed of orientation, realize the distance and bearing of two-dimensional signal to decoupling, then It is bidimensional FFT and obtains sub-aperture image result.
Fig. 5 is image mosaic display module figure, and correction function is multiplied by orientation data field to eliminate day in sub-aperture Influence of the line directional diagram to the point target echo amplitude of different azimuth position in sub-aperture, makes the point target of identical scattering strength exist Identical brightness is presented on image.
Seamlessly transitting for overlapping region is realized using weighted average stitching algorithm.
The validity further illustrated the present invention below by measured data, actual measurement radar parameter are as shown in table 2.
Table 2
Entirely Space-borne SAR Imaging system is realized on the monolithic Virtex7-XC7VX6907 chips of Xilinx companies, system Satellite-borne SAR data are done with distance to handle to the cascade sub-aperture of PCS and orientation Sinc interpolation, size of data 8192*8192, Sub-aperture size of data is 8192*512, can be compatible with anyon aperture quantity at this time and no more than the sub-aperture of 8192*512 sizes Footpath image procossing.
In the case where system frequency is 200MHz, 32 single-precision floating point sub-aperture numbers of 8192*512 pixels are handled According to, distance to processing time be 0.24 second, orientation is 0.08 second, and bidimensional FFT is 0.05 second, can be real within the 5.92s times Existing full aperture imaging.
Fig. 6 is MATLAB handling result figures, and Fig. 7 is FPGA handling result figures, it can be seen that two handling results substantially one It causes, but the processing time of FPGA is but greatly reduced.
Fig. 8 (a) (b) (c) is respectively the partial enlarged view on airport in imaging results, harbour and mountain range, it can be seen that this is Imaging resolution of uniting is high, and focusing effect is apparent and atural object is clearly demarcated, more can clearly show the important geographical distribution on island Information and the tendency on mountain range, contrast and noise are higher.
Based on system above, the present invention provides a kind of realization of the FPGA of satellite-borne SAR multi-modal imaging signal processing algorithm Method includes the following steps:
Step 1, satellite-borne SAR raw radar data and parameter profit are sent to FPGA cores by host computer and gigabit Ethernet Piece;
Step 2, echo data is subjected to sub-aperture division and calculates sub-aperture PFA processing procedures using parameter calculating module In the parameters such as required phase compensating factor;
In the step 2, echo data is first subjected to sub-aperture division first, the satellite-borne SAR data of big bandwidth are divided There is the sub-aperture of lap to be several and complete imaging, full aperture processing is avoided to cause the image aliasing of orientation;
Step 3, Estimation of Doppler central frequency is carried out according to echo data;
In the step 3, realize that dislocation is multiplied using shift register first and obtain auto-correlation function, then will be after multiplication Data deposit FIFO waits addition corresponding with next line result, and the data line finally obtained is carried out entirely using floating point accumulator Flowing water floating-point adds up, and the horizontal echo-signal for obtaining multiple range cells is averaged, then is ATAN by CORDIC cores and operates basis Formula obtains accurate doppler centroid.
Step 4, sub- pore size data is carried out based on the distance of dimensional variation to processing;
In the step 4, read satellite-borne SAR distance one by one first from DDR3SDRAM and carry out motion compensation to pulse, then Using PCS principles carry out it is a series of multiply again with FFT computings, often handled and waited in the write-in DDR of a pulse transposition Orientation processing;
Step 5, orientation high-precision Sinc interpolation is carried out;
In the step 5, after step 4 completion, by spaceborne data, Pulse by Pulse transposition is read from DDR first, will SAR64 double-precision floating point echo datas turn fixed point IP kernel by floating-point and are converted to fixed-point data, then will pinpoint echo data Integer part completes the weighting of echo data by addressing operation as address ram, piecemeal storage echo data and fractional part Summation obtains SINC interpolation results;
Step 6, bidimensional FFT imagings are carried out and data are sent back into host computer by Ethernet;
In the step 6, by the son of orientation treated data do bidimensional FFT imagings obtain satellite-borne SAR target scene Data are sent back host computer by aperture imaging result data by Ethernet;
Step 7, radiant correction processing and weighting splicing are done to sub- subaperture image in host computer, do after histogram equalization Host computer is shown.
In the step 7, the inverse for being multiplied by orientation antenna radiation pattern to orientation data field first in sub-aperture comes Influence of the antenna radiation pattern to the point target echo amplitude of different azimuth position in sub-aperture is eliminated, is spliced using weighted average and calculated Method realizes seamlessly transitting for overlapping region, obtains full aperture image, obtained full aperture image is converted into gray-scale map, does straight Last imaging results figure is shown in host computer after side's figure is balanced.
Above example is merely illustrative of the invention's technical idea, it is impossible to protection scope of the present invention is limited with this, it is every According to technological thought proposed by the present invention, any change done on the basis of technical solution each falls within the scope of the present invention Within.

Claims (6)

1. a kind of implementation method of the FPGA of satellite-borne SAR multi-modal imaging signal processing algorithm, it is characterised in that including walking as follows Suddenly:
Step 1, satellite-borne SAR raw radar data and parameter profit are sent to fpga chip by host computer and gigabit Ethernet;
Step 2, echo data is subjected to sub-aperture division, and calculates required parameter in sub-aperture PFA processing procedures;
Step 3, Estimation of Doppler central frequency is carried out according to echo data;
Step 4, sub- pore size data is carried out based on the distance of dimensional variation to processing;
Step 5, orientation high-precision Sinc interpolation is carried out;
Step 6, bidimensional FFT imagings are carried out and data are sent back into host computer by Ethernet;
Step 7, radiant correction processing and weighting splicing are done to sub- subaperture image in host computer, is done after histogram equalization upper Machine is shown.
2. a kind of implementation method of the FPGA of satellite-borne SAR multi-modal imaging signal processing algorithm as described in claim 1, special Sign is:In the step 2, when echo data is carried out sub-aperture division, the satellite-borne SAR data of big bandwidth are divided into several There is the sub-aperture of lap and complete imaging.
3. a kind of implementation method of the FPGA of satellite-borne SAR multi-modal imaging signal processing algorithm as described in claim 1, special Sign is:The detailed content of the step 3 is:It realizes that dislocation is multiplied using shift register first and obtains auto-correlation function, then Data deposit FIFO after multiplication is waited into addition corresponding with next line result, the data line finally obtained is tired out using floating-point Device is added to carry out full flowing water floating-point to add up, the horizontal echo-signal for obtaining multiple range cells is averaged, then does by CORDIC cores ATAN operates to obtain accurate doppler centroid.
4. a kind of implementation method of the FPGA of satellite-borne SAR multi-modal imaging signal processing algorithm as described in claim 1, special Sign is:The detailed content of the step 4 is:Satellite-borne SAR distance is read one by one first from DDR3SDRAM to be transported to pulse Dynamic compensation, then using PCS principles multiplied again with FFT computings, often handled and waited in the write-in DDR of a pulse transposition Orientation processing.
5. a kind of implementation method of the FPGA of satellite-borne SAR multi-modal imaging signal processing algorithm as described in claim 1, special Sign is:The detailed content of the step 5 is:By spaceborne data, Pulse by Pulse transposition is read from DDR first, double by SAR64 Accuracy floating-point echo data turns fixed point IP kernel by floating-point and is converted to fixed-point data, then the integer part for pinpointing echo data is made For address ram, piecemeal storage echo data and fractional part, the weighted sum that echo data is completed by addressing operation obtain SINC interpolation results.
6. a kind of implementation method of the FPGA of satellite-borne SAR multi-modal imaging signal processing algorithm as described in claim 1, special Sign is:The detailed content of the step 7 is:Orientation antenna radiation pattern is multiplied by orientation data field first in sub-aperture Inverse eliminate influence of the antenna radiation pattern to the point target echo amplitude of different azimuth position in sub-aperture, it is flat using weighting Equal stitching algorithm realizes seamlessly transitting for overlapping region, obtains full aperture image, obtained full aperture image is converted into ash Degree figure, last imaging results figure is shown in host computer after doing histogram equalization.
CN201711325281.8A 2017-12-13 2017-12-13 A kind of implementation method of the FPGA of satellite-borne SAR multi-modal imaging signal processing algorithm Pending CN108120980A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711325281.8A CN108120980A (en) 2017-12-13 2017-12-13 A kind of implementation method of the FPGA of satellite-borne SAR multi-modal imaging signal processing algorithm

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711325281.8A CN108120980A (en) 2017-12-13 2017-12-13 A kind of implementation method of the FPGA of satellite-borne SAR multi-modal imaging signal processing algorithm

Publications (1)

Publication Number Publication Date
CN108120980A true CN108120980A (en) 2018-06-05

Family

ID=62229848

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711325281.8A Pending CN108120980A (en) 2017-12-13 2017-12-13 A kind of implementation method of the FPGA of satellite-borne SAR multi-modal imaging signal processing algorithm

Country Status (1)

Country Link
CN (1) CN108120980A (en)

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108508426A (en) * 2018-07-03 2018-09-07 上海志良电子科技有限公司 SAR echo signal generation method based on multi-core DSP and echo simulator
CN109471097A (en) * 2018-11-14 2019-03-15 湖南华诺星空电子技术有限公司 A kind of through-wall radar Signal optimum processing method and device
CN109886941A (en) * 2019-01-31 2019-06-14 天津大学 SAR flood remote sensing imagery change detection method based on FPGA
CN110018459A (en) * 2019-04-11 2019-07-16 上海卫星工程研究所 Satellite-borne SAR echo auxiliary data process of analysis design method
CN110531355A (en) * 2019-07-04 2019-12-03 西安电子科技大学 A kind of satellite-borne synthetic aperture radar unifies imaging method in real time
CN110542900A (en) * 2019-10-12 2019-12-06 南京航空航天大学 SAR imaging method and system
CN111538005A (en) * 2020-05-13 2020-08-14 西安电子科技大学 SAR front-side-looking imaging method based on FPGA and multiple multi-core DSPs
CN111640147A (en) * 2020-06-04 2020-09-08 北京无线电测量研究所 SAR subimage registration method, device and computer equipment for step-by-step frequency band splicing
CN112346057A (en) * 2019-08-06 2021-02-09 北京华航无线电测量研究所 Two-dimensional interpolation module in image processing process
CN112946649A (en) * 2021-04-08 2021-06-11 电子科技大学 PFA imaging method suitable for any sub-aperture length
CN113376632A (en) * 2021-05-18 2021-09-10 南京航空航天大学 Large squint airborne SAR imaging method based on pretreatment and improved PFA
CN113484860A (en) * 2021-07-06 2021-10-08 中国科学院空天信息创新研究院 SAR image ocean front detection method and system based on Doppler center anomaly
CN114488148A (en) * 2022-01-12 2022-05-13 南京航空航天大学 Sparse TOPS-SAR-based imaging mode implementation method
CN118011397A (en) * 2024-04-08 2024-05-10 中国科学院空天信息创新研究院 Echo abnormality detection and compensation method

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101441271A (en) * 2008-12-05 2009-05-27 航天恒星科技有限公司 SAR real time imaging processing device based on GPU
CN103197317A (en) * 2013-03-26 2013-07-10 西安电子科技大学 Synthetic aperture radar (SAR) imaging method based on field programmable gate array (FPGA)
CN104007437A (en) * 2014-05-21 2014-08-27 西安电子科技大学 SAR real-time imaging processing method based on FPGA and multiple DSPs
CN104391297A (en) * 2014-11-17 2015-03-04 南京航空航天大学 Sub-aperture partition PFA (Polar Format Algorithm) radar imaging method
CN104793196A (en) * 2015-04-28 2015-07-22 西安电子科技大学 Real-time SAR (synthetic aperture radar) imaging method based on improved range migration algorithm
CN105866774A (en) * 2016-03-23 2016-08-17 南京航空航天大学 FPGA implementation method for polar coordinate format imaging algorithm of chirp signal
CN106896360A (en) * 2017-04-21 2017-06-27 南京航空航天大学 A kind of FPGA implementation method of SAR signal processing algorithms

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101441271A (en) * 2008-12-05 2009-05-27 航天恒星科技有限公司 SAR real time imaging processing device based on GPU
CN103197317A (en) * 2013-03-26 2013-07-10 西安电子科技大学 Synthetic aperture radar (SAR) imaging method based on field programmable gate array (FPGA)
CN104007437A (en) * 2014-05-21 2014-08-27 西安电子科技大学 SAR real-time imaging processing method based on FPGA and multiple DSPs
CN104391297A (en) * 2014-11-17 2015-03-04 南京航空航天大学 Sub-aperture partition PFA (Polar Format Algorithm) radar imaging method
CN104793196A (en) * 2015-04-28 2015-07-22 西安电子科技大学 Real-time SAR (synthetic aperture radar) imaging method based on improved range migration algorithm
CN105866774A (en) * 2016-03-23 2016-08-17 南京航空航天大学 FPGA implementation method for polar coordinate format imaging algorithm of chirp signal
CN106896360A (en) * 2017-04-21 2017-06-27 南京航空航天大学 A kind of FPGA implementation method of SAR signal processing algorithms

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
张晓等: ""基于***噪声抑制的ScanSAR辐射校正方法"", 《国外电子测量技术》 *
王绍清: ""星载SAR多普勒中心实时估计技术研究"", 《中国优秀博硕士学位论文全文数据库信息科技辑》 *

Cited By (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108508426B (en) * 2018-07-03 2021-11-30 上海志良电子科技有限公司 SAR echo signal generation method based on multi-core DSP and echo simulator
CN108508426A (en) * 2018-07-03 2018-09-07 上海志良电子科技有限公司 SAR echo signal generation method based on multi-core DSP and echo simulator
CN109471097A (en) * 2018-11-14 2019-03-15 湖南华诺星空电子技术有限公司 A kind of through-wall radar Signal optimum processing method and device
CN109471097B (en) * 2018-11-14 2022-04-22 湖南华诺星空电子技术有限公司 Through-wall radar signal optimization processing method and device
CN109886941A (en) * 2019-01-31 2019-06-14 天津大学 SAR flood remote sensing imagery change detection method based on FPGA
CN110018459A (en) * 2019-04-11 2019-07-16 上海卫星工程研究所 Satellite-borne SAR echo auxiliary data process of analysis design method
CN110018459B (en) * 2019-04-11 2023-06-09 上海卫星工程研究所 Method for designing analysis flow of satellite-borne SAR echo auxiliary data
CN110531355A (en) * 2019-07-04 2019-12-03 西安电子科技大学 A kind of satellite-borne synthetic aperture radar unifies imaging method in real time
CN112346057B (en) * 2019-08-06 2023-07-14 北京华航无线电测量研究所 Two-dimensional interpolation module in image processing process
CN112346057A (en) * 2019-08-06 2021-02-09 北京华航无线电测量研究所 Two-dimensional interpolation module in image processing process
CN110542900A (en) * 2019-10-12 2019-12-06 南京航空航天大学 SAR imaging method and system
CN111538005B (en) * 2020-05-13 2023-03-24 西安电子科技大学 SAR front-side-looking imaging method based on FPGA and multiple multi-core DSPs
CN111538005A (en) * 2020-05-13 2020-08-14 西安电子科技大学 SAR front-side-looking imaging method based on FPGA and multiple multi-core DSPs
CN111640147A (en) * 2020-06-04 2020-09-08 北京无线电测量研究所 SAR subimage registration method, device and computer equipment for step-by-step frequency band splicing
CN111640147B (en) * 2020-06-04 2023-10-31 北京无线电测量研究所 SAR sub-image registration method, device and computer equipment for step frequency band stitching
CN112946649A (en) * 2021-04-08 2021-06-11 电子科技大学 PFA imaging method suitable for any sub-aperture length
CN113376632A (en) * 2021-05-18 2021-09-10 南京航空航天大学 Large squint airborne SAR imaging method based on pretreatment and improved PFA
CN113376632B (en) * 2021-05-18 2023-12-15 南京航空航天大学 Large strabismus airborne SAR imaging method based on pretreatment and improved PFA
CN113484860B (en) * 2021-07-06 2023-07-14 中国科学院空天信息创新研究院 SAR image ocean front detection method and system based on Doppler center abnormality
CN113484860A (en) * 2021-07-06 2021-10-08 中国科学院空天信息创新研究院 SAR image ocean front detection method and system based on Doppler center anomaly
CN114488148A (en) * 2022-01-12 2022-05-13 南京航空航天大学 Sparse TOPS-SAR-based imaging mode implementation method
CN114488148B (en) * 2022-01-12 2024-03-29 南京航空航天大学 Sparse TOPS-SAR imaging mode-based implementation method
CN118011397A (en) * 2024-04-08 2024-05-10 中国科学院空天信息创新研究院 Echo abnormality detection and compensation method
CN118011397B (en) * 2024-04-08 2024-06-11 中国科学院空天信息创新研究院 Echo abnormality detection and compensation method

Similar Documents

Publication Publication Date Title
CN108120980A (en) A kind of implementation method of the FPGA of satellite-borne SAR multi-modal imaging signal processing algorithm
CN106896360A (en) A kind of FPGA implementation method of SAR signal processing algorithms
CN103197317B (en) Synthetic aperture radar (SAR) imaging method based on field programmable gate array (FPGA)
CN109581352B (en) Super-resolution angle measurement system based on millimeter wave radar
CN110542900B (en) SAR imaging method and system
US20070188371A1 (en) Spotlight synthetic aperture radar (SAR) system and method for generating a SAR map in real-time using a modified polar format algorithm (PFA) approach
CN106842210A (en) A kind of new multiple submatrixes synthetic aperture sonar fast imaging algorithm
CN102508251B (en) Method for rapidly implementing sector conversion in multi-beam image sonar
CN108614266A (en) A kind of implementation method of the FPGA of video SAR high-speed processing technologies
CN109444901A (en) A kind of multiple submatrixes SAS sub-aperture image method under isomerous environment
CN103257341A (en) Fast autofocus algorithm implementation method based on FPGA
CN109298417B (en) Building internal structure detection method and device based on radar signal processing
CN102565772B (en) Marine dynamic information extraction method on basis of SAR (Synthetic Aperture Radar) sub-aperture sequence images
Linchen et al. FPGA implementation of polar format algorithm for airborne spotlight SAR processing
CN113156435B (en) Missile-borne SAR front side view time domain imaging method based on embedded GPU
EP3444632A1 (en) Signal processing device and radar apparatus
CN108427111B (en) Radar ranging method and device
Xin et al. A wide-field SAR polar format algorithm based on quadtree sub-image segmentation
CN114626006B (en) FPGA (field programmable Gate array) realization method for real-time generation of CS (Circuit switched) algorithm compensation factor in radar imaging
Wang et al. Video SAR high-speed processing technology based on FPGA
CN114185047B (en) Double-base SAR moving target refocusing method based on optimal polar coordinate transformation
CN102654574A (en) Doppler centroid estimation method based on field programmable gate array (FPGA)
CN104166140B (en) Method and device for realizing inverse synthetic aperture radar imaging
Yang et al. ISAR cross-range scaling for non-uniform rotating targets using joint affine invariant normalization and weighted average method
Zhou et al. A High-Precision Motion Errors Compensation Method Based on Sub-Image Reconstruction for HRWS SAR Imaging

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20180605

RJ01 Rejection of invention patent application after publication