Summary of the invention
The object of the present invention is to provide a kind of spraying emulation mode based on the simulation of droplet cloud, by simulating actual farmland operation complex environment, set up virtual spraying analogue system, can provide spraying operation multinomial visual operational indicator analysis, and to take the effective deposition of droplet be objective function, optimize every technical parameter in spraying system, thereby provide certain reference for actual plant protection operation process.
In order to achieve the above object, the technical solution used in the present invention is:
The step of the inventive method is as follows:
1) the applying three-dimensional stereoscopic camera obtains plant body three dimensional point cloud, and deletes method and coordinate transform processing by a cloud, simulates actual crop plant;
2), under definite shower nozzle model, shower nozzle coordinate and spray technique Parameter Conditions, describe operator by Lagrangian track and describe spraying droplet Particles Moving track;
3) design a kind of spraying droplet Particles Moving track and plant point cloud intersection algorithm, differentiate droplet and whether drop on the plant body and calculate its drop point site;
4), according to the spraying operation rate request, constantly change the shower nozzle volume coordinate and complete dynamic similation, until export the mist droplet deposition rate and optimize spray height after reaching stop condition.
AB+AC-BC<2AD
Described continuous variation shower nozzle volume coordinate completes dynamic similation, can simulate the Dynamic spraying process, nozzle position is changed according to the spraying operation rate request, after reaching the emulation step number Maxstep that stop condition complete requirement, calculate the sprayed deposit rate of this operation height condition lower nozzle, by repeatedly analog result comparison, realize the optimization of spray height parameter, the computing formula of droplet product rate is as follows:
Wherein C is the mist droplet deposition rate;
for the droplet mark vector, dropping on value on plant is 1, and the value on the ground that falls is 0, N
ijfor all droplet mark vector values are 1, m
ijfor Quality of fogdrop corresponding to droplet identification number.
The beneficial effect that the present invention has is:
The present invention can provide spraying operation visual operational indicator analysis, and to take the effective deposition of droplet be objective function, acquisition optimum spraying operation height, thereby provide certain reference for reducing actual plant protection operation process Pesticides use amount.
Embodiment
Below in conjunction with drawings and Examples, the invention will be further described.
As shown in Figure 1, realization flow of the present invention is as follows:
1) the applying three-dimensional stereoscopic camera obtains plant body three dimensional point cloud, and deletes method and coordinate transform processing by a cloud, simulates actual crop plant;
2), under definite shower nozzle model, shower nozzle coordinate and spray technique Parameter Conditions, describe operator by Lagrangian track and describe spraying droplet Particles Moving track;
3) design a kind of spraying droplet Particles Moving track and plant point cloud intersection algorithm, differentiate droplet and whether drop on the plant body and calculate its drop point site;
4), according to the spraying operation rate request, constantly change the shower nozzle volume coordinate and complete dynamic similation, until export the mist droplet deposition rate and optimize spray height after reaching stop condition.
S110: use the 3 D stereo camera to get and obtain the plant cloud data, utilize the distance samples method to complete a cloud and delete, sampled distance is r, and the some cloud after deleting is carried out to coordinate transform, and coordinate transformation algorithm is described below:
Choosing plant model root is some N, and the stem that chooses upwards growth is some M, vector
for the normal vector of plane Π, N is a bit in the Π plane, and the formula of establishing Π is Ax+By+Cz+D=0,
?
?
Bring the N point coordinate into and obtain the XOY plane formula, select two uneven vectors on XOY plane, and carry out orthogonalization unit by Schmidt orthogonalization, the Schmidt orthogonalization procedure is as follows:
If α
1, α
2..., α
slinear independence, construct β
1, β
2..., β
smake its pairwise orthogonal, and β
ionly α
1, α
2..., α
ilinear combination, i=1,2 ..., s, then β
iunit, note
γ
1, γ
2..., γ
sthe standard Orthogonal Vectors,
β wherein
1=α
1,
……
To two groups of bases that obtain, use transformation for mula to obtain the coordinate of millet cake cloud under new coordinate system, transformation for mula is as follows: if vectorial γ is the arbitrary vector in the n-dimensional vector SPACE V, it is at base α
1, α
2..., α
swith base β
1, β
2..., β
scoordinate be respectively X=(x
1, x
1..., x
n)
t, Y=(y
1, y
1..., y
n)
t,
Be γ=x
1α
1+ x
2α
2+ ... + x
nα
n=y
1β
1+ y
2β
2+ ... + y
nβ
n,
The vectorial coordinate transformation for mula is Y=C
-1x
Wherein C is from base α
1, α
2..., α
sto base β
1, β
2..., β
sexcessive matrix, excessively Matrix Formula is as follows:
C=(α
1,α
2,…,α
n)
T(β
1,β
2,…,β
n)。
S120: under definite shower nozzle model, shower nozzle coordinate and spray technique Parameter Conditions, take particIe system as basic description means, by Lagrangian method, calculate droplet Particles Moving track, it is as follows that Lagrange is described the Particles Moving equation:
V wherein
px, v
py, v
pzfor droplet particle rapidity vector
at x, y, the speed component on the z direction, t is the time, F
gfor droplet particle gravity
m wherein
ρfor the droplet mass particle, D is the droplet particle diameter, and ρ is the droplet particle density, and g is acceleration of gravity, F
fx, F
fy, F
fzfor the droplet particle is subject to air resistance
at x, y, the component on the z direction, air resistance
wherein μ is aerodynamic force viscosity,
for droplet Particles Moving velocity, F
wx, F
wy, F
wzfor the droplet particle is subject to the wind drag force
at x, y, the component on the z direction,
for the wind drag force
wherein
for wind vector, when indoor calm experiment, think
by time increment method, the Particles Moving in each time step is done to certain simplification, think according to every step just the motion state of particle determine the motion state of every step end particle, by the time, every step is progressively followed the tracks of, obtain the movement locus of particle in whole flight course, the droplet equation of motion is one group of nonlinear ordinary differential equation group, here adopt easy Local approximation solution, by time increment method, publicity is simplified, solution formula is as follows:
Wherein, the τ momentum relaxation time, Δ t time step,
for at n constantly
at x, y, the projection on the z direction,
for at n+1 constantly
at x, y, the projection on the z direction, n+1 Δ t constantly, in the locus of known shower nozzle coordinate (x0, y0, z0), single step spraying droplet number num situation, tried to achieve by following system of equations by the locus of each droplet point:
Wherein, x
n, y
n, z
nfor the volume coordinate position at n moment droplet particle, x
n+1, y
n+1, z
n+1for the volume coordinate position at n+1 moment droplet particle, droplet particle life cycle shows that this particle participates in the cycle of track emulation, and in the spraying simulation process, droplet particle life cycle ends up being following three kinds of situations:
(1) particle falls earthward;
(2) particle is fallen on target or on manipulating object;
(3) particle rapidity approaches 0.
S130: a kind of droplet track and plant point cloud intersection algorithm are to utilize plant point and directed line segment spatial relation, a Rule of judgment is set and realizes whether droplet drops on plant and drop point site is calculated, suppose: droplet track directed line segment BC is arranged in space, one plant point cloud A is arranged, if the vertical line that AD is triangle ABC, the D point is intersection point, by triangular nature, is obtained:
AD+BD>AB
AD+DC>AC
The addition of binomial formula obtains
2AD+DC+BD>AB+AC
That is: 2AD+BC>AB+AC
Judge whether droplet reaches in a cloud, and drop point site is thought the point of Z-direction numerical value maximum in all mid points that satisfy condition if therefore obtain meeting following formula, wherein by the distance samples method, learn that AD should be greater than 3r.
AB+AC-BC<2AD
AB+AC-BC<2AD
S140, S150, S160: constantly change the shower nozzle volume coordinate and complete dynamic similation, can simulate the Dynamic spraying process, nozzle position is changed according to the spraying operation rate request, after reaching the emulation step number Maxstep that stop condition complete requirement, calculate the sprayed deposit rate of this operation height condition lower nozzle, by repeatedly analog result comparison, realize the optimization of spray height parameter, the computing formula of droplet product rate is as follows:
Wherein C is the mist droplet deposition rate;
for the droplet mark vector, dropping on value on plant is 1, and the value on the ground that falls is 0, N
ijfor all droplet mark vector values are 1, m
ijfor Quality of fogdrop corresponding to droplet identification number.
As shown in Figure 2, be the rapeseed plants original point cloud atlas by the 3 D stereo collected by camera, have 45235 points.As shown in Figure 3, being over-sampling deletes and, through the plant point cloud chart of row-coordinate conversion, has 4574 points apart from the r=2mm sampling, and the spatial point number is reduced to original 10.1%, reduce the calculated amount of subsequent step, selecting the N point coordinate is (24.745010376 ,-158.550506592,1111.473999023), the N point of take is initial point, selecting the M point coordinate is that (32.594978333 ,-141.628906250,1091.062011719) carry out coordinate transform.As shown in Figure 4, be that shower nozzle is 0.2m in spray height, single step spraying cloud point number num is 50, and maximum emulation step number Maxstep is 100 steps, and the shower nozzle Changing Pattern is along the Dynamic spraying design sketch under x direction of principal axis 0.6m/s mobile condition.As shown in Figure 5, that shower nozzle is that 0.1m is to 0.3m in spray height, every the 0.025m altitude simulation once, single step spraying cloud point number num is 50, maximum emulation step number Maxstep is 100 steps, and the shower nozzle Changing Pattern is along the Dynamic spraying design sketch under x direction of principal axis 0.6m/s mobile condition, the relation between spray height and spraying droplet deposition, as can be seen from Figure 5, under current spraying environmental baseline, when spray height is 0.2m, the mist droplet deposition rate reaches maximal value.