Outline

From PyMOLWiki
Revision as of 13:08, 28 July 2014 by Jaredsampson (talk | contribs) (→‎Example: added "Automation" heading. Added `outfile` variable in imagemagick script, and adjusted filenames for consistency.)
Jump to navigation Jump to search

Outline

There is currently (as of PyMOL v.1.7) no way to outline a selection directly within PyMOL. However, it can be accomplished indirectly with a composite image.

Relevant Settings

Example


The following script was used to generate the first two images from PyMOL.

bg_color white

fetch 3cyt, async=0

as surface
color marine

select outline, br. chain I and resi 85 around 2

set_view (\
     0.061975956,   -0.950684488,    0.303902954,\
     0.703773856,    0.257531703,    0.662103057,\
    -0.707715809,    0.172843516,    0.685028315,\
     0.000000000,    0.000000000, -152.244812012,\
    25.274658203,    8.288025856,    9.110867500,\
    51.143974304,  253.345642090,  -20.000000000 )

png base.png, ray=1

hide everything
as surface, outline
set ray_trace_mode, 2
set ray_trace_color, yellow
set ray_opaque_background, 0

png overlay.png, ray=1

In an image editing program, position overlay.png in a layer exactly overlapping base.png. The extra outline lines can also (optionally) be cleaned up by carefully erasing them using the graphics program's eraser tool.

Automation

If you need to have a fully automated solution, you can use ImageMagick or GraphicMagick to do this. If you use ImageMagick in a Unix box, you can type this command on the terminal:

composite -gravity center overlay.png base.png composite.png

See this page on ImageMagick for more information. If you need to use it with GraphicMagick or in Windows or Mac, search for the relevant guidances.

If you need to erase the lines inside the border, you can use the following command to do it. The `color` variable is the desired final output color of the border.

infile="overlay.png"
outfile="overlay_cleaned.png"
color="black"
w2=`convert $infile -format "%[fx:w-2]" info:`
h2=`convert $infile -format "%[fx:h-2]" info:`
convert $infile -background white -flatten -fill black +opaque white -bordercolor black -border 1 \
-fill none -draw "matte 2,2 floodfill matte $w2,2 floodfill matte $w2,$h2 floodfill matte 2,$h2 floodfill" \
-fill white +opaque none -fill black -opaque none \
-alpha off -morphology edgein octagon:3 \
-channel rgba -fill none +opaque white -fill $color -opaque white -shave 1x1 $outfile

See this thread in ImageMagick forum for more detail.