Split stack using select file?

8 posts / 0 new
Last post
tarek
Split stack using select file?

Dear all,
I have completed refinement of a dataset with a large box and used the same particles with smaller box for classification into sub-sets.
Now I have alignment parameters at a large box and want to create a new sub-stack according to the class-results.
How would you accomplish this keeping alignment parameters?

Best,
tarek

timgrant
Hi Tarek, 

Hi Tarek, 

Was this done in cisTEM? In which case you can great the new refinement pakackage for the desired class.  Then export the parameters from the refined class and then import a new refinement package, selecting your exported parameters, and the stack with the larger box size.

Would that work?

Cheers,

Tim

tarek
Hi Tim,

Hi Tim,

thanks for your reply.

It was done with cistem but this way should not work. Particle numbers between the exported parameters file and the stack with larger box will differ.

Let's assume particle no. with large box are 200,000. Selected particles after classification with small box are 20,000.

This 20,000 are of course a fraction of 200,000 but I would have to selectively pick the parameters out of the 200k par file.

 

Cheers,

tarek

timgrant
Hi Tarek,

Hi Tarek,

I'm not sure i understand? When you created your new refinement package with 20,000 particles it should have carried over the parameters.  You can then export the parameters of these 20k and then do what I suggest above?  Will that not work?

Tim

tarek
Hi Tim,

Hi Tim,

the 20,000 particles were from decimated data and only refined at large pixel size.

I need to take over the alignment parameters of the large box 200k stack but only take the 20k particles which were identified using the small box.

The point is: the 20 k particles have too crude alignment parameters at small box. When just exchanging the particles with the large boxed ones, I need to re-align too liberately (interpolation errors...) what will cause even poorer alignment. It's a tricky dataset, few particles, poor signal  :(

It would be great to have something similar like "original image name" in the par files that would be as in the imported the star file ;)

timgrant
Hi Tarek,

Hi Tarek,

Aaah, I see the problem now.  Yes, there is no easy way to do this.  You could probably do it through a combination of database queries and scripting, but it would be somewhat complicated.

Tim

tarek
thanks for the chat, I will

thanks for the chat, I will have to go the hard way...

Daniel Asarnow
Hey tarek,

Hey tarek,

You can use a few programs from my github (master branch) to accomplish it pretty easily, though with some extra files to track.

  1. par2star.py Assets/Parameters/output_N_*.par output_N.star --apix 1.0 --cs 2.7 --kv 200 --ac 0.1 --stack refpkg.mrcs
  2. star.py --augment --merge-fields=rlnClassNumber --merge-key=ucsfImageIndex --merge-source output_N.star refpkg.star refpkg_classes.star OR star.py output_N.star refpkg.star refpkg_classes.star --class K --copy-alignments output_N.star
  3. stack.py <star> <mrcs out> --star <star out> [--class K]

Maximum occupancy is used to assign particle parameters and class number. You can choose how to merge the information using all the options in star.py, or expand the --merge-fields example. Let me know if it does/doesn't work for you.

Log in or register to post comments