A computer implemented method for selecting colors comprising: associating, in dependence on a mathematical equation derived from measuring psychophysical responses of a plurality of test subjects to a first plurality of test colors, a color score with each of a second plurality of colors that are numerically defined in a color library; displaying on a visual user interface a combination of at least three colors from the color library, wherein the color score includes at least a color emotion score in said combination; displaying on the visual user interface a combined color score for said color combination; wherein the combined color score is calculated for inviting-uninviting color emotion or color harmony using a controller as follows: (a) calculating for each unique pair of the color combination that are adjacent to each other an adjacent pair color score using a first equation in dependence on the color score of each color in the adjacent pair; (b) calculating for each unique pair of the color combination that are non-adjacent to each other a non-adjacent pair color score using a second equation in dependence of the color score of each color in the non-adjacent pair; (c) wherein the combined color score = 1/n • (Σadjacent pair color score(s) + ∑non-adjacent pair color score(s)), where n is the total number of unique adjacent and non-adjacent pairs in the color combination; wherein the color emotion scores comprise at least one of the inviting-uninviting, exciting-calming, warm-cool, light-dark, clean-dirty, happy-sad, or fun-serious color emotion score, according to one of the following equations: Inviting = − 1.6 + 0.009 L ∗ − 13 2 + a ∗ − 3 0.19 2 + b ∗ − 19 0.29 2 1 2 + 0.015 L ∗; Exciting = − 1.13 + 0.02 L ∗ − 40 2 + a ∗ 0.37 2 + b ∗ 0.60 2 1 2; Warm = 0.02 + 0.12 C ∗ ab 0.6 cos h ab − 45 ° − 1.5 × 10 − 8 L ∗ 3.9; Dark = − 2.1 + 0.05 100 − L ∗; Clean = − 1.8 + 0.043 L ∗ − 28 2 + a ∗ − 3 1.18 2 + b ∗ − 27 1.27 2 1 2; Happy = − 1.3 + 0.025 L ∗ − 14 2 + a ∗ 0.73 2 + b ∗ − 12 0.85 2 1 2; or Fun = − 1.8 + 0.027 L ∗ − 10 2 + a ∗ 0.45 2 + b ∗ − 10 0.61 2 1 2, where L*, a * and b * represent the three CIELAB coordinates of lightness, red/green coordinate and yellow/blue coordinate for the color, respectively, and where C * ab, h ab and L * represent CIELAB chroma, hue angle and lightness, respectively. associating, in dependence on a mathematical equation derived from measuring psychophysical responses of a plurality of test subjects to a first plurality of test colors, a color score with each of a second plurality of colors that are numerically defined in a color library; displaying on a visual user interface a combination of at least three colors from the color library, wherein the color score includes at least a color emotion score in said combination; displaying on the visual user interface a combined color score for said color combination; wherein the combined color score is calculated for inviting-uninviting color emotion or color harmony using a controller as follows: (a) calculating for each unique pair of the color combination that are adjacent to each other an adjacent pair color score using a first equation in dependence on the color score of each color in the adjacent pair; (b) calculating for each unique pair of the color combination that are non-adjacent to each other a non-adjacent pair color score using a second equation in dependence of the color score of each color in the non-adjacent pair; (c) wherein the combined color score = 1/n • (Σadjacent pair color score(s) + ∑non-adjacent pair color score(s)), where n is the total number of unique adjacent and non-adjacent pairs in the color combination; (a) calculating for each unique pair of the color combination that are adjacent to each other an adjacent pair color score using a first equation in dependence on the color score of each color in the adjacent pair; (b) calculating for each unique pair of the color combination that are non-adjacent to each other a non-adjacent pair color score using a second equation in dependence of the color score of each color in the non-adjacent pair; (c) wherein the combined color score = 1/n • (Σadjacent pair color score(s) + ∑non-adjacent pair color score(s)), where n is the total number of unique adjacent and non-adjacent pairs in the color combination; wherein the color emotion scores comprise at least one of the inviting-uninviting, exciting-calming, warm-cool, light-dark, clean-dirty, happy-sad, or fun-serious color emotion score, according to one of the following equations: Inviting = − 1.6 + 0.009 L ∗ − 13 2 + a ∗ − 3 0.19 2 + b ∗ − 19 0.29 2 1 2 + 0.015 L ∗; Exciting = − 1.13 + 0.02 L ∗ − 40 2 + a ∗ 0.37 2 + b ∗ 0.60 2 1 2; Warm = 0.02 + 0.12 C ∗ ab 0.6 cos h ab − 45 ° − 1.5 × 10 − 8 L ∗ 3.9; Dark = − 2.1 + 0.05 100 − L ∗; Clean = − 1.8 + 0.043 L ∗ − 28 2 + a ∗ − 3 1.18 2 + b ∗ − 27 1.27 2 1 2; Happy = − 1.3 + 0.025 L ∗ − 14 2 + a ∗ 0.73 2 + b ∗ − 12 0.85 2 1 2; or Fun = − 1.8 + 0.027 L ∗ − 10 2 + a ∗ 0.45 2 + b ∗ − 10 0.61 2 1 2, where L*, a * and b * represent the three CIELAB coordinates of lightness, red/green coordinate and yellow/blue coordinate for the color, respectively, and where C * ab, h ab and L * represent CIELAB chroma, hue angle and lightness, respectively.