CN116878400A - Inversion mountain soil layer thickness obtaining method and system - Google Patents
Inversion mountain soil layer thickness obtaining method and system Download PDFInfo
- Publication number
- CN116878400A CN116878400A CN202310775206.0A CN202310775206A CN116878400A CN 116878400 A CN116878400 A CN 116878400A CN 202310775206 A CN202310775206 A CN 202310775206A CN 116878400 A CN116878400 A CN 116878400A
- Authority
- CN
- China
- Prior art keywords
- data
- soil
- vegetation coverage
- thickness
- obtaining
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 239000002689 soil Substances 0.000 title claims abstract description 218
- 238000000034 method Methods 0.000 title claims abstract description 67
- 238000012545 processing Methods 0.000 claims abstract description 15
- 238000012216 screening Methods 0.000 claims abstract description 8
- 238000005070 sampling Methods 0.000 claims description 15
- 238000010586 diagram Methods 0.000 claims description 12
- 238000004590 computer program Methods 0.000 claims description 9
- 230000001186 cumulative effect Effects 0.000 claims description 5
- 238000011835 investigation Methods 0.000 abstract description 15
- 238000001514 detection method Methods 0.000 abstract description 10
- 239000000463 material Substances 0.000 description 9
- 238000004458 analytical method Methods 0.000 description 6
- 238000005527 soil sampling Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 230000003628 erosive effect Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 230000000704 physical effect Effects 0.000 description 3
- 238000012795 verification Methods 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 238000010276 construction Methods 0.000 description 2
- 230000035558 fertility Effects 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 238000013421 nuclear magnetic resonance imaging Methods 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000010291 electrical method Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000004162 soil erosion Methods 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
- 238000000179 transient infrared spectroscopy Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/02—Measuring arrangements characterised by the use of optical techniques for measuring length, width or thickness
- G01B11/06—Measuring arrangements characterised by the use of optical techniques for measuring length, width or thickness for measuring thickness ; e.g. of sheet material
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Image Processing (AREA)
Abstract
A method and a system for acquiring soil layer thickness of an inversion mountain area relate to the technical field of geological detection. The method solves the problems of large acquisition difficulty, large error and high cost caused by the combination of investigation and interpolation in the existing soil thickness acquisition method. The acquisition method comprises the following steps: acquiring a remote sensing image of a mountain area, and screening the remote sensing image to obtain high-precision remote sensing image data; interpreting the high-precision remote sensing image data to obtain a vegetation coverage map; extracting a vegetation coverage map to obtain vegetation coverage data; acquiring soil thickness data measured in field work, and processing the soil thickness data and vegetation coverage data to obtain relation data between the vegetation coverage data and the soil thickness; obtaining a mountain area soil thickness map according to the relation data of the vegetation coverage data and the soil thickness and the vegetation coverage map; and obtaining a soil layer thickness value according to the soil thickness graph. The invention is suitable for the technical field of soil layer thickness detection in mountain areas.
Description
Technical Field
The invention relates to the technical field of geological exploration.
Background
The soil thickness refers to the vertical depth from the soil matrix layer to the soil surface, generally tens of centimeters to two meters, is an important basic attribute of the soil, can directly reflect the development degree of the soil, influences the soil fertility condition, and influences the earth surface processes such as vegetation growth, surface hydrology, soil moisture and the like.
The soil layer thickness of the inversion mountain area is mastered, so that a basis can be provided for land utilization planning, and the soil layer thickness is an important mark for judging the fertility, the profile development degree and the erosion resistance of the soil. Meanwhile, the knowledge of the soil layer thickness of the mountain area can also provide basis for resource utilization, such as development of mountain area mineral resources, construction of mountain area water conservancy and electric power and the like.
Therefore, the soil thickness of the mountain area has important significance in mountain soil utilization and water and soil conservation measure planning and ecological construction.
Many students study and analyze the soil thickness distribution conditions of different areas, but the soil thickness distribution conditions are obtained by investigation based on a large number of soil sampling points, and in order to achieve accuracy, at least 1000 sampling points are arranged in an area with the area of 100 square kilometers, so that a large amount of manpower and material resources are consumed. For plain areas, the mode of combining investigation and interpolation can reflect the actual thickness of soil to a certain extent, but for mountain areas, the difficulty of obtaining the soil thickness is high, a large amount of manpower and material resources are spent, and the thickness of the soil is far from actual because of low accuracy of the investigation points and the interpolation.
The soil layer thickness of the existing mountain area is mostly obtained by adopting a soil-rock physical parameter inversion method, a topography elevation difference method, a nuclear magnetic resonance imaging method and a soil physical property measurement method; however, the inversion method of physical parameters of the earth and rock requires a large amount of seismic logging and electrical method data, and the required data are complex, so that the cost is high. The terrain elevation difference potential method is only suitable for inversion of soil layer thickness in a small range, and is only suitable for mountain areas with steeper terrains. The nuclear magnetic resonance imaging method has lower detection depth and resolution, is similar to seismic logging, and can only be suitable for inversion of the soil layer thickness of a relatively shallow layer. The soil physical property measurement method needs to carry out depth detection on the near surface, the soil physical property is greatly influenced by the environment, and the reliability of the obtained data is low.
In summary, it is particularly important and urgent to explore a soil thickness acquisition method for inverting mountain areas accurately and with low cost, and to provide references for ecological geological planning.
Disclosure of Invention
The invention provides a method for acquiring soil layer thickness of an inversion mountain area, which solves the problems of high acquisition difficulty, high error and high cost caused by the fact that the existing soil thickness acquisition method adopts a mode of combining investigation and interpolation.
In order to achieve the above object, the present invention provides the following solutions:
the invention provides a method for obtaining the soil layer thickness of an inversion mountain area, which comprises the following steps:
s1, acquiring a remote sensing image of a mountain area, and screening the remote sensing image to obtain high-precision remote sensing image data;
s2, interpreting the high-precision remote sensing image data to obtain a vegetation coverage map;
s3, extracting the vegetation coverage map to obtain vegetation coverage data;
s4, acquiring soil thickness data measured in field work, and processing the soil thickness data and the vegetation coverage data to obtain relation data between the vegetation coverage data and the soil thickness;
s5, obtaining a mountain area soil thickness map according to the relation data of the vegetation coverage data and the soil thickness and the vegetation coverage map;
s6, obtaining a soil layer thickness value according to the soil thickness map.
Further, in a preferred embodiment, the step S2 specifically includes:
s21, processing the NDVI value of any one pixel information in the linear mixed pixel model into green vegetation information and non-vegetation information to obtain new pixel information;
s22, obtaining a new linear mixed pixel model according to the new pixel information;
s23, adopting the new linear mixed pixel model to interpret the high-precision remote sensing image data, and obtaining a vegetation coverage map.
Further, there is a preferred embodiment, the new pixel information in the step S21 is expressed as:
NDVI＝NDVI veg +NDVI soil ；
wherein, NDVI veg NDVI is green vegetation information soil Is non-vegetation information.
Further, in a preferred embodiment, the vegetation coverage data in the step S23 is expressed as:
the NDVI veg Is 95%, the NDVI is obtained by soil The cumulative distribution frequency of (2) was 5%.
Further, in a preferred embodiment, the vegetation coverage data in the step S3 is normally distributed.
Further, in a preferred embodiment, the step S4 specifically includes:
soil thickness data measured by field 150 sampling points are obtained, and the soil thickness data and the vegetation coverage data are processed to obtain relation data of the vegetation coverage data and the soil thickness.
Further, in a preferred embodiment, the step S5 specifically includes:
s51, carrying out average treatment on vegetation coverage of the same soil thickness according to the relation data of the vegetation coverage data and the soil thickness to obtain average data;
s52, drawing a scatter diagram and a correlation curve according to the mean value data;
and S53, drawing a mountain area soil thickness map according to the vegetation coverage map, the scatter diagram and the correlation curve.
The invention also provides a system for obtaining the soil layer thickness of the inversion mountain area, which comprises the following steps:
the storage device is used for acquiring a remote sensing image of a mountain area, screening the remote sensing image and obtaining high-precision remote sensing image data;
the storage device is used for interpreting the high-precision remote sensing image data to obtain a vegetation coverage map;
the storage device is used for extracting the vegetation coverage map to obtain vegetation coverage data;
the storage device is used for acquiring soil thickness data measured in field work, processing the soil thickness data and the vegetation coverage data and obtaining relation data of the vegetation coverage data and the soil thickness;
the storage device is used for obtaining a mountain area soil thickness map according to the relation data of the vegetation coverage data and the soil thickness and the vegetation coverage map;
and the storage device is used for obtaining the soil layer thickness value according to the soil thickness map.
The invention also provides a computer readable storage medium, wherein the computer readable storage medium is stored with a computer program, and the computer program is executed by a processor to execute the method for acquiring the soil layer thickness of the inversion mountain area.
The invention also provides a computer device, which comprises a memory and a processor, wherein the memory stores a computer program, and when the processor runs the computer program stored in the memory, the processor executes the method for acquiring the soil layer thickness of the inversion mountain area.
The beneficial effects of the invention are as follows:
1. the soil layer thickness detected by the existing inversion mountain area soil layer thickness method is based on a large number of soil sampling points, and at least 1000 sampling points are arranged in an area with the area of 100 square kilometers for realizing accuracy, so that field investigation workload is greatly increased, and detection cost is increased. The invention provides a method for obtaining soil layer thickness in inversion mountain areas, which adopts remote sensing image data and processes the remote sensing image to obtain vegetation coverage data, and simultaneously combines a small number of field sampling points, namely, the soil thickness value of the whole area can be inverted by only setting 150 field sampling points in an area with the area of 100 square kilometers, thereby greatly saving field investigation workload and manpower and material resources.
2. The soil layer thickness detected by the existing inversion mountain area soil layer thickness method is affected by the environment, so that the accuracy of the acquired soil layer thickness is low. The invention provides a method for acquiring soil layer thickness of an inversion mountain area, which adopts interpretation processing to remote sensing data to acquire vegetation coverage data, and simultaneously can invert the soil thickness value of the whole area by combining a small number of field investigation points, so that the soil layer thickness of the inversion mountain area is accurately acquired without being limited by geographical environment, and the soil layer thickness obtained by the soil layer thickness method has the advantages of less than 0.5 m difference from the actual measured soil layer thickness and high accuracy through verification.
3. The invention provides a method for obtaining soil layer thickness of inversion mountain area, which adopts NDVI value with 95% accumulated distribution frequency veg Green vegetation information and NDVI value accumulationNDVI with 5% distribution frequency soil The non-vegetation information can avoid the influence of noise on the NDVI value analysis of the remote sensing image and improve the analysis accuracy.
4. The invention provides an inversion mountain soil layer thickness acquisition method, which can plan vegetation suitable for planting with different soil thicknesses according to the obtained soil layer thickness, can monitor the erosion degree of mountain soil and provides a basis for evaluating the soil erosion condition of mountain areas.
The invention is suitable for the technical field of soil layer thickness detection in mountain areas.
Drawings
FIG. 1 is a flow chart of a method for obtaining an inversion of soil layer thickness in mountainous areas according to an embodiment;
FIG. 2 is a graph of vegetation coverage data according to embodiments five and eleven;
FIG. 3 is a table of field vegetation coverage versus soil thickness as described in embodiments six and eleven;
FIG. 4 is a plot of soil thickness versus vegetation coverage for embodiments seven and eleven;
fig. 5 is a schematic view of soil thickness according to the seventh and eleventh embodiments.
Detailed Description
The following describes in further detail the embodiments of the present invention with reference to the drawings and examples. The following examples will assist those skilled in the art in further understanding the present invention, but are not intended to limit the invention in any way. It should be noted that variations and modifications can be made by those skilled in the art without departing from the spirit of the invention, which falls within the scope of the invention.
Referring to fig. 1, the present embodiment provides a method for obtaining a soil layer thickness of an inversion mountain area, where the method includes:
s1, acquiring a remote sensing image of a mountain area, and screening the remote sensing image to obtain high-precision remote sensing image data;
s2, interpreting the high-precision remote sensing image data to obtain a vegetation coverage map;
s3, extracting the vegetation coverage map to obtain vegetation coverage data;
s4, acquiring soil thickness data measured in field work, and processing the soil thickness data and the vegetation coverage data to obtain relation data between the vegetation coverage data and the soil thickness;
s5, obtaining a mountain area soil thickness map according to the relation data of the vegetation coverage data and the soil thickness and the vegetation coverage map;
s6, obtaining a soil layer thickness value according to the soil thickness map.
When the method is actually applied, the remote sensing image of a mountain area is firstly obtained, the remote sensing image is screened, winter remote sensing image data with cloud cover less than 5% and vegetation in a non-flourishing period are screened, the winter remote sensing image data at the moment can better reflect the interpretation condition of vegetation coverage after the winter remote sensing image data, and further high-precision remote sensing image data are obtained, so that the accuracy of the soil layer thickness of the final reverse performance is higher. Interpreting the high-precision remote sensing image data to obtain a vegetation coverage map; extracting the vegetation coverage map to obtain vegetation coverage data; then acquiring soil thickness data measured in field work, and processing the soil thickness data and the vegetation coverage data to obtain relation data between the vegetation coverage data and the soil thickness; only a small number of field acquisition points are needed, and a large number of data acquisition points are not needed; taking the vegetation coverage map as a base map, and combining the relation data of the vegetation coverage data and the soil thickness to obtain a mountain area soil thickness map; and finally inverting the soil layer thickness of the whole domain according to the mountain soil thickness map. The method solves the problems of large acquisition difficulty, large error and high cost caused by the combination of investigation and interpolation in the existing method for acquiring the soil thickness in the mountain area.
The soil layer thickness detected by the existing inversion mountain area soil layer thickness method is based on a large number of soil sampling points, and at least 1000 sampling points are arranged in an area with the area of 100 square kilometers for realizing accuracy, so that field investigation workload is greatly increased, and detection cost is increased. The invention provides a method for obtaining soil layer thickness in inversion mountain areas, which adopts remote sensing image data and processes the remote sensing image to obtain vegetation coverage data, and simultaneously combines a small number of field sampling points, namely, the soil thickness value of the whole area can be inverted by only setting 150 field sampling points in an area with the area of 100 square kilometers, thereby greatly saving field investigation workload and manpower and material resources.
The soil layer thickness detected by the existing inversion mountain area soil layer thickness method is affected by the environment, so that the accuracy of the acquired soil layer thickness is low. The embodiment provides a method for obtaining the soil layer thickness of an inversion mountain area, which adopts interpretation processing to remote sensing data to obtain vegetation coverage data, and simultaneously combines a small amount of field investigation points to invert the soil thickness value of the whole area, so that the soil layer thickness of the inversion mountain area can be accurately collected without being limited by geographical environment, and the soil layer thickness obtained by the soil layer thickness method in the embodiment has the difference from the actual measured soil layer thickness of less than 0.5 meter and high accuracy.
According to the method for obtaining the soil layer thickness of the inversion mountain area, vegetation suitable for planting with different soil thicknesses can be planned according to the obtained soil layer thickness, meanwhile, the erosion degree of the soil in the mountain area can be monitored, and a basis is provided for evaluating the water and soil loss condition of the mountain area.
In the second embodiment, the step S2 in the method for obtaining the soil layer thickness of the inversion mountain area according to the first embodiment is illustrated, where the step S2 specifically includes:
s21, processing the NDVI value of any one pixel information in the linear mixed pixel model into green vegetation information and non-vegetation information to obtain new pixel information;
s22, obtaining a new linear mixed pixel model according to the new pixel information;
s23, adopting the new linear mixed pixel model to interpret the high-precision remote sensing image data, and obtaining a vegetation coverage map.
In practical application, the embodiment adopts the pixel binary model in the linear mixed pixel model to interpret the high-precision remote sensing image data, and specifically comprises the following steps: firstly, the NDVI value of any one pixel information in the pixel bipartite model is processed into green vegetation information and non-vegetation information, and new pixel information is obtained; the new pixel information is put into the original pixel bipartite model to obtain a new linear mixed pixel model; and then adopting the new linear mixed pixel model to interpret the NDVI value of the vegetation coverage part and the NDVI value of the non-vegetation in the high-precision remote sensing image data to obtain a new NDVI value, and obtaining a vegetation coverage map according to the NDVI value of the vegetation coverage part, the NDVI value of the non-vegetation and the new NDVI value.
In the embodiment, the high-precision remote sensing image data is processed by adopting a linear mixed pixel model, so as to obtain a vegetation coverage map. And because any one of the pixel information in the pixel bipartite model and the remote sensing image data are paired, the accuracy of the obtained vegetation coverage map is higher.
In the third embodiment, the new pixel information in step S21 in the method for obtaining the soil layer thickness of the inversion mountain area according to the second embodiment is illustrated, where the new pixel information is expressed as:
NDVI＝NDVI veg +NDVI soil ；
wherein, NDVI veg NDVI is green vegetation information soil Is non-vegetation information.
In practical application, the information in any one of the pixels in the pixel bipartite model is composed of green vegetation information and non-vegetation information, and the NDVI value of each pixel can be regarded as a weighted average value of the NDVI value of the vegetation coverage part and the NDVI value of the non-vegetation, so that the accuracy of remote sensing data analyzed by the pixel bipartite model is higher.
In the fourth embodiment, the vegetation coverage data in step S23 in the method for obtaining the soil layer thickness of the inversion mountain area according to the first embodiment is illustrated, where the vegetation coverage data is expressed as:
the NDVI veg Is 95%, the NDVI is obtained by soil The cumulative distribution frequency of (2) was 5%.
In practical application of the present embodiment, FVC represents vegetation coverage, and the value range is generally [0,1 ]]。NDVI veg The NDVI value representing the pure vegetation cover pixels is close to 1.NDVI soil The NDVI value, representing the bare earth covered pixel, is close to 0. Due to noise influence, the NDVI value of the remote sensing image is too high or too low, and the NDVI with the integrated distribution frequency of 95% of the NDVI value is selected veg And an NDVI value cumulative distribution frequency of 5% soil The influence of noise on the NDVI value analysis of the remote sensing image can be avoided, and the analysis accuracy is improved.
The embodiment provides a method for obtaining the soil layer thickness of an inversion mountain area, which adopts NDVI with the accumulated distribution frequency of the NDVI value of 95 percent veg NDVI with 5% cumulative distribution frequency of green vegetation information and NDVI value soil The non-vegetation information can avoid the influence of noise on the NDVI value analysis of the remote sensing image and improve the analysis accuracy.
In the fifth embodiment, referring to fig. 2, in the present embodiment, the vegetation coverage data in step S3 in the method for obtaining the inversion of the soil layer thickness in the mountain area according to the first embodiment is illustrated, and the vegetation coverage data is normally distributed.
In practical application, the embodiment extracts the vegetation coverage map to obtain vegetation coverage data, as shown in fig. 2, so that the vegetation coverage data is normally distributed, and the normally distributed data is favorable for data analysis and has high accuracy.
Referring to fig. 3, the present embodiment is an example of step S4 in the method for obtaining the soil layer thickness of the inversion mountain area according to the first embodiment, where the step S4 specifically includes:
soil thickness data measured by field 150 sampling points are obtained, and the soil thickness data and the vegetation coverage data are processed to obtain relation data of the vegetation coverage data and the soil thickness.
In practical application, 150 field sampling points are uniformly distributed in an area with the area of 100 square kilometers, soil thickness data measured by the 150 sampling points are obtained, the soil thickness data and the vegetation coverage data are processed, and relationship data between the field vegetation coverage data and the soil thickness are obtained, as shown in fig. 3; compared with the existing method for acquiring the soil layer thickness in the mountain area, the method has the advantages of greatly saving field investigation workload and saving manpower and material resources.
The soil layer thickness detected by the existing inversion mountain area soil layer thickness method is based on a large number of soil sampling points, and at least 1000 sampling points are arranged in an area with the area of 100 square kilometers for realizing accuracy, so that field investigation workload is greatly increased, and detection cost is increased. The invention provides a method for obtaining soil layer thickness in inversion mountain areas, which adopts remote sensing image data and processes the remote sensing image to obtain vegetation coverage data, and simultaneously combines a small number of field sampling points, namely, the soil thickness value of the whole area can be inverted by only setting 150 field sampling points in an area with the area of 100 square kilometers, thereby greatly saving field investigation workload and manpower and material resources.
Referring to fig. 4, the present embodiment is an example of step S5 in the method for obtaining the soil layer thickness of the inversion mountain area according to the first embodiment, where the step S5 specifically includes:
s51, carrying out average treatment on vegetation coverage of the same soil thickness according to the relation data of the vegetation coverage data and the soil thickness to obtain average data;
s52, drawing a scatter diagram and a correlation curve according to the mean value data;
and S53, drawing a mountain area soil thickness map according to the vegetation coverage map, the scatter diagram and the correlation curve.
In practical application, the present embodiment processes the relationship data between the field vegetation coverage data and the soil thickness obtained in the first embodiment, specifically: and (5) taking an average value of the field vegetation coverage of the same soil thickness to obtain average value data. Drawing a scatter diagram and a correlation curve according to the mean value data, as shown in fig. 4; and drawing a mountain area soil thickness map according to the scatter diagram and the correlation curve, as shown in fig. 5. And obtaining the soil layer thickness according to the mountain area soil thickness map.
The embodiment provides a method for obtaining soil layer thickness of inversion mountain areas, which comprises the steps of processing the relation data between the field vegetation coverage data and the soil thickness obtained in the first embodiment to obtain average data, drawing a scatter diagram and a correlation curve according to the average data, and simultaneously combining the vegetation coverage diagram to obtain the soil thickness diagram, so that the detection of the soil thickness of the mountain areas can be realized without a large amount of manpower and material resources, and the detection accuracy is high.
The eighth embodiment provides a system for obtaining soil layer thickness of an inversion mountain area, wherein the system is as follows:
the storage device is used for acquiring a remote sensing image of a mountain area, screening the remote sensing image and obtaining high-precision remote sensing image data;
the storage device is used for interpreting the high-precision remote sensing image data to obtain a vegetation coverage map;
the storage device is used for extracting the vegetation coverage map to obtain vegetation coverage data;
the storage device is used for acquiring soil thickness data measured in field work, processing the soil thickness data and the vegetation coverage data and obtaining relation data of the vegetation coverage data and the soil thickness;
the storage device is used for obtaining a mountain area soil thickness map according to the relation data of the vegetation coverage data and the soil thickness and the vegetation coverage map;
and the storage device is used for obtaining the soil layer thickness value according to the soil thickness map.
In a third embodiment, a method for obtaining a soil layer thickness of an inversion mountain area according to any one of the first to third embodiments is provided.
The tenth embodiment provides a computer device, which includes a memory and a processor, where the memory stores a computer program, and when the processor runs the computer program stored in the memory, the processor executes an obtaining method for inverting a soil layer thickness of a mountain area according to any one of the above.
In an eleventh embodiment, the method for obtaining the soil layer thickness of the inversion mountain area according to the first to seventh embodiments is described in verification, and the southern town of the Changping area of Beijing city is used as the verification area, the southern town area is about 200 square kilometers, and the mountain area is about 140 square kilometers. And screening out winter remote sensing image data with cloud cover below 5% and vegetation non-flourishing period by adopting a 30 m-precision Landsat 8OLI_TIRS remote sensing image. Using ENVI software, adopting a vegetation coverage graph interpreted by a pixel binary model as a base graph, normally distributing the interpreted vegetation coverage data, combining 150 points of field work measurement soil thickness data, extracting soil thickness data corresponding to the 150 points, and processing the soil thickness data and the vegetation coverage data to obtain relation data of the vegetation coverage data and the soil thickness, wherein the figure 2 shows that; averaging vegetation coverage of the same soil thickness, wherein the result is shown in figure 3; according to the mean value data, a scatter diagram and a correlation curve are drawn, as shown in fig. 4, the correlation curve is linear positive correlation, the correlation coefficient is 0.9185, and the correlation fitting condition is good. And taking the vegetation coverage map as a base map, and combining the vegetation coverage data and the soil thickness relationship data to obtain a mountain area soil thickness map, as shown in fig. 5. According to the mountain area soil thickness map, inverting the soil layer thickness value, comparing the obtained soil layer thickness value with the soil layer thickness value actually collected in the field, finding that the error difference is less than 0.5 m, and having high accuracy.
In conclusion, the soil layer thickness obtained by the soil layer thickness method disclosed by the embodiment has a difference of less than 0.5 meter from the actually measured soil layer thickness, and the accuracy is high.
Furthermore, the terms "first," "second," and the like, are used for descriptive purposes only and are not to be construed as indicating or implying a relative importance or implicitly indicating the number of technical features indicated. Thus, a feature defining "a first" or "a second" may explicitly or implicitly include at least one such feature. In the description of the present invention, the meaning of "plurality" means at least two, for example, two, three, etc., unless specifically defined otherwise.
In the description of the present specification, a description referring to terms "one embodiment," "some embodiments," "examples," "specific examples," or "some examples," etc., means that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the present invention. In this specification, schematic representations of the above terms are not necessarily directed to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples. Furthermore, the different embodiments or examples described in this specification and the features of the different embodiments or examples may be combined and combined by those skilled in the art without contradiction.
The above description is only an example of the present invention and is not limited to the present invention, but various modifications and changes will be apparent to those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention should be included in the scope of the claims of the present invention.
Claims (10)
1. The method for obtaining the soil layer thickness of the inversion mountain area is characterized by comprising the following steps:
s1, acquiring a remote sensing image of a mountain area, and screening the remote sensing image to obtain high-precision remote sensing image data;
s2, interpreting the high-precision remote sensing image data to obtain a vegetation coverage map;
s3, extracting the vegetation coverage map to obtain vegetation coverage data;
s4, acquiring soil thickness data measured in field work, and processing the soil thickness data and the vegetation coverage data to obtain relation data between the vegetation coverage data and the soil thickness;
s5, obtaining a mountain area soil thickness map according to the relation data of the vegetation coverage data and the soil thickness and the vegetation coverage map;
s6, obtaining a soil layer thickness value according to the soil thickness map.
2. The method for obtaining the soil layer thickness of the inverted mountain area according to claim 1, wherein the step S2 is specifically:
s21, processing the NDVI value of any one pixel information in the linear mixed pixel model into green vegetation information and non-vegetation information to obtain new pixel information;
s22, obtaining a new linear mixed pixel model according to the new pixel information;
s23, adopting the new linear mixed pixel model to interpret the high-precision remote sensing image data, and obtaining a vegetation coverage map.
3. The method according to claim 2, wherein the new pixel information in the step S21 is represented as:
NDVI＝NDVI veg +NDVI soil ；
wherein, NDVI veg NDVI is green vegetation information soil Is non-vegetation information.
4. The method according to claim 1, wherein the vegetation coverage data in step S3 is represented as:
the NDVI veg Is 95%, the NDVI is obtained by soil The cumulative distribution frequency of (2) was 5%.
5. The method for obtaining the soil layer thickness of the inverted mountain area according to claim 1, wherein the vegetation coverage data in the step S3 is normally distributed.
6. The method for obtaining the soil layer thickness of the inverted mountain area according to claim 1, wherein the step S4 is specifically:
soil thickness data measured by field 150 sampling points are obtained, and the soil thickness data and the vegetation coverage data are processed to obtain relation data of the vegetation coverage data and the soil thickness.
7. The method for obtaining the soil layer thickness of the inverted mountain area according to claim 1, wherein the step S5 is specifically:
s51, carrying out average treatment on vegetation coverage of the same soil thickness according to the relation data of the vegetation coverage data and the soil thickness to obtain average data;
s52, drawing a scatter diagram and a correlation curve according to the mean value data;
and S53, drawing a mountain area soil thickness map according to the vegetation coverage map, the scatter diagram and the correlation curve.
8. An acquisition system for inverting the soil layer thickness of a mountain area, which is characterized in that:
the storage device is used for acquiring a remote sensing image of a mountain area, screening the remote sensing image and obtaining high-precision remote sensing image data;
the storage device is used for interpreting the high-precision remote sensing image data to obtain a vegetation coverage map;
the storage device is used for extracting the vegetation coverage map to obtain vegetation coverage data;
the storage device is used for acquiring soil thickness data measured in field work, processing the soil thickness data and the vegetation coverage data and obtaining relation data of the vegetation coverage data and the soil thickness;
the storage device is used for obtaining a mountain area soil thickness map according to the relation data of the vegetation coverage data and the soil thickness and the vegetation coverage map;
and the storage device is used for obtaining the soil layer thickness value according to the soil thickness map.
9. A computer readable storage medium, characterized in that the computer readable storage medium has stored thereon a computer program which, when executed by a processor, performs a method of obtaining an inverted mountain soil layer thickness as claimed in any one of claims 1-7.
10. A computer device, characterized by: the apparatus comprises a memory and a processor, wherein the memory stores a computer program, and when the processor runs the computer program stored in the memory, the processor executes the method for acquiring the soil layer thickness of the inversion mountain area according to any one of claims 1-7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310775206.0A CN116878400A (en) | 2023-06-28 | 2023-06-28 | Inversion mountain soil layer thickness obtaining method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310775206.0A CN116878400A (en) | 2023-06-28 | 2023-06-28 | Inversion mountain soil layer thickness obtaining method and system |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN116878400A true CN116878400A (en) | 2023-10-13 |
Family
ID=88254008
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310775206.0A Pending CN116878400A (en) | 2023-06-28 | 2023-06-28 | Inversion mountain soil layer thickness obtaining method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116878400A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118670277A (en) * | 2024-05-30 | 2024-09-20 | 中国科学院南京土壤研究所 | Method, device, medium and product for identifying thickness of thick-layer black soil |
-
2023
- 2023-06-28 CN CN202310775206.0A patent/CN116878400A/en active Pending
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118670277A (en) * | 2024-05-30 | 2024-09-20 | 中国科学院南京土壤研究所 | Method, device, medium and product for identifying thickness of thick-layer black soil |
| CN118670277B (en) * | 2024-05-30 | 2024-12-13 | 中国科学院南京土壤研究所 | A method, device, medium and product for identifying the thickness of thick black soil |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Watlet et al. | Imaging groundwater infiltration dynamics in the karst vadose zone with long-term ERT monitoring | |
| Watkins et al. | Appraisal of fracture sampling methods and a new workflow to characterise heterogeneous fracture networks at outcrop | |
| Goodale et al. | Mapping monthly precipitation, temperature, and solar radiation for Ireland with polynomial regression and a digital elevation model | |
| Babaee et al. | Land subsidence from interferometric SAR and groundwater patterns in the Qazvin plain, Iran | |
| CN103234920B (en) | Based on the underground water enriching appraisal procedure of sensor information | |
| CN103810376A (en) | Ground daily rainfall predicting method based on satellite remote sensing and regression Kriging | |
| Herring et al. | Best practices for using electrical resistivity tomography to investigate permafrost | |
| CN116384717B (en) | Land utilization planning investigation project mapping data acquisition and analysis system | |
| CN113591288A (en) | Soil humidity data prediction method and device based on kriging interpolation | |
| Yu et al. | Spatial and temporal scale effect in simulating hydrologic processes in a watershed | |
| Huang et al. | An error budget for soil salinity mapping using different ancillary data | |
| Cui et al. | A two-step fusion framework for quality improvement of a remotely sensed soil moisture product: A case study for the ECV product over the Tibetan Plateau | |
| CN114970934A (en) | A Soil Thickness Type Prediction Method Based on Feature Ensemble Learning | |
| CN113624716A (en) | Soil nitrogen estimation method based on vegetation coverage | |
| CN105528523A (en) | Soil thickness inversion method based on remote sensing data | |
| Schneider et al. | Inverse modelling with a genetic algorithm to derive hydraulic properties of a multi-layered forest soil | |
| CN116878400A (en) | Inversion mountain soil layer thickness obtaining method and system | |
| Nazari et al. | Estimation of groundwater storage change in the Helmand River Basin (Afghanistan) using GRACE satellite data | |
| Coopersmith et al. | Using similarity of soil texture and hydroclimate to enhance soil moisture estimation | |
| CN109033459B (en) | A spatial grid construction method for soil bulk density data | |
| Yu et al. | Spatiotemporal sampling strategy for characterization of hydraulic properties in heterogeneous soils | |
| CN116630676A (en) | Large-scale-range field classification processing method and device and electronic equipment | |
| Guo-an et al. | Evaluation on the accuracy of digital elevation models | |
| Illawathure et al. | Evaluating soil moisture estimation from ground‐penetrating radar hyperbola fitting with respect to a systematic time‐domain reflectometry data collection in a boreal podzolic agricultural field | |
| CN118858583A (en) | A high-resolution soil scanning method and soil scanner |
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 |