INVESTIGATION OF BICYCLE TRAVEL TIME ESTIMATION USING BLUETOOTH SENSORS FOR LOW SAMPLING RATES

Similar documents
Big Data and High Performance Computing

Mixed Frequency Allocation Strategy for GSM-R

Electrical Engineering

RHT03-SPI Temperature & Humidity Sensor

Chapter 6 Basics of Digital Audio

Slicing of Aspect-Oriented Programs Based on Reduced Dynamic Dependence Graph *

SVPWM 三相电压源型逆变器的死区时间效应可能会导致电压损失, 电流波形畸变和转矩脉动. 为了改善电流波形, 减少转矩脉动, 详细分析了死区时间对输出电压的影响, 并提出了 SVPWM A

Multi-Channel CMOS Front-End IC for Physiological Signal Acquisition 生理信号采集的多通道 CMOS 模拟前端集成电路设计 张金勇王磊于力

Table of Contents. DS177-ZH LUXEON V2 Product Datasheet Lumileds Holding B.V. All rights reserved.

WiFi 雷达 无所不在的感知 清华大学软件学院

WiMAX 技术培训系列 ( 一 ) OFDM 物理层技术初探 MCG 赵昕 2006 年 2 月. All rights reserved 2006, Alcatel Shanghai Bell

Ultrasonic thickness measurement for aluminum alloy irregular surface parts based on spectral analysis

A Self-Adaptive Low-Voltage Current Mode ASK Demodulator for RFID Tags

Analog Applications 模拟应用期刊 Journal

M030 and M040. Operating Manual Genelec M030 and M040 Active Monitoring Loudspeakers 操作手册

US SEALED SERIES. Structure. Features 特点. Connector M Connector F

3-Channel Constant Current LED Driver UCS1903

1 st 5G Algorithm Competition SCMA

Application of generalized estimating equations for crash frequency modeling with temporal correlation *

09-Architecture Design 游戏架构设计. Ref cornell

Chapter 3. Cellular Concept

Game playing. Chapter 6

Advanced Design System. Release 2005A New Features

DETERMINATION OF VEHICULAR TRAVEL PATTERNS IN AN URBAN LOCATION USING BLUETOOTH TECHNOLOGY

downloaded from

1. Travel time measurement using Bluetooth detectors 2. Travel times on arterials (characteristics & challenges) 3. Dealing with outliers 4.

Pantomimes 圣诞节上演的话剧. Look Behind You! 就在你身后!

US SERIES CONNECTOR FOR AUTOMOBILES

YP G-HSD Hall Current Sensor 确保测量准确精度高 安装方便 售价低. 执行标准 Standard 过载能力 Overload 精度 Accuracy 线性度 Linearity 失调电压 Voltage maladjustment

Long-Awaited Film 期待已久的影片

GSM/GPRS. By Mendy Ouzillou Silicon Laboratories Inc. Austin, TX GSM 3GPP (GSM) burst current) GSM900 DCS V

SPECIFICATION. 宁波凯普电子有限公司 Ningbo Kepo Electronics Co.,Ltd.

Brief Introduction to Degree Program in Control. Science and Engineering (0811)

Version:2.0 M1+/M2/M3

Automated impact device with non-synchronous impacts: a practical solution for modal testing during operation *

ADS for your RF Board Design Flow

Stars Snub Twilight Premiere 明星怠慢 暮光之城 首映式

Fabrication of Schottky Diode in Standard CMOS Process

SMART Manufacturing Technologies - A Chinese Approach to Sustainable Manufacturing-

3DD13005ED 主要参数 MAIN CHARACTERISTICS

Air Quality Monitoring with Mosaic: Revisiting sensornets for practical applications. Dr. Wei Dong Zhejiang University

GAMI Newsletter. Foreword. News. News Industry Projects Training Events Contact. News Industry Projects Training Events Contact

1 BGJ 型变频调速纵剪线 ( 液压 ) FOB60000 USD/SET. BGJ variable frequency speed regulating slitting line (hydraulic pressure)

A SMALL PASSIVE UHF RFID TAG FOR METALLIC ITEM IDENTIFICATION. Mun Leng Ng, Kin Seong Leong, and Peter H. Cole

Synchronization Algorithms and VLSI Implementation for DC-OFDM based UWB System

Qualanod Specifications. Edition 1 st January 年 1 月 1 日版

Ansoft HFSS Antenna Design Kit Design Parameters

YUTAKA 1000 SERIES YU-1100 P4-5. YUTAKA1000series

关于台词的备注 : 请注意这不是广播节目的逐字稿件 本文稿可能没有体现录制 编辑过程中对节目做出的改变

NX70 Installation Instructions High Performance Pulse Output Counter Module

High frequency ratio antenna for RFID tags

8 English P.20 한국의 P.21 中 文 P.22

Appendix 6-3: HFSS 3D Excitations

The Advanced Digital Microphone System for Acoustic Guitar USER MANUAL

Qualanod Specifications Qualanod 2017 版标准. QUALISINO AGM, Foshan, Josef Schoppig Secretary General

Si720x 霍尔效应磁性开关 / 锁存器位置传感器数据表

第 21 届中国国际工业博览会. China International Industry Fair 2019

HETERO JUNCTION FIELD EFFECT TRANSISTOR NE4210M01

Measuring the performance of Knowledge Transfer from Universities to Industry in China. ZHONG Wei Renmin Univ

¼'' TS 法律声明 10. MUSIC Tribe

Internal Compact Dual-Band Printed Loop Antenna for Mobile Phone Application

3DD4243D. 封装 Package IPAK

产品规格书 金龙机电股份有限公司 JINLONG MACHINERY & ELECTRONIC CO., LTD. PRODUCT SPECIFICATION. 客户名称 / Buyer Name 客户信息. 客户料号 /Buyer s Part No.

FDTD Modeling of Noise in Computer Packages

SPECIFICATION FOR APPROVAL

How to Write a Quality Technical Paper and Where to Publish within IEEE 如何撰写高质量科技论文以及在 IEEE 进行论文发表

Species identification of sympatric cuckoo nestlings in a multiple-cuckoo system, China

FTP10N40/FTA10N40 400V N 沟道 MOS 场效应管. 产品特点 低的导通电阻 低的栅极电荷 ( 典型值为 34nC) 开关速度快 100% 雪崩测试 符合 RoHS 标准 / 无铅封装 产品应用 高效开关电源 适配器 / 充电器 有源功率因数校正 液晶面板电源 订购代码

JCS1404 封装 Package. High efficiency switch mode power supplies UPS

The Accurate and Compliant Test Solution Provider. Product Catalog

The audience has given us a lot of positive feedback, some typical ones are listed below:

4/17/01. Design Seminar. Agilent EEsof Customer Education and Applications. MESFET Power Amplifier Design: Small Signal Approach

ctbuh.org/papers The City Re-Imagined Title: Paul Friedli, Head of Transit Management, Schindler Max Schwitalla, Founder, Studio Schwitalla

用户手册 HE-XE & HE-XT 100, 102, 103, 104, 105

型号 :ZLE12864B-FFSSWE-NAA

1. Things you need to know 基本词汇 1. Yes, it s great to begin with a word of agreement: Yes 好的 hǎo de (= yes, agree) 是的 shì de (= yes, confirm) 2.

A Comparative Assessment of Crowded Source Travel Time Estimates: A Case Study of Bluetooth vs INRIX on a Suburban Arterial

YDS Presentation. KRL series for Inverter Air-con 横浜電子精工株式会社

Technology. F. Tabatabai, Member, IEEE, H.S. Al-Raweshidy, Senior Member, IEEE. filter, normally do not possess high quality factor at the

V GG -V TO. rout = Saturation Region I OUT. Vout V SAT V GS V OUT. Sheet 1 of 10. Cascode Current Mirror

ASSESSING THE POTENTIAL FOR THE AUTOMATIC DETECTION OF INCIDENTS ON THE BASIS OF INFORMATION OBTAINED FROM ELECTRONIC TOLL TAGS

MPA250 Multi-Purpose Amplifier Operation Manual

Metal Oxide Semiconductor Field Effect Transistor: Additional Concepts

MBD for ANSYS Simulation

CPA 130 Power Amplifier Operation Manual

Antenna Characterization for Bluetooth-based Travel Time Data Collection

操作说明书 Operating Instructions

Design of an Analog Front-end for Ambulatory Biopotential Measurement Systems

Advanced Manufacturing Technology 先进制造技术

激光焊接装置系列. Laser Welder Series

MIPI 技术及物理层测试的挑战. 是德科技 (Keysight) 携手 MIPI 联盟和 Synopsys 共同推动 MIPI 技术发展. Page 1

National Standard of the People's Republic of China 中华人民共和国国家标准

USING BLUETOOTH TM TO MEASURE TRAVEL TIME ALONG ARTERIAL CORRIDORS

Vocia VAM-1. Manual. January Biamp Systems, 9300 SW Gemini Drive, Beaverton, Oregon U.S.A. (503)

Grohe AG v. Zhejiang Jianlong Sanitary Ware Co., Ltd., A Dispute over Infringement of an Exterior Design Patent CHINA GUIDING CASES PROJECT

King Arthur 亚瑟王. A Legendary King 一个具有传奇色彩的国王. Read the text below and do the activity that follows. 阅读下面的短文, 然后完成练习

A2 产品技术规格书 深圳市爱普泰科电子有限公司 ABTEC-A2 音频测试分析仪 前面板图 后面板图 ABTEC-WI099-A0 产品 型号 产品 功能 产品 版本 V1.0 产品 外观

Innovative mobility data collection tools for sustainable planning

Big data in Thessaloniki

EMF Safety and Actions in China Telecom

Transcription:

ZHENYU MEI, Ph.D. E-mail: meizhenyu@zju.edu.cn Zhejiang University Department of Civil Engineering Zijingang Campus, Hangzhou, 38, China DIANHAI WANG, Ph.D. E-mail: wangdianhai@zju.edu.cn Zhejiang University School of Traffic and Transportation Engineering Zijingang Campus, Hangzhou, 38, China JUN CHEN, Ph.D. E-mail: chenjun@seu.edu.cn WEI WANG, Ph.D. E-mail: wangwei@seu.edu.cn Southeast University School of Traffic and Transportation Engineering Sipailou Campus, Nanjing, 296, China Information and Communication Technology Preliminary Communication Submitted: July 21, 213 Approved: Sep. 23, 214 INVESTIGATION OF BICYCLE TRAVEL TIME ESTIMATION USING BLUETOOTH SENSORS FOR LOW SAMPLING RATES ABSTRACT Filtering the data for bicycle travel time using Bluetooth sensors is crucial to the estimation of link travel times on a corridor. The current paper describes an adaptive filtering algorithm for estimating bicycle travel times using Bluetooth data, with consideration of low sampling rates. The data for bicycle travel time using Bluetooth sensors has two characteristics. First, the bicycle flow contains stable and unstable conditions. Second, the collected data have low sampling rates (less than 1%). To avoid erroneous inference, filters are introduced to purify multiple time series. The valid data are identified within a dynamically varying validity window with the use of a robust data-filtering procedure. The size of the validity window varies based on the number of preceding sampling intervals without a Bluetooth record. Applications of the proposed algorithm to the dataset from Genshan East Road and Moganshan Road in Hangzhou demonstrate its ability to track typical variations in bicycle travel time efficiently, while suppressing high frequency noise signals. KEY WORDS Bluetooth sensor; bicycle travel time; low sampling rates; filtering algorithm; 1. INTRODUCTION Bicycles are the most popular travel means in China, which contribute to 3% to 7% of daily trips. Accurate travel time information acquisition is essential to the effective planning and management of bicycle travel conditions. Traditionally, video camera data have been used as the primary source for measuring the quality of bicycle travel time, which have been expensive and invasive to privacy. In recent years, a novel technology that uses bypassing Bluetooth devices to capture vehicle travel times has been proven to be effective [1-3]. With the increased penetration of Bluetooth on the roadways in China, nearly 2% to 3% of bicycles in Hangzhou had active Bluetooth devices in 2. These sensors became a viable alternative to other methods of bicycle travel time data collection [4-6]. Bluetooth is a short-range communications protocol developed by the Special Interests Group for interdevice communications. This protocol creates some challenges for its additional uses, which could be optimized for travel time collection based on Bluetooth MAC (Media Access Control) address matching [7]. Bluetooth sensor systems may provide time stamps at which bicycles pass successive monitoring sensors; however, the data gathered by these systems require some form of filtering to remove outlier readings. Frequently, the outlier data include bicycles making a stop or taking a detour between two locations. These bicycles experience atypical travel times; therefore, these readings should be removed from the dataset of valid readings to avoid producing erroneous travel time estimates. Filtering algorithms are often based on removing the readings with unreasonably low speeds [8-11]. Promet Traffic&Transportation, Vol. 26, 214, No., 383-391 383

However, this approach has several disadvantages. The algorithm does not handle both stable (constant mean) and unstable (varying mean) traffic conditions, especially for bicycles flow. Moreover, the algorithm cannot be applied to low levels of market penetration (less than 1%). In response to the above issues, the current research has been dedicated to the following aspects: 1. Analysis of bicycle travel time error using Bluetooth sensors and problem formulation; 2. Description of the proposed filtering algorithm for travel time estimation using Bluetooth sensors; 3. Development of a filtering algorithm that considers low sampling rates; 4. Verification and implementation of the proposed method to the data from Genshan East Road and Moganshan Road in Hangzhou. The advantages of the proposed method are: 1. The proposed algorithms are simple and do not need a lot of computer resources for their realization; 2. The proposed algorithms can handle both stable (constant mean) and unstable (varying mean) traffic conditions, and correct abrupt changes in the means of bicycle travel times and spurious spikes; 3. Using the obtained coefficients described in the section of considering low sampling rates, the proposed algorithms can be applied to low levels of market penetration (less than 1%); 4. The proposed algorithms can be used not only for offline, but also for online data filtering. 2. PROBLEM FORMULATION Obtaining travel time measurements using Bluetooth devices involves matching an observed MAC address between two locations. The difference in time between the two readings is the travel time [12]. Bluetooth readers are capable of detecting MACs within a specific range; hence, the travel times obtained can be considered zone-to-zone. This is illustrated in Figure 1 where the dashed lines represent the Bluetooth detection zone. The average travel times obtained from the Bluetooth sensors can be expressed as: M BT ^ j, Lh = nl / j = 1 ^t ^Lh- t ^Lhh ja, jb, n L where MBT ^j, Lh is the Bluetooth-based average travel time between nodes A and B during period L, respectively; n is the number of readings; t is the time stamp when a bicycle or Bluetooth device is detected. The Bluetooth-equipped bicycle could be detected any time within the detection zone. The timestamp, t, for the Bluetooth-based method is defined as the time of the first discovery of the Bluetooth device no matter how many times it is detected afterwards. (1) Zone A Point A The authors realized that the travel time reported by the Bluetooth device will be subject to the following sources of error: spatial error, temporal error and sampling error. Spatial error is caused by the variation of detection zone radius. Temporal error is caused by multiple times detection. Sampling error is caused by other mode of transportation being detached [13]. Bluetooth MAC address-based travel times are expected to contain outliers. These outliers are usually caused by drivers who detour or stop at certain places between detector locations rather than being delayed by the bicycle stream. The existing noise removal procedure can be relatively trivial because only bicycles comprise the relatively large travel time data set. However, the Bluetooth device still contains several errors. A median offline filtering algorithm is designed to extract ground truth from the pool of Bluetooth readings. A detailed comparison between the proposed filter and an alternative method by Dion and Rakha [14] can be found in Mei et al. []. 3. PROPOSED FILTERING ALGORITHM The proposed filtering algorithm is the Moving Median Filter (MMF). A window of size L is moving from the beginning to the end of the time series with unit step k and the median of L values falling within the window are continuously calculated. As shown in Eq. (2), the technique estimates the median of L successive travel time values. A similar process is used to estimate the standard deviation in Eq. (3). il = i- L/ 2 i+ L/ 2 M Median^V, f, V h (2) 2 2 v = a^v - h + ^1 - ah $ ^v - h (3) k k 1 Sk 1 Link Figure 1 - Segment composition Zone B Point B where, M il, is the Median of L successive travel time values, where time i is the middle point of window; Vi is the travel time of bicycles at moment i; v k is the Standard Deviation (SD) of time series in the sampling interval k; v Sk is the smoothed SD of time series that is estimated in the sampling interval k with the use of the filter; a is the smoothing factor (varying between and 1, based on the number of valid readings in the current sampling interval). The flowchart for MMF is presented in Figure 2, and a detailed description of the algorithm is introduced below. Step 1: Select time series from database with N elements. 384 Promet Traffic&Transportation, Vol. 26, 214, No., 383-391

Start Input V...V...V 1 i n L=L min, i=l/2 Yes L<L max No Output Failure at point i i=i+1, L=L min Yes i<n-l/2 No M i,l=median(v i-l/2,,v i+l/2) Output modified V...V...V 1 i n Yes M i,l! (P min, P max) No End L=L min L=2*L Yes Abs(V -M ) D i i,l max No V=M i i,l i=i+1 Figure 2 - Flowchart of the proposed MMF Step 2: Start algorithm operation, assign L = Lmin and i = L/ 2. (L is the width of the moving window; Lmin is the minimal width of the window). Step 3: If L < Lmax, continue; otherwise, report program failure at point i, increment i by 1, assign L = Lmin. ( Lmax is the maximal width of the window), Step 4: If i < ^N- L/ 2h, calculate median value M il, for interval 6 i- L/ 2, i+ L/ 2@; otherwise, write filtered time series into database, stop. Step : If MiL!, 6 Pmin, Pmax@, assign L = Lmin, continue; otherwise, enlarge L by factor of 2 and go to step 3.(Maximal and minimal limits of MiL-, Pmax, Pmin ;) Step 6: If absolute value of ^vi- Mi, Lh > Dmax, assign vi = M, i L. (Limit of maximal deviation of current count rate from median value Dmax.) Step 7: Increment i by 1, go to step 3. This filter is optimal for the time series containing abrupt changes in the mean, followed by recovery. The flowchart of the proposed MMF is presented in Figure 2. Figures 3 (a) and (b) illustrate the application of filtering algorithm to survey data with Bluetooth sensor on Genshan East Road, and Figures 3 (c) and (d) illustrate the application on Moganshan Road. Genshan East Road and Moganshan Road are both fourlane arterials with four inside vehicle lanes and two bicycle outside lanes. The average bicycle volumes of two roads exceed 1, bic/h which provides good conditions for detecting bicycle data. Figures 3 (a), (b), (c) and (d) illustrate how the search window enhances the performance of the proposed algorithm. The size of the validity window is defined by the number of standard deviations which is user-definable in the algorithm. However, the basic validity ranges encompassing two or three standard deviations are envisioned to be utilized. The use of two-standard-deviations wide search window would mean that all data points Promet Traffic&Transportation, Vol. 26, 214, No., 383-391 38

(a) (b) 4 4 3 3 2 2 4 4 3 3 2 2 (c) 4 4 3 3 2 2 (d) 4 4 3 3 2 2 Sunday, April, 211 1 2 3 4 6 7 8 9 11 12 13 14 16 17 18 19 2 21 22 23 24 Tuesday, April 12, 211 1 2 3 4 6 7 8 9 11 12 13 14 16 17 18 19 2 21 22 23 24 Moganshan Road, Dengyun Road to Daguan Road Wednesday, April 13, 211 1 2 3 4 6 7 8 9 11 12 13 14 16 17 18 19 2 21 22 23 24 Moganshan Road, Dengyun Road to Daguan Road Thursday, April 14, 211 1 2 3 4 6 7 8 9 11 12 13 14 16 17 18 19 2 21 22 23 24 Figure 3 - Application of MMF filtering algorithm to survey data with Bluetooth sensor 386 Promet Traffic&Transportation, Vol. 26, 214, No., 383-391

within a 9% lognormal confidence interval would be considered valid, and that all other points falling outside this range would be rejected from consideration when estimating average link travel times. The use of a validity window three-standard-deviations wide validity window would mean that all data points within a 99% confidence interval would be considered valid. Through using the four sets of data of two arterials, the figures illustrate the variation in the interval average travel time, as well as the lower and upper limits of the validity window, after the consideration of trends outside the basic validity range. 4. CONSIDERATION OF LOW SAMPLING RATES To further improve the filtering algorithm for the consideration of low sampling rates, additional enhancements are made. Given the variable nature of traffic flow, estimating the travel times during an interval while using data collected in the previous intervals does not ensure that the resulting estimates are truly representative of the real average trip time of the interval [, 11, 16]. For instance, if bicycle flow is slowly increasing during a given portion of the day, the average travel time measured in consecutive intervals could be expected to increase gradually. Secondly, the assumption that the median travel time and standard deviation of the validity window would remain constant during intervals with no readings, as defined in Eqs. (2) and (3), could potentially result in the algorithm using outdated average travel times to determine the validity window limits. As can be observed in Figure 4, the application of the filtering criterion defined by Eqs. (2) and (3) does not produce good results. The filtering criterion is unable to follow the underlying travel time increase between 6:3 and 8:3 a.m. This failure is caused by the combination of the abrupt increase in roadway travel times and no observations. Because travel times change rapidly at the onset of the period, a gap of a few minutes is sufficient in this case to raise the average link travel time beyond the validity bounds. The validity window would remain constant during intervals with no readings. Since all subsequent observations lie outside the validity window, they are considered invalid and rejected. In this case, a period with no recorded travel time does not mean that there are no bicycles passing through the pair of Bluetooth sensors, but simply that no bicycles equipped with Bluetooth are travelling. To increase the responsiveness of the algorithm to the changes in traffic conditions and to situations with low sampling rates Eqs. (4) and () are introduced to modify the search window limits. Eqs. (4) and () compute for the number of standard deviations nv k used in Eqs. (2) and (3) within each sampling interval k, based on the number of intervals with zero readings. D max = nv k $ vk (4) nk n = m+ m6 1-^1-b h @ () vk v where nk is a user-defined parameter; m represents the minimum number of standard deviations to consider, and a sensitivity parameter bv. Eq. () provides a model that dynamically adjusts the size of the validity window based on the number of preceding sampling intervals without Bluetooth record. For any interval in which at least one reading was made in the preceding interval, the equation defines a validity window that corresponds to the minimum size specified by the user, m. If no readings were made in the previous interval, the size of the validity window is increased to m+ m^bvh. If no readings were made in the two previous intervals, the validity window is further increased 2 to m+ m71-^1-b vh A. The size of the validity window will continue to increase with each preceding interval without readings, until the maximum size of 2m is reached. 4 4 3 3 2 2 Monday, May 14, 212 1 2 3 4 6 7 8 9 11 12 13 14 16 17 18 19 2 21 22 23 24 Figure 4 Application of MMF filtering algorithm to survey data - with Bluetooth sensor for low sampling rates Promet Traffic&Transportation, Vol. 26, 214, No., 383-391 387

4 4 3 3 2 2 Monday, May 14, 212 Previously rejected valid 1 2 3 4 6 7 8 9 11 12 13 14 16 17 18 19 2 21 22 23 24 Figure - Effect of low sampling search limits on the operation of the proposed filtering algorithm Figure illustrates the effect of Eqs. (4) and () on the performance of the proposed algorithm. The figure compares various algorithms, considering an arterial road segment with a relatively low travel time sampling rate by applying a version of the filtering algorithm with only Eqs. (2) and (3), with a version that also includes Eqs. (4) and (). There is a noticeable difference in the size of the validity window used by both versions of the filtering algorithm. Identical validity window limits are used by both versions of the filtering algorithm for all sampling intervals, for which travel time readings are made in the preceding interval. Increasing differences are observed for intervals preceded by an increasing number of intervals without readings. In Figure, the effect of the modified filtering algorithm is particularly apparent in the intervals between : and 8:3. Within this period, thirteen travel times were observed, which were concentrated between 24 min. and 3 min. After the detection of the first bicycle passage, the minimum and maximum search window limits for the subsequent sampling interval, : to 6:, are set at 8.66 min. and 24.44 min., respectively, by the MMF algorithm. As time passes without any additional readings, the filtering algorithm based only on Eqs. (2) and (3) maintains a fixed validity window, whereas the algorithm with Eqs. (4) and () gradually alters the size of the validity window to 4.2 min. and 3.62 min. by the time the : to 6: interval is reached. This thirteen travel times concentrated between 24 min. and 3 min. were considered as valid by the modified algorithm with Eqs. (4) and (). A visual analysis of the time series of the recorded travel times appears to indicate that the thirteen travel time readings should be considered according to the Video cameras observations on spot.. ALGORITHM TESTING AND VALIDATION To evaluate the ability of the proposed filtering algorithm to follow the fluctuations in the observed travel times correctly, the algorithms defined by Eqs. (2) and (3), and (4) and () were applied to two arterial roads in Hangzhou, Genshan East Road and Moganshan Road, as shown in Figure 6. The series consists of readings made on two arterial roads. The first is Genshan East Road, between Jianguo North Road and Jichang Road, over a period of two consecutive days in April 211 and one day in May 212. The second is Moganshan Road, between Daguan Road and Dengyun Road, over a period of two consecutive days in April 211. The multiple days readings of two arterial road segments are utilized to illustrate the performance of the proposed algorithm. The link of Genshan East Road is 2.28 km long and has an average bicycle volume of 1,12 bic/h, and the link of Moganshan Road is 2.34 km long and has an average bicycle volume of 1,26 A Jianguo North Road A Dengyun Road Study Route Genshan East Road Study Route Moganshan Road Jichang Road Figure 6 - Data collection sketch map B Daguan Road B 388 Promet Traffic&Transportation, Vol. 26, 214, No., 383-391

bic/h. Bicyclist behaviour was also captured through video cameras. Video cameras acquisition is one of the most reliable travel time measurement methods available and travel time data collected by video camera can be used for evaluating the performance of the Bluetooth-based approach. Figures 7 (a), (b), (c), (d) and (e) illustrate the application results of the proposed filtering algorithm to the travel time series described above. The results shown in each figure are those associated with the calibrated parameters producing the best results. The average confidence interval is about 99.6%. The proposed filtering algorithm is able to track all major changes in the observed travel times correctly. There may be arguments about the validity of some readings; however, these are relatively few and do not appear to affect the underlying trend significantly. In particular, the figures show that the proposed filtering algorithm could effectively respond to abrupt increases and decreases in observed travel times. The ability to track these changes is particularly apparent in the data for April 13 in Figures 7 (a) and (b), which illustrate a significant change in travel times possibly caused by the onset of congestion due to an incident or illegal parking on the road. To successfully demonstrate the ability to respond to the underlying trend, Figure 7 shows, more specifi- (a) (b) 4 4 3 3 2 2 4 4 3 3 2 2 (c) 4 4 3 3 2 2 Wednesday, April 13, 211 1 2 3 4 6 7 8 9 11 12 13 14 16 17 18 19 2 21 22 23 24 Moganshan Road, Dengyun Road to Daguan Road Wednesday, April 13, 211 1 2 3 4 6 7 8 9 11 12 13 14 16 17 18 19 2 21 22 23 24 Moganshan Road, Dengyun Road to Daguan Road Thursday, April 14, 211 1 2 3 4 6 7 8 9 11 12 13 14 16 17 18 19 2 21 22 23 24 Figure 7 - Sample application to corridors for low sampling rates ( m= 4, b v =.4) Promet Traffic&Transportation, Vol. 26, 214, No., 383-391 389

(d) (e) 4 4 3 3 2 2 4 4 3 3 2 2 Thursday, April 14, 211 1 2 3 4 6 7 8 9 11 12 13 14 16 17 18 19 2 21 22 23 24 Monday, May 14, 211 1 2 3 4 6 7 8 9 11 12 13 14 16 17 18 19 2 21 22 23 24 Figure 7 - Sample application to corridors for low sampling rates ( m = 4, b v =.4) (continued) cally, the ability of the proposed algorithm to follow general travel time fluctuations on links with very low sampling rates. As indicated earlier, this is an issue with the existing algorithms. As shown in the Figure, the total number of Bluetooth readings made throughout an entire day on the test arterial segments varies between 74 and 96, which are considered small numbers. As a result of these very low sampling rates, the diagrams in Figure 7 of Genshan East Road in May 14, 212 typically feature large intervals with no readings. In the diagram, there are typically no Bluetooth readings between : and 6:. Aside from these two time points, periods of more than 3 min. without any reading is not uncommon. The presence of these large gaps greatly increases the task of the proposed filtering algorithm to identify underlying low frequency trends from the high frequency white noise, as well as to determine whether travel times observed after a large gap truly represent the existing bicycle flow conditions or should be considered invalid. 6. CONCLUSION Filtering the data for bicycle travel time using Bluetooth sensors is crucial to the estimation of link travel times on a corridor. The current paper describes an adaptive filtering algorithm to estimate link travel times using Bluetooth data. The proposed algorithm overcomes a number of shortcomings in the existing algorithms by handling both stable and unstable bicycle flow conditions and functioning at different levels of market penetration (including less than 1% of the traffic volume). The proposed filtering algorithm determines the range of valid travel time readings within each sampling interval by considering the number of readings in the current and previous sampling intervals, as well as the number of consecutive readings outside the validity range. Applications of the proposed algorithm to the dataset of link travel times collected by Bluetooth sensors on Genshan East Road and Moganshan Road in Hangzhou demonstrate its ability to track the underlying roadway travel times correctly, while suppressing high frequency noise. Despite the successful application of the proposed algorithm to the Bluetooth dataset, further tests are required to validate the findings of the current study. The use of historical data as an additional validation criterion or a means to provide travel time information during long intervals without travel time readings from Bluetooth readers should also be investigated. 39 Promet Traffic&Transportation, Vol. 26, 214, No., 383-391

ACKNOWLEDGMENT The current work is supported by the Zhejiang Provincial Natural Science Foundation of China (No. LY12E82), the National Natural Science Foundation of China (No.982) and the National High- Tech Research and Development Program (863 Program) (No. 211AA132). 梅振宇, 博士, 副教授, 单位 : 中国浙江大学土木工程系王殿海, 博士, 教授, 单位 : 中国浙江大学土木工程系陈峻, 博士, 教授, 单位 : 中国东南大学交通学院王炜, 博士, 教授, 单位 : 中国东南大学交通学院 摘要 对蓝牙传感器采集的自行车旅行时间数据进行过滤对于主干路的连线间旅行时间预测非常重要 该论文提出一种可以适用于处理自行车低样本率的蓝牙数据来预测旅行时间的过滤算法 采集来的自行车蓝牙数据有两大特征, 一是自行车车流数据包括稳定情况和不稳定情况 二是采集的数据为低抽样率 ( 小于 1%) 为了排除错误数据的干扰, 首先提出一个过滤算法来处理多组数据, 该算法通过设计一个鲁棒数据过滤器来控制动态的有效窗口大小, 由此确定有效的数据 该动态窗口的大小变化是基于前面没有蓝牙数据的样本间隔数量 最后, 将方法应用到杭州市的艮山东路和莫干山路, 结果显示其能高效地应对自行车流的波动和数据的冗杂性 关键词 蓝牙传感器, 自行车旅行时间, 低样本率, 过滤算法 REFERENCES [1] Ahmed H, EL-Darieby M, Abdulhai B, & Morgan Y. Bluetooth and Wi-Fi-Based Mesh Network Platform for Traffic Monitoring. Transportation Research Board 87 th Annual Meeting. Transportation Research Board, Washington, D.C., 28. [2] Haghani A, Hamedi M, Sadabadi, KF, Yound S, & Tarnoff PJ. Freeway Travel Time Ground Truth Data Collection Using Bluetooth Sensors. Transportation Research Board 89th Annual Meeting. Transportation Research Board, Washington, D.C., 2. [3] Sharifi E, Hamedi M, & Haghani A. Vehicle Detection Rate for Bluetooth Travel Time Sensors: A Case Study in Maryland and Delaware. Transportation Research Board 89 th Annual Meeting. Transportation Research Board, Washington, D.C., 2. [4] Traffic Congestion and Reliability, FHWA (Federal Highway Administration), U.S. Department of Transportation, Sept. 2. [] Urban Traffic Report, China Transportation Institute. Sept. 2. [6] Shan X. A Research on Urban Bicycle Transportation Rational Ridership and Road Resource Allocation (Doctoral dissertation). Southeast University, China, 27. [7] What Is Traffic Message Channel? http://www.tmcforum.com. Accessed July 29. [8] Brennan TM, Day CM, Wasson JS, Sturdevant JR, & Bullock DM. Assessing Signal Timing Plans for Winter Conditions. ITE Learned Journal of Transportation,Was hington,dc,1(1), 211, p. 9-76. [9] Brennan TM, Ernst JM, Day CM, Bullock DM, Krogmeier JV, & Martchouk M. Influence of Vertical Sensor Placement on Data Collection Efficiency From Bluetooth MAC Address Collection Devices. Journal of Transportation Engineering, 136(12), 211, p. 14 19. [] Bullock D, Haseman R, Wasson J, & Spitler R. Anonymous Bluetooth Probes for Measuring Airport Security Screening Passage Time: The Indianapolis Pilot Deployment. Transportation Research Board 89th Annual Meeting. CD-ROM. Transportation Research Board, Washington D.C, 2. [11] Wasson JS, Sturdevant JR, & Bullock DM. Real-Time Travel Time Estimates Using Media Access Control Address Matching. Institute of Transportation Engineers Journal, 78(6), 28, p.2 23. [12] Malinovskiy Y, Wu Y, Wang Y, & Lee U. Field Experiments on Bluetooth-based Travel Time Data Collection. Transportation Research Board 87 th Annual Meeting. CD-ROM. Transportation Research Board, Washington, D.C., 2. [13] Mei Z, Wang D, Chen J. Investigation with Bluetooth Sensors of Bicycle Travel Time Estimation on a Short Corridor. International Journal of Distributed Sensor Networks, 212, 212. [14] Dion F, & Rakha H. Estimating Dynamic Roadway Travel Times Using Automatic Vehicle Identification Data for Low Sampling Rates. Transportation Research Part B, 4(9), 26, p.74 766. [] Mei Z, & Tian B. Real-Time Travel Time Estimation: Filtering Raw Data in an Automatic Vehicle Identification Setting. In Proceedings of the First International Conference on Transportation Engineering, Chengdu, China, 27, p.34-39. [16] Kothuri SM, Tufte KA, Fayed E, & Bertini RL. Toward Understanding and Reducing Errors in Real-Time Estimation of Travel Times. Transportation Research Record: Journal of the Transportation Research Board, 249, 28, p. 21-28. Promet Traffic&Transportation, Vol. 26, 214, No., 383-391 391