Вот простая реализация алгоритма Брезенхема для построения многоточия на вашем image
(адаптировано из https://gist.github.com/bert/1085538):
void set_pixel(image a, int x, int y, float r, float g, float b) {
if (x >= 0 && x < a.w && y >= 0 && y < a.h) {
a.data[x + y * a.w + 0 * a.w * a.h] = r;
a.data[x + y * a.w + 1 * a.w * a.h] = g;
a.data[x + y * a.w + 2 * a.w * a.h] = b;
}
}
void plot_ellipsis_rect(image a, int x0, int y0, int x1, int y1, float r, float g, float b) {
int a = abs(x1 - x0), b = abs(y1 - y0), b1 = b & 1; /* values of diameter */
long dx = 4 * (1 - a) * b * b, dy = 4 * (b1 + 1) * a * a; /* error increment */
long err = dx + dy + b1 * a * a, e2; /* error of 1.step */
y0 += (b + 1) / 2;
y1 = y0 - b1; /* starting pixel */
a = 8 * a * a;
b1 = 8 * b * b;
do {
set_pixel(a, x1, y0, r, g, b); /* I. Quadrant */
set_pixel(a, x0, y0, r, g, b); /* II. Quadrant */
set_pixel(a, x0, y1, r, g, b); /* III. Quadrant */
set_pixel(a, x1, y1, r, g, b); /* IV. Quadrant */
e2 = 2 * err;
if (e2 >= dx) {
x0++;
x1--;
err += dx += b1;
}
if (e2 <= dy) {
y0++;
y1--;
err += dy += a;
}
} while (x0 <= x1);
while (y0 - y1 < b) { /* too early stop of flat ellipses a=1 */
set_pixel(a, x0 - 1, y0, r, g, b); /* -> finish tip of ellipse */
set_pixel(a, x1 + 1, y0++, r, g, b);
set_pixel(a, x0 - 1, y1, r, g, b);
set_pixel(a, x1 + 1, y1--, r, g, b);
}
}