CN114942441A - Progressive scanning terrain observation mode scanning parameter determination method - Google Patents

Progressive scanning terrain observation mode scanning parameter determination method Download PDF

Info

Publication number
CN114942441A
CN114942441A CN202210856103.2A CN202210856103A CN114942441A CN 114942441 A CN114942441 A CN 114942441A CN 202210856103 A CN202210856103 A CN 202210856103A CN 114942441 A CN114942441 A CN 114942441A
Authority
CN
China
Prior art keywords
scanning
sub
swath
overlapping
calculating
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
Application number
CN202210856103.2A
Other languages
Chinese (zh)
Other versions
CN114942441B (en
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.)
Aerospace Information Research Institute of CAS
Original Assignee
Aerospace Information Research Institute of CAS
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 Aerospace Information Research Institute of CAS filed Critical Aerospace Information Research Institute of CAS
Priority to CN202210856103.2A priority Critical patent/CN114942441B/en
Publication of CN114942441A publication Critical patent/CN114942441A/en
Application granted granted Critical
Publication of CN114942441B publication Critical patent/CN114942441B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • 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
    • 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
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

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 invention provides a progressive scanning terrain observation mode scanning parameter determination method which specifically comprises the steps of inputting SAR system parameters, setting initial overlapping coefficients of azimuth direction adjacent observation bands in sub mapping bands, constructing an equation set containing the overlapping coefficients according to a star-earth geometric relation, and calculating parameters such as a scanning period, a scanning angle and scanning time of each sub mapping band according to the equation set. And then, according to the actual scanning capability, the stepping precision and the pulse repetition frequency of the antenna, practical adjustment is carried out on the calculation parameters to obtain actually usable working parameters such as a scanning angle, stepping times, stepping resident pulse number and the like. And then, calculating the overlapping degree of the azimuth direction of each sub mapping band to the actual observation band according to the actual working parameters, adjusting the overlapping coefficient when the overlapping degree does not meet the requirement, and recalculating the working parameters until the overlapping degree meets the requirement. The method is simple and easy to operate, and does not need to consider various factors and the difference between the working parameter approximate value and the theoretical calculated value.

Description

Progressive scanning terrain observation mode scanning parameter determination method
Technical Field
The invention relates to the technical field of Synthetic Aperture Radars (SAR), in particular to a method for determining a scanning parameter of a Progressive scanning Terrain Observation (TOPS) mode.
Background
The synthetic aperture radar has the characteristics of all-time and all-weather imaging, and has wide application in the fields of national economy and national defense. The TOPS mode is a high-resolution wide swath imaging mode, compared with a Scanning (SCAN) mode, the TOPS mode eliminates the scallop effect and obtains a uniform noise equivalent coefficient in the whole sub swath; the TOPS mode has stronger fuzzy inhibition capability, and the fuzzy degree in the whole sub-swath is more uniform; the loss of the TOPS mode edge effect is smaller due to the long duration of the azimuthal scan within the sub swath. These characteristics make TOPS mode obtain better image quality under the condition of not reducing resolution ratio, even improving resolution ratio, and become preferred scheme in high resolution ratio wide swath application.
The TOPS mode working parameter design is complex, and the influence of various factors, such as system indexes of resolution, swath width, sub swath number and the like, needs to be comprehensively considered. After the factors are determined, working parameters such as the scanning period, the scanning angle of each sub measuring and drawing strip, the stepping times, the number of resident pulses and the like are further determined. The working parameters obtained by theoretical calculation often have multi-digit decimal values, even infinite decimal values, and cannot be directly used in practice. For example, in practice, the maximum scanning angle of an antenna is limited by the scanning capability of the antenna, and the scanning precision is limited; the number of steps is an integer, and the number of resident pulses per step is also an integer. Therefore, practical adjustment of these operating parameters one by one is required, and such adjustment may cause incomplete splicing of the adjacent observation bands in the azimuth direction in the sub-swath. If the splicing requirement is satisfied, the factors need to be considered one by one, and the difference between the practical adjustment value and the theoretical calculation value is used as a time error for recalculation. And repeating for multiple times until the splicing requirement is met. This method is not only complex in calculation, but also increases the scanning angle, increases the antenna load, and increases the system cost. In addition, too large a scan angle may cause an azimuthal scallop effect, resulting in a degradation of image quality.
Disclosure of Invention
In order to solve the technical problems, the invention provides a method for determining scanning parameters of a progressive scanning terrain observation mode, which does not need to consider time difference caused by practical adjustment of multiple factors, only takes the overlapping coefficient of the azimuth direction of each sub-swath to the adjacent observation band as a design basis, minimizes the scanning angle while meeting image splicing, and ensures that the image quality is not reduced. Although the invention is designed by taking a satellite-borne platform as an example, the invention is also applicable to an airborne platform.
In order to achieve the purpose, the invention adopts the technical scheme that:
a method for determining scanning parameters of a progressive scanning terrain observation mode comprises the following steps:
step 1, inputting system parameters, wherein the system parameters comprise: the antenna scanning angular speed, the number of the sub-swaths, the beam width of each sub-swath, the beam center slant distance and the beam ground projection speed;
step 2, setting an initial overlapping coefficient of the azimuth sub mapping band;
step 3, constructing an equation set containing an overlapping coefficient according to the geometric relation of the satellite and the ground, and calculating a scanning period, a scanning angle of each sub-swath and scanning time according to the equation set;
and 4, according to the actual scanning capability, the stepping precision and the pulse repetition frequency of the antenna, practically adjusting the scanning period, the scanning angle and the scanning time of each sub-swath in the step 3 to obtain actually usable working parameters, wherein the working parameters comprise the scanning angle, the stepping times, the stepping resident pulse number, the actual scanning time, the effective scanning track of the wave beam and the overlapping length of the wave beam centers of two continuous scanning.
And 5, calculating the near-end azimuth actual scanning overlapping degree of each sub mapping band beam according to the actually usable working parameters obtained in the step 4, adjusting the overlapping coefficient when the overlapping degree does not meet the requirement, calculating the scanning period, the scanning angle and the scanning time of each sub mapping band according to the step 3, and adjusting the actually usable working parameters according to the step 4 until the overlapping degree meets the requirement.
Further, in the step 2, the initial overlapping coefficient is set to 1, i.e. not overlapped, and then the overlapping coefficient is increased to meet the requirement; or the initial overlap coefficient is set to a number greater than 1 and less than 2, and then the overlap coefficient is reduced to satisfy the requirement.
Further, the step 3 specifically includes:
introducing overlap factork gap And obtaining a sub swath expression:
(|ωT b -θ c R c +V s ·T b =k gap ·V s ·T R
wherein the content of the first and second substances,ωis the sub swath antenna scan angular velocity,T b is the sub swath azimuth scan duration,θ c is the sub-swath beam width,V s is the sub-swath beam ground projection velocity,R c is the center slant distance of the sub-swath,T R is the scan period;
determining the number of the distance direction sub mapping bands according to the resolution and the breadth required by the system, and assuming that the number of the sub mapping bands is n, obtaining the following equation set:
Figure 854919DEST_PATH_IMAGE001
wherein the superscript (-) denotes the child swath index,T R represents a scanning period;
solving the equation set to determine the azimuth scanning duration of the sub mapping zone m
Figure 369077DEST_PATH_IMAGE002
M represents the mth sub-swath, m =1,2, …, n;
by the formula
Figure 498707DEST_PATH_IMAGE003
Calculating the maximum scanning angle of the m-th sub-swath wave beam azimuth to further obtain the antenna scanning angle range
Figure 781921DEST_PATH_IMAGE004
Wherein, in the step (A),
Figure 604383DEST_PATH_IMAGE005
is the maximum scan angle of the antenna beam in the azimuth direction within each sub swath.
Further, the step 4 specifically includes: and determining parameters of each sub-swath, wherein the parameters comprise a scanning angle, stepping times, a resident pulse number, actual scanning time, an effective beam scanning track and the beam center overlapping length of two continuous scanning times.
Assuming an antenna stepping angle ofCalculating the actual maximum scanning angle of the antenna by adopting the following formula:
Figure 340258DEST_PATH_IMAGE006
ceil(. cndot.) represents an upward integer;
by the formula
Figure 8000DEST_PATH_IMAGE007
Calculating the beam stepping times in the sub-mapping band;
by the formula
Figure 411299DEST_PATH_IMAGE008
Counting the number of resident pulses per step, whereinprf (m) Is the pulse repetition frequency of the sub swath;
by the formula
Figure 404663DEST_PATH_IMAGE009
Calculating the actual scanning time of each sub measuring and drawing band;
by the formula
Figure 627834DEST_PATH_IMAGE010
Calculating an actual scanning period;
by the formula
Figure 364846DEST_PATH_IMAGE011
Calculating the effective scanning track of the mth sub swath wave beam, namely the effective imaging length in the azimuth direction;
by the formula
Figure 622652DEST_PATH_IMAGE012
And calculating the overlapping length of the observed region of the beam center of the mth sub-swath in two continuous scanning periods.
Further, the step 5 specifically includes:
by the formula
Figure 521338DEST_PATH_IMAGE013
Calculating the difference of the synthetic aperture lengths corresponding to the central view angle and the near view angle of the antenna beam in the mth sub mapping band;
by the formula
Figure 497384DEST_PATH_IMAGE014
Calculating the overlapping length of the observation region of the beam near end of the mth sub-swath in two continuous scanning periods;
judgment of
Figure 506928DEST_PATH_IMAGE015
Whether the overlap requirement is satisfied, i.e., whether the following equation is satisfied:
Figure 619241DEST_PATH_IMAGE016
wherein (A), (B), (C) and CS 1 ,S 2 ) Indicating the overlap range of the system requirements;
if it is not
Figure 954407DEST_PATH_IMAGE017
If the overlapping degree is not in the range, adjusting the overlapping coefficient, and repeating the steps 3 to 5 until the overlapping degree meets the requirement.
The invention has the beneficial effects that:
the invention provides a high-efficiency and quick scanning parameter determination method, which is used for controlling a beam to carry out azimuth scanning. The invention can minimize the beam azimuth scanning angle and improve the image quality while meeting the overlapping degree of azimuth adjacent observation bands in the sub-swath. Different from the conventional TOPS mode scanning parameter determination method, the method does not consider the time difference caused by practical adjustment of parameters such as the stepping times, the resident pulse number and the like, only takes the overlapping coefficient as a design basis, and can adjust the coefficient at any time according to the requirement of the overlapping degree of the system.
Drawings
Fig. 1 is a flow chart of the method for determining scanning parameters in a progressive scanning topography observation mode according to the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention. In addition, the technical features involved in the respective embodiments of the present invention described below may be combined with each other as long as they do not conflict with each other.
As shown in fig. 1, the method for determining scanning parameters in a progressive scanning terrain observation mode specifically includes:
step 1, inputting system parameters, wherein the system parameters comprise: the antenna scanning angular speed, the number of the sub-swaths, the beam width of each sub-swath, the beam center slant distance and the beam ground projection speed;
step 2, setting an initial overlapping coefficient of the azimuth direction sub mapping band;
step 3, constructing an equation set containing an overlapping coefficient according to the geometric relation of the satellite and the ground, and calculating a scanning period, a scanning angle of each sub-swath and scanning time according to the equation set;
step 4, according to the actual scanning capability, the stepping precision and the Pulse Repetition Frequency (PRF) of the antenna, the scanning period, the scanning angle and the scanning time of each sub-swath in the step 3 are practically adjusted to obtain actually usable working parameters, wherein the actually usable working parameters comprise the scanning angle, the stepping times, the stepping resident pulse number, the actual scanning time, the effective scanning track of the wave beam and the overlapping length of the central tracks of the wave beam scanned for two times continuously;
and step 5, calculating the overlapping degree of the near-end azimuth of each sub mapping band wave beam to the adjacent actual observation band according to the actually usable working parameters obtained in the step 4, when the overlapping degree does not meet the requirement, adjusting the overlapping coefficient, calculating the scanning period, the scanning angle and the scanning time of each sub mapping band again according to the step 3, and adjusting the actually usable working parameters according to the step 4 until the overlapping degree meets the requirement.
Specifically, the invention comprises the following steps:
step 1, inputting system parameters, and inputting system parameters required by calculation, wherein the system parameters comprise: the antenna scanning angular velocity, the number of the sub-swaths, the beam width of each sub-swath, the beam center slant distance, the beam near-end slant distance, the beam ground projection speed and the like.
And 2, setting an initial overlapping coefficient of the azimuth sub mapping band.
The initial overlap factor may be set to 1, i.e. no overlap, and then the requirement is met by increasing the overlap factor; it can also be set to a larger value, such as 1.1, and then the overlap factor is reduced to meet the requirement, and the specific value can be determined according to the actual requirement.
And 3, constructing an equation set containing the overlapping coefficient according to the satellite-ground geometric relation, and calculating the scanning period, the scanning angle of each sub measuring and drawing band and the scanning time according to the equation set.
The theoretical calculation equation of the conventional TOPS mode has no overlapping coefficient, so that the invention introduces the overlapping coefficient on the basis of the conventional equation systemk gap From the geometrical relationships, by strict derivation, we obtain:
(|ωT b -θ c R c +V s ·T b =k gap ·V s ·T R
wherein the content of the first and second substances,ωis the sub swath antenna scan angular velocity,T b is the sub swath azimuth scan duration,θ c is the sub-swath beam width,V s is the sub-swath beam ground projection velocity,R c is the center slant distance of the sub-swath,T R is the scan period.
Determining the number of the sub mapping bands of the distance direction according to the resolution and the width required by the system, and assuming that the number of the sub mapping bands is n, obtaining an equation set as follows:
Figure 152170DEST_PATH_IMAGE018
wherein the superscript (-) denotes the child swath index,T R representing the scan period.
Solving the equation set can determine the m azimuth scanning duration of the sub mapping band
Figure 230985DEST_PATH_IMAGE019
M denotes the mth subpologram, m =1,2, …, n.
Then by the formula
Figure 463383DEST_PATH_IMAGE020
Calculating the azimuth maximum scanning angle of the mth sub-swath to further obtain the scanning angle range of the antenna
Figure 703872DEST_PATH_IMAGE021
. Wherein the content of the first and second substances,
Figure 654510DEST_PATH_IMAGE005
is the maximum scan angle of the antenna beam in the azimuth direction within each sub swath.
Step 4, according to the actual scanning capability, the stepping precision and the Pulse Repetition Frequency (PRF) of the antenna, the scanning period, the scanning angle and the scanning time of each sub-swath in the step 3 are practically adjusted to obtain actually usable working parameters, wherein the actually usable working parameters comprise the scanning angle, the stepping times, the stepping resident pulse number, the actual scanning time, the effective scanning track of the wave beam and the overlapping length of the central tracks of the wave beam scanned for two times continuously;
the phased array antenna adopts a step scanning mode, beam continuous scanning cannot be realized, the step precision is limited by antenna beam control precision, any angle cannot be scanned, the number of the stepping resident pulses at each time is an integer value, and any value cannot be obtained, so that practical adjustment needs to be carried out on the calculation result.
Assuming an antenna stepping angle ofCalculating the actual maximum scanning angle of the mth sub swath antenna by adopting the following formula:
Figure 271436DEST_PATH_IMAGE022
ceil(. cndot.) denotes an upward integer.
By the formula
Figure 358341DEST_PATH_IMAGE023
And calculating the beam stepping times in the mth mapping sub-band.
By the formula
Figure 766801DEST_PATH_IMAGE024
Calculating the number of resident pulses per step, whereinprf (m) Is the pulse repetition frequency of the sub swath.
By the formula
Figure 939156DEST_PATH_IMAGE025
And calculating the actual scanning time of the mth sub-swath.
By the formula
Figure 359773DEST_PATH_IMAGE026
The actual scan period is calculated.
By the formula
Figure 301185DEST_PATH_IMAGE027
Calculating the effective scanning track of the mth sub swath wave beam, namely the effective imaging length in the azimuth direction; wherein the content of the first and second substances,L m() is the m-th sub-mapping band, the antenna beam is along the azimuth-wise effective imaging length.
By the formula
Figure 883476DEST_PATH_IMAGE028
And calculating the overlapping length of the observed region of the beam center of the mth sub-swath in two continuous scanning periods.
And 5, calculating the actual overlapping degree of the near-end azimuth of each sub mapping band beam according to the actually usable working parameters obtained in the step 4, and judging the overlapping degree.
And judging whether the overlapping degree of each sub mapping band in the azimuth direction to the adjacent observation band meets the requirement. Because the near-end synthetic aperture length of the beam is shorter than the synthetic aperture lengths of the center and the far end of the beam, in order to ensure that the near end has no splicing gap, the overlapping degree of the sub mapping bands in the near-end direction to the adjacent observation bands needs to be calculated.
By the formula
Figure 543127DEST_PATH_IMAGE029
And calculating the difference of the synthetic aperture lengths corresponding to the central visual angle and the near visual angle of the antenna beam in the mth sub mapping band.
By the formula
Figure 501856DEST_PATH_IMAGE030
And calculating the overlapping length of the observation region in the m-th sub mapping band beam near end in two continuous scans.
Judgment of
Figure 563353DEST_PATH_IMAGE031
Whether the overlap requirement is satisfied, i.e., whether the following equation is satisfied:
Figure 582124DEST_PATH_IMAGE032
wherein (A), (B), (C), (D), (C), (B), (C)S 1 ,S 2 ) Indicating the range of overlap required by the system.
If it is not
Figure 463493DEST_PATH_IMAGE031
If the overlapping degree is not in the range, adjusting the overlapping coefficient, and repeating the steps 3 to 5 until the overlapping degree meets the requirement.
The present invention relates to the overlapping coefficientk gap And when the current parameter is more than or equal to 2, the method is equivalent to multi-view processing, so that the method is suitable for calculating the multi-view parameters in the TOPS mode.
And determining actual working parameters according to the calculation result, controlling the SAR system and realizing a TOPS mode of high-resolution repeated observation.
It will be understood by those skilled in the art that the foregoing is only a preferred embodiment of the present invention, and is not intended to limit the invention, and that any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the scope of the present invention.

Claims (5)

1. A method for determining scanning parameters of a progressive scanning terrain observation mode is characterized by comprising the following steps:
step 1, inputting system parameters, wherein the system parameters comprise: the antenna scanning angular speed, the number of the sub-swaths, the beam width of each sub-swath, the beam center slant distance and the beam ground projection speed;
step 2, setting initial overlapping coefficients of adjacent azimuth observation zones;
step 3, constructing an equation set containing an overlapping coefficient according to the geometric relation of the satellite and the ground, and calculating a scanning period, a scanning angle of each sub-swath and scanning time according to the equation set;
step 4, according to the actual scanning capability, the stepping precision and the pulse repetition frequency of the antenna, carrying out practical adjustment on the scanning period, the scanning angle and the scanning time of each sub-surveying and mapping band in the step 3 to obtain actually usable working parameters, wherein the working parameters comprise the scanning angle, the stepping times, the stepping resident pulse number, the actual scanning time, the effective scanning track of the wave beam and the near-end overlapping length of the wave beam of two adjacent times of scanning;
and 5, calculating the overlapping degree of the near-end direction of the wave beam of each sub mapping band to the actual observation band according to the actually usable working parameters obtained in the step 4, adjusting the overlapping coefficient when the overlapping degree does not meet the requirement, calculating the scanning period, the scanning angle and the scanning time of each sub mapping band according to the step 3, and adjusting the actually usable working parameters according to the step 4 until the overlapping degree meets the requirement.
2. The progressive scan topography observation mode scan parameter determination method according to claim 1, wherein: in the step 2, the initial overlapping coefficient is set to be 1, namely, the overlapping coefficient is not overlapped, and then the overlapping coefficient is increased to meet the requirement; or the initial overlap coefficient is set to a number greater than 1 and less than 2, and then the overlap coefficient is reduced to satisfy the requirement.
3. A progressive scan topography observation mode scan parameter determination method according to claim 2, wherein: the step 3 specifically includes:
introducing overlap factork gap And obtaining:
(|ωT b -θ c R c +V s ·T b =k gap ·V s ·T R
wherein the content of the first and second substances,ωis the sub swath antenna scan angular velocity,T b is the sub swath azimuth scan duration,θ c is the sub-swath beam width,V s is the sub-swath beam ground projection velocity,R c is the center slant distance of the sub-swath,T R is the scan period;
determining the number of the distance direction sub mapping bands according to the resolution and the breadth required by the system, and assuming that the number of the sub mapping bands is n, obtaining the following equation set:
Figure 307937DEST_PATH_IMAGE001
wherein the superscript (-) denotes the child swath index,T R represents a scanning period;
solving the above equation set to determine the azimuth scanning duration of the sub mapping band m
Figure 744734DEST_PATH_IMAGE002
M represents the mth sub-swath, m =1,2, …, n;
by the formula
Figure 575287DEST_PATH_IMAGE003
Calculating the maximum scanning angle of the m-th sub-swath antenna azimuth to further obtain the antenna scanning angle range
Figure 755733DEST_PATH_IMAGE004
Wherein, in the step (A),
Figure 620921DEST_PATH_IMAGE005
is the maximum scan angle of the antenna beam in the azimuth direction within each sub swath.
4. A progressive scan topographical observation mode scan parameter determination method as claimed in claim 3, wherein: the step 4 specifically includes:
assuming an antenna stepping angle ofCalculating the actual maximum scanning angle of the mth sub-swath antenna by adopting the following formula:
Figure 494199DEST_PATH_IMAGE006
ceil(. cndot.) represents an upward integer;
by the formula
Figure 812048DEST_PATH_IMAGE007
Calculating the scanning stepping times of the wave beam in the mth sub-mapping band;
by the formula
Figure 530605DEST_PATH_IMAGE008
Calculating the number of resident pulses of the mth sub-swath in each step, whereinprf (m) Is the pulse repetition frequency of the mth sub swath;
by the formula
Figure 250299DEST_PATH_IMAGE009
Calculating the actual scanning time of the mth sub-swath;
by the formula
Figure 294479DEST_PATH_IMAGE010
Calculating an actual scanning period;
by the formula
Figure 99624DEST_PATH_IMAGE011
Calculating the effective scanning track of the mth sub swath wave beam, namely the effective imaging length in the azimuth direction;
by the formula
Figure 887451DEST_PATH_IMAGE012
And calculating the azimuth overlapping length of the mth sub-swath beam center in the observation swath in two continuous scanning periods.
5. The progressive scan topography observation mode scan parameter determination method according to claim 4, wherein: the step 5 specifically includes:
by the formula
Figure 458722DEST_PATH_IMAGE013
Calculating the difference of the synthetic aperture lengths corresponding to the central view angle and the near view angle of the antenna beam in the mth sub mapping band;
by the formula
Figure 673803DEST_PATH_IMAGE014
Calculating the overlapping length of the proximal azimuth direction of the observation region scanned twice in the mth sub-mapping band;
judgment of
Figure 700665DEST_PATH_IMAGE015
Whether the overlap requirement is met, i.e. whether the following equation is met:
Figure 292183DEST_PATH_IMAGE016
wherein (A), (B), (C), (D), (C), (B), (C)S 1 ,S 2 ) Indicating the overlap range of the system requirements;
if it is not
Figure 720890DEST_PATH_IMAGE015
If the overlapping degree is not in the range, adjusting the overlapping coefficient, and repeating the steps 3 to 5 until the overlapping degree meets the requirement.
CN202210856103.2A 2022-07-21 2022-07-21 Progressive scanning terrain observation mode scanning parameter determination method Active CN114942441B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210856103.2A CN114942441B (en) 2022-07-21 2022-07-21 Progressive scanning terrain observation mode scanning parameter determination method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210856103.2A CN114942441B (en) 2022-07-21 2022-07-21 Progressive scanning terrain observation mode scanning parameter determination method

Publications (2)

Publication Number Publication Date
CN114942441A true CN114942441A (en) 2022-08-26
CN114942441B CN114942441B (en) 2022-10-21

Family

ID=82910933

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210856103.2A Active CN114942441B (en) 2022-07-21 2022-07-21 Progressive scanning terrain observation mode scanning parameter determination method

Country Status (1)

Country Link
CN (1) CN114942441B (en)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101344594A (en) * 2008-05-14 2009-01-14 中国科学院上海光学精密机械研究所 Scanning synthetic aperture laser imaging radar
US20100207808A1 (en) * 2007-07-04 2010-08-19 Deutsches Zentrum für Luft- und Raumfahrt e.V. Method for processing TOPS (Terrain Observation by Progressive Scan)-SAR (Synthetic Aperture Radar)-Raw Data
CN102346249A (en) * 2010-07-28 2012-02-08 中国科学院电子学研究所 Implementation method for wide swath earth observation step scanning mode of synthetic aperture radar
CN102967861A (en) * 2012-10-17 2013-03-13 中国人民解放军国防科学技术大学 Terrain observation by progressive scans synthetic aperture radar (TOPSAR) system parameter engineering design method
CN104635230A (en) * 2015-03-05 2015-05-20 北京航空航天大学 Method for MIMO (multi-input multi-output)-SAR (synthetic aperture radar) near field measurement imaging azimuth side lobe suppression
CN105068070A (en) * 2015-07-08 2015-11-18 中国科学院电子学研究所 Method, device and system of scanning ground moving target indication by satellite-borne SAR
CN107783156A (en) * 2017-10-12 2018-03-09 深圳市华讯方舟空间信息产业科技有限公司 Satellite tracking method
CN114488148A (en) * 2022-01-12 2022-05-13 南京航空航天大学 Sparse TOPS-SAR-based imaging mode implementation method
CN114527443A (en) * 2022-01-25 2022-05-24 上海卫星工程研究所 Ka-band DBF-SAR real-time processing system and method

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100207808A1 (en) * 2007-07-04 2010-08-19 Deutsches Zentrum für Luft- und Raumfahrt e.V. Method for processing TOPS (Terrain Observation by Progressive Scan)-SAR (Synthetic Aperture Radar)-Raw Data
CN101344594A (en) * 2008-05-14 2009-01-14 中国科学院上海光学精密机械研究所 Scanning synthetic aperture laser imaging radar
CN102346249A (en) * 2010-07-28 2012-02-08 中国科学院电子学研究所 Implementation method for wide swath earth observation step scanning mode of synthetic aperture radar
CN102967861A (en) * 2012-10-17 2013-03-13 中国人民解放军国防科学技术大学 Terrain observation by progressive scans synthetic aperture radar (TOPSAR) system parameter engineering design method
CN104635230A (en) * 2015-03-05 2015-05-20 北京航空航天大学 Method for MIMO (multi-input multi-output)-SAR (synthetic aperture radar) near field measurement imaging azimuth side lobe suppression
CN105068070A (en) * 2015-07-08 2015-11-18 中国科学院电子学研究所 Method, device and system of scanning ground moving target indication by satellite-borne SAR
CN107783156A (en) * 2017-10-12 2018-03-09 深圳市华讯方舟空间信息产业科技有限公司 Satellite tracking method
CN114488148A (en) * 2022-01-12 2022-05-13 南京航空航天大学 Sparse TOPS-SAR-based imaging mode implementation method
CN114527443A (en) * 2022-01-25 2022-05-24 上海卫星工程研究所 Ka-band DBF-SAR real-time processing system and method

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
CHARLOTTE BAY HASAGER 等: "Mapping Offshore Winds Around Iceland Using Satellite Synthetic Aperture Radar and Mesoscale Model Simulations", 《IEEE JOURNAL OF SELECTED TOPICS IN APPLIED EARTH OBSERVATIONS AND REMOTE SENSING》 *
GANG ZHANG 等: "Accurate Elimination Method of 3-D Acceleration for Highly Squinted TOPS SAR with Curved Trajectory", 《 2021 2ND CHINA INTERNATIONAL SAR SYMPOSIUM (CISS)》 *
张群等: "步进频率逆合成孔径雷达成像的一种运动补偿方法", 《西安电子科技大学学报》 *
陈粤 等: "星载SAR方位多通道TOPS成像模式参数设计", 《雷达科学与技术》 *
高贺利 等: "一种宽幅星载TOPSAR数据高效成像算法", 《南京信息工程大学学报》 *

Also Published As

Publication number Publication date
CN114942441B (en) 2022-10-21

Similar Documents

Publication Publication Date Title
CN102346249B (en) Implementation method for wide swath earth observation step scanning mode of synthetic aperture radar
CN110632594B (en) Long-wavelength spaceborne SAR imaging method
CN108459321B (en) Large squint high-resolution SAR imaging method based on distance-azimuth circle model
CN109521424B (en) Satellite-borne sliding spotlight SAR satellite attitude and PRF sequence design method
CN111257878B (en) Wave form design method based on pitching dimensional frequency intra-pulse scanning high-resolution wide-range SAR
CN112034459B (en) Linear synthetic aperture radar imaging method and system based on antenna pattern compensation
CN114545411B (en) Polar coordinate format multimode high-resolution SAR imaging method based on engineering realization
US5334980A (en) Method and system for sharpening impulse response in a synthetic aperture radar
CN110673144A (en) Sub-aperture large squint SAR imaging processing method based on time-varying labels
CN108226925A (en) A kind of overlapped sub-aperture algorithm suitable for missile-borne strabismus time-varying parameter SAR before big
CN108107432B (en) High-low orbit bistatic SAR phase-preserving imaging method based on time domain disturbance
CN110879391B (en) Radar image data set manufacturing method based on electromagnetic simulation and missile-borne echo simulation
CN114942441B (en) Progressive scanning terrain observation mode scanning parameter determination method
US6650274B1 (en) Radar imaging system and method
CN113030970A (en) Two-dimensional scanning high-resolution wide-range SAR waveform design method based on azimuth frequency scanning
CN113156437A (en) Method for evaluating influence of high-orbit SAR on radio frequency interference of low-orbit SAR imaging
CN114924269B (en) Distance ambiguity analysis method based on spaceborne F-SCAN SAR
CN112198490A (en) Ground clutter suppression method and device for airborne millimeter wave cloud detection radar and radar
CN115436940B (en) Sparse sliding bunching SAR imaging mode implementation method and device
CN107607950B (en) Satellite-borne Ka InSAR beam forming method based on DBF-TOPS weighting
CN110579762A (en) Terahertz circular track SAR rapid back projection imaging method
CN112986989B (en) Method for restraining distance ambiguity of orthogonal phase coding signal based on genetic algorithm
KR102408991B1 (en) High resolution image decoding system based on squint sar and method of decoding image using the same
Zhang et al. Multi-beam Doppler beam sharpening approach for airborne forward-looking radar imaging
CN110297240B (en) Imaging method and device of azimuth wide-beam 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