MATLAB: I am encountering this error on below mention code Undefined function or variable ‘greenchannel’. Error in bidemo (line 47) [Gout]=bilateral_each_channel(w,sigma_r,sigma_d,greenchannel);

bilateralMATLAB

%%Clear everything
clear all;
close all;
clc;
v = VideoReader('peng.mp4');
tic
numFrames = v.NumberOfFrames;
k=1;
while k<=numFrames
  frame=read(v,k);
  frame=imresize(frame,[256 256]);
  frame=double(frame);
  [a b z]=size(frame);
    frame = frame/255.0;
% Extract the individual red, green, and blue color channels.
redChannel = frame(:, :, 1);
greenChannel = frame(:, :, 2);
blueChannel = frame(:, :, 3);
sigma = 0.1; %noise standard deviation
xr = redChannel + sigma*randn(size(redChannel));
xg = greenChannel + sigma*randn(size(greenChannel));
xb = blueChannel + sigma*randn(size(blueChannel));
noisyRGB = cat(3,xr,xg,xb);
w=9; % window size
sigma_d=1;
sigma_r=30;
tic;
[Rout]=bilateral_each_channel(w,sigma_r,sigma_d,redChannel);
[Gout]=bilateral_each_channel(w,sigma_r,sigma_d,greenchannel);
[Bout]=bilateral_each_channel(w,sigma_r,sigma_d,bluechannel);
frame(:,:,1)=(Rout);
frame(:,:,2)=(Gout);
frame(:,:,3)=(Bout);
rgbFixed = cat(3, redMF,greenMF,blueMF);
imshow([frame,noisyRGB,rgbFixed],[0,1])
  k=k+1;
end
finalpsnr1=avgPSNR1/numFrames
% finalpsnr2=avgPSNR2/numFrames
toc

Best Answer

greenchannel is not the same as greenChannel

Related Solutions

MATLAB: Reading video file to Matlab

See if this demo helps you any:

% Demo macro to extract frames and get frame means from an avi movie
% and save individual frames to separate image files.
% Also computes the mean gray value of the color channels.
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 = 14;
% Change the current folder to the folder of this m-file.
% (The line of code below is from Brett Shoelson of The Mathworks.)
if(~isdeployed)
  cd(fileparts(which(mfilename)));
end
% Check that user has the Image Processing Toolbox installed.
hasIPT = license('test', 'image_toolbox');
if ~hasIPT
  % User does not have the toolbox installed.
  message = sprintf('Sorry, but you do not seem to have the Image Processing Toolbox.\nDo you want to try to continue anyway?');
  reply = questdlg(message, 'Toolbox missing', 'Yes', 'No', 'Yes');
  if strcmpi(reply, 'No')
    % User said No, so exit.
    return;
  end
end
% Open the rhino.avi demo movie that ships with MATLAB.
folder = fullfile(matlabroot, '\toolbox\images\imdemos');
movieFullFileName = fullfile(folder, 'rhinos.avi');
% Check to see that it exists.
if ~exist(movieFullFileName, 'file')
  strErrorMessage = sprintf('File not found:\n%s\nYou can choose a new one, or cancel', movieFullFileName);
  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)
      movieFullFileName = fullfile(folderName, baseFileName);
    else
      return;
    end
  else
    return;
  end
end
try
  videoObject = VideoReader(movieFullFileName)
  % Determine how many frames there are.
  numberOfFrames = videoObject.NumberOfFrames;
  vidHeight = videoObject.Height;
  vidWidth = videoObject.Width;
  
  numberOfFramesWritten = 0;
  % Prepare a figure to show the images in the upper half of the screen.
  figure;
  %   screenSize = get(0, 'ScreenSize');
  % Enlarge figure to full screen.
  set(gcf, 'units','normalized','outerposition',[0 0 1 1]);
  
  % Loop through the movie, writing all frames out.
  % Each frame will be in a separate file with unique name.
  meanGrayLevels = zeros(numberOfFrames, 1);
  meanRedLevels = zeros(numberOfFrames, 1);
  meanGreenLevels = zeros(numberOfFrames, 1);
  meanBlueLevels = zeros(numberOfFrames, 1);
  for frame = 1 : numberOfFrames
    % Extract the frame from the movie structure.
    thisFrame = read(videoObject, frame);
    
    % Display it
    hImage = subplot(2, 2, 1);
    image(thisFrame);
    caption = sprintf('Frame %4d of %d.', frame, numberOfFrames);
    title(caption, 'FontSize', fontSize);
    drawnow; % Force it to refresh the window.
    
    % Calculate the mean gray level.
    grayImage = rgb2gray(thisFrame);
    meanGrayLevels(frame) = mean(grayImage(:));
    
    % Calculate the mean R, G, and B levels.
    meanRedLevels(frame) = mean(mean(thisFrame(:, :, 1)));
    meanGreenLevels(frame) = mean(mean(thisFrame(:, :, 2)));
    meanBlueLevels(frame) = mean(mean(thisFrame(:, :, 3)));
    
    % Plot the mean gray levels.
    hPlot = subplot(2, 2, 2);
    hold off;
    plot(meanGrayLevels, 'k-', 'LineWidth', 2);
    hold on;
    plot(meanRedLevels, 'r-');
    plot(meanGreenLevels, 'g-');
    plot(meanBlueLevels, 'b-');
    grid on;
    
    % Put title back because plot() erases the existing title.
    title('Mean Gray Levels', 'FontSize', fontSize);
    if frame == 1
      xlabel('Frame Number');
      ylabel('Gray Level');
      % Get size data later for preallocation if we read
      % the movie back in from disk.
      [rows columns numberOfColorChannels] = size(thisFrame);
    end
    
    % 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;
    
    % Now let's do the entropy filter.
    % First get a gray level image
    grayImage = rgb2gray(thisFrame);
    subplot(2, 2, 3);
    imshow(grayImage);
    caption = sprintf('Gray Scale Image, frame #%d', frame);
    title(caption, 'FontSize', fontSize);
    % Display the filtered image.
    entropyFilteredImage = entropyfilt(grayImage);
    subplot(2, 2, 4);
    imshow(entropyFilteredImage, []);
    caption = sprintf('Entropy Filtered Image, frame #%d', frame);
    title(caption, 'FontSize', fontSize);
  end
  
  % Alert user that we're done.
  finishedMessage = sprintf('Done!  It processed %d frames of\n"%s"', numberOfFramesWritten, movieFullFileName);
  disp(finishedMessage); % Write to command window.
  uiwait(msgbox(finishedMessage)); % Also pop up a message box.
  
catch ME
  % Some error happened if you get here.
  stError = lasterror;
  strErrorMessage = sprintf('Error extracting movie frames from:\n\n%s\n\nError: %s\n\n)', movieFullFileName, stError.message);
  uiwait(msgbox(strErrorMessage));
end

MATLAB: Read three different RGB band and Swap the Band

Try this:

% Extract the individual red, green, and blue color channels.
redChannel = rgbImage(:, :, 1);
greenChannel = rgbImage(:, :, 2);
blueChannel = rgbImage(:, :, 3);
% Now swap around in any way you want.
newRedChannel = greenChannel;
newGreenChannel = blueChannel; 
newBlueChannel = redChannel;
% Recombine separate color channels into a single, true color RGB image.
NewRgbImage = cat(3, newRedChannel, newGreenChannel, newBlueChannel);

Best Answer

Related Solutions

MATLAB: Reading video file to Matlab

MATLAB: Read three different RGB band and Swap the Band

Related Question