#include <vtkSmartPointer.h>
#include <vtkSimplePointsReader.h>
#include <vtkPolyDataMapper.h>
#include <vtkActor.h>
#include <vtkProperty.h>
#include <vtkRenderer.h>
#include <vtkRenderWindow.h>
#include <vtkRenderWindowInteractor.h>
#include <iostream>
#include 
 
<fstream>
#include <sstream>
#include <string>
#include <vector>
#include <vtkVertexGlyphFilter.h>


//不新增以下程式碼，不能執行
#ifndef INITIAL_OPENGL
#define INITIAL_OPENGL

#include <vtkAutoInit.h>
VTK_MODULE_INIT(vtkRenderingOpenGL2)
VTK_MODULE_INIT(vtkInteractionStyle)
#endif
 


using namespace std;
//刪除字串中空格，製表符tab等無效字元
string Trim(string& str)
{
    //str.find_first_not_of(" \t\r\n"),在字串str中從索引0開始，返回首次不匹配"\t\r\n"的位置
    str.erase(0, str.find_first_not_of(" \t\r\n"));
    str.erase(str.find_last_not_of(" \t\r\n") + 1);
    return str;
}

int main()
{
    vtkSmartPointer<vtkPoints> points = vtkSmartPointer<vtkPoints>::New();
    
 
double x, y, z;
    ifstream fin("D:\\QTData\\Vtk-new\\ceshi1.csv"); //開啟檔案流操作-絕對路徑
    string line;
    while (getline(fin, line))   //整行讀取，換行符“\n”區分，遇到檔案尾標誌eof終止讀取
    {
        //cout <<"原始字串："<< line << endl; //整行輸出
        istringstream sin(line); //將整行字串line讀入到字串流istringstream中
        vector<string> fields; //宣告一個字串向量
        string field;
        while (getline(sin, field, ',')) //將字串流sin中的字元讀入到field字串中，以逗號為分隔符
        {
            fields.push_back(field); //將剛剛讀取的字串新增到向量fields中
        }

        string sx = Trim(fields[0]);
        string sy = Trim(fields[1]);
        string sz = Trim(fields[2]);

        //字串轉換成double型數字
        istringstream streamx, streamy, streamz;
        streamx.str(sx), streamy.str(sy), streamz.str(sz);
        streamx >> x, streamy >> y, streamz >> z;
        points->InsertNextPoint(x, y, z); //新讀取的資料賦予點的幾何結構  
    }
    fin.close();  //關閉檔案
    vtkSmartPointer<vtkPolyData> polyData = vtkSmartPointer<vtkPolyData>::New();
    polyData->SetPoints(points);

    //定義二維點剖分
    vtkDelaunay2D* delaunay = vtkDelaunay2D::New();
    delaunay->SetInputData(polyData);
    delaunay->Update();

    vtkSmartPointer<vtkVertexGlyphFilter> glyphFilter = vtkSmartPointer<vtkVertexGlyphFilter>::New();
#if VTK_MAJOR_VERSION <= 5  
    glyphFilter->SetInputConnection(polyData->GetProducerPort());
#else  
    glyphFilter->SetInputData(polyData);
#endif  
    glyphFilter->Update();

    vtkVertexGlyphFilter* vertexGlyphFilter = vtkVertexGlyphFilter::New();
    vertexGlyphFilter->SetInputData(polyData);
    vertexGlyphFilter->Update();

    // Visualize
    vtkSmartPointer<vtkPolyDataMapper> mapper = vtkSmartPointer<vtkPolyDataMapper>::New();
    mapper->SetInputConnection(vertexGlyphFilter->GetOutputPort());

    vtkSmartPointer<vtkActor> actor = vtkSmartPointer<vtkActor>::New();
    actor->SetMapper(mapper);
    actor->GetProperty()->SetPointSize(4);
    actor->GetProperty()->SetColor(0.0, 0.0, 0.0);

    //三角渲染
    vtkPolyDataMapper* triangleMapper = vtkPolyDataMapper::New();
    triangleMapper->SetInputConnection(delaunay->GetOutputPort());
    triangleMapper->Update();
    vtkActor* triangleActor = vtkActor::New();
    triangleActor->SetMapper(triangleMapper);

    vtkSmartPointer<vtkRenderer> renderer = vtkSmartPointer<vtkRenderer>::New();
    renderer->AddActor(actor);
    renderer->AddActor(triangleActor);

    renderer->SetBackground(1.3, 1.6, 1.3); // Background color green

    vtkSmartPointer<vtkRenderWindow> renderWindow = vtkSmartPointer<vtkRenderWindow>::New();
    renderWindow->AddRenderer(renderer);


    vtkSmartPointer<vtkRenderWindowInteractor> renderWindowInteractor = vtkSmartPointer<vtkRenderWindowInteractor>::New();
    renderWindowInteractor->SetRenderWindow(renderWindow);

    renderWindow->Render();
    renderWindowInteractor->Start();

    return EXIT_SUCCESS;
}