#pragma once

#include <cstdint>
#include <cmath>
#include <array>
#include <vector>
#include "image/image.h"
#include "common/calibration.h"
#include "nlohmann/json.hpp"

// #include <fstream>

#include "common/ccsds/ccsds.h"
#include "common/ccsds/ccsds_time.h"

#define SCI_PACKET_SIZE 1286
#define MIU_BYTE_OFFSET 48
#define MHS_OFFSET 49
#define MHS_WIDTH 90
#define PRT_OFFSET 1225
#define HKTH_offset 9

#define DAY_OFFSET 17453
#define SEC_OFFSET 3600 * 9 - 60 * 10

#define CAL_LIMIT 5000

namespace noaa_metop
{
    namespace mhs
    {
        class MHSReader
        {
        private:
            struct calib_line
            {
                std::array<uint16_t, 3> PRT_calib;
                std::array<uint16_t, 5> PRT_readings;
                std::array<uint8_t, 39> HK;
                std::array<std::array<uint16_t, 2>, 5> calibration_views; // space 0, warm 1
            };

        private:
            std::array<std::vector<std::array<uint16_t, MHS_WIDTH>>, 5> channels;
            std::array<std::vector<std::array<double, MHS_WIDTH>>, 5> calibrated_channels;
            void work(uint8_t *buffer);
            // std::ofstream deb_out;

            // NOAA specific stuff
            unsigned int last = 0;
            std::array<uint8_t, 50> MIU_data[80];
            uint32_t major_cycle_count = 0;
            uint32_t last_major_cycle = 0;

            // calib values
            nlohmann::json calib;
            std::vector<calib_line> calib_lines;
            std::vector<uint8_t> PIE_buff;
            double last_Tw = 0;

            uint16_t limits[2][2] = {{1500, 28000}, {20000, 64881}};

            // calib functions
            std::array<uint8_t, SCI_PACKET_SIZE> get_SCI_packet(int PKT);
            double get_u(double temp, int ch);
            double interpolate(double a1x, double a1y, double a2x, double a2y, double bx, int mode);
            double get_timestamp(int pkt, int offset, int ms_scale = 1000);

        public:
            MHSReader();
            int line = 0;
            std::vector<double> timestamps;
            void work_NOAA(uint8_t *buffer);
            void work_metop(ccsds::CCSDSPacket &packet);
            image::Image getChannel(int channel);
            nlohmann::json calib_out;
            void calibrate(nlohmann::json calib_coefs);
            nlohmann::json dump_telemetry(nlohmann::json calib_coefs);
        };
    } // namespace hirs
} // namespace noaa