Next Article in Journal
The Response of Cloud-Precipitation Recycling in China to Global Warming
Next Article in Special Issue
Multi-Task Fusion Deep Learning Model for Short-Term Intersection Operation Performance Forecasting
Previous Article in Journal
Development of Novel Classification Algorithms for Detection of Floating Plastic Debris in Coastal Waterbodies Using Multispectral Sentinel-2 Remote Sensing Imagery
Previous Article in Special Issue
An Analytical Framework for Accurate Traffic Flow Parameter Calculation from UAV Aerial Videos
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Higher-Order Graph Convolutional Network for Location Recommendation of an Air-Quality-Monitoring Station

1
Department of Automation, University of Science and Technology of China Hefei, Hefei 230026, China
2
Hefei Comprehensive National Science Center, Institute of Artificial Intelligence, Hefei 230088, China
3
Institute of Advanced Technology, University of Science and Technology of China, Hefei 230088, China
4
State Key Laboratory of Fire Science, University of Science and Technology of China, Hefei 230027, China
*
Author to whom correspondence should be addressed.
Submission received: 9 February 2021 / Revised: 10 April 2021 / Accepted: 16 April 2021 / Published: 20 April 2021

Abstract

:
The location recommendation of an air-quality-monitoring station is a prerequisite for inferring the air-quality distribution in urban areas. How to use a limited number of monitoring equipment to accurately infer air quality depends on the location of the monitoring equipment. In this paper, our main objective was how to recommend optimal monitoring-station locations based on existing ones to maximize the accuracy of a air-quality inference model for inferring the air-quality distribution of an entire urban area. This task is challenging for the following main reasons: (1) air-quality distribution has spatiotemporal interactions and is affected by many complex external influential factors, such as weather and points of interest (POIs), and (2) how to effectively correlate the air-quality inference model with the monitoring station location recommendation model so that the recommended station can maximize the accuracy of the air-quality inference model. To solve the aforementioned challenges, we formulate the monitoring station location as an urban spatiotemporal graph (USTG) node recommendation problem in which each node represents a region with time-varying air-quality values. We design an effective air-quality inference model-based proposed high-order graph convolution (HGCNInf) that could capture the spatiotemporal interaction of air-quality distribution and could extract external influential factor features. Furthermore, HGCNInf can learn the correlation degree between the nodes in USTG that reflects the spatiotemporal changes in air quality. Based on the correlation degree, we design a greedy algorithm for minimizing information entropy (GMIE) that aims to mark the recommendation priority of unlabeled nodes according to the ability to improve the inference accuracy of HGCNInf through the node incremental learning method. Finally, we recommend the node with the highest priority as the new monitoring station location, which could bring about the greatest accuracy improvement to HGCNInf.

1. Introduction

In recent years, with economic growth, environmental problems have become increasingly prominent and air pollution is receiving unprecedented attention [1,2,3,4]. The air-quality index (AQI) provides a number used by government agencies to communicate to the public how polluted the air is currently. As the AQI increases, an increasingly large percentage of the population is likely to experience increasingly severe adverse health effects [5,6,7]. To compute this AQI, an air-pollutant concentration from a monitor or model is required, such as carbon monoxide (CO), carbon dioxide (CO 2 ), hydrocarbons (HC), nitrogen oxides (NO x ), solid particulate matter (PM 2.5 and PM 10 ), etc. In order to reflect the air quality and its development trends in a timely and accurate manner, we need accurate air-quality-monitoring equipment.
However, it is unrealistic to establish a large number of monitoring equipment in different areas of the city due to the expensive construction and later maintenance costs [8]. Based on the premise of a limited number of monitoring equipment, a reasonable monitoring station location layout is the first basic link to accurately infer air quality. On the other hand, with the rapid development of urbanization and the continuous growth in the number of motor vehicles, the number of established monitoring equipment is too few to accurately reflect the air-quality distribution in an entire urban area.
Thus far, there is few research on the optimization methods of regional air-quality-monitoring networks and there is also no uniform standard for recommending optimal monitoring equipment locations. Therefore, we need to design a framework that could recommend optimal station locations to establish new monitoring equipment on the basis of the existing ones, which can bring about the greatest accuracy improvement to the inference model. In this way, it not only meets the accuracy requirements for air-quality inference but also saves economic costs.
In recent years, the research on the location selection is mainly divided into the following two categories: knowledge-driven approaches and data-driven approaches.
The knowledge-driven approaches use mathematical models [9,10] and physics knowledge [11] to solve the location recommendation problem through computational simulation. In addition, the United States, Europe, and Japan have also established regional air-quality-monitoring networks around photochemical smog pollution and aerosol pollution [12,13,14]. These methods all adopted a series of optimization techniques to determine the layout of network sites, including statistical analysis methods such as correlation analysis and cluster analysis [15,16] or mathematical methods such as multi-objective optimization [17,18]. In order to reach a stable state, the simulation process not only requires complex system programming but also consumes a lot of computing power. Simplifications and stationarity assumptions that may be unrealistic in modeling further degrade the model efficiency.
Recently, data-driven approaches for location recommendation have been developed based on fixed monitoring stations [19] or Taxi GPS trajectory [20,21] but without considering the road network spatial structure, such as road segment length, road types (highways, main roads, and streets), POI density, etc. Kang et al. [22] transformed the station recommendation into a graph problem, which is committed to covering the urban area with the least number of monitoring stations but without considering the influences of complex external factors, which might lead to geographically non-smooth values in air-quality distribution. Hsieh et al. [23] studied label propagation [24] on graphs and considered external influential factors but failed to capture the spatiotemporal correlation between nodes in graph.
In this paper, we want to solve a practical problem: how to recommend optimal station locations to establish air-quality-monitoring equipment based on existing ones to maximize the accuracy of the air-quality-inference model, so as to reflect the air-quality distribution and its development trend in a timely manner and accurately. This task is challenging for the following reasons:
(1) The air-quality distribution is affected by many complex external influential factors (such as weather, traffic volume, land use, etc). For commercial centers with heavy traffic flow, the air quality is often worse than that in parks or lakes. There will be geographical non-smooth data due to the influence of external influential factors, and it is difficult to obtain accurate air-quality distribution in unobserved areas through interpolation-based methods. For example, Figure 1 is a real-time record of a certain day at the Beijing Air-Quality-Monitoring Station. We can find that the PM 10 concentration data are not smooth. The monitoring stations in the red circle are very geographically close, but the PM 10 concentration varies greatly throughout the year. The most likely reason for this result is that the monitoring stations with low PM 10 concentration are close to parks or lakes and that the stations with high PM 10 concentration are located near commercial centers or main roads with heavy traffic.
(2) There is spatiotemporal interaction between different nodes in the urban spatiotemporal graph: At the same time, the air-quality value of a node is affected by other nodes; in the same space, the air-quality values of the same node at different times are correlated.
(3) Complex correlation between the monitoring station location recommendation model and the air-quality inference model. The task of the monitoring station location recommendation model is to select optimal locations to establish new monitoring stations to maximize the accuracy of an air-quality-inference model for inferring the air-quality distribution of an entire urban area. There is a close relationship between two models.
Aiming to solve the main research problems of this paper, the current research methods either fail to consider the urban spatial structure or the influence of external factors on air-quality distribution; the first problem is how to accurately infer air-quality distribution. However, the performance of an air-quality-inference model has a crucial impact on the subsequent monitoring-station location recommendation. The next most direct question is how to associate the monitoring-station location recommendation model with the air-quality-inference model to recommend optimal station locations. The existing methods cannot capture the correlation between nodes because of the performance limitation of the model, so the recommended station location is not optimal.
In order to solve the aforementioned challenges, we formulate monitoring station locations as an urban spatiotemporal graph (USTG) node recommendation problem in which each node represents a region with time-varying air-quality values and propose a two-step learning framework.
The main contributions of this research are as follows:
(1) We propose a variant of graph convolutional network (GCN) called higher-order graph convolutional network (HGCN). The ordinary GCN can capture the one-hop neighbor node’s spatiotemporal correlation. We improve the ordinary graph convolution network and design a high-order graph convolution network. Compared with an ordinary (one-hop) GCN, HGCN has a larger receptive field, which could capture more and higher-order neighbor node information.
(2) We designed an accurate air-quality-inference model (HGCNInf) based on the proposed HGCN to infer urban air-quality distribution. We applied the graph convolutional network to the field of air-quality monitoring for the first time. By modeling the urban area as an USTG, the HGCNInf can accurately infer the air quality of an entire urban area: the HGCN can effectively capture the spatiotemporal interactions of the air-quality distribution, and a fully connected neural network is used to extract the external influential factor features.
(3) We analyze and correlate the air-quality-inference model with the monitoring-station location recommendation model. By using the convolution weight parameters of HGCNInf, we design a GMIE based on the correlation degree between nodes in USTG reflecting the air-quality spatiotemporal changes, marking the recommendation priority of unlabeled nodes according to the ability to improve the inference accuracy of HGCNInf iteratively through the node incremental learning method. The recommended optimal node could bring about the greatest accuracy improvement to HGCNInf.
(4) We evaluate our model using Beijing air-quality data from 1 January 2015 to 31 December 2015, and the experimental results show that our approach far outperforms the state-of-the-art baseline methods.

2. Study Area and Data Requirement

2.1. Study Area

Beijing is China’s capital city (Figure 2), located at 40.250583 N and 116.462611 E, on the northwestern edge of the North China Plain. The average altitude of Beijing is 43.5 m above sea level. Beijing is located in a warm temperate and semi-humid area. The climate is a warm temperate and characterized by semi-humid continental monsoons. The average annual precipitation in the plain area is about 600 mm, and the annual average temperature is about 12.9 C.
The main sources of air pollution in Beijing are industrial exhaust gas and automobile exhaust emissions. Our study area covers the Haidian, Shijingshan, Fengtai, Xicheng, and Chaoyang districts, which are the main urban districts in Beijing. The air-quality data collected from the above areas can reflect the overall air-quality condition in Beijing.

2.2. Data Requirement

We evaluated our proposed model with Beijing air-quality data, which comes from the Beijing Municipal Ecological and Environment Monitoring Center. Air-quality data are collected every hour, and the time span is from 1 January 2015 to 31 December 2015. The data are from 17 monitoring stations. and a total of 17 × 8760 data points were collected. In the experiment, we divided the main urban area of Beijing into 30 × 38 sub-regions (1 km × 1 km), in which 17 sub-regions have established monitoring stations, while the air-quality distribution of the remaining 1123 sub-regions are unknown. The external factors include urban meteorological data, the number of points of interest (POIs) in the sub-region, and the geographic road network information in the sub-region (such as the length of highways and main roads, the number of intersections, etc.) [25]. These data features are represented as F = f 1 , f 2 , , f m , where f m represents the mth data feature. The types of POI are shown in Table 1.
Based on urban geographic information, the urban area is divided into disjoint grids. Each grid is a sub-area of 1 km × 1 km and is associated with time-varying air quality values. As shown in Figure 3.

3. Preliminary and Methodology

3.1. Preliminary

3.1.1. Air-Quality Index

The air-quality index (AQI) is a number used by government agencies to communicate to the public how polluted the air is currently. As the AQI increases, an increasingly large percentage of the population is likely to experience increasingly severe adverse health effects. To compute the AQI, an air-pollutant concentration from a monitor or model is required. The function used to convert the air pollutant concentration to AQI varies by pollutants and is different in different countries. The air-quality index values were divided into ranges, and each range was assigned a descriptor and a color code.

3.1.2. Graph Convolutional Networks

Convolutional Neural Networks (CNNs) are mainly effective on Euclidean domain data with regular spatial structures. Generalizing convolutional neural networks on arbitrary graph-structured data is a hot topic in recent researches. Bruna et al. [26] generalized the convolution operator to a general graph via the graph-Laplacian spectrum. By transforming both the filter and the signal to the Fourier domain, convolutions on the graphs can be calculated by multiplying the filter and the signal and then by transforming the result back into the discrete domain. As a meaningful translation operator cannot be expressed in the node domain [26], Defferrard et al. [27] defined the convolution operator G on graph in the Fourier domain. A signal x filtered by g w = d i a g ( w ) is as follows:
g w G x = g w ( L ) x = g w U Λ U τ x = U g w ( Λ ) U τ x .
where U is a matrix of eigenvectors of the normalized graph Laplace matrix L and Λ is a diagonal matrix of eigenvalues of L.
However, evaluating Equation (1) is computationally expensive, as its computational complexity is O N 2 . To circumvent this problem, the filter g w ( Λ ) is well approximated by a truncated expansion of K-order Chebyshev polynomials suggested in [28]. Based on [28], Defferrard et al. [27] designed fast localized convolutional filters on graphs to reduce the computational complexity.
g w G x = g w ( L ) x U k = 0 K 1 W k T k ( A ˜ ) U T x = k = 0 K 1 W k T k ( L ˜ ) x
where T k ( L ˜ ) is a K-order polynomial about L ˜ , L ˜ = 2 λ max L I N , a n d λ max is the largest eigenvalue of the normalized graph Laplacian matrix L. The computational complexity is reduced from O N 2 in Equation (2) to O N .
Furthermore, Kipf et al. [29] designed a multi-layer graph convolution network for semi-supervised learning via a first-order approximation of localized spectral filters on graphs. The definition is as follows:
Z = D ˜ 1 2 A ˜ D ˜ 1 2 X W
where Z is the convolved signal matrix, A ˜ = A + I N is the graph adjacency matrix with added self-connections, D ˜ i i = Σ j A ˜ i j , and W is a trainable weight matrix. The computational complexity of this convolution operation is the same as A ˜ X .
As a typical type of graph, based on the spectral graph convolutions in [29], Yu et al. [30] designed a spatiotemporal graph convolutional network for traffic forecasting. Xu et al. [31] designed a spatiotemporal graph multi-fusion network for urban mobile source emission prediction. Sun et al. [32] designed multi-view graph convolutional networks to predict citywide crowd flows in irregular regions. Especially, Cui et al. [33] designed a traffic higher-order graph convolutional recurrent neural network for network-scale traffic learning and forecasting.

3.1.3. Urban Spatiotemporal Graph

We take each grid element as a node and to construct an urban spatiotemporal graph G = ( V ; ε ; A ) (USTG), where V is the set of nodes associated with the time-varying air-quality value; V = N denotes the number of grid nodes; ε is the set of edges, indicating the connectivity between the nodes; and A R N × N is the graph adjacency matrix. In particular, each vertex v i V in the USTG is a three-dimensional representation; the first two dimensions represent the geographic location information of the node, and the third dimension represents the air-quality value that changes with time. In USTG, V can be represented as V = L U , where L represents the node set with air-quality-monitoring stations and U represents the node set without monitoring stations. We take the nodes in L as labeled nodes and the nodes in U as unlabeled nodes. Generally speaking, L     U . The spatial correlation of the spatiotemporal graph is reflected by the geographical distance characteristics between the nodes, and the temporal correlation is reflected by the air-quality value of different times associated with the nodes.

3.1.4. Inference Problem of Air Quality

Problem 1.
At time T i i = 1 , 2 , , n , based on USTG G = ( V ; ε ; A ) , the air-quality distribution f l R | L | × P of the labeled nodes L is used to infer the air-pollution distribution f u R | U | × P of the unlabeled nodes U. The task of air-quality inference is to learn a function F ( · ) that maps f l R | L | × P to f u * R | U | × P :
F f l ( T i ) , G ( V ; ε ; A ) = f u * ( T i ) .
where | L | is the number of labeled nodes, | U | is the number of unlabeled nodes, and P is the number of features associated with each node after one-hot encoding processing.

3.1.5. Location Recommendation for the Monitoring Station

Problem 2.
The main objective of monitoring equipment site location recommendation is to choose optimal practical locations to establish new monitoring equipment at candidate locations. The optimal location can bring about the greatest accuracy improvement to the inference model. There is a relatively close relationship between the monitoring station recommendation model and the air-quality-inference model. We need to analyze and correlate these two models and to design an effective algorithm based on a correlation degree to complete the monitoring-station location recommendation.

3.2. Methodology

In this section, we introduce the air-quality-inference model and the monitoring-station location recommendation model. We first introduce the conception of higher-order graph convolutional networks. Then, we design an air-quality-inference model based on the improved HGCN, which is applied effectively to the field of atmospheric environment monitoring. Next, we analyze the air-quality inference model and correlate it with the air-quality-monitoring station location recommendation. We design a greedy algorithm for minimizing information entropy based on the correlation degree between nodes in USTG, marking the recommendation priorities of unlabeled nodes according to the ability to improve the inference model accuracies using the incremental learning method to complete the monitoring-equipment station location recommendation.
The proposed two-step air-quality-inference learning framework is shown in Figure 4.
(1) Air-quality inference model. We propose an effective air-quality inference model based on the proposed higher-order graph convolution network (HGCNInf) to infer the air-quality distribution, consisting of a higher-order graph convolution network (HGCN) and several fully connected neural networks. By inputing historical monitoring data and by fetching global information, the HGCN is used to capture the spatiotemporal interactions of air-quality distribution and the fully connected neural network is used to extract external influential factor features, so that HGCNInf can accurately infer the air-quality distribution of an entire urban area.
(2) Monitoring-station location recommendation model. HGCNInf takes the spatiotemporal interactions of pollutant distribution and external influential factors into account, in which the model convolution weight parameters actually represent the correlation degree between the nodes in USTG that reflect the spatiotemporal changes in air quality. The high correlation between a node and other nodes indicates that the air quality of the node is easier to infer. Therefore, for the station location recommendation, it should be chosen in such a node that has a low correlation with other nodes. In this paper, we introduce the information entropy index to evaluate the correlation degree between nodes. The node with the lower information entropy indicates that it has a higher correlation degree with the remaining nodes. Next, we designed a greedy algorithm for minimizing information entropy (GMIE) based on the correlation degree between nodes in USTG, marking the recommendation priorities of unlabeled nodes according to the ability to improve the inference accuracies of HGCNInf through the node incremental learning method to complete the monitoring station location recommendation.

3.2.1. Higher-Order Graph Convolutional Network

Some recent researches [34,35,36] have shown that, for graph node semi-supervised learning, random walk statistics can be used to effectively learn feature representations because it could preserve the completeness of the graph structure. Under special conditions, if the activation function of the first layer in the graph convolutional network is an identity function, then the multi-layer GCN model can also learn random walks. For the 2-layer GCN model proposed in [29], we assume that the activation function of the first layer is an identity function:
σ A ^ X W ( 0 ) = A ^ X W ( 0 )
In this way, the GCN model can be further simplified to the following form:
Z = f ( X , A ) = sof tmax A ^ A ^ X W ( 0 ) W ( 1 ) = softmax A ^ 2 X W ( 0 ) W ( 1 ) = softmax A ^ 2 X W ( 2 )
where W ( 2 ) = W ( 0 ) W ( 1 ) is expressed as the weight parameter of A ^ 2 , A ^ = D ˜ 1 2 A ˜ D ˜ 1 2 . The two-layer one-hop graph convolutional network is converted into a one-layer two-hop convolutional network. Furthermore, A ^ 2 can be expressed as follows:
A ^ 2 = D ˜ 1 2 A ˜ D ˜ 1 2 D ˜ 1 2 A ˜ D ˜ 1 2 = D ˜ 1 2 A ˜ D ˜ 1 A ˜ D ˜ 1 2 D ˜ 1 A ˜
We can find that the 2-hop GCN network under certain conditions actually learns a one-step random walk, which can transfer the information of the node to the higher-order neighbor nodes. For the 2-layer GCN model proposed in [29], the flow of the input signal through the graph edges is reduced because the information features have gone through matrix multiplication (average neighbor node information) and a nonlinear activation function σ ( · ) . Therefore, it is not possible to directly take higher-order operations on A ^ based on the original multi-layer GCN network; otherwise, it causes network performance degradation. In order to make better use of the information of high-hop neighbor nodes, we increase the order of the graph convolution network while reducing the layer numbers of the network and then fuse the features extracted from each order network.
A ^ k = D ˜ 1 2 A ˜ D ˜ 1 A ˜ ( k 1 ) D ˜ 1 2
G C k = GCN A ^ k , X = A ^ k X W ( k )
Z = f ( X , A ) = G C { k } = softmax G C 1 , G C 2 , , G C k
Thus, the newly designed higher-order graph convolutional network uses fewer layers and parameters but can achieve better results.

3.2.2. The Air-Quality Graph Convolution

We first define the graph higher-order adjacency matrix.The one-hop neighbor matrix of the USTG is the adjacency matrix A, and the k-hop adjacency matrix can be obtained by the kth power of A. As shown in Figure 5, the multi-hop neighbor with respect to a node (red star) is the nodes that can be reached by the black edges in the graph.
We apply the proposed higher-order graph convolutional network to the field of air-quality monitoring. The problem of air-quality inference is specifically described in Section 3.1.4 and the proposed HGCNInf framework is shown in Figure 6.
In order to make the subsequent air-quality monitoring station location recommendation model more reasonable and interpretable, we modified the graph convolution network G C k described in Equation (9) to the following form.
G C k = GCN A ^ k , X = A ^ k W ( k ) X
where A ˜ = D 1 2 A D 1 2 R N × N , D is the degree matrix of A, and W ( k ) is the convolution weight parameter. The graph input signal X = f l f u R N × P , where f l R L × P . is the air-quality distribution of the labeled nodes, initializing the air-quality distribution of unlabeled nodes to f u = 0 0 0 0 R | U | × P .
Then, we concatenated the K-hop features extracted from the graph signal, defined as follows:
G C { k } = G C 1 , G C 2 , , G C k
where G C { k } R N × ( K × P ) is a set of K-hop features extracted by graph convolutional network. We fed G C { k } into a fully connected network f f c 1 to jointly train
f f c 1 : R N × ( K × P ) R N × P
G C = f f c 1 G C { k } , W f c 1
where air-distribution feature G C R N × P and the trainable parameter matrix W f c 1 R P × ( K × P ) .
We use the fully connected network f f c 2 to extract external factor features, defined as follows:
f f c 2 : R N × M R N × P
G C f = f f c 2 f m , W f c 2
where the trainable parameter matrix W f c 2 R P × M ; external factor features G C f R N × P ; and f m R N × M , which is the M external factors.
We concatenated the air-quality distribution feature G C R N × P with the external factor feature G C f R N × P and then fed them into the fully connected network f f c 3 , defined as follows:
f f c 3 : R N × P R N × P
GC a l l = f f c 3 GC + G C f , W f c 3
where the trainable parameter matrix W f c 3 R P × P and general feature G C a l l R N × P .
Further, the air quality distribution F u R | U | × P of unlabeled nodes is obtained, defined as follows:
F = softmax GC a l l
F = F l F u R N × P
where F l R | L | × P is the air-quality distribution output of the labeled nodes and F u R | U | × P is the air-quality distribution output of the unlabeled node.
We minimized the cross-entropy loss of the output F l and the known training label Y. The cross-entropy loss is defined as follows:
H Y , F l = diag V l Y log F l
where ∗ represents the Hadamard product and where d i a g ( V l ) is a diagonal matrix where elements are set to 1 when ( i , i ) V l and 0 otherwise. In order to make the graph convolution features more stable and interpretable, we added regularization terms to the loss function.
(1) We added the L1-norm of the graph convolution weight matrix to the loss function as a regular term to make these weight matrices as sparse as possible. L1 regularization is defined as follows:
R 1 = W g c i 1 = i = 1 k W g c i
(2) In order to limit the differences between adjacent features extracted from graph convolution, the convolution feature regularization term based on the L2 norm is added to the loss function. L2 regularization is defined as follows:
R 2 = G C { k } 2 = i = 1 k 1 G C i G C i + 1 2
Thus, the total loss function is defined as follows:
Loss = H Y , F l + α · R 1 + β · R 2
where the parameters α and β are used to control the weight of the regularization term.
Next, we minimized the model loss function to train the network as follows:
min Loss = min [ diag V l Y log F l + α i = 1 k W g c i + β i = 1 k 1 G C i G C i + 1 2 ]
When the network training was completed, we found the value with the highest quantization probability from the air-quality distribution F u and defined f u * as the inferred air-quality value. The definition is as follows:
f u = argmax F u

3.2.3. Air-Quality-Monitoring Station Location Recommendation

The main objective of this section is use the proposed HGCNInf to complete the location recommendation for a air-quality-monitoring station. The ultimate goal is to recommend an optimal station location from candidate locations (this paper defaults to all unlabeled nodes) that could bring about the greatest accuracy improvement to the inference model. Therefore, we need to correlate the air-quality-inference model with the monitoring-station location recommendation model. For the trained higher-order graph convolutional network, its weight parameters reflect the correlation degree between nodes. Since the designed HGCNInf takes into account both the spatiotemporal interactions of air pollutant distribution and the complex external factors, its convolution weight parameter represents the correlation degree between the nodes reflecting the spatiotemporal changes in air-quality distribution. This is of great significance to the monitoring-equipment location selection.
In the USTG, a higher correlation degree between a node and the remaining nodes indicate that the air quality of the node is easily inferred. Therefore, the monitoring station should be selected at a node that has the lowest correlation with the remaining nodes. We introduced the information entropy index to evaluate the correlation degree between nodes. The smaller the information entropy value, the higher the correlation degree between nodes. We used the convolution weight parameters of the trained HGCNInf to recommend optimal monitoring station locations that could also bring about the greatest accuracy improvement to HGCNInf. In this way, the two problems are well correlated.
Greedy algorithm based on information entropy minimization
We designed a greedy algorithm for minimizing information entropy (GMIE) based on the correlation degree between nodes in USTG to recommend optimal station locations. HGCNInf is used as the driving model through the incremental learning method, iteratively marking the recommendation priority of unlabeled nodes.
The basic framework is shown in Figure 7, and the specific steps are as follows from step 1 to step 6:
1. Calculate the average convolution weight parameter of the trained HGCNInf. The average convolution weight parameter is defined as follows:
W ¯ = 1 k i = 1 k A ˜ i W g c i
where k is the order of the higher-order graph convolutional network and A ˜ is the standardized adjacency matrix. Next, calculate the information entropy of the average convolution weight parameter of each node in the unlabeled node U, which is defined as follows:
H W ¯ u = u U W ¯ u log W ¯ u
2. Find the unlabeled node u U with the smallest information entropy value in H W ¯ u , and mark its recommendation priority as the lowest. Based on the air-quality distribution f u inferred from HGCNInf, associate the inferred air-quality value with the node u and then turn the node u from the unlabeled nodes into the labeled nodes to constitute a new graph input signal data f l . Remove the node u U from the unlabeled node U, and the remaining unlabeled nodes are represented as U .
3. Based on the new graph input data f l , drive HGCNInf to retrain and update the network convolution weight parameters, that is, update the correlation degree between the nodes in USTG reflecting the spatiotemporal changes in air-quality distribution and then derive the new air-quality distribution f u of unlabeled nodes. As described in step 1, new average convolution weight parameters W ¯ and information entropy H W ¯ u can be calculated.
4. Find the unlabeled node u U with the smallest information entropy value in H W ¯ u , and mark its priority as the second lowest. Based on the new air-quality distribution f u inferred from the air-quality inference model, associate the inferred air-quality value with the node u , and then, turn the node u from the unlabeled nodes into the labeled nodes to constitute a new graph input signal data f l . Remove the node u U from the unlabeled node U , and the remaining unlabeled nodes are represented as U .
5. Repeat step 1 to step 4 above until the recommendation priority of all unlabeled nodes U is marked.
6. Furthermore, input graph signal data at other moments and repeat step 1 to step 5.
Finally, analyze the results at multiple times and average the recommendation priority of unlabeled nodes, and then recommend the Top-N nodes with the highest average recommendation priority as the new monitoring station location. These Top-N sites with the highest average recommendation priority have the highest average information entropy, that is, the lowest correlation degree with the remaining nodes. Therefore, newly established monitoring equipment at these Top-N sites could bring about the greatest accuracy improvement to HGCNInf, so that it can reflect the air quality and its development trends in a timely and accurate manner.

4. Experiments and Results

4.1. Data Processing

The graph input data in this paper isare the PM 10 values of Beijing and the complex external influential factor data. We used one-hot encoding to process PM 10 data and used Z-Score Normalization to process the external influential factor data, the factor data were linearly transformed to be mapped into the interval [0, 1].
The adjacency matrix of the Beijing urban spatiotemporal graph was calculated by using the threshold Gaussian kernel weighting function based on the distance between nodes [30,37]. The calculation formula is as follows:
W i j = exp dist p i , p j 2 σ 2 , exp dist p i , p j 2 σ 2 ε 0 , otherwise
where dist p i , p j 2 represents the Euclidean distance between the nodes p i and p j , and σ 2 and ε are the thresholds used to control the distribution and sparsity of the adjacency matrix. In this experiment, σ 2 and ε are set to 10 and 0.05 respectively.

4.2. Experimental Settings

4.2.1. Evaluation

In order to evaluate and compare the performance between different models, we used two commonly used metrics: (1) Mean Absolute Percentage Error (MAPE); (2) Root Mean Squared Error (RMSE). They are defined as follows:
MAPE = 1 n t = 1 n Y i Y ^ i Y i   100 %
RMSE = 1 n i = 1 N ( Y i Y ^ i ) 2
where Y ^ i and Y i are the predicted values and the real values from selecting labeled nodes for model validation, and n is the number of all predicted nodes.

4.2.2. Baselines

We compared HGCNInf and GMIE with the following baseline models.
(1) HA: Historical average. HA takes the average value of air-quality data as the predicted value.
(2) SVR: Support vector regression [38]. SVR uses the historical AQI data of each monitoring station as training data to infer the air-quality value of unlabeled nodes.
(3) AQInf: Affinity-based AQI inference model [23]. For air-quality inference problem, this is the current state-of-the-art method. It is used to infer the air-quality distribution of the unlabeled nodes by using the AQI data of the labeled ones and by combining the external influential factors.
(4) GCNInf: An air-quality inference model based on graph convolutional networks (GCNInf). We built a two-layer graph convolutional networks as described in [29]. Similar to AQInf, GCNInf was used to infer the air-quality distribution of unlabeled nodes by using the AQI data of the labeled ones and by combining with the external influential factors.
All of the neural network model implementations are based on Pytorch version: 1.1.0, and evaluation and training are based on CPU Intel(R) Core(TM) i5-4210U CPU at 1.70GHz, NVIDIA GeForce 840M.

4.3. Experimental Results

4.3.1. The Effectiveness of HCNInf

In the experiment, since we only have real data of the labeled nodes (the sub-regions with monitoring stations), for the unlabeled nodes (the sub-regions without monitoring stations), there are no real data to verify the inference accuracy of each model. Therefore, we proposed to conduct cross validation by randomly selecting 12 nodes as labeled nodes from the 17 labeled ones, which generated 12 × 8760 training data, and the remaining 5 nodes were used to verify the inference accuracy of each model, which generated 5 × 8760 verification data.
(1) The performance comparison of air-quality-inference models
We trained the above air-quality-inference models separately three times and then took the average values as the final results. Table 2 and Figure 8 show the performance comparison of our proposed HGCNInf and other baseline models on the inferred Beijing air-quality distribution in two metrics. It can be found that HGCNInf is far superior to other baseline models both in the metric MAPE and RMSE.
In addition, comparing the results of the air-quality-inference models with external factors and the models without external factors, we can find that, for the same model, after considering complex external influential factors, the error in the air-quality inference decreased significantly. Therefore, the external influential factors play a significant role in improving the inference accuracy of air-quality distribution. We can conclude that the complex external influential factors are closely related to the air-quality distribution, which also explains why the monitoring stations that are geographically close to each other in the red circle differ greatly in the monitoring values throughout the year.
Furthermore, to test the robustness of each air-quality-inference model, we randomly discarded the same labeled nodes and observed the change degree in the air-quality-inference performance. Table 3 and Figure 9 show the inference results of HGCNInf and other baseline models on two evaluation metrics when we gradually randomly discarded the same labeled nodes. From the above chart, we can find that our proposed air-quality-inference model HGCNInf is far more robust than other baseline models.
(2) Training efficiency of HGCNInf We considered the influences of different order K of the graph convolution in HGCNInf on inference performance. As shown in the histogram of Figure 10a,b and Table 4, they showed the performance comparison of different orders of HGCNInf. For Beijing’s air-quality-distribution data, when we chose the order K of the higher-order convolutional network to be 3, the model HGCNInf had the best performance. Figure 8c and Table 5 show the average training times per 1000 times at each moment when the HGCNInf model takes different orders. The training time increases as the order of the higher-order convolutional network increases. In addition, when the HGCNInf model considers the external influential factors, the time consumed by training increases significantly. Comprehensively considering the model performance and training efficiency, we chose K = 3 as the order of the high-order convolutional network in the HGCNInf model, and the results are shown in Table 2.
(3) Weight visualization
To better understand the HGCNInf, we visualized a proportion of the average graph convolution weight W, graph adjacency matrix A, and their product of A*W, as shown in Figure 11.

4.3.2. The Effectiveness of GMIE

The main objective of this section is to recommend several optimal locations from the candidate locations to establish new monitoring stations. In this experiment, we verify whether the nodes recommended by our proposed monitoring-station location recommendation model (GMIE) can maximize the inference accuracy of HGCNInf for inferring the air-quality distribution of the entire city. Since there is no ground real air-quality data for unlabeled nodes, it is impossible to verify the specific accuracy improvement that the recommended nodes bring to the air-quality-inference models.
We adopt a cross-validation method to randomly select 7 nodes from 17 labeled nodes as training datas, and then, we reserve 5 labeled nodes as candidate sites to establish new monitoring equipment. The remaining 5 labeled nodes are used to verify the performance improvement that the recommended nodes bring to the air-quality-inference models.
Baselines. Since to our knowledge, AQInf+GEM is currently the state-of-the-art method for monitoring-station location recommendation [18], we compared our proposed HGCNInf+GMIE with GCNInf+GMIE and AQInf+GEM. The recommendation priorities of the candidate locations marked by each monitoring-station location recommendation model are shown in the Figure 12.
We use MAPE improvement and RMSE improvement to evaluate the effectiveness of monitoring-station location recommendation models. In theory, as the number of recommended station locations increases, the accuracy of the corresponding air-quality inference model should be greatly improved because more station locations are added into the labeled nodes. That is, MAPE and RMSE should be greatly reduced. Since the candidate locations provide the recommendation priority, the node with a higher recommendation priority can bring about greater improvement to the air-quality-inference model. Therefore, when gradually increasing the number of recommended nodes according to the recommendation priority, the declining slopes of MAPE and RMSE should become smaller and smaller or the rising slopes of MAPE improvement and RMSE improvement should become smaller and smaller.
We display how much both MAPE improvement and RMSE improvement models can be obtained by gradually adding more recommended station locations shown in Table 6, and Figure 13 and Figure 14. We can find that the performance of each station location recommendation model is consistent with the theory and that our proposed model’s performance has far exceeded the state-of-the-art method, which can bring about the greatest accuracy improvement to HGCNInf for inferring the air quality.

5. Discussion

5.1. Air-Quality Inference

The knowledge-driven method is based on the formation mechanism of air pollutants and the actual physical diffusion and migration laws, and then a suitable model is established, and the air-quality distribution in urban areas not covered by the monitoring station is predicted through calculation and simulation. Gibson et al. [39] used a Gaussian plume diffusion model to predict the PM 2.5 distribution of main points and routes. Rakowska et al. [40] constructed an urban street canyon model to monitor the impact of traffic flow on air quality. These methods are based on a complex mechanism knowledge, and it is difficult to establish a universal model for various environments due to the adopted empirical assumptions. For example, the Gaussian plume diffusion model that is widely used in the field of air pollution assumes that the distribution of air pollutants conform to the Gaussian distribution. However, in an urban street environment, the diffusion of mobile source pollutants on the road is affected by the buildings on both sides of the road, which does not satisfy this model.
Data-driven methods are also constantly developed and have achieved good results in fine-grained air-quality prediction in space. Zheng et al. [8] considered the influence of external factors related to air-quality distribution and proposed a semi-supervised learning method based on the co-training collaborative training framework [41] to infer the fine-grained air-quality distribution in urban space. However, the framework consists of two independent classifiers, and it is difficult to capture the complex interactions between different features. Chen et al. [42] used the K-nearest neighbor strategy to construct the relationship between external data sources (such as road network structure, traffic, etc.) and air-quality distribution. However, there is no exact theoretical standard for the specific value of K. The characteristics of the latest K sites may not be the most effective, and the degree of influence of different sites may change over time, which may lead to inconsistency problems. Hsieh et al. [23] designed an affinity-based air-quality-inference model AQInf (the affinity-based AQI inference model (AQInf)) based on graph label node propagation. However, due to the limitation of its own model performance, it is difficult for AQInf to capture the correlation between graph nodes that reflect the trend for air-quality temporal and spatial changes.
Our proposed HGCNInf consists of a higher-order graph convolution network (HGCN) and several fully connected neural networks (FNN). By inputing historical monitoring data and by fetching global information, the HGCN is used to capture the spatiotemporal interactions of air-quality distribution and the fully connected neural network is used to extract external influential factor features, so that HGCNInf could accurately infer the air-quality distribution of an entire urban area. As shown in Figure 8 and Figure 9, the prediction accuracy and robustness of HGCNInf are far better than those of each benchmark model. The experimental results show that our proposed model is also applicable even when there are very few label nodes available. In addition, the difference in the two evaluation indicators between the air quality inference model without external influencing factors and the external influencing factors added was compared. The same model after considering the influence of external factors improves the accuracy of urban air-quality distribution prediction and has a significant effect. Therefore, we can conclude that complex external influencing factors are closely related to the distribution of air quality, which also explains the phenomenon why the monitoring sites that are very close geographically in Figure 1 have large differences in values throughout the year.

5.2. Air-Quality-Monitoring Station Location Recommendation

HGCNInf takes the spatiotemporal interactions of pollutant distribution and external influential factors into account, for which the model convolution weight parameters actually represent the correlation degree between the nodes in USTG that reflects the spatiotemporal changes in air quality. Based on the correlation degree, we designed a GMIE that aimed to mark the recommendation priority of unlabeled nodes according to the ability to improve the inference accuracy of HGCNInf to complete the monitoring-station location recommendation. According to Figure 13 and Figure 14, our proposed model’s performance far exceeded the state-of-the-art method, which can bring about the greatest accuracy improvement to HGCNInf for inferring the air quality.

6. Conclusions

In this paper, we formulated the monitoring station location as a urban spatiotemporal graph (USTG) node recommendation problem. We designed an effective air-quality inference model HGCNInf based on a proposed higher-order graph convolution network for inferring the air-quality distribution throughout the city. HGCNInf can not only capture the spatiotemporal interactions of atmospheric pollutant distribution but also extract complex external influential factor features, so that HGCNInf could accurately infer the air-quality distribution of the entire urban area. In addition, HGCNInf can also learn the correlation between the nodes in USTG reflecting the air-quality distribution. Then, we designed a greedy algorithm for minimizing information entropy (GMIE) based on the correlation degree, marking the recommendation priority of unlabeled nodes according to the ability to improve the inference accuracy of HGCNInf to complete the task of recommending the optimal station location. The final recommended optimal station can bring HGCNInf the greatest accuracy improvement. By evaluating real-world Beijing air quality data, HGCNInf and GMIE far outperform the baselines, including the state-of-the-art method.
In the future, we will consider transfer learning methods to improve the applicability of the model to different regions and apply the proposed model to other areas of urban computing, such as traffic flow monitoring, urban regional population flow prediction, etc.

Author Contributions

Conceptualization, Y.K. and J.C.; methodology, Y.K. and J.C.; software, J.C.; validation, Y.C.; formal analysis, J.C. and Z.X.; investigation, J.C.; resources, Y.C. and Z.X.; data curation, J.C.; writing—original draft preparation, Y.K.; writing—review and editing, Z.X.; visualization, J.C.; supervision, Y.C. and Z.X.; project administration, Y.C.; funding acquisition, Y.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by the National Key R&D Program of China under grants 2018YFE0106800, 2018AAA0100800, and 2018YFC0213104; National Natural Science Foundation of China (62033012, 61725304, 61673361, 61773360, and 61872327); Major Special Science and Technology Project of Anhui, China (912198698036, 012223665049); Special Fund for Transformation of Scientific and Technological Achievements of Smart City Research Institute (Wuhu) of University of Science and Technology of China(2019ZX004); as well as the Fundamental Research Funds for the Central Universities under grant (WK2100000013, WK2380000001).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Li, W.; Shao, L.; Wang, W.; Li, H.; Wang, X.; Li, Y.; Li, W.; Jones, T.; Zhang, D. Air quality improvement in response to intensified control strategies in Beijing during 2013–2019. Sci. Total. Environ. 2020, 744, 140776. [Google Scholar] [CrossRef]
  2. Shi, K.; Wu, L. Forecasting air quality considering the socio-economic development in Xingtai. Sustain. Cities Soc. 2020, 61, 102337. [Google Scholar] [CrossRef]
  3. Suman. Air quality indices: A review of methods to interpret air quality status. Mater. Today Proc. 2021, 34, 863–868. [Google Scholar] [CrossRef]
  4. Xu, Z.; Kang, Y.; Cao, Y.; Li, Z. Deep amended COPERT model for regional vehicle emission prediction. Sci. China Inf. Sci. 2021, 64, 1–3. [Google Scholar] [CrossRef]
  5. Reames, T.G.; Bravo, M.A. People, place and pollution: Investigating relationships between air quality perceptions, health concerns, exposure, and individual-and area-level characteristics. Environ. Int. 2019, 122, 244–255. [Google Scholar] [CrossRef]
  6. Yang, X. The relationship between air pollution and human health. Green Build. Mater. 2019, 154, 53–56. [Google Scholar]
  7. Bai, L.; Shin, S.; Burnett, R.T.; Burnett, R.T.; Kwong, J.C.; Hystad, P.; Donkelaar, A.V.; Goldberg, M.S.; Lacigne, E.; Copes, R.; et al. Exposure to ambient air pollution and the incidence of congestive heart failure and acute myocardial infarction: A population-based study of 5.1 million Canadian adults living in Ontario. Environ. Int. 2019, 132, 105004. [Google Scholar] [CrossRef] [PubMed]
  8. Zheng, Y.; Liu, F.; Hsieh, H.P. U-air: When urban air quality inference meets big data. In Proceedings of the 19th ACM SIGKDD International Conference on Knowledge Discovery and Data Mining, Chicago, IL, USA, 11–14 August 2013; pp. 1436–1444. [Google Scholar]
  9. Qin, Y.; Qian, Y.; Rong, T. Research on the optimization of monitoring site selection based on atmospheric pollutants. China Environ. Sci. 2015, 35, 1056–1064. [Google Scholar]
  10. Liu, P.; Zheng, J.; Li, Z.; Zhong, L.; Wang, X. Research on Optimal Distribution Method of Regional Air Quality Monitoring Network. China Environ. Sci. 2010, 30, 907–913. [Google Scholar]
  11. Luo, H.; Feng, H.; Yan, X. Application of air quality prediction model for monitoring site selection. In Proceedings of the 2012 Annual Conference of Chinese Society for Environmental Sciences, Sapporo, Japan, 1–4 August 2012. [Google Scholar]
  12. USEPA. Fresno Supersite Final Report; Desert Research Institute for the OAQPS, U.S. Environmental Protection Agency: Research Triangle Park, NC, USA, 2005.
  13. Europe Environment Agency. The European Environment. State and Outlook 2005 [EB/OL]. Available online: http://www.eea.europa.eu/highlights/20051122115248 (accessed on 29 November 2005).
  14. Fukushima, H. Air pollution monitoring in East Asia. Sci. Technol. Trend 2006, 18, 55–63. [Google Scholar]
  15. Raffuse, S.M.; Sullivan, D.C.; McCarthy, M.C. Ambient Air Monitoring Network Assessment Guidance; U.S. Environmental Protection Agency: Washington, DC, USA, 2006.
  16. Arbeloa, S.F.J.; Caseiras, P.C.; Andres, L.P.M. Air quality monitoring: Optimization of a network around a hypothetical potash plant in open countryside. Atmos. Environ. 1993, 27A, 729–738. [Google Scholar] [CrossRef]
  17. Kao, J.J.; Hsieh, M.R. Utilizing multiobjective analysis to determine an air quality monitoring network in an industrial district. Atmos. Environ. 2006, 40, 1092–1103. [Google Scholar] [CrossRef]
  18. Sarigiannis, D.A.; Saisana, M. Multi-objective optimization of air quality monitoring. Environ. Monit. Assess. 2008, 136, 87–99. [Google Scholar] [CrossRef] [PubMed]
  19. Xu, Z.; Kang, Y.; Cao, Y. Emission stations location selection based on conditional measurement GAN data. Neurocomputing 2020, 388, 170–180. [Google Scholar] [CrossRef]
  20. Liu, D.; Weng, D.; Li, Y.; Bao, J.; Zheng, Y.; Qu, H.; Wu, Y. Smartadp: Visual analytics of large-scale taxi trajectories for selecting billboard locations. IEEE Trans. Vis. Comput. Graph. 2016, 23, 1–10. [Google Scholar] [CrossRef]
  21. Li, Y.; Zheng, Y.; Ji, S.; Wang, W.; U, L.H.; Gong, Z. Location selection for ambulance stations: A data-driven approach. In Proceedings of the 23rd SIGSPATIAL International Conference on Advances in Geographic Information Systems, Seattle, WA, USA, 3–6 November 2015; pp. 1–4. [Google Scholar]
  22. Kang, Y.; Li, Z.; Zhao, Y.; Qin, J.; Song, W. A novel location strategy for minimizing monitors in vehicle emission remote sensing system. IEEE Trans. Syst. Man Cybern. Syst. 2017, 48, 500–510. [Google Scholar] [CrossRef]
  23. Hsieh, H.P.; Lin, S.D.; Zheng, Y. Inferring air quality for station location recommendation based on urban big data. In Proceedings of the 21th ACM SIGKDD International Conference on Knowledge Discovery and Data Mining, Sydney, NSW, Australia, 10–13 August 2015; pp. 437–446. [Google Scholar]
  24. Zhu, X.; Ghahramani, Z.; Lafferty, J.D. Semi-supervised learning using gaussian fields and harmonic functions. In Proceedings of the 20th International Conference on Machine Learning (ICML-03), Washington, DC, USA, 21–24 August 2003; pp. 912–919. [Google Scholar]
  25. Vardoulakis, S.; Fisher, B.E.A.; Pericleous, K.; Gonzalez-Flesca, N. Modelling air quality in street canyons: A review. Atmos. Environ. 2003, 37, 155–182. [Google Scholar] [CrossRef] [Green Version]
  26. Bruna, J.; Zaremba, W.; Szlam, A.; LeCun, Y. Spectral networks and locally connected networks on graphs. arXiv 2013, arXiv:1312.6203. [Google Scholar]
  27. Defferrard, M.; Bresson, X.; Vandergheynst, P. Convolutional neural networks on graphs with fast localized spectral filtering. In Proceedings of the Advances in Neural Information Processing Systems, Bracelona, Spain, 9–10 December 2016; pp. 3844–3852. [Google Scholar]
  28. Hammond, D.K.; Vandergheynst, P.; Gribonval, R. Wavelets on graphs via spectral graph theory. Appl. Comput. Harmon. Anal. 2011, 30, 129–150. [Google Scholar] [CrossRef] [Green Version]
  29. Kipf, T.N.; Welling, M. Semi-supervised classification with graph convolutional networks. arXiv 2016, arXiv:1609.02907. [Google Scholar]
  30. Yu, B.; Yin, H.; Zhu, Z. Spatio-temporal graph convolutional networks: A deep learning framework for traffic forecasting. arXiv 2017, arXiv:1709.04875. [Google Scholar]
  31. Xu, Z.; Kang, Y.; Cao, Y.; Li, Z. Spatiotemporal Graph Convolution Multifusion Network for Urban Vehicle Emission Prediction. IEEE Trans. Neural Networks Learn. Syst. 2021. [Google Scholar] [CrossRef]
  32. Sun, J.; Zhang, J.; Li, Q.; Yi, X.; Liang, Y.; Zheng, Y. Predicting citywide crowd flows in irregular regions using multi-view graph convolutional networks. arXiv 2019, arXiv:1903.07789. [Google Scholar] [CrossRef]
  33. Cui, Z.; Henrickson, K.; Ke, R.; Wang, Y. Traffic graph convolutional recurrent neural network: A deep learning framework for network-scale traffic learning and forecasting. IEEE Trans. Intell. Transp. Syst. 2019, 21, 4883–4894. [Google Scholar] [CrossRef] [Green Version]
  34. Abu-El-Haija, S.; Perozzi, B.; Al-Rfou, R.; Alemi, A. Watch your step: Learning graph embeddings through attention. arXiv 2017, arXiv:1710.09599. [Google Scholar]
  35. Grover, A.; Leskovec, J. node2vec: Scalable feature learning for networks. In Proceedings of the 22nd ACM SIGKDD International Conference on Knowledge Discovery and Data Mining, San Francisco, CA, USA, 13–17 August 2016; pp. 855–864. [Google Scholar]
  36. Perozzi, B.; Al-Rfou, R.; Skiena, S. DeepWalk: Online Learning of Social Representations; ACM: New York, NY, USA, 2014; pp. 701–710. [Google Scholar]
  37. Shuman, D.I.; Narang, S.K.; Frossard, P.; Ortega, A.; Vandergheynst, P. The Emerging Field of Signal Processing on Graphs: Extending High-Dimensional Data Analysis to Networks and Other Irregular Domains. IEEE Signal Process. Mag. 2013, 30, 83–98. [Google Scholar] [CrossRef] [Green Version]
  38. Smola, A.J.; Schölkopf, B. A tutorial on support vector regression. Stat. Comput. 2004, 14, 199–222. [Google Scholar] [CrossRef] [Green Version]
  39. Gibson, M.D.; Kundu, S.; Satish, M. Dispersion model evaluation of PM2.5, NOx and SO2 from point and major line sources in Nova Scotia, Canada using AERMOD Gaussian plume air dispersion model. Atmos. Pollut. Res. 2013, 4, 157–167. [Google Scholar] [CrossRef] [Green Version]
  40. Rakowska, A.; Wong, K.C.; Townsend, T.; Chan, K.L.; Westerdahl, D.; Ng, S.; Močnik, G.; Drinovec, L.; Ning, Z. Impact of traffic volume and composition on the air quality and pedestrian exposure in urban street canyon. Atmos. Environ. 2014, 98, 260–270. [Google Scholar] [CrossRef]
  41. Blum, A.; Mitchell, T. Combining labeled and unlabeled data with co-training. In Proceedings of the 11th Annual Conference on Computational Learning Theory, Madison, WI, USA, 24–26 July 1998; pp. 92–100. [Google Scholar]
  42. Chen, L.; Cai, Y.; Ding, Y.; Lv, M.; Yuan, C.; Chen, G. Spatially fine-grained urban air quality estimation using ensemble semi-supervised learning and pruning. In Proceedings of the Acm International Joint Conference on Pervasive & Ubiquitous Computing, Heidelberg, Germany, 12–16 September 2016; pp. 1076–1087. [Google Scholar]
Figure 1. Beijing air-quality distribution on a certain day. The monitoring stations in the red circle are very geographically close, but the PM 10 concentration varies greatly throughout the year.
Figure 1. Beijing air-quality distribution on a certain day. The monitoring stations in the red circle are very geographically close, but the PM 10 concentration varies greatly throughout the year.
Remotesensing 13 01600 g001
Figure 2. The location of the Beijing air-quality experimental sites.
Figure 2. The location of the Beijing air-quality experimental sites.
Remotesensing 13 01600 g002
Figure 3. Division of urban areas based on geographic information.
Figure 3. Division of urban areas based on geographic information.
Remotesensing 13 01600 g003
Figure 4. The general framework. (1) Air-quality-inference model: HGCNInf consists of a higher -order graph convolution network (HGCN) and several fully connected neural networks (FNN). By inputing historical monitoring data and by fetching global information, the HGCN is used to capture the spatiotemporal interactions of air-quality distribution and the FNN is used to extract the features of external factors. (2) Monitoring-station location recommendation model: by using the convolution weight parameters of HGCNInf, we designed a GMIE based on the correlation degree between nodes in USTG reflecting the air-quality spatiotemporal changes, marking the recommendation priorities of unlabeled nodes according to the ability to improve the inference accuracy of HGCNInf via the node incremental learning method.
Figure 4. The general framework. (1) Air-quality-inference model: HGCNInf consists of a higher -order graph convolution network (HGCN) and several fully connected neural networks (FNN). By inputing historical monitoring data and by fetching global information, the HGCN is used to capture the spatiotemporal interactions of air-quality distribution and the FNN is used to extract the features of external factors. (2) Monitoring-station location recommendation model: by using the convolution weight parameters of HGCNInf, we designed a GMIE based on the correlation degree between nodes in USTG reflecting the air-quality spatiotemporal changes, marking the recommendation priorities of unlabeled nodes according to the ability to improve the inference accuracy of HGCNInf via the node incremental learning method.
Remotesensing 13 01600 g004
Figure 5. The architecture of proposed higher-order convolutional network. The multi-hop neighbor with respect to a node (the red star) is the nodes that can be reached by the black edges.
Figure 5. The architecture of proposed higher-order convolutional network. The multi-hop neighbor with respect to a node (the red star) is the nodes that can be reached by the black edges.
Remotesensing 13 01600 g005
Figure 6. The framework of the air-quality-inference model. The specific structure of the higher-order grapth convolutional network is shown in Figure 4.
Figure 6. The framework of the air-quality-inference model. The specific structure of the higher-order grapth convolutional network is shown in Figure 4.
Remotesensing 13 01600 g006
Figure 7. The basic framework of the monitoring-station location recommendation model.
Figure 7. The basic framework of the monitoring-station location recommendation model.
Remotesensing 13 01600 g007
Figure 8. The performance comparison of each inference model.
Figure 8. The performance comparison of each inference model.
Remotesensing 13 01600 g008
Figure 9. The robustness of each air-quality-inference model when gradually reducing the number of labeled nodes.
Figure 9. The robustness of each air-quality-inference model when gradually reducing the number of labeled nodes.
Remotesensing 13 01600 g009
Figure 10. The HGCNInf training efficiency of different orders.
Figure 10. The HGCNInf training efficiency of different orders.
Remotesensing 13 01600 g010
Figure 11. The HGCNInf weight visualization.
Figure 11. The HGCNInf weight visualization.
Remotesensing 13 01600 g011
Figure 12. The recommendation results of each model. The nodes marked with orange symbols in the figure are candidate nodes. Each candidate node has a recommendation priority marked by Arabic numerals (the smaller the number, the higher the recommendation priority).
Figure 12. The recommendation results of each model. The nodes marked with orange symbols in the figure are candidate nodes. Each candidate node has a recommendation priority marked by Arabic numerals (the smaller the number, the higher the recommendation priority).
Remotesensing 13 01600 g012
Figure 13. The MAPE results when gradually increase the number of recommended stations according to its priority.
Figure 13. The MAPE results when gradually increase the number of recommended stations according to its priority.
Remotesensing 13 01600 g013
Figure 14. The RMSE results when gradually increase the number of recommended stations according to its priority.
Figure 14. The RMSE results when gradually increase the number of recommended stations according to its priority.
Remotesensing 13 01600 g014
Table 1. The types of POIs.
Table 1. The types of POIs.
1.Parks6.Companies
2.Schools7.Hotels
3.Stadium8.Supermarkets
4.Restaurants9.Vehicle repair stations
5.Business center10.Gas stations
Table 2. The performance comparison of inference models.
Table 2. The performance comparison of inference models.
ModelBeijing Air-Quality Data Evaluation
MAPE (%)RMSE
Without External FactorsWith External FactorsWithout External FactorsWith External Factors
HA56.35/30.01/
SVR49.94/23.00/
AQInf35.2521.3116.9210.35
GCNInf28.5214.1013.937.24
HGCNInf19.187.719.993.90
Table 3. The robustness comparison of inference models.
Table 3. The robustness comparison of inference models.
The Robustness of Each Air Quality Inference Model
ModelHASVRAQInfGCNInfHGCNInf
MetricMAPE (%)RMSEMAPE (%)RMSEMAPE (%)RMSEMAPE (%)RMSEMAPE (%)RMSE
Number
056.3530.0149.9423.0021.3110.3514.107.247.713.90
159.4531.5051.7524.4023.7511.4115.827.749.013.98
263.3834.0252.1526.424.1512.0116.137.8211.636.62
369.6737.9955.1829.0027.1813.3918.318.6112.076.66
495.0542.9662.0033.0040.6819.5220.2710.6915.187.52
5131.9250.1873.0140.1861.0627.1826.4212.3422.9411.67
Table 4. The HGCNInf performance of different orders.
Table 4. The HGCNInf performance of different orders.
OrderK = 1K = 2K = 3K = 4K = 5
Metric
MAPE (%)Without external factors23.7523.0919.1820.2720.27
With external factors15.1815.187.7111.6311.63
RMSEWithout external factors11.4111.759.9910.6910.69
With external factors7.527.523.906.626.62
Table 5. The HGCNInf training time of different orders.
Table 5. The HGCNInf training time of different orders.
Training Time/(s)K = 1K = 2K = 3K = 4K = 5
Without external factors1.792.623.504.615.71
With external factors5.839.4915.2617.3920.75
Table 6. The improved MAPE and RMSE results with varying numbers of recommended locations.
Table 6. The improved MAPE and RMSE results with varying numbers of recommended locations.
The Improved MAPE and RMSE Results with the Increase in the Number of Recommended Nodes
ModelHGCNInf+GMIEGCNInf+GMIEAQInf+GEM
MetricMAPE (%)RMSEMAPE (%)RMSEMAPE (%)RMSE
Number
021.5710.8323.3111.2624.1512.01
112.946.8420.4810.1323.7511.41
29.014.2718.308.6122.8011.22
38.003.9816.377.8222.1010.69
47.713.9014.707.4021.6010.47
57.713.9014.107.4221.3110.35
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kang, Y.; Chen, J.; Cao, Y.; Xu, Z. A Higher-Order Graph Convolutional Network for Location Recommendation of an Air-Quality-Monitoring Station. Remote Sens. 2021, 13, 1600. https://0-doi-org.brum.beds.ac.uk/10.3390/rs13081600

AMA Style

Kang Y, Chen J, Cao Y, Xu Z. A Higher-Order Graph Convolutional Network for Location Recommendation of an Air-Quality-Monitoring Station. Remote Sensing. 2021; 13(8):1600. https://0-doi-org.brum.beds.ac.uk/10.3390/rs13081600

Chicago/Turabian Style

Kang, Yu, Jie Chen, Yang Cao, and Zhenyi Xu. 2021. "A Higher-Order Graph Convolutional Network for Location Recommendation of an Air-Quality-Monitoring Station" Remote Sensing 13, no. 8: 1600. https://0-doi-org.brum.beds.ac.uk/10.3390/rs13081600

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop