Content Based Image Retrieval Using Embedded Neural Networks with Bandletized Regions

One of the major requirements of content based image retrieval (CBIR) systems is to ensure meaningful image retrieval against query images. The performance of these systems is severely degraded by the inclusion of image content which does not contain the objects of interest in an image during the image representation phase. Segmentation of the images is considered as a solution but there is no technique that can guarantee the object extraction in a robust way. Another limitation of the segmentation is that most of the image segmentation techniques are slow and their results are not reliable. To overcome these problems, a bandelet transform based image representation technique is presented in this paper, which reliably returns the information about the major objects found in an image. For image retrieval purposes, artificial neural networks (ANN) are applied and the performance of the system and achievement is evaluated on three standard data sets used in the domain of CBIR.


Introduction
The number of digital images in the form of personalized and enterprise collections have grown immensely.Hence, there is a growing demand for powerful image indexing and retrieval in an automatic way.However, with such widespread use and availability of images, the textual annotation of images (using keywords) is becoming impracticable and unsuitable for the representation of images and their retrieval.This is the reason that content based image retrieval (CBIR) has become a great research interest amongst research communities [1][2][3].
Content based image retrieval systems generate meaningful image representations by considering the visual characteristics of images, i.e., color [4][5][6], texture [7][8][9][10], shape [11,12] and salient points [13,14] and bring closely resembling images in terms of distance as the semantic response.One of the major challenges faced in this regard is semantic gap, i.e., features at low level are not sufficient to characterize the high level image semantics [15].To bridge this gap to some extent, an important focus of research is on the enhancement of these features, so that the machine learning algorithms can make significant improvements to bridge this gap.Features can be split into two main categories, i.e., global image features and local image features of the region representations.Global image features are generated by considering the whole image as a single entity and then feature representations are obtained.Color histograms [6,16], texture [10,17,18], color layout [4,5] etc. are all examples of the global image features.Region representation or the local features are obtained either by the segmented regions [11,19,20] or by detecting the interest points, e.g., features obtained by scale-invariant feature transform (SIFT) [21] and sped up robust features (SURF) [14].Segmentation allows to generate the representations that take into account the actual objects of interest and avoid image regions that are less informative.This is the reason that the image representations obtained by the image segmentation are much more powerful than the global image representations.However, the drawback of the segmentation based image features is that currently there is no technique available that can perform the image segmentation in a better way, also the associated cost for image segmentation in terms of computational time makes it of no practical use [11,20,22,23].
To reap the benefits of the segmentation based image representations and to overcome the associated drawbacks, the focus of the current manuscript is on the identification of the image segments that contain major image objects by applying bandlet transform.Bandlet transform returns the geometric representation of the texture of the object regions, which can be used to discriminate the objects of interest in a fast way.The detailed procedure will be described in the Section 3.1.The major problem with the geometric output is that its empathy is complicated due to the close resemblance of the connected regions.In order to ensure the actual association, artificial neural networks are applied and correct texture classification is performed.We then apply the Gabor filter and generate the texture representation from it based on the classification output.To further enhance the image representation capabilities, we have also estimated the color content in the YCbCr color domain and defused it with texture.
As described by the Irtaza et al., major drawbacks faced by a CBIR system or query by image content which severely impact the retrieval performance [24] are: (1) the lack of output verification and (2) neighborhood similarity avoidance for the semantic association purpose.Therefore, we have followed their findings and also included the neighborhood in the semantic association process.Content based image retrieval is then performed by the artificial neural networks after training them with these obtained features.Bandelet transform is used for medical image retrieval [25].However, their approach for feature extraction is different than our approach.Before this, researchers have utilized the bandeletization property for image compression [26], image enhancement [27] and gender classification [28].
The remainder of paper is presented as follows: In Section 2 related work in the domain of CBIR are introduced.In Section 3 the proposed techniques is described in detail.Experimentation results are provided in Section 4. Finally, we have concluded our findings in Section 5.

Related Work
Numerous CBIR systems [1,29] are proposed so far to focus on the image searching problem and efficient image retrieval system in a more effective way.For this, the research intent is on the exploration of the new signature types [18] and powerful image similarity detection measures.In CBIR, image signature plays an imperative role to fabricate an efficient image search.Query image and images found in the repository are qualified as a collection of feature vectors and ranking of the relevant results occur on the basis of common norms, i.e., distance or semantic association by a machine learning technique [15].Signature development is usually performed through the analysis of color [4,16], texture [17], or shape [19] or by generating any of these combinations and representing them mathematically [8].Color features are extensively used in CBIR, which may be endorsed to the better potentiality in three dimensional domains over the gray level images which is single dimensional domain.Texture features as powerful visual features are used to capture repetitive patterns of a surface in the images.The formation of human identity and recognize objects in the real world is known as an important cue of the shape.Shape features form have been used for the purpose of retrieving images in many applications [30].Shape features extraction techniques are classified into contour based and region based methods.The classification of contour method based on shape boundary which is consist of contours and extracted features from the contour, while region based methods extract the features from the entire region.
A color estimation method in images through color difference histograms (CDH) is presented [16].Their work has diminished the reliance on the frequency of pixels in an image.The unique characteristic of CDH is that they count the perceptually uniform color differences between the two points with respect to color and edge orientations in L × a × b color space.Dominant colors and edge orientations are then used for image representation purposes.In the work of [31] a color-texture and dominant color based image retrieval system (CTDCIRS) is proposed and three different features from the images, i.e., motif co-occurrence matrix (MCM), dynamic dominant color (DDC), and difference between pixels of scan pattern (DBPSP) are offered.In their first step, using a color quantization algorithm, the image is divided into an eight coarse partition and the dominant colors are obtained from every partition.In their second step, MCM and DBPSP are applied for texture representation of an image using motif transform.The above main three features are amalgamated to expedite the efficiency of the system.In the work of Yue et al. [32], they introduced a method for feature extraction which is based on texture-color features.Combination of texture and color features are used for automatic retrieval purpose.In the initial step hue-saturation-value (HSV) color space is measured and co-occurrence matrices are then exerted to build the texture features.
The work of [33] presented a feature extraction approach by generating the curvelet representation of the images.The Curvelet transformations are combined with a vector codebook of region based sub-band clustering for extraction of dominant color.In this approach, the user defined image and target images are compared by using the principle of the most similar highest priority (MSHP) and evaluated for the retrieval performance.In the work of [34], a color histogram based on the wavelet is introduced, which also considers the texture and color component of the images for image retrieval.In the work of [35] the Micro Structure Descriptor (MSD) was proposed, which depend on the orientation of the edges.It depends on the primary colors in the small structures with similar edge orientation, which are simulated according to the human visual processing estimation.Hejazi and Ho [17] presented an image retrieval approach, which is based on textural information of an image and their approach considers the directionality, orientation and regularity of the texture using nonlinear modified discrete Radon transform, which are further used for image retrieval purposes.The Vector Quantization (VQ) is applied in the work of [18] for feature extraction.According to their approach, an image is cleaved into pixel blocks of 2 × 2 size in RGB domain.These blocks are then used for training vectors development.For the initial training sets purpose, linde buzo gray (LBG) and kekre proportionate error (KPE) algorithms are applied.After training, the same approach is used for testing and decent results are produced.
In the work of [36], they transformed the RGB image in the opponent chromaticity space.In the chromaticity space, Zernike chromaticity distribution moments are used to capture the characteristic of color contents.In the contour domain, texture features are extracted using a scale-invariant and rotation-invariant image descriptor.The combination of texture-color features of the human visual system provides efficient and flexible approximation processing.In the work of [13], they proposed a CBIR system using color and texture features of an image sub region.This is to identify regions of interest (ROI) almost by segmenting the image in the fixed sections after the application of morphological dilation.The color of the ROIs are calculated from the quantized histogram HSV color space and Texture are computed from gray level co-occurrence matrix (GLCM).The Query and target Images are compared with same numbers of Region of interests.Lin et al. [37], proposed three image features for efficient automatically retrieval of images.To extract the texture feature, difference between pixels of scan pattern (DBPSP)are used while color co-occurrence matrix (CCM)are used to get the color features.The last image feature depend on color distribution, called color histogram for K-mean (CHKM).
Irtaza and Jaffar [24] presented a potential solution of retrieving semantically similar images from large image repositories against any query image.Genetic algorithm and support vector machines are used to minimize the existing gap between high-level and low level features.To avoid the risk of dissociation, relevance feedback is also incorporated in their work.In the work of [38], presented a gender classification technique which is based on the feature selection through Eigen values of Gabor filters and Haar based wavelet packets.As shown the experiment result, Gabor features combined with random mean feature values they improve the performance of image retrieval and classification accuracy.Jhanwar et al. [39], presented a new technique for CBIR using motif co-occurence matrix (MCM).MCM is obtained by using motif transformed image and is similar to color co-occurence matrix.
In the work of [40], texture is the main difficulty in the way of segmentation.It is very difficult to estimate the parameters of texture model when texture segmentation are applied.These problem are controlled by using J-image segmentation (JSEG).JSEG consists of two independent steps: spatial segmentation and color quantization.Wang et al. [22] introduced a semantics classification method, which uses a wavelet-based approach for feature extraction and then for the image comparison.They used image-segmentation-based region matching approach.integrated region matching (IRM) proposed in this work is not efficient for texture classification due to uncertain modeling.Therefore, to address this issue, their idea was further processed by Chen et al. [23].In the work of [23], they proposed a method cluster-based retrieval of images by unsupervised learning (CLUE).They proposed an unsupervised clustering-based technique, which generate multiple clusters of retrieved results and give more accurate results as compared to the previous work, but their method suffers from issues such as numbers of clusters identification and segmentation uncertainty, due to which the results of this technique are not reliable.ElAlami [41] proposed a model which is based on three different techniques: the first one is concerned with the features extraction from images database.For this purpose, Color histogram and Gabor filter are used to extract the color and texture features.While the second technique depends on genetic algorithm and obtained the optimal boundaries of these discrete values, in the last technique, the selection of features consist of two successive function which are called preliminary and intensely reduction for extracting the most similar features from the original feature repository sets.
Some other visual features are also proposed for CBIR, Important points such as salient points and spatial features.SIFT [21] and SURF [32] based on salient points found in an image are the familiar visual features.Researchers have done a lot of work by using these salient points in content based image retrieval.Velmurugan and Baboo [14] applied SURF features by combining them with the color features to improve the retrieval accuracy.Mallat and Peyré [42] introduced bandelet approaches to geometric image representations.Orthogonal Bandelets using an adaptive segmentation are well appropriate solution to capture the regularity of edge structure.They applied wavelet coefficient of orthogonal transformation based on bandeletization.
Qu et al. [43] proposed a system based on bandelet transform and represented sharp image transitions such as edges by taking the advantage of geometric regularity of the image structure in image fusion.To create the fused image max rule is applied to select the geometric flow of the source image and bandelet coefficient.
The technique that we present in the current paper considers the most prominent objects which exist in an image using the object geometric representation obtained by bandelet transform in a precise manner.The texture information found in object boundaries are then utilized for use as the component of the feature vectors after applying the targeted parameters to the Gabor transform based on the Artificial Neural Network suggestions.The features are further improved by incorporating the color information in YCbCr domain.Image semantics are then obtained by the Artificial Neural Networks.

Proposed Method
The most important capability of the proposed method is its attribute for identifying the most prominent objects in an image.These objects are then considered as the core outcomes which are used for the generation of feature vectors.To achieve this, first of all, image transformations are generated through bandlet transform, which return the geometric boundaries of the major objects found in an image.We apply Gabor filter with targeted parameters (as will be described) to estimate the texture content around these boundaries.These geometric boundaries are vague in a sense that they could easily be deceived to be associated with unwanted texture classes as all of them closely resemble to one another, and if not carefully considered can result in the form of wrong parameter estimation, which will be consequence of the form of unsatisfactory image retrieval output.Therefore, to avoid this situation, geometric classification is performed through the backpropagation neural networks, which make it certain that the texture estimation parameters to apply Gabor filter should be approximated with maximum accuracy.
To increase the power of feature vectors, color components are also induced in the YCbCr domain after approximating it through wavelet decomposition over the color histograms.The proposed features are applied on all images present in the image repository, and their semantic classes are determined through ground truth training with Artificial Neural Networks and the finer neighborhood of every image.We generated inverted index over the semantic sets, which guarantees the fast image retrieval after determining the semantic class of query image.The Complete process of the proposed method is represented in Figure 1 and the detail of the process will come in the subsequent sections.

Bandelet Transform
The issue with Wavelet bases was that the same values of texture have different directions in an image.To overcome this limitation Le Pennec and Mallat et al. [42,44] proposed the Geometric regularity in an anisotropic way by eliminating the redundancy of wavelet transform using the concept of bandeletization.Bandelet transform is a major self adaptive multiscale geometry analysis method which exploits the recognized geometric information of images as compared to the non adaptive algorithms such as curvelet [9,33] and Contourlet transforms [7].Bandelet transform not only has the uniqueness of multiscale analysis, directionality and anisotropy but it also presents particular possessions of severe sampling and adaptability for image representation.Bandelet basis rules accumulate carriers extend in the direction perpendicular to the regularity of the maximum of the function as shown in Figure 2. Alpert transform is used for bandeletization that closely follows the geometry of underlying images.The main objective is to take the advantage of sharp image transitions by computing the geometric flow to form bandelet bases which capture the grayscale images constantly changing direction.As shown in Figure 3, bandelet transform divides the image into square blocks and obtains one contour (Ω i ) from it.If a small image region does not contain any contour, it means that the image intensity is uniform and regular in that region therefore the flow of line is not defined.

Alpert bases in bandelet transform
As per the work of [45], Alpert transform was applied to compute the bandelet bases to approximate images having some geomatric regularity.For this, the image is divided into square blocks and the geometric flow is estimated in every block.These square block are represented by S. In our implementation the block size is 8 × 8.
As we elaborated in Figure 3, if we use smaller blocks, i.e., 4×4, then it will divide the image in more chunks.However, the drawback is that in this case, the bandlet transform would not be able to capture the sharp edges and similarly, if the block size is larger, i.e., 16 × 16 or 32 × 32, then the geometric flow exceeds the object boundaries.Hence, through the experimental observations we used the block size of 8 × 8 for appropriate object estimation.Alpert transform parallel to the geometric flow which is constructed over the space l 2 (S) of wavelet coefficients in S with piecewise polynomials over bands of dyadic widths.A geometric flow direction γ is assumed to be known over S and warping operator W warps S into derivative of S. In the theory of warping function, any point x n = 2 j n is warped into xn = W (2 j n).Similarly l 2 (S) represents the function sampled in warped domain, i.e., ḡ(x n )(2 j n) ∈ S. To explain the multiresolution for each scale 2 l , the warped square S is recursively subdivided into 2 −l horizontal bands through: In Equation (1), 2nd term depends upon βl,i = βl−1,2i ∪ βl−1,2i+1 .Now the value of band is calculated in original Square "S" using Alpert multiresolution space Vl ⊂ l 2 ( S) i.e. βl,i = def W −1 ( βl,i ) ∈ S that has the width roughly equal to λ2 l and sampling points equals to 2 l (λ2 l ) 2 .Alpert vector is obtained through Equations ( 2) and (3): According to multiresolution space orthogonal bases ( hl,i,k ) i,k of each space are obtained by Gram Schmidt orthogonalization and resulting vector is obtained as: Alpert wavelets ( Ψl,i,k ) i.k are the orthogonal bases of orthogonal complement ( wl ) of ( Ṽl ).Therefore we compute Alpert wavelet vectors ( Ψl,i,k ) k after applying Gram-Schmidt orthogonalization of the family.
The consequential multiwavelet vectors ( Ψl,i,k ) have vanishing moments over the warped domain which is orthogonal to The above equation satisfies the condition, where (x) k = (x 1 ) k1 for each point xn = (x 1 , x2 ) in the warped domain.The orthogonal bases ( Ψl,i,k ) l,i,k of l 2 ( SS) describes an orhogonal alpert bases in l 2 (S) domain.
In square block at scale 2 l we calculate the orthogonal alpert bases β(S, γ) of l 2 (S) by: Bandelet transform provides the information of each square block and the flow of each square S is undefined then the projection onto β(S, γ) leaves the wavelet coefficient in S as unchanged.Bandelets have to provide the dyadic segmentation with bandeletization bases β(Γ j ) of the whole space of wavelet coefficients at a scale 2 j .
After applying alpert transform, we get a vector for each square, i.e., In the equation above ψv [n] are the coordinates of bandelet function and fit in the space of L 2 ([0, 1] 2 ).These coordinates are further used to calculate the bandelet bases which is called bandelization.
Bandelet bases are important factor to calculate the geometric images.Therefore in bandelet transform best bases are obtained by minimizing the Lagrangian function.
Bandelet transform pursues the above equation to generate the geometries of an image.In the above equation, "T" is the threshold.The value of threshold has impacts on the diversity of the image estimation.Different values can be adopted, i.e., 32, 48, 56, ..., etc.Therefore, in our implementation, we used a threshold value of 70 after detailed experimentation, which is able to estimate the theme object in an image.In the work of [43], each block is approximated in separable wavelet bases of L 2 (Ω) domain i.e., where I(Ω) is index set which depends upon the geometry of the boundary of (Ω) and x 1 , x 2 denotes the location of pixel in the image.Above equation represents the customized wavelets at the boundary and geometry flow is calculated in region (Ω).These wavelet bases are replaced by bandelet orthonormal bases of L 2 (Ω).Then, In the above equation c(x) defines the line flow associated to fix translation parameter x 2 , (x 1 , x 2 +c(x 1 )) be in the right place to (Ω) and is the direction of geometric flow which is extended.Then, c(x) is obtained as: This flow is parallel and c (x) is calculated as an expansion over translated function dilated by a scale factor 2 l .Then, the flow at this scale is characterized by: The bandeletization of wavelet coefficient use alpert transform to define a set of bandelet coefficients and by using these coefficients combined as inner product of original image f with bandelets: Local geometric flow depends upon these coefficient and scales.Therefore, for each scale 2 j of the orientation k a different geometry is obtained.After bandeletization process, we have achieved multiscale low and high pass filtering structure similar to wavelet transform.The above Equations ( 12) and ( 17) are used to calculate the geometry of the images as shown in Figure 4.The regions having contours are further used for texture classification to compute the features using Artificial Neural Network structure.Object categorization on the base of Geometric flow obtained through Bandletization.

Texture Feature Extraction using Bandelets
Texture is significant module of human visual perception and many researchers have done a lot of work to determine effectively characterize it in images.In this paper, we have proposed a new method to figure out the most prominent texture areas in the image that constitutes the major image objects.In the proposed method, first of all image transformations are generated through bandelet transform, which gives back the geometric boundaries of the major objects found in an image.We apply a Gabor filter with targeted parameters to estimate the texture content around these boundaries.These geometric boundaries are indefinite in the sense that they could easily be duped to be associated with undesired texture classes as all of them closely resemble to one another, and if not carefully considered can result in the form of wrong parameter estimation, which will be in consequence of the unsatisfactory image retrieval output.Therefore, to avoid this situation, geometric classification is performed through the back propagation neural networks, which makes it certain that the texture estimation parameters to apply Gabor filter should be approximated with maximum accuracy.The following are the main steps of texture feature extraction: (1) Convert input RGB image (I) of size M × N into gray scale image.
(2) Apply bandelet transform to calculate the geometry of an image and obtain the directional edges.
(3) Artificial Neural Network is used to classify these blocks having directional edges after training on the sample edge set as described in Figure 5. Once the network is trained, every geometric shape obtained in the step 2 is classified for parameter estimation.These parameters will further be described in the Gabor filter section.(4) After parameter estimation, the blocks with geometric contents are passed to the Gabor filter to estimate the texture.( 5) Steps 1 to 4 is repeated for whole image repository.

Artificial Neural Network
Neural networks (NN) are renounced as powerful and dominant tools in the area of pattern recognition and are inspired by biological neurons found in the human brain.The Least mean square rule and the gradient search method are used to minimize the average difference between input and target values on neural network [46].
The backpropagation neural networks are applied to classify the texture on the base of geometry returned by the bandelet transform.In this regard, we classify the texture directions in either horizontal, vertical, right/left diagonal or no contour blocks, by training on a small set developed manually.For this, we have placed 14 block samples representing the mentioned geometric types in every category, as described in Figure 5.To generate these samples, we consider only the image geometry and suppress the original image part.Once the network is trained, we apply it to classify every block present in the image.The reason to perform this task with the help of ANN instead of kernel (Window based operations used in image processing) is that, the geometry is not fixed and has different variations for same category.In this situation, the performance of the kernel based operations is miserable.Therefore, the ANN is applied and it classified the texture with maximum accuracy.Figure 6 shows the structure of neural network.The neural networks structure is defined with one hidden layer having 20 neurons and four output units.The sigmoid function is used in hidden layer and output layer as transfer function, i.e., Details of the neural network structure is summarized in Table 1.
After training of the neural networks, all blocks of an image are tested against the neural network, and their texture type is determined as: where m represents the total number of output units in the neural network structure, ȳf m returns the association factor of a particular output unit.The texture type m ↓ of the particular class is based on the output the neural network with highest association factor.

Gabor Feature
Gabor filters are extensively used in the vicinity of computer vision and pattern recognition.Several conquering applications of Gabor wavelet filter include feature extraction, segmentation of texture, face recognition, identification of finger prints , edge detection, contours detection, directional image enhancement, and image hierarchical representation, compression, recognition.Gabor is a strong technique to reduce the noise and can easily reduce image redundancy and repetition [47].Gabor filters can be convolved on a small portion of an image or can be convolved on the full image.An image region is expressed by the different Gabor responses generated through different orientations, different frequencies and angles [24,48].For an image I(x, y) with size M×N, its discrete Gabor wavelet transform is given by convolution: where s and t are filter mask size and ψ * mn is the complex conjugate of ψ mn which is the self similar function created from rotation and dilation of following wavelets, i.e.: where λ is the modulation frequency.The generating function is used to obtain self similar Gabor wavelets.
where m and n specify the scale and orientation of wavelet,with m = 0, 1, ..., M − 1 and n = 0, 1, ..., N − 1.In the above equation, we calculate the term x, ỹ, i.e.: where a > 1 and θ = nπ/N .In Gabor filter σ is the standard deviation of the Gaussian function, λ is the wavelength of harmonic function, and θ is the orientation.In our implementation, blocks having bandelet based geometric response are passed to the Gabor filter and based on the neural network classification, we select the parameters for the application of the Gabor filter [48].
Energy computation is performed using following equation: where F v is the feature vector.λ E are the eigen values, X is the eigen vector, and A is the Gabor response on a particular block.

Color Feature Extraction
In CBIR, color is the most imperative and significant visual attribute.It has been extensively studied and the motivation is that: the color estimation is not sensitive to rotation translation, and scale changes.Varieties of color spaces are available and serve effectively for different applications [6,49].In this paper, color features are extracted on the base of edge detection in YCbCr color space.Edges are extracted by applying canny edge detector on Y luminance component.The main steps of color feature extraction are as under:  6) To improve the feature performance, we apply wavelet transform at each histogram obtained in the previous step.We apply the discrete wavelet transform of H R at level 2, H G , and H B are applied at level 3.After this step, we have 128 bins i.e. 64 bins from H R , 32 bins from H G , and 32 bins from H B .(7) Calculate feature vector for every image in the repository.
Figure 7 shows that the color features are obtain as describe in above mentioned step.

Fusion Vector
Application of the above mentioned procedure generate two feature vectors, representing texture features obtained from Bandelet transform and color features obtained from YCbCr color space.Aggregation of these feature vectors in a single vector represents feature vector features against any image.

Content based image retrieval
Once the images present in the image repository are represented in the form of low level features, we can determine their semantic class.To determine the actual semantic class, a sub repository of the images is generated having the representation of "M" known classes, and every class contains R ≥ 2 images.In our implementation the value of "R" is set on 30, which means that 30 images from every semantic class from a ground truth image repository are used for the development of training repository.On this sub repository neural networks with specific class association parameters are trained.One against all classes (OAA) association rule is used for the network development, with the target to decrease the mean squared error between actual association and the NN based obtained association.The class specific training set can be defined as: Ω tr = Ω pos Ω neg where Ω pos are representing the R images from a particular class, and Ω neg are all other images in the training repository.Once the training is complete, semantic class of all images present in the image repository is determined on the base of decision function and the association rules.
Due to the object composition present in the images, many images may tend to show the association with some other classes, i.e., in the case of Mountain images which sometimes associate with the beach images.Therefore, a mechanism is required to reduce such associations.This is the reason, that the class finalization process also involves top "K" neighbors in the semantic association process using the majority voting rule (MVR) [50]: where C i (X) is the class wise association of input image and its top neighbors.
where l = {1, 2, 3..., n} represents the total number of neural network structures, and ȳf l returns the association factor of a particular neural network structure with a specific class.MVR counts the largest number of classifiers that agrees with each other [50].Therefore, the class association can be determined by: Once the semantic association of all images present in the image repository is determined, we store the semantic association values in a file that serves for us as the semantic association database.Therefore, after determining the semantic class of the query image through trained neural networks, we compute the Euclidean distance of the query image with all images of the same semantic class by taking into an account the values of the semantic association database, and generate the output on the basis of the feature similarities.

Performance Evaluation
To elaborate the retrieval capabilities of the proposed method numerous experiments are performed on three different image datasets.For implementation purposes, we have used Matlab 2010 in the Windows 7 environment using a core i3 machine by Dell.The detail of the experiments will be presented in the following subsections.Section 4.1 describes the datasets used for image retrieval purposes.Section 4.2 is about the retrieval precision and recall on randomly selected queries.Section 4.3 and 4.4 describe the comparison results of the proposed method with some state of the art works in CBIR.

Image Datasets
For our experiments we have used three image datasets namely: Corel, Coil, and Caltech 101.The Corel dataset contains 10,908 images of the size 384 × 256 or 256 × 384 each.For this dataset, we have reported the results on ten semantic categories having 100 images in each category.These semantic classes are namely: Africa, Buses, Beach, Dinosaurs, Buildings, Elephants, Horses, Mountains, Flowers, and Food.The reason for our choice to report the result on these categories is that: these categories are the same semantic groups used by most of the researchers who are working in the domain of CBIR to report the effectiveness of their work [31,33,37,39,41], so a clear performance comparison is possible in term of the reported results.To further elaborate the performance of the proposed system, experiments are also performed on Columbia object image library (COIL) [33] having 7200 images from 100 different categories.Finally, we have used the Caltech 101 image set.This dataset consists of 101 image categories and every category has different number of images in it.For the simplification purposes we have manually selected 30 categories which contain at least 100 images from every semantic category.

Retrieval Precision/Recall Evaluation
For our experiments we have written a computer simulation which randomly selects 300 images from image repository and uses them as the query image.As already described above, we are using image datasets having the images grouped in the form of semantic concepts, and on the base of their labels we can automatically determine their semantic association.We run this simulation on all three datasets mentioned previously and determine the performance by counting how many correct results are obtained against each query image.In proposed work, we reported the average result after running the five time experiments from each image category.For our experiments, a reverted index mechanism is proposed which after determining the semantic class of the query image, returns the relevant images against it, i.e., a method followed by the Google for text document search.According to the proposed method, we apply the trained neural networks on every image present in the image repository and determine its semantic class.The class association information is stored in a file which serves for us the semantic association database.The usability of this approach is that, after determining the semantic information for one time, we only need to determine the semantic class of the query image and relevance information about the semantic cluster is predetermined.Overall, the class association accuracy is determined in terms of the precision and recall using following formulas: where N A(q) represents the relevant images matching to the query image, N R(q) represents the retrieved images against the query image, and N t the total number of relevant images available in the database.Precision or specificity determines the ability of system to retrieve only those images which are relevant for any query image amongst all of the retrieved images, while Recall rate also known as sensitivity or true positive rate, determines the ability of classifier systems in terms of model association with their actual class.For the elaboration of results, top 20 retrieved images against any query image are used to compute the precision and recall.We have reported the average of the results as we mentioned previously after running our technique for five times.To elaborate the performance of the proposed system, we have randomly selected three images from 10 previously mentioned image categories from Corel image set and displayed their results in Figure 8.Here we have displayed the retrieval results which represent the precision obtained by our method against these query images.In this regard, we have displayed the results from the top 10 to top 40 retrievals against these randomly selected query images.The quantitative analysis of the proposed method suggests that the quality of the system is good in terms of precision as reliable results are appearing against these random selections.The most reliable results appear in the range of 10 to 30 images against query image as there are 100 images in a single category.The important thing to note here is that these are the results we have achieved without the involvement of any kind of the external supervision done by the user, as most of the relevance feedback based CBIR techniques do.In Figures 9 and 10, the same experiment is performed on Caltech 101 image set by randomly selecting 4 images.Precision and Recall are reported on top 10 to top 60 retrieval rates.Hence, on the base of retrieval accuracy we can say the proposed method is quite efficient.Another important thing to note is that the results reported here represent the retrieval against three random queries, while overall accuracy is reported on the average of 100 query images and performing the experiments for five times.

Comparison on Corel Image Set
To determine the usability of proposed method, it is compared with some state of the art methods in CBIR.In this regard, the technique is compared with [31,33,37,39,41].The reason for our choice to compare with these techniques is that these systems have reported their results on the common denomination of the ten semantic categories of Corel dataset as described earlier.Therefore, a clear performance comparison is possible.Table 2 explains the class wise comparison of the proposed system with other comparative systems in terms of precision.The results show that proposed system has performed better than all other systems in terms of average precision obtained.Table 3 describes the performance comparison in terms of Recall rates with the same systems.From the results, it could be easily observed that the proposed system has the highest recall rates.Figure 11 shows the performance of proposed method in term of precision with other state of art system.Figure 12 describes the performance comparison in terms of Recall rates with the same systems.From results it could be easily observed that the proposed system has the highest recall rates as show in Table 3.

Comparison on Coil Image Set
From the precision and recall results described for the Corel dataset, we can observe that integrating curvelet transform with enhanced dominant colors extraction and texture (ICTEDCT) has second highest rates in terms of precision and recall.Therefore we have reported the performance comparison on Coil dataset on different retrieval rates against ICTEDCT [33].For this experiment five images are selected from each image category and then performance of both systems is compared against each category.From the results elaborated in Figure 13, it can be clearly observed that proposed method is giving higher recall and precision rates as compare to ICTEDCT [33].Hence from the results of proposed method on Coil and Corel datasets, we can say that proposed method is much more precise and effective as compare to other CBIR systems.

Comparison with State-of-the-Art Methods
The retrieval results are compared with the state-of-art image retrieval methods, including the methods of Efficient content-based image retrieval using multiple support vector machines ensemble (EMSVM) [51], Simplicity [22], CLUE [23], patch based histogram of oriented gradients-local binary pattern ( Patch based HOG-LBP) [52], and Edge orientation difference histogram and color-SIFT (EODH and Color-SIFT) [53].The reason of our choice for comparison with these techniques is that: these systems have reported their results on the common denomination of the ten semantic categories of Corel dataset as described earlier.Hence, a clear performance comparison is possible with above state of art methods.Table 4 presents comparison of the proposed system with other comparative systems in terms of average precision.The results show that proposed system has performed better than all other systems in terms of average precision obtained.Same results are graphically illustrated in Figure 14.Table 5 describes the performance comparison in terms of Recall rates with the same systems.From the results, it can easily observed that the proposed system has the better recall rates.Figure 15 describes the performance comparison in terms of recall rates with the same systems.

Conclusions
With several application benefits, content based image retrieval has gained a lot of research attention.The paper has introduced a mechanism for automatic image retrieval.The major consideration of the paper was that the most prominent image results can appear if we generate the image representations that emphasize the core image objects instead of considering every image patch.Therefore, we have applied the bandelet transform for feature extraction, which considers the core objects found in an image.To further enhance the image representation capabilities, color features are also incorporated.Semantic association is performed through Artificial Neural Networks, and an inverted index mechanism is used to return the images against queries to assure the fast retrieval.The results of the proposed method are reported on three image data sets namely: Corel, Coil, and Caltech-101.The comparison results with other standard CBIR systems have revealed that the proposed system has outperforms other systems in terms of average precision and recall values.

Figure 2 .
Figure 2. Bandelet Transform [35,37].(a) Dyadic segmentation based on local directionality of the image; (b) A sample bandelet segmentation square that contains a strong regularity function shown by the red dash; (c) Geometric flow and sampling position; (d) Sampling position adapted to the warped geometric flow; (e) Illustration of a warping example.

Figure 4 .
Figure 4.Object categorization on the base of Geometric flow obtained through Bandletization.

Figure 6 .
Figure 6.The structure of neural network.

( 1 )
RGB image (I) is converted into Y C b C r color space.(2) After conversion, we separate the Y C b C r components and apply canny edge detector on Y component of the image.(3) In the next step, we combine the edges obtained in the previous step with unchanged C b and C r .(4) After step (3), convert the combined image back into single RGB image.(5) Now separate the individual R, G, and B components and calculate the histogram of each component.256 bins are obtained from H R , H G , and H B . (

Figure 8 .
Figure 8. Query Performance on Corel image dataset with top 10 to top 40 Retrievals.

Figure 9 .
Figure 9. Query Performance on Caltech image dataset with top 10 to top 60 Retrievals in terms of Precision.

Figure 10 .
Figure 10.Query Performance on Caltech image dataset with top 10 to top 60 Retrievals in terms of Recall.

Figure 12 .
Figure 12.Comparison of mean recall obtained by proposed method with other standard retrieval systems.

Figure 13 .
Figure 13.Comparison of precision and recall obtained by proposed method with integrating curvelet transform with enhanced dominant colors extraction and texture (ICTEDCT).

Figure 15 .
Figure 15.Comparison of mean recall obtained by proposed method with state of art retrieval systems.

Table 1 .
Summary of Neural network structure for every image category used in this work.

Table 2 .
Comparison of mean precision obtained by proposed method with other standard retrieval systems on top 20 retrievals.
Figure 11.Comparison of mean precision obtained by proposed method with other standard retrieval systems.

Table 3 .
Comparison of mean recall obtained by proposed method with other standard retrieval systems on top 20 retrievals.