Commit 569acab4 authored by myron's avatar myron

changed the name to the double flat field

parent ff90cd1c
......@@ -5863,7 +5863,7 @@ void CCspace_InterleavedReadingPreProcessing_chunk(CCspace * self ,int sn,int
// Write this in edf
{
char name[10000];
sprintf(name,"%s/projectionsmean_%20.10f.edf", dirname(self->params.OUTPUT_FILE) , self->params.FF2_SIGMA ) ;
sprintf(name,"%s/pyhst_dff_%.2f.edf", dirname(self->params.OUTPUT_FILE) , self->params.FF2_SIGMA ) ;
write_data_to_edf(self->ff2_globalmean, Size0, Size1, name);
}
......@@ -5968,7 +5968,7 @@ void CCspace_InterleavedReadingPreProcessing_chunk(CCspace * self ,int sn,int
{
char name[10000];
sprintf(name,"%s/projectionsmean_%20.10f.edf", dirname(self->params.OUTPUT_FILE) , self->params.FF2_SIGMA ) ;
sprintf(name,"%s/pyhst_dff_%.2f.edf", dirname(self->params.OUTPUT_FILE) , self->params.FF2_SIGMA ) ;
write_data_to_edf(self->ff2_globalmean, Size0, Size1, name);
if (self->params.verbosity >= 10) puts("wrote projectionsmean.edf");
......
......@@ -940,7 +940,7 @@ class Parameters:
several causes, for example a detector non-linearity.
If let to zero, no correction is done. If set to 1 then
one image is read from the double-ff EDF file name projectionsmean.edf residing in the same directory as the output file.
one image is read from the double-ff EDF file name pyhst_dff_1.23.edf residing in the same directory as the output file.
If the file is not there, it is generated first and then used for the reconstruction
The image must contain the logarithm of the correction factor. The data is
subtracted from the treated projection.
......@@ -952,12 +952,10 @@ class Parameters:
Mind that if DOUBLEFFCORRECTION_ONTHEFLY > 0, then we should have DOUBLEFFCORRECTION = 0.
In the current version, this parameter does not work for helical tomography.
If DOUBLEFFCORRECTION_ONTHEFLY is not zero, all the projections are averaged.
A file "projectionsmean.edf" is generated.
A file "pyhst_dff_0.12.edf" is generated.
This file is used as a double flat-field correction.
If FF2_SIGMA > 0, the average is high-pass filtered, where the filter is a complementary Gaussian function of STD FF2_SIGMA.
A file "projectionsmean_filt.edf" is generated.
This filtered average is used as a double flat field instead of the plain average.
"""
FF2_SIGMA = 0.0
......@@ -2052,7 +2050,7 @@ Units are pixels per projection.
DOUBLERUN4DFF=0
if DOUBLEFFCORRECTION:
dname = os.path.dirname(OUTPUT_FILE)
dd_path = os.path.join( dname,"projectionsmean_%20.10f.edf"%FF2_SIGMA)
dd_path = os.path.join( dname,"pyhst_dff_%.2f.edf"%FF2_SIGMA)
print( " CONTROLLO ", dd_path)
if not os.path.exists(dd_path):
DOUBLEFFCORRECTION = 0
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment