390 lines
16 KiB
C++
390 lines
16 KiB
C++
#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);
|
|
|
|
enableComputeInFloat(true);
|
|
}
|
|
//------------------------------------------------------------
|
|
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()) &&
|
|
!((TRasterFP)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;
|
|
}
|
|
}
|