Vectors [and more on masks] Vector space theory applies directly to several image processing/ representation problems

Similar documents
Image Gradients and Gradient Filtering Computer Vision

Filtering and Edge Detection

Edge Detection. Image Processing - Computer Vision

Edge Detection. Introduction to Computer Vision. Useful Mathematics Funcs. The bad news

Machine vision. Summary # 4. The mask for Laplacian is given

CS 4495 Computer Vision. Frequency and Fourier Transforms. Aaron Bobick School of Interactive Computing. Frequency and Fourier Transform

Image Enhancement in the frequency domain. GZ Chapter 4

Chapter 7 Wavelets and Multiresolution Processing

Machine vision, spring 2018 Summary 4

Edge Detection. CS 650: Computer Vision

Local enhancement. Local Enhancement. Local histogram equalized. Histogram equalized. Local Contrast Enhancement. Fig 3.23: Another example

Advanced Edge Detection 1

Lecture 7: Edge Detection

Local Enhancement. Local enhancement

GBS765 Electron microscopy

Math Real Analysis II

Convolution Spatial Aliasing Frequency domain filtering fundamentals Applications Image smoothing Image sharpening

Frequency2: Sampling and Aliasing

DISCRETE FOURIER TRANSFORM

Empirical Mean and Variance!

Announcements. Filtering. Image Filtering. Linear Filters. Example: Smoothing by Averaging. Homework 2 is due Apr 26, 11:59 PM Reading:

Laplacian Filters. Sobel Filters. Laplacian Filters. Laplacian Filters. Laplacian Filters. Laplacian Filters

Lecture 6: Edge Detection. CAP 5415: Computer Vision Fall 2008

ITK Filters. Thresholding Edge Detection Gradients Second Order Derivatives Neighborhood Filters Smoothing Filters Distance Map Image Transforms

Fourier series: Any periodic signals can be viewed as weighted sum. different frequencies. view frequency as an

Discrete Fourier Transform

ECE Digital Image Processing and Introduction to Computer Vision. Outline

I Chen Lin, Assistant Professor Dept. of CS, National Chiao Tung University. Computer Vision: 4. Filtering

Lecture 04 Image Filtering

Suppose we have a one-dimensional discrete signal f(x). The convolution of f with a kernel g is defined as h(x) = f(x)*g(x) = w 1. f(i)g(x i), (3.

Edge Detection in Computer Vision Systems

Fourier transforms and convolution

Edge Detection. Computer Vision P. Schrater Spring 2003

Image enhancement. Why image enhancement? Why image enhancement? Why image enhancement? Example of artifacts caused by image encoding

Roadmap. Introduction to image analysis (computer vision) Theory of edge detection. Applications

COMP344 Digital Image Processing Fall 2007 Final Examination

&& && F( u)! "{ f (x)} = f (x)e # j 2$ u x. f (x)! " #1. F(u,v) = f (x, y) e. f (x, y) = 2D Fourier Transform. Fourier Transform - review.

There are two things that are particularly nice about the first basis

Multiresolution schemes

Multiresolution schemes

Computer Vision. Filtering in the Frequency Domain

Review of Linear System Theory

The Frequency Domain, without tears. Many slides borrowed from Steve Seitz

Introduction to Computer Vision. 2D Linear Systems

Quadrature Formula for Computed Tomography

Nonlinear reverse-correlation with synthesized naturalistic noise

The Frequency Domain : Computational Photography Alexei Efros, CMU, Fall Many slides borrowed from Steve Seitz

Syllabus for IMGS-616 Fourier Methods in Imaging (RIT #11857) Week 1: 8/26, 8/28 Week 2: 9/2, 9/4

Scattered Interpolation Survey

CSE 473/573 Computer Vision and Image Processing (CVIP)

Math 2331 Linear Algebra

Spatial Enhancement Region operations: k'(x,y) = F( k(x-m, y-n), k(x,y), k(x+m,y+n) ]

MODULE 8 Topics: Null space, range, column space, row space and rank of a matrix

Math 24 Spring 2012 Sample Homework Solutions Week 8

(Refer Slide Time: 1:09)

Intro to harmonic analysis on groups Risi Kondor

Math 1180, Notes, 14 1 C. v 1 v n v 2. C A ; w n. A and w = v i w i : v w = i=1

Lecture 8: Interest Point Detection. Saad J Bedros

Invariant Scattering Convolution Networks

Additional Pointers. Introduction to Computer Vision. Convolution. Area operations: Linear filtering

arxiv: v1 [cs.cv] 10 Feb 2016

Computer Vision Lecture 3

Differential Operators for Edge Detection

A Possible Model of Noise EnhancedVisual Perception in Human Vision

Edge Detection PSY 5018H: Math Models Hum Behavior, Prof. Paul Schrater, Spring 2005

Wavelet-based Salient Points with Scale Information for Classification

These slides follow closely the (English) course textbook Pattern Recognition and Machine Learning by Christopher Bishop

EXAMINATION QUESTION PAPER

Image as a signal. Luc Brun. January 25, 2018

Image Filtering, Edges and Image Representation

Νεςπο-Ασαυήρ Υπολογιστική Neuro-Fuzzy Computing

CONTINUOUS IMAGE MATHEMATICAL CHARACTERIZATION

Reconnaissance d objetsd et vision artificielle

Reproducing Kernel Hilbert Spaces Class 03, 15 February 2006 Andrea Caponnetto

Digital Image Processing COSC 6380/4393

Image Filtering. Slides, adapted from. Steve Seitz and Rick Szeliski, U.Washington

Why does a lower resolution image still make sense to us? What do we lose? Image:

The Laplacian PDF Distance: A Cost Function for Clustering in a Kernel Feature Space

Lecture Notes 5: Multiresolution Analysis

MIT 2.71/2.710 Optics 10/31/05 wk9-a-1. The spatial frequency domain

Limit Theorems. MATH 464/506, Real Analysis. J. Robert Buchanan. Summer Department of Mathematics. J. Robert Buchanan Limit Theorems

Introduction to Computer Vision

Lagrange Multipliers

Sparse linear models

Wavelet Transform And Principal Component Analysis Based Feature Extraction

Vlad Estivill-Castro (2016) Robots for People --- A project for intelligent integrated systems

Computer Vision & Digital Image Processing. Periodicity of the Fourier transform

CITS 4402 Computer Vision

CHAPTER 6 : LITERATURE REVIEW

P E R E N C O - C H R I S T M A S P A R T Y

Focus was on solving matrix inversion problems Now we look at other properties of matrices Useful when A represents a transformations.

Colorado School of Mines Image and Multidimensional Signal Processing

Luminance, disparity and range statistics in 3D natural scenes

Lecture 8: Interest Point Detection. Saad J Bedros

Artifical Neural Networks

Linear Operators and Fourier Transform

Multiscale Image Transforms

Wavelets and multiresolution representations. Time meets frequency

f(s) e -i n π s/l d s

MTH 309Y 37. Inner product spaces. = a 1 b 1 + a 2 b a n b n

Transcription:

Vectors [and more on masks] Vector space theory applies directly to several image processing/ representation problems 1

Image as a sum of basic images What if every person s portrait photo could be expressed as a sum of 20 special images? è We would only need 20 numbers to model any photo è sparse rep on our Smart card. 2

Efaces 100 x 100 images of faces are approximated by a subspace of only 4 100 x 100 images, the mean image plus a linear combination of the 3 most important eigenimages 3

The image as an expansion 4

Different bases, different properties revealed 5

Fundamental expansion 6

Basis gives structural parts 7

Vector space review, part 1 8

Vector space review, Part 2 2 9

A space of images in a vector space n M x N image of real intensity values has dimension D = M x N n Can concatenate all M rows to interpret an image as a D dimensional 1D vector n The vector space properties apply n The 2D structure of the image is NOT lost 10

Orthonormal basis vectors help 11

Represent S = [10, 15, 20] 12

Projection of vector U onto V 13

Normalized dot product Can now think about the angle between two signals, two faces, two text documents, 14

Every 2x2 neighborhood has some constant, some edge, and some line component Confirm that basis vectors are orthonormal 15

Roberts basis cont. If a neighborhood N has large dot product with a basis vector (image), then N is similar to that basis image. 16

Standard 3x3 image basis Structureless and relatively useless! 17

Frie-Chen basis Confirm that bases vectors are orthonormal 18

Structure from Frie-Chen expansion Expand N using Frie- Chen basis 19

Sinusoids provide a good basis 20

Sinusoids also model well in images 21

Operations using the Fourier basis 22

A few properties of 1D sinusoids They are orthogonal Are they orthonormal? 23

F(x,y) as a sum of sinusoids 24

Continuous 2D Fourier Transform To compute F(u,v) we do a dot product of our image f(x,y) with a specific sinusoid with frequencies u and v 26

Power spectrum from FT 27

Examples from images Done with HIPS in 1997 28

Descriptions of former spectra 29

Discrete Fourier Transform Do N x N dot products and determine where the energy is. High energy in parameters u and v means original image has similarity to those sinusoids. 30

Bandpass filtering 31

Convolution of two functions in the spatial domain is equivalent to pointwise multiplication in the frequency domain 32

LOG or DOG filter Laplacian of Gaussian Approx Difference of Gaussians 33

LOG filter properties 34

Mathematical model 35

1D model; rotate to create 2D model 36

1D Gaussian and 1 st derivative 37

2 nd derivative; then all 3 curves 38

Laplacian of Gaussian as 3x3 39

G(x,y): Mexican hat filter 40

Convolving LOG with region boundary creates a zero-crossing Mask h(x,y) Input f(x,y) Output f(x,y) * h(x,y) 41

42

LOG relates to animal vision 43

1D EX. Artificial Neural Network (ANN) for computing g(x) = f(x) * h(x) level 1 cells feed 3 level 2 cells level 2 cells integrate 3 level 1 input cells using weights [-1,2,-1] 44

Experience the Mach band effect 45

Simple model of a neuron 46

Output conditioning: threshold versus smoother output signal 47

Canny edge detector 51

Summary of LOG filter n Convenient filter shape n Boundaries detected as 0-crossings n Psychophysical evidence that animal visual systems might work this way (your testimony) n Physiological evidence that real NNs work as the ANNs 53