static void
sinus (void)
{
params p;
gint bytes;
GimpPixelRgn dest_rgn;
gint x1, y1, x2, y2;
gpointer pr;
gint progress, max_progress;
prepare_coef(&p);
gimp_drawable_mask_bounds(drawable->drawable_id, &x1, &y1, &x2, &y2);
p.width = drawable->width;
p.height = drawable->height;
bytes = drawable->bpp;
gimp_pixel_rgn_init (&dest_rgn, drawable,
x1, y1, x2 - x1, y2 - y1, TRUE,TRUE);
progress = 0;
max_progress = (x2 - x1) * (y2 - y1);
for (pr = gimp_pixel_rgns_register (1, &dest_rgn);
pr != NULL;
pr = gimp_pixel_rgns_process (pr))
{
switch (bytes)
{
case 4:
compute_block_x (dest_rgn.data, dest_rgn.rowstride,
dest_rgn.x, dest_rgn.y, dest_rgn.w, dest_rgn.h,
4, assign_block_4, &p);
break;
case 3:
compute_block_x (dest_rgn.data, dest_rgn.rowstride,
dest_rgn.x, dest_rgn.y, dest_rgn.w, dest_rgn.h,
3, assign_block_3, &p);
break;
case 2:
compute_block_x (dest_rgn.data, dest_rgn.rowstride,
dest_rgn.x, dest_rgn.y, dest_rgn.w, dest_rgn.h,
2, assign_block_2, &p);
break;
case 1:
compute_block_x (dest_rgn.data, dest_rgn.rowstride,
dest_rgn.x, dest_rgn.y, dest_rgn.w, dest_rgn.h,
1, assign_block_1, &p);
break;
}
progress += dest_rgn.w * dest_rgn.h;
gimp_progress_update ((double) progress / (double) max_progress);
}
gimp_drawable_flush (drawable);
gimp_drawable_merge_shadow (drawable->drawable_id, TRUE);
gimp_drawable_update (drawable->drawable_id, x1, y1, x2 - x1, y2 - y1);
}
void
sinus_do_preview (GtkWidget *widget)
{
static GtkWidget *theWidget = NULL;
gint rowsize;
guchar *buf;
params p;
if (!do_preview)
return;
if (theWidget == NULL)
{
theWidget = widget;
}
rowsize = thePreview->width * thePreview->bpp;
buf = g_new (guchar, thePreview->width*thePreview->height*thePreview->bpp);
p.height = thePreview->height;
p.width = thePreview->width;
prepare_coef (&p);
if (thePreview->bpp == 3) /* [dindinx]: it looks to me that this is always true... */
compute_block_x (buf, rowsize, 0, 0, thePreview->width, thePreview->height,
3, assign_block_3, &p);
else if (thePreview->bpp == 1)
{
compute_block_x (buf, rowsize, 0, 0, thePreview->width,
thePreview->height, 1, assign_block_1, &p);
}
gimp_preview_area_draw (GIMP_PREVIEW_AREA (theWidget),
0, 0, thePreview->width, thePreview->height,
GIMP_RGB_IMAGE,
buf,
rowsize);
g_free (buf);
}
Matrix::Matrix (const ::VSMap &in, ::VSMap &out, void * /*user_data_ptr*/, ::VSCore &core, const ::VSAPI &vsapi)
: vsutl::FilterBase (vsapi, "matrix", ::fmParallel, 0)
, _clip_src_sptr (vsapi.propGetNode (&in, "clip", 0, 0), vsapi)
, _vi_in (*_vsapi.getVideoInfo (_clip_src_sptr.get ()))
, _vi_out (_vi_in)
, _sse_flag (false)
, _sse2_flag (false)
, _avx_flag (false)
, _avx2_flag (false)
, _range_set_src_flag (false)
, _range_set_dst_flag (false)
, _full_range_src_flag (false)
, _full_range_dst_flag (false)
/*, _mat_main ()*/
, _csp_out (fmtcl::ColorSpaceH265_UNSPECIFIED)
, _plane_out (get_arg_int (in, out, "singleout", -1))
, _proc_uptr ()
{
assert (&in != 0);
assert (&out != 0);
assert (&core != 0);
assert (&vsapi != 0);
vsutl::CpuOpt cpu_opt (*this, in, out);
_sse_flag = cpu_opt.has_sse ();
_sse2_flag = cpu_opt.has_sse2 ();
_avx_flag = cpu_opt.has_avx ();
_avx2_flag = cpu_opt.has_avx2 ();
_proc_uptr = std::unique_ptr <fmtcl::MatrixProc> (new fmtcl::MatrixProc (
_sse_flag, _sse2_flag, _avx_flag, _avx2_flag
));
// Checks the input clip
if (_vi_in.format == 0)
{
throw_inval_arg ("only constant pixel formats are supported.");
}
const ::VSFormat & fmt_src = *_vi_in.format;
if (fmt_src.subSamplingW != 0 || fmt_src.subSamplingH != 0)
{
throw_inval_arg ("input must be 4:4:4.");
}
if (fmt_src.numPlanes != NBR_PLANES)
{
throw_inval_arg ("greyscale format not supported as input.");
}
if ( ( fmt_src.sampleType == ::stInteger
&& ( fmt_src.bitsPerSample < 8
|| fmt_src.bitsPerSample > 12)
&& fmt_src.bitsPerSample != 16)
|| ( fmt_src.sampleType == ::stFloat
&& fmt_src.bitsPerSample != 32))
{
throw_inval_arg ("pixel bitdepth not supported.");
}
if (_plane_out >= NBR_PLANES)
{
throw_inval_arg (
"singleout is a plane index and must be -1 or ranging from 0 to 3."
);
}
// Destination colorspace
bool force_col_fam_flag;
const ::VSFormat * fmt_dst_ptr = get_output_colorspace (
in, out, core, fmt_src, _plane_out, force_col_fam_flag
);
if ( fmt_dst_ptr->colorFamily != ::cmGray
&& fmt_dst_ptr->colorFamily != ::cmRGB
&& fmt_dst_ptr->colorFamily != ::cmYUV
&& fmt_dst_ptr->colorFamily != ::cmYCoCg)
{
throw_inval_arg ("unsupported color family for output.");
}
if ( ( fmt_dst_ptr->sampleType == ::stInteger
&& ( fmt_dst_ptr->bitsPerSample < 8
|| fmt_dst_ptr->bitsPerSample > 12)
&& fmt_dst_ptr->bitsPerSample != 16)
|| ( fmt_dst_ptr->sampleType == ::stFloat
&& fmt_dst_ptr->bitsPerSample != 32))
{
throw_inval_arg ("output bitdepth not supported.");
}
if ( fmt_dst_ptr->sampleType != fmt_src.sampleType
|| fmt_dst_ptr->bitsPerSample < fmt_src.bitsPerSample
|| fmt_dst_ptr->subSamplingW != fmt_src.subSamplingW
|| fmt_dst_ptr->subSamplingH != fmt_src.subSamplingH)
{
throw_inval_arg (
"specified output colorspace is not compatible with the input."
);
}
// Preliminary matrix test: deduce the target color family if unspecified
if ( ! force_col_fam_flag
&& fmt_dst_ptr->colorFamily != ::cmGray)
{
int def_count = 0;
def_count += is_arg_defined (in, "mat" ) ? 1 : 0;
def_count += is_arg_defined (in, "mats") ? 1 : 0;
def_count += is_arg_defined (in, "matd") ? 1 : 0;
if (def_count == 1)
{
std::string tmp_mat (get_arg_str (in, out, "mat", ""));
tmp_mat = get_arg_str (in, out, "mats", tmp_mat);
tmp_mat = get_arg_str (in, out, "matd", tmp_mat);
fmtcl::ColorSpaceH265 tmp_csp =
find_cs_from_mat_str (*this, tmp_mat, false);
fmt_dst_ptr = find_dst_col_fam (tmp_csp, fmt_dst_ptr, fmt_src, core);
}
}
// Output format is validated.
_vi_out.format = fmt_dst_ptr;
const ::VSFormat &fmt_dst = *fmt_dst_ptr;
const int nbr_expected_coef = NBR_PLANES * (NBR_PLANES + 1);
bool mat_init_flag = false;
// Matrix presets
std::string mat (get_arg_str (in, out, "mat", ""));
std::string mats (( fmt_src.colorFamily == ::cmYUV ) ? mat : "");
std::string matd (( fmt_dst.colorFamily == ::cmYUV
|| fmt_dst.colorFamily == ::cmGray) ? mat : "");
mats = get_arg_str (in, out, "mats", mats);
matd = get_arg_str (in, out, "matd", matd);
if (! mats.empty () || ! matd.empty ())
{
fstb::conv_to_lower_case (mats);
fstb::conv_to_lower_case (matd);
select_def_mat (mats, fmt_src);
select_def_mat (matd, fmt_dst);
fmtcl::Mat4 m2s;
fmtcl::Mat4 m2d;
make_mat_from_str (m2s, mats, true);
make_mat_from_str (m2d, matd, false);
_csp_out = find_cs_from_mat_str (*this, matd, false);
_mat_main = m2d * m2s;
mat_init_flag = true;
}
// Range
_full_range_src_flag = (get_arg_int (
in, out, "fulls" ,
vsutl::is_full_range_default (fmt_src) ? 1 : 0,
0, &_range_set_src_flag
) != 0);
_full_range_dst_flag = (get_arg_int (
in, out, "fulld",
vsutl::is_full_range_default (fmt_dst) ? 1 : 0,
0, &_range_set_dst_flag
) != 0);
// Custom coefficients
const int nbr_coef = _vsapi.propNumElements (&in, "coef");
const bool custom_mat_flag = (nbr_coef > 0);
if (custom_mat_flag)
{
if (nbr_coef != nbr_expected_coef)
{
throw_inval_arg ("coef has a wrong number of elements.");
}
for (int y = 0; y < NBR_PLANES + 1; ++y)
{
for (int x = 0; x < NBR_PLANES + 1; ++x)
{
_mat_main [y] [x] = (x == y) ? 1 : 0;
if ( (x < fmt_src.numPlanes || x == NBR_PLANES)
&& y < fmt_dst.numPlanes)
{
int err = 0;
const int index = y * (fmt_src.numPlanes + 1) + x;
const double c = _vsapi.propGetFloat (&in, "coef", index, &err);
if (err != 0)
{
throw_rt_err ("error while reading the matrix coefficients.");
}
_mat_main [y] [x] = c;
}
}
}
mat_init_flag = true;
}
if (! mat_init_flag)
{
throw_inval_arg (
"you must specify a matrix preset or a custom coefficient list."
);
}
prepare_coef (fmt_dst, fmt_src);
if (_vsapi.getError (&out) != 0)
{
throw -1;
}
}
