%IM_SELECT_BLOB Fixed mapping selecting largest blob in binary images (DIP_Image)
%
% B = IM_SELECT_BLOB(A)
% B = A*IM_SELECT_BLOB
%
% Just the largest object in the image is returned.
%
% SEE ALSO (<a href="http://37steps.com/prtools">PRTools Guide</a>)
% DATASETS, DATAFILES, DIP_IMAGE
% Copyright: R.P.W. Duin, r.p.w.duin@37steps.com
% Faculty EWI, Delft University of Technology
% P.O. Box 5031, 2600 GA Delft, The Netherlands
function b = im_select_blob(a)
if nargin < 1 | isempty(a)
b = prmapping(mfilename,'fixed');
b = setname(b,'Select largest blob');
elseif isa(a,'prdataset') % allows datafiles too
isobjim(a);
b = filtim(a,mfilename);
b = setfeatsize(b,getfeatsize(a));
elseif isa(a,'double') | isa(a,'dip_image') % here we have a single image
if ~isa(a,'dip_image')
a = dip_image(a,'bin');
end;
labim = label(a);
sz = measure(labim,labim,{'size'});
[cc,ind] = max(double(sz));
I = double(measure(labim,labim,{'mean'}));
b = a.*(labim==round(I(ind)));
% c = measure(labim,labim,{'size','mean'});
% c = double(c);
% [cc,ind] = max(c(:,1));
% b = a.*(labim==round(c(ind,2)));
end
return