Three, summary of the invention
The drawbacks of in order to overcome traditional biological amount field observation mode and existing contactless remote sensing technology biometric measurementExisting limitation is calculated, the object of the present invention is to provide a kind of Biomass Models based on high-spectrum remote-sensing and photogrammetric technologyGrind construction method.
BROAD SUMMARY:
A kind of Biomass Models based on high-spectrum remote-sensing and photogrammetric technology grind construction method, it is characterized in that: benefit respectivelyVegetation characteristics parameters are extracted with high spectrum resolution remote sensing technique, Unmanned Aerial Vehicle Photogrammetric Technique and terrestrial photogrammetric survey technology,Combined ground sample ground measured data, establishes the inverse model between vegetation characteristics parameters and biomass respectively, and to grinding the mould builtType carries out precision analysis, determines optimal Biomass retrieval model.
Wherein, the step of extracting vegetation characteristics parameters using high spectrum resolution remote sensing technique are as follows:
(1) spectroscopic data pre-processes, and selects the average value of the measured value of a certain range before and after a measuring point on the curve of spectrum, makeesFor the value of the point, the sequence { R of N number of measuring point is providedi, i=1,2,3 ... N }, wherein the value of i point includes each K point in front and backAverage value, utilize mathematical model (1)Wherein, the new value of i point isR′i, substituted with the average value of the 2k+1 point including this point, and then determine the value of the point on the curve of spectrum;
(2) high-spectral data characteristic parameter extraction utilizes mathematical model (2) FDRλ(r)=(Rλ(j+1)-Rλ(j))/Δ λ is to bloomModal data is handled, wherein FDRλ(r)For the reflectivity first derivative at the wavelength i of the midpoint wave band j and j+1, Rλ(j+1)Wave bandValue is the reflectivity of j+1, Rλ(j)Band value is that the reflectivity Δ λ of j is the interval of wavelength j+1 to j, depending on band wavelength, and it is relatedCoefficient is expressed as mathematical model (3)And then obtain the differential curves of multiple orders, can effective green vegetation spectral signature;
(3) mathematical model (4) are utilizedEO-1 hyperion vegetation index model is established,In, R900、R1050The wavelength of reflectivity maximum respectively within the scope of near-infrared and short infrared wave band, R955、R1220RespectivelyThe wavelength at reflectivity minimum within the scope of near-infrared and short infrared wave band, it is quick between enhanced spectrum feature and vegetation biomassPerception.
Wherein, the detailed process of vegetation characteristics parameters is extracted using Unmanned Aerial Vehicle Photogrammetric Technique are as follows: need surrounding targetThe spiral mode of rising of trees carries out oblique photograph measurement, by oblique photograph measuring technique collect the trees different location,The image data of different angle rebuilds the three-dimensional point cloud model of ground vegetation, wherein between the adjacent image of sustained height acquisitionReach 90% or more Duplication, reaches 60% or more Duplication between the image of different height acquisition, phase when oblique photographMachine best angle is 45 degree;According to feature in three dimensional point cloud successively selected characteristic point, using Related Mathematical Models respectively intoRow tree height and the measuring and calculating of trunk any place diameter, crown mapping area and the measuring and calculating of hat width, the measuring and calculating of Crown surface area, Tree Crown VolumeMeasuring and calculating.
Wherein, the step of extracting vegetation characteristics parameters using terrestrial photogrammetric survey technology are as follows: with common optics or numberCamera is tool, establishes digital close view photogrammetric system by computer, and system photographs to permanent sample plot by general camera,Within foundation, elements of exterior orientation (x0, y0, f, xs, ys, zs, X, K, U) and it is the photogrammetric collinearity equation of unknown number, become by DLTIt changes, undetermined parameter is acquired, and then image space and object space are connected, by computer system, so that it may determine image spaceThe corresponding object coordinates of any point coordinate (x, y) (X, Y, Z) calculate tree height, the diameter of a cross-section of a tree trunk 1.3 meters above the ground, tree crown to restore object space solidVolume surveys tree index.
Wherein, Biomass Models, which are ground, builds detailed process are as follows:
(1) common Biomass Models W=aD is utilizedb, W=a (D2H)b, W=aDbHc、Wherein, W indicates that biomass, D indicate the diameter of a cross-section of a tree trunk 1.3 meters above the ground, and H indicates that tree is high, CwIndicate hat width, a and b are estimation parameter, by EO-1 hyperion, nobodyThe individual tree information that machine, terrestrial photogrammetric survey technology are extracted establishes Biomass retrieval model in conjunction with the measured data of ground;
(2) common model-evaluation index coefficient R is utilized2, the standard error SEE of estimated value, total relative errorTRE, average systenatic error MSE, average predictor error MPE need to evaluate its model to the Biomass Models of foundation,To select optimal models.
This invention has the advantage that compared with the conventional method
(1) change traditional biological amount and calculate mode on the spot, reduce field process amount, improve measurement efficiency;
(2) in terms of measurement range, big region is promoted to by original small range measuring and calculating and quickly, is efficiently calculated;
(3) evade the limitation of the prior art, mutually use for reference and utilize, realize accurate, the efficient measuring and calculating of biomass.
Five, specific embodiment:
The present invention is further illustrated With reference to embodiment, specifically:
Vegetation characteristics parameters are extracted using high spectrum resolution remote sensing technique, grind the step of building Biomass Models are as follows:
(1) spectroscopic data pre-processes, and selects the average value of the measured value of a certain range before and after a measuring point on the curve of spectrum, makeesFor the value of the point, the sequence { R of N number of measuring point is providedi, i=1,2,3 ... N }, wherein the value of i point includes each K point in front and backAverage value, utilize mathematical model (1)Wherein, the new value of i point isR′i, substituted with the average value of the 2k+1 point including this point, and then determine the value of the point on the curve of spectrum;
(2) high-spectral data characteristic parameter extraction utilizes mathematical model (2) FDRλ(r)=(Rλ(j+1)-Rλ(j))/Δ λ is to bloomModal data is handled, wherein FDRλ(r)For the reflectivity first derivative at the wavelength i of the midpoint wave band j and j+1, Rλ(j+1)Wave bandValue is the reflectivity of j+1, Rλ(j)Band value is that the reflectivity Δ λ of j is the interval of wavelength j+1 to j, depending on band wavelength, and it is relatedCoefficient is expressed as mathematical model (3)And then obtain the differential curves of multiple orders, can effective green vegetation spectral signature;
(3) mathematical model (4) are utilizedEO-1 hyperion vegetation index model is established,In, R900、R1050The wavelength of reflectivity maximum respectively within the scope of near-infrared and short infrared wave band, R955、R1220RespectivelyThe wavelength at reflectivity minimum within the scope of near-infrared and short infrared wave band, it is quick between enhanced spectrum feature and vegetation biomassPerception;
(4) common Biomass Models W=aD is utilizedb, W=a (D2H)b, W=aDbHc、Wherein, W indicates that biomass, D indicate the diameter of a cross-section of a tree trunk 1.3 meters above the ground, and H indicates that tree is high, CwIndicate hat width, a and b are estimation parameter, by EO-1 hyperion, nobodyThe individual tree information that machine, terrestrial photogrammetric survey technology are extracted establishes Biomass retrieval model in conjunction with the measured data of ground;
(5) common model-evaluation index coefficient R is utilized2, the standard error SEE of estimated value, total relative errorTRE, average systenatic error MSE, average predictor error MPE need to evaluate its model to the Biomass Models of foundation,To select optimal models.
Vegetation characteristics parameters are extracted using Unmanned Aerial Vehicle Photogrammetric Technique, grind the step of building Biomass Models are as follows:
(1) as shown in Fig. 2, surrounding target trees spiral, the mode of rising carries out oblique photograph measurement, passes through oblique photographMeasuring technique collects the image data of the trees different location, different angle, rebuilds the three-dimensional point cloud model of ground vegetation,In, 90% or more Duplication is reached between the adjacent image of sustained height acquisition, is reached between the image of different height acquisition60% or more Duplication, camera best angle when oblique photograph are 45 degree;
(2) tree height and the detailed process of trunk any place diameter measuring and calculating are as follows:
1. i indicates point number, value 1,2 ..., n, if certain point i on trunk on the three dimensional point cloud of target treesCoordinate be (Xi, Yi, Zi), the coordinate of another point i+1 is (Xi+1, Yi+1, Zi+1), it utilizes mathematical model (5)Acquire each Duan Shugao, wherein H indicates each section of tree high level, and i and i+1 is respectively about two points on the same vertical line of trunk at this time;
2. utilizing mathematical model (6)Acquire any place tree-walk value, wherein D is indicatedTree-walk value, i and i+1 is respectively two left-right dots in trunk same level at this time.
(3) detailed process of the measuring and calculating of crown mapping area and hat width are as follows:
1. the three dimensional point cloud of target trees is carried out hierarchy slicing, as shown in Fig. 3 a, 3b, if i-th on a certain sliceThe coordinate of a point is (Xi, Yi), the coordinate of i+1 point is (Xi+1, Yi+1), then utilize mathematical model (7)Calculate slice area, wherein i indicate point number, value 1,2 ...,n;
2., when slice area reaches maximum, maximum at this time is sliced area S as shown in Fig. 4 a, 4bmaxMean that tree crownProjected area;
3. utilizing mathematical model (8)Acquire the hat width of tree crown, wherein K indicates hatWidth, (Xa, Ya)、(Xb, Yb) respectively indicate lie farthest away on crown mapping section two o'clock a and b coordinate.
(4) detailed process of the measuring and calculating of Crown surface area are as follows:
1. utilizing mathematical model (9)It is respectively sliced after calculating tree crown layeringPerimeter, wherein (Xi, Yi) it is to be sliced upper i-th point of coordinate, (Xi+1, Yi+1) be i+1 point coordinate, i expression point number, takeValue is 1,2 ..., n;
It is made of 2. regarding as the tree crown after layering several intermediate stage bodies and top and the cone of bottom two, stage bodySide expansion be approximately it is trapezoidal, the expansion of the side of centrum is approximately triangle, is utilized respectively mathematical model (10)(11)Calculate the lateral area of approximate trapezoid and triangle, wherein j expression is cutPiece number, value are 1 ..., m, and when indicating stage body slice, value is 2 ..., m-1, and when indicating cone, value is 1 or m, Lj、Lj+1Respectively indicate the slice perimeter of jth layer and+1 layer of jth, hjIndicate that the side of folded rule body between adjacent two layers are sliced is high;
3. the lateral area of each rule body acquired is added up, utilize mathematical model (12)Acquire tree crown surfaceProduct, wherein Q is Crown surface area.
(5) detailed process of the measuring and calculating of Tree Crown Volume are as follows:
1. tree crown is layered, the tree crown after layering is regarded as by several intermediate stage bodies and top and the cone of bottom twoIt constitutes, utilizes mathematical model (13)(14)Calculate approximate platformThe volume of body and cone, wherein j indicates that slice number, value are 1 ..., m, and value is 2 ..., m- when indicating stage body slice1, when indicating cone, value is 1 or m, Sj、Sj+1Respectively indicate the slice area of jth layer and+1 layer of jth, HjIndicate adjacent two layersSegmentation height between slice;
2. the volume of each rule body acquired is added up, utilize mathematical model (15)Crown surface area is acquired,Wherein, V is Crown surface area;
(6) common Biomass Models W=aD is utilizedb, W=a (D2H)b, W=aDbHc、Wherein, W indicates that biomass, D indicate the diameter of a cross-section of a tree trunk 1.3 meters above the ground, and H indicates that tree is high, CwIndicate hat width, a and b are estimation parameter, by EO-1 hyperion, nobodyThe individual tree information that machine, terrestrial photogrammetric survey technology are extracted establishes Biomass retrieval model in conjunction with the measured data of ground;
(7) common model-evaluation index coefficient R is utilized2, the standard error SEE of estimated value, total relative errorTRE, average systenatic error MSE, average predictor error MPE need to evaluate its model to the Biomass Models of foundation,To select optimal models.
Vegetation characteristics parameters are extracted using terrestrial photogrammetric survey technology, grind the step of building Biomass Models are as follows:
(1) using common optics or digital camera as tool, digital close view photogrammetric system is established by computer, isSystem photographs to permanent sample plot by general camera, within foundation, elements of exterior orientation (x0, y0, f, xs, ys, zs, X, K, U) and it is unknownThe photogrammetric collinearity equation of number, is converted by DLT, acquires undetermined parameter, and then image space and object space are connected,By computer system, so that it may the corresponding object coordinates of image space any point coordinate (x, y) (X, Y, Z) are determined, to restore object spaceSpace multistory;
(2) specific steps of the high measuring and calculating of tree are as follows:
1. the flat position apart from the high left and right of 1 times of object to be measured tree tree is selected to guarantee goal tree and instrument as observation pointBetween intervisibility it is unobstructed;
2. frame level device, and measure the high h of instrument at this time;
3. as shown in figure 5, alignment target tree, carries out horizontal honest shooting.
(3) specific steps of diameter of a cross-section of a tree trunk 1.3 meters above the ground measuring and calculating are as follows:
1. selecting suitable observation point, frame level device measures height of instrument at this time;
2. as shown in fig. 6, screen center alignment target tree is carried out horizontal honest shooting, dragging photo horizontal line V and groundUpper thread VgIt is overlapped, the anti-diameter of a cross-section of a tree trunk 1.3 meters above the ground horizontal line released at 1.3m;
3. it is left along coordinate (U to choose the 1.3m diameter of a cross-section of a tree trunk 1.3 meters above the ground manuallyl, V1.3) and it is right along coordinate (Ur, V1.3);
4. utilizing mathematical model (16)Calculate the diameter of a cross-section of a tree trunk 1.3 meters above the ground d of goal tree.
(4) specific steps of Tree Crown Volume measuring and calculating are as follows:
1. setting height-age curve domain as D, tree crown is split with bis- straight line of X, Y is parallel to, tree crown is divided into n notEquitant daughter, as shown in Figure 7;
2. calculating the Tree Crown Volume that each daughter is controlled, and add up as n → ∞, mathematical model (17) V=∫ ∫ is utilizedHidxdy, acquire entire Tree Crown Volume, wherein dxdyFor the floor space of daughter,
(5) common Biomass Models W=aD is utilizedb, W=a (D2H)b, W=aDbHc、Wherein, W indicates that biomass, D indicate the diameter of a cross-section of a tree trunk 1.3 meters above the ground, and H indicates that tree is high, CwIndicate hat width, a and b are estimation parameter, by EO-1 hyperion, nobodyThe individual tree information that machine, terrestrial photogrammetric survey technology are extracted establishes Biomass retrieval model in conjunction with the measured data of ground;
(6) common model-evaluation index coefficient R is utilized2, the standard error SEE of estimated value, total relative errorTRE, average systenatic error MSE, average predictor error MPE need to evaluate its model to the Biomass Models of foundation,To select optimal models.