How to remove pixels in an image

Hi,

I need to remove pixels with zero pixel value in an image(dead pixels)

I have to fill the pixel value with zero with the average of the neighbourhood pixel values.

Can anyone suggest me something about this

Thanks in advance

Answer 1

[Image Analyst]

Detect the zero pixels to create a mask. Then use conv2() to get the average. Then make an output image the same as your input image. Then assign the "zero pixel locations" of your output image to the zero pixel locations of your convolved image. Is that good enough explanation for you to figure it out, or do you require the full code (which I don't have time to give you right now)?

[Jim]

Hi,

Thank you for ur reply

How to detect the dead pixels in an image

with what pixel value I have to create mask

Thanks in advance

[Image Analyst]

How about

deadPixelMap = grayImage == 0;

[Sean de Wolski]

Very good explanation IA, exactly the first thing to come to my mind. Convolve the image with a kernel, K, such that, K = ones(3)./8;K(5)=0; and then replace the dead pixels in each. zombies = ~I; (I(zombies) = Iconved(zombies);

[Jim]

Hi,

I find the position of dead pixels by using the command (not(find))

find will give the list of pixels with non zero value in the image

here I am using (not(find)) to give the list of pixels with zero

value in the image

After that by using conv2() one of the surrounding pixel value wiil replace the dead pixel(pixel value with zero value)

Can u suggest any other method

Thanks in advance

[Jim]

Hi sean de,

Thank you for ur reply

How to take the kernel matrix K=[1 1 1;1 0 1;1 1 1] or

K=[1/8 1/8 1/8;1/8 0 1/8;1/8 1/8 1/8]

Thanks in advance

[Imran Riaz]

@Image Analyst @Jim Can u share the code ? I have same pbm and task.

[Image Analyst]

Not sure what you mean but see if my salt and pepper demo will do what you want.

Answer 2

[Image Analyst]

Here is some code I made up for another question. It's basically the same except that it used the median filter instead of the conv2 filter. It's trivial to modify:

clc;  % Clear command window.
clear;  % Delete all variables.
close all;  % Close all figure windows except those created by imtool.
imtool close all;  % Close all figure windows created by imtool.
workspace;  % Make sure the workspace panel is showing.
fontSize = 15;

% Read in a standard MATLAB color demo image.
folder = fullfile(matlabroot, '\toolbox\images\imdemos');
baseFileName = 'peppers.png';
fullFileName = fullfile(folder, baseFileName);
% Get the full filename, with path prepended.
fullFileName = fullfile(folder, baseFileName);
if ~exist(fullFileName, 'file')
  % Didn't find it there.  Check the search path for it.
  fullFileName = baseFileName; % No path this time.
  if ~exist(fullFileName, 'file')
    % Still didn't find it.  Alert user.
    errorMessage = sprintf('Error: %s does not exist.', fullFileName);
    uiwait(warndlg(errorMessage));
    return;
  end
end
rgbImage = imread(fullFileName);
% Get the dimensions of the image.  numberOfColorBands should be = 3.
[rows columns numberOfColorBands] = size(rgbImage);
% Display the original color image.
subplot(3, 4, 1);
imshow(rgbImage);
title('Original color Image', 'FontSize', fontSize);
% Enlarge figure to full screen.
set(gcf, 'Position', get(0,'Screensize')); 

% Extract the individual red, green, and blue color channels.
redChannel = rgbImage(:, :, 1);
greenChannel = rgbImage(:, :, 2);
blueChannel = rgbImage(:, :, 3);

% Display the individual red, green, and blue color channels.
subplot(3, 4, 2);
imshow(redChannel);
title('Red Channel', 'FontSize', fontSize);
subplot(3, 4, 3);
imshow(greenChannel);
title('Green Channel', 'FontSize', fontSize);
subplot(3, 4, 4);
imshow(blueChannel);
title('Blue Channel', 'FontSize', fontSize);

% Generate a noisy image.  This has salt and pepper noise independently on
% each color channel so the noise may be colored.
noisyRGB = imnoise(rgbImage,'salt & pepper', 0.05);
subplot(3, 4, 5);
imshow(noisyRGB);
title('Image with Salt and Pepper Noise', 'FontSize', fontSize);

% Extract the individual red, green, and blue color channels.
redChannel = noisyRGB(:, :, 1);
greenChannel = noisyRGB(:, :, 2);
blueChannel = noisyRGB(:, :, 3);

% Display the noisy channel images.
subplot(3, 4, 6);
imshow(redChannel);
title('Noisy Red Channel', 'FontSize', fontSize);
subplot(3, 4, 7);
imshow(greenChannel);
title('Noisy Green Channel', 'FontSize', fontSize);
subplot(3, 4, 8);
imshow(blueChannel);
title('Noisy Blue Channel', 'FontSize', fontSize);

% Median Filter the channels:
redMF = medfilt2(redChannel, [3 3]);
greenMF = medfilt2(greenChannel, [3 3]);
blueMF = medfilt2(blueChannel, [3 3]);

% Find the noise in the red.
noiseImage = (redChannel == 0 | redChannel == 255);
% Get rid of the noise in the red by replacing with median.
noiseFreeRed = redChannel;
noiseFreeRed(noiseImage) = redMF(noiseImage);

% Find the noise in the green.
noiseImage = (greenChannel == 0 | greenChannel == 255);
% Get rid of the noise in the green by replacing with median.
noiseFreeGreen = greenChannel;
noiseFreeGreen(noiseImage) = greenMF(noiseImage);

% Find the noise in the blue.
noiseImage = (blueChannel == 0 | blueChannel == 255);
% Get rid of the noise in the blue by replacing with median.
noiseFreeBlue = blueChannel;
noiseFreeBlue(noiseImage) = blueMF(noiseImage);

% Reconstruct the noise free RGB image
rgbFixed = cat(3, noiseFreeRed, noiseFreeGreen, noiseFreeBlue);
subplot(3, 4, 9);
imshow(rgbFixed);
title('Restored Image', 'FontSize', fontSize);