python:浮点数格式

如何解决python:浮点数格式

我正在编写 python 脚本,该脚本将一个点的坐标(表示 3D 结构的几何中心)保存在多模型文件 (result.pdb) 中,该文件包含 240 个相同原子的模型,其中 240 个不同的坐标 (XYZ) 分开由 ENDMDL 术语。这是我的python脚本的一个片段:

# this function calculates XYZ of the given atomic coordinates
def lig_center(nested_array_list):
    a = numpy.array(nested_array_list)
    mean = numpy.mean(a,axis=0)
    return mean[0],mean[1],mean[2]  # x,y,z


# make the list of the centres for each of the conformation
conf_xyz=[]
i=1
# loop each of the 240 conformations
for conf in d.ch.conformations:
    number=i
    # get XYZ of all atoms of each conformations
    xyz=conf.getCoords()
    # get XYZ of the geometrical center of these coordinates
    col=lig_center(xyz)
    # store geometrical center as in a list
    col1 = list(col)
    # append it to the existing list called conf_xyz
    conf_xyz.append(col1)
    # store coordinates of the centers in the output pdb file:
    with open(outputfilename+".pdb","a") as f:
    # for each of the model first print its index
        print >>f,'MODEL'+str(i)
    # then print the name of the atom (always C and its XYZ taken from the list col:
        print >>f,'ATOM      1  C   COM     1      '+'%.2f  %.2f  %.2f' %(col[0],col[1],col[2])+' 1.00 0.00'
        print >>f,'ENDMDL'
    i=i+1

这个方法在输出格式上的问题。这是我在 result.pdb 输出文件中的内容:

MODEL79
ATOM      1  C   COM     1      22.06  -14.69  6.02 1.00 0.00
ENDMDL
MODEL80
ATOM      1  C   COM     1      20.01  -17.91  4.69 1.00 0.00
ENDMDL
MODEL81
ATOM      1  C   COM     1      04.38  -9.05  -7.02 1.00 0.00
ENDMDL
MODEL82
ATOM      1  C   COM     1      3.12  -3.63  -9.41 1.00 0.00
ENDMDL
MODEL83
ATOM      1  C   COM     1      4.13  -6.49  -7.51 1.00 0.00
ENDMDL
MODEL84
ATOM      1  C   COM     1      3.22  -3.67  -9.70 1.00 0.00
ENDMDL
MODEL85
ATOM      1  C   COM     1      4.22  -4.74  -8.54 1.00 0.00
ENDMDL
MODEL86
ATOM      1  C   COM     1      1.70  -2.19  -12.49 1.00 0.00
ENDMDL
MODEL87
ATOM      1  C   COM     1      3.32  -3.98  -10.74 1.00 0.00
ENDMDL

这是我应该拥有的:

MODEL79
ATOM      1  C   COM     1      22.06  -14.69  06.02 1.00 0.00
ENDMDL
MODEL80
ATOM      1  C   COM     1      20.01  -17.91  04.69 1.00 0.00
ENDMDL
MODEL81
ATOM      1  C   COM     1      04.38  -09.05  -07.02 1.00 0.00
ENDMDL
MODEL82
ATOM      1  C   COM     1      03.12  -03.63  -09.41 1.00 0.00
ENDMDL
MODEL83
ATOM      1  C   COM     1      04.13  -06.49  -07.51 1.00 0.00
ENDMDL
MODEL84
ATOM      1  C   COM     1      03.22  -03.67  -09.70 1.00 0.00
ENDMDL
MODEL85
ATOM      1  C   COM     1      04.22  -04.74  -08.54 1.00 0.00
ENDMDL
MODEL86
ATOM      1  C   COM     1      01.70  -02.19  -12.49 1.00 0.00
ENDMDL
MODEL87
ATOM      1  C   COM     1      03.32  -03.98  -10.74 1.00 0.00

注意在第二种情况下,由于每个值的 xx.xx 格式,XYZ (6-8) 列完美对齐(我在 x.xx 中添加了零,例如将 4.22 -4.74 -8.54 转换为 04.22 - 04.74 -08.54 等)。

假设在我的python脚本末尾定义了输出文件的格式,通过

print >>f,col[2])+' 1.00 0.00' 

解决方案之一是将其替换为

 print >>f,'ATOM      1  C   COM     1      '+'%05.2f  %05.2f  %05.2f' %(col_B[0],col_B[1],col_B[2])+' 1.00 0.00'

它对正数做正确的工作,但对负数没有做任何事情:

MODEL7 
ATOM 1 C COM 1 06.09 -19.59 -0.90 1.00 0.00 
ENDMDL 
MODEL8 
ATOM 1 C COM 1 06.51 -19.86 -1.27 1.00 0.00
ENDMDL 

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。

相关推荐


依赖报错 idea导入项目后依赖报错,解决方案:https://blog.csdn.net/weixin_42420249/article/details/81191861 依赖版本报错:更换其他版本 无法下载依赖可参考:https://blog.csdn.net/weixin_42628809/a
错误1:代码生成器依赖和mybatis依赖冲突 启动项目时报错如下 2021-12-03 13:33:33.927 ERROR 7228 [ main] o.s.b.d.LoggingFailureAnalysisReporter : *************************** APPL
错误1:gradle项目控制台输出为乱码 # 解决方案:https://blog.csdn.net/weixin_43501566/article/details/112482302 # 在gradle-wrapper.properties 添加以下内容 org.gradle.jvmargs=-Df
错误还原:在查询的过程中,传入的workType为0时,该条件不起作用 <select id="xxx"> SELECT di.id, di.name, di.work_type, di.updated... <where> <if test=&qu
报错如下,gcc版本太低 ^ server.c:5346:31: 错误:‘struct redisServer’没有名为‘server_cpulist’的成员 redisSetCpuAffinity(server.server_cpulist); ^ server.c: 在函数‘hasActiveC
解决方案1 1、改项目中.idea/workspace.xml配置文件,增加dynamic.classpath参数 2、搜索PropertiesComponent,添加如下 <property name="dynamic.classpath" value="tru
删除根组件app.vue中的默认代码后报错:Module Error (from ./node_modules/eslint-loader/index.js): 解决方案:关闭ESlint代码检测,在项目根目录创建vue.config.js,在文件中添加 module.exports = { lin
查看spark默认的python版本 [root@master day27]# pyspark /home/software/spark-2.3.4-bin-hadoop2.7/conf/spark-env.sh: line 2: /usr/local/hadoop/bin/hadoop: No s
使用本地python环境可以成功执行 import pandas as pd import matplotlib.pyplot as plt # 设置字体 plt.rcParams['font.sans-serif'] = ['SimHei'] # 能正确显示负号 p
错误1:Request method ‘DELETE‘ not supported 错误还原:controller层有一个接口,访问该接口时报错:Request method ‘DELETE‘ not supported 错误原因:没有接收到前端传入的参数,修改为如下 参考 错误2:cannot r
错误1:启动docker镜像时报错:Error response from daemon: driver failed programming external connectivity on endpoint quirky_allen 解决方法:重启docker -> systemctl r
错误1:private field ‘xxx‘ is never assigned 按Altʾnter快捷键,选择第2项 参考:https://blog.csdn.net/shi_hong_fei_hei/article/details/88814070 错误2:启动时报错,不能找到主启动类 #
报错如下,通过源不能下载,最后警告pip需升级版本 Requirement already satisfied: pip in c:\users\ychen\appdata\local\programs\python\python310\lib\site-packages (22.0.4) Coll
错误1:maven打包报错 错误还原:使用maven打包项目时报错如下 [ERROR] Failed to execute goal org.apache.maven.plugins:maven-resources-plugin:3.2.0:resources (default-resources)
错误1:服务调用时报错 服务消费者模块assess通过openFeign调用服务提供者模块hires 如下为服务提供者模块hires的控制层接口 @RestController @RequestMapping("/hires") public class FeignControl
错误1:运行项目后报如下错误 解决方案 报错2:Failed to execute goal org.apache.maven.plugins:maven-compiler-plugin:3.8.1:compile (default-compile) on project sb 解决方案:在pom.
参考 错误原因 过滤器或拦截器在生效时,redisTemplate还没有注入 解决方案:在注入容器时就生效 @Component //项目运行时就注入Spring容器 public class RedisBean { @Resource private RedisTemplate<String
使用vite构建项目报错 C:\Users\ychen\work>npm init @vitejs/app @vitejs/create-app is deprecated, use npm init vite instead C:\Users\ychen\AppData\Local\npm-