CN111896954A - Corner reflector coordinate positioning method for shipborne SAR image - Google Patents
Corner reflector coordinate positioning method for shipborne SAR image Download PDFInfo
- Publication number
- CN111896954A CN111896954A CN202010781923.0A CN202010781923A CN111896954A CN 111896954 A CN111896954 A CN 111896954A CN 202010781923 A CN202010781923 A CN 202010781923A CN 111896954 A CN111896954 A CN 111896954A
- Authority
- CN
- China
- Prior art keywords
- corner reflector
- sar image
- shipborne
- triangular conical
- coordinate positioning
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 60
- 238000012937 correction Methods 0.000 claims abstract description 14
- 230000005855 radiation Effects 0.000 claims abstract description 9
- 238000003384 imaging method Methods 0.000 claims description 38
- 230000010287 polarization Effects 0.000 claims description 8
- 238000005316 response function Methods 0.000 claims description 6
- 239000000463 material Substances 0.000 claims description 4
- 238000012216 screening Methods 0.000 claims description 4
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 230000003287 optical effect Effects 0.000 claims description 3
- 239000011148 porous material Substances 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 230000009467 reduction Effects 0.000 claims description 3
- 238000003786 synthesis reaction Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 5
- 238000012795 verification Methods 0.000 abstract description 5
- 230000006872 improvement Effects 0.000 description 6
- 238000004088 simulation Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 3
- 230000000052 comparative effect Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- PXFBZOLANLWPMH-UHFFFAOYSA-N 16-Epiaffinine Natural products C1C(C2=CC=CC=C2N2)=C2C(=O)CC2C(=CC)CN(C)C1C2CO PXFBZOLANLWPMH-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000005305 interferometry Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/02—Affine transformations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/80—Geometric correction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/70—Determining position or orientation of objects or cameras
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20016—Hierarchical, coarse-to-fine, multiscale or multiresolution image processing; Pyramid transform
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses a corner reflector coordinate positioning method of a shipborne SAR image, which comprises the following steps: selecting a corner reflector, calculating SAR image characteristics, calculating spatial characteristics, correcting characteristic points and deducing a central point of the corner reflector; based on the scale space, the invention considers the radiation characteristic and the space characteristic of the corner reflector in the SAR image, and after verification, the invention has equivalent precision and identical result for the average relative error and the root mean square error of the corner reflector positioning result of the high-resolution shipborne SAR image and the corner reflector coordinate positioning result after geometric correction based on the control point in the distance direction and the azimuth direction, thereby proving that the invention does not depend on the geographic coordinate information of the corner reflector, the corner reflector coordinate positioning precision can reach the sub-pixel level, and the azimuth coordinate positioning precision is superior to the distance direction coordinate because the focusing effect of the corner reflector in the azimuth direction is superior to the distance direction coordinate.
Description
Technical Field
The invention relates to the technical field of coordinate positioning, in particular to a corner reflector coordinate positioning method of a shipborne SAR image.
Background
The corner reflector can strongly reflect radar beams due to the particularity of the design of the corner reflector, particularly, echoes of the corner reflector are overlapped for many times in the imaging process of the synthetic aperture radar, a high-brightness area is presented in an SAR image, and then the scattering characteristic can be widely applied to external calibration of a radar system, indirect representation of ground micro-deformation, image registration basis, image restoration, CRINSAR technology and the like;
corner reflector positioning, which refers to a process of obtaining accurate SAR image coordinates corresponding to a centroid of a corner reflector, is fundamental work of corner reflector related application, and about the research of the corner reflector on the coordinate positioning in an SAR image, the method mainly focuses on satellite-borne SAR images, because the resolution of the satellite-borne SAR images is low, the images of the corner reflectors generally show highlight points or areas and have no obvious cross-hair characteristics, the commonly used methods for corner reflector positioning are mostly based on the radiation intensity characteristics of the corner reflectors, such as a template matching method based on intensity characteristics, a peak zero filling method, a method of a spectral correlation diagram and the like, other types of methods also comprise a method based on change detection, a dual threshold method based on PSInSAR, a comprehensive index multiple method based on intensity and coherence coefficients, a detection method based on prior knowledge and the like, and most of the methods can only achieve positioning accuracy at a level, such as template matching, spectral correlation, etc.; a small part of methods which can reach the accuracy of sub-pixel levels, such as a peak zero filling method, are influenced by factors such as background and noise, the positioning effect is difficult to achieve expectation, the characteristics of a cross wire formed by imaging a corner reflector in a high-resolution shipborne SAR image are obvious, the image coordinate positioning accuracy of the corner reflector has the potential of further improvement, and the method has important significance for obtaining the image coordinates of the sub-pixel levels with high accuracy by taking the corner reflector as a ground control point in the subsequent operation of using the shipborne SAR for interferometry and DEM extraction;
the identification of the corner reflector can be combined with the characteristics of gray scale, shape and the like expressed in the SAR image, and for a high-resolution shipborne SAR image, the characteristics of the corner reflector are rich, so that the identification strategies and the coordinate positioning method are more, according to the existing research, in the coordinate positioning method of the common corner reflector in the SAR image, for example, the method based on the radiation intensity characteristic is easily influenced by noise, the method based on the characteristic parameter can not reach the ideal precision, and the method based on the radiation calibration and the geometric correction needs to carry out field measurement on the corner reflector to obtain the geographic coordinate, which wastes time and labor, so the invention provides the corner reflector coordinate positioning method of the shipborne SAR image to solve the problems in the prior art.
Disclosure of Invention
In view of the above problems, an object of the present invention is to provide a corner reflector coordinate positioning method for a shipborne SAR image, which does not depend on geographic coordinate information of a corner reflector, and the corner reflector coordinate positioning accuracy can reach a sub-pixel level.
In order to realize the purpose of the invention, the invention is realized by the following technical scheme: a corner reflector coordinate positioning method of a shipborne SAR image comprises the following steps:
the method comprises the following steps: selective angle reflector
The method is characterized in that the whole RCS curve corresponding to the triangular conical corner reflector changes slowly, the RCS value reduction speed is slowest within the range of deviating from the maximum incidence direction by 40 degrees, and the triangular conical corner reflector has large tolerance to radar incident wave deviation, so that the triangular conical corner reflector is used as a common corner reflector, SAR is used for recording the amplitude and phase information of backscattered electromagnetic waves of the triangular conical corner reflector in an image form, an intensity map of the triangular conical corner reflector is obtained after the amplitude value is calibrated and squared, and the RCS of the triangular conical corner reflector is obtained on the intensity map, wherein the unit is m2;
Step two: SAR image feature computation
Regarding the backscattering signal of the triangular conical corner reflector to the SAR as an impulse response function, in a focused SAR image, the envelope surface of the expressed signal is as follows:
wherein h (x, t) represents the impulse response function of the t-x plane, LSFor synthesis of the pore time, τρIs the pulse width, t is the ramp time, x is the azimuth angle, ρrFor distance resolution, paFor azimuth resolution, the triangular conical corner reflector is obtained by the formula (1), and in an ideal state, the intensity is expressed in a sinc function form in the azimuth direction and the distance direction, so that the triangular conical corner reflector is expressed in a cross-hair shape in a high-resolution SAR image;
step three: spatial feature computation
Expressing the Gaussian and Gaussian difference scale space obtained by the triangular conical corner reflector imaging area as follows:
in the formula, E (t, x) represents E (h (t, x)), G (t, x, σ) represents a gaussian convolution kernel whose spatial scale is σ, L (t, x, σ), D (t, x, σ) represent a gaussian scale space and a gaussian difference scale space, respectively, and k is a ratio between spatial scales of two layers of gaussian scale spaces, as can be seen from formula (2), D (t, x, σ), spatial scale σ, and image resolution ρa,ρrCorrelation;
step four: feature point correction
Establishing a Gaussian pyramid and a Gaussian difference pyramid by the formula (2), detecting extreme points, obtaining a characteristic point X in a scale space, and correcting the X:
step five: derivative corner reflector center point
In actual processing, the number of feature points detected in the triangular pyramidal retroreflector imaging region is large, and the centroid of the retroreflector is detected as a feature point at a certain scale σ based on imaging feature formula (1) of the retroreflector, and therefore the centroid of the retroreflector corresponds to the imaging point in the feature point set P ═ X1,X2,...,XnAnd (4) screening, wherein according to the existing corner reflector coordinate positioning method: the template matching method and the peak zero-filling method firstly acquire an approximate coordinate X' of the corner reflector, and when the formula (4) is satisfied, the characteristic point is considered as the center point of the corner reflector:
|Xi-X′|<R (4)
wherein R is a distance threshold, i is 1,2, …, n.
The further improvement lies in that: in the first step, RCS is a radar scattering cross section and represents the echo intensity and the backscattering energy of the target generated after the radar transmits the electromagnetic waves, and RCS is a function consisting of frequency, incident wave polarization, receiving antenna polarization, target shape, structure and material characteristic factors and has no direct relation with the area of the target displayed on radar beams.
The further improvement lies in that: in the second step, in the focused SAR image, the echo energy distribution of the triangular conical corner reflector is similar to the point spread function PSF of the optical system.
The further improvement lies in that: in the third step, the imaging of the triangular conical corner reflector not only has obvious radiation characteristics, but also has spatial characteristics, and under the condition of different image resolutions, the imaging blur degree of the triangular conical corner reflector is different and is similar to the scale space of the image.
The further improvement lies in that: in the fourth step of the method, the first step of the method,and the corrected result of X comprises the precise coordinates of the characteristic points and the scale parameters.
The further improvement lies in that: in the fifth step, the size of the sigma is related to the size of the triangular conical corner reflector and the resolution of the shipborne SAR image.
The invention has the beneficial effects that: based on a scale space, the radiation characteristic and the space characteristic of a corner reflector in an SAR image are considered at the same time, and verification proves that the positioning accuracy of the corner reflector of a high-resolution shipborne SAR image is equivalent to the accuracy of the average relative error and the root mean square error of the coordinate positioning result of the corner reflector after geometric correction based on a control point in the distance direction and the azimuth direction, the results are identical, and the result proves that the positioning accuracy of the corner reflector can reach the sub-pixel level without depending on the geographic coordinate information of the corner reflector, and the focusing effect of the corner reflector in the azimuth direction is superior to that of the distance direction, so the positioning accuracy of the azimuth coordinate is superior to that of the distance direction coordinate, and compared with the distance direction, the strong reflection pixel in the azimuth direction is more, the details are more abundant, and the positioning accuracy is closer to the ideal imaging situation of the corner reflector.
Drawings
FIG. 1 is a RCS graph of a triangular pyramidal corner reflector of the present invention;
FIG. 2 is a three-dimensional perspective of the sinc function of the present invention;
FIG. 3 is a diagram of the triangular pyramidal corner reflector imaging of the present invention;
FIG. 4 is an imaging feature map of the corner reflector at a resolution of 241 × 241 images of 0.25m in accordance with the present invention;
FIG. 5 is an imaging feature map of the corner reflector at 1m, 61X 61 image resolution according to the present invention;
FIG. 6 is an imaging feature map of the corner reflector at 2m, 31X 31 image resolution according to the present invention;
FIG. 7 is an image feature map of the corner reflector at 4m, 16X 16 image resolution in accordance with the present invention;
FIG. 8 is a graph of corner reflector coordinate positioning error of the present invention;
FIG. 9 is a spatial scale view of a corner reflector centering feature of the present invention;
FIG. 10 is a graph of the coordinate difference between the present invention and the verification data.
Detailed Description
In order to further understand the present invention, the following detailed description will be made with reference to the following examples, which are only used for explaining the present invention and are not to be construed as limiting the scope of the present invention.
As shown in fig. 1,2 and 3, the present embodiment provides a method for positioning coordinates of a corner reflector of a ship-borne SAR image, including the following steps:
the method comprises the following steps: selective angle reflector
As the whole RCS curve corresponding to the triangular conical corner reflector changes slowly and the RCS value reduction speed is slowest within the range of deviating from the maximum incidence direction by 40 degrees, as shown in figure 1, the triangular conical corner reflector has large tolerance to the radar incident wave deviation, so the triangular conical corner reflector is taken as a commonly used corner reflector, the RCS is a radar scattering sectional area and represents the echo intensity and the backscattering energy of a target generated after the radar transmits an electromagnetic wave, and the RCS is a function consisting of frequency, incident wave polarization, receiving antenna polarization, target shape, structure and material characteristic factors and is a function consisting of frequency, incident wave polarization, receiving antenna polarization, target shape, structure and material characteristic factorsThere is no direct relationship to the area that the target exhibits on the radar beam; then, utilizing SAR to record the amplitude and phase information of the backscattering electromagnetic wave of the triangular conical corner reflector in the form of image, calibrating and squaring the amplitude value to obtain the intensity map of the triangular conical corner reflector, and obtaining the RCS of the triangular conical corner reflector on the intensity map, wherein the unit is m2;
Step two: SAR image feature computation
Regarding a backscattering signal of the triangular conical corner reflector to the SAR as an impulse response function, in a focused SAR image, the echo energy distribution of the triangular conical corner reflector is similar to a point spread function PSF of an optical system, and the envelope surface of the expressed signal is as follows:
wherein h (x, t) represents the impulse response function of the t-x plane, LSFor synthesis of the pore time, τρIs the pulse width, t is the ramp time, x is the azimuth angle, ρrFor distance resolution, paFor the azimuth resolution, the triangular conical corner reflector is obtained by the formula (1), and in an ideal state, the intensity is expressed in a sinc function form in the azimuth direction and the distance direction, as shown in fig. 2, so that the triangular conical corner reflector is expressed in a cross-hair shape in a high-resolution SAR image; as shown in FIG. 3;
step three: spatial feature computation
The imaging of the triangular conical corner reflector not only has obvious radiation characteristics, but also has spatial characteristics, under the condition of different image resolutions, the imaging blur degree of the triangular conical corner reflector is different and is similar to the scale space of an image, and the Gaussian difference scale space obtained by the imaging area of the triangular conical corner reflector is expressed as follows:
in the formula, E (t, x) represents E (h (t, x)), G (t, x, σ) represents a Gaussian convolution kernel having a spatial scale σ, and L (t, x, σ),D (t, x, sigma) represents a Gaussian scale space and a Gaussian difference scale space respectively, k is the ratio of the spatial scales of two layers of Gaussian scale spaces, and D (t, x, sigma) and the spatial scale sigma as well as the image resolution rho can be known from the formula (2)a,ρrCorrelation;
step four: feature point correction
Establishing a Gaussian pyramid and a Gaussian difference pyramid by the formula (2), detecting extreme points, obtaining a characteristic point X in a scale space, and correcting the X:
where X ═ t, X, σ, X is the imaging target point,the result after X correction comprises the precise coordinates and the scale parameters of the characteristic points;
step five: derivative corner reflector center point
In actual processing, the number of feature points detected in the triangular conical corner reflector imaging area is large, according to the imaging feature formula (1) of the corner reflector, the centroid of the corner reflector is detected as a feature point in a certain scale sigma, and the size of sigma is related to the size of the triangular conical corner reflector and the resolution of the ship-borne SAR image, so that the centroid of the corner reflector corresponds to the imaging point in a feature point set P { X ═ X { (X) }1,X2,...,XnAnd (4) screening, wherein according to the existing corner reflector coordinate positioning method: the template matching method and the peak zero-filling method firstly acquire an approximate coordinate X' of the corner reflector, and when the formula (4) is satisfied, the characteristic point is considered as the center point of the corner reflector:
|Xi-X′|<R (4)
wherein R is a distance threshold, i is 1,2, …, n.
Simulation test:
in order to verify the feasibility of corner reflector coordinate positioning based on a scale space, corner reflector imaging with different resolutions is simulated, and coordinate positioning is performed on a corner reflector central point based on the scale space.
Assuming that the resolutions of the azimuth direction and the range direction of the SAR image are the same, at the moment, only the positioning accuracy of the range direction coordinate needs to be considered, the simulation parameters comprise range direction bandwidth, range direction sampling frequency, range direction resolution and the like, the range direction bandwidth is 600MHz, under the condition of different range direction resolutions, the imaging simulation result of the corner reflector is shown in figures 4-7, the higher the image resolution is, the larger the number of pixels occupied by the bright area of the corner reflector is, and the more the details are.
And establishing a Gaussian and Gaussian difference pyramid for the corner reflector region images with different resolutions by using a siftDemo V4 program, detecting the characteristic points, screening out the characteristic points close to the central point of the image, and comparing the characteristic points with the coordinates of the central point (the central point of the image) of the real corner reflector. Under the condition of different distance-direction resolutions, the coordinate positioning error of the central point of the corner reflector and the corresponding spatial scale are respectively shown in fig. 8 and fig. 9. In fig. 8, under the condition of different resolutions (0 to 6m), the positioning accuracy of the central point of the corner reflector is extremely high (within 0.02 pixel error), and in fig. 9, the image resolution and the spatial scale of the feature point corresponding to the center of the corner reflector exhibit a strong inverse proportion relationship, which indicates that the imaging area of the corner reflector has a strong spatial scale characteristic, and the higher the image resolution is, the larger the spatial scale of the imaging area of the corner reflector is. Fig. 8 and 9 demonstrate the feasibility of corner reflector coordinate positioning based on scale space.
And (3) testing:
further verification is carried out by adopting ship-borne SAR measured data, the main parameters of the ship-borne SAR image are shown in a table 1, and a double-antenna rail-crossing imaging mode is adopted. The corner reflector has good imaging condition in the shipborne SAR image, and is in a more obvious cross-hair shape.
TABLE 1 Main parameters of shipborne SAR images
The control net is arranged at the central position of the corner reflector, and the accurate geographic coordinates of the central position of the corner reflector are obtained by using the traditional measuring means (total station forward intersection measurement), wherein errors in the measurement of the point coordinates are less than 6 mm.
And interpolating the imaging area of the corner reflector to 0.01 pixel by adopting a peak zero filling method to obtain image coordinates of the center points of all corner reflectors, selecting the corner reflectors CR7 and CR9 as geometric correction control points, and performing geometric correction by combining a Doppler imaging model and an affine transformation optimization model to obtain a high-precision image coordinate result of the corner reflector as shown in table 2. The results of table 2 can be made as comparative data for the corner reflector coordinate positioning method herein.
TABLE 2 geometric correction of relief reflector image coordinates
The specific coordinate positioning result of the corner reflector based on the scale space is shown in table 3, and by analyzing table 3, we know that the feature point scale values corresponding to all corner reflectors are close, and the sizes of all the triangular pyramid type corner reflectors arranged are the same, which is consistent with the simulation result of fig. 9, that is, for the corner reflectors with the same size, the feature points are detected on the same scale.
TABLE 3 corner Reflector coordinate positioning results based on scale space
The difference between the coordinate positioning result of the corner reflector based on the spatial scale and the comparison data obtained by combining the tables 2 and 3 is shown in fig. 10, and the average relative error and the root mean square error result of the coordinate positioning result and the comparison data obtained by the method are shown in table 4.
Table 4 error of the method and comparative data herein
Analyzing the above results in conclusion:
(1) the method has the advantages that the average relative error and the root mean square error of the corner reflector positioning result of the high-resolution shipborne SAR image and the corner reflector coordinate positioning result after geometric correction based on the control point in the distance direction and the azimuth direction are less than 0.31 pixel, and the accuracy is equivalent, which shows that the positioning of the corner reflector in the shipborne SAR image in the method can reach sub-pixel level accuracy under the condition of no ground control point;
(2) the results in table 3 indicate that the positioning accuracy of the azimuth coordinate is better than that of the distance coordinate because the focusing effect of the corner reflector in the azimuth direction is better than that of the distance direction, and compared with the distance direction, the strong reflection pixels in the azimuth direction are more, the details are richer and are closer to the ideal imaging situation of the corner reflector.
The corner reflector coordinate positioning method of the shipborne SAR image is based on a scale space, simultaneously considers the radiation characteristic and the space characteristic of a corner reflector in the SAR image, and after verification, the invention proves that the average relative error and the root mean square error of the corner reflector positioning result of the high-resolution shipborne SAR image and the corner reflector coordinate positioning result after geometric correction based on a control point in the distance direction and the azimuth direction are less than 0.31 pixel, the precision is equivalent, the result is identical, the result proves that the method does not depend on the geographic coordinate information of the corner reflector, the corner reflector coordinate positioning precision can reach the sub-pixel level, and the focusing effect of the corner reflector in the azimuth direction is superior to that in the distance direction, therefore, the positioning accuracy of the azimuth coordinate is superior to that of the distance coordinate, and compared with the distance coordinate, the azimuth upward strong reflection image elements are more, the details are richer, and the azimuth upward strong reflection image elements are closer to the ideal imaging situation of the corner reflector.
The foregoing illustrates and describes the principles, general features, and advantages of the present invention. It will be understood by those skilled in the art that the present invention is not limited to the embodiments described above, which are described in the specification and illustrated only to illustrate the principle of the present invention, but that various changes and modifications may be made therein without departing from the spirit and scope of the present invention, which fall within the scope of the invention as claimed. The scope of the invention is defined by the appended claims and equivalents thereof.
Claims (6)
1. A corner reflector coordinate positioning method of a shipborne SAR image is characterized by comprising the following steps:
the method comprises the following steps: selective angle reflector
The method is characterized in that the whole RCS curve corresponding to the triangular conical corner reflector changes slowly, the RCS value reduction speed is slowest within the range of deviating from the maximum incidence direction by 40 degrees, and the triangular conical corner reflector has large tolerance to radar incident wave deviation, so that the triangular conical corner reflector is used as a common corner reflector, SAR is used for recording the amplitude and phase information of backscattered electromagnetic waves of the triangular conical corner reflector in an image form, an intensity map of the triangular conical corner reflector is obtained after the amplitude value is calibrated and squared, and the RCS of the triangular conical corner reflector is obtained on the intensity map, wherein the unit is m2;
Step two: SAR image feature computation
Regarding the backscattering signal of the triangular conical corner reflector to the SAR as an impulse response function, in a focused SAR image, the envelope surface of the expressed signal is as follows:
wherein h (x, t) represents the impulse response function of the t-x plane, LSFor synthesis of the pore time, τρIs the pulse width, t is the ramp time, x is the azimuth angle, ρrFor distance resolution, paFor azimuth resolution, the triangular conical corner reflector is obtained by the formula (1), and in an ideal state, the intensity is expressed in a sinc function form in the azimuth direction and the distance direction, so that the triangular conical corner reflector is expressed in a cross-hair shape in a high-resolution SAR image;
step three: spatial feature computation
Expressing the Gaussian and Gaussian difference scale space obtained by the triangular conical corner reflector imaging area as follows:
in the formula, E (t, x) represents E (h (t, x)), G (t, x, σ) represents a gaussian convolution kernel whose spatial scale is σ, L (t, x, σ), D (t, x, σ) represent a gaussian scale space and a gaussian difference scale space, respectively, and k is a ratio between spatial scales of two layers of gaussian scale spaces, as can be seen from formula (2), D (t, x, σ), spatial scale σ, and image resolution ρa,ρrCorrelation;
step four: feature point correction
Establishing a Gaussian pyramid and a Gaussian difference pyramid by the formula (2), detecting extreme points, obtaining a characteristic point X in a scale space, and correcting the X:
step five: derivative corner reflector center point
In actual processing, the number of feature points detected in the triangular pyramidal retroreflector imaging region is large, and the centroid of the retroreflector is detected as a feature point at a certain scale σ based on imaging feature formula (1) of the retroreflector, and therefore the centroid of the retroreflector corresponds to the imaging point in the feature point set P ═ X1,X2,...,XnAnd (4) screening, wherein according to the existing corner reflector coordinate positioning method: the template matching method and the peak zero-filling method firstly acquire an approximate coordinate X' of the corner reflector, and when the formula (4) is satisfied, the characteristic point is considered as the center point of the corner reflector:
|Xi-X′|<R (4)
wherein R is a distance threshold, i is 1,2, …, n.
2. The corner reflector coordinate positioning method for the shipborne SAR image according to claim 1, characterized in that: in the first step, RCS is a radar scattering cross section and represents the echo intensity and the backscattering energy of the target generated after the radar transmits the electromagnetic waves, and RCS is a function consisting of frequency, incident wave polarization, receiving antenna polarization, target shape, structure and material characteristic factors and has no direct relation with the area of the target displayed on radar beams.
3. The corner reflector coordinate positioning method for the shipborne SAR image according to claim 1, characterized in that: in the second step, in the focused SAR image, the echo energy distribution of the triangular conical corner reflector is similar to the point spread function PSF of the optical system.
4. The corner reflector coordinate positioning method for the shipborne SAR image according to claim 1, characterized in that: in the third step, the imaging of the triangular conical corner reflector not only has obvious radiation characteristics, but also has spatial characteristics, and under the condition of different image resolutions, the imaging blur degree of the triangular conical corner reflector is different and is similar to the scale space of the image.
5. The corner reflector coordinate positioning method for the shipborne SAR image according to claim 1, characterized in that: in the fourth step of the method, the first step of the method,and the corrected result of X comprises the precise coordinates of the characteristic points and the scale parameters.
6. The corner reflector coordinate positioning method for the shipborne SAR image according to claim 1, characterized in that: in the fifth step, the size of the sigma is related to the size of the triangular conical corner reflector and the resolution of the shipborne SAR image.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010781923.0A CN111896954A (en) | 2020-08-06 | 2020-08-06 | Corner reflector coordinate positioning method for shipborne SAR image |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010781923.0A CN111896954A (en) | 2020-08-06 | 2020-08-06 | Corner reflector coordinate positioning method for shipborne SAR image |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111896954A true CN111896954A (en) | 2020-11-06 |
Family
ID=73246520
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010781923.0A Pending CN111896954A (en) | 2020-08-06 | 2020-08-06 | Corner reflector coordinate positioning method for shipborne SAR image |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111896954A (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112529945A (en) * | 2020-11-17 | 2021-03-19 | 西安电子科技大学 | Registration method for multi-view three-dimensional ISAR scattering point set |
CN113640758A (en) * | 2021-08-23 | 2021-11-12 | 中国科学院空天信息创新研究院 | SAR image scaler placement method and system under urban complex environment |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101376687B1 (en) * | 2012-12-06 | 2014-03-20 | 국방과학연구소 | Terrain analysis method based on grid map using quadtree |
CN107367716A (en) * | 2017-07-04 | 2017-11-21 | 武汉大学 | A kind of high-precision satellite-borne SAR geometric calibration method |
US10042048B1 (en) * | 2014-02-20 | 2018-08-07 | National Technology & Engineering Solutions Of Sandia, Llc | Superpixels for improved structure and terrain classification using multiple synthetic aperture radar image products |
CN108508439A (en) * | 2018-05-01 | 2018-09-07 | 南京理工大学 | The method that double carried SARs position target cooperative imaging volume |
CN109188433A (en) * | 2018-08-20 | 2019-01-11 | 南京理工大学 | The method of two-shipper borne SAR image target positioning based on no control point |
CN110363758A (en) * | 2019-07-22 | 2019-10-22 | 中国科学院合肥物质科学研究院 | A kind of Optical remote satellite image quality determines method and system |
US20190324133A1 (en) * | 2018-01-30 | 2019-10-24 | Oculii Corp. | Systems and methods for interpolated virtual aperature radar tracking |
-
2020
- 2020-08-06 CN CN202010781923.0A patent/CN111896954A/en active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101376687B1 (en) * | 2012-12-06 | 2014-03-20 | 국방과학연구소 | Terrain analysis method based on grid map using quadtree |
US10042048B1 (en) * | 2014-02-20 | 2018-08-07 | National Technology & Engineering Solutions Of Sandia, Llc | Superpixels for improved structure and terrain classification using multiple synthetic aperture radar image products |
CN107367716A (en) * | 2017-07-04 | 2017-11-21 | 武汉大学 | A kind of high-precision satellite-borne SAR geometric calibration method |
US20190324133A1 (en) * | 2018-01-30 | 2019-10-24 | Oculii Corp. | Systems and methods for interpolated virtual aperature radar tracking |
CN108508439A (en) * | 2018-05-01 | 2018-09-07 | 南京理工大学 | The method that double carried SARs position target cooperative imaging volume |
CN109188433A (en) * | 2018-08-20 | 2019-01-11 | 南京理工大学 | The method of two-shipper borne SAR image target positioning based on no control point |
CN110363758A (en) * | 2019-07-22 | 2019-10-22 | 中国科学院合肥物质科学研究院 | A kind of Optical remote satellite image quality determines method and system |
Non-Patent Citations (3)
Title |
---|
LOWE等: "Distinctive Image Features from Scale-Invariant Keypoints", 《INTERNATIONAL JOURNAL OF COMPUTER VISION》 * |
张桂芳等: "人工角反射器辐射特性及其像素级精定位参数", 《地震地质》 * |
戴国梦等: "基于尺度空间的角反射器车载SAR 影像坐标定位", 《测绘通报》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112529945A (en) * | 2020-11-17 | 2021-03-19 | 西安电子科技大学 | Registration method for multi-view three-dimensional ISAR scattering point set |
CN112529945B (en) * | 2020-11-17 | 2023-02-21 | 西安电子科技大学 | Multi-view three-dimensional ISAR scattering point set registration method |
CN113640758A (en) * | 2021-08-23 | 2021-11-12 | 中国科学院空天信息创新研究院 | SAR image scaler placement method and system under urban complex environment |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP3364212A1 (en) | A method and an apparatus for computer-assisted processing of sar raw data | |
CN112305510B (en) | DEM matching-based synthetic aperture radar image geometric calibration method | |
CN107238824B (en) | Satellite-borne SAR image geometric accurate correction method based on priori dem data | |
CN109100718A (en) | Sparse aperture ISAR self-focusing and transverse calibration method based on Bayesian learning | |
CN110865346B (en) | Satellite-borne SAR time parameter calibration method based on direct positioning algorithm | |
CN107340434B (en) | Satellite-borne synthetic aperture radar antenna direction picture capturing method based on internal calibration data | |
CN111896954A (en) | Corner reflector coordinate positioning method for shipborne SAR image | |
Chisham | Calibrating SuperDARN interferometers using meteor backscatter | |
CN110687508A (en) | Method for correcting monitoring data of micro-varying radar | |
CN115079109A (en) | Synthetic aperture radar imaging quality index evaluation system | |
CN115865142B (en) | Radar antenna beam pointing calibration method, system, device and storage medium | |
CN109959933B (en) | Multi-baseline circular synthetic aperture radar imaging method based on compressed sensing | |
CN112415512A (en) | SAR moving target focusing method based on advance and retreat method and golden section method | |
CN109061645B (en) | Virtual correlation sparse imaging method and system for synthetic aperture radiometer | |
CN116500648A (en) | Wind profile inversion method for foundation laser radar target area | |
CN102495408B (en) | Method for automatically searching point of synthetic aperture radar dot-matrix target image data | |
RU2723706C1 (en) | Method of obtaining a two-dimensional radar image of an object in multi-frequency pulse sounding and inverse synthesis of an aperture with determination of the third coordinate of the elements of the formed image | |
RU2327190C1 (en) | Method of obtaining three-dimensional radar image of rotating, pitching and banking object, with multi-frequency pulsed probing | |
Naething et al. | Analysis of SAR autofocus performance | |
CN114488142A (en) | Radar two-dimensional angle imaging method and system based on difference-sum beam | |
Leifer et al. | Computing RMS and Integrated Array Sidelobes | |
US12025693B2 (en) | Method for correcting a synthetic aperture radar antenna beam image | |
Daoud et al. | Simulated Data for SAR Tomographic Reconstruction in Urban Areas | |
CN114924239A (en) | Rapid double-station SAR echo simulation method based on electromagnetic scattering model | |
CN109581377B (en) | Resolution evaluation method and device for point echo imaging quality and storage medium |
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: 20201106 |
|
RJ01 | Rejection of invention patent application after publication |