Achim J. Lilienthal Mobile Robotics and Olfaction Lab, AASS, Örebro University

Size: px
Start display at page:

Download "Achim J. Lilienthal Mobile Robotics and Olfaction Lab, AASS, Örebro University"

Transcription

1 Achim J. Lilienthal Mobile Robotics and Olfaction Lab, Room T29, Mo, -2 o'clock AASS, Örebro University (please drop me an in advance)

2 4.!!!!!!!!! Pre-Class Reading!!!!!!!!! Pre-Class Reading Schedule o Class "Course Introduction" (Mar 29, 22) o Class 2 "Introduction" (Apr 2, 22)» Gonzalez/Woods Chapter "Introduction"» Gonzalez/Woods Chapter 2 "Fundamentals"» Lecture Notes from last year o Class 3 "" (Apr 2, 22)» Gonzalez/Woods Chapter 3 "Intensity Transformations and "» Lecture Notes from last year o Class 4 "Bilateral Filtering/Fourier Domain" (Apr 7, 22)» "A Gentle Introduction to Bilateral Filtering and its Applications", Sylvain Paris, Pierre Kornprobst, Jack Tumblin, and Frédo Durand, SIGGRAPH 28» "Bilateral Filtering for Gray and Color Images", C. Tomasi, R. Manduchi, Proc. Int. Conf. Computer Vision» Gonzalez/Woods Chapter 4 "Filtering in the Frequency Domain"» Lecture Notes from last year DIP'2 A. J. Lilienthal (Apr 2, 22) 2

3 2. General Introduction Schedule Lectures. week 3: Thu, Mar 29, 22, :5-2: o'clock, T27 2. week 4: Tue, Apr 3, 22, 5:5-7: o'clock, T 3. week 5: Thu, Apr 2, 22, :5-2: o'clock, T27 4. week 6: Tue, Apr 7, 22, 5:5-7: o'clock, T week 6: Thu, Apr 9, 22, :5-2: o'clock, T27 6. week 7: Tue, Apr 24, 22, 5:5-7: o'clock, T 7. week 7: Thu, Apr 26, 22, :5-2: o'clock, T27 8. week 8: Thu, May 3, 22, :5-2: o'clock, T27 9. week 9: Tue, May 8, 22, 5:5-7: o'clock, T. week 9: Thu, May, 22, :5-2: o'clock, T27 Exercises. week 4: Thu, Apr 5, 22, 3:5-7: o'clock, T2 2. week 5: Thu, Apr 2, 22, 3:5-7: o'clock, T2 3. week 6: Thu, Apr 9, 22, 3:5-7: o'clock, T 3. week 8: Thu, May 3, 22, 3:5-7: o'clock, T2 DIP'2 A. J. Lilienthal (Apr 2, 22) 3

4 Contents. Image Enhancement in the Spatial Domain 2. Grey Level Transformations 3. Histogram Processing 4. Operations Involving Multiple Images Applications People Tracking DIP'2 A. J. Lilienthal (Apr 2, 22) 4

5 Image Enhancement in the Spatial Domain DIP'2 A. J. Lilienthal (Apr 2, 22) 5

6 . Image Enhancement in the Spatial Domain Image Enhancement o image processing o the result is supposed to be "more suitable"» "more suitable" according to a certain application more suitable for visual interpretation DIP'2 A. J. Lilienthal (Apr 2, 22) 6

7 . Image Enhancement in the Spatial Domain We want to create an image which is "better" in some sense o helps visual interpretation (brightening, sharpening ) subjective o pre-processing for a subsequent image analysis algorithm performance metric (performance of a task) o make the image more "specific" application dependent T f(x,y) g(x,y) original image (or set of images) new image DIP'2 A. J. Lilienthal (Apr 2, 22) 7

8 . Image Enhancement in the Spatial Domain Spatial Domain versus Frequency Domain o spatial domain» direct manipulation of the pixels discussed in this lecture» two types of transformations in the spatial domain: pixel brightness transformations, point processing (depend only on the pixel value itself) spatial filters, local transformations or local processing (depend on a small neighbourhood around the pixel) o frequency domain: modifications of the Fourier transform» discussed in coming lectures DIP'2 A. J. Lilienthal (Apr 2, 22) 9

9 . Image Enhancement in the Spatial Domain Transformations in the Spatial Domain g ( x, y) = T[ f ( x, y)] o standard approach: T is applied to a sub-image centred at (x,y) o sub-image is called mask (kernel, filter, template, window) o mask processing or filtering o T can operate on a set of images DIP'2 A. J. Lilienthal (Apr 2, 22)

10 . Image Enhancement in the Spatial Domain Transformations in the Spatial Domain g ( x, y) = T[ f ( x, y)] o fill new array with weighted sum of pixel values from the locations surrounding the corresponding location in the image using the same set of weights each time DIP'2 A. J. Lilienthal (Apr 2, 22)

11 2 Gray Level Transformations DIP'2 A. J. Lilienthal (Apr 2, 22) 2

12 2. Grey Level Transformations Grey Level Transformations o simplest case: each pixel in the output image depends only on the corresponding pixel in the input image o x neighbourhood (point processing) o example: contrast stretching s = T (r) s = T (r) DIP'2 A. J. Lilienthal (Apr 2, 22) 3

13 2. Grey Level Transformations Grey Level Transformations contrast stretching thresholding DIP'2 A. J. Lilienthal (Apr 2, 22) 4

14 2. Grey Level Transformations Grey Level Transformations f = imread('bubbles.tif'); fp = imadjust(f, [..9], [..],.5); imshow(fp); o imadjust» parameters always specified in [,]» values below. clipped to.» values above.9 clipped to.» image intensity reversed since. <. DIP'2 A. J. Lilienthal (Apr 2, 22) 5

15 2. Grey Level Transformations Grey Level Transformations f = imread('bubbles.tif'); fp = imadjust(f, [..9], [..],.5); imshow(fp); o imadjust» parameters always specified in [,]» values below. clipped to.» values above.9 clipped to.» image intensity reversed since. <.» gamma function parameter < g = f γ DIP'2 A. J. Lilienthal (Apr 2, 22) 6

16 2. Grey Level Transformations Grey Level Transformations f = imread('bubbles.tif'); fp = imadjust(f, [..9], [..],.5); imshow(fp); fp = imadjust(f, [.55.9], [..], 3); DIP'2 A. J. Lilienthal (Apr 2, 22) 7

17 2. Grey Level Transformations Grey Level Transformations f = imread('bubbles.tif'); fp = imadjust(f, [..9], [..],.5); imshow(fp); fp = imadjust(f, [.55.9], [..], 3); o imadjust» gamma function parameter > g = f γ DIP'2 A. J. Lilienthal (Apr 2, 22) 8

18 2. Grey Level Transformations Contrast Stretching o piecewise linear function o power law transformation (gamma transformation) γ s = cr DIP'2 A. J. Lilienthal (Apr 2, 22) 9

19 2. Grey Level Transformations Common Grey Level Transformations (Single Image) o linear» identity» inverse (negative) o power law» n. power» n. root o logarithmic DIP'2 A. J. Lilienthal (Apr 2, 22) 2

20 2. Grey Level Transformations Common Grey Level Transformations (Single Image) o inverse transform DIP'2 A. J. Lilienthal (Apr 2, 22) 2

21 2. Grey Level Transformations Common Grey Level Transformations (Single Image) o linear» identity» inverse o piecewise linear o power law (gamma)» n. power» n. root o logarithmic... with more than one input image o sum, mean o transformation based on statistical operations (variance, t-test ) DIP'2 A. J. Lilienthal (Apr 2, 22) 23

22 3 Histogram Processing DIP'2 A. J. Lilienthal (Apr 2, 22) 24

23 3. Histogram Processing Grey Scale Histogram o shows the number of pixels per grey level f = imread('bubbles.tif'); imhist(f); % displays the histogram % histogram display default DIP'2 A. J. Lilienthal (Apr 2, 22) 26

24 3. Histogram Processing Grey Scale Histogram o shows the number of pixels per grey level f = imread('bubbles.tif'); h = imhist(f); % default number of bins = 256 imhist(f,8); % number of bins = 8 DIP'2 A. J. Lilienthal (Apr 2, 22) 27

25 3. Histogram Processing Grey Scale Histogram o shows the number of pixels per grey level f = imread('bubbles.tif'); h = imhist(f); % default number of bins = 256 h = imhist(f,6); % number of bins = 6 hn = h/numel(f); % normalized histogram % numel num. of elements (pixels) bar(hn) % normalized histogram DIP'2 A. J. Lilienthal (Apr 2, 22) 28

26 3. Histogram Processing Grey Scale Histogram o shows the number of pixels per grey level f = imread('bubbles.tif'); h = imhist(f); % default number of bins = 256 h = imhist(f,6); % number of bins = 6 hn = h/numel(f); % normalized histogram % numel num. of elements (pixels) bar(hn) % normalized histogram DIP'2 A. J. Lilienthal (Apr 2, 22) 29

27 3. Histogram Processing Grey Scale Histogram o neutral transform DIP'2 A. J. Lilienthal (Apr 2, 22) 3

28 3. Histogram Processing Grey Scale Histogram o neutral transform o inverse transform DIP'2 A. J. Lilienthal (Apr 2, 22) 3

29 3. Histogram Processing Grey Scale Histogram o neutral transform o inverse transform o logarithmic transform DIP'2 A. J. Lilienthal (Apr 2, 22) 32

30 3. Histogram Processing Histogram Equalization o contrast / brightness adjustments sometimes need to be automatised o "optimal" contrast for an image? flat histogram DIP'2 A. J. Lilienthal (Apr 2, 22) 36

31 3. Histogram Processing Histogram Equalization o consider the continuous case: s, r [,] o probability density functions (PDFs) of s and r are related by gray levels as random variables! s = T (r) p s ( s) = p r ( r) dr ds = p r ( r) T ( r) o transformation function = cumulative density function (CDF) ds dr r T ( r) p r ( ω) dω r d = T ( r) = pr ( ω) dω = pr ( r) p s ( s) = dr DIP'2 A. J. Lilienthal (Apr 2, 22) 37

32 3. Histogram Processing Histogram Equalization o discrete case pr rk ) = nk n ( s = = = k T ( rk ) pr ( rj ) j= o does not generally produce a uniform PDF o tends to spread the histogram o enables automatic contrast stretching k k j= n j n DIP'2 A. J. Lilienthal (Apr 2, 22) 38

33 3. Histogram Processing Histogram Equalization CDF DIP'2 A. J. Lilienthal (Apr 2, 22) 39

34 3. Histogram Processing Histogram Equalization DIP'2 A. J. Lilienthal (Apr 2, 22) 4

35 3. Histogram Processing Histogram Equalization f = imread('bubbles.tif'); g = histeq(f, 256); imshow(g); f = imread('bubbles.tif'); g = histeq(f, 4); % 4 output levels imshow(g); DIP'2 A. J. Lilienthal (Apr 2, 22) 4

36 5 DIP'2 A. J. Lilienthal (Apr 2, 22) 69

37 Neighbourhood Relations Between Pixels o a pixel has 4 or 8 neighbours in 2D depending on the neighbour definition:» 4-neighborhood each neighbor must share an edge with the pixel» 8- neighborhood each neighbor must share an edge or a corner with the pixel DIP'2 A. J. Lilienthal (Apr 2, 22) 7

38 Basics of o the pixel value in the output image is calculated from a local neighbourhood in the input image o the local neighbourhood is described by a mask with a typical size of 3x3, 5x5, 7x7, pixels o filtering is performed by moving the mask over the image o the centre pixel in the output image is given a value that depends on the input image and the weights of the mask DIP'2 A. J. Lilienthal (Apr 2, 22) 7

39 Basics of o filter subimage defines coefficients w(s,t) o used to update pixel at (x,y) DIP'2 A. J. Lilienthal (Apr 2, 22) 72

40 Linear o filter subimage defines coefficients w(s,t) o response of the filter at point (x,y) is given by a sum of products g ( x, y) = w( s, t) f ( x + s, y + t) s= at= b o also called convolution (convolution kernel) a b (-,-) (,-) (,-) (-,) (,) (,) (-,) (,) (,) DIP'2 A. J. Lilienthal (Apr 2, 22) 73

41 Linear Implementation o generic code: for P(x,y) in image for F(u,v) in filter Q(x,y) += F(u,v) P(x-u,y-v) end end How to Deal With the Border? o limit excursion of the centre of the mask smaller image o set outside pixel value zero border effects o mirroring border pixel values border effects o modify filter size along the border slower DIP'2 A. J. Lilienthal (Apr 2, 22) 75

42 Smoothing Spatial Filters (Averaging Filters) o for blurring» removal of small (irrelevant) details, bridging small gaps o for noise reduction» but: edges are also blurred DIP'2 A. J. Lilienthal (Apr 2, 22) 79

43 Smoothing Spatial Filters (Averaging Filters) o for blurring» removal of small (irrelevant) details, bridging small gaps o for noise reduction Smoothing Spatial Filters 3x3 Mean Filter / Box Filter o need for normalization to conserve the total energy of the image x /9 (sum of all greylevels) o can cause "ringing" o no good model of blurring in a defocused camera» turns a single "point" into a "box" DIP'2 A. J. Lilienthal (Apr 2, 22) 8

44 Smoothing Spatial Filters Mean Filter original Mean 5x5 Mean x DIP'2 A. J. Lilienthal (Apr 2, 22) 8

45 Linear in Matlab f = imread('bubbles.tif'); g = imfilter(f, w, filtering_mode, boundary_options, size_options); o filter matrix w o filtering modes» 'corr' or 'conv'» only important in the case of asymmetric filters» 'corr' (no mirroring) is the default DIP'2 A. J. Lilienthal (Apr 2, 22) 82

46 Linear in Matlab f = imread('bubbles.tif'); g = imfilter(f, w, filtering_mode, boundary_options, size_options); o boundary options» P padding with (default)» 'replicate' replicate values at the outer border» 'symmetric' mirror reflecting across the outer border» 'circular' repeating the image like a periodic function DIP'2 A. J. Lilienthal (Apr 2, 22) 83

47 Linear in Matlab f = imread('bubbles.tif'); g = imfilter(f, w, filtering_mode, boundary_options, size_options); o size options» 'same' same size as the input image (cropped padded image)» 'full' full size of the padded image» default is 'same' DIP'2 A. J. Lilienthal (Apr 2, 22) 84

48 Linear in Matlab f = imread('bubbles.tif'); g = imfilter(f, ones(8)/64, 'replicate'); imshow(g); 72 px DIP'2 A. J. Lilienthal (Apr 2, 22) 85

49 Linear in Matlab f = imread('bubbles.tif'); g = imfilter(f, fspecial('average',32), 'replicate'); imshow(g); DIP'2 A. J. Lilienthal (Apr 2, 22) 86

50 Smoothing Spatial Filters Mean Filter o square box filter generates defects» axis aligned streaks» blocky results output input example from "A Gentle Introduction to Bilateral Filtering and its Applications", Sylvain Paris, Pierre Kornprobst, Jack Tumblin, and Frédo Durand, SIGGRAPH 28 DIP'2 A. J. Lilienthal (Apr 2, 22) 87

51 box profile pixel weight pixel position unrelated pixels related pixels unrelated pixels from "A Gentle Introduction to Bilateral Filtering and its Applications", S. Paris, P. Kornprobst, J. Tumblin, and F. Durand, SIGGRAPH 28 DIP'2 A. J. Lilienthal (Apr 2, 22) 88

52 strategy to solve problems with box filters o use an isotropic (i.e. circular) window o use a window with a smooth falloff box kernel Gaussian kernel from "A Gentle Introduction to Bilateral Filtering and its Applications", S. Paris, P. Kornprobst, J. Tumblin, and F. Durand, SIGGRAPH 28 DIP'2 A. J. Lilienthal (Apr 2, 22) 89

53 Smoothing Spatial Filters Gaussian Filter o weighted average o 2D Gaussian kernel o higher weight in the centre to decrease blurring Why a Gaussian? o simple model of blurring in optical systems o smooth o also a Gaussian in the frequency domain /6 x DIP'2 A. J. Lilienthal (Apr 2, 22) 9

54 5 input from "A Gentle Introduction to Bilateral Filtering and its Applications", S. Paris, P. Kornprobst, J. Tumblin, and F. Durand, SIGGRAPH 28 DIP'2 A. J. Lilienthal (Apr 2, 22) 9

55 5 box average input from "A Gentle Introduction to Bilateral Filtering and its Applications", S. Paris, P. Kornprobst, J. Tumblin, and F. Durand, SIGGRAPH 28 DIP'2 A. J. Lilienthal (Apr 2, 22) 92

56 5 Gaussian box average input blur from "A Gentle Introduction to Bilateral Filtering and its Applications", S. Paris, P. Kornprobst, J. Tumblin, and F. Durand, SIGGRAPH 28 DIP'2 A. J. Lilienthal (Apr 2, 22) 93

57 Gaussian profile pixel weight pixel position unrelated pixels uncertain pixels related pixels uncertain pixels unrelated pixels from "A Gentle Introduction to Bilateral Filtering and its Applications", S. Paris, P. Kornprobst, J. Tumblin, and F. Durand, SIGGRAPH 28 DIP'2 A. J. Lilienthal (Apr 2, 22) 94

58 Gaussian profile o spatial parameter σ input small σ large σ limited smoothing strong smoothing from "A Gentle Introduction to Bilateral Filtering and its Applications", S. Paris, P. Kornprobst, J. Tumblin, and F. Durand, SIGGRAPH 28 DIP'2 A. J. Lilienthal (Apr 2, 22) 95

59 Gaussian profile o spatial parameter σ o how to set σ?» depends on the application» common strategy: proportional to image size e.g. 2% of the image diagonal property: independent of image resolution» depends on image content smooth "object areas" larger σ but don't smooth edges smaller σ from "A Gentle Introduction to Bilateral Filtering and its Applications", S. Paris, P. Kornprobst, J. Tumblin, and F. Durand, SIGGRAPH 28 DIP'2 A. J. Lilienthal (Apr 2, 22) 96

60 Smoothing Spatial Filters Median Filter o take the values of the input image corresponding to the desired sub-window (3x3, 5x5, ) o sort them o take the middle value (example: 3x3 the 5th largest) o forces pixels with distinct grey levels to be more like their neighbours o very good at reducing salt-and-pepper noise o less blurring than linear filters of the same size DIP'2 A. J. Lilienthal (Apr 2, 22) 97

61 Smoothing Spatial Filters Median Filter o take the median value over the sub-window X ray image of a circuit board Average 3x3 Median 3x3 DIP'2 A. J. Lilienthal (Apr 2, 22) 98

62 Median Filter in Matlab f = imread('bubbles.tif'); g = medfilt2(f, [32 32]); imshow(g); original image median 8x8 median 32x32 median 32x2 DIP'2 A. J. Lilienthal (Apr 2, 22) 99

63 Smoothing Spatial Filters Median Filter o take the values of the input image corresponding to the desired sub-window (3x3, 5x5, ) o sort them o take the middle value (example: 3x3 the 5th largest) o forces pixels with distinct grey levels to be more like their neighbours o very good at reducing salt-and-pepper noise o less blurring than linear filters of the same size o nonlinear filter (order statistics filter) o no equivalent in the frequency domain (order statistics filter) DIP'2 A. J. Lilienthal (Apr 2, 22)

64 Order Statistics Filters (Fractile Filters) o median o min, max» useful in mathematical morphology o percentile» generalization of median, min, max % percentile min (%) median (5%) max (%) DIP'2 A. J. Lilienthal (Apr 2, 22)

65 Sharpening Spatial Filters o highlight fine detail (also noise) o enhance edges o use image differentiation ( st order) f f( x+ ε, y) f( x, y) = lim x ε ε DIP'2 A. J. Lilienthal (Apr 2, 22) 3

66 Sharpening Spatial Filters o highlight fine detail (also noise) o enhance edges o use image differentiation ( st order) f f( x+ ε, y) f( x, y) = lim x ε ε f f x f i+, j i, j DIP'2 A. J. Lilienthal (Apr 2, 22) 4

67 Sharpening Spatial Filters o highlight fine detail (also noise) o enhance edges o use image differentiation Sharpening Spatial Filters D o approximation to st order derivation» equivalent to the D convolution mask f x f ( x + ) f ( x) - DIP'2 A. J. Lilienthal (Apr 2, 22) 5

68 Sharpening Spatial Filters o highlight fine detail (also noise) o enhance edges o use image differentiation Sharpening Spatial Filters D o approximation to st order derivation» equivalent to the D convolution mask f f( x + ) f( x ) x - DIP'2 A. J. Lilienthal (Apr 2, 22) 6

69 7 DIP'2 A. J. Lilienthal (Apr 2, 22) Gradient and Magnitude of the Gradient Sharpening Spatial Filters Based on the Gradient o Prewitt o Sobel o y f x f y f x f f + + = 2 2 = y f x f f,

70 Sharpening Spatial Filters o Prewitt gradient edge detector» 2 masks approximate G x and G y in f f f + = G x + G y x y - - DIP'2 A. J. Lilienthal (Apr 2, 22) 9

71 Sharpening Spatial Filters o Sobel Operators» 2 masks approximate G x and G y in» detects horizontal and vertical edges f f f + = G x + G y x y DIP'2 A. J. Lilienthal (Apr 2, 22)

72 Sharpening Spatial Filters Sobel Operators o weight 2 is supposed to smooth by emphasizing the centre DIP'2 A. J. Lilienthal (Apr 2, 22)

73 Sharpening Spatial Filters Sobel Operators o detection of vertical dark-light edges DIP'2 A. J. Lilienthal (Apr 2, 22) 2

74 Sharpening Spatial Filters Sobel Operators o combination of all the directional responses DIP'2 A. J. Lilienthal (Apr 2, 22) 3

75 Sharpening Spatial Filters o comparison between Sobel and Prewitt operator Sobel (~ G x + G y ) Prewitt (~ G x + G y ) DIP'2 A. J. Lilienthal (Apr 2, 22) 4

76 Sharpening Spatial Filters o highlight fine detail (also noise) o enhance edges o uses image differentiation Sharpening Spatial Filters D o approximation to st order derivation o approximation to 2 nd order derivation 2 f = f ( x + ) + f ( x ) 2 f ( x) 2 x» equivalent to the D convolution mask -2 DIP'2 A. J. Lilienthal (Apr 2, 22) 6

77 7 DIP'2 A. J. Lilienthal (Apr 2, 22) Sharpening Spatial Filters Laplace Filter o Laplacian (second order derivative) y f x f + = ) ( 2 ) ( ) ( 2 2 x f x f x f x f + + =

78 Sharpening Spatial Filters Laplace Filter o Laplacian (second order derivative) f f f = + = f ( x + ) + f ( x ) x y x o filter masks to implement the Laplacian» add the "digital implementation" of the two terms in the Laplacian (9 rotation symmetry) f ( x) -4 DIP'2 A. J. Lilienthal (Apr 2, 22) 8

79 Sharpening Spatial Filters Laplace Filter o Laplacian (second order derivative) f f f = + = f ( x + ) + f ( x ) x y x o filter masks to implement the Laplacian» add the "digital implementation" of the two terms in the Laplacian (9 rotation symmetry)» add also diagonal terms (45 rotation symmetry) f ( x) -4-8 DIP'2 A. J. Lilienthal (Apr 2, 22) 9

80 Sharpening Spatial Filters Laplace Filter o Laplacian (second order derivative) f f f = + = f ( x + ) + f ( x ) x y x o filter masks to implement the Laplacian» add the "digital implementation" of the two terms in the Laplacian (9 rotation symmetry)» add also diagonal terms (45 rotation symmetry) o negative values re-scale f ( x) -4-8 DIP'2 A. J. Lilienthal (Apr 2, 22) 2

81 Sharpening Spatial Filters Laplace Filter o detection of edges independent of direction o isotropic with respect to 9 rotations DIP'2 A. J. Lilienthal (Apr 2, 22) 2

82 Sharpening Spatial Filters Laplace Filter o Laplace filter + original image sharpening DIP'2 A. J. Lilienthal (Apr 2, 22) 22

83 Sharpening Spatial Filters o Laplace filter + original image sharpening DIP'2 A. J. Lilienthal (Apr 2, 22) 23

84 Sharpening Spatial Filters 2 nd order vs. st order o Laplacian (second order derivative) f f f = + = f ( x + ) + f ( x ) 2 f x y x o thinner edges o not so strong response to a step o better response to fine details o double response to edges o rotation independent one mask for all edges ( x) DIP'2 A. J. Lilienthal (Apr 2, 22) 24

85 Sharpening Spatial Filters Unsharp Masking o analog equivalent used in publishing industry o basic idea: subtract blurred version of an image from original image to generate the edges DIP'2 A. J. Lilienthal (Apr 2, 22) 25

86 Sharpening Spatial Filters in Matlab f = imread('bubbles.tif'); g = imfilter(f, fspecial('laplacian',.5)); g2 = imfilter(f, fspecial('unsharp',.5)); %... DIP'2 A. J. Lilienthal (Apr 2, 22) 26

87 6 Preparation for Next Class DIP'2 A. J. Lilienthal (Apr 2, 22) 29

88 6.!!!!!!!!! Preparation for Next Class!!!!!!!!! Pre-Class Reading Schedule o Class "Course Introduction" (Mar 29, 22) o Class 2 "Introduction" (Apr 2, 22)» Gonzalez/Woods Chapter "Introduction"» Gonzalez/Woods Chapter 2 "Fundamentals"» Lecture Notes from last year o Class 3 "" (Apr 2, 22)» Gonzalez/Woods Chapter 3 "Intensity Transformations and "» Lecture Notes from last year o Class 4 "Bilateral Filtering/Fourier Domain" (Apr 7, 22)» "A Gentle Introduction to Bilateral Filtering and its Applications", Sylvain Paris, Pierre Kornprobst, Jack Tumblin, and Frédo Durand, SIGGRAPH 28» "Bilateral Filtering for Gray and Color Images", C. Tomasi, R. Manduchi, Proc. Int. Conf. Computer Vision» Gonzalez/Woods Chapter 4 "Filtering in the Frequency Domain"» Lecture Notes from last year DIP'2 A. J. Lilienthal (Apr 2, 22) 3

89 Achim J. Lilienthal Mobile Robotics and Olfaction Lab, Room T29, Mo, -2 o'clock AASS, Örebro University (please drop me an in advance) 3

Digital Image Processing

Digital Image Processing Digital Image Processing Part : Image Enhancement in the Spatial Domain AASS Learning Systems Lab, Dep. Teknik Room T9 (Fr, - o'clock) achim.lilienthal@oru.se Course Book Chapter 3-4- Contents. Image Enhancement

More information

Digital Image Processing

Digital Image Processing Digital Image Processing Part 2: Image Enhancement Digital Image Processing Course Introduction in the Spatial Domain Lecture AASS Learning Systems Lab, Teknik Room T26 achim.lilienthal@tech.oru.se Course

More information

CoE4TN4 Image Processing. Chapter 3: Intensity Transformation and Spatial Filtering

CoE4TN4 Image Processing. Chapter 3: Intensity Transformation and Spatial Filtering CoE4TN4 Image Processing Chapter 3: Intensity Transformation and Spatial Filtering Image Enhancement Enhancement techniques: to process an image so that the result is more suitable than the original image

More information

Practical Image and Video Processing Using MATLAB

Practical Image and Video Processing Using MATLAB Practical Image and Video Processing Using MATLAB Chapter 10 Neighborhood processing What will we learn? What is neighborhood processing and how does it differ from point processing? What is convolution

More information

Chapter 3 Image Enhancement in the Spatial Domain. Chapter 3 Image Enhancement in the Spatial Domain

Chapter 3 Image Enhancement in the Spatial Domain. Chapter 3 Image Enhancement in the Spatial Domain It makes all the difference whether one sees darkness through the light or brightness through the shadows. - David Lindsay 3.1 Background 76 3.2 Some Basic Gray Level Transformations 78 3.3 Histogram Processing

More information

IMAGE ENHANCEMENT IN SPATIAL DOMAIN

IMAGE ENHANCEMENT IN SPATIAL DOMAIN A First Course in Machine Vision IMAGE ENHANCEMENT IN SPATIAL DOMAIN By: Ehsan Khoramshahi Definitions The principal objective of enhancement is to process an image so that the result is more suitable

More information

TDI2131 Digital Image Processing

TDI2131 Digital Image Processing TDI2131 Digital Image Processing Image Enhancement in Spatial Domain Lecture 3 John See Faculty of Information Technology Multimedia University Some portions of content adapted from Zhu Liu, AT&T Labs.

More information

EE482: Digital Signal Processing Applications

EE482: Digital Signal Processing Applications Professor Brendan Morris, SEB 3216, brendan.morris@unlv.edu EE482: Digital Signal Processing Applications Spring 2014 TTh 14:30-15:45 CBC C222 Lecture 15 Image Processing 14/04/15 http://www.ee.unlv.edu/~b1morris/ee482/

More information

Image Filtering. Median Filtering

Image Filtering. Median Filtering Image Filtering Image filtering is used to: Remove noise Sharpen contrast Highlight contours Detect edges Other uses? Image filters can be classified as linear or nonlinear. Linear filters are also know

More information

Prof. Feng Liu. Winter /10/2019

Prof. Feng Liu. Winter /10/2019 Prof. Feng Liu Winter 29 http://www.cs.pdx.edu/~fliu/courses/cs4/ //29 Last Time Course overview Admin. Info Computer Vision Computer Vision at PSU Image representation Color 2 Today Filter 3 Today Filters

More information

Prof. Vidya Manian Dept. of Electrical and Comptuer Engineering

Prof. Vidya Manian Dept. of Electrical and Comptuer Engineering Image Processing Intensity Transformations Chapter 3 Prof. Vidya Manian Dept. of Electrical and Comptuer Engineering INEL 5327 ECE, UPRM Intensity Transformations 1 Overview Background Basic intensity

More information

Image Enhancement in spatial domain. Digital Image Processing GW Chapter 3 from Section (pag 110) Part 2: Filtering in spatial domain

Image Enhancement in spatial domain. Digital Image Processing GW Chapter 3 from Section (pag 110) Part 2: Filtering in spatial domain Image Enhancement in spatial domain Digital Image Processing GW Chapter 3 from Section 3.4.1 (pag 110) Part 2: Filtering in spatial domain Mask mode radiography Image subtraction in medical imaging 2 Range

More information

Filtering in the spatial domain (Spatial Filtering)

Filtering in the spatial domain (Spatial Filtering) Filtering in the spatial domain (Spatial Filtering) refers to image operators that change the gray value at any pixel (x,y) depending on the pixel values in a square neighborhood centered at (x,y) using

More information

Filip Malmberg 1TD396 fall 2018 Today s lecture

Filip Malmberg 1TD396 fall 2018 Today s lecture Today s lecture Local neighbourhood processing Convolution smoothing an image sharpening an image And more What is it? What is it useful for? How can I compute it? Removing uncorrelated noise from an image

More information

Digital Image Processing

Digital Image Processing Digital Image Processing Lecture # 5 Image Enhancement in Spatial Domain- I ALI JAVED Lecturer SOFTWARE ENGINEERING DEPARTMENT U.E.T TAXILA Email:: ali.javed@uettaxila.edu.pk Office Room #:: 7 Presentation

More information

Table of contents. Vision industrielle 2002/2003. Local and semi-local smoothing. Linear noise filtering: example. Convolution: introduction

Table of contents. Vision industrielle 2002/2003. Local and semi-local smoothing. Linear noise filtering: example. Convolution: introduction Table of contents Vision industrielle 2002/2003 Session - Image Processing Département Génie Productique INSA de Lyon Christian Wolf wolf@rfv.insa-lyon.fr Introduction Motivation, human vision, history,

More information

CS534 Introduction to Computer Vision. Linear Filters. Ahmed Elgammal Dept. of Computer Science Rutgers University

CS534 Introduction to Computer Vision. Linear Filters. Ahmed Elgammal Dept. of Computer Science Rutgers University CS534 Introduction to Computer Vision Linear Filters Ahmed Elgammal Dept. of Computer Science Rutgers University Outlines What are Filters Linear Filters Convolution operation Properties of Linear Filters

More information

Image Enhancement in the Spatial Domain

Image Enhancement in the Spatial Domain Image Enhancement in the Spatial Domain Algorithms for improving the visual appearance of images Gamma correction Contrast improvements Histogram equalization Noise reduction Image sharpening Optimality

More information

Image Enhancement. DD2423 Image Analysis and Computer Vision. Computational Vision and Active Perception School of Computer Science and Communication

Image Enhancement. DD2423 Image Analysis and Computer Vision. Computational Vision and Active Perception School of Computer Science and Communication Image Enhancement DD2423 Image Analysis and Computer Vision Mårten Björkman Computational Vision and Active Perception School of Computer Science and Communication November 15, 2013 Mårten Björkman (CVAP)

More information

Midterm Review. Image Processing CSE 166 Lecture 10

Midterm Review. Image Processing CSE 166 Lecture 10 Midterm Review Image Processing CSE 166 Lecture 10 Topics covered Image acquisition, geometric transformations, and image interpolation Intensity transformations Spatial filtering Fourier transform and

More information

Motion illusion, rotating snakes

Motion illusion, rotating snakes Motion illusion, rotating snakes Image Filtering 9/4/2 Computer Vision James Hays, Brown Graphic: unsharp mask Many slides by Derek Hoiem Next three classes: three views of filtering Image filters in spatial

More information

Image Enhancement using Histogram Equalization and Spatial Filtering

Image Enhancement using Histogram Equalization and Spatial Filtering Image Enhancement using Histogram Equalization and Spatial Filtering Fari Muhammad Abubakar 1 1 Department of Electronics Engineering Tianjin University of Technology and Education (TUTE) Tianjin, P.R.

More information

Non Linear Image Enhancement

Non Linear Image Enhancement Non Linear Image Enhancement SAIYAM TAKKAR Jaypee University of information technology, 2013 SIMANDEEP SINGH Jaypee University of information technology, 2013 Abstract An image enhancement algorithm based

More information

Motivation: Image denoising. How can we reduce noise in a photograph?

Motivation: Image denoising. How can we reduce noise in a photograph? Linear filtering Motivation: Image denoising How can we reduce noise in a photograph? Moving average Let s replace each pixel with a weighted average of its neighborhood The weights are called the filter

More information

CSE 564: Scientific Visualization

CSE 564: Scientific Visualization CSE 564: Scientific Visualization Lecture 5: Image Processing Klaus Mueller Stony Brook University Computer Science Department Klaus Mueller, Stony Brook 2003 Image Processing Definitions Purpose: - enhance

More information

Vision Review: Image Processing. Course web page:

Vision Review: Image Processing. Course web page: Vision Review: Image Processing Course web page: www.cis.udel.edu/~cer/arv September 7, Announcements Homework and paper presentation guidelines are up on web page Readings for next Tuesday: Chapters 6,.,

More information

Image Processing. Adam Finkelstein Princeton University COS 426, Spring 2019

Image Processing. Adam Finkelstein Princeton University COS 426, Spring 2019 Image Processing Adam Finkelstein Princeton University COS 426, Spring 2019 Image Processing Operations Luminance Brightness Contrast Gamma Histogram equalization Color Grayscale Saturation White balance

More information

Image acquisition. Midterm Review. Digitization, line of image. Digitization, whole image. Geometric transformations. Interpolation 10/26/2016

Image acquisition. Midterm Review. Digitization, line of image. Digitization, whole image. Geometric transformations. Interpolation 10/26/2016 Image acquisition Midterm Review Image Processing CSE 166 Lecture 10 2 Digitization, line of image Digitization, whole image 3 4 Geometric transformations Interpolation CSE 166 Transpose these matrices

More information

IMAGE PROCESSING: AREA OPERATIONS (FILTERING)

IMAGE PROCESSING: AREA OPERATIONS (FILTERING) IMAGE PROCESSING: AREA OPERATIONS (FILTERING) N. C. State University CSC557 Multimedia Computing and Networking Fall 2001 Lecture # 13 IMAGE PROCESSING: AREA OPERATIONS (FILTERING) N. C. State University

More information

Image Filtering in Spatial domain. Computer Vision Jia-Bin Huang, Virginia Tech

Image Filtering in Spatial domain. Computer Vision Jia-Bin Huang, Virginia Tech Image Filtering in Spatial domain Computer Vision Jia-Bin Huang, Virginia Tech Administrative stuffs Lecture schedule changes Office hours - Jia-Bin (44 Whittemore Hall) Friday at : AM 2: PM Office hours

More information

CS/ECE 545 (Digital Image Processing) Midterm Review

CS/ECE 545 (Digital Image Processing) Midterm Review CS/ECE 545 (Digital Image Processing) Midterm Review Prof Emmanuel Agu Computer Science Dept. Worcester Polytechnic Institute (WPI) Exam Overview Wednesday, March 5, 2014 in class Will cover up to lecture

More information

Achim J. Lilienthal Mobile Robotics and Olfaction Lab, AASS, Örebro University

Achim J. Lilienthal Mobile Robotics and Olfaction Lab, AASS, Örebro University Achim J. Lilienthal Mobile Robotics and Olfaction Lab, Room T1227, Mo, 11-12 o'clock AASS, Örebro University (please drop me an email in advance) achim.lilienthal@oru.se 1 2. General Introduction Schedule

More information

Computer Vision. Intensity transformations

Computer Vision. Intensity transformations Computer Vision Intensity transformations Filippo Bergamasco (filippo.bergamasco@unive.it) http://www.dais.unive.it/~bergamasco DAIS, Ca Foscari University of Venice Academic year 2016/2017 Introduction

More information

Sensors and Sensing Cameras and Camera Calibration

Sensors and Sensing Cameras and Camera Calibration Sensors and Sensing Cameras and Camera Calibration Todor Stoyanov Mobile Robotics and Olfaction Lab Center for Applied Autonomous Sensor Systems Örebro University, Sweden todor.stoyanov@oru.se 20.11.2014

More information

Filtering Images in the Spatial Domain Chapter 3b G&W. Ross Whitaker (modified by Guido Gerig) School of Computing University of Utah

Filtering Images in the Spatial Domain Chapter 3b G&W. Ross Whitaker (modified by Guido Gerig) School of Computing University of Utah Filtering Images in the Spatial Domain Chapter 3b G&W Ross Whitaker (modified by Guido Gerig) School of Computing University of Utah 1 Overview Correlation and convolution Linear filtering Smoothing, kernels,

More information

Image Filtering. Reading Today s Lecture. Reading for Next Time. What would be the result? Some Questions from Last Lecture

Image Filtering. Reading Today s Lecture. Reading for Next Time. What would be the result? Some Questions from Last Lecture Image Filtering HCI/ComS 575X: Computational Perception Instructor: Alexander Stoytchev http://www.cs.iastate.edu/~alex/classes/2007_spring_575x/ January 24, 2007 HCI/ComS 575X: Computational Perception

More information

Midterm is on Thursday!

Midterm is on Thursday! Midterm is on Thursday! Project presentations are May 17th, 22nd and 24th Next week there is a strike on campus. Class is therefore cancelled on Tuesday. Please work on your presentations instead! REVIEW

More information

Image Enhancement in the Spatial Domain (Part 1)

Image Enhancement in the Spatial Domain (Part 1) Image Enhancement in the Spatial Domain (Part 1) Lecturer: Dr. Hossam Hassan Email : hossameldin.hassan@eng.asu.edu.eg Computers and Systems Engineering Principle Objective of Enhancement Process an image

More information

Images and Filters. EE/CSE 576 Linda Shapiro

Images and Filters. EE/CSE 576 Linda Shapiro Images and Filters EE/CSE 576 Linda Shapiro What is an image? 2 3 . We sample the image to get a discrete set of pixels with quantized values. 2. For a gray tone image there is one band F(r,c), with values

More information

Motivation: Image denoising. How can we reduce noise in a photograph?

Motivation: Image denoising. How can we reduce noise in a photograph? Linear filtering Motivation: Image denoising How can we reduce noise in a photograph? Moving average Let s replace each pixel with a weighted average of its neighborhood The weights are called the filter

More information

Image Processing for feature extraction

Image Processing for feature extraction Image Processing for feature extraction 1 Outline Rationale for image pre-processing Gray-scale transformations Geometric transformations Local preprocessing Reading: Sonka et al 5.1, 5.2, 5.3 2 Image

More information

LAB MANUAL SUBJECT: IMAGE PROCESSING BE (COMPUTER) SEM VII

LAB MANUAL SUBJECT: IMAGE PROCESSING BE (COMPUTER) SEM VII LAB MANUAL SUBJECT: IMAGE PROCESSING BE (COMPUTER) SEM VII IMAGE PROCESSING INDEX CLASS: B.E(COMPUTER) SR. NO SEMESTER:VII TITLE OF THE EXPERIMENT. 1 Point processing in spatial domain a. Negation of an

More information

Digital Image Processing. Lecture # 3 Image Enhancement

Digital Image Processing. Lecture # 3 Image Enhancement Digital Image Processing Lecture # 3 Image Enhancement 1 Image Enhancement Image Enhancement 3 Image Enhancement 4 Image Enhancement Process an image so that the result is more suitable than the original

More information

Robert Collins CSE486, Penn State. Lecture 3: Linear Operators

Robert Collins CSE486, Penn State. Lecture 3: Linear Operators Lecture : Linear Operators Administrivia I have put some Matlab image tutorials on Angel. Please take a look if you are unfamiliar with Matlab or the image toolbox. I have posted Homework on Angel. It

More information

Image analysis. CS/CME/BIOPHYS/BMI 279 Fall 2015 Ron Dror

Image analysis. CS/CME/BIOPHYS/BMI 279 Fall 2015 Ron Dror Image analysis CS/CME/BIOPHYS/BMI 279 Fall 2015 Ron Dror A two- dimensional image can be described as a function of two variables f(x,y). For a grayscale image, the value of f(x,y) specifies the brightness

More information

Spatial Domain Processing and Image Enhancement

Spatial Domain Processing and Image Enhancement Spatial Domain Processing and Image Enhancement Lecture 4, Feb 18 th, 2008 Lexing Xie EE4830 Digital Image Processing http://www.ee.columbia.edu/~xlx/ee4830/ thanks to Shahram Ebadollahi and Min Wu for

More information

Chapter 2 Image Enhancement in the Spatial Domain

Chapter 2 Image Enhancement in the Spatial Domain Chapter 2 Image Enhancement in the Spatial Domain Abstract Although the transform domain processing is essential, as the images naturally occur in the spatial domain, image enhancement in the spatial domain

More information

Image filtering, image operations. Jana Kosecka

Image filtering, image operations. Jana Kosecka Image filtering, image operations Jana Kosecka - photometric aspects of image formation - gray level images - point-wise operations - linear filtering Image Brightness values I(x,y) Images Images contain

More information

CEE598 - Visual Sensing for Civil Infrastructure Eng. & Mgmt.

CEE598 - Visual Sensing for Civil Infrastructure Eng. & Mgmt. CEE598 - Visual Sensing for Civil Infrastructure Eng. & Mgmt. Session 7 Pixels and Image Filtering Mani Golparvar-Fard Department of Civil and Environmental Engineering 329D, Newmark Civil Engineering

More information

Image restoration and color image processing

Image restoration and color image processing 1 Enabling Technologies for Sports (5XSF0) Image restoration and color image processing Sveta Zinger ( s.zinger@tue.nl ) What is image restoration? 2 Reconstructing or recovering an image that has been

More information

1.Discuss the frequency domain techniques of image enhancement in detail.

1.Discuss the frequency domain techniques of image enhancement in detail. 1.Discuss the frequency domain techniques of image enhancement in detail. Enhancement In Frequency Domain: The frequency domain methods of image enhancement are based on convolution theorem. This is represented

More information

Last Lecture. Lecture 2, Point Processing GW , & , Ida-Maria Which image is wich channel?

Last Lecture. Lecture 2, Point Processing GW , & , Ida-Maria Which image is wich channel? Last Lecture Lecture 2, Point Processing GW 2.6-2.6.4, & 3.1-3.4, Ida-Maria Ida.sintorn@it.uu.se Digitization -sampling in space (x,y) -sampling in amplitude (intensity) How often should you sample in

More information

Image Processing. 2. Point Processes. Computer Engineering, Sejong University Dongil Han. Spatial domain processing

Image Processing. 2. Point Processes. Computer Engineering, Sejong University Dongil Han. Spatial domain processing Image Processing 2. Point Processes Computer Engineering, Sejong University Dongil Han Spatial domain processing g(x,y) = T[f(x,y)] f(x,y) : input image g(x,y) : processed image T[.] : operator on f, defined

More information

8.2 IMAGE PROCESSING VERSUS IMAGE ANALYSIS Image processing: The collection of routines and

8.2 IMAGE PROCESSING VERSUS IMAGE ANALYSIS Image processing: The collection of routines and 8.1 INTRODUCTION In this chapter, we will study and discuss some fundamental techniques for image processing and image analysis, with a few examples of routines developed for certain purposes. 8.2 IMAGE

More information

Image Enhancement II: Neighborhood Operations

Image Enhancement II: Neighborhood Operations Image Enhancement II: Neighborhood Operations Image Enhancement:Spatial Filtering Operation Idea: Use a mask to alter piel values according to local operation Aim: De)-Emphasize some spatial requencies

More information

A.V.C. COLLEGE OF ENGINEERING DEPARTEMENT OF CSE CP7004- IMAGE PROCESSING AND ANALYSIS UNIT 1- QUESTION BANK

A.V.C. COLLEGE OF ENGINEERING DEPARTEMENT OF CSE CP7004- IMAGE PROCESSING AND ANALYSIS UNIT 1- QUESTION BANK A.V.C. COLLEGE OF ENGINEERING DEPARTEMENT OF CSE CP7004- IMAGE PROCESSING AND ANALYSIS UNIT 1- QUESTION BANK STAFF NAME: TAMILSELVAN K UNIT I SPATIAL DOMAIN PROCESSING Introduction to image processing

More information

Preparing Remote Sensing Data for Natural Resources Mapping (image enhancement, rectifications )

Preparing Remote Sensing Data for Natural Resources Mapping (image enhancement, rectifications ) Preparing Remote Sensing Data for Natural Resources Mapping (image enhancement, rectifications ) Why is this important What are the major approaches Examples of digital image enhancement Follow up exercises

More information

More image filtering , , Computational Photography Fall 2017, Lecture 4

More image filtering , , Computational Photography Fall 2017, Lecture 4 More image filtering http://graphics.cs.cmu.edu/courses/15-463 15-463, 15-663, 15-862 Computational Photography Fall 2017, Lecture 4 Course announcements Any questions about Homework 1? - How many of you

More information

Part I Feature Extraction (1) Image Enhancement. CSc I6716 Spring Local, meaningful, detectable parts of the image.

Part I Feature Extraction (1) Image Enhancement. CSc I6716 Spring Local, meaningful, detectable parts of the image. CSc I6716 Spring 211 Introduction Part I Feature Extraction (1) Zhigang Zhu, City College of New York zhu@cs.ccny.cuny.edu Image Enhancement What are Image Features? Local, meaningful, detectable parts

More information

Lecture 3: Linear Filters

Lecture 3: Linear Filters Signal Denoising Lecture 3: Linear Filters Math 490 Prof. Todd Wittman The Citadel Suppose we have a noisy 1D signal f(x). For example, it could represent a company's stock price over time. In order to

More information

>>> from numpy import random as r >>> I = r.rand(256,256);

>>> from numpy import random as r >>> I = r.rand(256,256); WHAT IS AN IMAGE? >>> from numpy import random as r >>> I = r.rand(256,256); Think-Pair-Share: - What is this? What does it look like? - Which values does it take? - How many values can it take? - Is it

More information

Digital Image Processing. Digital Image Fundamentals II 12 th June, 2017

Digital Image Processing. Digital Image Fundamentals II 12 th June, 2017 Digital Image Processing Digital Image Fundamentals II 12 th June, 2017 Image Enhancement Image Enhancement Types of Image Enhancement Operations Neighborhood Operations on Images Spatial Filtering Filtering

More information

Computer Vision, Lecture 3

Computer Vision, Lecture 3 Computer Vision, Lecture 3 Professor Hager http://www.cs.jhu.edu/~hager /4/200 CS 46, Copyright G.D. Hager Outline for Today Image noise Filtering by Convolution Properties of Convolution /4/200 CS 46,

More information

Sharpening Spatial Filters ( high pass)

Sharpening Spatial Filters ( high pass) Sharpening Spatial Filters ( high pass) Previously we have looked at smoothing filters which remove fine detail Sharpening spatial filters seek to highlight fine detail Remove blurring from images Highlight

More information

Image Processing COS 426

Image Processing COS 426 Image Processing COS 426 What is a Digital Image? A digital image is a discrete array of samples representing a continuous 2D function Continuous function Discrete samples Limitations on Digital Images

More information

Introduction. Computer Vision. CSc I6716 Fall Part I. Image Enhancement. Zhigang Zhu, City College of New York

Introduction. Computer Vision. CSc I6716 Fall Part I. Image Enhancement. Zhigang Zhu, City College of New York CSc I6716 Fall 21 Introduction Part I Feature Extraction ti (1) Zhigang Zhu, City College of New York zhu@cs.ccny.cuny.edu Image Enhancement What are Image Features? Local, meaningful, detectable parts

More information

DIGITAL IMAGE DE-NOISING FILTERS A COMPREHENSIVE STUDY

DIGITAL IMAGE DE-NOISING FILTERS A COMPREHENSIVE STUDY INTERNATIONAL JOURNAL OF RESEARCH IN COMPUTER APPLICATIONS AND ROBOTICS ISSN 2320-7345 DIGITAL IMAGE DE-NOISING FILTERS A COMPREHENSIVE STUDY Jaskaranjit Kaur 1, Ranjeet Kaur 2 1 M.Tech (CSE) Student,

More information

Computing for Engineers in Python

Computing for Engineers in Python Computing for Engineers in Python Lecture 10: Signal (Image) Processing Autumn 2011-12 Some slides incorporated from Benny Chor s course 1 Lecture 9: Highlights Sorting, searching and time complexity Preprocessing

More information

Image Enhancement contd. An example of low pass filters is:

Image Enhancement contd. An example of low pass filters is: Image Enhancement contd. An example of low pass filters is: We saw: unsharp masking is just a method to emphasize high spatial frequencies. We get a similar effect using high pass filters (for instance,

More information

TIRF, geometric operators

TIRF, geometric operators TIRF, geometric operators Last class FRET TIRF This class Finish up of TIRF Geometric image processing TIRF light confinement II(zz) = II 0 ee zz/dd 1 TIRF Intensity for d = 300 nm 0.9 0.8 0.7 0.6 Relative

More information

DIGITAL IMAGE PROCESSING (COM-3371) Week 2 - January 14, 2002

DIGITAL IMAGE PROCESSING (COM-3371) Week 2 - January 14, 2002 DIGITAL IMAGE PROCESSING (COM-3371) Week 2 - January 14, 22 Topics: Human eye Visual phenomena Simple image model Image enhancement Point processes Histogram Lookup tables Contrast compression and stretching

More information

Digital Image Fundamentals and Image Enhancement in the Spatial Domain

Digital Image Fundamentals and Image Enhancement in the Spatial Domain Digital Image Fundamentals and Image Enhancement in the Spatial Domain Mohamed N. Ahmed, Ph.D. Introduction An image may be defined as 2D function f(x,y), where x and y are spatial coordinates. The amplitude

More information

Digital Image Processing 3/e

Digital Image Processing 3/e Laboratory Projects for Digital Image Processing 3/e by Gonzalez and Woods 2008 Prentice Hall Upper Saddle River, NJ 07458 USA www.imageprocessingplace.com The following sample laboratory projects are

More information

Filtering. Image Enhancement Spatial and Frequency Based

Filtering. Image Enhancement Spatial and Frequency Based Filtering Image Enhancement Spatial and Frequency Based Brent M. Dingle, Ph.D. 2015 Game Design and Development Program Mathematics, Statistics and Computer Science University of Wisconsin - Stout Lecture

More information

Computer Graphics Fundamentals

Computer Graphics Fundamentals Computer Graphics Fundamentals Jacek Kęsik, PhD Simple converts Rotations Translations Flips Resizing Geometry Rotation n * 90 degrees other Geometry Rotation n * 90 degrees other Geometry Translations

More information

Image Filtering Josef Pelikán & Alexander Wilkie CGG MFF UK Praha

Image Filtering Josef Pelikán & Alexander Wilkie CGG MFF UK Praha Image Filtering 1995-216 Josef Pelikán & Alexander Wilkie CGG MFF UK Praha pepca@cgg.mff.cuni.cz http://cgg.mff.cuni.cz/~pepca/ 1 / 32 Image Histograms Frequency table of individual brightness (and sometimes

More information

02/02/10. Image Filtering. Computer Vision CS 543 / ECE 549 University of Illinois. Derek Hoiem

02/02/10. Image Filtering. Computer Vision CS 543 / ECE 549 University of Illinois. Derek Hoiem 2/2/ Image Filtering Computer Vision CS 543 / ECE 549 University of Illinois Derek Hoiem Questions about HW? Questions about class? Room change starting thursday: Everitt 63, same time Key ideas from last

More information

Chapter 6. [6]Preprocessing

Chapter 6. [6]Preprocessing Chapter 6 [6]Preprocessing As mentioned in chapter 4, the first stage in the HCR pipeline is preprocessing of the image. We have seen in earlier chapters why this is very important and at the same time

More information

Enhancement Techniques for True Color Images in Spatial Domain

Enhancement Techniques for True Color Images in Spatial Domain Enhancement Techniques for True Color Images in Spatial Domain 1 I. Suneetha, 2 Dr. T. Venkateswarlu 1 Dept. of ECE, AITS, Tirupati, India 2 Dept. of ECE, S.V.University College of Engineering, Tirupati,

More information

Digital Image Processing. Lecture 5 (Enhancement) Bu-Ali Sina University Computer Engineering Dep. Fall 2009

Digital Image Processing. Lecture 5 (Enhancement) Bu-Ali Sina University Computer Engineering Dep. Fall 2009 Digital Image Processing Lecture 5 (Enhancement) Bu-Ali Sina University Computer Engineering Dep. Fall 2009 Outline Image Enhancement in Spatial Domain Histogram based methods Histogram Equalization Local

More information

To process an image so that the result is more suitable than the original image for a specific application.

To process an image so that the result is more suitable than the original image for a specific application. by Shahid Farid 1 To process an image so that the result is more suitable than the original image for a specific application. Categories: Spatial domain methods and Frequency domain methods 2 Procedures

More information

Announcements. Image Processing. What s an image? Images as functions. Image processing. What s a digital image?

Announcements. Image Processing. What s an image? Images as functions. Image processing. What s a digital image? Image Processing Images by Pawan Sinha Today s readings Forsyth & Ponce, chapters 8.-8. http://www.cs.washington.edu/education/courses/49cv/wi/readings/book-7-revised-a-indx.pdf For Monday Watt,.3-.4 (handout)

More information

DIGITAL IMAGE PROCESSING Quiz exercises preparation for the midterm exam

DIGITAL IMAGE PROCESSING Quiz exercises preparation for the midterm exam DIGITAL IMAGE PROCESSING Quiz exercises preparation for the midterm exam In the following set of questions, there are, possibly, multiple correct answers (1, 2, 3 or 4). Mark the answers you consider correct.

More information

INSTITUTE OF AERONAUTICAL ENGINEERING Dundigal, Hyderabad

INSTITUTE OF AERONAUTICAL ENGINEERING Dundigal, Hyderabad INSTITUTE OF AERONAUTICAL ENGINEERING Dundigal, Hyderabad - 500 043 ELECTRONICS AND COMMUNICATION ENGINEERING QUESTION BANK Course Title Course Code Class Branch DIGITAL IMAGE PROCESSING A70436 IV B. Tech.

More information

Subband coring for image noise reduction. Edward H. Adelson Internal Report, RCA David Sarnoff Research Center, Nov

Subband coring for image noise reduction. Edward H. Adelson Internal Report, RCA David Sarnoff Research Center, Nov Subband coring for image noise reduction. dward H. Adelson Internal Report, RCA David Sarnoff Research Center, Nov. 26 1986. Let an image consisting of the array of pixels, (x,y), be denoted (the boldface

More information

What is an image? Bernd Girod: EE368 Digital Image Processing Pixel Operations no. 1. A digital image can be written as a matrix

What is an image? Bernd Girod: EE368 Digital Image Processing Pixel Operations no. 1. A digital image can be written as a matrix What is an image? Definition: An image is a 2-dimensional light intensity function, f(x,y), where x and y are spatial coordinates, and f at (x,y) is related to the brightness of the image at that point.

More information

CS 4501: Introduction to Computer Vision. Filtering and Edge Detection

CS 4501: Introduction to Computer Vision. Filtering and Edge Detection CS 451: Introduction to Computer Vision Filtering and Edge Detection Connelly Barnes Slides from Jason Lawrence, Fei Fei Li, Juan Carlos Niebles, Misha Kazhdan, Allison Klein, Tom Funkhouser, Adam Finkelstein,

More information

Image analysis. CS/CME/BioE/Biophys/BMI 279 Oct. 31 and Nov. 2, 2017 Ron Dror

Image analysis. CS/CME/BioE/Biophys/BMI 279 Oct. 31 and Nov. 2, 2017 Ron Dror Image analysis CS/CME/BioE/Biophys/BMI 279 Oct. 31 and Nov. 2, 2017 Ron Dror 1 Outline Images in molecular and cellular biology Reducing image noise Mean and Gaussian filters Frequency domain interpretation

More information

Midterm Examination CS 534: Computational Photography

Midterm Examination CS 534: Computational Photography Midterm Examination CS 534: Computational Photography November 3, 2015 NAME: SOLUTIONS Problem Score Max Score 1 8 2 8 3 9 4 4 5 3 6 4 7 6 8 13 9 7 10 4 11 7 12 10 13 9 14 8 Total 100 1 1. [8] What are

More information

Analysis of infrared images in integrated-circuit techniques by mathematical filtering

Analysis of infrared images in integrated-circuit techniques by mathematical filtering 10 th International Conference on Quantitative InfraRed Thermography July 27-30, 2010, Québec (Canada) Analysis of infrared images in integrated-circuit techniques by mathematical filtering by I. Benkö

More information

ECC419 IMAGE PROCESSING

ECC419 IMAGE PROCESSING ECC419 IMAGE PROCESSING INTRODUCTION Image Processing Image processing is a subclass of signal processing concerned specifically with pictures. Digital Image Processing, process digital images by means

More information

CSE 564: Visualization. Image Operations. Motivation. Provide the user (scientist, t doctor, ) with some means to: Global operations:

CSE 564: Visualization. Image Operations. Motivation. Provide the user (scientist, t doctor, ) with some means to: Global operations: Motivation CSE 564: Visualization mage Operations Klaus Mueller Computer Science Department Stony Brook University Provide the user (scientist, t doctor, ) with some means to: enhance contrast of local

More information

Chapter 3. Study and Analysis of Different Noise Reduction Filters

Chapter 3. Study and Analysis of Different Noise Reduction Filters Chapter 3 Study and Analysis of Different Noise Reduction Filters Noise is considered to be any measurement that is not part of the phenomena of interest. Departure of ideal signal is generally referred

More information

Image analysis. CS/CME/BioE/Biophys/BMI 279 Oct. 31 and Nov. 2, 2017 Ron Dror

Image analysis. CS/CME/BioE/Biophys/BMI 279 Oct. 31 and Nov. 2, 2017 Ron Dror Image analysis CS/CME/BioE/Biophys/BMI 279 Oct. 31 and Nov. 2, 2017 Ron Dror 1 Outline Images in molecular and cellular biology Reducing image noise Mean and Gaussian filters Frequency domain interpretation

More information

Image Enhancement. Image Enhancement

Image Enhancement. Image Enhancement SPATIAL FILTERING g h * h g FREQUENCY DOMAIN FILTERING G H. F F H G Copright RMR / RDL - 999. PEE53 - Processamento Digital de Imagens LOW PASS FILTERING attenuate or eliminate high-requenc components

More information

Color Space 1: RGB Color Space. Color Space 2: HSV. RGB Cube Easy for devices But not perceptual Where do the grays live? Where is hue and saturation?

Color Space 1: RGB Color Space. Color Space 2: HSV. RGB Cube Easy for devices But not perceptual Where do the grays live? Where is hue and saturation? Color Space : RGB Color Space Color Space 2: HSV RGB Cube Easy for devices But not perceptual Where do the grays live? Where is hue and saturation? Hue, Saturation, Value (Intensity) RBG cube on its vertex

More information

SECTION I - CHAPTER 2 DIGITAL IMAGING PROCESSING CONCEPTS

SECTION I - CHAPTER 2 DIGITAL IMAGING PROCESSING CONCEPTS RADT 3463 - COMPUTERIZED IMAGING Section I: Chapter 2 RADT 3463 Computerized Imaging 1 SECTION I - CHAPTER 2 DIGITAL IMAGING PROCESSING CONCEPTS RADT 3463 COMPUTERIZED IMAGING Section I: Chapter 2 RADT

More information

Digital Image Processing. Lecture # 4 Image Enhancement (Histogram)

Digital Image Processing. Lecture # 4 Image Enhancement (Histogram) Digital Image Processing Lecture # 4 Image Enhancement (Histogram) 1 Histogram of a Grayscale Image Let I be a 1-band (grayscale) image. I(r,c) is an 8-bit integer between 0 and 255. Histogram, h I, of

More information

A Study On Preprocessing A Mammogram Image Using Adaptive Median Filter

A Study On Preprocessing A Mammogram Image Using Adaptive Median Filter A Study On Preprocessing A Mammogram Image Using Adaptive Median Filter Dr.K.Meenakshi Sundaram 1, D.Sasikala 2, P.Aarthi Rani 3 Associate Professor, Department of Computer Science, Erode Arts and Science

More information

Denoising and Effective Contrast Enhancement for Dynamic Range Mapping

Denoising and Effective Contrast Enhancement for Dynamic Range Mapping Denoising and Effective Contrast Enhancement for Dynamic Range Mapping G. Kiruthiga Department of Electronics and Communication Adithya Institute of Technology Coimbatore B. Hakkem Department of Electronics

More information