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

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

Info

Publication number
CN112860832B
CN112860832B CN202110129852.0A CN202110129852A CN112860832B CN 112860832 B CN112860832 B CN 112860832B CN 202110129852 A CN202110129852 A CN 202110129852A CN 112860832 B CN112860832 B CN 112860832B
Authority
CN
China
Prior art keywords
cable
image data
data
module
information
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.)
Active
Application number
CN202110129852.0A
Other languages
Chinese (zh)
Other versions
CN112860832A (en
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/CN112860832B/en
Publication of CN112860832A publication Critical patent/CN112860832A/en
Application granted granted Critical
Publication of CN112860832B publication Critical patent/CN112860832B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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

Landscapes

  • Engineering & Computer Science (AREA)
  • Databases & Information Systems (AREA)
  • Theoretical Computer Science (AREA)
  • Remote Sensing (AREA)
  • Data Mining & Analysis (AREA)
  • Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Radar Systems Or Details Thereof (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

The embodiment of the invention provides a cable display method, device and equipment for a three-dimensional map and a storage medium, wherein the method comprises the following steps: loading three-dimensional shape data of a designated 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 inflection points in the geologic image data; sequentially connecting inflection points to the cable image data to represent the cable; writing cable information of the cable into three-dimensional shape data of the region with reference to the cable image data to generate a semantic layer; the three-dimensional shape data is loaded 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 error of cable identification.

Description

Cable display method, device and equipment for three-dimensional map and storage medium
Technical Field
The embodiment of the invention relates to the technical field of power engineering, in particular to a cable display method, device and equipment for a three-dimensional map and a storage medium.
Background
Nowadays, with the rapid development of urban construction, after entering modern society, underground cable transmission modes are commonly adopted in large cities due to the reasons of shortage of urban land, high traffic pressure, urban construction and the like. Especially in coastal areas, typhoons are frequent, overhead lines are very easy to be injured, and the problem is solved perfectly by burying cables. Compared with overhead lines, the cable has the advantages of small occupied area, reliable power transmission, strong anti-interference capability and the like.
In the existing scheme, the identification plates are marked in the cable area by related staff to serve as the positions of the cables mainly by means of manual marking and manual statistics one by one, and the method is easy to damage the identification plates or missing the identification plates by external factors, so that the correctness of the positions of the cables is doubtful finally.
Disclosure of Invention
The embodiment of the invention provides a cable display method, device, equipment and storage medium of a three-dimensional map, which are used for solving 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 shape data of a designated area;
Drawing geological image data for an underground space of the region according to radar data generated when a ground penetrating radar detects a cable in the region;
Detecting inflection points in the geologic image data;
sequentially connecting the inflection points to the cable image data to represent the cable;
Writing cable information of the cable into three-dimensional shape data of the region with reference to the cable image data to generate a semantic layer;
And loading the three-dimensional shape data to display the semantic layer of the cable in a three-dimensional map of the region.
Optionally, the mapping geological image data to the underground space of the region according to the radar data generated when the ground penetrating radar detects the cable in the region includes:
Acquiring radar data obtained by detecting the ground penetrating radar in the region in which the cable is embedded;
adjusting the radar data to two-dimensional geological image data;
determining a zero line in the geological image data, wherein the zero line is a junction position of a cable and a measuring line which are placed on the ground;
converting the geologic image data from a time domain to a frequency spectrum domain;
And executing preset image processing on the geological image data.
Optionally, before detecting the inflection point in the geologic image data, further comprising:
and performing at least two levels of wavelet denoising on the geological image data.
Optionally, performing at least two levels of wavelet denoising on the geologic image data, including:
performing wavelet decomposition on the address image data to obtain at least two levels of first candidate signals and second candidate signals, wherein the frequency of the second candidate signals is higher than that of the first candidate signals;
removing signals smaller than a threshold value from the second candidate signals of each level to obtain third candidate signals, wherein the threshold value corresponds to each level;
Reconstructing the first candidate signal and the third candidate signal with the highest level into the geological image data.
Optionally, detecting an inflection point in the geologic image data includes:
translating the geological image data with a preset window, and calculating gray level change information according to the following formula:
wherein E (u, v) is gray level change information, w (x, y) is a window, I (x+u, y+v) is gray level information after the window is translated, and I (x, y) is gray level information of the geological image data;
And determining an inflection point based on the gradation change information.
Optionally, the image processing includes at least one of:
Gain processing, IIR filtering, FIR filtering, arithmetic operations, deconvolution, offset processing, static correction.
Optionally, writing cable information of the cable into three-dimensional shape data of the region with reference to the cable image data to generate a semantic layer, including:
Creating an index in the cable image data;
According to the index value meeting the query condition, extracting the cable information of the cable;
writing the cable information into three-dimensional shape 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 device for a three-dimensional map, including:
The three-dimensional shape data loading module is used for loading three-dimensional shape data of a designated 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 inflection points in the geological image data;
A cable generating module for sequentially connecting the inflection points to the cable image data to represent the cable;
The semantic layer generation module is used for writing cable information of the cable into the three-dimensional shape 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 shape data so as to display the semantic layer of the cable in the three-dimensional map of the region.
Optionally, the geological image data rendering module includes:
The radar data acquisition sub-module is used for acquiring radar data obtained by detecting the ground penetrating radar in the region embedded with the cable;
The two-dimensional address image data acquisition sub-module is used for adjusting the radar data into two-dimensional geological image data;
the zero line sub-module is used for determining a zero line in the geological image data, wherein the zero line is a position where a cable placed on the ground is intersected with a measuring line;
a spectrum domain conversion sub-module, configured to convert the geological image data from a time domain to a spectrum domain;
and the image processing sub-module is used for executing preset image processing on the geological image data.
Optionally, before detecting the inflection point in the geologic image data, further comprising:
And the wavelet denoising module is used for performing at least two levels of wavelet denoising on the geological image data.
Optionally, the wavelet denoising module includes:
the wavelet decomposition sub-module is used for carrying out wavelet decomposition on the address image data to obtain at least two levels of first candidate signals and second candidate signals, and the frequency of the second candidate signals is higher than that of the first candidate signals;
A third candidate signal obtaining sub-module, configured to remove, for the second candidate signal of each level, a signal smaller than a threshold value, to obtain a third candidate signal, where the threshold value 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 level change information calculation sub-module is used for translating the geological image data with a preset window, and calculating gray level change information through the following formula:
wherein E (u, v) is gray level change information, w (x, y) is a window, I (x+u, y+v) is gray level information after the window is translated, and I (x, y) is gray level information of the geological image data;
And the inflection point determining submodule is used for determining an inflection point based on the gray change information.
Optionally, the image processing includes at least one of:
Gain processing, IIR filtering, FIR filtering, arithmetic operations, deconvolution, offset processing, static correction.
Optionally, the semantic layer generating module includes:
an index creation sub-module for creating an index in the cable image data;
the cable information extraction sub-module is used for inquiring index values meeting inquiry conditions according to the indexes and extracting cable information of the cable;
The cable information writing sub-module is used for writing the cable information into the three-dimensional shape data of the region;
And the semantic layer generation sub-module 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 apparatus for implementing cable display of a three-dimensional map, the computer apparatus including:
One or more processors;
A memory for storing one or more programs,
The 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 in any of the first aspects.
In a fourth aspect, embodiments of the present invention further provide a computer readable storage medium for implementing cable display of a three-dimensional map, the computer readable storage medium having stored thereon a computer program which, when executed by a processor, implements a cable display method of a three-dimensional map as set forth in any one of the first parties.
The embodiment of the invention provides a cable display method, device and equipment for a three-dimensional map and a storage medium, wherein the method comprises the following steps: loading three-dimensional shape data of a designated 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 inflection points in the geologic image data; sequentially connecting inflection points to the cable image data to represent the cable; writing cable information of the cable into three-dimensional shape data of the region with reference to the cable image data to generate a semantic layer; the three-dimensional shape data is loaded 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 error of cable identification.
Drawings
Fig. 1 is a flowchart of a cable display method of a three-dimensional map according to an embodiment of the present invention;
FIG. 2A is a schematic illustration of an underground hyperbola according to an embodiment of the invention;
FIG. 2B is a schematic diagram of inflection point detection according to an embodiment of the present invention;
FIG. 2C is a schematic diagram of a certain inflection point according to an embodiment of the present invention;
fig. 2D is a schematic diagram showing a cable trend 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 of 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 invention is described in further detail below with reference to the drawings and examples. It is to be understood that the specific embodiments described herein are merely illustrative of the invention and are not limiting thereof. It should be further noted that, for convenience of description, only some, but not all of the structures related to the present invention are shown in the drawings. Before discussing exemplary embodiments in more detail, it should be mentioned that some exemplary embodiments are described as processes or methods depicted as flowcharts. Although a flowchart depicts operations (or steps) as a sequential process, many of the operations can be performed in parallel, concurrently, or at the same time. Furthermore, the order of the operations may be rearranged. The process may be terminated when its operations are completed, but may have additional steps not included in the figures. The processes may correspond to methods, functions, procedures, subroutines, and the like.
Example 1
Fig. 1 is a flowchart of a cable display method of a three-dimensional map according to a first embodiment of the present invention, where the method may be applied to find a cable position of an underground cable, and the method may be performed by a cable display device of the three-dimensional map, where the cable display device 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 so on, and specifically includes the following steps:
S101, loading three-dimensional shape data of a designated region.
In the embodiment of the invention, the three-dimensional map for displaying the underground cable adopts WebGIS (network geographic information system), so that a user can conveniently view the position image of the underground cable at the client terminal. The WebGIS (geographic information system) refers to a GIS working on the Web, can realize basic functions of GIS such as space data retrieval, inquiry, drawing output, editing and the like, and is also a basis of geographic information release, sharing and communication collaboration on the internet, and the WebGIS is a technology for expanding and perfecting the geographic information system by using a Web technology.
Further, webGIS in this embodiment is a network-based client/server system; information exchange between the client and the server is performed using the internet, and when a user uploads underground cable data to the server terminal, other users can acquire 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, three-dimensional shape data of the region where the underground cable is located is loaded in the platform, the three-dimensional shape data in the embodiment is obtained from the GIS (Geographic Information System geographic information system), and the obtained 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 the underground cable. The ground penetrating radar judges the underground target body on the basis of the characteristic analysis of the reflected wave, and judges the space position, structure and distribution of the underground target body according to the travel time, amplitude and wave form of the reflected wave. Wherein, the electrical difference between the target body and the surrounding medium is the basic condition of the ground penetrating radar detection.
S102, 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.
In the present embodiment, geological image data is drawn for an underground space of a region from radar data generated when a ground penetrating radar detects an underground cable 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 electromagnetic wave energy 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 signals during data acquisition, the later processing and recognition of the images are facilitated, and the amplitude of the signals is improved by using a gain function, so that the subtle changes of the signals are easier to display and recognize. Meanwhile, the variable gain function can be adjusted according to feedback of the acquired signal after the initial value is set, and clipping phenomenon caused by overlarge gain is avoided.
In the way of drawing geologic image data, S102 includes the steps of:
s1021, radar data obtained by detecting the ground penetrating radar in the area with the embedded cable is obtained.
In this embodiment, a ground penetrating radar detection technique is used to locate and identify the underground cable. The ground penetrating radar judges the underground target body on the basis of the characteristic analysis of the reflected wave, and judges the space position, structure and distribution of the underground target body according to the travel time, amplitude and 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. Other possibilities include positioning devices such as GPS, odometer or Marker (MARK), power supplies, and carts. The transmitting and receiving antennas are present in pairs for transmitting and receiving radar waves reflected from the subsurface into the subsurface. The host is an acquisition system for sending transmit and receive control commands (including start-stop time, transmit frequency, repetition number, etc.) to the transmitter. The transmitter transmits radar waves into the subsurface according to a host command, and the receiver begins data acquisition according to a control command. The received reflected signal is converted into a digital signal by sampling and a/D conversion, and displayed and stored.
And S1022, adjusting the radar data into two-dimensional geological image data.
The electromagnetic wave energy of the ground penetrating radar determines the depth of the underground cable, i.e. positions the cable. Knowing the surface dielectric constant, the propagation velocity of an electromagnetic pulse of a ground penetrating radar in a medium can be calculated according to the following formula:
where c is the propagation velocity of an electromagnetic wave in air and ε 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 may be determined by drilling, measuring its depth, determining the propagation velocity over time, and then back calculating the dielectric constant according to the equation above. After the dielectric constant is obtained, the ground penetrating radar can be used for positioning the underground cable.
After the data acquisition is completed by adopting the ground penetrating radar, the radar data are adjusted into two-dimensional geological image data according to the position of the underground cable.
S1023, determining a zero line in the geological image data, wherein the zero line is the junction position of a cable and a measuring line which are placed on the ground.
Further, a zero line is determined in the two-dimensional geological image data, wherein the zero line is the junction position of a cable and a measuring line which are placed on the ground, the cable is placed on the ground to be orthogonal to the measuring line when the zero line is detected, the position of the cable can be recorded on the section of the antenna when the antenna passes through the cable, and the position of the zero line can be determined by identifying the cable.
S1024, converting the geological image data from the time domain to the spectrum domain.
Converting from a time domain to a spectral domain in the geologic image data upon completion of the zero line determination, wherein the time domain represents a change in a signal of the geologic image data over time; the spectrum domain represents a coordinate system used when signals of the geological image data are characterized in terms of frequency, wherein a sine wave is the only waveform existing in the frequency domain, and it can be considered that the change rule of the signal intensity of the geological image data with time is time domain characteristic, and the signals of the geological image data are synthesized by signals of single frequencies and are frequency domain characteristic.
In the process of converting the time domain into the frequency domain and carrying out spectrum analysis, the geological image data is converted into the frequency domain from the time domain, the amplitude distribution of various harmonic frequencies is represented, different detection mediums have different amplitude spectrum characteristics, and the mediums can be distinguished by the frequency domain characteristics.
Since time domain analysis and frequency domain analysis are two observation planes for a signal. The time domain analysis is to use a time axis as a coordinate to represent the relation of dynamic signals; the frequency domain analysis is to change the signal to be represented by coordinates on the frequency axis. Generally, the time domain is more vivid and visual, the frequency domain analysis is more concise, the analysis problem is more profound and convenient, however, the time domain and the frequency domain are mutually connected, and the time domain and the frequency domain are inexhaustible and supplement each other.
Specifically, the frequency spectrum analysis adopts Fourier spectrum transformation, and the Fourier spectrum analysis shifts the frequency spectrum to the circle center so that the image frequency distribution can be clearly seen, and interference signals with periodicity regularity can be separated. For example, if a sinusoidal disturbance is included, it can be seen that there is a set of symmetrically distributed bright spots centered at a certain point in addition to the center on the spectrogram shifted to the origin, and this set is generated by the disturbance noise, and this can be intuitively eliminated by placing a band-stop filter at this position.
The fourier transform formula is as follows:
Wherein F (w) is called an image function of F (t), F (t) is called an image primitive function of F (w), F (w) is an image of F (t), and F (t) is an image of F (w).
Further, the spectrum analysis also adopts hilbert transformation, which is a technology of directly converting the information of the record into instantaneous amplitude, instantaneous phase and instantaneous frequency in the time domain. The instantaneous amplitude is a measure of the intensity of the reflection, which is proportional to the square root of the total energy of the geological radar signal at that moment, and this feature is used to determine the change in a particular formation. When the stratum has obvious medium layering, sliding crack zone or subsurface water interface, the instantaneous amplitude will change strongly, and the apparent amplitude change is reflected in the instantaneous amplitude profile, i.e. the interface position. Instantaneous phase is a measure of the continuity of the same phase axis across the geological radar profile. 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 abnormality exists, its phase will vary significantly at the location of the abnormality, being significantly discontinuous in the cross-sectional view. Thus, the instantaneous phase is utilized to better distinguish the underground layering and the underground abnormality. When phase discontinuity occurs in the instantaneous phase image profile, it can be judged that there is layering or abnormality. The instantaneous frequency is the time rate of change of phase, which reflects the lithology change of the formation, and helps to identify the formation, and when the electromagnetic wave passes through different medium interfaces, the electromagnetic wave frequency will change obviously, and the change can be displayed more clearly in the instantaneous frequency image section. The magnitude and stability of the instantaneous frequency can be used to determine the stability and lithology changes of the subsurface medium.
For the same detected object, the three instantaneous messages can reflect the physical property change of the detected object at the same position when the three instantaneous messages are obviously changed at the same position. Because the resolution of the instantaneous phase spectrum is highest and the variation of the instantaneous frequency spectrum and the instantaneous amplitude spectrum is more visual, the approximate position of the underground cable or the abnormality can be determined according to the instantaneous frequency spectrum and the instantaneous amplitude spectrum, and then the cable contour line and the abnormality position can be accurately determined by utilizing the instantaneous phase spectrum.
The Hilbert transform can solve the problem that the data processing is difficult greatly because reflected waves, refracted waves, diffracted waves and scattered waves generated in the underground propagation process of high-frequency pulse electromagnetic waves emitted by the geological radar are mutually overlapped due to the nonuniformity of the earth medium. Meanwhile, various interference noises caused by the fact that the geological radar records the reflected wave characteristics by utilizing the broadband can be reduced.
S1025, performing preset image processing on the geological image data.
In the present embodiment, after the conversion of the geologic image data from the time domain to the frequency spectrum domain is completed, a preset image process is performed on the geologic image data. The image processing comprises IIR filtering, FIR filtering, arithmetic operation, deconvolution, offset processing and static correction. Specifically, removing horizontal noise by adopting an IIR horizontal high-pass filtering and FIR background denoising method, and removing low-frequency horizontal band-shaped interference; removing high-frequency interference by adopting IIR horizontal low-pass filtering and FIR horizontal superposition, and carrying out moving average filtering; the elimination of multiple reflections adopts predictive deconvolution; removing diffraction and correcting a layer with a larger inclination angle by adopting offset processing; the compensation phase change adopts static correction.
Optionally, according to the evaluation processing effect, if the effect is good, performing result interpretation and report writing, and if the effect is bad, performing data analysis and processing again until a better effect is achieved.
S103, detecting inflection points in the geological image data.
Further, an inflection point of the cable is detected in the geological image data after the image processing, and in order to detect the inflection point in the geological image data more conveniently, wavelet denoising of at least two levels is performed on the geological image data before the inflection point is detected in the geological image data, and denoising processing is performed. The concrete mode is as follows:
As shown in fig. 2A, at least two levels of wavelet denoising are performed on the geologic image data, so that geologic image data with higher recognition degree is obtained. The wavelet transformation mainly uses the characteristic multi-resolution, decorrelation and basis selection flexibility of the wavelet transformation, so that the wavelet transformation is quite useful in the aspect of denoising geological image data, and the geological image data is clear. After wavelet transformation, different rules are presented under different resolutions, a threshold is set, and wavelet coefficients are adjusted, so that the purpose of wavelet denoising can be achieved.
In the implementation process of the embodiment, wavelet decomposition is performed on the geological image data to obtain at least two levels of first candidate signals and second candidate signals, wherein the frequency of the second candidate signals is higher than that of the first candidate signals. Specifically, when the geological image data is subjected to wavelet decomposition, a first candidate signal with a wavelet coefficient smaller than the 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 noise signals from the second candidate signals of each level to obtain third candidate signals, wherein each level has a corresponding threshold value, and removing signals smaller than the corresponding threshold value from the second candidate signals of each level to obtain third candidate signals. Taking a secondary wavelet as an example, a second threshold is set in a second candidate signal of the secondary 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, and the second candidate signal larger than the second threshold is reserved, so that a third candidate signal is obtained.
Further, the first candidate signal and the third candidate signal with the highest level are reconstructed into geological image data. Specifically, after wavelet decomposition is completed, wavelet is approximately decomposed, wavelet detail decomposition is added after denoising, a noise-removed signal can be obtained, and the wavelet reconstruction of the signal can be performed according to the first candidate signal with the highest level and the third candidate signal with the high frequency coefficients of the first layer to the third layer after quantization treatment in the low frequency coefficient of the wavelet decomposition, so that the original geological image data after removing the noise can be recovered.
Alternatively, the method used for wavelet denoising is a Matlab (matrix & laboratory) algorithm.
Of course, the above two-level wavelet denoising method is merely an example, and when the denoising effect is not satisfied in implementing the embodiment of the present invention, multiple wavelet decomposition 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 wavelet denoising is completed, inflection points are detected in the geological image data, the cable image data is identified, and the position of the underground cable is determined. Wherein the geologic image data is represented in the computer as matrices, one value in each matrix corresponding to one pixel of the geologic image data. In the geological image data matrix, the numerical values in the matrix are traversed through a moving window, and the inflection point can be determined when a place with obvious pixel change is found. Because the data section of the cable pipeline in the ground penetrating radar detection imaging generally shows a hyperbola shape, and the vertex of the hyperbola reflects the cable wiring and the geological image data information which needs 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 with a preset window to detect gray level change information, and the gray level change information is calculated by the following formula:
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 window translation, and I (x, y) is gray scale information of geological image data;
further, determining an inflection point based on the gradation change information includes:
taylor series simplification is carried out on gray level change information:
the equivalent operation of the local tiny movement amount [ u, v ] and the gray level change information can be obtained:
Wherein [ u, v ] is the displacement of the window, M is the 2x2 partial derivative matrix;
the geological image data change caused by window movement is the eigenvalue analysis of the real symmetric matrix M, a corner response function R is set, and whether the pixel is a corner is judged by judging the size of R:
R=detM-k(traceM)2
Wherein detM is determinant of M matrix, traceM is characteristic value of trace operation R on M matrix depending on M.
If the value of the response function is larger than a preset threshold value, determining that an inflection point exists in the window. Specifically, in geologic image data, for corner points, |r| is large; for flat areas, |R| is very small; for edges, R is negative. The position where the inflection point exists can be determined by calculating the value of R in the window.
S104, sequentially connecting inflection points to the cable image data to represent the cable.
As shown in fig. 2C, when the detection of the inflection point is completed, the inflection point is extracted in the cable image data, and the inflection points are sequentially connected to represent the cable. The position of the inflection point is the position of the cable, so that the cable can be positioned by connecting the inflection point.
S105, cable information of the cable is written into the three-dimensional shape data of the region by referring to the cable image data so as to generate a semantic layer.
And writing the cable information of the cable into the three-dimensional shape data of the region according to the cable image data obtained by connecting the inflection points to generate a semantic layer. The cable information of the cable comprises pipeline characteristics, pipeline attribute data of properties and space data taking space positions as parameters, wherein the pipeline attribute data comprise channel names, communication relations, pipe hole numbers, cable lengths, states, trend and the like, and the space data comprise cable (starting points, inflection points and ending points) coordinates, cable middle joint coordinates, installation positions, top burial depths (or elevations) and the like.
Specifically, three-dimensional shape data are stored in cable image data of a three-dimensional shape body surface, indexes are created in the cable image data, index values meeting query conditions are queried according to the indexes, cable information of the cable is extracted, the cable information is written into the three-dimensional shape data of the region, and a three-dimensional semantic layer is generated according to the cable information in the three-dimensional shape data.
S106, loading the three-dimensional shape data to display a semantic layer of the cable in the three-dimensional map of the region.
As shown in fig. 2D, three-dimensional shape data is loaded in the three-dimensional semantic layer to display the semantic layer of the cable in the three-dimensional map of the region. Specifically, a coordinate system is established in the three-dimensional semantic layer, and three-dimensional shape data corresponding to the cable image data is loaded in the Web page so as to display the semantic layer of the cable in the three-dimensional map of the region and represent the three-dimensional trend of the cable.
It should be noted that, webGIS adopted in the drawing of the three-dimensional map and the cable trend display in the embodiment of the invention mainly comprises the steps of realizing 3D digitization of the underground cable by utilizing a3D modeling technology, a computer technology, a network technology, a database technology and the like, referring to the principles of 3D digitization, networking, realisation, simple visualization and high performance, adopting international Web3D standard VRML language and network cross-platform language Java3D for development, providing three-dimensional visual information of the underground cable, and realizing view management of different directions of the underground cable.
The embodiment of the invention provides a cable display method, device and equipment for a three-dimensional map and a storage medium, wherein the method comprises the following steps: loading three-dimensional shape data of a designated 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 inflection points in the geologic image data; sequentially connecting inflection points to the cable image data to represent the cable; writing cable information of the cable into three-dimensional shape data of the region with reference to the cable image data to generate a semantic layer; the three-dimensional shape data is loaded 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 error of cable identification.
It should be noted that, for simplicity of description, the method embodiments are shown as a series of acts, but it should be understood by those skilled in the art that the embodiments are not limited by the order of acts, as some steps may occur in other orders or concurrently in accordance with the embodiments. Further, those skilled in the art will appreciate that the embodiments described in the specification are presently preferred embodiments, and that the acts are not necessarily required by the embodiments of the invention.
Example two
Fig. 3 is a block diagram of a cable display device for a three-dimensional map according to a second embodiment of the present invention, which may specifically include the following modules:
the three-dimensional shape data loading module 301 is used for loading three-dimensional shape data of a designated area;
A geological image data drawing module 302 for drawing geological image data for an underground space of the region from radar data generated when a ground penetrating radar detects a cable in the region;
an inflection point detection module 303 for detecting an inflection point in the geologic image data;
a cable generating module 304, configured to sequentially connect the inflection points to the cable image data to represent the cable;
a semantic layer generating module 305, configured to write cable information of the cable into three-dimensional shape data of the region with reference to the cable image data, so as to generate a semantic layer;
A cable display module 306 for loading the three-dimensional shape data to display the semantic layer of the cable in a three-dimensional map of the region.
Optionally, the geological image data rendering module includes:
The radar data acquisition sub-module is used for acquiring radar data obtained by detecting the ground penetrating radar in the region embedded with the cable;
The two-dimensional address image data acquisition sub-module is used for adjusting the radar data into two-dimensional geological image data;
the zero line sub-module is used for determining a zero line in the geological image data, wherein the zero line is a position where a cable placed on the ground is intersected with a measuring line;
a spectrum domain conversion sub-module, configured to convert the geological image data from a time domain to a spectrum domain;
and the image processing sub-module is used for executing preset image processing on the geological image data.
Optionally, before detecting the inflection point in the geologic image data, further comprising:
And the wavelet denoising module is used for performing at least two levels of wavelet denoising on the geological image data.
Optionally, the wavelet denoising module includes:
the wavelet decomposition sub-module is used for carrying out wavelet decomposition on the address image data to obtain at least two levels of first candidate signals and second candidate signals, and the frequency of the second candidate signals is higher than that of the first candidate signals;
A third candidate signal obtaining sub-module, configured to remove, for the second candidate signal of each level, a signal smaller than a threshold value, to obtain a third candidate signal, where the threshold value 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 level change information calculation sub-module is used for translating the geological image data with a preset window, and calculating gray level change information through the following formula:
wherein E (u, v) is gray level change information, w (x, y) is a window, I (x+u, y+v) is gray level information after the window is translated, and I (x, y) is gray level information of the geological image data;
And the inflection point determining submodule is used for determining an inflection point based on the gray change information.
Optionally, the image processing includes at least one of:
Gain processing, IIR filtering, FIR filtering, arithmetic operations, deconvolution, offset processing, static correction.
Optionally, the semantic layer generating module includes:
an index creation sub-module for creating an index in the cable image data;
the cable information extraction sub-module is used for inquiring index values meeting inquiry conditions according to the indexes and extracting cable information of the cable;
The cable information writing sub-module is used for writing the cable information into the three-dimensional shape data of the region;
And the semantic layer generation sub-module 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 the 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 merely an example and should not be construed as limiting the functionality and scope of use of embodiments of the present invention.
As shown in FIG. 4, the computer device 12 is in the form of a general purpose computing device. Components of computer device 12 may include, but are not limited to: one or more processors or processing units 16, a system memory 28, a bus 18 that connects the various system components, including the system memory 28 and the processing units 16.
Bus 18 represents one or more of several types of bus structures, including a memory bus or memory controller, a peripheral bus, an accelerated graphics port, a processor, and a local bus using any of a variety of bus architectures. By way of example, and not limitation, such architectures include 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 can 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. The 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 or write to non-removable, nonvolatile magnetic media (not shown in FIG. 4, commonly referred to as a "hard disk drive"). Although not shown in fig. 4, a magnetic disk drive for reading from and writing to a removable non-volatile magnetic disk (e.g., a "floppy disk"), and an optical disk drive for reading from or writing to a removable non-volatile optical disk (e.g., a CD-ROM, DVD-ROM, or other optical media) may be provided. In such cases, each drive may be coupled to bus 18 through one or more data medium interfaces. Memory 28 may include at least one program product having a set (e.g., at least one) of program modules 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 in, for example, 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 or some combination of which may include an implementation of a network environment. Program modules 42 generally perform the functions and/or methods of the embodiments described herein.
The computer device 12 may also communicate with one or more external devices 14 (e.g., keyboard, pointing device, display 24, etc.), one or more devices that enable a user to interact with the computer device 12, and/or any devices (e.g., network card, modem, etc.) that enable the computer device 12 to communicate with one or more other computing devices. Such communication may occur through an input/output (I/O) interface 22. Moreover, computer device 12 may also communicate with one or more networks such as a Local Area Network (LAN), a Wide Area Network (WAN) and/or a public network, such as the Internet, through network adapter 20. As shown, network adapter 20 communicates with other modules of computer device 12 via bus 18. It should be appreciated that although not shown, other hardware and/or software modules may be used in connection with computer device 12, including, but not limited to: microcode, device drivers, redundant processing units, external disk drive arrays, RAID systems, tape drives, data backup storage systems, and the like.
The processing unit 16 executes various functional applications and data processing by running programs stored in the system memory 28, for example, implementing a cable display method of a three-dimensional map provided by an embodiment of the present invention.
Example IV
The fourth embodiment of the present invention further provides a computer readable storage medium, on which a computer program is stored, where the computer program when executed by a processor implements each process of the cable display method of the three-dimensional map, and the same technical effects can be achieved, so that repetition is avoided, and no further description is given here.
The 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 a combination of any 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 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.
Note that the above is only a preferred embodiment of the present invention and the technical principle applied. 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, while the invention has been described in connection with the above embodiments, the invention is not limited to the embodiments, but may be embodied in many other equivalent forms without departing from the spirit or scope of the invention, which is set forth in the following claims.

Claims (7)

1. A cable display method of a three-dimensional map, comprising:
loading three-dimensional shape data of a designated area; the three-dimensional shape data are acquired from a geographic information system, and the appointed area is an appointed area of a city;
Drawing geological image data for an underground space of the region according to radar data generated when a ground penetrating radar detects a cable in the region;
Detecting inflection points in the geologic image data;
Sequentially connecting the inflection points to the cable image data to represent the cable;
Writing cable information of the cable into three-dimensional shape data of the region with reference to the cable image data to generate a semantic layer;
Loading the three-dimensional shape data to display the semantic layer of the cable in a three-dimensional map of the region, wherein the three-dimensional map is a map constructed by adopting a network geographic information system, and the network geographic information system is a network-based client/server system;
The drawing of geological image data for an underground space of the region from radar data generated when a ground penetrating radar detects a cable in the region includes:
Acquiring radar data obtained by detecting the ground penetrating radar in the region in which the cable is embedded;
adjusting the radar data to two-dimensional geological image data;
determining a zero line in the geological image data, wherein the zero line is a junction position of a cable and a measuring line which are placed on the ground;
converting the geologic image data from a time domain to a frequency spectrum domain;
performing preset image processing on the geological image data;
detecting an inflection point in the geologic image data, comprising:
translating the geological image data with a preset window, and calculating gray level change information according to the following formula:
wherein E (u, v) is gray level change information, w (x, y) is a window, I (x+u, y+v) is gray level information after the window is translated, and I (x, y) is gray level information of the geological image data;
determining an inflection point based on the gray scale variation information;
Writing cable information of the cable into three-dimensional shape data of the region with reference to the cable image data to generate a semantic layer, comprising:
Creating an index in the cable image data;
According to the index value meeting the query condition, extracting the cable information of the cable;
writing the cable information into three-dimensional shape data of the region;
And generating a semantic layer according to the cable information in the three-dimensional body data.
2. The method of claim 1, further comprising, prior to detecting an inflection point in the geologic image data:
and performing at least two levels of wavelet denoising on the geological image data.
3. The method of claim 2, wherein performing at least two levels of wavelet denoising on the geologic image data comprises:
performing wavelet decomposition on the address image data to obtain at least two levels of first candidate signals and second candidate signals, wherein the frequency of the second candidate signals is higher than that of the first candidate signals;
removing signals smaller than a threshold value from the second candidate signals of each level to obtain third candidate signals, wherein the threshold value corresponds to each level;
Reconstructing the first candidate signal and the third candidate signal with the highest level into the geological image data.
4. The method of claim 1, wherein the image processing comprises at least one of:
Gain processing, IIR filtering, FIR filtering, arithmetic operations, deconvolution, offset processing, static correction.
5. A cable display device of a three-dimensional map, comprising:
The three-dimensional shape data loading module is used for loading three-dimensional shape data of a designated area; the three-dimensional shape data are acquired from a geographic information system, and the appointed area is an appointed area of a city;
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 inflection points in the geological image data;
The cable generating module is used for sequentially connecting the inflection points in the cable image data to represent the cable;
The semantic layer generation module is used for writing cable information of the cable into the three-dimensional shape data of the region by referring to the cable image data so as to generate a semantic layer;
The cable display module is used for loading the three-dimensional shape data to display the semantic layer of the cable in a three-dimensional map of the region, wherein the three-dimensional map is a map constructed by adopting a network geographic information system, and the network geographic information system is a client/server system based on a network;
A geologic image data rendering module, comprising:
The radar data acquisition sub-module is used for acquiring radar data obtained by detecting the ground penetrating radar in the region embedded with the cable;
The two-dimensional address image data acquisition sub-module is used for adjusting the radar data into two-dimensional geological image data;
the zero line sub-module is used for determining a zero line in the geological image data, wherein the zero line is a position where a cable placed on the ground is intersected with a measuring line;
a spectrum domain conversion sub-module, configured to convert the geological image data from a time domain to a spectrum domain;
an image processing sub-module, configured to perform preset image processing on the geological image data;
An inflection point detection module comprising:
The gray level change information calculation sub-module is used for translating the geological image data with a preset window, and calculating gray level change information through the following formula:
wherein E (u, v) is gray level change information, w (x, y) is a window, I (x+u, y+v) is gray level information after the window is translated, and I (x, y) is gray level information of the geological image data;
an inflection point determining sub-module for determining an inflection point based on the gray scale variation information;
the semantic layer generation module comprises:
an index creation sub-module for creating an index in the cable image data;
the cable information extraction sub-module is used for inquiring index values meeting inquiry conditions according to the indexes and extracting cable information of the cable;
The cable information writing sub-module is used for writing the cable information into the three-dimensional shape data of the region;
And the semantic layer generation sub-module is used for generating a semantic layer according to the cable information in the three-dimensional body data.
6. 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,
The 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 recited in any of claims 1-4.
7. A computer-readable storage medium for realizing 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, realizes the cable display method of a three-dimensional map as claimed in any one of claims 1 to 4.
CN202110129852.0A 2021-01-29 2021-01-29 Cable display method, device and equipment for three-dimensional map and storage medium Active CN112860832B (en)

Priority Applications (1)

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

Applications Claiming Priority (1)

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

Publications (2)

Publication Number Publication Date
CN112860832A CN112860832A (en) 2021-05-28
CN112860832B true CN112860832B (en) 2024-06-11

Family

ID=75987090

Family Applications (1)

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

Country Status (1)

Country Link
CN (1) CN112860832B (en)

Families Citing this family (1)

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

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102944818A (en) * 2012-11-27 2013-02-27 广东电网公司佛山供电局 Geographic information system (GIS) terminal cable fault on-line location method
CN109031464A (en) * 2018-07-05 2018-12-18 国网福建省电力有限公司 Buried cable 3-dimensional digital visable positioning method based on AR and BIM
JP2020085871A (en) * 2018-11-30 2020-06-04 すみえ 池田 Underground radar survey device
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 (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102944818A (en) * 2012-11-27 2013-02-27 广东电网公司佛山供电局 Geographic information system (GIS) terminal cable fault on-line location method
CN109031464A (en) * 2018-07-05 2018-12-18 国网福建省电力有限公司 Buried cable 3-dimensional digital visable positioning method based on AR and BIM
JP2020085871A (en) * 2018-11-30 2020-06-04 すみえ 池田 Underground radar survey device
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
基于探地雷达***的地下管线自动识别技术;涂碧君;中国优秀博硕士学位论文全文数据库(硕士)信息科技辑(第01期);第1、10-13、40-42页 *
海底电缆工程多源勘察数据管理及三维可视化研究;何旭涛等;南方能源建设;第7卷;第2-5页 *

Also Published As

Publication number Publication date
CN112860832A (en) 2021-05-28

Similar Documents

Publication Publication Date Title
US8094515B2 (en) Seismic data visualizations
CN112862888A (en) Cable positioning method and device, computer equipment and storage medium
CA2818790C (en) Seismic trace attribute
CN102305944A (en) Underwater sedimentation amount detecting analysis processing system
US20120320712A1 (en) Dip seismic attribute
US9335434B2 (en) Method and system for processing data acquired in an electromagnetic survey
CN115100363B (en) Underground abnormal body three-dimensional modeling method and device based on ground penetrating radar
CN112860832B (en) Cable display method, device and equipment for three-dimensional map and storage medium
CN113075706A (en) GNSS-R based snow depth inversion method and application thereof
MX2015002330A (en) Measurement transformation apparatus, methods, and systems.
CN111582114A (en) Seismic fault identification method, device, equipment and storage medium
AU2012397811A1 (en) Measurement correction apparatus, methods, and systems
CN115859081A (en) Visual detection method and device for pipeline
CN113568041B (en) Repeatability analysis method and system for time-lapse seismic three-dimensional towing cable acquired data
Sun et al. Research on detection and visualization of underground pipelines
JP3409001B2 (en) 3D voxel data display method and apparatus
JP3409002B2 (en) Three-dimensional exploration method and apparatus
WO2014107257A1 (en) Dip seismic attributes
CN115128675A (en) Method, device, equipment and storage medium for suppressing interbed multiples
Wang et al. Characteristic parameters extraction method of hidden Karst Cave from Borehole radar signal
CN111880223B (en) Attribute analysis method for well-side sound wave remote detection imaging fracture-cave carving
CHOI et al. Integrating GIS and 3D geostatistical methods for geotechnical characterization of soil properties
US11592589B2 (en) Seismic attribute map for gas detection
JP4059609B2 (en) 3D voxel data display method
Zhang et al. Migration imaging processing of subgrade cavity GPR detection based on TUFK method

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
GR01 Patent grant
GR01 Patent grant