quantized_mesh_tile.h

// Copyright (c) 2018 Coronis Computing S.L. (Spain)
// All rights reserved.
//
// This file is part of EMODnet Quantized Mesh Generator for Cesium.
//
// This program is free software: you can redistribute it and/or modify
// it under the terms of the GNU General Public License as published by
// the Free Software Foundation, either version 3 of the License, or
// (at your option) any later version.
//
// This program is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
// GNU General Public License for more details.
//
// You should have received a copy of the GNU General Public License
// along with this program. If not, see <https://www.gnu.org/licenses/>.
//
// Author: Ricard Campos (ricardcd@gmail.com)

#ifndef EMODNET_QMGC_QUANTIZED_MESH_TILE_H
#define EMODNET_QMGC_QUANTIZED_MESH_TILE_H

#include "quantized_mesh.h"
#include "ellipsoid.h"
#include "tin_creation/tin_creation_cgal_types.h"


class QuantizedMeshTile :
public QuantizedMesh, public ctb::Tile
{
    friend class TerrainTiler;
    friend class QuantizedMeshTiler;

public:
    // --- Typedefs ---
    typedef TinCreation::Point_3   Point_3 ;
    typedef TinCreation::Vector_3  Vector_3 ;

    QuantizedMeshTile(const ctb::TileCoordinate &coord,
                      const Ellipsoid& e = WGS84Ellipsoid())
            : QuantizedMesh()
            , Tile(coord)
            , m_ellipsoid(e) {}

    QuantizedMeshTile(const char *fileName,
                      const ctb::TileCoordinate &coord,
                      const Ellipsoid& e = WGS84Ellipsoid())
            : QuantizedMesh(fileName)
            , Tile(coord)
            , m_ellipsoid(e) {}

    QuantizedMeshTile(const QuantizedMesh &qm,
                      const ctb::TileCoordinate &coord,
                      const Ellipsoid& e = WGS84Ellipsoid())
            : QuantizedMesh(qm)
            , Tile(coord)
            , m_ellipsoid(e) {}

    void convertUVHToLonLatHeight( const unsigned short &u, const unsigned short &v, const unsigned short &h,
                                   double &lon, double &lat, double &height ) const ;

    bool exportToOFF( const std::string &outFilePath, const bool& useRealWorldValues = false ) ;

    Point_3 horizonOcclusionPoint( const std::vector<Point_3> &pts, const Point_3 &center ) ;

private:
    // Function as described in https://cesium.com/blog/2013/05/09/computing-the-horizon-occlusion-point/
    // We ommit the ellipsoid here for simplicity (hard-coded)
    double computeHorizonOcclusionPointMagnitude( const Point_3 &position, const Vector_3 &scaledSpaceDirectionToPoint ) ;

    // --- Attributes ---
    Ellipsoid m_ellipsoid ; 
};

#endif //EMODNET_QMGC_QUANTIZED_MESH_TILE_H