void process(struct dt_iop_module_t *self, dt_dev_pixelpipe_iop_t *piece, const void *const ivoid,
void *const ovoid, const dt_iop_roi_t *const roi_in, const dt_iop_roi_t *const roi_out)
{
dt_iop_hotpixels_gui_data_t *g = (dt_iop_hotpixels_gui_data_t *)self->gui_data;
const dt_iop_hotpixels_data_t *data = (dt_iop_hotpixels_data_t *)piece->data;
// The processing loop should output only a few pixels, so just copy everything first
memcpy(ovoid, ivoid, (size_t)roi_out->width * roi_out->height * sizeof(float));
int fixed;
if(piece->pipe->filters == 9u)
{
fixed = process_xtrans(data, ivoid, ovoid, roi_out, (const uint8_t(*const)[6])piece->pipe->xtrans);
}
else
{
fixed = process_bayer(data, ivoid, ovoid, roi_out);
}
if(g != NULL && self->dev->gui_attached && piece->pipe->type == DT_DEV_PIXELPIPE_FULL)
{
g->pixels_fixed = fixed;
}
}
/* Detect hot sensor pixels based on the 4 surrounding sites. Pixels
* having 3 or 4 (depending on permissive setting) surrounding pixels that
* than value*multiplier are considered "hot", and are replaced by the maximum of
* the neighbour pixels. The permissive variant allows for
* correcting pairs of hot pixels in adjacent sites. Replacement using
* the maximum produces fewer artifacts when inadvertently replacing
* non-hot pixels. */
void process(struct dt_iop_module_t *self, dt_dev_pixelpipe_iop_t *piece, void *i, void *o,
const dt_iop_roi_t *roi_in, const dt_iop_roi_t *roi_out)
{
const dt_image_t *img = &self->dev->image_storage;
dt_iop_hotpixels_gui_data_t *g = (dt_iop_hotpixels_gui_data_t *)self->gui_data;
const dt_iop_hotpixels_data_t *data = (dt_iop_hotpixels_data_t *)piece->data;
const float threshold = data->threshold;
const float multiplier = data->multiplier;
const int width = roi_out->width;
const int height = roi_out->height;
const int widthx2 = width * 2;
const gboolean markfixed = data->markfixed;
const int min_neighbours = data->permissive ? 3 : 4;
// The loop should output only a few pixels, so just copy everything first
memcpy(o, i, (size_t)roi_out->width * roi_out->height * sizeof(float));
int fixed = 0;
if(img->filters == 9u)
{
fixed = process_xtrans(i, o, roi_in, width, height, img->xtrans, threshold, multiplier, markfixed,
min_neighbours);
goto processed;
}
// Bayer sensor array
#ifdef _OPENMP
#pragma omp parallel for default(none) shared(roi_out, i, o) reduction(+ : fixed) schedule(static)
#endif
for(int row = 2; row < roi_out->height - 2; row++)
{
const float *in = (float *)i + (size_t)width * row + 2;
float *out = (float *)o + (size_t)width * row + 2;
for(int col = 2; col < width - 1; col++, in++, out++)
{
float mid = *in * multiplier;
if(*in > threshold)
{
int count = 0;
float maxin = 0.0;
float other;
#define TESTONE(OFFSET) \
other = in[OFFSET]; \
if(mid > other) \
{ \
count++; \
if(other > maxin) maxin = other; \
}
TESTONE(-2);
TESTONE(-widthx2);
TESTONE(+2);
TESTONE(+widthx2);
#undef TESTONE
if(count >= min_neighbours)
{
*out = maxin;
fixed++;
if(markfixed)
{
for(int i = -2; i >= -10 && i >= -col; i -= 2) out[i] = *in;
for(int i = 2; i <= 10 && i < width - col; i += 2) out[i] = *in;
}
}
}
}
}
processed:
if(g != NULL && self->dev->gui_attached && piece->pipe->type == DT_DEV_PIXELPIPE_FULL)
{
g->pixels_fixed = fixed;
}
}
