Remote sensing estimation method for felling time point, growth duration and aboveground biomass of timber bamboo forest

Abstract

The invention provides a remote sensing estimation method and device for felling time point, growth duration and aboveground biomass of bamboo forest, electronic equipment and a storage medium, and belongs to the technical field of forestry. The method comprises the following steps: respectively performing patch segmentation on microwave remote sensing data and optical remote sensing data of a monitoring area; correspondingly obtaining a plurality of microwave plaques and a plurality of optical plaques; respectively calculating pixel mean values of the plurality of microwave plaques and the plurality of optical plaques, and correspondingly obtaining a microwave plaque mean value and an optical plaque mean value; based on the microwave plaque mean value and the optical plaque mean value, determining a wood bamboo forest felling area; based on the average value of the optical plaques, determining a time point of felling the wood bamboo forest; determining the growth duration of the wood bamboo forest; and determining the aboveground biomass based on the vegetation coverage of the timber bamboo forest and the growth duration of the timber bamboo forest. According to the method, the difference between the microwave remote sensing data and the optical remote sensing data is combined, so that the felling time point, the growth duration and the aboveground biomass estimation precision of the bamboo forest are improved.

Description

Technical Field

[0001] The invention relates to the technical field of forestry, and in particular to a remote sensing estimation method for the felling time point, growth duration and above-ground biomass of a timber bamboo forest. Background Art

[0002] Bamboo forests are a material used in industry and agriculture, and they grow fast and have strong carbon sequestration capabilities, which are of great significance in the context of "dual carbon". However, due to the lack of unified planning and management of bamboo forests, it is difficult to determine the felling area of ​​timber bamboo forests.

[0003] At present, remote sensing technology is usually used to identify different bamboo felling areas. Specifically, the bamboo felling area is obtained by land use classification or visual interpretation through spectral information or combined with texture information. In many bamboo forest growing areas, the topography and landforms are quite different. For example, the influence of shadows will lead to the similarity of the spectra of other forests and grasses with bamboo forests, resulting in uncertainty in the extracted information of bamboo forests, grasses and their mixtures, resulting in low accuracy of remote sensing estimation of the felling time, growth duration and aboveground biomass of timber bamboo forests.

[0004] Therefore, how to improve the accuracy of remote sensing estimation of the felling time, growth duration and aboveground biomass of timber bamboo forests has become a technical problem that needs to be solved urgently. Summary of the invention

[0005] The present invention provides a remote sensing estimation method for the felling time point, growth duration and aboveground biomass of a timber bamboo forest, so as to solve the defect of low accuracy of remote sensing estimation of the felling time point, growth duration and aboveground biomass of a timber bamboo forest in the prior art, and improve the accuracy of remote sensing estimation of the felling time point, growth duration and aboveground biomass of a timber bamboo forest.

[0006] The present invention provides a remote sensing estimation method for felling time, growth duration and aboveground biomass of a timber bamboo forest, comprising the following steps: Performing patch segmentation on the microwave remote sensing data and the optical remote sensing data of the monitoring area respectively, and obtaining a plurality of microwave patches and a plurality of optical patches accordingly; Calculating pixel means of the plurality of microwave patches and the plurality of optical patches respectively, and obtaining microwave patch means and optical patch means correspondingly; Determining the felling area of ​​the timber bamboo forest based on the microwave patch mean and the optical patch mean; Determining the time point of felling the timber bamboo forest based on the optical patch mean value; Determine the growth period of timber bamboo forest; The aboveground biomass is determined based on the vegetation coverage of the timber bamboo forest and the growth duration of the timber bamboo forest.

[0007] According to a remote sensing estimation method for felling time point, growth duration and aboveground biomass of a timber bamboo forest provided by the present invention, the method of determining the felling area of ​​the timber bamboo forest based on the microwave patch mean and the optical patch mean comprises: Determining a timber forest felling area identification index based on the microwave patch mean and the optical patch mean; Determining the timber bamboo forest felling area based on the timber forest felling area identification index; The timber forest felling area identification index is calculated based on the following formula: ; in, represents the timber forest felling area identification index, express Year The average optical patch value in 2018 is express The mean of the annual optical patch mean, express Year The average microwave patch value in 2018 is express The mean of the annual microwave patch means.

[0008] According to a remote sensing estimation method for felling time point, growth duration and aboveground biomass of a timber bamboo forest provided by the present invention, the method of determining the felling area of ​​the timber bamboo forest based on the timber forest felling area identification index comprises: If the timber forest felling area identification index is greater than a preset threshold, the monitoring area is determined to be a timber bamboo forest felling area; If the timber forest felling area identification index is less than or equal to the preset threshold, it is determined that the monitoring area is not a timber bamboo forest felling area.

[0009] According to a remote sensing estimation method for felling time, growth duration and aboveground biomass of a timber bamboo forest provided by the present invention, the pixel mean of the plurality of microwave patches is calculated based on the following steps: Calculating the plaque mean of each microwave plaque respectively; Determining a pixel mean based on the plaque mean of each of the microwave plaques; Among the plurality of microwave plaques, The mean plaque value of each microwave plaque is calculated based on the following formula: ; in, Indicates The mean value of the microwave plaques, represents the product of the image element values ​​of the microwave, represents the pixel value of VV polarization, Pixel value representing VH polarization.

[0010] According to a remote sensing estimation method for felling time of a timber bamboo forest, growth duration and aboveground biomass provided by the present invention, the time of felling a timber bamboo forest is determined based on the optical patch mean, comprising: Calculate the absolute difference of the annual optical patch mean values ​​of any two adjacent years within m years, where the absolute difference of the annual optical patch mean value is the absolute value of the difference of the optical patch mean values ​​of the two years, and m is a positive integer; The year corresponding to the maximum value of the annual optical patch mean absolute difference is taken as the time point for felling the timber bamboo forest, and the maximum value of the annual optical patch mean absolute difference is the maximum value among multiple annual optical patch mean absolute differences.

[0011] According to a remote sensing estimation method for felling time of a timber bamboo forest, growth duration and aboveground biomass provided by the present invention, the time of felling a timber bamboo forest is determined based on the optical patch mean, comprising: Calculate the absolute difference of the monthly optical patch mean values ​​between any two adjacent months within month o, where the absolute difference of the monthly optical patch mean values ​​is the absolute value of the difference between the optical patch mean values ​​of the two months, and o is a positive integer; The month corresponding to the maximum value of the absolute difference of the monthly optical patch means is taken as the time point for felling the timber bamboo forest, and the maximum value of the absolute difference of the monthly optical patch means is the maximum value among the multiple absolute differences of the monthly optical patch means.

[0012] According to a remote sensing estimation method for felling time, growth duration and aboveground biomass of a timber bamboo forest provided by the present invention, the vegetation coverage is calculated based on the following formula: ; in, Indicates the vegetation coverage of timber bamboo forest, Represents the Normalized Difference Vegetation Index NDVI for a specific pixel, represents the minimum value of NDVI in the monitoring area, Indicates the maximum value of NDVI in the monitoring area.

[0013] The present invention also provides a remote sensing estimation device for the felling time point, growth duration and aboveground biomass of a timber bamboo forest, comprising the following modules: The patch segmentation module is used to perform patch segmentation on the microwave remote sensing data and the optical remote sensing data of the monitoring area, and obtain a plurality of microwave patches and a plurality of optical patches accordingly; A mean value calculation module, used to calculate the pixel means of the plurality of microwave patches and the plurality of optical patches respectively, and obtain the microwave patch mean value and the optical patch mean value accordingly; A felling area determination module, used for determining a felling area of ​​a timber bamboo forest based on the microwave patch mean and the optical patch mean; A time point determination module, used to determine the time point of felling the timber bamboo forest based on the optical patch mean value; A growth duration determination module is used to determine the growth duration of the timber bamboo forest; The biomass determination module is used to determine the aboveground biomass based on the vegetation coverage of the timber bamboo forest and the growth time of the timber bamboo forest.

[0014] The present invention also provides an electronic device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the computer program, a remote sensing estimation method for the felling time point, growth duration, and aboveground biomass of a timber bamboo forest as described in any one of the above-mentioned methods is implemented.

[0015] The present invention also provides a non-transitory computer-readable storage medium having a computer program stored thereon. When the computer program is executed by a processor, the remote sensing estimation method for the felling time point, growth duration and aboveground biomass of a timber bamboo forest as described in any one of the above methods is implemented.

[0016] The present invention also provides a computer program product, including a computer program, which, when executed by a processor, implements any of the above-mentioned remote sensing estimation methods for the felling time point, growth duration and above-ground biomass of a timber bamboo forest.

[0017] The remote sensing estimation method for the felling time point, growth duration and aboveground biomass of a timber bamboo forest provided by the present invention performs patch segmentation on the microwave remote sensing data and the optical remote sensing data of the monitoring area respectively, and obtains a plurality of microwave patches and a plurality of optical patches accordingly; calculates the pixel means of the plurality of microwave patches and the plurality of optical patches respectively, and obtains the microwave patch mean and the optical patch mean accordingly; determines the felling area of ​​the timber bamboo forest based on the microwave patch mean and the optical patch mean; determines the felling time point of the timber bamboo forest based on the optical patch mean; and determines the aboveground biomass based on the vegetation coverage of the timber bamboo forest. Optical remote sensing can invert the vegetation structure of the bamboo forest, such as leaf area index, vegetation coverage, etc., and microwave remote sensing can identify the geometric information of the bamboo forest, and comprehensively reflect the surface undulation, soil moisture, geometric information, etc. through the backscattering coefficient. This scheme combines the differences between microwave remote sensing data and optical remote sensing data to improve the remote sensing estimation accuracy of the felling time point, growth duration and aboveground biomass of the timber bamboo forest. BRIEF DESCRIPTION OF THE DRAWINGS

[0018] In order to more clearly illustrate the technical solutions in the present invention or the prior art, the following briefly introduces the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0019] Figure 1 It is a flow chart of a remote sensing estimation method for felling time point, growth duration and above-ground biomass of a timber bamboo forest provided by the present invention.

[0020] Figure 2 It is a structural schematic diagram of a remote sensing estimation device for felling time point, growth duration and above-ground biomass of a timber bamboo forest provided by the present invention.

[0021] Figure 3 It is a structural schematic diagram of the electronic device provided by the present invention. DETAILED DESCRIPTION

[0022] In order to make the purpose, technical solution and advantages of the present invention clearer, the technical solution of the present invention will be clearly and completely described below in conjunction with the drawings of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0023] It should be noted that in the description of the embodiments of the present invention, the terms "include", "comprise" or any other variants thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or device. In the absence of further restrictions, the elements defined by the sentence "include one..." do not exclude the existence of other identical elements in the process, method, article or device including the elements. The orientation or position relationship indicated by the terms "upper", "lower" and the like is based on the orientation or position relationship shown in the drawings, which is only for the convenience of describing the present invention and simplifying the description, rather than indicating or implying that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and therefore cannot be understood as a limitation on the present invention. Unless otherwise clearly specified and limited, the terms "installed", "connected" and "connected" should be understood in a broad sense, for example, it can be a fixed connection, a detachable connection, or an integral connection; it can be a mechanical connection or an electrical connection; it can be directly connected, or indirectly connected through an intermediate medium, or it can be a connection between two elements. For those skilled in the art, the specific meanings of the above terms in the present invention can be understood according to specific circumstances.

[0024] The terms "first", "second", etc. in the present invention are used to distinguish similar objects, and are not used to describe a specific order or sequence. It should be understood that the data used in this way can be interchanged under appropriate circumstances, so that the embodiments of the present invention can be implemented in an order other than those illustrated or described herein, and the objects distinguished by "first", "second", etc. are generally of the same type, and the number of objects is not limited. For example, the first object can be one or more. In addition, "and / or" means at least one of the connected objects, and the character " / " generally indicates that the objects associated with each other are in an "or" relationship.

[0025] Bamboo forests are a material used in industry and agriculture, and they grow fast and have strong carbon sequestration capacity, which is of great significance in the context of "dual carbon". However, due to the lack of unified planning and management of bamboo forests, it is difficult to determine the felling time and carbon sequestration estimation of bamboo forests in different locations.

[0026] Remote sensing technology has the characteristics of multi-time, multi-space and multi-spectral, and can effectively obtain relevant information on surface cover and changes. At present, the application of remote sensing technology to extract information such as different bamboo forest growth time and felling time points and aboveground biomass estimation is mostly obtained through land use classification or visual interpretation through spectral information or combined with texture information. The problem with this is that the representation of bamboo forests in optical remote sensing images is prone to "same object different spectrum" and "same spectrum different object" phenomena, resulting in large differences between information such as bamboo forest growth time and felling and the actual situation, and there is great uncertainty in the results of information extraction such as bamboo forest growth time; second, for many bamboo forest growth areas, the terrain and landforms are quite different, such as shadows, which will lead to the similarity of other forest and grass spectra with bamboo forests, which may lead to uncertainty in the extraction of bamboo forests, grass and their mixed information. In the actual application process, it is necessary to extract the spatiotemporal distribution information of the widespread forest and grass mixture, and conventional classification or extraction methods are difficult to effectively obtain relevant information such as bamboo forests. Since bamboo forests mostly grow in the south, different spectra of the same object and different objects with the same spectrum are more obvious. In addition, since it is cloudy and rainy in the south, the quality of optical data has great problems.

[0027] In view of the above problems, the present invention proposes the following embodiments. Figure 1-Figure 3 The invention describes a remote sensing estimation method for felling time, growth duration and aboveground biomass of a timber bamboo forest provided by an embodiment of the present invention.

[0028] Figure 1 : is a flow chart of the remote sensing estimation method for the felling time point, growth duration and aboveground biomass of a timber bamboo forest provided by the present invention, such as Figure 1 As shown, the method includes the following: S110, performing patch segmentation on the microwave remote sensing data and the optical remote sensing data of the monitoring area respectively, and correspondingly obtaining a plurality of microwave patches and a plurality of optical patches; S120, respectively calculating pixel means of the plurality of microwave patches and the plurality of optical patches, and obtaining a microwave patch mean and an optical patch mean accordingly; S130, determining a felling area of ​​the timber bamboo forest based on the microwave patch mean and the optical patch mean; S140, determining a time point for felling the timber bamboo forest based on the optical patch mean value; S150, determine the growth period of the timber bamboo forest; S160: Determine the aboveground biomass based on the vegetation coverage of the timber bamboo forest and the growth duration of the timber bamboo forest.

[0029] It should be noted that the execution subject of the remote sensing estimation method for the felling time, growth duration and aboveground biomass of timber bamboo forest provided in the embodiment of the present invention can be a server, computer equipment, such as a mobile phone, tablet computer, laptop computer, PDA, ultra-mobile personal computer (UMPC), netbook, etc.

[0030] In S110, the monitoring area is the area where the timber bamboo forest monitoring is required, and during the specific implementation process, it can be set according to the actual situation.

[0031] In S110, the microwave remote sensing data and optical remote sensing data of the monitoring area are high-resolution microwave remote sensing data and optical remote sensing data of the time series covering the current year and previous years of the monitoring area, such as Landsat data, Sentinel data, high-resolution series data, etc.

[0032] Optionally, in S110, segmentation is performed by integrating spectrum and texture, and vector patches are obtained after segmentation.

[0033] In S120, the pixel averages of the plurality of microwave patches and the plurality of optical patches are calculated by the following formula: ; in, represents the microwave patch mean / optical patch mean, Indicates The average value of microwave patches / optical patches, The number of plaques representing microwave plaques / optical plaques.

[0034] It should be noted that there is no strict timing restriction between S130, S140 and S160, that is, some of the steps can be executed first, and then other steps, or S130, S140 and S150 can be executed at the same time; in addition, it can be understood that S150 must be executed before executing S160.

[0035] The remote sensing estimation method for the felling time point, growth duration and aboveground biomass of a timber bamboo forest provided by the embodiment of the present invention performs patch segmentation on the microwave remote sensing data and the optical remote sensing data of the monitoring area respectively, and obtains a plurality of microwave patches and a plurality of optical patches accordingly; calculates the pixel means of the plurality of microwave patches and the plurality of optical patches respectively, and obtains the microwave patch mean and the optical patch mean accordingly; determines the felling area of ​​the timber bamboo forest based on the microwave patch mean and the optical patch mean; determines the felling time point of the timber bamboo forest based on the optical patch mean; and determines the aboveground biomass based on the vegetation coverage of the timber bamboo forest. Optical remote sensing can invert the vegetation structure of the bamboo forest, such as leaf area index, vegetation coverage, etc., and microwave remote sensing can identify the geometric information of the bamboo forest, and comprehensively reflect the surface undulation, soil moisture, geometric information, etc. through the backscattering coefficient. This scheme combines the differences between microwave remote sensing data and optical remote sensing data to improve the remote sensing estimation accuracy of the felling time point, growth duration and aboveground biomass of the timber bamboo forest.

[0036] In an optional embodiment, the determining of the timber bamboo forest felling area based on the microwave patch mean and the optical patch mean includes: Determining a timber forest felling area identification index based on the microwave patch mean and the optical patch mean; Determining the timber bamboo forest felling area based on the timber forest felling area identification index; The timber forest felling area identification index is calculated based on the following formula: ; in, represents the timber forest felling area identification index, express Year The average optical patch value in 2018 is express The mean of the annual optical patch mean, express Year The average microwave patch value in 2018 is express The mean of the annual microwave patch means.

[0037] Furthermore, the determining of the timber bamboo forest felling area based on the timber forest felling area identification index includes: If the timber forest felling area identification index is greater than a preset threshold, the monitoring area is determined to be a timber bamboo forest felling area; If the timber forest felling area identification index is less than or equal to the preset threshold, it is determined that the monitoring area is not a timber bamboo forest felling area.

[0038] Specifically, if the timber forest felling area identification index in the monitoring area is greater than the preset threshold, the monitoring area is determined to be a timber bamboo forest felling point; if the timber forest felling area identification index in the monitoring area is not greater than the preset threshold, it is determined to be a non-feeding area.

[0039] During the specific implementation process, the preset threshold can be adaptively set according to actual usage requirements.

[0040] In an optional embodiment, the pixel means of the plurality of microwave patches are calculated based on the following steps: Calculating the plaque mean of each microwave plaque respectively; Determining a pixel mean based on the plaque mean of each of the microwave plaques; Among the plurality of microwave plaques, The mean plaque value of each microwave plaque is calculated based on the following formula: ; in, Indicates The mean value of the microwave plaques, represents the product of the image element values ​​of the microwave, represents the pixel value of VV polarization, Represents the pixel value of VH polarization.

[0041] Here, VV polarization means that both the transmitted radar wave and the received radar wave adopt vertical polarization, and VH polarization means that the transmitted radar wave adopts vertical polarization and the received radar wave adopts horizontal polarization.

[0042] In an optional embodiment, the pixel means of the plurality of optical patches are calculated based on the following steps: Calculating the plaque mean of each optical plaque respectively; Determining a pixel mean based on the patch mean of each of the optical patches; Among the plurality of optical spots, The mean value of each optical patch is calculated based on the following formula: ; in, Indicates The patch mean of the light wave patches, represents the normalized difference vegetation index, represents the reflectivity in the near-infrared band, Indicates the reflectivity of the red light band.

[0043] Based on any of the above embodiments, determining the time point of felling the timber bamboo forest based on the optical patch mean value includes: Calculate the absolute difference of the annual optical patch mean values ​​of any two adjacent years within m years, where the absolute difference of the annual optical patch mean value is the absolute value of the difference of the optical patch mean values ​​of the two years, and m is a positive integer; The year corresponding to the maximum value of the annual optical patch mean absolute difference is taken as the time point for felling the timber bamboo forest, and the maximum value of the annual optical patch mean absolute difference is the maximum value among multiple annual optical patch mean absolute differences.

[0044] In the embodiment of the present invention, the year corresponding to the maximum value of the annual optical patch mean absolute difference is the later year of the two adjacent years corresponding to the maximum value of the annual optical patch mean absolute difference. Taking m=4 as an example, the annual optical patch mean absolute difference between the first year and the second year, the annual optical patch mean absolute difference between the second year and the third year, and the annual optical patch mean absolute difference between the third year and the fourth year are calculated. If the annual optical patch mean absolute difference between the third year and the fourth year is the largest, the fourth year is determined to be the time point for felling the timber bamboo forest.

[0045] Furthermore, the sub-annual time point of felling of timber bamboo forest is determined by the following formula: ; like The largest, then determine the Year is the felling year; in, It is the index for determining the felling time (year) of timber bamboo forest. , , …, Respectively , , …, The average optical patch value for the year.

[0046] Based on any of the above embodiments, determining the time point of felling the timber bamboo forest based on the optical patch mean value includes: Calculate the absolute difference of the monthly optical patch mean values ​​between any two adjacent months within month o, where the absolute difference of the monthly optical patch mean values ​​is the absolute value of the difference between the optical patch mean values ​​of the two months, and o is a positive integer; The month corresponding to the maximum value of the absolute difference of the monthly optical patch means is taken as the time point for felling the timber bamboo forest, and the maximum value of the absolute difference of the monthly optical patch means is the maximum value among the multiple absolute differences of the monthly optical patch means.

[0047] In the embodiment of the present invention, the month corresponding to the maximum value of the monthly optical patch mean absolute difference is the later month of the two adjacent months corresponding to the maximum value of the monthly optical patch mean absolute difference. Taking o=4 as an example, the monthly optical patch mean absolute difference between the first month and the second month, the monthly optical patch mean absolute difference between the second month and the third month, and the monthly optical patch mean absolute difference between the third month and the fourth month are calculated. The monthly optical patch mean absolute difference between the third month and the fourth month is the largest, and the fourth month is determined to be the time point for felling the timber bamboo forest.

[0048] Furthermore, the sub-monthly time point for felling timber bamboo forests is determined by the following formula: ; in, It is the index for determining the time (month) of felling of timber bamboo forest. , , …, Respectively , , …, The average optical patch value for the month.

[0049] It should be understood that the time point for felling bamboo forests can also be determined on a smaller scale based on the satellite data revisit period, such as 5 days, 10 days, etc. Take 10 days as an example: ; in, It is the index for determining the time (month) of felling of timber bamboo forest. , , …, Respectively , , …, The average optical patch value of the month, Is a positive integer.

[0050] Based on any of the above embodiments, the vegetation coverage is calculated based on the following formula: ; in, Indicates the vegetation coverage of timber bamboo forest, Represents the Normalized Difference Vegetation Index NDVI for a specific pixel, represents the minimum value of NDVI in the monitoring area, Indicates the maximum value of NDVI in the monitoring area.

[0051] Specifically, we selected optical images of the dry season (season with less precipitation, such as the dry season from October to March of the following year in southern China) and combined them with ground survey data to estimate the aboveground biomass of bamboo forests, and established a function of aboveground biomass and growth duration of timber bamboo forests: ; ; in, Indicates the growth time of the timber bamboo forest The estimated aboveground biomass at The growth time of timber bamboo forest estimated by combining ground data When the dependent variable (aboveground biomass ) and the independent variable (coverage) function, For the coverage of bamboo forest, is the soil normalized difference vegetation index value.

[0052] The following is a description of the remote sensing estimation device for the felling time, growth duration and aboveground biomass of a timber bamboo forest provided by an embodiment of the present invention. The remote sensing estimation device for the felling time, growth duration and aboveground biomass of a timber bamboo forest described below and the remote sensing estimation method for the felling time, growth duration and aboveground biomass of a timber bamboo forest described above can be referenced to each other.

[0053] Figure 2 Schematic diagram of the structure of the remote sensing estimation device for the felling time point, growth duration and aboveground biomass of a timber bamboo forest provided by the present invention. Figure 2 As shown, the remote sensing estimation device for the felling time point, growth duration and aboveground biomass of the timber bamboo forest may include but is not limited to; The patch segmentation module 210 is used to perform patch segmentation on the microwave remote sensing data and the optical remote sensing data of the monitoring area, respectively, to obtain a plurality of microwave patches and a plurality of optical patches; A mean value calculation module 220 is used to calculate the pixel means of the plurality of microwave patches and the plurality of optical patches respectively, and obtain the microwave patch mean value and the optical patch mean value accordingly; A felling area determination module 230, configured to determine a felling area of ​​the timber bamboo forest based on the microwave patch mean and the optical patch mean; A time point determination module 240, for determining a time point for felling the timber bamboo forest based on the optical patch mean value; A growth duration determination module 250 is used to determine the growth duration of the timber bamboo forest; The biomass determination module 260 is used to determine the aboveground biomass based on the vegetation coverage of the timber bamboo forest and the growth time of the timber bamboo forest.

[0054] The remote sensing estimation device for the felling time point, growth duration and aboveground biomass of a timber bamboo forest provided by the embodiment of the present invention performs patch segmentation on the microwave remote sensing data and the optical remote sensing data of the monitoring area, respectively, to obtain a plurality of microwave patches and a plurality of optical patches; calculates the pixel means of the plurality of microwave patches and the plurality of optical patches, respectively, to obtain a microwave patch mean and an optical patch mean; determines the felling area of ​​the timber bamboo forest based on the microwave patch mean and the optical patch mean. Optical remote sensing can invert the vegetation structure of the bamboo forest, such as leaf area index, vegetation coverage, etc. Microwave remote sensing can identify the geometric information of the bamboo forest, and comprehensively reflect the surface undulation, soil moisture, geometric information, etc. through the backscattering coefficient. This scheme combines the differences between microwave remote sensing data and optical remote sensing data to improve the remote sensing estimation accuracy of the felling time point, growth duration and aboveground biomass of the timber bamboo forest.

[0055] It should be noted that the remote sensing estimation device for the felling time, growth duration and aboveground biomass of a timber bamboo forest provided in an embodiment of the present invention can, during specific operation, execute the remote sensing estimation method for the felling time, growth duration and aboveground biomass of a timber bamboo forest described in any of the above embodiments, which will not be elaborated in this embodiment.

[0056] Figure 3 An example of a physical structure diagram of an electronic device is shown in FIG. Figure 3 As shown, the electronic device may include: a processor 310 , a communications interface 320 , a memory 330 and a communication bus 340 , wherein the processor 310 , the communications interface 320 , and the memory 330 communicate with each other via the communication bus 340 . The processor 310 can call the logic instructions in the memory 330 to execute a remote sensing estimation method for the time point of felling, growth duration and aboveground biomass of a timber bamboo forest, the method comprising: performing patch segmentation on the microwave remote sensing data and the optical remote sensing data of the monitoring area respectively, and correspondingly obtaining a plurality of microwave patches and a plurality of optical patches; calculating the pixel means of the plurality of microwave patches and the plurality of optical patches respectively, and correspondingly obtaining a microwave patch mean and an optical patch mean; determining the timber bamboo forest felling area based on the microwave patch mean and the optical patch mean; determining the time point of felling the timber bamboo forest based on the optical patch mean; determining the growth duration of the timber bamboo forest; and determining the aboveground biomass based on the vegetation coverage of the timber bamboo forest and the growth duration of the timber bamboo forest.

[0057] In addition, the logic instructions in the above-mentioned memory 330 can be implemented in the form of a software functional unit and can be stored in a computer-readable storage medium when it is sold or used as an independent product. Based on such an understanding, the technical solution of the present invention, in essence, or the part that contributes to the prior art or the part of the technical solution, can be embodied in the form of a software product, and the computer software product is stored in a storage medium, including a number of instructions for a computer device (which can be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the method described in each embodiment of the present invention. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), disk or optical disk, etc. Various media that can store program codes.

[0058] On the other hand, the present invention also provides a computer program product, which includes a computer program. The computer program can be stored on a non-transitory computer-readable storage medium. When the computer program is executed by a processor, the computer can execute the remote sensing estimation method for the felling time point, growth time and aboveground biomass of the timber bamboo forest provided by the above-mentioned methods. The method includes: performing patch segmentation on the microwave remote sensing data and the optical remote sensing data of the monitoring area, respectively, to obtain a plurality of microwave patches and a plurality of optical patches; respectively calculating the pixel means of the plurality of microwave patches and the plurality of optical patches, to obtain a microwave patch mean and an optical patch mean; determining the timber bamboo forest felling area based on the microwave patch mean and the optical patch mean; determining the felling time point of the timber bamboo forest based on the optical patch mean; determining the growth time of the timber bamboo forest; and determining the aboveground biomass based on the vegetation coverage of the timber bamboo forest and the growth time of the timber bamboo forest.

[0059] On the other hand, the present invention also provides a non-transitory computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, is implemented to execute the remote sensing estimation method for the felling time point, growth duration and aboveground biomass of a timber bamboo forest provided by the above-mentioned methods, the method comprising: performing patch segmentation on the microwave remote sensing data and the optical remote sensing data of the monitoring area, respectively, to obtain a plurality of microwave patches and a plurality of optical patches correspondingly; calculating the pixel means of the plurality of microwave patches and the plurality of optical patches, respectively, to obtain a microwave patch mean and an optical patch mean correspondingly; determining the felling area of ​​the timber bamboo forest based on the microwave patch mean and the optical patch mean; determining the felling time point of the timber bamboo forest based on the optical patch mean; determining the growth duration of the timber bamboo forest; and determining the aboveground biomass based on the vegetation coverage of the timber bamboo forest and the growth duration of the timber bamboo forest.

[0060] The device embodiments described above are merely illustrative, wherein the units described as separate components may or may not be physically separated, and the components displayed as units may or may not be physical units, that is, they may be located in one place, or they may be distributed on multiple network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the scheme of this embodiment. Ordinary technicians in this field can understand and implement it without paying creative labor.

[0061] Through the description of the above implementation methods, those skilled in the art can clearly understand that each implementation method can be implemented by means of software plus a necessary general hardware platform, and of course, can also be implemented by hardware. Based on this understanding, the above technical solution is essentially or the part that contributes to the prior art can be embodied in the form of a software product, and the computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, a disk, an optical disk, etc., including a number of instructions for a computer device (which can be a personal computer, a server, or a network device, etc.) to execute the methods described in each embodiment or some parts of the embodiments.

[0062] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A remote sensing estimation method for the felling time, growth duration and aboveground biomass of a timber bamboo forest, characterized in that: include: Performing patch segmentation on the microwave remote sensing data and the optical remote sensing data of the monitoring area respectively, and obtaining a plurality of microwave patches and a plurality of optical patches accordingly; Calculating pixel means of the plurality of microwave patches and the plurality of optical patches respectively, and obtaining microwave patch means and optical patch means correspondingly; Determining the felling area of ​​the timber bamboo forest based on the microwave patch mean and the optical patch mean; Determining the time point of felling the timber bamboo forest based on the optical patch mean value; Determine the growth period of timber bamboo forest; The aboveground biomass is determined based on the vegetation coverage of the timber bamboo forest and the growth duration of the timber bamboo forest.

2. The remote sensing estimation method for felling time, growth duration and aboveground biomass of timber bamboo forest according to claim 1 is characterized in that: The determining of the timber bamboo forest felling area based on the microwave patch mean and the optical patch mean comprises: Determining a timber forest felling area identification index based on the microwave patch mean and the optical patch mean; Determining the timber bamboo forest felling area based on the timber forest felling area identification index; The timber forest felling area identification index is calculated based on the following formula: ； in, represents the timber forest felling area identification index, express Year The average optical patch value in 2018 is express The mean of the annual optical patch mean, express Year The average microwave patch value in 2018 is express The mean of the annual microwave patch means.

3. The remote sensing estimation method for felling time, growth duration and aboveground biomass of timber bamboo forest according to claim 2, characterized in that: The step of determining the timber bamboo forest felling area based on the timber forest felling area identification index comprises: If the timber forest felling area identification index is greater than a preset threshold, the monitoring area is determined to be a timber bamboo forest felling area; If the timber forest felling area identification index is less than or equal to the preset threshold, it is determined that the monitoring area is not a timber bamboo forest felling area.

4. The remote sensing estimation method for felling time, growth duration and aboveground biomass of timber bamboo forest according to claim 1, characterized in that: The pixel mean of the plurality of microwave patches is calculated based on the following steps: Calculating the plaque mean of each microwave plaque respectively; Determining a pixel mean based on the plaque mean of each of the microwave plaques; Among the plurality of microwave plaques, The mean plaque value of each microwave plaque is calculated based on the following formula: ； in, Indicates The mean value of the microwave plaques, represents the product of the image element values ​​of the microwave, represents the pixel value of VV polarization, Pixel value representing VH polarization.

5. The remote sensing estimation method for felling time, growth duration and aboveground biomass of a timber bamboo forest according to any one of claims 1 to 4, characterized in that: The step of determining the time point of felling the timber bamboo forest based on the optical patch mean value comprises: Calculate the absolute difference of the annual optical patch mean values ​​of any two adjacent years within m years, where the absolute difference of the annual optical patch mean value is the absolute value of the difference of the optical patch mean values ​​of the two years, and m is a positive integer; The year corresponding to the maximum value of the annual optical patch mean absolute difference is taken as the time point for felling the timber bamboo forest, and the maximum value of the annual optical patch mean absolute difference is the maximum value among multiple annual optical patch mean absolute differences.

6. The remote sensing estimation method for felling time, growth duration and aboveground biomass of a timber bamboo forest according to any one of claims 1 to 4, characterized in that: The step of determining the time point of felling the timber bamboo forest based on the optical patch mean value comprises: Calculate the absolute difference of the monthly optical patch mean values ​​between any two adjacent months within month o, where the absolute difference of the monthly optical patch mean values ​​is the absolute value of the difference between the optical patch mean values ​​of the two months, and o is a positive integer; The month corresponding to the maximum value of the absolute difference of the monthly optical patch means is taken as the time point for felling the timber bamboo forest, and the maximum value of the absolute difference of the monthly optical patch means is the maximum value among the multiple absolute differences of the monthly optical patch means.

7. The remote sensing estimation method for felling time, growth duration and aboveground biomass of a timber bamboo forest according to any one of claims 1 to 4, characterized in that: The vegetation coverage is calculated based on the following formula: ； in, Indicates the vegetation coverage of timber bamboo forest, Represents the Normalized Difference Vegetation Index NDVI for a specific pixel, represents the minimum value of NDVI in the monitoring area, Indicates the maximum value of NDVI in the monitoring area.

8. A remote sensing estimation device for the felling time, growth duration and aboveground biomass of a timber bamboo forest, characterized in that: include: The patch segmentation module is used to perform patch segmentation on the microwave remote sensing data and the optical remote sensing data of the monitoring area, and obtain a plurality of microwave patches and a plurality of optical patches accordingly; A mean value calculation module, used to calculate the pixel means of the plurality of microwave patches and the plurality of optical patches respectively, and obtain the microwave patch mean value and the optical patch mean value accordingly; A felling area determination module, used for determining a felling area of ​​a timber bamboo forest based on the microwave patch mean and the optical patch mean; A time point determination module, used to determine the time point of felling the timber bamboo forest based on the optical patch mean value; A growth duration determination module is used to determine the growth duration of the timber bamboo forest; The biomass determination module is used to determine the aboveground biomass based on the vegetation coverage of the timber bamboo forest and the growth time of the timber bamboo forest.

9. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein: When the processor executes the computer program, the remote sensing estimation method for the felling time point, growth duration and aboveground biomass of the timber bamboo forest as claimed in any one of claims 1 to 7 is implemented.

10. A non-transitory computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the remote sensing estimation method for the felling time point, growth duration and aboveground biomass of a timber bamboo forest as claimed in any one of claims 1 to 7 is implemented.

11. A computer program product, comprising a computer program, characterized in that When the computer program is executed by a processor, the remote sensing estimation method for the felling time point, growth duration and aboveground biomass of a timber bamboo forest as claimed in any one of claims 1 to 7 is implemented.

Images