#include <sstream> /* std::ostringstream */
#include "tfxparam.h"
#include "tparamset.h"  // TPointParamP
#include "stdfx.h"
#include "tparamuiconcept.h"  //add 20140130

#include "ino_common.h"
#include "igs_radial_blur.h"

//------------------------------------------------------------
namespace {
const double ino_param_range = 100.0;
}
class ino_radial_blur final : public TStandardRasterFx {
  FX_PLUGIN_DECLARATION(ino_radial_blur)
  TRasterFxPort m_input;
  TRasterFxPort m_refer;

  TPointParamP m_center;
  TDoubleParamP m_blur;
  TDoubleParamP m_radius;
  TDoubleParamP m_twist;
  // TDoubleParamP m_twist_radius;
  TBoolParamP m_alpha_rendering;
  TBoolParamP m_anti_alias;
  TIntEnumParamP m_ref_mode;
  TIntEnumParamP m_type;
  // elliptical shape
  TDoubleParamP m_ellipse_aspect_ratio;
  TDoubleParamP m_ellipse_angle;
  TDoubleParamP m_intensity_correlation_with_ellipse;

public:
  ino_radial_blur()
      : m_center(TPointD(0.0, 0.0))
      , m_blur(0.01 * ino_param_range)
      , m_radius(0.0)
      , m_twist(0.0)
      // ,m_twist_radius(0.0)
      , m_alpha_rendering(true)
      , m_anti_alias(false)
      , m_ref_mode(new TIntEnumParam(0, "Red"))
      , m_type(new TIntEnumParam(0, "Accelerator"))
      , m_ellipse_aspect_ratio(1.0)
      , m_ellipse_angle(0.0)
      , m_intensity_correlation_with_ellipse(0.0) {
    this->m_center->getX()->setMeasureName("fxLength");
    this->m_center->getY()->setMeasureName("fxLength");
    this->m_radius->setMeasureName("fxLength");

    addInputPort("Source", this->m_input);
    addInputPort("Reference", this->m_refer);

    bindParam(this, "center", this->m_center);
    bindParam(this, "blur", this->m_blur);
    bindParam(this, "radius", this->m_radius);
    bindParam(this, "twist", this->m_twist);
    // bindParam(this,"twist_radius", this->m_twist_radius);
    bindParam(this, "alpha_rendering", this->m_alpha_rendering);
    bindParam(this, "anti_alias", this->m_anti_alias);
    bindParam(this, "reference", this->m_ref_mode);
    bindParam(this, "type", this->m_type);
    bindParam(this, "ellipse_aspect_ratio", this->m_ellipse_aspect_ratio);
    bindParam(this, "ellipse_angle", this->m_ellipse_angle);
    bindParam(this, "intensity_correlation_with_ellipse",
              this->m_intensity_correlation_with_ellipse);

    this->m_radius->setValueRange(0, (std::numeric_limits<double>::max)());
    /* 拡大のしすぎを防ぐためにMaxを制限する */
    this->m_blur->setValueRange(0.0 * ino_param_range, 1.0 * ino_param_range);
    this->m_ellipse_aspect_ratio->setValueRange(0.1, 10.0);
    this->m_ellipse_angle->setValueRange(-180.0, 180.0);
    this->m_twist->setValueRange(-180.0, 180.0);
    // this->m_twist_radius->setValueRange(0.0,1000.0);
    this->m_ref_mode->addItem(1, "Green");
    this->m_ref_mode->addItem(2, "Blue");
    this->m_ref_mode->addItem(3, "Alpha");
    this->m_ref_mode->addItem(4, "Luminance");
    this->m_ref_mode->addItem(-1, "Nothing");
    this->m_type->addItem(1, "Uniform Length");
    this->m_intensity_correlation_with_ellipse->setValueRange(-1.0, 1.0);
  }
  //------------------------------------------------------------
  TPointD get_render_center(const double frame, const TPointD &pos,
                            const TAffine affine) {
    /*--- 表示画像の中心をゼロとしたミリメートル単位の位置 ---*/
    TPointD center = this->m_center->getValue(frame);
    /*--- 単位変換(mm --> render用pixel) ---*/
    center = center * ino::pixel_per_mm();
    /*--- Affine変換(Geometry&Shrink&原点を中心から左下へ) ---*/
    center = affine * center;
    /*--- カメラ座標をrender画像座標の原点に ---*/
    return center - pos;
  }
  int get_render_int_margin(const double frame, const TRectD &bBox,
                            const TAffine affine, TPointD center,
                            double blur_radius, double pivot_radius) {
    /*--- 単位変換(mm-->render用pixel)と長さ(scale)変換 ---*/
    const double scale = ino::pixel_per_mm() * sqrt(fabs(affine.det()));
    /*--- margin計算...Twist時正確でない ---*/
    return igs::radial_blur::reference_margin(
        static_cast<int>(ceil(bBox.getLy())),
        static_cast<int>(ceil(bBox.getLx())), center,
        this->m_twist->getValue(frame) * M_PI_180, pivot_radius,
        this->m_blur->getValue(frame) / ino_param_range, blur_radius,
        this->m_type->getValue(), this->m_ellipse_aspect_ratio->getValue(frame),
        this->m_ellipse_angle->getValue(frame),
        this->m_intensity_correlation_with_ellipse->getValue(frame));
  }
  void get_render_enlarge(const double frame, const TAffine affine,
                          TRectD &bBox, double blur_radius,
                          double pivot_radius) {
    TPointD center(this->get_render_center(frame, bBox.getP00(), affine));
    int margin = this->get_render_int_margin(frame, bBox, affine, center,
                                             blur_radius, pivot_radius);
    if (0 < margin) {
      /* 拡大のしすぎを防ぐテキトーな制限 */
      if (4096 < margin) {
        margin = 4096;
      }
      bBox = bBox.enlarge(margin);
    }
  }
  //------------------------------------------------------------
  bool doGetBBox(double frame, TRectD &bBox,
                 const TRenderSettings &info) override {
    if (!this->m_input.isConnected()) {
      bBox = TRectD();
      return false;
    }
    const bool ret      = this->m_input->doGetBBox(frame, bBox, info);
    const double scale  = ino::pixel_per_mm() * sqrt(fabs(info.m_affine.det()));
    const double radius = this->m_radius->getValue(frame) * scale;
    const double pivot_radius = info.m_cameraBox.getLy() / 2.0;
    this->get_render_enlarge(frame, info.m_affine, bBox, radius, pivot_radius);
    return ret;
  }
  int getMemoryRequirement(const TRectD &rect, double frame,
                           const TRenderSettings &info) override {
    TRectD bBox(rect);
    const double scale  = ino::pixel_per_mm() * sqrt(fabs(info.m_affine.det()));
    const double radius = this->m_radius->getValue(frame) * scale;
    const double pivot_radius = info.m_cameraBox.getLy() / 2.0;
    this->get_render_enlarge(frame, info.m_affine, bBox, radius, pivot_radius);
    return TRasterFx::memorySize(bBox, info.m_bpp);
  }
  void transform(double frame, int port, const TRectD &rectOnOutput,
                 const TRenderSettings &infoOnOutput, TRectD &rectOnInput,
                 TRenderSettings &infoOnInput) override {
    rectOnInput = rectOnOutput;
    infoOnInput = infoOnOutput;
    const double scale =
        ino::pixel_per_mm() * sqrt(fabs(infoOnInput.m_affine.det()));
    const double radius       = this->m_radius->getValue(frame) * scale;
    const double pivot_radius = infoOnInput.m_cameraBox.getLy() / 2.0;
    this->get_render_enlarge(frame, infoOnOutput.m_affine, rectOnInput, radius,
                             pivot_radius);
  }
  bool canHandle(const TRenderSettings &info, double frame) override {
    // return false; // toonz has geometry control
    // return true;
    /* 2012-11-14:
     * trueのとき、Fx処理がジオメトリ変換の最後にまとめられてしまう。 */
    return m_blur->getValue(frame) == 0 ? true
                                        : isAlmostIsotropic(info.m_affine);
  }
  void doCompute(TTile &tile, double frame, const TRenderSettings &ri) override;

  // add 20140130
  void getParamUIs(TParamUIConcept *&concepts, int &length) override {
    concepts            = new TParamUIConcept[length = 2];
    concepts[0].m_type  = TParamUIConcept::ELLIPSE;
    concepts[0].m_label = "Radius";
    concepts[0].m_params.push_back(m_radius);
    concepts[0].m_params.push_back(m_center);
    concepts[0].m_params.push_back(m_ellipse_aspect_ratio);
    concepts[0].m_params.push_back(m_ellipse_angle);
    concepts[0].m_params.push_back(m_twist);

    concepts[1].m_type = TParamUIConcept::COMPASS;
    concepts[1].m_params.push_back(m_center);
    concepts[1].m_params.push_back(m_ellipse_aspect_ratio);
    concepts[1].m_params.push_back(m_ellipse_angle);
    concepts[1].m_params.push_back(m_twist);
  }
  // add 20140130
};
FX_PLUGIN_IDENTIFIER(ino_radial_blur, "inoRadialBlurFx");
//--------------------------------------------------------------------
namespace {
void fx_(const TRasterP in_ras,  // with margin
         const int margin, /* ここではinとrefは同サイズで使用している */
         const TRasterP refer_ras,  // no margin
         const int refer_mode,
         TRasterP out_ras,  // no margin
         const TPointD center, const double twist, const double blur,
         const double radius, const double pivot_radius,
         const bool alpha_rendering_sw, const bool anti_alias_sw,
         const int type, const double ellipse_aspect_ratio,
         const double ellipse_angle,
         const double intensity_correlation_with_ellipse) {
  TRasterGR8P ref_gr8;
  if ((refer_ras != nullptr) && (0 <= refer_mode)) {
    ref_gr8 = TRasterGR8P(out_ras->getLy(), out_ras->getLx() * sizeof(float));
    ref_gr8->lock();
    ino::ras_to_ref_float_arr(refer_ras,
                              reinterpret_cast<float *>(ref_gr8->getRawData()),
                              refer_mode);
  }

  TRasterGR8P in_gr8(in_ras->getLy(),
                     in_ras->getLx() * ino::channels() * sizeof(float));
  in_gr8->lock();
  ino::ras_to_float_arr(in_ras, ino::channels(),
                        reinterpret_cast<float *>(in_gr8->getRawData()));

  TRasterGR8P out_buffer(out_ras->getLy(),
                         out_ras->getLx() * ino::channels() * sizeof(float));
  out_buffer->lock();

  igs::radial_blur::convert(
      reinterpret_cast<float *>(in_gr8->getRawData()),
      reinterpret_cast<float *>(out_buffer->getRawData()), margin,
      out_ras->getSize(), ino::channels(),
      (ref_gr8) ? reinterpret_cast<float *>(ref_gr8->getRawData()) : nullptr,
      center + TPointD(margin, margin), twist, pivot_radius, blur, radius, type,
      anti_alias_sw, alpha_rendering_sw, ellipse_aspect_ratio, ellipse_angle,
      intensity_correlation_with_ellipse);

  in_gr8->unlock();
  ino::float_arr_to_ras(out_buffer->getRawData(), ino::channels(), out_ras, 0);
  out_buffer->unlock();
  if (ref_gr8) ref_gr8->unlock();
}
}  // namespace
//------------------------------------------------------------
void ino_radial_blur::doCompute(TTile &tile, double frame,
                                const TRenderSettings &ri) {
  /*------ 接続していなければ処理しない ----------------------*/
  if (!this->m_input.isConnected()) {
    tile.getRaster()->clear(); /* 塗りつぶしクリア */
    return;
  }
  /*------ サポートしていないPixelタイプはエラーを投げる -----*/
  if (!((TRaster32P)tile.getRaster()) && !((TRaster64P)tile.getRaster())) {
    throw TRopException("unsupported input pixel type");
  }
  /*------ パラメータを得る ----------------------------------*/
  const double scale  = ino::pixel_per_mm() * sqrt(fabs(ri.m_affine.det()));
  const double blur   = this->m_blur->getValue(frame) / ino_param_range;
  const double radius = this->m_radius->getValue(frame) * scale;
  const double twist =
      this->m_twist->getValue(frame) * 3.14159265358979 / 180.0;
  // const double twist_radius = this->m_twist_radius->getValue(frame)*scale;
  const bool alpha_rend_sw = this->m_alpha_rendering->getValue();
  const bool anti_alias_sw = this->m_anti_alias->getValue();
  const int refer_mode     = this->m_ref_mode->getValue();
  const int type           = this->m_type->getValue();
  const double ellipse_aspect_ratio =
      this->m_ellipse_aspect_ratio->getValue(frame);
  const double ellipse_angle = this->m_ellipse_angle->getValue(frame);
  const double pivot_radius  = ri.m_cameraBox.getLy() / 2.0;
  const double intensity_correlation_with_ellipse =
      this->m_intensity_correlation_with_ellipse->getValue(frame);

  TPointD center = this->m_center->getValue(frame);
  TPointD render_center(
      this->get_render_center(frame, tile.m_pos, ri.m_affine));
  /*------ 入力画像の範囲を得る ------------------------------*/
  TRectD tile_with_margin =
      TRectD(tile.m_pos /* Render画像上(Pixel単位)の位置 */
             ,
             TDimensionD(/* Render画像上(Pixel単位)のサイズ */
                         tile.getRaster()->getLx(), tile.getRaster()->getLy()));
  /*------ 入力画像のエリアの拡大 ----------------------------*/
  int margin = this->get_render_int_margin(frame, tile_with_margin, ri.m_affine,
                                           render_center, radius, pivot_radius);
  if (0 < margin) {
    /* 拡大のしすぎを防ぐテキトーな制限 */
    if (4096 < margin) {
      margin = 4096;
    }
    tile_with_margin = tile_with_margin.enlarge(margin);
  }
  /*------ 入力画像生成 --------------------------------------*/
  TTile enlarge_tile;
  this->m_input->allocateAndCompute(
      enlarge_tile, tile_with_margin.getP00(),
      TDimensionI(/* Pixel単位で四捨五入 */
                  static_cast<int>(tile_with_margin.getLx() + 0.5),
                  static_cast<int>(tile_with_margin.getLy() + 0.5)),
      tile.getRaster(), frame, ri);
  /*------ 参照画像生成 --------------------------------------*/
  TTile refer_tile;
  bool refer_sw = false;
  if (this->m_refer.isConnected()) {
    refer_sw = true;
    this->m_refer->allocateAndCompute(
        refer_tile, tile.m_pos,
        TDimension(tile.getRaster()->getLx(), tile.getRaster()->getLy()),
        tile.getRaster(), frame, ri);
  }
  /* ------ 保存すべき画像メモリを塗りつぶしクリア ---------- */
  tile.getRaster()->clear(); /* 塗りつぶしクリア */
  /* ------ (app_begin)log記憶 ------------------------------ */
  const bool log_sw = ino::log_enable_sw();

  if (log_sw) {
    std::ostringstream os;
    os << "params"
       << "  cx " << center.x << "  cy " << center.y << "  blur " << blur
       << "  radius " << radius << "  twist "
       << twist
       // << "  twist_radius " << twist_radius
       << "  reference " << refer_mode << "  alpha_rendering " << alpha_rend_sw
       << "  anti_alias " << anti_alias_sw << "  render_center "
       << render_center << "  frame " << frame << "  pixbits "
       << ino::pixel_bits(tile.getRaster()) << "  tile.m_pos " << tile.m_pos
       << "  tile_getLx " << tile.getRaster()->getLx() << "  y "
       << tile.getRaster()->getLy() << "  ri.m_cameraBox " << ri.m_cameraBox
       << "  ri.m_affine " << ri.m_affine << "  ri_m_affine_det "
       << ri.m_affine.det() << "  shrink_x " << ri.m_shrinkX << "  shrink_y "
       << ri.m_shrinkY << "  tile_with_margin " << tile_with_margin
       << "  enlarge_tile.m_pos " << enlarge_tile.m_pos
       << "  tile_allocated_getLx " << enlarge_tile.getRaster()->getLx()
       << "  y " << enlarge_tile.getRaster()->getLy();
    if (refer_sw) {
      os << "  refer_tile.m_pos " << refer_tile.m_pos << "  refer_tile_getLx "
         << refer_tile.getRaster()->getLx() << "  y "
         << refer_tile.getRaster()->getLy();
    }
  }
  /* ------ fx処理 ------------------------------------------ */
  try {
    tile.getRaster()->lock();
    enlarge_tile.getRaster()->lock();
    if (refer_tile.getRaster() != nullptr) {
      refer_tile.getRaster()->lock();
    }
    fx_(enlarge_tile.getRaster(), margin, refer_tile.getRaster(), refer_mode,
        tile.getRaster(), render_center, twist, blur, radius, pivot_radius,
        alpha_rend_sw, anti_alias_sw, type, ellipse_aspect_ratio, ellipse_angle,
        intensity_correlation_with_ellipse);
    if (refer_tile.getRaster() != nullptr) {
      refer_tile.getRaster()->unlock();
    }
    enlarge_tile.getRaster()->unlock();
    tile.getRaster()->unlock();
  }
  /* ------ error処理 --------------------------------------- */
  catch (std::bad_alloc &e) {
    if (refer_tile.getRaster() != nullptr) {
      refer_tile.getRaster()->unlock();
    }
    enlarge_tile.getRaster()->unlock();
    tile.getRaster()->unlock();
    if (log_sw) {
      std::string str("std::bad_alloc <");
      str += e.what();
      str += '>';
    }
    throw;
  } catch (std::exception &e) {
    if (refer_tile.getRaster() != nullptr) {
      refer_tile.getRaster()->unlock();
    }
    enlarge_tile.getRaster()->unlock();
    tile.getRaster()->unlock();
    if (log_sw) {
      std::string str("exception <");
      str += e.what();
      str += '>';
    }
    throw;
  } catch (...) {
    if (refer_tile.getRaster() != nullptr) {
      refer_tile.getRaster()->unlock();
    }
    enlarge_tile.getRaster()->unlock();
    tile.getRaster()->unlock();
    if (log_sw) {
      std::string str("other exception");
    }
    throw;
  }
}