#include "trandom.h"
#include "tfxparam.h"
#include "stdfx.h"

class NoiseFx : public TStandardRasterFx {
  FX_PLUGIN_DECLARATION(NoiseFx)

  TRasterFxPort m_input;
  TDoubleParamP m_value;
  TBoolParamP m_Red;
  TBoolParamP m_Green;
  TBoolParamP m_Blue;
  TBoolParamP m_BW;
  TBoolParamP m_Animate;

public:
  NoiseFx()
      : m_value(100.0)
      , m_Red(1.0)
      , m_Green(1.0)
      , m_Blue(1.0)
      , m_BW(0.0)
      , m_Animate(0.0) {
    bindParam(this, "Intensity", m_value);
    bindParam(this, "Red", m_Red);
    bindParam(this, "Green", m_Green);
    bindParam(this, "Blue", m_Blue);
    bindParam(this, "Black_White", m_BW);
    bindParam(this, "Animate", m_Animate);
    addInputPort("Source", m_input);
    m_value->setValueRange(0, (std::numeric_limits<double>::max)());
  }

  ~NoiseFx(){};

  bool doGetBBox(double frame, TRectD &bBox, const TRenderSettings &info) {
    if (m_input.isConnected())
      return m_input->doGetBBox(frame, bBox, info);
    else {
      bBox = TRectD();
      return false;
    }
  }

  void doCompute(TTile &tile, double frame, const TRenderSettings &ri);
  std::string getAlias(double frame, const TRenderSettings &info) const;

  bool canHandle(const TRenderSettings &info, double frame) { return true; }
};

//------------------------------------------------------------------------------

namespace {
template <class PIXEL>
class RayleighNoise {
  std::vector<double> noise_buf;

public:
  RayleighNoise(double sigma, int seed = 0) {
    noise_buf.resize(PIXEL::maxChannelValue + 1);
    TRandom random(seed);
    sigma = log(sigma * 0.07 + 1);
    if (PIXEL::maxChannelValue == 255)
      sigma = 2.0 * sigma * sigma * sigma * sigma;
    else
      sigma = 257.0 * 2.0 * sigma * sigma * sigma * sigma;
    for (int i = 0; i < PIXEL::maxChannelValue + 1; i++) {
      noise_buf[i] = sigma * sqrt(-log(1.0 - random.getFloat())) *
                     cos(M_PI * (2.0 * random.getFloat() - 1.0));
    }
  }
  double get_value(int index) {
    tcrop(index, 0, PIXEL::maxChannelValue);
    return noise_buf[index];
  }
};
}

template <typename PIXEL, typename PIXELGRAY, typename CHANNEL_TYPE>
void doNoise(TRasterPT<PIXEL> &ras, double sigma, bool bw, bool red, bool green,
             bool blue, bool animate, double frame) {
  RayleighNoise<PIXEL> noise(sigma);

  TRandom random;
  if (animate) random.setSeed(frame);
  int j;
  ras->lock();
  for (j = 0; j < ras->getLy(); j++) {
    PIXEL *pix    = ras->pixels(j);
    PIXEL *endPix = pix + ras->getLx();
    while (pix < endPix) {
      if (bw) {
        double index;
        double value;
        index = random.getFloat() * pix->m;  // PIXEL::maxChannelValue;
        value = noise.get_value(tfloor(index));

        // value =
        // (value<pix->m)?((value>0)?(int)value:0):PIXEL::maxChannelValue;
        // if(value)
        //{
        int reference;
        reference = PIXELGRAY::from(*pix).value;
        pix->r    = tcrop<int>(value + reference, 0, pix->m);
        pix->g    = tcrop<int>(value + reference, 0, pix->m);
        pix->b    = tcrop<int>(value + reference, 0, pix->m);
        //}
      } else {
        double r, g, b, index;
        if (red) {
          index  = random.getFloat() * pix->m;  // PIXEL::maxChannelValue;
          r      = noise.get_value(tfloor(index)) + pix->r;
          pix->r = tcrop<int>(r, 0, pix->m);
        }
        if (green) {
          index  = random.getFloat() * pix->m;  // PIXEL::maxChannelValue;
          g      = noise.get_value(tfloor(index)) + pix->g;
          pix->g = tcrop<int>(g, 0, pix->m);
        }
        if (blue) {
          index  = random.getFloat() * pix->m;  // PIXEL::maxChannelValue;
          b      = noise.get_value(tfloor(index)) + pix->b;
          pix->b = tcrop<int>(b, 0, pix->m);
        }
      }
      *pix++;
    }
  }

  ras->unlock();
}

//-------------------------------------------------------------------

void NoiseFx::doCompute(TTile &tile, double frame, const TRenderSettings &ri) {
  if (!m_input.isConnected()) return;

  m_input->compute(tile, frame, ri);

  double sigma = m_value->getValue(frame);
  bool bw      = m_BW->getValue();
  bool red     = m_Red->getValue();
  bool green   = m_Green->getValue();
  bool blue    = m_Blue->getValue();
  bool animate = m_Animate->getValue();
  if (sigma == 0.0) return;

  TRaster32P raster32 = tile.getRaster();
  if (raster32)
    doNoise<TPixel32, TPixelGR8, UCHAR>(raster32, sigma, bw, red, green, blue,
                                        animate, frame);
  else {
    TRaster64P raster64 = tile.getRaster();
    if (raster64)
      doNoise<TPixel64, TPixelGR16, USHORT>(raster64, sigma, bw, red, green,
                                            blue, animate, frame);
    else
      throw TException("Brightness&Contrast: unsupported Pixel Type");
  }
}
//------------------------------------------------------------------

std::string NoiseFx::getAlias(double frame, const TRenderSettings &info) const {
  std::string alias = getFxType();
  alias += "[";

  // alias degli effetti connessi alle porte di input separati da virgole
  // una porta non connessa da luogo a un alias vuoto (stringa vuota)
  int i = 0;

  for (i = 0; i < getInputPortCount(); i++) {
    TFxPort *port = getInputPort(i);
    if (port->isConnected()) {
      TRasterFxP ifx = port->getFx();
      assert(ifx);
      alias += ifx->getAlias(frame, info);
    }
    alias += ",";
  }

  bool addframe = 0;
  std::string paramalias("");
  for (i = 0; i < getParams()->getParamCount(); i++) {
    TParam *param = getParams()->getParam(i);
    paramalias += param->getName() + "=" + param->getValueAlias(frame, 3);
    if (param->getName() == "Animate" && param->getValueAlias(frame, 0) == "1")
      addframe = true;
  }

  if (addframe) alias += std::to_string(frame) + ",";
  alias += paramalias + "]";

  return alias;
}

//------------------------------------------------------------------

FX_PLUGIN_IDENTIFIER(NoiseFx, "noiseFx")