Color Image Evaluation for Small Space Based on FA and GEP

Aiming at the problem that color image is difficult to quantify, this paper proposes an evaluation method of color image for small space based on factor analysis (FA) and gene expression programming (GEP) and constructs a correlation model between color image factors and comprehensive color image. The basic color samples of small space and color images are evaluated by semantic differential method (SDmethod), color image factors are selected via dimension reduction in FA, factor score function is established, and by combining the entropy weight method to determine each factor weights then the comprehensive color image score is calculated finally.The best fitting function between color image factors and comprehensive color image is obtained by GEP algorithm, which can predict the users’ color image values. A color image evaluation system for small space is developed based on thismodel.The color evaluation of a control roomonAC frequency conversion rig is taken as an example, verifying the effectiveness of the proposed method. It also can assist the designers in other color designs and provide a fast evaluation tool for testing users’ color image.


Introduction
Rerferring to the topic of color image, color image processing has to be reviewed.It has become an important research field in computer graphics and image processing.For instance, Wang et al. proposed an effective color image retrieval method based on texture [1] and later presented a novel method for image retrieval based on structure elements' descriptor [2].The method of image authentication [3] and restoration [4] was also been put forward.The results are plentiful and of guiding significance.However, in this paper, color image would be discussed from another perspective of Kansei engineering.Known from color ergonomics research, color affects people's emotions, feelings, thoughts, and behavior unconsciously, either direct visual stimulation or indirect association.As a kind of important visual language, different color produces different stimulus to people's physical and mental; moreover psychological effect of color affects person's mood and behavior [5].In traditional color design, usually according to the designers' subjective experience and a lot of trials and errors to make a decision, now the color design tends to satisfy users' preferences, emphasizing feeling and experience of users.Color image is a kind of color characteristics that consider color property and color psychology together.In color design, designers must study the emotional reactions, such as the feeling associated with emotional cognitive in users' heart, and find out tacit knowledge of users [6].
As one kind of association produced by users' sense, color image relates to factors such as users' vision, experience of life, and cultural background, belongs to the category of mental activity, and can be quantitatively studied.At present, many methods that such as grey theory, fuzzy theory [7], neural network [8], genetic algorithm [9,10], and Kansei engineering [11,12] have been applied to product's form and color design.They attempted to quantify people's psychological feelings of color [13] and acquired consumers' color image finally.Hsiao and Tasi [14] used grey theory and back-propagation nerve network for children walker's color quantitative evaluation research.Hsiao et al. [15] put forward a kind of computer aided color selection system, based on the color harmony, fuzzy set theory, and so on, which could assist consumer to choose the suitable clothing color matching according to skin color and image preferences.Tsai et al. [16] proposed a kind of fast concept design method to predict  image evaluation based on fuzzy neural network and grey theory and applied it in electronic lock design.Tsai and Chou [17] applied grey theory and genetic algorithm to construct an automatic design support system, which could generate color matching scheme of two colors according to the image or evaluate color scheme automatically, and it was used to evaluate the color image of kettle.Ma et al. [18] constructed a design decision-making support model for custom color combinations, which would be applied to the sofa color design.
The above research provides an efficient image evaluation criterion and auxiliary support for product color design.However, these studies mainly focus on product color design, although there are a few studies about color image and preferences of interior space [19] that hardly consider color image problem of small space.Color image is a comprehensive psychological reaction for color that is produced by a variety of sensory organs.The color of small space is different from general product appearance.Suppose a person is surrounded by environmental color, the color of small space will have strong influence on people's perception.How to validate the fact that the color image designed by designers matched with users' expectations?Users have numerous color images, and the image is fuzzy and uncertain; how to quickly and accurately reflect the nonlinear relationship between images?How to change the perceptual color design to rational calculation and then return to the understanding of perceptual once again after design?However, very few studies have discussed these issues.
Therefore, the study aims to quantitatively study the psychological attribute and emotional characteristics of color and solve the problem that users' color image perception is fuzzy.From the perspective of cognitive psychology [20], the methods such as Kansei engineering are used to make the users' implicit color preference image externalization.Aiming at the complicated nonlinear relation between the comprehensive color image and color factors, meanwhile, and considering their intrinsic relevance, excellent GEP algorithm [21] was applied to construct the relationship model between color image and color scheme, and a color optimization design method for small space was built.Finally color image evaluation system for small space was developed based on the model and assisted designers to conduct color evaluation and optimization by users' preference.

Problem Description and Outline of the Method
2.1.Problem Description.Working space [22] refers to the needed activity space when people are operating machine, plus the total space, which is occupied by machinery, equipment, appliances, tools, processing objects, and so forth.
Combining with the definition of working space, the concept of small space is defined.It is a closed or semiclosed space; the range of activity for person is narrow and easily leads to oppression and other uncomfortable feelings.Such as aircraft cockpit, manned submersibles, engineering machinery cab, nuclear power plant control room, oil rig driller control room, and family car cab, in which the persons operate in sitting or standing posture.The space is restricted by various factors; besides the space of machinery and equipment layout, the activity space for person is very limited.The appropriate color design for small space will be able to adjust the space environment's atmosphere, color temperature, space size, and so on.Under the background of customization, designers need to conduct special custom to satisfy the demands of users.
Color is a sense of perception formed by the brain's visual cortex transmitted from nerve impulse, generated by light stimulate eye.The formation of users' perceptual image is a process from the divergent thinking to convergent thinking [23].Color information seen directly and color aesthetic psychology are comprehensively processed by users; the most similar images close to the space color are searched in the memory, and users will combine several subimages to form an overall image.On the contrary, the deconstruction of the users' perceptual image is the inverse process.In order to get the users' overall image, we need to decompose and know its subimages.Because of the users' different life background, experience, culture, preferences, age, gender, and so on, perceptual images produced by different users will form a large number of subimage sets.And the formation of perceptions often accompanied by complicated mental process leads to the fact that the subjective information of images is fuzzy and uncertain.So, in order to build the mapping relationship between the perceptual information and the color design of small space, according to the users' perceptual cognitive structure of thinking (shown in Figure 1), via FA to reduce dimensionality and form the users' perceptual image space, this paper extracted users' implicit knowledge, explored the Step 1 Step 2 Step 3 Step 4 relationship between the color factors and the overall image, and thus mastered users' color psychology.

Outline of the Method.
As shown in Figure 2, the implementation procedures of evaluation method mainly comprise the following steps.
Step 1 (choose the basic color samples of small space and color images).The test color samples were created by regularly adjusting the RGB parameters with constant units within the range of 0-255.Human, material, and financial resources were comprehensively considered; finally, the simulation experiment method of drawing space in plane was adopted.
The extensive color images were acquired by open-ended questionnaire.
Step 2 (analyze the color images qualitatively and quantitatively).The color images were qualitatively analyzed by designers and the representative adjectives were selected.Then, the perceptual image evaluation experiment was carried on by SD method [24].Through FA to reduce dimensionality of subimage sets, factor score function was established, and the weight of each factor was determined by the entropy weight method; finally, the comprehensive perceptual image score was calculated.
Step 3 (establish the prediction model of color image based on GEP).By GEP algorithm's strong ability of function finding, the best fitting function between the color factors and comprehensive perceptual image was obtained.The function could predict the users' color image value and grasp the users' psychology.
Step 4 (construct the color image evaluation system for small space).Based on above-mentioned model, the color image evaluation system for small space was developed, assisted designers to design the color, and provided effective tools to test users' color image.

Prediction Model of Color Image for Small
Space Based on FA and GEP where   ( = 1, 2, . . ., ;  = 1, 2, . . ., ) represents factor loading and   ( = 1, 2, . . ., ) represents special factor.The factor loading is deduced with the method of principal component analysis (PCA).By means of coordinate transformation, after standardization, the original relevant variables  1 ,  2 ,  3 , . . .,   were linearly combined, converted into another set of uncorrelated variables   : where  1 ,  2 , . . .,   , respectively, represent the first, second,. . ., th principal component of the original variables. 1 has the largest proportion in total variance, so it has the strongest ability to synthesize the information of original variables.The rest of principal components' ability to synthesize information declines in turn.FA can reduce the variable dimension and eliminate the correlation between the variables, so as to reduce the input of GEP independent variable in the late.

Entropy Weight Method.
German physicist R. Clausius put forward entropy and applied it in thermodynamics.C. E. Shannon introduced the concept of entropy in information theory later, as a measurement of uncertainty or information of random events [26].In information theory, entropy is used to measure the disorder degree of system.When the system is more orderly, the information entropy is smaller; oppositely, the information entropy is larger.The system has a variety of different state and the probability of each state to appear with the   ( = 1, 2, . . ., ), and then the entropy of the system is defined as When   = 1/ ( = 1, 2, . . ., ), that is to say, the presence of the probability to each state is equal, the maximum entropy is  max = ln .
This paper introduced the information entropy to determine the weights in the color image evaluation for small space.If there are  color samples of small space and  color factors, the data matrix of original evaluation  = (  ) × is as follows: A color factor   has information entropy   = − ∑  =1   ⋅ ln   and   =   / ∑  =1   ( = 1, 2, . . ., ;  = 1, 2, . . ., ).Calculation steps of entropy weights are as follows.
Step 1. Calculate the th color factor of the th color sample of small space's index weight   :   =   / ∑  =1   .
Step 3. Calculate the deviation degree of the th color factor (difference coefficient)   :   = 1 −   .

GEP Algorithm.
GEP is a new random search and optimization algorithm, invented and put forward by the Portuguese scientist Ferreira [27].GEP combines the advantages of genetic algorithm and genetic programming algorithm.It is a new kind of efficient and ideal data mining technique, taking the biotic genetic gene expression patterns as reference.Using adaptive random search method, and without relying on any prior knowledge, GEP can discover formula which can describe data inherent law from the data and shows strong universality and accuracy.In GEP, a computer program is coded into fixed length of linear symbol strings (chromosome), while in the individual fitness calculation, it will be expressed to different shapes and sizes of expression trees, so as to solve the problem.As shown in Figure 2, the main steps of using GEP for function mining are as follows [28].
Step 1. Code the individual and create the initial population.
Population contains a number of individuals (chromosomes).
Chromosome is comprised of more than one gene, and the genes are connected through linking function.Gene is comprised of head and tail.Head consists of function set or terminal set, and end consists of terminal set.
Step 2. Calculate the fitness.The fitness of each individual is calculated through fitness function.Fitness reflects the extent of excellence of individual to achieve the optimal solution in the course of evolution.If the optimal individual satisfies end condition, it will be transferred to the output.If not, transferred to genetic operation steps and then produces offspring with new characteristics.
Step 3. Do a series of genetic operation to produce a new generation of population, using principle of survival of the fittest to guide the evolution, including (1) keep the best individual, (2) selection, (3) replication, (4) mutation, (5) transposition, and (6) recombination.

Prediction Model of Color
Image for Small Space.The core technology of GEP is that the mutation process is completely separated from the evaluation process.Fixed linear string is used in mutation process, and expression tree is used in evaluation process.First of all, determine the symbols representing chromosome; namely, select function set and terminal set which suited the solution.GEP = (, ).Suppose function set  = {+, −, * , /, sin, cos} and common function include arithmetic operators, elementary mathematics function, Boolean operator, and relational operators; terminal set  = { 1 ,  2 ,  3 , . . .,   },  1 ,  2 ,  3 , . . .,   represents color image factor set.
Secondly, determine the genetic structure and the length of gene head.According to the complexity of the problem to define the length of gene head ℎ, the length of tail  and the length of head ℎ satisfy the relationship that  = ℎ(−1)+1.In the formula,  represents the maximum number of operation numbers in function set, namely, the number involved in the operation variable.It can guarantee that the genes are legitimate in this way.Suppose that a gene is composed of the elements in {+, −, * , /, sin, cos,  1 ,  2 ,  3 ,  4 ,  5 ,  6 }, so  = 2.If ℎ = 6, then  = 6(2 − 1) + 1 = 7. Randomly generate a legitimate GEP gene as follows: The italicized part represents the tail.The corresponding expression tree is shown in Figure 3.
The length of each gene is fixed in GEP coding.Gene is made up of two parts: the front effective K-expression and the back of filler components.From top to bottom, left to right sequence, the expression trees in Figure 3 are traversed, and the K-expression of expression trees is obtained.The effective length of the gene has 10 characters: sin − 5 * * +  4  1  2  3 .So the mathematical expression is as follows: Again, determine the number of genes in each chromosome and the linking function to connect the genes.Fourth, select the fitness function.There are three kinds of classical fitness function in GEP algorithm.According to the characteristics of the problem and reference [29], the fitness function is defined as follows: In (7),   = (1/) ∑  =1 (  −   ) 2 is the mean square error (MSE) of experimental samples.m is the total number of training set samples;   is the output value of the th sample of the th individual, calculated by the mathematical expressions acquired by GEP algorithm modeling;   is the target value of the th sample.
Finally, determine the genetic control parameters before the algorithm running, including the size of the population, evolution generation limit, and the probability of each genetic operators.

Implementation Procedures
The effectiveness and feasibility of the proposed method are demonstrated by taking the case of color image evaluation of a control room on AC frequency conversion rig for illustration purposes.The steps involved in implementing the evaluation system are described in the following sections.

Basic Color Samples. The amount of basic color samples
should not only reduce the workload of experiment but also not affect the measure conditions to express image.According to [17], to set color intervals, RGB parameters change with 64 units within the range of 0-255.Each primary axis gets 5 intervals; eventually 5 3 = 125 sets of basic color samples are obtained (shown in Figure 4), defined as  = { 1 ,  2 , . . .,   }.

Small Space Samples.
In order to avoid the mutual influence of multivariable, this paper only discusses the main color of small space and does not consider the impact of factors such as material and form for the moment.As shown in Figure 5, the space sample with transparent background (1) Modern ( 6) Gorgeous ( 11) Elegant ( 16) Quiet ( 21) Warm ( 26) Lively ( 2) Dynamic ( 7) Noble ( 12) Graceful ( 17) Solemn ( 22) Soft ( 27) Lightful (3) Simple ( 8) Male ( 13) Fresh ( 18) Classical ( 23) Joyful ( 28) Bright ( 4) Technological ( 9) Handsome ( 14) Cool ( 19) Plain ( 24) Romantic ( 29) Roomy ( 5) Fashionable    is made as a master page.In test process, the random superposition of the 125 color samples, which can form the color samples of small space to meet the experimental demand is shown in Figure 6.

Qualitative Analysis.
Through open-ended questionnaire [30] to acquire extensive color images, they were qualitatively analyzed by designers, according to their practical work experience.The repeated meaning words were deleted; and the similar meaning words were merged.Ultimately the selected representative adjectives covered color semantic cognitive space as much as possible.29 color images are preliminary selected (shown in Table 1).

Quantitative
Analysis.40 subjects (half male and half female) are invited to take part in the experiment.All of them are students, aged between 20 and 28 years old, and both eyesight and color vision are normal.As shown in Figures 7 and 5 points psychology scale is used to represent the different dimensions of psychological variation.The evaluation value of color image for small space is set between 0 and 1; the larger value represents the more close to the feeling of color image.The color samples of small space randomly are displayed on the computer screen and the subjects evaluate them according to the 29 color images, respectively.The experimental results are analyzed with SPSS Statistics19.0software for FA.First by calculating the correlation coefficient matrix, anti-image correlation matrix, Bartlett's test of sphericity, and Kaiser-Meyer-Olkin (KMO) test, the relationship between variables is tested.The statistics observed value is 5456.468 in Bartlett's test of sphericity; since the corresponding probability of  value is close to 0, less than the significance level of  ( equal to 0.05), it can be regarded as that there is significant difference between the correlation coefficient matrix and unit matrix.The value of KMO is 0.886; according to the KMO metrics provided by Kaiser, the original variables are suitable to conduct FA.
According to the principle [25] that "usually select the number of eigenvalues as factor number when cumulative variance contribution rate is greater than 0.85, " the factors are extracted by the method of PCA.Six factors of color image are extracted; the corresponding cumulative variance contribution rate reaches 87.289% (shown in Table 2) and meets the above principle.
In Table 3, the data shown in bold in each column represent the color images as having high loading on the 6 factors, respectively.For example, the fifth factor mainly explains four color images: the lightful, bright, lively, and roomy.Therefore, we can summarize some conclusions that the first factor mainly explains the feeling of being modern and fashionable, the second factor mainly explains the feeling of being relaxed and calm, the third factor mainly explains the feeling of being sweet and elegant, the fourth factor mainly   The factor score of each color sample of small space can be calculated by (8), then the method of calculating total score by factor weight is adopted, and the overall image is assessed comprehensively.Determining the weight of the factors is the key.Based on the principle of information entropy and the original evaluation data information, this paper used the differences between the values which reflected "information value" of factor to determine the weight.SD method is used to quantify the color image of subjects' preferences; evaluation data is sorted and described as follows: where  represents the number of color samples of small space;  represents the number of factors;    represents the th subject's preference to the th factor of the th color sample of small space.
According to the above-mentioned calculation steps of entropy weights, the objective weight of each factor can be calculated, denoted as   = ( 1 ,  2 , . . .,   ).Entropy weight method is a kind of objective method.In order to retain experts and policymakers' subjective opinions, namely, balance subjective and objective aspect, the entropy weight method is combined with each factor variance contribution rate in the above-mentioned FA.Suppose the variance contribution rate calculated in FA is subjective weight   and entropy   is the objective weight; the combination weight   of the th factor can be calculated as follows: Table 4 shows the calculation results of weights.The comprehensive color image evaluation of small space is as follows:

Color Image Prediction
Based on GEP.84 groups of data were randomly selected from 125 groups of experimental data as the training set and the rest as a validation set.In (13),    ,   , and  avg represent the predicted value, the actual value, and the actual average value, respectively.The range of -square is [0, 1], the closer to 1, showing that the equation's 6 variables have stronger ability to explain the comprehensive color image ().The calculation result shows that the MSE and -square of training set are 0.0006 and 0.9614, respectively.In order to verify the validity of the GEP model, 41 sets of data are taken into the model, and the curve fitting of validation set is shown in Figure 10.The MSE and -square are 0.0007 and 0.9453, respectively, achieving the ideal effect.4.4.Evaluation System.According to the above-mentioned method, VB is used to develop a color image evaluation system for small space.The system can assist designers in color design.In the system, designers input color value and users' factor score of color scheme.According to the function The calculation and analysis results of application example showed that the prediction result was close to the target and it had high prediction accuracy.The example of color image evaluation of driller control room proved that GEP had strong nonlinear and global search ability to find function and provided a solution way for rapid evaluation of color image.In the next step of research work, different algorithms will be used to compare the results and analysis.At the same time, the color collocation and the coupling relationship between color and shape will be the research content in next step.

Figure 1 :
Figure 1: The mapping of users' perceptual cognitive.

Figure 2 :
Figure 2: The color image evaluation method for small space.

Figure 5 :
Figure 5: Master page of the space sample.

Figure 6 :Figure 7 :
Figure 6: An example of the color sample of small space.

Figure 9 2 ∑Figure 9 :
Figure 9 shows the curve fitting of GEP algorithm in the training set.MSE and -square are used to verify the validity of the algorithm and the ability of the prediction.The formulas are as follows: MSE = 1   ∑ =1

Figure 10 :
Figure 10: The fitting curve of validation set.
[25]5]is a method of multivariate statistical analysis, which studies how to make numerous original variables condense into a few factors with the least amount of information loss and how to make the factors have certain named explanatory.Suppose that there are  original variables  1 ,  2 ,  3 , . . .,   , after standardized treatment, variables are expressed by  factors ( < ), by means of linear combination of  1 ,  2 ,  3 , . . .,   : 3.1.Theoretical Background 3.1.1.FA.

Table 1 :
Color images obtained from questionnaire.

Table 4 :
The weight of factors.

Table 5
shows the main operation parameters in GEP model.In function set, Sqrt represents square root; 2 represents square; Avg. 2 represents the mean of two variables;

Table 5 :
GEP parameter set.,  2 ,  3 ,  4 ,  5 ,  6 }  1 ,  2 ,  3 ,  4 ,  5 ,  6 , respectively, represent the 6 factors obtained from FA.The visual basic programming is used to express the GEP algorithm, and the program runs multiple times to get the best individual.The best individual's expression tree is shown in Figure8.Each expression tree represents a gene, and genes are connected by linking function "+" to form a chromosome.The best individual translated into mathematical expression is as follows: 2⋅  4 .