A radar monitoring data plane domain deformation determination method, storage medium and system

Abstract

The present invention provides a method for determining surface deformation of radar monitoring data, which uses a dichotomy method to eliminate noise data, sets different weights for the eliminated data according to image intensity, and finally calculates a comprehensive deformation value based on the cumulative deformation value of each cell and the weight value, and regards the comprehensive deformation value as the cumulative deformation value of the surface. The method for determining surface deformation of radar monitoring data of the present invention utilizes all deformation variables in the surface and radar reflection image intensity data so that the determined deformation value can fully reflect the overall deformation state of the surface. The present invention also provides a storage medium and a radar monitoring system, which use the method for determining surface deformation of radar monitoring data for monitoring, and the storage medium is used to store the method for determining surface deformation of radar monitoring data.

Description

Technical Field

[0001] The present invention relates to the field of surveying, mapping and safety monitoring technology, and in particular to a method for determining surface deformation of radar monitoring data. Background Art

[0002] Arc-type synthetic aperture deformation monitoring radars are capable of large-scale deformation monitoring. Conventional deformation analysis methods select one or several PS points for deformation data analysis. However, this approach remains within the framework of conventional GNSS point monitoring and fails to fully exploit the advantages of radar surface monitoring. Furthermore, due to the inherent characteristics of radar, individual points may contain singularities due to factors such as atmospheric errors, making the data unreliable. Therefore, selecting a surface area for radar data and determining deformation trends within that area can better leverage the advantages of radar.

[0003] However, considering that the deformation variables of a cell in a single scan of a radar surface range from a dozen to hundreds or even thousands, it is difficult to fully utilize them during data analysis. Therefore, it is necessary to determine a deformation variable to represent the degree of deformation of the surface during a single scan.

[0004] In summary, there is an urgent need for a method for determining surface deformation of radar monitoring data to solve the problems existing in the existing technology. Summary of the Invention

[0005] The present invention aims to provide a method for determining surface deformation using radar monitoring data, aiming to leverage the advantages of radar surface monitoring and implement surface deformation processing and analysis. The specific technical solution is as follows:

[0006] A method for determining surface deformation of radar monitoring data comprises the following steps:

[0007] Step S1: Select the analysis area, the cell set covered by the area is U{u1, u2, u3, ...u n}, the cumulative deformation value set corresponding to the cell set covered by the face region is D{d1, d2, d3, ...d n}, the image intensity value set corresponding to the cell set covered by the area is M{m1, m2, m3, ...m n}, where n is a natural number greater than 0, and n is less than or equal to the total number of cells within the radar monitoring range;

[0008] Step S2: Count the minimum value d in the set D min and the maximum value d max , using the dichotomy method to obtain the screening index d′ cmin and d′ cmax ;

[0009] Step S3, obtain the cumulative deformation value in [d′cmin , d′ cmax ] between the cell set U′{u′1,u′2,u′3,…u′ k}, the cumulative deformation value set corresponding to the set U′ is D c {d′1, d′2, d′3,...d′ k}, the corresponding image intensity value set is M′{m′1, m′2, m′3, ...m ′ k}, where k is a natural number and is less than n;

[0010] Step S4: Set q weight levels according to the image intensity values ​​in the set M′. The cells in the set U′ are divided into corresponding weight levels according to their image intensity values. The cumulative deformation value set of the cell set in the jth weight level is Where j is less than or equal to a;

[0011] Step S5: Calculate the comprehensive deformation value D Cv , the comprehensive deformation value is regarded as the cumulative deformation value of the surface area:

[0012]

[0013] in, Representative Set The number of elements in Representative Set The cumulative deformation value of the i-th cell in Q j and Q m are the weight values ​​of the j-th and m-th weight levels respectively.

[0014] In the above technical solution, preferably, in step S2, the screening index d' is obtained cmin The specific steps are as follows:

[0015] Step A1: Set the deformation value index d c ,shilling

[0016] Step A2: For each value in set D, when d i ≥d c When d i Add to set D g ; Otherwise, put it into set D s ; where 1≤i≤n;

[0017] Step A3: Statistical set D g The number of data in C Dg and set D s The number of data in C Ds ;

[0018] Step A4, when d' c = d c , recalculate where P is a set threshold value, 0

[0019] Step A5, repeat steps A2-A4 until the first time d c is assigned to d' cmin . In the above technical solution, preferably, in step S2, the specific steps of obtaining the screening index d' cmax are as follows:

[0020] Step B1, set the deformation value index d c , first let

[0021] Step B2, for each value in set D, when d i ≥ d c , then put d i into set D g ; otherwise, put it into set D s ; where 1≤i≤n;

[0022] Step B3, count the number of data C g in set D Dg and the number of data C s in set D Ds ;

[0023] Step B4, when d' c = d c , recalculate where P is a set threshold value, 0

[0024] Step B5, repeat steps B2-B4 until the first time d c is assigned to d' cmax .

[0025] In the above technical solution, preferably, -50≤ the image intensity value m of the surface domain covering unit cell <0; weight levels Q1, Q2, Q3, Q4 and Q5 are set; the image intensity value range corresponding to each weight level is as follows:

[0026] -50≤m<40, Q1=1;

[0027] -40≤m<30, Q2=2;

[0028] -30≤m<20, Q3=3;

[0029] -20≤m<10, Q4=4;

[0030] -10≤m<0, Q5=5.

[0031] Preferably in the above technical solution, the analysis surface domain comprises a plurality of continuous cells.

[0032] The application further provides a storage medium for storing the radar monitoring data surface domain deformation determination method.

[0033] The application further provides a radar monitoring system, which adopts the radar monitoring data surface domain deformation determination method.

[0034] The application has the following beneficial effects:

[0035] The application provides a radar monitoring data surface domain deformation determination method, which first screens out a part of maximum and minimum data, because considering the influence of environmental factors on radar electromagnetic waves, the possibility of the part of data having larger noise is relatively large.Secondly, different weights are set according to different image intensities, and the signal quality and calculation data reliability are considered to be higher when the intensity of electromagnetic wave reflection is stronger.Therefore, the comprehensive deformation value D Cv calculated by the method of the application can fully represent the deformation condition of the surface domain.

[0036] The radar monitoring surface domain deformation determination method of the application uses all deformation variables and radar reflection image intensity data in the surface domain, so that the determined deformation value can fully reflect the deformation state of the whole surface domain.

[0037] In addition to the purposes, features and advantages described above, the application has other purposes, features and advantages. The application will be further described below with reference to the drawings. BRIEF DESCRIPTION OF DRAWINGS

[0038] The accompanying drawings, which form a part of this application, are included to provide a further understanding of the application, and are incorporated herein for explanation by illustrating a preferred embodiment of the present application and to explain the principles of the application, and are not intended to limit the application. In the drawings:

[0039] Figure 1 is a schematic diagram of radar monitoring range monitoring. DETAILED DESCRIPTION

[0040] In order to facilitate the understanding of the application, the application will be described more fully below, and the preferred embodiments of the application are given. However, the application can be realized in many different forms, and is not limited to the embodiments described herein. On the contrary, the purpose of providing these embodiments is to make the disclosure of the application more thorough and comprehensive.

[0041] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application belongs. The terminology used in the description herein is for describing particular embodiments only and is not intended to be limiting of the application.

[0042] Embodiment 1:

[0043] The embodiment provides a radar monitoring data surface domain deformation determination method.

[0044] The circular arc synthetic aperture deformation monitoring radar can realize large range deformation monitoring. The monitoring range is a sector. The radar scanning result area is divided into a plurality of small sector cells by the resolution of the distance direction and the angle direction. The monitoring result is the cumulative deformation value of the small sector cells relative to the initial monitoring time. When the surface domain deformation is analyzed, the surface domain contains a plurality of continuous cells.

[0045] The determination method comprises the following steps:

[0046] Step S1, selecting an analysis surface domain, the cell set covered by the surface domain is U {u1, u2, u3,... u n}, the cumulative deformation value set corresponding to the cell set covered by the surface domain is D {d1, d2, d3,... d n}, and the image intensity value set corresponding to the cell set covered by the surface domain is M {m1, m2, m3,... m n}, wherein n is a natural number greater than 0, and n is less than or equal to the total number of cells in the radar monitoring range.

[0047] For the data of a scan, the cumulative deformation value of the cell u i is d i , and the image intensity value is m i ; that is, the elements in the sets U, D and M are in one-to-one correspondence.

[0048] Step S2, counting the minimum value d min and the maximum value d max in the set D, and obtaining the screening indexes d' cmin and d' cmax by using the dichotomy method.

[0049] Further, in step S2, the specific steps of obtaining the screening index d' cmin are as follows:

[0050] Step A1, setting a deformation value index d c , and first setting

[0051] Step A2, for each value in set D, when d i ≥ d c , then put d i into set D g ; otherwise, put d s into set D g ; where 1≤i≤n;

[0052] Step A3, count the number of data C Dg in set D s and the number of data C Ds in set D c ;

[0053] Step A4, when , let d' c = d ， recalculate where P is a set threshold, 0

[0054] Step A5, repeat steps A2-A4 until the first time to meet assign d c to d' cmin .

[0055] Further, in step S2, the specific steps of obtaining the screening index d' cmax are as follows:

[0056] Step B1, set the deformation value index d c , first let

[0057] Step B2, for each value in set D, when d i ≥ d c , then put d i into set D g ; otherwise, put d s into set D g ; where 1≤i≤n;

[0058] Step B3, count the number of data C Dg in set D s and the number of data C Ds in set D c ;

[0059] Step B4, when , let d' c = d c , recalculate where P is a set threshold, 0

[0060] Step B5, repeat steps B2-B4 until the first time to meet At this time, d c Assign to d′ cmax .

[0061] Step S3: Obtain the cumulative deformation value in the set U in [d′ cmin , d′ cmax ] between the cell set U′{u′1,u′2,u′3,…u′ k}, the cumulative deformation value set corresponding to the set U′ is D c {d′1, d′2, d′3,...d′ k}, the corresponding image intensity value set is M′{m′1, m′2, m′3, ...m′ k}, where k is a natural number and is less than n; those skilled in the art will appreciate that U′, D c and M′ are in one-to-one correspondence.

[0062] Step S4: Set a weight level according to the image intensity value in the set M′. The cells in the set U′ are divided into corresponding weight levels according to their image intensity values. The cumulative deformation value set of the cell set in the jth weight level is D j c , where j is less than or equal to a;

[0063] Step S5: Calculate the comprehensive deformation value D Cv , the comprehensive deformation value is regarded as the cumulative deformation value of the surface area:

[0064]

[0065] in, Representative Set The number of elements in Representative Set The cumulative deformation value of the i-th cell in Q j and Q m are the weight values ​​of the j-th and m-th weight levels respectively.

[0066] Further preferably, since generally -50≤the image intensity value of the area coverage cell m<0, this embodiment sets weight levels Q1, Q2, Q3, Q4 and Q5; the image intensity value range corresponding to each weight level is as follows:

[0067] When -50≤m<40, Q1=1;

[0068] When -40≤m<30, Q2=2;

[0069] When -30≤m<20, Q3=3;

[0070] -20≤m<10, Q4=4;

[0071] -10≤m<0, Q5=5.

[0072] For the cells in the set U', according to the above classification method, the cumulative deformation value set of the cells with image intensity m∈[-50, 40) is The cumulative deformation value set of the cells with image intensity m∈[-40, 30) is The cumulative deformation value set of the cells with image intensity m∈[-30, 20) is The cumulative deformation value set of the cells with image intensity m∈[-20, 10) is The cumulative deformation value set of the cells with image intensity m∈[-10, 0) is

[0073] The embodiment also provides a radar monitoring system which adopts the radar monitoring data surface domain deformation determination method.

[0074] In order to verify the effectiveness of the method of the embodiment, the method of the embodiment is compared with the existing conventional method.

[0075] Select a certain surface domain coverage area, the grid is 12*27, a total of 324 grids, and the cumulative deformation value set D formed by the grids is:

[0076] {18.19，11.23，5.64，4.09，1.47，1.41，3.84，8.00，18.75，29.50，38.19，45.41，22.91，14.49，7.11，5.14，1.53，1.63，4.45，8.59，19.23，27.84，35.16，41.41，23.41，16.14，8.30，6.04，2.38，2.57，5.06，9.07，16.72，22.75，27.27，31.36，29.67，23.14,15.20,13.09,8.03,7.11，8.71，10.74，14.70,17.72,19.84，21.52，37.03，34.47，30.00,28.30，23.36,20.41,21.63，18.06,16.86，17.27,15.69,14.03，38.47,39.38，38.72,38.28,32.06，28.03,25.09,18.28,13.65,11.98,8.59,5.54,49.56,52.13,54.09,53.88,47.41,41.78,35.47,25.41,16.61，13.01,7.46,5.00,61.13,64.00,66.31,65.25,57.69,50.38,44.97,35.84,28.52,23.58,17.25,14.29,69.31,73.38,76.81,75.38,68.63,62.31,56.38,47.16,39.44,32.31,25.34,22.66,74.13,76.44,78.56,77.25，71.81,65.50,56.28,46.28,37.66,26.64,20.25,18.45,75.13,75.25,76.50,76.50,75.06，72.75,66.75,59.13,49.25,36.22,26.84,22.20,69.81,71.00,73.94,75.50,74.56,72.69,68.25,62.28,54.31,42.56，33.19，27.81，66.56，68.06，72.13，73.63，73.50，73.13，67.69，61.41，53.84，43.16，35.38，31.94,61.16,64.63，67.94，69.13，69.19，68.69,64.81，60.59，56.47,48.66，40.09,35.56，54.09，59.38，62.81, 63.81, 63.56, 63.81, 62.34, 61.47, 59.38, 55.09, 45.88, 41.03, 57.84, 62.13, 64.06, 63.59, 61.94, 60.38, 58.44, 59.06, 60.16, 57.47, 50.19, 49.09, 63.13, 67.06, 66.25, 63.78, 62.19, 60.56, 58.88, 60.72, 61.44, 59.88, 55.13, 54.75, 66.88, 72.50, 71.00, 68.81, 66.63, 62.81, 58.59, 57.16, 55.22, 54.13, 51.03, 51.56, 64.00, 67.81, 65.44, 62.59, 58.78, 55.13, 52.00, 50.44, 46.56, 45.72, 43.81, 44.94, 60.09, 61.91, 60.16, 57.66, 53.72, 50.75, 49.06, 49.75, 49.41, 51.16, 51.47, 54.16, 57.22, 61.75, 59.91, 56.78, 53.63, 51.13, 48.50, 49.06, 48.84, 52.34, 55.06, 59.03, 52.94, 55.00, 52.59, 49.91, 45.69, 43.16, 43.84, 47.34, 49.94, 55.69, 59.03, 61.31, 49.88, 49.19, 46.09, 42.63, 39.59, 42.50, 47.72, 54.88, 59.22, 64.94, 65.38, 63.34, 58.59, 57.00, 54.38, 51.31, 48.03, 50.00, 53.03, 58.53, 62.09, 67.75, 68.88, 66.75, 65.44, 64.13, 62.28, 58.97, 54.16, 52.50, 51.91, 54.94, 57.84, 63.50, 67.50, 68.13, 67.69, 64.06, 63.50, 61.84, 57.56, 57.09, 57.41, 59.09, 60.53, 63.88, 66.06, 65.75, 70.69, 67.69, 68.50, 66.69, 63.81, 65.69, 64.88, 62.16, 60.50, 60.72, 60.63, 60.81.

[0077] The corresponding set of image intensity data M is:

[0078] {-37.16,-41.44,-51.02,-46.25,-49.70,-45.62,-39.81,-40.94,-44.31,-41.40,-42.71,-44.12,-33.99,-37.30,-49.96,-38.73,-40.59,-49.67,-38.58,-41.61,-39.81,-36.44,-36.98,-35.80,-35.61,-39.87,-48.54,-42.26,-47.84,-41.63,-42.29,-56.42,-41.98,-37.05,-32.49,-32.04,-43.72,-41.48,-40.29,-39.33,-36.42,-37.62,-49.16,-46.38,-47.69,-38.71,-35.33,-37.84,-30.16,-30.93,-32.80,-35.18,-37.24,-40.69,-42.05,-44.92,-47.87,-41.25,-45.98,-46.85,-41.60,-37.07,-35.81,-36.34,-35.72,-36.30,-42.86,-43.97,-37.35,-38.91,-43.87,-43.03,-20.17,-21.02,-22.93,-25.57,-27.29,-29.49,-33.41,-34.13,-36.29,-43.77,-36.65,-38.57,-20.87,-23.44,-27.26,-31.56,-32.39,-33.82,-38.47,-44.55,-51.94,-38.71,-37.41,-43.72,-21.84,-24.65,-28.19,-32.74,-38.04,-41.18,-40.45,-39.38,-39.28,-38.65,-40.80,-46.52,-24.57,-27.63,-30.40,-31.22,-29.99,-28.29,-27.80,-29.50,-31.80,-32.24,-34.29,-36.28,-30.61,-31.99,-31.26,-30.06,-29.62,-28.99,-29.58,-33.28,-34.72,-33.74,-37.14,-35.35,-26.15,-25.35,-25.08,-25.69,-27.94,-32.14,-36.72,-41.82,-54.68,-53.41,-44.88,-41.94,-33.38,-29.36,-27.67,-27.88,-29.49,-32.97,-37.09,-37.00,-38.90,-42.62,-41.14,-42.59,-19.74,-21.42,-25.22,-34.67,-38.01,-30.05,-27.79,-27.63,-29.54,-31.17,-31.22,-34.85,-22.36,-25.82,-31.64,-30.87,-26.93,-24.84,-23.60,-23.27,-24.00,-24.59,-25.11,-28.08,-32.03,-28.09,-26.19,-25.60,-25.90,-26.32,-26.15,-25.89,-26.59,-27.97,-28.66,-30.32,-30.88,-29.39,-29.32,-30.56,-31.96,-32.36,-32.72,-34.27,-35.67,-37.26,-43.93,-42.48,-16.27,-14.15,-13.12,-13.10,-13.98,-15.51,-17.37,-19.71,-23.25,-27.52,-29.09,-31.30,-14.39,-12.33,-11.21,-11.07,-11.90,-13.61,-16.14,-19.67,-25.27,-33.43,-29.68,-27.25,-19.79,-17.83,-16.78,-17.03,-18.41,-20.35,-21.82,-22.41,-23.12,-24.64,-26.62,-28.29,-26.18,-25.97,-21.27,-18.82,-17.63,-17.13,-17.36,-18.78,-22.18,-29.19,-36.02,-29.70,-32.63,-32.12,-24.84,-21.20,-19.28,-18.60,-18.97,-20.36,-22.58,-24.85,-26.19,-26.79,-49.75,-47.23,-38.81,-31.49,-26.01,-22.56,-20.62,-19.73,-19.46,-19.60,-20.21,-21.21,-25.23,-22.89,-23.19,-26.47,-30.05,-25.64,-22.18,-20.68,-20.55,-21.48,-23.02,-24.33,-21.77,-21.56,-23.20,-26.61,-32.00,-41.11,-43.19,-35.26,-31.90,-29.69,-27.86,-26.39,-21.75,-24.47,-31.20,-43.41,-30.16,-27.68,-28.81,-32.76,-37.38,-36.01,-32.77,-29.93,-23.62,-28.81,-33.21,-23.77,-19.51,-17.57,-17.05,-17.46,-18.36,-19.58,-21.40,-24.37}..

[0079] Using the method of the embodiment, 5 weight levels are set, wherein the cumulative deformation value set of the cell with image intensity m∈[-50, 40) is the cumulative deformation value set of the cell with image intensity m∈[-40, 30) is the cumulative deformation value set of the cell with image intensity m∈[-30, 20) is the cumulative deformation value set of the cell with image intensity m∈[-20, 10) is the cumulative deformation value set of the cell with image intensity m∈[-10, 0) is

[0080] The final calculated comprehensive deformation value D Cv = 52.31.

[0081] Using the conventional average method, the average value of the cumulative deformation values of all cells is 47.56, and the deformation difference is 4.75 mm. According to the data in set D and set M, the cells with relatively large deformation account for the majority in the surface domain analyzed this time, and therefore, when selecting a value to represent the cumulative deformation value of the surface domain, the comprehensive deformation value 52.31 is more appropriate than the average value 47.56.

[0082] The above only describes the preferred embodiments of the present application and is not used to limit the present application. For those skilled in the art, the present application can have various modifications and changes. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application shall be included in the protection scope of the present application.

Claims

1. A method for determining surface deformation of radar monitoring data, characterized in that: The following steps are involved: Step S1: Select the analysis area, the cell set covered by the area is U{u1, u2, u3, ...u n }, the cumulative deformation value set corresponding to the cell set covered by the face region is D{d1, d2, d3, ...d n }, the image intensity value set corresponding to the cell set covered by the area is M{m1, m2, m3, ...m n }, where n is a natural number greater than 0, and n is less than or equal to the total number of cells within the radar monitoring range; Step S2: Count the minimum value d in the set D min and the maximum value d max , using the dichotomy method to obtain the screening index d′ cmin and d′ cmqx ; Step S3, obtain the cumulative deformation value in [d′ cmin , d′ cmax ] between the cell set U′{u′1,u′2,u′3,…u′ k }, the cumulative deformation value set corresponding to the set U′ is D c {d′1, d′2, d′3,...d′ k }, the corresponding image intensity value set is M′(m′1, m′2, m′3, ...m′ k }, where k is a natural number and is less than n; Step S4: Set a weight level according to the image intensity value in the set M′. The cells in the set U′ are divided into corresponding weight levels according to their image intensity values. The cumulative deformation value set of the cell set in the jth weight level is Where j is less than or equal to a; Step S5: Calculate the comprehensive deformation value D Cv , the comprehensive deformation value is regarded as the cumulative deformation value of the surface area: in, Representative Set The number of elements in Representative Set The cumulative deformation value of the i-th cell in Q j and Q m are the weight values ​​of the j-th and m-th weight levels respectively.

2. The method for determining surface deformation based on radar monitoring data according to claim 1, characterized in that: In step S2, the screening index d' is obtained cmin The specific steps are as follows: Step A1: Set the deformation value index d c ,shilling Step A2: For each value in set D, when d i ≥d c When d i Put into set D g ; Otherwise, put it into set D s ; where 1≤i≤n; Step A3: Statistical set D g The number of data in C Dg and set D s The number of data in C Ds ; Step A4: When d′ c =d c , recalculate Where P is the set threshold, 0<P<1: Step A5: Repeat steps A2-A4 until the first At this time, d c Assign to d′ cmin .

3. The method for determining surface deformation based on radar monitoring data according to claim 1, characterized in that: In step S2, the screening index d' is obtained cmax The specific steps are as follows: Step B1: Set the deformation value index d c ,shilling Step B2: For each value in set D, when d i ≥d c When d i Put into set D g ; Otherwise, put it into set D s ; where 1≤i≤n; Step B3: Statistical set D g The number of data in C Dg and set D s The number of data in C Ds ; Step B4: When d′ c =d c , recalculate Where P is the set threshold, 0<P<1: Step B5: Repeat steps B2-B4 until the first At this time, d c Assign to d′ cmax .

4. The method for determining surface deformation based on radar monitoring data according to claim 1, wherein: -50≤The image intensity value of the area coverage cell m<0; Set the weight levels Q1, Q2, Q3, Q4 and Q5; The image intensity value range corresponding to each weight level is as follows: When -50≤m<40, Q1=1; When -40≤m<30, Q2=2; When -30≤m<20, Q3=3; When -20≤m<10, Q4=4; When -10≤m<0, Q5=5.

5. The method for determining surface deformation based on radar monitoring data according to claim 1, characterized in that: The analysis area consists of several continuous cells.

6. A storage medium, characterized in that The storage medium is used to store the radar monitoring data surface deformation determination method according to any one of claims 1 to 5.

7. A radar monitoring system, characterized in that: The system adopts the method for determining surface deformation of radar monitoring data as described in any one of claims 1 to 5.

Images