CN115147746B - Saline-alkali geological identification method based on unmanned aerial vehicle remote sensing image - Google Patents
Saline-alkali geological identification method based on unmanned aerial vehicle remote sensing image Download PDFInfo
- Publication number
- CN115147746B CN115147746B CN202211068025.6A CN202211068025A CN115147746B CN 115147746 B CN115147746 B CN 115147746B CN 202211068025 A CN202211068025 A CN 202211068025A CN 115147746 B CN115147746 B CN 115147746B
- Authority
- CN
- China
- Prior art keywords
- target
- saline
- alkali
- region
- preliminary
- 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
Links
- 239000003513 alkali Substances 0.000 title claims abstract description 340
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000012545 processing Methods 0.000 claims abstract description 38
- 238000005457 optimization Methods 0.000 claims abstract description 24
- 238000000605 extraction Methods 0.000 claims abstract description 17
- 238000006243 chemical reaction Methods 0.000 claims abstract description 9
- 238000005520 cutting process Methods 0.000 claims abstract description 9
- 239000002585 base Substances 0.000 claims description 9
- 238000012549 training Methods 0.000 claims description 7
- 230000007704 transition Effects 0.000 claims description 4
- 238000013507 mapping Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- 239000002689 soil Substances 0.000 abstract description 22
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical compound [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 description 7
- 230000011218 segmentation Effects 0.000 description 7
- 239000011780 sodium chloride Substances 0.000 description 7
- 238000004364 calculation method Methods 0.000 description 6
- 238000013528 artificial neural network Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 238000010606 normalization Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 150000003839 salts Chemical class 0.000 description 2
- 239000000758 substrate Substances 0.000 description 2
- 238000012800 visualization Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000001704 evaporation Methods 0.000 description 1
- 230000008020 evaporation Effects 0.000 description 1
- 239000003673 groundwater Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000003911 water pollution 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/17—Terrestrial scenes taken from planes or by drones
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/90—Dynamic range modification of images or parts thereof
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/90—Determination of colour characteristics
-
- 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/26—Segmentation of patterns in the image field; Cutting or merging of image elements to establish the pattern region, e.g. clustering-based techniques; Detection of occlusion
- G06V10/267—Segmentation of patterns in the image field; Cutting or merging of image elements to establish the pattern region, e.g. clustering-based techniques; Detection of occlusion by performing operations on regions, e.g. growing, shrinking or watersheds
-
- 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/56—Extraction of image or video features relating to colour
-
- 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/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
- G06V10/764—Arrangements for image or video recognition or understanding using pattern recognition or machine learning using classification, e.g. of video objects
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20212—Image combination
- G06T2207/20221—Image fusion; Image merging
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Multimedia (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- Health & Medical Sciences (AREA)
- Remote Sensing (AREA)
- Computing Systems (AREA)
- Databases & Information Systems (AREA)
- Evolutionary Computation (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Image Analysis (AREA)
Abstract
The invention relates to the technical field of image data processing, in particular to a saline-alkali geological identification method based on an unmanned aerial vehicle remote sensing image, which comprises the following steps: acquiring a remote sensing image to be detected of a target ground area under the saline-alkali condition to be detected; carrying out image enhancement extraction processing on a remote sensing image to be detected; carrying out color space conversion processing on the target remote sensing image; carrying out region division optimization on the color remote sensing image; merging and dividing the initial optimal saline-alkali area; cutting each target saline-alkali area; inputting each target saline-alkali image into a saline-alkali degree classification network, and determining the saline-alkali degree corresponding to the target saline-alkali image; and generating target saline-alkali information representing the saline-alkali condition of the target ground area. The method solves the technical problem of low accuracy of saline-alkali identification of the soil by carrying out image processing on the remote sensing image to be detected, improves the accuracy of saline-alkali identification of the soil, and is mainly applied to identification of the saline-alkali degree of the soil.
Description
Technical Field
The invention relates to the technical field of image data processing, in particular to a saline-alkali geological identification method based on an unmanned aerial vehicle remote sensing image.
Background
Under natural environment influence or artificial movement interference, the salinity in soil and groundwater rises to the earth's surface through the moisture of capillary along with the bottom, and waits to the moisture evaporation after, the salinity often gathers in earth's surface soil, forms saline and alkaline land. Large-scale geological salinization not only affects the land quality in the range, but also often causes yield reduction in cultivated land and water pollution, and even causes reduction in biodiversity. In addition, geological salinization may also cause many other serious land problems, with very serious direct or indirect economic losses to human society. Therefore, the saline-alkali identification of the soil is very important. At present, when saline-alkali identification is carried out on soil, the mode that usually adopts is: firstly, acquiring a soil ground image, inputting the soil ground image into a semantic segmentation network, extracting and segmenting the soil ground image through the semantic segmentation network, finally, respectively inputting a plurality of extracted and segmented images obtained after extraction and segmentation into a classification neural network, and detecting the saline-alkali degree of the plurality of extracted and segmented images through the classification neural network to obtain the saline-alkali degree corresponding to each extracted and segmented image.
However, when the above-described manner is adopted, there are often technical problems as follows:
due to the characteristics of the semantic segmentation network, the relation between pixels is not considered when the soil ground image is extracted and segmented through the semantic segmentation network, the space consistency is lacked, and the extraction and segmentation are not sensitive to details, so that the soil ground image is not accurate enough in extraction and segmentation, and therefore, the extracted and segmented images of regions with various saline-alkali degrees possibly exist in the obtained plurality of extracted and segmented images, the saline-alkali degrees corresponding to the segmented images cannot be extracted through accurate determination of a subsequent classification neural network, and further the saline-alkali identification of soil is low in accuracy.
Disclosure of Invention
This summary is provided to introduce a selection of concepts in a simplified form that are further described below in the detailed description. This summary is not intended to identify key features or essential features of the claimed subject matter, nor is it intended to be used to limit the scope of the claimed subject matter.
In order to solve the technical problem of low accuracy of saline-alkali identification of soil, the invention provides a saline-alkali geological identification method based on an unmanned aerial vehicle remote sensing image.
The invention provides a saline-alkali geological identification method based on unmanned aerial vehicle remote sensing images, which comprises the following steps:
acquiring a remote sensing image to be detected of a target ground area of a saline-alkali condition to be detected through a target camera installed on a target unmanned aerial vehicle;
carrying out image enhancement extraction processing on the remote sensing image to be detected to obtain a target remote sensing image;
carrying out color space conversion processing on the target remote sensing image to obtain a color remote sensing image;
carrying out region division optimization on the color remote sensing image to obtain a preliminary optimal saline-alkali region set;
merging and dividing the initial optimal saline-alkali regions in the initial optimal saline-alkali region set to obtain a target saline-alkali region set;
cutting each target saline-alkali area in the target saline-alkali area set, determining a target saline-alkali image corresponding to the target saline-alkali area to obtain a target saline-alkali image set;
inputting each target saline-alkali image in the target saline-alkali image set into a trained saline-alkali degree classification network, and determining the saline-alkali degree corresponding to the target saline-alkali image through the saline-alkali degree classification network;
and generating target saline-alkali information representing the saline-alkali condition of the target ground area according to the saline-alkali degree corresponding to each target saline-alkali image in the target saline-alkali image set.
Further, the color remote sensing image is subjected to region division optimization to obtain a preliminary optimal saline-alkali region set, which includes:
mapping the color remote sensing image into a target undirected graph;
and determining the initial optimal saline-alkali area set through an optimization algorithm according to the target undirected graph.
Further, the objective function of the optimization algorithm is:
wherein,is the objective function of the optimization algorithm,Fis the value of the objective function of the optimization algorithm,is a first objective function of the first set of functions,His the value of the first objective function,,so as to makeCounting at the bottomThe number of the pairs is logarithmic,is the second in the target undirected graphiPixel point and the secondjThe probability of a transition between individual pixel points,,,is the second in the target undirected graphiThe sum of the edge weights between each pixel point and each pixel point in the target undirected graph,Iis the number of pixel points in the target undirected graph, the targetSecond in the Infinite drawingiPixel point and the secondjThe edge weight between each pixel point is,exp() Is an exponential function with a natural constant as the base,is the second in the target undirected graphiPixel point and the secondjThe euclidean distance between the individual pixel points,is the second in the target undirected graphiHue value corresponding to each pixel point and the secondjThe difference in hue value corresponding to each pixel point,is the second in the target undirected graphiSaturation value corresponding to each pixel point and the secondjThe difference in saturation values corresponding to each pixel point,is the second in the target undirected graphiBrightness value corresponding to each pixel point and the secondjThe difference in luminance values corresponding to each pixel point,Uis the parameter of the model and is,is the coefficient of the model that is,is the second objective function of the first function,Objis the value of the second objective function,is the second in the target undirected graphkThe number of pixel points within a sub-region, a sub-region being a region in the target undirected graph,andare respectively the second in the target undirected graphkThe length and width of the smallest circumscribed rectangle corresponding to a sub-region,Kis the number of sub-regions into which the target undirected graph is divided.
Further, the preliminary optimal saline-alkali area in the preliminary optimal saline-alkali area set is merged and divided to obtain a target saline-alkali area set, and the method comprises the following steps:
normalizing the target remote sensing image to obtain a target normalized image;
determining a target preliminary region corresponding to each preliminary optimal saline-alkali region in the preliminary optimal saline-alkali region set according to the position of each preliminary optimal saline-alkali region in the color remote sensing image to obtain a target preliminary region set, wherein the target preliminary region in the target preliminary region set is a region in the target normalized image;
dividing a channel value of each channel in a target number of channels of the target normalized image into a preset number of channel levels to obtain a channel level set, wherein each channel of the target normalized image corresponds to the preset number of channel levels, and the number of the channel levels in the channel level set is the power of the preset number of target number;
determining a color channel histogram corresponding to each target preliminary region in the target preliminary region set;
determining a first region merging index between every two target preliminary regions according to the channel level set and color channel histograms corresponding to the two target preliminary regions in the target preliminary region set;
determining a second region merging index between every two target preliminary regions according to every two target preliminary regions in the target preliminary region set;
determining an overall region merging index between every two target preliminary regions according to a first region merging index and a second region merging index between every two target preliminary regions in the target preliminary region set;
when the overall area merging index between two target preliminary areas in the target preliminary area set is larger than a preset judgment threshold value, dividing the two target preliminary areas into the same type of target preliminary area;
and determining the target preliminary region in the same target preliminary region as a target saline-alkali region in the target saline-alkali region set.
Further, the formula for determining the correspondence between the first region merging indicators of the two target preliminary regions is as follows:
wherein,is the first in the target preliminary region setpA target preliminary region andqa first region merging indicator between the target preliminary regions,Cis the number of channel levels in the set of channel levels,is the first in the target preliminary region setpThe first of the channel level sets included in the corresponding color channel histogram for each target preliminary regioncThe histogram distribution values at the level of one channel,is the first in the target preliminary region setqA first one of the channel level sets included in the corresponding color channel histogram for each target preliminary regioncThe histogram distribution values at the level of one channel,is presetA number greater than 0.
Further, the determining a second region merging indicator between every two target preliminary regions according to every two target preliminary regions in the target preliminary region set includes:
translating the two target preliminary regions to enable the distance between the two target preliminary regions to be zero, and obtaining fitting regions corresponding to the two target preliminary regions;
and determining a second region merging index between the two target preliminary regions according to the two target preliminary regions and the corresponding fitting regions of the two target preliminary regions.
Further, the formula for determining the second region merging indicator between the two target preliminary regions is:
wherein,is the first in the target preliminary region setpA target preliminary region andqsecond region merging indicators between the target preliminary regions,exp() Is an exponential function with a natural constant as the base,is the first in the target preliminary region setpA target preliminary region andqthe number of pixel points in the fitting and region corresponding to each target preliminary region,is the first in the target preliminary region setpA target preliminary region andqthe number of edge pixel points included in the fitting and region corresponding to each target preliminary region,is the first in the target preliminary region setpA target preliminary region andqthe perimeter of the minimum bounding rectangle corresponding to the fitting and region corresponding to each target preliminary region,is the first in the target preliminary region setpThe number of pixel points within the individual target initialization region,is the first in the target preliminary region setpThe number of edge pixels included in each target preliminary region,is the first in the target preliminary region setpThe perimeter of the minimum bounding rectangle corresponding to each preliminary region of interest,is the first in the target preliminary region setqThe number of pixel points within the individual target initialization region,is the first in the target preliminary region setqThe number of edge pixels included in each target preliminary region,is the first in the target preliminary region setqThe perimeter of the minimum bounding rectangle corresponding to each target preliminary region.
Further, the formula for determining the correspondence of the overall region merging indicator between the two target preliminary regions is as follows:
wherein,is the first in the target preliminary region setpA target preliminary region andqthe overall region merging index between the target preliminary regions,exp() Is an exponential function with a natural constant as the base,is the first in the target preliminary region setpA target preliminary region andqa first region merging indicator between the target preliminary regions,is the first in the target preliminary region setpA target preliminary region andqsecond region merging indicators between the target preliminary regions.
Further, the image enhancement extraction processing is performed on the remote sensing image to be detected to obtain a target remote sensing image, and the image enhancement extraction processing includes:
carrying out image enhancement processing on the remote sensing image to be detected through an image enhancement algorithm to obtain an enhanced image to be detected;
and extracting a target ground area of the enhanced image to be detected to obtain the target remote sensing image.
Further, the training process of the saline-alkali degree classification network comprises the following steps:
constructing a saline-alkali degree classification network;
obtaining a sample saline-alkali image set, wherein a label corresponding to a sample saline-alkali image in the sample saline-alkali image set is the saline-alkali degree;
and training the saline-alkali degree classification network by utilizing the sample saline-alkali image set and the label corresponding to each sample saline-alkali image in the sample saline-alkali image set to obtain the trained saline-alkali degree classification network.
The invention has the following beneficial effects:
according to the saline-alkali geological identification method based on the unmanned aerial vehicle remote sensing image, the technical problem that the accuracy of saline-alkali identification of soil is low is solved by image processing of the remote sensing image to be detected, and the accuracy of saline-alkali identification of the soil is improved. Firstly, a target camera installed on a target unmanned aerial vehicle is used for acquiring a remote sensing image to be detected of a target ground area to be detected for saline-alkali conditions. The remote sensing image to be detected contains the information of the target ground area, so that the subsequent image processing of the remote sensing image to be detected can be facilitated, and the saline-alkali condition corresponding to the remote sensing image to be detected can be determined, so that the saline-alkali condition of the target ground area can be determined, and the accuracy of determining the saline-alkali condition of the target ground area can be improved. And secondly, carrying out image enhancement extraction processing on the remote sensing image to be detected to obtain a target remote sensing image. The image enhancement processing is carried out on the remote sensing image to be detected, the image contrast of the remote sensing image to be detected can be improved, the visualization degree of the remote sensing image to be detected is higher, and the subsequent image processing can be conveniently carried out on the remote sensing image to be detected. Secondly, in actual conditions, through the target camera, not only the target ground area but also an area other than the target ground area is often shot on the obtained remote sensing image to be detected. However, the regions other than the target ground region do not need to be subjected to saline-alkali condition judgment, so that the enhanced image to be detected is extracted to obtain the target remote sensing image only shooting the target ground region, subsequent steps can be omitted in the regions other than the target ground region, the calculation amount can be reduced, the occupation of calculation resources can be reduced, and the efficiency of saline-alkali identification on the target ground region can be improved. And then, carrying out color space conversion processing on the target remote sensing image to obtain a color remote sensing image. Since the color remote sensing image may be an HSV image. HSV images tend to have less correlation between the three channels of H (Hue), S (Saturation) and V (Value), and tend to be more visually consistent. Subsequent analysis of H, S and V three channels can be facilitated. And then, carrying out region division optimization on the color remote sensing image to obtain a primary optimal saline-alkali region set. And continuing to merge and divide the initial optimal saline-alkali regions in the initial optimal saline-alkali region set to obtain a target saline-alkali region set. In practical situations, the saline-alkali degree corresponding to each position in the target ground area is often more than one, that is, the target ground area may often include a plurality of areas with different saline-alkali degrees. Therefore, the color remote sensing images are subjected to area division, optimization and combination, a plurality of target saline-alkali areas with different corresponding saline-alkali degrees can be obtained, and the saline-alkali condition of the target ground area can be conveniently and accurately judged subsequently. And then, cutting each target saline-alkali area in the target saline-alkali area set, determining a target saline-alkali image corresponding to the target saline-alkali area, and obtaining a target saline-alkali image set. The different target saline-alkali regions of the corresponding saline-alkali conditions are cut off to obtain the target saline-alkali images in one-to-one correspondence with the saline-alkali degrees, so that the subsequent analysis of the multiple saline-alkali degrees corresponding to the target ground regions and the position of each saline-alkali degree in the target ground regions can be facilitated. Then, inputting each target saline-alkali image in the target saline-alkali image set into a trained saline-alkali degree classification network, and determining the saline-alkali degree corresponding to the target saline-alkali image through the saline-alkali degree classification network. And finally, generating target saline-alkali information representing the saline-alkali condition of the target ground area according to the saline-alkali degree corresponding to each target saline-alkali image in the target saline-alkali image set. Therefore, the method and the device solve the technical problem of low accuracy of saline-alkali identification of the soil by processing the image of the remote sensing image to be detected, and improve the accuracy of the saline-alkali identification of the soil.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions and advantages of the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and other drawings can be obtained by those skilled in the art without creative efforts.
Fig. 1 is a flow diagram of some embodiments of a method for saline-alkali geological identification based on unmanned aerial vehicle remote sensing images according to the invention;
FIG. 2 is a schematic diagram of an enhanced image to be detected and a target remote sensing image according to the present invention;
FIG. 3 is a schematic diagram of a distance between two target preliminary regions being zero according to the present invention;
fig. 4 is a schematic diagram of a target normalized image, a target saline-alkali region and a target saline-alkali image according to the invention.
Wherein the reference numerals in fig. 2 include: the method comprises the steps of an enhanced image to be detected 201, an image area 202 and a target remote sensing image 203.
The reference numerals in fig. 3 include: a first target preliminary region 301 and a second target preliminary region 302.
The reference numerals in fig. 4 include: the target normalized image 401, the first target saline-alkali region 402, the second target saline-alkali region 403, the first target saline-alkali image 404 and the second target saline-alkali image 405.
Detailed Description
To further explain the technical means and effects of the present invention adopted to achieve the predetermined objects, the following detailed description of the embodiments, structures, features and effects of the technical solutions according to the present invention will be given with reference to the accompanying drawings and preferred embodiments. In the following description, different references to "one embodiment" or "another embodiment" do not necessarily refer to the same embodiment. Furthermore, the particular features, structures, or characteristics may be combined in any suitable manner in one or more embodiments.
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.
The invention provides a saline-alkali geological identification method based on an unmanned aerial vehicle remote sensing image, which comprises the following steps:
acquiring a remote sensing image to be detected of a target ground area of a saline-alkali condition to be detected through a target camera installed on a target unmanned aerial vehicle;
carrying out image enhancement extraction processing on the remote sensing image to be detected to obtain a target remote sensing image;
carrying out color space conversion processing on the target remote sensing image to obtain a color remote sensing image;
carrying out region division optimization on the color remote sensing image to obtain a preliminary optimal saline-alkali region set;
merging and dividing the initial optimal saline-alkali regions in the initial optimal saline-alkali region set to obtain a target saline-alkali region set;
cutting each target saline-alkali area in the target saline-alkali area set, determining a target saline-alkali image corresponding to the target saline-alkali area, and obtaining a target saline-alkali image set;
inputting each target saline-alkali image in the target saline-alkali image set into a trained saline-alkali degree classification network, and determining the saline-alkali degree corresponding to the target saline-alkali image through the saline-alkali degree classification network;
according to the saline-alkali degree corresponding to each target saline-alkali image in the target saline-alkali image set, and generating target saline-alkali information representing the saline-alkali condition of the target ground area.
The following steps are detailed:
referring to fig. 1, a flow of some embodiments of a method for identifying saline-alkali geology based on unmanned aerial vehicle remote sensing images according to the present invention is shown. The saline-alkali geological identification method based on the unmanned aerial vehicle remote sensing image comprises the following steps:
s1, acquiring a remote sensing image to be detected of a target ground area to be detected for saline-alkali conditions through a target camera installed on a target unmanned aerial vehicle.
In some embodiments, the remote sensing image to be detected of the target ground area to be detected for saline-alkali condition can be acquired by a target camera installed on the target unmanned aerial vehicle.
Wherein, above-mentioned target unmanned aerial vehicle is the unmanned aerial vehicle that can normally fly. The target camera may be a remote sensing camera mounted on the target drone. The target ground area may be a ground area to be detected for saline-alkali conditions. The remote sensing image to be detected can be a remote sensing image of a target ground area to be detected for saline-alkali conditions.
As an example, the target unmanned aerial vehicle may be controlled to move to a target position, and the target camera is used to obtain a remote sensing image to be detected of a target ground area where saline-alkali conditions are to be detected. Wherein the target location may be a location directly above the target ground area. The distance between the target position and the target ground area may be a preset height, and may be set according to actual conditions. The shooting angle of the target camera may be a preset angle, and may be set according to actual conditions. The data set according to the actual situation needs to ensure that the remote sensing image to be detected can be acquired through the target camera.
And S2, carrying out image enhancement extraction processing on the remote sensing image to be detected to obtain a target remote sensing image.
In some embodiments, the remote sensing image to be detected may be subjected to image enhancement extraction processing to obtain a target remote sensing image.
The target remote sensing image can be a remote sensing image to be detected after image enhancement extraction processing.
As an example, this step may include the steps of:
firstly, carrying out image enhancement processing on the remote sensing image to be detected through an image enhancement algorithm to obtain an enhanced image to be detected.
The enhanced image to be detected can be a remote sensing image to be detected after image enhancement processing. Image enhancement algorithms may include, but are not limited to: histogram equalization algorithms, laplacian algorithms, and gamma transformation algorithms.
And secondly, extracting a target ground area of the enhanced image to be detected to obtain the target remote sensing image.
The target remote sensing image can be a region image corresponding to a target ground region.
As an example, as shown in fig. 2, a target ground area extraction process may be performed on the enhanced image to be detected 201 to obtain a target remote sensing image 203. The image area 202 may be an area where the target ground area corresponds to the enhanced image 201 to be detected. And cutting the image area 202 from the enhanced image 201 to be detected to obtain a target remote sensing image 203.
In practical situations, through the target camera, not only the target ground area but also an area other than the target ground area is often shot on the obtained remote sensing image to be detected. However, the regions other than the target ground region do not need to be subjected to saline-alkali condition judgment, so that the target ground region extraction processing is performed on the to-be-detected enhanced image to obtain a target remote sensing image only shooting the target ground region, the subsequent steps are not performed on the regions other than the target ground region, the calculation amount can be reduced, the occupation of calculation resources can be reduced, and the efficiency of saline-alkali identification on the target ground region can be improved.
And S3, performing color space conversion processing on the target remote sensing image to obtain a color remote sensing image.
In some embodiments, the target remote sensing image may be subjected to color space conversion processing to obtain a color remote sensing image.
The color remote sensing image can be an HSV image converted from a target remote sensing image.
As an example, the target remote sensing image may be converted into a color remote sensing image by a conventional technique.
And S4, carrying out region division optimization on the color remote sensing image to obtain a preliminary optimal saline-alkali region set.
In some embodiments, the color remote sensing image may be subjected to region division optimization to obtain a preliminary optimal saline-alkali region set.
As an example, this step may include the steps of:
and step one, mapping the color remote sensing image into a target undirected graph.
The target undirected graph can be an undirected graph mapped by the color remote sensing image. The nodes of the target undirected graph can be pixel points in the color remote sensing image. The edges of the target undirected graph can be connecting lines among pixel points in the color remote sensing image.
And secondly, determining the initial optimal saline-alkali area set through an optimization algorithm according to the target undirected graph.
Wherein, the saline and alkaline degree that preliminary optimum saline and alkaline area in the above-mentioned preliminary optimum saline and alkaline area set corresponds can be different. The degree of salt/alkali may be, but is not limited to: normal area, light saline-alkali area, moderate saline-alkali area and severe saline-alkali area.
For example, the objective function of the above optimization algorithm may be:
wherein,is the objective function of the optimization algorithm described above.FAre the objective function values of the above-described optimization algorithm.Is the first objective function.HIs the first objective function value.。So as to makeCounting at the bottomLogarithm. For example,。is the second in the above target undirected graphiA pixel point andfirst, thejThe transition probability between individual pixel points.。。Is the second in the above target undirected graphiAnd the sum of the edge weight value between each pixel point and each pixel point in the target undirected graph.IIs the number of pixel points in the target undirected graph. The second in the above objective undirected graphiPixel point and the secondjThe edge weight between each pixel point is。exp() Is an exponential function with a natural constant as the base.Is the second in the above target undirected graphiPixel point and the secondjThe Euclidean distance between pixel points.Is the second in the above target undirected graphiHue value (H, hue) corresponding to each pixel point and the secondjDifference values of hue values corresponding to the individual pixel points.Is the second in the above target undirected graphiSaturation value (S) corresponding to each pixel point and the second pixel pointjThe difference in saturation values corresponding to each pixel point.Is the second in the above target undirected graphiBrightness Value (V, value) corresponding to each pixel point and the second pixel pointjDifference in luminance values corresponding to each pixel point.UAre the model parameters.Are the model coefficients.Is the second objective function.ObjIs the second objective function value.Is the second in the above target undirected graphkThe number of pixels in each sub-region. The sub-region is a region in the above-described target undirected graph.Andare respectively the second in the above-mentioned target undirected graphkThe length and width of the minimum bounding rectangle corresponding to a sub-region.KIs the number of sub-regions obtained by dividing the target undirected graph. Wherein the model parametersUAnd coefficient of modelMay be preset. For example,Uand (2). Model coefficientsMay be positive real numbers.
When the value of the objective functionFAnd when the maximum value is obtained, dividing the obtained area to obtain a primary optimal saline-alkali area.
Due to the fact thatIs the second in the above target undirected graphiPixel point and the secondjThe transition probability between individual pixel points.。。Is the second in the above target undirected graphiAnd the sum of the edge weight value between each pixel point and each pixel point in the target undirected graph.IIs the number of pixel points in the target undirected graph. The second in the above objective undirected graphiPixel point and the secondjThe edge weight between each pixel point is. So that the first objective function valueHEntropy rate values corresponding to the target undirected graph can be characterized. The larger the entropy value is, the higher the structural similarity of each sub-region obtained after the color remote sensing image is divided is, that is, the greater the similarity between pixel points in the sub-region is, and the target function value isFThe larger the division result, the better the division result. Model coefficientsCan be set according to actual conditions, so that the target functionAdding model coefficientsThe result can be more in line with the actual situation. Due to the fact thatIs the second in the above-mentioned target undirected graphkThe number of pixels in each sub-region. The sub-region is a region in the above-described target undirected graph.Andare respectively the second in the above-mentioned target undirected graphkIndividual subareaThe field corresponds to the length and width of the minimum bounding rectangle. Therefore, it is not only easy to useCan characterize thekMorphological characteristics of the sub-regions. When the morphological characteristics of the sub-regions are larger, the total number of the divided sub-regionsKThe smaller the number of the division, the more stable the form of the divided sub-region, the more compact the structure in the sub-region, and the objective function valueFThe larger the division result, the better the division result.
And S5, merging and dividing the initial optimal saline-alkali regions in the initial optimal saline-alkali region set to obtain a target saline-alkali region set.
In some embodiments, the preliminary optimal saline-alkali regions in the preliminary optimal saline-alkali region set may be merged and divided to obtain a target saline-alkali region set.
The saline-alkali degree corresponding to each target saline-alkali area in the target saline-alkali area set can be different.
As an example, this step may include the steps of:
firstly, normalizing the target remote sensing image to obtain a target normalized image.
The target normalized image can be a normalized target remote sensing image.
And secondly, determining a target preliminary region corresponding to each preliminary optimal saline-alkali region in the preliminary optimal saline-alkali region set according to the position of each preliminary optimal saline-alkali region in the color remote sensing image, and obtaining a target preliminary region set.
The target preliminary region in the target preliminary region set may be a region in the target normalized image.
For example, the position of the preliminary optimal saline-alkali region in the color remote sensing image can be determined as the target position. And determining the region where the target position in the target normalized image is located as a target preliminary region corresponding to the preliminary optimal saline-alkali region. Wherein, the shape and size of the preliminary optimal saline-alkali region may be the same as the shape and size of the target preliminary region corresponding to the preliminary optimal saline-alkali region.
Because the color remote sensing image is the target remote sensing image after color space conversion, actual scenes corresponding to pixel points at the same position in the color remote sensing image and the target remote sensing image can be the same. The initial optimal saline-alkali region in the initial optimal saline-alkali region set can be a region in the color remote sensing image. The target preliminary region in the set of target preliminary regions may be a region in the target normalized image. The target normalized image may be a normalized target remote sensing image. Therefore, the actual scenes corresponding to the pixel points at the same position in the color remote sensing image, the target remote sensing image and the target normalized image can be the same. Therefore, the actual scenes corresponding to the image areas at the same position in the color remote sensing image, the target remote sensing image and the target normalized image can be the same. Therefore, the first position may be the same as the second position. The first location may be a location of the preliminary optimal saline-alkali region in the color remote sensing image. The second location may be a location of a preliminary region of interest in the normalized image of interest corresponding to the preliminary optimal saline-alkali region.
And thirdly, dividing the channel value of each channel in the target quantity of channels of the target normalized image into a preset quantity of channel grades to obtain a channel grade set.
The target number may be the number of channels of the target normalized image. For example, the channels of the target normalized image may be an R (Red), G (Green), and B (Blue) channel in RGB mode, respectively. The target number may be 3. Each channel of the target normalized image may correspond to a preset number of channel levels. The preset number may be a preset number. For example, the preset number may be 16. The number of channel levels in the channel level set may be raised to the power of a preset number of target numbers. For example, the number of channel levels in the set of channel levels described above may be 16 to the power of 3.
For example, the preset number may be 2. The channel values of the R channel of the target normalized image may include: 0.1,0.2,0.6 and 0.7. The channel values of the R channel of the target normalized image may be divided into 2 channel levels. Where 0.1 and 0.2 may be one channel level. 0.6 and 0.7 may be another channel rank.
And fourthly, determining a color channel histogram corresponding to each target preliminary region in the target preliminary region set.
The color channel histogram may be a histogram in which the channel level is used as a horizontal axis value and the number of target pixel points is used as a vertical axis value. The target pixel point can be a pixel point of which the pixel value in the target preliminary region belongs to the channel level.
And fifthly, determining a first region merging index between the two target preliminary regions according to the color channel histograms corresponding to each two target preliminary regions in the target preliminary region set and the channel level set.
For example, the above formula for determining the correspondence between the first region merging indicators of the two target preliminary regions may be:
wherein,is the first in the set of target preliminary regionspA target preliminary region andqand merging indexes of the first region between the target preliminary regions.CIs the number of channel levels in the set of channel levels described above.Is the first in the set of target preliminary regionspThe first one of the target preliminary regions in the above channel level set included in the corresponding color channel histogramcHistogram distribution values at each channel level.Is the first in the set of target preliminary regionsqThe first one of the target preliminary regions in the above channel level set included in the corresponding color channel histogramcHistogram distribution values at each channel level. The histogram distribution value may be a vertical axis value of the color channel histogram.Is a preset number greater than 0.May be a very small number.The effect of (3) is mainly to prevent the denominator from being 0.
Due to the fact thatIs the first in the set of target preliminary regionspThe first one of the target preliminary regions in the above channel level set included in the corresponding color channel histogramcHistogram distribution values at each channel level.Is the first in the set of target preliminary regionsqThe first one of the target preliminary regions in the above channel level set included in the corresponding color channel histogramcHistogram distribution values at each channel level. Therefore, the first and second electrodes are formed on the substrate,can characterize thepA target preliminary region andqthe difference in histogram distribution values of the c-th channel level of the respective target preliminary regions. Since only the difference between the two needs to be determined here, and the positive and negative need not be considered, the method has the advantages of simple structure, low cost and high reliabilityThe difference between the two can be characterized, anThe larger the difference between the two. Therefore, the first and second electrodes are formed on the substrate,the difference between the p-th target preliminary region and the q-th target preliminary region can be characterized and then squared to obtainThe difference between the p-th target preliminary region and the q-th target preliminary region may be made more conspicuous.The denominator can be prevented from being 0.The smaller the size of the tube is,the larger. Therefore it is firstpA target preliminary region andqfirst region merging indicator between target preliminary regionsThe larger, thepA target preliminary region andqthe more similar the individual target preliminary regions, the more likely they can be merged together.
And sixthly, determining a second region merging index between the two target preliminary regions according to every two target preliminary regions in the target preliminary region set.
For example, this step may include the following sub-steps:
and a first substep of translating the two target preliminary regions to make the distance between the two target preliminary regions zero, and obtaining fitting regions corresponding to the two target preliminary regions.
The fitting union region corresponding to the two target preliminary regions may be a region obtained by merging the two target preliminary regions.
For example, when the distance between the two target preliminary regions is greater than zero, the two target preliminary regions are translated to make the distance between the two target preliminary regions zero, and the fitting regions corresponding to the two target preliminary regions are obtained.
For another example, when the distance between the two target preliminary regions is zero, the two target preliminary regions do not need to be translated, and the two target preliminary regions can be directly merged together to obtain a fit region corresponding to the two target preliminary regions. As shown in fig. 3, the distance between the first target preliminary region 301 and the second target preliminary region 302 is zero. The first target preliminary region 301 and the second target preliminary region 302 may be two target preliminary regions.
And a second sub-step of determining a second region merging indicator between the two target preliminary regions according to the two target preliminary regions and the corresponding fitted regions of the two target preliminary regions.
For example, the formula for determining the correspondence between the second region merging indicators of the two target preliminary regions may be:
wherein,is the first in the set of target preliminary regionspA target preliminary region andqsecond region merging indicators between the target preliminary regions.exp() Is an exponential function with a natural constant as the base.Is the first in the set of target preliminary regionspA target preliminary region andqand the number of pixel points in the fitting and region corresponding to each target preliminary region.Is the first in the set of target preliminary regionspA target preliminary region andqand the number of edge pixel points included in the fitting and region corresponding to each target preliminary region.Is the first in the set of target preliminary regionspA target preliminary region andqand the perimeter of the minimum bounding rectangle corresponding to the fitting and region corresponding to each target preliminary region.Is the first in the set of target preliminary regionspThe number of pixel points in each target initial region.Is the first in the set of target preliminary regionspThe number of edge pixel points included in each target preliminary region.Is the first in the set of target preliminary regionspThe perimeter of the minimum bounding rectangle corresponding to each target preliminary region.Is the first in the set of target preliminary regionsqThe number of pixel points in each target initial region.Is the first in the set of target preliminary regionsqThe number of edge pixel points included in each target preliminary region.Is the first in the set of target preliminary regionsqThe perimeter of the minimum bounding rectangle corresponding to each target preliminary region.
The closer to 0 the value of (A) is, thepA target preliminary region andqthe more similar the respective target preliminary regions.Can realize the pairNormalization is performed to facilitate comparison, anThe larger the size of the tube, the larger the tube,the smaller, thepA target preliminary region andqsecond region merging indicator between target preliminary regionsThe larger, thepA target preliminary region andqthe more similar the individual target preliminary regions, the more likely they can be merged together.
And seventhly, determining an overall region merging index between every two target preliminary regions according to the first region merging index and the second region merging index between every two target preliminary regions in the target preliminary region set.
For example, the above formula for determining the correspondence between the two target preliminary regions and the overall region merging indicator may be:
wherein,is the first in the set of target preliminary regionspA target preliminary region andqintegral region merging finger between target preliminary regionsAnd (4) marking.exp() Is an exponential function with a natural constant as the base.Is the first in the set of target preliminary regionspA target preliminary region andqand merging indexes of the first region between the target preliminary regions.Is the first in the set of target preliminary regionspA target preliminary region andqand second region merging indexes between the target preliminary regions.
Due to the fact thatpA target preliminary region andqfirst region merging indicator between target preliminary regionsThe larger, thepA target preliminary region andqthe more similar the individual target preliminary regions, the more likely they can be merged together. First, thepA target preliminary region andqsecond region merging indicator between target preliminary regionsThe larger, thepA target preliminary region andqthe more similar the individual target preliminary regions, the more likely they can be merged together. Therefore it is firstpA target preliminary region andqthe larger the overall region merging index between the target preliminary regions is, thepA target preliminary region andqthe more similar the individual target preliminary regions, the more likely they can be merged together. And isRealize the pairNormalization is performed to facilitate comparison, anThe larger the size of the tube is,the larger. In addition, the overall region merging index integrates the first region merging index and the second region merging index, and the accuracy of determining the overall region merging index is improved.
And eighthly, dividing the two target preliminary regions into the same target preliminary region when the overall region merging index between the two target preliminary regions in the target preliminary region set is greater than a preset judgment threshold value.
The determination threshold may be a minimum overall region merging indicator that considers that the two target preliminary regions are not the same type of target preliminary region. For example, the decision threshold may be 0.75. The same kind of target preliminary region may include target preliminary regions having the same degree of salt and alkali.
And ninthly, determining the target preliminary region in the same target preliminary region as the target saline-alkali region in the target saline-alkali region set.
And S6, cutting each target saline-alkali area in the target saline-alkali area set, determining a target saline-alkali image corresponding to the target saline-alkali area, and obtaining a target saline-alkali image set.
In some embodiments, each target saline-alkali region in the target saline-alkali region set may be cut, and a target saline-alkali image corresponding to the target saline-alkali region is determined, so as to obtain a target saline-alkali image set.
As an example, as shown in fig. 4, a first target saline-alkali region 402 and a second target saline-alkali region 403 in a target normalized image 401 may be subjected to a cutting process, so as to obtain a first target saline-alkali image 404 and a second target saline-alkali image 405. Wherein the first target saline alkali area 402 and the second target saline alkali area 403 can be two target saline alkali areas. The first target saline-alkali image 404 and the second target saline-alkali image 405 may be two target saline-alkali images.
And S7, inputting each target saline-alkali image in the target saline-alkali image set into the trained saline-alkali degree classification network, and determining the saline-alkali degree corresponding to the target saline-alkali image through the saline-alkali degree classification network.
In some embodiments, each target saline-alkali image in the target saline-alkali image set may be input into a trained saline-alkali degree classification network, and the saline-alkali degree corresponding to the target saline-alkali image is determined through the saline-alkali degree classification network.
The saline-alkali degree classification network can be used for judging the saline-alkali degree corresponding to the target saline-alkali image. The saline-alkali degree classification network can be a classification recognition neural network.
As an example, the training process of the salinity classification network may include the following steps:
firstly, constructing a saline-alkali degree classification network.
This step can be implemented by the prior art, and is not described herein again.
And secondly, acquiring a sample saline-alkali image set.
The sample saline-alkali image in the sample saline-alkali image set can be a ground image with known saline-alkali degree. The label that sample saline-alkali image in above-mentioned sample saline-alkali image set corresponds can be saline-alkali degree.
And thirdly, training a saline-alkali degree classification network by using the sample saline-alkali image set and labels corresponding to the sample saline-alkali images in the sample saline-alkali image set to obtain the trained saline-alkali degree classification network.
Wherein, the loss function of the training saline-alkali degree classification network can be a cross entropy loss function.
This step can be implemented by the prior art, and is not described herein again.
And S8, generating target saline-alkali information representing the saline-alkali condition of the target ground area according to the saline-alkali degree corresponding to each target saline-alkali image in the target saline-alkali image set.
In some embodiments, the target saline-alkali information representing the saline-alkali condition of the target ground area may be generated according to the saline-alkali degree corresponding to each target saline-alkali image in the target saline-alkali image set.
As an example, the saline-alkali degrees corresponding to the target saline-alkali image in the target saline-alkali image set may be a moderate saline-alkali area and a severe saline-alkali area, respectively. The generated target saline-alkali information representing the saline-alkali condition of the target ground area can be 'the middle part of the target ground area is a severe saline-alkali area, the rest areas are moderate saline-alkali areas, and the saline-alkali degree is relatively severe'.
According to the saline-alkali geological identification method based on the unmanned aerial vehicle remote sensing image, the technical problem that the accuracy of saline-alkali identification of soil is low is solved by image processing of the remote sensing image to be detected, and the accuracy of saline-alkali identification of the soil is improved. Firstly, a target camera installed on a target unmanned aerial vehicle is used for acquiring a remote sensing image to be detected of a target ground area to be detected for saline-alkali conditions. The remote sensing image to be detected contains the information of the target ground area, so that the subsequent image processing of the remote sensing image to be detected can be facilitated, and the saline-alkali condition corresponding to the remote sensing image to be detected can be determined, so that the saline-alkali condition of the target ground area can be determined, and the accuracy of determining the saline-alkali condition of the target ground area can be improved. And secondly, carrying out image enhancement extraction processing on the remote sensing image to be detected to obtain a target remote sensing image. The image enhancement processing is carried out on the remote sensing image to be detected, the image contrast of the remote sensing image to be detected can be improved, the visualization degree of the remote sensing image to be detected is higher, and the subsequent image processing can be conveniently carried out on the remote sensing image to be detected. In practical situations, the remote sensing image to be detected obtained by the target camera often shoots not only a target ground area but also areas except the target ground area. However, the regions other than the target ground region do not need to be subjected to saline-alkali condition judgment, so that the enhanced image to be detected is extracted to obtain the target remote sensing image only shooting the target ground region, subsequent steps can be omitted in the regions other than the target ground region, the calculation amount can be reduced, the occupation of calculation resources can be reduced, and the efficiency of saline-alkali identification on the target ground region can be improved. And then, carrying out color space conversion processing on the target remote sensing image to obtain a color remote sensing image. Since the color remote sensing image may be an HSV image. HSV images tend to have less correlation between the three channels of H (Hue), S (Saturation) and V (Value), and tend to be more visually consistent. Subsequent analysis of H, S and V three channels can be facilitated. And then, carrying out region division optimization on the color remote sensing image to obtain a primary optimal saline-alkali region set. And continuing to merge and divide the initial optimal saline-alkali regions in the initial optimal saline-alkali region set to obtain a target saline-alkali region set. In practical situations, the saline-alkali degree corresponding to each position in the target ground area is often more than one, that is, the target ground area may often include a plurality of areas with different saline-alkali degrees. Therefore, the color remote sensing images are subjected to area division, optimization and combination, a plurality of target saline-alkali areas with different corresponding saline-alkali degrees can be obtained, and the saline-alkali condition of the target ground area can be conveniently and accurately judged subsequently. And then, cutting each target saline-alkali area in the target saline-alkali area set, determining a target saline-alkali image corresponding to the target saline-alkali area, and obtaining a target saline-alkali image set. The different target saline-alkali regions of the corresponding saline-alkali conditions are cut off to obtain the target saline-alkali images in one-to-one correspondence with the saline-alkali degrees, so that the subsequent analysis of the multiple saline-alkali degrees corresponding to the target ground regions and the position of each saline-alkali degree in the target ground regions can be facilitated. Then, each target saline-alkali image in the target saline-alkali image set is input to a trained saline-alkali degree classification network, and the saline-alkali degree corresponding to the target saline-alkali image is determined through the saline-alkali degree classification network. And finally, generating target saline-alkali information representing the saline-alkali condition of the target ground area according to the saline-alkali degree corresponding to each target saline-alkali image in the target saline-alkali image set. Therefore, the method and the device solve the technical problem of low accuracy of saline-alkali identification of the soil by image processing of the remote sensing image to be detected, and improve the accuracy of the saline-alkali identification of the soil.
The above embodiments are only used to illustrate the technical solutions of the present application, and not to limit the same; although the present application has been described in detail with reference to the foregoing embodiments, it should be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; such modifications and substitutions do not substantially depart from the spirit and scope of the embodiments of the present application and are intended to be included within the scope of the present application.
Claims (7)
1. A saline-alkali geological identification method based on unmanned aerial vehicle remote sensing images is characterized by comprising the following steps:
acquiring a remote sensing image to be detected of a target ground area of a saline-alkali condition to be detected through a target camera installed on a target unmanned aerial vehicle;
carrying out image enhancement extraction processing on the remote sensing image to be detected to obtain a target remote sensing image;
carrying out color space conversion processing on the target remote sensing image to obtain a color remote sensing image;
carrying out region division optimization on the color remote sensing image to obtain a preliminary optimal saline-alkali region set;
merging and dividing the initial optimal saline-alkali regions in the initial optimal saline-alkali region set to obtain a target saline-alkali region set;
cutting each target saline-alkali area in the target saline-alkali area set, and determining a target saline-alkali image corresponding to the target saline-alkali area to obtain a target saline-alkali image set;
inputting each target saline-alkali image in the target saline-alkali image set into a trained saline-alkali degree classification network, and determining the saline-alkali degree corresponding to the target saline-alkali image through the saline-alkali degree classification network;
generating target saline-alkali information representing the saline-alkali condition of the target ground area according to the saline-alkali degree corresponding to each target saline-alkali image in the target saline-alkali image set;
the method comprises the following steps of carrying out region division optimization on the color remote sensing image to obtain a primary optimal saline-alkali region set, wherein the method comprises the following steps:
mapping the color remote sensing image into a target undirected graph;
determining the initial optimal saline-alkali area set through an optimization algorithm according to the target undirected graph;
the objective function of the optimization algorithm is:
wherein,is the objective function of the optimization algorithm,Fis the value of the objective function of the optimization algorithm,is a first objective function of the first set of functions,His the value of the first objective function,,so as to makeCounting at the bottomThe number of the pairs is logarithmic,is the second in the target undirected graphiPixel point and the secondjThe probability of a transition between individual pixel points,,,is the second in the target undirected graphiThe sum of the edge weights between each pixel point and each pixel point in the target undirected graph,Iis the number of pixel points in the target undirected graphiPixel point and the secondjThe edge weight between each pixel point is,exp() Is an exponential function with a natural constant as the base,is the second in the target undirected graphiA pixel point and a secondjThe euclidean distance between the individual pixel points,is the second in the target undirected graphiHue value corresponding to each pixel point and the secondjThe difference in hue value corresponding to each pixel point,is the second in the target undirected graphiSaturation value corresponding to each pixel point and the secondjThe difference of the saturation values corresponding to each pixel point,is the second in the target undirected graphiBrightness value corresponding to each pixel point and the secondjThe difference in luminance values corresponding to each pixel point,Uis the parameter of the model and is,is the coefficient of the model that is,is the second objective function of the first function,Objis the value of the second objective function,is the second in the target undirected graphkThe number of pixel points within a sub-region, a sub-region being a region in the target undirected graph,andare respectively the second in the target undirected graphkThe length and width of the smallest circumscribed rectangle corresponding to a sub-region,Kthe number of sub-regions obtained by dividing the target undirected graph is obtained;
the preliminary optimal saline-alkali area in the preliminary optimal saline-alkali area set is merged and divided to obtain a target saline-alkali area set, and the method comprises the following steps:
normalizing the target remote sensing image to obtain a target normalized image;
determining a target preliminary region corresponding to each preliminary optimal saline-alkali region in the preliminary optimal saline-alkali region set according to the position of each preliminary optimal saline-alkali region in the color remote sensing image to obtain a target preliminary region set, wherein the target preliminary region in the target preliminary region set is a region in the target normalized image;
dividing a channel value of each channel in a target number of channels of the target normalized image into a preset number of channel levels to obtain a channel level set, wherein each channel of the target normalized image corresponds to the preset number of channel levels, and the number of the channel levels in the channel level set is the power of the preset number of target number;
determining a color channel histogram corresponding to each target preliminary region in the target preliminary region set;
determining a first region merging index between every two target preliminary regions according to the channel level set and color channel histograms corresponding to the two target preliminary regions in the target preliminary region set;
determining a second region merging index between every two target preliminary regions according to every two target preliminary regions in the target preliminary region set;
determining an overall region merging index between every two target preliminary regions according to a first region merging index and a second region merging index between every two target preliminary regions in the target preliminary region set;
when the overall area merging index between two target preliminary areas in the target preliminary area set is larger than a preset judgment threshold value, dividing the two target preliminary areas into the same type of target preliminary area;
and determining the target preliminary region in the same target preliminary region as a target saline-alkali region in the target saline-alkali region set.
2. The method for identifying the saline-alkali geology based on the unmanned aerial vehicle remote sensing image according to claim 1, wherein the formula for determining the correspondence of the first region merging index between the two target preliminary regions is as follows:
wherein,is the first in the target preliminary region setpA target preliminary region andqa first region merging indicator between the target preliminary regions,Cis the number of channel levels in the set of channel levels,is the first in the target preliminary region setpA first one of the channel level sets included in the corresponding color channel histogram for each target preliminary regioncThe histogram distribution values at the level of one channel,is the first in the target preliminary region setqA first one of the channel level sets included in the corresponding color channel histogram for each target preliminary regioncThe histogram distribution values at the level of one channel,is a preset number greater than 0.
3. The method for identifying the saline-alkali geology based on the unmanned aerial vehicle remote sensing image according to claim 1, wherein the determining a second region merging index between two target preliminary regions according to each two target preliminary regions in the target preliminary region set comprises:
translating the two target preliminary regions to enable the distance between the two target preliminary regions to be zero, and obtaining fitting regions corresponding to the two target preliminary regions;
and determining a second region merging index between the two target preliminary regions according to the two target preliminary regions and the corresponding fitting regions of the two target preliminary regions.
4. The method for identifying the saline-alkali geology based on the unmanned aerial vehicle remote sensing image according to claim 3, wherein the formula for determining the correspondence of the second region merging index between the two target preliminary regions is as follows:
wherein,is the first in the target preliminary region setpA target preliminary region andqsecond region merging indicators between the target preliminary regions,exp() Is an exponential function with a natural constant as the base,is the first in the target preliminary region setpA target preliminary region andqthe number of pixel points in the fitting and region corresponding to each target preliminary region,is the first in the target preliminary region setpA target preliminary region andqthe number of edge pixel points included in the fitting and region corresponding to each target preliminary region,is the first in the target preliminary region setpA target preliminary region andqthe perimeter of the minimum bounding rectangle corresponding to the fitting and region corresponding to each target preliminary region,is the first in the target preliminary region setpThe number of pixel points within the individual target initialization region,is the first in the target preliminary region setpThe number of edge pixels included in each target preliminary region,is the first in the target preliminary region setpThe circumference of the minimum bounding rectangle corresponding to each target preliminary regionThe length of the utility model is long,is the first in the target preliminary region setqThe number of pixel points within the individual target initialization region,is the first in the target preliminary region setqThe number of edge pixels included in each target preliminary region,is the first in the target preliminary region setqThe perimeter of the minimum bounding rectangle corresponding to each target preliminary region.
5. The method for identifying the saline-alkali geology based on the unmanned aerial vehicle remote sensing image according to claim 1, wherein the formula for determining the correspondence of the overall region merging index between the two target preliminary regions is as follows:
wherein,is the first in the target preliminary region setpA target preliminary region andqthe overall region merging indicators between the target preliminary regions,exp() Is an exponential function with a natural constant as the base,is the first in the target preliminary region setpA target preliminary region andqa first region merging indicator between the target preliminary regions,is the first in the target preliminary region setpA target preliminary region andqand second region merging indexes between the target preliminary regions.
6. The method for identifying the saline-alkali geology based on the unmanned aerial vehicle remote sensing image according to claim 1, wherein the step of performing image enhancement extraction processing on the remote sensing image to be detected to obtain a target remote sensing image comprises the following steps:
carrying out image enhancement processing on the remote sensing image to be detected through an image enhancement algorithm to obtain an enhanced image to be detected;
and extracting a target ground area of the enhanced image to be detected to obtain the target remote sensing image.
7. The method for identifying the saline-alkali geology based on the unmanned aerial vehicle remote sensing image according to claim 1, wherein the training process of the saline-alkali degree classification network comprises the following steps:
constructing a saline-alkali degree classification network;
obtaining a sample saline-alkali image set, wherein a label corresponding to a sample saline-alkali image in the sample saline-alkali image set is the saline-alkali degree;
and training the saline-alkali degree classification network by utilizing the sample saline-alkali image set and the label corresponding to each sample saline-alkali image in the sample saline-alkali image set to obtain the trained saline-alkali degree classification network.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211068025.6A CN115147746B (en) | 2022-09-02 | 2022-09-02 | Saline-alkali geological identification method based on unmanned aerial vehicle remote sensing image |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211068025.6A CN115147746B (en) | 2022-09-02 | 2022-09-02 | Saline-alkali geological identification method based on unmanned aerial vehicle remote sensing image |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115147746A CN115147746A (en) | 2022-10-04 |
CN115147746B true CN115147746B (en) | 2022-11-29 |
Family
ID=83415144
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211068025.6A Active CN115147746B (en) | 2022-09-02 | 2022-09-02 | Saline-alkali geological identification method based on unmanned aerial vehicle remote sensing image |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115147746B (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116310842B (en) * | 2023-05-15 | 2023-08-04 | 菏泽市国土综合整治服务中心 | Soil saline-alkali area identification and division method based on remote sensing image |
CN117630337A (en) * | 2024-01-04 | 2024-03-01 | 中国科学院华南植物园 | Coral sand saline-alkali monitoring system based on unmanned aerial vehicle |
CN117953400A (en) * | 2024-02-23 | 2024-04-30 | 广州市市政集团设计院有限公司 | Green micro-reconstruction system for old and old urban communities |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101614818A (en) * | 2009-07-09 | 2009-12-30 | 中国科学院遥感应用研究所 | A kind of radar remote sensing monitoring method of salting of soil |
CN107784657A (en) * | 2017-09-29 | 2018-03-09 | 西安因诺航空科技有限公司 | A kind of unmanned aerial vehicle remote sensing image partition method based on color space classification |
CN108256621A (en) * | 2016-12-28 | 2018-07-06 | 北京天诚同创电气有限公司 | Bee colony optimization method and device |
CN108680509A (en) * | 2018-08-17 | 2018-10-19 | 山东农业大学 | A kind of strand salt marsh area soil salt content evaluation method |
CN109447111A (en) * | 2018-09-20 | 2019-03-08 | 杭州师范大学 | A kind of remote sensing supervised classification method based on subclass training sample |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108647703B (en) * | 2018-04-19 | 2021-11-02 | 北京联合大学 | Saliency-based classification image library type judgment method |
CN114494081B (en) * | 2022-04-01 | 2022-07-05 | 武汉大学 | Unmanned aerial vehicle remote sensing mapping image enhancement method |
-
2022
- 2022-09-02 CN CN202211068025.6A patent/CN115147746B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101614818A (en) * | 2009-07-09 | 2009-12-30 | 中国科学院遥感应用研究所 | A kind of radar remote sensing monitoring method of salting of soil |
CN108256621A (en) * | 2016-12-28 | 2018-07-06 | 北京天诚同创电气有限公司 | Bee colony optimization method and device |
CN107784657A (en) * | 2017-09-29 | 2018-03-09 | 西安因诺航空科技有限公司 | A kind of unmanned aerial vehicle remote sensing image partition method based on color space classification |
CN108680509A (en) * | 2018-08-17 | 2018-10-19 | 山东农业大学 | A kind of strand salt marsh area soil salt content evaluation method |
CN109447111A (en) * | 2018-09-20 | 2019-03-08 | 杭州师范大学 | A kind of remote sensing supervised classification method based on subclass training sample |
Non-Patent Citations (3)
Title |
---|
Land Cover Classification by using of Multi-source Remote Sensing Image based on ELM;Qing Ding et al.;《2019 SAR in Big Data Era (BIGSARDATA)》;20191007;1-4 * |
基于图像分割和密度聚类的遥感动目标分块提取;林翊钧等;《北京航空航天大学学报》;20180828;第44卷(第12期);2510-2520 * |
基于改进的遗传算法的多目标优化问题研究;孔德剑;《计算机仿真》;20120215;第29卷(第2期);213-215 * |
Also Published As
Publication number | Publication date |
---|---|
CN115147746A (en) | 2022-10-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN115147746B (en) | Saline-alkali geological identification method based on unmanned aerial vehicle remote sensing image | |
Jia et al. | Detection and segmentation of overlapped fruits based on optimized mask R-CNN application in apple harvesting robot | |
CN109871902B (en) | SAR small sample identification method based on super-resolution countermeasure generation cascade network | |
CN107563413B (en) | Unmanned aerial vehicle aerial image farmland block object accurate extraction method | |
CN106909902B (en) | Remote sensing target detection method based on improved hierarchical significant model | |
CN110309781B (en) | House damage remote sensing identification method based on multi-scale spectrum texture self-adaptive fusion | |
CN109063754B (en) | Remote sensing image multi-feature joint classification method based on OpenStreetMap | |
CN111797712B (en) | Remote sensing image cloud and cloud shadow detection method based on multi-scale feature fusion network | |
CN103049763B (en) | Context-constraint-based target identification method | |
CN111738064B (en) | Haze concentration identification method for haze image | |
CN111027446B (en) | Coastline automatic extraction method of high-resolution image | |
Mai et al. | Semi-supervised fuzzy C-means clustering for change detection from multispectral satellite image | |
CN109886146B (en) | Flood information remote sensing intelligent acquisition method and device based on machine vision detection | |
CN110176005B (en) | Remote sensing image segmentation method based on normalized index and multi-scale model | |
CN109685045A (en) | A kind of Moving Targets Based on Video Streams tracking and system | |
CN107464249A (en) | A kind of sheep only contactless body measurement method | |
CN106960182A (en) | A kind of pedestrian integrated based on multiple features recognition methods again | |
CN115410096B (en) | Satellite remote sensing image multi-scale fusion change detection method, medium and electronic device | |
CN110889840A (en) | Effectiveness detection method of high-resolution 6 # remote sensing satellite data for ground object target | |
Zhang et al. | Superpixel generation for SAR imagery based on fast DBSCAN clustering with edge penalty | |
CN109447111A (en) | A kind of remote sensing supervised classification method based on subclass training sample | |
CN113435254A (en) | Sentinel second image-based farmland deep learning extraction method | |
Dong et al. | Cloud detection method for high resolution remote sensing imagery based on the spectrum and texture of superpixels | |
CN107103295A (en) | Optical remote sensing image cloud detection method of optic | |
CN110120030A (en) | Processing method, application, computer-readable medium and the disease occurring area measuring method of wheat diseases generation image |
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 |