min_atom_idx = min(target_atom)#顶部原子index
max_atom_idx = max(target_atom)#如果两端连接,这表示另一端index
link_num_min = mol_t.GetAtomWithIdx(min_atom_idx).GetDegree()#头部原子连接数
link_num_max = mol_t.GetAtomWithIdx(max_atom_idx).GetDegree()#最大index连接数,如果是两边都连,就是另一端
if link_num_min < 4 and link_num_max < 3:#前面限制链上取代,后面限制作为骨架的中间情况,弊端:无法识别环一边只连脂肪链的情况
if min(bondring) != 0:
bondring_1 = min(bondring) - 1#前面的甲基bond
bondring = list(bondring)
bondring.append(bondring_1)