Difference between revisions of "Cross correlation matrix"
Jump to navigation
Jump to search
| Line 7: | Line 7: | ||
This similarity of particles ''i'' and ''j'' measured in terms of the normalized cross correlation of the the two aligned particles, filtered to their common fourier components, and restricted to a region in direct space (indicated by a classification mask). The pseudo code will run as: | This similarity of particles ''i'' and ''j'' measured in terms of the normalized cross correlation of the the two aligned particles, filtered to their common fourier components, and restricted to a region in direct space (indicated by a classification mask). The pseudo code will run as: | ||
<code> | <code> | ||
| − | # read | + | # read particles ''i'' and ''j'' -> p<sub>i</sub>, p<sub>j</sub> |
| − | # | + | # align particles ''i'', and ''j'' -> A<sub>i</sub>p<sub>i</sub>, A<sub>j</sub>p<sub>j</sub> |
| − | # | + | # compute missing wedges for particles ''i'' and ''j'' -> W<sub>i</sub>, W<sub>j</sub> |
| − | + | # rotate missing wedges -> R<sub>i</sub>W<sub>i</sub>, R<sub>j</sub>W<sub>j</sub> | |
| − | # compute Fourier coefficients | + | # compute Fourier coefficients common to bCoth rotated missing wedges |
| + | C<sub>ij</sub> = intersection(R<sub>i</sub>W<sub>i</sub>, R<sub>j</sub>W<sub>j</sub>) | ||
| + | # filter p<sub>i</sub> with C<sub>ij</sub> | ||
| + | # compare the filtered particles inside the classification mask. | ||
</code> | </code> | ||
Revision as of 10:42, 19 April 2016
The cross correlation matrix (often called ccmatrix in Dynamo jargon) of a set of N particles is an N X N matrix. Each entry (i,j) represents the similarity of particles i and j in the data set.
Contents
Definition of similarity
This similarity of particles i and j measured in terms of the normalized cross correlation of the the two aligned particles, filtered to their common fourier components, and restricted to a region in direct space (indicated by a classification mask). The pseudo code will run as:
- read particles i and j -> pi, pj
- align particles i, and j -> Aipi, Ajpj
- compute missing wedges for particles i and j -> Wi, Wj
- rotate missing wedges -> RiWi, RjWj
- compute Fourier coefficients common to bCoth rotated missing wedges
Cij = intersection(RiWi, RjWj)
- filter pi with Cij
- compare the filtered particles inside the classification mask.
