Method and device for constructing three-dimensional evaluation model of urban underground space resources

Abstract

The present invention discloses a method and device for constructing a three-dimensional evaluation model of urban underground space resources. The method comprises: (1) generating a two-dimensional plane unit set A and a stratum set D according to a two-dimensional vector plane unit layer and a three-dimensional stratum model of the urban underground space resources respectively; (2) generating a three-dimensional space unit set U based on the preset depth of each layer and the set A; (3) assigning stratum attributes to each three-dimensional space unit in U based on the set U and the set D to obtain a stratum attribute set T; (4) assigning index attributes to each three-dimensional space unit in U based on the index attribute value of A to obtain an index attribute set Z; (5) merging the sets T and Z to obtain an attribute set S, and calculating the appropriate development degree value of each three-dimensional space unit based on S; (6) classifying and merging U based on the set S to obtain a three-dimensional evaluation model of urban underground space resources. The present invention has high processing efficiency and a wide range of applications.

Description

Technical Field

[0001] The present invention relates to geographic information technology, and in particular to a method and device for constructing a three-dimensional evaluation model of urban underground space resources. Background Art

[0002] As an important reserve resource for the expansion of national space, the scientific development and sustainable utilization of underground space is a major national demand and an important indicator of the process of urban modernization. For the construction of a three-dimensional evaluation model for urban underground space resources, current research is mainly limited to two-dimensional plane layered evaluation methods with different depth layers or three-dimensional evaluation methods based on regular discrete grids. The former easily separates the three-dimensional characteristics of underground space resources at the vertical level, and the processing process is cumbersome and inefficient; the latter is only applicable to key geological survey areas with abundant geological survey data and small scope due to its high data requirements. Summary of the invention

[0003] Purpose of the invention: In view of the problems existing in the prior art, the present invention provides a method and device for constructing a three-dimensional evaluation model of urban underground space resources with high efficiency and wide applicability.

[0004] Technical solution: The method for constructing a three-dimensional evaluation model of urban underground space resources described in the present invention includes:

[0005] (1) Generate a two-dimensional plane unit set A and a stratum set D based on the two-dimensional vector plane unit layer and the three-dimensional stratum model of the urban underground space resources respectively;

[0006] (2) Based on the preset depths of each layer and the two-dimensional plane unit set A, a three-dimensional space unit set U is generated;

[0007] (3) Based on the three-dimensional space unit set U and the stratum set D, a stratum attribute is assigned to each three-dimensional space unit in the three-dimensional space unit set U to obtain a stratum attribute set T;

[0008] (4) According to the index attribute value of the two-dimensional plane unit set A, the index attribute value is assigned to each three-dimensional space unit in the three-dimensional space unit set U to obtain the index attribute set Z;

[0009] (5) The formation attribute set T and the indicator attribute set Z are combined to obtain the attribute set S. Based on the attribute set S, the appropriate development degree value of each three-dimensional space unit is calculated and stored in the attribute set S;

[0010] (6) Based on the attribute set S, the three-dimensional space unit set U is classified and merged to obtain a three-dimensional evaluation model of urban underground space resources.

[0011] Furthermore, step (1) comprises:

[0012] (1-1) Load the two-dimensional vector plane unit layer and store it in the two-dimensional plane unit set A = {a i |i=1,2,…,AN}, where a i represents the i-th two-dimensional plane unit, AN represents the number of two-dimensional plane units;

[0013] (1-2) Load the 3D stratigraphic model data and store it in stratigraphic set D = {d k |k=1,2,…,DN}, where d k represents the kth stratum after sorting from shallow to deep, and DN represents the number of strata.

[0014] Furthermore, step (2) includes:

[0015] (2-1) Store the preset depths of each layer into the depth set DEP = {dep j |j=1,2,…,n}, where dep j represents the depth of the jth depth layer from shallow to deep, and n represents the number of layers;

[0016] (2-2) Select any element a in the two-dimensional plane unit set A i , and obtain a i Point set P2D i ;

[0017] (2-3) Select any depth dep in the depth set DEP j ;

[0018] (2-4) a i Point set P2D i Give depth dep j , and get the three-dimensional point set P3D ij ;

[0019] (2-5) The three-dimensional point set P3D ij Convert to 3D surface feature u ij , and stored as a three-dimensional space unit in the three-dimensional space unit set U;

[0020] (2-6) Loop through steps (2-3) to (2-5) until the depth set DEP is traversed and a is obtained. i The three-dimensional space unit corresponding to each depth;

[0021] (2-7) Loop through steps (2-2) to (2-6) until the two-dimensional plane unit set A is traversed and the three-dimensional space unit set U = {u ij |i=1,2,…,AN,j=1,2,…,n}, AN represents the number of elements in the two-dimensional plane unit set A.

[0022] Furthermore, step (3) includes:

[0023] (3-1) Select any unit u from the three-dimensional space unit set U ij ;

[0024] (3-2) Select any stratum d in the stratum set D k ;

[0025] (3-3) Based on u ij and d k , the formation d is calculated by the intersection Boolean operation method of three-dimensional polyhedrons. k in u ij The volume v of the intersecting part ij,k ,;

[0026] (3-4) Loop through steps (3-2)-(3-3) until all stratum sets are traversed and u is obtained. ij The set of intersection volumes V with all layers ij = {v ij,k |k=1,2,…,DN}, DN represents the number of strata;

[0027] (3-5) From the set V ij Select the stratum d corresponding to the largest intersecting volume;

[0028] (3-6) Assign the attributes of layer d to u ij The stratigraphic properties of ij ;

[0029] (3-7) Loop through steps (3-1) to (3-6) until the three-dimensional space unit set U is completely traversed, and obtain the stratigraphic attributes of all three-dimensional space units, which are stored in the stratigraphic attribute set T = {t ij |i=1,2,…,AN,j=1,2,…,n}.

[0030] Further, step (4) includes:

[0031] (4-1) Get any element a from the two-dimensional plane unit set A i ;

[0032] (4-2) Select the cells that match a from the three-dimensional space unit set U i The projections on the xy plane are exactly the same, but the three-dimensional space unit subsets U belong to different depth layers i = {u ij |j=1,2,…,n};

[0033] (4-3) From the three-dimensional space unit subset U i Read any three-dimensional space unit uij ;

[0034] (4-4) According to a i The sequence number i and the depth layer j to which it belongs generate a three-dimensional space unit u ij The corresponding unit number ij, and a i The index attribute is assigned to u ij The indicator attribute z ij ;

[0035] (4-5) Loop through steps (4-3) to (4-4) until U i After being traversed, we get U i The index properties of all three-dimensional space units in;

[0036] (4-6) Loop through steps (4-1) to (4-5) until the two-dimensional plane unit set A is traversed, and the index attributes of all three-dimensional space units are obtained and stored in the index attribute set Z = {z ij |i=1,2,…,AN,j=1,2,…,n};

[0037] Further, step (5) includes:

[0038] (5-1) Merge the formation attribute set T and the indicator attribute set Z to obtain the attribute set S = {s ij ={t ij ,z ij}|i=1,2,…,AN,j=1,2,…,n},s ij Represents the three-dimensional space unit u ij The indicator attributes of

[0039] (5-2) Select any u from the attribute set S ij The indicator attributes of ij ;

[0040] (5-2) According to the attribute s ij Various index data and formation information in the calculation unit u ij The appropriate development degree value is stored in the attribute s ij middle;

[0041] (5-3) Repeat steps (5-2)-(5-3) until the calculation of the appropriate development degree values ​​of all three-dimensional space units is completed.

[0042] Further, step (6) comprises:

[0043] (6-1) Read the preset thresholds of appropriate development levels and store them in the threshold set E = {e m |m=1,2,…,EN}, e mrepresents the threshold value of the mth level of appropriate development, and EN represents the level number;

[0044] (6-2) Select any threshold e from the threshold set E m ;

[0045] (6-3) From the three-dimensional space unit set U and the attribute set S, select the suitable development degree value in the interval (e m-1 ,e m ] and the corresponding attributes, and store them in the m-th level three-dimensional space unit set U m = {u ml |l=1,2,…,MN m} and attribute set S m ={s ml |l=1,2,…,MN m}, where u ml Represents the lth three-dimensional space unit in the mth level, MN m Represents the number of three-dimensional space units at the mth level;

[0046] (6-4) The three-dimensional space unit set U m All the elements in are combined into a three-dimensional solid element u m , the attribute set S m All attributes in are merged into one attribute s m ;

[0047] (6-5) Loop through steps (6-2) to (6-4) until the three-dimensional space unit set U is completely traversed, and all three-dimensional space units are classified to obtain the classified three-dimensional unit set U' = {u m |m=1,2,…,EN} and the corresponding attribute set S'={s m |m=1,2,…,EN};

[0048] (6-6) The classified three-dimensional unit set U' is converted into an obj format model file, and the attribute set S' is converted into a json format file to complete the three-dimensional evaluation model of underground space resources.

[0049] The device for constructing a three-dimensional evaluation model of urban underground space resources described in the present invention includes a processor and a computer program stored in a memory and executable on the processor, and the processor implements the above method when executing the program.

[0050] Beneficial effects: Compared with the prior art, the present invention has the following significant advantages: based on the two-dimensional index data affecting the development and utilization of underground space resources and the three-dimensional stratigraphic model in the region, the present invention generates three-dimensional space units with different depth layers, and then evaluates each three-dimensional space unit through the urban and rural land use evaluation standards, thereby obtaining a three-dimensional evaluation model for urban underground space resources. The present invention has high efficiency, simple processing, and low data requirements. Therefore, it has a wide range of applications. The present invention is the basis and premise for the development and utilization of urban underground space, and has important research significance and application value. BRIEF DESCRIPTION OF THE DRAWINGS

[0051] Figure 1 is the two-dimensional plane unit data used in this embodiment;

[0052] Figure 2 is the three-dimensional stratum model data used in this embodiment;

[0053] Figure 3 is a flow chart of the present invention;

[0054] Figure 4 This is the display effect of the three-dimensional space unit in this embodiment;

[0055] Figure 5 This is the final display effect after the three-dimensional space units are classified and merged in this embodiment. DETAILED DESCRIPTION

[0056] like Figure 1 and Figure 2 As shown, this embodiment selects two-dimensional evaluation unit vector data of a certain test area and three-dimensional model data of loose layer as experimental data, and the projection coordinate system used by the experimental data is the local coordinate system of the test area. The following is further explained by combining the accompanying drawings and describing a specific embodiment.

[0057] like Figure 3 As shown, this embodiment provides a method for constructing a three-dimensional evaluation model of urban underground space resources, including:

[0058] (1) Generate a two-dimensional plane unit set A and a stratigraphic set D based on the two-dimensional vector plane unit layer and the three-dimensional stratigraphic model of the urban underground space resources.

[0059] This step specifically includes:

[0060] (1-1) Load the two-dimensional vector plane unit layer and store the two-dimensional plane unit set A = {a i |i=1,2,…,AN}, where a i represents the i-th two-dimensional plane unit, and AN represents the number of two-dimensional plane units;

[0061] (1-2) Load the 3D stratigraphic model data and store it in stratigraphic set D = {d k |k=1,2,…,DN}, where d k represents the kth stratum after sorting from shallow to deep, and DN represents the number of strata.

[0062] In this embodiment, AN is 25, DN is 7, and the attributes of each unit are shown in Table 1;

[0063] Table 1 Two-dimensional plane element attribute table

[0064]

[0065]

[0066] (2) Based on the preset depths of each layer and the two-dimensional plane unit set A, a three-dimensional space unit set U is generated.

[0067] This step specifically includes:

[0068] (2-1) Store the preset depths of each layer into the depth set DEP = {dep j |j=1,2,…,n}, where dep j represents the depth of the jth depth layer from shallow to deep, and n represents the number of layers; in this embodiment, referring to the "Urban and Rural Land Assessment Standard CJJ132-2009", the depth set is set as DEP = {10m, 30m, 50m, 100m}, and n is 4;

[0069] (2-2) Select any element a in the two-dimensional plane unit set A i , and obtain a i Point set P2D i ;

[0070] (2-3) Select any depth dep in the depth set DEP j ;

[0071] (2-4) a i Point set P2D i Give depth dep j , and get the three-dimensional point set P3D ij ;

[0072] (2-5) The three-dimensional point set P3D ij Convert to 3D surface feature u ij , and stored as a three-dimensional space unit in the three-dimensional space unit set U;

[0073] (2-6) Loop through steps (2-3) to (2-5) until the depth set DEP is traversed and a is obtained. iThe three-dimensional space unit corresponding to each depth;

[0074] (2-7) Loop through steps (2-2) to (2-6) until the two-dimensional plane unit set A is traversed and the three-dimensional space unit set U = {u ij |i=1,2,…,AN,j=1,2,…,n}, such as Figure 4 shown.

[0075] (3) Based on the three-dimensional space unit set U and the stratum set D, the stratum attribute is assigned to each three-dimensional space unit in the three-dimensional space unit set U to obtain the stratum attribute set T.

[0076] This step specifically includes:

[0077] (3-1) Select any unit u from the three-dimensional space unit set U ij ;

[0078] (3-2) Select any stratum d in the stratum set D k ;

[0079] (3-3) Based on u ij and d k , the formation d is calculated by the intersection Boolean operation method of three-dimensional polyhedrons. k in u ij The volume v of the intersecting part ij,k ,;

[0080] (3-4) Loop through steps (3-2)-(3-3) until all stratum sets are traversed and u is obtained. ij The set of intersecting volumes V with all layers ij = {v ij,k |k=1,2,…,DN}, DN represents the number of strata;

[0081] (3-5) From the set V ij Select the stratum d corresponding to the largest intersecting volume;

[0082] (3-6) Assign the attributes of layer d to u ij The stratigraphic properties of ij ;

[0083] (3-7) Loop through steps (3-1) to (3-6) until the three-dimensional space unit set U is completely traversed, and obtain the stratigraphic attributes of all three-dimensional space units, which are stored in the stratigraphic attribute set T = {t ij |i=1,2,…,AN,j=1,2,…,n}.

[0084] (4) According to the index attribute value of the two-dimensional plane unit set A, the index attribute is assigned to each three-dimensional space unit in the three-dimensional space unit set U to obtain the index attribute set Z.

[0085] This step specifically includes:

[0086] (4-1) Get any element a from the two-dimensional plane unit set A i ;

[0087] (4-2) Select the cells that match a from the three-dimensional space unit set U i The projections on the xy plane are exactly the same, but the three-dimensional space unit subsets U belong to different depth layers i = {u ij |j=1,2,…,n};

[0088] (4-3) From the three-dimensional space unit subset U i Read any three-dimensional space unit u ij ;

[0089] (4-4) According to a i The sequence number i and the depth layer j to which it belongs generate a three-dimensional space unit u ij The corresponding unit number ij, and a i The index attribute is assigned to u ij The indicator attribute z ij ;

[0090] (4-5) Loop through steps (4-3) to (4-4) until U i After being traversed, we get U i The index properties of all three-dimensional space units in;

[0091] (4-6) Loop through steps (4-1) to (4-5) until the two-dimensional plane unit set A is traversed, and the index attributes of all three-dimensional space units are obtained and stored in the index attribute set Z = {z ij |i=1,2,…,AN,j=1,2,…,n};

[0092] (5) The stratum attribute set T and the indicator attribute set Z are combined to obtain the attribute set S. Based on the attribute set S, the appropriate development degree value of each three-dimensional space unit is calculated and stored in the attribute set S.

[0093] (5-1) Merge the formation attribute set T and the indicator attribute set Z to obtain the attribute set S = {s ij ={t ij ,z ij}|i=1,2,…,AN,j=1,2,…,n},s ij Represents the three-dimensional space unit uij The indicator attributes of

[0094] (5-2) Select any u from the attribute set S ij The indicator attributes of ij ;

[0095] (5-2) According to the attribute s ij Based on the various index data and stratum information in the land use assessment model, unit u is calculated ij The appropriate development degree value is stored in the attribute s ij The land use assessment model is an existing model and will not be introduced here.

[0096] (5-3) Repeat steps (5-2)-(5-3) until the calculation of the appropriate development degree values ​​of all three-dimensional space units is completed.

[0097] (6) Based on the attribute set S, the three-dimensional space unit set U is classified and merged to obtain a three-dimensional evaluation model of urban underground space resources, such as Figure 5 shown.

[0098] This step specifically includes:

[0099] (6-1) Read the preset thresholds of appropriate development levels and store them in the threshold set E = {e m |m=1,2,…,EN}, e m represents the threshold value of the mth level of suitable development, and EN represents the level number; in this embodiment, referring to the "Urban and Rural Land Assessment Standard CJJ132-2009", the threshold set is set as E = {1, 3, 6, ∞}, e0 = 0, and the corresponding suitable development types are non-construction land, unsuitable construction land, constructible land and suitable construction land;

[0100] (6-2) Select any threshold e from the threshold set E m ;

[0101] (6-3) From the three-dimensional space unit set U and the attribute set S, select the suitable development degree value in the interval (e m-1 ,e m ] and the corresponding attributes, and store them in the m-th level three-dimensional space unit set U m = {u ml |l=1,2,…,MN m} and attribute set S m ={s ml |l=1,2,…,MN m}, where u ml Represents the lth three-dimensional space unit in the mth level, MN mRepresents the number of three-dimensional space units at the mth level;

[0102] (6-4) The three-dimensional space unit set U m All the elements in are combined into a three-dimensional solid element u m , the attribute set S m All attributes in are merged into one attribute s m ;

[0103] (6-5) Loop through steps (6-2) to (6-4) until the three-dimensional space unit set U is completely traversed, and all three-dimensional space units are classified to obtain the classified three-dimensional unit set U' = {u m |m=1,2,…,EN} and the corresponding attribute set S'={s m |m=1,2,…,EN};

[0104] (6-6) The classified three-dimensional unit set U' is converted into an obj format model file, and the attribute set S' is converted into a json format file to complete the three-dimensional evaluation model of underground space resources.

[0105] In this embodiment, shp data is read only based on the feature processing interface provided by the GDAL open source code, and other libraries such as ArcEngine can also be used to read shp data. Model intersection and merging processing is performed only based on the three-dimensional model Boolean operation interface provided by the CGAL open source code, and three-dimensional Boolean operations can also be performed using interfaces provided by other libraries such as meshlab. The three-dimensional space unit model is only exported in OBJ format, and the three-dimensional stratum model in other formats such as FBX can also be exported. In this embodiment, the attribute data of the three-dimensional space unit is only exported in json format, and it can also be exported in other formats or integrated into the OBJ file for export together.

[0106] This embodiment also provides a device for constructing a three-dimensional evaluation model of urban underground space resources, including a processor and a computer program stored in a memory and executable on the processor, and the processor implements the above method when executing the program.

[0107] The above disclosure is only a preferred embodiment of the present invention, which cannot be used to limit the scope of the present invention. Therefore, equivalent changes made according to the claims of the present invention are still within the scope of the present invention.

Claims

1. A method for constructing a three-dimensional evaluation model for urban underground space resources, characterized in that The method includes: (1) Generate a two-dimensional plane unit set A and a stratum set D based on the two-dimensional vector plane unit layer and the three-dimensional stratum model of the urban underground space resources respectively; (2) Based on the preset depths of each layer and the two-dimensional plane unit set A, a three-dimensional space unit set U is generated; (3) Based on the three-dimensional space unit set U and the stratum set D, a stratum attribute is assigned to each three-dimensional space unit in the three-dimensional space unit set U to obtain a stratum attribute set T; (4) According to the index attribute value of the two-dimensional plane unit set A, the index attribute value is assigned to each three-dimensional space unit in the three-dimensional space unit set U to obtain the index attribute set Z; (5) The formation attribute set T and the indicator attribute set Z are combined to obtain the attribute set S. Based on the attribute set S, the appropriate development degree value of each three-dimensional space unit is calculated and stored in the attribute set S; (6) Based on the attribute set S, the three-dimensional space unit set U is classified and merged to obtain a three-dimensional evaluation model of urban underground space resources; Step (1) comprises: (1-1) Load the two-dimensional vector plane unit layer and store the two-dimensional plane unit set A = {a i |i=1,2,…,AN}, where a i represents the i-th two-dimensional plane unit, and AN represents the number of two-dimensional plane units; (1-2) Load the 3D stratigraphic model data and store it in stratigraphic set D = {d k |k=1,2,…,DN}, where d k represents the kth stratigraphic layer after sorting from shallow to deep, and DN represents the number of stratigraphic layers; Step (2) comprises: (2-1) Store the preset depths of each layer into the depth set DEP = {dep j |j=1,2,…,n}, where dep j represents the depth of the jth depth layer from shallow to deep, and n represents the number of layers; (2-2) Select any element a in the two-dimensional plane unit set A i , and obtain a i Point set P2D i ; (2-3) Select any depth dep in the depth set DEP j ; (2-4) a i Point set P2D i Give depth dep j , and get the three-dimensional point set P3D ij ; (2-5) The three-dimensional point set P3D ij Convert to 3D surface feature u ij , and stored as a three-dimensional space unit in the three-dimensional space unit set U; (2-6) Loop through steps (2-3) to (2-5) until the depth set DEP is traversed and a is obtained. i The three-dimensional space unit corresponding to each depth; (2-7) Loop through steps (2-2) to (2-6) until the two-dimensional plane unit set A is traversed and the three-dimensional space unit set U = {u ij |i=1,2,…,AN,j=1,2,…,n}, AN represents the number of elements in the two-dimensional plane unit set A; Step (3) includes: (3-1) Select any unit u from the three-dimensional space unit set U ij ; (3-2) Select any stratum d in the stratum set D k ; (3-3) Based on u ij and d k , the formation d is calculated by the intersection Boolean operation method of three-dimensional polyhedrons. k in u ij The volume v of the intersecting part ij,k ; (3-4) Loop through steps (3-2)-(3-3) until all stratum sets are traversed and u is obtained. ij The set of intersection volumes V with all layers ij = {v ij,k |k=1,2,…,DN}, DN represents the number of strata; (3-5) From the set V ij Select the stratum d corresponding to the largest intersecting volume; (3-6) Assign the attributes of layer d to u ij The stratigraphic properties of ij ; (3-7) Loop through steps (3-1) to (3-6) until the three-dimensional space unit set U is completely traversed, and obtain the stratigraphic attributes of all three-dimensional space units, which are stored in the stratigraphic attribute set T = {t ij |i=1,2,…,AN,j=1,2,…,n}; Step (4) comprises: (4-1) Get any element a from the two-dimensional plane unit set A i ; (4-2) Select the cells that match a from the three-dimensional space unit set U i The projections on the xy plane are exactly the same, but the three-dimensional space unit subsets U belong to different depth layers i = {u ij |j=1,2,…,n}; (4-3) From the three-dimensional space unit subset U i Read any three-dimensional space unit u ij ; (4-4) According to a i The sequence number i and the depth layer j to which it belongs generate a three-dimensional space unit u ij The corresponding unit number ij, and a i The index attribute is assigned to u ij The indicator attribute z ij ; (4-5) Loop through steps (4-3) to (4-4) until U i After traversing, we get U i The index properties of all three-dimensional space units in; (4-6) Loop through steps (4-1) to (4-5) until the two-dimensional plane unit set A is traversed, and the index attributes of all three-dimensional space units are obtained and stored in the index attribute set Z = {z ij |i=1,2,…,AN,j=1,2,…,n}; Step (5) comprises: (5-1) Merge the formation attribute set T and the indicator attribute set Z to obtain the attribute set S = {s ij ={t ij ,z ij }|i=1,2,…,AN,j=1,2,…,n},s ij Represents the three-dimensional space unit u ij The indicator attributes of (5-2) Select any u from the attribute set S ij The indicator attributes of ij ; (5-2) According to the attribute s ij Various index data and formation information in the calculation unit u ij The appropriate development degree value is stored in the attribute s ij middle; (5-3) Execute steps (5-2) to (5-3) repeatedly until the calculation of the appropriate development degree values ​​of all three-dimensional space units is completed; Step (6) comprises: (6-1) Read the preset thresholds of appropriate development levels and store them in the threshold set E = {e m |m=1,2,…,EN}, e m represents the threshold value of the mth level of appropriate development, and EN represents the level number; (6-2) Select any threshold e from the threshold set E m ; (6-3) From the three-dimensional space unit set U and the attribute set S, select the suitable development degree value in the interval (e m-1 ,e m ] and the corresponding attributes, and store them in the m-th level three-dimensional space unit set U m = {u ml |l=1,2,…,MN m } and attribute set S m ={s ml |l=1,2,…,MN m }, where u ml Represents the lth three-dimensional space unit in the mth level, MN m Represents the number of three-dimensional space units at the mth level; (6-4) The three-dimensional space unit set U m All the elements in are combined into a three-dimensional solid element u m , the attribute set S m All attributes in s are combined into one attribute s m ; (6-5) Loop through steps (6-2) to (6-4) until the three-dimensional space unit set U is completely traversed, and all three-dimensional space units are classified to obtain the classified three-dimensional unit set U' = {u m |m=1,2,…,EN} and the corresponding attribute set S'={s m |m=1,2,…,EN}; (6-6) The classified three-dimensional unit set U' is converted into an obj format model file, and the attribute set S' is converted into a json format file to complete the three-dimensional evaluation model of underground space resources.

2. A device for constructing a three-dimensional evaluation model of urban underground space resources, comprising a processor and a computer program stored in a memory and executable on the processor, characterized in that: When the processor executes the program, the method according to claim 1 is implemented.

Images