Accuracy Assessment and Comparative Analysis of IDW, Spline and Kriging in Spatial Interpolation of Landform (Topography): An Experimental Study

It is practically impossible and unnecessary to obtain spatial-temporal information of any given continuous phenomenon at every point within a given geographic area. The most practical approach has always been to obtain information about the phenomenon as in many sample points as possible within the given geographic area and estimate the values of the unobserved points from the values of the observed points through spatial interpolation. However, it is important that users understand that different interpolation methods have their strength and weaknesses on different datasets. It is not correct to generalize that a given interpolation method (e.g. Kriging, Inverse Distance Weighting (IDW), Spline etc.) does better than the other without taking into cognizance, the type and nature of the dataset and phenomenon involved. In this paper, we theoretically, mathematically and experimentally evaluate the performance of Kriging, IDW and Spline interpolation methods respectively in estimating unobserved elevation values and modeling landform. This paper undertakes a comparative analysis based on the prediction mean error, prediction root mean square error and cross validation outputs of these interpolation methods. Experimental results for each of the method on both biased and normalized data show that Spline provided a better and more accurate interpolation within the sample space than the IDW and Kriging methods. The choice of an interpolation method should be phenomenon and data set structure dependent.

Study of Quintic Spline Interpolation and Generated Velocity Profile for High Speed Machining

Modern high speed machining (HSM) machine tools often operates at high speed and high feedrate with high ac- celerations,in order to deliver the rapid feed motion.This paper presents an interpolation algorithm to generate continuous quintic spline toolpaths,with a constant travel increment at each step,while the smoother accelerations and jerks of two-order curve are obtained.Then an approach for reducing the feedrate fluctuation in high speed spline interpolation is presented.The presented ap- proach has been validated to quickly,reliably and effective with the simulation.

The Cubic Trigonometric Automatic Interpolation Spline

A class of cubic trigonometric interpolation spline curves with two parameters is presented in this paper. The spline curves can automatically interpolate the given data points and become C^2 interpolation curves without solving equations system even if the interpolation conditions are fixed. Moreover, shape of the interpolation spline curves can be globally adjusted by the two parameters. By selecting proper values of the two parameters,the optimal interpolation spline curves can be obtained.

Smooth cutting pattern generation technique for membrane structures using geodesic line on subplane and spline interpolation

Practical techniques for smooth geodesic patterning of membrane structures were investigated.For the geodesic search,adjustment of the subplane of the extracted elements series was proposed,and various spline approximation methods were used to flatten the strip for the generation of a smooth pattern.This search approach is very simple,and the geodesic line could be easily attained by the proposed method without the need for a difficult computation method.Smooth cutting patterning can also be generated by spline approximation without the noise in discrete nodal information.Additionally,the geodesic cutting pattern saved about 21%of the required area for the catenary model due to the reduction of the curvature of the planar pattern seam line.

Performance of global look-up table strategy in digital image correlation with cubic B-spline interpolation and bicubic interpolation

Global look-up table strategy proposed recently has been proven to be an efficient method to accelerate the interpolation, which is the most time-consuming part in the iterative sub-pixel digital image correlation （DIC） algorithms. In this paper, a global look-up table strategy with cubic B-spline interpolation is developed for the DIC method based on the inverse compositional Gauss-Newton （IC-GN） algorithm. The performance of this strategy, including accuracy, precision, and computation efficiency, is evaluated through a theoretical and experimental study, using the one with widely employed bicubic interpolation as a benchmark. The global look-up table strategy with cubic B-spline interpolation improves significantly the accuracy of the IC-GN algorithm-based DIC method compared with the one using the bicubic interpolation, at a trivial price of computation efficiency.

A new method to build boundary conditions for nonuniform B-splines interpolation

Discusses a new method to build boundary conditions for nonuniform B splines interpolation based on the curvature parameters with two advantages: no derivative of curve end is required and zero curvature at curve end is avoided, so that the shapes of the two end segments of curve can be controlled easily by adjusting the curvature parameters to meet the designer’s requirements.

Non-overshooting and Non-undershooting Cubic Spline Interpolation for Empirical Mode Decomposition

To suppress the overshoots and undershoots in the envelope fitting for empirical mode decomposition （EMD）, an alternative cubic spline interpolation method without overshooting and undershooting is proposed. On the basis of the derived slope constraints of knots of a non-overshooting and non-undershooting cubic interpolant, together with ＂not-a-knot＂ conditions the cubic spline interpolants are constructed by replacing the requirement for equal second order derivatives at every knot with Brodlie＇ s derivative formula. Analysis and simulation experiments show that this approach can effectively avoid generating new extrema, shifting or exaggerating the existing ones in a signal, and thus significantly improve the decomposition performance of EMD.

Trivariate Polynomial Natural Spline for 3D Scattered Data Hermit Interpolation

Consider a kind of Hermit interpolation for scattered data of 3D by trivariate polynomial natural spline, such that the objective energy functional （with natural boundary conditions） is minimal. By the spline function methods in Hilbert space and variational theory of splines, the characters of the interpolation solution and how to construct it are studied. One can easily find that the interpolation solution is a trivariate polynomial natural spline. Its expression is simple and the coefficients can be decided by a linear system. Some numerical examples are presented to demonstrate our methods.

Cubic Spline Interpolation on a Class of Triangulations

In this paper, we consider spaces of cubic C^1-spline on a class of triangulations. By using the inductive algorithm, the posed Lagrange interpolation sets are constructed for cubic spline space. It is shown that the class of triangulations considered in this paper are nonsingular for S1/3 spaces. Moreover, the dimensions of those spaces exactly equal to L. L. Schuraaker＇s low bounds of the dimensions. At the end of this paper, we present an approach to construct triangulations from any scattered planar points, which ensures that the obtained triangulations for S1/3 space are nonsingular.

Improvement of Orbit Prediction Algorithm for Spacecraft Through Simplified Precession-Nutation Model Using Cubic Spline Interpolation Method

For the on-orbit flight missions,the model of orbit prediction is critical for the tasks with high accuracy requirement and limited computing resources of spacecraft.The precession-nutation model,as the main part of extended orbit prediction,affects the efficiency and accuracy of on-board operation.In this paper,the previous research about the conversion between the Geocentric Celestial Reference System and International Terrestrial Reference System is briefly summarized,and a practical concise precession-nutation model is proposed for coordinate transformation computation based on Celestial Intermediate Pole(CIP).The idea that simplifying the CIP-based model with interpolation method is driven by characteristics of precession-nutation parameters changing with time.A cubic spline interpolation algorithm is applied to obtain the required CIP coordinates and Celestial Intermediate Origin locator.The complete precession nutation model containing more than 4000 parameters is simplified to the calculation of a cubic polynomial,which greatly reduces the computational load.In addition,for evaluating the actual performance,an orbit propagator is built with the proposed simplified precession-nutationmodel.Compared with the orbit prediction results obtained by the truncated series of IAU2000/2006 precession-nutation model,the simplified precession-nutation model with cubic spline interpolation can significantly improve the accuracy of orbit prediction,which implicates great practical application value in further on-orbit missions of spacecraft.

ERROR BOUNDS IN PERIODIC QUARTIC SPLINE INTERPOLATION

In this paper we develop periodic quartic spline interpolation theory which,in general,gives better fus to continuous functions than does the existing quintic spline interpolation theory.The main theorem of the paper is to establish that r=0,1,2,3.Also,the nanperiodic cases cannot be constructed empoly-ing the methodology of this paper because that will involve several other end conditions entirely different than(1,10).

A polynomial smooth epsilon-support vector regression based on cubic spline interpolation

Regression analysis is often formulated as an optimization problem with squared loss functions. Facing the challenge of the selection of the proper function class with polynomial smooth techniques applied to support vector regression models, this study takes cubic spline interpolation to generate a new polynomial smooth function ｜×｜ε^ 2, in g-insensitive support vector regression. Theoretical analysis shows that Sε^2 -function is better than pε^2 -function in properties, and the approximation accuracy of the proposed smoothing function is two order higher than that of classical pε^2 -function. The experimental data shows the efficiency of the new approach.

Combining Cubic Spline Interpolation and Fast Fourier Transform to Extend Measuring Range of Reflectometry

The reflectometry is a common method used to measure the thickness of thin films. Using a conventional method,its measurable range is limited due to the low resolution of the current spectrometer embedded in the reflectometer.We present a simple method, using cubic spline interpolation to resample the spectrum with a high resolution,to extend the measurable transparent film thickness. A large measuring range up to 385 m in optical thickness is achieved with the commonly used system. The numerical calculation and experimental results demonstrate that using the FFT method combined with cubic spline interpolation resampling in reflectrometry, a simple,easy-to-operate, economic measuring system can be achieved with high measuring accuracy and replicability.

OPTIMAL ERROR BOUNDS FOR THE CUBIC SPLINE INTERPOLATION OF LOWER SMOOTH FUNCTIONS(1)

In this paper,the kernel of the cubic spline interpolation is given.An optimal error bound for the cu- bic spline interpolation of lower smooth functions is obtained.

Inverse Spline Interpolation for All-time Resistivity of Central-Loop TEM

A convenient numerical calculation method (inverse spline interpolation) for all-time apparent resistivity intransient electromagnetic method (TEM) is proposed in this paper. Characteristic of early and late normalized inductiveelectromotive force was investigated. According to the turning point, the transient process is divided into the earlyphase, the turning point, and the late phase. Afterwards, apparent resistivity is obtained through inverse spline interpo-lation in the early and the late phases, respectively. Finally, the resistivities of the early-time and the late-time wereconnected together by the turning point. The result shows that the inverse spline method is feasible and the method alsolays a foundation for initial model construction in the TEM automatic inversion.

Downhole microseismic data reconstruction and imaging based on combination of spline interpolation and curveletsparse constrained interpolation

When cause of the aliasing lack probl using borehole sensors and microseimic events to image, spatial aliasing often occurred be- of sensors underground and the distance between the sensors which were too large. To solve em, data reconstruction is often needed. Curvelet transform sparsity constrained inversion was widely used in the seismic data reconstruction field for its anisotropic, muhiscale and local basis. However, for the downhole ease, because the number of sampling point is mueh larger than the number of the sensors, the advantage of the cnrvelet basis can＇t perform very well. To mitigate the problem, the method that joints spline and curvlet-based compressive sensing was proposed. First, we applied the spline interpolation to the first arri- vals that to be interpolated. And the events are moved to a certain direction, such as horizontal, which can be represented by the curvelet basis sparsely. Under the spasity condition, curvelet-based compressive sensing was applied for the data, and directional filter was also used to mute the near vertical noises. After that, the events were shifted to the spline line to finish the interpolation workflow. The method was applied to a synthetic mod- el, and better result was presented than using curvelet transform interpolation directly. We applied the method to a real dataset, a mieroseismic downhole observation field data in Nanyang, using Kirchhoff migration method to image the microseimic event. Compared with the origin data, artifacts were suppressed on a certain degree.

C^n-CONTINUOUS B-TYPE SPLINE CURVES WITH ITSLOCALIZATION INTERPOLATION AND APPROXIMATION

This paper presents a class of Cn- continuous B- type spline curves with some paramet- ric factors.The length of their local support is equal to4.Taking the different values of the parametric factors,the curves can become free- type curves or interpolate a set of given points even mix the both cases.When the parametric factors satisfy the certain conditions,the degrees of the curves can be decreased as low as possible.Besides,when all the parametric factors tend to zero,the curves globally approximate to the control polygon.

Modeling the Dynamic Gravity Variations of Northeastern Margin of Qinghai-Xizang (Tibet) Plateau by Using Bicubic Spline Interpolation Function

In this paper, the spatial-temporal gravity variation patterns of the northeastern margin of Qinghal-Xizang （Tibet） Plateau in 1992 - 2001 are modeled using bicubic spline interpolation functions and the relations of gravity change with seismicity and tectonic movement are discussed preliminarily. The results show as follows： ① Regional gravitational field changes regularly and the gravity abnormity zone or gravity concentration zone appears in the earthquake preparation process; ②In the significant time period, the gravity variation shows different features in the northwest, southeast and northeast parts of the surveyed region respectively, with Lanzhou as its boundary;③The gravity variation distribution is basically identical to the strike of tectonic fault zone of the region, and the contour of gravity variation is closely related to the fault distribution.

Thin-Plate Spline Analysis of Diversity in Field Populations of Codling Moth, Cydia pomonella L. （Lep., Tortricidae） in North West of Iran

In order to study morphological diversity of codling moth, Cydia pomonella （L.） using thin-plate spline analysis, nine geographical populations from four north western provinces of Iran namely East Azarbayjan, West Azarbayjan, Ardebil and Zandjan were collected during 2003 and 2004. 575 and 564 images were prepared from fore and hind wings, respectively. Then 15 and 11 landmarks were determined from fore and hind wings, respectively. With transforming of landmark＇s two dimensional coordinate data into partial warp scores, 26 and 18 scores were generated for fore and hind wings, respectively. Cluster analysis based on wing shape variables using Ward＇s algorithm assigned nine geographical populations into two groups. The pattern of grouping based on fore and hind wings was different in both sexes. Principal component analysis revealed discrimination between geographic populations and confirmed the result of cluster analysis. Among environmental parameters, wind speed showed the highest correlation with wing shape variables. Non significant correlation was observed between geographic and morphological distance matrices as revealed by Mantel test.