summaryrefslogtreecommitdiff
diff options
context:
space:
mode:
-rw-r--r--striperemoval.m24
1 files changed, 14 insertions, 10 deletions
diff --git a/striperemoval.m b/striperemoval.m
index 7563bf0..8c2f623 100644
--- a/striperemoval.m
+++ b/striperemoval.m
@@ -1,13 +1,17 @@
function [img_cleaned]=striperemoval(img, varargin)
-% Cleans an image from the stripes produced by interference fringes.
-% img - image with stripes to be removed
-% varargin - not used, it is here for compatibility
+%STRIPEREMOVAL removes interference fringes/stripes from an image.
+% [img_cleaned]=STRIPEREMOVAL(img, varargin)
+% img - image with stripes to be removed
+% varargin - not used, it is here for compatibility
% By Eugeniy E. Mikhailov 2014/11/25
%
-% be aware of "slow" spatial stripes their Fourier components will overlap with Gaussian.
-% Such parallel fringes appears like strong isolated peaks in the Fourier transformed image.
-% What is even nice they appear symmetrical with respect to origin but not symmetrical to
-% fold
+% Be aware of "slow" spatial stripes
+% i.e. stripes with characteristic width comparable to the Gaussian beam.
+% Their Fourier components will overlap with Gaussian.
+
+% Fringes appear like strong isolated peaks in the Fourier transformed image.
+% What is even nicer they appear symmetrically with respect to origin but not symmetrical to
+% the fold with respect to X or Y axes. See figure below
%
% ^
% |
@@ -15,13 +19,13 @@ function [img_cleaned]=striperemoval(img, varargin)
% |
% oOo
% -------------O0O-------------->
-% oOo <- gaussian beam fft
+% oOo <- Gaussian beam fft
% |
-% | * <- stripe fft peak
+% | * <- stripes fft peak
% |
%
%
-% so it easy to clean them, especially if their spatial frequency is high.
+% So it easy to clean them, especially if their spatial frequency is high.
%
img_fourier=fftshift(fft2(img)); % move to Fourier space