--- /dev/null
+
+/*
+ * CINELERRA
+ * Copyright (C) 2008 Adam Williams <broadcast at earthling dot net>
+ *
+ * This program is free software; you can redistribute it and/or modify
+ * it under the terms of the GNU General Public License as published by
+ * the Free Software Foundation; either version 2 of the License, or
+ * (at your option) any later version.
+ *
+ * This program is distributed in the hope that it will be useful,
+ * but WITHOUT ANY WARRANTY; without even the implied warranty of
+ * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+ * GNU General Public License for more details.
+ *
+ * You should have received a copy of the GNU General Public License
+ * along with this program; if not, write to the Free Software
+ * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
+ *
+ */
+
+#ifndef DENOISE_H
+#define DENOISE_H
+
+
+
+#include "bchash.inc"
+#include "guicast.h"
+#include "mutex.h"
+#include "pluginaclient.h"
+#include "vframe.inc"
+
+class DenoiseEffect;
+typedef enum { DECOMP, RECON } wavetype;
+
+class DenoiseLevel : public BC_FPot
+{
+public:
+ DenoiseLevel(DenoiseEffect *plugin, int x, int y);
+ int handle_event();
+ DenoiseEffect *plugin;
+};
+
+class DenoiseWindow : public PluginClientWindow
+{
+public:
+ DenoiseWindow(DenoiseEffect *plugin);
+ void create_objects();
+ void update();
+ DenoiseLevel *scale;
+ DenoiseEffect *plugin;
+};
+
+
+class DenoiseConfig
+{
+public:
+ DenoiseConfig();
+ void copy_from(DenoiseConfig &that);
+ int equivalent(DenoiseConfig &that);
+ void interpolate(DenoiseConfig &prev,
+ DenoiseConfig &next,
+ int64_t prev_frame,
+ int64_t next_frame,
+ int64_t current_frame);
+ double level;
+};
+
+class Tree
+{
+public:
+ Tree(int input_length, int levels);
+ ~Tree();
+
+ int input_length;
+ int levels;
+ double **values;
+};
+
+class WaveletCoeffs
+{
+public:
+ WaveletCoeffs(double alpha, double beta);
+ ~WaveletCoeffs();
+
+ double values[6];
+ int length;
+};
+
+class WaveletFilters
+{
+public:
+ WaveletFilters(WaveletCoeffs *wave_coeffs, wavetype transform);
+ ~WaveletFilters();
+
+ double g[6], h[6];
+ int length;
+};
+
+class DenoiseEffect : public PluginAClient
+{
+public:
+ DenoiseEffect(PluginServer *server);
+ ~DenoiseEffect();
+
+ int is_realtime();
+ void read_data(KeyFrame *keyframe);
+ void save_data(KeyFrame *keyframe);
+ int process_realtime(int64_t size, Samples *input_ptr, Samples *output_ptr);
+
+
+
+
+ void reset();
+ void update_gui();
+ void delete_dsp();
+
+ void process_window();
+ double dot_product(double *data, double *filter, char filtlen);
+ int convolve_dec_2(double *input_sequence,
+ int64_t length,
+ double *filter,
+ int filtlen,
+ double *output_sequence);
+ int64_t decompose_branches(double *in_data,
+ int64_t length,
+ WaveletFilters *decomp_filter,
+ double *out_low,
+ double *out_high);
+ int wavelet_decomposition(double *in_data,
+ int64_t in_length,
+ double **out_data);
+ int tree_copy(double **output,
+ double **input,
+ int length,
+ int levels);
+ int threshold(int window_size, double gammas, int levels);
+ double dot_product_even(double *data, double *filter, int filtlen);
+ double dot_product_odd(double *data, double *filter, int filtlen);
+ int convolve_int_2(double *input_sequence,
+ int64_t length,
+ double *filter,
+ int filtlen,
+ int sum_output,
+ double *output_sequence);
+ int64_t reconstruct_branches(double *in_low,
+ double *in_high,
+ int64_t in_length,
+ WaveletFilters *recon_filter,
+ double *output);
+ int wavelet_reconstruction(double **in_data,
+ int64_t in_length,
+ double *out_data);
+
+
+ PLUGIN_CLASS_MEMBERS(DenoiseConfig)
+
+// buffer for storing fragments until a complete window size is armed
+ double *input_buffer;
+ int64_t input_size;
+ int64_t input_allocation;
+// buffer for storing fragments until a fragment is ready to be read
+ double *output_buffer;
+ int64_t output_size;
+ int64_t output_allocation;
+ double *dsp_in;
+ double *dsp_out;
+// buffer for capturing output of a single iteration
+ double *dsp_iteration;
+ Tree *ex_coeff_d, *ex_coeff_r, *ex_coeff_rn;
+ WaveletCoeffs *wave_coeff_d, *wave_coeff_r;
+ WaveletFilters *decomp_filter, *recon_filter;
+// scaling factor for transferring from input_buffer
+ double in_scale;
+// power converted to scaling factor
+ double out_scale;
+
+// depends on the type of music
+ int64_t levels;
+// higher number reduces aliasing due to a high noise_level
+// also increases high end
+ int64_t iterations;
+// daub6 coeffs
+ double alpha;
+ double beta;
+// power
+ float output_level;
+// higher number kills more noise at the expense of more aliasing
+ float noise_level;
+ int64_t window_size;
+ int first_window;
+ int initialized;
+};
+
+
+#endif