Color normalization algorithm under various lighting conditions

2020 4 24
2 个回答
回答已采纳
更新时间：2020 4 28
5 次查看（30 天）

0 个投票

I am trying to take pictures using cellphones under various lighting conditions which would cause color distortions. Hence some kind of white balance algorithm is needed to correct the lighting of those photos so we get the standardized color. This can be commonly done if you have a "gray card" which many photographers use to tell their cameras what is the true white color under any light. With the "gray card", a person can take a picture with a standard color item by the side. With the item's color as a reference, we hope we can standardize the picture's color. We are doing this for an excellent color normalization.
Does anybody know any readily available algorithms (or something we can modify) for the above purpose? I really appreciate it if you can offer any suggestions.

4 个评论
显示 2更早的评论 隐藏 2更早的评论

darova
darova 2020-4-24
Can you take a picture with some identificator? Maybe slice of white paper?
LANG Wu
LANG Wu 2020-4-24
Yes, we can. After we take the picture with identificator, what algorithms can I use to normalize the color?
darova
darova 2020-4-24
Is it possible to make that identificator to be in the same place?
LANG Wu
LANG Wu 2020-4-24
That's hard to say. But maybe we can make it in the top right or left. Or we can take two pictures, one for the object, the other for the identificator.

请先登录，再进行评论。

请先登录，再回答此问题。

 采纳的回答

Image Analyst
Image Analyst 2020-4-24

0 个投票

There is a whole family of functions to do this. In the Image Processing Toolbox, check out illumgray(), illumpca(), illumwhite(), and chromadapt().

6 个评论
显示 4更早的评论 隐藏 4更早的评论

LANG Wu
LANG Wu 2020-4-24
Okay, thank you. I will check them out.
LANG Wu
LANG Wu 2020-4-25
It seems they are all auto white balance algorithms, which are not that helpful.
Image Analyst
Image Analyst 2020-4-25
I don't know why not. If you want to do it yourself, this is what I'd do
  1. Convert the image into lab colorspace with rgb2lab()
  2. Take the histogram of the L channel and threshold the channel to get your white/gray card. This will get you a mask
  3. Extract color channels: redChannel = rgbImage(:,:,1);, etc.
  4. Get the mean gray levels for every color channel using meanR = mean(redChannel(mask)), etc.
  5. Compute the max intensity maxGL = max([meanR, meanG, meanB])
  6. Compute the deltas : deltaR = uint8(maxGL - meanR), etc.
  7. Add the deltas to each color channel: newR = redChannel + deltaR, newG = greenChannel + deltaG, etc.
  8. Recombine into new RGB image newRGB = cat(3, newR, newG, newB)
This will color correct the image such that the brightest color channel of the white card does not change.
LANG Wu
LANG Wu 2020-4-28
Thanks for your answer. But can you explain more on how to implement step 2?
Image Analyst
Image Analyst 2020-4-28
Actually you can do steps 1 and 2 using rgb2gray()
grayImage = rgb2gray(rgbImage);
mask = grayImage > 180; % or whatever gray level your white card is.
% Extract the individual red, green, and blue color channels.
redChannel = rgbImage(:, :, 1);
greenChannel = rgbImage(:, :, 2);
blueChannel = rgbImage(:, :, 3);
meanR = mean(redChannel(mask))
meanG = mean(greenChannel(mask))
meanB = mean(blueChannel(mask))
maxR = max(redChannel(mask))
maxG = max(greenChannel(mask))
maxB = max(blueChannel(mask))
deltaR = uint8(maxR - meanR)
deltaG = uint8(maxG - meanG)
deltaB = uint8(maxB - meanB)
newR = redChannel + deltaR;
newG = greenChannel + deltaG;
newB = blueChannel + deltaB;
newRGB = cat(3, newR, newG, newB);
imshow(newRGB);
See attached color standardization and calibration tutorial from my course I teach.
LANG Wu
LANG Wu 2020-4-28
Thanks again. Actually what I want to do is select the white card area, get its illumination and tell the algorithm this is 'white', and then normalize the image color considering the white card's illumination as standard. I don't see how you define the 'white' color in your code. Why do you choose 'grayscale>180" and what does 'delta' mean?

请先登录，再进行评论。

更多回答（1 个）

darova
darova 2020-4-25

1 个投票

Try this trick
I = imread('peppers.png');
imshow(I)
p = round(ginput(1));
for i = 1:3
I(:,:,i) = I(:,:,i) + (255 - I(p(2),p(1),i));
end
imshow(I)

类别

在 帮助中心 和 File Exchange 中查找有关 Creating, Deleting, and Querying Graphics Objects 的更多信息

产品

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by