Skip to Content.
Sympa Menu

cgal-discuss - Re: [cgal-discuss] Atoms contributing to the Connolly surface

Subject: CGAL users discussion list

List archive

Re: [cgal-discuss] Atoms contributing to the Connolly surface


Chronological Thread 
  • From: Nico Kruithof <>
  • To:
  • Subject: Re: [cgal-discuss] Atoms contributing to the Connolly surface
  • Date: Wed, 29 Dec 2010 03:44:58 -0800
  • Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=mime-version:sender:in-reply-to:references:date :x-google-sender-auth:message-id:subject:from:to:content-type; b=YlypFyUPKb8FDQInSzNjQt98NK1J22Ctd4yjT/uHToifr5feDuaVmoBfiZzQ7pKJ7F y0Bsxa65UHj27bEp5yx0cdr+c1TgT/R4b9mhqmV0dj7L0UP4i9M0LkVYUMGRszS4GhTk w/rmrCJIsqPH/o2P/yQV7yNU5JeXnh2TokEtc=

Then, you should compute the Regular triangulation of the balls and extract the balls that contribute to a dual vertex lying outside the union of balls. To do that, you could iterate over all finite cells of the regular triangulation, check whether the dual point lies outside one (and therefore all) ball defining the cell. In that case all four vertices intersect the union of balls.

You could also use the skin surface package to mesh the union of balls.

Bests,
Nico

On Sun, Dec 26, 2010 at 6:18 PM, Francois Berenger <> wrote:
Nico Kruithof wrote:
Not exactly, the Connolly surface is the surface obtained by rolling a ball of radius r over the union of balls. If r==0, you are right, but not if r>0.

The rolling probe sphere should represent a solvent molecule in my case (water).


The easiest way is to compute the regular triangulation of the balls (weighted points) with the radius increased by r. Then, if a Voronoi cell of a weighted point (ball) is not entirely contained within that ball, the ball contributes to the Connoly surface. This shouldn't be too hard to implement.

I see. In fact, the dot surface should be OK for me, I don't really need the Connoly surface in fact.
I think I will do this outside of CGAL if it does not include such
high level constructs as molecule.get_dot_suface() or molecule.get_Connoly_suface().

Regards,
F.

Bests,
Nico


On Tue, Dec 21, 2010 at 10:27 AM, Laurent Rineau (GeometryFactory) < <mailto:>> wrote:

   On Tuesday 21 December 2010 09:29:41 Francois Berenger wrote:
    > Hi,
    >
    > Using CGAL, is it easy to extract the list of atoms
    > contributing to the Connolly surface for a globular
    > protein?

   Unless I am wrong, the Connolly surface is the skin surface with scale
   factor 1.

   See this chapter if you can find what you need:
    http://www.cgal.org/Manual/latest/doc_html/cgal_manual/Skin_surface_3/Chapter_main.html

   --
   Laurent Rineau, PhD
   R&D Engineer at GeometryFactory              http://www.geometryfactory.com/
   Release Manager of the CGAL Project       http://www.cgal.org/

   --
   You are currently subscribed to cgal-discuss.
   To unsubscribe or access the archives, go to
   https://lists-sop.inria.fr/wws/info/cgal-discuss




--
You are currently subscribed to cgal-discuss.
To unsubscribe or access the archives, go to
https://lists-sop.inria.fr/wws/info/cgal-discuss





Archive powered by MHonArc 2.6.16.

Top of Page