Mapping and characterizing thermal dilation of civil infrastructures with multi-temporal X-band synthetic Aperture Radar interferometry

: Temperature variation plays a signiﬁcant role in the long-term structural behaviour of civil infrastructures, but very few quantitative studies have measured and analysed the infrastructures’ global thermal dilation because of their large sizes and geometric complexities. The modern Differential Synthetic Aperture Radar Interferometry (DInSAR) technique has great potential in applications of their thermal dilation mapping and characterization due to the techniques’ unique capabilities for use in large areas, with high-resolution, and at low-costs for deformation measurements. However, the practical application of DInSAR in thermal dilation estimation is limited by difﬁculty in the precise separation from the residual topographic phase and the trend deformation phase. Moreover, due to a lack of thermal dilation characteristics analyses in previous studies, the thermal dilation mechanisms are still unclear to users, which restricts the accurate understanding of the thermal dilation evolution process. Given the above challenges, an advanced multi-temporal DInSAR approach is proposed in this study, and the effectiveness of this approach was presented using three cases studies concerning different infrastructure types. In this method, the coherent, incoherent, and semantic information of structures were combined in order to reﬁne the detection of point-like targets (PTs). Interferometric subsets with small temporal baselines and temperature differences were used for high-resolution topography estimation. A pre-analysis was adopted to determine the transmission direction, spatial pattern, and temporal variation of the thermal dilation. Then, both the traditional least squares estimation and our robust coherence-weighted least squares regression analysis were performed between the time series displacements and the corresponding temperatures to quantitatively estimate the thermal dilation model. The results were veriﬁed in terms of the estimated linear thermal dilation coefﬁcient, which indicates the greater reliability of our method. Furthermore, the thermal dilation characteristics of different civil infrastructure types were analysed, which facilitates a greater understanding of the thermal dilation evolution process of civil infrastructures.


Introduction
The widespread deterioration and the inadequate upkeep of existing civil infrastructures, such as buildings, highways, and bridges, has accelerated the growing costs of their replacements and repairs, and also brought significant risks to public safety and quality of life [1,2]. To ensure the safe operations of these infrastructures and to reduce the maintenance costs, the consecutive stability monitoring of such manmade structures in an economical manner is important, especially in densely inhabited areas, to detect the deformation at an early stage and to prevent the associated risks by taking timely measures [3,4]. However, it is difficult to achieve this goal using current in situ methods for civil infrastructure monitoring because most in situ techniques (e.g., leveling and total station surveys) are labour-intensive and not frequently performed [5,6]. Currently, professionals and technical officers are seeking sustainable solutions for the routine inspections of civil infrastructures with reduced human resources and allocated funds.
According to previous studies, the temperature variation is considered to have a more significant effect on the deformation of civil infrastructures than the effects of operation loads or structural damage [19,26,27]. Given the high sensitivity of the X-band SAR data on thermal expansion, the expanding or contracting of materials owing to temperature variation would cause a strong seasonal variation in the DInSAR time series displacements, especially when metallic components comprise most of the structure [19,[28][29][30][31][32][33][34]. Moreover, during the multi-temporal DInSAR monitoring, when the observation period does not correspond to exactly one year or an integral multiple of one year, this periodical displacement would also lead to a non-zero deformation signal of point-like targets (PTs) [30][31][32]. As a result, the thermal dilation of civil infrastructures must be incorporated into DInSAR processing and interpretation in order to reliably detect infrastructural trend deformation.
Although many studies have dealt with thermal dilation estimation during time series DInSAR monitoring [19,[28][29][30][31][32][33][34], a number of technical challenges also appeared. First, only a few studies have considered the resolution gap between the X-band SAR images (approximately 3 m) and the old Shuttle Radar Topography Mission (SRTM) (approximately 30 m) [20,35], which would cause uncertainties in the estimated topography, especially for infrastructures newly built after 2000, which affects the accuracy of the thermal dilation and trend deformation mapping. In addition, most of the previous studies ceased at the separation of the thermal dilation, rather than going deep into the investigation of thermal dilation characteristics, such as the transmission direction, spatial pattern, and magnitude, which means that the thermal dilation mechanisms remain unclear to users. Thus, challenges in better understanding the thermal dilation evolution process of civil infrastructures continue.
To address the existing issues that are illustrated above, a multi-temporal DInSAR analysis approach that is aimed at mapping and characterizing the thermal dilation of different infrastructure types, using denser PTs and more accurate topographic phase estimation, is proposed and applied in this study. In this method, the coherent, incoherent, and semantic information of structures were comprehensively considered to improve the density and accuracy of detectable PTs upon structures. Moreover, a high-resolution topography was estimated from the interferometric subsets with small temporal baselines and temperature differences, which eliminates the possible uncertainties of topography estimation and the effects of uncertainties on deformation mapping. Then, a two-step analysis method, which consists of a model analysis phase and a model generation phase, Remote Sens. 2018, 10, 941 3 of 22 was implemented to estimate the thermal dilation. First, the time series differential interferometric phases were visually pre-analysed in order to qualitatively investigate the possible transmission direction, spatial pattern, and temporal variation of the thermal dilation. Then, the pre-analysed results were used in the model generation phase. With the assumption that temperature is homogeneous along the structure, both the traditional least squares regression analysis and the robust coherence-weighted least squares regression analysis that is proposed in this study were implemented to establish the thermal dilation model. Finally, the thermal dilation models estimated from the two methods were verified by comparing the estimated linear thermal dilation coefficients against the actual physical properties of materials, which indicates the greater reliability of our method.
Based on the above approach, three case studies concerning different infrastructures types, including the skyscraper buildings, low buildings, and linear highways, were carried out and are discussed in this article. The experimental results demonstrated the effectiveness of the proposed method for the thermal dilation mapping of different infrastructure types, and the deformation maps of thermal dilation and trend deformation were separately produced. The thermal dilation characteristics, such as the transmission direction, spatial pattern, and magnitude were captured by comparing the thermal dilation results with the different structural information, including specific structural shapes, mechanical properties, construction materials, and sizes. Our results demonstrated that the transmission direction of the thermal dilation is highly dependent on the specific structural shape of the target, which is usually along the direction of the longest side of the geometric shape. Moreover, the static structural characteristics of infrastructures, such as the positions of structural key points, can be approximately inferred from the spatial pattern of thermal dilation. Finally, the magnitude of thermal dilation is proven to be highly variable for each infrastructure, which is associated with the material properties and sizes of the objects. These analytical findings provide conclusive evidence of the complex thermal dilation of different infrastructure types and aids in understanding the infrastructures' thermal dilation evolution process. This paper is structured as follows: Section 2 describes the multi-temporal DInSAR analysis method that is applied in this study. The basic information of study objects, their datasets, and preliminary results are introduced in Section 3. Then, a discussion of the thermal dilation characteristics on different civil infrastructure types is carried out in Section 4. The conclusions are given in Section 5.

Methodology
The main workflow of the multi-temporal DInSAR approach that is applied in this paper is illustrated in Figure 1, in which some steps have been improved to refine the thermal dilation monitoring of civil infrastructure using X-band high-resolution SAR images. The main idea of the improved steps marked by the coloured rectangles, including the refined point-like target detection, topographic phase estimation and displacement retrieval, as well as the thermal dilation mapping and characterization, are described in detail in the following.

Refined Point-Like Target Detection
In practical applications, it is important to maintain the high level of density and accuracy of the selected PTs when monitoring the detailed deformation of ground targets [16][17][18]. This is because detecting as many PTs as possible on the structures is conducive to reducing the calculation error, as well as improving the estimation accuracy in DInSAR monitoring [7,8,12]. Moreover, since the backscattering signals of different ground targets in the complex urban scenery are usually mixed together due to oblique scene illumination, it is difficult to separate these PTs using the traditional method. Therefore, new methods that are aimed at improving the density and accuracy of detectable PTs upon the civil infrastructures are required for this research. In this study, the coherent and incoherent information derived from the high-resolution SAR images, as well as the detailed semantic information about specific structures are fully integrated into the multi-temporal DInSAR processing to refine detection of PTs. Since the higher coherence of SAR images usually indicates a higher quality of the interferometric phase, the temporal coherence, which is an indicator of phase stability, was used to select a set of PT candidates [12,13]. On the other hand, the statistical values of backscattering signals, including the mean intensity and the amplitude deviation of the time series SAR images, were also used as independent incoherent information to identify another set of PT candidates. Then, the PT candidates selected from both the coherent and incoherent information are merged to maximize the number of detectable PTs, which is greater than the number of PTs selected by each method alone. The object cognitive abilities of the coherent and incoherent information in the high-resolution SAR images are complemented in this study.
In this study, the specific semantic information of the infrastructures such as geometric shape, geolocation buffer zone, and elevation were integrated together into the PTs refinement [31,36]. Since the higher coherence of SAR images usually indicates a higher quality of the interferometric phase, the temporal coherence, which is an indicator of phase stability, was used to select a set of PT candidates [12,13]. On the other hand, the statistical values of backscattering signals, including the mean intensity and the amplitude deviation of the time series SAR images, were also used as independent incoherent information to identify another set of PT candidates. Then, the PT candidates selected from both the coherent and incoherent information are merged to maximize the number of detectable PTs, which is greater than the number of PTs selected by each method alone. The object cognitive abilities of the coherent and incoherent information in the high-resolution SAR images are complemented in this study.
In this study, the specific semantic information of the infrastructures such as geometric shape, geolocation buffer zone, and elevation were integrated together into the PTs refinement [31,36]. First, the PTs located within a local buffer zone of the structures' geolocation in the horizontal plane were chosen as PTs candidates. Then, the outlier points and target points were distinguished according to their local elevation analysis, which improves the accuracy of the structural PTs. Since the PTs Remote Sens. 2018, 10, 941 5 of 22 elevation selected on the civil structures should be consecutive, if the difference between the mean elevation of surrounding PTs and the central PT is less than double of their root means squares, the PT is thought to be the target point. Otherwise, the PT would be removed as an outlier point.

Topographic Phase Estimation and Displacement Retrieval
After PTs refinement, the observed phase of the ith PT in the jth interferogram can be expressed as Equation (1) [7,9,13]: In the above equation, ∅ int stands for the interferogram phase of the PT; ∅ de f indicates the phase change caused by displacement in the line-of-sight (LOS) direction, ∅ top represents the topographic phase caused by the height of target; ∅ atm accounts for the phase contributed by the difference in atmospheric retardation between image pair; and ∅ orbit corresponds to the phase due to orbit inaccuracy. The noise term ∅ noise is small enough for a PT and does not significantly obscure the signal [9,10].
The difference in resolution between the SRTM Digital Elevation Model (DEM) and TerraSAR-X data might lead to the inaccurate estimation of the topographic phase, especially for modern civil infrastructures. Therefore, the SRTM DEM must be used cautiously, especially in rapidly developing cities like Shanghai and Hong Kong, and methods of topographic phase estimation should be further improved when aiming at the accurate thermal dilation monitoring of civil infrastructures [20,35]. In this study, the interferometric subsets with temperature differences of less than 2.5 • C and temporal baselines of less than 100 days were used for the high-resolution topography estimation. Although the temperature difference of 2.5 • C may cause thermal dilation, when the temperature is close to the reference temperature (without thermal dilation), the temperature of the infrastructures does not change much with such a small temperature difference and the corresponding thermal dilation can be ignored. Therefore, the deformation phase is supposed to be small, and the atmosphere and noise phases can be separated during the time series processing. Then, the high-resolution topography of the study area can be obtained and used as a high-resolution DEM in order to generate a set of differential interferograms from the total interferograms.
Then, a three-dimensional (3D) time-space unwrapping algorithm is applied to estimate the final deformation of each PT [9,11,13]. The observed deformation, including the linear and nonlinear displacements, can be decomposed into the vertical and horizontal directions for further analysis. The deformation measurement of DInSAR (D Los ) is the projection of the original displacement D with the components D Ver and D Hor in the vertical and horizontal directions, respectively, on the LOS direction. Assuming that the local incidence angle of the satellite is θ, the LOS displacement can be decomposed by the following equation.
Actually, high-resolution SAR datasets, including both ascending and descending tracks, to model a reliable 3D deformation are not always available. Therefore, based on the assumption that the civil infrastructure deformation is mainly caused by the thermal dilation, and the thermal dilation was always transmitted along one direction (either vertically or horizontally), which was also proven in many previous studies [26,37], a simplified conversion of D Los from the LOS direction into the vertical or horizontal directions displacements can be implemented by the following equations: Furthermore, for horizontal civil infrastructures, such as highways, the displacement along the longitudinal direction (D Lon ) can be calculated by the following equation: Here, α is the horizontal angle between the structure orientation and the LOS direction, as illustrated in Figure 2.
Remote Sens. 2018, 10, x FOR PEER REVIEW 6 of 22 (4) Furthermore, for horizontal civil infrastructures, such as highways, the displacement along the longitudinal direction ( ) can be calculated by the following equation: Here, is the horizontal angle between the structure orientation and the LOS direction, as illustrated in Figure 2.

Thermal Dilation Mapping and Characterization
Since most modern civil infrastructures are constructed by steel or reinforced concrete material, the change in ambient temperature would cause the thermal expansion and contraction of the structure and material, which would even exceed the magnitude of other phase components [20][21][22]37]. This phenomenon, if not properly handled during the DInSAR data processing and analysis, would have a significant impact on the interferometric products, particularly on the deformation velocity map, and would challenge the widespread application of DInSAR technology for civil infrastructure monitoring [19,[28][29][30][31][32][33][34]37].
Most previous studies modelled thermal dilation by adding a temperature dependent phase term into the traditional two-dimensional (linear velocity and topographic error) analysis model, or by using a periodic seasonal phase term [19,28,29]. In this study, the deformation signals were estimated through temporal-spatial filtering using the StamPS method [12,13], rather than through calculation with a prior deformation model hypothesis, such as a linear model [7]. Therefore, an additional two-step thermal dilation analysis method was applied to the deformation signals after the estimation of other phase components [31,37].
The two-step thermal dilation analysis method consists of a model analysis phase and a model generation phase. The model analysis phase qualitatively estimates the characteristics of thermal dilation for structures. In our experiments, this step is carried out by visually observing and pre-analysing the characteristics of the time series differential interferometric phases. Critical features that can be detected from the time series differential interferograms are the possible transmission direction, spatial pattern, and temporal variation of thermal dilation.
Then, the above qualitative analysis results were used in the model generation phase. The correlation pattern between the measured time series deformation along the thermal dilation transmission direction and the corresponding temperature variation was quantitatively established with a robust coherence-weighted least squares linear regression fitting [31,37]. In most of the previous empirical linear models, the regression parameters were commonly computed in order to minimize an error function based on least squares estimates, as in Equation (6):

Thermal Dilation Mapping and Characterization
Since most modern civil infrastructures are constructed by steel or reinforced concrete material, the change in ambient temperature would cause the thermal expansion and contraction of the structure and material, which would even exceed the magnitude of other phase components [20][21][22]37]. This phenomenon, if not properly handled during the DInSAR data processing and analysis, would have a significant impact on the interferometric products, particularly on the deformation velocity map, and would challenge the widespread application of DInSAR technology for civil infrastructure monitoring [19,[28][29][30][31][32][33][34]37].
Most previous studies modelled thermal dilation by adding a temperature dependent phase term into the traditional two-dimensional (linear velocity and topographic error) analysis model, or by using a periodic seasonal phase term [19,28,29]. In this study, the deformation signals were estimated through temporal-spatial filtering using the StamPS method [12,13], rather than through calculation with a prior deformation model hypothesis, such as a linear model [7]. Therefore, an additional two-step thermal dilation analysis method was applied to the deformation signals after the estimation of other phase components [31,37].
The two-step thermal dilation analysis method consists of a model analysis phase and a model generation phase. The model analysis phase qualitatively estimates the characteristics of thermal dilation for structures. In our experiments, this step is carried out by visually observing and pre-analysing the characteristics of the time series differential interferometric phases. Critical features that can be detected from the time series differential interferograms are the possible transmission direction, spatial pattern, and temporal variation of thermal dilation.
Then, the above qualitative analysis results were used in the model generation phase. The correlation pattern between the measured time series deformation along the thermal dilation transmission direction and the corresponding temperature variation was quantitatively established with a robust coherence-weighted least squares linear regression fitting [31,37]. In most of the previous empirical linear models, the regression parameters were commonly computed in order to minimize an error function based on least squares estimates, as in Equation (6): In the above equation, n is the number of observations, E i indicates the measurement, andÊ i represents the displacement that is predicted by the estimated model. However, this criterion is known to be sensitive to the presence of outliers in the datasets and therefore may lead to unreliable models [38,39]. In this study, we mitigate this problem by employing a robust weighted least squares function in which weights are assigned to the time series measurements. Since the interferograms with higher coherence would produce more reliable measurements, the use of the interferometric coherence of interferometric pairs to determine the weight of the corresponding time series observations (given by Equation (7)) is simple and effective.γ Here,γ i indicates the average coherence of the ith interferogram. The coherence of a pixel in the ith interferogram is calculated by the following equation: In the Equation (8), N indicates the number of pixels in an estimation window; y In this study, the estimated thermal dilations were validated in two ways. Qualitatively, their spatial patterns were compared with the patterns predicted from the specific structural mechanical properties of infrastructures. Quantitatively, the thermal dilation coefficients that were estimated from both the traditional method and our method using Equation (9) were compared with the actual physical properties of materials, which indicates the greater reliability of our method.
where D T indicates the magnitude of the thermal dilation along its transmission direction; α T is the thermal dilation coefficient of the material; and, L is the effective transmission length of the thermal dilation. Moreover, the thermal dilation characteristics, in terms of their transmission direction, spatial pattern, and magnitude were analysed to identify their potential correlation with the specific structural characteristics of different infrastructure types. These analytical findings aid in better understanding the thermal dilation evolution process of different civil infrastructure types and are valuable in optimizing their future maintenance strategies.

Case Studies and Preliminary Results
The method described above has been applied to three case studies concerning different infrastructure types, including the Shanghai World Financial Center (SWFC) building, the Shanghai New International Expo Center (SNIEC) buildings, and the Hong Kong Tsing Kwai Highway (TKH). For each case, we briefly introduced the SAR data and interpreted the thermal dilation result in detail.

The SWFC and SAR Data
The SWFC is a skyscraper building with the world's highest viewing platform of approximately 492 m, which ranks as the seventh highest building in the world. The optical picture and the SAR amplitude image of this building are shown in Figure 3a,b, respectively. In this case, a stack of 24 descending TerraSAR-X Stripmap mode images with approximately 3 m resolution was collected from 2014 to 2016 over Shanghai. The coverage of the TerraSAR-X images and the location of the SWFC are illustrated in Figure 3c. Leveling points F59 and F16 near the building (see triangles in Figure 3c) were used to validate the DInSAR results [21]. The basic SAR image information is described in Table 1, and the temperatures were acquired from the China Meteorological Data Sharing Service Website.

The SWFC and SAR Data
The SWFC is a skyscraper building with the world's highest viewing platform of approximately 492 m, which ranks as the seventh highest building in the world. The optical picture and the SAR amplitude image of this building are shown in Figure 3a,b, respectively. In this case, a stack of 24 descending TerraSAR-X Stripmap mode images with approximately 3 m resolution was collected from 2014 to 2016 over Shanghai. The coverage of the TerraSAR-X images and the location of the SWFC are illustrated in Figure 3c. Leveling points F59 and F16 near the building (see triangles in Figure 3c) were used to validate the DInSAR results [21]. The basic SAR image information is described in Table 1, and the temperatures were acquired from the China Meteorological Data Sharing Service Website.

Thermal Dilation Results
Based on our method, the estimated topography is close to the real height of 492 m and it is considered to be reliable. The time series differential interferometric phases of the SWFC are obtained and shown in Figure 4. Interferometric fringes with different densities along the building were observed in the differential interferograms with various temperature differences (∆T). Specifically, when the ∆T is large, several fringes can be observed in the interferograms, which are shown by the red rectangles in Figure 4. This periodical phase changing disappeared or became slower in the differential interferograms with a smaller ∆T, implying that these periodical fringes are mainly affected by the change in temperature. Moreover, the thermal dilation of the SWFC is most likely to propagate in the vertical direction because the fringes evidently transmitted along the vertical direction of the building. Based on our method, the estimated topography is close to the real height of 492 m and it is considered to be reliable. The time series differential interferometric phases of the SWFC are obtained and shown in Figure 4. Interferometric fringes with different densities along the building were observed in the differential interferograms with various temperature differences ( ∆ ). Specifically, when the ∆ is large, several fringes can be observed in the interferograms, which are shown by the red rectangles in Figure 4. This periodical phase changing disappeared or became slower in the differential interferograms with a smaller ∆ , implying that these periodical fringes are mainly affected by the change in temperature. Moreover, the thermal dilation of the SWFC is most likely to propagate in the vertical direction because the fringes evidently transmitted along the vertical direction of the building. . The time series differential interferograms of the SWFC, the red rectangles indicate the differential interferograms with larger temperature differences and denser fringes.
Since the thermal dilation of the SWFC is mainly along the vertical direction, the time series vertical deformations of this building were calculated and are shown in  . The time series differential interferograms of the SWFC, the red rectangles indicate the differential interferograms with larger temperature differences and denser fringes.
Since the thermal dilation of the SWFC is mainly along the vertical direction, the time series vertical deformations of this building were calculated and are shown in Figure 5. The PTs distribution achieved an excellent coverage of the entire building, thus ensuring a robust estimation of the deformation parameters. From the structural point of view, the bottom of the building is supposed to be a fixed point because the displacement at the bottom of the building is very small and gradually increases to the top. Periodical displacements of the building can be observed during the observation period, especially on the top half of it, with values ranging up to approximately 90 mm. The positive values indicate that the rising temperature induces an increase in the building's height, leading to the displacement towards the satellite, and vice versa.
We obtained the vertical length variation of the building by calculating the average deformation difference (∆ ) between the bottom and top PTs (blue dots in Figure 6), and found that it is highly correlated with the acquisition temperature (T). The standard deviation of the measurements calculated at each acquisition time is shown as an error bar that is marked by an orange line. Both the traditional least squares estimation and the robust coherence-weighted least squares regression analysis proposed in this study were implemented to identify the relationship between these two variables. The two linear fitting models established by the traditional method (green dotted line) and our method (red solid line), which are shown in Figure 6, clearly highlighted the evident linear correlation between the ∆ and T.  According to these two models, when ∆ equals to zero, the corresponding is 20.48 °C for the traditional method and 19.30 °C for our method, which are both close to the temperature of the master image (21.9 °C). Although the fitted reference temperature (without thermal dilation) of the traditional method is closer to 21.9 °C than our method, the accuracy of only one measurement was considered. However, deformation measurements that are close to the reference temperature would be less accurate than those with larger temperature differences. Actually, when the temperature is approximately 21.9 °C (like 3 °C), the temperature of the building does not change that much. Thus, the corresponding thermal dilation will be small, and the accuracy would be low. When the temperature differences are larger than 10 °C, the building will undergo significant temperature changes, and the thermal dilation is evident and not susceptible to noise. Therefore, we use the linear thermal dilation coefficients of the construction material calculated from the fitted models We obtained the vertical length variation of the building by calculating the average deformation difference (∆D) between the bottom and top PTs (blue dots in Figure 6), and found that it is highly correlated with the acquisition temperature (T). The standard deviation of the measurements calculated at each acquisition time is shown as an error bar that is marked by an orange line. Both the traditional least squares estimation and the robust coherence-weighted least squares regression analysis proposed in this study were implemented to identify the relationship between these two variables. The two linear fitting models established by the traditional method (green dotted line) and our method (red solid line), which are shown in Figure 6, clearly highlighted the evident linear correlation between the ∆D and T. We obtained the vertical length variation of the building by calculating the average deformation difference (∆ ) between the bottom and top PTs (blue dots in Figure 6), and found that it is highly correlated with the acquisition temperature (T). The standard deviation of the measurements calculated at each acquisition time is shown as an error bar that is marked by an orange line. Both the traditional least squares estimation and the robust coherence-weighted least squares regression analysis proposed in this study were implemented to identify the relationship between these two variables. The two linear fitting models established by the traditional method (green dotted line) and our method (red solid line), which are shown in Figure 6, clearly highlighted the evident linear correlation between the ∆ and T.  According to these two models, when ∆ equals to zero, the corresponding is 20.48 °C for the traditional method and 19.30 °C for our method, which are both close to the temperature of the master image (21.9 °C). Although the fitted reference temperature (without thermal dilation) of the traditional method is closer to 21.9 °C than our method, the accuracy of only one measurement was considered. However, deformation measurements that are close to the reference temperature would be less accurate than those with larger temperature differences. Actually, when the temperature is approximately 21.9 °C (like 3 °C), the temperature of the building does not change that much. Thus, the corresponding thermal dilation will be small, and the accuracy would be low. When the temperature differences are larger than 10 °C, the building will undergo significant temperature changes, and the thermal dilation is evident and not susceptible to noise. Therefore, we use the linear thermal dilation coefficients of the construction material calculated from the fitted models According to these two models, when ∆D equals to zero, the corresponding T is 20.48 • C for the traditional method and 19.30 • C for our method, which are both close to the temperature of the master image (21.9 • C). Although the fitted reference temperature (without thermal dilation) of the traditional method is closer to 21.9 • C than our method, the accuracy of only one measurement was considered. However, deformation measurements that are close to the reference temperature would be less accurate than those with larger temperature differences. Actually, when the temperature is approximately 21.9 • C (like ±3 • C), the temperature of the building does not change that much. Thus, the corresponding thermal dilation will be small, and the accuracy would be low. When the temperature differences are larger than 10 • C, the building will undergo significant temperature changes, and the thermal dilation is evident and not susceptible to noise. Therefore, we use the linear thermal dilation coefficients of the construction material calculated from the fitted models (8.64 × 10 −6 per • C for the traditional method and 10.10 × 10 −6 per • C for our method), which consider the total accuracy of all measurements, to evaluate the results. Since the actual thermal dilation coefficient is taken as 12 × 10 −6 per • C for pure steel and 9 × 10 −6 per • C for pure concrete [40], the coefficient derived from our method (10.10 × 10 −6 per • C) is within the reasonable range of the physical properties of typical steel and reinforced concrete material (9 ∼ 12 × 10 −6 per • C), while the coefficient derived from the traditional method is lower than the theoretical value. According to the thermal dilation coefficient validation, our method is proven to be more effective than the traditional method.
After calibrating and validating based on the local leveling points, the trend deformation and thermal dilation of the SWFC are obtained and are shown in Figure 7a,b, respectively. After separating the thermal dilation, this building is quite stable with a slight and uniform subsidence velocity of less than 5 mm/year, which indicates that the observed time series displacement is mainly caused by thermal dilation. The temperature deformation parameters are around zero at the bottom of the building and gradually increased from the bottom to approximately 5 mm/ • C at the top. consider the total accuracy of all measurements, to evaluate the results. Since the actual thermal dilation coefficient is taken as 12 10 per °C for pure steel and 9 10 per °C for pure concrete [40], the coefficient derived from our method ( 10.10 10 per °C) is within the reasonable range of the physical properties of typical steel and reinforced concrete material (9 ∼ 12 10 per °C), while the coefficient derived from the traditional method is lower than the theoretical value. According to the thermal dilation coefficient validation, our method is proven to be more effective than the traditional method.
After calibrating and validating based on the local leveling points, the trend deformation and thermal dilation of the SWFC are obtained and are shown in Figure 7a,b, respectively. After separating the thermal dilation, this building is quite stable with a slight and uniform subsidence velocity of less than 5 mm/year, which indicates that the observed time series displacement is mainly caused by thermal dilation. The temperature deformation parameters are around zero at the bottom of the building and gradually increased from the bottom to approximately 5 mm/°C at the top.

The SNIEC and SAR Data
The SNIEC is the international exhibition venue in the heart of Shanghai, which consists of 17 exhibition halls, including W1 to W5 in the west, N1 to N5 in the north, and E1 to E7 in the east (see Figure 8b). Each exhibition hall is approximately 185 m long and 70 m wide, with an entire indoor area of 200,000 square meters and an outdoor area of 130,000 square meters. The amplitude image and the optical image of these buildings are shown in Figure 8a,b, respectively. Each exhibition hall is clearly visible in the high-resolution SAR image. The same SAR dataset (see Figure 3c) with the SWFC was used for the thermal dilation monitoring, and the location of the SNIEC is represented as the green star in Figure 3c.

The SNIEC and SAR Data
The SNIEC is the international exhibition venue in the heart of Shanghai, which consists of 17 exhibition halls, including W1 to W5 in the west, N1 to N5 in the north, and E1 to E7 in the east (see Figure 8b). Each exhibition hall is approximately 185 m long and 70 m wide, with an entire indoor area of 200,000 square meters and an outdoor area of 130,000 square meters. The amplitude image and the optical image of these buildings are shown in Figure 8a,b, respectively. Each exhibition hall is clearly visible in the high-resolution SAR image. The same SAR dataset (see Figure 3c) with the SWFC was used for the thermal dilation monitoring, and the location of the SNIEC is represented as the green star in Figure 3c.

Thermal Dilation Results
The time series differential interferograms of SNIEC are illustrated in Figure 9, where the denser periodical fringes appeared at the top of exhibition halls in the differential interferograms with larger temperature differences (remarked by red rectangles). Since there are no obvious elevation differences between the two sides of the exhibition halls, and the distance between the two sides is less than 200 m, these fringes cannot be caused by the topographic phase or the atmosphere phase. According to the temperature correlation analysis of each differential interferogram, the thermal dilation is probably the main cause of these fringes. If the thermal dilation of the exhibition halls was mainly transmitted along the vertical direction, the thermal dilation magnitude at the top of buildings with almost the same height should be the same. However, the magnitudes of the thermal dilation at the top of the exhibition halls are different from one side to the other side in the differential interferograms, which indicates that the thermal dilation is most likely to propagate along the horizontal direction. Figure 9. The time series differential interferograms of the SNIEC, the red rectangles indicate the differential interferograms with larger temperature differences and denser fringes.

Thermal Dilation Results
The time series differential interferograms of SNIEC are illustrated in Figure 9, where the denser periodical fringes appeared at the top of exhibition halls in the differential interferograms with larger temperature differences (remarked by red rectangles). Since there are no obvious elevation differences between the two sides of the exhibition halls, and the distance between the two sides is less than 200 m, these fringes cannot be caused by the topographic phase or the atmosphere phase. According to the temperature correlation analysis of each differential interferogram, the thermal dilation is probably the main cause of these fringes. If the thermal dilation of the exhibition halls was mainly transmitted along the vertical direction, the thermal dilation magnitude at the top of buildings with almost the same height should be the same. However, the magnitudes of the thermal dilation at the top of the exhibition halls are different from one side to the other side in the differential interferograms, which indicates that the thermal dilation is most likely to propagate along the horizontal direction.
Based on the above analysis, we assumed that the thermal dilation of the exhibition halls was mainly transmitted in the horizontal direction. The time series east-west displacements of the SNIEC were then calculated and they are shown in Figure 10. First, we can see that the number of detectable PTs depends on their projection lengths in the north-south direction, which is related to the satellite's flying direction (nearly north-south). The most PTs were selected from the exhibition halls N1 to N5 since they were nearly along the north-south direction, while fewer PTs were identified at the tops of W1 to W5 and E1 to E7 due to their shorter projection lengths along the flying direction of the satellite, especially for the five exhibition halls in the west.
According to Figure 10, the deformation in the middle of the building is small and gradually increased to both sides, with values ranging up to approximately 20 mm, but in opposing directions. Since the transmission distances of the thermal dilation in the east-west direction are different for different exhibition halls, their thermal dilation magnitudes are also slightly varied. The thermal dilation of the north exhibition halls (N1 to N5) is smaller when compared with the west (W1 to W5) and east (E1 to E7) halls, due to the smaller scales along the east-west direction. When we further converted the east-west deformation to the length direction of each exhibition hall, a linear relationship between the average displacements in the length direction of exhibition halls (see blue dots in Figure 11) and the corresponding temperature variation was found.
interferogram, the thermal dilation is probably the main cause of these fringes. If the thermal dilation of the exhibition halls was mainly transmitted along the vertical direction, the thermal dilation magnitude at the top of buildings with almost the same height should be the same. However, the magnitudes of the thermal dilation at the top of the exhibition halls are different from one side to the other side in the differential interferograms, which indicates that the thermal dilation is most likely to propagate along the horizontal direction. Figure 9. The time series differential interferograms of the SNIEC, the red rectangles indicate the differential interferograms with larger temperature differences and denser fringes. Figure 9. The time series differential interferograms of the SNIEC, the red rectangles indicate the differential interferograms with larger temperature differences and denser fringes.
Based on the above analysis, we assumed that the thermal dilation of the exhibition halls was mainly transmitted in the horizontal direction. The time series east-west displacements of the SNIEC were then calculated and they are shown in Figure 10. First, we can see that the number of detectable PTs depends on their projection lengths in the north-south direction, which is related to the satellite's flying direction (nearly north-south). The most PTs were selected from the exhibition halls N1 to N5 since they were nearly along the north-south direction, while fewer PTs were identified at the tops of W1 to W5 and E1 to E7 due to their shorter projection lengths along the flying direction of the satellite, especially for the five exhibition halls in the west.
According to Figure 10, the deformation in the middle of the building is small and gradually increased to both sides, with values ranging up to approximately 20 mm, but in opposing directions. Since the transmission distances of the thermal dilation in the east-west direction are different for different exhibition halls, their thermal dilation magnitudes are also slightly varied. The thermal dilation of the north exhibition halls (N1 to N5) is smaller when compared with the west (W1 to W5) and east (E1 to E7) halls, due to the smaller scales along the east-west direction. When we further converted the east-west deformation to the length direction of each exhibition hall, a linear relationship between the average displacements in the length direction of exhibition halls (see blue dots in Figure 11) and the corresponding temperature variation was found.
The estimated thermal dilation models of the SNIEC are shown in Figure 11, where the green dotted line indicates the model derived from the traditional method and the red solid line indicates the model that was calculated using our method. The error bars, which represented by orange lines, indicate the standard deviation of measurements at each acquisition. According to these models, when the air temperature is approximately 21.52 °C in the traditional model or 21.32 °C in our model, there is no thermal dilation detected. Both of the models are close to the measured value of 21.9 °C. Based on this temperature, when the temperature changed by 1 °C, the thermal dilation of 8.77 10 (traditional model) and 9.80 10 (our model) of the building's length will be detected. The coefficient estimated from our method (9.80 10 per °C) is also in accordance with the linear thermal dilation coefficient of steel and reinforced concrete material, while that of the traditional method (8.77 10 per °C) is still slightly smaller than the theoretical value, which implies the greater reliability of our method.  The estimated thermal dilation models of the SNIEC are shown in Figure 11, where the green dotted line indicates the model derived from the traditional method and the red solid line indicates the model that was calculated using our method. The error bars, which represented by orange lines, indicate the standard deviation of measurements at each acquisition. According to these models, when the air temperature is approximately 21.52 • C in the traditional model or 21.32 • C in our model, there is no thermal dilation detected. Both of the models are close to the measured value of 21.9 • C. Based on this temperature, when the temperature changed by 1 • C, the thermal dilation of 8.77 × 10 −6 (traditional model) and 9.80 × 10 −6 (our model) of the building's length will be detected. The coefficient estimated from our method (9.80 × 10 −6 per • C) is also in accordance with the linear thermal dilation coefficient of steel and reinforced concrete material, while that of the traditional method (8.77 × 10 −6 per • C) is still slightly smaller than the theoretical value, which implies the greater reliability of our method. After calibrating and validating based on the local leveling points, the trend deformation and thermal dilation of the SNIEC buildings are illustrated in Figure 12a,b, respectively. The magnitude of the temperature deformation parameters on the exhibition halls increased from the middle of the buildings to approximately 1 mm/°C at both sides, but with opposing directions. From such a thermal dilation spatial pattern, we can infer that, structurally, the fixed points of the buildings are probably located in the middle of the roofs, and the thermal dilation transmitted from the middle to both sides reached the peak at both ends. The trend deformations of these buildings are small, ranging from −4 to 4 mm/year, which may be caused by the surface subsidence of Shanghai. Therefore, the thermal dilation is the main cause of the observed time series deformation. The slightly different trend deformation of the exhibition halls may be due to their different orientations.

The TKH and SAR Data
The TKH in Hong Kong is an elevated road, which is approximately 20 m above the ground. The optical image and SAR amplitude image of TKH are shown in Figure 13a,b, respectively. The blue rectangles in Figure 13a indicate the locations of the expansion joints of the TKH. In this study, we collected 28 TerraSAR-X Stripmap images from 2013 to 2014 in order to monitor the highway's After calibrating and validating based on the local leveling points, the trend deformation and thermal dilation of the SNIEC buildings are illustrated in Figure 12a,b, respectively. The magnitude of the temperature deformation parameters on the exhibition halls increased from the middle of the buildings to approximately 1 mm/ • C at both sides, but with opposing directions. From such a thermal dilation spatial pattern, we can infer that, structurally, the fixed points of the buildings are probably located in the middle of the roofs, and the thermal dilation transmitted from the middle to both sides reached the peak at both ends. The trend deformations of these buildings are small, ranging from −4 to 4 mm/year, which may be caused by the surface subsidence of Shanghai. Therefore, the thermal dilation is the main cause of the observed time series deformation. The slightly different trend deformation of the exhibition halls may be due to their different orientations. After calibrating and validating based on the local leveling points, the trend deformation and thermal dilation of the SNIEC buildings are illustrated in Figure 12a,b, respectively. The magnitude of the temperature deformation parameters on the exhibition halls increased from the middle of the buildings to approximately 1 mm/°C at both sides, but with opposing directions. From such a thermal dilation spatial pattern, we can infer that, structurally, the fixed points of the buildings are probably located in the middle of the roofs, and the thermal dilation transmitted from the middle to both sides reached the peak at both ends. The trend deformations of these buildings are small, ranging from −4 to 4 mm/year, which may be caused by the surface subsidence of Shanghai. Therefore, the thermal dilation is the main cause of the observed time series deformation. The slightly different trend deformation of the exhibition halls may be due to their different orientations.

The TKH and SAR Data
The TKH in Hong Kong is an elevated road, which is approximately 20 m above the ground. The optical image and SAR amplitude image of TKH are shown in Figure 13a,b, respectively. The blue rectangles in Figure 13a indicate the locations of the expansion joints of the TKH. In this study, we collected 28 TerraSAR-X Stripmap images from 2013 to 2014 in order to monitor the highway's

The TKH and SAR Data
The TKH in Hong Kong is an elevated road, which is approximately 20 m above the ground. The optical image and SAR amplitude image of TKH are shown in Figure 13a,b, respectively. The blue rectangles in Figure 13a indicate the locations of the expansion joints of the TKH. In this study, we collected 28 TerraSAR-X Stripmap images from 2013 to 2014 in order to monitor the highway's thermal dilation. The coverage of the SAR images and the location of the highway are illustrated in Figure 13c. The basic information of the SAR images is described in Table 2.
Remote Sens. 2018, 10, x FOR PEER REVIEW 15 of 22 thermal dilation. The coverage of the SAR images and the location of the highway are illustrated in Figure 13c. The basic information of the SAR images is described in Table 2.

Thermal Dilation Results
The time series differential interferometric phases of the TKH are shown in Figure 14. Although the annual temperature difference in Hong Kong is smaller than that of Shanghai, the temperature-dependent deformation can also be detected by the X-Band SAR images. The differential interferograms with temperature differences that are larger than 10 • C showed denser fringes than others, and those fringes changed more slowly in the differential interferograms with smaller temperature differences. Therefore, these fringes are most likely caused by the thermal dilation. According to their spatial pattern, the thermal dilation is probably transmitted along the longitudinal direction of the TKH.
The longitudinal time series displacement maps of the TKH are shown in Figure 15. Dense PTs were selected on the TKH, especially on the northern outline, which is very bright in the SAR amplitude image. In this case, several sudden discontinuities along the highway were observed and are marked by the red lines in Figure 15. Based on further analysis, these discontinuities appear to be located at the expansion joints of the TKH, which are represented in Figure 13a. By referring to the red lines, this highway can be divided into three segments, which correspond to the Cheung Tsing Bridge and two segments of viaducts from the west to the east. The accumulation of thermal dilation in the two terminations of the Cheung Tsing Bridge is evident. The thermal dilations of the two highway segments transmitted from one expansion joint to another. From a structural viewpoint, this structure is defined in the civil engineering framework as a hyperstatic structure, meaning that the high number of constraints (such as expansion joints) tend to hinder the deformation of the bridge and highway [18,28,29]. Therefore, such a spatial pattern of thermal dilation agreed well with the pattern predicted from the positions of expansion joints along the highway, which indicates the reliability of our results.

Thermal Dilation Results
The time series differential interferometric phases of the TKH are shown in Figure 14. Although the annual temperature difference in Hong Kong is smaller than that of Shanghai, the temperature-dependent deformation can also be detected by the X-Band SAR images. The differential interferograms with temperature differences that are larger than 10 °C showed denser fringes than others, and those fringes changed more slowly in the differential interferograms with smaller temperature differences. Therefore, these fringes are most likely caused by the thermal dilation. According to their spatial pattern, the thermal dilation is probably transmitted along the longitudinal direction of the TKH.
The longitudinal time series displacement maps of the TKH are shown in Figure 15. Dense PTs were selected on the TKH, especially on the northern outline, which is very bright in the SAR amplitude image. In this case, several sudden discontinuities along the highway were observed and are marked by the red lines in Figure 15. Based on further analysis, these discontinuities appear to be located at the expansion joints of the TKH, which are represented in Figure 13a. By referring to the red lines, this highway can be divided into three segments, which correspond to the Cheung Tsing Bridge and two segments of viaducts from the west to the east. The accumulation of thermal dilation in the two terminations of the Cheung Tsing Bridge is evident. The thermal dilations of the two highway segments transmitted from one expansion joint to another. From a structural viewpoint, this structure is defined in the civil engineering framework as a hyperstatic structure, meaning that the high number of constraints (such as expansion joints) tend to hinder the deformation of the bridge and highway [18,28,29]. Therefore, such a spatial pattern of thermal dilation agreed well with the pattern predicted from the positions of expansion joints along the highway, which indicates the reliability of our results. Figure 14. The time series differential interferograms of the TKH, the red rectangles indicate the differential interferorgrams with larger temperature differences and denser fringes.
By quantitatively analysing the time series average longitudinal displacements of the Cheung Tsing Bridge (blue dots in Figure 16) and the temperature variation, the thermal dilation model that was derived from the traditional least squares method and our coherence-weighted method for the TKH are illustrated in Figure 16 as the green dotted line and red solid line, respectively. According to these models, when the air temperature is approximately 17.24 °C for the traditional model and 16.92 °C for our model, no thermal dilation is observed on the highway. Both the traditional Figure 14. The time series differential interferograms of the TKH, the red rectangles indicate the differential interferorgrams with larger temperature differences and denser fringes.
By quantitatively analysing the time series average longitudinal displacements of the Cheung Tsing Bridge (blue dots in Figure 16) and the temperature variation, the thermal dilation model that was derived from the traditional least squares method and our coherence-weighted method for the TKH are illustrated in Figure 16 as the green dotted line and red solid line, respectively. According to these models, when the air temperature is approximately 17.24 • C for the traditional model and 16.92 • C for our model, no thermal dilation is observed on the highway. Both the traditional method and our method comprehensively consider all of the measurements at different temperatures to minimize the total errors, rather than only considering the measurement with no thermal dilation. Therefore, the measurements are not always exactly the same as the fitted values, and the slight differences in the no thermal dilation temperature are within a reasonable range. Based on the standard temperature, a temperature variation of 1 • C may lead to a change of 8.15 × 10 −6 and 9.37 × 10 −6 of the bridge's length, according to the traditional model and our model. The linear thermal dilation coefficient estimated from our method (9.37 × 10 −6 per • C) is within the reasonable range of the construction material physical properties, while that of the traditional model (8.15 × 10 −6 per • C) is again slightly smaller than the theoretical model, which indicates the reliability of our method.
Remote Sens. 2018, 10, x FOR PEER REVIEW 17 of 22 method and our method comprehensively consider all of the measurements at different temperatures to minimize the total errors, rather than only considering the measurement with no thermal dilation. Therefore, the measurements are not always exactly the same as the fitted values, and the slight differences in the no thermal dilation temperature are within a reasonable range. Based on the standard temperature, a temperature variation of 1 °C may lead to a change of 8.15 10 and 9.37 10 of the bridge's length, according to the traditional model and our model.
The linear thermal dilation coefficient estimated from our method (9.37 10 per °C) is within the reasonable range of the construction material physical properties, while that of the traditional model (8.15 10 per °C) is again slightly smaller than the theoretical model, which indicates the reliability of our method.  After modelling and separating the thermal dilation using our method, the trend deformation and thermal dilation of the TKH are calculated and are shown in Figure 17a,b, respectively. Clearly, the TKH is quite stable during the observation period, especially on the highway segment, and a method and our method comprehensively consider all of the measurements at different temperatures to minimize the total errors, rather than only considering the measurement with no thermal dilation. Therefore, the measurements are not always exactly the same as the fitted values, and the slight differences in the no thermal dilation temperature are within a reasonable range. Based on the standard temperature, a temperature variation of 1 °C may lead to a change of 8.15 10 and 9.37 10 of the bridge's length, according to the traditional model and our model.
The linear thermal dilation coefficient estimated from our method (9.37 10 per °C) is within the reasonable range of the construction material physical properties, while that of the traditional model (8.15 10 per °C) is again slightly smaller than the theoretical model, which indicates the reliability of our method.  After modelling and separating the thermal dilation using our method, the trend deformation and thermal dilation of the TKH are calculated and are shown in Figure 17a,b, respectively. Clearly, the TKH is quite stable during the observation period, especially on the highway segment, and a After modelling and separating the thermal dilation using our method, the trend deformation and thermal dilation of the TKH are calculated and are shown in Figure 17a,b, respectively. Clearly, the TKH is quite stable during the observation period, especially on the highway segment, and a slight trend deformation of less than 5 mm/year was observed on the Cheung Tsing Bridge. The temperature deformation parameters varied from −2.5 to 2.5 mm/ • C and propagated longitudinally on each sub-segment. The thermal dilation of the Cheung Tsing Bridge gradually increased from the middle to both ends with opposing directions. The thermal dilation stopped at the locations of the expansion joints and started to transmit again in the next segment because the expansion joints in the TKH hindered the thermal dilation transmission.
Remote Sens. 2018, 10, x FOR PEER REVIEW 18 of 22 slight trend deformation of less than 5 mm/year was observed on the Cheung Tsing Bridge. The temperature deformation parameters varied from −2.5 to 2.5 mm/°C and propagated longitudinally on each sub-segment. The thermal dilation of the Cheung Tsing Bridge gradually increased from the middle to both ends with opposing directions. The thermal dilation stopped at the locations of the expansion joints and started to transmit again in the next segment because the expansion joints in the TKH hindered the thermal dilation transmission.

Discussion of Thermal Dilation Mechanisms
To further investigate the thermal dilation evolution process, an in-depth analysis was developed by summarizing the thermal dilation results of different civil infrastructures types. Then, the common characteristics of civil infrastructure thermal dilations were concluded and are illustrated in Table 3, which provides valuable information for the engineers to conduct thermal dilation related research in the future.
First, according to the three case studies, we found that the transmission direction of thermal dilation is highly dependent on the specific structural shape of the target, and it tends to transmit along the direction of the longest side of the structures' geometric shape. Since the shape of the SWFC is very long in the vertical direction but is much shorter in the horizontal direction, the thermal dilation fringes mainly changed in the building's vertical direction. However, as far as the SNIEC and TKH are concerned, they are longer in the horizontal direction, but are much shorter in the vertical direction. Therefore, their thermal dilations are mainly accumulated along the horizontal direction, especially in the longitudinal directions. In general, the thermal dilations of vertical infrastructures, such as skyscrapers, tend to transmit along the vertical direction, while those of the horizontal infrastructures, such as low buildings and highways, are proven to propagate in the

Discussion of Thermal Dilation Mechanisms
To further investigate the thermal dilation evolution process, an in-depth analysis was developed by summarizing the thermal dilation results of different civil infrastructures types. Then, the common characteristics of civil infrastructure thermal dilations were concluded and are illustrated in Table 3, which provides valuable information for the engineers to conduct thermal dilation related research in the future. First, according to the three case studies, we found that the transmission direction of thermal dilation is highly dependent on the specific structural shape of the target, and it tends to transmit along the direction of the longest side of the structures' geometric shape. Since the shape of the SWFC is very long in the vertical direction but is much shorter in the horizontal direction, the thermal dilation fringes mainly changed in the building's vertical direction. However, as far as the SNIEC and TKH are concerned, they are longer in the horizontal direction, but are much shorter in the vertical direction. Therefore, their thermal dilations are mainly accumulated along the horizontal direction, especially in the longitudinal directions. In general, the thermal dilations of vertical infrastructures, such as skyscrapers, tend to transmit along the vertical direction, while those of the horizontal infrastructures, such as low buildings and highways, are proven to propagate in the horizontal direction.
Second, the spatial pattern of thermal dilation can provide valuable information regarding the static structural characteristics of infrastructures, such as the positions of structural key points. As for the SWFC, the temperature deformation parameter approximately equals zero at the bottom and increased gradually to the maximum value at the top, indicating that the bottom on the ground is probably the fixed point of the building and the top of the building is structurally free. In terms of the SNIEC, the temperature deformation parameter remains zero in the middle of the exhibition halls and it increased gradually to both sides with opposing values, which indicates that the fixed points are most likely located in the middle of the buildings, while both sides are structurally free. Regarding the TKH, the discontinuities of the temperature deformation parameter indicated the locations of the expansion joints along the highway, because these constraints hindered the accumulation of thermal dilation [18,27,28]. Therefore, the locations of structural key points along infrastructures can be approximately inferred from the estimated thermal dilation patterns.
Finally, the magnitude of thermal dilation is associated with the materials and sizes of infrastructures. Since the thermal dilation showed a highly linear correction with the temperature variation, it can be modelled while using a linear regression model. Based on the analysis of the thermal dilation model, two types of information would affect the thermal dilation magnitude. On the one hand, the thermal dilation coefficients of different material specifications are slightly different, which would affect the magnitude of the accumulated thermal dilation within a certain temperature range. On the other hand, the size of the structures also affects the magnitude of the accumulated thermal dilation. The thermal dilation can increase to a large value when there are no extra constraints on the structures, while it would be limited to a smaller value when extra constraints shorten the transmission size. Therefore, with the same temperature difference and similar building material, such as the SWFC and SNIEC, the accumulated thermal dilation of the SWFC is larger than that of SNIEC because the longer effective thermal dilation transmission distance would lead to a larger accumulation of thermal dilation in the structures.

Conclusions
Given the advantages of richer ground target information and a higher sensitivity to the small deformations in the high-resolution SAR images, a multi-temporal DInSAR approach, which estimates the civil infrastructure thermal dilations with an increasing number of PTs and high-resolution topography estimation based on TerraSAR-X images, was applied in this study. Both the coherent and incoherent information derived from the high-resolution SAR images, and the semantic information of the specific structures were integrated into the DInSAR method in order to refine the PTs identification. Moreover, the interferometric subsets with small temporal baselines and temperature differences were used for the high-resolution topography estimation of modern civil infrastructures. Then, a two-step analysis method consisting of a model analysis phase and a model generation phase was applied to estimate the thermal dilation. First, a pre-analysis was adopted to determine the transmission direction, spatial pattern, and temporal variation of thermal dilation. Then, both the traditional least squares regression analysis and the robust coherence-weighted least squares regression analysis proposed in this study were performed to quantitatively estimate the thermal dilation model. In all cases studies, the thermal dilation models that were estimated from the two methods were validated in terms of the estimated linear thermal dilation coefficients, which indicate the greater reliability of our method. Furthermore, the general thermal dilation characteristics of different civil infrastructure types were explored by analysing the potential relationship between the thermal dilation results and their specific structural characteristics. Our results demonstrate that the thermal dilation transmission direction is highly dependent on the specific structural shape of the target, and it tends to transmit along the longest side direction of the structures' geometric shape. Moreover, the static structural characteristics of infrastructures, such as the locations of structural key points, can be approximately inferred from the thermal dilation spatial patterns. Finally, the thermal dilation magnitude is associated with the materials and the sizes of the objects. The thermal dilation characteristics of different civil infrastructure types that are explored in this study can provide valuable information for better understanding the thermal dilation evolution process, and these findings are instructive for the future related work.
Author Contributions: All five authors contributed to this work. X.Q. processed the datasets, implemented the methodology, and finished the manuscript. L.Z. guided the results discussion and interpretation. X.D. and M.L. supervised the research and provided valuable suggestions for the revision. M.Y. collected and preprocessed the SAR images and acquisition temperatures.