CscCartography.h

Go to the documentation of this file.

#ifndef CscCartography_H_
#define CscCartography_H_
 
#include "Util/conscience_log.h"
 
#include "Axiomes/Csc3dTypes.h"
 
#include <string>
 
using std::string, std::vector, std::exception;
 
using namespace conscience_core::axiomes;
using namespace conscience_utils;
 
namespace conscience_utils::Cartography {
 
enum CscProjectionLambertType {
    CC42 = 1,
    CC43 = 2,
    CC44 = 3,
    CC45 = 4,
    CC46 = 5,
    CC47 = 6,
    CC48 = 7,
    CC49 = 8,
    CC50 = 9,
};
 
CscPoint3d *projectWGS84toLambertCC(const CscPoint3d *position, CscProjectionLambertType projection);
 
double isometricLatitude(double phi, double e);
 
double getDistanceBetweenTwoWGS84Points(const CscPoint3d *firstPoint, const CscPoint3d *secondPoint);
 
CscPoint3d *projectWGS84PointWithDistanceAndAngle(const CscPoint3d *initialPoint, double distance, double angle);
 
CscPoint3d *calculatePositionFromGps(double lat, double lon,
                                     double latCornerBLScene, double longCornerBLScene,
                                     double latCornerBRScene, double longCornerBRScene,
                                     double latCornerTLScene, double longCornerTLScene,
                                     double latCornerTRScene, double longCornerTRScene);
 
 
CscPoint3d findPointAtDistanceFrom(
    double latitude, double longitude, double rotationFromNorthRad, double distanceMeters);
 
CscPoint3d pictureCoordToGpsPositionWithFOV(
    double latitude, double longitude, double altitudeInMeters,
    double widthFOV, double heightFOV, double heading, double cameraOrientation,
    double photoX, double photoY, double imageHeight, double imageWidth);
 
    
enum CscCartographicProjectionType {
    EPSG9842, // CC42
    EPSG9843, // CC43
    EPSG9844, // CC44
    EPSG9845, // CC45
    EPSG9846, // CC46
    EPSG9847, // CC47
    EPSG9848, // CC48
    EPSG9849, // CC49
    EPSG9850, // CC50
    EPSG3395, // mercator
};
 
class CscCartographicProjection {
public:
    virtual ~CscCartographicProjection();
    
    CscPoint3d *projectPointFromWGS84(const string &target, const CscPoint3d *position);
    CscPoint3d *projectPointToWGS84(const string &src, const CscPoint3d *position);
    CscPoint3d *projectPointFromWGS84ToConsciencePosition(const CscPoint3d *position);
    CscPoint3d *projectPointFromConsciencePositionToWGS84(const CscPoint3d *position);
    CscCartographicProjectionType getEpsgProjectionFromString(string epsgCode);
    CscPoint3d *projectPointFromWGS84WithReferencePoint(const string &target, const CscPoint3d *position, const CscPoint3d *referencePoint);
    CscPoint3d *projectPointToWGS84WithReferencePoint(const string &src, const CscPoint3d *position, const CscPoint3d *referencePoint);
 
    static CscCartographicProjection *getInstance() {
        if (instance == nullptr) {
            instance = new CscCartographicProjection();
        }
        return instance;
    }
 
    inline static constexpr const char* EPSG3395_DEF = "EPSG3395";
    inline static constexpr const char* EPSG9849_DEF = "EPSG9849";
    inline static constexpr const char* ESRI54032_DEF = "ESRI54032";
 
    inline static constexpr const char* DEFAULT_PROJECTION_DEF = EPSG3395_DEF;
    inline static const CscCartographicProjectionType DEFAULT_PROJECTION = CscCartographicProjectionType::EPSG3395;
 
 
    inline static constexpr const char* FRANCE_PROJECTION_DEF = EPSG9849_DEF;
    inline static const CscCartographicProjectionType FRANCE_PROJECTION = CscCartographicProjectionType::EPSG9849;
 
private:
    CscCartographicProjection();
    // Position in Mercator of ConscienceLabs, its use for reduce number of digits in the position
    const long double latRef = 627905946.399036;
    const long double lngRef = -3945574.63582955;
    inline static CscCartographicProjection *instance = nullptr;
    map<CscCartographicProjectionType, string> cartographicProjectionType;
    map<string, CscCartographicProjectionType> cartographicProjectionTypeName;
    const string getProjectionDefinition(const string &projectionType);
    const string wgs84 = "+proj=longlat +datum=WGS84 +no_defs +type=crs";
    mutex projMutex;
};
 
}
 
#endif