CN116754847B - Method and device for estimating electromagnetic scattering intensity of far-region of sea surface composite target - Google Patents
Method and device for estimating electromagnetic scattering intensity of far-region of sea surface composite target Download PDFInfo
- Publication number
- CN116754847B CN116754847B CN202310668652.1A CN202310668652A CN116754847B CN 116754847 B CN116754847 B CN 116754847B CN 202310668652 A CN202310668652 A CN 202310668652A CN 116754847 B CN116754847 B CN 116754847B
- Authority
- CN
- China
- Prior art keywords
- target
- virtual
- intersection
- vector
- information
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 36
- 239000002131 composite material Substances 0.000 title claims abstract description 29
- 239000013598 vector Substances 0.000 claims abstract description 133
- 238000012545 processing Methods 0.000 claims abstract description 42
- 238000005259 measurement Methods 0.000 claims abstract description 17
- 230000014509 gene expression Effects 0.000 claims description 33
- 230000006870 function Effects 0.000 claims description 10
- 230000011218 segmentation Effects 0.000 claims description 10
- 238000012937 correction Methods 0.000 claims description 9
- 230000010354 integration Effects 0.000 claims description 3
- 238000003079 width control Methods 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 5
- 238000000638 solvent extraction Methods 0.000 description 4
- PCTMTFRHKVHKIS-BMFZQQSSSA-N (1s,3r,4e,6e,8e,10e,12e,14e,16e,18s,19r,20r,21s,25r,27r,30r,31r,33s,35r,37s,38r)-3-[(2r,3s,4s,5s,6r)-4-amino-3,5-dihydroxy-6-methyloxan-2-yl]oxy-19,25,27,30,31,33,35,37-octahydroxy-18,20,21-trimethyl-23-oxo-22,39-dioxabicyclo[33.3.1]nonatriaconta-4,6,8,10 Chemical compound C1C=C2C[C@@H](OS(O)(=O)=O)CC[C@]2(C)[C@@H]2[C@@H]1[C@@H]1CC[C@H]([C@H](C)CCCC(C)C)[C@@]1(C)CC2.O[C@H]1[C@@H](N)[C@H](O)[C@@H](C)O[C@H]1O[C@H]1/C=C/C=C/C=C/C=C/C=C/C=C/C=C/[C@H](C)[C@@H](O)[C@@H](C)[C@H](C)OC(=O)C[C@H](O)C[C@H](O)CC[C@@H](O)[C@H](O)C[C@H](O)C[C@](O)(C[C@H](O)[C@H]2C(O)=O)O[C@H]2C1 PCTMTFRHKVHKIS-BMFZQQSSSA-N 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 239000010432 diamond Substances 0.000 description 1
- 230000005672 electromagnetic field Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R29/00—Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
- G01R29/08—Measuring electromagnetic field characteristics
- G01R29/0864—Measuring electromagnetic field characteristics characterised by constructional or functional features
- G01R29/0892—Details related to signal analysis or treatment; presenting results, e.g. displays; measuring specific signal features other than field strength, e.g. polarisation, field modes, phase, envelope, maximum value
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- General Physics & Mathematics (AREA)
- Computer Graphics (AREA)
- Theoretical Computer Science (AREA)
- Remote Sensing (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses a method and a device for estimating electromagnetic scattering intensity of a sea surface composite target remote zone, wherein the method comprises the following steps: acquiring target range position information and incident point source information; determining a virtual incident wave set according to the target range position information and the incident point source information; performing subdivision and measurement processing on the target three-dimensional model to obtain a target triangular surface element model; processing the virtual incident wave set and the target triangular surface element model by using an intersection discrimination model to obtain an intersection discrimination result; and processing the intersection discrimination result by using a marine remote-zone scattered field calculation model to obtain the electromagnetic scattering intensity information of the sea surface composite target remote zone. The invention improves the calculation efficiency and saves the memory by reducing the number of the target triangle surface elements; by establishing the intersection point discrimination vector and introducing the target position measurement variance vector information in the intersection point discrimination process, the images of sea surface stormy waves are effectively reduced, and the calculation accuracy is improved.
Description
Technical Field
The invention relates to the technical field of electromagnetic fields and microwaves, in particular to a method and a device for estimating electromagnetic scattering intensity of a far-zone of a sea surface composite target.
Background
Accurate estimation of electromagnetic scattering intensity of a far region of a sea surface and sea composite target is always limited by factors such as application conditions, calculation efficiency, memory and the like. The traditional method has the defects of high calculation cost and huge consumption of calculation resources when calculating the composite scattering problem of the super-electric large-size target and the sea surface. On the premise of keeping the far-region electromagnetic scattering intensity solving precision of the sea surface and the offshore composite target, the calculating efficiency of the solving method is improved, the calculated amount is effectively reduced, and the method is a problem to be solved at present.
Disclosure of Invention
Aiming at the problem of improving the calculation efficiency of the solving method and effectively reducing the calculation amount on the premise of keeping the far-region electromagnetic scattering intensity solving precision of the sea surface and the offshore composite target, the invention discloses a method and a device for estimating the far-region electromagnetic scattering intensity of the sea surface composite target.
The invention discloses a sea surface composite target far-zone electromagnetic scattering intensity estimation method, which comprises the following steps:
s1, acquiring target range position information and incident point source information; the incident point source information comprises incident point source position information and incident point source waveform information; the target range position information comprises target three-dimensional model information, target boundary position information, target center position information and target boundary point position information;
s2, determining a virtual incident wave set according to the target range position information and the incident point source information;
s3, performing subdivision and measurement processing on the target three-dimensional model to obtain a target triangular surface element model; the target triangular surface element model comprises target position measurement variance vector information, triangular surface element quantity information contained in a target, normal vector information of each triangular surface element, vertex coordinate information of each triangular surface element and medium information of each triangular surface element;
s4, processing the virtual incident wave set and the target triangular surface element model by using an intersection judgment model to obtain an intersection judgment result;
and S5, processing the intersection discrimination result by using a remote sea area scattered field calculation model to obtain the remote sea area electromagnetic scattering intensity information of the sea surface composite target.
The determining a virtual incident wave set according to the target range position information and the incident point source information comprises the following steps:
s21, processing the target range position information and the incident point source position information by using a height coverage model to obtain a sea surface preset height h 0 ；
The expression of the height coverage model is as follows:
wherein, (x) s ,y s ,z s ) Is the position coordinate value of the incident point source in the target rectangular coordinate system, (x) a ,y a ,z a ) H is the position coordinate value of the geometric center point of the target in the rectangular coordinate system of the target 1 Altitude, which is the highest point of the target; the target rectangular coordinate system is a three-dimensional rectangular coordinate system established by taking the projection of a geometric center point of a target on the sea level as an origin and taking the sea level as an XOY axis;
s22, presetting the sea surface to be high h 0 Constructing a virtual horizontal plane;
s23, intersecting the virtual horizontal plane by utilizing a connecting line of the incident point source and the target boundary point to obtain a virtual intersection point;
s24, constructing a virtual parallelogram surrounding all virtual intersection points on the virtual horizontal plane;
s25, uniformly dividing the virtual parallelogram by utilizing the basic shape to obtain a virtual divided quadrilateral; the virtual segmentation quadrangle comprises a plurality of basic shapes;
s26, determining the incident point source as a starting point, determining the vertexes of the basic shapes contained in the virtual segmentation quadrangle as path points, and determining rays pointing from the starting point to the path points as virtual incident waves;
s27, constructing and obtaining a virtual incident wave set by utilizing all virtual incident waves.
The processing the virtual incident wave set and the target triangle surface element model by utilizing the intersection discrimination model to obtain an intersection discrimination result comprises the following steps:
s41, measuring variance vector information by using a direction vector and a target position of each virtual incident wave, and calculating to obtain a correction vector corresponding to the virtual incident wave; correction vector corresponding to the virtual incident wave The calculation expression is as follows:
wherein,for the direction vector of the virtual incident wave, < >>The variance vector is measured for the target location,
s42, correcting the vertex coordinate information of the triangular surface element by utilizing the target position measurement variance vector information to obtain triangular surface element vertex vector information;
and correcting the vertex coordinate information of the triangular surface element by utilizing the target position measurement variance vector information to obtain triangular surface element vertex vector information, wherein the calculation expression is as follows:
wherein,vectors corresponding to three vertex coordinates of the triangular surface element, respectively, < ->Three vertex vectors of the triangular surface element respectively;
s43, processing the virtual incident wave and the triangular surface element vertex vector information by utilizing an intersection judgment model to obtain an intersection judgment result of the virtual incident wave;
the processing the virtual incident wave and the triangle primitive vertex vector information by utilizing an intersection discrimination model to obtain an intersection discrimination result of the virtual incident wave comprises the following steps:
s431, utilizing the incident point source vectorCorrection vector corresponding to virtual incident wave>And triangular surface element vertex vector information, calculating to obtain an intersection judgment vector set; the intersection judgment vector set comprises a first intersection judgment vector T 1 Second intersection judgment vector T 2 Third intersection judgment vector T 3 Fourth intersection judgment vector T 4 And a fifth intersection judgment vector T 5 The method comprises the steps of carrying out a first treatment on the surface of the The incident point source vector->
The calculation expression of the intersection judgment vector set is as follows:
T 5 ＝T 3 ×T 1 ；
s432, vector calculation processing is carried out on the intersecting discrimination vector set to obtain a first intersecting discrimination vector (p 1 ,p 2 ,p 3 ) The method comprises the steps of carrying out a first treatment on the surface of the The first intersection discrimination vector (p 1 ,p 2 ,p 3 ) The calculated expression of (2) is:
s433, the first intersection discrimination vector is subjected to discrimination processing, and the intersection discrimination result of the virtual incident wave is obtained.
The step of performing a discrimination process on the first intersection discrimination vector to obtain an intersection discrimination result of the virtual incident wave includes:
judging the first intersection judgment vector (p 1 ,p 2 ,p 3 ) Whether the intersecting condition is met or not, and a condition judgment result is obtained; the intersecting conditions include:
p 1 ≥0,p 2 ≥β 2 ,p 3 ≥β 3 ,1-p 2 -p 3 ≥β 1 ；
when the first intersection discrimination vector (p 1 ,p 2 ,p 3 ) When the intersecting condition is met, determining that the virtual incident wave intersects with the triangular surface element, and calculating to obtain an intersection point coordinate W of the virtual incident wave and the triangular surface element, wherein the calculation expression is as follows:
when the house is atThe first intersection discrimination vector (p 1 ,p 2 ,p 3 ) And when the intersection condition is not met, determining that the virtual incident wave is not intersected with the triangular surface element.
S44, integrating all the intersection discrimination results of the virtual incident waves to obtain the intersection discrimination result of the virtual incident wave set. The intersection discrimination result of the virtual incident wave set comprises intersecting triangular face element information and intersection point coordinate information corresponding to each virtual incident wave.
The processing of the intersecting discrimination result by using the offshore remote region scattered field calculation model to obtain the electromagnetic scattering intensity information of the offshore composite target remote region comprises the following steps:
calculating the corresponding far-zone electromagnetic scattering intensity information of each virtual incident wave by using a marine far-zone scattering field calculation model;
and superposing the electromagnetic scattering intensity of the far-zone of each virtual incident wave to obtain the electromagnetic scattering intensity information of the far-zone of the sea surface composite target.
The offshore remote region fringe field computation model has the expression:
wherein t represents a time variable, N represents an integrated electromagnetic constant vector, τ represents an incident pulse width of an incident point source, and Q 1 And Q 2 The first integral function and the second integral function are respectively expressed, and the expression is as follows:
wherein, [ L ] 1 ,L 2 ]When the triangle surface element is projected onto the XOY plane of the target rectangular coordinate system, the value range of the vertex of the triangle surface element on the X axis is also the upper limit and the lower limit of the integral range of the X axis; obtaining three virtual incident waves and an XOY plane by using the virtual incident waves and two virtual incident waves nearest to the virtual incident wavesThree intersection points, two straight line equations corresponding to the two intersection points are determined by utilizing two connecting lines among the three intersection points, and the expressions of the straight line equations are y=q respectively 1 x+w 1 And y=q 2 x+w 2 ，q 1 、w 1 、q 2 、w 2 Coefficients of two linear equations respectively; connecting the intersection point of the virtual incident wave and the triangular surface element with an integral point on the triangular surface element to obtain a first linear equation, wherein the expression of the first linear equation is y=qx+w, and q and w represent coefficients of the first linear equation; τ 1 Representing the time taken for the virtual incident wave to reach the integration point of the triangular surface element from the incident point source through the intersection point of the virtual incident wave and the triangular surface element, t 0 The control coefficient is the incident pulse time width of the incident point source, and tau is the incident point source waveform time width; the first integral function and the second integral function integrate the variable t; the expression of the integral electromagnetic constant vector N is as follows:
wherein,an extra-unit normal vector representing a triangle bin, < >> A magnetic field unit vector, Z, in a virtual incident wave representing an incident point source 0 Is free space wave impedance, c is vacuum wave velocity, < ->Representing the incident point source vector +.>Is a modulus of the model.
The invention discloses a data processing device, which comprises:
a memory storing executable program code;
a processor coupled to the memory;
and the processor calls the executable program codes stored in the memory to execute the sea surface composite target far-zone electromagnetic scattering intensity estimation method.
The invention discloses a computer storage medium which stores computer instructions for executing the sea surface composite target far-zone electromagnetic scattering intensity estimation method when the computer instructions are called.
The beneficial effects of the invention are as follows:
1. the invention eliminates the limit that the traditional composite target far-field has to meet the eighth wavelength of the target bin size, and improves the calculation efficiency and saves the memory by reducing the number of target triangle bins.
2. By establishing the intersection point discrimination vector, the complexity of the intersection point discrimination process is effectively reduced, and the calculation efficiency is improved. Meanwhile, in the intersection point judging process, the image of sea surface wind wave factors is effectively reduced and the calculation accuracy is improved by introducing the target position measurement variance vector information.
3. According to the invention, when the electromagnetic scattering intensity of the far-region of the sea surface and the offshore composite target is calculated, the integral calculation region is divided into two regions, so that the calculation precision is effectively improved, and the calculated electromagnetic scattering intensity of the far-region is more accurate.
Drawings
FIG. 1 is a flow chart of the method of the present invention;
FIG. 2 is a schematic view of the virtual horizontal plane arrangement of the present invention;
FIG. 3 is a schematic diagram of vertices of a virtual partitioned quadrilateral obtained by uniformly partitioning the virtual graph with triangles;
FIG. 4 is a schematic diagram of vertices of a virtual partitioned quadrilateral obtained by uniformly partitioning the virtual graph with the quadrilateral;
FIG. 5 is a schematic diagram of a virtual split quadrilateral according to the present invention.
Detailed Description
For a better understanding of the present disclosure, an embodiment is presented herein.
FIG. 1 is a flow chart of the method of the present invention; FIG. 2 is a schematic view of the virtual horizontal plane arrangement of the present invention; FIG. 3 is a schematic diagram of vertices of a virtual partitioned quadrilateral obtained by uniformly partitioning the virtual graph with triangles; fig. 4 is a schematic diagram of vertices of a virtual partitioned quadrilateral obtained by uniformly partitioning the virtual graphic with the quadrilateral.
The invention discloses a sea surface composite target far-zone electromagnetic scattering intensity estimation method, which comprises the following steps:
s1, acquiring target range position information and incident point source information; the incident point source information comprises incident point source position information and incident point source waveform information; the target range position information comprises target three-dimensional model information, target boundary position information, target center position information and target boundary point position information;
the target three-dimensional model is established by acquiring target geometric structure information and utilizing the acquired target geometric structure information;
the incident point source information may be obtained by a preset method according to target range information.
S2, determining a virtual incident wave set according to the target range position information and the incident point source information;
s3, performing subdivision and measurement processing on the target three-dimensional model to obtain a target triangular surface element model; the target triangular surface element model comprises target position measurement variance vector information, triangular surface element quantity information contained in a target, normal vector information of each triangular surface element, vertex coordinate information of each triangular surface element and medium information of each triangular surface element;
the step S3 may be implemented by using HyperMesh software.
The target position measurement variance vector information is obtained by calculating variances on three coordinate axes of a target rectangular coordinate system based on a target three-dimensional model on the basis of a target center position measurement result obtained on the sea surface.
S4, processing the virtual incident wave set and the target triangular surface element model by using an intersection judgment model to obtain an intersection judgment result;
and S5, processing the intersection discrimination result by using a remote sea area scattered field calculation model to obtain the remote sea area electromagnetic scattering intensity information of the sea surface composite target.
The step S2 comprises the following steps:
s21, processing the target range position information and the incident point source position information by using a height coverage model to obtain a sea surface preset height h 0 ；
The expression of the height coverage model is as follows:
wherein, (x) s ,y s ,z s ) Is the position coordinate value of the incident point source in the target rectangular coordinate system, (x) a ,y a ,z a ) H is the position coordinate value of the geometric center point of the target in the rectangular coordinate system of the target 1 Altitude, which is the highest point of the target; the target rectangular coordinate system is a three-dimensional rectangular coordinate system established by taking the projection of a geometric center point of a target on the sea level as an origin and taking the sea level as an XOY axis;
through the height coverage model, the height of the virtual horizontal plane can be realized, the height of the virtual horizontal plane is lower when the target is lower along with the self-adaptive adjustment of the relative height relation between the target and the incident point source, and the height of the virtual horizontal plane is higher when the target is higher, so that the complexity of subsequent grid division and intersection point judgment is reduced.
S22, presetting the sea surface to be high h 0 Constructing a virtual horizontal plane;
s23, intersecting the virtual horizontal plane by utilizing a connecting line of the incident point source and the target boundary point to obtain a virtual intersection point;
s24, constructing a virtual parallelogram surrounding all virtual intersection points on the virtual horizontal plane;
s25, uniformly dividing the virtual parallelogram by utilizing the basic shape to obtain a virtual divided quadrilateral; the virtual segmentation quadrangle comprises a plurality of basic shapes;
the basic shape may be triangular or parallelogram;
s26, determining the incident point source as a starting point, determining the vertexes of the basic shapes contained in the virtual segmentation quadrangle as path points, and determining rays pointing from the starting point to the path points as virtual incident waves;
s27, constructing and obtaining a virtual incident wave set by utilizing all virtual incident waves.
The step S2 may be:
the incident wave passing through a virtual horizontal plane, such as a 'B' C 'D' in fig. 2, is discretized into a number of independent virtual incident waves. The direction of the ray from the incident point source to the target center point is then the main direction of the beam, while the virtual horizontal plane is a planar quadrilateral perpendicular to the main direction of the beam and at a distance from the target center point. The target is projected onto the ground (i.e., XOY plane) along the main beam direction and a rectangle such as ABCD in fig. 2 is found on the projection plane to just cover the projection. Four rays can be formed between the incident source point and four vertexes of the rectangle, the four rays intersect with a plane where a virtual horizontal plane is located, a connecting line of the four intersection points is a plane quadrangle, namely a virtual parallelogram, and then the plane quadrangle is divided into virtual incident wave emergent ports independent of each other, namely the required virtual horizontal plane. Two kinds of virtual incident wave dividing modes are introduced here, namely triangle virtual incident wave dividing and quadrilateral virtual incident wave dividing, wherein the two kinds of methods are to uniformly divide a plane quadrilateral by utilizing triangles and diamonds respectively, and are shown in fig. 3 and 4 respectively.
As shown in FIG. 5, for a virtual split quadrilateral ABCD as a space plane quadrilateral, segments AB and CD are equally divided into N-1 segments, AD and BC are equally divided into M-1 segmentsThen the points of the AB and the DC are correspondingly connected, the equal division points corresponding to the AD and the BC are connected to form (M-1) x (N-1) space quadrilaterals, each space quadrilaterals is divided into two triangles, and 2 x (M-1) x (N-1) space triangles are formed in total, namely the space quadrilaterals on the same plane can be divided into 2 * (M-1) triangle virtual incident waves.
All vertexes and triangular surface elements on the virtual horizontal plane are numbered, and as the coordinates of four vertexes A, B, C, D of the space quadrangle are available, after the virtual incident wave is divided, corner point and center point information corresponding to all virtual incident waves can be obtained. Assuming that the triangle virtual incident wave with the bin number 2*i is i, i+1, i+n, respectively, and the triangle virtual incident wave with the bin number 2 x i+1 is i+n, i+1, i+1+n, respectively. The virtual incident wave has three angular point rays and a central ray, and the incident source point is connected with the angular point of each triangle virtual incident wave to form the angular point rays and connected with the central point of the triangle virtual incident wave to form the central ray. The incident wave passes through the virtual horizontal plane to generate a plurality of virtual incident waves which are independent of each other.
The processing the virtual incident wave set and the target triangle surface element model by utilizing the intersection discrimination model to obtain an intersection discrimination result comprises the following steps:
s41, measuring variance vector information by using a direction vector and a target position of each virtual incident wave, and calculating to obtain a correction vector corresponding to the virtual incident wave; correction vector corresponding to the virtual incident wave The calculation expression is as follows:
wherein the method comprises the steps of，For the direction vector of the virtual incident wave, < >>The variance vector is measured for the target location,
s42, correcting the vertex coordinate information of the triangular surface element by utilizing the target position measurement variance vector information to obtain triangular surface element vertex vector information;
and correcting the vertex coordinate information of the triangular surface element by utilizing the target position measurement variance vector information to obtain triangular surface element vertex vector information, wherein the calculation expression is as follows:
wherein,vectors corresponding to three vertex coordinates of the triangular surface element, respectively, < ->Three vertex vectors of the triangular surface element respectively;
s43, processing the virtual incident wave and the triangular surface element vertex vector information by utilizing an intersection judgment model to obtain an intersection judgment result of the virtual incident wave;
the processing the virtual incident wave and the triangle primitive vertex vector information by utilizing an intersection discrimination model to obtain an intersection discrimination result of the virtual incident wave comprises the following steps:
s431, utilizing the incident point source vectorCorrection vector corresponding to virtual incident wave>And triangular surface element vertex vector information, calculating to obtain an intersection judgment vector set; the intersection judgment vector set comprises a first intersection judgment vector T 1 Second intersection judgment vector T 2 Third intersection judgment vector T 3 Fourth intersection judgment vector T 4 And a fifth intersection judgment vector T 5 The method comprises the steps of carrying out a first treatment on the surface of the The incident point source vector->
The calculation expression of the intersection judgment vector set is as follows:
T 5 ＝T 3 ×T 1 ；
s432, vector calculation processing is carried out on the intersecting discrimination vector set to obtain a first intersecting discrimination vector (p 1 ,p 2 ,p 3 ) The method comprises the steps of carrying out a first treatment on the surface of the The first intersection discrimination vector (p 1 ,p 2 ,p 3 ) The calculated expression of (2) is:
s433, the first intersection discrimination vector is subjected to discrimination processing, and the intersection discrimination result of the virtual incident wave is obtained.
The step of performing a discrimination process on the first intersection discrimination vector to obtain an intersection discrimination result of the virtual incident wave includes:
judging the first intersection judgment vector (p 1 ,p 2 ,p 3 ) Whether the intersecting condition is met or not, and a condition judgment result is obtained; the intersecting conditions include:
p 1 ≥0,p 2 ≥β 2 ,p 3 ≥β 3 ,1-p 2 -p 3 ≥β 1 ；
when the first intersection discrimination vector (p 1 ,p 2 ,p 3 ) When the intersecting condition is met, determining that the virtual incident wave intersects with the triangular surface element, and calculating to obtain an intersection point coordinate W of the virtual incident wave and the triangular surface element, wherein the calculation expression is as follows:
when the first intersection discrimination vector (p 1 ,p 2 ,p 3 ) And when the intersection condition is not met, determining that the virtual incident wave is not intersected with the triangular surface element.
S44, integrating all the intersection discrimination results of the virtual incident waves to obtain the intersection discrimination result of the virtual incident wave set. The intersection discrimination result of the virtual incident wave set comprises intersecting triangular face element information and intersection point coordinate information corresponding to each virtual incident wave.
The processing of the intersecting discrimination result by using the offshore remote region scattered field calculation model to obtain the electromagnetic scattering intensity information of the offshore composite target remote region comprises the following steps:
calculating the corresponding far-zone electromagnetic scattering intensity information of each virtual incident wave by using a marine far-zone scattering field calculation model;
and superposing the electromagnetic scattering intensity of the far-zone of each virtual incident wave to obtain the electromagnetic scattering intensity information of the far-zone of the sea surface composite target.
The offshore remote region fringe field computation model has the expression:
wherein t represents a time variable, N represents an integrated electromagnetic constant vector, τ represents an incident pulse width of an incident point source, and Q 1 And Q 2 The first integral function and the second integral function are respectively expressed, and the expression is as follows:
wherein, [ L ] 1 ,L 2 ]When the triangle surface element is projected onto the XOY plane of the target rectangular coordinate system, the value range of the vertex of the triangle surface element on the X axis is also the upper limit and the lower limit of the integral range of the X axis; obtaining three intersection points of the three virtual incident waves and an XOY plane by using the virtual incident waves and two virtual incident waves closest to the virtual incident waves, and determining two corresponding linear equations by using two connecting lines between the three intersection points, wherein the expressions of the two linear equations are y=q respectively 1 x+w 1 And y=q 2 x+w 2 ，q 1 、w 1 、q 2 、w 2 Coefficients of two linear equations respectively; connecting the intersection point of the virtual incident wave and the triangular surface element with an integral point on the triangular surface element to obtain a first linear equation, wherein the expression of the first linear equation is y=qx+w, and q and w represent coefficients of the first linear equation; τ 1 Representing the time taken for the virtual incident wave to reach the integration point of the triangular surface element from the incident point source through the intersection point of the virtual incident wave and the triangular surface element, t 0 For the control factor of the time width of the incident pulse of the incident point source, t is usually taken 0 =0.8τ, τ is the incident point source waveform time width; the first integral function and the second integral function integrate the variable t; the expression of the integral electromagnetic constant vector N is as follows:
wherein,an extra-unit normal vector representing a triangle bin, < >> A magnetic field unit vector, Z, in a virtual incident wave representing an incident point source 0 Is free space wave impedance, c is vacuum wave velocity, < ->Representing the incident point source vector +.>Is a modulus of the model. The offshore remote scattering field calculation model performs integral calculation by taking a triangular surface element intersected with the virtual incident wave as an integral area. The two connecting lines between the three intersecting points are utilized, and two connecting lines are selected from all connecting lines of the three intersecting points. The incident point source waveform may be a gaussian pulse.
The processing of the virtual incident wave set and the target triangle surface element model by using the intersection discrimination model to obtain an intersection discrimination result may be:
performing intersection judgment processing on the virtual incident wave set and the target triangle surface element model to obtain intersection point information of the incident rays, wherein the intersection point information of the incident rays is an intersection judgment result, and specifically comprises the following steps:
surrounding and dividing the target triangular face element model to obtain a target multi-stage cube dividing model;
performing intersection judgment processing on each virtual incident wave in the virtual incident wave set by using the target multi-stage cuboid segmentation model to obtain incident ray intersection point information;
and merging the incident ray intersection point information of all the virtual incident waves to obtain the incident ray intersection point information of the virtual incident wave set.
The surrounding and dividing treatment is carried out on the target triangular surface element model to obtain a target multi-stage cube dividing model, which comprises the following steps:
s401, constructing a first-stage cuboid, wherein the first-stage cuboid surrounds the target triangular surface element model; determining the first-stage cuboid as a current-stage cuboid;
s402, determining a corresponding cutting surface for the current cuboid, and cutting the current cuboid by utilizing the cutting surface to obtain a next cuboid; the number of the next-stage cuboid is increased by 1 for the number of the current cuboid;
the cutting face is used for cutting the current cuboid to obtain the next cuboid, and the method comprises the following steps:
cutting the current-stage cuboid by utilizing the cutting surface to obtain a first cutting cuboid and a second cutting cuboid;
acquiring the whole position range information of the triangular surface element cut by the cutting surface;
based on the whole position range information, performing volume expansion on the first cutting cuboid and the second cutting cuboid, so that the ranges of the first cutting cuboid and the second cutting cuboid both comprise the whole position range;
and constructing a next-stage cuboid of the current-stage cuboid by utilizing the first cut cuboid and the second cut cuboid after the volume expansion.
S403, judging whether the next-stage cuboid meets a cutting stop condition, if not, determining that the next-stage cuboid is the current-stage cuboid, and executing the step S402;
if the cutting stopping condition is met, constructing and obtaining a target multi-stage cube segmentation model by utilizing cuboids of all stages;
the cutting stop condition refers to that the number of triangular surface elements included in the series of cuboid is lower than a preset surface element threshold value, or the number of triangular surface elements cut by the cut surface included in the series of cuboid is lower than a preset cutting threshold value.
The determining the corresponding cutting surface for the current cuboid comprises the following steps:
obtaining vertex coordinate information of all triangular surface elements included in the current cuboid;
carrying out variance calculation processing on the vertex coordinate information of all triangular surface elements according to coordinate axes to obtain coordinate variance values of the triangular surface elements on each coordinate axis of a target rectangular coordinate system;
determining the coordinate axis with the largest coordinate variance value as a cutting coordinate axis, and acquiring the vertex coordinate values of all triangular surface elements on the cutting coordinate axis;
constructing and obtaining a coordinate value sequence according to the value increasing sequence by utilizing the vertex coordinates;
calculating to obtain a median value of the coordinate value sequence, and determining a corresponding cutting surface by using a plane corresponding to the median value on the cutting coordinate axis as the current cuboid; for example, if the cutting coordinate axis is the z axis and the median value is 5, the cutting plane is z=5.
The determination of the cutting surface ensures that the quantity of triangular surface elements contained in two cuboid bodies obtained after cutting is the same as much as possible, thereby providing cutting efficiency.
And performing intersection discrimination processing on each virtual incident wave in the virtual incident wave set by using the target multi-stage cuboid segmentation model to obtain incident ray intersection point information, wherein the method comprises the following steps of:
s404, constructing and obtaining an intersecting judging incident wave set by utilizing virtual incident waves intersecting the first-stage cuboid in the virtual incident wave set; the intersection judging incident wave set comprises a plurality of incident waves; initializing reflection frequency information of all incident waves;
s405, judging whether the incident waves intersect with each level of cuboid according to the order of the cuboid in the target multi-level cuboid segmentation model from small to large, so as to obtain a least-level intersected cuboid;
s406, intersecting judgment is carried out on triangular surface elements contained in the cube intersected by the incident wave and the minimum level number, and triangular surface element information intersected with the incident wave is determined;
s407, processing the incident wave and the intersected triangular face element information vertex vector information to obtain intersection point information of the incident wave; the intersection point information of the incident wave comprises intersection point coordinate information, reflection times information, reflection direction vector information and reflection intensity information;
judging whether the reflection times and the reflection intensity of the incident wave meet the intersection point judging and stopping result or not; if the intersection judgment stopping result is met, judging whether intersection judgment processing is completed for all the incident waves of the intersection judgment incident wave set, if yes, executing S408, and if not, executing S405 for the next incident wave of the intersection judgment incident wave set;
if the intersection point judging and stopping result is not met, replacing the starting point information of the incident wave by utilizing intersection point coordinate information of the incident wave, respectively replacing the direction vector information and the radiation intensity information of the incident wave by utilizing the reflection direction vector information and the reflection intensity information of the incident wave to obtain updated incident wave, and executing S405;
s408, combining the intersection point information of all the incident waves in the virtual incident wave set to obtain the intersection point information of the incident rays;
the incident wave information comprises starting point information of the incident wave, direction vector information of the incident wave and radiation intensity information of the incident wave;
the intersection point judging and stopping result is that the reflection times of the current incident wave is larger than a preset reflection times threshold, and the reflection intensity of the current incident wave is smaller than a preset reflection intensity threshold;
the S406 and S407 may be implemented by S431, S432, S433.
The invention discloses a data processing device, which comprises:
a memory storing executable program code;
a processor coupled to the memory;
and the processor calls the executable program codes stored in the memory to execute the sea surface composite target far-zone electromagnetic scattering intensity estimation method.
The invention discloses a computer storage medium which stores computer instructions for executing the sea surface composite target far-zone electromagnetic scattering intensity estimation method when the computer instructions are called.
The foregoing is merely exemplary of the present application and is not intended to limit the present application. Various modifications and changes may be made to the present application by those skilled in the art. Any modifications, equivalent substitutions, improvements, etc. which are within the spirit and principles of the present application are intended to be included within the scope of the claims of the present application.
Claims (2)
1. The method for estimating the electromagnetic scattering intensity of the far-zone of the sea surface composite target is characterized by comprising the following steps of:
s1, acquiring target range position information and incident point source information; the incident point source information comprises incident point source position information and incident point source waveform information; the target range position information comprises target three-dimensional model information, target boundary position information, target center position information and target boundary point position information;
s2, determining a virtual incident wave set according to the target range position information and the incident point source information;
s3, performing subdivision and measurement processing on the target three-dimensional model to obtain a target triangular surface element model; the target triangular surface element model comprises target position measurement variance vector information, triangular surface element quantity information contained in a target, normal vector information of each triangular surface element, vertex coordinate information of each triangular surface element and medium information of each triangular surface element;
s4, processing the virtual incident wave set and the target triangular surface element model by using an intersection judgment model to obtain an intersection judgment result;
s5, processing the intersection discrimination result by using an offshore remote zone scattered field calculation model to obtain the electromagnetic scattering intensity information of the offshore composite target remote zone;
the determining a virtual incident wave set according to the target range position information and the incident point source information comprises the following steps:
s21, processing the target range position information and the incident point source position information by using a height coverage model to obtain a sea surface preset height h 0 ；
The expression of the height coverage model is as follows:
wherein, (x) s ,y s ,z s ) Is the position coordinate value of the incident point source in the target rectangular coordinate system, (x) a ,y a ,z a ) H is the position coordinate value of the geometric center point of the target in the rectangular coordinate system of the target 1 Altitude, which is the highest point of the target; the target rectangular coordinate system is a three-dimensional rectangular coordinate system established by taking the projection of a geometric center point of a target on the sea level as an origin and taking the sea level as an XOY axis;
s22, presetting the sea surface to be high h 0 Constructing a virtual horizontal plane;
s23, intersecting the virtual horizontal plane by utilizing a connecting line of the incident point source and the target boundary point to obtain a virtual intersection point;
s24, constructing a virtual parallelogram surrounding all virtual intersection points on the virtual horizontal plane;
s25, uniformly dividing the virtual parallelogram by utilizing the basic shape to obtain a virtual divided quadrilateral; the virtual segmentation quadrangle comprises a plurality of basic shapes;
s26, determining the incident point source as a starting point, determining the vertexes of the basic shapes contained in the virtual segmentation quadrangle as path points, and determining rays pointing from the starting point to the path points as virtual incident waves;
s27, constructing and obtaining a virtual incident wave set by utilizing all virtual incident waves;
the processing the virtual incident wave set and the target triangle surface element model by utilizing the intersection discrimination model to obtain an intersection discrimination result comprises the following steps:
s41, measuring variance vector information by using a direction vector and a target position of each virtual incident wave, and calculating to obtain a correction vector corresponding to the virtual incident wave; correction vector corresponding to the virtual incident waveThe calculation expression is as follows:
wherein,for the direction vector of the virtual incident wave, < >>The variance vector is measured for the target location,
s42, correcting the vertex coordinate information of the triangular surface element by utilizing the target position measurement variance vector information to obtain triangular surface element vertex vector information;
s43, processing the virtual incident wave and the triangular surface element vertex vector information by utilizing an intersection judgment model to obtain an intersection judgment result of the virtual incident wave;
s44, integrating all the intersection discrimination results of the virtual incident waves to obtain the intersection discrimination results of the virtual incident wave set; the intersection discrimination result of the virtual incident wave set comprises intersecting triangular face element information and intersection point coordinate information corresponding to each virtual incident wave;
the processing the virtual incident wave and the triangle primitive vertex vector information by utilizing an intersection discrimination model to obtain an intersection discrimination result of the virtual incident wave comprises the following steps:
s431, utilizing the incident point source vectorCorrection vector corresponding to virtual incident wave>And triangular surface element vertex vector information, calculating to obtain an intersection judgment vector set; the intersection judgment vector set comprises a first intersection judgment vector T 1 Second intersection judgment vector T 2 Third intersection judgment vector T 3 Fourth intersection judgment vector T 4 And a fifth intersection judgment vector T 5 The method comprises the steps of carrying out a first treatment on the surface of the The incident point source vector->
The calculation expression of the intersection judgment vector set is as follows:
T 5 ＝T 3 ×T 1 ；
s432, vector calculation processing is carried out on the intersecting judgment vector set to obtain a first intersecting judgment vector (p 1 ,p 2 ,p 3 ) The method comprises the steps of carrying out a first treatment on the surface of the The first intersection discrimination vector (p 1 ,p 2 ,p 3 ) The calculated expression of (2) is:
s433, judging the first intersection judgment vector to obtain an intersection judgment result of the virtual incident wave;
the step of judging the intersection judgment vector to obtain the intersection judgment result of the virtual incident wave comprises the following steps:
judging the first intersection judgment vector (p 1 ,p 2 ,p 3 ) Whether the intersecting condition is met or not, and a condition judgment result is obtained; the intersecting conditions include:
p 1 ≥0,p 2 ≥β 2 ,p 3 ≥β 3 ,1-p 2 -p 3 ≥β 1 ；
when the first intersection discrimination vector (p 1 ,p 2 ,p 3 ) When the intersecting condition is met, determining that the virtual incident wave intersects with the triangular surface element, and calculating to obtain an intersection point coordinate W of the virtual incident wave and the triangular surface element, wherein the calculation expression is as follows:
when the first intersection discrimination vector (p 1 ,p 2 ,p 3 ) When the intersection condition is not satisfied, determining that the virtual incident wave is not intersected with the triangular surface element;
the processing of the intersecting discrimination result by using the offshore remote region scattered field calculation model to obtain the electromagnetic scattering intensity information of the offshore composite target remote region comprises the following steps:
calculating the corresponding far-zone electromagnetic scattering intensity information of each virtual incident wave by using a marine far-zone scattering field calculation model;
superposing the electromagnetic scattering intensity of the far-zone of each virtual incident wave to obtain the electromagnetic scattering intensity information of the far-zone of the sea surface composite target;
the offshore remote region fringe field computation model has the expression:
wherein t represents a time variable, N represents an integrated electromagnetic constant vector, τ represents an incident pulse width of an incident point source, and Q 1 And Q 2 The first integral function and the second integral function are respectively expressed, and the expression is as follows:
wherein, [ L ] 1 ,L 2 ]When the triangle surface element is projected onto the XOY plane of the target rectangular coordinate system, the value range of the vertex of the triangle surface element on the X axis is also the upper limit and the lower limit of the integral range of the X axis; obtaining three intersection points of the three virtual incident waves and an XOY plane by using the virtual incident waves and two virtual incident waves closest to the virtual incident waves, and determining two corresponding linear equations by using two connecting lines between the three intersection points, wherein the expressions of the two linear equations are y=q respectively 1 x+w 1 And y=q 2 x+w 2 ，q 1 、w 1 、q 2 、w 2 Coefficients of two linear equations respectively; connecting the intersection point of the virtual incident wave and the triangular surface element with an integral point on the triangular surface element to obtain a first linear equation, wherein the expression of the first linear equation is y=qx+w, and q and w represent coefficients of the first linear equation; τ 1 Representing the time taken for the virtual incident wave to reach the integration point of the triangular surface element from the incident point source through the intersection point of the virtual incident wave and the triangular surface element, t 0 An incident pulse time width control coefficient of an incident point source; the first and second integral functionsIntegrating the variable t; the expression of the integral electromagnetic constant vector N is as follows:
wherein,an extra-unit normal vector representing a triangle bin, < >>A magnetic field unit vector, Z, in a virtual incident wave representing an incident point source 0 Is free space wave impedance, c is vacuum wave velocity, < ->Representing the incident point source vector +.>Is a modulus of (2);
and correcting the vertex coordinate information of the triangular surface element by utilizing the target position measurement variance vector information to obtain triangular surface element vertex vector information, wherein the calculation expression is as follows:
wherein,vectors corresponding to three vertex coordinates of the triangular surface element, respectively, < ->Three vertex vectors of the triangular surface element respectively.
2. A data processing apparatus, the apparatus comprising:
a memory storing executable program code;
a processor coupled to the memory;
the processor invokes the executable program code stored in the memory to perform the sea surface composite target far zone electromagnetic scattering intensity estimation method of claim 1.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310668652.1A CN116754847B (en) | 2023-06-07 | 2023-06-07 | Method and device for estimating electromagnetic scattering intensity of far-region of sea surface composite target |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310668652.1A CN116754847B (en) | 2023-06-07 | 2023-06-07 | Method and device for estimating electromagnetic scattering intensity of far-region of sea surface composite target |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116754847A CN116754847A (en) | 2023-09-15 |
| CN116754847B true CN116754847B (en) | 2024-01-23 |
Family
ID=87952585
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310668652.1A Active CN116754847B (en) | 2023-06-07 | 2023-06-07 | Method and device for estimating electromagnetic scattering intensity of far-region of sea surface composite target |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116754847B (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104656078A (en) * | 2013-11-25 | 2015-05-27 | 北京环境特性研究所 | Sea surface ship method based on ray tracing |
| CN110705058A (en) * | 2019-09-19 | 2020-01-17 | 电子科技大学 | Near-field electromagnetic scattering simulation method for ultra-electric large-scale target |
| KR102201256B1 (en) * | 2019-12-30 | 2021-01-12 | 인천대학교 산학협력단 | Measurement Method of Radar Electromagnetic Field |
| CN114167379A (en) * | 2021-11-19 | 2022-03-11 | 南京航空航天大学 | Complex target scattering characteristic rapid calculation method for SAR imaging simulation |
| CN115267720A (en) * | 2022-08-18 | 2022-11-01 | 中国人民解放军91977部队 | Method for calculating composite electromagnetic scattering RCS (radar cross section) of marine ship target |
-
2023
- 2023-06-07 CN CN202310668652.1A patent/CN116754847B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104656078A (en) * | 2013-11-25 | 2015-05-27 | 北京环境特性研究所 | Sea surface ship method based on ray tracing |
| CN110705058A (en) * | 2019-09-19 | 2020-01-17 | 电子科技大学 | Near-field electromagnetic scattering simulation method for ultra-electric large-scale target |
| KR102201256B1 (en) * | 2019-12-30 | 2021-01-12 | 인천대학교 산학협력단 | Measurement Method of Radar Electromagnetic Field |
| CN114167379A (en) * | 2021-11-19 | 2022-03-11 | 南京航空航天大学 | Complex target scattering characteristic rapid calculation method for SAR imaging simulation |
| CN115267720A (en) * | 2022-08-18 | 2022-11-01 | 中国人民解放军91977部队 | Method for calculating composite electromagnetic scattering RCS (radar cross section) of marine ship target |
Non-Patent Citations (2)
| Title |
|---|
| 海面上电大尺寸目标电磁散射特性研究;姬伟杰;刘平;关晓东;王炯;;现代导航(第05期);362-370 * |
| 海面环境下舰船目标的RCS测试建模分析;王雅芬 等;数字海洋与水下攻防;第第4卷卷(第第3期期);183-187 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116754847A (en) | 2023-09-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Von Herzen et al. | Accurate triangulations of deformed, intersecting surfaces | |
| JP5111638B2 (en) | Apparatus and method for dividing a parametric curve into smaller subpatches | |
| Sigg et al. | Signed distance transform using graphics hardware | |
| US20100231583A1 (en) | Image processing apparatus, method and program | |
| US20020082800A1 (en) | Method of composing three-dimensional multi-viewpoints data | |
| Yguel et al. | Efficient GPU-based construction of occupancy grids using several laser range-finders | |
| Chen et al. | Denoising of point cloud data for computer-aided design, engineering, and manufacturing | |
| US20240185503A1 (en) | Displacement-centric acceleration for ray tracing | |
| CN112154394A (en) | Terrain detection method, movable platform, control device, system and storage medium | |
| Meerits et al. | Real-time scene reconstruction and triangle mesh generation using multiple RGB-D cameras | |
| CN111540063A (en) | Full-automatic high-precision splicing method based on multi-station laser point cloud data | |
| CN110796735A (en) | NURBS curved surface finite element plate shell mesh division method and computer implementation system | |
| CN116754847B (en) | Method and device for estimating electromagnetic scattering intensity of far-region of sea surface composite target | |
| Sagawa et al. | Parallel processing of range data merging | |
| Maurya et al. | Performance of greedy triangulation algorithm on reconstruction of coastal dune surface | |
| CN117115239B (en) | Entrance ray intersection point acquisition method for remote electromagnetic scattering intensity estimation | |
| Bronsvoort et al. | Display of profiled sweep objects | |
| Zhu et al. | Triangulation of well-defined points as a constraint for reliable image matching | |
| Splietker et al. | Directional TSDF: Modeling surface orientation for coherent meshes | |
| Chen et al. | A point-based offsetting method of polygonal meshes | |
| Wiemann et al. | An evaluation of open source surface reconstruction software for robotic applications | |
| CN103761416B (en) | Method for rapidly forecasting scattering sound field of underwater triangle corner reflector | |
| Sipos et al. | Creating good quality meshes from smooth implicit surfaces | |
| Teutsch et al. | Evaluation and correction of laser-scanned point clouds | |
| CN118334263B (en) | High-precision modeling method for fusion laser point cloud based on truncated symbol distance function |
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 |