Xmipp  v3.23.11-Nereus
Public Member Functions | Static Public Member Functions | Public Attributes | Static Public Attributes | List of all members
ProgCTFBasicParams Class Reference

#include <ctf_estimate_from_psd_base.h>

Inheritance diagram for ProgCTFBasicParams:
Inheritance graph
[legend]
Collaboration diagram for ProgCTFBasicParams:
Collaboration graph
[legend]

Public Member Functions

 ProgCTFBasicParams ()
 
void readParams ()
 Read parameters. More...
 
void readBasicParams (XmippProgram *program)
 Read parameters. More...
 
void show ()
 Show parameters. More...
 
void defineParams ()
 Define Parameters. More...
 
void produceSideInfo ()
 Produce side information. More...
 
void generate_model_halfplane (int Ydim, int Xdim, MultidimArray< double > &model)
 
void generate_model_quadrant (int Ydim, int Xdim, MultidimArray< double > &model)
 
- Public Member Functions inherited from XmippProgram
const char * getParam (const char *param, int arg=0)
 
const char * getParam (const char *param, const char *subparam, int arg=0)
 
int getIntParam (const char *param, int arg=0)
 
int getIntParam (const char *param, const char *subparam, int arg=0)
 
double getDoubleParam (const char *param, int arg=0)
 
double getDoubleParam (const char *param, const char *subparam, int arg=0)
 
float getFloatParam (const char *param, int arg=0)
 
float getFloatParam (const char *param, const char *subparam, int arg=0)
 
void getListParam (const char *param, StringVector &list)
 
int getCountParam (const char *param)
 
bool checkParam (const char *param)
 
bool existsParam (const char *param)
 
void addParamsLine (const String &line)
 
void addParamsLine (const char *line)
 
ParamDefgetParamDef (const char *param) const
 
virtual void run ()
 
virtual void quit (int exit_code=0) const
 
virtual int tryRun ()
 
void initProgress (size_t total, size_t stepBin=60)
 
void setProgress (size_t value=0)
 
void endProgress ()
 
void processDefaultComment (const char *param, const char *left)
 
void setDefaultComment (const char *param, const char *comment)
 
virtual void initComments ()
 
void setProgramName (const char *name)
 
void addUsageLine (const char *line, bool verbatim=false)
 
void clearUsage ()
 
void addExampleLine (const char *example, bool verbatim=true)
 
void addSeeAlsoLine (const char *seeAlso)
 
void addKeywords (const char *keywords)
 
const char * name () const
 
virtual void usage (int verb=0) const
 
virtual void usage (const String &param, int verb=2)
 
int version () const
 
virtual void show () const
 
virtual void read (int argc, const char **argv, bool reportErrors=true)
 
virtual void read (int argc, char **argv, bool reportErrors=true)
 
void read (const String &argumentsLine)
 
 XmippProgram ()
 
 XmippProgram (int argc, const char **argv)
 
virtual ~XmippProgram ()
 

Static Public Member Functions

static void defineBasicParams (XmippProgram *program)
 Define basic parameters. More...
 

Public Attributes

FileName fn_psd
 CTF filename. More...
 
double downsampleFactor
 Downsample performed. More...
 
Image< double > ctftomodel
 CTF amplitude to model. More...
 
Image< double > enhanced_ctftomodel
 CTF amplitude to model. More...
 
Image< double > enhanced_ctftomodel_fullsize
 CTF amplitude to model. More...
 
bool show_optimization
 Show convergence values. More...
 
bool selfEstimation
 
int ctfmodelSize
 X dimension of particle projections (-1=the same as the psd) More...
 
bool bootstrap
 Bootstrap estimation. More...
 
bool refineAmplitudeContrast
 Refine amplitude contrast. More...
 
bool fastDefocusEstimate
 Fast defocus estimate. More...
 
bool noDefocusEstimate
 No defocus estimate. More...
 
double lambdaPhase
 Regularization factor for the phase direction and unwrapping estimates (used in Zernike estimate) More...
 
int sizeWindowPhase
 Size of the average window used during phase direction and unwrapping estimates (used in Zernike estimate) More...
 
double min_freq
 Minimum frequency to adjust. More...
 
double max_freq
 Maximum frequency to adjust. More...
 
double Tm
 Sampling rate. More...
 
double defocus_range
 Defocus range. More...
 
double f1
 Enhancement filter low freq. More...
 
double f2
 Enhancement filter high freq. More...
 
double enhanced_weight
 Weight of the enhanced image. More...
 
Matrix1D< double > adjust
 Set of parameters for the complete adjustment of the CTF. More...
 
int modelSimplification
 Model simplification. More...
 
MultidimArray< double > x_contfreq
 Frequencies in axes. More...
 
MultidimArray< double > y_contfreq
 
MultidimArray< double > w_contfreq
 
MultidimArray< double > x_digfreq
 
MultidimArray< double > y_digfreq
 
MultidimArray< double > w_digfreq
 
MultidimArray< double > psd_exp_radial
 PSD data. More...
 
MultidimArray< double > psd_exp_enhanced_radial
 
MultidimArray< double > psd_exp_enhanced_radial_2
 
MultidimArray< double > psd_exp_enhanced_radial_derivative
 
MultidimArray< double > psd_theo_radial_derivative
 
MultidimArray< double > psd_exp_radial_derivative
 
MultidimArray< double > psd_theo_radial
 
MultidimArray< double > w_digfreq_r_iN
 
MultidimArray< double > w_digfreq_r
 
MultidimArray< std::complex< double > > psd_fft
 
std::vector< double > amplitud
 
MultidimArray< double > mask
 Masks. More...
 
MultidimArray< double > mask_between_zeroes
 
MultidimArray< double > w_count
 
double value_th
 
double min_freq_psd
 
double max_freq_psd
 
double max_gauss_freq
 
int show_inf
 
int action
 
double corr13
 
bool penalize
 
int evaluation_reduction
 
double heavy_penalization
 
double current_penalty
 
MultidimArray< double > * f
 
Matrix1D< double > * adjust_params
 
- Public Attributes inherited from XmippProgram
bool doRun
 
bool runWithoutArgs
 
int verbose
 Verbosity level. More...
 
int debug
 

Static Public Attributes

static constexpr double penalty = 32.0
 

Additional Inherited Members

- Protected Member Functions inherited from XmippProgram
void defineCommons ()
 
- Protected Attributes inherited from XmippProgram
int errorCode
 
ProgramDefprogDef
 Program definition and arguments parser. More...
 
std::map< String, CommentListdefaultComments
 
int argc
 Original command line arguments. More...
 
const char ** argv
 

Detailed Description

Definition at line 17 of file ctf_estimate_from_psd_base.h.

Constructor & Destructor Documentation

◆ ProgCTFBasicParams()

ProgCTFBasicParams::ProgCTFBasicParams ( )
inline

Definition at line 128 of file ctf_estimate_from_psd_base.h.

129  {
130  }

Member Function Documentation

◆ defineBasicParams()

void ProgCTFBasicParams::defineBasicParams ( XmippProgram program)
static

Define basic parameters.

Definition at line 99 of file ctf_estimate_from_psd_base.cpp.

100 {
101  program->addSeeAlsoLine("ctf_enhance_psd");
102  program->addParamsLine("== Downsampling");
103  program->addParamsLine(
104  " [--downSamplingPerformed <F=1>] : Downsampling performed to produce this PSD");
105  program->addParamsLine(
106  " : Note that the output CTF model is referred to the original sampling rate");
107  program->addParamsLine(
108  " : not the one of the downsampled image.");
109  program->addParamsLine("== CTF fit: Optimization constraints");
110  program->addParamsLine(
111  " [--min_freq <fmin=0.03>] : Minimum digital frequency (<0.5) to use in adjust. Its value");
112  program->addParamsLine(
113  " : should be a little lower than the dig. freq. of the first ");
114  program->addParamsLine(" : CTF zero.");
115  program->addParamsLine(
116  " [--max_freq <fmax=0.35>] : Maximum digital frequency (<0.5) to use in adjust.");
117  program->addParamsLine(
118  " : It should be higher than the last zero of the CTF.");
119  program->addParamsLine(
120  " [--fastDefocus <lambda=2> <size=10>] : Estimate first defocus with Zernike polynomials");
121  program->addParamsLine(
122  " :+Lambda is a regularization factor used during the estimation of the CTF phase");
123  program->addParamsLine(
124  " :+During the estimation, the phase values are averaged within a window of this size");
125  program->addParamsLine(
126  " [--noDefocus] : No defocus estimation");
127  program->addParamsLine(
128  " [--defocus_range <D=8000>] : Defocus range in Angstroms");
129  program->addParamsLine(
130  " [--refine_amplitude_contrast] : Refine amplitude contrast with respect to the input one");
131  program->addParamsLine(
132  " [--show_optimization+] : Show optimization process");
133  program->addParamsLine(
134  " [--radial_noise++] : By default, noise is astigmatic");
135  program->addParamsLine(
136  " [--enhance_weight++ <w=1>] : Weight of the enhanced term");
137  program->addParamsLine(
138  " [--model_simplification++ <s=0>]: 0 (no simplification)");
139  program->addParamsLine(
140  " : 1 (simplified envelope)");
141  program->addParamsLine(
142  " : 2 (last Gaussian removal)");
143  program->addParamsLine(
144  " : 3 (symmetric intermediate Gaussian)");
145  program->addParamsLine(
146  " [--bootstrapFit++ <N=-1>] : Perform bootstrap fit (Fourier pixels are randomly chosen)");
147  program->addParamsLine(
148  " : This is used to test the variability of the fit");
149  program->addParamsLine(
150  " : N defines the number of times the fit is repeated");
151  program->addParamsLine("==+ CTF fit: Output CTF models");
152  program->addParamsLine(
153  " [--ctfmodelSize <size=256>] : Size for the ctfmodel thumbnails");
154  program->addParamsLine("==+ PSD enhancement");
155  program->addParamsLine(
156  " [--enhance_min_freq <f1>] : Bandpass cutoff. Normalized to 0.5");
157  program->addParamsLine(
158  " : If fmax>0.35, f1 default=0.01");
159  program->addParamsLine(
160  " : If fmax<0.35, f1 default=0.02");
161  program->addParamsLine(
162  " [--enhance_max_freq <f2>] : Bandpass cutoff. Normalized to 0.5.");
163  program->addParamsLine(
164  " [--selfEstimation] : Estimate defocus without previous estimation");
165  program->addParamsLine(
166  " : If fmax>0.35, f2 default=0.08");
167  program->addParamsLine(
168  " : If fmax<0.35, f2 default=0.15");
169 }
void addSeeAlsoLine(const char *seeAlso)
void addParamsLine(const String &line)

◆ defineParams()

void ProgCTFBasicParams::defineParams ( )
virtual

Define Parameters.

Reimplemented from XmippProgram.

Reimplemented in ProgCTFEstimateFromPSDFast.

Definition at line 171 of file ctf_estimate_from_psd_base.cpp.

172 {
173  addUsageLine("Adjust a parametric model to a PSD file.");
174  addUsageLine(
175  "The PSD is enhanced ([[http://www.ncbi.nlm.nih.gov/pubmed/16987671][See article]]). ");
176  addUsageLine(
177  "And finally, the CTF is fitted to the PSD, being guided by the enhanced PSD ");
178  addUsageLine(
179  "([[http://www.ncbi.nlm.nih.gov/pubmed/17911028][See article]]).");
180  addParamsLine(" --psd <PSDfile> : PSD file");
181  addSeeAlsoLine("ctf_estimate_from_micrograph");
182  defineBasicParams(this);
183 }
void addSeeAlsoLine(const char *seeAlso)
static void defineBasicParams(XmippProgram *program)
Define basic parameters.
void addUsageLine(const char *line, bool verbatim=false)
void addParamsLine(const String &line)

◆ generate_model_halfplane()

void ProgCTFBasicParams::generate_model_halfplane ( int  Ydim,
int  Xdim,
MultidimArray< double > &  model 
)

Generate half-plane model at a given size. It is assumed that ROUT_Adjust_CTF has been already run

◆ generate_model_quadrant()

void ProgCTFBasicParams::generate_model_quadrant ( int  Ydim,
int  Xdim,
MultidimArray< double > &  model 
)

Generate quadrant model at a given size. It is assumed that ROUT_Adjust_CTF has been already run

◆ produceSideInfo()

void ProgCTFBasicParams::produceSideInfo ( )

Produce side information.

Definition at line 186 of file ctf_estimate_from_psd_base.cpp.

187 {
188  // Resize the frequency
189  x_digfreq.initZeros(YSIZE(*f), XSIZE(*f) / 2);
190  y_digfreq.initZeros(YSIZE(*f), XSIZE(*f) / 2);
191  w_digfreq.initZeros(YSIZE(*f), XSIZE(*f) / 2);
192  w_digfreq_r.initZeros(YSIZE(*f), XSIZE(*f) / 2);
193  x_contfreq.initZeros(YSIZE(*f), XSIZE(*f) / 2);
194  y_contfreq.initZeros(YSIZE(*f), XSIZE(*f) / 2);
195  w_contfreq.initZeros(YSIZE(*f), XSIZE(*f) / 2);
196 
197  Matrix1D<int> idx(2); // Indexes for Fourier plane
198  Matrix1D<double> freq(2); // Frequencies for Fourier plane
199 
201  {
202  XX(idx) = j;
203  YY(idx) = i;
204 
205  // Digital frequency
206  FFT_idx2digfreq(*f, idx, freq);
207  x_digfreq(i, j) = XX(freq);
208  y_digfreq(i, j) = YY(freq);
209  w_digfreq(i, j) = freq.module();
210  w_digfreq_r(i, j) = (int)(w_digfreq(i,j) * (double)YSIZE(w_digfreq));
211 
212  // Continuous frequency
213  digfreq2contfreq(freq, freq, Tm);
214  x_contfreq(i, j) = XX(freq);
215  y_contfreq(i, j) = YY(freq);
216  w_contfreq(i, j) = freq.module();
217  }
218 
219  // Build frequency mask
223  {
224  if (w_digfreq(i, j) >= max_freq
225  || w_digfreq(i, j) <= min_freq)
226  continue;
227  mask(i, j) = 1;
228  w_count(w_digfreq_r(i, j))++;
229  }
230 
231  // Enhance PSD for ctfmodels
233  prm.filter_w1 = 0.02;
234  prm.filter_w2 = 0.2;
235  prm.decay_width = 0.02;
236  prm.mask_w1 = 0.01;
237  prm.mask_w2 = 0.5;
240  if (fn_psd.find('@')==std::string::npos)
241  enhanced_ctftomodel.write(fn_psd.withoutExtension() + "_enhanced_psd.xmp");
242  else
243  enhanced_ctftomodel.write(fn_psd.withoutExtension() + "_enhanced_psd.stk");
244  CenterFFT(enhanced_ctftomodel(), false);
246 
247  // Enhance PSD for optimization
248  prm.filter_w1 = f1;
249  prm.filter_w2 = f2;
252  CenterFFT(enhanced_ctftomodel(), false);
253  enhanced_ctftomodel().resize(w_digfreq);
254 
255  // Divide by the number of count at each frequency
256  // and mask between min_freq and max_freq
257  double min_val = enhanced_ctftomodel().computeMin();
259  if (mask(i, j) <= 0)
260  enhanced_ctftomodel(i, j) = min_val;
261 
264 
265  enhanced_ctftomodel() = aux;
266  enhanced_ctftomodel().rangeAdjust(0, 1);
267 
268  FourierFilter Filter;
269  Filter.FilterShape = RAISED_COSINE;
270  Filter.FilterBand = HIGHPASS;
271  Filter.w1 = 0.01;
272  Filter.raised_w = 0.005;
273  enhanced_ctftomodel().setXmippOrigin();
277 
278  // Compute now radial average of the enhanced_ctftomodel
286  {
288  {
289  int r = w_digfreq_r(i, j);
290  w_digfreq_r_iN(r)+=1;
292  psd_exp_radial(r) += ctftomodel(i, j);
293  }
294  }
296  if (w_digfreq_r_iN(i)>0)
297  {
301  }
302 
303  // Compute its derivative
304  int state=0;
305  double maxDiff=0;
307  {
308  switch (state)
309  {
310  case 0:
311  if (w_digfreq(i,0)>min_freq)
312  state=1;
313  break;
314  case 1:
315  state=2; // Skip first sample
316  break;
317  case 2:
318  if (w_digfreq(i,0)>max_freq)
319  state=3;
320  else
321  {
324  maxDiff=std::max(maxDiff,fabs(diff));
325  }
326  break;
327  }
328  }
330  if (show_optimization)
331  {
332  psd_exp_enhanced_radial.write("PPPexpEnhanced_fast.txt");
333  psd_exp_radial.write("PPPexp_fast.txt");
334  }
335 
336 }
#define YSIZE(v)
MultidimArray< double > * f
Image< double > ctftomodel
CTF amplitude to model.
double filter_w2
Bandpass filter high frequency (in Fourier space, max 0.5)
MultidimArray< double > w_digfreq
double mask_w1
Lower frequency for the mask (in Fourier space, max 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)
MultidimArray< double > y_digfreq
double f1
Enhancement filter low freq.
double Tm
Sampling rate.
MultidimArray< double > w_digfreq_r_iN
MultidimArray< double > psd_exp_enhanced_radial_derivative
void FFT_idx2digfreq(T &v, const Matrix1D< int > &idx, Matrix1D< double > &freq)
Definition: xmipp_fft.h:80
FileName fn_psd
CTF filename.
#define STARTINGX(v)
MultidimArray< double > w_digfreq_r
#define i
#define FOR_ALL_ELEMENTS_IN_ARRAY2D(m)
double max_freq
Maximum frequency to adjust.
double mask_w2
Higher frequency for the mask (in Fourier space, max 0.5)
#define STARTINGY(v)
void CenterFFT(MultidimArray< T > &v, bool forward)
Definition: xmipp_fft.h:291
Image< double > enhanced_ctftomodel_fullsize
CTF amplitude to model.
MultidimArray< double > w_contfreq
MultidimArray< double > psd_theo_radial
MultidimArray< double > psd_exp_radial
PSD data.
#define XX(v)
Definition: matrix1d.h:85
#define XSIZE(v)
void write(const FileName &fn) const
#define HIGHPASS
#define RAISED_COSINE
void max(Image< double > &op1, const Image< double > &op2)
double decay_width
Decay width (raised cosine)
double f2
Enhancement filter high freq.
bool show_optimization
Show convergence values.
MultidimArray< double > psd_theo_radial_derivative
#define j
#define YY(v)
Definition: matrix1d.h:93
MultidimArray< double > x_digfreq
#define FOR_ALL_ELEMENTS_IN_ARRAY1D(v)
FileName withoutExtension() const
double filter_w1
Bandpass filter low frequency (in Fourier space, max 0.5)
void medianFilter3x3(MultidimArray< T > &m, MultidimArray< T > &out)
Definition: filters.h:1088
void applyFilter(MultidimArray< double > &PSD)
double min_freq
Minimum frequency to adjust.
Image< double > enhanced_ctftomodel
CTF amplitude to model.
MultidimArray< double > x_contfreq
Frequencies in axes.
MultidimArray< double > psd_exp_enhanced_radial
ProgClassifyCL2D * prm
void initZeros(const MultidimArray< T1 > &op)
MultidimArray< double > mask
Masks.
void generateMask(MultidimArray< double > &v)
MultidimArray< double > w_count
MultidimArray< double > y_contfreq
void applyMaskSpace(MultidimArray< double > &v)
void digfreq2contfreq(const Matrix1D< double > &digfreq, Matrix1D< double > &contfreq, double pixel_size)
Definition: xmipp_fft.h:125

◆ readBasicParams()

void ProgCTFBasicParams::readBasicParams ( XmippProgram program)

Read parameters.

Definition at line 33 of file ctf_estimate_from_psd_base.cpp.

34 {
35  downsampleFactor = program->getDoubleParam("--downSamplingPerformed");
36  show_optimization = program->checkParam("--show_optimization");
37  min_freq = program->getDoubleParam("--min_freq");
38  max_freq = program->getDoubleParam("--max_freq");
39  defocus_range = program->getDoubleParam("--defocus_range");
40  modelSimplification = program->getIntParam("--model_simplification");
41  bootstrap = program->checkParam("--bootstrapFit");
42  refineAmplitudeContrast = program->checkParam("--refine_amplitude_contrast");
43  fastDefocusEstimate = program->checkParam("--fastDefocus");
44  noDefocusEstimate = program->checkParam("--noDefocus");
45  selfEstimation = program->checkParam("--selfEstimation");
47  {
48  lambdaPhase=program->getDoubleParam("--fastDefocus",0);
49  sizeWindowPhase=program->getIntParam("--fastDefocus",1);
50  }
51  ctfmodelSize = program->getIntParam("--ctfmodelSize");
52  enhanced_weight = program->getDoubleParam("--enhance_weight");
53  if (!program->checkParam("--enhance_min_freq"))
54  f1 = (max_freq > 0.35) ? 0.01 : 0.02;
55  else
56  f1 = program->getDoubleParam("--enhance_min_freq");
57  if (!program->checkParam("--enhance_max_freq"))
58  f2 = (max_freq > 0.35) ? 0.08 : 0.15;
59  else
60  f2 = program->getDoubleParam("--enhance_max_freq");
61 }
double defocus_range
Defocus range.
double enhanced_weight
Weight of the enhanced image.
double getDoubleParam(const char *param, int arg=0)
bool fastDefocusEstimate
Fast defocus estimate.
int sizeWindowPhase
Size of the average window used during phase direction and unwrapping estimates (used in Zernike esti...
double f1
Enhancement filter low freq.
double downsampleFactor
Downsample performed.
double max_freq
Maximum frequency to adjust.
int ctfmodelSize
X dimension of particle projections (-1=the same as the psd)
bool bootstrap
Bootstrap estimation.
bool noDefocusEstimate
No defocus estimate.
bool refineAmplitudeContrast
Refine amplitude contrast.
double f2
Enhancement filter high freq.
bool show_optimization
Show convergence values.
int modelSimplification
Model simplification.
double min_freq
Minimum frequency to adjust.
bool checkParam(const char *param)
double lambdaPhase
Regularization factor for the phase direction and unwrapping estimates (used in Zernike estimate) ...
int getIntParam(const char *param, int arg=0)

◆ readParams()

void ProgCTFBasicParams::readParams ( )
virtual

Read parameters.

Reimplemented from XmippProgram.

Reimplemented in ProgCTFEstimateFromPSDFast.

Definition at line 63 of file ctf_estimate_from_psd_base.cpp.

64 {
65  fn_psd = getParam("--psd");
66  readBasicParams(this);
67 }
void readBasicParams(XmippProgram *program)
Read parameters.
FileName fn_psd
CTF filename.
const char * getParam(const char *param, int arg=0)

◆ show()

void ProgCTFBasicParams::show ( )

Show parameters.

Definition at line 70 of file ctf_estimate_from_psd_base.cpp.

71 {
72  if (verbose==0)
73  return;
74  std::cout
75  << "PSD file: " << fn_psd << std::endl
76  << "Downsampling: " << downsampleFactor << std::endl
77  << "Min Freq.: " << min_freq << std::endl
78  << "Max Freq.: " << max_freq << std::endl
79  << "Sampling: " << Tm << std::endl
80  << "Defocus range: " << defocus_range << std::endl
81  << "ctfmodelSize: " << ctfmodelSize << std::endl
82  << "Enhance min freq: " << f1 << std::endl
83  << "Enhance max freq: " << f2 << std::endl
84  << "Enhance weight: " << enhanced_weight << std::endl
85  << "Model simplification: " << modelSimplification << std::endl
86  << "Bootstrap: " << bootstrap << std::endl
87  << "Fast defocus: " << fastDefocusEstimate << std::endl
88  << "No defocus: " << noDefocusEstimate << std::endl
89  << "Refine amplitude contrast: " << refineAmplitudeContrast << std::endl
90  ;
92  std::cout
93  << "Regularization factor: " << lambdaPhase << std::endl
94  << "Window size: " << sizeWindowPhase << std::endl;
95 
96 }
double defocus_range
Defocus range.
double enhanced_weight
Weight of the enhanced image.
bool fastDefocusEstimate
Fast defocus estimate.
int sizeWindowPhase
Size of the average window used during phase direction and unwrapping estimates (used in Zernike esti...
double f1
Enhancement filter low freq.
double Tm
Sampling rate.
double downsampleFactor
Downsample performed.
FileName fn_psd
CTF filename.
double max_freq
Maximum frequency to adjust.
int ctfmodelSize
X dimension of particle projections (-1=the same as the psd)
bool bootstrap
Bootstrap estimation.
bool noDefocusEstimate
No defocus estimate.
bool refineAmplitudeContrast
Refine amplitude contrast.
int verbose
Verbosity level.
double f2
Enhancement filter high freq.
int modelSimplification
Model simplification.
double min_freq
Minimum frequency to adjust.
double lambdaPhase
Regularization factor for the phase direction and unwrapping estimates (used in Zernike estimate) ...

Member Data Documentation

◆ action

int ProgCTFBasicParams::action

Definition at line 102 of file ctf_estimate_from_psd_base.h.

◆ adjust

Matrix1D<double> ProgCTFBasicParams::adjust

Set of parameters for the complete adjustment of the CTF.

Definition at line 62 of file ctf_estimate_from_psd_base.h.

◆ adjust_params

Matrix1D<double>* ProgCTFBasicParams::adjust_params

Definition at line 126 of file ctf_estimate_from_psd_base.h.

◆ amplitud

std::vector<double> ProgCTFBasicParams::amplitud

Definition at line 83 of file ctf_estimate_from_psd_base.h.

◆ bootstrap

bool ProgCTFBasicParams::bootstrap

Bootstrap estimation.

Definition at line 36 of file ctf_estimate_from_psd_base.h.

◆ corr13

double ProgCTFBasicParams::corr13

Definition at line 112 of file ctf_estimate_from_psd_base.h.

◆ ctfmodelSize

int ProgCTFBasicParams::ctfmodelSize

X dimension of particle projections (-1=the same as the psd)

Definition at line 34 of file ctf_estimate_from_psd_base.h.

◆ ctftomodel

Image<double> ProgCTFBasicParams::ctftomodel

CTF amplitude to model.

Definition at line 25 of file ctf_estimate_from_psd_base.h.

◆ current_penalty

double ProgCTFBasicParams::current_penalty

Definition at line 122 of file ctf_estimate_from_psd_base.h.

◆ defocus_range

double ProgCTFBasicParams::defocus_range

Defocus range.

Definition at line 54 of file ctf_estimate_from_psd_base.h.

◆ downsampleFactor

double ProgCTFBasicParams::downsampleFactor

Downsample performed.

Definition at line 23 of file ctf_estimate_from_psd_base.h.

◆ enhanced_ctftomodel

Image<double> ProgCTFBasicParams::enhanced_ctftomodel

CTF amplitude to model.

Definition at line 27 of file ctf_estimate_from_psd_base.h.

◆ enhanced_ctftomodel_fullsize

Image<double> ProgCTFBasicParams::enhanced_ctftomodel_fullsize

CTF amplitude to model.

Definition at line 29 of file ctf_estimate_from_psd_base.h.

◆ enhanced_weight

double ProgCTFBasicParams::enhanced_weight

Weight of the enhanced image.

Definition at line 60 of file ctf_estimate_from_psd_base.h.

◆ evaluation_reduction

int ProgCTFBasicParams::evaluation_reduction

Definition at line 118 of file ctf_estimate_from_psd_base.h.

◆ f

MultidimArray<double>* ProgCTFBasicParams::f

Definition at line 125 of file ctf_estimate_from_psd_base.h.

◆ f1

double ProgCTFBasicParams::f1

Enhancement filter low freq.

Definition at line 56 of file ctf_estimate_from_psd_base.h.

◆ f2

double ProgCTFBasicParams::f2

Enhancement filter high freq.

Definition at line 58 of file ctf_estimate_from_psd_base.h.

◆ fastDefocusEstimate

bool ProgCTFBasicParams::fastDefocusEstimate

Fast defocus estimate.

Definition at line 40 of file ctf_estimate_from_psd_base.h.

◆ fn_psd

FileName ProgCTFBasicParams::fn_psd

CTF filename.

Definition at line 21 of file ctf_estimate_from_psd_base.h.

◆ heavy_penalization

double ProgCTFBasicParams::heavy_penalization

Definition at line 121 of file ctf_estimate_from_psd_base.h.

◆ lambdaPhase

double ProgCTFBasicParams::lambdaPhase

Regularization factor for the phase direction and unwrapping estimates (used in Zernike estimate)

Definition at line 44 of file ctf_estimate_from_psd_base.h.

◆ mask

MultidimArray<double> ProgCTFBasicParams::mask

Masks.

Definition at line 85 of file ctf_estimate_from_psd_base.h.

◆ mask_between_zeroes

MultidimArray<double> ProgCTFBasicParams::mask_between_zeroes

Definition at line 86 of file ctf_estimate_from_psd_base.h.

◆ max_freq

double ProgCTFBasicParams::max_freq

Maximum frequency to adjust.

Definition at line 50 of file ctf_estimate_from_psd_base.h.

◆ max_freq_psd

double ProgCTFBasicParams::max_freq_psd

Definition at line 92 of file ctf_estimate_from_psd_base.h.

◆ max_gauss_freq

double ProgCTFBasicParams::max_gauss_freq

Definition at line 95 of file ctf_estimate_from_psd_base.h.

◆ min_freq

double ProgCTFBasicParams::min_freq

Minimum frequency to adjust.

Definition at line 48 of file ctf_estimate_from_psd_base.h.

◆ min_freq_psd

double ProgCTFBasicParams::min_freq_psd

Definition at line 91 of file ctf_estimate_from_psd_base.h.

◆ modelSimplification

int ProgCTFBasicParams::modelSimplification

Model simplification.

Definition at line 64 of file ctf_estimate_from_psd_base.h.

◆ noDefocusEstimate

bool ProgCTFBasicParams::noDefocusEstimate

No defocus estimate.

Definition at line 42 of file ctf_estimate_from_psd_base.h.

◆ penalize

bool ProgCTFBasicParams::penalize

Definition at line 115 of file ctf_estimate_from_psd_base.h.

◆ penalty

constexpr double ProgCTFBasicParams::penalty = 32.0
static

Definition at line 123 of file ctf_estimate_from_psd_base.h.

◆ psd_exp_enhanced_radial

MultidimArray<double> ProgCTFBasicParams::psd_exp_enhanced_radial

Definition at line 74 of file ctf_estimate_from_psd_base.h.

◆ psd_exp_enhanced_radial_2

MultidimArray<double> ProgCTFBasicParams::psd_exp_enhanced_radial_2

Definition at line 75 of file ctf_estimate_from_psd_base.h.

◆ psd_exp_enhanced_radial_derivative

MultidimArray<double> ProgCTFBasicParams::psd_exp_enhanced_radial_derivative

Definition at line 76 of file ctf_estimate_from_psd_base.h.

◆ psd_exp_radial

MultidimArray<double> ProgCTFBasicParams::psd_exp_radial

PSD data.

Definition at line 73 of file ctf_estimate_from_psd_base.h.

◆ psd_exp_radial_derivative

MultidimArray<double> ProgCTFBasicParams::psd_exp_radial_derivative

Definition at line 78 of file ctf_estimate_from_psd_base.h.

◆ psd_fft

MultidimArray<std::complex<double> > ProgCTFBasicParams::psd_fft

Definition at line 82 of file ctf_estimate_from_psd_base.h.

◆ psd_theo_radial

MultidimArray<double> ProgCTFBasicParams::psd_theo_radial

Definition at line 79 of file ctf_estimate_from_psd_base.h.

◆ psd_theo_radial_derivative

MultidimArray<double> ProgCTFBasicParams::psd_theo_radial_derivative

Definition at line 77 of file ctf_estimate_from_psd_base.h.

◆ refineAmplitudeContrast

bool ProgCTFBasicParams::refineAmplitudeContrast

Refine amplitude contrast.

Definition at line 38 of file ctf_estimate_from_psd_base.h.

◆ selfEstimation

bool ProgCTFBasicParams::selfEstimation

Definition at line 32 of file ctf_estimate_from_psd_base.h.

◆ show_inf

int ProgCTFBasicParams::show_inf

Definition at line 97 of file ctf_estimate_from_psd_base.h.

◆ show_optimization

bool ProgCTFBasicParams::show_optimization

Show convergence values.

Definition at line 31 of file ctf_estimate_from_psd_base.h.

◆ sizeWindowPhase

int ProgCTFBasicParams::sizeWindowPhase

Size of the average window used during phase direction and unwrapping estimates (used in Zernike estimate)

Definition at line 46 of file ctf_estimate_from_psd_base.h.

◆ Tm

double ProgCTFBasicParams::Tm

Sampling rate.

Definition at line 52 of file ctf_estimate_from_psd_base.h.

◆ value_th

double ProgCTFBasicParams::value_th

Definition at line 90 of file ctf_estimate_from_psd_base.h.

◆ w_contfreq

MultidimArray<double> ProgCTFBasicParams::w_contfreq

Definition at line 68 of file ctf_estimate_from_psd_base.h.

◆ w_count

MultidimArray<double> ProgCTFBasicParams::w_count

Definition at line 87 of file ctf_estimate_from_psd_base.h.

◆ w_digfreq

MultidimArray<double> ProgCTFBasicParams::w_digfreq

Definition at line 71 of file ctf_estimate_from_psd_base.h.

◆ w_digfreq_r

MultidimArray<double> ProgCTFBasicParams::w_digfreq_r

Definition at line 81 of file ctf_estimate_from_psd_base.h.

◆ w_digfreq_r_iN

MultidimArray<double> ProgCTFBasicParams::w_digfreq_r_iN

Definition at line 80 of file ctf_estimate_from_psd_base.h.

◆ x_contfreq

MultidimArray<double> ProgCTFBasicParams::x_contfreq

Frequencies in axes.

Definition at line 66 of file ctf_estimate_from_psd_base.h.

◆ x_digfreq

MultidimArray<double> ProgCTFBasicParams::x_digfreq

Definition at line 69 of file ctf_estimate_from_psd_base.h.

◆ y_contfreq

MultidimArray<double> ProgCTFBasicParams::y_contfreq

Definition at line 67 of file ctf_estimate_from_psd_base.h.

◆ y_digfreq

MultidimArray<double> ProgCTFBasicParams::y_digfreq

Definition at line 70 of file ctf_estimate_from_psd_base.h.


The documentation for this class was generated from the following files: