40 addUsageLine(
"Reduce the number of pseudoatoms in a volume.");
41 addExampleLine(
" xmipp_pdb_reduce_pseudoatoms -i 1o7d.vol -o 1o7dreduced.vol --threshold 0.15");
45 addParamsLine(
" [--number <num=-1.0>] : Number of pseudoatoms to keep based on intensity");
46 addParamsLine(
" [--threshold <thresh=0.0>] : Intensity threshold for removing pseudoatoms");
61 std::cout <<
"PDB file: " << fn_volume << std::endl
62 <<
"Number: " << num << std::endl
63 <<
"Threshold: " << thresh << std::endl;
70 pdb.
read(fn_volume,
false, thresh);
73 pdb.
read(fn_volume,
true);
78 pdb.
read(fn_volume,
false, thresh);
80 pdb.
write(fn_out,
true);
void read(const FileName &fnPDB, const bool pseudoatoms=false, const double threshold=0.0)
Read rich phantom from either a PDB of CIF file.
double getDoubleParam(const char *param, int arg=0)
void write(const FileName &fnPDB, const bool renumber=false)
Write rich phantom to PDB or CIF file.
std::vector< double > intensities
void reduceNumberPseudoatoms()
const char * getParam(const char *param, int arg=0)
void addExampleLine(const char *example, bool verbatim=true)
void sort(struct DCEL_T *dcel)
std::vector< RichAtom > atomList
List of atoms.
std::vector< std::string > remarks
List of remarks.
FileName withoutExtension() const
bool checkParam(const char *param)
void addUsageLine(const char *line, bool verbatim=false)
void addParamsLine(const String &line)