Project

General

Profile

Creating Hole Template » History » Revision 2

Revision 1 (Amber Herold, 04/23/2010 08:34 AM) → Revision 2/5 (Amber Herold, 04/27/2010 11:55 AM)

h1. Creating Template 


 Any image can be used as a template for the hole finder. If you create your own, you just need to point the node preference to the right mrc file. Our philosophy for the template distributed with the Leginon package (your leginon python source path/holetemplate.mrc) was to make a hole template that had the average features of a 2 micron hole covered in vitreous ice. To accomplish this, we took fifty images of holes acquired by Leginon at 5000X magnification. Holes were then boxed out of these images using the boxer utility in the EMAN package making sure that the holes were in the center of the box. 


 !http://emg.nysbc.org/software/leginon/images/images/JAHCHole.png! Figuer 1-Boxing a hole 


 The hole images were then assembled into a single Imagic stack. Since the hole images were well centered, they were then averaged to produce an image of an average hole (fig. 2). Finally, the average hole image was rotationally averaged to produce a symmetrical generalized hole image (fig. 3). 


 !http://emg.nysbc.org/software/leginon/images/images/JAHCAvg.png!Figuer 2-Averager of 50 holes !http://emg.nysbc.org/software/leginon/images/images/JAHCRotAvg.png! Figure 3-Rotaionally averaged final template. 


 The following is the EMAN program and parameters associated with each step. 


 *    For 10-50 hole images, run boxer to box out holes and save them as a pair of IMAGIC file (hole.hed and hold.img). 


 <pre>boxer <hole_image.mrc></pre> 


 *    Average hole images and save the output as an MRC image file. 


 <pre>proc2d hole.hed average.mrc average</pre> 


 *    Rotationally average the average hole image and save as an MRC image file. 


 <pre>proc2d average.mrc rotav.mrc rotav</pre> 


 *    The result is the new hole template. 



 ______ 

 [[Summary of this application - MSI-T|< Summary of this application]] | [[JAHC finder nodes set-up|JAHC finder nodes set-up >]] 

 ______