• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.28 no.1
• /
• pp.22-32
• /
• 2024

In this paper we present a method of conic section approximation by hexic Bézier curves. The hexic Bézier approximants are G³ Hermite interpolations of conic sections. We show that there exists at least one hexic Bézier approximant for each weight of the conic section The hexic Bézier approximant depends one parameter and it can be obtained by solving a quartic polynomial, which is solvable algebraically. We present the explicit upper bound of the Hausdorff distance between the conic section and the hexic Bézier approximant. We also prove that our approximation method has the maximal order of approximation. The numerical examples for conic section approximation by hexic Bézier curves are given and illustrate our assertions.

• Proceedings of the KOSOMBE Conference
• /
• v.1997 no.11
• /
• pp.459-463
• /
• 1997

In this paper, we present the elimination of ECG artifact from the polysomnographic EEG and EOG. The idea of this method is that the ECG synchronized EEG segment is detected from ECG and regard samples of that segment a missing signal. After this, we used two interpolation methods to recover the missing segment. One is the Lagrange Polynomial Interpolation Method and the other is the Least Square Error AR Interpolation method. We tested those methods by applying to simulated signals. AR methods works well enough to reject the artifact about 10% of the main artifact level. We practically applied to real EEG and EOG signals. We also developed the algorithm to detect whether the artifact level is high or not. If the artifact level is high, then the interpolations are applied.

• Economic and Environmental Geology
• /
• v.25 no.1
• /
• pp.73-77
• /
• 1992

This paper compares numerical quadrature techniques for computing an inverse Fourier transform integral in two-dimensional resistivity modeling. The quadrature techniques using exponential and cubic spline interpolations are examined for the case of a homogeneous earth model. In both methods the integral over the interval from 0 to ${\lambda}_{min}$, where ${\lambda}_{min}$, is the minimum sampling spatial wavenumber, is calculated by approximating Fourier transformed potentials to a logarithmic function. This scheme greatly reduces the inverse Fourier transform error associated with the logarithmic discontinuity at ${\lambda}=0$. Numrical results show that, if the sampling intervals are adequate, the cubic spline interpolation method is more accurate than the exponential interpolation method.

• Structural Engineering and Mechanics
• /
• v.38 no.3
• /
• pp.261-281
• /
• 2011

A new method of formulation of a class of elements that are immune to mesh distortion effects is proposed here. The simple three-noded bar element with an offset of the internal node from the element center is employed here to demonstrate the method and the principles on which it is founded upon. Using the function space approach, the modified formulation is shown here to be superior to the conventional isoparametric version of the element since it satisfies the completeness requirement as the metric formulation, and yet it is in agreement with the best-fit paradigm in both the metric and the parametric domains. Furthermore, the element error is limited to only those that are permissible by the classical projection theorem of strains and stresses. Unlike its conventional counterpart, the modified element is thus not prone to any errors from mesh distortion. The element formulation is symmetric and thus satisfies the requirement of the conservative nature of problems associated with all self-adjoint differential operators. The present paper indicates that a proper mapping set for distortion immune elements constitutes geometric and displacement interpolations through parametric and metric shape functions respectively, with the metric components in the displacement/strain replaced by the equivalent geometric interpolation in parametric co-ordinates.

• Journal of the Korean Society for Precision Engineering
• /
• v.16 no.9
• /
• pp.172-178
• /
• 1999

3D linear and circular interpolations are a basic part for the machining of complex shapes. Until now, because of the absence of appropriate algorithms for the generation of 3D lines and circles, a full accomplishment for available machine tool resolution is difficult. this paper presents new algorithms for 3D linear and circular interpolation in the reference pulse technique. In 3D space, the line or circle is not expressed as an implicit function, it is only defined as the intersection of two surfaces. A 3D line is defined as the intersection of two planes, and a 3D circle is defined as the intersection of a plane and the surface of a sphere. Based on these concepts, interpolation algorithms are designed to follow intersection curves in 3D space, and a real-time 3D linear and circular interpolator was developed in software using a PC. The algorithm implemented in a PC showed promising results in interpolation error and speed performance. It is expected that it can be applied to the next generation computerized numerical control systems for the machining of 3D lines, circles and some other complex shapes.

• Journal of the Korean Institute of Intelligent Systems
• /
• v.19 no.5
• /
• pp.683-688
• /
• 2009

The purposes of this paper is to propose new techniques to reconstruct measured optical images by using the template matching and B-Spline interpolation method based on image processing technology. To do this, we detect the matching template and non-matching template from each optical image. And then, we match the overlaped images from base level by correcting roll, pitch, and yaw error of images. At last, the matching image is interpolated by B-Spline interpolation. Finally, we show the effectiveness and feasibility of the proposed method through some experiments.

• Korean Journal of Agricultural and Forest Meteorology
• /
• v.2 no.4
• /
• pp.175-182
• /
• 2000

Spatial interpolation of daily temperature forecasts and observations issued by public weather services is frequently required to make them applicable to agricultural activities and modeling tasks. In contrast to the long term averages like monthly normals, terrain effects are not considered in most spatial interpolations for short term temperatures. This may cause erroneous results in mountainous regions where the observation network hardly covers full features of the complicated terrain. We developed a spatial interpolation model for daily minimum temperature which combines inverse distance squared weighting and elevation difference correction. This model uses a time dependent function for 'mountain slope lapse rate', which can be derived from regression analyses of the station observations with respect to the geographical and topographical features of the surroundings including the station elevation. We applied this model to interpolation of daily minimum temperature over the mountainous Korean Peninsula using 63 standard weather station data. For the first step, a primitive temperature surface was interpolated by inverse distance squared weighting of the 63 point data. Next, a virtual elevation surface was reconstructed by spatially interpolating the 63 station elevation data and subtracted from the elevation surface of a digital elevation model with 1 km grid spacing to obtain the elevation difference at each grid cell. Final estimates of daily minimum temperature at all the grid cells were obtained by applying the calculated daily lapse rate to the elevation difference and adjusting the inverse distance weighted estimates. Independent, measured data sets from 267 automated weather station locations were used to calculate the estimation errors on 12 dates, randomly selected one for each month in 1999. Analysis of 3 terms of estimation errors (mean error, mean absolute error, and root mean squared error) indicates a substantial improvement over the inverse distance squared weighting.

• Journal of the Korea Society for Simulation
• /
• v.16 no.3
• /
• pp.57-63
• /
• 2007

In real time rendering of fetus the empty space leaping while traversing a ray is most frequently used accelerating technique. The main idea is to skip empty voxel samples which do not contribute the result image and it speeds up the rendering time by avoiding sampling data while traversing a ray in the empty region, saving a substantial number of interpolations. Calculating the distance from the nearest object boundary for every yokel can reduce the sampling operation. Among widely-well-known distance maps, those estimates the true distance, such as euclidean distance, takes a long time to compute because of the complicated floating-point operations, and others which uses approximated distance functions, such as city-block and chessboard, provides faster computation time but sampling error may can occur. In this paper, therefore, we analyze the characteristics of several distance maps and compare the number of samples and rendering time. And we aim to suggest the most appropriate distance map for rendering of fetus in ultrasound image.

• Korean Journal of Agricultural and Forest Meteorology
• /
• v.4 no.2
• /
• pp.95-102
• /
• 2002

Spatial interpolation has become a common procedure in converting temperature forecasts and observations at irregular points for use in regional scale ecosystem modeling and the model based decision support systems for resource management. Neglection of terrain effects in most spatial interpolations for short term temperatures may cause erroneous results in mountainous regions, where the observation network hardly covers full features of the complicated terrain. A spatial interpolation model for daytime hourly temperature was formulated based on error analysis of unsampled site with respect to the site topography. The model has a solar irradiance correction scheme in addition to the common backbone of the lapse rate - corrected inverse distance weighting. The solar irradiance scheme calculates the direct, diffuse and reflected components of shortwave radiation over any surfaces based on the sun-slope geometry and compares the sum with that over a reference surface. The deviation from the reference radiation is used to calculate the temperature correction term by an empirical conversion formula between the solar energy and the air temperature on any sloped surfaces at an hourly time scale, which can be prepared seasonally for each land cover type. When this model was applied to a 14 km by 22 km mountainous region at a 10 m horizontal resolution, the estimated hourly temperature surfaces showed a better agreement with the observed distribution than those by a conventional method.

• The Korean Journal of Nuclear Medicine
• /
• v.39 no.6
• /
• pp.464-472
• /
• 2005

Purpose: SPECT using a fan-beam collimator improves spatial resolution and sensitivity. For the reconstruction from fan-beam projections, it is necessary to implement direct fan-beam reconstruction methods without transforming the data into the parallel geometry. In this study, various fan-beam reconstruction algorithms were implemented and their performances were compared. Materials and Methods: The projector for fan-beam SPECT was implemented using a ray-tracing method. The direct reconstruction algorithms implemented for fan-beam projection data were FBP (filtered backprojection), EM (expectation maximization), OS-EM (ordered subsets EM) and MAP-EM OSL (maximum a posteriori EM using the one-step late method) with membrane and thin-plate models as priors. For comparison, the fan-beam protection data were also rebinned into the parallel data using various interpolation methods, such as the nearest neighbor, bilinear and bicubic interpolations, and reconstructed using the conventional EM algorithm for parallel data. Noiseless and noisy projection data from the digital Hoffman brain and Shepp/Logan phantoms were reconstructed using the above algorithms. The reconstructed images were compared in terms of a percent error metric. Results: for the fan-beam data with Poisson noise, the MAP-EM OSL algorithm with the thin-plate prior showed the best result in both percent error and stability. Bilinear interpolation was the most effective method for rebinning from the fan-beam to parallel geometry when the accuracy and computation load were considered. Direct fan-beam EM reconstructions were more accurate than the standard EM reconstructions obtained from rebinned parallel data. Conclusion: Direct fan-beam reconstruction algorithms were implemented, which provided significantly improved reconstructions.