Next Article in Journal
On the Asymptotic Capacity of Information-Theoretic Privacy-Preserving Epidemiological Data Collection
Next Article in Special Issue
Modeling Categorical Variables by Mutual Information Decomposition
Previous Article in Journal
Integrated Multiresonator Quantum Memory
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Systematic Classification of Curvature and Feature Descriptor of 3D Shape and Its Application to “Complexity” Quantification Methods

1
School of Integrated Design Engineering, Graduate School of Keio University, Yokohama 223-8522, Japan
2
Department of Mechanical Engineering, Keio University, Yokohama 223-8522, Japan
*
Author to whom correspondence should be addressed.
Entropy 2023, 25(4), 624; https://doi.org/10.3390/e25040624
Submission received: 7 March 2023 / Revised: 31 March 2023 / Accepted: 31 March 2023 / Published: 6 April 2023
(This article belongs to the Special Issue Information Geometry for Data Analysis)

Abstract

:
Generative design is a system that automates part of the design process, but it cannot evaluate psychological issues related to shapes, such as “beauty” and “liking”. Designers therefore evaluate and choose the generated shapes based on their experience. Among the design features, “complexity” is considered to influence “aesthetic preference”. Although feature descriptors calculated from curvature can be used to quantify “complexity”, the selection guidelines for curvature and feature descriptors have not been adequately discussed. Therefore, this study aimed to conduct a systematic classification of curvature and a feature descriptor of 3D shapes and to apply the results to the “complexity” quantification. First, we surveyed the literature on curvature and feature descriptors and conducted a systematic classification. To quantify “complexity”, we used five curvatures (Gaussian curvature, mean curvature, Casorati curvature, shape index, and curvature index) and a feature descriptor (entropy of occurrence probability) obtained from the classification and compared them with the sensory evaluation values of “complexity”. The results showed that the determination coefficient between the quantified and sensory evaluation values of “complexity” was highest when the mean curvature was used. In addition, the Casorati curvature tended to show the highest signal-to-noise ratio (i.e., a high determination coefficient irrespective of the parameters set in the entropy calculation). These results will foster the development of generative design of 3D shapes using psychological evaluation.

1. Introduction

In recent years, computer-aided automatic design technology, referred to as “generative design”, has been put into practical use and applied to product design. Generative design is defined as a design system that automates part of the design process [1] or utilizes computing and manufacturing abilities to create novel and efficient designs [2]. Combined with processing technologies such as three-dimensional (3D) printers and additive manufacturing, generative design has enabled greater freedom in product design [3].
The advantage of generative design is that it enables the generation of shapes that satisfy physical characteristics such as stiffness and mass. However, its shortcoming is that psychological characteristics such as “beauty” and “preference” of the shape, which are important in the early process of a product design, are not evaluated. Therefore, designers need to evaluate and select shapes on the basis of their experience and intuition [4]. Solving this problem requires quantifying the shape features that designers focus on when evaluating a design.
When evaluating a design, designers tend to focus on macroscopic shape features, which are defined as those that appear in the entire shape as a result of combining points, lines, surfaces, and other modeling elements [5]. The main macroscopic shape features include “complexity”, “symmetry”, “order”, “proportion”, and “novelty” [6,7,8,9,10,11]. Among them, complexity and order are considered to influence “aesthetic preference” [12,13]. For example, in product design, high complexity is considered to be a factor that leads people to perceive a product as having high quality [14]. Moreover, complexity is considered a determinant of shape preference [15]. Thus, complexity is assumed to play an important role in the cognitive characteristics of shapes. Quantifying complexity will enable designers to design shapes that humans prefer.
Recent studies have indicated that a relationship exists between curvature and complexity [16,17,18]. Therefore, complexity can be quantified using a shape-feature value (hereafter referred to as a feature descriptor [19]) calculated from curvature. However, guidelines for selecting a curvature and feature descriptor to evaluate the cognition of a product shape’s complexity have not been adequately discussed. Accurately quantifying complexity requires studies on curvatures and feature descriptors of 3D shapes so that their characteristics can be understood and selection guidelines can be considered.
This study aims to perform a systematic classification of curvature and feature descriptors of 3D shapes and to apply the results to a method for quantifying complexity. This paper is organized as follows. Section 2 describes a systematic classification of curvature and feature descriptors of 3D shapes. Section 3 explains the quantification of complexity using the curvature and feature descriptor obtained from their classification. Section 4 illustrates sensory evaluation experiments that verify the applicability of the quantification. Section 5 summarizes the results and challenges of the study.

2. Systematic Classification of Curvature and Feature Descriptors

This section describes the systematic classification of curvature and feature descriptors of 3D shapes. We conducted a systematic classification according to a protocol known as preferred reporting items for systematic reviews and meta-analysis (PRISMA) [20]. One feature of this protocol, which is also used in the field of engineering, is that the selected references can be represented by flowcharts and tree diagrams [21,22]. The classification was conducted by three researchers from October 2021 to May 2022. The last access date to the database was 21 March 2022.

2.1. Classification Method

2.1.1. Keyword Setting

A paper search was conducted using the “Web of Science” database. Specifically, the search criteria were set to include all keywords related to “curvature”, “3D”, “shape”, and “shape features”, which were relevant to our study. Table 1 shows the keywords and search criteria.

2.1.2. Paper Selection Based on Title and Abstract

To narrow the papers obtained in the previous subsubsection to the literature highly relevant to the present study, we carried out a paper selection process based on titles and abstracts.
  • Papers evaluating 3D shapes were extracted. That is, we excluded both papers based on two-dimensional (2D) shapes and papers based on silhouettes or images of 3D shapes.
  • Among the references extracted in procedure 1, we selected those that evaluate shapes via, for example, classification, retrieval, and recognition based on the features extracted from the shapes. Our reasoning is that these references are highly relevant to the objective of the present study, which is to evaluate the complexity of shapes. Other references involving the generation/deformation of shapes, such as by meshing, rendering, or noise detection, were eliminated.
  • Among papers extracted in procedure 2, we selected those that evaluate macroscopic features by calculating a feature descriptor for each shape, such as similarity or complexity, because these references are highly related to evaluations of the degree of complexity.

2.1.3. Paper Selection Based on Content

From the papers selected as described in Section 2.1.1, we selected those based on textual content. Specifically, we selected papers that satisfy both of the following conditions:
  • Curvature is calculated from a 3D shape;
  • A feature descriptor is calculated on the basis of the curvature.
Papers that did not satisfy both criteria, such as those involving calculations of a curvature or feature descriptor, were excluded from our survey.

2.2. Classification Results

2.2.1. Paper Selection

We obtained 1310 papers as a result of keyword settings. Figure 1 shows an overview of the paper selection procedure and the number of references obtained in each procedure.
The details of the paper selection procedure are described below. First, 46 papers were selected on the basis of titles and abstracts. Figure 2 shows a tree diagram of papers classified on the basis of their titles and abstracts. The classification results are summarized as follows.
  • The shapes to be evaluated were classified into 3D and 2D shapes (generated from silhouettes or images of 3D shapes).
  • The objectives of the studies on 3D shapes were divided into three main categories: (1) shape evaluation, (2) shape generation, and (3) shape deformation. Shape evaluation is the process of extracting features from shapes and evaluating their characteristics (i.e., similarity evaluation classification, retrieval, recognition, posture estimation, registration, abnormal detection, simulation, and segmentation). Shape generation refers to the generation of 3D shapes by rendering or mesh segmentation of data such as 2D images or point clouds. Shape deformation denotes the process of removing points that deviate from the actual shape by detecting noise in the point cloud acquired by 3D scanning.
  • The tasks for shape evaluation were further classified into multiple shape evaluations (e.g., similarity calculations, classification, retrieval, and recognition) and single shape evaluations (e.g., posture estimation, registration, abnormal detection, simulation, and segmentation).
In the present study, following the criteria described in Section 2.1.2, we selected papers that targeted 3D shapes and aimed to evaluate shapes, especially by comparing them among numerous shapes. The criteria of selected papers are indicated by the bold lines in Figure 2. Finally, 46 papers were checked on the basis of the content of their text, and those without curvatures or feature descriptors were excluded, resulting in 18 papers. Table 2 lists the reference numbers, curvatures, feature descriptors, authors, and publication years of the papers.
In these papers, polygon models were used to represent 3D shapes. A polygon model is a method of representing 3D shapes using a set of minute polygons, such as triangles and quadrilaterals. Because polygon models contain information about the shape surface necessary for calculating curvature and are commonly used in the field of industrial design [23,24], polygon models are considered a useful representation method in the present study. We therefore employed polygon models.
Table 2. Target papers.
Table 2. Target papers.
ReferenceCurvatureType of Feature Descriptor
(Information Considered)
AuthorYear
[25]Cone curvatureCurvatures of vertices and their surroundingsAdán, M2003
[26]Cone curvatureCurvatures of vertices and their surroundingsAdán, M2008
[27]Cone curvatureOccurrence probabilities of curvaturesAdán, M2014
[28]Casorati curvatureOccurrence probabilities of curvaturesWang, J2020
[29]Gaussian curvature
Mean curvature
Clusters of curvaturesBiasotti, S2015
[30]Shape index
Curvature index
Occurrence probabilities of curvaturesLi, B2014
[31]Gaussian curvatureCurvatures of vertices and their surroundingsFu, J2008
[32]Shape indexOccurrence probabilities of curvaturesJeong-Jun, S2003
[33]Shape indexCurvatures of verticesJunli, Z2014
[34]Shape indexCurvatures of vertices and their surroundingsLi, B2011
[35]Curvature tensorCurvatures of verticesMuzahid, M2021
[36]Cone curvatureTransition probabilityYuJie, L2013
[37]Mean curvatureClusters of curvaturesZou, K2015
[38]Gaussian curvatureTransition probabilities of curvaturesGuo, K2010
[39]Gaussian curvatureOccurrence probabilities of curvaturesSukumar, S2006
[40]Gaussian curvature
Mean curvature
Transition probabilities of curvaturesGuo, K2014
[41]Gaussian curvatureOccurrence probabilities of curvaturesSukumar, S2006
[42]Gaussian curvatureTransition probabilities of curvaturesMatsumoto, T2018

2.2.2. Systematic Classification and Overview of Curvature

In the 18 extracted papers, seven types of curvature (i.e., mean curvature, Gaussian curvature, shape index, curvature index, Casorati curvature, curvature tensor, and cone curvature) were used. Table 3 shows the calculation method and characteristics of each curvature. Most curvatures can be calculated using two curvatures: the maximum principal curvature k 1 and the minimum principal curvature k 2 ( k 1 > k 2 ) of the normal curvature created when the surface is cut at an arbitrary point (Figure 3).

2.2.3. Systematic Classification and Overview of Feature Descriptor

In the 18 extracted papers, five types of feature descriptor were used. Table 4 shows the advantages and disadvantages of each feature descriptor.
  • A feature descriptor that uses the curvature of each vertex of a polygonized shape is a vector whose elements are the curvatures at all vertices [33]. Compared with the other feature descriptors, a feature descriptor that uses the curvature of each vertex of a polygonized shape has the advantage of being computationally efficient but has the disadvantage of being unable to compare shapes with different numbers of vertices unless calculating statistic values (e.g., average and probability) [35].
  • A feature descriptor that uses the curvature of each vertex and its surrounding vertices is a vector whose elements are calculated from the curvature of each vertex and its surrounding vertices, or it is a scalar obtained by their calculation (integration, weighted addition, etc.) [31]. Compared with other feature descriptors without the surrounding vertices, a feature descriptor that uses the curvature of each vertex and its surrounding vertices has the advantage of enabling an accurate evaluation of imperceptible surface changes but has the disadvantage of low computational efficiency [25].
  • A feature descriptor that uses occurrence probabilities is a vector whose elements are the occurrence probabilities of curvature within a set range of curvatures, or entropy calculated from the vector [39]. Compared with other feature descriptors that use curvature values, a feature descriptor that uses occurrence probabilities has the advantage of allowing the user to set the range of curvatures of interest; however, it has the disadvantage that judging the validity of the curvature range setting and setting it appropriately are difficult [28].
  • A feature descriptor that uses transition probability is a vector of transition probabilities (probabilities of changing from one curvature to another) between two adjacent vertices within a set curvature range, or the entropy calculated from this vector, similar to occurrence probability [39]. Compared with feature descriptors that use the occurrence probabilities, a feature descriptor that uses transition probability has the advantage of being able to accurately evaluate surfaces within a localized range because it can take into account the proportion of changes from the curvature at surrounding vertices; however, it has the disadvantage of being computationally inefficient [42].
  • A feature descriptor that uses clustering methods (e.g., K-means and C-means methods) can reduce the information amount by clustering vector-type curvatures such as the cone curvature and curvature tensor to clarify the remarkable vertices of the shape. This descriptor cannot be used for scalar curvatures.

3. Application to a Method for Quantifying Complexity of 3D Shapes

This section describes a method for quantifying the complexity of 3D shapes. Section 3.1 describes how to calculate the five curvature indices: Gaussian curvature, mean curvature, Casorati curvature, shape index, and curvature index. Section 3.2 describes how to calculate feature descriptors using occurrence probability.

3.1. Methods for Calculating Curvature

This subsection describes the calculation methods of the curvatures in Table 3 from the polygon model. This study extracted five curvatures: Gaussian curvature, mean curvature, Casorati curvature, shape index, and curvature index. The curvature tensor (Figure 4a) and cone curvature (Figure 4b) were excluded for the following reasons:
  • They have a large computational load to calculate the curvature. Both curvatures are expressed as an array of curvatures for a single vertex. Therefore, they are computationally more expensive than other curvatures where a single value is calculated for a single vertex. In addition, the computation required to compute feature descriptors from curvatures can also be substantial.
  • They are difficult to use to evaluate curvatures of different sizes. A region (surrounding the vertex for which the curvature is calculated) must be set as a parameter. However, in actual design situations, many scenarios arise where shapes of different sizes are evaluated, and setting a common region is difficult.
The following illustration of the curvature calculation is in accordance with the polygon mesh in Figure 5. In this figure, v i is the vertex whose curvature is to be calculated, v i k ( k = 1 ,   2 , , n i ) are vertices adjacent to vertex v i , f i k ( k = 1 ,   2 , , n i ) are polygons with v i as a vertex, and a i k = ( v i k , v i , v i ( k + 1 ) ) is the angle of f i k at vertex v i .

3.1.1. Gaussian Curvature

Using Gauss–Bonnet’s law and paraboloid fitting, Gaussian curvature K i at the ith vertex v i is calculated via the following equation [43]:
K i = 2 π k = 1 n i α i k 1 3 A
where A is the sum of the polygon area surrounding vertex v i .

3.1.2. Mean Curvature

Using Gauss–Bonnet’s law and paraboloid fitting, mean curvature H i at the ith vertex v i is calculated as follows [43]:
H i = 1 4 k = 1 n e i k β i k 1 3 A
where e i k is the distance between vertex v i and its neighboring vertices v i k ( k = 1 ,   2 , , n i ) and β i k is the angle between adjacent polygons.

3.1.3. Casorati Curvature

The mean curvature is the curvature defined by the average of the principal curvatures, whereas the Gaussian curvature is the curvature defined by the product of the principal curvatures. Therefore, principal curvatures k 1 and k 2 can be calculated by combining Equations (1) and (2). When the maximum principal curvature is k 1 and the minimum principal curvature is k 2 , both are expressed by the following equations:
k 1 = H + H 2 K
k 2 = H H 2 K
The Casorati curvature C can be calculated using the principal curvatures obtained from Equations (3) and (4):
C = k 1 2 + k 2 2 2

3.1.4. Shape Index

The shape index can be calculated using the principal curvatures obtained from Equations (3) and (4). However, such an equation cannot be used to calculate the shape index in the plane ( k 1 = k 2 = 0 ) . Therefore, in the present study, the shape index SI can be calculated even for a flat surface by dividing the cases according to whether k 1 = k 2 = 0 or not (Equation (6)):
S I = { 1 2 1 π a r c t a n ( k 1 + k 2 k 1 k 2 )   ( k 1 0 ,   k 2 0 ) 1 2   ( k 1 = k 2 = 0 )

3.1.5. Curvature Index

The curvature index can be calculated by using the principal curvatures obtained in Equations (3) and (4). However, the resultant equation cannot be used to calculate the curvature index in the plane ( k 1 = k 2 = 0 ) . Therefore, in the present study, the curvature index CI can be calculated even for a flat surface by dividing the cases according to whether k 1 = k 2 = 0 or not (Equation (7)):
C I = { 2 π l o g ( k 1 2 + k 2 2 2 )   ( k 1 0 ,   k 2 0 ) 2 π l o g ( 1 e 15 2 )   ( k 1 = k 2 = 0 )

3.2. Method for Calculating Feature Descriptor

Sukumar et al. [41], who investigated systematic classification, suggested that the complexity of a shape could be quantified by considering its occurrence probability. In this approach, the Gaussian curvature at each vertex of a 3D shape is first calculated and the curvature is subsequently discretized by setting the number of states and deviation to express the continuous Gaussian curvature as discrete values. Finally, the information (occurrence probability) entropy of the discretized curvature is calculated as a feature descriptor to quantify complexity. Such discretization can be performed according to the curvature changes that can be discriminated by humans and the magnitude of the curvature that attracts their attention. Therefore, it is considered useful for evaluating human cognition. However, feature descriptors that use curvature values might be strongly affected when outliers (especially large or small curvatures) occur in a small area or when minute curvature changes that cannot be discriminated by humans occur in a wide area. Therefore, we decided to adopt a feature descriptor based on discretized occurrence probability.
Feature descriptors that consider surrounding curvature information, those that consider transition probabilities, and those that consider clustering were excluded from our calculation. Our reasoning is described below.
  • The feature descriptors that consider the surrounding information are excluded because the computational load to integrate the curvature of each vertex and its neighboring vertices into a single value by calculation is larger than that of the other descriptors. In addition, the purpose of considering surrounding information is to evaluate the similarity of shapes accurately, not to evaluate human cognition.
  • Feature descriptors that consider transition probabilities are excluded because evaluating the surface features of the entire shape is difficult. Although the feature descriptor can evaluate the transition of the curvature from one vertex to another, people are unlikely to focus only on local changes between vertices when evaluating complexity. In addition, checking and considering the results based on graphs such as contour plots are difficult when using transition probabilities.
  • Feature descriptors that consider clustering are excluded because evaluating the curved surface features of the entire shape is difficult. The feature descriptor extracts salient points (vertices with characteristic curvature) from the entire shape as local features to evaluate the similarity among shapes deformed in only one part. However, curvatures not considered salient (not complex) are not reflected in the feature descriptors. In addition, as mentioned in the previous section, clustering can only be applied to a vector curvature (e.g., cone curvature and curvature tensors), where multiple values are calculated from a single vertex.
In this subsection, we first provide an overview of information entropy and its calculation method. We then outline the method used to calculate feature descriptors using the probability of occurrence in the present study.

3.2.1. Information Entropy

Information entropy E expresses information clutter by the following equation [44]:
E = i = 1 n p i log p i
where n is the number of events and p i is the probability that event i is likely to occur. The information entropy increases as the probability of occurrence of each event is equally certain.

3.2.2. Calculation Methods

  • The curvature at each vertex of the 3D shape to be evaluated is calculated as described in Section 3.1.
  • The curvature at each vertex is discretized. Specifically, the minimum area value E min , the maximum area value E max , and the number of states V are first set as parameters. The range from the minimum to the maximum area is then equally divided by the number of states to create V states s i   ( i = 1 , 2 , , V ) . Finally, the curvature of each vertex is assigned to a state.
    s 1   ( E min     K   < E min + Δ E ) ,   s 2   ( E min + Δ E     K   < E min + 2 Δ E ) ,   ,   s V   ( E max     Δ E     K   <   E max ) ( Δ E = ( E max E min ) V )
  • After discretization, occurrence probability p i is calculated using the following equation:
    p i = N i N
    where N is the total number of vertices in the 3D shape to be evaluated and N i is the number of vertices assigned to state s i .
  • The entropy of the probability of occurrence is calculated using Equation (8).

4. Sensory Evaluation Experiment

This section describes the sensory evaluation experiments conducted to analyze the relationship between the feature descriptors calculated from five curvatures and sensory evaluation values of complexity. We first describe the experimental method, then present the results and discuss the sensory evaluation experiments conducted on three types of shapes.

4.1. Experimental Methods

4.1.1. Sample Shapes

In this experiment, extruded and rotated shapes, which are typical CAD shape deformation methods, were used, along with the shapes of an actual product as the “complexity” evaluation sample. The details of the shapes are described below.
  • Extruded shape
Extruded geometry is a 3D curved surface created from closed planar curves using extrusion, which is a 3D CAD geometry creation method. The extruded geometry was chosen for evaluation in this study because it enables the curvature characteristics of the planes, contours, and sides of an extruded surface to be compared. In this experiment, 50 2D closed curves used by Ujiie et al. [16] were edited in 3D CAD and the 3D shapes to be evaluated were created by extruding the surfaces in the vertical direction of the same shape (Figure 6a). The size of the closed planar curves was adjusted so that their maximum radius vector was equal, and the length of the extrusion was set to twice the maximum radius vector of the closed planar curves.
  • Rotated shape
A rotated shape is a three-dimensional curved surface created from closed planar curves using rotation, which is a shape creation method of 3D CAD. Rotational geometry was chosen as the target of evaluation in the present study because it enables the differences in curvature among planar, convex, and concave surfaces to be considered.
The rotated shape used in this experiment was created by editing 50 closed planar curves similar to the extruded shape in 3D CAD and rotating the right side of the shape 360 degrees using the vertical direction passing through the center of gravity of the 2D shape as the rotation axis (Figure 6b). The size of the closed planar curves was adjusted so that the maximum radius vectors of the curves were equal.
  • Vase shape
A vase is a product with a high degree of freedom in its shape, and various types of curved surfaces can appear. Therefore, using the vase shape as the evaluation target, we can consider the influence of various curved surfaces on the perception of complexity. In addition to its functional role of holding cut flowers, a vase also plays a role as an interior decoration; thus, it is important to take aesthetic preferences into account when modeling [45]. A vase is therefore considered an appropriate target for evaluation in the present study, which aims at aesthetic preferences in product shape.
A total of 50 vase shapes were used in this experiment: 25 vase shapes provided by Free3D [46] and 25 vase shapes provided by CGTrader [47]; the top of each vase was closed with a flat surface to prevent curvature from occurring inside the vase, which is not visible (Figure 6c).

4.1.2. Experimental Conditions

  • Evaluation method: A 5-point Likert scale was used for complexity: “low complexity: 1”, “somewhat low complexity: 2”, “undecided: 3”, “somewhat high complexity: 4”, and “high complexity: 5”.
  • Presentation method: An online questionnaire was used, with the sample shape colored gray and displayed on a white background. In addition, 10 of each sample shape were randomly selected and presented in duplicate to check the accuracy of the participants’ evaluation. Each shape was rotated at a constant speed (z-axis: 9.0 s per rotation).
  • Participants: 110 male and 110 female participants. To exclude participants with inaccurate evaluations, we calculated the absolute value of the difference in sensory evaluation values among the 10 shapes presented in duplicate and excluded participants whose sum was 11 or more. As a result, data from 81 participants for the extruded shape, 89 participants for the rotated shape, and 61 participants for the vase shape were used.
To calculate the curvature and feature descriptors from the geometry under evaluation, we performed the following additional procedures in the present study.
  • Equal division of polygons
The shapes created by the method described in Section 4.1.1 had polygons with nonuniform shapes and sizes. Because participants were expected to evaluate complexity from the entire shape, each polygon vertex should be uniformly distributed. In addition, because the area of polygons was involved in the formulas for calculating Gaussian curvature and mean curvature, the curvature value might be affected by the size of the area.
Therefore, in the present study, polygons were equally divided by applying the advancing-front method to the created geometry. An advantage of the advancing-front method is that the polygons are closer to equilateral triangles than those used in other methods of polygon division, such as the Delaunay division and bubble mesh [48]. In the present study, equal division was performed using software [49] that implements the method, and the results were analyzed for 50 extruded shapes, 46 rotated shapes, and 38 vase shapes for which equal division was possible.
2.
Determination of discretization parameters
When calculating the entropy of occurrence probability, discretization is required to divide the continuous curvature into multiple states. In the present study, suitable ranges were set for each shape and curvature by the following methods:
  • The curvature was calculated at all vertices of the shape to be evaluated.
  • A percentile range [ E 0.003 % , E 99.997 % ] that included curvatures between 0.003% and 99.997% of the calculated curvatures was calculated. This range was determined with reference to the range of (mean) ± 4 (standard deviation) for the normal distribution.
  • The curvatures included in the determined percentile range [ E 0.003 % , E 99.997 % ] were divided by V , as in the conventional method [39]. The curvatures above the maximum and below the minimum in the same range were assigned to the state with the largest and smallest curvatures, respectively, as in the method of Matsumoto et al. [42].
  • Integers from 2 to 20 were used as candidates for the number of states.

4.2. Experimental Results

4.2.1. Extruded Shape

Table 5 shows the results obtained when calculating the entropy of occurrence probability using five different curvatures for extruded shapes. The table lists as discretization parameters the maximum range E max , the minimum range E min , and the number of states V when the determination coefficient of the logarithmic approximation between the sensory evaluation values and the entropy is maximal. Note that the logarithmic approximation is applied based on Fechner’s law, which indicates the relationship between human sensitivity and stimuli using the logarithmic function and used in shape cognition studies [17,18,42,50]. In addition, to evaluate the stability of the discretization, the table lists the average, standard deviation, and the larger-the-better signal-to-noise (SN) ratio of the determination coefficients for the range of 5 to 20 states. Figure 7 shows the relationship of the determination coefficients with complexity versus the number of states. Table 5 and Figure 7 indicate that the highest correlation is obtained when the mean curvature is used, with a determination coefficient of 0.704. However, the determination coefficient is stable with respect to the number of states (the SN ratio is highest) when the Casorati curvature is used.

4.2.2. Rotated Shape

Table 6 shows the results of calculations of the entropy of occurrence probability using five different curvatures for rotated shapes. Figure 8 shows the determination coefficient between complexity and the number of states. Table 6 and Figure 8 indicate that the highest correlation is obtained when the mean curvature is used, with a determination coefficient of 0.551. However, the determination coefficient is stable with respect to the number of states (the SN ratio is highest) when the Casorati curvature is used.

4.2.3. Vase Shape

Table 7 shows the results of calculations of the entropy of occurrence probability using five types of curvature for vase shapes. Figure 9 shows the determination coefficient between the number of states and complexity. Table 7 and Figure 9 indicate that the highest correlation is obtained when the mean curvature is used, with a determination coefficient of 0.471. However, the determination coefficient is stable with respect to the number of states (the SN ratio is highest) when the Casorati curvature is used.

4.3. Experimental Discussion

4.3.1. Variation in Determination Coefficients with Different Numbers of States

When the Gaussian curvature and mean curvature were used, the determination coefficient with complexity tended to vary with the number of states V . This variation was observed in all three sample shapes, suggesting that the characteristics of the two curvatures affect the result. In the following, we focus on the number of states with large and small determination coefficients for the Gaussian curvature and mean curvature and discuss the reasons for the variation.
  • Gaussian curvature
As an example, we compare the number of states 6 with 5, where the change in the determination coefficient is large in rotated shapes. In the rotated shape, where the entropy variation is particularly large, the surface is discretized for each number of states as follows (Figure 10):
  • For V = 6 with a high determination coefficient, flat surfaces and surfaces with a principal curvature of 0 at one side are discretized into gray states, whereas convex surfaces are discretized into red states (Figure 10a).
  • For V = 5 with a low determination coefficient, most of the surfaces are in the gray states (Figure 10b).
In rotational shapes, steep convex surfaces tend to influence the evaluation of complexity. Therefore, the number of states 6 that can discriminate complex surfaces (convex surface) from uncomplex surfaces (i.e., the surfaces whose maximum and/or minimum principal curvature is almost 0) is consistent with the human perception of complexity and is considered to have a high determination coefficient. Similarly, in the extruded and vase shapes, the number of states that discriminate complex and uncomplex surfaces tends to have a high determination coefficient (Figure 10c,d).
2.
Mean curvature
As an example, we compare the number of states 7 with 6, where the change in the determination coefficient is large in the extruded shapes. In the extruded shapes, where the entropy variation is particularly large, the surface is discretized for each number of states as follows (Figure 11):
  • For V = 7 with a high determination coefficient, flat surfaces and moderate surfaces with a principal curvature of 0 at one side are discretized into gray states. However, steep surfaces with a principal curvature of 0 at one side are discretized into red areas (Figure 11a).
  • For V = 7 with a high determination coefficient, flat surfaces and steep surfaces with a principal curvature of 0 at one side are discretized into gray states. However, moderate surfaces with a principal curvature of 0 at one side are discretized into blue areas (Figure 11b).

4.3.2. Variation in Determination Coefficients with Different Curvatures

The experimental results show that the Casorati curvatures are stable and highly correlated for all the investigated shapes. However, depending on the number of states, the determination coefficients of the Gaussian curvature and mean curvature can exceed the determination coefficient of the Casorati curvature. Therefore, we discuss the factors by comparing the determination coefficients of the Gaussian curvature and mean curvature with those of the Casorati curvature at the number of states where the determination coefficients of the Gaussian curvature and mean curvature are maximal.
  • Gaussian curvature
As an example, in the extruded geometry, the highest determination coefficient when using the Gaussian curvature is found for a number of states 18. We focus on the extruded shape, where the entropy change between the two curvatures is large. At each curvature, the surfaces of extruded shape are discretized as follows (Figure 12).
  • For the Gaussian curvature, moderate surfaces with a principal curvature of 0 at one side and flat surfaces are discretized into gray states.
  • For the Casorati curvature, moderate surfaces with a principal curvature of 0 at one side are discretized into different states: blue, gray, and red.
Unlike the Gaussian curvature, the Casorati curvature can take positive values in the column plane. In addition, compared with the Gaussian curvature, the Casorati curvature tends to be larger for small (near zero) changes in principal curvature. Therefore, for extruded shapes, the Gaussian curvature is expected to quantify complexity more accurately than the Casorati curvature.
2.
Mean curvature
As an example, for the vase shape, the highest determination coefficient when using the mean curvature is found when the number of states is 19. We focus on the vase shape, where the entropy change between the two curvatures is large. For each curvature, the surfaces of the vase shape are discretized as follows (Figure 13):
  • For the mean curvature, flat and principal curvatures of 0 at one side are discretized into gray states, whereas convex edges are assigned to the red area and concave edges to the blue area.
  • For the Casorati curvature, flat surfaces and concave surfaces are discretized into the same blue area. Convex surfaces as edges are assigned to gray or red states.
Vase shapes tend to have a higher complexity evaluation when both convex and concave surfaces are included. Therefore, mean curvatures with different states for convex and concave surfaces are considered appropriate.

4.3.3. Variation in Determination Coefficients with Different Feature Descriptors

In the present study, the entropy of occurrence probability was calculated as a feature descriptor because it can be used to evaluate the diversity of curvature and because the discretization parameters can be set according to human cognition. To verify the validity of this feature descriptor, we calculated the curvature moment as a comparison target.
The curvature moment is a characteristic value calculated by considering the curvature group calculated from the shape as a probability distribution. Because some studies have evaluated the similarity of shapes using moments of curvature [51], it is applicable to the quantification of complexity in the present study. Here, the average, standard deviation, coefficient of variation, kurtosis, and skewness were calculated as moments.
Table 8 shows the determination coefficients when the entropy of occurrence probability is used and when the curvature moments are used. The table shows the average, standard deviation, and desired characteristic SN ratio of the 15 determination coefficients (3 shapes × 5 curvatures) calculated for each feature descriptor. The results in this table show that the occurrence probability has large and stable determination coefficients with complexity.
The lower determination coefficient when using the curvature moment is attributed to the effect of mesh irregularities. In this subsection, we focus on the standard deviation that exhibits the highest desirability characteristic SN ratio among the curvature moments. The shape shown in Figure 14 is evaluated to be of medium complexity; however, the standard deviation is large. This shape is uneven in the blue area in Figure 14, and the standard deviation is considered to have been increased because of the effect of this area. Because the occurrence probability entropy is discretized, the curvature value does not affect the entropy value. Therefore, we speculate that the high determination coefficient is attributable to the fact that the effect of the curvature of such outliers can be reduced.
One limitation of this study is that only the determination coefficient of the logarithmic approximation between the quantified and sensory evaluation values of complexity was used to evaluate their correlation. Although we confirmed some scatter diagrams (Appendix A) and this procedure is assumed to be valid, other procedures can be more suitable.

5. Conclusions

We conducted a systematic classification of curvature and feature descriptors in 3D shapes. The results showed that Gaussian curvature, mean curvature, Casorati curvature, curvature index, shape index, curvature tensor, and cone curvature can be used as curvatures calculated for 3D shapes. In addition, we confirmed that “feature descriptors considering surrounding information”, “feature descriptors considering occurrence probability”, “feature descriptors considering transition probability”, and “feature descriptors considering clustering” could be used as feature descriptors calculated for the 3D shapes.
We applied the results of the systematic classification to a method for quantifying the “complexity” of 3D shapes. Specifically, Gaussian curvature, mean curvature, Casorati curvature, curvature index, and shape index were calculated for the polygon model. In addition, the probability of occurrence (information) entropy of these five types of curvatures was calculated as a feature descriptor.
Finally, we conducted sensory evaluation experiments to verify the validity of the quantification method; the results showed that the highest determination coefficient for all three shapes was obtained when the mean curvature was used. To investigate the effect of discretization when calculating entropy, we calculated the mean value, standard deviation, and the SN ratio of the determination coefficient when the number of states was varied. The results showed that the Casorati curvature was stable for almost all of the investigated shapes.
Our future work will include (1) conducting sensory evaluation experiments using other sample shapes to verify whether the same trend can be confirmed and (2) constructing a comparative criterion to enable a comparison that includes other types of curvatures and feature descriptors excluded in the present study.

Author Contributions

Conceptualization, K.M. and T.K.; formal analysis, K.M. and T.S.; investigation, K.M. and T.S.; data curation, K.M.; investigation, K.M. and T.S.; resources, T.K.; writing—original draft preparation, K.M. and T.K.; writing—review and editing, T.K.; supervision, T.K.; project administration, T.K.; funding acquisition, T.K. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by the Japan Society for the Promotion of Science Grant-in-Aid for Scientific Research (C) (20K12514) and Scientific Research (B) (21H03528).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

This study was partly supported by the Japan Society for the Promotion of Science Grant-in-Aid for Scientific Research (C) (20K12514) and Scientific Research (B) (21H03528).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Figure A1 shows the relationship between the entropy using mean curvature and sensory evaluation values.
Figure A1. Scatter diagrams with entropy values using mean curvature on the x axis and the sensory evaluation values on the y axis: (a) the number of states V = 7; (b) V = 6.
Figure A1. Scatter diagrams with entropy values using mean curvature on the x axis and the sensory evaluation values on the y axis: (a) the number of states V = 7; (b) V = 6.
Entropy 25 00624 g0a1

References

  1. Singh, V.; Gu, N. Towards an integrated generative design framework. Des. Stud. 2012, 33, 185–207. [Google Scholar] [CrossRef]
  2. Shea, K.; Aish, R.; Gourtovavia, M. Towards integrated performance-driven generative design tools. Autom. Constr. 2005, 14, 253–264. [Google Scholar] [CrossRef]
  3. Lu, P.; Shih-Wen, H.; Fan, W. A Product Shape Design and Evaluation Model Based on Morphology Preference and Macroscopic Shape Information. Entropy 2021, 23, 639. [Google Scholar] [CrossRef]
  4. Krish, S. A practical generative design method. Comput. Aided Des. 2011, 43, 88–100. [Google Scholar] [CrossRef]
  5. Lin, M.H.; Lee, C.L. An Experimental Study for Applying Generative Design to Electronic Consumer Products. In Proceedings of the International Conference of Design, User Experience, and Usability, Las Vegas, NV, USA, 21–26 July 2013. [Google Scholar]
  6. Arnheim, R. Toward a Psychology of Art; University of California Press: Berkeley, CA, USA, 1966. [Google Scholar]
  7. Eysenck, H.J. The experimental study of the good Gestalt. Psychol. Rev. 1942, 49, 344–364. [Google Scholar] [CrossRef]
  8. Garner, W.R. The Processing of Information Structure; Psychology Press: New York, NY, USA, 1975. [Google Scholar]
  9. Gombrich, E.H. Art and Illusion: A Study in the Psychology of Pictorial Representation; Phaidon: London, UK, 1960. [Google Scholar]
  10. Gombrich, E.H. The Sense of Order: A Study in the Psychology of Decorative Art; Cornell University Press: New York, NY, USA, 1979. [Google Scholar]
  11. Solso, L. Cognition and the Visual Arts; MIT Press: Cambridge, MA, USA, 1994. [Google Scholar]
  12. Berlyne, D.E. Aesthetics and Psychobiology; Appleton Century Crofts: New York, NY, USA, 1971. [Google Scholar]
  13. Birkhoff, G.D. Aesthetic Measure; Harvard University Press: Cambridge, UK, 1933. [Google Scholar]
  14. Creusen, M.E.H.; Veryzer, R.W.; Schoormans, J.P.L. Product value importance and consumer preference for visual complexity and symmetry. Eur. J. Mark. 2010, 44, 1437–1452. [Google Scholar] [CrossRef]
  15. Blijlevens, J.; Creusen, M.E.H.; Schoormans, J.P.L. How consumers perceive product appearance: The identification of three product appearance attributes. Int. J. Des. 2009, 3, 27–35. [Google Scholar]
  16. Ujiie, Y.; Matsuoka, Y. Total absolute curvature to represent the complexity of diverse curved profiles. In Proceedings of the 6th Asian Design Conference, Tsukuba, Ibaraki, Japan, 14–17 October 2003. [Google Scholar]
  17. Ujiie, Y.; Kato, T.; Sato, K.; Matsuoka, Y. Curvature entropy for curved profile generation. Entropy 2012, 14, 533–558. [Google Scholar] [CrossRef] [Green Version]
  18. Okano, A.; Matsumoto, T.; Kato, T. Gaussian Curvature Entropy for Curved Surface Shape Generation. Entropy 2020, 22, 353. [Google Scholar] [CrossRef] [Green Version]
  19. Farlin, M.; Ranga, R. Curvature Based Robust Descriptors. In Proceedings of the British Machine Vision Conference, Guildford, UK, 2–7 September 2012. [Google Scholar]
  20. Page, J.M. The PRISMA 2020 statement: An updated guideline for reporting systematic reviews. Syst. Rev. 2021, 10, 89. [Google Scholar] [CrossRef]
  21. James, I.N.; Andrew, N. Is Additive Manufacturing Performance in Sports? A Systematic Review. Proceeding Inst. Mech. Eng. Part P J. Sport. Eng. Technol. 2020, 235, 163–175. [Google Scholar]
  22. Farashi, S.; Vafaee, F. Effect of printing parameters on the tensile strength of FDM 3D samples: A meta-analysis focusing on layer thickness and sample orientation. Prog. Addit. Manuf. 2022, 7, 565–582. [Google Scholar] [CrossRef]
  23. Antonelli, M.; Beccari, C.V.; Casciola, G.; Ciarloni, R. Subdivision surfaces integrated in a CAD system. Comput. Aided Des. 2013, 45, 1294–1305. [Google Scholar] [CrossRef]
  24. Ma, W. Subdivision surfaces for CAD—An overview. Comput. Aided Des. 2004, 37, 693–709. [Google Scholar] [CrossRef]
  25. Adan, M.; Adan, A.; Cerrada, C.; Merchan, P.; Salamanca, S. Weighted Cone-Curvature: Applications for 3D Shapes Similarity. In Proceedings of the Fourth International Conference on 3-D Digital Imaging and Modeling, Banff, AB, Canada, 6–10 October 2003. [Google Scholar]
  26. Adán, A.; Adán, M.; Santiago, S.; Pilar, M. Using Nonlocal Features for 3D Shape Grouping. In Planning Perspectives; IEEE: Piscataway Township, NJ, USA, 2008; pp. 644–653. [Google Scholar]
  27. Adán, A.; Adán, M. Consensus Strategy for Clustering Using RC-Images. Pattern Recognit. 2014, 47, 402–417. [Google Scholar] [CrossRef]
  28. Wang, J.; Yan, W.; Huang, C. Surface shape-based clustering for B-rep models. Multimed. Tools and Appl. 2020, 79, 25747–25761. [Google Scholar]
  29. Biasotti, S.; Andrea, C.; Bianca, F.; Michela, S. 3D Artifacts Similarity Based on the Concurrent Evaluation of Heterogeneous Properties. J. Comput. Cult. Herit. 2015, 8, 1–19. [Google Scholar] [CrossRef]
  30. Bo, L.; Godil, A.; Johan, A. Hybrid Shape Descriptor and Meta Similarity Generation for Non-Rigid and Partial 3D Model Retrieval. Multimed. Tools Appl. 2014, 72, 1531–1560. [Google Scholar]
  31. Fu, J.; Sanjay, B.J.; Timothy, W.S. Shape differentiation of freeform surfaces using a similarity measure based on an integral of Gaussian curvature. Comput. Aided Des. Appl. 2008, 40, 311–623. [Google Scholar] [CrossRef]
  32. Song, J.; Golshani, F. Shape-Based 3D Model Retrieval. Multimed. Tools Appl. 2020, 79, 25747–25761. [Google Scholar]
  33. Zhao, J.; Cuiting, L.; Wu, Z.; Duan, F.; Zhang, M.; Wang, K.; Taorui, J. 3D Facial Similarity Measure Based on Geodesic Network and Curvatures. Math. Probl. Eng. 2014, 2014, 832837. [Google Scholar] [CrossRef] [Green Version]
  34. González, J.A.; Navarro, M.; Hernández, H. Shape Descriptor Based on Curvature. Open Access Libr. J. 2022, 9, 1–13. [Google Scholar]
  35. Muzahid, A.A.M.; Wanggen, W.; Ferdous, S.; Lianyao, W.; Li, H. CurveNet: Curvature-Based Multitask Learning Deep Networks for 3D Object Recognition. IEEE/CAA J. Autom. Sin. 2021, 8, 1177–1187. [Google Scholar] [CrossRef]
  36. Yujie, L.; Xiao-Dong, Z.; Zongmin, L.; Hua, L. Extended cone-curvature based salient points detection and 3D model retrieval. Multimed. Tools Appl. 2013, 64, 671–693. [Google Scholar]
  37. Zou, K.; Zhaojun, Z.; Jianhua, Z.; Qian, Z. A 3D Model Feature Extraction Method Using Curvature-based Shape Distribution. In Proceedings of the 12th International Conference on Fuzzy Systems and Knowledge Discovery (FSKD), Zhangjiajie, China, 15–17 August 2015. [Google Scholar]
  38. Guo, K. 3D Shape Representation Using Gaussian Curvature Co-occurrence Matrix. Lect. Notes Comput. Sci. 2010, 6319, 373–380. [Google Scholar]
  39. Sukumar, S.; Page, D.; Koschan, A.F.; Roui-Abidi, B.; Abidi, M. Shape analysis algorithm based on information theory. In Proceedings of the International Conference on Image Processing, Barcelona, Spain, 14–17 September 2003. [Google Scholar]
  40. Guo, K.; Duan, G. 3D image retrieval based on differential geometry and co-occurrence matrix. Neural Comput. Appl. 2014, 24, 715–721. [Google Scholar] [CrossRef]
  41. Sukumar, S.; Page, D.; Gribok, A.; Koschan, A.; Abidi, M. Shape Measure for Identifying Perceptually Informative Parts of 3D Objects. In Proceedings of the International Symposium on 3D Data Processing, Visualization, and Transmission, Chapel Hill, NC, USA, 14–16 June 2006. [Google Scholar]
  42. Matsumoto, T.; Sato, K.; Matsuoka, Y.; Kato, T. Quantification of “complexity” in curved surface shape using total absolute curvature. Comput. Graph. 2019, 78, 108–115. [Google Scholar] [CrossRef]
  43. Magid, E.; Soldea, O.; Rivlin, E. A comparison of Gaussian and mean curvature estimation methods on triangular meshes of range image data. Comput. Vis. Image Underst. 2007, 107, 139–159. [Google Scholar] [CrossRef]
  44. Shannon, C.E.; Weaver, W. The Mathematical Theory of Communication; The University of Illinois Press: Urbana, IL, USA, 1964. [Google Scholar]
  45. Mata, P.M.; Kristensen, S.; Shea, K. Perception into a Tool for Three-Dimensional Shape Generation Using a Shape Grammar and a Parametric Model. J. Mech. Des. 2019, 141, 1–11. [Google Scholar] [CrossRef] [Green Version]
  46. Free3D. Available online: https://free3d.com/ja/ (accessed on 28 February 2023).
  47. cgtrader. Available online: https://www.cgtrader.com/ (accessed on 28 February 2023).
  48. Owen, J.S. A Survey of Unstructured Mesh Generation Technology, In Proceedings of the 7th International Meshing Roundtable, Dearborn, MI, USA, 26–28 October 1998.
  49. Ansys. Available online: https://www.ansys.com/ja-jp (accessed on 28 February 2023).
  50. Kato, T.; Matsumoto, T. Morphological evaluation of Closed Planar Curves and Its Application to Aesthetic Evaluation. Graph. Model. 2020, 109, 101064. [Google Scholar] [CrossRef]
  51. Cao, W.; Ping, H.; Yujie, L.; Ming, G.; Hua, L. Gaussian-curvature-derived Invariants for Isometry. Sci. China Inf. Sci. 2013, 56, 1–12. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Extraction procedure.
Figure 1. Extraction procedure.
Entropy 25 00624 g001
Figure 2. Classification of papers.
Figure 2. Classification of papers.
Entropy 25 00624 g002
Figure 3. Principal curvature at point p.
Figure 3. Principal curvature at point p.
Entropy 25 00624 g003
Figure 4. Calculation of curvatures.
Figure 4. Calculation of curvatures.
Entropy 25 00624 g004
Figure 5. Polygons around vertex v i .
Figure 5. Polygons around vertex v i .
Entropy 25 00624 g005
Figure 6. Sample shapes: (a) Extruded shapes; (b) Rotated shapes; (c) Vase shapes.
Figure 6. Sample shapes: (a) Extruded shapes; (b) Rotated shapes; (c) Vase shapes.
Entropy 25 00624 g006
Figure 7. Relationship between the number of states and determination coefficient in extruded shapes.
Figure 7. Relationship between the number of states and determination coefficient in extruded shapes.
Entropy 25 00624 g007
Figure 8. Relationship between the number of states and determination coefficient in rotated shapes.
Figure 8. Relationship between the number of states and determination coefficient in rotated shapes.
Entropy 25 00624 g008
Figure 9. Relationship between the number of states and determination coefficient in extruded shapes.
Figure 9. Relationship between the number of states and determination coefficient in extruded shapes.
Entropy 25 00624 g009
Figure 10. Shapes with large entropy changes between the number of states: (a) contour diagram of the Gaussian curvature with scatter plots of the principal curvatures of a rotated shape (V = 6); (b) ditto (V = 5); (c) an extruded shape example whose entropy variation is large between the number of states; (d) a vase shape example whose entropy variation is large between the number of states.
Figure 10. Shapes with large entropy changes between the number of states: (a) contour diagram of the Gaussian curvature with scatter plots of the principal curvatures of a rotated shape (V = 6); (b) ditto (V = 5); (c) an extruded shape example whose entropy variation is large between the number of states; (d) a vase shape example whose entropy variation is large between the number of states.
Entropy 25 00624 g010
Figure 11. Shape with large entropy changes between the number of states: (a) contour diagram of the mean curvature with scatter plots of the principal curvatures of an extruded shape (V = 7); (b) ditto (V = 6); (c) a rotated shape example whose entropy variation is large between the number of states; (d) a vase shape example whose entropy variation is large between the number of states.
Figure 11. Shape with large entropy changes between the number of states: (a) contour diagram of the mean curvature with scatter plots of the principal curvatures of an extruded shape (V = 7); (b) ditto (V = 6); (c) a rotated shape example whose entropy variation is large between the number of states; (d) a vase shape example whose entropy variation is large between the number of states.
Entropy 25 00624 g011
Figure 12. Shape with large entropy changes between curvatures: (a) contour diagram of the Gaussian curvature with scatter plots of the principal curvatures of an extruded shape; (b) contour diagram of the Casorati curvature with scatter plots of the principal curvatures of an extruded shape; (c) an extruded shape example whose curvature variation is large between the Gaussian and Casorati curvatures; (d) a vase shape example whose curvature variation is large between the Gaussian and Casorati curvature.
Figure 12. Shape with large entropy changes between curvatures: (a) contour diagram of the Gaussian curvature with scatter plots of the principal curvatures of an extruded shape; (b) contour diagram of the Casorati curvature with scatter plots of the principal curvatures of an extruded shape; (c) an extruded shape example whose curvature variation is large between the Gaussian and Casorati curvatures; (d) a vase shape example whose curvature variation is large between the Gaussian and Casorati curvature.
Entropy 25 00624 g012
Figure 13. Shape with large entropy changes between curvatures: (a) contour diagram of the mean curvature with scatter plots of the principal curvatures of a vase shape; (b) contour diagram of the Casorati curvature with scatter plots of the principal curvatures of an extruded shape; (c) an extruded shape example whose curvature variation is large between the Gaussian and Casorati curvatures; (d) a rotated shape example whose curvature variation is large between the Gaussian and Casorati curvature.
Figure 13. Shape with large entropy changes between curvatures: (a) contour diagram of the mean curvature with scatter plots of the principal curvatures of a vase shape; (b) contour diagram of the Casorati curvature with scatter plots of the principal curvatures of an extruded shape; (c) an extruded shape example whose curvature variation is large between the Gaussian and Casorati curvatures; (d) a rotated shape example whose curvature variation is large between the Gaussian and Casorati curvature.
Entropy 25 00624 g013
Figure 14. Shape with large difference between the sensory evaluation value and the curvature moment (standard deviation).
Figure 14. Shape with large difference between the sensory evaluation value and the curvature moment (standard deviation).
Entropy 25 00624 g014
Table 1. Keywords and search criteria.
Table 1. Keywords and search criteria.
CurvatureAND/OR3DAND/ORShapeAND/ORFeature
“Curvature”AND“3d”AND“Shape”AND“Preference”
OROROR
“Three-dimensions”“Form”“Aesthetic”
OROR
“Appearance”“Likeness”
OROR
“Outline”“Liking”
OROR
“Silhouette”“Complex”
OR
“Complexity”
OR
“Novelty”
OR
“Similarity”
OR
“Order”
Table 3. Classification of curvatures.
Table 3. Classification of curvatures.
CurvatureCalculation MethodPropertiesReferences
Mean curvature H = k 1 + k 2 2
  • Identify convex and concave surfaces
[29,37,40]
Gaussian curvature K = k 1 k 2
  • Identify saddle surface as negative values
  • Calculate 0 even on cylindrical surfaces
[29,31,38,39,40,41,42]
Casorati curvature C = k 1 2 + k 2 2 2
  • Identify saddle and uneven surfaces
  • Cannot be calculated on a flat surface
[28]
Shape index S I = 1 2 1 π arctan ( k 1 + k 2 k 1 k 2 )
  • Expressed in absolute value
  • Cannot be calculated on a flat surface
[30,32,33,34]
Curvature index C I = 2 π log ( k 1 2 + k 2 2 2 )
  • Identify flat and curved surface
  • Expressed in absolute value
[30]
Curvature tensor τ v = 1 G A N e θ · L G B · e ¯ · e ¯ t
  • Identify flat and curved surface
  • Appropriate area setting is necessary
[35]
Cone curvature α j = sign ( F j ) | π 2 1 t j i = 1 t j r i j |
  • Identify flat and curved surface
  • Smoothly calculated
[25,26,27,36]
G B : circular region around a vertex; G A   : area of G B ; N e : number of edges; θ : angle between two triangles; L G B : length of the part of the edge inside G B ; e : edge of triangles; e ¯ : unit vector of e ; e ¯ t : transpose vector of e ; F j : jth modeling wave sign ( F j ) : sign function depends on convex–concave at F j ; t j : number of vertices at F j ; r i j : angle of N i N C j .
Table 4. Classification of curvatures.
Table 4. Classification of curvatures.
Information ConsideredAdvantages/DisadvantagesReferences
Curvatures of vertices
  • Feature descriptors can be calculated from curvatures with a small amount of computation.
  • Comparisons between shapes require matching the number of vertices and the order in which curvatures are stored in the array.
[33,35]
Curvatures of vertices and their surroundings
  • Using the curvature of a vertex and its surrounding vertices, this feature descriptor takes the surface features of adjacent vertices into account.
  • This consideration of surface features increases the computational complexity because the curvature must be calculated for multiple vertices.
[24,25,26,31,34,36]
Occurrence probabilities of curvatures
  • Using the distribution of curvatures that appear, the diversity of curvatures is considered.
  • The curvature must be discretized by setting appropriate parameters in advance.
[27,28,30,32,34,37,38,39,40,41]
Transition probabilities of curvatures
  • Using the curvature of two vertices that are adjacent or arbitrarily far apart, surface features between the vertices are considered.
  • The curvature must be discretized by setting appropriate parameters in advance.
[42]
Clusters of curvatures
  • By classifying curvatures based on the similarity of values and using only some groups of curvatures as feature descriptors, shape features are expressed in a short array.
  • Because only some vertices and curvatures are used, detailed features of the surface might not be represented.
[29,30,34,36,37]
Table 5. Experimental results for extruded shapes.
Table 5. Experimental results for extruded shapes.
CurvatureParameterDetermination Coefficient
VEmaxEminMaximumAverageStandard DeviationSN Ratio
Gaussian curvature181.943−0.9490.6660.5880.068−0.804
Mean curvature71.734−0.8940.7040.5910.079−4.826
Casorati curvature101.4770.0000.6520.6100.041−4.350
Shape index81.0000.0000.0.1290.0760.030−26.167
Curvature index50.737−11.7850.0220.0130.016−90.085
Table 6. Experimental results for rotated shapes.
Table 6. Experimental results for rotated shapes.
CurvatureParameterDetermination Coefficient
VEmaxEminMaximumAverageStandard DeviationSN Ratio
Gaussian curvature190.959−0.4280.5030.3870.086−9.144
Mean curvature71.245−1.1880.5510.3540.144−14.918
Casorati curvature91.2000.0000.5010.4210.078−8.297
Shape index71.0000.0000.0050.0010.002−82.639
Curvature index90.340−11.2150.0880.0420.030−61.944
Table 7. Experimental results for vase shapes.
Table 7. Experimental results for vase shapes.
CurvatureParameterDetermination Coefficient
VEmaxEminMaximumAverageStandard DeviationSN Ratio
Gaussian curvature1417.652−4.9350.3510.1700.087−18.629
Mean curvature196.364−3.7910.4710.3210.070−10.731
Casorati curvature66.7990.0000.3980.3550.024−9.054
Shape index71.0000.0000.3560.2900.038−11.093
Curvature index131.441−11.2140.0660.0360.020−67.962
Table 8. Comparison between entropy of occurrence probability and curvature moment.
Table 8. Comparison between entropy of occurrence probability and curvature moment.
Feature DescriptorAverageStandard DeviationSN Ratio
MomentMean0.1420.307−44.786
Standard deviation0.3070.247−38.357
Coefficient of variation0.2100.183−43.974
Kurtosis0.1650.241−50.722
Skewness0.1530.211−58.771
Entropy of occurrence probability0.3930.222−34.310
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

Matsuyama, K.; Shimizu, T.; Kato, T. Systematic Classification of Curvature and Feature Descriptor of 3D Shape and Its Application to “Complexity” Quantification Methods. Entropy 2023, 25, 624. https://doi.org/10.3390/e25040624

AMA Style

Matsuyama K, Shimizu T, Kato T. Systematic Classification of Curvature and Feature Descriptor of 3D Shape and Its Application to “Complexity” Quantification Methods. Entropy. 2023; 25(4):624. https://doi.org/10.3390/e25040624

Chicago/Turabian Style

Matsuyama, Kazuma, Takahiro Shimizu, and Takeo Kato. 2023. "Systematic Classification of Curvature and Feature Descriptor of 3D Shape and Its Application to “Complexity” Quantification Methods" Entropy 25, no. 4: 624. https://doi.org/10.3390/e25040624

APA Style

Matsuyama, K., Shimizu, T., & Kato, T. (2023). Systematic Classification of Curvature and Feature Descriptor of 3D Shape and Its Application to “Complexity” Quantification Methods. Entropy, 25(4), 624. https://doi.org/10.3390/e25040624

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