User Tools

Site Tools


pipeline:meridien:sxmeridien
This version is outdated by a newer approved version.DiffThis version (2018/08/22 11:02) is a draft.
Approvals: 0/1

This is an old revision of the document!


sxmeridien

3D Refinement: Performs 3D structure refinement using a quasi-Maximum Likelihood approach.


Usage

Usage in command line

sxmeridien.py  stack  output_directory  initial_volume  --do_final=MERIDIEN_ITERATION_ID  --local_refinement  --radius=particle_radius  --mask3D=MASK3D  --symmetry=SYMMETRY  --inires=INITIAL_RESOLUTION  --delta=DELTA  --xr=XR  --ts=TS  --initialshifts  --skip_prealignment  --memory_per_node=MEMORY_PER_NODE  --center_method=CENTER_METHOD  --target_radius=TARGET_RADIUS  --an=ANGULAR_NEIGHBORHOOD  --shake=SHAKE  --small_memory  --ccfpercentage=CCFPERCENTAGE  --nonorm  --function=USER_FUNCTION


Typical usage

sxmeridien exists only in MPI version (Running MPI version does not require –MPI flag).

There are five ways to run the program:


1. Standard default run:
The standard refinement starts from exhaustive searches, uses initial reference structure

mpirun -np 64 --hostfile four_nodes.txt  sxmeridien.py  bdb:sparx_stack vton1 mask15.hdf --sym=c5  --initialshifts  --radius=120  --mask3D=mask15.hdf    >1ovotn &


2. Restart after the last fully finished iteration: One can change some parameters (MPI settings have to be the same)

mpirun -np 64 --hostfile four_nodes.txt  sxmeridien.py  vton1 --radius=100 >2ovotn &


3. Local refinement:
Local refinement starts from user-provided orientation parameters, delta has to be ⇐ 3.75

mpirun -np 64 --hostfile four_nodes.txt sxmeridien.py --local_refinement bdb:sparx_stack   vton3 --delta=1.875 --xr=2.0  --inires=5.5  --sym=c5  --radius=120  --mask3D=mask15.hdf >5ovotn &


4. Restart of local refinement after the last fully finished iteration.:
One can change some parameters (MPI settings have to be the same)

mpirun -np 64 --hostfile four_nodes.txt  sxmeridien.py --local_refinement  vton3  --xr=0.6 >6ovotn &


5. Final reconstruction only:
Do only final reconstruction using a fully finished iteration of meridien (here number 21).

mpirun -np 64 sxmeridien.py --do_final=21  meridien_outdir


Input

Main Parameters

stack
Input image stack: Input image stack. Required for new runs, not for restarts. (default none)
output_directory
Output directory: The results will be written here. If not given, the program will use name master_DATA_AND_TIME. For standard continuation run, local refinement from iteration, and final reconstruction only, the directory must exist and further iterations will be written to this directory. (default none)
initial_volume
Initial 3D reference: The 3D reference used in the initial iteration of 3D refinement. Required only for new runs. (default none)
--do_final
Do only final reconstruction: Specify the iteration which to perform final reconstruction. By setting to 0, program searches for the iteration that had best resolution, then performs correponding final reconstruction. Value must be zero or positive. (default -1)
--local_refinement
Perform local refinement: Perform local refinement starting from user-provided orientation parameters stored in the header of input image stack. (default False)
--radius
Particle radius [Pixels]: Outer particle radius in pixels < int(boxsize/2)-1. Ignored in final reconstruction. (default -1)
--mask3D
3D mask file: A mask applied to half-map during iterations of the program. If not given, a hard sphere of radius boxsize/2-1 will be used. Ignored in final reconstruction. (default none)
--symmetry
Point-group symmetry: Point-group symmetry of the refined structure. Supported point groups symmetries are: cn and dn, where n is multiplicity, oct, tet, and icos. Ignored in final reconstruction. (default c1)
--inires
Starting resolution [A]: Resolution of the initial map used to start the refinement. Ignored in final reconstruction. (default 25.0)
--delta
Initial angular sampling step [Degrees]: Initial angular sampling step. Ignored in final reconstruction. (default 7.5)
--xr
Search range [Pixels]: Range for translation search in both directions. Search is +/-xr. It can be fractional. Ignored in final reconstruction. (default 5.0)
--ts
Search step size [Pixels]: Step size of translation search in both directions. Search is within a circle of radius xr on a grid with steps ts. It can be fractional. (default 1.0)
--initialshifts
Read shifts from header: Start refinement using translation parameters located in the input file header to jumpstart the procedure. Specific to standard run mode. (default False value reversed in GUI)
--skip_prealignment
Skip the 2D pre-alignment step: Indicate if pre-alignment should be omitted. By default it is performed in order to find initial translation parameters. This accelerates convergence of the program. Do not use 2D pre-alignment if images are already centered. Specific to standard run modes. (default False question reversed in GUI)
--memory_per_node
Memory per node [GB]: User provided information about memory per node in GB (NOT per CPU). By default, it uses 2GB * (number of CPUs per node). Used in all modes. (default -1.0)


Advanced Parameters

--an
Angular neighborhood: Angular neighborhood for local search. Used only in Local Refinement mode. Ignored in final reconstruction. (default -1.0)
--center_method
Centering method: Method for centering averages during initial 2D prealignment of data (0: no centering; -1: average shift method; For 1-7, see center_2D in utilities.py). Specific to standard run modes. (default -1)
--target_radius
Target particle radius [Pixels]: For 2D prealignment, images will be shrank or enlarged to this radius. Specific to standard run modes. (default 29)
--shake
Shake: Shake randomizes grid searches by a factor of (shake x grid step). Ignored in final reconstruction. (default 0.5)
--small_memory
Keep data in memory: Indicate if data should be kept in memory or not. By default, data will be kept in memory. Ignored in final reconstruction. (default False question reversed in GUI)
--ccfpercentage
Correlation peaks to be included [%]: Percentage of correlation peaks to be included. 0.0 corresponds to hard matching. Ignored in final reconstruction. (default 99.9)
--nonorm
Apply image norm correction: Indicate if image norm correction should be applied or not. By default, apply image norm correction. Ignored in final reconstruction. (default False question reversed in GUI)
--function
Reference preparation function: Specify name of function that program should use to prepare the reference structure after each iteration. Ignored in final reconstruction. (default do_volume_mask)


Output


Description


Method


Reference


Developer Notes


Author / Maintainer

Pawel A. Penczek


Keywords

Category 1:: APPLICATIONS


Files

sparx/bin/sxmeridien.py


See also

References to relevant other routines.


Maturity

Alpha:: Fully developed.


Bugs

Glitches possible.


pipeline/meridien/sxmeridien.1534928538.txt.gz · Last modified: 2018/08/22 11:02 by fmerino