24小时热门版块排行榜    

查看: 3679  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿北京化工大学材料与化工(085600)296求调剂 +5 稻妻小编 2026-03-26 5/250 2026-03-26 13:32 by yujianx
[考研] 284求调剂 +4 junqihahaha 2026-03-26 4/200 2026-03-26 12:36 by 张凯十八号
[考研] 环境专硕324分求调剂推荐 +5 轩小宁—— 2026-03-26 5/250 2026-03-26 12:05 by i_cooler
[考研] 316求调剂 +9 梁茜雯 2026-03-19 9/450 2026-03-25 23:01 by peike
[考研] 303求调剂 +6 蓝山月 2026-03-25 6/300 2026-03-25 22:47 by 418490947
[考研] 材料专硕 335 分求调剂 +4 拒绝冷暴力 2026-03-25 4/200 2026-03-25 18:45 by haxia
[考研] 材料277求调剂 +4 min3 2026-03-24 4/200 2026-03-25 15:29 by fch1983
[考研] 347求调剂 +4 L when 2026-03-25 4/200 2026-03-25 13:37 by cocolv
[考研] 281求调剂 +4 Koxui 2026-03-24 5/250 2026-03-25 11:38 by userper
[考研] 0854电子信息求调剂 324 +4 Promise-jyl 2026-03-23 4/200 2026-03-25 11:36 by Sugarlight
[考研] B区考研调剂 +4 yqdszhdap- 2026-03-22 5/250 2026-03-25 08:51 by baoball
[考研] 307求调剂 +5 超级伊昂大王 2026-03-24 5/250 2026-03-24 15:46 by 星空星月
[考研] 321求调剂 +4 Ymlll 2026-03-24 4/200 2026-03-24 14:44 by sprinining
[考研] 335求调剂 +4 yuyu宇 2026-03-23 5/250 2026-03-23 23:49 by Txy@872106
[考研] 一志愿北京化工大学 070300 学硕 336分 求调剂 +7 vv迷 2026-03-22 7/350 2026-03-23 23:44 by Txy@872106
[考研] 一志愿陕师大生物学071000,298分,求调剂 +3 SYA! 2026-03-23 3/150 2026-03-23 19:09 by macy2011
[考研] 328求调剂 +4 LHHL66 2026-03-23 4/200 2026-03-23 14:55 by lbsjt
[考研] 一志愿北京化工大学070300 学硕336求调剂 +5 vv迷 2026-03-21 8/400 2026-03-22 14:20 by ColorlessPI
[考研] 0703化学调剂 +4 妮妮ninicgb 2026-03-21 4/200 2026-03-21 18:39 by 学员8dgXkO
[考研] 一志愿 南京航空航天大学大学 ,080500材料科学与工程学硕 +5 @taotao 2026-03-20 5/250 2026-03-20 20:16 by JourneyLucky
信息提示
请填处理意见