From 6309c880257d13b01f91e9c69097b683a2f87e88 Mon Sep 17 00:00:00 2001 From: jklebes Date: Thu, 19 Jan 2023 15:31:46 +0000 Subject: [PATCH] bug in hard stripe filter --- imbandpass.m | 8 +++++--- 1 file changed, 5 insertions(+), 3 deletions(-) diff --git a/imbandpass.m b/imbandpass.m index e672417..71830b4 100644 --- a/imbandpass.m +++ b/imbandpass.m @@ -85,22 +85,24 @@ %potentially add stripe to Fourier space mask switch p.Results.stripeFilter case "hard" + %the stripe Filter has this Fourier space width + stripeWidth=(100-p.Results.stripeTolerance)/100; %add a stripe of 0s to the mask switch p.Results.stripes case 'Horizontal' %Fourier space stripe goes other way than in real space! - for col = floor(center_coord_x-width/2):ceil(center_coord_x+p.Results.stripeTolerance/2) + for col = floor(center_coord_x-stripeWidth/2):ceil(center_coord_x+p.Results.stripeTolerance/2) mask(:,col)=0; end case 'Vertical' - for row = floor(center_coord_y-width/2):ceil(center_coord_y+p.Results.stripeToleranceh/2) + for row = floor(center_coord_y-stripeWidth/2):ceil(center_coord_y+p.Results.stripeTolerance/2) mask(row,:)=0; end end case "gaussian" - %add a gaussian stripe filter to the mask %the stripe Filter has this Fourier space width stripeWidth=(100-p.Results.stripeTolerance)/100; + %add a gaussian stripe filter to the mask switch p.Results.stripes case 'Horizontal' %Fourier space stripe goes other way than in real space!