Index: trunk/src/transform_field/ima_remove_background.m
===================================================================
--- trunk/src/transform_field/ima_remove_background.m	(revision 883)
+++ trunk/src/transform_field/ima_remove_background.m	(revision 897)
@@ -54,4 +54,5 @@
 B=imresize(Backg,size(A),'bilinear');% interpolate to the initial size image
 ImPart=(A-B);
-DataOut.A=ImPart.*(ImPart>threshold);
+DataOut.A=ImPart;
+%DataOut.A=ImPart.*(ImPart>threshold);
 DataOut.A=feval(Atype,DataOut.A);
Index: trunk/src/transform_field/ima_remove_particles.m
===================================================================
--- trunk/src/transform_field/ima_remove_particles.m	(revision 883)
+++ trunk/src/transform_field/ima_remove_particles.m	(revision 897)
@@ -32,10 +32,11 @@
 if strcmp(DataIn,'*')
     return
+    
 end
 
 %parameters
 threshold=200;
-nblock_x=30;%size of image subblocks for analysis
-nblock_y=30;
+nblock_x=10;%size of image subblocks for analysis
+nblock_y=10;
 %--------------------------------------------------------- 
 DataOut=DataIn;%default
@@ -49,5 +50,5 @@
 Atype=class(DataIn.A);
 A=double(DataIn.A);
-Backg=zeros(size(A));
+% Backg=zeros(size(A));
 Aflagmin=sparse(imregionalmin(A));%Amin=1 for local image minima
 Amin=A.*Aflagmin;%values of A at local minima
@@ -58,7 +59,8 @@
 Backg=Backg./Bmin; % find the average of minima in blocks
 B=imresize(Backg,size(A),'bilinear');% interpolate to the initial size image
-ImPart=(A-B);
-ImPart=ImPart.*(ImPart>threshold);
-DataOut.A=A-ImPart;%
+DataOut.A=B;
+% ImPart=(A-B);
+% ImPart=ImPart.*(ImPart>threshold);
+% DataOut.A=A-ImPart;%
 DataOut.A=feval(Atype,DataOut.A);