Skip to content
homurachan edited this page Oct 21, 2022 · 5 revisions

How to combine blocks to a whole map

It should be noted that a block after refinement by RELION or cryoSparc could have a small rotation to original maps.

Check https://github.com/homurachan/Block-based-recontruction/wiki/How-to-fit-block-into-original-map for how to correct small rotation.

There are 2 kinds of situation.

First is to combine blocks into a asymmetric unit. In this case you should download this below which is a standalone program.

https://github.com/homurachan/Block-based-recontruction/blob/master/addup_many_part_v5_diffsize.exe

Input and output are shown in the figure below.

Generally speaking, the max radius of interpretation should be significant smaller than boxsize/2 to avoid fading away of density at the edge.

The min radius means we keep all density without interpretation within the radius in pixels.

////////////////////////////////////////

The second situation is to expand blocks into a icosahedral map, e.g. https://github.com/homurachan/Block-based-recontruction/issues/1 . In this case you should download these below which are also standalone programs.

https://github.com/homurachan/Block-based-recontruction/blob/master/addup_many_part_into_full_icos_v2.exe and https://github.com/homurachan/Block-based-recontruction/blob/master/make_mrc_slice_into_3d_mrc.exe

The outputs of addup_many_part_into_full_icos_v2.exe are slices in Z-axis. The reason is that the program runs quite slow. By producing slices you can run the program multiple times to greatly save time.

Then using make_mrc_slice_into_3d_mrc.exe to read those slices and generating full-size map.

"Enter images root name:

asdf

Enter ny of full size output image:

1024

Enter output

asdf_out.mrc"

The output map asdf_out.mrc is what you need. Remember to check edge of each blocks, there might be some artifacts if the radius is not carefully chosen.

It is suggested to run "e2proc3d.py blockN.mrc blockN_norm.mrc --process=normalize.circlemean:radius=MAX_RADIUS" before combining blocks.