allow specifying a region for ringRemove

functions/tomography/ringremove.m

@@ -30,6 +30,9 @@ function sinoCorrected = ringremove(sino,settings)
 % additive_Strength : Default = 1
 %     Strength of the subtraction of the identified stripe artifacts from the
 %     original sinogram when using method 'additive'.
+% additive_Region : Default = []
+%     Region (of the radius) of the sinogram to be corrected when using method 
+%     'additive'. [] corresponds to 1:size(sino,1).
 % wavelet_Wname : Default = 'sym8'
 %     Wavelet base to be used for the decomposition of the sinogram when using method
 %     'wavelet'.
@@ -106,7 +109,8 @@ defaults.method = 'additive';
 defaults.additive_SmoothMethod = 'mean'; 
 defaults.additive_WindowSize = 10;      
 defaults.additive_GaussSigma = 3;        
-defaults.additive_Strength = 1;         
+defaults.additive_Strength = 1;
+defaults.additive_Region = [];
 
 defaults.wavelet_Wname = 'sym8';
 defaults.wavelet_DecNum = 3;    
@@ -136,8 +140,16 @@ switch settings.method
             otherwise
                 error('Choose between "mean", "median" and "gauss" for setting "additive_SmoothMethod"');
         end
+
+		if ~isempty(settings.additive_Region)
+			region = settings.additive_Region;
+		else
+			region = 1:size(sino,1);
+		end
+
         % identification of horizontal stripes
-        correctionProfile = sinoMean - sinoMeanSmoothed;
+        correctionProfile = zeros(size(sinoMean));
+        correctionProfile(region) = sinoMean(region) - sinoMeanSmoothed(region);
         % sinogram correction
         sinoCorrected = sino - settings.additive_Strength * correctionProfile;