| CPC G06Q 30/0203 (2013.01) [G01N 33/0098 (2013.01)] | 6 Claims |
|
1. A continuous monitoring method for forest stock, comprising the following steps:
(1) sampling design of sample plots: completing layout of monitoring sample plots, and determining a sample plot population, a sampling method, and spatial locations of sample plots;
(2) layout of intelligent sample plots: completing layout of first measurement and monitoring devices for the sample plots; the layout of the monitoring devices comprising:
installing tree diameter measurement sensors to measure diameters and perimeters of sample trees;
connecting the tree diameter measurement sensors to data collection terminals through wireless an ad hoc network technology to complete data collection and input;
networking the tree diameter measurement sensors and the data collection terminals with a mobile communication gateway or a Beidou short message gateway through the wireless ad hoc network technology to complete data summarization and transmission;
sending, by the mobile communication gateway, data back to a communication server through a mobile communication base station of a communication operator, or sending, by the Beidou short message gateway, the data to a Beidou director through a Beidou satellite; and
transmitting, by the communication server or the Beidou director, the data to a continuous monitoring application system for forest stock through an optical network;
(3) automatic collection of sample plot data:
(4) dynamic update of stock: detecting plot type change subclasses through remote sensing, and updating graphic and attribute information of forest resource change maps simultaneously; building a dynamic forest stand update model through intelligent sample plot data for plot type unchanged subclasses, and then updating attribute information of forest subclasses, wherein specific steps are as follows:
(4.1) computing the stock of each intelligent sample plot according to the data collected in step (3);
(4.2) computing forest stock and sampling precision of current monitored regions based on the intelligent sample plot, wherein the forest stock is computed by the following formula:
 where Vall_Plot is the forest stock of the current monitored regions based on the intelligent sample plot, vij is the stock of the ith sample plot of the jth population, sij is the area of the ith sample plot of the jth population, Sj is the total area of the jth population, n is the number of sample plots of the jth population, and m is the total number of the current monitored regions;
wherein the sampling precision is computed by the following formula:
 where Pvj is sampling precision of the current sample plot population, ta is a reliability index, Svj is an arithmetic square root of a sample variance of the jth population, and Vj is a sample mean of the jth population;
(4.3) determining whether the sampling precision meets the sampling design, and if so, performing step (4.4); otherwise, performing step (1) to adjust the sampling design and complement intelligent sample plots;
(4.4) determining whether to combine with first class survey, and if so, computing forest stock of the current monitored regions based on the intelligent sample plot as current period forest stock, and performing step (4.17); otherwise, performing step (4.5);
(4.5) determining whether remote sensing images are obtained in a monitoring period, and if so, performing step (4.6); otherwise, only updating the dynamic forest stand model and performing step (4.11);
(4.6) carrying out remote sensing change detection and update, with remote sensing change detection as the main approach, supplemented by on-site survey and file update, and zoning a spatial scope of plot type change subclasses to form a remote sensing interpretation map database;
(4.7) filling in on-site survey factors for remote sensing interpretation maps based on on-site survey and file update to form an on-site survey database;
(4.8) performing spatial update analysis on the on-site survey database and a base period forest resource subclass database, and performing spatial and attribute updates of the on-site survey factors on the base period forest resource subclass database to generate a current period forest resource subclass database;
(4.9) performing spatial joint analysis on the current period forest resource subclass database and the base period forest resource subclass database, and only retaining previous and subsequent plot type change subclasses for the joint results as a forest resource change database;
(4.10) summarizing differences between the current period subclass stock and the base period subclass stock of the forest resource change database to obtain a subclass stock variation of the plot type change subclasses, wherein a computation formula is as follows:
 Where ΔVArea_change is the subclass stock variation of the plot type change subclasses, Vcur_areai is the subclass stock of the ith subclass in the forest resource change database, Vbase_areai is base period subclass stock of the ith subclass in the forest resource change database, and 0 is the total number of subclasses in the forest resource change database;
(4.11) determining whether to combine with second class survey, and if so, performing step (4.12); otherwise, performing step (4.14);
(4.12) computing the total stock of sample plots in the plot type unchanged subclasses, and solving the difference between the ratio of the total stock to the total stock of base period sample plots and 1 as a dynamic forest stand model update parameter, wherein a computation formula is as follows:
 where k is the dynamic forest stand model update parameter, Vcur_Plot is the total stock of sample plots in the plot type unchanged subclasses, and Vbase_Plot is the total stock of the base period sample plots;
(4.13) multiplying the stock of the plot type unchanged subclasses by the dynamic forest stand model update parameter, and obtaining a forest stock variation of the forest stand structure change subclasses after summarization, wherein a computation formula is as follows:
 where ΔVStruct_vary is the forest stock variation of the forest stand structure change subclasses, ΔVbase_Struct is the stock of the plot type unchanged subclasses, Vcur_Plot is the total stock of sample plots in the plot type unchanged subclasses, and Vbase_Plot is the total stock of the base period sample plots;
(4.14) computing the total stock of sample plots in the plot type unchanged subclasses in each stratum or quota, and solving the difference between the ratio of the total stock to the total stock of base period sample plots and 1 as a dynamic forest stand model update parameter in each stratum or quota, wherein a computation formula is as follows:
 where kj is the dynamic forest stand model update parameter in the jth stratum or quota, Vcur_Plotj is the total stock of sample plots in the plot type unchanged subclasses in the jth stratum or quota, and Vbase_Plotj is the total stock of base period sample plots in the jth stratum or quota;
(4.15) multiplying base subclass stock of the plot type unchanged subclasses in each stratum or quota by the dynamic forest stand model update parameter, and obtaining a forest stock variation of the forest stand structure change subclasses in each stratum or quota after summarization, wherein a computation formula is as follows:
 where ΔVStruct_vary is the forest stock variation of the forest stand structure change subclasses, Vcur_Plotj is the total stock of sample plots in the plot type unchanged subclasses in the jth stratum or quota, Vbase_Plotj is the total stock of base period sample plots in the jth stratum or quota, Vbase_Structij is base period subclass stock of the ith plot type unchanged subclass in the jth stratum or quota, Vcur_Plotj is the stock of sample plots in the plot type unchanged subclasses in the jth stratum or quota, VbasePlot is the stock of base period sample plots in the jth stratum or quota, s is a quantity of plot type unchanged subclasses in the jth stratum or quota, and t is a quantity of strata or quotas;
(4.16) computing a forest stock variation of forest resource subclasses by the following formula:
ΔV=ΔVArea_change+ΔVStruct_vary
where ΔV is the forest stock variation of the forest resource subclasses, ΔVArea_change is the subclass stock variation of the plot type change subclasses, and ΔVStruct_vary is the forest stock variation of the forest stand structure change subclasses;
(4.17) ending the process;
(6) monitoring output: outputting current period stock monitoring data;
(7) determining whether a monitoring period arrives, and if so, ending the process; otherwise, returning to step (3) for continuous monitoring.
|