geomodel_builder_2d_from_3d.h

Go to the documentation of this file.

/*
 * Copyright (c) 2012-2017, Association Scientifique pour la Geologie et ses
 * Applications (ASGA). All rights reserved.
 *
 * Redistribution and use in source and binary forms, with or without
 * modification, are permitted provided that the following conditions are met:
 *     * Redistributions of source code must retain the above copyright
 *       notice, this list of conditions and the following disclaimer.
 *     * Redistributions in binary form must reproduce the above copyright
 *       notice, this list of conditions and the following disclaimer in the
 *       documentation and/or other materials provided with the distribution.
 *     * Neither the name of ASGA nor the
 *       names of its contributors may be used to endorse or promote products
 *       derived from this software without specific prior written permission.
 *
 * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
 * AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO,
 * THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
 * PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL ASGA BE LIABLE FOR ANY DIRECT,
 * INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES
 * (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
 * LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND
 * ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
 * (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF
 * THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
 *
 *     http://www.ring-team.org
 *
 *     RING Project
 *     Ecole Nationale Superieure de Geologie - GeoRessources
 *     2 Rue du Doyen Marcel Roubault - TSA 70605
 *     54518 VANDOEUVRE-LES-NANCY
 *     FRANCE
 */

#pragma once

#include <ringmesh/basic/common.h>

#include <ringmesh/basic/geometry.h>

#include <ringmesh/geomodel/geomodel_builder.h>

namespace RINGMesh
{
    class RINGMESH_API GeoModelBuilder2DFrom3D : public GeoModelBuilder< 2 >
    {
    public:
        GeoModelBuilder2DFrom3D( GeoModel2D& geomodel2d,
            const GeoModel3D& geomodel3d_from,
            const Geometry::Plane& plane )
            : GeoModelBuilder( geomodel2d ),
              geomodel3d_from_( geomodel3d_from ),
              plane_( plane )
        {
            // Definition of v axis (upward) of the 2D space (along the plane)
            auto upward_point = plane_.origin + vec3{ 0., 0., 1. };
            vec3 v_axis_point_direction;
            std::tie( std::ignore, v_axis_point_direction ) =
                Distance::point_to_plane( upward_point, plane_ );
            if( ( plane_.origin - v_axis_point_direction ).length()
                < geomodel3d_from.epsilon() )
            {
                // Case where plane is sub-horizontal
                // (v axis is set towards 3D x direction)
                auto towards_x_point = plane_.origin + vec3{ 1., 0., 0. };
                std::tie( std::ignore, v_axis_point_direction ) =
                    Distance::point_to_plane( towards_x_point, plane_ );
                ringmesh_assert(
                    ( plane_.origin - v_axis_point_direction ).length()
                    < geomodel3d_from.epsilon() );
            }
            v_axis = normalize( v_axis_point_direction - plane_.origin );
            u_axis = cross( v_axis, plane_.normal );
        }

    protected:
        vec2 get_2d_coord( const vec3& coord3d )
        {
            return { dot( coord3d, u_axis ), dot( coord3d, v_axis ) };
        }

    protected:
        const GeoModel3D& geomodel3d_from_;
        const Geometry::Plane& plane_;
        vec3 u_axis{};
        vec3 v_axis{};
    };

    class RINGMESH_API GeoModelBuilder2DProjection
        : public GeoModelBuilder2DFrom3D
    {
    public:
        GeoModelBuilder2DProjection( GeoModel2D& geomodel2d,
            const GeoModel3D& geomodel3d_from,
            const Geometry::Plane& plane )
            : GeoModelBuilder2DFrom3D( geomodel2d, geomodel3d_from, plane )
        {
            info.set_geomodel_name( geomodel3d_from_.name() + "_projected" );
        }

        void build_geomodel();

    private:
        void copy_geomodel_3d_topology();

        void copy_geomodel_3d_geological_informations();

        void project_geomodel_3d_mesh_entities();

        std::vector< vec2 > compute_projected_vertices(
            const GeoModelMeshEntity3D& entity );
    };
}