1 28994509 2004-04-21 devnull #include <u.h>
2 28994509 2004-04-21 devnull #include <libc.h>
3 28994509 2004-04-21 devnull #include <bio.h>
4 28994509 2004-04-21 devnull #include <draw.h>
5 28994509 2004-04-21 devnull #include "imagefile.h"
8 28994509 2004-04-21 devnull c1 = 2871, /* 1.402 * 2048 */
9 28994509 2004-04-21 devnull c2 = 705, /* 0.34414 * 2048 */
10 28994509 2004-04-21 devnull c3 = 1463, /* 0.71414 * 2048 */
11 cbeb0b26 2006-04-01 devnull c4 = 3629 /* 1.772 * 2048 */
14 28994509 2004-04-21 devnull Rawimage*
15 28994509 2004-04-21 devnull totruecolor(Rawimage *i, int chandesc)
17 28994509 2004-04-21 devnull int j, k;
18 28994509 2004-04-21 devnull Rawimage *im;
19 28994509 2004-04-21 devnull char err[ERRMAX];
20 28994509 2004-04-21 devnull uchar *rp, *gp, *bp, *cmap, *inp, *outp, cmap1[3*256];
21 28994509 2004-04-21 devnull int r, g, b, Y, Cr, Cb;
23 28994509 2004-04-21 devnull if(chandesc!=CY && chandesc!=CRGB24)
24 28994509 2004-04-21 devnull return _remaperror("remap: can't convert to chandesc %d", chandesc);
26 28994509 2004-04-21 devnull err[0] = '\0';
27 28994509 2004-04-21 devnull errstr(err, sizeof err); /* throw it away */
28 28994509 2004-04-21 devnull im = malloc(sizeof(Rawimage));
29 28994509 2004-04-21 devnull if(im == nil)
30 28994509 2004-04-21 devnull return nil;
31 28994509 2004-04-21 devnull memset(im, 0, sizeof(Rawimage));
32 28994509 2004-04-21 devnull if(chandesc == CY)
33 28994509 2004-04-21 devnull im->chanlen = i->chanlen;
35 28994509 2004-04-21 devnull im->chanlen = 3*i->chanlen;
36 28994509 2004-04-21 devnull im->chandesc = chandesc;
37 28994509 2004-04-21 devnull im->chans[0] = malloc(im->chanlen);
38 28994509 2004-04-21 devnull if(im->chans[0] == nil){
39 28994509 2004-04-21 devnull free(im);
40 28994509 2004-04-21 devnull return nil;
42 28994509 2004-04-21 devnull im->r = i->r;
43 28994509 2004-04-21 devnull im->nchans = 1;
45 28994509 2004-04-21 devnull cmap = i->cmap;
47 28994509 2004-04-21 devnull outp = im->chans[0];
49 28994509 2004-04-21 devnull switch(i->chandesc){
51 28994509 2004-04-21 devnull return _remaperror("remap: can't recognize channel type %d", i->chandesc);
53 28994509 2004-04-21 devnull if(i->nchans != 1)
54 28994509 2004-04-21 devnull return _remaperror("remap: Y image has %d chans", i->nchans);
55 28994509 2004-04-21 devnull if(chandesc == CY){
56 28994509 2004-04-21 devnull memmove(im->chans[0], i->chans[0], i->chanlen);
59 28994509 2004-04-21 devnull /* convert to three color */
60 28994509 2004-04-21 devnull inp = i->chans[0];
61 28994509 2004-04-21 devnull for(j=0; j<i->chanlen; j++){
62 28994509 2004-04-21 devnull k = *inp++;
63 28994509 2004-04-21 devnull *outp++ = k;
64 28994509 2004-04-21 devnull *outp++ = k;
65 28994509 2004-04-21 devnull *outp++ = k;
69 28994509 2004-04-21 devnull case CRGB1:
70 28994509 2004-04-21 devnull if(cmap == nil)
71 28994509 2004-04-21 devnull return _remaperror("remap: image has no color map");
72 28994509 2004-04-21 devnull if(i->nchans != 1)
73 28994509 2004-04-21 devnull return _remaperror("remap: can't handle nchans %d", i->nchans);
74 28994509 2004-04-21 devnull for(j=1; j<=8; j++)
75 28994509 2004-04-21 devnull if(i->cmaplen == 3*(1<<j))
77 28994509 2004-04-21 devnull if(j > 8)
78 28994509 2004-04-21 devnull return _remaperror("remap: can't do colormap size 3*%d", i->cmaplen/3);
79 28994509 2004-04-21 devnull if(i->cmaplen != 3*256){
80 28994509 2004-04-21 devnull /* to avoid a range check in loop below, make a full-size cmap */
81 28994509 2004-04-21 devnull memmove(cmap1, cmap, i->cmaplen);
82 28994509 2004-04-21 devnull cmap = cmap1;
84 28994509 2004-04-21 devnull inp = i->chans[0];
85 28994509 2004-04-21 devnull if(chandesc == CY){
86 28994509 2004-04-21 devnull for(j=0; j<i->chanlen; j++){
87 28994509 2004-04-21 devnull k = *inp++;
88 28994509 2004-04-21 devnull r = cmap[3*k+2];
89 28994509 2004-04-21 devnull g = cmap[3*k+1];
90 28994509 2004-04-21 devnull b = cmap[3*k+0];
91 28994509 2004-04-21 devnull r = (2125*r + 7154*g + 721*b)/10000; /* Poynton page 84 */
92 28994509 2004-04-21 devnull *outp++ = r;
95 28994509 2004-04-21 devnull for(j=0; j<i->chanlen; j++){
96 28994509 2004-04-21 devnull k = *inp++;
97 28994509 2004-04-21 devnull *outp++ = cmap[3*k+2];
98 28994509 2004-04-21 devnull *outp++ = cmap[3*k+1];
99 28994509 2004-04-21 devnull *outp++ = cmap[3*k+0];
104 28994509 2004-04-21 devnull case CRGB:
105 28994509 2004-04-21 devnull if(i->nchans != 3)
106 28994509 2004-04-21 devnull return _remaperror("remap: can't handle nchans %d", i->nchans);
107 28994509 2004-04-21 devnull rp = i->chans[0];
108 28994509 2004-04-21 devnull gp = i->chans[1];
109 28994509 2004-04-21 devnull bp = i->chans[2];
110 28994509 2004-04-21 devnull if(chandesc == CY){
111 28994509 2004-04-21 devnull for(j=0; j<i->chanlen; j++){
112 28994509 2004-04-21 devnull r = *bp++;
113 28994509 2004-04-21 devnull g = *gp++;
114 28994509 2004-04-21 devnull b = *rp++;
115 28994509 2004-04-21 devnull r = (2125*r + 7154*g + 721*b)/10000; /* Poynton page 84 */
116 28994509 2004-04-21 devnull *outp++ = r;
119 28994509 2004-04-21 devnull for(j=0; j<i->chanlen; j++){
120 28994509 2004-04-21 devnull *outp++ = *bp++;
121 28994509 2004-04-21 devnull *outp++ = *gp++;
122 28994509 2004-04-21 devnull *outp++ = *rp++;
126 28994509 2004-04-21 devnull case CYCbCr:
127 28994509 2004-04-21 devnull if(i->nchans != 3)
128 28994509 2004-04-21 devnull return _remaperror("remap: can't handle nchans %d", i->nchans);
129 28994509 2004-04-21 devnull rp = i->chans[0];
130 28994509 2004-04-21 devnull gp = i->chans[1];
131 28994509 2004-04-21 devnull bp = i->chans[2];
132 28994509 2004-04-21 devnull for(j=0; j<i->chanlen; j++){
133 28994509 2004-04-21 devnull Y = *rp++ << 11;
134 28994509 2004-04-21 devnull Cb = *gp++ - 128;
135 28994509 2004-04-21 devnull Cr = *bp++ - 128;
136 28994509 2004-04-21 devnull r = (Y+c1*Cr) >> 11;
137 28994509 2004-04-21 devnull g = (Y-c2*Cb-c3*Cr) >> 11;
138 28994509 2004-04-21 devnull b = (Y+c4*Cb) >> 11;
139 28994509 2004-04-21 devnull if(r < 0)
141 28994509 2004-04-21 devnull if(r > 255)
142 28994509 2004-04-21 devnull r = 255;
143 28994509 2004-04-21 devnull if(g < 0)
145 28994509 2004-04-21 devnull if(g > 255)
146 28994509 2004-04-21 devnull g = 255;
147 28994509 2004-04-21 devnull if(b < 0)
149 28994509 2004-04-21 devnull if(b > 255)
150 28994509 2004-04-21 devnull b = 255;
151 28994509 2004-04-21 devnull if(chandesc == CY){
152 28994509 2004-04-21 devnull r = (2125*r + 7154*g + 721*b)/10000;
153 28994509 2004-04-21 devnull *outp++ = r;
155 28994509 2004-04-21 devnull *outp++ = b;
156 28994509 2004-04-21 devnull *outp++ = g;
157 28994509 2004-04-21 devnull *outp++ = r;
162 28994509 2004-04-21 devnull return im;
