Hi, Esben Jannik Bjerrum,
    I have studied how to do it, though it may not be the best way. In 
`rdMolDraw` there are `getAtomCoords` and `getDrawCoords` in C++ APIs    
http://www.rdkit.org/docs/cppapi/classRDKit_1_1MolDraw2D.html#abd327050dfa838543103d1f2c5f28f23
     But these APIs are not wrapped for python.  So you may have to add the 
wrapper into the source and compile it manually.
Here are some codes I made in older version:### in 
`rdkit/Code/GraphMol/MolDraw2D/Wrap` ###namespace RDKit {....Point2D 
getDrawCoordsViaId(MolDraw2D &self,  int at_num) {

  return  self.getDrawCoords(at_num);

}

Point2D getAtomCoordsViaId(MolDraw2D &self, int at_num) {

  return self.getAtomCoords(at_num);

} 

}...BOOST_PYTHON_MODULE(rdMolDraw2D) 
{...python::class_<RDKit::MolDraw2D,boost::noncopyable>("MolDraw2D",docString.c_str(),python::no_init)
    .def ...)    .def("getDrawCoords", RDKit::getDrawCoordsViaId, "a")     
.def("getAtomCoords", RDKit::getAtomCoordsViaId,"b")

    ;...

(I don't know why I added "a" and "b" and they may not necessary. I am not good 
at C++ and boost-python)        BTW, I wonder is it possible to open these APIs 
officially ?

Hongbin Yang 

 From: Esben Jannik Bjerrum via Rdkit-discussDate: 2017-06-21 22:38To: 
rdkit-discuss@lists.sourceforge.netSubject: [Rdkit-discuss] Reverse scale of 
svg coordinates to atom coordinates
Hi RDkitters,  
 I'm experimenting a bit with an application with some user interactivity. I 
get the 
SVG coordinates from the Mol SVG drawing from the user interaction and 
need to get back to the atom coordinates with the goal of identifying 
the atom nearest the selected coordinates (or is there a smarter way to 
achieve that goal?).
Is this possible from Python currently?
I see that there is a MolDraw2D::getAtomCoords function in the cpp code for 
MolDraw2D, but I can't see it from the Python side, and there don't seem to be 
a way to get the scaling from the MolDraw2DSVG object.

As
 I understand it, the coordinates from the molecule are offset and 
scaled (and flipped for y) to fit the drawing canvas of the specified 
size. To get back to the original atom coordinates I must somehow 
reverse the transformation. Here's some script snippets illustrating 
what I try to achieve.

#Get som SVG depiction of a mol
mol = Chem.MolFromSmiles('CCCC')
mc = Chem.Mol(mol.ToBinary())
rdDepictor.Compute2DCoords(mc)
drawer = rdMolDraw2D.MolDraw2DSVG(300,300)
drawer.DrawMolecule(mc)
drawer.FinishDrawing()
svg = drawer.GetDrawingText().replace('svg:','')

#....#Visualization and User interaction code here gives SVG coordinates
#....
svg_x = 271.0svg_y = 237.0

#Is there a function to scale back the coordinates? alternatively get the 
scaling and the offset from drawer and handle it manually
atomcoords = drawer.getAtomCoords((svg_x, svg_y))
#But this function doesn't exist:-(
#...#Continue working with the rdkit mol

I would welcome some hints or suggestions.

Esben Jannik Bjerrum
cand.pharm, Ph.D
/Sent from my Ubuntu Touch Phone

Phone +45 2823 8009
http://dk.linkedin.com/in/esbenbjerrum
http://www.wildcardconsulting.dk
------------------------------------------------------------------------------
Check out the vibrant tech community on one of the world's most
engaging tech sites, Slashdot.org! http://sdm.link/slashdot
_______________________________________________
Rdkit-discuss mailing list
Rdkit-discuss@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/rdkit-discuss

Reply via email to