Multi-Orientation Estimation: Selectivity and Localization
|
|
- Clyde Cobb
- 5 years ago
- Views:
Transcription
1 in: M. van Ginkel, P.W. Verbeek, and L.J. van Vliet, Multi-orientation estimation: Selectivity and localization, in: H.E. Bal, H. Corporaal, P.P. Jonker, J.F.M. Tonino (eds.), ASCI 97, Proc. 3rd Annual Conference of the Advanced School for Computing and Imaging (Heijen, NL, June 2-4), ASCI, Delft, 1997, Multi-Orientation Estimation: Selectivity and Localization M. van Ginkel, P.W. Verbeek, L.J. van Vliet Pattern Recognition Group Faculty of Applied Physics, Delft University of Technology, Lorentzweg 1, 2628 CJ Delft, The Netherlands, Keywords: lines, edges, orientation, teture Abstract Filtering of an image with rotated versions of an orientation selective filter yields a set of images which can be stacked to form an orientation space. Orientation space provides a means of analyzing overlapping and touching anisotropic tetures. A set of rotated k th order directional derivatives yields a discrete orientation space, which allows interpolation. Net we apply a deconvolution scheme that results in improved orientation selectivity. This scheme allows decomposition of noisy multi-orientation patterns. 1 Introduction Images are composed of various types of structures, such as lines, edges and tetures. These in turn can be characterized by various properties, the most important being intensity, scale and orientation. In this paper we focus on improved orientation selectivity for multi-orientation estimation. An image can often locally be modeled as a weighted sum of translation invariant patterns or (paintbrush) strokes. Each stroke has a onedimensional intensity profile and a typical orientation: the profile orientation across the stroke, see figure 1. We define directions as angles in the interval (0, 2π), thus making a distinction between an arrow pointing to the left or right. When we refer to orientation (0, π) we make no such distinction. For many applications a reliable estimator for a single locally dominant orientation suffices, Kass & Witkin [4], Haglund [3] and Van Vliet & Verbeek [9]. There are a number of applications for Figure 1: An oriented pattern which an estimate of the local orientation will not be sufficient. Tetures can often be characterized by a number of overlapping patterns with a different orientation. At boundaries between two single orientation regions, an estimator for a locally dominant orientation will give the wrong answer. In such cases a more advanced analysis is required. For another important property, scale, an etensive framework has been created over the last decade, initiated by Witkin [12] and Koenderink [6]. In the scale space paradigm, scale is eplicitly dealt with by embedding an image in a new image with one etra dimension, the scale dimension. This image is generated from the original by iteratively applying an isotropic (e.g. Gaussian) filter. In order to deal with anisotropy an image or scale space can be replaced by a stack of directionally filtered versions of the original. The resulting orientation space or orientation+scale space is periodic along the orientation ais with a period of π. Orientation space and scale space differ in two aspects; the original image is part of scale space, but not of orientation space; scale space has a hierarchical structure, whereas orientation space does not. As shown in Van Vliet and Verbeek [10], such an orientation space can also be used to per- 99
2 form a segmentation of overlapping objects. The main objective of this paper is to establish how to generate an orientation space representation of an image with a high orientation selectivity. Interesting related work on orientation estimation and junction classification has been done by Andersson [1] and Michaelis [7]. 2 Constructing an orientation space Orientation space is created by applying rotated versions of some orientation selective filter to the image. Initially the only constraint on the choice of filter is that it allows proper sampling of orientation space. The filters we use are all two dimensional and operate in the (, y) plane. Yet, most of the ensuing (Fourier) analysis takes place in the orientation dimension and will be onedimensional. 2.1 Directional derivatives We will consider the possibility of using a Gaussian directional derivative of order k as orientation selective filter. In section 2.3 we will use the results of this section to derive a filter with a better orientation selectivity. The Gaussian is only used for regularization of the derivative operators. Therefore the σ, which specifies the width of the Gaussian, is not of interest in what follows. The Gaussian derivative filter of order k in the direction is given by: G (k) (, y, ) = (cos + sin y )k ep( 2 + y 2 2σ 2 ) (1) Net consider an image I(, y) that consists of a sum of strokes f i, with corresponding orientations θ i. Thus I(, y) can locally be described by: I(, y) = i f i ( cos θ i + y sin θ i ) (2) The response of the k th order directional Gaussian derivative to the stroke f i when = θ i, is denoted by f k,ma i. The response of the directional derivative to the image I(, y) can then be written as: I (k) (, y, ) = G (k) (, y, ) (,y) I(, y) = i cos k ( θ i )f k,ma i (, y) (3) Where (,y) denotes convolution in the (, y) plane. This can be rewritten as: I (k) (, y, ) = i cos k () f k,ma i (, y)δ( θ i ) (4) Convolution in the -dimension is denoted by (). The last equation has a simple interpretation. In orientation space each stroke gives an impulse in the -dimension convolved with the following kernel: a (k) () = cos k (5) A property of this kernel is that the orientation selectivity increases with increasing order k. This means that a specific selectivity can be chosen for a particular application. A less desirable property is that not only the orientation selectivity changes with the order k, but the radial frequency sensitivity as well. Therefore we will follow Knutsson [5] and decompose the filter into an angular and a radial part. The Fourier transform of a directional Gaussian derivative using polar coordinates ω and θ is given by: G (k) (ω, θ, ) = (jω) k ep( 1 2 σ2 ω 2 ) cos k ( θ) (6) Although separable, both the angular response, cos k ( θ), and the radial response depend on k. Since we wish to vary only the orientation selectivity of the filter, we will keep the radial response fied at k = 2 by multiplying G (k) by ω k+2. This is a rather arbitrary choice, but will suffice for the purpose of this paper. The actual choice of radial function depends on the application. The Fourier transform of the resulting filter A is as follows: Ã (k) (ω, θ, ) = ω k+2 G(k) (ω, θ, ) = j k ω 2 ep( 1 2 σ2 ω 2 ) cos k ( θ) (7) 2.2 Sampling the orientation space In practice only a finite number of filters can be used. Therefore the orientation space must be sampled. In this section we derive how many filters are needed for a given order k to allow reconstruction. Given a certain filter, Freeman & Adelson [2] and Perona [8] address the problem of creating a set of filters that allows interpolation and minimizes the error between this set and the filter response obtained by rotation of the original filter. They also present the filter constraints required to allow interpolation. The following analysis is essentially identical to the derivations in the references given above and shows that our filter allows interpolation. The Fourier coeffients of cos() are zero, ecept for ω = ±1. Repeated convolution (k 1 times) of 100
3 the Fourier series representation of cos by itself, yields the Fourier coefficients ã (k) (ω ) of a (k) (), which are then given by Pascal s triangle: ã (k) (ω ) = ( ) 1 k k k odd, ω odd k (k + ω ) k even, ω even k(8) 0 elsewhere The coefficients are depicted in figure 2a for order 10. The filter is bandlimited and the coefficients that correspond to frequencies higher than kω are zero. The Nyquist sampling theorem states that we need 2k + 1 samples on the interval (0, 2π) to allow reconstruction, or equivalently 2k + 1 filters. The number of filters needed can be reduced to k + 1, by noting the following symmetry: a (k) ( + π) = ( 1) k a (k) () (9) 2.3 Improving the angular resolution The kernel a (k) () can be well approimated by a Gaussian shape with a variance that decreases as 1 k. Therefore the peak width decreases 1 as k. Since the kernels are periodic the notion of peak width only makes sense on the interval ( π 2, π 2 ). The 1/ k behaviour is disappointing, because one would wish that doubling the number of filters would halve the peak width. We will now show how this can be achieved. Equation 4 shows that each stroke gives an impulse convolved with the kernel a (k) (). To improve the orientation selectivity we can simply deconvolve the resulting signal, because the blurring kernel is eactly known. As equation 8 shows, some Fourier series coefficients are zero, meaning that only a partial deconvolution can be performed. For order k the amplitude of the Fourier coefficients can be flattened (deconvolved) to: 1 k odd, ω odd k b(k) (ω ) = k+1 k even, ω even k 0 elsewhere (10) In figure 2b the coefficients b (k) (ω ) are depicted for k = 10. The resulting response in orientation space is given by: b (k) () = sin((k + 1)) sin (11) In figure 2c the kernels a (k) () and b (k) () are depicted for k = 10. It is clear that the peak width is indeed reduced, but at the cost of adding some side lobes. These lobes will cause problems if the amplitude of the different strokes differs too much. The first zerocrossing of b (k) () lies at π/(k + 1). This shows that the peak width indeed decreases approimately as 1/k. The response along the ais in orientation space to a stroke is equivalent to the angular response of the filter that is used. So instead of performing the deconvolution, it is also possible to use a filter that uses equation 11 as its angular response. The Fourier transform of the filter created by replacing the angular part of equation 7 by b (k) () is given by: B (k) (ω, θ, ) = j k ω 2 ep( 1 2 σ2 ω 2 sin((k + 1)( θ)) ) sin( θ) (12) In figure 3 the A and B filters are shown in the frequency domain as well as in the spatial domain. 2.4 Measuring the rms orientation intensity The filters A and B do not measure the rms orientation intensity, but respond to either locally symmetric or antisymmetric structures. For teture analysis we do not wish to discriminate between the two. The rms intensity is measured by the following procedure; First the filter result is squared, yielding an estimate of the local orientation energy, but with an image containing high frequency signals superimposed on it. By adding the square of its quadrature counterpart [5] the high frequency signals will cancel. Creating a quadrature counterpart for filter B results in a filter that is discontinuous for even k and nondifferentiable for odd k. Therefore it is not possible to use a quadrature filter approach. Instead the image will be smoothed by a Gaussian to remove the high frequency signals. The square root is taken of the smoothed image to yield the final estimate of the rms orientation intensity. Any interpolation must take place before the squaring operation, otherwise the necessary conditions to allow interpolation will no longer hold. All these operations take place in the (, y) planes of orientation space. 3 Single Peak detection In order to perform an analysis in orientation space, it is necessary to detect peaks in the dimension at every (, y) position. We assume that the peaks are non-overlapping and in the eperiments we have ensured this by choosing a sufficiently high order k. The shape of the peaks is 101
4 ω ω -0.2 π π (a) (b) 2 (c) 2 Figure 2: a.) Fourier coefficients ã (10) (ω ). The numbers show relative amplitude. b.) Fourier coefficients b(10) (ω ). c.) The kernels a (10) (), dashed, and b (10) (), solid line (a) A (9) (, y) (b) Ã(9) (ω, ω y) (c) B (9) (, y) (d) B (9) (ω, ω y) Figure 3: Spatial and frequency domain versions of the filters A and B known and given by b k (). Since we are not interested in the side lobes of this function, a simpler Gaussian shaped function will be used in the fit procedure. Because a fit has to be performed for every (, y) position, traditional iterative methods such as Levenberg-Marquardt require too much computation time. In this section we describe a fast alternative for the fit procedure. The σ of the Gaussian can be determined in advance since the width of the peaks is known. The Gaussian was fitted to a peak for order 9 yielding σ = radians. The fitted Gaussian coincides almost perfectly with the main lobe. Since the peak width decreases approimately as 1/(k + 1), the following relation gives the σ for any k: σ = (radians) (13) k + 1 To estimate the location and amplitude of the Gaussian the following error measure will be minimized: f Ag(ψ) 2 (14) Where f is the orientation response ( dimension), g a Gaussian positioned at ψ with amplitude factor A and the σ as given by equation 13. The inner product between functions m and n is defined as: (a) Figure 4: a.) Signal and inner product between shifted copies of the template and the signal (Gaussian filtering). b.) Signal and the detected peak. m() n() = (b) m()n()d (15) The integral is performed over one period of. Equation 14 can be written as: f 2 + Ag 2 2Af g (16) Only the last term is a function of ψ. Minimizing equation 14 with respect to ψ is equivalent to finding ma ψ (f g). Due to the symmetry of the Gaussian, convolving f with the Gaussian kernel yields f g for every value of ψ. The position of the maimum of the resulting signal is the ψ that maimizes f g. A typical f is shown in figure 4a. In the same figure f g is shown for every ψ. The net step is to fit the A parameter. Setting the deritivative of equation 16 with respect to A to zero and solving the resulting equation yields: 102
5 y orientation space y C slice CD slice AB A original image D (a) y A B C B (b) D (c) k=19 (d) (e) (f) (g) (h) (i) Figure 5: a) Schematic representation of orientation space. An, slice as indicated is used in (d-f). b) A noise free image containing two superimposed patterns. The intersection of the slice with the (, y) plane is indicated by line AB. The intersection with a second slice is indicated by line CD. c) Noisy version of (b). d) The slice indicated in (a) for the rms orientation intensity measured by filter A in image (b). e) Same as (d), but using filter B. f) Same as (e), but applied to image (c). (g-i) Same as (d-f) but using the other slice as indicated by line CD. Note the structures net to (f) indicating where in the slice each pattern gives a response. A = f g g 2 (17) In figure 4b the peak as detected by the algorithm, Ag, is shown plotted over the original signal f. To improve the fit it is necessary to remove the background offset before using the algorithm as discussed above. Since each parameter is optimized individually, there is no guarantee that a global optimum will be found. Maimizing f g will first yield the peaks with the largest volume. Such a peak is either a true peak or may consist of several overlapping peaks (This does not happen in the test images). Therefore a quality measure is needed to assess whether our assumption of isolated peaks holds. We will use the following quality measure: Q = f W g f W (18) Where f W is a windowed version of f with the window s center at the position of the detected peak. This is to prevent other peaks from influencing the quality measure. The window should be about as wide as the Gaussian, for eample 4σ. Q lies in the interval [ 1, 1], with 1 indicating a 103 high quality. By subtracting the peak that was found from f, a second peak can now be detected by repeating the procedure. Since Gaussian filtering can be implemented very efficiently [13], the entire procedure is very fast. 4 Eperiments In figure 5 we show the results of applying the filters proposed in this paper to an image consisting of two overlapping patterns. Each of the patterns has been created by calculating the Euclidean distance to some point in the image and computing the sine of the result. The amplitude of each individual pattern is one. The second test image was created by adding Gaussian noise with a variance of one to the first test image (SNR=6dB). We have computed the orientation space for both images using both filter A and filter B for k = 19. We show two slices of the orientation space for the following combinations: noise free image and filter A, noise free image and filter B, noisy image and filter B. The slices show the rms orientation intensity, using σ = 5 for the smoothing. The resulting slices show that the orientation selectivity is indeed improved. It is also clear that the scheme is quite robust with respect to noise. Furthermore, close eamination of the slices shows that there is indeed sub-piel infor-
6 (a) (b) (c) (d) Figure 6: a) Three touching patterns: left 51 degrees with respect to ais, middle 36, right 21. Image size 224 by 224 piels, diameter of circle 128 piels. Gaussian noise with variance 0.2 was added to the image. b), slice of orientation space halfway along the y ais. Parameters used: k = 27, σ filter = 1.5, σ smoothing = 5. c) Magnitude of second detected peak divided by magnitude of first detected peak for each, y position. d) Detected boundaries superimposed on (a). mation present along the orientation ais, and it is therefore possible to interpolate in (the original) orientation space to get improved accuracy. The net eperiment concerns touching orientation fields. A test image was created containing a circular region in the middle and two regions around it separated by a second vertical boundary in the middle of the image. The circular region has an orientation of 36 degrees with respect to the ais, the left region 51 degrees, and the right region 21 degrees. The amplitude of the patterns is again one, and Gaussian noise with a variance of 0.2 (SNR=13dB) was added. Using the method described in section 3, the dominant peak in orientation space was found for each, y position. In figure 6c the magnitude of the second peak divided by the magnitude of the first is shown. When two orientation fields are present, as at the boundaries of the regions, the two peaks found will have approimately the same magnitude and hence the quotient will approach one. Indeed, the quotient is maimal at the boundaries as the results in figure 6c show. Using a watershed [11] based technique, the ridges in figure 6c are detected and shown superimposed on the original image in figure 6d. They coincide almost perfectly with the actual boundaries. Since the parameters were chosen to ensure that there are no overlapping peaks in orientation space, it was not necessary to use the quality measure as defined in section 3. The location of the maima is not determined with subpiel accuracy in the current implementation. Subtracting a slightly misaligned peak results in larger residual flanks. The orientation of the center region in figure 6a corresponds eactly with a piel position, whereas the other patterns do not. As can be seen in figure 6c, the residual peaks are indeed larger for the other two regions. The location of the boundaries found in the second eample is very accurate. However, as the difference in angle between the different regions decreases, a higher order filter will be required to distinguish between the regions. But, as the 104
7 narrow wide [2] W.T. Freeman and E.H. Adelson, The Design and Use of Steerable filters, IEEE transactions on Pattern Analysis and Machine Intelligence, vol. 13, no. 9, September 1991, pp [3] L. Haglund, Adaptive Multidimensional Filtering, PhD thesis, Linköping University, Sweden, Figure 7: Localization depends on orientation with respect to the boundary filter becomes narrower in the frequency domain, the filter will become wider in the spatial domain. This will have a corresponding adverse effect on localization. Due to the shape of the filter as shown in figure 3c, localization depends on the orientation of the boundary with respect to the oriented pattern, see figure 7. Fortunately, the situation where the pattern(s) and the boundary have a similar orientation (allowing good localization) occurs most frequently in real world images. 5 Conclusions We have shown that the peak width of the response in orientation space can be made to depend linearly on the number of filters k. This occurs at the cost of introducing some side lobes. Simple eperiments with this scheme show that it works well for noisy images. The method used to detect the boundaries is still rather crude and more advanced techniques to analyze orientation space must be developed. Real orientation fields may not only differ in orientation, but also in amplitude and frequency content. Furthermore the interplay between filter order, frequency content, attainable orientation selectivity and localization and its consequences must be investigated. Acknowledgements This work was partially supported by the Rolling Grants program 94RG12 of the Netherlands Organization for Fundamental Research of Matter (FOM) and by the Royal Dutch Academy of Sciences (KNAW). References [1] M. Andersson, Controllable Multidimensional Filters and Models in Low Level Computer Vision, PhD thesis, Linköping University, Sweden, [4] M. Kass and A. Witkin, Analyzing Oriented Patterns, Computer Vision, Graphics and Image Processing, vol. 37, 1987, pp [5] H. Knutsson, Filtering and Reconstruction in Image Processing, PhD thesis, Linköping University, Sweden, [6] J.J. Koenderink, The Structure of Images, Biological Cybernetics, vol. 50, 1984, pp [7] M. Michaelis and G. Sommer, Junction Classification by multiple orientation detection, in: Jan-Olof Eklundh (ed.) ECCV 94, Third European Conference on Computer Vision, Stockholm, Sweden, May 1994, pp [8] P. Perona, Deformable Kernels for Early Vision, IEEE transactions on Pattern Analysis and Machine Intelligence, vol. 17, no. 5, May 1995, pp [9] L.J. van Vliet and P.W. Verbeek, Estimators for Orientation and Anisotropy in Digitized Images, in: J. van Katwijk, J.J. Gerbrands, M.R. van Steen, J.F.M. Tonino (eds.), ASCI 95, Proc. First Annual Conference of the Advanced School for Computing and Imaging (Heijen, NL, May 16-18), ASCI, Delft, 1995, pp [10] L.J. van Vliet and P.W. Verbeek, Segmentation of overlapping objects, in: L.J. van Vliet, I.T. Young (eds.), Abstracts of the ASCI Imaging Workshop 1995, Venray, The Netherlands, October 1995, pp [11] B.J.H. Verwer, L.J. van Vliet and P.W. Verbeek, Binary and Grey-value skeletons: Metrics and Algorithms, International Journal of Pattern Recognition and Artificial Intelligence, vol. 7, no. 5, 1993, pp [12] A. Witkin, Scale space filtering, Proc. Int. Joint Conf. on Artif. Intell., Karlsruhe, Germany, 1983, pp [13] I.T. Young, L.J. van Vliet, Recursive implementation of the Gaussian filter, Signal Processing, vol. 44, no. 2, 1995,
3D-Orientation Space; Filters and Sampling
3D-Orientation Space; Filters and Sampling Frank G.A. Faas and Lucas J. van Vliet Pattern Recognition Group, Delft University of Technology, Lorentzweg 1, 2628CJ, Delft, The Netherlands, {faas,lucas}@ph.tn.tudelft.nl
More informationRobust Curve Detection using a Radon Transform in Orientation Space
Robust Curve Detection using a Radon Transform in Orientation Space M. van Ginkel 1, M.A. Kraaijveld 2, L.J. van Vliet 1, E.P. Reding 2, P.W. Verbeek 1, and H.J. Lammers 2 1 Pattern Recognition Group,
More informationHow and what do we see? Segmentation and Grouping. Fundamental Problems. Polyhedral objects. Reducing the combinatorics of pose estimation
Segmentation and Grouping Fundamental Problems ' Focus of attention, or grouping ' What subsets of piels do we consider as possible objects? ' All connected subsets? ' Representation ' How do we model
More informationNormalized averaging using adaptive applicability functions with applications in image reconstruction from sparsely and randomly sampled data
Normalized averaging using adaptive applicability functions with applications in image reconstruction from sparsely and randomly sampled data Tuan Q. Pham, Lucas J. van Vliet Pattern Recognition Group,
More informationWeichuan Yu, Kostas Daniilidis, Gerald Sommer. Christian Albrechts University. Preusserstrasse 1-9, D Kiel, Germany
Rotated Wedge Averaging Method for Junction Characterization Weichuan Yu, Kostas Daniilidis, Gerald Sommer Institute of Computer Science Christian Albrechts University Preusserstrasse -9, D-45 Kiel, Germany
More informationINVARIANT CORNER DETECTION USING STEERABLE FILTERS AND HARRIS ALGORITHM
INVARIANT CORNER DETECTION USING STEERABLE FILTERS AND HARRIS ALGORITHM ABSTRACT Mahesh 1 and Dr.M.V.Subramanyam 2 1 Research scholar, Department of ECE, MITS, Madanapalle, AP, India vka4mahesh@gmail.com
More information(10) Image Segmentation
(0) Image Segmentation - Image analysis Low-level image processing: inputs and outputs are all images Mid-/High-level image processing: inputs are images; outputs are information or attributes of the images
More informationKey properties of local features
Key properties of local features Locality, robust against occlusions Must be highly distinctive, a good feature should allow for correct object identification with low probability of mismatch Easy to etract
More informationFast Anisotropic Gauss Filtering
Fast Anisotropic Gauss Filtering Jan-Mark Geusebroek, Arnold W. M. Smeulders, and Joost van de Weijer Intelligent Sensory Information Systems, Department of Computer Science, University of Amsterdam, Kruislaan
More informationFourier Transform and Texture Filtering
Fourier Transform and Texture Filtering Lucas J. van Vliet www.ph.tn.tudelft.nl/~lucas Image Analysis Paradigm scene Image formation sensor pre-processing Image enhancement Image restoration Texture filtering
More informationEffects of multi-scale velocity heterogeneities on wave-equation migration Yong Ma and Paul Sava, Center for Wave Phenomena, Colorado School of Mines
Effects of multi-scale velocity heterogeneities on wave-equation migration Yong Ma and Paul Sava, Center for Wave Phenomena, Colorado School of Mines SUMMARY Velocity models used for wavefield-based seismic
More informationAnnouncements. Edge Detection. An Isotropic Gaussian. Filters are templates. Assignment 2 on tracking due this Friday Midterm: Tuesday, May 3.
Announcements Edge Detection Introduction to Computer Vision CSE 152 Lecture 9 Assignment 2 on tracking due this Friday Midterm: Tuesday, May 3. Reading from textbook An Isotropic Gaussian The picture
More informationImage Sampling and Quantisation
Image Sampling and Quantisation Introduction to Signal and Image Processing Prof. Dr. Philippe Cattin MIAC, University of Basel 1 of 46 22.02.2016 09:17 Contents Contents 1 Motivation 2 Sampling Introduction
More informationImage Sampling & Quantisation
Image Sampling & Quantisation Biomedical Image Analysis Prof. Dr. Philippe Cattin MIAC, University of Basel Contents 1 Motivation 2 Sampling Introduction and Motivation Sampling Example Quantisation Example
More informationCircular Analysis based Line Detection Filters for Watermark Extraction in X-ray images of Etchings
Circular Analysis based Line Detection Filters for Watermark Extraction in X-ray images of Etchings M. van Staalduinen, J. C. A. van der Lubbe, E. Backer {M.vanStaalduinen, J.C.A.vanderLubbe, E.Backer}@ewi.tudelft.nl
More informationThe Generalised Radon Transform: Sampling and Memory Considerations
The Generalised Radon Transform: Sampling and Memory Considerations C.L. Luengo Hendriks, M. van Ginkel, P.W. Verbeek and L.J. van Vliet Pattern Recognition Group, Delft University of Technology, Lorentzweg
More informationPrecalculus Notes Unit 1 Day 1
Precalculus Notes Unit Day Rules For Domain: When the domain is not specified, it consists of (all real numbers) for which the corresponding values in the range are also real numbers.. If is in the numerator
More informationu u 1 u (c) Distributive property of multiplication over subtraction
ADDITIONAL ANSWERS 89 Additional Answers Eercises P.. ; All real numbers less than or equal to 4 0 4 6. ; All real numbers greater than or equal to and less than 4 0 4 6 7. ; All real numbers less than
More informationMulti-dimensional Image Analysis
Multi-dimensional Image Analysis Lucas J. van Vliet www.ph.tn.tudelft.nl/~lucas Image Analysis Paradigm scene Image formation sensor pre-processing Image enhancement Image restoration Texture filtering
More informationTexture Segmentation Using Multichannel Gabor Filtering
IOSR Journal of Electronics and Communication Engineering (IOSRJECE) ISSN : 2278-2834 Volume 2, Issue 6 (Sep-Oct 2012), PP 22-26 Texture Segmentation Using Multichannel Gabor Filtering M. Sivalingamaiah
More informationUsing a multipoint interferometer to measure the orbital angular momentum of light
CHAPTER 3 Using a multipoint interferometer to measure the orbital angular momentum of light Recently it was shown that the orbital angular momentum of light can be measured using a multipoint interferometer,
More informationTHE preceding chapters were all devoted to the analysis of images and signals which
Chapter 5 Segmentation of Color, Texture, and Orientation Images THE preceding chapters were all devoted to the analysis of images and signals which take values in IR. It is often necessary, however, to
More informationHOUGH TRANSFORM CS 6350 C V
HOUGH TRANSFORM CS 6350 C V HOUGH TRANSFORM The problem: Given a set of points in 2-D, find if a sub-set of these points, fall on a LINE. Hough Transform One powerful global method for detecting edges
More informationLecture Image Enhancement and Spatial Filtering
Lecture Image Enhancement and Spatial Filtering Harvey Rhody Chester F. Carlson Center for Imaging Science Rochester Institute of Technology rhody@cis.rit.edu September 29, 2005 Abstract Applications of
More informationFiltering and Edge Detection. Computer Vision I. CSE252A Lecture 10. Announcement
Filtering and Edge Detection CSE252A Lecture 10 Announcement HW1graded, will be released later today HW2 assigned, due Wed. Nov. 7 1 Image formation: Color Channel k " $ $ # $ I r I g I b % " ' $ ' = (
More informationEdge and local feature detection - 2. Importance of edge detection in computer vision
Edge and local feature detection Gradient based edge detection Edge detection by function fitting Second derivative edge detectors Edge linking and the construction of the chain graph Edge and local feature
More informationSegmentation-free estimation of length distributions using sieves and RIA morphology
Segmentation-free estimation of length distributions using sieves and RIA morphology Cris L. Luengo Hendriks and Lucas J. van Vliet Pattern Recognition Group, Delft University of Technology, Lorentzweg
More informationExtracting Wavefront Error From Shack-Hartmann Images Using Spatial Demodulation
Etracting Wavefront Error From Shack-Hartmann Images Using Spatial Demodulation Edwin J. Sarver, PhD; Jim Schwiegerling, PhD; Raymond A. Applegate, OD, PhD ABSTRACT PURPOSE: To determine whether the spatial
More informationTexture. Outline. Image representations: spatial and frequency Fourier transform Frequency filtering Oriented pyramids Texture representation
Texture Outline Image representations: spatial and frequency Fourier transform Frequency filtering Oriented pyramids Texture representation 1 Image Representation The standard basis for images is the set
More informationAliasing. Can t draw smooth lines on discrete raster device get staircased lines ( jaggies ):
(Anti)Aliasing and Image Manipulation for (y = 0; y < Size; y++) { for (x = 0; x < Size; x++) { Image[x][y] = 7 + 8 * sin((sqr(x Size) + SQR(y Size)) / 3.0); } } // Size = Size / ; Aliasing Can t draw
More informationThe Steerable Pyramid: A Flexible Architecture for Multi-Scale Derivative Computation
MITSUBISHI ELECTRIC RESEARCH LABORATORIES http://www.merl.com The Steerable Pyramid: A Flexible Architecture for Multi-Scale Derivative Computation Eero P. Simoncelli, William T. Freeman TR95-15 December
More informationMultimedia Computing: Algorithms, Systems, and Applications: Edge Detection
Multimedia Computing: Algorithms, Systems, and Applications: Edge Detection By Dr. Yu Cao Department of Computer Science The University of Massachusetts Lowell Lowell, MA 01854, USA Part of the slides
More informationFractional Discrimination for Texture Image Segmentation
Fractional Discrimination for Texture Image Segmentation Author You, Jia, Sattar, Abdul Published 1997 Conference Title IEEE 1997 International Conference on Image Processing, Proceedings of the DOI https://doi.org/10.1109/icip.1997.647743
More informationCoarse-to-fine image registration
Today we will look at a few important topics in scale space in computer vision, in particular, coarseto-fine approaches, and the SIFT feature descriptor. I will present only the main ideas here to give
More informationconvolution shift invariant linear system Fourier Transform Aliasing and sampling scale representation edge detection corner detection
COS 429: COMPUTER VISON Linear Filters and Edge Detection convolution shift invariant linear system Fourier Transform Aliasing and sampling scale representation edge detection corner detection Reading:
More informationCHAPTER 8 COMPOUND CHARACTER RECOGNITION USING VARIOUS MODELS
CHAPTER 8 COMPOUND CHARACTER RECOGNITION USING VARIOUS MODELS 8.1 Introduction The recognition systems developed so far were for simple characters comprising of consonants and vowels. But there is one
More informationRobust Lip Contour Extraction using Separability of Multi-Dimensional Distributions
Robust Lip Contour Extraction using Separability of Multi-Dimensional Distributions Tomokazu Wakasugi, Masahide Nishiura and Kazuhiro Fukui Corporate Research and Development Center, Toshiba Corporation
More informationAlgorithms for Recognition of Low Quality Iris Images. Li Peng Xie University of Ottawa
Algorithms for Recognition of Low Quality Iris Images Li Peng Xie University of Ottawa Overview Iris Recognition Eyelash detection Accurate circular localization Covariance feature with LDA Fourier magnitude
More informationAnnouncements, schedule. Lecture 8: Fitting. Weighted graph representation. Outline. Segmentation by Graph Cuts. Images as graphs
Announcements, schedule Lecture 8: Fitting Tuesday, Sept 25 Grad student etensions Due of term Data sets, suggestions Reminder: Midterm Tuesday 10/9 Problem set 2 out Thursday, due 10/11 Outline Review
More information3D Object Scanning to Support Computer-Aided Conceptual Design
ABSTRACT 3D Object Scanning to Support Computer-Aided Conceptual Design J.S.M. Vergeest and I. Horváth Delft University of Technology Faculty of Design, Engineering and Production Jaffalaan 9, NL-2628
More informationAnno accademico 2006/2007. Davide Migliore
Robotica Anno accademico 6/7 Davide Migliore migliore@elet.polimi.it Today What is a feature? Some useful information The world of features: Detectors Edges detection Corners/Points detection Descriptors?!?!?
More informationImage Enhancement Techniques for Fingerprint Identification
March 2013 1 Image Enhancement Techniques for Fingerprint Identification Pankaj Deshmukh, Siraj Pathan, Riyaz Pathan Abstract The aim of this paper is to propose a new method in fingerprint enhancement
More informationused to describe all aspects of imaging process: input scene f, imaging system O, and output image g.
SIMG-716 Linear Imaging Mathematics I 03 1 Representations of Functions 1.1 Functions Function: rule that assigns one numerical value (dependent variable) to each of a set of other numerical values (independent
More informationDigital Image Processing. Prof. P. K. Biswas. Department of Electronic & Electrical Communication Engineering
Digital Image Processing Prof. P. K. Biswas Department of Electronic & Electrical Communication Engineering Indian Institute of Technology, Kharagpur Lecture - 21 Image Enhancement Frequency Domain Processing
More informationELEC Dr Reji Mathew Electrical Engineering UNSW
ELEC 4622 Dr Reji Mathew Electrical Engineering UNSW Review of Motion Modelling and Estimation Introduction to Motion Modelling & Estimation Forward Motion Backward Motion Block Motion Estimation Motion
More informationComputer Vision. Fourier Transform. 20 January Copyright by NHL Hogeschool and Van de Loosdrecht Machine Vision BV All rights reserved
Van de Loosdrecht Machine Vision Computer Vision Fourier Transform 20 January 2017 Copyright 2001 2017 by NHL Hogeschool and Van de Loosdrecht Machine Vision BV All rights reserved j.van.de.loosdrecht@nhl.nl,
More informationA DH-parameter based condition for 3R orthogonal manipulators to have 4 distinct inverse kinematic solutions
Wenger P., Chablat D. et Baili M., A DH-parameter based condition for R orthogonal manipulators to have 4 distinct inverse kinematic solutions, Journal of Mechanical Design, Volume 17, pp. 150-155, Janvier
More informationSteerable Filters from Erlang Functions
Steerable Filters from Erlang Functions Anil A. Bharath Department of Biological & Medical Systems, Imperial College of Science, Technology & Medicine, Exhibition Road, London SW7 2BX a.bharath@ic.ac.uk
More informationImage features. Image Features
Image features Image features, such as edges and interest points, provide rich information on the image content. They correspond to local regions in the image and are fundamental in many applications in
More informationused to describe all aspects of imaging process: input scene f, imaging system O, and output image g.
SIMG-716 Linear Imaging Mathematics I 03 1 Representations of Functions 1.1 Functions Function: rule that assigns one numerical value (dependent variable) to each of a set of other numerical values (independent
More informationDisparity from Monogenic Phase
Disparity from Monogenic Phase Michael Felsberg Department of Electrical Engineering, Linköping University, SE-58183 Linköping, Sweden mfe@isy.liu.se, WWW home page: http://www.isy.liu.se/~mfe Abstract.
More informationEdge Detection. CS664 Computer Vision. 3. Edges. Several Causes of Edges. Detecting Edges. Finite Differences. The Gradient
Edge Detection CS664 Computer Vision. Edges Convert a gray or color image into set of curves Represented as binary image Capture properties of shapes Dan Huttenlocher Several Causes of Edges Sudden changes
More informationAdaptative Elimination of False Edges for First Order Detectors
Adaptative Elimination of False Edges for First Order Detectors Djemel ZIOU and Salvatore TABBONE D~partement de math~matiques et d'informatique, universit~ de Sherbrooke, Qc, Canada, J1K 2R1 Crin/Cnrs
More informationJournal of Asian Scientific Research FEATURES COMPOSITION FOR PROFICIENT AND REAL TIME RETRIEVAL IN CBIR SYSTEM. Tohid Sedghi
Journal of Asian Scientific Research, 013, 3(1):68-74 Journal of Asian Scientific Research journal homepage: http://aessweb.com/journal-detail.php?id=5003 FEATURES COMPOSTON FOR PROFCENT AND REAL TME RETREVAL
More informationImage Segmentation Techniques for Object-Based Coding
Image Techniques for Object-Based Coding Junaid Ahmed, Joseph Bosworth, and Scott T. Acton The Oklahoma Imaging Laboratory School of Electrical and Computer Engineering Oklahoma State University {ajunaid,bosworj,sacton}@okstate.edu
More informationSynthesis Imaging. Claire Chandler, Sanjay Bhatnagar NRAO/Socorro
Synthesis Imaging Claire Chandler, Sanjay Bhatnagar NRAO/Socorro Michelson Summer Workshop Caltech, July 24-28, 2006 Synthesis Imaging 2 Based on the van Cittert-Zernike theorem: The complex visibility
More informationLucas-Kanade Image Registration Using Camera Parameters
Lucas-Kanade Image Registration Using Camera Parameters Sunghyun Cho a, Hojin Cho a, Yu-Wing Tai b, Young Su Moon c, Junguk Cho c, Shihwa Lee c, and Seungyong Lee a a POSTECH, Pohang, Korea b KAIST, Daejeon,
More informationA Survey of Light Source Detection Methods
A Survey of Light Source Detection Methods Nathan Funk University of Alberta Mini-Project for CMPUT 603 November 30, 2003 Abstract This paper provides an overview of the most prominent techniques for light
More informationScanner Parameter Estimation Using Bilevel Scans of Star Charts
ICDAR, Seattle WA September Scanner Parameter Estimation Using Bilevel Scans of Star Charts Elisa H. Barney Smith Electrical and Computer Engineering Department Boise State University, Boise, Idaho 8375
More informationSchedule for Rest of Semester
Schedule for Rest of Semester Date Lecture Topic 11/20 24 Texture 11/27 25 Review of Statistics & Linear Algebra, Eigenvectors 11/29 26 Eigenvector expansions, Pattern Recognition 12/4 27 Cameras & calibration
More informationFourier transforms and convolution
Fourier transforms and convolution (without the agonizing pain) CS/CME/BioE/Biophys/BMI 279 Oct. 26, 2017 Ron Dror 1 Why do we care? Fourier transforms Outline Writing functions as sums of sinusoids The
More informationAnalysis of Planar Anisotropy of Fibre Systems by Using 2D Fourier Transform
Maroš Tunák, Aleš Linka Technical University in Liberec Faculty of Textile Engineering Department of Textile Materials Studentská 2, 461 17 Liberec 1, Czech Republic E-mail: maros.tunak@tul.cz ales.linka@tul.cz
More informationComputer Vision 2. SS 18 Dr. Benjamin Guthier Professur für Bildverarbeitung. Computer Vision 2 Dr. Benjamin Guthier
Computer Vision 2 SS 18 Dr. Benjamin Guthier Professur für Bildverarbeitung Computer Vision 2 Dr. Benjamin Guthier 1. IMAGE PROCESSING Computer Vision 2 Dr. Benjamin Guthier Content of this Chapter Non-linear
More informationComparative Study of Dual-Tree Complex Wavelet Transform and Double Density Complex Wavelet Transform for Image Denoising Using Wavelet-Domain
International Journal of Scientific and Research Publications, Volume 2, Issue 7, July 2012 1 Comparative Study of Dual-Tree Complex Wavelet Transform and Double Density Complex Wavelet Transform for Image
More informationOutlines. Medical Image Processing Using Transforms. 4. Transform in image space
Medical Image Processing Using Transforms Hongmei Zhu, Ph.D Department of Mathematics & Statistics York University hmzhu@yorku.ca Outlines Image Quality Gray value transforms Histogram processing Transforms
More informationIllumination-Robust Face Recognition based on Gabor Feature Face Intrinsic Identity PCA Model
Illumination-Robust Face Recognition based on Gabor Feature Face Intrinsic Identity PCA Model TAE IN SEOL*, SUN-TAE CHUNG*, SUNHO KI**, SEONGWON CHO**, YUN-KWANG HONG*** *School of Electronic Engineering
More informationFiltering Images. Contents
Image Processing and Data Visualization with MATLAB Filtering Images Hansrudi Noser June 8-9, 010 UZH, Multimedia and Robotics Summer School Noise Smoothing Filters Sigmoid Filters Gradient Filters Contents
More informationAccurate and Dense Wide-Baseline Stereo Matching Using SW-POC
Accurate and Dense Wide-Baseline Stereo Matching Using SW-POC Shuji Sakai, Koichi Ito, Takafumi Aoki Graduate School of Information Sciences, Tohoku University, Sendai, 980 8579, Japan Email: sakai@aoki.ecei.tohoku.ac.jp
More informationMultiresolution Image Processing
Multiresolution Image Processing 2 Processing and Analysis of Images at Multiple Scales What is Multiscale Decompostion? Why use Multiscale Processing? How to use Multiscale Processing? Related Concepts:
More informationDigital Image Processing (CS/ECE 545) Lecture 5: Edge Detection (Part 2) & Corner Detection
Digital Image Processing (CS/ECE 545) Lecture 5: Edge Detection (Part 2) & Corner Detection Prof Emmanuel Agu Computer Science Dept. Worcester Polytechnic Institute (WPI) Recall: Edge Detection Image processing
More informationSnakes reparameterization for noisy images segmentation and targets tracking
Snakes reparameterization for noisy images segmentation and targets tracking Idrissi Sidi Yassine, Samir Belfkih. Lycée Tawfik Elhakim Zawiya de Noaceur, route de Marrakech, Casablanca, maroc. Laboratoire
More informationTexture Segmentation Using Gabor Filters
TEXTURE SEGMENTATION USING GABOR FILTERS 1 Texture Segmentation Using Gabor Filters Khaled Hammouda Prof. Ed Jernigan University of Waterloo, Ontario, Canada Abstract Texture segmentation is the process
More informationFACE RECOGNITION USING INDEPENDENT COMPONENT
Chapter 5 FACE RECOGNITION USING INDEPENDENT COMPONENT ANALYSIS OF GABORJET (GABORJET-ICA) 5.1 INTRODUCTION PCA is probably the most widely used subspace projection technique for face recognition. A major
More informationComparison of Feature Detection and Matching Approaches: SIFT and SURF
GRD Journals- Global Research and Development Journal for Engineering Volume 2 Issue 4 March 2017 ISSN: 2455-5703 Comparison of Detection and Matching Approaches: SIFT and SURF Darshana Mistry PhD student
More informationTEXTURE ANALYSIS USING GABOR FILTERS
TEXTURE ANALYSIS USING GABOR FILTERS Texture Types Definition of Texture Texture types Synthetic Natural Stochastic < Prev Next > Texture Definition Texture: the regular repetition of an element or pattern
More informationAutomated Drill Design Software
Automated Drill Design Software Athulan Vijayaraghavan March 19, 2006 Abstract This section of the report discusses a tool which can create automated 3D CAD drill models based on geometric as well as manufacturing
More information4 Parametrization of closed curves and surfaces
4 Parametrization of closed curves and surfaces Parametrically deformable models give rise to the question of obtaining parametrical descriptions of given pixel or voxel based object contours or surfaces,
More informationTexture Segmentation by using Haar Wavelets and K-means Algorithm
Texture Segmentation by using Haar Wavelets and K-means Algorithm P. Ashok Babu Associate Professor, Narsimha Reddy Engineering College, Hyderabad, A.P., INDIA, ashokbabup2@gmail.com Dr. K. V. S. V. R.
More informationImplementing the Scale Invariant Feature Transform(SIFT) Method
Implementing the Scale Invariant Feature Transform(SIFT) Method YU MENG and Dr. Bernard Tiddeman(supervisor) Department of Computer Science University of St. Andrews yumeng@dcs.st-and.ac.uk Abstract The
More informationERROR RECOGNITION and IMAGE ANALYSIS
PREAMBLE TO ERROR RECOGNITION and IMAGE ANALYSIS 2 Why are these two topics in the same lecture? ERROR RECOGNITION and IMAGE ANALYSIS Ed Fomalont Error recognition is used to determine defects in the data
More informationBlob-Representation of Multidimensional Objects and Surfaces
Blob-Representation of Multidimensional Objects and Surfaces Edgar Garduño and Gabor T. Herman Department of Computer Science The Graduate Center City University of New York Presentation Outline Reconstruction
More informationSUMMARY: DISTINCTIVE IMAGE FEATURES FROM SCALE- INVARIANT KEYPOINTS
SUMMARY: DISTINCTIVE IMAGE FEATURES FROM SCALE- INVARIANT KEYPOINTS Cognitive Robotics Original: David G. Lowe, 004 Summary: Coen van Leeuwen, s1460919 Abstract: This article presents a method to extract
More informationObject Recognition Using Pictorial Structures. Daniel Huttenlocher Computer Science Department. In This Talk. Object recognition in computer vision
Object Recognition Using Pictorial Structures Daniel Huttenlocher Computer Science Department Joint work with Pedro Felzenszwalb, MIT AI Lab In This Talk Object recognition in computer vision Brief definition
More informationFundamentals of Digital Image Processing
\L\.6 Gw.i Fundamentals of Digital Image Processing A Practical Approach with Examples in Matlab Chris Solomon School of Physical Sciences, University of Kent, Canterbury, UK Toby Breckon School of Engineering,
More informationSimultaneous surface texture classification and illumination tilt angle prediction
Simultaneous surface texture classification and illumination tilt angle prediction X. Lladó, A. Oliver, M. Petrou, J. Freixenet, and J. Martí Computer Vision and Robotics Group - IIiA. University of Girona
More informationTEXTURE ANALYSIS USING GABOR FILTERS FIL
TEXTURE ANALYSIS USING GABOR FILTERS Texture Types Definition of Texture Texture types Synthetic ti Natural Stochastic < Prev Next > Texture Definition Texture: the regular repetition of an element or
More informationarxiv: v1 [cs.cv] 29 Sep 2007
Graph rigidity, Cyclic Belief Propagation and Point Pattern Matching Julian J. McAuley, Tibério S. Caetano and Marconi S. Barbosa October 27, 28 arxiv:7.43v [cs.cv] 29 Sep 27 Abstract A recent paper []
More informationAccurately measuring 2D position using a composed moiré grid pattern and DTFT
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 Accurately measuring 2D position using a composed moiré grid pattern and DTFT S. Van
More informationLecture 2: 2D Fourier transforms and applications
Lecture 2: 2D Fourier transforms and applications B14 Image Analysis Michaelmas 2017 Dr. M. Fallon Fourier transforms and spatial frequencies in 2D Definition and meaning The Convolution Theorem Applications
More informationQuantifying Three-Dimensional Deformations of Migrating Fibroblasts
45 Chapter 4 Quantifying Three-Dimensional Deformations of Migrating Fibroblasts This chapter presents the full-field displacements and tractions of 3T3 fibroblast cells during migration on polyacrylamide
More informationBinary Image Processing. Introduction to Computer Vision CSE 152 Lecture 5
Binary Image Processing CSE 152 Lecture 5 Announcements Homework 2 is due Apr 25, 11:59 PM Reading: Szeliski, Chapter 3 Image processing, Section 3.3 More neighborhood operators Binary System Summary 1.
More informationMixture Models and EM
Mixture Models and EM Goal: Introduction to probabilistic mixture models and the expectationmaximization (EM) algorithm. Motivation: simultaneous fitting of multiple model instances unsupervised clustering
More informationAnnouncements. Edges. Last Lecture. Gradients: Numerical Derivatives f(x) Edge Detection, Lines. Intro Computer Vision. CSE 152 Lecture 10
Announcements Assignment 2 due Tuesday, May 4. Edge Detection, Lines Midterm: Thursday, May 6. Introduction to Computer Vision CSE 152 Lecture 10 Edges Last Lecture 1. Object boundaries 2. Surface normal
More informationSegmentation and Grouping
Segmentation and Grouping How and what do we see? Fundamental Problems ' Focus of attention, or grouping ' What subsets of pixels do we consider as possible objects? ' All connected subsets? ' Representation
More informationA Robust Method for Circle / Ellipse Extraction Based Canny Edge Detection
International Journal of Research Studies in Science, Engineering and Technology Volume 2, Issue 5, May 2015, PP 49-57 ISSN 2349-4751 (Print) & ISSN 2349-476X (Online) A Robust Method for Circle / Ellipse
More informationExperiments with Edge Detection using One-dimensional Surface Fitting
Experiments with Edge Detection using One-dimensional Surface Fitting Gabor Terei, Jorge Luis Nunes e Silva Brito The Ohio State University, Department of Geodetic Science and Surveying 1958 Neil Avenue,
More informationBiomedical Image Analysis. Point, Edge and Line Detection
Biomedical Image Analysis Point, Edge and Line Detection Contents: Point and line detection Advanced edge detection: Canny Local/regional edge processing Global processing: Hough transform BMIA 15 V. Roth
More informationWEINER FILTER AND SUB-BLOCK DECOMPOSITION BASED IMAGE RESTORATION FOR MEDICAL APPLICATIONS
WEINER FILTER AND SUB-BLOCK DECOMPOSITION BASED IMAGE RESTORATION FOR MEDICAL APPLICATIONS ARIFA SULTANA 1 & KANDARPA KUMAR SARMA 2 1,2 Department of Electronics and Communication Engineering, Gauhati
More informationGrade 4 CCSS Pacing Guide: Math Expressions
Grade 4 CCSS Pacing Guide: Math Expressions Green : Areas of intensive focus, where students need fluent understanding and application of the core concepts (Approximately 70%) Blue : Rethinking and linking;
More informationParametric Texture Model based on Joint Statistics
Parametric Texture Model based on Joint Statistics Gowtham Bellala, Kumar Sricharan, Jayanth Srinivasa Department of Electrical Engineering, University of Michigan, Ann Arbor 1. INTRODUCTION Texture images
More information