Changeset 897 for trunk/src/transform_field

Timestamp:

May 25, 2015, 8:48:58 PM (10 years ago)

Author:

sommeria

Message:

sub_background corrected for volume and multilevel cases

Location:

trunk/src/transform_field

Files:

• ima_remove_background.m (modified) (1 diff)
• ima_remove_particles.m (modified) (3 diffs)

trunk/src/transform_field/ima_remove_background.m

@@ -54 +54 @@
 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);

trunk/src/transform_field/ima_remove_particles.m

@@ -32 +32 @@
 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 +50 @@
 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 +59 @@
 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);