A Heuristic Approach for Resolving Spatial Conﬂicts of Buildings in Urban Villages

: Building displacement is a common operation to resolve the spatial conﬂicts between map features, and it has important theoretical value and practical application signiﬁcance for multi-scale mapping. The prerequisite for a successful displacement operation is that there is extra space around the conﬂicting buildings into which they can be displaced. Otherwise, additional generalization operators need to be combined to resolve spatial conﬂicts. Based on this idea, this study proposes a new heuristic spatial conﬂict resolution framework that mainly resolves the spatial conﬂicts between buildings and other features in urban villages by combining three cartographic generalization operators: selection, displacement, and aggregation. This method ﬁrst reduces the density of buildings in the block through selection operation, then resolves the spatial conﬂicts between buildings and other features through displacement operation, and ﬁnally, the aggregation operation is performed to eliminate any remaining conﬂicts and newly generated conﬂicts. Experiments were carried out using real urban village data, and visual inspection and quantitative analysis were used to evaluate the experimental results. The evaluation results show that the proposed framework can not only resolve spatial conﬂicts well, but also maintain the spatial distribution and area balance of the buildings in urban villages.


Introduction
As one of the representative features in vector polygonal data, buildings have become the most common research objects in the field of map generalization.Research on map generalization related to buildings mainly includes building outline simplification [1,2], building selection [3], building group typification [4][5][6], and building arrangement pattern recognition [7][8][9][10].In the process of building generalization, due to the reduction in the map scale or the application of related generalization operators, buildings and other features (such as roads) are covered, or the distance between the features is less than the minimum visible distance specified by the cartographic constraints, which generates spatial overlap or proximity conflicts [11].To ensure the clarity of the generalized map and correctly convey knowledge and information, spatial conflicts must be resolved.Cartographers have designed multiple cartographic generalization operators to resolve spatial conflicts between map features, such as displacement [12], selection [3], aggregation [13], and typification [14].Among them, displacement is the most widely used operator because it does not significantly reduce the information content in maps [15,16].
The purpose of displacement operation is to resolve the spatial conflicts between map features while avoiding secondary conflicts.However, the prerequisite for the successful execution of the displacement operation on features is that there is extra displacement space around the conflicting features.When there is no excess displacement space around the conflicting features (e.g., urban villages), only relying on the displacement operation cannot adequately resolve the spatial conflicts, requiring a combination with other generalization operators.
Addressing the aforementioned questions, this study provides an applicable framework for resolving spatial conflicts in urban village buildings using a heuristic approach.First, constrained Delaunay triangulation is constructed after node interpolation of the building's outline.Then, the Voronoi-like polygon surrounding each building is constructed by extracting the skeleton lines of each constrained Delaunay triangle.Next, the density of the corresponding building is calculated by the area of each Voronoi-like polygon.Then, according to the calculated density value of buildings in each block, it is judged whether to perform the selection operation.This selection strategy comprehensively considers the size, shape, and orientation of the building as its importance index.On the one hand, by deleting unimportant buildings, the selection operation reduces the density of buildings in the block, and the generated space can be used for subsequent displacement operations.Moreover, spatial conflicts associated with the deleted buildings are also resolved (see Section 3.1).If the density of buildings in the block meets the given density threshold, the displacement operation is adopted to eliminate the spatial conflict between the buildings and other features.Based on the idea of combinatorial optimization, we propose a novel displacing method called Building Displacement based on Genetic Simulated Annealing Algorithm (BDGSA), which combines the advantages of a genetic algorithm (GA) and a simulated annealing algorithm (SA).It effectively overcomes the precocious phenomenon of traditional GA, and has strong spatial conflict resolution ability.Finally, the aggregation operation is used to resolve the remaining spatial conflicts after the displacement operation, mainly including unresolved spatial conflicts and newly generated spatial conflicts.
Therefore, the main contributions of this study are as follows: • A new heuristic framework for resolving spatial conflicts of building features that combines the three cartographic generalization operators of selection, displacement, and aggregation.

•
An efficient and intelligent algorithm for displacing buildings in urban villages.This algorithm is called Building Displacement based on Genetic Simulated Annealing Algorithm (BDGSA), and can be used to solve spatial conflicts between buildings and other features.
The rest of this paper is organized as follows.Section 2 presents a review of previous work on feature displacement, mainly the displacement of buildings.Section 3 describes the core steps and details of the proposed algorithm.Section 4 details the experimental results and evaluations of the method.Section 5 presents the conclusions and some suggestions for future research.

Related Works
Different displacement methods have been proposed to resolve the spatial conflicts between buildings and other features.These approaches can be divided into two categories: mechanical methods and optimization methods.
Using mechanical methods, conflicting features move sequentially according to a pre-calculated displacement order.Mackaness [17] proposed a proportional radial displacement algorithm for point features, which considers the effect of distance attenuation.At the same time, the algorithm can be extended to include other data types by sensitizing linear and polygonal features.RUAS [18] presented a displacement method for buildings that considers multiple cartographic constraints and combines multiple cartographic generalization operators.Lonergan and Jones [19] proposed an iterative building displacement method, which improves the map quality measure and search strategy.Basaraner [20] proposed a zone-based iterative building displacement method that uses spatial analysis to detect conflicting buildings in the generalization zones, and then the displacement candidate and vector are decided by means of Voronoi tessellation, spatial analysis techniques and combined multiple criteria.Sahbaz and Basaraner [21] proposed a more practical method of zonal displacement for buildings.This method first creates building groups and zones through Voronoi tessellation and buffering, then grid points are generated and weighted through kernel density estimation and buffer analysis to find appropriate positions, and finally, visual assessment and quantitative analysis are employed to measure the displacement results.Chen and Qian [22] proposed a method for regularizing buildings by combining skeleton lines and Minkowski addition.In this method, the angular bisector method is used to dissipate conflicts among the redundant vertices in the building outlines.Generally, mechanical methods solve spatial conflicts from a local perspective and will repeatedly cause secondary conflicts.
Due to the limitations of mechanical methods, optimization methods have become the focus of research.In contrast to mechanical methods, optimization methods consider multiple conflicting objects and resolve spatial conflicts simultaneously.Optimization methods can be classified into two categories: functional optimization methods and combinational optimization methods.
Functional optimization methods establish mathematical equations by combining various cartographic constraints with models in the fields of material mechanics, mathematics, and physics, and then realize the displacement of features in map generalization.Højholt [23] proposed a method for displacement based on the finite element method (FEM).This method treats displacement as the deformation of an elastic body.Some scholars have presented displacement using least-squares adjustment (LSA).Among them, Harrie and Sarjakoski [24] constructed an overdetermined equation system using generalization constraints, solving it based on least-squares adjustment.The displacement process did not end until a desired degree of accuracy was reached through the conjugate gradient method.Sester [25] regarded displacement as an optimization problem, and thus, interior constraints and exterior constraints can be formulated in terms of a least-squares adjustment approach.Some studies applied displacement based on the beams model [26] and the snake algorithm [27].However, as stated by Liu et al. [27], parameters in these approaches are difficult to determine and need to be identified using intelligent methods.From a global view, Ai et al. [16] presented an approach using a vector field and isoline model, which is effective in the maintenance of spatial relationships between buildings but may perform poorly in high-density areas.Maruyama et al. [28] addressed building displacement as a constrained optimization problem utilizing linear programming.They identified spatial relationships among map objects as constraints and optimized the cost function that penalizes the excessive displacement of buildings according to the map scale.
Combinatorial optimization methods directly use various combinatorial optimization algorithms to find an optimal solution among many candidate solutions.This kind of method has achieved satisfactory displacement effects.Ware and Jones [29] used steepest gradient descent (SGD) and simulated annealing (SA) approaches to reduce conflicts; the comparison results showed that the SA has a better ability to resolve spatial conflicts.Subsequently, Ware, Jones, and Thomas [30] improved this approach by taking generalization operations such as displacement, deletion, exaggeration, and reduction as candidates in the SA.Wilson, Ware, and Ware [31] used the genetic algorithm (GA) to realize the displacement of buildings and achieved a better displacement result.Sun et al. [32] proposed an immune genetic algorithm (IGA) that combines the self-adjustment of the individual concentration and the elite preservation methods with GA.Compared with the GA displacement results, IGA has satisfactory spatial conflict resolution results.Huang et al. [33] proposed an improved particle swarm optimization (PSO) algorithm for the displacement of buildings, which combines PSO with the cultural algorithm computational framework.These optimization algorithms all use a single population for evolution.To overcome the shortcomings of single-population evolution, Li et al. [34] proposed a building displacement method based on the multiple-population genetic algorithm (BDMPGA).This method uses an immigration operator to establish connections between the various populations to realize information exchange between various populations.The optimal displacement scheme is a result of the comprehensive evolution of multiple populations.
Through the literature summarized above, we can conclude the following: (1) Mechanical methods only use the displacement operation to displace features sequentially to resolve spatial conflicts.These kinds of methods are relatively simple and lack high intelligence.
(2) Optimization methods, such as functional optimization methods and combinatorial optimization methods, use mathematical concepts and heuristic search strategies to obtain the best displacement scheme.However, traditional optimization algorithms such as GA and SA cannot resolve all possible spatial conflicts and even generate secondary conflicts.
This study focused on resolving spatial conflicts between buildings and other features in urban villages.Its novelty lies in the proposed heuristic spatial conflict resolution framework, which overcomes previous defects using only a single displacement operator to solve spatial conflicts but adopts three cartographic generalization operators sequentially.In addition, for the ideal displacement effect, we propose the BDGSA for displacing buildings after the selection operation, which improves the deficiencies of the traditional GA and SA.The results of our experiment demonstrate the potential of this approach for the resolution of spatial conflict.

Methodologies
When the displacement space around conflicting features is insufficient, it is unrealistic to resolve spatial conflicts between the features only through displacement operation.Therefore, this study exploited a heuristic spatial conflict resolution framework by combining multiple generalization operators.Thus, the spatial conflicts between buildings and other features (i.e., buildings and roads) in urban villages are resolved by sequentially applying the three generalization operators of selection, displacement, and aggregation.Figure 1 presents the proposed framework for spatial conflict resolution of buildings in urban villages.This paper details the main steps of the framework in the following sections.

Building Selection Operation
The first step of our framework is to perform a selection operation of original buildings in the urban villages.The purpose is to reduce the building density and maintain the building density distribution in the block.The key steps of the selection operation include two parts: calculation of building density and determination of the selection strategy.
Building distribution density draws on the idea of point cluster distribution density and is defined as the distribution space occupied by each building.Based on the idea of 'spatial competition' [16], it is believed that each building exists to obtain its own living space, and the result of competition between adjacent buildings is a division of the space.This idea is consistent with the Voronoi diagram of point features; thus, the Voronoi polygons formed by constructing the Voronoi diagram for buildings are used to define the distribution density of buildings.Ai and Zhang [35] proposed the specific process of constructing Voronoi-like polygons of buildings.Constrained Delaunay triangulation is a powerful tool for detecting spatial proximity between vector features.The Voronoi-like polygon constructed based on the constrained Delaunay triangulation comprehensively considers the geometric characteristics of a building, such as the shape, size, and orientation, so that the detected spatial distance is consistent with the principle of human spatial cognition.Employing Voronoi-like polygons of buildings, the density of buildings can be defined as the reciprocal of the area of the Voronoi-like polygon surrounding each building.If the buildings are sparsely distributed, the area of the corresponding Voronoi-like polygon is larger, which means that the living space is larger and the density value is smaller; otherwise, the density value is larger.Equation (1) details the specific calculation method: where Density poly_i represents the density of the ith building, and VoronoiArea i represents the area of the Voronoi-like polygon corresponding to the ith building.

Building Selection Operation
The first step of our framework is to perform a selection operation of original buildings in the urban villages.The purpose is to reduce the building density and maintain the building density distribution in the block.The key steps of the selection operation include two parts: calculation of building density and determination of the selection strategy.
Building distribution density draws on the idea of point cluster distribution density and is defined as the distribution space occupied by each building.Based on the idea of 'spatial competition' [16], it is believed that each building exists to obtain its own living space, and the result of competition between adjacent buildings is a division of the space.This idea is consistent with the Voronoi diagram of point features; thus, the Voronoi polygons formed by constructing the Voronoi diagram for buildings are used to define the distribution density of buildings.Ai and Zhang [35] proposed the specific process of constructing Voronoi-like polygons of buildings.Constrained Delaunay triangulation is a powerful tool for detecting spatial proximity between vector features.The Voronoi-like polygon constructed based on the constrained Delaunay triangulation comprehensively considers the geometric characteristics of a building, such as the shape, size, and orientation, so that the detected spatial distance is consistent with the principle of human spatial cognition.Employing Voronoi-like polygons of buildings, the density of buildings can be defined as the reciprocal of the area of the Voronoi-like polygon surrounding each building.If the buildings are sparsely distributed, the area of the corresponding Voronoi-like This study used the density of buildings in the block as the termination condition for selection, i.e., after selection, if the density of buildings is less than the given density threshold, then the selection operation is stopped.The basic idea of the building selection operation is as follows.First, the constrained Delaunay triangulation is constructed, as presented in Figure 2b.Then, the Voronoi-like polygon is constructed based on the constructed constrained Delaunay triangulation, as presented in Figure 2c.Next, the density of each building surrounded by a Voronoi-like polygon is calculated.Then, the building with the highest density value is selected for deletion.If the building has the highest density, the area of the Voronoi-like polygon corresponding to the building is the smallest, which means that this building is less important.Next, to ensure the spatial distribution of buildings, adjacent buildings are not deleted continuously, i.e., after each building is deleted, other buildings in the first-order neighborhood are not deleted, and the building with the highest density is identified from the remaining buildings and deleted.The space occupied by the deleted building is equally divided by other nearby buildings.Therefore, after each round of deletion operation is completed, the Voronoi polygon is reconstructed, its density is calculated for the remaining buildings, and the next round of deletion operation is performed.
When the density of buildings in the block is less than a given density threshold DensityT, the selection operation is stopped.Figure 2d shows the result after the first round of the selection operation.The buildings that have been selected are shown in green.Figure 2e shows the result after the second round of the selection operation.The buildings that have been selected are shown in purple.Figure 2f shows the result after the third round of the selection operation.The buildings that have been selected are shown in cyan.The specific operation process is shown in Figure 2.
sity, the area of the Voronoi-like polygon corresponding to the building is the smallest, which means that this building is less important.Next, to ensure the spatial distribution of buildings, adjacent buildings are not deleted continuously, i.e., after each building is deleted, other buildings in the first-order neighborhood are not deleted, and the building with the highest density is identified from the remaining buildings and deleted.The space occupied by the deleted building is equally divided by other nearby buildings.Therefore, after each round of deletion operation is completed, the Voronoi polygon is reconstructed, its density is calculated for the remaining buildings, and the next round of deletion operation is performed.When the density of buildings in the block is less than a given density threshold , the selection operation is stopped.Figure 2d shows the result after the first round of the selection operation.The buildings that have been selected are shown in green.Figure 2e shows the result after the second round of the selection operation.The buildings that have been selected are shown in purple.Figure 2f shows the result after the third round of the selection operation.The buildings that have been selected are shown in cyan.The specific operation process is shown in Figure 2.

Building Displacement Operation
Our method first reduces the density of buildings by building selection operations and generates displacement spaces in each block.Then, the building displacement operation (i.e., BDGSA) is applied to resolve spatial conflicts for each building block separately.The following is a detailed introduction to the specific implementation of BDGSA from several aspects.

Population Initialization
BDGSA uses real number coding to generate chromosomes; specifically, the displacement distances of each building in the block in the X direction and Y direction are expressed as a chromosome.First, the displacement distance of each building in the X direction and the Y direction are randomly generated according to the determined maximum displacement range of the building.Then, the displacement distances of all the buildings in the block are spliced together to form a chromosome, which represents a candidate displacement scheme for the buildings in the block.Finally, randomly generated N chro-mosomes form the initial population of BDGSA for subsequent evolution.Figure 3 shows a representation of one chromosome and one population.

Population Initialization
BDGSA uses real number coding to generate chromosomes; specifically, the displacement distances of each building in the block in the X direction and Y direction are expressed as a chromosome.First, the displacement distance of each building in the X direction and the Y direction are randomly generated according to the determined maximum displacement range of the building.Then, the displacement distances of all the buildings in the block are spliced together to form a chromosome, which represents a candidate displacement scheme for the buildings in the block.Finally, randomly generated chromosomes form the initial population of BDGSA for subsequent evolution.Figure 3 shows a representation of one chromosome and one population.

Genetic Operations
The BDGSA consists of three genetic operations: chromosome selection, chromosome crossover and chromosome mutation.This study adopted the roule e selection method, the algorithmic crossover method, and the non-uniform mutation method, as detailed in this section.
The roule e selection method refers to the selection of an appropriate chromosome test strategy according to the chromosome's fitness value.Assuming that the population

Genetic Operations
The BDGSA consists of three genetic operations: chromosome selection, chromosome crossover and chromosome mutation.This study adopted the roulette selection method, the algorithmic crossover method, and the non-uniform mutation method, as detailed in this section.
The roulette selection method refers to the selection of an appropriate chromosome test strategy according to the chromosome's fitness value.Assuming that the population size is N and the fitness value of chromosome X i is f X i , the probability of chromosome X i being selected can be calculated using Equation (2).
It can be seen from the selection formula that if the fitness value of a chromosome is higher, the probability of it being selected is also greater; conversely, a chromosome with a lower fitness value is less likely to be selected.
The algorithmic crossover method is a crossover strategy adapted to real numbercoded chromosomes.Two new chromosomes are obtained through linear operations on two paired chromosomes.If the algorithmic crossover is performed between chromosome X i and chromosome X j , the two new chromosomes, X i and Y j , generated after the algorithmic crossover operation, can be calculated using Equation (3).
When performing a non-uniform mutation operation on the k-th gene x ik of the i-th chromosome, the new gene value, x ik , is determined using Equation (4).
where U k max is the maximum value of the gene, and U k min is the minimum value of the gene, f (gen) = r 2 × (1 − gen/MaxGen) 2 , where r 2 is a random number in the range of [0, 1] that meets the uniform probability distribution, gen is the current iteration number.MaxGen is the maximum iteration number, and r represents a random number in the range of [0, 1].

Simulated Annealing Operation
The purpose of introducing the idea of SA into GA is to improve the local search ability of traditional GA and then to ensure that BDGSA has strong local and global search abilities.The specific operation is after the fitness value of the chromosome in the offspring population is calculated by the fitness function, and the fitness value of each chromosome in the offspring population is compared with the corresponding chromosome in the parent population.If the offspring chromosomes are better, replace the corresponding chromosomes in the parent population with offspring chromosomes; otherwise, Metropolis sampling stability criteria are used to determine whether to update the corresponding chromosomes in the parent population with new chromosomes in the offspring population [29].

Objective Function
The basic idea of BDGSA is to integrate SA in the evolution process of GA.Therefore, similarly to the construction of the objective function and fitness function of the GA, BDGSA also considers three constraints, namely the number of remaining conflicts between buildings and buildings, the number of remaining conflicts between buildings and roads, and the total displacement distance of all buildings.Specific definitions are presented subsequently.
Regarding the number of remaining conflicts between buildings and buildings, for two buildings, B i and B j , MinDist B i , B j is the minimum distance between buildings B i and B j .con f B i , B j represents a spatial conflict between buildings B i and B j .If MinDist B i , B j is less than the predefined distance threshold, BB tol , then con f B i , B j = 1; otherwise, con f B i , B j = 0. Thus, the number of remaining conflicts between all buildings can be expressed as: Regarding the number of remaining conflicts between buildings and roads, for building B i and road R j , MinDist B i , R j is the minimum distance between building B i and road R j , and con f B i , R j represents a spatial conflict between building B i and road R j .If MinDist B i , R j is less than the predefined distance threshold, BR tol , then con f B i , R j = 1; otherwise, con f B i , R j = 0. Thus, the number of conflicts between buildings and roads can be expressed as: The total displacement distance of all buildings is measured to evaluate the positional accuracy after displacement.The displacement must be within a certain range of the positional accuracy, and the smaller the total displacement distances are, the better the positional accuracy is maintained.The total displacement distances of all the buildings are calculated by summing the displacement distances in the X and Y directions of all the buildings.The expression is as follows: Finally, the calculation formula of the objective function is as shown in Equation (8).
where w 1 , w 2 and w 3 are weights for f BB_con f , f BR_con f , and f TotalDist , respectively.The weight-setting principle of the three parameters is as follows: w 2 > w 1 > w 3 .The fitness function is obtained after the objective function undergoes a linear transformation.In this study, the reciprocal of the objective function value was taken as the fitness value, as shown in Equation (9).

Algorithm Flow
The specific implementation process of the BDGSA is shown in Figure 4.

1.
Initialize control parameters: population size, Num p , maximum number of iterations MaxGen, chromosome selection probability, P s , crossover probability, P c , mutation probability, P m , annealing initial temperature, T 0 , temperature cooling coefficient, k, and termination temperature, T end .2.
According to the coding method, randomly generate Num p chromosome s to form the initial population and use Equations ( 5) and ( 9) to calculate the fitness value, f i , of each chromosome in the population, where i = 1, 2, . . ., N.
Perform genetic operations on the population to generate a new offspring population and calculate the fitness value, f i , of each chromosome in the offspring population.

5.
Perform simulated annealing operation on the chromosomes in the offspring population, i.e., the fitness value, f i , of each chromosome in the offspring population is compared with the fitness value, f i , of the corresponding chromosome in the parent population.If f i > f i , then replace the existing chromosome with the new chromosome.
Otherwise, accept the new chromosome with probability P = exp f i − f i /T .6.
If the current temperature, T i , is less than the end temperature, T end , i.e., T i < T end , the algorithm ends and returns to the global optimal solution.Otherwise, the cooling operation is performed, i.e., T i+1 = k × T i , and transpose step (3).8.
According to the returned global optimal solution, i.e., the amount of displacement in the X direction and the amount of displacement amount in the Y direction corresponding to each building, perform the displacement operation on the building to obtain the result of the displacement of the building.

Building Aggregation Operation
Due to the limited displacement space around conflicting buildings in urban villages, the displacement operation still cannot eliminate spatial conflicts between the buildings.Moreover, it is possible that secondary conflicts will occur in local high-density areas.At this point, the aggregation operation can be used to merge the conflicting building pairs into a whole to eliminate spatial conflicts between the buildings.
According to the spatial relationship between remaining conflict buildings, the conflicts are divided into three types: intersecting relationship, tangent relationship, and separation relationship.Among them, the spatial conflict between a building with an intersecting relationship and a tangent relationship refers to the spatial overwhelming conflict, which can be resolved using the 'aggregation' tool within ArcGIS software.The spatial conflict between buildings with a separated relationship belongs to the spatial proximity conflict.For this type of spatial conflict, the conflicting buildings are first moved, and then aggregation is adopted to eliminate the spatial conflicts.The movement of conflicting buildings can squeeze out redundant blank gaps, thereby ensuring the balance of the building area of the entire block.
When moving conflicting buildings, two key parameters must be determined, namely the movement direction and the movement distance.For the movement direction, first construct a constrained Delaunay triangle network for the adjacent conflicting building pairs to obtain the visibility area.Then, extract the skeleton line of the visibility area by tracing each triangle, i.e., solid blue line in Figure 5. Finally, use least-squares multiplication to fit the extracted skeleton line, and take the normal direction of the fitted straight line as the movement direction of the conflicting building.For the movement distance, the area of conflicting buildings should be considered as the constraint condition.First, calculate the shortest distance between a pair of conflicting buildings.Then, determine the moving distance of conflicting buildings according to the area ratio.This method makes the conflicting buildings with large areas move a small distance, thus ensuring the position accuracy.The specific process is shown in Figure 5, where the red polygons represent the buildings after the movement.

Building Aggregation Operation
Due to the limited displacement space around conflicting buildings in urban villages, the displacement operation still cannot eliminate spatial conflicts between the buildings.Moreover, it is possible that secondary conflicts will occur in local high-density areas.At this point, the aggregation operation can be used to merge the conflicting building pairs into a whole to eliminate spatial conflicts between the buildings.
According to the spatial relationship between remaining conflict buildings, the conflicts are divided into three types: intersecting relationship, tangent relationship, and separation relationship.Among them, the spatial conflict between a building with an intersecting relationship and a tangent relationship refers to the spatial overwhelming conflict, which can be resolved using the 'aggregation' tool within ArcGIS software.The spatial conflict between buildings with a separated relationship belongs to the spatial proximity conflict.For this type of spatial conflict, the conflicting buildings are first moved, and then aggregation is adopted to eliminate the spatial conflicts.The movement of conflicting buildings can squeeze out redundant blank gaps, thereby ensuring the balance of the building area of the entire block.
When moving conflicting buildings, two key parameters must be determined, namely the movement direction and the movement distance.For the movement direction, first construct a constrained Delaunay triangle network for the adjacent conflicting building pairs to obtain the visibility area.Then, extract the skeleton line of the visibility area by tracing each triangle, i.e., solid blue line in Figure 5. Finally, use least-squares multiplication to fit the extracted skeleton line, and take the normal direction of the fi ed straight line as the movement direction of the conflicting building.For the movement distance, the area of conflicting buildings should be considered as the constraint condition.First, calculate the shortest distance between a pair of conflicting buildings.Then, determine the moving distance of conflicting buildings according to the area ratio.This method makes the conflicting buildings with large areas move a small distance, thus ensuring the position accuracy.The specific process is shown in Figure 5, where the red polygons represent the buildings after the movement.

Experimental Data
In this experiment, the proposed method was used to deal with spatial conflicts between buildings and roads and between buildings and buildings in urban villages.The

Experimental Data
In this experiment, the proposed method was used to deal with spatial conflicts between buildings and roads and between buildings and buildings in urban villages.The experimental data are shown in Figure 6.The experimental area was in the Chinese city of Guangzhou.The scale of the original data was 1:2000.There were 170 buildings in 4 blocks.The four blocks are marked with the numbers 1 , 2 , 3 , and 4 .It can be seen from Figure 6a that the density of buildings in each block is relatively high.When the road was widened, there were many spatial overlaps and proximity conflicts between the objects.The experimental framework was implemented by using ArcGIS 10.2 software assisted by MATLAB 2022a.

Se ing Parameters
Before the algorithm runs, the parameter values need to be set.Reasonable parameter values have a significant impact on the results of the algorithm.The relevant parameters involved in this algorithm are as follows: the selection density threshold, , the conflict distance threshold between buildings, , the conflict distance threshold between buildings and roads, , population size, , the maximum number of iterations, , selection probability, , crossover probability, , mutation probability, , initial temperature, , cooling rate, , termination temperature, , the weights of each index in the objective function, , , and .The building selection density, , was the empirical value obtained after many experiments.In this study, = 0.6 was used as the building selection density threshold, i.e., when the density of buildings in the block was less than 0.6, the selection operation was stopped.The conflict distance threshold between buildings, , was determined according to the scale with clear and legible constraints.Generally, the minimum distance threshold between objects is required to be 0.2mm according to the graphic limits used in cartography; therefore, the conflict distance threshold between buildings was determined in this experiment to be = 0.4m.The conflict distance threshold between a building and a road, , was determined by three factors: scale, legibility constraints, and the value of road widening.The finally determined conflict distance threshold between buildings and roads was = 0.9m.Population size, , refers to the number of chromosomes contained in the population.If the value of the parameter is small, the population will contain fewer chromosomes, which will shorten the running time of DBGSA, but may also lead to insufficient evolution and an inability to find the optimal solution.If the parameter takes a larger value, it will increase the computational cost of the algorithm and reduce the computational efficiency.Wilson, Ware, and Ware [31] used heuristic experiments to determine that there is a certain linear proportional relationship between the reasonable value of the parameter and the number of spatial conflicts in a block, , as detailed in Equation (10).

Setting Parameters
Before the algorithm runs, the parameter values need to be set.Reasonable parameter values have a significant impact on the results of the algorithm.The relevant parameters involved in this algorithm are as follows: the selection density threshold, DensityT, the conflict distance threshold between buildings, BB tol , the conflict distance threshold between buildings and roads, BR tol , population size, Num p , the maximum number of iterations, MaxGen, selection probability, P s , crossover probability, P c , mutation probability, P m , initial temperature, T 0 , cooling rate, k, termination temperature, T end , the weights of each index in the objective function, w 1 , w 2 , and w 3 .
The building selection density, DensityT, was the empirical value obtained after many experiments.In this study, DensityT = 0.6 was used as the building selection density threshold, i.e., when the density of buildings in the block was less than 0.6, the selection operation was stopped.The conflict distance threshold between buildings, BB tol , was determined according to the scale with clear and legible constraints.Generally, the minimum distance threshold between objects is required to be 0.2 mm according to the graphic limits used in cartography; therefore, the conflict distance threshold between buildings was determined in this experiment to be BB tol = 0.4 m.The conflict distance threshold between a building and a road, BR tol , was determined by three factors: scale, legibility constraints, and the value of road widening.The finally determined conflict distance threshold between buildings and roads was BR tol = 0.9 m.
Population size, Num p , refers to the number of chromosomes contained in the population.If the value of the parameter Num p is small, the population will contain fewer chromosomes, which will shorten the running time of DBGSA, but may also lead to insufficient evolution and an inability to find the optimal solution.If the parameter Num p takes a larger value, it will increase the computational cost of the algorithm and reduce the computational efficiency.Wilson, Ware, and Ware [31] used heuristic experiments to determine that there is a certain linear proportional relationship between the reasonable value of the parameter Num p and the number of spatial conflicts in a block, C, as detailed in Equation (10).
The maximum number of iterations, MaxGen, is the judgement condition for the termination of the algorithm.Setting a larger MaxGen parameter can obtain a more accurate optimal solution, but it will also increase the computational cost of the algorithm.At the same time, as the iteration progresses, the improvement range of the optimal solution will become smaller and smaller.Wilson, Ware, and Ware [31] used heuristic experiments to summarize the empirical formula between the maximum number of iterations, MaxGen, and the number of buildings in a block, S, as shown in Equation (11).
In the BDGSA, the crossover operation is used to generate new individuals; thus, the crossover probability, P c , should generally be a larger value.Scholars suggest that the value range of P c is [0.7, 0.9] [36].To compare with the standard GA [31], the value of P c in this experiment was 0.8.The mutation operation is also a method for generating new individuals, and its frequency is controlled by the mutation probability, P m .Scholars suggest that the value range of mutation probability, P m , is [0.001, 0.05] [36].To compare with the standard GA [31], the value of P m in the experiment was 0.008.
The initial temperature, T 0 , the temperature cooling coefficient, k, and the end temperature, T end , are the three parameters involved in the simulated annealing operation.Experiments show that the greater the initial temperature, the greater the probability of obtaining the optimal solution, but the time spent will also increase.Similarly, the smaller the temperature cooling coefficient, the slower the cooling, and although better optimization results will be obtained, it will increase the running time of the algorithm.Therefore, comprehensively considering the optimization quality and optimization efficiency, the values of the three parameters were determined using the empirical formula: initial temperature T 0 = 3, cooling rate k = 0.1, and termination temperature T end = 1.
Due to the small distance between objects in urban villages, there is a certain spatial conflict in itself.Additionally, road widening operations will aggravate the spatial conflict between roads and buildings.The principle of displacement is to maintain the position accuracy of the objects as much as possible on the premise of resolving all spatial conflicts.Therefore, the weight values in the objective function in this algorithm are as follows: w 1 = 50, w 2 = 100, and w 3 = 1.To facilitate the comparison of experimental results, the weight of the objective function in the standard GA was also the same.

Spatial Conflict Detection
In this study, the constructed constrained Delaunay triangulation and the extracted skeleton lines were used to detect the conflict area between the objects.The detection result is shown in Figure 7, where Figure 7a shows the constructed constrained Delaunay triangulation and the extracted skeleton line in the gap space between the objects, which is represented by solid red line.The red area in Figure 7b represents the conflict area between the detected objects.
It can be seen from the detected result that due to the crowded buildings in blocks, the distances between buildings and between buildings and roads are small, which cannot meet the requirements of cartographic constraints.In particular, the spatial conflict between buildings and roads is more serious, which affects the correct expression of the spatial relationship between the objects.

Experimental Results and Evaluation of the Building Selection Operation
The results of building selection in each block are shown in Figure 8. From the results, it can be observed that the proposed selection method effectively selected important buildings, reduced the density of buildings, and maintained the distribution of buildings in the block.The selection results at different densities could be obtained according to the density threshold.Therefore, this selection method is suitable for multi-scale applications.Table 1 shows the statistical information before and after the selection of each block.It can be seen from the statistical information that the total number of buildings after the selection operation is reduced to 68% of the original data, and the reduction in building density in each block is basically the same, which provides extra space for subsequent displacement operation.Due to the reduction in the number of buildings, certain spatial conflicts have been eliminated, and the number of spatial conflicts has been reduced from 131 to 80.

Experimental Results and Evaluation of the Building Selection Operation
The results of building selection in each block are shown in Figure 8. From the results, it can be observed that the proposed selection method effectively selected important buildings, reduced the density of buildings, and maintained the distribution of buildings in the block.The selection results at different densities could be obtained according to the density threshold.Therefore, this selection method is suitable for multi-scale applications.Table 1 shows the statistical information before and after the selection of each block.It can be seen from the statistical information that the total number of buildings after the selection operation is reduced to 68% of the original data, and the reduction in building density in each block is basically the same, which provides extra space for subsequent displacement operation.Due to the reduction in the number of buildings, certain spatial conflicts have been eliminated, and the number of spatial conflicts has been reduced from 131 to 80.

Experimental Results and Evaluation of the Building Selection Operation
The results of building selection in each block are shown in Figure 8. From the results, it can be observed that the proposed selection method effectively selected important buildings, reduced the density of buildings, and maintained the distribution of buildings in the block.The selection results at different densities could be obtained according to the density threshold.Therefore, this selection method is suitable for multi-scale applications.Table 1 shows the statistical information before and after the selection of each block.It can be seen from the statistical information that the total number of buildings after the selection operation is reduced to 68% of the original data, and the reduction in building density in each block is basically the same, which provides extra space for subsequent displacement operation.Due to the reduction in the number of buildings, certain spatial conflicts have been eliminated, and the number of spatial conflicts has been reduced from 131 to 80.To determine the visualization effect of the building selection results, a linear relationship was established between the building density and the grey value.Figure 9a presents a grayscale image of the building before the selection; Figure 9b shows a grayscale image of the building after the selection.It can be seen from the figure that the grey distribution of the corresponding regions before and after the selection is roughly the same.Figure 10 shows the result of spatial conflict detection on selected buildings.Figure 10a is the constructed constraint Delaunay triangulation (i.e., gray triangles) and the extracted skeleton line of the gap space (i.e., solid red line).Figure 10b shows the detected spatial conflict area, which are represented in red.Compared with the spatial conflict detection results before the building selection operation, some spatial conflicts have been eliminated.To determine the visualization effect of the building selection results, a linear relationship was established between the building density and the grey value.Figure 9a presents a grayscale image of the building before the selection; Figure 9b shows a grayscale image of the building after the selection.It can be seen from the figure that the grey distribution of the corresponding regions before and after the selection is roughly the same.Figure 10 shows the result of spatial conflict detection on selected buildings.Figure 10a is the constructed constraint Delaunay triangulation (i.e., gray triangles) and the extracted skeleton line of the gap space (i.e., solid red line).Figure 10b shows the detected spatial conflict area, which are represented in red.Compared with the spatial conflict detection results before the building selection operation, some spatial conflicts have been eliminated.To determine the visualization effect of the building selection results, a linear relationship was established between the building density and the grey value.Figure 9a presents a grayscale image of the building before the selection; Figure 9b shows a grayscale image of the building after the selection.It can be seen from the figure that the grey distribution of the corresponding regions before and after the selection is roughly the same.Figure 10 shows the result of spatial conflict detection on selected buildings.Figure 10a is the constructed constraint Delaunay triangulation (i.e., gray triangles) and the extracted skeleton line of the gap space (i.e., solid red line).Figure 10b shows the detected spatial conflict area, which are represented in red.Compared with the spatial conflict detection results before the building selection operation, some spatial conflicts have been eliminated.To determine the visualization effect of the building selection results, a linear relationship was established between the building density and the grey value.Figure 9a presents a grayscale image of the building before the selection; Figure 9b shows a grayscale image of the building after the selection.It can be seen from the figure that the grey distribution of the corresponding regions before and after the selection is roughly the same.Figure 10 shows the result of spatial conflict detection on selected buildings.Figure 10a is the constructed constraint Delaunay triangulation (i.e., gray triangles) and the extracted skeleton line of the gap space (i.e., solid red line).Figure 10b shows the detected spatial conflict area, which are represented in red.Compared with the spatial conflict detection results before the building selection operation, some spatial conflicts have been eliminated.To determine the visualization effect of the building selection results, a linear relationship was established between the building density and the grey value.Figure 9a presents a grayscale image of the building before the selection; Figure 9b shows a grayscale image of the building after the selection.It can be seen from the figure that the grey distribution of the corresponding regions before and after the selection is roughly the same.
Figure 10 shows the result of spatial conflict detection on selected buildings.Figure 10a is the constructed constraint Delaunay triangulation (i.e., gray triangles) and the extracted skeleton line of the gap space (i.e., solid red line).Figure 10b shows the detected spatial conflict area, which are represented in red.Compared with the spatial conflict detection results before the building selection operation, some spatial conflicts have been eliminated.

Experimental Results and Evaluation of the Building Displacement Operation
Figure 11a shows the experimental results of using the BDGSA to displace the buildings in each block; Figure 11b shows the superimposed effect of the BDGSA displacement results and the original building data.The blue-outlined polygons represent the buildings after displacement, and the grey polygons represent the buildings before displacement.Table 2 details the statistical information of the BDGSA displacement results, including the number of spatial conflicts between buildings before and after the displacement, the number of spatial conflicts between buildings and roads, and the total displacement distance of buildings in each block.The statistical information can be utilized to evaluate the displacement effect of the BDGSA.To determine the visualization effect of the building selection results, a linear relationship was established between the building density and the grey value.Figure 9a presents a grayscale image of the building before the selection; Figure 9b shows a grayscale image of the building after the selection.It can be seen from the figure that the grey distribution of the corresponding regions before and after the selection is roughly the same.Figure 10 shows the result of spatial conflict detection on selected buildings.Figure 10a is the constructed constraint Delaunay triangulation (i.e., gray triangles) and the extracted skeleton line of the gap space (i.e., solid red line).Figure 10b shows the detected spatial conflict area, which are represented in red.Compared with the spatial conflict detection results before the building selection operation, some spatial conflicts have been eliminated.

Experimental Results and Evaluation of the Building Displacement Operation
Figure 11a shows the experimental results of using the BDGSA to displace the buildings in each block; Figure 11b shows the superimposed effect of the BDGSA displacement results and the original building data.The blue-outlined polygons represent the buildings after displacement, and the grey polygons represent the buildings before displacement.Table 2 details the statistical information of the BDGSA displacement results, including the number of spatial conflicts between buildings before and after the displacement, the number of spatial conflicts between buildings and roads, and the total displacement distance of buildings in each block.The statistical information can be utilized to evaluate the displacement effect of the BDGSA.

Experimental Results and Evaluation of the Building Displacement Operation
Figure 11a shows the experimental results of using the BDGSA to displace the buildings in each block; Figure 11b shows the superimposed effect of the BDGSA displacement results and the original building data.The blue-outlined polygons represent the buildings after displacement, and the grey polygons represent the buildings before displacement.Table 2 details the statistical information of the BDGSA displacement results, including the number of spatial conflicts between buildings before and after the displacement, the number of spatial conflicts between buildings and roads, and the total displacement distance of buildings in each block.The statistical information can be utilized to evaluate the displacement effect of the BDGSA.To determine the visualization effect of the building selection results, a linear relationship was established between the building density and the grey value.Figure 9a presents a grayscale image of the building before the selection; Figure 9b shows a grayscale image of the building after the selection.It can be seen from the figure that the grey distribution of the corresponding regions before and after the selection is roughly the same.To determine the visualization effect of the building selection results, a linear relationship was established between the building density and the grey value.Figure 9a presents a grayscale image of the building before the selection; Figure 9b shows a grayscale image of the building after the selection.It can be seen from the figure that the grey distribution of the corresponding regions before and after the selection is roughly the same.To determine the visualization effect of the building selection results, a linear relationship was established between the building density and the grey value.Figure 9a presents a grayscale image of the building before the selection; Figure 9b shows a grayscale image of the building after the selection.It can be seen from the figure that the grey distribution of the corresponding regions before and after the selection is roughly the same.To determine the visualization effect of the building selection results, a linear relationship was established between the building density and the grey value.Figure 9a presents a grayscale image of the building before the selection; Figure 9b shows a grayscale image of the building after the selection.It can be seen from the figure that the grey distribution of the corresponding regions before and after the selection is roughly the same.To better illustrate the displacement results of the BDGSA, we selected some examples of residual conflicts between buildings in each block, which are marked with black circles A-I in Figure 11. Figure 12 shows the detailed information of examples of residual conflicts, which the blue polygons represent the buildings after displacement.Statistical analysis revealed that there were three main spatial relationships between the remaining conflict objects: separation, tangency, and intersection.Figure 12A-C show the remaining conflict types with detached relationships; Figure 12D-F show the remaining conflict types with tangent relationships; and Figure 12G-I show the remaining conflict types with intersecting relationships.Spatial conflicts with detached relationships belong to spatial proximity conflicts, while spatial conflicts with tangent and intersecting relationships belong to spatial overlap conflicts.To be er illustrate the displacement results of the BDGSA, we selected some examples of residual conflicts between buildings in each block, which are marked with black circles A-I in Figure 11. Figure 12 shows the detailed information of examples of residual conflicts, which the blue polygons represent the buildings after displacement.Statistical analysis revealed that there were three main spatial relationships between the remaining conflict objects: separation, tangency, and intersection.Figure 12A-C show the remaining conflict types with detached relationships; Figure 12D-F show the remaining conflict types with tangent relationships; and Figure 12G-I show the remaining conflict types with intersecting relationships.Spatial conflicts with detached relationships belong to spatial proximity conflicts, while spatial conflicts with tangent and intersecting relationships belong to spatial overlap conflicts.

Experimental Results and Evaluation of the Building Aggregation Operation
The aggregation operation was used to resolve the remaining spatial conflicts in this study.The steps used to achieve the aggregation operation were also different according to the type of spatial conflict; there was excess space between conflicting buildings with separate relationships.To maintain balance in the area of the buildings in the block, the conflicting building pairs were moved first and then aggregated.Figure 13 shows the movement effect of conflicting building pairs before they are aggregated, which red polygons represent the buildings after displacement and blue lines represent the extracted skeleton lines.The movement squeezed out any excess space between the conflicting building pairs.to the type of spatial conflict; there was excess space between conflicting buildings with separate relationships.To maintain balance in the area of the buildings in the block, the conflicting building pairs were moved first and then aggregated.Figure 13 shows the movement effect of conflicting building pairs before they are aggregated, which red polygons represent the buildings after displacement and blue lines represent the extracted skeleton lines.The movement squeezed out any excess space between the conflicting building pairs.Figure 14 shows the effect of first moving and then aggregating the remaining spatial conflicts with separate relationships in each block, and the effect of directly using the aggregation operation within ArcGIS software for the remaining spatial conflicts with tangent and intersecting relationships.In comparison with the result after displacement, it was found that after the aggregation operation, the remaining spatial conflicts between buildings in each block were all resolved.Figure 14 shows the effect of first moving and then aggregating the remaining spatial conflicts with separate relationships in each block, and the effect of directly using the aggregation operation within ArcGIS software for the remaining spatial conflicts with tangent and intersecting relationships.In comparison with the result after displacement, it was found that after the aggregation operation, the remaining spatial conflicts between buildings in each block were all resolved.
The aggregation operation was used to resolve the remaining spatial conflicts in this study.The steps used to achieve the aggregation operation were also different according to the type of spatial conflict; there was excess space between conflicting buildings with separate relationships.To maintain balance in the area of the buildings in the block, the conflicting building pairs were moved first and then aggregated.Figure 13 shows the movement effect of conflicting building pairs before they are aggregated, which red polygons represent the buildings after displacement and blue lines represent the extracted skeleton lines.The movement squeezed out any excess space between the conflicting building pairs.Figure 14 shows the effect of first moving and then aggregating the remaining spatial conflicts with separate relationships in each block, and the effect of directly using the aggregation operation within ArcGIS software for the remaining spatial conflicts with tangent and intersecting relationships.In comparison with the result after displacement, it was found that after the aggregation operation, the remaining spatial conflicts between buildings in each block were all resolved.

Discussion
To better evaluate the spatial conflict resolution ability of the BDGSA proposed in this study, we implemented two other intelligent optimization algorithms to perform comparative experiments.Because the BDGSA incorporates the idea of SA into the traditional GA, the first comparative method is the standard GA, which was used by Wilson [31] to displace buildings.The other comparative method is the SA, which was proposed by Ware [29] to reduce conflicts.The experimental results of the BDGSA, GA, and SA are compared from the two aspects of the legibility constraint and the positional accuracy constraint.
The legibility constraint uses the number of remaining spatial conflicts as the evaluation index, which indicates the ability of the algorithm to resolve spatial conflicts.The smaller this value is, the stronger the ability of the algorithm to resolve spatial conflicts is.The positional accuracy constraint uses the total displacement distance of all buildings as the evaluation index, which indicates the disruption degree of the positional accuracy.The smaller the value is, the better the positional accuracy is.
Figure 15a shows the displacement results based on the standard GA; Figure 15b shows the overlay display effect before and after building displacement using the GA.In the figure, the cyan-outlined polygons represent the displaced buildings using GA, and the gray polygons represent the buildings before displacement.The black circles A-I represent some examples of residual conflicts, which correspond to the one in Figure 11.The statistical information of the displacement results based on the GA is listed in Table 3.
From Table 3, it can be seen that the GA cannot resolve all the spatial conflicts between buildings and other features.In particular, all spatial conflicts have been resolved in block 1.There are three unresolved spatial conflicts between buildings and buildings and one spatial conflict between buildings and roads in block 2. The number and type of remaining spatial conflicts in block 3 are the same as those in block 2. There are five unresolved spatial conflicts between buildings in block 4.

Discussion
To be er evaluate the spatial conflict resolution ability of the BDGSA proposed in this study, we implemented two other intelligent optimization algorithms to perform comparative experiments.Because the BDGSA incorporates the idea of SA into the traditional GA, the first comparative method is the standard GA, which was used by Wilson [31] to displace buildings.The other comparative method is the SA, which was proposed by Ware [29] to reduce conflicts.The experimental results of the BDGSA, GA, and SA are compared from the two aspects of the legibility constraint and the positional accuracy constraint.The legibility constraint uses the number of remaining spatial conflicts as the evaluation index, which indicates the ability of the algorithm to resolve spatial conflicts.The smaller this value is, the stronger the ability of the algorithm to resolve spatial conflicts is.The positional accuracy constraint uses the total displacement distance of all buildings as the evaluation index, which indicates the disruption degree of the positional accuracy.The smaller the value is, the be er the positional accuracy is.
Figure 15a shows the displacement results based on the standard GA; Figure 15b shows the overlay display effect before and after building displacement using the GA.In the figure, the cyan-outlined polygons represent the displaced buildings using GA, and the gray polygons represent the buildings before displacement.The black circles A-I represent some examples of residual conflicts, which correspond to the one in Figure 11.The basic statistical information of the displacement results based on the GA is listed in Table 3. From Table 3, it can be seen that the GA cannot resolve all the spatial conflicts between buildings and other features.In particular, all spatial conflicts have been resolved in block 1.There are three unresolved spatial conflicts between buildings and buildings and one spatial conflict between buildings and roads in block 2. The number and type of remaining spatial conflicts in block 3 are the same as those in block 2. There are five unresolved spatial conflicts between buildings in block 4.    To determine the visualization effect of the building selection results, a linear relationship was established between the building density and the grey value.Figure 9a presents a grayscale image of the building before the selection; Figure 9b shows a grayscale image of the building after the selection.It can be seen from the figure that the grey distribution of the corresponding regions before and after the selection is roughly the same.To determine the visualization effect of the building selection results, a linear relationship was established between the building density and the grey value.Figure 9a presents a grayscale image of the building before the selection; Figure 9b shows a grayscale image of the building after the selection.It can be seen from the figure that the grey distribution of the corresponding regions before and after the selection is roughly the same.To determine the visualization effect of the building selection results, a linear relationship was established between the building density and the grey value.Figure 9a presents a grayscale image of the building before the selection; Figure 9b shows a grayscale image of the building after the selection.It can be seen from the figure that the grey distribution of the corresponding regions before and after the selection is roughly the same.To determine the visualization effect of the building selection results, a linear relationship was established between the building density and the grey value.Figure 9a presents a grayscale image of the building before the selection; Figure 9b shows a grayscale image of the building after the selection.It can be seen from the figure that the grey distribution of the corresponding regions before and after the selection is roughly the same.Figure 16a shows the displacement results based on SA; Figure 16b shows the overlay display effect before and after building displacement using SA.In the figure, the pinkoutlined polygons represent the displaced buildings using SA; the grey polygons represent the buildings before the displacement.The black circles A-I represent some examples of residual conflicts, which correspond to the one in Figure 11.The basic statistical information of the displacement results based on the SA is listed in Table 4. Table 4 shows that the SA had one unresolved spatial conflict between buildings located in block 2. Figure 16a shows the displacement results based on SA; Figure 16b shows the overlay display effect before and after building displacement using SA.In the figure, the pinkoutlined polygons represent the displaced buildings using SA; the grey polygons represent the buildings before the displacement.The black circles A-I represent some examples of residual conflicts, which correspond to the one in Figure 11.The basic statistical information of the displacement results based on the SA is listed in Table 4. Table 4 shows that the SA had one unresolved spatial conflict between buildings located in block 2.  Although these methods can be used to displace buildings and resolve spatial conflicts, they yield different displacement results.The statistical information above indicates the following: 1.For the legibility constraint, the total numbers of remaining spatial conflicts after building displacement using the BDGSA, GA, and SA are 10, 13, and 4, respectively, and the proportions of conflict resolution are 87.5%, 83.75%, and 95%, respectively.Therefore, in terms of the ability to resolve spatial conflicts independently, SA is the best, followed by the BDGSA.The reason for this situation is that the objective function of the BDGSA is constructed with the number of spatial conflicts and the total displacement distance as constraints, whereas the objective function of SA is constructed only with the number of spatial conflicts as a constraint.The different constraints considered in the objective function led to the much larger displacement distance of SA than that of the BDGSA.To determine the visualization effect of the building selection results, a linear relationship was established between the building density and the grey value.Figure 9a presents a grayscale image of the building before the selection; Figure 9b shows a grayscale image of the building after the selection.It can be seen from the figure that the grey distribution of the corresponding regions before and after the selection is roughly the same.To determine the visualization effect of the building selection results, a linear relationship was established between the building density and the grey value.Figure 9a presents a grayscale image of the building before the selection; Figure 9b shows a grayscale image of the building after the selection.It can be seen from the figure that the grey distribution of the corresponding regions before and after the selection is roughly the same.To determine the visualization effect of the building selection results, a linear relationship was established between the building density and the grey value.Figure 9a presents a grayscale image of the building before the selection; Figure 9b shows a grayscale image of the building after the selection.It can be seen from the figure that the grey distribution of the corresponding regions before and after the selection is roughly the same.To determine the visualization effect of the building selection results, a linear relationship was established between the building density and the grey value.Figure 9a presents a grayscale image of the building before the selection; Figure 9b shows a grayscale image of the building after the selection.It can be seen from the figure that the grey distribution of the corresponding regions before and after the selection is roughly the same.Although these methods can be used to displace buildings and resolve spatial conflicts, they yield different displacement results.The statistical information above indicates the following: 1.
For the legibility constraint, the total numbers of remaining spatial conflicts after building displacement using the BDGSA, GA, and SA are 10, 13, and 4, respectively, and the proportions of conflict resolution are 87.5%, 83.75%, and 95%, respectively.Therefore, in terms of the ability to resolve spatial conflicts independently, SA is the best, followed by the BDGSA.The reason for this situation is that the objective function of the BDGSA is constructed with the number of spatial conflicts and the total displacement distance as constraints, whereas the objective function of SA is constructed only with the number of spatial conflicts as a constraint.The different constraints considered in the objective function led to the much larger displacement distance of SA than that of the BDGSA.

2.
For the positional accuracy constraint, the total displacement distances of the BDGSA, GA and SA are 20.895m, 28.974 m, and 95.334 m, respectively.It can be seen from the statistical information that the total displacement distance of the BDGSA proposed in this study is the shortest, while the total displacement distance of SA is the largest.Therefore, the proposed BDGSA can better maintain the positional accuracy constraints of the displaced buildings, while the SA cannot.
Therefore, considering the two constraints of legibility constraint and positional accuracy constraint, the BDGSA achieves a better displacement effect.The main reason is that BDGSA introduces the idea of SA into the traditional GA and realizes the complementary advantages of GA and SA.BDGSA has strong local research ability and global search ability.Therefore, as the experimental results show, compared with GA and SA, BDGSA has a stronger comprehensive ability to solve spatial conflicts.
Successful displacement should not only solve all possible spatial conflicts after displacing some buildings but also ensure the positional accuracy of the displaced buildings.
To better demonstrate the difference between the proposed BDGSA and the other two comparative methods in maintaining the positional accuracy of displaced buildings, we calculated the detailed statistical information, and the relevant statistical indicators of the building displacement distances in each block of the three methods of BDGSA, GA, and SA, as shown in Tables 5-7, which indicate the following: 1.
The statistical information that the minimum amount of displacement of buildings using the BDGSA in each block is 0.000 m, while the minimum displacements of the GA and SA are greater than 0 m.This suggests that if some buildings in the BDGSA do not conflict with any neighboring objects, they will not be moved to maintain position accuracy.If the area of a conflicting building is larger, it indicates that the building is relatively important, and its location accuracy should be ensured as much as possible without displacement.

2.
The average displacement distance of buildings in each block of the BDGSA is smaller than the corresponding blocks of the GA and SA, which shows that the BDGSA displacement method has a shorter displacement distance.Compared with the BDGSA and GA, the SA exhibits the largest average displacement value in each block, which shows that the SA is weak in maintaining the accuracy of the building displacement.
The main reason is that SA tries to designate displacement candidate positions in a continuous displacement space; therefore, the size of the displacement distance is greatly affected by the subsequent positions.To determine the visualization effect of the building selecti tionship was established between the building density and the g sents a grayscale image of the building before the selection; Figu image of the building after the selection.It can be seen from the tribution of the corresponding regions before and after the select To determine the visualization effect of the building selecti tionship was established between the building density and the g sents a grayscale image of the building before the selection; Figu image of the building after the selection.It can be seen from the tribution of the corresponding regions before and after the select To determine the visualization effect of the building selecti tionship was established between the building density and the g sents a grayscale image of the building before the selection; Figu image of the building after the selection.It can be seen from the tribution of the corresponding regions before and after the select To determine the visualization effect of the building selecti tionship was established between the building density and the g sents a grayscale image of the building before the selection; Figu image of the building after the selection.It can be seen from the tribution of the corresponding regions before and after the select To determine the visualization effect of the building selecti tionship was established between the building density and the g sents a grayscale image of the building before the selection; Figu image of the building after the selection.It can be seen from the tribution of the corresponding regions before and after the select To determine the visualization effect of the building selecti tionship was established between the building density and the g sents a grayscale image of the building before the selection; Figu image of the building after the selection.It can be seen from the tribution of the corresponding regions before and after the select   To determine the visualization effect of the building selecti tionship was established between the building density and the g sents a grayscale image of the building before the selection; Figu image of the building after the selection.It can be seen from the tribution of the corresponding regions before and after the select To objectively evaluate the maintenance of the building's relative position, this study compares the changes in the area of Voronoi polygons before and after building displacement.The Voronoi polygon of one building represents its competition with the surrounding space.When the relative position between adjacent buildings changes, the size of the corresponding Voronoi polygons will also change.Therefore, by analyzing the size of changes in Voronoi polygons before and after buildings displacement, the effect of maintaining the spatial distribution and relative positional relationship of each building can be quantitatively evaluated.Figure 17a depicts the Voronoi polygons of the buildings in each block before the displacement, i.e., red polygons, and Figure 17b is the Voronoi polygons of buildings in each block after the displacement, i.e., red polygons.Visual observation revealed that the Voronoi polygons of the buildings before and after the displacement are very similar in structure.This indicates that the BDGSA can sufficiently maintain the spatial distribution of buildings after displacement.To determine the visualization effect of the building selecti tionship was established between the building density and the g sents a grayscale image of the building before the selection; Figu image of the building after the selection.It can be seen from the tribution of the corresponding regions before and after the select  To determine the visualization effect of the building selecti tionship was established between the building density and the g sents a grayscale image of the building before the selection; Figu image of the building after the selection.It can be seen from the tribution of the corresponding regions before and after the select  To determine the visualization effect of the building selecti tionship was established between the building density and the g sents a grayscale image of the building before the selection; Figu image of the building after the selection.It can be seen from the tribution of the corresponding regions before and after the select   To determine the visualization effect of the building selecti tionship was established between the building density and the g sents a grayscale image of the building before the selection; Figu image of the building after the selection.It can be seen from the tribution of the corresponding regions before and after the select   To objectively evaluate the maintenance of the building's relative position, this study compares the changes in the area of Voronoi polygons before and after building displacement.The Voronoi polygon of one building represents its competition with the surrounding space.When the relative position between adjacent buildings changes, the size of the corresponding Voronoi polygons will also change.Therefore, by analyzing the size of changes in Voronoi polygons before and after buildings displacement, the effect of maintaining the spatial distribution and relative positional relationship of each building can be quantitatively evaluated.Figure 17a depicts the Voronoi polygons of the buildings in each block before the displacement, i.e., red polygons, and Figure 17b is the Voronoi polygons of buildings in each block after the displacement, i.e., red polygons.Visual observation revealed that the Voronoi polygons of the buildings before and after the displacement are very similar in structure.This indicates that the BDGSA can sufficiently maintain the spatial distribution of buildings after displacement.For quantitative analysis, Figure 18 presents the areas of Voronoi polygons before and after the displacement of the buildings of block ①, block ②, block ③, and block ④ .The abscissa represents the area of the Voronoi polygon of the building before the displacement; the ordinate represents the area of the Voronoi polygon of the building after the displacement.The correlation coefficients, R 2 , of each block are 0.8623, 0.9747, 0.9771, and 0.9744, indicating that the area of the Voronoi polygons before and after the building displacement is highly correlated.The quantitative analysis results show that the displacement operation can adequately maintain the spatial distribution and relative position of the building after displacement.For quantitative analysis, Figure 18 presents the areas of Voronoi polygons before and after the displacement of the buildings of block 1 , block 2 , block 3 , and block 4 .The abscissa represents the area of the Voronoi polygon of the building before the displacement; the ordinate represents the area of the Voronoi polygon of the building after the displacement.The correlation coefficients, R 2 , of each block are 0.8623, 0.9747, 0.9771, and 0.9744, indicating that the area of the Voronoi polygons before and after the building displacement is highly correlated.The quantitative analysis results show that the displacement operation can adequately maintain the spatial distribution and relative position of the building after displacement.

Conclusions
Generalizing building data at different scales plays an important role in the application of network map services, automatic drive, and urban planning.In these applications, map generalization approaches are widely used to abstract the geometric representations of building data.However, due to the reduction in map scale or the application of related generalization operators, spatial conflicts will inevitably be generated between buildings and between buildings and other features, especially in blocks with high building density.Traditional approaches mainly adopt displacement operations to resolve spatial conflicts between features.The experiments showed that when there is no redundant displacement space around the conflicting buildings, spatial conflicts cannot be resolved using only the displacement operation, except with other generalization operators.In this sense, this research has presented a heuristic spatial conflict resolution framework combining the three cartographic generalization operators of selection, displacement, and aggregation.This solution was effective for our cases because the building density in urban villages is very high, which provides a new way to solve spatial conflict between buildings and other features during map generalization.
For the core displacement operation, we have proposed a new BDGSA based on the combinatorial optimization concept.This displacement method combines the advantages of the traditional GA and SA and has a stronger spatial conflict resolution ability.We evaluated the proposed approach with high-density block data derived from Guangzhou topographic datasets.The experimental results demonstrate the feasibility of the proposed  To determine the visualization effect of the tionship was established between the building de sents a grayscale image of the building before the image of the building after the selection.It can be tribution of the corresponding regions before and

Conclusions
Generalizing building data at different scales plays an important role in the application of network map services, automatic drive, and urban planning.In these applications, map generalization approaches are widely used to abstract the geometric representations of building data.However, due to the reduction in map scale or the application of related generalization operators, spatial conflicts will inevitably be generated between buildings and between buildings and other features, especially in blocks with high building density.Traditional approaches mainly adopt displacement operations to resolve spatial conflicts between features.The experiments showed that when there is no redundant displacement space around the conflicting buildings, spatial conflicts cannot be resolved using only the displacement operation, except with other generalization operators.In this sense, this research has presented a heuristic spatial conflict resolution framework combining the three cartographic generalization operators of selection, displacement, and aggregation.This solution was effective for our cases because the building density in urban villages is very high, which provides a new way to solve spatial conflict between buildings and other features during map generalization.
For the core displacement operation, we have proposed a new BDGSA based on the combinatorial optimization concept.This displacement method combines the advantages of the traditional GA and SA and has a stronger spatial conflict resolution ability.We evaluated the proposed approach with high-density block data derived from Guangzhou topographic datasets.The experimental results demonstrate the feasibility of the proposed approach from both qualitative and quantitative aspects, indicating that this research has strong theoretical significance and practical application value.
The heuristic framework proposed in this study sequentially adopted three cartographic generalization operators.The application order of the three cartographic generalization operators was fixed.For high-density blocks in urban villages, the sequential use of three cartographic generalization operators can solve the spatial conflicts between building features and other features.However, the sequential use of the three cartographic generalization operators also has some shortcomings.In practice, it is necessary to select the optimal application order of the generalization operators according to the specific spatial conflict types and conflict features.Future work will combine relevant models in machine learning and other fields to study the optimal application order of generalization operators.At the same time, the framework only combined three cartographic operators: selection, displacement, and aggregation.However, if the density of buildings in the test area meets the requirements, it is not necessary to reduce the density through the selection operator.Therefore, it is necessary to adopt other generalization operators according to different spatial conflict types and different conflict features to resolve spatial conflicts, such as the typification operator.Follow-up research will explore the combination of more generalization operators to resolve more complex types of spatial conflicts.

ISPRS 25 Figure 1 .
Figure 1.The framework for the spatial conflict resolution of buildings in urban villages.

Figure 1 .
Figure 1.The framework for the spatial conflict resolution of buildings in urban villages.

Figure 2 .
Figure 2. Schematic diagram of building selection operation to maintain spatial distribution: (a) experimental data; (b) the result of constrained Delaunay triangulation construction; (c) the result of Voronoi-like polygon construction; (d) the result after the first round of the selection operation; (e) the result after the second round of the selection operation; (f) the result after the third round of the selection operation.

Figure 3 .
Figure 3. Representation of one chromosome and one population: (a) one chromosome, ( = 1,2, … , ) represents the displacement distance of the th building in the X and Y directions; (b) one population, ( = 1,2, … , ) represents the th chromosome in a population.

Figure 3 .
Figure 3. Representation of one chromosome and one population: (a) one chromosome, B i (i = 1, 2, . . ., n) represents the displacement distance of the ith building in the X and Y directions; (b) one population, C j (i = 1, 2, . . ., m) represents the jth chromosome in a population.
ISPRS Int.J. Geo-Inf.2023, 12, x FOR PEER REVIEW 10 of 25 corresponding to each building, perform the displacement operation on the building to obtain the result of the displacement of the building.

Figure 5 .
Figure 5.The aggregation operation of conflicting building pairs: (a) the result of not considering the area factor; (b) the result of considering the area factor.

Figure 5 .
Figure 5.The aggregation operation of conflicting building pairs: (a) the result of not considering the area factor; (b) the result of considering the area factor.

Figure 7 .
Figure 7. Result of the detected spatial conflict: (a) constructed constrained Delaunay triangulation and extracted skeleton lines; (b) detected spatial conflict area.

Figure 8 .
Figure 8.Comparison of results before and after the building selection operation: (a) building before selection operation; (b) building after selection operation.

Figure 7 .
Figure 7. Result of the detected spatial conflict: (a) constructed constrained Delaunay triangulation and extracted skeleton lines; (b) detected spatial conflict area.

25 Figure 7 .
Figure 7. Result of the detected spatial conflict: (a) constructed constrained Delaunay triangulation and extracted skeleton lines; (b) detected spatial conflict area.

Figure 8 .
Figure 8.Comparison of results before and after the building selection operation: (a) building before selection operation; (b) building after selection operation.

Figure 8 .
Figure 8.Comparison of results before and after the building selection operation: (a) building before selection operation; (b) building after selection operation.

Figure 9 .
Figure 9. Grayscale comparison of building density before and after selection: (a) grayscale image of the building before the selection; (b) grayscale image of the building after the selection.

Figure 9 .
Figure 9. Grayscale comparison of building density before and after selection: (a) grayscale image of the building before the selection; (b) grayscale image of the building after the selection.

Figure 9 .
Figure 9. Grayscale comparison of building density before and after selection: (a) grayscale image of the building before the selection; (b) grayscale image of the building after the selection.

Figure 9 .
Figure 9. Grayscale comparison of building density before and after selection: (a) grayscale image of the building before the selection; (b) grayscale image of the building after the selection.

Figure 9 .
Figure 9. Grayscale comparison of building density before and after selection: (a) grayscale image of the building before the selection; (b) grayscale image of the building after the selection.

Figure 9 .
Figure 9. Grayscale comparison of building density before and after selection: (a) grayscale image of the building before the selection; (b) grayscale image of the building after the selection.

Figure 10 .
Figure 10.The spatial conflict detection results after the building selection operation: (a) constructed constrained Delaunay triangulation and extracted skeleton lines; (b) detected spatial conflict area.

Figure 11 .
Figure 11.Experimental results of the BDGSA: (a) displacement result of the BDGSA; (b) results of the overlay with the buildings before displacement.

Figure 10 .
Figure 10.The spatial conflict detection results after the building selection operation: (a) constructed constrained Delaunay triangulation and extracted skeleton lines; (b) detected spatial conflict area.

Figure 10 .
Figure 10.The spatial conflict detection results after the building selection operation: (a) constructed constrained Delaunay triangulation and extracted skeleton lines; (b) detected spatial conflict area.

Figure 11 .
Figure 11.Experimental results of the BDGSA: (a) displacement result of the BDGSA; (b) results of the overlay with the buildings before displacement.

Figure 11 .
Figure 11.Experimental results of the BDGSA: (a) displacement result of the BDGSA; (b) results of the overlay with the buildings before displacement.

Figure 12 .
Figure 12.Remaining spatial conflict types after the displacement operation using the BDGSA.Figure 12. Remaining spatial conflict types after the displacement operation using the BDGSA.

Figure 12 .
Figure 12.Remaining spatial conflict types after the displacement operation using the BDGSA.Figure 12. Remaining spatial conflict types after the displacement operation using the BDGSA.

Figure 13 .
Figure 13.Movement effect of the conflicting building pairs before aggregation: (a) movement effect of the conflicting building pairs with tangent relationships; (b,c) movement effect of the conflicting building pairs with detached relationships.

Figure 14 .
Figure 14.Result after the aggregation operation: (a) the aggregate result of the remaining spatial conflict examples; (b) final experimental result.

Figure 13 .
Figure 13.Movement effect of the conflicting building pairs before aggregation: (a) movement effect of the conflicting building pairs with tangent relationships; (b,c) movement effect of the conflicting building pairs with detached relationships.

Figure 13 .
Figure 13.Movement effect of the conflicting building pairs before aggregation: (a) movement effect of the conflicting building pairs with tangent relationships; (b,c) movement effect of the conflicting building pairs with detached relationships.

Figure 14 .
Figure 14.Result after the aggregation operation: (a) the aggregate result of the remaining spatial conflict examples; (b) final experimental result.

Figure 14 .
Figure 14.Result after the aggregation operation: (a) the aggregate result of the remaining spatial conflict examples; (b) final experimental result.

Figure 15 .
Figure 15.Experimental results of the GA: (a) displacement results of GA; (b) results of the overlay with the buildings before displacement.

Figure 15 .Table 3 .
Figure 15.Experimental results of the GA: (a) displacement results of GA; (b) results of the overlay with the buildings before displacement.

Figure 16 .
Figure 16.Experimental results of SA: (a) displacement results of SA; (b) results of the overlay with buildings before displacement.

Figure 16 .Table 4 .
Figure 16.Experimental results of SA: (a) displacement results of SA; (b) results of the overlay with buildings before displacement.

Figure 9 .
Figure 9. Grayscale comparison of building density before and after selection: (a) grayscale image the visualization effect of the building selecti tionship was established between the building density and the g sents a grayscale image of the building before the selection; Figu image of the building after the selection.It can be seen from the tribution of the corresponding regions before and after the select 48 0.010 0.828 10.272 0.214 0.191 ISPRS Int.J. Geo-Inf.2023, 12, x FOR PEER REVIEW

Figure 9 .
Figure 9. Grayscale comparison of building density before and after sele of the building before the selection; (b) grayscale image of the building a Figure 10 shows the result of spatial conflict detection on se 10a is the constructed constraint Delaunay triangulation (i.e., gra tracted skeleton line of the gap space (i.e., solid red line).Figure spatial conflict area, which are represented in red.Compared wit tection results before the building selection operation, some spa eliminated.

Figure 9 .
Figure 9. Grayscale comparison of building density before and after sele of the building before the selection; (b) grayscale image of the building a Figure 10 shows the result of spatial conflict detection on se 10a is the constructed constraint Delaunay triangulation (i.e., gra tracted skeleton line of the gap space (i.e., solid red line).Figure spatial conflict area, which are represented in red.Compared wit tection results before the building selection operation, some spa eliminated.

Figure 9 .
Figure 9. Grayscale comparison of building density before and after sele of the building before the selection; (b) grayscale image of the building a Figure 10 shows the result of spatial conflict detection on se 10a is the constructed constraint Delaunay triangulation (i.e., gra tracted skeleton line of the gap space (i.e., solid red line).Figure spatial conflict area, which are represented in red.Compared wit tection results before the building selection operation, some spa eliminated. 48

Figure 9 .
Figure 9. Grayscale comparison of building density before and after sele of the building before the selection; (b) grayscale image of the building a Figure 10 shows the result of spatial conflict detection on se 10a is the constructed constraint Delaunay triangulation (i.e., gra tracted skeleton line of the gap space (i.e., solid red line).Figure spatial conflict area, which are represented in red.Compared wit tection results before the building selection operation, some spa eliminated.

Figure 17 .
Figure 17.Voronoi polygons of buildings before and after the displacement: (a) Voronoi polygons of buildings before the displacement; (b) Voronoi polygons of buildings after the displacement.

Figure 17 .
Figure 17.Voronoi polygons of buildings before and after the displacement: (a) Voronoi polygons of buildings before the displacement; (b) Voronoi polygons of buildings after the displacement.

Figure 18 .
Figure 18.Relationship between the areas of the Voronoi polygons of buildings before and after displacement: (a) relationship between the areas of the Voronoi polygons of buildings before and after displacement of block ①; (b) relationship between the areas of the Voronoi polygons of buildings before and after displacement of block ②; (c) relationship between the areas of the Voronoi polygons of buildings before and after displacement of block ③; (d) relationship between the areas of the Voronoi polygons of buildings before and after displacement of block ④.

Figure 18 .
Figure 18.Relationship between the areas of the Voronoi polygons of buildings before and after displacement: (a) relationship between the areas of the Voronoi polygons of buildings before and after displacement of block

Figure 9 .
Figure 9. Grayscale comparison of building density bef of the building before the selection; (b) grayscale image

Table 1 .
Statistical information before and after the selection operation of each block.

Table 1 .
Statistical information before and after the selection operation of each block.

Table 1 .
Statistical information before and after the selection operation of each block.

Table 1 .
Statistical information before and after the selection operation of each block.

Table 1 .
Statistical information before and after the selection operation of each block.

Table 1 .
Statistical information before and after the selection operation of each block.

Table 2 .
Statistical information of the BDGSA displacement results.

Table 1 .
Statistical information before and after the selection operation of each block.

Table 1 .
Statistical information before and after the selection operation of each block.

Table 1 .
Statistical information before and after the selection operation of each block.

Table 1 .
Statistical information before and after the selection operation of each block.

Table 3 .
Statistical information of the GA displacement results.

Table 1 .
Statistical information before and after the selection operation of each block.

Table 1 .
Statistical information before and after the selection operation of each block.

Table 1 .
Statistical information before and after the selection operation of each block.

Table 1 .
Statistical information before and after the selection operation of each block.

Table 4 .
Statistical information of the SA displacement results.

Table 1 .
Statistical information before and after the selection operation of each block.

Table 1 .
Statistical information before and after the selection operation of each block.

Table 1 .
Statistical information before and after the selection operation of each block.

Table 1 .
Statistical information before and after the selection operation of each block.

Table 5 .
Statistical information of the BDGSA displacement distance.

Table 1 .
Statistical information before and after the selection operation o

Table 1 .
Statistical information before and after the selection operation o

Table 1 .
Statistical information before and after the selection operation o

Table 1 .
Statistical information before and after the selection operation o

Table 6 .
Statistical information of the GA displacement distance.

Table 1 .
Statistical information before and after the selection operation o

Table 1 .
Statistical information before and after the selection operation o

Table 1 .
Statistical information before and after the selection operation o

Table 1 .
Statistical information before and after the selection operation o

Table 7 .
Statistical information of the SA displacement distance.

Table 1 .
Statistical information before and after the selection operation o

Table 1 .
Statistical information before and after the selection operation o

Table 1 .
Statistical information before and after the selection operation o

Table 1 .
Statistical information before and after the selection operation o

Table 7 .
Statistical information of the SA displacement distance.

Table 1 .
Statistical information before and after the sele