24小时热门版块排行榜    

查看: 3647  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料类求调剂 +11 wana_kiko 2026-02-28 14/700 2026-03-02 08:46 by 聪明的大松鼠
[考研] 085600 英一数二272求调剂 5+3 vida_a 2026-03-01 7/350 2026-03-02 07:51 by ms629
[考研] 哈工大计算机刘劼团队招生 +3 hit_aiot 2026-03-01 5/250 2026-03-02 07:48 by 得鹿梦鱼111
[考研] 0854复试调剂 276 +4 wmm9 2026-03-01 5/250 2026-03-02 06:41 by 汪!?!
[考研] 材料学硕318求调剂 +11 February_Feb 2026-03-01 13/650 2026-03-01 23:53 by ccp273206157
[考研] 材料化工调剂 +12 今夏不夏 2026-03-01 13/650 2026-03-01 23:32 by L135790
[考研] 材料学硕318求调剂 +5 February_Feb 2026-03-01 5/250 2026-03-01 23:31 by L135790
[考研] 292求调剂 +6 yhk_819 2026-02-28 6/300 2026-03-01 23:23 by 向上的胖东
[基金申请] 成果系统访问量大,请一小时后再尝试。---NSFC啥时候好哦,已经两天这样了 +4 NSFC2026我来了 2026-02-28 4/200 2026-03-01 22:37 by 铁门栓
[考研] 0856调剂 +5 刘梦微 2026-02-28 5/250 2026-03-01 22:30 by wang_dand
[硕博家园] 博士自荐 +7 科研狗111 2026-02-26 11/550 2026-03-01 22:24 by 哲平L
[考研] 298求调剂 +6 axyz3 2026-02-28 6/300 2026-03-01 19:00 by 18137688336
[考研] 328求调剂 +3 aaadim 2026-03-01 5/250 2026-03-01 17:29 by njzyff
[考研] 290求调剂 +9 材料专硕调剂; 2026-02-28 11/550 2026-03-01 17:21 by sunny81
[考研] 281求调剂 +4 2026计算机_诚心 2026-03-01 7/350 2026-03-01 17:20 by 2026计算机_诚心
[考研] 321求调剂一志愿东北林业大学材料与化工英二数二 +4 虫虫虫虫虫7 2026-03-01 7/350 2026-03-01 16:52 by caszguilin
[考研] 0856材料求调剂 +4 麻辣鱿鱼 2026-02-28 4/200 2026-03-01 16:51 by caszguilin
[考研] 295复试调剂 +3 简木ChuFront 2026-03-01 3/150 2026-03-01 14:27 by zzxw520th
[考研] 材料调剂 +4 爱擦汗的可乐冰 2026-02-28 4/200 2026-03-01 00:38 by 猫猫球alter
[硕博家园] 【博士招生】太原理工大学2026化工博士 +4 N1ce_try 2026-02-24 8/400 2026-02-26 08:40 by N1ce_try
信息提示
请填处理意见