Tecnical features of the dabase.
As is defined by UNESCO : "Biological diversity or biodiversity is defined as the diversity of all living forms at different levels of complexity: genes, species, ecosystems and even landscapes and seascapes. Biodiversity is shaped by climatic conditions, the properties of soils and sediments, evolutionary processes and human action. Biodiversity can be greatly enhanced by human activities; however, it can also be adversely impacted by such activities due to unsustainable use or by more profound causes linked to our development models."
It is clear that climate change and biodiversity are interconnected. Biodiversity is impacted by climate change but it also makes an important contribution to both climate-change mitigation and adaptation through the ecosystem services that it supports. Therefore, conserving and sustainable managing biodiversity is crucial to meet the clime change.
Biodiversity conservation is an urgent environmental issue that must be specially attended. It is as critical to humans as it is to the other lifeforms on Earth. Countries of the world acknowledge that species research is crucial in order to obtain and develop the right methods and tools to understand and protect biodiversity. Thus, biodiversity conservation has became a top priority for researchers 
In this sense, a big effort is being carried out by the scientific community in order to study the huge biodiversity present on the planet. Sadly, spiders have been one of most unattended groups in conservation biology . These arachnids are plentiful and ecologically crucial in almost every terrestrial and semi-terrestrial habitats   . Moreover, they present a series of extraordinary qualities, such as the ability to react to environmental changes and anthropogenic impacts  .
Several works have studied the spiders’ behavior. Some of them focuses on the use of the way spiders build their webs as a source of information for species identification . Artificial intelligent systems have been proven to be of incalculable value for these systems.  proposed a model for spider behavior modeling, which provides simulations of how a specific spider specie builds its web.  recorded how spiders build their webs in a controlled scenario for further spatiotemporal analysis.
Because spider webs carry an incredibly amount of information, this chapter presents an study about its usage for automatic classification of spider species. In particular, computer vision and artificial intelligence techniques will be used for this aim. Moreover, the amount of information will be such, that it would be enough to perform the spider specie classification. This is, to authors extend, a novel approach on this problem.
The remainder of this paper is organized as follow. First, the database is briefly presented. Section 3 explains how images were preprocessed in order to extract the spiderwebs from the background. The feature extraction and classification techniques are introduced in sections 4 and 5. Next, experiments and results are shown in detail. Finally, the conclusions derived from the results are presented.
The database contains spider web images of four different species named Allocyclosa, Anapisona Simoni, Micrathena Duodecimspinosa and Zosis Geniculata. Each class has respectively 28, 41, 39 and 49 images, which makes a total of 150 images. Some examples can be seen in Figure 1. Since the webs imageswere taken in both controlled and uncontrolled environments the lightness condition and background differ between classes.
The images that correspond to
|Class||Number of samples||Size (pixels)||Bits number (color/gray)|
|1||28||1024x768||24 (True color)|
|2||41||2240x1488||8 (Gray scale)|
|3||39||2240x1488||24 (True color)|
|4||42||2216x2112||24 (True color)|
Class 2 (C2) corresponds to the
Class 3 (C3) is composed by
Finally, Class 4 (C4) corresponds to
The technique features of images from each class are summarized in Table 1
As can be seen in Figure 1, spider web images were taken in both controlled and uncontrolled environments. Thus, the preprocessing step was vital in order to isolate the spider webs and remove possible effects of background in the system.
Image processing techniques were employed in order to isolate the spider webs from light reflections and elements of the background of the image such as leaves or tree branches. Once it was applied, a new normalized database was obtained.
3.1. Spider web selection
Since the image collection was not taken for this research, there are information that does not provide valid data for the spiders’ study. This information corresponds to any external element of the spiderwebs. This is why the region of interest of the image was manually selected. This is represented in figure 2.
Once the spider web has been select, an adjustment of the proportional ratio was necessary in order to obtain a proportional square image. This will be explained in detail in the following section.
3.2. Image contrast
To enhance the contour of cobweb’s threads an increase of the color contrast was first applied. A spacial filtering was applied to enhance or attenuate details in order to obtain a better visual interpretation and prepare the data for the next preprocessing step. By using this filtering, the value of each pixel is modified according to neighbors’ values, transforming the original gray levels so that they become more similar or different to the corresponding neighboring pixels.
In general, the convolution of a image
3.3. Image binarization
The binarization process transforms the image to a black and white format in a way that it does not change the esential properties of the image. Equation 2 defines the binarization process, where
This threshold is computed by using the very know Otsu’s method , which assigns the membership of each pixel to a determined group by computing the optimal value from which carrying out that assignment.
3.4. Image denoising
Once the image has been binarized, a denoising process was used aiming to eliminate any irrelevant information. To achieve this goal, two specific techniques were applied:
The Wiener Filtering applied a spatial filtering using statistical methods in order to reduce noise and smooth shapes. It gradually smooths the image by changing the areas where the noise is very apparent, but keeping the areas where the details are present and the noise is less apparent. The Wiener filter is adapted to the local image variance.
In this work, the algorithm
If the noise variance is not given,
On the other hand, morphological operations are those transformations that modify the structure or shape of the objects in the image based on the their geometry and shape, simplifying the images. These techniques can be used to denoise an image, for feature extraction or processing specific regions.
The resulting image after image denoise can be observed in figure 6.
3.5. Center of the spiderwebs
Finally, the center of the spiderwebs was used as the source of discriminative information in order to classify spiders. Thus, once the images were preprocessed this center area was selected to conform the experimentation database. Figures 7 and 8 show the center of the spiderwebs for each specie.
4. Features extractors
In general, feature extraction refers to the process of obtaining some numerical measures of images such as area, radio, perimeter, etc. Also it concerns to the process of transforming a set of original features; with dimension
Two well known techniques were used: Discrete Wavelet Transform (DWT) and Discrete Cosine Transform (DCT). These were selected as they have been successfully used in other biometric studies. Besides being able to reduce the dimensionality of data and reduce the computational requirements of the classifier stage, these techniques can improve the generalization of the information and the system’s success rate.
4.1. Wavelet transform
The discrete wavelet transform (DWT) is based on the idea of decomposing a signal in terms of displaced and dilated versions of a finite wave called mother wavelet. The Wavelet transform is a preprocessing and feature extraction technique which can be directly applied to the image of spiderwebs. The DWT is defined in  as follows:
where is the transform function:
In the wavelet analysis is common evaluate the results as approximations and details. The approximations are the low frequency components of the signal and the details are the high frequency components. For many signals the most important information is content in the low frequencies. This content is what gives identity to the signal. The following is the diagram that has been used in this work, which would be one-dimensional DWT. The filtering process to obtain the approximations and detail of the discrete wavelet transform shown in the following figure 9:
The application of different mother families on pre-processing (artifacts elimination) and on the feature extraction has a set of good and discriminate parameters. Unlike the Fourier transform, the wavelet transform can be implemented on many bases. The different categories of wavelets (continuous, discrete, orthogonal, etc..) and various types of wavelet functions within each category provide a wide number of options for analyzing a signal. This allows selection of the base functions whose shape better approximates the characteristics of the signal to be represented or analyzed. On this work, the families Daubechies1 (db1), Biorthogonal 3.7 (bior3.7) and Discrete Meyer (dmey) were used.
4.2. Discrete Cosine Transform
Discrete Cosine Transform (DCT) was applied for noise and details of high frequency elimination . Besides, this transform has a good energy compaction property that produces uncorrelated coefficients, where the base vectors of the DCT depend only on the order of the transformation selected, and not of the statistical properties of the input data.
Another important aspect of the DCT is its capacity to quantify the coefficients utilizing quantification values, which are chosen of visual way. This transformation has had a great acceptance inside the image digital processing, as there is a high correlation among elements for the data of a conventional image.
5. Classification: Support Vector Machine
Once the images were transformed to a set of features, the classification stage tried to produce an answer to the spider identification problem. In this work, the well known Support Vector Machine (SVM) technique has been used.
The SVM is a method of structural risk minimization (SRM) derived from the statistical learning theory developed by Vapnik and Chervonenkis . It is enclosed in the group of supervised learning methods of pattern recognition, and it is used for classification and regression analysis.
Based on characteristic points called Support Vectors (SVs), the SVM uses an hyperplane or a set of hyperplanes to divide the space in zones enclosing a common class. Labeling these zones the system is able to identify the membership of a testing sample. The interesting aspect of the SVM is that it is able to do so even when the problem is not linearly separable. This is achieved by projecting the problem into a higher dimensional space where the classes are linearly separable. The projection is performed by an operator known as kernel, and this technique is called the kernel trick  . The use of hyperplanes to divide the space gives rise to margins as shown in figure 10.
In this work, the Suykens’ et. al. LS-SVM  was used along with the Radial Basis Function kernel (RBF-kernel). The regularization parameter and the bandwidth of the RBF function were automatically optimized by the validation results obtained from 10 iterations of a Hold-Out cross-validation process. Two samples from each class (from the training set) were used for testing and the remaining for training as we saw that the number of training samples has a big impact in the LS-SVM optimal parameters. Once the optimal parameters were found, they were used to retrain the LS-SVM using all available training samples.
6. Experiments and results
To sum up, the proposed system normalized all the images to 10x10 pixels. This system used the first M features obtained from the DCT projection of the spider webs images and the outcome of the DWT transformation of the spider webs images as inputs for a RBF-kernel LS-SVM with regularization and kernel parameters. The former parameter (the number of features) was varied during experimentation, while the later two parameters (the regularization and the kernel parameters) were automatically optimized by iteration using validation results. To obtain more reliable results, the available samples were divided into training and test sets, so that the system was trained and tested with totally different samples.
The well known K-Folds cross-validation and Hold-Out cross-validation techniques were used to obtain the final results. In particular, experiments with K equal 3, 5, 7, and 10 were run. The percent of training samples in Hold-Out cross-validation was 50, 40, 30, 20, 10 respectively. It is worth it to mention that the training and testing sets were computed for each class individually, having into account that each class has different number of samples. These experiments were performed for both datasets, i.e. using the whole spiderwebs and only the center area.
6.1. DCT results
In order to obtain the optimal number of coefficients, 30 experiments were performed using the Hold-Out cross validation technique. As the size of image was normalized to 10x10 pixels the total number of characteristics corresponded to 100, therefore, in this phase, the number of coefficients was swept from 1 to 100 coefficients. Figure 11 represents the mean of those 30 experiments. It can be observed that 60 is the optimal number of coefficients.
Table 2 shows the results reached for K-Fold cross-validation and Hold-Out cross-validation using the optimal number of coefficients.
|Success rate (%)||Hold-Out cross-validation Success rate (%) (% of training samples)|
|10||98.75% ± 0.18||50||98.69% ± 1.20|
|7||98.56% ± 0.58||40||98.45% ± 1.38|
|5||98.44% ± 0.76||30||96.89% ± 2.06|
|3||98.07% ± 0.52||20||94.29% ± 4.07|
|-||-||10||79.75% ± 6.54|
6.2. DWT results
In this case, the length of the feature vector depends on the wavelet family. Windows
The Table 3, 4, 5 shows the results reached for K-Fold cross-validation and Hold-Out cross-validation for each DWT family.
|Success rate (%)||Hold-Out cross-validation|
(% of training samples)
|Success rate (%)|
|10||99.40% ||50||98.42% |
|7||99.31% ||40||96.89% |
|5||99.18% ||30||96.26% |
|3||98.47% ||20||94.38% |
|Success rate (%)||Hold-out cross-validation|
(% of training samples)
|Success rate (%)|
|10||99.40% ± 0.53||50||98.42% ± 1.74|
|7||99.31% ± 0.27||40||96.89% ± 1.93|
|5||99.18% ± 0.31||30||96.26% ± 2.40|
|3||98.47% ± 0.12||20||94.38% ± 2.89|
|-||-||10||91.08% ± 4.28|
|Success rate (%)||Hold-Out cross-validation|
(% of training samples)
|Success rate (%)|
|10||98.62% ± 0.36||50||97.70% ± 2.10|
|7||98.30% ± 0.50||40||97.12% ± 1.54|
|5||98.55% ± 0.25||30||95.66% ± 2.29|
|3||97.89% ± 0.35||20||94.77% ± 2.28|
|-||-||10||89.14% ± 6.28|
7. Discussion and conclusions
This work has faced the problem of spider web recognition improving the results obtained by the previous work . It is important to note that, to the authors extend, these are the only published works using the proposed technique.
Images were preprocessed to isolate the center of the spider web and remove the effects of the background in the system. The resulting images were then transformed by using DCT and DWT. For the former, the optimal number of DCT coefficients
The results confirmed the improvement compared to , where only three species (versus the four species used in thius work) were classified with a maximum success rate of 95%. Thus, tables 2, 3, 4 and 5 show that the new system reached performance of around 99% on K-Fold cross-validation and 98% on Hold-Out cross- validation. Moreover, the obtained standard deviation was significantly low, although, as expected, slightly higher on Hold-Out as the number of samples for training lowered. All in all, the standard deviation achieved in both K-Fold and Hold-Out procedures are smaller than those obtained on .
When comparing DCT and DWT, the DWT provided a better behavior for this problem. It is worth it to emphasize that the images have been normalized to a size of 10x10, this is, quite compressed, considering the spatial distribution of the threads in the spider webs.
The results achieved by this work support the conclusions derived from  stating that the center of the spiderwebs provide enough discriminative information to recognize different species of spiders. However, it is still necessary to run more experiments with a larger database and execute a more detailed study on which parts of the spiderweb provide the most discriminative information before make stronger conclusions. On the other hand, this will allow to test the system’s performance with larger training sets, which will be interesting having into account that the system clearly improved when the number of training samples increased.
Sytnik, K.M., Preservation of biological diversity: Top-priority tasks of society and state 2010Ukrainian Journal of Physical Optics, 11 (SUPPL. 1), S2 S10
- 2. Carvalho, J. C., Cardoso, P., Crespo, L.C., Henriques, S., Carvalho,R., Gomes, P., "Biogeographic patterns of spiders in coastal dunes along a gradient of mediterraneity." Biodiversity and conservation (2011):1-22.
Johnston J. M. 2000 The contribution of microarthropods to aboveground food webs: A review and model of belowground transfer in a coniferous forest 143 226 238
Peterson A. T. Osborne D. R. Taylor D. H. 1989Tree trunk arthropod faunas as food resources for birds. Ohio Journal of Science 89 1 23 25
Cardoso P. MA Arnedo Triantis. K. A. Borges P. A. V. 2010 Drivers of diversity in Macaronesian spiders and the role of species extinctionsJ Biogeogr 37: 1034 EOF 1046 EOF
Finch-D O. Blick T. Schuldt A. 2008 Macroecological patterns of spider species richness across EuropeBiodivers Conserv 17: 2849 EOF 2868 EOF
Eberhard,W.G., Behavioral Characters for the Higher Classification of Orb-Weaving Spiders, Evolution, Vol. 5Sep., 1982 1067 1095Society for the Study of Evolution
Eberhard,W.G., Early Stages of Orb Construction by Philoponella Vicina, Leucauge Mariana, and Nephila Clavipes (Araneae, Uloboridae and Tetragnathidae), and Their Phylogenetic Implications, Journal of Arachnology, Vol. 2Summer, 1990 205 234American Arachnological Society
Eberhard W. G. Computer Simulation. of-Web Orb. Construction J. American Zoologist. pp 229 238February 1, 1969
th European colloquium of arachnology, Aarhus, Denmark, Suresh P. B. Zschokke S. computerised A. method to. observe spider. web building. behaviour in. a. semi-natural light. environment 2000
IWINAC 2011, Part II, LNCS 6687, Ticay-Rivas Jaime. R. del Pozo-Baños Marcos. Eberhard William. G. Alonso Jesús. B. Travieso Carlos. Spider Recognition. by Biometric. Web Analysis. I. W. I. N. A. 409 417 2011
- 12. Jing Hu; Si, J.; Olson, B.P.; Jiping He;, "Feature detection in motor cortical spikes by principal component analysis," Neural Systems and Rehabilitation Engineering, IEEE Transactions on, vol.13, no.3, pp.256-262, Sept. 2005.
- 13. Qingfu Zhang; YiuWing Leung;, "Aclass of learning algorithms for principal component analysis and minor component analysis," Neural Networks, IEEE Transactions on, vol.11, no.1, pp.200-204, Jan 2000.
- 14. Langley, P.; Bowers, E.J.; Murray, A.;, "Principal Component Analysis as a Tool for Analyzing Beat-to-Beat Changes in ECG Features: Application to ECG-Derived Respiration," Biomedical Engineering, IEEE Transactions on, vol.57, no.4, pp.821-829, April 2010.
Haibo Yao; Lei Tian;, "A genetic-algorithm-based selective principal component analysis (GA-SPCA) method for high-dimensional data feature extraction," Geoscience and Remote Sensing, IEEE Transactions on, vol.41, no.6, pp. 1469 1478June 2003
Nan Liu; Han Wang;, "Feature Extraction with Genetic Algorithms Based Nonlinear Principal Component Analysis for Face Recognition," Pattern Recognition, 2006ICPR 2006. 18th International Conference on, 3no., 461 464
- 17. V. Vapnik, “The Nature of Statistical learning Theory.” Springer Verlag, New York, 1995.
Vojislav Kevman. Learning and Soft Computing: Support Vector Machines, Neural Networks, and Fuzzy Logic models”Puiblished by The MIT Press, 2001
- 19. B. Schölkopf y A.J. Smola. “Learning with Kernels. Support Vector Machines, Regularization, Optimization, and Beyond“, Published by The MIT Press, 2002.
Suykens J. A. K. Van Gestel T. De Brabanter J. De Moor B. Vandewalle J. Least “. Squares Support. Vector Machines”. World Scientific. Singapore 2002 9-81238-151-1
Lim, Jae S., Two-Dimensional Signal and Image Processing, Englewood Cliffs, NJ, Prentice Hall, 1990 548
http://www.unesco.org/new/en/natural-sciences/special-themes/biodiversity-initiative/.Las visit in March 2011
Ahmed N. Natarajan T. Rao K. R. “. Discrete Cosine. Transform” I. E. E. E. transactions on. Computes 90 93 1974
Mallat, S., “A theory for multiresolution signal decomposition: the wavelet representation”, IEEE Pattern Analysis and Machine Intelligence, Vol. 7 674 693 1989
Otsu, N.; “A threshold selection method from gray-level histograms”.IEEE Trans. Sys., Man., Cyber. 9 1 62 66 1979