CN113109832A - Bamboo degree judgment method based on laser echo intensity - Google Patents
Bamboo degree judgment method based on laser echo intensity Download PDFInfo
- Publication number
- CN113109832A CN113109832A CN202110291159.3A CN202110291159A CN113109832A CN 113109832 A CN113109832 A CN 113109832A CN 202110291159 A CN202110291159 A CN 202110291159A CN 113109832 A CN113109832 A CN 113109832A
- Authority
- CN
- China
- Prior art keywords
- bamboo
- echo intensity
- laser echo
- laser
- value
- 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
- 235000017166 Bambusa arundinacea Nutrition 0.000 title claims abstract description 78
- 235000017491 Bambusa tulda Nutrition 0.000 title claims abstract description 78
- 241001330002 Bambuseae Species 0.000 title claims abstract description 78
- 235000015334 Phyllostachys viridis Nutrition 0.000 title claims abstract description 78
- 239000011425 bamboo Substances 0.000 title claims abstract description 78
- 238000000034 method Methods 0.000 title claims abstract description 39
- 244000302661 Phyllostachys pubescens Species 0.000 claims abstract description 46
- 235000003570 Phyllostachys pubescens Nutrition 0.000 claims abstract description 46
- 238000012937 correction Methods 0.000 claims abstract description 24
- 238000005259 measurement Methods 0.000 claims abstract description 13
- 238000002310 reflectometry Methods 0.000 claims description 7
- 241000196324 Embryophyta Species 0.000 claims description 6
- 230000008569 process Effects 0.000 description 5
- 238000001228 spectrum Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000035772 mutation Effects 0.000 description 2
- 241001330024 Bambusoideae Species 0.000 description 1
- 241000544043 Blyxa aubertii Species 0.000 description 1
- 241000218631 Coniferophyta Species 0.000 description 1
- 241000745988 Phyllostachys Species 0.000 description 1
- 241000209504 Poaceae Species 0.000 description 1
- 206010048245 Yellow skin Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 210000000481 breast Anatomy 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 241001233957 eudicotyledons Species 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 230000002503 metabolic effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003071 parasitic effect Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 241000894007 species Species 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 239000002023 wood Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/88—Lidar systems specially adapted for specific applications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/48—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
- G01S7/4802—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Electromagnetism (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
The invention discloses a moso bamboo degree judgment method based on laser echo intensity, and relates to the technical field of laser moso bamboo degree judgment. The method comprises the steps of establishing a sample plot of a rectangular area, and erecting a three-dimensional laser scanner by taking the periphery and the center of the sample plot as a measuring station; extracting point cloud data from the middle of each section of bamboo joint in a window of 2cm by 2 cm; correcting the distance measurement value and the incident angle of the laser echo strength value; and obtaining a fitting model of bamboo joint numbers of bamboo stalks with different bamboo degrees and the laser echo intensity value. The method comprises the steps of obtaining point cloud data of the middle parts of all bamboo sections of the moso bamboo through a three-dimensional laser scanner, and correcting distance measurement values and incidence angles of laser echo intensity values; calculating the average value of the laser echo intensity values after correction as the laser echo intensity value of the bamboo joint section, and obtaining a fitting model of bamboo joint numbers of bamboo stalks with different bamboo degrees and the laser echo intensity value; the bamboo degree of the moso bamboo is judged by calculating whether the laser echo intensity value of each section of the moso bamboo meets the fitting model, and the discrimination efficiency and accuracy of the bamboo degree of the moso bamboo are improved.
Description
Technical Field
The invention belongs to the technical field of laser moso bamboo degree judgment, and particularly relates to a moso bamboo degree judgment method based on laser echo intensity.
Background
The moso bamboo as a large-scale scattered bamboo seed of Phyllostachys of Bambusoideae of Gramineae has the characteristics of fast growth, high yield, good quality, wide application, short felling period, early yield and high economic value, and is the bamboo seed which is widely distributed in the world, has the largest cultivation area and has obvious economic benefit at present. The xylem of the plant is different from that of a common arbor tree species in that the xylem of the plant has a traditional annual ring to determine the age unlike gymnosperms or dicotyledons. The age of the bamboo plant of the moso bamboo (bamboo age for short) has a close relationship with the living strength, the metabolic capability, the renewal capability and the mechanical properties of the bamboo, so that the accurate judgment of the age of the moso bamboo has great significance for cultivation management, felling utilization, scientific research and the like of the moso bamboo, and is one of the essential bases.
The age of bamboo is commonly expressed by degree, one degree represents two years, the moso bamboo grown from one to two years is called first degree bamboo, the moso bamboo grown from three to four years is called second degree bamboo, and so on. The currently used bamboo age determination methods mainly comprise three methods:
1) the bamboo number marking method comprises the following steps: the method can accurately judge the age of the moso bamboo by marking the bamboo forming age on the bamboo stalk when the moso bamboo is formed by bamboo shoots, but the workload is large and the condition of label leakage is easy to cause;
2) bamboo stalk color skin method: as the moso bamboo culm grows older, its surface color, wax layer and parasitic lichen characteristics of the moso bamboo change. The first-degree bamboo stalks are cyan and have white fluff adhered to the surfaces, the second-degree bamboo stalks are green yellow, the third-degree bamboo stalks are yellow skins, and the fourth-degree bamboo stalks are over 6 years old and have yellow-red colors. The method is simple and easy to identify, but the age of the moso bamboo is easy to be judged by mistake;
3) identifying branch leaf-changing habit: the twig of the first generation moso bamboo is an unshaped twig, new branches and leaves grow with the age of the mother bamboo, old branches fall off in the process of the growth and development of the new branches, traces are left on the branches, and the formula for calculating the age of the moso bamboo can be obtained by induction and arrangement, wherein the age of the moso bamboo is (2 multiplied by the number of the traces of the branches) + the age of the small branches with leaves. The method can accurately calculate the real age of the moso bamboos, but requires the patience analysis, needs to have rich plant knowledge and practical experience, is easily influenced by external environment mutation, and cannot be used for accurately judging the age of the moso bamboos under the influence of the external environment mutation such as locust disasters, fire drought and waterlogging and the like.
The invention provides a moso bamboo degree judging method based on laser echo intensity, aiming at efficiently judging the moso bamboo degree.
Disclosure of Invention
The invention aims to provide a moso bamboo degree judging method based on laser echo intensity, which is used for solving the problems in the background technology.
In order to solve the technical problems, the invention is realized by the following technical scheme:
the invention relates to a moso bamboo degree judging method based on laser echo intensity, which comprises the following steps:
a00: establishing a sample plot of a rectangular area, and erecting a three-dimensional laser scanner by taking the periphery and the center of the sample plot as a measuring station;
a01: selecting 8 plants of one to seven-year-old moso bamboos, cutting the parts from the lower branches of the moso bamboos to the bamboo stems on the ground, and extracting point cloud data from the middle part of each section of bamboo joint in a window of 2cm to 2 cm; the point cloud data consists of a plurality of laser echo intensity values;
a02: correcting the distance measurement value and the incident angle of the laser echo strength value;
a03: calculating the average value of the corrected laser echo intensity values as the laser echo intensity value of the section of bamboo joint;
a04: obtaining fitting models of bamboo joint numbers of bamboo stalks with different bamboo degrees and laser echo intensity values:
y1=0.48x2-17.937x-1014.5
y2=-0.4425x2+12.983x-1157.8
y3=-0.0547x2-4.2539x-1050.8
y4=0.0016x6-0.0849x5-19.105x3+109.77x2-315.42x-680.66
wherein, y1、y2、y3、y4Respectively fitting laser echo strength values of the bamboo stalks of one to four degrees; and x is a bamboo knot number.
As a preferred technical solution, a00 further includes placing 3 plane targets in the sample chamber and performing station-finding stitching through the plane targets.
As a preferred technical solution, a02 specifically includes the following processes:
b00: neglecting the influence of the reflectivity of the object and the atmospheric attenuation on the laser echo intensity value, and obtaining a simplified laser echo intensity empirical model: i ═ Gf (θ) f (r);
wherein, I is the original laser echo intensity value, G is a constant, f (theta) is an incidence angle correction function, and f (R) is a distance correction function;
b01: fixing the incidence angle theta of the three-dimensional laser scanner, selecting a 3In plane target as a scanning object, changing the distance between the target and the scanner by a fixed step length, keeping the center of an electric detector of the three-dimensional laser scanner and the center of the target at the same height for scanning, and obtaining the relation between different distance measurement and laser echo intensity values under the fixed incidence angle;
b02: fixing the incident distance R of the three-dimensional laser scanner, rotating the target plane by a fixed step length, changing the incident angle of the scanned target, and acquiring the relationship between different incident angles and laser echo intensity values at a fixed distance;
b03: and establishing a piecewise fitting function between the laser echo intensity value and the ranging and incidence angles, calculating model parameters, and correcting the laser echo intensity value.
As a preferable technical scheme, a target plane in B01 is vertical to the ground, and an incidence angle theta is kept to be 0; when R is less than or equal to 10m, the step length is 1 m; when R is more than 10m, the step length is 5 m.
As a preferable technical scheme, the distance between the three-dimensional laser scanner in B02 and the target is 10 m; when theta is less than or equal to 40 degrees, the step length is 10 degrees; theta > 40 deg. and the step size is 5 deg..
As a preferred technical solution, B03 specifically includes the following processes:
c00: selecting point cloud data with the distance of 8-10m and the incident angle of 40-50 degrees to perform polynomial fitting, and acquiring an extreme point of a piecewise fitting function as a critical point; wherein the distance correction critical point is 9.9m, and the angle correction critical point is 45 degrees;
c01: obtaining a laser echo intensity value correction formula containing distance and incident angle:
wherein, IsFor the corrected laser echo intensity value, (X, Y, Z) is the three-dimensional coordinate of the point cloud space, RsTo reference the distance, θsAs a reference angle, αiAnd betaiAre all polynomial coefficients, | n | is the modulus of the normal vector; i is the original laser echo intensity; k is a system constant, theta is an incident angle, and R is an incident distance;
c02: obtaining a point cloud normal vector n ═ (n) according to a K nearest neighbor domain method1,n2,n3) Solving the cosine value of the point incidence angle according to the cosine law, and solving the inverse cosine function of the cosine value to obtain the corresponding incidence angle theta;
c03: performing polynomial fitting on f (theta), f (R) by using a least square method, selecting different indexes N, and determining N by using the root mean square error RMSE of points participating in fitting and fitting functions as a selection basis of f (theta), f (R)1=4,N2=3,N3=3,N4=4。
As a preferable technical solution, when calculating an average value of a plurality of corrected laser echo intensity values in a03, a laser echo intensity queue is formed by arranging the plurality of corrected laser echo intensity values according to intensity, and the head 10% and the tail 10% of the queue are removed.
The invention has the following beneficial effects:
the method comprises the steps of obtaining point cloud data of the middle parts of all bamboo sections of the moso bamboo through a three-dimensional laser scanner, and correcting distance measurement values and incidence angles of laser echo intensity values; calculating the average value of the laser echo intensity values after correction as the laser echo intensity value of the bamboo joint section, and obtaining a fitting model of bamboo joint numbers of bamboo stalks with different bamboo degrees and the laser echo intensity value; the bamboo degree of the moso bamboo is judged by calculating whether the laser echo intensity value of each section of the moso bamboo meets the fitting model, and the discrimination efficiency and accuracy of the bamboo degree of the moso bamboo are improved.
Of course, it is not necessary for any product in which the invention is practiced to achieve all of the above-described advantages at the same time.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings used in the description of the embodiments will be briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art that other drawings can be obtained according to the drawings without creative efforts.
FIG. 1 is a flowchart of a method for determining bamboo degree based on laser echo intensity according to the present invention;
FIG. 2 is a flow chart illustrating the distance measurement and incident angle correction of the laser echo intensity value in A02 according to the present invention;
FIG. 3 is a flow chart of the present invention B03 for calculating model parameters and correcting the laser echo intensity values;
FIG. 4 is a distribution diagram of the three-dimensional laser scanner in a rectangular plot in accordance with the present invention;
FIG. 5 is a schematic diagram of extracting point cloud data from a window of 2cm by 2cm at the middle of a bamboo joint in the present invention;
FIG. 6 is a graph showing the fitting results of the fitting model of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Referring to fig. 1, the present invention is a method for determining bamboo degree of moso bamboo based on laser echo intensity, including the following steps:
a00: referring to fig. 4, a sample plot of a rectangular area is established, a three-dimensional laser scanner is erected with the periphery and the center of the sample plot as a measuring station, and overlapping parts are uniformly scanned during scanning; specifically, the above-mentioned pattern is 20m × 20m, and 188 moso bamboos are cumulatively scanned; arranging 3 plane targets in the sample plot, and carrying out station measurement splicing through the plane targets, wherein the splicing error is controlled within 1 mm; in addition, each wood ruler is manually checked, and the information of the age, the diameter of breast height (1.3 m) and the height of bamboo of each moso bamboo is recorded;
a01: selecting 8 plants of one to seven-year-old moso bamboos, cutting the parts from the lower branches of the moso bamboos to the bamboo stems on the ground, and extracting point cloud data from the middle part of each section of bamboo joint in a window of 2cm to 2cm, as shown in figure 5; the point cloud data consists of a plurality of laser echo intensity values; when the laser echo intensity value correcting device is actually used, when the average value of a plurality of corrected laser echo intensity values is calculated, a laser echo intensity queue is formed by arranging the plurality of corrected laser echo intensity values according to the intensity, and 10% of the head of the queue and 10% of the tail of the queue are removed so as to reduce errors;
a02: acquiring point cloud data by adopting a Fence function in Leica matching software Cyclone, and correcting a distance measurement value and an incident angle of a laser echo intensity value;
a03: calculating the average value of the corrected laser echo intensity values as the laser echo intensity value of the section of bamboo joint;
a04: obtaining fitting models of bamboo joint numbers of bamboo stalks with different bamboo degrees and laser echo intensity values:
y1=0.48x2-17.937x-1014.5,R2=0.9611
y2=-0.4425x2+12.983x-1157.8,R2=0.9231
y3=-0.0547x2-4.2539x-1050.8,R2=0.9546
y4=0.0016x6-0.0849x5-19.105x3+109.77x2-315.42x-680.66,R2=0.9466
wherein, y1、y2、y3、y4Respectively fitting laser echo strength values of the bamboo stalks of one to four degrees; x is the bamboo knot number; furthermore, the fitting order is determined while ensuring the fitting model correlation coefficient R2Above 0.95, the fitting result of the function model is accurate, and the fitting result is shown in figure 6.
As a preferred embodiment, please refer to fig. 2, a02 specifically includes the following steps:
actually, the point cloud laser echo intensity represents the reflection spectrum information of the target to the laser, and the phenomenon of same object and different spectrum or same spectrum of foreign objects can occur under the influence of various factors such as the system characteristics of a scanner, atmospheric attenuation, the physical properties of the surface of the target, the scanning environment and the like. Therefore, various factors influencing the point cloud laser echo intensity data need to be corrected aiming at the phenomenon; the reflectivity of the plane target is high, so that the reflectivity of the target can be regarded as 100%, the influence of the reflectivity of an object on the echo intensity of the point cloud laser is ignored, and the influence of the reflectivity rho on the echo intensity of the laser is regarded as a fixed constant; in the ground scanning process, the influence of atmospheric attenuation on the laser echo intensity can be ignored due to the fact that the laser measuring distance is small.
B00: neglecting the influence of the reflectivity of the object and the atmospheric attenuation on the laser echo intensity value, and obtaining a simplified laser echo intensity empirical model: i ═ Gf (θ) f (r);
wherein, I is the original laser echo intensity value, G is a constant, f (theta) is an incidence angle correction function, and f (R) is a distance correction function;
b01: fixing the incidence angle theta of the three-dimensional laser scanner, selecting a 3In plane target as a scanning object, changing the distance between the target and the scanner by a fixed step length, keeping the center of an electric detector of the three-dimensional laser scanner and the center of the target at the same height for scanning, and obtaining the relation between different distance measurement and laser echo intensity values under the fixed incidence angle; specifically, the target plane is perpendicular to the ground, and the incidence angle theta is kept to be 0; when R is less than or equal to 10m, the step length is 1 m; when R is more than 10m, the step length is 5 m;
b02: fixing the incident distance R of the three-dimensional laser scanner, rotating the target plane by a fixed step length, changing the incident angle of the scanned target, and acquiring the relationship between different incident angles and laser echo intensity values at a fixed distance; specifically, the distance between the three-dimensional laser scanner and the target is 10 m; when theta is less than or equal to 40 degrees, the step length is 10 degrees; theta is more than 40 degrees, and the step length is 5 degrees;
b03: and establishing a piecewise fitting function between the laser echo intensity value and the ranging and incidence angles, calculating model parameters, and correcting the laser echo intensity value.
As a preferred embodiment, please refer to fig. 3, B03 specifically includes the following steps:
in fact, for the Leica C05 ground three-dimensional laser scanner, the laser echo intensity of the point cloud appears as two distinct peak regions in the range of 0m-30m, roughly in a linear relationship within 0 ° -50 °, drops significantly after 50 °, and after 50 °, the laser echo intensity rapidly decreases with increasing incidence angle. Compared with the ranging factor, the influence of large-angle measurement on the laser echo intensity is more obvious. Therefore, in the test process, the relation between the laser echo intensity and two factors, namely the laser ranging value and the incident angle, is considered to establish a piecewise fitting function relation, so that model parameters are calculated, and the point cloud laser echo intensity value is corrected;
c00: selecting point cloud data with the distance of 8-10m and the incident angle of 40-50 degrees to perform polynomial fitting, and acquiring an extreme point of a piecewise fitting function as a critical point; wherein the distance correction critical point is 9.9m, and the angle correction critical point is 45 degrees;
c01: obtaining a laser echo intensity value correction formula containing distance and incident angle:
wherein, IsFor the corrected laser echo intensity value, (X, Y, Z) is the three-dimensional coordinate of the point cloud space, RsFor a reference distance (10 m in the present embodiment), θsIs a reference angle (0 degree in the present embodiment), alphaiAnd betaiAre all polynomial coefficients, | n | is the modulus of the normal vector; i is the original laser echo intensity; k is a system constant, theta is an incident angle, and R is an incident distance;
c02: obtaining a point cloud normal vector n ═ (n) according to a K nearest neighbor domain method1,n2,n3) Solving the cosine value of the point incidence angle according to the cosine law, and solving the inverse cosine function of the cosine value to obtain the corresponding incidence angle theta;
c03: carrying out polynomial fitting on f (theta), f (R) by using a least square method, selecting different indexes N, and using root mean square errors RMSE of points participating in fitting and fitting functions as selection basis of f (theta), f (R); see Table 1 for details, determination of N1=4,N2=3,N3=3,N 44, avoiding overfitting while ensuring f (theta), f (r) to have higher fitting accuracy; in particular, the method comprises the following steps of,
table 1 shows the piecewise polynomial function order selection and fitting function root mean square error:
table 2 shows the standard deviation of the target laser echo intensities before and after correction:
the average laser echo intensities of the target after the distance and the incidence angle are corrected tend to be consistent, the standard deviation of the average value of the laser echo intensities of the point cloud before and after the correction is calculated is shown in table 2, the standard deviation before and after the correction is obviously reduced, and the correction result is proved to be practical and effective.
When the device is actually used, point cloud data of the middle parts of all bamboo sections of the moso bamboo are obtained through the three-dimensional laser scanner, and distance measurement values and incident angle correction are carried out on laser echo intensity values; calculating the average value of the laser echo intensity values after correction as the laser echo intensity value of the bamboo joint section, and obtaining a fitting model of bamboo joint numbers of bamboo stalks with different bamboo degrees and the laser echo intensity value; the bamboo degree of the moso bamboo is judged by calculating whether the laser echo intensity value of each section of the moso bamboo meets the fitting model, and the discrimination efficiency and accuracy of the bamboo degree of the moso bamboo are improved.
It should be noted that, in the above system embodiment, each included unit is only divided according to functional logic, but is not limited to the above division as long as the corresponding function can be implemented; in addition, specific names of the functional units are only for convenience of distinguishing from each other, and are not used for limiting the protection scope of the present invention.
In addition, it is understood by those skilled in the art that all or part of the steps in the method for implementing the embodiments described above may be implemented by a program instructing associated hardware, and the corresponding program may be stored in a computer-readable storage medium.
The preferred embodiments of the invention disclosed above are intended to be illustrative only. The preferred embodiments are not intended to be exhaustive or to limit the invention to the precise embodiments disclosed. Obviously, many modifications and variations are possible in light of the above teaching. The embodiments were chosen and described in order to best explain the principles of the invention and the practical application, to thereby enable others skilled in the art to best utilize the invention. The invention is limited only by the claims and their full scope and equivalents.
Claims (7)
1. The moso bamboo degree judging method based on the laser echo intensity is characterized by comprising the following steps of:
a00: establishing a sample plot of a rectangular area, and erecting a three-dimensional laser scanner by taking the periphery and the center of the sample plot as a measuring station;
a01: selecting 8 plants of one to seven-year-old moso bamboos, cutting the parts from the lower branches of the moso bamboos to the bamboo stems on the ground, and extracting point cloud data from the middle part of each section of bamboo joint in a window of 2cm to 2 cm; the point cloud data consists of a plurality of laser echo intensity values;
a02: correcting the distance measurement value and the incident angle of the laser echo strength value;
a03: calculating the average value of the corrected laser echo intensity values as the laser echo intensity value of the section of bamboo joint;
a04: obtaining fitting models of bamboo joint numbers of bamboo stalks with different bamboo degrees and laser echo intensity values:
y1=0.48x2-17.937x-1014.5
y2=-0.4425x2+12.983x-1157.8
y3=-0.0547x2-4.2539x-1050.8
y4=0.0016x6-0.0849x5-19.105x3+109.77x2-315.42x-680.66
wherein, y1、y2、y3、y4Respectively fitting laser echo strength values of the bamboo stalks of one to four degrees; and x is a bamboo knot number.
2. The method of claim 1, wherein a00 further comprises placing a 3-plane target in the sample chamber and performing station-finding stitching through the plane target.
3. The method for determining the bamboo degree based on the laser echo intensity of claim 2, wherein a02 specifically comprises the following steps:
b00: neglecting the influence of the reflectivity of the object and the atmospheric attenuation on the laser echo intensity value, and obtaining a simplified laser echo intensity empirical model: i ═ Gf (θ) f (r);
wherein, I is the original laser echo intensity value, G is a constant, f (theta) is an incidence angle correction function, and f (R) is a distance correction function;
b01: fixing the incidence angle theta of the three-dimensional laser scanner, selecting a 3In plane target as a scanning object, changing the distance between the target and the scanner by a fixed step length, keeping the center of an electric detector of the three-dimensional laser scanner and the center of the target at the same height for scanning, and obtaining the relation between different distance measurement and laser echo intensity values under the fixed incidence angle;
b02: fixing the incident distance R of the three-dimensional laser scanner, rotating the target plane by a fixed step length, changing the incident angle of the scanned target, and acquiring the relationship between different incident angles and laser echo intensity values at a fixed distance;
b03: and establishing a piecewise fitting function between the laser echo intensity value and the ranging and incidence angles, calculating model parameters, and correcting the laser echo intensity value.
4. The method for determining the degree of bamboo by using the echo intensity of laser according to claim 3, wherein the target plane in B01 is perpendicular to the ground, and the incident angle θ is kept to be 0; when R is less than or equal to 10m, the step length is 1 m; when R is more than 10m, the step length is 5 m.
5. The method of claim 3, wherein the three-dimensional laser scanner in B02 is mounted at a distance of 10m from the target; when theta is less than or equal to 40 degrees, the step length is 10 degrees; theta > 40 deg. and the step size is 5 deg..
6. The method for determining the bamboo degree based on the laser echo intensity of claim 3, wherein B03 specifically comprises the following steps:
c00: selecting point cloud data with the distance of 8-10m and the incident angle of 40-50 degrees to perform polynomial fitting, and acquiring an extreme point of a piecewise fitting function as a critical point; wherein the distance correction critical point is 9.9m, and the angle correction critical point is 45 degrees;
c01: obtaining a laser echo intensity value correction formula containing distance and incident angle:
wherein, IsFor the corrected laser echo intensity value, (X, Y, Z) is the three-dimensional coordinate of the point cloud space, RsTo reference the distance, θsAs a reference angle, αiAnd betaiAre all polynomial coefficients, | n | is the modulus of the normal vector; i is the original laser echo intensity; k is a system constant, theta is an incident angle, and R is an incident distance;
c02: obtaining a point cloud normal vector n ═ (n) according to a K nearest neighbor domain method1,n2,n3) Solving the cosine value of the point incidence angle according to the cosine law, and solving the inverse cosine function of the cosine value to obtain the corresponding incidence angle theta;
c03: performing polynomial fitting on f (theta), f (R) by using a least square method, selecting different indexes N, and determining N by using the root mean square error RMSE of points participating in fitting and fitting functions as a selection basis of f (theta), f (R)1=4,N2=3,N3=3,N4=4。
7. The method for determining the degree of bamboo or bamboo on the basis of laser echo intensities according to claim 1, 3, 4, 5 or 6, wherein when the average value of the plurality of corrected laser echo intensity values is calculated in A03, the plurality of corrected laser echo intensity values are arranged according to the intensity to form a laser echo intensity queue, and the head 10% and the tail 10% of the queue are removed.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110291159.3A CN113109832A (en) | 2021-03-18 | 2021-03-18 | Bamboo degree judgment method based on laser echo intensity |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110291159.3A CN113109832A (en) | 2021-03-18 | 2021-03-18 | Bamboo degree judgment method based on laser echo intensity |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113109832A true CN113109832A (en) | 2021-07-13 |
Family
ID=76711854
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110291159.3A Pending CN113109832A (en) | 2021-03-18 | 2021-03-18 | Bamboo degree judgment method based on laser echo intensity |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113109832A (en) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7411688B1 (en) * | 2006-03-17 | 2008-08-12 | Arius3D Inc. | Method and system for laser intensity calibration in a three-dimensional multi-color laser scanning system |
CN109945838A (en) * | 2019-04-19 | 2019-06-28 | 福建农林大学 | Three-dimensional laser point cloud acquisition methods and auxiliary device suitable for mao bamboo woods |
CN110346782A (en) * | 2019-05-31 | 2019-10-18 | 华东师范大学 | A kind of correcting method of long range ground three-dimensional laser radar echo strength data |
CN110415259A (en) * | 2019-07-30 | 2019-11-05 | 南京林业大学 | A kind of shade tree point cloud recognition methods based on laser reflection intensity |
-
2021
- 2021-03-18 CN CN202110291159.3A patent/CN113109832A/en active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7411688B1 (en) * | 2006-03-17 | 2008-08-12 | Arius3D Inc. | Method and system for laser intensity calibration in a three-dimensional multi-color laser scanning system |
CN109945838A (en) * | 2019-04-19 | 2019-06-28 | 福建农林大学 | Three-dimensional laser point cloud acquisition methods and auxiliary device suitable for mao bamboo woods |
CN110346782A (en) * | 2019-05-31 | 2019-10-18 | 华东师范大学 | A kind of correcting method of long range ground three-dimensional laser radar echo strength data |
CN110415259A (en) * | 2019-07-30 | 2019-11-05 | 南京林业大学 | A kind of shade tree point cloud recognition methods based on laser reflection intensity |
Non-Patent Citations (10)
Title |
---|
XINYUE TAO 等: "Phytolith sizes and assemblages differentiate genera and ecotypes of woody bamboos in subtropical Southwest China", 《REVIEW OF PALAEOBOTANY AND PALYNOLOGY》, 31 December 2020 (2020-12-31), pages 1 - 7 * |
廖晓和 等: "基于车载点云数据的树木提取与分", 《测绘通报》, no. 11, 31 December 2020 (2020-12-31), pages 163 - 166 * |
张仓皓 等: "毛竹立竹度无人机遥感识别有效高度的研究", 《遥感技术与应用》, vol. 35, no. 6, 31 December 2020 (2020-12-31), pages 1436 - 1446 * |
程小龙 等: "基于分段多项式模型的地面三维激光扫描激光强度改正", 《激光与光电子学进展》, 31 December 2017 (2017-12-31), pages 112802 - 1 * |
苏德添 等: "基于毛竹冠层LiDAR点云密度测算其蓄积量", 《中国激光》, vol. 47, no. 4, 30 April 2020 (2020-04-30), pages 1 - 8 * |
蔡越等: "基于激光回波强度判别毛竹年龄", 《中国激光》, no. 01, 31 January 2018 (2018-01-31), pages 1 - 9 * |
谭凯 等: "TLS强度数据的入射角及距离效应改正方法", 《武汉大学学报·信息科学版》, vol. 42, no. 2, 28 February 2017 (2017-02-28), pages 223 - 228 * |
赵谦,西安电子科技大学出版社, vol. 1, pages: 169 * |
邓英英 等: "天目山近自然毛竹纯林的竹秆空间结构特征", 《浙江农林大学学报》, vol. 28, no. 2, 31 December 2011 (2011-12-31), pages 173 - 179 * |
陈锦 等: "基于地面三维激光扫描强度数据的潮滩表层含水量估算", 《地球信息科学》, vol. 22, no. 2, 29 February 2020 (2020-02-29), pages 290 - 297 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wu et al. | An accurate skeleton extraction approach from 3D point clouds of maize plants | |
CN111104858B (en) | Large-scale crop phenology extraction method based on morphological model method | |
CN110969654A (en) | Corn high-throughput phenotype measurement method and device based on harvester and harvester | |
CN102331244A (en) | Method of acquiring leaf area index (LAI) of rice canopy by using common digital camera | |
CN110598532B (en) | Tree pest and disease damage monitoring system and method | |
CN105806266B (en) | Trees canopy leaf area computational methods based on laser scanning data | |
CN114694047A (en) | Corn sowing quality evaluation method and device | |
CN106683092B (en) | Device and method for measuring and calculating crown canopy density of blueberries | |
CN108280410A (en) | One kind being based on binary-coded crops recognition methods and system | |
Li et al. | Three-dimensional reconstruction and phenotype measurement of maize seedlings based on multi-view image sequences | |
CN107292957B (en) | Three-dimensional reconstruction method and device for crop canopy | |
CN113109832A (en) | Bamboo degree judgment method based on laser echo intensity | |
CN109631821A (en) | A kind of method of quick nondestructive measurement bamboo leaves area | |
Ma et al. | A method for calculating and simulating phenotype of soybean based on 3D reconstruction | |
CN117178769A (en) | Automatic garden plant maintenance method and system | |
CN105403177B (en) | A kind of detection method of the leaf area index based on video detection | |
CN116402879A (en) | Corn ear position measuring method, device, system and storage medium | |
Guo et al. | Apple tree canopy leaf spatial location automated extraction based on point cloud data | |
CN116258844A (en) | Rapid and accurate identification method for phenotype character of cotton leaf | |
CN115808668A (en) | Method for extracting number of field wheat tillers based on mean shift algorithm of voxel interpolation | |
CN112215184B (en) | Method for detecting yield of camellia oleifera fruit trees based on three-dimensional laser scanner | |
CN108804854A (en) | The tree management method for parameter estimation and its system that one planting fruit-trees machinery is pruned | |
CN114494689A (en) | Identification method of tomato drought stress | |
CN114049390A (en) | Wheat seedling planting density measuring device and method based on machine vision | |
CN106767717A (en) | A kind of rice maturity GVG based on unmanned plane is than sampling and detection method |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210713 |
|
RJ01 | Rejection of invention patent application after publication |