When going through the filter code converting to stacks I noticed that the code may be adjusting for downsampling in two ways and Im not sure it should do both.
One scaling is via the pixel_size, which scales bandwidth applied to radial grid inside phase flipping. The radial grid is derived from the omega passed from evaluate.
Another scaling occurs via filter.scale which scales omega directly (before passing to filter.evaluate).
We generally get best results with phase flipping at full image size. I think it is worth reviewing the code for correctness in case this is why the phase flipping after downsample under performs expectations.
When going through the filter code converting to stacks I noticed that the code may be adjusting for downsampling in two ways and Im not sure it should do both.
One scaling is via the
pixel_size, which scales bandwidth applied to radial grid inside phase flipping. The radial grid is derived from theomegapassed fromevaluate.Another scaling occurs via
filter.scalewhich scalesomegadirectly (before passing tofilter.evaluate).We generally get best results with phase flipping at full image size. I think it is worth reviewing the code for correctness in case this is why the phase flipping after downsample under performs expectations.