1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
|
function blurshader = Create2DGaussianBlurShader(stddev, width)
% blurshader = Create2DGaussianBlurShader(stddev, width)
%
% Create a GLSL shader which uses a 2D gaussian filter kernel to
% apply gaussian blur onto a texture. Returns the handle 'blurshader',
% which can be used with Screen('DrawTexture') to apply the shader during
% texture drawing.
%
% stddev = Standard deviation of kernel to use. Defaults to 2.5.
% width = Width of filter-kernel. Currently only a width of 5 is
% supported, so this argument is optional and defaults to 5.
%
% Currently you will need the Matlab imageprocessing toolbox for this to
% work, or the Octave image package, as we use the toolbox function 'fspecial'
% to compute the filter kernel.
%
% Example of use:
%
% Somewhere during initialization of your script you call:
% blurshader = Create2DGaussianBlurShader(2.5);
%
% When you want to draw an image blurred with a gaussian lowpass, you do:
%
% Screen('DrawTexture', win, mytexture, [], [], [], [], [], [], blurshader);
%
% History:
% 15.04.2006 written by Mario Kleiner.
if nargin < 1
stddev = 2.5;
end;
if nargin < 2 || isempty(width)
width = 11;
end;
if width~=11
% error('Create2DGaussianBlurShader: Invalid width specified. Currently only width=5 supported.');
end;
% Build and initialize gaussian blur shader (5x5, stddev=2.5):
try
% Use the imageprocessing toolbox function fspecial for creation of
% kernel:
if IsOctave
pkg load image;
end
kernel = fspecial('gaussian', width, stddev);
catch
error('Create2DGaussianBlurShader failed: You do not seem to have the image processing toolbox installed.');
end;
% Create shader from convolution kernel and return handle to it:
blurshader = EXPCreateStatic2DConvolutionShader(kernel);
return;
|
