edem生成颗粒后，后处理导出颗粒中颗粒团形貌图片的c++代码

以下是一个简单的C++程序，可以将edem生成的颗粒团的形貌导出为图片。这个程序使用了OpenCV库来处理图像。

#include <iostream>
#include <iomanip>
#include <fstream>
#include <string>
#include <opencv2/core/core.hpp>
#include <opencv2/highgui/highgui.hpp>

using namespace std;
using namespace cv;

int main(int argc, char** argv)
{
    if (argc < 2) {
        cout << "Usage: " << argv[0] << " <input_file> [output_file]" << endl;
        return 1;
    }

    string input_file = argv[1];
    string output_file = "output.png";
    if (argc >= 3) {
        output_file = argv[2];
    }

    ifstream file(input_file);
    if (!file.is_open()) {
        cout << "Failed to open input file: " << input_file << endl;
        return 1;
    }

    int num_particles, num_dimensions;
    file >> num_particles >> num_dimensions;

    Mat image(512, 512, CV_8UC3, Scalar(255, 255, 255));

    for (int i = 0; i < num_particles; i++) {
        double x, y, z, radius;
        file >> x >> y >> z >> radius;

        Point center(x * 10 + 256, y * 10 + 256);
        int r = radius * 10;
        circle(image, center, r, Scalar(0, 0, 0), -1);
    }

    imwrite(output_file, image);
    cout << "Saved image to: " << output_file << endl;

    return 0;
}

此程序的输入文件应该是.edat文件，例如：

10 3
-0.016854 0.036321 0.001108 0.001000
-0.015609 0.038271 0.000822 0.001000
-0.013704 0.038652 0.000833 0.001000
-0.017669 0.034667 0.001151 0.001000
-0.014751 0.031810 0.000728 0.001000
-0.013601 0.036400 0.000760 0.001000
-0.011202 0.038575 0.000878 0.001000
-0.017399 0.032178 0.000546 0.001000
-0.015687 0.029238 0.001089 0.001000
-0.012709 0.034580 0.001112 0.001000

其中第一行是颗粒数和维度，接下来的每一行是一个颗粒的位置和半径。程序将这些颗粒绘制为一个黑色的圆圈，并将结果保存为PNG图像。你可以将其编译并运行，如下所示：

g++ -o export_particles export_particles.cpp `pkg-config --cflags --libs opencv`
./export_particles example.edat example.png

这将生成一个名为example.png的文件，其中包含颗粒团的形貌。