CN113239721A - Road ground collapse analysis display method - Google Patents

Abstract

The invention discloses a road ground collapse analysis display method. The method includes step S1: making a shape format file; s2: performing spatial operation; s3: storing data: storing the geoelectrical section distribution graphic results corresponding to each detection line, namely an equal apparent resistivity section diagram and an electrical geological section diagram on a hard disk medium; s4: tiling and dividing storage: dividing each ground electric profile map into a plurality of tile grades; s5: and (4) displaying the ground fault section result diagram. The road ground collapse analysis display method disclosed by the invention is used for realizing rapid and convenient dynamic query analysis and display from the urban road ground collapse survey grating result graph at the client by combining the support of a Java programming language and a PostgreSQL database on geospatial vector data and expanding an SQL support function realized on a DEM-9IM spatial algorithm by utilizing the PostGIS database.

Description

Road ground collapse analysis display method

Technical Field

The invention relates to the technical field of electronic geographic information processing, in particular to a road ground collapse analysis display method.

Background

In the raster result map of urban road surface collapse survey, because the result map is usually very long, it is difficult for a computer to conveniently and quickly display a picture with such a size, and the result map is even more displayed on a webpage end or a mobile end.

Disclosure of Invention

The invention aims to overcome the defects in the prior art and provide a road ground collapse analysis and display method. The user access passage ground collapse analysis display method can check details of the geoelectricity section result diagram at any position, and is convenient and rapid.

The technical scheme of the invention is as follows:

the road ground collapse analysis display method comprises the following steps: the method comprises the following steps:

s1: making a shape format file: sequentially acquiring data on site according to a line acquisition plan by adopting a high-density resistivity method; processing the geophysical prospecting engineering detection circuit plane graph into a linear shape shapefile format file by ArcGIS software;

s2: spatial operation: importing the shape format file into a PostGIS database, and carrying out spatial query operation and search by a spatial SQL operation function of the PostGIS according to conditions;

s3: storing data: storing the geoelectrical section distribution graphic results corresponding to each detection line, namely an equal apparent resistivity section diagram and an electrical geological section diagram on a hard disk medium;

s4: tiling and dividing storage: dividing each ground electric profile map into a plurality of tile grades;

s5: and (4) displaying the ground fault section result diagram.

In the step S4, dividing each ground electrical section distribution map into six tile grades, which are 0,1, 2, 3, 4 and 5 in sequence; and each tile level is sliced by vertical pixel length/256 x horizontal pixel length/256 tiles count.

Further, when 256px is less, the white fill is performed.

Further, step S5 specifically includes the following steps:

s51: taking a geophysical prospecting engineering detection circuit plan as a base map, and superposing a layer of geoelectricity section result slice display map at the right lower part of the geophysical prospecting engineering detection circuit plan;

s52: and a user can check the ground electric fracture result diagrams of other positions in the geophysical prospecting engineering control circuit plan by dragging the transverse slider.

The invention has the beneficial effects that: the road ground collapse analysis display method disclosed by the invention is used for realizing rapid and convenient dynamic query analysis and display from the urban road ground collapse survey grating result graph at the client by combining the support of a Java programming language and a PostgreSQL database on geospatial vector data and expanding an SQL support function realized on a DEM-9IM spatial algorithm by utilizing the PostGIS database.

Drawings

FIG. 1 is a plan view of a geophysical prospecting engineering survey circuit of the present invention;

FIG. 2 is a schematic view of a tile cut of the present invention;

FIG. 3 is a plan view of a geophysical prospecting engineering detection circuit and a ground electric section result slice comprehensive display effect diagram;

FIG. 4 is a slice display effect diagram of ground fault section result of the present invention.

Detailed Description

For a better understanding of the invention, reference will now be made to the following examples and accompanying drawings.

As shown in fig. 1 to 4, the method for analyzing and displaying the road surface collapse comprises the following steps:

(1) sequentially acquiring data on site according to a line acquisition plan by adopting a high-density resistivity method; processing the geophysical prospecting engineering detection circuit plane graph into a linear shape shapefile format file by ArcGIS software; each line is formed by connecting a plurality of points (x, y) in sequence, as shown in fig. 1 below: "Geophysic #1, Geophysic #2, Geophysic #3, Geophysic # 4" indicates that a total of 4 lines are collected, wherein "Geophysic # 1" is formed by sequentially connecting points [ (10,10) (20,10) (30,10) ]: the property fields of the shape file at least include but are not limited to the following fields: id. code, name, geo.

Table 1: exploration drilling record point layer field table

id	code	name	geom
				1	gmddz1001	Geophysical prospecting wire-Geogaphic #1	LineString(10 10,20 10,30 10)
2	gmddz1002	Geophysical prospecting wire-Geogaphic #2	LineString(10 20,20 20,30 20)
				3	gmddz1003	Geophysical prospecting wire-Geogaphic #3	LineString(40 50,50 50,60 50)
4	gmddz1004	Geophysical prospecting wire-Geogaphic #4	LineString(40 60,50 60,60 60)

Description of the main fields:

id-self-increment id;

code-a geophysical unique code;

name-geophysical prospecting circuit logical name;

a geo-geophysical route geometry field, describing location information of the route on a two-dimensional map.

(2) And importing the shapefile format file into a PostGIS database, and conveniently and quickly performing space query operation and search according to various conditions by means of the PostGIS space SQL operation function capability.

(3) Storing the geological profile distribution graphic results (namely an iso-apparent resistivity profile and an electrical geological profile) corresponding to each detection line on a computer hard disk medium by taking the corresponding code as a name, wherein each graphic result corresponds to a folder named by the code; since each geographical profile graphic result picture is very large (about 50MB or so) in most cases, and it is difficult to load and display quickly at a network client with very limited bandwidth, it is necessary to tile and store this graphic result by levels, and to load the result by levels and demands in a tiled manner to increase the loading speed. Therefore, it is specified that each map of the terrain profile is divided into 6 tile levels, which are 0,1, 2, 3, 4, and 5 in turn, and each tile level is sliced by the number of "vertical pixel length/256 × horizontal pixel length/256" tiles, and those with less than 256px are filled with white, as shown in fig. 2 (3,0), (3,1) two slices. For convenience of description, we can use two-dimensional rows (rows) and columns (cols), i.e. a total of rows × cols slices: it is specified that, from level 0 to level 5, the resolution size of the ground electric profile distribution result graph is correspondingly doubled every last level, and assuming that the size of level 0 is 896 × 512 pixels, the size of level 1 is 1792 × 1024pixels, so that the size calculation formula of any level can be obtained based on the size conversion of level 0, that is: (896 × (1+ level)) × (512 × (1+ level)) pixels. As a class 5 sizing formula:

(896 × (1+5)) × (512 × (1+5)) ═ 5376 × 3072 pixcels; for the convenience of analysis during client imaging, a description file which is specified by a standard JSON format and has a file name of tiles is established in a corresponding folder which stores the tile slicing rules by taking code as a name and recorded, and the content of the file is detailed as follows:

description of the meaning of each field:

tilename-tile name, used to describe which geophysical prospecting line this is;

tileformat-slice format, MIME type, fixed as image/jpg;

basesize-geoelectrical section achievement diagram 0 level resolution reference size, unit is pixel;

tilesize-slice resolution size, fixed at 256 × 256 pixels;

tiles-hierarchical description, describing the number of columns and rows per level;

level-level of each level;

cols-number of columns per stage;

rows-the number of rows per stage;

(4) the relationship between the actual length of the geophysical prospecting wire and the ground fault profile distribution result diagram takes a code of a 'gmddz 1001' circuit as an example, and the actual length of the circuit is 896 meters, and the transverse length of the corresponding 0-level ground fault profile distribution result diagram is 896 pixels. Then, the actual length represented per pixel at level 0 is: 896 m/896 m pixels 1 m, from which the formula representing the actual length per pixel is given by: actual mileage divided by pixel length; the actual length represented per pixel at 5 levels is: 896 m ÷ 5376 pixels ≈ 0.1667 m; loading the geophysical links stored in the PostGIS database onto a client map, and then enabling a user to select any geophysical link (taking 'Geographic # 1' as an example) and quickly and conveniently view a corresponding geoelectrical section result diagram according to the following method:

a. as shown in fig. 3, a geophysical prospecting engineering detection circuit plan is used as a base map, and then a layer of result slice display map is superimposed on the lower right;

b. as shown in fig. 4, the production slice display area superimposed on the lower right is fixed as a three-column 256 × 256-pixel slice, the number of rows is obtained according to the rows field parameter of the corresponding level in the tiles json record file, but the displayable area is fixed as only two rows can be seen, and other rows that cannot be seen can be viewed by scrolling through the vertical scroll bar and are defaulted as 0 level. A sliding bar capable of sliding transversely is arranged below the probe line, the sliding bar represents the length of the geophysical line (for example, the geophysical line 'Geogaphic # 1' is 896 meters), and the sliding bar is positioned at the beginning of the geophysical line, namely 0 meter by default;

c. the result slice display image area superposed at the lower right is fixed into a three-column 256 × 256-pixel slice, the row number is obtained according to the rows field parameter of the corresponding level in the tiles json record file, but the displayable area is fixed to be only capable of seeing two rows, other rows which cannot be seen can be viewed through rolling of a vertical scroll bar, and the default is 0 level. A sliding bar capable of sliding transversely is arranged below the position, and the sliding bar represents the length of the geophysical line (for example, the geophysical line 'Geogaphic # 1' is 896 meters), and is positioned at the beginning of the geophysical line by default, namely 0 meter (as shown in figure 4);

d. when the user drags the horizontal slider, for example, to 10 meters, the corresponding geophysical line "Geographic # 1" on the map also synchronously shows a red dot representing the current viewing position at 10 meters. Meanwhile, according to parameters such as the horizontal axis 896pixels of the resolution reference dimension corresponding to the current level 0 and the level 0, the actual length 896 meters of the object detection line, the size 256 × 256pixels of the slice and the like, the number of the corresponding slice to be displayed is calculated, and the calculation formula is as follows: 896 pixels/896 m × 10 m 10 pixels, i.e. the slice at 10 m is at 10 pixels transverse to the production map; 10 pixels/256 is 0 or more than 10, i.e. the pixels in the result map 10 are located in the slice with 0 rows; according to the above specification of the area of the result slice display image, the result slice display image is fixed into three columns, three columns of 0,1 and 2 should be taken out respectively, and then according to the rows field parameter corresponding to level 0 of the tiles json record file, the number of rows is 2, so that the number of the last slice to be displayed should be: (0,0), (0,1), (1,0), (1,1), (2,0), (2,1), after calculating the serial number of the slice to be displayed, the slice of the corresponding serial number can be put into the corresponding position of the display picture area of the result slice; therefore, no matter where the user wants to see, the user only needs to drag the slider to the position where the user wants to see, and then the slice corresponding to the result picture can be called without calling the whole result picture at one time, and the method is convenient and fast.

The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the invention, so any modifications, equivalents or improvements made within the spirit and principles of the present invention should be included in the scope of the present invention.

Claims (4)

1. The road ground collapse analysis display method comprises the following steps: the method is characterized in that: the method comprises the following steps:

s1: making a shape format file: sequentially acquiring data on site according to a line acquisition plan by adopting a high-density resistivity method; processing the geophysical prospecting engineering detection circuit plane graph into a linear shape shapefile format file by ArcGIS software;

s2: spatial operation: importing the shape format file into a PostGIS database, and carrying out spatial query operation and search by a spatial SQL operation function of the PostGIS according to conditions;

s3: storing data: storing the geoelectrical section distribution graphic results corresponding to each detection line, namely an equal apparent resistivity section diagram and an electrical geological section diagram on a hard disk medium;

s4: tiling and dividing storage: dividing each ground electric profile map into a plurality of tile grades;

s5: and (4) displaying the ground fault section result diagram.

2. The road ground collapse analysis display method according to claim 1: the method is characterized in that: in the step S4, dividing each ground electric section distribution map into six tile grades, which are 0,1, 2, 3, 4 and 5 in sequence; and each tile level is sliced by vertical pixel length/256 x horizontal pixel length/256 tiles count.

3. The road ground collapse analysis and display method according to claim 2: the method is characterized in that: when the amount of the additive is less than 256px, the additive is performed.

4. The road ground collapse analysis display method according to claim 1: the method is characterized in that: the step S5 includes the following steps:

s51: taking a geophysical prospecting engineering detection circuit plan as a base map, and superposing a layer of geoelectricity section result slice display map at the right lower part of the geophysical prospecting engineering detection circuit plan;

s52: and a user can check the ground electric fracture result diagrams of other positions in the geophysical prospecting engineering control circuit plan by dragging the transverse slider.

Images