Subversion Repositories fpga-median

% +----------------------------------------------------------------------------
% Universidade Federal da Bahia
% ------------------------------------------------------------------------------
% PROJECT: FPGA Median Filter
% ------------------------------------------------------------------------------
% FILE NAME            : median_ref.m
% AUTHOR               : João Carlos Bittencourt
% AUTHOR'S E-MAIL      : joaocarlos@ieee.org
% -----------------------------------------------------------------------------
% RELEASE HISTORY
% VERSION  DATE        AUTHOR        DESCRIPTION
% 1.0      2013-08-27  joao.nunes    initial version
% -----------------------------------------------------------------------------
% KEYWORDS: median, filter, image processing
% -----------------------------------------------------------------------------
% PURPOSE: Reference model for Median Filter.
% -----------------------------------------------------------------------------
 
% This is only a simple verification reference model based on default Median Filter
 
% Clear environment
clc
clear
% Set-up vectors
% I have set the vector to size 12 in order to perform a basic testbench.
% result = zeros(227,227);
img_ref = imread('images/image22.jpg');
 
img_ref = rgb2gray(img_ref);
[width, height] = size(img_ref);
% Add one column with zeros for pipelining verification purpose.
% The hardware assumes that, in the first round in a row, previous values are zeros.
img = zeros(width+1,width+1);
for i = 1 : height,
     img(i,2:width+1) = img_ref(i, 1:width);
end
 
 
% Default Median Algorithm
window_width = 3;
window_height = 3;
edgex = floor(window_width/2);
edgey = floor(window_height/2);
tic
for x = edgex : width - edgex,
    for y = edgey : height - edgey,
        temp = zeros(edgex,edgey);
        for fx = 1 : window_width,
            for fy = 1 : window_height,
                temp(fx,fy) = img(x + fx - edgex, y + fy - edgey);
            end       
        end
        temp = reshape(temp.',1,[]);
        srt = sort(temp); % remove comma to view step by step outputs
        result(x,y) = srt(5);
    end
end
toc
wtime = toc
fprintf ( 1, '  MY_PROGRAM took %f seconds to run.\n', wtime );
 
%result
imshow(mat2gray(result));
imwrite(mat2gray(result), 'images/image22_median.jpg', 'jpg');