% wiener filtering problem
load hw7image;
[sx sy] = size(ho);
% adding noise
hon = ho + 10*randn([sx sy]);

% autocorrelation and noise estimates

% derive and apply wiener filter

% determine MSE's

% work on 48x48 blocks
blkfiltim = zeros([sx sy]);
for bx = 1:15
    for by = 1:10
        subim = hon((by-1).*48+[1:48],(bx-1).*48+[1:48]);
        % do filtering here
        subimfilt = subim;  % place holder
        blkfiltim((by-1).*48+[1:48],(bx-1).*48+[1:48]) = subimfilt;
    end
end

% determine MSE's
