TIN2VTK.cpp

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#include <tclap/CmdLine.h>
 
#include "BaseLib/TCLAPArguments.h"
 
// STL
#include <memory>
#include <string>
#include <vector>
 
#include "BaseLib/FileTools.h"
#include "BaseLib/MPI.h"
#include "InfoLib/GitInfo.h"
 
// GeoLib
#include "GeoLib/IO/TINInterface.h"
#include "GeoLib/Point.h"
#include "GeoLib/PointVec.h"
#include "GeoLib/Surface.h"
#include "MeshLib/IO/VtkIO/VtuInterface.h"
 
// MeshLib
#include "BaseLib/Logging.h"
#include "MeshLib/Mesh.h"
#include "MeshToolsLib/convertMeshToGeo.h"
 
int main(int argc, char* argv[])
{
    TCLAP::CmdLine cmd(
        "Converts TIN file into VTU file.\n\n"
        "OpenGeoSys-6 software, version " +
            GitInfoLib::GitInfo::ogs_version +
            ".\n"
            "Copyright (c) 2012-2026, OpenGeoSys Community "
            "(http://www.opengeosys.org)",
        ' ', GitInfoLib::GitInfo::ogs_version);
    TCLAP::ValueArg<std::string> inArg("i", "input-tin-file",
                                       "Input (.tin). The name "
                                       "of the input file "
                                       "containing the input TIN ",
                                       true, "", "INPUT_FILE");
    cmd.add(inArg);
    TCLAP::ValueArg<std::string> outArg(
        "o", "output-vtu-file",
        "Output (.vtu). The name of the output file "
        "the mesh will be written to",
        true, "", "OUTPUT_FILE");
    cmd.add(outArg);
    auto log_level_arg = BaseLib::makeLogLevelArg();
    cmd.add(log_level_arg);
    cmd.parse(argc, argv);
 
    BaseLib::MPI::Setup mpi_setup(argc, argv);
    BaseLib::initOGSLogger(log_level_arg.getValue());
 
    INFO("reading the TIN file...");
    const std::string tinFileName(inArg.getValue());
    std::string point_vec_name{"SurfacePoints"};
    GeoLib::PointVec point_vec{point_vec_name, std::vector<GeoLib::Point*>{},
                               GeoLib::PointVec::NameIdMap{}};
    std::unique_ptr<GeoLib::Surface> sfc(
        GeoLib::IO::TINInterface::readTIN(tinFileName, point_vec));
    if (!sfc)
    {
        return EXIT_FAILURE;
    }
    INFO("TIN read:  {:d} points, {:d} triangles", point_vec.size(),
         sfc->getNumberOfTriangles());
 
    INFO("converting to mesh data");
    std::unique_ptr<MeshLib::Mesh> mesh(MeshToolsLib::convertSurfaceToMesh(
        *sfc, BaseLib::extractBaseNameWithoutExtension(tinFileName),
        std::numeric_limits<double>::epsilon()));
    INFO("Mesh created: {:d} nodes, {:d} elements.", mesh->getNumberOfNodes(),
         mesh->getNumberOfElements());
 
    INFO("Write it into VTU");
    MeshLib::IO::VtuInterface writer(mesh.get());
    writer.writeToFile(outArg.getValue());
 
    return EXIT_SUCCESS;
}