Download All

Code covered by the BSD License  

Highlights from
Free-knot spline approximation

image thumbnail
from Free-knot spline approximation by Bruno Luong
Least squares approximation of 1D data using free-knots spline

KramerBruckner(A, window) 
function B = KramerBruckner(A, window)
%
% function B = KramerBruckner(A, window)
%
% Kramer & Bruckner filtering
%
% AUTHOR: Bruno Luong <brunoluong@yahoo.com>
% HISTORY
%   Original: 12-Jul-2009
%

if nargin<2
    window = [];
end

% cast
classA = class(A);

% Min/Max filtering
A = double(A);
[minimg maximg] = minmaxfilt(A, window, 'both', 'same');

% Kramer & Brunckner
meanimg = (minimg+maximg)/2;
B = minimg;
B(A>=meanimg) = maximg(A>=meanimg);

if ~strcmp(classA,'double')
    B = feval(classA,B);
end

Contact us at files@mathworks.com