24小时热门版块排行榜    

查看: 3664  |  回复: 9

meatball1982

铜虫 (小有名气)

[求助] 用C读取xyz文件的坐标信息

有一个由pdb转成的xyz文件。
如图,我想用C将其中的坐标信息提取出来,存在一个新的文件或是数组当中。

7.688102   -2.282595    2.994640
7.303011   -1.469669    3.647740
…………
7.158374    3.922640   -0.572473

只能用C。其实用awk或其它的很容易实现。但这部分最后可能会在一个大的C当中,所以只能用C。

我的思路是每次读入一行,进行判断。但如何判断我还没有想好。因为每行的信息长度不同。

现在有其它的工作,比较着急。非常感谢。

要求,程序要有程序说明和相应的注释。用C读取xyz文件的坐标信息
捕获.GIF
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : diala.xyz
  • 2013-06-26 13:59:40, 1.41 K

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

唉。还是学吧。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

绿遍山原

铜虫 (小有名气)

【答案】应助回帖


感谢参与,应助指数 +1
jjdg: 金币+1, 感谢参与 2013-06-27 00:33:56
只要坐标就只读每行前四列:
比如
fgetl(line, LINE_MAX, fp);
int n = sscanf(line, "%s %lg %lg %lg", name, &x, &y, &z);
if (n < 4) error("Invalid input"
要夢遊,不要催眠。
2楼2013-06-26 15:17:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meatball1982

铜虫 (小有名气)

引用回帖:
2楼: Originally posted by 绿遍山原 at 2013-06-26 15:17:39
只要坐标就只读每行前四列:
比如
fgetl(line, LINE_MAX, fp);
int n = sscanf(line, "%s %lg %lg %lg", name, &x, &y, &z);
if (n < 4) error("Invalid input"

绿遍山原:
您好。
1.能把整个程序写出来吗?
2.我查了一下,查到的fgetl都是matlab中的函数,这是一个C函数吗?
非常感谢。
唉。还是学吧。
3楼2013-06-26 17:41:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaochong8693

木虫之王 (文坛精英)

祝楼主诸事顺利,万事如意
4楼2013-06-26 18:11:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

绿遍山原

铜虫 (小有名气)

引用回帖:
3楼: Originally posted by meatball1982 at 2013-06-26 17:41:10
绿遍山原:
您好。
1.能把整个程序写出来吗?
2.我查了一下,查到的fgetl都是matlab中的函数,这是一个C函数吗?
非常感谢。...

是c的函数
要夢遊,不要催眠。
5楼2013-06-26 22:30:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jackie1179

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与,应助指数 +1
jjdg: 金币+2, 感谢参与 2013-06-27 00:34:06
meatball1982: 金币+50, ★★★很有帮助, 程序完整,有说明。尝试过,好用。 2013-06-27 15:56:29
CODE:
这个或许可以满足你的要求
#include <stdio.h>
#include <stdlib.h>
#include <malloc.h>
#include <string.h>

#define BUFSIZE 512

int main(int argc,char *argv[])
{
  FILE *ifp,*ofp;
  char buf[BUFSIZE];
  int number;       // record entry number
  int i;
  double *x,*y,*z;

  if((ifp = fopen("diala.xyz","r+")) == NULL)   // input file handle
  {
    printf("Cannot open this file!\n");
        exit(1);
  }
  if((ofp = fopen(argv[1],"w+")) == NULL)   // output filename obtained from first commandline parameter
  {
    printf("Cannot write to this file!\n");
        exit(1);
  }
  fgets(buf,BUFSIZE,ifp);   // read record entry number
  number = atoi(buf);

  // memory allocation for three arrays
  x = (double *)malloc(sizeof(double) * number);
  y = (double *)malloc(sizeof(double) * number);
  z = (double *)malloc(sizeof(double) * number);


  i = 0;
  while(fgets(buf,BUFSIZE,ifp) != NULL)  // read info line by line
  {
    sscanf(buf,"%*d%*s%lf%lf%lf",&x[i],&y[i],&z[i]);
        fprintf(ofp,"%lf  %lf  %lf\n",x[i],y[i],z[i]);  // write into output file
        i++;
  }

  fclose(ifp);
  fclose(ofp);

  // release memory
  free(x);
  free(y);
  free(z);
  

  return 0;

}

6楼2013-06-26 22:47:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meatball1982

铜虫 (小有名气)

引用回帖:
6楼: Originally posted by jackie1179 at 2013-06-26 22:47:21
这个或许可以满足你的要求
#include <stdio.h>
#include <stdlib.h>
#include <malloc.h>
#include <string.h>

#define BUFSIZE 512

int main(int argc,char *argv[])
{
  F ...

稍加修改,正是我想要的。非常感谢!
唉。还是学吧。
7楼2013-06-27 15:55:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

迷失的心

铜虫 (初入文坛)

引用回帖:
6楼: Originally posted by jackie1179 at 2013-06-26 22:47:21
这个或许可以满足你的要求
#include <stdio.h>
#include <stdlib.h>
#include <malloc.h>
#include <string.h>

#define BUFSIZE 512

int main(int argc,char *argv[])
{
  F ...

想知道buf的大小为什么是512呢?别的数不行吗?
努力才能做到最好
8楼2013-06-30 10:58:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jackie1179

木虫 (正式写手)

引用回帖:
8楼: Originally posted by 迷失的心 at 2013-06-30 10:58:20
想知道buf的大小为什么是512呢?别的数不行吗?...

buf的大小可以是其他数,此值大小不固定,但应尽量设的大一点,否则有可能出错。
9楼2013-07-01 14:45:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

迷失的心

铜虫 (初入文坛)

引用回帖:
9楼: Originally posted by jackie1179 at 2013-07-01 14:45:14
buf的大小可以是其他数,此值大小不固定,但应尽量设的大一点,否则有可能出错。...

多谢,只是512这个数太常见了....
努力才能做到最好
10楼2013-07-03 22:36:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 meatball1982 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料292调剂 +8 橘颂思美人 2026-03-23 8/400 2026-03-24 16:33 by laoshidan
[考博] 申博26年 +4 八6八68 2026-03-19 4/200 2026-03-24 15:49 by 小Ben呵呵
[考研] 085404电子信息284分求调剂 +4 13659058978 2026-03-24 4/200 2026-03-24 12:15 by syl20081243
[考研] 【双一流院校新能源、环境材料,材料加工与模拟招收大量调剂】 +4 Higraduate 2026-03-22 7/350 2026-03-24 11:23 by 种大树
[考研] 一志愿山东大学药学学硕求调剂 +3 开开心心没烦恼 2026-03-23 4/200 2026-03-24 00:06 by 开开心心没烦恼
[考研] 求材料,环境专业调剂 +3 18567500178 2026-03-18 3/150 2026-03-23 23:50 by 热情沙漠
[考研] 一志愿北京化工大学 070300 学硕 336分 求调剂 +7 vv迷 2026-03-22 7/350 2026-03-23 23:44 by Txy@872106
[考研] 化学308分求调剂 +3 你好明天你好 2026-03-23 3/150 2026-03-23 20:11 by macy2011
[考研] 0703化学求调剂 +4 奶油草莓. 2026-03-22 5/250 2026-03-23 19:37 by pswait
[考研] 336求调剂 +4 收到VS 2026-03-20 4/200 2026-03-23 19:02 by macy2011
[考研] 324求调剂 +6 lucky呀呀呀鸭 2026-03-20 6/300 2026-03-22 16:01 by ColorlessPI
[考研] 求调剂 +5 Zhangbod 2026-03-21 7/350 2026-03-22 13:13 by Zhangbod
[考研] 考研调剂 +4 来好运来来来 2026-03-21 4/200 2026-03-22 12:15 by 星空星月
[考研] 286求调剂 +10 Faune 2026-03-21 10/500 2026-03-21 23:34 by 314126402
[考研] 初试 317 +7 半拉月丙 2026-03-20 7/350 2026-03-21 22:26 by peike
[考研] 求调剂 +4 要好好无聊 2026-03-21 4/200 2026-03-21 18:57 by 学员8dgXkO
[考研] 311求调剂 +3 勇敢的小吴 2026-03-20 3/150 2026-03-21 17:40 by ColorlessPI
[考研] 一志愿华中科技大学,080502,354分求调剂 +5 守候夕阳CF 2026-03-18 5/250 2026-03-21 01:06 by JourneyLucky
[考研] 一志愿华南师大 070300(化学)304分求调剂 +3 0703武芊慧雪304 2026-03-18 3/150 2026-03-21 00:48 by JourneyLucky
[考研] 材料学硕318求调剂 +5 February_Feb 2026-03-19 5/250 2026-03-19 23:51 by 23Postgrad
信息提示
请填处理意见