[Chimera-users] how to select atom within a rectangular region

Eric Pettersen pett at cgl.ucsf.edu
Tue Nov 1 10:58:30 PDT 2016


	Well, it certainly could be done by a custom Python script, e.g.:

from chimera import openModels, Molecule, selection
sel_atoms = []
for m in openModels.list(modelTypes=[Molecule]):
	for a in m.atoms:
		x, y, z = a.coord().data()
		if x0 <= x <= x1 and y0 <= y <= y1 and z0 <= z <= z1:
			sel_atoms.append(a)
selection.setCurrent(sel_atoms)

This selection is working on the original coordinates, i.e. as provided by the input source.  If you want to include rotations/translations that occurred in Chimera, you would use the .xformCoord() method instead of the .coord() method.
	But perhaps if you described why you need to select atoms within a certain Cartesian box, we could provide better help.

—Eric

	Eric Pettersen
	UCSF Computer Graphics Lab


> On Oct 31, 2016, at 11:37 AM, Cherry Dons <cherrydons at gmail.com> wrote:
> 
> Hi Chimera Developers,
> 
> This is Dengpan, one of chimera users. I have a question about how to select atoms within a region? For example, I I want to select atoms within a rectangular region  
> x0 to x1,
> y0 to y1,
> z0 to z1.
> 
> Could you please give me a hint how to do it? Thanks.
> 
> Dengpan Dong
> 
> -- 
> Department of Materials Science and Engineering,
> University of Utah
> _______________________________________________
> Chimera-users mailing list: Chimera-users at cgl.ucsf.edu
> Manage subscription: http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://plato.cgl.ucsf.edu/pipermail/chimera-users/attachments/20161101/ab88c1b2/attachment.html>


More information about the Chimera-users mailing list