A site vertical design earthwork volume calculation method and system based on GIS

Abstract

The application discloses a kind of based on GIS's site vertical design earthwork quantity calculation method and system, it is related to earthwork quantity calculation technical field, including: according to topographic map, construct site present situation DEM;According to vertical design drawing, construct site design DEM;Site design DEM and present situation DEM are carried out grid analysis, obtain surface fill and dig earthwork grid data;Combined with control regulation data, calculate the surface fill and dig earthwork quantity of each plot of planning land;Calculate the underground engineering excavation quantity of each plot of planning land, building underground space, subway, underground comprehensive pipe gallery;The total amount of fill and dig of each plot of planning land is calculated, and fill and dig special map is generated.The system includes present situation DEM construction, design DEM construction, surface dig and fill quantity calculation, underground engineering excavation quantity calculation, total amount of dig and fill calculation 5 modules.The method is simple, efficient and accurate, low cost, wide application range, considering the underground engineering earth volume, the result is more reasonable and accurate, provide more reliable basis for site vertical design and construction.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of earthwork calculation, in particular to a site vertical design earthwork calculation method and system based on GIS. BACKGROUND

[0002] Earthwork calculation is an important part of site vertical design. Accurate and rapid earthwork calculation can more scientifically and reasonably determine the elevation relationship of each part of the construction site and arrange the construction sequence, which is of great significance to reducing engineering cost and accelerating construction progress. Traditional site vertical design earthwork calculation is mainly realized based on CAD earthwork calculation software, but this method is slow in calculation and is suitable for small-scale earthwork calculation. For large-scale earthwork calculation, the efficiency is not high, and it is necessary to use excel to calculate the earthwork of each plot in the project area. It is not convenient to use, and it does not consider and calculate the earthwork generated by underground space development and construction, resulting in inaccurate and unreasonable site earthwork calculation. Moreover, this method has poor visualization and is not intuitive, which is not conducive to guiding earthwork allocation.

[0003] In recent years, with the rapid development of spatial information technology, GIS, unmanned aerial photography and BIM technology have been applied to earthwork calculation.

[0004] For example, Chinese patent 202410199703.5 discloses a method and system for calculating earthwork based on unmanned aerial vehicles, which includes planning the flight path of the unmanned aerial vehicle, obtaining image data and point cloud data of the earthwork terrain, fusing the image data and point cloud data of the earthwork terrain and removing the occlusion, and calculating the earthwork using the fused data. Chinese patent 202311455540.4 discloses a method and system for earthwork allocation based on fusion of BIM and GIS, which establishes a GIS model by obtaining images of the project site through topographic survey, establishes a BIM model through construction drawings, fuses the two models, and considers the soil that does not meet the building requirements to calculate the earthwork. These earthwork calculation methods based on GIS technology, unmanned aerial photography and BIM technology have high calculation accuracy, but the data collection and production of laser point cloud and BIM model takes a long time and costs a lot, and the earthwork calculation process is relatively complicated. The most critical problem is that these methods only consider and calculate the changes of the ground surface of the project area, and do not consider the earthwork of underground space, subway and underground comprehensive pipe gallery, resulting in inaccurate and unreasonable earthwork calculation in site vertical design, which affects the scientificity of earthwork allocation and construction cost.

[0005] At present, there is no effective solution to the problems in the related art. SUMMARY

[0006] In view of the deficiencies of the prior art, the application provides a GIS-based site vertical design earthwork volume calculation method and system to overcome the above technical problems of the prior art, and more accurately, low-cost and simply calculate the earthwork volume in site vertical design.

[0007] In order to achieve the above-mentioned purpose, the application adopts the following technical solutions:

[0008] A GIS-based site vertical design earthwork volume calculation method, which comprises the following steps:

[0009] Step 1: Obtain a site status topographic map file, and construct a site status digital elevation model (DEM) according to the elevation points in the topographic map;

[0010] Step 2: Construct a site design digital elevation model (DEM) according to the design elevations in the site vertical design map, wherein the design elevations are derived from the control points marked in the site design map, and the control points include road intersection control points, road section control points, road outside land control points and control points at non-construction areas of the site;

[0011] Step 3: Perform grid analysis on the site design digital elevation model (DEM) and the status digital elevation model (DEM) to obtain a ground filling and digging grid data with the earthwork volume as the pixel value, and generate a site ground filling and digging thematic map;

[0012] Step 4: According to the control detailed planning data, make planning land surface vector data composed of road surfaces and land plots, and calculate the ground filling and digging earthwork volume of each plot in the planning land surface vector data based on the ground filling and digging grid data and the planning land surface vector data;

[0013] Step 5: According to the planning index data of the underground space of buildings, subways and underground comprehensive pipe galleries in the control detailed planning, calculate the underground engineering excavation volume of each plot in the planning land surface vector data; the underground engineering excavation volume includes the underground space excavation volume of buildings, the subway excavation volume and the underground comprehensive pipe gallery excavation volume; wherein the underground space excavation volume of buildings is calculated based on the plot area, the volume rate and the parking space index, and the subway and underground comprehensive pipe gallery excavation volumes are calculated based on the planning length and the cross section size;

[0014] Step 6: Based on the ground filling and digging volume and the underground engineering excavation volume of each plot, calculate the filling and digging volume of each plot, and generate a site filling and digging thematic map;

[0015] Step 7: Accumulate the excavation volume and the filling volume of all plots in the site respectively to obtain the total excavation volume and the total filling volume of the site, and continuously adjust the design elevations to re-enter step 2 for recalculation until the difference between the total excavation volume and the total filling volume of the site is less than a predetermined threshold.

[0016] The building underground space excavation amount calculation formula of each plot in step 5 is:

[0017] V u = S u × H u ;

[0018] S u = N t × S t ;

[0019] N t = S j / N j ;

[0020] S j = S d × R;

[0021] In the formula, V u is the building underground space excavation amount of the plot, S u is the building underground space land area of the plot, H u is the building underground space excavation thickness, N t is the number of underground parking spaces of the plot, S t is the land area per vehicle, S j is the above-ground building area of the plot, N j is the number of parking spaces per 100 square meters of building area, S d is the plot area, and R is the plot volume rate.

[0022] Among them, the building underground space excavation thickness is obtained by calculating the sum of the underground space beam height, net height, bottom plate thickness and pipeline reserved space height, and is calculated respectively according to the types of residential and commercial land.

[0023] The subway excavation amount calculation formula of each road surface in step 5 is:

[0024] V u = D c × S j ;

[0025] S j = π × (D z / 2) 2 ;

[0026] In the formula, V u is the subway excavation amount of the road surface, D c is the planned subway length of the road surface, S j is the cross-sectional area of the planned subway, and D z is the diameter of the planned subway.

[0027] The underground comprehensive pipe gallery excavation amount calculation formula of each road surface is:

[0028] V u = D c × S j ;

[0029] S j = W j × H j ;

[0030] In the formula, V u is the underground comprehensive pipe gallery excavation amount of the road surface, D c is the planned underground comprehensive pipe gallery length of the road surface, S j is the section size of the planned underground comprehensive pipe gallery, W j is the width of the planned underground comprehensive pipe gallery section, and H j is the height of the planned underground comprehensive pipe gallery section.

[0031] Preferably, in step 2, the control point generation step is:

[0032] Step 201: mark 1 control point on each road red line turning point at the intersection and the intersection center, and generate road intersection control points;

[0033] Step 202: mark 1 control point on the road center line and the road red line on both sides at the road section slope change point, and generate road section control points;

[0034] Step 203: mark control points at the edges and slope change points inside the plot, and generate road outside plot control points;

[0035] Step 204: take the existing terrain elevation points of the water system and reserved structures in the existing terrain map as control points in the non-construction area of the site.

[0036] Preferably, in step 2, the site design digital elevation model DEM and the site existing digital elevation model DEM have the same coordinate system and ground resolution.

[0037] Preferably, the step 3 comprises:

[0038] Step 301: calculate the difference between the pixel values in the site existing digital elevation model DEM and the design digital elevation model DEM by grid calculation, and obtain site terrain elevation difference grid data;

[0039] Step 302: set the pixel values in the non-construction area in the site terrain elevation difference grid data to 0 by grid calculation;

[0040] Step 303: multiply the pixel value in the terrain elevation grid data by the ground area represented by each pixel, to obtain site surface filling and digging grid data with pixel value as excavation earthwork volume;

[0041] Step 304: the positive value, negative value and 0 of the pixel in the site surface filling and digging grid data represent excavation, filling and unchanged terrain respectively, and a site surface filling and digging thematic map is made by setting different color symbols.

[0042] Preferably, the step 4 comprises:

[0043] Step 401: according to the control detailed planning data, a planning land surface vector data composed of road surfaces and land plots within the site is made, and the road surfaces are segmented at intersections;

[0044] Step 402: the site surface filling and digging grid data and the planning land surface vector data are subjected to spatial statistical analysis, and the site surface filling and digging earthwork volume of each plot in the planning land surface vector data is calculated.

[0045] Preferably, the step 6 comprises:

[0046] Step 601: the site surface filling and digging volume and the underground engineering excavation volume of each plot are summed up to calculate

[0047] the filling and digging volume of each plot;

[0048] Step 602: a site filling and digging thematic map is made by taking the digging and filling type and the land use type as characteristic values and by using the layer color setting mode, the corresponding digging and filling volume value is marked in the center of the plot, and the title, the compass, the scale, and the legend elements are added to the map.

[0049] The application discloses a site vertical design earthwork volume calculation method and system based on GIS, which is based on a site current terrain map, a site design map and control detailed planning data, calculates site surface terrain excavation and filling volume and underground engineering excavation volume, and obtains excavation volume and filling volume of the site.

[0050] 1. The application is based on GIS technology, can realize simple, rapid and high-precision calculation of earthwork volume in large-area site vertical design, and can rapidly calculate excavation and filling earthwork volume of any region in the site.

[0051] 2. The application takes the terrain map in site design as a data source, does not need to additionally collect laser point cloud, current status data such as oblique photography, saves current status data acquisition investment while ensuring high-precision calculation of earthwork volume, reduces the complexity of earthwork volume calculation, improves the efficiency of earthwork volume calculation, and expands the applicability.

[0052] 3. In addition to calculating the volume of excavated and filled earthwork on the surface, this invention also calculates the volume of excavated earthwork for underground spaces of buildings, subways, and underground integrated pipe corridors within the site. This makes the earthwork volume calculation more accurate, scientific, and reasonable, providing more reliable data support for the vertical design of the site. At the same time, combined with the development sequence of the plots, it details the earthwork allocation routes for each plot, effectively guiding the construction and development of the site, helping to reduce construction costs and accelerate the construction progress. Attached Figure Description

[0053] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are only embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on the provided drawings without creative effort.

[0054] Figure 1 This is a flowchart of a GIS-based method and system for calculating earthwork volume in site vertical design, according to an embodiment of the present invention. Detailed Implementation

[0055] To make the objectives, technical solutions, and advantages of this invention clearer, the invention will be further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative and not intended to limit the invention. All other embodiments obtained by those skilled in the art based on the embodiments of this invention without inventive effort are within the scope of protection of this invention.

[0056] According to an embodiment of the present invention, a method and system for calculating the earthwork volume of a site based on GIS is provided.

[0057] The present invention will now be further described in conjunction with the accompanying drawings and specific embodiments, such as... Figure 1 As shown, according to an embodiment of the present invention, a method for calculating earthwork volume in site vertical planning and design based on GIS is provided. The method includes the following steps:

[0058] Step 1: Obtain the existing topographic map file of the site, and construct a digital elevation model (DEM) of the existing site based on the elevation points in the topographic map, including:

[0059] Step 101: Collect 1:500 scale, CAD format topographic map files, perform data format conversion on the topographic map files by layer, convert the water system and building layers in the topographic map into area vector data, and convert the elevation annotations in the topographic map into elevation point vector data;

[0060] Step 102: After removing the abnormal points in the elevation point vector data, use the Kriging method to perform spatial interpolation processing on the elevation point vector data to generate an initial existing digital elevation model DEM raster file;

[0061] Step 103: Make a site range boundary vector data, and use the site range boundary to cut the initial existing digital elevation model DEM raster file to obtain a site existing digital elevation model DEM.

[0062] Step 2: According to the design elevation in the site vertical design drawing, a site design digital elevation model DEM is constructed, including:

[0063] Step 201: Collect CAD format site vertical design drawings with the same coordinate system as the existing topographic map file, and convert the design elevation annotation points in the site vertical design drawing into GIS format design elevation vector point data;

[0064] Step 201: Check and process the design elevation vector point data to ensure that the design elevation vector point data has no abnormal data and contains the following control points:

[0065] At each road red line turning point of the intersection and the intersection center, one control point is marked to generate a road intersection control point;

[0066] At the road red line and road center line at the slope change point of the road section, one control point is marked on each side to generate a road section control point;

[0067] Control points are marked at the edges and slope change points inside the plot to generate road outside plot control points;

[0068] The existing topographic elevation points in the existing topographic map in the non-construction area of the water system, reserved structures and the like are taken as control points in the non-construction area of the site.

[0069] Step 203: Use the Kriging method to perform spatial interpolation processing on the design elevation vector point data, set the same pixel size as the site existing digital elevation model DEM, and generate an initial design digital elevation model DEM raster file;

[0070] Step 204: Use the site range boundary to cut the initial design digital elevation model DEM raster file to obtain a site design digital elevation model DEM.

[0071] Step 3: Perform raster analysis on the site design digital elevation model DEM and the existing digital elevation model DEM to obtain a ground filling and excavation raster data with earthwork volume as the pixel value, and generate a site ground filling and excavation thematic map, including:

[0072] Step 301: Subtract the pixel values in the site status digital elevation model DEM and the design digital elevation model DEM by grid calculation to obtain site terrain elevation difference grid data;

[0073] Step 302: Set the pixel values in the site terrain elevation difference grid data in the non-construction area to 0 by grid calculation;

[0074] Step 303: Multiply the pixel values in the terrain elevation difference grid data by the ground area represented by each pixel to generate site surface filling and digging grid data with pixel values as excavation volume;

[0075] Step 304: The positive value, negative value and 0 of the pixel in the site surface filling and digging grid data represent excavation, filling and unchanged terrain respectively, and a site surface filling and digging thematic map is made by setting different color symbols.

[0076] Step 4: According to the control detailed planning data, a planning land surface vector data composed of road surfaces and land plots is made, and the surface filling and digging earthwork volume of each plot in the planning land surface vector data is calculated based on the surface filling and digging grid data and the planning land surface vector data, including:

[0077] Step 401: According to the control detailed planning data, a planning land surface vector data composed of road surfaces and land plots within the site range is made, and the road surfaces at intersections are segmented;

[0078] Step 402: Spatial statistical analysis is performed on the surface filling and digging grid data and the planning land surface vector data to calculate the surface filling and digging earthwork volume of each plot in the planning land surface vector data.

[0079] Step 5: According to the planning index data of building underground space, subway and underground comprehensive pipe gallery in the control detailed planning, the underground engineering excavation volume of each plot in the planning land surface vector data is calculated, including:

[0080] Step 501: According to the plot area, volume rate and allocated parking space data in the control detailed planning, the building underground space excavation volume of each plot is calculated, and the calculation formula is:

[0081] V u =S u ×H u ;

[0082] S u =N t ×S t ;

[0083] N t =S j / N j ;

[0084] S j =S d ×R;

[0085] In the formula, V u S represents the excavation volume for the underground space of the buildings on the plot. u H represents the area of ​​the underground space of the buildings on the plot. u N represents the thickness of the excavated soil for the underground space of a building. t S represents the number of underground parking spaces on the plot. t S is the area of ​​one underground parking space. j N represents the above-ground building area of ​​the plot. j S is the number of parking spaces to be provided per 100 square meters of building area. d R represents the plot area, and R represents the plot's floor area ratio.

[0086] The thickness of the excavated soil in the underground space of the building is obtained by calculating the sum of the beam height, net height, bottom slab thickness, and reserved space height for pipelines in the underground space, and is calculated separately for residential and commercial land types.

[0087] Step 502: Using road surfaces as the calculation unit, calculate the excavation volume for the subway and the underground utility tunnel for each road surface based on the planned length and cross-sectional dimensions of the subway and underground utility tunnel. The formula for calculating the subway excavation volume for each road surface is as follows:

[0088] V u =D c ×Sj;

[0089] S j =π×(D) z / 2) 2 ;

[0090] In the formula, V u D represents the volume of subway excavation for the road surface. c S represents the planned length of the subway line across the road surface. j To plan the cross-sectional area of ​​the subway, D z To determine the diameter of the planned subway line;

[0091] The formula for calculating the excavation volume of underground utility tunnels on each road surface is as follows:

[0092] V u =D c ×S j ;

[0093] S j =W j ×H j ;

[0094] In the formula, Vu D is the underground comprehensive pipe gallery excavation quantity of the road surface c S is the planned underground comprehensive pipe gallery length of the road surface j W is the section size of the planned underground comprehensive pipe gallery j H is the width of the section of the planned underground comprehensive pipe gallery j H is the height of the section of the planned underground comprehensive pipe gallery.

[0095] Step 503: summing up the subway excavation quantity and the underground comprehensive pipe gallery excavation quantity of each road surface to obtain the underground engineering excavation quantity of each plot in the planned land surface vector data.

[0096] Step 6: calculating the fill and excavation quantity of each plot based on the ground surface fill and excavation quantity and the underground engineering excavation quantity of each plot, and generating a site fill and excavation thematic map, including:

[0097] Step 601: summing up the ground surface fill and excavation quantity and the underground engineering excavation quantity of each plot by attribute calculation to calculate the fill and excavation quantity of each plot;

[0098] Step 602: taking the excavation and filling type as a characteristic value, and making a site fill and excavation thematic map by layer color setting, marking the corresponding fill and excavation quantity value in the center of the plot, and adding title, compass, scale, and legend elements.

[0099] Step 603: taking the land use type as a characteristic value, and making a site fill and excavation thematic map by layer color setting, marking the corresponding fill and excavation quantity value in the center of the plot, and adding title, compass, scale, and legend elements.

[0100] Step 7: respectively accumulating the excavation quantity and the fill quantity of all plots in the site to obtain the total excavation quantity and the total fill quantity of the site, and continuously adjusting the design elevation to re-enter step 2 for recalculation until the difference between the total excavation quantity and the total fill quantity of the site is not more than a predetermined threshold.

[0101] According to another embodiment of the application, a GIS-based site vertical design earthwork quantity calculation system is provided, comprising:

[0102] A current DEM construction module is configured to perform spatial interpolation and clipping processing on the elevation points in the site current topographic map to generate a site current digital elevation model DEM.

[0103] A design DEM construction module is configured to perform spatial interpolation and clipping processing on the design elevation data in the site design map to generate a site design digital elevation model DEM.

[0104] The ground excavation and filling amount calculation module is configured to perform grid analysis on the site status digital elevation model DEM and the site design digital elevation model DEM, obtain ground excavation and filling grid data of the site by combining a ground area represented by each pixel, and calculate ground excavation and filling earthwork amount of each plot in the planning land surface vector data through spatial statistical analysis.

[0105] The underground engineering excavation amount calculation module is configured to calculate underground engineering excavation amount of each plot in the planning land surface vector data according to planning index data of underground space of buildings, subways and underground comprehensive pipe galleries in the control detailed planning.

[0106] The total excavation and filling amount calculation module is configured to sum up the ground excavation and filling amount and the underground engineering excavation amount of each plot in the planning land surface vector data, obtain excavation and filling amount of each plot, generate a site excavation and filling thematic map, and accumulate excavation amount and filling amount of all plots in the site to obtain total excavation amount and total filling amount of the site.

Claims

1. A GIS-based site vertical design earthwork volume calculation method, characterized in that, The method comprises the following steps: Step 1: Obtain a site current topographic map file, and construct a site current digital elevation model (DEM) according to elevation points in the topographic map; Step 2: Construct a site design digital elevation model (DEM) according to design elevations in a site vertical design drawing, wherein the design elevations are derived from control points marked in the site design drawing, and the control points include road intersection control points, road section control points, road outer plot control points, and control points at non-construction areas of the site; Step 3: Perform grid analysis on the site design digital elevation model (DEM) and the current digital elevation model (DEM) to obtain a ground surface filling and digging grid data with earthwork volume as pixel value, and generate a site ground surface filling and digging thematic map; Step 4: Produce planning land surface vector data composed of road surfaces and plots according to control detailed planning data, and calculate earthwork volume of each plot in the planning land surface vector data based on the ground surface filling and digging grid data and the planning land surface vector data; Step 5: Calculate underground engineering excavation volume of each plot in the planning land surface vector data according to planning index data of underground space of buildings, subways, and underground comprehensive pipe galleries in the control detailed planning, wherein the underground engineering excavation volume includes building underground space excavation volume, subway excavation volume, and underground comprehensive pipe gallery excavation volume; The building underground space excavation volume is calculated based on plot area, volume rate, and parking space index, and the subway and underground comprehensive pipe gallery excavation volume is calculated based on planning length and cross section size; Step 6: Calculate filling and digging volume of each plot based on the ground surface filling and digging volume and the underground engineering excavation volume of each plot, and generate a site filling and digging volume thematic map; Step 7: Accumulate excavation volume and filling volume of all plots in the site respectively to obtain total excavation volume and total filling volume of the site, and continuously adjust the design elevations to re-enter step 2 for recalculation until the difference between the total excavation volume and the total filling volume of the site is less than a predetermined threshold; The building underground space excavation volume calculation formula of each plot in step 5 is as follows: V u = Su x H u ; S u = N t x S t ; N t = S j / N j ; S j = S d x R; In the formula, V u is the amount of excavation of the underground space of the building of the plot, S u is the land area of the underground space of the building of the plot, H u is the thickness of the excavated soil of the underground space of the building, N t is the number of underground parking spaces of the plot, S t is the land area per car, S j is the above-ground building area of the plot, N j is the number of parking spaces per 100 square meters of building area, S d is the plot area, R is the plot ratio; The building underground space excavation thickness is obtained by calculating the sum of beam height, net height, bottom plate thickness, and pipe reserved space height, and is calculated respectively for residential and commercial land types; The subway excavation volume calculation formula of each road surface in step 5 is as follows: V u = D c x S j ; S j = π x (D z / 2) 2 ; wherein V u is the amount of subway excavation of the road surface, D c is the planned subway length of the road surface, S j is the cross-sectional area of the planned subway, D z is the diameter of the planned subway; The underground comprehensive pipe gallery excavation volume calculation formula of each road surface in step 5 is as follows: V u = D c × S j ; S j = W j x H j ; In the formula, V u is the underground comprehensive pipe gallery excavation amount of the road surface, D c is the planned underground comprehensive pipe gallery length of the road surface, S j is the cross-sectional dimension of the planned underground comprehensive pipe gallery, W j is the width of the cross section of the planned underground comprehensive pipe gallery, H j is the height of the cross section of the planned underground comprehensive pipe gallery.

2. The GIS-based site vertical design earthwork volume calculation method of claim 1, wherein, In step 2, the control point generation steps are as follows: Step 201: Mark one control point at each road red line turning point of the intersection and the intersection center to generate road intersection control points; Step 202: Mark one control point on the road center line and the road red line on both sides of the road section slope change point to generate road section control points; Step 203: Mark control points at the edges and slope change points inside the plot to generate road outer plot control points; Step 204: Take the water system and the current topographic elevation points at non-construction areas of reserved structures in the current topographic map as control points at non-construction areas of the site. 3.The GIS-based site vertical design earthwork volume calculation method of claim 1, wherein, In step 2, the site design digital elevation model DEM and the site status digital elevation model DEM have the same coordinate system and ground resolution. 4.The GIS-based site vertical design earthwork volume calculation method of claim 1, wherein, The step 3 comprises: Step 301: calculating the difference between the site status digital elevation model DEM and the design digital elevation model DEM by grid calculation, and obtaining site terrain elevation difference grid data; Step 302: setting the pixel value in the site terrain elevation difference grid data in the non-construction area to 0 by grid calculation; Step 303: multiplying the pixel value in the terrain elevation difference grid data by the ground area represented by each pixel, and obtaining site surface filling and digging grid data with pixel value as the earthwork volume; Step 304: the positive value, negative value and 0 of the site surface filling and digging grid data represent excavation, filling and unchanged terrain respectively, and a site surface filling and digging thematic map is made by setting different color symbols.

5. The GIS-based site vertical design earthwork volume calculation method of claim 1, wherein, The step 4 comprises: Step 401: making planning land surface vector data composed of road surfaces and land plots within the site range according to the control detailed planning data, and dividing the road surfaces at intersections; Step 402: performing spatial statistical analysis on the site surface filling and digging grid data and the planning land surface vector data, and calculating the surface filling and digging earthwork volume of each plot in the planning land surface vector data.

6. The GIS-based site vertical design earthwork volume calculation method of claim 1, wherein, The step 6 comprises: Step 601: summing the surface filling and digging volume and the underground engineering excavation volume of each plot to calculate the filling and digging volume of each plot; Step 602: making a site filling and digging thematic map by layer color setting with the filling and digging type and land use type as characteristic values, marking the corresponding filling and digging volume values in the plot center, and adding title, compass, scale and legend elements in the map.

Images