37 std::string
const& id_area_fname, std::string
const& csv_fname,
38 std::vector<std::pair<std::size_t, double>>
const& ids_and_areas,
39 std::vector<MeshLib::Node*>
const& mesh_nodes)
41 std::ofstream ids_and_area_out(id_area_fname);
42 if (!ids_and_area_out)
44 OGS_FATAL(
"Unable to open the file '{:s}' - aborting.", id_area_fname);
46 std::ofstream csv_out(csv_fname);
49 OGS_FATAL(
"Unable to open the file '{:s}' - aborting.", csv_fname);
52 ids_and_area_out << std::setprecision(20);
53 csv_out << std::setprecision(20);
55 ids_and_area_out << ids_and_areas[0].first <<
" "
56 << ids_and_areas[0].second;
57 csv_out <<
"ID x y z area node_id\n";
58 csv_out << 0 <<
" " << *mesh_nodes[ids_and_areas[0].first]
59 << ids_and_areas[0].second <<
" " << ids_and_areas[0].first;
60 for (std::size_t k(1); k < ids_and_areas.size(); k++)
62 ids_and_area_out <<
"\n"
63 << ids_and_areas[k].first <<
" "
64 << ids_and_areas[k].second;
66 << k <<
" " << *mesh_nodes[ids_and_areas[k].first]
67 << ids_and_areas[k].second <<
" " << ids_and_areas[k].first;
69 ids_and_area_out <<
"\n";
73int main(
int argc,
char* argv[])
76 "Computes ids of mesh nodes that are in polygonal regions and resides "
77 "on the top surface. The polygonal regions have to be given in a gml- "
78 "or gli-file. The found mesh nodes and the associated area are written "
79 "as txt and csv data. The documentation is available at "
80 "https://docs.opengeosys.org/docs/tools/model-preparation/"
81 "computesurfacenodeidsinpolygonalregion.\n\n"
82 "OpenGeoSys-6 software, version " +
85 "Copyright (c) 2012-2025, OpenGeoSys Community "
86 "(http://www.opengeosys.org)",
88 TCLAP::ValueArg<std::string> mesh_in(
89 "m",
"mesh-input-file",
90 "Input (.vtu). The name of the file containing the input mesh",
true,
93 TCLAP::ValueArg<std::string> geo_in(
95 "Input (.gml). The name of the input file containing the polygons",
96 true,
"",
"INPUT_FILE");
99 TCLAP::ValueArg<std::string> gmsh_path_arg(
100 "",
"gmsh-path",
"Input (.msh). The path to the input binary file",
101 false,
"",
"INPUT_FILE");
102 cmd.add(gmsh_path_arg);
104 cmd.add(log_level_arg);
105 cmd.parse(argc, argv);
110 std::unique_ptr<MeshLib::Mesh const> mesh(
112 INFO(
"Mesh read: {:d} nodes, {:d} elements.", mesh->getNumberOfNodes(),
113 mesh->getNumberOfElements());
117 gmsh_path_arg.getValue());
119 INFO(
"Geometry '{:s}' read: {:d} points, {:d} polylines.",
124 Eigen::Vector3d
const dir({0.0, 0.0, -1.0});
127 auto computeElementTopSurfaceAreas =
128 [](
MeshLib::Mesh const& mesh, Eigen::Vector3d
const& d,
double angle)
130 std::unique_ptr<MeshLib::Mesh> surface_mesh(
137 std::vector<double> areas(computeElementTopSurfaceAreas(*mesh, dir, angle));
138 std::vector<MeshLib::Node*> all_sfc_nodes(
142 std::for_each(all_sfc_nodes.begin(), all_sfc_nodes.end(),
145 std::vector<MeshLib::Node*>
const& mesh_nodes(mesh->getNodes());
149 for (std::size_t j(0); j < plys.size(); j++)
151 if (!plys[j]->isClosed())
155 std::string polygon_name;
157 if (polygon_name.empty())
159 polygon_name =
"Polygon-" + std::to_string(j);
164 std::vector<std::pair<std::size_t, double>> ids_and_areas;
165 for (std::size_t k(0); k < all_sfc_nodes.size(); k++)
168 if (polygon.isPntInPolygon(surface_node))
170 ids_and_areas.emplace_back(surface_node.
getID(), areas[k]);
173 if (ids_and_areas.empty())
175 ERR(
"Polygonal part of surface '{:s}' doesn't contains nodes. No "
176 "output will be generated.",
182 std::string id_and_area_fname(out_path + polygon_name);
183 std::string csv_fname(out_path + polygon_name);
184 id_and_area_fname += std::to_string(j) +
".txt";
185 csv_fname += std::to_string(j) +
".csv";
187 "Polygonal part of surface '{}' contains {} nodes. Writing to "
188 "files '{}' and '{}'.",
190 ids_and_areas.size(),
193 writeToFile(id_and_area_fname, csv_fname, ids_and_areas, mesh_nodes);
196 std::for_each(all_sfc_nodes.begin(), all_sfc_nodes.end(),
197 std::default_delete<MeshLib::Node>());
void writeToFile(std::string const &id_area_fname, std::string const &csv_fname, std::vector< std::pair< std::size_t, double > > const &ids_and_areas, std::vector< MeshLib::Node * > const &mesh_nodes)