• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Parking Availability Prediction with Coarse-Grained Human Mobility Data

    2022-08-23 02:15:42AuroraGonzalezVidalFernandoTerrosoenzandAntonioSkarmeta
    Computers Materials&Continua 2022年6期

    Aurora Gonzalez-Vidal,Fernando Terroso-Sáenzand Antonio Skarmeta

    1University of Murcia,Murcia,30003,Spain

    2Universidad Católica de Murcia(UCAM),Murcia,30107,Spain

    Abstract: Nowadays,the anticipation of parking-space demand is an instrumental service in order to reduce traffic congestion levels in urban spaces.The purpose of our work is to study,design and develop a parking-availability predictor that extracts the knowledge from human mobility data,based on the anonymized human displacements of an urban area, and also from weather conditions.Most of the existing solutions for this prediction take as contextual data the current road-traffic state defined at very high temporal or spatial resolution.However,access to this type of fine-grained location data is usually quite limited due to several economic or privacy-related restrictions.To overcome this limitation, our proposal uses urban areas that are defined at very low spatial and temporal resolution.We conducted several experiments using three Artificial Neural Networks: Multilayer Perceptron, Gated Recurrent Units and bidirectional Long Short Term Memory networks and we tested their suitability using different combinations of inputs.Several metrics are provided for the sake of comparison within our study and between other studies.The solution has been evaluated in a real-world testbed in the city of Murcia(Spain)integrating an open human-mobility dataset showing high accuracy.A MAPE between 4%and 10%was reported in horizons of 1 to 3 h.

    Keywords:Parking space;human mobility mining;recurrent neural networks;prediction

    1 Introduction

    At the dawn of the smart cities era,municipalities are still facing a lack of proper transportation policies for the efficient management of human mobility.As a result,problems related to road traffic congestion have become challenging issues within urban settlements.In economic terms,it is estimated that traffic jams cost$87 millions in the USA in 20191https://inrix.com/scorecard/.Also,traffic congestions have a direct impact on the air quality of a region since petrol vehicles are one of the most important sources of air pollutants[1].

    The development of services for parking space management is a sensible solution to reduce trafficcongestion,since roughly 40%of road traffic in cities is caused by drivers searching for parking space2https://www-03.ibm.com/press/us/en/pressrelease/35515.wss.In that sense,the anticipation of future parking-space demand would allow operators to develop shortterm and long-term counteractions to avoid mobility problems in urban regions.

    For that prediction,a paramount contextual factor to be considered is the road-traffic conditions of the urban areas.In that sense,a large course of action for parking availability prediction leverages the endless deployment of new and more powerful sensors across cities thanks to the Internet of Things(IoT)[2].More in detail,mechanisms based on static infrastructure sensors[3],mobile ones[4]or the crowdsensing paradigm[5]have already been proposed.In general terms,these solutions rely on highresolution mobility data that captures human or vehicle displacements at quite fine-grained temporal or spatial scales.For example,a solution based on mobile sensors mounted on floating vehicles(e.g.,taxis or buses)would allow monitoring the state of each road segment.

    However, high-resolution human mobility data is quite sensitive in terms of privacy.Its access is rather limited because of regulatory and economic policies [6].At the same time, the open data movement has promoted the release of an increasing number of human-mobility datasets [7–9].However, they are usually filtered and pre-processed for anonymization.Hence, open mobility data provides a much lower temporal and spatial resolution than ad-hoc mobility feeds.The development of intelligent systems leveraging such coarse-grained mobility data is still scarce in the domain.

    In this context,the present work introduces a novel mechanism to forecast the number of available spaces of a set of urban parking lots in the city of Murcia(Spain).Our system relies on two datasources.

    On the one hand, the predictor is timely fed with the sheer number of free spaces in the target parking lots.This is provided by a set of counter sensors installed in the parking lots which are,in turn,part of a wide IoT infrastructure deployed in the city.On the other hand,the system also considers the anonymized human trips co-occurring within the city.Such trips are extracted from a nationwide open mobility dataset released by the Spanish Ministry of Transportation(SMT)considering mobile-phone location data from several telephone operating carriers3https://www.mitma.es/ministerio/covid-19/evolucion-movilidad-big-data/opendata-movilidad.To ensure users’privacy,trips are defined on a spatial tessellation of the country based on large geographical areas.

    These two sources are used to train an Artificial Neural Network (ANN) able to anticipate the values from the IoT counters several hours ahead.Given the sequential nature of the inputs,two types of Recurrent Neural Networks(RNNs)have been evaluated,a LSTM[10]and a GRU[11]one.Both types of architectures can learn long and short-term patterns from a given input sequence but some slight differences in terms of complexity exist between them.Nevertheless,both have been successfully applied to the parking availability prediction problem [12].Hence, the present work also relies on the hypothesis that the latent human movement of a region provides valuable information about its parking occupancy level so such knowledge should be integrated as part of the prediction system.

    It is important to note that the inputs are defined at different spatial scales and sources.The occupancy of the parking lots is extracted from an IoT infrastructure,and,in spatial terms,such data is related to very particular spots within the city.However, the human trips dataset only provides a region-based view of the actual human mobility.Therefore, our solution studies whether human mobility data defined at a very high level gives insight into road-traffic behaviour related to city locations captured by the parking lots’sensors.Some parking-space predictors incorporate in their pipelines data defined at a coarse spatial granularity such as weather conditions.However,to the best of the authors’knowledge,there is a scarcity of proposals making use of coarse-grained feeds related to human movement information.This targets an important research gap in the parking availability prediction domain.

    Furthermore,our work also brings three important operational benefits with respect to existing solutions.Firstly,we can deploy it in cities without an IoT infrastructure able to capture fine-grained mobility data.Secondly,as a side effect,it can be regarded as a cost-effective mechanism.Finally,our approach relies on pre-anonymized human data so it will limit the privacy concerns among end-users.

    The rest of the paper is structured as follows.Section 2 provides an overview of current approaches for parking space prediction from different perspectives.Next,Section 3 puts forward the proposed forecasting mechanism and the main results of the performed experiments are described in Section 4.Finally,the main conclusions and the future work are summed up in Section 5.

    2 Related Work

    This section provides an overview of existing approaches for parking space prediction.We gathered three dimensions of the solutions, their input sources, applied algorithms and the urban service that they provide.Last, we also briefly review the main findings within the human mobility mining discipline.

    2.1 Input Sources

    Sensors,parking meters and crowdsensing are common sources for collecting real-time parking availability.

    Regarding static sensors, the mechanisms in [4,13,14] use parking meter data for its prediction.Counting data from parking lots at each moment is used in [15–21] where the available slots are collected every 5–15 min.Infrared and ultrasonic sensors can be used to identify a car in a particular spot[2], as well as Radio-Frequency Identification(RFID) [22] for single vehicles that are provided with a radio frequency tag.Furthermore,ferromagnetic sensors installed under the asphalt are used in[23].Others works focus on the ticket data collected from a smart-parking system[24].In terms of communication protocols,some proposals use an energy efficient LoRaWAN channel to transmit the collected data to a centralized repository[25].

    Concerning mobile sensors,most works rely on vehicle-mounted devices sensing traffic contextual features[13].A similar approach is followed in[4]where the GPS traces of the fleet of taxis are used for parking availability monitoring.More generalist ways to collect data use machine vision[26]and car count in order to know the number of free spaces in real-time.In this scenario, crowdsensing focuses mainly on developing smartphone apps to ensure the accuracy of parking information and location[5].

    Lastly, weather conditions might affect parking demand at different degrees.For that reason,several works have incorporated this data as part of their input [13,15,20].Other contextual urban feed that has also used in some works has been Point of Interest(POI)data indicating the density of certain venues like restaurants or bars or schools[20,21],since the POI distribution of a region strongly affects its human mobility patterns[27].

    The present work takes as input count-based parking data and weather conditions as other proposals in the literature.Our approach also considersmobiledata as input in the form of human mobility flows but its spatial and temporal granularity is much larger than the sources of previous works [4,13].As a matter of fact, the temporal resolution of the traffic dataset in [4] was fixed to 10 min instead of 1 h as in our case as it is pointed out in Sec.3.2.1.Furthermore,whilst works based on mobile data rely on street-based reports [4,13] our solution considers the incoming or outgoing human flows of large-scale geographical urban areas.For estimating human flows,other works rely on social media publications(geo-tagged information)[28,29].Unlike previous proposals,our approach combines datasources of different nature.Those are raw sensor data from the counters installed in the parking lots and transmitted through a IoT infrastructure(see Sec.3.1).and open data capturing the human flows in the target city at very low resolution.As we have seen,this fusion of raw sensor data and an open mobility feed for parking availability prediction has not been fully explored in the literature.

    2.2 Applied Methods

    Different ANN models have been applied for parking demand forecasting.Since this can be regarded as a time-series prediction problem,most works have used different types of Recurrent Neural Networks(RNNs),such as Long short-term memory(LSTM)[14,17,18]models.Moreover,the work in [13] combines a Graph Neural Network (GNN) and a LSTM model so as to also consider the underlying road topology of the target area.Hierarchical GNNs have been also used to predict the availability of different parking lots since they capture well the spatial correlations among different entities[21].Apart from that,Multi-layer Perceptrons(MLPs)have been also proposed and used as a baseline to compare with other algorithms,especially deep learning algorithms[23,30,31].GRU and LSTM have been proven to be the best RNN architectures for solving a wide set of sequential data problems[32].In[12],the GRU architecture achieves better results in nearly all scenarios(4 countries and different exogenous variables) compared to the LSTM version for 6 h-ahead prediction.GRU have also been successfully applied to capture dynamic temporal autocorrelations between several parking lots [33].In [25], authors proposed a network architecture combining three types of layers,1D convolutional, a LSTM and three fully connected ones for the prediction task.Also, generative adversarial networks (GANs) that merge spatio-temporal and weather data from different urban datasources have been used[20].

    Furthermore,researchers in[15]make use of Markov chain analysis as they assume the Markov property in the context of parking space prediction.Similarly,a continuous-time Markov queue is put forward in[16].A likelihood maximization approach is applied to deal with the heterogeneity of crowdsensed data is presented in[5].Other methods like Kalman filters have been also studied[4].In addition to that,some proposals have used ensemble-learning algorithms like Random Undersampling Boosts to perform the prediction task [24].In a similar way, [19] proposes an ensemble model where three different algorithms,Random Forest regressor,Gradient boosting regressor and Adaboost regressor are combined by means of bagging and boosting techniques to provide a final prediction outcome.

    In terms of applied models,our proposal has evaluated three different models,MLP,GRU and bi-LSTM ones in order to find the most suitable hand for the task in hand as it is described in Sec.4.4.2.

    2.3 Provided Urban Services

    Different parking scenarios:research works deal with improving efficiency at closed parking lots which are paid parking lots(both indoors and outdoors),andpen parking lots or ticketing.

    As long as the first scenario is concerned,a mechanism to predict the occupancy of three private and public parking lots in the next 60 min in Geneva (Switzerland) is proposed in [15].A set of 14 parking lots in San Francisco (USA) are used as a use case in [16].In this case, the future available parking places area provided a time horizon of up to 12 min.In the same city, the proposal in[4] controls on-street parking spaces of 579 road segments.The work in [17] provides a prediction mechanism for a single parking lot in Shenzhen (China) with a prediction horizon of 30 min.In[21]the proposed system was able to provide predictions for a set of different parking lots at Beijing and Shenzhen 15, 30 and 45 min ahead.A holistic IoT infrastructure that able to make predictions of 30 parking locations of Birmingham (UK) is proposed in [18].Unlike previous approaches, the predictions are given in hour intervals instead of minutes and the same dataset is used for 30-minute time horizons in[19].

    Regarding the second group of on-street parking spaces,the work in[13]applies its predictor to a downtown area in Pittsburgh(USA).In[20],the prediction targets an outdoor regulated parking zone in Santander (Spain).In both cases, the prediction horizon is up to 30 min.In [5], the crowd-based prediction mechanism is tested in seven on-street public parking spaces in London(UK).In[14],3000 on-street parking slots were controlled by the solution to provide predictions up to 30 min ahead.The problema can be shifted from regression to classification[24].This way,the solution is able to classify the level of service at different on-street parking areas in Madrid(Spain).In[25],the predictor focuses on a street in the city of Montova(Italy)providing parking lots prediction from 1 to 8 h ahead.

    Our work forecasts the available parking places of two public garages with a larger time granularity(hourvs.minutes)and a wider time horizon of the prediction,as is stated in Sec.4.4.2.

    2.4 Human Mobility Mining

    Human mobility repeats over time(e.g.,people tend to go to work at the same time every working day)so it is usually modeled as a time pattern highly predictable at some extend[34].The prediction of where and when people are going to move is an instrumental tool in many domains[35].

    Nowadays it is possible to capture location data at multiple spatial and temporal scales from a large palette of technologies like GPS, Bluetooth or WiFi.Based on these technologies, a variety of machine learning models have been proposed for human mobility forecasting [36].These solutions basically rely on the raw spatio-temporal trajectories from movingentities, like taxis or individuals,equipped with the location technologies.This allows collecting their current location every few minutes or even seconds.

    In this context,human mobility patterns have already been successfully extracted from such open feeds and used in different predictive problems[37].However,there is a scarcity of proposals leveraging such mobility coarse-grained feeds as reliable input sources to improve other predictive tasks.

    3 Setting Overview

    This section describes the urban scenario for testing our approach.In that sense,we describe the available data sources and their suitability as input features for parking availability prediction.

    3.1 General Description

    The feasibility of our solution has been evaluated in Murcia, a city in the South-East of Spain.Fig.1 depicts the target urban area within this city.The city of Murcia is involved in a holistic smartcity project called MiMurcia4https://mapamurcia.inf.um.es/that is developing a platform that interconnects data and services to provide citizens with personalised and updated information about the status of the city and to make it smart.Murcia is a city with a great trajectory in the use of IoT data for providing smart services such as water [38] and energy management [39].Within this platform, several data is published:irrigation, problems reported by citizens (with pictures) and the parking and weather information that our solution uses.

    Figure 1:Map of the target city.Each colored polygon is a particular Mobility area from the mobility study whereas the dots indicate the locations of the two target parking lots

    Our solution would eventually enrich this ecosystem with a decision-support system for road traffic management.This way,it would provide drivers with real-time information about the available parking spaces in their destination areas on the basis of their time of arrival to such destinations.

    3.2 Target Datasets

    Given the urban infrastructure described in the previous section,we focused on the following data sources related to the city to perform the prediction task.

    3.2.1 Low-Resolution Human Mobility

    As it was mentioned in Sec.1, we used the nationwide human mobility dataset released by the SMT in December 2020.This dataset covers a 9-month period from February 29th to November 30th,2020,and it indicates the number of trips among 3216 ad-hoc administrative areas(hereby Mobility Areas,MA)per hour in Spain both in its peninsular and insular extension.In that sense,a single trip stands for the spatial displacement of an individual with a distance above 500 meters.Consequently,this dataset can be regarded as a set of tuples taking the form〈date,hour,morigin,mdest,ntrp〉,reporting that there was ntrphuman trips from the MA moriginto mdestduring the indicateddateandhour.

    This data has been collected through Call Detail Records (CDRs) from 13 million users of an unspecified mobile-phone carrier [40].Once anonymised, representative mobility statistics at the nation-level of the population of Spain were inferred and made publicly available open data.In its raw form, the dataset comprises 830, 450, 300 trips among MAs.To filter the trips related to the target city,we focused on the 3 particular MAs that spatially cover its geographical extension.Fig.1 shows this set of three MAs,M=〈mA,mB,mC〉with areas of 5.08,3.04 and 3.66 km2 and populations of 54,279 people,47,249 and 46,738respectively.Next,we extracted the time series for each target area,comprising the time-ordered sequence of incoming trips ih towards m ∈M per hour h.This way,the mobility profile of each MA of the city was composed in terms of its incoming flows of people.Given the geographical and population magnitudes of these MAs,their profiles provide much more coarse-grained mobility information than one based on road segments as much of the actual state of the art relies on.

    3.2.2 Parking Space

    We are collecting also the occupancy information from two underground parking lots in the area that operate 24/7,La Vega and Libertad parking.The actual locations of both of them are depicted in Fig.1 as colored dots.In particular,the number of free spaces is collected every two minutes for each of them.This will constitute the target variable of our system.

    The two parking lots are very close toThe Circular Squareof Murcia,which is one of the most important squares in the city.Very crowded streets link to this place, including theGran Vía, full of commerce and bustle.La Vega Hospitalis also in the area.The data was aggregated into hourly intervals in order to make it coincide with that of mobility flows and the tables were joined.Adding the number of parking slots from both parking we count with 642 spaces.Last,the quality of the data from this source has already been tested in previous studies[11].

    3.2.3 Weather Conditions

    Wind (strength and direction), rainfall, daylight, relative humidity and temperature have been proven to have an impact on the choice of travel modes(walk,public transportation and vehicles)[41]and therefore,we have also studied the suitability of weather data as input for the target predictor.Since people would normally make a decision based on current weather,we decided not to incorporate future weather information into the system.We used the IMIDA (The Research Institute of Agriculture and Food Development of Murcia) service in order to obtain the historical weather information5https://www.imida.es/.They provided us with an hourly historical set of data including the following variables:temperature(mean, min and max) (?C), humidity (mean, min and max) (%), radiation (mean and max) (w/m2),wind speed (mean and max) (m/s2), wind direction (mean)(degrees), precipitations (mm), dew point(?C)and vapour pressure deficit(kPa)from the closest weather station to the studied zone.

    3.3 Correlation Study

    Given these input sources, we carried out a correlation study to assess whether the extracted mobility flows and the weather parameters are reliable sources for parking-space prediction.In that sense, statistical correlation measures the relationship between two random variables, whether it is causal or not.Although it could refer to any statistical association, is commonly related to linear correlation.

    In our scenario, we can assume that our data is normal, which means that the distribution of each variable fits a Gaussian distribution.As the sample becomes larger, the sample’s distribution approximates to a normal distribution according to the Central Limit Theroem[10.Therefore,given our scenarios,we can use Pearson’s correlation coefficient to summarize the correlation between the variables.

    Pearson’s correlation coefficient is the covariance of the two variables divided by the product of their standard deviations.Given paired data{(x1,y1).,...,(xn,yn)}consisting of n pairs,it is written as follows,

    Pearson’s correlation coefficient is a number between-1 and 1 that describes a negative or positive correlation respectively.A value of zero indicates no correlation.

    Fig.2 shows the correlation heatmap between the free parking spaces (free), the considered mobility flows〈ImA,ImB,ImC〉namedflow1,flow2andflow3in the figure and the weather data.From this,we chose the first two flows since they show a greater negative correlation with the output variable than the third.We also chose the mean temperature(tmed),mean radiation(radmed)and mean wind speed(vvmed).The rest of the weather variables are highly correlated amongst themselves and would add noise to the model.Also, even though from this correlation analysis we find that precipitations(prec)might not be crucial,this could be due to the lack of rain in the studied region.In that sense,we decided to include it nevertheless because it is an important factor when choosing transportation.

    Figure 2:Correlation plot of all considered input variables

    Regarding the mobility flows, we can see from Fig.2 that the three of them have a negative correlation with the parking-space variable.However,flows ImAand ImBhad a slightly higher negative correlation(-0.642,-0.719)than the flow for the area mC(-0.437)so we opted for just including the two former flows to the predictor.This difference in the coefficients makes sense as area mCis the one more distant to the target parking lots location(see Fig.1).Hence,it associated latent mobility had a lower impact on the actual demand of the parking lots.

    For the sake of clarity,Fig.3 shows the two mobility flows selected from the analysis,ImAand ImB.Both flows are depicted along with the evolution of available parking spaces during the whole period of study.

    Figure 3: The two human mobility flows finally used for the prediction task along with number of available parking spaces during the period of study (a) Available parking spaces and the ImA flow during the period of study(b)Available parking spaces and the ImB flow during the period of study

    Furthermore,there are many phenomena in which the past influences the present in the sense that they can be used to foretell what will happen.When such phenomena are represented as a time series,they are said to have an auto-regressive property.The correlation for univariate time series observations with previous time steps,called lags,is an autocorrelation.A plot of the autocorrelation of a time series by lag is called the AutoCorrelation Function(ACF).In this case,instead of correlation between two different variables,the autocorrelationat lagkis between two values of the same variable at timesxiandxi+k.The previous formula is then transformed to the following,

    Second, the heroine is accomplishing a difficult task, removing a settled stain from clothing. Psyche, in Cupid and Psyche, has to perform three impossible tasks to prove her devotion to Cupid.

    Finally,partial autocorrelation is a summary of the relationship between an observation in a time series with observations at prior time steps with the relationships of intervening observations removed.The autocorrelation for an observation and some prior time step of itself is comprised of both the direct correlation and indirect correlations.These indirect correlations are a linear function of the correlation of the observation,with observations at intervening time steps.

    Fig.4 shows that parking spaces have a positive autocorrelation in the first 7 lags, reaches its maximum negative autocorrelation in lag 12 and then again the best positive correlation is present 24 h later.Regarding the partial autocorrelation, we see that there are many values statistically significant.From this we can conclude that the free parking presents an autorregresive behaviour and that including lagged values as inputs of the predictive models will likely imply a greater accuracy.

    Figure 4:The visualization of the autocorrelation and partial autocorrelation of free parking spaces(a)Autocorrelation parking availability(b)Partial autocorrelation parking availability

    4 Proposed Solution for Parking Space Prediction

    In the light of the available datasources in the urban scenario and their suitability as independent variables for the prediction task at hand, we are now able to properly define a solution for parking space prediction based on them.

    Our system takes as input three sources:historical data of available spaces in the target parking lots, the human mobility flows at low resolution and the weather conditions in the parkings’surrounding areas.

    4.1 Problem Formulation

    The parking space prediction problem that we focus on can be formulated as follows:

    Given the hour h ∈〈0,...,23〉,the number of available parking spaces during the last hprev hours from the two target parking lots,Ph=(ph,ph-1,..,ph-hprev),the number of incoming trips of a set of MAsand the weather conditions of the region Wh=〈wh,wh-1,..,wh-hprev〉,Find a mapping function F(Ph,Ih,Wh)→ph+T,where ph+Tis the sheer number of available parking spaces in the parking lots in the next T hours.

    4.2 Candidate Prediction Methods

    4.2.1 MLP Model

    This is one of the most popular and successful ANNs architectures that for a wide range of regression problems[42].Its inner architecture basically comprises an input layer,one or more hidden layers and an output one where each layer is connected to another one by non-linear functions.Each layer comprises a set of neurons connected to all neurons of the next layer but not to each other.In computational terms,the output of a MLP layer be defined as h(W,b)(X)=φ(XW+b),where W is the connections’weight matrix,X is the matrix with the input features,b is the vector with the bias terms andφis the activation function.

    4.2.2 GRU Model

    Given the sequence nature of the two validated input sources of our approach, we have used a GRU model for the prediction task [43].This is a foremost variant of Recurrent Neural Networks(RNN)able to learn short-term and long-term patterns in sequences of data.Unlike other RNNs,like LSTM,a GRU model has a slightly simpler structure which makes it faster to train[44].

    Fig.5a shows the general structure of a GRU network whereas Fig.5b depicts the inner structure of a GRU cell.As we can see,a GRU model just follows the composition of a regular RNN.Concerning the structure of the cells,they make use of a gated mechanism to memorize long-term patterns in the target sequence.Thus,a cell receives as input the current input vector x(t)and the previous state vector h(t- 1).Then,the cell generates the associated output y(t)which is also the state vector h(t)of the next cell.

    Figure 5: Illustrative architecture of the GRU model (a) A GRU model with three layers unrolled thought time(b)Gated structure of a GRU cell.FC stands for Fully Connected

    More in detail, the cell comprises three different gates, the update z(t), the reset r(t) and the memory-content g(t).The computations of each gate are as follows,

    where W{z,r,g}are the weight matrices for the input xt,U{z,r,g}are the weight matrices for the connections to the previous short-term state h(t-1)and b{z,r,g}are the bias terms of each layer.

    4.2.3 Bi-LSTM Model

    Last, we have evaluated a Long-Short Term Memory (LSTM) model.This is also a wellestablished RNN variant[45].We have evaluated a bidirectional LSTM model[46].Unlike atraditionalLSTM network, a bidirectional approach allows the model not only to learn “causal” patterns by means of its past and present inputs but also to“l(fā)ook ahead”into the future.As Fig.6a depicts,this is done by concatenating two RNNs one reading the input from left to right and the other from right to left.Then,the output of each cell is combined with its counterpart.

    Figure 6:Illustrative architecture of the bidirectional LSTM model.(a)A bidirectional LSTM layer with six cells(b)Gated structure of a LSTM cell.FC stands for Fully Connected

    The inner structure of a LSTM cell is depicted in Fig.5b,where h(t-1)indicates the short-term state at time instant and x(t)is the current input vector at instant t.Concerning the cell outputs,y(t)is the predicted availability at instant t whereas c(t)is the long-term state that traverses the network from left to right.Focusing on the four inner gates of the cell that modulates the outputs of the model,they can be formulated as follows:

    where x(t)is the parking availability at time instant t,W{f,i,o,g}are the weight matrices for this feature,U{f,i,o,g}are the weight matrices for the connections to the previous short-term state h(t-1)and b{f,i,o,g}are the bias terms of the four gates.

    As we can see,the key difference between the GRU and the LSTM cell is that the former has three layers that are reset r(t),update z(t)and memory-content g(t)gates while the latter has four layers that are the input i(t),output o(t)and forget f(t)gates along with an output layer g(t).Thus,GRU is less complex than LSTM because it has a smaller number of gates.

    4.3 Generation of the Models

    4.3.1 Data Pre-Processing

    The chosen period of analysis consisted of two ranges: 13th of July-31st of July and 29th of September-19th of October.We avoided August and September because in Murcia activities are reduced given the temperatures and holidays.Also,COVID-19 restrictions were more severe in those periods,so we decided to let them out of the analysis in order not to contaminate the sample.In the end,we considered 927 data points in roughly 2 months.

    The raw parking data and the raw weather data are collected with more granularity than the mobility flows,that is 4 and 15 min respectively.However,we homogenised all data in hourly intervals,which was the mobility flows’granularity.Given that we had to resample to a wider time frame,we used down-sampling and the aggregation function was the mean.After downsampling,due to failures in the communication system in the parking set,a 6%of the sample was missing.In order to impute the values,we used linear interpolation.

    The performance of Deep Learning models and regression models is improved when scaling[47].Scaling the inputs and outputs used in the training of the model is meaningful since small weights and errors in prediction values are used.Using unscaled inputs may slow or destabilize the learning process,while un-scaled target variables in regression problems may cause gradient explosions,leading to failure of the learning process.For that purpose, we have used Standardization (Z-score Normalization) in all our inputs.Feature standardization makes the values of each feature in the data have zero-mean(when subtracting the mean in the numerator)and unit variance.The formula is

    4.3.2 Models’Configuration

    All models had a training rate of 70 %, batch size 32 and learning factor 0.01, Mean Squared Error (MSE) as loss function and the Adam optimizer was used (α= 0.001,β1=0.9,β2=0.999,epsilon=10-7).The hyperparameters of the models were selected using a grid search with Tab.1 configuration.MLP used the ReLU activation function with 3 layers,50 cells per layer and 50 epochs.GRU used the Hyberbolic tangent function for activation with 2 layers, 50–20 cells in each of them and 100 epochs.Finally,LSTM used the same activation function with 3 layers,50–20–20 cells in each of them and 400 epochs.

    Table 1: Configuration of the grid search for the models’hyperparameters

    Fig.7 shows the learning curves of the models with the provided configuration.As we can see,we avoid overfitting all models with the provided datasets.

    Figure 7:Learning curves of the three target models in the experiment.(a)MLP model(b)GRU model(c)Bi-LSTM model

    Last,the number of parameters of these three models when they are fed with the incoming human flows,available parking spaces and weather conditions are 8801 in the case of the MLP,13611 for the GRU model and 33321 for the bi-LSTM model.

    4.4 Evaluation of the Approach

    4.4.1 Metrics

    4.4.2 Obtained Results

    Tab.2 shows the metric values obtained by the models.This table also includes, as a baseline, a naive method that just returns the last value phin each series as the predicted value p∧h+T.Surprisingly,this type of straightforward mechanism is sometimes quite difficult to outperform.

    Table 2:Metrics values for the evaluated models for different time horizons T(h).The parameters in brackets of each model indicate its input sources.The best value per metric is shown in bold

    Table 2:Continued

    The table includes the results for 3 different hourly horizons (T): 1, 2 and 3 h ahead parking prediction.Obviously, a decrease in accuracy is observed the further we try to estimate, especially in the combinations(model and input)that provide better results.As was highlighted in bold in the table,GRU model is clearly better than the other alternatives in all cases.However,there are certain differences between metrics and scenarios when including or not the weather in the inputs.

    In 1 h-ahead forecasting, when penalizing the points further away from the mean, the best combination of inputs includes previous parking observations and mobility flows: RMSE=20.36 parking spaces and CVRMSE=4.39%.However, when not penalizing it, the best combination includes also the weather information:MAE=15.32 parking spaces and MAPE=4.04%.The choice between the two depends on the objectives of parking spaces.We consider that if we are mistaken by 10 spaces should count more than just two times being mistaken by 5 spaces since precision is very important.Therefore,we would choose not to include the weather.

    In 2 h-ahead and 3 h-ahead the better results are achieved when not considering the weather as input for all metrics,having a CVRMSE=8.5%and CVRMSE=10%respectively.For the 3 h-ahead scenario,the Bi-LSTM model achieved slightly better results than the GRU model.Nevertheless,it is important to remark that this improvement comes at the cost of using a much complex neural network as a bi-LSTM not only uses more complex cells than a GRU model but also its layers include extra connections(parameters)for the in-reverse flow of information as it was put forward in Sec.4.2.3.

    Fig.8 shows the actual available spaces in the target parking lots in the test interval along with the predictions of the generated models for a time horizon T=1 h.As we can see,the actual availability of the target parking lots follows a quite strong daily seasonal pattern.However, this figure also shows that only the GRU model was able to correctly predict the abnormal availability of the parking lots during the interval between the 80 and 120 h.This clearly shows a better matching between the prediction and the actual free parking when using GRU compared to MLP.

    Figure 8: Time series of the available parking spaces (ph) and the prediction outcome of the models(.(a)MLP(I,P,W)(b)GRU(I,P,W)(c)Bi-LSTM(I,P,W)

    Finally, Tab.3 shows a comparison between the results obtained by our solution and previous approaches in the literature for availability prediction in closed parking lots(Sec.2.3).

    From this table, it is important to remark that results in [13–17] are given with respect the occupancy rate of the car parks.If we compare such values with the computed rate-based metrics of our approach(0.0233 MAE and 0.0311 RMSE for 1 h time horizon),we can see that our approach obtained slightly more accurate results.Furthermore,the time horizon is much larger than the baseline approaches as they range from a few minutes to 1 h.Regarding the proposal in[18],the absolute values of the metrics are also quite similar with a much larger time horizon.

    All in all, results state that the usage of the human mobility dataset has clearly improved the accuracy of the predictor.However,such a dataset reflects all types of human displacements without considering the adopted means of transport.Since the parking lots are only used by motorized vehicles,the impact of using the SMT feed also reflects that the human mobility contained in this source is mainly defined by private-car trips.Consequently, this also confirms the fact that such vehicles are still the dominant means of transport for regular trips in most urban settlements.

    Table 3:Comparison of our approach with existing approaches of the literature.In brackets the time horizon of each solution

    5 Conclusions

    The management of road traffic in large cities is one of the most challenging problems for authorities and operators.Since many urban traffic jams are caused by drivers searching for parking spaces,the accurate prediction of available spaces is an instrumental tool to enable efficient road-traffic policies.

    However, most of the existing solutions generally assume that such contextual data is easily available through a large palette of location sensors.This limits the usability of current solutions because the access to location data in most real-world scenarios is restricted due to many regulations.

    In this context,the present work proposes a novel mechanism for parking space prediction that leverages area-based human mobility data as the primary contextual source.The key goal was to evaluate whether such a coarse-grained mobility feed was useful to detect latent movement behaviors in particular locations within a city.To do so,we tested our approach in a real-world setting by making use of different ANNsto predict the available spaces of two downtown parking lots.Results have shown that enriching the models with the flows from the open dataset meaningfully improved their prediction accuracy.

    In the future, with regards to research, our system can be used to feed fine-grained systems to improve their generalisation and also to study general parking patterns between cities to develop transfer the extracted knowledge between the similar ones using transfer learning.With regards to the development of new services,our parking availability prediction can be part of a system that includes route suggestions to drivers according to users’preferences with regards to time,pollution,money,and availability to use other transport services in conjunction with their private car.Also,a new business model can emerge,such as the automatic,hourly,and renting of private parking spaces according to real-time needs.

    Code availability:The source code of the project is available at https://github.com/fterroso/mobility_murcia_umu_ucam.

    Funding Statement:This work has been sponsored by UMU-CAMPUS LIVING LAB EQC2019-006176-P funded by ERDF funds,by the European Commission through the H2020 PHOENIX(grant agreement 893079)and DEMETER(grant agreement 857202)EU Projects.It was also co-financed by the European Social Fund(ESF)and the Youth European Initiative(YEI)under the Spanish Seneca Foundation(CARM).

    Conflicts of Interest:The authors declare that they have no conflicts of interest to report regarding the present study.

    av天堂中文字幕网| 夜夜看夜夜爽夜夜摸| 免费av不卡在线播放| 成年人黄色毛片网站| 偷拍熟女少妇极品色| 久久久国产成人精品二区| 热99re8久久精品国产| 亚洲欧美激情综合另类| 欧美成狂野欧美在线观看| 香蕉av资源在线| 最好的美女福利视频网| 久久久久国产精品人妻aⅴ院| 午夜精品一区二区三区免费看| 99热这里只有是精品在线观看 | av在线老鸭窝| 极品教师在线免费播放| 91麻豆精品激情在线观看国产| 久久久久久九九精品二区国产| 午夜福利在线观看免费完整高清在 | 国产av不卡久久| 大型黄色视频在线免费观看| 一级黄色大片毛片| 成人精品一区二区免费| 可以在线观看的亚洲视频| 99热只有精品国产| 宅男免费午夜| 免费在线观看影片大全网站| 好男人电影高清在线观看| 中文字幕av在线有码专区| 成人高潮视频无遮挡免费网站| 国产亚洲精品久久久久久毛片| 亚洲欧美精品综合久久99| 99热6这里只有精品| 蜜桃久久精品国产亚洲av| 欧美日韩福利视频一区二区| 亚洲不卡免费看| 观看免费一级毛片| 欧美高清成人免费视频www| 亚洲,欧美,日韩| 精品一区二区三区视频在线观看免费| 亚洲欧美日韩无卡精品| 97超视频在线观看视频| 国内毛片毛片毛片毛片毛片| 桃红色精品国产亚洲av| 国产精品伦人一区二区| 十八禁国产超污无遮挡网站| 国产 一区 欧美 日韩| 国内精品久久久久精免费| 亚洲,欧美,日韩| 亚洲一区二区三区色噜噜| 国产真实伦视频高清在线观看 | 女人被狂操c到高潮| 直男gayav资源| 人人妻人人看人人澡| 国产毛片a区久久久久| 婷婷丁香在线五月| 国产伦精品一区二区三区四那| 免费无遮挡裸体视频| 国产精品久久久久久久电影| 搡老岳熟女国产| 国产精品一区二区三区四区久久| 啪啪无遮挡十八禁网站| 日本撒尿小便嘘嘘汇集6| 99在线视频只有这里精品首页| 午夜激情欧美在线| 亚洲av成人精品一区久久| 欧美绝顶高潮抽搐喷水| 成人av一区二区三区在线看| 日本免费一区二区三区高清不卡| 欧美午夜高清在线| 两个人视频免费观看高清| 精品熟女少妇八av免费久了| 色尼玛亚洲综合影院| 好看av亚洲va欧美ⅴa在| 熟女电影av网| 悠悠久久av| 国产免费男女视频| 波多野结衣高清无吗| .国产精品久久| 国产精品亚洲美女久久久| 99久久成人亚洲精品观看| 国产主播在线观看一区二区| 看十八女毛片水多多多| 亚洲自偷自拍三级| 我的女老师完整版在线观看| 99久久久亚洲精品蜜臀av| 成年女人看的毛片在线观看| 国产精品亚洲一级av第二区| 看十八女毛片水多多多| 成熟少妇高潮喷水视频| 51国产日韩欧美| 日韩中字成人| 91av网一区二区| 午夜影院日韩av| 国产亚洲精品av在线| 国产色爽女视频免费观看| 欧美成人一区二区免费高清观看| 免费电影在线观看免费观看| 免费黄网站久久成人精品 | 午夜日韩欧美国产| 亚洲午夜理论影院| 国产亚洲精品av在线| 日日摸夜夜添夜夜添小说| 国产av麻豆久久久久久久| 麻豆一二三区av精品| 亚洲在线自拍视频| 中文字幕人妻熟人妻熟丝袜美| 亚洲第一电影网av| 久久久久久久久久成人| 欧美乱色亚洲激情| 成人三级黄色视频| 91在线精品国自产拍蜜月| 一卡2卡三卡四卡精品乱码亚洲| 久久久久久久久中文| 日本黄大片高清| 变态另类丝袜制服| 51国产日韩欧美| 免费观看精品视频网站| 日本五十路高清| 久久精品影院6| 午夜免费激情av| 免费看日本二区| 精品福利观看| 搡老岳熟女国产| av中文乱码字幕在线| 亚洲熟妇中文字幕五十中出| 亚洲av熟女| 99热这里只有是精品在线观看 | 美女高潮喷水抽搐中文字幕| 毛片一级片免费看久久久久 | 亚洲欧美精品综合久久99| 国产精品野战在线观看| 97人妻精品一区二区三区麻豆| 国产精品影院久久| 亚洲av中文字字幕乱码综合| 亚洲无线观看免费| 欧美激情国产日韩精品一区| 18禁在线播放成人免费| 日本熟妇午夜| 99久国产av精品| 亚洲av电影不卡..在线观看| 波多野结衣高清无吗| 人人妻,人人澡人人爽秒播| 国产乱人伦免费视频| 成人特级黄色片久久久久久久| 亚洲av二区三区四区| 给我免费播放毛片高清在线观看| 超碰av人人做人人爽久久| 搡女人真爽免费视频火全软件 | 成人精品一区二区免费| 俄罗斯特黄特色一大片| 老鸭窝网址在线观看| 国产美女午夜福利| 乱码一卡2卡4卡精品| 成人国产综合亚洲| 此物有八面人人有两片| 一个人看视频在线观看www免费| 国产精品1区2区在线观看.| 亚洲一区二区三区色噜噜| 一进一出好大好爽视频| 男女视频在线观看网站免费| 久久草成人影院| 久久久久久久亚洲中文字幕 | 免费人成视频x8x8入口观看| 成年免费大片在线观看| 在线国产一区二区在线| 波多野结衣巨乳人妻| 免费在线观看影片大全网站| 天堂√8在线中文| 一级a爱片免费观看的视频| 成年女人永久免费观看视频| 亚洲国产色片| 国产又黄又爽又无遮挡在线| 亚洲中文字幕日韩| or卡值多少钱| 亚洲第一区二区三区不卡| 免费看光身美女| 美女 人体艺术 gogo| 国产美女午夜福利| 亚洲av免费高清在线观看| 久久6这里有精品| 美女cb高潮喷水在线观看| 成年免费大片在线观看| 国产av一区在线观看免费| 观看免费一级毛片| 脱女人内裤的视频| 亚洲男人的天堂狠狠| 久久国产精品人妻蜜桃| 午夜视频国产福利| 久久久久国产精品人妻aⅴ院| 国产伦精品一区二区三区四那| 91狼人影院| 日日摸夜夜添夜夜添av毛片 | 老熟妇仑乱视频hdxx| 国产亚洲精品av在线| 啦啦啦观看免费观看视频高清| 一区福利在线观看| 国产极品精品免费视频能看的| av天堂中文字幕网| 动漫黄色视频在线观看| 日韩免费av在线播放| 国产精品影院久久| 国产一区二区在线观看日韩| 亚洲中文字幕一区二区三区有码在线看| 日韩免费av在线播放| 五月伊人婷婷丁香| 久久这里只有精品中国| 久久精品人妻少妇| 欧美日本亚洲视频在线播放| 中文字幕久久专区| av黄色大香蕉| 亚洲成av人片在线播放无| 国产高清三级在线| 亚洲va日本ⅴa欧美va伊人久久| ponron亚洲| 一个人免费在线观看的高清视频| 国产精品亚洲av一区麻豆| 国产精品亚洲美女久久久| 999久久久精品免费观看国产| 亚洲 欧美 日韩 在线 免费| 久久久久久久久久黄片| 午夜精品在线福利| 午夜视频国产福利| av在线观看视频网站免费| 超碰av人人做人人爽久久| 欧美极品一区二区三区四区| 99久久九九国产精品国产免费| 精品一区二区三区人妻视频| 国产成人福利小说| 国产单亲对白刺激| 国产人妻一区二区三区在| 丁香欧美五月| 国产野战对白在线观看| 淫秽高清视频在线观看| 美女免费视频网站| 久久99热这里只有精品18| 久久久久久久久大av| 国产国拍精品亚洲av在线观看| 一个人免费在线观看的高清视频| 一级av片app| 国产高潮美女av| 国产精品自产拍在线观看55亚洲| 夜夜夜夜夜久久久久| 人妻制服诱惑在线中文字幕| 丁香欧美五月| 给我免费播放毛片高清在线观看| 亚洲国产日韩欧美精品在线观看| 天堂√8在线中文| 国产中年淑女户外野战色| 欧美日韩国产亚洲二区| 男人舔女人下体高潮全视频| 非洲黑人性xxxx精品又粗又长| 国产精品电影一区二区三区| 一级a爱片免费观看的视频| 色精品久久人妻99蜜桃| 久久久久免费精品人妻一区二区| 久久精品国产亚洲av香蕉五月| 国产精华一区二区三区| 国产精品电影一区二区三区| 精品久久久久久久久亚洲 | 精品久久久久久久久av| 国产精品久久久久久久久免 | av欧美777| 午夜老司机福利剧场| 高潮久久久久久久久久久不卡| 国产视频一区二区在线看| 国产精品女同一区二区软件 | 波多野结衣巨乳人妻| 国产爱豆传媒在线观看| 午夜免费激情av| 波多野结衣巨乳人妻| 欧美黄色片欧美黄色片| 午夜激情福利司机影院| 久久精品国产亚洲av涩爱 | 午夜精品久久久久久毛片777| 亚洲欧美日韩东京热| 女生性感内裤真人,穿戴方法视频| 男人和女人高潮做爰伦理| 亚洲经典国产精华液单 | 亚洲av电影不卡..在线观看| 成人毛片a级毛片在线播放| 国内毛片毛片毛片毛片毛片| 十八禁国产超污无遮挡网站| 色综合亚洲欧美另类图片| 午夜精品久久久久久毛片777| 亚洲一区二区三区不卡视频| 婷婷六月久久综合丁香| 午夜激情福利司机影院| 脱女人内裤的视频| 欧美激情久久久久久爽电影| 全区人妻精品视频| 国产蜜桃级精品一区二区三区| 欧美性感艳星| 亚洲成av人片免费观看| 国产人妻一区二区三区在| 成熟少妇高潮喷水视频| 午夜免费男女啪啪视频观看 | 欧美日韩国产亚洲二区| www.色视频.com| 国产精品亚洲av一区麻豆| 精品国内亚洲2022精品成人| 亚洲欧美日韩高清在线视频| 久久伊人香网站| 精品一区二区三区av网在线观看| 黄色日韩在线| 神马国产精品三级电影在线观看| 久久久久久久午夜电影| 嫩草影院新地址| 男人狂女人下面高潮的视频| 在线十欧美十亚洲十日本专区| 亚洲自拍偷在线| 亚洲最大成人av| 狂野欧美白嫩少妇大欣赏| 日韩欧美在线乱码| 全区人妻精品视频| 少妇人妻一区二区三区视频| 亚洲狠狠婷婷综合久久图片| 欧美中文日本在线观看视频| 免费看光身美女| 国产一区二区在线观看日韩| 免费人成视频x8x8入口观看| 欧美成人免费av一区二区三区| 在线十欧美十亚洲十日本专区| 搡老岳熟女国产| 欧美在线一区亚洲| 成人亚洲精品av一区二区| 欧美日韩综合久久久久久 | 亚洲熟妇熟女久久| 亚洲欧美清纯卡通| 国产精品亚洲美女久久久| 欧美日本亚洲视频在线播放| 啪啪无遮挡十八禁网站| 亚洲第一电影网av| 麻豆久久精品国产亚洲av| 亚洲av成人av| x7x7x7水蜜桃| 91九色精品人成在线观看| eeuss影院久久| 亚洲第一电影网av| 成人高潮视频无遮挡免费网站| 永久网站在线| 好男人电影高清在线观看| 国产野战对白在线观看| eeuss影院久久| 3wmmmm亚洲av在线观看| 亚洲内射少妇av| 能在线免费观看的黄片| 中文字幕人妻熟人妻熟丝袜美| 国产亚洲精品久久久久久毛片| 精品国产三级普通话版| 国产精品久久久久久人妻精品电影| 亚洲av第一区精品v没综合| 女人被狂操c到高潮| 在线观看舔阴道视频| 老熟妇乱子伦视频在线观看| 亚洲欧美精品综合久久99| 亚洲第一电影网av| 亚洲,欧美,日韩| 久久久成人免费电影| 欧美xxxx性猛交bbbb| 久久精品国产亚洲av涩爱 | 日日摸夜夜添夜夜添av毛片 | 国产亚洲精品av在线| 精品久久久久久,| 国产亚洲精品久久久com| 看片在线看免费视频| 午夜免费成人在线视频| 99国产综合亚洲精品| 免费看美女性在线毛片视频| 最近视频中文字幕2019在线8| 一进一出抽搐gif免费好疼| 久久久久久久午夜电影| 热99re8久久精品国产| 午夜日韩欧美国产| 村上凉子中文字幕在线| 美女xxoo啪啪120秒动态图 | 免费无遮挡裸体视频| 一a级毛片在线观看| 99国产精品一区二区蜜桃av| 国产成人啪精品午夜网站| 变态另类成人亚洲欧美熟女| 国产精品久久久久久精品电影| 夜夜爽天天搞| 亚洲,欧美,日韩| 哪里可以看免费的av片| 看片在线看免费视频| 男女做爰动态图高潮gif福利片| 亚洲人成电影免费在线| 欧美最新免费一区二区三区 | 99热这里只有是精品50| 精品国产三级普通话版| 欧美丝袜亚洲另类 | 啦啦啦观看免费观看视频高清| 亚洲国产欧洲综合997久久,| 99热这里只有精品一区| 亚洲激情在线av| 国产高清三级在线| 日韩欧美一区二区三区在线观看| www.熟女人妻精品国产| 精品久久久久久久久亚洲 | 精品久久久久久久久av| 欧美3d第一页| 91狼人影院| 婷婷色综合大香蕉| 狠狠狠狠99中文字幕| 国产精品永久免费网站| 精品不卡国产一区二区三区| 国产91精品成人一区二区三区| 欧美xxxx黑人xx丫x性爽| 九色国产91popny在线| 在线观看免费视频日本深夜| 麻豆国产97在线/欧美| 美女高潮的动态| 亚洲色图av天堂| 国产精品,欧美在线| 午夜精品久久久久久毛片777| 老司机福利观看| 婷婷亚洲欧美| 欧美黄色片欧美黄色片| 国产中年淑女户外野战色| 757午夜福利合集在线观看| 99热只有精品国产| 亚洲av美国av| 1024手机看黄色片| 国产精品1区2区在线观看.| 日本一二三区视频观看| 午夜精品一区二区三区免费看| 国产乱人视频| 99国产综合亚洲精品| 99热这里只有精品一区| 国产亚洲精品久久久久久毛片| 亚洲内射少妇av| 欧美日韩亚洲国产一区二区在线观看| 此物有八面人人有两片| 久久久色成人| 亚洲专区国产一区二区| 成人午夜高清在线视频| 日日摸夜夜添夜夜添av毛片 | 国产精品免费一区二区三区在线| 成年女人毛片免费观看观看9| 欧美又色又爽又黄视频| 国内精品一区二区在线观看| 午夜亚洲福利在线播放| 白带黄色成豆腐渣| 18禁裸乳无遮挡免费网站照片| 两性午夜刺激爽爽歪歪视频在线观看| 嫩草影院精品99| 国产精品99久久久久久久久| 免费看美女性在线毛片视频| 91在线观看av| 欧美成人性av电影在线观看| 老司机深夜福利视频在线观看| 国产精品自产拍在线观看55亚洲| 日韩大尺度精品在线看网址| 国产一区二区三区在线臀色熟女| 高潮久久久久久久久久久不卡| 亚洲精品日韩av片在线观看| 国产av一区在线观看免费| 久久精品国产自在天天线| 国产视频内射| 婷婷六月久久综合丁香| 蜜桃久久精品国产亚洲av| 美女高潮喷水抽搐中文字幕| 久久国产精品人妻蜜桃| 少妇丰满av| 十八禁国产超污无遮挡网站| 91麻豆精品激情在线观看国产| 最近在线观看免费完整版| 国产精品98久久久久久宅男小说| 亚洲熟妇熟女久久| 丰满人妻熟妇乱又伦精品不卡| 男女床上黄色一级片免费看| 久久久国产成人免费| 中文字幕久久专区| 亚洲真实伦在线观看| 久久天躁狠狠躁夜夜2o2o| 老司机福利观看| 12—13女人毛片做爰片一| 国产伦人伦偷精品视频| 亚洲美女搞黄在线观看 | 女生性感内裤真人,穿戴方法视频| 可以在线观看的亚洲视频| 熟妇人妻久久中文字幕3abv| 欧美黄色片欧美黄色片| 99久久精品国产亚洲精品| 日本成人三级电影网站| 欧美性猛交黑人性爽| 亚洲人成网站在线播| 亚洲中文字幕一区二区三区有码在线看| av在线天堂中文字幕| 日韩精品青青久久久久久| 精品一区二区三区视频在线| 熟女电影av网| 午夜亚洲福利在线播放| 99热精品在线国产| 波多野结衣高清作品| 欧美黄色淫秽网站| 国产午夜福利久久久久久| 男女之事视频高清在线观看| 丁香六月欧美| 午夜a级毛片| 天堂影院成人在线观看| 婷婷色综合大香蕉| 精品无人区乱码1区二区| 精品不卡国产一区二区三区| 最近最新免费中文字幕在线| 小蜜桃在线观看免费完整版高清| 久久精品国产亚洲av香蕉五月| 好男人在线观看高清免费视频| 国内精品一区二区在线观看| 色哟哟·www| www.熟女人妻精品国产| 成人毛片a级毛片在线播放| 日韩欧美免费精品| 中文在线观看免费www的网站| 人人妻人人澡欧美一区二区| 一本久久中文字幕| 看片在线看免费视频| 热99在线观看视频| www.www免费av| 国产精品国产高清国产av| 久久久久久大精品| 九色国产91popny在线| 美女免费视频网站| 亚洲国产色片| 成年版毛片免费区| 51午夜福利影视在线观看| 欧美乱色亚洲激情| 少妇的逼好多水| 午夜福利成人在线免费观看| 午夜福利高清视频| 午夜福利在线观看吧| 亚洲,欧美,日韩| 久久精品影院6| 国产伦精品一区二区三区四那| 变态另类丝袜制服| 婷婷亚洲欧美| 国产麻豆成人av免费视频| 久久精品久久久久久噜噜老黄 | 69av精品久久久久久| av在线蜜桃| 成熟少妇高潮喷水视频| 精品久久国产蜜桃| 国产白丝娇喘喷水9色精品| av视频在线观看入口| 一区二区三区高清视频在线| 99热6这里只有精品| 成人三级黄色视频| 人妻久久中文字幕网| 精品人妻视频免费看| 国产 一区 欧美 日韩| 床上黄色一级片| 午夜免费男女啪啪视频观看 | 国产乱人伦免费视频| 精品人妻熟女av久视频| 亚洲精品一区av在线观看| 一本精品99久久精品77| 精品久久久久久久久av| netflix在线观看网站| 99riav亚洲国产免费| 国产精品一区二区性色av| 最近在线观看免费完整版| 精品人妻视频免费看| 成熟少妇高潮喷水视频| 男女做爰动态图高潮gif福利片| 18+在线观看网站| 天美传媒精品一区二区| 亚洲欧美日韩东京热| 日韩有码中文字幕| av中文乱码字幕在线| 18禁黄网站禁片午夜丰满| 一进一出好大好爽视频| 国产 一区 欧美 日韩| 国产高清激情床上av| 国产高清三级在线| 精品久久久久久久久亚洲 | 精品国产三级普通话版| 97超级碰碰碰精品色视频在线观看| 人妻夜夜爽99麻豆av| 老鸭窝网址在线观看| 免费av观看视频| 国产欧美日韩精品亚洲av| 99热这里只有精品一区| 亚洲美女黄片视频| 日韩 亚洲 欧美在线| 女同久久另类99精品国产91| 我要看日韩黄色一级片| 男女视频在线观看网站免费| 亚洲avbb在线观看| 97人妻精品一区二区三区麻豆| 欧美极品一区二区三区四区| 91九色精品人成在线观看| 波多野结衣高清无吗| 如何舔出高潮| 97超视频在线观看视频| 怎么达到女性高潮| 欧美性猛交黑人性爽| 中文字幕高清在线视频| 国产精品三级大全| 内射极品少妇av片p| 黄色视频,在线免费观看| 久久国产精品人妻蜜桃| 久久精品国产清高在天天线| 一本精品99久久精品77| 丰满乱子伦码专区| 久久久久久久亚洲中文字幕 | 国产精华一区二区三区| 国产极品精品免费视频能看的| 熟女人妻精品中文字幕| 1000部很黄的大片| av女优亚洲男人天堂| 身体一侧抽搐| 青草久久国产| 亚洲 欧美 日韩 在线 免费| 亚洲人成电影免费在线|