A curved surface reconstruction modeling method and system based on section editing

Abstract

The application provides a curved surface reconstruction modeling method and system based on a section editing, and belongs to the technical field of geological curved surface reconstruction. The application utilizes a two-dimensional section line of an original geological curved surface and a two-dimensional section line after updating, calculates a characteristic line F path and a corresponding elevation according to linear interpolation, sets a sampling range radius, a calculation point selection range radius and a reference elevation point range radius, sets a re-interpolation calculation interval, splits a single long characteristic line, and for each characteristic line, the following steps are cyclically executed to determine a sampling point set, a calculation point set and a reference elevation point set of each characteristic line, the reference elevation point set is merged into the sampling point set, a new elevation of the calculation point set is recalculated by Kriging interpolation, after traversing all the characteristic lines, a grid curved surface is updated according to the recalculated new elevation, and a section drawing is regenerated. The interpolation result accuracy and reliability make the curved surface correction more accurate.

Description

Technical Field

[0001] The present invention relates to the technical field of geological surface reconstruction, and in particular to a surface reconstruction modeling method and system based on section plane editing. Background Art

[0002] Conventional 3D geological surface modeling is usually based on geological exploration data, such as drill hole data, geological profile data, etc. Discrete data points are constructed into a continuous 3D surface model through interpolation algorithms (such as Kriging interpolation and triangulation interpolation).

[0003] In actual use, there is a certain error between theoretical interpolation and actual measurement. After more detailed on-site measurements and tests, the surface needs to be modified. The existing technology is to create a similar two-dimensional polyline based on the original geological two-dimensional profile line according to user needs, and calculate the corresponding three-dimensional point information of the newly created two-dimensional polyline through the profile line. This three-dimensional point information is added to the surface modeling to realize the surface update, and the geological profile map is updated after the surface is updated.

[0004] During the modification process, due to issues such as data accuracy, calculation efficiency, model consistency maintenance, and user interactivity, the user's manually drawn two-dimensional polyline cannot accurately reflect the actual measured value. At the same time, due to too many sampling points, the interpolation calculation fails, resulting in unsuccessful surface updates. Summary of the Invention

[0005] In order to solve the technical problems existing in the prior art, the present invention provides a surface reconstruction modeling method and system based on section plane editing. By selecting an appropriate number of sampling points to participate in the calculation, the calculation points that need to be recalculated and the reference points involved in the calculation are determined, which not only ensures the data quantity, but also prevents the calculation failure caused by excessive sampling data volume. At the same time, the calculation points that need to be recalculated in elevation will not be missed. The accuracy and reliability of the Kriging interpolation results make the surface correction more accurate. The calculation is performed in segments and the surface is locally reconstructed in blocks, which will not cause data overflow. For the endpoints of repeated calculations, the average value or the value of the first calculation is used as the calculation value of the point, which will not cause calculation errors.

[0006] The present invention provides a surface reconstruction modeling method based on section editing, comprising the following steps:

[0007] Step S1, selecting the original two-dimensional section line L1 and the updated two-dimensional section line L2 of the geological surface, wherein the updated two-dimensional section line is obtained by directly drawing a polyline on the section diagram to replace the original geological section line;

[0008] Step S2, using the two-dimensional profile line L1 of the original geological surface and the updated two-dimensional profile line L2, calculate the path of the characteristic line F and the corresponding elevation according to linear interpolation;

[0009] Step S3, setting the sampling range radius R1 of the sampling points around the characteristic line F that need to participate in the interpolation calculation, the selection range radius R2 of the points around the characteristic line F that need to re-interpolate the elevation calculation, and the reference elevation point range radius R3 around the characteristic line F that need to participate in the interpolation calculation; wherein R1>R2, R3>R2;

[0010] Step S4, setting the re-interpolation calculation spacing S of the feature line F;

[0011] Step S5: determine whether the single-segment feature line F is too long. If the length of the single-segment feature line exceeds a preset length threshold, the single-segment feature line is split;

[0012] Step S6, sequentially obtain the feature line outer bounding box BOX of a feature line segment, calculate the feature line sampling point outer bounding box BOX1 of the sampling point set P1 of the feature line segment F, and the feature line recalculation point outer bounding box BOX2 of the calculation point set P2 of the feature line segment F that needs to recalculate the elevation;

[0013] Step S7, determining the sampling point set P1, calculation point set P2 and reference elevation point set P3 of the characteristic line segment;

[0014] Step S8, merging the reference elevation point set P3 of the feature line segment into its sampling point set P1;

[0015] Step S9, using Kriging interpolation to recalculate the new elevations of the calculation points in the calculation point set P2 of the segment of the characteristic line according to the sampling point set P1 of the segment of the characteristic line and the set interpolation calculation interval S;

[0016] Step S10, looping through steps S6 to S9 until all feature lines are traversed;

[0017] Step S11: update the mesh surface according to the recalculated new elevation and regenerate the profile.

[0018] Preferably, step S1 specifically includes the following steps:

[0019] Select a sub-entity from the geological two-dimensional profile, determine the original geological surface to be modified, the profile line L1 corresponding to the geological surface, the corresponding characteristic line F, and the elevation set E1 corresponding to each key point on the characteristic line F that meets the preset rules;

[0020] A polyline L2 is drawn on the geological two-dimensional section map to replace the section line in the original geological two-dimensional section map, to obtain an updated two-dimensional section line L2, and the updated two-dimensional section line L2 is selected.

[0021] Preferably, step S2 specifically includes the following steps: restoring the plane coordinates of the characteristic line F from the geological two-dimensional profile, and establishing a first key point set based on the key points on the characteristic line F that meet the preset rules; projecting the selected updated two-dimensional profile line L2 into the geological two-dimensional profile map, establishing a second key point set based on the key points on the updated two-dimensional profile line L2 that meet the preset rules, and using the projection relationship to calculate the plane coordinates and new elevation set E2 corresponding to the second key point set; inserting the key points in the first key point set that do not have a corresponding relationship with the key points of the second key point set into the second key point set, and using linear interpolation to calculate the corresponding elevation and add it to the corresponding position of the new elevation set E2.

[0022] Preferably, a Boolean AND operation is performed on the first key point set and the second key point set, and key points in the first key point set that do not correspond to key points in the second key point set are inserted into the second key point set. The specific method is as follows:

[0023] The key points in the first key point set are sequentially found in the second key point set according to their positions on the feature lines. If they exist, they are ignored. If they do not exist, the key points in the first key point set are inserted into the second key point set.

[0024] Preferably, step S9 specifically includes the following steps: adding the reference elevation point set P3 to the sampling point set P1, using the points in the reference elevation point set P3 to cover the elevations of the points in the sampling point set P1 for points with repeated plane coordinates, and adding the new elevation point set E2 corresponding to the characteristic line to the sampling point set P1 as interpolation sampling points.

[0025] Preferably, the preset length threshold in step S5 is greater than or equal to 4 times the sampling range radius R1. When the length of a single feature line segment exceeds the preset length threshold, a new key point is directly inserted into the feature line segment using linear interpolation calculation method as the endpoint of the split segment.

[0026] Preferably, step S7 includes the following steps:

[0027] According to the two vertices of the single-segment feature line, calculate an original feature line outer bounding box BOX;

[0028] According to the setting value of R1, the original bounding box BOX is expanded to a range with a radius of R1 to form the outer bounding box BOX1 of the feature line sampling point;

[0029] According to the setting value of R2, the original bounding box is expanded to a range with a radius of R2 to form a feature line and recalculate the outer bounding box BOX2 of the point.

[0030] Preferably, step S8 specifically includes the following steps:

[0031] Filter out the points within the bounding box BOX1 outside the feature line sampling point, and then filter out the sampling point set P1 in this range by the distance of the feature line;

[0032] Filter out the feature lines and recalculate the points within the bounding box BOX2, and then filter out the calculation point set P2 in this range by the distance of the feature lines;

[0033] The reference elevation points within the bounding box BOX1 are filtered out through the feature line sampling point, and then the reference elevation point set P3 within the range is filtered out through the distance to the feature line.

[0034] Preferably, step S11 specifically includes the following steps: updating the local surface of the grid according to the recalculated calculation point set P2, and regenerating the geological section line in the cross-section diagram at the original position.

[0035] The present invention also provides a surface reconstruction modeling system based on section plane editing, comprising a processor capable of executing a computer program, wherein the computer program can implement the above-mentioned surface reconstruction modeling method based on section plane editing.

[0036] Compared with the prior art, the present invention has the following beneficial effects:

[0037] 1. The present invention selects an appropriate number of sampling points to participate in the calculation, determines the calculation points that need to be recalculated and the important feature reference points involved in the calculation, which not only ensures the data quantity and eliminates the interference of too far points, but also prevents the calculation failure caused by excessive sampling data volume. At the same time, it also does not omit the calculation points that need to be recalculated, and improves the accuracy and reliability of the Kriging interpolation results, making the surface correction more accurate.

[0038] 2. The present invention adopts segmented calculation and blocks to locally reconstruct the surface, which greatly improves the calculation efficiency, reduces the loss of data accuracy, and does not cause data overflow. For the endpoints of repeated calculations, the average value or the value of the first calculation is used as the calculated value of the point, which will not cause calculation errors. BRIEF DESCRIPTION OF THE DRAWINGS

[0039] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the following briefly introduces the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on the structures shown in these drawings without paying any creative work.

[0040] Figure 1 A schematic diagram of a characteristic line calculation method according to an embodiment of the present invention;

[0041] Figure 2A schematic diagram of sampling points and calculation points in characteristic line calculation according to an embodiment of the present invention;

[0042] Figure 3 A schematic diagram of points in a common area that may be repeatedly calculated in characteristic line calculation according to an embodiment of the present invention;

[0043] Figure 4 Another schematic diagram of points in a common area that may be repeatedly calculated in characteristic line calculation according to an embodiment of the present invention;

[0044] Figure 5 A schematic diagram of a geological two-dimensional cross-section according to an embodiment of the present invention;

[0045] Figure 6 A schematic diagram of an original geological surface according to an embodiment of the present invention;

[0046] Figure 7 This is an effect diagram of reconstructed geological surface according to an embodiment of the present invention; R1=100R2=50;

[0047] Figure 8 This is an effect diagram of reconstructed geological surface according to an embodiment of the present invention; R1=200R2=100;

[0048] Figure 9 This is an effect diagram of reconstructed geological surface according to an embodiment of the present invention; R1=400R2=100;

[0049] Figure 10 This is an effect diagram of reconstructed geological surface according to an embodiment of the present invention; R1=400R2=200;

[0050] Figure 11 This is a flow chart of a surface reconstruction modeling method based on section plane editing according to an embodiment of the present invention. DETAILED DESCRIPTION

[0051] The specific embodiments of the present invention are described in detail below.

[0052] like Figure 11 As shown, the present invention provides a surface reconstruction modeling method based on section plane editing, comprising the following steps:

[0053] Step S1, selecting the original two-dimensional section line L1 and the updated two-dimensional section line L2 of the geological surface, wherein the updated two-dimensional section line is obtained by directly drawing a polyline on the section diagram to replace the original geological section line;

[0054] Step S2, using the two-dimensional profile line L1 of the original geological surface and the updated two-dimensional profile line L2, calculate the path of the characteristic line F and the corresponding elevation according to linear interpolation;

[0055] Step S3, setting the sampling range radius R1 of the sampling points around the characteristic line F that need to participate in the interpolation calculation, the selection range radius R2 of the points around the characteristic line F that need to re-interpolate the elevation calculation, and the reference elevation point range radius R3 around the characteristic line F that need to participate in the interpolation calculation; wherein R1>R2, R3>R2;

[0056] Step S4, setting the re-interpolation calculation spacing S of the feature line F;

[0057] Step S5: determine whether the single-segment feature line F is too long. If the length of the single-segment feature line exceeds a preset length threshold, the single-segment feature line is split;

[0058] Step S6, sequentially obtain the feature line outer bounding box BOX of a feature line segment, calculate the feature line sampling point outer bounding box BOX1 of the sampling point set P1 of the feature line segment F, and the feature line recalculation point outer bounding box BOX2 of the calculation point set P2 of the feature line segment F that needs to recalculate the elevation;

[0059] Step S7, determining the sampling point set P1, calculation point set P2 and reference elevation point set P3 of the characteristic line segment;

[0060] Step S8, merging the reference elevation point set P3 of the feature line segment into its sampling point set P1;

[0061] Step S9, using Kriging interpolation to recalculate the new elevations of the calculation points in the calculation point set P2 of the segment of the characteristic line according to the sampling point set P1 of the segment of the characteristic line and the set interpolation calculation interval S;

[0062] Step S10, looping through steps S6 to S9 until all feature lines are traversed;

[0063] Step S11: update the mesh surface according to the recalculated new elevation and regenerate the profile.

[0064] According to a specific embodiment of the present invention, step S1 specifically includes the following steps:

[0065] Select a sub-entity from the geological two-dimensional profile, determine the original geological surface to be modified, the profile line L1 corresponding to the geological surface, the corresponding characteristic line F, and the elevation set E1 corresponding to each key point on the characteristic line F that meets the preset rules;

[0066] A polyline L2 is drawn on the geological two-dimensional section map to replace the section line in the original geological two-dimensional section map, to obtain an updated two-dimensional section line L2, and the updated two-dimensional section line L2 is selected.

[0067] According to a specific embodiment of the present invention, step S2 specifically includes the following steps: restoring the plane coordinates of the characteristic line F from the geological two-dimensional profile, and establishing a first key point set based on the key points on the characteristic line F that meet the preset rules; projecting the selected updated two-dimensional profile line L2 into the geological two-dimensional profile map, establishing a second key point set based on the key points on the updated two-dimensional profile line L2 that meet the preset rules, and using the projection relationship to calculate the plane coordinates and new elevation set E2 corresponding to the second key point set; inserting the key points in the first key point set that do not have a corresponding relationship with the key points of the second key point set into the second key point set, and using linear interpolation to calculate the corresponding elevation and add it to the corresponding position of the new elevation set E2.

[0068] According to a specific embodiment of the present invention, a Boolean AND operation is performed on the first key point set and the second key point set, and key points in the first key point set that do not correspond to key points in the second key point set are inserted into the second key point set. The specific method is as follows:

[0069] The key points in the first key point set are sequentially found in the second key point set according to their positions on the feature lines. If they exist, they are ignored. If they do not exist, the key points in the first key point set are inserted into the second key point set.

[0070] According to a specific embodiment of the present invention, step S9 specifically includes the following steps: adding the reference elevation point set P3 to the sampling point set P1, using the points in the reference elevation point set P3 to cover the elevations of the points in the sampling point set P1 for points with repeated plane coordinates, and adding the new elevation point set E2 corresponding to the characteristic line to the sampling point set P1 as interpolation sampling points.

[0071] According to a specific embodiment of the present invention, the preset length threshold in step S5 is greater than or equal to 4 times the sampling range radius R1. When the length of a single-segment feature line exceeds the preset length threshold, a new key point is directly inserted into the feature line segment using a linear interpolation method as the endpoint of the split line segment.

[0072] According to a specific embodiment of the present invention, step S7 includes the following steps:

[0073] According to the two vertices of the single-segment feature line, calculate an original feature line outer bounding box BOX;

[0074] According to the setting value of R1, the original bounding box BOX is expanded to a range with a radius of R1 to form the outer bounding box BOX1 of the feature line sampling point;

[0075] According to the setting value of R2, the original bounding box is expanded to a range with a radius of R2 to form a feature line and recalculate the outer bounding box BOX2 of the point.

[0076] According to a specific embodiment of the present invention, step S8 specifically includes the following steps:

[0077] Filter out the points within the bounding box BOX1 outside the feature line sampling point, and then filter out the sampling point set P1 in this range by the distance of the feature line;

[0078] Filter out the feature lines and recalculate the points within the bounding box BOX2, and then filter out the calculation point set P2 in this range by the distance of the feature lines;

[0079] The reference elevation points within the bounding box BOX1 are filtered out through the feature line sampling point, and then the reference elevation point set P3 within the range is filtered out through the distance to the feature line.

[0080] According to a specific embodiment of the present invention, step S11 specifically includes the following steps: updating the local surface of the grid according to the recalculated calculation point set P2, and regenerating the geological section line in the cross-section diagram at the original position.

[0081] The present invention also provides a surface reconstruction modeling system based on section plane editing, comprising a processor capable of executing a computer program, wherein the computer program can implement the above-mentioned surface reconstruction modeling method based on section plane editing.

[0082] Example 1

[0083] like Figure 11 As shown, in step S1, the two-dimensional section line L1 of the original geological surface and the updated two-dimensional section line L2 are selected, and the updated two-dimensional section line is obtained by directly drawing a polyline on the section diagram to replace the original geological section line;

[0084] Step S2, using the two-dimensional profile line L1 of the original geological surface and the updated two-dimensional profile line L2, calculate the path of the characteristic line F and the corresponding elevation according to linear interpolation;

[0085] Step S3, setting the sampling range radius R1 of the sampling points around the characteristic line F that need to participate in the interpolation calculation, the selection range radius R2 of the points around the characteristic line F that need to re-interpolate the elevation calculation, and the reference elevation point range radius R3 around the characteristic line F that need to participate in the interpolation calculation; wherein R1>R2, R3>R2;

[0086] Step S4, setting the re-interpolation calculation spacing S of the feature line F;

[0087] Step S5: determine whether the single-segment feature line F is too long. If the length of the single-segment feature line exceeds a preset length threshold, the single-segment feature line is split;

[0088] Step S6, sequentially obtain the feature line outer bounding box BOX of a feature line segment, calculate the feature line sampling point outer bounding box BOX1 of the sampling point set P1 of the feature line segment F, and the feature line recalculation point outer bounding box BOX2 of the calculation point set P2 of the feature line segment F that needs to recalculate the elevation;

[0089] Step S7, determining the sampling point set P1, calculation point set P2 and reference elevation point set P3 of the characteristic line segment;

[0090] Step S8, merging the reference elevation point set P3 of the feature line segment into its sampling point set P1;

[0091] Step S9, using Kriging interpolation to recalculate the new elevations of the calculation points in the calculation point set P2 of the segment of the characteristic line according to the sampling point set P1 of the segment of the characteristic line and the set interpolation calculation interval S;

[0092] Step S10, looping through steps S6 to S9 until all feature lines are traversed;

[0093] Step S11: update the mesh surface according to the recalculated new elevation and regenerate the profile.

[0094] Example 2

[0095] like Figure 1-11 As shown, in step S1, select the original geological surface (such as Figure 6 The two-dimensional section line L1 (as shown) Figure 5 As shown) and the updated two-dimensional section line L2, the updated two-dimensional section line is obtained by directly drawing a polyline on the section diagram to replace the original geological section line;

[0096] Step S2, using the two-dimensional profile line L1 of the original geological surface and the updated two-dimensional profile line L2, calculate the path of the characteristic line F and the corresponding elevation according to linear interpolation;

[0097] Step S3, set the sampling range radius R1 of the sampling points around the characteristic line F that need to participate in the interpolation calculation, the selection range radius R2 of the points around the characteristic line F that need to re-interpolate the elevation, and the range radius R3 of the reference elevation points around the characteristic line F that need to participate in the interpolation calculation; wherein R1>R2, R3>R2; wherein in this embodiment, the reference elevation point is a drill hole, a flat hole, etc. Figure 1-4 As shown;

[0098] Step S4, setting the re-interpolation calculation spacing S of the feature line F;

[0099] Step S5: determine whether the single-segment feature line F is too long. If the length of the single-segment feature line exceeds a preset length threshold, the single-segment feature line is split;

[0100] Step S6, sequentially obtain the feature line outer bounding box BOX of a feature line segment, calculate the feature line sampling point outer bounding box BOX1 of the sampling point set P1 of the feature line segment F, and the feature line recalculation point outer bounding box BOX2 of the calculation point set P2 of the feature line segment F that needs to recalculate the elevation;

[0101] Step S7, determining the sampling point set P1, calculation point set P2 and reference elevation point set P3 of the characteristic line segment;

[0102] Step S8, merging the reference elevation point set P3 of the feature line segment into its sampling point set P1;

[0103] Step S9, using Kriging interpolation to recalculate the new elevations of the calculation points in the calculation point set P2 of the segment of the characteristic line according to the sampling point set P1 of the segment of the characteristic line and the set interpolation calculation interval S;

[0104] Step S10, looping through steps S6 to S9 until all feature lines are traversed;

[0105] Step S11, based on the recalculated new elevation, update the mesh surface and regenerate the profile. Figure 7-10 The following figure shows the resulting profiles after selecting different R1 and R2 values. The dark gray areas in the figure represent points within the R2 range, while the white areas represent the characteristic lines. Larger values ​​increase the range of variation and allow for greater scope for modification. You can select different values ​​as needed.

[0106] In this embodiment, step S1 specifically includes the following steps:

[0107] Select a sub-entity from the geological two-dimensional profile, determine the original geological surface to be modified, the profile line L1 corresponding to the geological surface, the corresponding characteristic line F, and the elevation set E1 corresponding to each key point on the characteristic line F that meets the preset rules;

[0108] A polyline L2 is drawn on the geological two-dimensional section map to replace the section line in the original geological two-dimensional section map, to obtain an updated two-dimensional section line L2, and the updated two-dimensional section line L2 is selected.

[0109] In this embodiment, step S2 specifically includes the following steps: restoring the plane coordinates of the characteristic line F from the geological two-dimensional profile, and establishing a first key point set based on the key points on the characteristic line F that meet the preset rules; projecting the selected updated two-dimensional profile line L2 into the geological two-dimensional profile, establishing a second key point set based on the key points on the updated two-dimensional profile line L2 that meet the preset rules, and using the projection relationship to calculate the plane coordinates and new elevation set E2 corresponding to the second key point set; inserting the key points in the first key point set that do not have a corresponding relationship with the key points of the second key point set into the second key point set, and using linear interpolation to calculate the corresponding elevation and add it to the corresponding position of the new elevation set E2.

[0110] In this embodiment, a Boolean AND operation is performed on the first key point set and the second key point set, and key points in the first key point set that do not correspond to key points in the second key point set are inserted into the second key point set. The specific method is as follows:

[0111] The key points in the first key point set are sequentially found in the second key point set according to their positions on the feature lines. If they exist, they are ignored. If they do not exist, the key points in the first key point set are inserted into the second key point set.

[0112] In this embodiment, step S9 specifically includes the following steps: adding the reference elevation point set P3 to the sampling point set P1; for points with repeated plane coordinates, using the points in the reference elevation point set P3 to cover the elevations of the points in the sampling point set P1; and adding the new elevation point set E2 corresponding to the characteristic line to the sampling point set P1 as interpolation sampling points.

[0113] In this embodiment, the preset length threshold in step S5 is greater than or equal to 4 times the sampling range radius R1. When the length of a single feature line segment exceeds the preset length threshold, a new key point is directly inserted into the feature line segment using the linear interpolation method as the endpoint of the split segment.

[0114] In this embodiment, step S7 includes the following steps:

[0115] According to the two vertices of the single-segment feature line, calculate an original feature line outer bounding box BOX;

[0116] According to the setting value of R1, the original bounding box BOX is expanded to a range with a radius of R1 to form the outer bounding box BOX1 of the feature line sampling point;

[0117] According to the setting value of R2, the original bounding box is expanded to a range with a radius of R2 to form a feature line and recalculate the outer bounding box BOX2 of the point.

[0118] In this embodiment, step S8 specifically includes the following steps:

[0119] Filter out the points within the bounding box BOX1 outside the feature line sampling point, and then filter out the sampling point set P1 in this range by the distance of the feature line;

[0120] Filter out the feature lines and recalculate the points within the bounding box BOX2, and then filter out the calculation point set P2 in this range by the distance of the feature lines;

[0121] The reference elevation points within the bounding box BOX1 are filtered out through the feature line sampling point, and then the reference elevation point set P3 within the range is filtered out through the distance to the feature line.

[0122] In this embodiment, step S11 specifically includes the following steps: updating the local surface of the grid according to the recalculated calculation point set P2, and regenerating the geological section line in the cross-section diagram at the original position.

[0123] The foregoing description is merely a preferred embodiment of the present invention and is not intended to limit the present invention. Those skilled in the art will readily appreciate that various modifications and variations of the present invention are possible. Any modifications, equivalent substitutions, or improvements made within the spirit and principles of the present invention are intended to be within the scope of protection of the present invention.

Claims

1. A surface reconstruction modeling method based on section editing, characterized in that: The steps include: Step S1, selecting the two-dimensional section line L1 of the original geological surface and the updated two-dimensional section line L2, wherein the updated two-dimensional section line L2 is obtained by directly drawing a polyline on the section diagram to replace the original geological section line; Step S2, using the two-dimensional profile line L1 of the original geological surface and the updated two-dimensional profile line L2, calculate the path of the characteristic line F and the corresponding elevation according to linear interpolation; Step S3, setting the sampling range radius R1 of the sampling points around the characteristic line F that need to participate in the interpolation calculation, the selection range radius R2 of the points around the characteristic line F that need to re-interpolate the elevation calculation, and the reference elevation point range radius R3 around the characteristic line F that need to participate in the interpolation calculation; wherein R1>R2, R3>R2; Step S4, setting the re-interpolation calculation spacing S of the feature line F; Step S5: determine whether the single-segment feature line F is too long. If the length of the single-segment feature line exceeds a preset length threshold, the single-segment feature line is split; Step S6, sequentially obtain the feature line outer bounding box BOX of a feature line segment, calculate the feature line sampling point outer bounding box BOX1 of the sampling point set P1 of the feature line segment F, and the feature line recalculation point outer bounding box BOX2 of the calculation point set P2 of the feature line segment F that needs to recalculate the elevation; Step S7, determining the sampling point set P1, calculation point set P2 and reference elevation point set P3 of the characteristic line segment; Step S8, merging the reference elevation point set P3 of the feature line segment into its sampling point set P1; Step S9, using Kriging interpolation to recalculate the new elevations of the calculation points in the calculation point set P2 of the segment of the characteristic line according to the sampling point set P1 of the segment of the characteristic line and the set interpolation calculation interval S; Step S10, looping through steps S6 to S9 until all feature lines are traversed; Step S11: update the mesh surface according to the recalculated new elevation and regenerate the profile.

2. The surface reconstruction modeling method based on section editing according to claim 1, characterized in that: Step S1 specifically includes the following steps: Select a sub-entity from the geological two-dimensional profile, determine the original geological surface to be modified, the profile line L1 corresponding to the geological surface, the corresponding characteristic line F, and the elevation set E1 corresponding to each key point on the characteristic line F that meets the preset rules; A polyline L2 is drawn on the geological two-dimensional section map to replace the section line in the original geological two-dimensional section map, to obtain an updated two-dimensional section line L2, and the updated two-dimensional section line L2 is selected.

3. The surface reconstruction modeling method based on section editing according to claim 2, characterized in that: Step S2 specifically includes the following steps: restoring the plane coordinates of the characteristic line F from the geological two-dimensional profile, establishing a first key point set based on the key points on the characteristic line F that meet the preset rules; projecting the selected updated two-dimensional profile line L2 onto the geological two-dimensional profile, establishing a second key point set based on the key points on the updated two-dimensional profile line L2 that meet the preset rules, and calculating the plane coordinates and new elevation set E2 corresponding to the second key point set using the projection relationship; The key points in the first key point set that do not correspond to the key points in the second key point set are inserted into the second key point set, and the corresponding elevations are calculated using linear interpolation and added to the corresponding positions in the new elevation set E2.

4. The surface reconstruction modeling method based on section editing according to claim 3 is characterized in that: Perform a Boolean AND operation on the first key point set and the second key point set, and insert the key points in the first key point set that do not correspond to the key points in the second key point set into the second key point set. The specific method is as follows: The key points in the first key point set are sequentially found in the second key point set according to their positions on the feature lines. If they exist, they are ignored. If they do not exist, the key points in the first key point set are inserted into the second key point set.

5. The surface reconstruction modeling method based on section editing according to claim 3 is characterized in that: Step S9 specifically includes the following steps: adding the reference elevation point set P3 to the sampling point set P1; for points with repeated plane coordinates, using the points in the reference elevation point set P3 to cover the elevations of the points in the sampling point set P1; and adding the new elevation point set E2 corresponding to the characteristic line to the sampling point set P1 as interpolation sampling points.

6. The surface reconstruction modeling method based on section editing according to claim 1, characterized in that: In step S5, the preset length threshold is greater than or equal to 4 times the sampling range radius R1. When the length of a single feature line segment exceeds the preset length threshold, a new key point is directly inserted into the feature line segment using the linear interpolation method as the endpoint of the split segment.

7. The surface reconstruction modeling method based on section editing according to any one of claims 1 to 6, characterized in that: Step S7 includes the following steps: According to the two vertices of the single-segment feature line, calculate an original feature line outer bounding box BOX; According to the setting value of R1, the original bounding box BOX is expanded to a range with a radius of R1 to form the outer bounding box BOX1 of the feature line sampling point; According to the setting value of R2, the original bounding box is expanded to a range with a radius of R2 to form a feature line and recalculate the outer bounding box BOX2 of the point.

8. The surface reconstruction modeling method based on section editing according to claim 1, characterized in that: Step S8 specifically includes the following steps: Filter out the points within the bounding box BOX1 outside the feature line sampling point, and then filter out the sampling point set P1 in this range by the distance of the feature line; Filter out the feature lines and recalculate the points within the bounding box BOX2, and then filter out the calculation point set P2 in this range by the distance of the feature lines; The reference elevation points within the bounding box BOX1 are filtered out through the feature line sampling point, and then the reference elevation point set P3 within the range is filtered out through the distance to the feature line.

9. The surface reconstruction modeling method based on section editing according to claim 1, characterized in that: Step S11 specifically includes the following steps: updating the local surface of the grid according to the recalculated calculation point set P2, and regenerating the geological section line in the cross-section diagram at the original position.

10. A surface reconstruction modeling system based on section editing, characterized in that: The invention comprises a processor, wherein the processor is capable of executing a computer program, and the computer program can implement the surface reconstruction modeling method based on section plane editing according to any one of claims 1 to 9.

Images