GeoLib::IO::TINInterface Class Reference

Detailed Description

Interface for reading and writing Triangulated Irregular Network (TIN) file formats.

Definition at line 28 of file TINInterface.h.

#include <TINInterface.h>

Static Public Member Functions
static GeoLib::Surface *	readTIN (std::string const &fname, GeoLib::PointVec &pnt_vec, std::vector< std::string > *errors=nullptr)
static void	writeSurfaceAsTIN (GeoLib::Surface const &surface, std::string const &file_name)

Member Function Documentation

readTIN()

GeoLib::Surface * GeoLib::IO::TINInterface::readTIN ( std::string const & fname, GeoLib::PointVec & pnt_vec, std::vector< std::string > * errors = nullptr )

static

Reads TIN file

Parameters

fname	TIN file name
pnt_vec	a point vector to which triangle vertices are added
errors	a vector of error messages

Returns

a pointer to a GeoLib::Surface object created from TIN data. nullptr is returned if it fails to read the file.

Definition at line 27 of file TINInterface.cpp.

{
    // open file
    std::ifstream in(fname.c_str());
    if (!in)
    {
        WARN("readTIN(): could not open stream from {:s}.", fname);
        if (errors)
        {
            errors->push_back("readTINFile error opening stream from " + fname);
        }
        return nullptr;
    }
 
    auto* sfc = new GeoLib::Surface(pnt_vec.getVector());
    std::size_t id;
    MathLib::Point3d p0;
    MathLib::Point3d p1;
    MathLib::Point3d p2;
    std::string line;
    while (std::getline(in, line).good())
    {
        // allow empty lines
        if (line.empty())
        {
            continue;
        }
 
        // parse line
        std::stringstream input(line);
        // read id
        if (!(input >> id))
        {
            in.close();
            delete sfc;
            return nullptr;
        }
        // read first point
        if (!(input >> p0[0] >> p0[1] >> p0[2]))
        {
            ERR("Could not read coords of 1st point of triangle {:d}.", id);
            if (errors)
            {
                errors->push_back(
                    std::string("readTIN error: ") +
                    std::string(
                        "Could not read coords of 1st point in triangle ") +
                    std::to_string(id));
            }
            in.close();
            delete sfc;
            return nullptr;
        }
        // read second point
        if (!(input >> p1[0] >> p1[1] >> p1[2]))
        {
            ERR("Could not read coords of 2nd point of triangle {:d}.", id);
            if (errors)
            {
                errors->push_back(
                    std::string("readTIN error: ") +
                    std::string(
                        "Could not read coords of 2nd point in triangle ") +
                    std::to_string(id));
            }
            in.close();
            delete sfc;
            return nullptr;
        }
        // read third point
        if (!(input >> p2[0] >> p2[1] >> p2[2]))
        {
            ERR("Could not read coords of 3rd point of triangle {:d}.", id);
            if (errors)
            {
                errors->push_back(
                    std::string("readTIN error: ") +
                    std::string(
                        "Could not read coords of 3rd point in triangle ") +
                    std::to_string(id));
            }
            in.close();
            delete sfc;
            return nullptr;
        }
 
        // check area of triangle
        double const d_eps(std::numeric_limits<double>::epsilon());
        if (MathLib::calcTriangleArea(p0, p1, p2) < d_eps)
        {
            ERR("readTIN: Triangle {:d} has zero area.", id);
            if (errors)
            {
                errors->push_back(std::string("readTIN: Triangle ") +
                                  std::to_string(id) +
                                  std::string(" has zero area."));
            }
            delete sfc;
            return nullptr;
        }
 
        // determine size pnt_vec to insert the correct ids
        std::size_t const s(pnt_vec.getVector().size());
 
        std::size_t const pnt_pos_0(
            pnt_vec.push_back(new GeoLib::Point(p0, s)));
        std::size_t const pnt_pos_1(
            pnt_vec.push_back(new GeoLib::Point(p1, s + 1)));
        std::size_t const pnt_pos_2(
            pnt_vec.push_back(new GeoLib::Point(p2, s + 2)));
        // create new Triangle
        if (pnt_pos_0 != std::numeric_limits<std::size_t>::max() &&
            pnt_pos_1 != std::numeric_limits<std::size_t>::max() &&
            pnt_pos_2 != std::numeric_limits<std::size_t>::max())
        {
            sfc->addTriangle(pnt_pos_0, pnt_pos_1, pnt_pos_2);
        }
    }
 
    if (sfc->getNumberOfTriangles() == 0)
    {
        WARN("readTIN(): No triangle found.", fname);
        if (errors)
        {
            errors->push_back("readTIN error because of no triangle found");
        }
        delete sfc;
        return nullptr;
    }
 
    return sfc;
}

References MathLib::calcTriangleArea(), ERR(), GeoLib::TemplateVec< T >::getVector(), GeoLib::PointVec::push_back(), and WARN().

Referenced by main(), and FileIO::Legacy::readSurface().

writeSurfaceAsTIN()

void GeoLib::IO::TINInterface::writeSurfaceAsTIN ( GeoLib::Surface const & surface, std::string const & file_name )

static

Writes surface data into TIN file

Parameters

surface	surface object
file_name	TIN file name

Definition at line 162 of file TINInterface.cpp.

{
    std::ofstream os(file_name.c_str());
    if (!os)
    {
        WARN("writeSurfaceAsTIN(): could not open stream to {:s}.", file_name);
        return;
    }
    os.precision(std::numeric_limits<double>::max_digits10);
    const std::size_t n_tris(surface.getNumberOfTriangles());
    for (std::size_t l(0); l < n_tris; l++)
    {
        GeoLib::Triangle const& tri(*(surface[l]));
        os << l << " " << *(tri.getPoint(0)) << " " << *(tri.getPoint(1)) << " "
           << *(tri.getPoint(2)) << "\n";
    }
    os.close();
}

References GeoLib::Surface::getNumberOfTriangles(), GeoLib::Triangle::getPoint(), and WARN().

Referenced by main(), and FileIO::Legacy::writeTINSurfaces().

---

The documentation for this class was generated from the following files:

• GeoLib/IO/TINInterface.h
• GeoLib/IO/TINInterface.cpp