Spatial Data Analysis for Deformation Monitoring of Bridge Structures

: Weather conditions and di ﬀ erent operational loads often cause changes in essential parts of engineering structures, and this a ﬀ ects the static and dynamic behavior and reliability of these structures. Therefore, geodetic monitoring is an integral part of the diagnosis of engineering structures and provides essential information about the current state (condition) of the structure. The development of measuring instruments enables deformation analyses of engineering structures using non-conventional surveying methods. Nowadays, one of the most e ﬀ ective techniques for spatial data collection is terrestrial laser scanning (TLS). TLS is frequently used for data acquisition in cases where three-dimensional (3D) data with high resolution is needed. Using suitable data processing, TLS can be used for static deformation analysis of the structure being monitored. For dynamic deformation measurements (structural health monitoring) of bridge structures, ground-based radar interferometry and accelerometers are often used for vibration mode determination using spectral analysis of frequencies. This paper describes experimental deformation monitoring of structures performed using TLS and ground-based radar interferometry. The procedure of measurement, the analysis of the acquired spatial data, and the results of deformation monitoring are explained and described.


Introduction
Bridge structures are integral parts of the transport infrastructure in the Slovak Republic, and their number in recent decades has increased. The modernization of the transport infrastructure has caused an increase in traffic intensity, which is also reflected in the increased operating load of the bridges. This causes changes in the spatial position and the shape of the structures, which affect their static and dynamic functions and reliability. Due to these facts, taking geodetic measurements is an integral part of bridge structure diagnosis; the geometry of the structures can be measured, along with their static and dynamic responses.
Because terrestrial laser scanning (TLS) is used by the wider community of surveyors and structural and civil engineers, it is used for different surveying jobs, e.g., [1][2][3][4][5][6][7]. The benefit of TLS over the more usual surveying methods is the three-dimensional (3D) data collection's efficiency. TLS enables non-contact measurement of the 3D coordinates of points lying on the surface of the scanned object. The state-of-the-art instruments are able to measure with a scan rate of up to 2 million p/s, which results in an expressive reduction of the time necessary for the measurements and an increases in the amount of information about the structures being measured. Using suitable spatial data analysis approaches, the resulting point cloud can be used for deformation analysis [8].

Deformation Analysis Using TLS
Laser scanning is a non-selective technique for 3D data collection; therefore, with the exception of special cases of repetition of measurements from a single position of the instrument, it is practically impossible to reproduce the identical position of measured points in point clouds from different epochs of deformation measurements. This fact makes it impossible to directly determine the displacements, therefore the allocation of identical parts of the object is performed in a separate step, which usually involves modeling parts of the point cloud [8]. The approaches for deformation analysis from terrestrial laser scanning results can be divided according to [15] into five categories: • Point-based approaches; • Point cloud-based approaches; • Surface-based approaches; • Geometry-based approaches; • Parameter-based approaches.
Point-based approaches are based on determination of displacements, which are based on the difference in the coordinates of repeatedly measured points in the clouds [16]. Such an approach can be used only if the scanner allows repeated measurement of points without changing the orientation of the instrument. In practice, this means that, from one position of the instrument, the surface of the monitored object is repeatedly scanned at certain time intervals using identical scanning parameters in each epoch while the stability of the instrument must be ensured during the whole monitoring process.
Point cloud-based approaches determine the displacements and deformations from the difference of two clouds. The relationship between clouds is determined by their transformation (registration) into a common coordinate system, most often using the ICP algorithm [17]. For example, according to [18], the clouds are divided into clusters after the registration process using hierarchical clustering methods, e.g., octant tree. The displacements between the clouds can then be determined in different ways. One way is to calculate the distance between the nearest neighbors. Another way is to calculate the average distance between corresponding clusters, which is calculated from the distances between identical pairs of points from both clouds, or by calculating the Hausdorff distance (the maximum distance from the set of the distances of nearest neighbors). The angles between the normal vectors of the planar surfaces approximating the points of the cloud in the corresponding clusters can also indicates deformations. A simple determination of the coordinate differences between the corresponding clusters' centroids also belongs to this category. The surface-based approaches are based on modeling point clouds by surfaces [19]. The displacements are determined by calculating the distances between the cloud points in the current epoch and the triangular network created in the basic epoch [8]. The second approach is to determine the distances between two surfaces in a defined regular grid of points [20].
Geometry-based approaches determine changes in the spatial position and orientation of the geometric primitives (planar surface, cylinder, sphere, etc.) modeled from point clouds. Parameter-based approaches are essentially an extension of geometry-based approaches [8]. In this case, displacements and deformations are determined by the statistical testing of changes in the characteristic parameters of geometric primitives in individual epochs [1].

Geometry-Based Approach for Deformation Analysis of Brgidge Strucutres
In most cases, the structural elements of bridge structures consist of regular geometric shapes and their combinations. This fact enables us to model the chosen part of the monitored structure by geometric primitives and therefore leads to the utilization of a geometry-based deformation analysis approach to determine the displacements of the monitored object. To improve the quality (accuracy) of the deformation analysis, the monitored parts (points) can be modelled by regression planes fitted to segmented parts of the point cloud. The advantage of regression planes is improved precision due to the reduced measurement noise, which, in the end, results in an improved deformation analysis.
In general, displacements of monitored points can be identified as coordinate differences between their coordinates in individual measurement epochs. For vertical displacements, this will be the difference between the heights (Z coordinates) of the monitored points. In the following section, an approach based on regression plane modelling is described. The heights of the monitored points are determined by modelling small planar surfaces using orthogonal regression. The monitored points are fixed in the horizontal plane (XY coordinates) in a local horizon, and only their heights are calculated (Z coordinates) ( Figure 1). The segmentation of the point cloud is done by RANSAC [21] from a pre-segmented subset of points using fencing boxes. The fencing boxes are defined around the fixed position of the monitored points by their dimensions along the axis of a local coordinate system of the monitored structure. The parameters of the planes are calculated by orthogonal regression from the general equation of a plane: where a, b, and c are the normal vector's parameters of the plane; X, Y, and Z are the coordinates of the point lying in the plane; and d is the scalar product of the normal vector of the plane and the position vector of any point of the plane. The parameters of normal vectors are calculated using singular value decomposition (SVD) [22]: where the design matrix A has the dimensions nx3 where n is the number of points included in the plane's parameter estimation. The column vectors of U nxn are normalized eigenvectors of matrix AA T . The column vectors of V 3×3 are normalized eigenvectors of A T A. The matrix Σ nx3 contains eigenvalues on the diagonals. Then the regression plane's normal vector is the column vector of V, which corresponds to the smallest eigenvalue from Σ.
reference points after the registration of point clouds into a common coordinate system. From the orthogonal distances of the subset of points from the plane, the standard deviation of the regression plane is calculated. The dispersion of the points around the plane reflects the noise of the point cloud, which is mainly affected by the rangefinder (distance measurement) of the instrument used. To minimize the effects of systematic errors, the measurements should be performed under equal conditions in each epoch (position of the scanner, temperature, etc.) if possible. The vertical displacements of the monitored points are calculated as the Z coordinate differences of these points (calculated by (4)) in each measurement epoch. Finally, the vertical displacements can be transformed into the normal direction to the object's surface based on the estimated normal vectors. The benefit of this procedure is that the monitored object itself defines the direction of the displacement and therefore it is not dependent on the scanner position or on the coordinate system axis orientation.

Structural Health Monitoring Using Ground-Based Radar
Ground-based radar interferometry is an innovative and rapidly evolving technique for the dynamic deformation analysis of engineering structures, even those of large investments such as dams and bridge structures [23][24][25]. The radar measurements use the stepped frequency continuous wave (SF-CW) approach in microwave bandwidth [26]. This approach allows determination of the monitored object's displacements in the radar's field of view in a radial direction (from the instrument to the object). The measurement is based on the transmission of a set of pulses created by a beam of frequency modulated electromagnetic signals (waves). To get radial resolution (range resolution), the instrument emits short pulses of beams. The relationship between the pulse duration and the radial resolution is defined by the formula: The design matrix has the form: where (X i − X 0 ), (Y i − Y 0 ), and (Z i − Z 0 ) are the reduced coordinates (reduced to the subset's centroid). The advantage of using SVD is that no initial parameters are needed for estimation. The above-mentioned procedure is performed by RANSAC, while the seed points are the points nearest to the defined positions of the monitored points.
The heights (Z coordinates) of the monitored points are calculated by projecting the predefined position (XY) of these points onto the estimated regression planes ( Figure 1) using the formula: The benefit of this approach is that the position of the monitored points does not change with the thermal expansion of the structure; the uncertainty caused by different scan densities (e.g., by different positions of the instrument) in the single epochs are also eliminated. Using uncertainty propagation law, the standard deviations of the results are calculated from the standard deviation of the vertical component of the registration error and the standard deviation of the regression plane. The vertical component is calculated as a quadratic mean of differences, ∆Z, between the identical reference points after the registration of point clouds into a common coordinate system. From the orthogonal distances of the subset of points from the plane, the standard deviation of the regression plane is calculated. The dispersion of the points around the plane reflects the noise of the point cloud, which is mainly affected by the rangefinder (distance measurement) of the instrument used. To minimize the effects of systematic errors, the measurements should be performed under equal conditions in each epoch (position of the scanner, temperature, etc.) if possible.
The vertical displacements of the monitored points are calculated as the Z coordinate differences of these points (calculated by (4)) in each measurement epoch. Finally, the vertical displacements can be transformed into the normal direction to the object's surface based on the estimated normal vectors.
The benefit of this procedure is that the monitored object itself defines the direction of the displacement and therefore it is not dependent on the scanner position or on the coordinate system axis orientation.

Structural Health Monitoring Using Ground-Based Radar
Ground-based radar interferometry is an innovative and rapidly evolving technique for the dynamic deformation analysis of engineering structures, even those of large investments such as dams and bridge structures [23][24][25]. The radar measurements use the stepped frequency continuous wave (SF-CW) approach in microwave bandwidth [26]. This approach allows determination of the monitored object's displacements in the radar's field of view in a radial direction (from the instrument to the object). The measurement is based on the transmission of a set of pulses created by a beam of frequency modulated electromagnetic signals (waves). To get radial resolution (range resolution), the instrument emits short pulses of beams. The relationship between the pulse duration and the radial resolution is defined by the formula: where c is the light's speed in vacuum (empty space) and τ is the time of flight (ToF) of the measuring signal. The frequency response is then determined by comparison of the emitted and received signal. By the application of inverse Fourier transformation (IFT), this response can be transformed from the frequency domain to the time domain. After this procedure, the radar creates a so-called range bin profile, which is a one-dimensional (1D) image representing the signal to noise ratio (SNR) at individual radial distances from the instrument ( Figure 2). In this range bin profile, the peaks of the reflected signal are detected. Each peak represents an average reflection from the whole range resolution area, e.g., when the range resolution is 0.5 m then the corresponding signal peak represents the average reflection from the parts of the structure monitored, which is in this 0.5 m range interval (radial resolution area).
Appl. Sci. 2020, 10, x FOR PEER REVIEW 5 of 13 where c is the light's speed in vacuum (empty space) and τ is the time of flight (ToF) of the measuring signal. The frequency response is then determined by comparison of the emitted and received signal. By the application of inverse Fourier transformation (IFT), this response can be transformed from the frequency domain to the time domain. After this procedure, the radar creates a so-called range bin profile, which is a one-dimensional (1D) image representing the signal to noise ratio (SNR) at individual radial distances from the instrument ( Figure 2). In this range bin profile, the peaks of the reflected signal are detected. Each peak represents an average reflection from the whole range resolution area, e.g., when the range resolution is 0.5 m then the corresponding signal peak represents the average reflection from the parts of the structure monitored, which is in this 0.5 m range interval (radial resolution area). After the range bin profile is created and the targets are identified, the changes in the position of the monitored parts are determined by differential interferometry (DI). The DI determines the radial displacements by comparison of the phase shift of the emitted and received measuring signal on the base of the formula: where λ is the wavelength of the signal and φ1 and φ2 are the phase shifts of the signal in first and subsequent epoch.
As the basic principle of the radar interferometry indicates, the displacements of the monitored structure are determined in the radial direction (from the instrument to the structure being monitored). In practice, it is almost impossible to ensure the direct deformation monitoring in a After the range bin profile is created and the targets are identified, the changes in the position of the monitored parts are determined by differential interferometry (DI). The DI determines the radial displacements by comparison of the phase shift of the emitted and received measuring signal on the base of the formula: where λ is the wavelength of the signal and ϕ 1 and ϕ 2 are the phase shifts of the signal in first and subsequent epoch.
As the basic principle of the radar interferometry indicates, the displacements of the monitored structure are determined in the radial direction (from the instrument to the structure being monitored). In practice, it is almost impossible to ensure the direct deformation monitoring in a required direction, e.g., vertical or horizontal. The solution can be the transformation of the displacement from the radial direction to the required one based on the known geometry of the measurements (relationship between the instrument and the object monitored); e.g., the transformation to the vertical direction is illustrated in Figure 3.

Spectral Analysis of the Data
Different methods of spectral analysis can be used when the identification of the bridge's vibration modes is required. Fourier transform (FT) is frequently applied for this purpose. It expresses the signal in the form of time series by continuously differentiable functions. The signal can then be transformed between time and frequency domains, and it can be represented as a continuous or discrete signal.
In practice, the signal represented by a finite number of measured values is analyzed by FT, known as discrete Fourier transformation (DFT). The mathematical definition of DFT depends on the data analyzed. When the goal is the analysis of the dynamic deformations of bridge structures, the fast Fourier transformation (FFT) is used in most cases. The FFT is described by the formula: where γx(k) represents the autocorrelation function and w(k) is the spectral window function [27].
An alternative to the approach described in the previous section is the so-called Welch method. Using the Welch method, overlapping segments of the signal are used for the determination of the spectral density. FFT is then executed using these overlapping segments. The approach ensures smoother periodograms and improved precision of the estimated frequencies, though the magnitude spectrum's resolution is decreased [28].
The estimation of the mutual spectral density of synchronized data can be performed by their spectral analysis. By the calculation of an average value of the normalized periodograms, the average normalized spectral density (ANSPD) can be computed. The ANPSD characterizes the spectral density of every time series processed in the data analysis. It allows us to create a complex picture about the static and dynamic behavior of the structure being monitored. Time synchronized measurement of the entire structure by a G-B radar gives the possibility to calculate the average periodograms for the whole structure (Figure 4). In practice, the signal represented by a finite number of measured values is analyzed by FT, known as discrete Fourier transformation (DFT). The mathematical definition of DFT depends on the data analyzed. When the goal is the analysis of the dynamic deformations of bridge structures, the fast Fourier transformation (FFT) is used in most cases. The FFT is described by the formula: where γ x (k) represents the autocorrelation function and w(k) is the spectral window function [27]. An alternative to the approach described in the previous section is the so-called Welch method. Using the Welch method, overlapping segments of the signal are used for the determination of the spectral density. FFT is then executed using these overlapping segments. The approach ensures smoother periodograms and improved precision of the estimated frequencies, though the magnitude spectrum's resolution is decreased [28]. The estimation of the mutual spectral density of synchronized data can be performed by their spectral analysis. By the calculation of an average value of the normalized periodograms, the average normalized spectral density (ANSPD) can be computed. The ANPSD characterizes the spectral density of every time series processed in the data analysis. It allows us to create a complex picture about the static and dynamic behavior of the structure being monitored. Time synchronized measurement of the entire structure by a G-B radar gives the possibility to calculate the average periodograms for the whole structure (Figure 4). spectral analysis. By the calculation of an average value of the normalized periodograms, the average normalized spectral density (ANSPD) can be computed. The ANPSD characterizes the spectral density of every time series processed in the data analysis. It allows us to create a complex picture about the static and dynamic behavior of the structure being monitored. Time synchronized measurement of the entire structure by a G-B radar gives the possibility to calculate the average periodograms for the whole structure (Figure 4).

The Liberty Bridge-A Case Study
The Liberty Bridge is part of the cycling route between the Bratislava district of Devínska Nová Ves (the Slovak Republic) and Schlosshof (Austria). It spans the river Morava and a stagnant pool of associated flows with a total length of 525.0 m [29]. The bridge is located in an inundation area on both sides in a protected floodplain forest ( Figure 5). The substructure consists of reinforced-concrete pillars in which the supports are anchored. The superstructure of the bridge over the river consists of three bridge sections with a total span of 180.0 m. The main part of the superstructure is a steel triangular truss beam with an orthotropic deck. The suspended bridge section has a radius of 376.35 m and is suspended on four pylons, which create dual-hinged rectangular frames. The pylons' diameter is 0.914 m, while the height of each is 17.7 m.

The Liberty Bridge-A Case Study
The Liberty Bridge is part of the cycling route between the Bratislava district of Devínska Nová Ves (the Slovak Republic) and Schlosshof (Austria). It spans the river Morava and a stagnant pool of associated flows with a total length of 525.0 m [29]. The bridge is located in an inundation area on both sides in a protected floodplain forest ( Figure 5). The substructure consists of reinforced-concrete pillars in which the supports are anchored. The superstructure of the bridge over the river consists of three bridge sections with a total span of 180.0 m. The main part of the superstructure is a steel triangular truss beam with an orthotropic deck. The suspended bridge section has a radius of 376.35 m and is suspended on four pylons, which create dual-hinged rectangular frames. The pylons' diameter is 0.914 m, while the height of each is 17.7 m.

Deformation Analysis Usign Terrestrial Laser Scanning
Monitoring using TLS was executed in three epochs of measurement using a Leica ScanStation2 instrument. The surface of the bottom of the main suspended bridge section was scanned from a single instrument's position. The instrument was set up on the longitudinal axis of the structure on the Slovakian riverbank to ensure that the whole structure's rear section could be scanned ( Figure 6). The minimum density of the scanned points on the structure was 3 mm. The next two epochs of the measurements were performed using the same Leica ScanStation2 with approximately the same conditions as in the first (initial) epoch.

Deformation Analysis Usign Terrestrial Laser Scanning
Monitoring using TLS was executed in three epochs of measurement using a Leica ScanStation2 instrument. The surface of the bottom of the main suspended bridge section was scanned from a single instrument's position. The instrument was set up on the longitudinal axis of the structure on the Slovakian riverbank to ensure that the whole structure's rear section could be scanned ( Figure 6). The minimum density of the scanned points on the structure was 3 mm. The next two epochs of the measurements were performed using the same Leica ScanStation2 with approximately the same conditions as in the first (initial) epoch.
instrument. The surface of the bottom of the main suspended bridge section was scanned from a single instrument's position. The instrument was set up on the longitudinal axis of the structure on the Slovakian riverbank to ensure that the whole structure's rear section could be scanned ( Figure 6). The minimum density of the scanned points on the structure was 3 mm. The next two epochs of the measurements were performed using the same Leica ScanStation2 with approximately the same conditions as in the first (initial) epoch.  The reference point network consisted of four control points (VB1-VB4), realized by Leica HDS type targets. Two of the reference points were the points of the reference network established for the construction of the bridge, stabilized with concrete pillars. Due to the vegetation and the fact that the whole structure is situated in a protected area of floodplain forest, where cutting out the vegetation (for lines of sight) is not allowed, two of the reference points were stabilized by metallic fasteners on the base of the pillars.
The scanned data was transformed to a local coordinate system, which was defined for the bridge for the purpose of geodetic monitoring with a local horizon. The reference points defined the coordinate system in each measurement epoch. The aim of the data processing was the identification and quantification of the bridge's main section in discrete points positioned on the bottom of the bridge structure. The vertical displacements were determined using the geometry-based approach of deformation analysis described in Section 2.1. During the data processing, square fences of 75 mm × 75 mm were defined around the fixed position of the monitored points (these created fencing boxes for the pre-segmentation of the points in each epoch). The monitored points were situated on the bottom of the transverse girders between the diagonal reinforcements of the supporting girders on both sides of the bridge; there were 23 points on each side (left and right), giving 46 in total.
The deformation analysis shows the changed position of every monitored point, except for the points on the ends of the bridge (Figure 7). The bridge is anchored to the pillars on both ends, which means there is no displacement is these parts. The standard deviation of the displacements, calculated using the uncertainty propagation law, varied from 1.3 mm to 1.8 mm. The vertical displacements increasing from the ends to the center of the bridge where reach negative values at the level of −13 mm (left side) and −10 mm (right side). The displacements were caused by the structure's temperature difference between the monitoring epochs and by the weight of the layer of fallen snow (approx. 100 mm).
Appl. Sci. 2020, 10, 8731 9 of 14 means there is no displacement is these parts. The standard deviation of the displacements, calculated using the uncertainty propagation law, varied from 1.3 mm to 1.8 mm. The vertical displacements increasing from the ends to the center of the bridge where reach negative values at the level of −13 mm (left side) and −10 mm (right side). The displacements were caused by the structure's temperature difference between the monitoring epochs and by the weight of the layer of fallen snow (approx. 100 mm).

Deformation Analysis Usign Ground-Based Radar Interferometry
The dynamic measurements were performed during different types of structural load, which were proposed based on the finite element method (FEM) model of the bridge using the structure's technical documentation [29]. Four loading epochs were defined as follows: 1st-without loading; 2nd-one person walking; 3rd-one person running; 4th-one person jumping in the middle of the bridge's span. Each epoch was performed in three phases. The first phase was the no load state (but the data was recorded). The second phase represented the load (4 types) and the last phase was again

Deformation Analysis Usign Ground-Based Radar Interferometry
The dynamic measurements were performed during different types of structural load, which were proposed based on the finite element method (FEM) model of the bridge using the structure's technical documentation [29]. Four loading epochs were defined as follows: 1st-without loading; 2nd-one person walking; 3rd-one person running; 4th-one person jumping in the middle of the bridge's span. Each epoch was performed in three phases. The first phase was the no load state (but the data was recorded). The second phase represented the load (4 types) and the last phase was again the no load state with the aim of determining the structure's damping. The individual loading epochs each took approximately 60 s.
The main suspended bridge section was measured from a single position of the interferometric radar IDS IBIS-S (Figure 8). Because the FEM model of the structure assumed significant frequencies of the structure's vibration on the level of 10 Hz, the data was registered with a frequency of 100 Hz to ensure higher relevance of the results (relative displacements), although it is sufficient to measure with a frequency that is twice as high as the accepted frequencies for their proper identification. Despite the fact that the G-B radar measurement allows monitoring of the entire structure, the measurement and data processing was focused specifically on two points of the structure. The first point was at the center of the suspended structure and the second one was positioned at the anchorage of the suspension cable.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 9 of 13 the no load state with the aim of determining the structure's damping. The individual loading epochs each took approximately 60 s. The main suspended bridge section was measured from a single position of the interferometric radar IDS IBIS-S (Figure 8). Because the FEM model of the structure assumed significant frequencies of the structure's vibration on the level of 10 Hz, the data was registered with a frequency of 100 Hz to ensure higher relevance of the results (relative displacements), although it is sufficient to measure with a frequency that is twice as high as the accepted frequencies for their proper identification. Despite the fact that the G-B radar measurement allows monitoring of the entire structure, the measurement and data processing was focused specifically on two points of the structure. The first point was at the center of the suspended structure and the second one was positioned at the anchorage of the suspension cable. The data analysis was realized in several steps, and its aim was the determination of the amplitude and the frequency of the vibration of the structure's selected parts. The processing was arranged into following steps: 1. Radial displacements' transformation to vertical; 2. Identification of the monitored parts (P01, P02); 3. Modal analysis.
The displacements were transformed from the radial to the vertical direction on the basis of the position and orientation of the radar in relation to the structure being monitored. This geometry was The data analysis was realized in several steps, and its aim was the determination of the amplitude and the frequency of the vibration of the structure's selected parts. The processing was arranged into following steps:
The displacements were transformed from the radial to the vertical direction on the basis of the position and orientation of the radar in relation to the structure being monitored. This geometry was obtained from the laser scanning data. The second step was the identification of the peaks of the signal that corresponded to the chosen monitored parts of the structure from the range bin profile (Figure 9). The figure shows the estimated signal to ratio (SNR) of the signal's reflection, depending on the structure's range. The last step was the determination of the relative displacements, auto-spectral analysis of the radar data using the approaches described in the chapter 3.1, and the comparison of the results with the FEM model of the structure. The data analysis was realized in several steps, and its aim was the determination of the amplitude and the frequency of the vibration of the structure's selected parts. The processing was arranged into following steps: 1. Radial displacements' transformation to vertical; 2. Identification of the monitored parts (P01, P02); 3. Modal analysis.
The displacements were transformed from the radial to the vertical direction on the basis of the position and orientation of the radar in relation to the structure being monitored. This geometry was obtained from the laser scanning data. The second step was the identification of the peaks of the signal that corresponded to the chosen monitored parts of the structure from the range bin profile (Figure 9). The figure shows the estimated signal to ratio (SNR) of the signal's reflection, depending on the structure's range. The last step was the determination of the relative displacements, autospectral analysis of the radar data using the approaches described in the chapter 3.1, and the comparison of the results with the FEM model of the structure.  The vertical displacements are affected by a only very small extent by the pedestrian walking ( Table 1). The fast movement (running) of the pedestrian has almost the same effect on the maximum vertical displacements as that of the walking pedestrian. One person jumping in the center of the span results in the maximum vertical displacement, with values of 2.55 mm (P01) and 4.63 mm (P02). During the measurement epoch without any load on the structure, a frequency of deformation at the level of 1.53 Hz was determined ( Table 2). Signals from both monitored points were of low coherence with a relatively high phase delay of around 65.0 • . The dominant frequency of vibration determined during the epoch with one person walking corresponds to the structure's 22nd vibration mode (from the FEM model). The estimated 2.01 Hz corresponds to the frequency of steps of standard walking. A phase shift of 23.0 • is caused by a short delay in the structure's response at points P01 and P02. The 3rd loading epoch shows the vibration of the structure caused by the running person to be in the range of 1.53-4.56 Hz, which corresponds to the 2nd, 23rd, and 59th vibration modes. The 4th load epoch was executed when one person was jumping synchronously at the middle of the bridge's main span. The vibrations with a frequency of 1.82 Hz with minimum phase shift are caused purely by a person jumping. The identified frequency is close to the 22nd vibration mode defined by FEM. To verify the results, measurements using accelerometers were also performed. Two HBM B12/200 one-axial accelerometers were installed at the measured points, P01 and P02. The accelerometers were levelled in a vertical position using Zeiss tribraches, as shown in Figure 10, to ensure the measurement of the acceleration in the vertical direction. The accelerometers used were inductive sensors with a relative error of up to ±2%. The frequency of the data registration was 100 Hz, as in the case of G-B radar, and it was performed using a HBM Spider 8 A/D transducer and Catman Easy software.
coherence with a relatively high phase delay of around 65.0°. The dominant frequency of vibration determined during the epoch with one person walking corresponds to the structure's 22nd vibration mode (from the FEM model). The estimated 2.01 Hz corresponds to the frequency of steps of standard walking. A phase shift of 23.0° is caused by a short delay in the structure's response at points P01 and P02. The 3rd loading epoch shows the vibration of the structure caused by the running person to be in the range of 1.53-4.56 Hz, which corresponds to the 2nd, 23rd, and 59th vibration modes. The 4th load epoch was executed when one person was jumping synchronously at the middle of the bridge's main span. The vibrations with a frequency of 1.82 Hz with minimum phase shift are caused purely by a person jumping. The identified frequency is close to the 22nd vibration mode defined by FEM. To verify the results, measurements using accelerometers were also performed. Two HBM B12/200 one-axial accelerometers were installed at the measured points, P01 and P02. The accelerometers were levelled in a vertical position using Zeiss tribraches, as shown in Figure 10, to ensure the measurement of the acceleration in the vertical direction. The accelerometers used were inductive sensors with a relative error of up to ±2%. The frequency of the data registration was 100 Hz, as in the case of G-B radar, and it was performed using a HBM Spider 8 A/D transducer and Catman Easy software. The relative displacements from the data measured were calculated by double integration. The drift of the accelerometers and the errors caused by the integration were minimized by a Butterworth high-pass filter with a cut-off frequency of 0.5 Hz. The reduction of the spectrum's magnitude by the The relative displacements from the data measured were calculated by double integration. The drift of the accelerometers and the errors caused by the integration were minimized by a Butterworth high-pass filter with a cut-off frequency of 0.5 Hz. The reduction of the spectrum's magnitude by the filter (0.7% at 1 Hz) had no significant influence on the determined displacements. Table 3 shows the results of the measurements by the accelerometers. The values of the vertical displacements are very close to the displacements obtained by the G-B radar measurements in all four epochs of the experiment ( Table 3). The maximum vertical displacements were indicated by one person jumping at the middle of the superstructure, with values of 2.34 mm (P01) and 4.82 mm (P02).
In the 1st measurement epoch (without any load), no significant frequencies were identified from the measured data (Table 4). During the epoch with one person walking, the structure's 22nd and 44th vibration modes with frequencies 2.05 Hz and 3.99 Hz were determined. The results of the 3rd epoch (running person), 2.90 Hz and 4.59 Hz, correspond to the 23rd and 59th vibration modes and are in accordance with the results of the G-B radar measurements.

Conclusions
This paper deals with the deformation monitoring of bridge structures by non-conventional surveying methods: terrestrial laser scanning and ground-based radar interferometry. The basic principles of deformation monitoring using the mentioned methods are described with a special focus on the data analysis approaches. The paper describes the experimental deformation monitoring of the Liberty Bridge (Bratislava, Slovak Republic).
The results of the case study presented in this paper show that both methods described are highly effective for the deformation analysis of engineering structures. The accuracy of 3D coordinate measurement of points by state-of the-art laser scanners is in the range of a few millimeters. The precision can be increased using the data processing approach described in the paper. In addition to the described geometry-based approach for deformation analysis, information concerning the behavior of the reflected laser beam is needed, as shown in [30].
The G-B radar interferometry is an innovative contactless technique that allows the monitoring of the entire structure's deformation synchronously with high frequency. The accuracy of the results depends on the system configuration and on the quality of the reflected signal. The results obtained by G-B radar interferometry were verified by accelerometer measurements. Comparison of the result shows close compliance between the results of both methods. The measured deformations allow for the prediction of possible structural failure, often before an accident occurs.
The results of the measurements performed and described in the paper confirmed the assumptions about the static and dynamic responses of the structure being monitored, as defined in the technical documentation of the bridge construction project.