A method and device for reverse recognition of objects in synthetic aperture radar satellite echo domain data

Abstract

The invention discloses a method and device for reverse cognition of synthetic aperture radar satellite echo domain data objects. The method comprises: obtaining a target shape structure information set; the target shape structure information set comprises target shape structure information; the target shape structure information comprises three-dimensional size information, surface material information and shape information; performing grid discretization processing on the target shape structure information set to obtain a target discrete grid information set; obtaining synthetic aperture radar information and detection echo data of a synthetic aperture radar satellite; the synthetic aperture radar information comprises synthetic aperture radar position information and synthetic aperture radar waveform information; the detection echo data is radar echo data obtained after the synthetic aperture radar of the synthetic aperture radar satellite detects a ground object; the radar echo data is two-dimensional matrix data; performing comparison and prediction processing on the target discrete grid information set, the synthetic aperture radar information and the detection echo data of the synthetic aperture radar to obtain synthetic aperture radar echo object shape structure information.

Description

Technical Field

[0001] The present invention relates to the technical field of satellite data interpretation, and in particular to a method and device for reverse cognition of synthetic aperture radar satellite echo domain data objects. Background Art

[0002] As an active microwave remote sensing detection method, synthetic aperture radar (SAR) has the advantages of long range, wide coverage, long continuous working time, and diverse imaging modes. How to achieve rapid and accurate interpretation of targets based on SAR satellite images is an urgent problem to be solved in the field of remote sensing. Due to the poor visual readability and difficulty of visual interpretation of SAR satellite images, extensive research has been carried out in the field of SAR image interpretation at home and abroad, but most methods are based on the image domain and are limited to the understanding and matching of image features. They do not fully utilize the electromagnetic scattering mechanism information contained in SAR satellite radar imaging data, resulting in the inability to reconstruct the three-dimensional size and structure of objects. Summary of the invention

[0003] In order to improve the accuracy of SAR satellite image object interpretation, the present invention proposes a synthetic aperture radar satellite echo domain data object inverse cognition method and device.

[0004] In a first aspect of the embodiments of the present application, a method for inverse recognition of objects in synthetic aperture radar satellite echo domain data is disclosed, comprising:

[0005] S1, obtaining a target shape structure information set; the target shape structure information set includes target shape structure information; the target shape structure information includes three-dimensional size information, surface material information and shape information;

[0006] S2, performing grid discretization processing on the target shape structure information set to obtain a target discrete grid information set;

[0007] S3, obtaining synthetic aperture radar information and detection echo data of synthetic aperture radar satellite; the synthetic aperture radar information includes synthetic aperture radar position information and synthetic aperture radar waveform information; the detection echo data is radar echo data obtained after the synthetic aperture radar of the synthetic aperture radar satellite detects the ground object; the radar echo data is two-dimensional matrix data;

[0008] S4, performing comparison and prediction processing on the target discrete grid information set, synthetic aperture radar information and synthetic aperture radar detection echo data to obtain synthetic aperture radar echo object shape structure information.

[0009] The grid discretization processing is performed on the target shape structure information set to obtain the target discrete grid information set, including:

[0010] S21, performing grid discretization processing on each target shape structure information in the target shape structure information set to obtain corresponding target discrete grid information;

[0011] S22, constructing a target discrete grid information set using the target discrete grid information corresponding to all the target shape structure information of the target shape structure information set; the target discrete grid information set includes the target discrete grid information and the corresponding target shape structure information.

[0012] The step of performing grid discretization processing on each target shape structure information in the target shape structure information set to obtain corresponding target discrete grid information includes:

[0013] S211, using the acquired target shape and structure information, establishing a target three-dimensional model; the target three-dimensional model includes three-dimensional size information, structure distribution information and medium information of the target;

[0014] S212, processing the target three-dimensional model to obtain target discrete grid information.

[0015] The processing of the target three-dimensional model to obtain target discrete grid information includes:

[0016] S2121, performing segmentation processing on the target three-dimensional model to obtain basic facet information of the target and calculation area grid information of the target;

[0017] S2122, using the target calculation area grid information, perform grid generation processing on the basic surface element to obtain the target discrete grid information.

[0018] The target discrete grid information set, synthetic aperture radar information and synthetic aperture radar detection echo data are compared and predicted to obtain synthetic aperture radar echo object shape structure information, including:

[0019] S41, performing radar echo calculation processing using the target discrete grid information set and the synthetic aperture radar information to obtain a target radar echo set; the target radar echo set includes radar echo calculation data corresponding to each target shape structure information;

[0020] S42, for each radar echo calculation data in the target radar echo set, respectively calculating a difference value between the radar echo calculation data and the detection echo data;

[0021] S43, determining the target shape structure information corresponding to the radar echo calculation data having the smallest difference value with the detection echo data, as the synthetic aperture radar echo object shape structure information.

[0022] The radar echo calculation process is performed using the target discrete grid information set and the synthetic aperture radar information to obtain the target radar echo set, including:

[0023] S411, performing virtual incident wave construction processing on each target discrete grid information of the target discrete grid information set and synthetic aperture radar information, respectively, to obtain a virtual incident wave set corresponding to the target discrete grid information;

[0024] S412, for each target discrete grid information of the target discrete grid information set, using a far-zone scattering field calculation model, calculating and processing the target discrete grid information and the corresponding virtual incident wave set to obtain a target radar echo of the target shape structure information corresponding to the target discrete grid information;

[0025] S413, constructing a target radar echo set by using the target radar echoes with all target shape and structure information.

[0026] The step of performing virtual incident wave construction processing on each target discrete grid information of the target discrete grid information set and synthetic aperture radar information to obtain a virtual incident wave set corresponding to the target discrete grid information comprises:

[0027] S4111, for each target discrete grid information of the target discrete grid information set, using synthetic aperture radar position information, connects a line with the coordinates of each grid center point in the target discrete grid information to obtain a corresponding virtual incident wave;

[0028] S4112, using all virtual incident waves, construct a virtual incident wave set corresponding to the target discrete grid information.

[0029] In a second aspect of the embodiment of the present application, a synthetic aperture radar satellite echo domain data object inverse recognition device is disclosed, characterized in that the device comprises:

[0030] A memory storing executable program code;

[0031] a processor coupled to the memory;

[0032] The processor calls the executable program code stored in the memory to execute the synthetic aperture radar satellite echo domain data object inverse cognition method.

[0033] In a third aspect of an embodiment of the present application, a computer storable medium is disclosed, wherein the computer storable medium stores computer instructions, and when the computer instructions are called by a computer, they are used to execute the synthetic aperture radar satellite echo domain data object inverse cognition method.

[0034] In a fourth aspect of the embodiments of the present application, an information data processing terminal is disclosed, which is used to implement the synthetic aperture radar satellite echo domain data object inverse cognition method.

[0035] The beneficial effects of the present invention are:

[0036] The invention discloses a synthetic aperture radar satellite echo domain data object inverse cognition method, which is aimed at the field of SAR satellite data interpretation. The invention aims at the fact that the current SAR image interpretation method is mainly limited to image domain data and does not fully utilize the electromagnetic scattering mechanism information of the object. The invention proposes an electromagnetic inverse cognition method based on echo domain data, which can directly invert the size, main structure and other information of the object from the echo domain data, thereby helping to achieve high-precision interpretation of SAR satellite data. The effectiveness of the invention is verified by inverting the measured SAR images of typical regular scatterers. BRIEF DESCRIPTION OF THE DRAWINGS

[0037] Figure 1 It is a flow chart for implementing the method of the present invention. DETAILED DESCRIPTION

[0038] In order to better understand the content of the present invention, an embodiment is given here.

[0039] Figure 1 It is a flow chart for implementing the method of the present invention.

[0040] In a first aspect of the embodiments of the present application, a method for inverse recognition of objects in synthetic aperture radar satellite echo domain data is disclosed, comprising:

[0041] S1, obtaining a target shape structure information set; the target shape structure information set includes target shape structure information; the target shape structure information includes three-dimensional size information, surface material information and shape information;

[0042] S2, performing grid discretization processing on the target shape structure information set to obtain a target discrete grid information set;

[0043] S3, obtaining synthetic aperture radar information and detection echo data of synthetic aperture radar satellite; the synthetic aperture radar information includes synthetic aperture radar position information and synthetic aperture radar waveform information; the detection echo data is radar echo data obtained after the synthetic aperture radar of the synthetic aperture radar satellite detects the ground object; the radar echo data is two-dimensional matrix data;

[0044] S4, comparing and predicting the target discrete grid information set, the synthetic aperture radar information and the detection echo data of the synthetic aperture radar to obtain the shape and structure information of the synthetic aperture radar echo object;

[0045] The S2 comprises:

[0046] S21, performing grid discretization processing on each target shape structure information in the target shape structure information set to obtain corresponding target discrete grid information;

[0047] S22, constructing a target discrete grid information set using target discrete grid information corresponding to all target shape structure information of the target shape structure information set; the target discrete grid information set includes target discrete grid information and corresponding target shape structure information;

[0048] The S21 includes:

[0049] S211, using the acquired target shape and structure information, establish a target three-dimensional model; the target three-dimensional model includes three-dimensional size information, structure distribution information, medium information, etc. of the target.

[0050] S212, processing the target three-dimensional model to obtain target discrete grid information.

[0051] The processing of the target three-dimensional model to obtain target discrete grid information includes:

[0052] S2121, performing segmentation processing on the target three-dimensional model to obtain basic facet information of the target and calculation area grid information of the target;

[0053] S2122, using the target calculation area grid information, perform grid generation processing on the basic surface element to obtain the target discrete grid information.

[0054] The method of comparing and predicting the target discrete grid information set, the synthetic aperture radar information and the detection echo data of the synthetic aperture radar to obtain the shape and structure information of the synthetic aperture radar echo object includes:

[0055] S41, performing radar echo calculation processing using the target discrete grid information set and the synthetic aperture radar information to obtain a target radar echo set; the target radar echo set includes radar echo calculation data corresponding to each target shape structure information;

[0056] S42, for each radar echo calculation data in the target radar echo set, respectively calculating a difference value between the radar echo calculation data and the detection echo data;

[0057] S43, determining the target shape structure information corresponding to the radar echo calculation data having the smallest difference value with the detection echo data, as the synthetic aperture radar echo object shape structure information.

[0058] The radar echo calculation process is performed using the target discrete grid information set and the synthetic aperture radar information to obtain the target radar echo set, including:

[0059] S411, performing virtual incident wave construction processing on each target discrete grid information of the target discrete grid information set and synthetic aperture radar information, respectively, to obtain a virtual incident wave set corresponding to the target discrete grid information;

[0060] S412, for each target discrete grid information of the target discrete grid information set, using a far-zone scattering field calculation model, calculating and processing the target discrete grid information and the corresponding virtual incident wave set to obtain a target radar echo of the target shape structure information corresponding to the target discrete grid information;

[0061] S413, constructing a target radar echo set by using the target radar echoes with all target shape and structure information.

[0062] The step of performing virtual incident wave construction processing on each target discrete grid information of the target discrete grid information set and synthetic aperture radar information to obtain a virtual incident wave set corresponding to the target discrete grid information comprises:

[0063] S4111, for each target discrete grid information of the target discrete grid information set, using synthetic aperture radar position information, connects a line with the coordinates of each grid center point in the target discrete grid information to obtain a corresponding virtual incident wave;

[0064] S4112, using all virtual incident waves, construct a virtual incident wave set corresponding to the target discrete grid information.

[0065] The expression of the far-zone scattering field calculation model is:

[0066]

[0067] in, represents the virtual incident wave Far-field electromagnetic scattering intensity information at time t, t represents the time variable, N represents the integrated electromagnetic constant vector, τ represents the incident pulse width in the synthetic aperture radar waveform information, Q 1 and Q 2 Represent the first integral function and the second integral function respectively, and their expressions are:

[0068]

[0069] Among them, [L 1 ,L 2] are the upper and lower limits of the integration range of the x-axis; using the virtual incident wave and the two nearest virtual incident waves, three intersection points of the three virtual incident waves and the grid in the target discrete grid information are obtained, and using the two connecting lines between the three intersection points, the corresponding two straight line equations are determined, and their expressions are y=q 1 x+w 1 and y=q 2 x+w 2 ,q 1 、w 1 ,q 2 、w 2 are the coefficients of two straight line equations respectively; the intersection of the virtual incident wave and the grid is connected with the integral point on the grid to obtain the first straight line equation, the expression of the first straight line equation is y=qx+w, q and w represent the coefficients of the first straight line equation; τ 1 represents the time taken for the virtual incident wave to pass through the intersection of the virtual incident wave and the grid from the synthetic aperture radar to the integration point of the grid, t 0 is the incident pulse width control coefficient of the synthetic aperture radar; the first integral function and the second integral function both integrate the variable t; the expression of the integral electromagnetic constant vector N is:

[0070]

[0071] in, Represents the unit external normal vector of the basic surface element corresponding to the mesh, represents the magnetic field unit vector in the virtual incident wave, Z 0 is the free space wave impedance, c is the vacuum wave speed, Represents the virtual incident wave corresponding vector The modulus of Modulo value, we get is the virtual incident wave Far-field electromagnetic scattering intensity information of the corresponding grid;

[0072] By using the far-field electromagnetic scattering intensity information of the grids corresponding to all virtual incident waves, a target radar echo is constructed to obtain the target shape structure information corresponding to the target discrete grid information.

[0073] Among them, [L 1 ,L 2 ] also represents the value range of the vertex on the X-axis when the local discrete grid is projected onto the XOY plane of the target rectangular coordinate system, usually t 0 =0.8τ. The intersection point between the virtual incident wave and the grid can be determined by using the method from S212211 to S212218.

[0074] The calculation expression of the difference value is:

[0075]

[0076] Among them, e i The difference between the i-th radar echo calculation data in the target radar echo set and the detection echo data, A i ij The i-th row and j-th column element of the matrix corresponding to the i-th radar echo calculation data, M and N are the row and column dimensions respectively, B ij is the element of the i-th row and j-th column of the matrix corresponding to the detection echo data.

[0077] The synthetic aperture radar echo object shape structure information is used to characterize the prediction result of the three-dimensional size information, surface material information and shape information of the target object in the synthetic aperture radar echo.

[0078] The detection echo data of the synthetic aperture radar is a target SAR image obtained at the position of the synthetic aperture radar.

[0079] The far-field scattering field calculation model uses the basic surface element intersecting with the virtual incident wave as the integral area for integral calculation. The two lines between the three intersections can be selected from any two lines among all the lines of the three intersections. The synthetic aperture radar waveform can be a Gaussian pulse.

[0080] The S2121 includes:

[0081] S21211, establish a three-dimensional rectangular coordinate system with the center point of the target three-dimensional model as the origin; use basic surface elements to segment the target three-dimensional model to obtain basic surface element information of the target; the basic surface element information includes the number of basic surface elements, the normal vector information of each basic surface element, the vertex coordinate information of each basic surface element, and the medium information of each basic surface element; the medium information includes dielectric constant, electrical conductivity, magnetic permeability and magnetic permeability.

[0082] The step S21211 may be to use 3ds Max software to implement segmentation processing on the target 3D model. From 3ds Max software, export the basic facet metadata file of the target 3D model in stl format; set the number of processes n, use the main process of the computing server to read the stl file, and determine the number of basic facets N by counting the number of times "endfacet" appears in the stl file. tri , get the normal vector information and vertex coordinate information of each basic surface element from the statistical stl file.

[0083] The step S21211 may also be implemented using HyperMesh software.

[0084] S21212, using vertex coordinate information of the basic facets to determine the region range of the basic facets, and using the region ranges of all the basic facets to determine the region range of the target three-dimensional model;

[0085] S21213, using the preset highest calculation frequency and the target fine structure size L min , determine the grid discrete size value;

[0086] The S21213 includes:

[0087] S212131, based on the highest calculation frequency a max and the conductivity of air σ, the minimum wavelength b of the detection wave in air is calculated min , its calculation expression is:

[0088]

[0089] Among them, μ 0 is the magnetic permeability coefficient in vacuum; the detection wave may be an electromagnetic wave.

[0090] S212132, select L min and λ min The minimum value among / 20 is the grid discrete size value;

[0091] S21214, determining the target calculation area grid information according to the grid discrete size value and the area range of the target three-dimensional model;

[0092] The S21214 includes:

[0093] S212141, using the area range of the target three-dimensional model, divided by the grid discrete size value, to obtain a calculation area of ​​the target three-dimensional model;

[0094] S212142, perform unit discretization processing on the calculation area of ​​the target three-dimensional model to obtain the calculation area grid information of the target; the calculation area grid includes a plurality of area grid lines; the calculation area grid information includes equations of the area grid lines and range information of the calculation area grid; the range information of the calculation area grid is expressed in the three-dimensional rectangular coordinate system as (Xmin: 1: Xmax, Ymin: 1: Ymax, Zmin: 1: Zmax); the area grid line is a combination of two variables selected from the three variables of the X-axis, Y-axis and Z-axis of the three-dimensional rectangular coordinate system, and the two variables are selected from the area range (Xmin: 1: Xmax, Ymin: 1: Ymax, Zmin: 1: Zmax) according to the coordinate axes where they are located, and the resulting straight line is constructed; the value set of the endpoints of the area grid line on the three axes constitutes the range information of the calculation area grid.

[0095] The S212, after completing S21214 and before performing the grid generation process, further includes:

[0096] S21215, dividing the calculation area of ​​the target three-dimensional model to obtain three-dimensional virtual topological structure information of the target characteristics; the three-dimensional virtual topological structure information includes range information of the sub-areas obtained by dividing; according to the three-dimensional virtual topological structure information, each sub-area is allocated to a corresponding calculation process of the calculation server to achieve optimal communication and load balancing;

[0097] S21216, according to the sub-region corresponding to each calculation process and the vertex coordinate information of the basic surface element, the basic surface elements included in each calculation process are statistically obtained.

[0098] The S21216 includes: calculating the center coordinate information of the basic face element according to the vertex coordinate information of the basic face element; judging the calculation process to which the basic face element belongs according to the center coordinate information of the basic face element and the corresponding sub-area of ​​each calculation process, and obtaining the basic face element information contained in each calculation process by counting, which is specifically the position and quantity.

[0099] The S2122 includes:

[0100] S21221, using the target calculation area grid information, performing intersection judgment processing on the basic face elements to obtain the intersection points of all basic face elements and the regional grid lines of the calculation area grid;

[0101] S21222, using the intersection of the basic facet and the regional grid line of the calculation regional grid, numbering each regional grid line to obtain local discrete grid information of each regional grid line; the local discrete grid information includes grid position information and grid center point coordinate information;

[0102] The intersection point between the virtual incident wave and the grid may be determined by methods from S212211 to S212218.

[0103] Specifically, step S2122 may be, in each calculation process, performing mesh generation processing on the included basic surface elements to obtain a local discrete mesh of the target;

[0104] S21223, integrating the local discrete grid information of all regional grid lines to obtain target discrete grid information.

[0105] The S21223 may be to use the local discrete grids of all regional grid lines as the target discrete grid information.

[0106] The S21223 may be that each computing process sends the local discrete grid of the target to the main process of the computing server, and the main process integrates all collected grid data to achieve rapid generation of the target discrete grid.

[0107] The method of using the target calculation area grid information to perform intersection judgment processing on the basic face elements to obtain the intersection points of all basic face elements and the regional grid lines of the calculation area grid includes:

[0108] S212211, determining the absolute values ​​of the components of the normal vector of the basic face element on the three coordinate axes of the three-dimensional rectangular coordinate system, and determining the coordinate axes corresponding to the two components with the smallest absolute values; constructing a projection plane using the coordinate axes corresponding to the two components; and projecting the basic face element onto the projection plane of the three-dimensional rectangular coordinate system to obtain a projection basic shape;

[0109] S212212, constructing a projected rectangular area using the vertex coordinates of the projected basic shape within the value range of the coordinate axes corresponding to the two components;

[0110] S212213, judging whether the projected rectangular area contains regional grid lines, and obtaining a judgment result;

[0111] If the result of the judgment is yes, execute S212214; if the result of the judgment is no, execute S212218;

[0112] The projection rectangular area contains regional grid lines, which means that the regional grid lines pass through the projection rectangular area.

[0113] S212214, determining the intersection of the regional grid lines contained in the projected rectangular area with the projection plane;

[0114] S212215, using the intersection point as a vector starting point and using the vertices of the projection basic shape as vector end points respectively, to construct a projection vector set; the projection vector set includes a plurality of projection vectors;

[0115] S212216, constructing a projection vector combination set by using all combinations of two projection vectors in the projection vector set;

[0116] S212217, for each projection vector combination in the projection vector combination set, respectively calculate the corresponding vector cross product value to obtain a vector cross product value set;

[0117] Determine whether the vector cross product value set meets the cross discrimination condition, and obtain a cross discrimination result;

[0118] When the intersection determination result satisfies the intersection determination condition, determining that the intersection point of the regional grid line and the basic bin is inside the basic bin;

[0119] When the intersection determination result is that the intersection determination condition is not satisfied, determining that the intersection point between the regional grid line and the basic bin is not inside the basic bin;

[0120] The crossover discrimination condition is that each vector cross product value in the vector cross product value set is a positive number, or each vector cross product value in the vector cross product value set is a negative number, or one vector cross product value in the vector cross product value set is 0, or two vector cross product values ​​in the vector cross product value set are 0.

[0121] S212218, if the basic bin is not the last basic bin in the preset sequence, execute S212211 for the next basic bin in the preset sequence;

[0122] If the basic surface element is the last basic surface element according to a preset sequence, the intersection points of all basic surface elements and the regional grid lines of the calculation regional grid are obtained.

[0123] The preset order may be the spatial distribution order of the basic face elements in the target three-dimensional model, for example, in an order from top to bottom and from left to right.

[0124] The method of using the intersection of the basic facet and the regional grid line of the calculation regional grid to perform numbering processing on each regional grid line to obtain a local discrete grid of each regional grid line includes:

[0125] S212221, for each regional grid line, numbering the intersections in order from small to large absolute values ​​of the coordinates of the intersections contained therein to obtain numbered intersections;

[0126] S212222, construct a directed line segment using all two numbered intersections with adjacent numbers; in the directed line segment, the intersection with a smaller number among the two numbered intersections is used as the starting point of the directed line segment, and the intersection with a larger number among the two numbered intersections is used as the end point of the directed line segment;

[0127] S212223, determining all directed line segments whose starting points are odd-numbered and whose end points are even-numbered as local discrete grids of the regional grid lines.

[0128] The S21222 may be: for each regional grid line, in the order of the absolute values ​​of the coordinates of the intersections it contains from small to large, n numbered intersections are obtained, 1, 2, 3, 4, ..., n, the directed line segment between intersection 1 and intersection 2 is Ⅰ, the directed line segment between intersection 2 and intersection 3 is Ⅱ, the directed line segment between intersection 3 and intersection 4 is Ⅲ, and so on. The directed line segments numbered Ⅰ, Ⅲ, Ⅴ, ... to the last odd number less than or equal to n are the local discrete grids of the regional grid line; for the local discrete grids, the corresponding medium parameters are marked.

[0129] The acquisition of the target shape and structure information may be achieved through an Internet database, data collection on the target, or a method of generating simulation data.

[0130] The target data collection may be achieved by radar measurement, laser measurement or acoustic measurement.

[0131] The method for generating the simulation data may be to use mechanical simulation software.

[0132] The three-dimensional size information of the target may be the length, width and height; the component information includes the power part, load part, peripheral part and fine structure part; the surface material information includes the type, parameters and thickness of the surface material, etc. The shape information includes the specific three-dimensional size information of each component.

[0133] The fine structure may be a mast, an antenna, a sensor, etc.

[0134] The step S2 can be implemented by three-dimensional mechanical simulation software; the three-dimensional mechanical simulation software can be 3ds Max software.

[0135] The target electromagnetic echo image is obtained by actual measurement using SAR radar.

[0136] The basic face element may be a triangle face element, a rhombus face element, a regular pentagon face element, etc.

[0137] In a second aspect of the embodiment of the present application, a synthetic aperture radar satellite echo domain data object inverse recognition device is disclosed, characterized in that the device comprises:

[0138] A memory storing executable program code;

[0139] a processor coupled to the memory;

[0140] The processor calls the executable program code stored in the memory to execute the synthetic aperture radar satellite echo domain data object inverse cognition method.

[0141] In a third aspect of an embodiment of the present application, a computer storable medium is disclosed, wherein the computer storable medium stores computer instructions, and when the computer instructions are called by a computer, they are used to execute the synthetic aperture radar satellite echo domain data object inverse cognition method.

[0142] In a fourth aspect of the embodiments of the present application, an information data processing terminal is disclosed, which is used to implement the synthetic aperture radar satellite echo domain data object inverse cognition method.

[0143] The above is only an embodiment of the present application and is not intended to limit the present application. For those skilled in the art, the present application may have various changes and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application should be included in the scope of the claims of the present application.

Claims

1. A method for inverse recognition of objects in synthetic aperture radar satellite echo domain data, characterized in that: include: S1, obtain the target shape structure information set; The target shape structure information set includes target shape structure information; The target shape structure information includes three-dimensional size information, surface material information and shape information; S2, performing grid discretization processing on the target shape structure information set to obtain a target discrete grid information set; S3, obtaining synthetic aperture radar information and detection echo data of synthetic aperture radar satellites; the synthetic aperture radar information includes synthetic aperture radar position information and synthetic aperture radar waveform information; The detection echo data is radar echo data obtained after the synthetic aperture radar of the synthetic aperture radar satellite detects the ground object; The radar echo data is two-dimensional matrix data; S4, performing comparison and prediction processing on the target discrete grid information set, synthetic aperture radar information and synthetic aperture radar detection echo data to obtain synthetic aperture radar echo object shape structure information; The method of comparing and predicting the target discrete grid information set, the synthetic aperture radar information and the detection echo data of the synthetic aperture radar to obtain the synthetic aperture radar echo object shape structure information includes: S41, performing radar echo calculation processing using the target discrete grid information set and the synthetic aperture radar information to obtain a target radar echo set; the target radar echo set includes radar echo calculation data corresponding to each target shape structure information; S42, for each radar echo calculation data in the target radar echo set, respectively calculating a difference value between the radar echo calculation data and the detection echo data; S43, determining the target shape structure information corresponding to the radar echo calculation data having the smallest difference value with the detection echo data as the synthetic aperture radar echo object shape structure information; The radar echo calculation process is performed using the target discrete grid information set and the synthetic aperture radar information to obtain the target radar echo set, including: S411, performing virtual incident wave construction processing on each target discrete grid information of the target discrete grid information set and synthetic aperture radar information, respectively, to obtain a virtual incident wave set corresponding to the target discrete grid information; S412, for each target discrete grid information of the target discrete grid information set, using a far-zone scattering field calculation model, calculating and processing the target discrete grid information and the corresponding virtual incident wave set to obtain a target radar echo of the target shape structure information corresponding to the target discrete grid information; S413, constructing a target radar echo set by using target radar echoes with all target shape and structure information; The step of performing virtual incident wave construction processing on each target discrete grid information of the target discrete grid information set and synthetic aperture radar information to obtain a virtual incident wave set corresponding to the target discrete grid information comprises: S4111, for each target discrete grid information of the target discrete grid information set, using synthetic aperture radar position information, connects a line with the coordinates of each grid center point in the target discrete grid information to obtain a corresponding virtual incident wave; S4112, constructing a virtual incident wave set corresponding to the target discrete grid information using all virtual incident waves; The calculation expression of the difference value is: Wherein, ei is the difference between the calculated data of the i-th radar echo in the target radar echo set and the detection echo data, A i ij The i-th row and j-th column element of the matrix corresponding to the i-th radar echo calculation data, M and N are the row and column dimensions respectively, B ij is the element of the i-th row and j-th column of the matrix corresponding to the detection echo data.

2. The synthetic aperture radar satellite echo domain data object inverse recognition method according to claim 1, characterized in that: The grid discretization processing is performed on the target shape structure information set to obtain the target discrete grid information set, including: S21, performing grid discretization processing on each target shape structure information in the target shape structure information set to obtain corresponding target discrete grid information; S22, constructing a target discrete grid information set using the target discrete grid information corresponding to all the target shape structure information of the target shape structure information set; the target discrete grid information set includes the target discrete grid information and the corresponding target shape structure information.

3. The method for inverse recognition of objects in synthetic aperture radar satellite echo domain data according to claim 2, characterized in that: The step of performing grid discretization processing on each target shape structure information in the target shape structure information set to obtain corresponding target discrete grid information includes: S211, using the acquired target shape and structure information, establishing a target three-dimensional model; the target three-dimensional model includes three-dimensional size information, structure distribution information and medium information of the target; S212, processing the target three-dimensional model to obtain target discrete grid information.

4. The synthetic aperture radar satellite echo domain data object inverse recognition method according to claim 3, characterized in that: The processing of the target three-dimensional model to obtain target discrete grid information includes: S2121, performing segmentation processing on the target three-dimensional model to obtain basic facet information of the target and calculation area grid information of the target; S2122, using the target calculation area grid information, perform grid generation processing on the basic surface element to obtain the target discrete grid information.

5. A synthetic aperture radar satellite echo domain data object reverse recognition device, characterized in that: The device comprises: A memory storing executable program code; a processor coupled to the memory; The processor calls the executable program code stored in the memory to execute the synthetic aperture radar satellite echo domain data object inverse cognition method according to any one of claims 1 to 4.

6. A computer storable medium, characterized in that: The computer storable medium stores computer instructions, and when the computer instructions are called by a computer, they are used to execute the synthetic aperture radar satellite echo domain data object inverse recognition method according to any one of claims 1 to 4.

7. An information data processing terminal, characterized in that: The information data processing terminal is used to implement the synthetic aperture radar satellite echo domain data object inverse recognition method as described in any one of claims 1 to 4.

Images