CN111160204A - Geological radar image identification method and system based on principal component analysis BP neural network - Google Patents
Geological radar image identification method and system based on principal component analysis BP neural network Download PDFInfo
- Publication number
- CN111160204A CN111160204A CN201911342790.0A CN201911342790A CN111160204A CN 111160204 A CN111160204 A CN 111160204A CN 201911342790 A CN201911342790 A CN 201911342790A CN 111160204 A CN111160204 A CN 111160204A
- Authority
- CN
- China
- Prior art keywords
- geological radar
- radar image
- sample data
- data set
- neural network
- 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.)
- Granted
Links
- 238000013528 artificial neural network Methods 0.000 title claims abstract description 58
- 238000000034 method Methods 0.000 title claims abstract description 46
- 238000000513 principal component analysis Methods 0.000 title claims abstract description 39
- 238000003708 edge detection Methods 0.000 claims abstract description 36
- 230000009467 reduction Effects 0.000 claims abstract description 28
- 238000012549 training Methods 0.000 claims abstract description 21
- 230000000877 morphologic effect Effects 0.000 claims abstract description 20
- 230000008030 elimination Effects 0.000 claims abstract description 18
- 238000003379 elimination reaction Methods 0.000 claims abstract description 18
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 17
- 125000004122 cyclic group Chemical group 0.000 claims abstract description 14
- 239000011435 rock Substances 0.000 claims abstract description 13
- 238000012545 processing Methods 0.000 claims description 18
- 230000008569 process Effects 0.000 claims description 17
- 238000004590 computer program Methods 0.000 claims description 14
- 238000001514 detection method Methods 0.000 claims description 14
- 239000011159 matrix material Substances 0.000 claims description 12
- 238000001914 filtration Methods 0.000 claims description 8
- 238000010276 construction Methods 0.000 claims description 6
- 230000009977 dual effect Effects 0.000 claims description 6
- 230000002708 enhancing effect Effects 0.000 claims description 5
- 230000005540 biological transmission Effects 0.000 claims description 4
- 238000000605 extraction Methods 0.000 claims description 4
- 238000009499 grossing Methods 0.000 claims description 4
- 230000004044 response Effects 0.000 claims description 4
- 230000000007 visual effect Effects 0.000 claims description 4
- 235000019994 cava Nutrition 0.000 claims description 3
- 238000002372 labelling Methods 0.000 claims description 3
- 230000002401 inhibitory effect Effects 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 9
- 239000010410 layer Substances 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 230000007547 defect Effects 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 230000002159 abnormal effect Effects 0.000 description 3
- 230000005856 abnormality Effects 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 2
- 230000001066 destructive effect Effects 0.000 description 2
- 239000000945 filler Substances 0.000 description 2
- 239000011229 interlayer Substances 0.000 description 2
- 238000011835 investigation Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000000149 penetrating effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/10—Terrestrial scenes
- G06V20/13—Satellite images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
- G06N3/084—Backpropagation, e.g. using gradient descent
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/20—Image preprocessing
- G06V10/28—Quantising the image, e.g. histogram thresholding for discrimination between background and foreground patterns
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/20—Image preprocessing
- G06V10/30—Noise filtering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/40—Extraction of image or video features
- G06V10/44—Local feature extraction by analysis of parts of the pattern, e.g. by detecting edges, contours, loops, corners, strokes or intersections; Connectivity analysis, e.g. of connected components
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Artificial Intelligence (AREA)
- Multimedia (AREA)
- Evolutionary Computation (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Health & Medical Sciences (AREA)
- Computing Systems (AREA)
- Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Molecular Biology (AREA)
- Evolutionary Biology (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Remote Sensing (AREA)
- Astronomy & Astrophysics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention provides a geological radar image identification method and system based on a principal component analysis (BP) neural network. The geological radar image identification method based on the principal component analysis BP neural network comprises the steps of marking a label of a geological radar image, wherein the label comprises complete rocks, a fault fracture zone, a fissure zone, a water-rich zone and a karst cave; sequentially carrying out noise elimination, binarization and morphological edge detection on the geological radar image marked with the label to obtain a digital image and form a sample data set; reducing the dimensionality of the sample data set by using a principal component analysis algorithm, and simultaneously keeping the characteristics of the sample data set which have the maximum contribution to the variance; performing cyclic training on the BP neural network by using the sample data set after dimensionality reduction; and receiving a geological radar image in real time, sequentially performing noise elimination, binarization, morphological edge detection and dimension reduction, inputting the image into a trained BP neural network, and outputting a geological radar image recognition result.
Description
Technical Field
The invention belongs to the field of geological radar image processing, and particularly relates to a geological radar image identification method and system based on principal component analysis (BP) neural network.
Background
The statements in this section merely provide background information related to the present disclosure and may not necessarily constitute prior art.
In the prior tunnel advanced geological prediction methods, the geological radar has the characteristics of high efficiency, no damage to targets, high detection data resolution, strong anti-interference capability and the like, and is widely applied. The principle of geological radar detection is that when an abnormal medium and a medium around the abnormal medium are different in electrical property, pulse electromagnetic waves emitted by the geological radar are reflected when being transmitted to an abnormal interface, and reflected signals are received and recorded by a receiving antenna. By analyzing the reflected signal, information such as the spatial position and the buried depth of the anomaly can be inferred.
The inventor finds that in the actual detection process of the geological radar, the working environment is usually disordered, and the influence of noise is added, so that the imaging quality of the radar is low, and the abnormality is difficult to accurately identify, so that the accuracy of abnormality detection excessively depends on the experience and level of technicians; and when the tunnel is long or the data volume is large, a large amount of time and manpower are consumed for abnormity identification, so that the application and popularization of the geological radar technology are limited to a certain extent.
Disclosure of Invention
In order to solve the problems, the invention provides a geological radar image recognition method and system based on a principal component analysis BP neural network, which realize intelligent recognition of geological radar image abnormity, greatly improve recognition accuracy, avoid the defect of recognition depending on experience, save time and improve efficiency.
In order to achieve the purpose, the invention adopts the following technical scheme:
the invention provides a geological radar image recognition method based on principal component analysis (BP) neural network, which comprises the following steps:
labeling a label of a geological radar image, wherein the label comprises complete rocks, a fault fracture zone, a fissure zone, a water-rich zone and a karst cave;
sequentially carrying out noise elimination, binarization and morphological edge detection on the geological radar image marked with the label to obtain a digital image and form a sample data set;
reducing the dimensionality of the sample data set by using a principal component analysis algorithm, and simultaneously keeping the characteristics of the sample data set which have the maximum contribution to the variance;
performing cyclic training on the BP neural network by using the sample data set after dimensionality reduction;
and receiving a geological radar image in real time, sequentially performing noise elimination, binarization, morphological edge detection and dimension reduction, inputting the image into a trained BP neural network, and outputting a geological radar image recognition result.
A second aspect of the present invention provides a geological radar image recognition system based on principal component analysis BP neural network, comprising:
the image label marking module is used for marking labels of the geological radar image, and the labels comprise complete rocks, fault fracture zones, water-rich zones and karst caves;
the sample data set construction module is used for sequentially carrying out noise elimination, binarization and morphological edge detection on the geological radar image marked with the label to obtain a digital image and form a sample data set;
a dimensionality reduction module for reducing the dimensionality of the sample data set using a principal component analysis algorithm while preserving features in the sample data set that contribute most to variance;
the BP neural network training module is used for carrying out cyclic training on the BP neural network by using the sample data set after dimensionality reduction;
and the image real-time identification module is used for receiving the geological radar image in real time, sequentially performing noise elimination, binarization, morphological edge detection and dimension reduction processing, inputting the processed image into a trained BP neural network, and outputting a geological radar image identification result.
A third aspect of the present invention provides a computer readable storage medium, on which a computer program is stored, which when executed by a processor, implements the steps in the method for geological radar image recognition based on principal component analysis, BP, neural network as described above.
A fourth aspect of the present invention provides a computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor when executing the program implements the steps of the geological radar image recognition method based on principal component analysis BP neural network as described above.
The invention has the beneficial effects that:
sequentially carrying out noise elimination, binarization and morphological edge detection on a geological radar image marked with a label to obtain a digital image and form a sample data set; reducing the dimensionality of the sample data set by using a principal component analysis algorithm, and simultaneously keeping the characteristics of the sample data set which have the maximum contribution to the variance; performing cyclic training on the BP neural network by using the sample data set after dimensionality reduction; and finally, receiving a geological radar image in real time, sequentially performing noise elimination, binarization, morphological edge detection and dimension reduction, inputting the image into a BP neural network after training, and outputting a geological radar image recognition result, so that the intelligent recognition of geological radar image abnormity is realized, the recognition accuracy is greatly improved, the defect of recognition by depending on experience is avoided, the time is saved, and the efficiency is improved.
Drawings
The accompanying drawings, which are incorporated in and constitute a part of this specification, are included to provide a further understanding of the invention, and are incorporated in and constitute a part of this specification, illustrate exemplary embodiments of the invention and together with the description serve to explain the invention and not to limit the invention.
FIG. 1 is a flow chart of a geological radar image recognition method based on principal component analysis (BP) neural network provided by the embodiment of the invention;
fig. 2 is a schematic diagram of a BP neural network with a 5-layer structure according to an embodiment of the present invention.
Detailed Description
The invention is further described with reference to the following figures and examples.
It is to be understood that the following detailed description is exemplary and is intended to provide further explanation of the invention as claimed. Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs.
It is noted that the terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of exemplary embodiments according to the invention. As used herein, the singular forms "a", "an" and "the" are intended to include the plural forms as well, and it should be understood that when the terms "comprises" and/or "comprising" are used in this specification, they specify the presence of stated features, steps, operations, devices, components, and/or combinations thereof, unless the context clearly indicates otherwise.
Example one
As shown in fig. 1, the present embodiment provides a geological radar image recognition method based on principal component analysis BP neural network, which includes:
step S101: and labeling a label of the geological radar image, wherein the label comprises complete rocks, a fault fracture zone, a fissure zone, a water-rich zone and a karst cave.
In specific implementation, the waveform characteristics of the geological radar image and the change rule of the frequency, the amplitude, the phase and the electromagnetic wave energy absorption condition and other detailed characteristics can represent different geological phenomena. Among them, in engineering investigation, the common unfavorable geological phenomena are: fault fracture zone, fissure zone, water-rich zone, karst cave and lithology change zone.
The general medium of intact rock mass is relatively even, and the electrical property difference is very little, does not have obvious reflection interface, and radar image and wave form characteristic usually show: the energy clusters are uniformly distributed or strong reflection fine and bright stripes exist only in local parts; generally, a low-amplitude reflection wave group is formed, the waveform is uniform, disordered reflection is avoided, and the automatic gain gradient is relatively small.
The fault is a destructive geological structure, broken rock mass, mud or underground water and the like usually develop in the fault, the medium is extremely uneven, the electric difference is large, and the rock masses at two sides of the fault usually develop joints and folds and have poor medium uniformity. The fracture zone usually exists in fault affected zone, dike and weak interlayer, and the fracture also has various non-uniform fillings and large dielectric difference. The geological radar image shows that the fault and crack interface reflection is strong, the amplitude near the reflecting surface is obviously enhanced and greatly changed, the energy mass is unevenly distributed, diffraction and scattering are often generated in a broken zone and a crack zone, the waveform is disordered, the in-phase axis is broken, and the deep part is even blurred.
The water-rich zone appears in the geological radar image as follows: the geological radar wave generates strong amplitude reflection on the surface of the aquifer, and the electromagnetic wave generates multiple strong reflection with a certain rule when penetrating the aquifer, so as to generate diffraction and scattering in the water-rich zone.
The karst cave appears in geological radar images as follows: the cavity is composed of a plurality of hyperbolic strong reflection waves, and multiple reflection wave groups with high amplitude, low frequency and equal spacing are generally arranged on the side wall of the cavity, and particularly the reflection waves are stronger when no filler is filled or water is filled.
Step S102: and sequentially carrying out noise elimination, binarization and morphological edge detection on the geological radar image marked with the label to obtain a digital image and form a sample data set.
Specifically, before denoising the geological radar image, the method further comprises: and carrying out image enhancement processing on the original geological radar image, enhancing useful information in the image and improving the visual effect of the image.
In specific implementation, noise in the original geological radar image is eliminated by adopting low-pass filtering, and the influence of noise generated by interference in the acquisition and transmission processes on the radar image is reduced; and high-pass filtering is adopted to enhance the high-frequency signal information of the target body outline in the geological radar image, so that the useful image characteristics are highlighted.
Specifically, edge extraction is carried out on the geological radar image by using a Canny edge detection algorithm; the process is as follows:
reducing noise by using Gaussian smoothing processing;
calculating the gradient strength and direction of each pixel point in the image, and reserving the maximum gradient value and direction at each point;
utilizing non-maximum value to inhibit and eliminate spurious response caused by edge detection;
determining true and potential edges using dual threshold detection; the specific process of determining true and potential edges using dual threshold detection is to use a large threshold to detect more confident edge points. The tracking is then performed along the previously derived gradient direction, using a smaller threshold in the tracking, until the original starting point is returned. Obtaining a binary image, wherein each point represents whether an edge point exists or not;
edge detection is finally accomplished by suppressing isolated weak edges.
Among them, the Canny edge detection algorithm is optimal for step-type edges affected by white noise. The purpose of the Canny edge detection algorithm is to return a binary image, with non-zero values indicating the presence of edges in the image, and scale and direction information relating to the edges.
Step S103: the dimensionality of the sample data set is reduced using a principal component analysis algorithm while preserving features in the sample data set that contribute most to variance.
In this embodiment, the process of reducing the dimensionality of the sample data set by using the principal component analysis algorithm while maintaining the features in the sample data set that contribute most to the variance is as follows:
normalizing the samples in the sample data set;
solving a covariance matrix of the sample characteristics;
selecting k maximum eigenvalues to form an eigenvector matrix; wherein k is a positive integer greater than or equal to 2;
and projecting the sample data onto the eigenvector matrix to determine principal components.
Thus, the original multidimensional problem is reduced in dimensionality and greatly simplified.
Step S104: and performing cyclic training on the BP neural network by using the sample data set after dimensionality reduction.
For example: the BP neural network of this embodiment is a 5-layer BP neural network structure, as shown in fig. 2.
Specifically, the process of performing cyclic training on the BP neural network by using the sample data set after dimensionality reduction comprises the following steps:
initializing a weight value and a threshold value, wherein the weight value and the threshold value are random values in a (-1, 1) interval;
the input signal is transmitted in the forward direction, the square of the detection error is calculated to correct the weight and the threshold, the error information is transmitted in the reverse direction from the output layer, and each weight is corrected to reduce the error;
and when the square error is smaller than the preset target error value, finishing iteration and outputting a weight vector, otherwise, continuously transmitting the input signal in the forward direction until the square error is smaller than the preset target error value or the preset iteration times are reached.
Step S105: and receiving a geological radar image in real time, sequentially performing noise elimination, binarization, morphological edge detection and dimension reduction, inputting the image into a trained BP neural network, and outputting a geological radar image recognition result.
In the embodiment, denoising, binarization and morphological edge detection processing are sequentially carried out on a geological radar image marked with a label to obtain a digital image and form a sample data set; reducing the dimensionality of the sample data set by using a principal component analysis algorithm, and simultaneously keeping the characteristics of the sample data set which have the maximum contribution to the variance; performing cyclic training on the BP neural network by using the sample data set after dimensionality reduction; and finally, receiving a geological radar image in real time, sequentially performing noise elimination, binarization, morphological edge detection and dimension reduction, inputting the image into a BP neural network after training, and outputting a geological radar image recognition result, so that the intelligent recognition of geological radar image abnormity is realized, the recognition accuracy is greatly improved, the defect of recognition by depending on experience is avoided, the time is saved, and the efficiency is improved.
Example two
The embodiment provides a geological radar image recognition system based on principal component analysis (BP) neural network, which comprises:
(1) the image label marking module is used for marking labels of the geological radar image, and the labels comprise complete rocks, fault fracture zones, water-rich zones and karst caves;
in specific implementation, the waveform characteristics of the geological radar image and the change rule of the frequency, the amplitude, the phase and the electromagnetic wave energy absorption condition and other detailed characteristics can represent different geological phenomena. Among them, in engineering investigation, the common unfavorable geological phenomena are: fault fracture zone, fissure zone, water-rich zone, karst cave and lithology change zone.
The general medium of intact rock mass is relatively even, and the electrical property difference is very little, does not have obvious reflection interface, and radar image and wave form characteristic usually show: the energy clusters are uniformly distributed or strong reflection fine and bright stripes exist only in local parts; generally, a low-amplitude reflection wave group is formed, the waveform is uniform, disordered reflection is avoided, and the automatic gain gradient is relatively small.
The fault is a destructive geological structure, broken rock mass, mud or underground water and the like usually develop in the fault, the medium is extremely uneven, the electric difference is large, and the rock masses at two sides of the fault usually develop joints and folds and have poor medium uniformity. The fracture zone usually exists in fault affected zone, dike and weak interlayer, and the fracture also has various non-uniform fillings and large dielectric difference. The geological radar image shows that the fault and crack interface reflection is strong, the amplitude near the reflecting surface is obviously enhanced and greatly changed, the energy mass is unevenly distributed, diffraction and scattering are often generated in a broken zone and a crack zone, the waveform is disordered, the in-phase axis is broken, and the deep part is even blurred.
The water-rich zone appears in the geological radar image as follows: the geological radar wave generates strong amplitude reflection on the surface of the aquifer, and the electromagnetic wave generates multiple strong reflection with a certain rule when penetrating the aquifer, so as to generate diffraction and scattering in the water-rich zone.
The karst cave appears in geological radar images as follows: the cavity is composed of a plurality of hyperbolic strong reflection waves, and multiple reflection wave groups with high amplitude, low frequency and equal spacing are generally arranged on the side wall of the cavity, and particularly the reflection waves are stronger when no filler is filled or water is filled.
(2) The sample data set construction module is used for sequentially carrying out noise elimination, binarization and morphological edge detection on the geological radar image marked with the label to obtain a digital image and form a sample data set;
specifically, before denoising the geological radar image, the method further comprises: and carrying out image enhancement processing on the original geological radar image, enhancing useful information in the image and improving the visual effect of the image.
In specific implementation, noise in the original geological radar image is eliminated by adopting low-pass filtering, and the influence of noise generated by interference in the acquisition and transmission processes on the radar image is reduced; and high-pass filtering is adopted to enhance the high-frequency signal information of the target body outline in the geological radar image, so that the useful image characteristics are highlighted.
Specifically, edge extraction is carried out on the geological radar image by using a Canny edge detection algorithm; the process is as follows:
reducing noise by using Gaussian smoothing processing;
calculating the gradient strength and direction of each pixel point in the image, and reserving the maximum gradient value and direction at each point;
utilizing non-maximum value to inhibit and eliminate spurious response caused by edge detection;
determining true and potential edges using dual threshold detection; the specific process of determining true and potential edges using dual threshold detection is to use a large threshold to detect more confident edge points. The tracking is then performed along the previously derived gradient direction, using a smaller threshold in the tracking, until the original starting point is returned. Obtaining a binary image, wherein each point represents whether an edge point exists or not;
edge detection is finally accomplished by suppressing isolated weak edges.
Among them, the Canny edge detection algorithm is optimal for step-type edges affected by white noise. The purpose of the Canny edge detection algorithm is to return a binary image, with non-zero values indicating the presence of edges in the image, and scale and direction information relating to the edges.
(3) A dimensionality reduction module for reducing the dimensionality of the sample data set using a principal component analysis algorithm while preserving features in the sample data set that contribute most to variance;
in this embodiment, the process of reducing the dimensionality of the sample data set by using the principal component analysis algorithm while maintaining the features in the sample data set that contribute most to the variance is as follows:
normalizing the samples in the sample data set;
solving a covariance matrix of the sample characteristics;
selecting k maximum eigenvalues to form an eigenvector matrix; wherein k is a positive integer greater than or equal to 2;
and projecting the sample data onto the eigenvector matrix to determine principal components.
Thus, the original multidimensional problem is reduced in dimensionality and greatly simplified.
(4) The BP neural network training module is used for carrying out cyclic training on the BP neural network by using the sample data set after dimensionality reduction;
for example: the BP neural network of this embodiment is a 5-layer BP neural network structure, as shown in fig. 2.
Specifically, the process of performing cyclic training on the BP neural network by using the sample data set after dimensionality reduction comprises the following steps:
initializing a weight value and a threshold value, wherein the weight value and the threshold value are random values in a (-1, 1) interval;
the input signal is transmitted in the forward direction, the square of the detection error is calculated to correct the weight and the threshold, the error information is transmitted in the reverse direction from the output layer, and each weight is corrected to reduce the error;
and when the square error is smaller than the preset target error value, finishing iteration and outputting a weight vector, otherwise, continuously transmitting the input signal in the forward direction until the square error is smaller than the preset target error value or the preset iteration times are reached.
(5) And the image real-time identification module is used for receiving the geological radar image in real time, sequentially performing noise elimination, binarization, morphological edge detection and dimension reduction processing, inputting the processed image into a trained BP neural network, and outputting a geological radar image identification result.
In the embodiment, denoising, binarization and morphological edge detection processing are sequentially carried out on a geological radar image marked with a label to obtain a digital image and form a sample data set; reducing the dimensionality of the sample data set by using a principal component analysis algorithm, and simultaneously keeping the characteristics of the sample data set which have the maximum contribution to the variance; performing cyclic training on the BP neural network by using the sample data set after dimensionality reduction; and finally, receiving a geological radar image in real time, sequentially performing noise elimination, binarization, morphological edge detection and dimension reduction, inputting the image into a BP neural network after training, and outputting a geological radar image recognition result, so that the intelligent recognition of geological radar image abnormity is realized, the recognition accuracy is greatly improved, the defect of recognition by depending on experience is avoided, the time is saved, and the efficiency is improved.
EXAMPLE III
The present embodiment provides a computer-readable storage medium, on which a computer program is stored, which when executed by a processor implements the steps in the principal component analysis BP neural network-based geological radar image recognition method as described above.
Example four
The embodiment provides a computer device, which comprises a memory, a processor and a computer program stored on the memory and capable of running on the processor, wherein the processor executes the program to realize the steps in the geological radar image recognition method based on the principal component analysis (BP) neural network.
As will be appreciated by one skilled in the art, embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of a hardware embodiment, a software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, optical storage, and the like) having computer-usable program code embodied therein.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
It will be understood by those skilled in the art that all or part of the processes of the methods of the embodiments described above can be implemented by a computer program, which can be stored in a computer-readable storage medium, and when executed, can include the processes of the embodiments of the methods described above. The storage medium may be a magnetic disk, an optical disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), or the like.
The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention, and various modifications and changes may be made by those skilled in the art. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Claims (10)
1. A geological radar image recognition method based on principal component analysis (BP) neural network is characterized by comprising the following steps:
labeling a label of a geological radar image, wherein the label comprises complete rocks, a fault fracture zone, a fissure zone, a water-rich zone and a karst cave;
sequentially carrying out noise elimination, binarization and morphological edge detection on the geological radar image marked with the label to obtain a digital image and form a sample data set;
reducing the dimensionality of the sample data set by using a principal component analysis algorithm, and simultaneously keeping the characteristics of the sample data set which have the maximum contribution to the variance;
performing cyclic training on the BP neural network by using the sample data set after dimensionality reduction;
and receiving a geological radar image in real time, sequentially performing noise elimination, binarization, morphological edge detection and dimension reduction, inputting the image into a trained BP neural network, and outputting a geological radar image recognition result.
2. The geological radar image recognition method based on principal component analysis (BP) neural network as claimed in claim 1, wherein before denoising the geological radar image, further comprising: and carrying out image enhancement processing on the original geological radar image, enhancing useful information in the image and improving the visual effect of the image.
3. The geological radar image recognition method based on principal component analysis (BP) neural network as claimed in claim 1, characterized in that the noise in the original geological radar image is eliminated by low-pass filtering, and the influence of the noise generated by interference in the acquisition and transmission processes on the radar image is reduced; and high-pass filtering is adopted to enhance the high-frequency signal information of the target body outline in the geological radar image, so that the useful image characteristics are highlighted.
4. The geological radar image recognition method based on principal component analysis (BP) neural network as claimed in claim 1, characterized in that edge extraction is performed on the geological radar image by using Canny edge detection algorithm; the process is as follows:
reducing noise by using Gaussian smoothing processing;
calculating the gradient strength and direction of each pixel point in the image, and reserving the maximum gradient value and direction at each point;
utilizing non-maximum value to inhibit and eliminate spurious response caused by edge detection;
determining true and potential edges using dual threshold detection;
edge detection is finally accomplished by suppressing isolated weak edges.
5. The method of claim 1, wherein the principal component analysis (BP) neural network-based geological radar image recognition method is characterized in that the principal component analysis algorithm is used to reduce the dimensionality of the sample data set while maintaining the features in the sample data set that contribute most to the variance by:
normalizing the samples in the sample data set;
solving a covariance matrix of the sample characteristics;
selecting k maximum eigenvalues to form an eigenvector matrix; wherein k is a positive integer greater than or equal to 2;
and projecting the sample data onto the eigenvector matrix to determine principal components.
6. The geological radar image recognition method based on principal component analysis (BP) neural network as claimed in claim 1, wherein the process of performing cyclic training on the BP neural network by using the sample data set after dimensionality reduction comprises:
initializing a weight value and a threshold value, wherein the weight value and the threshold value are random values in a (-1, 1) interval;
the input signal is transmitted in the forward direction, the square of the detection error is calculated to correct the weight and the threshold, the error information is transmitted in the reverse direction from the output layer, and each weight is corrected to reduce the error;
and when the square error is smaller than the preset target error value, finishing iteration and outputting a weight vector, otherwise, continuously transmitting the input signal in the forward direction until the square error is smaller than the preset target error value or the preset iteration times are reached.
7. A geological radar image recognition system based on principal component analysis (BP) neural network is characterized by comprising:
the image label marking module is used for marking labels of the geological radar image, and the labels comprise complete rocks, fault fracture zones, water-rich zones and karst caves;
the sample data set construction module is used for sequentially carrying out noise elimination, binarization and morphological edge detection on the geological radar image marked with the label to obtain a digital image and form a sample data set;
a dimensionality reduction module for reducing the dimensionality of the sample data set using a principal component analysis algorithm while preserving features in the sample data set that contribute most to variance;
the BP neural network training module is used for carrying out cyclic training on the BP neural network by using the sample data set after dimensionality reduction;
and the image real-time identification module is used for receiving the geological radar image in real time, sequentially performing noise elimination, binarization, morphological edge detection and dimension reduction processing, inputting the processed image into a trained BP neural network, and outputting a geological radar image identification result.
8. The geological radar image recognition system based on principal component analysis (BP) neural network of claim 7, wherein before denoising geological radar images in the sample data set construction module, further comprising: carrying out image enhancement processing on an original geological radar image, enhancing useful information in the image and improving the visual effect of the image;
or
In the sample data set construction module, noise in the original geological radar image is eliminated by adopting low-pass filtering, and the influence of noise generated by interference in the acquisition and transmission processes on the radar image is reduced; enhancing high-frequency signal information of a target body contour in the geological radar image by adopting high-pass filtering, and highlighting useful image characteristics;
or
In the sample data set construction module, performing edge extraction on a geological radar image by using a Canny edge detection algorithm; the process is as follows:
reducing noise by using Gaussian smoothing processing;
calculating the gradient strength and direction of each pixel point in the image, and reserving the maximum gradient value and direction at each point;
utilizing non-maximum value to inhibit and eliminate spurious response caused by edge detection;
determining true and potential edges using dual threshold detection;
finally completing edge detection by inhibiting isolated weak edges;
or
In the dimensionality reduction module, the dimensionality of the sample data set is reduced by using a principal component analysis algorithm, and meanwhile, the process of keeping the characteristics which have the maximum contribution to the variance in the sample data set is as follows:
normalizing the samples in the sample data set;
solving a covariance matrix of the sample characteristics;
selecting k maximum eigenvalues to form an eigenvector matrix; wherein k is a positive integer greater than or equal to 2;
projecting the sample data to a feature vector matrix to determine principal components;
or
In the BP neural network training module, the process of performing cyclic training on the BP neural network by using the sample data set after dimensionality reduction comprises the following steps:
initializing a weight value and a threshold value, wherein the weight value and the threshold value are random values in a (-1, 1) interval;
the input signal is transmitted in the forward direction, the square of the detection error is calculated to correct the weight and the threshold, the error information is transmitted in the reverse direction from the output layer, and each weight is corrected to reduce the error;
and when the square error is smaller than the preset target error value, finishing iteration and outputting a weight vector, otherwise, continuously transmitting the input signal in the forward direction until the square error is smaller than the preset target error value or the preset iteration times are reached.
9. A computer readable storage medium, having stored thereon a computer program, characterized in that the program, when being executed by a processor, carries out the steps of the method for geological radar image recognition based on principal component analysis, BP, neural network according to any of claims 1-6.
10. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, characterized in that the processor implements the steps in the method for geological radar image recognition based on principal component analysis, BP, neural network according to any of claims 1-6 when executing the program.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911342790.0A CN111160204B (en) | 2019-12-23 | 2019-12-23 | Geological radar image recognition method and system based on principal component analysis BP neural network |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911342790.0A CN111160204B (en) | 2019-12-23 | 2019-12-23 | Geological radar image recognition method and system based on principal component analysis BP neural network |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111160204A true CN111160204A (en) | 2020-05-15 |
| CN111160204B CN111160204B (en) | 2024-01-30 |
Family
ID=70558127
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911342790.0A Active CN111160204B (en) | 2019-12-23 | 2019-12-23 | Geological radar image recognition method and system based on principal component analysis BP neural network |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111160204B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112001885A (en) * | 2020-07-14 | 2020-11-27 | 武汉理工大学 | PCB fault detection method, storage medium and system based on machine vision |
| WO2021068848A1 (en) * | 2019-10-09 | 2021-04-15 | 山东大学 | Tunnel structure disease multi-scale measurement and intelligent diagnosis system and method |
| CN112801035A (en) * | 2021-02-24 | 2021-05-14 | 山东大学 | Knowledge and data dual-drive-based carrying type lithology intelligent identification method and system |
| CN113486971A (en) * | 2021-07-19 | 2021-10-08 | 国网山东省电力公司日照供电公司 | User state identification method and system based on principal component analysis and neural network |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102722759A (en) * | 2012-05-17 | 2012-10-10 | 河海大学 | Method for predicting power supply reliability of power grid based on BP neural network |
| CN104200284A (en) * | 2014-09-01 | 2014-12-10 | 山东科技大学 | Mine fault structure prediction method |
| CN107092990A (en) * | 2017-05-03 | 2017-08-25 | 西安电子科技大学 | The shield construction ground settlement forecast system and method analyzed based on big data |
| US20180060398A1 (en) * | 2016-08-23 | 2018-03-01 | Interana, Inc. | Methods for stratified sampling-based query execution |
| CN108182427A (en) * | 2018-01-30 | 2018-06-19 | 电子科技大学 | A kind of face identification method based on deep learning model and transfer learning |
-
2019
- 2019-12-23 CN CN201911342790.0A patent/CN111160204B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102722759A (en) * | 2012-05-17 | 2012-10-10 | 河海大学 | Method for predicting power supply reliability of power grid based on BP neural network |
| CN104200284A (en) * | 2014-09-01 | 2014-12-10 | 山东科技大学 | Mine fault structure prediction method |
| US20180060398A1 (en) * | 2016-08-23 | 2018-03-01 | Interana, Inc. | Methods for stratified sampling-based query execution |
| CN107092990A (en) * | 2017-05-03 | 2017-08-25 | 西安电子科技大学 | The shield construction ground settlement forecast system and method analyzed based on big data |
| CN108182427A (en) * | 2018-01-30 | 2018-06-19 | 电子科技大学 | A kind of face identification method based on deep learning model and transfer learning |
Non-Patent Citations (1)
| Title |
|---|
| 蒋林权: ""基于机器视觉的砀山梨表面缺陷识别"" * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2021068848A1 (en) * | 2019-10-09 | 2021-04-15 | 山东大学 | Tunnel structure disease multi-scale measurement and intelligent diagnosis system and method |
| US12013485B2 (en) | 2019-10-09 | 2024-06-18 | Shandong University | Multi-scale inspection and intelligent diagnosis system and method for tunnel structural defects |
| CN112001885A (en) * | 2020-07-14 | 2020-11-27 | 武汉理工大学 | PCB fault detection method, storage medium and system based on machine vision |
| CN112801035A (en) * | 2021-02-24 | 2021-05-14 | 山东大学 | Knowledge and data dual-drive-based carrying type lithology intelligent identification method and system |
| CN113486971A (en) * | 2021-07-19 | 2021-10-08 | 国网山东省电力公司日照供电公司 | User state identification method and system based on principal component analysis and neural network |
| CN113486971B (en) * | 2021-07-19 | 2023-10-27 | 国网山东省电力公司日照供电公司 | User state identification method and system based on principal component analysis and neural network |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111160204B (en) | 2024-01-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111160204B (en) | Geological radar image recognition method and system based on principal component analysis BP neural network | |
| CN104991269B (en) | A kind of margin guide and the full waveform inversion fast method of structural constraint | |
| CN102879766B (en) | Method and apparatus for detecting and tracking faint target of high frequency ground wave radar | |
| CN110133643B (en) | Plant root system detection method and device | |
| CN107705313B (en) | Remote sensing image ship target segmentation method | |
| CN107678019A (en) | Radar signal multi-target detection method and device based on CA CFAR | |
| CN109146803B (en) | SAR image radiation resolution improving method and device based on multi-angle image | |
| CN114677525B (en) | Edge detection method based on binary image processing | |
| CN104535973A (en) | Target detection method by use of airborne early warning radar | |
| CN107479037B (en) | PD radar clutter area distinguishing method | |
| CN115343685A (en) | Multi-dimensional ground penetrating radar detection method, device and equipment applied to disease identification | |
| He et al. | GPR image denoising with NSST-UNET and an improved BM3D | |
| CN105652256A (en) | High-frequency surface wave radar TBD method based on polarization information | |
| CN112666552B (en) | Ground penetrating radar data background clutter self-adaptive removing method | |
| CN112713907B (en) | Marine CSEM noise reduction method and system based on dictionary learning | |
| CN115542278A (en) | Ground penetrating radar original data-based pavement damage area judgment method | |
| CN113705516A (en) | Sound source depth discrimination method based on radiation noise interference fringe characteristics | |
| CN111399038B (en) | Slope parameter extraction method and device and computer readable storage medium | |
| CN113506312A (en) | Ultraviolet discharge image segmentation method and computer readable medium | |
| CN110596709A (en) | Method for self-adaptively detecting multiple targets by vehicle-mounted millimeter wave radar | |
| CN112379427B (en) | Method and device for suppressing surface wave noise | |
| CN117152622B (en) | Boundary optimization model training, boundary optimization method, device, equipment and medium | |
| CN114170200B (en) | Metal pitting defect degree evaluation method and system based on artificial intelligence | |
| CN116051426B (en) | Synthetic aperture radar image processing method | |
| CN107607931A (en) | A kind of laser radar echo image processing 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 |