
 visual assessment of different skin pathologies is a result of ambient light that enters the skin and is scattered and diffused within it .
it is a combination of selective absorption and scattering of specific light wavelengths due to the physical properties of chromophores composing the skin .
well - trained dermatologists analyze the skin color and interpret the clinical pathologies based on their knowledge and experience .
dermatologists evaluate lesion conditions based on the distribution , size , shape , border , and symmetry but mostly on the color aspect .
diagnoses based on colour are subjective as colour perception depends on human visual response to light .
the human eye does not have the same sensitivity for all wavelengths   and between individuals .
colour is sensed by the human eye over the visible wave range and is subjectively interpreted as a unique sensation while it is a combination of wavelengths .
this lack of spectral discrimination means that the eye can be affected by metamerism which potentially affects the analysis .
imaging systems for skin analysis often try to mimic the eye analysis . nowadays , digital imaging systems are more and more available to clinicians , but imaging systems are mostly restricted to colour cameras .
such systems are limited in terms of spectral information as it is based on the trichromatic model .
however , it does not take advantage of the skin / light interaction which occurs over the whole spectrum range .
there is a growing awareness of skin disease condition worldwide . to improve the subjective assessment made on colour information
several studies [ 68 ] point out the possible differentiation of skin variegation above specific wavelength values as compared to healthy skin , meaning that spectral information is an important tool of assessment .
a reflectance spectrum provides precise objective physical information compared to subjective colour measurement . in vivo
optical spectroscopy is based on the study of light interaction of molecules with electromagnetic radiation .
spectroscopy measures the light intensity as a function of wavelength in form of a spectrum .
this type of measurement is linked to the optical property of the skin and is a result of the absorption , scattering , and emission properties of the skin .
currently , there is no system that can replace the diagnosis abilities of experienced clinicians .
however , the use of optical instruments increases the amount of complementary information to the dermatologist
. it can potentially provide information not detectable by the human eye and can lead to objective skin chromophore quantification .
this can be obtained by combining advantages of both spectrophotometer ( spectral resolution ) and digital camera ( spatial resolution ) .
multispectral imaging ( msi ) systems overcome their respective limitations ( lack of spatial variation and lack of spectral resolution ) .
the development of methods to assist diagnosis of skin pathology is based on objective assessment of skin characteristics .
the study of skin reflectance can be correlated with its biochemical and morphological composition to reveal information about its condition .
several researches base their skin parameters retrieval by multivariate methods such as partial least squares regressions , support vector machine ( svm ) , or blind source separation ( bss )   such as independent component analysis ( ica )   or principal component analysis ( pca )   to determine the concentration of skin chromophores .
these techniques assumed that skin reflectance is a linear combination of different source component spectra weighed by their mixing quantities .
generally , melanin and haemoglobin are assumed to be the two main components of the skin .
this category of analysis can be affected when the skin composition is different from the assumption .
these methods are based on strong hypothesis of the skin composition and the results might be influenced when one or more hypotheses are not satisfied .
another category refers to the analysis of reflectance spectra by means of physical models of light transportation that are based on optical skin properties .
different light propagation models have been developed such as model based on monte carlo simulations [ 15 , 16 ] , kubelka - munk [ 17 , 18 ] , and the modified beer - lambert law . the main motivation is to retrieve skin parameters by inversing the model to match a measured reflectance spectrum .
this category of analysis is based on a priori physical knowledge of skin absorption and scattering properties and thickness .
these methods tend to be more flexible regarding the skin composition and are less likely to be affected by unexpected data and to output incoherent results .
a research by shi and dimarzio   presents a hyperspectral imaging system applied to foot wound care and analyzed data with two methods ( beer - lambert law and two - layer optical model ) ; however , the system is limited to the visible range and only three parameters ( oxy / deoxyhaemoglobin concentration and epidermis thickness ) are retrieved .
have developed a multispectral imaging system covering visible and near - infrared light combine with an algorithm retrieving five parameters ( total melanin content of the epidermis and papillary dermis , collagen and haemoglobin content , and the presence of melanin in the papillary dermis ) but does not consider skin layer thickness . as a result ,
our development focuses on a technology combining a large amount of skin information ( both spectrally and spatially ) and analysis of the skin parameters quantification .
we address this problem by developing an exploratory imaging system retrieving skin parameter maps dedicated to dermatology / cosmetology . in this contribution
the msi reconstructs reflectance cubes which contain only the spectral reflectance of cutaneous data by means of an artificial neural - network algorithm .
hypercubes are analyzed by using an inverse model of light propagation which retrieves skin parameter maps based on the analysis of reflectance spectra by means of a physical model of light transportation in skin .
our motivation is to provide a system for prospective study of different skin lesions . the developed multispectral imaging ( msi ) system covers visible and near - infrared wave range and is suitable for dermatological and cosmetological constraints ( ergonomics and fast acquisition time ) .
the system presented is called asclepios , standing for analysis of skin characteristics by light emission and processing of images of spectrum .
an msi is generally composed of elements similar to colour acquisition system ; the main difference is the increased number of channels . for ergonomic purpose ,
the system is separated into two parts , an illumination system and a hand - held acquisition device ( see figure 2 ) .
this setting has the advantage of reducing the weight and size of the hand - held device and minimizes the number of calibration steps .
the illumination system is composed of a light source and a spectral selective device which is positioned in front of the illumination .
it is a 175-watts light bulb with spectral range from 380 to 1000  nm .
the spectral selective device is a filter wheel ( lambda 10 - 3 , model lb10-nw , sutter instrument , usa ) holding a set of ten interference filters ( 25  mm diameter ) which is placed in front of the light source in the light compartment .
unlike moncrieff et al .   who studied the optimal selection of spectral filters to recover human skin information in the visible range
, we consider equally divided filter in both the visible and near - infrared range .
the ten medium bandpass interference filters ( cvi melles griot , usa ) have been chosen from 420 to 960  nm .
the full width at half maximum of each filter is 80  nm ; there are overlays between the filters to compensate the central wavelength tolerance guarantees by the manufacturer ( 16  nm ) ( see figure 3 ) .
the choice of interference filters for the spectral decomposition is motivated by the large commercial offer available in terms of spectral band and transmittance .
light at specific waveband is transmitted from the illumination system to the back of the hand - held device by a liquid light guide .
the hand - held device is a light opaque device which protects the acquisition from external light ; thus , the skin is only illuminated with light at specific wavebands .
the device is composed of a wireless trigger , a camera , a lens , and liquid light guide .
the extremity of the hand - held device is composed of a nozzle which sets a constant focus distance of 10  cm between the skin area and the camera .
the selected camera consists of a 12  bits cmos digital sensor ( model mv1-d1312l-160-cl , photonfocus , switzerland ) with a resolution of 1312  1082 pixels and extended sensitivity in the visible and near - infrared ( 370 to 1100  nm ) .
a c - mount magnifying lens ( model mevis - c , linos , germany ) provides a useful area of 32  38  mm with a depth of field of 5  mm , yielding a spatial resolution of 33  pixelsmm . after positioning the hand - held device extremity on the skin area under study
the pressure triggers the acquisition which is performed in less than two seconds . if the patient is properly installed during the acquisition process , we consider the scene completely static
the system is operated using a personal computer running on windows xp operating system ( microsoft , seattle , wa , usa ) .
an in - house software , developed using visual c++ ( microsoft , seattle , washington ) , controls the filter wheel rotation , the image acquisitions , the reflectance cube reconstruction , and the database .
the interest of asclepios system is its capacity to reconstruct reflectance cubes from acquired multispectral images . to ensure reproducibility of the reconstruction
this specificity is linked to the different exposure time of each monoband image but also because dust or scratch might affect each filter .
a raw monoband image can be modelled as follows : 
 ( 1)[r]=[o]+[us ] , 
 where [ r ] is the raw image , [ o ] is the offset , [ u ] is the useful signal we want to extract , and [ s ] is the sensor response . using the cited modelled
, it is necessary to calibrate first the offset noise and then pixel gain . 

the offset corresponds to the intensity values that the camera acquires when observing a black target .
the offset frame is obtained by averaging multispectral images of a spectralon gray - scale standard target with 2% reflectance and is denoted [ o ] in ( 1 ) ( it is noted o which stands for offset frame ) .
the [ s ] is in most cases characterized by uneven intensities in an image of a uniform surface .
this effect is generally caused by nonuniform illumination , blemish on optical elements , difference of sensitivity within the sensor , and vignetting .
its goal is to compensate any variation within the detector for a given amount of light .
the acquisition of flat - field multispectral image is obtained by averaging multispectral image of a spectralon gray - scale standard target with 99% reflectance across the visible and is denoted f and stands for the flat - field image .
the gain correction image is composed of specific coefficient for each pixel of each monoband image .
the coefficients are obtained using the following formula : 
 ( 2)[s]=[f][o]dr , 
 where f is the flat - field monoband image and dr is the dynamic range of the camera . in the case of asclepios , dr is equal to 4.096 , corresponding to 12  bits images .
the final correction of the raw image is applied using ( 2 ) : 
 ( 3)[u]=[r][o][s ] . the cube reconstruction aims to extract only the skin reflectance spectra in each pixel of the data acquired . in order to reconstruct spectra linked to the physical properties of the skin element , the proposed method takes into account a model of light propagation in which the spectral response of all the elements involved in the acquisition process is employed ( figure 4 ) where i( ) is the spectral radiance of the illuminant , k( ) is the spectral transmittance of the kth filter , r( ) is the spectral reflectance of the surface , o( ) is the spectral transmittance of the optical system , ( ) is the spectral sensitivity of the camera and dk is the acquired image by the kth filter .
such model aims at separating each element of the acquisition chain to only keep the reflectance information r( ) .
there exist different methods to solve the inverse problem of spectral reflectance estimation of an object from the camera response .
there are pseudoinverse calculus methods and least - squares ones ; however , these methods are affected by noise amplification which is the major drawback .
another method is based on the a priori spectral reflectance information from the surface to be imaged .
it is this method that we select to perform the reconstruction of reflectance cubes because this technique is robust to noise and moreover it has generalization capabilities , meaning that it can reconstruct a wide range of reflectance spectrum , even ones that it has not learnt .
the reconstruction of reflectance cube is based on an artificial neural network- ( ann ) based algorithm proposed by mansouri et al . .
ann technique is a two - step process , a learning and a reconstruction part .
the neural network , employed in the reflectance cube reconstruction , uses heteroassociative memories due to its modularity with regard to the different size of the input and output vectors .
the ann generates a coefficient matrix which is obtained upon association of the multispectral image of the patches with their respective known spectral values . from figure 4 , the reflectance spectra r( ) of each pixel is reconstructed using the camera response d( ) and the result obtained from the ann learning step ( called coefficient matrix ) .
the reconstruction requires a multispectral image as input in order to output a reflectance cube of cutaneous data .
the reconstruction provides a 3-dimensional volume ( x , y ,  ) called a reflectance cube of cutaneous data where x and y represent the spatial dimensions and  the spectral one . the reconstruction of a reflectance cube is fast and simple as the operation is a product between the coefficient matrix and the camera response .
the use of heteroassociative memory allows reconstruction with different wave ranges within the capability ( 420 to 780  nm ) with a current step equal to ten nanometers .
the limitation of the 420780  nm window is a result of the learning protocol which is based on the gretag macbeth colorchart .
a study involving 150 healthy participants covering five of the six skin phototypes ( spt ) ( according to the fitzpatrick scale ) was performed using asclepios .
three acquisitions are performed at three different body locations , two facultative skin colour areas and a constitutive skin colour one .
constitutive skin colour is the natural , genetically determined colour of the skin whereas facultative skin colour is skin area affected by the environment ( sun , hormones , etc . )
the reconstruction validation is performed by comparing data acquired using a commercial spectrophotometer with reflectance cube reconstructed by asclepios .
the asclepios data used for comparison is of the same size and at similar location to that the data acquired by the spectrophotometer .
the results reveal that our system reconstructs reflectance cube with an average goodness of fit coefficient ( gfc ) greater than of 0.997 which , according to hernndez - andrs et al . , considers that the reconstruction is good if the gfc is above 0.99 .
the clinical data were acquired at the department of dermatology of hospital kuala lumpur under the supervision of dr .
r. baba ( head of dermatology department ) and hospital serdang ( malaysia ) with the collaboration of dr .
the data collection was performed with the collaboration of the department of electrical and electronic engineering of the university teknologi petronas ( malaysia ) .
the research study was registered to the national medical research register which supports the implementation of the national institute of health nih guideline on the conduct of research in the ministry of health malaysia ( moh ) .
we acquired a total of 22 melasma data from 10 patients and a total of 110 vitiligo data from 32 patients .
the acquisition process required the patient consent and all procedures were performed following the dermatology guideline .
the interest of the previously reconstructed reflectance cube is the retrieval of skin component parameters from each spatial area by mean of spectral analysis .
analysis of reflectance cubes can provide noninvasive evaluation of skin condition through the 2d mapping of relative meaningful skin chromophore concentrations .
as previously mentioned , we propose to extract skin parameter maps by using spectral analysis based on a light propagation model in skin .
the presentation of the method is divided into two sections , light propagation model description and genetic algorithm inversion method . the developed method for skin parameter maps by spectral analysis
the structure of human skin can be seen as a three - layer medium : epidermis , dermis , and subcutaneous fat .
several skin models have been developed with various number of skin parameters and different numbers of layers ( 2 , 3 , 5 , 7 , and even 22 layers ) .
however , regardless of the sophistication of the models , the numerous chromophores , and layers composing the skin , it is accepted that the human skin appearance , in terms of colour and reflectance spectrum in the visible domain , is mostly a result of melanin and haemoglobin concentration . following this assessment , we selected a skin model with two layers , epidermis and dermis .
the epidermis layer main component is the melanin   while the dermis main component is the haemoglobin .
we consider the epidermal thickness to be the effective melanin layer thickness and the dermal thickness as the effective haemoglobin layer thickness .
the physical and optical properties of the epidermis and dermis have been selected from the literature [ 1 , 30 , 3639 ] and are summarized in table 1 . the lower value of the epidermis thickness aims to take into account the epidermis thinning related to melasma disease . the light propagation model selected in this study is based on kubelka - munk ( k - m )   model .
it is originally based on a simple relationship between scattering and absorption coefficient of layers paint and its overall reflectance .
the k - m theory describes the radiation transfer in diffuse scattering media by applying energy transport equations .
k - m equations allow quantitative studies of absorption , scattering , and luminescence in diffuse scattering media .
k - m model has been extended to skin analysis by different groups [ 1 , 17 , 18 , 41 ] .
k - m model has an analytical form and it offers rapid skin optical parameters determination using inversion procedure .
we detail the principle of the kubelka - munk model in terms of reflectance and transmittance for a single layer .
it is based on the thickness of the layer dlayer , the layer absorption coefficient a.layer , and the layer scattering coefficient s.layer .
the equations related the melanin and haemoglobin absorption with the layer absorption coefficient a.layer and the layer scattering coefficient s.layer are detailed , respectively , for the epidermis and the dermis .
the epidermis absorption property is mostly due to the melanin and lightly by the baseline absorption coefficient .
the melanin spectral absorption coefficient   is approximated by 
 ( 4)a.melanin()=6.610113.33  [ cm1 ] .
the combined absorption of the different negligible skin components ( carotene , keratin , and collagen ) is taken into account by jacques which defines an equation for the baseline absorption coefficient free of the major chromophores which is similar for these two layers : 
 ( 5)a.baseline()=0.244 + 85.3exp(164)66.2[cm1 ] . for simplification purpose ,
epidermis and stratum corneum are regarded as a single layer because the stratum corneum light absorption is low and transmits light uniformly in the visible wave range .
the optical absorption coefficient of the epidermis a.epidermis is expressed as a function of the wavelength and depends mostly on the volume fraction of melanosome and lightly on the baseline skin absorption coefficient : 
 ( 6)a.epidermis()=fmela.melanin()+(1fmel)a.baseline()[cm1 ] , 
 where  is the wavelength in nanometres and a.melanin is in cm .
the fmel refers to the concentration of melanin in % and a.baseline is the baseline absorption in cm . in the dermis
, the absorption is performed by the main chromophore , the blood . within the red blood cell ,
the haemoglobin is a major absorber and it is decomposed into oxyhaemoglobin ( hbo2 ) and deoxyhaemoglobin ( hb ) components .
the values of absorption coefficient for the deoxy- and oxyhaemoglobin were obtained from the oregon medical laser centre website .
oxyhaemoglobin absorption spectra have two absorption peaks at around 542 and 578  nm , revealing the characteristic  w  shape . for this study ,
the oxygenated haemoglobin can be estimated using the concentration of oxy- and deoxyhaemoglobin by the following equation:(8)rso2=chbo2chbo2+chb[% ] .
the dermal absorption coefficient a.dermis is expressed by 
 ( 9)a.dermis()=fblood(chbo2a.oxy())+fblood(1chbo2)a.deoxy()+(1fblood)a.baseline( )  [ cm1 ] , 
 where  is the wavelength in nanometres .
chbo2 is the concentration of oxyhaemoglobin in blood . a.oxy and a.deoxy refer , respectively , to the absorption coefficient of the oxyhaemoglobin and deoxyhaemoglobin in cm .
subcutaneous layer is ignored because its main function is light absorption and mostly contains fat   and limited visible light reaches this layer . 

the equations to calculate the reflectance r and transmittance t for one layer defined by kubelka - munk are expressed by 
 ( 10)rlayer()=(1()2)(exp(k()dlayer )  exp(k()dlayer))((1+())2exp(k()dlayer )  (1())2exp(k()dlayer))1,tlayer()=4((1+)2exp(kdlayer )  (1)2exp(kdlayer))1 , 
 where k is the backward flux variable of one layer expressed by 
 ( 11)klayer = klayer(klayer+2slayer ) , 
 and  is the forward flux variable of one layer expressed by 
 ( 12)layer = klayerklayer+2slayer , 
 where k and s are expressed by 
 ( 13)klayer=2alayer , slayer=2slayer .
the total reflection of a two - layer medium is defined by the following equation : 
 ( 14)rtotal = rl1l2=r1+tl12rl21rl1rl2,ttotal = tl1l2=tl1tl21rl1rl2 , 
 where l1 and l2 , respectively , refer to the epidermis layer and the dermis layer .
the model allows fast computation of the total reflectance and transmittance of a medium based on the absorption and scattering coefficient and thickness of each layer .
the five biological parameters used to simulate a reflectance spectra can be represented by an input function denoted by 
 ( 15)p=(fmel , depi,fblood , coxy , ddermis ) .
we presented a model of light propagation based on a two - layer skin medium using the kubelka - munk theory .
the model generates reflectance spectra using five physiological parameters that can vary . from our research
this inverse problem is nonlinear due to the complex structure , in terms of scattering and absorption properties , of the two - layer skin model .
the inversion means to reverse the model of light propagation in which parameters are input to generate a reflectance spectrum .
the inversion of our model retrieves the five parameters by matching simulated spectra generated by the model with measured one . among the different existing techniques of search and optimization , our interest focused on genetic algorithm ( ga ) .
ga is a metaheuristic method that optimizes a problem through iterative improvement of a candidate solution with a quality measurement .
the search of a global solution using ga is done on a population of candidates rather than a single candidate .
the strengths of ga are its capacity to explore large search space ( through mimic of natural evolution ) , and most importantly , its strength against becoming trapped into a local minimum ( by introduction of random search ) .
ga requires a representation of the solution space and a fitness function to evaluate the population of individuals .
one individual is defined in real value of the five skin parameters of the k - m model and is bounded by its physical limit defined in table 1 from the literature .
attempt using ga to retrieve skin parameters of a monte carlo for multilayered ( mcml ) media has been tested by zhang et al .   and a hybrid version of it was developed by choi .
both prove the usefulness of the ga , but these techniques also retrieve unrealistic values mainly due to the nonconstraint of the search space . in our version ,
the boundary condition of each element ( minimum and maximum parameter values ) is restricted to their physical limit .
this nonrestriction aims to be applied on any data taken wherever on the body , however ; this might potentially influence the results .
the evolution of the population is roughly based on an iterative three - step process : selection , crossover , and mutation ( see figure 5 ) .
it is probably the most important aspect of the genetic algorithm as it requires classifying the best individual which is the most similar to the measured spectrum . 

the genes of a chromosome are input into the forward km model to generate a simulated spectrum .
the fitness function calculates the similarity between a measured spectrum and a simulated one ( using the forward km model ) . the calculated fitness value is used to classify the population from the best to worst one for the selection process . in the search of an optimal fitness function for our ga application ,
the following five different metric scales have been tested : the root mean squared error ( rmse ) , the goodness of fit coefficient ( gfc ) , the reconstruction percentage ( recp ) , the modified spectral angle similarity ( msas ) , and the spectral similarity value ( ssv ) .
the first stage of the evolution process is to apply genetic operation to the population t which will generate an intermediate population t ( called mating pool ) . for the evolution ,
then , most of the individuals from the intermediate population will be kept as part of the new population for the next generation . only few individuals from the intermediate population mutate before being added to the new population . during the next step of the evolution , the random operator selects individuals from the population and places them in the t population .
crossover operation consists of swapping one randomly selected parameter between two randomly chosen parents to generate two offsprings .
the mutation alters randomly one parameter of the child chromosome by generating a new random value for the parameter .
there are different criteria to stop the evolution of the algorithm . due to the difficulty to reach a specific fitness function value for different spectra
this technique has the advantage of setting a finite computation time for the algorithm and , if properly selected , one can reach optimum result .
however , at the final iteration , the algorithm may not have yet converged . as there is no golden rule to define the value of the different parameters of the ga ,
rmse was selected because it performed better than four other metrics ( reconstruction percentage , goodness of fit coefficient , modified spectral angle similarity , and spectral similarity value ) tested in the optimization process of the ga .
the proposed method of combining kubelka - munk model ( with a layered skin structure ) with ga optimization can retrieve accurately simulated skin parameters as suggested by the results .
the accuracy is defined by the retrieval error of two characteristic sets of parameters ( one lightly and one darkly pigmented skin ) with a white gaussian noise with amplitude of 0.1 added to the simulated spectra .
the average root mean squared error of the fitness function is 0.25  10 with an average error of less than 1.5 percent for each parameter .
this algorithm outputs the concentration of melanin , the concentration of haemoglobin , the thickness of both epidermis ( considered as the melanin layer ) and dermis ( considered as the haemoglobin layer ) , and the oxygenated haemoglobin .
the genetic algorithm routine and the kubelka - munk model have been developed on matlab ( the mathworks , inc .
, ma , usa ) running on a personal computer running on windows xp operating system ( microsoft , seattle , wa , usa ) .
quantitative and objective assessment of skin lesions is critical in the detection of variety of skin conditions . in order to use our method for clinically relevant skin analysis
, we first validated our technique using a population of healthy participants and then applied our method to the measurement of parameters from patients with vitiligo or melasma skin lesions .
the choice to study these diseases is based on the generic target of the system .
these diseases are easily visually assessed by dermatologist with known parameter variation ( melanin ) .
another experiment to validate the oxygenated haemoglobin percentage is performed on healthy participants using images of normal tissue and tissue with temporary induced ischemia .
we present the results of our method based on participants from a healthy population and compare the results obtained with data from the literature .
the volume fraction of melanosome increases with the spt type which is in accordance with the literature [ 33 , 45 , 49 ] .
currently , the only method to classify and compare our system is the relation between spt and melanin concentration in skin and relation between healthy skin and skin with temporary induced ischemia .
the data for our method of skin parameters retrieval are the same as the data used for the validation of the reflectance cube reconstruction . for each reflectance cube ,
the distribution of the skin is fairly uniform as the data were taken from healthy patient and care was taken to acquire with no visible lesions or nevus .
the average value of each parameter was calculated for each entire retrieved map because of low disparity within values of each map .
the aim of this test is to corroborate the skin parameter values / concentration from the literature with the parameters retrieved using our method on skin sample from various spt .
our results , presented by body location , show the different retrieved parameters for different spt in tables 3 , 4 , and 5 .
the results are fairly consistent with the finding of robertson and rees   which reveal a thicker epidermis layer for the back of the hand which agrees with our retrieved parameters .
the decrease of thickness measured by reflectance confocal microscopy ( rcm ) again matches the decrease of thickness calculated by our method .
the major difference between the body locations is found for the melanin content of the back of the hand which is twice the one of face and lower back for spt ii , iii , but the trends , while slightly lower , are similar for the remaining spt .
the oxygenated haemoglobin variation between skin types ii and iv can be mostly noticed for the face and back location . this effect may be linked to the increase of melanin concentration . as the melanin concentration increases , the characteristic  w  shape of the oxygenated haemoglobin is more difficult to notice on the skin reflectance spectrum . the oxygenated haemoglobin algorithm estimation might be affected when there is a high concentration of melanin , hence the difference in the results .
also the overall values of oxygenated haemoglobin are on average 20% lower than the one reported by zonios et al .   but similar to the one from tseng et al . .
the difference may be a result of the interrogation depth limit of our system and requires further investigations .
the volume fraction of haemoglobin is very low for spt vi compared to other skin types ; this might be a limitation of the system .
it seems to be affected by the low variation of the reflectance spectrum for spt vi .
the accuracy of the method may need to be improved for this specific phototype . a possible solution to overcome this limitation is to select specific spectrum range ( possibly around the characteristic  w  shape ) to increase its accuracy .
further investigation is required to analyze the variation of the volume fraction of haemoglobin for the back site .
the model seems to be affected by the relative flat spectrum of spt vi , especially with regard to the dermis thickness and volume fraction of haemoglobin .
first , using ultrasound technology , the determination of the accuracy of the layer thicknesses may be verified .
we applied our method to quantify and compare data acquired from two populations with skin diseases that have characteristic effect on the melanin composition .
melasma ( also called chloasma ) is a hyperpigmentation skin disease that is characterized by an increase level of melanin concentration released by melanocytes .
the symptoms are characterized by dark , irregular , and well - demarcated skin lesions .
vitiligo ( also called leukoderma ) is a common genetic autoimmune skin disease caused by the disappearance of melanocytes in the epidermis resulting in hypopigmentation .
the average parameter values are obtained from the two different areas ( healthy and diseased ) which are manually selected on each map .
the relative difference between healthy and hypo / hyperpigmented skin area is calculated by the following formula : 
 ( 16)rdiff = valrefvalmeasured , 
 where valref is the retrieved healthy skin parameter value and valmeasured is the retrieved hypo / hyperpigmented skin parameter value .
tables 6 and 7 summarize the difference of parameters between healthy and disease affected reflectance spectra .
the variation of the other parameters ( epidermis and dermis thickness , haemoglobin concentration , and oxygenated haemoglobin ) is not significant .
the standard deviations of the average melanin difference for both melasma and vitiligo are high .
this is a result of the inclusion of all types of vitiligo ( mild to severe ) and all spt types . for spt ii and iii
, the difference between melasma and healthy skin is very high , sometimes with a relative difference of 4 , whereas for spt v and vi , the relative difference is around 1.3 to 1.5 as the original melanin concentration is already high . 

figure 7 presents the results of the five parameters retrieved using our method on a vitiligo lesion .
the melanin map is clearly affected in the area of vitiligo lesions ( fair patch on the subfigure ( a ) ) . in order to highlight the interest of the result obtained by the proposed method
it is to be noted that the code has not been optimized as the aim is to show the relevance of the method for skin analysis .
the variation of epidermis / dermis thickness , haemoglobin concentration , and oxygenated haemoglobin parameters is only of around 10% which we do not consider to be characteristic for melasma and vitiligo diseases .
this finding is consistent with the medical expectation which states that vitiligo only modifies the melanin content of the skin . as previously noted for healthy data , a difference of dermis thickness
is observed between the different spts which does not correspond to a histological characteristic as spt difference is due to a change of melanin concentration .
this is particularly noticeable for a specific lesion ( figure 7 ) where the cross - correlation between the melanin concentration map and the dermis thickness map is 0.41 .
the algorithm might be affected by the limited information in the nir to accurately retrieve dermis thickness in some case meaning a possible crosstalk in the model .
this will be investigated further by adding layer to the model to try to decorrelate these parameters .
the experiment aims to induce temporary ischemia to the lower arms part of healthy volunteers .
ischemia is defined as a restriction in blood supply to tissues , leading to a decrease of oxygen ( and glucose ) .
the first acquisition is a baseline data acquired with the subject at rest for 5 minutes .
temporary induced ischemia is simulated by inflating a pressure cuff to the participant left upper arm .
the pressure cuff blocks the blood flow and is maintained for two minutes before performing the second acquisition .
finally , the third acquisition is taken one minute after reperfusion ( release of the cuff pressure ) .
the three multispectral images of the lower inner arm , acquired during the experiment , are reconstructed into reflectance cube and analyzed by the previously defined algorithm .
the data are processed by the algorithm and the average parameter values are obtained from identical areas manually selected on each map . 

the values of melanin concentration , epidermis thickness , and dermis thickness do not show significant variation between the three acquisitions .
the occlusion does not have effect on these parameters , which is the result expected .
the percentage of difference between the baseline and acquisition 2 and acquisition 3 is calculated by the following formula : 
 ( 17)pdiff=(valacqvalbaseline)valbaseline , 
 where valacq is the retrieved skin parameter value for acquisition ( either number 2 or 3 ) and valbaseline is the retrieved baseline skin parameter value .
the occlusion effect ( see table 9 ) shows an increase in volume fraction of haemoglobin and a decrease in oxygenated haemoglobin concentration compared to baseline .
it is consistent with finding by vogel et al .   and the decreased value of oxygenated haemoglobin of 22.7% is similar to the decrease reported by zuzak et al .  
the retrieve values of both oxygenated blood concentration and volume of blood fraction for the third acquisition compared to baseline reveal only a slight increase ( + 6% ) .
this potentially means that one minute after the release of the pressure cuff , the hyperemia ( increase of blood flow and oxygenated haemoglobin ) is nearly over and that the tissues are returning to homeostasis .
the temporary induced ischemia experiment demonstrates the potential of our system to record oxygenated haemoglobin and volume fraction of haemoglobin .
further studies might be required to fully assess the applicability and accuracy of our system for both dermatological and cosmetological applications .
we presented a multispectral camera which has the capacity to reconstruct reflectance cubes only linked to the skin characteristics .
the reflectance cube allows quantitative measure of cutaneous tissue based on a kubelka - munk model combined with evolutionary algorithm . using this approach , quantification maps of five skin parameters ( melanin concentration , epidermis / dermis thickness , haemoglobin concentration , and the oxygenated haemoglobin )
the developed algorithm was tested on a set of healthy skin data acquired using our system .
the results retrieved by the algorithm are in good agreement with the data from the literature .
finally , the usefulness of the developed technique was tested during a clinical study based on healthy skin , vitiligo , and melasma skin lesions .
the results show that our research method can retrieve skin parameters in accordance with the expected skin composition for each lesion .
for example , with regard to the vitiligo , where the main characteristic is a lack of melanin , the method clearly shows a decrease in melanin concentration . in the case of the melasma
the first step consisting in developing a generic system for the multispectral analysis of skin data has been achieved .
because the first tests give promising results , our in vivo imaging system may be useful for prospective study for a use in dermatological and cosmetological studies .
specific pathologies with strong socioeconomic impacts could be chosen to reveal possible lesion spectral signature .
the work perspectives involve modifying the system to be disease - specific ( in terms of wavelength selection , field of view , contact / noncontact  ) and to further validate the system for these specific purposes . to that regard ,
the perspective of validation for all parameters includes the use of histology ( although it is difficult to set up ) or measurements of skin layer depth by ultrasound imaging system .
 in vivo quantitative assessment of skin lesions is an important step in the 
 evaluation of skin condition . 
 an objective measurement device can help as 
 a valuable tool for skin analysis . 
 we propose an explorative new multispectral camera specifically developed for dermatology / cosmetology applications . 
 the multispectral imaging system provides images of skin reflectance at different wavebands covering visible and near - infrared domain . 
 it is coupled 
 with a neural network - based algorithm for the reconstruction of reflectance 
 cube of cutaneous data . 
 this cube contains only skin optical reflectance 
 spectrum in each pixel of the bidimensional spatial information . 
 the reflectance cube is analyzed by an algorithm based on a kubelka - munk model 
 combined with evolutionary algorithm . 
 the technique allows quantitative 
 measure of cutaneous tissue and retrieves five skin parameter maps : melanin 
 concentration , epidermis / dermis thickness , haemoglobin concentration , and 
 the oxygenated hemoglobin . 
 the results retrieved on healthy participants 
 by the algorithm are in good accordance with the data from the literature . 
 
 the usefulness of the developed technique was proved during two experiments : a clinical study based on vitiligo and melasma skin lesions and a skin oxygenation experiment ( induced ischemia ) with healthy participant where 
 normal tissues are recorded at normal state and when temporary ischemia is 
 induced . 
1. Introduction
2. Materials and Methods
3. Results and Discussion
4. Conclusion