Method for calculating the effort of an airborne geophysical survey based on the vector product

Abstract

The application discloses a kind of methods for the statistical work of aerial geophysical survey line based on vector product, it is related to aerial geophysical survey line technical field, comprising: reading the aerial geophysical survey line database of survey area, obtains the survey line number, first and last point plane coordinates and survey line length of each survey line;According to the preset survey line number rule, if the current survey line is original survey line, then the current survey line original length is included in survey area workload;If the current survey line is a connecting survey line, then the first and last point coordinates of the current survey line are extracted, and the first and last point coordinates of the corresponding original survey line are obtained based on the survey line number rule index, the slope of the original survey line is calculated;The projection coordinates of the first and last point coordinates of the current connecting line on the straight line where the original survey line is located are calculated, and the position of the first and last point projection coordinates on the original survey line is judged by constructing vector product;The first and last node of the current connecting line is reconstructed, and the length of the deduplicated connecting line is calculated;The length of all original survey lines and deduplicated connecting lines is accumulated.The application solves the problem of low efficiency in the prior art.

Description

Technical Field

[0001] This invention relates to the field of airborne object detection line technology, and in particular to a method for calculating the workload of airborne object detection lines based on vector product. Background Technology

[0002] During airborne geophysical exploration, due to various reasons such as weather and flight routes, it is sometimes impossible to complete the entire survey line in one go. In order to ensure the accuracy of subsequent data processing, the connection head must cross the cutting line during the connection process. Therefore, when there is a connection in the survey line, the length of the survey line will be repeated.

[0003] When counting the workload of aerial geophysical surveys, it is necessary to deduct redundant work. The usual practice is to draw flight path charts, then calculate the distance, and then subtract the calculated work from the total workload. This purely manual method is prone to missing lines, and the accuracy of the calculation will also bring a large error to the statistical workload, affecting the statistical accuracy. Most importantly, it is inefficient.

[0004] The key to improving the efficiency of field data processing is to quickly and effectively obtain the survey line workload.

[0005] The information disclosed in this background section is intended only to enhance the understanding of the overall background of the invention and should not be construed as an admission or in any way implying that the information constitutes prior art known to those skilled in the art. Summary of the Invention

[0006] The purpose of this invention is to provide a method for calculating the workload of airborne object detection lines based on vector product, which solves the problems of low statistical accuracy and low efficiency in existing technologies.

[0007] To achieve the above objectives, this invention provides a method for calculating the workload of airborne object detection lines based on vector product, comprising the following steps:

[0008] A method for calculating the workload of airborne object detection lines based on vector product, characterized by comprising:

[0009] S1: Read the airborne geophysical survey line database of the survey area, obtain the survey line number, the plane coordinates of the first and last points and the survey line length of each survey line. The survey line includes the original survey line and the connecting survey line generated by the connecting flight. The connecting head crosses the cutting line, resulting in the connecting line and the corresponding original survey line having inter-segment overlap.

[0010] S2: Based on the preset line numbering rules, determine whether the current line is an original line or a connecting line;

[0011] If the current survey line is the original survey line, then the original length of the current survey line will be included in the workload of the survey area;

[0012] If the current survey line is a connecting survey line, extract the coordinates of the first and last points of the current survey line, and obtain the coordinates of the first and last points of the corresponding original survey line based on the survey line number rule index;

[0013] S3: Calculate the slope of the original survey line based on the coordinates of the first and last points of the original survey line;

[0014] S4: Based on the slope of the original survey line, calculate the projected coordinates of the first and last points of the current connection on the straight line of the original survey line to obtain the projected coordinates of the first and last points.

[0015] S5: Construct a vector product based on the projected coordinates of the first and last points and the coordinates of the first and last points of the original survey line; determine the position of the projected coordinates of the first and last points on the original survey line based on the numerical sign of the vector product.

[0016] S6: Based on the position of the first and last point projection coordinates on the original survey line, reconstruct the first and last nodes of the current wiring and calculate the wiring length after deduplication;

[0017] S7: Sum the lengths of all original survey lines and the deduplicated connection lines to obtain the total workload of the survey area.

[0018] Furthermore, in S2, the preset survey line numbering rule is as follows:

[0019] When the same survey line number corresponds to multiple survey line records, the first record is the original survey line, and the rest are connection survey lines according to the record order.

[0020] The test line number is composed of "original test line unique identifier + connection sequence number".

[0021] Furthermore, in S2, if the current survey line is a connecting survey line, then the coordinates of the first point in the current connecting survey line record are read. , ) and the coordinates of the final point ( , ), and according to the survey line number rules, query the survey line database to obtain the coordinates of the first point of the original survey line corresponding to the current connection survey line ( , ) and the coordinates of the final point ( , ).

[0022] Furthermore, S3: Based on the coordinates of the first and last points of the original survey line, calculate the slope of the original survey line, including;

[0023] Based on the coordinates of the first and last points of the original survey line ( , ), ( , ), calculate the slope k of the original survey line:

[0024]

[0025] in, The slope of the original survey line, ( , ) represents the coordinates of the first point of the original survey line, ( , ) represents the coordinates of the end point of the original survey line.

[0026] Furthermore, S4: Based on the slope of the original survey line, calculate the projected coordinates of the first and last points of the current connection on the straight line of the original survey line, and obtain the projected coordinates of the first and last points, including:

[0027] Based on the slope of the original survey line and the coordinates of the first and last points of the current connection, the projected coordinates of the first point are obtained. , ) and the final point projection coordinates ( , The formula is:

[0028] ,

[0029] ,

[0030] in, The slope of the original survey line, ( , ) represents the coordinates of the first point of the original survey line, ( , ) represents the coordinates of the starting point of the current connection. , ) represents the coordinates of the end point of the current connection. , ) are the projection coordinates of the first point, ( , ) represents the projected coordinates of the final point.

[0031] Further, S5: Based on the projected coordinates of the first and last points and the coordinates of the first and last points of the original survey line, construct a vector product, and determine the position of the projected coordinates of the first and last points on the original survey line based on the numerical sign of the vector product, including:

[0032] Based on the projection coordinates of the first and last points ( , ), ( , ) and the coordinates of the beginning and end points of the original survey line ( , ), ( , Construct the first vector product. Second vector product ;

[0033] like >0, determine the projection coordinates of the first point ( , Located at the end of the original survey line ( , One side;

[0034] like <0, determine the projection coordinates of the first point ( , Located at the starting point of the original survey line ( , One side;

[0035] like =0, determine the projection coordinates of the first point ( , ) and the first point of the original survey line ( , )coincide;

[0036] like >0, determine the projection coordinates of the final point ( , Located at the end of the original survey line ( , One side;

[0037] like <0, determine the projection coordinates of the final point ( , Located at the starting point of the original survey line ( , One side;

[0038] like =0, determine the projection coordinates of the final point ( , ) and the end point of the original survey line ( , )coincide.

[0039] Furthermore, S6: Based on the projected coordinates of the first and last points on the original survey line, reconstruct the first and last nodes of the current wiring, and calculate the deduplicated wiring length, including:

[0040] like <0 and >0, and ( , )and( , If the connection between the two lines covers the original survey line segment, it is determined that the current connection completely covers the original survey line, and the deduplication length of the current connection is set to 0.

[0041] like <0 and >0, and ( , )and( , The line connecting the two points covers the original survey line segment, then ( , ) as the first node of the reconstructed wiring, retaining the original last node of the wiring. , )constant;

[0042] like >0 and >0, and ( , If it is located within the original survey line segment, then it is determined by ( , ) as the reconstructed end node of the wiring, retaining the original start node of the wiring ( , )constant;

[0043] like <0 and If <0, it is determined that the current wiring is located in the extension direction of the first point of the original survey line and has no overlap, and the first and last nodes of the original wiring are retained unchanged.

[0044] like >0 and If the value is >0, it is determined that the current wiring is located in the extension direction of the end point of the original measurement line and there is no overlap. The start and end nodes of the original wiring are retained unchanged.

[0045] Furthermore, it also includes procedures for handling anomalies where the original survey line is vertical:

[0046] When detected - When the slope k does not exist due to the value of 0, the projected coordinates are set. = , = ; = , = Overlap detection is performed directly based on the vertical coordinate.

[0047] Furthermore, in S7, the accumulation operation is performed after all survey line records in the survey area database have been read, or it is immediately accumulated into the total workload after each survey line has been read and processed.

[0048] On the other hand, the present invention also provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to implement the method for calculating the workload of airborne object detection lines based on vector product as described above.

[0049] On the other hand, the present invention also provides a non-transitory computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the method for calculating the workload of airborne object detection lines based on vector product as described above.

[0050] On the other hand, the present invention also provides a computer program product, including a computer program that, when executed by a processor, implements the method for calculating the workload of airborne object detection lines based on vector product as described above.

[0051] Compared with existing technologies, the present invention provides a method for calculating the workload of aerial object detection lines based on vector product. Traditional methods rely on manually drawing flight path maps and calculating distances, which are easily affected by human factors, leading to missed lines or measurement errors. The present invention introduces a vector product (dot product) sign determination rule. By calculating the projected coordinates of the beginning and end points of the connection line on the original survey line and constructing the vector product, the spatial position relationship of the projection point relative to the original survey line segment (first-side extension, intra-segment extension, or last-side extension) is accurately determined by the positive or negative sign of the value. This determination logic based on analytical geometry eliminates the visual error of manual map reading and improves the accuracy of workload statistics to the level of the original data accuracy.

[0052] Existing survey line index reconstruction technologies mainly focus on the repetitive judgment of the entire survey line (index level), which is difficult to handle the problem of local segment overlap caused by multiple wiring within a single survey line. This invention addresses the industry characteristic of airborne geophysical exploration that "the wiring head must cross the cutting line" and focuses on segment-level deduplication within a single survey line. By identifying the coverage area of ​​the projection point, it dynamically adjusts the first and last nodes of the wiring and accurately deducts the length of the physical overlap caused by "crossing the cutting line", avoiding the crude drawback of the traditional method of "one-size-fits-all" deduction of the entire wiring length.

[0053] This invention achieves automatic association and matching of original survey lines and connections through preset survey line number rules (such as "original identifier + connection sequence number" or "group sorting"), eliminating the need for manual screening. Combined with the fast calculation logic of vector product, it can compress the manual statistical work that originally required several hours to be completed in seconds. In particular, for complex survey areas with large areas, multiple sorties, and frequent connection, the efficiency improvement is particularly significant, meeting the urgent needs of modern airborne geophysical survey field rapid evaluation and data processing.

[0054] This invention not only provides a standard algorithm based on slope projection, but also specially designs an anomaly handling mechanism for when the original survey line is vertical. In this scenario, the algorithm automatically switches to a direct judgment mode based on the ordinate, avoiding division by zero errors in the calculation process. This fault-tolerant design ensures the stability and reliability of the algorithm when processing survey lines at any azimuth angle, and guarantees the smooth execution of the entire automated process. Attached Figure Description

[0055] Figure 1 This is a flowchart of a method for calculating the workload of airborne object detection lines based on vector product according to an embodiment of the present invention;

[0056] Figure 2 This is a schematic diagram of an electronic device based on a method for calculating the workload of airborne object detection lines according to an embodiment of the present invention. Detailed Implementation

[0057] The specific embodiments of the present invention will now be described in detail with reference to the accompanying drawings, but it should be understood that the scope of protection of the present invention is not limited to the specific embodiments.

[0058] Unless otherwise expressly stated, throughout the specification and claims, the term "comprising" or its variations such as "including" or "comprises" shall be understood to include the stated elements or components without excluding other elements or other components.

[0059] like Figure 1 As shown, a method for calculating the workload of airborne object detection lines based on vector product according to a preferred embodiment of the present invention is characterized by comprising:

[0060] S1: Read the airborne geophysical survey line database of the survey area, obtain the survey line number, the plane coordinates of the first and last points and the survey line length of each survey line. The survey line includes the original survey line and the connecting survey line generated by the connecting flight. The connecting head crosses the cutting line, resulting in the connecting line and the corresponding original survey line having inter-segment overlap.

[0061] S2: Based on the preset line numbering rules, determine whether the current line is an original line or a connecting line;

[0062] If the current survey line is the original survey line, then the original length of the current survey line will be included in the workload of the survey area;

[0063] If the current survey line is a connecting survey line, extract the coordinates of the first and last points of the current survey line, and obtain the coordinates of the first and last points of the corresponding original survey line based on the survey line number rule index;

[0064] S3: Calculate the slope of the original survey line based on the coordinates of the first and last points of the original survey line;

[0065] S4: Based on the slope of the original survey line, calculate the projected coordinates of the first and last points of the current connection on the straight line of the original survey line to obtain the projected coordinates of the first and last points.

[0066] S5: Construct a vector product based on the projected coordinates of the first and last points and the coordinates of the first and last points of the original survey line; determine the position of the projected coordinates of the first and last points on the original survey line based on the numerical sign of the vector product.

[0067] S6: Based on the position of the first and last point projection coordinates on the original survey line, reconstruct the first and last nodes of the current wiring and calculate the wiring length after deduplication;

[0068] S7: Sum the lengths of all original survey lines and the deduplicated connection lines to obtain the total workload of the survey area.

[0069] In this embodiment of the invention, relying on the preset survey line numbering rules (the first line in the same survey line number group is the original survey line, or the naming rule of "original identifier + connection sequence number"), the system can automatically complete the association matching of the original survey line and the connection line. Taking the measured data of a certain aeromagnetic survey area in a northern basin as an example, the survey area contains a total of 1278 survey line records, of which 213 are connection survey lines. The system completes the type determination and corresponding relationship binding of all survey lines in only 0.2 seconds, without the need for manual screening, and completely solves the problem of "missed connection line judgment and mismatched original survey line" that is easy to occur in the traditional mode of manual verification.

[0070] The technique of using slope projection and vector product to determine the location of the projection point can be accurately located simply by relying on the sign of the vector product: for example, for a certain original survey line L 203 Corresponding wiring L 203_J02 The first vector product V1 = −12.3 < 0 (located in the direction of the extension of the first point of the original survey line) is calculated for the first projection point, and the second vector product V2 = 45.6 > 0 (located in the direction of the extension of the last point of the original survey line) is calculated for the last projection point. The system automatically determines that the connection completely covers the original survey line segment and directly sets the length after deduplication to 0. For extreme cases where the original survey line is in the vertical direction, the abnormal handling logic automatically switches to the vertical coordinate judgment. No manual adjustment of the algorithm is required throughout the process, and the judgment accuracy reaches 100%, which is far higher than the error level of 3% to 8% of traditional manual measurement.

[0071] After completing the reconstruction of the first and last nodes and the deduplication length calculation for a single connection, the system supports two modes: real-time accumulation and batch accumulation. Taking the above-mentioned 1278 test lines as an example, the total time for the entire statistical process is only 1.7 seconds, and the deviation between the statistical results and the results of manual calculation is only 0.03%, which is far below the industry's allowable error threshold of 1%. This is more than 400 times more efficient than traditional manual statistics (average time of 2 hours per test area), and can meet the needs of rapid acceptance and real-time calculation of workload at field operation sites.

[0072] This invention strictly adheres to the industry standard that "the connector must cross the cutting line". It only deducts the effective length of the overlap between segments and will not delete the effective flight data by mistake. For example, if a connector only crosses the cutting line by 5 meters, the system can accurately identify and deduct the 5-meter overlap segment, and retain the remaining effective flight length. This makes up for the defect of existing index-level deduplication technology that directly deletes the entire connector. The statistical results are completely consistent with the actual situation of field operations.

[0073] The method directly interfaces with the GeoProbe format GDB database and can complete the entire calculation process by reading the native fields of the survey line index. Unlike traditional algorithms, it does not require complex data format conversion and can be directly embedded into existing airborne geophysical data processing pipelines, resulting in extremely low engineering implementation costs.

[0074] like Figure 1 As shown, in S2, the preset survey line numbering rule is as follows:

[0075] When the same survey line number corresponds to multiple survey line records, the first record is the original survey line, and the rest are connection survey lines according to the record order.

[0076] The test line number is composed of "original test line unique identifier + connection sequence number".

[0077] In this embodiment of the invention, traditional line identification relies on manual viewing of flight track charts or reviewing flight logs, which is not only inefficient but also prone to omissions and misjudgments due to personnel fatigue or lack of experience. This embodiment, by solidifying the rules of "the first line of the same survey line number is the original, and the rest are connections" and "original identifier + connection sequence number," endows the system with automatic parsing capabilities. Regardless of whether the survey area contains thousands or tens of thousands of survey line records, the system can instantly complete the classification according to the rules, achieving a 100% identification accuracy. For example, when the system reads survey line number L... 100 When three records are found, the first one is automatically identified as the original survey line, and the second one (L) is identified as the original survey line. 100_J01 ) and Article 3 (L 100_J02 This is a wiring method that requires no human intervention, completely solving the subjectivity problem of manual verification;

[0078] This rule is flexibly designed and compatible with the two mainstream paradigms of current airborne geophysical data storage: For compact data tables: when multiple records of the same survey line number are stored consecutively, the "record order" can be used to determine the data without modifying the database structure; For standardized data tables: when original survey lines and connecting survey lines are stored together, the "sponge naming method" (such as L) is used. 100_J01 The original survey line identifier can be deduced by parsing the string; this invention can be adapted to both old historical archived data and newly collected standardized data, greatly reducing the technical promotion threshold.

[0079] This rule is not only a classification basis, but also an "index key" for data association. By parsing the "unique identifier of the original survey line" in the connection survey line number, the system can quickly locate and retrieve the coordinates of the first and last points of the corresponding original survey line in the database. This strong association based on naming rules is faster and more unambiguous than algorithms based on spatial distance fuzzy matching, ensuring that the input data for projection calculation in step S4 and deduplication calculation in step S6 is absolutely accurate.

[0080] like Figure 1 As shown, in S2, if the current survey line is a connecting survey line, then the coordinates of the first point in the current connecting survey line record are read. , ) and the coordinates of the final point ( , ), and according to the survey line number rules, query the survey line database to obtain the coordinates of the first point of the original survey line corresponding to the current connection survey line ( , ) and the coordinates of the final point ( , ).

[0081] In this embodiment of the invention, the source logic of the coordinate data is clearly defined: the coordinates of the beginning and end points of the current connection ( , ), ( , The reading method reads directly from its own records, rather than through the original survey line or other complex transformations. This "self-contained" reading method ensures that the geometry of the line itself (such as bending and offset) is completely preserved. It will not lose the actual flight trajectory information of the line due to correlation query, and provides the most authentic original data for subsequent accurate calculation of overlap length.

[0082] Existing technologies often use fuzzy matching based on spatial distance to associate tracks and survey lines, which can easily lead to ambiguity in areas with densely packed survey lines. This invention utilizes survey line number rules as a "precise index" to directly locate the corresponding original survey line record in the database and extract its first and last points. , ), ( , This deterministic query based on naming rules is like finding a person by their ID number. It ensures that even in complex survey areas with crisscrossing survey lines and extremely small spacing, the original survey line can be accurately located in milliseconds, fundamentally eliminating the "misattribution" phenomenon that may occur in spatial matching algorithms.

[0083] This step logically completes the "pre-pruning" of the data. The system only needs to load the coordinates of the four key points (a total of eight values) of the current connection and its corresponding original survey line to carry out subsequent calculations. Compared with algorithms such as global search or buffer analysis that require loading a large amount of surrounding survey line data, this method controls the data scale at the constant level (O(1)), making the processing time of a single connection negligible and ensuring that the overall statistical process still maintains a second-level response when facing tens of thousands of connections.

[0084] The original survey line coordinates obtained in this step provide a unified geometric benchmark for subsequent projection calculation (S4) and vector product determination (S5). Since the coordinates of the original survey line and the connection line both come from the same survey line database (such as GeoProbeGDB) and are in the same coordinate system, the cumulative errors caused by different data sources or coordinate transformations are avoided, ensuring that the calculation results of the vector product truly reflect the spatial position relationship, thereby guaranteeing the accuracy of the final deduplication length.

[0085] like Figure 1 As shown, S3: Calculate the slope of the original survey line based on the coordinates of the first and last points of the original survey line, including;

[0086] Based on the coordinates of the first and last points of the original survey line ( , ), ( , ), calculate the slope k of the original survey line:

[0087]

[0088] in, The slope of the original survey line, ( , ) represents the coordinates of the first point of the original survey line, , ) represents the coordinates of the end point of the original survey line.

[0089] In this embodiment of the invention, the slope k quantitatively describes the inclination and direction of the original survey line (such as due south, due north, or southeast, northwest), and is the mathematical cornerstone for defining the "straight line where the original survey line is located". Through this calculation, the system abstracts the discrete coordinates of the first and last points into a continuous linear equation. This transformation makes it possible to calculate the projected coordinates of the first and last points of the connection in the subsequent S4 step, which is equivalent to establishing a precise "reference axis" for overlap detection on a two-dimensional plane.

[0090] This step transforms complex spatial geometric positional relationships into simple algebraic operations. Compared to manual measurement or complex spatial topology operations that rely on graphical interfaces, slope-based calculations have extremely high execution efficiency. Calculating a slope only requires one subtraction and one division, with a computational complexity of O(1). This enables the algorithm to easily handle the processing needs of tens of thousands of wiring lines in large-scale survey areas, ensuring a second-level response throughout the entire process.

[0091] Regardless of whether the original survey line is approximately horizontal, vertical, or tilted at an arbitrary angle, the slope calculation formula provided in this step is applicable. This means that the algorithm is not limited by the direction of the survey line layout. Whether it is a regular grid-like survey line or an irregular survey line along the terrain, it can accurately capture its geometric features, greatly expanding the applicability of the method.

[0092] like Figure 1 As shown, S4: Based on the slope of the original survey line, calculate the projected coordinates of the first and last points of the current connection on the straight line of the original survey line, and obtain the projected coordinates of the first and last points, including:

[0093] Based on the slope of the original survey line and the coordinates of the first and last points of the current connection, the projected coordinates of the first point are obtained. , ) and the final point projection coordinates ( , The formula is:

[0094] ,

[0095] ,

[0096] in, The slope of the original survey line, ( , ) represents the coordinates of the first point of the original survey line, , ) represents the coordinates of the starting point of the current connection. , ) represents the coordinates of the end point of the current connection. , ) are the projection coordinates of the first point, ( , ) represents the projected coordinates of the final point.

[0097] In this embodiment of the invention, the principle of vertical projection is used to visualize the start and end points of the wiring ( , ), ( , Mapping the original survey line onto an infinitely long straight line, we obtain the projected coordinates. , ), ( , The essence of this operation is to reduce the problem of line segment overlap detection on a two-dimensional plane to the problem of interval coverage detection on a one-dimensional straight line. Through projection, the system can ignore the small offset of the line in the direction perpendicular to the measurement line (such as instrument drift or wind direction), accurately capture its true overlap position in the flight direction of the measurement line, and provide an accurate linear measurement benchmark for subsequent determination of "crossing the cutting line".

[0098] The projection formula used in this invention is not a simple approximation calculation, but a precise derivation of the "coordinates of the foot of the perpendicular from a point to a line" in analytical geometry. This formula directly uses the slope k of the original survey line for calculation, avoiding the accuracy loss caused by converting it to an angle system. Especially when dealing with long-distance survey lines, this formula can effectively suppress error accumulation and ensure that the accuracy of the coordinate calculation of the projected point is consistent with the accuracy of the original GPS measurement data, thus ensuring the accuracy of the workload statistics from the source.

[0099] Projected coordinates ( , ), ( , The projected coordinates are the core variables for constructing the vector products V1 and V2 in the subsequent S5 step. Without the precise location of the projection points, it is impossible to determine whether the connection is located at the beginning, end, or middle of the original survey line by using the sign of the vector product. Furthermore, once overlap is determined, these projected coordinates themselves become the "anchor points" for reconstructing the beginning and end nodes of the connection (for example, when it is determined that the overlap at the beginning needs to be deducted, the coordinates are directly selected). , (This serves as a new starting point), giving the deduplication operation a clear geometric basis.

[0100] like Figure 1 As shown, S5: Based on the projected coordinates of the first and last points and the coordinates of the first and last points of the original survey line, construct a vector product, and determine the position of the projected coordinates of the first and last points on the original survey line based on the numerical sign of the vector product, including:

[0101] Based on the projection coordinates of the first and last points ( , ), ( , ) and the coordinates of the beginning and end points of the original survey line ( , ), ( , Construct the first vector product. Second vector product ;

[0102] like >0, determine the projection coordinates of the first point ( , Located at the end of the original survey line ( , One side;

[0103] like <0, determine the projection coordinates of the first point ( , Located at the starting point of the original survey line ( , One side;

[0104] like =0, determine the projection coordinates of the first point ( , ) and the first point of the original survey line ( , )coincide;

[0105] like >0, determine the projection coordinates of the final point ( , Located at the end of the original survey line ( , One side;

[0106] like <0, determine the projection coordinates of the final point ( , Located at the starting point of the original survey line ( , One side;

[0107] like =0, determine the projection coordinates of the final point ( , ) and the end point of the original survey line ( , )coincide.

[0108] In this embodiment of the invention, the sign of the vector product is used as the basis for judgment, and the geometric position judgment is transformed into a simple algebraic sign comparison. Compared with the commonly used "calculation of vertical distance", "angle discrimination" or "complex conditional branch statements" in traditional technology, this method only requires one multiplication and one sign judgment, and the computational complexity is extremely low (O(1)). This "branchless" algorithm design greatly improves the pipeline execution efficiency of the CPU, enabling the system to maintain a microsecond-level single-point judgment speed when processing massive survey line data, which is the key to realizing full-process second-level statistics.

[0109] By combining the signs of V1 and V2, the system can accurately distinguish three critical states of the projection point: negative half-axis (V<0): the projection point is located outside the first point of the original survey line (the direction of the first point's extension); positive half-axis (V>0): the projection point is located outside the last point of the original survey line (the direction of the last point's extension); zero point (V=0): the projection point coincides with the first / last point. This refined judgment capability allows the algorithm to strictly follow the airborne geophysical survey specification that "the connection point must cross the cutting line"—only when the projection point crosses from V<0 to V>0 (i.e., crosses the first or last point) can it be judged as a valid coverage or overlap, thus avoiding misjudgment caused by the ambiguity of the threshold setting in the existing technology.

[0110] In conventional geometric calculations, determining whether a point lies on a line segment often uses distance formulas (such as calculating whether the sum of the distances from the point to the endpoint equals the length of the line segment). This method is highly susceptible to floating-point rounding errors, causing points that should coincide to be judged as non-coincident. This invention uses the sign of the vector product to determine the sign, focusing only on the positive or negative attribute of the value. It is not sensitive to small numerical fluctuations. Unless the projected point happens to fall on the endpoint (within the allowable error range), the sign determination result has extremely strong stability. It effectively shields the computational noise caused by rounding multiple decimal places in GPS coordinates, ensuring the robustness of the algorithm in complex data environments.

[0111] like Figure 1 As shown, S6: Based on the projected coordinates of the first and last points on the original survey line, reconstruct the first and last nodes of the current wiring, and calculate the deduplicated wiring length, including:

[0112] like <0 and >0, and ( , )and( , If the connection between the two lines covers the original survey line segment, it is determined that the current connection completely covers the original survey line, and the deduplication length of the current connection is set to 0.

[0113] like <0 and >0, and ( , )and( , The line connecting the two points covers the original survey line segment, then ( , ) as the first node of the reconstructed wiring, retaining the original last node of the wiring. , )constant;

[0114] like >0 and >0, and ( , If it is located within the original survey line segment, then it is determined by ( , ) as the reconstructed end node of the wiring, retaining the original start node of the wiring ( , )constant;

[0115] like <0 and If <0, it is determined that the current wiring is located in the extension direction of the first point of the original survey line and has no overlap, and the first and last nodes of the original wiring are retained unchanged.

[0116] like >0 and If the value is >0, it is determined that the current wiring is located in the extension direction of the end point of the original measurement line and there is no overlap. The start and end nodes of the original wiring are retained unchanged.

[0117] In this embodiment of the invention, differentiated node reconstruction strategies are formulated for different overlap patterns, realizing a closed loop from "qualitative judgment" to "quantitative deduction": Complete coverage is zeroed out: when V1<0 and V2>0 and the connection covers the entire segment, the length is decisively set to 0, avoiding the problem in traditional methods where the entire connection is mistakenly included in the effective workload due to "the connection head crossing the cutting line"; Local overlap is trimmed: when it is determined to be partial coverage (such as V1<0, V2>0 and the connection partially covers), the first node after reconstruction is accurately positioned to the projection point (x6, y6), only the overlapping segment is deducted, and the effective segment is completely retained; this refined processing capability makes the statistical results no longer rely on the "empirical value" estimated by humans, but accurate to the "true value" of the coordinate data, reducing the workload statistical error from a few percent as common in the industry to a few ten-thousandths;

[0118] The five decision branches designed in this step constitute a complete truth table, covering all possible spatial relationships between the wiring and the original survey line: cross-end coverage (extending from the first side to the last side); embedded at the first end (projection point within the segment, first node unchanged); embedded at the last end (projection point within the segment, last node unchanged); no overlap at the first end (all outside the first point); no overlap at the last end (all outside the last point). Regardless of how the survey line is laid out or how the wiring flies, the algorithm can find the corresponding processing logic, eliminating the risk of "missed deductions" or "erroneous deductions" caused by logical loopholes in existing technologies, and ensuring the rigor of the statistical results.

[0119] By "reconstructing the first and last nodes," the system actually simulates the physical process of "cutting off overlapping segments" at the mathematical level. For example, reconstructing the first node of the connection from (x2, y2) to (x6, y6) means that the system recognizes that only flights starting from the end of the original survey line are valid. This dynamic correction not only outputs an accurate length value, but also logically restores an ideal flight trajectory without redundancy. This provides important basic support for subsequent data quality assessment, flight cost settlement, and reserve calculation.

[0120] When determining overlap, this invention preserves the effective length of the non-overlapping portion to the maximum extent. For example, if only the first 10 kilometers of a certain line overlaps with the original survey line, this invention only deducts 10 kilometers and retains the remaining 90 kilometers. In contrast, traditional index deduplication might delete all 100 kilometers. This "keep as much as possible" strategy not only conforms to the "better to be too much than too little" operational standard of airborne geophysical surveys, but also strives to secure the maximum reasonable rights for enterprises in project settlement.

[0121] like Figure 1 As shown, the method for calculating the workload of airborne object survey lines based on vector product also includes an anomaly handling step for original survey lines that are vertical:

[0122] When detected - When the slope k does not exist due to the value of 0, the projected coordinates are set. = , = ; = , = Overlap determination is performed directly based on the vertical coordinate.

[0123] In this embodiment of the invention, when the original survey line is a vertical line in the standard projection formula, the denominator is zero, resulting in the absence of slope k. Without intervention, the computer will crash due to a "division by zero error". This step uses pre-detection logic to automatically bypass the complex slope calculation and projection formula when x1−x0=0 is detected, and instead executes a simplified assignment operation. This mechanism transforms potentially fatal errors into controllable logical branches, ensuring that the system can operate stably without interruption when facing survey lines of any azimuth angle.

[0124] Given the unique characteristics of vertical survey lines, this step abandons the common projection formula and simplifies the process by directly utilizing vertical geometric features: the horizontal coordinate of the projection is set to x5=x0 (the projection point falls on the vertical line), and the vertical coordinate is directly taken as the original value of the connection point, y5=y2. This processing simplifies the originally complex floating-point multiplication and division operations into simple assignment operations, with computational overhead being almost zero. In data containing a large number of north-south survey lines, this optimization can significantly improve the overall processing speed, demonstrating the ingenuity of the algorithm design.

[0125] Figure 2 This is a schematic diagram of the structure of the electronic device provided in an embodiment of the present invention.

[0126] like Figure 2 As shown, the electronic device may include a processor 610, a communications interface 620, a memory 630, and a communication bus 640. The processor 610, communications interface 620, and memory 630 communicate with each other via the communication bus 640. The processor 610 can call logical instructions from the memory 630 to execute a method for calculating the workload of airborne object detection lines based on vector product statistics.

[0127] Furthermore, the logical instructions in the aforementioned memory 630 can be implemented as software functional units and, when sold or used as independent products, can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, in essence, or the part that contributes to the prior art, or a part of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of the present invention. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0128] On the other hand, the present invention also provides a computer program product, which includes a computer program that can be stored on a non-transitory computer-readable storage medium. When the computer program is executed by a processor, the computer is able to perform the methods provided above for calculating the workload of airborne object detection lines based on vector product statistics.

[0129] In another aspect, the present invention also provides a non-transitory computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, is implemented to perform the methods provided above for calculating the workload of airborne object detection lines based on vector product statistics.

[0130] The device embodiments described above are merely illustrative. The units described as separate components may or may not be physically separate. The components shown as units may or may not be physical units; that is, they may be located in one place or distributed across multiple network units. Some or all of the modules can be selected to achieve the purpose of this embodiment according to actual needs. Those skilled in the art can understand and implement this without any creative effort.

[0131] Through the above description of the embodiments, those skilled in the art can clearly understand that each embodiment can be implemented by means of software plus necessary general-purpose hardware platforms, and of course, it can also be implemented by hardware. Based on this understanding, the above technical solutions, in essence or the part that contributes to the prior art, can be embodied in the form of a software product. This computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, magnetic disk, optical disk, etc., including several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute the methods described in the various embodiments or some parts of the embodiments.

[0132] The foregoing description of specific exemplary embodiments of the invention is for illustrative and explanatory purposes. These descriptions are not intended to limit the invention to the precise forms disclosed, and it will be apparent that many changes and variations can be made in accordance with the foregoing teachings. The exemplary embodiments were chosen and described in order to explain the specific principles of the invention and its practical application, thereby enabling those skilled in the art to implement and utilize various different exemplary embodiments of the invention, as well as various different choices and variations. The scope of the invention is intended to be defined by the claims and their equivalents.

Claims

1. A method for calculating the workload of airborne object detection lines based on vector product, characterized in that, include: S1: Read the airborne geophysical survey line database of the survey area, obtain the survey line number, the plane coordinates of the first and last points and the survey line length of each survey line. The survey line includes the original survey line and the connecting survey line generated by the connecting flight. The connecting head crosses the cutting line, resulting in the connecting line and the corresponding original survey line having inter-segment overlap. S2: Based on the preset line numbering rules, determine whether the current line is an original line or a connecting line; If the current survey line is the original survey line, then the original length of the current survey line will be included in the workload of the survey area; If the current survey line is a connecting survey line, extract the coordinates of the first and last points of the current survey line, and obtain the coordinates of the first and last points of the corresponding original survey line based on the survey line number rule index; S3: Calculate the slope of the original survey line based on the coordinates of the first and last points of the original survey line; S4: Based on the slope of the original survey line, calculate the projected coordinates of the first and last points of the current connection on the straight line of the original survey line to obtain the projected coordinates of the first and last points. S5: Construct a vector product based on the projected coordinates of the first and last points and the coordinates of the first and last points of the original survey line; determine the position of the projected coordinates of the first and last points on the original survey line based on the numerical sign of the vector product. S6: Based on the position of the first and last point projection coordinates on the original survey line, reconstruct the first and last nodes of the current wiring and calculate the wiring length after deduplication; S7: Sum the lengths of all original survey lines and the deduplicated connection lines to obtain the total workload of the survey area.

2. The method for calculating the workload of airborne object detection lines based on vector product statistics according to claim 1, characterized in that, In S2, the preset survey line numbering rule is as follows: When the same survey line number corresponds to multiple survey line records, the first record is the original survey line, and the rest are connection survey lines according to the record order. The test line number is composed of "original test line unique identifier + connection sequence number".

3. The method for calculating the workload of airborne object detection lines based on vector product statistics according to claim 1, characterized in that, In S2, if the current survey line is a connecting survey line, then read the coordinates of the first point in the current connecting survey line record. , ) and the coordinates of the final point ( , ), and according to the survey line number rules, query the survey line database to obtain the coordinates of the first point of the original survey line corresponding to the current connection survey line ( , ) and the coordinates of the final point ( , ).

4. The method for calculating the workload of airborne object detection lines based on vector product statistics according to claim 1, characterized in that, S3: Calculate the slope of the original survey line based on the coordinates of the first and last points of the original survey line, including: Based on the coordinates of the first and last points of the original survey line ( , ), ( , ), calculate the slope k of the original survey line: in, The slope of the original survey line, ( , ) represents the coordinates of the first point of the original survey line, , ) represents the coordinates of the end point of the original survey line.

5. The method for calculating the workload of airborne object detection lines based on vector product statistics according to claim 1, characterized in that, S4: Based on the slope of the original survey line, calculate the projected coordinates of the first and last points of the current connection on the straight line of the original survey line, and obtain the projected coordinates of the first and last points, including: Based on the slope of the original survey line and the coordinates of the first and last points of the current connection, the projected coordinates of the first point are obtained. , ) and the final point projection coordinates ( , The formula is: ， ， in, The slope of the original survey line, ( , ) represents the coordinates of the first point of the original survey line, , ) represents the coordinates of the starting point of the current connection. , ) represents the coordinates of the end point of the current connection. , ) are the projection coordinates of the first point, ( , ) represents the projected coordinates of the final point.

6. The method for calculating the workload of airborne object detection lines based on vector product statistics according to claim 1, characterized in that, S5: Based on the projected coordinates of the first and last points and the coordinates of the first and last points of the original survey line, construct a vector product. Determine the position of the projected coordinates of the first and last points on the original survey line based on the sign of the vector product, including: Based on the projection coordinates of the first and last points ( , ), ( , ) and the coordinates of the beginning and end points of the original survey line ( , ), ( , Construct the first vector product. Second vector product ; like >0, determine the projection coordinates of the first point ( , Located at the end of the original survey line ( , One side; like <0, determine the projection coordinates of the first point ( , Located at the starting point of the original survey line ( , One side; like =0, determine the projection coordinates of the first point ( , ) and the first point of the original survey line ( , )coincide; like >0, determine the projection coordinates of the final point ( , Located at the end of the original survey line ( , One side; like <0, determine the projection coordinates of the final point ( , Located at the starting point of the original survey line ( , One side; like =0, determine the projection coordinates of the final point ( , ) and the end point of the original survey line ( , )coincide.

7. The method for calculating the workload of airborne object detection lines based on vector product statistics according to claim 1, characterized in that, S6: Based on the projected coordinates of the first and last points on the original survey line, reconstruct the first and last nodes of the current wiring, and calculate the deduplicated wiring length, including: like <0 and >0, and ( , )and( , If the connection between the two lines covers the original survey line segment, it is determined that the current connection completely covers the original survey line, and the deduplication length of the current connection is set to 0. like <0 and >0, and ( , )and( , The line connecting the two points covers the original survey line segment, then ( , ) as the first node of the reconstructed wiring, retaining the original last node of the wiring. , )constant; like >0 and >0, and ( , If it is located within the original survey line segment, then it is determined by ( , ) as the reconstructed end node of the wiring, retaining the original start node of the wiring ( , )constant; like <0 and If <0, it is determined that the current wiring is located in the extension direction of the first point of the original survey line and has no overlap, and the first and last nodes of the original wiring are retained unchanged. like >0 and If the value is >0, it is determined that the current wiring is located in the extension direction of the end point of the original measurement line and there is no overlap. The start and end nodes of the original wiring are retained unchanged.

8. The method for calculating the workload of airborne object detection lines based on vector product statistics according to claim 1, characterized in that, It also includes procedures for handling anomalies where the original survey line is vertical: When detected - When the slope k does not exist due to the value of 0, the projected coordinates are set. = , = ; = , = Overlap determination is performed directly based on the vertical coordinate.

9. The method for calculating the workload of airborne object detection lines based on vector product statistics according to claim 1, characterized in that, In S7, the accumulation operation is performed after all survey line records in the survey area database have been read, or it is immediately added to the total workload after each survey line has been read and processed.

10. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes the program, it implements the method for calculating the workload of airborne object detection lines based on vector product as described in any one of claims 1 to 9.

Images