function T = WlsToT(wls,sd)
% T = WlsToT(wls,[sd])
%
% Produce the identity color matching matrix associated
% with the passed wavelength sampling.
%
% Arguments may be either [start delta n] description or
% a list of wavelengths.
%
% The second input sd is the standard deviation of a Gaussian
% blur function which may be incorporated into the T matrix
% if desired.
%
% This second feature should be used with some
% caution, as I haven't thought it completely through.  The
% current implementation normalizes each row of T so that the
% entries sum to 1.  This seems as if it is what we want to
% get answers consistent with what we get when we don't
% incorporate wavelength blurring.
%
% 7/11/03  dhb  A few comments added.

% Force into wavelength representation.
wls = MakeItWls(wls);

[m, n] = size(wls);
if nargin == 1
    T = eye(m,m);
else
    T = zeros(m,m);
    for i = 1:m
        T(i,:) = NormalPDF(wls,wls(i),sd^2)';
        T(i,:) = T(i,:) / sum(T(i,:)');
    end
end