| 273 | } |
| 274 | |
| 275 | void div_image(const imagef &src1, const imagef &src2, imagef &dst, const vec4F &scale) |
| 276 | { |
| 277 | dst.resize(src1); |
| 278 | |
| 279 | //#pragma omp parallel for |
| 280 | for (int y = 0; y < (int)dst.get_height(); y++) |
| 281 | { |
| 282 | for (uint32_t x = 0; x < dst.get_width(); x++) |
| 283 | { |
| 284 | const vec4F &s1 = src1(x, y); |
| 285 | const vec4F &s2 = src2(x, y); |
| 286 | |
| 287 | vec4F d; |
| 288 | |
| 289 | for (uint32_t c = 0; c < 4; c++) |
| 290 | { |
| 291 | float v = s2[c]; |
| 292 | if (v == 0.0f) |
| 293 | d[c] = 0.0f; |
| 294 | else |
| 295 | d[c] = (s1[c] * scale[c]) / v; |
| 296 | } |
| 297 | |
| 298 | dst(x, y) = d; |
| 299 | } |
| 300 | } |
| 301 | } |
| 302 | |
| 303 | vec4F avg_image(const imagef &src) |
| 304 | { |
no test coverage detected