Difference between revisions of "Operations on volumes"

From Dynamo
Jump to navigation Jump to search
Line 50: Line 50:
  
 
===Flipping  along a direction===
 
===Flipping  along a direction===
You can use the Matlab command {{t|flip}}. This requires using {{t|dread}} <tt>dread</tt> to bring the contents of a file into the Matlab workspace.
+
You can use the Matlab command {{t|flip}}. This requires using {{t|dread}} <tt>dread</tt> dread to bring the contents of a file into the Matlab workspace.
  
 
Also, you can use the ''Dynamo'' command {{t|dpktomo.tools.flipYZ}} to flip the contains of a file controlling the maximum amount of memory used as a given time.
 
Also, you can use the ''Dynamo'' command {{t|dpktomo.tools.flipYZ}} to flip the contains of a file controlling the maximum amount of memory used as a given time.

Revision as of 12:05, 15 February 2017

There are different operations that are frequently required when processing 3d data, both for small volumes (as templates, averages and particle) or tomograms. Dynamo and Matlab itself provide several commands.

Using the functions

In general, volumes need to be read into the worskpace, creating a variable. The operator can then be used onto the variable and then the result is written back into a file. For instance, to bin the contents of the file myFile.mrc and write the result into the file binned.mrc

 a = dread('myFile.mrc'); 
 b = dbin(a,2);  
 dwrite(b,'binned.mrc'); 

Most processing options accept as input indistinctly a filename or a variable. Note also that dread accepts the alternative command line syntax:

 dread myFile.mrc a 

Dealing with large tomograms

Reading full tomograms into the worskpace might be a bad idea, leading to overcrowding you RAM memory and freezing your system. In these situations, it is frequently a good approach to let Dynamo divide automatically the tomogram and operate separately on the chunks (we call this "chunk partition" in the documentation). There are two different devices for this task:

Commands oriented to chunk processing

Some tasks already integrate the operation through a chunk partition operated on the fly, for instance the binning of tomograms, which allows the command dpktomo.partition.bin

Partition objects

Technical information

For operations that you want to define yourself, you can use the class behind dpktomo.partition.bin, i.e. mbgeom.partition.GridPartition

Functionalities

Binning

Command dbin.

Normalization

Command dynamo_normalize_roi. You can also pass a mask, so that the normalization will ensure zero mean and std=1 inside the mask.

Filters

The bandpass filter can be performed with the command dynamo_bandpass .html dbandpass. The extent of the band is given in fourier pixels.

Median filtering is available through dynamo_median.html dmedian.

Contrast inversion

You can just use the minus operator in Matlab:

 dread file.mrc a 
 b = -a;

possibly normalizing again.

Flipping along a direction

You can use the Matlab command flip. This requires using dread dread dread to bring the contents of a file into the Matlab workspace.

Also, you can use the Dynamo command dpktomo.tools.flipYZ to flip the contains of a file controlling the maximum amount of memory used as a given time.

Flipping two directions

You can use the Matlab commad permute to interchange two directions in a volume.

Rotations and shifts

dynamo_rot.html drot,dynamo_move .html dmove

Symmetrization

The general function is dynamo_sym.html dsym. The notation for the symmetry operators in this command is the one used when defining project parameters.

Advanced functions

A (somewhat) GUI version is available as the command dvolume. The command line counterpart is dynamo_particle.html dparticle, which accepts several operators (binning, symmetrization, etc). The notation of the operators in dparticle is the same used in other instances of Dynamo (normally under the name actions or particle actions) when volumetric transformations are operated on the fly. For instance when depicting the contents of data folders with ddbrowse.