Чтение и отображение изображения с сильным желтым отливом.

A = imread('hallway.jpg');
imshow(A)
title('Original Image')

Выберите пиксель на изображении, который должен выглядеть белым или серым, например точку на столбе. Не выбирайте насыщенный пиксель, например точку на потолке. Выбранная точка отображается зеленым цветом.

x = 2800;
y = 1000;
gray_val = impixel(A,x,y);
drawpoint('Position',[x y],'Color','g');

Используйте выбранный цвет в качестве ссылки для освещения сцены и исправьте баланс белого изображения.

B = chromadapt(A,gray_val);

Отображение скорректированного изображения. Столбы теперь выглядят белыми, как и ожидалось, а остальная часть изображения не имеет желтого оттенка.

imshow(B)
title('White-Balanced Image')

Откройте файл изображения, содержащий минимально обработанные линейные интенсивности RGB.

A = imread('foosballraw.tiff');

Данные изображения являются необработанными данными датчика после коррекции уровня черного и масштабирования до 16 бит на пиксель. Интерполяция интенсивностей для восстановления цвета. Массив цветовых фильтров имеет формат RGGB.

A = demosaic(A,'rggb');

Отображение изображения. Поскольку изображение находится в линейном цветовом пространстве RGB, примените гамма-коррекцию, чтобы изображение отображалось на экране правильно.

A_sRGB = lin2rgb(A);
imshow(A_sRGB)
title('Original Image')

Изображение имеет диаграмму X-Rite ® ColorChecker ® в сцене. Чтобы получить цвет окружающего света, отобразите значения RGB пикселя в одном из нейтральных фрагментов диаграммы. Интенсивность красного канала ниже, чем интенсивность двух других каналов, что указывает на то, что свет голубовато-зелёный.

x = 1510;
y = 1250;
light_color = [A(y,x,1) A(y,x,2) A(y,x,3)]

light_color = 1x3 uint16 row vector

    7361   14968   10258

Сбалансировать цветовые каналы изображения. Используйте 'ColorSpace' чтобы указать, что изображение и светильник выражены в линейном RGB.

B = chromadapt(A,light_color,'ColorSpace','linear-rgb');

Отображение цветосбалансированного изображения с гамма-коррекцией.

B_sRGB = lin2rgb(B);
imshow(B_sRGB)
title('Color-Balanced Image')

Убедитесь, что серый фрагмент сбалансирован по цвету. Три цветовых канала в цветовом сером пятне имеют одинаковую интенсивность, как и ожидалось.

patch_color = [B(y,x,1) B(y,x,2) B(y,x,3)]

patch_color = 1x3 uint16 row vector

   13010   13010   13010