singleLED = load('hdrCalGamma_1x1');
LCDvalue=unique(singleLED.inputs(:,1));
%lowlevel=find(LEDvalue<=10);
%LEDvalue(lowlevel)=[];

figure;
hold on;
for ii=1:length(LCDvalue);
x = singleLED.inputs(find(singleLED.inputs(:,1)==LCDvalue(ii)),2);
y = singleLED.XYZ(find(singleLED.inputs(:,1)==LCDvalue(ii)),2);
plot(x,y/max(y),'ro'); 
plot(sort(x),sort(x/max(x)),'b-');
end
hold off;

xlabel('LED levels');
ylabel('luminance scaled to 1');
legend('measured','linear');
%legend('measured');
title('LED response with a single LED')
grid on;