vtk文件编写

　　 在paraview中加载vtk文件，可以很好的显示三维空间图像，如下cpp代码：

#include <iostream>
#include <fstream>
#include <windows.h>
#include <ctime>

using namespace std;
int main() {
    cout << "this is beginning!" << endl;
    int Nx = , Ny = ;
    fstream outfile;
    char filename[] = "output.vtk";
    outfile.open(filename, ios::out);
    if (!outfile) {
        cout << "open failed" << endl;
    }
    outfile << "# vtk DataFile Version 2.0" << endl;
    outfile << filename << endl;
    outfile << "ASCII " << endl;
    outfile << "DATASET STRUCTURED_GRID" << endl;
    outfile << "DIMENSIONS " << Nx << " " << Ny << " " <<  << endl;
    outfile << "POINTS " << Nx * Ny *  << " float" << endl;
    for (int i = ; i < ; i++) {
        for (int j = ; j < ; j++) {
            outfile << i << " " << j << " " <<  << endl;
        }
    }

    outfile << "POINT_DATA " << Nx * Ny *  << endl;
    outfile << "SCALARS CON float 1" << endl;
    outfile << "LOOKUP_TABLE default" << endl;

    // 获取随机数
    srand(time());
    for (int i = ; i < ; i++)
    {
        for (int j = ; j < ; j++)
        {
            outfile << rand() %  << "\t";
        }
    }
    outfile.close();
    system("pause");
    return ;
}

格式是固定的，将输出的文件在paraview中打开即可。如下：