1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
|
/* This file is an image processing operation for GEGL
*
* GEGL is free software; you can redistribute it and/or
* modify it under the terms of the GNU Lesser General Public
* License as published by the Free Software Foundation; either
* version 3 of the License, or (at your option) any later version.
*
* GEGL is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
* Lesser General Public License for more details.
*
* You should have received a copy of the GNU Lesser General Public
* License along with GEGL; if not, see <http://www.gnu.org/licenses/>.
*
* Copyright 2013 Victor Oliveira <victormatheus@gmail.com>
* Copyright 2013 Carlos Zubieta <czubieta.dev@gmail.com>
*/
int
block_index (int pos,
int size)
{
return pos < 0 ? ((pos + 1) / size - 1) : (pos / size);
}
__kernel void calc_block_color(__global float4 *in,
__global float4 *out,
int xsize,
int ysize,
int roi_x,
int roi_y,
int4 bbox,
int line_width,
int block_count_x )
{
int gidx = get_global_id(0);
int gidy = get_global_id(1);
int cx = block_index (roi_x, xsize) + gidx;
int cy = block_index (roi_y, ysize) + gidy;
int px0 = max (bbox.s0, cx * xsize) - roi_x + xsize;
int py0 = max (bbox.s1, cy * ysize) - roi_y + ysize;
int px1 = min (bbox.s2, cx * xsize + xsize) - roi_x + xsize;
int py1 = min (bbox.s3, cy * ysize + ysize) - roi_y + ysize;
int i, j;
float4 col = (float4)(0.0f, 0.0f, 0.0f, 0.0f);
int real_xsize = px1 - px0;
int real_ysize = py1 - py0;
float weight = 1.0f / (real_xsize * real_ysize);
for (j = py0; j < py1; ++j)
{
for (i = px0; i < px1; ++i)
{
col += in[j * line_width + i];
}
}
out[gidy * block_count_x + gidx] = col * weight;
}
#define NORM_MANHATTAN 0
#define NORM_EUCLIDEAN 1
#define NORM_INFINITY 2
#define SQR(x) ((x)*(x))
__kernel void kernel_pixelize(__global float4 *in,
__global float4 *out,
int xsize,
int ysize,
float xratio,
float yratio,
int roi_x,
int roi_y,
float4 bg_color,
int norm,
int block_count_x)
{
int gidx = get_global_id(0);
int gidy = get_global_id(1);
int src_width = get_global_size(0);
int cx = block_index (gidx + roi_x, xsize) - block_index (roi_x, xsize);
int cy = block_index (gidy + roi_y, ysize) - block_index (roi_y, ysize);
float4 grid_color = in[cx + cy * block_count_x];
float4 out_color = bg_color;
int x_pos = gidx + roi_x;
int y_pos = gidy + roi_y;
int rect_shape_width = ceil (xsize * xratio);
int rect_shape_height = ceil (ysize * yratio);
int off_shape_x = floor ((xsize - xratio * xsize) / 2.0f);
int off_shape_y = floor ((ysize - yratio * ysize) / 2.0f);
int start_x = block_index (x_pos, xsize) * xsize - roi_x;
int start_y = block_index (y_pos, ysize) * ysize - roi_y;
float shape_area = rect_shape_width * rect_shape_height;
float center_x = start_x + off_shape_x + (float)(rect_shape_width) / 2.0f;
float center_y = start_y + off_shape_y + (float)(rect_shape_height) / 2.0f;
if (norm == NORM_MANHATTAN &&
(fabs (gidx - center_x) * rect_shape_height +
fabs (gidy - center_y) * rect_shape_width
< shape_area))
out_color = grid_color;
if (norm == NORM_EUCLIDEAN &&
SQR ((gidx - center_x) / (float)rect_shape_width) +
SQR ((gidy - center_y) / (float)rect_shape_height) <= 1.0f)
out_color = grid_color;
if (norm == NORM_INFINITY &&
(gidx >= start_x + off_shape_x &&
gidy >= start_y + off_shape_y &&
gidx < start_x + off_shape_x + rect_shape_width &&
gidy < start_y + off_shape_y + rect_shape_height))
out_color = grid_color;
out[gidx + gidy * src_width] = out_color;
}
|
