Blame src/utils/rescaler_utils.c

Packit 9c6abc
// Copyright 2012 Google Inc. All Rights Reserved.
Packit 9c6abc
//
Packit 9c6abc
// Use of this source code is governed by a BSD-style license
Packit 9c6abc
// that can be found in the COPYING file in the root of the source
Packit 9c6abc
// tree. An additional intellectual property rights grant can be found
Packit 9c6abc
// in the file PATENTS. All contributing project authors may
Packit 9c6abc
// be found in the AUTHORS file in the root of the source tree.
Packit 9c6abc
// -----------------------------------------------------------------------------
Packit 9c6abc
//
Packit 9c6abc
// Rescaling functions
Packit 9c6abc
//
Packit 9c6abc
// Author: Skal (pascal.massimino@gmail.com)
Packit 9c6abc
Packit 9c6abc
#include <assert.h>
Packit 9c6abc
#include <stdlib.h>
Packit 9c6abc
#include <string.h>
Packit 9c6abc
#include "src/dsp/dsp.h"
Packit 9c6abc
#include "src/utils/rescaler_utils.h"
Packit 9c6abc
Packit 9c6abc
//------------------------------------------------------------------------------
Packit 9c6abc
Packit 9c6abc
void WebPRescalerInit(WebPRescaler* const wrk, int src_width, int src_height,
Packit 9c6abc
                      uint8_t* const dst,
Packit 9c6abc
                      int dst_width, int dst_height, int dst_stride,
Packit 9c6abc
                      int num_channels, rescaler_t* const work) {
Packit 9c6abc
  const int x_add = src_width, x_sub = dst_width;
Packit 9c6abc
  const int y_add = src_height, y_sub = dst_height;
Packit 9c6abc
  wrk->x_expand = (src_width < dst_width);
Packit 9c6abc
  wrk->y_expand = (src_height < dst_height);
Packit 9c6abc
  wrk->src_width = src_width;
Packit 9c6abc
  wrk->src_height = src_height;
Packit 9c6abc
  wrk->dst_width = dst_width;
Packit 9c6abc
  wrk->dst_height = dst_height;
Packit 9c6abc
  wrk->src_y = 0;
Packit 9c6abc
  wrk->dst_y = 0;
Packit 9c6abc
  wrk->dst = dst;
Packit 9c6abc
  wrk->dst_stride = dst_stride;
Packit 9c6abc
  wrk->num_channels = num_channels;
Packit 9c6abc
Packit 9c6abc
  // for 'x_expand', we use bilinear interpolation
Packit 9c6abc
  wrk->x_add = wrk->x_expand ? (x_sub - 1) : x_add;
Packit 9c6abc
  wrk->x_sub = wrk->x_expand ? (x_add - 1) : x_sub;
Packit 9c6abc
  if (!wrk->x_expand) {  // fx_scale is not used otherwise
Packit 9c6abc
    wrk->fx_scale = WEBP_RESCALER_FRAC(1, wrk->x_sub);
Packit 9c6abc
  }
Packit 9c6abc
  // vertical scaling parameters
Packit 9c6abc
  wrk->y_add = wrk->y_expand ? y_add - 1 : y_add;
Packit 9c6abc
  wrk->y_sub = wrk->y_expand ? y_sub - 1 : y_sub;
Packit 9c6abc
  wrk->y_accum = wrk->y_expand ? wrk->y_sub : wrk->y_add;
Packit 9c6abc
  if (!wrk->y_expand) {
Packit 9c6abc
    // This is WEBP_RESCALER_FRAC(dst_height, x_add * y_add) without the cast.
Packit 9c6abc
    // Its value is <= WEBP_RESCALER_ONE, because dst_height <= wrk->y_add, and
Packit 9c6abc
    // wrk->x_add >= 1;
Packit 9c6abc
    const uint64_t ratio =
Packit 9c6abc
        (uint64_t)dst_height * WEBP_RESCALER_ONE / (wrk->x_add * wrk->y_add);
Packit 9c6abc
    if (ratio != (uint32_t)ratio) {
Packit 9c6abc
      // When ratio == WEBP_RESCALER_ONE, we can't represent the ratio with the
Packit 9c6abc
      // current fixed-point precision. This happens when src_height ==
Packit 9c6abc
      // wrk->y_add (which == src_height), and wrk->x_add == 1.
Packit 9c6abc
      // => We special-case fxy_scale = 0, in WebPRescalerExportRow().
Packit 9c6abc
      wrk->fxy_scale = 0;
Packit 9c6abc
    } else {
Packit 9c6abc
      wrk->fxy_scale = (uint32_t)ratio;
Packit 9c6abc
    }
Packit 9c6abc
    wrk->fy_scale = WEBP_RESCALER_FRAC(1, wrk->y_sub);
Packit 9c6abc
  } else {
Packit 9c6abc
    wrk->fy_scale = WEBP_RESCALER_FRAC(1, wrk->x_add);
Packit 9c6abc
    // wrk->fxy_scale is unused here.
Packit 9c6abc
  }
Packit 9c6abc
  wrk->irow = work;
Packit 9c6abc
  wrk->frow = work + num_channels * dst_width;
Packit 9c6abc
  memset(work, 0, 2 * dst_width * num_channels * sizeof(*work));
Packit 9c6abc
Packit 9c6abc
  WebPRescalerDspInit();
Packit 9c6abc
}
Packit 9c6abc
Packit 9c6abc
int WebPRescalerGetScaledDimensions(int src_width, int src_height,
Packit 9c6abc
                                    int* const scaled_width,
Packit 9c6abc
                                    int* const scaled_height) {
Packit 9c6abc
  assert(scaled_width != NULL);
Packit 9c6abc
  assert(scaled_height != NULL);
Packit 9c6abc
  {
Packit 9c6abc
    int width = *scaled_width;
Packit 9c6abc
    int height = *scaled_height;
Packit 9c6abc
Packit 9c6abc
    // if width is unspecified, scale original proportionally to height ratio.
Packit 9c6abc
    if (width == 0) {
Packit 9c6abc
      width =
Packit 9c6abc
          (int)(((uint64_t)src_width * height + src_height / 2) / src_height);
Packit 9c6abc
    }
Packit 9c6abc
    // if height is unspecified, scale original proportionally to width ratio.
Packit 9c6abc
    if (height == 0) {
Packit 9c6abc
      height =
Packit 9c6abc
          (int)(((uint64_t)src_height * width + src_width / 2) / src_width);
Packit 9c6abc
    }
Packit 9c6abc
    // Check if the overall dimensions still make sense.
Packit 9c6abc
    if (width <= 0 || height <= 0) {
Packit 9c6abc
      return 0;
Packit 9c6abc
    }
Packit 9c6abc
Packit 9c6abc
    *scaled_width = width;
Packit 9c6abc
    *scaled_height = height;
Packit 9c6abc
    return 1;
Packit 9c6abc
  }
Packit 9c6abc
}
Packit 9c6abc
Packit 9c6abc
//------------------------------------------------------------------------------
Packit 9c6abc
// all-in-one calls
Packit 9c6abc
Packit 9c6abc
int WebPRescaleNeededLines(const WebPRescaler* const wrk, int max_num_lines) {
Packit 9c6abc
  const int num_lines = (wrk->y_accum + wrk->y_sub - 1) / wrk->y_sub;
Packit 9c6abc
  return (num_lines > max_num_lines) ? max_num_lines : num_lines;
Packit 9c6abc
}
Packit 9c6abc
Packit 9c6abc
int WebPRescalerImport(WebPRescaler* const wrk, int num_lines,
Packit 9c6abc
                       const uint8_t* src, int src_stride) {
Packit 9c6abc
  int total_imported = 0;
Packit 9c6abc
  while (total_imported < num_lines && !WebPRescalerHasPendingOutput(wrk)) {
Packit 9c6abc
    if (wrk->y_expand) {
Packit 9c6abc
      rescaler_t* const tmp = wrk->irow;
Packit 9c6abc
      wrk->irow = wrk->frow;
Packit 9c6abc
      wrk->frow = tmp;
Packit 9c6abc
    }
Packit 9c6abc
    WebPRescalerImportRow(wrk, src);
Packit 9c6abc
    if (!wrk->y_expand) {     // Accumulate the contribution of the new row.
Packit 9c6abc
      int x;
Packit 9c6abc
      for (x = 0; x < wrk->num_channels * wrk->dst_width; ++x) {
Packit 9c6abc
        wrk->irow[x] += wrk->frow[x];
Packit 9c6abc
      }
Packit 9c6abc
    }
Packit 9c6abc
    ++wrk->src_y;
Packit 9c6abc
    src += src_stride;
Packit 9c6abc
    ++total_imported;
Packit 9c6abc
    wrk->y_accum -= wrk->y_sub;
Packit 9c6abc
  }
Packit 9c6abc
  return total_imported;
Packit 9c6abc
}
Packit 9c6abc
Packit 9c6abc
int WebPRescalerExport(WebPRescaler* const rescaler) {
Packit 9c6abc
  int total_exported = 0;
Packit 9c6abc
  while (WebPRescalerHasPendingOutput(rescaler)) {
Packit 9c6abc
    WebPRescalerExportRow(rescaler);
Packit 9c6abc
    ++total_exported;
Packit 9c6abc
  }
Packit 9c6abc
  return total_exported;
Packit 9c6abc
}
Packit 9c6abc
Packit 9c6abc
//------------------------------------------------------------------------------