Showing 1–33 of 33 results for author: Skala, V

Line-Torus Intersection for Ray Tracing: Alternative Formulations

Authors: Vaclav Skala

Abstract: Intersection algorithms are very important in computation of geometrical problems. Algorithms for a line intersection with linear or quadratic surfaces are quite efficient. However, algorithms for a line intersection with other surfaces are more complex and time consuming. In this case the object is usually closed into a simple bounding volume to speed up the cases when the given line cannot inter… ▽ More Intersection algorithms are very important in computation of geometrical problems. Algorithms for a line intersection with linear or quadratic surfaces are quite efficient. However, algorithms for a line intersection with other surfaces are more complex and time consuming. In this case the object is usually closed into a simple bounding volume to speed up the cases when the given line cannot intersect the given object. In this paper new formulations of the line-torus intersection problem are given and new specification of the bounding volume for a torus is given as well. The presented approach is based on an idea of a line intersection with an envelope of rotating sphere that forms a torus. Due to this approach new bounding volume can be formulated which is more effective as it enables to detect cases when the line passes the "hole" of a torus, too. △ Less

Submitted 9 January, 2023; originally announced January 2023.

Comments: Draft of the paper published Line-Torus Intersection: Alternative Formulations, WSEAS Trans. on Computers, ISSN 2224-2872, Vol7., No.12, pp.288-297, 2013

ACM Class: I.3.5

Journal ref: WSEAS Trans. on Computers, ISSN 2224-2872, Vol7., No.12, pp.288-297, 2013

A Comparison of Fundamental Methods for Iso-surface Extraction

Authors: Jan Patera, Vaclav Skala

Abstract: In this paper four fundamental methods for an iso-surface extraction are compared, based on cell decomposition to tetrahedra. The methods are compared both on mathematically generated data sets as well as on real data sets. The comparison using mathematical data is made from different points of view such as area approximation, volume approximation. On the other hand, the Hausdorff distance and roo… ▽ More In this paper four fundamental methods for an iso-surface extraction are compared, based on cell decomposition to tetrahedra. The methods are compared both on mathematically generated data sets as well as on real data sets. The comparison using mathematical data is made from different points of view such as area approximation, volume approximation. On the other hand, the Hausdorff distance and root mean square are used to compare methods on real data sets. The presented comparison can be helpful when deciding among tested methods which one to choose, as well as when we need to compare a newly developed method with other existing approaches. △ Less

Submitted 4 January, 2023; originally announced January 2023.

Comments: Draft of the paper submitted to Machine Graphics and Vision, Polish Academy of Sciences, Vol.13, No.4., pp.329-344, ISSN 1230-0535, 2004

Journal ref: Machine Graphics and Vision, Polish Academy of Sciences, Vol.13, No.4., pp.329-344, ISSN 1230-0535, 2004

Robust Surface Reconstruction from Orthogonal Slices

Authors: Radek Svitak, Vaclav Skala

Abstract: The surface reconstruction problem from sets of planar parallel slices representing cross sections through 3D objects is presented. The final result of surface reconstruction is always based on the correct estimation of the structure of the original object. This paper is a case study of the problem of the structure determination. We present a new approach, which is based on considering mutually or… ▽ More The surface reconstruction problem from sets of planar parallel slices representing cross sections through 3D objects is presented. The final result of surface reconstruction is always based on the correct estimation of the structure of the original object. This paper is a case study of the problem of the structure determination. We present a new approach, which is based on considering mutually orthogonal sets of slices. A new method for surface reconstruction from orthogonal slices is described and the benefit of orthogonal slices is discussed too. △ Less

Submitted 4 January, 2023; originally announced January 2023.

Comments: Draft of the papers published in Machine Graphics and Vision, Vol.13, No.3, Polish Academy of Sciences, Vol.13, No.3, pp.221-233, ISSN 1230-0535, 2004

Journal ref: Machine Graphics and Vision, Vol.13, No.3, Polish Academy of Sciences, Vol.13, No.3, pp.221-233, ISSN 1230-0535, 2004

Poseidon: Non-server WEB Forms Off-line Processing System

Authors: Vaclav Skala

Abstract: The proposed Poseidon system is based on email services of filled forms instead of WEB server based services. This approach is convenient especially for small applications or small-medium companies. It is based on PDF forms that are available on a WEB page. PDF forms can be downloaded, off-line filled in, printed and finally sent by email for final processing. Data are actually stored in the local… ▽ More The proposed Poseidon system is based on email services of filled forms instead of WEB server based services. This approach is convenient especially for small applications or small-medium companies. It is based on PDF forms that are available on a WEB page. PDF forms can be downloaded, off-line filled in, printed and finally sent by email for final processing. Data are actually stored in the local outbox waiting for a connection to a mail server. This follows an idea of the standard "paper" letter sending. Filled in data are actually sent when a user is on-line, therefore a user is "free" of being on-line when filling the forms. When the PDF form is processed on the recipient side, answer is sent back via email as well. Typical application is e.g. in conference management systems, systems for submission to journals etc. The great advantage of the PDF forms use is that they can be easily made or modified by a non-specialized administrative person easily. △ Less

Submitted 27 December, 2022; originally announced December 2022.

Comments: Draft of the paper submitted to International Journal of Computers, ISSN: 2367-8895

ACM Class: H.4

HS-Patch: A New Hermite Smart Bicubic Patch Modification

Authors: Vaclav Skala, Michal Smolik, Lukas Karlicek

Abstract: Bicubic four-sided patches are widely used in computer graphics, CAD/CAM systems etc. Their flexibility is high and enables to compress a surface description before final rendering. However, computer graphics hardware supports only triangular meshes. Therefore, four-sided bicubic patches are approximated by a triangular mesh. The border curves of a bicubic patch are of degree 3, while diagonal and… ▽ More Bicubic four-sided patches are widely used in computer graphics, CAD/CAM systems etc. Their flexibility is high and enables to compress a surface description before final rendering. However, computer graphics hardware supports only triangular meshes. Therefore, four-sided bicubic patches are approximated by a triangular mesh. The border curves of a bicubic patch are of degree 3, while diagonal and anti-diagonal curves are of degree 6. Therefore the resulting shape and texturing depend on the actual mapping, i.e. how the tessellation of a bicubic patch is made. The proposed new modification of the Hermite bicubic patch, the HS-patch, is a result of additional restriction put on the Hermite bicubic patch formulation - the diagonal and anti-diagonal curves are of degree 3. This requirement leads to a new Hermite based bicubic four-sided patch with 12 control points and another 4 control points, i.e. twist vectors, are computed from those 12 control points. △ Less

Submitted 22 December, 2022; originally announced December 2022.

Comments: Draft of the paper: NAUN Journal International Journal of Mathematics and Computers in Simulation, Vol.8, pp.292-299, ISSN: 1998-0159, 2014. arXiv admin note: substantial text overlap with arXiv:2212.11986, arXiv:2212.11875

ACM Class: I.3

Journal ref: NAUN Journal International Journal of Mathematics and Computers in Simulation, Vol.8, pp.292-299, ISSN: 1998-0159, 2014

BS-Patch: Constrained Bezier Parametric Patch

Authors: Vaclav Skala, Vit Ondracka

Abstract: Bezier parametric patches are used in engineering practice quite often, especially in CAD/CAM systems oriented to mechanical design. In many cases quadrilateral meshes are used for tessellation of parameters domain. We propose a new modification of the Bezier cubic rectangular patch, the BS-patch, which is based on the requirement that diagonal curves must be of degree 3 instead of degree 6 as it… ▽ More Bezier parametric patches are used in engineering practice quite often, especially in CAD/CAM systems oriented to mechanical design. In many cases quadrilateral meshes are used for tessellation of parameters domain. We propose a new modification of the Bezier cubic rectangular patch, the BS-patch, which is based on the requirement that diagonal curves must be of degree 3 instead of degree 6 as it is in the case of the Bezier patch. Theoretical derivation of conditions is presented and some experimental results as well. The BS-Patch is convenient for applications where for different tessellation of the u-v domain different degrees of diagonal curves are not acceptable. △ Less

Submitted 22 December, 2022; originally announced December 2022.

Comments: Draft of the paper: WSEAS Trans.on Mathematics, E-ISSN 2224-2880, No.5, Vol.12, pp.598-607, 2013. arXiv admin note: text overlap with arXiv:2212.11875

ACM Class: I.3

Journal ref: WSEAS Trans.on Mathematics, E-ISSN 2224-2880, No.5, Vol.12, pp.598-607, 2013

S-patch: Modification of the Hermite parametric patch

Authors: Vaclav Skala, Vit Ondracka

Abstract: A new modification of the Hermite cubic rectangular patch is proposed: the S-Patch, which is based on the requirement that diagonal curves must be of degree 3 instead of degree 6 as it is in the case of the Hermite patch. Theoretical derivation of conditions is presented and some experimental results as well. The S-Patch is convenient for applications, where different tessellation of the u-v domai… ▽ More A new modification of the Hermite cubic rectangular patch is proposed: the S-Patch, which is based on the requirement that diagonal curves must be of degree 3 instead of degree 6 as it is in the case of the Hermite patch. Theoretical derivation of conditions is presented and some experimental results as well. The S-Patch is convenient for applications, where different tessellation of the u-v domain is needed, boundary and diagonal curves of different degrees are not acceptable. △ Less

Submitted 22 December, 2022; originally announced December 2022.

Comments: Draft of the paper: S-Patch: Modification of the Hermite Parametric Patch, ICGG 2010 conference, Kyoto, Japan, 2010

ACM Class: I.3

Summation Problem Revisited -- More Robust Computation

Authors: Vaclav Skala

Abstract: Numerical data processing is a key task across different fields of computer technology use. However, even simple summation of values is not precise due to the floating point representation use. This paper presents a practical algorithm for summation of values convenient for medium and large data sets. The proposed algorithm is simple, easy to implement. Its computational complexity is O(N) in the… ▽ More Numerical data processing is a key task across different fields of computer technology use. However, even simple summation of values is not precise due to the floating point representation use. This paper presents a practical algorithm for summation of values convenient for medium and large data sets. The proposed algorithm is simple, easy to implement. Its computational complexity is O(N) in the contrary of the Exact Sign Summation Algorithm (ESSA) approach with O(N^2) run-time complexity. The proposed algorithm is especially convenient for cases when exponent data differ significantly and many small values are summed with higher values △ Less

Submitted 8 November, 2022; originally announced November 2022.

Comments: 9 pages, 3 Figs, 3 Tabs. Presented at Recent Advances in Computer Science Conf, 2013

ACM Class: F.2.1

Algorithm for Finding an Exact Maximum Distance in E2 with Oexp(N) Complexity: Analysis and Experimental Results

Authors: Vaclav Skala

Abstract: This paper describes a novel and fast, simple and robust algorithm with O(N) expected complexity which enables to decrease run time needed to find the maximum distance of two points in E2. It can be easily modified for the E3 case in general. The proposed algorithm has been evaluated experimentally on larger different datasets in order to verify it and prove expected properties of it. Experiments… ▽ More This paper describes a novel and fast, simple and robust algorithm with O(N) expected complexity which enables to decrease run time needed to find the maximum distance of two points in E2. It can be easily modified for the E3 case in general. The proposed algorithm has been evaluated experimentally on larger different datasets in order to verify it and prove expected properties of it. Experiments proved the advantages of the proposed algorithm over the standard algorithms based on the Brute force, convex hull or convex hull diameters approaches. The proposed algorithm gives a significant speed-up to applications, when medium and large data sets are processed. It is over 10 000 times faster than the standard Brute force algorithm for 10 mil. points randomly distributed points in E2 and over 4 times faster than convex hull diameter computation. The speed-up of the proposed algorithm grows with the number of points processed. △ Less

Submitted 4 September, 2022; originally announced September 2022.

Comments: This contribution was partially presented as an extended abstract at the ICNAAM 2013 conference, DOI: 10.1063/1.4826047, arXiv DOI:10.48550/arXiv.2208.04730

MSC Class: 02.60.-x; 02.30.Jr; 02.60 Dc; 89.20.Ff

Polar, Spherical and Orthogonal Space Subdivisions for an Algorithm Acceleration: O(1) Point-in-Polygon/Polyhedron Test

Authors: Vaclav Skala

Abstract: Acceleration of algorithms is becoming a crucial problem, if larger data sets are to be processed. Evaluation of algorithms is mostly done by using computational geometry approach and evaluation of computational complexity. However in todays engineering problems this approach does not respect that number of processed items is always limited and a significant role plays also speed of read/write ope… ▽ More Acceleration of algorithms is becoming a crucial problem, if larger data sets are to be processed. Evaluation of algorithms is mostly done by using computational geometry approach and evaluation of computational complexity. However in todays engineering problems this approach does not respect that number of processed items is always limited and a significant role plays also speed of read/write operations. One general method how to speed up an algorithm is application of space subdivision technique and usually the orthogonal space subdivision is used. In this paper non-orthogonal subdivisions are described. The proposed approach can significantly improve memory consumption and run-time complexity. The proposed modified space subdivision techniques are demonstrated on two simple problems Point-in-Convex Polygon and Point-in-Convex Polyhedron tests. △ Less

Submitted 26 August, 2022; originally announced August 2022.

Comments: 4 pages, 8 figures

MSC Class: 65D17; 65D18 ACM Class: I.3.5

Fast O_{expected}(N) Algorithm for Finding Exact Maximum Distance in E^2 Instead of O(N^2) or O(N lg N)

Authors: Vaclav Skala

Abstract: This paper describes novel and fast, simple and robust algorithm with O(N) expected complexity which enables to decrease run-time needed to find an exact maximum distance of two points in E2. The proposed algorithm has been evaluated experimentally on larger different datasets. The proposed algorithm gives a significant speed-up to applications, when medium and large data sets are processed. It is… ▽ More This paper describes novel and fast, simple and robust algorithm with O(N) expected complexity which enables to decrease run-time needed to find an exact maximum distance of two points in E2. The proposed algorithm has been evaluated experimentally on larger different datasets. The proposed algorithm gives a significant speed-up to applications, when medium and large data sets are processed. It is over 10 000 times faster than the standard algorithm for 10^6 points randomly distributed points in E2. Experiments proved the advantages of the proposed algorithm over the standard algorithm and convex hull diameters approaches. △ Less

Submitted 7 August, 2022; originally announced August 2022.

MSC Class: 68xx; 68U05

New Geometric Continuity Solution of Parametric Surfaces

Authors: Vaclav Skala

Abstract: This paper presents a new approach to computation of geometric continuity for parametric bi-cubic patches, based on a simple mathematical reformulation which leads to simple additional conditions to be applied in the patching computation. The paper presents an Hermite formulation of a bicubic parametric patch, but reformulations can be made also for Bézier and B-Spline patches as well. The present… ▽ More This paper presents a new approach to computation of geometric continuity for parametric bi-cubic patches, based on a simple mathematical reformulation which leads to simple additional conditions to be applied in the patching computation. The paper presents an Hermite formulation of a bicubic parametric patch, but reformulations can be made also for Bézier and B-Spline patches as well. The presented approach is convenient for the cases when valencies of corners are different from the value 4, in general. △ Less

Submitted 7 August, 2022; originally announced August 2022.

MSC Class: 68xx; 68U05

A New Approach to Line-Sphere and Line-Quadrics Intersection Detection and Computation

Authors: Vaclav Skala

Abstract: Line intersection with convex and un-convex polygons or polyhedron algorithms are well known as line clipping algorithms and very often used in computer graphics. Rendering of geometrical problems often leads to ray tracing techniques, when an intersection of many lines with spheres or quadrics is a critical issue due to ray-tracing algorithm complexity. A new formulation of detection and computat… ▽ More Line intersection with convex and un-convex polygons or polyhedron algorithms are well known as line clipping algorithms and very often used in computer graphics. Rendering of geometrical problems often leads to ray tracing techniques, when an intersection of many lines with spheres or quadrics is a critical issue due to ray-tracing algorithm complexity. A new formulation of detection and computation of the intersection of line (ray) with a quadric surface is presented, which separates geometric properties of the line and quadrics that enables pre-computation. The presented approach is especially convenient for implementation with SSE instructions or on GPU △ Less

Submitted 7 August, 2022; originally announced August 2022.

MSC Class: 68xx; 68U05

Point-in-Convex Polygon and Point-in-Convex Polyhedron Algorithms with O(1) Complexity using Space Subdivision

Authors: Vaclav Skala

Abstract: There are many space subdivision and space partitioning techniques used in many algorithms to speed up computations. They mostly rely on orthogonal space subdivision, resp. using hierarchical data structures, e.g. BSP trees, quadtrees, octrees, kd-trees, bounding volume hierarchies, etc. However in some applications a non-orthogonal space subdivision can offer new ways for actual speed up. In the… ▽ More There are many space subdivision and space partitioning techniques used in many algorithms to speed up computations. They mostly rely on orthogonal space subdivision, resp. using hierarchical data structures, e.g. BSP trees, quadtrees, octrees, kd-trees, bounding volume hierarchies, etc. However in some applications a non-orthogonal space subdivision can offer new ways for actual speed up. In the case of convex polygon in E3 a simple Point-in-Polygon test is of the O(N) complexity and the optimal algorithm is of O(lg N) computational complexity. In the E3 case, the complexity is O(N) even for the convex polyhedron as no ordering is defined. New Point-in-Convex Polygon and Point-in-Convex Polyhedron algorithms are presented based on space subdivision in the preprocessing stage resulting to O(1) run-time complexity. The presented approach is simple to implement. Due to the principle of duality, dual problems, e.g. line-convex polygon, line clipping, can be solved in a similarly. △ Less

Submitted 7 August, 2022; originally announced August 2022.

MSC Class: 68xx; 68U05

A New Robust Algorithm for Computation of a Triangle Circumscribed Sphere in E3 and a Hypersphere Simplex

Authors: Vaclav Skala

Abstract: There are many applications in which a bounding sphere containing the given triangle E3 is needed, e.g. fast collision detection, ray-triangle intersecting in raytracing etc. This is a typical geometrical problem in E3 and it has also applications in computational problems in general. In this paper a new fast and robust algorithm of circumscribed sphere computation in the -dimensional space is pre… ▽ More There are many applications in which a bounding sphere containing the given triangle E3 is needed, e.g. fast collision detection, ray-triangle intersecting in raytracing etc. This is a typical geometrical problem in E3 and it has also applications in computational problems in general. In this paper a new fast and robust algorithm of circumscribed sphere computation in the -dimensional space is presented and specification for the E3 space is given, too. The presented method is convenient for use on GPU or with SSE or Intel AVX instructions on a standard CPU △ Less

Submitted 7 August, 2022; originally announced August 2022.

MSC Class: 68xx; 68U05

Projective Geometry, Duality and Plucker Coordinates for Geometric Computations with Determinants on GPUs

Authors: Vaclav Skala

Abstract: Many algorithms used are based on geometrical computation. There are several criteria in selecting appropriate algorithm from already known. Recently, the fastest algorithms have been preferred. Nowadays, algorithms with a high stability are preferred. Also technology and computer architecture, like GPU etc., plays a significant role for large data processing. However, some algorithms are ill-cond… ▽ More Many algorithms used are based on geometrical computation. There are several criteria in selecting appropriate algorithm from already known. Recently, the fastest algorithms have been preferred. Nowadays, algorithms with a high stability are preferred. Also technology and computer architecture, like GPU etc., plays a significant role for large data processing. However, some algorithms are ill-conditioned due to numerical representation used; result of the floating point representation. In this paper, relations between projective representation, duality and Plucker coordinates will be explored with demonstration on simple geometric examples. The presented approach is convenient especially for application on GPUs or vector-vector computational architectures △ Less

Submitted 7 August, 2022; originally announced August 2022.

Comments: arXiv admin note: substantial text overlap with arXiv:1708.06684

MSC Class: 68xx; 68U05

Clipping and Intersection Algorithms: Short Survey and References

Authors: Vaclav Skala

Abstract: This contribution presents a brief survey of clipping and intersection algorithms in E2 and E3 with a nearly complete list of relevant references. Some algorithms use the projective extension of the Euclidean space and vector-vector operations, which supports GPU and SSE use. This contribution presents a brief survey of clipping and intersection algorithms in E2 and E3 with a nearly complete list of relevant references. Some algorithms use the projective extension of the Euclidean space and vector-vector operations, which supports GPU and SSE use. △ Less

Submitted 27 June, 2022; originally announced June 2022.

Comments: 22 pages

MSC Class: 68xx; 68U05

Two Methods for Iso-Surface Extraction from Volumetric Data and Their Comparison

Authors: Vaclav Skala, Alex Brusi

Abstract: There are various methods for extracting iso-surfaces from volumetric data. Marching cubes or tetrahedra or raytracing methods are mostly used. There are many specific techniques to increase speed of computation and decrease memory requirements. Although a precision of iso-surface extraction is very important, too, it is not mentioned usually. A comparison of the selected methods was made in diffe… ▽ More There are various methods for extracting iso-surfaces from volumetric data. Marching cubes or tetrahedra or raytracing methods are mostly used. There are many specific techniques to increase speed of computation and decrease memory requirements. Although a precision of iso-surface extraction is very important, too, it is not mentioned usually. A comparison of the selected methods was made in different aspects: iso-surface extraction process time, number of triangles generated and estimation of radius, area and volume errors based on approximation of a sphere. Surprisingly, experiments proved that there is no direct relation between precision of extracted and human perception of the extracted iso-surface △ Less

Submitted 3 January, 2022; originally announced January 2022.

MSC Class: 68U05 ACM Class: I.3

Journal ref: Machine Graphics & Vision, No.1/2, Vol.9, pp.149-166, Poland Academy of Sciences, Poland, ISSN 1230-0535, 2000

Line Clipping in E3 with Expected Complexity O(1)

Authors: Vaclav Skala

Abstract: A new line clipping algorithm against convex polyhedron in E3 with an expected complexity O(1) is presented. The suggested approach is based on two orthogonal projections to E2 co-ordinate system and on pre-processing of the given polyhedron. The pre-processing enables to speed up solution significantly. The proposed method is convenient for those applications when many lines are clipped against c… ▽ More A new line clipping algorithm against convex polyhedron in E3 with an expected complexity O(1) is presented. The suggested approach is based on two orthogonal projections to E2 co-ordinate system and on pre-processing of the given polyhedron. The pre-processing enables to speed up solution significantly. The proposed method is convenient for those applications when many lines are clipped against constant convex polyhedron. Theoretical considerations and experimental results are also presented. △ Less

Submitted 3 January, 2022; originally announced January 2022.

Comments: arXiv admin note: text overlap with arXiv:2111.07987

MSC Class: 68U05 ACM Class: I.3

Journal ref: Machine Graphics and Vision, Poland Academy of Sciences, Vol.5, No.4, pp.551-562, ISSN 1230-0535, 1996

Two New Algorithms for Line Clipping in E2 and Their Comparison

Authors: Vaclav Skala, Duc Huy Bui

Abstract: Many algorithms for clipping a line by a rectangular area or a convex polygon in E2 or by a non-convex or convex polyhedron in E3 have been published. The line segment clipping by the rectangular window in E2 is often restricted to the use of the Cohen-Sutherland (CS) algorithm or its modifications based on some presumptions like small clipping window or more sophisticated coding technique, etc. T… ▽ More Many algorithms for clipping a line by a rectangular area or a convex polygon in E2 or by a non-convex or convex polyhedron in E3 have been published. The line segment clipping by the rectangular window in E2 is often restricted to the use of the Cohen-Sutherland (CS) algorithm or its modifications based on some presumptions like small clipping window or more sophisticated coding technique, etc. The line clipping problem solution is a bottleneck of many packages and applications and, therefore, it would be desirable to use the fastest algorithm even though it is more complex. △ Less

Submitted 3 January, 2022; originally announced January 2022.

MSC Class: 68U05 ACM Class: I.3

Journal ref: Machine Graphics & Vision, Vol. 9, no. 1/2, pp. 297-306, 2000

A New Algorithm for Pyramidal Clipping of Line Segments in E3

Authors: Vaclav Skala, Duc Huy Bui

Abstract: A new algorithm for clipping a line segment against a pyramid in E3 is presented. This algorithm avoids computation of intersection points which are not end-points of the output line segment. It also allows solving all cases more effectively. The performance of this algorithm is shown to be consistently better than existing algorithms, including the Cohen-Sutherland, Liang-Barsky and Cyrus-Beck al… ▽ More A new algorithm for clipping a line segment against a pyramid in E3 is presented. This algorithm avoids computation of intersection points which are not end-points of the output line segment. It also allows solving all cases more effectively. The performance of this algorithm is shown to be consistently better than existing algorithms, including the Cohen-Sutherland, Liang-Barsky and Cyrus-Beck algorithms. △ Less

Submitted 3 January, 2022; originally announced January 2022.

MSC Class: 68U05 ACM Class: I.3

Journal ref: Machine GRAPHICS & VISION, Poland Academy of Sciences, Vol. 9, No. 4, 2000, pp. 841-850, ISSN 1230-0535, 2000

A Comparison of O(1) and Cyrus-Beck Line Clipping Algorithms in E2 and E3

Authors: Vaclav Skala, Pavel Lederbuch, Bohumir Sup

Abstract: A comparison of a new algorithm for line clipping in E2 and E3 by convex polygon and/or polyhedron with O(1) processing complexity and Cyrus- Beck algorithm is presented. The new algorithm in E2 is based on dual space representation and space subdivision technique. The principle of algorithm in E3 is based on the projection of polyhedron to three orthogonal E2 coordinate systems. Algorithms have o… ▽ More A comparison of a new algorithm for line clipping in E2 and E3 by convex polygon and/or polyhedron with O(1) processing complexity and Cyrus- Beck algorithm is presented. The new algorithm in E2 is based on dual space representation and space subdivision technique. The principle of algorithm in E3 is based on the projection of polyhedron to three orthogonal E2 coordinate systems. Algorithms have optimal complexities O(1) and demonstrates that preprocessing can be used to speed up the line clipping significantly. Obvious applications are for one polygon and/or polyhedron and many clipped lines. Detailed theoretical estimations and experimental results are also presented. △ Less

Submitted 15 November, 2021; originally announced November 2021.

Big Geo Data Surface Approximation using Radial Basis Functions: A Comparative Study

Authors: Zuzana Majdisova, Vaclav Skala

Abstract: Approximation of scattered data is often a task in many engineering problems. The Radial Basis Function (RBF) approximation is appropriate for big scattered datasets in $n-$dimensional space. It is a non-separable approximation, as it is based on the distance between two points. This method leads to the solution of an overdetermined linear system of equations. In this paper the RBF approximation… ▽ More Approximation of scattered data is often a task in many engineering problems. The Radial Basis Function (RBF) approximation is appropriate for big scattered datasets in $n-$dimensional space. It is a non-separable approximation, as it is based on the distance between two points. This method leads to the solution of an overdetermined linear system of equations. In this paper the RBF approximation methods are briefly described, a new approach to the RBF approximation of big datasets is presented, and a comparison for different Compactly Supported RBFs (CS-RBFs) is made with respect to the accuracy of the computation. The proposed approach uses symmetry of a matrix, partitioning the matrix into blocks and data structures for storage of the sparse matrix. The experiments are performed for synthetic and real datasets. △ Less

Submitted 20 June, 2018; originally announced June 2018.

Comments: arXiv admin note: text overlap with arXiv:1806.04243

Journal ref: Computers and Geosciences, Vol.109, pp.51-58, ISSN 0098-3004, Elsevier, December 2017

A New Radial Basis Function Approximation with Reproduction

Authors: Zuzana Majdisova, Vaclav Skala

Abstract: Approximation of scattered geometric data is often a task in many engineering problems. The Radial Basis Function (RBF) approximation is appropriate for large scattered (unordered) datasets in d-dimensional space. This method is useful for a higher dimension d>=2, because the other methods require a conversion of a scattered dataset to a semi-regular mesh using some tessellation techniques, which… ▽ More Approximation of scattered geometric data is often a task in many engineering problems. The Radial Basis Function (RBF) approximation is appropriate for large scattered (unordered) datasets in d-dimensional space. This method is useful for a higher dimension d>=2, because the other methods require a conversion of a scattered dataset to a semi-regular mesh using some tessellation techniques, which is computationally expensive. The RBF approximation is non-separable, as it is based on a distance of two points. It leads to a solution of overdetermined Linear System of Equations (LSE). In this paper a new RBF approximation method is derived and presented. The presented approach is applicable for d dimensional cases in general. △ Less

Submitted 18 April, 2018; originally announced April 2018.

Journal ref: Proceedings of IHCI 2016; GET 2016; and CGVCVIP 2016, Portugal, pp.215-222, ISBN 978-989-8533-52-4, IADIS Press, 2016

Least Square Error Method Robustness of Computation: What is not usually considered and taught

Authors: Vaclav Skala

Abstract: There are many practical applications based on the Least Square Error (LSE) approximation. It is based on a square error minimization 'on a vertical' axis. The LSE method is simple and easy also for analytical purposes. However, if data span is large over several magnitudes or non-linear LSE is used, severe numerical instability can be expected. The presented contribution describes a simple method… ▽ More There are many practical applications based on the Least Square Error (LSE) approximation. It is based on a square error minimization 'on a vertical' axis. The LSE method is simple and easy also for analytical purposes. However, if data span is large over several magnitudes or non-linear LSE is used, severe numerical instability can be expected. The presented contribution describes a simple method for large span of data LSE computation. It is especially convenient if large span of data are to be processed, when the 'standard' pseudoinverse matrix is ill conditioned. It is actually based on a LSE solution using orthogonal basis vectors instead of orthonormal basis vectors. The presented approach has been used for a linear regression as well as for approximation using radial basis functions. △ Less

Submitted 1 January, 2018; originally announced February 2018.

O(lgN) Line Clipping Algorithm in E2

Authors: Vaclav Skala

Abstract: A new O(lg N) line clipping algorithm in E2 against a convex window is presented. The main advantage of the presented algorithm is the principal acceleration of the line clipping problem solution. A comparison of the proposed algorithm with others shows a significant improvement in run-time. Experimental results for selected known algorithms are also shown. A new O(lg N) line clipping algorithm in E2 against a convex window is presented. The main advantage of the presented algorithm is the principal acceleration of the line clipping problem solution. A comparison of the proposed algorithm with others shows a significant improvement in run-time. Experimental results for selected known algorithms are also shown. △ Less

Submitted 1 January, 2018; originally announced January 2018.

A Fast Algorithm for Line Clipping by Convex Polyhedron in E3

Authors: Vaclav Skala

Abstract: A new algorithm for line clipping against convex polyhedron is given. The suggested algorithm is faster for higher number of facets of the given polyhedron than the traditional Cyrus-Beck's and others algorithms with complexity O(N) . The suggested algorithm has O(N) complexity in the worst N case and expected O(sqrt(N))) complexity. The speed up is achieved because of 'known order' of triangles.… ▽ More A new algorithm for line clipping against convex polyhedron is given. The suggested algorithm is faster for higher number of facets of the given polyhedron than the traditional Cyrus-Beck's and others algorithms with complexity O(N) . The suggested algorithm has O(N) complexity in the worst N case and expected O(sqrt(N))) complexity. The speed up is achieved because of 'known order' of triangles. Some principal results of comparisons of selected algorithms are presented and give some imagination how the proposed algorithm could be used effectively. △ Less

Submitted 1 January, 2018; originally announced January 2018.

A Comparative Study of LOWESS and RBF Approximations for Visualization

Authors: Michal Smolik, Vaclav Skala, Ondrej Nedved

Abstract: Approximation methods are widely used in many fields and many techniques have been published already. This comparative study presents a comparison of LOWESS (Locally weighted scatterplot smoothing) and RBF (Radial Basis Functions) approximation methods on noisy data as they use different approaches. The RBF approach is generally convenient for high dimensional scattered data sets. The LOWESS metho… ▽ More Approximation methods are widely used in many fields and many techniques have been published already. This comparative study presents a comparison of LOWESS (Locally weighted scatterplot smoothing) and RBF (Radial Basis Functions) approximation methods on noisy data as they use different approaches. The RBF approach is generally convenient for high dimensional scattered data sets. The LOWESS method needs finding a subset of nearest points if data are scattered. The experiments proved that LOWESS approximation gives slightly better results than RBF in the case of lower dimension, while in the higher dimensional case △ Less

Submitted 1 January, 2018; originally announced January 2018.

A two-level approach to implicit surface modeling with compactly supported radial basis functions

Authors: Rongjiang Pan, Vaclav Skala

Abstract: We describe a two-level method for computing a function whose zero-level set is the surface reconstructed from given points scattered over the surface and associated with surface normal vectors. The function is defined as a linear combination of compactly supported radial basis functions (CSRBFs). The method preserves the simplicity and efficiency of implicit surface interpolation with CSRBFs and… ▽ More We describe a two-level method for computing a function whose zero-level set is the surface reconstructed from given points scattered over the surface and associated with surface normal vectors. The function is defined as a linear combination of compactly supported radial basis functions (CSRBFs). The method preserves the simplicity and efficiency of implicit surface interpolation with CSRBFs and the reconstructed implicit surface owns the attributes, which are previously only associated with globally supported or globally regularized radial basis functions, such as exhibiting less extra zero-level sets, suitable for inside and outside tests. First, in the coarse scale approximation, we choose basis function centers on a grid that covers the enlarged bounding box of the given point set and compute their signed distances to the underlying surface using local quadratic approximations of the nearest surface points. Then a fitting to the residual errors on the surface points and additional off-surface points is performed with fine scale basis functions. The final function is the sum of the two intermediate functions and is a good approximation of the signed distance field to the surface in the bounding box. Examples of surface reconstruction and set operations between shapes are provided. △ Less

Submitted 28 July, 2017; originally announced August 2017.

Journal ref: Engineering with Computers 2011 27:299-307

Fractions, Projective Representation, Duality, Linear Algebra and Geometry

Authors: Vaclav Skala

Abstract: This contribution describes relationship between fractions, projective representation, duality, linear algebra and geometry. Many problems lead to a system of linear equations. This paper presents equivalence of the Cross-product operation and solution of a system of linear equations Ax=0 or Ax=b using projective space representation and homogeneous coordinates. It leads to conclusion that divisio… ▽ More This contribution describes relationship between fractions, projective representation, duality, linear algebra and geometry. Many problems lead to a system of linear equations. This paper presents equivalence of the Cross-product operation and solution of a system of linear equations Ax=0 or Ax=b using projective space representation and homogeneous coordinates. It leads to conclusion that division operation is not required for a solution of a system of linear equations, if the projective representation and homogeneous coordinates are used. An efficient solution on CPU and GPU based architectures is presented with an application to barycentric coordinates computation as well. △ Less

Submitted 28 July, 2017; originally announced August 2017.

Journal ref: CZECH-SLOVAK CONFERENCE ON GEOMETRY AND GRAPHICS 2016, ISBN 978-80-7464-874-8

Space Subdivision to Speed-up Convex Hull Construction in E3

Authors: Vaclav Skala, Zuzana Majdisova, Michal Smolik

Abstract: Convex hulls are fundamental geometric tools used in a number of algorithms. This paper presents a fast, simple to implement and robust Smart Convex Hull (S-CH) algorithm for computing the convex hull of a set of points in E3. This algorithm is based on "spherical" space subdivision. The main idea of the S-CH algorithm is to eliminate as many input points as possible before the convex hull constru… ▽ More Convex hulls are fundamental geometric tools used in a number of algorithms. This paper presents a fast, simple to implement and robust Smart Convex Hull (S-CH) algorithm for computing the convex hull of a set of points in E3. This algorithm is based on "spherical" space subdivision. The main idea of the S-CH algorithm is to eliminate as many input points as possible before the convex hull construction. The experimental results show that only a very small number of points are used for the final convex hull calculation. Experiments made also proved that the proposed S-CH algorithm achieves a better time complexity in comparison with other algorithms in E3. △ Less

Submitted 9 August, 2017; originally announced August 2017.

Journal ref: Advances in Engineering Software, Vol.91, pp.12-22, ISSN 0965-9978, Elsevier, January 2016

Fast Algorithm for Finding Maximum Distance with Space Subdivision in E2

Authors: Vaclav Skala, Zuzana Majdisova

Abstract: Finding an exact maximum distance of two points in the given set is a fundamental computational problem which is solved in many applications. This paper presents a fast, simple to implement and robust algorithm for finding this maximum distance of two points in E2. This algorithm is based on a polar subdivision followed by division of remaining points into uniform grid. The main idea of the algori… ▽ More Finding an exact maximum distance of two points in the given set is a fundamental computational problem which is solved in many applications. This paper presents a fast, simple to implement and robust algorithm for finding this maximum distance of two points in E2. This algorithm is based on a polar subdivision followed by division of remaining points into uniform grid. The main idea of the algorithm is to eliminate as many input points as possible before finding the maximum distance. The proposed algorithm gives the significant speed up compared to the standard algorithm. △ Less

Submitted 9 August, 2017; originally announced August 2017.

Journal ref: ICIG 2015 Proceedings Part II, China, pp.261-274, Springer, 2015

Continuous Global Optimization in Surface Reconstruction from an Oriented Point Cloud

Authors: Rongjiang Pan, Vaclav Skala

Abstract: We introduce a continuous global optimization method to the field of surface reconstruction from discrete noisy cloud of points with weak information on orientation. The proposed method uses an energy functional combining flux-based data-fit measures and a regularization term. A continuous convex relaxation scheme assures the global minima of the geometric surface functional. The reconstructed sur… ▽ More We introduce a continuous global optimization method to the field of surface reconstruction from discrete noisy cloud of points with weak information on orientation. The proposed method uses an energy functional combining flux-based data-fit measures and a regularization term. A continuous convex relaxation scheme assures the global minima of the geometric surface functional. The reconstructed surface is implicitly represented by the binary segmentation of vertices of a 3D uniform grid and a triangulated surface can be obtained by extracting an appropriate isosurface. Unlike the discrete graph-cut solution, the continuous global optimization entails advantages like memory requirements, reduction of metrication errors for geometric quantities, allowing globally optimal surface reconstruction at higher grid resolutions. We demonstrate the performance of the proposed method on several oriented point clouds captured by laser scanners. Experimental results confirm that our approach is robust to noise, large holes and non-uniform sampling density under the condition of very coarse orientation information. △ Less

Submitted 28 July, 2017; originally announced July 2017.

Journal ref: Computer Aided Design, Vol.43, No.8, pp.896-901, Elsevier, ISSN 0010-4485, 2011