Methodology for selecting measurement points that optimize information gain for model updating

Information collected through sensor measurements has the potential to improve knowledge of complex-system behavior, leading to better decisions related to system management. In this situation, and particularly when using digital twins, the quality of sensor data determines the improvement that sensors have on decision-making. The choice of the monitoring system, including sensor types and their configuration, is typically made using engineering judgement alone. As the price of sensor devices is usually low, large sensor networks have been implemented. As sensors are often used to monitor at high frequencies over long periods, very large data sets are collected. However, model predictions of system behavior are often influenced by only a few parameters. Informative data sets are thus difficult to extract as they are often hidden amid redundant and other types of irrelevant data when updating key parameter values. This study presents a methodology for selecting informative measurements within large data sets for a given model-updating task. By selecting the smallest set that maximizes the information gain, data sets can be significantly refined, leading to increased data-interpretation efficiency. Results of an excavation case study show that the information gains with refined measurement sets that are much smaller than the entire data set are better than using the data set prior to refinement for the same probability of identification, while the computational time of model updating is significantly reduced. This methodology thus supports engineers for significant data filtering to improve model-updating performance.


Introduction
The Internet of Things (IoT) is a system that consists of connecting devices, objects, and people through the exchange of data and information [1].It is one of the key aspects guiding the development of responsive cities around the world.However, there are challenges.When scaling up the IoT concept to a city level, "Things" (physical objects) become much more complicated than simple household items [2].For example, when performance degrades, replacement is usually neither an easy option nor desirable for city elements.
Monitoring complex infrastructure systems, such as bridges, requires the recognition of several special characteristics.Effects, rather than causes, are generally measured.Data-only methods cannot support decisions related to activities such as retrofitting and replacement avoidance.Most measurement interpretation tasks require model-based diagnosis methodologies to infer causes from effects and then predict future behavior (prognosis), particularly when predictions involve extrapolation [3].Additionally, uncertainty magnitudes and bias are much larger as city systems are strongly affected by boundary conditions such as the environment.At best, behavior models produce approximate predictions [4].
Eventually, when infrastructure "Things" are upgraded due to, for example, a change in demand, replacement should be the last option [5].More sustainable options usually involve retrofitting and improvement rather than replacement [6,7].These actions need to be well-designed to control costs and minimize the impact on users.
The availability of inexpensive sensing devices [8,9] and computing tools [10] has made it feasible to monitor civil engineering systems.Many bridges are monitored under a structural health monitoring framework [11][12][13].
A particular characteristic of civil engineering systems, compared with other engineering fields, is that civil structures are often designed using conservative, rather than accurate, behavior models.Also, construction practices often reflect the assumption that design specifications are minimum requirements.Most civil infrastructure thus has reserve capacity that is greater than the reserve which was intended by safety factors [14].Decisions may also be conservative when asset management is made without quantifying this reserve capacity, leading to unsustainable and uneconomical maintenance activities.
Monitoring has the potential to improve understanding of structural behavior if it is followed by appropriate data interpretation methodologies.This process is called structural identification.Much research has been carried out to develop model-based data interpretation methodologies, see for example [15][16][17] among many others.Common data interpretation methodologies are residual minimization (RM) [18,19], Bayesian model updating (BMU) [20][21][22], and more recently, a special implementation of Bayesian model updating called error domain model falsification (EDMF) [23,24].EDMF is an easy-to-use methodology that has been developed to be compatible with typically available engineering knowledge as well as standard civil-structure assessment concepts [25,26].These methodologies mostly differ in the criteria used for model updating and in the assumptions related to the quantification of uncertainties.Since each complex system has a unique combination of geometry and material properties as well as specific monitoring goals, the most appropriate data interpretation methodology is case-dependent.A methodology map to provide guidance for the appropriate use of these methodologies has been recently proposed by [27].
The performance of structural identification directly depends on the choice of sensor configuration [28,29].To maximize the information gain with a minimum number of sensors, researchers have developed sensor placement algorithms such as those proposed by [30][31][32], among others.Since the general task has exponential computational complexity with respect to the number of sensor types and locations, most researchers have used greedy algorithms to reduce the computational effort [33][34][35].
Several studies have focused on determining the best objective function for optimal sensor placement.Multiple approaches have been proposed that involve either minimizing the information entropy in either posterior model parameter distributions [36][37][38] or maximizing information joint entropy in multiple-model predictions have been proposed [39,40].Most approaches provide ranking of potential sensor locations based on their ability to identify either bridge parameter values or structural damage [41,42].In addition to structural monitoring, optimal sensor placement studies have involved wind predictions [43], wind turbine condition monitoring [44], and steel pipeline damage detection [45], among others.
Recently, a novel strategy called the hierarchical algorithm has been introduced that explicitly accounts for the shared information within sensor data sets through a joint entropy metric [39,40].This methodology has then been adapted to structural identification contexts such as multiple static tests and dynamic excitations [46] and multiple objective criteria [47], and has been validated using field measurements [48].
Evaluation of these sensor placement methodologies has shown that most sensor locations provide redundant information [49].Measurement system design (also called optimal sensor placement when only one type of sensor is used) is usually carried out by engineers using only qualitative rules of thumb and estimations of signal-to-noise ratios, leading to sensor networks involving unnecessarily large numbers of sensors.As sensor placement algorithms are used for model updating, generation of model classes prior to monitoring as well as defining monitoring goals are necessary.In practice, such activities often take place only after monitoring is carried out.This is unfortunate since these methodologies can be used to filter data sets, thus improving data interpretation.In such cases, algorithms for optimal sensor placement should be adapted for the task of measurement point selection, which refers to the selection of the most informative data among previously collected data sets from multiple sensors.A measurement point is the data collected by a sensor at a given timestep.Each sensor thus collects multiple measurements throughout the duration of the measurement event.
The hierarchical algorithm has already been used to reduce sets of measurements in two applications: an excavation case study [49] and a wheel-flat identification [50].The aim was to enhance data interpretation by either reducing the computation time or improving the precision of parameter value identification.Efficient methodologies for measurement point selection should involve four stages to avoid removing informative data points.Stages include model generation, data cleansing, data-point selection, and validation.Previous studies [49,50] proposed stepwise processes for measurement point selection and did not account for the influence of the number of measurements on the selection process, leading to non-optimal data sets.
This paper contains a proposal for a comprehensive methodology for the last three stages mentioned above to improve structural identification.In this methodology, a strategy is proposed for data cleansing that involves both an outlier detection methodology and filtering of measurement points that provide redundant information.A result validation scheme is also included involving an iterative process to define the optimal number of measurements.
The paper is organized as follows.The error domain model falsification is presented in the background section.The following section presents the methodology for measurement point selection.Then, the methodology is applied to a full-scale case study.

Background-error domain model falsification
Error domain model falsification (EDMF) is a recently developed probabilistic methodology for data interpretation that has been introduced by [23].The particularity of this methodology is its robustness of the model-updating accuracy to variations in correlations between measurement points for most forms of uncertainty.The methodology involves several steps.First, the model class is chosen.A model class involves the creation of a parametric model through selection of its most critical behavior-sensitive characteristics, such as material properties, geometry, boundary conditions, and excitations, as well as the quantification of nonparametric model uncertainties ( U i,g ) and measurement errors ( U i,y ).Additionally, possible ranges of model parameters that can be identified during monitoring are defined.Then, a set of model instances is generated, where each instance represents a unique combination of model parameter values.The output of EDMF is a set of plausible model instances (i.e., plausible combinations of parameter values) among the initial population through comparisons of model instance predictions with field measurements.
For a measurement point i , the predictions of a model instance g i (Θ) are generated by assigning a unique set of parameter values Θ .The true structural response R i (unknown in practice) is linked to the field measurements y i and model prediction g i (Θ) using Eq. (1),where n y is the number of measurement points.
Distributions of model U i,g and measurement U i,y uncer- tainties are combined in a unique source of uncertainty U i,c (1) using Monte Carlo sampling.Equation ( 1) is then transformed in Eq. (2) as recommended by [51].
Plausible instances are selected through the falsification of instances that have residuals larger than defined threshold bounds.These thresholds are calculated using a probability target of identification ϕ on the combined uncertainty distribution U i,c .This target ϕ is usually set to 95% [23].When multiple measurements are considered, the Šidák correction [52] is used to maintain a constant level of confidence.Then, threshold bounds, t i,low and t i,high are evaluated (Eq. 3) for measurement point i.These bounds are calculated to be the shortest intervals in the probability density function of combined uncertainties f U i u i .
The candidate model set (CMS),Ω �� , contains unfalsified model instances defined using Eq. ( 4) from the initial model set Ω .These instances are set to be equally likely since lit- tle reliable information is usually available to describe the combined uncertainty distribution [51].They are assigned an equal probability, while falsified model instances are assigned a null probability.
It may happen that all model instances are falsified.In such situations, model predictions are incompatible with measurements given the estimations of nonparametric uncertainty sources.This result means that the current model does not adequately reflect the true structural behavior.Using EDMF thus leads to the re-evaluation of assumptions related to the choice of model class and uncertainty quantification [53].This situation highlights an important advantage of EDMF compared with other structural identification approaches, such as traditional implementations of BMU, particularly when there are few measurement points.

Overview
In this section, the methodology for selecting measurement points for data interpretation is presented.This methodology aims to improve the performance of data interpretation by removing measurements that do not provide information.By including only informative (i.e., useful) data, model updating requires less computational time than required for the entire set.A potential increase in precision is linked to the target probability that is corrected (such as using the Šidák correction [52]) to remain constant, independent of the number of measurements.
Figure 1 presents the framework of the methodology for measurement point selection.This methodology is subdivided into two sides, the model side and the measurement side.Inputs on the model side generate a numerical model (such as a finite-element model) and define a model class.
A model class includes the most relevant parameters (i.e., parameters that have the largest influence on model predictions for a given test condition based on a sensitivity analysis).On the measurement side, the inputs are the selection of sensor types and locations.
The first stage involves the collection of data under test conditions, constituting the initial data set.On the model side, the model class is used to generate model instances with a unique combination of model parameter values that Fig. 1 Flowchart of the methodology for measurement point selection.Relevant paper section numbers are provided for parts of the flowchart are within ranges.These parameter ranges are defined using engineering knowledge.These model instances are used to obtain predictions for each measurement collected.On the measurement side, field measurements are collected.
In the second stage, the initial data set is filtered using two model-based methodologies.First, outliers are removed.Measurements are defined as outliers if the value is outside the range of predictions, including nonparametric uncertainties.Then, non-informative measurements are removed.Non-informative measurements are defined using information entropy that evaluates the variability of model instance predictions for each measurement.On the measurement side, field measurements are filtered.
The third stage of the methodology involves selecting an optimal set of measurements using the hierarchical algorithm, following an iterative process.This algorithm evaluates the total information gain of a set of measurements, and the smallest set of measurements that maximizes joint entropy is chosen, and the remaining measurement points are discarded.Field measurements are then selected based on hierarchical-algorithm results.
In Stage 4, the model updating is made using the optimal set of measurements, reducing model parameter value ranges.The final stage (Stage 5) involves a validation procedure to ensure that the model updating is enhanced by comparing information gains using the entire data set and the selected measurement set.Using a smaller set of data, the precision of the model parameter value identification must be either equal or better in order to validate the results of the methodology.

Model class and model instance predictions
In this stage, the numerical model built for the inverse problem is instantiated in multiple model instances.Each instance is a distinct numerical model which has a unique combination of model parameter values.Each model instance has unique prediction values (Fig. 2).Then, the numerical model is used to generate model instance predictions for each measurement collected.These model instance predictions are the main input required for the measurement point selection in the next stages.

Outlier detection
This stage is divided into two steps.The first step involves using an outlier detection methodology to remove inconsistent data.When continuous data sets are involved, traditional outlier detection methods are designed to detect anomalies in measurements, see for example [54,55] among others.When traditional methods are appropriate for the case study, they can also be used in this step.
Traditional outlier detection methods are not suitable for examining data sets that consist of non-time-dependent measurements, such as in some civil engineering situations.Also, conventional outlier detection methodologies based on statistical methods and signal processing are not suitable in the context of sparse measurements due to the relatively low number of measurements.Therefore, threshold bound solutions, as proposed below, are recommended [56].
At each measurement point, a distribution is generated by combining the model instance prediction distribution with the combined uncertainty distribution using Monte Carlo sampling.This new distribution covers expected measurement values for a given current model class.A measurement point is judged to be an outlier if the measured value lies outside this prediction range (Fig. 3a).When an outlier is detected, it is then removed from the data set for the next stages.
An important assumption of the outlier detection methodology is that the model class that has been selected is correct.If a large number of the measurements (i.e., most measurements on multiple sensors) does not lie in the distribution of model instance predictions, including uncertainties, it is possible that the selected model class is not appropriate.In such situations, the model class must be modified, either through changing the model parameters included in the set or correcting the nonparametric uncertainty estimations, among other measures [57].
When it is established that a sensor has provided several outlier measurements (i.e., at least 20% of outliers), all data collected by this device is excluded, even though some of the data falls within the range of predictions.For example, this situation may happen either when the device has malfunctioned or if it has not been properly installed.

Non-informative measurement point filtering
The second stage involves removing data points that do not provide information.In computer science, information gain is defined as a reduction in uncertainty due to the new availability of data.In this methodology, the information gain is evaluated by the ability of a measurement set to discriminate between model instances based on their predictions and uncertainties.
To evaluate this variability, a model instance prediction histogram is generated at each measurement point i.
The range of model instance predictions is subdivided into N I,i intervals where the interval width is given by the dif- ference between threshold bounds taken from combined uncertainty U i,c (Eq.2).In the present study, this definition is adapted to account for a unique uncertainty distribution for each model instance.Details are provided in Sect.4.3.
The probability that the model instance prediction g i,j falls inside the j th interval is equal to P � g i,j � = m i,j ∕ ∑N I,i j=1 m i,j with m i,j is the number of model instances falling inside this specific interval.
The information entropy H g i of a measurement point i is evaluated using Eq. ( 5).As the information entropy measures the disorder in predictions, measurement points associated with large information entropy values represent a high potential for model instance discrimination.In other words, when the value of information entropy is larger than zero, this measurement may help reduce parametric uncertainties and be informative.
By evaluating measurement point information entropy, non-informative data are removed from the set (Fig. 3b).This filtering removes data points where model instance predictions have no variabilities compared with the uncertainties that do not provide additional information.This filtering helps increase the speed of the measurement point selection in the next stage.This step has two goals.First, it reduces the computational time of the measurement point selection as these optimization algorithms, as its computational complexity is O(n 2 ) with respect to the number of measurement points.Secondly, it reveals data points that are not useful, thus improving the understanding of the information collected by the data set.

Measurement point selection
An optimization algorithm is used to rationally choose a minimum of measurement points that maximize the information gain.The aim is to find the measurement set that will minimize parameter value ranges after monitoring.This strategy involves two principal components: an objective function to evaluate the sets of measurements and an optimization scheme.

Assessing information gain using joint entropy
When systems are monitored, such as bridges, measurements at nearby sensor locations are typically correlated.Therefore, selecting measurement points based on their information entropy values alone leads to measurement sets that provide significant redundant information [40].To account for the mutual information between measurement points, joint entropy has been introduced as a new objective function for sensor placement [39].Joint entropy H g i,i+1 quantifies the information entropy between sets of predictions at measurement points.For a set of two measurement points, joint entropy is calculated following Eq.( 6), where P g i,j , g i+1,k is the joint probability that model instance predictions fall inside the j th interval at measurement point i and the k th interval at measurement i + 1 .In this equation,k ∈ 1, … , N I,i+1 and N I,i+1 is the maximum number of prediction intervals at the i + 1 measurement and i + 1 ∈ 1, … , n y with the number of measurement points n y.
(5) H g i = − ∑ N I,i j=1 P g i,j log 2 P g i,j Due to the redundancy in information gain between measurement points, the joint entropy is less than or equal to the sum of individual information entropies at measurements i and i + 1 .Equation ( 6) can be changed to Eq. (7), where I g i,i+1 is the shared information between measurement points i and i + 1.
For a given number of measurements, the set of measurements associated with the largest joint entropy will provide the largest reduction in model parameter values.Therefore, joint entropy is used as the objective function of the optimization scheme.

Optimal set of measurement points using the hierarchical algorithm
Due to a large number of possible combinations of measurement points, a greedy search algorithm is recommended to reduce the computational time, similarly to other sensor placement methodologies [36].The hierarchical algorithm is a greedy search algorithm that efficiently evaluates joint entropy between sensor locations [39].This algorithm is used for ranking measurement points based on the joint entropy of sets of measurements.The ranking of measurement points is then used to select the appropriate number of measurements in the data set (Fig. 4).As long as the joint entropy is increasing, this means that additional measurements will increase the information gain (i.e., reducing model parameter ranges) after data interpretation.When the maximum joint entropy value is reached, this means that the remaining measurements are not helpful and should thus be removed from the data sets.The first set of selected measurements is defined as the minimum number of measurements to reach the maximum joint entropy value.( 6) Fig. 4 Measurement point selection using the hierarchical algorithm 123 An important assumption of the hierarchical algorithm is that histogram intervals are calculated based on the number of measurements in the data sets.Intervals taken in the evaluation of the joint entropy value is directly related to the threshold bounds defined in EDMF (Sect.2).These threshold bounds are corrected using the Šidák correction [52] which is used to maintain a constant level of probability with respect to the number of measurements.Thus, the number of measurements slightly influences the joint entropy values.This means that joint entropy values depend on this number (Eq. 6).When no results of the data filtering exist, such as during the first run, this number is equal to the initial size of measurement sets.An iterative process is then necessary until the initial assumption of the size of the selected measurement sets is correct.In each iteration, the hierarchical algorithm is run with the number of measurements selected in the previous iteration.Measurements selected in the last iteration constitute the data sets used for model updating.

Data interpretation
Once measurement points have been selected, these data are then used for model updating.Several methodologies exist.In this methodology, EDMF is recommended as it is an easy-to-use method that provides robust results in the presence of correlated and systematic uncertainties (Fig. 5).In this figure, the process of model instance falsification using EDMF is shown.Model instances in red are falsified as their predictions significantly differ from the measured values.All non-falsified model instances are included in the CMS.Nevertheless, the methodology for measurement point selection presented in this paper is also compatible with other data interpretation methodologies such as traditionally implemented BMU and residual minimization.

Validation
The last step of the measurement point selection methodology involves a validation procedure for the point selection.Such validation helps guarantee that the refinement of the data sets does not affect the quality of model updating.It is made using the following two checks: 1. Falsification is performed with each point that has a zero entropy value (Fig. 3).These measurements must not alter the falsification result.2. Falsification is performed with all data.Using reduced data sets, the precision of the data interpretation results should be equal to or better than using entire data sets.
If both criteria are satisfied, the measurement point selection is validated.This means that the selection of measurement points enhances the model-updating process since the computational time of the data interpretation method is reduced, and the precision of parameter value identification is equal to or higher than using all data sets.When these criteria are not satisfied, it means that either outliers have not been detected or that the model class selection was not appropriate.In such situations, an iteration of the entire methodology is required.
To compare data interpretation results, the falsification performance metric is introduced in Eq. ( 8) [48].This metric is bounded between 0 and 1.When the falsification value is close to 1, information gain is significant, while a value of 0 shows that the measurements did not increase knowledge.

Excavation numerical modeling
The case study is the excavation site of the Bukit Panjang Hawker Centre in Singapore.In this application, monitoring outcomes are used to improve the extrapolation of the retaining wall at subsequent excavation phases.Following (8) Falsif icationperformance = 1 − Sizeof thecandidatemodelset Sizeof theinitialmodelset code requirements in Singapore, several inclinometers must be placed on the excavation site.Additionally, each sensor measures the wall-deflection profile at 15 heights during the four excavation phases.In this study, measurements involve 10-15 individual measurements made by the inclinometers during the four excavation phases.The 10-m deep excavation is approximately 60 by 40 m large and is situated on the Bukit Timah Granite formation (Fig. 6) [58].Six boreholes were drilled on the excavation site.The 3-m top layer consists of sandy silt and man-made backfill materials.This is underlain by a 10-to 13-m-thick layer of sandy silt, followed by the granitic rock layer at approximately 15 m below the ground surface.On the east side of the excavation site, a 5-m-thick layer of coarse sand is found between the sandy silt and the granitic rock.

Excavation numerical modeling
Two-dimensional numerical models are traditionally used to predict retaining wall behavior.However, the prediction accuracy is affected by 3D effects such as corner constraints [59].To improve the quality of the model predictions, 2D and 3D numerical models are built on this case study [60].The 2D model (Fig. 7a) is used to generate model instance predictions, while the 3D effects are quantified using interpolation of the prediction discrepancies between both models.Inclinometer locations, named S1 to S10, are displayed on the 3D model (Fig. 7b).In both models, soil layers are described using the hardening soil with small strain stiffness model [61], while the rock layer is depicted using the Hoek-Brown model [62].The initial water table is 2 m below the ground level.
To retain earth, the support system includes diaphragm walls, soldier pile walls, toe pins, and two layers of steel struts and waler beams.The 800-mm-thick diaphragm walls are modeled as elastic plate members with reduced lateral stiffness due to the presence of construction joints.Toe pins and soldier pile walls are modeled as elastic plate members, while struts and waler beams are modeled as node-to-node anchors and beam elements, respectively.Soil-wall interactions are then modeled using an interface element without thickness.More information about the excavation modeling can be found in [49,60].
The construction sequence modeled in the finite element analysis involves six phases that are presented in Table 1.The aim of phase 0A is to generate initial ground stresses.The diaphragm wall is "wished-in-place" in phase 0B, and installation effects are assumed to be negligible.Fully coupled flow deformation calculations are performed to account for the combined effects of groundwater flow and time-dependent consolidation.Phases 1 to 4 correspond to the four-stage analysis, where model predictions are compared with field measurements for model updating.
Based on a composite scale sensitivity analysis [63], four model parameters are identified as having the most influence on model predictions, and they are presented in Table 2. Parameter bounds are estimated using engineering judgment.On the four-dimensional space, 1000 model instances are generated using random sampling.Uncertainty sources present in this case study are listed in Table 3.As the 2D model does not consider the corner effects, its predictions are underestimated when compared with 3D model predictions.Therefore, the uncertainty source related to the 2D model simplification has only positive magnitudes.Other sources are estimated based on literature [23,26,53] and quantification methods reported in [60].

Adaptation of the hierarchical algorithm for excavation
The uncertainty magnitude of 3D effects involved in an excavation case study is influenced by model parameter values [60].Therefore, falsification thresholds differ for each model instance.Consequently, the methodology requires the evaluation of the subset widths at sensor locations in the hierarchical algorithm.For each model instance, the combined uncertainty distribution at each measurement point location is computed.Then, the average of standard deviations mean,i of model instance combined uncertainty distributions at a measurement point i is calculated.The subset width at this measurement point is equal to 4 × mean,i .This value rep- resents the difference between upper and lower threshold bounds that are defined as two standard deviations from the mean value of predictions, similar to bridge case studies [40].Finally, the information entropy of a single measurement point is evaluated by Eq. ( 5).Results show that 30 measurements have been detected as outliers.These outliers fall into two categories.The first category involves measurements at the bottom of the excavation.In this case, the model predicts zero value inclination and has a distribution of combined uncertainties with a small standard deviation.A small deviation (in absolute values) of the measurements may result in the measurement being qualified as an outlier.This category includes most of the suspicious measurements, representing 20 out of the 28 measurements.The second category includes faulty measurement that occurs randomly.
Figure 8 presents the results of the outlier detection for two measurements.In the first case (Fig. 8a), the measurement falls within the distributions of the model predictions and combined uncertainties, while the second measurement value is clearly outside the expected range (Fig. 8b).Measurement 237 is thus taken to be a reliable measurement,  while Measurement 435 is an outlier and is removed from the measurement set.

Information entropy
The second step involves removing data points that are noninformative.This data cleansing is made by evaluating the information entropy of model instance predictions for each measurement point individually.
When the measurement point is associated with an information entropy value equal to zero (Eq.6), this means that the variability of model instance predictions is very low compared to the combined uncertainties for this data point, meaning that uncertainties have much larger magnitudes than prediction ranges.In such cases, no model instances will be falsified by this measurement, and this measurement is thus not providing any information.This is typically the case when the signal-to-noise ratio of a measurement point is low.Nonetheless, the physical justification of a noninformative measurement point may not be trivial.
Figures 8c, d present the distribution of information entropy values for the 498 measurement points after outlier detection (y-axis bounded to 50).The 320 measurement points have a zero entropy value, meaning that they can be removed from the measurement sets without affecting the information gain during data interpretation.The remaining 148 measurements have non-zero information entropy and are potentially providing useful information.These data points constitute the initial sets of potential measurements for the next step of the process.

Joint entropy assessment
In this stage, the hierarchical algorithm is used to evaluate the joint entropy values of sets of measurements.The aim is to find the smallest number of measurements that maximize the joint entropy (i.e., the expected information gain).As joint entropy values and selected measurement points are affected by the hypothesis of the number of measurements that will be needed (Eq.6), it is an iterative process.
Figure 9 presents the joint entropy values with respect to the number of measurements.Results of two hypotheses on the number of measurements are shown: 528 measurements (initial number of measurements), and 36 after the iterative process that is proposed in this paper.When the hypothesis involves a smaller number of sensors, the number of useful measurements and the maximum joint entropy values increase.As the Šidák correction increases threshold bound values with the number of measurements (Eq.5), smaller information gain is possible with more measurements for this hypothesis.

Optimal set of measurements
On the initial 528 and 148 measurements after data cleansing, only 36 measurements provide the most useful information (Table 4).These measurements constitute the optimal set of measurements for data interpretation.This corresponds to a reduction of 93% of the initial set and a reduction of 76% compared with the set after data cleansing.This selection must not compromise the information gain, and this will be assessed in the last step of the analysis.

Result validation
The last step involves validating the measurement point selection by comparing the true information of the selected set of measurements with the entire set after the outlier detection methodology.This step ensures that no information is lost during the data filtering.Field measurements from inclinometers (except I1 and I7) and all phases of excavation are shown in Fig. 10.I1 and I7 are not included due to the small deflection measured at all excavation phases.
Table 5 presents the comparison of the falsification performance and CMS of sets of measurements.Five sets are shown, and they include various numbers of measurements.The falsification performance significantly increases with the data selection methodology proposed in this paper.From the initial 20,000 model instances, only 55% are falsified when all measurements are considered in the falsification process.This falsification performance increases to 89% with the proposed methodology, showing a significant increase in information gain.As the Šidák correction also increases threshold bounds for falsification (Eq.3), adding more data in the measurement set decreases the identification performance as most data do not provide any information.
Compared to the previous methodology [49], a similar data set is selected (36 instead of 34 measurements).However, the information gain is improved as an additional 781 (26%) model instances have been falsified as not the same measurements are included in the set.This result is explained by the iterative process proposed in the methodology that  helps evaluate the information gain from each measurement point more accurately.A two-step process has been introduced in the methodology section.The first step involves the falsification performance of all non-informative data points.Following the results of the data cleansing section, this set includes 148 measurement points.When the falsification is performed with this set, all initial model instances remain in the candidate model set, confirming the usefulness of these measurements.The second step of the validation process involves comparing the falsification performance of the entire data set after the outlier detection with the set after the measurement point selection.In the present case study, the latter set enables a better falsification performance than the entire data set.Results of the measurement point selection are thus validated.Explanations for this better falsification performance are provided in the next section.

Discussion
For a constant probability of identification, the measurement point selection algorithm helps significantly reduce data sets required for structural identification without compromising the falsification performance.In the Bukit Panjang case study, such significant data selection has even improved the falsification performance.This result contradicts the paradigm in information theory "information never hurts" [64,65].
This contradiction can be partially explained by the data interpretation methodology used.EDMF [23], used in the study, involves the Šidák correction to maintain a given probability of identification in the data interpretation.By using this correction, threshold bounds increase when the number of measurements increases.This increase leads to a reduction in the information gain of measurements with an increase in the number of measurements.The contradiction with the information-never-hurts paradigm has also been observed for the residual minimization method [49].Additional work should be performed to evaluate the influence of the number of measurements on the information gain when traditionally implemented BMU [30,32,37] is used for data interpretation.
Figure 11 presents the influence of the number of measurements on the measurement point performance in terms of information gain for the full-scale case study.Results are shown for the expected information gain, calculated using information entropy (Eq.5) and observed information gain (Eq.8) using the best measurement of the case study.The falsification performance of the measurement point decreases with an increase in the number of measurements, showing that adding measurements to the data set reduces the individual information gained by each measurement point.The optimal set of measurements is thus a tradeoff between adding a new measurement with unique information and reducing the information gain of other measurements in the set.
The expected information gain (Fig. 11a), measured with information entropy, also decreases with increasing the number of measurements.As the information gain per measurement may vary with respect to the size of the data set, an iterative process to reach the optimal set of measurement points is necessary.This iterative process also leads to a new selection of measurement points as the unique information gain provided by each measurement is more accurately assessed.Even if the total number of measurement points is similar, the falsification performance of the measurement sets may be significantly improved by the iterative process, as shown in Table 5.
The following limitations of the work are recognized.In the case study, the combined uncertainties used in the hierarchical algorithm are calculated as the mean values of combined uncertainties associated with the 1000 model instances in order to reduce the computational time.This simplification may affect the evaluation of the expected measurement point performance, especially in terms of joint entropy evaluation.Also, this simplification may explain small discrepancies between the results of the hierarchical algorithm and the falsification performance of measurement points.Additionally, Monte Carlo simulations (using 1,000,000 samples) are used to generate threshold bounds which may slightly affect evaluations of expected and observed performances.
This study shows that algorithms that were initially built for measurement system design algorithms, such as the hierarchical algorithm, have the potential to be used for new applications such as measurement point selection.In the context of smart cities and the Internet of Things (IoT), design frameworks are increasingly becoming data-dependent.However, when inexpensive sensor devices are available, the amount of data increases exponentially.In these situations, large data sets are difficult to store, manage and analyze.This often results in weak data interpretation.A systematic and rational method for measurement point selection supports the extraction of valuable information among data sets and leads to the potential for better and more timely decision-making than that occurs in current practice.

Conclusions
In this paper, a methodology is proposed to select an optimal set of measurements in large data sets.A full-scale case study has been used to illustrate the proposed methodology.This case study involves updating the numerical model of the excavation of the Bukit Panjang Hawker Centre in Singapore through a model falsification approach.Results have shown that refined data sets provide more information than full data sets for the same probability of identification.Specific conclusions are as follows: • The methodology for measurement point selection supports engineers in the task of selecting informative measurement points to improve model-updating performance.
• Reducing the number of measurements using the proposed methodology helps avoid over-instrumentation without reducing the quality of the identification.
• By comparing identification performance using all measurement points and the subsets, methodologies for measurement point selection are effectively validated using field measurements.

Fig. 2
Fig. 2 Generation of model instance predictions

Fig. 3
Fig. 3 Data cleansing step.a Outlier detection methodology; b data filtering using information entropy

Fig. 5
Fig. 5 Model updating using error domain model falsification.It is common for model predictions to cluster either at the upper or lower threshold due to the bias in the uncertainties

123 4 . 4
Data filtering4.4.1 Outlier detectionThis step involves data cleansing to remove doubtful measurements.The outlier detection methodology removes data from the measurement sets if the measurements fall outside the 99% interval of the distribution of the combined uncertainties sum with the distribution of model instance predictions.

Fig. 8
Fig. 8 Results of the data cleansing.a Outlier detection: Measurement 237; b Outlier detection: Measurement 435.c Distribution of information entropy values of measurement points; d Zoom on the distribution of the information entropy values of measurement points (y-axis bounded to 50)

Fig. 9
Fig.9 Joint entropy as a function of number of measurement points for a constant probability of identification.Both data sets show a diminishing increase in joint entropy to nearly zero after 15 measurements.However, the highest entropy is achieved with the smaller data set

Fig. 10
Fig. 10 Deflection measurements by the inclinometers at the four phases of the excavation

Fig. 11
Fig. 11 Information gain of the best measurement point as a function of the number of measurements.a Expected information gain (information entropy); b observed information gain (falsification performance)

Table 4
Size of measurement sets with respect to the step of the methodology