Measuring Accessibility to Various ASFs from Public Transit using Spatial Distance Measures in Indian Cities

Nowadays, accessibility to facilities is one of the most discussed issues in sustainable urban planning. In the current research, two spatial distance accessibility measures were applied to evaluate the accessibility to amenities, services, and facilities (ASFs) from public transit (PT) by walking distance in six Indian cities. The first stage accounts for distance measures using the Euclidean distance with a new methodical approach derived from the built-up area with a spatial resolution of 30 m from Landsat data, and for the network distance method, the actual road distances using OpenStreetMap (OSM) for different threshold ranges of distances were derived. Meanwhile, in the second stage, indicators such as built-up area, network connectivity, and network density with the percentage of ASFs are evaluated and combined for normalization process for ranking the city. The present study assesses the accessibility to various ASFs from PT at city level and explores whether the actual road network access (by measuring distance) in Indian cities is contributing to a high level of accessibility. It adopts a unique approach using statistical tools while assessing both Euclidean and network distances. It models a framework for overall benchmarking in all six cities by ranking them for their accessibility. The results show various scenarios in terms of the rank of cities, which had been strongly affected by distance metrics (Euclidean vs. network) and thus emphasize the careful use of these measures as supporting tools for planning. This facilitates the identification of the local barriers and problems with network access that affect the actual distance. This unique approach can help policymakers to identify the gaps in PT coverage for reaching ASFs. Furthermore, it helps in crucial implementation by strategic planning that can be achieved using these distance criteria.


Introduction
For the past five decades, accessibility has been central to physical planning and spatial modelling. Hansen (1959) defined accessibility as the ease of access to a desirable destination by linking land use and activities with the transport network [1]. Later on, researchers called for a paradigm shift from auto mobility-oriented planning to accessibility-oriented planning [2]. Evidence shows that better access to services by feasible modes with less impact on the environment is an influencing factor in achieving sustainable accessibility [3].
In recent years, public transport has become increasingly focused on improvements in sustainability [4]. Improving access to public transit (PT) can be considered as an effective way to reduce the external cost and negative impacts of motorized vehicles [5]. Better access will encourage people to use more environment friendly transportation modes, including walking and cycling. These features that promote various forms of physical activity (such as walking) can be referred to as

Measuring accessibility: Methods and Techniques:
Measuring accessibility is an essential component in many transport and accessibility plans, as it gives evidence for preparing guidelines and policies while also creating accessibility index [30]. Accessibility measures often focus primarily on one specific mode or a limited set of destination types, such as jobs. Accessibility measures are classified based on two main aspects, i.e., travel behavior and potential availability of opportunities [31]. These are divided into four important categories of active measures: i) gravity-based, ii) topology-based, iii) distance-based, and iv) utilitybased [32]. All these measures are now currently focusing on the growing interest in spatial data technologies and are expected to have many precise results [33]. Studies are now widely carried out using geographical information system (GIS)-based mapping of accessibility to public transit [13,32,34]. Moreover, studies focusing on access to health care and travel time to hospitals are now widespread, and the distance to shops and cultural facilities were found to be an important factor [35].
Measures are described by a set of parameters concerning 1) a population-based spatial unit, 2) an accessibility measure, and 3) a distance measure chosen for measuring accessibility [36]. While assessing geographic accessibility, the use of these criteria may produce different results depending on the distance measure selected. Some researchers explored the geographical accessibility of health facilities by calculating the accessibility to the nearest facility using the shortest Euclidean distance, (extension of Minkowski distance measures), network distance or a combination of distance types [35,36]. Whereas few studies used the choice of Manhattan distance and Euclidean distances [37], but investigations later concluded that comparing Euclidean distance and network distance is more reliable [38,39]. Moreover, considering network distance measures further helps in the decisionmaking process for resource location and construction of new facilities [40].
Furthermore, in various countries, accessibility is tested using multiple accessibility instruments (AI). One such is Spatial Network Analysis for Multimodal Urban Transport Systems (SNAMUTS) in Australia, where it uses geographical information systems (GIS) tools to assess the linkage between the public transport network by clustering land use and visualizing its strengths and weaknesses across the metropolitan areas by spatial mapping. Spatial Network Analysis of Public Transport Accessibility Systems (SNAPTA) in the United Kingdom contributes to identifying barriers and creates an accessibility index to integrate NMT with PT. It was developed with an idea to identify "whether, at a reasonable cost and with reasonable ease, one can reach to services and activities". A few studies, such as Isochrones Maps to Facilities (IMaFa) in Spain and Method for Arriving at Maximus Recommendable Size of Shopping Centre (MaReSi) in Norway, measure accessibility to activities such as shopping centres. The IMaFa tool calculates services in a buffer area using Euclidean distance and travel times along a network performed in GIS. MaReSi calculates services by using walking and bicycle distances from dwelling to shopping center within 1-2 km by road distance to the site [41].
Land Use and Public Transport Accessibility Index (LUPTAI), developed in Australia, aims to quantify access by walking and the public transport network to destinations such as education, health

Measuring Accessibility: Methods and Techniques
Measuring accessibility is an essential component in many transport and accessibility plans, as it gives evidence for preparing guidelines and policies while also creating accessibility index [30]. Accessibility measures often focus primarily on one specific mode or a limited set of destination types, such as jobs. Accessibility measures are classified based on two main aspects, i.e., travel behavior and potential availability of opportunities [31]. These are divided into four important categories of active measures: (i) gravity-based, (ii) topology-based, (iii) distance-based, and (iv) utility-based [32]. All these measures are now currently focusing on the growing interest in spatial data technologies and are expected to have many precise results [33]. Studies are now widely carried out using geographical information system (GIS)-based mapping of accessibility to public transit [13,32,34]. Moreover, studies focusing on access to health care and travel time to hospitals are now widespread, and the distance to shops and cultural facilities were found to be an important factor [35].
Measures are described by a set of parameters concerning (1) a population-based spatial unit, (2) an accessibility measure, and (3) a distance measure chosen for measuring accessibility [36]. While assessing geographic accessibility, the use of these criteria may produce different results depending on the distance measure selected. Some researchers explored the geographical accessibility of health facilities by calculating the accessibility to the nearest facility using the shortest Euclidean distance, (extension of Minkowski distance measures), network distance or a combination of distance types [35,36]. Whereas few studies used the choice of Manhattan distance and Euclidean distances [37], but investigations later concluded that comparing Euclidean distance and network distance is more reliable [38,39]. Moreover, considering network distance measures further helps in the decision-making process for resource location and construction of new facilities [40].
Furthermore, in various countries, accessibility is tested using multiple accessibility instruments (AI). One such is Spatial Network Analysis for Multimodal Urban Transport Systems (SNAMUTS) in Australia, where it uses geographical information systems (GIS) tools to assess the linkage between the public transport network by clustering land use and visualizing its strengths and weaknesses across the metropolitan areas by spatial mapping. Spatial Network Analysis of Public Transport Accessibility Systems (SNAPTA) in the United Kingdom contributes to identifying barriers and creates an accessibility index to integrate NMT with PT. It was developed with an idea to identify "whether, at a reasonable cost and with reasonable ease, one can reach to services and activities". A few studies, such as Isochrones Maps to Facilities (IMaFa) in Spain and Method for Arriving at Maximus Recommendable Size of Shopping Centre (MaReSi) in Norway, measure accessibility to activities such as shopping centres. The IMaFa tool calculates services in a buffer area using Euclidean distance and travel times along a network performed in GIS. MaReSi calculates services by using walking and bicycle distances from dwelling to shopping center within 1-2 km by road distance to the site [41].
Land Use and Public Transport Accessibility Index (LUPTAI), developed in Australia, aims to quantify access by walking and the public transport network to destinations such as education, health shopping, jobs. LUPTAI considers walking in two ways: whether it is a single mode of access for a destination considering walking distances, i.e., (400, 600, 800, 1000 and 1200 m) or a mode of access to public transport i.e., (0-20, 20-40 and 40 plus min) [42]. Tools like Walk Score ® (www.walkscore.com) developed in the U.S.A tests access to the nearest public transit and services such as shopping and educational resources (schools and libraries) by walk scores and this tool is available in the public domain [6]. All these instruments are known for interventions that cover the aspects of physical distance, travel time and travel cost.
By comparing all these AI tools irrespective of their strength and weakness, a study showed that instruments using spatial separation measured by an existing network rather than gravity measures have a high level of communication value on ground level, meaning that even non-expert users can easily interpret network flows or link connectivity [43]. Therefore, the presented study adopts spatial accessibility measures used by researchers in the past [36,37,44]. However, these studies have not often been conducted in developing nations like India. Hence, there is little quantitative work on the accessibility of Indian cites done previously. This study uses spatial data for developing an index to compare accessibility in Indian cities, by excluding cost and time factors, as they create complexity for comparing various cities. Here, the distance from ASFs is investigated using both Euclidean and network distance, and we compare both techniques within and among the cities in the Indian context.

Selection of Indicators
In a few studies, the macro-scale evaluation considered land use and land cover, job density and population density [44,45]. However, Indian demographic information corresponds to the 2011 census, and authoritarian land use plans are more than a decade old. This study, therefore, only includes built-up area extracted from Landsat 8 images, and thus the ASFs mapped are from year 2018-2019.
A comprehensive spatial accessibility assessment would help policymakers to take appropriate actions for improving overall accessibility in cities. Relevant indicators are selected from various literature regarding mobility/walkability [46][47][48], for land cover [49], and for network density, connectivity [50,51], as shown in Table 1. The selection of diverse and multifold opportunities will provide a viable option to assess accessibility plans in a much better way. Education, health care, leisure, entertainment activities, and commuting capabilities are commonly seen as the leading demands of urban contemporary life. Unlike economic development infrastructures, social infrastructure includes the facilities and places that primarily serve the specific needs of everyday life [52]. With this understanding, various utilities in developed nations were identified as relating to accessibility to a specific infrastructure, such as health [39,52], parks and green spaces [38,53], schools [54,55], transit [32,[56][57][58]. A common goal for these studies was to limit the usage of private vehicles and encourage access to public transit by walking and bicycling. Based on this, various sectors are categorized as ASFs (Table 2). In India, as per guidelines, policies and a few research studies have identified recreation activities as "Amenities" in ISPRS Int. J. Geo-Inf. 2020, 9, 446 6 of 36 the urban local bodies act, (Government of India, 1992) [59], while police stations and post offices are considered in the category of "Services" [60]. Education and health are considered in the category of "Facilities" in the 83rd Constitutional Amendment Act (Government of India, 2002), which focuses on universal access [61].

Identifying Reliable Walking Distance
To measure accessibility using distance, researchers used standard walking distances similar to the planners' thumb rule, i.e., 400 and 800 m [62,63]. Even the World Bank suggests that bus services should be made available for less than 400 m distance to encourage pedestrians [64]. A few other recommended distances from different literature reviews are shown in Table 3. Japan 800, 1600 m Resident (neighbourhood) [66] Glasgow, Scotland 400 m (high-level access) Resident [67] Athens, Greece 500 m < (high-level access) Resident [68] Melbourne and Perth Australia 400 m (high-level bicycle accessibility) Resident [69] Catania, south Italy 300-600 m Recreational area [39] Nagpur, India 400 m (walking Distance) PT [70] Ahmedabad, India 400 m (walking Distance) PT [71] Indian policies which encourage accessibility specify ideal distances to access urban facilities and services. The Urban and Regional Development Plans Formulation and Implementation-2015 (URDPFI) guidelines [72] dictate that the ideal distance for facilities and transit must be available within 600 to 800 m. The colour code used for mapping and identifying Euclidean and network distances is shown in Table 4.

Methodology
This paper introduces a framework to assess accessibility at a macro level with the help of derived indicators from literature studies for various ASFs from PT, as shown in Figure 2.

Methodology
This paper introduces a framework to assess accessibility at a macro level with the help of derived indicators from literature studies for various ASFs from PT, as shown in Figure 2.

Selection of Cities in India
India's urban population continues to grow at an unprecedented rate. The majority of new development is occurring in major cities such as the metropolitan (over 4 million inhabitants) and million-plus cities (over 1 million inhabitants), popularly known in India as tier-2 cities [73]. As per Wilber Smith Associates' report, two types of accessibility indices are developed: (i) Public Transport Accessibility Index and (ii) Service Accessibility Index. Public Transport Accessibility Index is formulated as the inverse of the average distance (in km) to the nearest bus stop/railway station (suburban/metro). The Service Accessibility Index is computed as the percentage of work trips accessible within 15 to 30 minutes for each city. For both the indices, medium category cities, i.e., those that fall in Tier 2 cities, have better accessibility compared to low and high category cities [74,75]. However, even this study only addressed a few tier-2 cities, as shown in Figure 4. Moreover, studies in India are very rarely conducted specifically for accessibility to each ASF in a broader perspective and how it influences the city's accessibility.
Previously in India, studies conducted at city level focused on accessibility assessment to jobs through public transport. These studies are mostly conducted in tier-1 cities, or cities with a population over 3 million as per the 2011 census, such as Mumbai, Delhi, Kolkata, Chennai, Hyderabad, Bengaluru, Ahmedabad, Surat, and Pune. Only a limited number of studies have focused on issues of access by walking to public transit at the city level to educational facilities [70,76]. However, these studies were limited for tier-2 cities. Most of the Indian tier-2 cities with a population of 1.5 to 3.0 million have an average trip length of less than 6km and the urban radius (which is determined as population distribution by distance from the center of the city) is within 15 km from the center, as shown in Figure 3a,b [77,78].
perspective and how it influences the city's accessibility.
Previously in India, studies conducted at city level focused on accessibility assessment to jobs through public transport. These studies are mostly conducted in tier-1 cities, or cities with a population over 3 million as per the 2011 census, such as Mumbai, Delhi, Kolkata, Chennai, Hyderabad, Bengaluru, Ahmedabad, Surat, and Pune. Only a limited number of studies have focused on issues of access by walking to public transit at the city level to educational facilities [70,76]. However, these studies were limited for tier-2 cities. Most of the Indian tier-2 cities with a population of 1.5 to 3.0 million have an average trip length of less than 6km and the urban radius (which is determined as population distribution by distance from the center of the city) is within 15 km from the center, as shown in Figure 3 (a, b) [77,78].  [77,78]. .
For the current study, large cities with a population of over 3.1 million were ignored, as these cities exhibit unique characteristics in comparison with the million-plus cities in India. Moreover, limited studies were conducted for the shortlisted cities, on the accessibility of PT (PT in this study considers only bus stops). Therefore, confining this scope, the present study attempts to assess the accessibility to various ASFs from bus stops at the city level. All these cities were shortlisted based on population size and presence of a public transportation facility in the city as shown in Figure 4. For this paper, out of the 11 shortlisted cities, only 6 were considered (Visakhapatnam, Nagpur, Bhopal, Jaipur, Lucknow, and Indore), as shown in Figure 4 with blue colour bar, and they are mapped accordingly in Figure 5. For the current study, large cities with a population of over 3.1 million were ignored, as these cities exhibit unique characteristics in comparison with the million-plus cities in India. Moreover, limited studies were conducted for the shortlisted cities, on the accessibility of PT (PT in this study considers only bus stops). Therefore, confining this scope, the present study attempts to assess the accessibility to various ASFs from bus stops at the city level. All these cities were shortlisted based on population size and presence of a public transportation facility in the city as shown in Figure 4. For this paper, out of the 11 shortlisted cities, only 6 were considered (Visakhapatnam, Nagpur, Bhopal, Jaipur, Lucknow, and Indore), as shown in Figure 4 with blue colour bar, and they are mapped accordingly in Figure 5.
The average mean distance between bus stops in all the selected six cities was almost below 500 m, as calculated using the average neighborhood tool in ArcGIS 10.5.1. This is the reliable distance between bus stops (<500 m) as per the World Banking Group [64]. The shape of the city is identified from a study wherein it considers dynamic growth, density and development patterns along with the road network [77], and is mentioned in Table 5.  The average mean distance between bus stops in all the selected six cities was almost below 500 m, as calculated using the average neighborhood tool in ArcGIS 10.5.1. This is the reliable distance between bus stops (<500m) as per the World Banking Group [64]. The shape of the city is identified from a study wherein it considers dynamic growth, density and development patterns along with the road network [77], and is mentioned in Table 5.

Data Collection:
Secondary and spatial data were collected from various sources, as shown in Table A1 in Appendix A. Data are mapped using OpenStreetMap (OSM) in Google Earth (GE) and imported into ArcGIS 10.5.1 for analysis, as shown in Figure 6. For built-up area, the Landsat 8 image of datum WGS84 was used shown in Table 6. Bands 4, 3 and 2 were combined to create true-colour composite images. For analyzing the land cover, supervised classification was performed by a pixel-based supervised classifier, in which other land covers such as forest, agriculture, and water bodies were excluded from this study and only built-up area was delineated and used. The data from OSM in GE were overlapped with built-up area, as accessibility is crucial for identifying geo-located features and labels. The combination of earth observations with OSM features is a tool for large maps and offers significant opportunities for the creation of ideas [81]. Table 6. Extracted Landsat 8 imagery data.

Data Collection
Secondary and spatial data were collected from various sources, as shown in Table A1 in Appendix A. Data are mapped using OpenStreetMap (OSM) in Google Earth (GE) and imported into ArcGIS 10.5.1 for analysis, as shown in Figure 6. For built-up area, the Landsat 8 image of datum WGS84 was used shown in Table 6. Bands 4, 3 and 2 were combined to create true-colour composite images. For analyzing the land cover, supervised classification was performed by a pixel-based supervised classifier, in which other land covers such as forest, agriculture, and water bodies were excluded from this study and only built-up area was delineated and used. The data from OSM in GE were overlapped with built-up area, as accessibility is crucial for identifying geo-located features and labels. The combination of earth observations with OSM features is a tool for large maps and offers significant opportunities for the creation of ideas [81].

Data Analysis
The analysis aims at developing a step by step assessment method. Initially, the quality and accuracy of data were checked by testing them for one city and secondly data processing and mapping were carried out by measuring accessibility by Euclidean and network distance for various indicators in six cities. From the literature, network connectivity, network density and detour index were calculated using the relevant equations. Lastly, the cumulative accessibility score (index) was calculated by combining the relevant indicators through the process of normalization.

Data Quality and Accuracy
Research quality may be defined by the reliability and credibility of information and data. Precision and error-free data are the key characteristics of data accuracy [82]. This study uses builtup area extracted from Landsat 8. Initially, the study examines the accuracy assessment of land cover classification using GE with Nagpur city for the year 2018. The supervised classification was used to classify the image of Nagpur city for projected coordinate system of "WGS_1984_UTM_Zone_44N". Under land cover categories, built-up and non-built-up were classified (refer to Figure A1 in Appendix A). After classification of land cover types, 20 random points were generated in ArcGIS 10.5.1 and random points were converted to KML for opening in GE. Each random point's value was verified from GE for accuracy assessment. GE was used to measure the number of correctly classified ground truth pixels, as shown in Equations (1) and (2). Since GE has a high resolution, it is an important source of information, especially in an urban area where the pattern of land cover is a complex mosaic of different land uses [83]. The source of data is open, and therefore, a manual visual interpretation approach on the GE interface was carried out for 26 November 2018. The result shows that the overall accuracy of land cover for 2018 is 80% and Kappa (K) is 0.84, which is acceptable in both overall accuracy and Kappa accuracy (Refer figure A2 in Appendix A).
Secondly, data acquired from various government sources were mapped using OSM in GE. To check accuracy, ground truthing was conducted randomly in the city of Nagpur for three areas in the year 2018 with a buffer ring of 800 m from the main bus stop in each area, as shown in Figure A3 in Appendix A. In general, the Kappa coefficient is performed to understand the plots' mismatch and is mostly used in land cover classification, which was done initially. In this case, the data point for schools, playgrounds and recreation etc. was vpresumed in OSM, and ground truthing was done by directly visiting the site and checking the existence of ASFs by using GPS coordinates. Mismatch is mainly seen in playground and recreational areas, where it is verified by reaching out to the precise

Data Analysis
The analysis aims at developing a step by step assessment method. Initially, the quality and accuracy of data were checked by testing them for one city and secondly data processing and mapping were carried out by measuring accessibility by Euclidean and network distance for various indicators in six cities. From the literature, network connectivity, network density and detour index were calculated using the relevant equations. Lastly, the cumulative accessibility score (index) was calculated by combining the relevant indicators through the process of normalization.

Data Quality and Accuracy
Research quality may be defined by the reliability and credibility of information and data. Precision and error-free data are the key characteristics of data accuracy [82]. This study uses built-up area extracted from Landsat 8. Initially, the study examines the accuracy assessment of land cover classification using GE with Nagpur city for the year 2018. The supervised classification was used to classify the image of Nagpur city for projected coordinate system of "WGS_1984_UTM_Zone_44N". Under land cover categories, built-up and non-built-up were classified (refer to Figure A1 in Appendix A). After classification of land cover types, 20 random points were generated in ArcGIS 10.5.1 and random points were converted to KML for opening in GE. Each random point's value was verified from GE for accuracy assessment. GE was used to measure the number of correctly classified ground truth pixels, as shown in Equations (1) and (2). Since GE has a high resolution, it is an important source of information, especially in an urban area where the pattern of land cover is a complex mosaic of different land uses [83]. The source of data is open, and therefore, a manual visual interpretation approach on the GE interface was carried out for 26 November 2018. The result shows that the overall accuracy of land cover for 2018 is 80% and Kappa (K) is 0.84, which is acceptable in both overall accuracy and Kappa accuracy (Refer Figure A2 in Appendix A).
Secondly, data acquired from various government sources were mapped using OSM in GE. To check accuracy, ground truthing was conducted randomly in the city of Nagpur for three areas in the year 2018 with a buffer ring of 800 m from the main bus stop in each area, as shown in Figure A3 in Appendix A. In general, the Kappa coefficient is performed to understand the plots' mismatch and is mostly used in land cover classification, which was done initially. In this case, the data point for schools, playgrounds and recreation etc. was vpresumed in OSM, and ground truthing was done by directly visiting the site and checking the existence of ASFs by using GPS coordinates. Mismatch is mainly seen in playground and recreational areas, where it is verified by reaching out to the precise location. Where a few locations are rarely inconsistent, and when one or two locations are found, these are drawn into other services in the evaluation process. The results are shown in Table A2 in Appendix A for the individual area, and the overall accuracy within all the areas is shown in Figure A4 in Appendix A. The accuracy results confirm that the overall accuracy is above 90.0% and Kappa (K) is above 0.85 in all the cases. Studies showed that accuracy higher than 75% is reliable and acceptable [83]. Previously, a study in Nagpur conducted on validation georeferenced Google Earth imagery (GEI) with ground truthing for recreation and playgrounds had a confirmed 95% accuracy, with a kappa value of 0.93 (i.e., 93%), which followed the same method [84]. Other studies in Indian cities compared Landsat 8 images with GE and OSM data, resulting in 87 percent total accuracy in large scale land utilization [85]. This shows that OSM data in Indian cities are sufficiently reliable for analysis.
Total Accuracy = Number of correct plots total number of plots (1) K observed accuracy = observed accuracy − chance agreement 1 − Chance agreement (2)

Data Processing
Initially, the distance from bus stops to the grids is considered and as per various levels of accessibility indices distance, colour-coding is assigned to the grid. Furthermore, this grid is overlapped with the existing built-up area (This is done to categorize levels of accessibility as per distance), as shown in Figure 7, for the city of Bhopal. A similar process was adopted for the other five cities and the percentage of built-up area is observed in Figure 8. A high association of built-up area to PT will help in reducing travel distance and vice versa [4].
Initially, the distance from bus stops to the grids is considered and as per various levels of accessibility indices distance, colour-coding is assigned to the grid. Furthermore, this grid is overlapped with the existing built-up area (This is done to categorize levels of accessibility as per distance), as shown in Figure 7, for the city of Bhopal. A similar process was adopted for the other five cities and the percentage of built-up area is observed in figure 8. A high association of built-up area to PT will help in reducing travel distance and vice versa [4].

Measuring Accessibility:
Two measures of accessibility are selected (i) Euclidean distance and (ii) network distance-both of these measures have a single component. The first measure is proximity distance, while the latter one is actual distance calculated using the existing ground network. Furthermore, network connectivity, density and detour index are analyzed accordingly.

Measuring Accessibility
Two measures of accessibility are selected (i) Euclidean distance and (ii) network distance-both of these measures have a single component. The first measure is proximity distance, while the latter one is actual distance calculated using the existing ground network. Furthermore, network connectivity, density and detour index are analyzed accordingly.

Euclidean Distance Measure Using a Near Tool
Near-tool is used to calculate the distance in ArcGIS 10.5.1. This tool calculates the distance from each point in the class to the next point in a different feature class (ESRI, 2018). Here, the origins are ASFs and destinations are bus stops. The process is based on the Pythagoras theorem Equation (3).
from Equation (3) a = shortest distance between two points (ASF and bus stop) x1, y1 are the geographic coordinates of the centroid of each ASFs. x2, y2 are the geographic coordinates of each bus stop. Euclidean distance accessibility levels are considered up to 2400 m; above this range of distance, bus stop facilities are considered not available. Initially, the distance from the bus stop is measured to the nearest ASFs. For example, in Figure 9, for hospitals and clinics, the distance is measured from the bus stops. Depending on the range of distance, every hospital and clinic will be designated a color and symbol, intersecting the built-up area with the same colour ( Figure 7). Similarly, this process is repeated for the remaining sectors of ASFs in all the six cities. As all the ASFs are represented as a point, moreover for recreation and playgrounds the points are considered at the entrance. Here, the ASFs percentage for very high access range (0-400 m) is shown in Table 7, and the total percentage of ASFs combined (average of recreational places, playgrounds, cinema halls, other services, hospitals and clinics, schools and colleges), noted for Euclidean distance for entire range (0-2400 m) in all the six cities, is shown Table 8.   . Spatial evaluation of hospitals and clinics from bus stops using Euclidean distance in six Indian cities.   . Spatial evaluation of hospitals and clinics from bus stops using Euclidean distance in six Indian cities.

Network Distance
A network layer was first created in ArcGIS 10.5.1. Using the Network Analyst tool, the origins and destinations are determined. The shortest part to reach the bus stops from ASFs is measured in network analysts. The ArcGIS network analyst tool makes use of Dijkstra's algorithm to find the shortest path between two points along a network. One-way restrictions are ignored, while U-turns are permitted at junctions because no mode was considered, and no barriers were assumed restricting movement along the roads. The ArcGIS 10.5.1 Network Analyst allows all points that are used in the analysis to snap to the nearest orthogonal position and measure the distance. Different colour codes are given based on distances from 400 to 2400 m.
In this case, the OSM's network data are more accurate and up-to-date relative to the available urban development plans in selected cities. A topology analysis tool in ArcGIS 10.5.1 is used for necessary inspections and modifications to OSM's road system. Accordingly, it modified pseudo nodes and dangles errors.
The distance was calculated from the closest bus stop to ASFs. Here, for example, hospitals and clinics are mapped by considering the entire network distance to a bus stop and assigned a different colour according to the range of distance, as shown in Figure 10. Similarly, distance is measured for the remaining sectors of ASFs in all the six cities. Here, the ASFs percentage for very high access range is shown in Table 9, and the total percentage of ASFs combined is noted for network distance in all the six cities, as shown in Table 10.
In all the cases of Tables 8 and 10 for Euclidean and network distance, it is observed that over 68% of ASFs are in the moderately high level of accessibility. Thus, network connectivity and density have been examined for the very high accessibility range 0-400 m and moderate-high range 0-1200 m. The entire city network is examined to identify the network condition affecting the percentage of accessibility.
The distance was calculated from the closest bus stop to ASFs. Here, for example, hospitals and clinics are mapped by considering the entire network distance to a bus stop and assigned a different colour according to the range of distance, as shown in Figure 10. Similarly, distance is measured for the remaining sectors of ASFs in all the six cities. Here, the ASFs percentage for very high access range is shown in Table 9, and the total percentage of ASFs combined is noted for network distance in all the six cities, as shown in Table 10.     9. Spatial evaluation of hospitals and clinics from bus stops using Euclidean distance in six Indian cities.

Network Connectivity and Density
"The connectivity of a network may be defined as the degree of completeness of the links between nodes" [86]. The greater the degree of connectivity within a road network, the more efficient that system is. While examining any road network, one should investigate connectivity associations with links and nodes. Studies have shown that in response to the addition of any new links in the future, or the improvement of the existing links, changes to the structure of a road network will have a greater impact on physical accessibility. If there are a greater number of links connecting nodes, the chances of having the shortest part influences accessibility in a better way [51]. Considering the same criteria, this study uses a beta index to identify which city has an efficient ratio of links to nodes that influences connections between ASFs and bus stops.
The beta index is expressed as shown in Equation (4), where E is the total number of edges and V is the total number of vertices in the network. A highly connected transport network will have a high value, and a simple network will have a low value [50]. As the number of edges (links) increases, the connectivity between the vertices (nodes) rises, and the beta index changes progressively.
Here, in this study, the entire city's edges and vertices are calculated for distance range of 0-400 m and 0-1200 m from bus stops to the network directly. This is done using buffer tool in ArcGIS 10.5.1. The ratio of total edges to vertices is calculated in the segments of 0-400 m and 0-1200 m and likewise for the entire city within the flat buffer area. The edges to vertices ratio is calculated overall and not at a single junction, as shown in Figure 11, for Bhopal city, and similar analysis is done for other cities as summarized in Table 11.  Similarly, network density is calculated using Equation (5), for all six cities at various levels, as shown in Figure 12. The density of the road network is the total amount of the road network length to the size of the city multiplied by 100 [87,88]. Here, in this case, the total area in the range of 0-400m, 0-1200m and the entire city is calculated accordingly, the road lengths are considered and summed.
where D is the the density of the roads network, ∑L is the summation of the total road's length, and S is the city area (square kilometers). 15  Similarly, network density is calculated using Equation (5), for all six cities at various levels, as shown in Figure 12. The density of the road network is the total amount of the road network length to the size of the city multiplied by 100 [87,88]. Here, in this case, the total area in the range of 0-400 m, 0-1200 m and the entire city is calculated accordingly, the road lengths are considered and summed.
where D is the the density of the roads network, L is the summation of the total road's length, and S is the city area (square kilometers).
Similarly, network density is calculated using Equation (5), for all six cities at various levels, as shown in Figure 12. The density of the road network is the total amount of the road network length to the size of the city multiplied by 100 [87,88]. Here, in this case, the total area in the range of 0-400m, 0-1200m and the entire city is calculated accordingly, the road lengths are considered and summed. = (∑ / ) * 100% (5) where D is the the density of the roads network, ∑L is the summation of the total road's length, and S is the city area (square kilometers).

Detour Index
The detour index (DI) is a measure of the amount of detour from the shortest route that connects two points, as shown in Equation (6) [89]. The detour index sums up a relation in a single index between network distance and the Euclidean distance. A higher detour index value means that the network distance deviates to a greater extent from the Euclidean distance and vice versa.
(Here, road network distance of all ASFs to PT/Euclidean distance of all ASFs to PT and then multiplying by 100).
As predicted, calculated average network distances are longer than the average Euclidean distances for all the six cities, as shown in Table 6. The detour index values for the six cities are presented in Figure 13. The detour index (DI) is a measure of the amount of detour from the shortest route that connects two points, as shown in Equation (6) [89]. The detour index sums up a relation in a single index between network distance and the Euclidean distance. A higher detour index value means that the network distance deviates to a greater extent from the Euclidean distance and vice versa.

=
Network Distance Euclidean Distance × 100 (Here, road network distance of all ASFs to PT / Euclidean distance of all ASFs to PT and then multiplying by 100). As predicted, calculated average network distances are longer than the average Euclidean distances for all the six cities, as shown in Table 6. The detour index values for the six cities are presented in Figure 13.

Normalization and Benchmarking for Comparing Cities
Different units of indicator values such as distances, density percentages have been standardized as a common unit in a standard form. Such normalization is necessary to prevent inconsistencies, as one indicator value may be between 0 and 200 and other indicator value may be within the range of 0 and 50. In this study, vector normalization is used based on Equation (7). The different dimensions of attributes are transformed into non-dimensional attributes that allow comparisons between criteria [90]. Since different criteria are generally measured in different units, thus it brings into a single unit ranging from 0-1. It creates an evaluation matrix consisting of m alternatives and n criteria, with the intersection of each alternative and criteria given as X ij , therefore, have a matrix (X ij )m × n [91]. Generally, this normalization technique is used in studies of Multiple-Criteria Decision Analysis such as TOPSIS, AHP [92,93] The matrix (X ij )m × n is normalized to form R = (r ij )m × n where i = 1, 2, . . . . . . , m and j = 1, 2, . . . , n.
Benchmarking is used to compare the level of performance of the indicators that are used to measure the performance of cities. This is further comparable between countries and allows verification for others in the targeted level [21,93].

Sensitivity Analysis
Sensitivity analysis is a study of how uncertainty in the output of a mathematical method or system (numeric or otherwise) can be divided and assigned to unique sources of uncertainty in its inputs. Sensitivity analysis tests the robustness of the results, i.e., whether the results produced using a certain method are accurate and identical to those obtained by other methods, not skewed or affected by a certain phase involved in the above normalization process [94]. Here, this analysis is performed only to check whether the ranks produced by both normalization methods are similar. Although in this study the vector method is adopted to rescale the indicators score, the resultant for all ASFs is also compared with results obtained using another rescaling method (minimum-maximum) of normalization. Here the transformation is based on a range of values Equation (8)

Results
The outcomes drawn from the analysis are put into two stages. Initially, in the first stage, the study compares the difference between the Euclidean and network distance for all ASFs from bus stops by evaluating it for all six cities. The first stage includes four assessment criteria-(i) the percentage of ASFs available within the various range of distances 0-400, 401-800 and 801-1200 m, (ii) correlation between Euclidean and network distances for entire city, (iii) mean of Euclidean and network distances for entire city, (iv) difference between Euclidean and network distances within various range of distances for ASFs sector-wise (0-400, 401-800, 801-1200 m). In the second stage, the discussion is carried out by relating the criteria (v) built-up area, network connectivity, network density with the percentage of ASFs in the range of 0-400 m, and the detour index. Finally, (vi) normalization is done by merging the results of the two stages. The comparison between Euclidean and network distance (0-400 m) for very high-level ASFs is shown in Figure 14. The result shows 50% to 69% of combined ASFs are easily accessible through Euclidean distance of 0-400 m range, but through network distance, the combined ASFs percentage for the same range is nearly half of the Euclidean distance in all the six cities. This observation suggests the presence of local natural and humanmade physical barriers; undesirable location of bus stops and bus routes; and network issues, which tend to increase the distance and lower the accessibility to ASFs. (ii) Correlation between Euclidean and network distance The overall Euclidean and network distance correlation is performed for entire ASFs in all the cities using Pearson's correlation method to find the relationship. First, evidence shows high correlation at very similar rates in all the cities except Indore city as shown in Table A3 in Appendix A. It has been verified at a 0.01 (2-tailed) significant level. Since, the R 2 variance only shows the relationship between Euclidean and network distance, hence descriptive statistics analysis is performed to identify the difference between Euclidean and network distances in meters. (iii) Mean of Euclidean and network distance Descriptive statistics are performed in the Social Sciences Statistical Package (SPSS 16.0) for the selected six cities, and the combined mean average of ASFs for every city is given in Table 12. It provides scope for identifying the average difference in the city. The difference is very less in Nagpur city with 212.4m compared to other cities while it was higher in Jaipur city.  The percentage of ASFs across the identified individual sectors (example hospital and clinics) is observed by comparing Table 7 with Table 9 for the (0-400 m) very high accessible range. In Visakhapatnam city, for hospitals and clinics, the low difference (9.77%) between Euclidean and network distance indicates high accessibility by using the network. In Nagpur city, for playgrounds the high difference (26.31%) between Euclidean and network distance indicates low accessibility by using network in 0-400 m range. For other services, the percentage difference is 44.44% between Euclidean and network at a 0-400 m range which is high when compared to other sectors in Bhopal city.
The results from Tables 8 and 10 are compared for the range of 0-400 m combined ASFs. It is observed that Indore city is having the highest percentage by Euclidean distance while Jaipur city is having the lowest. This scenario is entirely different in case of network distance, where almost 30% of combined ASFs are inaccessible in the case of Indore city as shown in Figure 14. It is observed that almost 70% of combined ASFs are accessible within 1200 m range of network in all the cities (Table 10). Above all, hospitals and clinics, other services have high accessibility to bus stops in all the cities in the range of 0-400 m compared to the rest of the ASFs through Euclidean as well as network distance.

(ii) Correlation between Euclidean and network distance
The overall Euclidean and network distance correlation is performed for entire ASFs in all the cities using Pearson's correlation method to find the relationship. First, evidence shows high correlation at very similar rates in all the cities except Indore city as shown in Table A3 in Appendix A. It has been verified at a 0.01 (2-tailed) significant level. Since, the R 2 variance only shows the relationship between Euclidean and network distance, hence descriptive statistics analysis is performed to identify the difference between Euclidean and network distances in meters.

(iii) Mean of Euclidean and network distance
Descriptive statistics are performed in the Social Sciences Statistical Package (SPSS 16.0) for the selected six cities, and the combined mean average of ASFs for every city is given in Table 12. It provides scope for identifying the average difference in the city. The difference is very less in Nagpur city with 212.4m compared to other cities while it was higher in Jaipur city. (iv) Difference between Euclidean and network distances within various range of distances for ASFs sector-wise (0-400, 401-800, 801-1200 m) The correlation, mean difference between Euclidean and network distance for ASFs were observed, as they influence the identification of network connections in the city for reaching bus stops. While considering the entire city, high correlations were found among all the cities as shown in Table 11. Now, by considering it for different threshold ranges 0-400, 401-800, 801-1200 m for ASFs by sector wise in all the six cities, as shown in Table A4 in Appendix A. For the range 0-400 m, the observed correlation is 0.02 to 0.86, p < 0.001, for 401-800 it is 0.09 to 0.79, p < 0.001; and for 801-1200 m it is 0.25 to 0.91, p < 0.001. These values show significant correlations when comparing all sectors of ASFs, and the range varies as high as R 2 = 0.91 and as low as R 2 = 0.02 in all the six cities. Even in this case, very low correlation is observed in Indore city compared to the other five cities. The average difference between the Euclidean and network distances within each range in all six cities is as high as 50% in each sector for all ASFs. The resultant difference between Euclidean and network distances vary which is very logical as it is based on ground reality.
(v) Relating built-up area, network connectivity and network density with a percentage of AFSs in 0-400 m range: The relationship between built-up and network connectivity and density is one of the important characteristics for identifying the city's accessibility [56]. In all six cities, as shown in Figure 15, it is observed that Indore city has the highest built-up area percentage and Visakhapatnam city has the lowest.
Despite less network density and built-up area, Visakhapatnam city has the highest percentage of combined ASFs accessible using network distance from bus stops compared to other five cities. ASFs are accessible using network distance from bus stops. In the case of Jaipur city, even though it has both high network density and better connectivity, as shown in Table 11 but has a low percentage of combined ASFs present in this range which in turn affects accessibility. For Lucknow city, the network density is very high compared to the other five cities.

lowest.
Despite less network density and built-up area, Visakhapatnam city has the highest percentage of combined ASFs accessible using network distance from bus stops compared to other five cities. ASFs are accessible using network distance from bus stops. In the case of Jaipur city, even though it has both high network density and better connectivity, as shown in Table 11 but has a low percentage of combined ASFs present in this range which in turn affects accessibility. For Lucknow city, the network density is very high compared to the other five cities. Visakhapatnam has the lowest detour index for all combined ASFs 0-400m range of accessibility, as well as for the entire city network, shown in Figure 13. The length of edges within the road network connecting bus stops from ASFs is lower compared to the rest of the cities. While it is upturned in the case of the city of Jaipur, where the edge length is 0-400 m high. In the case of Indore, the detour index is very high, which indicates that the edge length is substantially large in the city connecting bus stops from all combined ASFs compared to rest of the city. It is evident from the findings that a few cities are strong in one dimension and weak in another. To generate consistent results from the different analyses, normalization is done. Visakhapatnam has the lowest detour index for all combined ASFs 0-400 m range of accessibility, as well as for the entire city network, shown in Figure 13. The length of edges within the road network connecting bus stops from ASFs is lower compared to the rest of the cities. While it is upturned in the case of the city of Jaipur, where the edge length is 0-400 m high. In the case of Indore, the detour index is very high, which indicates that the edge length is substantially large in the city connecting bus stops from all combined ASFs compared to rest of the city. It is evident from the findings that a few cities are strong in one dimension and weak in another. To generate consistent results from the different analyses, normalization is done.

(vi) Results of Normalization and Sensitivity analysis
The results of Euclidean and network distances are separately evaluated, and ranks were given to the cities accordingly. The performance was evaluated for the high range of accessibility considering 0-400 m and 0-1200 m for both Euclidean and network distance. For Euclidean distance, the percentage of ASFs accessible (Table 7), built-up area ( Figure 8) and the average distance from Appendix A Table A4 were considered and normalized as shown in Tables 13 and 14. Similarly, for network distance, percentage ASFs accessible (Table 10), network connectivity (Table 11) and the average distance from Appendix A Table A4 were considered and normalized and the results are shown in Tables 15 and 16. The results of vector method are compared to the minimum-maximum method for the same indicators as those produced by vector method from Table 13. By rearranging them rank wise, it is observed that scores of both ways vary, but the order of ranks was quite similar, as shown in Table 17. In the vector method, the division of metric units is more precise compared to min-max method. From the overall results of normalization, it is observed that within 0-400 m and 0-1200 m range of Euclidean distance, Indore city ranks first compared to other cities, as shown in Table 13. However, when considered through the network distance, it is ranked lower. When considered through the network distance, Bhopal city ranks first for network connectivity, and the mean distance to reach ASFs from bus stops is high. In the case of Visakhapatnam city for Euclidean distance, it is poorer because the built-up area is very low in the 0-400 m range. However, when compared by network distance, it ranks second. This is because a large amount of combined ASFs are accessible in this range refer (Table 10) compared to cities of Lucknow, Jaipur, and Indore.

Discussions
In this section, the choice of apt indicators to measure accessibility is discussed elaborately. The distance approaches used to evaluate accessibility levels and the ranking of the cities are considered in the following section. Limitations and data constraints for the present study are debated. Overall, how the results can be applied for planning is discussed in detail as follows:

(i) Choice of indicators and methodological concerns
The prime focus is the selection of the most appropriate indicators for defining the accessibility to ASFs from PT. This may depend on several factors, such as the aim of the research or application, the number, and types of data available, the methodological assumptions and the choices (type of distance) made. The evaluation of ASFs is performed as the study aims to understand variation in accessibility in various cities from PT, and particularly which ASFs (such as hospitals and clinics) are more accessible Using simple distance measures such as Euclidean and network may be most appropriate [38], as they are localized and easy to evaluate ASFs from PT. The choice of built-up area for evaluating the distance in this study becomes crucial for identifying geo-located features and labels. It acts as a background to categorize levels of accessibility as per distance using the grid-based method. In several developed nations, the importance of accessibility instruments in planning practice has increased [41]. But with developing nations, this is still in the commencement stage. However, nowadays, availability of data in OSM has facilitated GIS implementation with the transport module (e.g., ESRI Network Analyst Extension) [39]. The possibility of focusing on city level accessibility using various parameters for spatial assessment in this research became easier in comparing various cities. Hence, the assumption made in the hypothesis was proved, as this study consists of various indicators (based on physical elements such as built up, connectivity and density of ASFs). Although the accessibility of each ASFs differs, grouping as per various accessibility level by evaluating combination of influencing indicators, various scores and ranks are assigned.
(ii) The distance approaches Euclidean distances are easy to measure in ArcGIS 10.5.1 software and are used in planning applications traditionally [39]. They do not need a network street layer when measuring the distance between ASFs and PT. However, while measuring by network measures, some additional work is normally needed when setting, such as manually snapping points to the network. Euclidean distances do not adequately represent the actual distances. Accessibility can be measured by different methods as discussed in the literature but measuring accessibility by network distance method is an essential tool for any study and at any scale (city, area and neighborhoods) as it focuses on ground reality [48]. This helps to identify the problems at a deeper level, making it easier for both municipal and local authorities to provide recommendations for improvement in preparing a master or regional plan. The network distances used are more accurate and closer to the actual distances. Since they are calculated using road segments, they represent the actual distances that the user has to walk/drive to reach ASFs, which is very well understood from the literature [40]. However, network distances could minimize this effect because, when generating the network in the ArcGIS 10.5.1, ASFs centroids are moved to the nearest element of the street network before the distance is calculated. Although there may be opportunities to use network distance and accessibility points, researchers may prefer Euclidean distance to reduce computation time, when accessibility measures are considered if there is no availability of accurate network data.

(iii) Limitations and data constraints
Data in research play an enormous role while choosing elements for assessment usability, mentioned in the literature analysis, this study uses secondary data. In such cases, there is no choice but to use what one has. However, the quality and accuracy of the data are determined. In all the cities, indicators were measured within the municipal area. Since no spatial data of ASFs outside the city are considered; for this reason, the 'border effect' may exist and influence the result of Euclidean and network distances which is reported by similar studies of Apparicio [36,37].
(iv) Implications of results in planning The findings of this study encourage that almost 70% of ASFs are present within 1200 m of bus stops through network distance in all the six cities. However, it is also suggested to increase density of the ASF's within 400 m in areas nearer to bus stops. This will lead to a higher percentage of ridership by bus. The method of statistical assessment to perform linear regression is easy to understand. The analysis shows a significant correlation for various indicators. In particular, a strong correlation has been found between indicators measured with Euclidean and network distances at the entire overall city level (higher than 0.55, p < 0.001) in all the cities. Significant spatial variations in geographical accessibility were observed across six cities. The analysis demonstrates that regardless of the usage of travel mode and population density, the majority of ASFs in the cities of Lucknow and Jaipur have low access to PT through network distance. Meanwhile, it is high in cities of Visakhapatnam and Indore which is significantly different from results of Euclidean distance. Moreover, the entire data are represented in percentages, which helps in decision making and can be used in further studies that encourage public transportation, NMT infrastructure and high priority bus lanes. The red colour marked on grids and network maps stands out as a challenge, showing the implementation of upcoming bus stops in the future.
There is an uneven distribution of ASFs, where most are concentrated only at the center of the city like for Visakhapatnam. This hampers the availability of the nearest bus stop to a few ASFs at other locations, particularly in newly developed areas. Consideration of built-up area and network density will help to improve the number of additional ASFs within the high access level that are closer to the bus stops and targets in improving the network condition where development would occur in the future. This is provided by the spatial metrics studied in this research.
The location of bus stops used to assess the reliability of bus services have a direct impact on the performance of the transit system and services. However, this study does not evaluate errors in the choice of distance types with aggregation method [36], cost and time [58], and population [38]. Such limited exactitudes could be combined for future studies to automate demand-based bus stop locations which emphasize land use distribution and demand for the service. The results suggest policy and design implications in all the six cities, such that in the initial stage, new regulations of existing guidelines must be implemented at the policy level. The upcoming infrastructure of ASFs must be planned as per consideration by measuring the distance to any nearest PT. Pedestrian planning guidelines in Indian code IRC: 103-2012 should be revised to include distance factor along with the quality of pedestrian infrastructure [95]. The new network should be planned by providing a more significant number of links connecting the nodes. Comparing cities develops a common framework to identify the level of access to each town and issues in those cities [9]. Application of this kind in similar studies will help in comparing cities with benchmarked cities globally and will surely pave a new way. This study focuses on city level assessment, due to its scale/density of activities, it will further hamper the accessibility of the city in future. The ranks are assigned to understand the position of the city by comparing its performance with other cities and to compete to stand as a benchmark and become a target for other cities to achieve and improve PT. Summarizing the nature of ranking in this study shows city of Bhopal, Indore, Visakhapatnam has high rank which stands out as a benchmark, whereas Nagpur, Jaipur, and Lucknow scores low rank for network distance 0-400 m range. In general, the higher the rank, the better it encourages accessibility for ASFs from bus stops.

Conclusions
The study adopts a unique approach as it integrates several ASFs and evaluates spatial geographic access to ASFs from bus stops by using both Euclidean and network distance measures in six Indian cities. The study tries to identify the level of accessibility to every ASFs sector-wise connecting the bus stops. Some ASFs from a place of residence may have high accessibility at the neighborhood level but still they do not have a decent bus stop connectivity. Sometimes, people going to the nearest ASFs in a neighborhood may be debated, as they may prefer to move to other ASFs further from home. This study considers the links for the walkable distance between bus stops and ASFs to encourage people to use public transport. The grid map overlapped with the built-up area helps in identifying areas with no transit at all and areas with only a few bus stops. The present work proposes and discusses a set of indicators that can be used in urban planning further. Measures are divided into two major categories: Euclidean and network. First accounts using straight line distance while the second calculates distance by using road network. More precise and accurate results are obtained when using network distance, although using network distance requires more data and analysis (e.g., a network layer). The results obtained by using Euclidean distance overestimate the distance. The results of this research indicate the necessity of land use planning by increasing the density of ASFs in areas close to bus stops, concentrated within 400 m of transit, which in turn leads to a higher percentage of ridership by bus.
In India, current policies on accessibility target on building up new infrastructure, but they do not define how accessibility is achieved in cities and assessment measures required for accessibility planning. In this research, various ASFs are identified based on levels of accessibility and combined with other indicators for ranking. Overall, the percentages at various levels and scores can be considered for further studies and will help researchers and planners in terms of the methodologies used, accepted assumptions, and limitations.
In any context, indicators used for the assessment of the accessibility of ASFs from PT represent very important tool for planning of urban areas in India. They may provide a framework for developing different public transit policies in municipalities and in metropolitan areas (regions, metropolitan areas).
Further research might include the refinement and addition of indicators so that they can further strengthen the integration of various land uses connecting to PT with relevant factors that affect accessibility. A similar methodology can also be adopted in other cities in India for evaluating accessibility. This research explores the possibility of including some of the factors such as slope, elevation and road condition that affects accessibility to reach PT irrespective of shortest distance. These methods can be further applied to other transit such as metros in future, as the choice of origins and destinations may influence interpretation differently for a city.  Acknowledgments: This research is funded by the scholarship awarded to Yenisetty Pavan Teja from Ministry of Human Resource Development (MHRD), Government of India. We also like to thank the editor and two anonymous reviewers for their constructive comments, which helped us to improve the manuscript.

Conflicts of Interest:
The authors declare no conflict of interest.  Figure A3: Overlapping ASFs on Google Earth; for accuracy assessment in 3 areas by ground truthing Figure A3. Overlapping ASFs on Google Earth; for accuracy assessment in 3 areas by ground truthing. Table A2. Accuracy assessment and Kappa coefficient for Trimurti Nagar area (ASFs).

Trimurti Nagar Shankar Nagar Dhantoli
Overall Accuracy Kappa coefficient Figure A4. Overall accuracy assessment and Kappa coefficient results for 3 areas.