rkttest.m

%% This is a script for displaying
%% rotational kernel transformation (RKT)
%% in comparison to the median filter.
%% (Parameters are arbitrary)
%%
%% Requires to run:
%% - from Image processing toolbox: medfilt2
%% - following files:   rktfilt.m       ..the function for filtering
%%                      rktkern.m       ..the function for creating a kernel
%%                      rktexample.mat  ..the sample image

load ('rktexample.mat')

figure
subplot(1,3,1)
imagesc(example);
title('Original image');

subplot(1,3,2)
imagesc(medfilt2(example, [5 5]));
title('Median filtered image');

subplot(1,3,3)
imagesc(rktfilt(example, rktkern(9)));
title('RKT filtered image');