Changeset 2811 for libpipi/trunk/pipi/dither/ordered.c

Timestamp:

09/01/08 01:05:50 (5 years ago)

Author:

sam

Message:

• ordered.c: allow to rotate the dither pattern, using nearest-neighbour rotation interpolation.

File:

• libpipi/trunk/pipi/dither/ordered.c (modified) (4 diffs)

libpipi/trunk/pipi/dither/ordered.c

@@ -21 +21 @@
 
 #include <stdlib.h>
+#include <math.h>
 
 #include "pipi.h"
@@ -27 +28 @@
 pipi_image_t *pipi_dither_ordered(pipi_image_t *img, pipi_image_t *kernel)
 {
+    return pipi_dither_ordered_ext(img, kernel, 0.0);
+}
+
+pipi_image_t *pipi_dither_ordered_ext(pipi_image_t *img, pipi_image_t *kernel,
+                                      double angle)
+{
+    double sint, cost;
     pipi_image_t *dst;
     pipi_pixels_t *dstp, *kernelp;
     float *dstdata, *kerneldata;
-    int x, y, w, h, kw, kh;
+    int x, y, w, h, kx, ky, kw, kh;
 
     w = img->w;
@@ -36 +44 @@
     kw = kernel->w;
     kh = kernel->h;
+
+    cost = cos(angle * (M_PI / 180));
+    sint = sin(angle * (M_PI / 180));
 
     dst = pipi_copy(img);
@@ -50 +61 @@
             float p, q;
 
+            kx = (int)(cost * x - sint * y + 2 * w * h) % kw;
+            ky = (int)(cost * y + sint * x + 2 * w * h) % kh;
+
             p = dstdata[y * w + x];
-            q = p > kerneldata[(y % kh) * kw + (x % kw)] ? 1. : 0.;
+            q = p > kerneldata[ky * kw + kx] ? 1. : 0.;
             dstdata[y * w + x] = q;
         }