Exp_Shift_Computation

Purpose

This program computes experimental shifts at the unit cell center from the experimental correlation peaks stored in[filein]. Only points with crosscorrelation coefficient above[cc_peak_factor]*Maximum-Correlation are considered. The former experimental are used to linearly interpolate shifts to the each lattice point. The result is stored in[fileout]

Usage


$ exp_shift_computation ...

Parameters

`` Input correlation peaks file (MRC quadsearch output)
`` Output experimental shifts file
`` Crosscorrelation thershold (0-1) (as in MRC ccunbend)

Examples and notes

Computation of shifts using peaks in correlation file0001PASS1.cor with correlation above0.13*Maximum-Correlation. The computed experimental shifts are stored in[[ExpShifts]].doc


$ xmipp_exp_shift_computation -i 0001PASS1.cor -o ExpShifts.doc -cc_peak_factor 0.13

The file format for experimental shifts ouput is as follows:


; Dimensions: dimx dimy Ox Oy
; 4096 4096 2048.000000 2048.000000
; Lat. Vectors: ax ay bx by
; 47.874779  91.586548  91.066162  -48.395161
; NumFila NumElem  h   k  X   Y IncrementoX  IncrementoY
      1  6   3 -24     6.04  3484.24   -3.686   -9.694
      2  6  -2 -21    39.86  2881.13   -0.461   -6.025
      3  6   2 -23    49.23  3344.26   -2.728   -8.412
      4  6 -17 -13    50.27  1120.17   -2.719   -9.666
      5  6   6 -25    58.59  3807.40   -5.045   -9.048
      6  6   8 -26    63.28  4038.97   -4.578   -6.517
      7  6  -1 -21    87.74  2972.71   -0.436   -6.212
      8  6   1 -22    92.42  3204.28   -1.569   -7.980

whereh,k are the laticce indexes,X,Y are the positions of the unit cell center (correlation peaks, as suplied by MRC quadserch program) and[[IncrementoX]],[[IncrementoY]] are the unit cell deviation from the regular position (as computed by MRC quadserch program).

This output file is the input for the phantom generation program project

--Main.AlfredoSolano - 30 Jan 2007

Exp_Shift_Computation