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.

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. subject to

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 planting

To 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 probability , with a slight bias towards local pollination

From 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:where is the pollen or solution vector at iteration and is the current best solution found among all solutions at the current generation/iteration. 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:where is the standard gamma function and this distribution is valid for large step .

Rule 2 and Rule 3 can be represented as follows:where and are pollen from different flowers of the same plant species. This essentially mimics the flower constancy in a limited neighborhood. Mathematically, if and come from the flower type or selected from the same population, this equivalently becomes a local random walk if we draw 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 to switch between common global pollination and intensive local pollination.

2.1. Pollination Intelligence Algorithm

Flower Pollination Algorithm (Algorithm 1) is limited to switching probabilities, thereby incapable of making a proper balance to search globally or locally. Experimentally, setting the probability 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 at different stages of finding the best solution:where and , . The essence of introducing this scheme is to get rid of having a constant value of throughout the iteration compared to what was presented in Rule 4 of Flower Pollinated Algorithm.

Minimize,
Initialize a population of n flowers/pollen gametes with random solutions
Find the best solutionin the initial population
Define a switch probability
Define a stopping criterion
while
  i=1 : n (all n flowers in the population)
   if,
    Draw a (d-dimensional) step vector L which obeys a Lévy distribution
    Global pollination via
   else
    Drawfrom a uniform distribution in [0, 1]
    Do local pollination via
   end if
   Evaluate new solutions
   If new solutions are better, update them in the population
  end for
   Find the current best solution
end while
Output the best solution found

Lé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:

Here , , , and are parameters that are fitted to have a promising search region. Parameter 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 , where is the current iteration, is the maximum number of iteration, and is a constant. The parameter defines how far the movement should be towards or outwards the destination. The parameter 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 is greater than 1, the current solution will be away from the best solution, whereas if has a value less than 1, the best solution will be much nearer to the current solution. Finally, the parameter 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).

Minimize,
Initialize a population of n flowers/pollen gametes with random solutions
Find the best solutionin the initial population
Perform an iterative process using equation (8) this will set values of
Define a stopping criterion
while
  i= 1 : n (all n flowers in the population)
   if,
    Define the value for,,and.
    Perform global pollination using equations (9) and (10)
   else
    Drawfrom a uniform distribution in [0, 1]
    Do local pollination via
end if
    Evaluate new solutions
    If new solutions are better, update them in the population
  end for
    Find the current best solution
end while
Output the best solution found

To 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 Tables 13. 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 13.

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 Table 4 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.

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.

Abbreviations

I:Number of selectable crops
Pi:Price of crops i
Yi:Yield of crop i
:Projected yield of crop i
Ci:Variable cost of planting crop i
:Fixed cost of cultivating the whole land
Ai:Area of land assigned to crop i
L:Total amount of land available for cultivation.

Data Availability

The data used to support the findings of the study are included within the article.

Conflicts of Interest

This research work is self-sponsored by the authors with reference materials well cited accordingly. There are no conflicts of interest regarding the publication of this paper.