Nothing Special   »   [go: up one dir, main page]

Next Article in Journal
Active-Optical Sensors Using Red NDVI Compared to Red Edge NDVI for Prediction of Corn Grain Yield in North Dakota, U.S.A.
Next Article in Special Issue
Quantitative Ethylene Measurements with MOx Chemiresistive Sensors at Different Relative Air Humidities
Previous Article in Journal
Multi-UAV Routing for Area Coverage and Remote Sensing with Minimum Time
Previous Article in Special Issue
Quantification Method for Electrolytic Sensors in Long-Term Monitoring of Ambient Air Quality
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Electronic Nose Feature Extraction Methods: A Review

College of Electronics and Information Engineering, Southwest University, Chongqing 400715, China
*
Author to whom correspondence should be addressed.
Sensors 2015, 15(11), 27804-27831; https://doi.org/10.3390/s151127804
Submission received: 30 August 2015 / Revised: 10 October 2015 / Accepted: 27 October 2015 / Published: 2 November 2015
(This article belongs to the Special Issue Gas Sensors—Designs and Applications)

Abstract

:
Many research groups in academia and industry are focusing on the performance improvement of electronic nose (E-nose) systems mainly involving three optimizations, which are sensitive material selection and sensor array optimization, enhanced feature extraction methods and pattern recognition method selection. For a specific application, the feature extraction method is a basic part of these three optimizations and a key point in E-nose system performance improvement. The aim of a feature extraction method is to extract robust information from the sensor response with less redundancy to ensure the effectiveness of the subsequent pattern recognition algorithm. Many kinds of feature extraction methods have been used in E-nose applications, such as extraction from the original response curves, curve fitting parameters, transform domains, phase space (PS) and dynamic moments (DM), parallel factor analysis (PARAFAC), energy vector (EV), power density spectrum (PSD), window time slicing (WTS) and moving window time slicing (MWTS), moving window function capture (MWFC), etc. The object of this review is to provide a summary of the various feature extraction methods used in E-noses in recent years, as well as to give some suggestions and new inspiration to propose more effective feature extraction methods for the development of E-nose technology.

1. Introduction

An electronic nose (E-nose) is an instrument which comprises an array of electronic chemical sensors with partial specificity and an appropriate pattern-recognition system, capable of recognising simple or complex odors. It is specifically used to sense odorant molecules in analogy to the human nose. However, the architecture of an E-nose also applies in gas sensing for the detection of individual components or mixtures of gases/vapours [1], which is playing an increasing role in general purpose detection of gases in many applications such as odor analysis [2,3,4], quality control of food industry [5,6,7,8,9,10], environment protection [11,12,13], public health [14,15,16,17,18,19], explosives detection [20] and spaceflight applications [21]. The main hardware component of an E-nose is an array of non-specific gas sensors, i.e., sensors that interact with a broad range of chemicals with varying strengths. Correspondingly, the analyte stimulates sensors in the array, which elicits a characteristic response called “fingerprint”. The main software component of an E-nose is its feature extraction and pattern recognition algorithms, which process the sensor characteristic response, extract and select useful information and realize the pattern recognition.
According to the components of an E-nose system, the performance improvement of an E-nose system mainly contains three kinds of optimization: sensitive material selection and sensor array optimization, feature extraction and selection method, and pattern recognition method. Many research groups are focusing on improving the performance of E-noses [22,23,24]. In fact, they mainly emphasize the three types of optimizations. Sensitive material selection and sensor array optimization refers to the hardware structure of the E-nose. If the E-nose is homemade, a more typical approach is to first research the chemical composition of the samples and then select sensors which are responsive to those chemical groups with cross sensitivity to form the sensor array. Then from this, you would tune and adjust the sensors in the array according to the responses, which is called sensitive material selection and sensor array optimization. The optimization is all-round, as the type of sensors, the number of sensors, even the position, order and topological structure of sensors etc., can affect the kinetic processes of the sensor response. However, if a commercial E-nose is used, we just get what is in the array and cannot change the hardware structure of the E-nose. This means that we cannot perform sensor array optimization from the hardware perspective. If we wish to improve the performance of an E-nose we must emphasize the feature extraction and selection method as well as the pattern recognition method. Although performance improvement of an homemade E-nose system can be done through sensitive material selection and sensor array optimization, the feature extraction method is one of the key points of performance improvement of E-nose systems because feature extraction is the first step of the sensor signal processing and the feature selection [25] and pattern recognition are performed on the basis of features which play a significant role in affecting the effect of the subsequent feature selection and pattern recognition algorithms. No matter which pattern recognition method is been selected its recognition rate has a strong relationship with the sample distribution in the feature space. The aim of feature extraction is to extract robust information from the characteristic sensor response with less redundancy, which can represent the different “fingerprint” patterns well, to ensure the effectiveness of the subsequent pattern recognition algorithm. However, it is still a great challenge for people to extract features from the “fingerprint” to further improve the pattern recognition accuracy.
According to the research on feature extraction during the most recent twenty years, there are many feature extraction methods which have been used in E-nose applications. Most of these feature extraction methods can be roughly divided into three types: the first type is to extract piecemeal signal features [26] from the original response curves of sensors, such as maximum values, integrals, differences, primary derivatives, secondary derivatives, the adsorption slope, and the maximum adsorption slope at a specific interval from the response curves, etc. The second type is based on curve fitting [27], which fits the response curves based on a specific model and extracts a set of fitting parameters as the features. The third one is based on some transforms, and very often the fast Fourier transform (FFT) [28,29,30,31] and the discrete wavelet transform (DWT) have been used as feature extraction tools for this purpose [6,28,32].
Besides the above conventional feature extraction methods, many new original methods have been proposed in recent years. Phase space (PS) [33] and dynamic moments (DM) [34], which are usually used in dynamical systems, are also applied to E-nose signals for extraction. Parallel factor analysis (PARAFAC) [35,36] is one of the most popular multi-way data decomposition methods and is attracting increasing interest because it is a processing technique that simultaneously determines the pure contributions to the dataset, optimizing each factor as a time, in trilinear systems. Energy vector (EV) [37] is a vector of energy, which contains the energy of each sensor and all the mutual energies. It is suitable to study the relations between signals of sensors belonging to the same array. In addition, the power density spectrum (PSD) method [38] is also applied in E-nose feature extraction. Another type of methods are based of window functions, especially windowed time slicing (WTS) [39,40,41], wherein the time response of each sensor is multiplied by time windows to obtain the area values and these values are further used as features.
This review provides a summary of the main methods of feature extraction used in E-noses in recently years, which are conducive to analysis and research on E-nose technology, by describing and comparing the basic types of feature extraction methods which differ as the application and E-nose experiments change and by providing examples of research in which E-noses have been utilized to analyze and detect different materials. In addition, we also hope that more potential and effective feature extraction methods can be proposed in further work to enhance the performance of E-noses with the inspirations proposed in this paper.

2. Feature Extraction Methods

2.1. Feature Extraction from Original Response Curves

The first feature extraction method is to extract piecemeal signal features from the original response curves of sensors, including steady-state response and transient responses such as maximum value, integrals, derivatives, area values, rising time, falling time, rising slopes, falling slopes, etc. Because the maximum value represents the final steady-state feature of the entire dynamic response process in the final balance, which reflects the maximum reaction degree change of sensors responding to odors, it is usually used as the most common and simple E-nose feature. In the early stages of the published history of E-noses, transient information was not used in sensor signal feature extraction. A variety of steady-state models (maximum value) have been used as feature extraction method for gas sensor signals, as illustrated in Table 1.
Table 1. Some maximum value feature models.
Table 1. Some maximum value feature models.
ModelDescriptionReferences
Difference x i j = V i j max V i j min [42]
Relative difference x i j = V i j max / V i j min [43]
Fractional difference x i j = ( V i j max V i j min ) / V i j min [44]
Logarithm difference log ( V i j max / V i j min ) [45]
Sensor normalization x i j = x i j / ( x i max x i min ) [45,46]
Array normalization x i j = x i j / i x i j 2 [47,48]
x i j and V i j max represent the processed maximum value feature and the original maximum response value of the of the i-th sensor to the j-th gas, respectively. The difference method can usually eliminate the additive errors, which are added both to the baseline and the steady-state response (meaning the response of the i-th sensor to the j-th gas). There is some evidence that the relative and fractional difference are helpful to compensate for the temperature influence on the sensors and the fractional difference linearizes the mechanism that generates a concentration dependence in metal oxide chemiresistor sensors [47]. The log difference is more suitable when the variation of the concentration of the odor is very large because it is able to linearize the highly nonlinear relationship between odor concentration and the sensor output. The normalization models, which limit each sensor output between 0 and 1, thereby keep each element of the response vector at the same magnitude. Not only can it reduce the calculation error of stoichiometric recognition, but also is very effective when it is not the odor concentration which is of interest, but rather the precise identification of the odor. This method for example provides good recognition results for the discrimination of several types of coffee using neural networks [48].
Besides the steady-state features of response curves, various other transient features such as derivatives and integrals of original response curves were taken into research in many special applications. Balasubramanian et al. used a MOS-based E-nose system, which was designed and fabricated at the Bio-Imaging and Sensing Center, North Dakota State University (Fargo, ND, USA) to analyze the headspace from beef strip loins stored at 4 °C and 10 °C. The designed system consisted of a cylindrical sampling chamber made of Teflon on to which seven thick film MOS (Figaro USA Inc., Glenview, IL, USA) and one integrated sensor sensitive to temperature and relative humidity (Sensirion AG, Zurich, Switzerland) were mounted at regular intervals. The headspace gas from the meat packages was drawn into the chamber using a diaphragm-type air pump, and a fan and valve assembly fixed at the bottom of the sample chamber facilitated purging the sensor array after analysis. Six areas features under the response curves were extracted from six gas sensors (TGS 812, TGS 822, TGS 880, TGS 2602, TGS 2611-1 and TGS 2611-2). Another four features were extracted from the remaining two sensors related to the relative humidity (RH), temperature and carbon dioxide readings. The results showed that the MOS-based electronic system integrated with radial basis function neural networks (RBFNN) obtained above 90% total maximum classification accuracy in identifying spoiled meat samples from the unspoiled meat samples with six area features and four additional environmental features for the meat samples stored at 4 °C and 10 °C [49].
Distante et al., used a homemade E-nose consisting of five SnO2 sensors prepared by means of sol-gel technology with Pd, Pt, Os, and Ni as doping elements. The film thickness was 100 nm and the films were deposited on alumina substrates supplied with interdigitated electrodes and platinum heater, by the spin coating technique at 3000 rpm, dried at 80 °C and heat treated in air at 600 °C. After deposition, the sensors were mounted onto a TO8 socket and inserted in a test chamber. Three volatile organic compounds (VOCs) were investigated: acetone, hexanal and pentanone in 50% relative humidity (RH) and dry air. Several feature extraction methods, both steady-state feature maximum value and transient feature, were considered. The results showed that integral and derivatives methods have higher recognition rates than that of maximum values and they also concluded that the desorption stage is more informative than the adsorption stage of the original response curves [50].
Roussel et al., applied a homemade E-nose to discriminate satisfactory wines and unsatisfactory vinegar flavour wines. The measurement device is a metal oxide gas sensor array containing five Figaro TGS series SnO2 sensors (TGS 800, 822, 824, 825 and 842) and modified to host a thermo-hygrometer and to acquire the sensor voltage on a 12 bit A/D converter. Twenty-nine features are extracted from each sensor response curve including seven mean response values on different time intervals, six slopes on different time intervals, eight primary derivatives and eight secondary derivatives. And then these features have been validated and evaluated qualitatively by three specific indexes: the repeatability, discrimination and redundancy. However, the authors just proposed the three indexes to qualitatively analyze the features but did not give any classification results for specific data using different features, which were considered better than others according to their proposed three evaluating indexes [51].
Llobet et al. established an array to discriminate between and quantify ethanol, toluene and o-xylene at concentration ranges from 25 ppm to 100 ppm, which consisted of four commercially available Taguchi gas sensors: two TGS822, one TGS813 and one TGS800. The model TGS822 exhibits high sensitivity to alcohol and organic solvents. The model TGS813 is mainly sensitive to general combustible gases and the model TGS800 is sensitive to gasoline exhaust and other air contaminants. The steady-state feature (conductance change Δ G ) and conductance rise time (Tr, measured from 20% to 60% of Δ G ) were extracted as features and they found that the using Tr as feature gave much better qualitative recognition results and using Δ G and Tr as features results in equal quantitative recognition rates, as shown in Table 2, which meant that Tr was concentration-independent [26].
Table 2. Classification rates (CR) of different features in [26].
Table 2. Classification rates (CR) of different features in [26].
FeaturesQualitative CR (%)Quantitative CR (%)
Difference maximum value Δ G 66.795
Conductance rise time Tr9295
Eklöv et al., investigated the concentration of gas mixtures of hydrogen and ethanol with a homemade E-nose with four platinum MOSFET sensors, which had a thin discontinuous platinum film covered gate oxide. The gate voltages of each MOSFET sensor were detected and the sampling frequency was set at 1 Hz. All measurements consisted of an absorption stage for 2 min and a desorption stage for 8 min with a flow rate of 100 mL/min. A sequence of different concentrations of hydrogen and ethanol that ranged between 0 and 50 ppm were tested. The extracted features could be divided into two types: simple parameters from the original response curves and coefficients from different curve fitting models. It is necessary to select these parameters systematically and take into consideration the noise properties as well as the relationship between the features and the recognition ability of the features in a specific problem [52]. The features extracted from original response curves are shown in Table 3. Several combinations of features were proposed and gave satisfactory concentration prediction results.
Table 3. Description of features extracted from original response curves in [52].
Table 3. Description of features extracted from original response curves in [52].
ParametersDescription
Baseline 1 5 T = gason 4 s gasOn (sensor value)
Final response, responseSensor value (averaged over 5 s) at gasOff baseline
30/90 s on/off responseSensor value (averaged over 5 s) 30/90 s after gasOn/Off—baseline
Maximum responseMax(sensor value)—baseline
Min/max derivativeMin/max difference between two samples during measurement
On/off derivative(Sensor value 10 s after gasOn/Off−baseline)/10
Plateau derivative(Response 90 s on response)/30
Derivative T gason gasOff (sensor value baseline)
Off integral T gasoff gasOff + 119 s (sensor value baseline)
Short on/off integral T = gason/gasoff gasOn / gasoff + 9 s (sensor value baseline)
Response/on integralResponse/on integral
T0-90%Time from gasOn for sensor value to reach baseline + 0.9 × response
T0-60%Time from gasOn for sensor value to reach baseline + 0.6 × response
T100-10%Time from gasOff for sensor value to reach baseline + 0.1 × response
T100-40%Time from gasOff for sensor value to reach baseline + 0.4 × response
These parameters can be assessed according to the signal to noise ratio (SNR), correlation and significance in a 2-dimensional PCA plot. The signal to standard deviation ratio ( S / σ - r a t i o n ) [52], defined in Equation (1), was used to estimate the property of the features. It can provide similar information as a SNR:
S / σ - r a t i o n = stdev all concentration steps ( average all measurments at each concentration step ( f e a t u r e ) ) stdev all concentration steps ( f e a t u r e average all measurments at this concentration step ( f e a t u r e ) )
In addition, principal component analysis (PCA) was used to analyze the features. The PCA result gave score plots and loading plots. The trends, groups and outliers of the samples can be observed in score plots and the correlation and similarity between the features can be observed. A longer distance from the origin in the loading plot and a larger SNR implied that the feature contained more information. The artificial neural networks (ANN) classification results showed that: (1) a large error is obtained only with the final response features; (2) it was possible to improve the prediction accuracy to use features according to the S / σ - r a t i o n calculation; (3) there was information which can improve the recognition ability in the loading plot; (4) features extracted from original response curves may gave comparable results as the complex curve fitting coefficients; (5) there was more information about the detected materials in the shape of the whole response than only the response value.
Zhang et al. used an array of six TGS gas sensors to detect four common flammable liquids (gasoline, kerosene, diesel oil and ethanol) and three incombustible drinks (ice black tea, orange juice and Coca Cola) in less than 10 s for each recognition. Because they needed to realize fast recognition with an E-nose, they must extract robust information from the first few seconds of the sensor response curves as less redundant as possible. They researched four features (integrals, differences, primary derivatives and secondary derivatives) at a specific interval from the original response curves to classification [53]. Because integrals, differences, primary derivatives and secondary derivatives of the response curves are continuous with the interval changes, different positions of the curves can contain different information, which may consecutively change with the interval changes. Therefore, a calculating process was applied to fix at which interval of the response curves robust information can be extracted St shown in Equation (2), was use to evaluated the information quality of the feature extracted at the time t of the response curves from the sensor array, which is a criterion to measure the separative degree among different classes of data based on the scatter matrix [54]. The expression of St is:
S t = T r { S w t 1 S m t }
where Swt is the scatter matrix within a class, while S m t is the scatter matrix among different classes. PCA and discriminant function analysis (DFA) results showed that 85.7% classification rate was obtained with the four features in optimal time t, while only 57.1% classification rate was obtained with the maximum value feature.
Šetkus et al., applied a homemade E-nose with s metal oxide (MOX) sensors array that contained commercial- and laboratory-made sensors to investigate some volatile organic compounds (VOCs), namely acetone, acetic acid, butyric acid, and acetaldehyde. The sensor array included two TGS type sensors (TGS2602, TGS2610). The laboratory-made sensors were based on tin oxide thin films, which were modified by catalytic metals such as Au, Pt, Pd and Ru, and were grown by rheotaxial (fused metal) growth and the thermal oxidation technique (RGTO). The authors researched three kinetic features: (1) the time the maximum value (τmax); (2) the time interval between the beginning of the adsorption and the end of the transient signal; (3) the slope of the response calculated in the desorption stage. PCA results showed that the dynamic features added significant information and allowed a better discrimination of the VOCs than that of the relative maximum response value [55].
Table 4. Description of features extracted from original response curves in [56].
Table 4. Description of features extracted from original response curves in [56].
Extracted FeatureRepresent Meaning
Max slopeThe respond rate of sensor to different vinegar gas
MaximumThe maximum respond value
Average of last 20 pointsThe stationary phase of equilibrium between reversible adsorption and desorption
Average of whole pointsSensor respond value during the whole process
Zou et al., used an E-nose based on five tin-oxide sensors from Figaro Co. Ltd. (TGS 813, TGS 880, TGS 822, TGS 825, TGS 812), one humidity sensor (HS-01) and one temperature sensor (Pt100) to identify two different types of vinegars. Four features are extracted from each curve, which are shown in Table 4. They are the slope max, maximum, average of the last 20 points and the average of whole points of curve [56]. In addition, an evaluation index, called “distinguish index” (DI), was proposed to evaluate the features.
Wei et al., extracted the maximum values, area values and 70th s values as the features from the E-nose responses and the PCA results indicated that the 70th s values provided the most accurate results when classifying unshelled peanuts and peanut kernels, which could be distinguished reliably by using a two-dimensional PCA [57].
From the above, many features can be extracted from original response curves (Table 5). Some of the features have specific physical meanings and reflect different information of the reaction kinetics at different aspects. Maximum value represents the final steady-state feature of the entire dynamic response process in the final balance, reflecting the maximum reaction degree change of sensors responding to odors. Intervals may represent the accumulative total of the reaction degree changing and differences may represent the degree of the reaction. Primary derivatives (slopes) may represent the rate of the reaction of sensors responding to odors and secondary derivatives may represent the acceleration of the reaction, etc. [53,58].
Table 5. Summary of commonly used features extracted from original response curves
Table 5. Summary of commonly used features extracted from original response curves
FeatureDescriptionReferences
Maximum responseMax(sensor value)[42,43,44,45,46,47,51,52,55,56,59]
Responses of special timeResponse value of special time in the whole response curve[51,52,56,57]
Time of special responsesTime of special response value in the whole response curve[26,51,52,55]
AreaArea values of sensor response curve and time axis surrounded[49,57]
Integral I = a b x ( t ) d t [50,52,53,59]
Derivative D = d x ( t ) d t [50,51,52,53,56]
Difference x ( t j ) x ( t i ) [52,53,55,59]
Second derivative D = d 2 x ( t ) d 2 t [51,53]
By extracting these features we can get information about different aspects from the original reactions between the sensors and the odors. However, it is evident that which extraction methods from the original response curves are appropriate depends on the different sensors and applications and for different applications and sensors, the feature extraction methods which give the best features are different. Moreover, researchers are likely to use composite features because these types of features which usually contain steady-state and transient features can obtain better performance. It is not easy to draw a conclusion about which method is the best one. On the whole, for maximum features, the normalization models give better performance than the other models. The transient features can receive better recognition accuracy than that of the steady-state feature and especially integral methods generally give the best features for many different E-nose applications. Based on the previous researches, a summary is listed in Table 5, which provides a brief description of the features commonly used in previous researches.

2.2. Features Extraction from Curve Fitting Parameters

Instead of extracting characteristics from the raw data curve, an alternative method is to model sensor response curves and then extract features based on the model coefficients. Curve fitting is a data processing method which approximately characterizes or analogizes the functional relationship of discrete points using continuous curves. The curve fitting method approximates discrete data using analytical expressions. The basic idea is to determine the fitting function in accordance with the trend of the discrete points. It utilizes the discrete points, but is not limited to these discrete points. In curve fitting, the critical and complex problem is how to select and envision the specific form of the unknown function.
There are two strict mathematical laws on theories to follow for curve fitting. However, two approaches are generally followed: (1) determine the basic type of function through the study of the physical concepts between the variables and deep understanding of professional knowledge and (2) determine the type of function through the observation of general trends of the curve of experimental data. The root mean square error (RMSE) for the validation data can be calculated for each curve fitting model and used to evaluate the effect of the curve fitting.
Polynomial function models are usually applied to fit a number of kinds of curves. The curve fitting and features extraction can be carried out on not only the adsorption or desorption parts of the response curves respectively, but also the entire response process. In polynomial functions model, n-degree polynomials are used to fit the response curves:
y = A 0 + A 1 x + A 2 x 2 + A 3 x 3 + + A n x n
During the adsorption part the response curves, y represents the sensor response value, x represents the time from gas on to gas off. During the desorption part of the response curves, y represents the final response value subtracting the sensor value, x represents the time from gas off to the end. The fitting model parameters A0, A1, A2, A3, ···, An are used as features.
Exponential function models which are expressed in Equation (4) are another common model in studies used to model E-nose response curves [59]. When the number N equals 1 or 2, the exponential function model is called single-exponential function model and double-exponential functions, which are the most commonly used exponential function models:
y = A 0 + i N A i exp ( x T i ) , i = 1 , 2 , 3 ,
where y and x are defined as the same in the polynomial function model and the fitting model parameters A0, T1, ···, TN are used as features.
Auto-Regressive with eXogenous variables (ARX) models are widely built to represent dynamic behavior, used for time series prediction and system identification, and used a different way from the former two types of methods. The combination of sensor and target gas can be regarded as a dynamic system, which gives a response when exposed to a step function. This type of system can be modeled using an ARX model. In [52], not only are simple parameters from original response curves extracted as the features shown in Table 3, but also coefficients from different curve fitting models are taken into consideration. Four models, third-order polynomial function, single-exponential functions, double-exponential functions and ARX model, shown in Table 6, were applied to fit the response curves of E-nose and the corresponding fitting model parameters were extracted as features for recognition.
Table 6. Description of features extracted from curve fitting models in [52].
Table 6. Description of features extracted from curve fitting models in [52].
ParametersDescription
Polynomial on/off y = A 0 + A 1 x + A 2 x 2 + A 3 x 3
On : y =(senor value baseline) and  x =time from gasOn to gasOff
Off : y =(maximum value senor value) and  x =time from gasOff to end
Exponential on/off y = A ( 1 exp ( x T ) ) , where y and x are defined as in the polynomial fit.
Exponential on/off y = A 0 + A 1 exp ( x T 1 ) + A 2 exp ( x T 2 ) , where y and x are defined as in the polynomial fit.
ARX on/off y ( t ) = a 1 y ( t 1 ) + a 2 y ( t 2 ) + b u ( t 1 ) ,
On: y(t) = (sensor value − baseline), t time from gas on—5 s to gas off
u(t) = 0 if test gas off and 1 if test gas on
Off: y(t) = (maximum value − sensor value), t = time from gas off—5 s to end
u(t) = 1 if test gas off and 0 if test gas on
Table 7 showed neural network prediction results for different features in [52]. As expected, only using one feature, the final response gave a large error. It is possible to improve the prediction accuracies according to the information obtained from the S / σ - r a t i o n . On/off derivatives have higher S/a-ratios than max/min derivatives and short on/off integrals, and the prediction errors showed the same tendency than using final responses, so combined with on/off derivatives as features they can obtain lower errors than for the other two combinations (models 2–4 in Table 7). Information from the PCA loading plot can be used to improve the prediction ability too. Final response, 30 s and 90 s on response, which all are adjacent in the loading plot gave worse prediction ability than that of final response together with 30 s and 90 s off response (models 5 and 6 in Table 7), which are located farther. On the whole, the curve fitting model parameters for the entire curves gave worse results compared with the features extracted from original response curves, except for the polynomial coefficients, which can obtain equally good results. Hence, curve fitting parameter features, which refer to the whole response, cannot guarantee better recognition results than those of the simple and piecemeal features combination, which gives comparable results as the more complex curve fitting parameters if chosen carefully.
Carmel et al. researched a large dataset composed of 30 volatile odorous substances using a MOSES II E-nose equipped with two sensor modules: an eight-sensor quartz-microbalance (QMB) module and, and an eight-sensor metal-oxide (MOX) module. The samples were put in 20 mL vials in a headspace sampler and the headspace contents were injected into the E-nose, first introduced into the QMB chamber and then the MOX chamber. The injection lasts for 30 s and then a 15 min desorption stage follows. In total, 300 measurements were collected with an average of 10 per chemical. The exponential model, Lorentzian model and double-sigmoid model were applied to fir a curve the response for feature extraction [27].
Table 7. Results of neural network estimation from different features in [52].
Table 7. Results of neural network estimation from different features in [52].
ModelFeaturesANN ArchitectureRMSE for Validation (ppm)
1Final response1-2-213.4
2Final response3-3-22.2
On/off derivative
3Final response3-3-22.6
Min/max derivative
4Final response3-3-23.6
Short on/off integral
5Final response3-3-24.5
30 s on response
90 s on response
6Final response3-3-22.4
30 s off response
90 s off response
7Final response5-4-21.6
30 s on/off response
90 s on/off response
8Final response5-4-21.7
On/off derivative
on/off integral
9Final response6-5-21.7
On/off derivative
Plateau derivative
Response/on integral
10Polynomial on/off8-6-21.6
111. Exponential on/off4-4-22.0
122. Exponential on/off10-6-22.1
13ARX on/off6-5-22.3
The ANN architecture column presents the number of inputs, number of neurons in the hidden layer and number of outputs.
The Lorentzian model is derived from a simple physical description of the measurement process. It applies four parameters, with respective specific physical meaning, which can be obtained from a curve fitting process. The double-sigmoid model was also used for curve fitting. The model with nine parameters function had been constructed by multiplying a monotonically decreasing sigmoid function by a monotonically increasing one. The explicit form of Lorentzian model and double-sigmoid model are expressed in Equations (5) and (6), respectively, and the details of the parameters in the two equations are presented in [27]:
x i ( t ) = { 0 ,                                                              t < t i β i τ i tan 1 ( t t i τ i ) ,                                      t i < t < t i + T β i τ i [ tan 1 ( t t i τ i ) tan 1 ( t t i T τ i ) ] ,       t > t i + T
where t is the response time and t [ 0 , T ] , T is the whole response time, ti is the time how long a particle of the investigated chemical makes its way between the inlet and sensor i, which is just the time when the signal starts to rise. The parameters ti, T, β i  and  τ i can be used as features.
x i ( t ) = α i π [ 1 exp ( ( t β i γ i + e i ) δ i ) ] η i × [ π 2 tan 1 ( t μ i v i ) ] λ i
where the parameters α i , β i γ i e i η i , μ i , v i  and  λ i can be used as features.
Table 8. Classification rates (CR) of different feature sets for QMB or MOX modules in [27].
Table 8. Classification rates (CR) of different feature sets for QMB or MOX modules in [27].
ClassifiersModulesFeature SetsCR (%)
k-Nearest Neighbor (with k = 3)QMB m a j ( ψ max , β Lor , τ Lor , t Lor , T Lor ) 97.8
( β Lor , τ Lor , T Lor ) 90
( ψ max ) 84.4
m a j ( ψ max , β Exp , τ Exp , t Exp , T Exp ) 96.7
Mahalanobis distance discrimination in four-dimensional PCA spaceQMB m a j ( ψ max , β Lor , τ Lor , t Lor , T Lor ) 83.3
( β Lor , τ Lor , T Lor ) 96.7
( ψ max ) 42.2
m a j ( ψ max , β Exp , τ Exp , t Exp , T Exp ) 60
k-Nearest Neighbor (with k = 5)MOX m a j ( ψ max , β Lor , τ Lor , t Lor , T Lor ) 93.1
( β Lor , τ Lor , T Lor ) 90
( ψ max ) 86.9
m a j ( ψ max , β Exp , τ Exp , t Exp , T Exp ) 93.5
maj(feature 1,feature 2,...,feature n) means the following: use each of the n features separately for the purpose of classification, and then decide on the final classification by a majority rule.
Fifty-one feature sets were tested and applied to 101 datasets with nine classifiers. The nine classifiers were based on two types of classifiers. A k-nearest-neighbors (k-NN) classifier [60] with k = 1 , 3 , 5 , 7 , and classifier in the base of the shortest Mahalanobis distance discrimination [60] in the 1 to 5 dimensional principal components spaces. Average classification performance was calculated and ranked for every feature sets and Table 8 showed the classification performance of several feature sets in [27]. From the ranking of the fifty-one feature sets, the best feature set was maj ( ψ max , β Lor , τ Lor , t Lor , T Lor ) , which meant classify by majority rule among ψ max , which was the difference between the maximum value and its baseline, and the parameters of the Lorentzian model. The best non-majority feature set was ( β Lor , τ Lor , T Lor ) . On the whole, the Lorentzian model was better than the exponential model with respect to classification rate and they were both much better than the standard features, which classified by majority rule among ψ max and the area beneath the curve, the area beneath the curve left of the peak, and the time it took for the signal to reach its maximum value. Among the three models investigated in [27], the analysis results showed that among the three models (exponential function models, Lorentzian model and double-sigmoid model) with respect to goodness of fit, computation speed, robustness and classification rate, the Lorentzian model was the best model for feature extraction and the most potent set of features resulting in higher classification rates than other combinations of features was m a j ( ψ max , β Lor , τ Lor , t Lor , T Lor ) . Some other models are also used to characterize E-nose measurements, such as models based on sigmoid function [61], fractional function, arc tangent function, and hyperbolic tangent function [62]. A summary of different curve fitting models is shown in Table 9.
Table 9. Summary of common used curve fitting models.
Table 9. Summary of common used curve fitting models.
ModelDescriptionReferences
Third-order polynomial function Y = A 0 + A 1 x + A 2 x 2 + A 3 x 3 [52,62]
Single-exponential function Y = A ( 1 e x p ( x T ) ) [52,62]
Double-exponential function Y = A 0 + A 1 exp ( x T 1 ) + A 2 exp ( x T 2 ) [52,62]
ARX model y ( t ) = a 1 y ( t 1 ) + a 2 y ( t 2 ) + b x ( t 1 ) [52]
Lorentzian modelEquation (5)[27]
Double-sigmoid modelEquation (6)[27]
Sigmoid function f ( k ; θ ) = θ 1 1 1 + e θ 2 ( θ 3 k ) ( 1 1 1 + e θ 4 ( θ 5 k ) ) [61]
Fractional function Y = x / A x + B [62]
Arctangent function Y = A × arctan ( x / B ) [62]
Hyperbolic tangent function Y = A × tanh ( x / B ) [62]

2.3. Features Extraction from Transform Domain

Another method to extract features for E-noses is based on some transforms and then the transform coefficients are used as features to discriminate materials. The most commonly used transforms in E-nose signal processing are the Fourier transform and wavelet transform.
The widely used Fourier transform, for which the basis functions are sine and cosine, maps the original data into a new space. It decomposes the original response into the superposition of the dc component and different harmonic components, and the feature characterized by amplitude of each component can be used for qualitative and quantitative analysis. However, the problem of using the Fourier transform is the impossibility of localizing frequencies in the time domain. This problem can be overcome if we can add the information about the localization of these frequencies in the time domain that can be addressed with the theory of wavelet transform. Wavelet transform [63] is an extension of the Fourier transform. It maps the signals into a new space with basis functions quite localizable in time and frequency space. The wavelet transform decomposes the original response into the approximation (low frequencies) and details (high frequencies). It displays good anti-interference ability [64] for the subsequent pattern recognition if one uses the wavelet coefficients of certain sub-bands as features. Unlike the Fourier transform, the shapes of components of the decomposed signal are different according to the different shapes of the using mother wavelet. The common types of wavelet are Haar wavelet, Daubechies wavelet, Symlet wavelet, Coiflets wavelet, Biorthogonal wavelet, Morlet wavelet, Gaussian wavelet, Mexican hat wavelet, Meyer wavelet, Morlet wavelet, etc.
Jia et al. [62] constructed a gas sensor array with six metal oxide sensors and one electrochemical sensor and used it to detect the seven species of pathogen most common in wound infections: P. aeruginosa, E. coli, Acinetobacter spp., S. aureus, S. epidermidis, K. pneumoniae, and S. pyogenes. They researched and compared the effect of different features of sensors including maximum values, integrals, derivatives, six types of curve fitting model parameters, FFT coefficients and DWT coefficients. In this research, a Daubechies family wavelet (db N) was used and the variant N was the order of the wavelet. For the Daubechies family wavelet, the N was the vanishing moment. A larger N can result in stronger localization ability in the frequency domain and better frequency band allocation results, and then, the energy is more concentrated. The results showed that the selected FFT feature, the amplitudes of the DC component, and DWT feature, the mean of the approximating coefficients after decomposition of 12 levels with the sixth order Daubechies family wavelet (db6), can both achieve a 100% correct classification rate, shown in Table 10, which is generally higher than that of the other features except for integral features, single-exponential function and hyperbolic tangent function fitting parameters, which can provide the same 100% correct classification rate.
Table 10. Classification rates (CR) of three types of features in [62].
Table 10. Classification rates (CR) of three types of features in [62].
Feature TypeFeatureCR (%)
Original response curveMaximum 94.29
Integrals 100
Derivatives 97.14
Curve fitting parametersThree-order polynomial function 91.43
Fractional function 97.14
Single-exponential function 100
Double-exponential function 68.57
Arctangent function 91.43
Hyperbolic tangent function 100
Transform domainFFT100
DWT100
Distante et al. [50] also analyzed and researched the effect of features of E-noses from original response curves and transform domains for the application to volatile organic compound detection. The results are shown in Table 11. It is evident that the wavelet transform coefficients can obtain the highest classification rate and it is interesting to note that integrals method produces very informative features as compared with the results obtained with the wavelet descriptors.
Table 11. Classification rates (CR) of different features in [50].
Table 11. Classification rates (CR) of different features in [50].
Feature TypeFeatureCR (%)
Original response curve Difference maximum83
Relative maximum81
Fractional maximum82
Log maximum81
Derivatives95.45
Integrals99.5
Transform domainFourier coefficient 96
Wavelet transform coefficient100
Ratton et al. researched not only the Fourier transform and wavelet transform but also the Gram Schmidt orthogonalization approach in E-nose feature extraction for the investigation of four test gases (methanol, ethanol, formaldehyde, and acetone). For this study, the wavelet transform with Haar wavelet showed the highest overall performance [65]. In [28], a 1-D discrete wavelet transform and a fuzzy adaptive resonant theory map (ARTMAP) neural network are employed as a novel feature extraction and pattern classification method. The result shows that the wavelet technique is more effective than FFT in terms of data compression and is highly tolerant of the presence of additive noise and drift in the sensor responses. Tian et al. [32] proposed a new method of detecting wound pathogens by selecting the wavelet transform coefficients preferentially with a scatter matrix and using the mean of the selected coefficients as the feature. The new feature extraction method showed high performance in identifying seven species of pathogens. Without the influence of drift, the classification rate of the new features using a probabilistic neural network classifier can reached 100% while the maximum value feature obtained only 88.57%. Under the effect of a strong drift, although both the classification rates decreased, the anti-drift ability of the new features was evidently stronger than that of the maximum value feature. Ionescu et al. [66] demonstrated that a single, thermally modulated tungsten oxide-based resistive sensor can discriminate between different vapours based on DWT. It can be found that DWT outperformed FFT in the extraction of important feature from the sensor response and allowed for straightforward gas recognition in feature space according to the 2-dimension PCA score plots.

2.4. Feature Extraction from Phase Space

The phase space (PS), which can be used for signal processing in many fields [67,68,69], is a key concept in the field of dynamic systems research. We assume that the state of a system is completely described by m scalar variables. The different states are represented as different points in an m-dimension vector space defined by an orthonormal basis where each direction is corresponding to one of the scalar variables. The basic characteristic of the PS is the correspondence between each point and the transient state of the system. A general PS can be defined according to Taken’s embedding theorem [70]. In PS, the time course of any system is described by time parametric trajectories, which contains the dynamic characteristics of the system. The trajectories present a large variety of shapes depending on the nature of the phenomena. Neglecting the scale effects, the shapes of the trajectories are expected to be associated to the properties of the phenomena. From this point of view, it is interesting to define some morphological descriptors which are able to encode the shape of the trajectories. These morphological descriptors can then be used to obtain information about the system dynamics and used as the features [71].
The first attempts to introduce the PS and dynamic moments (DM) [72,73] to represent the temporal evolution of chemical sensor signals (QCM) was presented in Refs. [33,34]. Martinelli et al. described the sensor signal changes in a PS by the response values and its first derivative like Figure 1 and then extracted the area, called phase space integral (PSI), spanned by the trajectory during the evolution either in adsorption or in desorption phase as feature, which gave a substantial improvement in terms of both error of estimation and classification [33].
Figure 1. Phase space of sensor with f ( t ) and f ( t ) as variables.
Figure 1. Phase space of sensor with f ( t ) and f ( t ) as variables.
Sensors 15 27804 g001
Vergara et al. took some morphological descriptors DM in PS analogous to the second moments of the area of a geometrical figure in a 2D space into consideration. The usefulness of the method of PS and DM is assessed by analyzing the transient response of metal oxide gas sensors either to a step change in gas concentration or to thermal modulation [71].
In [74], a gas sensor array with six TGS sensors, which were TGS813, TGS2600, TGS2602, TGS2610, TGS2611 and TGS2620, was used in the investigation on 11 kinds of gases at four different concentrations. The test gases were benzene, toluene, xylene, acetone, butanone, methanol, ethanol, formaldehyde, acetaldehyde, pentane and cyclohexane, while the four concentrations were 100, 200, 300 and 400 ppm. The phase space entire (PSE) feature extraction method extracted six features from PS (Figure 2), which can describe the response curves completely. Suppose S t to be the response of sensors, which is a function of time. Six parameters which can be obtained from the response curve in PS are S 0 (the maximum of S t ), max ( d S t / d t ) / S 0 (the ratio of height and length in adsorption process), a / b (the location of ma ( d S t / d t ) ), min ( d S t / d t ) / S 0 (the ratio of height and length in desorption process), a / b (the location of min ( d S t / d t ) ) and c is wrap value, which details are expressed in [74]. The difference of PSE from other methods based on PS is that the sensor response curves could be reconstructed well with the six features in PSE. PSE was compared with the PSI and maximum value methods and the recognition results of three feature methods with FDA are shown in Table 12.
Figure 2. Six parameters extracted from PS.
Figure 2. Six parameters extracted from PS.
Sensors 15 27804 g002
Table 12. Classification rates (CR) of different features in [74].
Table 12. Classification rates (CR) of different features in [74].
FeaturesNumber of FeaturesCR (%)
Relative difference maximum 190.9
Phase space integral (PSI)181.1
Phase space entire (PSE)6100

2.5. Other Methods for Feature Extraction

Besides the above feature extraction methods, many new original methods have been proposed in recent years. The “fingerprint” obtained from an E-nose is strongly dependent on the kind and concentration of volatiles to which it are exposed, the duration of this exposure and the temperature and humidity of the system. In some cases, it can be appropriate to increase the amount of data in order to capture as much information as possible which properly describes the system, if more reliable models are desired.
A novel method of feature extraction is the use of the three-way structure of the signals obtained from the recording of the different sensor values during a period of time, for each sample. Following this strategy, the contribution of the inner relationship between the three ways can be exploited to obtain more robust information about the system. Parallel factor analysis (PARAFAC) is one of the most popular multi-way data decomposition methods. Originated from psychometrics [75,76], PARAFAC is gaining interest because it is a processing technique which simultaneously determines the pure contributions to the dataset, optimizing each factor as a time, in trilinear systems. In [35], the use of PARAFAC as a feature extraction technique from a specific three-way dataset (samples × sensors × time) for food investigation by E-noses was researched. PARAFAC performed a suitable feature extraction task, incorporating the most of the information of this complex data system. The methodology of PARAFAC also was applied in feature extraction of gas sensor data of monitoring the indoor air pollutants. The results showed that PARAFAC could more accurately monitor and predict the indoor air quality variations than multiway principal component analysis (MPCA), since it is able to capture both the time-correlations and the variable-correlation between indoor air pollutants [36].
Due to the fact that the sensitivities of the gas sensors are different and the signals of the measured samples also are different, the energy characteristics of the signals can be used for feature extraction. An approach called global feature extraction, which focuses on the whole sensors array instead of individual sensors, based on the mutual energy of sensor signal called energy vector (EV) was introduced by Vergara [37]. It is very useful to study the behavior of individual sensors and the theory of signals about the mutual energy is suitable to study the relations between signals of sensors belonging to the same array. The sensors of an array simultaneously interact with the same chemical pattern and as a consequence the sensor signals share a certain degree of correlation. The amount of correlation may change according to the quality and quantity of the materials to which the array is exposed. The EV was defined as:
E = [ ξ 1 , 1 , ξ 1 , 2 , ... , ξ 1 , n , ξ 2 , 1 , ξ 2 , 2 , ... , ξ 2 , n , ... , ξ n 1 , n 1 , ... , ξ n , n ]
where ξ x y = t 0 t 1 x ( t ) y ( t ) d t was the mutual energy between signals x(t) and y(t) in the interval (t0, t1) and n is the number of sensors that compose the array. The method was assessed by solving a practical problem of identification and prediction of pollutant species by building and validating a linear discriminant analysis classifier method performed by partial least squares (PLS-DA). Different validation techniques have been implemented, which have shown that the EV provides sufficient accuracy, and then it is of interest for simple, small size and real time gas analyzers. The results showed that this feature extraction methods based on EV outperformed the FFT method [37].
Kish et al. presented a new feature extraction method for E-noses based on a resistance noise power density spectrum (PDS), which demonstrates that even a single sensor can be sufficient to distinguish between many different chemical species [38]. In this method, it is first assumed that only one sensor is used and its response is independent for each investigated gas. If the PDS of the resistance fluctuations in the sensor has K different frequency ranges, in which the dependence of the response on the concentration of the gases is different from the response in the other ranges, one can write:
d S ( f 1 ) = B 1 , 1 C 1 + B 1 , 2 C 2 + + B 1 , N C N d S ( f K ) = B K , 1 C 1 + B K , 2 C 2 + + B K , N C N
where dS(fi) is the change of the PDS of resistance fluctuations at the i-th characteristic frequency (or frequency range), and the Bi,j quantities are calibration constants in the linear response. Thus, a single sensor is able to provide a set of independent equations to determine the gas composition around the sensor. If we have P different sensors, and if we can use the same characteristic frequency ranges for all sensors, then we have:
d S ( 1 ) ( f 1 ) = B 1 , 1 C 1 + B 1 , 2 C 2 + + B 1 , N C N d S ( 1 ) ( f K ) = B K , 1 C 1 + B K , 2 C 2 + + B K , N C N d S ( 2 ) ( f 1 ) = B 1 , 1 C 1 + B 1 , 2 C 2 + + B 1 , N C N d S ( 2 ) ( f K ) = B K , 1 C 1 + B K , 2 C 2 + + B K , N C N d S ( P ) ( f K ) = B K , 1 C 1 + B K , 2 C 2 + + B K , N C N
where dS(m)(fi) is the change of the PDS of resistance fluctuations at the i-th characteristic frequency range in the m-th sensor and the number of independent equations is P × K .
Wang et al. compared three E-nose features for tomato detection: kurtosis coefficient features, similitude entropy features and the energy features. The results showed that according to the separability measurement, the similitude entropy feature extraction method had its advantages in electronic nose detection [77].
Figure 3. Temporal representation of WTS.
Figure 3. Temporal representation of WTS.
Sensors 15 27804 g003
Gutierrez-Osuna et al. [39] suggested a method called window time slicing (WTS) for creating pseudosensors based upon time slicing of sensor responses, wherein, the time response of each sensor was multiplied by four time windows to obtain the area values shown in Figure 3. First, the response signals are multiplied by four smooth, bell-shaped windowing functions as shown in Figure 3, and then the resulting time traces are integrated with respect to time as shown in Equation (10). In fact, the integrals are the areas surrounded by the response signal curves and the four windowing function curves. Therefore, for each sensor, four “areas” are obtained, which represent the area under the two curves. The windows are designed to overlap and span the entire response interval. These areas were further used as features and were selected using GA for classification of fragrances, hog farm air and cola beverage samples [78]:
A = [ W 1 , W 2 , , W 4 ] W i = k = 1 N R ( t k ) K i ( t k ) Δ t K i ( t k ) = 1 1 + ( t k c i a i ) 2 b i
where A is the feature vector extracted from one sensor which consist of four elements, W i is the area surrounded by the response signal curve and the windowing function curves, R ( t k ) is the response value at the time t k , K i ( t k ) is the value of the i-th windowing function at the time t k , k = 1 , 2 , , N denotes the sampling points of each curve, Δ t is the sampling interval between two sampling points, and the parameters ( a i , b i , c i ) define the width, shape and center of the different windowing functions K i ( t k ) .
However, the drawback of this method is the loss of information owing to the gap between each window. In [79], an attempt has been made to identify the optimum positions of the time windows in order to maximize the classification rate. To avoid the drawback of WTS, a bell shaped window defined by normal probability distribution function with mean at its position and standard deviation as 10 s was moved along the time response (shown in Figure 4) of the sensors and simultaneously multiplied with the sensor response to obtain the area values at different positions. It can be proved that moving window time slicing (MWTS) as an approach to extract area values as features can enhance the performance of an E-nose [80]. Comparing MWTS with other techniques viz. peak value, rise time, fall time, DWT, FFT and WTS, it can be seen that, in Table 13, the MWTS method performed better than the others for classification of Kangra orthodox black tea.
Table 13. Classification rates (CR, %) of different features in [80].
Table 13. Classification rates (CR, %) of different features in [80].
FeaturesAll SensorsSensors 2–4Sensors 2 and 3Sensors 2 and 4Sensors 3 and 4
Maximum value72.9575.7275.2274.974.83
Rise time75.553.4341.2534.9536.18
Fall time61.8887.7085.1283.1793.43
DWT86.3294.6795.0794.7094.90
FFT84.7092.1792.0592.0291.97
WTS89.4092.9292.5492.0293.20
SITO-WTS a94.52
SITO-MWTS b97.3596.2093.3096.35
a Considering position 4 of sensor 3, and positions 1 and 2 of sensor; b Considering bin 3 of sensors 2 and 4, and bin 10 of sensor 3.
Guo et al. [58] proposed a novel feature extraction method also based on window functions called moving window function capturing (MWFC) shown in Figure 5. A 64 points window was placed around the peak value and then moved 64 points to the left and right along with the time axis, respectively. Thus three area values surrounded by two curves can be obtained during the moving process and the three area values were chosen as features simultaneously. The width, position, shape of the window were compared and discussed, and the final results showed that the MWFC feature performed better than the contrast features such as peak value, rising slope, descending slope, FFT, DWT and WFC, showing in Table 14. The same conclusion can be drawn for another two E-nose datasets [81,82].
Figure 4. Schematic diagram of MWTS technique.
Figure 4. Schematic diagram of MWTS technique.
Sensors 15 27804 g004
Figure 5. The schematic diagram of MWFC.
Figure 5. The schematic diagram of MWFC.
Sensors 15 27804 g005
Table 14. Classification rates (CR) of different features in [58].
Table 14. Classification rates (CR) of different features in [58].
Feature ExtractionAccuracy Rate
Peak value87.5
Rising slope87.5
Descending slope85
FFT90.0
DWT92.5
WFC95.0
MWFC97.5

3. Conclusions and Outlook

E-noses have been utilized in a wide and diverse range of applications such as odor analysis, quality control in the food industry, environmental protection detection, public health, explosives detection and spaceflight applications to improve effectiveness, efficiency and safety. Due to the complexity of biological olfaction, the artificial olfactory system has non-linearity characteristics. The characteristic response of the sensor array is a whole cognition of the odors called “fingerprint”. It is crucial to extract useful and robust information from the sensor characteristic response with less redundancy. This paper has presented an overview of the methods of feature extraction used in E-noses for different applications in recent years, such as extraction from the original response curves, curve fitting parameters, transform domains, phase space and dynamic moments, parallel factor analysis, energy vector, power density spectrum and window functions, etc. A few additional approaches, not covered in this review due to space constraints, also offer promising results. Feature extraction from the sensor signals is a key procedure to further improve the performance of an E-nose, but the evaluation of a feature extraction method is influenced by the type of sensors, parameters of experiments, detection targets, demands of specific application, and so on. Different methods are suitable for different situations, and we have to choose the method according to the actual conditions. Although there are no widely approved evaluation criteria for various features, we can give some advice on feature extraction according to the previous research. On the whole, for steady-state features, normalization models give better performance than other models. Transient features have more information than steady-state features. Integrals usually give better performance than maximum values, and the same conclusion can be drawn for derivatives, especially when one sensor extracts several derivatives. Curve fitting parameter features cannot guarantee better recognition results than those of the simple and piecemeal features combination, which gives comparable results to more complex curve fitting parameters if chosen carefully. The single-exponential model is a fine choice, which is simple and easy to deploy and can usually obtain better performance than steady-state features. DWT generally obtains better results than any other feature extraction. In addition, we also hope that more effective feature extraction strategies such as feature fusion can be proposed in further work for enhance the performance of E-noses with the inspiration provided in this paper. Due to the differences of selectivity, sensitivity and specificity of sensors, the optimal features of sensors are various. Usually, the same features are extracted from different sensors with the special methods which we reviewed above. However, it is possible that the features are not optimal for each sensor and better performance can be obtained if different features are extracted and fused for pattern recognition with different feature methods for different sensors.

Acknowledgments

The work was supported by National Natural Science Foundation of China (Grant Nos. 61372139, 61571372, 61101233, 60972155), Program for New Century Excellent Talents in University (Grant Nos. [2013]47), “Spring Sunshine Plan” Research Project of Ministry of Education of China (Grant No. z2011148), Fundamental Research Funds for the Central Universities (Grant Nos. SWU113068, XDJK2014A009, XDJK2013B011, XDJK2014C016, XDJK2016A001), Program for Excellent Talents in scientific and technological activities for Overseas Scholars, Ministry of Personnel in China (Grant No. 2012-186), University Excellent Talents Supporting Foundations in of Chongqing (Grant No. 2011-65), University Key Teacher Supporting Foundations of Chongqing (Grant No. 2011-65), Training Programs for Innovation and Entrepreneurship of China(Grant Nos. 201410635007, 201510635059).

Author Contributions

Shukai Duan is the group leader and he was responsible for the project management and in charge of revising this manuscript. Jia Yan was in charge of the preparation and revising of this manuscript. Xiuzhen Guo was in charge of data collection and revising of this manuscript. Pengfei Jia and Lidan Wang provided valuable advice about the revised manuscript. Chao Peng and Songlin Zhang were involved in discussions and data collection.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gardner, J.W.; Bartlett, P.N. A Brief History of Electronic Noses. Sens. Actuators B Chem. 1994, 18/19, 211–220. [Google Scholar] [CrossRef]
  2. Chen, Q.; Liu, A.; Zhao, J.; Ouyang, Q. Classification of tea category using a portable electronic nose based on an odor imaging sensor array. J. Pharm. Biomed. Anal. 2013, 84, 77–83. [Google Scholar] [CrossRef] [PubMed]
  3. Phaisangittisagul, E.; Nagle, H.T. Predicting odor mixture's responses on machine olfaction sensors. Sens. Actuators B Chem. 2011, 155, 473–482. [Google Scholar] [CrossRef]
  4. Zhang, B.; Huang, Y.; Zhang, Q.; Liu, X.; Li, F.; Chen, K. Fragrance discrimination of Chinese Cymbidium species and cultivars using an electronic nose. Sci. Hortic. 2014, 172, 271–277. [Google Scholar] [CrossRef]
  5. Konduru, T.; Rains, G.C.; Li, C. Detecting sour skin infected onions using a customized gas sensor array. J. Food Eng. 2015, 160, 19–27. [Google Scholar] [CrossRef]
  6. Dai, Y.; Zhi, R.; Zhao, L.; Gao, H.; Shi, B.; Wang, H. Longjing tea quality classification by fusion of features collected from E-nose. Chemometr. Intell. Lab. 2015, 144, 63–70. [Google Scholar] [CrossRef]
  7. Längkvist, M.; Coradeschi, S.; Loutfi, A.; Rayappan, J.B.B. Fast Classification of meat spoilage markers using nanostructured ZnO thin films and unsupervised feature learning. Sensors 2013, 13, 1578–1592. [Google Scholar] [CrossRef] [PubMed]
  8. Musatov, V.Y.; Sysoev, V.V.; Sommer, M.; Kiselev, I. Assessment of meat freshness with metal oxide sensor microarray electronic nose: A practical approach. Sens. Actuators B Chem. 2010, 144, 99–103. [Google Scholar] [CrossRef]
  9. Hartyani, P.; Dalmadi, I.; Knorr, D. Electronic nose investigation of Alicyclobacillus acidoterrestris inoculated apple and orange juice treated by high hydrostatic pressure. Food Control 2013, 32, 262–269. [Google Scholar] [CrossRef]
  10. Pan, L.; Zhang, W.; Zhu, N.; Mao, S.; Tu, K. Early detection and classification of pathogenic fungal disease in post-harvest strawberry fruit by electronic nose and gas chromatography—Mass spectrometry. Food Res. Int. 2014, 62, 162–168. [Google Scholar] [CrossRef]
  11. Romain, A.C.; Nicolas, J. Long term stability of metal oxide-based gas sensors for E-nose environmental applications: An overview. Sens. Actuators B Chem. 2010, 146, 502–506. [Google Scholar] [CrossRef]
  12. Cesare, F.D.; Pantalei, S.; Zampetti, E.; Macagnano, A. Electronic nose and SPME techniques to monitor phenanthrene biodegradation in soil. Sens. Actuators B Chem. 2008, 131, 63–70. [Google Scholar] [CrossRef] [Green Version]
  13. Becherb, C.; Kaulb, P.; Mitrovicsa, J.; Warmerb, J. The detection of evaporating hazardous material released from moving sources using a gas sensor network. Sens. Actuators B Chem. 2010, 146, 513–520. [Google Scholar] [CrossRef]
  14. Li, D.; Lei, T.; Zhang, S.; Shao, X.; Xie, C. A novel headspace integrated E-nose and its application in discrimination of Chinese medical herbs. Sens. Actuators B Chem. 2015, 221, 556–563. [Google Scholar] [CrossRef]
  15. Wilson, A.D. Review of Electronic-nose Technologies and algorithms to detect hazardous chemicals in the environment. Procedia Tech. 2012, 1, 453–463. [Google Scholar] [CrossRef]
  16. He, Q.; Yan, J.; Shen, Y.; Bi, Y.; Ye, G.; Tian, F.; Wang, Z. Classification of electronic nose data in wound infection detection based on PSO-SVM combined with wavelet transform. Int. Autom. Soft Comput. 2012, 18, 967–979. [Google Scholar] [CrossRef]
  17. Adiguzel, Y.; Kulah, H. Breath sensors for lung cancer diagnosis. Biosens. Bioelectron. 2015, 65, 121–138. [Google Scholar] [CrossRef] [PubMed]
  18. Bruins, M.; Rahim, Z.; Bos, A.; Sande, W.; Endtz, H.; Belkum, A. Diagnosis of active tuberculosis by E-nose analysis of exhaled air. Tuberculosis 2013, 93, 232–238. [Google Scholar] [CrossRef] [PubMed]
  19. Burfeind, O.; Bruins, M.; Bos, A.; Sannmann, I.; Voigtsberger, R.; Heuwieser, W. Diagnosis of acute puerperal metritis by electronic nose device analysis of vaginal discharge in dairy cows. Theriogenology 2014, 82, 64–70. [Google Scholar] [CrossRef] [PubMed]
  20. Norman, A.; Stam, F.; Morrissey, A.; Hirschfelder, M.; Enderlein, D. Packaging effects of a novel explosion-proof gas sensor. Sens. Actuators B Chem. 2003, 95, 287–290. [Google Scholar] [CrossRef]
  21. Young, R.C.; Buttner, W.J.; Linnell, B.R.; Ramesham, R. Electronic nose for space program applications. Sens. Actuators B Chem. 2003, 93, 7–16. [Google Scholar] [CrossRef]
  22. Stetter, J.R.; Penrose, W.R. Understanding chemical sensors and chemical sensor arrays (electronic noses): Past, present, and future. Sensors 2002, 10, 189–229. [Google Scholar] [CrossRef]
  23. Wilson, D.M.; Garrod, S.; Hoyt, S.; Kennoch, S.M. Array optimization and preprocessing techniques for chemical sensing microsystems. Sensors 2002, 10, 77–106. [Google Scholar] [CrossRef]
  24. Hierlemann, A.; Berberich, M.S.; Weimar, U.; Kraus, G.; Pfau, A.; Gopel, W. Pattern recognition and multicomponent analysis. Sensors 1996, 2, 119–180. [Google Scholar] [CrossRef]
  25. Gardner, J.W.; Boilot, P.; Hines, E.L. Enhancing electronic performance by sensor selection using a new integer-based genetic algorithm approach. Sens. Actuators B Chem. 2005, 106, 114–121. [Google Scholar] [CrossRef]
  26. Llobet, E.; Brezmes, J.; Vilanova, X.; Sueiras, J.E.; Correig, X. Qualitative and quantitative analysis of volatile organic compounds using transient and steady-state responses of a thick-film tin oxide gas sensor array. Sens. Actuators B Chem. 1997, 41, 13–21. [Google Scholar] [CrossRef]
  27. Carmel, L.; Levy, S.; Lancet, D.; Harel, D. A feature extraction method for chemical sensors in electronic noses. Sens. Actuators B Chem. 2003, 93, 67–76. [Google Scholar] [CrossRef]
  28. Llobet, E.; Brezmes, J.; Ionescu, R.; Vilanova, X.; Al-Khalifa, S.; Gardner, J.W.; Barsan, N.; Correig, X. Wavelet transform fuzzy ARTMAP based pattern recognition for fast gas identification using a micro-hotplate gas sensor. Sens. Actuators B Chem. 2002, 83, 238–244. [Google Scholar] [CrossRef]
  29. Sankaran, S.; Khot, L.R.; Panigrahi, S. Biology and applications of olfactory sensing system: A review. Sens. Actuators B Chem. 2012, 171/172, 1–17. [Google Scholar] [CrossRef]
  30. Lamagna, A.; Reich, S.; Rodriguez, D.; Scoccola, N.N. Performance of an E-nose in hops classification. Sens. Actuators B Chem. 2004, 102, 278–283. [Google Scholar] [CrossRef]
  31. Yang, S.; Larsen, P.K.; Alkjær, T.; Juul-Kristensen, B.; Simonsen, E.B.; Lynnerup, N. Height estimations based on eye measurements throughout a gait cycle. Forensic Sci. Int. 2014, 236, 170–174. [Google Scholar] [CrossRef] [PubMed]
  32. Tian, F.; Xu, X.; Shen, Y.; Yan, J.; Ma, J.; Liu, T. Detection of wound pathogen by an intelligent electronic nose. Sensor Mater. 2009, 21, 155–166. [Google Scholar]
  33. Martinelli, E.; Falconi, C.; Amico, A.D.; Natale, C.D. Feature extraction of chemical sensors in phase space. Sens. Actuators B Chem. 2003, 95, 132–139. [Google Scholar] [CrossRef]
  34. Martinelli, E.; Pennazza, G.; Di Natale, C.; D’Amico, A. Chemical sensors clustering with the dynamic moments approach. Sens. Actuators B Chem. 2004, 101, 346–352. [Google Scholar] [CrossRef]
  35. Padilla, M.; Montoliu, I.; Pardo, A.; Perera, A.; Marco, S. Feature extraction on three way E-nose signals. Sens. Actuators B Chem. 2006, 116, 145–150. [Google Scholar] [CrossRef]
  36. Lee, S.C.; Liu, H.; Kim, M.; Kim, J.T.; Yoo, C. Online monitoring and interpretation of periodic diurnal and seasonal variations of indoor air pollutants in a subway station using parallel factor analysis (PARAFAC). Energy Build. 2014, 68, 87–98. [Google Scholar] [CrossRef]
  37. Vergara, A.; Martinelli, E.; Llobet, E.; Giannini, F.; Amico, A.D.; Natale, C.D. An alternative global feature extraction of temperature modulated micro-hotplate gas sensors array using an energy vector approach. Sens. Actuators B Chem. 2007, 124, 352–359. [Google Scholar] [CrossRef]
  38. Kish, L.; Vajtai, R.; Granqvist, C. Extracting information from noise spectra of chemical sensors: single sensor electronic noses and tongues. Sens. Actuators B Chem. 2000, 71, 55–59. [Google Scholar] [CrossRef]
  39. Ricardo, G.; Troy, N.H. A method for evaluating data-preprocessing techniques for odor classification with an array of gas sensors. IEEE Trans. Syst. Man. Cybern. B Cybern. 1999, 29, 626–632. [Google Scholar] [CrossRef] [PubMed]
  40. Diego, L.; Garcia, A. Detection of defective virgin olive oils by metal-oxide sensors. Eur. Food Res. Technol. 2002, 215, 118–123. [Google Scholar] [CrossRef]
  41. Kermani, B.G.; Schiffman, S.S.; Troy, H. Performance of the Levenberg—Marquardt neural network training method in electronic nose applications. Sens. Actuators B Chem. 2005, 110, 13–22. [Google Scholar] [CrossRef]
  42. Ohnishi, M.; Ishibashi, T.; Kijima, Y.; Ishimoto, C.; Seto, J. A molecular recognition system for odorants incorporating biomimetic gas-sensitive devices using Langmuir-Blodgett films. Sens. Mater. 1992, 4, 53–60. [Google Scholar]
  43. Egashira, M.; Shimizu, Y.; Takao, Y. Trimethylamine sensor based on semiconductive metal oxides for detection of fish freshness. Sens. Actuators 1990, B1, 108–112. [Google Scholar] [CrossRef]
  44. Pearce, T.C.; Gardner, J.W.; Friel, S.; Bartlett, P.N.; Blair, N. Electronic nose for monitoring the flavour of beers. Analyst 1993, 118, 371–377. [Google Scholar] [CrossRef]
  45. Moore, S.W.; Gardner, J.W.; Hines, E.L.; Gopel, W.; Weimar, U. A modified multilayer perceptron model for gas mixture analysis. Sens. Actuators B Chem. 1993, 15/16, 344–348. [Google Scholar] [CrossRef]
  46. Gardner, J.W.; Hines, E.L.; Wilkinson, M. The application of artificial neural networks in an electronic nose. Meas. Sci. Technol. 1990, 1, 446–451. [Google Scholar] [CrossRef]
  47. Gardner, J.W. Detection of vapours and odours from a multisensor array using pattern recognition Part 1: Principal components and cluster analysis. Sens. Actuators B Chem. 1991, 4, 109–115. [Google Scholar] [CrossRef]
  48. Gardner, J.W.; Shurmer, H.V.; Tan, T.T. Application of an electronic nose to the discrimination of coffee. Sens. Actuators B Chem. 1992, 6, 71–75. [Google Scholar] [CrossRef]
  49. Balasubramanian, S.; Panigrahi, S.; Logue, C.; Gub, H.; Marchello, M. Neural networks-integrated metal oxide-based artificial olfactory system for meat spoilage identification. J. Food Eng. 2009, 91, 91–98. [Google Scholar] [CrossRef]
  50. Distantea, C.; Leo, M.; Sicilianoa, P.; Persaud, K. On the study of feature extraction methods for an electronic nose. Sens. Actuators B Chem. 2002, 87, 274–288. [Google Scholar] [CrossRef]
  51. Roussel, S.; Forsberg, G.; Steinmetz, V.; Grenier, P.; Maurel, B.V. Optimisation of Electronic Nose Measurements. Part I: Methodology of Output Feature Selection. J. Food Eng. 1998, 37, 207–222. [Google Scholar] [CrossRef]
  52. Eklov, T.; Martensson, P.; Lundstrom, I. Enhanced selectivity of MOSFET gas sensors by systematical analysis of transient parameters. Anal. Chim. Acta 1997, 353, 291–300. [Google Scholar] [CrossRef]
  53. Zhang, S.; Xie, C.; Zeng, D.; Zhang, Q.; Li, H.; Bi, Z. A feature extraction method and a sampling system for fast recognition of flammable liquids with a portable E-nose. Sens. Actuators B Chem. 2007, 124, 437–443. [Google Scholar] [CrossRef]
  54. Webb, A.R. Statistical Pattern Recognition, 2nd ed.; Wiley: New York, NY, USA, 2002; pp. 243–245. [Google Scholar]
  55. Šetkus, A.; Olekas, A.; Senuliene, D.; Falasconi, M.; Pardo, M.; Sberveglieri, G. Analysis of the dynamic features of metal oxide sensors in response to SPME fiber gas release. Sens. Actuators B Chem. 2010, 146, 539–544. [Google Scholar] [CrossRef]
  56. Zou, X.; Zhao, J.; Wu, S.; Huang, X. Vinegar Classification Based on Feature Extraction and Selection from Tin Oxide Gas Sensor Array Data. Sensors 2003, 3, 101–109. [Google Scholar] [CrossRef]
  57. Wei, Z.; Wang, J.; Zhang, W. Detecting internal quality of peanuts during storage using electronic nose responses combined with physicochemical methods. Food Chem. 2015, 177, 89–96. [Google Scholar] [CrossRef] [PubMed]
  58. Guo, X.; Peng, C.; Zhang, S.; Yan, Jia.; Duan, S.; Wang, L.; Jia, P.; Tian, F. A Novel Feature Extraction Approach Using Window Function Capturing and QPSO-SVM for Enhancing Electronic Nose Performance. Sensors 2015, 15, 15198–15217. [Google Scholar] [CrossRef] [PubMed]
  59. Galdikas, A.; Mironas, A.; Senuliene, D.; Strazdiene, V.; Setkus, A.; Zelenin, D. Response time based output of metal oxide gas sensors applied to evaluation of meat freshness with neural signal analysis. Sens. Actuators B Chem. 2000, 69, 258–265. [Google Scholar] [CrossRef]
  60. Duda, R.O.; Hart, P.E.; Stork, D.G. Pattern Classification, 2nd ed.; Wiley: New York, NY, USA, 2000. [Google Scholar]
  61. Holmberg, M.; Gustafsson, F.; Hornsten, E.G.; Winquist, F.; Nilsson, L.E.; Ljung, I.; Lundstro, L. Bacteria classification based on feature extraction from sensor data. Biotechnol. Tech. 1998, 12, 319–324. [Google Scholar] [CrossRef]
  62. Yan, J.; Tian, F.; He, Q.; Shen, Y.; Xu, S.; Feng, J.; Chaibou, K. Feature Extraction from Sensor Data for Detection of Wound Pathogen Based on Electronic Nose. Sens. Mater. 2012, 24, 57–73. [Google Scholar]
  63. Daubechies, I. Ten Lectures on Wavelets; Society for Industrial and Applied Mathematics: Philadelphia, PA, USA, 1992. [Google Scholar]
  64. Feng, J.; Tian, F.; Yan, J.; He, Q.; Shen, Y.; Pan, L. A background elimination method based on wavelet transform in wound infection detection by electronic nose. Sens. Actuators B Chem. 2011, 157, 395–400. [Google Scholar] [CrossRef]
  65. Ratton, L.; Kunt, T.; McAvoy, T.; Fuja, T.; Cavicchi, R.; Semancik, S. A comparative study of signal processing techniques for clustering microsensor data (a first step towards an artificial nose). Sens. Actuators B Chem. 1997, 41, 105–120. [Google Scholar] [CrossRef]
  66. Ionescu, R.; Llobet, E. Wavelet transform-based fast feature extraction from temperature modulated semiconductor gas sensors. Sens. Actuators B Chem. 2002, 81, 289–295. [Google Scholar] [CrossRef]
  67. Fang, Y.; Chen, M.; Zheng, X. Extracting features from phase space of EEG signals in brain—Computer interfaces. Neurocomputing 2015, 151, 1477–1485. [Google Scholar] [CrossRef]
  68. Samadani, M.; Kitio Kwuimy, C.A.; Nataraj, C. Model-based fault diagnostics of nonlinear systems using the features of the phase space response. Commun. Nonlinear Sci. 2015, 20, 583–593. [Google Scholar] [CrossRef]
  69. Su, Z.; Wu, T.; Yang, P.; Wang, Y. Dynamic analysis of heartbeat rate signals of epileptics using multidimensional phase space reconstruction approach. Phys. A 2008, 387, 2293–2305. [Google Scholar] [CrossRef]
  70. Takens, F. Detecting Strange Attractors in Turbulence. In Dynamical Systems and Turbulence; Rand, D., Young, L.S., Eds.; Springer: Warwick, UK, 1980; pp. 366–381. [Google Scholar]
  71. Vergara, A.; Llobet, E.; Martinelli, E.; Natale, C.D.; Amico, A.D.; Correig, X. Feature extraction of metal oxide gas sensors using dynamic moments. Sens. Actuators B Chem. 2007, 122, 219–226. [Google Scholar] [CrossRef]
  72. Fichera, A.; Losenno, C.; Pagano, A. Clustering of chaotic dynamics of a lean gas-turbine combustor. Appl. Energy 2001, 69, 101–117. [Google Scholar] [CrossRef]
  73. Llop, M.F.; Jand, N.; Gallucci, K.; Llauro, F.X. Characterizing gas—Solid fluidization by nonlinear tools: Chaotic invariants and dynamic moments. Chem. Eng. Sci. 2012, 71, 252–263. [Google Scholar] [CrossRef]
  74. Zhang, S.; Xi, C.; Hu, M.; Li, H.; Bai, Z.; Zeng, D. An entire feature extraction method of metal oxide gas sensors. Sens. Actuators B Chem. 2008, 132, 81–89. [Google Scholar] [CrossRef]
  75. Harshman, R.A. Foundations of the PARAFAC procedure: Model and conditions for an ‘explanatory’ multi-mode factor analysis. UCLA Work. Pap. Phon. 1970, 16, 1–84. [Google Scholar]
  76. Carroll, J.D.; Chang, I. Analysis of individual differences in multidimensional scaling via an N-way generalization of and Eckart-Young decomposition. Psychometrika 1970, 35, 283–319. [Google Scholar] [CrossRef]
  77. Wang, W.; Guan, X.; Wang, J.; Yang, Y.; Liu, H. The research of feature extraction methods in the tomatoes detection. In Proceedings of the IEEE International Conference on Instrumentation & Measurement, Computer, Communication and Control, Harbin, China, 23 December 2012; pp. 702–705.
  78. Kermani, B.; Schiffman, S.; Nagle, H. Using Neural Networks and Genetic Algorithms to Enhance Performance in an Electronic Nose. IEEE Tans. Biomed. Eng. 1999, 46, 429–439. [Google Scholar] [CrossRef]
  79. Phaisangittisagul, E. Signal Processing Using Wavelets for Enhancing Electronic Nose Performance. Ph.D. Thesis, North Carolina State University, Raleigh, NC, USA, 2007. [Google Scholar]
  80. Kaur, R.; Kumar, R.; Gulati, A.; Ghanshyam, C.; Kapur, P.; Bhondekar, A.P. Enhancing electronic nose performance: A novel feature selection approach using dynamic social impact theory and moving window time slicing for classification of Kangra orthodox black tea (Camellia sinensis (L.) O. Kuntze). Sens. Actuators B Chem. 2012, 166/167, 309–319. [Google Scholar] [CrossRef]
  81. Xu, X. Research on the Key Technology of Medical Electronic Nose. Ph.D. Thesis, Chong Qing University, Chongqing, China, 2009. [Google Scholar]
  82. Dang, L.; Tian, F.; Zhang, L.; Kadria, C.; Yin, X.; Peng, X.; Liu, S. A novel classifier ensemble for recognition of multiple indoor air contaminants by an electronic nose. Sens. Actuators B Chem. 2014, 207, 67–74. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Yan, J.; Guo, X.; Duan, S.; Jia, P.; Wang, L.; Peng, C.; Zhang, S. Electronic Nose Feature Extraction Methods: A Review. Sensors 2015, 15, 27804-27831. https://doi.org/10.3390/s151127804

AMA Style

Yan J, Guo X, Duan S, Jia P, Wang L, Peng C, Zhang S. Electronic Nose Feature Extraction Methods: A Review. Sensors. 2015; 15(11):27804-27831. https://doi.org/10.3390/s151127804

Chicago/Turabian Style

Yan, Jia, Xiuzhen Guo, Shukai Duan, Pengfei Jia, Lidan Wang, Chao Peng, and Songlin Zhang. 2015. "Electronic Nose Feature Extraction Methods: A Review" Sensors 15, no. 11: 27804-27831. https://doi.org/10.3390/s151127804

APA Style

Yan, J., Guo, X., Duan, S., Jia, P., Wang, L., Peng, C., & Zhang, S. (2015). Electronic Nose Feature Extraction Methods: A Review. Sensors, 15(11), 27804-27831. https://doi.org/10.3390/s151127804

Article Metrics

Back to TopTop