compare macbeth data
Contents
clear
load sess2
load macbeth
Macbeth data
Y = macbeth_cie(19:24,2);
g = macbeth_rgb(19:24,2);
xlog = log10(Y);
ylog = log10(g);
c = polyfit(xlog,ylog,1);
b = c(1);
A = 10^c(2);
fprintf('power law A = %g b = %g\n',A,b);
xr = linspace(0,100,200);
plot(Y,g,'ko');
hold on
plot(xr,A*xr.^b,'k-');
hold off
xlabel('Y');
ylabel('rgb');
grid;
power law A = 30.9576 b = 0.458003
my data
Y = macbeth_cie(19:24,2);
g = rgb2(19:24,2)*255;
xlog = log10(Y);
ylog = log10(g);
c = polyfit(xlog,ylog,1);
b2 = c(1);
A2 = 10^c(2);
fprintf('power law A = %g b = %g\n',A2,b2);
xr = linspace(0,100,200);
plot(Y,g,'ko');
hold on
plot(xr,A2*xr.^b2,'k-');
hold off
xlabel('Y');
ylabel('rgb');
grid;
power law A = 18.1175 b = 0.612007
g1 = macbeth_rgb(19:24,2);
plot(Y,g,'ko',Y,g1,'bx');
legend('my data','Macbeth standard');
hold on
plot(xr,A*xr.^b,'k-');
hold off
xlabel('Y');
ylabel('rgb');
grid;