Assessment of Presence Conditions of Pavement Markings with Image Processing

Similar documents
A Statistical Method for Eliminating False Counts Due to Debris, Using Automated Visual Inspection for Probe Marks

HOW TO IMPROVE HIGH-FREQUENCY BUS SERVICE RELIABILITY THROUGH SCHEDULING

Airspace Complexity Measurement: An Air Traffic Control Simulation Analysis

Kernel Principal Component Analysis and its Applications in Face Recognition and Active Shape Models

Motion 2. 1 Purpose. 2 Theory

A RECURSION EVENT-DRIVEN MODEL TO SOLVE THE SINGLE AIRPORT GROUND-HOLDING PROBLEM

American Airlines Next Top Model

Transfer Scheduling and Control to Reduce Passenger Waiting Time

Authors. Courtney Slavin Graduate Research Assistant Civil and Environmental Engineering Portland State University

Supplemental Information

According to FAA Advisory Circular 150/5060-5, Airport Capacity and Delay, the elements that affect airfield capacity include:

PREFACE. Service frequency; Hours of service; Service coverage; Passenger loading; Reliability, and Transit vs. auto travel time.

Discriminate Analysis of Synthetic Vision System Equivalent Safety Metric 4 (SVS-ESM-4)

Hydrological study for the operation of Aposelemis reservoir Extended abstract

UC Berkeley Working Papers

1. Introduction. 2.2 Surface Movement Radar Data. 2.3 Determining Spot from Radar Data. 2. Data Sources and Processing. 2.1 SMAP and ODAP Data

You Must Be At Least This Tall To Ride This Paper. Control 27

Clinical results using the Holladay 2 intraocular lens power formula

An Analytical Approach to the BFS vs. DFS Algorithm Selection Problem 1

ARRIVAL CHARACTERISTICS OF PASSENGERS INTENDING TO USE PUBLIC TRANSPORT

METROBUS SERVICE GUIDELINES

ANALYSIS OF THE CONTRIUBTION OF FLIGHTPLAN ROUTE SELECTION ON ENROUTE DELAYS USING RAMS

ANALYZING IMPACT FACTORS OF AIRPORT TAXIING DELAY BASED ON ADS-B DATA

A Multilayer and Time-varying Structural Analysis of the Brazilian Air Transportation Network

Cross-sectional time-series analysis of airspace capacity in Europe

Estimates of the Economic Importance of Tourism

Impact of Landing Fee Policy on Airlines Service Decisions, Financial Performance and Airport Congestion

Quantitative Analysis of the Adapted Physical Education Employment Market in Higher Education

Documentation of the Elevation Selected to Model Helicopter Noise at HTO

Delineation of Urban Freeway Gore Area Crash Cushions in Texas

Time Benefits of Free-Flight for a Commercial Aircraft

EA-12 Coupled Harmonic Oscillators

along a transportation corridor in

4. Serrated Trailing Edge Blade Designs and Tunnel Configuration

A Study on Berth Maneuvering Using Ship Handling Simulator

ScienceDirect. Prediction of Commercial Aircraft Price using the COC & Aircraft Design Factors

White Paper: Assessment of 1-to-Many matching in the airport departure process

CHAPTER 5 SIMULATION MODEL TO DETERMINE FREQUENCY OF A SINGLE BUS ROUTE WITH SINGLE AND MULTIPLE HEADWAYS

Remote Sensing into the Study of Ancient Beiting City in North-Western China

Effective Removal of Pavement Markings

PHY 133 Lab 6 - Conservation of Momentum

An Analysis of Dynamic Actions on the Big Long River

The Computerized Analysis of ATC Tracking Data for an Operational Evaluation of CDTI/ADS-B Technology

SLOPE CALCULATION. Wilderness Trekking School 1

Predicting a Dramatic Contraction in the 10-Year Passenger Demand

The Combination of Flight Count and Control Time as a New Metric of Air Traffic Control Activity

Optimization Model Integrated Flight Schedule and Maintenance Plans

NCUTCD Proposal for Changes to the Manual on Uniform Traffic Control Devices

Provincial Railway Technical Standards

Towards Improved Management of Pavement Markings and Markers. Maurice Masliah, Ph.D., itrans Consulting Inc.

Todsanai Chumwatana, and Ichayaporn Chuaychoo Rangsit University, Thailand, {todsanai.c;

A GEOGRAPHIC ANALYSIS OF OPTIMAL SIGNAGE LOCATION SELECTION IN SCENIC AREA

Northeast Stoney Trail In Calgary, Alberta

Runway Roughness Evaluation- Boeing Bump Methodology

J. Oerlemans - SIMPLE GLACIER MODELS

An Econometric Study of Flight Delay Causes at O Hare International Airport Nathan Daniel Boettcher, Dr. Don Thompson*

Empirical Studies on Strategic Alli Title Airline Industry.

Predicting Flight Delays Using Data Mining Techniques

Digital twin for life predictions in civil aerospace

NOTES ON COST AND COST ESTIMATION by D. Gillen

Pump Fillage Calculation (PFC) Algorithm for Well Control

Photopoint Monitoring in the Adirondack Alpine Zone

Evaluation of High-Occupancy-Vehicle

Sensitivity Analysis for the Integrated Safety Assessment Model (ISAM) John Shortle George Mason University May 28, 2015

OPTIMUM ANCHORING FOR LONG SPAN CABLE-STAYED BRIDGES

15:00 minutes of the scheduled arrival time. As a leader in aviation and air travel data insights, we are uniquely positioned to provide an

SHIP MANAGEMENT SURVEY. January June 2018

Validation of Runway Capacity Models

ScienceDirect. Aircraft parking stands: proposed model for Indonesian airports

Comparison of Arrival Tracks at Different Airports

AIS DATA ANALYSIS FOR REALISTIC SHIP TRAFFIC SIMULATION MODEL

SHIP MANAGEMENT SURVEY* July December 2015

Analysis of en-route vertical flight efficiency

Appendix B Ultimate Airport Capacity and Delay Simulation Modeling Analysis

USE OF 3D GIS IN ANALYSIS OF AIRSPACE OBSTRUCTIONS

1-Hub or 2-Hub networks?

Central Coast Origin-Destination Survey

Polynomial Roller Coaster

Special edition paper Development of a Crew Schedule Data Transfer System

Word Count: 3,565 Number of Tables: 4 Number of Figures: 6 Number of Photographs: 0. Word Limit: 7,500 Tables/Figures Word Count = 2,250

Schedule Compression by Fair Allocation Methods

Simulation of disturbances and modelling of expected train passenger delays

This report was prepared by the Lake Zurich Police Department Traffic Safety Division. Intersection location and RLR camera approaches identified:

HEATHROW COMMUNITY NOISE FORUM

ALLOMETRY: DETERMING IF DOLPHINS ARE SMARTER THAN HUMANS?

MEASURING ACCESSIBILITY TO PASSENGER FLIGHTS IN EUROPE: TOWARDS HARMONISED INDICATORS AT THE REGIONAL LEVEL. Regional Focus.

Unit 6: Probability Plotting

Dr. Dimitris P. Drakoulis THE REGIONAL ORGANIZATION OF THE EASTERN ROMAN EMPIRE IN THE EARLY BYZANTINE PERIOD (4TH-6TH CENTURY A.D.

Bird Strike Damage Rates for Selected Commercial Jet Aircraft Todd Curtis, The AirSafe.com Foundation

AIR TRANSPORT MANAGEMENT Universidade Lusofona January 2008

RE: Draft AC , titled Determining the Classification of a Change to Type Design

Interstate 90 and Mercer Island Mobility Study APRIL Commissioned by. Prepared by

Fixed-Route Operational and Financial Review

Study on impact force calculation formula. of ship lock gravity dolphin

Anomaly Detection in airlines schedules. Asmaa Fillatre Data Scientist, Amadeus

NORTH CASCADE SLACIER CLIMATE PROJECT Director: Dr. Mauri S. Pelto Department of Environmental Science Nichols College, Dudley MA 01571

Statistical Study of the Impact of. Adaptive Traffic Signal Control. Traffic and Transit Performance

Exemplar for Internal Achievement Standard Geography Level 1. Conduct geographic research, with direction

SAMTRANS TITLE VI STANDARDS AND POLICIES

PREFERENCES FOR NIGERIAN DOMESTIC PASSENGER AIRLINE INDUSTRY: A CONJOINT ANALYSIS

Transcription:

Assessment of Presence Conditions of Pavement Markings with Image Processing Yunlong Zhang and Hancheng Ge This paper presents a systematic approach that can automatically determine the presence conditions of pavement markings from marking images developed with digital image processing techniques. These techniques are used to correct geometric deformity, detect colors of pavement markings, segment images, enhance images, detect edge lines of pavement markings, and recognize the features of pavement markings appearing in the photographs. The performance of this system was investigated with the photograph data sets provided by the National Transportation Product Evaluation Program Mississippi test deck. The empirical results (when compared with those from the manual method and expert observations) show that this method produced accurate and reliable results. The techniques in this study can be used to develop an automated system for accurate and speedy assessments of pavement marking conditions and might help agencies make a better decision in the maintenance of pavement markings. Pavement markings play a critical role in maintaining a safe driving environment, especially at nighttime. The Manual on Uniform Traffic Control Devices (MUTCD) () states that pavement markings are commonly placed with paint, thermoplastic, or other materials. Every year state agencies spend more than $ billion in maintaining pavement markings in the United States and Canada. The National Transportation Product Evaluation Program (NTPEP) states that the durability-related performance of pavement markings is evaluated based on the overall percentage of marking material remaining, as well as the retained retroreflectivity. Currently, the most widely adopted method used to determine the presence conditions of pavement markings is expert observation, a subjective technique that may not provide consistent and reliable results for agencies to make cost-effective decisions while satisfying safety requirements. There is a need to develop a new method or system to estimate the presence conditions of pavement markings with higher accuracy and consistency and also reduce labor costs with an automated process. This research is designed to apply the image processing techniques for automatically determining the presence conditions of pavement markings with a higher accuracy, speed, and consistency. Photographs from NTPEP are used to develop components of the image processing process and also for performance assessment. The developed method, as a theoretical foundation, has the potential to be used in an automated system to assess the conditions of any pavement markings. Zachry Department of Civil Engineering, Texas A&M University, TAMU, College Station, TX 77-. Corresponding author: Y. Zhang, yzhang@ civil.tamu.edu. Transportation Research Record: Journal of the Transportation Research Board, No. 7, Transportation Research Board of the National Academies, Washington, D.C., 0, pp.. DOI:./7- Previous Studies In general, pavement marking performance is judged based on two criteria: retroreflectivity and durability. Durability is related to the percentage of pavement marking remaining on a surface and the ability to retain retroreflectivity over time. The study made by Montebello and Schroeder () in 000 focused on providing guidelines for pavement markings in county and city highways. The study found that for roadways with high annual average daily traffic, a more durable product may be a better alternative than paint. Bead application plays an important role in the retroreflectivity of all pavement marking materials. Proper application can lead to increased nighttime visibility and greater line durability. In a study by Cottrell and Hanson (), the effectiveness of pavement markings was investigated. The study examined the durability and the cost-effectiveness of of Virginia s pavement materials. The durability (service life) was estimated by Virginia Department of Transportation s (DOT s) district staff responsible for maintaining the various kinds of markings on Virginia DOT maintained roads. The conclusion was that the large paint contract was the most costeffective for two-lane roads under most traffic volume conditions and for four- and six-lane low-volume roads. In 00, Migletz et al. () evaluated the service life of durable, longer-lasting pavement markings sponsored by FHWA. Six different kinds of markings were evaluated for the durable performance in this study. They established criteria for the minimum retroreflectivity value used to define the end of a pavement marking s service life. Finally, statistical modeling was used to determine the relationship between decreasing retroreflectivity with time (in months) and cumulative traffic passages in freeway and nonfreeway. Gates et al. () studied the determination of the durability of various pavement materials on concrete pavements in Texas. For NTPEP testing, the authors invited Texas A&M Transportation Institute research staff to rate each of the total materials as good, marginal, or poor, based on a combination of the retroreflectivity and durability performance on concrete pavements. They concluded that epoxy materials and preformed tapes should be used on portland cement concrete roadways and thermoplastic should be used only for short-term applications with low to medium traffic. The measurement of the performance of pavement markings is important for agencies around the world to guarantee the effectiveness of markings and determine when the old pavement markings need to be replaced or restored. Visual rating systems, such as the one developed by NTPEP, have been widely used to measure the performance of pavement markings. More recently, advances in image processing techniques have provided an opportunity for the use of this technology in the measurement of pavement markings. Burrow et al. () developed a method to automatically measure road-marking erosion via digital video image analysis techniques.

Zhang and Ge Techniques such as image digitization and image segmentation were used to analyze the image to ascertain whether the picture elements are clearly road markings or not. Furthermore, the authors used the technique of image enhancement to remove misclassified picture elements. Finally, the determination of road-marking erosion can be produced by comparing the ideal road markings based on the feature description. The test results demonstrated that the method developed in this paper was more accurate and reliable than the manual method in the determination of road-marking erosion. A paper written by McCall and Trivedi (7) presented a solution to the problem of detecting road lanes with the steerable filter, which can provide robustness to various environments such as lighting changes, shadows, and road marking variation. According to the test results for a,000-frame image sequence that includes various road conditions, the method developed was shown to have a strong ability of robustness and high accuracy for the detection. In 00, Collado et al. () studied road detection and classification for driver assistance systems, which track several road lanes and identify the type of lane boundaries. From the study, an algorithm was developed that used an edge filter to extract the longitudinal road markings to which a straight line model could be fitted. On the basis of this information, the position, orientation, and type of road lanes can be automatically detected from the camera image. In 00, Noda et al. () developed a method for the recognition of pavement markings on the generated road surface from in-vehicle camera images. The method combined the techniques of projection transformation, image blurring, image clipping, and a patternlearning algorithm. According to the experimental results, the method had a good performance in the recognition of pavement markings and overcame the effect of various environmental factors. Methodology This study used Hough transformation to detect the edges of pavement marking. Because images could be taken with different angles, height, and other settings, camera calibration was used to recover camera settings. Color detection, image segmentation and enhancement, and feature recognition were among the components of signal processing used to determine the presence value of a pavement marking. Hough Transformation for Edge Detection The Hough transformation [first introduced by Hough as cited by Leavers () and expanded to the field of computer vision by Ballard ()] is a powerful technique used to extract features of a particular shape from an image. The main advantage of Hough transformation is a good tolerance of image noise. In this study, Hough transformation was only applied in the detection of edge lines of pavement markings. The simplest Hough transformation is the linear transformation for detecting straight lines. The formula of a straight line (y = k * x + b in x y coordinate system) can be reexpressed as r = x * cos θ + y * sin θ in the polar coordinate system, where r is the distance from the origin to the straight line and θ [0, π] is the orientation of r with respect to the x-axis. Similarly, points (x i, y i ) of the straight line in the x y coordinate system correspond with sinusoidal curves that have an intersect point in the polar coordinate system. The test implementation finds that Hough transformation is capable of detecting edge lines of pavement markings, as shown in Figure. After the process of edge detection in the original image (Figure a), all intersection points (the number of edges that the system has (a) (b) (c) (d) FIGURE Test implementation of Hough transformation: (a) original picture, (b) intersection in polar space, (c) intersection points in -D, and (d) line detection.

Transportation Research Record 7 detected) are found and marked by rectangular boxes in the polar space as shown in Figure, b and c. In the three-dimensional view (Figure c), intersection points are represented by peaks. Figure d demonstrates all detected edges with green lines. Camera Calibration In this study, all photographs provided by NTPEP have oblique distortion because of the angle between the camera and the road plane. A part of marking could overlap with another one in the vertical direction, as the red line illustrates in Figure a. Without the correction of geometric deformity, it would be difficult to detect the exact edges of each pavement marking to segment them in images. Therefore, in order to make such image markings parallel, camera calibration is an essential process for such geometric deformity. A camera s calibration model, first proposed by Schoepflin and Dailey (), recovers the intrinsic and extrinsic camera parameters, such as camera height, focus length, angle between camera and the ground, the aspect ratio, the position of the camera center, and the camera s heading in world coordinates. To simplify the calibration process and enhance the efficiency, it is assumed that the intrinsic camera parameters are fixed (except for focus length). A vanishing point VP(u 0, v 0 ), a prerequisite of camera calibration, is formed by intersecting a pair of parallel edge lines of pavement markings in the image coordinate system. After the vanishing point is acquired, the parameters including focus length f, the pan angle θ, and the tilt angle ϕ for calibrating the camera can be determined. Finally, objective positions in the world coordinate system can be reconstructed by a transformation function. Statistical Color Model for Pavement Markings Color is considered as a visual perceptual property that is crucial in attracting a driver s attention. However, it is challenging to predict which color value belongs to the color class of pavement markings because the color values of pavement markings in images can be widely dispersed under different outdoor conditions. To overcome the effects from various outdoor conditions, the color model is developed to detect pavement marking colors with a backpropagation neural network (BPNN), with a highly accurate and robust detection for pavement markings using labeled pavement markings color samples from NTPEP images. For the basic concept of neural networks, Murino () and Zhang () have authored good introductions that can be referenced. The general structure of BPNN employed in this work is shown in Figure, in which one hidden layer is adopted with 0 neural nodes. H S V High Order Input Terms FIGURE BPNN for color model training (H hue; S saturation; V value; f focus length). f (color) The hue saturation value (HSV) color space, defined by a wide range of colors with representations of hue (H), saturation (S), and value (V), was adopted as the input for the color detection, because it was more perceptual than other color spaces and more appropriate to be used in linear color segmentation. Additionally, it should be noted that the inputs were expanded with the high-order terms to simplify the structure and acquire higher accuracy because multilayer neural networks can be substituted with one-layer networks by using the expanded inputs (). After testing various types of highorder input terms, the high-order polynomials were finally chosen as the expended inputs with a higher accuracy as following: H H H high-order polynomials = S S S V V V () Therefore, the total count of the inputs used in practice was, including three inputs H, S, V and nine inputs of high-order polynomials as shown in Equation. The output of BPNN is a switch value that represents whether a pixel belongs to the color of pavement markings in the image ( means yes, 0 means no). Image Segmentation After color detection, pixels of pavement markings in the image need to be identified and separated from the road surface as the prerequisite for the further determination. Image segmentation is crucial for the whole system, and it makes further analysis of the image only focus on the separated road-marking pixels and hence measures the performance characteristics of the pavement marking itself. Because of such contrast of the brightness between pavement markings and the background in an image, a gray-level value is expected to distinguish these two regions as a threshold value. A variety of techniques have been successfully implemented to segment the images, such as the average method, variance method, and minimum error thresholding method. After investigation, the minimum error thresholding method, first introduced by Kittler and Illingworth (), was adopted because its performance is better than that of the average method and the variance method in image segmentation, even though those methods have a similar perforation in the image that has two obviously distinct peaks (composed of objective and background) in the histogram. Image Enhancement During the process of segmentation, some pixels outside of the region of pavement markings may be incorrectly classified into those of pavement markings. Such incorrectly segmented pixels are known as noise and could significantly affect further analysis of the presence conditions of pavement markings. Because of this potential, there is a need to apply filtering techniques to eliminate such noises. The filtering techniques improve the quality of the image by using such techniques as low-pass filter, histogram modification, median filter, and frequency domain processing, which can further reduce noise and increase the accuracy for further analysis. After testing, the median filter was employed to enhance the image by removing noise. As its name implies, the median filter, the best-known filter in nonlinear spatial filters, replaces the value of a pixel by the median of values in its neighborhood. In this case, the number of white pixels

Zhang and Ge 7 in the neighborhood of a colored pixel will determine whether this pixel belongs to pavement markings or the noise. If the number of white pixels within a neighborhood is lower than a threshold value, this pixel will be considered as a noise and then removed. The opposite holds true, as a pixel will remain if its number of neighboring white pixels is higher than a threshold value. On the basis of Burrow et al. s research (), the optimized size of a neighborhood is * pixels as a window, and the threshold number of white pixels for the determination of noise pixels is set as 0 to acquire maximum efficiency. Feature Recognition The feature recognition in which the quantification of pavement marking erosion will be involved is based on two results: the computation of the number of remaining white pixels within the area between the ideal edge lines of pavement markings in a segmented and enhanced image and the correlation of histographic data with criteria photographs provided by NTPEP as the supplementation. A method is developed to detect ideal edge lines of pavement markings in an image according to the properties of the remaining pavement marking pixels in the image. The key to feature recognition is to estimate the likely position and size of the ideal pavement markings, which are commonly shown as parallelograms in images. The resulting transformation from Hough transformation is usually composed of discontinuous lines that are difficult to be distinguished and unsuitable to determine directly the likely position and size of ideal pavement markings. Fortunately, the slopes of detected edge lines can be provided through Hough transformation, which is useful in correcting parallelogram-shaped pavement markings to rectangular-shaped ones for easily determining the ideal edge lines of pavement marking with higher accuracy. Each pixel in a segmented and enhanced image with the x y coordinate has a corresponding point with the v w coordinate under a spatial relationship. Such a spatial relationship can be mathematically described in Equations and : v = x y tan α ( ) w = y () where α is the angle between the parallelogram-shaped pavement marking in the image and the horizontal line. To investigate further the ideal edge lines of pavement markings, it is necessary to calculate the centroid (known as the barycenter of a plane figure) of remaining pavement marking pixels with the following equations, which could be used to estimate the centerline of the ideal pavement marking: C C x y ( ) xdy x dx = ( ) A ( ) ydx ydy = () A where (C x, C y ) = position of the centroid, A = area of remaining pavement marking pixels (total number of white pixels in image), D y (x) = distance between original point O and pavement marking pixel of segmented and enhanced image in x-axis, and D x (y) = distance in y-axis. Hence the centerline can be located as a horizontal line going through the centroid of the remaining pavement marking pixels in the image. After the centerline is located, the ideal edge lines of pavement markings as boundaries can be estimated based on the position of the detected centerline. Because the pavement marking is separated by the centerline into two symmetric parts, the centerline of each part can be calculated by the same method stated previously. The double distance between the centerline of each part and that of the whole image can be calculated and considered as the distance between the ideal edge lines and the centerline of the whole image. Therefore, two ideal edge lines, as the ideal width of the pavement marking, can be finally determined at an equivalent distance apart from the centerline of the whole pavement marking. The percentage of remaining pavement markings could be further determined from the ratio of the number of remaining white pixels and the total pixels between the two ideal edge lines of the pavement marking. The percentage corresponds to the level as follows: 0% to % is for, % to % is for, % to 7% is for, and so forth. Each % decrease corresponds to one decreased level after. However, before determining the final presence condition of pavement markings, one must account for the correlation of histograms between test and criteria photographs provided by NTPEP. The Pearson correlation coefficient r shown in the following equation is most widely used because it is mainly sensitive to a linear relationship between two variables: r XY Cov( X, Y) E X X Y Y = ( X Y )= = [( µ )( µ )] corr, σxσy σ X σ Y where = n xiyi xi yi ( ) ( () n x x n y y ) i i i i X and Y = histograms of test and criteria images; E = expected value; Cov = covariance; µ x and µ y = mean of the histograms X and Y, respectively; σ X and σ Y = standard deviation of the histograms X and Y; n = number of bins of the histogram; and i = bin index. Each test image will be evaluated by the correlation analysis with all criteria photographs from s to. The level of the test image will correspond to the criteria photograph with the highest correlation value. Thereafter, the final determination of the presence conditions of pavement markings will be conducted by taking into account both the results of the percentage and the correlation analysis. If the highest correlation value is below 0., the percentage is only considered as the determination of the presence conditions without the correlation result. If the highest correlation value is more than 0., the final presence conditions (L) of pavement

Transportation Research Record 7 markings will be determined by a linear regression equation as follows: L = 0. P+ 0. C ( 7) where P is the level corresponding to the percentage of remaining white pixels within the ideal pavement marking and C is the level corresponding to the criteria photograph with the highest correlation. The final presence condition L will be rounded to an integer value as the final determination of the presence conditions of pavement markings. Description of Data The photograph data sets were provided by the Mississippi NTPEP pavement marking test deck located along U.S.-7. There are two sites located in a flat area, including a concrete pavement site and an asphalt site. A total of different pavement marking materials (0 of them were adopted in this study) provided by eight vendors were installed at this test deck in 00. The performance of the pavement markings, including retroreflectivity, durability rating, and color, is periodically monitored from the initial set of readings over years. At the time of measurement, the photograph for each pavement marking was recorded by digital camera as well. These photographs are used to test the effectiveness and robustness of the system developed in this study. At each measurement, the presence conditions of pavement markings were also produced by averaging several experts grades and recorded in an Excel file. An example image of pavement markings transversely installed at the Mississippi test deck is illustrated in Figure. During the evaluation period of more than years, a total of,0 photographs were collected from July 00 to July 00. More than,0 pavement markings (including repeated counts) appeared in these photographs. The photograph resolution is,0 7 in JPEG format. To evaluate the proposed algorithm, photographs including a total of,0 pavement markings (, on asphalt surfaces and,7 on concrete surfaces) were chosen, covering two different types of pavement surfaces and various lighting conditions. FIGURE Example image of Mississippi test deck. Three hundred ninety-one photographs were taken from asphalt surfaces and photographs from concrete surfaces. Additionally, to further test the capability of determining the presence conditions of pavement markings for the proposed system, a total of 0 criteria photographs from NTPEP were adopted in this study as well. The criteria photographs labeled from s to are produced with a variety of pavement surfaces and various outdoor conditions such as lighting, weather, and so forth. Examples of criteria photographs are illustrated in Figure. Results and Analysis To assess the system, with respect to the percentage of remaining pavement markings, the manual method is adopted in this study for comparison with the developed system by means of an estimate of the number of remaining pixels of remaining pavement markings within the ideal area of pavement markings using external observations and Matlab image processing tools. With visual observation, the pixels of remaining pavement markings in the image are summed up. Thereafter, the percentage of remaining pavement markings in the image can be determined by the ratio of the actual area and the ideal area of pavement markings manually. This percentage is considered as accurate for further comparison. The comparison with the percentage of pavement markings determined manually for 0 criteria photographs is demonstrated in Table. A paired t-test is used to investigate the difference between the proposed algorithm and the manual method for each criteria photograph at the % confidence level. The p-value of., which is larger than.0, indicated that there is no statistically significant difference between these two methods. In NTPEP Mississippi test deck photographs, each pavement marking has a grade of the presence condition by averaging grades from multiple experts observations, considered as the ground-truth values for comparison with grades from the developed system. To investigate how well the developed system performed, it is better to start from analyzing the pavement markings for which the developed system in this study cannot work well. Failures happened in some pavement markings that were difficult to extract from the pavement surface because of similar coloring with the background (pavement surface). According to the overall observation in this study as shown in Table, 0.0% of presence condition levels of pavement markings obtained by the developed system are consistent with ones conducted by experts observations, which means that the proposed algorithm does not yield significantly different performance from expert observations. Apparently, the developed system in this study performed worse on pavement markings applied on concrete surfaces than on those applied on asphalt surfaces. However, the difference in success rate did not appear to be significant. Although there are 70 ratings inconsistent with ratings provided by experts observations, these differences are ± in most cases, as shown in Figure. This demonstrates that the error distributions for pavement markings applied to both asphalt and concrete surfaces. In general, tests for pavement markings on concrete surfaces yielded more and larger errors than tests for pavement markings on asphalt surfaces. The error distributions for each level of pavement markings are specifically illustrated in Table, supporting the conclusions stated. Clearly, a ± error accounts for a large proportion of tested errors for

Zhang and Ge FIGURE Examples of criteria photographs. pavement markings graded from to. Sometimes, it is difficult to visually distinguish pavement markings graded at s and for experts as well. On the basis of visual inspection, most of these pavement markings that are inconsistent with the results from the expert observations are graded at or. This is because experts will focus on the horizontal erosion of pavement markings as well as the vertical erosion, such as the change of thickness, the wear of glass beads, and so forth. Nevertheless, because the developed system cannot evaluate the presence conditions of pavement markings in the vertical direction because of the limitations of twodimensional photographs adopted as the data source in this study, some inconsistency will occur. Moreover, it is expected that more errors appeared in the pavement markings on concrete surfaces than asphalt surfaces because of more samples in the former, in addition to the previously stated similarity of colors. The maximum error is ±, which only occurred in two pavement markings applied on the concrete surface with s and 7 graded by experts. Compared with the amount of pavement markings evaluated here, such errors are only a small part of the overall results and are therefore insignificant to the broader picture. Additionally, the visual inspection did not completely support the field ratings. Generally, the error distributions in this case study demonstrate that the determinations produced by the developed system are not greatly dispersed as compared with those made by experts. In this study, as mentioned earlier, there are a total of 0 different products of pavement markings, such as temporary tape, thermo plastic, methyl methacrylate (MMA), waterborne, epoxy, and solvent-borne, which are all used to test the system s performance. According to the levels graded by the expert observations for each pavement marking displayed in test deck photographs, the results revealed that the developed system did not perform well for pavement markings made from materials such as MMA on both asphalt and concrete surfaces. Figure shows the comparison between the levels of presence conditions noted by expert observations and the developed system for four MMA markings located in the NTPEP Mississippi test deck. For all measurements, the levels graded by the developed system are equal to or lower than ones made by the expert observations during the whole period of evaluation. Additionally, visual inspection of the curves indicates that the differences of levels between two algorithms become smaller over time, which implies that the developed system is capable of performing better over time. Meanwhile, after the results of pavement markings on asphalt and concrete surfaces have been examined, it can be shown that the developed system may not be appropriate for evaluating thermoplastic pavement markings in the first years because of overevaluation, such as grading the pavement markings as. Thermoplastic pavement markings, as the most common pavement markings used on roadways in Texas for years, have been

0 Transportation Research Record 7 TABLE Comparison of Percentage of Pavement Markings with Criteria Images Criteria Photograph () Pixels of Remaining Pavement Markings Percentage of Remaining Pavement Markings Manual Method Proposed System,0 7.. 0,..,77..,7. 0.,777..,7..,7.0.,. 7.0,7.70.7,.7.,..0 7 7, 7. 7.,.7 7.,..,0. 0.7,7.7.,7..,..0,.0.,7.. Mean 0.70 0. Variance 0.07 0.07 known to last up to to years depending on traffic volumes. At to years, the erosion of thermo plastic pavement markings occurred mostly in the vertical direction and slightly in the horizontal direction. As stated previously, the developed system in this study only focused on the horizontal erosion of pavement markings, which is the main reason for overgrading thermoplastic pavement markings. The same problem of overgrading the presence conditions can be found in performed thermoplastic markings. However, this may not be a problem because the system is only needed when significant marking wear has occurred, so the issue of grading markings as may not occur or does not matter to maintenance decision making. As shown in the table below, there are several pavement markings, such as durable tape, waterborne, and temporary tape, for which the developed system is capable of better determining the presence conditions because of the high rate of consistency with the levels made by the expert observations. Consistency with s Graded by Expert Observation Pavement Marking Type Quantity/Total (markings) Percentage Durable tape 0/ Waterborne / Temporary tape 7/ Overall the developed system can perform well in the determination of presence conditions of pavement markings such as durable tape, waterborne, and temporary thermoplastic. For MMA and thermo plastic markings, the developed system is capable of determining the presence conditions well or years after the application. Conclusions The goal of this research is to automatically determine the presence conditions of pavement markings with digital image processing techniques. To achieve this goal, a system was developed in this study that consists of five components that can be used to correct the geometric deformity shown in the NTPEP test deck photographs, detect colors of pavement markings, segment images, enhance images, detect the edge lines of pavement markings, and recognize features of pavement markings. To validate the performance of the developed system, different photographs provided by the NTPEP are tested for effectiveness and reliability. With 0 criteria photographs, the developed system is capable of grading the level of presence conditions of pavement markings in the criteria photographs, consistent with the results from the manual method (considered as reliable, but time consuming and labor intensive). With the NTPEP Mississippi test deck photo graphs, the developed system is capable of determining the presence conditions of pavement markings such as durable tape, waterborne, and temporary thermoplastic with high rates of consistency with TABLE Results for Testing Developed System with Test Deck Photographs Asphalt Concrete Total Variable Quantity Percentage Quantity Percentage Quantity Percentage Total images 0.00 0.00 0.00 Total pavement, 0.00,7 0.00,0 0.00 markings (lines) Inconsistency with 7.0. 70. observation results Consistency with observation results,7.0, 77.,70 0.0 Note: Consistency means that the level of presence condition of pavement markings produced by the developed system is exactly equal to the one conducted by expert observation and vice versa.

Zhang and Ge FIGURE Error distributions for testing the developed system. the corresponding levels made by expert observations. Recommendations for future research include extending the work to measure other important characteristics of pavement markings, such as luminance, retroreflectivity, and so on, with two-dimensional photographs or video logs under various outdoor conditions. Moreover, it is recommended that a geographic information system be integrated into the system, which should be able to help agencies manage pavement markings with better decisions. In addition, the developed system can be integrated into a mobile pavement marking retroreflectometer van for more comprehensive and faster evaluation of pavement markings. The system can also be integrated as a component into an image-based system that detects and analyzes other types of traffic control devices such as signs and markings. TABLE Error Distributions for Each of Pavement Markings Error Distribution of Presence Condition of Pavement Marking Asphalt Concrete ± ± ± Total ± ± ± ± Total 7 0 0 7 0 0 7 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 Total 7 7

Transportation Research Record 7 No. 7th Line: MMA 7 made by expert observation made by the developed system Section 7-No. th Line: MMA 7 made by expert observation made by the developed system No. th Line: MMA 7 Section 7-No. 0th Line: MMA 7 No. th Line: MMA 7 Section -No. th Line: MMA 7 No. th Line: MMA 7 Data Collection Intervals MMA Markings on Asphalt Surface Section -No. 0th Line: MMA 7 Data Collection Intervals MMA Markings on Concrete Surface FIGURE Comparisons between two methods for MMA markings. References. Manual on Uniform Traffic Control Devices. FHWA, U.S. Department of Transportation, 00.. Montebello, D., and J. Schroeder. Cost of Pavement Marking Materials. Publication 000-. Minnesota Department of Transportation, Minneapolis, 000.. Cottrell, B. H., and R. A. Hanson. Determining the Effectiveness of Pavement Marking Materials. Publication VTRC-0-R. Virginia Transportation Research Council and FHWA, U.S. Department of Transportation, Charlottesville, 00.. Migletz, J., J. L. Graham, D. W. Harwood, and K. M. Bauer. Service Life of Durable Pavement Markings. In Transportation Research Record: Journal of the Transportation Research Board, No. 7, TRB, National Research Council, Washington, D.C., 00 pp... Gates, T., H. Hawkins, and E. Rose. Effective Pavement Marking Materials and Applications for Portland Cement Concrete Roadways. Report FHWA/TX-0/0-. Texas Transportation Institute, Texas A&M University System, College Station, 00.. Burrow, M. P. N., H. T. Evdorides, and M. S. Snaith. Road Marking Assessment Using Digital Image Analysis. Proc., Institution of Civil Engineers: Transport, Vol., No., 000, pp. 7. 7. McCall, J. C., and M. M. Trivedi. An Integrated, Robust Approach to Lane Marking Detection and Lane Tracking. Proc., IEEE Intelligent Vehicles Symposium, Parma, Italy, IEEE, New York, 00, pp. 7.. Collado, J. M., C. Hilario, A. De La Escalera, and J. M. Armingol. Adaptive Road Lanes Detection and Classification. Lecture Notes in Computer Science, Vol. 7, 00, pp... Noda, M., T. Takahashi, D. Deguchi, I. Ide, H. Murase, Y. Kojima, and T. Naito. Recognition of Road Markings from In-Vehicle Camera Images by a Generative Learning Method. Proceedings of the Institute of Electronics, Information and Communication Engineers (IEICE), Vol., No., 00, pp... Leavers, V. F. Which Hough Transform? CVGIP: Image Understanding, Vol., No.,, pp. 0.. Ballard, D. H. Generalizing the Hough Transform to Detect Arbitrary Shapes. Pattern Recognition, Vol., No.,, pp... Schoepflin, T. N., and D. J. Dailey. Dynamic Camera Calibration of Roadside Traffic Management Cameras for Vehicle Speed Estimation. IEEE Transactions on Intelligent Transportation System, Vol., No., 00, pp. 0.. Murino, V. Structured Neural Networks for Pattern Recognition. IEEE Transactions on Systems, Man, and Cybernetics, Part B: Cybernetics, Vol., No.,, pp... Zhang, Y. X. Artificial Neural Networks Based on Principal Component Analysis Input Selection for Clinical Pattern Recognition Analysis. Talanta, Vol. 7, No., 007, pp. 7.. Pao, Y. H., and Y. Takefuji. Functional-Link Net Computing: Theory, System Architecture, and Functionalities. Computer, Vol., No.,, pp. 7 7.. Kittler, J., and J. Illingworth. Minimum Error Thresholding. Pattern Recognition, Vol., No.,, pp. 7. The Signing and Marking Materials Committee peer-reviewed this paper.