<html><head><meta http-equiv="Content-Type" content="text/html; charset=us-ascii"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">Hi Rakesh,<div class=""><br class=""></div><div class=""> Here are the commands you were running (I would usually write those to a file with suffix ".cxc" and have ChimeraX open that file, makes it easier to debug). The error is because the volume surface is not created until it is displayed. And in a script nothing will be displayed until the end unless you put it a command to display it. So add a "wait 1" command before the "color zone #3..." command. The "wait 1" command means draw one frame and that will cause the surface to be computed.</div><div class=""><br class=""></div><div class=""> Maybe the "wait 1" is not needed. It looks like your script on line 8 is showing the volume #1 as a surface when really you want show volume #3 as a surface. Maybe if you change that to #3 it will make the surface without rendering.</div><div class=""><div class=""><br class=""></div><div class=""><div class="">open ../emmaps/emd_5584.mrc</div><div class="">volume #1 originIndex 0,0,0</div><div class="">save ../emmaps/emd_5584_cen.mrc</div><div class="">open 3j31_A.pdb</div><div class="">move 701.4400024414062,701.4400634765625,701.4400024414062 model #2</div><div class="">save 3j31_A_cen.pdb model #2</div><div class="">open ../emmaps/emd_5584_cen.mrc</div><div class="">volume #1 region all style surface level 0.05</div><div class="">open 3j31_A_cen.pdb</div><div class="">color #4/A blue</div><div class="">color zone #3 near #4/A distance 5.0 sharpEdges true</div><div class="">volume splitbyzone #3</div><div class="">save ../emmaps/emd_5584_seg.mrc #5.2</div><div class=""><br class=""></div><div class=""> Tom</div><div class=""><br class=""></div><div><br class=""><blockquote type="cite" class=""><div class="">On Feb 21, 2022, at 8:23 PM, Rakesh Ramachandran via ChimeraX-users <<a href="mailto:chimerax-users@cgl.ucsf.edu" class="">chimerax-users@cgl.ucsf.edu</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div dir="ltr" class="">Hi Tom,<br class=""><div class=""><br class=""></div><div class=""> I tried your suggestion however I still get this error message.</div><div class=""><br class=""></div><div class="">ChimeraX --nogui --exit --silent --cmd 'open ../emmaps/emd_5584.mrc; volume #1 originIndex 0,0,0; save ../emmaps/emd_5584_cen.mrc; open 3j31_A.pdb; move 701.4400024414062,701.4400634765625,701.4400024414062 model #2; save 3j31_A_cen.pdb model #2; open ../emmaps/emd_5584_cen.mrc; volume #1 region all style surface level 0.05; open 3j31_A_cen.pdb; color #4/A blue; color zone #3 near #4/A distance 5.0 sharpEdges true; volume splitbyzone #3; save ../emmaps/emd_5584_seg.mrc #5.2'<br class="">Traceback (most recent call last):<br class=""> File "/home/rakesh/WORK/Softwares/Molecular_Visualizers/chimerax/lib/python3.9/site-packages/ChimeraX_main.py", line 713, in init<br class=""> run(sess, cmd)<br class=""> File "/home/rakesh/WORK/Softwares/Molecular_Visualizers/chimerax/lib/python3.9/site-packages/chimerax/core/commands/run.py", line 38, in run<br class=""> results = command.run(text, log=log, return_json=return_json)<br class=""> File "/home/rakesh/WORK/Softwares/Molecular_Visualizers/chimerax/lib/python3.9/site-packages/chimerax/core/commands/cli.py", line 2897, in run<br class=""> result = ci.function(session, **kw_args)<br class=""> File "/home/rakesh/WORK/Softwares/Molecular_Visualizers/chimerax/lib/python3.9/site-packages/chimerax/std_commands/color.py", line 1319, in color_zone<br class=""> raise UserError('color zone: No surfaces specified.')<br class="">chimerax.core.errors.UserError: color zone: No surfaces specified.<br class=""></div><div class=""><br class=""></div><div class="">What I am trying to do here is to change the originIndex to 0,0,0 for the original map, then move the original fitted pdb structure according to the shifted map, color the fitted pdb and use volume splitbyzone to segment the map. Let me know if I am doing something wrong here, the whole chimera command is printed above.</div><div class=""><br class=""></div><div class="">Also when I segment the map and save it, the box size is equal to the original map, is there a way I can resize it and keep it cubic, probably with the box size equal to the longest dimension of the segmented map. </div><br class=""><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Feb 18, 2022 at 8:20 PM Tom Goddard <<a href="mailto:goddard@sonic.net" target="_blank" class="">goddard@sonic.net</a>> wrote:<br class=""></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi Rakesh,<br class="">
<br class="">
ChimeraX displays just a single plane when you open a map file larger than 256 Mvoxels -- often the situation for electron microscopy tomograms, for example a 2048 by 2048 by 100 map (400 Mvoxels), but also for large electron microscopy single particle reconstructions, e.g. 800 x 800 x 800 (512 Mvoxels). Maps larger than 16 Mvoxels will not be shown at full resolution (step 1), instead the step size will be chosen to only display at most 16 Mvoxels when the map is first opened, again to improve load speed for common interactive use cases.<br class="">
<br class="">
So if you are writing a script and need the full map shown always include a volume command to assure it is shown how you want, for example to show all planes,<br class="">
<br class="">
volume #1 region all style surface<br class="">
<br class="">
or to assure it is shown at full resolution<br class="">
<br class="">
volume #1 step 1 region all style surface<br class="">
<br class="">
and to make sure the threshold level is what you want you can specify it too<br class="">
<br class="">
volume #1 level 4.5 step 1 region all style surface<br class="">
<br class="">
There are volume command options that set the 256 MVoxel limit for planes (option name "voxelLimitForPlane") and 16 MVoxel full resolution limit (option name "voxelLimit") that you can use to change the default for the current session, but unfortunately ChimeraX currently does not offer a way to save new default values.<br class="">
<br class="">
These are all documented<br class="">
<br class="">
<a href="https://www.cgl.ucsf.edu/chimerax/docs/user/commands/volume.html#sampling" rel="noreferrer" target="_blank" class="">https://www.cgl.ucsf.edu/chimerax/docs/user/commands/volume.html#sampling</a><br class="">
<br class="">
Tom<br class="">
<br class="">
<br class="">
> On Feb 18, 2022, at 12:47 PM, Rakesh Ramachandran <<a href="mailto:rakeshramachandran10@gmail.com" target="_blank" class="">rakeshramachandran10@gmail.com</a>> wrote:<br class="">
> <br class="">
> Hi Tom,<br class="">
> <br class="">
> I am trying to segment volume using 'volume splitbyzone' in ChimeraX for a large number of maps. But what I observed is that some of the maps are displayed as planes instead of surfaces and this affects the volume splitbyzone command. <br class="">
> <br class="">
> I tried using 'volume style surface' but it doesn't seem to work unless I manually set it in the GUI. What option do I use to generate surfaces for map volumes automatically from the command line in such cases?<br class="">
> <br class="">
> -- <br class="">
> Regards<br class="">
> Rakesh Ramachandran<br class="">
> Postdoctoral Scholar<br class="">
> Department of Bioengineering and Therapeutic Sciences<br class="">
> University of California, San Francisco<br class="">
> Byers Hall Room 501<br class="">
> 1700 4th Street<br class="">
> San Francisco, CA 94143, USA<br class="">
<br class="">
</blockquote></div><br clear="all" class=""><div class=""><br class=""></div>-- <br class=""><div dir="ltr" class=""><div dir="ltr" class="">Regards<br class="">Rakesh Ramachandran<br class="">Postdoctoral Scholar<br class="">Department of Bioengineering and Therapeutic Sciences<br class="">University of California, San Francisco<br class="">Byers Hall Room 501<br class="">1700 4th Street<br class="">San Francisco, CA 94143, USA<br class="">Mob No: +1 415-489-8017<br class=""></div></div></div>
_______________________________________________<br class="">ChimeraX-users mailing list<br class=""><a href="mailto:ChimeraX-users@cgl.ucsf.edu" class="">ChimeraX-users@cgl.ucsf.edu</a><br class="">Manage subscription:<br class="">https://www.rbvi.ucsf.edu/mailman/listinfo/chimerax-users<br class=""></div></blockquote></div><br class=""></div></div></body></html>