CN112860832A - Cable display method, device, equipment and storage medium for three-dimensional map - Google Patents

Cable display method, device, equipment and storage medium for three-dimensional map Download PDF

Info

Publication number
CN112860832A
CN112860832A CN202110129852.0A CN202110129852A CN112860832A CN 112860832 A CN112860832 A CN 112860832A CN 202110129852 A CN202110129852 A CN 202110129852A CN 112860832 A CN112860832 A CN 112860832A
Authority
CN
China
Prior art keywords
cable
image data
data
geological
dimensional
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202110129852.0A
Other languages
Chinese (zh)
Inventor
刘玉宝
曾振达
汪进锋
刘培杰
黄智鹏
王东芳
廖镇尧
黄杨珏
赖舷
孙新宇
叶杭
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guangdong Power Grid Co Ltd
Heyuan Power Supply Bureau of Guangdong Power Grid Co Ltd
Original Assignee
Guangdong Power Grid Co Ltd
Heyuan Power Supply Bureau of Guangdong Power Grid Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Guangdong Power Grid Co Ltd, Heyuan Power Supply Bureau of Guangdong Power Grid Co Ltd filed Critical Guangdong Power Grid Co Ltd
Priority to CN202110129852.0A priority Critical patent/CN112860832A/en
Publication of CN112860832A publication Critical patent/CN112860832A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/29Geographical information databases

Abstract

The embodiment of the invention provides a cable display method, a cable display device, equipment and a storage medium of a three-dimensional map, wherein the method comprises the following steps: loading three-dimensional body data of a specified area; drawing geological image data for the underground space of the area according to radar data generated when the ground penetrating radar detects the cable in the area; detecting an inflection point in the geological image data; sequentially connecting inflection points to the cable image data to represent the cable; writing cable information of the cable into three-dimensional body data of the region by referring to the cable image data to generate a semantic layer; and loading the three-dimensional body data to display the semantic layer of the cable in a three-dimensional map of the region. Compared with the traditional manual marking and manual statistics, the embodiment of the invention provides a more accurate and visual cable display method, and reduces the cable identification error.

Description

Cable display method, device, equipment and storage medium for three-dimensional map
Technical Field
The embodiment of the invention relates to the technical field of power engineering, in particular to a cable display method, device, equipment and storage medium for a three-dimensional map.
Background
Nowadays, with the rapid development of city construction, after entering modern society, due to the reasons of shortage of urban land, high traffic pressure, city appearance construction and the like, an underground cable power transmission mode is generally adopted in large cities. Especially in coastal areas, typhoons are frequently generated, overhead lines are easily damaged, and buried cables perfectly solve the problem. Compared with an overhead line, the cable has the advantages of small occupied area, reliable power transmission, strong anti-interference capability and the like.
In the existing scheme, related workers mark the signboard in the cable area as the position of the cable mainly by means of manual marking and manual statistics one by one, and the method is easy to be affected by external factors, the signboard is damaged manually or omitted and lost, and finally the correctness of the cable position is questioned.
Disclosure of Invention
The embodiment of the invention provides a cable display method, a cable display device, cable display equipment and a cable display storage medium of a three-dimensional map, and aims to solve the problem of realizing cable visualization in the three-dimensional map.
In a first aspect, an embodiment of the present invention provides a cable display method for a three-dimensional map, including:
loading three-dimensional body data of a specified area;
drawing geological image data for the underground space of the region according to radar data generated when the ground penetrating radar detects the cable in the region;
detecting an inflection point in the geological image data;
sequentially connecting the inflection points to the cable image data to represent the cable;
writing the cable information of the cable into the three-dimensional body data of the region by referring to the cable image data to generate a semantic layer;
and loading the three-dimensional body data to display the semantic layer of the cable in a three-dimensional map of the region.
Optionally, the mapping geological image data of the underground space of the region according to radar data generated when the ground penetrating radar detects the cable in the region includes:
acquiring radar data obtained by detecting the area pre-buried with the cable by the ground penetrating radar;
adjusting the radar data into two-dimensional geological image data;
determining a zero line in the geological image data, wherein the zero line is the joint position of a cable placed on the ground and a measuring line;
converting the geological image data from a time domain to a frequency spectrum domain;
and executing preset image processing on the geological image data.
Optionally, before detecting an inflection point in the geological image data, the method further includes:
and performing wavelet denoising of at least two levels on the geological image data.
Optionally, performing at least two levels of wavelet denoising on the geological image data, including:
performing wavelet decomposition on the address image data to obtain a first candidate signal and a second candidate signal of at least two levels, wherein the frequency of the second candidate signal is higher than that of the first candidate signal;
removing signals smaller than a threshold value from the second candidate signal of each level to obtain a third candidate signal, wherein the threshold value corresponds to each level;
reconstructing the first candidate signal and the third candidate signal of the highest rank as the geological image data.
Optionally, detecting an inflection point in the geological image data comprises:
translating the geological image data by a preset window, and calculating gray change information by the following formula:
Figure BDA0002925083210000031
wherein E (u, v) is gray scale change information, w (x, y) is a window, I (x + u, y + v) is gray scale information after the window is translated, and I (x, y) is gray scale information of the geological image data;
an inflection point is determined based on the gray-scale variation information.
Optionally, the image processing comprises at least one of:
gain processing, IIR filtering, FIR filtering, arithmetic operation, deconvolution, offset processing, static correction.
Optionally, writing the cable information of the cable into the three-dimensional shape data of the area with reference to the cable image data to generate a semantic layer, including:
creating an index in the cable image data;
inquiring an index value meeting the inquiry condition according to the index, and extracting the cable information of the cable;
writing the cable information into the three-dimensional body data of the region;
and generating a semantic layer according to the cable information in the three-dimensional body data.
In a second aspect, an embodiment of the present invention further provides a cable display apparatus for a three-dimensional map, including:
the three-dimensional body data loading module is used for loading three-dimensional body data of a specified area;
the geological image data drawing module is used for drawing geological image data for the underground space of the region according to radar data generated when the ground penetrating radar detects the cable in the region;
the inflection point detection module is used for detecting an inflection point in the geological image data;
a cable generation module for sequentially connecting the inflection points to the cable image data to represent the cable;
the semantic layer generating module is used for writing the cable information of the cable into the three-dimensional body data of the region by referring to the cable image data so as to generate a semantic layer;
and the cable display module is used for loading the three-dimensional body data so as to display the semantic map layer of the cable in a three-dimensional map of the region.
Optionally, the geological image data drawing module includes:
the radar data acquisition submodule is used for acquiring radar data obtained by detecting the ground penetrating radar to the area pre-buried with the cable;
the two-dimensional address image data acquisition submodule is used for adjusting the radar data into two-dimensional geological image data;
the zero line submodule is used for determining a zero line in the geological image data, and the zero line is the position of the connection between a cable placed on the ground and a measuring line;
the frequency spectrum domain conversion sub-module is used for converting the geological image data from a time domain into a frequency spectrum domain;
and the image processing submodule is used for executing preset image processing on the geological image data.
Optionally, before detecting an inflection point in the geological image data, the method further includes:
and the wavelet denoising module is used for performing wavelet denoising of at least two levels on the geological image data.
Optionally, the wavelet denoising module includes:
the wavelet decomposition sub-module is used for performing wavelet decomposition on the address image data to obtain a first candidate signal and a second candidate signal of at least two levels, wherein the frequency of the second candidate signal is higher than that of the first candidate signal;
a third candidate signal obtaining sub-module, configured to remove, for the second candidate signal of each level, a signal smaller than a threshold to obtain a third candidate signal, where the threshold corresponds to each level;
and the geological image data reconstruction sub-module is used for reconstructing the first candidate signal and the third candidate signal with the highest level into the geological image data.
Optionally, the inflection point detecting module includes:
the gray scale change information calculation submodule is used for translating the geological image data by a preset window and calculating gray scale change information by the following formula:
Figure BDA0002925083210000041
wherein E (u, v) is gray scale change information, w (x, y) is a window, I (x + u, y + v) is gray scale information after the window is translated, and I (x, y) is gray scale information of the geological image data;
and the inflection point determining submodule is used for determining an inflection point based on the gray scale change information.
Optionally, the image processing comprises at least one of:
gain processing, IIR filtering, FIR filtering, arithmetic operation, deconvolution, offset processing, static correction.
Optionally, the semantic layer generating module includes:
an index creating sub-module for creating an index in the cable image data;
the cable information extraction sub-module is used for inquiring the index value meeting the inquiry condition according to the index and extracting the cable information of the cable;
the cable information writing submodule is used for writing the cable information into the three-dimensional body data of the region;
and the semantic layer generating submodule is used for generating a semantic layer according to the cable information in the three-dimensional body data.
In a third aspect, an embodiment of the present invention further provides a computer device for implementing cable display of a three-dimensional map, where the computer device includes:
one or more processors;
a memory for storing one or more programs,
when executed by the one or more processors, cause the one or more processors to implement the cable display method of a three-dimensional map as described in any one of the first aspects.
In a fourth aspect, the embodiment of the present invention further provides a computer-readable storage medium for implementing cable display of a three-dimensional map, where the computer-readable storage medium stores thereon a computer program, and the computer program, when executed by a processor, implements the cable display method of the three-dimensional map according to any one of the first aspects.
The embodiment of the invention provides a cable display method, a cable display device, equipment and a storage medium of a three-dimensional map, wherein the method comprises the following steps: loading three-dimensional body data of a specified area; drawing geological image data for the underground space of the area according to radar data generated when the ground penetrating radar detects the cable in the area; detecting an inflection point in the geological image data; sequentially connecting inflection points to the cable image data to represent the cable; writing cable information of the cable into three-dimensional body data of the region by referring to the cable image data to generate a semantic layer; and loading the three-dimensional body data to display the semantic layer of the cable in a three-dimensional map of the region. Compared with the traditional manual marking and manual statistics, the embodiment of the invention provides a more accurate and visual cable display method, and reduces the cable identification error.
Drawings
Fig. 1 is a flowchart of a cable display method for a three-dimensional map according to an embodiment of the present invention;
FIG. 2A is a schematic representation of a hyperbola of the subsurface provided by one embodiment of the present invention;
FIG. 2B is a schematic diagram of inflection point detection according to a first embodiment of the present invention;
FIG. 2C is a schematic diagram of determining an inflection point according to an embodiment of the present invention;
fig. 2D is a schematic diagram illustrating cable directions displayed in a three-dimensional map according to an embodiment of the present invention;
fig. 3 is a schematic structural diagram of a cable display device for a three-dimensional map according to a second embodiment of the present invention;
fig. 4 is a schematic structural diagram of a computer device according to a third embodiment of the present invention.
Detailed Description
The present invention will be described in further detail with reference to the accompanying drawings and examples. It is to be understood that the specific embodiments described herein are merely illustrative of the invention and are not limiting of the invention. It should be further noted that, for the convenience of description, only some of the structures related to the present invention are shown in the drawings, not all of the structures. Before discussing exemplary embodiments in more detail, it should be noted that some exemplary embodiments are described as processes or methods depicted as flowcharts. Although a flowchart may describe the operations (or steps) as a sequential process, many of the operations can be performed in parallel, concurrently or simultaneously. In addition, the order of the operations may be re-arranged. The process may be terminated when its operations are completed, but may have additional steps not included in the figure. The processes may correspond to methods, functions, procedures, subroutines, and the like.
Example one
Fig. 1 is a flowchart of a cable display method of a three-dimensional map according to an embodiment of the present invention, where the present embodiment is applicable to finding a cable location of an underground cable, and the method may be executed by a cable display apparatus of the three-dimensional map, where the cable display apparatus of the three-dimensional map may be implemented by software and/or hardware, and may be configured in a computer device, for example, a personal computer, a server, a workstation, an application program, and the like, and specifically includes the following steps:
s101, loading three-dimensional body data of a specified area.
In the embodiment of the invention, the WebGIS (Web geographic information System) is adopted for displaying the three-dimensional map of the underground cable, so that a user can conveniently view the position image of the underground cable at a client terminal. The WebGIS (network geographic information system) is a GIS working on a Web network, can realize GIS basic functions of retrieval, query, drawing output, editing and the like of spatial data, and is also a basis of geographic information publishing, sharing and communication cooperation on the Internet, and the WebGIS utilizes Web technology to expand and perfect the geographic information system.
Further, in this embodiment, the WebGIS is a web-based client/server system; the information exchange between the client and the server is performed by using the internet, and when a user uploads underground cable data to the server terminal, other users can obtain the underground cable data distributed on different sites and different computer platforms from the server terminal.
In order to display the position of the underground cable in the WebGIS, the three-dimensional shape data of the area where the underground cable is located is loaded in the platform, the three-dimensional shape data in the embodiment is acquired from a GIS (Geographic Information System), and the acquired three-dimensional shape data is stored in the server side, so that the client side can conveniently view and read the data.
In this embodiment, a ground penetrating radar detection technique is used to locate and identify underground cables. The ground penetrating radar judges the underground target body on the basis of the characteristic analysis of the reflected wave, and judges the spatial position, the structure and the distribution of the underground target body according to the travel time, the amplitude and the wave form of the reflected wave. The electrical difference between the target and the surrounding medium is a basic condition for the ground penetrating radar detection.
And S102, drawing geological image data of the underground space of the region according to radar data generated when the ground penetrating radar detects the cable in the region.
In the embodiment, geological image data is drawn on the underground space of the region according to radar data generated when the ground penetrating radar detects underground cables in the region. Because the energy of the electromagnetic wave detected by the radar can be absorbed by the medium in the process of propagating in the medium, the energy of the electromagnetic wave is weakened along with the increase of the depth, and the signal amplitude is correspondingly reduced, so that the signal identification and recognition are not facilitated. In order to improve the visibility of the acquired signal during the data acquisition process and facilitate the later processing and identification of the image, a variable gain (gain) function is adopted to improve the amplitude of the signal, so that the slight change of the signal is easier to display and identify. Meanwhile, the variable gain function can adjust the size according to the feedback of the collected signal after setting an initial value, and the clipping phenomenon caused by overlarge gain is avoided.
In the method for drawing the geological image data, S102 includes the steps of:
and S1021, acquiring radar data obtained by detecting the ground penetrating radar to an area with a pre-buried cable.
In this embodiment, a ground penetrating radar detection technique is used to locate and identify underground cables. The ground penetrating radar judges the underground target body on the basis of the characteristic analysis of the reflected wave, and judges the spatial position, the structure and the distribution of the underground target body according to the travel time, the amplitude and the wave form of the reflected wave.
The ground penetrating radar mainly comprises a host (main control unit), a transmitter, a transmitting antenna, a receiver and a receiving antenna. Others are also possible including positioning devices such as GPS, odometers or Markers (MARKs), power supplies, and carts, among others. The transmitting and receiving antennas are present in pairs for transmitting and receiving radar waves from underground reflections into the ground. The host is an acquisition system for sending transmit and receive control commands (including parameters such as start and stop time, transmit frequency, number of repetitions, etc.) to the transmitter. The transmitter transmits radar waves into the ground according to the host command, and the receiver starts data acquisition according to the control command. After sampling and A/D conversion, the received reflected signal is converted into digital signal for display and storage.
And S1022, adjusting the radar data into two-dimensional geological image data.
The electromagnetic wave of the ground penetrating radar can determine the depth of the underground cable, namely the cable is positioned. With the known ground dielectric constant, the propagation speed of the electromagnetic pulse of the ground penetrating radar in the medium can be calculated according to the following formula:
Figure BDA0002925083210000081
where c is the propagation velocity of electromagnetic waves in air and epsilon is the dielectric constant of the medium.
Further, the depth s of the target cable is calculated according to the propagation speed c and the propagation time t.
s=vt/2
Alternatively, if the dielectric constant is unknown, a target volume can be determined by drilling, then measuring its depth, finding the propagation velocity over time, and then back-calculating the dielectric constant according to the above equation. After the dielectric constant is obtained, the underground cable can be positioned by using the ground penetrating radar.
After the data acquisition is completed by adopting the ground penetrating radar, the radar data is adjusted into two-dimensional geological image data according to the position of the underground cable.
And S1023, determining a zero line in the geological image data, wherein the zero line is the position of the connection between the cable placed on the ground and the survey line.
Further, a zero line is determined in two-dimensional geological image data, wherein the zero line is the position of the connection of a cable placed on the ground and a measuring line, the cable is placed on the ground to be orthogonal to the measuring line during detection of the zero line, the position of the cable is recorded on a section when an antenna passes through the cable, and the position of the zero line can be determined by identifying the cable.
And S1024, converting the geological image data from a time domain to a frequency spectrum domain.
After the zero line is determined, converting the geological image data from a time domain into a frequency spectrum domain, wherein the time domain represents the change of a signal of the geological image data along with time; the spectrum domain represents a coordinate system used for the frequency-aspect characteristic of the signal of the geological image data, wherein a sine wave is the only existing waveform in the frequency domain, and the time-domain characteristic can be considered as the change rule of the signal intensity of the geological image data along with time, and the frequency-domain characteristic is synthesized by the signal of which single frequency is the signal of the geological image data.
In the process of converting time domain into frequency spectrum domain for spectrum analysis, geological image data is converted from time domain into frequency domain, representing the amplitude distribution of various harmonic frequencies, different detection media have different amplitude spectrum characteristics, and the media can be distinguished by the frequency domain characteristics.
Since the time domain analysis and the frequency domain analysis are two observation planes for the signal. The time domain analysis is to represent the relation of dynamic signals by taking a time axis as a coordinate; the frequency domain analysis is to transform the signal into a coordinate representation with the frequency axis. Generally, the representation of the time domain is more visual and visual, the frequency domain analysis is more concise, and the analysis problem is more profound and convenient, but they are interrelated, none of which is indispensable, and complement each other.
Specifically, Fourier spectrum transformation is adopted in the frequency spectrum analysis, the frequency spectrum is shifted to the circle center through the Fourier spectrum analysis, the image frequency distribution can be clearly seen, and interference signals with periodic regularity can be separated. For example, if there is a sinusoidal interference, it can be seen on the spectrogram shifted to the origin that there exists a symmetrically distributed set of bright spots centered around a certain point in addition to the center, and this set is caused by the interference noise, and it is very intuitive to eliminate the interference by placing a band-stop filter at this position.
The fourier transform equation is as follows:
Figure BDA0002925083210000101
wherein, F (w) is the image function called f (t), f (t) is the image primitive function called F (w), F (w) is the image of f (t), and f (t) is the image of F (w).
Furthermore, the frequency spectrum analysis also adopts Hilbert transform, which is a technology for directly converting the recorded information into instantaneous amplitude, instantaneous phase and instantaneous frequency in a time domain. The instantaneous amplitude is a measure of the reflection intensity, which is proportional to the square root of the total energy of the geological radar signal at that time, and this feature is used to facilitate the determination of changes in the particular formation. When obvious medium layering, a slip zone or a groundwater interface exist in a stratum, the instantaneous amplitude can be changed strongly and is reflected in an instantaneous amplitude profile, namely, the interface position has obvious amplitude change. Instantaneous phase is a measure of the continuity of the in-phase axis across the geological radar cross-section. When an electromagnetic wave propagates in an isotropic homogeneous medium, its phase is continuous; when an electromagnetic wave propagates in a medium in which an anomaly exists, the phase of the electromagnetic wave changes significantly at the anomaly position, and is obviously discontinuous in a cross-sectional view. Therefore, the instantaneous phase can be used for better distinguishing the underground layering from the underground abnormity. When phase discontinuity occurs in the instantaneous phase image section, the layering or abnormity can be judged to exist at the position. The instantaneous frequency is the time change rate of the phase, reflects the lithological change of the formation, is helpful for identifying the formation, and when the electromagnetic wave passes through different medium interfaces, the frequency of the electromagnetic wave is obviously changed, and the change can be more clearly displayed in the instantaneous frequency image section. The stability and lithology changes of the underground medium can be judged by using the magnitude and the stable condition of the instantaneous frequency.
For the same detection object, the three kinds of instant information obviously change at the same position, and the physical property change of the detection object at the position can be reflected. Because the resolution of the instantaneous phase spectrum is the highest among the three parameters, and the change of the instantaneous frequency spectrum and the instantaneous amplitude spectrum is more intuitive, the approximate position of the underground cable or the abnormity can be determined according to the instantaneous frequency spectrum and the instantaneous amplitude spectrum, and then the contour line of the cable and the position of the abnormity can be accurately determined by using the instantaneous phase spectrum.
The problem that due to the nonuniformity of the earth medium, reflected waves, refracted waves, diffracted waves and scattered waves generated in the underground propagation process of high-frequency pulse electromagnetic waves transmitted by a geological radar are superposed with each other, so that the great difficulty of data processing can be solved through Hilbert transformation. Meanwhile, various interference noises caused by the fact that the geological radar records the characteristics of reflected waves by utilizing a wide frequency band can be reduced.
And S1025, executing preset image processing on the geological image data.
In this embodiment, after the conversion of the geological image data from the time domain to the frequency spectrum domain is completed, a preset image process is performed on the geological image data. The image processing comprises IIR filtering, FIR filtering, arithmetic operation, deconvolution, offset processing and static correction. Specifically, an IIR horizontal high-pass filtering and FIR background denoising method is adopted to remove horizontal band interference of low frequency; removing high-frequency interference by adopting IIR horizontal low-pass filtering, FIR horizontal superposition and moving average filtering; removing multiple reflections and adopting prediction deconvolution; removing diffraction and correcting the layer with larger inclination angle and adopting offset treatment; the compensation phase change employs static correction.
Optionally, according to the evaluation processing effect, if the effect is good, the result explanation and the report compiling are performed, and if the effect is not good, the data analysis and processing are performed again until a better effect is achieved.
And S103, detecting an inflection point in the geological image data.
Furthermore, an inflection point of the cable is detected in the geological image data after the image processing, and in order to more conveniently detect the inflection point in the geological image data, at least two levels of wavelet denoising are carried out on the geological image data before the inflection point is detected in the geological image data, so that the denoising processing is carried out. The specific mode is as follows:
as shown in fig. 2A, at least two levels of wavelet denoising are performed on the geological image data to obtain geological image data with higher recognition degree. The wavelet transform mainly utilizes the characteristics of multiple resolutions, decorrelation and base selection flexibility, so that the wavelet transform is very useful in the aspect of geological image data denoising, and the geological image data is clear. After wavelet transformation, different rules are presented under different resolutions, threshold value thresholds are set, and wavelet coefficients are adjusted, so that the purpose of wavelet denoising can be achieved.
In the implementation process of the embodiment, the wavelet decomposition is performed on the geological image data to obtain at least two levels of a first candidate signal and a second candidate signal, wherein the frequency of the second candidate signal is higher than that of the first candidate signal. Specifically, when the geological image data is subjected to wavelet decomposition, a first candidate signal with a wavelet coefficient smaller than a first threshold and a second candidate signal with a wavelet coefficient larger than the first threshold are obtained according to the first threshold of the geological image data.
Further, removing a noise signal from the second candidate signal of each level to obtain a third candidate signal, wherein each level has a corresponding threshold, and the second candidate signal of each level removes signals smaller than the corresponding threshold to obtain the third candidate signal. Illustratively, taking a two-level wavelet as an example, a second threshold is set in a second candidate signal of the two-level wavelet, wavelet decomposition is performed again according to the second threshold, the second candidate signal smaller than the second threshold is considered as noise, the noise in the second candidate signal is removed, the second candidate signal larger than the second threshold is retained, and a third candidate signal is obtained.
Further, the first candidate signal and the third candidate signal with the highest rank are reconstructed into geological image data. Specifically, after wavelet decomposition is completed, wavelet approximate decomposition is performed, wavelet detail decomposition after denoising is performed, and a signal without noise can be obtained.
Optionally, the method adopted by the wavelet denoising is a Matlab (matrix & laboratory) algorithm.
Of course, the above-mentioned two-level wavelet denoising method is only used as an example, and in the implementation of the embodiment of the present invention, if the denoising effect is not satisfactory, multiple wavelet decompositions may be performed on the geological image data according to the actual situation, that is, N-level wavelet denoising is adopted until the satisfactory denoising effect is achieved.
As shown in fig. 2B, after the wavelet denoising is completed, an inflection point is detected in the geological image data, the cable image data is identified, and the position of the underground cable is determined. The geological image data is expressed as matrixes in a computer, and one value in each matrix corresponds to one pixel point of the geological image data. In the geological image data matrix, traversing the numerical values in the matrix through a moving window, and finding a place with obvious pixel change, namely determining an inflection point. As the data section of the cable pipeline generally presents a hyperbolic shape in the ground penetrating radar detection imaging, and the vertex of the hyperbolic curve just reflects the cable routing and geological image data information required to be extracted, the inflection point detection method is adopted to identify the underground cable radar image.
In this embodiment, the geological image data is translated by a preset window to detect the gray-scale change information, and the calculation of the gray-scale change information is implemented by the following formula:
Figure BDA0002925083210000131
wherein E (u, v) is gray scale change information, w (x, y) is a window, I (x + u, y + v) is gray scale information after the window is translated, and I (x, y) is gray scale information of geological image data;
further, determining an inflection point based on the gray-scale variation information includes:
and carrying out Taylor series simplification on the gray change information:
Figure BDA0002925083210000132
the local minute movement amount [ u, v ] and the gradation change information are equivalently operated to obtain:
Figure BDA0002925083210000133
wherein [ u, v ] is the displacement of the window, and M is a 2x2 partial derivative matrix;
Figure BDA0002925083210000134
the geological image data change caused by window movement is the characteristic value analysis of the real symmetric matrix M, an angular point response function R is set, and whether a pixel is an angular point is judged by judging the size of R:
R=detM-k(traceM)2
where detM is a determinant of M matrix, and traceM is a characteristic value that a trajectory operation R on the M matrix depends on M.
And if the numerical value of the response function is larger than a preset threshold value, determining that an inflection point exists in the window. Specifically, in geological image data, | R | is large for the corner points; for flat areas, | R | is small; for edges, R is negative. From this, the value of R in the window can be calculated to determine where the inflection point exists.
And S104, sequentially connecting inflection points to the cable image data to represent the cable.
As shown in fig. 2C, after the detection of the inflection point is completed, the inflection points are extracted from the cable image data and are sequentially connected to represent the cable. The cable is positioned at the inflection point, so that the cable can be positioned by connecting the inflection point.
And S105, writing the cable information of the cable into the three-dimensional body data of the region by referring to the cable image data to generate a semantic layer.
And writing the cable information of the cable into the three-dimensional body data of the region according to the cable image data obtained by the connection inflection point to generate a semantic layer. The cable information of the cable comprises pipeline attribute data of pipeline characteristics and properties and spatial data taking spatial positions as parameters, the pipeline attribute data comprises channel names, communication relations, pipe hole numbers, cable lengths, states, trends and the like, and the spatial data comprises cable (starting point, inflection point and end point) coordinates, cable middle joint coordinates, installation positions, upper jacking buried depth (or elevation) and the like.
Specifically, the three-dimensional body data is stored in cable image data of a three-dimensional body surface, an index is created in the cable image data, an index value meeting a query condition is inquired according to the index, cable information of the cable is extracted, the cable information is written into the three-dimensional body data of the area, and a three-dimensional semantic layer is generated according to the cable information in the three-dimensional body data.
And S106, loading the three-dimensional body data to display a semantic layer of the cable in the three-dimensional map of the area.
As shown in fig. 2D, the three-dimensional semantic layer is loaded with three-dimensional shape data to display the semantic layer of the cable in the three-dimensional map of the area. Specifically, a coordinate system is established in the three-dimensional semantic layer, and three-dimensional body data corresponding to the cable image data is loaded in a Web page so as to display the semantic layer of the cable in a three-dimensional map of the area and express the three-dimensional trend of the cable.
It should be noted that, the WebGIS used for drawing the three-dimensional map and displaying the cable trend in the embodiments of the present invention mainly includes the steps of implementing 3D digitization of the underground cable by using a3D modeling technology, a computer technology, a network technology, a database technology, etc., developing the underground cable by using an international Web3D standard VRML language and a network cross-platform language Java3D according to the principles of 3D digitization, networking, reality, simplicity, visualization, and high performance efficiency, and providing three-dimensional visualization information of the underground cable to implement view management of the underground cable in different directions.
The embodiment of the invention provides a cable display method, a cable display device, equipment and a storage medium of a three-dimensional map, wherein the method comprises the following steps: loading three-dimensional body data of a specified area; drawing geological image data for the underground space of the area according to radar data generated when the ground penetrating radar detects the cable in the area; detecting an inflection point in the geological image data; sequentially connecting inflection points to the cable image data to represent the cable; writing cable information of the cable into three-dimensional body data of the region by referring to the cable image data to generate a semantic layer; and loading the three-dimensional body data to display the semantic layer of the cable in a three-dimensional map of the region. Compared with the traditional manual marking and manual statistics, the embodiment of the invention provides a more accurate and visual cable display method, and reduces the cable identification error.
It should be noted that, for simplicity of description, the method embodiments are described as a series of acts or combination of acts, but those skilled in the art will recognize that the present invention is not limited by the illustrated order of acts, as some steps may occur in other orders or concurrently in accordance with the embodiments of the present invention. Further, those skilled in the art will appreciate that the embodiments described in the specification are presently preferred and that no particular act is required to implement the invention.
Example two
Fig. 3 is a block diagram of a cable display device of a three-dimensional map according to a second embodiment of the present invention, which may specifically include the following modules:
a three-dimensional body data loading module 301, configured to load three-dimensional body data of a specified area;
the geological image data drawing module 302 is used for drawing geological image data for the underground space of the region according to radar data generated when the ground penetrating radar detects cables in the region;
an inflection point detection module 303, configured to detect an inflection point in the geological image data;
a cable generation module 304, configured to sequentially connect the inflection points in the cable image data to represent the cable;
a semantic layer generating module 305, configured to write the cable information of the cable into the three-dimensional body data of the area by referring to the cable image data, so as to generate a semantic layer;
a cable display module 306, configured to load the three-dimensional shape data, so as to display the semantic map layer of the cable in a three-dimensional map of the area.
Optionally, the geological image data drawing module includes:
the radar data acquisition submodule is used for acquiring radar data obtained by detecting the ground penetrating radar to the area pre-buried with the cable;
the two-dimensional address image data acquisition submodule is used for adjusting the radar data into two-dimensional geological image data;
the zero line submodule is used for determining a zero line in the geological image data, and the zero line is the position of the connection between a cable placed on the ground and a measuring line;
the frequency spectrum domain conversion sub-module is used for converting the geological image data from a time domain into a frequency spectrum domain;
and the image processing submodule is used for executing preset image processing on the geological image data.
Optionally, before detecting an inflection point in the geological image data, the method further includes:
and the wavelet denoising module is used for performing wavelet denoising of at least two levels on the geological image data.
Optionally, the wavelet denoising module includes:
the wavelet decomposition sub-module is used for performing wavelet decomposition on the address image data to obtain a first candidate signal and a second candidate signal of at least two levels, wherein the frequency of the second candidate signal is higher than that of the first candidate signal;
a third candidate signal obtaining sub-module, configured to remove, for the second candidate signal of each level, a signal smaller than a threshold to obtain a third candidate signal, where the threshold corresponds to each level;
and the geological image data reconstruction sub-module is used for reconstructing the first candidate signal and the third candidate signal with the highest level into the geological image data.
Optionally, the inflection point detecting module includes:
the gray scale change information calculation submodule is used for translating the geological image data by a preset window and calculating gray scale change information by the following formula:
Figure BDA0002925083210000171
wherein E (u, v) is gray scale change information, w (x, y) is a window, I (x + u, y + v) is gray scale information after the window is translated, and I (x, y) is gray scale information of the geological image data;
and the inflection point determining submodule is used for determining an inflection point based on the gray scale change information.
Optionally, the image processing comprises at least one of:
gain processing, IIR filtering, FIR filtering, arithmetic operation, deconvolution, offset processing, static correction.
Optionally, the semantic layer generating module includes:
an index creating sub-module for creating an index in the cable image data;
the cable information extraction sub-module is used for inquiring the index value meeting the inquiry condition according to the index and extracting the cable information of the cable;
the cable information writing submodule is used for writing the cable information into the three-dimensional body data of the region;
and the semantic layer generating submodule is used for generating a semantic layer according to the cable information in the three-dimensional body data.
The cable display device of the three-dimensional map provided by the embodiment of the invention can execute the cable display method of the three-dimensional map provided by any embodiment of the invention, and has corresponding functional modules and beneficial effects of the execution method.
EXAMPLE III
Fig. 4 is a schematic structural diagram of a computer device according to a third embodiment of the present invention. FIG. 4 illustrates a block diagram of an exemplary computer device 12 suitable for use in implementing embodiments of the present invention. The computer device 12 shown in FIG. 4 is only one example and should not bring any limitations to the functionality or scope of use of embodiments of the present invention.
As shown in FIG. 4, computer device 12 is in the form of a general purpose computing device. The components of computer device 12 may include, but are not limited to: one or more processors or processing units 16, a system memory 28, and a bus 18 that couples various system components including the system memory 28 and the processing unit 16.
Bus 18 represents one or more of any of several types of bus structures, including a memory bus or memory controller, a peripheral bus, an accelerated graphics port, and a processor or local bus using any of a variety of bus architectures. By way of example, such architectures include, but are not limited to, Industry Standard Architecture (ISA) bus, micro-channel architecture (MAC) bus, enhanced ISA bus, Video Electronics Standards Association (VESA) local bus, and Peripheral Component Interconnect (PCI) bus.
Computer device 12 typically includes a variety of computer system readable media. Such media may be any available media that is accessible by computer device 12 and includes both volatile and nonvolatile media, removable and non-removable media.
The system memory 28 may include computer system readable media in the form of volatile memory, such as Random Access Memory (RAM)30 and/or cache memory 32. Computer device 12 may further include other removable/non-removable, volatile/nonvolatile computer system storage media. By way of example only, storage system 34 may be used to read from and write to non-removable, nonvolatile magnetic media (not shown in FIG. 4, and commonly referred to as a "hard drive"). Although not shown in FIG. 4, a magnetic disk drive for reading from and writing to a removable, nonvolatile magnetic disk (e.g., a "floppy disk") and an optical disk drive for reading from or writing to a removable, nonvolatile optical disk (e.g., a CD-ROM, DVD-ROM, or other optical media) may be provided. In these cases, each drive may be connected to bus 18 by one or more data media interfaces. Memory 28 may include at least one program product having a set (e.g., at least one) of program modules that are configured to carry out the functions of embodiments of the invention.
A program/utility 40 having a set (at least one) of program modules 42 may be stored, for example, in memory 28, such program modules 42 including, but not limited to, an operating system, one or more application programs, other program modules, and program data, each of which examples or some combination thereof may comprise an implementation of a network environment. Program modules 42 generally carry out the functions and/or methodologies of the described embodiments of the invention.
Computer device 12 may also communicate with one or more external devices 14 (e.g., keyboard, pointing device, display 24, etc.), with one or more devices that enable a user to interact with computer device 12, and/or with any devices (e.g., network card, modem, etc.) that enable computer device 12 to communicate with one or more other computing devices. Such communication may be through an input/output (I/O) interface 22. Also, computer device 12 may communicate with one or more networks (e.g., a Local Area Network (LAN), a Wide Area Network (WAN), and/or a public network such as the Internet) via network adapter 20. As shown, network adapter 20 communicates with the other modules of computer device 12 via bus 18. It should be understood that although not shown in the figures, other hardware and/or software modules may be used in conjunction with computer device 12, including but not limited to: microcode, device drivers, redundant processing units, external disk drive arrays, RAID systems, tape drives, and data backup storage systems, among others.
The processing unit 16 executes various functional applications and data processing, such as a cable display method of a three-dimensional map provided by an embodiment of the present invention, by executing a program stored in the system memory 28.
Example four
The fourth embodiment of the present invention further provides a computer-readable storage medium, where a computer program is stored on the computer-readable storage medium, and when the computer program is executed by a processor, the computer program implements each process of the cable display method for a three-dimensional map, and can achieve the same technical effect, and in order to avoid repetition, the details are not repeated here.
A computer readable storage medium may include, for example, but is not limited to, an electronic, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any combination of the foregoing. More specific examples (a non-exhaustive list) of the computer readable storage medium would include the following: an electrical connection having one or more wires, a portable computer diskette, a hard disk, a Random Access Memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disc read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing. In the context of this document, a computer readable storage medium may be any tangible medium that can contain, or store a program for use by or in connection with an instruction execution system, apparatus, or device.
It is to be noted that the foregoing is only illustrative of the preferred embodiments of the present invention and the technical principles employed. It will be understood by those skilled in the art that the present invention is not limited to the particular embodiments described herein, but is capable of various obvious changes, rearrangements and substitutions as will now become apparent to those skilled in the art without departing from the scope of the invention. Therefore, although the present invention has been described in greater detail by the above embodiments, the present invention is not limited to the above embodiments, and may include other equivalent embodiments without departing from the spirit of the present invention, and the scope of the present invention is determined by the scope of the appended claims.

Claims (10)

1. A cable display method of a three-dimensional map is characterized by comprising the following steps:
loading three-dimensional body data of a specified area;
drawing geological image data for the underground space of the region according to radar data generated when the ground penetrating radar detects the cable in the region;
detecting an inflection point in the geological image data;
sequentially connecting the inflection points to the cable image data to represent the cable;
writing the cable information of the cable into the three-dimensional body data of the region by referring to the cable image data to generate a semantic layer;
and loading the three-dimensional body data to display the semantic layer of the cable in a three-dimensional map of the region.
2. The method of claim 1, wherein mapping geological image data of the subsurface of the region from radar data generated by a ground penetrating radar during cable detection in the region comprises:
acquiring radar data obtained by detecting the area pre-buried with the cable by the ground penetrating radar;
adjusting the radar data into two-dimensional geological image data;
determining a zero line in the geological image data, wherein the zero line is the joint position of a cable placed on the ground and a measuring line;
converting the geological image data from a time domain to a frequency spectrum domain;
and executing preset image processing on the geological image data.
3. The method of claim 1, prior to detecting an inflection point in the geological image data, further comprising:
and performing wavelet denoising of at least two levels on the geological image data.
4. The method of claim 3, wherein performing at least two levels of wavelet de-noising on the geological image data comprises:
performing wavelet decomposition on the address image data to obtain a first candidate signal and a second candidate signal of at least two levels, wherein the frequency of the second candidate signal is higher than that of the first candidate signal;
removing signals smaller than a threshold value from the second candidate signal of each level to obtain a third candidate signal, wherein the threshold value corresponds to each level;
reconstructing the first candidate signal and the third candidate signal of the highest rank as the geological image data.
5. The method of claim 1, wherein detecting an inflection point in the geological image data comprises:
translating the geological image data by a preset window, and calculating gray change information by the following formula:
Figure FDA0002925083200000021
wherein E (u, v) is gray scale change information, w (x, y) is a window, I (x + u, y + v) is gray scale information after the window is translated, and I (x, y) is gray scale information of the geological image data;
an inflection point is determined based on the gray-scale variation information.
6. The method of claim 2, wherein the image processing comprises at least one of:
gain processing, IIR filtering, FIR filtering, arithmetic operation, deconvolution, offset processing, static correction.
7. The method according to claims 1-3, wherein writing cable information of the cable into the three-dimensional shape data of the region with reference to the cable image data to generate a semantic layer comprises:
creating an index in the cable image data;
inquiring an index value meeting the inquiry condition according to the index, and extracting the cable information of the cable;
writing the cable information into the three-dimensional body data of the region;
and generating a semantic layer according to the cable information in the three-dimensional body data.
8. A cable display apparatus for a three-dimensional map, comprising:
the three-dimensional body data loading module is used for loading three-dimensional body data of a specified area;
the geological image data drawing module is used for drawing geological image data for the underground space of the region according to radar data generated when the ground penetrating radar detects the cable in the region;
the inflection point detection module is used for detecting an inflection point in the geological image data;
a cable generation module for sequentially connecting the inflection points to the cable image data to represent the cable;
the semantic layer generating module is used for writing the cable information of the cable into the three-dimensional body data of the region by referring to the cable image data so as to generate a semantic layer;
and the cable display module is used for loading the three-dimensional body data so as to display the semantic map layer of the cable in a three-dimensional map of the region.
9. A computer device for enabling cable display of a three-dimensional map, the computer device comprising:
one or more processors;
a memory for storing one or more programs,
when executed by the one or more processors, cause the one or more processors to implement the cable display method of a three-dimensional map of any of claims 1-7.
10. A computer-readable storage medium for enabling cable display of a three-dimensional map, characterized in that the computer-readable storage medium has stored thereon a computer program which, when executed by a processor, implements a cable display method of a three-dimensional map according to any one of claims 1-7.
CN202110129852.0A 2021-01-29 2021-01-29 Cable display method, device, equipment and storage medium for three-dimensional map Pending CN112860832A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110129852.0A CN112860832A (en) 2021-01-29 2021-01-29 Cable display method, device, equipment and storage medium for three-dimensional map

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110129852.0A CN112860832A (en) 2021-01-29 2021-01-29 Cable display method, device, equipment and storage medium for three-dimensional map

Publications (1)

Publication Number Publication Date
CN112860832A true CN112860832A (en) 2021-05-28

Family

ID=75987090

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110129852.0A Pending CN112860832A (en) 2021-01-29 2021-01-29 Cable display method, device, equipment and storage medium for three-dimensional map

Country Status (1)

Country Link
CN (1) CN112860832A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115830263A (en) * 2023-02-16 2023-03-21 四川省地质调查研究院 Hydrogeology dynamic three-dimensional model construction method and system

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111429575A (en) * 2020-04-01 2020-07-17 中冶建筑研究总院(深圳)有限公司 Three-dimensional visual monitoring method, system, equipment and storage medium
CN111914627A (en) * 2020-06-18 2020-11-10 广州杰赛科技股份有限公司 Vehicle identification and tracking method and device

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111429575A (en) * 2020-04-01 2020-07-17 中冶建筑研究总院(深圳)有限公司 Three-dimensional visual monitoring method, system, equipment and storage medium
CN111914627A (en) * 2020-06-18 2020-11-10 广州杰赛科技股份有限公司 Vehicle identification and tracking method and device

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
何旭涛等: "海底电缆工程多源勘察数据管理及三维可视化研究", 南方能源建设, vol. 7, pages 2 - 5 *
涂碧君: "基于探地雷达系统的地下管线自动识别技术", 中国优秀博硕士学位论文全文数据库(硕士)信息科技辑, no. 01, pages 10 - 13 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115830263A (en) * 2023-02-16 2023-03-21 四川省地质调查研究院 Hydrogeology dynamic three-dimensional model construction method and system

Similar Documents

Publication Publication Date Title
US7881501B2 (en) Local dominant wave-vector analysis of seismic data
US8094515B2 (en) Seismic data visualizations
CN111164462B (en) Artificial source surface wave exploration method, surface wave exploration device and terminal equipment
US20090132170A1 (en) Seismic data processing and visualization
AU2019382288B2 (en) Passive seismic imaging
CA2818790C (en) Seismic trace attribute
CN114428324B (en) Pre-stack high-angle fast Fourier transform seismic imaging method, system and equipment
US5673191A (en) Method and apparatus for identifying geological structures using wavelet analysis of potential fields
US6490526B2 (en) Method for characterization of multi-scale geometric attributes
CN111337992A (en) Method for obtaining depth of field source based on downward continuation of bit field data
CN112862888A (en) Cable positioning method and device, computer equipment and storage medium
CN115100363B (en) Underground abnormal body three-dimensional modeling method and device based on ground penetrating radar
JP2000512385A (en) Sampling and reconstruction of propagating wave fields
CN112860832A (en) Cable display method, device, equipment and storage medium for three-dimensional map
CA2693247C (en) Seismic attributes for structural analysis
AU2012397811A1 (en) Measurement correction apparatus, methods, and systems
CN115859081A (en) Visual detection method and device for pipeline
Sun et al. Research on detection and visualization of underground pipelines
CN114167511A (en) Continuous-fraction expansion downward continuation-based bit field data rapid inversion method
CN111582114A (en) Seismic fault identification method, device, equipment and storage medium
CN111880223B (en) Attribute analysis method for well-side sound wave remote detection imaging fracture-cave carving
NO20181090A1 (en) Visualizing attributes of multiple fault surfaces in real time
CN115586527B (en) Cloud road deformation early warning system based on DS-InSAR technology
CA2796915C (en) Systems and methods for computing a default 3d variogram model
CN115128675A (en) Method, device, equipment and storage medium for suppressing interbed multiples

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination