An Urban Built Environment Analysis Approach for Street View Images Based on Graph Convolutional Neural Networks

Abstract

Traditionally, research in the field of traffic safety has predominantly focused on two key areas—the identification of traffic black spots and the analysis of accident causation. However, such research heavily relies on historical accident records obtained from the traffic management department, which often suffer from missing or incomplete information. Moreover, these records typically offer limited insight into the various attributes associated with accidents, thereby posing challenges to comprehensive analyses. Furthermore, the collection and management of such data incur substantial costs. Consequently, there is a pressing need to explore how the features of the urban built environment can effectively facilitate the accurate identification and analysis of traffic black spots, enabling the formulation of effective management strategies to support urban development. In this study, we research the Kowloon Peninsula in Hong Kong, with a specific focus on road intersections as the fundamental unit of our analysis. We propose leveraging street view images as a valuable source of data, enabling us to depict the urban built environment comprehensively. Through the utilization of models such as random forest approaches, we conduct research on traffic black spot identification, attaining an impressive accuracy rate of 87%. To account for the impact of the built environment surrounding adjacent road intersections on traffic black spot identification outcomes, we adopt a node-based approach, treating road intersections as nodes and establishing spatial relationships between them as edges. The features characterizing the built environment at these road intersections serve as node attributes, facilitating the construction of a graph structure representation. By employing a graph-based convolutional neural network, we enhance the traffic black spot identification methodology, resulting in an improved accuracy rate of 90%. Furthermore, based on the distinctive attributes of the urban built environment, we analyze the underlying causes of traffic black spots. Our findings highlight the significant influence of buildings, sky conditions, green spaces, and billboards on the formation of traffic black spots. Remarkably, we observe a clear negative correlation between buildings, sky conditions, and green spaces, while billboards and human presence exhibit a distinct positive correlation.

1. Introduction

With the rapid development of cities, people’s living standards gradually improve; the number of motor vehicles increases; and the problem of the unbalanced development of people, vehicles, and roads becomes more prominent, which can easily produce traffic hazards and lead to traffic accidents [1,2,3]. Traffic black spots are defined as locations with a high incidence of traffic accidents, which account for only 0.25% of the total road network mileage, although the number of traffic accidents occurring in them accounts for 25% of the total number of accidents [4,5,6,7]. Therefore, accurately identifying traffic black spots can reduce the occurrence of traffic accidents.

How to accurately identify traffic black spots on the basis of existing data has not yet been uniformly determined. In previous studies, traffic accident data have mostly come from the historical accident records of the traffic management department, and the number of traffic accidents is analyzed to determine the locations of the traffic black spots by using mathematical and statistical methods [8,9,10]. The data are relatively singular and not easy to collect and manage, and more importantly the number of accident features is small, which does not allow for a specific analysis of the causes of traffic black spots. The emergence of geographic big data and the rapid development of machine learning and deep learning make it possible to portray the built environment of urban roads from street view image data [11,12,13]. In real life, traffic accidents are also more likely to occur on densely populated streets with crowded and dilapidated buildings, without better road infrastructure, and with missing road markings and blurred and uneven road surfaces [14,15,16]. These features are all part of the built environment of urban roads, which can be described by street image data. Using the existing deep learning semantic segmentation technology, street image data can be segmented into more than 100 urban built environment features [13], which could be very useful as a supplement to a traffic management department’s historical accident record data feature factors, so as to carry out a specific analysis of the formation of the traffic black spots.

Traffic black spot identification is the focus of research in the field of traffic safety. By accurately identifying traffic black spots in urban roads, urban traffic management departments can make targeted improvements to reduce the incidence of traffic accidents to protect people’s right to life, property, and health. In this paper, on the basis of the above research, street view image data are assessed using deep learning semantic segmentation technology to describe the characteristics of the built environment of urban roads, so as to obtain the information containing the characteristics of the built environment of urban roads, which can be used to identify black spots in urban road traffic. At the same time, the data from traffic accidents and street view image data are processed and constructed into the structure of a graph network and a graph convolutional neural network (GCNN) model is used to identify black spots in traffic. This GCNN model is used for the recognition of traffic black spots to make the recognition results more accurate. Firstly, taking the road intersection as the research unit, the street view image data are extracted using a semantic segmentation model to portray the built environment of urban road intersections and machine learning models, such as support vector machine (SVM), gradient boosting (GBDT), and random forest (RF) models, which are used to identify the traffic black spots. Then, the traffic black spot recognition method is improved by considering that the interconnection of built environment features of adjacent road intersections will have an impact on the formation of traffic black spots. The traffic accident data and street scene image data are processed, the graph network structure is constructed, and the graph convolutional neural network (GCNN) model is used for the recognition of traffic black spots. In the analysis of the causes of traffic black spots, the importance of each feature of the built environment of urban road intersections on the formation of traffic black spots and the interconnections between each feature are analyzed to make the research results more reliable. The proposed method can effectively improve the urban road planning and design process and make human environment planning and design more scientific, more reasonable, and more in line with the standard of the ideal modern city. In addition, the street scene image data are easier to collect and manage, which can save a significant amount of manpower cost and time cost.

In this paper, we take the Kowloon Peninsula of Hong Kong as the study area and road intersections as the study unit, and traffic black spot identification and cause analysis are the research objectives. The built environment feature dataset of urban road intersections is constructed from street view images by the semantic segmentation model, which is used to identify and analyze the causes of traffic black spots at urban road intersections. The main contributions of this study are as follows:

(1)

By capitalizing on the characteristics of the built environment surrounding urban road intersections, and employing techniques such as random forest modeling, traffic black spots are identified for each individual road intersection;

(2)

To enhance the methodology for detecting traffic black spots, a graph network structure is proffered as a valuable approach. Recognizing the spatial correlation intrinsic to traffic black spots, this research adopts road intersections as nodes and establishes interconnections between adjacent road intersections as edges within the graph. By incorporating the built environment characteristics of the urban setting as node attributes, a graph network topology is constructed. Subsequently, a graph convolutional neural network is employed to carry out traffic black spot identification;

(3)

Starting from the characteristics of the built environment of urban road intersections, we analyze their correlation with the formation of traffic black spots, further expanding the research ideas in this direction.

2. Related Work

2.1. Previous Research on the Urban Built Environment

In early studies, the descriptions of the built environment were mainly based on site or community audits, which were time-consuming, labor-intensive, and could not be applied at large spatial scales [17,18,19]. Some scholars use GIS tools to describe the urban environment [20]. With the development of network mapping technology, urban streetscape images along road networks have been produced, giving rise to new ideas for the study of urban built-up area environments. Zhang et al. [13] used street view images to portray the urban visual environment, mining street view picture information with the help of the deep learning semantic segmentation model, creatively proposing scene semantic trees and scene expression vectors, which help to qualitatively understand and quantitatively analyze the cityscape. Through their research, Rundle et al. [20] demonstrated through research that Google Street View images can be used to evaluate the audit of the neighborhood environment, which has the characteristics of high efficiency and low cost. Bromm et al. [21] used Google Street View images to construct road environment features instead of collecting data in the field to evaluate the environmental features, and the results of the study proved that using street view images to evaluate the road environment is a reliable method. Badland et al. [22] assessed bicycling- and walking-related built environment features through a field survey and a Google Street View image analysis of 48 street segments in Oakland. The results show that the two methods demonstrated consistency across most of the assessed features, suggesting that Google Street View imagery can be a useful supplement to field surveys for assessing urban environmental features. Wilson et al. [23] showed that the average person can clearly interpret 84% of the features in the environmental factors represented by Google Street View; and Clarke [24] and others have shown that Street View images can provide some reliable indicators for evaluating local food environments and recreational facilities.

By evaluating the built environment of roadway segments and intersections, it was found that factors such as pedestrian and vehicular traffic flow, roadway width, and lighting have a definite impact on drivers [25], which can increase the probability of traffic accidents. In addition, factors such as medians and trees can have an impact on pedestrians on the roadway [26], which can lead to crashes between pedestrians and motor vehicles. Hanson et al. [27] used street view imagery to examine associations between pedestrians and roadway infrastructure features and crash severity, which provided recommendations for pedestrian safety. Johnson et al. [28] used Google Street View photos to assess the impacts of real-world car accidents and found that car drivers may have a greater likelihood of being injured. Wang et al. [29] then analyzed the built environment characteristics of traffic black spot areas and found that most of the areas with traffic black spots had bad roads, while in areas with fewer traffic accidents, the roads were wide and clean, and both motor vehicle and pedestrian densities were lower. In addition, Qin et al. [30] found that there is a connection between the built environments of neighboring road intersections based on the urban built environment using the graph convolutional neural network (GCNN) model, which can affect the probability of traffic congestion at the intersection. All of these studies have laid the theoretical and technical foundation for this paper to identify traffic black spots based on street view image data.

2.2. Previous Research on Semantic Segmentation of Urban Street View Images

Traditional image segmentation methods include two types of segmentation methods, namely edge-based segmentation methods and threshold-based segmentation methods. The threshold-based segmentation method consists in dividing each pixel in the image into different categories according to the gray level, which is suitable for images with different background and target gray levels. Edge detection-based methods, on the other hand, determine the edges of the region by certain detection methods. The detection principle lies in the fact that the gray values at the edges of an image are generally discontinuous, and such drastic changes can be detected by a simple derivation. Traditional image segmentation methods tend to use only the gradient histogram or the color space of the image, and these methods tend to be limited to specific segmentation tasks and are not very general.

There have been new breakthroughs in image semantic segmentation methods, and a seminal study is the fully convolutional neural network (FCN) model proposed by Long et al. [31], which is transformed into a fully convolutional network by deleting the fully connected layer in the convolutional neural network and then transforming it into a fully convolutional network. Based on the FCN model, the Seg Net model proposed by Badrinarayanan et al. [32] makes an innovation in the up-sampling operation; instead of using the hopping connection structure of the FCN model, the image detail information is recovered by means of maximally pooling the indices during the decoding process. The U-Net model proposed by Ronneberger et al. [33], on the other hand, is based on the shape of the network results. In contrast to the FCN model, the U-Net model supplements the detailed information by using more shallow features. Another more important approach is Deeplab-v2 proposed by Chen et al. [34]. The PSP Net model proposed by Zhao et al. [35] introduces more contextual information features through feature fusion operations and global average pooling operations based on the FCN model. Fu et al. [36] proposed the dual attention network (DANet) as a way to adaptively integrate local features and their global dependencies, which relatively significantly improves the classification results.

2.3. Study of Black Spots on Urban Roads

Since 1940, Western countries have carried out research on the identification of traffic accident black spots according to their own road safety conditions and objectives [36,37,38,39,40]. As shown in

Table 1. Definition of road black spots in Western countries.

Country (Region)	Definition of Road Black Spots
Austria	Locations with at least three traffic fatalities and injuries in a 3-year period and an accident factor of not less than 0.8
Denmark	Locations where the number of accidents exceeds the number predicted by statistical laws
Belgium	Locations with three or more traffic accidents in a 3-year period, or with a priority factor S exceeding 15
Germany	Accident black spots are defined by the number of accidents on a road section or at a location that exceeds the identification threshold within a certain period of time
Hungary	Road segments within the urban area with at least four traffic accidents in 3 years and a length of up to 1000 m, or outside the urban area with at least four traffic accidents in 3 years and a length of up to 100 m
Norway	Any point with at least four traffic fatalities and injuries within 5 years and not exceeding 100 m in length, and road sections with at least ten fatalities and injuries within 5 years and not exceeding 1000 m in length
Portugal	Intersections are road segments with at least five crashes with a severity index of more than 20 in a given year and a length of no more than 200 m
Switzerland	Locations where the accident rate exceeds this range by calculating the average accident rate for road segments and intersections and setting a maximum probability of deviation

, although the economic development and traffic operation environment of each country is different, the definition and understanding of accident black spots is also different, but at present, most of the countries are still adopting the accident number method or the accident rate method as the only means of accident black spot identification, i.e., through statistics and other methods of processing and analyzing the data and combining them with the actual road safety conditions and other factors to set up the accident black spot identification thresholds, and when the number of accidents or the accident rate exceeds the threshold, it will be defined as an accident black spot [41,42].

However, the accident number method and the accident rate method have greater limitations; not only is the threshold setting highly subjective, but the identification cycle is long and the required accident data are huge, so they cannot be widely used in urban traffic accident black spot identification.

The research on black spot identification of road traffic accidents has made remarkable progress, constantly promoting the development of theory and practice. The research teams used a variety of technical tools, such as a hierarchical Bayesian model [43], safety potential evaluation method [44,45,46], an innovative method combining dynamic segmentation and the DBSCAN algorithm [17], exhaustive statistical analysis based on Poisson and negative exponential distributions [18], classical Bayesian method [19], gray clustering model [20], kernel density estimation method [21], and a combination of Geographic Information System (GIS) and firefly algorithms and black spot identification algorithms based on support vector machines with maximum classification boundaries and deep learning networks [22] to effectively identify and manage accident black spots. In addition, the research also involves accident distribution patterns and spatio-temporal characterization, driver behavior studies, and recommendations to enhance road safety by improving traffic infrastructure, rationally planning road signs and markings, and optimizing guardrail design [23]. These results not only improve the efficiency of black spot identification, but also provide practical optimization paths for accident risk assessment and safety management.

In summary, despite the continuous theoretical improvement in the field of road traffic accident black spot identification, there are challenges in practical applications, especially in meeting the demand for short-term rapid identification of accident black spots on urban expressways. In addition, the subjectivity of threshold setting and segmentation of roadway unit length in the recognition process limits the ability of these techniques to accurately map the objective spatial distribution of accident black spots.

3. Study Area and Data

3.1. Study Area

The study area selected for this paper is the Kowloon Peninsula of Hong Kong, located in the northern part of the Hong Kong Special Administrative Region, a major administrative region and commercial center of Hong Kong, with a total area of about 46.93 square kilometers, including five districts, namely Yau Tsim Mong, Sham Shui Po, Kowloon City District, Wong Tai Sin, and Kwun Tong, which is one of the three major districts of Hong Kong (Figure 1). With a more developed economic and transportation network, and many famous scenic spots as well as major roads, it is an important convergence of logistics and people flow. As a typical urban area in Hong Kong, the Kowloon Peninsula has a very high traffic density and traffic congestion is very common, especially during the morning and evening peak periods.

3.2. Research Data Acquisition and Processing

The dataset used in this study is a self-constructed dataset, which consists of three parts:

Road network data: The road network data need to be obtained on the OpenStreetMap (OSM) website by itself, and the original data type obtained is in the OSM format, which needs to be converted into the SHP format and then processed.

Traffic Accident Data: The traffic accident data include a total of 10,392 data points with latitude and longitude coordinates in the Kowloon Peninsula area of Hong Kong for two years in 2014 and 2015, and the distribution is shown in

Table 2. Number of traffic accidents distributed in various districts of Kowloon Peninsula.

District	Number of Traffic Accidents
Wong Tai Sin District	1330
Kowloon City District	2158
Kwun Tong District	2402
Sham Shui Po District	1766
Yau Tsim Mong District	2736

.

Street View Image Data: The street view image data were collected using the Tencent Street View platform with Python 3.8 crawler technology, and a total of 8044 street view sampling points were selected and 32,176 street view images were acquired.

3.2.1. Road Data Acquisition and Processing

The road network map of the Kowloon area of Hong Kong obtained from the OpenStreetMap (https://openmaptiles.org/; Access Date: 1 December 2023) website contains geographic information about the roads, streets, highways, bridges, tunnels, sidewalks, and other transportation facilities and buildings in the area. The obtained road network of the Kowloon Peninsula area of Hong Kong is shown in Figure 2.

The obtained road network attributes of the Kowloon Peninsula of Hong Kong are shown in

Table 3. OpenStreetMap road network attributes.

Property Name	Hidden Meaning
fid	The “fid” can be used to uniquely identify an element on the map.
shape	The “shape” attribute is usually used to indicate the shape of a line element.
length	The “length” attribute is usually used to indicate the length of a road line element.
highway	The “highway” attribute is usually used to indicate the type of road.
lanes	The “lanes” attribute is usually used to indicate the number of lanes in a roadway.
name	The “name” attribute is usually used to identify the name of a location.

.

The value of the “lanes” attribute indicates the number of lanes in the road, including all lanes in a two-way road. For example, a two-way road with two lanes can be labeled as “lanes = 4”, where two lanes are used in each direction.

The value of the “shape” attribute indicates the shape of the line elements; in this study, the road network is folded lines.

The value of the “highway” attribute indicates the class and type of the road, and there are mainly the following categories:

(1)

Motorway: highway;

(2)

Trunk: trunk highway;

(3)

Primary: major roads;

(4)

Secondary: secondary roads;

(5)

Tertiary: third-class roads;

(6)

Residential: residential roads;

(7)

Unclassified: unclassified roads.

The contents of the road network data attribute table are shown in

Table 4. Table of road network data attribute values.

fid	Shape	Length	Highway	Lanes	Name
1	broken line	715.644	trunk	3	Long Xiang Street
2	broken line	832.543	primary	4	Cheung Sha Bend Road
3	broken line	126.559	unclassified	1	Lung Ping Road
4	broken line	66.248	residential	1	Fuk Wah Street
5	broken line	483.324	secondary	2	Cornwall Street

.

OpenStreetMap data are contributed by community volunteers, so there may be an incomplete coverage of the map data, and at the same time, according to the research objectives of this paper, it is also necessary to carry out a secondary processing of the road network, so as to ensure that the adopted road network can be used to carry out the research on the identification of traffic black spots at road intersections.

The process of secondary processing of the road network mainly includes the following steps:

Road network connectivity processing: in some cases, the road network may have the problem of discontinuity, that is, some areas cannot be reached or can only be reached through longer paths, which will affect the recognition results of traffic black spots, so it is necessary to carry out connectivity processing.

Road network segmentation: road network segmentation consists in cutting the road network into small segments, which can be nodes, edges, or regions, which can improve the computational efficiency and accuracy of the traffic black spot recognition algorithm. In this study, the road network segmentation is performed at road network intersections.

Remove redundant road network lines. Turn two-way roads into unidirectional lines, circular intersections into unidirectional cross-connection-type lines, etc. Such treatment can simplify the road network, avoid repeated calculations, and unify the road representation.

Using MapGIS 10.6 software to process the road network, verify the road connectivity, and eliminate the disconnected roads, unconnected roads, and branch lines. Secondly, the redundant road network lines are deleted, so that the road network map used in this study can be initially obtained, and a road network processing example is shown in Figure 3, Figure 4 and Figure 5.

Deleting redundant branch roads in the road network, as shown in Figure 3, simplifies the road network and does not interfere with the research experiments.

As shown in Figure 4, the complexity of roadway intersection routes in the roadway network is generally simplified to one adjacent route, with the roadway intersection becoming a simpler cross-connection type.

As shown in Figure 5, the two-way lanes in the road network are simplified into one-way lanes, while the lines are connected to form road intersections.

After processing the above steps, the road network map used in the experiment can be initially obtained, as shown in Figure 6.

The road network shown in Figure 7 is processed using MapGIS software to interrupt intersectionality from road intersections to form road intersection nodes, and Figure 7 demonstrates the results of localized road intersection node processing in the Kowloon Peninsula area.

The road network road intersection nodes in the entire Kowloon Peninsula area were processed to obtain a total of 2519 intersections, and the results are shown in Figure 8.

After the above steps, the final road network map and road intersection node location map used in this experiment are obtained.

3.2.2. Traffic Accident Data Acquisition and Processing

(1)

Traffic accident data acquisition

The traffic accident data used from 2014 to 2015 were obtained from the Hong Kong Traffic Management Department. The data attributes include latitude and longitude coordinates, weather conditions, and accident severity, etc., and there are a total of 10,392 entries.

The acquired traffic accident data attributes are shown in

Table 5. Table of attributes of traffic accidents.

Property Name	Hidden Meaning
fid	The “fid” can be used to uniquely identify an element on the map.
shape	The “shape” attribute is usually used to indicate the shape of a line element.
year	The “year” attribute is usually used to indicate the year of a traffic accident.
severity	The “severity” attribute is usually used to indicate the severity of the crash.
latitude	The “latitude” attribute is usually used to indicate the longitude of the location of the accident.
longitude	The “longitude” attribute is usually used to indicate the latitude of the location of the accident.
weather	The “weather” attribute is usually used to indicate the weather conditions of the day.

.

According to the research objectives of this experiment, focusing on selecting fid attribute, year attribute, latitude attribute, and longitude attribute, the content of the attribute table of traffic accident data is displayed, as shown in

Table 6. Table of values of traffic accident data attributes.

Fid	Shape	Year	Latitude	Longitude
1	point	2014	22.300664	114.167181
2	point	2014	22.252020	114.135996
3	point	2014	22.331050	114.205178
4	point	2015	22.306150	114.224390
5	point	2015	22.320042	114.208088

.

(2)

Traffic accident data processing

The main purpose of the experiment is to remove the traffic accident data with missing latitude and longitude attributes, and it also includes some steps such as filling the missing values and processing the outliers (the latitude and longitude coordinates are not on the road network, etc.). After processing, a total of 8374 pieces of experimentally available traffic accident data were obtained.

A buffer zone will be established using MapGIS 10.6 with the road intersection node as the center of the circle, as shown by the black dot in Figure 9, with a radius of 100 m.

Then, combining the 8374 traffic accident data points, using MapGIS to import them into the established road intersection buffer and unify the coordinate system, the final results are displayed as shown in Figure 10.

If the number of traffic accidents in this buffer zone exceeds twenty or more, this defines the road intersection as a traffic black spot. After counting, there are 584 traffic black spots in total, and the final results are displayed as shown in Figure 11, where the red nodes denote the traffic black spots, which are replaced by the number 1, and the blue nodes denote the non-traffic black spots, which are replaced by the number 0.

The MapGIS platform was also used to count the distribution of traffic black spots in the five districts of Kowloon Peninsula, as shown in

Table 7. Table showing the number of traffic black spots distributed in various districts of the Kowloon Peninsula.

District	Number of Traffic Black Spots
Wong Tai Sin District	24
Kowloon City District	140
Kwun Tong District	49
Sham Shui Po District	80
Yau Tsim Mong District	294

.

3.2.3. Streetscape Image Data Acquisition

In this study, we utilize the street view image data provided by the Tencent Street View platform to extract street-level road environment features. By analyzing the road network data, we selected a total of 8044 street view sampling points with an average sampling interval of about 100 m. Four images were captured at each sampling point with orientations of 0°, 90°, 180°, and 270°, constituting a 360° panoramic view to fully characterize the environment at each sampling point. We acquired a total of 32,176 street view images of Hong Kong’s Kowloon Peninsula, each with a size of 1024 × 1024 pixels. The final sampling point location map and the acquired street view images are shown in Figure 12, with the green nodes representing the street view sampling points.

3.2.4. Streetscape Image Data Processing

In previous research, most scholars would apply traditional machine learning models (e.g., decision trees, convolutional neural networks) to classify street view images at the image level until the emergence of FCN [18] in 2015, which transforms the last layers of convolutional networks from the last few fully connected layers to the fully convolutional layers, realizing a breakthrough from image-level to pixel-level classification. In this study, we will use FCN, a classical semantic segmentation model, to recognize the pixels in street view images, including more than 100 features such as trees, buildings, billboards, roads, sidewalks, signals, and vehicles, and on the basis of which we will calculate the percentage of each feature in each group of street view images.

The basic structure of the entire FCN network is shown in Figure 13.

In this study, a deep learning fully convolutional network trained on the ADE_20K dataset, which has a pixel contrast accuracy of 81% in the training dataset [47], is used to perform semantic segmentation on the pixel level of the acquired street view images to obtain more than 100 street-level environmental features such as green trees, buildings, roads, and sky. The final input result image of the semantic segmentation model is at 512 × 512 pixels and the semantic segmentation result image is shown in Figure 14.

As can be seen from the result of semantic segmentation, in Figure 14, the model has a better recognition effect, but there is still a slight error. For large-scale street scene recognition, such an error is acceptable. At the same time, the model will also obtain a picture element feature percentage data table, the attributes of which are shown in

Table 8. Semantic segmentation model element feature shared attributes.

Property Name	Meaning
svid	The “svid” is an integer that uniquely identifies the sampling point.
direct	The “direct” attribute is usually used to indicate the location of the angle at which the sample point was taken.
road	The “road” attribute is usually used to indicate the percentage of roads in an image.
building	The “building” attribute is usually used to represent the percentage of buildings in an image.
sky	The “sky” attribute is usually used to represent the percentage of sky in a picture.
tree	The “tree” attribute is usually used to represent the percentage of trees in an image.
streetlight	The “streetlight” attribute is usually used to represent the percentage of streetlights in an image.

.

The results of the feature occupancy of some elements of the semantic segmentation model for street scene images are displayed in

Table 9. Feature share of some elements of the semantic segmentation model for streetscape images.

Svid	Direct	Road	Building	Tree
37011007140313112032600	Left	36.90986633	25.16593933	16.90559387
37011007140313112032600	Right	36.37313843	20.44830322	12.4294281
37011007140313112032600	Front	13.66424561	6.915283203	28.41835022
37011007140313112032600	Back	23.45466614	32.78198242	21.04797363
37011007140313112045500	Left	34.24072266	35.12573242	7.035064697
37011007140313112045500	Right	33.47587585	38.87634277	6.537246704
37011007140313112045500	Front	20.10383606	66.08657837	0
37011007140313112045500	Back	26.7665863	17.46864319	36.40022278

.

As shown in Table 9, the svid attribute is the unique number of the street scene sampling point; while each sampling point has four pictures, according to the different direct attributes to differentiate the direction, in the process of data processing, it is necessary to calculate the average of the four-picture occupancy of the same sampling point, to obtain the unique elemental feature occupancy data of each sampling point. At the same time, zero is assigned to the picture occupancy of the feature to avoid null values. The final elemental feature ratio of the semantic segmentation data of street view images for each sampling point is shown in

Table 10. Feature share of data elements for semantic segmentation of streetscape images.

Svid	Road	Building	Tree	Sky	Bus	Sidewalk
37011007140313112009700	27.600479	21.327877	19.700336	11.957073	3.814697	4.408550
37011007140313112032600	28.646755	39.389324	12.493134	4.594326	0.020790	2.093601
37011007140313112045500	36.336517	42.506981	2.013302	11.781406	0.026512	2.521706
37011007140313112153000	21.691227	45.594501	1.884079	12.111092	0.000000	7.178211
37011007140313113823900	27.237415	51.236629	0.212097	3.229713	0.012398	10.075188
37011013140704150800500	31.939507	27.535725	11.226654	6.665039	0.362587	1.970863
37011013140704150802500	28.536606	40.852642	4.955864	7.218075	0.582981	9.016323
37011013140704152537000	25.959492	52.578068	0.099945	5.836773	0.000000	7.768345

.

The coordinates of each sampling point are imported into the road intersection buffer established in the previous section, and the number of sampling points falling into the buffer and the percentage of each element’s characteristics are counted as a means of depicting the built environment of the road intersection, and Figure 15 shows the distribution of the sampling points of the streetscape in the road intersection buffer.

The number of sampling points distributed in the buffer zone of roadway intersections was counted and the results are shown in Figure 16.

Statistical results show that each buffer has a sufficient number of sampling points, indicating that the results of semantic segmentation of street view images can be used to represent the built environment of road intersections in the study area, which lays the foundation for the smooth progress of the experiments on traffic black spot recognition.

Finally, the results of semantic segmentation of street scene images in each buffer are processed by data connection, filling, etc., and the final part of the built environment of some road intersections is shown in

Table 11. Characteristic percentage of road built environment elements at road intersections.

fid	Road	Building	Tree	Sky	Bus	Sidewalk
0	36.90986633	25.16593933	16.90559387	10.83641052	4.64630127	2.752304077
1	36.37313843	20.44830322	12.4294281	18.67103577	0	3.379058838
2	13.66424561	6.915283203	28.41835022	14.63470459	10.61248779	7.374191284
3	23.45466614	32.78198242	21.04797363	3.686141968	0	4.128646851
4	34.24072266	35.12573242	7.035064697	9.188842773	0.025939941	0
5	33.47587585	38.87634277	6.537246704	9.12437439	0	0
6	26.7665863	17.46864319	36.40022278	0.064086914	0.057220459	0.422286987
7	38.45748901	40.20271301	6.348419189	8.185195923	0.006103516	3.675079346
8	28.87039185	43.71452332	0.653076172	15.32402039	0.099945068	3.146362305
9	41.0949707	33.22982788	0.941085815	18.69354248	0	1.778030396
10	42.93022156	3.246688843	0.149536133	38.7336731	0.001144409	0.415802002
11	45.43418884	0.122070313	1.889419556	39.11361694	0.030136108	0.311660767

.

4. Methodology

4.1. Machine Learning Models

4.1.1. SVM

SVM is an algorithm under a supervised learning model, commonly used in classification problems. The earliest SVM model is a linear classifier dedicated to binary classification, with a statistical theory foundation. The SVM method maps the sample space into a high-dimensional or even infinite-dimensional eigenspace (Hilbert space) through a nonlinear mapping p, so that a problem that is nonlinearly differentiable in the original sample space is transformed into a problem that is linearly differentiable in the eigenspace. Simply put, it is lift and dimension linearization. Due to classification, regression, and other problems, it is likely that the sample set in the low-dimensional sample space cannot be linearly processed, but in the high-dimensional feature space this can be achieved through a linear hyperplane linear division (or regression). Generally, while upgrading will bring the complexity of computation, the SVM method applies the expansion theorem of the kernel function and establishes a linear learning machine in high-dimensional feature space; therefore, compared with the linear model, not only does it almost not increase the complexity of the computation, but to a certain extent it also avoids the “dimensionality catastrophe”. Different SVMs can be generated by choosing different kernel functions, and the following four kernel functions are commonly used: (1) linear kernel function $K\left(x,y\right)=x·y$; (2) polynomial kernel function $K\left(x,y\right)={\left[\left(x·y\right)+1\right]}^d$; (3) radial basis functions $K\left(x,y\right)=\sqrt{\frac{-{\left|x-y\right|}^2}{d^2}}$; (4) two-layer neural network kernel function $K\left(x,y\right)=tanh\left(a\left(x,y\right)+b\right)$. Figure 17 is the construction flowchart of the SVM model for binary classification of traffic black spot recognition.

The SVM model is trained using data with selected features, and during the training process, appropriate parameters and kernel functions need to be selected to optimize the performance of the model. In this process, the model parameters need to be selected, and the more important parameters are selected as follows:

(1)

The choice of kernel function. Depending on the different SVM model kernel functions, the accuracy of the SVM model will be different; generally, it will be a linear kernel function, Gaussian kernel function, Sigmoid kernel function, and so on;

(2)

Optimization of kernel parameters. The main parameters are the kernel parameter g and the penalty factor c. The penalty factor indicates the importance of misclassification, and the parameter g indicates the size of the impact of a single sample, and the smaller the value, the greater the impact.

The specific description of the model parameters is shown in

Table 12. Important parameters of the SVM algorithm model (described in words).

Parameter Name	Description
kernel	Linear kernel function
c	Penalty coefficients for error terms
gamma	Kernel function coefficients

.

The SVM function package in the machine learning Sklearn toolkit was used to model the SVM algorithm, and the model accuracy was higher when the kernel function was chosen to be linear, in addition to the other parameters of the algorithmic model (

Table 13. Selection of model kernel function for SVM algorithm.

Kernel Value	Accuracy
Linear	81.75%
Rbf	76.83%
Sigmoid	71.58%

), which were appropriately tuned using the grid search method.

According to

Table 14. Values of important parameters of the SVM algorithm model (direct textual description).

Parameter	Retrieve a Value
kernel	linear
c	10
gamma	auto

, the accuracy of the model is high when the kernel function is chosen to be linear. The appropriate kernel function is chosen to find the optimal classification hyperplane, and the support vector set of each sample feature is obtained to form the discriminant function of classification.

4.1.2. GBDT Model

GBDT is a common basic machine learning algorithm based on boosting iterative decision trees that can be used for classification, regression, and other problems. The steps of the algorithm first train a decision tree with an initial sample, and then train the next decision tree based on the residuals of the previous decision tree, thus gradually reducing the iterative residuals, thus improving the iterative results and finally constructing a new combinatorial model. When all the convergence criteria are satisfied, all the base models are obtained. In the GBDT classification algorithm, the final output is discrete categories rather than continuous values, so we cannot fit the category output error directly from the output categories. To solve this difficulty, we fit the loss function using a method similar to the log-likelihood loss function used in logistic regression. The dataset of this experiment has a relatively large number of features, and the traditional linear model cannot achieve a better classification effect. GBDT, as one of the representative models of integrated learning, consists in training multiple weak learners and synthesizing their classification results to form a strong learner. Figure 18 shows the flowchart of the construction of the GBDT model for the binary classification of the traffic black spot recognition.

Figure 18 consists of multiple weak learners. X represents the feature attributes, Y represents the prediction results of the weak classifiers, and the results of multiple weak learners are combined to obtain the final classification results.

The traffic black spot recognition model is constructed using the GBDT model, which is an integrated learning method based on a decision tree, and its main idea is to construct a strong classifier by combining multiple weak classifiers. When building the GBDT model, it is necessary to select the hyperparameters such as the depth of the decision tree and the number of leaf nodes, as well as to adjust the parameters such as the learning rate and the loss function. The selection and description of model parameters are shown in

Table 15. Important parameters of the GBDT algorithm model.

Parameter Name	Hidden Meaning	Usage
loss	Loss function	Deviance log-likelihood function
max_depth	Growth depth of each decision tree	Exponential loss function
n_estimators	Number of iterations of the weak learner	A positive integer
learning_rate	Learner learning step size	A positive integer

.

Gradient Boosting Classifier in the machine learning Sklearn toolkit was used to model the GBDT algorithm, with loss as the loss function, and then 2519 data points were randomly divided into test and validation sets, and the parameters of the algorithmic model were appropriately tuned using the grid search method, including the parameters of n_estimators, max_depth, and so on.

As shown in

Table 16. Comparison table of loss fetch value tuning reference.

Loss Value	Model Accuracy
Deviance	82.51%
Exponential	81.87%

, the values of the loss parameter are set to deviance and exponential, respectively, and it can be seen that the model is slightly more accurate when the value of the loss parameter is set to deviance.

According to Figure 19, the accuracy of the model is higher when the n_estimators parameter in the GBDT model is set to 23.

Moreover, the results of the grid search method show that the accuracy of the model is better when the other parameter, max_depth, is set to 7 and when min_samples_leaf is set to 1.

Table 17. Values of some parameters of the GBDT model.

Parameter Name	Value	Clarification
max_depth	7	Growth depth of each decision tree
n_estimators	23	Number of decision trees
learning_rate	0.05	Learning rate
min_samples_leaf	1	Minimum number of samples for leaf nodes

shows some of the parameter settings of the GBDT model, on the basis of which the GBDT model is used for the identification study of traffic black spots.

4.1.3. Random Forest

Random forest is one kind of cluster classification model, which consists in building a forest in a random way. The forest consists of a significant number of decision trees, and there is no association between each decision tree. After obtaining the random forest model, when a new sample enters the random forest for each decision tree to judge, the bagging collection strategy is relatively simple, for classification problems usually use the voting method to obtain the category with most votes for the final model output. Firstly, from the original training dataset, the bootstrap method is applied to randomly draw k new self-help sample sets with put-back, and from this, k decision trees are constructed, and the samples that are not drawn each time make up K out-of-bag data. Second, there are n features, and mtry features are randomly drawn at each node of each tree, and by calculating the amount of information embedded in each feature, one of the features with the most classification ability is selected for node splitting. Third, each tree is maximally grown without any clipping. Fourth, the generated multiple trees are formed into a random forest, which is used to classify the new data, and the classification result is based on how many tree classifiers are voted for. For classification problems, the random forest model can naturally handle most classification tasks, and it is a natural and nonlinear classification model. Figure 20 shows the flowchart of the construction of the random forest model for binary classification of traffic black spot identification.

Next, we discuss model training using the random forest algorithm. Random forest is an integrated learning method that combines multiple decision trees to form a more accurate classifier. Some parameters need to be set during training, such as the number of trees, the maximum depth of each tree, and the minimum number of samples for leaf nodes, as shown in

Table 18. Important parameters of random forest.

Parameter Name	Clarification
bootstrap	Indicates sampling with put-back
max_depth	Growth depth of each decision tree
n_estimators	Number of decision trees
learning_rate	Learning rate

.

In the tuning process of random forest models, there are several principles:

Training set random sampling selection. In this study, the bagging method is chosen for sampling, which can prevent overfitting and at the same time has a relatively fast training speed.

Selection of training objective function. There are mainly two kinds, information gain and Gini coefficient; the Gini index tends to select features with more values as split features, while the information gain index easily causes overfitting when using high latitude data. According to the data of this experiment, the Gini coefficient index is selected.

Parameter selection. The max_depth parameter indicates the maximum depth of the tree, generally set as max_depth = 3, and then to try to increase the depth, the min_samples_leaf parameter indicates the minimum number of samples specified to be contained in each leaf node. If this parameter is set too low it will cause overfitting, so generally it is set as min_samples_leaf = 5, and followed by up and down adjustment.

In this paper, we use the Random Forest Classifier in the machine learning toolkit to construct the random forest algorithm model, which can set a series of parameters to tune the model, we use the logloss function for the calculation of the loss function, and in this process, randomly classify the test set and validation set on 2519 pieces of data, and we use the grid search method. The algorithm model is tuned, after which the max_depth parameter is set to 13 and the max_features parameter is set to 46. Figure 21 shows the tuning diagram of the random forest model parameter n_estimators as set to 15, and

Table 19. Values of some important parameters of random forest.

Parameter Name	Retrieve a Value	Clarification
bootstrap	True	Indicates sampling with put-back
max_depth	13	Growth depth of each decision tree
n_estimators	15	Number of decision trees
learning_rate	0.05	Learning rate

shows the values of some parameter settings of random forest.

4.1.4. Fully Connected Neural Network

Figure 22 is a flowchart of the construction of a fully connected neural network model for the binary classification of traffic black spot identification.

Figure 22 consists of an input layer, a hidden layer, and an output layer, where x represents the feature attributes and y represents the classification results of traffic black spots.

A fully connected neural network model is constructed, and data are fed into the model to train the model for the classification task. A fully connected neural network can consist of multiple hidden layers, each layer containing multiple neurons and each neuron being connected to all the neurons in the previous layer. Some parameters need to be set during training, and after network tuning, some of the hyperparameters of the fully connected neural network proposed in this study are set at the following values:

(1)

Input layer vector dimension: 139 features of each intersection semantic segmentation will be input;

(2)

Number of hidden layers: the number of hidden layers is three;

(3)

Neurons: the neurons of the three hidden layers are 79, 17, and four, respectively;

(4)

Activation function: the activation function is chosen as ReLU, which is less likely to have gradient explosion or gradient disappearance;

(5)

Learning rate: the learning rate is between 0.005 and 0.01;

(6)

Back-propagation algorithm: after the neural network calculates the error according to the error function in the process of iteration. It needs to return the result to the end of input weight adjustment, which needs to consider the model effect and speed. In this study, Adam’s algorithm is chosen as the back-propagation algorithm for this model.

The values of some parameter settings of the fully connected neural network are shown in

Table 20. Hyperparameter settings for fully connected neural networks.

Parameter Name	Grade	Value
Output Layer Vector Dimension	n	139
Number of hidden layers	l	3
Hidden layer neurons	o	79, 17, 4
Activation function	F(x)	ReLu
Learning rate	$\eta$	0.0001
Algorithm	A	Adam

.

4.2. Traffic Black Spot Recognition Model Construction Based on Graph Convolutional Neural Network Model (GCNN)

In this experiment, roads can be considered as nodes of the graph, and connections between neighboring roads are considered as edges of the graph. Then, the road features from the street view image data source can be input into the GCNN as node features. With the GCNN model, the following operations can be performed:

(1)

Neighbor node feature aggregation: the GCNN model can fuse the information between adjacent roads by aggregating the neighbor node features of each node. This allows for the features of neighboring roads to be used to better portray the urban built environment;

(2)

Contextual information learning: the GCNN model can learn the contextual information of each node, including the attributes, shapes, and traffic flows of neighboring roads. This helps to synthesize the relationship between roads to further portray the urban built environment;

(3)

Feature representation learning: the GCNN model can learn the representation of road features through multi-layer convolutional operations. Each GCNN layer can extract more advanced feature representations to capture more complex information about the urban built environment.

In the following, specific steps are described for the study of traffic black spot recognition at road intersections using the GCNN model.

(1)

Introduction of the dataset

This experiment requires the use of street view image semantic segmentation data and the corresponding road intersection data. There are 2519 intersections in total. The intersections are mainly divided into two categories: one is traffic black spot intersections and the other is non-traffic black spot intersections. Each intersection is connected to at least one other intersection, and the connectivity information between the 2519 intersection data points is obtained through Python data analysis. Each intersection is treated as a graph node, the connection of adjacent road intersections is the edge of the graph, the street view semantic feature information of this intersection is treated as the node information of the corresponding node, and there are 139 node features, and the value is the percentage of this feature in the image of this intersection, and the neighbor matrix is constructed by obtaining the adjacency table based on the connectivity between nodes to obtain the structure of the graph network, as shown in Figure 23, where the red color represents the intersections, and the gray color represents street view image sampling points. The final data form is as follows:

(1)

x: node features with dimension 2519 × 139 and type np.ndarray;

(2)

y: node corresponding labels, including two categories (traffic black spots and non-traffic black spots), np.ndarray;

(3)

Adjacency: adjacency matrix, dimension 2519 × 2519, type scipy.sparse.coo_matrix;

(4)

train_mask, val_mask, test_mask: the same mask as the number of nodes, which is used to classify the training set, validation set, and test set. When the node belongs to the training set or validation set or test set, the corresponding position is True, and otherwise it is False.

In GCNN node classification, how to choose and construct the feature matrix and adjacency matrix is crucial, because they directly affect the performance of the model.

(2)

Construction of the GCNN model

First of all, the data should be processed to obtain the node features and labels, adjacency matrix, training set, validation set, and test set. To normalize the data, to increase the self-connections, and to perform the definition of the graph convolution layer, the node input dimensions are 139 dimensions, the output is two dimensions, no bias is used, and at the same time, since the adjacency matrix is a sparse matrix, the sparse matrix multiplication method is used. Defining the graph convolution model, the hidden layer has a total of two layers, as shown in Figure 24.

(3)

Selection of GCNN model parameters

PyTorch-Geometric is a PyTorch-based graph neural network library specialized in processing graph network structures. It provides a series of tools for processing graph data, including data loading, data preprocessing, graph convolutional neural network model construction, model training, and model evaluation.

In this experiment, PyTorch, as well as the graph neural network library PyTorch-Geometric (PYG), are used. The Adam function is used for gradient descent, the loss function is used for the computation of the loss value, and the activation function is used for the calculation of the loss value by using ReLu. The hyperparameters are optimized, and the definitions of the hyperparameters are shown in

Table 21. Graph convolutional neural network (GCNN) hyperparameter values.

Parameter Name	Grade	Value
Output layer vector dimension	n	139
Number of hidden layers	l	2
Hidden layer neurons	o	40, 6
Activation function	F(x)	ReLu
Learning rate	$\eta$	0.01

.

5. Experimental Results and Analysis

5.1. Indicators for Model Evaluation

The data partitioning is kept consistent for all machine learning models and is divided into a training set and a test set. This ensures the comparability of the experimental results. When constructing the algorithmic model, all the sample data (a total of 2519 entries, each with 139 features) in the two categories (traffic black spots and non-traffic black spots) are randomly divided into 80% as the training set and 20% as the test set.

The main evaluation metrics commonly used in machine learning classification models are the following: (i) Accuracy; (ii) Precision; (iii) Recall; and (iv) F1_score; these metrics can help us evaluate the performance of the algorithm in different situations as well as provide valuable feedback when improving the algorithm. A confusion matrix is a matrix used to evaluate the performance of classification algorithms. In the binary classification problem of traffic black spot recognition at road intersections, the confusion matrix is a 2 × 2 matrix that contains the values True Positive (TP), False Positive (FP), True Negative (TN), and False Negative (False Negative, FN).

(1)

Definition of “Accuracy”: the ratio of the number of correctly categorized samples to the total number of samples. Accuracy is one of the most commonly used evaluation indices for classification models, but for the case of data imbalance, the accuracy may be misleading, and the expression is:

$$\mathrm{Accuracy}=\frac{\mathrm{T}\mathrm{P}+\mathrm{T}\mathrm{N}}{\mathrm{T}\mathrm{P}+\mathrm{T}\mathrm{N}+\mathrm{F}\mathrm{P}+\mathrm{F}\mathrm{N}}$$

(1)

(2)

Definition of “Precision”: the proportion of samples classified as positive that are truly positive. Precision is the evaluation of the ability of the classifier, i.e., the proportion of truly positive samples among all the samples determined as positive by the classifier, expressed as:

$$\mathrm{Precision}=\frac{\mathrm{T}\mathrm{P}}{\mathrm{T}\mathrm{P}+\mathrm{F}\mathrm{P}}$$

(2)

(3)

Definition of “Recall”: the proportion of all positive samples that are correctly categorized as positive samples. Recall is an index for evaluating the ability of a classifier to detect positive samples, i.e., the percentage of positive samples that can be detected by the classifier out of all real positive samples, expressed as:

$$\mathrm{Recall}=\frac{\mathrm{T}\mathrm{P}}{\mathrm{T}\mathrm{P}+\mathrm{F}\mathrm{N}}$$

(3)

(4)

Definition of F1 Score (F1$\_\mathrm{score}$): the F1 value is a metric that combines precision and recall, i.e., the reconciled average of precision and recall. The higher the F1 value, the better the classification performance of the model, and the expression is:

$$\mathrm{F}1\_\mathrm{score}=\frac{2\ast \mathrm{Precision}\ast \mathrm{Recall}}{\mathrm{Precision}+\mathrm{Recall}}$$

(4)

In the experiments of this paper, the precision rate represents the proportion of real traffic black spots among the number of all judged traffic black spots; the recall rate represents the ability of the maximization model to find all traffic black spots in the dataset, i.e., the proportion of the number of correctly judged traffic black spots to the proportion of the actual number of all traffic black spots; and the F1 score is a combination of the two metrics of the precision rate and the recall rate, and the higher the value of the F1 score, the better the classification performance of the model.

The data are randomly assigned and the ratio of data in the training set and test set is 4:1. SVM, GBDT algorithm, random forest algorithm, and fully CNN methods are utilized for training and evaluation of the data.

5.1.1. Evaluation of SVM Model

Table 22. Different model evaluation metrics and loss function values.

	Accuracy	Precision	Recall	F1_Score	Loss (Mean)	Loss (Max)	Loss (Min)	Loss (Std)
SVM	81.14%	74.12%	63.08%	41.86%	6.513	7.346	5.658	0.269
GBDT algorithm	85.58%	85.02%	74.32%	63.99%	4.636	5.795	3.741	0.430
Random forest algorithm	87.19%	86.12%	80.07%	75.06%	3.725	4.928	2.738	0.335
Fully CNN	84.54%	83.04%	73.57%	68.09%	4.884	5.973	3.634	0.421

demonstrates that the average Accuracy was 81.14%, the average Precision was 74.12, the average Recall was 63.08%, and the average F1 score was 41.86%. The average Loss value was 6.513, the maximum Loss value was 7.346, the minimum Loss value was 5.658, and the standard deviation of Loss value std was 0.269.

5.1.2. Evaluation of GBDT Model

Table 22 demonstrates that the average Accuracy was 85.58%, the average Precision was 85.02%, the average Recall was 74.32%, and the average F1 score was 663.99%. The average Loss value was 4.636, the maximum Loss value was 5.795, the minimum Loss value was 3.741, and the standard deviation of Loss value std was 0.430.

The values of the loss function for the GBDT model, with loss values reaching an average of 4.6 per round, represent a reduction in loss values relative to the SVM model.

5.1.3. Evaluation of Random Forest Model

Table 22 demonstrates that the average Accuracy was 87.19%, the average Precision was 86.12%, the average Recall was 79.07%, and the average F1 score was 73.06%. The average Loss value was 3.725, the maximum Loss value was 4.928, the minimum Loss value was 2.738, and the standard deviation of Loss value std was 0.335.

The values of the loss function for the random forest model, with the loss values reaching an average of 3.7 per round, represent a significant decrease relative to the SVM model.

5.1.4. Evaluation of Fully CNN Models

Table 22 demonstrates that the average Accuracy was 84.54%, the average Precision was 83.04%, the average Recall was 73.57%, and the average F1 score was 68.09%. The average Loss value was 4.884, the maximum Loss value was 5.973, the minimum Loss value was 3.634, and the standard deviation of Loss value std was 0.421.

From the specific values of the evaluation index of the fully connected neural network model, it can be found that the accuracy of the fully connected neural network model reaches 84%, and there is a decrease in the accuracy relative to that of the random forest model, because the deep learning model is more prone to overfitting due to the limitation of the amount of data, and the value of the F1 score index reaches 68%, and in general, as the most basic aspect of the deep learning field of the network model, the fully connected neural network model also shows better performance.

5.2. Comparative Results and Analysis of Machine Learning Model Performance

In the study of traffic black spot recognition for each road intersection, it was found that, among the four models, namely SVM, GBDT, RF, and fully connected neural network, the random forest model was the most effective for traffic black spot recognition, with an accuracy rate of 87%, and at the same time the value of the loss was relatively small, and the results of the evaluation of the various models are shown in

Table 23. Comparison of the results of the four model constructions.

Evaluation Indicators	Model Algorithm Name
	SVM	GBDT	Random Forest	Fully CNN
Accuracy	81.14%	85.58%	87.19%	84.54%
Precision Rate	74.12%	85.02%	86.12%	83.04%
Recall	63.08%	74.32%	80.07%	73.54%
F1_score	41.86%	63.99%	75.06%	68.09%
Loss Mean	6.513	4.636	3.725	4.884
Loss (max)	7.346	5.795	4.928	5.973
Loss (min)	5.658	3.741	2.738	3.634
Loss (std)	0.269	0.430	0.335	0.421

, which demonstrates that it is possible to utilize the characteristics of the built environment of the roads at each road intersection to conduct the study of traffic black spot recognition by means of the results of the above experiments.

The results of the above models show that, in the case of only considering the built environment of urban roads with only a single intersection, the random forest model is the most effective, with an accuracy rate of about 87%, while the precision rate and the recall rate both reach more than 80% at the same time, and the accuracy rates of the other models reach more than 80% as well. The above information fully exemplifies the feasibility of using street view images to identify the traffic black spots at road intersections.

However, all of the above machine learning models only recognize the road intersections in each study area as an independent individual, and have not yet considered the influence of the interconnections between the road intersections on the recognition of traffic black spots; therefore, in the next section, we will analyze the spatial correlation of the traffic black spots to show that there is a clustering effect between the traffic black spots.

5.3. Spatial Correlation Analysis of Traffic Black Spots

Spatial correlation is the correlation between phenomena or variables in geographic space. It involves the factor of geographic location and examines whether there is a correlation or similarity between neighboring regions in geographic space.

In previous studies, the spatial distributions of traffic black spots are generally not considered to be spatially independent of each other, and the probability of them being the same traffic black spot is higher when the physical environments are very similar. Therefore, this paper will analyze the spatial correlation of traffic black spots at road intersections in the Kowloon Peninsula region of Hong Kong through the use of spatial analysis methods such as the global Moran’s I.

Spatial correlation analysis of traffic black spots at road intersections in the Kowloon Peninsula region of Hong Kong can reveal their geospatial distribution patterns, clustering, and the degree of association with other factors. This is of great significance for transportation planning, traffic safety, and traffic management, and helps to improve the efficiency and safety of the transportation system.

5.3.1. Global Moran’s I

Global Moran’s I is a commonly used statistical index that can be used to measure the overall spatial correlation of traffic black spots at road intersections in Kowloon Peninsula, Hong Kong.

In this section, the study takes the Kowloon Peninsula of Hong Kong as the study area and the road intersections as the study units, and the MapGIS spatial statistics module is used to calculate spatial statistics and construct the Global Moran’s I to analyze the global spatial correlation of the distribution of traffic black spots. The global Moran’s I ranges from −1 to 1. Positive values indicate positive correlation, which means that there is a positive correlation in the distribution of traffic black spots at road intersections in the Kowloon Peninsula region of Hong Kong, and the closer the value is to 1, the stronger the correlation is. Negative values indicate negative correlation, and the closer the value is to −1, the stronger the correlation. A value close to 0 indicates no spatial correlation.

In the process of constructing the global Moran’s I for the traffic black spots of road intersections in the Kowloon Peninsula region of Hong Kong, the most crucial aspect is the setting of the parameter value of the conceptualization of spatial relationship, which mainly reflects the inherent relationship between the traffic black spots, and the common conceptualizations include the inverse distance, the undifferentiated area, and the connecting edges only, and so on. Based on the characteristics of road intersections, this study decides to adopt the inverse distance method: it shows that the closer the road intersections are spatially, the more likely they are to interact with each other and cause each other to become traffic black spots.

In MapGIS, analysis based on Moran’s index, z-value, and p-value can help us to assess the spatial autocorrelation of geospatial data.

Table 24. p-value, z-value, and confidence level.

z (Standard Deviation)	p (Probability)	Confidence Level
Absolute value of z is greater than 1.65	less than 0.1	90%
Absolute value of z is greater than 1.96	less than 0.05	95%
Absolute value of z is greater than 2.58	less than 0.01	99%

illustrates the relationship between z-value, p-value, and confidence level.

When p is less than 0.1 and the absolute value of Z is greater than 1.65, this indicates that there is a 90% probability that the distribution of traffic black spots at road intersections in the Kowloon Peninsula region of Hong Kong is correlated.

When p is less than 0.05 and the absolute value of Z is greater than 1.96, this means that there is a 95% probability that the distribution of traffic black spots at road intersections in the Kowloon Peninsula region of Hong Kong is correlated.

When p is less than 0.01 and the absolute value of Z is greater than 2.58, this means that there is a 99% probability that there is a correlation between the distribution of traffic black spots at road intersections in the Kowloon Peninsula region of Hong Kong.

Other cases can be recognized as accepting the null hypothesis, i.e., it is a random state.

From

Table 25. Global Moran’s I for traffic black spots.

Moran’s I	Z-Value	p-Value
0.55	38.08%	0.00

, from the results of the construction of the global Moran’s I of traffic black spots at road intersections in the Kowloon Peninsula region of Hong Kong, it can be seen that the global Moran’s I is 0.55, which indicates that there is a positive correlation in the distribution of traffic black spots, and that this correlation may be due to the built-up environments of adjacent road intersections interacting with each other. A p-value of 0 and a z-value of 38.08 indicate that there is a positive correlation in the global Moran’s I, and that it is statistically significant, supporting the existence of spatial autocorrelation in the distribution of traffic black spots at road intersections in the Kowloon Peninsula region of Hong Kong.

The results in Figure 25 also indicate that there is a global correlation between traffic black spots at road intersections in the Kowloon Peninsula region of Hong Kong.

5.3.2. Local Moran’s I

The local Moran’s I is an index used to measure the local spatial autocorrelation of each geographic unit. The local Moran’s I is constructed for the traffic black spots at road intersections in the Kowloon Peninsula region of Hong Kong, and it identifies the spatial aggregation characteristics by comparing the attributes of each road intersection and its neighboring road intersections. Unlike the global Moran’s I, the local Moran’s I can help to identify the geographic units with local aggregation of the traffic black spots at the road intersections in the Kowloon Peninsula region of Hong Kong and reveal the spatial distribution of the traffic black spots at the road intersection patterns and heterogeneity.

The results in

Table 26. Number of distributions of each cluster set of localized local Moran’s I for traffic black spots.

Each Clustering Set	Number
High–high clustering	599
High–low clustering	20
Low–low clustering	819
Low–high clustering	340
Not significant	741

indicate that there are four types of results, such as this is a traffic black spot aggregation area, a non-traffic black spot aggregation area, and so on, in the local distribution of traffic black spots at road intersections in the Kowloon Peninsula region of Hong Kong. Among them, the non-traffic black spot aggregation area is the most significant.

Figure 26 shows the spatial clustering distribution of traffic black spots at road intersections in the Kowloon Peninsula region of Hong Kong. The pink part represents the traffic black spot aggregation area, the light blue part represents the non-traffic black spot aggregation area, the red part represents the area where the traffic black spot aggregation area is surrounded by the non-traffic black spot aggregation area, and the dark blue part represents the area where the non-traffic black spot aggregation area is surrounded by the traffic black spot aggregation area. The non-traffic black spot agglomeration area is mainly concentrated on highways, where there are fewer people and buildings, and this physical environment leads to fewer traffic accidents and less possibility of becoming a traffic black spot, while the high traffic black spot agglomeration area is mainly concentrated in the Yau Tsim Mong District, where there is a dense road network, with a transportation junction, and a dense population and buildings, which have an agglomeration effect and lead to a greater probability of traffic accidents, leading it to become a traffic black spot. The probability of traffic accidents becomes greater, resulting in a higher probability of it becoming a traffic black spot. From the perspective of spatial correlation, it can be understood that the clustering area of a traffic black spot also affects the traffic safety of the surrounding area, resulting in a higher likelihood of traffic accidents.

5.3.3. Summary of Spatial Correlation of Traffic Black Spots

The global Moran index analysis of traffic black spots at road intersections indicates that the distribution of traffic black spots is visibly spatially aggregated and positively correlated in the whole Kowloon Peninsula area. The local Moran index analysis of traffic black spots at road intersections indicates that there are visible traffic black spot aggregation areas and non-traffic black spot aggregation areas, which also indicates that adjacent road intersections will influence each other, increasing or decreasing the likelihood of each other becoming traffic black spots.

The above experiments show that there exists a spatial correlation of traffic black spots at road intersections in the Kowloon Peninsula region of Hong Kong, i.e., adjacent road intersections will affect each other in terms of their probability of becoming traffic black spots. Based on this situation, this paper improves the recognition method of road intersection traffic black spots, constructs a traffic black spot network topology, takes road intersections as nodes, uses the connecting lines between nodes as edges of the network, and the built environment of road intersections as node features, and uses the GCNN graphical convolutional neural network for traffic black spot recognition, which is a better way to perform contextual learning and feature aggregation of neighboring node information, and achieve better traffic black spot recognition.

5.4. GCNN-Based Model Experiment and Result Analysis

(1)

Experimental results

The data are divided randomly and the ratio of data in the training set and test set is 8:1:1. The accuracy of the test set is close to 90% and the Loss value is 0.43. This is shown in Figure 27.

The results of the graphical convolutional neural network for traffic black spot identification at road intersections are visualized in 2D space, as shown in Figure 28.

According to the experimental results of the graph convolutional neural network, it can be found that the GCNN model can be better used in the identification of traffic black spots. Through the extraction of semantic information between the built environment features of adjacent road intersections, the accuracy of the identification of traffic black spots has been improved to a certain extent, and the model value Loss has decreased significantly, which also confirms that, in daily life, the built environment is more similar between the neighboring roads and that, when a traffic accident occurs in a certain place, it brings about the negative effects of vehicular congestion, crowd gathering, and other negative effects which will be transmitted to the neighboring roads, which will make the probability of the occurrence of traffic accidents greatly increased.

6. Discussion

In this paper, based on the street view image data, we use the semantic segmentation model to obtain the built environment features of urban roads, take the road intersection as the research unit, carry out the research idea of traffic black spot identification and its cause analysis, construct a model for identifying traffic black spots based on the built environment of adjacent road intersections, and broaden the research idea of traffic black spot identification. The research results of this paper are summarized as follows: (1) In this paper, the Tencent Street View platform and Python crawler technology are used to obtain the street view image data of Hong Kong’s Kowloon area, and deep learning semantic segmentation technology is used to construct the built environment features of urban road intersections and define the concept of traffic black spots. (2) In this paper, based on the set of built environment features of road intersections, the support vector machine model, gradient boosting model, random forest model, and fully connected neural network model are applied to identify traffic black spots at road intersections, and the experimental results show the feasibility of identifying traffic black spots at road intersections based on the street view image data, and that their experimental accuracy can reach 87%. (3) In this paper, based on the above experiments, considering that the interconnection between the built environment features of adjacent road intersections will have an impact on the identification of traffic black spots, the traffic accident data and street view image data are processed, and the graph convolutional neural network model is used to identify traffic black spots at intersections and achieves better experimental performance. (4) This paper uses principal component analysis to analyze in detail several features that have the greatest impact on the formation of traffic black spots at urban road intersections, such as buildings, billboards, sky, greenery, etc., so that the traffic management department and the urban planning department can regulate the density of building development, pedestrian and vehicular flow, billboards, and greenery according to these findings.

We acknowledge three major limitations of this work and point out three directions for future work. (1) The scope of this paper is limited to the study of only one region, Kowloon Peninsula, Hong Kong, and only road intersections were selected as the study unit, resulting in a small amount of data. In the future, we will continue to expand the size of the whole dataset so that the results will be more reliable and convincing. (2) Secondly, this paper did not consider the time division in the process of identifying traffic black spots. Since most people have different travelling purposes on weekdays and weekends, and the situation of the morning and evening rush hours on weekdays is very different from other times of the day when people and traffic gather, all these factors may affect the distribution of traffic black spots; therefore, in the next stage of this study, the identification of traffic black spots by time period can be more refined to propose corresponding preventive measures to achieve “one point one policy”. (3) This paper is not comprehensive enough in the research and analysis of data. In the process of data collection, the current traffic management department is not able to collect all the information of the accident correctly, there is a part of the data which cannot be recorded by the precise geographic coordinates, and the attributes of the existing data on traffic accidents are not perfect, such as the lack of records on the characteristics of the vehicles and people involved, which will have an impact on the analysis of the causes of the traffic black spots.

7. Conclusions and Future Work

This study focuses on the recognition and causal analysis of black spots in road intersections using street view images. The research methodology revolves around utilizing street view image data as a foundation and employing a semantic segmentation model to extract the characteristics of the urban roads’ built environment. Specifically, this study adopts road intersections as the primary research unit to conduct the identification and cause analysis of traffic black spots. It innovatively constructs a model that identifies traffic black spots based on the built environment of adjacent road intersections, thereby expanding the current understanding of traffic black spot identification. The proposed model and accompanying research offer a novel perspective for traffic management authorities to effectively address traffic safety concerns. By leveraging insights from this study, it becomes possible to mitigate the occurrence of traffic accidents and reduce the probability of developing traffic black spots within the urban built environment.