%?function?exampleFilter
%
%?Example?of?how?the?Kalman?filter?performs?for?real?data
obj=VideoReader（‘julius.avi‘）;

aviobj=VideoWriter（‘outputvideo.avi‘）
open（aviobj）
nframes=obj.NumberOfframes;
for?k=1:nframes
????img=read（objk）;
????im2?=?imnoise（img‘gaussian‘0.10.2）
????f=im2frame（im2）
????writeVideo（aviobjf）
end
????
close（aviobj）
%--------------------------------------------------------------------------
%---------------------End?of?Gaussian?noise?adding-------------------------
%--------------------------------------------------------------------------
%implay（‘outputvideo.avi‘）;

filename?=?‘outputvideo.avi‘;
mov?=?VideoReader（filename）;

%?Output?folder
%getting?no?of?frames
videoframes=[];
newframes=[];
realframe=[];
wienerframe=[];
psnrk=[];
filteredframe=[];
filteredframe1=[];
numberOfframes?=?mov.NumberOfframes;
numberOfframesWritten?=?0;

for?frame?=?1?:numberOfframes
???thisframe?=?read（mov?frame）;
??%?outputbaseFileName?=?sprintf（‘%3.3d.bmp‘?frame）;
???
???%outputFullFileName?=?fullfile（outputFolder?outputbaseFileName）;
???%imwrite（thisframe?outputFullFileName?‘bmp‘）;
???%progressIndication?=?sprintf（‘Wrote?frame?%4d?of?%d.‘?framenumberOfframes）;
???videoframes?=?cat（3videoframesthisframe）;
???%disp（progressIndication）;
???%numberOfframesWritten?=?numberOfframesWritten?+?1;
end
for?frame=1:size（videoframes3）
????x=videoframes（::frame）;
????y=im2single（x）;
????newframes?=?cat（3newframesy）;
end

%Apply?filter
k=Kalman_Stack_Filter（newframes）;
k75=Kalman_Stack_Filter（newframes0.75）;
outputFolder=?fullfile（cd?‘filteredframe‘）;
if?~exist（outputFolder?‘dir‘）
???mkdir（outputFolder）;
end
y=size（k3）;
disp（y）;
for?frame=1:y
????thisframe=k（::frame）;
????thisframe1=k75（::frame）;
????y=im2single（thisframe）;
????y1=im2single（thisframe1）;
????filteredframe?=?cat（3filteredframey）;
????filteredframe1=cat（3filteredframe1y1）;

????%outputbaseFileName1?=?sprintf（‘%3.3d.bmp‘?frame）;
????%disp（x）;
???%disp（outputbaseFileName1）;
????%outputFullFileName1?=?fullfile（outputFolder?outputbaseFileName1）;
????%disp（outputFullFileName1）;
???%disp（y）;
???%imwrite（thisframe?outputFullFileName1?‘bmp‘）;
end
%--------------------------------------------------------------------------
%---------------------Wiener?Filering--------------------------------------
%--------------------------------------------------------------------------
%implay（videoframes）;

outputFolder3?=?fullfile（cd?‘wienerframes‘）;
if?~exist（outputFolder3?‘dir‘）
???mkdir（outputFolder3）;
end
disp（outputFolder3）;

for?frame=1:size（newframes3）
???I=videoframes（::frame）;
???I=im2single（I）;
????%x=frame+5;
????%outputbaseFileName1?=?sprintf（‘%3.3d.bmp‘x）;
????%outputFullFileName1?=?fullfile（outputFolder3?outputbaseFileName1）;
????%disp（outputFullFileName1）;
????J=wiener2（I[5?5]）;
????wienerframe=cat（3wienerframeJ）;
end
%--------------------------------------------------------------------------
%---------------------End?of?Wiener?Filering---------------------

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2016-09-12?16:53??Comparing-and-combining-Kalman-and-Wiener-Filter-for-video-denoising-master\
?????文件??????613275??2016-09-12?16:53??Comparing-and-combining-Kalman-and-Wiener-Filter-for-video-denoising-master\FinalYearProjectPaper-2.docx
?????文件????????1088??2016-09-12?16:53??Comparing-and-combining-Kalman-and-Wiener-Filter-for-video-denoising-master\Kalman_Stack_Filter.m
?????文件?????????742??2016-09-12?16:53??Comparing-and-combining-Kalman-and-Wiener-Filter-for-video-denoising-master\README.md
?????文件????????6065??2016-09-12?16:53??Comparing-and-combining-Kalman-and-Wiener-Filter-for-video-denoising-master\exampleFilter.m
?????文件????????4035??2016-09-12?16:53??Comparing-and-combining-Kalman-and-Wiener-Filter-for-video-denoising-master\gaussiannoise.m
?????文件??????820224??2016-09-12?16:53??Comparing-and-combining-Kalman-and-Wiener-Filter-for-video-denoising-master\julius.avi
?????文件?????????370??2016-09-12?16:53??Comparing-and-combining-Kalman-and-Wiener-Filter-for-video-denoising-master\psnrcalculationnew.m
?????文件?????????532??2016-09-12?16:53??Comparing-and-combining-Kalman-and-Wiener-Filter-for-video-denoising-master\terminal.m
?????文件?????????249??2016-09-12?16:53??Comparing-and-combining-Kalman-and-Wiener-Filter-for-video-denoising-master\wienerfilter.m