Yes good idea but you have a number of errors in there. Mainly that your histogram equalization is not inside your loop over all frames. See my attached demo.
% Demo to extract frames from a input video, and build a new movie where each frame is histogram equalized.
% Illustrates the use of the VideoReader and VideoWriter classes.
% A Mathworks demo (different than mine) is located here http://www.mathworks.com/help/matlab/examples/convert-between-image-sequences-and-video.html
%-------------------------------------------------------------------------------------------------------------------------------------------
% Initialization steps
clc; % Clear the command window.
close all; % Close all figures (except those of imtool.)
imtool close all; % Close all imtool figures.
clear; % Erase all existing variables.
workspace; % Make sure the workspace panel is showing.
fontSize = 15;
%-------------------------------------------------------------------------------------------------------------------------------------------
% Open the rhino.avi demo movie that ships with MATLAB.
% First get the folder that it lives in.
folder = fileparts(which('rhinos.avi')); % Determine where demo folder is (works with all versions).
% Pick one of the two demo movies shipped with the Image Processing Toolbox.
% Comment out the other one.
inputFullFileName = fullfile(folder, 'rhinos.avi');
% movieFullFileName = fullfile(folder, 'traffic.avi');
% Check to see that it exists.
if ~exist(inputFullFileName, 'file')
strErrorMessage = sprintf('File not found:\n%s\nYou can choose a new one, or cancel', inputFullFileName);
response = questdlg(strErrorMessage, 'File not found', 'OK - choose a new movie.', 'Cancel', 'OK - choose a new movie.');
if strcmpi(response, 'OK - choose a new movie.')
[baseFileName, folderName, FilterIndex] = uigetfile('*.avi');
if ~isequal(baseFileName, 0)
inputFullFileName = fullfile(folderName, baseFileName);
else
return;
end
else
return;
end
end
%-------------------------------------------------------------------------------------------------------------------------------------------
% Open up the VideoReader for reading an input video file.
inputVideoReaderObject = VideoReader(inputFullFileName)
% Determine how many frames there are.
numberOfFrames = inputVideoReaderObject.NumberOfFrames;
inputVideoRows = inputVideoReaderObject.Height
inputVideoColumns = inputVideoReaderObject.Width
%-------------------------------------------------------------------------------------------------------------------------------------------
% Create an output VideoWriter object to write the video out to a new, different file.
outputFullFileName = fullfile(pwd, 'NewRhinos.avi');
outputVideoWriterObject = VideoWriter(outputFullFileName);
open(outputVideoWriterObject);
% Specify the output video size.
shrinkFactor = 4; % Shrink by a factor of 4 in both directions.
outputVideoRows = round(inputVideoRows / shrinkFactor)
outputVideoColumns = round(inputVideoColumns / shrinkFactor)
numberOfFramesWritten = 0;
% Prepare a figure to show the images.
hFig = figure('Name', 'Demo by Image Analyst', 'NumberTitle', 'off');
%-------------------------------------------------------------------------------------------------------------------------------------------
% Loop through the movie, writing all frames out.
for frame = 1 : numberOfFrames
% Extract the frame from the movie structure.
thisInputFrame = read(inputVideoReaderObject, frame);
% Histogram equalize the image.
% First convert the image to LAB color space because to make it look natural
% we want to only filter the overall brightness, not each individual color channel.
LAB = rgb2lab(thisInputFrame);
L = LAB(:, :, 1); % Get the brightness or luminance channel.
% L is floating point in the range 0-100. If we are to use adapthisteq,
% it must be in the range 0-1.
L = L / 100;
% Do the histogram equalization on this frame.
L = adapthisteq(L, 'NumTiles', [8 8], 'ClipLimit', 0.005);
% Undo the scaling to get it back into its original range.
L = L * 100;
% Put the filtered result back into the LAB array.
LAB(:, :, 1) = L;
% Turn it from LAB color space back into RGB colorspace.
outputFrame = 255 * lab2rgb(LAB);
% Make it uint8 so we can put side by side with the original frame.
outputFrame = uint8(outputFrame);
imshowpair(thisInputFrame, outputFrame,'montage');
progressIndication = sprintf('Processed frame %4d of %d.', frame, numberOfFrames);
title(progressIndication, 'FontSize', fontSize);
drawnow;
% Write this new, resized frame out to the new video file.
writeVideo(outputVideoWriterObject, outputFrame);
% Update user with the progress. Display in the command window.
progressIndication = sprintf('Processed frame %4d of %d.', frame, numberOfFrames);
disp(progressIndication);
% Increment frame count (should eventually = numberOfFrames
% unless an error happens).
numberOfFramesWritten = numberOfFramesWritten + 1;
end
% Close the output video object. You don't need to close the input video reader.
close(outputVideoWriterObject);
%================================================================================
% NOW WE'RE DONE!
%================================================================================
% Alert user that we're done.
finishedMessage = sprintf('Done! It processed %d frames of\n"%s" and created output video\n%s.\nClick OK to see the output video.', numberOfFramesWritten, inputFullFileName, outputFullFileName);
fprintf('%s\n', finishedMessage); % Write to command window.
uiwait(msgbox(finishedMessage)); % Also pop up a message box.
% OPTIONAL: Play the movie in the default Windows program that plays videos.
winopen(outputFullFileName);
