ICM Manual v.3.9
by Ruben Abagyan,Eugene Raush and Max Totrov
Copyright © 2020, Molsoft LLC
Nov 14 2024

Contents
 
Introduction
Reference Guide
Command Line User's Guide
 ICM-shell
 ICM graphics
 Str.Analysis
 Sequence
 Molcart
 Pharmacophores
 Energy
 Molecules
 Animation
  Movie rotate view
  Mc movie
 Symmetry
 X-ray
 Plotting
 Docking/VLS
 Examples
 _chemSuper
 PROTAC Modeling
 Chemical Conformation Generator
 RIDGE
 RIDE
 MD
References
Glossary
 
Index
PrevICM Language Reference
Animation
Next

[ Movie rotate view | Mc movie ]

How to rotate and zoom in a script


First, if you intend to save images as movie frames from your script, start ICM with the -24 option (see above). This will increase the image quality.

The simplest script will use the View(v1,v2,r) function to interpolate between views with different rotation and zoom. Example:

 
 build string "ACDF" # tetrapeptide 
   # zoom out and create a small remote view 
 v1=View()  # the 1st view 
   # zoom in and rotate 
 v2=View()  # the last view 
 for i=1,100 
   set view View(v1,v2,i*0.01) # interpolation 
 endfor 
Now you can add a specific display ( cpk, skin or anything else) and add the write image command after set view .
The following script and macro make a .mov file for a specified molecule.
 
# makemovie -f qt -o tamox.mov -c qt_anim -r 15 `ls -tr ../f*.tif` 
 print "DO NOT FORGET icmt -24;  set view 10 10 640 480" 
 macro molkino i_tZoom (20) i_tRotate (20) l_preview (yes) R_startView R_endView s_Files ("/icmdata/movies/tamox/f") 
  GRAPHICS.ballQuality = 15 
  nfr=30 
  for i=0,i_tZoom*nfr 
    set view View(R_startView ,R_endView, 0.5-0.5*Cos(180.*i/Real(i_tZoom*nfr))) 
    if(!l_preview) write image s_Files + i 
  endfor 
  B = 0. 
  for i=0,i_tRotate*nfr 
    A = 360./Real(nfr*i_tRotate) * i 
    A = 180. - 180.*Cos(A/2.) 
    rotate view Rot({0. 1. 0.} A-B) 
    B = A 
    if(!l_preview) write image s_Files + (i_tZoom*nfr+ i) 
  endfor 
 endmacro 
#read object "tamoxifen" 
 read object "propecia.ob" 
 set window 10 10 640 480 
 GRAPHICS.ballRadius=0.4 
 display xstick 
 color a_//c* green 
 color a_//h* white 
 read rarray "v0" 
 read rarray "v1" 
 set view v0 

Saving a movie to a movie file (in addition to the ICM animation), see also write movie To export an extermal movie file you need to use the following set of commands


write movie "/tmp/forutube.avi"  # give a name to your movie file. Option heavy for high quality
# your script
write movie on # option smooth for anti-aliasing
# some actions
write movie off
# preparations for the next scene
write movie exit


How to make a molecular movie from a Monte Carlo trajectory


[ Image collection | Image frame by frame ]

We assume that you have a Monte Carlo trajectory saved as an ICM trajectory file. Displaying this trajectory interactively on your graphical screen is the simplest form of animation. You just load the object, the ICM trajectory file and run the display trajectory command. Adjust your window and view and the allowed types of representation. This method will not work in two situations.
  • First, skin and surface are not automatically recalculated and redisplayed upon the conformational change.
  • Second, if you would like to display computationally heavy forms of graphical representation the interactive trajectory may become too slow.
In these two situations you may resort to external tools generating movie from individual image files.
In this case, the process of making a molecular movie animating the trajectory can be split into two steps.
The first step is a preparation of a series of images stored as separate files. Watch the ICM trajectory by the display trajectory command and select a range of frames you are going to include to the movie. There are two ways how a series of the molecular graphics images can be saved as tif or rgb formatted files.

Image collection (type I)


The following molecular representations are allowed in this case: wire, cpk, ball, stick, xstick and ribbon . The image option of the ICM display trajectory command is required. An example:
 
 build       "alpha" 
 read trajectory  "alpha" 
 display ribbon 
 color ribbon a_/1:9   magenta 
 color ribbon a_/10:18 green 
 print "Adjust the view, press ENTER to continue" 
 pause             # Pause to adjust a view of the molecule 
                   # start of the collection of the image files 
                   # (note, tif files are to be saved) 
 nframe = Nof(frame) 
 print "Total number of frames = ", nframe 
                   # specify the range you need, for example: 
 ifrom = 10 
 ito   = 20 
                   # otherwise, for all frames: 
#  ifrom = 1 
#  ito   = nframe 
 
 display trajectory ifrom ito image  
Computationally more expensive skin and surface molecular representations (or a combination of any of these two with those mentioned above) require an ICM script, as in the following example.

Image collection (type II)


 
 build       "alpha" 
 read trajectory  "alpha" 
 display xstick a_//n,ca,c 
 set window 600 30 640 480  # good to convert it later into NTSC format 
# change   640 480  to appropriate dimensions for PAL 
 lineWidth = 3. 
 print "Adjust the view, press ENTER to continue" 
 pause 
 
 nframe = Nof(frame) 
 print "Total number of frames = ", nframe 
 ifrom = 10 
 ito   = 20 
 icount = 0 
 l_confirm = no 
 IMAGE.compress = yes 
 for iframe = ifrom, ito 
   load frame iframe 
   icount = icount + 1 
   imgname = s_tempDir + "alpha_" + String(icount-1) 
   display xstick a_//n,ca,c green 
   display surface a_/7:12/!n,ca,c a_/7:12/!n,ca,c magenta 
   write image imgname 
   undisplay surface       # should be removed and redisplayed at the next step 
 endfor 

Whatever protocol you make use of, a series of files should appear in the s_tempDir directory. List their names in a file, for example, alpha.lst:
 
/usr/tmp/alpha_0.tif 
/usr/tmp/alpha_1.tif 
/usr/tmp/alpha_2.tif 
/usr/tmp/alpha_3.tif 
/usr/tmp/alpha_4.tif 
/usr/tmp/alpha_5.tif 
/usr/tmp/alpha_6.tif 
/usr/tmp/alpha_7.tif 
/usr/tmp/alpha_8.tif 
/usr/tmp/alpha_9.tif 
(Note, image numbers are started from zero, and "_" (underscore) is inserted before the current frame number.) Now, you can put these image files into a movie file. We suggest to use the makemovie command for Silicon Graphics machines. The following is a sample UNIX-shell command line which should be entered to make a non-interlaced MPEG-formatted movie with the frame size 640 by 480 saved in the Apple QuickTime movie file alpha.qt:
 
makemovie -f qt -o alpha.qt -c qt_anim -s 640,480 -r 30 ` cat alpha.lst` 

Upon completion, use the movieplayer IRIX command to see the movie. UNIX man on-line help should allow you to learn more about these (and other) helpful commands. And let us give you several helpful hints:
  • While the collection of the image files is in progress, do not hide any portion or the whole area of the graphics window by opening another window which can cover the first one, and do not move the mouse cursor through the graphics window area. This procedure typically requires quite a time, so be ready to be off your terminal/workstation for the whole time required for the image collection to complete.
  • Do not forget to deactivate your screen saver before you start!
  • Make sure you have enough disk storage. All images are assumed to be stored in the directory specified by s_tempDir string variable. The easiest way to estimate the required disk storage is to store the content of a single image by the write image command, and multiply the size of the resulting tif or rgb file by the number of frames. You could do it without leaving your ICM-session. For example, if you have a molecular image in the graphics window, then:
     
     write image "testsize" 
     unix ls -l testsize.tif | awk '{print $5}' > testsize.rar 
     read rarray "testsize.rar" 
     sizeValue = testsize[1] 
     nFrame = 0 
     nFrame = Ask("Number of frames?", nFrame) 
     printf "Required storage for %d frames is %8.3f Mbytes\n", \ 
             nFrame, Real(nFrame)*sizeValue/1048576. 
     unix rm testsize.tif testsize.rar 
     delete sizeValue testsize nFrames 
    
(Note that in the above example the file is stored in the current directory, not in that defined by the s_tempDir variable).

Prev
Stereo reconstruction
Home
Up
Next
Symmetry

Copyright© 1989-2024, Molsoft,LLC - All Rights Reserved. Copyright© 1989-2024, Molsoft,LLC - All Rights Reserved. This document contains proprietary and confidential information of Molsoft, LLC. The content of this document may not be disclosed to third parties, copied or duplicated in any form, in whole or in part, without the prior written permission from Molsoft, LLC.