Index: trunk/src/transform_field/ima_erode_particles.m
===================================================================
--- trunk/src/transform_field/ima_erode_particles.m	(revision 906)
+++ 	(revision )
@@ -1,62 +1,0 @@
-% 'ima_remove_particles': removes particles from an image (keeping the local minimum)
-% requires the Matlab image processing toolbox
-%------------------------------------------------------------------------
-%%%%  Use the general syntax for transform fields with a single input %%%%
-% OUTPUT: 
-% DataOut:   output field structure 
-%
-%INPUT:
-% DataIn:  first input field structure
-
-%=======================================================================
-% Copyright 2008-2014, LEGI UMR 5519 / CNRS UJF G-INP, Grenoble, France
-%   http://www.legi.grenoble-inp.fr
-%   Joel.Sommeria - Joel.Sommeria (A) legi.cnrs.fr
-%
-%     This file is part of the toolbox UVMAT.
-%
-%     UVMAT is free software; you can redistribute it and/or modify
-%     it under the terms of the GNU General Public License as published
-%     by the Free Software Foundation; either version 2 of the license,
-%     or (at your option) any later version.
-%
-%     UVMAT is distributed in the hope that it will be useful,
-%     but WITHOUT ANY WARRANTY; without even the implied warranty of
-%     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
-%     GNU General Public License (see LICENSE.txt) for more details.
-%=======================================================================
-
-function DataOut=ima_remove_particles(DataIn)
-%------------------------------------------------------------------------
-DataOut=[];  %default  output field
-if strcmp(DataIn,'*')
-    return
-end
-
-%parameters
-radius=4;
-SE=strel('disk',radius);
-%--------------------------------------------------------- 
-DataOut=DataIn;%default
-
-if ~isfield(DataIn,'A')
-    DataOut.Txt='remove_particles only valid for input images';
-    return
-end
-[npy,npx]=size(DataIn.A);
-[X,Y]=meshgrid(1:npx,1:npy);
-%BACKGROUND LEVEL
-Atype=class(DataIn.A);
-%SE=ones(4);
-Aerode=imerode(DataIn.A,SE);
-Aflagmin=find(DataIn.A==Aerode);
-% Backg=zeros(size(A));
-%Aflagmin=imregionalmin(DataIn.A);%Amin=1 for local image minima
-Xmin=X(Aflagmin);
-Ymin=Y(Aflagmin);
-Amin=double(DataIn.A(Aflagmin));
-F = TriScatteredInterp([Xmin Ymin], Amin);
-DataOut.A=F(X,Y);
-
-% DataOut.A=feval(Atype,DataOut.A);
-