CN112415537A - Model for depicting forest photosynthetic effective radiation distribution by using three-dimensional point cloud data - Google Patents
Model for depicting forest photosynthetic effective radiation distribution by using three-dimensional point cloud data Download PDFInfo
- Publication number
- CN112415537A CN112415537A CN202010396373.0A CN202010396373A CN112415537A CN 112415537 A CN112415537 A CN 112415537A CN 202010396373 A CN202010396373 A CN 202010396373A CN 112415537 A CN112415537 A CN 112415537A
- Authority
- CN
- China
- Prior art keywords
- par
- point
- forest
- canopy
- dimensional
- 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.)
- Granted
Links
- 238000009826 distribution Methods 0.000 title claims abstract description 40
- 230000005855 radiation Effects 0.000 title claims abstract description 21
- 230000000243 photosynthetic effect Effects 0.000 title description 5
- 230000035515 penetration Effects 0.000 claims description 17
- 238000000034 method Methods 0.000 claims description 16
- 230000005540 biological transmission Effects 0.000 claims description 6
- 238000011160 research Methods 0.000 abstract description 6
- 230000035790 physiological processes and functions Effects 0.000 abstract 1
- 238000005259 measurement Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 2
- 238000012821 model calculation Methods 0.000 description 2
- 230000029553 photosynthesis Effects 0.000 description 2
- 238000010672 photosynthesis Methods 0.000 description 2
- 239000011148 porous material Substances 0.000 description 2
- 101100136092 Drosophila melanogaster peng gene Proteins 0.000 description 1
- 230000002457 bidirectional effect Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000011888 foil Substances 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000002310 reflectometry Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 230000005068 transpiration Effects 0.000 description 1
- 238000009827 uniform distribution Methods 0.000 description 1
- 239000011800 void material Substances 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
- 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
- G01S17/89—Lidar systems specially adapted for specific applications for mapping or imaging
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Electromagnetism (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Optical Radar Systems And Details Thereof (AREA)
- Image Processing (AREA)
Abstract
本发明提供了一种利用三维点云数据刻画森林光合有效辐射(PAR)分布的模型，属于植被遥感建模的研究领域。针对目前未能在精细尺度上利用激光雷达数据有效估算森林冠层PAR三维分布的难题，本发明以多站获取的森林地面激光雷达数据为基础，提出了一种森林PAR三维时空分布估算模型，模型分开考虑太阳直射和散射PAR在森林中的传输过程，同时顾及元素角度的影响。本发明所构建的模型可在精细尺度上定量刻画森林中PAR的三维时空分布，精度较高，有助于森林生理过程的研究。
The invention provides a model for depicting the distribution of photosynthetically active radiation (PAR) in a forest by using three-dimensional point cloud data, which belongs to the research field of vegetation remote sensing modeling. Aiming at the problem that the three-dimensional distribution of forest canopy PAR cannot be effectively estimated by using lidar data on a fine scale at present, the present invention proposes a forest PAR three-dimensional spatiotemporal distribution estimation model based on the forest ground lidar data obtained from multiple stations. The model separately considers direct solar and scattered PAR transport in the forest, while taking into account the effects of element angles. The model constructed by the invention can quantitatively describe the three-dimensional space-time distribution of PAR in the forest on a fine scale, with high precision, which is helpful for the research on the physiological process of the forest.
Description
One, the technical field
The invention relates to a model for depicting three-dimensional space-time distribution of forest photosynthetically active radiation on a fine scale through three-dimensional point cloud data, in particular to a method for accurately estimating three-dimensional distribution of forest photosynthetically active radiation at any time on the fine scale by acquiring three-dimensional point cloud data of a forest canopy through a ground laser radar technology, belonging to the research field of vegetation remote sensing modeling.
Second, background Art
The Photosynthetically Active Radiation (PAR) refers to the part of solar radiation with the wavelength of 400-700nm which can drive vegetation to carry out photosynthesis, and is an important environmental factor influencing the photosynthesis and transpiration process of the vegetation. The difference of the PAR received by different positions of the forest canopy is obvious, the difference changes along with the change of the position of the sun, and meanwhile, the photosynthetic speed of the leaves has a nonlinear characteristic on the PAR response, so that the three-dimensional space-time distribution of the photosynthetically active radiation of the forest canopy has important significance for researching the material and energy exchange process between the forest and the atmosphere.
The PAR of the forest canopy at any time and specific position can be measured by using a quantum sensor or an optical instrument, however, when the random distribution degree of canopy elements is high, the cost for obtaining the three-dimensional distribution of the PAR is high, and meanwhile, the distribution of the canopy elements can be influenced in the measurement process, so that the method is usually used for measuring the distribution of the PAR below the canopy and is less used for obtaining the three-dimensional space-time distribution of the PAR of the forest canopy. The three-dimensional space-time distribution of the PAR of the forest canopy is closely related to the canopy structure, and in order to overcome the limitation of ground measurement, some researches describe the three-dimensional space-time distribution of light on the forest canopy by constructing a forest canopy three-dimensional model and further utilizing methods such as ray tracing and the like. The three-dimensional model of the forest canopy can be divided into four types according to the complexity. The simplest type is that a forest canopy is assumed to be a mixed medium layer with elements distributed randomly, and the penetration coefficient of light is calculated by utilizing a Leaf Area Index (LAI) through a Bell's law theory, so that the energy value of the corresponding position of the canopy is obtained. For example, Myneni et al, in "A procedure for interpreting the radiation region of defined and truncated foil sites, part I. The mathematical expression of bell's law is relatively simple and widely used, but only for forest canopies with a uniform distribution of elements. The second method is to assume the forest canopy as a geometric body, e.g. a sphere, ellipsoid. For example, Chen et al in "A four-scale bidirectional reflectance model based on selectivity architecture" describes the reflectivity information of the canopy by assuming the canopy is a cone. The third method is to consider information such as leaf area density and leaf angle distribution in the voxel and the surface element, assuming that the forest canopy is composed of the voxel or the surface element. For example, in Huanghua and other countries, the RAPID principle and application of a three-dimensional remote sensing mechanism model are that a forest canopy is assumed to be composed of surface elements, and the radiation distribution in the canopy is described by using a radiometric method. The fourth method is to simulate the three-dimensional scene of the canopy with a fine structure by a computer, and the model can obtain more accurate PAR three-dimensional distribution. For example, Zhongjie in "simulation study of the three-dimensional distribution of direct photosynthetic active radiation based on virtual plant canopy", i.e. the three-dimensional distribution of the direct PAR of a single wood is calculated by simulating its real three-dimensional structure. However, as the forest canopy is complex and changeable, high cost and time are required for constructing a real and detailed forest three-dimensional scene. The three-dimensional forest structure model is still different from the real forest canopy structure in the nature.
The laser radar technology (LiDAR) with the three-dimensional drawing capability can acquire the real three-dimensional structure of the natural canopy and provides a foundation for describing the three-dimensional space-time distribution of the PAR of the forest canopy. The most common method is by voxelizing the canopy point cloud, each voxel calculating its PAR value using bell's law. To simplify the calculation, the PAR value is usually calculated using the path length instead of the LAI. The penetration factor of light at different locations in the canopy is calculated using the path length of light in the canopy by Peng in a "Modeling spatial patterns of understory light intensity using air laser scanner (LiDAR)". Bell's Law assumes that the canopy elements are randomly distributed, so voxel size is generally more difficult to determine. The second method is to estimate the penetration coefficient of the light at different locations of the canopy using the number of points or energy information of the points of the LiDAR data. For example, Bode in "Subcanopy Solar Radiation model: the method comprises the steps of Predicting the specific radiation of a crown layer, estimating the penetration coefficient of light at different positions of the crown layer by using the difference between the total number of emission points and the number of interception points of the crown layer, and estimating the penetration coefficient of the light by Chen at the position of "Sensitivity of direct radiation gap fraction from air turbine for laser to laser and reflection characteristics" by using the energy information emitted by LiDAR and the energy information returned.
However, existing studies are primarily based on Airborne Lidar (ALS) data to estimate the three-dimensional spatiotemporal distribution of the sub-canopy PAR using bell's law theory. First, bell's law only applies to the case where elements are randomly distributed, with limitations. Secondly, the ALS data points are of low density and do not accurately describe the three-dimensional spatiotemporal distribution of the forest canopy PAR. Meanwhile, because the photosynthetic efficiency of the vegetation is higher under the condition of scattering PAR than under the condition of direct incidence PAR, the transmission process of the direct solar radiation and the scattering PAR in the canopy needs to be considered separately, but many researches do not consider the factor. In addition, the PAR of the element surface of the canopy is influenced by the included angle between the normal vector of the element and the incident direction of the sun, the smaller the included angle is, the larger the element light-intercepting area is, and the larger the surface PAR is, but the research for describing the three-dimensional distribution of the PAR of the canopy by using LiDAR data at present cannot consider the influence of the real element angle distribution of the canopy.
Third, the invention
The invention aims to provide a model for estimating the three-dimensional space-time distribution of the PAR of a forest canopy on a fine scale by using ground laser radar (TLS) data. The principle of the invention is as follows:
the solar incident PAR is composed primarily of direct and diffuse components. The model will consider the transmission of the direct and scattered PAR inside the canopy separately. The penetration coefficients of direct incidence and scattering PAR at different positions of the canopy are the key for estimating the three-dimensional space-time distribution of the forest PAR, the direct incidence and scattering PAR of each point can be obtained by calculating the penetration coefficient of the direct incidence and scattering PAR of each point in TLS data and further utilizing the measured incident direct incidence and scattering PAR values of the sun, and finally the three-dimensional space-time distribution result of the forest canopy PAR is obtained.
The technical scheme of the invention mainly comprises the following steps:
(1) firstly, a TLS scanner is used for acquiring complete point cloud data of a forest canopy in a multi-station mode. Namely, the TLS and the target are erected at different positions of the canopy, and the station number of the scanner is adjusted according to the density of the sample. And finally, splicing the TLS data of different stations through the target to obtain complete three-dimensional point cloud data of the canopy. The values of direct solar radiation and the scattered PAR were measured with a BF5 radiometer.
(2) The direct PAR for each point is estimated. Assuming that each point in the canopy corresponds to a micro-plane, when the micro-plane is perpendicular to the incident direction of the sun, the PAR value of the surface can be expressed as the product of the transmission coefficient of the incident direct solar PAR and the direct solar PAR of the point. Therefore, the determination of the direct PAR penetration coefficient of each point is very critical. Assuming direct incident solar radiation PAR as a parallel light source, depending on the direction of the sunThe point cloud is rotated so that the solar incident rays are parallel to the-z axis. The rotated point cloud data is subjected to three-dimensional meshing, and the length and the width (V) of the mesh0) Greater than the canopy average dot spacing (fig. 1 a). Searching for a point q to be solvedj(x0，y0，z0) The point set of the grid (FIG. 1b) is higher than z0Is projected onto the XY plane and rasterized (fig. 1 c). The grid size is the mean point spacing (V)1). The grid without points is identified as pores, the proportion of which is the penetration coefficient of the direct PAR at that point. And then obtaining the PAR value of the micro-plane perpendicular to the incident direction of the sun by utilizing the penetration coefficient of the incident direct solar radiation PAR and the direct solar radiation PAR at the point. When an included angle exists between the micro plane of each point and the incident direction of the sun, the normal vector of each point is calculated by using a method of a near point structural surface, namely, 6 points around the point to be solved are searched to calculate the covariance matrix of the point to be solved, and then the eigenvalue and the eigenvector of the point are calculated, wherein the vector corresponding to the minimum eigenvalue is the normal vector of the point. And calculating a cosine value of an included angle between a normal vector of the point to be solved and the incident direction of the sun, and multiplying the cosine value by the calculated direct incidence PAR to obtain a final direct incidence PAR result of each point.
(3) The scatter PAR for each spot was estimated. The scattering PAR for each point in the canopy can be expressed as the product of the solar incident scattering PAR and the transmission coefficient of the scattering PAR for that point. Therefore, the determination of the scattering PAR penetration coefficient of each point is very critical. The solar incident scattering PAR is assumed to be isotropic. And translating the TLS data to enable the point to be solved to be a central point, converting a coordinate system into a polar coordinate system, and setting an inclination angle interval and an azimuth angle interval to be V to divide the upper hemisphere of the canopy into [90/V ] × [360/V ] conical bodies. When there is no point in the cone, it is identified as a void. The ratio of the pores is used as the transmission coefficient of the scattering PAR of the point to be determined. The product of the penetration coefficients of the solar incident scattering PAR and the scattering PAR is the scattering PAR of the point to be determined.
(4) And estimating the three-dimensional distribution of the PAR of the forest at a specific moment. The sum of the direct PAR and the scattered PAR of each point in the TLS data is the total PAR of each point, and thus the three-dimensional spatial distribution of the corresponding temporal canopy total PAR, direct PAR, and scattered PAR can be obtained.
Description of the drawings
FIG. 1 direct PAR penetration coefficient calculation
(a) For gridding TLS data after rotating according to solar direction
(b) To select the grid where the point to be solved is located
(c) Projecting the points on an XY plane and rasterizing
FIG. 2 forest canopy TLS data and results of three-dimensional distribution of PAR every two hours
(a) TLS data for forest canopy
(b) - (f) canopy in 2017, 7, 27, 8: 00-16: 00 results of PAR three-dimensional distribution every two hours, the brighter the color, the higher the PAR value, the minimum value of 0. mu. mol. m-2·s-1Maximum value of 2036.6 μmol · m-2·s-1
Fifth, detailed description of the invention
The invention is further described below by way of specific examples:
and (2) acquiring forest sample point cloud data (figure 2a) by using a ground laser radar scanner according to the technical scheme step (1). And (3) setting the sampling interval of TLS to be 1cm at a position of 10m according to the density of the sample, and totally setting 5 stations, one station at the center and one station at each of four corners of the sample. Three positions are selected to place targets for data stitching processing. And after the data are acquired, splicing by using Cyclone software, and cutting the center area of the sample. Noise points in the data were manually removed. The TLS data is rotated according to the reference point such that the y-axis is in the due north direction. The values of the direct solar incident PAR and the scattered PAR were measured with a BF5 radiometer.
Calculating direct PAR of each point according to the step (2) of the technical scheme, calculating scattering PAR of each point according to the step (3) of the technical scheme, and finally obtaining three-dimensional space-time distribution of total PAR of the canopy according to the step (4) of the technical scheme. The accuracy of the model was verified by comparing the model calculation results with the radiometer measurement results for 9 specific locations (table 1). The results of the three-dimensional distribution of the PAR every two hours of the canopy are finally obtained (fig. 2b-2 f).
TABLE 1 comparison of model calculations with radiometer measurements
Claims (2)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010396373.0A CN112415537B (en) | 2020-05-08 | 2020-05-08 | Model for depicting forest photosynthetic effective radiation distribution by using three-dimensional point cloud data |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010396373.0A CN112415537B (en) | 2020-05-08 | 2020-05-08 | Model for depicting forest photosynthetic effective radiation distribution by using three-dimensional point cloud data |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112415537A true CN112415537A (en) | 2021-02-26 |
| CN112415537B CN112415537B (en) | 2022-07-12 |
Family
ID=74844058
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010396373.0A Active CN112415537B (en) | 2020-05-08 | 2020-05-08 | Model for depicting forest photosynthetic effective radiation distribution by using three-dimensional point cloud data |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112415537B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116758006A (en) * | 2023-05-18 | 2023-09-15 | 广州广检建设工程检测中心有限公司 | Scaffold quality detection method and device |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102708254A (en) * | 2012-05-18 | 2012-10-03 | 福州大学 | Method for analyzing intercept ability of canopy on photosynthetic active radiation based on virtual plant |
| CN107831497A (en) * | 2017-09-26 | 2018-03-23 | 南京大学 | A kind of method that forest building-up effect is quantitatively portrayed using three dimensional point cloud |
| CN108871561A (en) * | 2018-04-25 | 2018-11-23 | 南京大学 | A kind of method of continuous observation rice canopy absorbed photosynthtic active radiation |
| CN110006900A (en) * | 2019-04-23 | 2019-07-12 | 中国农业科学院棉花研究所 | Method, device and equipment for obtaining the interception rate of photosynthetically active radiation in plant canopy |
-
2020
- 2020-05-08 CN CN202010396373.0A patent/CN112415537B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102708254A (en) * | 2012-05-18 | 2012-10-03 | 福州大学 | Method for analyzing intercept ability of canopy on photosynthetic active radiation based on virtual plant |
| CN107831497A (en) * | 2017-09-26 | 2018-03-23 | 南京大学 | A kind of method that forest building-up effect is quantitatively portrayed using three dimensional point cloud |
| CN108871561A (en) * | 2018-04-25 | 2018-11-23 | 南京大学 | A kind of method of continuous observation rice canopy absorbed photosynthtic active radiation |
| CN110006900A (en) * | 2019-04-23 | 2019-07-12 | 中国农业科学院棉花研究所 | Method, device and equipment for obtaining the interception rate of photosynthetically active radiation in plant canopy |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116758006A (en) * | 2023-05-18 | 2023-09-15 | 广州广检建设工程检测中心有限公司 | Scaffold quality detection method and device |
| CN116758006B (en) * | 2023-05-18 | 2024-02-06 | 广州广检建设工程检测中心有限公司 | Scaffold quality detection method and device |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112415537B (en) | 2022-07-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Gastellu-Etchegorry et al. | Simulation of satellite, airborne and terrestrial LiDAR with DART (I): Waveform simulation with quasi-Monte Carlo ray tracing | |
| CN102914501B (en) | Method for calculating extinction coefficients of three-dimensional forest canopy by using laser-point cloud | |
| Li et al. | VBRT: A novel voxel-based radiative transfer model for heterogeneous three-dimensional forest scenes | |
| CN111563962A (en) | A Remote Sensing Image Simulation Method Based on Geometric Radiation Integrated Sampling | |
| CN104049259A (en) | Lidar three-dimensional imaging system based on virtual instrument | |
| CN108414998A (en) | A kind of laser satellite altitude meter echo waveform analog simulation method and equipment | |
| CN107993281B (en) | Method for simulating optical characteristics of space target visible light | |
| CN112068153A (en) | A method for estimating canopy gap ratio based on ground-based lidar point cloud | |
| CN105371789A (en) | Method for utilizing aviation laser point cloud to calculate effective leaf area index | |
| Ma et al. | Characterizing the three-dimensional spatiotemporal variation of forest photosynthetically active radiation using terrestrial laser scanning data | |
| Ahmadibeni et al. | Aerial and ground vehicles synthetic SAR dataset generation for automatic target recognition | |
| Lin et al. | Reconstruction of a large-scale realistic three-dimensional (3-D) mountain forest scene for radiative transfer simulations | |
| JP5137039B2 (en) | Thermal radiation energy simulation apparatus and method | |
| CN114663786B (en) | Stand radiation flux calculating method based on point cloud data and computer graphics | |
| CN112415537A (en) | Model for depicting forest photosynthetic effective radiation distribution by using three-dimensional point cloud data | |
| Kawata et al. | Multiple scattering analysis of airborne POLDER image data over the sea | |
| Yin et al. | Simulating satellite waveform Lidar with DART model | |
| CN108318458B (en) | Method for measuring outdoor typical feature pBRDF (binary RDF) suitable for different weather conditions | |
| Li et al. | Achieve accurate recognition of 3D point cloud images by studying the scattering characteristics of typical targets | |
| Henriksson et al. | Optical reflectance tomography using TCSPC laser radar | |
| Wang et al. | Research on method of high-precision 3D scene optical remote sensing imaging simulation | |
| CN104933302A (en) | Method and system for realizing aviation remote sensing simulation of unmanned aerial vehicle load | |
| Wang et al. | Analysis on the inversion accuracy of LAI based on simulated point clouds of terrestrial LiDAR of tree by ray tracing algorithm | |
| Labarre et al. | MATISSE-v2. 0: new functionalities and comparison with MODIS satellite images | |
| CN115964837A (en) | Underwater complex target laser scanning point cloud simulation method based on ray tracing |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |