Research on the Evolution Mechanism of Congestion in the Entrances and Exits of Parking Facilities Based on the Improved Spatial Autoregressive Model

,e entrance and exit area of parking facilities has the characteristics of high concentration of urban traffic and prominent traffic intertwining phenomenon, which easily induces rapid congestion of mixed heterogeneous traffic at specific times and local locations and quickly spreads to the entire road section or even a larger area. In order to better understand the congestion distribution characteristics and propagation effects of access section of the parking entrance and exit from the mid and microperspective, a 5m ∗ lane width pixel grid is used to divide the frontage road research. It also proposes a spatially robust autoregressive model and complex network tools suitable for analysis of local traffic flow to analyze it. ,e results show that as spatial scale increases, the congestion propagation decreases sharply and spatial adjacency within the fourth order can account for more than 90% of the propagation; the frontage road to the entrance and exit is the place where the congestion first happens, and the congestion gradually attenuates as it propagates to the inner lane and the upstream of the road segments; the lateral congestion propagation attenuates faster, so the area affected by congestion is mainly distributed in the outermost lane. ,is paper can provide theoretical guidance for alleviating traffic congestion in the entrance and exit areas of parking facilities and has theoretical and empirical significance.


Introduction
Increasing car park [1], which exacerbates the contradiction between the supply and demand of urban parking, along with the lag of parking facilities planning [2][3][4] and other factors, causes more severe traffic congestion problems in the entrances and exits of parking facilities. Since the entrances and exits of parking facilities are generally set up on secondary arterial roads or branch roads with low capacity, the interweaving of incoming and outgoing vehicles with external traffic will reduce the traffic capacity [5,6], increase the delay time of main road traffic [7,8], and even cause a bottleneck and congestion on the upstream road.
At present, research focused on the traffic conditions in entrances and exits of the off-street parking is still very rare. Previous research mainly focuses on the complex impact of the bus stop [9][10][11][12][13][14] and on-street parking [15][16][17][18] on the traffic flow. Compared with other road facilities, the settings of entrances and exits of parking facilities can bring more complicated impact: on the one hand, due to the uncertainty of the occurrence and attraction of parking facilities, the temporal and spatial distribution of parked vehicles are more disorderly; the intermittently concentrated entering and departing of vehicles enhance the occasional disturbance to the road; on the other hand, the interwoven flow of incoming and outgoing vehicles, the traffic flow of the main road, and the flow of motor vehicles and nonmotorized vehicles [19] increase the heterogeneity of the traffic flow and the complexity of the occurrence of congestion in this area to a certain extent.
In the analysis of traffic congestion, scholars have built various models to reveal the essence of congestion. At the macrolevel [20][21][22][23][24], traffic congestion is generally regarded as a process of compression, blockage, and spread of traffic flow. At the microlevel [25][26][27][28][29], more attention is paid to conflicts, queuing, and easing of vehicles. Spatial autocorrelation is a geographic concept, which refers to the potential interdependence between observed data of some variables in the same distribution area [30,31]. e development of this theory in the field of transportation has resulted in its wide use in traffic feature recognition [32] at the macrolevel and mid level and accident analysis [33][34][35][36]. In addition, complex network theory, as a new scientific theory for revealing network complexity phenomena, has a wide range of applications for the organization and optimization of transportation logistics [37][38][39][40] and the macrointerpretation of transportation systems [41][42][43][44][45][46], but it is rarely used in the study of traffic flow characteristics at the meso-microlevels.
To sum up, although there are rich research outcomes on congestion propagation, research on the frontage roads of parking facilities is still relatively rare; besides, although the spatial autocorrelation theory has developed in the field of transportation, due to its higher requirements on spatial stability, it is not suitable for the research of transportation theory at the mid and microlevel. In response to the abovementioned problems, improving the spatial autocorrelation model, building a model of the spatial congestion propagation in the entrances and exits of the parking facilities from the mid and microperspective under the guidance of the complex network theory and other theories, and exploring the attenuation effect and key nodes in the process of congestion propagation based on the observed data are of great significance for alleviating traffic congestion in the entrances and exits of parking facilities.

Data Description
In order to determine the objective of the study, the frontage road of this article is defined as follows: the free road segments from about 140 m upstream to about 20 m downstream from the entrance and exit access points. is paper selects the entrance and exit of the parking lot of the West Affiliated Hospital of Yangzhou University as a typical off-street parking area. e frontage road is the one-way three-lane road (a part of the two-way six-lane road with a separation zone in the middle), and the upstream area is far away from the intersection. e traffic volume and entry rate vary significantly in different periods, and the road segment is free of vegetation, billboards, and other structures, which are conducive for data collection. e data set in this paper is obtained by aerial filming, with a wind-proof drone MAVIC_AIR2 recording the traffic conditions in the experimental area at an altitude of 100-150 m (each video lasts 10-15 minutes, 100 videos in total). en, the open-source vehicle-tracking algorithm based on YOLOV5 and Deep-Sort [47] is used to monitor the location parameters of the moving vehicle [48,49] in real time (as shown in Figure 1).
In order to describe the location of the lanes in an easy way, the lanes in the experimental area are numbered, where the outermost lane is Lane 3, the middle lane is Lane 2, and the innermost lane is Lane 1. In addition, in order to study the relationship of the traffic state of each spatial unit in the entrance and exit and the spatial congestion propagation, this paper divides the road into three lanes and then further divides these three lanes into 5-meter wide spatial units, which can ensure the spatial accuracy of the model and meet the requirements of velocity measurement, as shown in Figure 2. e investigation finds that the time for vehicles to enter the parking lot is generally within the range of 3 to 13 seconds. In order to extract the traffic information of each spatial unit to the greatest extent, this paper uses a time interval of 3 seconds to calculate the average speed of vehicles in each spatial unit. e sample data are shown in Table 1.

Improved Spatial Autoregressive Model.
e traditional spatial autoregressive model requires the data to have spatial stability. However, due to the influence of external factors such as the actual environment and traffic conditions, the spatial stability of the measured data is often poor. In order to reduce the influence of external factors on the modeling results, this paper improves the traditional spatial autoregressive model and builds a robust spatial autoregressive model based on repeated observation data that consider spatial instability and multiple correlations.
e traditional spatial autoregressive model adds a spatial lag term to the general regression model, that is, the dependent variable of a spatial object is related to the independent variable on the same object, and it is also related to the independent variable and dependent variable of adjacent objects. Its general form is shown in Here, z represents the spatial attributes of the research object, or the dependent variable; X represents the spatial attribute of the adjacent objects, or the independent variable; ρ is the spatial autoregressive coefficient; β is the regression coefficient; μ represents the residual of the model; λ represents the residual spatial regression coefficient; ε represents the random error; and W 1 and W 2 represent the spatial adjacency matrix of the spatial units of frontage road.
Observing the influencing factors of the attribute value of the spatial unit in the study area, it is not difficult to find that the change of the attribute value of the spatial unit is mainly affected by the adjacent spatial unit, while the explanatory variable X has almost no contribution. erefore, X and W 2 in formula (1) can be assigned to 0, and the firstorder spatial autoregressive model can be obtained, as shown in In order to study the propagation mechanism of congestion at different spatial scales, this paper divides adjacent spatial units according to "distance" and takes the influence of spatial units under each adjacent number on the research object as an independent variable. Taking into account the correlation between the independent variables, the direct use of traditional parameter estimation methods based on the least squares method will cause larger errors. erefore, this paper proposes an improved first-order spatial autoregressive model that takes into account spatial instability and multiple correlations.
In the improved spatial autoregressive model, the disturbance of the research object is the dependent variable and the disturbances of other spatial units are the independent variables to describe the relationship between the spatial units in the research attributes, as shown in formula Here, y ′ represents the velocity disturbance sequence of the research object; W k represents the k-order spatial adjacency matrix considering spatial instability; ε represents the error term of the improved spatial autoregressive model; l is the constant term; ρ k represents the spatial lag term; the regression coefficient of W k ; c represents the regression coefficient of the constant term; k represents the number of spatial adjacencies; and K is the maximum number of spatial adjacencies considered. For the space unit under study, the historical sequence of the previous c cycle formula (4) can be constructed to obtain the velocity sequence of each space unit: Y � y 1,1 y 2,1 y 1,2 y 2,2 · · · y n−1,1 y n,1 y n−1,2 y n,2 ⋮ ⋱ ⋮ y 1,c y 2,c · · · y n−1,c y n,c Here, y i represents the data of space unit i in c cycles and y i � [y i,1 , y i,2 , y i,3 , . . . , y i,c ] T is the speed history data sequence of space unit i.   erefore, the speed perturbation sequence of the unit of the space is calculated according to the speed history data sequence. As shown in formulas (5) and (6), y i is the historical velocity sequence of space unit i and y is the historical velocity perturbation sequence:

Considering the Spatial Adjacency Matrix at Is Not
Stationary. e traditional spatial adjacency matrix based on the boundary adjacency method or the center-of-gravity distance method assumes that all adjacent objects have the same weight value, which is inconsistent with the actual situation of congestion propagation. On the one hand, the number of adjacencies in different spatial units is different. On the other hand, the traffic flow is not homogenous and the effects of adjacent units in different directions are not the same. e spatial adjacency matrix determined based on the boundary adjacency method is shown in formulas (5) and (6). In the formula, w k i,j indicates that there is a k-order adjacency relationship between space unit i and space unit j. Determining the spatial adjacency matrix based on the boundary adjacency method is a commonly used method to determine the adjacency relationship between area units; if two spatial units have the same boundary, they have a firstorder spatial adjacency relationship and the k-order adjacency relationship can be determined according to the transitivity of the spatial adjacency relationship.

(7)
In order to accurately describe the relationship between the spatial units, first, according to the spatial adjacency matrix W k , determine the set of spatial units that have a k-order adjacency relationship with each spatial unit, denoted as f k (i, j), j � 1, 2, 3, . . . , r, . . . , R; then, the rth space unit that has an adjacency relationship of order k with the space unit i is denoted as f k (i, r).
Let vector y i � [y i,1 , y i,2 , y i,3 , . . . , y i,r , . . . , y i,R ] T be the velocity history perturbation sequence of space unit i, and H k i is the matrix of velocity history perturbation sequences of all space units that have a k-order adjacency relationship with space unit i, as shown in In order to preserve the features in H k i to the greatest extent and reduce their dimensionality, the principal component analysis method is used to extract the first principal component that is most similar to the velocity history disturbance sequence ] T be the weight value of each spatial unit that has a kth order adjacency with the spatial unit i, in order to ensure each weight value can accurately describe the degree of association between spatial units, and the covariance between H k i u k i and y i should be the largest, as shown in Use the Lagrangian method to solve the model, as in formula (9): Calculate the partial derivative of F with u k i and λ, as shown in According to formula (12), e value u k f k (i,r) in u k i represents the spatial association relationship of the rth spatial unit that has a k-order adjacency with the spatial unit i; the sign of u k f k (i,r) indicates the positive or negative correlation between the space unit i and the rth space unit that has a k-order adjacency relationship, and the absolute value indicates the relationship between the space units. e greater the correlation, the greater the absolute value, and vice versa. In actual road conditions, the propagation of congestion is the same direction, so only the positive numbers in the vector are retained to construct the smoothed spatial adjacency matrix W k , as shown in 4 Journal of Advanced Transportation · · · · · · · · · · · · · · · · · · w k n,1 w k n,2 · · · w k n,j · · · w k n,n

Parameter Calculation Method.
In order to study the spatial propagation process of traffic congestion in the entrance and exit area of parking facilities, in the improved spatial autoregressive model, the disturbance of the research unit is used as the dependent variable and the spatial lag term W k y ′ of the disturbance under different spatial adjacency k is used as an independent variable; there is a linear correlation between the spatial lag terms W k y ′ under different spatial adjacency numbers k.
According to the space lag term W 1 y ′ of the disturbance under different number of space delay periods, the independent variable matrix is constructed, as shown in

Improved Spatial Autoregressive Model Parameter
Estimation Process. For the convenience of calculation, the dependent variable matrix y ′ of the improved spatial autoregressive model is marked as F 0 , the independent variable matrix is marked as E 0 , and the parameter estimation is performed. Extract the first principal component t 1 from the independent variable matrix, where t 1 is the linear combination of E 0 of the independent variables as shown in formulas (17) and (18). Since the dependent variable contains only one variable, the first principal component of the dependent variable is F 0 .
For the needs of regression analysis, the extracted principal components must simultaneously meet the requirements of the first principal components t 1 and F 0 of the variables to reflect the characteristics of the variable group to the greatest extent and the correlation between the two first principal components to reach the maximum. Considering the requirements of modeling comprehensively, covariance is selected to characterize the principal components, as shown in According to the Lagrangian multiplier method, the following formula can be obtained: Respectively, find the partial derivatives of u 1 and λ to get the following formulas: According to formula (22), we can get According to the principle of principal component analysis, the sample is expressed as a linear combination of components, that is, regression on t 1 is performed on E 0 and F 0 , respectively, as shown in formulas (24) and (25), where E 1 and F 1 are residual matrixes.
According to (25), the following formulas can be obtained from the least square estimation: Replace E 0 and F 0 with the residual terms E 1 and F 1 . Repeat the above steps to obtain the weight vector u 2 (formula (28)) of the second principal component t 2 of E 0 and the regression coefficients p 2 (formula (29)) and r of E 1 and F 1 with respect to r 2 (formula (30)).
Journal of Advanced Transportation 5 Assuming that a total of h principal components are extracted, the regression model between F 0 and E 0 can be obtained, as shown in F 0 � e 1 r 1 + e 2 r 2 + e 3 r 3 + · · · + e H r H + F H , where B 0 is the perturbation sequence F 0 of the dependent variable y ′ of the improved spatial autoregressive model and the regression coefficient vector between the perturbation matrix E 0 of the independent variable. Determine the value of the regression coefficient according to the position of the respective variable in Z (K). For example, the kth element in B 0 represents the estimated value of ρ k .

Determination of the Optimal Number of Principal
Components. e partial least squares method is a method of establishing a regression model by extracting the principal components of the independent variables. If the number of principal components is too small, the relationship between the variables cannot be reflected; otherwise, overfitting will occur. In order to ensure the validity of the model, the crossvalidity criterion is selected to determine the number of principal components to be extracted. e cross-validity criterion is mainly to test whether the extracted hth principal component can significantly improve the prediction performance of the model compared to the (h − 1) th principal component. For the improved spatial autoregressive model, the cross-validity criterion is selected to determine the number of principal components, which is mainly calculated by using all the spatial unit disturbance value sequences to extract the difference between the predicted value and the true value of the spatial unit disturbance when the number of principal components is (h − 1). e sum of squared deviations SS h−1 is shown in formula (34). And, extract the hth principal component based on the sequence after removing the perturbation value of its own spatial unit and calculate the sum of squared deviations of the predicted value and the true value based on the new sequence PRESS h , as shown in According to SS h and PRESS h , calculate the effectiveness of extracting the hth principal component Q 2 h (formula (36)). When Q 2 h ≥ 1 − 0.95 2 , it indicates that the hth principal component has significantly improved the prediction performance of the model, and it is necessary to continue to extract the principal components; when Q 2 h < 1 − 0.95 2 , indicating the introduction after the hth principal component, the model has no obvious improvement, and there is no need to continue to extract principal components.

Improving the Validation of the Spatial Autoregressive
Model. For the speed disturbance sequence studied, the improved spatial autoregressive model is used to fit and the model's fitting effect is tested from two aspects: the model's goodness-of-fit (R 2 ) index and visual observation. When R 2 is closer to 1, it means that the fitting effect is good. On the contrary, when it is close to 0, it means that the fitting effect is poor. When the change trend of the measured value curve and the predicted curve is close, it indicates that the model can reasonably describe the space average speed at each position. On the contrary, it means that the model cannot reasonably reflect the real situation.

Analysis Method of Congestion Evolution in the Entrance and Exit
Area. e improved spatial autoregressive model takes the velocity disturbance of each spatial unit as the dependent variable and independent variable and uses the partial least square method to estimate the parameters of the improved spatial autoregressive model and obtain the parameter estimated value ρ k under the number of adjacent spaces in each space. e estimated value ρ k represents the degree of correlation between the degree of velocity disturbance of the k-order adjacent unit of each spatial unit and the degree of correlation between the spatial unit. erefore, the parameters in the improved spatial autoregressive model are defined as the propagation structure of each spatial unit velocity disturbance in space. e influence of the disturbance of the space unit on its adjacent space unit of order k is denoted as V(W k ), as shown in V W k � ρ k W k (k � 1, 2, . . . , k, . . . , K).
e spatial propagation process of congestion can be summarized as a process in which the influence of local spatial units gradually spreads outwards through neighboring units. Congestion first spreads from spontaneously congested spatial units to neighboring spatial units and then from neighboring spatial units to second-order units. e continuous spread to the surroundings eventually leads to the proliferation and transformation of congestion. erefore, the spatial propagation effects of local congestion can be calculated by formula (37), and the distribution of congestion spatial propagation effects under different spatial adjacency numbers can be investigated and then the spatial process of congestion propagation can be studied.

Propagation Effect of Congestion.
Suppose that the space unit is congested and the speed of its first-order adjacent unit is disturbed through the adjacency relationship. e first-order propagation effect generated by the space 6 Journal of Advanced Transportation unit, that is, the influence on each adjacent unit can be represented by V(W 1 ) i,j (j is the first-order adjacent space unit of i). Similarly, after congestion is propagated k times, the influence of spatial unit i on its k-order adjacent unit is expressed as V(W k ) i,j . en, the ratio of the kth propagation effect to the total propagation effect is shown in

Interaction Degree Matrix between Spatial
Units. e change of velocity perturbation of each space unit in the study area is the result of the combined effect of multiple space units. In order to quantify the degree of mutual influence between spatial units, the disturbance formed by the influence of spatial unit j by i is recorded as v i,j and then the degree of influence of spatial unit i on spatial unit j is calculated: W � w 1,1 w 1,2 · · · w 2,1 w 2,2 · · · w i,1 w i,2 · · · w 1,j · · · w 1,n w 2,j · · · w 2,n w i,j · · · w i,n ⋮ ⋮ ⋱ w n,1 w n,2 · · · ⋮ ⋱ ⋮ w n,j · · · w n,n

Identification of Key Nodes.
In the discipline of complex networks, nodes that are easily damaged and have a greater impact on other nodes are called key nodes. With reference to the definition of key nodes in a complex network, this paper proposes a method for identifying key nodes in the congestion evolution process of the entrance and exit areas: the space where the speed of the space unit is susceptible to influence and has a greater impact on the speed of adjacent space units. e unit is defined as the key node in the congestion evolution process of the entrance and exit area.  Table 2. e first of the above five traffic states is the congestion formed under the influence of its own factors, while the other four are under the comprehensive influence of the congestion propagation effect of neighboring units. When the disturbance rises sharply and reaches a certain threshold, the traffic state gradually changes from unblocked to congested, indicating that the propagation effect is the main reason for the decline in service level and even congestion in the entrance and exit areas.
Based on the method of judging the traffic state, this paper defines the space unit with a speed greater than 10 km/h and an impact on adjacent space units greater than 30% as a key node in the congested traffic state/process.

Goodness-of-Fit Test.
e goodness-of-fit test is conducted on the velocity disturbance sequence in the experimental area under the improved spatial autoregressive model. According to the parameter fitting results, the goodness-of-fit (R 2 ) value of the improved spatial autoregressive model is 0.7423, indicating that the model has a good fitting effect.
According to the test shown in Figure 3, divergence between the actual observed value and predicted value in each spatial unit is small and the actual observed value curve and the predicted value curve are very close, indicating the improved spatial autoregressive model can reasonably describe the trend of speed at various locations. In a nutshell, the model performs well in predicting the speed disturbance of each spatial unit and can truly and reasonably reflect the spatial correlation of the traffic of various units in the area.

Analysis on the Evolution Law of Congestion in the Entrance and Exit.
Congestion evolution analysis, the method proposed in this paper, is used to study the evolution of congestion in the entrance and exit. Twenty sets of observed data whose traffic flow is 2000 pcu/h and entry rate is 20% are collected for research. By calculating the congestion propagation in the entrance and exit and the influence matrix of different spatial units, the key nodes of the evolution of congestion in the area can be perceived and its formation mechanism can be analyzed.

Analysis of Congestion Propagation in the Entrance and
Exit. Based on the criteria in Table 2, the traffic of each spatial unit in the entrance and exit can be determined and the spatial distribution is shown in In order to further explore the interaction between spatial units and understand the characteristics of congestion propagation, five spatial units with different traffic are selected and the distributions of congestion of spatial units of different orders are compared. e result is shown in Figure 5. It can be seen from the figure that the propagation whose spatial adjacency is within the fourth order accounts for more than 90% of the total propagation. Propagation gradually decreases as the spatial adjacency increases. is shows that when congestion happens in a spatial unit, its impact on the most adjacent spatial unit (first-order adjacent unit) is the greatest, and as the spatial scale increases, the impact decreases sharply.

Spatial Evolution Mechanism of Congestion in the
Entrance and Exit. In order to further analyze the evolution mechanism of congestion in the entrance and exit, this paper calculates the influence matrix of spatial units based on the improved spatial autoregressive model and quantifies the degree of mutual influence between spatial units under the influence of congestion propagation. e spatial unit is regarded as the node, the mutual influence between the spatial units as the edge, and the degree of mutual influence as the weight, and a directed weighted network of the spatial units can be constructed, as shown in Figure 6.
It can be seen from Figure 6 that the directed weighted network of spatial units in the experimental area takes the shape of a "band," which is consistent with the abovementioned attenuation characteristics of congestion propagation. To better analyze and study the evolution mechanism of congestion in the area, this paper divides the congestion evolution process into entrance and exit section area, transition area, and upstream area based on different spatial locations where the congestion happens. Figure 7 draws a directed weighted network for the three areas and uses the abstract network to analyze the evolution process of congestion.
(1) Section Area of the Entrance and Exit. e directed weighted network near the entrance and exit shows that affected by the parking behavior, there are significant differences in the weight values of the spatial units near the section. Based on the directed weighted network near the section (Figure 8), Figure 9 is drawn. All units in Figure 9 are spatial units in the section area of the entrance and exit; A, B, and C indicate Lane 1, Lane 2, and Lane 3, respectively, 1, 2, 3, 4, and 5 indicate adjacent spatial units from upstream to downstream of the road; C3, C4, and C5 are spontaneous congestion units, and other units are general units. It can be seen from the figure that the lateral interference on the spatial unit accounts for about 85% of the total in Lane 1 (outermost lane), and the lateral interference on the spatial unit accounts for 60% of the total in Lane 2 (middle lane). e main reason for the spontaneous congestion units in Lane 3 (innermost lane) is that the service capacity of the entrance and exit of the parking lot cannot meet the needs of parking and the entry time is prolonged due to the conflict between motor vehicles and nonmotorized vehicles. ese ultimately make it impossible to ease the vehicles and cause congestion to propagate upstream.
(2) Transition Area. From the directed weighted network in the transition area ( Figure 10) and the diagram of spatial units ( Figure 11) (all the units in the figure are units in the transition area and all are general units, and the meanings of all the other symbols are the same as the preceding text), it can be seen that under the influence of the vehicle stagnation, the ratio of longitudinal interference to the total interference the spatial units receive in Lane 3 varies within the range of [0.5, 0.65], and the ratio of longitudinal interference to the total interference and the distance between the spatial units and the entrance and exit is inversely correlated. e ratio of the longitudinal interference to the total interference the spatial units receive on Lane 2 varies within the range of [0.35, 0.6], and the ratio of the lateral interference to the total interference the spatial units receive on Lane 2 varies within the range of [0.1, 0.45]. Moreover, the ratio of longitudinal interference to the total interference is positively correlated with the distance between the spatial units and the entrance and exit, and the ratio of lateral interference to the total interference is negatively correlated   with the distance between the spatial units and the entrance and exit. is is because the spatial units in Lane 2 is affected by the combination of stagnation of vehicles in the same lane and lane changing of driving vehicles in Lane 3, and the speed is disturbed. As the distance increases, the impact of lane changing on Lane 2 gradually decreases, and the impact of stagnation on the upstream units gradually becomes the main factor, which is manifested by the gradual increase in the weight value of the units in Lane 2. e longitudinal interference the spatial units receive on Lane 1 accounts for about 65% of the total, which is mainly caused by the vehicles changing lanes.
(3) Upstream Area. e directed weighted network of the spatial units in the upstream area of the entrance and exit is shown in Figure 12. e weight value between the spatial units in the upstream area tends to be stable, which means that the congestion is propagated upstream in the same pattern. e relationship between spatial units is illustrated by Figure 13 (all the units in the figure are units in the upstream area and all are general units, and the meanings of all the other symbols are the same as the preceding text). It can be seen from the figure that the lateral interference the spatial units receive in Lane 1 accounts for more than 90% of the total. e main reason for this phenomenon is that the vehicles in Lane 1 rarely stagnate, and the spatial units are mainly affected by "friction effect"; in Lane 2, the lateral interference the spatial units receive in Lane 2 is about 60%.
is is because the intermittent queuing of vehicles parked near the entrance and exit causes the vehicles in the upstream area to stagnate. erefore, the mutual influence between the spatial units in the same lane is relatively high. Lane 3 is the lane directly connected to the entrance and exit. Due to the combination of frequent queuing, following the traffic flow and lane changing of parked vehicles, the longitudinal and lateral interferences the spatial units receive in Lane 3 are both about 50%.
To sum up, from the perspective of the impact value, the congestion propagation is the strongest in the frontage road to entrance and exit of the parking lot, where the congestion first happens; there is moderate propagation in the transition area, the main area where congested vehicles stagnate; congestion propagation is weak in the downstream area, stagnation of vehicles is barely seen, and there is a relatively fast-flowing traffic. From the perspective of the proportion, from the frontage road to the transition area and then to the upstream area, the proportion of longitudinal interference gradually increases, while the proportion of lateral interference gradually decreases.
is shows that the lateral propagation of congestion attenuates faster than longitudinal propagation.

Key Congestion Nodes.
In order to further analyze the relationship between the location and number of key nodes in the process of congestion evolution and the traffic volume of the road segments and the entry rate, this paper calculates the location of key nodes of congestion propagation when the traffic volume is 1500 pcu/h, 2000 pcu/h, and 2500 pcu/h, respectively, and when the entry rate of parked vehicles stands at 15%, 20%, and 25%, respectively (as shown in Table 3). Obviously, the traffic volume and the entry rate are positively correlated with the number of key nodes. As the traffic volume and the entry rate increase, key nodes spread from the outer lane to the inner lane and from the frontage road to the entrance and exit of the parking lot to the upstream area. Further comparative analysis shows that when the traffic volume is 1500 pcu/h, the change of the entry rate does not have a significant impact on the number of key nodes; when the traffic volume increases to 2000 pcu/ h (especially when it increases to 2500 pcu/h), the number of  Immune unit   65 66 67 68  44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60  62 63 64  61  35 36 37 38 39 40 41 42 43   78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94  96 97 98  95  99 100 101 102  69 70 71 72 73 74    key nodes increases significantly along with the increase in the entry rate, which shows that the impact of the entry rate on the key nodes is enhanced as the traffic volume of the road segments increases. e analysis suggests that the degree, distribution, and number of congestion of spatial units in the entrance and exit are constantly changing. When the traffic volume of the road segments and entry rate are high, the degree and number of congestion of the spatial units of the entrance and exit will increase significantly, and vice versa. However, for each specific unit, whether it will become a congestion unit is not determined solely by the traffic volume or the entry rate, but by the combination of the traffic volume, the entry rate, and the location of the spatial unit. In the directed weighted network in the entrance and exit, some spatial units (key nodes under the effect of congestion propagation) quickly evolve into congestion units because of greater turbulence caused by the increase in traffic volume or entry rate. However, some units are less or not susceptible to the increase in traffic volume or entry rate. erefore, in daily traffic management, the key nodes in the evolution process should be eased; for spontaneous congestion units, by improving the service capacity of entrances and exits and reducing the entry time of parked vehicles, the impact of spontaneous congestion units on adjacent units can be reduced; in terms of the key nodes under the influence of propagation, reducing stagnation of vehicles can alleviate the longitudinal propagation of congestion; in addition, guiding vehicles to enter the parking lot from other entrances with facilities such as parking grading guidance screens to reduce entry rate improve the frontage roads to the entrances and exits and enhance the quality of travel for citizens.

Conclusions and Discussion
is paper starts with studying the spatial evolution pattern of congestion propagation in the entrance and exit of a parking lot. First, it discusses the shortcomings of traditional spatial autoregressive model and proposes an improved robust spatial autoregressive model based on the actual characteristics of the traffic flow in the entrance and exit. Second, it constructs a smooth spatial adjacency matrix based on the principal component analysis and designs a parameter estimation method based on partial least squares.
ird, it selects the frontage road to entrance and exit of a typical parking lot for research, and the fitting result (R 2 is 0.7423) shows that the model has a good fitting performance for the velocity disturbance sequence of the spatial units. Fourth, by identifying and observing the distribution of different types of congestion units in the research section, the mechanism of congestion formation and the main propagation direction are analyzed. Finally, it explores and analyzes the congestion propagation and the spatial evolution mechanism of the congestion in the experimental area. e results show that as the spatial scale increases, the congestion propagation decreases sharply, and spatial adjacency within the fourth order can account for more than 90% of the propagation; the frontage road to the entrance and exit is the place where the congestion first happens; the congestion gradually attenuates as it propagates to the inner lane and the upstream of the road segments; the lateral congestion propagation attenuates faster, so the area affected by congestion is mainly distributed in the outermost lane. In summary, the proposed method for analyzing the spatial propagation characteristics of traffic congestion can be better applied to the frontage roads of parking facilities. e results of the case analysis show that for the local traffic congestion at the entrance and exit, the relevant departments should proceed from a holistic perspective to accurately monitor and orderly channel the key space nodes. e conclusions are of great theoretical and practical significance for improving the frontage roads to the entrances and exits of parking facilities and the quality of travel for citizens.

Data Availability
e measured data used to support the findings of this study are included within the article and the supplementary files.