--- /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
+ *
+ */
+
+#include "clip.h"
+#include "bistogramconfig.h"
+#include "units.h"
+
+#include <math.h>
+
+
+
+
+HistogramPoint::HistogramPoint()
+ : ListItem<HistogramPoint>()
+{
+}
+
+HistogramPoint::~HistogramPoint()
+{
+}
+
+int HistogramPoint::equivalent(HistogramPoint *src)
+{
+ return EQUIV(x, src->x) && EQUIV(y, src->y) && EQUIV(gradient, src->gradient);
+}
+
+
+
+
+HistogramPoints::HistogramPoints()
+ : List<HistogramPoint>()
+{
+ insert(0.0,0.0);
+ first->gradient = 1.0;
+ first->xoffset_left = 0.0;
+ first->xoffset_right = 0.05;
+ insert(1.0,1.0);
+ last->gradient = 1.0;
+ last->xoffset_left = -0.05;
+ last->xoffset_right = 0.0;
+
+}
+
+HistogramPoints::~HistogramPoints()
+{
+}
+
+HistogramPoint* HistogramPoints::insert(float x, float y)
+{
+ HistogramPoint *current = first;
+
+// Get existing point after new point
+ while(current)
+ {
+ if(current->x > x)
+ break;
+ else
+ current = NEXT;
+ }
+
+// Insert new point before current point
+ HistogramPoint *new_point = new HistogramPoint;
+ if(current)
+ {
+ insert_before(current, new_point);
+ }
+ else
+// Append new point to list
+ {
+ append(new_point);
+ }
+
+ new_point->x = x;
+ new_point->y = y;
+ new_point->xoffset_left = -0.05;
+ new_point->xoffset_right = 0.05;
+
+
+ return new_point;
+}
+
+void HistogramPoints::boundaries()
+{
+ HistogramPoint *current = first;
+ while(current)
+ {
+ CLAMP(current->x, 0.0, 1.0);
+ CLAMP(current->y, 0.0, 1.0);
+ current = NEXT;
+ }
+}
+
+int HistogramPoints::equivalent(HistogramPoints *src)
+{
+ HistogramPoint *current_this = first;
+ HistogramPoint *current_src = src->first;
+ while(current_this && current_src)
+ {
+ if(!current_this->equivalent(current_src)) return 0;
+ current_this = current_this->next;
+ current_src = current_src->next;
+ }
+
+ return !current_this ^ !current_src ? 0 : 1;
+}
+
+void HistogramPoints::copy_from(HistogramPoints *src)
+{
+ while(last)
+ delete last;
+ HistogramPoint *current = src->first;
+ while(current)
+ {
+ HistogramPoint *new_point = new HistogramPoint;
+ new_point->x = current->x;
+ new_point->y = current->y;
+ append(new_point);
+ current = NEXT;
+ }
+}
+
+void HistogramPoints::interpolate(HistogramPoints *prev,
+ HistogramPoints *next,
+ double prev_scale,
+ double next_scale)
+{
+ HistogramPoint *current = first;
+ HistogramPoint *current_prev = prev->first;
+ HistogramPoint *current_next = next->first;
+
+ while(current && current_prev && current_next)
+ {
+ current->x = current_prev->x * prev_scale +
+ current_next->x * next_scale;
+ current->y = current_prev->y * prev_scale +
+ current_next->y * next_scale;
+ current = NEXT;
+ current_prev = current_prev->next;
+ current_next = current_next->next;
+ }
+}
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+HistogramConfig::HistogramConfig()
+{
+ reset(1);
+}
+
+void HistogramConfig::reset(int do_mode)
+{
+ reset_points();
+
+
+ for(int i = 0; i < HISTOGRAM_MODES; i++)
+ {
+ output_min[i] = 0.0;
+ output_max[i] = 1.0;
+ }
+
+ if(do_mode)
+ {
+ automatic = 0;
+ threshold = 0.1;
+ split = 0;
+ smoothMode = 0;
+ }
+}
+
+void HistogramConfig::reset_points()
+{
+ for(int i = 0; i < HISTOGRAM_MODES; i++)
+ {
+ while(points[i].last) delete points[i].last;
+ points[i].insert(0.0,0.0);
+ points[i].last->gradient = 1.0;
+ points[i].insert(1.0,1.0);
+ points[i].last->gradient = 1.0;
+ }
+}
+
+
+void HistogramConfig::boundaries()
+{
+ for(int i = 0; i < HISTOGRAM_MODES; i++)
+ {
+ points[i].boundaries();
+ CLAMP(output_min[i], HIST_MIN_INPUT, HIST_MAX_INPUT);
+ CLAMP(output_max[i], HIST_MIN_INPUT, HIST_MAX_INPUT);
+ output_min[i] = Units::quantize(output_min[i], PRECISION);
+ output_max[i] = Units::quantize(output_max[i], PRECISION);
+ }
+ CLAMP(threshold, 0, 1);
+}
+
+int HistogramConfig::equivalent(HistogramConfig &that)
+{
+ for(int i = 0; i < HISTOGRAM_MODES; i++)
+ {
+ if(!points[i].equivalent(&that.points[i]) ||
+ !EQUIV(output_min[i], that.output_min[i]) ||
+ !EQUIV(output_max[i], that.output_max[i])) return 0;
+ }
+
+ if(automatic != that.automatic ||
+ !EQUIV(threshold, that.threshold)) return 0;
+
+ return 1;
+}
+
+void HistogramConfig::copy_from(HistogramConfig &that)
+{
+ for(int i = 0; i < HISTOGRAM_MODES; i++)
+ {
+ points[i].copy_from(&that.points[i]);
+ output_min[i] = that.output_min[i];
+ output_max[i] = that.output_max[i];
+ }
+
+ automatic = that.automatic;
+ threshold = that.threshold;
+ split = that.split;
+ smoothMode = that.smoothMode;
+}
+
+void HistogramConfig::interpolate(HistogramConfig &prev,
+ HistogramConfig &next,
+ int64_t prev_frame,
+ int64_t next_frame,
+ int64_t current_frame)
+{
+ double next_scale = (double)(current_frame - prev_frame) / (next_frame - prev_frame);
+ double prev_scale = 1.0 - next_scale;
+
+ for(int i = 0; i < HISTOGRAM_MODES; i++)
+ {
+ points[i].interpolate(&prev.points[i], &next.points[i], next_scale, prev_scale);
+ output_min[i] = prev.output_min[i] * prev_scale + next.output_min[i] * next_scale;
+ output_max[i] = prev.output_max[i] * prev_scale + next.output_max[i] * next_scale;
+ }
+
+ threshold = prev.threshold * prev_scale + next.threshold * next_scale;
+ automatic = prev.automatic;
+ split = prev.split;
+ smoothMode = prev.smoothMode;
+}
+
+
+void HistogramConfig::dump()
+{
+ for(int j = 0; j < HISTOGRAM_MODES; j++)
+ {
+ printf("HistogramConfig::dump mode=%d\n", j);
+ HistogramPoints *points = &this->points[j];
+ HistogramPoint *current = points->first;
+ while(current)
+ {
+ printf("%f,%f ", current->x, current->y);
+ fflush(stdout);
+ current = NEXT;
+ }
+ printf("\n");
+ }
+}
+
+
+