41 addUsageLine(
"Remove the flickering in a tilt series. This is a phenomenon rather ");
42 addUsageLine(
"common in X-ray microscopy. For doing so, a low-pass filter of a given ");
43 addUsageLine(
"frequency (normalized to 0.5) is applied across the time series, i.e., a ");
44 addUsageLine(
"line is formed with all the lines at the same position along the tilt series. ");
45 addUsageLine(
"This time series is the filtered and put back into the tilt series. ");
46 addParamsLine(
" -i <file> : Input images (selfile or stack)");
57 size_t Zdim,
dummy, Ydim, Xdim, Ndim;
61 V().initZeros(Zdim,Ydim,Xdim);
65 for (
size_t objId :
SF.
ids())
79 std::cout <<
"Removing fluctuations from a tilt series\n";
80 std::cout <<
"Input series: " <<
fnIn << std::endl
81 <<
"Output root: " <<
fnOut << std::endl
82 <<
"Max freq (LPF): " <<
maxFreq << std::endl
99 std::complex<double>
zero=0;
115 for (
size_t k=maxPixel;
k<
XSIZE(lineFourier);
k++)
void init_progress_bar(long total)
double getDoubleParam(const char *param, int arg=0)
double maxFreq
Cutoff frequency of the lowpass filter (<0.5)
void write(const FileName &name="", size_t select_img=ALL_IMAGES, bool isStack=false, int mode=WRITE_OVERWRITE, CastWriteMode castMode=CW_CAST, int _swapWrite=0)
void readParams()
Read parameters from command line.
FileName fnOut
Rootname for the output.
ql0001_ & k(htemp+1),(cvec+1),(atemp+1),(bj+1),(bl+1),(bu+1),(x+1),(clamda+1), &iout, infoqp, &zero,(w+1), &lenw,(iw+1), &leniw, &glob_grd.epsmac
#define DIRECT_A1D_ELEM(v, i)
const char * getParam(const char *param, int arg=0)
void show() const
Show parameters.
void produceSideInfo()
Produce side info.
void progress_bar(long rlen)
int verbose
Verbosity level.
#define DIRECT_A3D_ELEM(v, k, i, j)
void setDimensions(int Xdim, int Ydim, int Zdim, size_t Ndim)
int read(const FileName &name, DataMode datamode=DATA, size_t select_img=ALL_IMAGES, bool mapData=false, int mode=WRITE_READONLY)
void defineParams()
Define params.
void addUsageLine(const char *line, bool verbatim=false)
void initZeros(const MultidimArray< T1 > &op)
void addParamsLine(const String &line)
ql0001_ & zero(ctemp+1),(cvec+1),(a+1),(b+1),(bl+1),(bu+1),(x+1),(w+1), &iout, ifail, &zero,(w+3), &lwar2,(iw+1), &leniw, &glob_grd.epsmac
FileName fnIn
Input images.