Convolution Script for BIOL 319

convolution.m

clear all; close all;
A = imread('BrownDischer079.tiff');
A = single(A);
scale = max(max(A))/256;
A = A/scale;
figure(1);
co=[0:255]/256; colormap([co;co;co]');
image(A); axis off tight; title('original image');

[x,y] = meshgrid(-25:25);
gFilter=exp(-.5*(x.^2/2+y.^2/45));
laplaceFilter=[0 -1 0; -1 4 -1; 0 -1 0];
K=conv2(laplaceFilter,gFilter);
combinedFilter=K(2:(end-1),2:(end-1));
filteredImage = conv2(combinedFilter, A);
figure(2);
colormap([co;co;co]');
image(filteredImage); axis off tight; title('filtered image');