User Tools

Site Tools

This version (2021/10/11 17:15) is a draft.
Approvals: 0/1
The Previously approved version (2021/02/01 17:22) is available.Diff

sp_signalsubtract centershift

Center Map : Center map of remaining density.


Usage in command line centershift --cvol1=map_center_1 --cvol2=map_center_2 --shiftparams=shift_parameters --diffimgs=images_to_subtract --outdir=output_subtraction_directory --volradius=map_radius --apix=pixel_size --verbosity=verbosity_level

Typical usage

sp_signalsubtract comprises the steps needed to subtract masked regions from a map in order to focus on the remaining features from that map. There are six steps for this process (five within SPHIRE):

1. Average and low-pass filter a map for segmentation

2. Locate regions to be subtracted

3. Mask a soft mask of the region to be excluded

4. Compute re-projections of map to be subtracted, and subtract them from the original images

5. Compute reconstruction of density remaining after signal-subtraction

6. Center map of the remaining density : --mode centershift --cvol1=centered_halfmap_1 --cvol2=centered_halfmap_2 --shiftparams=meridien_parameters --diffimgs bdb:output_subtraction_directory#all_subtracted --outdir=output_subtraction_directory --apix=pixel_size --volradius=structure_radius

As inputs “–cvol1” and “–cvol2”, use the half-maps generated by the previous reconstruction step. For “–shiftparams”, use the Meridien parameters for the iteration used for those reconstructions. e.g., of the form “final_params_015.txt”. For “–diffimgs”, use the signal-subtracted image stack from above.

The reconstruction will be roughly centered by alignment to a Gaussian blob of radius specified by the flag “–volradius”. (Additional centering options will be provided in the future.) If the pixel size is provided (–apix), the units of radius will be assumed to be in Angstroms. If the pixel size is not provided, the radius will be assumed to be in units of voxels. If your structure is not globular, try the shortest dimension.

The output centered parameter file will be combine the Meridien parameters with the centering parameters. The output reconstruction will not use smear information, and may be noisier than the reconstruction from the previous step. Simply make sure that the centering was adequate. The better the centering, the smaller the required shifts will be in subsequent refinements.


Main Parameters

First map to average and center. (default required string)
Second map to average and center. (default None)
Meridien parameters to combine with centering parameters. (default None)
Images from signal subtraction to reconstruct. (default None)
Directory where outputs will be written. (default required string)
Radius to use for centering reconstruction. If the structure is not globular, try the shortest dimension. (default required string)
Pixel size in Angstroms. (default None)

Advanced Parameters

Verbosity level
Controls how much information will be written to the screen. (default 2)


Reconstruction of subtracted particles
Negative of above reconstruction, subtracted region shouldn't have significant densities
Centered reconstruction
Original parameters
Parameters after centering
Reconstruction (filtered) before centering
Reconstruction after centering




Developer Notes

Author / Maintainer

Tapu Shaikh





See also


Beta:: Under evaluation and testing. Please let us know if there are any bugs.


pipeline/subtract/centershift.txt · Last modified: 2021/10/11 17:15 by shaikh