Next Article in Journal
Effects of Canopy and Multi-Epoch Observations on Single-Point Positioning Errors of a GNSS in Coniferous and Broadleaved Forests
Next Article in Special Issue
Impervious Surfaces Mapping at City Scale by Fusion of Radar and Optical Data through a Random Forest Classifier
Previous Article in Journal
Sentinel-1 and 2 Time-Series for Vegetation Mapping Using Random Forest Classification: A Case Study of Northern Croatia
Previous Article in Special Issue
An Improved Index for Urban Population Distribution Mapping Based on Nighttime Lights (DMSP-OLS) Data: An Experiment in Riyadh Province, Saudi Arabia
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Remote Sensing-Based Analysis of Urban Landscape Change in the City of Bucharest, Romania

1
Faculty of Geography, University of Bucharest, 010041 Bucharest, Romania
2
GISBOX, 060586 Bucharest, Romania
*
Author to whom correspondence should be addressed.
Remote Sens. 2021, 13(12), 2323; https://doi.org/10.3390/rs13122323
Submission received: 7 May 2021 / Revised: 1 June 2021 / Accepted: 10 June 2021 / Published: 13 June 2021
(This article belongs to the Special Issue Data Fusion for Urban Applications)

Abstract

:
The paper investigates the urban landscape changes for the last 50 years in Bucharest, the capital city of Romania. Bucharest shows a complex structural transformation driven by the socialist urban policy, followed by an intensive real-estate market development. Our analysis is based on a diachronic set of high-resolution satellite imagery: declassified CORONA KH-4B from 1968, SPOT-1 from 1989, and multisensor stacked layers from Sentinel-1 SAR together with Sentinel-2MSI from 2018. Three different datasets of land cover/use are extracted for the reference years. Each dataset reveals its own urban structure pattern. The first one illustrates a radiography of the city in the second part of the 20th century, where rural patterns meet the modern ones, while the second one reveals the frame of a city in a full process of transformation with multiple constructions sites, based on the socialist model. The third one presents an image of a cosmopolitan city during an expansion process, with a high degree of landscape heterogeneity. All the datasets are included in a built-up change analysis in order to map and assess the spatial transformations of the city pattern over 5 decades. In order to quantify and map the changes, the Built-up Change Index (BCI) is introduced. The results highlight a particular situation linked to the policy development visions for each decade, with major changes of about 50% for different built-up classes. The GIS analysis illustrates two major landscape transformations: from the old semirural structures with houses surrounded by gardens from 1968, to a compact pattern with large districts of blocks of flats in 1989, and a contemporary city defined by an uncontrolled urban sprawl process in 2018.

1. Introduction

Urban landscape change analysis is a necessary approach for the Eastern European cities [1,2], being related to quality of life [3] and socio-ecologic configuration [4,5]. The landscape change analysis plays a key role in understanding the present-day pattern of the city, useful for decision-making, urban planning [6,7], and urban modelling [8,9]. Land cover and land use features are the first issue in urban remote sensing analyses [10] and in the study of city evolution and configuration, while the thematic classification of higher spatial resolution satellite data for the urban area is essential in mapping process for practical uses [11].
Bucharest is the largest urban area of Romania and the sixth biggest European Union capital city with more than 2.1 million inhabitants [12]. The transformation of Bucharest urban structures was massive, as the socialist urban policies from the second half of the 20th century were followed by the real-estate market-driven policies during the last 20–25 years [13,14,15]. The city of Bucharest represents an example of intensive and complex urban landscape change [16]. It shows the effect of a remarkable demographic and industrial growth between 1950–1989, followed by an uncertain transition period between 1990 and 2001 and finally by a real-estate and private investment boom after [17,18,19,20]. During the socialist times, one of the urban policies was implemented around the “working class” people, which materialized with the construction of blocks of flats and factories. People were concentrated in districts of blocks of flats [21]. This process resulted in a new urban pattern with a compact city featured by small parks, narrow streets, and a lack of open spaces and parking places [15]. Moreover, a particularity for Bucharest was the clear delineation of the built-up area by the ring road beyond which dominated the agricultural land [22].
Similar to other former socialist cities, Bucharest followed the same development trajectory in its transition from socialism to democracy [23]. The collapse of communism and the transition to new real-estate market rules was abrupt and led to the development of new and uncontrolled process like urban sprawl or peri-urbanization [24]. Urban sprawl is defined as the unplanned extension of the city by a mix of low density built-up spaces into the surrounding agricultural area [25]. In particular, Eastern European countries are susceptible to urban sprawl. The main consequences of this process are the inefficient use of land resources [26] and the decrease of life quality [27].
Thus, the large temporal frame of about 5 decades was marked by major changes of the political system, followed by a total transformation of the urban planning policies in the context of the transition from the centralized economy to the present-day market driven economy. This process can be modelled with the help of remote sensing data.
The declassified data archive hosted by the USGS Earth Explorer portal [28] is used in many contributions focused on land cover changes and urban structures analysis. Previous contributions mostly combined image interpretation of CORONA KH-4B images with thematic classification of the multispectral Landsat imagery [29,30,31,32,33,34]. Some authors integrated KH-4B data with 2–3 m spatial resolution with Very High Resolution (VHR) imagery in urban area change detection approaches, applied in different regional studies from Kazakhstan [35], Hungary [36], China [37], or Iraq [38].
The latest developments of the Earth Observation (EO) tools in Europe, after the initialization of the European Space Agency (ESA) and Copernicus Programme [39,40,41], opened up new opportunities for urban landscape change research. The launch of Sentinel-1 Synthetic Aperture Radar (S1 SAR) and Sentinel-2 Multispectral Images (S2 MSI) satellites with enhanced resolutions offer new possibilities for LCLU (land cover land use) change detection analysis, using the synergy of the imagery data in order to explore in detail selected areas [42,43,44].
Different authors focused on medium to high resolution satellite data for urban land cover/land use change mapping, but the main limit was the archive data of multispectral images starting from the early 1970s, although intensive transformation of urban structure started earlier, during the 1960s. In this respect, mainly large urban regions or big cities were modelled in terms of LCLU dynamics.
Topaloğlu et al. [45] developed a case study focused on LCLU changes in Istanbul, Turkey based on the integration of Sentinel-2 MSI with Landsat 8 OLI imagery. The thematic classification process applied on Sentinel-2 imagery was assessed using pixel-based approaches [46] or OBIA techniques [47]. Ghanea et al. [48] and Lefebvre et al. [49] focused on urban area LCLU classification.
Due to the improved spatial resolution, SPOT images were used worldwide for urban mapping [50] urban sprawl monitoring [51], change detection analysis [52], and classification of urban land use [53].
Sentinel-1 delivers adequate information for discriminating the features on the ground thanks to the C-band Synthetic Aperture Radar (C-SAR), useful for land cover classification [54,55,56] and mainly for urban footprint delineation [57,58]. There are some studies using the multisensor approach by developing a synergy of S2 MSI and S1 SAR images for urban feature extraction [59,60,61,62] or crop classification [63,64,65,66].
The city of Bucharest was the subject of some studies regarding remote sensing change detection. Most of them are based on Landsat data for land abandonment mapping in the urban sprawl area [67], for the identification of the change pattern of built-up area [68] and for land cover change detection [69]. Sandric et al. [70] compared CORONA KH 7 and Ikonos imagery for mapping and describing urban structural changes in some critical areas. Noaje and Sion [71] integrated CORONA and VHR imagery to retrieve a qualitative analysis of urban landscape change for a sample area in the city. Most of these studies are based on image classification followed by change detection analysis and only indicate the area of change for a few land cover classes. They lack the quantification of classes involved [72].
The aim of this study is to understand and reveal urban dynamics in the capital-city of Bucharest, like super-compaction and artificiality with increase of built-up category. Moreover, qualitative data about urban sprawl and suburbanization can be extracted, as the information is stored in archived earth observation images.
The approach proposes a case study of Bucharest urban landscape change over the last 5 decades, with an emphasis on built-up land cover classes. The basement of the starting point of the study consists of three datasets from different reference years: CORONA KH-4B declassified imagery from 1968, SPOT-1 multispectral and panchromatic imagery from 1989, and a stacked dataset containing Sentinel-1A SAR together with Sentinel-2A MSI images from 2018. The study focuses on the key LCLU classes for the urban landscape change analysis, extracted from diachronic imagery by using complementary approaches: (1) a digital photogrammetric processing of declassified data CORONA KH-4B, followed by (2) visual interpretation and features digitization for 1968 and 1989 years and (3) a digital processing and automatic thematic classification of Sentinel-1 SAR and Sentinel-2 MSI multisensor stacked data for multiple land cover features mapping. Finally, the resulted datasets were integrated into a metric and statistical analysis. Three detailed case studies illustrate the major changes that occurred in Bucharest.

2. Materials and Methods

The city of Bucharest is situated in the south-eastern part of the country, in the Romanian Plain, at 55–97 m above sea level (a.s.l.) (Figure 1), developed along Dâmbovița River floodplain and terrace fields.
The analysis is performed within the official limits of the administrative unit, divided into six sectors, covering 228 km2. Bucharest has a population density of 8,200 inhabitants per km2 and a building density of about 500 units / km2 [12]. The urban configuration is radial and concentric at first sight [73], but it is more complex as an effect of the multi-coredevelopment within each sector. Along with the population increase during the communist time, the city faced the necessity of constructing new living spaces.
The current study proposes a methodology for the integration of historical imagery dataset with recent EO images. Therefore, this approach is based on three datasets taken from different spaceborne sensors presented in Table 1.
The workflow follows four levels corresponding to the main stages of data processing and analysis (Figure 2). The first part is related to the identification of the available Data Input for the analysis: CORONA KH-4B, SPOT-1, Sentinel-2A MSI, and Sentinel-1A SAR C band satellite imagery. The second level corresponds to data processing, including a complete digital photogrammetric processing approach for CORONA images [74] and data correction and enhancement for multispectral and radar images, using ENVI and SNAP softwares. The third level refers to the preliminary results of corrected data. The fourth level corresponds to GIS data processing for LCLU feature delineation, based on image interpretation techniques on CORONA 4B and SPOT 1 images and semi-automatic classification of Sentinel-1 and Sentinel-2 decorrelated data stack. The final step refers to the integration of all datasets generated from previous steps, resized on the official limits of the municipality of Bucharest.
The photogrammetric processing of the historical CORONA KH-4 B images from 1968 started from two digitized film frames, available on Earth Explorer geoportal [28]. The non-metric camera model was set according to the technical parameters, while the exterior orientation was based on 31 Ground Control Points (GCP) collected on available orthophotos at 1:5,000 scale [75] and Digital Elevation Model (DEM) based on 1:25,000 Military Topographic Maps [76]. The final bundle-adjustment of Aerial Triangulation (AT) subset was accomplished by using Brown model, integrated in MATCH-AT Trimble and endorsing an accuracy of less than 0.5 meters [77]. The 16-parameter model describes five trapezoidal distortions in X and Y, three corrections for lens eccentricity, and three for radial distortion [78,79]. The resulted orthophotomosaic with 2 m spatial resolution was the support of photointerpretation of the LCLU classes for the single-date mapping from 1968. A manual feature extraction was preferred due to the limited spectral resolution and because the scanned film frame grain size strongly reduces the quality of the automatic image segmentation accuracy.
Satellite image processing focused on the preparation of SPOT-1, Sentinel-1, and Sentinel-2 imagery for data extraction.
SPOT-1 images were freely delivered for scientific use by CNES ISIS (Centre National D’Etudes Spatiales, Initiative for Space Innovative Standard) [80]. Images processed at level 2A came in two scenes panchromatic at 10 m spatial resolution and multispectral with Near Infrared, Red, and Green bands at 20 m spatial resolution. In order to generate a continuous coverage, both scenes were mosaicked for panchromatic and multispectral data. To ensure the resolution for multispectral image, we produced a Gram–Schmidt pansharpened combination [81]. As a result, the false color image at 10 m spatial resolution was the subject of visual interpretation for LCLU feature extraction.
The Sentinel-1 or S1 radar imagery was radiometrically calibrated in order to generate the backscatter intensity component image [82]. This followed a general workflow dealing with radar images processing [83], including the identification of radar data with the orbit and the terrain features (Figure 2). After the application of the Range-Doppler Terrain Correction followed by a speckle filtering, the backscatter image of the VV-VH polarization components was prepared for co-registering with the optical imagery S2 MSI at L2A processing level in terms of surface reflectance, at 10 m spatial resolution [84].
The image classification process based on the recent satellite Sentinel-1 and Sentinel-2 images from 2018 was performed after building a hybrid and highly decorrelated image stack dataset [63]. Training reference data for eight LCLU classes generated on the image dataset was analysed using the separability tool [85]. This explains the limited number of classes that can be extracted due to the similar spectral signature of some of the classes. For an increased accuracy of the classification, Support Vector Machine (SVM) algorithm with a radial basis function was used for the LCLU classes delineation [86,87].
Urban landscape change analysis and mapping used a GIS application for the harmonization of the three datasets to the built-up area coverage. For statistic approach and metric analysis, we introduced the Built-up Change Index (BCI) as a percentage difference of built-up area quantified for a reference surface of 1 hectare (1).
B C I ( % ) = B u i l t u p   c h a n g e R e f e r e n c e   s u r f a c e × 1 100
The analysis was performed in detail for three complementary case studies from the city of Bucharest, which illustrate the direction and magnitude of changes in each key area.

3. Results

The first results of this analysis are represented by the LCLU maps for the three reference periods: 1968, 1989, and 2018. The thematic classification of the S1-S2 layer was statistically compared with a set of random sample points representing ground truth data through an error matrix. The assessment of S1-S2 classification raster returned a high overall accuracy of 96% with a kappa coefficient of 0.942 for built-up area, which validates the result for the integration in the urban landscape analysis. These new products offer support for conducting spatial and statistical analysis for changes that occurred in the urban landscape.
The thematic map representing the past date LCLU (Figure 3) has a legend of 24 classes where the built-up area is the most detailed one, composed of nine different sub-classes, with a clear separation based on functional criteria. The map from 1968 shows an accurate radiography of a city at the beginning of a demographic explosion, with large parts covered by districts of houses with gardens, surrounded by large agricultural plots and forest patches. This past landscape preserves large semi-rural districts with houses and gardens to the border of the built-up area, and an emerging block of flat districts built on empty grounds, near intensive and systematic industrial areas developments [88,89].
LCLU map from 1989 reveals the image of a compact city with clearly delineated neighbourhoods of blocks of flats, most dominant in the western and eastern parts. The presence of small houses surrounded by gardens was rare in 1989. Moreover, large industrial areas appeared at the border of the city. The arable land at the outskirts decreased in the southern part. Also, there are many construction sites in the city center. Major hydrological works appeared, such as the construction of Morii Lake and Dîmbovița river regulation.
The most recent LCLU map (Figure 4), from 2018, highlights the residential areas of the block of flat districts or residential apartment buildings that surround the entire historical core and modern time development belts, partly replacing the former residential house with garden districts.
Overall, the recent LCLU map is the image of a compact city, with a regressive dynamic of green areas and empty grounds, with brownfields in the process of transformation after the real-estate boom of 2002–2008 [17]. A second set of results is represented by the maps showing the BCI for each time interval: 1968–1989 and 1989–2018 (Figure 5). This index is introduced in order to quantify the city’s fragmentation for each time period, with a focus on built-up classes. It offers a visual support for understanding the urban sprawl process, which occurred mostly at the periphery.
The period 1968–1989 illustrates the large-scale development of blocks of flats and factories replacing the house districts, marked by an increase of BCI index between +75 and +100%. The increase of built-up areas is the direct consequence of population growth from 1.36 million in 1966 up to 2.06 million in 1992 [12]. Population density almost doubled from 5900 inhabitants/km2 in 1966 to 9052 inhabitants/km2 in 1992.
On the other side, the map also reveals decreases of built-up area as a trend of a sustainable policy of urban development. The Palace of Parliament, also known as The House of the People, was constructed after 1984 and it is composed of a very large building and large open spaces around it [17]. It replaced two districts of houses with gardens—Uranus and partly Izvor [90] with a negative BCI of −75% to −100%.
The BCI map for the period 1989–2018 (Figure 5) highlights the conversion from arable land and gardens into large commercial or retail areas. Moreover, the map spatially explains the continuity of urban growth through the sprawl process, featuring the transition stage to the market-driven economy but mainly the real-estate development after 2001 [91]. This is related to the highest built-up area development rates between +75% and +100%. Here, the arable land or gardens and the large brownfields were replaced by large commercial or retail areas, after they were abandoned due to a high rate of bankruptcy at the beginning of 1990s [92].
On the other hand, a decrease of built-up area of −50% is reported in the eastern part of the city: This is the place of the declared protected area of Văcărești Nature Reserve of 190 hectares [93]. Similar trends are observed along the lakes of Colentina River, a special protection area of patrimony. Giulești district, in the north-western part, is another example of brownfields that were located in less attractive places for real-estate development due to the large railway transportation infrastructure in a railway junction area, degraded and partly dismantled after 1989 (marshalling yards, locomotive depots, rolling stock maintenance and other specific facilities).
The transformation of the urban structures is presented from a quantitative perspective in Figure 6 and Figure 7a,b, showing the intensity of changes.
The change diagrams focus on the key urban LCLU classes, selected after simplification of the first datasets (1968 and 1989). These were in fact the classes used for supervised classification training on the recent image stack (2018). Table 2 contains a description of these classes according to the ground truth data and providing an argument of their selection. For example, it is easy to notice the replacement of the green zones by built-up areas, in the context of the economic development of Bucharest and the demographic increase, while industrial areas combine with retail areas for the 1989 and 2018 datasets in terms of spectral signatures.
According to the description provided, it was essential to limit the number of key classes in order to facilitate the production of more relevant results and to map the change in a simple and accurate cartographic formula. Different uncertainties in class separation and description were taken into account but they were investigated when evaluating and interpreting the results.

4. Discussion

Our study brings an added value that refers to the increased number of classes used for the thematic classification. Also, the analysis covers a larger temporal frame of 5 decades compared to other similar studies that were conducted for Bucharest [94,95,96].
The approach started from the four satellite images collection featured by a variable spatial resolution of 2 meters, 10 meters and 20 meters. First of all, the resolution of the multispectral SPOT image was resampled based on panchromatic band at 10 meters in a pansharpening product, and also the 20-meter Sentinel-2 bands were resampled at 10 meters. The target was also to use a similar mapping unit size when extracting the LCLU classes with a minimum mapping unit of 1 ha, for each reference year: 1968, 1989 and 2018. We used a semantic generalization by joining of land use classes that refer to the same land cover class, in order to integrate the different number of classes extracted for this diachronic dataset.
Multi-date satellite imagery synergy is a reality built on a complex methodology for data calibration and geographical information production. A combination of photogrammetric and digital image processing workflows opened an opportunity to generate three LCLU data coverages for the Bucharest urban area.
Two aspects are essential in this context. First, spatial and spectral resolution differences between the images did not allow a unique approach for time series data production. It is the reason why we combined the visual interpretation of 1968 and 1989 images with the thematic classification of a decorrelated image stack, where multispectral and radar signatures helped the mapping process of the key LCLU classes. The second aspect is the topological models building and their evaluation. In this context, the RMSE computation for the 1968 orthoimage mosaic was done and a map of errors distribution was interpreted before the interpretation of the dataset for LCLU layer production.
In order to assess the accuracy of 2018 image stack classification, the results were compared with the ground truth data systematically extracted from the 0.5 m resolution color orthophotos (ANCPI—National Agency for Cadastre and Land Registration, 2012) using a special confusion matrix (Table 3). An independent set of 346 polygons was produced and equally distributed to the evaluated LCLU classes. The overall accuracy of 94% and the kappa coefficient of 0.92 indicate a very good result in classification. The confusion matrix indicates the accuracy per each class and the same low results can be noticed for street-related class and for the retail area class, imposed by the almost similar spectral response with other classes.
A noticeable aspect that resulted from the comparison between 1968 and 1989 thematic datasets is the disappearance of the high fragmentation of gardening plots around the built-up area, especially in the eastern and southern parts. The problem of land use classes still occurred for the year 1968. It is true that it is not a semi-urban, semi-rural or semi-wild area. This is reflected by the classes introduced. We use the term “mixed class” in order to define a particular pattern featured by small houses surrounded by vegetable gardens, vineyards, and orchards, where people develop a survival-like small agriculture, which is a feature of Balkan cities. It is not a “class like category”, because this was a spatial projection of a traditional way of life, adapted to local conditions that survived over the centuries and was systematically destroyed during the communism period by house and farm demolition, replaced by block ensembles/industrial plants and later by the real-estate market pressures (an intensive one after the year 2001, when ownership recovery was helped by a new special legislation).
Other important changes are related to hydrotechnical works, in order to protect and control water level during floods [96], such as the channelization of Dâmbovița River and the building of Morii Lake reservoir in the western part of the city.
Among the mapped classes, the built-up class was selected to highlight the relevant urban changes. City development policies introduced new zonation of the functional areas. For example, in 1968, the limits of commercial and retail areas were not clearly defined, while residential areas were not combined with other structures like houses and gardens or industrial districts [97].
For an adequate understanding of the urban structural changes in the city of Bucharest, it is necessary to explore some complementary case studies to reveal specific transformations. Three selected areas are presented in the context of built-up area change rate levels, together with a comparative survey of urban configuration in 1968, 1989, and 2018. From this point of view, the selection focuses on:
1. The Palace of Parliament area, downtown of Bucharest, with a dramatic replacement of an entire urban area by an extremely large building that was supposed to be the emblem of the communist regime in Romania.
2. Drumul Taberei neighbourhood, in the western part of the city, where rural-like structures were replaced by a totally new district with blocks of flats and industrial enterprises, and later by new developments related to the market-driven economy and real-estate investments.
3. Băneasa area, where arable land and orchards were replaced by office parks, logistics, and retail areas.
These case studies were selected because of their complementary characteristics. They summarise the most common key aspects that Bucharest faced in the last decades: (i) city reshape for the Parliament House area, and (ii) massive conversion from garden houses into blocks of flats in Drumul Taberei district, (iii) urban sprawl and sub-urbanisation in Băneasa district.
The Palace of Parliament area (Figure 8) is probably the most spectacular urban landscape change not only in Romania, but in Central Eastern Europe, well-known in literature as the “Vanished Bucharest” [98]. Here, a rate of BCI between 50 and 100% is visible along the axis, which splits the historical centre into two parts. It also illustrates the total replacement of the old traditional urban structures, composed of small houses and gardens, together with a dense and complex street network, which resulted from the unplanned development before the 19th century (Figure 8a,b), with socialist-style public building of the Palace of Parliament and the blocks of flats along Dâmbovița River (Figure 8c) and some green areas [15,99].
According to the urban planners’ documentations [100], a total surface of 580 hectares was affected by demolition works after 1984, although the projects started after the earthquake of 4 March 1977. The qualitative change of these structure can be explained by the rising from the ground of urban identity elements: 1 monastery, 10 churches, 5 moved churches to new sites, 1 stadium, and 1 inn, while the new block-of-flats type building along the avenue hide other remained historical monuments of an outstanding value [101]. This case study is an example for erasing the city’s legacy and the forced introduction of a totalitarian urban policy, which is contrary to the adequate urban LCLU features [98]. The area still has a real problem related to urban re-integration, although the palace is now part of the administrative buildings ensemble, as it is the house of the Romanian Parliament [18,102]. In the vicinity of the Palace of Parliament, there is the emerging building of the new Romanian Orthodox Church cathedral (started 2010), which is to be finished in the coming years.
Drumul Taberei district area is a typical urban landscape change that occurred in the city of Bucharest at the end of 1960s and the beginning of 1970s. It is an example of totally new urban structures superposed on former rural structures, as the town’s built-up area was for a long time in a close relationship with the neighbouring agricultural area. This finding is also confirmed by other studies [95] and is related to the unplanned development until the 19th century in the context of an unfortified urban centre. Figure 9 highlights the total replacement of houses and farms with gardens, orchards, and vineyards, which resulted from an unplanned configuration (Figure 9a), with structures of a particular physiognomy, after the relaxation of the urban policies of the 50s and the increasing respect to the new urban regulations emphasized by Athens Charter [103]. This urban district with block-of-flats and green areas, with service spaces and other, was projected as a new township for 100,000 inhabitants in the context of demographic increase and intensive industrialization. However, nowadays, it hosts about 400,000 people (Figure 9b). In the last 10 years, new urban structures with the same physiognomy appeared here even if the real-estate market value decreased (Figure 9c), as a result of emerging superficial investments and lack of modernization of urban transport.
Băneasa area (Figure 10), in the northern part of the city, is a typical example of urban structures in the globalization trend. Here, the urban sprawl process can be easily noticed and driven by economic-market rules and lack of planning regulation. This situation is common in the majority of the former socialist countries [14]. This area also integrates the Aurel Vlaicu International Airport, which has been operating since 1920s and nowadays is enclosed in the built-up area. Besides the airport, this area is crossed by the main national and European highway connecting Bucharest to Ploiești and Brașov. This main route is used daily by residents who commute between Bucharest and Ploiesti and by people who travel during the weekend. At the same time, it represents the main road access to the main international airport of Romania at Otopeni (Henri Coandă). It is the reason for a lot of infrastructure investments for road traffic improvement during the last two decades. Therefore, this area is characterized by a spectacular urban sprawl and development, mainly after the 1990s and more intensively after 2001 [104]. The urban land ownership recovery starting from 1995, as well as the release of the latest Bucharest Urban Master Plan, contributed to the current landscape as it is today: complex built-up area and Băneasa forest zone. This area is the largest complex urban development in Romania after 1990, where 221 hectares were covered by a combination of business parks and residential districts with villas, green areas, and a large commercial space. This represents an increase of the urban landscape change rate between +75 and +100%, around Băneasa airport (Figure 10c).
Bański [105], found that the greatest loss of agricultural lands took place in the periphery of large urban area, and we can confirm this situation for Bucharest city. Grigorescu et al. [106] underline the conversion of arable land to urban, industrial, and commercial uses after 1990, especially around the capital cities such as Bucharest.
These three case studies confirm the complex pattern of the urban landscape change in Bucharest. This is an illustration of the incoherent urban development policy of a city and the chaotic transition from a communist regime to a capitalist market driven economy. The planning strategies occurred after 1831 legislation were based on radial and concentric configuration of the built-up area within the city. However, the systematic planning was possible only after the first master plans from 1921 and 1935, which tried to draw the main street network and to differentiate between the urban structures and the development areas [107,108]. This normal historical evolution was interrupted by the socialist urban planning after 1948, where private planning strategies were oriented towards the demographic and the industrial growth of the town between 1950 and 1980, followed by a severe transformation of the city center and the loss of local identity between 1980 and 1989 [92]. New financial opportunities were exploited as soon as the ownership regime changed after the 1990s. These changes are visible especially after 2000 when a new Master Plan was released and investments increased for expansion of built-up area based on conversion of brownfields.

5. Conclusions

Mapping the urban landscape changes with a span of about 50 years was accomplished using the synergy between CORONA, Spot and Sentinel remote sensing imagery. The declassified dataset from 1968 allowed the manual extraction of the former structure of the urban landscape, same as for the reference year 1989. The multisensor approach based on the recent radar and optical Sentinel images enabled the semi-automatic production of the urban land cover features. The data is integrated and harmonized, thus obtaining a very accurate result for a city with a complex and heterogeneous urban pattern.
The main progress element of the analysis is the detailed GIS mapping and the statistical quantification of the key urban structures change for the reference 5-decade period, which is a step forward from the existing approaches. Its dynamics draw a projection of the urban landscape transformation and urban system dynamics with all the background systemic connection, from land to social life and economy. The urban imprint is actually a key feature strictly connected to the LCLU data layers.
The metric analysis using BCI index quantifies the changes that occurred in a synthetic way. The analysis emphasizes the urban sprawl of Bucharest within the conversion of the agricultural lands to built-up areas. The systematic replacement of the traditional land use pattern composed by small houses and agriculture gardens with the new block of flats met the necessities of the growing population. However, the uncontrolled real-estate developments have many social effects including life quality in the city of Bucharest. New problems were generated, such as low living comfort, lack of green space, lack of parking spaces, and deficiency of the urban and transport infrastructures. The landscape change is also related to the replacement of historical districts with socialist buildings. This evolution led to the fragmentation of the historical area and the decreasing of the architectural patrimony density.
Bucharest, unlike other Eastern European capital cities, had a transformation based on the destruction of a lot of semi-rural zones but especially on urban areas with historical architecture of a real cultural value, in order to free the land for the construction of new government buildings that correspond to the ideology of communist propaganda (after 1984). This is also reflected in the city evolution of the post-communist period that failed to integrate these ensembles in a typical urban structure. We believe that not only the shape, but especially the model of spatial distribution by extending to the periphery of the new ensembles, is a distinctive feature for Eastern European cities.

Author Contributions

Conceptualization, C.N. and M.V.; methodology, C.N., M.V., I.C.; software, C.N., M.V., I.C., L.T., B.O.; validation, C.N., M.V., I.C., B.-A.M.; writing—original draft preparation, C.N., M.V., I.C., B.-A.M.; writing—review and editing, C.N., M.V., I.C., B.-A.M., L.T., B.O.; supervision, C.N.; funding acquisition, C.N. All authors have read and agreed to the published version of the manuscript.

Funding

The APC was funded by University of Bucharest.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

We thanks to the Centre National d’Etudes Spatiales (CNES) for providing us for free the SPOT 1 satellite image for scientific use under the program Incitation à l’utilisation Scientifique des Images Spot (ISIS).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cegielska, K.; Noszczyk, T.; Kukulska, A.; Szylar, M.; Hernik, J.; Dixon-Gough, R.; Jombach, S.; Valánszki, I.; Kovács, K.F. Land use and land cover changes in post-socialist countries: Some observations from Hungary and Poland. Land Use Policy 2018, 78, 1–18. [Google Scholar] [CrossRef]
  2. Fan, C.; Myint, S. A comparison of spatial autocorrelation indices and landscape metrics in measuring urban landscape fragmentation. Landsc. Urban Plan. 2014, 121, 117–128. [Google Scholar] [CrossRef]
  3. Schwarz, N. Urban form revisited—Selecting indicators for characterising European cities. Landsc. Urban Plan. 2010, 96, 29–47. [Google Scholar] [CrossRef]
  4. Norton, B.A.; Evans, K.L.; Warren, P.H. Urban biodiversity and landscape ecology: Patterns, processes and planning. Curr. Landsc. Ecol. Rep. 2016, 1, 178–192. [Google Scholar] [CrossRef] [Green Version]
  5. Mascarenhas, A.; Haase, D.; Ramos, T.; Santos, R. Pathways of demographic and urban development and their effects on land take and ecosystem services: The case of Lisbon Metropolitan Area, Portugal. Land Use Policy 2019, 82, 181–194. [Google Scholar] [CrossRef]
  6. Jepsen, M.R.; Kuemmerle, T.; Müller, D.; Erb, K.; Verburg, P.H.; Haberl, H.; Vesterager, J.P.; Andric, M.; Antrop, M.; Austrheim, G.; et al. Transitions in European land-management regimes between 1800 and 2010. Land Use Policy 2015, 49, 53–64. [Google Scholar] [CrossRef]
  7. Hersperger, A.M.; Bürgi, M.; Wende, W.; Bacău, S.; Grădinaru, S.R. Does landscape play a role in strategic spatial planning of European urban regions? Landsc. Urban Plan. 2020, 194, 103702. [Google Scholar] [CrossRef]
  8. Herold, M.; Couclelis, H.; Clarke, K. The role of spatial metrics in the analysis and modeling of urban land use change. Comput. Environ. Urban Syst. 2005, 29, 369–399. [Google Scholar] [CrossRef]
  9. Ngie, A.; Abutaleb, K.; Fethi, A.; Taiwo, O.-J. Spatial modelling of urban change using satellite remote sensing. In Life in a Changing Urban Landscape, Proceedings of the IGU Urban Geography Commission (Urban Challenges in a Complex World), Krakow, Poland, 15–22 August 2014; Nico, K., Ronnie, D., Gustav, V., Eds.; University of Johannesburg: Johannesburg, South Africa, 2014; pp. 3–12. [Google Scholar]
  10. Mesev, V. Remotely-Sensed Cities; CRC Press: Boca Raton, FL, USA, 2003. [Google Scholar]
  11. Aplin, P. Comparison of simulated IKONOS and SPOT HRV imagery for classifying urban areas. In Remotely-Sensed Cities; CRC Press: Boca Raton, FL, USA, 2020. [Google Scholar]
  12. The Site of Romanian National Institute for Statistics. Available online: http://statistici.insse.ro/ (accessed on 12 June 2020).
  13. Stanilov, K. Housing trends in Central and Eastern European cities during and after the period of transition. In Cities between Competitiveness and Cohesion; Springer Science and Business Media LLC: Berlin/Heidelberg, Germany, 2007; Volume 92, pp. 173–190. [Google Scholar]
  14. Garcia-Ayllon, S. Urban transformations as indicators of economic change in post-communist Eastern Europe: Territorial diagnosis through five case studies. Habitat Int. 2018, 71, 29–37. [Google Scholar] [CrossRef]
  15. Badiu, D.L.; Onose, D.A.; Niță, M.R.; Lafortezza, R. From “red” to green? A look into the evolution of green spaces in a post-socialist city. Landsc. Urban Plan. 2019, 187, 156–164. [Google Scholar] [CrossRef]
  16. Marin, V.; Chelcea, L. The many (still) functional housing estates of Bucharest, Romania: A viable housing provider in Europe’s densest capital city. In The Life and Afterlife of Gay Neighborhoods; Springer Science and Business Media LLC: Berlin/Heidelberg, Germany, 2018; pp. 167–190. [Google Scholar]
  17. Mihai, B.; Nistor, C.; Simion, G. Post-socialist urban growth of Bucharest, Romania—a change detection analysis on Landsat imagery (1984–2010). Acta Geogr. Slov. 2015, 55. [Google Scholar] [CrossRef] [Green Version]
  18. Ianoş, I.; Sorensen, A.; Merciu, C. Incoherence of urban planning policy in Bucharest: Its potential for land use conflict. Land Use Policy 2017, 60, 101–112. [Google Scholar] [CrossRef]
  19. Tudor, C.A.; Iojă, I.C.; Rozylowicz, L.; Pǎtru-Stupariu, I.; Hersperger, A.M. Similarities and differences in the assessment of land-use associations by local people and experts. Land Use Policy 2015, 49, 341–351. [Google Scholar] [CrossRef]
  20. Boentje, J.P.; Blinnikov, M. Post-Soviet forest fragmentation and loss in the Green Belt around Moscow, Russia (1991–2001): A remote sensing perspective. Landsc. Urban Plan. 2007, 82, 208–221. [Google Scholar] [CrossRef]
  21. Rufat, S.; Marcińczak, S. The equalising mirage? Socioeconomic segregation and environmental justice in post-socialist Bucharest. Neth. J. Hous. Environ. Res. 2020, 35, 917–938. [Google Scholar] [CrossRef]
  22. Suditu, B. Bucureștiul în locuințe și locuitori (Bucharest in housing and inhabitant); Editura Compania: Bucharest, Romania, 2016; ISBN 6066800167. [Google Scholar]
  23. Ianos, I.; Sirodoev, I.; Pascariu, G. Land-use conflicts and environmental policies in two post-socialist urban agglomerations: Bucharest and Chişinǎu. Carpathian J. Earth Environ. Sci. 2012, 7, 125–136. [Google Scholar]
  24. Shaw, B.J.; van Vliet, J.; Verburg, P.H. The peri-urbanization of Europe: A systematic review of a multifaceted process. Landsc. Urban Plan. 2020, 196, 103733. [Google Scholar] [CrossRef]
  25. EEA. European Environment Agency Urban Sprawl in Europe—The Ignored Challenge; Wiley: Hoboken, NJ, USA, 2006; ISBN 9781405139175. [Google Scholar]
  26. Oueslati, W.; Alvanides, S.; Garrod, G. Determinants of urban sprawl in European cities. Urban Stud. 2015, 52, 1594–1614. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Hamidi, S.; Ewing, R.; Tatalovich, Z.; Grace, J.B.; Berrigan, D. Associations between Urban Sprawl and Life Expectancy in the United States. Int. J. Environ. Res. Public Health 2018, 15, 861. [Google Scholar] [CrossRef] [Green Version]
  28. USGS. USGS EarthExplorer; USGS: Reston, VI, USA, 2015.
  29. Yang, L.; Huang, C.; Homer, C.G.; Wylie, B.; Coan, M.J. An approach for mapping large-area impervious surfaces: Synergistic use of Landsat-7 ETM+ and high spatial resolution imagery. Can. J. Remote Sens. 2003, 29, 230–240. [Google Scholar] [CrossRef]
  30. Dittrich, A.; Buerkert, A.; Brinkmann, K. Assessment of land use and land cover changes during the last 50 years in oases and surrounding rangelands of Xinjiang, NW China. J. Agric. Rural Dev. Trop. Subtrop. 2010, 111, 129–142. [Google Scholar]
  31. Li, M.; Mao, L.; Shen, W.; Liu, S.; Wei, A. Change and fragmentation trends of Zhanjiang mangrove forests in southern China using multi-temporal Landsat imagery (1977–2010). Estuarine Coast. Shelf Sci. 2013, 130, 111–120. [Google Scholar] [CrossRef]
  32. Brinkmann, K.; Schumacher, J.; Dittrich, A.; Kadaore, I.; Buerkert, A. Analysis of landscape transformation processes in and around four West African cities over the last 50 years. Landsc. Urban Plan. 2012, 105, 94–105. [Google Scholar] [CrossRef]
  33. Fekete, A. Urban and rural landslide hazard and exposure mapping using landsat and corona satellite imagery for tehran and the Alborz Mountains, Iran. AIMS Geosci. 2017, 3, 37–66. [Google Scholar] [CrossRef]
  34. Saleem, A.; Corner, R.; Awange, J. On the possibility of using CORONA and Landsat data for evaluating and mapping long-term LULC: Case study of Iraqi Kurdistan. Appl. Geogr. 2018, 90, 145–154. [Google Scholar] [CrossRef]
  35. Ratcliffe, I.; Henebry, G. Using Declassification Intelligence Satellite Pictures with Quickbird Imagery to Study Urban Land Cover Dynamics: A Case Study from Kazakhstan. Annu. Proc. ASPRS 2004, 198, 1–10. [Google Scholar]
  36. Mucsi, L.; Liska, C.M.; Henits, L.; Tobak, Z.; Csendes, B.; Nagy, Z. The evaluation and application of an urban land cover map with image data fusion and laboratory measurements. Hung. Geogr. Bull. 2017, 66, 145–156. [Google Scholar] [CrossRef] [Green Version]
  37. Dong, J.-J.; Chen, N.-H.; Ma, Y.-H.; Chen, J.-Y.; Jin-Jin, D.; Ning-Hua, C.; Yi-Hang, M.; Jian-Yu, C. Land use change and information extraction of rural residential land based on Corona KH-4B Imagery. In Proceedings of the 2012 2nd International Conference on Remote Sensing, Environment and Transportation Engineering; Institute of Electrical and Electronics Engineers (IEEE), Nanjing, China, 1–3 June 2012; pp. 1–3. [Google Scholar]
  38. Scardozzi, G. Multitemporal satellite images for knowledge of the assyrian capital cities and for monitoring landscape transformations in the upper course of Tigris River. Int. J. Geophys. 2011, 2011, 1–17. [Google Scholar] [CrossRef] [Green Version]
  39. ESA Copernicus Programme. Copernicus 2017; ESA: Paris, France, 2017. [Google Scholar]
  40. Aschbacher, J.; Milagro-Pérez, M.P. The European Earth monitoring (GMES) programme: Status and perspectives. Remote Sens. Environ. 2012, 120, 3–8. [Google Scholar] [CrossRef]
  41. Copernicus. Copernicus Open Access Hub; Copernicus: Brussels, Belgium, 2018. [Google Scholar]
  42. Shahtahmassebi, A.R.; Song, J.; Zheng, Q.; Blackburn, G.A.; Wang, K.; Huang, L.Y.; Pan, Y.; Moore, N.; Shahtahmassebi, G.; Haghighi, R.S.; et al. Remote sensing of impervious surface growth: A framework for quantifying urban expansion and re-densification mechanisms. Int. J. Appl. Earth Obs. Geoinf. 2016, 46, 94–112. [Google Scholar] [CrossRef] [Green Version]
  43. Santos, N.D.; Gonçalves, G. Remote Sensing Applications Based on Satellite Open Data (Landsat8 and Sentinel-2). In Proceedings of the Conferência Nacional de Geodecisão, Barreiro, Portugal, 15–16 May 2014; pp. 1–10. [Google Scholar]
  44. Drusch, M.; Del Bello, U.; Carlier, S.; Colin, O.; Fernandez, V.; Gascon, F.; Hoersch, B.; Isola, C.; Laberinti, P.; Martimort, P.; et al. Sentinel-2: ESA’s optical high-resolution mission for GMES operational services. Remote Sens. Environ. 2012, 120, 25–36. [Google Scholar] [CrossRef]
  45. Topaloğlu, R.; Sertel, E.; Musaoglu, N. Assessment of classification accuracies of Sentinel-2 And Landsat-8 data for land cover/use mapping. In Proceedings of the XXIII ISPRS Congress, Prague, Czech Republic, 12–19 July 2016; Volume XLI-B8. [Google Scholar]
  46. Rujoiu-Mare, M.-R.; Olariu, B.; Mihai, B.-A.; Nistor, C.; Săvulescu, I. Land cover classification in Romanian Carpathians and Subcarpathians using multi-date Sentinel-2 remote sensing imagery. Eur. J. Remote Sens. 2017, 50, 496–508. [Google Scholar] [CrossRef] [Green Version]
  47. Novelli, A.; Aguilar, M.A.; Nemmaoui, A.; Aguilar, F.J.; Tarantino, E. Performance evaluation of object based greenhouse detection from Sentinel-2 MSI and Landsat 8 OLI data: A case study from Almería (Spain). Int. J. Appl. Earth Obs. Geoinf. 2016, 52, 403–411. [Google Scholar] [CrossRef] [Green Version]
  48. Ghanea, M.; Moallem, P.; Momeni, M. Building extraction from high-resolution satellite images in urban areas: Recent methods and strategies against significant challenges. Int. J. Remote Sens. 2016, 37, 5234–5248. [Google Scholar] [CrossRef]
  49. Lefebvre, A.; Sannier, C.; Corpetti, T. Monitoring urban areas with Sentinel-2A Data: Application to the update of the copernicus high resolution layer imperviousness degree. Remote Sens. 2016, 8, 606. [Google Scholar] [CrossRef] [Green Version]
  50. Sertel, E.; Akay, S.S. High resolution mapping of urban areas using SPOT-5 images and ancillary data. Int. J. Environ. Geoinf. 2015, 2, 63–76. [Google Scholar] [CrossRef] [Green Version]
  51. Skupinski, G.; BinhTran, D.; Weber, C. Les images satellites Spot multi-dates et la métrique spatiale dans l’étude du changement urbain et suburbain—Le cas de la basse vallée de la Bruche (Bas-Rhin, France). Cybergeo 2009. [Google Scholar] [CrossRef]
  52. Martín, L.; Howarth, P. Change-detection accuracy assessment using SPOT multispectral imagery of the rural-urban fringe. Remote Sens. Environ. 1989, 30, 55–66. [Google Scholar] [CrossRef]
  53. da Costa, S.M. Spot imagery for classification of urban land use: A comparison with Landsat TM imagery—A study of Belo Horizonte area. In Proceedings of the ISPRS Congress 17, Washington, DC, USA, 2–14 August 1992; pp. 575–582. [Google Scholar]
  54. Abdikan, S.; Sanli, F.B.; Ustuner, M.; Calò, F. Land cover mapping using Sentinel-1 sar data. ISPRS Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2016, B7, 757–761. [Google Scholar] [CrossRef]
  55. Dostálová, A.; Hollaus, M.; Milenković, M.; Wagner, W. Forest area derivation from Sentinel-1 data. ISPRS Ann. Photogramm. Remote Sens. Spat. Inf. Sci. 2016, 7, 227–233. [Google Scholar] [CrossRef] [Green Version]
  56. Georgescu, F.A.; Tanase, R.; Datcu, M.; Raducanu, D. Patch-Based Image Classification for Sentinel-1 and Sentinel-2 Earth Observation Image Data Products; European Space Agency: Paris, France, 2016; Volume SP-740. [Google Scholar]
  57. Badea, A.C.; Badea, G. Advanced of Identifying Urban Footprint using Sentinel-1. In Embracing Our Smart World Where the Continents Connect: Enhancing The Geospatial Maturity of Societies, Proceeding of the FIG Congress 2018, Istambul, Turkey, 6–11 May 2018; FIG: Copenhagen, Denmark, 2018. [Google Scholar]
  58. Serco Italia. Urban Clasification with Sentinel-1. Available online: https://rus-copernicus.eu/portal/ (accessed on 1 December 2018).
  59. Jacob, A.; Ban, Y. Sentinel-1A SAR data for global urban mapping: Preliminary results. In Proceedings of the 2015 IEEE International Geoscience and Remote Sensing Symposium (IGARSS), Milan, Italy, 26–31 July 2015; pp. 1179–1182. [Google Scholar]
  60. Vigiak, O.; Malagó, A.; Bouraoui, F.; Vanmaercke, M.; Obreja, F.; Poesen, J.; Habersack, H.; Fehér, J.; Grošelj, S. Modelling sediment fluxes in the Danube River Basin with SWAT. Sci. Total. Environ. 2017, 599–600, 992–1012. [Google Scholar] [CrossRef] [PubMed]
  61. Ban, Y.; Webber, L.; Gamba, P.; Paganini, M. EO4Urban: Sentinel-1A SAR and Sentinel-2A MSI data for global urban services. In Proceedings of the 2017 Joint Urban Remote Sensing Event (JURSE), Dubai, United Arab Emirates, 6–8 March 2017; pp. 1–4. [Google Scholar] [CrossRef]
  62. Pesaresi, M.; Corbane, C.; Julea, A.; Florczyk, A.J.; Syrris, V.; Soille, P. Assessment of the added-value of Sentinel-2 for detecting built-up areas. Remote Sens. 2016, 8, 299. [Google Scholar] [CrossRef] [Green Version]
  63. Clerici, N.; Calderón, C.A.V.; Posada, J.M. Fusion of Sentinel-1A and Sentinel-2A data for land cover mapping: A case study in the lower Magdalena region, Colombia. J. Maps 2017, 13, 718–726. [Google Scholar] [CrossRef] [Green Version]
  64. Wenbo, X.; Bingfang, W.; Yichen, T.; Jianxi, H.; Yong, Z. Synergy of multitemporal Radarsat SAR and Landsat ETM data for extracting agricultural crops structure. In Proceedings of the IEEE International Geoscience and Remote Sensing Symposium, 2004. IGARSS ’04, Anchorage, AK, USA, 20–24 September 2004; Volume 6, pp. 4073–4076. [Google Scholar]
  65. Peters, J.; Van Coillie, F.; Westra, T.; De Wulf, R. Synergy of very high resolution optical and radar data for object-based olive grove mapping. Int. J. Geogr. Inf. Sci. 2011, 25, 971–989. [Google Scholar] [CrossRef]
  66. Van Tricht, K.; Gobin, A.; Gilliams, S.; Piccard, I. Synergistic Use of Radar Sentinel-1 and Optical Sentinel-2 Imagery for Crop Mapping: A Case Study for Belgium. Remote Sens. 2018, 10, 1642. [Google Scholar] [CrossRef] [Green Version]
  67. Grădinaru, S.R.; Kienast, F.; Psomas, A. Using multi-seasonal Landsat imagery for rapid identification of abandoned land in areas affected by urban sprawl. Ecol. Indic. 2019, 96, 79–86. [Google Scholar] [CrossRef]
  68. Ianoş, I.; Pascariu, G.; Sîrodoev, I.; Henebry, G. Divergent patterns of built-up urban space growth following post-socialist changes. Urban Stud. 2016, 53, 3172–3188. [Google Scholar] [CrossRef]
  69. Grivei, A.-C.; Radoi, A.; Datcu, M. Land cover change detection in Satellite Image Time Series using an active learning method. In Proceedings of the 2017 9th International Workshop on the Analysis of Multitemporal Remote Sensing Images (MultiTemp), Brugge, Belgium, 27–29 May 2017; pp. 1–4. [Google Scholar]
  70. Sandric, I.; Mihai, B.; Savulescu, I.; Suditu, B.; Chitu, Z. Change detection analysis for urban development in Bucharest-Romania, using high resolution satellite imagery. In Proceedings of the Urban Remote Sensing Joint Event 2007, Paris, France, 11–13 April 2007; pp. 1–8. [Google Scholar]
  71. Noaje, I.; Sion, I.G. Environmental changes analysis in Bucharest city using Corona, spot HRV and IKONOS images. ISPRS Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2012, B7, 329–334. [Google Scholar] [CrossRef] [Green Version]
  72. Tang, J.; Wang, L.; Yao, Z. Analyses of urban landscape dynamics using multi-temporal satellite images: A comparison of two petroleum-oriented cities. Landsc. Urban Plan. 2008, 87, 269–278. [Google Scholar] [CrossRef]
  73. Mihailescu, V. Evoluția geografică a unui oraș—București (The geographical evolution of a city—Bucharest); PAIDEIA: Bu-charest, Romania, 2003; ISBN 9789735961237. [Google Scholar]
  74. Mihai, B.; Nistor, C.; Toma, L.; Săvulescu, I. High resolution landscape change analysis with CORONA KH-4B Imagery. A case study from iron gates reservoir area. Procedia Environ. Sci. 2016, 32, 200–210. [Google Scholar] [CrossRef] [Green Version]
  75. ANCPI National Agency for Cadaster and Land Registration. Geoportal INIS Viewer; ANCPI: Bucharest, Romania, 2018. [Google Scholar]
  76. DTM. Military Mapping Directorate. Available online: www.geomil.ro (accessed on 20 August 2003).
  77. Brown, D. The bundle adjustment: Progress and prospects. Int. Arch. Photogramm. 1976, 21, 33. [Google Scholar]
  78. Sigle, M.; Heuchel, T. Match-At: Recent developments and performance. Photogrammetric Week 2001, 1, 189–194. [Google Scholar]
  79. Pang, S.; Sun, M.; Hu, X.; Zhang, Z. SGM-based seamline determination for urban orthophoto mosaicking. ISPRS J. Photogramm. Remote Sens. 2016, 112, 1–12. [Google Scholar] [CrossRef]
  80. Centre National d’Etudes Spatiales. Available online: https://isis-mission.cnes.fr/en/ISIS/index.htm (accessed on 12 November 2019).
  81. Alparone, L.; Aiazzi, B.; Baronti, S.; Garzelli, A. Remote Sensing Image Fusion; CRC Press: Boca Raton, FL, USA, 2015. [Google Scholar]
  82. Bourbigot, M.; Johnsen, H.; Piantanida, R. Sentinel-1 Product Definition; ESA: Paris, France, 2016. [Google Scholar]
  83. Stepanian, P.M.; Chilson, P.B.; Kelly, J.F. An introduction to radar image processing in ecology. Methods Ecol. Evol. 2014, 5, 730–738. [Google Scholar] [CrossRef] [Green Version]
  84. Zuhlke, M.; Fomferra, N.; Brockmann, C.; Peters, M.; Veci, L.; Malik, J.; Regner, P. SNAP (Sentinel Application Platform) and the ESA Sentinel 3 Toolbox. In Sentinel-3 for Science Workshop; The European Space Agency: Paris, France, 2015; Volume 12, p. 21. ISBN 9789292212988. [Google Scholar]
  85. Aplin, P. Image Analysis, Classification and Change Detection in Remote Sensing, with algorithms for ENVI/IDL. Int. J. Geogr. Inf. Sci. 2009, 23, 129–130. [Google Scholar] [CrossRef]
  86. Wessel, M.; Brandmeier, M.; Tiede, D. Evaluation of different machine learning algorithms for scalable classification of tree types and tree species based on sentinel-2 Data. Remote Sens. 2018, 10, 1419. [Google Scholar] [CrossRef] [Green Version]
  87. Noi, P.T.; Kappas, M. Comparison of Random Forest, k-Nearest Neighbor, and Support Vector Machine Classifiers for Land Cover Classification Using Sentinel-2 Imagery. Sensors 2017, 18, 18. [Google Scholar] [CrossRef] [Green Version]
  88. Purcaru, G. The Vrancea, Romania, earthquake of March 4, 1977—A quite successful prediction. Phys. Earth Planet. Inter. 1979, 18, 274–287. [Google Scholar] [CrossRef]
  89. Armaş, I.; Gavriş, A. Census-based social vulnerability assessment for Bucharest. Procedia Environ. Sci. 2016, 32, 138–146. [Google Scholar] [CrossRef] [Green Version]
  90. Danta, D. Ceausescu’s Bucharest. Geogr. Rev. 1993, 83, 170. [Google Scholar] [CrossRef]
  91. O’Neill, B. The political agency of cityscapes. J. Soc. Archaeol. 2009, 9, 92–109. [Google Scholar] [CrossRef]
  92. Stan, A. Urban expansion in Bucharest, after 1990: Errors and benefits. In Planning Capital Cities: Belgrade, Bucharest, Sofia; Doytchinov, G., Dukić, A., Ioniță, C., Eds.; Verlag der Technischen Universität Graz: Graz, Austria, 2015; pp. 224–234. ISBN 9783851253986. [Google Scholar]
  93. Manea, G.; Matei, E.; Vijulie, I.; Tîrlă, L.; Cuculici, R.; Cocoş, O.; Tişcovschi, A. Arguments for Integrative Management of Protected Areas in the Cities—Case Study in Bucharest City. Procedia Environ. Sci. 2016, 32, 80–96. [Google Scholar] [CrossRef] [Green Version]
  94. Nae, M.; Turnock, D. The new Bucharest: Two decades of restructuring. Cities 2011, 28, 206–219. [Google Scholar] [CrossRef]
  95. Aldea, M.; Petrescu, F. Urban growth patterns for Bucharest, Romania: Analysis of Landsat imagery. In Proceedings of the 34th EARSeL Symposium 2014, Prague, Czech Republic, 16–20 June 2014; pp. 65–70. [Google Scholar] [CrossRef]
  96. Zaharia, L.; Ioana-Toroimac, G.; Cocoş, O.; Ghiţă, F.A.; Mailat, E. Urbanization effects on the river systems in the Bucharest city region (Romania). Ecosyst. Health Sustain. 2016, 2, e01247. [Google Scholar] [CrossRef] [Green Version]
  97. Sebestyen, M. Urban image and national representation: Bucharest in the 19th and the beginning of the 20th century. In Planning Capital Cities: Belgrade, Bucharest, Sofia; Doytchinov, G., Dukić, A., Ioniță, C., Eds.; Verlag der Technischen Universität Graz: Graz, Austria, 2015; pp. 44–62. ISBN 973851253986. [Google Scholar]
  98. Cina’, G. Bucharest: From Village to Metropolis: Urban Identity and New Trends; Capitel: Bucharest, Romania, 2010; ISBN 9789731885001. [Google Scholar]
  99. Stroe, M. Bucharest’s urban planning instruments during the communist regime: Systematization sketches, plans, projects and interventions. In Planning Capital Cities: Belgrade, Bucharest, Sofia; Doytchinov, G., Dukić, A., Ioniță, C., Eds.; Verlag der Tech-nischen Universität Graz: Graz, Austria, 2015; pp. 116–140. ISBN 9783851253986. [Google Scholar]
  100. UIA. I.U. of A; Bucharest: Kristiansand, Norway, 2000. [Google Scholar]
  101. Otoiu, D. National (ist) Ideology and Urban Planning: Building the Victory of Socialism in Bucharest, Romania. In Art and Politics: Case-Studies from Eastern Europe; Levandauskas, V., Ed.; Vytautas Magnus University, Art Institute: Kaunas, Lithuania, 2007; pp. 119–129. [Google Scholar]
  102. Alexandru, M. Urban planning through major planning documents after 1999: Urban centrality between vision and reality. In Planning Capital Cities: Belgrade, Bucharest, Sofia; Doytchinov, G., Dukić, A., Ioniță, C., Eds.; Verlag der Technischen Universität Graz: Graz, Austria, 2015; pp. 234–248. ISBN 9783851253986. [Google Scholar]
  103. Ioana, T. Improving the Quality of Life in the Neighbourhood. COMMUNITY Centers as Means of Rehabilitation; New Europe College—Institute for Advanced Studies: Paris, France, 2004; pp. 309–346. [Google Scholar]
  104. Stoica, I.-V.; Vîrghileanu, M.; Zamfir, D.; Mihai, B.-A.; Săvulescu, I. Comparative Assessment of the Built-Up Area Expansion Based on Corine Land Cover and Landsat Datasets: A Case Study of a Post-Socialist City. Remote Sens. 2020, 12, 2137. [Google Scholar] [CrossRef]
  105. Bański, J. The consequences of changes of ownership for agricultural land use in Central European countries following the collapse of the Eastern Bloc. Land Use Policy 2017, 66, 120–130. [Google Scholar] [CrossRef]
  106. Grigorescu, I.; Mitrică, B.; Mocanu, I.; Ticană, N. Urban Sprawl and Residential Development in The Romanian Metropolitan Areas; Institute of Geography: Bucharest, Romania, 2012; pp. 43–59. [Google Scholar]
  107. Lascu, N. Legislation and Urban Development. Bucharest 1831–1952. (Legislaţie si Dezvoltarea Urbană. Bucureşti 1831–1952); Ion Mincu, University of Architecture: Bucharest, Romania, 1997. [Google Scholar]
  108. Udrea, A. The first urban plans of Bucharest in the rise of the 20th century. In Planning Capital Cities: Belgrade, Bucharest, Sofia; Doytchinov, G., Dukić, A., Ioniță, C., Eds.; Verlag der Technischen Universität Graz: Graz, Austria, 2015; pp. 62–80. ISBN 9783851253986. [Google Scholar]
Figure 1. Location map and overview of the city of Bucharest. Basemap: Natural colour composite of Sentinel-2 MSI from 14 August 2018 (ESA Copernicus).
Figure 1. Location map and overview of the city of Bucharest. Basemap: Natural colour composite of Sentinel-2 MSI from 14 August 2018 (ESA Copernicus).
Remotesensing 13 02323 g001
Figure 2. Flow chart of data processing for urban landscape change analysis.
Figure 2. Flow chart of data processing for urban landscape change analysis.
Remotesensing 13 02323 g002
Figure 3. LCLU maps of the City of Bucharest from 1968, with CORONA KH-4B orthophoto in background (left) and from 1989 with SPOT 1 panchromatic in background (right).
Figure 3. LCLU maps of the City of Bucharest from 1968, with CORONA KH-4B orthophoto in background (left) and from 1989 with SPOT 1 panchromatic in background (right).
Remotesensing 13 02323 g003
Figure 4. LCLU map of the City of Bucharest in 2018, based on Sentinel-1 and Sentinel-2 data.
Figure 4. LCLU map of the City of Bucharest in 2018, based on Sentinel-1 and Sentinel-2 data.
Remotesensing 13 02323 g004
Figure 5. Maps of the built-up area change 1968–1989 and 1989–2018 in the City of Bucharest.
Figure 5. Maps of the built-up area change 1968–1989 and 1989–2018 in the City of Bucharest.
Remotesensing 13 02323 g005
Figure 6. The percentage of LCLU classes in the city of Bucharest for each reference year.
Figure 6. The percentage of LCLU classes in the city of Bucharest for each reference year.
Remotesensing 13 02323 g006
Figure 7. Statistical approach representing the percentages: (a) of the 1968 LCLU features to the composition of the 1989 LCLU classes; (b) of the 1989 LCLU features to the composition of 2018 LULC classes.
Figure 7. Statistical approach representing the percentages: (a) of the 1968 LCLU features to the composition of the 1989 LCLU classes; (b) of the 1989 LCLU features to the composition of 2018 LULC classes.
Remotesensing 13 02323 g007
Figure 8. The Palace of Parliament area in the downtown of the city of Bucharest. Percent rates of the LCLU change (1968–2016) and qualitative features.
Figure 8. The Palace of Parliament area in the downtown of the city of Bucharest. Percent rates of the LCLU change (1968–2016) and qualitative features.
Remotesensing 13 02323 g008
Figure 9. Drumul Taberei district to the western edge of built-up area, an example of rural structures and agricultural land, replaced by a totally new urban district.
Figure 9. Drumul Taberei district to the western edge of built-up area, an example of rural structures and agricultural land, replaced by a totally new urban district.
Remotesensing 13 02323 g009
Figure 10. LCLU change in Băneasa zone, with increasing rate to the north of the Băneasa International Airport.
Figure 10. LCLU change in Băneasa zone, with increasing rate to the north of the Băneasa International Airport.
Remotesensing 13 02323 g010
Table 1. Remote sensing imagery used in the urban landscape change analysis of Bucharest City.
Table 1. Remote sensing imagery used in the urban landscape change analysis of Bucharest City.
DatasetsDateData SourceSpatial ResolutionSpectral ResolutionProcessing Level
CORONA KH-4B5 May 1968DECLASS 1 from USGS Earth Explorer2 mPanchromatic imageScanned film frames
SPOT 18 July 1989CNES10 mPanchromatic
image
Level 2A
20 m Multispectral
Image, 3 bands
Sentinel-1A SAR17 August 2018Copernicus ESA Sentinel Scientific Hub10 mC-bandIW-SLC, dual polarisation VV-VH
Sentinel-2A MSI14 August 2018Copernicus ESA Sentinel Scientific Hub10 m, 20 mMultispectral 13 bandsL1C processing level
Table 2. Urban LCLU key classes description for the city of Bucharest.
Table 2. Urban LCLU key classes description for the city of Bucharest.
Class NameDescriptionObservations
Block of flatsCompact areas with blocks of flats of 4–12 floors in big ensembles (1960s and 1970c) and linear block of flats pattern along boulevards (mainly 1980s).Building grouped with very narrow green areas, parking places, and playgrounds around with a typical pattern of planned areas.
Houses and other small buildingsHouses and gardens with one to two floor surrounded by gardens, in different patterns, from non-organized to rectangular. Usually small surface building, with gardens, remained from the traditional urban fabric (compact in the centre and more disperse to the periphery)—specific to the old urban landscape before the communist era.
Retail and industrial areasIndustrial building with big surfaces of halls and production facilities.
Commercial or retail areas with large buildings surrounded by parking places.
Difficult to separate industrial and commercial facilities with semi-automatic image supervised classification. Interpretation was the only solution to extract them from 1968 and 1989 images. Also includes the Băneasa international airport area and transportation facilities along railways.
StreetsStreet network including the entire paved/asphalt covered streets and boulevard/avenues of the city.In the 1968 image, the interpretation of these features made it possible to extract paved and non-paved streets. Semi-automatic approaches focus mainly on paved/ asphalt covered street areas (polygons).
Bare soil and construction sitesBarren land on the place of former demolished houses and industrial facilities, around railways and motorways as well as on the place of some former agricultural land/gardens at the periphery.Still remained around the Palace of Parliament (downtown) for a long time after 1984, largely extended on industrial areas before the conversion of land use to residential and commercial.
High vegetationForested grounds to the city periphery (ex. Băneasa, northern from airport) on large parks with compact configuration and linear patterns along alleys.Typical forested grounds from the regional forested landscape and adapted to the urban parks landscape in different periods.
Low vegetationPastures and other grass-covered areas, together with some agricultural land at the periphery.It corresponds to parks and natural reserves and also to open air stadium and sport facilities, as well as to some of the largest gardens around recently built leisure facilities around residential villa ensembles (after 2001).
Water bodiesLakes and rivers in built-up areasAnthropogenic lakes along Colentina and Dâmbovița rivers and the hydrotechnical works along Dâmbovița River floodplain.
Table 3. Classification confusion matrix for 2018 S1-S2 image stack.
Table 3. Classification confusion matrix for 2018 S1-S2 image stack.
Ground Truth ( %)
ClassArable LandWater BodiesHigh VegetationLow VegetationStreetsBuilt-Up HousesBuilt-Up Block of FlatsRetail Area
Arable land99.980.000.000.000.000.070.000.15
Water bodies0.0099.780.000.000.480.000.580.00
High vegetation0.000.0099.750.460.000.000.000.00
Low vegetation0.000.000.0199.480.000.080.000.00
Street0.000.060.000.0066.592.4412.7910.97
Built-up house0.010.000.000.0011.0695.363.1736.95
Built-up block of flats0.000.160.250.0321.751.6582.350.03
Retail area0.000.000.000.000.120.400.0051.59
Total100100100100100100100100
Overall accuracy: 94.0686%
Kappa coefficient: 0.9255
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Nistor, C.; Vîrghileanu, M.; Cârlan, I.; Mihai, B.-A.; Toma, L.; Olariu, B. Remote Sensing-Based Analysis of Urban Landscape Change in the City of Bucharest, Romania. Remote Sens. 2021, 13, 2323. https://doi.org/10.3390/rs13122323

AMA Style

Nistor C, Vîrghileanu M, Cârlan I, Mihai B-A, Toma L, Olariu B. Remote Sensing-Based Analysis of Urban Landscape Change in the City of Bucharest, Romania. Remote Sensing. 2021; 13(12):2323. https://doi.org/10.3390/rs13122323

Chicago/Turabian Style

Nistor, Constantin, Marina Vîrghileanu, Irina Cârlan, Bogdan-Andrei Mihai, Liviu Toma, and Bogdan Olariu. 2021. "Remote Sensing-Based Analysis of Urban Landscape Change in the City of Bucharest, Romania" Remote Sensing 13, no. 12: 2323. https://doi.org/10.3390/rs13122323

APA Style

Nistor, C., Vîrghileanu, M., Cârlan, I., Mihai, B. -A., Toma, L., & Olariu, B. (2021). Remote Sensing-Based Analysis of Urban Landscape Change in the City of Bucharest, Romania. Remote Sensing, 13(12), 2323. https://doi.org/10.3390/rs13122323

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop