Scattered Data Interpolation Using Quartic Triangular Patch for Shape-Preserving Interpolation and Comparison with Mesh-Free Methods

Scattered data interpolation is important in sciences, engineering, and medical-based problems. Quartic Bézier triangular patches with 15 control points (ordinates) can also be used for scattered data interpolation. However, this method has a weakness; that is, in order to achieve <inline-formula...

Full description

Bibliographic Details
Main Authors: Samsul Ariffin Abdul Karim, Azizan Saaban, Van Thien Nguyen
Format: Article
Language:English
Published: MDPI AG 2020-06-01
Series:Symmetry
Subjects:
Online Access:https://www.mdpi.com/2073-8994/12/7/1071
_version_ 1797563747402776576
author Samsul Ariffin Abdul Karim
Azizan Saaban
Van Thien Nguyen
author_facet Samsul Ariffin Abdul Karim
Azizan Saaban
Van Thien Nguyen
author_sort Samsul Ariffin Abdul Karim
collection DOAJ
description Scattered data interpolation is important in sciences, engineering, and medical-based problems. Quartic Bézier triangular patches with 15 control points (ordinates) can also be used for scattered data interpolation. However, this method has a weakness; that is, in order to achieve <inline-formula> <math display="inline"> <semantics> <mrow> <msup> <mi>C</mi> <mn>1</mn> </msup> </mrow> </semantics> </math> </inline-formula> continuity, the three inner points can only be determined using an optimization method. Thus, we cannot obtain the exact Bézier ordinates, and the quartic scheme is global and not local. Therefore, the quartic Bézier triangular has received less attention. In this work, we use Zhu and Han’s quartic spline with ten control points (ordinates). Since there are only ten control points (as for cubic Bézier triangular cases), all control points can be determined exactly, and the optimization problem can be avoided. This will improve the presentation of the surface, and the process to construct the scattered surface is local. We also apply the proposed scheme for the purpose of positivity-preserving scattered data interpolation. The sufficient conditions for the positivity of the quartic triangular patches are derived on seven ordinates. We obtain nonlinear equations that can be solved using the regula-falsi method. To produce the interpolated surface for scattered data, we employ four stages of an algorithm: (a) triangulate the scattered data using Delaunay triangulation; (b) assign the first derivative at the respective data; (c) form a triangular surface via convex combination from three local schemes with <inline-formula> <math display="inline"> <semantics> <mrow> <msup> <mi>C</mi> <mn>1</mn> </msup> </mrow> </semantics> </math> </inline-formula> continuity along all adjacent triangles; and (d) construct the scattered data surface using the proposed quartic spline. Numerical results, including some comparisons with some existing mesh-free schemes, are presented in detail. Overall, the proposed quartic triangular spline scheme gives good results in terms of a higher coefficient of determination (R<sup>2</sup>) and smaller maximum error (Max Error), requires about 12.5% of the CPU time of the quartic Bézier triangular, and is on par with Shepard triangular-based schemes. Therefore, the proposed scheme is significant for use in visualizing large and irregular scattered data sets. Finally, we tested the proposed positivity-preserving interpolation scheme to visualize coronavirus disease 2019 (COVID-19) cases in Malaysia.
first_indexed 2024-03-10T18:47:46Z
format Article
id doaj.art-bcf2a2c13fce46fc8706b54133f9a54d
institution Directory Open Access Journal
issn 2073-8994
language English
last_indexed 2024-03-10T18:47:46Z
publishDate 2020-06-01
publisher MDPI AG
record_format Article
series Symmetry
spelling doaj.art-bcf2a2c13fce46fc8706b54133f9a54d2023-11-20T05:24:49ZengMDPI AGSymmetry2073-89942020-06-01127107110.3390/sym12071071Scattered Data Interpolation Using Quartic Triangular Patch for Shape-Preserving Interpolation and Comparison with Mesh-Free MethodsSamsul Ariffin Abdul Karim0Azizan Saaban1Van Thien Nguyen2Fundamental and Applied Sciences Department and Centre for Smart Grid Energy Research (CSMER), Institute of Autonomous System, Universiti Teknologi PETRONAS, Bandar Seri Iskandar, Seri Iskandar 32610, Perak Darul Ridzuan, MalaysiaSchool of Quantitative Sciences, UUMCAS, Universiti Utara Malaysia, Kedah 06010, MalaysiaFPT University, Education Zone, Hoa Lac High Tech Park, Km29 Thang Long Highway, Thach That Ward, Hanoi 10000, VietnamScattered data interpolation is important in sciences, engineering, and medical-based problems. Quartic Bézier triangular patches with 15 control points (ordinates) can also be used for scattered data interpolation. However, this method has a weakness; that is, in order to achieve <inline-formula> <math display="inline"> <semantics> <mrow> <msup> <mi>C</mi> <mn>1</mn> </msup> </mrow> </semantics> </math> </inline-formula> continuity, the three inner points can only be determined using an optimization method. Thus, we cannot obtain the exact Bézier ordinates, and the quartic scheme is global and not local. Therefore, the quartic Bézier triangular has received less attention. In this work, we use Zhu and Han’s quartic spline with ten control points (ordinates). Since there are only ten control points (as for cubic Bézier triangular cases), all control points can be determined exactly, and the optimization problem can be avoided. This will improve the presentation of the surface, and the process to construct the scattered surface is local. We also apply the proposed scheme for the purpose of positivity-preserving scattered data interpolation. The sufficient conditions for the positivity of the quartic triangular patches are derived on seven ordinates. We obtain nonlinear equations that can be solved using the regula-falsi method. To produce the interpolated surface for scattered data, we employ four stages of an algorithm: (a) triangulate the scattered data using Delaunay triangulation; (b) assign the first derivative at the respective data; (c) form a triangular surface via convex combination from three local schemes with <inline-formula> <math display="inline"> <semantics> <mrow> <msup> <mi>C</mi> <mn>1</mn> </msup> </mrow> </semantics> </math> </inline-formula> continuity along all adjacent triangles; and (d) construct the scattered data surface using the proposed quartic spline. Numerical results, including some comparisons with some existing mesh-free schemes, are presented in detail. Overall, the proposed quartic triangular spline scheme gives good results in terms of a higher coefficient of determination (R<sup>2</sup>) and smaller maximum error (Max Error), requires about 12.5% of the CPU time of the quartic Bézier triangular, and is on par with Shepard triangular-based schemes. Therefore, the proposed scheme is significant for use in visualizing large and irregular scattered data sets. Finally, we tested the proposed positivity-preserving interpolation scheme to visualize coronavirus disease 2019 (COVID-19) cases in Malaysia.https://www.mdpi.com/2073-8994/12/7/1071quartic splinetriangulationscattered datacontinuitysurface reconstructionpositivity-preserving
spellingShingle Samsul Ariffin Abdul Karim
Azizan Saaban
Van Thien Nguyen
Scattered Data Interpolation Using Quartic Triangular Patch for Shape-Preserving Interpolation and Comparison with Mesh-Free Methods
Symmetry
quartic spline
triangulation
scattered data
continuity
surface reconstruction
positivity-preserving
title Scattered Data Interpolation Using Quartic Triangular Patch for Shape-Preserving Interpolation and Comparison with Mesh-Free Methods
title_full Scattered Data Interpolation Using Quartic Triangular Patch for Shape-Preserving Interpolation and Comparison with Mesh-Free Methods
title_fullStr Scattered Data Interpolation Using Quartic Triangular Patch for Shape-Preserving Interpolation and Comparison with Mesh-Free Methods
title_full_unstemmed Scattered Data Interpolation Using Quartic Triangular Patch for Shape-Preserving Interpolation and Comparison with Mesh-Free Methods
title_short Scattered Data Interpolation Using Quartic Triangular Patch for Shape-Preserving Interpolation and Comparison with Mesh-Free Methods
title_sort scattered data interpolation using quartic triangular patch for shape preserving interpolation and comparison with mesh free methods
topic quartic spline
triangulation
scattered data
continuity
surface reconstruction
positivity-preserving
url https://www.mdpi.com/2073-8994/12/7/1071
work_keys_str_mv AT samsulariffinabdulkarim scattereddatainterpolationusingquartictriangularpatchforshapepreservinginterpolationandcomparisonwithmeshfreemethods
AT azizansaaban scattereddatainterpolationusingquartictriangularpatchforshapepreservinginterpolationandcomparisonwithmeshfreemethods
AT vanthiennguyen scattereddatainterpolationusingquartictriangularpatchforshapepreservinginterpolationandcomparisonwithmeshfreemethods