{"title":"On the BIC for determining the number of control points in B-spline surface approximation in case of correlated observations","authors":"G. Kermarrec, H. Alkhatib","doi":"10.1515/jogs-2020-0110","DOIUrl":null,"url":null,"abstract":"Abstract B-spline curves are a linear combination of control points (CP) and B-spline basis functions. They satisfy the strong convex hull property and have a fine and local shape control as changing one CP affects the curve locally, whereas the total number of CP has a more general effect on the control polygon of the spline. Information criteria (IC), such as Akaike IC (AIC) and Bayesian IC (BIC), provide a way to determine an optimal number of CP so that the B-spline approximation fits optimally in a least-squares (LS) sense with scattered and noisy observations. These criteria are based on the log-likelihood of the models and assume often that the error term is independent and identically distributed. This assumption is strong and accounts neither for heteroscedasticity nor for correlations. Thus, such effects have to be considered to avoid under-or overfitting of the observations in the LS adjustment, i.e. bad approximation or noise approximation, respectively. In this contribution, we introduce generalized versions of the BIC derived using the concept of quasi- likelihood estimator (QLE). Our own extensions of the generalized BIC criteria account (i) explicitly for model misspecifications and complexity (ii) and additionally for the correlations of the residuals. To that aim, the correlation model of the residuals is assumed to correspond to a first order autoregressive process AR(1). We apply our general derivations to the specific case of B-spline approximations of curves and surfaces, and couple the information given by the different IC together. Consecutively, a didactical yet simple procedure to interpret the results given by the IC is provided in order to identify an optimal number of parameters to estimate in case of correlated observations. A concrete case study using observations from a bridge scanned with a Terrestrial Laser Scanner (TLS) highlights the proposed procedure.","PeriodicalId":44569,"journal":{"name":"Journal of Geodetic Science","volume":null,"pages":null},"PeriodicalIF":0.9000,"publicationDate":"2020-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"2","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Journal of Geodetic Science","FirstCategoryId":"1085","ListUrlMain":"https://doi.org/10.1515/jogs-2020-0110","RegionNum":0,"RegionCategory":null,"ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q4","JCRName":"REMOTE SENSING","Score":null,"Total":0}
引用次数: 2
Abstract
Abstract B-spline curves are a linear combination of control points (CP) and B-spline basis functions. They satisfy the strong convex hull property and have a fine and local shape control as changing one CP affects the curve locally, whereas the total number of CP has a more general effect on the control polygon of the spline. Information criteria (IC), such as Akaike IC (AIC) and Bayesian IC (BIC), provide a way to determine an optimal number of CP so that the B-spline approximation fits optimally in a least-squares (LS) sense with scattered and noisy observations. These criteria are based on the log-likelihood of the models and assume often that the error term is independent and identically distributed. This assumption is strong and accounts neither for heteroscedasticity nor for correlations. Thus, such effects have to be considered to avoid under-or overfitting of the observations in the LS adjustment, i.e. bad approximation or noise approximation, respectively. In this contribution, we introduce generalized versions of the BIC derived using the concept of quasi- likelihood estimator (QLE). Our own extensions of the generalized BIC criteria account (i) explicitly for model misspecifications and complexity (ii) and additionally for the correlations of the residuals. To that aim, the correlation model of the residuals is assumed to correspond to a first order autoregressive process AR(1). We apply our general derivations to the specific case of B-spline approximations of curves and surfaces, and couple the information given by the different IC together. Consecutively, a didactical yet simple procedure to interpret the results given by the IC is provided in order to identify an optimal number of parameters to estimate in case of correlated observations. A concrete case study using observations from a bridge scanned with a Terrestrial Laser Scanner (TLS) highlights the proposed procedure.
b样条曲线是控制点(CP)与b样条基函数的线性组合。它们满足强凸包特性,具有精细的局部形状控制,改变一个凸包只会局部影响曲线,而凸包总数对样条控制多边形的影响更为普遍。信息准则(IC),如Akaike IC (AIC)和Bayesian IC (BIC),提供了一种确定最佳CP数的方法,使b样条近似在最小二乘(LS)意义上与分散和噪声观测值最优拟合。这些标准是基于模型的对数似然,并且通常假设误差项是独立的和同分布的。这种假设是强有力的,既没有考虑到异方差,也没有考虑到相关性。因此,必须考虑这些影响,以避免LS平差中观测值的欠拟合或过拟合,即不良近似或噪声近似。在这篇文章中,我们介绍了利用拟似然估计量(QLE)的概念推导出的BIC的广义版本。我们自己对广义BIC标准的扩展(i)明确地考虑了模型的错误规范和复杂性(ii),另外还考虑了残差的相关性。为此,假设残差的相关模型对应于一阶自回归过程AR(1)。我们将我们的一般推导应用于曲线和曲面的b样条近似的具体情况,并将不同IC给出的信息耦合在一起。随后,为了确定在相关观测的情况下要估计的参数的最优数量,提供了一个说明性但简单的程序来解释由集成电路给出的结果。使用地面激光扫描仪(TLS)扫描的桥梁观测结果的具体案例研究突出了所提出的程序。