From: Esben J. B. <esb...@ro...> - 2014-01-29 14:15:38
|
Hi All-knowing Pymol'ers. I wonder if anyone had a couple of tips/ pointers / script snippets. I have a list of coordinates on a regular grid with associated floating point values from a calculation that I want to visualise in PyMOL. Is it nescessary to write this to an Xplor/CNS/CCP4 file format and subsequently load it, or can it be loaded directly into a map object via some fancy python code? Best Regards Esben Jannik Bjerrum |
From: Jason V. <jas...@sc...> - 2014-01-30 21:24:35
|
Hi Esben, I have a list of coordinates on a regular grid with associated floating > point values from a calculation that I want to visualise in PyMOL. Is it > nescessary to write this to an Xplor/CNS/CCP4 file format and subsequently > load it, or can it be loaded directly into a map object via some fancy > python code? It depends on how you want to visualize the data. If dots/points or other discrete representations are okay then you can use a XYZ or a modified structure file. If you want the benefit of a field (volume, gradient, isosurface) then I suggest the map. If you already have the data in a well-formatted brick, check out the asCCP4 function in http://www.pymolwiki.org/index.php/Tiff2ccp4. That writes the CCP4 header and data (of a TIFF file, but you can modify that for your needs). If you want, send me a copy of your data. Or, I also have a very poorly written script that will convert raw volumes of data to CCP4 maps, similar to the aforementioned script. Email me personally and I can send you a copy if you prefer (but it's really ugly code). Cheers, -- Jason On Wed, Jan 29, 2014 at 8:02 AM, Esben Jannik Bjerrum < esb...@ro...> wrote: > > Hi All-knowing Pymol'ers. > > I wonder if anyone had a couple of tips/ pointers / script snippets. > > I have a list of coordinates on a regular grid with associated floating point values from a calculation that I want to visualise in PyMOL. Is it nescessary to write this to an Xplor/CNS/CCP4 file format and subsequently load it, or can it be loaded directly into a map object via some fancy python code? > > Best Regards > Esben Jannik Bjerrum > > ------------------------------------------------------------------------------ > WatchGuard Dimension instantly turns raw network data into actionable > security intelligence. It gives you real-time visual feedback on key > security issues and trends. Skip the complicated setup - simply import > a virtual appliance and go from zero to informed in seconds. > http://pubads.g.doubleclick.net/gampad/clk?id=123612991&iu=/4140/ostg.clktrk > _______________________________________________ > PyMOL-users mailing list (PyM...@li...) > Info Page: https://lists.sourceforge.net/lists/listinfo/pymol-users > Archives: http://www.mail-archive.com/pym...@li... -- Jason Vertrees, PhD Director of Core Modeling Products Schrödinger, Inc. (e) Jas...@sc... (o) +1 (603) 374-7120 |
From: Thomas H. <tho...@sc...> - 2014-01-30 22:45:03
|
Hi Esben, I just pushed a change to the open-source SVN repo which makes loading maps with a script a little easier. If you have your map in a numpy array, you can do: yourdata = ... # array with shape (a,b,c) yourgridspacing = (0.2, 0.2, 0.2) # grid spacing from chempy.brick import Brick b = Brick.from_numpy(yourdata, yourgridspacing) cmd.load_brick(b, 'yourmap') Hope that helps. Cheers, Thomas On 30 Jan 2014, at 16:24, Jason Vertrees <jas...@sc...> wrote: > Hi Esben, > > I have a list of coordinates on a regular grid with associated floating point values from a calculation that I want to visualise in PyMOL. Is it nescessary to write this to an Xplor/CNS/CCP4 file format and subsequently load it, or can it be loaded directly into a map object via some fancy python code? > > It depends on how you want to visualize the data. If dots/points or other discrete representations are okay then you can use a XYZ or a modified structure file. If you want the benefit of a field (volume, gradient, isosurface) then I suggest the map. > > If you already have the data in a well-formatted brick, check out the asCCP4 function in http://www.pymolwiki.org/index.php/Tiff2ccp4. That writes the CCP4 header and data (of a TIFF file, but you can modify that for your needs). > > If you want, send me a copy of your data. Or, I also have a very poorly written script that will convert raw volumes of data to CCP4 maps, similar to the aforementioned script. Email me personally and I can send you a copy if you prefer (but it's really ugly code). > > Cheers, > > -- Jason > > > On Wed, Jan 29, 2014 at 8:02 AM, Esben Jannik Bjerrum <esb...@ro...> wrote: > > > > Hi All-knowing Pymol'ers. > > > > I wonder if anyone had a couple of tips/ pointers / script snippets. > > > > I have a list of coordinates on a regular grid with associated floating point values from a calculation that I want to visualise in PyMOL. Is it nescessary to write this to an Xplor/CNS/CCP4 file format and subsequently load it, or can it be loaded directly into a map object via some fancy python code? > > > > Best Regards > > Esben Jannik Bjerrum -- Thomas Holder PyMOL Developer Schrödinger, Inc. |
From: Esben J. B. <esb...@ro...> - 2014-02-02 17:40:34
|
Hi Jason and Thomas, Thx for you suggestions, Maps are definitely the way to go, and I already have the data it in a numpy array. I also found a nifty example from brick01.py in the examples directory, but to get it to render in linux, numpy needs to be compiled into PyMOL, does this also hold for Your solution Thomas? Loading directly from Python is the way to go in the future. In the meantime I solved my problem by writing a small script that dumps the data in dx file format, which can subsequently be loaded into PyMOL. #!/usr/bin/python # Esben Jannik Bjerrum 2014 #Write DX format import math nx,ny,nz = 100,100,100 # Number of grid point in X,y,z dim xmin,ymin,zmin = 0,0,0 # coordinate of grid lowest corner. hx = hy = hz = 0.2 #Spacing between grid points f = open("tester.dx","w") f.write("#Comment\n") f.write("object 1 class gridpositions counts %s %s %s\n"% (nx,ny,nz)) f.write("origin %s %s %s\n"% (xmin,ymin,zmin)) f.write("delta %s 0.0 0.0\ndelta 0.0 %s 0.0\ndelta 0.0 0.0 %s\n"%(hx,hy,hz)) f.write("object 2 class gridconnections counts %s %s %s\n"%(nx,ny,nz)) f.write("object 3 class array type double rank 0 times %s\n"%(nx*ny*nz)) #The data with Z increasing first, then Y, then X. i = 0 for x in range(nx): for y in range(ny): for z in range(nz): #Generation of sample data mx=1.5*(x-nx/2.) my=2*(y-nx/2.) mz=3*(z-nz/2.) #Write values of point number (x,y,z) f.write("%e "%(1.-2./math.exp(math.sqrt(mx*mx+my*my+mz*mz)/10.0))) #f.write("(%s,%s,%s)"%(x,y,z)) if i == 2: f.write("\n") i = 0 else: i += 1 f.write("\n") f.close() Best Regards Esben On Thursday, January 30, 2014 11:45 PM, Thomas Holder <tho...@sc...> wrote: Hi Esben, I just pushed a change to the open-source SVN repo which makes loading maps with a script a little easier. If you have your map in a numpy array, you can do: yourdata = ... # array with shape (a,b,c) yourgridspacing = (0.2, 0.2, 0.2) # grid spacing from chempy.brick import Brick b = Brick.from_numpy(yourdata, yourgridspacing) cmd.load_brick(b, 'yourmap') Hope that helps. Cheers, Thomas On 30 Jan 2014, at 16:24, Jason Vertrees <jas...@sc...> wrote: > Hi Esben, > > I have a list of coordinates on a regular grid with associated floating point values from a calculation that I want to visualise in PyMOL. Is it nescessary to write this to an Xplor/CNS/CCP4 file format and subsequently load it, or can it be loaded directly into a map object via some fancy python code? > > It depends on how you want to visualize the data. If dots/points or other discrete representations are okay then you can use a XYZ or a modified structure file. If you want the benefit of a field (volume, gradient, isosurface) then I suggest the map. > > If you already have the data in a well-formatted brick, check out the asCCP4 function in http://www.pymolwiki.org/index.php/Tiff2ccp4. That writes the CCP4 header and data (of a TIFF file, but you can modify that for your needs). > > If you want, send me a copy of your data. Or, I also have a very poorly written script that will convert raw volumes of data to CCP4 maps, similar to the aforementioned script. Email me personally and I can send you a copy if you prefer (but it's really ugly code). > > Cheers, > > -- Jason > > > On Wed, Jan 29, 2014 at 8:02 AM, Esben Jannik Bjerrum <esb...@ro...> wrote: > > > > Hi All-knowing Pymol'ers. > > > > I wonder if anyone had a couple of tips/ pointers / script snippets. > > > > I have a list of coordinates on a regular grid with associated floating point values from a calculation that I want to visualise in PyMOL. Is it nescessary to write this to an Xplor/CNS/CCP4 file format and subsequently load it, or can it be loaded directly into a map object via some fancy python code? > > > > Best Regards > > Esben Jannik Bjerrum -- Thomas Holder PyMOL Developer Schrödinger, Inc. |
From: Thomas H. <tho...@sc...> - 2014-02-03 02:01:38
|
Hi Esben, On 02 Feb 2014, at 12:27, Esben Jannik Bjerrum <esb...@ro...> wrote: > ... I also found a nifty example from brick01.py in the examples directory, but to get it to render in linux, numpy needs to be compiled into PyMOL, does this also hold for Your solution Thomas? yes, PyMOL needs to be compiled with numpy for chempy.brick support. Incentive PyMOL ships with numpy support, and open-source PyMOL has numpy support if it can "import numpy" during compilation. Cheers, Thomas -- Thomas Holder PyMOL Developer Schrödinger, Inc. |
From: Gianluca S. <gia...@ib...> - 2014-06-12 14:58:41
|
Hi everyone, I was following this instructions to load a nparray as a map, but I get the error: class Brick has no attribute 'from_numpy' Is it a new feature in chempy, or am I missing something? Cheers, Gian On 2/3/14 3:01 AM, Thomas Holder wrote: > Hi Esben, > > On 02 Feb 2014, at 12:27, Esben Jannik Bjerrum <esb...@ro...> wrote: >> ... I also found a nifty example from brick01.py in the examples directory, but to get it to render in linux, numpy needs to be compiled into PyMOL, does this also hold for Your solution Thomas? > > yes, PyMOL needs to be compiled with numpy for chempy.brick support. Incentive PyMOL ships with numpy support, and open-source PyMOL has numpy support if it can "import numpy" during compilation. > > Cheers, > Thomas > -- Gianluca Santoni, Dynamop Group Institut de Biologie Structurale 6 rue Jules Horowitz 38027 Grenoble Cedex 1 France _________________________________________________________ Please avoid sending me Word or PowerPoint attachments. See http://www.gnu.org/philosophy/no-word-attachments.html |
From: Thomas H. <tho...@sc...> - 2014-06-12 16:05:39
|
Hi Gian, yes, the from_numpy method is new in PyMOL 1.7. But the underlying import functionality is not new, so you could basically copy the code from that method and create a valid brick instance yourself. Have a look at: https://sourceforge.net/p/pymol/code/HEAD/tree/trunk/pymol/modules/chempy/brick.py#l28 Cheers, Thomas On 12 Jun 2014, at 16:40, Gianluca Santoni <gia...@ib...> wrote: > Hi everyone, > I was following this instructions to load a nparray as a map, but I get > the error: > class Brick has no attribute 'from_numpy' > > Is it a new feature in chempy, or am I missing something? > Cheers, > Gian > > On 2/3/14 3:01 AM, Thomas Holder wrote: >> Hi Esben, >> >> On 02 Feb 2014, at 12:27, Esben Jannik Bjerrum <esb...@ro...> wrote: >>> ... I also found a nifty example from brick01.py in the examples directory, but to get it to render in linux, numpy needs to be compiled into PyMOL, does this also hold for Your solution Thomas? >> >> yes, PyMOL needs to be compiled with numpy for chempy.brick support. Incentive PyMOL ships with numpy support, and open-source PyMOL has numpy support if it can "import numpy" during compilation. >> >> Cheers, >> Thomas > > -- > Gianluca Santoni, > Dynamop Group > Institut de Biologie Structurale > 6 rue Jules Horowitz > 38027 Grenoble Cedex 1 > France -- Thomas Holder PyMOL Developer Schrödinger, Inc. |
From: Shiven S. <shi...@gm...> - 2014-06-12 21:49:14
|
Hi Martin, I guess you can do this using modified GLSL shaders … Two examples I posted on the PyMOL wiki a while ago: http://www.pymolwiki.org:8888/index.php/GLSL_Shaders Best, Shiven On Jun 12, 2014, at 5:10 PM, Martin Hediger <ma...@bl...> wrote: > Hi all > > I wonder if it's possible to have something like ray_trace_mode, 2 (i.e. > black and white with outlines) but in the interactive view. Also, is it > possible to have part of the structure in color (say some active site > amino acids) and the rest in black and white? > > Thanks for help > Martin > > ------------------------------------------------------------------------------ > HPCC Systems Open Source Big Data Platform from LexisNexis Risk Solutions > Find What Matters Most in Your Big Data with HPCC Systems > Open Source. Fast. Scalable. Simple. Ideal for Dirty Data. > Leverages Graph Analysis for Fast Processing & Easy Data Exploration > http://p.sf.net/sfu/hpccsystems > _______________________________________________ > PyMOL-users mailing list (PyM...@li...) > Info Page: https://lists.sourceforge.net/lists/listinfo/pymol-users > Archives: http://www.mail-archive.com/pym...@li... |
From: Martin H. <ma...@bl...> - 2014-06-12 22:10:24
|
Impressive, how did you figure this out? Martin On 12.06.14 23:18, Shiven Shandilya wrote: > Hi Martin, > > I guess you can do this using modified GLSL shaders … > > Two examples I posted on the PyMOL wiki a while ago: > http://www.pymolwiki.org:8888/index.php/GLSL_Shaders > > Best, > Shiven > > > On Jun 12, 2014, at 5:10 PM, Martin Hediger <ma...@bl... > <mailto:ma...@bl...>> wrote: > >> Hi all >> >> I wonder if it's possible to have something like ray_trace_mode, 2 (i.e. >> black and white with outlines) but in the interactive view. Also, is it >> possible to have part of the structure in color (say some active site >> amino acids) and the rest in black and white? >> >> Thanks for help >> Martin >> >> ------------------------------------------------------------------------------ >> HPCC Systems Open Source Big Data Platform from LexisNexis Risk Solutions >> Find What Matters Most in Your Big Data with HPCC Systems >> Open Source. Fast. Scalable. Simple. Ideal for Dirty Data. >> Leverages Graph Analysis for Fast Processing & Easy Data Exploration >> http://p.sf.net/sfu/hpccsystems >> _______________________________________________ >> PyMOL-users mailing list (PyM...@li... >> <mailto:PyM...@li...>) >> Info Page: https://lists.sourceforge.net/lists/listinfo/pymol-users >> Archives: http://www.mail-archive.com/pym...@li... > -- Dr. Martin R. Hediger Schaffhauserstrasse 104 8057 Zürich |
From: Martin H. <ma...@bl...> - 2014-06-12 21:10:57
|
Hi all I wonder if it's possible to have something like ray_trace_mode, 2 (i.e. black and white with outlines) but in the interactive view. Also, is it possible to have part of the structure in color (say some active site amino acids) and the rest in black and white? Thanks for help Martin |
From: Marcelo M. <MM...@na...> - 2014-07-01 06:37:23
|
Hi Everyone, If one has a functional Nvidia Nvision 3D setup, can you export a video that retains 3D? How is this done? Thanks in advance! Marcelo |