Polygon Morphing and Its Application in Orebody Modeling

Three different polygon morphing methods are examined. The first one is based on the utilization of the trimmed skeleton of the symmetric difference of the source and target polygons as an intermediate polygon. The second one reduces the problem to the problem of morphing compatible planar triangulations and utilizes the representation of planar triangulations as a matrix constructed using barycentric coordinates of the planar triangulation’s vertices relative to their neighbors. The third and last one describes the polygon by the parametric curve representation based on estimated Fourier parameters and thus transfers the morphing process to Fourier parametric space. The different features and comparative results of these methods are shown by the tests with different examples. Thesemethods are used for generating a set of polygonal sections from two nonplanar polygonal sections which are nearly planar in 3D before constructing a threedimensional object from these nonplanar sections.


Introduction
Morphing or metamorphosis is usually defined as the gradual, smooth, and continuous transformation of a source object into a target object 1, 2 . Morphing is used in a variety of applications such as computer graphics, computer animation, industrial design, geometric modelling, and medicine.
The morphing problem has been examined in many domains, such as two-dimensional images 3-5 , polygons and polylines 1, 6-12 , free-form curves 13 , polyhedra and volumetric representations 6, 14, 15 . The morphing process for polygons and polyhedras consists of solving two main problems. The first one, known as the vertex correspondence problem, is to generate a correspondence between geometric elements of the source and target objects. The second one, known as the vertex path problem, is to find paths that corresponding geometric elements travel during the morphing process. In this work, three different morphing methods for polygons with some correspondence between their vertices are compared, and a morphing method for nonplanar polygonal sections in 3D without correspondence between their vertices is generated. Several approaches have been described on the polygon morphing problem. Shapira and Rappoport 11 presented a polygon morphing technique that takes into account both the interiors and the boundaries of the polygons. This is achieved by representing the interiors of the polygons with compatible skeletons which are called star skeletons and applying the morphing to the parametric description of these skeletons. The generation of corresponding skeletons of two polygons is used, and though theoretically possible for all shapes, it seems natural for similar shapes, but uncertain for rather different shapes. Alexa et al. 1 morphed the interiors of the polygons rather than their boundaries to achieve locally least distorting intermediate polygons. They compatibly triangulated the source and target polygons and developed the vertex paths during the morph in an attempt to maintain compatibility and preserve shape. Firstly they apply an algorithm to divide the source and target polygons into separate isomorphic triangles. Then, they develop a method for interpolating the locations of corresponding boundary and interior vertices along their paths from the source to the target polygons. Blanding et al. 6 described the skeleton-based polygon morphing and then generalized the same procedure directly to three dimensions. They used the trimmed skeleton of the symmetric difference of the source and target polygons as intermediate polygons. Gotsman and Surazhsky 9 described polygon morphing using compatible triangulation to reduce the polygon morphing problem to the problem of morphing compatible planar triangulations and utilized the representation of planar triangulations as a matrix constructed using barycentric coordinates of the planar triangulation's vertices relative to their neighbors. Chen and Sun 7 described polygon morphing using Fourier parametrization to define the polygon by the parametric curve representation based on estimated Fourier parameters and thus transfer the morphing process to Fourier parametric space. In this work, the last three polygon morphing methods mentioned above are examined, and the last one is used for morphing nonplanar sections which are nearly planar in 3D.

Polygon Morphing
Given a source polygon P 0 and a target polygon P 1 with a correspondence between their vertices, the morphing procedure can be defined by the continuous function P : t → P t , 0 ≤ t ≤ 1, where P 0 P 0 and P 1 P 1 are the source and target polygons, and for all 0 < t < 1, P t P t are intermediate polygons. Additionally, the intermediate polygon P t P t is similar to P 0 P 0 as t → 0 and similar to P 1 P 1 as t → 1.

Skeleton-Based Polygon Morphing
Skeleton-based polygon morphing method 6 is based on utilization of the trimmed skeleton of the symmetric difference of the source and target polygons as intermediate polygons.
Each boundary point of the intermediate polygon generated by the morphing process can be defined as points equidistant from the nearest points of the source and target polygons. Similarly, skeletons consist of points minimally equidistant from two boundary points of a given shape. Therefore the skeleton of a shape that incorporates the boundaries of both the source and target polygons can be used for constructing middle polygon generated on the morphing process. Figure 1: The main steps of skeleton-based polygon morphing. Two rectangles P 0 and P 1 are combined by taking their symmetric difference a and the resulting shape is skeletonized b and trimmed c to obtain the intermediate polygon P 1/2 .

Mathematical Problems in Engineering
In this method, firstly the intermediate polygon P 1/2 P 1/2 is computed by using the source polygon P 0 P 0 and the target polygon P 1 P 1 as shown in Figure 1. To achieve this, first the skeleton of the symmetric difference of a combination of the source and target polygons is computed. Then, the intermediate polygon P 1/2 P 1/2 is generated from the skeleton by applying a trimming operation to the skeleton. Similiary, the intermediate polygon P 1/4 P 1/4 is generated by using the polygons P 0 P 0 and P 1/2 P 1/2 , and the intermediate polygon P 3/4 P 3/4 is generated by using the polygons P 1/2 P 1/2 and P 1 P 1 . The same process is applied recursively until a desired sequence of intermediate polygons is produced.

Polygon Morphing Using Compatible Triangulation
Polygon morphing method using compatible triangulation 9 reduces the problem to the problem of morphing compatible planar triangulations.
The first task in polygon morphing using compatible triangulation is to compatibly triangulate the interiors of the given source polygon P 0 and the target polygon P 1 with the matching vertices. Then the polygons P 0 and P 1 are enclosed by identical convex polygons Q 0 and Q 1 , and subsequently the annuli between Q 0 and P 0 and between Q 1 and P 1 are compatibly triangulated as shown in Figure 2. These two compatible triangulations with common convex boundary are represented by matrices constructed using barycentric coordinates of the these triangulation's vertices relative to their neighbors. The next step is to generate the convex combination morph between these two compatible triangulations with common convex boundary to produce desired intermediate compatible triangulation with common convex boundary. Finally, all edges of these intermediate triangulations which do not belong to polygons P 0 and P 1 are discarded to obtain the desired intermediate polygon.
Given polygons P 0 and P 1 with matching vertices, an overview of all the steps of the method of polygon morphing using compatible triangulation is shown in Table 1.

Polygon Morphing Using Fourier Parametrization
Polygon morphing method using Fourier parametrization 7 describes the polygon by the parametric curve representation based on estimated Fourier parameters and thus transfers the morphing process to Fourier parametric space.
In this method, firstly the polygon coordinates are transformed to local coordinates, and then polygons are divided into matching segments using the correspondence between their vertices. The Fourier parametrization for each matching segment is computed, and the  The main steps of the method of morphing polygons P 0 and P 1 using compatible triangulation.
1 Compatibly triangulate the interiors of P 0 and P 1 with the matching vertices. 2 Enclose P 0 and P 1 by identical convex polygons Q 0 and Q 1 , respectively. 3 Compatibly triangulate the annuli between Q 0 and P 0 and between Q 1 and P 1 , respectively. 4 Morph the two compatible triangulations using the method of morphing compatible triangulation. 5 Discard all edges of the two triangulations not belonging to P 0 and P 1 .

Results and Discussion
The morphing problem can be solved by using a lot of different methods, but some constraints or desirable properties can be given for morphing procedure. Firstly, if the target polygon is a translation or rotation of the source polygon, then the intermediate polygons should also be an appropriate translation or rotation of the source polygon. Secondly, the area of the polygon generated by morphing should change smoothly as the source polygon turns into the target polygon. As Figure 3 shows, among the examined three methods, polygon morphing method using Fourier parametrization is the only one which guarantees these properties. Thirdly, if the source and target polygons are closed and simple, then the intermediate polygons should also be closed and simple. Figure 3 shows that all of the examined methods achieve this property. The comparative results of these methods are shown in the following experiments. The results of all methods seem natural for similar shapes, but the results of skeleton-based polygon morphing method seem ambiguous for rather different shapes.
As a distinction from the other two methods, skeleton-based polygon morphing method is capable of morphing between two polygons with different topologies, as shown in

Morphing of Nonplanar Sections
In mining, for 3D orebody modelling, three-dimensional objects are generated from nonplanar sections by building 3D triangulated or tetrahedral models. In the case where nonplanar sections are away from each other, three-dimensional models cannot be smooth as shown in Figure 5. To build higher quality 3D models, these nonplanar sections should be morphed to generate intermediate sections. Polygon morphing method using Fourier parametrization is utilized for morphing two nonplanar sections which are nearly planar. To morph nonplanar sections which are nearly planar, firstly the source and target nonplanar sections are triangulated separately. The nonplanar section is specified in 3D, even though the triangulation is 2D. Thus, the nonplanar section is transformed in a plane, the transformed section is triangulated, and the topology of this triangulation is used for the original nontransformed section. To find a plane for triangulation, the most appropriate face of the oriented bounding box of the the nonplanar section can be used. These steps of the triangulation of the nonplanar section are shown in Figure 6. · · · n t 1 /t. Thus, these transformed planar sections can be specified in 2D. Then, each transformed section in 2D is separately parametrized proportionally to arclength as g : 0, 1 → R 2 , t → g t x t , y t , and the Fourier parametrization for each curve is separately computed as where m 0 and m 1 are the center of the planes D 0 and D 1 , respectively. Thus the intermediate sections are obtained. An overview of all these steps to morph nonplanar sections is shown in Table 3. Figure 7 shows a set of intermediate sections generated from two nonplanar polygonal sections in 3d, and Figure 8 shows a three-dimensional object of the source, target and intermediate sections. the transformations of nonplanar sections can have an important effect on the results of morphing nonplanar section, namely this method relies on having data organized in nearly planar sections.

Conclusion
We have presented the different features and comparative results of three different morphing methods which are skeleton-based polygon morphing, polygon morphing using compatible triangulation, and polygon morphing using Fourier parametrization. According to results, we developed a generalization of the intermediate sections from two nonplanar polygonal  sections in 3D which are nearly planar, before constructing a three-dimensional object from these nonplanar polygonal sections to build higher quality 3D models. To achieve this, firstly we needed to solve the vertex path problem. In the skeleton-based polygon morphing method, the source and target objects are combined into a single object as desired instead of generating a correspondence between the geometric elements of the source and target objects. This process is important and has a major effect on the morph history generated. But in mining, there are nonplanar sections with many different initial relative positions and orientations in 3D. Therefore transformed combination of two nonplanar sections in 2D can be constructed, and polygon morphing method using Fourier parametrization is utilized for morphing these two 2D sections, and so intermediate sections are obtained in 2D. Then these intermediate sections are transformed to appropriate location between the source and target nonplanar polygonal sections. In this manner, intermediate polygonal sections in 3D are generated on the morphing process of the source and target nonplanar polygonal sections in 3D. Our method relies on having data organized in nearly planar sections, since transformation of folded sections in 2D cannot be simple, and the method is based on the transformations of nonplanar sections.