Hybrid Metaheuristic-Based Spatial Modeling and Analysis of Logistics Distribution Center

: The location analysis of logistics distribution centers is one of the most critical issues in large-scale supply chains. While a number of algorithms and applications have been provided for this end, comparatively fewer investigations have been made into the integration of geographical information. This study proposes logistic distribution center location analysis that considers current geographic and embedded information gathered from a geographic information system (GIS). After reviewing the GIS, the decision variables and parameters are estimated using spatial analysis. These variables and parameters are utilized during mathematical problem-based analysis stage. While a number of existing algorithms have been proposed, this study applies a hybrid metaheuristic algorithm integrating particle swarm optimization (PSO) and genetic algorithm (GA). Using the proposed method, a more realistic mathematical model is established and solved for accurate analysis of logistics performance. To demonstrate the effectiveness of the proposed method, Korea Post distribution centers were considered in South Korea. Through tests with several real-world scenarios, it is proven experimentally that the proposed solution is more effective than existing PSO variations.


Introduction
The efficiency of logistics networks and systems in decreasing travel time and reaching longer-distance markets is considered important for enhanced economic growth [1].By acknowledging the importance of logistics activities, the World Bank helps countries measure their improvement by providing a logistics performance index (LPI).The updated LPI methodology [2] features micro-level performance data and geospatial data in their qualitative assessments, demonstrating that geospatial data plays an important role in measuring logistics performance.
Most countries consider their transportation infrastructure and logistics centers as business generators [3].The location decision of logistics centers and their performances are crucial, as they help companies not only minimize costs, traffic congestion, and environmental pollution levels, but also improve the scheduling system and vehicle routing [4][5][6][7].This is especially the case with environmental concerns since it is not only the responsibility of logistics companies but also an issue for government and other stakeholders, as well as one of customers awareness [8].
Furthermore, Europlatforms [9] has explained that logistics centers not only operate a commercial basis to benefit transport, logistics and distribution but also offer key logistical functions, such as coordination, centralization, consolidation, collaboration, and integration in a specific area's hub.Then, in the process of formulating these logistical functions, the system will be not only working for the logistics companies but also has the ability to minimize or even overcome the general issues mentioned above [10].Creating specific strategies and human research evaluation are also needed to see the system's overall efficiency [11,12].While a number of research studies have been proposed for analyzing logistic facilities, efforts to incorporate real-world data, such as geographic information system (GIS) data including facility information, have more recently begun.The integration of GIS data and existing algorithms can help enhance the performance of logistics facilities.This study proposes a hybrid metaheuristic model based on spatial analysis for measuring distribution performance of logistic centers quantitatively.In particular, post-logistics in South Korea are examined in this research.
The methodology of this study mainly consists of two parts: (a) spatial and network analysis using GIS, and (b) mathematical model-based location analysis using hybrid PSO-GA model.Hybrid of PSO-GA is effective in incorporating very large-scale instances and can be developed using various parameters.
To validate our research findings, we performed logistic center location analysis on the GIS data of Korea Post locations around the country, including population density and other information.South Korea has a total land area of 97,230 km 2 , with 81.8% of the 2020 population being urban and uneven population density [9,13].
The remainder of this paper is organized as follows.Section 2 provides relevant background knowledge and a literature review.Section 3 proposes the model and methodology with a deeper theoretical background analysis.Finally, Section 4 presents the application, and results of the proposed methodology are provided in Section 4.

Background and Literature Reviews
Rodrigue et al. [14] examined the standardization of logistics centers of various sizes based on terminal size and logistics facilities hierarchy.Moreover, selecting and evaluating the size and hierarchy of logistics centers has been one of the essential elements of business operations.Well-planned logistics center allocation helps reduce logistics costs and improve the efficiency of distribution flows [15].
Furthermore, delivery options and distribution flows are improved to welcome changes in shopping behavior, such as extending conventional home delivery to car boot delivery, in-home delivery, store collection, and more [16].It is important to note that the ability to ship goods reliably at a low cost determines a country's participation strategy in global value chains [17].
In terms of logistics costs, there are five common components, as indicated by Pohit et al. [18]: transportation costs, warehousing costs for both in-house and outsource, administration costs, inventory-holding costs, and inventory-carrying costs.In addition to considering the aforementioned roles, logistics center locations should be suitable for integrated and other transportation modes.In fact, the fluidity of transport flows and profitability are closely interlinked by transport accessibility, market coverage, land availability, and other spatial factors.
Furthermore, in reaching the optimal decisions is needed in the least amount of time.The possibility in transitioning the predictive analytics and planning of just-in-time management to develop the dynamics physical distribution can be reached by implementing the mathematical model as well as artificial intelligence [19].
For this reason, a number of existing studies have mentioned the importance of building larger and more efficient facilities to meet regional and national demands [20].The larger the size of the logistics network, the greater the complexity of its functions [21][22][23].Thus, logistics centers should tackle this issue by adopting a spatial strategy to evaluate the optimal location for maximizing service coverage [24].
Identifying suitable logistics locations can be conducted in a number of ways.Quantitative analyses such as those performed by Hagino and Endo [25] have estimated the location of facilities and distribution centers using a multinomial logit model (MNL).Many metaheuristics, including adaptive PSO, have been applied to achieve the location of logistics centers as well [26].
High performance optimization of a large-scale logistic network is carried out using a modified discrete algorithm of PSO that handles binary decision variables (BPSO) [27].Satisfactory results are also achieved considering the updated inertia weight, best-fit solutions, and worst-fit solutions from a hybrid PSO and genetic algorithm [28].This study applied an integrated PSO-genetic algorithm approach for logistics network performance analysis.Table 1 shows the results of several existing studies on metaheuristic-based logistic centers.Alternatively, spatial density analysis highlights trip distance minimization as a major factor in choosing the location of facilities in the real world [37].Similar analyses have been noticed in that the logistics activities concentration or spatial clusters are counted while having logistics relocation [38].GISs handle all the necessary physical and environmental constraints and have a significant role as a decision support tool for optimizing locations [39].The implementation of GIS also improves service quality and employee performance by providing optimum routes [40].Furthermore, GIS can compute distances in the network location model to obtain the shortest or fastest routes and solve logistics network problems.For instance, Klose and Drexl [41] proposed the location of logistics facilities using the criteria of minimizing the sum of distances between nodes and the nearest facility.
Respecting the complexity of logistics, we considered a well-formulated geographicalarea-and-mathematical-model combination for evaluating and selecting the appropriate location of logistics centers [42].Optimal clustering of numerous logistics activities can be approached by employing GIS as a spatial analyzer and BPSO as a metaheuristic model [29].Therefore, the aim of this study was to integrate GIS and metaheuristic algorithms that will manage large-scale, real-world logistics instances.In this research, a comprehensive integrated solution for the selection of logistics center locations is presented using real GIS data.

Integration of GIS Data and PSO-GA Algorithm
As mentioned previously, there are two stages used in this study to achieve an efficient solution for analyzing the logistics center location problem.While most of the research studies focused on quantitative analysis of distribution network using mathematical frameworks, the proposed method provides how GIS information is used for the quantitative analyses and their mathematical models.These incorporations and the detailed mechanism are the contributions of this study.
The proposed method is provided in Figure 1.

Integration of GIS Data and PSO-GA Algorithm
As mentioned previously, there are two stages used in this study to achieve an efficient solution for analyzing the logistics center location problem.While most of the research studies focused on quantitative analysis of distribution network using mathematical frameworks, the proposed method provides how GIS information is used for the quantitative analyses and their mathematical models.These incorporations and the detailed mechanism are the contributions of this study.
The proposed method is provided in Figure 1.The first step is a geographic analysis that involves spatial analysis and network analysis.Spatial analysis calculates the logistics densities by mapping the smaller Korea Post locations (demands) and showing the priority values of the demands.Network analysis provides the exact distances between demands.In this stage, the calculated densities and distances among the demands are considered when choosing the optimal route.The next step is to analyze the influence of the location utilizing the hybrid PSO-GA as the metaheuristic tool.There are several existing research studies that evaluate the logistics center's performance using PSO and GA separately.Table 1 shows the performances of methods using PSO and GA, and the standalone tests and combination tests using these metaheuristics are performed.As the result, the PSO-GA algorithm has better performances in this study.

Spatial and Network Analysis Using GIS
Geospatial analysis provides effective and value-based decision making and the best potential analysis using the capabilities of GIS in geospatial data management [43].In general, the spatial network structure is constructed using the shortest paths among facilities as an adjacency matrix.Additionally, the characteristics of the network should be accounted for to build a hierarchical transportation system.Moreover, the logistics center is expressed by a candidate point (on the geographic data vector) and located on the vertex nodes of the road transportation system [44].
The prepared network will then be calculated to find its Origin-Destination (OD) distance (distance between any selected demands) with distance metrics to obtain the shortest paths or time metrics to obtain the fastest paths using Dijkstra's algorithm.In addition to finding the shortest and fastest paths, spatial analysis also applies powerful The first step is a geographic analysis that involves spatial analysis and network analysis.Spatial analysis calculates the logistics densities by mapping the smaller Korea Post locations (demands) and showing the priority values of the demands.Network analysis provides the exact distances between demands.In this stage, the calculated densities and distances among the demands are considered when choosing the optimal route.The next step is to analyze the influence of the location utilizing the hybrid PSO-GA as the metaheuristic tool.There are several existing research studies that evaluate the logistics center's performance using PSO and GA separately.Table 1 shows the performances of methods using PSO and GA, and the standalone tests and combination tests using these metaheuristics are performed.As the result, the PSO-GA algorithm has better performances in this study.

Spatial and Network Analysis Using GIS
Geospatial analysis provides effective and value-based decision making and the best potential analysis using the capabilities of GIS in geospatial data management [43].In general, the spatial network structure is constructed using the shortest paths among facilities as an adjacency matrix.Additionally, the characteristics of the network should be accounted for to build a hierarchical transportation system.Moreover, the logistics center is expressed by a candidate point (on the geographic data vector) and located on the vertex nodes of the road transportation system [44].
The prepared network will then be calculated to find its Origin-Destination (OD) distance (distance between any selected demands) with distance metrics to obtain the shortest paths or time metrics to obtain the fastest paths using Dijkstra's algorithm.In addition to finding the shortest and fastest paths, spatial analysis also applies powerful statistical tools to address dependence and heterogeneity through essential event frequencies at neighboring locations (hotspots identification) [45].
The statistical tools in GIS were assigned in the initial analysis to determine the autocorrelation using the global Moran's I calculation [46].Moran's I uses the position, distance, and value of neighbor characteristics.The index value ranges from −1 to 1, where −1 indicates dispersed activities, and 1 indicates clustered activities.If a value is approximately zero, it indicates a random pattern [46].
Mathematically, Moran's I statistic for spatial autocorrelation is given in Equation (1).
where n = total number of features/spatial units; X i = the index value at location j; X = the global mean value; and, W ij = the spatial weight among location "i" and location "j".
To exhibit significant clustering or dispersion statistically, the z-score and p-value are required.As shown in Equation ( 2), z-score is the standard deviation of a feature from its mean and is based on randomization null hypothesis computation.
Used in ( 2), E[I] and V[I] are defined in Equations ( 3) and ( 4), respectively. (3) When the p-value is very, the observed spatial pattern is the result of a random process, and the null hypothesis should be rejected.On the other hand, the higher (or lower) the z-score, the more intense the clustering.As with global Moran's I index value, a near-zero zscore indicates no apparent spatial clustering [39].The decision to reject the null hypothesis is determined by the confidence level as shown in Table 2. Finally, hot spot values are calculated using Getis-Ord tools as shown in Equation ( 5) and used in the mathematical model as h i parameters.

Mathematical Model for Logistics Network Analysis
Most location allocation problems involve selecting a certain number of logistics centers, and mapping to the demand or destination points within a maximized or minimized distance.One of the optimization problems allocates facilities (logistics centers) and assigns them to the demand or destination points.This optimization also minimizes the sum of the weighted distances between all destination points and related facilities [47].The minimized weighted demand distance between the destinations and facilities is given by a set of destinations: D = (d 1 , d 2 , . . . ,d n ).Each destination (d i ) has coordinates and demands.The number of required facilities is labeled as p and the p-facility locations (x) out of the n-destinations are selected to map each destination to a facility.A binary parameter is applied to denote the destination d i mapped to facility x j .For each destination, the binary parameter is expressed by a ij as shown in Equation ( 6) and followed by the objective function cost(x) in Equation (7).Table 3 provides the parameter descriptions.
a ij 1, i f the destination i is assigned to f acility at point j 0, otherwise As shown in Table 4, this mathematical model is solved using an integrated model, including PSO and GA.PSO is a stochastic optimization-based method that uses observed social behavior and initializes a random solution through a potential solution in the problem space.

Mathematical Models
x x Presented as dimensional vectors (D), each particle is described in Equation ( 9) using the randomly generated initial velocity in Equation (10).Equation (11) defines the best local position of each particle, while Equations ( 12) and ( 13) are their adjusted personal positions and velocities according to the best solution for each particle (pbest) and the best solution of all populations (gbest) among particles in their neighborhood.
This study uses PSO because of it is simple, easy to implement, requires few parameters and has an important characteristic that GA does not have [48].One drawback of PSO is premature convergence.
On the other hand, GA, developed by John Holland, uses recombination, mutation, and selection processes as its operators.GA has been applied to wide-range optimization problems and discrete and continuous systems.GA features enable the calculation of various type of optimizations because multiple offspring act independently in a population and explore the search space simultaneously in many directions [49].Despite the advantages of GA, population size, new population and other selection parameters should be chosen wisely.Otherwise, it will be difficult to converge and may produce meaningless results.
To obtain an objective function that is superior to PSO and GA individually, we added the capability to present pbest and gbest.We expected the potential to achieve a more robust solution in a reasonable time and the detailed procedures are provided in Table 5.

Procedures Description
Step 1. Initialization (for k = 0) Set the starting value of pbest and gbest as zero.
Step 2. For i = 1 to N Calculating the solution for the first iteration.
Step 3. Generate the initial solution randomly Make an initial random population/particles as basic starting solution to generate an initial population.
Step 4. Calculate the initial solution By using the objective function, the fitness for all solutions is calculated.The objective function will maximize the profit.
Step 5. Assign pbest Get the best solution for each particle/initial position and update the pbest with the best solution (pbest').
Step 6. Assign gbest Get the best position among all particles and update the gbest with the best solution (gbest').
Step 7. Generate initial velocities Velocities are generated randomly.
Step 8. Storing the solution Store the best solution, best cost, and worst cost.
Step 9. Crossover process 1 (for k = 1: nc/2) Children are generated from all current populations (half of current populations)-Group A.
Step 10.Generate the solution The fitness for the Group A solution.
Step 11.Crossover process 2 Children are generated from all current populations with its pbest solution-Group B.
Step 12. Generate the solution The fitness for Group B solution.
Step 13.Crossover process 3 Children are generated from all current populations with its gbest solution-Group C.
Step 14. Generate the solution The fitness for Group C solution.
Step 15.Mutation process (for k = 1: nm) Select the mutation operator by a limited predetermined random rate.
Step 16.Generate the solution The fitness for mutation solution.
Step 17. Merge all solutions All solutions from current populations, Crossover groups, and mutation populations.
Step 18. Sort the best solution Select the best ones to form the next iteration population.
Step 19.Update velocities Applying the velocity limits.
Step 20.Modify the current positions Using the updated velocities.
Step 21.Generate the solution Calculating the initial solutions.
Step 22. Update pbest The best position of the ith particle.
Step 23.Update gbest The best position of the particle group.
Step 24.Finalize the algorithm For k = itermax.
Step 25.Assign gbest and stop The best solution and stop.
During this process, there were three crossover actions.The first is to generate the offspring/children from all populations to be half of the population (Group A).The second crossover process generates the children with the pbest solution to obtain the solution (Group B).The last crossover process contains the solution generated from all populations with its gbest (Group C).After performing crossover three times, the process continues to carry out the mutation and obtain the fitness solution.Finally, Groups A, B, and C and the mutation solution are sorted to obtain the best solutions for the next iteration.

Initial Analysis
For the initial evaluation, the locations of logistics centers along with the Korea Post amenities dataset from the IGIS were assigned to the study area [50].The dataset includes 1999 nodes (centers and demands) from around the country, including place name, address, amenity type, and latitude-longitudinal locations.The dataset is presented on a layer of shape-file from IGIS and is processed using QGIS 3.20 Odense software [51].
In Figures 2 and 3, the red dots represent Korea Post Distribution (KPD) center locations while the blue dots represent the small retail Korea Post locations (demands).Using the OpenStreetMap layer, Figure 2 shows the detailed environment of the dot locations in the Seoul area, and Figure 3 shows all of the demands in South Korea.Because Seoul is the largest metropolitan city in South Korea, with more than 10 million people, its logistics demand is also extremely high.This is revealed by the number of blue dots on the Korea Post locations map in comparison to the population map [13], with clustered demands shown in Figure 3.

Spatial and Network Analysis for Location Section
The model parameters for analyzing spatial analysis involve the distance matrix and logistics densities as priority values in choosing alternative locations.The distance matri can be plotted using the GIS software distance matrix toolbox.The aim matrices calculat the exact length and clusters of the network used in the mathematical model.
In this research, we complete the clustering process by dividing the study area usin The dark red hexagons in Figure 4 represent the area with the highest number of demands (blue dots) within 10 km.The number of demand dots is important for setting the network and spatial characteristics of logistics systems.
the exact length and clusters of the network used in the mathematical model.
In this research, we complete the clustering process by dividing the study area using the hot spot analysis plugin in GIS software with a hexagon-shaped grid, as shown in Figure 4.Note that the hexagon can be replaced by other shapes.The coverage area for each hexagon was 10 km and the mapped points were counted within each hexagon region.After that, the number of facilities is counted to obtain the weighted points within the radius.

Spatial and Network Analysis for Location Section
The model parameters for analyzing spatial analysis involve the distance matrix and logistics densities as priority values in choosing alternative locations.The distance matrix can be plotted using the GIS software distance matrix toolbox.The aim matrices calculate the exact length and clusters of the network used in the mathematical model.
In this research, we complete the clustering process by dividing the study area using the hot spot analysis plugin in GIS software with a hexagon-shaped grid, as shown in Figure 4.Note that the hexagon can be replaced by other shapes.The coverage area for each hexagon was 10 km and the mapped points were counted within each hexagon region.After that, the number of facilities is counted to obtain the weighted points within the radius.
The facilities locations were analyzed using the Moran's I pattern analysis.The analysis contains three types of disparities: distributed, centralized, and clustered, as shown in Figure 5. Moran's I is calculated to ensure the logistics demands in South Korea are clustered (+1).By comparing the Figures 4 and 5, we can also obtain the general picture of logistics demands in South Korea that are clustered because the darker red hexagons are located in dense areas such as Seoul, Busan, and Daegu.The QGIS analysis revealed the disparities in logistics facilities around South Korea and confirmed that facilities are clustered in regions.Using the Getis-Ord toolbox, a weighted value is generated to provide the environmental difference between the highand low-demand areas as shown in Figure 6.The estimated weighted value is then used The QGIS analysis revealed the disparities in logistics facilities around South Korea and confirmed that facilities are clustered in regions.Using the Getis-Ord toolbox, a weighted value is generated to provide the environmental difference between the highand low-demand areas as shown in Figure 6.The estimated weighted value is then used with the PSO-GA method for location analysis.The QGIS analysis revealed the disparities in logistics facilities around South Korea and confirmed that facilities are clustered in regions.Using the Getis-Ord toolbox, a weighted value is generated to provide the environmental difference between the highand low-demand areas as shown in Figure 6.The estimated weighted value is then used with the PSO-GA method for location analysis.

Location Selection Using PSO-GA
As explained in the previous section, location selection is a part of the logistics network problem in this study where "p" centers on a network and "m" nodes are considered for "n" candidate facilities, with parameters being the minimum total weighted or unweighted distances of the related.
The logistics network analysis is calculated with the weighted and unweighted areas that were estimated in the previous section.The mathematical model is solved using the

Location Selection Using PSO-GA
As explained in the previous section, location selection is a part of the logistics network problem in this study where "p" centers on a network and "m" nodes are considered for "n" candidate facilities, with parameters being the minimum total weighted or unweighted distances of the related.
The logistics network analysis is calculated with the weighted and unweighted areas that were estimated in the previous section.The mathematical model is solved using the hybrid PSO-GA method, as provided in Section 3.2.Then, the results are compared with existing methods, such as the BPSO and PSO algorithms.
As mentioned previously, the size of a logistics network problem is variously based on the number of nodes and the number of related "p" centers.This study considers five demand-number scenarios to test the model: 100, 150, 200, 350, and 600.
The demand numbers reflect the solution time differences and help us choose the proper network structure.As another performance distances travelled between demands are counted and reported.Table 6 lists the proposed parameter.
The three algorithms (the proposed hybrid PSO-GA, PSO, and BPSO) were calculated using © MATLAB 2021Ra and run on an AMD Ryzen 3 4300U laptop with Radeon Graphics, 2.70 GHz with installed 16.00 GB RAM.
The proposed hybrid PSO-GA model uses an initial population size of 100 demands and the same number is applied for the BPSO and PSO algorithms, as provided in Table 4.The iteration number was set to 500 for PSO and 500 for GA in the PSO-GA process.For comparison, the iteration number set to 1000 for both BPSO and PSO.A small maximum velocity (Vmax) determines a value of 0 and promotes exploration as a pure random search.On the other hand, a large maximum velocity indicates limited exploration.In the latter case, the maximum velocity is set to 10 for maximum and the minimum velocity is set to −10 and is updated for each iteration.Choosing the inertia factor (w) to slow particle velocities is a challenging step.If the value is less than 1, convergence is prevented.If it is set to −1 < w < 1, V ij becomes zero over time.In this experiment, all of the algorithms used 0.9 for the inertia factor value.The learning factors (c 1 and c 2 ) are both 2.0.There were no hardware or software differences in the processing of the solver operations.

PSO-GA, BPSO and PSO Comparison
The idea of the proposed algorithm is to improve the solution generation with each iteration by selection, mutation, and crossover operators.GA operators run until significant improvements are provided to the next generation.The best solutions are selected to achieve better solutions with each iteration.Therefore, the proposed PSO-GA is expected to provide superior logistics in less time.
The three algorithms PSO-GA, BPSO, and PSO are tested using two types of logistics network problems: weighted and unweighted.The weighted logistics network problem means that an area has a demand density value from GIS and, the unweighted logistics network problem means that an area has demand density with only area size.Both of the logistics network problems are applied for 100, 150, 200, 350, and 600 nodes with various numbers of medians.The results of solving for the solution time and total cost are given in Tables 7 and 8 respectively.To add complexity, the p-number is involved in the calculation for each network size using one to five instances.
Table 7 presents the solution times of each method, separated mainly based on their demand densities (weighted and unweighted).Moreover, the total cost values differ among these performances.
Consider the first case, for instance, when the network size is 100, and p is 1: the center that considers a weighted priority area calculates the total cost as $6,784,846, while the unweighted one calculates the total cost as $6,185,364.Calculations were performed repeatedly using the PSO-GA, BPSO, and PSO algorithms to validate the superior effectiveness of PSO-GA.Detailed performances and analyses are provided in the following section.

Analysis and Discussion
This study analyzes the performance of logistics centers by providing the hybrid PSO-GA algorithm, and the result is compared with the existing techniques and approaches such as BPSO and PSO to evaluate the performance of proposed method.
Reviewing the Tables 7 and 8, the solution times and total costs of the three algorithms can be compared.Using the same parameters, the hybrid PSO-GA shows better performance in terms of solution time and total cost compared to the BPSO, and Tables 9 and 10 provide a comparison for the solution times and total costs respectively.The negative value indicates that the PSO-GA has better results than BPSO or PSO.Compared to BPSO, PSO-GA offers an 80% better solution-time performance among the repetition.On the same hand, PSO-GA provides approximately 88% better solutiontime performance compared to the PSO algorithm.At a higher p-number, based on the presented results, the proposed PSO-GA shows very small differences, at some points almost 0%.
The proposed PSO-GA also provides better performance compared to BPSO in weighted logistics networks in both solution time and total cost results.The PSO-GA hybrid is 80% more effective than BPSO in solving the cost calculation.However, at a higher level of complexity, the PSO-GA and BPSO provide the same quality of performance with almost zero difference between the two solvers.In PSO-GA and BPSO comparisons, PSO-GA provides a slightly faster solution time for the weighted logistics network compared to the unweighted one.On the other hand, PSO-GA provides 88% faster solution during repetition when handling unweighted pmedian problems compare to the original PSO.A noticeable characteristic is that the larger the network size of a problems, the longer PSO needs to solve it.On the other hand, PSO solves higher-complexity problems faster.
PSO-GA solved the problems within a significantly shorter period of time compared to PSO.The proposed algorithm sorts the offspring solution and takes only the best to be reviewed as a new population for the next iteration; PSO-GA searches for the optimum solution among the best solutions, so it saves time providing the final solution.The best solutions of PSO-GA result not only from crossing the population over its pbest and gbest, but also through crossover and population mutation.
Therefore, obtaining the best solution is achieved faster than by reviewing ordinary populations.Figure 7 shows the effectiveness of the proposed model.
Despite solving the weighted logistics network 78% faster than the unweighted network, the total cost calculation demonstrates that PSO-GA provides a slightly better management cost for unweighted logistics network problems.In particular, PSO-GA offers 4.71% better average cost-saving performance for weighted problems, and 4.79% for unweighted.

2
−0.90% −3.10% −0.17% −1.66% 3 −0.12%−1.94% −0.06% −0.68% 4 0.00% −0.90% 0.00% −0.33% 5 0.00% −0.77% 0.00% −0.27% Therefore, obtaining the best solution is achieved faster than by reviewing ordinary populations.Figure 7 shows the effectiveness of the proposed model.Despite solving the weighted logistics network 78% faster than the unweighted network, the total cost calculation demonstrates that PSO-GA provides a slightly better management cost for unweighted logistics network problems.In particular, PSO-GA offers 4.71% better average cost-saving performance for weighted problems, and 4.79% for unweighted.Despite the good result, this study has limitation in terms of distribution scale.Since logistics is not restricted to the distribution between one city and another but also includes the distribution between countries, this should be considered to obtain a broad perspective in distribution processes globally.Having a positive result compared to existing studies, the proposed model is expected to be improved in the future to signify the differences and help researchers in reviewing the performance of logistics centers.

Conclusions
The selection of a logistic center requires a solution that can connect transportation nodes and minimize the transportation cost between those nodes.This study proposes the integration of a geographic information system with a hybrid metaheuristic algorithm.The first step uses a GIS-based analysis, which provides spatial and network analysis using hotspot analysis.The calculation uses Moran's index considering z-score, p-value, and Getis-Ord statistics.These values are important because they provide network priorities.These priorities are determined using demand density values.
The next step is an analyzing process using the proposed PSO-GA method.The parameters are driven from the previous GIS analysis.The proposed method provides shorter solution times for both weighted and unweighted logistics networks compared to other existing metaheuristics, such as BPSO and PSO.However, in solving larger and higher levels of network complexity, the difference in solution times between PSO-GA and BPSO is almost zero.The total cost difference between PSO-GA and BPSO is also small, but the proposed PSO-GA does provide a comparatively lower total cost.PSO-GA also demonstrates better performance in both weighted and unweighted logistics network solutions compared to PSO alone.Overall, PSO-GA is proven to be an efficient method for logistics center location analysis compared to the existing optimization models.
In this study, Seoul, South Korea, was a test site for the analysis.In future studies, much broader regions for logistic selection should be considered.Modification of the proposed method is expected.

Figure 1 .
Figure 1.The schema for the proposed method.

Figure 1 .
Figure 1.The schema for the proposed method.

Figure 3 .
Figure 3. Mapped Korea Post and logistics center locations by QGIS 3.20 Odense software.

Figure 4 .
Figure 4.The clustering analysis of logistics demand.Figure 4. The clustering analysis of logistics demand.

Figure 4 .
Figure 4.The clustering analysis of logistics demand.Figure 4. The clustering analysis of logistics demand.

Figure 6 .
Figure 6.Disparity pattern in a region, Seoul using QGIS.

Figure 6 .
Figure 6.Disparity pattern in a region, Seoul using QGIS.

Figure 7 .
Figure 7. PSO-GA, BPSO, and PSO comparison in Solution Time and Total Cost.(a) Solution Time for Weighted Network, (b) Solution Time for Unweighted Network, (c) Total Cost for Weighted Network, and (d) Total Cost for Unweighted Network.

Figure 7 .
Figure 7. PSO-GA, BPSO, and PSO comparison in Solution Time and Total Cost.(a) Solution Time for Weighted Network, (b) Solution Time for Unweighted Network, (c) Total Cost for Weighted Network, and (d) Total Cost for Unweighted Network.

Table 2 .
The uncorrected critical z-scores and p-values for different confidence levels.

Table 5 .
The procedures of the proposed PSO-GA model.

Table 7 .
Solution time (Unit: s) for weighted and unweighted center locations analysis.

Table 8 .
Total cost for weighted and unweighted logistics network.

Table 9 .
The comparative solution time comparison for weighted and unweighted networks of each solver.

Table 10 .
The comparative total cost comparison for unweighted and weighted logistics network of each solver compared to PSO-GA solver.