JiscMail Logo
Email discussion lists for the UK Education and Research communities

Help for CCPEM Archives


CCPEM Archives

CCPEM Archives


CCPEM@JISCMAIL.AC.UK


View:

Message:

[

First

|

Previous

|

Next

|

Last

]

By Topic:

[

First

|

Previous

|

Next

|

Last

]

By Author:

[

First

|

Previous

|

Next

|

Last

]

Font:

Proportional Font

LISTSERV Archives

LISTSERV Archives

CCPEM Home

CCPEM Home

CCPEM  November 2016

CCPEM November 2016

Options

Subscribe or Unsubscribe

Subscribe or Unsubscribe

Log In

Log In

Get Password

Get Password

Subject:

Re: Relion/2.0 GCTF wrapper errors

From:

Kai Zhang <[log in to unmask]>

Reply-To:

[log in to unmask]

Date:

Wed, 30 Nov 2016 21:56:07 -0000

Content-Type:

text/plain

Parts/Attachments:

Parts/Attachments

text/plain (371 lines)

Thanks for reporting this anyway.
I agree the script is not so user-friendly and will probably update it.
But if you find the rows of the the partiles are in random order, don't
worry abouth that.
Each particle still points to exactly the correct .mrcs file, it just
looks a bit ugly in such a random order.

Best wishes,
Kai


> Sorry, I jumped the gun. I fixed it.
>
>
> Zachary Berndsen, PhD
>
> The Scripps Research Institute
>
> Hazen Building (HL 104)
>
> Ward Lab
>
> ________________________________
> From: Collaborative Computational Project in Electron cryo-Microscopy
> <[log in to unmask]> on behalf of Zack Berndsen <[log in to unmask]>
> Sent: Wednesday, November 30, 2016 1:04:58 PM
> To: [log in to unmask]
> Subject: Re: [ccpem] Relion/2.0 GCTF wrapper errors
>
>
> Hi Kai,
>
>
> star_replaceUVA.com is not working properly for me. It's mixing up
> columns. Particles were extracted from the mics_all_gctf.star file and no
> processing was performed on the extracted particles. I am attempting to
> replace ctf info in the extracted particles.star file with the the local
> ctf info.
>
>
> Here is the the original particles.star file with the correct columns:
>
>
> data_
>
>
> loop_
>
> _rlnCoordinateX #1
>
> _rlnCoordinateY #2
>
> _rlnImageName #3
>
> _rlnMicrographName #4
>
> _rlnVoltage #5
>
> _rlnDefocusU #6
>
> _rlnDefocusV #7
>
> _rlnDefocusAngle #8
>
> _rlnSphericalAberration #9
>
> _rlnCtfBfactor #10
>
> _rlnCtfScalefactor #11
>
> _rlnPhaseShift #12
>
> _rlnAmplitudeContrast #13
>
> _rlnMagnification #14
>
> _rlnDetectorPixelSize #15
>
> _rlnCtfFigureOfMerit #16
>
>  2199.450878   800.450878
> [log in to unmask]
> mc2_5x5_Bft300_ref0/16aug05c_00003sq_00004hl_00002esup-a_movie_DW.mrc
> 300.000000 26826.996094 26031.109375    86.949692     2.700000
> 0.000000     1.000000     0.000000     0.100000 38168.000000     5.000000
>     0.132110
>
>
>
> Here is the the particles_localUVA.star file generated by
> star_replaceUVA.com:
>
>
>
> data_
>
>
> loop_
>
> _rlnCoordinateX #1
>
> _rlnCoordinateY #2
>
> _rlnImageName #3
>
> _rlnMicrographName #4
>
> _rlnVoltage #5
>
> _rlnDefocusU #6
>
> _rlnDefocusV #7
>
> _rlnDefocusAngle #8
>
> _rlnSphericalAberration #9
>
> _rlnCtfBfactor #10
>
> _rlnCtfScalefactor #11
>
> _rlnPhaseShift #12
>
> _rlnAmplitudeContrast #13
>
> _rlnMagnification #14
>
> _rlnDetectorPixelSize #15
>
> _rlnCtfFigureOfMerit #16
>
> 3045.048009  1922.652009
> [log in to unmask]
>
> mc2_5x5_Bft300_ref0/16aug05c_a_00012gr_00112sq_00011hl_00003esup-a_movie_DW.mrc
>  300.000000        2.700000  0.000000  1.000000  0.000000  0.100000
> 38168.000000  5.000000  0.187550
>
> The columns are clearly incorrect. DU and DV are suppose to follow voltage
> but CS does instead.
>
> Here is one of the ctf_local.star files with the correct columns:
>
>
>
> data_
>
>
> loop_
>
> _rlnMicrographName #1
>
> _rlnCoordinateX #2
>
> _rlnCoordinateY #3
>
> _rlnVoltage #4
>
> _rlnDefocusU #5
>
> _rlnDefocusV #6
>
> _rlnDefocusAngle #7
>
> _rlnSphericalAberration #8
>
> _rlnDetectorPixelSize #9
>
> _rlnCtfFigureOfMerit #10
>
> _rlnMagnification #11
>
> _rlnAmplitudeContrast #12
>
> _rlnAutopickFigureOfMerit #13
>
> mc2_5x5_Bft300_ref0/16aug05c_00003sq_00010hl_00002esup-a_movie.mrc
> 1294.000000  2701.000000   300.000000 24318.199219 23865.933594
> 74.219482     2.700000    14.000000     0.065700 106870.234375
> 0.100000     0.000000
>
>
>
>
> I can edit my star_replaceUVA.com script to reference the correct columns
> if that is the problem, you'll just have to point me to the correct lines.
>
>
> Thanks,
>
>
> Zachary Berndsen, PhD
>
> The Scripps Research Institute
>
> Hazen Building (HL 104)
>
> Ward Lab
>
> ________________________________
> From: Zack Berndsen
> Sent: Tuesday, November 29, 2016 3:16:58 PM
> To: [log in to unmask]; Zack Berndsen
> Subject: Re: Relion/2.0 GCTF wrapper errors
>
>
> Kai,
>
>
> I attempted a local ctf run with a set of reentered
> coordinate_extract.star files from a previous relion/2.0 processing run
> and it seems to have worked properly. I will read up on the _UVA.com
> script but it seems straight forward enough. Thanks a lot for the help.
>
>
> Best,
>
>
> Zachary Berndsen, PhD
>
> The Scripps Research Institute
>
> Hazen Building (HL 104)
>
> Ward Lab
>
> ________________________________
> From: Collaborative Computational Project in Electron cryo-Microscopy
> <[log in to unmask]> on behalf of Zack Berndsen <[log in to unmask]>
> Sent: Monday, November 28, 2016 1:17:09 PM
> To: [log in to unmask]
> Subject: [ccpem] Relion/2.0 GCTF wrapper errors
>
>
> Hi,
>
>
> This is a two part question. First, I am encountering a strange error in
> the final steps of GCTF execution through relion/2.0. I have run the
> wrapper successfully several times without the error and the solution is
> not obvious to me. Here is the error:
>
>
>
> ...
>
> *
> MotionCor2_5x5_Bft100_ref0/15feb27c_1_00018gr_00004sq_00017hl_00004ed-a.mrc
>
> ...
>
> Estimating CTF parameters using Kai Zhang's Gctf ...
>
> 5.60/5.60 hrs
> ............................................................~~(,_,">
>
> +++ RELION: command line arguments (with defaults for optional ones
> between parantheses) +++
>
> ...
>
> ERROR: Unexpected number of words on Final values line in
> CtfFind/mc2_5x5_ref0_rn2/MotionCor2_5x5_Bft100_ref0/15feb27c_1_00003gr_00021sq_00026hl_00004ed-a_gctf.log
>
> File: /tmp/relion2-beta/src/ctffind_runner.cpp line: 659
>
>
> Final few lines of the log file in question (I don't see anything
> obviously wrong):
>
>
> **************************************   LAST CYCLE
> ************************************************************ *
>
>
>    Defocus_U   Defocus_V       Angle         CCC
>
>     21388.87    20984.57       87.63    0.120526  Final Values
>
>
> Resolution limit estimated by EPA: RES_LIMIT 2.827
>
> Estimated Bfactor: B_FACTOR  21.37
>
> Processing done successfully.
>
> 4.59       87.75    0.120526  Final Values
>
>
> Resolution limit estimated by EPA: RES_LIMIT 2.827
>
> Estimated Bfactor: B_FACTOR  21.37
>
> Processing done successfully.
>
>
>
>
> Has anyone encountered this before? If so do you know what is causing it
> and how to fix it?
>
>
> Second, I've noticed in the log files that the XMAG and DStep are not
> being passed to the executable properly (correct values were given in the
> command). I know the wrapper is suppose to input the microscopy options
> from the gui into the command and ignore gctf defaults but the fact that
> they appear incorrect in the log file makes me wonder if they are in fact
> being used the in calculations. The apix is correct however and I assume
> this is the most important option (does it overwrite the others?). Also,
> my previous reconstructions where GCTF was used are great and I can't see
> any issues with them so I am assuming I should just ignore these
> inconsistencies, but it would be nice to know for sure. Are these values
> used in the calculations? See section of log file below:
>
>
>
> ...
>
> GCTF , v0.50 2014-11-05
>
> All parameters to be used:
>
> General options:
>
>         --apix               1.310
>
>         --dstep              14.000
>
>         --kv                 300.00
>
>         --cs                 2.70
>
>         --ac                 0.10
>
> ...
>
> ******************************** CTFFIND3 parameter format for RELION:
> *********************************************
>
>
> NOTE: NOTHING related to CTFFIND3. Just use it for RELION.
>
>
>  Input image file name
>
> CtfFind/mc2_5x5_ref0_rn2/MotionCor2_5x5_Bft100_ref0/15feb27c_1_00003gr_00021sq_00026hl_00004ed-a.mrc
>
>
>  Output diagnostic file name
>
> CtfFind/mc2_5x5_ref0_rn2/MotionCor2_5x5_Bft100_ref0/15feb27c_1_00003gr_00021sq_00026hl_00004ed-a.ctf
>
>
>  CS[mm], HT[kV], AmpCnst, XMAG, DStep[um]
>
>   2.7    300.0    0.10  106870.2   14.000
>
>
> Thank you!
>
>
>
> Zachary Berndsen, PhD
>
> The Scripps Research Institute
>
> Hazen Building (HL 104)
>
> Ward Lab
>

Top of Message | Previous Page | Permalink

JiscMail Tools


RSS Feeds and Sharing


Advanced Options


Archives

May 2024
April 2024
March 2024
February 2024
January 2024
December 2023
November 2023
October 2023
September 2023
August 2023
July 2023
June 2023
May 2023
April 2023
March 2023
February 2023
January 2023
December 2022
November 2022
October 2022
September 2022
August 2022
July 2022
June 2022
May 2022
April 2022
March 2022
February 2022
January 2022
December 2021
November 2021
October 2021
September 2021
August 2021
July 2021
June 2021
May 2021
April 2021
March 2021
February 2021
January 2021
December 2020
November 2020
October 2020
September 2020
August 2020
July 2020
June 2020
May 2020
April 2020
March 2020
February 2020
January 2020
December 2019
November 2019
October 2019
September 2019
August 2019
July 2019
June 2019
May 2019
April 2019
March 2019
February 2019
January 2019
December 2018
November 2018
October 2018
September 2018
August 2018
July 2018
June 2018
May 2018
April 2018
March 2018
February 2018
January 2018
December 2017
November 2017
October 2017
September 2017
August 2017
July 2017
June 2017
May 2017
April 2017
March 2017
February 2017
January 2017
December 2016
November 2016
October 2016
September 2016
August 2016
July 2016
June 2016
May 2016
April 2016
March 2016
February 2016
January 2016
December 2015
November 2015
October 2015
September 2015
August 2015
July 2015
June 2015
May 2015
April 2015
March 2015
February 2015
January 2015
December 2014
November 2014
October 2014
September 2014
August 2014
July 2014
June 2014
May 2014
April 2014
March 2014
February 2014
January 2014
December 2013
November 2013
October 2013
September 2013
August 2013
July 2013
June 2013
May 2013
April 2013
March 2013
February 2013


JiscMail is a Jisc service.

View our service policies at https://www.jiscmail.ac.uk/policyandsecurity/ and Jisc's privacy policy at https://www.jisc.ac.uk/website/privacy-notice

For help and support help@jisc.ac.uk

Secured by F-Secure Anti-Virus CataList Email List Search Powered by the LISTSERV Email List Manager