annotate src/imgproc.c @ 853:2b6fa2618f76

[xemacs-hg @ 2002-05-28 08:44:22 by ben] merge my stderr-proc ws make-docfile.c: Fix places where we forget to check for EOF. code-init.el: Don't use CRLF conversion by default on process output. CMD.EXE and friends work both ways but Cygwin programs don't like the CRs. code-process.el, multicast.el, process.el: Removed. Improvements to call-process-internal: -- allows a buffer to be specified for input and stderr output -- use it on all systems -- implement C-g as documented -- clean up and comment call-process-region uses new call-process facilities; no temp file. remove duplicate funs in process.el. comment exactly how coding systems work and fix various problems. open-multicast-group now does similar coding-system frobbing to open-network-stream. dumped-lisp.el, faces.el, msw-faces.el: Fix some hidden errors due to code not being defined at the right time. xemacs.mak: Add -DSTRICT. ================================================================ ALLOW SEPARATION OF STDOUT AND STDERR IN PROCESSES ================================================================ Standard output and standard error can be processed separately in a process. Each can have its own buffer, its own mark in that buffer, and its filter function. You can specify a separate buffer for stderr in `start-process' to get things started, or use the new primitives: set-process-stderr-buffer process-stderr-buffer process-stderr-mark set-process-stderr-filter process-stderr-filter Also, process-send-region takes a 4th optional arg, a buffer. Currently always uses a pipe() under Unix to read the error output. (#### Would a PTY be better?) sysdep.h, sysproc.h, unexfreebsd.c, unexsunos4.c, nt.c, emacs.c, callproc.c, symsinit.h, sysdep.c, Makefile.in.in, process-unix.c: Delete callproc.c. Move child_setup() to process-unix.c. wait_for_termination() now only needed on a few really old systems. console-msw.h, event-Xt.c, event-msw.c, event-stream.c, event-tty.c, event-unixoid.c, events.h, process-nt.c, process-unix.c, process.c, process.h, procimpl.h: Rewrite the process methods to handle a separate channel for error input. Create Lstreams for reading in the error channel. Many process methods need change. In general the changes are fairly clear as they involve duplicating what's used for reading the normal stdout and changing for stderr -- although tedious, as such changes are required throughout the entire process code. Rewrote the code that reads process output to do two loops, one for stdout and one for stderr. gpmevent.c, tooltalk.c: set_process_filter takes an argument for stderr. ================================================================ NEW ERROR-TRAPPING MECHANISM ================================================================ Totally rewrite error trapping code to be unified and support more features. Basic function is call_trapping_problems(), which lets you specify, by means of flags, what sorts of problems you want trapped. these can include -- quit -- errors -- throws past the function -- creation of "display objects" (e.g. buffers) -- deletion of already-existing "display objects" (e.g. buffers) -- modification of already-existing buffers -- entering the debugger -- gc -- errors->warnings (ala suspended errors) etc. All other error funs rewritten in terms of this one. Various older mechanisms removed or rewritten. window.c, insdel.c, console.c, buffer.c, device.c, frame.c: When creating a display object, added call to note_object_created(), for use with trapping_problems mechanism. When deleting, call check_allowed_operation() and note_object deleted(). The trapping-problems code records the objects created since the call-trapping-problems began. Those objects can be deleted, but none others (i.e. previously existing ones). bytecode.c, cmdloop.c: internal_catch takes another arg. eval.c: Add long comments describing the "five lists" used to maintain state (backtrace, gcpro, specbind, etc.) in the Lisp engine. backtrace.h, eval.c: Implement trapping-problems mechanism, eliminate old mechanisms or redo in terms of new one. frame.c, gutter.c: Flush out the concept of "critical display section", defined by the in_display() var. Use an internal_bind() to get it reset, rather than just doing it at end, because there may be a non-local exit. event-msw.c, event-stream.c, console-msw.h, device.c, dialog-msw.c, frame.c, frame.h, intl.c, toolbar.c, menubar-msw.c, redisplay.c, alloc.c, menubar-x.c: Make use of new trapping-errors stuff and rewrite code based on old mechanisms. glyphs-widget.c, redisplay.h: Protect calling Lisp in redisplay. insdel.c: Protect hooks against deleting existing buffers. frame-msw.c: Use EQ, not EQUAL in hash tables whose keys are just numbers. Otherwise we run into stickiness in redisplay because internal_equal() can QUIT. ================================================================ SIGNAL, C-G CHANGES ================================================================ Here we change the way that C-g interacts with event reading. The idea is that a C-g occurring while we're reading a user event should be read as C-g, but elsewhere should be a QUIT. The former code did all sorts of bizarreness -- requiring that no QUIT occurs anywhere in event-reading code (impossible to enforce given the stuff called or Lisp code invoked), and having some weird system involving enqueue/dequeue of a C-g and interaction with Vquit_flag -- and it didn't work. Now, we simply enclose all code where we want C-g read as an event with {begin/end}_dont_check_for_quit(). This completely turns off the mechanism that checks (and may remove or alter) C-g in the read-ahead queues, so we just get the C-g normal. Signal.c documents this very carefully. cmdloop.c: Correct use of dont_check_for_quit to new scheme, remove old out-of-date comments. event-stream.c: Fix C-g handling to actually work. device-x.c: Disable quit checking when err out. signal.c: Cleanup. Add large descriptive comment. process-unix.c, process-nt.c, sysdep.c: Use QUIT instead of REALLY_QUIT. It's not necessary to use REALLY_QUIT and just confuses the issue. lisp.h: Comment quit handlers. ================================================================ CONS CHANGES ================================================================ free_cons() now takes a Lisp_Object not the result of XCONS(). car and cdr have been renamed so that they don't get used directly; go through XCAR(), XCDR() instead. alloc.c, dired.c, editfns.c, emodules.c, fns.c, glyphs-msw.c, glyphs-x.c, glyphs.c, keymap.c, minibuf.c, search.c, eval.c, lread.c, lisp.h: Correct free_cons calling convention: now takes Lisp_Object, not Lisp_Cons chartab.c: Eliminate direct use of ->car, ->cdr, should be black box. callint.c: Rewrote using EXTERNAL_LIST_LOOP to avoid use of Lisp_Cons. ================================================================ USE INTERNAL-BIND-* ================================================================ eval.c: Cleanups of these funs. alloc.c, fileio.c, undo.c, specifier.c, text.c, profile.c, lread.c, redisplay.c, menubar-x.c, macros.c: Rewrote to use internal_bind_int() and internal_bind_lisp_object() in place of whatever varied and cumbersome mechanisms were formerly there. ================================================================ SPECBIND SANITY ================================================================ backtrace.h: - Improved comments backtrace.h, bytecode.c, eval.c: Add new mechanism check_specbind_stack_sanity() for sanity checking code each time the catchlist or specbind stack change. Removed older prototype of same mechanism. ================================================================ MISC ================================================================ lisp.h, insdel.c, window.c, device.c, console.c, buffer.c: Fleshed out authorship. device-msw.c: Correct bad Unicode-ization. print.c: Be more careful when not initialized or in fatal error handling. search.c: Eliminate running_asynch_code, an FSF holdover. alloc.c: Added comments about gc-cons-threshold. dialog-x.c: Use begin_gc_forbidden() around code to build up a widget value tree, like in menubar-x.c. gui.c: Use Qunbound not Qnil as the default for gethash. lisp-disunion.h, lisp-union.h: Added warnings on use of VOID_TO_LISP(). lisp.h: Use ERROR_CHECK_STRUCTURES to turn on ERROR_CHECK_TRAPPING_PROBLEMS and ERROR_CHECK_TYPECHECK lisp.h: Add assert_with_message. lisp.h: Add macros for gcproing entire arrays. (You could do this before but it required manual twiddling the gcpro structure.) lisp.h: Add prototypes for new functions defined elsewhere.
author ben
date Tue, 28 May 2002 08:45:36 +0000
parents b39c14581166
children a8d8f419b459
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
428
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
1 /* Image processing functions
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
2 Copyright (C) 1998 Jareth Hein
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
3
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
4 This file is a part of XEmacs
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
5
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
6 XEmacs is free software; you can redistribute it and/or modify it
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
7 under the terms of the GNU General Public License as published by the
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
8 Free Software Foundation; either version 2, or (at your option) any
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
9 later version.
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
10
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
11 XEmacs is distributed in the hope that it will be useful, but WITHOUT
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
12 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
13 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
14 for more details.
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
15
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
16 You should have received a copy of the GNU General Public License
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
17 along with XEmacs; see the file COPYING. If not, write to
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
18 the Free Software Foundation, Inc., 59 Temple Place - Suite 330,
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
19 Boston, MA 02111-1307, USA. */
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
20
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
21 /* Synched up with: Not in FSF. */
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
22
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
23 /* Original author: Jareth Hein */
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
24
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
25 /* Parts of this file are based on code from Sam Leffler's tiff library,
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
26 with the original copyright displayed here:
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
27
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
28 Copyright (c) 1988-1997 Sam Leffler
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
29 Copyright (c) 1991-1997 Silicon Graphics, Inc.
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
30
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
31 Permission to use, copy, modify, distribute, and sell this software and
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
32 its documentation for any purpose is hereby granted without fee, provided
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
33 that (i) the above copyright notices and this permission notice appear in
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
34 all copies of the software and related documentation, and (ii) the names of
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
35 Sam Leffler and Silicon Graphics may not be used in any advertising or
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
36 publicity relating to the software without the specific, prior written
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
37 permission of Sam Leffler and Silicon Graphics. */
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
38
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
39 /* Quantizing code based off of the paper
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
40 Color Image Quantization for Frame Buffer Display, Paul Heckbert,
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
41 Siggraph '82 proceedings, pp. 297-307 */
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
42
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
43 #include <config.h>
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
44 #include "lisp.h"
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
45 #include "imgproc.h"
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
46
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
47 static void
647
b39c14581166 [xemacs-hg @ 2001-08-13 04:45:47 by ben]
ben
parents: 428
diff changeset
48 get_histogram(quant_table *qt, UChar_Binary *pic,
428
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
49 int width, int height, Colorbox* box)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
50 {
647
b39c14581166 [xemacs-hg @ 2001-08-13 04:45:47 by ben]
ben
parents: 428
diff changeset
51 register UChar_Binary *inptr;
428
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
52 register int red, green, blue;
647
b39c14581166 [xemacs-hg @ 2001-08-13 04:45:47 by ben]
ben
parents: 428
diff changeset
53 register int j, i;
428
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
54
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
55 box->rmin = box->gmin = box->bmin = 999;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
56 box->rmax = box->gmax = box->bmax = -1;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
57 box->total = width * height;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
58
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
59 inptr = pic;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
60 for (i = 0; i < height; i++)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
61 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
62 for (j = width; j-- > 0;)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
63 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
64 red = *inptr++ >> COLOR_SHIFT;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
65 green = *inptr++ >> COLOR_SHIFT;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
66 blue = *inptr++ >> COLOR_SHIFT;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
67 if (red < box->rmin)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
68 box->rmin = red;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
69 if (red > box->rmax)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
70 box->rmax = red;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
71 if (green < box->gmin)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
72 box->gmin = green;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
73 if (green > box->gmax)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
74 box->gmax = green;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
75 if (blue < box->bmin)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
76 box->bmin = blue;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
77 if (blue > box->bmax)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
78 box->bmax = blue;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
79 qt->histogram[red][green][blue]++;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
80 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
81 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
82 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
83
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
84 static Colorbox *
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
85 largest_box(quant_table *qt)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
86 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
87 register Colorbox *p, *b;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
88 register int size;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
89
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
90 b = NULL;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
91 size = -1;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
92 for (p = qt->usedboxes; p != NULL; p = p->next)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
93 if ((p->rmax > p->rmin || p->gmax > p->gmin ||
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
94 p->bmax > p->bmin) && p->total > size)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
95 size = (b = p)->total;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
96 return (b);
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
97 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
98
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
99 static void
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
100 shrinkbox(quant_table *qt, Colorbox* box)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
101 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
102 register int *histp, ir, ig, ib;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
103
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
104 if (box->rmax > box->rmin)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
105 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
106 for (ir = box->rmin; ir <= box->rmax; ++ir)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
107 for (ig = box->gmin; ig <= box->gmax; ++ig)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
108 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
109 histp = &(qt->histogram[ir][ig][box->bmin]);
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
110 for (ib = box->bmin; ib <= box->bmax; ++ib)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
111 if (*histp++ != 0)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
112 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
113 box->rmin = ir;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
114 goto have_rmin;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
115 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
116 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
117 have_rmin:
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
118 if (box->rmax > box->rmin)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
119 for (ir = box->rmax; ir >= box->rmin; --ir)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
120 for (ig = box->gmin; ig <= box->gmax; ++ig)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
121 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
122 histp = &(qt->histogram[ir][ig][box->bmin]);
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
123 ib = box->bmin;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
124 for (; ib <= box->bmax; ++ib)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
125 if (*histp++ != 0)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
126 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
127 box->rmax = ir;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
128 goto have_rmax;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
129 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
130 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
131 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
132 have_rmax:
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
133 if (box->gmax > box->gmin)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
134 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
135 for (ig = box->gmin; ig <= box->gmax; ++ig)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
136 for (ir = box->rmin; ir <= box->rmax; ++ir)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
137 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
138 histp = &(qt->histogram[ir][ig][box->bmin]);
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
139 for (ib = box->bmin; ib <= box->bmax; ++ib)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
140 if (*histp++ != 0)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
141 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
142 box->gmin = ig;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
143 goto have_gmin;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
144 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
145 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
146 have_gmin:
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
147 if (box->gmax > box->gmin)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
148 for (ig = box->gmax; ig >= box->gmin; --ig)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
149 for (ir = box->rmin; ir <= box->rmax; ++ir)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
150 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
151 histp = &(qt->histogram[ir][ig][box->bmin]);
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
152 ib = box->bmin;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
153 for (; ib <= box->bmax; ++ib)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
154 if (*histp++ != 0)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
155 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
156 box->gmax = ig;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
157 goto have_gmax;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
158 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
159 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
160 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
161 have_gmax:
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
162 if (box->bmax > box->bmin)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
163 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
164 for (ib = box->bmin; ib <= box->bmax; ++ib)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
165 for (ir = box->rmin; ir <= box->rmax; ++ir)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
166 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
167 histp = &(qt->histogram[ir][box->gmin][ib]);
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
168 for (ig = box->gmin; ig <= box->gmax; ++ig)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
169 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
170 if (*histp != 0)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
171 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
172 box->bmin = ib;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
173 goto have_bmin;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
174 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
175 histp += B_LEN;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
176 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
177 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
178 have_bmin:
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
179 if (box->bmax > box->bmin)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
180 for (ib = box->bmax; ib >= box->bmin; --ib)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
181 for (ir = box->rmin; ir <= box->rmax; ++ir)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
182 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
183 histp = &(qt->histogram[ir][box->gmin][ib]);
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
184 ig = box->gmin;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
185 for (; ig <= box->gmax; ++ig)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
186 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
187 if (*histp != 0)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
188 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
189 box->bmax = ib;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
190 goto have_bmax;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
191 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
192 histp += B_LEN;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
193 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
194 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
195 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
196 have_bmax:
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
197 ;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
198 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
199
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
200 static void
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
201 splitbox(quant_table *qt, Colorbox* ptr)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
202 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
203 int hist2[B_LEN];
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
204 int first = 0, last = 0;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
205 register Colorbox *new;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
206 register int *iptr, *histp;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
207 register int i, j;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
208 register int ir,ig,ib;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
209 register int sum, sum1, sum2;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
210 enum { RED, GREEN, BLUE } axis;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
211
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
212 /*
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
213 * See which axis is the largest, do a histogram along that
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
214 * axis. Split at median point. Contract both new boxes to
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
215 * fit points and return
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
216 */
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
217 i = ptr->rmax - ptr->rmin;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
218 if (i >= ptr->gmax - ptr->gmin && i >= ptr->bmax - ptr->bmin)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
219 axis = RED;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
220 else if (ptr->gmax - ptr->gmin >= ptr->bmax - ptr->bmin)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
221 axis = GREEN;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
222 else
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
223 axis = BLUE;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
224 /* get histogram along longest axis */
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
225 switch (axis)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
226 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
227 case RED:
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
228 histp = &hist2[ptr->rmin];
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
229 for (ir = ptr->rmin; ir <= ptr->rmax; ++ir)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
230 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
231 *histp = 0;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
232 for (ig = ptr->gmin; ig <= ptr->gmax; ++ig)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
233 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
234 iptr = &(qt->histogram[ir][ig][ptr->bmin]);
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
235 for (ib = ptr->bmin; ib <= ptr->bmax; ++ib)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
236 *histp += *iptr++;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
237 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
238 histp++;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
239 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
240 first = ptr->rmin;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
241 last = ptr->rmax;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
242 break;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
243 case GREEN:
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
244 histp = &hist2[ptr->gmin];
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
245 for (ig = ptr->gmin; ig <= ptr->gmax; ++ig)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
246 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
247 *histp = 0;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
248 for (ir = ptr->rmin; ir <= ptr->rmax; ++ir)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
249 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
250 iptr = &(qt->histogram[ir][ig][ptr->bmin]);
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
251 for (ib = ptr->bmin; ib <= ptr->bmax; ++ib)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
252 *histp += *iptr++;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
253 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
254 histp++;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
255 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
256 first = ptr->gmin;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
257 last = ptr->gmax;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
258 break;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
259 case BLUE:
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
260 histp = &hist2[ptr->bmin];
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
261 for (ib = ptr->bmin; ib <= ptr->bmax; ++ib)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
262 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
263 *histp = 0;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
264 for (ir = ptr->rmin; ir <= ptr->rmax; ++ir)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
265 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
266 iptr = &(qt->histogram[ir][ptr->gmin][ib]);
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
267 for (ig = ptr->gmin; ig <= ptr->gmax; ++ig)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
268 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
269 *histp += *iptr;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
270 iptr += B_LEN;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
271 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
272 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
273 histp++;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
274 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
275 first = ptr->bmin;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
276 last = ptr->bmax;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
277 break;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
278 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
279 /* find median point */
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
280 sum2 = ptr->total / 2;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
281 histp = &hist2[first];
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
282 sum = 0;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
283 for (i = first; i <= last && (sum += *histp++) < sum2; ++i)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
284 ;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
285 if (i == first)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
286 i++;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
287
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
288 /* Create new box, re-allocate points */
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
289 new = qt->freeboxes;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
290 qt->freeboxes = new->next;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
291 if (qt->freeboxes)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
292 qt->freeboxes->prev = NULL;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
293 if (qt->usedboxes)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
294 qt->usedboxes->prev = new;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
295 new->next = qt->usedboxes;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
296 qt->usedboxes = new;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
297
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
298 histp = &hist2[first];
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
299 for (sum1 = 0, j = first; j < i; j++)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
300 sum1 += *histp++;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
301 for (sum2 = 0, j = i; j <= last; j++)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
302 sum2 += *histp++;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
303 new->total = sum1;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
304 ptr->total = sum2;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
305
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
306 new->rmin = ptr->rmin;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
307 new->rmax = ptr->rmax;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
308 new->gmin = ptr->gmin;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
309 new->gmax = ptr->gmax;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
310 new->bmin = ptr->bmin;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
311 new->bmax = ptr->bmax;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
312 switch (axis)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
313 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
314 case RED:
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
315 new->rmax = i-1;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
316 ptr->rmin = i;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
317 break;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
318 case GREEN:
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
319 new->gmax = i-1;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
320 ptr->gmin = i;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
321 break;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
322 case BLUE:
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
323 new->bmax = i-1;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
324 ptr->bmin = i;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
325 break;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
326 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
327 shrinkbox (qt, new);
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
328 shrinkbox (qt, ptr);
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
329 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
330
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
331
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
332 static C_cell *
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
333 create_colorcell(quant_table *qt, int num_colors, int red, int green, int blue)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
334 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
335 register int ir, ig, ib, i;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
336 register C_cell *ptr;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
337 int mindist, next_n;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
338 register int tmp, dist, n;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
339
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
340 ir = red >> (COLOR_DEPTH-C_DEPTH);
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
341 ig = green >> (COLOR_DEPTH-C_DEPTH);
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
342 ib = blue >> (COLOR_DEPTH-C_DEPTH);
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
343 ptr = (C_cell *)xmalloc(sizeof (C_cell));
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
344 *(qt->ColorCells + ir*C_LEN*C_LEN + ig*C_LEN + ib) = ptr;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
345 ptr->num_ents = 0;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
346
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
347 /*
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
348 * Step 1: find all colors inside this cell, while we're at
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
349 * it, find distance of centermost point to furthest corner
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
350 */
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
351 mindist = 99999999;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
352 for (i = 0; i < num_colors; ++i)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
353 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
354 if (qt->rm[i]>>(COLOR_DEPTH-C_DEPTH) != ir ||
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
355 qt->gm[i]>>(COLOR_DEPTH-C_DEPTH) != ig ||
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
356 qt->bm[i]>>(COLOR_DEPTH-C_DEPTH) != ib)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
357 continue;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
358 ptr->entries[ptr->num_ents][0] = i;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
359 ptr->entries[ptr->num_ents][1] = 0;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
360 ++ptr->num_ents;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
361 tmp = qt->rm[i] - red;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
362 if (tmp < (MAX_COLOR/C_LEN/2))
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
363 tmp = MAX_COLOR/C_LEN-1 - tmp;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
364 dist = tmp*tmp;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
365 tmp = qt->gm[i] - green;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
366 if (tmp < (MAX_COLOR/C_LEN/2))
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
367 tmp = MAX_COLOR/C_LEN-1 - tmp;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
368 dist += tmp*tmp;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
369 tmp = qt->bm[i] - blue;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
370 if (tmp < (MAX_COLOR/C_LEN/2))
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
371 tmp = MAX_COLOR/C_LEN-1 - tmp;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
372 dist += tmp*tmp;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
373 if (dist < mindist)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
374 mindist = dist;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
375 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
376
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
377 /*
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
378 * Step 3: find all points within that distance to cell.
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
379 */
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
380 for (i = 0; i < num_colors; ++i)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
381 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
382 if (qt->rm[i] >> (COLOR_DEPTH-C_DEPTH) == ir &&
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
383 qt->gm[i] >> (COLOR_DEPTH-C_DEPTH) == ig &&
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
384 qt->bm[i] >> (COLOR_DEPTH-C_DEPTH) == ib)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
385 continue;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
386 dist = 0;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
387 if ((tmp = red - qt->rm[i]) > 0 ||
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
388 (tmp = qt->rm[i] - (red + MAX_COLOR/C_LEN-1)) > 0 )
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
389 dist += tmp*tmp;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
390 if ((tmp = green - qt->gm[i]) > 0 ||
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
391 (tmp = qt->gm[i] - (green + MAX_COLOR/C_LEN-1)) > 0 )
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
392 dist += tmp*tmp;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
393 if ((tmp = blue - qt->bm[i]) > 0 ||
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
394 (tmp = qt->bm[i] - (blue + MAX_COLOR/C_LEN-1)) > 0 )
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
395 dist += tmp*tmp;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
396 if (dist < mindist)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
397 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
398 ptr->entries[ptr->num_ents][0] = i;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
399 ptr->entries[ptr->num_ents][1] = dist;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
400 ++ptr->num_ents;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
401 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
402 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
403
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
404 /*
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
405 * Sort color cells by distance, use cheap exchange sort
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
406 */
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
407 for (n = ptr->num_ents - 1; n > 0; n = next_n)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
408 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
409 next_n = 0;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
410 for (i = 0; i < n; ++i)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
411 if (ptr->entries[i][1] > ptr->entries[i+1][1])
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
412 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
413 tmp = ptr->entries[i][0];
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
414 ptr->entries[i][0] = ptr->entries[i+1][0];
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
415 ptr->entries[i+1][0] = tmp;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
416 tmp = ptr->entries[i][1];
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
417 ptr->entries[i][1] = ptr->entries[i+1][1];
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
418 ptr->entries[i+1][1] = tmp;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
419 next_n = i;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
420 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
421 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
422 return (ptr);
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
423 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
424
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
425 static int
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
426 map_colortable(quant_table *qt, int num_colors)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
427 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
428 register int *histp = &(qt->histogram[0][0][0]);
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
429 register C_cell *cell;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
430 register int j, tmp, d2, dist;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
431 int ir, ig, ib, i;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
432
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
433 for (ir = 0; ir < B_LEN; ++ir)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
434 for (ig = 0; ig < B_LEN; ++ig)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
435 for (ib = 0; ib < B_LEN; ++ib, histp++)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
436 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
437 if (*histp == 0)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
438 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
439 *histp = -1;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
440 continue;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
441 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
442 cell = *(qt->ColorCells +
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
443 (((ir>>(B_DEPTH-C_DEPTH)) << C_DEPTH*2) +
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
444 ((ig>>(B_DEPTH-C_DEPTH)) << C_DEPTH) +
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
445 (ib>>(B_DEPTH-C_DEPTH))));
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
446 if (cell == NULL )
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
447 cell = create_colorcell (qt, num_colors,
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
448 ir << COLOR_SHIFT,
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
449 ig << COLOR_SHIFT,
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
450 ib << COLOR_SHIFT);
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
451 if (cell == NULL) /* memory exhausted! punt! */
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
452 return -1;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
453 dist = 9999999;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
454 for (i = 0; i < cell->num_ents &&
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
455 dist > cell->entries[i][1]; ++i)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
456 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
457 j = cell->entries[i][0];
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
458 d2 = qt->rm[j] - (ir << COLOR_SHIFT);
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
459 d2 *= d2;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
460 tmp = qt->gm[j] - (ig << COLOR_SHIFT);
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
461 d2 += tmp*tmp;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
462 tmp = qt->bm[j] - (ib << COLOR_SHIFT);
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
463 d2 += tmp*tmp;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
464 if (d2 < dist)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
465 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
466 dist = d2;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
467 *histp = j;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
468 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
469 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
470 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
471 return 0;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
472 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
473
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
474 quant_table *
647
b39c14581166 [xemacs-hg @ 2001-08-13 04:45:47 by ben]
ben
parents: 428
diff changeset
475 build_EImage_quantable(UChar_Binary *eimage, int width, int height, int num_colors)
428
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
476 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
477 quant_table *qt;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
478 Colorbox *box_list, *ptr;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
479 int i,res;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
480
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
481 qt = (quant_table*)xmalloc_and_zero (sizeof(quant_table));
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
482 if (qt == NULL) return NULL;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
483
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
484 assert (num_colors < 257 && num_colors > 2);
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
485 /*
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
486 * STEP 1: create empty boxes
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
487 */
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
488 qt->usedboxes = NULL;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
489 box_list = qt->freeboxes = (Colorbox *)xmalloc (num_colors*sizeof (Colorbox));
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
490 qt->freeboxes[0].next = &(qt->freeboxes[1]);
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
491 qt->freeboxes[0].prev = NULL;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
492 for (i = 1; i < num_colors-1; ++i)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
493 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
494 qt->freeboxes[i].next = &(qt->freeboxes[i+1]);
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
495 qt->freeboxes[i].prev = &(qt->freeboxes[i-1]);
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
496 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
497 qt->freeboxes[num_colors-1].next = NULL;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
498 qt->freeboxes[num_colors-1].prev = &(qt->freeboxes[num_colors-2]);
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
499
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
500 /*
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
501 * STEP 2: get histogram, initialize first box
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
502 */
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
503 ptr = qt->freeboxes;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
504 qt->freeboxes = ptr->next;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
505 if (qt->freeboxes)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
506 qt->freeboxes->prev = NULL;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
507 ptr->next = qt->usedboxes;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
508 qt->usedboxes = ptr;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
509 if (ptr->next)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
510 ptr->next->prev = ptr;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
511 get_histogram (qt, eimage, width, height, ptr);
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
512
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
513 /*
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
514 * STEP 3: continually subdivide boxes until no more free
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
515 * boxes remain or until all colors assigned.
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
516 */
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
517 while (qt->freeboxes != NULL)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
518 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
519 ptr = largest_box(qt);
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
520 if (ptr != NULL)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
521 splitbox (qt, ptr);
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
522 else
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
523 qt->freeboxes = NULL;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
524 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
525
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
526 /*
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
527 * STEP 4: assign colors to all boxes
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
528 */
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
529 for (i = 0, ptr = qt->usedboxes; ptr != NULL; ++i, ptr = ptr->next)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
530 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
531 qt->rm[i] = ((ptr->rmin + ptr->rmax) << COLOR_SHIFT) / 2;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
532 qt->gm[i] = ((ptr->gmin + ptr->gmax) << COLOR_SHIFT) / 2;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
533 qt->bm[i] = ((ptr->bmin + ptr->bmax) << COLOR_SHIFT) / 2;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
534 qt->um[i] = ptr->total;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
535 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
536 qt->num_active_colors = i;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
537
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
538 /* We're done with the boxes now */
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
539 xfree (box_list);
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
540 qt->freeboxes = qt->usedboxes = NULL;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
541
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
542 /*
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
543 * STEP 5: scan histogram and map all values to closest color
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
544 */
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
545 /* 5a: create cell list as described in Heckbert */
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
546 qt->ColorCells = (C_cell **)xmalloc_and_zero (C_LEN*C_LEN*C_LEN*sizeof (C_cell*));
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
547 /* 5b: create mapping from truncated pixel space to color
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
548 table entries */
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
549 res = map_colortable (qt, num_colors);
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
550
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
551 /* 5c: done with ColorCells */
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
552 for (i = 0; i < C_LEN*C_LEN*C_LEN; i++)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
553 if (qt->ColorCells[i]) xfree (qt->ColorCells[i]);
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
554 xfree (qt->ColorCells);
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
555
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
556 if (res)
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
557 {
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
558 /* we failed in memory allocation, so clean up an leave */
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
559 xfree(qt);
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
560 return NULL;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
561 }
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
562
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
563 return qt;
3ecd8885ac67 Import from CVS: tag r21-2-22
cvs
parents:
diff changeset
564 }