The dynamic optimal interpolation(DOI)method is a technique based on quasi-geostrophic dynamics for merging multi-satellite altimeter along-track observations to generate gridded absolute dynamic topography(ADT).Compa...The dynamic optimal interpolation(DOI)method is a technique based on quasi-geostrophic dynamics for merging multi-satellite altimeter along-track observations to generate gridded absolute dynamic topography(ADT).Compared with the linear optimal interpolation(LOI)method,the DOI method can improve the accuracy of gridded ADT locally but with low computational efficiency.Consequently,considering both computational efficiency and accuracy,the DOI method is more suitable to be used only for regional applications.In this study,we propose to evaluate the suitable region for applying the DOI method based on the correlation between the absolute value of the Jacobian operator of the geostrophic stream function and the improvement achieved by the DOI method.After verifying the LOI and DOI methods,the suitable region was investigated in three typical areas:the Gulf Stream(25°N-50°N,55°W-80°W),the Japanese Kuroshio(25°N-45°N,135°E-155°E),and the South China Sea(5°N-25°N,100°E-125°E).We propose to use the DOI method only in regions outside the equatorial region and where the absolute value of the Jacobian operator of the geostrophic stream function is higher than1×10^(-11).展开更多
This paper introduces the definition of the Orthogonal Type Node Configuration and discusses the corresponding multivariate Lagrange, Hermite and Birkhoff interpolation problems in high dimensional space R s(s>2). ...This paper introduces the definition of the Orthogonal Type Node Configuration and discusses the corresponding multivariate Lagrange, Hermite and Birkhoff interpolation problems in high dimensional space R s(s>2). This node configuration can be considered to be a kind of extension of the Cross Type Node Configuration , in R 2 to high dimensional spaces. And the Mixed Type Node Configuration in R s(s>2) is also discussed in this paper in an example.展开更多
We construct general structures of one and two variable interpolation function, without depending on the existence of divided difference or inverse differences, and we also discuss the block based osculatory interpola...We construct general structures of one and two variable interpolation function, without depending on the existence of divided difference or inverse differences, and we also discuss the block based osculatory interpolation in one variable case. Clearly, our method offers marly flexible interpolation schemes for choices. Error terms for the interpolation are determined and numerical examples are given to show the effectlveness of the results.展开更多
As a branch of quantum image processing,quantum image scaling has been widely studied.However,most of the existing quantum image scaling algorithms are based on nearest-neighbor interpolation and bilinear interpolatio...As a branch of quantum image processing,quantum image scaling has been widely studied.However,most of the existing quantum image scaling algorithms are based on nearest-neighbor interpolation and bilinear interpolation,the quantum version of bicubic interpolation has not yet been studied.In this work,we present the first quantum image scaling scheme for bicubic interpolation based on the novel enhanced quantum representation(NEQR).Our scheme can realize synchronous enlargement and reduction of the image with the size of 2^(n)×2^(n) by integral multiple.Firstly,the image is represented by NEQR and the original image coordinates are obtained through multiple CNOT modules.Then,16 neighborhood pixels are obtained by quantum operation circuits,and the corresponding weights of these pixels are calculated by quantum arithmetic modules.Finally,a quantum matrix operation,instead of a classical convolution operation,is used to realize the sum of convolution of these pixels.Through simulation experiments and complexity analysis,we demonstrate that our scheme achieves exponential speedup over the classical bicubic interpolation algorithm,and has better effect than the quantum version of bilinear interpolation.展开更多
In this study,a wavelet multi-resolution interpolation Galerkin method(WMIGM)is proposed to solve linear singularly perturbed boundary value problems.Unlike conventional wavelet schemes,the proposed algorithm can be r...In this study,a wavelet multi-resolution interpolation Galerkin method(WMIGM)is proposed to solve linear singularly perturbed boundary value problems.Unlike conventional wavelet schemes,the proposed algorithm can be readily extended to special node generation techniques,such as the Shishkin node.Such a wavelet method allows a high degree of local refinement of the nodal distribution to efficiently capture localized steep gradients.All the shape functions possess the Kronecker delta property,making the imposition of boundary conditions as easy as that in the finite element method.Four numerical examples are studied to demonstrate the validity and accuracy of the proposedwavelet method.The results showthat the use ofmodified Shishkin nodes can significantly reduce numerical oscillation near the boundary layer.Compared with many other methods,the proposed method possesses satisfactory accuracy and efficiency.The theoretical and numerical results demonstrate that the order of theε-uniform convergence of this wavelet method can reach 5.展开更多
Missing value is one of the main factors that cause dirty data.Without high-quality data,there will be no reliable analysis results and precise decision-making.Therefore,the data warehouse needs to integrate high-qual...Missing value is one of the main factors that cause dirty data.Without high-quality data,there will be no reliable analysis results and precise decision-making.Therefore,the data warehouse needs to integrate high-quality data consistently.In the power system,the electricity consumption data of some large users cannot be normally collected resulting in missing data,which affects the calculation of power supply and eventually leads to a large error in the daily power line loss rate.For the problem of missing electricity consumption data,this study proposes a group method of data handling(GMDH)based data interpolation method in distribution power networks and applies it in the analysis of actually collected electricity data.First,the dependent and independent variables are defined from the original data,and the upper and lower limits of missing values are determined according to prior knowledge or existing data information.All missing data are randomly interpolated within the upper and lower limits.Then,the GMDH network is established to obtain the optimal complexity model,which is used to predict the missing data to replace the last imputed electricity consumption data.At last,this process is implemented iteratively until the missing values do not change.Under a relatively small noise level(α=0.25),the proposed approach achieves a maximum error of no more than 0.605%.Experimental findings demonstrate the efficacy and feasibility of the proposed approach,which realizes the transformation from incomplete data to complete data.Also,this proposed data interpolation approach provides a strong basis for the electricity theft diagnosis and metering fault analysis of electricity enterprises.展开更多
The refinement of the 1cm geoid holds significant importance in engineering applications,where the accuracy of the geoid is frequently constrained by its resolution.However,there has been limited exploration into the ...The refinement of the 1cm geoid holds significant importance in engineering applications,where the accuracy of the geoid is frequently constrained by its resolution.However,there has been limited exploration into the specific relationship between geoid resolution and accuracy.This article aims to address this gap by thoroughly examining said relationship.This study employs the CapRCR modification to calculate the geoid of Gongzhuling City.The findings indicate that the accuracy can be enhanced by 2%to 9%through encryption of the geoid resolution from 2.5′×2.5′to 1.5′×1.5′.Furthermore,this improvement can be augmented by 15%to 21%through encryption of the gravity anomaly to the same resolution.The accuracy of the geoid exhibits a linear relationship with the resolution of gravity anomalies.The theoretical accuracy of the geoid,excluding integration errors,is determined to be 1.21 cm.In engineering applications,opting for a resolution of 30''×30''can result in the geoid containing only an integration error of 2 mm.However,to attain an accuracy of 1 cm,it becomes imperative to further mitigate data errors.Consequently,the final refined quasi-geoid accuracy is established at 0.56 mm.展开更多
The bivariate interpolation in two dimensional space R2 is more complicated than that in one dimensional space R, because there is no Haar space of continuous functions in R2. Therefore, the bivariate interpolation ha...The bivariate interpolation in two dimensional space R2 is more complicated than that in one dimensional space R, because there is no Haar space of continuous functions in R2. Therefore, the bivariate interpolation has not a unique solution for a set of arbitrary distinct pairwise points. In this work, we suggest a type of basis which depends on the points such that the bivariate interpolation has the unique solution for any set of distinct pairwise points. In this case, the matrix of bivariate interpolation has the semi inherited factorization.展开更多
The accuracy of interpolation models applied to groundwater depends, among other factors, on the interpolation method chosen. Therefore, it is necessary to compare different approaches. For this, different methods of ...The accuracy of interpolation models applied to groundwater depends, among other factors, on the interpolation method chosen. Therefore, it is necessary to compare different approaches. For this, different methods of interpolation of nitrate concentrations were contrasted in sixty-seven wells in an aquifer in Aguascalientes, Mexico. Four general interpolation methods were used in ArcGIS 10.5 to make the maps: IDW, Kriging, Natural Neighbor and Spline. In the modeling, only method type was varied. The input parameters (location, temporality, and nitrate concentration) were the same in the four interpolations;despite this, different maximum and minimum values were obtained for each interpolation method: for IDW, 0.2 to 22.0 mg/l, for Kriging, 3.5 to 16.5 mg/l, for Natural Neighbor, 0.3 to 21.7 mg/l and for Spline −30.8 to 37.2 mg/l. Finally, an assessment of the maps obtained was conducted by comparing them with the Official Mexican Standard (OMS), where 24 of the 67 wells were found outside the 10 mg/l that the OMS establishes as maximum permissible limit for human consumption. Taking as a starting point the measured values of nitrates (0.25 to 22.12 mg/l), as well as the spatial distribution of the interpolated values, it was determined that the Krigging method best fitted the data measured in the wells within the studied aquifer.展开更多
The purpose of this paper is to investigate the spatial interpolation of rainfall variability with deterministic and geostatic inspections in the Prefecture of Kilkis (Greece). The precipitation data where recorded fr...The purpose of this paper is to investigate the spatial interpolation of rainfall variability with deterministic and geostatic inspections in the Prefecture of Kilkis (Greece). The precipitation data where recorded from 12 meteorological stations in the Prefecture of Kilkis for 36 hydrological years (1973-2008). The cumulative monthly values of rainfall were studied on an annual and seasonal basis as well as during the arid-dry season. In the deterministic tests, the I.D.W. and R.B.F. checks were inspected, while in the geostatic tests, Ordinary Kriging and Universal Kriging respectively. The selection of the optimum method was made based on the least Root Mean Square Error (R.M.S.E.), as well as on the Mean Error (M.E.), as assessed by the cross validation analysis. The geostatical Kriging also considered the impact of isotropy and anisotropy across all time periods of data collection. Moreover, for Universal Kriging, the study explored spherical, exponential and Gaussian models in various combinations. Geostatistical techniques consistently demonstrated greater reliability than deterministic techniques across all time periods of data collection. Specifically, during the annual period, anisotropy was the prevailing characteristic in geostatistical techniques. Moreover, the results for the irrigation and seasonal periods were generally comparable, with few exceptions where isotropic methods yielded lower (R.M.S.E.) in some seasonal observations.展开更多
General interpolation formulae for barycentric interpolation and barycen- tric rational Hermite interpolation are established by introducing multiple parameters, which include many kinds of barycentric interpolation a...General interpolation formulae for barycentric interpolation and barycen- tric rational Hermite interpolation are established by introducing multiple parameters, which include many kinds of barycentric interpolation and barycentric rational Her- mite interpolation. We discussed the interpolation theorem, dual interpolation and special cases. Numerical example is given to show the effectiveness of the method.展开更多
In this paper, we obtain a properly posed set of nodes for interpolation on a sphere. Moreover it is applied to construct properly posed set of nodes for Lagrange interpolation on the trivariate polynomial space of to...In this paper, we obtain a properly posed set of nodes for interpolation on a sphere. Moreover it is applied to construct properly posed set of nodes for Lagrange interpolation on the trivariate polynomial space of total degree n.展开更多
The radial basis function (RBF) interpolation approach proposed by Freedman is used to solve inverse problems encountered in well-logging and other petrophysical issues. The approach is to predict petrophysical prop...The radial basis function (RBF) interpolation approach proposed by Freedman is used to solve inverse problems encountered in well-logging and other petrophysical issues. The approach is to predict petrophysical properties in the laboratory on the basis of physical rock datasets, which include the formation factor, viscosity, permeability, and molecular composition. However, this approach does not consider the effect of spatial distribution of the calibration data on the interpolation result. This study proposes a new RBF interpolation approach based on the Freedman's RBF interpolation approach, by which the unit basis functions are uniformly populated in the space domain. The inverse results of the two approaches are comparatively analyzed by using our datasets. We determine that although the interpolation effects of the two approaches are equivalent, the new approach is more flexible and beneficial for reducing the number of basis functions when the database is large, resulting in simplification of the interpolation function expression. However, the predicted results of the central data are not sufficiently satisfied when the data clusters are far apart.展开更多
In order to get the spatial grid data of monthly precipitation and monthly average temperature of Sanjiangyuan area, the Co-Kriging (COK) and thin plate smoothing splines(TPS) interpolation methods were applied by usi...In order to get the spatial grid data of monthly precipitation and monthly average temperature of Sanjiangyuan area, the Co-Kriging (COK) and thin plate smoothing splines(TPS) interpolation methods were applied by using the climate data during 1971-2000 of 58 meteorological stations around Qinghai Province and the 3 arc-second digital elevation model (DEM) data. The performance was evaluated by the smallest statistical errors by general cross validation (GCV). Root-mean-squared predicted errors (RMSE) and mean absolute errors (MAE) were used to compare the performance of the two methods. The results showed that: 1) After combing covariates into the models, both methods performed better; 2) The performance of TPS was significantly better than COK: for monthly average temperature, the RMSE derived from TPS was 69.48% higher than COK, as MAE increased by 70.56%. And for monthly precipitation, the RMSE derived from TPS was 28.07% higher than COK, as MAE increased by 29.06%.展开更多
To alleviate the distortion of XRII X-ray image intensifier images in the C-arm CT computer tomography imaging system an algorithm based on the Delaunay triangulation interpolation is proposed.First the causes of the ...To alleviate the distortion of XRII X-ray image intensifier images in the C-arm CT computer tomography imaging system an algorithm based on the Delaunay triangulation interpolation is proposed.First the causes of the phenomenon the classical correction algorithms and the Delaunay triangulation interpolation are analyzed.Then the algorithm procedure is explained using flow charts and illustrations. Finally experiments are described to demonstrate its effectiveness and feasibility. Experimental results demonstrate that the Delaunay triangulation interpolation can have the following effects.In the case of the same center the root mean square distances RMSD and standard deviation STD between the corrected image with Delaunay triangulation interpolation and the ideal image are 5.760 4 ×10 -14 and 5.354 2 ×10 -14 respectively.They increase to 1.790 3 2.388 8 2.338 8 and 1.262 0 1.268 1 1.202 6 after applying the quartic polynomial model L1 and model L2 to the distorted images respectively.The RMSDs and STDs between the corrected image with the Delaunay triangulation interpolation and the ideal image are 2.489 × 10 -13 and 2.449 8 ×10 -13 when their centers do not coincide. When the quartic polynomial model L1 and model L2 are applied to the distorted images they are 1.770 3 2.388 8 2.338 8 and 1.269 9 1.268 1 1.202 6 respectively.展开更多
Aim To study wind tunnel test data interpolation methods for flight vehicle with aerodynamic axial asymmetry. Methods For different body aerodynamic roll angles, proper wind tunnel test schemes were selected and ...Aim To study wind tunnel test data interpolation methods for flight vehicle with aerodynamic axial asymmetry. Methods For different body aerodynamic roll angles, proper wind tunnel test schemes were selected and trigonometric series were used for aerodynamic interpolation. Results and Conclusion A simple and effective scheme for wind tunnel test and an accurate aerodynamic interpolation method are developed with satisfactory results.展开更多
A novel spatial interpolation method based on integrated radial basis function artificial neural networks (IRBFANNs) is proposed to provide accurate and stable predictions of heavy metals concentrations in soil at u...A novel spatial interpolation method based on integrated radial basis function artificial neural networks (IRBFANNs) is proposed to provide accurate and stable predictions of heavy metals concentrations in soil at un- sampled sites in a mountain region. The IRBFANNs hybridize the advantages of the artificial neural networks and the neural networks integration approach. Three experimental projects under different sampling densities are carried out to study the performance of the proposed IRBFANNs-based interpolation method. This novel method is compared with six peer spatial interpolation methods based on the root mean square error and visual evaluation of the distribution maps of Mn elements. The experimental results show that the proposed method performs better in accuracy and stability. Moreover, the proposed method can provide more details in the spatial distribution maps than the compared interpolation methods in the cases of sparse sampling density.展开更多
In this research, we present a seismic trace interpolation method which uses seismic data with surface-related multiples. It is different from conventional seismic data interpolation using information transformation o...In this research, we present a seismic trace interpolation method which uses seismic data with surface-related multiples. It is different from conventional seismic data interpolation using information transformation or extrapolation of adjacent channels for reconstruction of missing seismic data. In this method there are two steps, first, we construct pseudo-primaries by cross-correlation of surface multiple data to extract the missing near- offset information in multiples, which are not displayed in the acquired seismic record. Second, we correct the pseudo-primaries by applying a Least-squares Matching Filter (LMF) and RMS amplitude correction method in time and space sliding windows. Then the corrected pseudo-primaries can be used to fill the data gaps. The method is easy to implement, without the need to separate multiples and primaries. It extracts the seismic information contained by multiples for filling missing traces. The method is suitable for seismic data with surfacerelated multiples.展开更多
A direct interpolation algorithm for spatial curves on a surface is proposed for pen-cutting of sculptured surfaces. The algorithm can carry out the direct interpolation for projective curves lying on the sculptured s...A direct interpolation algorithm for spatial curves on a surface is proposed for pen-cutting of sculptured surfaces. The algorithm can carry out the direct interpolation for projective curves lying on the sculptured surface. It is based on the geometric and kinetic relationships between drive curves and cutter-contact (C-C) curves. It evaluates the parameter of drive curves corresponding to interpolation points by the Taylor formula. Then it gains coordinates of interpolation points indirectly by inverse calculation. Finally, it generates the motion commands for machine tool controller. This method extends the locus-controlled function in the computer numerical control (CNC) system effectively and improves the efficiency for the numerical control (NC) machining of sculptured surfaces. The simulation shows that the proposed algorithm is feasible and practical. This algorithm can also be applied to the machining of the whole surface.展开更多
基金supported by National Natural Science Foundation of China under Grants 42192531 and 42192534the Special Fund of Hubei Luojia Laboratory(China)under Grant 220100001the Natural Science Foundation of Hubei Province for Distinguished Young Scholars(China)under Grant 2022CFA090。
文摘The dynamic optimal interpolation(DOI)method is a technique based on quasi-geostrophic dynamics for merging multi-satellite altimeter along-track observations to generate gridded absolute dynamic topography(ADT).Compared with the linear optimal interpolation(LOI)method,the DOI method can improve the accuracy of gridded ADT locally but with low computational efficiency.Consequently,considering both computational efficiency and accuracy,the DOI method is more suitable to be used only for regional applications.In this study,we propose to evaluate the suitable region for applying the DOI method based on the correlation between the absolute value of the Jacobian operator of the geostrophic stream function and the improvement achieved by the DOI method.After verifying the LOI and DOI methods,the suitable region was investigated in three typical areas:the Gulf Stream(25°N-50°N,55°W-80°W),the Japanese Kuroshio(25°N-45°N,135°E-155°E),and the South China Sea(5°N-25°N,100°E-125°E).We propose to use the DOI method only in regions outside the equatorial region and where the absolute value of the Jacobian operator of the geostrophic stream function is higher than1×10^(-11).
文摘This paper introduces the definition of the Orthogonal Type Node Configuration and discusses the corresponding multivariate Lagrange, Hermite and Birkhoff interpolation problems in high dimensional space R s(s>2). This node configuration can be considered to be a kind of extension of the Cross Type Node Configuration , in R 2 to high dimensional spaces. And the Mixed Type Node Configuration in R s(s>2) is also discussed in this paper in an example.
基金The Grant (11RC05) of Scienti/fic Research Foundation for Talents of Hefei Universitythe Grant (11KY06ZR) of Scientific Research Foundation Hefei University+1 种基金the Key Project Foundation (KJ2008A027) of the Department of Education of Anhui Provincethe Project Foundation (KJ2010B182,KJ2011B152, KJ2011B137) of the Department of Education of Anhui Province
文摘We construct general structures of one and two variable interpolation function, without depending on the existence of divided difference or inverse differences, and we also discuss the block based osculatory interpolation in one variable case. Clearly, our method offers marly flexible interpolation schemes for choices. Error terms for the interpolation are determined and numerical examples are given to show the effectlveness of the results.
基金Project supported by the Scientific Research Fund of Hunan Provincial Education Department,China (Grant No.21A0470)the Natural Science Foundation of Hunan Province,China (Grant No.2023JJ50268)+1 种基金the National Natural Science Foundation of China (Grant Nos.62172268 and 62302289)the Shanghai Science and Technology Project,China (Grant Nos.21JC1402800 and 23YF1416200)。
文摘As a branch of quantum image processing,quantum image scaling has been widely studied.However,most of the existing quantum image scaling algorithms are based on nearest-neighbor interpolation and bilinear interpolation,the quantum version of bicubic interpolation has not yet been studied.In this work,we present the first quantum image scaling scheme for bicubic interpolation based on the novel enhanced quantum representation(NEQR).Our scheme can realize synchronous enlargement and reduction of the image with the size of 2^(n)×2^(n) by integral multiple.Firstly,the image is represented by NEQR and the original image coordinates are obtained through multiple CNOT modules.Then,16 neighborhood pixels are obtained by quantum operation circuits,and the corresponding weights of these pixels are calculated by quantum arithmetic modules.Finally,a quantum matrix operation,instead of a classical convolution operation,is used to realize the sum of convolution of these pixels.Through simulation experiments and complexity analysis,we demonstrate that our scheme achieves exponential speedup over the classical bicubic interpolation algorithm,and has better effect than the quantum version of bilinear interpolation.
基金supported by the National Natural Science Foundation of China (No.12172154)the 111 Project (No.B14044)+1 种基金the Natural Science Foundation of Gansu Province (No.23JRRA1035)the Natural Science Foundation of Anhui University of Finance and Economics (No.ACKYC20043).
文摘In this study,a wavelet multi-resolution interpolation Galerkin method(WMIGM)is proposed to solve linear singularly perturbed boundary value problems.Unlike conventional wavelet schemes,the proposed algorithm can be readily extended to special node generation techniques,such as the Shishkin node.Such a wavelet method allows a high degree of local refinement of the nodal distribution to efficiently capture localized steep gradients.All the shape functions possess the Kronecker delta property,making the imposition of boundary conditions as easy as that in the finite element method.Four numerical examples are studied to demonstrate the validity and accuracy of the proposedwavelet method.The results showthat the use ofmodified Shishkin nodes can significantly reduce numerical oscillation near the boundary layer.Compared with many other methods,the proposed method possesses satisfactory accuracy and efficiency.The theoretical and numerical results demonstrate that the order of theε-uniform convergence of this wavelet method can reach 5.
基金This research was funded by the National Nature Sciences Foundation of China(Grant No.42250410321).
文摘Missing value is one of the main factors that cause dirty data.Without high-quality data,there will be no reliable analysis results and precise decision-making.Therefore,the data warehouse needs to integrate high-quality data consistently.In the power system,the electricity consumption data of some large users cannot be normally collected resulting in missing data,which affects the calculation of power supply and eventually leads to a large error in the daily power line loss rate.For the problem of missing electricity consumption data,this study proposes a group method of data handling(GMDH)based data interpolation method in distribution power networks and applies it in the analysis of actually collected electricity data.First,the dependent and independent variables are defined from the original data,and the upper and lower limits of missing values are determined according to prior knowledge or existing data information.All missing data are randomly interpolated within the upper and lower limits.Then,the GMDH network is established to obtain the optimal complexity model,which is used to predict the missing data to replace the last imputed electricity consumption data.At last,this process is implemented iteratively until the missing values do not change.Under a relatively small noise level(α=0.25),the proposed approach achieves a maximum error of no more than 0.605%.Experimental findings demonstrate the efficacy and feasibility of the proposed approach,which realizes the transformation from incomplete data to complete data.Also,this proposed data interpolation approach provides a strong basis for the electricity theft diagnosis and metering fault analysis of electricity enterprises.
文摘The refinement of the 1cm geoid holds significant importance in engineering applications,where the accuracy of the geoid is frequently constrained by its resolution.However,there has been limited exploration into the specific relationship between geoid resolution and accuracy.This article aims to address this gap by thoroughly examining said relationship.This study employs the CapRCR modification to calculate the geoid of Gongzhuling City.The findings indicate that the accuracy can be enhanced by 2%to 9%through encryption of the geoid resolution from 2.5′×2.5′to 1.5′×1.5′.Furthermore,this improvement can be augmented by 15%to 21%through encryption of the gravity anomaly to the same resolution.The accuracy of the geoid exhibits a linear relationship with the resolution of gravity anomalies.The theoretical accuracy of the geoid,excluding integration errors,is determined to be 1.21 cm.In engineering applications,opting for a resolution of 30''×30''can result in the geoid containing only an integration error of 2 mm.However,to attain an accuracy of 1 cm,it becomes imperative to further mitigate data errors.Consequently,the final refined quasi-geoid accuracy is established at 0.56 mm.
文摘The bivariate interpolation in two dimensional space R2 is more complicated than that in one dimensional space R, because there is no Haar space of continuous functions in R2. Therefore, the bivariate interpolation has not a unique solution for a set of arbitrary distinct pairwise points. In this work, we suggest a type of basis which depends on the points such that the bivariate interpolation has the unique solution for any set of distinct pairwise points. In this case, the matrix of bivariate interpolation has the semi inherited factorization.
文摘The accuracy of interpolation models applied to groundwater depends, among other factors, on the interpolation method chosen. Therefore, it is necessary to compare different approaches. For this, different methods of interpolation of nitrate concentrations were contrasted in sixty-seven wells in an aquifer in Aguascalientes, Mexico. Four general interpolation methods were used in ArcGIS 10.5 to make the maps: IDW, Kriging, Natural Neighbor and Spline. In the modeling, only method type was varied. The input parameters (location, temporality, and nitrate concentration) were the same in the four interpolations;despite this, different maximum and minimum values were obtained for each interpolation method: for IDW, 0.2 to 22.0 mg/l, for Kriging, 3.5 to 16.5 mg/l, for Natural Neighbor, 0.3 to 21.7 mg/l and for Spline −30.8 to 37.2 mg/l. Finally, an assessment of the maps obtained was conducted by comparing them with the Official Mexican Standard (OMS), where 24 of the 67 wells were found outside the 10 mg/l that the OMS establishes as maximum permissible limit for human consumption. Taking as a starting point the measured values of nitrates (0.25 to 22.12 mg/l), as well as the spatial distribution of the interpolated values, it was determined that the Krigging method best fitted the data measured in the wells within the studied aquifer.
文摘The purpose of this paper is to investigate the spatial interpolation of rainfall variability with deterministic and geostatic inspections in the Prefecture of Kilkis (Greece). The precipitation data where recorded from 12 meteorological stations in the Prefecture of Kilkis for 36 hydrological years (1973-2008). The cumulative monthly values of rainfall were studied on an annual and seasonal basis as well as during the arid-dry season. In the deterministic tests, the I.D.W. and R.B.F. checks were inspected, while in the geostatic tests, Ordinary Kriging and Universal Kriging respectively. The selection of the optimum method was made based on the least Root Mean Square Error (R.M.S.E.), as well as on the Mean Error (M.E.), as assessed by the cross validation analysis. The geostatical Kriging also considered the impact of isotropy and anisotropy across all time periods of data collection. Moreover, for Universal Kriging, the study explored spherical, exponential and Gaussian models in various combinations. Geostatistical techniques consistently demonstrated greater reliability than deterministic techniques across all time periods of data collection. Specifically, during the annual period, anisotropy was the prevailing characteristic in geostatistical techniques. Moreover, the results for the irrigation and seasonal periods were generally comparable, with few exceptions where isotropic methods yielded lower (R.M.S.E.) in some seasonal observations.
基金supported by the grant of Key Scientific Research Foundation of Education Department of Anhui Province, No. KJ2014A210
文摘General interpolation formulae for barycentric interpolation and barycen- tric rational Hermite interpolation are established by introducing multiple parameters, which include many kinds of barycentric interpolation and barycentric rational Her- mite interpolation. We discussed the interpolation theorem, dual interpolation and special cases. Numerical example is given to show the effectiveness of the method.
文摘In this paper, we obtain a properly posed set of nodes for interpolation on a sphere. Moreover it is applied to construct properly posed set of nodes for Lagrange interpolation on the trivariate polynomial space of total degree n.
基金supported by the National Science and Technology Major Projects(No.2011ZX05020-008)Well Logging Advanced Technique and Application Basis Research Project of Petrochina Company(No.2011A-3901)
文摘The radial basis function (RBF) interpolation approach proposed by Freedman is used to solve inverse problems encountered in well-logging and other petrophysical issues. The approach is to predict petrophysical properties in the laboratory on the basis of physical rock datasets, which include the formation factor, viscosity, permeability, and molecular composition. However, this approach does not consider the effect of spatial distribution of the calibration data on the interpolation result. This study proposes a new RBF interpolation approach based on the Freedman's RBF interpolation approach, by which the unit basis functions are uniformly populated in the space domain. The inverse results of the two approaches are comparatively analyzed by using our datasets. We determine that although the interpolation effects of the two approaches are equivalent, the new approach is more flexible and beneficial for reducing the number of basis functions when the database is large, resulting in simplification of the interpolation function expression. However, the predicted results of the central data are not sufficiently satisfied when the data clusters are far apart.
基金Supported by Forestry Science and Technology Support Project (2008BADB0B0203)National Technology Support Project (2007BAC03A08-5)
文摘In order to get the spatial grid data of monthly precipitation and monthly average temperature of Sanjiangyuan area, the Co-Kriging (COK) and thin plate smoothing splines(TPS) interpolation methods were applied by using the climate data during 1971-2000 of 58 meteorological stations around Qinghai Province and the 3 arc-second digital elevation model (DEM) data. The performance was evaluated by the smallest statistical errors by general cross validation (GCV). Root-mean-squared predicted errors (RMSE) and mean absolute errors (MAE) were used to compare the performance of the two methods. The results showed that: 1) After combing covariates into the models, both methods performed better; 2) The performance of TPS was significantly better than COK: for monthly average temperature, the RMSE derived from TPS was 69.48% higher than COK, as MAE increased by 70.56%. And for monthly precipitation, the RMSE derived from TPS was 28.07% higher than COK, as MAE increased by 29.06%.
基金The Natural Science Foundation of Anhui Province(No.1308085MF96)the Project of Chuzhou University(No.2012qd06,2011kj010B)+1 种基金the Scientific Research Foundation of Education Department of Anhui Province(No.KJ2014A186)the National Basic Research Program of China(973 Program)(No.2010CB732503)
文摘To alleviate the distortion of XRII X-ray image intensifier images in the C-arm CT computer tomography imaging system an algorithm based on the Delaunay triangulation interpolation is proposed.First the causes of the phenomenon the classical correction algorithms and the Delaunay triangulation interpolation are analyzed.Then the algorithm procedure is explained using flow charts and illustrations. Finally experiments are described to demonstrate its effectiveness and feasibility. Experimental results demonstrate that the Delaunay triangulation interpolation can have the following effects.In the case of the same center the root mean square distances RMSD and standard deviation STD between the corrected image with Delaunay triangulation interpolation and the ideal image are 5.760 4 ×10 -14 and 5.354 2 ×10 -14 respectively.They increase to 1.790 3 2.388 8 2.338 8 and 1.262 0 1.268 1 1.202 6 after applying the quartic polynomial model L1 and model L2 to the distorted images respectively.The RMSDs and STDs between the corrected image with the Delaunay triangulation interpolation and the ideal image are 2.489 × 10 -13 and 2.449 8 ×10 -13 when their centers do not coincide. When the quartic polynomial model L1 and model L2 are applied to the distorted images they are 1.770 3 2.388 8 2.338 8 and 1.269 9 1.268 1 1.202 6 respectively.
文摘Aim To study wind tunnel test data interpolation methods for flight vehicle with aerodynamic axial asymmetry. Methods For different body aerodynamic roll angles, proper wind tunnel test schemes were selected and trigonometric series were used for aerodynamic interpolation. Results and Conclusion A simple and effective scheme for wind tunnel test and an accurate aerodynamic interpolation method are developed with satisfactory results.
基金The National Natural Science Foundation of China(No.61261007,61062005)the Key Program of Yunnan Natural Science Foundation(No.2013FA008)
文摘A novel spatial interpolation method based on integrated radial basis function artificial neural networks (IRBFANNs) is proposed to provide accurate and stable predictions of heavy metals concentrations in soil at un- sampled sites in a mountain region. The IRBFANNs hybridize the advantages of the artificial neural networks and the neural networks integration approach. Three experimental projects under different sampling densities are carried out to study the performance of the proposed IRBFANNs-based interpolation method. This novel method is compared with six peer spatial interpolation methods based on the root mean square error and visual evaluation of the distribution maps of Mn elements. The experimental results show that the proposed method performs better in accuracy and stability. Moreover, the proposed method can provide more details in the spatial distribution maps than the compared interpolation methods in the cases of sparse sampling density.
基金sponsored by:the National Basic Research Program of China (973 Program) (2007CB209605)the National Natural Science Foundation of China (40974073)the National Hi-tech Research and Development Program of China (863 Program) (2009AA06Z206)
文摘In this research, we present a seismic trace interpolation method which uses seismic data with surface-related multiples. It is different from conventional seismic data interpolation using information transformation or extrapolation of adjacent channels for reconstruction of missing seismic data. In this method there are two steps, first, we construct pseudo-primaries by cross-correlation of surface multiple data to extract the missing near- offset information in multiples, which are not displayed in the acquired seismic record. Second, we correct the pseudo-primaries by applying a Least-squares Matching Filter (LMF) and RMS amplitude correction method in time and space sliding windows. Then the corrected pseudo-primaries can be used to fill the data gaps. The method is easy to implement, without the need to separate multiples and primaries. It extracts the seismic information contained by multiples for filling missing traces. The method is suitable for seismic data with surfacerelated multiples.
基金The Natural Science Foundation of Jiangsu Province(No.BK2003005).
文摘A direct interpolation algorithm for spatial curves on a surface is proposed for pen-cutting of sculptured surfaces. The algorithm can carry out the direct interpolation for projective curves lying on the sculptured surface. It is based on the geometric and kinetic relationships between drive curves and cutter-contact (C-C) curves. It evaluates the parameter of drive curves corresponding to interpolation points by the Taylor formula. Then it gains coordinates of interpolation points indirectly by inverse calculation. Finally, it generates the motion commands for machine tool controller. This method extends the locus-controlled function in the computer numerical control (CNC) system effectively and improves the efficiency for the numerical control (NC) machining of sculptured surfaces. The simulation shows that the proposed algorithm is feasible and practical. This algorithm can also be applied to the machining of the whole surface.