channel_transform.h

Go to the documentation of this file.

#pragma once

 
#include "core/exception.h"
#include "nlohmann/json.hpp"
 
#include "projection/thinplatespline.h" // TODOREWORK Move to here?
 
namespace satdump
{
    class ChannelTransform
    {
    public:
        enum transform_type_t
        {
            TYPE_INVALID = 0,
            TYPE_NONE = 1,
            TYPE_AFFINE = 2,
            TYPE_AFFINE_SLANTX = 3,
            TYPE_AFFINE_INTERPX = 4,
            TYPE_INTERP_XY = 5,
        };
 
    private:
        transform_type_t d_type = TYPE_INVALID;
 
        double affine_a_x;
        double affine_a_y;
        double affine_b_x;
        double affine_b_y;
 
        double slantx_center;
        double slantx_ampli;
 
        std::vector<std::pair<double, double>> interpx_points;
        std::vector<std::pair<std::pair<double, double>, std::pair<double, double>>> interp_xy_points;
        std::shared_ptr<projection::VizGeorefSpline2D> interp_fwd_interpolator;
        std::shared_ptr<projection::VizGeorefSpline2D> interp_rev_interpolator;
 
    public:
        ChannelTransform &init_none();

        ChannelTransform &init_affine(double a_x, double a_y, double b_x, double b_y);

        ChannelTransform &init_affine_slantx(double a_x, double a_y, double b_x, double b_y, double slantx_center, double slantx_ampli);

        ChannelTransform &init_affine_interpx(double a_x, double a_y, double b_x, double b_y, std::vector<std::pair<double, double>> interpx_points);
 
        // TODOREWORK
        ChannelTransform &init_interp_xy(std::vector<std::pair<std::pair<double, double>, std::pair<double, double>>> interp_xy_points);
 
    public:
        friend void to_json(nlohmann::json &j, const ChannelTransform &v)
        {
            j["type"] = (int)v.d_type;
 
            if (v.d_type == TYPE_AFFINE || v.d_type == TYPE_AFFINE_SLANTX || v.d_type == TYPE_AFFINE_INTERPX)
            {
                j["ax"] = v.affine_a_x;
                j["ay"] = v.affine_a_y;
                j["bx"] = v.affine_b_x;
                j["by"] = v.affine_b_y;
            }
 
            if (v.d_type == TYPE_AFFINE_SLANTX)
            {
                j["slantx_center"] = v.slantx_center;
                j["slantx_ampli"] = v.slantx_ampli;
            }
 
            if (v.d_type == TYPE_AFFINE_INTERPX)
            {
                j["interpx_points"] = v.interpx_points; // v.fwd_interpolator.get_points();
            }
 
            if (v.d_type == TYPE_INTERP_XY)
            {
                j["interp_xy_points"] = v.interp_xy_points; // v.fwd_interpolator.get_points();
            }
        }
 
        friend void from_json(const nlohmann::json &j, ChannelTransform &v)
        {
            v.d_type = (transform_type_t)j["type"].get<int>();
 
            if (v.d_type == TYPE_AFFINE || v.d_type == TYPE_AFFINE_SLANTX || v.d_type == TYPE_AFFINE_INTERPX)
            {
                v.affine_a_x = j["ax"];
                v.affine_a_y = j["ay"];
                v.affine_b_x = j["bx"];
                v.affine_b_y = j["by"];
            }
 
            if (v.d_type == TYPE_AFFINE_SLANTX)
            {
                v.slantx_center = j["slantx_center"];
                v.slantx_ampli = j["slantx_ampli"];
            }
 
            if (v.d_type == TYPE_AFFINE_INTERPX)
            {
                v.interpx_points = j["interpx_points"].get<std::vector<std::pair<double, double>>>();
 
                v.interp_fwd_interpolator = std::make_shared<projection::VizGeorefSpline2D>(1);
                v.interp_rev_interpolator = std::make_shared<projection::VizGeorefSpline2D>(1);
 
                for (auto &p : v.interpx_points)
                {
                    v.interp_fwd_interpolator->add_point(p.first, p.first, &p.second);
                    v.interp_rev_interpolator->add_point(p.second, p.second, &p.first);
                }
 
                v.interp_fwd_interpolator->solve();
                v.interp_rev_interpolator->solve();
            }
 
            if (v.d_type == TYPE_INTERP_XY)
            {
                v.interp_xy_points = j["interp_xy_points"].get<std::vector<std::pair<std::pair<double, double>, std::pair<double, double>>>>();
 
                v.interp_fwd_interpolator = std::make_shared<projection::VizGeorefSpline2D>(2);
                v.interp_rev_interpolator = std::make_shared<projection::VizGeorefSpline2D>(2);
 
                for (auto &p : v.interp_xy_points)
                {
                    printf("Point %f %f => %f %f\n", p.first.first, p.first.second, p.second.first, p.second.second);
                    double p1[2] = {p.second.first, p.second.second};
                    v.interp_fwd_interpolator->add_point(p.first.first, p.first.second, p1);
                    double p2[2] = {p.first.first, p.first.second};
                    v.interp_rev_interpolator->add_point(p.second.first, p.second.second, p2);
                }
 
                printf("%d %d\n", v.interp_fwd_interpolator->solve(), v.interp_rev_interpolator->solve());
            }
        }

        transform_type_t getType();

        void render();
 
    public: // Must be inline for performance reasons!
        inline void forward(double *x, double *y) const
        {
            if (d_type == TYPE_NONE)
                ;
            else if (d_type == TYPE_AFFINE)
            {
                *x = *x * affine_a_x + affine_b_x;
                *y = *y * affine_a_y + affine_b_y;
            }
            else if (d_type == TYPE_AFFINE_SLANTX)
            {
                *y -= (*x - slantx_center) * slantx_ampli;
                *x = *x * affine_a_x + affine_b_x;
                *y = *y * affine_a_y + affine_b_y;
            }
            else if (d_type == TYPE_AFFINE_INTERPX)
            {
                *x = *x * affine_a_x + affine_b_x;
                *y = *y * affine_a_y + affine_b_y;
                interp_fwd_interpolator->get_point(*x, *x, x);
            }
            else if (d_type == TYPE_INTERP_XY)
            {
                double out[2]; // TODOREWORK
              //  printf("%f %f FWD => ", *x, *y);
                interp_fwd_interpolator->get_point(*x, *y, out);
                *x = out[0];
                *y = out[1];
              //  printf("%f %f\n ", *x, *y);
            }
            else if (d_type == TYPE_INVALID)
                throw satdump_exception("Invalid Channel Transform!\n");
        }

        inline void reverse(double *x, double *y) const
        {
            if (d_type == TYPE_NONE)
                ;
            else if (d_type == TYPE_AFFINE)
            {
                *x = (*x - affine_b_x) / affine_a_x;
                *y = (*y - affine_b_y) / affine_a_y;
            }
            else if (d_type == TYPE_AFFINE_SLANTX)
            {
                *x = (*x - affine_b_x) / affine_a_x;
                *y = (*y - affine_b_y) / affine_a_y;
                *y += (*x - slantx_center) * slantx_ampli;
            }
            else if (d_type == TYPE_AFFINE_INTERPX)
            {
                interp_rev_interpolator->get_point(*x, *x, x);
                *x = (*x - affine_b_x) / affine_a_x;
                *y = (*y - affine_b_y) / affine_a_y;
            }
            else if (d_type == TYPE_INTERP_XY)
            {
                double out[2]; // TODOREWORK
                printf("%f %f REV => ", *x, *y);
                interp_rev_interpolator->get_point(*x, *y, out);
                *x = out[0];
                *y = out[1];
                printf("%f %f\n ", *x, *y);
            }
            else if (d_type == TYPE_INVALID)
                throw satdump_exception("Invalid Channel Transform!\n");
        }
    };
} // namespace satdump