-
Notifications
You must be signed in to change notification settings - Fork 205
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
relion_helix_toolbox doesn't work #122
Comments
Original comment by Jing_basel (Bitbucket: Jing_basel, GitHub: Unknown): Hi Sjors, I trired: "relion_helix_toolbox" I got the same error message: Only one option can be specified at one time! valid_options = 127
|
Original comment by Jing_basel (Bitbucket: Jing_basel, GitHub: Unknown): Hi Sjors, I also tried the examples copied from the wiki page: relion_helix_toolbox --cylinder --o ref.mrc --boxdim 300 --cyl_outer_diameter 200 --cyl_inner_diameter -1 --angpix 1.34 --sphere_percentage 0.9 --width 5 relion_helix_toolbox --simulate_helix --o VipASimRefCyl160Sub80.mrc --subunit_diameter 80 --cyl_outer_diameter 160 --angpix 1 --rise 21.8 --twist 29.41 --boxdim 450 --sym_Cn 6 --sphere_percentage 0.9 --width 5 Still got the same error message. From what I understand of the wiki page: relion_helix_toolbox [cylinder] --help And I also tried: relion_helix_toolbox cylinder --help relion_helix_toolbox --help cylinder Still not working : ( Do you mind just showing me what's the correct command to get this "help" option work? Thanks a lot! Jing |
Original comment by Jing_basel (Bitbucket: Jing_basel, GitHub: Unknown): Hi Sjors, Thanks for your reply! I tried the command: Still got the same error message: Only one option can be specified at one time! valid_options = 128 Maybe it's a c++ compiling issue? |
Original comment by Shaoda He (Bitbucket: sh802, GitHub: sh802): Dear Jing,
// valid_options += (do_impose_helical_symmetry_fourier_space) ? (1) : (0); Best, Shaoda |
Original comment by Jing_basel (Bitbucket: Jing_basel, GitHub: Unknown): Hi Shaoda, Thanks a lot for your suggestion! //valid_options += (do_remove_segments_with_bad_psi) ? (1) : (0); I'm not familiar with c++, but it seems like some of those options are set to true by default? Jing |
Original comment by Shaoda He (Bitbucket: sh802, GitHub: sh802): Hi Jing. I am happy that you can use the options you want. But I think the bug still exists and everything seems fine with our compiler. Could you do me a favour? Could you replace "src/apps/helix_toolbox.cpp" with this one: ftp://ftp.mrc-lmb.cam.ac.uk/pub/she/2016-10-13/helix_toolbox.cpp , recompile and type in the command "relion_helix_toolbox --help" ? There should be an additional line before the error message: Best, Shaoda |
Original comment by Jing_basel (Bitbucket: Jing_basel, GitHub: Unknown): Hi Shaoda, /scicore/home/basler/wangji/src/relion2.0beta/src/apps/helix_toolbox.cpp: In member function ‘void helix_bilder_parameters::run()’:
|
Original comment by Jing_basel (Bitbucket: Jing_basel, GitHub: Unknown): Alright, the culprit is option "do_remove_mics_with_bad_ctf". I added a line to print the value of valid_options after each valid_option += xxxxx and here's the output: [wangji@login11 contracted_2aa]$ relion_helix_toolbox --help valid_options_psi = 0 Once I commented out the ctf thingy, which is actually not listed in the helix_toolbox functions: Everything else seems fine: Jing |
Original comment by Shaoda He (Bitbucket: sh802, GitHub: sh802): Hi Jing, I am sorry for replying late. Good spot! That's my fault. Thanks for helping us solve the problem. Could you try compiling this file: ftp://ftp.mrc-lmb.cam.ac.uk/pub/she/2016-10-17/helix_toolbox.cpp ? I think everything should be fine now. Best, Shaoda |
Original comment by Jing_basel (Bitbucket: Jing_basel, GitHub: Unknown): Hi Shaoda, Thank you very much for fixing the code! The new helix_toolbox.cpp file works quite well except just one more problem that I found: the "remove_bad_ctf" doesn't take the arg of "--EPA_lowest_res", alone or after "--ctf_fom_min". I had a look through the input data.star file and the "_rlnCtfMaxResolution #21" ranges from 3-5 A. However no particle was removed by the relion_helix_tool --EPA_lowest_res criteria. Here's the result of a few tests: [wangji@login11 fourth_round]$ relion_helix_toolbox --remove_bad_ctf --i run_data.star --o out.star --ctf_fom_min 0.1 [wangji@login11 fourth_round]$ relion_helix_toolbox --remove_bad_ctf --i run_data.star --o out.star --EPA_lowest_res 5 [wangji@login11 fourth_round]$ relion_helix_toolbox --remove_bad_ctf --i run_data.star --o out.star --EPA_lowest_res 4 [wangji@login11 fourth_round]$ relion_helix_toolbox --remove_bad_ctf --i run_data.star --o out.star --EPA_lowest_res 3 [wangji@login11 fourth_round]$ relion_helix_toolbox --remove_bad_ctf --i run_data.star --o out.star --EPA_lowest_res 2 [wangji@login11 fourth_round]$ relion_helix_toolbox --remove_bad_ctf --i run_data.star --o [wangji@login11 fourth_round]$ relion_helix_toolbox --remove_bad_ctf --i run_data.star --o out.star --ctf_fom_min 0.1 --o out.star --EPA_lowest_res 4
|
Original comment by Shaoda He (Bitbucket: sh802, GitHub: sh802): Hi Jing. Thanks for testing this option. I haven't found the bug by myself because I always use "subset selection" (relion_display) to get rid of the bad micrographs after CTF correction. We encourage users to work inside pipelines whenever possible. Anyway, please copy the file: ftp://ftp.mrc-lmb.cam.ac.uk/pub/she/src/helix.cpp to your "src" directory and try again. Best, |
Original comment by Jing_basel (Bitbucket: Jing_basel, GitHub: Unknown): Hi Shaodao, Thanks for the update! The new helix.cpp caused an installation error: [ 48%] Building CXX object src/apps/CMakeFiles/relion_lib.dir//pipeliner.cpp.o Is it an easy fix? Otherwise I'll just use the subset selection to discard micrographs. Cheers, Jing |
Original comment by Shaoda He (Bitbucket: sh802, GitHub: sh802): Hi Jing, I am sorry that we have very different versions of RELION 2.0. Although this problem can be fixed by replacing "metadata_label.h", I don't want you to do that because it risks messing up other parts of the program. Could you try the old "helix.cpp" file with modified "void excludeLowCTFCCMicrographs()" function? void excludeLowCTFCCMicrographs(
} Best, Shaoda |
Original comment by Jing_basel (Bitbucket: Jing_basel, GitHub: Unknown): Hi Shaoda, [wangji@login11 fourth_round]$ relion_helix_toolbox --remove_bad_ctf --i run_data.star --o test.star --EPA_lowest_res 3 Sorry about this back and forth, maybe I'll just wait for the new version of Relion2. Cheers, Jing |
Original comment by Shaoda He (Bitbucket: sh802, GitHub: sh802): Hi Jing, I am sorry that with so many versions of RELION 2.0 around I don't know exactly what you have. Could you send the input file ("run_data.star"???, you may compress the file if it is too large) and the commands to me by email (sh802@cam.ac.uk or she@mrc-lmb.cam.ac.uk) ? If the option in my version works, I suggest that you just wait for the next public release of RELION 2.0 and select subsets with the GUI for the moment :) Shaoda |
Original comment by Shaoda He (Bitbucket: sh802, GitHub: sh802): Hi Jing. Thanks for your email. I have performed tests on your STAR file and the results are stored in: ftp://ftp.mrc-lmb.cam.ac.uk/pub/she/helix/jing (commands and the screen output are copied to "notes" file). I think my version works well. And additional tests using my STAR file have also been done: ftp://ftp.mrc-lmb.cam.ac.uk/pub/she/helix/shaoda (Now I have added new parameters "--df_min" and "--df_max"). I guess we might have different set of EMDL labels. If you are still interested in what happened, could you send "src/metadata_label.h" to me by email? Or you may just wait for the next release of RELION :) I also suggest that micrographs (with their CTF information estimated), not extracted particles, should be selected with "helix_toolbox --remove_bad_ctf". And it is encouraged to select subsets in the GUI if you don't have complicated selection criteria. Best, Shaoda |
Hi, ERROR: unknown error in /programs/x86_64-linux/relion/2.1b1_cu8.0/src/gpu_utils/cuda_projector.cu at line 126 (error-code 30) |
validate tomogram_ids against tilt image files during import
Originally reported by: Jing_basel (Bitbucket: Jing_basel, GitHub: Unknown)
Dear all,
I would like to use the region_helix_tool box to get rid of some particles with huge out of plane tilt angle and soon I realised no matter what I type after the "region_helix_toolbox" command, it always prints out a list of options followed by the same error message:
Only one option can be specified at one time! valid_options = 127
File: /scicore/home/basler/wangji/src/relion2.0beta/src/apps/helix_toolbox.cpp line: 294
Is there something specific need to be done to compile the helix_tool box? There are plenty nice options of this tool box, it would be great to have it work.
Thanks a lot,
Jing
The text was updated successfully, but these errors were encountered: