CN105608735A - Crop root space division method and system - Google Patents

Abstract

The present invention discloses a crop root space division method and system, and cylindrical distribution characteristic of the root growth around a plant can be described. The method comprises the step of (S1) carrying out modeling of a crop root to be divided, and obtaining the root model formed by ordered three-dimensional fold line sections, wherein the origin of the three-dimensional coordinate system used for the modeling is a root newborn position, a Z axis is a straight line which passes the root newborn position and is perpendicular to the ground; (S2) specifying the range defining cylinder of root division and the scale parameter of each voxel, and carrying out the division of voxel on a root growth space according to the range defining cylinder of the root division and the scale parameter, wherein the axis of the range defining cylinder is the straight line where the Z axis is located, (S3) distributing each level of root in the range of the range defining cylinder in each voxel, carrying out statistical calculation of the root parameter in each voxel, and obtaining the quantitative distribution of the crop root in each voxel according to a calculation result.

Description

Crop root space-division method and system

Technical field

The present invention relates to crop root configuration technical field, be specifically related to a kind of crop root space-division method and beSystem.

Background technology

Crop root configuration is the important component part of crop pattern structure, but at present to the research of its morphosis far notAnd aerial part, its main cause be plant growth in underground, obtain with Measurement and analysis difficulty very big to its Develop Data. But rootSystem is the important sources of Crop moisture and nutrient, and its morphological research is most important for agriculture water and fertilizer management.

Aspect crop root spatial division, there is scholar to adopt the method for voxel to divide space to root system at present, itsThe method that adopts cube voxel to divide is divided root system space, but the method is taking cube voxel as division unit, difficultyTo portray the cylinder distribution characteristics of crop root around plant strain growth.

Summary of the invention

The object of the invention is to, a kind of crop root space-division method and system are provided, can describe root system aroundThe 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, comprising:

S1, crop root to be divided is carried out to modeling, obtain the root system model that formed by orderly three-dimensional polyline section, itsIn, the initial point of modeling three-dimensional system of coordinate used is the nascent position of root system, Z axis was the nascent position and perpendicular to the ground of root systemStraight line;

S2, the scope definition cylinder of specifying root system division and the scale parameter of each voxel, according to described scope definition circleCylinder and scale parameter carry out Root growing space the division of voxel, and wherein, the cylindrical axial line of described scope definition isThe straight line at described Z axis place;

S3, at different levels within the scope of described scope definition cylinder are distributed in each voxel, in each voxelRoot Parameters carries out statistical computation, and obtains the quantization profile of described crop root in each voxel according to result of calculation.

On the other hand, the present invention proposes a kind of crop root spatial division system, comprising:

Modeling unit, carries out modeling for the crop root to be divided, obtains the root being made up of orderly three-dimensional polyline sectionBe model, wherein, the initial point of modeling three-dimensional system of coordinate used is the nascent position of root system, Z axis be the nascent position of root system and withThe straight line that ground is vertical;

Division unit, is used to specify the scope definition cylinder of root system division and the scale parameter of each voxel, according to describedScope definition cylinder and scale parameter carry out Root growing space the division of voxel, wherein, and described scope definition cylinderAxial line be the straight line at described Z axis place;

Computing unit, at different levels within the scope of described scope definition cylinder are distributed in each voxel, rightRoot Parameters in each voxel carries out statistical computation, and obtains the quantification of described crop root in each voxel according to result of calculationDistribute.

Crop root space-division method and system described in the embodiment of the present invention, based on actual measurement crop root 3-dimensional digitalChange data, the root system within the scope of scope definition cylinder is carried out to space voxel division, and calculate the root system in each voxel cellParameter, obtains the quantization profile of crop root in each voxel according to Root Parameters, thereby it is raw around plant to describe root systemLong cylindrical distribution feature, thus watering and fertilizing management in arable farming process can be instructed.

Brief description of the drawings

Fig. 1 is the schematic flow sheet of crop root space-division method one embodiment of the present invention;

Fig. 2 carries out the schematic diagram after voxel division to Maize Roots space;

Fig. 3 is the distribution schematic diagram of corn root section in a voxel;

Fig. 4 is the distribution schematic diagram of corn root section in another voxel;

Fig. 5 is the structural representation of crop root spatial division system of the present invention one embodiment.

Detailed description of the invention

For making object, technical scheme and the advantage of the embodiment of the present invention clearer, below in conjunction with the embodiment of the present inventionIn accompanying drawing, the technical scheme in the embodiment of the present invention is clearly described, obviously, described embodiment is the present inventionPart embodiment, instead of whole embodiment. Based on the embodiment in the present invention, those of ordinary skill in the art are not havingMake the every other embodiment obtaining under creative work prerequisite, all belong to the scope of protection of the invention.

As shown in Figure 1, the present embodiment discloses a kind of crop root space-division method, comprising:

S1, crop root to be divided is carried out to modeling, obtain the root system model that formed by orderly three-dimensional polyline section, itsIn, the initial point of modeling three-dimensional system of coordinate used is the nascent position of root system, Z axis was the nascent position and perpendicular to the ground of root systemStraight line;

S2, the scope definition cylinder of specifying root system division and the scale parameter of each voxel, according to described scope definition circleCylinder and scale parameter carry out Root growing space the division of voxel, and wherein, the cylindrical axial line of described scope definition isThe straight line at described Z axis place;

S3, at different levels within the scope of described scope definition cylinder (comprising root section and lateral root) are divided in each voxelJoin, the Root Parameters in each voxel is carried out to statistical computation, and obtain described crop root in each voxel according to result of calculationQuantization profile.

Crop root space-division method described in the embodiment of the present invention, based on actual measurement crop root three-dimensional digital numberAccording to, the root system within the scope of scope definition cylinder is carried out to space voxel division, and calculates the Root Parameters in each voxel cell,Obtain the quantization profile of crop root in each voxel according to Root Parameters, thereby can describe the circle of root system around plant strain growthCylindricality distribution characteristics, thus watering and fertilizing management in arable farming process can be instructed.

Alternatively, in another embodiment of crop root space-division method of the present invention, described appointment root system is dividedThe scale parameter of scope definition cylinder and each voxel, comprising:

Specify the cylindrical starting altitude value of described scope definition, cylinder height and cylinder radius;

Determine the division number that described scope definition cylinder is divided into respectively to voxel at height, radius and angle directionParameter.

Crop root spatial distribution division methods of the present invention can comprise five parts:

(1) crop root skeleton structure is carried out to three dimensional representation;

Crop root is multilayer, multilevel hierarchy, and secondary root growth, on upper level root, and is linear structure, by cropRoot system each single as a construction unit, and show each construction unit with orderly three-dimensional polyline segment table.

(2) specify the scope definition cylinder of root system division and the scale parameter of each voxel, according to scope and voxel yardstickParameter is carried out Root growing space the division of voxel;

Specify cylindrical 3 parameters of Root growing space scope definition: cylinder starting altitude value H₀, cylinder heightH, cylinder radius R, wherein, to cross the nascent position of root system straight line perpendicular to the ground as cylindrical axial line.

Determine after range circle cylinder, according to practical study demand, determine cylinder respectively in height, radius and angle sideDivision number parameter to being divided into voxel: N_h、N_rAnd N_a, the division step-length of three directions is respectively

$\begin{array}{ccc}h=\frac{H}{N_h}& r=\frac{R}{N_r}& \mathrm{\θ}=\frac{2\mathrm{\π}}{N_a}\end{array}.$

Cylinder is divided into N=N_h×N_r×N_aIndividual voxel, each voxel is segment, highly identical but shared volume is notWith, each voxel has unique coordinated indexing in cylindrical space.

(3) by extraneous scope definition cylinder root system excision;

May be only a part for Root growing space owing to calculating target zone, therefore, the root system outside target zone beJoin root system spatial division and calculate, need reject in advance. It is single initial that the outer root system excision of target zone likely comprises (1)Part target zone outside (2) single end portion outside target zone with (3) mid portion three kinds of situations outside target zone.Detect and upgrade single composition to each single excision by traveling through current single each broken line. Excision be calculate single withThe intersection point of dividing border, intersection point joins in this single broken line as end points, divides extraneous point deleted.

(4) at different levels within the scope of scope definition cylinder are distributed in each voxel;

(5) Root Parameters in each voxel is carried out to statistical computation, and obtain the quantization parameter in each voxel.

Alternatively, in another embodiment of crop root space-division method of the present invention, described by described scope definitionAt different levels within the scope of cylinder are distributed, comprising in each voxel:

S31, for each single A={P₀,P₁,…,P_m-1,P_m, calculate this single A={P₀,P₁,…,P_m-1,P_mInEach end points P_iCorresponding circular cylindrical coordinate (h_i,r_i,θ_i), wherein, P_i-1P_iFor forming this i article of single broken line, i ∈ (1,2,…,m)，h_iFor the coordinate in described scope definition cylinder height direction, r_iFor described scope definition cylinder radial directionOn coordinate, θ_iFor the coordinate on angle direction;

S32, for each broken line P_i-1P_i, calculate this broken line P_i-1P_iSegmentation Number N_refine, computing formula isWherein, $\stackrel{\&OverBar;}{|{\mathrm{\&theta;}}_{i-1}-{\mathrm{\&theta;}}_i|}=\left\{\begin{array}{ccc}|{\mathrm{\&theta;}}_{i-1}-{\mathrm{\&theta;}}_i|& if& |{\mathrm{\&theta;}}_{i-1}-{\mathrm{\&theta;}}_i|<\mathrm{\&pi;}\\ 2\mathrm{\&pi;}|{\mathrm{\&theta;}}_{i-1}-{\mathrm{\&theta;}}_i|& if& |{\mathrm{\&theta;}}_{i-1}-{\mathrm{\&theta;}}_i|>\mathrm{\&pi;}\end{array}\right.,h=\frac{H}{N_h},$ H is described cylinder height, N_hFor the division number parameter of described short transverse,R is described cylinder radius, N_rForThe division number parameter of described radial direction,N_aFor the division number parameter of described angle direction;

S33, determine the corresponding N of this each single broken line_refineAlong ent, and according to described N_refineAlong ent upgradesThis is single, calculates described N_refineThe circular cylindrical coordinate that Along ent is corresponding;

S34, for this each single broken line P_jP_j+1, according to a P_jCorresponding circular cylindrical coordinate is searched a P_jThe voxel at placeV_j, according to a P_j+1Corresponding circular cylindrical coordinate is searched a P_j+1The voxel V at place_j+1, and calculate voxel index distanceComputing formula is $d_V^j=|H_{je}-H_{js}|+|R_{je}-R_{js}|+\stackrel{\&OverBar;}{|A_{je}-A_{js}|},$ Wherein,For lower bracket function, $\stackrel{\&OverBar;}{|A_{je}-A_{js}|}=\{\begin{array}{ccc}|A_{je}-A_{js}|& if& |A_{je}-A_{js}|\&le;N_a/2\\ N_a-|A_{je}-A_{js}|& if& |A_{je}-A_{js}|>N_a/2\end{array},$ H_jeFor voxel V_jAt the coordinated indexing of described scope definition cylinder height direction, H_jsFor voxel V_j+1At described scope definition circleThe coordinated indexing of height of column direction, R_jeFor voxel V_jAt the coordinated indexing of described scope definition cylinder radial direction, R_jsForVoxel V_j+1At the coordinated indexing of described scope definition cylinder radial direction, A_jeFor voxel V_jIn the coordinated indexing of angle direction,A_jsFor voxel V_j+1In the coordinated indexing of angle direction;

S35, for this each single broken line P_jP_j+1, according to this broken line P_jP_j+1Corresponding voxel index distanceBy this broken line P_jP_j+1In each voxel, distribute.

Alternatively, in another embodiment of crop root space-division method of the present invention, described according to this broken lineP_jP_j+1Corresponding voxel index distanceBy this broken line P_jP_j+1In each voxel, distribute, comprising:

Judge described voxel index distanceWhether be 0;

If 0, by this broken line P_jP_j+1Be assigned to described voxel V_jIn, or

If not 0, judge described voxel index distanceWhether be greater than 0;

If be greater than 0, judge this broken line P_jP_j+1Whether crossing with described axial line, if intersect, determine intersection point P_t，By line segment P_jP_tBe assigned to described voxel V_jIn, by line segment P_tP_j+1Be assigned to described voxel V_j+1In, or

If non-intersect, judge described voxel index distanceWhether be 1;

If 1, described in judgement | H_je-H_js|、|R_je-R_js| andIn whether one be 1, other two is 0, ifDescribed | H_je-H_js|、|R_je-R_js| andIn one be 1, other two is 0, determines this broken line P_jP_j+1With instituteState voxel V_jAnd V_j+1The intersection point P of phase cross surface_m, by line segment P_jP_mBe assigned to described voxel V_jIn, by line segment P_mP_j+1Described in being assigned toVoxel V_j+1In, or

If described voxel index distanceBe not 1, judge described voxel index distanceWhether be 2;

If 2, described in judgement | H_je-H_js|、|R_je-R_js| andIn whether one be 0, other two is 1, ifDescribed | H_je-H_js|、|R_je-R_js| andIn one be 0, other two is 1, determines this broken line P_jP_j+1With instituteState voxel V_jIntersection point P_m1With this broken line P_jP_j+1With described voxel V_j+1Intersection point P_m2, by line segment P_jP_m1Be assigned to described voxelV_jIn, by line segment P_m2P_j+1Be assigned to described voxel V_j+1In, by line segment P_m1P_m2Be assigned in the voxel at its place, or

If not 2, judge described voxel index distanceWhether be 3;

If 3, described in judgement | H_je-H_js|、|R_je-R_js| andWhether be 1, if described in | H_je-H_js|、|R_je-R_js| andBe 1, according to step S32 to S35 to this broken line P_jP_j+1Carry out voxel distribution.

Taking the spatial division of maize root system as example, specific embodiment of the invention method is described below, specifically comprises following stepRapid:

(1) three dimensional representation of maize root system

Maize root system is multistage linear structure, withRepresent one single (comprising root section or lateral root), wherein n presentation layerLevel, is root section when n=0, represents the first order root hair when n=1; I represents single sequence number. Single growth relationship ∈ symbols at different levelsRepresent, asRepresentFor being grown in j one-level root hair(this one-level root hair is grown in i article of root sectionOn) onK article of secondary root hair.

Each single skeleton represents with broken line, and it represents by orderly discrete point set again, asItsIn, current single skeleton is by the continuously orderly broken line P of m bar₀P₁，P₁P₂，…，P_m-1P_mComposition.

(2) Maize Roots space voxel is divided

Be subject to the impact of soil nutrrient and water, the growth scope of maize root system in soil is different, therefore, and corn rootBe that space-division method is only considered the root system spatial distribution situation in range of interest. Due to most maize root systems in soil withCorn stem is that growth is expanded in axle center downwards, and root system spatial distribution is with apart from root pitch cone radius from presenting certain variation tendency, because ofThis, the quantization profile that adopts cylindrical space division methods to realize maize root system space is calculated.

Specify cylindrical 3 parameters of Root growing space scope definition: cylinder starting altitude value H₀, cylinder heightH, cylinder radius R. Wherein, taking the nascent position of root system straight line perpendicular to the ground as cylindrical axial line, H₀Represent earth's surface moreHighly, consider as the root system aerial part of aerial root if need, can be according to the height adjustment H of aerial root₀; H is distance H₀DownwardThe degree of depth; R defines apart from axial line and is less than the root system part within the scope of R.

Determine after range circle cylinder, according to practical study demand, determine cylinder respectively in height, radius and angle sideDivision number parameter to being divided into voxel: N_h、N_rAnd N_a, the division step-length of three directions is respectively

$\begin{array}{ccc}h=\frac{H}{N_h}& r=\frac{R}{N_r}& \mathrm{\&theta;}=\frac{2\mathrm{\&pi;}}{N_a}\end{array}.$

Cylinder is divided into N=N_h×N_r×N_aIndividual voxel, each voxel is segment, highly identical but shared volume is notWith, each voxel has unique coordinated indexing in cylindrical space, as (H_i,R_j,A_k)。

(3) the outer root system excision of target zone

May be only a part for Root growing space owing to calculating target zone, therefore, the root system outside target zone beJoin root system spatial division and calculate, need reject in advance. The outer root system of target zone is cut with that may to comprise (1) single initialPart target zone outside (2) single end portion outside target zone with (3) mid portion three kinds of situations outside target zone.Detect and upgrade single composition by traveling through current single each broken line, withFor example, suppose P₀,P₁,…,P_iIn target zone, and P_i+1，…，P_m-1，P_mOutside target zone, traversal each point P_i, when occurring P_iIn scope andP_i+1Outside scope, when situation, calculate P_iP_i+1Intersection point with cylinder border, is designated as P_s, after cutting

(4) Distribution Calculation of maize root system in each voxel

For root system being assigned in each voxel, with singleFor example, arthmetic statement is as follows:

(41) forEach some P on skeleton_i, calculate its corresponding circular cylindrical coordinate (h_i,r_i,θ_i)。

(42) segment the larger line segment of single middle span. In order to reduce computation complexity, make every line segment three sidesUpwards the span of (short transverse, radial direction and angle direction) is less than the step-length of this direction, to ensure that every line segment is the partyUpwards drop at most in two adjacent voxels. With singleOn line segment P_i-1P_iFor example, utilize the head calculating in step (41)Point counting method of cutting apart that the circular cylindrical coordinate of tail point is calculated this line segment is:

Wherein,For lower bracket function; In addition, because angle is a circle for sealing, therefore adopt neighbour in the time calculating angleProcess, that is:

$\stackrel{\&OverBar;}{|{\mathrm{\&theta;}}_{i-1}-{\mathrm{\&theta;}}_i|}=\left\{\begin{array}{ccc}|{\mathrm{\&theta;}}_{i-1}-{\mathrm{\&theta;}}_i|& if& |{\mathrm{\&theta;}}_{i-1}-{\mathrm{\&theta;}}_i|<\mathrm{\&pi;}\\ 2\mathrm{\&pi;}|{\mathrm{\&theta;}}_{i-1}-{\mathrm{\&theta;}}_i|& if& |{\mathrm{\&theta;}}_{i-1}-{\mathrm{\&theta;}}_i|>\mathrm{\&pi;}\end{array}\right.,$

According to calculated N_refineValue, at line segment P_iP_i+1In the equidistant N that inserts_refine-1 point, and upgradeCalculate simultaneouslyThe circular cylindrical coordinate of newly-increased point.

(43) calculate the intersection point of single and associated voxels and preserve the single data in each voxel. With singleOn line segmentP_jP_j+1For example, search a P according to the circular cylindrical coordinate of head and the tail end points_jAnd P_j+1The voxel at place, is designated as respectively V_jAnd V_j+1, two individualitiesThe cylinder coordinated indexing of element is designated as respectively (H_js,R_js,A_js) and (H_je,R_je,A_je), calculate voxel index distanceWherein $d_H^j=|H_{je}-H_{js}|,d_R^j=|R_{je}-R_{js}|,d_A^j=\stackrel{\&OverBar;}{|A_{je}-A_{js}|},$

$\stackrel{\&OverBar;}{|A_i-A_j|}=\{\begin{array}{ccc}|A_i-A_j|& if& |A_i-A_j|\&le;N_a/2\\ N_a-|A_i-A_j|& if& |A_i-A_j|>N_a/2\end{array}.$

According toValue, the distribution of line segment in voxel comprises multiple situation:

1) ifP_jWith at P_j+1In same voxel, without calculating line segment and voxel intersection point, current line segment is dividedBe fitted on voxel V_jIn.

2) ifAnd P_jP_j+1Crossing with the axial line of cylinder, calculate intersection point P_t, by line segment P_jP_tBe assigned to voxel V_jIn, line segment P_tP_j+1Be assigned to voxel V_j+1In.

3) ifAnd P_jP_j+1Non-intersect with the axial line of cylinder, now P_jAnd P_j+1In two adjacent voxels, if, in two adjacent voxels of short transverse, calculate P_jP_j+1With voxel V_jAnd V_j+1Phase cross surfaceIntersection point P_m, by line segment P_jP_mBe assigned to voxel V_jIn, line segment P_mP_j+1Be assigned to voxel V_j+1In. In like manner can solve radial direction phaseAdjacent $(d_R^j=1,$ $\left.\begin{array}{cc}{d}_H^j=0,& d_A^j=0,\end{array}\right)$ Or angle direction is adjacent $\left(\begin{array}{ccc}{d}_A^j=1,& d_H^j=0,& d_R^j=0\end{array}\right)$ Situation.

4) ifAnd P_jP_j+1Non-intersect with the axial line of cylinder, now P_jAnd P_j+1Place voxel and certain intermediateElement is adjacent, one of three kinds of situations below occurring:

Situation one: $\begin{array}{ccc}{d}_H^j=1,& d_R^j=1,& d_A^j=0\end{array};$

Situation two: $\begin{array}{ccc}{d}_H^j=1,& d_R^j=0,& d_A^j=1\end{array};$

Situation three: $\begin{array}{ccc}{d}_H^j=0,& d_R^j=1,& d_A^j=1\end{array}.$

Taking situation one as example, calculate respectively P_jP_j+1With V_jAnd V_j+1Intersection point, be designated as P_m1And P_m2, by line segment P_jP_m1Be assigned toVoxel V_jIn, line segment P_m2P_j+1Be assigned to voxel V_j+1In, simultaneously by line segment P_m1P_m2Be stored in its corresponding voxel. In like manner canThe situation of solution situation two and situation three.

5) ifAnd P_jP_j+1Crossing with the axial line of cylinder, whenTime, this situation is littleOccur, once occur, utilize the method for step (42) that line segment is divided again, and be converted into the method processing in (43).

Alternatively, in another embodiment of crop root space-division method of the present invention, described to the root in each voxelBe that parameter is carried out statistical computation, and obtain the quantization profile of described crop root in each voxel according to result of calculation, comprising:

To each voxel, calculate in this voxel quantity, length, volume and the surface area of at different levels, and according to described at different levelsQuantity, length, volume and the surface area of root obtains the quantization profile of described crop root in each voxel.

In addition, at present crop root configuration research is carried out mainly for topological structure, rare to root system spatial distributionQuantitative description method, can describe crop root in soil by the spatial distribution quantitative method in the embodiment of the present inventionDistribution situation.

Alternatively, in another embodiment of crop root space-division method of the present invention, before described S3, also bagDraw together:

By extraneous described scope definition cylinder root system excision.

Alternatively, in another embodiment of crop root space-division method of the present invention, described by described scope definitionThe extraneous root system excision of cylinder, comprising:

Each single for described crop root, travels through this single each broken line, by two end points all describedThe excision of the extraneous broken line of scope definition cylinder, by an end points at the extraneous broken line of described scope definition cylinderWith the intersection point on the cylindrical border of described scope definition, and should be at the line of the extraneous end points of described scope definition cylinderThe broken line excision of composition.

As shown in Figure 5, the present embodiment discloses a kind of crop root spatial division system, comprising:

Modeling unit 1, carries out modeling for the crop root to be divided, obtains the root being made up of orderly three-dimensional polyline sectionBe model, wherein, the initial point of modeling three-dimensional system of coordinate used is the nascent position of root system, Z axis be the nascent position of root system and withThe straight line that ground is vertical;

Division unit 2, is used to specify the scope definition cylinder of root system division and the scale parameter of each voxel, according to describedScope definition cylinder and scale parameter carry out Root growing space the division of voxel, wherein, and described scope definition cylinderAxial line be the straight line at described Z axis place;

Computing unit 3, at different levels within the scope of described scope definition cylinder are distributed in each voxel, rightRoot Parameters in each voxel carries out statistical computation, and obtains the quantification of described crop root in each voxel according to result of calculationDistribute.

Crop root spatial division system described in the embodiment of the present invention, based on actual measurement crop root three-dimensional digital numberAccording to, the root system within the scope of scope definition cylinder is carried out to space voxel division, and calculates the Root Parameters in each voxel cell,Obtain the quantization profile of crop root in each voxel according to Root Parameters, thereby can describe the circle of root system around plant strain growthCylindricality distribution characteristics, thus watering and fertilizing management in arable farming process can be instructed.

Alternatively, in another embodiment of crop root spatial division system of the present invention, described division unit, is used in reference toThe fixed cylindrical starting altitude value of described scope definition, cylinder height and cylinder radius, determine described scope definition circleCylinder is divided into respectively the division number parameter of voxel at height, radius and angle direction.

Alternatively, in another embodiment of crop root spatial division system of the present invention, described computing unit, comprising:

Circular cylindrical coordinate computation subunit, for for each single A={P₀,P₁,…,P_m-1,P_m, calculate this single A＝{P₀,P₁,…,P_m-1,P_mIn each end points P_iCorresponding circular cylindrical coordinate (h_i,r_i,θ_i), wherein, P_i-1P_iSingle for forming thisI article of broken line, i ∈ (1,2 ..., m), h_iFor the coordinate in described scope definition cylinder height direction, r_iFor described scopeDefine the coordinate on cylinder radial direction, θ_iFor the coordinate on angle direction;

Segmentation Number computation subunit, for for each broken line P_i-1P_i, calculate this broken line P_i-1P_iSegmentation NumberN_refine, computing formula isWherein,For under roundFunction, $\stackrel{\&OverBar;}{|{\mathrm{\&theta;}}_{i-1}-{\mathrm{\&theta;}}_i|}=\left\{\begin{array}{ccc}|{\mathrm{\&theta;}}_{i-1}-{\mathrm{\&theta;}}_i|& if& |{\mathrm{\&theta;}}_{i-1}-{\mathrm{\&theta;}}_i|<\mathrm{\&pi;}\\ 2\mathrm{\&pi;}-|{\mathrm{\&theta;}}_{i-1}-{\mathrm{\&theta;}}_i|& if& |{\mathrm{\&theta;}}_{i-1}-{\mathrm{\&theta;}}_i|>\mathrm{\&pi;}\end{array}\right.,h=\frac{H}{N_h},$ H is described cylinder height, N_hFor described height sideTo division number parameter,R is described cylinder radius, N_rFor the division number parameter of described radial direction,N_aFor the division number parameter of described angle direction;

Along ent is determined subelement, for determining the corresponding N of this each single broken line_refineAlong ent, and according to instituteState N_refineAlong ent upgrades that this is single, calculates described N_refineThe circular cylindrical coordinate that Along ent is corresponding;

Voxel index is apart from computation subunit, for for this each single broken line P_jP_j+1, according to a P_jRightThe circular cylindrical coordinate of answering is searched a P_jThe voxel V at place_j, according to a P_j+1Corresponding circular cylindrical coordinate is searched a P_j+1The body at placeElement V_j+1, and calculate voxel index distanceComputing formula is $d_V^j=|H_{je}-H_{js}|+|R_{je}-R_{js}|+\stackrel{\&OverBar;}{|A_{je}-A_{js}|},$ Wherein, $\stackrel{\&OverBar;}{|A_{je}-A_{js}|}=\{\begin{array}{ccc}|A_{je}-A_{js}|& if& |A_{je}-A_{js}|\&le;N_a/2\\ N_a-|A_{je}-A_{js}|& if& |A_{je}-A_{js}|>N_a/2\end{array},$ H_jeFor voxel V_jIn described scope definition cylinder height directionCoordinated indexing, H_jsFor voxel V_j+1At the coordinated indexing of described scope definition cylinder height direction, R_jeFor voxel V_jAt described modelEnclose the coordinated indexing that defines cylinder radial direction, R_jsFor voxel V_j+1At the coordinate of described scope definition cylinder radial directionIndex, A_jeFor voxel V_jAt the coordinated indexing of angle direction, A_jsFor voxel V_j+1In the coordinated indexing of angle direction;

Distribute subelement, for for this each single broken line P_jP_j+1, according to this broken line P_jP_j+1Corresponding bodyElement index distanceBy this broken line P_jP_j+1In each voxel, distribute.

Although described by reference to the accompanying drawings embodiments of the present invention, those skilled in the art can not depart from thisIn the situation of bright spirit and scope, make various modifications and variations, such amendment and modification all fall into by claimsWithin limited range.

Claims (10)

1. a crop root space-division method, is characterized in that, comprising:

S1, crop root to be divided is carried out to modeling, obtain the root system model that formed by orderly three-dimensional polyline section, wherein, buildThe initial point of mould three-dimensional system of coordinate used is the nascent position of root system, and Z axis was the nascent position of root system and straight line perpendicular to the ground;

S2, the scope definition cylinder of specifying root system division and the scale parameter of each voxel, according to described scope definition cylinderWith scale parameter, Root growing space is carried out to the division of voxel, wherein, the cylindrical axial line of described scope definition is described ZThe straight line at axle place;

S3, at different levels within the scope of described scope definition cylinder are distributed in each voxel, to the root system in each voxelParameter is carried out statistical computation, and obtains the quantization profile of described crop root in each voxel according to result of calculation.

2. crop root space-division method according to claim 1, is characterized in that, the model that described appointment root system is dividedEnclose the scale parameter that defines cylinder and each voxel, comprising:

Specify the cylindrical starting altitude value of described scope definition, cylinder height and cylinder radius;

Determine the division number ginseng that described scope definition cylinder is divided into respectively to voxel at height, radius and angle directionNumber.

3. crop root space-division method according to claim 2, is characterized in that, described by described scope definition circleAt different levels within the scope of cylinder are distributed, comprising in each voxel:

S31, for each single A={P₀,P₁,…,P_m-1,P_m, calculate this single A={P₀,P₁,…,P_m-1,P_mIn eachEnd points P_iCorresponding circular cylindrical coordinate (h_i,r_i,θ_i), wherein, P_i-1P_iFor forming this i article of single broken line, i ∈ (1,2 ...,m)，h_iFor the coordinate in described scope definition cylinder height direction, r_iFor the seat on described scope definition cylinder radial directionMark, θ_iFor the coordinate on angle direction;

S32, for each broken line P_i-1P_i, calculate this broken line P_i-1P_iSegmentation Number N_refine, computing formula isWherein,For lower bracket function, $\stackrel{\&OverBar;}{|{\mathrm{\&theta;}}_{i-1}-{\mathrm{\&theta;}}_i|}=\left\{\begin{array}{cc}|{\mathrm{\&theta;}}_{i-1}-{\mathrm{\&theta;}}_i|& if|{\mathrm{\&theta;}}_{i-1}-{\mathrm{\&theta;}}_i|<\mathrm{\&pi;}\\ 2\mathrm{\&pi;}-|{\mathrm{\&theta;}}_{i-1}-{\mathrm{\&theta;}}_i|& if|{\mathrm{\&theta;}}_{i-1}-{\mathrm{\&theta;}}_i|>\mathrm{\&pi;}\end{array}\right.,$ H is described cylinder height, N_hFor drawing of described short transverseDivide number parameter,R is described cylinder radius, N_rFor the division number parameter of described radial direction,N_aForThe division number parameter of described angle direction;

S33, determine the corresponding N of this each single broken line_refineAlong ent, and according to described N_refineAlong ent upgrades this listRoot, calculates described N_refineThe circular cylindrical coordinate that Along ent is corresponding;

S34, for this each single broken line P_jP_j+1, according to a P_jCorresponding circular cylindrical coordinate is searched a P_jThe voxel at placeV_j, according to a P_j+1Corresponding circular cylindrical coordinate is searched a P_j+1The voxel V at place_j+1, and calculate voxel index distanceCalculate publicFormula is $d_V^j=|H_{je}-H_{js}|+|R_{je}-R_{js}|+\stackrel{\&OverBar;}{|A_{je}-A_{js}|},$ Wherein, $\stackrel{\&OverBar;}{|A_{je}-A_{js}|}=\left\{\begin{array}{cc}|A_{je}-A_{js}|& if|A_{je}-A_{js}|\&le;N_a/2\\ N_a-|A_{je}-A_{js}|& if|A_{je}-A_{js}|>N_a/2\end{array}\right.,$ H_jeForVoxel V_jAt the coordinated indexing of described scope definition cylinder height direction, H_jsFor voxel V_j+1At described scope definition cylinderThe coordinated indexing of short transverse, R_jeFor voxel V_jAt the coordinated indexing of described scope definition cylinder radial direction, R_jsFor voxelV_j+1At the coordinated indexing of described scope definition cylinder radial direction, A_jeFor voxel V_jAt the coordinated indexing of angle direction, A_jsForVoxel V_j+1In the coordinated indexing of angle direction;

S35, for this each single broken line P_jP_j+1, according to this broken line P_jP_j+1Corresponding voxel index distanceWillThis broken line P_jP_j+1In each voxel, distribute.

4. crop root space-division method according to claim 3, is characterized in that, described according to this broken line P_jP_j+1Corresponding voxel index distanceBy this broken line P_jP_j+1In each voxel, distribute, comprising:

Judge described voxel index distanceWhether be 0;

If 0, by this broken line P_jP_j+1Be assigned to described voxel V_jIn, or

If not 0, judge described voxel index distanceWhether be greater than 0;

If be greater than 0, judge this broken line P_jP_j+1Whether crossing with described axial line, if intersect, determine intersection point P_t, by lineSection P_jP_tBe assigned to described voxel V_jIn, by line segment P_tP_j+1Be assigned to described voxel V_j+1In, or

If non-intersect, judge described voxel index distanceWhether be 1;

If 1, described in judgement | H_je-H_js|、|R_je-R_js| andIn whether one be 1, other two is 0, if described in|H_je-H_js|、|R_je-R_js| andIn one be 1, other two is 0, determines this broken line P_jP_j+1With described bodyElement V_jAnd V_j+1The intersection point P of phase cross surface_m, by line segment P_jP_mBe assigned to described voxel V_jIn, by line segment P_mP_j+1Be assigned to described voxelV_j+1In, or

If described voxel index distanceBe not 1, judge described voxel index distanceWhether be 2;

If 2, described in judgement | H_je-H_js|、|R_je-R_js| andIn whether one be 0, other two is 1, if described in|H_je-H_js|、|R_je-R_js| andIn one be 0, other two is 1, determines this broken line P_jP_j+1With described bodyElement V_jIntersection point P_m1With this broken line P_jP_j+1With described voxel V_j+1Intersection point P_m2, by line segment P_jP_m1Be assigned to described voxel V_jIn, by line segment P_m2P_j+1Be assigned to described voxel V_j+1In, by line segment P_m1P_m2Be assigned in the voxel at its place, or

If not 2, judge described voxel index distanceWhether be 3;

If 3, described in judgement | H_je-H_js|、|R_je-R_js| andWhether be 1, if described in | H_je-H_js|、|R_je-R_js|WithBe 1, according to step S32 to S35 to this broken line P_jP_j+1Carry out voxel distribution.

5. crop root space-division method according to claim 1, is characterized in that, described to the root system in each voxelParameter is carried out statistical computation, and obtains the quantization profile of described crop root in each voxel according to result of calculation, comprising:

To each voxel, calculate in this voxel quantity, length, volume and the surface area of at different levels, and according to described at different levelsQuantity, length, volume and surface area obtain the quantization profile of described crop root in each voxel.

6. crop root space-division method according to claim 1, is characterized in that, before described S3, also comprises:

By extraneous described scope definition cylinder root system excision.

7. crop root space-division method according to claim 6, is characterized in that, described by described scope definition circleThe extraneous root system excision of cylinder, comprising:

Each single for described crop root, travels through this single each broken line, by two end points all in described scopeDefine the excision of the extraneous broken line of cylinder, by an end points in the extraneous broken line of described scope definition cylinder and instituteState the intersection point on the cylindrical border of scope definition, and should be at the line composition of the extraneous end points of described scope definition cylinderBroken line excision.

8. a crop root spatial division system, is characterized in that, comprising:

Modeling unit, carries out modeling for the crop root to be divided, obtains the root system mould being made up of orderly three-dimensional polyline sectionType, wherein, the initial point of modeling three-dimensional system of coordinate used is the nascent position of root system, Z axis was the nascent position of root system and and groundVertical straight line;

Division unit, is used to specify the scope definition cylinder of root system division and the scale parameter of each voxel, according to described scopeDefine the division that cylinder and scale parameter carry out Root growing space voxel, wherein, the cylindrical axle of described scope definitionHeart line is the straight line at described Z axis place;

Computing unit, at different levels within the scope of described scope definition cylinder are distributed in each voxel, to each bodyRoot Parameters in element carries out statistical computation, and obtains the quantification of described crop root in each voxel according to result of calculation and divideCloth.

9. crop root spatial division system according to claim 8, is characterized in that, described division unit, is used in reference toThe fixed cylindrical starting altitude value of described scope definition, cylinder height and cylinder radius, determine described scope definition circleCylinder is divided into respectively the division number parameter of voxel at height, radius and angle direction.

10. crop root spatial division system according to claim 9, is characterized in that, described computing unit, comprising:

Circular cylindrical coordinate computation subunit, for for each single A={P₀,P₁,…,P_m-1,P_m, calculate this single A={P₀,P₁,…,P_m-1,P_mIn circular cylindrical coordinate (h corresponding to each end points Pi_i,r_i,θ_i), wherein, P_i-1P_iFor forming this single i articleBroken line, i ∈ (1,2 ..., m), h_iFor the coordinate in described scope definition cylinder height direction, r_iFor described scope definition circleCoordinate on cylinder radial direction, θ_iFor the coordinate on angle direction;

Segmentation Number computation subunit, for for each broken line P_i-1P_i, calculate this broken line P_i-1P_iSegmentation Number N_refine，Computing formula isWherein,For lower bracket function, $\stackrel{\&OverBar;}{|{\mathrm{\&theta;}}_{i-1}-{\mathrm{\&theta;}}_i|}=\left\{\begin{array}{cc}|{\mathrm{\&theta;}}_{i-1}-{\mathrm{\&theta;}}_i|& if|{\mathrm{\&theta;}}_{i-1}-{\mathrm{\&theta;}}_i|<\mathrm{\&pi;}\\ 2\mathrm{\&pi;}-|{\mathrm{\&theta;}}_{i-1}-{\mathrm{\&theta;}}_i|& if|{\mathrm{\&theta;}}_{i-1}-{\mathrm{\&theta;}}_i|>\mathrm{\&pi;}\end{array}\right.,$ H is described cylinder height, N_hFor drawing of described short transverseDivide number parameter,R is described cylinder radius, N_rFor the division number parameter of described radial direction,N_aForThe division number parameter of described angle direction;

Along ent is determined subelement, for determining the corresponding N of this each single broken line_refineAlong ent, and according to describedN_refineAlong ent upgrades that this is single, calculates described N_refineThe circular cylindrical coordinate that Along ent is corresponding;

Voxel index is apart from computation subunit, for for this each single broken line P_jP_j+1, according to a P_jCorresponding circular cylindrical coordinate is searchedPoint P_jThe voxel V at place_j, according to a P_j+1Corresponding circular cylindrical coordinate is searched a P_j+1The voxel V at place_j+1, and calculate voxel index distanceMeterCalculation formula is $d_V^j=|H_{je}-H_{js}|+|R_{je}-R_{js}|+\stackrel{\&OverBar;}{|A_{je}-A_{js}|},$ Wherein, $\stackrel{\&OverBar;}{|A_{je}-A_{js}|}=\left\{\begin{array}{cc}|A_{je}-A_{js}|& if|A_{je}-A_{js}|\&le;N_a/2\\ N_a-|A_{je}-A_{js}|& if|A_{je}-A_{js}|>N_a/2\end{array}\right.,$ H_jeFor voxel V_jAt the coordinated indexing of described scope definition cylinder height direction, H_jsFor voxel V_j+1At described scope definition circleThe coordinated indexing of height of column direction, R_jeFor voxel V_jAt the coordinated indexing of described scope definition cylinder radial direction, R_jsForVoxel V_j+1At the coordinated indexing of described scope definition cylinder radial direction, A_jeFor voxel V_jIn the coordinated indexing of angle direction,A_jsFor voxel V_j+1In the coordinated indexing of angle direction;

Distribute subelement, for for this each single broken line P_jP_j+1, according to this broken line P_jP_j+1Corresponding voxel ropeDraw distanceBy this broken line P_jP_j+1In each voxel, distribute.