Best apertures from entropy minimization
|
|
- Randall Cross
- 5 years ago
- Views:
Transcription
1 Best apertures from entropy minimization Bart M. ter Haar Romeny, PhD Eindhoven University of Technology Excerpt from the book "Front-End Vision & Multi-Scale Image Analysis", Kluwer Academic Publishers, << MathVisionTools`; Frame > FalseD; PlotRange AllD; MathVisionTools 2008 Biomedical Image-Analysis, Technische Universiteit Eindhoven, the Netherlands This Mathematica add-on, version (July 23rd, 2010) is for academic use only. Please, contact Dr. Markus van Almsick or Prof. Bart ter Haar Romeny for bug reports and commercial use. MathVisionTools` is a package written at Eindhoven University of Technology. It contains the functions gd and rasterplot, used below. It can be downloaded from Derivation The derivation is based on: M. Nielsen, "From paradigm to algorithms in computer vision". PhD thesis, Datalogisk Institut Kopenhagen University, Denmark, Dept. of Computer Science, Universitetsparken 1, DK-2100 Kopenhagen 0, Denmark, ISSN We do a measurement Ø finite aperture. All locations treated similar Ø translation invariance. Measurement is linear Ø superposition principle. Ø The observation must be a convolution: hhxl=ÿ - LHaL ghx -al a. LHxL is the luminance, ghxl is our aperture, hhxl the result of our measurement. A. The aperture function ghxl should be a normalized filter: Ÿ - ghxl x = 1. B. The mean of the filter ghxl is Ÿ - x ghxl x = x 0 = 0. C. The width is the variance: Ÿ - x 2 ghxl x=s 2. The information or entropy of our filter is by definition: H = Ÿ - -ghxl ln ghxl x. We look for the ghxl for which the entropy is minimal, given the constraints: Ÿ - ghxl x = 1, Ÿ - x ghxl x= 0 andÿ - x 2 ghxl x= s 2. When we want to find a minimum under a set of given constraints, we apply a standard mathematical technique named variational methods. The energy E becomes: E = Ÿ - -ghxl ln ghxl x + l 1Ÿ - ghxl x+l 2Ÿ - x ghxl x +l 3 Ÿ - x 2 ghxl x and is minimum when E = 0. The l's are the Lagrange multipliers. g
2 2 << VariationalMethods`; var = VariationalDA g@xdlog@g@xdd+λ 1 g@xd+λ 2 xg@xd+λ 3 x 2 g@xd, g@xd, xe 1 Log@g@xDD+λ 1 +x λ 2 +x 2 λ 3 Solve@var == 0, g@xdd ::g@xd 1+λ 1 +x λ 2 +x2 λ 3 >> g@x_d = First@g@xD ê. Solve@var == 0, g@xddd 1+λ 1 +x λ 2 +x2 λ 3 An important first finding is that g is an exponential function. l 3 must be negative, otherwise the function explodes, which is physically unrealistic. The constraint equations are now: eqn1 = SimplifyB g@xd x == 1, λ3 < 0F λ 3 λ 1 λ λ 3 π eqn2 = SimplifyB xg@xd x 0, λ3 < 0F λ 1 λ2 2 4 λ 3 λ2 0 eqn3 = SimplifyB x 2 g@xd x == σ 2, λ 3 < 0F 1+λ 1 λ2 2 4 λ 3 π Iλ2 2 2 λ 3 M 4 H λ 3 L 5ê2 σ 2 Now we can solve for all three l's: solution = Solve@8eqn1, eqn2, eqn3<, 8λ 1, λ 2, λ 3 <D ::λ LogB 4 4 π 2 σ 4F, λ 2 0, λ σ 2>> g@x_, σ_d = Simplify@g@xD ê. Flatten@solutionD, σ > 0D x 2 2 σ 2 2 π σ which is the Gaussian function. A beautiful result. We have found the Gaussian as the unique solution to the set of constraints, which in principle are a formal statement of the uncommittment of the observation.
3 3 There are 11 known axiomatic derivations of the Gaussian kernel as the optimal aperture for uncommitted observations [Weickert 2002]. Relaxing the constraints gives other families, e.g.: - separability: Poisson scale-space [Duits, Felsberg] - tuned to spatial frequency: Gabor kernels Plot@g@x, 1D, 8x, 4, 4<D Derivatives of sampled, observed data All partial derivatives of the Gaussian kernel are solutions too of the diffusion equation. g := 1 x2 +y2 ExpB F; σ = 1; 2 π σ2 2 σ 2 Plot3D@g, 8x, 3.5, 3.5<, 8y, 3.5, 3.5<D
4 4 x gd, 8x, 3.5, 3.5<, 8y, 3.5, 3.5<D Plot3DAEvaluateA y ge, 8x, 3.5, 3.5<, 8y, 3.5, 3.5<E Plot3DAEvaluateA x I y gme, 8x, 3.5, 3.5<, 8y, 3.5, 3.5<E
5 5 Plot3DAEvaluateA x,x,x,x I y,y gme, 8x, 3.5, 3.5<, 8y, 3.5, 3.5<E The derivative of the observed data is 8L 0Hx, yl GHx, y; sl<, which can be rewritten as L 0 Hx, yl GHx, y; sl. x x Differentiation of discrete data is convolution with a Gaussian derivative æ Differentiation and observation are done in a single step: convolution with a Gaussian derivative kernel. æ Differentiation is now done by integration, namely by the convolution integral. æ The Gaussian kernel is the physical analogue of a mathematical point, the Gaussian derivative kernels are the physical analogons of the mathematical differential operators. Equivalence is reached for the limit when the scale of the Gaussian goes to zero: The function gd[im,nx,ny,σ] implements a convolution with a Gaussian derivative: im = Table@If@80 < x < 170&&80 < y < 170, 1, 0D, 8y, 1, 256<, 8x, 1, 256<D; σ = 1; im x = gd@im, 1, 0, σd; im y = gd@im, 0, 1, σd; grad = im x 2 +im y 2 ; p1 = RasterPlotê@9im, im x, im y, grad=; GraphicsGrid@8p1<, ImageSize 300D FrontEndVision Version 2.0 for Mathematica 6 The first order derivative of an image gives edges. Left: original test image LHx, yl, resolution Second: the derivative with respect to x: L at scale x s = 1 pixel. Note the positive and negative edges. Third: the derivative with respect to y: L at scale s = 1 pixel. Right: the gradient J L scale of s = 1 pixel which gives all edges. y x N2 + J L x N2 at a im = ColorSeparate@Import@"mr256.gif"DDP1, 1T;
6 6 p1 = RasterPlot@imD Lxx = gd@im, 2, 0, 1D; Lyy = gd@im, 0, 2, 1D; RasterPlot@Lxx + LyyD
7 7 The Laplacian of the image, as seen by a center-surround cell of scale s=1 pixel. p2 = ListContourPlot@Reverse@Lxx + LyyD, Contours 80<, ContourStyle Red, ContourShading None, Frame NoneD Show@p1, p2d æ The Gaussian kernel is the unique kernel that generates no spurious resolution. It is the blown-up physical point operator, the Gaussian derivatives are the blown-up physical multi-scale derivative operators.
8 8 up ddx.jpg", ImageSize 300D æ Because convolving is an integration, the Gaussian kernel has by definition a strong regularizing effect. æ Recently some interesting papers have shown the complete equivalence of Gaussian scale space regularization with a number of other methods for regularization such as splines, thin plate splines, graduated convexity etc. [Scherzer2000a, Nielsen1996b, Nielsen1997b]. æ The set of Gaussian derivative kernels (including the zeroth order derivative: the Gaussian kernel itself) forms a complete set of derivatives. This set is sometimes referred to as the N-jet. A scale-space is a stack of 2D images, where scale is the third dimension. And here are - a scale-space of the intensity (no derivatives, only blurred); - a scale-space of the Laplacian ( 2 L x L y 2 ). im = ColorSeparate@Import@"mr256.gif"DDP1, 1T; GraphicsGrid@8Table@RasterPlot@gD@im, 0, 0, E τ DD, 8τ, 0, 2.1,.3<D, Table@RasterPlot@gD@im, 2, 0, E τ D+gD@im, 0, 2, E τ DD, 8τ, 0, 2.1,.3<D<, ImageSize 450D
2. Foundations of scale-space
2. Foundations of scale-space 13 2. Foundations of scale-space "There are many paths to the top of the mountain, but the view is always the same" -Chinese proverb. 2.1 Constraints for an uncommitted front-end
More information3. The Gaussian kernel
3. The Gaussian kernel Of all things, man is the measure. Protagoras the Sophist (480-4 B.C.)
More information8. Differentiation and regularization
43 8. Regularization 8. Differentiation and regularization 8. Regularization Regularization is the technique to make data behave well when an operator is applied to them. Such data could e.g. be functions,
More informationScale-Time Kernels and Models
Scale-Time Kernels and Models Bart M. ter Haar Romeny 1, Luc M.J. Florack, and Mads Nielsen 3 1 Image Sciences Institute, Utrecht University, Utrecht, The Netherlands B.terHaarRomeny@isi.uu.nl, http://www.isi.uu.nl
More informationGeometry-driven diffusion
Geometry-driven diffusion In[5]:= FalseD; SetOptions@ListPlot, Joined > TrueD; MathVisionTools 2008 Biomedical
More informationDifferential Geometry for Image Processing
MSc TU/e Course: Differential Geometry for Image Processing Teachers: R. Duits MF 7.072 (responsible teacher, lecturer) E.J. Bekkers MF 7.074 (teacher, instructor) Course Description: This MSc course aims
More informationEdge Detection. CS 650: Computer Vision
CS 650: Computer Vision Edges and Gradients Edge: local indication of an object transition Edge detection: local operators that find edges (usually involves convolution) Local intensity transitions are
More information16. Deblurring Gaussian blur
6. Deblurring Gaussian blur 277 6. Deblurring Gaussian blur 6. Deblurring To discuss an application where really high order Gaussian derivatives are applied, we study the deblurring of Gaussian blur by
More informationScale and Statistics in Variable Conductance Diffusion
Scale and Statistics in Variable Conductance Diffusion Terry S. Yoo Department of Computer Science University of North Carolina Chapel Hill, NC 7599-3175 Introduction Non-linear diffusion or Variable Conductance
More informationFeature Vector Similarity Based on Local Structure
Feature Vector Similarity Based on Local Structure Evgeniya Balmachnova, Luc Florack, and Bart ter Haar Romeny Eindhoven University of Technology, P.O. Box 53, 5600 MB Eindhoven, The Netherlands {E.Balmachnova,L.M.J.Florack,B.M.terHaarRomeny}@tue.nl
More informationA Riemannian Framework for Denoising Diffusion Tensor Images
A Riemannian Framework for Denoising Diffusion Tensor Images Manasi Datar No Institute Given Abstract. Diffusion Tensor Imaging (DTI) is a relatively new imaging modality that has been extensively used
More informationAn Efficient Method for Tensor Voting Using Steerable Filters
An Efficient Method for Tensor Voting Using Steerable Filters Erik Franken 1, Markus van Almsick 1, Peter Rongen 2, Luc Florack 1,, and Bart ter Haar Romeny 1 1 Department of Biomedical Engineering, Technische
More informationECE Digital Image Processing and Introduction to Computer Vision
ECE592-064 Digital Image Processing and Introduction to Computer Vision Depart. of ECE, NC State University Instructor: Tianfu (Matt) Wu Spring 2017 Outline Recap, image degradation / restoration Template
More informationComputer Vision & Digital Image Processing
Computer Vision & Digital Image Processing Image Restoration and Reconstruction I Dr. D. J. Jackson Lecture 11-1 Image restoration Restoration is an objective process that attempts to recover an image
More informationTRACKING and DETECTION in COMPUTER VISION Filtering and edge detection
Technischen Universität München Winter Semester 0/0 TRACKING and DETECTION in COMPUTER VISION Filtering and edge detection Slobodan Ilić Overview Image formation Convolution Non-liner filtering: Median
More informationFor continuous models, the so-called linear superposition operator for an image transformation can be expressed as:
Computational Vision U. Minn. Psy 5036 Linear systems, Convolutions and Optical blur ü Linear Systems Linear system models are important to vision for modeling: optical, retinal sampling, and neural transformations
More informationDetectors part II Descriptors
EECS 442 Computer vision Detectors part II Descriptors Blob detectors Invariance Descriptors Some slides of this lectures are courtesy of prof F. Li, prof S. Lazebnik, and various other lecturers Goal:
More informationLinear Operators and Fourier Transform
Linear Operators and Fourier Transform DD2423 Image Analysis and Computer Vision Mårten Björkman Computational Vision and Active Perception School of Computer Science and Communication November 13, 2013
More informationConvolution and Linear Systems
CS 450: Introduction to Digital Signal and Image Processing Bryan Morse BYU Computer Science Introduction Analyzing Systems Goal: analyze a device that turns one signal into another. Notation: f (t) g(t)
More informationITK Filters. Thresholding Edge Detection Gradients Second Order Derivatives Neighborhood Filters Smoothing Filters Distance Map Image Transforms
ITK Filters Thresholding Edge Detection Gradients Second Order Derivatives Neighborhood Filters Smoothing Filters Distance Map Image Transforms ITCS 6010:Biomedical Imaging and Visualization 1 ITK Filters:
More informationImage enhancement. Why image enhancement? Why image enhancement? Why image enhancement? Example of artifacts caused by image encoding
13 Why image enhancement? Image enhancement Example of artifacts caused by image encoding Computer Vision, Lecture 14 Michael Felsberg Computer Vision Laboratory Department of Electrical Engineering 12
More informationSURF Features. Jacky Baltes Dept. of Computer Science University of Manitoba WWW:
SURF Features Jacky Baltes Dept. of Computer Science University of Manitoba Email: jacky@cs.umanitoba.ca WWW: http://www.cs.umanitoba.ca/~jacky Salient Spatial Features Trying to find interest points Points
More informationFraunhofer Institute for Computer Graphics Research Interactive Graphics Systems Group, TU Darmstadt Fraunhoferstrasse 5, Darmstadt, Germany
Scale Space and PDE methods in image analysis and processing Arjan Kuijper Fraunhofer Institute for Computer Graphics Research Interactive Graphics Systems Group, TU Darmstadt Fraunhoferstrasse 5, 64283
More informationIntroduction to Nonlinear Image Processing
Introduction to Nonlinear Image Processing 1 IPAM Summer School on Computer Vision July 22, 2013 Iasonas Kokkinos Center for Visual Computing Ecole Centrale Paris / INRIA Saclay Mean and median 2 Observations
More informationInvertible Orientation Bundles on 2D Scalar Images
Invertible Orientation Bundles on 2D Scalar Images Stiliyan N. Kalitzin, Bart M ter Haar Romeny, Max A Viergever Image Science Institute, Utrecht University Abstract. A general approach for multiscale
More informationLinear Diffusion. E9 242 STIP- R. Venkatesh Babu IISc
Linear Diffusion Derivation of Heat equation Consider a 2D hot plate with Initial temperature profile I 0 (x, y) Uniform (isotropic) conduction coefficient c Unit thickness (along z) Problem: What is temperature
More informationConsistent Positive Directional Splitting of Anisotropic Diffusion
In: Boštjan Likar (ed.): Proc. of Computer Vision Winter Workshop, Bled, Slovenia, Feb. 7-9, 2001, pp. 37-48. Consistent Positive Directional Splitting of Anisotropic Diffusion Pavel Mrázek and Mirko Navara
More informationf(s) e -i n π s/l d s
Pointwise convergence of complex Fourier series Let f(x) be a periodic function with period l defined on the interval [,l]. The complex Fourier coefficients of f( x) are This leads to a Fourier series
More informationFiltering and Edge Detection
Filtering and Edge Detection Local Neighborhoods Hard to tell anything from a single pixel Example: you see a reddish pixel. Is this the object s color? Illumination? Noise? The next step in order of complexity
More informationGeneralized Laplacian as Focus Measure
Generalized Laplacian as Focus Measure Muhammad Riaz 1, Seungjin Park, Muhammad Bilal Ahmad 1, Waqas Rasheed 1, and Jongan Park 1 1 School of Information & Communications Engineering, Chosun University,
More informationMultiscale Image Transforms
Multiscale Image Transforms Goal: Develop filter-based representations to decompose images into component parts, to extract features/structures of interest, and to attenuate noise. Motivation: extract
More informationReading. 3. Image processing. Pixel movement. Image processing Y R I G Q
Reading Jain, Kasturi, Schunck, Machine Vision. McGraw-Hill, 1995. Sections 4.-4.4, 4.5(intro), 4.5.5, 4.5.6, 5.1-5.4. 3. Image processing 1 Image processing An image processing operation typically defines
More informationEdge Detection PSY 5018H: Math Models Hum Behavior, Prof. Paul Schrater, Spring 2005
Edge Detection PSY 5018H: Math Models Hum Behavior, Prof. Paul Schrater, Spring 2005 Gradients and edges Points of sharp change in an image are interesting: change in reflectance change in object change
More informationAbstract In this thesis I will examine the relations between turbulence as found in physical systems and turbulence in optic ow elds. I will record so
Turbulence in Optical Flow Fields Kim Steenstrup Pedersen Technical Report DIKU-TR-/3 Department of Computer Science University of Copenhagen Universitetsparken 1 DK-21 Kbenhavn DENMARK 3th January 21
More informationCS 179: LECTURE 16 MODEL COMPLEXITY, REGULARIZATION, AND CONVOLUTIONAL NETS
CS 179: LECTURE 16 MODEL COMPLEXITY, REGULARIZATION, AND CONVOLUTIONAL NETS LAST TIME Intro to cudnn Deep neural nets using cublas and cudnn TODAY Building a better model for image classification Overfitting
More informationIntroduction to Computer Vision. 2D Linear Systems
Introduction to Computer Vision D Linear Systems Review: Linear Systems We define a system as a unit that converts an input function into an output function Independent variable System operator or Transfer
More informationCalculus Overview. f(x) f (x) is slope. I. Single Variable. A. First Order Derivative : Concept : measures slope of curve at a point.
Calculus Overview I. Single Variable A. First Order Derivative : Concept : measures slope of curve at a point. Notation : Let y = f (x). First derivative denoted f ʹ (x), df dx, dy dx, f, etc. Example
More information6.4 Directional derivatives
103 6.5 First order gauge coordinates for rotation in this way means something different than a rotation of the image. There would be no way otherwise to recognize rotated images from non-rotated ones!
More informationTotal Variation Image Edge Detection
Total Variation Image Edge Detection PETER NDAJAH Graduate School of Science and Technology, Niigata University, 8050, Ikarashi 2-no-cho, Nishi-ku, Niigata, 950-28, JAPAN ndajah@telecom0.eng.niigata-u.ac.jp
More informationDeep Learning. Convolutional Neural Network (CNNs) Ali Ghodsi. October 30, Slides are partially based on Book in preparation, Deep Learning
Convolutional Neural Network (CNNs) University of Waterloo October 30, 2015 Slides are partially based on Book in preparation, by Bengio, Goodfellow, and Aaron Courville, 2015 Convolutional Networks Convolutional
More informationReview of Fourier transforms
Review of Fourier transforms Note - a lot of what we discuss here applies to Laplace transforms as well by a 90 degree rotation in the complex plane. I will focus on Fourier here as it this directly relevant
More informationEdge Detection. Computer Vision P. Schrater Spring 2003
Edge Detection Computer Vision P. Schrater Spring 2003 Simplest Model: (Canny) Edge(x) = a U(x) + n(x) U(x)? x=0 Convolve image with U and find points with high magnitude. Choose value by comparing with
More informationECE Digital Image Processing and Introduction to Computer Vision. Outline
2/9/7 ECE592-064 Digital Image Processing and Introduction to Computer Vision Depart. of ECE, NC State University Instructor: Tianfu (Matt) Wu Spring 207. Recap Outline 2. Sharpening Filtering Illustration
More informationEdge Detection in Computer Vision Systems
1 CS332 Visual Processing in Computer and Biological Vision Systems Edge Detection in Computer Vision Systems This handout summarizes much of the material on the detection and description of intensity
More informationFundamental Solutions and Green s functions. Simulation Methods in Acoustics
Fundamental Solutions and Green s functions Simulation Methods in Acoustics Definitions Fundamental solution The solution F (x, x 0 ) of the linear PDE L {F (x, x 0 )} = δ(x x 0 ) x R d Is called the fundamental
More informationToday s lecture. Local neighbourhood processing. The convolution. Removing uncorrelated noise from an image The Fourier transform
Cris Luengo TD396 fall 4 cris@cbuuse Today s lecture Local neighbourhood processing smoothing an image sharpening an image The convolution What is it? What is it useful for? How can I compute it? Removing
More informationLecture 6: Edge Detection. CAP 5415: Computer Vision Fall 2008
Lecture 6: Edge Detection CAP 5415: Computer Vision Fall 2008 Announcements PS 2 is available Please read it by Thursday During Thursday lecture, I will be going over it in some detail Monday - Computer
More informationSpatial Enhancement Region operations: k'(x,y) = F( k(x-m, y-n), k(x,y), k(x+m,y+n) ]
CEE 615: Digital Image Processing Spatial Enhancements 1 Spatial Enhancement Region operations: k'(x,y) = F( k(x-m, y-n), k(x,y), k(x+m,y+n) ] Template (Windowing) Operations Template (window, box, kernel)
More informationMATHEMATICAL MODEL OF IMAGE DEGRADATION. = s
MATHEMATICAL MODEL OF IMAGE DEGRADATION H s u v G s u v F s u v ^ F u v G u v H s u v Gaussian Kernel Source: C. Rasmussen Gaussian filters pixel 5 pixels 0 pixels 30 pixels Gaussian filter Removes high-frequency
More informationComputer Vision. Filtering in the Frequency Domain
Computer Vision Filtering in the Frequency Domain Filippo Bergamasco (filippo.bergamasco@unive.it) http://www.dais.unive.it/~bergamasco DAIS, Ca Foscari University of Venice Academic year 2016/2017 Introduction
More informationCSE 473/573 Computer Vision and Image Processing (CVIP)
CSE 473/573 Computer Vision and Image Processing (CVIP) Ifeoma Nwogu inwogu@buffalo.edu Lecture 11 Local Features 1 Schedule Last class We started local features Today More on local features Readings for
More informationLecture 7: Edge Detection
#1 Lecture 7: Edge Detection Saad J Bedros sbedros@umn.edu Review From Last Lecture Definition of an Edge First Order Derivative Approximation as Edge Detector #2 This Lecture Examples of Edge Detection
More informationMotion Estimation (I) Ce Liu Microsoft Research New England
Motion Estimation (I) Ce Liu celiu@microsoft.com Microsoft Research New England We live in a moving world Perceiving, understanding and predicting motion is an important part of our daily lives Motion
More informationCrossing Preserving Multi Scale Vesselness
Crossing Preserving Multi Scale Vesselness Julius Hannink, Remco Duits and Erik Bekkers Eindhoven University of Technology, Department of Biomedical Engineering and Department of Mathematics and Computer
More informationErkut Erdem. Hacettepe University February 24 th, Linear Diffusion 1. 2 Appendix - The Calculus of Variations 5.
LINEAR DIFFUSION Erkut Erdem Hacettepe University February 24 th, 2012 CONTENTS 1 Linear Diffusion 1 2 Appendix - The Calculus of Variations 5 References 6 1 LINEAR DIFFUSION The linear diffusion (heat)
More informationEECS564 Estimation, Filtering, and Detection Exam 2 Week of April 20, 2015
EECS564 Estimation, Filtering, and Detection Exam Week of April 0, 015 This is an open book takehome exam. You have 48 hours to complete the exam. All work on the exam should be your own. problems have
More informationRiesz-Transforms vs. Derivatives: On the Relationship Between the Boundary Tensor and the Energy Tensor
Riesz-Transforms vs. Derivatives: On the Relationship Between the Boundary Tensor and the Energy Tensor Ullrich Köthe 1 and Michael Felsberg 1 Cognitive Systems Group, University of Hamburg koethe@informatik.uni-hamburg.de
More informationFront End Vision: A Multiscale Geometry Engine
Front End Vision: A Multiscale Geometry Engine Scale-Space Theory in Computer Vision versus Front-End Biological Vision Bart ter Haar Romeny, PhD & Luc M.J.Florack 2, PhD Utrecht University : Image Sciences
More informationCS 556: Computer Vision. Lecture 21
CS 556: Computer Vision Lecture 21 Prof. Sinisa Todorovic sinisa@eecs.oregonstate.edu 1 Meanshift 2 Meanshift Clustering Assumption: There is an underlying pdf governing data properties in R Clustering
More informationCHAPTER 4 PRINCIPAL COMPONENT ANALYSIS-BASED FUSION
59 CHAPTER 4 PRINCIPAL COMPONENT ANALYSIS-BASED FUSION 4. INTRODUCTION Weighted average-based fusion algorithms are one of the widely used fusion methods for multi-sensor data integration. These methods
More informationSpatially adaptive alpha-rooting in BM3D sharpening
Spatially adaptive alpha-rooting in BM3D sharpening Markku Mäkitalo and Alessandro Foi Department of Signal Processing, Tampere University of Technology, P.O. Box FIN-553, 33101, Tampere, Finland e-mail:
More informationGalilean-diagonalized spatio-temporal interest operators
Galilean-diagonalized spatio-temporal interest operators Tony Lindeberg, Amir Akbarzadeh and Ivan Laptev Computational Vision and Active Perception Laboratory (CVAP) Department of Numerical Analysis and
More informationDiscrete Fourier Transform
Discrete Fourier Transform DD2423 Image Analysis and Computer Vision Mårten Björkman Computational Vision and Active Perception School of Computer Science and Communication November 13, 2013 Mårten Björkman
More informationSIFT: SCALE INVARIANT FEATURE TRANSFORM BY DAVID LOWE
SIFT: SCALE INVARIANT FEATURE TRANSFORM BY DAVID LOWE Overview Motivation of Work Overview of Algorithm Scale Space and Difference of Gaussian Keypoint Localization Orientation Assignment Descriptor Building
More informationA REDUCED COMPLEXITY TWO-DIMENSIONAL BCJR DETECTOR FOR HOLOGRAPHIC DATA STORAGE SYSTEMS WITH PIXEL MISALIGNMENT
A REDUCED COMPLEXITY TWO-DIMENSIONAL BCJR DETECTOR FOR HOLOGRAPHIC DATA STORAGE SYSTEMS WITH PIXEL MISALIGNMENT 1 S. Iman Mossavat, 2 J.W.M.Bergmans 1 iman@nus.edu.sg 1 National University of Singapore,
More informationEnergy Tensors: Quadratic, Phase Invariant Image Operators
Energy Tensors: Quadratic, Phase Invariant Image Operators Michael Felsberg and Erik Jonsson Linköping University, Computer Vision Laboratory, SE-58183 Linköping, Sweden Abstract. In this paper we briefly
More informationIntroduction to the Discrete Fourier Transform
Introduction to the Discrete ourier Transform Lucas J. van Vliet www.ph.tn.tudelft.nl/~lucas TNW: aculty of Applied Sciences IST: Imaging Science and technology PH: Linear Shift Invariant System A discrete
More informationPDEs in Image Processing, Tutorials
PDEs in Image Processing, Tutorials Markus Grasmair Vienna, Winter Term 2010 2011 Direct Methods Let X be a topological space and R: X R {+ } some functional. following definitions: The mapping R is lower
More informationWhy does a lower resolution image still make sense to us? What do we lose? Image:
2D FREQUENCY DOMAIN The slides are from several sources through James Hays (Brown); Srinivasa Narasimhan (CMU); Silvio Savarese (U. of Michigan); Bill Freeman and Antonio Torralba (MIT), including their
More informationarxiv: v1 [cs.cv] 10 Feb 2016
GABOR WAVELETS IN IMAGE PROCESSING David Bařina Doctoral Degree Programme (2), FIT BUT E-mail: xbarin2@stud.fit.vutbr.cz Supervised by: Pavel Zemčík E-mail: zemcik@fit.vutbr.cz arxiv:162.338v1 [cs.cv]
More informationSCALE-SPACE - Theory and Applications
SCALE-SPACE - Theory and Applications Scale is embedded in the task: do you want only coins or TREASURE? SCALE-SPACE Theory and Applications - Scale-space theory is a framework of multiscale image/signal
More informationTwo-Dimensional Clifford Windowed Fourier Transform
Two-Dimensional Clifford Windowed Fourier Transform Mawardi Bahri, Eckhard M. S. Hitzer and Sriwulan Adji Abstract Recently several generalizations to higher dimension of the classical Fourier transform
More informationLinear Diffusion and Image Processing. Outline
Outline Linear Diffusion and Image Processing Fourier Transform Convolution Image Restoration: Linear Filtering Diffusion Processes for Noise Filtering linear scale space theory Gauss-Laplace pyramid for
More informationReview Smoothing Spatial Filters Sharpening Spatial Filters. Spatial Filtering. Dr. Praveen Sankaran. Department of ECE NIT Calicut.
Spatial Filtering Dr. Praveen Sankaran Department of ECE NIT Calicut January 7, 203 Outline 2 Linear Nonlinear 3 Spatial Domain Refers to the image plane itself. Direct manipulation of image pixels. Figure:
More informationOptic Flow Computation from Cardiac MR Tagging Using a Multiscale Differential Method
Optic Flow Computation from Cardiac MR Tagging Using a Multiscale Differential Method A Comparative Study with Velocity-Encoded MRI A. Suinesiaputra 1,2, L.M.J. Florack 2, J.J.M. Westenberg 1, B.M. ter
More informationNON-LINEAR DIFFUSION FILTERING
NON-LINEAR DIFFUSION FILTERING Chalmers University of Technology Page 1 Summary Introduction Linear vs Nonlinear Diffusion Non-Linear Diffusion Theory Applications Implementation References Page 2 Introduction
More informationInvertible Apertured Orientation Filters in Image Analysis
International Journal of Computer Vision 31(2/3), 145 158 (1999) c 1999 Kluwer Academic Publishers. Manufactured in The Netherlands. Invertible Apertured Orientation Filters in Image Analysis STILIYAN
More informationPDE Based Image Diffusion and AOS
PDE Based Image Diffusion and AOS by Jarno Ralli April 27, 204 Abstract This technical paper shows how AOS Additive Operator Splitting scheme, together with TDMA Tri-Diagonal Matrix Algorithm, can be used
More informationAdditional Pointers. Introduction to Computer Vision. Convolution. Area operations: Linear filtering
Additional Pointers Introduction to Computer Vision CS / ECE 181B andout #4 : Available this afternoon Midterm: May 6, 2004 W #2 due tomorrow Ack: Prof. Matthew Turk for the lecture slides. See my ECE
More informationErrata for Robot Vision
Errata for Robot Vision This is a list of known nontrivial bugs in Robot Vision 1986) by B.K.P. Horn, MIT Press, Cambridge, MA ISBN 0-6-08159-8 and McGraw-Hill, New York, NY ISBN 0-07-030349-5. Thanks
More informationThe GET Operator. LiTH-ISY-R-2633
The GET Operator Michael Felsberg Linköping University, Computer Vision Laboratory, SE-58183 Linköping, Sweden, mfe@isy.liu.se, http://www.isy.liu.se/cvl/ October 4, 004 LiTH-ISY-R-633 Abstract In this
More informationProblem 1 (Exercise 2.2, Monograph)
MS&E 314/CME 336 Assignment 2 Conic Linear Programming January 3, 215 Prof. Yinyu Ye 6 Pages ASSIGNMENT 2 SOLUTIONS Problem 1 (Exercise 2.2, Monograph) We prove the part ii) of Theorem 2.1 (Farkas Lemma
More informationFront-End Vision and Multi-Scale Image Analysis
Front-End Vision and Multi-Scale Image Analysis Computational Imaging and Vision Managing Editor MAX A. VIERGEVER Utrecht University, Utrecht, The Netherlands Editorial Board GUNILLA BORGEFORS, Centre
More informationSparse linear models
Sparse linear models Optimization-Based Data Analysis http://www.cims.nyu.edu/~cfgranda/pages/obda_spring16 Carlos Fernandez-Granda 2/22/2016 Introduction Linear transforms Frequency representation Short-time
More informationA New Scheme for Anisotropic Diffusion: Performance Evaluation of Rotation Invariance, Dissipativity and Efficiency no author given no address given A
A New Scheme for Anisotropic Diffusion: Performance Evaluation of Rotation Invariance, Dissipativity and Efficiency Hanno Scharr 1 and Joachim Weickert 1 Interdisciplinary Center for Scientific Computing,
More informationVariations. ECE 6540, Lecture 10 Maximum Likelihood Estimation
Variations ECE 6540, Lecture 10 Last Time BLUE (Best Linear Unbiased Estimator) Formulation Advantages Disadvantages 2 The BLUE A simplification Assume the estimator is a linear system For a single parameter
More informationKAZE Features. 1 Introduction. Pablo Fernández Alcantarilla 1, Adrien Bartoli 1, and Andrew J. Davison 2
KAZE Features Pablo Fernández Alcantarilla 1, Adrien Bartoli 1, and Andrew J. Davison 2 1 ISIT-UMR 6284 CNRS, Université d Auvergne, Clermont Ferrand, France {pablo.alcantarilla,adrien.bartoli}@gmail.com
More informationPartial Differential Equations and Image Processing. Eshed Ohn-Bar
Partial Differential Equations and Image Processing Eshed Ohn-Bar OBJECTIVES In this presentation you will 1) Learn what partial differential equations are and where do they arise 2) Learn how to discretize
More informationNormative theory of visual receptive fields
Normative theory of visual receptive fields Tony Lindeberg Computational Brain Science Lab, Department of Computational Science and Technology, KTH Royal Institute of Technology, SE- 44 Stockholm, Sweden.
More informationLaplacian Eigenimages in Discrete Scale Space
Laplacian Eigenimages in Discrete Scale Space Martin Tschirsich and Arjan Kuijper,2 Technische Universität Darmstadt, Germany 2 Fraunhofer IGD, Darmstadt, Germany Abstract. Linear or Gaussian scale space
More informationMath 5588 Final Exam Solutions
Math 5588 Final Exam Solutions Prof. Jeff Calder May 9, 2017 1. Find the function u : [0, 1] R that minimizes I(u) = subject to u(0) = 0 and u(1) = 1. 1 0 e u(x) u (x) + u (x) 2 dx, Solution. Since the
More informationNotes on Regularization and Robust Estimation Psych 267/CS 348D/EE 365 Prof. David J. Heeger September 15, 1998
Notes on Regularization and Robust Estimation Psych 67/CS 348D/EE 365 Prof. David J. Heeger September 5, 998 Regularization. Regularization is a class of techniques that have been widely used to solve
More informationErrata and Updates for ASM Exam MAS-I (First Edition) Sorted by Date
Errata for ASM Exam MAS-I Study Manual (First Edition) Sorted by Date 1 Errata and Updates for ASM Exam MAS-I (First Edition) Sorted by Date Practice Exam 5 Question 6 is defective. See the correction
More informationErrata and Updates for ASM Exam MAS-I (First Edition) Sorted by Page
Errata for ASM Exam MAS-I Study Manual (First Edition) Sorted by Page 1 Errata and Updates for ASM Exam MAS-I (First Edition) Sorted by Page Practice Exam 5 Question 6 is defective. See the correction
More informationBlobs & Scale Invariance
Blobs & Scale Invariance Prof. Didier Stricker Doz. Gabriele Bleser Computer Vision: Object and People Tracking With slides from Bebis, S. Lazebnik & S. Seitz, D. Lowe, A. Efros 1 Apertizer: some videos
More informationAdvanced Machine Learning & Perception
Advanced Machine Learning & Perception Instructor: Tony Jebara Topic 1 Introduction, researchy course, latest papers Going beyond simple machine learning Perception, strange spaces, images, time, behavior
More informationOn the FPA infrared camera transfer function calculation
On the FPA infrared camera transfer function calculation (1) CERTES, Université Paris XII Val de Marne, Créteil, France (2) LTM, Université de Bourgogne, Le Creusot, France by S. Datcu 1, L. Ibos 1,Y.
More informationMore on Stochastics and the Phenomenon of Line-Edge Roughness
More on Stochastics and the Phenomenon of Line-Edge Roughness Chris A. Mack 34 th International Photopolymer Science and Technology Conference Chiba, Japan, June 28, Conclusions We need more than just
More informationLecture 04 Image Filtering
Institute of Informatics Institute of Neuroinformatics Lecture 04 Image Filtering Davide Scaramuzza 1 Lab Exercise 2 - Today afternoon Room ETH HG E 1.1 from 13:15 to 15:00 Work description: your first
More information