# Drawable Region of the Generalized Log Aesthetic Curves

Authors: R. U. Gobithaasan; R. Karpagavalli; Kenjiro T. Miura
Journal: Journal of Applied Mathematics
 (2013)

---

## Abstract

The main characteristic of visually pleasing curves used for product design is a monotonic curvature profile. Recently, a planar curve called Generalized Log Aesthetic Curve (GLAC) has been extended from the Log Aesthetic Curve (LAC), and it has an additional shape parameter,ν. This curve preserves the monotonicity of curvature and is said to produce visually pleasing curves. This paper delves on the drawable region of the GLAC segment which indicates the probable solutions of shape parameters from given interpolating points and the direction of travel at those points. The first section reviews the formulation of GLAC and its related bounds. The section describes the algorithm for identifying the drawable region. It is followed by the section describing how small changes of ν widen the drawable boundaries. The final section discusses the superiority of GLAC compared to LAC for use in industrial product design.


---

## Body

## 1. Introduction

The importance of aesthetic shapes to design industrial products leads to the studies of planar aesthetic curves [1]. High quality curves are also known as fair curves, minimal energy curves, aesthetic curves, and so forth. The main characteristic of these curves is that it has a monotonic curvature profile. The mathematical term used to describe planar curves with monotonic (either increase or decrease) curvature is called a spiral [2]. A curve is categorized as fair if its curvature plot consists of relatively minimum number of monotone pieces.In 1999, Harada et al. presented a quantitative method of investigating curves used in automobile design which is called Logarithmic Distribution Diagram of Curvature (LDDC). They defined aesthetic curves as curves with a constant LDDC gradient (denoted asα). Consequently in 2003, Kanaya et al. simplified the formulation of LDDC to Logarithmic Curvature Graph (LCG). In 2005, Miura analyzed the features of aesthetic curves and derived a general formula of aesthetic curves called Log Aesthetic Curve (LAC). Yoshida and Saito [3] further investigated LAC to identify the overall shapes of LAC. They proposed a method to draw LAC curve segment interactively by using two endpoints and their respective tangent vectors. In 2009, Levien and Sequin stated that the LA curve is the most promising curve for aesthetic design. In 2012, Yoshida and Saito derived a method to render the drawable boundary for LAC segments.The number of research papers produced since the introduction of LAC has been increasing exponentially such as the reformulation of 2D and 3D LAC in the form of variational principal [4], analytic representation of LA curve using Incomplete Gamma function for efficient computation [5], and the formulation of G2 LA spline for automobile design application [6].Recently, a planar curve called Generalized Log Aesthetic Curve (GLAC) has been developed by extending the formulation of Generalized Cornu Spiral (GCS) [7] in a similar manner to the LAC [8]. The family of GLAC comprises planar curves of high quality such as GCS, LAC, clothoid, Nielsen's spiral, logarithmic spiral, circle involute, and so forth. The GLAC segment has an additional parameter (denoted as ν) to determine its shape as compared to GCS and LAC segment. Hence, an extra constraint can be satisfied when shaping the GLAC segment. Recent advancements of GLAC include the extension of spatial GLAC [9] and the shape analysis of GLAC [10].Interactive control of GLAC segment has been developed similar to interactive LAC segment [11]. However, it is necessary to analyze the drawable region of interactive GLAC segment so that the users are aware of the flexibility of GLAC for shape design. The identification of drawable region of GLAC is similar to the method employed for LAC [12]. However, a slight modification has been carried out to determine the shape parameter Λ. We modify the arc length instead of modifying the θd directly. We will modify arc length in the range of 0<s<∞, which will give θd such that 0<θd<π.


### 1.1. The Formulation of GLAC

GLAC is derived via the curve synthesis process where a formulation of a curve is derived from defined curvature function. The advantage of GLAC is that it has an extra shape parameter (ν) which can be used to dictate the curvature radius (ρ) at the starting point of the GLAC. GLAC becomes LAC when ν becomes 0. The LCG and its gradient function, λ(s) are derived for GLAC as follows [8]:
(1)LCGGLAC(s)={1(Λαs+1)-1α+νlog[1(Λαs+1)-1/α+ν],log[1Λ(Λαs+1)(1+ν(Λαs+1)1/α)][1(Λαs+1)-1α+ν]},λ(s)GLAC=α+ν(Λαs+1)1/α(1+α).Arc length and tangential angle function of GLAC are shown below, respectively:(2)sGLAC(ρ)={1Λlog[1ρ-1-ν]ifα=0,1Λα((ρ-1-ν)-α-1)otherwise,θGLAC(s)={1Λ(1-e-Λs)+νsifα=0,1Λlog[Λs+1]+νsifα=1,1Λ(α-1)((Λαs+1)(α-1)/α-1)+νsotherwise.The parametric form of GLAC is(3)CGLAC(s)={∫0scos[θGLAC(u)]du,∫0ssin[θGLAC(u)]du},
where {Λ,α,ν}∈ℝ are the variables that can be used to shape the GLAC segment.


### 1.2. The Bounds of GLAC

The bounds are derived for both arc length and tangential angle function of GLACwhich are depicted in Table1 and Table 2 respectively. For the case of ν=0, the bounds of GLAC are similar to the bounds of LAC as stated in [10].Table 1
The lower bound and upper bound fors.
Arc length
ν
<
0
(
ν
≠
-
1
)
ν
>
0
ν
=
-
1
LB
UB
LB
UB
LB
UB
α
<
0
—
(
-
ν
)
-
α
-
1
Λ
α
—
-
1
Λ
α
—
0
α
=
0
—
-
log
[
-
ν
]
Λ
—
—
—
0
α
>
0
-
1
Λ
α
(
-
ν
)
-
α
-
1
Λ
α
-
1
Λ
α
—
-
1
Λ
α
0Table 2
The lower bound and upper bound forθ.
Tangential angle
ν
<
0 (ν≠-1)
ν
>
0
ν
=
-
1
LB
UB
LB
UB
LB
UB
α
<
0
—
A
—
B
—
0
α
=
0
—
1
Λ
(
1
+
ν
-
ν
log
[
-
ν
]
)
—
—
—
0
0
<
α
<
1
—
A
—
—
—
0
α
=
1
—
-
1
Λ
(
1
+
ν
+
log
[
-
ν
]
)
—
—
—
0
α
>
1
B
A
B
—
1
Λ
α
(
1
-
α
)
0
A and B represent1Λ(α-1)((-ν)1-α-1)+νΛα((-ν)-α-1) and1Λ(11-α-να).


## 1.1. The Formulation of GLAC

GLAC is derived via the curve synthesis process where a formulation of a curve is derived from defined curvature function. The advantage of GLAC is that it has an extra shape parameter (ν) which can be used to dictate the curvature radius (ρ) at the starting point of the GLAC. GLAC becomes LAC when ν becomes 0. The LCG and its gradient function, λ(s) are derived for GLAC as follows [8]:
(1)LCGGLAC(s)={1(Λαs+1)-1α+νlog[1(Λαs+1)-1/α+ν],log[1Λ(Λαs+1)(1+ν(Λαs+1)1/α)][1(Λαs+1)-1α+ν]},λ(s)GLAC=α+ν(Λαs+1)1/α(1+α).Arc length and tangential angle function of GLAC are shown below, respectively:(2)sGLAC(ρ)={1Λlog[1ρ-1-ν]ifα=0,1Λα((ρ-1-ν)-α-1)otherwise,θGLAC(s)={1Λ(1-e-Λs)+νsifα=0,1Λlog[Λs+1]+νsifα=1,1Λ(α-1)((Λαs+1)(α-1)/α-1)+νsotherwise.The parametric form of GLAC is(3)CGLAC(s)={∫0scos[θGLAC(u)]du,∫0ssin[θGLAC(u)]du},
where {Λ,α,ν}∈ℝ are the variables that can be used to shape the GLAC segment.


## 1.2. The Bounds of GLAC

The bounds are derived for both arc length and tangential angle function of GLACwhich are depicted in Table1 and Table 2 respectively. For the case of ν=0, the bounds of GLAC are similar to the bounds of LAC as stated in [10].Table 1
The lower bound and upper bound fors.
Arc length
ν
<
0
(
ν
≠
-
1
)
ν
>
0
ν
=
-
1
LB
UB
LB
UB
LB
UB
α
<
0
—
(
-
ν
)
-
α
-
1
Λ
α
—
-
1
Λ
α
—
0
α
=
0
—
-
log
[
-
ν
]
Λ
—
—
—
0
α
>
0
-
1
Λ
α
(
-
ν
)
-
α
-
1
Λ
α
-
1
Λ
α
—
-
1
Λ
α
0Table 2
The lower bound and upper bound forθ.
Tangential angle
ν
<
0 (ν≠-1)
ν
>
0
ν
=
-
1
LB
UB
LB
UB
LB
UB
α
<
0
—
A
—
B
—
0
α
=
0
—
1
Λ
(
1
+
ν
-
ν
log
[
-
ν
]
)
—
—
—
0
0
<
α
<
1
—
A
—
—
—
0
α
=
1
—
-
1
Λ
(
1
+
ν
+
log
[
-
ν
]
)
—
—
—
0
α
>
1
B
A
B
—
1
Λ
α
(
1
-
α
)
0
A and B represent1Λ(α-1)((-ν)1-α-1)+νΛα((-ν)-α-1) and1Λ(11-α-να).


## 2.  Drawable Region of GLAC

Interactive control of GLAC segment can be developed either by stating the value ofν or allowing the optimization techique to find a suitable value for ν. It is possible to construct the drawable region for interactive GLAC if the designer chooses to set the value for ν. The analysis of the drawable region can then be used to identify all possible solutions that exist for a given input.


### 2.1.  The Method

The method used is similar to LAC where three control points will be placed in a rectangle of corners(±1,±1). Values of α<0 are used to construct the drawable region for the case of ν>0. For the case of -1<ν<0, there is no restriction for α. Figure 1 illustrates the setup of control points in order to obtain the drawable region of GLAC.A drawable region of interactive GLAC segment.
(a)
GLAC segment with new endpoints
(b)
GLAC segment and a drawable boundaryLet pointsP0, P1, and P2 be the control points for drawing a GLAC segment. The points P0 and P2 are placed at (-1,0) and (1,0). The point P1 is the point on the drawable boundary in the rectangle. To note, P1 is unknown and will be determined by the tangential angle θd on the curve. GLAC segment with tangential angle from 0 to θd is used for all α. Point P1 can be found on the rectangle by performing transformation. By changing the arc length from 0 to a large number, the points of P1 will form the drawable boundary. Note that we are modifying arc length to obtain θd in the range of 0<θd<π. Hence, by applying reflection on the x-axis and y-axis, a complete drawable boundary can be found.


### 2.2. The Algorithm

The objective of this algorithm (Algorithm1) is to determine the second control point through the formulation of GLAC curve segment and by plotting it on the rectangle.Algorithm 1
Remark: Control points{P0,P2} are given with a preferred α. Via GLAC function, second
control point is determined. New points are transformed to the originals and plotted.s is modified
from 0 to a large number such thatθd is in range of 0<θd<π to get the complete drawable boundary.
Input: P0(-1,0),P2(1,0),α,ν,s.
Output: Complete drawable boundary of interactive GLAC.
Begin
Step 1. If -1<ν<0 then
ifα=0 then
Λ←-log[-ν]s.
else
Λ←(-ν)-α-1sα.
else ifν>0 and α<0 then
Λ←-1sα.
Step 2. Determine SP^(x0,y0) and EP^(x2,y2) using (3).
Step 3. Determine scaling factor r←∥P2P0∥∥SP^EP^∥ and angle θe←cos-1[x2x22+y22].
Step 4. Determine P1^ using the tangent line of both end points.
Step 5. Scale points to r, reflect through x-axis, rotate to -θe and translate to {P0,P2}.
Step 6. Plot point P1.
Step 7. Modify s from 0 to a large number and repeat Step 1. to 6.
Step 8. Reflect points via x-axis andy-axis.
Step 9. Output.
End.


## 2.1.  The Method

The method used is similar to LAC where three control points will be placed in a rectangle of corners(±1,±1). Values of α<0 are used to construct the drawable region for the case of ν>0. For the case of -1<ν<0, there is no restriction for α. Figure 1 illustrates the setup of control points in order to obtain the drawable region of GLAC.A drawable region of interactive GLAC segment.
(a)
GLAC segment with new endpoints
(b)
GLAC segment and a drawable boundaryLet pointsP0, P1, and P2 be the control points for drawing a GLAC segment. The points P0 and P2 are placed at (-1,0) and (1,0). The point P1 is the point on the drawable boundary in the rectangle. To note, P1 is unknown and will be determined by the tangential angle θd on the curve. GLAC segment with tangential angle from 0 to θd is used for all α. Point P1 can be found on the rectangle by performing transformation. By changing the arc length from 0 to a large number, the points of P1 will form the drawable boundary. Note that we are modifying arc length to obtain θd in the range of 0<θd<π. Hence, by applying reflection on the x-axis and y-axis, a complete drawable boundary can be found.


## 2.2. The Algorithm

The objective of this algorithm (Algorithm1) is to determine the second control point through the formulation of GLAC curve segment and by plotting it on the rectangle.Algorithm 1
Remark: Control points{P0,P2} are given with a preferred α. Via GLAC function, second
control point is determined. New points are transformed to the originals and plotted.s is modified
from 0 to a large number such thatθd is in range of 0<θd<π to get the complete drawable boundary.
Input: P0(-1,0),P2(1,0),α,ν,s.
Output: Complete drawable boundary of interactive GLAC.
Begin
Step 1. If -1<ν<0 then
ifα=0 then
Λ←-log[-ν]s.
else
Λ←(-ν)-α-1sα.
else ifν>0 and α<0 then
Λ←-1sα.
Step 2. Determine SP^(x0,y0) and EP^(x2,y2) using (3).
Step 3. Determine scaling factor r←∥P2P0∥∥SP^EP^∥ and angle θe←cos-1[x2x22+y22].
Step 4. Determine P1^ using the tangent line of both end points.
Step 5. Scale points to r, reflect through x-axis, rotate to -θe and translate to {P0,P2}.
Step 6. Plot point P1.
Step 7. Modify s from 0 to a large number and repeat Step 1. to 6.
Step 8. Reflect points via x-axis andy-axis.
Step 9. Output.
End.


## 3.  Analysis Result

Figures2, 3, 4, 5, 6, 7, 8, 9, 10, and 11 show the obtained result for the drawable region of the interactive GLAC segment depending on α andν. Figures 2, 3, 4, 5, 6, and 7 show the pattern of how the decrease and increase ofν affects the drawable region. Figures 8, 9, 10, and 11 show the drawable region for α=-1,0,1,2 as ν increases. The white region implies the drawable region of GLAC segment, while the gray region indicates no solution. Hence, more probable GLAC shapes can be drawn when the white region widens.The drawable region of GLAC whenν=-0.3.
(a)
α
=
-
3
(b)
α
=
-
2
(c)
α
=
-
1
(d)
α
=
-
0.5
(e)
α
=
-
0.1
(f)
α
=
0
(g)
α
=
0.5
(h)
α
=
0.8
(i)
α
=
1
(j)
α
=
1.1
(k)
α
=
2
(l)
α
=
3The drawable region of GLAC whenν=-0.1.
(a)
α
=
-
3
(b)
α
=
-
2
(c)
α
=
-
1
(d)
α
=
-
0.5
(e)
α
=
-
0.1
(f)
α
=
0
(g)
α
=
0.5
(h)
α
=
0.8
(i)
α
=
1
(j)
α
=
1.1
(k)
α
=
2
(l)
α
=
3The drawable region of GLAC whenν=-0.01.
(a)
α
=
-
3
(b)
α
=
-
2
(c)
α
=
-
1
(d)
α
=
-
0.5
(e)
α
=
-
0.1
(f)
α
=
0
(g)
α
=
0.5
(h)
α
=
0.8
(i)
α
=
1
(j)
α
=
1.1
(k)
α
=
2
(l)
α
=
3The drawable region of GLAC whenν=0.01.
(a)
α
=
-
3
(b)
α
=
-
2
(c)
α
=
-
1
(d)
α
=
-
0.5
(e)
α
=
-
0.1The drawable region of GLAC whenν=0.1.
(a)
α
=
-
3
(b)
α
=
-
2
(c)
α
=
-
1
(d)
α
=
-
0.5
(e)
α
=
-
0.1The drawable region of GLAC whenν=0.3.
(a)
α
=
-
3
(b)
α
=
-
2
(c)
α
=
-
1
(d)
α
=
0.5
(e)
α
=
-
0.1The drawable region forα=-1 as ν increases.
(a)
ν
=
-
0.3
(b)
ν
=
-
0.1
(c)
ν
=
-
0.01
(d)
ν
=
0.01
(e)
ν
=
0.1
(f)
ν
=
0.3The drawable region forα=0 as ν increases.
(a)
ν
=
-
0.3
(b)
ν
=
-
0.1
(c)
ν
=
-
0.01The drawable region forα=1 as ν increases.
(a)
ν
=
-
0.3
(b)
ν
=
-
0.1
(c)
ν
=
-
0.01The drawable region forα=2 as ν increases.
(a)
ν
=
-
0.3
(b)
ν
=
-
0.1
(c)
ν
=
-
0.01The drawable region of GLAC shows a better result when-1<ν<0 compared to when ν>0. This can be observed in Figure 4 where the drawable boundaries are larger compared to Figure 7. When ν is positive (small value, e.g., when ν=0.01) and -1<α<0, the drawable region is quite larger as compared to when -1<ν<0. However, whenν is positive, the drawable region gets smaller for all other α. This implies that the shape variable ν plays a major role in dictating the drawable region effectively. The drawable regions can be expended remarkably with a small change ofν. Generally, -1<ν<0 gives a wider drawable region when α>1 and α<-1. We can note that when α=2, the drawable region is very large compared to LAC. LAC has regions that have no solutions, whereas GLAC provides the solution.According to Yoshida and Saito [12], there are cases where the LAC segment is not drawable even when the second control point is within the drawable boundary. GLAC faces a similar situation where the second control point is within the drawable region yet the curve segment is not drawable.Figure12 is the comparison on the drawable region of LAC and GLAC when α=3. It portrays well the role of ν in GLAC. Note that GLAC becomes LAC when ν is 0. Since LAC is one of the cases in GLAC, the numerical example above clearly shows that for α>1, GLAC has a wider solution compared to LAC. Hence, GLAC provides greater flexibility where users would be able to create various shapes with a fixed start and end interpolating points.Comparison on the drawable region of LAC and GLAC forα=3.
(a)
GLAC(ν=0)
(b)
GLAC(ν=-0.01)


## 4. Conclusion and Future Works

An algorithm has been proposed to elucidate possible solutions that exist for interactive control of GLAC segment and the result is promising. The drawable region of GLAC proved to be far greater than LAC. This gives more freedom for designers to choose control points, tangent angle, and suitable shape parameters in order to design visually pleasing industrial products. The extra shape parameterν facilitates the ability to set the curvature radius at the origin which makes the GLAC more flexible.Future works include the development of aG2 algorithm for GLAC by manipulating the shapes variables similar to the work proposed in [13, 14]. A different kind of approach will be redefining the curvature function in such a way that a curvature continuous Log Aesthetic Spline (LAS) can be created.



---

*Source: 732457-2013-11-14.xml*