A real scene three-dimensional archive data application method and device and readable storage medium

Abstract

The application discloses a kind of real scene three-dimensional archive data claim method, device and readable storage medium, comprising: obtaining real scene three-dimensional data and the range of claim space of data claim;Based on Tile operation rule, real scene three-dimensional data is framed and handled;Query the all Tile map corresponding to the range of claim space in Tile map data, and obtain the space range of each Tile map;If the range of claim space contains the space range of current Tile map, then the Tile map data corresponding to current Tile map is converted into real scene three-dimensional data and output;If the range of claim space intersects with the space range of current Tile map, obtain the intersection area of the range of claim space and the space range of current Tile map, according to intersection area, current Tile map is cut, and the Tile map data of cutting region is reserved, the Tile map data in cutting region is converted into real scene three-dimensional data and output.The application realizes the adaptive response to different use roles, real-time demand, desensitization level.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of image processing, and in particular to a real scene three-dimensional archive data application method, device and readable storage medium. BACKGROUND

[0002] In the application scenario of archive data application, the rapid application and sharing of real scene three-dimensional data are hindered by the complexity of the data, which reduces the practical value. At present, the cutting technology in the field of real scene three-dimensional archives is generally based on static hierarchical pre-generation, and lacks adaptive response capability for different user roles, real-time needs and desensitization levels.

[0003] Therefore, the present application develops a real scene three-dimensional archive data application method, device and readable storage medium to solve the above problems. SUMMARY

[0004] The present application proposes a real scene three-dimensional archive data application method, device and readable storage medium to solve the problem of lack of adaptive response capability for different user roles, real-time needs and desensitization levels in the existing real scene three-dimensional archive data application method.

[0005] The present application achieves the above-mentioned purposes through the following technical solutions:

[0006] The present application provides a real scene three-dimensional archive data application method, which comprises:

[0007] Obtaining real scene three-dimensional data and an application space range of data application;

[0008] Performing framing processing on the real scene three-dimensional data based on a Tile operation rule to obtain Tile map data;

[0009] Querying all Tile maps corresponding to the application space range in the Tile map data, and obtaining the space range of each Tile map;

[0010] Judging the positional relationship between the space range of each Tile map and the application space range;

[0011] If the application space range contains the space range of the current Tile map, the Tile map data corresponding to the current Tile map is converted into real scene three-dimensional data and then outputted;

[0012] If the application space range intersects with the space range of the current Tile map, the intersection area of the space range of the current Tile map and the application space range is obtained, the current Tile map is cut according to the intersection area, the Tile map data of the cutting area is retained, the Tile map data in the cutting area is converted into real scene three-dimensional data and then outputted, and a real scene three-dimensional data result is obtained.

[0013] Further, the acquisition of real three-dimensional data and the acquisition of the spatial range of the data includes:

[0014] The real three-dimensional data includes digital elevation model data, digital orthophoto data, digital surface model data, and oblique photography data.

[0015] The real three-dimensional data is stored in the file server in a data scale framing manner.

[0016] The acquisition of the spatial range of the data includes a preset range and administrative division information, the preset range includes administrative region location information described by text, and the administrative division information includes administrative region name and coordinate information of administrative region spatial range.

[0017] The spatial range is converted into split sentences by a semantic tool and returned in geojson format.

[0018] Further, the spatial range of the data is queried in the Tile map data, and the spatial range of each Tile map is acquired, including:

[0019] All point string information in the spatial range of the data is acquired, and the all point string information is all coordinate information in the administrative region spatial range.

[0020] The point string information is converted by a coordinate conversion formula as follows:

[0021] ;

[0022] ;

[0023] In the formula, Tile column number represents the real three-dimensional column number result to be queried, Tile row number represents the real three-dimensional row number result to be queried, Tile width is the longitude range covered by a single Tile at this level, which is determined by the scale; Lng represents the longitude information of the point to be input, Lng origin is the real three-dimensional origin longitude information stored in the metadata file, Tile' width is the latitude range covered by a single Tile' at this level, which is determined by the scale, Lng' represents the latitude information of the point to be input, Lng' origin is the real three-dimensional origin latitude information stored in the metadata file, and Z represents the zoom ratio of the real three-dimensional data.

[0024] The administrative region spatial range (minLon, minLat, maxLon, maxLat) is converted into the Tile column number and row number of the corresponding level, and the conversion formula is as follows:

[0025] ;

[0026] ;

[0027] ;

[0028] ;

[0029] wherein StartX / Y represents the starting row number and column number, EndX / Y represents the ending row number and column number, minLon is the minimum longitude, maxLon is the maximum longitude; minLat is the minimum latitude, maxLat is the maximum latitude, and Z is the zoom level.

[0030] Further, the positional relationship between the spatial range of each Tile map and the claimed spatial range is determined, including:

[0031] If all the vertices in the spatial range of the current Tile map are within the claimed spatial range, and the edges of the spatial range of the current Tile map do not intersect with the edges of the claimed spatial range, it is determined that the claimed spatial range contains the spatial range of the current Tile map.

[0032] If the lower-left corner coordinate (x1, y1) and the upper-right corner coordinate (x2, y2) of the claimed spatial range satisfy:

[0033] max(x1, a1) < min(x2, a2) & max(y1, b1) < min(y2, b2)

[0034] it is determined that the spatial range of the current Tile map intersects with the claimed spatial range.

[0035] Further, the intersection region of the spatial range of the current Tile map and the claimed spatial range is obtained by performing a clipping operation on the spatial range of the current Tile map and the claimed spatial range.

[0036] Further, the current Tile map is clipped according to the intersection region, including:

[0037] Converting the geojson data corresponding to the intersection region in the claimed spatial range into GeoRegion format data;

[0038] Performing coordinate conversion on the data after conversion to ensure that the coordinate system is consistent with that of the current Tile map;

[0039] Constructing a three-dimensional clipping space plane group according to the GeoRegion format;

[0040] Traverse all model vertices of the current Tile map, determine whether the vertex is inside the three-dimensional clipping space plane group;

[0041] Reserve the vertices inside the three-dimensional clipping space plane group, and construct a clipping region according to the vertices.

[0042] Further, according to the GeoRegion format, construct the three-dimensional clipping space plane group, including:

[0043] According to the Delaunay triangulation tool, convert the GeoRegion format data into a convex polygon;

[0044] According to the vertices of the convex polygon, calculate the minimum X, Y value and the maximum X, Y value, generate four vertical planes, which are the right plane, the left plane, the front plane and the back plane;

[0045] According to the elevation data in the digital elevation model data, calculate the maximum Z value and the minimum Z value, and generate the upper plane and the lower plane;

[0046] According to the right plane, the left plane, the front plane, the back plane, the upper plane and the lower plane, construct the three-dimensional clipping space plane group.

[0047] Further, it also includes:

[0048] Output the real three-dimensional data result to the specified path;

[0049] Automatically package the output real three-dimensional data result;

[0050] Use GDAL tool to complete the projection conversion of real three-dimensional data result, generate real three-dimensional data result cache configuration file, real three-dimensional data result cache, and add cache to three-dimensional data scene;

[0051] Save the map and workspace, publish the vector map service and three-dimensional data service, and the steps are as follows:

[0052] Connect the data service publishing platform using the username and password;

[0053] Get the type and data workspace of the real three-dimensional data result, and copy the data workspace path address;

[0054] According to the configuration file, automatically generate the metadata information of the real three-dimensional;

[0055] According to the type and workspace path address of the real three-dimensional data, call the data service publishing interface;

[0056] Register the map service and three-dimensional data service to the real three-dimensional GIS data platform automatically, complete the three-dimensional service superposition and display of the real three-dimensional data result.

[0057] The application further discloses a real-scene three-dimensional archive data application device.

[0058] The application device comprises:

[0059] The application device comprises:

[0060] The application device comprises:

[0061] The application device comprises:

[0062] The application device comprises:

[0063] The application device comprises:

[0064] The application further provides a readable storage medium, wherein the readable storage medium stores a computer program, and the computer program is executed by a processor to realize the steps of the real-scene three-dimensional archive data application method.

[0065] The application has the following beneficial effects:

[0066] The real-scene three-dimensional archive data application method, device and readable storage medium can accurately locate the space range and object category corresponding to the data applied by a user, replace the traditional manual screening operation, improve the intelligent extraction capability and spatial accuracy, automatically execute a hierarchical clipping strategy according to different data levels, security requirements or user roles, control the precision and divide the area of the real-scene three-dimensional data, support flexible configuration and on-demand delivery of data, improve the system resource utilization rate and release efficiency, construct a real-scene three-dimensional data release system linked with the clipping result, realize a full-process controllable data visualization and delivery mechanism, prevent unauthorized access and data leakage, and enhance the system data security and sensitive area protection capability. Thus, adaptive responses to different use roles, real-time requirements and desensitization levels are realized. Attached Figure Description

[0067] Figure 1 This is a flowchart illustrating a method, apparatus, and readable storage medium for applying for real-scene 3D archive data according to this application. Detailed Implementation

[0068] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some, not all, of the embodiments of the present invention. The components of the embodiments of the present invention described and shown in the accompanying drawings can generally be arranged and designed in various different configurations.

[0069] Therefore, the following detailed description of the embodiments of the invention provided in the accompanying drawings is not intended to limit the scope of the claimed invention, but merely to illustrate selected embodiments of the invention. All other embodiments obtained by those skilled in the art based on the embodiments of the invention without inventive effort are within the scope of protection of the invention.

[0070] The specific embodiments of the present invention will now be described in detail with reference to the accompanying drawings.

[0071] like Figure 1 As shown, a method for applying for real-scene 3D archive data includes:

[0072] S1: Obtain real-scene 3D data and the application space range for data application;

[0073] S2: Based on the Tile operation rules, the real-world 3D data is divided into sections to obtain Tile map sheet data;

[0074] S3: Query all Tile sheets in the Tile sheet data corresponding to the claimed space range, and obtain the space range of each Tile sheet;

[0075] S4: Determine the positional relationship between the spatial extent of each Tile sheet and the claimed spatial extent;

[0076] S5: If the claimed spatial range includes the spatial range of the current Tile map sheet, then convert the Tile map sheet data corresponding to the current Tile map sheet into real-scene 3D data and output it.

[0077] S6: If the application space range intersects with the space range of the current Tile map, obtain the intersection area of the space range of the current Tile map and the application space range, crop the current Tile map according to the intersection area, retain the Tile map data in the cropped area, convert the Tile map data in the cropped area into real scene three-dimensional data and output, and obtain the real scene three-dimensional data result.

[0078] Specifically, the application discloses a real scene three-dimensional archive data application method, which comprises the following steps:

[0079] Step 1: data acquisition: acquiring DEM (digital elevation model), DOM (digital orthophoto map), DSM (digital surface model), oblique photography and other real scene three-dimensional data needing data application, and storing the data in a file server in a data scale and map sheet manner.

[0080] Step 2: data application semanticization: through the data self-defined range of data application and the administrative division information of application, a semantic tool (such as text2sql) is used to convert into a split statement, the space range of application is obtained, and the geojson format is returned.

[0081] Step 3: calculating the map sheet range of data: inputting the geojson space to return the map sheet range information of real scene three-dimensional data, and the calculation steps are as follows:

[0082] Tile operation rule: real scene three-dimensional data is stored in a quadtree tiling structure, each Tile is recursively divided by 4 to the power of N, the higher the level (the larger the value), the lower the precision, and the smaller the data amount (for example: L1 level (original precision) → L2 level (4 sub-blocks) → L3 level (16 sub-blocks)), the Tile is quickly positioned through a quadtree index, a parent node contains four child nodes, and the row and column numbers of the child nodes are twice the row and column numbers of the parent node (for example, the child nodes of the parent node Tile_1_1_1 are Tile_2_2_2, Tile_2_3_2 and the like).

[0083] According to the application space range of step 2, the Tile row and column numbers of the corresponding map sheet data are inquired.

[0084] 2.1, all point string information in the geojson range is obtained, and coordinate conversion is performed (taking longitude as an example)

[0085] ;

[0086] In the formula, Tile column number represents the real three-dimensional column number result to be queried, Tile width is the longitude range covered by a single Tile at this level, which is determined by the scale; Lng represents the longitude information of the input point, and the Lng origin is the real three-dimensional origin longitude information stored in the metadata file (such as metadata.xml); Z represents the zoom scale of the real three-dimensional data.

[0087] The latitude calculation coordinate conversion is the same, and the longitude parameter in the formula is replaced, and the calculation formula is as follows:

[0088] ;

[0089] 2.2, range calculation, input geographic range (minLon, minLat, maxLon, maxLat), zoom level Z, convert minLon / minLat and maxLon / maxLat into the starting Tile column number and row number and the ending Tile column number and row number of the corresponding level:

[0090] ;

[0091] ;

[0092] ;

[0093] ;

[0094] In the formula, EndX / Y represents the starting row number and column number, EndX / Y represents the ending row number and column number, minLon is the minimum longitude, maxLon is the maximum longitude; minLat is the minimum latitude, maxLat is the maximum latitude, and Z is the zoom level.

[0095] 2.3, traverse the row and column number range [StartX, EndX] × [StartY, EndY], and generate all Tile row and column number combinations.

[0096] Step 4: Calculate the map graph relationship: obtain all the row and column numbers, traverse the spatial range information of each Tile map obtained in step 3, calculate the position relationship with the current range, and the calculation steps are as follows:

[0097] 4.1, inclusion relationship judgment: the space range surface A contains the current Tile map range surface B, then the formula is as follows:

[0098] A⊇B⟺∀P∈B,P∈A;

[0099] Need to determine whether all vertices of B are in A, and B edge and A edge no intersection (or only in the boundary contact).

[0100] 4.2, intersection relationship judgment: if the left lower corner (x1, y1) and the right upper corner (x2, y2) of the range of the space of the application and the left lower corner (a1, b1) and the right upper corner (a2, b2) of the range of the Tile map meet:

[0101] max(x1, a1) < min(x2, a2) & max(y1, b1) < min(y2, b2);

[0102] Then judge that the range of the space of the application and the current Tile map range intersect.

[0103] Contain: directly output the corresponding real scene three-dimensional data of step 1 corresponding to the range data of the current single tile map.

[0104] 4.3, map intersection space range calculation, wherein the calculation process is as follows:

[0105] 1. Perform merge analysis between the range of the space of the application and the range of the current Tile map, calculate the intersection point (do the clipping operation);

[0106] 2, output the node and chain, and construct the coordinate string information of the viewable area surface;

[0107] 3, establish topology and new object, and generate surface object;

[0108] 4, remove a large number of broken polygons, and fuse similar polygons;

[0109] 5, repeat steps 1-4 to get the intersection area of the viewable area.

[0110] Step 5: Real scene three-dimensional automatic clipping: the map information of the containing relationship does not participate in the clipping calculation, and the intersection space range generates two-dimensional data participating in the real scene three-dimensional data clipping, and the calculation steps are as follows:

[0111] 5.1, convert the range geojson data generated by the intersection into GeoRegion format;

[0112] 5.2, use proj4js or GDAL tool for coordinate conversion, and keep consistent with the Tile map coordinate system.

[0113] 5.3, define a three-dimensional clipping space plane group, and the calculation steps are as follows:

[0114] Calculate the convex hull: use the Delaunay triangulation tool to convert the irregular polygon of GeoRegion into a convex polygon.

[0115] Bounding box plane generation: Calculate the minimum / maximum X, Y values according to the convex hull vertices, and generate vertical planes. The formula of the vertical plane is as follows:

[0116]

[0117] Where: (A, B, C) is the normal vector of the plane, which determines the orientation of the plane. D is the offset of the plane to the origin, which determines the position of the plane. If the result is positive, it means that the point is on the side of the plane normal vector (reserved); if it is negative, it means that the point is on the opposite side of the plane; if it is 0, it means that it is located on the plane, x, y, x are the longitude, latitude and elevation information of the convex hull vertices.

[0118] According to the minimum / maximum X / Y values calculated from the convex hull vertices, four vertical planes can be obtained:

[0119] ;

[0120] From top to bottom, they are the right plane, the left plane, the front plane, and the back plane.

[0121] Generate height constraint planes. According to the maximum and minimum values of Z in the metadata file of the real scene three-dimensional (such as metadata.xml), the upper and lower planes can be generated, and the results are as follows:

[0122] ;

[0123] 5.4, Traverse all model vertices of the Tile map, and use OpenSceneGraph (OSG) to determine whether they are inside the vertical clipping planes and the upper and lower planes.

[0124] 5.5, Face recombination: Rebuild the topological relationship of the retained faces to avoid breaking caused by clipping. The steps are as follows:

[0125] (1) Extract the inside vertices of the retained model;

[0126] (2) Regenerate the triangular mesh (use Delaunay triangulation or incremental method);

[0127] (3) Repair the texture coordinates (use UV mapping);

[0128] 5.6, Configure the output path of the clipping result;

[0129] 5.7, Configure the data name of the clipping result, and the naming rule is Tile map number + data name + clipping timestamp;

[0130] 5.8, Input the input file, clipping range, output path and other parameters, retain the data in the clipping area, and exclude external data;

[0131] 5.9, Real scene three-dimensional result data post-processing: use LOD optimization on the cropped model data to generate multiple levels of detail, balance rendering performance and accuracy;

[0132] 5.10, Convert the data to real scene three-dimensional data (such as osgb) format and output to the specified path.

[0133] Step 6: Automatic packaging of the data results, the steps are as follows:

[0134] 6.1, Create a Zip output stream: specify the path of the compressed ZIP file;

[0135] 6.2, Recursively traverse the real scene three-dimensional data folder: get all files and subfolders in the folder;

[0136] 6.3, Add files to the Zip stream: convert each file to a ZipEntry and write it to the compressed package.

[0137] Step 7: Push the path and basic information of the real scene three-dimensional data results to the data processing system.

[0138] Step 8: Use GDAL tools to complete the real scene three-dimensional data projection conversion.

[0139] Step 9: Generate real scene three-dimensional data cache configuration file.

[0140] Step 10: Generate real scene three-dimensional data cache and add cache to three-dimensional data scene.

[0141] Step 11: Save the map and workspace, publish vector map service and three-dimensional data service. The steps are as follows:

[0142] 11.1, Connect to the data service publishing platform using a username and password;

[0143] 11.2, Get the real scene three-dimensional data type workspace location and copy the data workspace path address;

[0144] 11.3, Automatically generate real scene three-dimensional metadata information according to the configuration file;

[0145] 11.4, According to the real scene three-dimensional data type and working path, call the data service publishing interface.

[0146] Step 12: Automatically register the service to the real scene three-dimensional GIS data platform, complete the three-dimensional service superposition and display of the real scene three-dimensional data results.

[0147] The real scene three-dimensional archive data application method, device and readable storage medium can accurately position the space range and object category corresponding to the data applied by the user, replace the traditional manual screening operation, improve the intelligent extraction capability and spatial accuracy; according to different data levels, security requirements or user roles, automatically execute a hierarchical clipping strategy, control the precision and region block of the real scene three-dimensional data, support flexible configuration and on-demand delivery of data, improve system resource utilization and release efficiency; construct a real scene three-dimensional data release system linked with the clipping result, realize a full-process controllable data visualization and delivery mechanism, prevent unauthorized access and data leakage, and enhance the system data security and sensitive area protection capability. The steps of user request submission, semantic processing, data clipping, release notification and the like are integrated to form an automatic processing chain without manual intervention, realize an automatic closed loop of "user initiation-system response-data delivery", and significantly improve the efficiency and service response speed.

[0148] The above only describes the preferred embodiments of the present application, and it should be noted that those skilled in the art can make several improvements and refinements without departing from the technical principles of the present application, and these improvements and refinements should also be considered within the protection scope of the present application.

Claims

1. A method for real-time three-dimensional archival data retrieval, characterized by, The application relates to a method for outputting real three-dimensional data based on a data request, and belongs to the field of real three-dimensional data processing. The method comprises the following steps: acquiring real three-dimensional data and a request space range of a data request; performing framing processing on the real three-dimensional data based on a Tile operation rule to obtain Tile map data; inquiring all Tile maps corresponding to the request space range in the Tile map data and acquiring the space range of each Tile map; judging the positional relationship between the space range of each Tile map and the request space range; if the request space range contains the space range of the current Tile map, the Tile map data corresponding to the current Tile map is converted into real three-dimensional data and then outputted; if the request space range intersects with the space range of the current Tile map, the intersection area of the space range of the current Tile map and the request space range is acquired, the current Tile map is cropped according to the intersection area, the Tile map data in the cropped area is reserved, the Tile map data in the cropped area is converted into real three-dimensional data and then outputted, and real three-dimensional data results are obtained; The method comprises the following steps: acquiring real three-dimensional data, wherein the real three-dimensional data comprises digital elevation model data, digital orthophoto map data, digital surface model data and oblique photography data; storing the real three-dimensional data in a file server in a framing manner according to a data scale; acquiring a request space range of a data request, wherein the request space range comprises a preset range and administrative division information, the preset range comprises administrative region position information described by text, and the administrative division information comprises administrative region names and coordinate information of administrative region space ranges; converting the request space range into split statements by using a semantic tool and returning the split statements in a geojson format; inquiring all Tile maps corresponding to the request space range in the Tile map data and acquiring the space range of each Tile map, which comprises the following steps: acquiring all point string information in the request space range, wherein the point string information is all coordinate information in the administrative region space range; ； ； In the formula, the Tile column number represents a real three-dimensional column number result to be queried, the Tile row number represents a real three-dimensional row number result to be queried, the Tile width is a longitude range covered by a single Tile at a corresponding level, and is determined by a scale; the Lng represents longitude information of an input point, the Lng origin is real three-dimensional origin longitude information stored in a metadata file, the Tile' width is a latitude range covered by a single Tile' at a corresponding level, and is determined by a scale, represents input point latitude information, is real three-dimensional origin latitude information stored in a metadata file, and Z represents a scaling ratio of real three-dimensional data. performing coordinate conversion on the point string information, and the coordinate conversion formula is as follows: ; ; ; ; converting the administrative region space range (minLon, minLat, maxLon, maxLat) into Tile column numbers and row numbers of a corresponding level, and the conversion formula is as follows: In the formula, StartX / Y represents the starting row number and column number, EndX / Y represents the ending row number and column number, minLon is the minimum longitude, maxLon is the maximum longitude, minLat is the minimum latitude, maxLat is the maximum latitude, and Z is the zoom level. cropping the current Tile map according to the intersection area, which comprises the following steps: converting geojson data corresponding to the intersection area in the request space range into GeoRegion format data; performing coordinate conversion on the converted data to make the converted data consistent with the coordinate system of the current Tile map; constructing a three-dimensional cropping space plane group according to the GeoRegion format; traversing all model vertices of the current Tile map and judging whether the vertices are inside the three-dimensional cropping space plane group. Reserve the vertices located inside the three-dimensional clipping space plane group, and form a clipping area according to the vertices; According to the GeoRegion format, build a three-dimensional clipping space plane group, including: According to the Delaunay triangulation tool, convert the GeoRegion format data into a convex polygon; According to the vertices of the convex polygon, calculate the minimum X, Y values and the maximum X, Y values, generate four vertical planes, which are the right plane, the left plane, the front plane and the back plane; According to the elevation data in the digital elevation model data, calculate the maximum Z value and the minimum Z value, and generate the upper plane and the lower plane; According to the right plane, the left plane, the front plane, the back plane, the upper plane and the lower plane, build a three-dimensional clipping space plane group.

2. A method of issuing a real scene three-dimensional archive data according to claim 1, characterized in that, Determine the spatial range of each Tile map and the positional relationship of the claimed spatial range, including: If all the vertices in the spatial range of the current Tile map are within the claimed spatial range, and the edges of the spatial range of the current Tile map do not intersect with the edges of the claimed spatial range, it means that the claimed spatial range contains the spatial range of the current Tile map; If the lower left corner coordinate (x1, y1) and the upper right corner coordinate (x2, y2) of the claimed spatial range satisfy: max(x1, a1) < min(x2, a2) & max(y1, b1) < min(y2, b2) Then determine that the spatial range of the current Tile map and the claimed spatial range intersect.

3. A method of issuing a real scene three-dimensional archive data according to claim 1, characterized in that, By performing a clipping operation on the spatial range of the current Tile map and the claimed spatial range, the intersection area of the spatial range of the current Tile map and the claimed spatial range is obtained.

4. A method of claiming real scene three-dimensional archival data according to claim 1, wherein, Further comprising: Output the real three-dimensional data result to a specified path; Automatically package the output real three-dimensional data result; Use the GDAL tool to complete the projection conversion of the real three-dimensional data result, generate a real three-dimensional data result cache configuration file, a real three-dimensional data result cache, and add the cache to the three-dimensional data scene; Save the map and workspace, publish the vector map service and three-dimensional data service, and the steps are as follows: Connect the data service publishing platform using a username and password; Get the type and data workspace of the real three-dimensional data result, and copy the data workspace path address; According to the configuration file, automatically generate the metadata information of the real three-dimensional; According to the real three-dimensional data type and the workspace path address, call the data service publishing interface; Automatically register the map service and the three-dimensional data service to the real three-dimensional GIS data platform, complete the three-dimensional service superposition and display of the real three-dimensional data result.

5. A live three-dimensional archival data retrieval device, characterized by, Including: An acquisition module for acquiring real three-dimensional data and a claimed spatial range of data; A framing module for performing framing processing on the real three-dimensional data based on the Tile operation rule to obtain Tile map data; A query module for querying all Tile maps corresponding to the claimed spatial range in the Tile map data, and obtaining the spatial range of each Tile map; A judging module is configured to judge the spatial range of each Tile map and the positional relationship between the spatial range and the claimed spatial range; A first output module is configured to convert the Tile map data corresponding to the current Tile map into real scene three-dimensional data and then output the real scene three-dimensional data if the claimed spatial range contains the spatial range of the current Tile map; A second output module is configured to obtain an intersection region between the spatial range of the current Tile map and the claimed spatial range, crop the current Tile map according to the intersection region, retain the Tile map data in the cropped region, convert the Tile map data in the cropped region into real scene three-dimensional data, and then output the real scene three-dimensional data to obtain real scene three-dimensional data results if the claimed spatial range intersects with the spatial range of the current Tile map; The real scene three-dimensional data and the claimed spatial range of the data claim are obtained, including: The real scene three-dimensional data is obtained, and the real scene three-dimensional data includes digital elevation model data, digital orthographic image data, digital surface model data, and oblique photography data; The real scene three-dimensional data is stored in the file server in a way of being divided into frames according to a data scale; The claimed spatial range of the data claim is obtained, and the claimed spatial range includes a preset range and administrative division information, the preset range includes administrative region position information described by text, and the administrative division information includes the name of the administrative region and the coordinate information of the spatial range of the administrative region; The claimed spatial range is converted into split statements by a semantic tool and returned in a geojson format; All Tile maps corresponding to the claimed spatial range in the Tile map data are queried, and the spatial range of each Tile map is obtained, including: All point string information in the claimed spatial range is obtained, and the all point string information is all coordinate information in the spatial range of the administrative region; The point string information is converted in coordinates, and the conversion formula is as follows: ； ； In the formula, the Tile column number represents a real three-dimensional column number result to be queried, the Tile row number represents a real three-dimensional row number result to be queried, the Tile width is a longitude range covered by a single Tile at a corresponding level, and is determined by a scale; the Lng represents longitude information of an input point, the Lng origin is real three-dimensional origin longitude information stored in a metadata file, the Tile' width is a latitude range covered by a single Tile' at a corresponding level, and is determined by a scale, represents latitude information of an input point, is real three-dimensional origin latitude information stored in a metadata file, and Z represents a scaling ratio of real three-dimensional data. The spatial range of the administrative region (minLon, minLat, maxLon, maxLat) is converted into the Tile column number and the row number of the corresponding level, and the conversion formula is as follows: ; ; ; ; In the formula, StartX / Y represents the starting row number and column number, EndX / Y represents the ending row number and column number, minLon represents the minimum longitude, maxLon represents the maximum longitude, minLat represents the minimum latitude, maxLat represents the maximum latitude, and Z represents the zoom level. The current Tile map is cropped according to the intersection region, including: The geojson data corresponding to the intersection region in the claimed spatial range is converted into GeoRegion format data; The converted data is converted in coordinates to be consistent with the coordinate system of the current Tile map; A three-dimensional cropping space plane group is constructed according to the GeoRegion format; All model vertices of the current Tile map are traversed, and it is judged whether the vertices are inside the three-dimensional cropping space plane group; The vertices inside the three-dimensional cropping space plane group are retained, and the cropping region is formed according to the vertices; The three-dimensional cropping space plane group is constructed according to the GeoRegion format, including: The GeoRegion format data is converted into a convex polygon according to a Delaunay triangulation tool; According to the vertexes of the convex polygon, minimum X and Y values and maximum X and Y values are calculated, and four vertical planes are generated, which are a right plane, a left plane, a front plane and a back plane; According to the elevation data in the digital elevation model data, maximum and minimum Z values are calculated, and an upper plane and a lower plane are generated; A three-dimensional clipping space plane group is constructed according to the right plane, the left plane, the front plane, the back plane, the upper plane and the lower plane.

6. A readable storage medium characterized by, The readable storage medium has a computer program stored thereon, and the computer program is executed by the processor to implement the steps of the real scene three-dimensional archive data application method according to any one of claims 1 to 4.

Images