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; } }