Changes in the Building Stock of Da Nang between 2015 and 2017

This descriptor introduces a novel dataset, which contains the number and types of buildings in the city of Da Nang in Central Vietnam. The buildings were classified into nine distinct types and initially extracted from a satellite image of the year 2015. Secondly, changes were identified based on a visual interpretation of an image of the year 2017, so that new buildings, demolished buildings and building upgrades can be quantitatively analyzed. The data was aggregated by administrative wards and a hexagonal grid with a diameter of 250 m to protect personal rights and to avoid the misuse of a single building’s information. The dataset shows an increase of 19,391 buildings between October 2015 and August 2017, with a variety of interesting spatial patterns. The center of the city is mostly dominated by building changes and upgrades, while most of the new buildings were constructed within a distance of five to six kilometers from the city center. Dataset: https://doi.org/10.5281/zenodo.3757710 Dataset License: CC-BY 4.0

Rapidly urbanizing cities require reliable and up-to-date information about their inhabitants to ensure the provision of food and energy, but also for the planning of logistics, public services and infrastructure. However, the collection of this information by surveys or an official census is time-consuming and expensive. As an alternative, methods of earth observation can be employed to extract the number of buildings as an indicator of the size and spatial distribution of the urban population [1,2]. The presented dataset provides building types and numbers for the city of Da Nang for the years 2015 and 2017. The initial number of building types was retrieved from a very high-resolution (VHR) Pléiades satellite image using an object-based image analysis (OBIA) approach, as was already described in two previous studies by the authors [3,4]. Results from these studies were then refined by the approach presented in this data descriptor to extract the changes between 2015 and 2017 based on visual inspection and the identification of changes, classified into five different categories: no change, new building, changed building, upgraded building and demolished building. The building types and changes were aggregated in two datasets: administrative boundaries and hexagons of 250 m in width. These two datasets are the first publicly available and most up to date assessments of the building stock of Da Nang. Such information is valuable for urban planners to estimate the future growth of the city, the demand of the population (food, energy, or water) and the persons potentially affected by rising sea levels under climate change.

Methods
Several steps were required to identify the changes in the urban area of Da Nang, including the collection of reference data, the retrieval of built-up areas and building types from satellite images and the visual identification of changes over the investigated period. These steps are summarized in Figure 1. The building types and changes were aggregated in two datasets: administrative boundaries and hexagons of 250 m in width. These two datasets are the first publicly available and most up to date assessments of the building stock of Da Nang. Such information is valuable for urban planners to estimate the future growth of the city, the demand of the population (food, energy, or water) and the persons potentially affected by rising sea levels under climate change.

Methods
Several steps were required to identify the changes in the urban area of Da Nang, including the collection of reference data, the retrieval of built-up areas and building types from satellite images and the visual identification of changes over the investigated period. These steps are summarized in Figure 1.

Study Area
The study area covers the urban and suburban parts of the city of Da Nang, located in Central Vietnam. It has a population of 1.125 million inhabitants, with an annual growth rate around 2.75% according to the estimates of the United Nations [5]. Because of its harbor, which is the fourth largest in the country, the entirety of Central Vietnam has become economically important and has experienced a considerable population increase since 2007, mostly due to migration from rural areas [6]. The administrative area covered by the dataset is divided into seven districts and 50 wards ( Figure 2, black lines), of which six districts belong to the city of Da Nang. Because the acquisition of VHR data is limited by financial and technical constraints, an area of interest (AOI) of 225 square kilometers was defined for the analysis (Figure 2, red line), which covers the urban and suburban settlement areas of the city and large parts of its rural surroundings. Because of the seamless transition of the urban area at the southern border of the study area, parts of the Điện Bàn district of the adjacent Quảng Nam province are also covered by the AOI and thus have been included in the presented dataset.

Study Area
The study area covers the urban and suburban parts of the city of Da Nang, located in Central Vietnam. It has a population of 1.125 million inhabitants, with an annual growth rate around 2.75% according to the estimates of the United Nations [5]. Because of its harbor, which is the fourth largest in the country, the entirety of Central Vietnam has become economically important and has experienced a considerable population increase since 2007, mostly due to migration from rural areas [6]. The administrative area covered by the dataset is divided into seven districts and 50 wards ( Figure 2, black lines), of which six districts belong to the city of Da Nang. Because the acquisition of VHR data is limited by financial and technical constraints, an area of interest (AOI) of 225 square kilometers was defined for the analysis (Figure 2, red line), which covers the urban and suburban settlement areas of the city and large parts of its rural surroundings. Because of the seamless transition of the urban area at the southern border of the study area, parts of the Điện Bàn district of the adjacent Quảng Nam province are also covered by the AOI and thus have been included in the presented dataset.

Satellite Data
Two image products from the Pléiades satellite constellation acquired on 24. 10.2015 and 13.08.2017 were used in this study. For this period, the United Nations estimated a growth of 3.5% for the city of Da Nang [5]. The sensor acquires four multispectral bands (blue, green, red, and infrared) at 2.8 m and one panchromatic band at 0.7 m spatial resolution (resampled to 50 cm by the data provider). Both products were pan-sharpened to a final spatial resolution of 50 cm and radiometrically and geometrically calibrated as described by Warth et al. [4]. The two images were acquired from slightly different off-nadir angles, which resulted in the geometric shift of high buildings towards the sensor. Furthermore, the different acquisition dates within the year lead to differences in illumination and the shadows that are cast from large objects. Lastly, there are phenologic differences between October and August, so urban greenspaces or trees can produce unequal patterns in both images. All these slight distortions can make the visual identification of changes challenging, especially in densely built-up areas. Figure 3 gives an example on the differences between both images for a small selected area in the An Hải Bắc ward of the Sơn Trà district in the eastern part of Da Nang. Although the building in the center has not changed between both image acquisitions, the red color of the roof looks different because of sunlight illumination. Moreover, the shadows at the northern edges of tall buildings are only present in the image from October 2015. These seasonal differences were reduced during the radiometric calibration, but could not entirely be removed. Accordingly, these radiometric differences can potentially increase the chance of falsely assigned changes at the building level during the latter visual interpretation of changes at the building level. The figure also demonstrates that actual changes in the morphologic structure of the city, such as the creation of paved roads (in the northern and western part of the image), or the construction of new buildings (two bright lines in the center north) or other infrastructure (blue sports court in the center) can clearly be identified by visual comparison. However, slight differences had to be kept in mind by the interpreters in order to avoid errors of commission during the identification of changes between 2015 and 2017.

Satellite Data
Two image products from the Pléiades satellite constellation acquired on 24. 10.2015 and 13.08.2017 were used in this study. For this period, the United Nations estimated a growth of 3.5% for the city of Da Nang [5]. The sensor acquires four multispectral bands (blue, green, red, and infrared) at 2.8 m and one panchromatic band at 0.7 m spatial resolution (resampled to 50 cm by the data provider). Both products were pan-sharpened to a final spatial resolution of 50 cm and radiometrically and geometrically calibrated as described by Warth et al. [4]. The two images were acquired from slightly different off-nadir angles, which resulted in the geometric shift of high buildings towards the sensor. Furthermore, the different acquisition dates within the year lead to differences in illumination and the shadows that are cast from large objects. Lastly, there are phenologic differences between October and August, so urban greenspaces or trees can produce unequal patterns in both images. All these slight distortions can make the visual identification of changes challenging, especially in densely built-up areas. Figure 3 gives an example on the differences between both images for a small selected area in the An Hải Bắc ward of the Sơn Trà district in the eastern part of Da Nang. Although the building in the center has not changed between both image acquisitions, the red color of the roof looks different because of sunlight illumination. Moreover, the shadows at the northern edges of tall buildings are only present in the image from October 2015. These seasonal differences were reduced during the radiometric calibration, but could not entirely be removed. Accordingly, these radiometric differences can potentially increase the chance of falsely assigned changes at the building level during the latter visual interpretation of changes at the building level. The figure also demonstrates that actual changes in the morphologic structure of the city, such as the creation of paved roads (in the northern and western part of the image), or the construction of new buildings (two bright lines in the center north) or other infrastructure (blue sports court in the center) can clearly be identified by visual comparison. As a qualitative historic reference, an image acquired on 15.08.1977 by the Keyhole-9 (KH-9) mission "Hexagon" [7] was used in this study. Its spatial resolution of 7 m allowed us to identify the historic city center (see Section 3.2) of Da Nang.

Reference Data
To gain a better understanding of the city's morphology and the predominant building types, three field campaigns have been carried out in 2015 and 2016. Using mobile handheld tablet devices, a total number of 803 geocoded questionnaires on single buildings were collected, including information on building size, shape and height, its usage, the color and materials of the roof, photographs and further parameters regarding the structure of the buildings. These datasets were collected in March 2015 (n = 233), March 2016 (n = 465) and December 2016 (n = 105). Based on the collected data, building archetypes were defined in close collaboration with local experts. The spatial distribution of these buildings and the assigned building type are shown in Figure 1. Typically, in Vietnamese cities, shop houses with a longish shape and three to five floors of mixed use are the most common in the city. They can be refined into sub-classes (such as traditional, new, commercial, rowhouse), as initially suggested by Moon et al. (2009) [8], but only with additional knowledge about the degree of commercial use and construction materials, which cannot be retrieved from remote sensing data for the entire city.

Building Typology
The building typology was developed based on an existing study by Downes et al. (2016) for Ho Chi Minh City [9]. It has already been applied to estimate the number of different fractions of household waste [3] and was successfully transferred to the Rwandan capital, Kigali, to demonstrate the spatial patterns of urbanization between 2008 and 2015 [10]. The criteria for the building types are based on height, size, shape and roof material to allow an assignment of building types for the entire urban area of Da Nang, solely based on satellite image information. Table 1 lists the nine building types, together with a short description of their identification from remote sensing images. As a qualitative historic reference, an image acquired on 15.08.1977 by the Keyhole-9 (KH-9) mission "Hexagon" [7] was used in this study. Its spatial resolution of 7 m allowed us to identify the historic city center (see Section 3.2) of Da Nang.

Reference Data
To gain a better understanding of the city's morphology and the predominant building types, three field campaigns have been carried out in 2015 and 2016. Using mobile handheld tablet devices, a total number of 803 geocoded questionnaires on single buildings were collected, including information on building size, shape and height, its usage, the color and materials of the roof, photographs and further parameters regarding the structure of the buildings. These datasets were collected in March 2015 (n = 233), March 2016 (n = 465) and December 2016 (n = 105). Based on the collected data, building archetypes were defined in close collaboration with local experts. The spatial distribution of these buildings and the assigned building type are shown in Figure 1. Typically, in Vietnamese cities, shop houses with a longish shape and three to five floors of mixed use are the most common in the city. They can be refined into sub-classes (such as traditional, new, commercial, rowhouse), as initially suggested by Moon et al. (2009) [8], but only with additional knowledge about the degree of commercial use and construction materials, which cannot be retrieved from remote sensing data for the entire city.

Building Typology
The building typology was developed based on an existing study by Downes et al. (2016) for Ho Chi Minh City [9]. It has already been applied to estimate the number of different fractions of household waste [3] and was successfully transferred to the Rwandan capital, Kigali, to demonstrate the spatial patterns of urbanization between 2008 and 2015 [10]. The criteria for the building types are based on height, size, shape and roof material to allow an assignment of building types for the entire urban area of Da Nang, solely based on satellite image information. Table 1 lists the nine building types, together with a short description of their identification from remote sensing images. Detached housing, with a mix of residential and commercial use. Low-rise with one to two floors. Comprised of wood, brickwork, and reinforced concrete, tin roof.
Often located along small alleyways or in peri-urban locations.
Data 2020, 5, x FOR PEER REVIEW 5 of 18

Typical building type in Da
Nang. Detached/semidetached/terraced shophouse. It is a two to five story urban building, which allows for a shop or other public activity at the street level, with residential accommodation on the upper floors.
Data 2020, 5, x FOR PEER REVIEW 5 of 18 Detached housing, with a mix of residential and commercial use. Low-rise with one to two floors. Comprised of wood, brickwork, and reinforced concrete, tin roof. Often located along small alleyways or in periurban locations.  Mostly a single-family detached dwelling, sometimes two to three attached multifamily buildings, low rise, two to four floors, built from brickwork or concrete, located in newly developed urban areas.    Detached housing, with a mix of residential and commercial use. Low-rise with one to two floors. Comprised of wood, brickwork, and reinforced concrete, tin roof. Often located along small alleyways or in periurban locations.  Mostly a single-family detached dwelling, sometimes two to three attached multifamily buildings, low rise, two to four floors, built from brickwork or concrete, located in newly developed urban areas.   Mostly a single-family detached dwelling, sometimes two to three attached multifamily buildings, low rise, two to four floors, built from brickwork or concrete, located in newly developed urban areas.
Data 2020, 5, x FOR PEER REVIEW 5 of 18 Detached housing, with a mix of residential and commercial use. Low-rise with one to two floors. Comprised of wood, brickwork, and reinforced concrete, tin roof. Often located along small alleyways or in periurban locations. Mostly a single-family detached dwelling, sometimes two to three attached multifamily buildings, low rise, two to four floors, built from brickwork or concrete, located in newly developed urban areas.    Detached housing, with a mix of residential and commercial use. Low-rise with one to two floors. Comprised of wood, brickwork, and reinforced concrete, tin roof. Often located along small alleyways or in periurban locations. Mostly a single-family detached dwelling, sometimes two to three attached multifamily buildings, low rise, two to four floors, built from brickwork or concrete, located in newly developed urban areas.    Detached housing, with a mix of residential and commercial use. Low-rise with one to two floors. Comprised of wood, brickwork, and reinforced concrete, tin roof. Often located along small alleyways or in periurban locations. Mostly a single-family detached dwelling, sometimes two to three attached multifamily buildings, low rise, two to four floors, built from brickwork or concrete, located in newly developed urban areas.

External Data
To add a spatial relation to the identified changes, administrative boundaries were retrieved from the database of global administrative areas (GADM, [11]) which divides the study area into eight districts and 51 wards, of which 50 belong to Da Nang and one belongs to the adjacent Quảng Nam province. We use the wards as an administrative division of the city into homogenous units of analysis.
As a measure which is independent from political divisions, hexagonal cells with a width of 250 meters were computed. Compared to squared measures, hexagons are more suitable to represent nearest neighbor areas, reduce the sampling bias at the edges and are reportedly easier to interpret by the human eye [12]. Both measures allow for a detailed and precise analyses of the city without violating the personal property rights of the residents or house owners.

Data Processing and Change Detection
The delineation of built-up areas and the extraction of building heights from the Pléiades image of 2015 was conducted using an object-based image analysis approach, as described by Warth et al. [4]. As described in that study, built-up areas were identified using a multiresolution segmentation. This segmentation was based on the four spectral bands of the acquisition, their three principal components, a Canny edge detection algorithm [13], and a digital surface model retrieved from the tri-stereoscopic information of the Pléiades product [14]. The image objects were then classified into "built-up" and "non-built-up" using rule-based thresholds applied to the derived object features, resulting in an accuracy of 88.4%. All built-up areas were divided into single buildings using spatial intersection with cadastral boundaries and visual identification, as described by Vetter-Gindele et al. (2019) [3]. The geometric accuracy of single buildings was assessed by the manual digitization of 50 randomly selected buildings, whose areas were then compared to the automatically derived building footprints. The comparison resulted in an overall agreement of 87.2%. While accuracies of over 95% were achieved for buildings of average sizes (50-160 m²), larger underestimations were observed for very small (<50 m 2 ) and very large (>200 m 2 ) buildings, which only make up a small fraction of the buildings in Da Nang.

External Data
To add a spatial relation to the identified changes, administrative boundaries were retrieved from the database of global administrative areas (GADM, [11]) which divides the study area into eight districts and 51 wards, of which 50 belong to Da Nang and one belongs to the adjacent Quảng Nam province. We use the wards as an administrative division of the city into homogenous units of analysis.
As a measure which is independent from political divisions, hexagonal cells with a width of 250 meters were computed. Compared to squared measures, hexagons are more suitable to represent nearest neighbor areas, reduce the sampling bias at the edges and are reportedly easier to interpret by the human eye [12]. Both measures allow for a detailed and precise analyses of the city without violating the personal property rights of the residents or house owners.

Data Processing and Change Detection
The delineation of built-up areas and the extraction of building heights from the Pléiades image of 2015 was conducted using an object-based image analysis approach, as described by Warth et al. [4]. As described in that study, built-up areas were identified using a multiresolution segmentation. This segmentation was based on the four spectral bands of the acquisition, their three principal components, a Canny edge detection algorithm [13], and a digital surface model retrieved from the tri-stereoscopic information of the Pléiades product [14]. The image objects were then classified into "built-up" and "non-built-up" using rule-based thresholds applied to the derived object features, resulting in an accuracy of 88.4%. All built-up areas were divided into single buildings using spatial intersection with cadastral boundaries and visual identification, as described by Vetter-Gindele et al. (2019) [3]. The geometric accuracy of single buildings was assessed by the manual digitization of 50 randomly selected buildings, whose areas were then compared to the automatically derived building footprints. The comparison resulted in an overall agreement of 87.2%. While accuracies of over 95% were achieved for buildings of average sizes (50-160 m²), larger underestimations were observed for very small (<50 m 2 ) and very large (>200 m 2 ) buildings, which only make up a small fraction of the buildings in Da Nang.

External Data
To add a spatial relation to the identified changes, administrative boundaries were retrieved from the database of global administrative areas (GADM, [11]) which divides the study area into eight districts and 51 wards, of which 50 belong to Da Nang and one belongs to the adjacent Quảng Nam province. We use the wards as an administrative division of the city into homogenous units of analysis.
As a measure which is independent from political divisions, hexagonal cells with a width of 250 meters were computed. Compared to squared measures, hexagons are more suitable to represent nearest neighbor areas, reduce the sampling bias at the edges and are reportedly easier to interpret by the human eye [12]. Both measures allow for a detailed and precise analyses of the city without violating the personal property rights of the residents or house owners.

Data Processing and Change Detection
The delineation of built-up areas and the extraction of building heights from the Pléiades image of 2015 was conducted using an object-based image analysis approach, as described by Warth et al. [4]. As described in that study, built-up areas were identified using a multiresolution segmentation. This segmentation was based on the four spectral bands of the acquisition, their three principal components, a Canny edge detection algorithm [13], and a digital surface model retrieved from the tri-stereoscopic information of the Pléiades product [14]. The image objects were then classified into "built-up" and "non-built-up" using rule-based thresholds applied to the derived object features, resulting in an accuracy of 88.4%. All built-up areas were divided into single buildings using spatial intersection with cadastral boundaries and visual identification, as described by Vetter-Gindele et al. (2019) [3]. The geometric accuracy of single buildings was assessed by the manual digitization of 50 randomly selected buildings, whose areas were then compared to the automatically derived building footprints. The comparison resulted in an overall agreement of 87.2%. While accuracies of over 95% were achieved for buildings of average sizes (50-160 m²), larger underestimations were observed for very small (<50 m 2 ) and very large (>200 m 2 ) buildings, which only make up a small fraction of the buildings in Da Nang.

External Data
To add a spatial relation to the identified changes, administrative boundaries were retrieved from the database of global administrative areas (GADM, [11]) which divides the study area into eight districts and 51 wards, of which 50 belong to Da Nang and one belongs to the adjacent Quảng Nam province. We use the wards as an administrative division of the city into homogenous units of analysis.
As a measure which is independent from political divisions, hexagonal cells with a width of 250 m were computed. Compared to squared measures, hexagons are more suitable to represent nearest neighbor areas, reduce the sampling bias at the edges and are reportedly easier to interpret by the human eye [12]. Both measures allow for a detailed and precise analyses of the city without violating the personal property rights of the residents or house owners.

Data Processing and Change Detection
The delineation of built-up areas and the extraction of building heights from the Pléiades image of 2015 was conducted using an object-based image analysis approach, as described by Warth et al. [4]. As described in that study, built-up areas were identified using a multiresolution segmentation. This segmentation was based on the four spectral bands of the acquisition, their three principal components, a Canny edge detection algorithm [13], and a digital surface model retrieved from the tri-stereoscopic information of the Pléiades product [14]. The image objects were then classified into "built-up" and "non-built-up" using rule-based thresholds applied to the derived object features, resulting in an accuracy of 88.4%. All built-up areas were divided into single buildings using spatial intersection with cadastral boundaries and visual identification, as described by Vetter-Gindele et al. (2019) [3]. The geometric accuracy of single buildings was assessed by the manual digitization of 50 randomly selected Data 2020, 5, 42 7 of 18 buildings, whose areas were then compared to the automatically derived building footprints. The comparison resulted in an overall agreement of 87.2%. While accuracies of over 95% were achieved for buildings of average sizes (50-160 m 2 ), larger underestimations were observed for very small (<50 m 2 ) and very large (>200 m 2 ) buildings, which only make up a small fraction of the buildings in Da Nang. In the end, the derived footprints resulted in building areas which were found to be reliable enough to serve as one of the variables for the classification of building types in the next step.
Accordingly, building types were assigned based on the criteria given in Table 1 and manual refinement. This classification typology resulted in an overall number of 241,217 building features in the study area for the year 2015. To understand the quality of the building type classification, an accuracy assessment was conducted based on the 803 building references which were collected during the field campaign, as described in Section 2.3 [3]. The confusion matrix between the reference data and the classified building types is shown in Table A1 in the Appendix A and discussed in more detail in Section 3.5. An overall classification accuracy of 90.9% was achieved.
In the next step, a typology of changes was developed, consisting of five classes ( Table 2): no change (1), construction of a new building (2), changes in the building, e.g., roof color or size, while maintaining the building type according to Table 1 (3), upgrading of the building to a new building type, mostly of higher quality (4), and demolition of the building (5). All buildings identified for the year 2015 were then examined for changes by a visual comparison of the two images and one type of change was assigned to each feature. In the case of a newly constructed building, a new feature was manually digitized at this location. This resulted in a total number of 260,608 buildings for 2017, which indicates an increase of 19,391 for the investigated area.
The process of visual interpretation was conducted by independent persons, with additional quality control by another person to ensure the highest level of objectivity and to minimize errors of omission and commission. However, as demonstrated by Figure 3, the identification of changes can still be subject to the interpreter's perception of the two images and their ability to identify changes. This was especially challenging in densely built-up areas where buildings stand extremely close (Figure 4, left part). This required a more time-consuming inspection and an independent crosscheck in those areas. For the more regular parts of Da Nang (Figure 4, right part) changes were clearly visible and unambiguous.
The selection of a manual assignment of building types for 2017 and the derivation of the corresponding changes was considerably faster than repeating the semi-automated building classification as it was conducted for the image of 2015. At the same time, it allowed us to achieve higher accuracies in the final data product-for example, by manually correcting smaller misclassifications of the building types of 2015 while searching for changes between both images. In the end, the derived footprints resulted in building areas which were found to be reliable enough to serve as one of the variables for the classification of building types in the next step. Accordingly, building types were assigned based on the criteria given in Table 1 and manual refinement. This classification typology resulted in an overall number of 241,217 building features in the study area for the year 2015. To understand the quality of the building type classification, an accuracy assessment was conducted based on the 803 building references which were collected during the field campaign, as described in section 2.3. [3]. The confusion matrix between the reference data and the classified building types is shown in Table A1 in the appendix and discussed in more detail in section 3.5. An overall classification accuracy of 90.9% was achieved.
In the next step, a typology of changes was developed, consisting of five classes ( Table 2): no change (1), construction of a new building (2), changes in the building, e.g., roof color or size, while maintaining the building type according to Table 1 (3), upgrading of the building to a new building type, mostly of higher quality (4), and demolition of the building (5). All buildings identified for the year 2015 were then examined for changes by a visual comparison of the two images and one type of change was assigned to each feature. In the case of a newly constructed building, a new feature was manually digitized at this location. This resulted in a total number of 260,608 buildings for 2017, which indicates an increase of 19,391 for the investigated area.
The process of visual interpretation was conducted by independent persons, with additional quality control by another person to ensure the highest level of objectivity and to minimize errors of omission and commission. However, as demonstrated by Figure 3, the identification of changes can still be subject to the interpreter's perception of the two images and their ability to identify changes. This was especially challenging in densely built-up areas where buildings stand extremely close (Figure 4, left part). This required a more time-consuming inspection and an independent crosscheck in those areas. For the more regular parts of Da Nang (Figure 4, right part) changes were clearly visible and unambiguous.
The selection of a manual assignment of building types for 2017 and the derivation of the corresponding changes was considerably faster than repeating the semi-automated building classification as it was conducted for the image of 2015. At the same time, it allowed us to achieve higher accuracies in the final data product-for example, by manually correcting smaller misclassifications of the building types of 2015 while searching for changes between both images. In the end, the derived footprints resulted in building areas which were found to be reliable enough to serve as one of the variables for the classification of building types in the next step. Accordingly, building types were assigned based on the criteria given in Table 1 and manual refinement. This classification typology resulted in an overall number of 241,217 building features in the study area for the year 2015. To understand the quality of the building type classification, an accuracy assessment was conducted based on the 803 building references which were collected during the field campaign, as described in section 2.3. [3]. The confusion matrix between the reference data and the classified building types is shown in Table A1 in the appendix and discussed in more detail in section 3.5. An overall classification accuracy of 90.9% was achieved.
In the next step, a typology of changes was developed, consisting of five classes ( Table 2): no change (1), construction of a new building (2), changes in the building, e.g., roof color or size, while maintaining the building type according to Table 1 (3), upgrading of the building to a new building type, mostly of higher quality (4), and demolition of the building (5). All buildings identified for the year 2015 were then examined for changes by a visual comparison of the two images and one type of change was assigned to each feature. In the case of a newly constructed building, a new feature was manually digitized at this location. This resulted in a total number of 260,608 buildings for 2017, which indicates an increase of 19,391 for the investigated area.
The process of visual interpretation was conducted by independent persons, with additional quality control by another person to ensure the highest level of objectivity and to minimize errors of omission and commission. However, as demonstrated by Figure 3, the identification of changes can still be subject to the interpreter's perception of the two images and their ability to identify changes. This was especially challenging in densely built-up areas where buildings stand extremely close (Figure 4, left part). This required a more time-consuming inspection and an independent crosscheck in those areas. For the more regular parts of Da Nang (Figure 4, right part) changes were clearly visible and unambiguous.
The selection of a manual assignment of building types for 2017 and the derivation of the corresponding changes was considerably faster than repeating the semi-automated building classification as it was conducted for the image of 2015. At the same time, it allowed us to achieve higher accuracies in the final data product-for example, by manually correcting smaller misclassifications of the building types of 2015 while searching for changes between both images. In the end, the derived footprints resulted in building areas which were found to be reliable enough to serve as one of the variables for the classification of building types in the next step. Accordingly, building types were assigned based on the criteria given in Table 1 and manual refinement. This classification typology resulted in an overall number of 241,217 building features in the study area for the year 2015. To understand the quality of the building type classification, an accuracy assessment was conducted based on the 803 building references which were collected during the field campaign, as described in section 2.3. [3]. The confusion matrix between the reference data and the classified building types is shown in Table A1 in the appendix and discussed in more detail in section 3.5. An overall classification accuracy of 90.9% was achieved.
In the next step, a typology of changes was developed, consisting of five classes ( Table 2): no change (1), construction of a new building (2), changes in the building, e.g., roof color or size, while maintaining the building type according to Table 1 (3), upgrading of the building to a new building type, mostly of higher quality (4), and demolition of the building (5). All buildings identified for the year 2015 were then examined for changes by a visual comparison of the two images and one type of change was assigned to each feature. In the case of a newly constructed building, a new feature was manually digitized at this location. This resulted in a total number of 260,608 buildings for 2017, which indicates an increase of 19,391 for the investigated area.
The process of visual interpretation was conducted by independent persons, with additional quality control by another person to ensure the highest level of objectivity and to minimize errors of omission and commission. However, as demonstrated by Figure 3, the identification of changes can still be subject to the interpreter's perception of the two images and their ability to identify changes. This was especially challenging in densely built-up areas where buildings stand extremely close (Figure 4, left part). This required a more time-consuming inspection and an independent crosscheck in those areas. For the more regular parts of Da Nang (Figure 4, right part) changes were clearly visible and unambiguous.
The selection of a manual assignment of building types for 2017 and the derivation of the corresponding changes was considerably faster than repeating the semi-automated building classification as it was conducted for the image of 2015. At the same time, it allowed us to achieve higher accuracies in the final data product-for example, by manually correcting smaller misclassifications of the building types of 2015 while searching for changes between both images. In the end, the derived footprints resulted in building areas which were found to be reliable enough to serve as one of the variables for the classification of building types in the next step. Accordingly, building types were assigned based on the criteria given in Table 1 and manual refinement. This classification typology resulted in an overall number of 241,217 building features in the study area for the year 2015. To understand the quality of the building type classification, an accuracy assessment was conducted based on the 803 building references which were collected during the field campaign, as described in section 2.3. [3]. The confusion matrix between the reference data and the classified building types is shown in Table A1 in the appendix and discussed in more detail in section 3.5. An overall classification accuracy of 90.9% was achieved.
In the next step, a typology of changes was developed, consisting of five classes ( Table 2): no change (1), construction of a new building (2), changes in the building, e.g., roof color or size, while maintaining the building type according to Table 1 (3), upgrading of the building to a new building type, mostly of higher quality (4), and demolition of the building (5). All buildings identified for the year 2015 were then examined for changes by a visual comparison of the two images and one type of change was assigned to each feature. In the case of a newly constructed building, a new feature was manually digitized at this location. This resulted in a total number of 260,608 buildings for 2017, which indicates an increase of 19,391 for the investigated area.
The process of visual interpretation was conducted by independent persons, with additional quality control by another person to ensure the highest level of objectivity and to minimize errors of omission and commission. However, as demonstrated by Figure 3, the identification of changes can still be subject to the interpreter's perception of the two images and their ability to identify changes. This was especially challenging in densely built-up areas where buildings stand extremely close (Figure 4, left part). This required a more time-consuming inspection and an independent crosscheck in those areas. For the more regular parts of Da Nang (Figure 4, right part) changes were clearly visible and unambiguous.
The selection of a manual assignment of building types for 2017 and the derivation of the corresponding changes was considerably faster than repeating the semi-automated building classification as it was conducted for the image of 2015. At the same time, it allowed us to achieve higher accuracies in the final data product-for example, by manually correcting smaller misclassifications of the building types of 2015 while searching for changes between both images.  In the final step, the building types and their changes between 2015 and 2017 were summarized for each ward (Figure 4, black dashed lines) and hexagon (Figure 4, blue lines). As described in Table  3, each polygon contains attributes on the absolute number of buildings for both years, as well as the absolute and relative amount of change.

Datatype
The two datasets presented in this descriptor (please see the zip file in the Supplementary Materials) will be described in the next sections and the changes in the building stock for the urban areas of Da Nang city, aggregated by administrative wards and a systematic grid of hexagon of 250 meters in width, will be explained. Both are stored as vector datasets in the ESRI shapefile format, which is compatible with common Geographic Information Systems (GIS), such as QGIS or ArcMap. Each dataset consists of several separate files with the same name, but different extensions: a vector geometry (*.shp), an attribute table (*.dbf), a projection file (*.prj), and an index linking the attribute In the final step, the building types and their changes between 2015 and 2017 were summarized for each ward (Figure 4, black dashed lines) and hexagon (Figure 4, blue lines). As described in Table  3, each polygon contains attributes on the absolute number of buildings for both years, as well as the absolute and relative amount of change.

Datatype
The two datasets presented in this descriptor (please see the zip file in the Supplementary Materials) will be described in the next sections and the changes in the building stock for the urban areas of Da Nang city, aggregated by administrative wards and a systematic grid of hexagon of 250 meters in width, will be explained. Both are stored as vector datasets in the ESRI shapefile format, which is compatible with common Geographic Information Systems (GIS), such as QGIS or ArcMap. Each dataset consists of several separate files with the same name, but different extensions: a vector geometry (*.shp), an attribute table (*.dbf), a projection file (*.prj), and an index linking the attribute In the final step, the building types and their changes between 2015 and 2017 were summarized for each ward (Figure 4, black dashed lines) and hexagon (Figure 4, blue lines). As described in Table  3, each polygon contains attributes on the absolute number of buildings for both years, as well as the absolute and relative amount of change.

Datatype
The two datasets presented in this descriptor (please see the zip file in the Supplementary Materials) will be described in the next sections and the changes in the building stock for the urban areas of Da Nang city, aggregated by administrative wards and a systematic grid of hexagon of 250 meters in width, will be explained. Both are stored as vector datasets in the ESRI shapefile format, which is compatible with common Geographic Information Systems (GIS), such as QGIS or ArcMap. Each dataset consists of several separate files with the same name, but different extensions: a vector geometry (*.shp), an attribute table (*.dbf), a projection file (*.prj), and an index linking the attribute table to the geometries (*.shx). Both datasets are stored in the UTM coordinate system, zone 49N (EPSG: 32649). The encoding of the attribute table (*.cpg) is UTF-8 to support the correct display of  (Figure 4, blue lines). As described in Table  3, each polygon contains attributes on the absolute number of buildings for both years, as well as the absolute and relative amount of change.

Datatype
The two datasets presented in this descriptor (please see the zip file in the Supplementary Materials) will be described in the next sections and the changes in the building stock for the urban areas of Da Nang city, aggregated by administrative wards and a systematic grid of hexagon of 250 meters in width, will be explained. Both are stored as vector datasets in the ESRI shapefile format, which is compatible with common Geographic Information Systems (GIS), such as QGIS or ArcMap. Each dataset consists of several separate files with the same name, but different extensions: a vector geometry (*.shp), an attribute table (*.dbf), a projection file (*.prj), and an index linking the attribute table to the geometries (*.shx). Both datasets are stored in the UTM coordinate system, zone 49N (EPSG: 32649). The encoding of the attribute table (*.cpg) is UTF-8 to support the correct display of In the final step, the building types and their changes between 2015 and 2017 were summarized for each ward (Figure 4, black dashed lines) and hexagon (Figure 4, blue lines). As described in Table 3, each polygon contains attributes on the absolute number of buildings for both years, as well as the absolute and relative amount of change.  (Figure 4, blue lines). As described in Table  3, each polygon contains attributes on the absolute number of buildings for both years, as well as the absolute and relative amount of change.

Datatype
The two datasets presented in this descriptor (please see the zip file in the Supplementary Materials) will be described in the next sections and the changes in the building stock for the urban areas of Da Nang city, aggregated by administrative wards and a systematic grid of hexagon of 250 m in width, will be explained. Both are stored as vector datasets in the ESRI shapefile format, which is compatible with common Geographic Information Systems (GIS), such as QGIS or ArcMap. Each dataset consists of several separate files with the same name, but different extensions: a vector geometry (*.shp), an attribute table (*.dbf), a projection file (*.prj), and an index linking the attribute

Datatype
The two datasets presented in this descriptor (please see the zip file in the Supplementary Materials) will be described in the next sections and the changes in the building stock for the urban areas of Da Nang city, aggregated by administrative wards and a systematic grid of hexagon of 250 m in width, Data 2020, 5, 42 9 of 18 will be explained. Both are stored as vector datasets in the ESRI shapefile format, which is compatible with common Geographic Information Systems (GIS), such as QGIS or ArcMap. Each dataset consists of several separate files with the same name, but different extensions: a vector geometry (*.shp), an attribute table (*.dbf), a projection file (*.prj), and an index linking the attribute table to the geometries (*.shx). Both datasets are stored in the UTM coordinate system, zone 49N (EPSG: 32649). The encoding of the attribute table (*.cpg) is UTF-8 to support the correct display of Vietnamese characters. Table 3 summarizes the attributes contained by the administrative (adm) and hexagonal (hex) datasets. Both datasets contain the absolute number of changes (c1 to c5, as defined in Table 2), the absolute number of buildings per year and type (t2015_1 to t2017_9, as defined in Table 1), as well as the absolute (nc) number of changes per building type in each ward and hexagon. The relative number of changes (pc) denotes the share of the changes in one building type compared to the total number of changes in the reference area. For the administrative dataset, the results can also be summarized by district, using the respective attributes (dist_ID or district). As one suggestion for further spatial analyses, the hexagon dataset additionally contains an attribute named dist_mean, which defines the Euclidean distance of each hexagon from the historic city center in Hải Châu district. This will be discussed in Section 3.4. Accordingly, the dataset does not contain any data related to individuals collected during the field campaigns. The attributes are purely a result of the geospatial analyses described in the methodology section. This aggregation grants that no user can calculate the original building or plot geometries, nor their exact type or location.

Administrative Level
The dataset <DaNang_buildings_adm_UTM49N> contains the changes in the building stock aggregated by wards. These represent the official administrative units, so this dataset is best compatible with other socio-demographic data, such as the official Vietnam Population and Housing Census released by the General Statistics Office of Vietnam [15]. However, it must be kept in mind that, in particular, the wards at the edges are not entirely covered by the satellite images, so absolute numbers have to be interpreted with care. For this reason, the attribute table of the administrative dataset contains an additional attribute, "per_cover", which tells the percentage of each ward's area inside the AOI. It is illustrated in Figure A1 in the Appendix B, which highlights the ten wards which are only covered by 70% or less of the satellite images. Figure 5 shows the dominance of new construction activities over demolitions in the observed time between 2015 and 2017. Except for the already very densely built-up wards of An Hải Đông, Bình Hiên, and Phước Ninh, which are located near the city center, and the rural wards of Hòa Vang district in the western outskirts of the city, all parts of Da Nang have experienced a considerable increase in building numbers, mostly exceeding 500 new buildings per ward. In contrast, the number of demolitions is visibly lower throughout the city. Over 90% of all wards range between 10 and 200 demolished buildings. Of course, these numbers are strongly dependent on the size of a ward. For this reason, the attribute "area" is contained in the dataset to calculate the absolute changes per square kilometer (c1/area or c5/area), as shown in Figure 6, which gives a more representative picture for each ward. It highlights the increase in Hòa Minh in the north and the relative high demolitions in Hải Châu district, representing the historic city center with its extreme building density. Another method of calculating a spatially uncorrelated representation of changes is given by the hexagons described in Section 3.4.   Figure 7 demonstrates further spatial patterns revealed in the data: It shows that multifamily apartments largely increased between 2015 and 2017, but with spatially distinct patterns. While localtype apartments (building type 5 in Table 2), which are predominantly inhabited by local people, show indifferent patterns and even slight decreases in the outer wards, there is a tremendous increase in hotels (modern-type apartments, building type 6 in Table 2) along the coasts in the east (Ngũ Hành Sơn district) and in the north (Liên Chiểu district), which underlines Da Nang's ambitions to increase its touristic sector [16].   Figure 7 demonstrates further spatial patterns revealed in the data: It shows that multifamily apartments largely increased between 2015 and 2017, but with spatially distinct patterns. While localtype apartments (building type 5 in Table 2), which are predominantly inhabited by local people, show indifferent patterns and even slight decreases in the outer wards, there is a tremendous increase in hotels (modern-type apartments, building type 6 in Table 2) along the coasts in the east (Ngũ Hành Sơn district) and in the north (Liên Chiểu district), which underlines Da Nang's ambitions to increase its touristic sector [16].  Figure 7 demonstrates further spatial patterns revealed in the data: It shows that multifamily apartments largely increased between 2015 and 2017, but with spatially distinct patterns. While local-type apartments (building type 5 in Table 2), which are predominantly inhabited by local people, show indifferent patterns and even slight decreases in the outer wards, there is a tremendous increase in hotels (modern-type apartments, building type 6 in Table 2) along the coasts in the east (Ngũ Hành Sơn district) and in the north (Liên Chiểu district), which underlines Da Nang's ambitions to increase its touristic sector [16].
However, as already stated above, the size of the wards largely differs and ranges between 0.23 and 47 square kilometers, with an average ward size of 7.1 square kilometers. Accordingly, absolute building numbers and their changes are spatially biased by the ward size. For this reason, a second dataset is introduced, which provides a regular representation of the aggregated results. Data 2020, 5, x FOR PEER REVIEW 12 of 18 However, as already stated above, the size of the wards largely differs and ranges between 0.23 and 47 square kilometers, with an average ward size of 7.1 square kilometers. Accordingly, absolute building numbers and their changes are spatially biased by the ward size. For this reason, a second dataset is introduced, which provides a regular representation of the aggregated results.

Spatial Aggregation
The dataset <DaNang_buildings_hex250m_UTM49N.shp> aggregates the changes at a spatially regular grid, where each hexagon has a width of 250 m and encloses an area of 72 hectares. This allows us to identify patterns that are not related to administrative boundaries. As an example, Figure  8 shows the absolute differences in bungalows and villa-type buildings calculated by the following equations: t2017_3 -t2015_3 and t2017_4 -t2015_4. It shows that these building types and their changes are not drastically shaping the morphology of the city, as most of the hexagons have values close to zero. However, it also shows that most of the bungalows were demolished in the new urban development areas of Hòa Xuân (Cẩm Lệ district) and Hòa Quy (Ngũ Hành Sơn district) in the southwest of the city center, which is one of Da Nang's fastest urbanizing areas. Large masses of sand were piled to create new space for development at a confluence area of the Hán River. In turn, most of the new villa-type buildings are constructed along the coastal area in the east to increase the number of holiday residences.

Spatial Aggregation
The dataset <DaNang_buildings_hex250m_UTM49N.shp> aggregates the changes at a spatially regular grid, where each hexagon has a width of 250 m and encloses an area of 72 hectares. This allows us to identify patterns that are not related to administrative boundaries. As an example, Figure 8 shows the absolute differences in bungalows and villa-type buildings calculated by the following equations: t2017_3-t2015_3 and t2017_4-t2015_4. It shows that these building types and their changes are not drastically shaping the morphology of the city, as most of the hexagons have values close to zero. However, it also shows that most of the bungalows were demolished in the new urban development areas of Hòa Xuân (Cẩm Lệ district) and Hòa Quy (Ngũ Hành Sơn district) in the southwest of the city center, which is one of Da Nang's fastest urbanizing areas. Large masses of sand were piled to create new space for development at a confluence area of the Hán River. In turn, most of the new villa-type buildings are constructed along the coastal area in the east to increase the number of holiday residences. However, as already stated above, the size of the wards largely differs and ranges between 0.23 and 47 square kilometers, with an average ward size of 7.1 square kilometers. Accordingly, absolute building numbers and their changes are spatially biased by the ward size. For this reason, a second dataset is introduced, which provides a regular representation of the aggregated results.

Spatial Aggregation
The dataset <DaNang_buildings_hex250m_UTM49N.shp> aggregates the changes at a spatially regular grid, where each hexagon has a width of 250 m and encloses an area of 72 hectares. This allows us to identify patterns that are not related to administrative boundaries. As an example, Figure  8 shows the absolute differences in bungalows and villa-type buildings calculated by the following equations: t2017_3 -t2015_3 and t2017_4 -t2015_4. It shows that these building types and their changes are not drastically shaping the morphology of the city, as most of the hexagons have values close to zero. However, it also shows that most of the bungalows were demolished in the new urban development areas of Hòa Xuân (Cẩm Lệ district) and Hòa Quy (Ngũ Hành Sơn district) in the southwest of the city center, which is one of Da Nang's fastest urbanizing areas. Large masses of sand were piled to create new space for development at a confluence area of the Hán River. In turn, most of the new villa-type buildings are constructed along the coastal area in the east to increase the number of holiday residences.  The notable increase in buildings in Cam Le Cẩm Lệ is furthermore underlined by Figure 9, where two hotspots of new buildings (change type 2 as defined in Table 1) are visible: these are the wards of An Khê, Hòa Phát, and Hòa Minh in the north and the new construction areas in Hòa Xuân ward of Cẩm Lệ district enclosed by the Hàn River in the southwest. In particular, the increase in built-up areas in Cẩm Lệ was already confirmed in 2016 by Bachofer and Rau [17] and Linh and Chuong [18]. This is especially alarming because these new residential areas are expected to lie below the tideline in 2050 [19] (Figure A3). In turn, most of the building changes (within the same type, change type 3 as defined in Table 1) occurred in the districts Hải Châu and Ngũ Hành Sơn, around the historic center of Da Nang. These are the parts of the city which are already extremely dense ( Figure A2) and probably do not allow large increases in the building stock. Instead, additional parts or floors were added to the existing ones.
The notable increase in buildings in Cam Le Cẩm Lệ is furthermore underlined by Figure 9, where two hotspots of new buildings (change type 2 as defined in Table 1) are visible: these are the wards of An Khê, Hòa Phát, and Hòa Minh in the north and the new construction areas in Hòa Xuân ward of Cẩm Lệ district enclosed by the Hàn River in the southwest. In particular, the increase in built-up areas in Cẩm Lệ was already confirmed in 2016 by Bachofer and Rau [17] and Linh and Chuong [18]. This is especially alarming because these new residential areas are expected to lie below the tideline in 2050 [19] (Figure A3). In turn, most of the building changes (within the same type, change type 3 as defined in Table 1) occurred in the districts Hải Châu and Ngũ Hành Sơn, around the historic center of Da Nang. These are the parts of the city which are already extremely dense ( Figure A2) and probably do not allow large increases in the building stock. Instead, additional parts or floors were added to the existing ones.  [4]. However, they used hexagons with a diameter of 500 m, so both datasets were normalized by calculating the demolitions per square kilometer. As shown in Figure 10, the presented dataset not only offers finer spatial details, which make it more suitable to describe urban changes, but is also less prone to systematic overestimation, as it was reported by Warth et al. (2019) [4].  To furthermore demonstrate the value of the presented dataset, we compared the number of demolitions displayed by the hexagonal aggregation to the demolitions automatically identified by , who derived digital surface models of the city from 2015 and 2017 to detect demolitions based on the threshold of the height difference [4]. However, they used hexagons with a diameter of 500 m, so both datasets were normalized by calculating the demolitions per square kilometer. As shown in Figure 10, the presented dataset not only offers finer spatial details, which make it more suitable to describe urban changes, but is also less prone to systematic overestimation, as it was reported by Warth et al. (2019) [4].
The notable increase in buildings in Cam Le Cẩm Lệ is furthermore underlined by Figure 9, where two hotspots of new buildings (change type 2 as defined in Table 1) are visible: these are the wards of An Khê, Hòa Phát, and Hòa Minh in the north and the new construction areas in Hòa Xuân ward of Cẩm Lệ district enclosed by the Hàn River in the southwest. In particular, the increase in built-up areas in Cẩm Lệ was already confirmed in 2016 by Bachofer and Rau [17] and Linh and Chuong [18]. This is especially alarming because these new residential areas are expected to lie below the tideline in 2050 [19] (Figure A3). In turn, most of the building changes (within the same type, change type 3 as defined in Table 1) occurred in the districts Hải Châu and Ngũ Hành Sơn, around the historic center of Da Nang. These are the parts of the city which are already extremely dense ( Figure A2) and probably do not allow large increases in the building stock. Instead, additional parts or floors were added to the existing ones.  [4]. However, they used hexagons with a diameter of 500 m, so both datasets were normalized by calculating the demolitions per square kilometer. As shown in Figure 10, the presented dataset not only offers finer spatial details, which make it more suitable to describe urban changes, but is also less prone to systematic overestimation, as it was reported by Warth et al. (2019) [4].

Statistical Description
Lastly, we highlight the potential for statistical analyses of the prepared data, which can be conducted based on the attribute tables of the shapefiles. Figure 11 aggregates the changes at the

Statistical Description
Lastly, we highlight the potential for statistical analyses of the prepared data, which can be conducted based on the attribute tables of the shapefiles. Figure 11 aggregates the changes at the building level for seven districts covered by the satellite images (please note that Điện Bàn district was excluded because of very low coverage and small values). It highlights the enormous share of changes in the total building stock of the city. Figure 11. Changed and unchanged buildings for the districts in Da Nang. Figure 12 shows the proportion of new and demolished buildings based on the hexagon dataset, in relation to their average distance to the historic city center of Da Nang. Interestingly, both types of changes show a similar pattern, even if they are not displayed in absolute values. Accordingly, 50% of all new buildings and 43% of all demolitions between 2015 and 2017 occurred within a distance of three to six kilometers from the historic city center. Within a radius of two kilometers, 12.9% of new buildings were built, but 21.2% of demolitions occurred. This again underlines that the central residential areas in Da Nang are already very dense and do not allow much more growth. Instead, the amount of changed buildings (change type 3, as defined in Table 1) is visibly higher in the more central areas. Of all changed buildings in Da Nang, 77% were located within the first six kilometers.   Figure 12 shows the proportion of new and demolished buildings based on the hexagon dataset, in relation to their average distance to the historic city center of Da Nang. Interestingly, both types of changes show a similar pattern, even if they are not displayed in absolute values. Accordingly, 50% of all new buildings and 43% of all demolitions between 2015 and 2017 occurred within a distance of three to six kilometers from the historic city center. Within a radius of two kilometers, 12.9% of new buildings were built, but 21.2% of demolitions occurred. This again underlines that the central residential areas in Da Nang are already very dense and do not allow much more growth. Instead, the amount of changed buildings (change type 3, as defined in Table 1) is visibly higher in the more central areas. Of all changed buildings in Da Nang, 77% were located within the first six kilometers.

Statistical Description
Lastly, we highlight the potential for statistical analyses of the prepared data, which can be conducted based on the attribute tables of the shapefiles. Figure 11 aggregates the changes at the building level for seven districts covered by the satellite images (please note that Điện Bàn district was excluded because of very low coverage and small values). It highlights the enormous share of changes in the total building stock of the city. Figure 11. Changed and unchanged buildings for the districts in Da Nang. Figure 12 shows the proportion of new and demolished buildings based on the hexagon dataset, in relation to their average distance to the historic city center of Da Nang. Interestingly, both types of changes show a similar pattern, even if they are not displayed in absolute values. Accordingly, 50% of all new buildings and 43% of all demolitions between 2015 and 2017 occurred within a distance of three to six kilometers from the historic city center. Within a radius of two kilometers, 12.9% of new buildings were built, but 21.2% of demolitions occurred. This again underlines that the central residential areas in Da Nang are already very dense and do not allow much more growth. Instead, the amount of changed buildings (change type 3, as defined in Table 1) is visibly higher in the more central areas. Of all changed buildings in Da Nang, 77% were located within the first six kilometers.

Discussions on Data Quality
This section summarizes the most important aspects concerning accuracy. It is obvious that the quality of the presented dataset depends upon two major error sources: the first is the accuracy of all input products, which determines how well the changes between 2015 and 2017 can be identified and classified. If buildings were not identified or classified for 2015, it is likely that an error is propagated into the change detection approach. As reported by Warth et al. [4], the built-up areas of 2015 have been assessed at an accuracy of 88.4%. Accordingly, the actual size of an individual building triggers potential errors. This affects small buildings in the first instance, because they are more sensitive to misclassifications of small image objects. We did not see a propagation of this error into later products, as the quality of building areas was assessed by Vetter-Gindele et al. [3] and reported with 87.2% accuracy. As the subsequent classification of building types was also based on the buildings' sizes, assessed from the satellite image, incorrect estimates could result in the assignment of a wrong type. This is partly disproved by the confusion matrix shown in Table A1, which shows that most of the misclassifications occurred between buildings of a similar size (e.g., single basic and single local, bungalows and villas, or also multifamily local and multifamily modern). Accordingly, we see a more plausible reason for the incorrect classifications, due to atypical buildings which do not fit into the typology presented in Table 1, rather than due to small errors in the size of their footprint.
One shortcoming of the of the presented accuracy assessment is the statistical dominance of single family local-type buildings. On the one hand, this type makes up 90% of all buildings in Da Nang and has producer and user accuracy measurements above 95%. On the other hand, the small representation of halls, outbuildings and special buildings allows few conclusions to be drawn about the actual accuracy of these building types, and their low producer accuracies (54.5%, 70.6% and 75.0%), indicate that they are rather underestimated in the dataset in favor of more accurate and complete residential buildings and an overall classification accuracy of 90.9%.
However, the chosen approach of visual identification of changes between 2015 and 2017 allowed the interpreters to correct errors in the building type classification of 2015, thereby minimizing misclassifications above the statistically reported accuracy. This leads to the second potential error source of this dataset, which is the subjective interpretation of building types for the image of 2017 and the resulting changes. We tried to minimize the impact of the single interpreter's perception by distributing this task over a broad range of individuals and by cross-checking random buildings to identify systematic errors. Despite this, interpretation errors can still arise from spectral differences in the two images (Figure 3), or in cases where actual changes are overlooked. Accordingly, these errors of omission and commission cannot entirely be avoided. Based on the statistical tests and the considerations above, their impact on the aggregated numbers at an administrative level and the hexagonal grids is considered very small.     [19].