The transportation network system is a special kind of general Linear Programming Problems LPs. Recently, the transportation network has been widely investigated in logistics, mathematics, transportation economics, and operations systems where distribution and transportation of commodities from town sources to final destinations is a vital issue . The job of transporter’s decisions can be optimized by reformulating the transportation Problem as generalization of the classical transportation network  .
The conventional transportation network can be formulated as a mathematical model which optimizes the cost objective function satisfying certain constraints. In two dimensional problems, transporting goods from m sources to n destinations is to be minimized. Multi-dimensional transportation network   is the large-scale extension of small classical transportation networks and is convenient for dealing with transportation problems with different supply nodes, different demand nodes as well as problems having multiple modes of transportation or delivering different kinds of goods. Thus, the proposed problem would be a more complex and large-scale problem than conventional transportation problems.
Nature based Inspired Algorithms (NIAs)   are the main source for improving intelligent systems. It also enables the decision maker to get the solutions for difficult, large-scale, and non-formulated mathematical problems. A function is optimized by Nature-based-Inspired-Algorithms by creating a candidate agent in terms of the measure of fitness (objective function). The distinguished characteristics of biogeography are maintained by Biogeography based optimization (BBO)  . Biological species can be classified according to space, and time by researching biogeography. In BBO, solution features are migrated between species which are motivated by the scientific model of biogeography.
This paper presents an improved Biogeography based optimization IBBO for solving multiobjective multi-dimensional transportation problem. The proposed approach was implemented to the flow of the commodities of the main roads to main nodes in the North Western Coastal Strip of Egypt, where there are four main roads and three nodes. Numerical results were reported in order to set up real numerical burden of the proposed system and to assess its convergence performances for solving real geographical transportation problem.
This paper is organized as follows: Optimal Transportation-Based Graph network is presented in Section 2. The Classical Transportation Problems are described in Section 3. Multiple Objective Transportation Problem is presented in Section 4. Multi-dimensional Bi-Criterion Transportation Problem is discussed in Section 5. Relaxed Multi-dimensional multi-objective transportation problem is introduced in Section 6. Section 7 presents a Modified BBO algorithm. Experimental numerical results and discussions are discussed in Section 8. Finally, the conclusion is investigated in Section 9.
2. Optimal Transportation-Based Graph Network
The analysis and solution of the multidimensional transportation problem  indicate whether a certain road should be used for transporting commodities from the perspective of optimizing transportation costs and, equivalently, if a connection between landmarks n1 and n2 should be included in the shape representation of the network. Figure 1 shows the main roads and main nodes in the north western coastline .
3. Classical Transportation Problem
The aim of the Classical Transportation network is to minimize the routing cost of goods from the supply node to the destination node. The general transportation problem can be represented as: each good is available at each of m sources and it is required that given quantities of the product be transported to each of n final destinations. The minimum cost of transporting a unit of each product from any source to any destination is known. The transportation schedule which minimizes the total cost of shipment is to be determined . The single objective, two-dimensional transportation problem can be formulated as:
: The quantity of the product to be transported from origin i to destination j.
: The available quantity at source i.
: The required quantity at destination j.
: The unit cost of transporting from source i to destination j.
The goal of transportation problem is to minimize the cost function which satisfies the constraints. The traditional objective cost is one of minimizing the total transportation cost.
Figure 1. The topological shape of the road networks in the north-western coastline .
4. Multiple Objective Transportation Problem
Multiple objective transportation optimization problem  can be defined as: “A set of decision variables which optimizes a vector cost functions which represent the cost functions, also satisfying the constraints”. Hence, the terminology “optimize” means to find a solution which would return the minimum cost values of all the cost functions acceptable to the decision maker Kasana and Kumar  formulated the multi-objective transportation two-dimensional problem as follows:
,is the required unit cost needed for shipping one unit of the product quantity from source i to the destination j satisfying l objective functions.
Let us consider m sources, and n destination, . Further, the cost of shipping one unit of goods from ith plant to jth destination is . The necessary and sufficient condition for the existence of feasible solution to balanced classical transportation problem is as follows:
The constraint & represent
equations in positive variables. Each variable appears in exactly two constraints, one is associated with source and other with destination.
5. Representation of Multi-Dimensional Bi-Criterion Transportation Problem
A conventional transportation network can be formulated in the form of a two-dimensional table for & as illustrated in Figure 2. In this figure, each cell illustrates one value of the ’s. When these values ’s are summed along the rows of the table they must equal bi, and when they are summed down the columns they must equal ai.
Figure 2. Two dimensional transportation problem.
5.1. Graphical and Mathematical Representation of Three-Dimensional Transportation Problem
Three index multi-dimensional transportation problem is a natural extension of the two-dimensional transportation type, which is explained by Haley . The multi-dimensional transportation problem in which there are m different sources, n different destinations and p of commodities type to be shipped can be represented as follows.
Three index multi-dimensional transportation problems may be represented as a block diagram, in which the layers in all directions form restricted transportation problem which is declared in Figure 3. The solid problem can be set out as a three-dimensional block for & & , in which, each cell represents one of the ’s. When these values ’s are summed along the rows (constant j and k) they equal , also, when these values are summed along the columns (constant k and i) they equal , finally, when these values ’s are summed down the heights (constant i and j) they equal .
5.2. Mathematical Representation of Four-Dimensional Transportation Problem
The mathematical representation of four-dimensional transportation problem is stated according to previous description. For easy representation, all basic notions will be illustrated and defined as follows:
Figure 3. Three index multidimensional transportation problem with m different sources and n different destinations and v different vehicles .
are m origins which represent the main sources sites.
are n destinations locations where the products are transported to them.
represent the type of vehicle (means of transportation).
represent the types of the transported product.
is the quantity of shipped goods from the origin i by vehicle k with a product l
is the quantity to be shipped to destination j by vehicle k with a product l
is the unit shipping cost from origin i to destination j by vehicle k with product l.
is the shipped quantity from source i to destination j by vehicle k with product l.
Figure 4 illustrates the representation of four-dimensional transportation problem with m different sources and n different destinations, v different vehicles and p different commodities.
5.3. Multiobjective Four-Dimensional Transportation Problem
The mathematical representation of the multiobjective four-dimensional transportation problem that has four indices, can be represented as follows:
Figure 4. Representation of four-dimensional transportation problem with m different sources and n different destinations, v different vehicles and p different commodities.
6. Relaxed Multi-Dimensional Multi-Objective Problem
The author considers some modification of the usual multi-index multi-objective transportation problem , by not allowing upper bounds for the admissible capacity for the total quantity to be sent from ith origin to the jth destination.
7. Modified BBO Algorithm
Biogeography Based Algorithm  is population-based evolutionary algorithm. It is based on investigation the biogeography science. The simulation of biogeography science declares the Migration of species (animals, birds, swarm), the development of species between habitat and Extinction of some species. In biogeography-based algorithm, habitat convenience Index is analogs to species island which is denoted as a candidate individual. This Habitat is known with high HSI are devoted to life. HSI corresponds to the BBO solution’s goodness. Rainfall, topographic diversity, temperature, land area, vegetation diversity, and others are some aspects which are included with HSI. Suitability index variables (SIVs) are known as aspects which identify habitability. Figure 5 demonstrates a species model in a single type of island.
Habitability is a term in which SIVs is the island’s independent variable and the dependent variable is HSI. High-HSI habitats are the habitat with a large number of species and Low habitats are habitats with few species. Facets of High HSI solution are given to low HSI, these facets are acquired by Low HSI provided by High HSI solution. Facets of High HSI solution is to emigrate to Low HSI solution. Emigration and Immigration tends to reform the solutions and thus emerging a solution to optimization problem. Like other evolutionary population dependent algorithms, BBO solution search procedure is an iterative procedure. After, BBO population’s initialization migration and mutation are the two type of procedures which necessitates the recited iterations. The schemes are described as.
1) Initialization of the population
The arbitrarily depressed population of habitat is created by IBBO, where every population habitat is a d-dimensional vector (the problem’s variables). represents the ith habitat in the population. Every habitat is created using the following equations.
is the lower bound of
is the upper bound of
is the rand number,
Figure 5. Species model .
In order to create the suitable design of habitat using IBBO, first, consider each solution consists of a sequence of p layers (each layer has a size of , where m is the number of origins and n is the number of different destinations). Each habitat (Figure 6) consists of p layers (p is the number of commodities). we generate each layer randomly such that
E.g., in Figure 7, there are two origins (m = 2) and we have three destination (n = 3), with 2 commodities (p = 2). In order to design the suitable structure of the solution using IBBO, first consider each individual consists of two layers (i.e., number of commodities p = 2).
Figure 6. Structure of individual for M_MOTP with m sources and n destinations and p commodities.
Figure 7. (a) and (b) Illustration of layer 1’s representation, (a) Transportation graph, (b) solution structure; (c) and (d) Illustration of layer 2's representation; (c) Transportation graph; (d) Solution structure.
We generate the 1st layer (i.e., 1st commodity) randomly such that, , , , , .
For the 2nd layer (i.e., 2nd commodity) randomly such that, , , , , .
From studying the transportation problem, it is obvious to note that all generated solution has the following distinguished features.
a) All of generated habitats are feasible.
b) The habitat length is only , that is , , for each k layer, where the problem has m sources and n destinations and p commodities.
By taking advantage of the emigration rate (µj) as well as immigration rate (λi) facets are probabilistically shared between the habitats this procedure is known as Migration . To bestowing the facets between candidate solution for modifying goodness the migration operator is liable. According to the probability of (µj) and (λi) emigration solution and immigration solutions are selected respectively.
In BBO algorithm solution’s variety are keeping the mutation  is culpable. For low and high HSI candidate solutions mutation renders a possibility for improving the solution’s goodness. It is able to intensify the solution’s quality even if they have more innumerable solutions already.
4) Pseudo-Code of the IBBO Algorithm
From the above discussion in Section 2, BBO’s pseudo-code is depicted in Figure 8.
5) Dominance Criteria
To improve the IBBO to deal with multiple objective optimization problems, the dominance criteria implemented   was in the proposed algorithm, which presented as follows:
For any optimization problem that has more than one objective function (i.e., ) any two candidate solution x1 and x2 could have one of two situations, “one dominates the other” or “none dominates the other”.
A solution x1 dominates the solution x2, if the following two conditions are verified .
“the operator refers to worse and the operator refers to better”.
a) The solution x1 is no worse than the solution x2 in all objectives.
b) The solution x1 is strictly better than the solution x2 in at least one objective function (i.e., for at least one ).
If any of these two conditions is violated, the solution x1 does not dominate the solution x2. This algorithm is repeatedly for all solutions, all solutions that
Figure 8. Pseudo code of the proposed algorithm.
are not dominated by any other solutions are constitute “the nondominated set”.
6) Archive Algorithm
The goal of the archive function  is to collect a new set of solutions in each iteration counter t, it uses the solutions in the old archive to update the solutions in the archive pool, see Figure 9. In general, the goal of this archive is to collect useful data about the underlying optimization problem during the run and update the content of the stored data (current archive content) see Figure 1.
7) Elitist Strategy
In order to guarantee rapid convergence to the nondominated-optimal solutions, elitist strategy could be implemented in the algorithm procedure. So, we propose an archiving algorithm that ensures the progress towards the nondominated-optimal solutions and at the same time covering of the whole range of the found nondominated solutions.
8. Experimental, Results and Discussions
The proposed algorithm was implemented using MATLAB version 7.0 (R14/2004). To ensure the effectiveness of the method on the multidimensional transportation problem, a numerical example was used in the computational studies.
Real application in the North Western Coastal Strip of Egypt
The following topological shape of the North Western Coastal Strip of Egypt (Figure 1) presents the geographical phenomena in a simplified manner, and does not mean the representation of all its characteristics and geographical relations; it abstracts them from some aspects that help to simplify them . The
Figure 9. Block structure of archive-based selection algorithm.
transformation of the north-west coast road network into a topological form is intended to facilitate its analysis and to make conclusions that can predict what its future status will be. The topological shape  consists of the vertices or nodes. The main nodes, the Edges or Arcs, are the roads or direct links between the nodes. The contract has gained importance for its privileged position on the network.
Table 1 gives the main roads which were investigated in our study. Figure 10 gives the daily cargo flows on the Alexandria International Highway Salloum and some major and regional roads in the Northwest Coast 2014 . The movement of goods transport and circulation is related to the study of commercial activity and its impact on transport, which mean that the role of transport and its impact on trade in the study area and the quality of goods and cargo are classified to:
Outbound trade: Carriages are transported by quarry ore, limestone, sand, gravel and child from quarries located in the province to the rest of the Republic. On the other hand, the trucks transport the production of dates, olives and fish in the province to the rest of the Republic.
Trade Imports: Products, petroleum products and water represent the most important commodities imported to the region from factories and factories. Imports also include cement. In addition, a large proportion of imports.
The nature of movable materials requires the use of a certain method, not the other, as in the case of liquid or flammable materials. Table 2 shows the movement rate of transport of goods, when examining the quality of transport between the destinations of the movement. Transport of goods on roads along the northwestern coastline.
The effectiveness of the proposed algorithm will be illustrated on the numerical example, taken from  for which .
The 1st objective cost matrix is as follows:
Figure 10. Daily cargo flows on the Alexandria International Highway Salloum and some major and regional roads in the Northwest Coast 2014 .
Table 1. The main roads in the northwestern coastal strip of Egypt.
Table 2. The type of goods transported on the roads.
The 2nd objective cost matrix is as follows:
The capacity of the ith plant (source) is and the requirement of the jth warehouse (destination) is for different commodity is as follows:
Figure 11. The detected Pareto optimal solution.
Figure 12. Curve fitted Pareto optimal solution.
Figure 11 shows the results obtained by proposed approaches versus the results obtained from I-GA , the results declare the validity of the proposed algorithm, the set of points that dominate the obtained solutions obtained by I-GA . IBBO dominate 60% of points of I-GA, on the other hand I-GA dominate 40% of points of IBBO, which illustrates the superiority of the proposed algorithm. Figure 12 represents curve fitted Pareto optimal solution, which helps the decision maker to make the right decision and select the best compromise alternative.
In this research, we present an improved evolutionary algorithm called improved biogeography-based optimization algorithm (IBBO) for solving multiobjective multi-dimensional transportation problem with an application to flow of the commodities of the main roads to main nodes in the North Western Coastal Strip of Egypt. Finally, we report numerical simulation on real application in order to assess its performances and convergence, also to set up the actual simulation burden of the proposed algorithm. The proposed method implements the concept of dominance criteria to deal with the multiple objective functions, which enable the decision maker to detect the set of all Pareto frontier of the problem which has a large-scale size. The main characteristics features of the proposed IBBO could be described as follows:
1) The detected Pareto optimal solutions are well distributed and have well satisfactory diversity characteristics.
2) Numerical results verified the superiority and the capability of the proposed algorithm.
3) This approach could also be used to treat many real applications in transportation network and urban planning.
In the future work the risk analysis will be investigated, including the best and safe road to transportation.
 Mousa, A.A. (2010) Using Genetic Algorithm and TOPSIS Technique for Multi-objective Transportation Problem: A Hybrid Approach. International Journal of Computer Mathematics, 87, 3017-3029.
 El-Shorbagy, M.A., Mousa, A.A.A., ALoraby, H. and Abo-Kila, T. (2020) Evolutionary Algorithm for Multi-Objective Multi-Index Transportation Problem under Fuzziness. Journal of Applied Research on Industrial Engineering, 7, 36-56.
 An, O., Chen, X., Li, Y., Han, Y., Zhang, J. and Shi, H. (2021) An Improved Non-Dominated Sorting Biogeography-Based Optimization Algorithm for the (Hybrid) Multi-Objective Flexible Job-Shop Scheduling Problem. Applied Soft Computing, 99, Article ID: 106869.
 Tu, J., Chen, H., Liu, J., Asghar Heidari, A., Zhang, X., Wang, M., Ruby, R. and Pham Q.-V. (2021) Evolutionary Biogeography-Based Whale Optimization Methods with Communication Structure: Towards Measuring the Balance. Knowledge-Based Systems, 212, Article ID: 106642.
 Mirjalili, S. (2019) Biogeography-Based Optimisation. In: Mirjalili, S., Ed., Evolutionary Algorithms and Neural Networks, Vol. 780, Springer, Cham, 57-72.
 Ma, H. and Simon, D. (2011) Blended Biogeography-Based Optimization for Constrained Optimization. Engineering Applications of Artificial Intelligence, 24, 517-525.
 Osman M.S., Abo-Sinna, M.A. and Mousa, A.A. (2006) IT-CEMOP: An Iterative Co-Evolutionary Algorithm for Multiobjective Optimization Problem with Nonlinear Constraints. Journal of Applied Mathematics & Computation, 183, 373-389.
 Abd Elazeem, A.E.M., Mousa, A.A.A., El-Shorbagy, M.A., Elagan, S.K. and Abo-Elnaga, Y. (2021) Detecting All Non-Dominated Points for Multi-Objective Multi-Index Transportation Problems. Sustainability, 13, Article No. 1372.
 Farrokh Ghatte, H. (2021) A Hybrid of Firefly and Biogeography-Based Optimization Algorithms for Optimal Design of Steel Frames. Arabian Journal for Science and Engineering, 46, 4703-4717.
 Farag, M.A., El-Shorbagy, M.A. Mousa, A.A. and El-Desoky, I.M. (2020) A New Hybrid Metaheuristic Algorithm for Multiobjective Optimization Problems. International Journal of Computational Intelligence Systems, 13, 920-940.
 Osman, M.S., Abo-Sinna, M.A. and Mousa, A.A. (2009) A Epsilon-Dominance Based Multiobjective Genetic Algorithm for Economic Emission Load Dispatch Optimization Problem. Electric Power Systems Research, 79, 1561-1567.
 Mousa, A.A. and Elattar, E.E. (2014) Best Compromise Alternative to EELD Problem Using Hybrid Multiobjective Quantum Genetic Algorithm. Applied Mathematics & Information Sciences, 8, 2889-2902.