24小时热门版块排行榜    

查看: 3648  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 276求调剂 +4 路lyh123 2026-02-28 5/250 2026-03-02 11:20 by yuchj
[考研] 295求调剂 +8 19171856320 2026-02-28 8/400 2026-03-02 11:19 by yuchj
[考研] 281求调剂 +5 2026计算机_诚心 2026-03-01 8/400 2026-03-02 11:05 by 汪!?!
[考研] 284求调剂 +10 天下熯 2026-02-28 11/550 2026-03-02 11:03 by 无际的草原
[考研] 0856材料与化工,270求调剂 +8 YXCT 2026-03-01 9/450 2026-03-02 11:01 by 无际的草原
[考研] 欢迎采矿、地质、岩土、计算机、人工智能等专业的同学报考 +5 pin8023 2026-02-28 7/350 2026-03-02 10:33 by ZY,先生
[考研] 化工270求调剂 +6 什么名字qwq 2026-03-02 6/300 2026-03-02 10:21 by 黑!在干嘛
[考研] 调剂 +3 13853210211 2026-03-02 4/200 2026-03-02 10:16 by 13853210211
[考研] 材料学调剂 +10 提神豆沙包 2026-02-28 12/600 2026-03-02 09:26 by 李老师!
[考研] 282求调剂 +3 2103240126 2026-03-02 4/200 2026-03-02 09:25 by 汪!?!
[考研] 材料复试调剂 +4 学材料的点 2026-03-01 5/250 2026-03-02 08:26 by houyaoxu
[考研] 085600 英一数二272求调剂 5+3 vida_a 2026-03-01 7/350 2026-03-02 07:51 by ms629
[考研] 材料化工调剂 +12 今夏不夏 2026-03-01 13/650 2026-03-01 23:32 by L135790
[硕博家园] 博士自荐 +7 科研狗111 2026-02-26 11/550 2026-03-01 22:24 by 哲平L
[考研] 0856化工专硕求调剂 +12 董boxing 2026-03-01 12/600 2026-03-01 19:45 by 材子momo
[考研] 291分工科求调剂 +9 science饿饿 2026-03-01 10/500 2026-03-01 18:55 by 18137688336
[考研] 321求调剂一志愿东北林业大学材料与化工英二数二 +4 虫虫虫虫虫7 2026-03-01 7/350 2026-03-01 16:52 by caszguilin
[考研] 285求调剂 +8 满头大汗的学生 2026-02-28 8/400 2026-03-01 16:47 by caszguilin
[考研] 寻找调剂 +4 LYidhsjabdj 2026-02-28 4/200 2026-03-01 10:56 by sunny81
[考研] 085600材料工程一志愿中科大总分312求调剂 +8 吃宵夜1 2026-02-28 10/500 2026-02-28 20:27 by L135790
信息提示
请填处理意见