<figure><img src="image_1.png" /></figure>

# Comparison between SVM and KNN classifiers for iris recognition using a new unsupervised neural approach in segmentation

Hicham Ohmaid¹, S. Eddarouich², A. Bourouhou³, M. Timouyas⁴

¹National School of Computer Science and Systems Analysis (ENSIAS), Mohammed V University, Rabat/10112, Morocco

²Regional Educational Center, Rabat/10112, Morocco

³,⁴Higher Normal School of Technical Education (ENSET), Mohammed V University, ENSET, Rabat/10112, Morocco

**Article Info**

**Article history:**

Received Dec 12, 2019

Revised Feb 17, 2020

Accepted Mar 19, 2020

**Keywords:**

Biometric

Discrete wavelet transformation

Hough transform

Iris segmentation

K-nearest-neighbor (KNN)

Support vector machine (SVM)

kernel function

Unsupervised neural approach

**ABSTRACT**

A biometric system of identification and authentication provides automatic recognition of an individual based on certain unique features or characteristics he or she possesses. Iris recognition is a biometric identification method that applies pattern recognition to images of the iris. Owing to the unique epigenetic patterns of the iris, iris recognition is considered one of the most accurate methods in the field of biometric identification. The segmentation algorithm proposed in this article starts with determining the regions of the eye using unsupervised neural approach, after the outline of the eye is found using the Canny edge, The Hough Transform is employed to determine the center and radius of the pupil and the iris. Then the normalization allows transforming the segmented circular iris region into a fixed-size rectangular shape using Daugman's rubber sheet model. A discrete wavelet transformation (DWT) is applied to the normalized iris to lower the size of iris models and to improve classifier accuracy. Finally, the URIBIS iris database is used for individual user verification by using the KNN classifier or support vector machine (SVM) which is based on the analysis of iris code as feature extraction is discussed.

This is an open access article under the **CC BY-SA** license.

<figure><img src="image_2.png" /></figure>

**Corresponding Author:**

Hicham Ohmaid,

National School of Computer Science and Systems Analysis (ENSIAS),

Mohammed V University, Rabat/10112, Morocco.

Email: hicham.ohmaid@gmail.com

## 1. INTRODUCTION

Nowadays, biometric recognition has become a reliable way to identify and recognize individuals. Conventional authentication methods like passwords and identification cards are not always credible because these methods can be forgotten or stolen. Identification based on physiological or behavioral characteristics of people (such as irises, fingerprints, faces, hand geometry, signatures, fingerprints, et cetera) can accurately authenticate people's identity [1].

Iris recognition offers the highest accuracy in identifying individuals compared to any other biometric approaches. The iris is defined as the circular part between the pupil and sclera of the eye. It is so unique that no two irises are alike, even among identical twins or between the left and right eye of the same person [2-3]. These visible characteristics are thought to be discriminatory and unique to each eye, as well as stable over an individual's lifetime. This makes the iris a very useful biometric identifier when it is possible to capture iris images of sufficient quality [4].
<figure><img src="image_3.png" /></figure>

The system is inspired by Daugman's work [5]. It is composed of a number of sub-systems, which are described in Figure 1 and Figure 2 the contains stages of an iris detection system.

The segmentation step enables the localization of the iris region in an eye image, which consists of finding two contours of the iris (pupil and iris). In other words, the internal contour (pupil/iris) and outer (iris/sclera) are located. Once the iris region is segmented, the next stage is to normalize this segmentation. The normalization transforms the iris texture from Cartesian to polar coordinates in order to limit the problems of variations in the eye (like the optical size of the iris or the position of the pupil in the iris). This step is based on Daugman's rubber sheet model [5]. Feature extraction is a technique used to extract information from the iris image, in order to lower the size of iris models and improve classifier accuracy. These features can not be used to reconstruct images, but these values are used in classification. Finally, matching is performed by comparing the features of a template iris with the feature vectors of templates in the database, and a decision is formulated [6-7].

<figure><img src="image_4.png" /><figcaption>Figure 1. Iris recognition system</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>Figure 2. Typical stages of iris recognition</figcaption></figure>

The main motivation in this research is to propose a new effective and robust algorithm to segment and classification the clear or noisy iris images. The proposed algorithm added a new pre-processing step using a new unsupervised neural approach to divide the iris image in two regions namely iris and eyelashes region, sclera and skin, in order to facilitate the determination of the iris contour in the phase Canny edge detection [8]. Unlike other iris segmentation algorithms, that process the whole image of the eye (which contains the non-iris regions that generate segmentation errors) such as Daugman algorithms. Then, discrete wavelet transformation (DWT) is used for extracting the optimum features of iris images and reducing the runtime of classification of these iris templates. Finally, a comparison was made between the classifiers SVM and KNN.

## 2. METHODOLOGY

The iris recognition system consists of the following steps:

- Image acquisition

- Iris segmentation

- Normalization

- Feature extraction

- Classification

In the approach presented here, segmentation was achieved using the Hough transform with the aim of locating the iris region of an eye image. The normalization enabled the transformation of the segmented
<figure><img src="image_6.png" /></figure>

circular iris region into a fixed-size rectangular shape using Daugman's rubber sheet model. Two level discrete wavelet transformation (DWT) was applied to the normalized iris for feature extraction. The support vector machine was used to classify the similarity between the iris templates. Figure 3 shows a diagram of the iris recognition system.

<figure><img src="image_7.png" /><figcaption>Figure 3. Diagram of iris recognition system</figcaption></figure>

## 2.1. Iris segmentation

In this section, we present a new iris segmentation approach used as a pre-processing step, based on neural competitive concepts [9]. As shown in Figure 4, this approach allows partitioning the iris image in two regions namely iris and eyelashes region, sclera and skin. Then, the outline of the eye is found using the Canny edge, The Hough Transform is employed to determine the center and radius of the pupil and the iris. Figure 4 shows the iris segmentation steps.

<figure><img src="image_8.png" /><figcaption>Figure 4. Algorithm used for iris recognition</figcaption></figure>
<figure><img src="image_9.png" /></figure>

### 2.1.1. Iris segmentation using the unsupervised neural approach

In an unsupervised context, that means when no prior information about the data sample, we involve the neural competitive clustering procedure to the Iris segmentation. The objective of these approaches to classifying the image pixels according to their distribution in the representation space and to assign them a label. As illustrated in Figure 5, it starts first by the organizationIt starts first, by the organization of the pixels of Iris image in an observation matrix (each row represents a pixel and each column represents an attribute) to estimate the underlying probability density function (pdf) of the pixels distribution using a non-parametric estimator. In the second step, the procedure uses an artificial neural network with competitive training (NNCT) to extract the local maxima of the pdf. Following a modes detection method using a technique to detect the existing interneural connection [10]. The last step is for affecting the remaining pixels to their classes.

<figure><img src="image_10.png" /><figcaption>Figure 5. Architecture of the segmentation procedure</figcaption></figure>

### 2.1.2. The estimation of underlying probability density function

After constructing the observation matrix of an Iris image pixels $\Gamma = \{X_1, X_2, ..., X_Q\}$ considering as a set of Q N-dimensional observations (in this case N=3) With $X_q = [x_{q,1}, x_{q,2}, ..., x_{q,n}, ..., x_{q,N}]$, q=1,2,...,Q and a probability density function P(X), anon-parametric method based on The Parzen [11] window using to estimate this underlying density function. The technique is a fast estimation algorithm that is proposed by Postaire and Vasseur [12]. First, the range of variation of each component of these observations is normalized to the interval [0, R], where R is an integer such as $R \ge 2$, by means of the transformation defined as:

$$y_{n,q} = \frac{\left( \frac{x_{n,q} - \min_{q} x_{n,q}}{\max_{q} x_{n,q} - \min_{q} x_{n,q}} \right)}{\max_{q} x_{n,q} - \min_{q} x_{n,q}} * R \quad (1)$$

Each axis of the so normalized data space is then partitioned into R exclusive intervals of unit width. This discretization defines a set of $R^N$ hypercube of unit side length. Each hypercube noted $H(X)$, is a site defined by its N coordinates $x_1, x_2, ..., x_n, ..., x_N$ which are the integer parts of the coordinates of its centre X. To be more specific, let $Y_q = [y_{1,q}, y_{2,q}, ..., y_{n,q}, ..., y_{N,q}]$, q=1, 2, Q be the Q observations in the normalized space. Each observation Yq is found inside a non-empty hypercube with the coordinates $x_n = \text{int}(y_{n,q})$, n=1, 2, N, where $\text{int}(y_{n,q})$ designates the integer parts of $y_{n,q}$. If several observations fall in the same hypercube, this one appears many times on the list of non-empty hypercubes. Furthermore, the number of times the hypercube H(X) appears in that list indicates the number of data points q[H(X)] which falls in this hypercube. Subsequently, the value of the local density estimated is:

$$p(X) = \frac{q[H(X)]}{Q} \quad (2)$$

Since the volume of H(X) is equal to unity.

So, this fast procedure allows only the estimation of the underlying probability density function at the centers of the non-empty hypercube whose number never exceeds the number Q of available observations (pixels). At the centers of the hypercube cells, which are not on that list, the density estimates are known to be null. At the end of this fast algorithm, all the available information for clustering is in the discrete set $\bar{X}$ of estimated values of the underlying probability density function P(X).
<figure><img src="image_11.png" /></figure>

### 2.1.3. The extraction of local maxima by neural network

Assimilating the modes to the local maxima of the pdf, the proposed approach uses the Neural Networks with Competitive Training (NNCT) [13]. In the training algorithm, we work only on the pdf by presenting sequentially, the centers of the non-empty hypercube of the set X to the network, instead of the Q observations.

The neural network is composed of two layers: the input layer and the output layer. The first one is made of N units $I_n$, n=1,2,...,N, such that unit $I_n$ is solicited by the attribute $X_N$ of the non-empty hypercube $H(X)$ when this one is presented to the network. However, each output neuron materializes a hypercube which represents the site of one local maximum of the pdf in the set X, and presents its weight by the mean vector $\mu_k(X)$, k=1,2,...,K. The number of the output neuron is first initialized arbitrarily.

During the training phase, the output neurons enter into competition with each other by comparing the distance $D[\mu_k(X), H(X)]$, k=1,2,...,K, between the input hypercube $H(X)$ and each output neuron $\mu_k(X)$, the winner is the closest one to the hypercube, then we compare the values of the pdf associated to the winner neuron $\mu_g(X)$ and to $H(X)$. The distance measure used in this training algorithm is Mahalanobis distance that gives the best results for the non-Gaussian distribution [14] instead of Euclidian distance as in the (NNCT).

### 2.1.4. Detection of significant modes of pdf

The aim of this phase is to connect each group of the closest modes in such a way that we get a map which preserves the shape and structure of the classes existing in the image, by an improved Competitive Hebbian Learning method (CHLim) [10] which allows eliminating the influence of the output neural number.

To generate the induced Delaunay triangulation, the CHLim, given the K modes detected by CNN as prototypes in RN, successively adds connections among them. The method does not change the weight of prototypes, but only generates topology according to these prototypes. For each mode $H(X_k)$, can be connected to the two closest prototypes by an edge using Mahalanobis distance as a measure of resemblance instead of Euclidian distance, it works as an activation function for competition between neurons. This leads to the induced Delaunay triangulation, which is limited to those regions of the input space $R^N$.

### 2.1.5. Classification procedure

Once the different modes are identified, the classification method [15] that we use in the present work consists first, in defining the pixels falling into any mode of a connected set as the prototype of one cluster. Then, the remaining pixels, which do not fall in one of the detected modes, are assigned to the clusters attached to their nearest neighbor among these prototypes by means of Mahalanobis distance. Figure 6 presents the procedure for classifying the iris image.

<figure><img src="image_12.png" /><figcaption>Figure 6. Illustration of the Iris segmentation procedure</figcaption></figure>

### 2.1.6. Hough transform

Hough transform is a standard image analysis tool for finding curves that can be defined in a parametrical form such as lines, circles, parabolas, and hyperbolas [16]. The edge map is then used in a
<figure><img src="image_13.png" /></figure>

voting process to maximize the defined Hough transform for the desired contour. Considering the obtained edge points as (xj; yj), j = 1; 2; ......, n, a Hough transform can be written as:

$$H(x_c, y_c, r) = \sum_{j=1}^{n} h(x_j, y_j, x_c, y_c, r) \qquad (3)$$

Where

$$h(x_j, y_j, x_c, y_c, r) = \begin{cases} 1 & \text{if } g(x_j, y_j, x_c, y_c, r) = 0 \\ 0 & \text{otherwise} \end{cases} \qquad (4)$$

The limbus and pupil are both modeled as circles and the parametric function g is defined as:

$$g(x_j, y_j, x_c, y_c, r) = (x_j - x_c)^2 + (y_j - y_c)^2 - r^2 \qquad (5)$$

Assuming a circle with center (xc; yc) and radius r, the edge points that are located over the circle result in a zero value for the function. The value of g is then transformed to 1 by the h function, which represents the local pattern of the contour. The local patterns are then used in a voting procedure using the Hough transform, H, in order to locate the proper pupil and limbus boundaries. In order to detect the limbus, only vertical edge information is used. The two eyelids usually cover the upper and lower parts; these parts contain horizontal edge information. The horizontal edge information is used for detecting the upper and lower eyelids, which are modeled as parabolic arcs [9].

After the Hough Transform process is complete, six parameters are stored: the radius and x and y (the center coordinates for both circles). Eyelids are isolated by first fitting a line to the upper and lower eyelid using the linear Hough transform. A second horizontal line is then drawn that intersects with the first line at the iris edge that is closest to the pupil. This process is completed for both the top and bottom eyelids. The second horizontal line allows maximum isolation of eyelid regions. Figure 7 shows an example of properly segmented iris.

<figure><img src="image_14.png" /><figcaption>Figure 7. Examples of the segmented iris</figcaption></figure>

## 2.2. Iris normalization

The size of the iris may change because of variation in the illumination, even for an iris from the same person. This elastic deformation in iris texture must be compensated for to achieve more accurate recognition results. Normalization is a technique to prepare a segmented iris for feature extraction and transforms the segmented circular iris region into a fixed-size rectangular shape. The Cartesian-to-polar transform of the iris region is based on Daugman's rubber sheet model [5], as illustrated in Figure 8.
<figure><img src="image_15.png" /><figcaption>Figure 8. Cartesian-to-polar coordinate transform</figcaption></figure>

The Cartesian-to-polar transform of the iris region is modeled by (4).

$$I(x(r, \theta), y(r, \theta)) \rightarrow I(r, \theta) \qquad (4)$$

with

$$\begin{aligned} x(r, \theta) &= (1-r)x_p(\theta) + rx_i(\theta) \\ y(r, \theta) &= (1-r)y_p(\theta) + ry_i(\theta) \end{aligned} \qquad (5)$$

The coordinates of the inner and outer boundaries along the θ direction are xp, yp and xi, yi, respectively. Image-normalized enhancement is required to obtain accurate features because the normalized iris presents many variations due to variations in light illuminations. By using histogram equalization, pixels intensities are spread over the entire iris template. Figure 9 shows an iris template with histogram equalization.

<figure><img src="image_16.png" /><figcaption>Figure 9. (a) Iris normalized into polar coordinates, (b) Enhancement of the iris normalized image</figcaption></figure>

## 2.3. Feature extraction

In the encoding stage, two level discrete wavelet transformation (DWT) is applied to the normalized iris region, in order to determine DWT coefficients, by passing the signal in different frequency ranges, namely low-low (LL), low-high (LH), high-low (HL) and high-high (HH), as indicated in Figure 10. The frequency range can be represented as LL < LH < HL < HH [17-18].

The LL sub-band represents the feature characteristics of the iris; this sub-band (LL) is used in the classification phase. Figure 11 indicates the resolution of the normalized and enhanced iris image (48x432). After applying DWT to an enhanced iris image, the resolution of sub-band is (24 × 216). This sub-band is used instead of the original normalized iris data to reduce the resolution of the iris model, and the runtime of the classification is similarly reduced.
<figure><img src="image_17.png" /></figure>

<figure><img src="image_18.png" /><figcaption>Figure 10. Iris template using DWT</figcaption></figure>

<figure><img src="image_19.png" /><figcaption>Figure 11. (a) Enhancement of the iris normalized image, (b) DWT transforms normalized iris image into LL sub-bands</figcaption></figure>

## 2.4. Support vector machine SVM

The support vector machine (SVM) performs pattern recognition according to the principle of structural risk minimization. The development of SVM as a classifier has two major aspects: the first aspect is to find the hyperplane that optimally separates the two classes, and the second aspect is the transformation of the linearly non-separable classification problem into a linearly separable problem [19].

Given an input entity vector learning set and the class tag (x_i,y_i), i=1,2,...,l, where $x_i \in R^n$ and $y \in \{1, -1\}$. In the case of a linearly separable problem, the separating hyperplane (which defines the boundary between class 1 and class 2) can be represented as follows:

$$ \min_{\omega,b,\xi} \frac{1}{2} \omega^T \omega + C \sum_{i=1}^{l} \xi_i \quad (6) $$

Subject to

$$ y_i(\omega^T\phi(x_i) + b) \geq 1 - \xi_i; \xi_i \geq 0 \quad (7) $$

In the case of a non-linearly separable problem, the training vectors x_i are mapped into higher dimensional space by the function $\emptyset$ in order to make the problem linearly separable. The SVM then finds a linear separating hyperplane with the maximal margin in this higher dimensional space. The penalty parameter of the error term is C > 0. Furthermore, $K(x_i, x_j) = \emptyset(x_i)^T(x_j)$ is called the kernel function [20].

- Linear:

$$ K(x_i, x_j) = x_i^T x_j \quad (8) $$

- Polynomial:

$$ K(x_i, x_j) = (\gamma x_i^T x_j + r)^d, \gamma > 0 \quad (9) $$

- Gaussian:

$$ K(x_i, x_j) = \exp(-\frac{\|x_i - x_j\|^2}{2\gamma^2}) \quad (10) $$

## 3. RESULTS

The Classification Learner App (included in the Statistics and Machine Learning Toolbox for MATLAB) is an application that can be used to train models to classify data using supervised machine learning methods. The SVM and KNN classifiers were used for our approach [20-21].
<figure><img src="image_20.png" /></figure>

By using Classification Learner, we can explore data, select features, specify validation schemes, and train and export classification models to the MATLAB workspace to generate predictions for new data.

To evaluate the effectiveness of the proposed system, a set of eye images obtained from the UBIRIS databases were used [22]. Figure 12 shows an example of the eye images contained in this database. The total number of iris classes is 100, and each iris class contains three images stored in the database for this experiment. In this work, 300 images were used for training, and 100 images were used for testing purposes.

<figure><img src="image_21.png" /><figcaption>Figure 12. Eye images in the UBIRIS database</figcaption></figure>

Four kernel functions (Linear, Quadratic, Cubic, and Medium Gaussian) were used, and the classification efficiency of each kernel was measured. Since the highest classification accuracy was obtained by a Medium Gaussian kernel, as shown in Table 1, this kernel was used in our system for classification and recognition purposes. As the simulation result shown in Table 2, the KNN algorithm assured better the accuracy of 95%, as well as the execution time is too low compared to that of the SVM algorithm.

<table><caption>Table 1. Efficiency of various kernel functions (SVM)</caption><thead><tr><td>Kernel Type</td><td>Classification Accuracy (%)</td><td>Training Time (sec)</td></tr></thead><tbody><tr><td>Linear</td><td>68.2</td><td>430.89</td></tr><tr><td>Quadratic</td><td>84.3</td><td>533.09</td></tr><tr><td>Cubic</td><td>88.6</td><td>580.83</td></tr><tr><td>Medium Gaussian</td><td>88.9</td><td>602.28</td></tr></tbody></table>

<table><caption>Table 2. Efficiency of KNN classifier</caption><thead><tr><td>Classification Accuracy (%)</td><td>Training Time (sec)</td></tr></thead><tbody><tr><td>KNN 95.00</td><td>19.94</td></tr></tbody></table>

The accuracy of our proposed technique is 95%, is shown in Table 3, which is better than the other methods that are used for runtime comparison.

<table><caption>Table 3. Accuracy comparison of our proposed technique with others</caption><thead><tr><td>Author</td><td>Classification Accuracy (%)</td></tr></thead><tbody><tr><td>H Ohmaid et al. (Our proposed)</td><td>95</td></tr><tr><td>MH Hamd and al. [23]</td><td>94</td></tr><tr><td>SG Firake and al. [24]</td><td>92.85</td></tr><tr><td>J poonia and al. [25]</td><td>82.5</td></tr></tbody></table>

## 4. CONCLUSIONS

In this paper, segmentation was achieved using pretreatment step based on an unsupervised neural approach. The normalization enabled the transformation of the segmented circular iris region into a fixed-size rectangular shape using Daugman's rubber sheet model. Discrete wavelet transformation (DWT) is used for extracting the optimum features of iris images and reducing the runtime of classification of these iris templates. Finally, a comparison was made between the classifiers SVM and KNN. Experimental evaluation using the URIBIS database clearly demonstrated the proposed technique performs in an efficient manner, especially for the KNN classifier that gave an accuracy of 95%.
<figure><img src="image_22.png" /></figure>

REFERENCES

[1] Fang, Bin Et Tang, Yuan Yan. Improved class statistics estimation for sparse data problems in offline signature verification. *IEEE Transactions on Systems, Man, and Cybernetics, Part C (Applications and Reviews)*. vol. 35, no 3, pp. 276-286, 2005.

[2] Gao, Xinbo, Zhong, Juanjuan, LI, Jie, et al. Face sketch synthesis algorithm based on E-HMM and selective ensemble. *IEEE Transactions on Circuits and Systems for Video Technology*, vol. 18, no 4, p. 487-496, 2008.

[3] Yu, Li, Zhang, David, Et Wang, Kuanquan. The relative distance of key point based iris recognition. *Pattern Recognition*, vol. 40, no 2, p. 423-430, 2007.

[4] Zhang, Taiping, Fang, Bin, Liu, Weining, et al. Total variation norm-based nonnegative matrix factorization for identifying discriminant representation of image patterns. *Neurocomputing*, vol. 71, no 10, p. 1824-1831, 2008.

[5] Daugman, J. Biometric personal identification system based on iris analysis, United States Patent, US 5291560. 1994.

[6] Daugman, J., "How iris recognition works," Image Processing. 2002. Proceedings. 2002 International Conference on, vol.1, no., pp.I-33, I-36 vol.1, 2002.

[7] Masek, Libor, et al. Recognition of human iris patterns for biometric identification. 2003.

[8] Ohmaid, Hicham, Eddarouich, S., Bourouhou, A., et al. "Iris segmentation using a new unsupervised neural approach". *Int J Artif Intell (IJAI)*, vol. 9, no 1, pp.58-64, 2020.

[9] H. Ohmaid, M. Timouyas and S. Eddarouich, "A comparative study of similarity measurement in a new neural model of Unsupervised clustering," 2019 Third International Conference on Intelligent Computing in Data Sciences (ICDS), Marrakech, Morocco, 2019, pp. 1-8, doi: 10.1109/ICDS47004.2019.8942242.

[10] Timouyas M., Eddarouich S. and Hammouch A. "Mode region detection using improved Competitive Hebbian Learning for unsupervised clustering," *Engineering Science and Technology: An International Journal (ESTIJ)*, ISSN: 2250-3498 Vol.7, No.4, pp. 26-35, 2017.

[11] Parzen, E. "An Estimation of a Probability Density Function and Mode," *Ann. Math. Stat.*, vol. 33, pp.1065-1076, 1962.

[12] Postaire, J.-G., & Vasseur, C. P. A. "A fast Algorithm for non Parametric Probability Density Estimation," *IEEE, Trans. on Pattern Anal. and Machine Intel*. PAMI-4, n 6, pp. 663-666, 1982.

[13] Eddarouich, S., & Sbihi, A. "Neural Network for Modes Detection in Pattern Classification." ICTIS'07, Morocco, Fez, 3-5 pp. 300-303, 2007.

[14] Timouyas, M., Eddarouich, S., and Hammouch, A. "A new approach of classification for non-Gaussian distribution upon competitive training," (ICCS'12), Agadir, Morocco, pp.1-6, 2012.

[15] Eddarouich, S., and Sbihi, A. "Détection des Modes par Approche Neuronale pour la Classification des Données d'un Mélange des Distributions Normales," *Proceeding of the 6th African Conference on Research in Computer Science (CARI'02)*, Yaoundé, Cameroon, pp. 61-68, 2002.

[16] Wildes, Richard P. Iris recognition: an emerging biometric technology. *Proceedings of the IEEE*, vol. 85, no 9, p. 1348-1363, 1997.

[17] Ma, Li, Tan, Tieniu, Wang, Yunhong, et al. Efficient iris recognition by characterizing key local variations. *IEEE Transactions on Image processing*, vol. 13, no 6, p. 739-750, 2004.

[18] Ganorkar, Sanjay Et Memane, Mayuri. "Iris recognition using discrete wavelet Transform." *International Journal of Advances in Engineering & Technology*, vol. 4, no 1, p. 356, 2012.

[19] Rana, Humayan Kabir, Azam, Md Shafiul, Akhtar, Mst Rashida, et al. "A fast iris recognition system through optimum feature extraction". *Peer J Computer Science*, vol. 5, p. e184, 2019.

[20] Czajka, A., Bowyer, K. W., Krumdick, M., & VidalMata, R. G. "Recognition of image-orientation-based iris spoofing". *IEEE Transactions on information Forensics and Security*, 12(9), 2184-2196, 2017.

[21] Zhang, Hui and Guan, Xiangfeng. Iris recognition based on grouping KNN and rectangle conversion. *In: 2012 IEEE International Conference on Computer Science and Automation Engineering*. IEEE. p. 131-134, 2012.

[22] "Ubiris: Noisy Visible Wavelength Iris Image Databases." [Online]. Available: URL:http://iris.di.ubi.pt/.

[23] Hamd, Muthana H. and Ahmed, Samah K. Biometric system design for iris recognition using intelligent algorithms. *International Journal of Modern Education and Computer Science*, vol. 11, no 3, p. 9, 2018.

[24] Firake, Sarita G. and Mahajan, P. M. Comparison of iris recognition using gabor wavelet, principal component analysis and independent component analysis. *International Journal of Innovative Research in Computer and Communication Engineering*, vol. 4, no 6, p. 12334-12342, 2016.

[25] Jyoti, P., Parvati, B., Sandeep, K. G., and al. New improved feature extraction approach of IRIS recognition. *International Journal of Computer Systems*, vol. 3, no 1, p. 1-3, 2016.