Next Article in Journal / Special Issue
Recycling-Oriented Characterization of Post-Earthquake Building Waste by Different Sensing Techniques
Previous Article in Journal
Incremental Learning for Dermatological Imaging Modality Classification
Previous Article in Special Issue
Low-Cost Hyperspectral Imaging with A Smartphone
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effective Recycling Solutions for the Production of High-Quality PET Flakes Based on Hyperspectral Imaging and Variable Selection

Department of Chemical Engineering, Materials & Environment, Sapienza, Rome University, Via Eudossiana 18, 00184 Rome, Italy
*
Author to whom correspondence should be addressed.
Submission received: 30 July 2021 / Revised: 1 September 2021 / Accepted: 6 September 2021 / Published: 8 September 2021
(This article belongs to the Special Issue Hyperspectral Imaging and Its Applications)

Abstract

:
In this study, effective solutions for polyethylene terephthalate (PET) recycling based on hyperspectral imaging (HSI) coupled with variable selection method, were developed and optimized. Hyperspectral images of post-consumer plastic flakes, composed by PET and small quantities of other polymers, considered as contaminants, were acquired in the short-wave infrared range (SWIR: 1000–2500 nm). Different combinations of preprocessing sets coupled with a variable selection method, called competitive adaptive reweighted sampling (CARS), were applied to reduce the number of spectral bands useful to detect the contaminants in the PET flow stream. Prediction models based on partial least squares-discriminant analysis (PLS-DA) for each preprocessing set, combined with CARS, were built and compared to evaluate their efficiency results. The best performance result was obtained by a PLS-DA model using multiplicative scatter correction + derivative + mean center preprocessing set and selecting only 14 wavelengths out of 240. Sensitivity and specificity values in calibration, cross-validation and prediction phases ranged from 0.986 to 0.998. HSI combined with CARS method can represent a valid tool for identification of plastic contaminants in a PET flakes stream increasing the processing speed as requested by sensor-based sorting devices working at industrial level.

1. Introduction

Plastics represent one of the most used materials, in daily life, in a wide range of applications, due to their peculiar characteristics and low production costs [1]. As a consequence, there has been an uncontrolled growth of large quantities of plastic waste [2], especially from packaging, still creating a series of challenges for industrialized countries at a political, economic, social, and environmental level [3]. In order to achieve circular economy and recycling targets, set by European and national legislation, to prevent the environmental impacts of plastic packaging waste, it is essential to implement efficient plastic waste recovery strategies [3,4,5]. Several actions can be taken to improve plastic recycling processes, thus allowing to bring high-quality recycled products to the market. In this context, the on-line sorting step of the mechanical recycling process plays a preeminent role in order to improve processing performance, increasing recycled plastic quality. Contaminants, i.e., other materials and other types of polymers, inside the post-consumer stream of a specific recycled polymer, can degrade the final properties of the secondary raw material [6,7,8]. A correct recognition and separation of materials in recycling plants is, thus, crucial.
Optical-based sorting of polymers is one of the key points in order to produce high-quality plastics as secondary raw materials [9,10]. Many spectroscopy- based approaches can be applied for plastic classification, including near-infrared (NIR) and short-wave infrared (SWIR) spectroscopy [10,11], Raman spectroscopy [12], and laser induced breakdown spectroscopy (LIBS) [13]. Compared with other spectroscopic techniques, NIR/SWIR spectroscopy has the advantage of rapid detection, little sample preparation and low cost [14]. NIR/SWIR spectroscopy can, thus, be profitably utilized to perform automatic plastic sorting or to implement quality control strategies in recycling plants [14]. In fact, useful information about polymers can be obtained in the NIR/SWIR range, as their molecules absorb light by overtone or combination vibrations [10].
Hyperspectral imaging (HSI), operating in the NIR or SWIR range, represents an attractive solution to characterize, identify, and classify various materials, thanks to its ability to provide information on spectral features and spatial distribution [15,16,17,18,19,20,21]. The selection of SWIR range, being characterized by a greater number of wavelengths, can allow to recognize many different polymers, even with slight spectral differences, reducing the misclassification errors [10,22,23,24,25]. HSI is based on the use of an integrated hardware and software architecture able to acquire and process data, useful to obtain spatial and spectral information of the investigated object. This information is contained in a three-dimensional dataset (i.e., two spatial dimensions and a spectral dimension), the so-called “hypercube” [9,10]. As the hyperspectral image is characterized by high-dimensional data, its spectral information is often affected by multicollinearity [26] and requires some data processing time. In fact, not all variables (i.e., wavelengths) are useful (i.e., presence of noise) or necessary (i.e., redundancy of information) to build the prediction models [27,28]. Thus, a variables selection approach could optimize classification logics of plastic waste, improving the efficiency and speed of optical sorting machines working at industrial level. The goal of variable selection is to obtain a small set of variables offering the best or at least comparable generalization or simplification capacity compared to the original set of variables [29,30]. Therefore, the variable selection can play an important role in HSI analysis before modeling in order to extract the most relevant and sensitive information.
The proposed study was carried out to develop efficient strategies for sensor-based sorting of plastic waste in recycling plants, based on hyperspectral imaging (HSI) and variable selection approach, in particular for the production of a high-quality recycled polyethylene terephthalate (PET) flakes stream. The use of variable selection methods to reduce processing time in optical sensors-based sorting is a current goal [31,32]. In more detail, the purpose is to reduce the number of variables to be used without losing quality in the recognition during the sorting process. For this reason, the use of new techniques for variable selection can be profitable improve the quality of sorting process.
In this study, HSI based analysis, working in the SWIR range (1000–2500 nm), was applied as a fast and non-destructive detection technique useful to obtain high predictive results. Different preprocessing strategies, among the most used in literature on NIR data [33,34,35,36,37,38,39,40,41,42,43], were evaluated to select the most efficient set or sequence, able to emphasize PET and contaminant spectral differences.
Competitive adaptive reweighted sampling (CARS) method, which performs by simulating the Darwinian “survival of the fittest” theory of evolution [44,45] was applied to eliminate the useless or irrelevant variables, and to select an optimal combination of effective wavelengths useful to recognize contaminants in a PET stream. Partial least squares discriminant analysis (PLS-DA) models were built for each preprocessing set, in order to evaluate the one showing the best efficiency to identify classes of polymers, i.e., PET and other polymers considered as a single class of contaminants.

2. Materials and Methods

2.1. Samples Overview

Plastic samples, representative of a flow-stream of PET flakes contaminated by other polymers, were collected from a recycling plant. In this scenario, the contaminants have limited and finite variability sources [46], allowing the possibility to create a representative prediction model with defined wavelengths.
Plastic flakes of PET and other polymers were selected and divided into calibration and prediction datasets for the evaluation of the PLS-DA models (Figure 1). In detail, the calibration dataset (CAL) was created from an individual image containing 36 samples divided into 18 PET and 18 contaminant flakes (Figure 1a). Principal component analysis (PCA) was used to set classes and defining the calibration set. The calibration dataset was pre-processed and cross-validated (CV) for building a PLS-DA model to detect the presence of contaminants on PET stream. The prediction image (PRED) was created from a set of plastic samples external to the model, characterized by 18 PET and 18 flakes of contaminants randomly arranged (Figure 1b).

2.2. Data Acquisition and Analysis

Hyperspectral images acquisition was carried out at the Raw Materials Laboratory (RawMaLab) of the Department of Chemical Engineering, Materials and Environment of Sapienza University of Rome by the Sisuchema XLTM Chemical Imaging Workstation (Specim Ltd., Oulu, Finland) (Figure 2). The HSI platform is based on a push-broom acquisition architecture, with a camera operating from 1000 to 2500 nm (SWIR range). The selected configuration of the device covers a maximum field of view of 20 cm with a pixel resolution of 625 µm. The HSI platform is equipped with a diffuse line illumination unit, consisting of quartz halogen lamps producing dual linear light, covering a spectrum range of 920 to 2514 nm, thus optimizing the imaging of various surfaces [47]. The working distance between the spectrograph lens and the sample tray plan was 30 cm. The device technical specifications are summarized in Table 1. Reflectance of hypercube was automatically set up by an internal standard reference target. A total of 240 wavelengths were collected and analyzed for each dataset. The number of pixels collected for the calibration dataset was 3420 for the PET class and 2089 for the contaminant class, while the number of pixels of the prediction dataset was 3594 for the PET class and 2836 for the contaminant class.
PLS_toolbox (ver. 8.8 Eigenvector Research, Inc., Wenatchee, WA, USA) running in the Matlab environment (version R2020a, The Mathworks, Inc., Nertick, MA, USA) was used to analyze the acquired hyperspectral images.

2.3. Data Preprocessing

Different preprocessing strategies, according to the most applied to infrared spectral data [33,34,35,36,37,38,39,40,41,42,43], including those related to plastic samples [10,15,17,48,49,50,51,52,53], were selected to build each pretreatment sequence, that is:
  • Standard Normal Variate (SNV): SNV was applied to reduce the scattering effects in the spectral data and to obtain a general linearization of the relationship between signal and concentration [10,33,34,35,36,37,38,39];
  • Savitzky–Golay (SG) derivative: Derivatives are a common method used to remove unimportant baseline signal from data. SG first derivative filter was applied to emphasize the spectral differences with second polynomial order and 15 points window [10,33,34,35,36,37,40];
  • Multiplicative Scatter Correction (MSC): MSC is widely used for infrared data (such as SNV and derivation). MSC was useful to remove artifacts or imperfections from data, such as undesirable scatter effect [32,33,34,35,36,37];
  • Smoothing: Smoothing (Savitzky–Golay routine) was used as low-pass filter (15 points) for removing high-frequency noise due to the derivation process [10,40,41,42,43];
  • Detrend: Detrend was applied on spectra to remove the effects of baseline shift and curvilinearity [34,35,39];
  • Mean Center (MC): Centering is one of the most common types of preprocessing, usually applied. MC has the effect to include an adjustable intercept in multivariate models [10,38,41,42,43].

2.4. Principal Component Analysis (PCA)

PCA is often applied for HSI data exploration, useful to provide an overview of multivariate data and to evaluate the selected preprocessing combinations [54]. PCA allows the decomposition of preprocessed spectral data into linear combinations of the original spectral data, called principal components (PCs), collecting the spectral variations in reduced set of factors. The first PCs were used to analyze the common characteristics of samples and their grouping, as the samples characterized by similar spectral signatures tend to aggregate in the score plot of the first two or three components [54].

2.5. Competitive Adaptive Reweighted Sampling (CARS)

CARS is an innovative and useful wavelength selection approach [44] used in NIR spectroscopy to select variables (i.e., significant wavelengths) [45]. CARS has the potential to select an optimal combination of the useful wavelengths from the full spectrum, combined with PLS regression [54]. In CARS method, regression coefficients (RC) absolute values of PLS model are used to evaluate the weight of each wavelength. Based on the importance of each wavelength, CARS sequentially selects N subsets of wavelengths by N Monte Carlo sampling run in an iterative and competitive manner. First, in each sampling run, samples are randomly selected in a fixed ratio (e.g., 80%) to build a calibration model. Then, based on RC, exponentially decreasing function (EDF) and adaptive reweighted sampling (ARS) procedures are applied to select the key wavelengths. Finally, the subset with the lowest root mean-square error of cross validation (RMSECV) is chosen.

2.6. Partial Least Square Discriminant Analysis (PLS-DA)

PLS-DA was used to identify predefined classes of materials (i.e., PET and contaminants), by forming discriminant functions from input variables (i.e., wavelengths) to produce a new set of transformed values useful to provide a more accurate discrimination than any single variable [55]. Venetian blind (number of data splits = 10) as cross-validation method was used, in order to evaluate the complexity of the models and to select the appropriate number of latent variables (LVs) (Figure S1). The optimal number of LVs was also determined by the smaller difference between RMSEC and RMSECV [56,57,58].

PLS-DA Performances

The classification performances obtained by PLS-DA models were evaluated in terms of statistical parameters: sensitivity, specificity and efficiency (Equations (1), (2), and (3)).
Sensitivity = True   Positive ( True   Positive   +   False   Negative )
Specificity = True   Negative ( True   Negative   +   False   Positive )
Efficiency = ( sensitivity × specificity )  

3. Experimental Results and Discussion

3.1. Average Raw Reflectance Spectra

The average raw reflectance spectra and the standard deviation of the two classes of polymers are shown in Figure 3. PET average spectrum was characterized by absorption bands of C-H2 and C-H of the third harmonic region (1131 and 1182 nm), C-H of the second harmonic region (1402, 1665, 1723, 1825, 1910, and 1960 nm) and C-H stretching vibrations + C-H deformation of first combination region (2090, 2136, 2161, 2186, and 2261 nm). The average spectrum of contaminants showed a complex fingerprint due to the presence of different types of polymers, with the main absorption bands located around 1220, 1402, 1735, and 2317 nm.

3.2. Preprocessing Sets and Variables Selection

In order to emphasize the spectral differences between PET and contaminants, three sets of preprocessing techniques were used, that is:
  • Set 1: Detrend + Smoothing + MC;
  • Set 2: SNV + MC;
  • Set 3: MSC + Derivative + MC.
The average reflectance spectra of the PET and contaminant classes resulting from the application of the aforementioned preprocessing sets, are shown in Figure 4.
Subsequently, for each preprocessing set, the CARS method [44] was applied, in order to reduce the number of wavelengths useful to discriminate the spectral characteristics between PET and contaminants. The selected wavelengths for each adopted preprocessing sets are shown in Table 2.

3.3. PCA Results of Preprocessing Set 1 (Detrend + Smoothing + MC)

PCA score and loadings plots are shown in Figure 5. Most of the variance was captured by the first two PCs, as shown in the score plot (Figure 5a), where PC1 and PC2 explained about the 74.09% and 21.07% of the variance, respectively. The PCA score plot showed two clouds corresponding to the two analyzed classes (i.e., PET and contaminant). The cluster separation was acceptable with a low overlapping of clouds in the fourth quadrant. In more detail, the PET scores, due to the low spectral variance and high uniformity detected in PET samples, were more grouped than the contaminant scores. The variance of the contaminant was greater than the PET class, as it was influenced by the spectral combination of different polymers, as shown in the PCA score plot. The loadings plot of PC1 and PC2 was shown in Figure 5b. The main PC1 variance was given by the wavelengths around 1240 and 1720 nm for positive values, while for negative values it was mainly given by the wavelengths around 1320 and 2100 nm. PC2 was mostly influenced by wavelengths around 1730, 1910, and 2100 nm for positive values, whereas negative values were highlighted by wavelengths about 1005 and 2450 nm.

3.4. PCA Results of the Preprocessing Set 2 (SNV + MC)

PCA score and loadings plots are shown in Figure 6. The PCA model showed a captured variance of 95.07% with 3 PCs. The best separation between PET and contaminant clusters was allowed by PC1 vs. PC2, as shown in the PCA score plot (Figure 6a). Cluster separation was very noticeable with few overlapping pixels compared to the previous preprocessing sets. PET cluster was mainly located in the first and second quadrant, while the distribution of the contaminant scores was mainly localized in the third and fourth quadrant. Both clusters showed a similar variance distribution. Therefore, the preprocessing Set 2 (SNV + MC) approach allowed to minimize the intra-class variance, emphasizing the differences between PET and contaminant classes. The loadings plots of PC1 and PC2 are shown in Figure 6b. The PC1 variance was mainly given for positive values by the wavelengths around 1720, 2250, and 2480 nm, and for negative values by the wavelengths around 1020, 1130, and 1326 nm. PC2 was mostly marked for positive values by wavelengths around 1206, 1720, and 1920 nm, and for negative values by wavelengths around 1650 and 2255 nm.

3.5. PCA Results of Preprocessing Set 3 (MSC + Derivative + MC)

PCA scores and loadings plots are shown in Figure 7. The PCA model showed a captured variance of 95.59% with 3 PCs. The best separation between PET and contaminants was allowed by PC1 vs. PC2. The PCA score plot showed two clusters related to PET and contaminant classes. The score plot showed a cluster separation and a low cluster overlap in the central zone (Figure 7a). PET cluster was mainly located in the first quadrant, while the class of contaminants was mainly localized in the second quadrant. Therefore, preprocessing Set 3 (MSC + Derivative + MC) allowed to minimize the intra-class variance, and to preserve the spectral differences between the two classes. The loadings plot of PC1 and PC2 is shown in Figure 7b. The PC1 variance was mainly given, for positive values, by the wavelengths around 2274 nm, and for negative values by the wavelengths around 2300 nm. PC2 was mainly influenced for positive values by wavelengths around 1050 nm, and for negative values by wavelengths around 1060 nm.

3.6. Classification Performances

3.6.1. PLS-DA Models Constructed for a Limited Set of Spectral Variables

Starting from the characteristics detected by the PCA of each preprocessing set with selected variables, a PLS-DA model was constructed. The correct number of LVs was chosen based on the smallest difference between the root mean square error for calibration (RMSEC) and cross-validation (RMSECV) values (Table 3). PLS-DA model of preprocessing Set 1 (Detrend + Smoothing + MC) showed a variance captured of 96.66% with 4 LVs; the PLS-DA model of preprocessing Set 2 (SNV + MC) showed a variance captured of 98.64% with 3 LVs and, finally, PLS-DA model of preprocessing Set 3 (MSC + Derivative + MC) showed a variance captured of 93.68% with 3 LVs. The PLS-DA models prediction results are shown in Figure 8. In general, in all models, PET and contaminant flakes samples were properly recognized using the PLS-DA model, with the presence of a few pixels not correctly classified. The only exception was related to the results achieved based on preprocessing Set 1 (Detrend + Smoothing + MC), showing a sample with multiple misclassified pixels (highlighted with a yellow circle in Figure 8) and border-effect in some flakes. The results obtained by PLS-DA models related to preprocessing Set 2 (SNV + MC) and 3 (MSC + Derivative + MC) showed a similar prediction quality, with few misclassification pixels mainly due to border-effect. However, the few pixels not correctly assigned, do not significantly affect the correct class recognition.
The classification performances obtained by the different preprocessing sets, shown in Table 4, revealed sensitivity and specificity values in calibration, cross-validation, and prediction ranging from 0.957 to 0.999. Efficiency values in prediction ranges from 0.969 to 0.991 confirming the positive quality of all PLS-DA models combined with variables selection. Based on the measured performance parameters, Set 2 (SNV + MC) and Set 3 (MSC + Derivative + MC) show a similar result in terms of specificity, sensitivity, and efficiency. However, Set 3 (MSC + Derivative + MC) is better because it uses a smaller number of wavelengths (14 for Set 3 vs. 29 for Set 2). Therefore, the best PLS-DA model was the one obtained starting from preprocessing Set 3 (MSC + Derivative + MC).

3.6.2. Comparison of Full Spectrum and Reduced Wavelength PLS-DA with Preprocessing Set 3 (MSC + Derivative + MC)

Finally, the performances of full spectrum PLS-DA using preprocessing Set 3 (MSC + Derivative + MC) were compared with those obtained in variables selection mode with the same preprocessing set. In details, the full spectrum PLS-DA model showed a captured variance of 99.39% with 5 LVs. The LVs number was chosen based on the smallest difference between the RMSEC and RMSECV values (Table 5). Full spectrum PLS-DA prediction results are shown in Figure 9. In particular, PET and contaminant classes were well predicted, with sensitivity and specificity values in calibration, cross-validation and prediction phases and efficiency (Table 6) ranging from 0.986 to 1.000 for both classes.
The comparison of the prediction results based on the PLS-DA-Set 3 (MSC + Derivative + MC) applied to the full spectrum hypercubes (Figure 9) and to the 14 selected wavelengths (Figure 8c) showed as they are similar. In detail, analyzing the sensitivity, specificity and efficiency values of the two models, it can be noticed a slight increase in misclassified pixel/spectra in the PLS-DA in variable selection model. However, the misclassified pixels were mainly located along the boundary of the samples, not affecting the correct attribution of the class.

4. Economic and Environmental Impact

The systematic implementation of the HSI detection and classification-based logic could have important effects both at commercial-industrial and at economic-environmental level. The proposed approach can produce not only a better separation efficiency, but also a product of better quality. The fulfilment of these two goals generates social, economic, and environmental benefits [59]. In fact, in an economic viability context, a stronger and widespread PET recycling sector generates employment and contributes to reduce the volume of municipal solid waste [60]. In addition, high-quality recycled PET contributes to reduce the consumption of energy and non-renewable raw materials, [61], according to the sustainable development goals (SDGs) of UN Agenda 2030, and in particular to SDG 12, and to the principles of circular economy.

5. Conclusions

The application of HSI in the SWIR region was investigated to evaluate the feasibility of a rapid and non-destructive method for the identification of plastic contaminants in a recycled PET flakes stream, producing a high-quality secondary raw material. CARS was tested as variable selection method after the application of three different preprocessing sequences to identify the best combination for the recognition of contaminants in PET stream with a limited number of wavelengths. The results of the variable selection obtained by CARS were evaluated by a PLS-DA model for each set of selected wavelengths. The best prediction results in calibration and cross-validation were provided by the combination of CARS and the preprocessing Set 3 (MSC + Derivative + MC), reducing the spectral dataset from 240 to 14 wavelengths. In addition, a comparison was made between the performances of the full spectrum PLS-DA model using preprocessing Set 3 (MSC + Derivative + MC) and those obtained in variable selection mode with the same preprocessing set. The results demonstrated that the correctness of the classification was similar, further highlighting the possibility to identify plastic contaminants in the recycled PET flakes stream using a limited number of key wavelengths, useful for online sorting applications.
The current study supplied an effective procedure for variable selection from hyperspectral images, reducing data redundancy and obtaining a prediction efficiency close to that obtained by the full spectrum PLS-DA model. The obtained results enable the possibility to build a multispectral detection system based on filters analyzing selected spectral regions, with a significant reduction in costs compared to a conventional full spectrum hyperspectral camera and ensuring a high quality of recycled PET stream.

Supplementary Materials

The following is available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/jimaging7090181/s1, Figure S1: image of cross-validation.

Author Contributions

Conceptualization G.B. and S.S.; methodology, G.B., G.C., P.C. and S.S.; software, G.C. and P.C.; validation, G.C. and P.C.; formal analysis, G.C. and P.C.; investigation, G.C. and P.C.; resources, G.B., S.S.; data curation, P.C. and G.C.; writing—original draft preparation, P.C. and G.C.; writing—review and editing, G.B. and S.S.; visualization, P.C. and G.C.; supervision, G.B. and S.S.; project administration, G.B. and S.S.; funding acquisition, G.B. and S.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors thank Antonio Uzzo of 3UVision Srl who provided the samples.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ragaert, K.; Delva, L.; Van Geem, K. Mechanical and chemical recycling of solid plastic waste. Waste Manag. 2017, 69, 24–58. [Google Scholar] [CrossRef]
  2. Vilaplana, F.; Karlsson, S. Quality concepts for the improved use of recycled polymeric materials: A review. Macromol. Mater. Eng. 2008, 293, 274–297. [Google Scholar] [CrossRef]
  3. Plastics Europe. The Circular Economy for Plastics—A European overview. 2019. Available online: https://www.plasticseurope.org/it/resources/publications/1899-circular-economy-plastics-european-overview (accessed on 23 June 2021).
  4. Schroeder, P.; Anggraeni, K.; Weber, U. The Relevance of Circular Economy Practices to the Sustainable Development Goals. J. Ind. Ecol. 2019, 23, 77–95. [Google Scholar] [CrossRef] [Green Version]
  5. Ellen MacArthur Foundation. The New Plastics Economy: Rethinking the Future of Plastics. Report Produced by World Economic Forum and Ellen MacArthur Foundation. 2016. Available online: https://www.ellenmacarthurfoundation.org/publications/the-new-plastics-economy-rethinking-the-future-of-plastics (accessed on 23 June 2021).
  6. Alsewailem, F.D.; Alrefaie, J.K. Effect of contaminants and processing regime on the mechanical properties and moldability of postconsumer polyethylene terephthalate bottles. Waste Manag. 2018, 81, 88–93. [Google Scholar] [CrossRef]
  7. Schyns, Z.O.G.; Shaver, M.P. Mechanical Recycling of Packaging Plastics: A Review. Macromol. Rapid Commun. 2021, 42, 2000415. [Google Scholar] [CrossRef]
  8. Lahtela, V.; Kärki, T. Mechanical Sorting Processing of Waste Material Before Composite Manufacturing—A Review. J. Eng. Sci. Technol. Rev. 2018, 11, 35–46. [Google Scholar] [CrossRef]
  9. Serranti, S.; Gargiulo, A.; Bonifazi, G. Hyperspectral imaging for process and quality control in recycling plants of polyolefin flakes. J. Near Infrared Spectrosc. 2012, 20, 573–581. [Google Scholar] [CrossRef]
  10. Bonifazi, G.; Capobianco, G.; Serranti, S. A hierarchical classification approach for recognition of low-density (LDPE) and high-density polyethylene (HDPE) in mixed plastic waste based on short-wave infrared (SWIR) hyperspectral imaging. Spectrochim. Acta Part A Mol. Biomol. Spectrosc. 2018, 198, 115–122. [Google Scholar] [CrossRef]
  11. Serranti, S.; Cucuzza, P.; Bonifazi, G. Hyperspectral imaging for VIS-SWIR classification of post-consumer plastic packaging products by polymer and color. In Proceedings of the SPIE 11525, SPIE Future Sensing Technologies, Online, 13 November 2020; Volume 1152510. [Google Scholar] [CrossRef]
  12. Allen, V.; Kalivas, J.H.; Rodriguez, R.G. Post-Consumer Plastic Identification Using Raman Spectroscopy. Appl. Spectrosc. 1999, 53, 672–681. [Google Scholar] [CrossRef]
  13. Gondal, M.A.; Siddiqui, M.N. Identification of different kinds of plastics using laser-induced breakdown spectroscopy for waste management. J. Environ. Sci. Health Part A 2007, 42, 1989–1997. [Google Scholar] [CrossRef] [PubMed]
  14. Wu, X.; Li, J.; Yao, L.; Xu, Z. Auto-sorting commonly recovered plastics from waste household appliances and electronics using near-infrared spectroscopy. J. Clean. Prod. 2020, 246, 11873. [Google Scholar] [CrossRef]
  15. Serranti, S.; Fiore, F.; Bonifazi, G.; Takeshima, A.; Takeuchi, H.; Kashiwada, S. Microplastics characterization by hyperspectral imaging in the SWIR range. In Proceedings of the SPIE 11197, SPIE Future Sensing Technologies, Online, 12 November 2019; Volume 1119710. [Google Scholar] [CrossRef]
  16. Hibbitts, C.A.; Bekker, D.; Hanson, T.; Knuth, A.; Goldberg, A.; Ryan, K.; Cantillo, D.; Daubon, D.; Morgan, F. Dual-band discrimination and imaging of plastic objects. In Proceedings of the SPIE 11012, Detection and Sensing of Mines, Explosive Objects, and Obscured Targets XXIV, Online, 22 May 2019; Volume 1101211. [Google Scholar] [CrossRef]
  17. Serranti, S.; Gargiulo, A.; Bonifazi, G. Characterization of post-consumer polyolefin wastes by hyperspectral imaging for quality control in recycling processes. Waste Manag. 2011, 31, 2217–2227. [Google Scholar] [CrossRef]
  18. Mehrubeoglu, M.; Zemlan, M.; Henry, S. Hyperspectral Imaging for Differentiation of Foreign Materials from Pinto Beans; Paper 96110A; SPIE: San Diego, CA, USA, 2015; Volume 9611. [Google Scholar] [CrossRef]
  19. Ferrari, C.; Foca, G.; Calvini, R.; Ulrici, A. Fast exploration and classification of large hyperspectral image datasets for early bruise detection on apples. Chemom. Intell. Lab. Syst. 2015, 146, 108–119. [Google Scholar] [CrossRef] [Green Version]
  20. Amigo, J.M.; Martì, I.; Gowen, A. Hyperspectral Imaging and Chemometrics. A Perfect Combination for the Analysis of Food Structure, Composition and Quality. Data Handl. Sci. Technol. 2013, 28, 343–370. [Google Scholar] [CrossRef]
  21. Ulrici, A.; Serranti, S.; Ferrari, C.; Cesare, D.; Foca, G.; Bonifazi, G. Efficient chemometric strategies for PET-PLA discrimination in recycling plants using hyperspectral imaging. Chemometr. Intell. Lab. 2013, 122, 31–39. [Google Scholar] [CrossRef]
  22. Singh, N.; Hui, D.; Singh, R.; Ahuja, I.P.S.; Feo, L.; Fraternali, F. Recycling of plastic solid waste: A state of art review and future applications. Compos. Part B Eng. 2017, 115, 409–422. [Google Scholar] [CrossRef]
  23. Caballero, D.; Bevilacqua, M.; Amigo, J.M. Application of hyperspectral imaging and chemometrics for classifying plastics with brominated flame retardants. J. Spectr. Imaging. 2019, 8. [Google Scholar] [CrossRef] [Green Version]
  24. Lorenzo-Navarro, J.; Serranti, S.; Bonifazi, G.; Capobianco, G. Performance Evaluation of Classical Classifiers and Deep Learning Approaches for Polymers Classification Based on Hyperspectral Images. In Advances in Computational Intelligence IWANN 2021. Lecture Notes in Computer Science; Rojas, I., Joya, G., Catala, A., Eds.; Springer: Cham, Switzerland, 2021; Volume 12862. [Google Scholar] [CrossRef]
  25. Serranti, S.; Palmieri, R.; Bonifazi, G.; Cózar, A. Characterization of microplastic litter from oceans by an innovative approach based on hyperspectral imaging. Waste Manag. 2018, 76, 117–125. [Google Scholar] [CrossRef]
  26. Yun, Y.H.; Li, H.D.; Deng, B.C.; Cao, D.S. An overview of variable selection methods in multivariate analysis of near-infrared spectra. Trends Anal. Chem. 2019, 113, 102–115. [Google Scholar] [CrossRef]
  27. Fan, S.; Zhang, B.; Li, J.; Huang, W.; Wang, C. Effect of spectrum measurement position variation on the robustness of NIR spectroscopy models for soluble solids content of apple. Biosyst. Eng. 2016, 143, 9–19. [Google Scholar] [CrossRef]
  28. Mehmood, T.; Liland, K.H.; Snipen, L.; Sæbø, S. A review of variable selection methods in Partial Least Squares Regression. Chemom. Intell. Lab. Sys. 2012, 118, 62–69. [Google Scholar] [CrossRef]
  29. Cheng, J.H.; Sun, D.W. Rapid and non-invasive detection of fish microbial spoilage by visible and near infrared hyperspectral imaging and multivariate analysis. LWT Food Sci. Technol. 2015, 62, 1060–1068. [Google Scholar] [CrossRef]
  30. Pierna, J.A.F.; Abbas, O.; Baeten, V.; Dardenne, P. A Backward Variable Selection method for PLS regression (BVSPLS). Anal. Chim. Acta 2009, 642, 89–93. [Google Scholar] [CrossRef] [PubMed]
  31. Belmerhnia, L.; Djermoune, E.H.; Carteret, C.; Brie, D. Simultaneous variable selection for the classification of near infrared spectra. Chemom. Intell. Lab. Sys. 2021, 211, 104268. [Google Scholar] [CrossRef]
  32. Bonifazi, G.; Capobianco, G.; Gasbarrone, R.; Serranti, S. Hazelnuts classification by hyperspectral imaging coupled with variable selection methods. In Sensing for Agriculture and Food Quality and Safety XIII; International Society for Optics and Photonics: Bellingham, WA, USA, 2021; Volume 11754, p. 117540Q. [Google Scholar]
  33. Camacho, W.; Karlsson, S. Quantification of antioxidants in polyethylene by near infrared (NIR) analysis and partial least squares (PLS) regression. Int. J. Polym. Anal. Charact. 2002, 7, 41–51. [Google Scholar] [CrossRef]
  34. González-Martín, I.; González-Pérez, C.; Hernández-Méndez, J.; Alvarez-García, N. Determination of fatty acids in the subcutaneous fat of Iberian breed swine by near infrared spectroscopy (NIRS) with a fibre-optic probe. Meat Sci. 2003, 65, 713–719. [Google Scholar] [CrossRef]
  35. Rinnan, Å.; van den Berg, F.; Engelsen, S.B. Review of the most common preprocessing techniques for near-infrared spectra. TrAC Trends Anal. Chem. 2009, 28, 1201–1222. [Google Scholar] [CrossRef]
  36. Esquerre, C.; Gowen, A.A.; Burger, J.; Downey, G.; O′Donnell, C.P. Suppressing sample morphology effects in near infrared spectral imaging using chemometric data pre-treatments. Chemom. Intell. Lab. Syst. 2012, 117, 129–137. [Google Scholar] [CrossRef]
  37. Vidal, M.; Amigo, J.M. Preprocessing of hyperspectral images. Essential steps before image analysis. Chemom. Intell. Lab. Syst. 2012, 117, 138–148. [Google Scholar] [CrossRef]
  38. Feng, Y.Z.; Sun, D.W. Near-infrared hyperspectral imaging in tandem with partial least squares regression and genetic algorithm for non-destructive determination and visualization of Pseudomonas loads in chicken fillets. Talanta 2013, 109, 74–83. [Google Scholar] [CrossRef]
  39. Calvini, R.; Ulrici, A.; Amigo, J.M. Practical comparison of sparse methods for classification of Arabica and Robusta coffee species using near infrared hyperspectral imaging. Chemom. Intell. Lab. Syst. 2015, 146, 503–511. [Google Scholar] [CrossRef] [Green Version]
  40. Barnes, R.J.; Dhanoa, M.S.; Lister, S.J. Standard normal variate transformation and de-trending of near-infrared diffuse reflectance spectra. Appl. Spectrosc. 1989, 43, 772–777. [Google Scholar] [CrossRef]
  41. Savitzky, A.; Golay, M.J.E. Smoothing and Differentiation of Data by Simplified Least Squares Procedures. Anal. Chem. 1964, 36, 1627–1639. [Google Scholar] [CrossRef]
  42. Amigo, J.M.; Babamoradi, H.; Elcoroaristizabal, S. Hyperspectral image analysis. A tutorial. Anal. Chim. Acta 2015, 896, 34–51. [Google Scholar] [CrossRef]
  43. Sun, D.W.; Rinnan, Å.; Nørgaard, L.; van den Berg, F.; Thygesen, J.; Bro, R.; Engelsen, S.B. Chapter 2–Data Preprocessing. In Infrared Spectroscopy for Food Quality, Analysis and Control; Academic Press: Cambridge, MA, USA, 2009; pp. 29–50. [Google Scholar] [CrossRef]
  44. Li, H.; Liang, Y.; Xu, Q.; Cao, D. Key wavelengths screening using competitive adaptive reweighted sampling method for multivariate calibration. Anal. Chim. Acta 2009, 648, 77–84. [Google Scholar] [CrossRef]
  45. Wang, Y.; Jiang, F.; Gupta, B.B.; Rho, S.; Liu, Q.; Hou, H.; Jing, D.; Shen, W. Variable Selection and Optimization in Rapid Detection of Soybean Straw Biomass Based on CARS. IEEE Access 2018, 6, 5290–5299. [Google Scholar] [CrossRef]
  46. Pieszczek, L.; Daszykowski, M. Improvement of recyclable plastic waste detection–A novel strategy for the construction of rigorous classifiers based on the hyperspectral images. Chemom. Intell. Lab. Sys. 2019, 187, 28–40. [Google Scholar] [CrossRef]
  47. Detection of Tumoral Epithelial Lesions Using Hyperspectral Imaging and Deep Learning. Available online: https://www.specim.fi/downloads/SisuCHEMA_2_2015.pdf (accessed on 23 August 2021).
  48. Calvini, R.; Orlandi, G.; Foca, G.; Ulrici, A. Development of a classification algorithm for efficient handling of multiple classes in sorting systems based on hyperspectral imaging. J. Spec. Imaging 2018, 7, 1–15. [Google Scholar] [CrossRef] [Green Version]
  49. Vidal, M.; Gowen, A.; Amigo, J.M. NIR hyperspectral imaging for plastics classification. NIR News 2012, 23, 13–15. [Google Scholar] [CrossRef]
  50. Serranti, S.; Gargiulo, A.; Bonifazi, G. Classification of polyolefins from building and construction waste using NIR hyperspectral imaging system. Resour. Conserv. Recycl. 2012, 61, 52–58. [Google Scholar] [CrossRef]
  51. Hu, B.; Serranti, S.; Fraunholcz, N.; Di Maio, F.; Bonifazi, G. Recycling-oriented characterization of polyolefin packaging waste. Waste Manag. 2013, 33, 574–584. [Google Scholar] [CrossRef]
  52. Luciani, V.; Bonifazi, G.; Hu, B.; Rem, P.; Serranti, S. Upgrading of PVC rich wastes by magnetic density separation and hyperspectral imaging quality control. Waste Manag. 2015, 45, 118–125. [Google Scholar] [CrossRef]
  53. Serranti, S.; Luciani, V.; Bonifazi, G.; Hu, B.; Rem, P. An innovative recycling process to obtain pure polyethylene and polypropylene from household waste. Waste Manag. 2015, 35, 12–20. [Google Scholar] [CrossRef]
  54. Bro, R.; Smilde, A.K. Principal component analysis. Anal. Methods 2014, 6, 2812–2831. [Google Scholar] [CrossRef] [Green Version]
  55. Ballabio, D.; Consonni, V. Classification tools in chemistry. Part 1: Linear models. PLSDA. Anal. Methods 2013, 5, 3790–3798. [Google Scholar] [CrossRef]
  56. Balage, J.M.; Amigo, J.M.; Antonelo, D.S.; Mazon, M.R.; Silva, S.d.L. Shear force analysis by core location in Longissimus steaks from Nellore cattle using hyperspectral images—A feasibility study. Meat Sci. 2018, 143, 30–38. [Google Scholar] [CrossRef] [PubMed]
  57. Currà, A.; Gasbarrone, R.; Cardillo, A.; Trompetto, C.; Fattapposta, F.; Pierelli, F.; Missori, P.; Bonifazi, G.; Serranti, S. Near-infrared spectroscopy as a tool for in vivo analysis of human muscles. Sci. Rep. 2019, 9, 8623. [Google Scholar] [CrossRef]
  58. Suhandy, D.; Yulia, M. Potential application of UV-visible spectroscopy and PLS-DA method to discriminate Indonesian CTC black tea according to grade levels. IOP Conf. Ser. Earth Environ. Sci. 2019, 258, 012042. [Google Scholar] [CrossRef]
  59. Barboza, E.S.; Lopez, D.R.; Amico, S.C.; Ferreira, C.A. Determination of a recyclability index for the PET glycolysis. Res. Conserv. Recyc. 2009, 53, 122–128. [Google Scholar] [CrossRef]
  60. Awaja, F.; Pavel, D. Recycling of PET. Eur. Polymer J. 2005, 41, 1453–1477. [Google Scholar] [CrossRef]
  61. Vollmer, I.; Jenks, M.J.; Roelands, M.C.; White, R.J.; van Harmelen, T.; de Wild, P.; Weckhuysen, B.M. Beyond mechanical recycling: Giving new life to plastic waste. Angew. Chem. Int. Ed. 2020, 59, 15402–15423. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Source images of calibration dataset showing PET (red square) and contaminant (green square) flakes (a) and source image of the prediction dataset showing PET and contaminant flakes (the latter marked by green circles) (b).
Figure 1. Source images of calibration dataset showing PET (red square) and contaminant (green square) flakes (a) and source image of the prediction dataset showing PET and contaminant flakes (the latter marked by green circles) (b).
Jimaging 07 00181 g001
Figure 2. Sisuchema XLTM chemical imaging workstation (Specim Ltd.).
Figure 2. Sisuchema XLTM chemical imaging workstation (Specim Ltd.).
Jimaging 07 00181 g002
Figure 3. Average raw spectra and standard deviation of PET (a) and contaminant (b) classes in the SWIR range.
Figure 3. Average raw spectra and standard deviation of PET (a) and contaminant (b) classes in the SWIR range.
Jimaging 07 00181 g003aJimaging 07 00181 g003b
Figure 4. Average pre-processed spectra of PET and contaminant classes in the SWIR range adopting three sets of preprocessing techniques: Set 1: Detrend + Smoothing + MC (a); Set 2: SNV + MC (b); and Set 3: MSC + Derivative + MC (c).
Figure 4. Average pre-processed spectra of PET and contaminant classes in the SWIR range adopting three sets of preprocessing techniques: Set 1: Detrend + Smoothing + MC (a); Set 2: SNV + MC (b); and Set 3: MSC + Derivative + MC (c).
Jimaging 07 00181 g004aJimaging 07 00181 g004b
Figure 5. PCA results for the preprocessing Set 1 (Detrend + Smoothing + MC): PCA score plot (PC1-PC2) (a) and loadings plot of PC1 and PC2 related to PET and contaminant classes (b).
Figure 5. PCA results for the preprocessing Set 1 (Detrend + Smoothing + MC): PCA score plot (PC1-PC2) (a) and loadings plot of PC1 and PC2 related to PET and contaminant classes (b).
Jimaging 07 00181 g005
Figure 6. PCA results for the preprocessing Set 2 (SNV + MC): PCA score plot (PC1-PC2) (a) and loadings plot of PC1 and PC2 related to PET and contaminant classes (b).
Figure 6. PCA results for the preprocessing Set 2 (SNV + MC): PCA score plot (PC1-PC2) (a) and loadings plot of PC1 and PC2 related to PET and contaminant classes (b).
Jimaging 07 00181 g006
Figure 7. PCA results for preprocessing Set 3 (MSC + Derivative + MC): PCA score plot (PC1-PC2) (a) and loadings plot of PC1 and PC2 related to PET and contaminant classes (b).
Figure 7. PCA results for preprocessing Set 3 (MSC + Derivative + MC): PCA score plot (PC1-PC2) (a) and loadings plot of PC1 and PC2 related to PET and contaminant classes (b).
Jimaging 07 00181 g007
Figure 8. Prediction maps as resulting from the preliminary utilization of the 3 different preprocessing strategies (i.e., Set 1: Detrend + Smoothing + MC, Set 2: SNV + MC and Set 3: MSC + Derivative + MC) applied to the reduced set of wavelengths, resulting from CARS processing, and the further PLS-DA modeling. Prediction maps related to the utilized wavelengths as resulting from preprocessing Set 1: Detrend + Smoothing + MC (a), Set 2: SNV + MC (b) and Set 3: MSC + Derivative + MC (c).
Figure 8. Prediction maps as resulting from the preliminary utilization of the 3 different preprocessing strategies (i.e., Set 1: Detrend + Smoothing + MC, Set 2: SNV + MC and Set 3: MSC + Derivative + MC) applied to the reduced set of wavelengths, resulting from CARS processing, and the further PLS-DA modeling. Prediction maps related to the utilized wavelengths as resulting from preprocessing Set 1: Detrend + Smoothing + MC (a), Set 2: SNV + MC (b) and Set 3: MSC + Derivative + MC (c).
Jimaging 07 00181 g008
Figure 9. Prediction maps for the full spectrum PLS-DA model of preprocessing Set 3 (MSC + Derivative + MC).
Figure 9. Prediction maps for the full spectrum PLS-DA model of preprocessing Set 3 (MSC + Derivative + MC).
Jimaging 07 00181 g009
Table 1. SISUCHEMA XLTM details (SWIR range).
Table 1. SISUCHEMA XLTM details (SWIR range).
Optical Characteristics
SpectrographImspector N25E
Spectral Range1000–2500 nm ±
Spectral resolution10 nm (30 µm slit)
Spectral sampling/pixel6.3 nm
Spatial resolutionRms spot radius <15 µm (320)
AberrationsInsignificant astigmatism, smile or keystone <5 µm
Numerical apertureF/2.0
Slit width options30 µm (50 or 80 µm optional)
Effective slit length9.6 mm
Total efficiency (typical)>50%, independent of polarization
Stray ligth<0.5% (halogen lamp, 1400 nm notch filter)
Field of view (mm)15 mm lens
200
Pixel dimension (mm)x0.625
y(y dimension in mm × 0.03)/9.6
Scanning speed (mm/s)72.50
Scanning rate100 hyperspectral line images/s (max), corresponding to
−60 mm/s with 600 micron pixel
Electrical Characteristics
CameraMCT camera
Pixels in full frame320 (spatial) × 256 (spectral)
Active pixels320 (spatial) × 240 (spectral)
Pixel size on sampleScalable from 30 to 300 µm
Cooling4-stage Peltier for detector array, additional Peltier for active cooling of the detector package
Camera output14-bit LVDS
Signal to noise ratio800:1 (at max signal level)
Frame grabberNational Instruments PCL-1422
Table 2. Different preprocessing sets and corresponding selected wavelengths, using CARS method.
Table 2. Different preprocessing sets and corresponding selected wavelengths, using CARS method.
SetPreprocessingSelected Wavelengths (nm)Number of
Wavelengths
1Detrend + Smoothing + MC1000, 1018, 1024, 1030, 1308, 1314, 1320, 1327, 1333, 1339, 1723, 1729, 1905, 1911, 1917, 2086, 2092, 2099, 2105, 2249, 2255, 2261, 2442, 2448 and 245425
2SNV + MC1000, 1018, 1024, 1030, 1131, 1207, 1308, 1314, 1320, 1327, 1333, 1339, 1346, 1346, 1654, 1723, 1911, 1917, 1923, 2249, 2255, 2261, 2448, 2454, 2479, 2486, 2492, 2498 and 250029
3MSC + Derivative + MC1049, 1055, 1062, 1119, 1291, 2217, 2224, 2274, 2280, 2286, 2292, 2299, 2411 and 241714
Table 3. Root mean square error for calibration (RMSEC) and cross-validation (RMSECV) for the three PLS-DA models constructed for preprocessing Set 1, 2, and 3.
Table 3. Root mean square error for calibration (RMSEC) and cross-validation (RMSECV) for the three PLS-DA models constructed for preprocessing Set 1, 2, and 3.
Preprocessing SetClassesRMSECRMSECVLVs Number
Set 1
(Detrend + Smoothing + MC)
PET0.2479650.2483364
Contaminant0.2479650.248336
Set 2
(SNV + MC)
PET0.2377050.2378033
Contaminant0.2377050.237803
Set 3
(MSC + Derivative + MC)
PET0.1264120.1266123
Contaminant0.1264120.126612
Table 4. PLS-DA classification performances in variables selection mode, for calibration (CAL), cross-validation (CV) and prediction (PRED) phases.
Table 4. PLS-DA classification performances in variables selection mode, for calibration (CAL), cross-validation (CV) and prediction (PRED) phases.
PLS-DA ModelClassesSensitivitySpecificityEfficiency (PRED)
Set 1
Detrend + Smoothing + MC
(4 LVs)
CALPET0.9740.9890.969
Contaminant0.9890.974
CVPET0.9740.988
Contaminant0.9880.974
PREDPET0.9830.957
Contaminant0.9570.983
Set 2
SNV + MC
(3 LVs)
CALPET0.9920.9990.987
Contaminant0.9990.992
CVPET0.9920.999
Contaminant0.9990.992
PREDPET0.9950.979
Contaminant0.9790.995
Set 3
MSC + Derivative + MC
(3 LVs)
CALPET0.9860.9980.991
Contaminant0.9980.986
CVPET0.9860.998
Contaminant0.9980.986
PREDPET0.9940.988
Contaminant0.9880.994
Table 5. Root mean square error for calibration (RMSEC) and cross-validation (RMSECV) for the full spectrum PLS-DA model using preprocessing Set 3 (MSC + Derivative + MC).
Table 5. Root mean square error for calibration (RMSEC) and cross-validation (RMSECV) for the full spectrum PLS-DA model using preprocessing Set 3 (MSC + Derivative + MC).
Preprocessing SetClassesRMSECRMSECVLVs Number
Set 3
(MSC + Derivative + MC)
PET0.1055490.1056955
Contaminant0.1055490.105695
Table 6. Full spectrum PLS-DA (Set 3: MSC + Derivative + MC) classification performances for calibration (CAL), cross-validation (CV), and prediction (PRED) phases.
Table 6. Full spectrum PLS-DA (Set 3: MSC + Derivative + MC) classification performances for calibration (CAL), cross-validation (CV), and prediction (PRED) phases.
PLS-DA ModelClassesSensitivitySpecificityEfficiency (PRED)
Full spectrum
PLS-DA
(Set 3: MSC + Derivative + MC)
CALPET0.9860.9981.000
Contaminant0.9980.986
CVPET0.9860.998
Contaminant0.9980.986
PREDPET1.0001.000
Contaminant1.0001.000
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cucuzza, P.; Serranti, S.; Bonifazi, G.; Capobianco, G. Effective Recycling Solutions for the Production of High-Quality PET Flakes Based on Hyperspectral Imaging and Variable Selection. J. Imaging 2021, 7, 181. https://0-doi-org.brum.beds.ac.uk/10.3390/jimaging7090181

AMA Style

Cucuzza P, Serranti S, Bonifazi G, Capobianco G. Effective Recycling Solutions for the Production of High-Quality PET Flakes Based on Hyperspectral Imaging and Variable Selection. Journal of Imaging. 2021; 7(9):181. https://0-doi-org.brum.beds.ac.uk/10.3390/jimaging7090181

Chicago/Turabian Style

Cucuzza, Paola, Silvia Serranti, Giuseppe Bonifazi, and Giuseppe Capobianco. 2021. "Effective Recycling Solutions for the Production of High-Quality PET Flakes Based on Hyperspectral Imaging and Variable Selection" Journal of Imaging 7, no. 9: 181. https://0-doi-org.brum.beds.ac.uk/10.3390/jimaging7090181

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop