20111ee114_1_CA4_sol

20111ee114_1_CA4_sol - x_i – 1, y_i – 1]); return;...

Info iconThis preview shows pages 1–4. Sign up to view the full content.

View Full Document Right Arrow Icon
UCLA, Electrical Engineering Department Winter Quarter 2011 EE114, Part 2: Image Processing Solutions to Computer assignment #4 1. m_files We need to write a couple of m_files. GenerateGaussian.m function h = GenerateGaussian (size, lambda) % this function generates a Gaussian matrix % generate the meshgrid and the Gaussian matrix [x,y]= meshgrid([- size / 2 : 1 : size / 2]); h = exp( - (x .* x + y .*y) / lambda); % find the sum and normalize the matrix sum_total = sum(sum(h)); h = h ./ sum_total; return; ConvolveImages.m function out = ConvolveImages (image, h) % this function convolves two images and % crop the result to the size of the first one % get the dimensions of the input images [y_i, x_i] = size(image); [y_h, x_h] = size (h); % convolve the images out = uint8(conv2(image, h)); %crop it to the correct size out = imcrop(out, [(x_h + 1)/ 2, (y_h + 1) / 2,
Background image of page 1

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon
Background image of page 2
Background image of page 3

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon
Background image of page 4
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: x_i – 1, y_i – 1]); return; PSNR.m function psnr = PSNR (i1, i2) % calculate PSNR of two images d = double(i1) – double(i2) ; d2 = d .* d; psnr = 10 * log10 (255 * 255/ mean(mean(d2))); return; 2. Results Under matlab, type image = imread(‘image.tif’); imagesc(image); colormap (gray); to load and display image. Type h5 = GenerateGaussian(40, 5); i5 = ConvolveImages(image, h5); imagesc(i5); PSNR5 = PSNR(image, i5); to generate Gaussian matrix with λ = 5, convolve image with it, display the result and calculate the PSNR against the original. PSNR5 = 25.3865 Type h10 = GenerateGaussian(40, 10); i10 = ConvolveImages(image, h10); imagesc(i10); PSNR10 = PSNR(image, i10); to generate Gaussian matrix with λ = 10, convolve image with it, display the result and calculate the PSNR against the original. PSNR10 = 23.1213...
View Full Document

This note was uploaded on 05/21/2011 for the course EE 114 taught by Professor Vanschaar during the Spring '11 term at UCLA.

Page1 / 4

20111ee114_1_CA4_sol - x_i – 1, y_i – 1]); return;...

This preview shows document pages 1 - 4. Sign up to view the full document.

View Full Document Right Arrow Icon
Ask a homework question - tutors are online