<rt id="bn8ez"></rt>
<label id="bn8ez"></label>

  • <span id="bn8ez"></span>

    <label id="bn8ez"><meter id="bn8ez"></meter></label>

    posts - 431,  comments - 344,  trackbacks - 0

    如果你的環境還沒準備好, 可以官方網站看如何配置環境:http://openbabel.org/wiki/Install_Python_bindings
    1. 通過使用OBMol, OBAtom和OBBond來創建原子和鍵
    import openbabel

    mol = openbabel.OBMol()
    print 'Should print 0 (atoms)'
    print mol.NumAtoms()

    a = mol.NewAtom()
    a.SetAtomicNum(6)   # carbon atom
    a.SetVector(0.0, 1.0, 2.0) # coordinates

    b = mol.NewAtom()
    mol.AddBond(1, 2, 1)   # atoms indexed from 1
    print 'Should print 2 (atoms)'
    print mol.NumAtoms()
    print 'Should print 1 (bond)'
    print mol.NumBonds()

    mol.Clear();

    2. 通過OBConversion來讀取分子, 并輸出不同格式文件或字符串值
    import openbabel

    obConversion = openbabel.OBConversion()
    obConversion.SetInAndOutFormats("smi", "mdl")    //讀取smiles值, 然后輸出mdl值

    mol = openbabel.OBMol()
    obConversion.ReadString(mol, "C1=CC=CS1")

    print 'Should print 5 (atoms)'
    print mol.NumAtoms()

    mol.AddHydrogens()
    print 'Should print 9 (atoms) after adding hydrogens'
    print mol.NumAtoms()      //輸出原子個數

    outMDL = obConversion.WriteString(mol)

    3. 計算fp值
    import pybel
    smiles = ['CCCC', 'CCCN']
    mols = [pybel.readstring("smi", x) for x in smiles]   # Create two molecules from the SMILES
    fps = [x.calcfp() for x in mols]   # Calculate their fingerprints
    print fps[0].bits, fps[1].bits
    print fps[0].fp[0]

    mol2 = pybel.readstring('smi', 'c2ccc1ccccc1c2')
    fp2 = mol2.calcfp("FP4")
    print fp2
    print fp2.bits


    mol3 = pybel.readstring('smi', 'C1CCCCC1')
    fp3 = mol3.calcfp()

    print fp3.__or__(fp2)  //計算相似度值

    4. 讀取sdf文件
    #encoding=utf-8
    import pybel
    for mymol in pybel.readfile("sdf", "structures_all.sdf"):
        fp = mymol.calcfp("FP2")
        print fp

    5. 輸出txt文件和sdf文件

    print mymol.write("smi")    //'CCCC'
    mymol.write("smi", "outputfile.txt")
    largeSDfile = Outputfile("sdf", "multipleSD.sdf")
    largeSDfile.write(mymol)
    largeSDfile.write(myothermol)
    largeSDfile.close()

    posted on 2009-10-25 12:37 周銳 閱讀(3747) 評論(0)  編輯  收藏 所屬分類: ChemistryJava 、Openbabel
    主站蜘蛛池模板: 久久精品国产亚洲精品| 国产成在线观看免费视频| 亚洲国产午夜福利在线播放| 亚洲欧美日韩中文高清www777| 114一级毛片免费| 国产午夜亚洲精品不卡| 在线观看免费大黄网站| 亚洲色偷偷色噜噜狠狠99网| 成年女人午夜毛片免费看| 亚洲一卡一卡二新区无人区| 成人a视频片在线观看免费| 亚洲成_人网站图片| 永久免费观看的毛片的网站| 亚洲成aⅴ人片久青草影院按摩| 免费黄色小视频网站| 久久久久亚洲精品无码网址色欲| 凹凸精品视频分类国产品免费| 美女免费精品高清毛片在线视| 国产亚洲精品无码拍拍拍色欲| 精品免费视在线观看| 亚洲日韩国产精品无码av| 成人AV免费网址在线观看| 蜜桃传媒一区二区亚洲AV| 亚洲午夜无码片在线观看影院猛| 91在线视频免费观看| 亚洲日本乱码一区二区在线二产线| 久久久久久久久免费看无码| 蜜桃传媒一区二区亚洲AV| 亚洲日本va在线视频观看| 国产高清免费视频| 在线观看亚洲免费| 亚洲av无码片在线播放| 日韩欧毛片免费视频| 特a级免费高清黄色片| 久久精品国产精品亚洲蜜月| 91免费人成网站在线观看18| 日本亚洲中午字幕乱码| 国产亚洲成AV人片在线观黄桃| 美女视频黄免费亚洲| 国产精品免费在线播放| 亚洲国产精品成人精品小说|