RB 32F - Rotate, Shift, Sampled, Back Projection in Fourier space |

PURPOSE

Changes the scale, rotates, and shifts image circularly. Then calculates two randomly selected sample recontruction and a total-3D reconstruction using interpolation in Fourier space.Rotates counter-clockwise around the center (NSAM/2 + 1, NROW/2 + 1). (Negative angles = clockwise. Note that the terms "clockwise" and "counter-clockwise" refer to the mirrored x-y system used for image display). This operation is the same as 'RT SQ' followed by 'BP 32F.

SEE ALSO

RT SQ	[Rotate and Shift - image, Quadratic interpolation ||]
BP 32F	[Back Projection - 3D, Sampled, Interpolated in Fourier space ||*]
BP 3F	[Back Projection - 3D, Interpolated in Fourier space ||*]

USAGE

.OPERATION: RB 32F

.TEMPLATE FOR IMAGE FILES (E.G. STK@****): INSTACK03@
[Enter template for the names of the file(s) to be rotated. This template can alternatively refer to a series of single images (e.g. img***), or to a series of stacked images (e.g stack3@***), or to a whole stack of images ending in '@' (e.g. stk@). Templates containing an asterisk will have the astericks substituted with a number for the file.]

If the template contains an asterisk the following question will appear:

.ENTER FILE NUMBERS OR SELECTION DOC. FILE NAME: 1-700 OR IMAGE_SEL
[Enter projection file numbers or the name of the document file containing projection file numbers in the first column.]

.ANGLES DOC FILE: ANGLES27
[Enter the name of the document file containing Eulerian angles for the projections used (psi, theta, phi).]

.SYMMETRIES DOC FILE: symmetries
[Enter the name of the document file containing Eulerian angles defining symmetries (psi, theta, phi). The angles should be such that when used in operation 'RT 3D' transform the volume into itself. The identity operation (angles 0,0,0) MUST be included.
The symmetries are enforced in two steps: projections related by the symmetry operations are added in Fourier space to the volume; in order to reduce the residual artifacts, additional real-space symmetrization is applied.
If there are no symmetries, enter symbol "*".]

.TRANSFORMED OUTPUT IMAGES TEMPLATE (E.G. ROT@****)~',
[Enter template for the names of the rotated and shifted output image files. This template can alternatively refer to a series of single images (e.g. img***), or to a series of stacked images (e.g stack3@***), or to a whole stack of images ending in '@' (e.g. stk@). Templates containing an asterisk will have the astericks substituted with a number for the file. However a single asterick will skip any output for these images.]

If this output template contains an asterick the following question appears:

.ENTER FILE NUMBERS OR SELECTION DOC. FILE NAME: select_03
[Enter the file numbers to be substituted into the output file name template or name of a selection document file that contains the file numbers in the first register.]

The following questions always appear:

.RECONSTRUCTED VOLUME: VOLUME001
[Name for the backprojected output volume calculated using all input projections specified in the selection document file.]

.FIRST SAMPLE VOLUME: VOLUME1001
[Name for the backprojected output volume calculated from a randomly selected set formed from half of the input projections.]

.SECOND SAMPLE VOLUME: VOLUME2001
[Name for the backprojected output volume calculated from remaining half of the input projections.] .ENTER REG. NUMBERS FOR ANGLE, SCALE, X, & Y SHIFT: 1 2 3 4 [Enter the column positions in the document file where these parameters are stored. If scale register is zero the scale defaults to 1.0.]

OTES: 1. A change of scale will not change the image dimensions.

All the parameters agree with the output of the marker program (operation 'MK').

2. The input file may be three-dimensional. In this case,

each of the slices is separately rotated and shifted with respect to the z-axis, and stored in the output file, which will be created with the same dimensions.

3. Shifts can be any real numbers.

5. The accuracy of interpolation is much higher than that of the

bilinear interpolation used in the 'RT' operation, but it is lower than the accuracy of bi-spline interpolation. In particular, the derivatives of the interpolated image are not continuous. The speed, however, is much higher than bi-spline interpolation.

6. Since quadratic interpolation is used, it is possible for

the output image min. & max. intensities to be outside the range of the input image.

7. 5x5x5 Keiser-Bessel interpolation scheme is used.

8. Sample volumes can be used to evaluate the resolution. The division of

projection data set is performed in a random manner. Therefore, the two sample volumes will be different if the command is repeated on the same data set.

9. If you do not have enough memory in your computer to use

this operation substitute 'RB 3F' instead and invoke it three times to create the various output volumes.

10. This operation is NOT parallelized for use with MPI.

11. Implemented by Paul Penczek and ArDean Leith.

SUBROUTINES

WIW32D, ROTQS, ROT2QS, QUADRI, GETDOCDAT

CALLER

VTIL2