Engineering construction abnormal feature analysis method based on multi-angle monitoring video image

Abstract

The invention provides an engineering construction abnormal feature analysis method based on a multi-angle monitoring video image. The method comprises the following steps: S1, entering a natural resource information system; s2, acquiring image data stored in a natural resource information system; and S3, analyzing and determining the abnormal change condition of engineering construction according to the acquired image data.

Description

Technical Field

[0001] The present invention relates to the technical field of engineering construction monitoring, and particularly to a method for analyzing abnormal characteristics of engineering construction based on multi-angle monitoring video images. Background Art

[0002] In the field of monitoring abnormal changes in engineering construction, traditional methods mainly rely on manual on-site shooting and collection, which have problems such as low monitoring efficiency, high false alarm rate, and poor dynamic adaptability. With the integrated application of geographic information system (GIS), remote sensing (RS) and global positioning system (GPS) technologies (i.e., 3S technology), it is bound to change the traditional way and improve aspects such as monitoring efficiency, accuracy, and dynamic adaptability. Summary of the Invention

[0003] The present invention aims to at least solve the technical problems existing in the prior art, and particularly innovatively proposes a method for analyzing abnormal characteristics of engineering construction based on multi-angle monitoring video images.

[0004] In order to achieve the above object of the present invention, the present invention provides a method for analyzing abnormal characteristics of engineering construction based on multi-angle monitoring video images, including the following steps:

[0005] S1, enter the natural resource information system;

[0006] S2, obtain the image data stored in the natural resource information system;

[0007] S3, determine the abnormal change situation of the engineering construction according to the obtained image data analysis.

[0008] In a preferred embodiment of the present invention, the method for determining the abnormal change situation of the engineering construction according to the engineering construction monitoring images stored in the natural resource information system in step S3 includes the following steps:

[0009] S31, obtain the selected date, and the selected date includes the selection of the reference date and the comparison date; the reference date includes year, month, and day, the comparison date includes year, month, and day, and if the reference date and the comparison date are not the same day, they are different dates;

[0010] S32, obtain the search string stored under the selected date according to the selected date;

[0011] S33, search for the stored image corresponding to the search string according to the search string, and use the search string and the stored image corresponding to the search string to obtain the engineering construction monitoring images of different dates;

[0012] S34, determine the abnormal change situation of the engineering construction according to the engineering construction monitoring images of different dates.

[0013] In a preferred embodiment of the present invention, the method for determining the abnormal changes in the engineering construction according to the engineering construction monitoring images on different dates in step S34 includes the following steps:

[0014] S341, obtain the pixel values of each channel in the engineering construction monitoring image:

[0015]

[0016] E j represents the pixel value of the j-th channel, where j = r, g, b;

[0017] (x″ p ,y″ q ) j represents the pixel value of the j-th channel at the position (x″ p ,y″ q ),

[0018] p = 1, 2, 3,..., P, where P is the number of horizontal pixel points;

[0019] q = 1, 2, 3,..., Q, where Q is the number of vertical pixel points;

[0020] S342, generate the pixel value of the merged channel according to the pixel values of each channel in the obtained engineering construction monitoring image:

[0021]

[0022] where, E rgb represents the pixel value of the merged channel;

[0023] δ r represents the pixel value coefficient of the red channel;

[0024] δ g represents the pixel value coefficient of the green channel;

[0025] δ b represents the pixel value coefficient of the blue channel;

[0026] δ r (x″ p ,y″ q ) r +δ g (x″ p ,y″ q ) g +δ r (x″ p ,y″ q ) b represents the pixel value of the merged channel at the position (x″ p ,y″ q );

[0027] p = 1, 2, 3, ……, P, where P is the number of horizontal pixel points;

[0028] q = 1, 2, 3, ……, Q, where Q is the number of vertical pixel points;

[0029] S343. Determine the pixel value for pixel values on the same date:

[0030] S3431. Arrange all the pixel values at the same pixel coordinate in ascending order, which are Pix1, Pix2, Pix3, ……, Pix H ;

[0031] S3432. Take as the pixel value participating in the calculation to obtain the first determined pixel value:

[0032]

[0033] where Pix1 represents the pixel value ranked 1st after arranging the pixel values at the same pixel coordinate in ascending order;

[0034] Pix2 represents the pixel value ranked 2nd after arranging the pixel values at the same pixel coordinate in ascending order;

[0035] Pix3 represents the pixel value ranked 3rd after arranging the pixel values at the same pixel coordinate in ascending order;

[0036] Pix H represents the pixel value ranked Hth after arranging the pixel values at the same pixel coordinate in ascending order;

[0037] represents the pixel value ranked h1th after arranging the pixel values at the same pixel coordinate in ascending order;

[0038] represents the pixel value ranked h2th after arranging the pixel values at the same pixel coordinate in ascending order;

[0039] Pix con1 represents the first determined pixel value;

[0040] H represents the number of images taken on the same date;

[0041]

[0042] S3433. Take all the pixel values Pix1, Pix2, Pix3, ……, Pix at the same pixel coordinate HAs the pixel value participating in the calculation, obtain the second determined pixel value:

[0043]

[0044] Among them, Pix con2 represents the second determined pixel value;

[0045] H represents the number of images taken on the same date;

[0046] Pix j represents the pixel value at the j-th position after arranging the pixel values of the same pixel coordinates in ascending order;

[0047] S3434, determine the size relationship between the first determined pixel value and the second determined pixel value:

[0048] If |Pix con1 - Pix con2 | ≤ Pix sd , then the second determined pixel value is the determined pixel value;

[0049] If |Pix con1 - Pix con2 | > Pix sd , then the first determined pixel value is the determined pixel value; Pix sd is the preset pixel value exceeding the threshold;

[0050] S344, perform pixel value abnormality determination on the pixel values of different dates:

[0051] If then there is no abnormal change in the engineering construction at the pixel point (x″′ p , y″′ q ); is the preset first abnormal change threshold,

[0052] If then there is an abnormal change in the engineering construction at the pixel point (x″′ p , y″′ q );

[0053] Among them, then there is a level III abnormal change in the engineering construction at the pixel point (x″′ p , y″′ q ); is the preset second abnormal change threshold;

[0054] then at the pixel point (x″′ p , y″′ q) there is a Class II abnormal change in the engineering construction at is the preset third abnormal change threshold;

[0055] then there is a Class I abnormal change in the engineering construction at the pixel point (x″′ p , y″′ q );

[0056] p = 1, 2, 3, ……, P, where P is the number of horizontal pixel points;

[0057] q = 1, 2, 3, ……, Q, where Q is the number of vertical pixel points; Pix dat (x″′ p , y″′ q ) is the pixel value of the reference date image at the pixel point (x″′ p , y″′ q ); Pix com (x″′ p , y″′ q ) is the pixel value of the comparison date image at the pixel point (x″′ p , y″′ q );

[0058] The Class III abnormal change is lower than the Class II abnormal change, and the Class I abnormal change is higher than the Class II abnormal change.

[0059] In a preferred embodiment of the present invention, the method for entering the natural resource information system in step S1 includes the following steps:

[0060] S11, input the account number and password on the natural resource information system login interface; the number of digits of the account number and password is not less than 6 digits;

[0061] S12, after inputting the account number and password on the natural resource information system login interface, perform security encryption processing on the input account number and password;

[0062] S13, after performing security encryption processing on the input account number and password, send the securely encrypted account number and password to the natural resource information system verification platform;

[0063] S14, after the natural resource information system verification platform receives the securely encrypted account number and password, decrypt the received securely encrypted account number and password;

[0064] S15, if the decrypted account number and password both exist in the natural resource information system verification platform, then enter the natural resource information system;

[0065] Otherwise, if the account number and password entered on the natural resource information system login interface are incorrect, the correct account number and password need to be re-entered.

[0066] In a preferred embodiment of the present invention, the method for securely encrypting the input account number and password in step S12 includes the following steps:

[0067] S121, obtain the current moment, convert the current moment into a ten-digit standard format, and the current moment includes year, month, day, and hour;

[0068] S122, determine the size relationship between the number of digits of the input account number and password and 10:

[0069] S1221, if the number of digits of the account number is less than 10, then remove e1 digits from the front of the ten-digit standard format moment, and the number of digits after removal is equal to the number of digits of the account number; perform the next step; e1 = 10 - e0, e0 is the number of digits of the account number, and e1 is the number of digits removed from the moment;

[0070] If the number of digits of the account number is greater than 10, then add e2 digits e3 to the front of the ten-digit standard format moment, and the number of digits after addition is equal to the number of digits of the account number; perform the next step; e2 = e0 - 10, e0 is the number of digits of the account number, e2 is the number of digits added to the moment, and e3 is the last digit value in the ten-digit standard format;

[0071] If the number of digits of the account number is equal to 10, perform the next step;

[0072] S1222, if the number of digits of the password is less than 10, then remove f1 digits from the front of the ten-digit standard format moment, and the number of digits after removal is equal to the number of digits of the password; perform the next step; f1 = 10 - f0, f0 is the number of digits of the password, and f1 is the number of digits removed from the moment;

[0073] If the number of digits of the password is greater than 10, then add f2 digits f3 to the front of the ten-digit standard format moment, and the number of digits after addition is equal to the number of digits of the password; perform the next step; f2 = f0 - 10, f0 is the number of digits of the password, f2 is the number of digits added to the moment, and f3 is the last digit value in the ten-digit standard format;

[0074] If the number of digits of the password is equal to 10, perform the next step;

[0075] At this time, the account number and password remain unchanged, but the moments corresponding to the account number and password may change;

[0076] S123, shift the first digits of the account number and password backward by f4 digits respectively, f4 = floor(f3 / 2), and f4 is the number of digits shifted backward;

[0077] S124. Convert each character in the account, password, and time into 7-bit binary ASCII form. Calculate the ASCII-formatted account with the ASCII-formatted time, and calculate the ASCII-formatted password with the ASCII-formatted time to obtain the encrypted account and the encrypted password respectively.

[0078]

[0079] g′2(φ) represents the φ-th bit value in the ASCII-formatted encrypted account;

[0080] g2(φ) represents the φ-th bit value in the ASCII-formatted account;

[0081] t2(φ) represents the φ-th bit value in the ASCII-formatted time;

[0082] φ = 1, 2, 3, ……, φ1; φ1 is the number of bits of the ASCII-formatted account;

[0083] φ1 = 7φ2, where φ2 is the number of bits of the account;

[0084] ◇ represents the swap value operator;

[0085]

[0086] v′2(θ) represents the θ-th bit value in the ASCII-formatted encrypted password;

[0087] v2(θ) represents the θ-th bit value in the ASCII-formatted password;

[0088] t′2(θ) represents the θ-th bit value in the ASCII-formatted time;

[0089] θ = 1, 2, 3, ……, θ1; θ1 is the number of bits of the ASCII-formatted password;

[0090] θ1 = 7θ2, where θ2 is the number of bits of the password;

[0091] Convert the ASCII-formatted encrypted account and the ASCII-formatted encrypted password into character-formatted encrypted account and encrypted password.

[0092] The present invention also discloses a computer system, including:

[0093] A processor;

[0094] A memory for storing processor-executable instructions;

[0095] Wherein, when the processor is configured to execute the executable instructions, it implements the engineering construction anomaly feature analysis method based on multi-angle monitoring video images described above.

[0096] The present invention also discloses a computer-readable storage medium, including:

[0097] a memory storing a computer program thereon;

[0098] a processor configured to execute the program in the memory to implement the method for analyzing abnormal features of engineering construction based on multi-angle monitoring video images as described above.

[0099] In summary, due to the adoption of the above technical solutions, the present invention can improve the monitoring efficiency, accuracy, dynamic adaptability, etc., and realize the rapid monitoring of abnormal changes in engineering construction by using satellite images, drone cruising, and on-site inspections.

[0100] Additional aspects and advantages of the present invention will be given in part in the following description, become apparent in part from the following description, or be learned through the practice of the present invention. BRIEF DESCRIPTION OF THE DRAWINGS

[0101] The above and / or additional aspects and advantages of the present invention will become apparent and be readily understood from the following description of the embodiments in conjunction with the accompanying drawings, where:

[0102] Figure 1 is a flowchart of the present invention.

[0103] Figure 2 is a schematic diagram of the present invention. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0104] The embodiments of the present invention will be described in detail below. The examples of the embodiments are shown in the accompanying drawings, where the same or similar reference numerals denote the same or similar elements or elements having the same or similar functions throughout. The embodiments described below by referring to the accompanying drawings are exemplary only for explaining the present invention and should not be construed as limiting the present invention.

[0105] The present invention discloses a method for analyzing abnormal features of engineering construction based on multi-angle monitoring video images, as Figures 1-2 shown, including the following steps:

[0106] S1. Enter the natural resources information system;

[0107] S2. Obtain the image data stored in the natural resources information system;

[0108] S3. Analyze and determine the abnormal changes in engineering construction based on the obtained image data.

[0109] In a preferred embodiment of the present invention, the method for entering the natural resources information system in step S1 includes the following steps:

[0110] S11. Enter the account number and password on the login interface of the natural resources information system; the number of digits of both the account number and the password is not less than 6 digits;

[0111] S12. After entering the account number and password on the login interface of the natural resources information system, perform security encryption processing on the entered account number and password;

[0112] S13. After performing security encryption processing on the entered account number and password, send the securely encrypted account number and password to the verification platform of the natural resources information system;

[0113] S14. After the verification platform of the natural resources information system receives the securely encrypted account number and password, decrypt the received securely encrypted account number and password;

[0114] S15. If both the decrypted account number and password exist in the verification platform of the natural resources information system, enter the natural resources information system;

[0115] Otherwise, the account number and password entered on the login interface of the natural resources information system are incorrect, and the correct account number and password need to be entered again.

[0116] In a preferred embodiment of the present invention, the method for performing security encryption processing on the entered account number and password in step S12 includes the following steps:

[0117] S121. Obtain the current moment, convert the current moment into a ten - digit standard format, where the current moment includes year, month, day, and hour. For example, if the obtained moment is 17:00 on March 26, 2025, the ten - digit standard format is 2025032617; if the obtained moment is 21:00 on April 16, 2024, the ten - digit standard format is 2024041621;

[0118] S122. Determine the size relationship between the number of digits of the entered account number and password and 10:

[0119] S1221. If the number of digits of the account number is less than 10, remove the first e1 digits from the front of the ten - digit standard format of the moment, and the number of removed digits is equal to the number of digits of the account number; proceed to the next step; e1 = 10 - e0, where e0 is the number of digits of the account number and e1 is the number of digits removed from the moment;

[0120] If the number of digits of the account number is greater than 10, add e2 digits e3 to the front of the ten - digit standard format of the moment, and the number of added digits is equal to the number of digits of the account number; proceed to the next step; e2 = e0 - 10, where e0 is the number of digits of the account number, e2 is the number of digits added to the moment, and e3 is the last digit value in the ten - digit standard format;

[0121] If the number of digits of the account number is equal to 10, proceed to the next step;

[0122] S1222, if the number of digits of the password is less than 10, then remove f1 digits from the front of the ten - digit standard - format time, and after removal, it is equal to the number of digits of the password; proceed to the next step; f1 = 10 - f0, where f0 is the number of digits of the password and f1 is the number of digits removed from the time.

[0123] If the number of digits of the password is greater than 10, then add f2 digits f3 to the front of the ten - digit standard - format time, and after addition, it is equal to the number of digits of the password; proceed to the next step; f2 = f0 - 10, where f0 is the number of digits of the password, f2 is the number of digits added to the time, and f3 is the last digit value in the ten - digit standard format.

[0124] If the number of digits of the password is equal to 10, proceed to the next step; for example, if the input account and password are Zhanghao123 and miMa3211 respectively, and the ten - digit standard - format time is 2025050515, then after processing, the times are 52025050515 and 25050515 respectively; if the input account and password are Zhanghao12332 and miMa3211222 respectively, and the ten - digit standard - format time is 2025050515, then after processing, the times are 5552025050515 and 552025050515 respectively.

[0125] At this time, the account and password remain unchanged, but the times corresponding to the account and password may change.

[0126] S123, move the first digits of the account and password backward by f4 digits respectively, where f4 = floor(f3 / 2), and f4 is the number of digits moved backward; for example, if the input account and password are Zhang1hao3 and mi2Ma2 respectively, and the ten - digit standard - format time is 2025060715, then after processing, the account and password are haZng1hao3 and i2mMa2 respectively; if the input account and password are Zhang1hao3 and mi2Ma2 respectively, and the ten - digit standard - format time is 2025060716, then the processed account and password are hanZg1hao3 and i2Mma2 respectively.

[0127] S124, convert each character in the account, password, and time into 7 - bit binary ASCII form, calculate the ASCII - form account with the ASCII - form time, and calculate the ASCII - form password with the ASCII - form time to obtain the encrypted account and encrypted password respectively.

[0128]

[0129] g′2(φ) represents the φ - th digit value in the ASCII - form encrypted account.

[0130] g2(φ) represents the φ-th digit value in the account number in ASCII form;

[0131] t2(φ) represents the φ-th digit value in the time in ASCII form;

[0132] φ = 1, 2, 3, ……, φ1; φ1 is the number of digits of the account number in ASCII form;

[0133] φ1 = 7φ2, where φ2 is the number of digits of the account number;

[0134] ◇ represents the swap value operator;

[0135]

[0136] v′2(θ) represents the θ-th digit value in the encrypted password in ASCII form;

[0137] v2(θ) represents the θ-th digit value in the password in ASCII form;

[0138] t′2(θ) represents the θ-th digit value in the time in ASCII form;

[0139] θ = 1, 2, 3, ……, θ1; θ1 is the number of digits of the password in ASCII form;

[0140] θ1 = 7θ2, where θ2 is the number of digits of the password;

[0141] Convert the encrypted account number in ASCII form and the encrypted password in ASCII form into the encrypted account number and encrypted password in character form.

[0142] For example, the input account number and password are ZHangHao12234 and MIma912 respectively, and the obtained time is 23:00 on January 22, 2025;

[0143] In step S521, the time in the ten-digit standard format is 2025012223;

[0144] In step S522, the account number is ZHangHao12234, and the time corresponding to the account number ZHangHao12234 is 3332025012223;

[0145] The password is MIma912, and the time corresponding to the password MIma912 is 5012223;

[0146] In step S523, the account number is HanZgHao12234, and the time corresponding to the account number HanZgHao12234 is 3332025012223;

[0147] The password is ImaM912, and the time corresponding to the password ImaM912 is 5012223;

[0148] In step S524, the ASCII form of the account number HanZgHao12234 is 1001000110000111011101011010110011110010001100001110111101100010110010011001001100110110100, and the ASCII form of the time corresponding to the time 3332025012223 corresponding to the account number HanZgHao12234 is 0110011011001101100110110010011000001100100110101011000001100010110010011001001100100110011; The encrypted account number in ASCII form is 1111011101001010111011101000101011111110101010100101111100000000000000000000000000010000111, and the corresponding encrypted account number is {R]hWzT_NUT NUT NUTSOH BEL, as shown in Table 1.

[0149] Table 1 ASCII corresponding to each character in the account number, time, and encrypted account number

[0150] Character ASCII Character ASCII Character ASCII H 1001000 3 0110011 { 1111011 a 1100001 3 0110011 R 1010010 n 1101110 3 0110011 ] 1011101 Z 1011010 2 0110010 h 1101000 g 1100111 0 0110000 W 1010111 H 1001000 2 0110010 z 1111010 a 1100001 5 0110101 T 1010100 o 1101111 0 0110000 _ 1011111 1 0110001 1 0110001 NUT 0000000 2 0110010 2 0110010 NUT 0000000 2 0110010 2 0110010 NUT 0000000 3 0110011 2 0110010 SOH 0000001 4 0110100 3 0110011 BEL 0000111

[0151] The ASCII form of the account number of the password ImaM912 is 1001001110110111000011001101011100101100010110010, and the ASCII form of the time corresponding to the time 5012223 corresponding to the password ImaM912 is 0110101011000001100010110010011001001100100110011; The encrypted password in ASCII form is 1111100101110110100001111111000101100000110000001, and the corresponding encrypted password is |]P DEL VT ETX SOH, as shown in Table 2.

[0152] Table 2 ASCII corresponding to each character in the password, time, and encrypted password

[0153] Character ASCII Character ASCII Character ASCII I 1001001 5 0110101 | 1111100 m 1101101 0 0110000 ] 1011101 a 1100001 1 0110001 P 1010000 M 1001101 2 0110010 DEL 1111111 9 0111001 2 0110010 VT 0001011 1 0110001 2 0110010 ETX 0000011 2 0110010 3 0110011 SOH 0000001

[0154] In a preferred embodiment of the present invention, the method for decrypting the received securely encrypted account and password in step S54 includes the following steps:

[0155] S541. Obtain the current moment and convert the current moment into a ten - digit standard format;

[0156] S542. Determine the magnitude relationship between the number of digits of the received encrypted account and the encrypted password and 10:

[0157] S5221. If the number of digits of the encrypted account is less than 10, then remove e′1 digits from the front of the ten - digit standard format moment, and after removal, it is equal to the number of digits of the encrypted account; proceed to the next step; e′1 = 10 - e′0, where e′0 is the number of digits of the encrypted account, and e1′ is the number of digits removed from the moment;

[0158] If the number of digits of the encrypted account is greater than 10, then add e′2 digits e′3 to the front of the ten - digit standard format moment, and after addition, it is equal to the number of digits of the encrypted account; proceed to the next step; e′2 = e′0 - 10, where e′0 is the number of digits of the encrypted account, e′2 is the number of digits added to the moment, and e′3 is the last digit value in the ten - digit standard format;

[0159] If the number of digits of the encrypted account is equal to 10, proceed to the next step;

[0160] S5222. If the number of digits of the encrypted password is less than 10, then remove f1′ digits from the front of the ten - digit standard format moment, and after removal, it is equal to the number of digits of the encrypted password; proceed to the next step; f1′ = 10 - f0′, where f0′ is the number of digits of the encrypted password, and f1′ is the number of digits removed from the moment;

[0161] If the number of digits of the encrypted password is greater than 10, then add f2′ digits f3′ to the front of the ten - digit standard format moment, and after addition, it is equal to the number of digits of the encrypted password; proceed to the next step; f2′ = f0′ - 10, where f0′ is the number of digits of the encrypted password, f2′ is the number of digits added to the moment, and f3′ is the last digit value in the ten - digit standard format;

[0162] If the number of digits of the encrypted password is equal to 10, proceed to the next step;

[0163] At this time, the encrypted account and the encrypted password remain unchanged, but the moments corresponding to the encrypted account and the encrypted password may change;

[0164] S543. Move the first f4′ digits of the encrypted account and the encrypted password one position backward as a whole, where f4′ = floor(f3′ / 2), and f4′ is the number of digits as a whole;

[0165] S544. Convert each character in the encrypted account, encrypted password, and time into 7-bit binary ASCII form. Calculate the ASCII-form encrypted account with the ASCII-form time, and calculate the ASCII-form encrypted password with the ASCII-form time to obtain the decrypted account and decrypted password respectively.

[0166]

[0167] g″′2(φ′) represents the φ′-th digit value in the decrypted account in ASCII form;

[0168] g″2(φ′) represents the φ′-th digit value in the encrypted account in ASCII form;

[0169] t″2(φ′) represents the φ′-th digit value in the time in ASCII form;

[0170] φ′ = 1, 2, 3, ……, φ′1; φ′1 is the number of digits of the encrypted account in ASCII form;

[0171] φ′1 = 7φ′2, where φ′2 is the number of digits of the encrypted account;

[0172] ◇ represents the swap value operator;

[0173]

[0174] v″′2(θ′) represents the θ′-th digit value in the decrypted password in ASCII form;

[0175] v″2(θ′) represents the θ′-th digit value in the encrypted password in ASCII form;

[0176] t″2(θ′) represents the θ′-th digit value in the time in ASCII form;

[0177] θ′ = 1, 2, 3, ……, θ′1; θ′1 is the number of digits of the encrypted password in ASCII form;

[0178] θ′1 = 7θ′2, where θ′2 is the number of digits of the encrypted password;

[0179] Convert the decrypted account and decrypted password in ASCII form into character-form decrypted account and decrypted password.

[0180] In a preferred embodiment of the present invention, in step S2, the image data stored in the natural resource information system is sourced from obtaining image data by one or any combination of satellite imagery, drone cruising, and on-site inspections. The method of storing it includes the following steps:

[0181] S21. Obtain the pixel values of each channel in the engineering construction monitoring image:

[0182]

[0183] C j represents the pixel value of the j-th channel, where j = r, g, b;

[0184] (x p , y q ) j represents the pixel value of the j-th channel at the position (x p , y q );

[0185] p = 1, 2, 3, ……, P, where P is the number of horizontal pixel points;

[0186] q = 1, 2, 3, ……, Q, where Q is the number of vertical pixel points;

[0187] When j = r, the pixel value of the red channel is:

[0188]

[0189] When j = g, the pixel value of the green channel is:

[0190]

[0191] When j = b, the pixel value of the blue channel is:

[0192]

[0193] C r represents the pixel value of the r-th channel (red channel), (x p , y q ) r represents the pixel value of the r-th channel (red channel) at the position (x p , y q );

[0194] C g represents the pixel value of the g-th channel (green channel), (x p , y q ) g represents the pixel value of the g-th channel (green channel) at the position (x p , y q );

[0195] C b represents the pixel value of the b-th channel (blue channel), (x p , y q ) b represents the pixel value of the b-th channel (blue channel) at the position (xp , y q ), the pixel value of

[0196] S22, determine whether the pixel values of each channel in the engineering construction monitoring image are binary pixel values that only contain 0 and 1:

[0197] If the pixel values of each channel in the engineering construction monitoring image are binary pixel values that only contain 0 and 1, then perform the next step;

[0198] If the pixel values of each channel in the engineering construction monitoring image are not binary pixel values that only contain 0 and 1, then change the pixel values of each channel in the engineering construction monitoring image to binary pixel values that only contain 0 and 1, and perform the next step;

[0199] S23, determine whether the number of bits of each pixel value in each channel of the engineering construction monitoring image is consistent with the preset bit number threshold:

[0200] If the number of bits of each pixel value in each channel of the engineering construction monitoring image is consistent with the preset bit number threshold, then perform the next step;

[0201] If the number of bits of the pixel values in each channel of the engineering construction monitoring image is inconsistent with the preset bit number threshold, then fill 0 at the front of the pixel values in each channel of the engineering construction monitoring image to make the number of bits of the pixel values in each channel of the engineering construction monitoring image consistent with the preset bit number threshold, and perform the next step;

[0202]

[0203] C′ j represents the pixel value whose number of bits of the pixel value of the j-th channel is consistent with the preset bit number threshold, j = r, g, b;

[0204] (x p , y q )′ j represents the pixel value whose number of bits of the pixel value at the position (x p , y q ) of the j-th channel is consistent with the preset bit number threshold,

[0205] p = 1, 2, 3,..., P, where P is the number of horizontal pixel points;

[0206] q = 1, 2, 3,..., Q, where Q is the number of vertical pixel points;

[0207] When j = r, the pixel value whose number of bits of the pixel value of the red channel is consistent with the preset bit number threshold is:

[0208]

[0209] When j = g, the pixel values of the green channel with the number of bits consistent with the preset bit threshold are:

[0210]

[0211] When j = b, the pixel values of the blue channel with the number of bits consistent with the preset bit threshold are:

[0212]

[0213] C′ r represents the pixel value of the r-th channel (red channel) with the number of bits consistent with the preset bit threshold, (x p , y q )′ r represents the pixel value of the r-th channel (red channel) at the position (x p , y q ) with the number of bits consistent with the preset bit threshold;

[0214] C′ g represents the pixel value of the g-th channel (green channel) with the number of bits consistent with the preset bit threshold, (x p , y q )′ g represents the pixel value of the g-th channel (green channel) at the position (x p , y q ) with the number of bits consistent with the preset bit threshold;

[0215] C′ b represents the pixel value of the b-th channel (blue channel) with the number of bits consistent with the preset bit threshold, (x p , y q )′ b represents the pixel value of the b-th channel (blue channel) at the position (x p , y q ) with the number of bits consistent with the preset bit threshold;

[0216] S24, connect the pixel values of each channel:

[0217]

[0218] is the channel-connected pixel value of the j-th channel;

[0219] (x p , y q )′ j represents the pixel value of the j-th channel at the position (x p , y q ) with the number of bits consistent with the preset bit threshold;

[0220] p = 1, 2, 3, ……, P, where P is the number of horizontal pixel points;

[0221] q = 1, 2, 3, ……, Q, where Q is the number of vertical pixel points;

[0222] When j = r, the channel - connected pixel value of the red channel is:

[0223]

[0224] is the channel - connected pixel value of the red channel;

[0225] (x p , y q )′ r represents the pixel value whose number of bits of the pixel value of the red channel at position (x p , y q ) is the same as the preset number - of - bits threshold;

[0226] p = 1, 2, 3, ……, P, where P is the number of horizontal pixel points;

[0227] q = 1, 2, 3, ……, Q, where Q is the number of vertical pixel points;

[0228] When j = g, the channel - connected pixel value of the green channel is:

[0229]

[0230] is the channel - connected pixel value of the green channel;

[0231] (x p , y q )′ g represents the pixel value whose number of bits of the pixel value of the green channel at position (x p , y q ) is the same as the preset number - of - bits threshold;

[0232] p = 1, 2, 3, ……, P, where P is the number of horizontal pixel points;

[0233] q = 1, 2, 3, ……, Q, where Q is the number of vertical pixel points;

[0234] When j = b, the channel - connected pixel value of the blue channel is:

[0235]

[0236] is the channel - connected pixel value of the blue channel;

[0237] (xp , y q )′ b represents the pixel value whose number of bits of the pixel value at position (x p , y q ) is consistent with the preset bit number threshold;

[0238] p = 1, 2, 3,..., P, where P is the number of horizontal pixel points;

[0239] q = 1, 2, 3,..., Q, where Q is the number of vertical pixel points;

[0240] S25. Fully connect the channel - connected pixel values of each channel to obtain the fully - connected pixel value:

[0241]

[0242] where F cen is the fully - connected pixel value;

[0243] is the channel - connected pixel value of the red channel;

[0244] is the channel - connected pixel value of the green channel;

[0245] is the channel - connected pixel value of the blue channel;

[0246] S26. Generate the fully - connected merged pixel value according to the fully - connected pixel value:

[0247] F mes = coalesceway(F cen ),

[0248] where F mes represents the merged pixel value;

[0249] F cen is the fully - connected pixel value;

[0250] coalesceway() represents the merged pixel value calculation function, preferably using the SHA - 3(384) function, and the merged pixel value is a 384 - bit value containing only 0 and 1;

[0251]

[0252] where F neu represents the fully - connected merged pixel value;

[0253] F mes represents the merged pixel value;

[0254] K represents the number of merged pixel values;

[0255]

[0256] Among them, K represents the number of merged pixel values;

[0257] floor() represents the rounding function;

[0258] P is the number of horizontal pixel points;

[0259] Q is the number of vertical pixel points;

[0260] represents the number of bits of the merged pixel value;

[0261] S27. Compare the pixel values of each channel in step S23 with the fully connected merged pixel values to obtain the swapped pixel values of each channel:

[0262]

[0263] Among them, represents the swapped value of;

[0264] represents the value of the p ,y q )′ j in the th bit;

[0265] F neu (τ) represents the value of the τ-th bit in the fully connected merged pixel value F neu ;

[0266] ◇ represents the swap value operator;

[0267] p, q, P, τ, η satisfy:

[0268]

[0269] Among them, τ is τ in F neu (τ);

[0270] η is the preset bit threshold, equal to the depth of the engineering construction monitoring image;

[0271] P is the number of horizontal pixel points;

[0272] p is in which p is the subscript value of the abscissa;

[0273] q is in which q is the subscript value of the ordinate;

[0274] is in

[0275] When j = r, the exchange value is:

[0276]

[0277] represents the value of the p -th q bit in the pixel value (x r , y );

[0278] (x p , y q )′ r represents the pixel value in the red channel at position (x p , y q ) whose number of bits is consistent with the preset bit threshold;

[0279] When j = g, the exchange value is:

[0280]

[0281] represents the value of the p -th q bit in the pixel value (x g , y );

[0282] (x p , y q )′ g represents the pixel value in the green channel at position (x p , y q ) whose number of bits is consistent with the preset bit threshold;

[0283] When j = b, the exchange value is:

[0284]

[0285] represents the value of the p -th q bit in the pixel value (x b , y );

[0286] (x p , y q )′ b represents the pixel value in the blue channel at position (x p , y qpixel values whose number of bits is the same as the preset bit threshold;

[0287] S28, the swapped pixel values of each channel are:

[0288]

[0289] C″ j indicating the swapped pixel value of the j-th channel, where j = r, g, b;

[0290] (x p ,y q )″ j indicating the swapped pixel value of the j-th channel at the position (x p ,y q );

[0291] p = 1, 2, 3,..., P, where P is the number of horizontal pixel points;

[0292] q = 1, 2, 3,..., Q, where Q is the number of vertical pixel points;

[0293] When j = r, the swapped pixel value of the red channel is:

[0294]

[0295] C″ r indicating the swapped pixel value of the red channel;

[0296] (x p ,y q )″ r indicating the swapped pixel value of the red channel at the position (x p ,y q );

[0297] When j = g, the swapped pixel value of the green channel is:

[0298]

[0299] C″ g indicating the swapped pixel value of the green channel;

[0300] (x p ,y q )″ g indicating the swapped pixel value of the green channel at the position (x p ,y q );

[0301] When j = b, the swapped pixel value of the blue channel is:

[0302]

[0303] C″b Represents the exchanged pixel value of the blue channel;

[0304] (x p , y q )″ b Represents the exchanged pixel value of the blue channel at the position (x p , y q );

[0305] S29. Generate an image based on the pixel values of each channel and store it in the natural resource information system. After generating the image, the combined pixel value F mes As the search string of the stored image, the stored image can be searched using the search string, and the search string is stored in the date storage set named after the date corresponding to the image shooting date corresponding to the search string. That is, if the shooting date is March 9, 2025, the search string is stored in the file named 20250309; if the shooting date is April 9, 2024, the search string is stored in the file named 20240409; if the shooting date is February 2, 2025, the search string is stored in the file named 20250202.

[0306] In a preferred embodiment of the present invention, the method for obtaining the image data stored in the natural resource information system in step S2 includes the following steps:

[0307] S2-1. Obtain the selected date, and the selected date includes the selection of the reference date and the comparison date; the reference date includes year, month, and day, the comparison date includes year, month, and day, and the reference date and the comparison date are not the same day, that is, different dates;

[0308] S2-2. Obtain the search string stored on that date according to the selected date;

[0309] S2-3. Search for the stored image corresponding to the search string according to the search string, and use the search string and the stored image corresponding to the search string to obtain the engineering construction monitoring images of different dates;

[0310] S2-4. Determine the abnormal changes in the engineering construction according to the engineering construction monitoring images of different dates.

[0311] In a preferred embodiment of the present invention, the method for obtaining the engineering construction monitoring images of different dates by using the search string and the image corresponding to the search string in step S2-3 includes the following steps:

[0312] S2-31. Obtain the pixel values of each channel in the stored image:

[0313]

[0314] D j represents the pixel value of the j-th channel in the stored image, where j = r, g, b;

[0315] (x′ p ,y′ q ) j represents the pixel value of the j-th channel in the stored image at the position (x′ p ,y′ q );

[0316] p = 1, 2, 3, ……, P, where P is the number of horizontal pixel points;

[0317] q = 1, 2, 3, ……, Q, where Q is the number of vertical pixel points;

[0318] When j = r, the pixel value of the red channel in the stored image is:

[0319]

[0320] When j = g, the pixel value of the green channel in the stored image is:

[0321]

[0322] When j = b, the pixel value of the blue channel in the stored image is:

[0323]

[0324] D r represents the pixel value of the r-th channel (red channel) in the stored image, (x′ p ,y′ q ) r represents the pixel value of the r-th channel (red channel) in the stored image at the position (x p ,y q );

[0325] D g represents the pixel value of the g-th channel (green channel) in the stored image, (x′ p ,y′ q ) g represents the pixel value of the g-th channel (green channel) in the stored image at the position (x p ,y q );

[0326] D b represents the pixel value of the b-th channel (blue channel) in the stored image, (x′ p ,y′ q ) b represents the pixel value of the b-th channel (blue channel) in the stored image at the position (xp , y q ); the pixel value of

[0327] S2-32. Determine whether the pixel values of each channel in the stored image are binary pixel values containing only 0 and 1:

[0328] If the pixel values of each channel in the stored image are binary pixel values containing only 0 and 1, then proceed to the next step;

[0329] If the pixel values of each channel in the stored image are not binary pixel values containing only 0 and 1, then change the pixel values of each channel in the stored image to binary pixel values containing only 0 and 1, and proceed to the next step;

[0330] S2-33. Determine whether the number of bits of each pixel value in each channel of the stored image is consistent with the preset bit number threshold:

[0331] If the number of bits of each pixel value in each channel of the stored image is consistent with the preset bit number threshold, then proceed to the next step;

[0332] If the number of bits of the pixel values in each channel of the stored image is inconsistent with the preset bit number threshold, then fill 0 at the front of the pixel values in each channel of the stored image to make the number of bits of the pixel values in each channel of the stored image consistent with the preset bit number threshold, and proceed to the next step;

[0333]

[0334] D′ j represents the pixel value whose number of bits of the pixel value of the j-th channel in the stored image is consistent with the preset bit number threshold, where j = r, g, b;

[0335] (x′ p , y′ q )′ j represents the pixel value whose number of bits of the pixel value at the position (x′ p , y′ q ) of the j-th channel in the stored image is consistent with the preset bit number threshold,

[0336] p = 1, 2, 3,..., P, where P is the number of horizontal pixel points;

[0337] q = 1, 2, 3,..., Q, where Q is the number of vertical pixel points;

[0338] When j = r, the pixel value whose number of bits of the pixel value of the red channel in the stored image is consistent with the preset bit number threshold is:

[0339]

[0340] When j = g, the pixel value whose number of bits of the pixel value of the green channel in the stored image is consistent with the preset bit number threshold is:

[0341]

[0342] When j = b, the pixel values in the blue channel of the stored image whose number of bits is the same as the preset bit threshold are:

[0343]

[0344] D′ r represents the pixel values in the r-channel (red channel) of the stored image whose number of bits is the same as the preset bit threshold, (x′ p ,y′ q )′ r represents the pixel values in the r-channel (red channel) of the stored image at the position (x′ p ,y′ q ) whose number of bits is the same as the preset bit threshold;

[0345] D′ g represents the pixel values in the g-channel (green channel) of the stored image whose number of bits is the same as the preset bit threshold, (x′ p ,y′ q )′ g represents the pixel values in the g-channel (green channel) of the stored image at the position (x′ p ,y′ q ) whose number of bits is the same as the preset bit threshold;

[0346] D′ b represents the pixel values in the b-channel (blue channel) of the stored image whose number of bits is the same as the preset bit threshold, (x′ p ,y′ q )′ b represents the pixel values in the b-channel (blue channel) of the stored image at the position (x′ p ,y′ q ) whose number of bits is the same as the preset bit threshold;

[0347] S2-34, concatenate all K′ search strings to obtain the fully connected search string:

[0348]

[0349] where, F str represents the fully connected search string;

[0350] F sea represents the search string;

[0351] K′ represents the number of search strings;

[0352]

[0353] Among them, K′ represents the number of search strings;

[0354] floor() represents the rounding function;

[0355] P is the number of horizontal pixel points;

[0356] Q is the number of vertical pixel points;

[0357] represents the number of digits of the search string;

[0358] S2-35. Compare the pixel values of each channel in step S2-33 with the fully connected search string to obtain the swapped pixel values of each channel:

[0359]

[0360] Among them, represents the swapped value of;

[0361] represents the value of the p ,y′ q )′ j in the th digit;

[0362] F str (τ) represents the value of the τth digit in the fully connected search string F str ; F str represents the fully connected search string;

[0363] ◇ represents the swap value operator;

[0364] p, q, P, τ, η satisfy:

[0365]

[0366] Among them, τ is the τ in F str (τ);

[0367] η is the preset digit threshold, equal to the depth of the stored image;

[0368] P is the number of horizontal pixel points;

[0369] p is the subscript value of the abscissa in;

[0370] q is the subscript value of the ordinate in;

[0371] For in

[0372] When j = r, the exchange value is:

[0373]

[0374] denotes the value of the p -th q bit in the pixel value (x′ r in ;

[0375] (x′ p , y′ q )′ r denotes the pixel value in the red channel at position (x′ p , y′ q ) whose number of bits is the same as the preset bit threshold;

[0376] When j = g, the exchange value is:

[0377]

[0378] denotes the value of the p -th q bit in the pixel value (x′ g in ;

[0379] (x′ p , y′ q )′ g denotes the pixel value in the green channel at position (x′ p , y′ q ) whose number of bits is the same as the preset bit threshold;

[0380] When j = b, the exchange value is:

[0381]

[0382] denotes the value of the p -th q bit in the pixel value (x′ b in ;

[0383] (x′ p , y′ q )′ bThe pixel value indicating the blue channel at position (x′ p , y′ q ) is consistent with the preset bit threshold;

[0384] S2-36, store the swapped pixel values of each channel in the image as:

[0385]

[0386] D″ j represents the swapped pixel value of the j-th channel in the stored image, where j = r, g, b;

[0387] (x′ p , y′ q )″ j represents the swapped pixel value of the j-th channel in the stored image at position (x′ p , y′ q );

[0388] p = 1, 2, 3, ……, P, where P is the number of horizontal pixel points;

[0389] q = 1, 2, 3, ……, Q, where Q is the number of vertical pixel points;

[0390] When j = r, the swapped pixel value of the red channel in the stored image is:

[0391]

[0392] D″ r represents the swapped pixel value of the red channel in the stored image;

[0393] (x′ p , y′ q )″ r represents the swapped pixel value of the red channel in the stored image at position (x′ p , y′ q );

[0394] When j = g, the swapped pixel value of the green channel in the stored image is:

[0395]

[0396] D″ g represents the swapped pixel value of the green channel in the stored image;

[0397] (x′ p , y′ q )″ g represents the swapped pixel value of the green channel in the stored image at position (x′ p , y′ q );

[0398] When j = b, the swapped pixel value of the blue channel in the stored image is:

[0399]

[0400] D″ b represents the swapped pixel value of the blue channel in the stored image;

[0401] (x′ p ,y′ q )″ b represents the swapped pixel value of the blue channel at the position (x′ p ,y′ q ) in the stored image;

[0402] S2-37. Generate an image based on the swapped pixel values of each channel, and thus an image is obtained, which is an image generated from a search string and the image retrieved by the search string. All other search strings are processed in the same way to obtain all the images for two periods.

[0403] It further includes

[0404] S2-371. In step S2-37, it further includes channel connecting the connected swapped pixel values of each channel in the stored image:

[0405]

[0406] is the channel-connected swapped pixel value of the j-th channel;

[0407] (x′ p ,y′ q )″ j represents the swapped pixel value of the j-th channel at the position (x′ p ,y′ q ) in the stored image;

[0408] p = 1, 2, 3,..., P, where P is the number of horizontal pixel points;

[0409] q = 1, 2, 3,..., Q, where Q is the number of vertical pixel points;

[0410] When j = r, the channel-connected swapped pixel value of the red channel is:

[0411]

[0412] is the channel-connected swapped pixel value of the red channel;

[0413] (x′ p ,y′q )′ r represents the swapped pixel value at position (x′ p , y′ q ) in the red channel of the stored image;

[0414] p = 1, 2, 3, ……, P, where P is the number of horizontal pixel points;

[0415] q = 1, 2, 3, ……, Q, where Q is the number of vertical pixel points;

[0416] When j = g, the channel - connected pixel value of the green channel is:

[0417]

[0418] is the channel - connected swapped pixel value of the green channel;

[0419] (x′ p , y′ q )′ g represents the swapped pixel value at position (x′ p , y′ q ) in the green channel of the stored image;

[0420] p = 1, 2, 3, ……, P, where P is the number of horizontal pixel points;

[0421] q = 1, 2, 3, ……, Q, where Q is the number of vertical pixel points;

[0422] When j = b, the channel - connected pixel value of the blue channel is:

[0423]

[0424] is the channel - connected swapped pixel value of the blue channel;

[0425] (x′ p , y′ q )′ b represents the swapped pixel value at position (x′ p , y′ q ) in the blue channel of the stored image;

[0426] p = 1, 2, 3, ……, P, where P is the number of horizontal pixel points;

[0427] q = 1, 2, 3, ……, Q, where Q is the number of vertical pixel points;

[0428] S2 - 372, fully connect the channel - connected swapped pixel values of each channel in the stored image to obtain the fully - connected swapped pixel value:

[0429]

[0430] Among them, F' cen is the fully connected exchange pixel value;

[0431] is the channel-connected exchange pixel value of the red channel;

[0432] is the channel-connected exchange pixel value of the green channel;

[0433] is the channel-connected exchange pixel value of the blue channel;

[0434] S2-373. Generate the fully connected merged exchange pixel value according to the fully connected exchange pixel value:

[0435] F' mes = coalesceway(F' cen ),

[0436] Among them, F' mes represents the merged exchange pixel value;

[0437] F' cen is the fully connected exchange pixel value;

[0438] coalesceway() represents the merged pixel value calculation function, preferably using the SHA-3(384) function, and the merged exchange pixel value is a 384-bit value containing only 0 and 1;

[0439] S2-374. Determine whether the merged exchange pixel value F' mes is consistent with the search string:

[0440] If the merged exchange pixel value F' mes is consistent with the search string, generate an image according to the exchange pixel values of each channel, and thus obtain a correct image, that is, a correct image generated by a search string and the image searched by the search string. All other search strings are also executed according to the corresponding steps to obtain all correct images in two periods;

[0441] If the merged exchange pixel value F' mes is inconsistent with the search string, delete the generated image and calculate the images corresponding to other search strings and the search string.

[0442] In a preferred embodiment of the present invention, the method for determining the abnormal changes in engineering construction according to the engineering construction monitoring images on different dates in step S2-4 includes the following steps:

[0443] S2-41, obtain the pixel values of each channel in the engineering construction monitoring image:

[0444]

[0445] E j represents the pixel value of the j-th channel, where j = r, g, b;

[0446] (x″ p ,y″ q ) j represents the pixel value of the j-th channel at the position (x″ p ,y″ q );

[0447] p = 1, 2, 3, ……, P, where P is the number of horizontal pixel points;

[0448] q = 1, 2, 3, ……, Q, where Q is the number of vertical pixel points;

[0449] When j = r, the pixel value of the red channel is:

[0450]

[0451] When j = g, the pixel value of the red channel is:

[0452]

[0453] When j = b, the pixel value of the red channel is:

[0454]

[0455] E r represents the pixel value of the r-th channel (red channel), (x″ p ,y″ q ) r represents the pixel value of the r-th channel (red channel) at the position (x″ p ,y″ q );

[0456] E g represents the pixel value of the g-th channel (green channel), (x″ p ,y″ q ) g represents the pixel value of the g-th channel (green channel) at the position (x″ p ,y″ q );

[0457] E b represents the pixel value of the b-th channel (blue channel), (x″ p ,y″ q ) bIndicates the pixel value at position (x″ p , y″ q ) in the b-th channel (blue channel);

[0458] S2-42. Generate the pixel values of the merged channel according to the pixel values of each channel in the obtained engineering construction monitoring image:

[0459]

[0460] That is:

[0461]

[0462] Where E rgb Indicates the pixel value of the merged channel;

[0463] δ r Indicates the pixel value coefficient of the red channel;

[0464] δ g Indicates the pixel value coefficient of the green channel;

[0465] δ b Indicates the pixel value coefficient of the blue channel;

[0466] δ r (x″ p , y″ q ) r + δ g (x″ p , y″ q ) g + δ r (x″ p , y″ q ) b Indicates the pixel value of the merged channel at position (x″ p , y″ q );

[0467] p = 1, 2, 3, ……, P, where P is the number of horizontal pixel points;

[0468] q = 1, 2, 3, ……, Q, where Q is the number of vertical pixel points;

[0469] S2-43. Determine the pixel values for the pixel values on the same date:

[0470] S2-431. Arrange all the pixel values at the same pixel coordinates in ascending order, which are Pix1, Pix2, Pix3, ……, Pix H ;

[0471] S2-432. Take As the pixel values participating in the calculation, the first determined pixel value is obtained:

[0472]

[0473] Among them, Pix1 represents the pixel value ranked first after arranging the pixel values of the same pixel coordinate in ascending order;

[0474] Pix2 represents the pixel value ranked second after arranging the pixel values of the same pixel coordinate in ascending order;

[0475] Pix3 represents the pixel value ranked third after arranging the pixel values of the same pixel coordinate in ascending order;

[0476] Pix H represents the pixel value ranked at the H-th after arranging the pixel values of the same pixel coordinate in ascending order;

[0477] represents the pixel value ranked at the h1-th after arranging the pixel values of the same pixel coordinate in ascending order;

[0478] represents the pixel value ranked at the h2-th after arranging the pixel values of the same pixel coordinate in ascending order;

[0479] Pix con1 represents the first determined pixel value;

[0480] H represents the number of images taken on the same date;

[0481]

[0482] S2-433, take all the pixel values Pix1, Pix2, Pix3,..., Pix of the same pixel coordinate H As the pixel values participating in the calculation, the second determined pixel value is obtained:

[0483]

[0484] Among them, Pix con2 represents the second determined pixel value;

[0485] H represents the number of images taken on the same date;

[0486] Pix j represents the pixel value ranked at the j-th after arranging the pixel values of the same pixel coordinate in ascending order;

[0487] S2-434, determine the magnitude relationship between the first determined pixel value and the second determined pixel value:

[0488] If |Pix con1 -Pix con2 | ≤ Pix sd , then the second determined pixel value is the determined pixel value;

[0489] If |Pix con1 -Pix con2 | > Pix sd , then the first determined pixel value is the determined pixel value; Pix sd is the preset pixel value exceeding the threshold;

[0490] The pixel value of a pixel coordinate can be determined through steps S2-431 to S2-434, and the pixel values of other pixel coordinates are determined in the same way; correspondingly, the image pixel values of the reference date and the comparison date can be determined;

[0491] S2-44, perform pixel value anomaly determination on the pixel values of different dates:

[0492] If then there is no abnormal change in the engineering construction at the pixel point (x″′ p , y″′ q ); is the preset first abnormal change threshold,

[0493] If then there is an abnormal change in the engineering construction at the pixel point (x″′ p , y″′ q );

[0494] Among them, then there is a level III abnormal change in the engineering construction at the pixel point (x″′ p , y″′ q ); is the preset second abnormal change threshold;

[0495] then there is a level II abnormal change in the engineering construction at the pixel point (x″′ p , y″′ q ); is the preset third abnormal change threshold;

[0496] then there is a level I abnormal change in the engineering construction at the pixel point (x″′ p , y″′ q );

[0497] p = 1, 2, 3,..., P, where P is the number of horizontal pixel points;

[0498] q = 1, 2, 3, ……, Q, where Q is the number of vertical pixel points; Pix dat (x″′ p , y″′ q ) is the pixel value of the reference date image at the pixel point (x″′ p , y″′ q ); Pix com (x″′ p , y″′ q ) is the pixel value of the comparison date image at the pixel point (x″′ p , y″′ q );

[0499] The abnormal change of level III is lower than that of level II, and the abnormal change of level I is higher than that of level II.

[0500] It also includes that if any one or any combination of level I abnormal change, level II abnormal change, and level III abnormal change occurs, the geographical location of the engineering construction where the image is located is obtained, and its occupied area in mu is:

[0501] S = nξ / ψ,

[0502] n is the number of abnormal change pixel points;

[0503] ξ is the proportionality coefficient of the actual area to the image area;

[0504] ψ is the area conversion unit, and its value is 666.7;

[0505] S is the abnormal change area in mu;

[0506]

[0507] R is the image resolution; pixels / cm;

[0508] P is the number of horizontal pixel points;

[0509] Q is the number of vertical pixel points;

[0510] s represents the actual area corresponding to the image.

[0511] The present invention also discloses a computer system, including:

[0512] A processor;

[0513] A memory for storing executable instructions of the processor;

[0514] Wherein, the processor is configured to implement the engineering construction abnormal feature analysis method based on multi-angle monitoring video images when executing the executable instructions.

[0515] The present invention also discloses a computer-readable storage medium, comprising:

[0516] a memory storing a computer program thereon;

[0517] a processor configured to execute the program in the memory to implement the method for analyzing abnormal features of engineering construction based on multi-angle monitored video images as described above.

[0518] Although embodiments of the present invention have been shown and described, those of ordinary skill in the art can understand that various changes, modifications, substitutions, and variations can be made to these embodiments without departing from the principles and spirit of the present invention. The scope of the present invention is defined by the claims and their equivalents.

Claims

1. A method for analyzing abnormal features in engineering construction based on multi-angle monitored video images, characterized in that, It includes the following steps: S1. Enter the natural resources information system; S2. Obtain the image data stored in the natural resources information system; S3. Analyze the obtained image data to determine the abnormal changes in the project construction.

2. The engineering construction anomaly feature analysis method based on multi-angle monitored video images according to claim 1, wherein The method for determining the abnormal changes in the project construction according to the project construction monitoring images stored in the natural resources information system in step S3 includes the following steps: S31. Obtain the selected date, where the selected date includes the selection of the reference date and the comparison date; the reference date includes year, month, and day, the comparison date includes year, month, and day, and the reference date and the comparison date are not the same day, i.e., different dates; S32. Obtain the search string stored on that date according to the selected date; S33. Search for the stored image corresponding to the search string according to the search string, and use the search string and the stored image corresponding to the search string to obtain the project construction monitoring images of different dates; S34. Determine the abnormal changes in the project construction according to the project construction monitoring images of different dates.

3. The engineering construction anomaly feature analysis method based on multi-angle monitoring video images according to claim 1, wherein The method for determining the abnormal changes in the project construction according to the project construction monitoring images of different dates in step S34 includes the following steps: S341. Obtain the pixel values of each channel in the project construction monitoring image: E j represents the pixel value of the j-th channel, where j = r, g, b; (x′ p ′,y′ q ′) j represents the pixel value of the j-th channel at the position (x′ p ′,y′ q ′), p = 1, 2, 3, ……, P, where P is the number of horizontal pixel points; q = 1, 2, 3, ……, Q, where Q is the number of vertical pixel points; S342. Generate the pixel values of the merged channel according to the pixel values of each channel in the obtained project construction monitoring image: Among them, E rgb represents the pixel value of the merged channel; δ r represents the pixel value coefficient of the red channel; δ g represents the pixel value coefficient of the green channel; δ b represents the pixel value coefficient of the blue channel; δ r (x′ p ′,y′ q ′) r +δ g (x′ p ′,y′ q ′) g +δ r (x′ p ′,y′ q ′) b represents the pixel value of the merged channel at the position (x′ p ′,y′ q ′); p = 1, 2, 3, ……, P, where P is the number of horizontal pixel points; q = 1, 2, 3, ……, Q, where Q is the number of vertical pixel points; S343. Determine the pixel values of the same date: S3431, arrange all pixel values at the same pixel coordinate in ascending order, which are Pix1, Pix2, Pix3, ……, Pix H ; S3432, take as the pixel value participating in the calculation to obtain the first determined pixel value: Among them, Pix1 represents the pixel value ranked first after arranging the pixel values of the same pixel coordinates in ascending order; Pix2 represents the pixel value ranked second after arranging the pixel values of the same pixel coordinates in ascending order; Pix3 represents the pixel value ranked third after arranging the pixel values of the same pixel coordinates in ascending order; Pix H represents the pixel value at the H-th position after arranging the pixel values at the same pixel coordinates in ascending order; Indicates the pixel value that is in the h1-th position after arranging the pixel values of the same pixel coordinates in ascending order; It represents the pixel value located at the h2nd position after arranging the pixel values of the same pixel coordinates in ascending order successively; Pix con1 represents a first determined pixel value; H represents the number of images taken on the same date; S3433, take all pixel values Pix1, Pix2, Pix3, ……, Pix at the same pixel coordinate H as the pixel values participating in the calculation to obtain the second determined pixel value: Among them, Pix con2 represents the second determined pixel value; H represents the number of images taken on the same date; Pix j represents the pixel value at the j-th position after arranging the pixel values at the same pixel coordinates in ascending order; S3434. Determine the magnitude relationship between the first determined pixel value and the second determined pixel value: If |Pix con1 -Pix con2 | ≤ Pix sd , then the second determined pixel value is the determined pixel value; If |Pix con1 -Pix con2 | > Pix sd , then the first determined pixel value is the determined pixel value; Pix sd is the preset pixel value exceeding the threshold; S344. Determine the pixel value anomalies of different dates: If then there is no abnormal change in the engineering construction at the pixel point (x′ p ″, y′ q ″); is a preset first abnormal change threshold, If then there is an abnormal change in the engineering construction at the pixel point (x′ p ″, y′ q ″); Among them, then a level III abnormal change occurs in the engineering construction at the pixel point (x′ p ″, y′ q ″); is the preset second abnormal change threshold; Then, a Class II abnormal change occurred in the engineering construction at the pixel point (x′ p ″, y′ q ″); is the preset third abnormal change threshold; Then a Class-I abnormal change occurred in the engineering construction at the pixel point (x′ p ″, y′ q ″); p = 1, 2, 3, ……, P, where P is the number of horizontal pixel points; q = 1, 2, 3, ……, Q, where Q is the number of vertical pixel points; Pix dat (x′ p ″, y′ q ″) is the pixel value of the reference date image at the pixel point (x′″ p , y′″ p ); Pix com (x′″ p , y′″ q ) is the pixel value of the comparison date image at the pixel point (x′″ p , y′″ q ); The abnormal changes of level III are lower than those of level II, and the abnormal changes of level I are higher than those of level II.

4. The engineering construction anomaly feature analysis method based on multi-angle monitoring video images according to claim 1, characterized in that The method for entering the natural resources information system in step S1 includes the following steps: S11. Enter the account number and password on the login interface of the natural resources information system; the number of digits of both the account number and the password is not less than 6 digits; S12. After entering the account number and password on the login interface of the natural resources information system, perform security encryption processing on the entered account number and password; S13. After performing security encryption processing on the entered account number and password, send the securely encrypted account number and password to the verification platform of the natural resources information system; S14. After the natural resource information system verification platform receives the securely encrypted account and password, it decrypts the received securely encrypted account and password. S15. If both the decrypted account and password exist in the natural resource information system verification platform, then enter the natural resource information system. Otherwise, the account and password entered on the natural resource information system login interface are incorrect, and the correct account and password need to be entered again.

5. The engineering construction anomaly feature analysis method based on multi-angle monitoring video images according to claim 1, wherein The method for securely encrypting the input account and password in step S12 includes the following steps: S121. Obtain the current moment, and convert the current moment into a ten - digit standard format, where the current moment includes year, month, day, and hour. S122. Determine the size relationship between the number of digits of the input account and password and 10: S1221. If the number of digits of the account is less than 10, then remove e1 digits from the front of the ten - digit standard format moment, and after removal, it is equal to the number of digits of the account; proceed to the next step; e1 = 10 - e0, where e0 is the number of digits of the account and e1 is the number of digits removed from the moment. If the number of digits of the account is greater than 10, then add e2 digits e3 to the front of the ten - digit standard format moment, and after addition, it is equal to the number of digits of the account; proceed to the next step; e2 = e0 - 10, where e0 is the number of digits of the account, e2 is the number of digits added to the moment, and e3 is the last digit value in the ten - digit standard format. If the number of digits of the account is equal to 10, proceed to the next step. S1222. If the number of digits of the password is less than 10, then remove f1 digits from the front of the ten - digit standard format moment, and after removal, it is equal to the number of digits of the password; proceed to the next step; f1 = 10 - f0, where f0 is the number of digits of the password and f1 is the number of digits removed from the moment. If the number of digits of the password is greater than 10, then add f2 digits f3 to the front of the ten - digit standard format moment, and after addition, it is equal to the number of digits of the password; proceed to the next step; f2 = f0 - 10, where f0 is the number of digits of the password, f2 is the number of digits added to the moment, and f3 is the last digit value in the ten - digit standard format. If the number of digits of the password is equal to 10, proceed to the next step. At this time, the account and password remain unchanged, but the moments corresponding to the account and password may change. S123. Move the first digits of the account and password backward by f4 digits respectively, where f4 = floor(f3 / 2), and f4 is the number of digits moved backward. S124. Convert each character in the account, password, and moment into 7 - bit binary ASCII form, calculate the ASCII - form account with the ASCII - form moment, and calculate the ASCII - form password with the ASCII - form moment, respectively obtaining the encrypted account and encrypted password.

6. A computer system, characterized in that, Including: A processor; A memory for storing processor - executable instructions; Wherein, when the processor is configured to execute the executable instructions, it implements the engineering construction anomaly feature analysis method based on multi - angle monitoring video images as described in any one of claims 1 - 5.

7. A computer-readable storage medium, characterized in that, Including: A memory on which a computer program is stored; A processor for executing the program in the memory to implement the engineering construction anomaly feature analysis method based on multi - angle monitoring video images as described in any one of claims 1 - 5.

Images