COM_van_vliet_gaussian_coefficients.hh

Go to the documentation of this file.

/* SPDX-FileCopyrightText: 2024 Blender Authors
 *
 * SPDX-License-Identifier: GPL-2.0-or-later */
 
#pragma once
 
#include <cstdint>
#include <memory>
 
#include "BLI_map.hh"
#include "BLI_math_vector_types.hh"
 
#include "COM_cached_resource.hh"
 
namespace blender::compositor {
 
class Context;
 
/* ------------------------------------------------------------------------------------------------
 * Van Vliet Gaussian Coefficients Key.
 */
class VanVlietGaussianCoefficientsKey {
 public:
  float sigma;
 
  VanVlietGaussianCoefficientsKey(float sigma);
 
  uint64_t hash() const;
};
 
bool operator==(const VanVlietGaussianCoefficientsKey &a,
                const VanVlietGaussianCoefficientsKey &b);
 
/* -------------------------------------------------------------------------------------------------
 * Van Vliet Gaussian Coefficients.
 *
 * A caches resource that computes and caches the coefficients of the fourth order IIR filter
 * approximating a Gaussian filter computed using Van Vliet's design method. This is based on the
 * following paper:
 *
 *   Van Vliet, Lucas J., Ian T. Young, and Piet W. Verbeek. "Recursive Gaussian derivative
 *   filters." Proceedings. Fourteenth International Conference on Pattern Recognition (Cat. No.
 *   98EX170). Vol. 1. IEEE, 1998.
 *
 * However, to improve the numerical stability of the filter, it is decomposed into a bank of
 * two parallel second order IIR filters, each having a causal and a non causal filter. */
class VanVlietGaussianCoefficients : public CachedResource {
 private:
  /* The causal and non causal feedforward coefficients for the first second-order filter. */
  double2 first_causal_feedforward_coefficients_;
  double2 first_non_causal_feedforward_coefficients_;
  /* The feedback coefficients for the first second order filter. This is the same for both the
   * causal and non causal filters. */
  double2 first_feedback_coefficients_;
 
  /* The causal and non causal feedforward coefficients for the second second-order filter. */
  double2 second_causal_feedforward_coefficients_;
  double2 second_non_causal_feedforward_coefficients_;
  /* The feedback coefficients for the second second-order filter. This is the same for both the
   * causal and non causal filters. */
  double2 second_feedback_coefficients_;
 
  /* The difference equation of the IIR filter rely on previous outputs to compute the new output,
   * and those previous outputs need to be properly initialized somehow. To do Neumann boundary
   * condition, we multiply the boundary value with this coefficient to simulate an infinite stream
   * of the boundary value. See the implementation for more information. */
  double first_causal_boundary_coefficient_;
  double first_non_causal_boundary_coefficient_;
  double second_causal_boundary_coefficient_;
  double second_non_causal_boundary_coefficient_;
 
 public:
  VanVlietGaussianCoefficients(Context &context, float sigma);
 
  const double2 &first_causal_feedforward_coefficients() const;
  const double2 &first_non_causal_feedforward_coefficients() const;
  const double2 &first_feedback_coefficients() const;
 
  const double2 &second_causal_feedforward_coefficients() const;
  const double2 &second_non_causal_feedforward_coefficients() const;
  const double2 &second_feedback_coefficients() const;
 
  double first_causal_boundary_coefficient() const;
  double first_non_causal_boundary_coefficient() const;
  double second_causal_boundary_coefficient() const;
  double second_non_causal_boundary_coefficient() const;
};
 
/* ------------------------------------------------------------------------------------------------
 * Van Vliet Gaussian Coefficients Container.
 */
class VanVlietGaussianCoefficientsContainer : CachedResourceContainer {
 private:
  Map<VanVlietGaussianCoefficientsKey, std::unique_ptr<VanVlietGaussianCoefficients>> map_;
 
 public:
  void reset() override;
 
  /* Check if there is an available VanVlietGaussianCoefficients cached resource with the given
   * parameters in the container, if one exists, return it, otherwise, return a newly created one
   * and add it to the container. In both cases, tag the cached resource as needed to keep it
   * cached for the next evaluation. */
  VanVlietGaussianCoefficients &get(Context &context, float sigma);
};
 
}  // namespace blender::compositor