CN106872979A - Ocean wave parameter acquisition methods based on shifting target SAR image refocusing - Google Patents
Ocean wave parameter acquisition methods based on shifting target SAR image refocusing Download PDFInfo
- Publication number
- CN106872979A CN106872979A CN201710120632.5A CN201710120632A CN106872979A CN 106872979 A CN106872979 A CN 106872979A CN 201710120632 A CN201710120632 A CN 201710120632A CN 106872979 A CN106872979 A CN 106872979A
- Authority
- CN
- China
- Prior art keywords
- motion
- phase error
- wave
- compensation
- curve
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 230000033001 locomotion Effects 0.000 claims abstract description 91
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 39
- 230000014509 gene expression Effects 0.000 claims abstract description 18
- 238000012545 processing Methods 0.000 claims abstract description 18
- 230000009466 transformation Effects 0.000 claims abstract description 9
- 238000004458 analytical method Methods 0.000 claims abstract description 3
- 238000012937 correction Methods 0.000 claims description 16
- 238000013508 migration Methods 0.000 claims description 11
- 230000005012 migration Effects 0.000 claims description 11
- 238000005096 rolling process Methods 0.000 claims description 5
- 238000011478 gradient descent method Methods 0.000 claims description 4
- 230000001133 acceleration Effects 0.000 claims description 2
- 230000006835 compression Effects 0.000 claims description 2
- 238000007906 compression Methods 0.000 claims description 2
- 230000005484 gravity Effects 0.000 claims description 2
- 238000012544 monitoring process Methods 0.000 abstract description 5
- 238000005259 measurement Methods 0.000 abstract description 3
- 238000007405 data analysis Methods 0.000 abstract description 2
- 230000008569 process Effects 0.000 description 6
- 238000001514 detection method Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 238000005457 optimization Methods 0.000 description 5
- 238000009795 derivation Methods 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 238000007667 floating Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000002592 echocardiography Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000005314 correlation function Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9021—SAR image post-processing techniques
- G01S13/9027—Pattern recognition for feature extraction
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9021—SAR image post-processing techniques
- G01S13/9029—SAR image post-processing techniques specially adapted for moving target detection within a single SAR image or within multiple SAR images taken at the same time
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention provides a kind of ocean wave parameter acquisition methods based on shifting target SAR image refocusing,By the SAR image for obtaining target area,The refocusing algorithm of inverse transformation and ISAR motion compensation is based on to the data application comprising moving target subgraph in the SAR image,The corresponding phase error compensation curve of SAR image and the image after being focused on,Phase error expression formula in the phase error compensation matched curve combination the heave of waves motion model of SIN function form is carried out into correspondence analysis,Obtain motion amplitude and the cycle of undulatory motion target,And the motion amplitude according to undulatory motion target deduces the parameters such as wave height and the cycle of wave with the cycle,Processing procedure is simple and quick,Cost is relatively low,Result is accurate,And above-mentioned Data Analysis Services advantage is had based on the method,Only need to place light keg buoy in measured zone,Monitoring mode is simply easily implemented,Can be with accurate measurement inshore information,Floor space is small,And do not limited by marine site.
Description
Technical Field
The invention belongs to the field of ocean remote sensing information extraction and the field of synthetic aperture radar signal processing, in particular to the field of ocean remote sensing sea wave information extraction and synthetic aperture radar data fine processing, and relates to a sea wave parameter acquisition method based on sea surface fluctuation target SAR image refocusing.
Background
In the research of ocean dynamic environment, sea waves are an important element and an important parameter in oceanographic research, so that the prediction and the acquisition of the information of the sea waves have a very close relationship with ocean military technology, safe sailing at sea and the like, and have very important scientific significance and practical value.
The existing ocean monitoring means mainly comprise the traditional modes of buoys, submerged buoys and the like and the emerging radio remote measuring modes of satellites, lasers, radars and the like. The traditional monitoring means has limited detection area range, can only detect data on a plurality of points and lines, is easily limited by severe climates such as high sea conditions, typhoons, tsunamis and the like, and has low working efficiency. Various emerging radio telemetry approaches have been rapidly developed to effectively avoid this limitation. Satellite telemetry has the advantages of global coverage, but the detection period is longer and the resolution is lower; the microwave shore-based radar has high resolution, can provide more comprehensive sea wave and ocean current information, has high precision, is widely deployed in places such as ship radars, shore-based ports, important navigation channels and the like at present, but has limited detection distance due to large propagation attenuation, can generally only detect within a plurality of kilometers, and is difficult to meet medium-long distance marine environment monitoring; the passive radar based on GNSS developed in recent years can detect sea surface wind field and wave field, but the detection distance and area are limited, so that the passive radar based on GNSS has certain limitation; the shore-based high-frequency ground wave radar can overcome the difficulties, can provide large-area data and continuous time series of sea wave, ocean current and sea surface wind speed data, but has high manufacturing cost and large antenna occupied area, and more importantly, has larger blind areas for near-shore areas and is difficult to obtain important information of the sea wave within several kilometers of the near-shore areas.
Therefore, in general, the following technical problems still exist in the field of ocean remote sensing sea wave information extraction at present: the processing means is complicated, the near-coast information cannot be accurately obtained, the manufacturing cost is high, the occupied area of the antenna is large, the detection distance is limited, and the like.
Disclosure of Invention
Technical problem to be solved
The invention provides a sea wave parameter acquisition method based on sea surface fluctuation target SAR image refocusing, which at least partially solves the technical problems.
(II) technical scheme
According to one aspect of the invention, a sea wave parameter acquisition method based on sea surface fluctuation target SAR image refocusing is provided, and comprises the following steps: acquiring an SAR image of a target area, and applying a refocusing algorithm based on inverse transformation and ISAR motion compensation to data containing a moving target sub-image in the SAR image to obtain a focused SAR image and a phase error compensation curve corresponding to the image; fitting a sine function to the phase error compensation curve to obtain a phase error compensation fitting curve in the form of the sine function; correspondingly analyzing the phase error compensation fitting curve in the form of the sine function by combining a phase error expression in a wave fluctuation motion model to obtain the motion amplitude and the cycle of a fluctuation motion target; and deducing parameters such as wave height, wave period and the like of the sea waves according to the motion amplitude and the motion period of the fluctuating motion target.
(III) advantageous effects
According to the technical scheme, the sea wave parameter acquisition method based on sea surface fluctuation target SAR image refocusing provided by the invention has at least one of the following beneficial effects:
1. the method establishes an error model of the fluctuating moving target, obtains an error result through quantitative analysis, innovatively adopts a refocusing method based on inverse transformation and ISAR motion compensation to realize sea wave parameter acquisition, and can more accurately and more simply obtain the sea wave parameter acquisition result while completing the refocusing of the moving target.
2. By means of fine processing of SAR image data, wave height and period parameter acquisition results of sea waves can be obtained, the processing process is simple and rapid, the cost is low, the results are accurate, the method has the advantages of data analysis and processing, small buoys only need to be placed in a measurement area, the monitoring mode is simple and easy to implement, and limitation of the measurement area is avoided.
Drawings
Fig. 1 is a technical flowchart of a sea wave parameter acquisition method based on sea surface fluctuation target SAR image refocusing.
Fig. 2 is a schematic diagram of a heave motion target point-satellite-earth geography.
Fig. 3 is a flowchart of a refocusing processing algorithm based on inverse transformation and ISAR motion compensation.
FIG. 4 is a flowchart of a distance alignment algorithm processing in the ISAR algorithm.
Fig. 5 is a flowchart of the range-phase correction process in the ISAR algorithm.
Fig. 6 is an SAR image of the offshore area of a certain island in China voyage under a Terra-SAR sliding bunching high-resolution mode.
Fig. 7 shows an optical image corresponding to fig. 6 for google earth.
FIG. 8 is a graph showing a comparison of the results of the experiment before and after the refocusing treatment.
Fig. 9 is a phase error compensation curve diagram corresponding to a good focusing effect.
Detailed Description
The invention provides a sea surface fluctuation target SAR image refocusing-based sea surface parameter acquisition method, which realizes accurate acquisition of sea surface parameters by adopting an algorithm based on inverse transformation and ISAR motion compensation and refocusing.
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 in conjunction with specific embodiments.
In one embodiment of the invention, a sea wave parameter acquisition method based on sea surface fluctuation target SAR image refocusing is provided. Fig. 1 is a technical flowchart of a sea wave parameter acquisition method based on sea surface fluctuation target SAR image refocusing. Referring to fig. 1, the sea wave parameter acquisition method based on sea surface fluctuation target SAR image refocusing comprises the following steps:
step A: establishing a model of wave fluctuation motion, and quantitatively analyzing an inclined distance error caused by the model, thereby obtaining a phase error expression caused by the inclined distance error;
the heave motion expression:
slope distance error:
phase error:
whereinA vector representing the heave motion of the target point, Δ R (η) representing the slope error,representing the modelled derived phase error, AvRepresenting the amplitude of the heave motion, λ representing the pulse wavelength, γ representing the sum of the downward viewing angle and the earth's center angle, ω representing the period of the heave motion, η representing any time instant corresponding to the heave motion,indicates the initial phase angle.
A small buoy is placed in an observation sea area to serve as a target point for the fluctuation motion of sea waves. Fig. 2 is a schematic view of observation target points-satellite-earth geography, as shown in fig. 2, T represents a stationary target point, S represents a satellite position of η at any azimuth, and S' represents a satellite subsatellite point at the time. Assuming that T has undulations in the form of sinusoidal vibrations along the radial direction of the earth, P represents the instantaneous position of the target point undulation motion.
Based on the above, the derivation process of the 3 formulas in step a is as follows:
vector representing the undulation motion of the target point:
the instantaneous slope distance vector between the static target and the satellite at the moment isTherefore, the instantaneous slope distance of the undulating target at the moment can be obtained as follows:
by applying the geographical relation and the vector knowledge, the approximate instantaneous slope distance of the fluctuating target can be obtained as follows:
then it can be deduced that the instantaneous slope error caused by the rolling object motion is:
the phase error caused by the instantaneous slope error is:
as can be seen from the phase error expression of equation (3): the phase error caused by the heave motion is in a simple harmonic vibration mode, and the phase error caused by the heave motion and the heave motion have the same period as seen by referring to a target point heave motion vector expression, and the amplitude of the phase error is determined by the heave motion amplitude, the downward viewing angle, the geocentric angle between the target and the satellite and the pulse wavelength.
After the sea wave motion target model is established in the step A to derive the phase error expression, the phase error expression can be directly called in the subsequent step D execution process without repeated derivation. Of course, step a may be omitted if the model is known and the phase error expression is also known.
Particularly, when analyzing the wave information of different areas, the wave fluctuation motion model established in the step a is established once, and modeling is not required to be carried out again each time, and the deduced error expression can be directly called when analyzing the wave information of different areas subsequently.
And B: acquiring an SAR image of a target area, and applying a refocusing algorithm based on inverse transformation and ISAR motion compensation to data containing a moving target sub-image in the SAR image to obtain a focused SAR image and a phase error compensation curve corresponding to the image;
fig. 3 is a flowchart of a refocusing processing algorithm based on Inverse transform and ISAR motion compensation, and as shown in fig. 3, the step B of applying an Inverse transform and ISAR (Inverse SAR/Inverse synthetic aperture radar) based motion compensation refocusing algorithm to the target data of the heave motion further includes:
substep B1: extracting a subimage containing a fluctuating moving target from an SAR image obtained by acquiring information in a target area;
substep B2: performing azimuth Fourier transform on the sub-image containing the fluctuating moving target to obtain equivalent original data before azimuth compression;
the data can be regarded as azimuth unfocused data after completely compensating the motion between the stationary target and the satellite, and can be used as raw data of input ISAR processing focusing, and the distance frequency domain can be equivalent to the following formula:
wherein f is0Representing the carrier frequency, B representing the transmission pulse bandwidth, TsRepresenting the imaging time, R (η) represents the instantaneous slope distance between the target and the satellite due to the heave motion of the target:
R(η)＝R0+ΔR(η) (7)
substep B3: performing refocusing processing on the equivalent original data obtained in the step B2 by adopting a distance alignment algorithm to correct a distance migration error caused by the fluctuation motion to obtain a distance migration error correction curve;
fig. 4 is a processing flow chart of a distance alignment algorithm in the ISAR algorithm, and as shown in fig. 4, the distance alignment algorithm adopts a global minimum entropy distance alignment algorithm, performs 32-fold interpolation on distance direction data to obtain distance alignment accuracy better than 32-fold distance resolution, and can completely correct distance migration errors caused by rolling motion.
The minimum entropy distance alignment algorithm applied to the orientation uncompressed data in fig. 4 achieves the distance-to-envelope alignment of the target, and eliminates the dislocation of adjacent echoes generated by motion in the distance direction without providing a priori knowledge. Establishing a model of the optimization problem, wherein the objective function is entropy values of all one-dimensional range image sums, and can be expressed as:
pave(r) represents the sum of all target one-dimensional range images:
where r is the range-oriented sample point and Δ r (n) represents the range offset of the nth echo.
paveThe degree of sharpening of (r) may be used to measure the degree of alignment of the echo envelopes, and when the degree of sharpening reaches the highest, the echoes are substantially aligned. And the Shannon entropy can be used to measure pave(r) degree of sharpening, the higher the waveform degree of sharpening, the smaller the entropy value, and therefore the optimization criterion is paveThe entropy value of (r) is minimum, and a global optimal solution of delta r (n) can be obtained through iteration, and the specific algorithm steps are described as follows:
step S401: initializing an offset Δ r (n);
step S402: calculate the current pave(r)；
Step S403: calculating pave(r), if the entropy value is not reduced, obtaining a final estimated value of the offset delta r (n), and exiting the loop algorithm, otherwise, executing the step S404;
step S404: computingIs calculated according to the following formula:
step S405: using FFT to convertCorrelating with each pulse;
step S406: taking the delta r (n) when the absolute value of the cross-correlation function reaches the maximum as the distance offset required by the iteration to obtain new pave(r), return to step 402.
Substep B4: further performing phase correction on the data of the range migration error correction curve obtained in the step B3, and correcting a phase error caused by the fluctuating motion by adopting a phase compensation algorithm to obtain a phase error compensation curve;
fig. 5 is a flowchart of distance phase correction processing in the ISAR algorithm, and as shown in fig. 5, the phase compensation algorithm adopts a minimum entropy phase correction algorithm based on a variable step gradient descent method to obtain optimal operation efficiency and result. Establishing an optimization problem model for solving the compensation phase to ensure that the image focusing effect is the best, wherein an objective function is an entropy value of an image and is expressed as follows:
where n is the azimuth sample point, g (r, n) is the phase compensated image, which can be expressed as:
wherein f (r, m) represents data of the azimuth frequency domain after the distance alignment,indicating the amount of phase compensation for the mth echo.
The degree of focus of an image can be measured by the entropy of the image, and the higher the degree of focus, the smaller the entropy of the image and the sharper the image. Therefore, the optimization criterion is that the value of entropy reaches the minimum, under the condition of giving an initial value of a compensation phase, a numerical value iterative algorithm is utilized to gradually approximate to finally solve a target function to obtainIn addition, the experiment searches for the global optimal solution by a gradient descent method with variable step lengthThe value of (2) further improves the iterative computation speed, and simultaneously ensures the final image quality.
The minimum entropy phase correction algorithm based on the variable step gradient descent method comprises the following specific implementation steps:
step S501: initializing a compensation phase to be zero;
step S502: calculating a gradient vector grad of the entropy function at the current position;
step S503: determining magnitude of negative direction change along gradientNamely, a Step length parameter Step needs to be determined;
the specific principle of the method is that the step length is halved each time to make the algorithm converge towards the minimum entropy quickly, and when the convergence is close to the minimum entropy, the optimization result is obtained by continuously fine-tuning the step length.
Step S504: will obtainAndand is added up, as a new acquisition value,and compensating the original image by using the new phase acquisition value to obtain a new image, and circulating to the Step 502 until Step is equal to 0.
Substep B5: and B4, performing azimuth focusing processing on the data of the phase error compensation curve obtained in the step B4 to finally obtain a focused SAR image.
And the azimuth focusing is realized by adopting Fourier transform or inverse Fourier transform.
In the algorithm execution process, the final output result is an SAR image with good focusing effect and a phase error compensation curve, and the SAR image and the phase error compensation curve have a one-to-one correspondence relationship.
And C: fitting a sine function to the phase error compensation curve obtained in the step B4 in combination with the first-order fitting term in the range migration correction curve obtained in the step B3 to obtain a phase error compensation fitting curve in the form of a sine function;
step D: c, correspondingly analyzing the phase error compensation fitting curve in the form of the sine function obtained in the step C in combination with the phase error expression in the sea wave fluctuating motion model established in the step A to obtain the motion amplitude and the motion period of the fluctuating motion target;
the phase error compensation curve represents the phase of each pulse corresponding to the compensationAnd the pulse number represents the azimuth sampling point, i.e. equivalent to the azimuth time:
the phase correction of the defocused target is performed to obtain an SAR image with a good focusing effect, and it can be considered that the phase error is completely compensated, that is, the compensation phase in the obtained phase error compensation curve is the error phase obtained by the analysis in the first step:
the period T of the function can be obtained by sinusoidal fitting of the compensation phase0And amplitude A0Then, the period T and the amplitude a of the rolling motion can be obtained by equation (14) as:
step E: and deducing parameters such as wave height, wave period and the like of the sea waves according to the motion amplitude and the motion period of the fluctuating motion target.
Specifically, the derivation process of obtaining the sea wave parameters according to the motion parameters of the undulating motion target is as follows:
assuming a target mass m, buoyancy and gravity are balanced at rest
Wave height of sea wave is h0Then the buoyancy force applied to the object is increased to:
Fs＝ρgS(h+h0) (17)
the object is acted by upward buoyancy to cause the fluctuation, the movement from the rest position to the highest point is a quarter of a cycle, and the height which can be reached is the amplitude of the fluctuation of the object.
Where A is the amplitude of the heave motion of the object and a is the acceleration.
It should be noted that the above formula (19) is a rough approximation, and the increase of the buoyancy of the object caused by the actual sea wave should be smaller than the above calculation, so a should also be smaller than the result in the above formula, and the ratio is about one fourth to one half, that is:
it can be seen that A is inversely proportional to h, and h reflects the draught of the object and is directly proportional to the weight of the object itself, so that the fluctuation range of the object is inversely proportional to the weight of the object itself, which is consistent with the common sense.
It is to be noted that, in the attached drawings or in the description, the implementation modes not shown or described are all the modes known by the ordinary skilled person in the field of technology, and are not described in detail. Furthermore, the above definition of the model for establishing the wave motion is not limited to the various specific expressions and specific observation models mentioned in the embodiments, and those skilled in the art can easily modify or replace them, for example:
(1) the SAR image can be acquired by using other forms except a satellite, such as a communication form carrying radar transmitting and receiving devices, such as an airplane, a hot air balloon, a microwave shore-based radar and the like, and then a corresponding heave motion model is established;
(2) the wave motion model can be replaced by other functions describing the wave motion besides simple harmonic motion.
By using the technical scheme of the embodiment, the wave parameters of the coast of a certain island in the sea area near the Chinese voyage are obtained. FIG. 6 is an SAR image of the offshore area of a certain island in China voyage under a Terra-SAR sliding bunching high-resolution mode; fig. 7 shows an optical image corresponding to fig. 6 for google earth. The data used for processing comes from a certain island coast of the sea area around the China voyage in the Terra-SAR satellite sliding bunching mode. Table 1 shows the parameters associated with the corresponding image of fig. 6. Combining fig. 6 and 7, it can be seen that the line of mustache defocusing of the island edge is tied to the floating barrel on the bank, because of the up-and-down fluctuating movement of the sea wave, a certain bright line is taken for experiment, a comparison graph before and after refocusing treatment is obtained, as shown in fig. 8, then a phase curve for compensating the line is obtained, as shown in fig. 9, by performing sine fitting on fig. 9 and combining a distance migration correction curve, considering that a certain error still exists, the result of obtaining parameters of the fluctuating movement is obtained as follows:
Ae＝0.14～0.25m
Te＝2.8～3.5s
wherein A iseRepresenting the magnitude of the acquisition, TeThe obtained period is represented, and considering that the floating bucket has weight, the wave height of the fluctuation motion of the floating bucket is about ten times of the obtained amplitude, and the period is equivalent to the period of the fluctuation motion, so the obtained sea wave parameter result is as follows:
As＝1.4～2.5m
Ts＝2.8～3.5s
wherein A issIndicating the wave height, TsRepresenting the wave period.
According to the data, the Chinese vogue belongs to the Bohai sea area, the average wave height of the sea wave is within the range of 0.5-2 m, the period is within 3-5 s, and the obtained wave parameters are within the range, so that the correctness of the method is demonstrated, and meanwhile, the accuracy is certain.
TABLE 1 Terra-SAR Tanshutan data corresponding parameters
| Radar parameter | Numerical value | Radar parameter | Numerical value |
| λ | 0.03112m | PRF | 42300 |
| Bandwidth of | 300MHz | Downward viewing angle | 38.799° |
| Sampling rate | 329.658MHz | Earth's center angle | 3.13° |
| Short distance | 613.981km | Distance resolution | 0.5m |
| Synthetic pore size time | 5.39s | Azimuthal resolution | 0.25m |
Thus, the embodiment of the present invention is described.
In summary, in the embodiment, an algorithm based on inverse transformation and ISAR motion compensation and refocusing is used for accurately obtaining the phase error of the SAR image data containing the moving target, and the obtained phase error is compared with an error expression in the established sea wave motion model, so that the obtained results of parameters such as the wave height and the period of the sea wave can be obtained.
Of course, according to actual needs, the sea wave parameter acquisition method based on sea surface fluctuation target SAR image refocusing provided by the invention further comprises other common algorithms and steps, and is not repeated herein because the method is irrelevant to the innovation of the invention.
The above-mentioned embodiments are intended to illustrate the objects, technical solutions and advantages of the present invention in further detail, and it should be understood that the above-mentioned embodiments are only illustrative of the present invention and are not intended to limit the present invention, and any modifications, equivalents, improvements and the like made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Claims (10)
1. A sea wave parameter obtaining method based on sea surface fluctuation target SAR image refocusing is characterized by comprising the following steps:
acquiring an SAR image of a target area, and applying a refocusing algorithm based on inverse transformation and ISAR motion compensation to data containing a moving target sub-image in the SAR image to obtain a focused SAR image and a phase error compensation curve corresponding to the image;
fitting a sine function to the phase error compensation curve to obtain a phase error compensation fitting curve in the form of the sine function;
correspondingly analyzing the phase error compensation fitting curve in the form of the sine function by combining a phase error expression in a wave fluctuation motion model to obtain the motion amplitude and the cycle of a fluctuation motion target; and
and deducing parameters such as wave height, wave period and the like of the sea waves according to the motion amplitude and the motion period of the fluctuating motion target.
2. A wave parameter acquisition method according to claim 1, wherein the phase error in the wave motion model is expressed as:
wherein,phase error deduced for modeling, AvRepresenting the amplitude of the heave motion, λ representing the pulse wavelength, γ representing the sum of the downward viewing angle and the earth's center angle, ω representing the period of the heave motion, η representing any time instant corresponding to the heave motion,indicates the initial phase angle.
3. A method of acquiring sea wave parameters according to claim 2, wherein the phase error compensation fit curve expression of the sine functional form is:
wherein,to fit the resulting phase error, a represents the amplitude of the heave motion.
4. A method of obtaining sea wave parameters according to claim 3, wherein the phase error compensation fit curve in the form of a sine function is combined with a phase error expression in a model of the heave motion of the sea wave to correspond as follows:
fitting to obtain period T of phase error compensation fitting curve in sine function form0And amplitude A0And the coefficient corresponding to the front of the sine function is equal to the period of the sine function, so that the period T and the amplitude A of the rolling motion are respectively as follows:
T＝T0
5. a sea wave parameter acquisition method according to claim 4, wherein the wave height of the sea wave is acquired according to the amplitude of the undulating moving object by adopting the following formula:
wherein h is0As waves of seaH is the draught of the object, T is the period of the fluctuation motion, A is the amplitude of the fluctuation motion, g is the acceleration of gravity, and a takes a value between 2 and 4.
6. A method for obtaining sea wave parameters according to any one of claims 1-5, wherein the step of performing corresponding analysis on the phase error compensation fitting curve in the form of a sine function in combination with the phase error expression in the sea wave heave motion model further comprises:
and establishing a model of the wave fluctuation motion, and quantitatively analyzing the slope distance error caused by the model, thereby obtaining the phase error expression caused by the slope distance error.
7. A wave parameter acquisition method according to any one of claims 1-5, characterized in that the ISAR motion compensated refocusing algorithm is divided into motion compensation and azimuth focusing, the motion compensation consisting of a distance alignment algorithm and a phase compensation algorithm;
the step of applying a refocusing algorithm based on inverse transformation and ISAR motion compensation to the data containing the moving target sub-image in the SAR image comprises the following steps:
extracting a subimage containing a fluctuating moving target from an SAR image obtained by acquiring information in a target area;
performing azimuth Fourier transform on the sub-image containing the fluctuating moving target to obtain equivalent original data before azimuth compression;
carrying out refocusing processing on the equivalent original data by adopting a distance alignment algorithm to correct a distance migration error caused by fluctuation to obtain a distance migration error correction curve;
performing phase correction on the data of the range migration error correction curve, and correcting a phase error caused by fluctuation motion by adopting a phase compensation algorithm to obtain a phase error compensation curve; and
and performing azimuth focusing processing on the data of the phase error compensation curve to finally obtain a focused SAR image.
8. A wave parameter acquisition method according to claim 7, characterized in that:
the distance alignment algorithm adopts a global minimum entropy distance alignment algorithm, 32-time interpolation is carried out on distance direction data to obtain distance alignment accuracy which is 32 times better than distance resolution, and distance migration errors caused by rolling motion are completely corrected; and/or
The phase compensation algorithm adopts a minimum entropy phase correction algorithm based on a variable step gradient descent method to obtain the optimal operation efficiency and a phase correction result.
9. A method of obtaining sea wave parameters according to claim 7, wherein said azimuthal focusing is performed using Fourier transform or inverse Fourier transform.
10. The ocean wave parameter acquisition method according to claim 7, wherein when the phase error compensation curve is fitted, fitting of a sine function is performed by combining the phase error compensation curve with a first-order fitting term in the distance migration correction curve to obtain the phase error compensation fitting curve in the form of the sine function.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710120632.5A CN106872979B (en) | 2017-03-02 | 2017-03-02 | Sea wave parameter acquisition method based on sea surface fluctuation moving target SAR image refocusing |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710120632.5A CN106872979B (en) | 2017-03-02 | 2017-03-02 | Sea wave parameter acquisition method based on sea surface fluctuation moving target SAR image refocusing |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106872979A true CN106872979A (en) | 2017-06-20 |
| CN106872979B CN106872979B (en) | 2019-12-20 |
Family
ID=59169273
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710120632.5A Active CN106872979B (en) | 2017-03-02 | 2017-03-02 | Sea wave parameter acquisition method based on sea surface fluctuation moving target SAR image refocusing |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106872979B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108460419A (en) * | 2018-03-07 | 2018-08-28 | 中国科学院武汉岩土力学研究所 | The fracture parameters extracting method of drilling optical imagery and the fusion of radar imagery information |
| CN111983612A (en) * | 2020-08-26 | 2020-11-24 | 中国科学院空天信息创新研究院 | SAR sliding bunching mode azimuth declivity method |
| CN113589281A (en) * | 2020-04-30 | 2021-11-02 | 北京理工大学重庆创新中心 | GEO SAR ship target imaging method based on micro-Doppler analysis |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3783058B2 (en) * | 2003-07-19 | 2006-06-07 | 独立行政法人情報通信研究機構 | Method and system for inverse estimation of wave direction spectrum from radar image |
| CN102034233A (en) * | 2010-10-21 | 2011-04-27 | 苏州科技学院 | Method for detecting SAR (stop and reveres) image wave group parameters based on contourlet conversion |
| CN102799770A (en) * | 2012-06-29 | 2012-11-28 | 哈尔滨工程大学 | Method for modeling sea wave significant wave height inversion model based on particle swarm optimization (PSO) self-adaptive piecewise linear fitting |
| CN103383455A (en) * | 2013-07-03 | 2013-11-06 | 南京信息工程大学 | Ocean wave parameter extraction method based on shape-from-shading (SFS) technology |
| CN104331588A (en) * | 2014-03-06 | 2015-02-04 | 浙江海洋学院 | A joint inversion device of sea surface wind field and sea waves |
-
2017
- 2017-03-02 CN CN201710120632.5A patent/CN106872979B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3783058B2 (en) * | 2003-07-19 | 2006-06-07 | 独立行政法人情報通信研究機構 | Method and system for inverse estimation of wave direction spectrum from radar image |
| CN102034233A (en) * | 2010-10-21 | 2011-04-27 | 苏州科技学院 | Method for detecting SAR (stop and reveres) image wave group parameters based on contourlet conversion |
| CN102799770A (en) * | 2012-06-29 | 2012-11-28 | 哈尔滨工程大学 | Method for modeling sea wave significant wave height inversion model based on particle swarm optimization (PSO) self-adaptive piecewise linear fitting |
| CN103383455A (en) * | 2013-07-03 | 2013-11-06 | 南京信息工程大学 | Ocean wave parameter extraction method based on shape-from-shading (SFS) technology |
| CN104331588A (en) * | 2014-03-06 | 2015-02-04 | 浙江海洋学院 | A joint inversion device of sea surface wind field and sea waves |
Non-Patent Citations (3)
| Title |
|---|
| DAIYIN ZHU，ET AL: "ISAR Range Alignment by Minimizing the Entropy of the Average Range Profile", 《IEEE》 * |
| J. WANG,ET AL: "Minimum-entropy phase adjustment for ISAR", 《IEE PROC.-RADAR SONAR NAVIG.》 * |
| XUEJIAO WEN,ET AL: "Focusing and Parameter Estimating of Fluctuating Target in High Resolution Spaceborne SAR", 《2016 CIE ITERNATIONAL CONFERENCE ON RADAR》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108460419A (en) * | 2018-03-07 | 2018-08-28 | 中国科学院武汉岩土力学研究所 | The fracture parameters extracting method of drilling optical imagery and the fusion of radar imagery information |
| CN108460419B (en) * | 2018-03-07 | 2020-05-08 | 中国科学院武汉岩土力学研究所 | Fracture parameter extraction method for fusion of drilling optical image and radar imaging information |
| CN113589281A (en) * | 2020-04-30 | 2021-11-02 | 北京理工大学重庆创新中心 | GEO SAR ship target imaging method based on micro-Doppler analysis |
| CN111983612A (en) * | 2020-08-26 | 2020-11-24 | 中国科学院空天信息创新研究院 | SAR sliding bunching mode azimuth declivity method |
| CN111983612B (en) * | 2020-08-26 | 2022-04-15 | 中国科学院空天信息创新研究院 | SAR sliding bunching mode azimuth declivity method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106872979B (en) | 2019-12-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111583214B (en) | Sea surface wind speed inversion method based on RBF neural network and based on marine radar image | |
| CN110609287B (en) | Double-frequency radar scatterometer and method for simultaneously measuring sea surface wind field and flow field | |
| Romeiser et al. | First analysis of TerraSAR-X along-track InSAR-derived current fields | |
| CN103926572B (en) | A kind of clutter suppression method of battle array radar self adaptation subspace, airborne anon-normal side | |
| CN103969643B (en) | One carries out X-band pathfinder inverting ocean wave parameter method based on novel wave dispersion relation band filter | |
| CN108007401A (en) | A kind of river and lake storehouse bank deformation detecting device and method based on boat-carrying InSAR platforms | |
| CN103941257B (en) | A kind of method of the pathfinder image inverting wind direction of ocean surface based on wave number energy spectrum | |
| CN104613944A (en) | Distributed water depth prediction method based on GWR (geographically weighted regression) and BP (back propagation) neural network | |
| CN109856635A (en) | CSAR ground moving target refocusing imaging method | |
| CN104035095A (en) | Low-altitude wind shear wind speed estimating method based on space-time optimal processor | |
| CN106872979B (en) | Sea wave parameter acquisition method based on sea surface fluctuation moving target SAR image refocusing | |
| CN104950307B (en) | Accurate locating method for onboard tri-channel SAR-GMTI (Synthetic Aperture Radar-Ground Moving Target Indication) | |
| CN105487074B (en) | A kind of double-base synthetic aperture radar numerical distance Doppler imaging method | |
| CN105891832A (en) | Sea surface wind speed inversion method based on cross polarization model and CMOD5N | |
| CN104020471A (en) | Partitioning processing-based SAR real-time imaging method and system thereof | |
| CN102565772B (en) | Marine dynamic information extraction method on basis of SAR (Synthetic Aperture Radar) sub-aperture sequence images | |
| CN108710111A (en) | A kind of two-dimentional space-variant bearing calibration of airborne biradical Forward-looking SAR orientation phase | |
| CN106768179A (en) | The measuring method of the tidal level based on continuous operation GNSS stations signal-to-noise ratio data | |
| CN105180852A (en) | GB-SAR deformation monitoring method based on triple stepping | |
| CN108646247A (en) | Inverse synthetic aperture radar imaging method based on Gamma process linear regression | |
| CN107356923A (en) | A kind of ISAR based on sub-aperture division is imaged envelope alignment method | |
| Li et al. | Mapping the bathymetry of shallow coastal water using single-frame fine-resolution optical remote sensing imagery | |
| CN110554377B (en) | Single-channel SAR two-dimensional flow field inversion method and system based on Doppler center offset | |
| CN114002720A (en) | Ship positioning and meteorological data inversion method based on ocean tide load | |
| CN108983235A (en) | A kind of cross polarization SAR Wind-field Retrieval method based on swath |
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 | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20211020 Address after: 250000 floor 24, block B, No. 188, Xiaoqinghe North Road, Tianqiao District, Jinan City, Shandong Province Patentee after: Zhongke satellite (Shandong) Technology Group Co.,Ltd. Address before: 100089 No.19, Beisihuan West Road, Haidian District, Beijing Patentee before: Aerospace Information Research Institute,Chinese Academy of Sciences Effective date of registration: 20211020 Address after: 100089 No.19, Beisihuan West Road, Haidian District, Beijing Patentee after: Aerospace Information Research Institute,Chinese Academy of Sciences Address before: 100190 No. 19 West North Fourth Ring Road, Haidian District, Beijing Patentee before: Institute of Electronics, Chinese Academy of Sciences |