# Agricultural Model for Allocation of Crops Using Pollination Intelligence Method

Authors: C. N. Ejieji; A. E. Akinsunmade
Journal: Applied Computational Intelligence and Soft Computing
 (2020)

---

## Abstract

An agricultural model for allocation of crops is considered in this work using Pollination Intelligence Method. The model was constructed to solve farmer’s decision making in allocating crops to a piece of land using market price, known yield of crops, cost incurred during planting, and the total amount of land available. A new class of metaheuristic method called Flower Pollinated Algorithm is also presented in this work to solve the designed model. An improved version of the Flower Pollinated Algorithm called Pollination Intelligence Algorithm using an iterative scheme to override the switch parameter in Flower Pollinated Algorithm is also presented and used in solving the designed model. A case study of a farmer in Ife, Osun State, Nigeria, was used to implement the model, and the results obtained suggested that instead of allocating crops to land randomly based on farmer’s intuition, cost of planting, yield of crops, and market price were factors that must be considered by farmers for optimal profit before planting crops.


---

## Body

## 1. Introduction

The roles of agriculture cannot be swept off in human life as agricultural practices have helped in provision of jobs, availability of raw materials for industrial uses, income generation for agrarian countries, to mention but a few. Some countries like Nigeria who have in time past diverted the nation’s economy from agriculture after the discovery of oil in 1956 are now in recent years revisiting and paying much interest to agricultural practices as a result of the recent decline in the price of oil. For countries like Saudi Arabia, agricultural projects are being sponsored majorly via irrigation process to help improve cropping activities in the country and this has tremendously helped with increasing the nation’s economy. It is worthwhile to note that agriculture is one of the factors for growing economy and this calls for proper attention. Land as a factor of production plays a vital role in agricultural practices, and farmers are always faced with the problem of decision making when it comes to choosing crops to be planted for a particular period and the amount of land needed for cultivation, intuition process although has been used in allocating crops to land, yet not optimally guaranteed. It is therefore still a big challenge for farmers to choose crops to be planted and the area of land to be assigned for cultivation. To achieve this, factors that enhance maximum profit such as crop yield, weather condition, price of crops influenced by market scale and demand, government policy on some crops, pest/diseases, and fertility of soil must be considered before planting crops on a piece of land, and this has to do with agricultural farm management as crop allocated to small cultivation space can eventually turn out to be the most productive, yielding huge profit amidst other selected crops planted by a farmer or group of farmers [1]. Farm management as carried out by farmers has been defined by Dillon [2] as the process by which resources and situations are manipulated by the farm manager in trying with less than full information, to achieve his or her goals. To tackle this problem, a number of optimization models have been developed by different authors, depending on their geographical zones. The first appearance of a model of such was done by Krishna [3] using known data from an Indian settlement, and this work has served as a basis and motivation for many research studies on land allocation and crop planning. Others include Keith [4] who suggested that the climate of a place is worth considering as a useful tool in maximizing techniques in farm planning, using linear programing approach. Campbell et al. [5] designed a cropland allocation model using linear programing technique and geographic information system. Sante and Crecente [6] and Wankhade and Lunge [7] also used linear programing approach as well as Ahmed et al. [8]. Meanwhile, most of these models were solved using traditional methods of optimization with few attention on the modern method of optimization. Modern methods of optimization are the most trending methods used in solving optimization problems and have been applied to virtually all fields of life including engineering, science, and finance to mention but a few. A significant number of these methods are algorithms that are derived from studying biological and environmental factors, human and animal activities, and chemical processes. Some of these methods include Genetic Algorithm (GA) developed by Holand [9]; Simulated Annealing (SA), developed by Kirkpatrick et al. [10]; Differential Evolution Method (DE) by Storn and Price [11]; and recently Firefly Algorithm (FA), an algorithm developed from studying the nature of fireflies and their light production by Yang [12]. Others include Wind Driven Optimization (WDO) by Bayraktar et al. [13], Brain Storm Optimization (BSO) by Shi [14], and Flower Pollinated Algorithm by Yang [15].In recent years, models have been developed for solving problems relating to crop selection and land management, like the work of Sante and Crecente [6], Hassan et al. [16], Tanko et al. [17], Nedunchezhian and Thirunavukkarasu [18], and Manos et al. [19]. Moreover, linear programing technique has been highlighted by Igwe et al. [20], as an efficient method of resource allocation in production planning, particularly in achieving increased agricultural productivity, and this has been applied to determine optimum enterprise combination in the work of Bamiro et al. [21]. Researchers like Wankhade and Lunge [7] have likewise made use of linear programing techniques to determine optimum allocation of crops using Simplex Algorithm and Ahmed et al. [8] to optimize the cropping pattern in Saudi Arabia, Majeke et al. [22] to mention but a few. As said earlier, the modern method of optimization has not been fully explored in solving land allocation problems. Badarudin et al. [23] gave a valid conclusion that metaheuristic approaches are good methods for solving problems relating to agriculture and land allocation. Based on this conclusion, we seek further study on the work of the abovementioned authors. In this work, agricultural model was developed for optimal crop allocation using Pollination Intelligence Algorithm which was constructed by modifying an existing metaheuristic method and used to solve the crop allocation model.


## 2. Agricultural Model for Crop Allocation

We consider a farmer or group of farmers in a region having a large portion of land say 20 hectares which is suitable for planting different types of crops. Available land must be maximized in such a way that maximum profit is obtained from all crops planted, while considering market price, history of crop yield in the region, cost of planting, and harvesting, with a set target of meeting the demand of the region. The farmer’s decision is to optimally allocate his crop of interest to a portion of land for cultivation.(1)Maximize∑i=1IPiYi−CiAi−FcL,subject to(2)∑iIAi≤L,(3)PiAi≥Yip,(4)Ai≥0.From the above model, the objective function equation (1) maximizes the production of the farmer considering market price of crops, yield of crops, cost of planting the crops, and total cost of cultivating the land. The constraint equations (2)–(4) are designed to meet the land requirement in such a way that all crops’ acreage must not be above the available land for cultivation, the crop acreage must be allocated in such a way that the farmer’s production target is met, and crop acreage must be nonnegative. Before implementing the above model, the following assumptions are made:(1)
Available land was suitable for planting farmer’s selected crop of interest(2)
Selected crops were known to be planted in farmer’s region(3)
Cost for pesticides, herbicides, fertilizer, and irrigation was treated under variable cost throughout the period of plantingTo solve this model, a new class of metaheuristic optimization method called Flower Pollinated Algorithm by Yang [15] was presented. Flower Pollinated Algorithm (FPA) as proposed by Yang [15] has the ability of solving different types of optimization problems, but has not been applied to solve land allocation models. The algorithm is based on the principle of pollination. The primary purpose of flower is ultimately reproduction. Pollen agents such as insects, birds, bats, and other animals tend to visit flowering plants having been attracted by its nature (bright colour and scent). Apart from this, abiotic factor can also be responsible for the transfer of pollen grain in flowering plants, and ten percent of the pollination process involves this process. There are over $2000$ varieties of pollinators, and they tend to behave by moving randomly. Based on this characteristic of flower pollination, Yang [15] developed Flower Pollination Algorithm (FPA). For simplicity, the following rules were used:(1)
Biotic and cross-pollination can be considered as a process of global pollination process, and pollen-carrying pollinators move in a random direction which obeys Lévy flight(2)
Self-pollination and abiotic factors were used to generalize local pollination(3)
Pollinators such as insects can develop flower constancy which is equivalent to reproduction probability that is proportional to the similarity of two flowers involved(4)
The interaction of switching of local pollination and global pollination can be controlled by a switch probabilityp∈0,1, with a slight bias towards local pollinationFrom implementation point of view, a set of updating formulas are needed to convert the rules into updating equations. In global pollination, flower pollen gametes are carried by pollinators such as insects, and pollen can travel over a long distance because insects can often fly and move in a much longer range. Rule 1 and flower constancy can be represented as follows:(5)xit+1=xit+Lλg∗−xit,where xit is the pollen or solution vector xi at iteration t and g∗ is the current best solution found among all solutions at the current generation/iteration. Lλ is the parameter that corresponds to the strength of the pollinators which is essentially the step size. Since insects may move over a long distance, Lévy flight is designed to mimic this characteristic efficiently.L > 0 is drawn from a Lévy distribution:(6)L≈λΓλsinπλ/2π1s1+λs≥s0>0,where Γλ is the standard gamma function and this distribution is valid for large step s>0.Rule 2 and Rule 3 can be represented as follows:(7)xit+1=xit+exjt−xkt,where xjt and xkt are pollen from different flowers of the same plant species. This essentially mimics the flower constancy in a limited neighborhood. Mathematically, if xjt and xkt come from the flower type or selected from the same population, this equivalently becomes a local random walk if we draw e from a uniform distribution in [0, 1]. Though flower pollination activities can occur at all scales, both local and global, adjacent flower patches or flowers in the not so far away neighborhood are more likely to be pollinated by locally. In order to mimic this, we can effectively use a switch probability (Rule 4) or proximity probability p to switch between common global pollination and intensive local pollination.


### 2.1. Pollination Intelligence Algorithm

Flower Pollination Algorithm (Algorithm1) is limited to switching probabilities, thereby incapable of making a proper balance to search globally or locally. Experimentally, setting the probability p=0.8 works well for most application, and there is therefore need of increasing the exploration in such a way that the probability goes along with the random behavior of pollinating agents. To achieve this, an iterative scheme is defined to generate the value of p at different stages of finding the best solution:(8)pn+1=θpn1−pn,where p∈0,1 and 0<θ≤4, n=0,1,2,…. The essence of introducing this scheme is to get rid of having a constant value of p throughout the iteration compared to what was presented in Rule 4 of Flower Pollinated Algorithm.Algorithm 1: Flower Pollinated Algorithm by Yang [15].
Minimizefx, x=x1,x2,…,xdInitialize a population of n flowers/pollen gametes with random solutionsFind the best solutiong∗in the initial populationDefine a switch probabilityp∈0,1Define a stopping criterionwhilet<MaxGenerationfori = 1 : n (all n flowers in the population)ifrand<p,Draw a (d-dimensional) step vector L which obeys a Lévy distributionGlobal pollination viaxit+1=xit+Lλg∗−xitelseDrawefrom a uniform distribution in [0, 1]Do local pollination viaxit+1=xit+exjt−xktend ifEvaluate new solutionsIf new solutions are better, update them in the populationend forFind the current best solutiong∗end whileOutput the best solution foundLévy flight has randomness in its step size which leads to huge variation in the position of the solution, at once it can have a long-range step and at another time it can have short steps. Due to this, Flower Pollinated Algorithm starts searching in a totally different search range and it may skip the global minima. To overcome this drawback and improve the exploration ability, Lévy flight is replaced by Sine Cosine operator. Sine Cosine Algorithm (SCA) as presented by Mirjalili [24] helps in searching promising areas of search space and we can deduce that pollinating agents tend to be distributed evenly in the search space:(9)xit+1=xit+r1sinr2r3g∗−xit,r4<0.5,(10)xit+1=xit+r1cosr2r3g∗−xit,r4≤0.5.Herer1, r2, r3, and r4 are parameters that are fitted to have a promising search region. Parameter r1 dictates the next position regions (or movement direction) which could be either in the space between the solution and destination or outside it, set as r1=T−ta/T, where t is the current iteration, T is the maximum number of iteration, and a is a constant. The parameter r2=a∗π∗rand defines how far the movement should be towards or outwards the destination. The parameter r3=a∗rand gives random weights that are being introduced for the purpose of bringing the solution either towards the best or away from the best. If the value of r3 is greater than 1, the current solution will be away from the best solution, whereas if r3 has a value less than 1, the best solution will be much nearer to the current solution. Finally, the parameter r4 is a random number between [0, 1] and it is used to switch between the Cosine and Sine components.Using equation (8) to replace the switch probability in Rule 4 of Flower Pollination Algorithm (FPA), and equation (9) and (10) in place of equation (5), we have an improved version of the algorithm called Pollination Intelligence Algorithm (PIA) which will be used together with Flower Pollinated Algorithm (FPA) and LINGO optimization modeling tool to solve the Agricultural Land Model in equations (1)–(4) (Algorithm 2).Algorithm 2: Pollination Intelligence Algorithm.
Minimizefx, x=x1,x2,…,xdInitialize a population of n flowers/pollen gametes with random solutionsFind the best solutiong∗in the initial populationPerform an iterative process using equation (8) this will set values ofpDefine a stopping criterionwhilet<MaxGenerationfori = 1: n (all n flowers in the population)ifrand<p,Define the value forr1,r2,r3andr4.Perform global pollination using equations (9) and (10)elseDrawefrom a uniform distribution in [0, 1]Do local pollination viaxit+1=xit+exjt−xktend ifEvaluate new solutionsIf new solutions are better, update them in the populationend forFind the current best solutiong∗end whileOutput the best solution foundTo implement the model, a case study of a farmer in Ife, Osun State, Nigeria, was used, data showing cost of production were collected from the farmer’s farm diary and record book, and it is presented in Tables1–3. The farmer is known to plant crops that can be grown in the region. Yield data for the selected crops were collected from http://www.nigeria.opendataforafrica.org, and the average yields (kg/acre) of the selected crops (maize, cassava, and yam) in Nigeria are 233.695, 1358.81, and 1282.305, respectively, and the market price (per kg) of the selected crops is 160, 120 and 200 which was selected by visiting Ile-Ife Market in Nigeria. The farmer’s records showing cost of planting for the selected crops are presented in Tables 1–3.Table 1
Dominion farm, Ife, Osun State, Nigeria (cost of planting maize per acre).
ItemCost (Naira)Clearing of land13,000.00Maize seedlings2,500.00Insecticide3,000.00Workmanship5,000.00Total cost of planting23,500.00Total income85,400.00Profit61,900.00Table 2
Dominion farm, Ife, Osun State, Nigeria (cost of planting cassava per acre).
ItemCost (Naira)Clearing of land13,000.00Making of ridges10,000.00Cassava stem4,000.00Workmanship6,500.00Herbicide2,500.00Total cost of planting36,000.00Total income146,000.00Profit110,000.00Table 3
Dominion farm, Ife, Osun State, Nigeria (cost of planting yam per acre).
ItemCost (Naira)Clearing of land13,000.00Making of ridges10,000.00Yam seedlings8,000.00Workmanship10,000.00Herbicide3,000.00Pesticide3,000.00Total cost of planting47,000.00Total income203,000.00Profit156,000.00


## 2.1. Pollination Intelligence Algorithm

Flower Pollination Algorithm (Algorithm1) is limited to switching probabilities, thereby incapable of making a proper balance to search globally or locally. Experimentally, setting the probability p=0.8 works well for most application, and there is therefore need of increasing the exploration in such a way that the probability goes along with the random behavior of pollinating agents. To achieve this, an iterative scheme is defined to generate the value of p at different stages of finding the best solution:(8)pn+1=θpn1−pn,where p∈0,1 and 0<θ≤4, n=0,1,2,…. The essence of introducing this scheme is to get rid of having a constant value of p throughout the iteration compared to what was presented in Rule 4 of Flower Pollinated Algorithm.Algorithm 1: Flower Pollinated Algorithm by Yang [15].
Minimizefx, x=x1,x2,…,xdInitialize a population of n flowers/pollen gametes with random solutionsFind the best solutiong∗in the initial populationDefine a switch probabilityp∈0,1Define a stopping criterionwhilet<MaxGenerationfori = 1 : n (all n flowers in the population)ifrand<p,Draw a (d-dimensional) step vector L which obeys a Lévy distributionGlobal pollination viaxit+1=xit+Lλg∗−xitelseDrawefrom a uniform distribution in [0, 1]Do local pollination viaxit+1=xit+exjt−xktend ifEvaluate new solutionsIf new solutions are better, update them in the populationend forFind the current best solutiong∗end whileOutput the best solution foundLévy flight has randomness in its step size which leads to huge variation in the position of the solution, at once it can have a long-range step and at another time it can have short steps. Due to this, Flower Pollinated Algorithm starts searching in a totally different search range and it may skip the global minima. To overcome this drawback and improve the exploration ability, Lévy flight is replaced by Sine Cosine operator. Sine Cosine Algorithm (SCA) as presented by Mirjalili [24] helps in searching promising areas of search space and we can deduce that pollinating agents tend to be distributed evenly in the search space:(9)xit+1=xit+r1sinr2r3g∗−xit,r4<0.5,(10)xit+1=xit+r1cosr2r3g∗−xit,r4≤0.5.Herer1, r2, r3, and r4 are parameters that are fitted to have a promising search region. Parameter r1 dictates the next position regions (or movement direction) which could be either in the space between the solution and destination or outside it, set as r1=T−ta/T, where t is the current iteration, T is the maximum number of iteration, and a is a constant. The parameter r2=a∗π∗rand defines how far the movement should be towards or outwards the destination. The parameter r3=a∗rand gives random weights that are being introduced for the purpose of bringing the solution either towards the best or away from the best. If the value of r3 is greater than 1, the current solution will be away from the best solution, whereas if r3 has a value less than 1, the best solution will be much nearer to the current solution. Finally, the parameter r4 is a random number between [0, 1] and it is used to switch between the Cosine and Sine components.Using equation (8) to replace the switch probability in Rule 4 of Flower Pollination Algorithm (FPA), and equation (9) and (10) in place of equation (5), we have an improved version of the algorithm called Pollination Intelligence Algorithm (PIA) which will be used together with Flower Pollinated Algorithm (FPA) and LINGO optimization modeling tool to solve the Agricultural Land Model in equations (1)–(4) (Algorithm 2).Algorithm 2: Pollination Intelligence Algorithm.
Minimizefx, x=x1,x2,…,xdInitialize a population of n flowers/pollen gametes with random solutionsFind the best solutiong∗in the initial populationPerform an iterative process using equation (8) this will set values ofpDefine a stopping criterionwhilet<MaxGenerationfori = 1: n (all n flowers in the population)ifrand<p,Define the value forr1,r2,r3andr4.Perform global pollination using equations (9) and (10)elseDrawefrom a uniform distribution in [0, 1]Do local pollination viaxit+1=xit+exjt−xktend ifEvaluate new solutionsIf new solutions are better, update them in the populationend forFind the current best solutiong∗end whileOutput the best solution foundTo implement the model, a case study of a farmer in Ife, Osun State, Nigeria, was used, data showing cost of production were collected from the farmer’s farm diary and record book, and it is presented in Tables1–3. The farmer is known to plant crops that can be grown in the region. Yield data for the selected crops were collected from http://www.nigeria.opendataforafrica.org, and the average yields (kg/acre) of the selected crops (maize, cassava, and yam) in Nigeria are 233.695, 1358.81, and 1282.305, respectively, and the market price (per kg) of the selected crops is 160, 120 and 200 which was selected by visiting Ile-Ife Market in Nigeria. The farmer’s records showing cost of planting for the selected crops are presented in Tables 1–3.Table 1
Dominion farm, Ife, Osun State, Nigeria (cost of planting maize per acre).
ItemCost (Naira)Clearing of land13,000.00Maize seedlings2,500.00Insecticide3,000.00Workmanship5,000.00Total cost of planting23,500.00Total income85,400.00Profit61,900.00Table 2
Dominion farm, Ife, Osun State, Nigeria (cost of planting cassava per acre).
ItemCost (Naira)Clearing of land13,000.00Making of ridges10,000.00Cassava stem4,000.00Workmanship6,500.00Herbicide2,500.00Total cost of planting36,000.00Total income146,000.00Profit110,000.00Table 3
Dominion farm, Ife, Osun State, Nigeria (cost of planting yam per acre).
ItemCost (Naira)Clearing of land13,000.00Making of ridges10,000.00Yam seedlings8,000.00Workmanship10,000.00Herbicide3,000.00Pesticide3,000.00Total cost of planting47,000.00Total income203,000.00Profit156,000.00


## 3. Results and Discussion

We used LINGO optimizer, Flower Pollinated Algorithm (FPA), and the new method Pollination Intelligence Algorithm (PIA) to solve the model, and the obtained results were compared with farmer’s intuition of allocating crops and profit earned from maize, cassava, and yam plantation, respectively, using 3 acres as used by the farmer. Results from Table4 show the optimal allocation of crops to agricultural land for cultivation and the profit obtained by solving the model equation using the solvers and the methods. Meanwhile, results obtained in Table 4 show that the farmland has been wasted following the farmers’ intuition process. The amount of land allocated to maize and cassava was out of hand to not guarantee the desired profit. The land would have best been maximized if yam has been given a priority following the nature of its yield and market price demand. Figures 1 and 2 show the allocation of crops and profit obtained from the model using the available number of lands used by the farmer for different crops selected based on his intuition and the optimal allocation obtained from the model using LINGO, Flower Pollinated Algorithm (FPA), and Pollination Intelligence Algorithm (PIA). Figures 1 and 2 establish that metaheuristic method of optimization can best be used in solving land allocation problems as crops were optimally allocated to the available land space for profit generation as against the result obtained from LINGO software which has the characteristics of the traditional method of optimization. The further result from these figures shows that optimal profit and allocation was achieved using the newly developed Pollination Intelligence Algorithm and solution obtained was very close to that of Flower Pollination Algorithm with a profit of 751,970 Naira and 751,300 Naira, respectively. Pollination Intelligence Algorithm (PIA) from this result displays superiority over the other methods used in solving the model.Table 4
Result of crop acreage and profit obtained from the model.
Land (acres)FarmerLingoFPAPIA3333Maize (acre)10.14605940.0123990.006299Cassava (acre)11.3688100.0101500.016452Yam (acre)11.4851312.9774502.977242Profit (Naira)327,900.00561,452.40751,300.00751,970.00Figure 1
Chart showing allocation of crops.Figure 2
Chart showing profits obtained.


## 4. Conclusion

An agricultural model for allocating crops is presented in this work. The developed model was solved using LINGO, Flower Pollination Algorithm by Yang [15], and a newly designed method called Pollination Intelligence Algorithm. Results from Table 4 gave an insight into the farmer’s decision in assigning crops of interest to his portion of land for cultivation. The Improved Flower Pollinated Algorithm called Pollination Intelligence Algorithm had a result relatively close to the solution obtained from Flower Pollinated Algorithm by Yang [15] and produces more profit from the assigned crop acreage. It is worthwhile to note that more iterative schemes can be developed to extend Flower Pollinated Algorithm and a more distributed function having the same properties with Sine Cosine operator and Lévy flight can be used to mimic the movement of pollinating agents. Furthermore, it can be deduced from the result that metaheuristic method of optimization has the ability of solving crop allocation problems, and instead of farmers using the intuition process in farm management, necessary factors must be put in place before embarking on farming activities.



---

*Source: 4830359-2020-06-01.xml*