Abstract
Standard Offset surfaces are defined as locus of the points which are at constant distance along the unit normal direction from the generator surfaces. Offset are widely used in various practical applications, such as tolerance analysis, geometric optics and robot path-planning. In some of the engineering applications, we need to extend the concept of standard offset to the generalized offset where distance offset is not necessarily constant and offset direction are not necessarily along the normal direction. Normally, a generalized offset is functionally more complex than its progenitor because of the square root appears in the expression of the unit normal vector. For this, an approximation method of its construction is necessary. In many situation it is necessary to fill or reconstruct certain function defined in a domain in which there is a lack of information inside one or several sub-domains (holes). In some practical cases, we may have some specific geometrical constrains, of industrial or design type, for example, the case of a specified volume inside each one of these holes. The problem of filling holes or completing a 3D surface arises in all sorts of computational graphics areas, like CAGD, CAD-CAM, Earth Sciences, computer vision in robotics, image reconstruction from satellite and radar information, etc. In this work we present an approximation method of filling holes of the generalized offset of a surface when there is a lack information in a sub-domain of the function that define it. We prove the existence and uniqueness of solution of this problem, we show how to compute it and we establish a convergence result of this approximation method. Finally, we give some graphical and numerical examples.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
1 Introduction
In the field of computer aided geometric design (CAGD), offset curves and surfaces have got considerable attention since they are widely used in various practical applications such as tolerance analysis, geometric optics and robot path-planning.
This work can be considered as a continuation or a note of our recently published work [2]. Then, for some applications of this kind of problem in the domain of chemical sciences, one can consult such applications in [2], for example the authors in [1] investigate the electrocoagulation/flotation (ECF) treatment efficiency for the removal of copper, turbidity and organic substances from the waste offset printing developer (WOPD). Meanwhile, the authors in [17] show the ability to recover the subsurface of surface-enhanced Raman scattering (SERS) signal. In [19], the authors reviews the interplay of weak non-covalent interactions involved in the formation of self-assembled mono-layers of organic molecules and the strong chemical binding in directed-assembly of organic molecules on solid surfaces.
The authors in [4] have elucidated the relationship between the coverage of cadmium carboxylate on CdSe surfaces and the trapping dynamics. Reducing the metal enrichment increases the rate of hole trapping and reduces the PLQY. In [16] the authors have employed transient absorption spectroscopy and transient photocurrent measurements to characterize the dynamics of photogenerated holes in CoPi-modified BiVO\(_4\) photoanodes. The competition between water oxidation and electron/hole recombination was quantified using a simple kinetic model. In [6] transient absorption spectroscopy was used to analyse the trapping of electrons and holes generated by band gap excitation of NiO nanoparticle films. The analysis revealed that after a fs pulse the hole is trapped on a sub-picosecond time scale in a “Ni\(^{4+}\)” state. Meanwhile in [13] the authors have presented direct time-resolved spectroscopic evidence for quasi-type-II band alignment in graded alloy CdS\(_x\)Se\(_{1-x}\) NCs through analysis of the ultrafast charge carrier dynamics as a function of chemical composition. The fast \(\thicksim 3\) ps hole-trapping process disappears as sulfur composition increases to form a graded CdS-rich shell with VB offset such that excited holes are confined to the CdSe-rich core. The authors in [12] demonstrate that measuring the sub-30-pm displacements of atoms from high-symmetry positions in the atomically resolved scanning tunnelling microscopy allows the physical order parameter fields to be visualized in real space on the single-atom level. This local crystallographic analysis is applied to the in-situ-grown manganite surfaces. In particular, using direct bond-angle mapping the authors report direct observation of structural domains on manganite surfaces, and trace their origin to surface-chemistry-induced stabilization of ordered Jahn-Teller displacements. In recent decades, the alkali metals have been measured by X-ray photoelectron spectroscopy (XPS), which is a powerful technique to detect elemental compositions of the surface of the material and determine the related chemical states. In [21] the XPS profiles of the alkali metals are analyzed via the ZPS analysis method, the BOLS notation, tight-binding (TB) theory, and density functional theory (DFT) method, which not only facilitates the reproduction of thermal and size trend of energy changes, but also clarifies the physical origin of their core-electron BE shifts in alkali metals.
Offset curves and surfaces are defined as locus of the points which distance along the normal direction from the original ones.
The parametric representation of the offset of a given parametric surface \(\varvec{s}(u,v)\) is another parametric surface given by
where d is the signed offset distance and \(\varvec{n}(u,v)\) is the unit normal vector at the point \(\varvec{s}(u,v)\).
Except a very surfaces as planes, spheres, cones, torus or cylinders, offset surfaces cannot be expressed in the same forms as their original ones. Therefore, approximation methods are needed.
In some of the engendering applications, we need to extend the concept of the standard offset surface, such as geodesic offset, where the constant distance is replaced by geodesic distance, or generalized offset, where the offset distance is not necessarily along the normal direction.
Generalized offset curves and surfaces were first introduced by Bréchner[3] and have been extended further, from the differential geometric as well as algebraic point of view by Pottman[18].
For a regular parametric surface \(\varvec{r}(u,v)=(x(u,v),y(u,v),z(u,v))\), its two unit tangent vectors in the directions of u and v and its normal vector are given by
respectively, where \(\langle \,\cdot \rangle _3\) is the usual Euclidean norm in \(\mathbb {R}^3\).
In this situation, the generalized surface \(\varvec{r}^0(u,v)\) with the variable offset distance and direction determined by the vectors \(d_1(u,v)\varvec{e}_1(u,v)\), \(d_2(u,v)\varvec{e}_2(u,v)\) and \(d_3(u,v)\varvec{n}(u,v)\) is defined by
Hole-filling techniques appear in many different real applications, like surface reconstruction in engineering [20], 3D human body scanning, dental reconstruction, reverse engineering, etc. In 3D scanning applications for example, data can be missing due to accessibility limitations, occlusion, reflecting spaces or surfaces perpendicular to camera/scanner.
In [11] the procedure considered presents a surface reconstruction and hole-filling scheme, based on a network of parallel and/or orthogonal curves, called also wireframes, where previously some specific one-dimensional curve reconstructions have been accomplished inside the hole that they want to fill. Under this perspective, the authors face the problem of filling surface by mixing the models presented in some of their references in an efficient way. Instead of addressing the problem as a whole, which yields a very complicated computational frame, the authors generate the reconstructed filling surface through one-dimensional filling curves subsequently used to generate the global filling smooth surface. The method proposed is developed just for explicit surface.
In [10] the methods and the procedures studied manage to find a smooth bivariate spline function that minimizes certain functional, and that properly takes into account all the features that the authors want this approximate spline to fulfill, this means technique of minimizing functionals in order to obtain splines verifying certain required conditions has been used in the literature in the last years. Among these possible conditions, there are the volume restrictions, that can be included as an interpolation condition inside the corresponding finite element space constructed, or a just as another approximation term to be included in the associated quadratic functional. Also, the authors has included some other approximation or interpolation conditions at some specific points outside the holes.
In [9] the authors propose to generalize the filling method previously developed in other works in order to fill holes with some shape conditions, i.e. in such a way that the filling path “inherits” as much as possible the shape of the original surface where it is known. In other words, in real single-variable calculus, there are definitions for increasing, decreasing, concave and convex functions. These definitions involve derivatives and make reference to shape characteristics of the functions. Following this idea, given a data function f with a hole, the method proposed in their paper consists of estimating the shape of f inside the hole by constructing functions that estimate the unknown derivatives of f whose derivatives inside the hole be as close as possible to the estimates ones of f, pretending in this way that the shapes of the function f and its filing path be close.
In [8] the authors treat the problem of adequately filling the holes of 3D-surfaces, not necessarily explicit, and even closed in parametric form, as spheroid type ones: that is, such obtained by a radial function of the spherical coordinates.
In [7] the authors develop two different approaches of a method to fill polygonal holes in a given surfaces by using smoothing variational splines: discontinuous filling and regular filling. In the discontinuous case, they fill the holes with spline functions in a finite element space that minimizes an energy functional. Such filling are chosen to be smooth and as close as possible to the original surface in the neighborhoods of the hole. In this approach the global reconstructed surface will not be continuous. On the contrary, in the regular case, they do not only fill the holes, but also they replace the known surface with another very similar in such a way that the global reconstructed surface will have the desired regularity.
In some practical cases we also know some specific geometrical constrains, of industrial and design type, as the special case of a specified volume inside each on of this sub-domains. The studied method in this work manages to find a function of a biquadratic spline space that minimizes certain quadratic functional that includes the usual semi-norm of order two in a Sobolev space. This minimizing functional technique has proven to be effective on some approximation methods in different functional spaces (see [14, 15] and the references therein).
We study here the problem in which we have surface points (e.g. from a given function, or from some data acquisition procedure) where there is a lack of information inside some sub-domains. As we said this work is a continuation of [2] let’s give a very clear observation between the two works, in fact first this work it is about filling holes and second the approximation is made with the biquadratic splines functions that are of class 1, while in [2] it is done with the bi-cubic splines that are of class 2. This clearly shows that the degree of approximation in this work has been significantly reduced.
The remainder of the manuscript is organized as follows. In sect. 2 we give some necessary formulations and preliminaries. In sect. 3 we present an approximation method of the generalized offset surface with the hole, in the first step we formulate the problem, in the second step we show how to compute such surface and in the last step we prove a convergence result. Section 4 is devoted to study an approximation method of the filling hole, in the first step we formulate the problem, in the second step we show a computation result. In sect. 5 we define the smoothing parametric biquadratic spline filling the hole. Finally, we present some numeric and graphical examples in order to prove the effectiveness and the useful of our method.
2 Notations and preliminaries
We denote by \(\langle \,\cdot \,\rangle _k\) and \(\langle \, , \, \rangle _k\) respectively, the Euclidean norm and the inner product in \(\mathbb {R}^r\).
For any interval (a, b) ((c, d)) with \(a<b\) (\(c<d\)) we consider the rectangle \(R=(a,b)\times (c,d)\).
For any set \(\omega \subset \mathbb {R}^2\) let \(H^3(\omega ;\mathbb {R}^3)\) be the usual Sobolev space of (classes of) functions \(\varvec{u}\) belong to \(L^2(\omega ;\mathbb {R}^3)\), together with their partial derivatives \(D^{\varvec{\beta }}\varvec{u}\), with \(\varvec{\beta }=(\beta _1,\beta _2)\), in the distribution sense, of order \(\vert \varvec{\beta }\vert =\beta _1+\beta _2\le 3\). This space is equipped with the inner products
the semi-norms \(\displaystyle \vert \varvec{u}\vert _{\ell ,\omega }=(\varvec{u},\varvec{u})_\ell ^\frac{1}{2}\), for \(0\le \ell \le 3\), and the norm \(\displaystyle \Vert \varvec{u}\Vert _\omega =\left( \sum _{\ell =0}^3 \vert \varvec{u}\vert _{\ell ,\omega }^2\right) ^\frac{1}{2}\).
Let \(T_n=\{a=x_0<x_1<\cdots <x_n=b\}\) and \(T_m=\{c=y_0<y_1<\cdots <y_m=d\}\) be two partitions of [a, b] and [c, d] respectively, and \(S_2(T_n)\), \(S_2(T_m)\) the spaces of quadratic spline functions of class \(C^1\) constructed from \(T_n\) and \(T_m\) respectively. It is verified that \(\dim (T_n)=n+2\) and \(\dim (T_m)=m+2\)
Let \(\{\phi _1(x),\ldots ,\phi _{n+2}(x)\}\) and \(\{\psi _1(y),\ldots ,\psi _{m+2}(y)\}\) be respectively the quadratic B-spline bases functions of \(S_2(T_n)\) and \(S_2(T_m)\), we consider the space of biquadratic functions
Let \(M=\dim S_2(T_n,T_m)=(n+2)(m+2)\).
For \(i=1,\ldots ,n+2\), \(j=1,\ldots ,m+2\) and for any \((x,y)\in R\), we define \(B_{(m+2)(i-1)+j}(x,y)=\phi _i(x)\psi _j(y)\). Then \(\{B_1,\ldots ,B_M\}\) is the biquadratic B-spline basis functions of \(S_2(T_n,T_m)\).
Being \(\displaystyle h=\max \{\frac{b-a}{n}, \frac{d-c}{m}\}\), let \(V_h=(S_2(T_n,T_m))^3\) be the parametric space of biquadratic splines functions. It is verified that \(V_h\subset H^2(R;\mathbb {R}^3)\cap C^1(\bar{R};\mathbb {R}^3)\).
Let \(\mathbb {R}^{k,3}\) be the space of real matrices of k rows and 3 columns equipped with the inner product \(\displaystyle \langle A,B\rangle _{k,3}=\sum _{i=1}^k\sum _{j=1}^3 a_{ij} b_{ij}\), with \(A=\mathop {(a_{ij})}\nolimits _{\begin{array}{c} 1\le i\le k\\ 1\le j\le 3 \end{array}}\), \(B=\mathop {(b_{ij})}\nolimits _{\begin{array}{c} 1\le i\le k\\ 1\le j\le 3 \end{array}}\) and the corresponding norm \(\displaystyle \langle A\rangle _{k,3}=\langle A,A\rangle _{k,3}^{\frac{1}{2}}\).
Let \(\mathcal {T}_h\) be the set of partition rectangles \(T_n\times T_m\) of R, we consider a closed set \(H\subset R\) (the hole) and let \(\displaystyle H_h=\mathop {\bigcup }\limits _{\begin{array}{c} K\in \mathcal {T}_h\\ K\cap H\ne \emptyset \end{array}} K.\)
Now, given a surface \(\mathcal {S}\) parameterized by the function \(\varvec{f}:\; R\rightarrow \mathbb {R}^3\) and its generalized offset surface \(\mathcal {O}\) parameterized by the function
for all \((x,y)\in R.\) Suppose we can know the points of \(\mathcal {S}\) by the parametrization \(\varvec{f}\) only over \(R-H\).
We wish to approximate the surface \(\mathcal {O}\) parameterized by \(\varvec{s}_f\) from a finite set of points of \(\varvec{f} (R-H)\) and filling the holes by a constructed function of type
such that \(\varvec{\sigma }^h\in V_h\).
3 Approximating \(\varvec{s}_f\vert _{R-H_h}\)
3.1 Formulating of the problem
For each \(N\in \mathbb {N}^\star\), let \(A_N=\{\varvec{a}_1,\ldots ,\varvec{a}_k\}\) be a finite set of \(k=k(N)\) points of \(R-H\) such that
Hence, we can obtain, for N sufficiently large, that \(k(N)>N\).
Now, let \(d_1(u,v)\varvec{e}_1(u,v)\), \(d_2(u,v)\varvec{e}_2(u,v)\) and \(d_3(u,v)\varvec{n}(u,v)\) be some given offset variables distances and directions.
Let \(\rho :\, H^2(R;\mathbb {R}^3)\rightarrow \mathbb {R}^{k,3}\) be the Lagrangian operator defined by \(\rho \varvec{v}=(v(\varvec{a}_i))_{1\le i\le k}\) and suppose that
where \(\mathbb {P}_2(R;\mathbb {R}^3)\) is the linear space of the restrictions to R of the bivariate polynomials of degree less than or equal to 2 with vectorial coefficients in \(\mathbb {R}^3\).
For any real number \(\varepsilon >0\), let \(J_1:\, H^2(R;\mathbb {R}^3)\rightarrow \mathbb {R}\) be the functional given by
Problem 1 We consider the following minimization problem: Find \(\varvec{\sigma }_1^h\in V_h\) such that
Proposition 1
Problem (6) has a unique solution, called the generalized offset biquadratic spline in \(V_h\) relative to N, \(A_N\), \(\rho (\varvec{s}_{\varvec{f}})\), \(d_1\), \(d_2, d_3\) and \(\varepsilon\), that is characterized as the unique solution of the following variational problem: Find \(\varvec{\sigma }_1^h\in V_h\) such that
Proof
Let define the following application
Using (4) we obtain that the application \(\, [[\,\,\cdot \,\,]]\,\) defines a norm in \(V_h\) which is equivalent to the usual norm \(\Vert \,\cdot \,\Vert _R\). Meanwhile, it is easily to prove that the symmetric bilinear form \(\tilde{a}:\, V_h\times V_h \rightarrow \mathbb {R}\) defined by
is continuous and \(V_h\)-elliptic.
Then, the result is obtained applying Lax-Milgram Lemma (see [5]) taking into account that the linear application \(\Phi : H^2(R;\mathbb {R}^3) \rightarrow \mathbb {R}\) defined by
is continuous. \(\square\)
3.2 Computing the solution \(\varvec{\sigma }_1^h\)
Now, we are going to see how to compute in practice the generalized offset biquadratic spline \(\varvec{\sigma }_1^h\in V_h\) relative to N, \(A_N\), \(\rho (\varvec{s}_{\varvec{f}})\), \(d_1\), \(d_2\), \(d_3\) and \(\varepsilon\). Let \(\{\omega _1,\ldots ,\omega _M\}\) be the biquadratic B-spline basis functions of \(S_2(T_n,T_m)\) and \(\{\varvec{j}_1,\varvec{j}_2,\varvec{j}_3\}\) the canonical one of \(\mathbb {R}^3\). We define
Then \(\{\varvec{v}_1,\ldots ,\varvec{v}_{3M}\}\) is a basis of \(V_h\).
Thus \(\displaystyle \varvec{\sigma }_1^h=\sum _{i=1}^{3M} \alpha _i \varvec{v}_i\), with \(\alpha _1,\ldots .\alpha _{3M}\in \mathbb {R}\) are the solution of the linear system
where
Remark 1
We can observe that for computing the solution of the linear system (8) only certain values of \(\varvec{s}_{\varvec{f}}\) on \(R-H\) are necessary.
3.3 Convergence result
The objective of this subsection is to prove, under adequate conditions, that the generalized offset biquadratic spline in \(V_h\) relative to N, \(A_N\), \(\rho (\varvec{s}_{\varvec{f}})\), \(d_1\), \(d_2\), \(d_3\) and \(\varepsilon\) converges to \(\varvec{s}_{\varvec{f}}\) locally in \(R-H\) as N tends to \(+\infty\) and h tends to 0.
For this, we can reason as [2, sect. 5] by using the unique parametric biquadratic spline function \(\varvec{S}_M\in V_h\) which solves the interpolation problem
and considering its restriction \(\left. \varvec{S}_M\right| _{R-H}\).
The interpolation function \(\varvec{S}_M\) verifies that there exist a positive real constant C such that
Hence, reasoning as in [2, Theorem 7] we can prove the following result.
Theorem 1
Suppose the hypotheses (4) and the following ones
and
hold. Then, for any \(\varvec{x}\in R-H\) there exists an open subset \(\omega _{\varvec{x}}\subset R-H\) such that \(\varvec{x}\in \omega _{\varvec{x}}\) and
Proof
From (9), (10) one has that \(h\rightarrow 0\) as \(N\rightarrow +\infty\). Then, given \(\varvec{x}\in R-H\), there exists \(h_1>0\) such that, for \(h\le h_1\), \(\varvec{x}\in R-H_h\) and thus there exists an open subset \(\omega _{\varvec{x}} \subset R-H\) such that \(\varvec{x}\in \omega _{\varvec{x}}\) and \(\omega _{\varvec{x}}\subset R-H_h\), for all \(h\le h_1\).
From here, and reasoning as in [2, Theorem 7] we can obtained the result.\(\square\)
4 Filling holes
4.1 Formulating the problem
We denote by \(B_h=A_N\cap H_h=\{\varvec{b}_1,\ldots ,\varvec{b}_{M_1}\}\), \(C_h^N=\{\varvec{c}_1,\ldots , \varvec{c}_{M_2} \}\) the set of all the knots of \(T_n\times T_m\) that they are on some side of the boundary of \(H_h\) and \(D_h^N=\{ \varvec{d}_1,\ldots ,\varvec{d}_{M_3}\}\) the set of all knots of \(C_h^N\) that are not vertices.
For \(i=1,2,3\) we consider the operators \(\rho _i^h: H^2(R;\mathbb {R}^3)\rightarrow \mathbb {R}\) given by
where \(\displaystyle \frac{\partial \varvec{v}}{\partial \varvec{n}}\) indicates the normal derivative of the restriction of \(\varvec{v}\) to \(H_h\).
Let \(V_h^*\) be the restriction set of the functions of \(V_h\) to \(H_h\) and let consider the subsets
Now, let \(J_2^h: X_h\rightarrow \mathbb {R}\) be the functional given by
Problem 2 We consider the following minimization problem: Find \(\varvec{\sigma }_2^h\in X_h\) such that
Proposition 2
Problem (11) has a unique solution that is characterized as the unique solution of the following variational problem: Find \(\varvec{\sigma }_2^h\in X_h\) such that
Proof
We consider the application \(\tilde{a}: X_h\times X_h\rightarrow \mathbb {R}\) given by
Obviously this is a bilinear and symmetric form in \(X_h\). It is easily to prove that \(\tilde{a}\) is continuous and \(X_h\)-coercive taking into account that the application
defines a norm in \(X_h\) which is equivalent to the usual norm \(\Vert \,\cdot \,\Vert _{H_h}\).
Let \(\varphi (\varvec{v})=2\left( \langle \rho _1^h(\varvec{v}),\rho _1^h(\varvec{\sigma }_1^h)\rangle _{M_1,3}\right)\), which is clearly a continuous and linear application.
So, by applying Stampacchia Theorem (see [5]), we conclude that there exists a unique \(\varvec{\sigma }_2^h\in X_h\) such that \(\tilde{a}(\varvec{\sigma }_2^h,\varvec{\omega }-\varvec{\sigma }_2^h)\ge \varphi (\varvec{\omega }-\varvec{\sigma }_2^h)\), for all \(\varvec{\omega }\in X_h\), which implies that \(\tilde{a}(\varvec{\sigma }_2^h,\varvec{v})\ge \varphi (\varvec{v})\), for all \(\varvec{v}\in X_h^0\),
Taking into account that \(X_h^0\) is a linear subspace, if \(\varvec{v}\in X_0^h\) then \(-\varvec{v}\in X_0^h\) and thus \(\tilde{a}(\varvec{\sigma }_2^h,-\varvec{v})\ge \varphi (-\varvec{v})\). From this we obtain that \(\tilde{a}(\varvec{\sigma }_2^h,\varvec{v})=\varphi (\varvec{v})\), for any \(\varvec{v}\in X_h^0\).
Furthermore, \(\varvec{\sigma }_2^h\) is the minimum in \(X_h\) of the functional \(\displaystyle \Phi (\varvec{v})=\frac{1}{2} \tilde{a}(\varvec{v},\varvec{v})-\varphi (\varvec{v})\), which is the minimum of \(J_2\) in \(X_h\) since \(\Phi (\varvec{v})=J_2(\varvec{v})-\langle \rho _1^h(\varvec{\sigma }_1^h)\rangle _{M_1,3}\).
Hence, we conclude the result. \(\square\)
Proposition 3
There exists \((\varvec{\sigma }_2^h,\varvec{\lambda }_1,\varvec{\lambda }_2)\in X_h\times \mathbb {R}^{M_2,3} \times \mathbb {R}^{M_3,3}\) such that
Proof
Let us denote by \(\{\varvec{\omega }_1^h,\ldots ,\varvec{\omega }_{M_2}^h\}\) the basis functions of \(X_h\) associate to the degrees of freedom \(\{v(\varvec{c}_i)\, \vert \, i=1,\ldots , M_2\}\) and by \(\{\varvec{\omega }_{M_2+1}^h,\ldots ,\varvec{\omega }_{M_2+M_3}^h\}\) the basis functions of \(X_h\) associate to the degrees of freedom \(\displaystyle \{\frac{\partial \varvec{v}}{\partial \varvec{n}}(\varvec{d}_{j})\,\vert \, j=1,\ldots ,M_3\}\).
For each \(\varvec{v}\in X_h\), let \(\displaystyle \varvec{\omega }=\varvec{v}-\sum _{j=1}^{M_2} \varvec{v}(\varvec{c}_j) \varvec{\omega }_j^h-\sum _{j=1}^{M_3} \frac{\partial \varvec{v}}{\partial \varvec{n}}(\varvec{d}_j)\varvec{\omega }_{M_2+j}^h\). Then \(\omega \in X_h\). Moreover, for any \(k=1,\ldots ,M_2\),
and, for \(k=1,\ldots M_3\), one has
So \(\rho _2^h(\varvec{\omega })=\rho _3^h(\varvec{\omega })=0\) and consequently \(\varvec{\omega }\in X_h^0\).
Let \(\varvec{\sigma }_2^h\) the solution of (12). Then we have \(\varvec{\sigma }_2^h \in X_h\) and
By linearity we obtain
If we denote by
then we conclude that (12) holds.
The uniqueness of \(\varvec{\lambda }_1\) and \(\varvec{\lambda }_2\) is immediate. \(\square\)
4.2 Computing the solution \(\varvec{\sigma }_2^h\)
From the basis \(\{\omega _1,\ldots ,\omega _M\}\) of \(S_2(T_n,T_m)\), let \(\{\omega _1^h,\ldots ,\omega _\mu ^h\}\) be the set of the functions of this basis such that
and we define
Then \(\{\varvec{v}_1^h,\ldots ,\varvec{v}_\mu ^h\}\) is a basis of \(V_h^*\).
Thus \(\displaystyle \varvec{\sigma }_2^h=\sum _{i=1}^{\mu } \beta _i \varvec{v}_i^h\), where \(\beta _1,\ldots ,\beta _{\mu }\in \mathbb {R}\) are de coefficients of \(\varvec{\sigma }_2^h\).
Applying Proposition 3, by linearity, this coefficients can be obtained from the unique solution of the linear system
where
Remark 2
Observe that for computing the solution of (14) only the values of \(\varvec{\sigma }_1^h\) on \(H_h\) and its normal derivative on the boundary of \(H_h\) are neccesary.
5 The solution of the problem
Definition 1
We define the smoothing parametric biquadratic spline function filling holes of \({\varvec{f}}\vert _{R-H_h}\) in \(V_h\) associated with \(A_N\) and \(\varepsilon\) as the function
By construction it is verified that \(\varvec{\sigma }\in V_h\).
6 Numerical and graphical examples
To test the effectiveness of our method we suppose that the initial surface parameterized by the function \(\varvec{f}:\; R=(0,1)\times (0,1)\rightarrow \mathbb {R}^3\) is given by
We consider a hole H that is the close of the interior of the ellipses \(25(x-0.35)^2+150(y-0.75)^2=1\) and \(150(x-0.75)^2+35(y-0.3)^2=1\) and the surface with these holes parameterized by \({\varvec{f}}\vert _{R-H}\) (see Fig. 1 left side).
Given the variable distances defined by \(d_1(x,y)=0.2\sin (0.5(x-y)^2+0.3)\), \(d_2(x,y)=0.2(\sin (0.5(x-y)^2+0.3)\) and \(d_3(x,y)=-0.2\cos (0.2\cos (0.5(x-y)^2+0.3)\), we consider the corresponding generalized offset surface with holes \(\displaystyle \varvec{s}_{\varvec{f}}\vert _{R-H}\) (see Fig. 1 right side).
First, we start with a graphical example. For this, from a partition \(\mathcal {T}_h\) of \(8\times 8\) equal squares, i.e. \(9\times 9\) knots, we construct the approximated hole \(H_h\) (see Fig. 2 left side) and we take a set \(A_N\) of \(k(N)=500\) random points in the closed set of \(R-H\) (see Fig. 2 right side).
Second, we compute the smoothing parametric biquadratic spline in \(V_h\) associated with \(A_N\), \(\rho ({\varvec{s}}_{\varvec{f}})\) and \(\varepsilon =10^{-9}\), called \(\varvec{\sigma }_1^h\), and an estimation of the relative \(L^2\)-error given by the number
where \(\{\xi _1,\ldots ,\xi _{5000}\}\) is a set of 5000 random points of the closed set of \(R-H_h\).
Fig. 3 shows the original surfaces parameterized by the functions \(\varvec{f}\vert _{R-H_h}\) and \(\varvec{s}_{\varvec{f}}\vert _{R-H_h}\), on the left side, and the approximated generalized offset surfaces parameterized by \(\varvec{\sigma }_1^h\), on the right side, from bottom to top in both cases. The obtained estimated relative error is \(Error_1=2.7695\times 10^{-4}\).
Third, we compute the smoothing parametric biquadratic spline function in \(V_h\vert _H\) associated with \(H_h\), \(\rho _i\varvec{\sigma }_1^h\), for \(i=1,2,3\), and \(\varepsilon =10^{-9}\), denoted by \(\varvec{\sigma }_2^h\).
Fourth, we compute the smoothing parametric biquadratic spline function filling holes of \(\varvec{f}\) in \(V_h\) associated with \(A_N\) and \(\varepsilon =10^{-9}\), denoted by \(\varvec{\sigma }^h\) and an estimation of the relative \(L^2\)-error given by the number
where \(\{\xi _1,\ldots ,\xi _{5000}\}\) is a set of 5000 random points of R.
Fig. 4 shows the original surface with the approximated hole parameterized by \(\varvec{f}\vert _{R-H}\) and its approximated generalized offset surface filling hole, parameterized by the smoothing parametric biquadratic spline filling hole function \(\varvec{\sigma }^h\), from bottom to top. The obtained estimated relative error is \(Error=2.7216\times 10^{-4}\).
Finally, Table 1 shows the estimated relative errors \(Error_1\) and Error from the different values of the parameters of the problem \(n=m\), k(N) and \(\varepsilon\).
Reflection and conclusion.
In the introduction we have cited and talked a bit about various references in the literature that have dealt with the problem of holes with various methods and extended techniques to fill these holes, our idea here was to present a comparison with these works, but due to the originality of our work we believe that such a comparison cannot be made. In fact, no paper presents the offset problem, even our extended method for filling holes is different. The complexity of our work translates into mixing two approximation methods, the first is offset and the second is filling holes. Studying these criteria together is not so easy.
It should be noted that approach approximation method of this paper is made with the biquadratic splines functions that are of class 1, while for example in [2] it is done with the bi-cubic splines that are of class 2. This clearly shows that the degree of approximation in this work has been significantly reduced. Moreover, from Table 1, one can observe that as the data points increases the estimation of the error decreases, and also when the value of the parameter \(\varepsilon\) decreases the estimation of the error decreases, hence we can conclude the compatibility between the theory of the convergence result and the numerical ones. In short, we can conclude form the study of the table, where several data of the problem appear, and figures the effectiveness of our algorithm as an approximation method.
As an open research for the future is the possibility of extending the manuscript to find some numerical application directly or indirectly with chemistry.
References
S. Adamovic, M. Prica, B. Dalmacija, S. Rapajic, D. Novakovic, Z. Pavlovic, S. Maletic, Feasibility of electrocoagulation/flotation treatment of waste offset printing developer based on the response surface analysis. Arab. J. Chem. 9, 152–162 (2016)
R. Akhrif, A. Kouibia, M. Pasadas, Approximation of generalized offset surfaces by bicubic splines. J. Math. Chem. 58, 647–662 (2020)
E. Bréchner, General tool Offset Curves and Surfaces, in Geometry Processing for Design and Manufacturing (SAM, Philadelphia, 1992)
E. Busby, N.C. Anderson, J.S. Owen, M.Y. Sfeir, Effect of surface stoichiometry on blinking and hole trapping dynamics in CdSe nanocrystals. J. Phys. Chem. 119, 27797–27803 (2015)
P.G. Ciarlet, The Finite Element Method for Elliptic Problems (North-Amsterdam, 1978)
L. Damario, J. Fohlinger, G. Boschloo, L. Hammarstrom, Unveiling hole trapping and surface dynamics of NiO nanoparticles. Chem. Sci. 9, 223–230 (2018)
M.A. Fortes, P. González, M. Pasadas, M.L. Rodríguez, Hole fillig on surfaces by discrete variational splines. Appl. Numer. Math. 62(9), 1050–1060 (2012)
M.A. Fortes, P. González, M. Pasadas, M.L. Rodríguez, A hole filling method for explicit and parametric surfaces by using \(C^1\)-Powell-Sabin splines. Math. Comput. Simul. 99, 71–81 (2014)
M.A. Fortes, P. González, A. Palomares, M. Pasadas, Filling holes with shape preserving conditions. Math. Comput. Simul. 118, 198–212 (2015)
M.A. Fortes, P. González, A. Palomares, M. Pasadas, Filling holes with geometric and volumetric constraints. Comput. Math. Appl. 73(4), 671–683 (2017)
M.A. Fortes, P. González, A. Palomares, M.L. Rodríguez, Filling holes using a mesh of filled curves. Math. Comput. Simul. 164, 78–93 (2019)
Z. Gai, W. Lin, J.D. Burton, K. Fuchigami, P.C. Snijders, T.Z. Ward, E.Y. Tsymbal, J. Shen, S. Jesse, S.V. Kalinin, A.P. Baddorf, Chemically induced Jahn-Teller ordering on manganite surfaces. Nat. Commun. (2014). https://doi.org/10.1038/ncomms5528
J.D. Keene, J.R. McBride, N.J. Orfield, S.J. Rosenthal, Elimination of hole-surface overlap in graded CdS\(_x\)Se\(_{1-x}\) nanocrystals revealed by ultrafast fluorescence upconversion spectroscopy. Am. Chem. Soc. 8(10), 10665–10673 (2014)
A. Kouibia, M. Pasadas, Approximation by discrete variational splines. J. Comput. Appl. Math. 116, 145–156 (2000)
A. Kouibia, M. Pasadas, Approximation of surfaces by fairness bicubic splines. Adv. Comput. Math. 20, 87–103 (2004)
Y. Ma, F. Le Formal, A. Kafizas, S.R. Pendleburya, J.R. Durrant, Efficient suppression of back electron/hole recombination in cobalt phosphate surface modified undoped bismuth vanadate photoanodes. J. Math. Chem. 3, 20649–20657 (2015)
R. Odion, P. Strobbia, T. Vo-Dinh, Surface-enhanced patially offset Raman spectroscopy (SESORS) for biomedical applications, Proc. SPIE 10484, Advanced Biomedical and Clinical Diagnostic and Surgical Guidance Systems XVI, p. 1048404 (2018). https://doi.org/10.1117/12.2300088
H. Pottman, General offset surfaces. Neural parallel Sci. Comput. 5, 55–80 (1997)
F. Tao, Nanoscale surface chemistry in self-and directed-assembly of organic molecules on solid surfaces and synthesis of nanostructed organic architectures. Pure Appl. Chem. 80(1), 45–57 (2008)
C. Wang, P. Hu, A hole-filling of triangular meshes in engineering. Int. J. Comput. Methods Eng. Sci. Mech. 14(5), 465–571 (2013)
W. Zhu, Z. Huang, M. Bo, J. Huang, C. Peng, H. Liu, Surface, size and thermal effects in Alkali metal with core-electron binding-energy shifts. Chin. J. Chem. Phys. 34(5), 628–638 (2021)
Acknowledgements
This work has been supported by FEDER/Junta de Andalucía-Consejería de Transformaciín Económica, Industria, Conocimiento y Universidades (Research Project A-FQM-76-UGR20, University of Granada) and by the Junta de Andalucía (Research Group FQM191).
Funding
Funding for open access charge: Universidad de Granada / CBUA.
Author information
Authors and Affiliations
Corresponding author
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Kouibia, A., Pasadas, M. Reconstruction approximating method by biquadratic splines of offset surfaces holes. J Math Chem 60, 423–439 (2022). https://doi.org/10.1007/s10910-021-01322-7
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10910-021-01322-7