Project

General

Profile

Fitting the tilt axis model » History » Version 6

Anchi Cheng, 07/15/2010 10:52 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 5 Anchi Cheng
The easiest way to find the tilt axis uses the fact that the projected feature moves according to a cosine function when it is off from CCD or optical axis but not off from eucentric height.
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 6 Anchi Cheng
# Perform the quick test run with this custom *Fixed Tilt Axis Model"
18 2 Anchi Cheng
19 6 Anchi Cheng
h2. Quick fixed model test run
20 1 Anchi Cheng
21 6 Anchi Cheng
Instead of a full range tilt series, much can be learned at low tilts.  Again, this uses the fact that defocus changes in a sine curve when the tilt axis is off.
22 2 Anchi Cheng
23 6 Anchi Cheng
Use Leginon Tomography to do a small range tilt since sine curve has the highest slope around zero.
24
# Set the tilt series to cover max 8, min 8 and step at 2 deg.
25
# Take tilt series at these low tilts and at an area where defocus change can be observed easily.  Increase dose if needed.
26
# Compare the defocus of the resulting images.  Optimize the tilt axis offset so that defocus differences among the tilts in the same tilt direction is minimized.
27 1 Anchi Cheng
28 6 Anchi Cheng
h2. Full Length  Tilt-Axis Model Testing and Fitting
29 1 Anchi Cheng
30 6 Anchi Cheng
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
31
adjusted to eucentric height. Tilt steps should be made to allow at least 15 images taken per tilt direction.  Note that this can refine the above rough results but may give
32
worse results at times.
33 1 Anchi Cheng
34 6 Anchi Cheng
h3. First Test Run: Fixed Model
35 1 Anchi Cheng
36
The first time the tomography application is used on a microscope, all model
37
parameters are default to zero. A test run should be done with the specimen at eucentric
38
height and all default settings including "keep the tilt axis parameters fixed" activated.
39
If the feature tracking is good, and the defocii of the images in the tilt series do not
40
significantly changed, the goniometer behavior and camera alignment are close to ideal,
41
and no improvement of the model is required.
42
43
44
If the xy tracking deteriorates quickly in the first few tilts, the magnification
45
should be dropped so that tracking is possible through out the tilt series. Once a rough
46
model is established, refinement can then be done at higher magnifications.
47
48
49
50
51
52
h3. Runs to Improve the Model-Fitting
53
54
55
56
A complete tilt model fitting is difficult with small number of data points or at only
57
low tilts. Therefore, two tilt series need to be taken in order to get a good averaged
58
model that is later fixed.
59
60
61
# tomography/settings/model>deactivate "keep the tilt axis parameters
62
fixed".
63
# tomography/settings/model>initialize the model with "this preset only" unless
64
this is a refinement of <link linkend="Tomo_low_mag_fitting">a rough model obtained at
65
a lower magnification</link>, in which case, initialize the model with "this preset
66
and lower mags".
67
# tomography/toolbar>click on "reset learning" !http://emg.nysbc.org/software/leginon/images/icons/refresh.png! to have
68
a fresh start of tilt series data included in the fitting.
69
# Presets Manager>send the "tomo" preset to scope.
70
# tomography>Start the tilt series image collection by clicking on "Simulated
71
Target" tool.
72
# repeat 4 and 5 to dynamically fit the model the second time if the first one goes
73
to completion.
74
 If the target shifts away during the run, reduce the magnification, repeat 1
75
through 9 in this procedure and then return to the higher, intended mag and refined
76
the model using the same procedure and proper setting for refinement as mentioned in
77
2.
78
# Check on web image viewer to see if the result of the second run in step 6
79
produces images of constant defocus through out each half of the tilt seriese.
80
# If the model is good, activate "keep the tilt axis parameters fixed" in the
81
tomography settings.
82
# Repeat step 2. Other than reset the learning, the tool produces an output of image
83
shift you should copy down and apply to the "tomo" preset in the future runs.
84
85
86
87
88
89
90
91
92
93
h2. Understand the results of the model fitting
94
95
96
97
If the above procedure does not yield a good model, it is necessary to study the graphs
98
output by the run to determine possible causes and derive the best fixed model that at least
99
work for most angles.
100
101
See [[Fitting the tilt axis model]] for definition of the displayed values
102
103
h3. First fitting tilt series
104
105
Figure 1 shows a typical result from the first dynamically fitted tilt series run. x and
106
107
*Figure 1* 
108
!http://emg.nysbc.org/software/leginon/images/images/tomoxyz1.png!
109
110
According to the Zheng et. al. model, the tilt-axis can be characterized by three
111
parameters: phi, the angle between the tilt axis and the detector x axis (column), offset,
112
the distance between the center of the detector and the tilt axis on the xy plane, and z0,
113
the distance between the specimen from the xy plane that contains the tilt axis.
114
115
As a first tilt series used in the model fitting, the data is not fitted until it
116
accumulate enough data points and spread over more tilts. We arbitrarily select 30 degrees as
117
the starting angle where the fit is started.
118
119
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
120
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.
121
122
By activate "show model parameters" check box in the viewer, more plots are shown
123
(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.
124
125
*Figure 2* 
126
!http://emg.nysbc.org/software/leginon/images/images/tomomodel1.png!
127
128
h3. Second fitting tilt series
129
130
Figure 3 shows the z-axis plot of the second run since the x and y does not change
131
significantly. Figure 4 shows its corresponding model parameters. The zero and first tilt
132
uses the initial model and the rest are fitted dynamically. In this run, because more than
133
one tilt series is found in the memory since last "reset learning", the model fitting starts
134
after the first tilt image is acquired. model parameters are rather stable over the whole
135
tilt series and suggests that the tilt axis is tilted from x-axis by about -1.3 degree, and
136
offset from the center of the detector by -0.6 um, and the specimen is off from eucentric
137
height by about -0.7 um
138
139
*Figure 3* 
140
!http://emg.nysbc.org/software/leginon/images/images/tomoz2.png!
141
142
*Figure 4*
143
!http://emg.nysbc.org/software/leginon/images/images/tomomodel2.png!
144
145
146
Although it is not possible to maintain perfect defocus prediction in the full range of
147
the tilt by using a fixed model, we found that the overall performance is better if an
148
average tilt axis model that works well in the mid-range tilt is used as a fixed model. For
149
example, for the behavior in the above figure, we choose:
150
151
152
h3. Fixed axis tilt series
153
154
Use these custom values as initial model and turn on "Keep the tilt axis parameters
155
fixed". If the result is good, judging by consistent defocus and target tracking through out
156
the tilt series, this model will be saved in the database as best model automatically, and
157
you can revert back to initialize with the model of "only this preset".
158 3 Anchi Cheng
159
h2. Low Magnification Model Fitting
160
161
The fitting of optical axis offset does not always works if the offset is so large that
162
the feature moves out of view with even a small tilt. In such a case, it is worth first
163
collect a tomography series at a lower magnification to define roughly the model.
164
165
166
At beginning of each session, or forced by the user, the model is initialized. By
167
default, at the initialization, Tomography node uses past fitting results that show good
168
agreement with the experimental data at the magnification of the preset used. If a good
169
model is not found, that from lower magnifications will be used. It is possible to force the
170
node to use a model fitted at a particular magnification by selecting it in
171
Tomography/Settings/Model.
172
173
174
Therefore, we recommend that, in case of fitting failure on good contrast images, the
175
followings should be done:
176
177
178
#  Tomography/Settings/Image Acquisition> change the preset to "hl".
179
#  Tomography/Settings> adjust Tilt and Exposure parameters to match.
180
#  Acquire the tiltseries images.
181
#  If the tracking is good, change the preset back.
182
#  Tomography/Settings/Model> Initialize with the model of (the mag of "hl"
183
preset).
184
#  Acquire the tomo-series.
185
#  If tracking is good, change back to Initialize with the model of "this preset and
186
lower" mag.
187
188 1 Anchi Cheng
See [[Troubles with Tomography]] section for various modeling problem solving examples
189 5 Anchi Cheng
190
______
191
192
[[Reading the tomography graphs|< Reading the tomography graphs]] | [[Full Protocol on a F30 with an energy filter| Full Protocol on a F30 with an energy filter >]]