版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機(jī)號驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

meatball1982

銅蟲 (小有名氣)

應(yīng)助: 8 (幼兒園)
金幣: -3.4
紅花: 8
帖子: 297
在線: 108.3小時(shí)
蟲號: 1530793
注冊: 2011-12-09
性別: GG
專業(yè): 理論和計(jì)算化學(xué)

[求助] 用C讀取xyz文件的坐標(biāo)信息

有一個(gè)由pdb轉(zhuǎn)成的xyz文件。
如圖，我想用C將其中的坐標(biāo)信息提取出來，存在一個(gè)新的文件或是數(shù)組當(dāng)中。

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

只能用C。其實(shí)用awk或其它的很容易實(shí)現(xiàn)。但這部分最后可能會(huì)在一個(gè)大的C當(dāng)中，所以只能用C。

我的思路是每次讀入一行，進(jìn)行判斷。但如何判斷我還沒有想好。因?yàn)槊啃械男畔㈤L度不同。

現(xiàn)在有其它的工作，比較著急。非常感謝。

要求，程序要有程序說明和相應(yīng)的注釋。

捕獲.GIF

回復(fù)此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺(tái)，不對該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對版權(quán)有異議，請聯(lián)系郵箱：xiaomuchong@tal.com
附件 1 : diala.xyz

2013-06-26 13:59:40, 1.41 K

» 猜你喜歡

本科新能源科學(xué)與工程，一志愿華理能動(dòng)285求調(diào)劑已經(jīng)有5人回復(fù)
一志愿華北電力大學(xué)能動(dòng)專碩，293，求調(diào)劑已經(jīng)有5人回復(fù)
調(diào)劑310 已經(jīng)有7人回復(fù)
食品工程專碩一志愿中海洋309求調(diào)劑已經(jīng)有6人回復(fù)
275求調(diào)劑已經(jīng)有14人回復(fù)
283求調(diào)劑已經(jīng)有3人回復(fù)
085602 化工專碩 338分求調(diào)劑已經(jīng)有12人回復(fù)
312，生物學(xué)求調(diào)劑已經(jīng)有3人回復(fù)
化學(xué)0703 調(diào)劑 306分一志愿211 已經(jīng)有4人回復(fù)
311（085601）求調(diào)劑已經(jīng)有3人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:

需要將dem轉(zhuǎn)化成xyz坐標(biāo)的文本形式，自己經(jīng)驗(yàn)有限，求流程。已經(jīng)有5人回復(fù)
xyz文件打開方法已經(jīng)有4人回復(fù)
讀取Google截圖坐標(biāo)生成TXT坐標(biāo)文件的程序已經(jīng)有14人回復(fù)
發(fā)布一個(gè)自己寫的色坐標(biāo)計(jì)算軟件CIE1931xy.V1.6 已經(jīng)有15人回復(fù)
【求助】怎么導(dǎo)出xyz坐標(biāo) 已經(jīng)有4人回復(fù)
【求助】求將vasp的輸出文件轉(zhuǎn)化成xyz坐標(biāo)的腳本已經(jīng)有11人回復(fù)

唉。還是學(xué)吧。

1樓 2013-06-26 13:59:40

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

綠遍山原

銅蟲 (小有名氣)

應(yīng)助: 64 (初中生)
金幣: 707.7
紅花: 5
帖子: 159
在線: 201.2小時(shí)
蟲號: 1128371
注冊: 2010-10-21
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

【答案】應(yīng)助回帖

★
感謝參與，應(yīng)助指數(shù) +1
jjdg: 金幣+1, 感謝參與 2013-06-27 00:33:56

只要坐標(biāo)就只讀每行前四列:
比如
fgetl(line, LINE_MAX, fp);
int n = sscanf(line, "%s %lg %lg %lg", name, &x, &y, &z);
if (n < 4) error("Invalid input"

贊一下

回復(fù)此樓

要夢遊，不要催眠。

2樓2013-06-26 15:17:39

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

meatball1982

銅蟲 (小有名氣)

應(yīng)助: 8 (幼兒園)
金幣: -3.4
紅花: 8
帖子: 297
在線: 108.3小時(shí)
蟲號: 1530793
注冊: 2011-12-09
性別: GG
專業(yè): 理論和計(jì)算化學(xué)

引用回帖:

2樓: Originally posted by 綠遍山原 at 2013-06-26 15:17:39
只要坐標(biāo)就只讀每行前四列:
比如
fgetl(line, LINE_MAX, fp);
int n = sscanf(line, "%s %lg %lg %lg", name, &x, &y, &z);
if (n < 4) error("Invalid input"

綠遍山原:
您好。
1.能把整個(gè)程序?qū)懗鰜韱幔?br /> 2.我查了一下，查到的fgetl都是matlab中的函數(shù)，這是一個(gè)C函數(shù)嗎？
非常感謝。

贊一下

回復(fù)此樓

唉。還是學(xué)吧。

3樓2013-06-26 17:41:10

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

xiaochong8693

木蟲之王 (文壇精英)

應(yīng)助: 58 (初中生)
金幣: 74351.8
散金: 3788
紅花: 233
帖子: 47068
在線: 2659.7小時(shí)
蟲號: 752823
注冊: 2009-04-20
專業(yè): 動(dòng)物學(xué)

祝樓主諸事順利，萬事如意

贊一下

回復(fù)此樓

4樓2013-06-26 18:11:58

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

綠遍山原

銅蟲 (小有名氣)

應(yīng)助: 64 (初中生)
金幣: 707.7
紅花: 5
帖子: 159
在線: 201.2小時(shí)
蟲號: 1128371
注冊: 2010-10-21
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

引用回帖:

3樓: Originally posted by meatball1982 at 2013-06-26 17:41:10
綠遍山原:
您好。
1.能把整個(gè)程序?qū)懗鰜韱幔?br /> 2.我查了一下，查到的fgetl都是matlab中的函數(shù)，這是一個(gè)C函數(shù)嗎？
非常感謝。...

是c的函數(shù)

回復(fù)此樓

要夢遊，不要催眠。

5樓2013-06-26 22:30:34

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

jackie1179

木蟲 (正式寫手)

應(yīng)助: 11 (小學(xué)生)
金幣: 1929.2
紅花: 1
帖子: 460
在線: 47.2小時(shí)
蟲號: 1119284
注冊: 2010-10-11
專業(yè): 生物信息學(xué)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感謝參與，應(yīng)助指數(shù) +1
jjdg: 金幣+2, 感謝參與 2013-06-27 00:34:06
meatball1982: 金幣+50, ★★★很有幫助, 程序完整，有說明。嘗試過，好用。 2013-06-27 15:56:29

CODE:

這個(gè)或許可以滿足你的要求

#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;

}

贊一下

回復(fù)此樓

6樓2013-06-26 22:47:21

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

meatball1982

銅蟲 (小有名氣)

應(yīng)助: 8 (幼兒園)
金幣: -3.4
紅花: 8
帖子: 297
在線: 108.3小時(shí)
蟲號: 1530793
注冊: 2011-12-09
性別: GG
專業(yè): 理論和計(jì)算化學(xué)

引用回帖:

6樓: Originally posted by jackie1179 at 2013-06-26 22:47:21
這個(gè)或許可以滿足你的要求
#include <stdio.h>
#include <stdlib.h>
#include <malloc.h>
#include <string.h>

#define BUFSIZE 512

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

稍加修改，正是我想要的。非常感謝！

贊一下

回復(fù)此樓

唉。還是學(xué)吧。

7樓2013-06-27 15:55:15

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

迷失的心

銅蟲 (初入文壇)

應(yīng)助: 3 (幼兒園)
金幣: 146.9
帖子: 37
在線: 35.5小時(shí)
蟲號: 2526371
注冊: 2013-06-29
專業(yè): 計(jì)算數(shù)學(xué)與科學(xué)工程計(jì)算

引用回帖:

想知道buf的大小為什么是512呢？別的數(shù)不行嗎？

贊一下

回復(fù)此樓

努力才能做到最好

8樓2013-06-30 10:58:20

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

jackie1179

木蟲 (正式寫手)

應(yīng)助: 11 (小學(xué)生)
金幣: 1929.2
紅花: 1
帖子: 460
在線: 47.2小時(shí)
蟲號: 1119284
注冊: 2010-10-11
專業(yè): 生物信息學(xué)

引用回帖:

8樓: Originally posted by 迷失的心 at 2013-06-30 10:58:20
想知道buf的大小為什么是512呢？別的數(shù)不行嗎？...

buf的大小可以是其他數(shù)，此值大小不固定，但應(yīng)盡量設(shè)的大一點(diǎn)，否則有可能出錯(cuò)。

贊一下

回復(fù)此樓

9樓2013-07-01 14:45:14

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

迷失的心

銅蟲 (初入文壇)

應(yīng)助: 3 (幼兒園)
金幣: 146.9
帖子: 37
在線: 35.5小時(shí)
蟲號: 2526371
注冊: 2013-06-29
專業(yè): 計(jì)算數(shù)學(xué)與科學(xué)工程計(jì)算

引用回帖:

9樓: Originally posted by jackie1179 at 2013-07-01 14:45:14
buf的大小可以是其他數(shù)，此值大小不固定，但應(yīng)盡量設(shè)的大一點(diǎn)，否則有可能出錯(cuò)。...

多謝，只是512這個(gè)數(shù)太常見了....

贊一下

回復(fù)此樓

努力才能做到最好

10樓2013-07-03 22:36:29

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 meatball1982 的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 食品工程專碩一志愿中海洋309求調(diào)劑 +3	小張zxy張 2026-03-26	6/300	2026-03-28 15:55 by 小張zxy張
[考研] 311（085601）求調(diào)劑 +3	liziyeyeye 2026-03-28	3/150	2026-03-28 15:29 by 1018329917
[考研] 0856，材料與化工321分求調(diào)劑 +12	大饞小子 2026-03-27	13/650	2026-03-28 10:56 by self2008
[考研] 295求調(diào)劑 +5	1428151015 2026-03-27	6/300	2026-03-28 04:04 by fmesaito
[考研] 材料292調(diào)劑 +12	橘頌思美人 2026-03-23	12/600	2026-03-27 15:44 by caszguilin
[考博] 26申博 +3	加油沖��！ 2026-03-26	3/150	2026-03-27 15:38 by cls512
[考研] 307求調(diào)劑 +8	超級伊昂大王 2026-03-24	9/450	2026-03-27 15:34 by 超級伊昂大王
[考研] 316求調(diào)劑 +5	Pigcasso 2026-03-24	5/250	2026-03-27 12:10 by zhshch
[考研] 315調(diào)劑 +4	0860求調(diào)劑 2026-03-26	5/250	2026-03-27 11:23 by wangjy2002
[考研] 材料學(xué)碩，求調(diào)劑 6+5	糖葫蘆888ll 2026-03-22	10/500	2026-03-27 08:18 by hypershenger
[考研] 071000生物學(xué)求調(diào)劑，初試成績343 +6	小小甜面團(tuán) 2026-03-25	6/300	2026-03-26 23:01 by 不吃魚的貓
[考研] 總分322求生物學(xué)/生化與分子/生物信息學(xué)相關(guān)調(diào)劑 +5	星沉uu 2026-03-26	6/300	2026-03-26 19:02 by macy2011
[考研] 一志愿南京郵電大學(xué) 288分材料考研求調(diào)劑 +3	jl0720 2026-03-26	3/150	2026-03-26 13:39 by zzll406
[考研] 299求調(diào)劑 +4	15188958825 2026-03-25	4/200	2026-03-25 22:56 by 418490947
[考研] 303求調(diào)劑 +6	藍(lán)山月 2026-03-25	6/300	2026-03-25 22:47 by 418490947
[考研] 機(jī)械學(xué)碩總分317求調(diào)劑�。。�！ +4	Acaciad 2026-03-25	4/200	2026-03-25 19:59 by hanserlol
[考研] 285求調(diào)劑 +3	AZMK 2026-03-24	3/150	2026-03-25 12:23 by userper
[考研] 材料專碩找調(diào)劑 +5	哈哈哈吼吼吼哈 2026-03-23	5/250	2026-03-24 19:07 by 了了了了。。
[考研] 344求調(diào)劑 +3	desto 2026-03-24	3/150	2026-03-24 10:09 by 搏擊518
[考研] 285求調(diào)劑 +6	ytter 2026-03-22	6/300	2026-03-22 12:09 by 星空星月