|
|
Appendix
% fifth, calculate the FSI (full spectrum index)
% interpolate SPD from 380nm to 730nm with 1nm internal
t=380:1:730;
spd=interp1(wavelength_spd,spd,t','spline');
spd(isnan(spd)) = 0.0;
spd=spd/max(spd);
% Equal energy accumulative spd
EEspd=(1/351:1/351:1);
% circular calculation of FSI
for j=1:351
for i=1:351
cum(i,:)=sum(spd(1:i,:),1)./sum(spd,1);
leastwo(i,:)=(cum(i,:)-EEspd(i)).^2;
end
sumleastwo(j,:)=sum(leastwo,1);
spd=circshift(spd,1);
end
meanleastwo=mean(sumleastwo,1);
FSI=100-4*meanleastwo;
|