A 1064 nm Dispersive Raman Spectral Imaging System for Food Safety and Quality Evaluation

: Raman spectral imaging is an effective method to analyze and evaluate the chemical composition and structure of a sample, and has many applications for food safety and quality research. This study developed a 1064 nm dispersive Raman spectral imaging system for surface and subsurface analysis of food samples. A 1064 nm laser module is used for sample excitation. A bifurcated optical ﬁber coupled with Raman probe is used to focus excitation laser on the sample and carry scattering signal to the spectrograph. A high throughput volume phase grating disperses the incoming Raman signal. A 512 pixels Indium-Gallium-Arsenide (InGaAs) detector receives the dispersed light signal. A motorized positioning table moves the sample in two-axis directions, accumulating hyperspectral image of the sample by the point-scan method. An interface software was developed in-house for parameterization, data acquisition, and data transfer. The system was spectrally calibrated using naphthalene and polystyrene. It has the Raman shift range of 142 to 1820 cm − 1 , the spectral resolution of 12 cm − 1 at full width half maximum (FWHM). The spatial resolution of the system was evaluated using a standard resolution glass test chart. It has the spatial resolution of 0.1 mm. The application of the system was demonstrated by surface and subsurface detection of metanil yellow contamination in turmeric powder. Results indicate that the 1064 nm dispersive Raman spectral imaging system is a useful tool for food safety and quality evaluation. Increasing the offset distance enhanced the 1631 cm the layer The 1631 cm − 1 peak can be used to identify the turmeric packed inside gelatin capsule. The results demonstrate this SORS measurement technique has a potential for non-destructive subsurface detection and the identification of gelatin-encapsulated pure and mixed samples of turmeric and metanil yellow powder. Increasing the offset distance enhanced the 1631 cm − 1 peak in the multiple layer capsule samples. The 1631 cm − 1 peak can be used to identify the turmeric packed inside gelatin capsule. The results demonstrate this SORS measurement technique has a potential for non-destructive subsurface detection and the identiﬁcation of gelatin-encapsulated pure and mixed samples of turmeric and metanil yellow powder.


Introduction
Incidents of foodborne illness outbreaks have necessitated the development of a reliable technique for food safety and quality evaluation. Advanced sensing techniques, such as optical methods, hold promise for food authentication [1][2][3]. Among the optical techniques, Raman spectroscopy has an advantage due to its higher specificity and sensitivity. Simple spectral measurement method and insensitivity to the presence of water in the sample (eliminating the need for the removal of water from the sample) are additional advantages of Raman spectroscopic technique. Unique fingerprints for chemicals in Raman spectra allow analyses of a wide range of chemicals. It has been widely used to evaluate the chemical properties of food samples [4][5][6].
Raman spectra can be used to identify and evaluate one or more chemicals in a mixture. Raman spectroscopy has a broad application in food processing [7][8][9][10]. It has also been applied to food quality and safety evaluation [11][12][13][14][15][16][17][18]. Commercially available backscattering Raman spectroscopy, transmission Raman spectroscopy, and Fourier transform Raman (FT-Raman) spectroscopy have been used to evaluate food samples. Backscattering Raman spectroscopy is suitable for surface analysis, 1.
develop a 1064 nm dispersive Raman spectral imaging system for food safety and quality evaluation; 2.
present the detailed description of the system, hardware components and interface software; 3. demonstrate application example of the 1064 nm dispersive Raman spectral imaging system for detection of metanil yellow contamination in turmeric powder at different concentrations; and, 4. demonstrate application for non-destructive subsurface detection and identification of gelatin-encapsulated pure and mixed samples of metanil yellow and turmeric powder by SORS method. Figure 1 shows a schematic of the 1064 nm dispersive Raman spectral imaging system. The Raman spectrograph uses a high throughput volume phase grating (VPG) (BaySpec, Inc., San Jose, CA, USA) optimized for 1064 nm laser excitation. A concave lens (Lens 1) in the spectrograph guides the incoming light to the VPG, where the light is diffracted into different angular output paths. The concave lens (Lens 2) reflects the dispersed light to the 512 pixels Indium-Gallium-Arsenide (InGaAs) detector (Nunavut, BaySpec, Inc., San Jose, CA, USA). The detector is thermoelectrically deep cooled to −55 • C during spectral acquisition to minimize the dark current. The detector is connected to the computer using a USB cable for detector control and data transfer.

System Design
A 1064 nm laser module (MiniLite, BaySpec, Inc., San Jose, CA, USA) with a power of 600 mW is used for sample excitation. A fiber-optic Raman probe (RPB, InPhotonics, Inc., Norwood, MA, USA) is used to focus the laser on the sample surface and acquire the scattering signal. A long-pass filter assembled in the probe eliminates the light at and below 1064 nm. The probe is coupled to the Raman spectrograph via a bifurcated optical fiber that consists of excitation fiber and collection fiber. The excitation fiber (diameter 105 µm) carries the laser from the source to the Raman probe, is connected to the output of laser source by angled physical contact (APC) connector. The laser is focused on the sample surface by an extension tip of the probe. The end face polished at an 8-degree angle in the APC connector allows any light that is reflected back to the source to reflect out into the fiber cladding, achieving a better return loss. The collection fiber (diameter 200 µm) transfers the acquired scattering Raman signal to the spectrometer, is connected to the spectroscopy by fiber-optic connector (FC). The Raman probe uses bandpass, diachronic and edge filters for separating the excitation light and the scattered signal. The probe has a focal length of 7.5 mm and spot size of 158 µm. A two-axis motorized positioning table (MAXY4009W1-S4, Velmex, Bloomfield, NY, USA) is used to hold and move the sample in two perpendicular directions. The positioning table can move through an area of 127 mm × 127 mm with a displacement resolution of 6.35 µm. Stepper motor controller (Velmex, Bloomfield, NY, USA) is programmed to control its movement. Sample movement at an adjustable step size below the fixed position Raman probe allows accumulation of Raman spectral image of the entire surface area of the sample, which can be analyzed both spectrally and spatially. 2. present the detailed description of the system, hardware components and interface software; 3. demonstrate application example of the 1064 nm dispersive Raman spectral imaging system for detection of metanil yellow contamination in turmeric powder at different concentrations; and, 4. demonstrate application for non-destructive subsurface detection and identification of gelatin-encapsulated pure and mixed samples of metanil yellow and turmeric powder by SORS method. Figure 1 shows a schematic of the 1064 nm dispersive Raman spectral imaging system. The Raman spectrograph uses a high throughput volume phase grating (VPG) (BaySpec, Inc., San Jose, CA, USA) optimized for 1064 nm laser excitation. A concave lens (Lens 1) in the spectrograph guides the incoming light to the VPG, where the light is diffracted into different angular output paths. The concave lens (Lens 2) reflects the dispersed light to the 512 pixels Indium-Gallium-Arsenide (InGaAs) detector (Nunavut, BaySpec, Inc., San Jose, CA, USA). The detector is thermoelectrically deep cooled to −55 °C during spectral acquisition to minimize the dark current. The detector is connected to the computer using a USB cable for detector control and data transfer.

System Design
A 1064 nm laser module (MiniLite, BaySpec, Inc., San Jose, CA, USA) with a power of 600 mW is used for sample excitation. A fiber-optic Raman probe (RPB, InPhotonics, Inc., Norwood, MA, USA) is used to focus the laser on the sample surface and acquire the scattering signal. A long-pass filter assembled in the probe eliminates the light at and below 1064 nm. The probe is coupled to the Raman spectrograph via a bifurcated optical fiber that consists of excitation fiber and collection fiber. The excitation fiber (diameter 105 µm) carries the laser from the source to the Raman probe, is connected to the output of laser source by angled physical contact (APC) connector. The laser is focused on the sample surface by an extension tip of the probe. The end face polished at an 8-degree angle in the APC connector allows any light that is reflected back to the source to reflect out into the fiber cladding, achieving a better return loss. The collection fiber (diameter 200 µm) transfers the acquired scattering Raman signal to the spectrometer, is connected to the spectroscopy by fiber-optic connector (FC). The Raman probe uses bandpass, diachronic and edge filters for separating the excitation light and the scattered signal. The probe has a focal length of 7.5 mm and spot size of 158 µm. A two-axis motorized positioning table (MAXY4009W1-S4, Velmex, Bloomfield, NY, USA) is used to hold and move the sample in two perpendicular directions. The positioning table can move through an area of 127 mm × 127 mm with a displacement resolution of 6.35 µm. Stepper motor controller (Velmex, Bloomfield, NY, USA) is programmed to control its movement. Sample movement at an adjustable step size below the fixed position Raman probe allows accumulation of Raman spectral image of the entire surface area of the sample, which can be analyzed both spectrally and spatially.

Software
Interface software of the Raman system was developed on a LabView platform (National Instruments, Austin, TX, USA) for parameter setup and data transfer. The interface software controls system operational parameters, such as initialization, exposure time, spectral acquisition and display, sample movement, and data transfer. The software development kits (SDKs) that were provided by manufacturers of InGaAs detector and positioning table were used to develop the interface software. The acquired hyperspectral data are stored in the band interleaved by pixel (BIP) format, which can be commercially analyzed by available software, such as ENVI (4.5, ITT Visual Information Solutions, Boulder, CO, USA) and Matlab (R2013a, MathWorks, Natick, MA, USA).

Spectral Calibration and Resolution
The X-axis in the Raman spectral measurement, termed as Raman shift (wavenumber, cm −1 ), is the change in the energy state from the excitation source. Typically, chemicals with known wavenumber shifts are used to assign a wavenumber to the corresponding pixels along the spectral dimension. The process is called spectral calibration. In this study, spectral calibration was done using polystyrene and naphthalene. American Society for Testing and Materials (ASTM) International enlists six chemicals (polystyrene, naphthalene, 1,4-bis(2-methylstyryl) benzene, sulfur, toluene-acetonitrile (concentration 50/50 by volume), 4-acetamidophenol, benzonitrile, and cyclohexane) as Raman shift standards for spectrometer calibration [46]. Polystyrene and naphthalene were packaged separately in two Petri-dishes and scanned covering the surface area of 10 mm × 10 mm with a step size of 1 mm for X-and Y-axes, collecting a total of 100 spectra from each Petri-dish. Figure 2a shows an average Raman spectrum of naphthalene (polystyrene spectrum not shown). Seven spectral peaks were identified for naphthalene between 513 to 1576 cm −1 . Similarly, five peaks were identified for polystyrene in the wavenumber range of 620 to 1602 cm −1 . Spectral calibration was performed using the twelve peaks covering the wavenumber range of 513 to 1602 cm −1 . Quadratic fitting was used to develop spectral calibration model. The quadratic model was used to determine the wavenumbers of corresponding pixels. The Raman system covered the wavenumber range of 142 to 1820 cm −1 (Figure 2b) after spectral calibration. The spectral resolution of the Raman spectrometer is 12 cm −1 at full width half maximum (FWHM).

Software
Interface software of the Raman system was developed on a LabView platform (National Instruments, Austin, TX, USA) for parameter setup and data transfer. The interface software controls system operational parameters, such as initialization, exposure time, spectral acquisition and display, sample movement, and data transfer. The software development kits (SDKs) that were provided by manufacturers of InGaAs detector and positioning table were used to develop the interface software. The acquired hyperspectral data are stored in the band interleaved by pixel (BIP) format, which can be commercially analyzed by available software, such as ENVI (4.5, ITT Visual Information Solutions, Boulder, CO, USA) and Matlab (R2013a, MathWorks, Natick, MA, USA).

Spectral Calibration and Resolution
The X-axis in the Raman spectral measurement, termed as Raman shift (wavenumber, cm −1 ), is the change in the energy state from the excitation source. Typically, chemicals with known wavenumber shifts are used to assign a wavenumber to the corresponding pixels along the spectral dimension. The process is called spectral calibration. In this study, spectral calibration was done using polystyrene and naphthalene. American Society for Testing and Materials (ASTM) International enlists six chemicals (polystyrene, naphthalene, 1,4-bis(2-methylstyryl) benzene, sulfur, toluene-acetonitrile (concentration 50/50 by volume), 4-acetamidophenol, benzonitrile, and cyclohexane) as Raman shift standards for spectrometer calibration [46]. Polystyrene and naphthalene were packaged separately in two Petri-dishes and scanned covering the surface area of 10 mm × 10 mm with a step size of 1 mm for X-and Y-axes, collecting a total of 100 spectra from each Petri-dish. Figure 2a shows an average Raman spectrum of naphthalene (polystyrene spectrum not shown). Seven spectral peaks were identified for naphthalene between 513 to 1576 cm −1 . Similarly, five peaks were identified for polystyrene in the wavenumber range of 620 to 1602 cm −1 . Spectral calibration was performed using the twelve peaks covering the wavenumber range of 513 to 1602 cm −1 . Quadratic fitting was used to develop spectral calibration model. The quadratic model was used to determine the wavenumbers of corresponding pixels. The Raman system covered the wavenumber range of 142 to 1820 cm −1 (Figure 2b) after spectral calibration. The spectral resolution of the Raman spectrometer is 12 cm −1 at full width half maximum (FWHM).

Spatial Resolution
The spatial resolution of the 1064 nm Raman spectral imaging system was evaluated using a standard resolution glass test chart (Edmund Optics Inc., Barrington, NJ, USA). The test chart was placed on the top of a Petri-dish containing naphthalene. The 1064 nm laser source illuminated the test chart. A step size of 0.1 mm was used to scan X and Y directions. Figure 3 shows the single band image of the resolution test chart. The laser spot on the white space penetrated through the glass test chart to the naphthalene, while the illumination on the black dots suppressed the Raman spectra of naphthalene, creating a contrast between white space and the black dot. In the Figure 3, the small inner dots have the diameter of 0.25 mm with 0.5 mm distance between each dot. The outer larger dots have the diameter of 0.5 mm with 1 mm dot spacing. Both the dots and their dot spacing can be distinguished in the image due to the small step size used to scan the chart. The step size used to scan the X and Y direction is the spatial resolution of the Raman image. Figure 3 shows that the system can achieve the spatial resolution of 0.1 mm.

Spatial Resolution
The spatial resolution of the 1064 nm Raman spectral imaging system was evaluated using a standard resolution glass test chart (Edmund Optics Inc., Barrington, NJ, USA). The test chart was placed on the top of a Petri-dish containing naphthalene. The 1064 nm laser source illuminated the test chart. A step size of 0.1 mm was used to scan X and Y directions. Figure 3 shows the single band image of the resolution test chart. The laser spot on the white space penetrated through the glass test chart to the naphthalene, while the illumination on the black dots suppressed the Raman spectra of naphthalene, creating a contrast between white space and the black dot. In the Figure 3, the small inner dots have the diameter of 0.25 mm with 0.5 mm distance between each dot. The outer larger dots have the diameter of 0.5 mm with 1 mm dot spacing. Both the dots and their dot spacing can be distinguished in the image due to the small step size used to scan the chart. The step size used to scan the X and Y direction is the spatial resolution of the Raman image. Figure 3 shows that the system can achieve the spatial resolution of 0.1 mm.

Application to Authenticate Turmeric Powder
Turmeric (Curcuma long L.), which is a yellow colored herbaceous root from ginger family, is commonly used for food seasoning, for medicinal purpose and as dietary supplements [32]. Anti-inflammatory, anticarcinogenic, antioxidant, wound-healing effects are some of the medicinal value that is attributed to turmeric [47][48][49][50]. The medicinal value of turmeric is associated with the yellow color pigment "curcumin" (diferuloyol methane) content in it. Curcumin content varies in the range of 0.3% to 8.6% [47,[51][52][53][54], due to nutrients and acidity content in soil [55,56], fertilizer, soil type and cultivar [57][58][59]. Due to the isolation of curcumin from turmeric for medicinal and cosmetic purposes, it is adulterated with color dyes as a substitute for curcumin [60]. Metanil yellow (C18H14N3NaO3S), a yellow color toxic azo dye, is added to turmeric for brighter color and appealing appearance [32]. Toxicologically, metanil yellow is classified as a CII category substance by the Joint Food and Agricultural Organization of the United Nations (FAO) and World Health Organization (WHO) expert committee on Food Additives [61]. Study on rats shows that the consumption of metanil yellow causes neurotoxicity [62], hepatocellular carcinoma [63], tumor development [64], deleterious effect on gastric mucin [65], and lymphocytic leukemia [66]. This study demonstrates the application of the 1064 nm dispersive Raman spectral imaging system for the detection of metanil yellow contamination in turmeric powder and subsurface detection and identification of gelatin-encapsulated pure and mixed samples of metanil yellow and turmeric powder. For each sample, a dark current spectrum was acquired with the laser off and a cap covering the probe and subtracted from the Raman spectrum at each pixel during the spectral acquisition. Five-point moving average filter was used to smooth the spectra.

Application to Authenticate Turmeric Powder
Turmeric (Curcuma long L.), which is a yellow colored herbaceous root from ginger family, is commonly used for food seasoning, for medicinal purpose and as dietary supplements [32]. Anti-inflammatory, anticarcinogenic, antioxidant, wound-healing effects are some of the medicinal value that is attributed to turmeric [47][48][49][50]. The medicinal value of turmeric is associated with the yellow color pigment "curcumin" (diferuloyol methane) content in it. Curcumin content varies in the range of 0.3% to 8.6% [47,[51][52][53][54], due to nutrients and acidity content in soil [55,56], fertilizer, soil type and cultivar [57][58][59]. Due to the isolation of curcumin from turmeric for medicinal and cosmetic purposes, it is adulterated with color dyes as a substitute for curcumin [60]. Metanil yellow (C 18 H 14 N 3 NaO 3 S), a yellow color toxic azo dye, is added to turmeric for brighter color and appealing appearance [32]. Toxicologically, metanil yellow is classified as a CII category substance by the Joint Food and Agricultural Organization of the United Nations (FAO) and World Health Organization (WHO) expert committee on Food Additives [61]. Study on rats shows that the consumption of metanil yellow causes neurotoxicity [62], hepatocellular carcinoma [63], tumor development [64], deleterious effect on gastric mucin [65], and lymphocytic leukemia [66]. This study demonstrates the application of the 1064 nm dispersive Raman spectral imaging system for the detection of metanil yellow contamination in turmeric powder and subsurface detection and identification of gelatin-encapsulated pure and mixed samples of metanil yellow and turmeric powder. For each sample, a dark current spectrum was acquired with the laser off and a cap covering the probe and subtracted from the Raman spectrum at each pixel during the spectral acquisition. Five-point moving average filter was used to smooth the spectra.

Detection of Metanil Yellow Contamination in Turmeric Powder
Metanil yellow (70% dye, Aldrich, Carson City, NV, USA) and organic turmeric powder (Frontier Natural Products CO-OP, Norway, IA, USA) were mixed together in a vortex mixer (Scientific Industries Inc., Bohemia, NY, USA) for ten minutes to prepare sample mixture at 1%, 3%, 5%, 7%, and 10% metanil yellow concentration (w/w). The weight of each mixture sample was 0.27 g. Each mixture sample was packed into a shallow nickel-plated sample container (internal volume 25 mm × 25 mm × 1 mm), and its surface was leveled flush with the edge of the container. A Raman spectral image of each sample was obtained at 1 s exposure time and 80 mW laser power, covering a spatial area of 25 mm × 25 mm with a step size of 0.25 mm, collecting a total of 10,000 spectral pixels. Figure 4 shows the Raman spectra of turmeric powder and metanil yellow. The yellow color appearance of turmeric and metanil yellow is due to the similarly extended conjugation in both the chemicals [32]. However, a different chemical composition results in different vibrational modes specific to their chemical structure. The sharp peaks of metanil yellow at 1147 cm −1 and 1433 cm −1 , due to the N = N site are most definitive to its identification and quantification [67][68][69]. The 1631 cm −1 peak in turmeric spectrum is due to the carbonyl group. The carbonyl group is not present in metanil yellow [32]. Therefore, the 1631 cm −1 peak is most definitive for turmeric identification.

Detection of Metanil Yellow Contamination in Turmeric Powder
Metanil yellow (70% dye, Aldrich, Carson City, NV, USA) and organic turmeric powder (Frontier Natural Products CO-OP, Norway, IA, USA) were mixed together in a vortex mixer (Scientific Industries Inc., Bohemia, NY, USA) for ten minutes to prepare sample mixture at 1%, 3%, 5%, 7%, and 10% metanil yellow concentration (w/w). The weight of each mixture sample was 0.27 g. Each mixture sample was packed into a shallow nickel-plated sample container (internal volume 25 mm × 25 mm × 1 mm), and its surface was leveled flush with the edge of the container. A Raman spectral image of each sample was obtained at 1 s exposure time and 80 mW laser power, covering a spatial area of 25 mm × 25 mm with a step size of 0.25 mm, collecting a total of 10,000 spectral pixels. Figure 4 shows the Raman spectra of turmeric powder and metanil yellow. The yellow color appearance of turmeric and metanil yellow is due to the similarly extended conjugation in both the chemicals [32]. However, a different chemical composition results in different vibrational modes specific to their chemical structure. The sharp peaks of metanil yellow at 1147 cm −1 and 1433 cm −1 , due to the N = N site are most definitive to its identification and quantification [67][68][69]. The 1631 cm −1 peak in turmeric spectrum is due to the carbonyl group. The carbonyl group is not present in metanil yellow [32]. Therefore, the 1631 cm −1 peak is most definitive for turmeric identification.   Figure 5a, the dark pixels in the spectral images of turmeric show low spectral intensity of turmeric at 1147 cm −1 and 1433 cm −1 . However, spectral images of metanil yellow have white pixels revealing high spectral intensity at 1147 cm −1 and 1433 cm −1 (Figure 5a). The 1147 cm −1 and 1433 cm −1 are the highest intensity metanil yellow peaks. Figure 5b is the raw spectral images of turmeric-metanil yellow samples at different concentrations. The white pixels in the spectral images are due to the high spectral peak intensities at 1147 cm −1 and 1433 cm −1 .   Figure 5a, the dark pixels in the spectral images of turmeric show low spectral intensity of turmeric at 1147 cm −1 and 1433 cm −1 . However, spectral images of metanil yellow have white pixels revealing high spectral intensity at 1147 cm −1 and 1433 cm −1 (Figure 5a). The 1147 cm −1 and 1433 cm −1 are the highest intensity metanil yellow peaks. Figure 5b is the raw spectral images of turmeric-metanil yellow samples at different concentrations. The white pixels in the spectral images are due to the high spectral peak intensities at 1147 cm −1 and 1433 cm −1 . For identification of metanil yellow pixels, binary detection images were created. The 1147 cm −1 and 1433 cm −1 peaks were used to convert the Raman spectral images into two separate single-band binary images. An intensity threshold was used to convert all the pixels below the threshold value into background pixels of turmeric. It was observed that the intensities of some of the metanil yellow pixels were lower than the set threshold value (false negative), while the intensities of some turmeric pixels were higher than the set threshold value (false positive). To avoid the false-positive and false-negative cases, a pixel-to-pixel analysis was done to select the most appropriate threshold value. Each pixel in the binary image was evaluated with its corresponding spectrum. A final intensity threshold value of 550 was set. All the pixels below the threshold value were converted to background pixels, and the remaining pixels represented metanil yellow. Figure 6 shows the 1147 cm −1 and 1433 cm −1 binary detection images of the samples, in which the white pixels represent the detected metanil yellow particles. A gradual increase in the number of white pixels at higher concentration shows more metanil yellow particles were detected at increasing concentration. A total of 8, 19, 35, 55, and 80 metanil yellow pixels were detected in the 1147 cm −1 binary images of 1%, 3%, 5%, 7%, and 10% concentration samples ( Figure 6). A similar number of metanil yellow pixels were detected in the 1433 cm −1 binary images. The metanil yellow pixels have a similar spatial distribution in both the binary images. For identification of metanil yellow pixels, binary detection images were created. The 1147 cm −1 and 1433 cm −1 peaks were used to convert the Raman spectral images into two separate single-band binary images. An intensity threshold was used to convert all the pixels below the threshold value into background pixels of turmeric. It was observed that the intensities of some of the metanil yellow pixels were lower than the set threshold value (false negative), while the intensities of some turmeric pixels were higher than the set threshold value (false positive). To avoid the false-positive and false-negative cases, a pixel-to-pixel analysis was done to select the most appropriate threshold value. Each pixel in the binary image was evaluated with its corresponding spectrum. A final intensity threshold value of 550 was set. All the pixels below the threshold value were converted to background pixels, and the remaining pixels represented metanil yellow. Figure 6 shows the 1147 cm −1 and 1433 cm −1 binary detection images of the samples, in which the white pixels represent the detected metanil yellow particles. A gradual increase in the number of white pixels at higher concentration shows more metanil yellow particles were detected at increasing concentration. A total of 8, 19, 35, 55, and 80 metanil yellow pixels were detected in the 1147 cm −1 binary images of 1%, 3%, 5%, 7%, and 10% concentration samples ( Figure 6). A similar number of metanil yellow pixels were detected in the 1433 cm −1 binary images. The metanil yellow pixels have a similar spatial distribution in both the binary images. Out of ten thousand pixels that were acquired from 10% concentration sample, 80 detected metanil yellow pixels correspond to 0.80% of total pixels. Similarly, in 1% concentration sample, 0.08% of total pixels were detected as metanil yellow. The percentage of detected metanil yellow pixels were correlated with its actual concentration in the sample with a correlation coefficient of 0.99 (Figure 7). The result shows a linear relationship between the metanil yellow concentration in the sample and the percentage of its detected pixels.

Subsurface Detection of Turmeric Powder and Metanil Yellow
Spatially offset Raman spectroscopic (SORS) technique retrieves Raman information of the inner content of the sample from diffusely scattering media [36]. The SORS technique collects scattering signal from a series of offset distance from the laser excitation spot [70]. It uses a laser focus unit at a fixed 45° incident angle to deliver laser light onto the sample surface. An optical probe is fixed perpendicular to the sample to collect Raman scattering signal. Figure 8 shows the SORS technique to acquire Raman information of the inner content. First, Raman spectrum of sample at the no-offset position is measured by overlapping the focal point of laser and the optical probe. Then, the sample and the laser focus unit is moved at an incremental step away from the probe creating an increased offset distance between the probe and the laser focus point on the sample. As the offset distance is increased, the laser penetrates deeper through the thick top surface layer. At an increased offset distance, the signal contribution from the deep subsurface layer increases, gradually outweighing the signal contribution of the top surface layer, enabling the detection of subsurface Out of ten thousand pixels that were acquired from 10% concentration sample, 80 detected metanil yellow pixels correspond to 0.80% of total pixels. Similarly, in 1% concentration sample, 0.08% of total pixels were detected as metanil yellow. The percentage of detected metanil yellow pixels were correlated with its actual concentration in the sample with a correlation coefficient of 0.99 (Figure 7). The result shows a linear relationship between the metanil yellow concentration in the sample and the percentage of its detected pixels. Out of ten thousand pixels that were acquired from 10% concentration sample, 80 detected metanil yellow pixels correspond to 0.80% of total pixels. Similarly, in 1% concentration sample, 0.08% of total pixels were detected as metanil yellow. The percentage of detected metanil yellow pixels were correlated with its actual concentration in the sample with a correlation coefficient of 0.99 (Figure 7). The result shows a linear relationship between the metanil yellow concentration in the sample and the percentage of its detected pixels.

Subsurface Detection of Turmeric Powder and Metanil Yellow
Spatially offset Raman spectroscopic (SORS) technique retrieves Raman information of the inner content of the sample from diffusely scattering media [36]. The SORS technique collects scattering signal from a series of offset distance from the laser excitation spot [70]. It uses a laser focus unit at a fixed 45° incident angle to deliver laser light onto the sample surface. An optical probe is fixed perpendicular to the sample to collect Raman scattering signal. Figure 8 shows the SORS technique to acquire Raman information of the inner content. First, Raman spectrum of sample at the no-offset position is measured by overlapping the focal point of laser and the optical probe. Then, the sample and the laser focus unit is moved at an incremental step away from the probe creating an increased offset distance between the probe and the laser focus point on the sample. As the offset distance is increased, the laser penetrates deeper through the thick top surface layer. At an increased offset distance, the signal contribution from the deep subsurface layer increases, gradually outweighing the signal contribution of the top surface layer, enabling the detection of subsurface

Subsurface Detection of Turmeric Powder and Metanil Yellow
Spatially offset Raman spectroscopic (SORS) technique retrieves Raman information of the inner content of the sample from diffusely scattering media [36]. The SORS technique collects scattering signal from a series of offset distance from the laser excitation spot [70]. It uses a laser focus unit at a fixed 45 • incident angle to deliver laser light onto the sample surface. An optical probe is fixed perpendicular to the sample to collect Raman scattering signal. Figure 8 shows the SORS technique to acquire Raman information of the inner content. First, Raman spectrum of sample at the no-offset position is measured by overlapping the focal point of laser and the optical probe. Then, the sample and the laser focus unit is moved at an incremental step away from the probe creating an increased offset distance between the probe and the laser focus point on the sample. As the offset distance is increased, the laser penetrates deeper through the thick top surface layer. At an increased offset distance, the signal contribution from the deep subsurface layer increases, gradually outweighing the signal contribution of the top surface layer, enabling the detection of subsurface layer material. The SORS signal can be used to retrieve Raman information of multiple layers from the same sample. The SORS method has been applied for diagnosis of breast cancer by measuring calcification composition through varying tissue thickness [71]. It has been applied for evaluation of internal maturity of tomatoes by evaluating the carotenoids change in them [72]. Although, the SORS method can be used for detection and analysis of subsurface layer material, it cannot be used to measure the depth of the sample.
Empty gelatin capsules of size 000 (length 26.1 mm, single-wall thickness 0.11 mm) were obtained from Capsuline (Pompano Beach, FL, USA) to prepare capsule-turmeric samples. A one-layer capsule sample was prepared by packing turmeric tightly in a capsule and capping it. A two-layer capsule sample was prepared by inserting the turmeric packed capsule inside an empty capsule and capping the outer capsule only. Similarly, three-, four-, and five-layer capsule samples were prepared by inserting the turmeric packed capsule inside empty capsules and capping the outermost capsule. The wall thickness of the capsule layers ranged from 0.11 mm (one-layer capsule) to 0.55 mm (five-layer capsule).
Each capsule sample was held immobile in the positioning platform. The 1064 nm laser module focused a laser spot on the capsule body near the edge of its cap. First, Raman spectrum of capsule-turmeric sample was obtained. Then, the positioning table moved the sample and the laser focus unit in incremental steps in the right direction away from the Raman probe, creating an increased offset distance for each subsequent spectral acquisition. A total of 31 spectra along the longitudinal flat surface body of each capsule sample proceeding away from the capped end were collected using 20 s exposure time, 400 mW laser power, and a step size of 0.1 mm from no offset to 3 mm offset range. The spatially offset Raman spectra of samples were subtracted with the empty capsule spectrum to eliminate the influence of capsule signal in the SORS data.
Appl. Sci. 2018, 8, x 9 of 17 layer material. The SORS signal can be used to retrieve Raman information of multiple layers from the same sample. The SORS method has been applied for diagnosis of breast cancer by measuring calcification composition through varying tissue thickness [71]. It has been applied for evaluation of internal maturity of tomatoes by evaluating the carotenoids change in them [72]. Although, the SORS method can be used for detection and analysis of subsurface layer material, it cannot be used to measure the depth of the sample. Empty gelatin capsules of size 000 (length 26.1 mm, single-wall thickness 0.11 mm) were obtained from Capsuline (Pompano Beach, FL, USA) to prepare capsule-turmeric samples. A one-layer capsule sample was prepared by packing turmeric tightly in a capsule and capping it. A two-layer capsule sample was prepared by inserting the turmeric packed capsule inside an empty capsule and capping the outer capsule only. Similarly, three-, four-, and five-layer capsule samples were prepared by inserting the turmeric packed capsule inside empty capsules and capping the outermost capsule. The wall thickness of the capsule layers ranged from 0.11 mm (one-layer capsule) to 0.55 mm (five-layer capsule).
Each capsule sample was held immobile in the positioning platform. The 1064 nm laser module focused a laser spot on the capsule body near the edge of its cap. First, Raman spectrum of capsule-turmeric sample was obtained. Then, the positioning table moved the sample and the laser focus unit in incremental steps in the right direction away from the Raman probe, creating an increased offset distance for each subsequent spectral acquisition. A total of 31 spectra along the longitudinal flat surface body of each capsule sample proceeding away from the capped end were collected using 20 s exposure time, 400 mW laser power, and a step size of 0.1 mm from no offset to 3 mm offset range. The spatially offset Raman spectra of samples were subtracted with the empty capsule spectrum to eliminate the influence of capsule signal in the SORS data.       Figure 10a for two-layer capsule-turmeric sample. Similarly, for three-, four-, and five-layer capsule samples, enhancement of 1599 cm −1 and 1631 cm −1 peak is evident as the offset distance is increased (Figure 11b-d).
At no offset, the turmeric peaks evident in one-and two-layer capsule samples are attenuated as the capsule layer is increased (Figure 10c,d) because the capsule layer suppresses the deep subsurface turmeric signal. In the five-layer capsule sample (Figure 10d), 1599 cm −1 and 1631 cm −1 peaks are not evident at no offset. This prevents identification of turmeric inside the five-layer capsule. The no offset measurement measures backscattering Raman signal. This indicates that the backscattering Raman measurement technique is not useful to analyze the subsurface material through the thick top surface layer. However, increasing the offset distance enhanced the 1599 cm −1 and 1631 cm −1 peaks (Figure 10d), allowing for the identification of subsurface turmeric despite the presence of thick gelatin capsule layers.  Figure 10 shows the spatially offset Raman spectra from capsule-turmeric sample at two-, three-, four-, and five-capsule layer at no offset, 0.1 mm, 0.5 mm, 1.0 mm, 1.5 mm, 2.0 mm, 2.5 mm, and 3.0 mm offset. Steep increase in the 1599 cm −1 and 1631 cm −1 intensity with increasing offset distance is readily visible in Figure 10a for two-layer capsule-turmeric sample. Similarly, for three-, four-, and five-layer capsule samples, enhancement of 1599 cm −1 and 1631 cm −1 peak is evident as the offset distance is increased (Figure 11b-d).
At no offset, the turmeric peaks evident in one-and two-layer capsule samples are attenuated as the capsule layer is increased (Figure 10c  The effectiveness of the SORS method was further validated by identifying metanil yellow and turmeric powder packed inside one up to five-layer of capsule. Turmeric-metanil yellow mixture sample at 50% concentration (w/w) was prepared and packed inside gelatin capsule to obtain one-, two-, three-, four-, and five-layer samples. A series of spatially offset Raman spectra were collected from no offset to 3 mm offset using a 0.1 mm increment along the longitudinal flat surface of capsule sample. Figure 11 shows the SORS spectra from the one-layer capsule sample. The 994 cm −1 , 1147 cm −1 , 1189 cm −1 , 1433 cm −1 , and 1602 cm −1 are metanil yellow peaks. The 1631 cm −1 is a turmeric peak. The turmeric peak at 1599 cm −1 is not resolved from the 1602 cm −1 peak of metanil yellow. Metanil yellow peaks at 1147 cm −1 , 1189 cm −1 , and 1433 cm −1 are evident at no offset, and are enhanced by increasing the offset distance. The 994 cm −1 metanil yellow peak is not evident at no offset. Increasing The effectiveness of the SORS method was further validated by identifying metanil yellow and turmeric powder packed inside one up to five-layer of capsule. Turmeric-metanil yellow mixture sample at 50% concentration (w/w) was prepared and packed inside gelatin capsule to obtain one-, two-, three-, four-, and five-layer samples. A series of spatially offset Raman spectra were collected from no offset to 3 mm offset using a 0.1 mm increment along the longitudinal flat surface of capsule sample. Figure 11 shows the SORS spectra from the one-layer capsule sample. The 994 cm −1 , 1147 cm −1 , 1189 cm −1 , 1433 cm −1 , and 1602 cm −1 are metanil yellow peaks. The 1631 cm −1 is a turmeric peak. The turmeric peak at 1599 cm −1 is not resolved from the 1602 cm −1 peak of metanil yellow. Metanil yellow peaks at 1147 cm −1 , 1189 cm −1 , and 1433 cm −1 are evident at no offset, and are enhanced by increasing the offset distance. The 994 cm −1 metanil yellow peak is not evident at no offset. Increasing the offset distance enhanced the 994 cm −1 peak. The 1631 cm −1 turmeric peak is not evident at no offset. Increasing the offset distance gradually increased the 1631 cm −1 peak intensity.
the offset distance enhanced the 994 cm −1 peak. The 1631 cm −1 turmeric peak is not evident at no offset. Increasing the offset distance gradually increased the 1631 cm −1 peak intensity. Figure 11. Spatially offset Raman spectra of turmeric-metanil yellow mixture packaged inside one-layer capsule. Figure 12 shows the SORS spectra from two-, three-, four-, and five-layer capsule samples at no offset, 0.1 mm, 0.5 mm, 1.0 mm, 1.5 mm, 2.0 mm, 2.5 mm, and 3.0 mm offset. Increasing the capsule layer attenuated the metanil yellow and turmeric peaks at no offset. The metanil yellow peaks at 1147 cm −1 and 1433 cm −1 were evident at no offset in the two-and three-layer capsule samples ( Figure  12a,b). The 1147 cm −1 peak is not evident at no offset in the four-and five-layer capsule samples (Figure 12c,d). Increasing the offset distance enhanced the 1147 cm −1 peak (Figure 12c,d). The 1433 cm −1 peak is evident at no offset in the two-, three-, and four-layer capsule samples (Figure 12a-c). An increase in the peak intensity at 1433 cm −1 was observed from no offset to 1.0 mm offset in two-, three-, and four-layer capsule samples (Figure 12a-c). After 1.0 mm offset, the 1433 cm −1 peak intensity gradually decreased. In the five-layer capsule sample, the 1433 cm −1 peak is not evident at no offset (Figure 12d). Increasing the offset distance enhanced the 1433 cm −1 peak in the five-layer capsule sample.
The turmeric peak at 1631 cm −1 is not evident at no offset in the two-, three-, four-, and five-layer capsule samples (Figure 12a-d). Increasing the offset distance enhanced the 1631 cm −1 peak in the multiple layer capsule samples. The 1631 cm −1 peak can be used to identify the turmeric packed inside gelatin capsule. The results demonstrate this SORS measurement technique has a potential for non-destructive subsurface detection and the identification of gelatin-encapsulated pure and mixed samples of turmeric and metanil yellow powder. Figure 11. Spatially offset Raman spectra of turmeric-metanil yellow mixture packaged inside one-layer capsule. Figure 12 shows the SORS spectra from two-, three-, four-, and five-layer capsule samples at no offset, 0.1 mm, 0.5 mm, 1.0 mm, 1.5 mm, 2.0 mm, 2.5 mm, and 3.0 mm offset. Increasing the capsule layer attenuated the metanil yellow and turmeric peaks at no offset. The metanil yellow peaks at 1147 cm −1 and 1433 cm −1 were evident at no offset in the two-and three-layer capsule samples (Figure 12a,b). The 1147 cm −1 peak is not evident at no offset in the four-and five-layer capsule samples (Figure 12c,d). Increasing the offset distance enhanced the 1147 cm −1 peak (Figure 12c,d). The 1433 cm −1 peak is evident at no offset in the two-, three-, and four-layer capsule samples (Figure 12a-c). An increase in the peak intensity at 1433 cm −1 was observed from no offset to 1.0 mm offset in two-, three-, and four-layer capsule samples (Figure 12a-c). After 1.0 mm offset, the 1433 cm −1 peak intensity gradually decreased. In the five-layer capsule sample, the 1433 cm −1 peak is not evident at no offset (Figure 12d). Increasing the offset distance enhanced the 1433 cm −1 peak in the five-layer capsule sample.
The turmeric peak at 1631 cm −1 is not evident at no offset in the two-, three-, four-, and five-layer capsule samples (Figure 12a-d). Increasing the offset distance enhanced the 1631 cm −1 peak in the multiple layer capsule samples. The 1631 cm −1 peak can be used to identify the turmeric packed inside gelatin capsule. The results demonstrate this SORS measurement technique has a potential for non-destructive subsurface detection and the identification of gelatin-encapsulated pure and mixed samples of turmeric and metanil yellow powder.

Conclusions
A 1064 nm dispersive Raman spectral imaging system was developed for food safety and quality evaluation. The system is capable of hyperspectral Raman imaging and spatially offset Raman spectroscopic (SORS) measurement of the food samples. The application of the hyperspectral Raman imaging system was demonstrated by detecting metanil yellow contamination in turmeric powder at five concentration levels. A hyperspectral Raman image of the mixture sample at each concentration was collected covering the sample surface area of 25 mm × 25 mm, using a step size of 0.25 mm along the X and Y direction. The system could detect metanil yellow at concentrations as low as 1%. The detected metanil yellow pixels linearly correlated with its concentration in the

Conclusions
A 1064 nm dispersive Raman spectral imaging system was developed for food safety and quality evaluation. The system is capable of hyperspectral Raman imaging and spatially offset Raman spectroscopic (SORS) measurement of the food samples. The application of the hyperspectral Raman imaging system was demonstrated by detecting metanil yellow contamination in turmeric powder at five concentration levels. A hyperspectral Raman image of the mixture sample at each concentration was collected covering the sample surface area of 25 mm × 25 mm, using a step size of 0.25 mm along the X and Y direction. The system could detect metanil yellow at concentrations as low as 1%. The detected metanil yellow pixels linearly correlated with its concentration in the mixture sample (R 2 = 0.99). Application of the system for SORS measurement was demonstrated by subsurface detection of gelatin-encapsulated pure and mixed samples of turmeric and metanil yellow powder. This system is a versatile platform capable of 1064 nm hyperspectral Raman imaging and SORS measurement for a wide variety of food and pharmaceutical samples.