Urban Planning Using a Geospatial Approach: A Case Study of Libya

Large scale developmental projects firstly require the selection of one or more cities to be developed. In Libya, the selection process is done by selected organiza-tions, which is highly influenced by human judgement that can be inconsiderate of socioeconomic and environmental factors. In this study, we propose an automated selection process, which takes into consideration only the important factors for city (cities) selection. Specifically, a geospatial decision-making tool, free of human bias, is proposed based on the fuzzy overlay (FO) and technique for order performance by similarity to ideal solution (TOPSIS) techniques for development projects in Libya. In this work, a dataset of 17 evaluation criteria (GIS factors) across five urban conditioning factors were prepared. The dataset served as input to the FO model to calculate weights (importance) for each criterion. A support vector machine (SVM) classifier was then trained to refine weights from the FO model. TOPSIS was then applied on the refined results to rank the cities for development. Experimental results indicate promising overall accuracy and kappa statistics. Our findings also show that highest and lowest success rates are 0.94 and 0.79, respectively, while highest and lowest prediction rates are 0.884 and 0.673, respectively. The results suggest that demography and vegetation are the two most influential factors contributing to the selection of city for development in Libya. This study is limited to analysis of six cities; the procedure developed through this study can be extended to other cities. It is of the opinion that evaluated criteria can be adjusted according to the environment and the current development of the cities.


Introduction
Rapid population growth and urbanization have caused many problems in the implementation of developmental projects in cities. Haphazard infrastructural project execution that includes disregard in prioritizing city (or cities) selection is also a factor hampering sustainable development practices. Development projects that rely on selected organizations, which in turn rely on human judgment, can lead to unrealistic criteria evaluation, causing delays in project execution [1]. However, the fact remains that continued infrastructure development is unavoidable, especially since urban cities constantly need to evolve and grow to keep up with the times [2].
The selection of a city (or a group of cities) is one of the most important steps for sustainable development. The selection criteria must ensure that the city (cities) has high priority for development and is (are) in line with the needs of the local citizens.
Moreover, timely selection requires effective planning and analysis and must consider multiple conflicting and disproportionate factors (such as those that have critical socioeconomic and environmental implications to different stakeholders). Urban planning application using remote sensing (RS) and geographical information systems (GIS) is one of the many areas that can be explored for city selection. Such applications would not only eliminate human bias but would also be able to make more objective decisions based on data.

Remote sensing and GIS applications in urban planning
Remote sensing can be applied in different aspects of urban planning such as (but not limited to) urban traffic analysis, urban environment analysis (air and water pollutions), and urban expansion. With recent developments in remote sensing technologies, remote sensing data can be exploited for urban studies. One example is the classification of land use based on high spatial and spectral resolution data such as orthomosaic and elevation images. Multidimensional spatiotemporal data can now be reliably obtained by sensors in different scale ranges and with flexible repetition rates [3].

Remote sensing sensors for urban planning
Medium-to high-resolution satellite imagery can be used by urban planners and land managers to monitor land conditions to support decision-making for sustainable urban development. Remote sensors are able to provide voluminous amounts of data, which can be exploited to produce/update GIS maps or for detection changes in urban land covers. High-resolution satellite sensors available on IKONOS, for example, can collect diverse geospatial data for studying vegetation. The sensors can sense 4 m resolution multispectral and 1 m resolution panchromatic, Quickbird imageries with 2.4 m resolution multispectral and 61 cm resolution panchromatic, and Worldview-4 imageries with 1.24 m resolution multispectral and 31 cm panchromatic. Medium-resolution satellite sensors, available on Landsat-8, Sentinel-2, and SPOT, are also valuable data sources for urban and vegetation change detection from various time periods during the same season, which further supports analyzing any past changes. Analysis of such data can then be used for decision-making and planning for further development of a particular urban area [4].

Integration of GIS and remote sensing for urban planning
Remote sensing data can be integrated with other spatial data to perform various types of full-fledged assessments. GIS techniques can be utilized to integrate the required spatial data and analytic data from various sources, such as field survey data, topographic maps, aerial photographs, and also archived data. The data can be represented as location (i.e., latitude and longitudes) or even as tabular attributes. GIS techniques play a substantial role in the data integration process of multilayer spatial information along with statistical information in various developmental scenarios [5].

Methods and approaches for urban planning
Multi-criteria decision-making (MCDM) is concerned with making a decision by evaluating multiple conflicting criteria. It embodies various methods and procedures where the gist is the formal incorporation of multiple conflicting criteria in the analytical process [6]. In the context of GIS, this refers to the spatial decisionmaking process based on GIS data with geolocation tags. Spatial decision-making techniques have been used to solve many GIS problems such as locating solar plants, urban planning, and project construction optimization [7]. Advanced MCDM methods include simple additive weighting (SAW) [8], analytic hierarchy process (AHP) [9], and TOPSIS [10]. Fuzzy set theory and random set theory are also MCDM techniques that incorporate sophisticated algorithms to resolve uncertainty in data [11][12][13][14].
TOPSIS is a MCDM technique that deals with real-world problems. It basically ranks criteria on the basis of the shortest distance from the positive ideal solution (PIS) and the farthest distance from the negative ideal solution (NIS) [15]. The work in [2] illustrates the application of a GIS-based MCDM tool for urban infrastructural planning. Awasthi et al. [16] presented a fuzzy TOPSIS method for selecting the best location for an urban distribution center in Canada. Uysal and Tosun [17] proposed a fuzzy TOPSIS-based maintenance management system using 17 criteria categorized under 5 contending parameters. The criteria were deduced from questionnaire feedbacks and interviews administered to company maintenance managers. In addition, Momeni et al. [18] presented a fuzzy TOPSIS-based method for maintenance strategy selection. Baysal et al. [19] developed a two-stage fuzzy method to determine the best sub-municipal projects among a set of proposed projects. The method simplifies the selection process and provides an objective decision outcome for stakeholders. Shelton and Medina [20] presented an integrated method to prioritize transportation projects in Wilmington Area, USA, based on multi-criteria decision support systems, AHP and TOPSIS methods. The process optimally selects the important routes that best serve the interest of the general public.
Based on the literature, TOPSIS has been successfully applied in many fields, producing reasonably accurate results. This study proposes an automated TOPSISbased solution for prioritizing urban projects based on criteria that meet sustainable development. Specifically, this work addresses the following questions on the value of remote sensing (and GIS) to urban planning: 1. Which remotely sensed dataset(s) is (are) useful for urban planning? 2. Which criteria can be derived from remotely sensed data?
3. What are the major factors that need to be considered in urban developments?
From these questions, this study further looks at the automated prioritization of urban projects based on criteria that meet sustainable development practices. The specific objectives are (i) to identify factors that play major roles in urban development and (ii) to develop a geospatial solution based on TOPSIS for prioritizing projects for urban development.

Study area and preparation of the conditioning factors
This study focuses on Libya, a country in the Maghreb region of North Africa ( Figure 1). Libya borders the Mediterranean Sea to the north and Egypt to the east. Along the southeast of Libya is Sudan, Chad. To the south is Niger. Algeria and Tunisia constitute the western border. Libya is the 17th largest nation in the world with a landmass of over 1,759,540 km 2 . The study area in the northern part of Libya covers six districts, namely, Darnah, Al Jabal Al Akhdar, Benghazi, Al Marj, Al Qubbah, and Al Hizam Al Akhdar ( Figure 1). Libya is geographically bounded between 20°00 0 00″ E and 23°30 0 00″ E and 31°00 0 00″ N and 33°00 0 00″ N. The climate in Libya is categorized by hot and dry summers with high temperatures. The mean annual temperature in the coastal region ranges from 14.2°C (Shahat) to 21.0°C (Tripoli Airport) and at stations in the interior region (inland) between 21.3°C (Al Qaryat) and (Ghat) 23.4°C (1945-2009). Libya is one of the driest countries in the world with mean annual rainfall along the Libyan coast ranging between 140 and 550 mm and rarely exceeding 50 mm in the interior regions . December and January are the wettest months with 6 months (October-March) receiving 87.1% of the total annual precipitation. The majority of rainfall occurs in the winter season with the rainy season beginning in September-October and ends in March-April [21].

Data preparation
The data used in this study include Landsat satellite imagery acquired in the year 2017 with 15 m resolution panchromatic and 30 m resolution multispectral, Shuttle Radar Topography Mission (SRTM) digital elevation model (DEM) downloaded from USGS data archive with 30 m resolution, population density map obtained from GHSL with 250 m resolution, road network map from Diva-GIS, and MODIS satellite imagery from where the land surface temperature with 0.25°resolution was derived. Other data include rainfall data at 0.25°, net primary productivity (NPP) at 0.1°, NDVI at 0.1°, and air quality (CO, NO 2 ) at 0.25°resolution ( Table 1). Details of the Landsat data are presented in Table 2. Seven set of images with overlapping areas were acquired between 4 February 2017 and 1 March 2017. In addition, the highest cloud cover was 1.03%, which does not pose a problem for land use information extraction from the study area. Since the images have overlapping areas, they were preprocessed and mosaicked to create one seamless image of the area for effective and efficient processing ( Figure 1d).

Data preprocessing
Four preprocessing steps were performed on the Landsat satellite images: (i) Pan-sharpening using a fusion of the panchromatic and multispectral bands for the enhancement of the spatial resolution of multispectral band; (ii) atmospheric correction, which is applied to correct the atmospheric distortion by retrieving surface reflectance and engage topographic correction as well as adjacency effect correction; (iii) radiometric correction, which converts radiance values to the pure surface reflectance to enhance image capability and contrast; and (iv) mosaicking to create one seamless image coverage of the area for effective and efficient processing [22]. The MODIS source data was preprocessed using MODIS Conversion Tool Kit (MCTK). Note that the spatial resolution of the MODIS dataset varied according to the source. However, during the preprocessing, they were resampled to 30 m to match the DEM and Landsat resolutions.

Urban conditioning factor dataset
In this study, we considered 17 critical urban conditioning factors for selecting the most suitable city or cities for sustainable urban development. The factors are grouped into five main categories: (i) topography, (ii) land use and infrastructure, (iii) demography and climate, (iv) vegetation, and (v) air quality.

Topography
Topography is a very important consideration for urban development projects [23]. For this study, altitude and slope are the two main factors related to topography. Altitude is important for citing facility because it affects the living conditions as well as breathing behavior. The collected DEM shows that the study area is between À4 and 865 m above mean sea level (Figure 2a). The slope factor, which ranges from 0 to 14°(can be classified as almost flat), was also generated ( Figure 2b). Such data is important when estimating cost. For example, any increase in slope will increase the cost of facility installation and maintenance since moving workers, transport vehicles, and machineries will be more difficult (i.e., up and down a slope). Low slope areas also may incur undesirable cost, in the instance of weather anomalies such as dust/sand storms.

Land use and infrastructure
Land use and infrastructure are also important considerations for urban projects. Land use information can show human activity patterns, whereas infrastructure can indicate development status in a particular city. The land use of the study area was derived from Landsat images and refined based on Google maps ( Figure 2c). In this study we applied the SVM classifier, which was based on object-based image analysis (OBIA) using the ENVI 5.3 software. Training sites for the SVM were collected form all land use classes by stratified random method (i.e., at least 80 sample points for each class) [24]. The area contains five main land use types: (i) irrigated crops, (ii) vegetation, (iii) artificial areas, (iv) bare lands, and (v) waterbodies. Most parts of the study area were bare land (desert), which were predominantly located in the middle and southern parts of the study area. The northern part mostly comprised of irrigated crops and artificial areas, specifically man-made features and urban areas.
Road networks play an important role in the country's economy, serving the people by linking main cities to industrial and commercial sites. In this study, three types of roads, namely, main routes, secondary routes, and trail routes (Figure 2d-f), were considered as factors for city selection. The northern part of the city is supported by main routes (Figure 2d). These roads mainly link other cities to Benghazi, which support Benghazi city itself. Main routes span a significant number of kilometers within the study area. The study area also contains several kilometers of secondary routes (Figure 2e). Unlike main routes, secondary routes are found in most parts of the study area and in different cities including Benghazi. Secondary routes mainly support transportation of goods and are used for civil construction projects. Finally, trail routes support the rural areas, mainly for transportation of agricultural produce to the markets. The class of roads plays a vital role in selecting a city for development according to the available budget. Cities that can support more people will normally be prioritized for development projects.
Another important factor related to infrastructure is the percentage of built-up areas in a particular city. This is important because cities with many built-up areas indicate little or no space for new projects. On the contrary, cities with fewer builtup areas mean that they are more suitable for new developmental projects. The Normalized Difference Built-up Index (NDBI), which is a quantity of the intensity of urban area from satellite images, was used in this study [25]. The NDBI was initially established regarding the ratio of bands 4 and 5 of TM sensor. However, the NDBI can be adopted on Landsat-8 data or even any multispectral sensor data [26]. It basically extracts the urban areas where there is an upper reflectance in the short-wave infrared band associated to the near-infrared band. The accuracy of built-up areas extracted using NDBI is reported to be around 93% [27,28]. We calculate the NDBI (Eq. (1)) based on the work in [27]: SWIR is the short-wave infrared band ranging from 1.57 to 1.65 μm. NIR is the near-infrared band in the range of 0.85-0.88 μm. Figure 2h shows the percentage of built-up areas calculated for each district considered in this work. The largest built-up area is found in Benghazi city (>45%). Some cities though, such as Al Marj and Darnah, have very small of built-up areas that are not detectable through satellite data.
Besides infrastructure factors, distance to the city is also a critical factor. Preferably, a city's location should be as near as possible to the capital or large cities such as Benghazi. This is because it facilitates ease of access to better business opportunities and education. Therefore, in this work, the distance to Benghazi city is one of the important parameters. Specifically, the desired distance to the city should range from 0 to 600 km (Figure 2g) so that cities such as Darnah and Al Qubbah (which are as far as 400 km away from Benghazi) are also covered.

Demography and climate
An increase in a city's population often leads to an increase in urbanization. Moreover, if the population increase is rapid, urbanization often happens randomly. This can be a major problem for most cities in a developing country. However, proper planning and effective decision-making can mitigate this problem. In this study, the population density (Figure 2i) was analyzed. The analysis results indicate that the northern part (mostly around Benghazi city) is most populated with a density of 669 people per 250 m cell of raster data.
Climate is a factor that also influences the selection of cities for developmental projects. Land surface temperature (LST) and rainfall are two factors considered in this work (Figures 2j and 3k). The LST map shows that the southern part of the study area (mostly desert with no vegetation) has higher surface temperature compared to the northern parts. Another observation from the map reveals that Benghazi has slightly higher temperature than other urbanized areas. In arid regions, people often prefer to settle in areas with low temperature. The average day-night temperature ranges from 11°to 20°Centigrade. High temperatures are also observed in the west-southern part, whereas the lowest temperature is found in the northern part of Al Jabal Al Akhdar cities.
Rainfall, which is another climate factor, is also considered in deciding the location of settlement and development. This is because rainfall frequency and intensity affect the dryness of the cities, the local climate system, as well as agriculture activities. Figure 2k presents the rainfall intensity of the study area for year 2016 where minimum and maximum rainfall intensities were 48 and 1119 mm per month, respectively. The central part of the area has less amount of rainfall compared to other areas.

Vegetation
Normalized difference vegetation index (NDVI) is an indicator derived from remote sensing satellite data. It is mostly used to monitor vegetation cover over any area on the planet. It serves as a good indicator for vegetation cover of the study area. The presence of abundant vegetation is able to lower the local temperature as well as reduces the negative effects of noise and air pollutants. In the study area, the northern part has higher NDVI compared to the south (Figure 2l).
Larger amounts of vegetation can indicate higher vegetation productivity. Having higher vegetation productivity helps asses the net primary productivity (Figure 2m). Plant productivity plays a major role in the global carbon cycle by absorbing some of the carbon dioxide released through coal, oil, and other fossilfuel burning. Large NPP values are found in the southern part of the study area.

Air quality
Air quality directly affects the environment and consequently people's health. In this work, we have considered the CO and NO 2 (Figure 2n and o) air quality indicators. In 2016, higher CO and NO 2 levels were measured in the southern part of Benghazi. Benghazi city also recorded high levels of these gasses for the year under investigation. The air quality data was extracted from the MODIS source with a resolution of 0.25°. We utilized the ENVI 5.3 software to process the MODIS imagery. However, in order to prepare unprocessed MODIS satellite images for analysis, they must firstly be converted into ENVI format. This was done using the MODIS Conversion Tool Kit.

Modeling process
This section describes the modeling process, specifically the application of TOPSIS for scheduling and prioritizing the cities for urban development ( Figure 3). First, a medium-resolution Landsat-8 satellite image from the study area was acquired and preprocessed. Then, the image was segmented using a multiresolution segmentation algorithm and classified into several classes using object-based image classification. The multiresolution algorithm has three main parameters, namely, scale, shape, and compactness. Since these parameters are data and application dependent, in this study, we had to select them empirically via trial and error. This meant that the best values were determined via visual examination of the segmentation results. After the segmentation process, several attributes were selected and used as class predictors in the classification algorithm. From the spectral attributes, the five bands of the Landsat-8 image were selected. For spatial attributes, shape index, roundness, compactness, and density were used [29,30]. In the classification step, the support vector machine (SVM) algorithm was used. Although the SVM is a relatively simple binary classifier, it has very good generalization capabilities if properly trained [31,32].
Several other digital data such as DEM and population density were also obtained from various online sources. The factors used as described in the previous section are widely reported in the literature for selecting urban projects or relevant projects. Fuzzy overlay (FO) analysis was carried out to determine the importance of each parameter to achieving the goal of the study. The SVM classifier was further applied to refine the results obtained from the FO model. Finally, the cities were sorted according to their importance by applying the TOPSIS model on the results of the SVM.

Fuzzy overlay and TOPSIS models
Fuzzy overlay analysis is based on the fuzzy set theory that relies on membership relationship of events to define specific sets or classes [33]. Operationally, FO is similar to overlay analysis but differs in the reclassified values and results from the combination of multiple criteria. It involves problem definition, partitioning into sub-models and determining the significant layers. FO transforms the data to a common scale and defines the likelihood of the data belonging to a specific class, for example, slope values being transformed into the probability of fitting into the favorable suitability set based on a scale of 0 to 1, expressed in terms of membership [34]. Input raster are not weighted in FO since the transformed values indicate the possibility of membership rather than using ratio scale as with weighted overlay and weighted sum. The equation using fuzzy Gaussian function can be given as [35] μ The inputs f i and f 2 are the spread and the midpoint, respectively. Midpoint can be a user-defined value with a fuzzy membership of 1. The default is the midpoint of the range of values of the input raster. Spread defines the membership of the Gaussian function. It generally ranges from 0.01 to 1. Increasing the spread causes the fuzzy membership curve to become steeper. Fuzzy overlay analysis quantifies the possibilities of each cell or location to a specified set based on membership value.
As previously mentioned, the results of FO are refined using the SVM, which develops a linear regression between suitability status and criteria factors. SVM aims to determine an optimal separating hyperplane (maximizing the margin width) between two classes in feature space [36]. The training points near the hyperplane are called support vectors and are utilized for classification once the decision line/surface is obtained. The separating hyperplane is found as follows: where w is the cofficienct vector that defines the hyperplane orientation in the feature space, b is the offset of the hyperplane from the origin, and ε i is the positive slack variables. The optimal hyperplane is found by solving the following optimization problem [36,37]: where α i is the Lagrange multiplier and C is the penalty for data classification. The following decision function is applied as follows: Developed in [38][39][40], TOPSIS is a multi-criteria decision tool based on the intuition that a selected alternative has the shortest possible geometric distance from the PIS. In other words, the alternative has the longest geometric distance from the NIS [41]. The analysis compares a set of alternatives by assigning weightage to each criterion to compute the geometric distance between possible alternatives to determine the ideal alternative based on the assumption that the criteria uniformly increases or decreases. TOPSIS allows trade-offs between criteria; a poor result in one criterion can be compensated by a good result in another criterion. TOPSIS provides a more realistic model than non-compensatory methods by including or excluding alternative solutions using hard cutoffs. Consider X ij as the inputs for matrix of priorities where there are i ¼ 1, …, m alternatives and j ¼ 1, …, n criteria. There are six steps associated with the implementation of TOPSIS as follows [42]: Step 1: Construct the normalized decision matrix calculated using Eq. (6): Step 2: Construct the weighted normalized decision matrix using Eq. (7): Step 3: The positive and negative ideal solutions are determined by Step 4: Calculation of separation (positive and negative) measurement using Euclidean distance. Eq. (9) is used to calculate the distance.
Step 5: Closeness to the ideal solution is calculated using Eq. (10): Step 6: Ranking alternatives based on closeness to the ideal solution. TOPSIS has been used in different circumstances (e.g., individual and grouping). By applying the TOPSIS model using the results of the FO as input, the cities were sorted according to their importance for proposed urban development projects.

Results and discussion
According to [43], the most contributing factors to urban suitability are topography, land use and infrastructure, vegetation, demography and climate, and air quality. Therefore, these factors should be thoroughly analyzed to discover the most (and the least) suitable area for urbanization. Hence, in this work, 17 detailed factors were analyzed in order to rank each's importance (via weight assignment) for the selection of a city (or cities) for sustainable urban development. Subsequently, a suitability map was generated based on the FO (Figure 4a) method. A continuous scale was used for suitability weightage, which ranges from 0 (less suitable) to 1 (highly suitable). From the generated map, areas indicated as most suitable are located in the northern parts, especially the areas surrounding Benghazi and the northern parts of Al Marj and Al Jabal Al Akhdar. Sole reliance on the generated map, however, does not help much in deciding city development prioritization. As a result, the map was further refined to make it much more distinct for decision-makers. To do this, the map was firstly reclassified into three categorical classes: (i) not suitable, (ii) less suitable, and (iii) highly suitable. This was done using the natural break classification method (Figure 4b) where several samples were selected from the not suitable and highly suitable areas (results of FO) to generate training and testing data. These datasets were then used to train a SVM to classify between the two classes. Table 3 presents the estimated factors and their coefficients. The result indicates that land use, distance to primary route,   Table 3.
List of criteria, estimated coefficient, and their code.
distance to capital city, rainfall, NPP, and NO 2 have negative effects on the suitability level of the selection process. The remaining factors have positive effects. Among the positive factors, population density has the highest effects on the selection process. Based on the estimated coefficients, the suitability map in Figure 5 was produced. It can be seen that the map reflects the same thing as in the previous suitability map. However, it is clearly more informative for decision-makers. Based on this, the cities were ranked according to their importance using TOPSIS method. Table 4 presents the positive ideal, negative ideal, closeness coefficient, and TOPSIS rank for each of the cities being analyzed. According to the closeness coefficients, the ranking order for the cities is as follows: The ranking results were then used to generate a map for final decision-making ( Figure 6). Cities with green and light green colors are suggested to be prioritized first for development. More details about the TOPSIS analysis can be found in the Appendices.

TOPSIS analysis
Based on the results, the importance of each group of factors was evaluated ( Figure 7). The bar chart shows the importance of the standardized factor weights    Table 4.
Ranking of cities on the basis of importance for urban development. in each group. Demography and vegetation are the two most influential factors with positive contribution, followed by vegetation and topography. The other factors have negative contribution.

Accuracy assessment
Recall that the refined suitability map was recategorized into the three classes of "not suitable," "less suitable," and "highly suitable," These classes reflect the degree of urban development suitability in the study area. The categorized suitability map can be validated accurately through each class. The continuous refined is suitability map ranging from 0 to 1. The most suitable areas that range from 0.751 to 1 fall into high suitable class, while the moderate suitable areas for urban development were extracted from 0.401 to 0.751 from the continuous refined suitability map.

Accuracy metric Value
Correctly classified instances 1178 (78.5%) Incorrectly classified instances 322 (21.4%) Kappa statistic 0.67  Table 6. Accuracy assessment of SVM modeling based on ROC. Consequently, the least feasible areas were classified out of 0.001 to 0.40. The categorized suitability map was then validated based on the same randomly selected samples ( Table 5). The SVM model accurately classified 1178 samples, which is about 78.5% of the total samples tested, which produced kappa index of o.67. The kappa index is calculated using Eq. (11) [44]: where P c is the proportion of number of pixels that are correctly classified and is calculated as

Conclusion
An automated geospatial solution for selecting and ranking cities in Libya for urban development is proposed in this chapter. The suitability map showed that most areas indicated to be suitable are in the northern part of Libya. The results indicate that land use, distance to primary route, distance to large city, rainfall, NPP, and NO 2 have negative effects on the level of suitability for the selection process, whereas the other factors have positive effects with population density taking the lead. It is revealed that SVM model accurately classifies 1178 samples, about 78.5% of the total samples tested which produced kappa statistic of 0.67. The high-priority city was selected as Benghazi that is followed by Al Jabal Al Akhdar. The results suggest that demography and vegetation are the two most influential factors contributing to the selection of city for development in Libya. This study is limited to analysis of six cities; the procedure developed through this study can be extended to other cities. It is of the opinion that evaluated criteria can be adjusted according to the environment and the current development of the cities.