Multimodal Prediction and Personalization of Photo Editswith Deep Generative Models

Ardavan SaeediCSAIL, MIT

Matthew D. HoffmanAdobe Research

Stephen J. DiVerdiAdobe Research

Asma GhandehariounMedia Lab, MIT

Matthew J. JohnsonGoogle Brain

Ryan P. AdamsHarvard and Google Brain

Abstract

Professional-grade software applications arepowerful but complicated—expert users canachieve impressive results, but novices oftenstruggle to complete even basic tasks. Photoediting is a prime example: after loading aphoto, the user is confronted with an arrayof cryptic sliders like “clarity”, “temp”, and“highlights.” An automatically generated sug-gestion could help, but there is no single “cor-rect” edit for a given image—different expertsmay make very different aesthetic decisionswhen faced with the same image, and a sin-gle expert may make different choices depend-ing on the intended use of the image (or ona whim). We therefore want a system thatcan propose multiple diverse, high-quality ed-its while also learning from and adapting toa user’s aesthetic preferences. In this work,we develop a statistical model that meets theseobjectives. Our model builds on recent ad-vances in neural network generative modelingand scalable inference, and uses hierarchicalstructure to learn editing patterns across manydiverse users. Empirically, we find that ourmodel outperforms other approaches on thischallenging multimodal prediction task.

1 Introduction

Many office workers spend most of their working daysusing pro-oriented software applications. These applica-tions are often powerful, but complicated. This complex-ity may overwhelm and confuse novice users, and evenexpert users may find some tasks time-consuming andrepetitive. We want to use machine learning and statisti-cal modeling to help users manage this complexity.

Fortunately, modern software applications collect largeamounts of data from users with the aim of providingthem with better guidance and more personalized ex-periences. A photo-editing application, for example,could use data about how users edit images to learn whatkinds of adjustments are appropriate for what images,and could learn to tailor its suggestions to the aestheticpreferences of individual users. Such suggestions canhelp both experts and novices: experts can use them asa starting point, speeding up tedious parts of the editingprocess, and novices can quickly get results they couldnot have otherwise achieved.

Several models have been proposed for predicting andpersonalizing user interaction in different software ap-plications.

These existing models are limited in that they only pro-pose a single prediction or are not readily personalized.Multimodal predictions1 are important in cases where,given an input from the user, there could be multiple pos-sible suggestions from the application. For instance, inphoto editing/enhancement, a user might want to applydifferent kinds of edits to the same photo depending onthe effect he or she wants to achieve. A model shouldtherefore be able to recommend multiple enhancementsthat cover a diverse range of styles.

In this paper, we introduce a framework for multimodalprediction and personalization in software applications.We focus on photo-enhancement applications, thoughour framework is also applicable to other domains wheremultimodal prediction and personalization is valuable.Figure 1 demonstrates our high-level goals: we want tolearn to propose diverse, high-quality edits, and we wantto be able to personalize those proposals based on users’historical behavior.

1We mean “multimodal” in the statistical sense (i.e., comingfrom a distribution with multiple maxima), rather than in thehuman-computer-interaction sense (i.e., having multiple modesof input or output).

arX

iv:1

704.

0499

7v1

[st

at.M

L]

17

Apr

201

7

Our modeling and inference approach is based on thevariational autoencoder (VAE) Kingma and Welling[2013] and a recent extension of it, the structured varia-tional autoencoder (SVAE) Johnson et al. [2016]. Alongwith our new models, we develop approximate inferencearchitectures that are adapted to our model structures.

We apply our framework to three different datasets (col-lected from novice, semi-expert, and expert users) of im-age features and user edits from a photo-enhancementapplication and compare its performance qualitativelyand quantitatively to various baselines. We demonstratethat our model outperforms other approaches.

2 Background and related work

In this section, we first briefly review the frameworks(VAEs and SVAEs) that our model is built upon; next,we provide an overview of the available models for pre-dicting photo edits and summarize their pros and cons.

2.1 Variational autoencoder (VAE)

The VAE, introduced in Kingma and Welling [2013], hasbeen successfully applied to various models with contin-uous latent variables and a complicated likelihood func-tion (e.g., a neural network with nonlinear hidden lay-ers). In these settings, posterior inference is typicallyintractable, and even approximate inference may be pro-hibitively expensive to run in the inner loop of a learn-ing algorithm. The VAE allows this difficult inferenceto be amortized over many learning updates, makingeach learning update cheap even with complex likelihoodmodels.

As an instance of such models, consider modeling aset of N i.i.d. observations y = {yn}Nn=1 with the fol-lowing generative process: zn

iid∼ h and yn ∼ f(gθ(zn)),where zn is a latent variable generated from aprior h(z) (e.g., N (0, I )) and the likelihood func-tion pθ(yn|zn) = f(yn; gθ(zn)) is a simple distribu-tion f whose parameters gθ(zn) can be a compli-cated function of zn. For example, pθ(yn|zn) mightbeN (yn;µ(zn; θ),Σ(zn; θ)) where the mean and the co-variance depend on zn through a multi-layer perceptron(MLP) richly parameterized by weights and biases θ. SeeFigure 2(a) for the graphical model representation of thisgenerative process.

In the VAE framework, the posterior density pθ(z|y) isapproximated by a recognition network qφ(z|y), whichcan take the form of a flexible conditional density modelsuch as an MLP parameterized by φ. To learn the pa-rameters of the likelihood function θ and the recognitionnetwork φ, the following lower bound on the marginal

likelihood is maximized with respect to θ and φ:

LVAE(φ, θ) , Eqφ(z|y)[log pθ(y|z)]− KL(qφ(z|y)||p(z)).(1)

To compute a Monte Carlo estimate of the gradient ofthis objective with respect to φ, Kingma and Welling[2013] propose a reparameterization trick for samplingfrom qφ(z|y) by first sampling from an auxiliary noisevariable and then applying a differentiable map to thesampled noise. This yields a differentiable Monte Carloestimate of the expectation with respect to φ. Given thegradients, the parameters are updated by stochastic gra-dient ascent.

2.2 Structured variational autoencoder (SVAE)

Johnson et al. [2016] extend the VAE inference schemeto latent graphical models with neural network observa-tion distributions. This SVAE framework combines theinterpretability of graphical models with the flexible rep-resentations found by deep learning. For example, con-sider a latent Gaussian mixture model (GMM) with non-linear observations:

π ∼ Dir(α), (µk,Σk) ∼ NIW(λ), sniid∼ π,

zn|sn, {µk,Σk}Kk=1 ∼ N (µsn ,Σsn), (2)yn|zn, θ ∼ N (µ(zn; θ),Σ(zn; θ)) .

Note that the nonlinear observation model for each datapoint yn|zn, θ resembles that of the VAE, while the la-tent variable model for zn is a GMM (see Figure 2(b)).This latent GMM can represent explicit latent cluster as-signments while also capturing complex non-Gaussiancluster shapes in the observations.

To simplify the SVAE notation, we consider a generalsetting in which we denote the global parameters of agraphical model by γ and the local latent variables by z.Furthermore, we assume that p(z|γ) and p(γ) are a con-jugate pair of exponential family densities with sufficientstatistics tz(z) and tγ(γ). We continue to use θ to de-note the parameters of a potentially complex, nonlinearobservation likelihood. Using a mean field family dis-tribution for approximating the posterior the variationallower bound (VLB) can be written as:

L(ηγ , θ, ηz) , Eq(γ)q(z)

[log

p(γ)p(z|γ)pθ(y|z)q(γ)q(z)

](3)

where ηγ and ηz are the parameters of the variationaldistributions q(γ) and q(z) respectively.

Due to the non-conjugate likelihood function pθ(y|z),standard variational inference methods cannot be appliedto the latent GMM. To solve this problem, the SVAE re-places the non-conjugate likelihood with a recognition

User Group I User Group IIOriginal Image

(a) (b)Figure 1: The main goals of our proposed models: (a) Multimodal photo edits: For a given photo, there may be multiple validaesthetic choices that are quite different from one another. (b) User categorization: A synthetic example where different userclusters tend to prefer different slider values. Group 1 users prefer to increase the exposure and temperature for the baby images;group 2 users reduce clarity and saturation for similar images.

model r(y;φ) that generates conjugate evidence poten-tials. We can then define a surrogate objective L̂ withconjugacy structure:

L̂(ηγ , ηz, φ) , (4)

Eq(γ)q(z)

[log

p(γ)p(z|γ) exp{〈r(y;φ), tz(z)〉}q(γ)q(z)

]where the potentials exp{〈r(y;φ), tz(z)〉} havea conjugate form to p(z|γ). By choosing ηzto optimize this surrogate objective, writingη∗z(ηγ , φ) , arg maxηz L̂(ηγ , ηz, φ), the SVAE ob-jective is then LSVAE(ηγ , θ, φ) , L(ηγ , θ, η

∗z(ηγ , φ))

which can be shown to lower bound the variational infer-ence objective in eq. 3. As in the stochastic variationalinference (SVI) algorithm Hoffman et al. [2013], thereis a simple expression for the natural gradient of thisobjective with respect to the variational parameters withconjugate priors; the gradients w.r.t. other variational pa-rameters, such as those parameterizing neural networks,can be computed using the reparameterization trick.

2.3 Related work on the prediction of photo edits

There are two main categories of models, parametricand nonparametric, that have been used for prediction ofphoto edits:

Parametric methods These methods approximate aparametric function by minimizing a squared (or a simi-lar) loss. The loss is typically squared L2 distance in Labcolor space, which more closely approximates humanperception than RGB space [Sharma and Bala, 2002].This loss is reasonable if the goal is to learn from a setof consistent, relatively conservative edits. But when ap-plied to a dataset of more diverse edits, a model that min-imizes squared error will tend to predict the average edit.

At best, this will lead to conservative predictions; in theworst case, the average of several good edits may pro-duce a bad result.

Bychkovsky et al. [2011] collect a dataset of 5000 pho-tos enhanced by 5 different experts; they identify a setof features and learn to predict the user adjustments aftertraining on the collected dataset. They apply a numberof regression techniques such as LASSO and Gaussian-process regression and show their proposed adjustmentscan match the adjustments of one of the 5 experts. Theirmethod only proposes a single adjustment and the per-sonalization scheme that they suggest requires the userto edit a set of selected training photos.

Yan et al. [2016] use a deep neural network to learn amapping from an input photo to an enhanced one fol-lowing a particular style; their results show that the pro-posed model is able to capture the nonlinear and complexnature of this mapping. They also incorporate semanticawareness in their model, so their model can predict theadjustments based on the semantically meaningful ob-jects (e.g., human, animal, sky, etc.) in the photo. Thismethod also only proposes a single style of adjustment.

Jaroensri et al. [2015] propose a technique that can pre-dict an acceptable range of adjustments for a given photo.The authors crowd-sourced a dataset of photos with var-ious brightness and contrast adjustments, and asked theparticipants to label each edited image as “acceptable” or“unacceptable”.

From this labeled dataset they learn a support vector ma-chine classifier that can determine whether an adjust-ment is acceptable or not. They use this model to pre-dict the acceptable range of edits by first sampling fromthe parameter space and then using their learned modelto analyze each sample. Although their model is able to

propose a range of edits to the user, it requires a bal-anced, human-labeled training set of “acceptable” and“unacceptable” images. Since the number of bad editsmay grow exponentially with the dimensionality of theadjustment space, they mostly limit their study to two-dimensional brightness and contrast adjustments.

Nonparametric methods These methods are typicallyable to propose multiple edits or some uncertainty overthe range of adjustments.

Lee et al. [2015] propose a method that can generate adiverse set of edits for an input photograph. The authorshave a curated set of exemplar images in various styles.They use an example-based style-transfer algorithm totransfer the style from an exemplar image to an inputphotograph. To choose the right exemplar image, theydo a semantic similarity search using features that theyhave learned via a convolutional neural network (CNN).Although their approach can recommend multiple editsto a photo, their edits are destructive; that is, the user isnot able to customize the model’s edits.

Koyama et al. [2016] introduce a model for personaliz-ing photo edits only based on the history of edits by asingle user. The authors use a self-reinforcement proce-dure in which after every edit by a user they 1) updatethe distance metric between the user’s past photos 2) up-date a feature vector representation of the user’s photosand 3) update an enhancement preference model basedon the feature vectors and the user’s enhancement his-tory. This model requires data collection from a singleuser and does not benefit from other users’ information.

2.4 Related multimodal prediction models

Traditional neural networks using mean squared error(MSE) loss cannot naturally handle multimodal predic-tion problems, since MSE is minimized by predicting theaverage response. Neal [1992] introduces stochastic la-tent variables to the network and proposes training Sig-moid Belief Networks (SBN) with only binary stochasticvariables. However, this model is difficult to train, andit can only make piecewise-constant predictions and istherefore not a natural fit to continuous-response predic-tion problems.

Bishop [1994] proposes mixture density networks(MDN), which are more suitable for continuous data. In-stead of using stochastic units, the model directly outputsthe parameters of a Gaussian mixture model. That is, asome of the network outputs are used as mixing weightsand the rest provide the means and variances of the mix-ture components. The complexity of MDNs’ predictivedistributions is limited by the number of mixture compo-

nents if the optimal predictive distribution cannot be wellapproximated by a relatively small number of Gaussians,then an MDN may not be an ideal choice.

Tang and Salakhutdinov [2013] add deterministic hid-den variables to SBNs in order to model continuous dis-tributions. The authors showed improvements over theSBN; nevertheless, training the stochastic units remaineda challenge due to the difficulty of doing approximate in-ference on a large number of discrete variables.

Dauphin and Grangier [2015] propose a new class ofstochastic networks called linearizing belief networks(LBN). LBN combines deterministic units with stochas-tic binary units multiplicatively. The model uses deter-ministic linear units which act as multiplicative skip con-nections and allow the gradient to flow without diffusion.The empirical results show that this model can outper-form standard SBNs.

3 Models

Given the limitations of the available methods for pre-dicting photo edits (described in Section 2.3), our goalis to propose a framework in which we can: 1) recom-mend a set of diverse, parametric edits based on a labeleddataset of photos and their enhancements, 2) categorizethe users based on their style and type of edits they apply,and finally 3) personalize the enhancements based on theuser category. We focus on the photo-editing applicationin this paper, but the proposed framework is applicable toother domains where users must make a selection froma large, richly parameterized design space where thereis no single right answer (for example, many audio pro-cessing algorithms have large numbers of user-tunableparameters).

Our framework is based on VAEs and their extensionSVAEs, and follows a mixture-of-experts design [Mur-phy, 2012, Section 11.2.4]. We first introduce a condi-tional VAE that can generate diverse set of enhancementsto a given photo. Next, we extend the model to catego-rize the users based on their adjustment style. Our modelcan provide interpretable clusters of users with similarstyle. Furthermore, the model can provide personalizedsuggestions by first estimating a user’s category and thensuggesting likely enhancements conditioned on that cat-egory.

3.1 Multimodal prediction with conditionalGaussian mixture variational autoencoder(CGM-VAE)

Given a photo, we are interested in predicting a set ofedits. Each photo is represented by a feature vector xn

N

✓yn

zn

N

✓yn

zn

sn

K

µk,⌃k

⇡

↵

�

(a) (b)Figure 2: (a) VAE with Gaussian latent variables (Section 2.1)(b) SVAE with a GMM latent graphical model (Section 2.2)

and its corresponding edits yn are represented by a vec-tor of slider values (e.g. contrast, exposure, saturation,etc.). We assume that there are L clusters of possible ed-its for each image. To generate the sliders yn for a givenimage xn, we first sample a cluster assignment sn anda set of latent features zn from its corresponding mix-ture component N (µsn ,Σsn). Next, conditioned on theimage and zn, we sample the slider values. The overallgenerative process for the slider values {yn}Nn=1 condi-tioned on the input images {xn}Nn=1 is

sn|π iid∼ π, zn|sn,{µ`,Σ`}L`=1 ∼ N (µsn ,Σsn),

yn|xn, zn, θ ∼ N (µ(zn, xn; θ),Σ(zn, xn; θ)), (5)

where µ(zn, xn; θ) and Σ(zn, xn; θ) are flexible para-metric functions, such as MLPs, of the input imagefeatures xn concatenated with the latent features zn.Summing over all possible values for the latent vari-ables sn and zn, the marginal likelihood p(yn|xn) =∑sn

∫znp(yn, sn, zn|xn)dzn yields complex, multi-

modal densities for the image edits yn.

The posterior p(s, z|x, y) is intractable. We approximateit with variational recognition models as

pθ(s, z|x, y) ≈ qφs(s|x, y)qφz (z|x, y, s). (6)

Note that this variational distribution does not model sand z as independent. For qφs(s|x, y), we use an MLPwith a final softmax layer, and for qφz (z|x, y, s), we usea Gaussian whose mean and covariance are the output ofan MLP that takes s, x, and y as input.

Given this generative model and variational family, toperform inference we maximize a variational lowerbound on log pθ(y|x), writing the objective as

L(θ, φ) ,Eqφ(s,z|x,y)[log pθ(y|z, x)]

− KL(qφ(s, z|x, y)||pθ(s, z)).By marginalizing over the latent cluster assignments s,the CGM-VAE objective can be optimized using stochas-tic gradient methods and the reparameterization trick as

in Section 2.1. Marginalizing out the discrete latent vari-ables is not computationally intensive since s and y areconditionally independent given z, pθ(s, z) is cheap tocompute relative to pθ(y|x, z), and we use a relativelysmall number of clusters. However, with a very largediscrete latent space, one could use alternate approachessuch as the Gumbel-Max trick Maddison et al. [2016].

Figure 3 (parts a and b) outlines the graphical modelstructures of the CGM-VAE and its variational distribu-tions qφ.

3.2 Categorization and personalization

In order to categorize the users based on their adjustmentstyle, we extend the basic CGM-VAE model to a hierar-chical model that clusters users based on the edits theymake. While the model in the previous section consid-ered each image-edit pair xn, yn in isolation, we now or-ganize the data according to U distinct users, using xunto denote the nth image of user u and yun to denote thecorresponding slider values. Nu denotes the number ofphotos edited by user u. As before, we assume a GMMwith L components {µ`,Σ`}L`=1 and mixing weights πto model the user categories.

For each user u we sample a cluster index su to indicatethe user’s category, then for each photo n ∈ {1, . . . , Nu}we sample the latent attribute vector zun from the corre-sponding mixture component:

su|π iid∼ π, zun|su, {(µ`,Σ`)}L`=1iid∼ N (µsu ,Σsu).

Finally, we use the latent features zun to generate thevector of suggested slider values yun. As before, we usea multivariate normal distribution with mean and vari-ance generated from an MLP parameterized by θ:

yun|xun, zun, θ iid∼ N (µ(zun,xun; θ),Σ(zun, xun; θ)).

For inference in the CGM-SVAE model, our goal is tomaximize the following VLB:

L , Eq

[log

p(y, z, s|x)

q(s, z)

]. (7)

To optimize this objective, we follow a similar approachto the SVAE inference framework described in Sec-tion 2.2. In the following we define the variational factorsand the recognition networks that we use.

Variational factors For the local variables z and s, werestrict q(z|s) to be normal with natural parameters ηzand we have q(s) in the categorical form with naturalparameter ηs. As in the CGM-VAE, we marginalize overcluster assignments at the user level.

N

✓yn

znxn

sn ⇡

{µ`,⌃`}L`=1

N

yn

znxn

sn

✓

⇡

U

zun

su

xun

yun

{µ`,⌃`}L`=1

Nu

U

zun

su

xun

yunNu

(a) (b) (c) (d)Figure 3: (a) The graphical model for CGM-VAE introduced in Section 3.1 (b) The dependency structure for the variational ap-proximations qφs(s|x, y) and qφz (z|x, y, s) in CGM-VAE (c) The CGM-SVAE model introduced in Section 3.2) for categorizationand personalization. There are U users and each user u has Nu photos. (d) The dependency structure in the variational distributionsfor the CGM-SVAE model. Note that the recognition network for su depends on all the images and their corresponding slidervalues of user u.

For a dataset of U users, the VLB factorizes as follows:

L(θ, φ) ,1

U

∑u

[ Nu∑n=1

Eqφ(z,s|x,y)[log pθ(yun|xun, zun)]

− KL(qφ(zun|xun, yun, su)||pθ(zun|su))]

− KL(qφ(su|{xun, yun}Nun=1)||p(su)).

Figure 3 (parts c and d) outlines the graphical modelstructures of the CGM-SVAE and its variational distri-butions qφ.

To adapt the recognition network r((x, y), φ) used in thelocal inference objective (eq. 4) to our model structure,we write

q(su|{xun, yun}Nun=1;φ) ∝

p(su) exp

{⟨log π+

∑Nun=1 log r(yun, xun;φ), ts(su)

⟩},

where ts(su) denotes the one-hot vector encoding of themixture component index su. That is, for each user im-age xun and corresponding set of slider values yun, therecognition network produces a potential over the user’slatent mixture component su. These image-by-imageguesses are then combined with each other and with theprior to produce the inferred variational factor on su.

This recognition network architecture is both natural andconvenient. It is natural because a powerful enough r canset rk(yun, xun;φ) ∝ pθ(yun|xun, su = k), in whichcase qφ(su|{xun, yun}Nun=1) ≡ p(su|{xun, yun}Nun=1)and there is no approximation error. It is convenientbecause it analyzes image-edit pair independently, andthese evidence potentials are combined in a symmetric,exchangeable way that extends to any number of user im-ages Nu.

4 Experiments

We evaluate our models and several strong baselines onthree datasets. We focus on the photo editing softwareAdobe Lightroom. The datasets that we use cover threedifferent types of users that can be roughly described as1) casual users who do not use the application regularly,2) frequent users who have more familiarity with the ap-plication and use it more frequently 3) experts who havemore experience in editing photos than the other twogroups. We split all three datasets into 10% test, 10%validation, and 80% train set.

Datasets The casual users dataset consists of 345000images along with the slider values that a user has ap-plied to the image in Lightroom. There are 3200 users inthis dataset. Due to privacy concerns, we only have ac-cess to the extracted features from a convolutional neuralnetwork (CNN) applied to the images. Hence, each im-age in the dataset is represented by a 1024-dimensionalvector. For the possible edits to the image, we only focuson 11 basic sliders in Lightroom. Many common editingtasks boil down to adjusting these sliders. The 11 basicsliders have different ranges of values, so we standardizethem to all have a range between−1 and 1 when trainingthe model.

The frequent users dataset contains 45000 images (in theform of CNN features) and their corresponding slidervalues. There are 230 users in this dataset. These usersgenerally apply more changes to their photos comparedto the users in the casual group.

Finally, the expert users dataset (Adobe-MIT5k, col-lected by Bychkovsky et al. [2011]) contains 5000 im-ages and edits applied to these images by 5 different ex-perts, for a total of 25000 edits.

We augment this dataset by creating new images after ap-plying random edits to the original images. To generate

0.0

0.2M

LP

0.0

0.2

MD

N

0.0

0.2

LBN

0.0

0.2

CG

M-

VAE

Temperature0.0

0.2

True

Tint Exposure Contrast Highlights

Figure 4: Marginal statistics for the prediction of the slidersin the casual users dataset (test set). Due to space limitation,We only display the top 5 mostly used sliders in the dataset.LBN has limited success compared to CGM-VAE. MLP mostlyconcentrates around the mean edit. The quantitative compari-son between different methods in terms of the distance betweennormalized histograms is provided in Table 1.

a random edit from a slider, we add uniform noise from arange of±10% of the total range of that slider. Given theaugmented set of images, we extract the “FC7” featuresof a VGG-16 Simonyan and Zisserman [2014] pretrainednetwork and use the 4096-dimensional feature vector asa representation of each image in the dataset. After aug-menting the dataset, we have 15000 images and 75000edits in total. Similar to other datasets, we only focus onthe basic sliders in Adobe Lightroom.

Baselines We compare our model for multimodal pre-diction with several models: a multilayer perceptron(MLP), mixture density network (MDN), and lineariz-ing belief network (LBN). The MLP is trained to predictthe mean and variance of a multivariate Gaussian dis-tribution; this model will demonstrate the limitations ofeven a strong model that makes unimodal predictions.The MDN and LBN, which are specifically designed formultimodal prediction, are other baselines for predictingmultimodal densities. Table 1 summarizes our quantita-tive results.

We use three different evaluation metrics to compare themodels. The first metric is the predictive log-likelihoodcomputed over a held-out test set of different datasets.Another metric is the Jensen-Shannon divergence (JSD)between normalized histograms of marginal statistics ofthe true sliders and the predicted sliders. Figure 4 showssome histograms of these marginal statistics for the ca-sual users.

Finally, we use the mean squared error in the CIE-LABcolor space between the expert-retouched image and themodel-proposed image. We use the CIE-LAB color

space as it is more perpetually linear compared to RGB.We only calculate this error for the experts dataset (testset) since that is the only dataset with available retouchedimages. To compute this metric, we first apply the pre-dicted sliders from the models to the original image andthen convert the generated RGB image to a LAB image.For reference the difference between white and black inCIE-Lab is 100 and photos with no adjustments result inan error of 10.2 . Table 1, shows that our model outper-forms the baselines across all these metrics.

Hyperparameters For the CGM-VAE model, wechoose the dimension of the latent variable from {2, 20}and the number of mixture components from the set {3,5, 10}. For the remaining hyperparameters see the sup-plementary materials.

Tasks In addition to computing the predictive log-likelihood and JSD over the held-out test sets for all threedatasets, we consider the following two tasks:

1. Multimodal prediction: We predict different editsapplied to the same image by the users in the expertsdataset. Our goal is to show that CGM-VAE is ableto capture different styles from the experts.

2. Categorizing the users and adapting the predictionsbased on users’ categories: We show that the CGM-SVAE model, by clustering the users, makes betterpredictions for each user. We also illustrate how in-ferred user clusters differ in terms of edits they ap-ply to similar images.

4.1 Multimodal predictions

To show that the model is capable of multimodal pre-dictions, we propose different edits for a given image inthe test subset of the experts dataset. To generate theseedits, we sample from different cluster components ofour CGM-VAE model trained on the experts dataset. Foreach image we generate 20 different samples and alignthese samples to the experts’ sliders. From the 5 expertsin the dataset, 3 propose a more diverse set of edits com-pared to the others; hence, we only align our results tothose three to show that the model can reasonably cap-ture a diverse set of styles.

For each image in the test set, we compare the predic-tions of MLP, LBN, MDN and the CGM-VAE with theedits from the 3 experts. In MLP (and also MDN), wedraw 20 samples from the Gaussian (mixture) distribu-tion with parameters generated from the MLP (MDN).For the LBN, since the network has stochastic units, wedirectly sample 20 times from the network. We align

Pred

icti

ons

Expe

rts

Original imageAvg. LAB error: 4.75

Tem

pera

ture

Tint

Expo

sure

Con

tras

tH

ighl

ight

sSh

adow

sW

hite

sBl

acks

Cla

rity

Vib

ranc

eSa

tura

tion

Brig

htne

ss

P

E

Tem

pera

ture

Tint

Expo

sure

Con

tras

tH

ighl

ight

sSh

adow

sW

hite

sBl

acks

Cla

rity

Vib

ranc

eSa

tura

tion

Brig

htne

ss

Tem

pera

ture

Tint

Expo

sure

Con

tras

tH

ighl

ight

sSh

adow

sW

hite

sBl

acks

Cla

rity

Vib

ranc

eSa

tura

tion

Brig

htne

ss�1.5�1.2�0.9�0.6�0.30.00.30.60.91.21.5

Pred

icti

ons

Expe

rts

Original imageAvg. LAB error: 3.22

Tem

pera

ture

Tint

Expo

sure

Con

tras

tH

ighl

ight

sSh

adow

sW

hite

sBl

acks

Cla

rity

Vib

ranc

eSa

tura

tion

Brig

htne

ss

P

E

Tem

pera

ture

Tint

Expo

sure

Con

tras

tH

ighl

ight

sSh

adow

sW

hite

sBl

acks

Cla

rity

Vib

ranc

eSa

tura

tion

Brig

htne

ss

Tem

pera

ture

Tint

Expo

sure

Con

tras

tH

ighl

ight

sSh

adow

sW

hite

sBl

acks

Cla

rity

Vib

ranc

eSa

tura

tion

Brig

htne

ss

�1.5�1.2�0.9�0.6�0.30.00.30.60.91.21.5

(a) (b)Figure 5: Multimodal photo edits: Sample slider predictions from the CGM-VAE model (denoted by P in the figure) comparedto the edits of 3 most active experts in the expert users dataset (denoted by E). The images are selected from the test subset of thedataset; the 3 samples are selected from a set of 10 proposals from the CGM-VAE model such that they align with the experts. Toshow the difference between the model and experts, we apply their sliders to the original image. For more examples, refer to thesupplementary material.

0 10 20 30 40 50 60 70 80 90 100

% of a user’s history

−4−2

02468

101214

(Nor

mal

ized

)Pr

edic

itiv

eLo

g-Li

kelih

ood Casual users

0 10 20 30 40 50 60 70 80 90 100

% of a user’s history

−20−15−10−5

05

101520 Frequent users

0 10 20 30 40 50 60 70 80 90 100

% of a user’s history

−10

0

10

20

30

40 Expert users

0 10 20 30 40 50 60 70 80 90 100

% of a user’s history

−5

0

5

10

15

20 All users

CasualFrequentExpert

Figure 6: Predictive log-likelihood for users in the test set of different datasets. For each user in the test set, we compute thepredictive log-likelihood of 20 images, given 0 to 30 images and their corresponding sliders from the same user. 30 sampletrajectories and the overall average± s.e. is shown for casual, frequent and expert users. The figure shows that knowing more aboutthe user (up to around 10 images) can increase the predictive log-likelihood. The log-likelihood is normalized by subtracting offthe predictive log-likelihood computed given zero images. Note the different y-axis in the plots. The rightmost plot is provided forcomparing the average predictive log-likelihood across datasets.

these samples to the experts’ edits and find the LAB er-ror between the expert-retouched image and the model-proposed image.

To report the results we average across the 3 experts andacross all the test images. The LAB error in Table 1 indi-cates that CGM-VAE model outperforms other baselinesin terms of predicting expert edits. Some sample editproposals and their corresponding LAB errors are pro-vided in Figure 5. This figure shows that the CGM-VAEmodel can propose a diverse set of edits that is reason-ably close to those of experts. For further examples seethe supplementary material.

4.2 Categorization and personalization

Next, we demonstrate how the CGM-SVAE model canleverage the knowledge from a user’s previous edits andpropose better future edits. For the users in the test sets ofall three datasets, we use between 0 and 30 image-sliderpairs to estimate the posterior of each user’s cluster mem-bership. We then evaluate the predictive log-likelihoodfor 20 other slider values conditioned on the images andthe inferred cluster memberships.

Figure 6 depicts how adding more image-slider com-binations can generally improve the predictive log-likelihood. The log-likelihood is normalized by subtract-

(a) (b)Figure 7: User categorization: Two examples of sample edits for three different user groups which the CGM-SVAE model hasidentified (in the experts dataset). (a) For similar flower photos, users in group I prefer to use low contrast and vibrance, whereasgroup II users tend to increase the exposure and vibrance from their default values. There is also group III users which do not showany specific preference for similar flower photos. (b) The same user groups for another set of similar photos with dominant bluecolors. For more examples, see the supplementary materials.

Dataset Casual Frequent ExpertEval. Metric LL JSD LL JSD LL JSD LABMLP −15.71± 0.21 0.26± 0.04 −2.72± 0.31 0.11± 0.02 −4.28± 0.12 0.22± 0.06 7.81± 0.26

LBN −7.12± 0.15 0.14± 0.02 −3.7± 0.43 0.13± 0.02 −4.89± 0.24 0.17± 0.04 7.44± 0.29

MDN −14.53± 0.25 0.31± 0.06 −1.67± 0.47 0.24± 0.08 −4.91± 0.07 0.28± 0.11 8.41± 0.27

CGM-VAE −6.39± 0.11 0.10± 0.02 −1.42± 0.18 0.08± 0.02 −2.6± 0.15 0.12± 0.05 6.72± 0.27

Table 1: Quantitative results: LL: Predictive log-likelihood for our model CGM-VAE and the three baselines. The predictive log-likelihood is calculated over the test sets from all three datasets. JSD: Jensen-Shannon divergence between normalized histogramsof the true sliders and our model predictions over the test sets (lower is better). See Figure 4 for an example of these histograms.LAB: LAB error between the images retouched by the experts and the images retouched by the model predictions. For each image wegenerate 3 proposals and compare that with the images generated by the top 3 active experts in the experts dataset.

ing off the predictive log-likelihood computed given zeroimages. The effect of adding more images is shown for30 different sampled users; the overall average for thetest dataset is also shown in the figure. To compare howvarious datasets benefit from this model, the average val-ues from the 3 datasets are overlaid. According to thisfigure, the frequent users benefit more than the casualusers and the expert users benefit the most. 2

To illustrate how the trained CGM-SVAE model pro-poses edits for different user groups, we use a set ofsimilar images in the experts dataset and show the pre-dicted slider values for those images. Figure 7 showshow the inferred user groups edit two groups of similarimages. This figure provides further evidence that the

2To apply the CGM-SVAE model to the experts dataset,we split the image-slider combinations from each of the 5 ex-perts into groups of 50 image-sliders and pretend that eachgroup belongs to a different user. This way we can have moreusers to train the CGM-SVAE model. However, this meansthe same expert may have some image-sliders in both train andtest datasets. The significant advantage gained in the expertsdataset might be due in part to this way of splitting the experts.Note that there are still no images shared across train and testsets.

model is able to propose a diverse set of edits across dif-ferent groups; moreover, it shows each user group mayhave a preference over which slider to use. For more ex-amples see the supplementary material.

5 Conclusion

We proposed a framework for multimodal predictionof photo edits and extend the model to make person-alized suggestions based on each user’s previous ed-its. Our framework outperforms several strong baselinesand demonstrates the benefit of having interpretable la-tent structure in VAEs. Although we only applied ourframework to the data from photo editing applications, itcan be applied to other domains where multimodal pre-diction, categorization and personalization are essential.Our proposed models could be extended further by as-suming more complicated graphical model structure suchas admixture models instead of the Gaussian mixturemodel that we used. Furthermore, the categories learnedby our model can be utilized to gain insights about thetypes of the users in the dataset.

ReferencesC. M. Bishop. Mixture density networks. 1994.

V. Bychkovsky, S. Paris, E. Chan, and F. Durand. Learningphotographic global tonal adjustment with a database of in-put/output image pairs. In Computer Vision and PatternRecognition (CVPR), 2011 IEEE Conference on, pages 97–104. IEEE, 2011.

Y. N. Dauphin and D. Grangier. Predicting distribu-tions with linearizing belief networks. arXiv preprintarXiv:1511.05622, 2015.

M. D. Hoffman, D. M. Blei, C. Wang, and J. Paisley. Stochas-tic variational inference. The Journal of Machine LearningResearch, 14(1):1303–1347, 2013.

R. Jaroensri, S. Paris, A. Hertzmann, V. Bychkovsky, andF. Durand. Predicting range of acceptable photographictonal adjustments. In Computational Photography (ICCP),2015 IEEE International Conference on, pages 1–9. IEEE,2015.

M. J. Johnson, D. Duvenaud, A. B. Wiltschko, S. R. Datta, andR. P. Adams. Composing graphical models with neural net-works for structured representations and fast inference. InNeural Information Processing Systems, 2016.

D. P. Kingma and M. Welling. Auto-encoding variationalbayes. arXiv preprint arXiv:1312.6114, 2013.

Y. Koyama, D. Sakamoto, and T. Igarashi. Selph: Progressivelearning and support of manual photo color enhancement. InProceedings of the 2016 CHI Conference on Human Factorsin Computing Systems, pages 2520–2532. ACM, 2016.

J.-Y. Lee, K. Sunkavalli, Z. Lin, X. Shen, and I. S. Kweon.Automatic content-aware color and tone stylization. arXivpreprint arXiv:1511.03748, 2015.

C. J. Maddison, A. Mnih, and Y. W. Teh. The concrete distribu-tion: A continuous relaxation of discrete random variables.arXiv preprint arXiv:1611.00712, 2016.

K. P. Murphy. Machine learning: a probabilistic perspective.MIT press, 2012.

R. M. Neal. Connectionist learning of belief networks. Artifi-cial intelligence, 56(1):71–113, 1992.

G. Sharma and R. Bala. Digital color imaging handbook. CRCpress, 2002.

K. Simonyan and A. Zisserman. Very deep convolu-tional networks for large-scale image recognition. CoRR,abs/1409.1556, 2014.

Y. Tang and R. R. Salakhutdinov. Learning stochastic feedfor-ward neural networks. In Advances in Neural InformationProcessing Systems, pages 530–538, 2013.

Z. Yan, H. Zhang, B. Wang, S. Paris, and Y. Yu. Automaticphoto adjustment using deep neural networks. ACM Trans-actions on Graphics (TOG), 35(2):11, 2016.