%CMFXYZ	Color matching function%%	XYZ = CMFXYZ(lambda)%	XYZ = CMFXYZ(lambda, spect)%%	Return the CIE color matching function, where each row of XYZ%	corresponds to the elements of lambda.%% SEE ALSO:  CCXYZ%%%	  By Peter Corke, Copyright (c) Peter Corke, 1999  Machine Vision Toolbox for Matlab%			  Peter Corke, 3/95, pic@brb.dmt.csiro.au% $Header: /home/autom/pic/cvsroot/image-toolbox/cmfxyz.m,v 1.2 2005/11/24 11:14:35 pic Exp $% $Log: cmfxyz.m,v $% Revision 1.2  2005/11/24 11:14:35  pic% Change to copyright notices.%% Revision 1.1.1.1  2002/05/26 10:50:20  pic% initial import%function xyz = cmfxyz(lambda, spect)	ciedat = [	0	0	0	0	380	0.0014	0.0000	0.0065	385	0.0022	0.0001	0.0105	390	0.0042	0.0001	0.0201	395	0.0076	0.0002	0.0362	400	0.0143	0.0004	0.0679	405	0.0232	0.0006	0.1102	410	0.0435	0.0012	0.2074	415	0.0776	0.0022	0.3713	420	0.1344	0.0040	0.6456	425	0.2148	0.0073	1.0391	430	0.2839	0.0116	1.3856	435	0.3285	0.0168	1.6230	440	0.3483	0.0230	1.7471	445	0.3481	0.0298	1.7826	450	0.3362	0.0380	1.7721	455	0.3187	0.0480	1.7441	460	0.2908	0.0600	1.6692	465	0.2511	0.0739	1.5281	470	0.1954	0.0910	1.2876	475	0.1421	0.1126	1.0419	480	0.0956	0.1390	0.8130	485	0.0580	0.1693	0.6162	490	0.0320	0.2080	0.4652	495	0.0147	0.2586	0.3533	500	0.0049	0.3230	0.2720	505	0.0024	0.4073	0.2123	510	0.0093	0.5030	0.1582	515	0.0291	0.6082	0.1117	520	0.0633	0.7100	0.0782	525	0.1096	0.7932	0.0573	530	0.1655	0.8620	0.0422	535	0.2257	0.9149	0.0298	540	0.2904	0.9540	0.0203	545	0.3597	0.9803	0.0134	550	0.4334	0.9950	0.0087	555	0.5121	1.0000	0.0057	560	0.5945	0.9950	0.0039	565	0.6784	0.9786	0.0027	570	0.7621	0.9520	0.0021	575	0.8425	0.9154	0.0018	580	0.9163	0.8700	0.0017	585	0.9786	0.8163	0.0014	590	1.0263	0.7570	0.0011	595	1.0567	0.6949	0.0010	600	1.0622	0.6310	0.0008	605	1.0456	0.5668	0.0006	610	1.0026	0.5030	0.0003	615	0.9384	0.4412	0.0002	620	0.8544	0.3810	0.0002	625	0.7514	0.3210	0.0001	630	0.6424	0.2650	0.0000	635	0.5419	0.2170	0.0000	640	0.4479	0.1750	0.0000	645	0.3608	0.1382	0.0000	650	0.2835	0.1070	0.0000	655	0.2187	0.0816	0.0000	660	0.1649	0.0610	0.0000	665	0.1212	0.0466	0.0000	670	0.0874	0.0320	0.0000	675	0.0636	0.0232	0.0000	680	0.0468	0.0170	0.0000	685	0.0329	0.0119	0.0000	690	0.0227	0.0082	0.0000	695	0.0158	0.0057	0.0000	700	0.0114	0.0041	0.0000	705	0.0081	0.0029	0.0000	710	0.0058	0.0021	0.0000	715	0.0041	0.0015	0.0000	720	0.0029	0.0010	0.0000	725	0.0020	0.0007	0.0000	730	0.0014	0.0005	0.0000	735	0.0010	0.0004	0.0000	740	0.0007	0.0002	0.0000	745	0.0005	0.0002	0.0000	750	0.0003	0.0001	0.0000	755	0.0002	0.0001	0.0000	760	0.0002	0.0001	0.0000	765	0.0001	0.0000	0.0000	770	0.0001	0.0000	0.0000	775	0.0001	0.0000	0.0000	780	0.0000	0.0000	0.0000	1e9	0	0	0];	xyz = interp1(ciedat(:,1)*1e-9, ciedat(:,2:4), lambda);	if nargin == 2,		xyz = spect(:)' * xyz;	end