Abra una imagen de prueba. La imagen es el dato sin procesar capturado con una cámara digital Canon EOS 30D después de corregir el nivel de negro y escalar las intensidades a 16 bits por píxel. No se ha aplicado demosaicing, balance de blancos, mejora de color, filtrado de ruido o corrección gamma.

A = imread('foosballraw.tiff');

Interpolar Utilizodo laFunciónpara Obener Una Imagen De Color。demosaicEl patrón de matriz de filtro de color es RGGB.

a_demosaiced = demosaic（a，'rggb');

La imagen contiene un gráfico ColorChecker. Especifique el iluminador de la verdad del suelo, que se calculó de antemano utilizando los parches neutros del gráfico.

illuminant_groundtruth = [0.0717 0.1472 0.0975];

Para evitar sesgar la estimación del iluminante, excluya el gráfico ColorChecker creando una máscara.

mask = true(size(A_demosaiced,1), size(A_demosaiced,2)); mask(920:1330,1360:1900) = false;

Ejecute tres algoritmos de estimación de iluminante diferentes: , , y .illumwhiteIllumgray.illumpca

Illuminant_Whitepatch = IllumWhite（A_Demosaiced，'Mask',mask); illuminant_grayworld = illumgray(A_demosaiced,'Mask',mask); illuminant_pca = illumpca(A_demosaiced,'Mask',mask);

Compare cada estimación con la verdad del suelo calculando el ángulo entre cada iluminador estimado y la verdad del suelo utilizando la función.罗缎Cuanto Menor SeaElÁngulo，MejorSeráLaiveración。La Magnitud de laErgeraciónnodoctaporquesóloladireccióndeliLuinante se Utiliza para equilibrar Una Imagen ConAdaptaciónCromática。

angle_whitepatch = colorangle(illuminant_whitepatch, illuminant_groundtruth)

angle_whitepatch = 5.0921

angle_grayworld = colorangle(illuminant_grayworld, illuminant_groundtruth)

angle_grayworld = 5.1036

Angle_pca = icorarangle（Illuminant_pca，Illiminant_groundtruth）

angle_pca = 5.0134

El valor de es más pequeño, lo que indica que el algoritmo de estimación del iluminante PCA está más cerca de la iluminación de la verdad del suelo para esta imagen.angle_pca