A surface twin scene modeling method and system based on exploration digital results

Abstract

The application provides a kind of surface twin scene modeling method and system based on exploration digital result, belong to engineering geological exploration and geographic information system technical field.The application generates real three-dimensional terrain grid model using exploration CAD terrain surface data and orthographic image;Corresponding river, road and building three-dimensional grid model is generated using exploration CAD river, road and building boundary data;Polygonal area generated by river, road and building boundary data is divided using ear cutting method;River and road grid model are refined to adapt to the grid spacing of terrain grid model;For the part of terrain grid model protruding river and road grid model, appropriate correction is carried out;Finally, map material is given to river, road and building model, and finally three-dimensional digital twin scene is generated.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of engineering geological exploration and geographic information system, and particularly relates to a surface twin scene modeling method and system based on exploration digital results. BACKGROUND

[0002] In today's many industry fields, such as water conservancy and hydropower, natural resource monitoring, precision agriculture, ecological restoration engineering, and large infrastructure surrounding environment simulation, the demand for constructing a digital twin scene according to field exploration information is increasingly urgent. However, the traditional modeling method has many shortcomings when facing these practical application scenarios:

[0003] Insufficient use of exploration result data: Exploration result data is mostly saved in the form of two-dimensional CAD drawings, and the delivered results are mostly in the form of two-dimensional results. The form of geological engineering exploration results is relatively backward: it is not possible to quickly build a twin scene of the surface in the exploration area according to the completed digital exploration results. In view of the above difficulties, it is urgent to research a surface twin scene modeling method based on exploration digital results to fill the gap in existing technology and meet the needs of the industry. SUMMARY

[0004] In view of the deficiencies in the prior art, the present application provides a surface twin scene modeling method and system based on exploration digital results, aiming to provide a more realistic and accurate surface twin scene that meets the exploration results for the water conservancy, water conservancy, highway, railway, petroleum and other industries as a simulation simulation, auxiliary design base. The present application comprises: generating a real three-dimensional terrain mesh model using exploration CAD terrain surface data and orthophoto; generating corresponding river, road and building three-dimensional mesh models using exploration CAD river, road and building boundary data; using ear cutting method to subdivide the polygon area generated by the river, road and building boundary data; refining the river and road mesh models to adapt to the grid spacing of the terrain mesh model; correcting the part of the terrain mesh model protruding the river and road mesh model; finally, assigning a map material to the river, road and building model, and finally generating a three-dimensional digital twin scene. The modeling process of the present application is convenient for programmatic automatic generation, suitable for the workflow of mainstream rendering engines, and the data form is adaptive.

[0005] To achieve the above technical purposes, the present application adopts the following technical solutions:

[0006] A surface twin scene modeling method based on exploration digital results, comprising the following steps:

[0007] Step S1, generating a three-dimensional terrain surface model using exploration CAD terrain surface data and aerial orthographic image data, reconstructing the terrain surface using a spatial interpolation algorithm, adopting vertex-triangle face-grid description, and orthographic image as texture mapping to the terrain grid model;

[0008] Step S2, generating and refining a three-dimensional river model using exploration CAD river system boundary line data, adopting a polygon triangulation algorithm, and adopting point-triangle face-grid description;

[0009] Step S3, generating and refining a three-dimensional road model using exploration CAD road boundary line data, adopting a polygon triangulation algorithm, and adopting point-triangle face-grid description;

[0010] Step S4, generating a three-dimensional building model using exploration CAD building boundary data, adopting a polygon triangulation algorithm, and adopting point-triangle face-grid description;

[0011] Step S5, modifying the terrain grid model generated in step S1 to adapt to the three-dimensional river model generated in step S2 and the three-dimensional road model generated in step S3;

[0012] Step S6, making river material, road material and building material, respectively assigning to the three-dimensional river model, the three-dimensional road model and the three-dimensional building model, and obtaining a three-dimensional digital twin model.

[0013] Preferably, the specific steps of generating the three-dimensional terrain surface model include:

[0014] Step S1.1, collecting exploration digital result terrain surface vertex-face-grid data;

[0015] Step S1.2, calculating the terrain surface range, setting the grid spacing, and forming a plane grid under X-Y coordinates;

[0016] Step S1.3, calculating the elevation value of the grid point using a spatial interpolation method;

[0017] Step S1.4, cutting the orthographic image according to the terrain surface range;

[0018] Step S1.5, cutting the orthographic image as texture mapping to the terrain grid model;

[0019] Step S1.6, generating a three-dimensional terrain surface model.

[0020] Preferably, the spatial interpolation method is inverse distance weighted average method.

[0021] Preferably, the three-dimensional terrain surface model grid adopts exploration digital result terrain surface data, and the texture adopts aerial orthographic image.

[0022] Preferably, the river, the road and the building adopt the river, road and building boundary data corresponding to the exploration digital results.

[0023] Preferably, the specific steps of the polygonal triangulation algorithm for generating and refining the three-dimensional river model include:

[0024] Step S2.1, collecting the boundary vertex three-dimensional coordinate data of any river in the study area;

[0025] Step S2.2, for the river, using the polygonal triangulation algorithm to triangulate the boundary vertex three-dimensional coordinate data of the river to generate a triangular mesh model;

[0026] Step S2.3, using a mesh subdivision algorithm to encrypt the generated triangular mesh model, dividing each triangular face into four, repeating the subdivision, so that the edges of the triangular faces are smaller than the edges of the terrain grid triangular faces;

[0027] Step S2.4, obtaining the three-dimensional mesh model of the river, if the three-dimensional river model of each river has been generated, then end, otherwise go to step S2.1.

[0028] Preferably, the specific steps of the polygonal triangulation algorithm for generating and refining the three-dimensional road model include:

[0029] Step S3.1, collecting the boundary vertex three-dimensional coordinate data of any road in the study area;

[0030] Step S3.2, for the road, using the polygonal triangulation algorithm to triangulate the boundary vertex three-dimensional coordinate data of the road to generate a triangular mesh model;

[0031] Step S3.3, using a mesh subdivision algorithm to encrypt the generated triangular mesh model, dividing each triangular face into four, repeating the subdivision, so that the edges of the triangular faces are smaller than the terrain grid spacing;

[0032] Step S3.4, obtaining the three-dimensional mesh model of the road, if the three-dimensional road model of each road has been generated, then end, otherwise go to step S3.1.

[0033] Preferably, the specific steps of the polygonal triangulation algorithm for generating a three-dimensional building model are:

[0034] Step S4.1, collecting the boundary vertex three-dimensional coordinate data and building height data of any building in the study area;

[0035] Step S4.2, for the building, generating a building facade polygon along the building boundary point according to the height, and finally generating a top surface polygon of the modeling object according to the height;

[0036] Step S4.3, using a polygon triangulation algorithm to divide the facade and roof of the current building to generate a triangular mesh model;

[0037] Step S4.4, obtain the three-dimensional mesh model of the building, if the three-dimensional building model of each building has been generated, end, otherwise go to step S4.1.

[0038] Preferably, the step S5 integrates and corrects the terrain mesh model to adapt to the specific steps of the three-dimensional river model and the three-dimensional road model:

[0039] Step S5.1, collect the three-dimensional river model and the three-dimensional road model generated in steps S2 and S3;

[0040] Step S5.2, calculate the bounding box of each three-dimensional river model and three-dimensional road model;

[0041] Step S5.3, for the terrain mesh model generated in step S1, calculate the terrain mesh points contained in the bounding box generated in step S5.2;

[0042] Step S5.4, from the terrain mesh point as the starting point, the direction is downward, judge whether it collides with the mesh in the current three-dimensional river model or three-dimensional road model, if it collides, the corresponding vertex of the terrain mesh is concave downward by a specified depth;

[0043] Step S5.5, repeat steps 5.2 to 5.4 until the correction of each three-dimensional road model and three-dimensional river model to the terrain mesh model generated in step S1 is completed.

[0044] Preferably, the polygon triangulation algorithm used in steps S2, S3 and S4 is ear cutting method, and the subdivision algorithm is Loop subdivision algorithm.

[0045] Preferably, the step S6 makes river material, road material and building material, and the specific steps of assigning three-dimensional river model, three-dimensional road model and three-dimensional building model are as follows:

[0046] Step S6.1, collect river, road and building surface texture map images;

[0047] Step S6.2, make river, road and building material shaders according to the map images;

[0048] Step S6.3, assign river, road and building material shaders to three-dimensional river model, three-dimensional road model and three-dimensional building model to obtain a three-dimensional digital twin model.

[0049] The application further provides a surface twin scene modeling system based on exploration digital results.

[0050] Compared with the prior art, the application has the following beneficial effects:

[0051] (1) The application restores the real scene with high precision: through the high-precision two-dimensional drawing data of geological engineering exploration, the terrain, water system, road and building model highly consistent with the real environment can be generated, thereby providing reliable digital twin scene for various engineering field applications.

[0052] (2) The application generates a twin scene model with high efficiency and automation: the traditional manual modeling method is abandoned, and the programmatic automatic generation method is adopted, so that the digital twin scene model can be quickly generated according to the engineering exploration data in a short time, the project cycle is greatly shortened, the labor cost is reduced, and the application is especially suitable for the modeling task of large-area and complex scene in the engineering field. BRIEF DESCRIPTION OF DRAWINGS

[0053] Figure 1 The figure is a flowchart of the surface twin scene modeling method based on exploration digital results of the preferred embodiment of the application.

[0054] Figure 2 The figure is a generation schematic diagram of the three-dimensional terrain surface model of the preferred embodiment of the application.

[0055] Figure 3 The figure is a generation schematic diagram of the three-dimensional river model and the three-dimensional road model of the preferred embodiment of the application.

[0056] Figure 4 The figure is a generation schematic diagram of the three-dimensional building model of the preferred embodiment of the application.

[0057] Figure 5 The figure is a process schematic diagram of the three-dimensional terrain correction to adapt to the river, road and building of the preferred embodiment of the application.

[0058] Figure 6 The figure is a schematic diagram of ear cutting method triangulation of the embodiment of the application.

[0059] Figure 7 The figure is a schematic diagram of Loop subdivision algorithm refining the triangular mesh model of the embodiment of the application.

[0060] Figure 8 The figure is a result map of the generation of the three-dimensional terrain surface model of the embodiment of the application. DETAILED DESCRIPTION

[0061] With reference to the drawings and the embodiments of the present application, the technical solutions in the embodiments of the present application will be described clearly and completely. Obviously, the described embodiments are only a part of the embodiments of the present application, rather than all the embodiments of the present application. Based on the embodiments of the present application, all other embodiments obtained by those of ordinary skill in the art without creative work fall within the scope of the present application.

[0062] The present application provides a surface twin scene modeling method based on exploration digital results, comprising the following steps:

[0063] Step S1, generating a three-dimensional terrain surface model using exploration CAD terrain surface data and aerial orthographic image data, reconstructing the terrain surface using a spatial interpolation algorithm, and adopting vertex-triangle face-grid description, with the orthographic image being mapped as texture onto the terrain grid model;

[0064] Step S2, generating and refining a three-dimensional river model using exploration CAD river system boundary line data, adopting a polygon triangulation algorithm, and adopting point-triangle face-grid description;

[0065] Step S3, generating and refining a three-dimensional road model using exploration CAD road boundary line data, adopting a polygon triangulation algorithm, and adopting point-triangle face-grid description;

[0066] Step S4, generating a three-dimensional building model using exploration CAD building boundary data, adopting a polygon triangulation algorithm, and adopting point-triangle face-grid description;

[0067] Step S5, modifying the terrain grid model generated in step S1 to adapt to the three-dimensional river model generated in step S2 and the three-dimensional road model generated in step S3;

[0068] Step S6, making river material, road material and building material, respectively giving the three-dimensional river model, the three-dimensional road model and the three-dimensional building model, and obtaining a three-dimensional digital twin model.

[0069] According to a specific embodiment of the present application, the specific steps of generating the three-dimensional terrain surface model comprise:

[0070] Step S1.1, collecting exploration digital result terrain surface vertex-face-grid data;

[0071] Step S1.2, calculating the terrain surface range, setting the grid spacing, and forming a plane grid under X-Y coordinates;

[0072] Step S1.3, calculating the elevation value of the grid point using a spatial interpolation method;

[0073] Step S1.4, clipping the orthographic image according to the terrain surface range;

[0074] Step S1.5, mapping the clipped orthographic image as texture onto the terrain mesh model;

[0075] Step S1.6, generating the three-dimensional terrain surface model.

[0076] According to one specific embodiment of the present application, the spatial interpolation method is inverse distance weighted average method.

[0077] According to one specific embodiment of the present application, the three-dimensional terrain surface model mesh adopts the exploration digital result terrain surface data, and the texture adopts the aerial orthographic image.

[0078] According to one specific embodiment of the present application, the river, the road and the building adopt the corresponding river, road and building boundary data of the exploration digital result.

[0079] According to one specific embodiment of the present application, the specific steps of the polygonal triangulation algorithm for generating and refining the three-dimensional river model include:

[0080] Step S2.1, collecting the boundary vertex three-dimensional coordinate data of any river in the study area;

[0081] Step S2.2, for the river, using the polygonal triangulation algorithm to triangulate the boundary vertex three-dimensional coordinate data of the river to generate a triangular mesh model;

[0082] Step S2.3, using a mesh subdivision algorithm to encrypt the generated triangular mesh model, dividing each triangular face into four, repeating the subdivision, so that the edges of the triangular face are smaller than the edges of the terrain mesh triangular face;

[0083] Step S2.4, obtaining the three-dimensional mesh model of the river, if the three-dimensional river model of each river has been generated, then ending, otherwise turning to step S2.1.

[0084] According to one specific embodiment of the present application, the specific steps of the polygonal triangulation algorithm for generating and refining the three-dimensional road model include:

[0085] Step S3.1, collecting the boundary vertex three-dimensional coordinate data of any road in the study area;

[0086] Step S3.2, for the road, using the polygonal triangulation algorithm to triangulate the boundary vertex three-dimensional coordinate data of the road to generate a triangular mesh model;

[0087] Step S3.3, using a mesh subdivision algorithm to encrypt the generated triangular mesh model, dividing each triangle into four, repeating the subdivision, so that the edge of the triangle is smaller than the terrain grid spacing;

[0088] Step S3.4, obtaining the three-dimensional grid model of the road, if the three-dimensional road model of each road has been generated, then end, otherwise go to step S3.1.

[0089] According to an embodiment of the present application, the specific steps of the polygon triangulation algorithm for generating a three-dimensional building model are:

[0090] Step S4.1, collecting the boundary vertex three-dimensional coordinate data and building height data of any building in the study area;

[0091] Step S4.2, for the building, generating building facade polygons along the building boundary points according to height, and finally generating top surface polygons of the modeling object according to height;

[0092] Step S4.3, using a polygon triangulation algorithm to subdivide the facade and top surface of the current building to generate a triangular mesh model;

[0093] Step S4.4, obtaining the three-dimensional grid model of the building, if the three-dimensional building model of each building has been generated, then end, otherwise go to step S4.1.

[0094] According to an embodiment of the present application, the specific steps of step S5 for integrating and correcting the terrain grid model to adapt to the three-dimensional river model and the three-dimensional road model are:

[0095] Step S5.1, collecting the three-dimensional river model and the three-dimensional road model generated in steps S2 and S3;

[0096] Step S5.2, calculating the bounding box of the boundary of each three-dimensional river model and three-dimensional road model;

[0097] Step S5.3, for the terrain grid model generated in step S1, calculating the terrain grid points contained in the bounding box generated in step S5.2;

[0098] Step S5.4, starting from the terrain grid point, emitting a ray in the downward direction, and judging whether it collides with the grid in the current three-dimensional river model or three-dimensional road model, if it collides, then the corresponding vertex of the terrain grid is depressed downward by a specified depth;

[0099] Step S5.5, repeating steps 5.2 to 5.4 until the correction of the terrain grid model generated in step S1 by each three-dimensional road model and three-dimensional river model is completed.

[0100] According to a specific embodiment of the present application, the polygon triangulation algorithm used in the steps S2, S3 and S4 is ear cutting method, and the subdivision algorithm is Loop subdivision algorithm.

[0101] According to a specific embodiment of the present application, the specific steps of manufacturing the river material, the road material and the building material in the step S6 to give the three-dimensional river model, the three-dimensional road model and the three-dimensional building model are as follows:

[0102] Step S6.1, collecting the surface texture map images of the river, the road and the building;

[0103] Step S6.2, manufacturing the river, road and building material shaders according to the map images;

[0104] Step S6.3, giving the river, road and building material shaders to the three-dimensional river model, the three-dimensional road model and the three-dimensional building model to obtain the three-dimensional digital twin model.

[0105] The present application also provides a surface twin scene modeling system based on exploration digital results, comprising a processor capable of executing a computer program capable of realizing the above-mentioned surface twin scene modeling method based on exploration digital results.

[0106] Embodiment 1

[0107] The present application provides a surface twin scene modeling method based on exploration digital results, comprising the following steps:

[0108] Step S1, generating a three-dimensional terrain surface model by using exploration CAD terrain surface data and aerial orthographic image data, reconstructing the terrain surface by using a spatial interpolation algorithm, and adopting vertex-triangle face-grid description, and mapping the orthographic image as texture to the terrain grid model;

[0109] Step S2, generating and refining a three-dimensional river model by using exploration CAD river system boundary line data and adopting a polygon triangulation algorithm, and adopting point-triangle face-grid description;

[0110] Step S3, generating and refining a three-dimensional road model by using exploration CAD road boundary line data and adopting a polygon triangulation algorithm, and adopting point-triangle face-grid description;

[0111] Step S4, generating a three-dimensional building model by using exploration CAD building boundary data and adopting a polygon triangulation algorithm, and adopting point-triangle face-grid description;

[0112] Step S5, modifying the terrain grid model generated in step S1 to adapt to the three-dimensional river model generated in step S2 and the three-dimensional road model generated in step S3.

[0113] Step S6, making river material, road material and building material, respectively giving three-dimensional river model, three-dimensional road model and three-dimensional building model, obtaining three-dimensional digital twin model.

[0114] Embodiment 2

[0115] The application provides a surface twin scene modeling method based on exploration digital results, including the following steps:

[0116] (1) input exploration result CAD data and aerial orthographic image;

[0117] (2) read the terrain surface from the exploration result CAD data, reconstruct the terrain surface grid by using a spatial interpolation algorithm, adopt vertex-face-grid description, and cut the aerial orthographic image according to the range of the terrain surface grid model to serve as the terrain texture, Figure 2 The flow chart for establishing the three-dimensional terrain surface grid model is as follows:

[0118] A, read in the terrain surface data; read in the aerial orthographic image;

[0119] B, calculate the bounding box range of the terrain surface;

[0120] C, set the grid spacing to form a plane grid on the XY coordinate plane;

[0121] D, generate the elevation value at each plane grid intersection point by using a spatial interpolation algorithm, such as the inverse distance weighting interpolation method;

[0122] E, obtain the reconstructed terrain grid, and the expression form is vertex-face-grid;

[0123] F, cut the orthographic image according to the bounding box range of the terrain surface;

[0124] G, associate the orthographic image texture to the terrain grid to generate the three-dimensional terrain surface model, and the expression form is vertex-face-grid.

[0125] (3), read the river boundary data and road boundary data from the exploration result CAD data, generate the river grid model and road grid model by using the polygon triangulation ear cutting algorithm, and according to the set terrain grid spacing, the generated grid model is refined and encrypted, Figure 3 The flow chart for generating the three-dimensional river road grid model is as follows:

[0126] A, read in the river boundary data and road boundary data;

[0127] B, extract the three-dimensional coordinate data of the river and road boundary polygon vertexes;

[0128] C. using polygon ear clipping method to clip river and road boundary;

[0129] D. generating initial river and road grid model;

[0130] E. performing Loop grid subdivision on initial grid model until grid edge length is less than terrain grid spacing;

[0131] F. generating final river and road grid model, expressed as vertex-face-grid.

[0132] (4) reading building boundary data and building height data from exploration results CAD data, generating building side and top polygon according to height, and then using polygon ear clipping method to clip the generated building side and top polygon to generate building grid model, Figure 4 a flowchart for generating 3D building grid model, and the specific implementation process is as follows:

[0133] A. reading in building boundary and building height information;

[0134] B. extracting 3D coordinates of building boundary polygon vertices;

[0135] C. generating side elevation rectangle on each edge of building boundary polygon according to building height, and finally generating building top polygon;

[0136] D. using ear clipping method to triangulate the generated building side and top polygon;

[0137] E. generating building grid model, expressed as vertex-face-grid.

[0138] (5) correcting terrain grid model according to the generated river and road grid model, so that the river and road grid model is more suitable for the surface terrain grid model, Figure 5 a flowchart for 3D terrain correction to adapt to river, road and building, and the specific implementation process is as follows:

[0139] A. reading in the generated ground, river and road grid model;

[0140] B. starting to correct the ground grid model, if all river and road models have participated in terrain correction, then executing step H;

[0141] C. calculating the bounding box of the current river or road grid model;

[0142] D. calculating the terrain vertices contained in the current bounding box;

[0143] E. Shoot a ray downward from the containing terrain point, check if it collides with the current river and road model;

[0144] F. If collision, the ground point is depressed by a specified depth, otherwise, no modification is made to the terrain vertex;

[0145] G. Execute the next batch of river and road mesh model correction terrain mesh, execute step B;

[0146] H. End the correction of the ground terrain mesh model.

[0147] (5) Collect river texture map, road texture map and building texture, use 3D rendering engine editor to make river material, road material and building material, import the corrected terrain, river and building mesh model, specify the corresponding material, finally render into three-dimensional ground twin scene, the specific implementation process is as follows:

[0148] A. Collect river, road and building texture map;

[0149] B. Use image processing software or program to process the texture map, so that the map texture is clear without distortion and the like;

[0150] C. Make river, road and building material;

[0151] D. Import the corrected terrain, river and building mesh model;

[0152] E. Associate the material and the mesh model together;

[0153] F. Finally render into three-dimensional ground twin scene.

[0154] Figure 6 For the river, road and building boundary polygon ear cutting algorithm of the embodiment, the polygon area is divided into triangular grid; Figure 7 For the Loop subdivision mesh algorithm for refining the triangular face of the river and road mesh model, which aims to recursively encrypt the mesh model; Figure 8 For the reconstructed ground terrain mesh model, the schematic diagram of superimposed display of aerial orthographic image.

[0155] The above only describes the embodiments of the present application and does not limit the present application. Any modification, equivalent replacement and improvement within the scope of the application shall be included in the protection scope of the present application.

Claims

1. A method for modeling a surface twin scene based on exploration digital results, characterized in that, The method comprises the following steps: Step S1, generating a three-dimensional terrain surface model by using exploration CAD terrain curved surface data and aerial orthographic image data, reconstructing the terrain curved surface by using a spatial interpolation algorithm, adopting vertex-triangle face-grid description, and mapping the orthographic image as texture onto the terrain grid model; Step S2, generating and refining a three-dimensional river model by using exploration CAD river system boundary line data and a polygon triangle partitioning algorithm, adopting point-triangle face-grid description, and encrypting by using a grid subdivision algorithm during refinement, dividing each triangle face into four, and repeating the subdivision to make the edges of the triangle faces smaller than the edges of the terrain grid triangle faces; Step S3, generating and refining a three-dimensional road model by using exploration CAD road boundary line data and a polygon triangle partitioning algorithm, adopting point-triangle face-grid description, and encrypting by using a grid subdivision algorithm during refinement, dividing each triangle face into four, and repeating the subdivision to make the edges of the triangle faces smaller than the terrain grid spacing; Step S4, generating a three-dimensional building model by using exploration CAD building boundary data and a polygon triangle partitioning algorithm, adopting point-triangle face-grid description; Step S5, modifying the terrain grid model generated in step S1 to adapt to the three-dimensional river model generated in step S2 and the three-dimensional road model generated in step S3; Step S6, making river material, road material and building material, respectively assigning the three-dimensional river model, the three-dimensional road model and the three-dimensional building model, and obtaining a three-dimensional digital twin model; The specific steps of step S5 for integrating and modifying the terrain grid model to adapt to the three-dimensional river model and the three-dimensional road model are as follows: Step S5.1, collecting the three-dimensional river model and the three-dimensional road model generated in steps S2 and S3; Step S5.2, calculating the boundary bounding box of each three-dimensional river model and three-dimensional road model; Step S5.3, for the terrain grid model generated in step S1, calculating the terrain grid points contained in the bounding box generated in step S5.2; Step S5.4, starting from the terrain grid point, emitting a ray downward, and judging whether it collides with the grid in the current three-dimensional river model or three-dimensional road model, if it collides, the corresponding vertex of the terrain grid is depressed downward by a specified depth; Step S5.5, repeating steps 5.2 to 5.4 until the modification of the terrain grid model generated in step S1 by each three-dimensional road model and three-dimensional river model is completed.

2. The method of claim 1, wherein, The specific steps for generating the three-dimensional terrain surface model include: Step S1.1, collecting exploration digital result terrain curved surface vertex-face-grid data; Step S1.2, calculating the terrain curved surface range, setting the grid spacing, and forming a plane grid under X-Y coordinates; Step S1.3, calculating the elevation value of the grid point by using a spatial interpolation method; Step S1.4, cropping the orthographic image according to the terrain curved surface range; Step S1.5, mapping the cropped orthographic image as texture onto the terrain grid model; Step S1.6, generating a three-dimensional terrain surface model.

3. The method for modelling a surface twin scene based on exploration digital deliverables according to any one of claims 1 or 2, characterized in that, The spatial interpolation method is an inverse distance weighted average method.

4. The method of claim 1, wherein: The three-dimensional terrain surface model grid adopts exploration digital result terrain curved surface data, and the texture adopts aerial orthographic image.

5. The method of claim 1, wherein, The river, the road and the building adopt corresponding river, road and building boundary data of exploration digital result.

6. The method according to claim 5, wherein, The specific steps of the polygonal triangulation algorithm for generating and refining the three-dimensional river model include: Step S2.1, collecting boundary vertex three-dimensional coordinate data of any river in the study area; Step S2.2, for the river, using the polygonal triangulation algorithm to divide the boundary vertex three-dimensional coordinate data of the river to generate a triangular mesh model; Step S2.3, using a mesh subdivision algorithm to encrypt the generated triangular mesh model, dividing each triangular face into four, and repeating the subdivision so that the edges of the triangular faces are smaller than the edges of the terrain mesh triangles; Step S2.4, obtaining the three-dimensional mesh model of the river, if the three-dimensional river model of each river has been generated, then ending, otherwise turning to step S2.

1.

7. The method of claim 5, wherein, The specific steps of the polygonal triangulation algorithm for generating and refining the three-dimensional road model include: Step S3.1, collecting boundary vertex three-dimensional coordinate data of any road in the study area; Step S3.2, for the road, using the polygonal triangulation algorithm to divide the boundary vertex three-dimensional coordinate data of the road to generate a triangular mesh model; Step S3.3, using a mesh subdivision algorithm to encrypt the generated triangular mesh model, dividing each triangular face into four, and repeating the subdivision so that the edges of the triangular faces are smaller than the terrain mesh spacing; Step S3.4, obtaining the three-dimensional mesh model of the road, if the three-dimensional road model of each road has been generated, then ending, otherwise turning to step S3.

1.

8. The method of claim 5, wherein: The specific steps of the polygonal triangulation algorithm for generating the three-dimensional building model are: Step S4.1, collecting boundary vertex three-dimensional coordinate data and building height data of any building in the study area; Step S4.2, for the building, generating a building facade polygon along the building boundary point according to the height, and finally generating a top surface polygon of the modeling object according to the height; Step S4.3, using the polygonal triangulation algorithm to divide the facade and top surface of the current building to generate a triangular mesh model; Step S4.4, obtaining the three-dimensional mesh model of the building, if the three-dimensional building model of each building has been generated, then ending, otherwise turning to step S4.

1.

9. The method of claim 1, wherein: The polygonal triangulation algorithm used in steps S2, S3 and S4 is the ear cutting method, and the subdivision algorithm is the Loop subdivision algorithm.

10. The method of claim 1, wherein: The specific steps of the step S6 for making river material, road material and building material, and assigning the three-dimensional river model, the three-dimensional road model and the three-dimensional building model are: Step S6.1, collecting river, road and building surface texture map images; Step S6.2, making river, road and building material shaders according to the map images; Step S6.3, assigning the river, road and building material shaders to the three-dimensional river model, the three-dimensional road model and the three-dimensional building model to obtain a three-dimensional digital twin model.

11. A system for modeling a surface twin scene based on exploration digital deliverables, characterized in that, A computer program product comprising a computer readable medium having stored thereon the computer program of claim 11. A computer program product comprising a computer readable medium having stored thereon the computer program of claim 11.

Images