5

我有两个文件，第一个看起来像这样

125 6.678
435 9.084
234 8.874

依此类推，由我的程序生成大约 2,048,000 个条目。第二个文件是由 gnuplot 生成的文件，看起来有点像：

65 321456 985
78 98374  834
54 8493   848

依此类推，大约有 2,048,000 个条目。

现在我需要做的是使用 3D 中的 gnuplot 将第一个文件的第二列和第二个文件的两列相互绘制。我认为第一个任务是将它们全部放入同一个文件中，我只想编写一个简单的 c 程序来快速读取这两个文件并将相关列放入一个文件中，但我不知道该怎么做。我知道如何复制文件的全部内容并使用 ac 程序将它们写入另一个文件 - 例如，我有以下代码：

#include <stdio.h>
#include <stdlib.h>

int main (int argc, const char * argv[])
{
FILE *avalanche_size_BM;
FILE *avalanche_size_BM_2000;
char ch;


avalanche_size_BM = fopen("/Users/maheensiddiqui/Documents/MSc_Proj/avalanche_size_BM.dat","r");

if (avalanche_size_BM == NULL) 
{
    printf("I couldn't open.\n");
    exit(0);
}
avalanche_size_BM_2000 = fopen("/Users/maheensiddiqui/Desktop/avalanche_size_BM_2000.dat", "w");
if (avalanche_size_BM_2000 == NULL) 
{
    printf("I couldn't open.\n");
    exit(0);
}

printf("\n success!!");

while((ch=getc(avalanche_size_BM))!=EOF)
      putc(ch,avalanche_size_BM_2000);

fclose(avalanche_size_BM);
fclose(avalanche_size_BM_2000);



return(0);
}

但是我如何告诉它只读取第一个文件中的第二列和第二个文件中的前两列并将它们（而不是所有列）复制到第三个文件中，我可以用它来绘制我的 3-D 图.

任何帮助都感激不尽！

谢谢

4

2 回答 2

4

非 C 答案

这假设您只想要结果并且并不真正担心如何在 C 中执行它。

要使用简单的命令行工具生成文件，您可以使用paste并awk获取结果：

paste读取多个文件并将所有行组合在一起，由制表符分隔

因此，如果您有两个文件，例如 foo.txt 和 bar.txt，它们具有以下内容：

foo.txt：

a1 b1 c1 
a2 b2 c2

bar.txt：

d1 e1 f1
d2 e2 f2

paste foo.txt bar.txt输出：

a1 b1 c1    d1 e1 f1
a2 b2 c2    d2 e2 f2

然后，您可以使用 awk 过滤掉所需的列。

因此，例如，如果您想要列 2、4 和 5 (b1, d1, e1)，您可以使用以下内容：

paste foo.txt bar.txt | awk '{ print $2 " " $4 " " $5 }'

我将假设两个文件都有相同数量的条目在 OSX 中从 Finder -> Utilities 启动命令 shell

在 shell 提示下，您可以这样做：

$ cd Desktop
$ paste avalanche_size_BM.dat avalanche_size_BM_2000.dat | \
    awk '{ print $2 " " $4 " " $5 }' > avalanche_size_3d.dat

这将在桌面子目录中从第一个文件的第 2 列和第二个文件的前两列创建一个名为 avalanche_size_3d.dat 的文件（因为它们在粘贴在一起时最终位于位置 4 和 5）。

C-答案

这本质上就像编写一个小型粘贴实用程序

#include <stdio.h>
#include <stdlib.h>

int main()
{
     FILE *fp1 = fopen ("file1", "r");
     FILE *fp2 = fopen("file2", "r");
     FILE *out = fopen("fileout", "w");

     while (!feof(fp1) && !feof(fp2)) {
         char buf[256], col1[256], col2[256], col3[256];
         if ( fgets(buf, sizeof(buf), fp1) == 0 )
            break;
         sscanf(buf, "%s %s", col2, col1); /* we only need col1, col2 is reused later */
         if ( fgets(buf, sizeof(buf), fp2) == 0 ) 
            break;
         sscanf(buf, "%s %s", col2, col3); 
         fprintf(out, "%s %s %s\n", col1, col2, col3);
     }
     fclose(fp1);
     fclose(fp2);
     fclose(out);
}
于 2013-08-01T01:47:14.470 回答
2

在 C 语言中，最简单的方法是从两个文件中的每一个中读取一行，然后将这两行作为单行打印到第三个文件中。一旦你这样做了，你就可以解析你真正想要打印的每一行的哪些部分，但是将两个文件合并为一个可能足以让你用来gnuplot绘制你的数据。

要读入一行输入，您可以使用fgets(). 您可以为每行使用两个单独的缓冲区。从您从第一个文件获得的第一行，您删除最后的换行符，然后将其打印到新文件中。然后，打印第二个文件中的行。

char buf1[MAX_LINE];
char buf2[MAX_LINE];
FILE *infile1 = fopen(..., "r");
FILE *infile2 = fopen(..., "r");
FILE *outfile = fopen(..., "w");

while (fgets(buf1, sizeof(buf1), infile1) != 0
       && fgets(buf2, sizeof(buf2), infile2) != 0) {
    strchr(buf1, '\n')[0] = '\0';
    fprintf(outfile, "%s %s", buf1, buf2);
}

fclose(outfile);
fclose(infile1);
fclose(infile2);

读取一行输入后，您可以使用解析输入sscanf()来检索您有兴趣打印出的行部分：

int file1_column1;
double file1_column2;
sscanf(buf1, "%d %lf", &file1_column1, &file1_column2);

int file2_column1;
int file2_column2;
int file2_column3;
sscanf(buf2, "%d %d %d", &file2_column1, &file2_column2, &file2_column3);

您可以使用以下方法打印出值fprintf()：

fprintf(outfile, "%f %d %d\n", file1_column2, file2_column1, file2_column3);
于 2013-08-01T01:41:42.243 回答