%?Rpeak?=?...
%?Gpeak?=?...
%?Bpeak?=?...
%?Rthresh?=?...
%?Gthresh?=?...
%?Bthresh?=?...
I5?=?zeros(M,?N,?'uint8');
I1R?=?(I1(:,:,1)-Rpeak)<Rthresh;
I1G?=?(I1(:,:,2)-Gpeak)<Gthresh;
I1B?=?(I1(:,:,3)-Bpeak)<Bthresh;
I5(I1R&I1G&I1B)?=?255;