Next Article in Journal
EvolveNet: Evolving Networks by Learning Scale of Depth and Width
Previous Article in Journal
Preface to the Special Issue on “Mathematical Methods for Computer Science”
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Birational Quadratic Planar Maps with Generalized Complex Rational Representations †

1
Department of Mathematics, Hohai University, Nanjing 211100, China
2
School of Physical and Mathematical Sciences, Nanjing Tech University, Nanjing 211816, China
3
Department of Computer Science, Rice University, Houston, TX 77251, USA
*
Author to whom correspondence should be addressed.
This paper is an extension of our paper published in Wang, X.; Wu, M.; Liu, Y.; Ni, Q. Constructing quadratic birational maps via their complex rational representation. Comput. Aided Geom. Desi. 2021, 85, 101969.
Mathematics 2023, 11(16), 3609; https://doi.org/10.3390/math11163609
Submission received: 28 July 2023 / Revised: 18 August 2023 / Accepted: 19 August 2023 / Published: 21 August 2023
(This article belongs to the Section Computational and Applied Mathematics)

Abstract

:
Complex rational maps have been used to construct birational quadratic maps based on two special syzygies of degree one. Similar to complex rational curves, rational curves over generalized complex numbers have also been constructed by substituting the imaginary unit with a new independent quantity. We first establish the relationship between degree one, generalized, complex rational Bézier curves and quadratic rational Bézier curves. Then we provide conditions to determine when a quadratic rational planar map has a generalized complex rational representation. Thus, a rational quadratic planar map can be made birational by suitably choosing the middle Bézier control points and their corresponding weights. In contrast to the edges of complex rational maps of degree one, which are circular arcs, the edges of the planar maps can be generalized to hyperbolic and parabolic arcs by invoking the hyperbolic and parabolic numbers.

1. Introduction

A birational planar map is a rational function p ( s , t ) from R 2 to R 2 with a rational function p 1 ( x , y ) inverse to p ( s , t ) , i.e., p 1 p ( s , t ) = ( s , t ) . The research on birational maps traces its history back to 1863 when birational maps were known as Cremona transformations. H. Hudson published a monograph in 1927 to classify Cremona transformations in the plane and in 3D space [1]. However, [1] predates by over half a century Bézier and B-spline curves and surfaces.
Free-form deformation (FFD) is an artist-friendly tool, which is used to model deformations by warping the ambient space [2]. FFD has numerous applications such as image morphing [3,4], facial animation [5] and medical image registration [6].
A preimage of a given ( x , y ) is defined as any ( s , t ) for which p ( s , t ) = ( x , y ) . Computing FFD preimages is also required in many applications, e.g., extended FFD [7] and directly manipulated FFD [8]. However, computing preimages usually involves solving a system of nonlinear equations using either numerical methods or algebraic techniques. The computation of preimages can be simplified when the rational map used in FFD is birational. Hence, the notion of birational FFD is put forward in [9], i.e., an FFD in which the control points can be moved freely and the FFD can be made birational by adjusting the weights.
Birational FFDs of degree 1 × 1 , i.e., birational quadrilateral maps, are constructed by considering axial moving lines in [10]. The inverses of birational quadrilateral maps are also discussed from the viewpoint of generalized barycentric coordinates [11,12]. The relationship between the inverses of rational quadrilateral maps and their characteristic conic is presented in [13]. Hence, by setting the isoparametric curves to be monoid curves, Seberberg et al. [9] studied the construction of the birational 2D FFD of degree 1 × n . To characterize the birational maps, Botbol et al. [14] provides a birationality criterion for bigraded birational maps. Note that the current methods to construct birational maps focus on the discussion of 2D rational maps related to bidegree cases [9,10].
Complex rational curves were introduced by J. Sánchez-Reyes by allowing complex weights in rational Bézier curves [15]. Based on the observation that complex rational curves have two special syzygies of low degrees [16], we can construct birational quadratic planar maps by exploiting their complex rational representations [17], which are complex rational maps of degree one. Complex rational planar maps of degree one represent exactly Möbius transformations. Based on Möbius transformations, Lipman et al. also introduced the four-point interpolant (FPI) maps, where four endpoints are free to move. Since Möbius transformations map lines to circles, the boundaries of the complex rational planar maps of degree one and FPI are limited circular arcs. Recently, the planar and bent corner-operated and transimilar (COTS) maps have been proposed to take the unit square to a region bounded by four log-spiral edges [18,19]. Based on the self-similar property of the log-spiral, the tiles of the COTS map of a regular pattern are similar, which will reduce the query cost of point inclusion testing (PIT) [20] and total area calculation (TAC). To increase the flexibility to construct birational planar maps and inspired by recent work on generalized complex rational curves [21], we shall discuss the construction of birational quadratic planar maps by probing their generalized complex representations. The main contributions of this study are as follows:
  • We introduce generalized complex rational curves and planar maps, and we discuss their properties.
  • We prove that degree one, generalized, complex rational maps correspond to quadratic birational planar maps.
  • We propose a new method to construct quadratic birational planar maps via generalized complex rational planar maps.
This paper is organized in the following fashion. We begin in Section 2 by recalling some preliminary results about rational curves over generalized complex numbers. In Section 3, we establish the relationship between quadratic rational Bézier curves and degree one complex rational curves. By considering the bivariate case of generalized complex rational curves, we construct birational quadratic planar maps by degree one, generalized, complex rational maps in Section 4. We close in Section 5 with a brief summary of our work.

2. Rational Curves over Generalized Complex Numbers

2.1. Complex, Hyperbolic and Parabolic Numbers

Hyperbolic and parabolic numbers generalize complex numbers by introducing a new quantity α : either α 2 = 1 (complex numbers), or α 2 = + 1 (hyperbolic numbers), or α 2 = 0 (parabolic numbers). The quantity α here is not a real number but an independent quantity. In addition, α commutes with every real number, and multiplication is associative and distributes through addition. Thus, we have the following multiplication rules:
( a + b α ) ( x + y α ) = ( a x + b y α 2 ) + ( a y + b x ) α = ( a x b y ) + ( a y + b x ) α , α 2 = 1 ; ( a x + b y ) + ( a y + b x ) α , α 2 = 1 , α ± 1 ; ( a x ) + ( a y + b x ) α , α 2 = 0 , α 0 .
The conjugate of a + b α is defined by
a + b α ¯ = a b α .
Multiplying a generalized complex number by its conjugate yields a real number, i.e., ( a + b α ¯ ) ( a + b α ) = a 2 b 2 α 2 . The modulus (or absolute value) of a + b α is defined by
| a + b α | = | ( a + b α ¯ ) ( a + b α ) | = | a 2 b 2 α 2 | .
For generalized complex numbers, the analog of Euler’s formula is also true [22]. We have the following exponential formulas:
  • (Complex numbers) e α t = cos ( t ) + α sin ( t ) ;
  • (Hyperbolic numbers) e α t = cosh ( t ) + α sinh ( t ) ;
  • (Parabolic numbers) e α t = 1 + α t .
The exponential of generalized complex numbers represents either a rotation for complex numbers, a scissor shear for hyperbolic numbers [23] or a classical shear for parabolic numbers [24].
A generalized complex number can be identified with a point in the plane R 2 . Using the standard basis { 1 , α } , we identify a + b α with the point or vector ( a , b ) . In analogy with complex numbers, we shall call the coefficient of 1 the real component and the coefficient of α the imaginary component of a hyperbolic or parabolic number.
The hyperbolic numbers can be used to help solve cubic Equations [25] and write Lorentz transformations in the special theory of relativity [22]. The parabolic numbers can be used to help analyze the mechanics of a mechanism [26].

2.2. Complex, Hyperbolic and Parabolic Rational Curves and Planar Maps

Generalized complex Bézier curves are defined by
c ( s ) = n ( s ) d ( s ) = i = 0 n p i w i B i n ( s ) i = 0 n w i B i n ( s ) , s [ 0 , 1 ] ,
where p i = x i + y i α are generalized complex control points, w i = u i + v i α are generalized complex weights and B i n ( s ) = n ! i ! ( n i ) ! s i ( 1 s ) n i are Bernstein polynomials. c ( s ) is called a complex, hyperbolic or parabolic rational Bézier curve when α 2 = 1 , + 1 , 0 , respectively. By separating the real and imaginary parts of (1), we obtain
c ( s ) = n 0 ( s ) + α n 1 ( s ) d 0 ( s ) + α d 1 ( s ) ,
where n 0 ( s ) , n 1 ( s ) , d 0 ( s ) , d 1 ( s ) R [ s ] . Multiplying the numerator and denominator in (1) by the conjugate of the denominator yields
c ( s ) = n 0 d 0 α 2 n 1 d 1 + α n 1 d 0 n 0 d 1 d 0 2 α 2 d 1 2 .
Thus, c(s) corresponds to a real rational curve in the homogeneous form:
C ( s ) = n 0 d 0 α 2 n 1 d 1 , n 1 d 0 n 0 d 1 , d 0 2 α 2 d 1 2 .
Note that when n 0 ( s ) n 1 ( s ) , d 0 ( s ) , d 1 ( s ) are polynomials of degree n, i.e., c(s) is a generalized complex rational curve of degree n, its counterpart C ( s ) = n 0 d 0 α 2 n 1 d 1 , n 1 d 0 n 0 d 1 , d 0 2 α 2 d 1 2 is a real rational curve of degree 2 n .
To discuss rational planar maps, we extend the generalized complex rational curves c ( s ) in 2D as follows:
c ( s , t ) = n 0 ( s , t ) + α n 1 ( s , t ) d 0 ( s , t ) + α d 1 ( s , t ) ,
where n 0 ( s , t ) , n 1 ( s , t ) , d 0 ( s , t ) , d 1 ( s , t ) R [ s , t ] . Recall that a syzygy of Q ( t 1 , t 2 , , t m ) R k [ t 1 , t 2 , , t m ] is a vector of polynomials
L ( t 1 , t 2 , , t m ) = l 0 ( t 1 , t 2 , , t m ) , l 1 ( t 1 , t 2 , , t m ) , , l k 1 ( t 1 , t 2 , , t m )
such that
L ( t 1 , t 2 , , t m ) · Q ( t 1 , t 2 , , t m ) 0 .
Following the results in [21], we have the following observation
Lemma 1.
Given a rational planar map C ( s , t ) = ( n 0 d 0 α 2 n 1 d 1 , n 1 d 0 n 0 d 1 , d 0 2 α 2 d 1 2 ) with generalized complex rational representation (4),
L 0 ( s , t ) = ( d 0 ( s , t ) , α 2 d 1 ( s , t ) , n 0 ( s , t ) ) , L 1 ( s , t ) = ( d 1 ( s , t ) , d 0 ( s , t ) , n 1 ( s , t ) ) .
are two syzygies of C ( s , t ) . Moreover, L 0 ( s , t ) × L 1 ( s , t ) = C ( s , t ) .
Proof. 
These results follow directly by substituting (5) into L i ( s , t ) · C ( s , t ) , i = 0 , 1 and L 0 ( s , t ) × L 1 ( s , t ) .    □
Note that the set of all syzygies of C ( s , t ) is also a free module of rank two over the ring R [ s , t ] (for more details, the interested reader is referred to the proof of Proposition 2.1 in [27]). Moreover, μ -bases of generalized complex rational curves can be calculated via these syzygies in the one-dimensional case [16,21].

3. Generalized Complex Rational Linear Curves

Let p l ( s ) denote the generalized complex rational Bézier curve of degree one:
p l ( s ) = w 0 p 0 ( 1 s ) + w 1 p 1 s w 0 ( 1 s ) + w 1 s , s [ 0 , 1 ] ,
where p i = x i + α y i , w i = u i + α v i , i = 0 , 1 . Note that when w 0 , w 1 R , p l ( s ) will degenerate to a line segment connecting p 0 and p 1 . Hence, in the rest of this paper, we assume at least one of w 0 , w 1 R .
Lemma 2.
Let
p ( s ) = p 0 w 0 B 0 2 ( s ) + p * w * B 1 2 ( s ) + p 1 w 1 B 2 2 ( s ) w 0 B 0 2 ( s ) + w * B 1 2 ( s ) + w 1 B 2 2 ( s )
be the corresponding real rational curve of p l ( s ) . Then, the control points and the weights in (7) have the following properties:
(1) 
w 0 = w 0 w ¯ 0 , w * = ( w 0 w ¯ 1 + w 1 w ¯ 0 ) / 2 , w 1 = w 1 w ¯ 1 , w 0 , w 1 , w 2 R ;
(2) 
p * = ( w 0 w ¯ 1 p 0 + w 1 w ¯ 0 p 1 ) / ( w 0 w ¯ 1 + w ¯ 0 w 1 ) = ( w 0 w ¯ 1 p 0 + w 1 w ¯ 0 p 1 ) / ( 2 w * ) ;
(3) 
p 1 p * = w 0 w ¯ 1 2 w * ( p 1 p 0 ) , p * p 0 = w ¯ 0 w 1 2 w * ( p 1 p 0 ) ;
(4) 
( p 1 p * ) / ( p 1 p 0 ) is the conjugate of ( p * p 0 ) / ( p 1 p 0 ) ;
(5) 
p * = ( p 0 + p 1 ) / 2 + w ¯ 0 w 1 w 0 w ¯ 1 4 w * ( p 1 p 0 ) = ( p 0 + p 1 ) / 2 + λ * ( p 1 p 0 ) α , where
λ * = u 1 v 0 u 0 v 1 2 ( u 0 u 1 + α 2 v 0 v 1 ) ;
(6) 
Geometrically, p * is on three different lines when α 2 is set to 1 , + 1 , 0 . All three lines pass through a common point ( p 0 + p 1 ) / 2 . More precisely,
(a) 
When α 2 = 1 , the line is the perpendicular bisector of the line segment between p 0 and p 1 ;
(b) 
When α 2 = 1 , the line is symmetric to the line p 0 p 1 about the axis ( p 0 + p 1 ) / 2 + ( 1 + α ) t , t R ;
(c) 
When α 2 = 0 , the line is parallel to the y-axis, which also passes through the midpoint ( p 0 + p 1 ) / 2 .
(7) 
The weight w * satisfies
( w * ) 2 = w 0 w 1 1 4 α 2 ( λ * ) 2 .
Proof. 
Items (1), (2), and (3) can be verified by multiplying the denominator and numerator of (6) by w ¯ 0 ( 1 s ) + w ¯ 1 s .
Since the conjugate of w ¯ 0 w 1 2 w * is w 0 w ¯ 1 2 w * , it follows from (3) that ( p 1 p * ) / ( p 1 p 0 ) is the conjugate of ( p * p 0 ) / ( p 1 p 0 ) . Hence, (4) holds.
Again, according to (3),
( p 1 p * ) ( p * p 0 ) = ( p 0 + p 1 ) 2 p * = w 0 w ¯ 1 w ¯ 0 w 1 2 w * ( p 1 p 0 ) .
By solving for p * in (9), we obtain p * = ( p 0 + p 1 ) / 2 + w ¯ 0 w 1 w 0 w ¯ 1 4 w * ( p 1 p 0 ) . Substituting w i = u i + α v i , i = 0 , 1 into w ¯ 0 w 1 w 0 w ¯ 1 4 w * yields
w ¯ 0 w 1 w 0 w ¯ 1 4 w * = ( u 1 v 0 u 0 v 1 ) 2 ( u 0 u 1 + α 2 v 0 v 1 ) α .
By setting λ * = ( u 1 v 0 u 0 v 1 ) / ( 2 ( u 0 u 1 + α 2 v 0 v 1 ) ) , it follows from (9) that p * = ( p 0 + p 1 ) / 2 + λ * ( p 1 p 0 ) α . Thus, (5) holds.
We shall prove (6) based on (5). Since
w ¯ 0 w 1 w 0 w ¯ 1 4 w * ( p 1 p 0 ) = λ * α ( x 1 x 0 ) + ( y 1 y 0 ) α = λ * ( ( y 1 y 0 ) α 2 + ( x 1 x 0 ) α ) ,
corresponds to the vector λ * ( y 1 y 0 ) α 2 , ( x 1 x 0 ) R 2 , it follows that
(a)
When α 2 = 1 , the vector λ * ( y 1 y 0 ) α 2 , x 1 x 0 = λ * ( y 1 y 0 ) , x 1 x 0 , which is perpendicular to the vector ( x 1 x 0 , y 1 y 0 ) . Thus, p * is on the perpendicular bisector of the line segment between p 0 and p 1 ;
(b)
When α 2 = 1 , the vector λ * ( y 1 y 0 ) α 2 , x 1 x 0 = λ * y 1 y 0 , x 1 x 0 . Hence, p * ( p 0 + p 1 ) / 2 = λ * y 1 y 0 , x 1 x 0 . Since the vector y 1 y 0 , x 1 x 0 is symmetric to the vector p 1 p 0 = ( x 1 x 0 , y 1 y 0 ) about the axis y = x , it follows that p * is on a line symmetric to the line p 0 p 1 about the axis ( p 0 + p 1 ) / 2 + t ( 1 , 1 ) , t R ;
(c)
When α 2 = 0 , the vector λ * ( y 1 y 0 ) α 2 , x 1 x 0 = λ * 0 , x 1 x 0 . Hence, p * = ( x 0 + x 1 , y 0 + y 1 ) / 2 + λ * 0 , ( x 1 x 0 ) , which is on a line parallel to the y-axis.
Next, we shall prove (7). Recall from (8) that λ * = u 1 v 0 u 0 v 1 2 ( u 0 u 1 + α 2 v 0 v 1 ) R . On the other hand,
w * = ( w 0 w ¯ 1 + w 1 w ¯ 0 ) / 2 = u 0 u 1 α 2 v 0 v 1 , w 0 w 1 = ( u 0 2 α 2 v 0 2 ) ( u 1 2 α 2 v 1 2 ) .
Hence, when u 0 u 1 α 2 v 0 v 1 0 ,
1 4 α 2 ( λ * ) 2 = 1 α 2 ( u 1 v 0 u 0 v 1 ) 2 ( u 0 u 1 α 2 v 0 v 1 ) 2 = w 0 w 1 ( w * ) 2 .
Thus, ( w * ) 2 = w 0 w 1 1 4 α 2 ( λ * ) 2 and (7) holds.    □
Based on the above analysis, we are now ready to derive conditions for a general real rational quadratic Bézier curve to have a generalized complex rational representation.
Theorem 1.
Suppose that the two endpoints p 0 , p 1 and their weights w 0 , w 1 > 0 are fixed for the rational quadratic curve (7). If the middle control point and its weight satisfy the following two conditions:
(1) 
p * can be expressed as p * = ( p 0 + p 1 ) / 2 + α ( p 1 p 0 ) λ for some λ R ;
(2) 
w * can be expressed as
w * = w 0 w 1 1 4 α 2 λ 2 ,
then p ( s ) has a generalized complex rational representation.
Proof. 
When the middle control point and its corresponding weight satisfy conditions (1) and (2), it can be verified that p ( s ) corresponds to a generalized complex rational representation w 0 p 0 ( 1 s ) + w 0 z 0 p 1 s w 0 ( 1 s ) + w 0 z 0 s if we choose w 0 such that w 0 w ¯ 0 = w 0 and set
w 1 = w 0 z 0 ,
where
z 0 = w * + α ( w * ) 2 w 0 w 1 α 2 / w 0 , α 2 = ± 1 , w * ( 1 + 2 λ α ) / w 0 , α 2 = 0 .
   □
Remark 1.
Based on Theorem 1, a rational quadratic Bézier curve p ( s ) in (7) with complex rational representation is fully determined by p 0 , p 1 , w 0 , w 1 and the parameter λ. To ensure that the curve segment is located in the convex hull of the control polygon p 0 p 1 p 2 , we shall adopt the positive square root for w * in (10). The complementary segment of the curve is generated by reversing of the sign of w * [28].
Remark 2.
When w 0 0 , we can divide the numerator and the denominator of p l ( s ) by w 0 . Hence, p l ( s ) can be rewritten as
p l ( s ) = p 0 ( 1 s ) + w 1 / w 0 p 1 s ( 1 s ) + w 1 / w 0 s = p 0 ( 1 s ) + z 0 p 1 s ( 1 s ) + z 0 s .
Thus, the choice of the ratio z 0 = w 1 / w 0 , instead of w 0 itself, affects the shape of the curve p ( s ) .
Remark 3.
z 0 in (12) can be expressed in exponential form:
z 0 = w 1 w 0 e α arctan ( 2 λ ) , α 2 = 1 ; w 1 w 0 e α artanh ( 2 λ ) , α 2 = 1 ; w 1 w 0 e α ( 2 λ ) , α 2 = 0 .
Remark 4.
In the case of hyperbolic rational curves, since α 2 = 1 , it follows from (10) that the parameter λ must satisfy: λ ( 1 / 2 , 1 / 2 ) .
Example 1.
( α 2 = 1 ) Consider two endpoints p 0 = ( 2 , 1 ) , p 1 = ( 5 , 0 ) and their corresponding weights w 0 = 4 , w 1 = 9 . Choose p * = ( p 0 + p 1 ) / 2 + α ( p 1 p 0 ) λ with λ = 1 / 2 , 1 , 3 / 2 . From (10), we need to set the weight w * :
w * = w 0 w 1 1 4 λ 2 α 2 = 3 2 , 6 5 , 3 2 5 .
Choose w 0 = 1 + α 3 . Then, w 1 = 3 1 + 3 2 2 + 3 4 2 + 6 α , 3 2 1 5 2 3 5 + 3 2 + 3 2 5 α , 3 1 3 3 2 10 + 3 3 + 3 2 10 α . The corresponding rational quadratic curves p ( s ) are shown in Figure 1a (blue curves).
Example 2.
( α 2 = 1 ) Consider two endpoints p 0 = ( 1 , 2 ) , p 1 = ( 5 , 3 ) and their corresponding weights w 0 = 3 , w 1 = 9 . Choose p * = ( p 0 + p 1 ) / 2 + α ( p 1 p 0 ) λ with λ = 3 / 16 , 5 / 16 , 7 / 16 . From (10), we need to set the weight w * :
w * = w 0 w 1 1 4 α 2 λ 2 = 8 3 7 , 6 , 24 5 .
Let w 0 = ( 2 , 1 ) . Then, w 1 = 17 21 + 10 21 α , 5 + 4 α , 23 5 + 22 5 α . The corresponding rational quadratic curves p ( s ) are shown in Figure 1b (blue curves).
Example 3.
( α 2 = 0 ) Consider two endpoints p 0 = ( 1 , 2 ) , p 1 = ( 5 , 3 ) and their corresponding weights w 0 = 4 , w 1 = 9 . Choose p * = ( p 0 + p 1 ) / 2 + α ( p 1 p 0 ) λ with λ = 1 / 3 , 1 , 5 / 3 . From (10), we need to set the weight w * :
w * = w 0 w 1 1 4 α 2 λ 2 = 6 .
Let w 0 = ( 2 , 1 ) . Then, w 1 = 3 + 7 2 α , 3 + 15 2 α , 3 + 23 2 α . The corresponding rational quadratic curves p ( s ) are shown in Figure 1c (blue curves).

4. Construction of Quadratic Birational Maps via Generalized Complex Representations

To generate planar rational maps, we generalize p l ( s ) in (6) to the bivariate setting as follows.
Consider three non-collinear endpoints p 100 , p 010 , p 001 . Let
p l ( s , t ) = p 100 w 100 u + p 010 w 010 s + p 001 w 001 t w 100 u + w 010 s + w 001 t , 0 u , s , t 1 , u = 1 s t ,
where w 100 , w 010 , w 001 are generalized complex weights, denoted by
w 100 = u 0 + α v 0 , w 010 = u 1 + α v 1 , w 001 = u 2 + α v 2 , p 100 w 100 = x ^ 0 + α y ^ 0 , p 010 w 010 = x ^ 1 + α y ^ 1 , p 001 w 001 = x ^ 2 + α y ^ 2 .
To avoid the degenerate case, we also assume w 100 , w 010 , w 001 0 . In addition, at least one of w 100 , w 010 , w 001 is not fully real throughout the rest of the paper. Note that the complex, hyperbolic and parabolic rational curves are invariant with respect to Möbius transformations, i.e., applying a Möbius transformation to the control points and weights is equivalent to applying the Möbius transformation to every point on the curve (for more details, the readers are referred to Theorem 4.5 in [21]).
Multiplying the numerator and denominator in (13) by the conjugate of the denominator yields the corresponding quadratic real rational planar map:
p ( s , t ) = i + j + k = 2 p i j k w i j k B i j k 2 ( u , s , t ) i + j + k = 2 w i j k B i j k 2 ( u , s , t ) ,
where B i j k n ( u , s , t ) are Bernstein polynomials defined by
B i j k n ( u , s , t ) = n ! i ! j ! k ! u i s j t k , i + j + k = n ,
and w i j k = ( w i 1 j 1 k 1 w i 2 j 2 k 2 ¯ + w i 1 j 1 k 1 ¯ w i 2 j 2 k 2 ) / 2 ,
p i j k = ( w i 1 j 1 k 1 w i 2 j 2 k 2 ¯ p i 1 j 1 k 1 + w i 1 j 1 k 1 ¯ w i 2 j 2 k 2 p i 2 j 2 k 2 ) / ( 2 w i j k ) ,
( i , j , k ) = ( i 1 + i 2 , j 1 + j 2 , k 1 + k 2 ) , i 1 + j 1 + k 1 = i 2 + j 2 + k 2 = 1 .
Denote by P ( s , t ) = a ( s , t ) , b ( s , t ) , c ( s , t ) , the homogeneous form of p ( s , t ) = ( a ( s , t ) c ( s , t ) , b ( s , t ) c ( s , t ) ) . Next, we shall derive the syzygies of a planar rational map P ( s , t ) that is endowed with a generalized complex representation (13).
From Lemma 1, there are two special syzygies L 0 ( s , t ) , L 1 ( s , t ) in the form (5) of P ( s , t ) , which are R [ s , t ] -independent. Let l 0 ( x , y ; s , t ) = L 0 ( s , t ) · ( x , y , 1 ) and l 1 ( x , y ; s , t ) = L 0 ( s , t ) · ( x , y , 1 ) . Then, after calculation, we obtain
l 0 ( x , y ; s , t ) l 1 ( x , y ; s , t ) = m s s + m t t + m 0 ,
where
m s = x ( u 1 u 0 ) + α 2 y ( v 1 v 0 ) + ( x ^ 1 x ^ 0 ) x ( v 1 v 0 ) + y ( u 1 u 0 ) + ( y ^ 1 y ^ 0 ) , m t = x ( u 2 u 0 ) + α 2 y ( v 2 v 0 ) + ( x ^ 2 x ^ 0 ) x ( v 2 v 0 ) + y ( v 2 v 0 ) + ( y ^ 2 y ^ 0 ) , m 0 = x u 0 + α 2 y v 0 x ^ 0 x v 0 + y u 0 y ^ 0 .
With the above notation in hand, we are now ready to give the inversion formulas for a rational quadratic map with a generalized complex rational representation.
Theorem 2.
For a real planar rational quadratic map p ( s , t ) with a generalized complex rational representation (13), the inversion equation is
s = det ( M 1 ) / det ( M 0 ) , t = det ( M 2 ) / det ( M 0 ) ,
where
M 0 = ( m s , m t ) , M 1 = ( m 0 , m t ) , M 2 = ( m s , m 0 ) .
Hence, the rational map (13) with generalized complex rational representation is birational.
Proof. 
The proof is similar to the proof of Theorem 4 in [17].    □
Given three endpoints p 200 , p 020 , p 002 and their corresponding weights w 200 , w 020 , w 002 > 0 , next we shall show how to obtain p 110 , p 101 , p 011 and w 110 , w 101 , w 011 such that the rational quadratic map is endowed with a generalized complex rational representation.
Note that not all three middle points are free to be chosen arbitrarily. Without loss of generality, we consider the two boundary curves p ( s , 0 ) , p ( 0 , t ) first. Choose two parameters λ 110 , λ 101 R to fix the middle points p 110 , p 101 , i.e.,
p 110 = ( p 200 + p 020 ) / 2 + α ( p 200 p 020 ) λ 110 , p 101 = ( p 200 + p 002 ) / 2 + α ( p 200 p 002 ) λ 101 .
Then from (10) in Theorem 1, we obtain
w 110 = w 200 w 020 1 4 α 2 λ 110 2 , w 101 = w 200 w 002 1 4 α 2 λ 101 2 .
By setting the generalized complex weight w 100 such that w 200 = w 100 ¯ w 100 , we can compute w 010 , w 001 from (11).
Let p 100 = p 200 , p 010 = p 020 , p 001 = p 002 , then
w 011 = ( w 010 w 001 ¯ + w 010 w 001 ¯ ) / 2 , p 011 = ( w 010 w 001 ¯ p 010 + w 001 w 010 ¯ p 001 ) / ( 2 w 011 ) .
In this way, we can construct a planar rational quadratic map (14) that is endowed with a generalized complex rational representation
p l ( s , t ) = p 100 w 100 u + p 010 w 010 s + p 001 w 001 t w 100 u + w 010 s + w 001 t .
Based on Theorem 2, the planar rational quadratic map p ( s , t ) is also birational.
In general, the method of constructing a planar rational quadratic map that is endowed with a generalized complex rational representation can be summarized in the following Algorithm 1.
Algorithm 1 To construct a birational quadratic planar map with a generalized complex rational representation
Input: Three endpoints p 200 , p 020 , p 002 ; weights w 200 , w 020 , w 002 > 0 ; two parameters λ 110 , λ 101 R
Output: A quadratic planar map p ( s , t ) that is endowed with generalized complex representation
  • Calculate two real weights w 110 = w 200 w 020 1 4 α 2 λ 110 2 , w 101 = w 200 w 002 1 4 α 2 λ 101 2
  • Choose a generalized complex weight w 100 such that w 100 ¯ w 100 = w 200 . Calculate the other two weights w 010 , w 001 as (11)
  • Calculate the real weight w 011 = ( w 010 w 001 ¯ + w 010 w 001 ¯ ) / 2
  • Set three middle points p 110 = ( p 200 + p 020 ) / 2 + α ( p 200 p 020 ) λ 110 , p 101 = ( p 200 + p 002 ) / 2 + α ( p 200 p 002 ) λ 101 , p 011 = ( w 010 w 001 ¯ p 010 + w 001 w 010 ¯ p 001 ) / ( 2 w 011 )
return p ( s , t ) = i + j + k = 2 p i j k w i j k B i j k 2 ( u , s , t ) i + j + k = 2 w i j k B i j k 2 ( u , s , t ) .
Example 4.
( α 2 = 1 ) Fix three endpoints p 200 = ( 0 , 0 ) , p 020 = ( 6 , 2 ) , p 002 = ( 3 , 4 ) . Set the weights w 200 = 36 , w 020 = 35 , w 002 = 8 . The control points p 110 , p 101 are fixed by choosing λ 110 = 1 / 12 , λ 101 = 1 / 6 , i.e.,
p 110 = ( p 200 + p 020 ) / 2 + α ( p 200 p 020 ) λ 110 = 19 / 6 + 3 α / 2 , p 101 = ( p 200 + p 002 ) / 2 + α ( p 200 p 002 ) λ 101 = 1 / 6 + 5 α / 2 ) .
Hence
w 110 = w 200 w 020 1 4 α 2 λ 110 2 = 36 , w 101 = w 200 w 002 1 4 α 2 λ 101 2 = 18 .
From (12), we obtain
z 010 = w 110 + α ( w 110 ) 2 w 200 w 020 α 2 / w 200 = 1 + α / 6 , z 001 = w 101 + α ( w 101 ) 2 w 200 w 002 α 2 / w 200 = 1 / 2 + α / 6 .
Choose w 100 = 6 + 0 α . Hence, from (11),
w 010 = w 100 z 010 = 6 + α , w 001 = w 100 z 001 = 3 + α .
Set p 100 = p 200 , p 010 = p 020 , p 001 = p 002 . w 011 and p 011 can be calculated as follows:
w 011 = ( w 010 w 001 ¯ + w 010 w 001 ¯ ) / 2 , = 17 , p 011 = ( w 010 w 001 ¯ p 010 + w 001 w 010 ¯ p 001 ) / ( 2 w 011 ) = 159 / 34 + 93 α / 34 .
Thus, we obtain a rational quadratic map p ( s , t ) with complex rational representation
p l ( s , t ) = p 100 w 100 u + p 010 w 010 s + p 001 w 001 t w 100 u + w 010 s + w 001 t = ( 38 s + 13 t ) + ( 18 s + 15 t ) α ( 6 3 t ) + ( s + t ) α .
See Figure 2 for an illustration of the birational map p ( s , t ) .
Let n 0 ( s , t ) = 38 s + 13 t , n 1 ( s , t ) = 18 s + 15 t , d 0 ( s , t ) = 6 3 t , d 1 ( s , t ) = s + t . Solving the system of equations
l 0 ( x , y ; s , t ) = d 0 ( s , t ) x + α 2 d 1 ( s , t ) y n 0 ( s , t ) = 0 , l 1 ( x , y ; s , t ) = d 1 ( s , t ) x + d 0 ( s , t ) y n 1 ( s , t ) = 0 ,
with respect to the variables s , t , we obtain the inversion equation:
s = 6 x 2 15 x y 2 + 13 y 3 x 2 79 x 3 y 2 + 117 y + 336 , t = 6 x 2 18 x y 2 + 38 y 3 x 2 79 x 3 y 2 + 117 y + 336 .
Example 5.
( α 2 = 0 ) Fix three endpoints p 200 = ( 0 , 0 ) , p 020 = ( 6 , 2 ) , p 002 = ( 3 , 4 ) . Set the weights w 200 = 9 , w 020 = 36 , w 002 = 1 . The control points p 110 , p 101 are fixed by choosing λ 110 = 1 / 6 , λ 101 = 1 / 6 , i.e.,
p 110 = ( p 200 + p 020 ) / 2 + α ( p 200 p 020 ) λ 110 = 3 + 0 α , p 101 = ( p 200 + p 002 ) / 2 + α ( p 200 p 002 ) λ 101 = 3 / 2 + 5 α / 2 .
Hence
w 110 = w 200 w 020 1 4 α 2 λ 110 2 = 18 , w 101 = w 200 w 002 1 4 α 2 λ 101 2 = 3 .
From (12), we obtain
z 010 = w 110 1 + 2 α λ 110 / w 200 = 2 2 α / 3 , z 001 = w 101 1 + 2 α λ 101 / w 200 = 1 / 3 + α / 9 .
Choose w 100 = 6 + 0 α . Hence, from (11),
w 010 = w 100 z 010 = 6 + 2 α , w 001 = w 100 z 001 = 1 + α .
Set p 100 = p 200 , p 010 = p 020 , p 001 = p 002 . w 011 and p 011 can be calculated as follows:
w 011 = ( w 010 w 001 ¯ + w 010 w 001 ¯ ) / 2 = 6 , p 011 = ( w 010 w 001 ¯ p 010 + w 001 w 010 ¯ p 001 ) / ( 2 w 011 ) = 9 / 2 + 2 α .
Thus, we obtain a rational quadratic map p ( s , t ) with complex rational representation
p l ( s , t ) = p 100 w 100 u + p 010 w 010 s + p 001 w 001 t w 100 u + w 010 s + w 001 t = ( 36 s + 3 t ) + ( 24 s + 7 t ) α ( 3 + 3 s 2 t ) + ( 2 t ) α .
See Figure 3 for an illustration of the birational map p ( s , t ) .
Let n 0 ( s , t ) = 36 s + 3 t , n 1 ( s , t ) = 24 s + 7 t , d 0 ( s , t ) = 3 + 3 s 2 t , d 1 ( s , t ) = 2 t , solving the system of equations
l 0 ( x , y ; s , t ) = d 0 ( s , t ) x + α 2 d 1 ( s , t ) y n 0 ( s , t ) = 0 , l 1 ( x , y ; s , t ) = d 1 ( s , t ) x + d 0 ( s , t ) y n 1 ( s , t ) = 0 ,
and with respect to the variables s , t , we obtain the inversion equation:
s = x 2 + 15 x 9 y 3 x 2 + 11 x 27 y 60 , t = 2 x 2 18 y x 2 + 11 x 27 y 60 .

5. Conclusions

We extended the construction of birational quadratic maps based on rational maps over generalized complex numbers. When we need to construct birational quadratic maps, the new method provides more options for choosing the middle control points and their corresponding weights. Hence, we increase the flexibility of the boundary curves of birational quadratic maps. The proof of the birational property is based on the construction of two special syzygies. The formulas for the inversion equation are also provided. Although we have generalized the construction of birational maps based on generalized complex representations, the choice of the middle control points is still limited, and we cannot construct quadratic birational planar maps by moving the three middle points freely. In the future, we hope to construct general birational maps by setting α 2 to the real values other than 1 , 0 , 1 , which will be used as a shape parameter to control the shape of the boundaries.

Author Contributions

Conceptualization, X.W.; Methodology, R.G.; Validation, Y.H., Q.N. and R.L.; Writing—original draft, X.W.; Writing—review & editing, R.G.; Funding acquisition, X.W. and Q.N. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China No. 122011292, No. 61772167, and No. 11771420; Natural Science Foundation of the Jiangsu Higher Education Institutions of China: 22KJB110015.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hudson, H.P. Cremona Transformations in Plane and Space; The University Press: Cambridge, UK, 1927; Volume 38. [Google Scholar]
  2. Sederberg, T.W.; Parry, S.R. Free-form deformation of solid geometric models. In Proceedings of the 13th Annual Conference on Computer Graphics and Interactive Techniques, Dallas, TX, USA, 18–22 August 1986; pp. 151–160. [Google Scholar]
  3. Lee, S.Y.; Chwa, K.Y.; Shin, S.Y. Image metamorphosis using snakes and free-form deformations. In Proceedings of the 22nd Annual Conference on Computer Graphics and Interactive Techniques, Los Angeles, CA, USA, 6–11 August 1995; pp. 439–448. [Google Scholar]
  4. Wolberg, G. Image morphing: A survey. Vis. Comput. 1998, 14, 360–372. [Google Scholar] [CrossRef]
  5. Parke, F.I.; Waters, K. Computer Facial Animation; CRC Press: Boca Raton, FL, USA, 2008. [Google Scholar]
  6. Rohlfing, T.; Maurer, C.R.; Bluemke, D.A.; Jacobs, M.A. Volume-preserving nonrigid registration of MR breast images using free-form deformation with an incompressibility constraint. IEEE Trans. Med. Imaging 2003, 22, 730–741. [Google Scholar] [CrossRef] [PubMed]
  7. Coquillart, S. Extended free-form deformation: A sculpturing tool for 3d geometric modeling. In Proceedings of the 17th Annual Conference on Computer Graphics and Interactive Techniques, Dallas, TX, USA, 6–10 August 1990; pp. 187–196. [Google Scholar]
  8. Gain, J.E.; Dodgson, N.A. Preventing self-intersection under free-form deformation. IEEE Trans. Vis. Comput. Graph. 2001, 7, 289–298. [Google Scholar] [CrossRef]
  9. Sederberg, T.W.; Goldman, R.N.; Wang, X. Birational 2D free-form deformation of degree 1× n. Comput. Aided Geom. Des. 2016, 44, 1–9. [Google Scholar] [CrossRef]
  10. Sederberg, T.W.; Zheng, J. Birational quadrilateral maps. Comput. Aided Geom. Des. 2015, 32, 1–4. [Google Scholar] [CrossRef]
  11. Deng, C.; Zhu, F.; Liu, J. The limit of a family of barycentric coordinates for quadrilaterals. Comput. Aided Geom. Des. 2015, 38, 38–39. [Google Scholar] [CrossRef]
  12. Floater, M.S. The inverse of a rational bilinear mapping. Comput. Aided Geom. Des. 2015, 33, 46–50. [Google Scholar] [CrossRef]
  13. Deng, C.; Li, Y.; Mu, X.; Zhao, Y. Characteristic conic of rational bilinear map. J. Comput. Appl. Math. 2019, 346, 277–283. [Google Scholar] [CrossRef]
  14. Botbol, N.; Busé, L.; Chardin, M.; Hassanzadeh, S.H.; Simis, A.; Tran, Q.H. Effective criteria for bigraded birational maps. J. Symb. Comput. 2017, 81, 69–87. [Google Scholar] [CrossRef]
  15. Sánchez-Reyes, J. Complex rational Bézier curves. Comput. Aided Geom. Des. 2009, 26, 865–876. [Google Scholar] [CrossRef]
  16. Wang, X.; Goldman, R. μ-bases for complex rational curves. Comput. Aided Geom. Des. 2013, 30, 623–635. [Google Scholar] [CrossRef]
  17. Wang, X.; Wu, M.; Liu, Y.; Ni, Q. Constructing quadratic birational maps via their complex rational representation. Comput. Aided Geom. Des. 2021, 85, 101969. [Google Scholar] [CrossRef]
  18. Kurzeja, K.; Rossignac, J. Becots: Bent corner-operated tran-similar maps and lattices. Comput.-Aided Des. 2020, 129, 102912. [Google Scholar] [CrossRef]
  19. Rossignac, J. Corner-operated tran-similar (cots) maps, patterns, and lattices. ACM Trans. Graph. 2020, 39, 1–14. [Google Scholar] [CrossRef]
  20. Tilove. Set membership classification: A unified approach to geometric intersection problems. IEEE Trans. Comput. 1980, C-29, 874–883. [Google Scholar] [CrossRef]
  21. Du, J.; Goldman, R.; Wang, X. Rational curves over generalized complex numbers. J. Symb. Comput. 2019, 93, 56–84. [Google Scholar] [CrossRef]
  22. Kulyabov, D.; Korolkova, A.; Gevorkyan, M. Hyperbolic numbers as einstein numbers. J. Phys. Conf. Ser. 2020, 1557, 012027. [Google Scholar] [CrossRef]
  23. Salomon, D. Computer Graphics and Geometric Modeling; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2012. [Google Scholar]
  24. VanArsdale, D. Homogeneous transformation matrices for computer graphics. Comput. Graph. 1994, 18, 177–191. [Google Scholar] [CrossRef]
  25. Sobczyk, G. The hyperbolic number plane. Coll. Math. J. 1995, 26, 268–280. [Google Scholar] [CrossRef]
  26. Fischer, I. Dual-Number Methods in Kinematics, Statics and Dynamics; CRC Press: Boca Raton, FL, USA, 1998. [Google Scholar]
  27. Chen, F.; Cox, D.; Liu, Y. The μ-basis and implicitization of a rational parametric surface. J. Symb. Comput. 2005, 39, 689–706. [Google Scholar] [CrossRef]
  28. Farin, G.E.; Farin, G. Curves and Surfaces for CAGD: A Practical Guide; Morgan Kaufmann: Burlington, MA, USA, 2002. [Google Scholar]
Figure 1. Generalized complex rational curves of degree one. The control polygon (red color) and the corresponding curve (blue color) vary as the parameter λ changes.
Figure 1. Generalized complex rational curves of degree one. The control polygon (red color) and the corresponding curve (blue color) vary as the parameter λ changes.
Mathematics 11 03609 g001
Figure 2. A birational quadratic planar map ( α 2 = 1 ). The points p 011 on the red line segments are fixed when the points p 101 , p 110 on the blue line segments are chosen.
Figure 2. A birational quadratic planar map ( α 2 = 1 ). The points p 011 on the red line segments are fixed when the points p 101 , p 110 on the blue line segments are chosen.
Mathematics 11 03609 g002
Figure 3. A birational quadratic planar map ( α 2 = 0 ). The points p 011 on the red line segments are fixed when the points p 101 , p 110 on the blue line segments are chosen.
Figure 3. A birational quadratic planar map ( α 2 = 0 ). The points p 011 on the red line segments are fixed when the points p 101 , p 110 on the blue line segments are chosen.
Mathematics 11 03609 g003
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wang, X.; Han, Y.; Ni, Q.; Li, R.; Goldman, R. Birational Quadratic Planar Maps with Generalized Complex Rational Representations. Mathematics 2023, 11, 3609. https://doi.org/10.3390/math11163609

AMA Style

Wang X, Han Y, Ni Q, Li R, Goldman R. Birational Quadratic Planar Maps with Generalized Complex Rational Representations. Mathematics. 2023; 11(16):3609. https://doi.org/10.3390/math11163609

Chicago/Turabian Style

Wang, Xuhui, Yuhao Han, Qian Ni, Rui Li, and Ron Goldman. 2023. "Birational Quadratic Planar Maps with Generalized Complex Rational Representations" Mathematics 11, no. 16: 3609. https://doi.org/10.3390/math11163609

APA Style

Wang, X., Han, Y., Ni, Q., Li, R., & Goldman, R. (2023). Birational Quadratic Planar Maps with Generalized Complex Rational Representations. Mathematics, 11(16), 3609. https://doi.org/10.3390/math11163609

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop