User Tools

Site Tools


pipeline:isac:sxisac2

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
Next revision Both sides next revision
pipeline:isac:sxisac2 [2018/02/21 16:57]
moriya
pipeline:isac:sxisac2 [2018/06/20 13:12]
127.0.0.1 external edit
Line 80: Line 80:
 \\ \\
 ===== Description ===== ===== Description =====
-=== Time and Memory ===+ 
 +\\ 
 +==== Method ==== 
 +The program will perform the following steps (to save computation time, in case of inadvertent termination, i.e. power failure or other causes, the program can be restarted from any saved step location, see options) 
 + 
 +  - The images in the input stacked will be phase-flipped. 
 +  - The data stack will be pre-aligned (output is in subdirectory 2dalignment, in particular it contains the overall 2D average aqfinal.hdf, it is advisable to confirm it is correctly centered). 
 +    * In case 2dalignment directory exists steps 1 and 2 are skipped.  
 +  - The alignment shift parameters will be applied to the input data. 
 +  - **IMPORTANT**: Input aligned images will be resized such that the original user-provided radius will be now target_radius and the box size target_nx + xr - 1.  The pixel size of the modified data is thus original_pixel_size * original_radius_size / target_radius. 
 +    * The pseudo-code for adjusting the size of the radius and the size of the images is as follows: 
 +    * shrink_ratio = target_radius / original_radius_size 
 +    * new_pixel_size = original_pixel_size * shrink_ratio 
 +    * if shrink_ratio is different than 1: resample images using shrink_ratio 
 +    * if new_pixel_size > target_nx : cut image to be target_nx in size 
 +    * if new_pixel_size < target_nx : pad image to be target_nx in size 
 +    * The target_radius and target_nx options allow the user to finely adjust the image so that it contains enough background information. 
 +  - The program will iterate through generations of ISAC2 by alternating two steps. The outcome of these two steps is in subdirectory generation_*** (stars replaced by the current generation number). 
 +    *  Calculation of candidate class averages. 
 +    *  Calculation of validated class averages.  
 +  - The program will terminate when it cannot find any more reproducible class averages. 
 +  - If no restart option is given the program will pick-up from the last saved point. 
 + 
 +Also see the reference below. 
 + 
 +\\ 
 +==== Time and Memory ====
  
 Unfortunately, ISAC2 is very time- and memory-consuming. (OBSLETE?!! For example, on my cluster, it takes 15 hours to process 50,000 64x64 particles on 256 processors. Therefore, before embarking on the big dataset, it is recommended to run a test dataset (about 2,000~5,000 particles) first to get a rough idea of timing.) Unfortunately, ISAC2 is very time- and memory-consuming. (OBSLETE?!! For example, on my cluster, it takes 15 hours to process 50,000 64x64 particles on 256 processors. Therefore, before embarking on the big dataset, it is recommended to run a test dataset (about 2,000~5,000 particles) first to get a rough idea of timing.)
  
 \\ \\
-=== Retrieval of images signed to selected group averages ===+==== Retrieval of images signed to selected group averages ====
   - Open in e2display.py file class_averages.hdf located in the main directory.   - Open in e2display.py file class_averages.hdf located in the main directory.
   - Delete averages whose member particles should not be included in the output.    - Delete averages whose member particles should not be included in the output. 
Line 97: Line 123:
  
 \\ \\
-=== RCT information retrieval ===+==== RCT information retrieval ====
 Let us assume we would want to generate a RCT reconstruction using as a basis group number 12 from ISAC2 generation number 3.  We have to do the following steps: Let us assume we would want to generate a RCT reconstruction using as a basis group number 12 from ISAC2 generation number 3.  We have to do the following steps:
  
Line 111: Line 137:
   - Extract the needed alignment parameters.  The order is phi,sx,sy,mirror.  sx and mirror are used to transfer to tilted images.   - Extract the needed alignment parameters.  The order is phi,sx,sy,mirror.  sx and mirror are used to transfer to tilted images.
     * $ sxheader.py  group3_12.12.hdf  --params=xform.align2d  --export=params_group3_12.txt     * $ sxheader.py  group3_12.12.hdf  --params=xform.align2d  --export=params_group3_12.txt
- 
-\\ 
-==== Method ==== 
-The program will perform the following steps (to save computation time, in case of inadvertent termination, i.e. power failure or other causes, the program can be restarted from any saved step location, see options)  : 
- 
-  - The images in the input stacked will be phase-flipped. 
-  - The data stack will be pre-aligned (output is in subdirectory 2dalignment, in particular it contains the overall 2D average aqfinal.hdf, it is advisable to confirm it is correctly centered). 
-    * In case 2dalignment directory exists steps 1 and 2 are skipped.  
-  - The alignment shift parameters will be applied to the input data. 
-  - **IMPORTANT**: Input aligned images will be resized such that the original user-provided radius will be now target_radius and the box size target_nx + xr - 1.  The pixel size of the modified data is thus original_pixel_size * original_radius_size / target_radius. 
-    * The pseudo-code for adjusting the size of the radius and the size of the images is as follows: 
-    * shrink_ratio = target_radius / original_radius_size 
-    * new_pixel_size = original_pixel_size * shrink_ratio 
-    * if shrink_ratio is different than 1: resample images using shrink_ratio 
-    * if new_pixel_size > target_nx : cut image to be target_nx in size 
-    * if new_pixel_size < target_nx : pad image to be target_nx in size 
-    * The target_radius and target_nx options allow the user to finely adjust the image so that it contains enough background information. 
-  - The program will iterate through generations of ISAC2 by alternating two steps. The outcome of these two steps is in subdirectory generation_*** (stars replaced by the current generation number). 
-    *  Calculation of candidate class averages. 
-    *  Calculation of validated class averages.  
-  - The program will terminate when it cannot find any more reproducible class averages. 
-  - If no restart option is given the program will pick-up from the last saved point. 
- 
-Also see the reference below. 
  
 \\ \\
pipeline/isac/sxisac2.txt · Last modified: 2020/08/19 16:40 by fschoenfeld