HOME GUIDE OPERATIONS DOCS ERRORS FORMATS INSTALL NEW TIPS WEB SITES

RM ORF - Orientational search between a 3d and 2d Radon transform.

(3/19/90)

PURPOSE

Calculates the angular crosscorelation function in three Euler angles between a 3d Radon transform and 2d Radon transform. Simultaneous translational alignment is possible. (This is the original program. It has been replaced by a set of more efficient programs: Combination of RM FT and RM ORF)

USAGE

.OPERATION: RM ORF X11,X12,X13,X14,X15,X16

[The optional registers will return: Phi,Theta,Psi,Value of CCF maximum,x-shift,y-shift of last aligned 2D-transform.]

.DEBUG OUTPUT: 0
[Answer: 0 otherwise very large amounts of output are produced]

.TEMPLATE FOR 2D RADON OUTPUT: *
[If * is answered no output is produced. Otherwise every Radon transform (i.e. one for every combination of theta and phi) that is extracted from the 3D Radon transform is written to disk.]

.TEMPLATE FOR CCF OUTPUT: *
[If * is answered no ouptu is produced. Otherwise for each combination of theta and phi one cross-correlation function between the input and the extracted Radon transform is written to disk.]

If a file name was given for CCF output:

.PSI NUMBER OF CCF: 11
[Anwer the number of the angle psi for which the cross-corelation function is written. Ths angle psi depends on the start value for the psi-search and the psi increment.]

.3D RADON FILE: 3DRAD001
[Enter the name of the 3D Radon transform.]

.ENTER FILE TEMPLATE: RAD2D***
[Enter example name of two-dimensional input Radon transform.]
.ENTER FILE NUMBERS: 1-188
[Enter file numbers of input 2D transforms.]

.SAMPLING RANGE AND INCREMENT IN EXTRACTED RADON TRANSFORM:

.ALPHA FROM,TO (DEF=2D): -90,88

.INCREMENT (DEF=2D: 2
[Enter parameter for the Radon transforms that are extracted from the 3D transform. The default is the same sampling as the 2D Radon transform.]

.PEAK FILE: peak001
[Enter name of file for the 3D angular cross-coraltion function.]

.(S)HIFTALIGN,(N)O SHIFTALIGN: S
[Enter S if a simultaneous translational alignment is desired.]

if "S" was answered:

.MAXIMUM SHIFT (IN INTEGER PIXELS): 5
[Enter radius for translational search.]

.(S)UBSEARCH OR (C)OMPLETE SEARCH?(=DEF): C
[For C a search of the ranges specified below is executed. If S is answered, the ranges specified below are centered around angles that are read from the file header.]

if C was answered:

.POSITION FOR WRITING ANGLES (1-3):
[Enter header position were the alignment results will be kept in each 2D file. (s. explanation for S-option below.]

if S was answerd:

.POSITION FOR READING, WRITING ANGLES (1-3): 1,2
[Enter from which file header location to read the current orienatation/translation if a subsearch is done, and where to write the results of the alignment. Position 1 is buffer location 15, position 2 is buffer location 31, position 3 is buffer location 36. The values stored are: phi,theta,psi,x-shift,y-shift.]

.THETA RANGE OF SEARCH (FROM,TO): -90,80

.INCREMENT: 10
[Enter search range and incement for theta (second Euler angle around y-axis.]

.PHI RANGE OF SEARCH (FROM,TO): -180,170

.INCRMENT: 10
[Enter search range and increment of first Euler Angle. (around z-axis.]

.PSI RANGE OF SEARCH (FROM,TO): -180,170

.INCREMENT: 10
[Enter search range of third Euler angle (around z'-axis, within the 2D projection plane.]

.DOCUMENT FILE: ALIDOC001
[Enter name of document file where the alignment parameters of each 2D-transform are stored. Key: counter, column 1:file number, column 2: phi, column 3: theta, column 4: psi, column 5:x-shift, column 6: y-shift.]

NOTES

  1. For more details see: M.Radermacher, Ultramicroscopy 53, 1994, 121-136

  2. Developed by M.Radermacher

  3. Not Distributed!

SUBROUTINES

MRRMOR, MRCC1F, MRFOUR1

CALLER

UTIL5

© Copyright Notice /       Enquiries: spider@wadsworth.org