learning overcomplete representations in an overcomplete basis the number of basis vectors is greater than the dimensionality of the input and the representation of an input is not a unique combination of basis vectors overcomplete representations have been advocated because they have greater robustness in the presence of noise can be sparser and can have greater flexibility in matching structure in the data overcomplete codes have also been proposed as a model of some of the response properties of neurons in primary visual cortex previous work has focused on finding the best representation of a signal using a fixed overcomplete basis or dictionary we present an algorithm for learning an overcomplete basis by viewing it as probabilistic model of the observed data we show that overcomplete bases can yield a better approximation of the underlying statistical distribution of the data and can thus lead to greater coding efficiency this can be viewed as a generalization of the technique of independent component analysis and provides a method for bayesian reconstruction of signals in the presence of noise and for blind source separation when there are more sources than mixtures

extrapolative spatial models for detecting perceptual boundaries in colour images we present a novel approach for automatically building image specific extrapolative spatial models of non boundary local energy that can be used to perform local statistical tests to detect perceptual boundaries the non boundary model consists of statistics of local energy that are spatially extrapolated by a non boundary confidence map and a scale adaptive normalised filtering algorithm we exploit the flexibility of steerable filters to both extract oriented local energy and to provide local statistics of the energy distribution in the orientation domain to compute the non boundary confidence map finally we apply our local thresholding technique separately to the three channels of colour images and adopt a max operator to combine the results we provide a qualitative and quantitative comparison on real images from a hand segmented natural image database against the best combination of the most widely cited colour edge detectors and automatic global thresholding methods

design of multiparameter steerable functions using cascade basis reduction an efficient method of computing the least squares optimal basis functions to steer any function locally is presented the method combines the lie group theoretic and the singular value decomposition approaches its efficiency is demonstrated with the design of basis functions to steer a gabor function under the four parameter linear transformation group

contour detection and hierarchical image segmentation this paper investigates two fundamental problems in computer vision contour detection and image segmentation we present state of the art algorithms for both of these tasks our contour detector combines multiple local cues into a globalization framework based on spectral clustering our segmentation algorithm consists of generic machinery for transforming the output of any contour detector into a hierarchical region tree in this manner we reduce the problem of image segmentation to that of contour detection extensive experimental evaluation demonstrates that both our contour detection and segmentation methods significantly outperform competing algorithms the automatically generated hierarchical segmentations can be interactively refined by user specified annotations computation at multiple image resolutions provides a means of coupling our system to recognition applications

sorted random projections for robust rotation invariant texture classification this paper presents a simple novel yet very powerful approach for robust rotation invariant texture classification based on random projection the proposed sorted random projection maintains the strengths of random projection in being computationally efficient and low dimensional with the addition of a straightforward sorting step to introduce rotation invariance at the feature extraction stage a small set of random measurements is extracted from sorted pixels or sorted pixel differences in local image patches the rotation invariant random features are embedded into a bag of words model to perform texture classification allowing us to achieve global rotation invariance the proposed unconventional and novel random features are very robust yet by leveraging the sparse nature of texture images our approach outperforms traditional feature extraction methods which involve careful design and complex steps we report extensive experiments comparing the proposed method to six state of the art methods rp patch lbp wmfs and the methods of lazebnik et al and zhang et al in texture classification on five databases curet brodatz uiuc umd and kth tips our approach leads to significant improvements in classification accuracy producing consistently good results on each database including what we believe to be the best reported results for brodatz umd and kth tips

analyzing weighted ell_1 minimization for sparse recovery with nonuniform sparse models in this paper we introduce a nonuniform sparsity model and analyze the performance of an optimized weighted l1 minimization over that sparsity model in particular we focus on a model where the entries of the unknown vector fall into two sets with entries of each set having a specific probability of being nonzero we propose a weighted l1 minimization recovery algorithm and analyze its performance using a grassmann angle approach we compute explicitly the relationship between the system parameters the weights the number of measurements the size of the two sets the probabilities of being nonzero so that when i i d random gaussian measurement matrices are used the weighted l1 minimization recovers a randomly selected signal drawn from the considered sparsity model with overwhelming probability as the problem dimension increases this allows us to compute the optimal weights we demonstrate through rigorous analysis and simulations that for the case when the support of the signal can be divided into two different subclasses with unequal sparsity fractions the weighted l1 minimization outperforms the regular l1 minimization substantially we also generalize our results to signal vectors with an arbitrary number of subclasses for sparsity

feature vector field and feature matching in this paper we propose a feature vector field for images which is built by the inner products and exterior products of image gradients the feature vector field effectively represents image edges and feature points including corners and edge points with big curvature using the feature vector field some novel descriptors are constructed for point matching and curve matching the descriptors are all invariant to image euclidean transformations and linear intensity changes the experiments show that the descriptors also have a good adaptability to small image affine transformation jpeg compression and nonlinear changes of intensity

certifying the restricted isometry property is hard this paper is concerned with an important matrix condition in compressed sensing known as the restricted isometry property rip we demonstrate that testing whether a matrix satisfies rip is np hard as a consequence of our result it is impossible to efficiently test for rip provided p np

image denoising using scale mixtures of gaussians in the wavelet domain we describe a method for removing noise from digital images based on a statistical model of the coefficients of an overcomplete multiscale oriented basis neighborhoods of coefficients at adjacent positions and scales are modeled as the product of two independent random variables a gaussian vector and a hidden positive scalar multiplier the latter modulates the local variance of the coefficients in the neighborhood and is thus able to account for the empirically observed correlation between the coefficient amplitudes under this model the bayesian least squares estimate of each coefficient reduces to a weighted average of the local linear estimates over all possible values of the hidden multiplier variable we demonstrate through simulations with images contaminated by additive white gaussian noise that the performance of this method substantially surpasses that of previously published methods both visually and in terms of mean squared error

optical flow using spatiotemporal filters a model is presented consonant with current views regarding the neurophysiology and psychophysics of motion perception that combines the outputs of a set of spatiotemporal motion energy filters to estimate image velocity a parallel implementation computes a distributed representation of image velocity a measure of image flow uncertainty is formulated preliminary results indicate that this uncertainty measure may be used to recognize ambiguity due to the aperture problem the model appears to deal with the aperture problem as well as the human visual system since it extracts the correct velocity for some patterns that have large differences in contrast at different spatial orientations

subband coding of images subband coding has become quite popular for the source encoding of speech this paper presents a simple yet efficient extension of this concept to the source coding of images we specify the constraints for a set of two dimensional quadrature mirror filters qmf s for a particular frequency domain partition and show that these constraints are satisfied by a separable combination of one dimensional qmf s bits are then optimally allocated among the subbands to minimize the mean squared error for dpcm coding of the subbands also an adaptive technique is developed to allocate the bits within each subband by means of a local variance mask optimum quantization is employed with quantizers matched to the laplacian distribution subband coded images are presented along with their signal to noise ratios snr s the snr performance of the subband coder is compared to that of the adaptive discrete cosine transform dct vector quantization and differential vector quantization for bit rates of 0 67 1 0 and 2 0 bits per pixel for 256 256 monochrome images the adaptive subband coder has the best snr performance

noise reduction procedures for gravity gradiometer data abstractnoise suppression of airborne gravity gradiometer data is a crucial part of the data processing stream we considered two approaches to removing noise kriging and directional filtering kriging is an estimation procedure for the interpolation of spatial data the estimator is calculated from the data variogram which characterizes the noise level and correlation length of the measurements directional filtering uses a user defined operator that is oriented to preferentially smooth the data along the strike but it leaves short wavelength components in the cross strike direction for definition of the trend edges both methods were applied to a recently collected offshore gravity gradient survey the kriging and directional filtering results revealed a similar level of smoothness but the main difference between them was the extra smoothing along the strike for the directionally filtered data because kriging is a data driven procedure it provides an objective estimate of the data noise level and

minimax estimation via wavelet shrinkage we attempt to recover an unknown function from noisy sampled data using orthonormal bases of compactly supported wavelets we develop a nonlinear method which works in the wavelet domain by simple nonlinear shrinkage of the empirical wavelet coefficients the shrinkage can be tuned to be nearly minimax over any member of a wide range of triebel and besov type smoothness constraints and asymptotically minimax over besov bodies with p q linear estimates cannot achieve even the minimax rates over triebel and besov classes with p 2 so the method can significantly outperform every linear method e g kernel smoothing spline sieve in a minimax sense variants of our method based on simple threshold nonlinear estimators are nearly minimax our method possesses the interpretation of spatial adaptivity it reconstructs using a kernel which may vary in shape and bandwidth from point to point depending on the data least favorable distributions for certain of the triebel and besov scales generate objects with sparse wavelet transforms many real objects have similarly sparse transforms which suggests that these minimax results are relevant for practical problems sequels to this paper which was first drafted in november 1990 discuss practical implementation spatial adaptation properties universal near minimaxity and applications to inverse problems

using basic image features for texture classification representing texture images statistically as histograms over a discrete vocabulary of local features has proven widely effective for texture classification tasks images are described locally by vectors of for example responses to some filter bank and a visual vocabulary is defined as a partition of this descriptor response space typically based on clustering in this paper we investigate the performance of an approach which represents textures as histograms over a visual vocabulary which is defined geometrically based on the basic image features of griffin and lillholm proc spie 6492 09 1 11 2007 rather than by clustering bifs provide a natural mathematical quantisation of a filter response space into qualitatively distinct types of local image structure we also extend our approach to deal with intra class variations in scale our algorithm is simple there is no need for a pre training step to learn a visual dictionary as in methods based on clustering and no tuning of parameters is required to deal with different datasets we have tested our implementation on three popular and challenging texture datasets and find that it produces consistently good classification results on each including what we believe to be the best reported for the kth tips and equal best reported for the uiuctex databases

on photometric issues in 3d visual recognition from asingle 2d image we describe the problem of recognition under changing illumination conditions and changing viewing positions from a computational and human vision perspective on the computational side we focus on the mathematical problems of creating an equivalence class for images of the same 3d object undergoing certain groups of transformations mostly those due to changing illumination and briefly discuss those due to changing viewing positions the computational treatment culminates in proposing a simple scheme for recognizing via alignment an image of a familiar object taken from a novel viewing position and a novel illumination condition on the human vision aspect the paper is motivated by empirical evidence inspired by mooney images of faces that suggest a relatively high level of visual processing is involved in compensating for photometric sources of variability and furthermore that certain limitations on the admissible representations of image information may exist the psychophysical observations and the computational results that follow agree in several important respects such as the same apparent limitations on image representations

error analysis of reweighted l 1 greedy algorithm for noisy reconstruction sparse solutions for an underdetermined system of linear equations x u can be found more accurately by l 1 minimization type algorithms such as the reweighted l 1 minimization and l 1 greedy algorithms than with analytical methods in particular in the presence of noisy data recently a generalized l 1 greedy algorithm was introduced and applied to signal and image recovery numerical experiments have demonstrated the convergence of the new algorithm and the superiority of the algorithm over the reweighted l 1 minimization and l 1 greedy algorithms although the convergence has not yet been proven theoretically in this paper we provide an error bound for the reweighted l 1 greedy algorithm a type of the generalized l 1 greedy algorithm in the noisy case and show its improvement over the reweighted l 1 minimization

steerable pyramids and tight wavelet frames in l_ 2 bbr d we present a functional framework for the design of tight steerable wavelet frames in any number of dimensions the 2 d version of the method can be viewed as a generalization of simoncelli s steerable pyramid that gives access to a larger palette of steerable wavelets via a suitable parametrization the backbone of our construction is a primal isotropic wavelet frame that provides the multiresolution decomposition of the signal the steerable wavelets are obtained by applying a one to many mapping n th order generalized riesz transform to the primal ones the shaping of the steerable wavelets is controlled by an m m unitary matrix where m is the number of wavelet channels that can be selected arbitrarily this allows for a much wider range of solutions than the traditional equiangular configuration steerable pyramid we give a complete functional description of these generalized wavelet transforms and derive their steering equations we describe some concrete examples of transforms including some built around a mallat type multiresolution analysis of l2 rd and provide a fast fourier transform based decomposition algorithm we also propose a principal component based method for signal adapted wavelet design finally we present some illustrative examples together with a comparison of the denoising performance of various brands of steerable transforms the results are in favor of an optimized wavelet design equalized principal component analysis which consistently performs best

the hermite transform theory the author introduces a scheme for the local processing of visual information called the hermite transform the problem is addressed from the point of view of image coding and therefore the scheme is presented as an analysis resynthesis system the objectives of the present work however are not restricted to coding the analysis part is designed so that it can also serve applications in the area of computer vision indeed derivatives of gaussians which have found widespread application in feature detection over the past few years play a central role in the hermite analysis it is also argued that the proposed processing scheme is in close agreement with current insight into the image processing that is carried out by the human visual system in particular it is demonstrated that the hermite transform is in better agreement with human visual modeling than gabor expansions

a theory for multiresolution signal decomposition the wavelet representation multiresolution representations are effective for analyzing the information content of images the properties of the operator which approximates a signal at a given resolution were studied it is shown that the difference of information between the approximation of a signal at the resolutions 2 sup j 1 and 2 sup j where j is an integer can be extracted by decomposing this signal on a wavelet orthonormal basis of l sup 2 r sup n the vector space of measurable square integrable n dimensional functions in l sup 2 r a wavelet orthonormal basis is a family of functions which is built by dilating and translating a unique function psi x this decomposition defines an orthogonal multiresolution representation called a wavelet representation it is computed with a pyramidal algorithm based on convolutions with quadrature mirror filters wavelet representation lies between the spatial and fourier domains for images the wavelet representation differentiates several spatial orientations the application of this representation to data compression in image coding texture discrimination and fractal analysis is discussed

nonlocal hierarchical dictionary learning using wavelets for image denoising exploiting the sparsity within representation models for images is critical for image denoising the best currently available denoising methods take advantage of the sparsity from image self similarity pre learned and fixed representations most of these methods however still have difficulties in tackling high noise levels or noise models other than gaussian in this paper the multiresolution structure and sparsity of wavelets are employed by nonlocal dictionary learning in each decomposition level of the wavelets experimental results show that our proposed method outperforms two state of the art image denoising algorithms on higher noise levels furthermore our approach is more adaptive to the less extensively researched uniform noise

image denoising via sparse and redundant representations over learned dictionaries we address the image denoising problem where zero mean white and homogeneous gaussian additive noise is to be removed from a given image the approach taken is based on sparse and redundant representations over trained dictionaries using the k svd algorithm we obtain a dictionary that describes the image content effectively two training options are considered using the corrupted image itself or training on a corpus of high quality image database since the k svd is limited in handling small image patches we extend its deployment to arbitrary image sizes by defining a global image prior that forces sparsity over patches in every location in the image we show how such bayesian treatment leads to a simple and effective denoising algorithm this leads to a state of the art denoising performance equivalent and sometimes surpassing recently published leading alternative denoising methods

a modified k means algorithm for circular invariant clustering several important pattern recognition applications are based on feature vector extraction and vector clustering directional patterns are commonly represented by rotation variant vectors f sub d formed from features uniformly extracted in m directions it is often desirable that pattern recognition algorithms are invariant under pattern rotation this paper introduces a distance measure and a k means based algorithm namely circular k means ck means to cluster vectors containing directional information such as f sub d in a circular shift invariant manner a circular shift of f sub d corresponds to pattern rotation thus the algorithm is rotation invariant an efficient fourier domain representation of the proposed measure is presented to reduce computational complexity a split and merge approach smck means suited to the proposed ck means technique is proposed to reduce the possibility of converging at local minima and to estimate the correct number of clusters experiments performed for textural images illustrate the superior performance of the proposed algorithm for clustering directional vectors f sub d compared to the alternative approach that uses the original k means and rotation invariant feature vectors transformed from f sub d

scale affine invariant interest point detectors in this paper we propose a novel approach for detecting interest points invariant to scale and affine transformations our scale and affine invariant detectors are based on the following recent results 1 interest points extracted with the harris detector can be adapted to affine transformations and give repeatable results geometrically stable 2 the characteristic scale of a local structure is indicated by a local extremum over scale of normalized derivatives the laplacian 3 the affine shape of a point neighborhood is estimated based on the second moment matrix our scale invariant detector computes a multi scale representation for the harris interest point detector and then selects points at which a local measure the laplacian is maximal over scales this provides a set of distinctive points which are invariant to scale rotation and translation as well as robust to illumination changes and limited changes of viewpoint the characteristic scale determines a scale invariant region for each point we extend the scale invariant detector to affine invariance by estimating the affine shape of a point neighborhood an iterative algorithm modifies location scale and neighborhood of each point and converges to affine invariant points this method can deal with significant affine transformations including large scale changes the characteristic scale and the affine shape of neighborhood determine an affine invariant region for each point we present a comparative evaluation of different detectors and show that our approach provides better results than existing methods the performance of our detector is also confirmed by excellent matching resultss the image is described by a set of scale affine invariant descriptors computed on the regions associated with our points

the design and use of steerable filters the authors present an efficient architecture to synthesize filters of arbitrary orientations from linear combinations of basis filters allowing one to adaptively steer a filter to any orientation and to determine analytically the filter output as a function of orientation steerable filters may be designed in quadrature pairs to allow adaptive control over phase as well as orientation the authors show how to design and steer the filters and present examples of their use in the analysis of orientation and phase angularly adaptive filtering edge detection and shape from shading one can also build a self similar steerable pyramid representation the same concepts can be generalized to the design of 3 d steerable filters

a generalized l1 greedy algorithm for image reconstruction in ct the sparse vector solutions for an underdetermined system of linear equations ax b have many applications in signal recovery and image reconstruction in tomography under certain conditions the sparsest solution can be found by solving a constrained l 1 minimization problem min x 1 subject to ax b recently the reweighted l 1 minimization and l 1 greedy algorithm have been introduced to improve the convergence of the l 1 minimization problem as an extension a generalized l 1 greedy algorithm for computerized tomography ct is proposed in this paper it is implemented as a generalized total variation minimization for images with sparse gradients in ct numerical experiments are also given to illustrate the advantage of the new algorithm

classification via sparse representation of steerable wavelet frames on grassmann manifold application to target recognition in sar image automatic target recognition has been widely studied over the years yet it is still an open problem the main obstacle consists in extended operating conditions e g depression angle change configuration variation articulation and occlusion to deal with them this paper proposes a new classification strategy we develop a new representation model via the steerable wavelet frames the proposed representation model is entirely viewed as an element on grassmann manifolds to achieve target classification we embed grassmann manifolds into an implicit reproducing kernel hilbert space rkhs where the kernel sparse learning can be applied specifically the mappings of training sample in rkhs are concatenated to form an overcomplete dictionary it is then used to encode the counterpart of query as a linear combination of its atoms by designed grassmann kernel function it is capable to obtain the sparse representation from which the inference can be reached the novelty of this paper comes from 1 the development of representation model by the set of directional components of riesz transform 2 the quantitative measure of similarity for proposed representation model by grassmann metric and 3 the generation of global kernel function by grassmann kernel extensive comparative studies are performed to demonstrate the advantage of proposed strategy

joint entropy of continuously differentiable ultrasonic waveforms this study is based on an extension of the concept of joint entropy of two random variables to continuous functions such as backscattered ultrasound for two continuous random variables x and y the joint probability density p x y is ordinarily a continuous function of x and y that takes on values in a two dimensional region of the real plane however in the case where x f t and y g t are both continuously differentiable functions x and y are concentrated exclusively on a curve t f t g t in the x y plane this concentration can only be represented using a mathematically singular object such as a schwartz distribution its use for imaging requires a coarse graining operation which is described in this study subsequently removal of the coarse graining parameter is accomplished using the ergodic theorem the resulting expression for joint entropy is applied to several data sets showing the utility of the concept for both materials characterization and detection of targeted liquid nanoparticle ultrasonic contrast agents in all cases the sensitivity of these techniques matches or exceeds sometimes by a factor of two that demonstrated in previous studies that employed signal energy or alternate entropic quantities

the laplacian pyramid as a compact image code we describe a technique for image encoding in which local operators of many scales but identical shape serve as the basis functions the representation differs from established techniques in that the code elements are localized in spatial frequency as well as in space pixel to pixel correlations are first removed by subtracting a lowpass filtered copy of the image from the image itself the result is a net data compression since the difference or error image has low variance and entropy and the low pass filtered image may represented at reduced sample density further data compression is achieved by quantizing the difference image these steps are then repeated to compress the low pass image iteration of the process at appropriately expanded scales generates a pyramid data structure the encoding process is equivalent to sampling the image with laplacian operators of many scales thus the code tends to enhance salient image features a further advantage of the present code is that it is well suited for many image analysis tasks as well as for image compression fast algorithms are described for coding and decoding

3d steerable wavelets in practice we introduce a systematic and practical design for steerable wavelet frames in 3d our steerable wavelets are obtained by applying a 3d version of the generalized riesz transform to a primary isotropic wavelet frame the novel transform is self reversible tight frame and its elementary constituents riesz wavelets can be efficiently rotated in any 3d direction by forming appropriate linear combinations moreover the basis functions at a given location can be linearly combined to design custom and adaptive steerable wavelets the features of the proposed method are illustrated with the processing and analysis of 3d biomedical data in particular we show how those wavelets can be used to characterize directional patterns and to detect edges by means of a 3d monogenic analysis we also propose a new inverse problem formalism along with an optimization algorithm for reconstructing 3d images from a sparse set of wavelet domain edges the scheme results in high quality image reconstructions which demonstrate the feature reduction ability of the steerable wavelets as well as their potential for solving inverse problems

analysis of digitized waveforms using shannon entropy ii high speed algorithms based on green s functions an improved approach for calculating continuous waveform entropy hc is described this approach is based on the use of the green s function for the one dimensional differential operator l u u with boundary conditions u 1 u 0 0 the method is applied to data acquired from a plexiglas block having artificial defects in order to produce an entropy image the defect contrast obtained with the new green s function computation is then compared with that obtained from a previously described fourier series approach it is found that the new approach produces the same or higher image contrast in a time that is roughly three orders of magnitude smaller than that required by the fourier series method this speedup arises from two sources the green s function approach has greater inherent immunity to noise and the number of calculations required in the green s function approach is much smaller than the number required in the fourier series approach

dictionary learning algorithms for sparse representation algorithms for data driven learning of domain specific overcomplete dictionaries are developed to obtain maximum likelihood and maximum a posteriori dictionary estimates based on the use of bayesian models with concave schur concave csc negative log priors such priors are appropriate for obtaining sparse representations of environmental signals within an appropriately chosen environmentally matched dictionary the elements of the dictionary can be interpreted as concepts features or words capable of succinct expression of events encountered in the environment the source of the measured signals this is a generalization of vector quantization in that one is interested in a description involving a few dictionary entries the proverbial 25 words or less but not necessarily as succinct as one entry to learn an environmentally adapted dictionary capable of concise expression of signals generated by the environment we develop algorithms that iterate between a representative set of sparse representations found by variants of focuss and an update of the dictionary using these sparse representations experiments were performed using synthetic data and natural images for complete dictionaries we demonstrate that our algorithms have improved performance over other independent component analysis ica methods measured in terms of signal to noise ratios of separated sources in the overcomplete case we show that the true underlying dictionary and sparse sources can be accurately recovered in tests with natural images learned overcomplete dictionaries are shown to have higher coding efficiency than complete dictionaries that is images encoded with an overcomplete dictionary have both higher compression fewer bits per pixel and higher accuracy lower mean square error

properties of an entropy based signal receiver with an application to ultrasonic molecular imaging qualitative and quantitative properties of the finite part hf of the shannon entropy of a continuous waveform f t in the continuum limit are derived in order to illuminate its use for waveform characterization simple upper and lower bounds on hf based on features of f t are defined quantitative criteria for a priori estimation of the average case variation of hf and logef where ef is the signal energy of f t are also derived these provide relative sensitivity estimates that could be used to prospectively choose optimal imaging strategies in real time ultrasonic imaging machines where system bandwidth is often pushed to its limits to demonstrate the utility of these sensitivity relations for this application a study designed to assess the feasibility of identification of angiogenic neovasculature targeted with perfluorocarbon nanoparticles that specifically bind to v 3 integrin expression in tumors was performed the outcome of this study agrees with the prospective sensitivity estimates tha

accurate image rotation using hermite expansions in this paper we propose an approach for the accurate rotation of a digital image using hermite expansions this exploits the fact that if a 2 d continuous bandlimited hermite expansion is rotated the resulting function can be expressed as a hermite expansion with the same bandlimit furthermore the hermite coefficients of the initial 2 d expansion and the rotated expansion are mapped through an invertible linear relationship two efficient methods to compute the mapping between hermite coefficients during rotation are proposed we also propose a method for connecting the hermite expansion and a discrete image using this method we can obtain the hermite expansion from a discrete image and vice versa combining these techniques we propose new methods for the rotation of discrete images we assess the accuracy of our methods and compare them with an existing fft based method implementing three shears we find that the method proposed here consistently has better accuracy than the fft based method

ideal spatial adaptation by wavelet shrinkage summary with ideal spatial adaptation an oracle furnishes information about how best to adapt a spatially variable estimator whether piecewise constant piecewise polynomial variable knot spline or variable bandwidth kernel to the unknown function estimation with the aid of an oracle offers dramatic advantages over traditional linear estimation by nonadaptive kernels however it is a priori unclear whether such performance can be obtained by a procedure relying on the data alone we describe a new principle for spatially adaptive estimation selective wavelet reconstruction we show that variable knot spline fits and piecewise polynomial fits when equipped with an oracle to select the knots are not dramatically more powerful than selective wavelet reconstruction with an oracle we develop a practical spatially adaptive method riskshrink which works by shrinkage of empirical wavelet coefficients riskshrink mimics the performance of an oracle for selective wavelet reconstruction as well as it is possible to do so a new inequality in multivariate normal decision theory which we call the oracle inequality shows that attained performance differs from ideal performance by at most a factor of approximately 2 log n where n is the sample size moreover no estimator can give a better guarantee than this within the class of spatially adaptive procedures riskshrink is essentially optimal relying only on the data it comes within a factor log 2 n of the performance of piecewise polynomial and variableknot spline methods equipped with an oracle in contrast it is unknown how or if piecewise polynomial methods could be made to function this well when denied access to an oracle and forced to rely on data alone

the hermite transform applications it is demonstrated how the hermite transform can be used for image coding and analysis hierarchical coding structures based on increasingly specified basic patterns i e general 2 d patterns general 1 d patterns and specific 1 d patterns such as edges and corners are presented in the image coding application the relation with existing pyramid coders is described a new coding scheme based on local one dimensional image approximations is introduced in the image analysis application the relation between the hermite transform and existing line edge detection schemes is described it is shown that by concentrating on more specific patterns the coding efficiency can be increased since fewer coefficients have to be coded meanwhile sufficient descriptive power can be maintained for approximating the most interesting features in natural images

complete discrete 2 d gabor transforms by neural networks for image analysis and compression a three layered neural network is described for transforming two dimensional discrete signals into generalized nonorthogonal 2 d gabor representations for image analysis segmentation and compression these transforms are conjoint spatial spectral representations which provide a complete image description in terms of locally windowed 2 d spectral coordinates embedded within global 2 d spatial coordinates in the present neural network approach based on interlaminar interactions involving two layers with fixed weights and one layer with adjustable weights the network finds coefficients for complete conjoint 2 d gabor transforms without restrictive conditions in wavelet expansions based on a biologically inspired log polar ensemble of dilations rotations and translations of a single underlying 2 d gabor wavelet template image compression is illustrated with ratios up to 20 1 also demonstrated is image segmentation based on the clustering of coefficients in the complete 2 d gabor transform

flexible spatial configuration of local image features local image features have been designed to be informative and repeatable under rigid transformations and illumination deformations even though current state of the art local image features present a high degree of repeatability their local appearance alone usually does not bring enough discriminative power to support a reliable matching resulting in a relatively high number of mismatches in the correspondence set formed during the data association procedure as a result geometric filters commonly based on global spatial configuration have been used to reduce this number of mismatches however this approach presents a trade off between the effectiveness to reject mismatches and the robustness to nonrigid deformations in this paper we propose two geometric filters based on a semilocal spatial configuration of local features that are designed to be robust to nonrigid deformations and to rigid transformations without compromising its efficacy to reject mismatches we compare our methods to the hough transform which is an efficient and effective mismatch rejection step based on the global spatial configuration of features in these comparisons our methods are shown to be more effective in the task of rejecting mismatches for rigid transformations and nonrigid deformations at comparable time complexity figures finally we demonstrate how we can integrate these methods in a probabilistic recognition system such that the final verification step uses not only the similarity between features but also their semilocal configuration

contrast context histogram an efficient discriminating local descriptor for object recognition and image matching in this paper we propose a new invariant local descriptor called the contrast context histogram cch for image matching and object recognition by representing the contrast distributions of a local region it serves as a distinctive local descriptor of the region our experiments demonstrate that contrast based local descriptors can represent local regions with more compact histogram bins because of its high matching accuracy and efficient computation the cch has the potential to be used in a number of real time applications

curvilinear structure analysis by ranking the orientation responses of path operators the analysis of thin curvilinear objects in 3d images is a complex and challenging task in this article we introduce a new non linear operator called rorpo ranking the orientation responses of path operators inspired by the multidirectional paradigm currently used in linear filtering for thin structure analysis rorpo is built upon the notion of path operator from mathematical morphology this operator unlike most operators commonly used for 3d curvilinear structure analysis is discrete non linear and non local from this new operator two main curvilinear structure characteristics can be estimated an intensity feature that can be assimilated to a quantitative measure of curvilinearity and a directional feature providing a quantitative measure of the structure s orientation we provide a full description of the structural and algorithmic details for computing these two features from rorpo and we discuss computational issues we experimentally assess rorpo by comparison with three of the most popular curvilinear structure analysis filters namely frangi vesselness optimally oriented flux and hybrid diffusion with continuous switch in particular we show that our method provides up to 8 percent more true positive and 50 percent less false positives than the next best method on synthetic and real 3d images

a tiling of phase space through self convolution a novel tiling for fourier half space is constructed using erlang functions this tiling is constructed by shifts in discrete time and self convolutions in the fourier domain we show that a subframe of this complex wavelet family can be used for the detection and characterization of events step and peaks in one dimensional 1 d discrete time signals

the sure let approach to image denoising we propose a new approach to image denoising based on the image domain minimization of an estimate of the mean squared error stein s unbiased risk estimate sure unlike most existing denoising algorithms using the sure makes it needless to hypothesize a statistical model for the noiseless image a key point of our approach is that although the nonlinear processing is performed in a transformed domain typically an undecimated discrete wavelet transform but we also address nonorthonormal transforms this minimization is performed in the image domain indeed we demonstrate that when the transform is a ldquotightrdquo frame an undecimated wavelet transform using orthonormal filters separate subband minimization yields substantially worse results in order for our approach to be viable we add another principle that the denoising process can be expressed as a linear combination of elementary denoising processes linear expansion of thresholds let armed with the sure and let principles we show that a denoising algorithm merely amounts to solving a linear system of equations which is obviously fast and efficient quite remarkably the very competitive results obtained by performing a simple threshold image domain sure optimized on the undecimated haar wavelet coefficients show that the sure let principle has a huge potential

maximum likelihood orientation estimation of 1 d patterns in laguerre gauss subspaces a method for measuring the orientation of linear 1 d patterns based on a local expansion with laguerre gauss circular harmonic lg ch functions is presented it lies on the property that the polar separable lg ch functions span the same space as the 2 d cartesian separable hermite gauss 2 d hg functions exploiting the simple steerability of the lg ch functions and the peculiar block linear relationship among the two expansion coefficients sets maximum likelihood ml estimates of orientation and cross section parameters of 1 d patterns are obtained projecting them in a proper subspace of the 2 d hg family it is shown in this paper that the conditional ml solution derived by elimination of the cross section parameters surprisingly yields the same asymptotic accuracy as the ml solution for known cross section parameters the accuracy of the conditional ml estimator is compared to the one of state of art solutions on a theoretical basis and via simulation trials a thorough proof of the key relationship between the lg ch and the 2 d hg expansions is also provided

the curvelet transform for image denoising we describe approximate digital implementations of two new mathematical transforms namely the ridgelet transform and the curvelet transform our implementations offer exact reconstruction stability against perturbations ease of implementation and low computational complexity a central tool is fourier domain computation of an approximate digital radon transform we introduce a very simple interpolation in the fourier space which takes cartesian samples and yields samples on a rectopolar grid which is a pseudo polar sampling set based on a concentric squares geometry despite the crudeness of our interpolation the visual performance is surprisingly good our ridgelet transform applies to the radon transform a special overcomplete wavelet pyramid whose wavelets have compact support in the frequency domain our curvelet transform uses our ridgelet transform as a component step and implements curvelet subbands using a filter bank of a spl grave trous wavelet filters our philosophy throughout is that transforms should be overcomplete rather than critically sampled we apply these digital transforms to the denoising of some standard images embedded in white noise in the tests reported here simple thresholding of the curvelet coefficients is very competitive with state of the art techniques based on wavelets including thresholding of decimated or undecimated wavelet transforms and also including tree based bayesian posterior mean methods moreover the curvelet reconstructions exhibit higher perceptual quality than wavelet based reconstructions offering visually sharper images and in particular higher quality recovery of edges and of faint linear and curvilinear features existing theory for curvelet and ridgelet transforms suggests that these new approaches can outperform wavelet methods in certain image reconstruction problems the empirical results reported here are in encouraging agreement

training less color object recognition for autonomous robotics abstract this paper looks at the challenge of object recognition from the perspective of achieving the final goals of practical real world autonomous robot vision applications that of identifying a target of interest in the robot s field of view and properly localizing its position in preparation for the higher level goal of tracking and navigation a unified framework is introduced that combines a multidimensional feature histogram approach with a multiscale pyramid approach for training less color object recognition and localization with direct application to autonomous robotic agents this framework addresses the high computational cost associated with multidimensional processing by deriving the multidimensional laplacian feature histogram pyramid a novel approach to a unified multidimensional multiscale histogram representation furthermore a taylor series formulation is employed to combine the multiscale levels of the multidimensional laplacian feature histogram pyramid into one efficient multidimensional multiscale laplacian taylor feature histogram for rapid object recognition and localization the paper describes the criteria for target detection and localization by autonomous robots and how this newly developed framework fits these needs comparative results demonstrate the robustness of this recognition framework to noise and localization of target objects in cluttered scenes

noise analysis of a new singularity index we analyze the noise sensitivity of a new singularity index that was designed to detect impulse singularities in signals of arbitrary dimensionality while rejecting step like singularities see muralidhar ieee signal process lett vol 20 no 1 pp 7 10 2013 and muralidhar proc ieee int conf image process 2012 for example the index responds strongly to curvilinear masses ridges in images while weakly to jump discontinuities edges we analyze the detection power of the index in the presence of noise our analysis is geared towards answering the following questions a in the presence of noise only what is the probability of falsely detecting an impulse given a threshold b given an impulse submerged in noise what is the probability of detecting it given a threshold and c since the index is designed to be edge suppressing what is the probability of incorrectly detecting an edge submerged in noise given a threshold we compare the detection power of the index with that of a nominal impulse detector the second derivative operator simulations and example applications in 1 d and 2 d reveal the efficacy of the new singularity index for correctly detecting impulses submerged in noise while suppressing edges a software version of the 2 d singularity index can be downloaded from http live ece utexas edu research singularityindex singularityindexcode zip

a computational approach to edge detection this paper describes a computational approach to edge detection the success of the approach depends on the definition of a comprehensive set of goals for the computation of edge points these goals must be precise enough to delimit the desired behavior of the detector while making minimal assumptions about the form of the solution we define detection and localization criteria for a class of edges and present mathematical forms for these criteria as functionals on the operator impulse response a third criterion is then added to ensure that the detector has only one response to a single edge we use the criteria in numerical optimization to derive detectors for several common image features including step edges on specializing the analysis to step edges we find that there is a natural uncertainty principle between detection and localization performance which are the two main goals with this principle we derive a single operator shape which is optimal at any scale the optimal detector has a simple approximate implementation in which edges are marked at maxima in gradient magnitude of a gaussian smoothed image we extend this simple detector using operators of several widths to cope with different signal to noise ratios in the image we present a general method called feature synthesis for the fine to coarse integration of information from operators at different scales finally we show that step edge detector performance improves considerably as the operator point spread function is extended along the edge

steerable wavelet frames based on the riesz transform we consider an extension of the 1 d concept of analytical wavelet to n d which is by construction compatible with rotations this extension called a monogenic wavelet yields a decomposition of the wavelet coefficients into amplitude phase and phase direction the monogenic wavelet is based on the hypercomplex monogenic signal which is defined using riesz transforms and perfectly isotropic wavelets frames employing the new concept of clifford frames we can show that the monogenic wavelet generates a wavelet frame furthermore this approach yields wavelet frames that are steerable with respect to direction applications to descreening and contrast enhancement illustrate the versatility of this approach to image analysis and reconstruction

local shading analysis local analysis of image shading in the absence of prior knowledge about the viewed scene may be used to provide information about the scene the following has been proved every image point has the same image intensity and first and second derivatives as the image of some point on a lambertian surface with principal curvatures of equal magnitude further if the principal curvatures are assumed to be equal there is a unique combination of image formation parameters up to a mirror reversal that will produce a particular set of image intensity and first and second derivatives a solution for the unique combination of surface orientation etc is presented this solution has been extended to natural imagery by using general position and regional constraints to obtain estimates of the following surface orientation at each image point the qualitative type of the surface i e whether the surface is planar cylindrical convex concave or saddle the illuminant direction within a region algorithms to recover illuminant direction and estimate surface orientation have been evaluated on both natural and synthesized images and have been found to produce useful information about the scene

signal denoising and ultrasonic flaw detection via overcomplete and sparse representations sparse signal representations from overcomplete dictionaries are the most recent technique in the signal processing community applications of this technique extend into many fields in this paper this technique is utilized to cope with ultrasonic flaw detection and noise suppression problem in particular a noisy ultrasonic signal is decomposed into sparse representations using a sparse bayesian learning algorithm and an overcomplete dictionary customized from a gabor dictionary by incorporating some a priori information of the transducer used nonlinear postprocessing including thresholding and pruning is then applied to the decomposed coefficients to reduce the noise contribution and extract the flaw information because of the high compact essence of sparse representations flaw echoes are packed into a few significant coefficients and noise energy is likely scattered all over the dictionary atoms generating insignificant coefficients this property greatly increases the efficiency of the pruning and thresholding operations and is extremely useful for detecting flaw echoes embedded in background noise the performance of the proposed approach is verified experimentally and compared with the wavelet transform signal processor experimental results to detect ultrasonic flaw echoes contaminated by white gaussian additive noise or correlated noise are presented in the paper

a differential model of the complex cell the receptive fields of simple cells in the visual cortex can be understood as linear filters these filters can be modeled by gabor functions or gaussian derivatives gabor functions can also be combined in an energy model of the complex cell response this letter proposes an alternative model of the complex cell based on gaussian derivatives it is most important to account for the insensitivity of the complex response to small shifts of the image the new model uses a linear combination of the first few derivative filters at a single position to approximate the first derivative filter at a series of adjacent positions the maximum response over all positions gives a signal that is insensitive to small shifts of the image this model unlike previous approaches is based on the scale space theory of visual processing in particular the complex cell is built from filters that respond to the 2d differential structure of the image the computational aspects of the new model are studied in one and two dimensions using the steerability of the gaussian derivatives the response of the model to basic images such as edges and gratings is derived formally the response to natural images is also evaluated using statistical measures of shift insensitivity the neural implementation and predictions of the model are discussed

focus of attention from local color symmetries in this paper a continuous valued measure for local color symmetry is introduced the new algorithm is an extension of the successful gray value based symmetry map proposed by reisfeld et al the use of color facilitates the detection of focus points fps on objects that are difficult to detect using gray value contrast only the detection of fps is aimed at guiding the attention of an object recognition system therefore fps have to fulfill three major requirements stability distinctiveness and usability the proposed algorithm is evaluated for these criteria and compared with the gray value based symmetry measure and two other methods from the literature stability is tested against noise object rotation and variations of lighting as a measure for the distinctiveness of fps the principal components of fp centered windows are compared with those of windows at randomly chosen points on a large database of natural images finally usability is evaluated in the context of an object recognition task

translation invariant contourlet transform and its application to image denoising most subsampled filter banks lack the feature of translation invariance which is an important characteristic in denoising applications in this paper we study and develop new methods to convert a general multichannel multidimensional filter bank to a corresponding translation invariant ti framework in particular we propose a generalized algorithme agrave trous which is an extension of the algorithme agrave trous introduced for 1 d wavelet transforms using the proposed algorithm as well as incorporating modified versions of directional filter banks we construct the ti contourlet transform tict to reduce the high redundancy and complexity of the tict we also introduce semi translation invariant contourlet transform stict then we employ an adapted bivariate shrinkage scheme to the stict to achieve an efficient image denoising approach our experimental results demonstrate the benefits and potential of the proposed denoising approach complexity analysis and efficient realization of the proposed ti schemes are also presented

the design of approximate hilbert transform pairs of wavelet bases several authors have demonstrated that significant improvements can be obtained in wavelet based signal processing by utilizing a pair of wavelet transforms where the wavelets form a hilbert transform pair this paper describes design procedures based on spectral factorization for the design of pairs of dyadic wavelet bases where the two wavelets form an approximate hilbert transform pair both orthogonal and biorthogonal fir solutions are presented as well as iir solutions in each case the solution depends on an allpass filter having a flat delay response the design procedure allows for an arbitrary number of vanishing wavelet moments to be specified a matlab program for the procedure is given and examples are also given to illustrate the results

image filtering edge detection and edge tracing using fuzzy reasoning we characterize the problem of detecting edges in images as a fuzzy reasoning problem the edge detection problem is divided into three stages filtering detection and tracing images are filtered by applying fuzzy reasoning based on local pixel characteristics to control the degree of gaussian smoothing filtered images are then subjected to a simple edge detection algorithm which evaluates the edge fuzzy membership value for each pixel based on local image characteristics finally pixels having high edge membership are traced and assembled into structures again using fuzzy reasoning to guide the tracing process the filtering detection and tracing algorithms are tested on several test images comparison is made with a standard edge detection technique

multiresolution monogenic signal analysis using the riesz laplace wavelet transform the monogenic signal is the natural 2d counterpart of the 1d analytic signal we propose to transpose the concept to the wavelet domain by considering a complexified version of the riesz transform which has the remarkable property of mapping a real valued primary wavelet basis of l2 r2 into a complex one the riesz operator is also steerable in the sense that it give access to the hilbert transform of the signal along any orientation having set those foundations we specify a primary polyharmonic spline wavelet basis of l2 r2 that involves a single mexican hat like mother wavelet laplacian of a b spline the important point is that our primary wavelets are quasi isotropic they behave like multiscale versions of the fractional laplace operator from which they are derived which ensures steerability we propose to pair these real valued basis functions with their complex riesz counterparts to specify a multiresolution monogenic signal analysis this yields a representation where each wavelet index is associated with a local orientation an amplitude and a phase we give a corresponding wavelet domain method for estimating the underlying instantaneous frequency we also provide a mechanism for improving the shift and rotation invariance of the wavelet decomposition and show how to implement the transform efficiently using perfect reconstruction filterbanks we illustrate the specific feature extraction capabilities of the representation and present novel examples of wavelet domain processing in particular a robust tensor based analysis of directional image patterns the demodulation of interferograms and the reconstruction of digital holograms

a morphological affine and galilean invariant scale space for movies we study a model of multiscale analysis or scale space applied to movies this model comes from a thorough formalization that has been done in the theory of scale space of static image this formulation has led one to associate with each multiscale analysis a partial differential equation pde we examine the case of movies and insist on the motion aspects more precisely it has been proved by alvarez guichard lions and morel 1993 that there exists a unique affine and morphological and galilean invariant scale space for movies the amg model this model is described by a partial differential equation we focus on terms appearing in that equation we show that this model provides a reliable definition of an optical multiscale acceleration at the practical level scale is interpreted as a way of characterizing reliable trajectories as we prove by experiments the amg model is a riddle for decimating spurious trajectories due to any kinds of nonadditive impurities and noise simple discrete formulae are given to implement the model

retinal vessel centerline extraction using multiscale matched filters confidence and edge measures motivated by the goals of improving detection of low contrast and narrow vessels and eliminating false detections at nonvascular structures a new technique is presented for extracting vessels in retinal images the core of the technique is a new likelihood ratio test that combines matched filter responses confidence measures and vessel boundary measures matched filter responses are derived in scale space to extract vessels of widely varying widths a vessel confidence measure is defined as a projection of a vector formed from a normalized pixel neighborhood onto a normalized ideal vessel profile vessel boundary measures and associated confidences are computed at potential vessel boundaries combined these responses form a six dimensional measurement vector at each pixel a training technique is used to develop a mapping of this vector to a likelihood ratio that measures the vesselness at each pixel results comparing this vesselness measure to matched filters alone and to measures based on the hessian of intensities show substantial improvements both qualitatively and quantitatively the hessian can be used in place of the matched filter to obtain similar but less substantial improvements or to steer the matched filter by preselecting kernel orientations finally the new vesselness likelihood ratio is embedded into a vessel tracing framework resulting in an efficient and effective vessel centerline extraction algorithm

recognition by symmetry derivatives and the generalized structure tensor we suggest a set of complex differential operators that can be used to produce and filter dense orientation tensor fields for feature extraction matching and pattern recognition we present results on the invariance properties of these operators that we call symmetry derivatives these show that in contrast to ordinary derivatives all orders of symmetry derivatives of gaussians yield a remarkable invariance they are obtained by replacing the original differential polynomial with the same polynomial but using ordinary coordinates x and y corresponding to partial derivatives moreover the symmetry derivatives of gaussians are closed under the convolution operator and they are invariant to the fourier transform the equivalent of the structure tensor representing and extracting orientations of curve patterns had previously been shown to hold in harmonic coordinates in a nearly identical manner as a result positions orientations and certainties of intricate patterns e g spirals crosses parabolic shapes can be modeled by use of symmetry derivatives of gaussians with greater analytical precision as well as computational efficiency since gaussians and their derivatives are utilized extensively in image processing the revealed properties have practical consequences for local orientation based feature extraction the usefulness of these results is demonstrated by two applications 1 tracking cross markers in long image sequences from vehicle crash tests and 2 alignment of noisy fingerprints

local energy pattern for texture classification using self adaptive quantization thresholds local energy pattern a statistical histogram based representation is proposed for texture classification first we use normalized local oriented energies to generate local feature vectors which describe the local structures distinctively and are less sensitive to imaging conditions then each local feature vector is quantized by self adaptive quantization thresholds determined in the learning stage using histogram specification and the quantized local feature vector is transformed to a number by n nary coding which helps to preserve more structure information during vector quantization finally the frequency histogram is used as the representation feature the performance is benchmarked by material categorization on kth tips and kth tips2 a databases our method is compared with typical statistical approaches such as basic image features local binary pattern lbp local ternary pattern completed lbp weber local descriptor and vz algorithms vz mr8 and vz joint the results show that our method is superior to other methods on the kth tips2 a database and achieving competitive performance on the kth tips database furthermore we extend the representation from static image to dynamic texture and achieve favorable recognition results on the university of california at los angeles ucla dynamic texture database

contrast sensitivity of the wavelet dual tree complex wavelet curvelet and steerable pyramid transforms accurate estimation of the contrast sensitivity of the human visual system is crucial for perceptually based image processing in applications such as compression fusion and denoising conventional contrast sensitivity functions csfs have been obtained using fixed sized gabor functions however the basis functions of multiresolution decompositions such as wavelets often resemble gabor functions but are of variable size and shape therefore to use the conventional csfs in such cases is not appropriate we have therefore conducted a set of psychophysical tests in order to obtain the csf for a range of multiresolution transforms the discrete wavelet transform the steerable pyramid the dual tree complex wavelet transform and the curvelet transform these measures were obtained using contrast variation of each transforms basis functions in a 2afc experiment combined with an adapted version of the quest psychometric function method the results enable future image processing applications that exploit these transforms such as signal fusion superresolution processing denoising and motion estimation to be perceptually optimized in a principled fashion the results are compared with an existing vision model hdr vdp2 and are used to show quantitative improvements within a denoising application compared with using conventional csf values

generalized sampling expansion for functions on the sphere functions on the sphere appear in several applications including geodesics imaging and acoustics sampling of these functions may result in aliasing if the sampling condition is not met the generalized sampling expansion introduced by papoulis enables the reconstruction of a band limited function sampled at a frequency lower than the nyquist frequency using the outputs of several linear time invariant systems this paper formulates the generalized sampling expansion for functions on the sphere using spherical harmonics decomposition facilitating sub nyquit sampling without aliasing error an analysis of linear systems on the sphere and the aliasing phenomenon in the spherical harmonics domain is presented examples demonstrating the performance of the method and its limitations are studied

two view line matching algorithm based on context and appearance in low textured images a novel approach for line detection and matching is proposed aimed at achieving good performance with low textured scenes under uncontrolled illumination conditions line detection is performed by means of phase based edge detector over gaussian scale space followed by a multi scale fusion stage which has been proven to be profitable in minimizing the number of fragmented and overlapped segments line matching is performed by an iterative process that uses structural information collected through the use of different line neighborhoods making the set of matched lines grow robustly at each iteration results show that this approach is suitable to deal with low textured scenes and also robust under a wide variety of image transformations highlightsa novel approach for line detection and matching is proposed lines are detected through a fusion process over the gaussian scale space lines are matched based on their appearance geometric properties and neighborhoods our proposal minimizes the number of fragmented and overlapped segments our proposal is suitable for low textured scenes under uncontrolled illumination

image representations using multiscale differential operators differential operators have been widely used for multiscale geometric descriptions of images the efficient computation of these differential operators is always desirable moreover it has not been clear whether such representations are invertible for certain applications it is usually required that such representations should be invertible so that one can facilitate the processing of information in the transform domain and then recover it in this paper such problems are studied we consider multiscale differential representations of images using different types of operators such as the directional derivative operators and laplacian operators in particular we provide a general approach to represent images by their multiscale and multidirectional derivative components for practical implementation efficient pyramid like algorithms are derived using the spline technique for both the the composition and reconstruction of images it is shown that using these representations various meaningful geometric information of images can be extracted at multiple scales therefore these representations can be used for edge based image processing purposes furthermore the intrinsic relationships of the proposed representations with the compact wavelet models and some classical multiscale approaches are also elucidated in the paper

discriminative learning of local image descriptors in this paper we explore methods for learning local image descriptors from training data we describe a set of building blocks for constructing descriptors which can be combined together and jointly optimized so as to minimize the error of a nearest neighbor classifier we consider both linear and nonlinear transforms with dimensionality reduction and make use of discriminant learning techniques such as linear discriminant analysis lda and powell minimization to solve for the parameters using these techniques we obtain descriptors that exceed state of the art performance with low dimensionality in addition to new experiments and recommendations for descriptor learning we are also making available a new and realistic ground truth data set based on multiview stereo data

lapb locally adaptive patch based wavelet domain edge preserving image denoising abstract image denoising is one of the most diversified research areas in the field of image processing and computer vision it is highly desirable for a denoising technique to preserve important image features such as edges corners and other sharp structures of an image after denoising wavelet transforms show excellent proficiency in providing efficient edge preserving image denoising due to their capability of suppressing noisy signals from an image this paper presents a novel edge preserving image denoising technique based on wavelet transforms the multi level decomposition of the noisy image is carried out to transform the data into the wavelet domain a locally adaptive patch based lapb thresholding scheme is used to effectively reduce noise while preserving relevant features of the original image experimental results on benchmark test images demonstrate that the proposed method achieves competitive denoising performance in comparison to various state of the art algorithms

analysis of focus measure operators for shape from focus shape from focus sff has widely been studied in computer vision as a passive depth recovery and 3d reconstruction method one of the main stages in sff is the computation of the focus level for every pixel of an image by means of a focus measure operator in this work a methodology to compare the performance of different focus measure operators for shape from focus is presented and applied the selected operators have been chosen from an extensive review of the state of the art the performance of the different operators has been assessed through experiments carried out under different conditions such as image noise level contrast saturation and window size such performance is discussed in terms of the working principles of the analyzed operators

a least squares framework for component analysis over the last century component analysis ca methods such as principal component analysis pca linear discriminant analysis lda canonical correlation analysis cca locality preserving projections lpp and spectral clustering sc have been extensively used as a feature extraction step for modeling classification visualization and clustering ca techniques are appealing because many can be formulated as eigen problems offering great potential for learning linear and nonlinear representations of data in closed form however the eigen formulation often conceals important analytic and computational drawbacks of ca techniques such as solving generalized eigen problems with rank deficient matrices e g small sample size problem lacking intuitive interpretation of normalization factors and understanding commonalities and differences between ca methods this paper proposes a unified least squares framework to formulate many ca methods we show how pca lda cca lpp sc and its kernel and regularized extensions correspond to a particular instance of least squares weighted kernel reduced rank regression ls wkrrr the ls wkrrr formulation of ca methods has several benefits 1 provides a clean connection between many ca techniques and an intuitive framework to understand normalization factors 2 yields efficient numerical schemes to solve ca techniques 3 overcomes the small sample size problem 4 provides a framework to easily extend ca methods we derive weighted generalizations of pca lda sc and cca and several new ca techniques

range image statistics can explain the anomalous perception of length a long standing puzzle in visual perception is that the apparent extent of a spatial interval e g the distance between two points or the length of a line does not simply accord with the length of the stimulus but varies as a function of orientation in the retinal image here we show that this anomaly can be explained by the statistical relationship between the length of retinal projections and the length of their real world sources using a laser range scanner we acquired a database of natural images that included the three dimensional location of every point in the scenes an analysis of these range images showed that the average length of a physical interval in three dimensional space changes systematically as a function of the orientation of the corresponding interval in the projected image the variation being in good agreement with perceived length this evidence implies that the perception of visual space is determined by the probability distribution of the possible real world sources of retinal images

rotation and illumination invariant interleaved intensity order based local descriptor the region descriptors using local intensity ordering patterns have become more popular recent years for image matching due to its enhanced discriminative ability however the dimension of these descriptors increases rapidly with the slight increase in the number of local neighbors under consideration and becomes unreasonable for image matching due to time constraint in this paper we reduce the dimension of the descriptor and matching time significantly while keeping up the comparable performance by considering the number of neighboring sample points in an interleaved manner the proposed interleaved order based local descriptor iold considers the local neighbors of a pixel as a set of interleaved neighbors and constructs the descriptor over each set separately and finally combines them to produce a single pattern we extract the local ordering pattern to cope up with the illumination effect in an inherent rotation invariant manner the novelty lies with using multiple neighboring sets in an interleaved fashion we also explored the local intensity order pattern in a multisupport region scenario results are compared over three challenging and widely adopted image matching data sets with other prominent descriptors under various image transformations results based on experiments suggest that the proposed iold descriptor outperforms in terms of both improved matching performance and reduced matching time we also found that the amount of improvement is significant under complex illumination difference while showing more robustness toward noise

asymmetric pulse modeling for fri sampling we consider sampling and reconstruction of finite rate of innovation fri signals such as a train of pulses where the pulses have varying degrees of asymmetry we address the problem of asymmetry modeling starting from a given symmetric prototype we show that among the class of unitary operators that are linear and invariant to translation and scale the fractional hilbert frh operator is unique for parametrically modeling pulse asymmetry the frh operator is obtained by a trigonometric interpolation between the standard hilbert and identity operators where the interpolation weights are determined by the degree of asymmetry the frh operators are also steerable which allows for estimation of the asymmetry factors in addition to the delays and amplitudes using the high resolution spectral estimation techniques that are used for solving standard fri problems we also develop the discrete counterpart using discrete frh operators and show that all the desirable properties carry over smoothly to the discrete setting as well we derive closed form expressions for the cramer rao bounds and hammersley chapman robbins bound on the variances of the estimators for continuous and discrete parameters respectively experimental results show that the proposed estimators have variances that meet the lower bounds we demonstrate an application of the proposed discrete frh methodology on real electrocardiogram ecg signals in the presence of noise specifically we show how the asymmetry of qrs complexes in various channels of an ecg signal could be modeled accurately

design of steerable filters for feature detection using canny like criteria we propose a general approach for the design of 2d feature detectors from a class of steerable functions based on the optimization of a canny like criterion in contrast with previous computational designs our approach is truly 2d and provides filters that have closed form expressions it also yields operators that have a better orientation selectivity than the classical gradient or hessian based detectors we illustrate the method with the design of operators for edge and ridge detection we present some experimental results that demonstrate the performance improvement of these new feature detectors we propose computationally efficient local optimization algorithms for the estimation of feature orientation we also introduce the notion of shape adaptable feature detection and use it for the detection of image corners

a shearlet approach to edge analysis and detection it is well known that the wavelet transform provides a very effective framework for analysis of multiscale edges in this paper we propose a novel approach based on the shearlet transform a multiscale directional transform with a greater ability to localize distributed discontinuities such as edges indeed unlike traditional wavelets shearlets are theoretically optimal in representing images with edges and in particular have the ability to fully capture directional and other geometrical features numerical examples demonstrate that the shearlet approach is highly effective at detecting both the location and orientation of edges and outperforms methods based on wavelets as well as other standard methods furthermore the shearlet approach is useful to design simple and effective algorithms for the detection of corners and junctions

spatially adaptive estimation via fitted local likelihood techniques this paper offers a new technique for spatially adaptive estimation the local likelihood is exploited for nonparametric modeling of observations and estimated signals the approach is based on the assumption of a local homogeneity of the signal for every point there exists a neighborhood in which the signal can be well approximated by a constant the fitted local likelihood statistics are used for selection of an adaptive size and shape of this neighborhood the algorithm is developed for a quite general class of observations subject to the exponential distribution the estimated signal can be uni and multivariable we demonstrate a good performance of the new algorithm for image denoising and compare the new method versus the intersection of confidence interval ici technique that also exploits a selection of an adaptive neighborhood for estimation

wavelet based fingerprint image retrieval this paper presents a novel approach for personal identification based on a wavelet based fingerprint retrieval system which encompasses three image retrieval tasks namely feature extraction similarity measurement and feature indexing we propose the use of different types of wavelets for representing and describing the textural information presented in fingerprint images in a compact way for that purpose the feature vectors used to characterize the fingerprints are obtained by computing the mean and the standard deviation of the decomposed images in the wavelet domain these feature vectors are used both to retrieve the most similar fingerprints given a query image and their indexation is used to reduce the search spaces of candidate images the different types of wavelets used in our study include gabor wavelets tree structured wavelet decomposition using both orthogonal and bi orthogonal filter banks as well as the steerable wavelets to evaluate the retrieval accuracy of the proposed approach a total number of eight different data sets were considered we also took into account different combinations of the above wavelets with six similarity measures the results show that the gabor wavelets combined with the square chord similarity measure achieves the best retrieval effectiveness

linear shape from shading in many situations the reflectance function of a surface is approximately linear and there is an effielent closed form solution to the shape from shading problem when boundary conditions e g edges singular points are not available good estimates of shape may still be extracted by using the assumption of general viewing position an improved method for estimating the illuminant direction is also presented

natural image statistics for natural image segmentation we integrate a model for filter response statistics of natural images into a variational framework for image segmentation incorporated in a sound probabilistic distance measure the model drives level sets toward meaningful segmentations of complex textures and natural scenes despite its enhanced descriptive power our approach preserves the efficiency of level set based segmentation since each region comprises two model parameters only analyzing thousands of natural images we select suitable filter banks validate the statistical basis of our model and demonstrate that it outperforms variational segmentation methods using second order statistics

differentiation of discrete multidimensional signals we describe the design of finite size linear phase separable kernels for differentiation of discrete multidimensional signals the problem is formulated as an optimization of the rotation invariance of the gradient operator which results in a simultaneous constraint on a set of one dimensional low pass prefilter and differentiator filters up to the desired order we also develop extensions of this formulation to both higher dimensions and higher order directional derivatives we develop a numerical procedure for optimizing the constraint and demonstrate its use in constructing a set of example filters the resulting filters are significantly more accurate than those commonly used in the image and multidimensional signal processing literature

a novel active contour model for medical images via the hessian matrix and eigenvalues abstract this paper presents a new level set formulation for active contour models acm we propose the idea of integrating the eigenvalue information of hessian matrix into the level set function by this new level set function the principal curvature information of images is used to enhance the ability of segmenting boundary regions the advantages of our model are as follows firstly the interior and exterior object boundaries can be segmented with the initial contour being anywhere in the input image secondly this method can work with heterogeneous images thirdly the proposed model can produce smooth and right boundaries of objects having vital importance in medical operations extensive experiments demonstrate that the proposed model can obtain better segmentation results