Remove off-by-.5 errors and previous attempts to fix them by moving the rotation center.
The bounding box is wrong (too large) but it does not seem worthwhile to fix this.

@@ -71,3 +70,1 @@
-  area.xmax = ceilf(scale_coord(area.xmax, center_x, relative_scale_x));
-  area.ymin = floorf(scale_coord(area.ymin, center_y, relative_scale_y));
-  area.ymax = ceilf(scale_coord(area.ymax, center_y, relative_scale_y));
+  area.xmin = floorf(scale_coord(area.xmin, center_x, relative_scale_x) + 0.1f);

@@ -542,0 +435,4 @@
+    const float w = get_width();
+    const float h = get_height();
+
+    if (sampler == PixelSampler::Nearest) {

@@ -546,2 +447,2 @@
-      this->wrap_pixel(u, v, extend_x, extend_y);
-      this->read_elem_sampled(u, v, sampler, result);
+      // compute (linear interpolation) intersection with Clip
+      float m = 1.0f;

@@ -548,0 +456,4 @@
+        clear_elem(result);
+        return;
+      }
+      if (is_a_single_elem_) {

@@ -548,0 +461,4 @@
+      }
+      else if (sampler == PixelSampler::Bilinear) {
+        // Sample using Extend or Repeat
+        math::interpolate_bilinear_wrap_fl(buffer_, result, w, h, num_channels_, x, y,

@@ -548,0 +467,4 @@
+      }
+      else { // PixelSampler::Bicubic
+        // Wrap is not implemented in interpolate_cubic_bspline, this uses Extend always.
+        // This will produce minor errors in pixels next to wrapped borders.

@@ -548,0 +454,4 @@
+    }
+
+    if (is_a_single_elem_) {
+      memcpy(result, buffer_, get_elem_bytes_len());