CN108492332A - Leaf area index real-time computing technique in a kind of forest three-dimensional scenic - Google Patents
Leaf area index real-time computing technique in a kind of forest three-dimensional scenic Download PDFInfo
- Publication number
- CN108492332A CN108492332A CN201810289091.3A CN201810289091A CN108492332A CN 108492332 A CN108492332 A CN 108492332A CN 201810289091 A CN201810289091 A CN 201810289091A CN 108492332 A CN108492332 A CN 108492332A
- Authority
- CN
- China
- Prior art keywords
- forest
- sample
- leaf area
- area index
- coordinate
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2200/00—Indexing scheme for image data processing or generation, in general
- G06T2200/04—Indexing scheme for image data processing or generation, in general involving 3D image data
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30181—Earth observation
- G06T2207/30188—Vegetation; Agriculture
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geometry (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Image Processing (AREA)
Abstract
Leaf area index real-time computing technique in a kind of forest three-dimensional scenic, belongs to forestry scientific research forest D visualized simulation technical field.Calculate sample ground forest quantity step in forest three-dimensional scenic;Calculate density of crop step;Build the density of crop and leaf area index model step;Leaf area index calculating operation step in forest three-dimensional scenic.It is an advantage of the invention that:Using this model for calculating leaf area index in arbitrary standing forest sample ground according to the density of crop in forest three-dimensional scenic, the shortcomings that overcoming the permanent sample plot of the density of crop and leaf area index in forest three-dimensional scenic only for load.The forest three-dimensional scenic density of crop being calculated can meet the real-time calculating and display of leaf area index in a wide range of scene forest.The present invention provides standing forest densitometer in forest three-dimensional scenic and calculates method, and calculate the numerical value of leaf area index in forest three-dimensional scenic by institute's structure leaf area index and the density of crop relational model on the basis of forest three-dimensional scenic is built.
Description
Technical field
The present invention relates to leaf area index real-time computing techniques in a kind of forest three-dimensional scenic, belong to forestry scientific research
Forest D visualized simulation technical field.
Background technology
With the development of computer graphics, forest D visualized simulation technology has been applied to current forest-science and has ground
The every field such as study carefully, and the importance that leaf area index is embodied as forest photosynthesis and upgrowth situation, for mould
Forest Growth under the influence of quasi- research luminous environment also has irreplaceable role.For only focusing on the forest of 3D visual effect
Illumination visual Simulation can be arranged with use environment light (such as the forest three-dimensional scenic simulation in game) and sun shade etc.
Property method simulate luminous environment in woods, but for the forestry for more focusing on stand structure and growth rhythm visualizes, in woods
The validity of leaf area index and real-time are then most important.On the other hand, in forest three-dimensional scenic leaf area index only for
It investigates and for the permanent sample plot of loading simulation, range is small, and leaf area index data are only to the sample effective, it is impossible to be used in mould
The quasi- leaf area index for calculating other sample ground in the scene, to forest three when can not meet forest three-dimensional scenic visual Simulation
Appoint the needs that leaf area index calculates in real time in equally in dimension scene.Therefore, the density of crop and leaf area index model construction
It is the basis that leaf area index calculates in real time in realization forest three-dimensional scenic.The density of crop is to influence Stand Growth in forest culture and management work
Several factors in the factor that can effectively control, largely determine the internal structure of standing forest, control and adjustment
The density of crop also becomes orest management person and researcher's questions of common interest.
D visualized simulation technology is completed by computer graphics techniques, is mainly used in virtual forest ring
Qualitative or quantitative comprehensive analysis is carried out to relevant issues such as Forest Growth and operation in border, to fast and efficiently solve complexity
Forest management planning problem, accelerate Forestry informationlization and modernization, Management offorestry pushed longer to be touched on the spot from the period
Rope develops to the good visualization business decision direction of science.
In forest D visualized simulation scene, realizes the high validity visual research of crown canopy structure, be conducive to section
The more intuitive efficient analyzing processing information of staff is ground, landscape planning and forestry scientific research can be also widely used in
Deng.
Important component of the leaf area index as crown canopy structural parameters can be to be studied on standing forest and large scale level
The structure of forest cover provides important evidence with growth.
Leaf area index calculates in forest three-dimensional scenic only obtains data by on-site inspection at present, then in forest three-dimensional
Loading section permanent sample plot in scene finally calls in permanent sample plot and surveys with showing the sample leaf area index of gained.
This method scene coverage area is small, cannot calculate the leaf area index in each piece of sample ground so that forest is three-dimensional
Leaf area index calculates critical constraints in real time in scene.
Invention content
For overcome the deficiencies in the prior art, leaf area index calculates in real time in a kind of forest three-dimensional scenic of present invention offer
Method.
The three-dimensional visible of forest illumination in the present invention is turned to the application field of the present invention;The density of crop and leaf area index
Basis of the model as this invention;Obtain target of the leaf area index as the present invention in forest three-dimensional scenic;Forest is three-dimensional
The density of crop calculates the core as this invention in real time in scene;Reality of the forest three-dimensional visualization software module as the present invention
Exampleization is applied.Leaf area index in forest three-dimensional scenic is calculated to the skill that the present invention is combined into density of crop model in real time
Art feature.
Leaf area index real-time computing technique in a kind of forest three-dimensional scenic, includes the following steps;
Calculate sample ground forest quantity step in forest three-dimensional scenic;
Calculate density of crop step;
Build the density of crop and leaf area index model step;
Leaf area index calculating operation step in forest three-dimensional scenic.
It is an advantage of the invention that:
Real-time calculate of the density of crop is the core realized leaf area index and calculated in forest three-dimensional scenic.That is,
It is the present invention that the density of crop, which is combined with leaf area index relational model with the real-time calculating of the density of crop in forest three-dimensional scenic,
Technology contents.
The relational models such as leaf area index and the density of crop or the age of stand are generally used for stand structure and light radiation rule
Qualitative and Quantitative research, the present invention is using this model for calculating arbitrary standing forest sample according to the density of crop in forest three-dimensional scenic
Leaf area index in ground overcomes the permanent sample plot of the density of crop and leaf area index in forest three-dimensional scenic only for load
Disadvantage.The forest three-dimensional scenic density of crop being calculated can meet the real-time meter of leaf area index in a wide range of scene forest
It calculates and shows.The present invention provides standing forest densitometer calculation side in forest three-dimensional scenic on the basis of forest three-dimensional scenic is built
Method, and calculate by institute's structure leaf area index and the density of crop relational model number of leaf area index in forest three-dimensional scenic
Value.
Description of the drawings
When considered in conjunction with the accompanying drawings, by referring to following detailed description, can more completely more fully understand the present invention with
And be easy to learn the advantage that many of which is adjoint, but attached drawing described herein is used to provide further understanding of the present invention,
The part of the present invention is constituted, the illustrative embodiments of the present invention and their descriptions are used to explain the present invention, does not constitute to this hair
Bright improper restriction, such as figure are wherein:
Fig. 1 is the flow diagram of the present invention.
Fig. 2 be the present invention forest three-dimensional scenic in sample ground centre coordinate transition diagram.
Fig. 3 is that forest quantity calculates schematic diagram in the forest three-dimensional scenic sample ground of the present invention.
Fig. 4 calculates example schematic for standing forest densitometer in the forest three-dimensional scenic of the present invention.
Fig. 5 is leaf area index effect diagram in the forest three-dimensional scenic of the present invention.
E, S, W, N of Fig. 2 is respectively the orientation normal vector of all directions four on sample ground;Xyz is world coordinate system;M ' is
Video camera position coordinate;M for sample central point position coordinates;M ' M normal vectors are video camera towards normal vector;d
For video camera point and sample central point distance d=20m;α is angle (0 ° of the video camera towards normal vector and north orientation (N) normal vector
≤α≤180°)。
E, S, W, N of Fig. 4 is respectively the azimuthal point of all directions four on sample ground;A, b, c, e are respectively sample map layer and woods
The intersection point of the wooden bounding box.
Present invention will be further explained below with reference to the attached drawings and examples.
Specific implementation mode
Obviously, those skilled in the art belong to the guarantor of the present invention based on many modifications and variations that spirit of the invention is done
Protect range.
Those skilled in the art of the present technique are appreciated that unless expressly stated, singulative " one " used herein, " one
It is a ", " described " and "the" may also comprise plural form.It is to be further understood that is used in the specification of the present invention arranges
It refers to there are the feature, integer, step, operation, element and/or component, but it is not excluded that presence or addition to take leave " comprising "
Other one or more features, integer, step, operation, element, component and/or their group.It should be understood that when claiming element, group
When part is "connected" to another element, component, it can be directly connected to other elements either component or there may also be in
Between element or component.Wording "and/or" used herein includes any cell of one or more associated list items
With whole combinations.
Those skilled in the art of the present technique are appreciated that unless otherwise defined, all terms used herein (including technology art
Language and scientific terminology) there is meaning identical with the general understanding of the those of ordinary skill in fields.
For ease of the understanding to embodiment, below in conjunction with being further explained explanation, and each embodiment not structure
At limitation of the invention.
Embodiment 1:As shown in Figure 1, Figure 2, shown in Fig. 3, Fig. 4 and Fig. 5, leaf area index is counted in real time in a kind of forest three-dimensional scenic
Calculation method, in forest three-dimensional scenic leaf area index calculating process, it is important to which the then calculating of the density of crop is further built
The density of crop calculates leaf area index with leaf area index relational model.
As shown in Figure 1, leaf area index real-time computing technique in a kind of forest three-dimensional scenic, includes the following steps;
Calculate sample ground forest quantity step in forest three-dimensional scenic;
Calculate density of crop step;
Build the density of crop and leaf area index model step;
Leaf area index in forest three-dimensional scenic calculates step;
It completes.
Wherein:
Calculate forest three-dimensional scenic in sample forest quantity step and calculate density of crop step;Include the following steps;
Standing forest densitometer calculates step in forest three-dimensional scenic,
The density of crop is by traversing using video camera position 20m as sample China firs in 20m × 20m samples ground at center forward
Number is calculated.
System three-dimensional rendering window can real-time display current interface center latitude and longitude coordinates and height, with DEM elevation maps
Coordinate transformation can obtain the x of the central point, y, z coordinate value, using the north orientation of DEM elevation maps as due north, respectively eastwards, west,
South, north respectively extend certain distance, the distance be sample the half of diagonal line length (be sample 20m × 20m, diagonal line length half is
14.14m), thus can respectively obtain the coordinate in four orientation, four azimuthal points be sequentially connected be 20m × 20m square-like
Ground.
A figure layer is built with this square, then each tree bounding box and figure layer crosspoint in the sample ground are inquired by ray
Number (each tree is generated by coordinate data table in system) calculates forest quantity and simultaneously preserves, sample in forest three-dimensional scenic
It is as shown in Figure 3 that interior forest quantity calculates schematic diagram.
The density of crop on the sample ground can be calculated finally by the ratio of forest quantity and sample area in sample ground.
Specific operating procedure is as follows:
Step 1) reads screen center's point coordinates.
Using video camera position, as sample central point, coordinate where reading video camera for the first time are 20m the present invention forward
Topographic map longitude and latitude and height coordinate, and it is P (λ, φ, h) to preserve screen center's point coordinates.
P is screen center's point coordinates,
λ is screen center's point coordinates (longitude),
φ is screen center's point coordinates (latitude),
H is screen center's point coordinates (height).
Step 2), coordinate conversion.As shown in Fig. 2,
After reading video camera position latitude and longitude coordinates and height coordinate, in order to subsequently calculate needs, by the coordinate and
DEM topographic map coordinate transformations obtain M ' (x ', y ', z ') coordinate of the video camera position, and preserve.X ', y ', z ' are respectively
For coordinate value.
Wherein,
A is the conversion coefficient of longitude;
B is the conversion coefficient of latitude;
C is the conversion coefficient of height coordinate;
P is screen center's point coordinates,
λ is screen center's point coordinates (longitude),
φ is screen center's point coordinates (latitude),
H is screen center's point coordinates (height).
Sample central point be located at video camera position 20m forward, sample center point coordinate be M (x, y, z).X, y, z point
It Wei not coordinate value.
1) video camera towards normal vector and east orientation (E) normal vector angle be 0 °~90 °,
2) video camera towards normal vector with west to (W) normal vector angle be 0 °~90 °,
Wherein, α is video camera towards the angle (0 °≤α≤180 °) of normal vector and north orientation (N) normal vector, d=20m,
ztmaxFor sample the maximum forest coordinate of inner height.
Step 3), with establishing sample.
Using the north orientation of DEM elevation maps as due north, respectively extend certain distance to East, West, South, North, northern four direction respectively, it should
Half, that is, 14.14m apart from for sample diagonal line length.
If the coordinate of four azimuthal points in all directions is respectively E (xE, yE, zE), S (xS, yS, zS), W (xW, yW, zW), N
(xN, yN, zN)。
xE、yE、zEThe respectively coordinate in east site, xS、yS、zSThe coordinate in respectively southern site, xW, yW, zWRespectively
The coordinate in west site, xW、yW、zWThe coordinate in respectively northern site.
Wherein, d1=14.14m.
According to the distance value, in conjunction with sample center point coordinate, sample can be respectively obtained the coordinate of four azimuthal points and guarantor
Deposit, four azimuthal points be sequentially connected be 20m × 20m sample.
E, S, W, N of Fig. 2 is respectively the orientation normal vector of all directions four on sample ground;
Xyz is world coordinate system;
M ' is video camera position coordinate;
M for sample central point position coordinates;
M ' M normal vectors are video camera towards normal vector;
D is video camera point with sample central point distance d=20m;
α is angle (0 °≤α≤180 °) of the video camera towards normal vector and north orientation (N) normal vector.
Step 4) calculates forest quantity.As shown in figure 3,
It is a square that the azimuthal point of all directions four on sample ground, which is sequentially connected, virtual with this square structure
Figure layer.
The intersection point number of the figure layer and forest bounding box in forest three-dimensional scenic is inquired by ray, and is preserved.
Step 5), forest inspection of quantity.
After obtaining the intersection point number of sample map layer and forest bounding box, because bounding box has a certain size, individual woodss
Although wood intersects with figure layer, it is nonetheless possible to except sample ground.
Pt (the x of every forest are obtained by intersection pointt, yt, zt) coordinate, the square region that with sample four azimuthal points are constituted
Constraints is established in domain, judges above-mentioned forest whether within the requirement of constraints according to constraints, final determining sample
Interior forest quantity simultaneously preserves.
Constraints is as follows:
Wherein, zESWNminFor sample the minimum z coordinate of four azimuthal points, zESWNmaxFor sample the maximum z of four azimuthal points is sat
Mark.
Step 6) calculates the density of crop.As shown in figure 4, E, S, W, N of Fig. 4 are respectively the side of all directions four on sample ground
Site;
A, b, c, e are respectively the intersection point of sample map layer and forest bounding box.
The density of crop SD on the sample ground can be obtained by the ratio of forest quantity and sample area in sample ground, and is preserved.
The density of crop and leaf area index model step are built, is included the following steps;
Leaf area index calculates step,
It has been embedded in the relational model of the density of crop and leaf area index in system, leaf area index can be calculated.
So far leaf area index has just been calculated, it is aobvious for use in the calling of leaf area index in forest three-dimensional scenic
Show.
Leaf area index in forest three-dimensional scenic calculates step, includes the following steps;
The density of crop in sample ground in known forest three-dimensional scenic, according to the density of crop and leaf area index relational model, certainly
It is dynamic to calculate leaf area index value, and preserve.
LAI=fSD2+gSD+j (6)
Wherein,
LAI is leaf area index,
SD is the density of crop,
SD2For square of the density of crop,
F is density of crop quadratic term,
G is Monomial coefficient,
J is constant.
Embodiment 2:As shown in Figure 1, Figure 2, shown in Fig. 3, Fig. 4 and Fig. 5, leaf area index is counted in real time in a kind of forest three-dimensional scenic
Calculation method, including walk as follows poly-:
Walking in poly- 1, forest three-dimensional scenic forest quantity in sample ground, to calculate step poly-,
Forest number is calculated by taking 20m × 20m samples ground as an example and preserve in forest three-dimensional scenic.
It walks poly- 1.1, set screen center's coordinate as P (750300,3000300,10).
Walk the conversion of poly- 1.2, coordinate.
Video camera position coordinate:
M ' (x ', y ', z ')=(750300,3000300,10)+(- 750000, -3000000,0)
M ' (x ', y ', z ')=(300,200,10),
Wherein,
X ' is camera coordinates (longitude),
Y ' is camera coordinates (latitude),
Z ' is camera coordinates (height),
M ' is camera coordinates,
A is the conversion coefficient of longitude,
B is the conversion coefficient of latitude;
C is the conversion coefficient of height coordinate;
λ is screen center's point coordinates (longitude),
φ is screen center's point coordinates (latitude),
H is screen center's point coordinates (height).
If video camera towards normal vector and east orientation (E) normal vector angle be 0 °~90 °, and the sample ground maximum forest of inner height
Coordinate is ztmax=15, then:
X=x '+dcos α=300+20 × 1/2=310,
Y=y '-dcos α=200-20 × 1/2=190,
Z=ztmax+5=15+5=20,
Therefore M (x, y, z)=(310,190,20),
Wherein,
X ' is camera coordinates (longitude),
Y ' is camera coordinates (latitude),
X for sample center point coordinate (longitudes),
Y for sample center point coordinate (latitudes),
Z for sample center point coordinate (height),
D is distance (20m),
Ztmax for sample maximum forest coordinates of inner height,
α is angle of the video camera towards normal vector and north orientation (N) normal vector.
Walk poly- 1.3, with establishing sample and computer azimuth point coordinates.
If the coordinate of four azimuthal points in all directions is respectively E (xE, yE, zE), S (xS, yS, zS), W (xW, yW, zW), N
(xN, yN, zN).Wherein, d1=14.14m, then:
xE=x+d1=310+14.14=324.14,
yS=y+d1=190+14.14=204.14,
xW=x-d1=310-14.14=285.86,
yN=y-d1=190-14.14=175.86,
Then the coordinate in four orientation in all directions is respectively:
E(xE, yE, zE)=(324.14,190,20),
S(xS, yS, zS)=(310,204.14,20),
W(xW, yW, zW)=(285.86,190,20),
N(xN, yN, zN)=(310,175.86,20).
It walks poly- 1.4, calculate forest quantity,
If with M (310,190,20) coordinate for sample central point when, share 30 plants of forests in 20 × 20m samples ground.
Poly- 2, forest inspection of quantity is walked,
Constraints is established, forest quantity is examined, calculates the density of crop.
The coordinate of 30 plants of forests in sample ground is respectively obtained by intersection point, whether 30 plants of forest coordinates of verification meet constraints
It is required that if 30 plants of forests all meet the requirements, then forest quantity is 30 plants in sample ground after examining.
Density of crop SD=forests quantity/sample area,
SD=30 plants/20m × 20m of the density of crop,
SD=750 plants/ha of the density of crop.
The calculating of poly- 3, leaf area index is walked,
After the density of crop is calculated, according to the embedded constructed density of crop and leaf area index model, it is calculated
Leaf area index, the calling for use in leaf area index in forest three-dimensional scenic are shown.
Leaf area index LAI=aSD2+ bSD+c,
That is LAI=-0.000 0006SD2+0.0015SD+1.423;
Then:
LAI=-0.000 0006 × 7502+0.0015×750+1.423;
:LAI=2.21.
As described above, being explained in detail to the embodiment of the present invention, as long as but essentially without this hair of disengaging
Bright inventive point and effect can have many deformations, this will be readily apparent to persons skilled in the art.Therefore, this
The variation of sample is also integrally incorporated within protection scope of the present invention.
Claims (4)
1. leaf area index real-time computing technique in a kind of forest three-dimensional scenic, it is characterised in that include the following steps;
Calculate sample ground forest quantity step in forest three-dimensional scenic;
Calculate density of crop step;
Build the density of crop and leaf area index model step;
Leaf area index calculating operation step in forest three-dimensional scenic.
2. leaf area index real-time computing technique in a kind of forest three-dimensional scenic according to claim 1, it is characterised in that
Calculate forest three-dimensional scenic in sample forest quantity step and calculate density of crop step;Include the following steps;
Standing forest densitometer calculates step in forest three-dimensional scenic,
The density of crop is by traversing using video camera position 20m as sample China fir numbers in 20m × 20m samples ground at center forward
It is calculated;
System three-dimensional rendering window can real-time display current interface center latitude and longitude coordinates and height, with DEM elevation map coordinates
Conversion can obtain the x of the central point, y, z coordinate value, using the north orientation of DEM elevation maps as due north, respectively to East, West, South, North, north
It is each to extend certain distance, the distance be sample the half of diagonal line length (be sample 20m × 20m, diagonal line length half is
14.14m), thus can respectively obtain the coordinate in four orientation, four azimuthal points be sequentially connected be 20m × 20m square-like
Ground;
With this square one figure layer of structure, then each tree bounding box in the sample ground is inquired by ray and intersects points with figure layer
(each tree is generated by coordinate data table in system) calculates forest quantity and preserves;
The density of crop on the sample ground can be calculated finally by the ratio of forest quantity and sample area in sample ground;
Specific operating procedure is as follows:
Step 1) reads screen center's point coordinates;
Using video camera position, as sample central point, coordinate where reading video camera for the first time are topographic map longitude and latitude to 20m forward
Degree and height coordinate, and it is P (λ, φ, h) to preserve screen center's point coordinates;
P is screen center's point coordinates,
λ is screen center's point coordinates (longitude),
φ is screen center's point coordinates (latitude),
H is screen center's point coordinates (height);
Step 2), coordinate conversion;After reading video camera position latitude and longitude coordinates and height coordinate, needed to subsequently calculate
It wants, by the coordinate and DEM topographic map coordinate transformations, obtains M ' (x ', y ', z ') coordinate of the video camera position, and preserve;
Wherein,
A is the conversion coefficient of longitude;
B is the conversion coefficient of latitude;
C is the conversion coefficient of height coordinate;
P is screen center's point coordinates,
λ is screen center's point coordinates (longitude),
φ is screen center's point coordinates (latitude),
H is screen center's point coordinates (height);
Sample central point be located at video camera position 20m forward, sample center point coordinate be M (x, y, z);
1), video camera towards normal vector and east orientation (E) normal vector angle be 0 °~90 °,
2), video camera towards normal vector with west to (W) normal vector angle be 0 °~90 °,
Wherein, α is video camera towards the angle (0 °≤α≤180 °) of normal vector and north orientation (N) normal vector, d=20m, ztmaxFor sample
The maximum forest coordinate of ground inner height;
Step 3), with establishing sample;
Using the north orientation of DEM elevation maps as due north, respectively extend certain distance to East, West, South, North, northern four direction respectively, the distance
For sample half, that is, 14.14m of diagonal line length;
If the coordinate of four azimuthal points in all directions is respectively E (xE, yE, zE), S (xS, yS, zS), W (xW, yW, zW), N (xN, yN,
zN);
Wherein, d1=14.14m;
According to the distance value, in conjunction with sample center point coordinate, with can respectively obtaining the sample coordinate of four azimuthal points and preservation is somebody's turn to do
Four azimuthal points be sequentially connected be 20m × 20m sample;
E, S, W, N are respectively the orientation normal vector of all directions four on sample ground;
Xyz is world coordinate system;
M ' is video camera position coordinate;
M for sample central point position coordinates;
M ' M normal vectors are video camera towards normal vector;
D is video camera point with sample central point distance d=20m;
α is angle (0 °≤α≤180 °) of the video camera towards normal vector and north orientation (N) normal vector;
Step 4) calculates forest quantity;
It is a square that the azimuthal point of all directions four on sample ground, which is sequentially connected, with this square one virtual graph of structure
Layer;
The intersection point number of the figure layer and forest bounding box in forest three-dimensional scenic is inquired by ray, and is preserved;
Step 5), forest inspection of quantity;
After obtaining the intersection point number of sample map layer and forest bounding box, because bounding box has a certain size, though individual forests
So intersect with figure layer, it is nonetheless possible to except sample ground;
Pt (the x of every forest are obtained by intersection pointt, yt, zt) coordinate, the square area that with sample four azimuthal points are constituted builds
Vertical constraints, judges above-mentioned forest whether within the requirement of constraints according to constraints, final to determine in sample ground
Forest quantity simultaneously preserves;
Constraints is as follows:
Wherein, zESWNminFor sample the minimum z coordinate of four azimuthal points, zESWNmaxFor sample the maximum z coordinate of four azimuthal points;
Step 6) calculates the density of crop;E, S, W, N are respectively the azimuthal point of all directions four on sample ground;
A, b, c, e are respectively the intersection point of sample map layer and forest bounding box;
The density of crop SD on the sample ground can be obtained by the ratio of forest quantity and sample area in sample ground, and is preserved.
3. leaf area index real-time computing technique in a kind of forest three-dimensional scenic according to claim 1, it is characterised in that
The density of crop and leaf area index model step are built, is included the following steps;
Leaf area index calculates step,
It has been embedded in the relational model of the density of crop and leaf area index in system, leaf area index can be calculated;
So far leaf area index has just been calculated, the calling for use in leaf area index in forest three-dimensional scenic is shown.
4. leaf area index real-time computing technique in a kind of forest three-dimensional scenic according to claim 1, it is characterised in that
Steps are as follows for leaf area index calculating operation:
The density of crop in sample ground in known forest three-dimensional scenic, it is automatic to count according to the density of crop and leaf area index relational model
Leaf area index value is calculated, and is preserved;
LAI=fSD2+gSD+j (6)
Wherein,
LAI is leaf area index,
SD is the density of crop,
SD2For square of the density of crop,
F is density of crop quadratic term,
G is Monomial coefficient,
J is constant.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810289091.3A CN108492332B (en) | 2018-04-03 | 2018-04-03 | Real-time calculation method for leaf area index in forest three-dimensional scene |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810289091.3A CN108492332B (en) | 2018-04-03 | 2018-04-03 | Real-time calculation method for leaf area index in forest three-dimensional scene |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108492332A true CN108492332A (en) | 2018-09-04 |
CN108492332B CN108492332B (en) | 2021-05-18 |
Family
ID=63318163
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810289091.3A Expired - Fee Related CN108492332B (en) | 2018-04-03 | 2018-04-03 | Real-time calculation method for leaf area index in forest three-dimensional scene |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108492332B (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111738119A (en) * | 2020-06-12 | 2020-10-02 | 云南电网有限责任公司电力科学研究院 | Forest fire area estimation method and device based on MODIS satellite monitoring |
CN112560661A (en) * | 2020-12-10 | 2021-03-26 | 首都师范大学 | Leaf area index calculation method and device, electronic equipment and readable storage medium |
CN112861336A (en) * | 2021-02-01 | 2021-05-28 | 中国林业科学研究院资源信息研究所 | Virtual simulation method for measuring forest stand section area by angle gauge |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011133451A (en) * | 2009-11-27 | 2011-07-07 | Kyushu Univ | Optical vegetation index sensor |
CN102997871A (en) * | 2012-11-23 | 2013-03-27 | 南京大学 | Method for inverting effective leaf area index by utilizing geometric projection and laser radar |
CN103398957A (en) * | 2013-08-12 | 2013-11-20 | 河海大学 | Hyperspectrum and laser radar-based method for extracting vertical distribution of leaf area |
CN103486928A (en) * | 2013-10-08 | 2014-01-01 | 北京林业大学 | Method for measuring canopy density, density and living vegetation volume under condition of troublesome vegetation |
CN103760565A (en) * | 2014-02-10 | 2014-04-30 | 中国科学院南京地理与湖泊研究所 | Regional scale forest canopy height remote sensing retrieval method |
CN103983230A (en) * | 2014-05-29 | 2014-08-13 | 福州大学 | Verification method for indirect measurement of ground leaf area index |
WO2014122800A1 (en) * | 2013-02-05 | 2014-08-14 | 国立大学法人千葉大学 | Fisheye image data creation program and lai computation program |
CN104462741A (en) * | 2014-09-03 | 2015-03-25 | 中国科学院遥感与数字地球研究所 | City green radiation benefit amount calculation model fused with crown three-dimensional structure |
CN104457626A (en) * | 2014-12-08 | 2015-03-25 | 中国科学院合肥物质科学研究院 | Plant leaf area index measurement method based on laser radar point cloud technology |
CN104778693A (en) * | 2015-04-08 | 2015-07-15 | 云挺 | Leaf area index calculation method based on projection algorithm and active contour model |
CN105654543A (en) * | 2014-09-25 | 2016-06-08 | 薛联凤 | Laser point cloud data-oriented broad-leaved tree real leaf modeling and deforming method |
CN107590540A (en) * | 2017-09-18 | 2018-01-16 | 中国林业科学研究院资源信息研究所 | A kind of forest hat width evaluation method dependent on neighboring trees feature |
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 |
-
2018
- 2018-04-03 CN CN201810289091.3A patent/CN108492332B/en not_active Expired - Fee Related
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011133451A (en) * | 2009-11-27 | 2011-07-07 | Kyushu Univ | Optical vegetation index sensor |
CN102997871A (en) * | 2012-11-23 | 2013-03-27 | 南京大学 | Method for inverting effective leaf area index by utilizing geometric projection and laser radar |
WO2014122800A1 (en) * | 2013-02-05 | 2014-08-14 | 国立大学法人千葉大学 | Fisheye image data creation program and lai computation program |
CN103398957A (en) * | 2013-08-12 | 2013-11-20 | 河海大学 | Hyperspectrum and laser radar-based method for extracting vertical distribution of leaf area |
CN103486928A (en) * | 2013-10-08 | 2014-01-01 | 北京林业大学 | Method for measuring canopy density, density and living vegetation volume under condition of troublesome vegetation |
CN103760565A (en) * | 2014-02-10 | 2014-04-30 | 中国科学院南京地理与湖泊研究所 | Regional scale forest canopy height remote sensing retrieval method |
CN103983230A (en) * | 2014-05-29 | 2014-08-13 | 福州大学 | Verification method for indirect measurement of ground leaf area index |
CN104462741A (en) * | 2014-09-03 | 2015-03-25 | 中国科学院遥感与数字地球研究所 | City green radiation benefit amount calculation model fused with crown three-dimensional structure |
CN105654543A (en) * | 2014-09-25 | 2016-06-08 | 薛联凤 | Laser point cloud data-oriented broad-leaved tree real leaf modeling and deforming method |
CN104457626A (en) * | 2014-12-08 | 2015-03-25 | 中国科学院合肥物质科学研究院 | Plant leaf area index measurement method based on laser radar point cloud technology |
CN104778693A (en) * | 2015-04-08 | 2015-07-15 | 云挺 | Leaf area index calculation method based on projection algorithm and active contour model |
CN107590540A (en) * | 2017-09-18 | 2018-01-16 | 中国林业科学研究院资源信息研究所 | A kind of forest hat width evaluation method dependent on neighboring trees feature |
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 |
Non-Patent Citations (6)
Title |
---|
SETH A. EX 等: "Stand Density Index Estimates Leaf Area Index in Uneven-Aged Ponderosa Pine Stands", 《WESTERN JOURNAL OF APPLIED FORESTRY》 * |
ZAIYANG MA 等: "3D Visual Simulation of Chinese Fir Based on the Influence of Different Stand Spatial Structures", 《2017 2ND INTERNATIONAL CONFERENCE ON IMAGE, VISION AND COMPUTING》 * |
孙华 等: "偏最小二乘回归在Hyperion影像叶面积指数反演中的应用", 《中国农学通报》 * |
林岳峰: "光学与雷达遥感数据联合反演森林生物量方法研究", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
王雄宾 等: "华北落叶松人工林叶面积指数与立地指数、密度、林龄关系模型研究", 《中国农学通报》 * |
黄瑞荣 等: "森林场景中太阳与阴影实时动态可视化模拟", 《南京林业大学学报(自然科学版)》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111738119A (en) * | 2020-06-12 | 2020-10-02 | 云南电网有限责任公司电力科学研究院 | Forest fire area estimation method and device based on MODIS satellite monitoring |
CN111738119B (en) * | 2020-06-12 | 2024-02-09 | 云南电网有限责任公司电力科学研究院 | Forest mountain fire area estimation method and device based on MODIS satellite monitoring |
CN112560661A (en) * | 2020-12-10 | 2021-03-26 | 首都师范大学 | Leaf area index calculation method and device, electronic equipment and readable storage medium |
CN112861336A (en) * | 2021-02-01 | 2021-05-28 | 中国林业科学研究院资源信息研究所 | Virtual simulation method for measuring forest stand section area by angle gauge |
CN112861336B (en) * | 2021-02-01 | 2023-07-25 | 中国林业科学研究院资源信息研究所 | Virtual simulation method for determining forest breaking area by angle gauge |
Also Published As
Publication number | Publication date |
---|---|
CN108492332B (en) | 2021-05-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Schrotter et al. | The digital twin of the city of Zurich for urban planning | |
Sargent et al. | Anthropogenic and biogenic CO2 fluxes in the Boston urban region | |
Zheng et al. | Computational-geometry-based retrieval of effective leaf area index using terrestrial laser scanning | |
CN108492332A (en) | Leaf area index real-time computing technique in a kind of forest three-dimensional scenic | |
CN112256815B (en) | Three-dimensional planning auxiliary examination and approval system and method | |
Zhang et al. | The Application of GIS 3D Modeling and Analysis Technology in Real Estate Mass Appraisal–Taking landscape and sunlight factors as the example | |
Lynch et al. | Perceptual values in the cultural landscape: A computer model for assessing and mapping perceived mystery in rural environments | |
Joly et al. | A quantitative approach to the visual evaluation of landscape | |
Zhang et al. | Digital protection and utilization of architectural heritage using knowledge visualization | |
Dissegna et al. | Modeling mean radiant temperature distribution in urban landscapes using DART | |
Gallardo‐Cruz et al. | Relating species richness to the structure of continuous landscapes: alternative methodological approaches | |
Maan et al. | Tree species biomass and carbon stock measurement using ground based-LiDAR | |
Beiranvand et al. | Evaluation of changes in per capita green space through remote sensing data | |
Zwolinski et al. | Computing and monitoring potential of public spaces by shading analysis using 3d Lidar data and advanced image analysis | |
Guo et al. | 3D Digital protection and representation of burial ruins based on LiDAR and UAV survey | |
CN111291725A (en) | Method for processing homeland utilization data | |
Teng et al. | Extraction and Analysis of Spatial Feature Data of Traditional Villages Based on the Unmanned Aerial Vehicle (UAV) Image | |
Adão et al. | Bringing together UAS-based land surveying and procedural modelling of buildings to set up enhanced VR environments for cultural heritage | |
Ma et al. | Low‐Altitude Photogrammetry and Remote Sensing in UAV for Improving Mapping Accuracy | |
Calvagna et al. | Innovative Multidisciplinary Methodology for the Analysis of Traditional Marginal Architecture | |
Ni et al. | A method for the registration of multiview range images acquired in forest areas using a terrestrial laser scanner | |
KR101089256B1 (en) | Geographic Information analysis method using of programing Invisible Depth Analysis | |
Yang et al. | Capturing spatial patterns of rural landscapes with point cloud | |
Banfi | Drone meets historic building information modelling (HBIM): Unmanned aerial vehicle (UAV) photogrammetry for multi-resolution semantic models | |
Parrinello et al. | Integration and modelling of 3D data as strategy for structural diagnosis in endangered sites. The study case of Church of the Annunciation in Pokcha (Russia) |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210518 |