ParallelOrdinaryKriging/XYZFile.cpp at master · mcleary/ParallelOrdinaryKriging · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
//
//  XYZFile.cpp
//  ParallelDTM
//
//  Created by Thales Sabino on 8/25/16.
//  Copyright © 2016 Thales Sabino. All rights reserved.
//

#include "XYZFile.h"

#include <iostream>
#include <fstream>
#include <sstream>
#include <iterator>

using namespace std;

PointVector ReadXYZFile(const std::string& FilePath)
{
    PointVector Data;

    ifstream InputFile(FilePath, ios::in);

    cout << "Reading file ... " << flush;
    string FileContents(istreambuf_iterator<char>(InputFile), (istreambuf_iterator<char>()));
    stringstream FileContentsStream(FileContents);
    cout << "done. Parsing ... " << flush;

    string Line;
    while (getline(FileContentsStream, Line))
    {
        istringstream iss(Line);
        vector<string> Tokens{ istream_iterator<string>{iss}, istream_iterator<string>{} };

        if (Tokens.size() == 3 || Tokens.size() == 6)
        {
            Data.emplace_back(stof(Tokens[0]), stof(Tokens[1]), stof(Tokens[2]));
        }
        else if (Tokens.size() == 8)
        {
            Data.emplace_back(stof(Tokens[2]), stof(Tokens[3]), stof(Tokens[4]));
        }
    }
    cout << "done" << endl;

    return Data;
}

void WriteXYZFile(const std::string& Filepath, PointVector Points)
{
    cout << "Writing to " << Filepath << " ... " << flush;
    ofstream OutputFile(Filepath);
    for(auto Point : Points)
    {
        OutputFile << Point.x << " " << Point.y << " " << Point.z << endl;
    }
    cout << "done" << endl;
}