Subtract average image from a series of images
Show older comments
I'm having some difficulty subtracting an averaged image from a series of images from which the average was obtained. The images have been obtained from a MP4 video file.
This is what i have done:
Ven23 = VideoReader('MyVideo23.MP4'); %Read in the video file
Ven23frame1 = read(Ven23,1); % Read in the first frame of the video
sumImage23 = im2double(Ven23frame1); %Double the accuracy of the frame
for i=2:30
Ven23Frames= read(Ven23,i);
sumImage23 = sumImage23 + im2double(Ven23Frames);
end
avgVen23= imshow(sumImage23/30);
This all seems to work; I get an image which looks like the average. However, whenever I try to subtract this average image from any particular frame from the video I get this error:
Undefined operator '-' for
input arguments of type
'matlab.graphics.primitive.Image'.
Then I tried to save the average image and frames as bmp files, load them back in and subtract them. Doing this just returns a black screen. I also tried doubling the accuracy of the reloaded images to no avail.
I'm obviously doing something wrong and I would be grateful if anyone could point it out to me. Thanks.
3 Comments
Image Analyst
on 2 Mar 2016
You forgot to show the line of code that generated the error! What does this show:
whos avgVen23
And did you use im2double on the image you're going to subtract from also?
N B
on 2 Mar 2016
N B
on 3 Mar 2016
Accepted Answer
Guillaume
on 2 Mar 2016
0 votes
In view of the error, I would suspect that you're trying to subtract your mean image from a variable called exactly image. Unfortunately, at the point where you're doing the subtraction, the variable image does not exist, so instead matlab calls the function image which returns an object of type 'matlab.graphics.primitive.Image'. Subtraction is not defined for that (a graphics object), hence the error message. At the same time, you probably have an upside down picture of a boy's head popping up in a figure.
In conclusion:
- you have a bug in that you're using a variable you've never defined. Can't help you more, since you didn't show that part of the code.
- do not name your variable image
10 Comments
N B
on 2 Mar 2016
Guillaume
on 2 Mar 2016
Oh! I've just spotted the problem. Ignore the above. The problem is this:
avgVen23 = imshow(sumImage23/30);
Which should be:
avgVen23 = sumImage23/30;
imshow(avgVen23);
You're assigning the return value of imshow to your avgVen23. imshow returns a handle to a graphics object, not the original image.
On the other hand, the whole code can be simplified to:
Ven23 = VideoReader('MyVideo23.MP4'); %Read in the video file
avgVen23 = mean(read(Ven23, [1 30]), 4, 'native')
N B
on 3 Mar 2016
Guillaume
on 3 Mar 2016
If you're subtracting images, make sure you convert them to double first.
Image Analyst
on 3 Mar 2016
If you used double() and the differences are in the range of 0 - 5 or 10, then it's possible it will be too dark. I think if you used im2double() it would probably be okay, but it depends on how imshow() is scaling it. Either way, use [] in imshow() to "fix" the problem:
imshow(avgVen23, []);
N B
on 3 Mar 2016
Guillaume
on 3 Mar 2016
What is the maximum and minimum of the difference between the 2 images:
max(double(Ven23frame2(:)) - double(avgVen23(:))
min(double(Ven23frame2(:)) - double(avgVen23(:))
N B
on 3 Mar 2016
Guillaume
on 3 Mar 2016
I suspect that
class(Ven23frame2)
class(avgVen23)
will return 'uint8' and 'double' respectively, since the min/max in the first case (with double) is not an integer. That means you didn't use
avgVen23 = mean(read(Ven23, [1 30]), 4, 'native')
to generate the mean image but some other way, where at some point you converted the image with im2double.
You need to learn the difference between the classes of arrays, what arithmetic operations are possible with them and the effect it has on display.
An array can be of class uint8. That means pixel values are in the range 0-255. When you imshow such an array (without using [] in the imshow call), matlab assumes that range and display 0 as black and 255 as white. It's not a good idea to perform arithmetic on uint8 arrays as any result below 0 gets clamped to 0 and any result above 255 gets clamped to 255.
If the array is of class double. imshow (and most image processing functions) assume that the intensities are in the range 0-1. 0 and anything below 0 gets displayed as black, 1 and anything higher gets displayed as white.
If you convert a uint8 array to double using double(img), the intensities do not change (they're still in the range 0-255) and if you try to display that, because it is double intensities 1 to 254 will show as white. im2double on the other hand rescale the intensities in the range 0-1.
To avoid all this trouble with display always use imshow with []. It will use whatever range is in the array:
imshow(img , [])
When you do the min/max difference with im2double the difference is 0.1, which is 10% of the intensity range (0-1), so you should see some colour.
If you don't get anywhere, attach save the output of
read(Ven23, [1 30])
in a mat file and post that.
N B
on 3 Mar 2016
More Answers (0)
Categories
Find more on Image Arithmetic in Help Center and File Exchange
Products
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
