Fitting the tilt axis model » History » Version 3
Anchi Cheng, 07/15/2010 10:25 PM
1 | 1 | Anchi Cheng | h1. Fitting the tilt axis model |
---|---|---|---|
2 | |||
3 | 2 | Anchi Cheng | h2. Find the Rough Tilt-Axis Location Manually |
4 | 1 | Anchi Cheng | |
5 | 2 | Anchi Cheng | The easiest way to find the tilt axis uses the fact that the projected feature moves in a cosine function when it is off from CCD or optical axis. |
6 | 1 | Anchi Cheng | |
7 | 2 | Anchi Cheng | Do this at the microscope: |
8 | 1 | Anchi Cheng | |
9 | 2 | Anchi Cheng | # Use stage alpha wobbler (symmetrical tilts) to put the specimen to eucentric height. |
10 | # Send to scope alternatively between 0 and 30 degree. Track the position of an object near the center of the viewing screen at the magnification you want use the model. If the tilt axis does not pass the center of the screen, the object will move upon this asymmetric movement. |
||
11 | # Move image shift x,y until the asymmetric tilt movement is minimized. |
||
12 | # Use a larger tilt to do the same if desired. |
||
13 | # Check the asymmetric movement between 0 and -30 degree to make sure it works reasonably in both way. |
||
14 | # Save the image shift as tomo preset's image shift. |
||
15 | # Estimate the combined image shift sqrt(x^2+y^2) and enter that as the custom, tilt axis-offset. |
||
16 | # From what is observed during the tilts, estimate the direction of the tilt axis and entered as the custom phi value. |
||
17 | # Perform test run with the custom *Fixed Tilt Axis Model" |
||
18 | |||
19 | h2. Checking and Improving Tilt-Axis Model |
||
20 | |||
21 | Test runs should be done at the magnification, defocus, and tilts up to the highest values planned for the real ones with a specimen with isotropic features in view and |
||
22 | 1 | Anchi Cheng | adjusted to eucentric height. Tilt steps should be made to allow at least 15 images taken |
23 | per tilt direction. |
||
24 | |||
25 | |||
26 | |||
27 | h3. First Test Run: Fixed Model |
||
28 | |||
29 | |||
30 | |||
31 | The first time the tomography application is used on a microscope, all model |
||
32 | parameters are default to zero. A test run should be done with the specimen at eucentric |
||
33 | height and all default settings including "keep the tilt axis parameters fixed" activated. |
||
34 | If the feature tracking is good, and the defocii of the images in the tilt series do not |
||
35 | significantly changed, the goniometer behavior and camera alignment are close to ideal, |
||
36 | and no improvement of the model is required. |
||
37 | |||
38 | |||
39 | If the xy tracking deteriorates quickly in the first few tilts, the magnification |
||
40 | should be dropped so that tracking is possible through out the tilt series. Once a rough |
||
41 | model is established, refinement can then be done at higher magnifications. |
||
42 | |||
43 | |||
44 | |||
45 | |||
46 | |||
47 | h3. Runs to Improve the Model-Fitting |
||
48 | |||
49 | |||
50 | |||
51 | A complete tilt model fitting is difficult with small number of data points or at only |
||
52 | low tilts. Therefore, two tilt series need to be taken in order to get a good averaged |
||
53 | model that is later fixed. |
||
54 | |||
55 | |||
56 | # tomography/settings/model>deactivate "keep the tilt axis parameters |
||
57 | fixed". |
||
58 | # tomography/settings/model>initialize the model with "this preset only" unless |
||
59 | this is a refinement of <link linkend="Tomo_low_mag_fitting">a rough model obtained at |
||
60 | a lower magnification</link>, in which case, initialize the model with "this preset |
||
61 | and lower mags". |
||
62 | # tomography/toolbar>click on "reset learning" !http://emg.nysbc.org/software/leginon/images/icons/refresh.png! to have |
||
63 | a fresh start of tilt series data included in the fitting. |
||
64 | # Presets Manager>send the "tomo" preset to scope. |
||
65 | # tomography>Start the tilt series image collection by clicking on "Simulated |
||
66 | Target" tool. |
||
67 | # repeat 4 and 5 to dynamically fit the model the second time if the first one goes |
||
68 | to completion. |
||
69 | If the target shifts away during the run, reduce the magnification, repeat 1 |
||
70 | through 9 in this procedure and then return to the higher, intended mag and refined |
||
71 | the model using the same procedure and proper setting for refinement as mentioned in |
||
72 | 2. |
||
73 | # Check on web image viewer to see if the result of the second run in step 6 |
||
74 | produces images of constant defocus through out each half of the tilt seriese. |
||
75 | # If the model is good, activate "keep the tilt axis parameters fixed" in the |
||
76 | tomography settings. |
||
77 | # Repeat step 2. Other than reset the learning, the tool produces an output of image |
||
78 | shift you should copy down and apply to the "tomo" preset in the future runs. |
||
79 | |||
80 | |||
81 | |||
82 | |||
83 | |||
84 | |||
85 | |||
86 | |||
87 | |||
88 | h2. Understand the results of the model fitting |
||
89 | |||
90 | |||
91 | |||
92 | If the above procedure does not yield a good model, it is necessary to study the graphs |
||
93 | output by the run to determine possible causes and derive the best fixed model that at least |
||
94 | work for most angles. |
||
95 | |||
96 | See [[Fitting the tilt axis model]] for definition of the displayed values |
||
97 | |||
98 | h3. First fitting tilt series |
||
99 | |||
100 | Figure 1 shows a typical result from the first dynamically fitted tilt series run. x and |
||
101 | |||
102 | *Figure 1* |
||
103 | !http://emg.nysbc.org/software/leginon/images/images/tomoxyz1.png! |
||
104 | |||
105 | According to the Zheng et. al. model, the tilt-axis can be characterized by three |
||
106 | parameters: phi, the angle between the tilt axis and the detector x axis (column), offset, |
||
107 | the distance between the center of the detector and the tilt axis on the xy plane, and z0, |
||
108 | the distance between the specimen from the xy plane that contains the tilt axis. |
||
109 | |||
110 | As a first tilt series used in the model fitting, the data is not fitted until it |
||
111 | accumulate enough data points and spread over more tilts. We arbitrarily select 30 degrees as |
||
112 | the starting angle where the fit is started. |
||
113 | |||
114 | As can be seen in Figure 1, the model is significantly different from the initial as z-axis "Prediction" jumps once the fitting is started. In particular, the model z0 which corresponds to the offset of the specimen from the tilt axis in z direction dominates deviation. Such problem is evident by the near-straight line feature shift in x-axis over more than 1 um. There is also some offset in |
||
115 | x direction since the "Prediction" and "Position" flattens somewhat as it approaches zero tilt. The tilt axis has minimal tilt from x-axis since in the y-direction, despite the bump at about 40 degrees, the range of variation in the whole tilt series is about 0.06 um. |
||
116 | |||
117 | By activate "show model parameters" check box in the viewer, more plots are shown |
||
118 | (Figure 2). phi and "optical axis" offset remain fixed below 30 deg while z0 is recalculated at each tilt. As you can see, the fitting of the former two parameters are not very stable in this first trial. Therefore, we need to do the second run. |
||
119 | |||
120 | *Figure 2* |
||
121 | !http://emg.nysbc.org/software/leginon/images/images/tomomodel1.png! |
||
122 | |||
123 | h3. Second fitting tilt series |
||
124 | |||
125 | Figure 3 shows the z-axis plot of the second run since the x and y does not change |
||
126 | significantly. Figure 4 shows its corresponding model parameters. The zero and first tilt |
||
127 | uses the initial model and the rest are fitted dynamically. In this run, because more than |
||
128 | one tilt series is found in the memory since last "reset learning", the model fitting starts |
||
129 | after the first tilt image is acquired. model parameters are rather stable over the whole |
||
130 | tilt series and suggests that the tilt axis is tilted from x-axis by about -1.3 degree, and |
||
131 | offset from the center of the detector by -0.6 um, and the specimen is off from eucentric |
||
132 | height by about -0.7 um |
||
133 | |||
134 | *Figure 3* |
||
135 | !http://emg.nysbc.org/software/leginon/images/images/tomoz2.png! |
||
136 | |||
137 | *Figure 4* |
||
138 | !http://emg.nysbc.org/software/leginon/images/images/tomomodel2.png! |
||
139 | |||
140 | |||
141 | Although it is not possible to maintain perfect defocus prediction in the full range of |
||
142 | the tilt by using a fixed model, we found that the overall performance is better if an |
||
143 | average tilt axis model that works well in the mid-range tilt is used as a fixed model. For |
||
144 | example, for the behavior in the above figure, we choose: |
||
145 | |||
146 | |||
147 | h3. Fixed axis tilt series |
||
148 | |||
149 | Use these custom values as initial model and turn on "Keep the tilt axis parameters |
||
150 | fixed". If the result is good, judging by consistent defocus and target tracking through out |
||
151 | the tilt series, this model will be saved in the database as best model automatically, and |
||
152 | you can revert back to initialize with the model of "only this preset". |
||
153 | 3 | Anchi Cheng | |
154 | h2. Low Magnification Model Fitting |
||
155 | |||
156 | The fitting of optical axis offset does not always works if the offset is so large that |
||
157 | the feature moves out of view with even a small tilt. In such a case, it is worth first |
||
158 | collect a tomography series at a lower magnification to define roughly the model. |
||
159 | |||
160 | |||
161 | At beginning of each session, or forced by the user, the model is initialized. By |
||
162 | default, at the initialization, Tomography node uses past fitting results that show good |
||
163 | agreement with the experimental data at the magnification of the preset used. If a good |
||
164 | model is not found, that from lower magnifications will be used. It is possible to force the |
||
165 | node to use a model fitted at a particular magnification by selecting it in |
||
166 | Tomography/Settings/Model. |
||
167 | |||
168 | |||
169 | Therefore, we recommend that, in case of fitting failure on good contrast images, the |
||
170 | followings should be done: |
||
171 | |||
172 | |||
173 | # Tomography/Settings/Image Acquisition> change the preset to "hl". |
||
174 | # Tomography/Settings> adjust Tilt and Exposure parameters to match. |
||
175 | # Acquire the tiltseries images. |
||
176 | # If the tracking is good, change the preset back. |
||
177 | # Tomography/Settings/Model> Initialize with the model of (the mag of "hl" |
||
178 | preset). |
||
179 | # Acquire the tomo-series. |
||
180 | # If tracking is good, change back to Initialize with the model of "this preset and |
||
181 | lower" mag. |
||
182 | |||
183 | See [[Troubles in Tomography]] section for various modeling problem solving examples |