CN105608735B - Crop root space-division method and system - Google Patents
Crop root space-division method and system Download PDFInfo
- Publication number
- CN105608735B CN105608735B CN201510989947.4A CN201510989947A CN105608735B CN 105608735 B CN105608735 B CN 105608735B CN 201510989947 A CN201510989947 A CN 201510989947A CN 105608735 B CN105608735 B CN 105608735B
- Authority
- CN
- China
- Prior art keywords
- voxel
- range
- cylinder
- root
- broken line
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computer Graphics (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Image Generation (AREA)
Abstract
The present invention discloses a kind of crop root space-division method and system, can describe the cylindrical distribution feature that root system surrounds plant strain growth.Method includes：S1, crop root to be divided is modeled, obtains the root model being made of ordered three-dimensional broken line, wherein, it models the origin of three-dimensional system of coordinate used and comes into being position for root system, Z axis was that root system is come into being position and straight line perpendicular to the ground；The range that S2, specified root system divide defines the scale parameter of cylinder and each voxel, the division of cylinder and scale parameter by Root growing space progress voxel is defined according to the range, wherein, the range defines straight line of the axial line of cylinder where the Z axis；S3, it the range is defined at different levels in the range of cylinder is allocated in each voxel, the Root Parameters in each voxel are carried out with statistics calculating, and quantization profile of the crop root in each voxel is obtained according to result of calculation.
Description
Technical field
The present invention relates to crop root construction techniques fields, and in particular to a kind of crop root space-division method and is
System.
Background technology
Crop root configuration is the important component of crop pattern structure, but at present to the research of its morphosis far not
And aerial part, the main reason is that plant growth obtains its Develop Data very big with measurement analysis difficulty in underground.But root
System is the important sources that crop absorbs moisture and nutrient, and morphological research is most important for the water and fertilizer management of agricultural.
At present in terms of the division of crop root space, there is scholar to divide root system to space using the method for voxel,
Root system space is divided using the method that cube voxel divides, but this method is using cube voxel as division unit, it is difficult
To portray the cylinder distribution characteristics that crop root surrounds plant strain growth.
Invention content
The object of the present invention is to provide a kind of crop root space-division method and systems, can describe root system and surround
The cylindrical distribution feature of plant strain growth.
For this purpose, on the one hand, the present invention proposes a kind of crop root space-division method, including：
S1, crop root to be divided is modeled, obtains the root model being made of ordered three-dimensional broken line,
In, it models the origin of three-dimensional system of coordinate used and comes into being position for root system, Z axis was that root system is come into being position and perpendicular to the ground
Straight line；
The range that S2, specified root system divide defines the scale parameter of cylinder and each voxel, and circle is defined according to the range
Root growing space is carried out the division of voxel with scale parameter by cylinder, wherein, the axial line that the range defines cylinder is
Straight line where the Z axis；
S3, it the range is defined at different levels in the range of cylinder is allocated in each voxel, in each voxel
Root Parameters carry out statistics calculating, and obtain quantization profile of the crop root in each voxel according to result of calculation.
On the other hand, the present invention proposes a kind of crop root space partitioning system, including：
Modeling unit for being modeled to crop root to be divided, obtains the root being made of ordered three-dimensional broken line
Be model, wherein, model the origin of three-dimensional system of coordinate used and come into being position for root system, Z axis be root system come into being position and with
The vertical straight line in ground；
Division unit, for the range that root system divides to be specified to define the scale parameter of cylinder and each voxel, according to described
Range defines the division that Root growing space is carried out voxel with scale parameter by cylinder, wherein, the range defines cylinder
Axial line be the straight line where the Z axis；
Computing unit is allocated for the range to be defined at different levels in the range of cylinder in each voxel, right
Root Parameters in each voxel carry out statistics calculating, and obtain quantization of the crop root in each voxel according to result of calculation
Distribution.
Crop root space-division method and system described in the embodiment of the present invention, based on actual measurement crop root 3-dimensional digital
Change data, the root system in the range of cylinder is defined to range and carries out spatial voxel division, and calculate the root system in each voxel cell
Parameter obtains quantization profile of the crop root in each voxel according to Root Parameters, it is thus possible to describe root system and be given birth to around plant
Long cylindrical distribution feature, so as to instruct the watering during arable farming and fertilizing management.
Description of the drawings
Fig. 1 is the flow diagram of one embodiment of crop root space-division method of the present invention；
Fig. 2 is that the schematic diagram after voxel division is carried out to Maize Roots space；
Fig. 3 is distribution schematic diagram of the corn root section in a voxel；
Fig. 4 is distribution schematic diagram of the corn root section in another voxel；
Fig. 5 is the structure diagram of one embodiment of crop root space partitioning system of the present invention.
Specific embodiment
Purpose, technical scheme and advantage to make the embodiment of the present invention are clearer, below in conjunction with the embodiment of the present invention
In attached drawing, the technical solution in the embodiment of the present invention is explicitly described, it is clear that described embodiment be the present invention
Part of the embodiment, instead of all the embodiments.Based on the embodiments of the present invention, those of ordinary skill in the art are not having
All other embodiments obtained under the premise of creative work are made, shall fall within the protection scope of the present invention.
As shown in Figure 1, the present embodiment discloses a kind of crop root space-division method, including：
S1, crop root to be divided is modeled, obtains the root model being made of ordered three-dimensional broken line,
In, it models the origin of three-dimensional system of coordinate used and comes into being position for root system, Z axis was that root system is come into being position and perpendicular to the ground
Straight line；
The range that S2, specified root system divide defines the scale parameter of cylinder and each voxel, and circle is defined according to the range
Root growing space is carried out the division of voxel with scale parameter by cylinder, wherein, the axial line that the range defines cylinder is
Straight line where the Z axis；
S3, at different levels (including root sections and lateral root) defining the range in the range of cylinder are divided in each voxel
Match, the Root Parameters in each voxel are carried out with statistics calculating, and the crop root is obtained in each voxel according to result of calculation
Quantization profile.
Crop root space-division method described in the embodiment of the present invention, based on actual measurement crop root three-dimensional digital number
According to, the root system in the range of cylinder is defined to range and carries out spatial voxel division, and calculate the Root Parameters in each voxel cell,
Quantization profile of the crop root in each voxel is obtained according to Root Parameters, it is thus possible to describe the circle that root system surrounds plant strain growth
Cylindricality distribution characteristics, so as to instruct the watering during arable farming and fertilizing management.
Optionally, in another embodiment of crop root space-division method of the present invention, what the specified root system divided
Range defines the scale parameter of cylinder and each voxel, including：
The range is specified to define starting altitude value, cylinder height and the cylinder radius of cylinder；
It determines the range defining the division number that cylinder is divided into voxel in height, radius and angle direction respectively
Parameter.
Crop root spatial distribution division methods of the present invention can include five parts：
(1) three dimensional representation is carried out to crop root skeleton structure；
Crop root is multilayer, multilevel hierarchy, and secondary root growth is linear structure on upper level root, by crop
One structural unit of each single conduct of root system, and each structural unit is represented with ordered three-dimensional broken line.
(2) range that root system is specified to divide defines the scale parameter of cylinder and each voxel, according to range and voxel scale
Root growing space is carried out the division of voxel by parameter；
Specified Root growing space range defines 3 parameters of cylinder：Cylinder starting altitude value H0, cylinder height
H, cylinder radius R, wherein, it is come into being axial line of the position straight line perpendicular to the ground as cylinder using crossing root system.
After determining range cylinder, according to practical study demand, determine cylinder respectively in height, radius and angle side
To the division number parameter for being divided into voxel：Nh、NrAnd Na, then the division step-length in three directions be respectively
Cylinder is divided into N=Nh×Nr×NaA voxel, each voxel are segment, and highly identical but shared volume is not
Together, each voxel has unique coordinated indexing in cylindrical space.
(3) root system defined range outside cylinder range is cut off；
It can be solely a part for Root growing space due to calculating target zone, the root system outside target zone is not
It is added to root system space and divides calculating, need rejected in advance.The outer root system excision of target zone is it is possible that include (1) single starting
Part outside target zone (2) single end section outside target zone and the three kinds of situations outside target zone of (3) middle section.
It is detected by traversing current single each broken line and updates single composition and single cut off to each.Excision be calculate it is single with
The intersection point on boundary is divided, intersection point is added to as endpoint in the single broken line, and is divided the point outside range and be then deleted.
(4) range is defined at different levels in the range of cylinder in each voxel to be allocated；
(5) statistics calculating is carried out to the Root Parameters in each voxel, and obtains the quantization parameter in each voxel.
Optionally, it is described to define the range in another embodiment of crop root space-division method of the present invention
At different levels in the range of cylinder are allocated in each voxel, including：
S31, for the single A={ P of each0,P1,…,Pm-1,Pm, calculate the single A={ P0,P1,…,Pm-1,PmIn
Each endpoint PiCorresponding circular cylindrical coordinate (hi,ri,θi), wherein, Pi-1PiTo form i-th single broken line, i ∈ (1,
2 ..., m), hiThe coordinate on cylinder height direction, r are defined for the rangeiCylinder radius direction is defined for the range
On coordinate, θiFor the coordinate on angle direction；
S33, the corresponding N for determining single each broken linerefineAlong ent, and according to the NrefineAlong ent updates
This is single, calculates the NrefineThe corresponding circular cylindrical coordinate of Along ent；
S35, for single each broken line PjPj+1, according to broken line PjPj+1Corresponding voxel index distanceBy broken line PjPj+1It is allocated in each voxel.
Optionally, it is described according to the broken line in another embodiment of crop root space-division method of the present invention
PjPj+1Corresponding voxel index distanceBy broken line PjPj+1It is allocated in each voxel, including：
Judge the voxel index distanceWhether it is 0；
If 0, then by broken line PjPj+1It is assigned to the voxel VjIn or
If not 0, then judge the voxel index distanceWhether 0 is more than；
If more than 0, then broken line P is judgedjPj+1Whether intersect with the axial line, if intersecting, it is determined that go out intersection point Pt,
By line segment PjPtIt is assigned to the voxel VjIn, by line segment PtPj+1It is assigned to the voxel Vj+1In or
If non-intersect, the voxel index distance is judgedWhether it is 1；
If 1, then judge described in | Hje-Hjs|、|Rje-Rjs| andIn whether one be 1, other two be 0, if
It is described | Hje-Hjs|、|Rje-Rjs| andIn one be 1, other two be 0, it is determined that go out broken line PjPj+1With institute
State voxel VjAnd Vj+1The intersection point P of intersectionm, by line segment PjPmIt is assigned to the voxel VjIn, by line segment PmPj+1It is assigned to described
Voxel Vj+1In or
If the voxel index distanceIt is not 1, then judges the voxel index distanceWhether it is 2；
If 2, then judge described in | Hje-Hjs|、|Rje-Rjs| andIn whether one be 0, other two be 1, if
It is described | Hje-Hjs|、|Rje-Rjs| andIn one be 0, other two be 1, it is determined that go out broken line PjPj+1With institute
State voxel VjIntersection point Pm1With broken line PjPj+1With the voxel Vj+1Intersection point Pm2, by line segment PjPm1It is assigned to the voxel
VjIn, by line segment Pm2Pj+1It is assigned to the voxel Vj+1In, by line segment Pm1Pm2Be assigned in the voxel where it or
If not 2, then judge the voxel index distanceWhether it is 3；
If 3, then judge described in | Hje-Hjs|、|Rje-Rjs| andWhether it is 1, if described | Hje-Hjs|、|
Rje-Rjs| andIt is 1, according to step S32 to S35 to broken line PjPj+1Carry out voxel distribution.
For being divided below by the space of maize root system, illustrate the specific implementation method of the present invention, specifically include following step
Suddenly：
(1) three dimensional representation of maize root system
Maize root system is in multistage linear structure, withRepresent one single (including root section or lateral root), wherein n expression layers
Grade is root section during n=0, and when n=1 represents first order root hair；I represents single serial number.Single growth relationship ∈ symbols at different levels
It represents, such asIt representsTo be grown in j-th of level-one root hair(the level-one root hair is grown in i-th root sectionOn) on
Kth two level root hair.
Each single skeleton is represented with broken line, is represented again by orderly discrete point set, such asIts
In, current single skeleton is by the continuous orderly broken line P of m items0P1, P1P2..., Pm-1PmComposition.
(2) Maize Roots spatial voxel divides
It being influenced by soil nutrrient and water, the growth scope of maize root system in the soil is different, therefore, corn root
It is that space-division method only considers root system space distribution situation in range of interest.Due to most maize root systems in the soil with
Corn stem extends downwards growth for axle center, and root system spatial distribution with away from root pitch cone radius from showing certain variation tendency, because
This, the quantization profile for realizing maize root system space using cylindrical space division methods calculates.
Specified Root growing space range defines 3 parameters of cylinder：Cylinder starting altitude value H0, cylinder height
H, cylinder radius R.Wherein, position straight line perpendicular to the ground is come into being as the axial line of cylinder, H using root system0Represent earth's surface more
Highly, if need to consider the root system aerial part such as aerial root, can H be adjusted according to the height of aerial root0；H is distance H0Downward
Depth；R defines the root system part being less than apart from axial line in the range of R.
After determining range cylinder, according to practical study demand, determine cylinder respectively in height, radius and angle side
To the division number parameter for being divided into voxel：Nh、NrAnd Na, then the division step-length in three directions be respectively
Cylinder is divided into N=Nh×Nr×NaA voxel, each voxel are segment, and highly identical but shared volume is not
Together, each voxel has unique coordinated indexing in cylindrical space, such as (Hi,Rj,Ak)。
(3) the outer root system excision of target zone
It can be solely a part for Root growing space due to calculating target zone, the root system outside target zone is not
It is added to root system space and divides calculating, need rejected in advance.The outer root system of target zone, which is cut with, may include (1) single starting
Part outside target zone (2) single end section outside target zone and the three kinds of situations outside target zone of (3) middle section.
It detects by traversing currently single each broken line and updates single composition, withFor, it is assumed that P0,
P1,…,PiIn target zone, and Pi+1..., Pm-1, PmOutside target zone, each point P is traversedi, when there is PiIn the range of and
Pi+1Outside range during situation, P is calculatediPi+1With the intersection point on cylinder boundary, it is denoted as Ps, then after cutting
(4) distribution of the maize root system in each voxel calculates
In order to which root system is assigned in each voxel, with singleFor, algorithm description is as follows：
(41) it is directed toEach point P on skeletoni, calculate its corresponding circular cylindrical coordinate (hi,ri,θi)。
(42) the larger line segment of single middle span is segmented.In order to reduce computation complexity so that every line segment is three sides
The span of (short transverse, radial direction and angle direction) is less than the step-length of the direction upwards, to ensure every line segment in the party
It at most falls in two adjacent voxels upwards.With singleOn line segment Pi-1PiFor, utilize the head calculated in step (41)
The segmentation fractal methods that the circular cylindrical coordinate of tail point calculates the line segment are：
Wherein,For lower bracket function；Further, since angle is a circle for closing, therefore neighbour is used when calculating angle
Processing, i.e.,：
According to the N calculatedrefineValue, in line segment PiPi+1In be equidistantly inserted into Nrefine- 1 point, and updateIt calculates simultaneously
The circular cylindrical coordinate of newly-increased point.
According toValue, distribution of the line segment in voxel include a variety of situations：
If 1)Then PjWith in Pj+1In same voxel, without calculating line segment and voxel intersection point, by current line segment point
It is fitted on voxel VjIn.
If 2)And PjPj+1Intersect with the axial line of cylinder, calculate intersection point Pt, by line segment PjPtIt is assigned to voxel Vj
In, line segment PtPj+1It is assigned to voxel Vj+1In.
If 4)And PjPj+1It is non-intersect with the axial line of cylinder, P at this timejAnd Pj+1Place voxel and some intermediate
Element is adjacent, when one of following three kinds of situations of appearance：
By taking situation one as an example, P is calculated respectivelyjPj+1With VjAnd Vj+1Intersection point, be denoted as Pm1And Pm2, by line segment PjPm1It is assigned to
Voxel VjIn, line segment Pm2Pj+1It is assigned to voxel Vj+1In, while by line segment Pm1Pm2It is stored in its corresponding voxel.Similarly may be used
The situation of solution situation two and situation three.
If 5)And PjPj+1Intersect with the axial line of cylinder, whenWhen, such case is seldom
Occur, once occurring, divided line segment again using the method for step (42), and the method processing being converted into (43).
Optionally, in another embodiment of crop root space-division method of the present invention, the root in each voxel
It is that parameter carries out statistics calculating, and obtain quantization profile of the crop root in each voxel according to result of calculation, including：
To each voxel, the quantity of at different levels, length, volume and surface area in the voxel are calculated, and according to described at different levels
Quantity, length, volume and the surface area of root obtain quantization profile of the crop root in each voxel.
It is rare to root system spatial distribution in addition, crop root configuration research at present is carried out mainly for topological structure
Quantitative description method, through the embodiment of the present invention in spatial distribution quantitative method crop root can be described in the soil
Distribution situation.
Optionally, it in another embodiment of crop root space-division method of the present invention, before the S3, also wraps
It includes：
The root system range defined outside cylinder range is cut off.
Optionally, it is described to define the range in another embodiment of crop root space-division method of the present invention
Root system excision outside cylinder range, including：
For each single of the crop root, single each broken line is traversed, by two endpoints described
Range defines the broken line excision outside cylinder range, and one endpoint is defined the broken line outside cylinder range in the range
The intersection point on the boundary of cylinder is defined with the range and the line of the endpoint outside cylinder range should be defined in the range
The broken line excision of composition.
As shown in figure 5, the present embodiment discloses a kind of crop root space partitioning system, including：
Modeling unit 1 for being modeled to crop root to be divided, obtains the root being made of ordered three-dimensional broken line
Be model, wherein, model the origin of three-dimensional system of coordinate used and come into being position for root system, Z axis be root system come into being position and with
The vertical straight line in ground；
Division unit 2, for the range that root system divides to be specified to define the scale parameter of cylinder and each voxel, according to described
Range defines the division that Root growing space is carried out voxel with scale parameter by cylinder, wherein, the range defines cylinder
Axial line be the straight line where the Z axis；
Computing unit 3 is allocated for the range to be defined at different levels in the range of cylinder in each voxel, right
Root Parameters in each voxel carry out statistics calculating, and obtain quantization of the crop root in each voxel according to result of calculation
Distribution.
Crop root space partitioning system described in the embodiment of the present invention, based on actual measurement crop root three-dimensional digital number
According to, the root system in the range of cylinder is defined to range and carries out spatial voxel division, and calculate the Root Parameters in each voxel cell,
Quantization profile of the crop root in each voxel is obtained according to Root Parameters, it is thus possible to describe the circle that root system surrounds plant strain growth
Cylindricality distribution characteristics, so as to instruct the watering during arable farming and fertilizing management.
Optionally, in another embodiment of crop root space partitioning system of the present invention, the division unit, for referring to
The fixed range defines starting altitude value, cylinder height and the cylinder radius of cylinder, determines the range defining circle
Cylinder is divided into the division number parameter of voxel in height, radius and angle direction respectively.
Optionally, in another embodiment of crop root space partitioning system of the present invention, the computing unit, including：
Circular cylindrical coordinate computation subunit, for for the single A={ P of each0,P1,…,Pm-1,Pm, calculate the single A
={ P0,P1,…,Pm-1,PmIn each endpoint PiCorresponding circular cylindrical coordinate (hi,ri,θi), wherein, Pi-1PiIt is single to form this
I-th broken line, i ∈ (1,2 ..., m), hiThe coordinate on cylinder height direction, r are defined for the rangeiFor the range
Define the coordinate on cylinder radius direction, θiFor the coordinate on angle direction；
Along ent determination subelement, for determining the corresponding N of single each broken linerefineAlong ent, and according to institute
State NrefineAlong ent updates that this is single, calculates the NrefineThe corresponding circular cylindrical coordinate of Along ent；
Subelement is distributed, for for single each broken line PjPj+1, according to broken line PjPj+1Corresponding body
Element index distanceBy broken line PjPj+1It is allocated in each voxel.
Although being described in conjunction with the accompanying embodiments of the present invention, those skilled in the art can not depart from this hair
Various modifications and variations are made in the case of bright spirit and scope, such modifications and variations are each fallen within by appended claims
Within limited range.
Claims (8)
1. a kind of crop root space-division method, which is characterized in that including：
S1, crop root to be divided is modeled, obtains the root model being made of ordered three-dimensional broken line, wherein, it builds
The origin of three-dimensional system of coordinate used in mould is come into being position for root system, and Z axis was that root system is come into being position and straight line perpendicular to the ground；
The range that S2, specified root system divide defines the scale parameter of cylinder and each voxel, and cylinder is defined according to the range
Root growing space is carried out to the division of voxel with scale parameter, wherein, the axial line that the range defines cylinder is the Z
Straight line where axis；
S3, it the range is defined at different levels in the range of cylinder is allocated in each voxel, to the root system in each voxel
Parameter carries out statistics calculating, and obtains quantization profile of the crop root in each voxel according to result of calculation；
It is described the range is defined at different levels in the range of cylinder in each voxel to be allocated, including：
S31, for the single A={ P of each0,P1,…,Pm-1,Pm, calculate this it is single in each endpoint PiCorresponding circular cylindrical coordinate
(hi,ri,θi), wherein, Pi-1PiTo form i-th single broken line, i ∈ (1,2 ..., m), hiCircle is defined for the range
Coordinate on height of column direction, riThe coordinate on cylinder radius direction, θ are defined for the rangeiFor the seat on angle direction
Mark；
S32, for each broken line Pi-1Pi, calculate broken line Pi-1PiSegmentation number Nrefine, calculation formula isWherein,For lower bracket function,H be the cylinder height, NhFor drawing for the short transverse
Divide number parameter,R be the cylinder radius, NrFor the division number parameter of the radial direction,NaFor
The division number parameter of the angle direction；
S33, the corresponding N for determining single each broken linerefineAlong ent, and according to the NrefineAlong ent updates the list
Root calculates the NrefineThe corresponding circular cylindrical coordinate of Along ent；
S34, for single each broken line PjPj+1, according to point PjCorresponding circular cylindrical coordinate searches point PjThe voxel at place
Vj, according to point Pj+1Corresponding circular cylindrical coordinate searches point Pj+1The voxel V at placej+1, and calculate voxel index distanceIt calculates public
Formula isWherein,
HjeFor voxel VjThe coordinated indexing in cylinder height direction, H are defined in the rangejsFor voxel Vj+1Circle is defined in the range
The coordinated indexing in height of column direction, RjeFor voxel VjThe coordinated indexing in cylinder radius direction, R are defined in the rangejsFor
Voxel Vj+1The coordinated indexing in cylinder radius direction, A are defined in the rangejeFor voxel VjIn the coordinated indexing of angle direction,
AjsFor voxel Vj+1In the coordinated indexing of angle direction；
S35, for single each broken line PjPj+1, according to broken line PjPj+1Corresponding voxel index distanceIt will
Broken line PjPj+1It is allocated in each voxel.
2. crop root space-division method according to claim 1, which is characterized in that the model that the specified root system divides
The scale parameter for defining cylinder and each voxel is enclosed, including：
The range is specified to define starting altitude value, cylinder height and the cylinder radius of cylinder；
It determines the range defining cylinder and is divided into the division number of voxel in height, radius and angle direction respectively and joins
Number.
3. crop root space-division method according to claim 1, which is characterized in that described according to broken line PjPj+1
Corresponding voxel index distanceBy broken line PjPj+1It is allocated in each voxel, including：
Judge the voxel index distanceWhether it is 0；
If 0, then by broken line PjPj+1It is assigned to the voxel VjIn or
If not 0, then judge the voxel index distance dVWhether j is more than 0；
If more than 0, then broken line P is judgedjPj+1Whether intersect with the axial line, if intersecting, it is determined that go out intersection point Pt, by line
Section PjPtIt is assigned to the voxel VjIn, by line segment PtPj+1It is assigned to the voxel Vj+1In or
If non-intersect, the voxel index distance is judgedWhether it is 1；
If 1, then judge described in | Hje-Hjs|、|Rje-Rjs| andIn whether one be 1, other two be 0, if described |
Hje-Hjs|、|Rje-Rjs| andIn one be 1, other two be 0, it is determined that go out broken line PjPj+1With the voxel
VjAnd Vj+1The intersection point P of intersectionm, by line segment PjPmIt is assigned to the voxel VjIn, by line segment PmPj+1It is assigned to the voxel Vj+1
In or
If the voxel index distanceIt is not 1, then judges the voxel index distanceWhether it is 2；
If 2, then judge described in | Hje-Hjs|、|Rje-Rjs| andIn whether one be 0, other two be 1, if described |
Hje-Hjs|、|Rje-Rjs| andIn one be 0, other two be 1, it is determined that go out broken line PjPj+1With the voxel
VjIntersection point Pm1With broken line PjPj+1With the voxel Vj+1Intersection point Pm2, by line segment PjPm1It is assigned to the voxel VjIn,
By line segment Pm2Pj+1It is assigned to the voxel Vj+1In, by line segment Pm1Pm2Be assigned in the voxel where it or
If not 2, then judge the voxel index distanceWhether it is 3；
If 3, then judge described in | Hje-Hjs|、|Rje-Rjs| andWhether it is 1, if described | Hje-Hjs|、|Rje-Rjs|
WithIt is 1, according to step S32 to S35 to broken line PjPj+1Carry out voxel distribution.
4. crop root space-division method according to claim 1, which is characterized in that the root system in each voxel
Parameter carries out statistics calculating, and obtains quantization profile of the crop root in each voxel according to result of calculation, including：
To each voxel, the quantity of at different levels, length, volume and surface area in the voxel are calculated, and according to described at different levels
Quantity, length, volume and surface area obtain quantization profile of the crop root in each voxel.
5. crop root space-division method according to claim 1, which is characterized in that before the S3, further include：
The root system range defined outside cylinder range is cut off.
6. crop root space-division method according to claim 5, which is characterized in that described that the range is defined into circle
Root system excision outside cylinder range, including：
For each single of the crop root, single each broken line is traversed, by two endpoints in the range
The broken line excision outside cylinder range is defined, an endpoint is defined in the range to broken line and the institute outside cylinder range
State range define cylinder boundary intersection point and should the range define the endpoint outside cylinder range line composition
Broken line excision.
7. a kind of crop root space partitioning system, which is characterized in that including：
Modeling unit for being modeled to crop root to be divided, obtains the root system mould being made of ordered three-dimensional broken line
Type, wherein, it models the origin of three-dimensional system of coordinate used and comes into being position for root system, Z axis was that root system is come into being position and and ground
Vertical straight line；
Division unit, for the range that root system divides to be specified to define the scale parameter of cylinder and each voxel, according to the range
The division that Root growing space is carried out voxel with scale parameter by cylinder is defined, wherein, the range defines the axis of cylinder
Heart line is the straight line where the Z axis；
Computing unit is allocated for the range to be defined at different levels in the range of cylinder in each voxel, to each body
Root Parameters in element carry out statistics calculating, and obtain quantization of the crop root in each voxel point according to result of calculation
Cloth；
The computing unit, including：
Circular cylindrical coordinate computation subunit, for for the single A={ P of each0,P1,…,Pm-1,Pm, calculate the single A={ P0,
P1,…,Pm-1,PmIn each endpoint PiCorresponding circular cylindrical coordinate (hi,ri,θi), wherein, Pi-1PiTo form this single i-th
Broken line, i ∈ (1,2 ..., m), hiThe coordinate on cylinder height direction, r are defined for the rangeiCircle is defined for the range
Coordinate on cylinder radius direction, θiFor the coordinate on angle direction；
Divide number computation subunit, for for each broken line Pi-1Pi, calculate broken line Pi-1PiSegmentation number Nrefine,
Calculation formula isWherein,For lower bracket function,H be the cylinder height, NhFor drawing for the short transverse
Divide number parameter,R be the cylinder radius, NrFor the division number parameter of the radial direction,NaFor
The division number parameter of the angle direction；
Along ent determination subelement, for determining the corresponding N of single each broken linerefineAlong ent, and according to described
NrefineAlong ent updates that this is single, calculates the NrefineThe corresponding circular cylindrical coordinate of Along ent；
Voxel index is apart from computation subunit, for for single each broken line PjPj+1, according to point PjCorresponding circle
Cylindrical coordinates searches point PjThe voxel V at placej, according to point Pj+1Corresponding circular cylindrical coordinate searches point Pj+1The voxel V at placej+1, and count
Calculate voxel index distanceCalculation formula isWherein,HjeFor voxel VjThe seat in cylinder height direction is defined in the range
Mark index, HjsFor voxel Vj+1The coordinated indexing in cylinder height direction, R are defined in the rangejeFor voxel VjIn the range
Define the coordinated indexing in cylinder radius direction, RjsFor voxel Vj+1The coordinate rope in cylinder radius direction is defined in the range
Draw, AjeFor voxel VjIn the coordinated indexing of angle direction, AjsFor voxel Vj+1In the coordinated indexing of angle direction；
Subelement is distributed, for for single each broken line PjPj+1, according to broken line PjPj+1Corresponding voxel rope
Draw distanceBy broken line PjPj+1It is allocated in each voxel.
8. crop root space partitioning system according to claim 7, which is characterized in that the division unit, for referring to
The fixed range defines starting altitude value, cylinder height and the cylinder radius of cylinder, determines the range defining circle
Cylinder is divided into the division number parameter of voxel in height, radius and angle direction respectively.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510989947.4A CN105608735B (en) | 2015-12-24 | 2015-12-24 | Crop root space-division method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510989947.4A CN105608735B (en) | 2015-12-24 | 2015-12-24 | Crop root space-division method and system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105608735A CN105608735A (en) | 2016-05-25 |
| CN105608735B true CN105608735B (en) | 2018-07-06 |
Family
ID=55988646
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510989947.4A Active CN105608735B (en) | 2015-12-24 | 2015-12-24 | Crop root space-division method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105608735B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109858180B (en) * | 2019-02-26 | 2020-11-10 | 中国农业大学 | Regional crop water consumption space pattern optimization design method |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2008263850A (en) * | 2007-04-19 | 2008-11-06 | Tottori Univ | Method and apparatus for assuming root system distribution |
| CN101650836A (en) * | 2009-09-10 | 2010-02-17 | 北京农业信息技术研究中心 | Self-adaptive gridding method and self-adaptive gridding system of geometric curved surfaces of three-dimensional plant organs |
| CN101706968A (en) * | 2009-12-10 | 2010-05-12 | 江苏大学 | Image-based method for reconstructing three-dimensional models of fruit tree limbs |
| CN101887596A (en) * | 2010-06-01 | 2010-11-17 | 中国科学院自动化研究所 | Three-dimensional model reconstruction method of tree point cloud data based on partition and automatic growth |
| CN103728418A (en) * | 2013-12-31 | 2014-04-16 | 马玮 | Analyzing method for researching parameters of architecture and spatial distribution of crop roots |
| CN104091012A (en) * | 2014-07-02 | 2014-10-08 | 河南工程学院 | Virtual root system construction method based on group behavior environment sensitivity |
-
2015
- 2015-12-24 CN CN201510989947.4A patent/CN105608735B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2008263850A (en) * | 2007-04-19 | 2008-11-06 | Tottori Univ | Method and apparatus for assuming root system distribution |
| CN101650836A (en) * | 2009-09-10 | 2010-02-17 | 北京农业信息技术研究中心 | Self-adaptive gridding method and self-adaptive gridding system of geometric curved surfaces of three-dimensional plant organs |
| CN101706968A (en) * | 2009-12-10 | 2010-05-12 | 江苏大学 | Image-based method for reconstructing three-dimensional models of fruit tree limbs |
| CN101887596A (en) * | 2010-06-01 | 2010-11-17 | 中国科学院自动化研究所 | Three-dimensional model reconstruction method of tree point cloud data based on partition and automatic growth |
| CN103728418A (en) * | 2013-12-31 | 2014-04-16 | 马玮 | Analyzing method for researching parameters of architecture and spatial distribution of crop roots |
| CN104091012A (en) * | 2014-07-02 | 2014-10-08 | 河南工程学院 | Virtual root system construction method based on group behavior environment sensitivity |
Non-Patent Citations (3)
| Title |
|---|
| Reconciling root plasticity and architectural ground rules in tree root growth models with voxel automata;Rachmat Mulia 等;《Plant & Soil》;20100812;第337卷(第1期);第77-92页 * |
| 作物根系构型三维探测与重建方法研究进展;温维亮 等;《中国农业科学》;20150130;第48卷(第3期);第436-448页 * |
| 均质壤土下玉米根系三维空间分布的模拟与参数分析;张吴平;《土壤学报》;20070930;第44卷(第5期);第949-955页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105608735A (en) | 2016-05-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103942838A (en) | Point cloud data based single tree three-dimensional modeling and morphological parameter extracting method | |
| CN108012608B (en) | A kind of prepartion of land method based on GNSS | |
| CN101887596B (en) | Three-dimensional model reconstruction method of tree point cloud data based on partition and automatic growth | |
| CN106327577B (en) | Dimensional topography optimal method based on local curvature's entropy and quad-tree structure | |
| CN104486435B (en) | Low energy consumption ecological environmental monitoring node deployment method based on sensor network | |
| CN101650836B (en) | Self-adaptive gridding method and self-adaptive gridding system of geometric curved surfaces of three-dimensional plant organs | |
| CN105654552B (en) | A kind of quick Delaunay network construction methods towards Arbitrary distribution large-scale point cloud data | |
| CN106643529A (en) | Rapid measuring method for growing height of agricultural crops in mountainous areas based on unmanned aerial vehicle photographed image | |
| CN106340061B (en) | A kind of mountain area point cloud filtering method | |
| CN107358579B (en) | Game war fog-lost realization method | |
| CN102930509A (en) | Intelligent filtering method for airborne laser point cloud data | |
| CN102915560B (en) | Threshold-irrelative point cloud filtering method and device for airborne laser radar | |
| CN105354335B (en) | For the time-varying tree graph layout method of time-varying hierarchical data and application | |
| CN104318611A (en) | Three-dimensional tree re-establishing method based on laser scanning point cloud | |
| Ersson et al. | Simulated productivity of one-and two-armed tree planting machines | |
| CN109002418A (en) | Tree breast-height diameter automatic calculating method based on Voxels growing and ground laser point cloud | |
| CN109146990B (en) | Building outline calculation method | |
| CN107103603A (en) | A kind of Objects extraction method of inclination measurement scene | |
| CN105608735B (en) | Crop root space-division method and system | |
| CN105785919B (en) | Point five axis of cloud is without global interference generating tool axis vector feasible zone computational methods | |
| CN103432754B (en) | A kind of corn maze and method for designing thereof | |
| CN110335329A (en) | A method of power circuit cross-section diagram is obtained based on DEM | |
| CN110472317A (en) | A method of the sprinkler tip spatial configuration optimal for Campus green space | |
| CN106447510A (en) | Terrace field information extracting method and terrace field information extracting system | |
| CN106373192A (en) | Non-topology-consistency three-dimensional grid block tracking algorithm |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |