Difference between revisions of "Walkthrough on GUI based tilt series alignment"

From Dynamo
Jump to navigation Jump to search
 
(87 intermediate revisions by 2 users not shown)
Line 1: Line 1:
Dynamo includes a package for automated alignment and reconstruction of tilt series. This walkthrough guides you through the steps on how to use it in GUI-based mode. A different [[Walkthrough on command line based tilt series alignment | page in this wiki]] shows how to operate this procedure through the command line.
+
''Dynamo'' includes a package for automated alignment and reconstruction of tilt series. This walkthrough guides users through a series of steps on how to use it in GUI-based mode. The [[Walkthrough on command line based tilt series alignment | walkthrough on command line based tilt series alignment]] shows how to operate this procedure through the command line.
  
 
= Data =
 
= Data =
  
We have prepared a version of a tilt series from the EMPIAR entry 10164, containing a set of virus like particles of VLP.
+
We use a binned version of the first tilt series from the EMPIAR entry 10164, depicting a set of virus like particles (VLP). Visit the course materials to access the tilt series or download it directly from EMPIAR. You can also follow this tutorial using your own data.
 
<tt>wget https://wiki.dynamo.biozentrum.unibas.ch/html/w/doc/data/hiv/align/b001ts001.mrc</tt>
 
 
 
In a Mac:
 
 
 
<tt>curl -O https://wiki.dynamo.biozentrum.unibas.ch/html/w/doc/data/hiv/align/b001ts001.mrc</tt>
 
 
 
In the workshop, a copy should be already available on your local data folder:
 
 
 
<tt>~/data/b001ts001.mrc</tt>
 
 
 
In this binned version, the pixel size is 2.7 Angstrom (original was 1.35, this is bin level 1 in ''Dynamo'' or bin level 2 in IMOD)
 
  
 
= Creating an alignment workflow=  
 
= Creating an alignment workflow=  
  
The basic function for invoking the GUI is <tt>dtsa</tt> (short form of <tt>dynamo_tilt_series_alignment</tt>:  
+
Start the alignment GUI by typing the command <tt>dtsa</tt> (short form of <tt>dynamo_tilt_series_alignment</tt>) into the ''Dynamo'' command line:  
  
   <tt>u = dtsa();</tt>
+
   <tt>dtsa</tt>
  
Here, <tt>u</tt> is an object that will reside in memory during the session and allows to interoperate with the workflow through the command line. It is not necessary when proceeding through the GUI.
+
A pop-up window appears. Click on "create a new workflow":
After running the command above, a pop-up window appears. Click on create a new workflow.
 
  
 
[[ File:alignGUI_creatingAWorkflow.png |thumb|center| 600px|creating a workflow]]
 
[[ File:alignGUI_creatingAWorkflow.png |thumb|center| 600px|creating a workflow]]
  
The data for the workflow (tilt series, angles) could have been introduced directly when creating the workflow, we however will introduce each item explicitly in this walkthrough.  
+
Here, the input data for the workflow (tilt series, angles, ...) is introduced:
 
 
The minimal data to initiate a workflow is
 
* tilt series matrix (use the file indicated above)
 
* tilt angles.
 
* discarded tilt angles (leave empty)
 
  
 
[[ File:alignGUI_inputWorkflowParameters.png |thumb|center| 600px|input workflow parameters]]
 
[[ File:alignGUI_inputWorkflowParameters.png |thumb|center| 600px|input workflow parameters]]
  
== Discarded tilt angles ==
+
The option <tt>configuration file</tt> is used to pass information on possible tilts that should be discarded. Here, we leave this option untouched. Press on ''create workflow'', to open the alignment GUI.
 
 
If you are unsure on the tilts to be discarded, you can do it in a later stage while you inspect the images in the  <tt>dmarkers</tt> GUI, using the <tt>x</tt> and <tt>shift + x</tt> buttons.  
 
 
 
  
 
= GUI description =
 
= GUI description =
 
+
The alignment GUI guides the process of generating gold bead coordinates, aligning the tilt series, correcting the CTF and to finally reconstruct the tomogram. It has 5 parts. Starting from top to bottom, they are:
The GUI guides the process of creating and editing a  set of '''working markers''' (in the areas of '''detection''', '''reindexing''' and '''refinement'''), and then using them to create aligned stacks ('''alignment''' tab), analyse and correct their CTF (through wrappers to <tt>ctfffind4</tt> and <tt>imod</tt>) to create 3D reconstructions of the tomographic data.  
+
# '''Menu:''' Different tools for visualization and customization of general settings. The option ''Settings'' > ''edit[general]'' for example lets you adjust the pixel size and different acquisition parameters such as Cs etc. Here we keep the default settings.
 +
# '''Icons:''' Quick access to tools such as ''dmarkers'', ''occupancy plot'', ''gold bead gallery'', ''gold bead targeting'', ''fitting'' and ''marker selection''. Their functionalities will be explained later.
 +
# '''Tabs:''' The tabs describe the 6 ''execution areas'' that are executed from left to right: ''detection'', ''indexing'', ''refinement'', ''alignment'', ''ctf'', ''reconstruction''. Each execution areas contains a set of processing steps.
 +
# '''Steps:''' The processing steps of each tab are run sequentually from top to bottom. However, after execution, the steps can be revisited to change or update information. Steps need to be active (letters highlighted in black, not grey) to be executed. They are automatically activated once the previous steps have been executed or the required input is given. Input parameters are edited by clicking the plus button (+) on the left of the step button. When a step is under execution, the text of the button becomes green. After execution the text becomes bold. At the right side of each of the steps, a tool icon is located. The tool icon can be right-clicked to provide a pop-up menu with specific tools for each step. Next to the tool icon there is the ''View'' button which provides specific tools to visualize and inspect the results of each step.
 +
# '''Run options:''' Options to run multiple steps.
  
 
[[ File:alignGUI_alignmentGUIBeforeCompletionOfAnyStep.png |thumb|center| 600px|alignment GUI before completion of any step]]
 
[[ File:alignGUI_alignmentGUIBeforeCompletionOfAnyStep.png |thumb|center| 600px|alignment GUI before completion of any step]]
== Settings==
 
 
=== Acquisition settings ===
 
 
Here you can enter pixel size (2.7 Angstroms), Cs (2.2), nominal defocus (-4.0 microns, used if the CTF area is connected).
 
 
== Alignment GUI Controls==
 
 
The GUI possesses different pre-defined areas. Starting from top to bottom, you can find a toolbox with the options: ''Tools'', ''Help'', ''Settings'', ''IO'', ''Working Markers'' and ''Fit''. These tools are extremely useful to visualize and customize general settings as tilt scheme and defocus value among others. Each of the tools are briefly explained:
 
 
* Tools
 
* Help
 
* Settings
 
* IO
 
* Working markers
 
* Fit
 
 
Under the toolbox, we find several icons: ''dmarkers'', ''occupancy plot'', ''gold bead gallery'', ''gold bead targeting'', ''fitting'' and ''marker selection''.
 
 
* dmarkers GUI: this icon opens a GUI that displays the tilt series data and the current working markers. It allows you to check the position of the working markers in each micrograph. If you want to make editions of the GUI permanent, you need to save the current markers as explained bellow.
 
 
* Occupancy plot: shows a scheme of which markers are present in which tilt numbers.
 
 
* Gold bead gallery: it crops the 2D patches which are currently contained in the working markers, and therefore, treated as gold beads, and display them in gallery mode.
 
 
= Execution Areas =
 
In the GUI, right below the icons, we can find six ''execution areas'': '''detection''', '''indexing''', '''refinement''', '''alignment''', '''ctf''', '''reconstruction'''. Each of them comprises a sets of steps that can be executed sequentially. Sequential execution is the natural way of proceeding; but steps can be visited in any order, i.e., a step can be re-executed in a later moment after updating its information.
 
 
[[ File:alignGUI_passingTheCursorOnTheGUIElementsShowsAHelpElement.png |thumb|center| 600px| By passing the cursor on the GUI elements shows a help element  ]]
 
 
=== Steps ===
 
Each ''execution area'' has one or more steps. Steps need to be ''open'' (letters highlighted in black, not grey) to be executed. They become open when the required input(s) are ready to be processed. The required input(s) can be provided by users or computed and stored in a previous step of this alignment workflow. When a step is under execution, its pushbutton becomes green.
 
 
* Step Toolbox
 
At the right side of each of the steps, a tool icon is located. The tool icon can be secondary clicked to provide a pop-up menu with specific tools for each step. There are also some common tools for all the steps which handle the inspection and listing of input and output items generated by the step.
 
 
Next to the tool icon, there is a <tt>View</tt> pushbutton. The <tt>View</tt> pushbutton provides specific tools to visualize and inspect the results of each step.
 
  
 +
= Processing =
 +
We are now ready to run all the processing steps. In the following, each step is described in detail.
 
== Detection ==
 
== Detection ==
 
The goal of this step is to find the positions of the gold bead projections in the tilt series.  To do that, ''Dynamo'' cross correlates each micrograph against a 2D model of the gold bead and then locate the cross correlation peaks, which provide putative positions for actual projections of gold beads. Further analysis  of these peaks will determine which ones correspond to actual gold bead projection (''peak feature'' and ''peak selection'').
 
The goal of this step is to find the positions of the gold bead projections in the tilt series.  To do that, ''Dynamo'' cross correlates each micrograph against a 2D model of the gold bead and then locate the cross correlation peaks, which provide putative positions for actual projections of gold beads. Further analysis  of these peaks will determine which ones correspond to actual gold bead projection (''peak feature'' and ''peak selection'').
Line 92: Line 40:
 
Click on ''creates a binned tilt series'' button to create a binned tilt series stack.
 
Click on ''creates a binned tilt series'' button to create a binned tilt series stack.
 
Rightclick on ''View'' button to visualize the newly binned tilt series stack in <tt>tmshow</tt>. In the <tt>tmshow</tt> window, visualization parameters as binning, band and thickness can be adjusted if needed.
 
Rightclick on ''View'' button to visualize the newly binned tilt series stack in <tt>tmshow</tt>. In the <tt>tmshow</tt> window, visualization parameters as binning, band and thickness can be adjusted if needed.
[[ File:alignGUI_completedStepsBecomeBoldFontSize.png |thumb|center| 600px|completed steps become bold font size]]
 
  
Now, we are going to measure the diameter of the gold beads precisely in <tt>tmshow</tt>. Zoom in a gold bead by holding the <tt>shift</tt> key and scrolling up your mouse wheel. Adjust the visualized region by holding the <tt>shift</tt> key again and dragging with right mouse click on <tt>tmshow</tt> window.
+
Now, we are going to measure the diameter of the gold beads precisely in <tt>tmshow</tt> in the original (unbinned) tilt series data. To do that, press the ''IO'' in the upper toolbox of the ''alignment workflow GUI'', then go to ''tilt series matrix > inspect > On-the-fly view [dtmshow]''. The <tt>tmshow</tt> window opens, adjust contrast if necessary for better visualization.
 +
[[ File:openOriginalTiltSeriesDataGUI.png |thumb|center| 600px| open unbinned tilt series data to measure gold bead diameter]]
 +
 
 +
Zoom in a gold bead by holding the <tt>shift</tt> key and scrolling up your mouse wheel (alternatively, use the icon with the red rectangle. Avoid using the lens icons). Adjust the visualized region by holding the <tt>shift</tt> key again and dragging with right mouse click on <tt>tmshow</tt> window.
  
 
[[ File:alignGUI_openTomoshowOnTheFullSizedMatrixToCheckTheSizeOfAGoldBead.png |thumb|center| 600px|open tomoshow on the full sized matrix to check the size of a gold bead]]
 
[[ File:alignGUI_openTomoshowOnTheFullSizedMatrixToCheckTheSizeOfAGoldBead.png |thumb|center| 600px|open tomoshow on the full sized matrix to check the size of a gold bead]]
Line 100: Line 50:
 
Once you have located the gold bead of interest, click on the ruler icon at the top toolbox to activate length measurements on the displayed image data. Draw a line across the gold bead, and the length of the line appears in pixels, which correspond to the gold bead diameter. In this case, the length of this line should be 32 pixel aprox.
 
Once you have located the gold bead of interest, click on the ruler icon at the top toolbox to activate length measurements on the displayed image data. Draw a line across the gold bead, and the length of the line appears in pixels, which correspond to the gold bead diameter. In this case, the length of this line should be 32 pixel aprox.
  
[[ File:alignGUI_theRulerToolAllowsToMeasureInPixelsTheDiameterOfAGoldBead.png |thumb|center| 600px|the ruler tool allows to measure in pixels the diameter of a gold bead]]
+
[[ File:OriginalGoldBeadDiameter.png |thumb|center| 600px|the ruler tool allows to measure in pixels the diameter of a gold bead]]
  
 
=== Detection of gold beads ===
 
=== Detection of gold beads ===
This area is the one that needed some design decision on the size of the gold bead. We will use <tt>dtmshow</tt>
+
Here we provide the size of the gold beads. The measured diameter of the gold beads is 32 pixels in the original image, therefore its radius is 16 pixels. We introduce the bead radius and some other parameters (see screenshot) to create an accurate initial gold bead template (here, the binning is just used internally for quick and robust bead detection).
 
 
[[ File:alignGUI_thePlusIconAtCornerShowsTheValuesOfTheParametersOfAGivenStep.png |thumb|center| 600px|the plus icon at corner shows the values of the parameters of a given step]]
 
  
 +
[[ File:alignGUI_DetectionParameters.png |thumb|center| 600px|the plus icon at corner reveals the values of the parameters of a given step.]]
  
 +
After pressing the button ''detect peaks'', ''Dynamo'' creates a raw cross-correlation matrix based on the initial template which is binary and extracts 300 cc-peaks by default. Then, the extracted peaks are averaged to create a refined template. The refined template is used to compute a refined cc-matrix from which 300 peaks are extracted and used as initial working markers. Once the ''detect peaks'' button turns black bold, the computation is finished. Right click on the <tt>view</tt> area of the step to access the different visualisation options.
  
Right click on the <tt>view</tt> area of the step to access the different visualisation options.
+
[[ File:alignGUI_markerSetsCanBeViewedAsCloudOfPoints.png |thumb|center| 600px| Secondary click the view panel to open a menu to access different results of a step. Marker sets can be viewed as cloud of points.]]
[[ File:alignGUI_secondaryClickInTheViewPanelOpensAMenuToAccessTheDifferentResultsOfAStep.png |thumb|center| 600px|secondary click in the view panel opens a menu to access the different results of a step ]]
 
[[ File:alignGUI_markerSetsCanBeViewedAsCloudOfPoints.png |thumb|center| 600px|marker sets can be viewed as cloud of points]]
 
  
This view shows the x,y coordinates of ''all'' cross correlation peaks on all the micrographs in the tilt series.
+
This view shows the overlaid sum of the x,y coordinates of ''all'' cc peaks from ''all'' micrographs in the tilt series.
 
[[ File:alignGUI_cloudOfUnindexedCcPeaks.png |thumb|center| 600px|cloud of unindexed cc peaks]]
 
[[ File:alignGUI_cloudOfUnindexedCcPeaks.png |thumb|center| 600px|cloud of unindexed cc peaks]]
 
We can also check the positions of all the detected spots on each micrograph:
 
We can also check the positions of all the detected spots on each micrograph:
 
[[ File:alignGUI_showsTheDetectedSpotsOnTheVideoGUI.png |thumb|center| 600px|shows the detected spots on the video GUI]]
 
[[ File:alignGUI_showsTheDetectedSpotsOnTheVideoGUI.png |thumb|center| 600px|shows the detected spots on the video GUI]]
This GUI does not allow for edition of the markers.  
+
This GUI does not allow for edition of the markers. '''Tip:''' enable ''auto'' option at the bottom center of the visualization GUI to facilitate tilt series visualization.
 
[[ File:alignGUI_staticVideoAnnotationsCannotBeEdited.png |thumb|center| 600px|static video annotations cannot be edited]]
 
[[ File:alignGUI_staticVideoAnnotationsCannotBeEdited.png |thumb|center| 600px|static video annotations cannot be edited]]
The neighbourhood of cross correlation peaks can be visualised individually:
+
The neighbourhood of cross correlation peaks can be visualised individually as 2D patches from the tilt series image data:
 
[[ File:alignGUI_EachPatchWillBeAPutativeGoldBead.png |thumb|center| 600px|Each patch will be a putative gold bead]]
 
[[ File:alignGUI_EachPatchWillBeAPutativeGoldBead.png |thumb|center| 600px|Each patch will be a putative gold bead]]
 
[[ File:alignGUI_theSetOfCcPeaksIsShownAsAMontage.png |thumb|center| 600px|the set of cc peaks is shown as a montage]]
 
[[ File:alignGUI_theSetOfCcPeaksIsShownAsAMontage.png |thumb|center| 600px|the set of cc peaks is shown as a montage]]
 
This visualization is normally a good gauge on how the detection procedure has worked. The example below would be an example of detection procedure carried with wrong search parameters
 
This visualization is normally a good gauge on how the detection procedure has worked. The example below would be an example of detection procedure carried with wrong search parameters
 
[[ File:alignGUI_exampleOfFailedSearch.png |thumb|center| 600px|example of failed gold bead detection]]
 
[[ File:alignGUI_exampleOfFailedSearch.png |thumb|center| 600px|example of failed gold bead detection]]
Note that you have the option of testing your parameters against a single image before launching a full detection procedure.
+
Note that you have the option of testing your parameters against a single image before launching a full detection procedure. Right-click on the tool icon on the right side of ''detect peaks'' pushbutton.
 
[[ File:alignGUI_toolboxOptionForTestingTheParametersOnACentralImage.png |thumb|center| 600px|toolbox option for testing the parameters on a central image]]
 
[[ File:alignGUI_toolboxOptionForTestingTheParametersOnACentralImage.png |thumb|center| 600px|toolbox option for testing the parameters on a central image]]
 
[[ File:alignGUI_resultsOfTestingOnACentralImage.png |thumb|center| 600px|results of testing on a central image]]
 
[[ File:alignGUI_resultsOfTestingOnACentralImage.png |thumb|center| 600px|results of testing on a central image]]
  
=== Computation of observation features s===
+
=== Computation of observation features ===
At this stage, the observed cross correlation peaks are analysed. In the current ''Dynamo''' version a merit figure is assigned to each peak based on a "rotational merit" that measures its similarity to a circle. The results can be depicted on patches representing individual gold beads:
+
At this stage, the observed cross correlation peaks are analysed. A merit figure is assigned to each peak based on a ''rotational merit'' that measures its rotational symmetry. Press the button ''Peak features'' to compute the rotational symmetry of the current working markers. The working markers that do not meet rotational symmetry will be discarded. The results can be depicted on patches representing individual gold beads:
 
[[ File:alignGUI_patchGalleryForFeatures.png |thumb|center| 600px|patch gallery for features]]
 
[[ File:alignGUI_patchGalleryForFeatures.png |thumb|center| 600px|patch gallery for features]]
[[ File:alignGUI_theGalleryAllowsRelatingAppearanceOfGoldBeadAnValueOfTheFeatures.png |thumb|center| 600px|the gallery allows relating appearance of gold bead an value of the features]]
+
[[ File:alignGUI_theGalleryAllowsRelatingAppearanceOfGoldBeadAnValueOfTheFeatures.png |thumb|center| 600px|the gallery shows the detected gold beads, the CC score (yellow) and the rotational merit (red). Tilt number is purple.]]
 
or extracting properties of the whole data set.
 
or extracting properties of the whole data set.
[[ File:alignGUI_featurePlotterAllowsCheckingOfTheFeatureValuesOnTheWholeSet.png |thumb|center| 600px|feature plotter allows checking of the feature values on the whole set ]]
+
[[ File:alignGUI_featurePlotterAllowsCheckingOfTheFeatureValuesOnTheWholeSet.png |thumb|center| 600px|feature plotter allows checking of the feature values on the whole set]]
[[ File:alignGUI_differentPlottingOptions.png |thumb|center| 600px|different plotting options ]]
+
[[ File:alignGUI_scatter2dOption.png |thumb|center| 600px| We select ''scatter'', set the y axis to ''rotationaMerit'' and then press ''Scatter 2d.]]
[[ File:alignGUI_scatter2dOption.png |thumb|center| 600px|scatter 2d option]]
+
[[ File:alignGUI_eachPointRepresentsACcPeak.png |thumb|center| 600px| We get a scatter plot of each peak and its relation to the CC and rotational merit score. Seeing two distinct populations means that the majority of the gold beads are distinguishable from other signal such as background.]]
[[ File:alignGUI_coloringOptionForScatteredPlots.png |thumb|center| 600px|coloring option for scattered plots]]
 
[[ File:alignGUI_eachPointRepresentsACcPeak.png |thumb|center| 600px|each point represents a cc peak]]
 
  
 
=== Selection of best gold beads===
 
=== Selection of best gold beads===
The cluster with the best markers is then selected. The individual patches can be inspected
+
The population with the best markers is automatically selected by clicking on ''Peak selection''. The results can be inspected again in a plot
 
[[ File:alignGUI_checkObservationCloudAfterSelectionOrClustering.png |thumb|center| 600px|check observation cloud after selection or clustering ]]
 
[[ File:alignGUI_checkObservationCloudAfterSelectionOrClustering.png |thumb|center| 600px|check observation cloud after selection or clustering ]]
or, again the cloud of observations
 
 
[[ File:alignGUI_observationCloudHasBeenThinned.png |thumb|center| 600px|observation cloud has been thinned]]
 
[[ File:alignGUI_observationCloudHasBeenThinned.png |thumb|center| 600px|observation cloud has been thinned]]
and the location of the observations in the micrographs:
+
or in the micrographs themselves:
 
[[ File:alignGUI_checkMarkerPositionsAfterClustering.png |thumb|center| 600px|check marker positions after clustering]]
 
[[ File:alignGUI_checkMarkerPositionsAfterClustering.png |thumb|center| 600px|check marker positions after clustering]]
 
[[ File:alignGUI_markerPositionsAfterClustering.png |thumb|center| 600px|marker positions after clustering]]
 
[[ File:alignGUI_markerPositionsAfterClustering.png |thumb|center| 600px|marker positions after clustering]]
  
 
== Indexing ==  
 
== Indexing ==  
 
+
The goal of this execution area is to find ''trails'' of gold beads, i.e., we connect the observations from different micrographs that correspond to the same gold bead. All the steps in this area will constantly update a set of ''working markers''.
The goal of this task is to find ''trails'' of gold beads, i.e., observations in different micrographs that correspond to the same 3d gold bead. All the steps in this area will update a set of "working markers".
 
  
 
=== Rough Alignment ===
 
=== Rough Alignment ===
Here all couples of micrographs are compared. For each couple, we look for the relative shifts that generates the most matches of cross-correlation peaks between one micrograph and the next one. This procedure induces trails of matched pairs of observations of variable length along the tilt series. The scheme that depicts which trails have representatives in which tilts is called ''occupancy graph''
+
In this step, micrographs are compared pair-wise. For each couple, we look for the relative shifts that generate the most matches of cross-correlation peaks between one micrograph and the next one. This procedure induces trails of matched pairs of observations of variable length along the tilt series. The scheme that depicts which trails have representatives in which tilts is called ''occupancy graph''
 
[[ File:alignGUI_occupancySchemeOfTheResultsOfAStep.png |thumb|center| 600px|occupancy scheme of the results of a step]]
 
[[ File:alignGUI_occupancySchemeOfTheResultsOfAStep.png |thumb|center| 600px|occupancy scheme of the results of a step]]
[[ File:alignGUI_roughAlignmentSelectsManyPotentialTrails.png |thumb|center| 600px|rough alignment selects many potential trails ]]
+
[[ File:alignGUI_roughAlignmentSelectsManyPotentialTrails.png |thumb|center| 600px| The rough alignment selects many potential trails. On the x-axis is the tilt number, on the y-axis the trail number and each red dot marks the presence of a marker.]]
Note that these trails are very lacunary and don't conserve the identity of a trail; whenever a trail that represents an actual 3d marker gets interrupted, the same 3d marker might generate a different trail in other tilts.
+
Note that these trails are very lacunary and do not conserve the identity of a trail; whenever a trail that represents an actual 3d marker gets interrupted, the same 3d marker might generate a different trail in other tilts.
  
=== Selection of stable trails===
+
=== Keep long chains ===
We select the longest trails, which will then be used to generate a 3d model of the markers:
+
Click on ''keep long chains'' to select the longest trails, which will then be used to generate a 3d model of the markers:
[[ File:alignGUI_trailsWithFewMarkersAreThenRejected.png |thumb|center| 600px|trails with few markers are then rejected]]
+
[[ File:alignGUI_trailsWithFewMarkersAreThenRejected.png |thumb|center| 600px|short trails have been rejected]]
 
Note that from this step onwards, the toolbox of each step offers the possibility to revert the current markers to the output generated at that state. Thus, if further processing of the markers in a later step leads you to loose a good marker set, you can always come back to this point.
 
Note that from this step onwards, the toolbox of each step offers the possibility to revert the current markers to the output generated at that state. Thus, if further processing of the markers in a later step leads you to loose a good marker set, you can always come back to this point.
[[ File:alignGUI_noteWorkingMarkersCanBeRestoredToTheResultOfAPreviousStep.png |thumb|center| 600px|note working markers can be restored to the result of a previous step]]
+
[[ File:alignGUI_noteWorkingMarkersCanBeRestoredToTheResultOfAPreviousStep.png |thumb|center| 600px|working markers can be restored to the result of a previous step]]
  
 
=== Iterative reindexing ===
 
=== Iterative reindexing ===
The 3d model is projected on each micrograph. Observations that are closest (and below a distance threshold) to each reprojection are assigned to the corresponding 3d marker. This refines the 3d model, and the process is iterated till no further improvement occurs (or till a maximum of ten iterations is reached).
+
The 3d model is projected on each micrograph. Observations that are closest (and below a distance threshold) to each reprojection are assigned to the corresponding 3d marker. This refines the 3d model, and the process is iterated till no further improvement occurs (or till a maximum of ten iterations is reached). Keep the default parameters and click on ''iterative reindexing''.
 
[[ File:alignGUI_parametersForIterativeReindexing.png |thumb|center| 600px|parameters for iterative reindexing]]
 
[[ File:alignGUI_parametersForIterativeReindexing.png |thumb|center| 600px|parameters for iterative reindexing]]
[[ File:alignGUI_resultOfIterativeReindexing.png |thumb|center| 600px|result of iterative reindexing]]
+
[[ File:alignGUI_resultOfIterativeReindexing.png |thumb|center| 600px|Result of iterative reindexing. Missing markers have been added to the trails]]
  
 
===Tilt gap filling ===
 
===Tilt gap filling ===
''Dynamo'' can be directed to try to analyse those  tilts that wind up not having any marker  (and that have not been explicitly rejected by the user). The cloud of observations found on that micrograph will be compared to the cloud or reproductions generated by the current model.
+
Here, potential tilts with too few (or zero) markers can be rescued. The cloud of observations found on that micrograph will be compared to the cloud or reprojections generated by the current model. Click on ''fill tilts without markers''.
 
[[ File:alignGUI_parametersForTiltGapFiller.png |thumb|center| 600px|parameters for tilt gap filler]]
 
[[ File:alignGUI_parametersForTiltGapFiller.png |thumb|center| 600px|parameters for tilt gap filler]]
[[ File:alignGUI_occupancyAfterTiltGapFiller.png |thumb|center| 600px|occupancy after tilt gap filler]]
+
[[ File:alignGUI_occupancyAfterTiltGapFiller.png |thumb|center| 600px|Results of tilt gap filler. Note how (in this example) a lot of markers have been rescued on the first tilt.]]
[[ File:alignGUI_weHave66TracesAfterGapFiller.png |thumb|center| 600px|we have 66 traces after gap filler ]]
+
[[ File:alignGUI_weHave66TracesAfterGapFiller.png |thumb|center| 600px|Click on ''working markers'' to see that (in this example) we have 66 traces after tilt gap filling.]]
[[ File:alignGUI_resultsOfTiltGapFiller.png |thumb|center| 600px|results of tilt gap filler]]
+
[[ File:alignGUI_resultsOfTiltGapFiller.png |thumb|center| 600px|results of tilt gap filler (''show markers created by this step on tilt series'')]]
  
 
===Trail extension===
 
===Trail extension===
Line 182: Line 126:
  
 
Note that newly added trails have empty spaces in the central area, they are likely to be gold beads located far away from the tilt axes and which get into the field of view only for high tilts.  
 
Note that newly added trails have empty spaces in the central area, they are likely to be gold beads located far away from the tilt axes and which get into the field of view only for high tilts.  
[[ File:alignGUI_someOfTheNewTrailsHaveEmptySpaces.png |thumb|center| 600px|some of the new trails have empty spaces]]
+
[[ File:alignGUI_someOfTheNewTrailsHaveEmptySpaces.png |thumb|center| 600px|some of the new trails may have empty spaces]]
 
[[ File:alignGUI_resultOfTrailExtensor.png |thumb|center| 600px|result of trail extensor]]
 
[[ File:alignGUI_resultOfTrailExtensor.png |thumb|center| 600px|result of trail extensor]]
  
 
==Refinement ==
 
==Refinement ==
This area refines the positions of the gold beads. We will run all steps of this area pressing on a single button (the ''run this tab'' option at the bottom of the GUI)
+
This area refines the positions of the gold beads. We will run all steps of this area by pressing on a single button (the ''run this tab'' option at the bottom of the GUI). The button will remain green while the different steps are computed.
 
[[ File:alignGUI_tabForRefinement.png |thumb|center| 600px|tab for refinement]]
 
[[ File:alignGUI_tabForRefinement.png |thumb|center| 600px|tab for refinement]]
The Botton remains green while the succession of different steps gets computed.
+
Check the current refinement error:
[[ File:alignGUI_allStepsOfRefinementAreaAreBeingComputed.png |thumb|center| 600px|all steps of refinement area are being computed]]
 
 
[[ File:alignGUI_betterFiitingIsAttained.png |thumb|center| 600px|better fitting is attained]]
 
[[ File:alignGUI_betterFiitingIsAttained.png |thumb|center| 600px|better fitting is attained]]
We can use the scissors icon in order to "trim" observations (or full markers) with too high residuals
+
Optionally, we can use the scissors icon to "trim" observations (or full markers) with too high residuals. Click on the scissor icon. A window pops up with options to set different thresholds. Click on ''Tools'' and then ''Show residuals per observation'':
 
[[ File:alignGUI_trimmingOfMarkersCanBeAccessedThroughTheScissorIcon.png |thumb|center| 600px|trimming of markers can be accessed through the scissor icon]]
 
[[ File:alignGUI_trimmingOfMarkersCanBeAccessedThroughTheScissorIcon.png |thumb|center| 600px|trimming of markers can be accessed through the scissor icon]]
 +
You get a plot with the residuals for each observation ordered by the value of the residual. You can use this plot to estimate the residual distribution and to manually adapt the residual thresholds, if necessary.
 
[[ File:alignGUI_residualsAtEachObservation.png |thumb|center| 600px|residuals at each observation]]
 
[[ File:alignGUI_residualsAtEachObservation.png |thumb|center| 600px|residuals at each observation]]
  
 
== Alignment ==
 
== Alignment ==
This area creates the aligned tilt series. It creates two version, a binned one and one with full resolution pixel size. Note that the binning level of this one is independent of the binning level used in the detection step.
+
The final step of the alignment procedure is to create the aligned tilt series using the coordinates of the markers that were optimized until now.
  
 
=== Fix the markers ===
 
=== Fix the markers ===
Here we just freeze the working markers and create a marker file that will define all subsequent steps of alignment and reconstruction.
+
Here we freeze the final version of the working markers and create a marker file that will define all subsequent steps of alignment and reconstruction. Click on ''fix alignment markers''.
  
 
[[ File:alignGUI_areaForCreationOfAnAlignedTiltSeries.png |thumb|center| 600px|area for creation of an aligned tilt series]]
 
[[ File:alignGUI_areaForCreationOfAnAlignedTiltSeries.png |thumb|center| 600px|area for creation of an aligned tilt series]]
  
 
=== Align the tilt series ===
 
=== Align the tilt series ===
We then proceed to the alignment task itself, generating the binned and the full resolution stacks.
+
We then proceed to the alignment task itself, generating full resolution aligned tilt series. We also have the option to generate a binned version of the aligned tilt series that can be used for quick visualization or to create binned reconstructions in the next steps. Click on ''Create full-size aligned stack'' and ''Create binned aligned stack''. For a quick assessment, open the binned stack in tomoshow:
 
[[ File:alignGUI_visualizationOfAlignedTiltStack.png |thumb|center| 600px|visualization of aligned tilt stack]]
 
[[ File:alignGUI_visualizationOfAlignedTiltStack.png |thumb|center| 600px|visualization of aligned tilt stack]]
 
[[ File:alignGUI_binnedAlignedTiltStack.png |thumb|center| 600px|binned aligned tilt stack]]
 
[[ File:alignGUI_binnedAlignedTiltStack.png |thumb|center| 600px|binned aligned tilt stack]]
 +
A quick way to gauge the quality of the alignment is by checking that the beads move nicely horizontally over all tilts.
  
 
== CTF  ==
 
== CTF  ==
In this GUI CTF estimation and correction are delegated to CTFFIND4 and Imod's ctfphaseflip programs respectively.
+
In this GUI, CTF estimation and correction are done using CTFFIND4 and IMOD's ctfphaseflip. This part will not be used in the workshop.
''This part will not be used in the workshop.'''
 
  
 
== Reconstruction ==
 
== Reconstruction ==
  
 
=== Binned reconstruction ===
 
=== Binned reconstruction ===
This version of ''Dynamo'' fixes the tomogram center at the 3d reconstruction center of the aligned stack.
+
Click on ''binned reconstruction'' to create a binned tomogram using the previously created binned aligned tilt series. Adapt the tomogram height if necessary (the tomogram center is automatically located at the center of the 3d reconstruction of the aligned stack). Select ''binned SIRT-like reconstruction'' to generate a high contrast tomogram. Run the step ''binned reconstruction'' and visualize the resulting binned SIRT-like tomogram:
[[ File:alignGUI_reconstructionArea.png |thumb|center| 600px|reconstruction area]]
+
[[ File:alignGUI_chooseAVisualization.png |thumb|center| 700px|choose a visualization]]
[[ File:alignGUI_chooseAVisualization.png |thumb|center| 600px|choose a visualization]]
+
[[ File:alignGUI_binnedWeightedBackProjectionReconstruction.png |thumb|center| 600px|binned SIRT-like reconstruction]]
[[ File:alignGUI_binnedWeightedBackProjectionReconstruction.png |thumb|center| 600px|binned weighted back projection reconstruction]]
 
  
 
=== Full sized reconstruction ===
 
=== Full sized reconstruction ===
The full sized reconstruction can be created related to a coordinate directly read in the binned tomogram, a feature useful for creating reconstructions localised on an area of interest. In this case, we will focus on one of the VLPS.
+
The full sized reconstruction is the one you will be using for later image processing steps that require the highest resolution possible (e.g., subtomogram averaging refinements). It can be created related to a coordinate directly read from the binned tomogram. This feature is useful to create reconstructions localised on an area of interest. In this tutorial, we will focus on one of the VLPs, resulting in a tomogram containing just one VLP. First, use the pin tool to select a point on the binned reconstruction. We use the previously opened binned SIRT-like tomogram for this:
 +
[[ File:alignGUI_thePinToolIsUsefulForLocatingIsolatedPoints.png |thumb|center| 600px|Click in the center of an object of interest using the pin tool.]]
  
Use the pin tool to select a point on the binned reconstruction.
+
We check the size of the structure by using the measuring tool:
[[ File:alignGUI_thePinToolIsUsefulForLocatingIsolatedPoints.png |thumb|center| 600px|the pin tool is useful for locating isolated points]]
 
 
 
[[ File:alignGUI_wePutAPinInTheCenterOfAnObjectOfInterest.png |thumb|center| 600px|we put a pin in the center of an object of interest]]
 
 
 
We also check the size of a structure of interest. We measure it in the binned tomogram.
 
 
[[ File:alignGUI_weCheckTheDistanceOfTheObjectOfInterest.png |thumb|center| 600px|we check the distance of the object of interest]]
 
[[ File:alignGUI_weCheckTheDistanceOfTheObjectOfInterest.png |thumb|center| 600px|we check the distance of the object of interest]]
[[ File:alignGUI_weUseTheOptionOfExpressingTheCenterOfTheFullSizedTomogramInTermsOfAPositionInTheBinnedTomogram.png |thumb|center| 600px|we use the option of expressing the center of the full sized tomogram in terms of a position in the binned tomogram]]
 
 
A viewer system is not yet connected to this step, we can however use the common inspection tools to fetch the outcome (the full sized WBP)
 
[[ File:alignGUI_asAnSpecificViewerIsNotAvailableForThisStepWeCanUseTheToolboxOfTheStepsItems.png |thumb|center| 600px|as an specific viewer is not available for this step we can use the toolbox of the steps items]]
 
[[ File:alignGUI_theFullSizeReconstructionIncludesOnlyTheObjectOfInterest.png |thumb|center| 600px|the full size reconstruction includes only the object of interest]]
 
 
As it is a full resolution reconstruction, it is much noisier than the corresponding area in the binned tomogram. You probably need to use a rather deep bandpass (around 0.2 of nyquist) in order to distinguish it.
 
[[ File:alignGUI_YouMightNeedToActivateTheBandpassToSeeTheVLP.png |thumb|center| 600px|You might need to activate the bandpass to see the VLP ]]
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
  
 +
Before generating the final localized tomogram, we need to set some parameters. We take the measured distance of 83 pixels and compensate for the binning factor of two, leading to 83*2*2=332 pixels. For the dimensions of the tomogram we add some spare space and define the x,y and z length to be 500 pixels. Select ''use center of binned coordniates'' and enter the values we measured before with the pin tool. Select ''WBP reconstruction'' and run the step ''full reconstruction''. Use the inspection tools to fetch the result:
 +
[[ File:alignGUI_asAnSpecificViewerIsNotAvailableForThisStepWeCanUseTheToolboxOfTheStepsItems.png |thumb|center| 700px|as an specific viewer is not available for this step we can use the toolbox of the steps items]]
  
]
+
As it is a WBP reconstruction, it is much noisier than the corresponding area in the SIRT-like tomogram. Adjust the contrast, lowpass filter and thickness to increase the signal:
 +
[[ File:alignGUI_theFullSizeReconstructionIncludesOnlyTheObjectOfInterest.png |thumb|center| 600px|the full size WBP reconstruction includes only the object of interest]]

Latest revision as of 23:05, 5 August 2021

Dynamo includes a package for automated alignment and reconstruction of tilt series. This walkthrough guides users through a series of steps on how to use it in GUI-based mode. The walkthrough on command line based tilt series alignment shows how to operate this procedure through the command line.

Data

We use a binned version of the first tilt series from the EMPIAR entry 10164, depicting a set of virus like particles (VLP). Visit the course materials to access the tilt series or download it directly from EMPIAR. You can also follow this tutorial using your own data.

Creating an alignment workflow

Start the alignment GUI by typing the command dtsa (short form of dynamo_tilt_series_alignment) into the Dynamo command line:

 dtsa

A pop-up window appears. Click on "create a new workflow":

creating a workflow

Here, the input data for the workflow (tilt series, angles, ...) is introduced:

input workflow parameters

The option configuration file is used to pass information on possible tilts that should be discarded. Here, we leave this option untouched. Press on create workflow, to open the alignment GUI.

GUI description

The alignment GUI guides the process of generating gold bead coordinates, aligning the tilt series, correcting the CTF and to finally reconstruct the tomogram. It has 5 parts. Starting from top to bottom, they are:

  1. Menu: Different tools for visualization and customization of general settings. The option Settings > edit[general] for example lets you adjust the pixel size and different acquisition parameters such as Cs etc. Here we keep the default settings.
  2. Icons: Quick access to tools such as dmarkers, occupancy plot, gold bead gallery, gold bead targeting, fitting and marker selection. Their functionalities will be explained later.
  3. Tabs: The tabs describe the 6 execution areas that are executed from left to right: detection, indexing, refinement, alignment, ctf, reconstruction. Each execution areas contains a set of processing steps.
  4. Steps: The processing steps of each tab are run sequentually from top to bottom. However, after execution, the steps can be revisited to change or update information. Steps need to be active (letters highlighted in black, not grey) to be executed. They are automatically activated once the previous steps have been executed or the required input is given. Input parameters are edited by clicking the plus button (+) on the left of the step button. When a step is under execution, the text of the button becomes green. After execution the text becomes bold. At the right side of each of the steps, a tool icon is located. The tool icon can be right-clicked to provide a pop-up menu with specific tools for each step. Next to the tool icon there is the View button which provides specific tools to visualize and inspect the results of each step.
  5. Run options: Options to run multiple steps.
alignment GUI before completion of any step

Processing

We are now ready to run all the processing steps. In the following, each step is described in detail.

Detection

The goal of this step is to find the positions of the gold bead projections in the tilt series. To do that, Dynamo cross correlates each micrograph against a 2D model of the gold bead and then locate the cross correlation peaks, which provide putative positions for actual projections of gold beads. Further analysis of these peaks will determine which ones correspond to actual gold bead projection (peak feature and peak selection).

Visualization matrix

This step just creates a binned version of the tilt series stack to accelerate visualizations that will be invoked in a later point. If you click on the box standing on the left of creates a binned tilt series button with a + symbol, an option box appears, indicating the default binning factor which is 2. This parameter can be changed by users. Binning factor 2 in Dynamo is equivalent to a binning factor of 4 in IMOD. Click on creates a binned tilt series button to create a binned tilt series stack. Rightclick on View button to visualize the newly binned tilt series stack in tmshow. In the tmshow window, visualization parameters as binning, band and thickness can be adjusted if needed.

Now, we are going to measure the diameter of the gold beads precisely in tmshow in the original (unbinned) tilt series data. To do that, press the IO in the upper toolbox of the alignment workflow GUI, then go to tilt series matrix > inspect > On-the-fly view [dtmshow]. The tmshow window opens, adjust contrast if necessary for better visualization.

open unbinned tilt series data to measure gold bead diameter

Zoom in a gold bead by holding the shift key and scrolling up your mouse wheel (alternatively, use the icon with the red rectangle. Avoid using the lens icons). Adjust the visualized region by holding the shift key again and dragging with right mouse click on tmshow window.

open tomoshow on the full sized matrix to check the size of a gold bead

Once you have located the gold bead of interest, click on the ruler icon at the top toolbox to activate length measurements on the displayed image data. Draw a line across the gold bead, and the length of the line appears in pixels, which correspond to the gold bead diameter. In this case, the length of this line should be 32 pixel aprox.

the ruler tool allows to measure in pixels the diameter of a gold bead

Detection of gold beads

Here we provide the size of the gold beads. The measured diameter of the gold beads is 32 pixels in the original image, therefore its radius is 16 pixels. We introduce the bead radius and some other parameters (see screenshot) to create an accurate initial gold bead template (here, the binning is just used internally for quick and robust bead detection).

the plus icon at corner reveals the values of the parameters of a given step.

After pressing the button detect peaks, Dynamo creates a raw cross-correlation matrix based on the initial template which is binary and extracts 300 cc-peaks by default. Then, the extracted peaks are averaged to create a refined template. The refined template is used to compute a refined cc-matrix from which 300 peaks are extracted and used as initial working markers. Once the detect peaks button turns black bold, the computation is finished. Right click on the view area of the step to access the different visualisation options.

Secondary click the view panel to open a menu to access different results of a step. Marker sets can be viewed as cloud of points.

This view shows the overlaid sum of the x,y coordinates of all cc peaks from all micrographs in the tilt series.

cloud of unindexed cc peaks

We can also check the positions of all the detected spots on each micrograph:

shows the detected spots on the video GUI

This GUI does not allow for edition of the markers. Tip: enable auto option at the bottom center of the visualization GUI to facilitate tilt series visualization.

static video annotations cannot be edited

The neighbourhood of cross correlation peaks can be visualised individually as 2D patches from the tilt series image data:

Each patch will be a putative gold bead
the set of cc peaks is shown as a montage

This visualization is normally a good gauge on how the detection procedure has worked. The example below would be an example of detection procedure carried with wrong search parameters

example of failed gold bead detection

Note that you have the option of testing your parameters against a single image before launching a full detection procedure. Right-click on the tool icon on the right side of detect peaks pushbutton.

toolbox option for testing the parameters on a central image
results of testing on a central image

Computation of observation features

At this stage, the observed cross correlation peaks are analysed. A merit figure is assigned to each peak based on a rotational merit that measures its rotational symmetry. Press the button Peak features to compute the rotational symmetry of the current working markers. The working markers that do not meet rotational symmetry will be discarded. The results can be depicted on patches representing individual gold beads:

patch gallery for features
the gallery shows the detected gold beads, the CC score (yellow) and the rotational merit (red). Tilt number is purple.

or extracting properties of the whole data set.

feature plotter allows checking of the feature values on the whole set
We select scatter, set the y axis to rotationaMerit and then press Scatter 2d.
We get a scatter plot of each peak and its relation to the CC and rotational merit score. Seeing two distinct populations means that the majority of the gold beads are distinguishable from other signal such as background.

Selection of best gold beads

The population with the best markers is automatically selected by clicking on Peak selection. The results can be inspected again in a plot

check observation cloud after selection or clustering
observation cloud has been thinned

or in the micrographs themselves:

check marker positions after clustering
marker positions after clustering

Indexing

The goal of this execution area is to find trails of gold beads, i.e., we connect the observations from different micrographs that correspond to the same gold bead. All the steps in this area will constantly update a set of working markers.

Rough Alignment

In this step, micrographs are compared pair-wise. For each couple, we look for the relative shifts that generate the most matches of cross-correlation peaks between one micrograph and the next one. This procedure induces trails of matched pairs of observations of variable length along the tilt series. The scheme that depicts which trails have representatives in which tilts is called occupancy graph

occupancy scheme of the results of a step
The rough alignment selects many potential trails. On the x-axis is the tilt number, on the y-axis the trail number and each red dot marks the presence of a marker.

Note that these trails are very lacunary and do not conserve the identity of a trail; whenever a trail that represents an actual 3d marker gets interrupted, the same 3d marker might generate a different trail in other tilts.

Keep long chains

Click on keep long chains to select the longest trails, which will then be used to generate a 3d model of the markers:

short trails have been rejected

Note that from this step onwards, the toolbox of each step offers the possibility to revert the current markers to the output generated at that state. Thus, if further processing of the markers in a later step leads you to loose a good marker set, you can always come back to this point.

working markers can be restored to the result of a previous step

Iterative reindexing

The 3d model is projected on each micrograph. Observations that are closest (and below a distance threshold) to each reprojection are assigned to the corresponding 3d marker. This refines the 3d model, and the process is iterated till no further improvement occurs (or till a maximum of ten iterations is reached). Keep the default parameters and click on iterative reindexing.

parameters for iterative reindexing
Result of iterative reindexing. Missing markers have been added to the trails

Tilt gap filling

Here, potential tilts with too few (or zero) markers can be rescued. The cloud of observations found on that micrograph will be compared to the cloud or reprojections generated by the current model. Click on fill tilts without markers.

parameters for tilt gap filler
Results of tilt gap filler. Note how (in this example) a lot of markers have been rescued on the first tilt.
Click on working markers to see that (in this example) we have 66 traces after tilt gap filling.
results of tilt gap filler (show markers created by this step on tilt series)

Trail extension

This step aims at locating trails that have been excluded by previous steps of the alignment procedure, or that were never detected as trails because previous, lesser quality stages of the 3d model were not able to analyse the observations correctly.

parameters of tilt extensor
trail extensor adds a few new trails

Note that newly added trails have empty spaces in the central area, they are likely to be gold beads located far away from the tilt axes and which get into the field of view only for high tilts.

some of the new trails may have empty spaces
result of trail extensor

Refinement

This area refines the positions of the gold beads. We will run all steps of this area by pressing on a single button (the run this tab option at the bottom of the GUI). The button will remain green while the different steps are computed.

tab for refinement

Check the current refinement error:

better fitting is attained

Optionally, we can use the scissors icon to "trim" observations (or full markers) with too high residuals. Click on the scissor icon. A window pops up with options to set different thresholds. Click on Tools and then Show residuals per observation:

trimming of markers can be accessed through the scissor icon

You get a plot with the residuals for each observation ordered by the value of the residual. You can use this plot to estimate the residual distribution and to manually adapt the residual thresholds, if necessary.

residuals at each observation

Alignment

The final step of the alignment procedure is to create the aligned tilt series using the coordinates of the markers that were optimized until now.

Fix the markers

Here we freeze the final version of the working markers and create a marker file that will define all subsequent steps of alignment and reconstruction. Click on fix alignment markers.

area for creation of an aligned tilt series

Align the tilt series

We then proceed to the alignment task itself, generating full resolution aligned tilt series. We also have the option to generate a binned version of the aligned tilt series that can be used for quick visualization or to create binned reconstructions in the next steps. Click on Create full-size aligned stack and Create binned aligned stack. For a quick assessment, open the binned stack in tomoshow:

visualization of aligned tilt stack
binned aligned tilt stack

A quick way to gauge the quality of the alignment is by checking that the beads move nicely horizontally over all tilts.

CTF

In this GUI, CTF estimation and correction are done using CTFFIND4 and IMOD's ctfphaseflip. This part will not be used in the workshop.

Reconstruction

Binned reconstruction

Click on binned reconstruction to create a binned tomogram using the previously created binned aligned tilt series. Adapt the tomogram height if necessary (the tomogram center is automatically located at the center of the 3d reconstruction of the aligned stack). Select binned SIRT-like reconstruction to generate a high contrast tomogram. Run the step binned reconstruction and visualize the resulting binned SIRT-like tomogram:

choose a visualization
binned SIRT-like reconstruction

Full sized reconstruction

The full sized reconstruction is the one you will be using for later image processing steps that require the highest resolution possible (e.g., subtomogram averaging refinements). It can be created related to a coordinate directly read from the binned tomogram. This feature is useful to create reconstructions localised on an area of interest. In this tutorial, we will focus on one of the VLPs, resulting in a tomogram containing just one VLP. First, use the pin tool to select a point on the binned reconstruction. We use the previously opened binned SIRT-like tomogram for this:

Click in the center of an object of interest using the pin tool.

We check the size of the structure by using the measuring tool:

we check the distance of the object of interest

Before generating the final localized tomogram, we need to set some parameters. We take the measured distance of 83 pixels and compensate for the binning factor of two, leading to 83*2*2=332 pixels. For the dimensions of the tomogram we add some spare space and define the x,y and z length to be 500 pixels. Select use center of binned coordniates and enter the values we measured before with the pin tool. Select WBP reconstruction and run the step full reconstruction. Use the inspection tools to fetch the result:

as an specific viewer is not available for this step we can use the toolbox of the steps items

As it is a WBP reconstruction, it is much noisier than the corresponding area in the SIRT-like tomogram. Adjust the contrast, lowpass filter and thickness to increase the signal:

the full size WBP reconstruction includes only the object of interest