* image.cpp (formerly sampled.cpp)

/*
* image.cpp (formerly sampled.cpp)
* by [email protected] at Wed Feb 27 09:26:05 CET 2002
*/

#ifdef __GNUC__
#ifndef __clang__
#pragma implementation
#endif
#endif

#include "image.hpp"
#include "error.hpp"
#include <string.h> /* strlen() */
#include "gensio.hpp"

/* --- 4-byte hashing */

#if SIZEOF_INT>=4
typedef unsigned int u32_t;
typedef signed int s32_t;
#else
typedef unsigned long u32_t;
typedef signed long s32_t;
#endif

/* vvv Dat: moved these out of Hash46 to pacify VC6.0 */
const unsigned M=1409;
/** Number of _value_ data bytes (they are not hashed) */
const unsigned D=2;
/** Size of each tuple in the array `t' */
const unsigned HD=4+D;
/** A tuple is considered free iff its first byte equals FREE */
const unsigned char FREE=255;

/**
* M=1409
* h(K)=K%1409
* h(i,K)=-i*(1+(K%1408)) (i in 0..1408)
* K=[a,r,g,b]=(a<<24)+(r<<16)+(g<<8)+b (4 bytes)
* h(K)=(253*a+722*r+(g<<8)+b)%1409
* h(i,K)=-i*(1+(896*a+768*r+(g<<8)+b)%1408)
*
* -- No dynamic growing or re-hashing. Useful for hashing colormap palettes
* with a maximum size of 256.
* -- Deleting not supported.
*
* Implementation: Use 32 bit integers for calculation.
*
* Imp: remove unused attr `size'
*/
class Hash46 {
public:
/** Creates an empty hash. */
Hash46();
inline unsigned getSize() const { return size; }
inline unsigned getLength() const { return size; }
inline unsigned getMaxSize() const { return M; }
inline bool isFull() const { return size==M; }
/** @return NULLP or the pointer to a tuple */
inline unsigned char* lookup(unsigned char k[4]) {
unsigned char *ret=walk(k);
return ret==NULLP || *ret==FREE ? (unsigned char*)NULLP : ret;
}
/** Can be called only if !isFull()
* @return NULL if isFull() and not found; otherwise: pointer to the tuple
* found or the place to which the insert can take place.
*/
unsigned char* walk(unsigned char k[4]);
protected:
/** Number of non-free tuples in the hash. */
unsigned size;
unsigned char t[M*HD];
};

Hash46::Hash46(): size(0) {
memset(t, FREE, sizeof(t));
}

unsigned char *Hash46::walk(unsigned char k[4]) {
u32_t hk, hik;
hk=HD*(((((u32_t)1<<24)%M)*k[0]+(((u32_t)1<<16)%M)*k[1]+
(((u32_t)1<< 8)%M)*k[2]+k[3])%M);
hik=HD*(1+((((u32_t)1<<24)%(M-1))*k[0]+(((u32_t)1<<16)%(M-1))*k[1]+
(((u32_t)1<< 8)%(M-1))*k[2]+k[3])%(M-1));
/* fprintf(stderr, "hk=%u hik=%u\n", hk, hik); */
register unsigned char *p=t+hk;
unsigned i=M;
/* fprintf(stderr, "?? %02x %02x %02x %02x\n", k[0], k[1], k[2], k[3]); */
do {
/* fprintf(stderr, "examining %02x %02x %02x %02x %d\n", p[0], p[1], p[2], p[3], (k[0]=p[0] && k[1]==p[1] && k[2]==p[2])); */
if (*p==FREE || (k[0]==p[0] && k[1]==p[1] && k[2]==p[2] && k[3]==p[3])) return p;
/* ^^^ huge == BUGFIX at Sun Apr 14 00:16:59 CEST 2002 */
if (hk>=hik) { hk-=hik; p-=hik; }
else { hk+=M*HD-hik; p+=M*HD-hik; }
} while (--i!=0);
/* fprintf(stderr, "full\n"); */
return (unsigned char*)NULLP;
}

/* --- */

const unsigned char Image::Sampled::cs2cpp[6]= { 0, 1, 3, 3, 4, 4 };
char const *Image::Sampled::cs2devcs(unsigned char cs) {
static const char *names[]={ (char*)NULLP, "Gray", "RGB", "RGB", "CMYK", "CMYK" };
return cs>=1 && cs<=5 ? names[cs] : (char*)NULLP;
}

static void fatal_image_too_large() {
Error::sev(Error::EERROR) << "Image: Image too large." << (Error*)0;
}

static slen_t multiply_check(slen_t a, slen_t b) {
const slen_t result = a * b;
/* Check for overflow. Works only if everything is unsigned. */
if (result / a != b) fatal_image_too_large();
return result;
}

static slen_t multiply_check(slen_t a, slen_t b, slen_t c) {
return multiply_check(multiply_check(a, b), c);
}

static slen_t add_check(slen_t a, slen_t b) {
/* Check for overflow. Works only if everything is unsigned. */
if (b > (slen_t)-1 - a) fatal_image_too_large();
return a + b;
}

#if 0
static slen_t add_check(slen_t a, slen_t b, slen_t c) {
return add_check(add_check(a, b), c);
}
#endif

static slen_t add_check(slen_t a, slen_t b, slen_t c, slen_t d) {
return add_check(add_check(a, b), add_check(c, d));
}

void Image::Sampled::init(slen_t l_comment, slen_t l_header, dimen_t wd_, dimen_t ht_,
/* ^^^ 24 is required for /Transparent in out_tiff_work */
unsigned char bpc_, unsigned char ty_, unsigned char cpp_) {
/* Even if we continue from here, most probably we'll reach
* ``sam2p.yes: Error: applyProfile: invalid combination, no applicable OutputRule''.
* So more work is needed to support output images of size 0.
*/
if (wd_ <= 0 || ht_ <= 0) Error::sev(Error::EERROR) << "Image: Image of size 0." << (Error*)0;
bpc=bpc_;
ty=ty_;
wd=wd_;
ht=ht_;
cpp=cpp_;
// pred=1;
transpc=0x1000000UL; /* Dat: this means: no transparent color */
const slen_t rlens = add_check(multiply_check(bpc_, cpp_, wd_), 7) >> 3;
rlen = rlens;
if (rlen != rlens) fatal_image_too_large();
beg=new char[len=add_check(l_comment, l_header, multiply_check(rlen, ht_), bpc)];
rowbeg=(headp=const_cast<char*>(beg)+l_comment)+l_header;
trail=const_cast<char*>(beg)+len-bpc;
memset(trail, 0, bpc);
}

Image::Gray* Image::Sampled::toGray0(unsigned char bpc_) {
unsigned char *crow=new unsigned char[wd*3+7*3], *p, *pend;
Image::Gray *img=new Image::Gray(wd, ht, bpc_);
unsigned char *outp=(unsigned char*)img->getRowbeg();
dimen_t htc;
memset(crow+wd*3, '\0', 7*3); /* *3 BUGFIX at Tue Jan 18 17:04:15 CET 2005 */
unsigned i;
/* Dat: not optimising for minimal rounding error since caller should ensure
* that there is no such error at all.
*/
if (bpc_==1) {
assert(img->getBpc()==1);
for (htc=0;htc<ht;htc++) {
copyRGBRow((char*)crow, htc);
for (p=crow-3,pend=p+wd*3; p<pend; ) {
i =(*(p+=3)!=0)<<7; i|=(*(p+=3)!=0)<<6; i|=(*(p+=3)!=0)<<5; i|=(*(p+=3)!=0)<<4;
i|=(*(p+=3)!=0)<<3; i|=(*(p+=3)!=0)<<2; i|=(*(p+=3)!=0)<<1; i|=(*(p+=3)!=0);
*outp++=i;
}
}
} else if (bpc_==2) {
for (htc=0;htc<ht;htc++) {
copyRGBRow((char*)crow, htc);
for (p=crow-3,pend=p+wd*3; p<pend; ) {
i =(*(p+=3)/85)<<6; i|=(*(p+=3)/85)<<4; i|=(*(p+=3)/85)<<2; i|=(*(p+=3)/85);
*outp++=i;
}
}
} else if (bpc_==4) {
for (htc=0;htc<ht;htc++) {
copyRGBRow((char*)crow, htc);
for (p=crow-3,pend=p+wd*3; p<pend; ) {
i =(*(p+=3)/17)<<4; i|=(*(p+=3)/17);
*outp++=i;
}
}
} else if (bpc_==8) {
for (htc=0;htc<ht;htc++) {
copyRGBRow((char*)crow, htc);
for (p=crow-3,pend=p+wd*3; p!=pend; ) {
*outp++=*(p+=3);
}
}
} else assert(0);
delete [] crow;
return img;
}
Image::RGB* Image::Sampled::toRGB0(unsigned char bpc_) {
unsigned char *crow=new unsigned char[wd*3+7], *p, *pend;
Image::RGB *img=new Image::RGB(wd, ht, bpc_);
unsigned char *outp=(unsigned char*)img->getRowbeg();
dimen_t htc;
memset(crow+wd*3, '\0', 7);
unsigned i;
/* Dat: not optimising for minimal rounding error since caller should ensure
* that there is no such error at all.
*/
if (bpc_==1) {
for (htc=0;htc<ht;htc++) {
copyRGBRow((char*)crow, htc);
for (p=crow,pend=crow+wd*3; p<pend; ) {
i =(*p++!=0)<<7; i|=(*p++!=0)<<6; i|=(*p++!=0)<<5; i|=(*p++!=0)<<4;
i|=(*p++!=0)<<3; i|=(*p++!=0)<<2; i|=(*p++!=0)<<1; i|=(*p++!=0);
*outp++=i;
}
}
} else if (bpc_==2) {
for (htc=0;htc<ht;htc++) {
copyRGBRow((char*)crow, htc);
for (p=crow,pend=crow+wd*3; p<pend; ) {
i =(*p++/85)<<6; i|=(*p++/85)<<4; i|=(*p++/85)<<2; i|=(*p++/85);
*outp++=i;
}
}
} else if (bpc_==4) {
for (htc=0;htc<ht;htc++) {
copyRGBRow((char*)crow, htc);
for (p=crow,pend=crow+wd*3; p<pend; ) {
i =(*p++/17)<<4; i|=(*p++/17);
*outp++=i;
}
}
} else if (bpc_==8) {
for (htc=0;htc<ht;htc++) {
copyRGBRow((char*)outp, htc);
outp+=wd*3;
}
} else assert(0);
delete [] crow;
return img;
}

Image::Indexed* Image::Sampled::toIndexed0()/* const*/ {
unsigned char *crow=new unsigned char[wd*3], k[6], *p, *pend, *w;
Image::Indexed *img=new Image::Indexed(wd, ht, 256, 8);
dimen_t htc;
unsigned char *pal=(unsigned char*)img->getHeadp(), *outp=(unsigned char*)img->getRowbeg();
unsigned ncols=0;
Hash46 h;
k[0]=0;
for (htc=0;htc<ht;htc++) {
copyRGBRow((char*)crow, htc);
for (p=crow,pend=crow+wd*3; p!=pend; p+=3) {
memcpy(k+1, p, 3);
w=h.walk(k);
assert(w!=NULL); /* Hash cannot be full since h.M>=256. */
/* fprintf(stderr, "w=%p\n", w); */
if (*w==/*h.*/FREE) {
if (ncols==256) { delete img; delete [] crow; return (Image::Indexed*)NULLP; }
/* ^^^ too many colors; cannot convert image to indexed */
memcpy(w,k,4);
memcpy(pal,k+1,3);
/* fprintf(stderr,"newcol=%02x #%02x%02x%02x\n", k[0], k[1], k[2], k[3]); */
/* fprintf(stderr,"newcol=pal #%02x%02x%02x\n", pal[0], pal[1], pal[2]); */
pal+=3;
*outp++=w[4]=ncols++;
} else { /* a color that we have already seen */
*outp++=w[4];
}
}
}
img->setNcolsMove(ncols);
delete [] crow;

/* Now img is ready. The user should call packPal() to make it even tighter. */
return img;
}

void Image::Sampled::to8mul() {
if (bpc==8) return;
if (wd==0 || ht==0) { bpc=8; return; }

unsigned oldBpc=bpc;
slen_t wdcpp=wd*cpp;
const char *oldBeg=beg;
unsigned char *p=(unsigned char*)rowbeg;

bpc=8;
rlen=wd;
beg=new char[len=rowbeg-oldBeg+rlen*ht+bpc];
headp= const_cast<char*>(beg)+(headp-oldBeg);
rowbeg=const_cast<char*>(beg)+(rowbeg-oldBeg);
trail= const_cast<char*>(beg)+len-bpc;
memcpy(const_cast<char*>(beg), oldBeg, rowbeg-beg);

unsigned char *to=(unsigned char*)rowbeg, *toend;
unsigned int i, j;
Image::Sampled::dimen_t htc;
if (oldBpc==1) {
htc=ht; while (htc--!=0) {
toend=to+((wdcpp)&~7);
while (to!=toend) {
i=*p++;
*to++=(i>>7)*255;
*to++=((i>>6)&1)*255;
*to++=((i>>5)&1)*255;
*to++=((i>>4)&1)*255;
*to++=((i>>3)&1)*255;
*to++=((i>>2)&1)*255;
*to++=((i>>1)&1)*255;
*to++=( i&1)*255;
}
if (0!=(j=(wdcpp)&7)) {
i=*p; /* No mem overrun, even if (wd&7)==0 */
while (j--!=0) { *to++=(i>>7)*255; i<<=1; }
}
}
} else if (oldBpc==2) {
htc=ht; while (htc--!=0) {
toend=to+((wdcpp)&~3);
while (to!=toend) {
i=*p++;
*to++=(i>>6)*85;
*to++=((i>>4)&3)*85;
*to++=((i>>2)&3)*85;
*to++=( i&3)*85;
}
if (0!=(j=(wdcpp)&3)) {
i=*p; /* No mem overrun, even if (wd&7)==0 */
while (j--!=0) { *to++=(i>>6)*85; i<<=2; }
}
}
} else if (oldBpc==4) {
htc=ht; while (htc--!=0) {
toend=to+((wdcpp)&~1);
while (to!=toend) {
i=*p++;
*to++=(i>>4)*17;
*to++=( i&15)*17;
}
if (0!=((wdcpp)&1)) *to++=(*p++>>4)*17;
}
} else assert(0 && "invalid bpc");

delete [] const_cast<char*>(oldBeg);
}

void Image::Sampled::to8nomul() {
if (bpc==8) return;
if (wd==0 || ht==0) { bpc=8; return; }

unsigned oldBpc=bpc;
slen_t wdcpp=wd*cpp;
const char *oldBeg=beg;
unsigned char *p=(unsigned char*)rowbeg;

bpc=8;
rlen=wd;
beg=new char[len=rowbeg-oldBeg+rlen*ht+bpc];
headp= const_cast<char*>(beg)+(headp-oldBeg);
rowbeg=const_cast<char*>(beg)+(rowbeg-oldBeg);
trail= const_cast<char*>(beg)+len-bpc;
memcpy(const_cast<char*>(beg), oldBeg, rowbeg-beg);

unsigned char *to=(unsigned char*)rowbeg, *toend;
unsigned int i, j;
Image::Sampled::dimen_t htc;
if (oldBpc==1) {
htc=ht; while (htc--!=0) {
toend=to+((wdcpp)&~7);
while (to!=toend) {
i=*p++;
*to++=(i>>7);
*to++=((i>>6)&1);
*to++=((i>>5)&1);
*to++=((i>>4)&1);
*to++=((i>>3)&1);
*to++=((i>>2)&1);
*to++=((i>>1)&1);
*to++=( i&1);
}
if (0!=(j=(wdcpp)&7)) {
i=*p++; /* No mem overrun, even if (wd&7)==0 */
while (j--!=0) { *to++=(i>>7); i<<=1; }
}
}
} else if (oldBpc==2) {
// assert(0);
htc=ht; while (htc--!=0) {
toend=to+((wdcpp)&~3);
while (to!=toend) {
i=*p++;
*to++=(i>>6);
*to++=((i>>4)&3);
*to++=((i>>2)&3);
*to++=( i&3);
}
if (0!=(j=(wdcpp)&3)) {
i=*p++;
// fprintf(stderr,"j=%d\n",j);
while (j--!=0) { *to++=(i>>6); i<<=2; }
}
}
assert((slen_t)((char*)to-rowbeg)==(slen_t)wd*cpp*ht);
} else if (oldBpc==4) {
htc=ht; while (htc--!=0) {
toend=to+((wdcpp)&~1);
while (to!=toend) {
i=*p++;
*to++=(i>>4);
*to++=( i&15);
}
if (0!=((wdcpp)&1)) *to++=(*p++>>4);
}
} else assert(0 && "invalid bpc");

delete [] const_cast<char*>(oldBeg);
}

unsigned char Image::Sampled::minRGBBpc() const {
unsigned char *crow=new unsigned char[wd*3], *p, *pend=crow+wd*3;
register unsigned minbpb=0;
dimen_t htc;
for (htc=0;htc<ht;htc++) {
copyRGBRow((char*)crow, htc);
for (p=crow; p!=pend; p++) {
if ((*p&15)*17!=*p) { delete [] crow; return 8; } /* 4 bits are not enough */
else if ((*p&3)*85!=*p) minbpb=3; /* 2 bits are not enough */
else if ((*p&1)*255!=*p) minbpb|=1; /* 1 bit is not enough */
}
}
delete [] crow;
return 1+minbpb;
}

bool Image::Sampled::hasPixelRGB(Image::Sampled::rgb_t rgb) const {
/* by [email protected] at Sat Jan 8 13:24:19 CET 2005 */
/* Dat: this dumb implementation will be overridden */
if (rgb>0xffffffUL) return false;
unsigned char *crow=new unsigned char[wd*3], *p, *pend=crow+wd*3, t[3];
dimen_t htc;
t[0]=(rgb>>16)&255; t[1]=(rgb>>8)&255; t[2]=rgb&255;
for (htc=0;htc<ht;htc++) {
copyRGBRow((char*)crow, htc);
for (p=crow; p!=pend; p+=3) {
if (t[0]==p[0] && t[1]==p[1] && t[2]==p[2]) { delete [] crow; return true; }
}
}
delete [] crow;
return false;
}

bool Image::Gray::hasPixelRGB(Image::Sampled::rgb_t rgb) const {
/* by [email protected] at Sat Jan 8 13:24:19 CET 2005 */
/* Dat: faster than Image::Sampled::hasPixelRGB */
if (rgb>0xffffffUL) return false;
unsigned char t[3];
t[0]=(rgb>>16)&255; t[1]=(rgb>>8)&255; t[2]=rgb&255;
if (t[0]!=t[1] || t[0]!=t[2]) return false;
if (bpc==8) {
unsigned char *p=(unsigned char*)rowbeg, *pend=p+wd*ht;
/* Imp: use memchr() if available */
while (p!=pend && t[0]!=p[0]) p++;
return p!=pend;
}
unsigned char *crow=new unsigned char[wd*3], *p, *pend=crow+wd*3;
dimen_t htc;
for (htc=0;htc<ht;htc++) {
copyRGBRow((char*)crow, htc); /* Imp: avoid this if bpp==8 */
p=crow; while (p!=pend && t[0]!=p[0]) p+=3;
if (p!=pend) { delete [] crow; return true; }
}
delete [] crow;
return false;
}

Image::Indexed* Image::Sampled::addAlpha0(Image::Indexed *iimg, Image::Gray *al) {
if (iimg==NULLP) Error::sev(Error::EERROR) << "addAlpha: too many colors, transparency impossible" << (Error*)0;
iimg->to8();
unsigned ncols;
if ((ncols=iimg->getNcols()) == 256) {
iimg->packPal(); // TODO: Do a more lightweight palette packing if there are only 2 colors (such as PNG import through PNM).
if ((ncols=iimg->getNcols())==256) Error::sev(Error::EERROR) << "addAlpha: too many colors, transparency impossible" << (Error*)0;
}
iimg->setNcolsMove(ncols+1);
/* fprintf(stderr,"old ncols=%u\n", ncols); */
iimg->setPal(ncols,0); /* black */
iimg->setTransp(ncols);
assert(iimg->getRlen()==iimg->getWd());
assert(iimg->getWd()==al->getWd());
char *p=iimg->getRowbeg(), *pend=p+iimg->getRlen()*iimg->getHt(), *alq=al->getRowbeg();
while (p!=pend) {
if ((unsigned char)*alq++!=255) *p=ncols; /* make it transparent */
p++;
}
return iimg;
}

/* --- */

Image::Indexed::Indexed(Image::Sampled::dimen_t wd_, Image::Sampled::dimen_t ht_, unsigned short ncols_, unsigned char bpc_) {
param_assert(ncols_<=256);
/* vvv Dat: `3' is here for an extra palette entry */
init(3,3*ncols_,wd_,ht_,bpc_,TY_INDEXED,1);
transp=-1;
cs=CS_Indexed_RGB;
}
void Image::Indexed::setNcols(unsigned short ncols_) {
headp=rowbeg-ncols_*3;
}
void Image::Indexed::setNcolsMove(unsigned short ncols_) {
param_assert(ncols_<=256);
unsigned ncols=getNcols();
if (ncols_==ncols) return;
if (ncols_<ncols || (slen_t)(headp-beg)>=(ncols_-ncols)*3) {
memmove(rowbeg-ncols_*3, headp, (ncols_<ncols ? ncols_ : ncols)*3);
/* ^^^ *3 BUGFIX at Sun Apr 14 00:50:34 CEST 2002 */
} else { /* Imp: test this routine */
/* Tue Jun 11 16:22:52 CEST 2002 */
assert(ncols_>ncols);
const char *oldBeg=beg, *oldHeadp=headp, *oldRowbeg=rowbeg, *oldEnd=beg+len;
slen_t delta=(ncols_-ncols)*3;
// substr_grow(headp-oldBeg, ncols*3, ncols_*3); /* no such method */
beg=new char[len+delta];
headp= const_cast<char*>(beg)+(headp-oldBeg);
rowbeg=const_cast<char*>(beg)+(rowbeg-oldBeg)+delta;
trail= const_cast<char*>(beg)+(trail-oldBeg)+delta;
assert(beg+(headp-oldBeg)==rowbeg-ncols_*3);
/* Dat: this->xoffs is left unchanged */
memcpy(headp, oldHeadp, oldRowbeg-oldHeadp);
memcpy(rowbeg, oldRowbeg, oldEnd-oldRowbeg);
delete [] const_cast<char*>(oldBeg);
}
headp=rowbeg-ncols_*3;
}
void Image::Indexed::setPal(unsigned char color, Image::Sampled::rgb_t rgb) {
assert(color<(rowbeg-headp)/3);
unsigned char *p=(unsigned char*)headp+3*color;
*p++=rgb>>16;
*p++=rgb>>8;
*p=rgb;
}
void Image::Indexed::setTransp(unsigned char color) {
// param_assert(color>=0); /* always */
assert(transp==-1);
transp=color;
unsigned char *p=(unsigned char*)headp+3*color;
transpc=((Image::Sampled::rgb_t)p[0]<<16)+(p[1]<<8)+p[2];
}

bool Image::Indexed::setTranspc(rgb_t color) {
if (color!=0x1000000UL && color!=transpc) {
char t[3];
t[0]=color>>16; t[1]=color>>8; t[2]=color;
char *p=headp, *pend=rowbeg;
while (p!=pend) { /* Examine the palette. */
if (p[0]==t[0] && p[1]==t[1] && p[2]==t[2]) {
transpc=color;
transp=(p-headp)/3; /* destroy old transparency */
}
p+=3;
}
}
return transp!=-1;
}

bool Image::Indexed::wouldSetTranspc(rgb_t color) const {
if (transp!=-1) return true;
if (color!=0x1000000UL && color!=transpc) {
char t[3];
t[0]=color>>16; t[1]=color>>8; t[2]=color;
char *p=headp, *pend=rowbeg;
while (p!=pend) { /* Examine the palette. */
if (p[0]==t[0] && p[1]==t[1] && p[2]==t[2]) return true;
p+=3;
}
}
return false;
}

void Image::Indexed::setTranspcAndRepack(rgb_t color) {
if (!(color!=0x1000000UL && color!=transpc)) return;
char t[3];
t[0]=color>>16; t[1]=color>>8; t[2]=color;
char *p=headp, *pend=rowbeg;
bool need_repack = false;
while (p!=pend) { /* Examine the palette. */
if (p[0]==t[0] && p[1]==t[1] && p[2]==t[2]) {
transpc=color;
transp=(p-headp)/3; /* destroy old transparency */
need_repack = true;
}
p+=3;
}
if (need_repack) {
const unsigned char old_bpc = bpc;
packPal(); /* May change bpc. */
setBpc(old_bpc);
}
}

void Image::Indexed::to8() { to8nomul(); }
Image::Indexed* Image::Indexed::toIndexed()/* const*/ { return this; }
Image::RGB* Image::Indexed::toRGB(unsigned char bpc_)/* const*/ { return toRGB0(bpc_); }
Image::Gray* Image::Indexed::toGray(unsigned char bpc_)/* const*/ { return toGray0(bpc_); }
bool Image::Indexed::canGray() const {
char *p=headp, *pend=rowbeg, *tp=p+transp*3;
/* ignore transparent color at Sat Jun 15 15:18:24 CEST 2002 */
if (transp!=-1 && tp!=pend-3) {
while (p!=pend) { /* Examine the palette. */
if (p!=tp && (p[0]!=p[1] || p[1]!=p[2])) return false; /* Found a non-gray color. */
p+=3;
}
} else {
if (transp!=-1 && tp==pend-3) pend-=3; /* both conditions are important */
while (p!=pend) { /* Examine the palette. */
if (p[0]!=p[1] || p[1]!=p[2]) return false; /* Found a non-gray color. */
p+=3;
}
}
return true;
}
unsigned char Image::Indexed::minRGBBpc() const {
unsigned char *p=(unsigned char*)headp, *pend=(unsigned char*)rowbeg;
unsigned char *tp=p+transp*3;
/* ignore transparent color at Sat Jun 15 15:18:24 CEST 2002 */
register unsigned minbpb=0;
while (p!=pend) { /* Examine the palette. */
if (p==tp) { p+=3; continue; } /* ignore transparent color */
if ((*p&15)*17!=*p) return 8; /* 4 bits are not enough */
else if ((*p&3)*85!=*p) minbpb=3; /* 2 bits are not enough */
else if ((*p&1)*255!=*p) minbpb|=1; /* 1 bit is not enough */
p++;
}
return 1+minbpb;
}
void Image::Indexed::copyRGBRow(char *to, Image::Sampled::dimen_t whichrow) const {
param_assert(whichrow<ht);
if (wd==0) return;
unsigned char *p=(unsigned char*)rowbeg+rlen*whichrow;
char *r, *toend=to+3*wd;
unsigned int i, j;

if (bpc==1) {
toend-=3*(wd&7);
while (to!=toend) {
i=*p++;
r=headp+3*(i>>7); *to++=*r++; *to++=*r++; *to++=*r++;
r=headp+3*((i>>6)&1); *to++=*r++; *to++=*r++; *to++=*r++;
r=headp+3*((i>>5)&1); *to++=*r++; *to++=*r++; *to++=*r++;
r=headp+3*((i>>4)&1); *to++=*r++; *to++=*r++; *to++=*r++;
r=headp+3*((i>>3)&1); *to++=*r++; *to++=*r++; *to++=*r++;
r=headp+3*((i>>2)&1); *to++=*r++; *to++=*r++; *to++=*r++;
r=headp+3*((i>>1)&1); *to++=*r++; *to++=*r++; *to++=*r++;
r=headp+3*( i&1); *to++=*r++; *to++=*r++; *to++=*r++;
}
i=*p; /* No mem overrun, even if (wd&7)==0 */
j=wd&7;
while (j--!=0) { r=headp+3*(i>>7); *to++=*r++; *to++=*r++; *to++=*r++; i<<=1; }
} else if (bpc==2) {
toend-=3*(wd&3);
while (to!=toend) {
i=*p++;
r=headp+3*(i>>6); *to++=*r++; *to++=*r++; *to++=*r++;
r=headp+3*((i>>4)&3); *to++=*r++; *to++=*r++; *to++=*r++;
r=headp+3*((i>>2)&3); *to++=*r++; *to++=*r++; *to++=*r++;
r=headp+3*( i&3); *to++=*r++; *to++=*r++; *to++=*r++;
}
i=*p; /* No mem overrun, even if (wd&7)==0 */
j=wd&3;
while (j--!=0) { r=headp+3*(i>>6); *to++=*r++; *to++=*r++; *to++=*r++; i<<=2; }
} else if (bpc==4) {
toend-=3*(wd&1);
while (to!=toend) {
i=*p++;
r=headp+3*(i>>4); *to++=*r++; *to++=*r++; *to++=*r++;
r=headp+3*( i&15); *to++=*r++; *to++=*r++; *to++=*r++;
}
if (0!=(wd&1)) { r=headp+3*(*p>>4); *to++=*r++; *to++=*r++; *to++=*r++; }
} else if (bpc==8) {
// fprintf(stderr, "p=%u pp=%u ppp=%u\n", p[0], p[1], p[2]);
while (to!=toend) {
r=headp+3**p++; *to++=*r++; *to++=*r++; *to++=*r++;
}
} else assert(0 && "invalid bpc");
}
void Image::Indexed::packPal() {
/* Convert samples, make bpc=8. */
to8();

unsigned oldNcols=getNcols();
unsigned char *p, *pend;
assert((rowbeg-headp)%3==0);
assert(transp>=-1);
assert(transp<(int)oldNcols);
if (oldNcols<=1) return; /* Cannot optimize further. */

/* Find unused colors. old2new[c]=(is c used at least once)?1:0 */
unsigned char old2new[256], newpal[768];
memset(old2new, 0, sizeof(old2new));
for (p=(unsigned char*)rowbeg, pend=p+wd*ht; p!=pend; p++) old2new[*p]=1;

/* Find and eliminate duplicate colors. Build the new palette in the
* beginning of newpal. Use a Hash46 for a quick lookup of colors already
* seen. Use the previously computed old2new, but also overwrite it.
*/
Hash46 h;
int newTransp=-1;
unsigned char *op=old2new, *opend=op+oldNcols, *w, k[6],
*ptransp=(unsigned char*)headp+3*transp; /* ==p-3 if no transparent color */
/* ^^^ headp BUGFIX at Fri Mar 22 18:02:18 CET 2002 */
p=(unsigned char*)headp;
// fprintf(stderr, "oldNcols=%d\n", (int)oldNcols);
unsigned newNcols=0;
while (op!=opend) {
// fprintf(stderr, "color=%d %d\n", (int)(op-old2new), p-ptransp);
if (0!=*op) { /* Map the color only if it is used in the image. */
// fprintf(stderr, "used=%d\n", (int)(op-old2new));
if (p==ptransp) { k[0]=1; k[1]=k[2]=k[3]=0; newTransp=newNcols; }
else { k[0]=0; memcpy(k+1,p,3); }
w=h.walk(k);
assert(w!=NULL); /* Hash cannot be full since h.M>=256. */
if (*w==/*h.*/FREE) {
memcpy(newpal+3*newNcols, p /* k+1 */, 3);
/* ^^^ side effect: make the transparent color black */
memcpy(w,k,4); w[4]=newNcols; *op=newNcols++;
} else *op=w[4];
}
p+=3; op++;
}
// fprintf(stderr,"newTransp=%d transp=%d\n", newTransp, transp);

// assert((newTransp==-1) == (transp==-1));
assert(newTransp==-1 || transp!=-1);
/* ^^^ BUGFIX: == not true, because image may have transparency, but no
* transparent pixels.
*/
assert((char*)p==headp+oldNcols*3);
if (newNcols==oldNcols && transp==newTransp) {
/* Could not change # colors. */
if (transp==-1) goto done;
if ((unsigned)transp==oldNcols-1) { setPal(transp, 0); goto done; }
}

/* Make the transparent color last. */
if (newTransp!=-1 && newTransp!=(int)newNcols-1) {
assert(transp!=-1);
unsigned newLast=newNcols-1;
memcpy(newpal+3*newTransp, newpal+3*newLast, 3);
memset(newpal+3*newLast, 0, 3); transpc=0; /* make it black */
for (op=old2new; op!=opend; op++) if (*op==newLast) *op=newTransp;
old2new[transp]=newLast;
transp=newTransp=newLast;
p=newpal+newTransp*3;
}

/* Update the image. */
for (p=(unsigned char*)rowbeg, pend=p+wd*ht; p!=pend; p++) {
assert(*p<oldNcols);
*p=old2new[*p];
}

/* Update the palette. */
headp=rowbeg-3*newNcols; /* public method getNcols() uses rowbeg-headp */
memcpy(headp, newpal, 3*newNcols);
transp=newTransp;
/* vvv BUGFIX at Tue May 21 13:10:30 CEST 2002 */
if (newTransp==-1) transpc=0x1000000UL; /* Dat: this means: no transparent color */
else { transp=-1; setTransp(newTransp); }

done:
sortPal();
}
void Image::Indexed::sortPal() {
/* Run packPal() first (but it's recursive!) if transp is in the middle. */
unsigned ncols = getNcols(), i;
assert(transp == -1 || transp + 0U == ncols - 1);
assert(ncols <= 256);
if (ncols == 0) return; /* Safe if ncols == 0 and transp == -1. */
if (transp + 0U == ncols - 1) --ncols;
if (ncols <= 1) return;
#if SIZEOF_SHORT>=4
typedef unsigned short d_t;
#elif SIZEOF_INT>=4
typedef unsigned d_t;
#else
typedef unsigned long d_t;
#endif
d_t d[256];
unsigned char *p, *pend;
for (i = 0, p = (unsigned char*)headp; i < ncols; ++i, p += 3) {
d[i] = (d_t)p[0] << 24 | (d_t)p[1] << 16 | (d_t)p[2] << 8 | i;
/*printf("c[%d]=0x%08x\n", i, d[i]);*/
}
for (i = 1; i < ncols; ++i) {
if (d[i] < d[i - 1]) break;
}
if (i >= ncols) return; /* Palette already sorted. */

/* Heap sort (unstable). Based on Knuth's TAOCP 5.2.3.H .
* Although heap sort is unstable, sortPal implements a stable sort, because
* the color index (i) is included in the sorted number (d[i]).
*/
{ unsigned k;
d_t tmp, *i, *j, *l=d+(ncols>>1), *r=d+ncols-1;
while (1) { /* h2: */
k=l-d;
if (k!=0) {
tmp=*--l;
} else {
tmp=*r; *r=d[0];
if (--r==d) { d[0]=tmp; break; }
k++;
}
i=j=l;
while ((j+=k)<=r) { /* h4: */
k<<=1;
if (j<r && j[0]<j[1]) {
if (!(tmp<*++j)) break;
k++;
} else if (!(tmp<*j)) break;
*i=*j; i=j;
}
/* h8: */
*i=tmp;
}
}

unsigned char old2new[256];
for (i = 0, p = (unsigned char*)headp; i < ncols; ++i) {
d_t di = d[i];
/*printf("d[%d]=0x%08x\n", i, di);*/
assert((i == 0 || di >= d[i - 1]) && "bug in sorting palette");
old2new[di & 255] = i;
*p++ = di >> 24; *p++ = di >> 16; *p++ = di >> 8;
}

/* Update the image. */
for (p=(unsigned char*)rowbeg, pend=p+wd*ht; p!=pend; p++) {
*p=old2new[*p];
}
}

void Image::Indexed::delete_separated(register Indexed **p) {
while (*p!=NULLP) delete *p++;
}
Image::Indexed **Image::Indexed::separate() {
assert(getNcols()>=1);
unsigned char ncols1=getNcols()-1;
signed nncols=getNcols()-(transp==-1 ? 0 : 1);
register unsigned char curcol;
Indexed **ret=new Indexed*[nncols+1], **curimg=ret;
Image::Sampled::dimen_t htc;
assert(cpp==1);
slen_t wdcpp=wd/* *cpp*/;
register unsigned char *p;
char *to, *toend;
register unsigned int i;

ret[nncols]=(Indexed*)NULLP;
to8();
for (curcol=0; curcol<=ncols1; curcol++) {
if (transp==(signed int)curcol) continue;
curimg[0]=new Indexed(wd, ht, /*ncols:*/2, /*bpc:*/1);
memcpy(curimg[0]->headp, headp+3*curcol, 3); /* copy the color value */
curimg[0]->setTransp(1);
to=curimg[0]->rowbeg; p=(unsigned char*)rowbeg;
htc=ht; while (htc--!=0) {
toend=to+((wdcpp+7)>>3);
while (to!=toend) {
i =(*p++!=curcol)<<7; i|=(*p++!=curcol)<<6;
i|=(*p++!=curcol)<<5; i|=(*p++!=curcol)<<4;
i|=(*p++!=curcol)<<3; i|=(*p++!=curcol)<<2;
i|=(*p++!=curcol)<<1; i|=(*p++!=curcol);
*to++=i;
}
if (0!=(wdcpp&7)) p+=(wdcpp&7)-8; /* negative */
}
curimg++;
}
assert(curimg==ret+nncols);
return ret;
}

Image::Indexed *Image::Indexed::calcAlpha() {
/* by [email protected] at Tue Jun 4 21:27:29 CEST 2002 */
assert(getNcols()>=1);
packPal(); /* removes transparency if no transparent pixel */
if (transp==-1) return (Image::Indexed*)NULLP;
to8();
Indexed *ret=new Indexed(wd, ht, /*ncols:*/2, /*bpc:*/1);
Image::Sampled::dimen_t htc;
assert(cpp==1);
slen_t wdcpp=wd/* *cpp*/;
register unsigned char *p;
char *to, *toend;
register unsigned int i, i8, i7;
unsigned char transpx=transp;

ret->headp[0]=ret->headp[1]=ret->headp[2]='\xFF'; /* white */
ret->headp[3]=ret->headp[4]=ret->headp[5]='\x00'; /* black, transparent */
ret->setTransp(1);
to=ret->rowbeg; p=(unsigned char*)rowbeg;
assert(transpx!=0);
#if 0
printf("tx=%u\n", transpx);
printf("%u %u %u\n", headp[0], headp[1], headp[2]);
#endif
htc=ht; while (htc--!=0) {
// putchar('.'); printf("mod=%d\n",(to-ret->rowbeg)%ret->rlen);
// assert((to-ret->rowbeg)%ret->rlen==0);
toend=to+(wdcpp>>3); /* ((wdcpp+7)>>3)-1; */
/* ^^^ BUGFIX at Tue Sep 17 11:08:46 CEST 2002 */
assert(toend>=to);
while (to!=toend) {
#if 1 /* add ->pal[0] funcitonality at Sat Jun 15 14:24:25 CEST 2002 */
i=0; i8=256;
/* vvv p[-1]=0 BUGFIX at Sun Dec 8 23:21:47 CET 2002 */
while ((i8>>=1)!=0) if (*p++==transpx) { p[-1]=0; i|=i8; }
#else
i =(*p++==transpx)<<7; i|=(*p++==transpx)<<6;
i|=(*p++==transpx)<<5; i|=(*p++==transpx)<<4;
i|=(*p++==transpx)<<3; i|=(*p++==transpx)<<2;
i|=(*p++==transpx)<<1; i|=(*p++==transpx);
#endif
*to++=i;
}
#if 1 /* This works even when p gets modified; this puts fixed 0 pads at EOLs */
if ((wdcpp&7)!=0) {
i7=1<<(7-(wdcpp&7)); i8=256; i=0;
/* vvv p[-1]=0 BUGFIX at Sun Dec 8 23:21:47 CET 2002 */
while ((i8>>=1)!=i7) if (*p++==transpx) { p[-1]=0; i|=i8; }
*to++=i;
}
#else
if (0!=(wdcpp&7)) p+=(wdcpp&7)-8; /* negative */
#endif
}
assert(ret->rlen==((wd+7)>>3));
/* printf("rlen=%d %d\n", ret->rlen, to-ret->rowbeg); */
assert(to==ret->rowbeg+ht*ret->rlen);
return ret;
}

void Image::Indexed::setBpc(unsigned char bpc_) {
unsigned ncols=getNcols();
if (bpc_==0) {
if (ncols<=2) bpc_=1;
else if (ncols<=4) bpc_=2;
else if (ncols<=16) bpc_=4;
else bpc_=8;
} else {
if (bpc_==1) assert(ncols<=2);
else if (bpc_==2) assert(ncols<=4);
else if (bpc_==4) assert(ncols<=16);
else if (bpc_!=8) param_assert(0 && "invalid bpc_");
}
// fprintf(stderr, "bpc: %u -> %u\n", bpc, bpc_);
if (bpc==bpc_) return;
to8(); /* Imp: make the transition without the intermediate 8-bits... */
if (bpc_==8) return;
if (ht==0 || wd==0) { bpc=bpc_; return; }

const char *oldBeg=beg;
unsigned char *p=(unsigned char*)rowbeg;
assert(cpp==1);
slen_t wdcpp=wd/* *cpp*/;
bpc=bpc_;
rlen=(((rlen_t)bpc_)*wd+7)>>3;
beg=new char[len=rowbeg-oldBeg+rlen*ht+bpc];
headp= const_cast<char*>(beg)+(headp-oldBeg);
rowbeg=const_cast<char*>(beg)+(rowbeg-oldBeg);
trail= const_cast<char*>(beg)+len-bpc;
memcpy(const_cast<char*>(beg), oldBeg, rowbeg-beg);

unsigned char *to=(unsigned char*)rowbeg, *toend;
unsigned int i;
Image::Sampled::dimen_t htc;
if (bpc_==1) {
// This reads bytes from trail.
htc=ht; while (htc--!=0) {
toend=to+((wdcpp+7)>>3);
while (to!=toend) {
i =*p++<<7; i|=*p++<<6; i|=*p++<<5; i|=*p++<<4;
i|=*p++<<3; i|=*p++<<2; i|=*p++<<1; i|=*p++;
*to++=i;
}
if (0!=(wdcpp&7)) p+=(wdcpp&7)-8; /* negative */
}
} else if (bpc_==2) {
// This reads bytes from trail.
htc=ht; while (htc--!=0) {
toend=to+((wdcpp+3)>>2);
while (to!=toend) {
i =*p++<<6; i|=*p++<<4; i|=*p++<<2; i|=*p++;
*to++=i;
}
if (0!=(wdcpp&3)) p+=(wdcpp&3)-4;
}
} else if (bpc_==4) {
// This reads bytes from trail.
htc=ht; while (htc--!=0) {
toend=to+((wdcpp+1)>>1);
while (to!=toend) {
i =*p++<<4; i|=*p++;
*to++=i;
}
if (0!=(wdcpp&1)) p--;
}
} else assert(0 && "invalid bpc");
delete [] const_cast<char*>(oldBeg);
}
Image::Sampled::rgb_t Image::Indexed::getPal(unsigned char color) const {
unsigned char *p=(unsigned char*)headp+3*color;
return ((Image::Sampled::rgb_t)p[0]<<16)+(p[1]<<8)+p[2];
}

void Image::Indexed::dumpDebug(GenBuffer::Writable& gw) {
gw <<"% ncols=" << getNcols()
<< " rlen=" << rlen
<< " ht=" << ht
<< " wd=" << wd
<< " cpp=" << (unsigned)cpp
<< " bpc=" << (unsigned)bpc
<< " transp=" << transp
<< " transpc=" << (transp==-1?"none":rgb2webhash(getPal(transp)))
<< " ty=" << (unsigned)ty
//<< " pred=" << (unsigned)pred
<< '\n';
unsigned char *p=(unsigned char*)headp;
while (p!=(unsigned char*)rowbeg) {
gw << rgb2webhash(((Image::Sampled::rgb_t)p[0]<<16)+(p[1]<<8)+p[2]) << '\n';
p+=3;
}
gw << '\n';
gw.vi_write(rowbeg,rlen*ht);
}

Image::Sampled* Image::Indexed::addAlpha(Image::Gray *al) {
// Error::sev(Error::WARNING) << "Indexed: alpha channel ignored" << (Error*)0; return this;
if (al->getHt()!=ht || al->getWd()!=wd) Error::sev(Error::EERROR) << "addAlpha: image dimension mismatch" << (Error*)0;
bool ign_mid=false;
unsigned char lightest, darkest;
al->to8();
al->calcExtrema(lightest, darkest);
if (darkest==255) return this; /* no transparent pixels at all */
to8();
if (transp>=0) { /* Already have a transparent index. Join. */
register char *p=rowbeg;
register unsigned char transp_=transp;
char *pend=rowbeg+rlen*ht, *alq=al->getRowbeg();
/* Imp: choose an image color instead of black... */
/* Dat: 0..254: transparent, 255: opaque */
while (p!=pend) {
if ((unsigned char)(*alq+1)>1) ign_mid=true;
/* fprintf(stderr,"alq=%u\n", (unsigned char)*alq); */
if ((unsigned char)*alq++!=255) p[0]=transp_; /* black out transparent-wannabe pixels */
p++;
}
} else { /* No transparent color yet. */
packPal();
unsigned ncols=getNcols();
char *p=rowbeg, *pend=rowbeg+rlen*ht, *alq=al->getRowbeg();
/* Imp: choose an image color instead of black... */
/* Dat: 0..254: transparent, 255: opaque */
while (p!=pend) {
if ((unsigned char)(*alq+1)>1) ign_mid=true;
/* fprintf(stderr,"alq=%u\n", (unsigned char)*alq); */
if ((unsigned char)*alq++!=255) p[0]=ncols; /* may set to 0 if ncols==256 */
p++;
}
if (ncols==256) { /* Try again, probably now we have fewer colors */
packPal();
if ((ncols=getNcols())==256) Error::sev(Error::EERROR) << "Indexed::addAlpha: too many colors, transparency impossible" << (Error*)0;
for (p=rowbeg,alq=al->getRowbeg(); p!=pend; p++)
if ((unsigned char)*alq++!=255) *p=ncols;
}
setNcolsMove(ncols+1);
setPal(ncols,0); /* black */
setTransp(ncols);
}
if (ign_mid) Error::sev(Error::WARNING) << "addAlpha: half-transparent pixels made transparent" << (Error*)0;
return this;
}
void Image::Indexed::makeTranspZero() {
if (transp<1) return; /* no transparency or already 0 */
unsigned char oldBpc=bpc;
register unsigned char transpch=transp;
/* Imp: make this faster by not converting to 8 bits */
if (oldBpc!=8) to8();

/* Update the image. */
register unsigned char *p; unsigned char *pend;
for (p=(unsigned char*)rowbeg, pend=p+wd*ht; p!=pend; p++) {
if (*p==0) *p=transp;
else if (*p==transpch) *p=0;
}
rgb_t rzero=getPal(0), rtransp=getPal(transp);
setPal(transp, rzero); setPal(0, rtransp);
transp=0;
if (oldBpc!=8) setBpc(oldBpc);
}

/* --- */

Image::Gray::Gray(Image::Sampled::dimen_t wd_, Image::Sampled::dimen_t ht_, unsigned char bpc_) {
init(0,0,wd_,ht_,bpc_,TY_GRAY,1);
cs=CS_GRAYSCALE;
}
void Image::Gray::to8() { to8mul(); }
Image::RGB * Image::Gray::toRGB(unsigned char bpc_)/* const*/ { return toRGB0(bpc_); }
Image::Gray* Image::Gray::toGray(unsigned char bpc_)/* const*/ { return bpc==bpc_ ? this : toGray0(bpc_); }
Image::Indexed* Image::Gray::toIndexed()/* const*/ {
Image::Indexed *img=new Image::Indexed(wd, ht, (1<<bpc), bpc);
unsigned i;
rgb_t rgb;
if (bpc==1) {
img->setPal(0,0);
img->setPal(1,0xffffffL);
} else if (bpc==2) {
img->setPal(0,0);
img->setPal(1,0x555555L);
img->setPal(2,0xaaaaaaL);
img->setPal(3,0xffffffL);
} else if (bpc==4) {
for (i=0,rgb=0;i<16;i++,rgb+=(rgb_t)0x111111L) img->setPal(i,rgb);
} else if (bpc==8) {
for (i=0,rgb=0;i<256;i++,rgb+=(rgb_t)0x010101L) img->setPal(i,rgb);
}
memcpy(img->getRowbeg(), rowbeg, rlen*ht);
return img;
}
bool Image::Gray::canGray() const { return true; }
void Image::Gray::copyRGBRow(char *to, Image::Sampled::dimen_t whichrow) const {
param_assert(whichrow<ht);
if (wd==0) return;
unsigned char *p=(unsigned char*)rowbeg+rlen*whichrow;
char *toend=to+3*wd;
unsigned int i, j, k;

if (bpc==1) {
toend-=3*(wd&7);
while (to!=toend) {
i=*p++;
k=(i>>7)*255; *to++=k; *to++=k; *to++=k;
k=((i>>6)&1)*255; *to++=k; *to++=k; *to++=k;
k=((i>>5)&1)*255; *to++=k; *to++=k; *to++=k;
k=((i>>4)&1)*255; *to++=k; *to++=k; *to++=k;
k=((i>>3)&1)*255; *to++=k; *to++=k; *to++=k;
k=((i>>2)&1)*255; *to++=k; *to++=k; *to++=k;
k=((i>>1)&1)*255; *to++=k; *to++=k; *to++=k;
k=( i&1)*255; *to++=k; *to++=k; *to++=k;
}
i=*p; /* No mem overrun, even if (wd&7)==0 */
j=wd&7;
while (j--!=0) { k=(i>>7)*255; *to++=k; *to++=k; *to++=k; i<<=1; }
} else if (bpc==2) {
toend-=3*(wd&3);
while (to!=toend) {
i=*p++;
k=(i>>6)*85; *to++=k; *to++=k; *to++=k;
k=((i>>4)&3)*85; *to++=k; *to++=k; *to++=k;
k=((i>>2)&3)*85; *to++=k; *to++=k; *to++=k;
k=( i&3)*85; *to++=k; *to++=k; *to++=k;
}
i=*p; /* No mem overrun, even if (wd&7)==0 */
j=wd&3;
while (j--!=0) { k=(i>>6)*85; *to++=k; *to++=k; *to++=k; i<<=2; }
} else if (bpc==4) {
toend-=3*(wd&1);
while (to!=toend) {
i=*p++;
k=(i>>4)*17; *to++=k; *to++=k; *to++=k;
k=( i&15)*17; *to++=k; *to++=k; *to++=k;
}
if (0!=(wd&1)) { k=(*p>>4)*17; *to++=k; *to++=k; *to++=k; }
} else if (bpc==8) {
while (to!=toend) { *to++=*p; *to++=*p; *to++=*p++; }
} else assert(0 && "invalid bpc");
}
#if 0
void Image::Gray::setBpc(unsigned char bpc_) {
(void)bpc_;
assert(0 && "unimplemented"); /* unimplemented */
}
#endif

Image::Sampled* Image::Gray::addAlpha(Image::Gray *al) {
// Error::sev(Error::WARNING) << "Gray: alpha channel ignored" << (Error*)0; return this;
if (al->getHt()!=ht || al->getWd()!=wd) Error::sev(Error::EERROR) << "addAlpha: image dimension mismatch" << (Error*)0;
bool ign_mid=false;
unsigned char lightest, darkest;
al->to8();
al->calcExtrema(lightest, darkest);
if (darkest==255) return this; /* no transparent pixels at all */
char *p=rowbeg, *pend=rowbeg+rlen*ht, *alq=al->getRowbeg();
/* Imp: choose an image color instead of black... */
/* Dat: 0..254: transparent, 255: opaque */
while (p!=pend) {
if ((unsigned char)(*alq+1)>1) ign_mid=true;
/* fprintf(stderr,"alq=%u\n", (unsigned char)*alq); */
if ((unsigned char)*alq++!=255) *p=0; /* black out transparent-wannabe pixels */
p++;
}
if (ign_mid) Error::sev(Error::WARNING) << "addAlpha: half-transparent pixels made transparent" << (Error*)0;
return addAlpha0(toIndexed(), al);
}

void Image::Gray::calcExtrema(unsigned char &lightest, unsigned char &darkest) {
to8();
unsigned l=0, d=255, val;
char *p=rowbeg, *pend=rowbeg+ht*wd;
while (p!=pend) {
val=*(unsigned char*)p++;
if (val>l) l=val;
if (val<d) d=val;
}
lightest=l; darkest=d;
}

/* --- */

Image::RGB::RGB(Image::Sampled::dimen_t wd_, Image::Sampled::dimen_t ht_, unsigned char bpc_) {
init(0,0,wd_,ht_,bpc_,TY_RGB,3);
cs=CS_RGB;
}
void Image::RGB::to8() { to8mul(); }
Image::Gray* Image::RGB::toGray(unsigned char bpc_)/* const*/ { return toGray0(bpc_); }
Image::RGB* Image::RGB::toRGB(unsigned char bpc_)/* const*/ { return bpc==bpc_ ? this : toRGB0(bpc_); }
Image::Indexed* Image::RGB::toIndexed()/* const*/ { return toIndexed0(); }
bool Image::RGB::canGray() const { return false; }
void Image::RGB::copyRGBRow(char *to, Image::Sampled::dimen_t whichrow) const {
param_assert(whichrow<ht);
if (wd==0) return;
unsigned char *p=(unsigned char*)rowbeg+rlen*whichrow;
char *toend=to+3*wd;
unsigned int i, j, k;

if (bpc==1) {
toend-=(wd*3)&7;
while (to!=toend) {
i=*p++;
k=(i>>7)*255; *to++=k;
k=((i>>6)&1)*255; *to++=k;
k=((i>>5)&1)*255; *to++=k;
k=((i>>4)&1)*255; *to++=k;
k=((i>>3)&1)*255; *to++=k;
k=((i>>2)&1)*255; *to++=k;
k=((i>>1)&1)*255; *to++=k;
k=( i&1)*255; *to++=k;
}
i=*p; /* No mem overrun, even if (wd&7)==0 */
j=(wd*3)&7;
while (j--!=0) { k=(i>>7)*255; *to++=k; i<<=1; }
} else if (bpc==2) {
toend-=(wd*3)&3;
while (to!=toend) {
i=*p++;
k=(i>>6)*85; *to++=k;
k=((i>>4)&3)*85; *to++=k;
k=((i>>2)&3)*85; *to++=k;
k=( i&3)*85; *to++=k;
}
i=*p; /* No mem overrun, even if (wd&7)==0 */
j=(wd*3)&3;
while (j--!=0) { k=(i>>6)*85; *to++=k; i<<=2; }
} else if (bpc==4) {
toend-=(wd*3)&7;
while (to!=toend) {
i=*p++;
k=(i>>4)*17; *to++=k;
k=( i&15)*17; *to++=k;
}
if (0!=((wd*3)&1)) *to++=(*p>>4)*17;
} else if (bpc==8) {
memcpy(to, p, 3*wd);
} else assert(0 && "invalid bpc");
}
#if 0
void Image::RGB::setBpc(unsigned char bpc_) {
(void)bpc_;
assert(0 && "unimplemented"); /* unimplemented */
}
#endif

Image::Sampled* Image::RGB::addAlpha(Image::Gray *al) {
if (al->getHt()!=ht || al->getWd()!=wd) Error::sev(Error::EERROR) << "addAlpha: image dimension mismatch" << (Error*)0;
bool ign_mid=false;
unsigned char lightest, darkest;
al->to8();
al->calcExtrema(lightest, darkest);
if (darkest==255) return this; /* no transparent pixels at all */
char *p=rowbeg, *pend=rowbeg+rlen*ht, *alq=al->getRowbeg();
/* Imp: choose an image color instead of black... */
/* Dat: 0..254: transparent, 255: opaque */
while (p!=pend) {
if ((unsigned char)(*alq+1)>1) ign_mid=true;
/* fprintf(stderr,"alq=%u\n", (unsigned char)*alq); */
if ((unsigned char)*alq++!=255) p[0]=p[1]=p[2]=0; /* black out transparent-wannabe pixels */
p+=3;
}
if (ign_mid) Error::sev(Error::WARNING) << "addAlpha: half-transparent pixels made transparent" << (Error*)0;
return addAlpha0(toIndexed(), al);
}

/* --- */

char *Image::Sampled::rgb2webhash(rgb_t rgb) {
static char tmp[8];
char *p=tmp;
*p='#';
*++p='0'+(rgb>>20); if (*p>'9') *p+='a'-'0'-10;
*++p='0'+((rgb>>16)&15); if (*p>'9') *p+='a'-'0'-10;
*++p='0'+((rgb>>12)&15); if (*p>'9') *p+='a'-'0'-10;
*++p='0'+((rgb>>8)&15); if (*p>'9') *p+='a'-'0'-10;
*++p='0'+((rgb>>4)&15); if (*p>'9') *p+='a'-'0'-10;
*++p='0'+((rgb )&15); if (*p>'9') *p+='a'-'0'-10;
*++p='\0';
return tmp;
}

GenBuffer::Writable& operator<<(GenBuffer::Writable& gw, Image::Sampled const& img) {
slen_t buflen=img.getWd()*3;
char *buf=new char[buflen];
Image::Sampled::dimen_t y, ht=img.getHt();
/* vvv in the xv program: image file must be >=30 bytes long to be treated as image */
gw << "P6\n###############\n" << img.getWd() << ' ' << ht;
if (img.getTranspc()>=0x1000000UL) gw << "\n#Opaque";
else gw << "\n#T" << img.rgb2webhash(img.getTranspc());
gw << "\n255\n";
for (y=0; y<ht; y++) {
img.copyRGBRow(buf, y);
gw.vi_write(buf, buflen);
}
delete [] buf;
return gw;
}

static Image::Loader *first_image_loader=(Image::Loader*)NULLP;

void Image::register0(Image::Loader *anew) {
param_assert(anew!=NULLP);
anew->next=first_image_loader;
first_image_loader=anew;
}

#if 0 /* removed by code refactoring */
Image::Sampled* Image::load(char const* format, filep_t f_, SimBuffer::Flat const& loadHints) {} /* Removed. */
#endif

// #include <unistd.h> /* sleep() */

#if 0
Rule::Sampled *Rule::load(char const* filename) {
static char buf[2*Applier::MAGIC_LEN];
FILE *f=fopen(filename, "rb");
unsigned got=0;
if (f==NULLP) Error::sev(Error::EERROR) << "Cannot open/read image file: " << FNQ(filename) << (Error*)0;
slen_t ret=fread(buf, 1, Applier::MAGIC_LEN, f);
/* vvv Imp: clarify error message: may be a read error */
if (ret==0) Error::sev(Error::EERROR) << "Zero-length image file: " << FNQ(filename) << (Error*)0;
if (ret<Applier::MAGIC_LEN) memset(buf+ret, '\0', Applier::MAGIC_LEN-ret);
#if 0
unsigned long pos=fseek(f, 0, SEEK_END);
pos=(pos<=Applier::MAGIC_LEN)?0:pos-Applier::MAGIC_LEN;
if (0!=fseek(f, pos, SEEK_SET)
|| (got=fread(buf+Applier::MAGIC_LEN, 1, Applier::MAGIC_LEN, f))==0
#else
if (0
#endif
|| (rewind(f), 0)
|| ferror(f))
Error::sev(Error::EERROR) << "I/O error in image file: " << FNQ(filename) << (Error*)0;
if (got!=0 && got!=Applier::MAGIC_LEN) memmove(buf+2*Applier::MAGIC_LEN-got, buf+Applier::MAGIC_LEN, got);
Applier *p=first_image_loader;
Applier::reader_t reader;
while (p!=NULLP) {
if (NULLP!=(reader=p->checker(buf,buf+Applier::MAGIC_LEN))) { return reader(f); }
p=p->next;
}
Error::sev(Error::EERROR) << "Unknown image format: " << FNQ(filename) << (Error*)0;
// Error::sev(Error::WARNING) << "Zero-length image1." << (Error*)0;
// Error::sev(Error::WARNING) << "Zero-length image2." << (Error*)0;
return 0; /*notreached*/
}
#endif

Image::Sampled *Image::load(Image::Loader::UFD* ufd0, SimBuffer::Flat const& loadHints, char const* format) {
Filter::UngetFILED &ufd=*(Filter::UngetFILED*)ufd0;
/* Dat: format arg used in in_pnm.cpp */
static char buf[Loader::MAGIC_LEN+1];
slen_t ret=ufd.vi_read(buf, Loader::MAGIC_LEN);
/* vvv Imp: clarify error message: may be a read error */
if (ufd.hadError()) Error::sev(Error::EERROR) << "I/O error pre in image file: " << FNQ(ufd.getFilenameDefault("-")) << (Error*)0;
if (ret==0) Error::sev(Error::EERROR) << "Zero-length image file: " << FNQ(ufd.getFilenameDefault("-")) << (Error*)0;
if (ret<Loader::MAGIC_LEN) memset(buf+ret, '\0', Loader::MAGIC_LEN-ret);
buf[Loader::MAGIC_LEN]='\0';
/* Dat: do not read the trailer onto buf+Loader::MAGIC_LEN, because no ->checker() uses it yet. */
Loader *p=first_image_loader;
Loader::reader_t reader;
ufd.unread(buf, ret); /* tries to seek back, on failure calls ufd.getUnget().vi_write() */
// ^^^ rewind(f); /* checker might have read */
/* ^^^ do this early for the checkers */
while (p!=NULLP) {
/* vvv each checker() must rewind ufd for itself */
if ((format==(char const*)NULLP || 0==strcmp(p->format, format))
&& (Loader::checker_t)0!=p->checker
&& (Loader::reader_t)0!=(reader=p->checker(buf,buf+Loader::MAGIC_LEN, loadHints, ufd0))
) {
// fprintf(stderr, "%p %p\n", ufd0, &ufd);
return reader(ufd0, loadHints);
}
p=p->next;
}
// sleep(1000);
Error::sev(Error::EERROR) << "Unknown input image format: " << FNQ(ufd.getFilenameDefault("-")) << (Error*)0;
return 0; /*notreached*/
}

#if 0 /* not used anywhere, except in test_main. */
Image::Sampled *Image::load(char const* filename, SimBuffer::Flat const& loadHints, filep_t stdin_f, char const* format) { /* Commented out. */
Filter::UngetFILED ufd(filename, stdin_f==NULLP ? stdin : (FILE*)stdin_f,
Filter::UngetFILED::CM_closep|Filter::UngetFILED::CM_keep_stdinp);
return load((Image::Loader::UFD*)&ufd, loadHints, format);
// Imp: better error message, something like: if (f==NULLP) Error::sev(Error::EERROR) << "Cannot open/read image file: " << FNQ(filename) << (Error*)0;
}
#endif

#if 0 /* before Sat Apr 19 13:42:04 CEST 2003 */
Image::Sampled *Image::load(char const* filename, SimBuffer::Flat const& loadHints, filep_t stdin_f, char const* format) {
/* Dat: format arg used in in_pnm.cpp */
static char buf[2*Loader::MAGIC_LEN];
bool stdin_p=filename[0]=='-' && filename[1]=='\0';
FILE *f=!stdin_p ? fopen(filename, "rb") : stdin_f!=NULLP ? (FILE*)stdin_f : stdin;
unsigned got=0;
if (f==NULLP) Error::sev(Error::EERROR) << "Cannot open/read image file: " << FNQ(filename) << (Error*)0;
slen_t ret=fread(buf, 1, Loader::MAGIC_LEN, f);
/* vvv Imp: clarify error message: may be a read error */
if (ret==0) Error::sev(Error::EERROR) << "Zero-length image file: " << FNQ(filename) << (Error*)0;
if (ret<Loader::MAGIC_LEN) memset(buf+ret, '\0', Loader::MAGIC_LEN-ret);
#if 0
/* Dat: do not read the trailer, because no ->checker() uses it yet. */
unsigned long pos=fseek(f, 0, SEEK_END);
pos=(pos<=Loader::MAGIC_LEN)?0:pos-Loader::MAGIC_LEN;
if (0!=fseek(f, pos, SEEK_SET)
|| (got=fread(buf+Loader::MAGIC_LEN, 1, Loader::MAGIC_LEN, f))==0
#else
if (0
#endif
|| (rewind(f), 0)
|| ferror(f))
Error::sev(Error::EERROR) << "I/O error pre in image file: " << FNQ(filename) << (Error*)0;
if (got!=0 && got!=Loader::MAGIC_LEN) memmove(buf+2*Loader::MAGIC_LEN-got, buf+Loader::MAGIC_LEN, got);
Loader *p=first_image_loader;
Loader::reader_t reader;
while (p!=NULLP) {
if ((format==(char const*)NULLP || 0==strcmp(p->format, format))
&& (Loader::checker_t)NULLP!=p->checker
&& (Loader::reader_t)NULLP!=(reader=p->checker(buf,buf+Loader::MAGIC_LEN, loadHints, f))
) {
rewind(f); /* checker might have read */
Image::Sampled *ret=reader(f, loadHints);
if (ferror(f) || (!stdin_p && 0!=fclose(f))) /* don't close stdin */
Error::sev(Error::EERROR) << "I/O error post in image file: " << FNQ(filename) << (Error*)0;
return ret;
}
p=p->next;
}
Error::sev(Error::EERROR) << "Unknown input image format: " << FNQ(filename) << (Error*)0;
// Error::sev(Error::WARNING) << "Zero-length image1." << (Error*)0;
// Error::sev(Error::WARNING) << "Zero-length image2." << (Error*)0;
return 0; /*notreached*/
}
#endif

unsigned Image::printLoaders(GenBuffer::Writable &out) {
unsigned num=0;
Loader *p=first_image_loader;
while (p!=NULLP) {
if (p->checker!=(Loader::checker_t)0 && p->format!=(char const*)NULLP) { num++; out << ' ' << p->format; }
p=p->next;
}
return num;
}

/* --- */

Image::SampledInfo::SampledInfo(Sampled *img_)
:hasTransp(false)
,nncols(257)
,canGray(false)
,minRGBBpc(8)
,img(img_)
,imgs((Indexed**)NULLP) {
param_assert(img_!=NULLP);
Sampled *bak=img;
if ((img=img->toIndexed())==NULLP) {
img=bak;
} else {
if (bak!=img) delete bak;
assert(img->getTy()==img->TY_INDEXED);
Indexed *iimg=PTS_dynamic_cast(Image::Indexed*,img);
/* This packPal() contains a call to sortPal(), which converts the indexed
* image to canonical form.
*/
iimg->packPal();
nncols=iimg->getNcols();
if (true==(hasTransp=iimg->hasTransp())) nncols--;
}
minRGBBpc=img->minRGBBpc();
if ((canGray=img->canGray())==true && nncols==257) nncols=256;
sf=(sf_t)((img->getTy()==img->TY_BLACKBOX) ? 0+SF_Asis : 0+SF_None);
/* Dat: 0+: pacify gcc-3.1 */
}

Image::SampledInfo::~SampledInfo() {
delete img;
if (imgs!=NULLP) { Image::Indexed::delete_separated(imgs); delete imgs; }
}

void Image::SampledInfo::separate() {
// bool ok;
// ASSERT_SIDE(ok=(sf!=SF_Transparent2 && sf!=SF_Transparent4 && sf!=SF_Transparent8));
// if (!ok) return;
if (sf!=SF_Transparent2 && sf!=SF_Transparent4 && sf!=SF_Transparent8) return;
Indexed *iimg=PTS_dynamic_cast(Indexed*,img);
imgs=iimg->separate();
}

bool Image::SampledInfo::setSampleFormat(sf_t sf_, bool WarningOK, bool TryOnly, Sampled::rgb_t Transparent) {
/* at Sat Jun 15 11:48:51 CEST 2002: added transparency warnings */
/* fprintf(stderr, "sf=%u sf_=%u transparent=0x%lx\n", sf, sf_, Transparent+0UL); */
// assert(sf_==SF_Asis);
Indexed *iimg;
Sampled *bak=img;
param_assert(sf_!=SF_None);
param_assert(sf==SF_None || sf==SF_Asis || sf==sf_);
if (sf_==sf) return true; /* already converted */
bool zero=img->getWd()==0 || img->getHt()==0;
if (sf==SF_Asis && sf_!=SF_Asis && sf_!=SF_Bbox) {
Error::sev(Error::WARNING) << "SampleFormat: can't convert image loaded as Asis to other" << (Error*)0;
return false;
}
switch (sf_) {
case SF_Bbox:
sf=SF_Bbox; return true;
case SF_Opaque:
if (!hasTransp && nncols==1 && PTS_dynamic_cast(Indexed*,img)->wouldSetTranspc(Transparent)) return false;
if (hasTransp || nncols!=1) return false;
assert(img->getTy()==img->TY_INDEXED);
/* The color can be calculated: PTS_dynamic_cast(Indexed*>(img)->getPal(0); */
/* Conversion is not necessary. */
sf=SF_Opaque; return true;
case SF_Transparent:
if (!hasTransp && nncols==1 && PTS_dynamic_cast(Indexed*,img)->setTranspc(Transparent)) {
/* Must return true eventually, because setTranspc has modified the image. */
hasTransp=true; nncols=0;
} else if (!hasTransp || nncols!=0) {
return false;
}
assert(img->getTy()==img->TY_INDEXED);
/* Conversion is not necessary. */
sf=SF_Transparent; return true;
case SF_Gray1:
/* vvv strict hasTransp added at Mon Sep 9 22:53:24 CEST 2002 */
if (nncols>2 || !canGray || minRGBBpc>1 || zero || hasTransp || img->hasPixelRGB(Transparent)) return false;
/* ^^^ Imp: !! make the hasPixelRFB() check a lot faster (cache results) */
if (TryOnly) return WarningOK || (nncols>=2 && !hasTransp);
if (hasTransp) {
if (!WarningOK) return false;
/* Dat: nncols may be 1 or 2 ! (both transparent and opaque black) */
Error::sev(Error::WARNING) <<
(nncols<=1 ? "SampleFormat: Mask would be better than "
: "SampleFormat: Transparent2 would be better than ")
<< "Gray1" << (Error*)0;
}
if (nncols<2) {
if (!WarningOK) return false;
Error::sev(Error::WARNING) << "SampleFormat: Opaque would be better than Gray1" << (Error*)0;
}
img=img->toGray(1); if (bak!=img) delete bak;
assert(img!=NULLP);
assert(img->getBpc()==1);
sf=SF_Gray1; return true;
case SF_Indexed1:
if (nncols>2 || zero || hasTransp || img->hasPixelRGB(Transparent)) return false;
if (TryOnly) return WarningOK || nncols>=2;
if (nncols<2) {
if (!WarningOK) return false;
Error::sev(Error::WARNING) << "SampleFormat: Opaque would be better than Indexed1" << (Error*)0;
}
if (canGray && minRGBBpc==1) {
Error::sev(Error::NOTICE) << "SampleFormat: Gray1 would be better than Indexed1" << (Error*)0;
}
assert(img->getTy()==img->TY_INDEXED);
// img=img->toIndexed(); /* should be a no-op */
// assert(img!=NULLP); if (bak!=img) delete bak;
{ iimg=PTS_dynamic_cast(Indexed*,img);
iimg->setBpc(1);
if (iimg->wouldSetTranspc(Transparent)) return false; /* Dat: false if must be changed to become transparent. */
iimg->setTranspcAndRepack(Transparent);
}
sf=SF_Indexed1; return true;
case SF_Mask:
if (!hasTransp && nncols==2 && PTS_dynamic_cast(Indexed*,img)->wouldSetTranspc(Transparent)) { hasTransp=true; --nncols; }
if (nncols>1 || zero) return false;
if (TryOnly) return WarningOK || nncols+(hasTransp?1:0)==2;
if (nncols==1 && !hasTransp) {
if (!WarningOK) return false;
Error::sev(Error::WARNING) << "SampleFormat: Opaque would be better than Mask" << (Error*)0;
}
if (nncols==0) {
assert(hasTransp);
Error::sev(Error::WARNING) << "SampleFormat: Transparent would be better than Mask" << (Error*)0;
}
assert(img->getTy()==img->TY_INDEXED);
{ iimg=PTS_dynamic_cast(Indexed*,img);
iimg->setBpc(1);
if (!iimg->wouldSetTranspc(Transparent)) return false;
iimg->setTranspcAndRepack(Transparent);
}
/* printf("gett=%d\n", PTS_dynamic_cast(Indexed*,img)->getTransp()); */
/* vvv BUGFIX: <1U -> <2U */
assert(PTS_dynamic_cast(Indexed*,img)->getTransp()==-1 || PTS_dynamic_cast(Indexed*,img)->getTransp()+0U<2U);
/* ^^^ color 0 is opaque, color 1 is transparent, thanks to
* img->packPal() called in SampleInfo() -- but setTranspc may have changed this
*/
sf=SF_Mask; return true;
case SF_Transparent2:
if (!hasTransp && nncols==4 && PTS_dynamic_cast(Indexed*,img)->wouldSetTranspc(Transparent)) { hasTransp=true; --nncols; }
if (nncols>3 || zero) return false;
if (TryOnly) return WarningOK || (hasTransp && nncols>=2);
Error::sev(Error::NOTICE) << "SampleFormat: Transparent2 separates colors" << (Error*)0;
if (!hasTransp) {
if (!WarningOK) return false;
Error::sev(Error::WARNING) << "SampleFormat: Indexed2 would be better than Transparent2" << (Error*)0;
}
if (nncols<2) {
if (!WarningOK) return false;
Error::sev(Error::WARNING) << "SampleFormat: Mask would be better than Transparent2" << (Error*)0;
}
assert(img->getTy()==img->TY_INDEXED);
{ iimg=PTS_dynamic_cast(Indexed*,img);
iimg->setBpc(2); /* BUGFIX at Sat Jun 15 13:55:25 CEST 2002 */
iimg->setTranspcAndRepack(Transparent);
// imgs=iimg->separate(); /* postponed because of GIF89a output */
}
sf=SF_Transparent2; return true;
case SF_Gray2:
if (nncols>4 || !canGray || minRGBBpc>2 || zero || hasTransp || img->hasPixelRGB(Transparent)) return false;
if (TryOnly) return WarningOK || (nncols>2 && !hasTransp);
if (hasTransp) {
if (!WarningOK) return false;
Error::sev(Error::WARNING) <<
(nncols<=3 ? "SampleFormat: Transparent2 would be better than "
: "SampleFormat: Transparent4 would be better than ")
<< "Gray2" << (Error*)0;
}
if (nncols<=2) {
if (!WarningOK) return false;
Error::sev(Error::WARNING) << "SampleFormat: Gray1 would be better than Gray2" << (Error*)0;
}
img=img->toGray(2); if (bak!=img) delete bak;
assert(img!=NULLP);
sf=SF_Gray2; return true;
case SF_Indexed2:
if (nncols>4 || zero || hasTransp || img->hasPixelRGB(Transparent)) return false;
if (TryOnly) return WarningOK || nncols>2;
if (nncols<=2) {
if (!WarningOK) return false;
Error::sev(Error::WARNING) << "SampleFormat: Indexed1 would be better than Indexed2" << (Error*)0;
}
if (canGray && minRGBBpc<=2) {
Error::sev(Error::NOTICE) << "SampleFormat: Gray2 would be better than Indexed2" << (Error*)0;
}
assert(img->getTy()==img->TY_INDEXED);
{ iimg=PTS_dynamic_cast(Indexed*,img);
iimg->setBpc(2);
if (iimg->wouldSetTranspc(Transparent)) return false;
iimg->setTranspcAndRepack(Transparent);
}
sf=SF_Indexed2; return true;
case SF_Transparent4:
if (!hasTransp && nncols==16 && PTS_dynamic_cast(Indexed*,img)->wouldSetTranspc(Transparent)) { hasTransp=true; --nncols; }
if (nncols>15 || zero) return false;
if (TryOnly) return WarningOK || (hasTransp && nncols>=4);
Error::sev(Error::NOTICE) << "SampleFormat: Transparent4 separates colors" << (Error*)0;
if (!hasTransp) {
if (!WarningOK) return false;
Error::sev(Error::WARNING) << "SampleFormat: Indexed4 would be better than Transparent4" << (Error*)0;
}
if (nncols<4) {
if (!WarningOK) return false;
Error::sev(Error::WARNING) << "SampleFormat: Transparent2 would be better than Transparent4" << (Error*)0;
}
assert(img->getTy()==img->TY_INDEXED);
{ iimg=PTS_dynamic_cast(Indexed*,img);
iimg->setBpc(4);
iimg->setTranspcAndRepack(Transparent);
// imgs=iimg->separate(); /* postponed because of GIF89a output */
}
sf=SF_Transparent4; return true;
case SF_Rgb1:
if (nncols>8 || minRGBBpc>1 || zero || hasTransp || img->hasPixelRGB(Transparent)) return false;
if (TryOnly) return WarningOK || (nncols>4 && !canGray && !hasTransp);
if (hasTransp) {
if (!WarningOK) return false;
Error::sev(Error::WARNING) << "SampleFormat: Transparent4 would be better than " << "Rgb1" << (Error*)0;
}
if (canGray) {
if (!WarningOK) return false;
Error::sev(Error::WARNING) << "SampleFormat: Gray1 would be better than Rgb1" << (Error*)0;
}
if (nncols<=4) {
if (!WarningOK) return false;
Error::sev(Error::WARNING) << "SampleFormat: "
<< (hasTransp ? "Transparent2 may" : "Indexed2 would")
<< " be better than Rgb1" << (Error*)0;
}
img=img->toRGB(1); if (bak!=img) delete bak;
assert(img!=NULLP);
sf=SF_Rgb1; return true;
case SF_Gray4:
if (nncols>16 || !canGray || minRGBBpc>4 || zero || hasTransp || img->hasPixelRGB(Transparent)) return false;
/* ^^^ BUGFIX at Sat Jun 1 18:27:10 CEST 2002 */
if (TryOnly) return WarningOK || (nncols>4 && !hasTransp);
if (hasTransp) {
if (!WarningOK) return false;
Error::sev(Error::WARNING) <<
(nncols<=15 ? "SampleFormat: Transparent4 would be better than "
: "SampleFormat: Transparent8 may be better than ")
<< "Gray4" << (Error*)0;
}
if (nncols<=4) {
if (!WarningOK) return false;
Error::sev(Error::WARNING) << "SampleFormat: Gray2 would be better than Gray4" << (Error*)0;
}
img=img->toGray(4); if (bak!=img) delete bak;
assert(img!=NULLP);
sf=SF_Gray4; return true;
case SF_Indexed4:
if (nncols>16 || zero || hasTransp || img->hasPixelRGB(Transparent)) return false;
if (TryOnly) return WarningOK || (nncols>4 && minRGBBpc>=4);
if (nncols<=4) {
if (!WarningOK) return false;
Error::sev(Error::WARNING) << "SampleFormat: Indexed2 would be better than Indexed4" << (Error*)0;
}
if (minRGBBpc<=1) {
if (!WarningOK) return false;
Error::sev(Error::WARNING) << "SampleFormat: RGB1 would be better than Indexed4" << (Error*)0;
}
if (canGray && minRGBBpc<=4) {
if (!WarningOK) return false;
Error::sev(Error::NOTICE) << "SampleFormat: Gray4 would be better than Indexed4" << (Error*)0;
}
assert(img->getTy()==img->TY_INDEXED);
{ iimg=PTS_dynamic_cast(Indexed*,img);
iimg->setBpc(4);
if (iimg->wouldSetTranspc(Transparent)) return false;
iimg->setTranspcAndRepack(Transparent);
}
sf=SF_Indexed4; return true;
case SF_Transparent8:
if (!hasTransp && nncols==256 && PTS_dynamic_cast(Indexed*,img)->wouldSetTranspc(Transparent)) { hasTransp=true; --nncols; }
if (nncols>255 || zero) return false;
if (!WarningOK) return false;
if (TryOnly) return true;
Error::sev(Error::WARNING) << "SampleFormat: Transparent8 separates too many colors" << (Error*)0;
if (!hasTransp) {
Error::sev(Error::WARNING) << "SampleFormat: Indexed8 would be much better than Transparent8" << (Error*)0;
}
if (nncols<16) {
Error::sev(Error::WARNING) << "SampleFormat: Transparent4 would be better than Transparent8" << (Error*)0;
}
assert(img->getTy()==img->TY_INDEXED);
{ iimg=PTS_dynamic_cast(Indexed*,img);
iimg->setBpc(8); /* should be a no-op */
iimg->setTranspcAndRepack(Transparent);
// imgs=iimg->separate(); /* postponed because of GIF89a output */
}
sf=SF_Transparent8; return true;
case SF_Rgb2:
if (nncols>64 || minRGBBpc>2 || zero || hasTransp || img->hasPixelRGB(Transparent)) return false;
if (TryOnly) return WarningOK || (nncols>16 && !canGray && !hasTransp);
if (hasTransp) {
if (!WarningOK) return false;
Error::sev(Error::WARNING) << "SampleFormat: Transparent8 would be better than " << "Rgb2" << (Error*)0;
}
if (canGray) {
if (!WarningOK) return false;
Error::sev(Error::WARNING) << "SampleFormat: Gray2 would be better than Rgb2" << (Error*)0;
}
if (nncols<=16) {
if (!WarningOK) return false;
Error::sev(Error::WARNING) << "SampleFormat: "
<< (hasTransp ? "Transparent4 may" : "Indexed4 would")
<< " be better than Rgb2" << (Error*)0;
}
img=img->toRGB(2); if (bak!=img) delete bak;
assert(img!=NULLP);
sf=SF_Rgb2; return true;
case SF_Gray8:
if (nncols>256 || !canGray || zero || hasTransp || img->hasPixelRGB(Transparent)) return false;
if (TryOnly) return WarningOK || (nncols>16 && !hasTransp);
if (hasTransp) {
if (!WarningOK) return false;
Error::sev(Error::WARNING) <<
(nncols<=255 ? "SampleFormat: Transparent8 may be better than "
: "SampleFormat: ignoring transparency for ")
<< "Gray8" << (Error*)0;
}
if (minRGBBpc<=4) {
/* BUGFIX at Wed Jul 3 01:07:44 CEST 2002 */
if (!WarningOK) return false;
Error::sev(Error::WARNING) << "SampleFormat: Gray4 would be better than Gray8" << (Error*)0;
}
img=img->toGray(8); if (bak!=img) delete bak;
assert(img!=NULLP);
sf=SF_Gray8; return true;
case SF_Indexed8:
// fprintf(stderr, "nncols=%u hasTransp=%u zero=%u\n", nncols, hasTransp, zero);
if (nncols>256 || zero || hasTransp || img->hasPixelRGB(Transparent)) return false;
// assert(0);
if (TryOnly) return WarningOK || (nncols>16 && minRGBBpc>=8);
if (nncols<=16) {
if (!WarningOK) return false;
Error::sev(Error::WARNING) << "SampleFormat: Indexed4 would be better than Indexed8" << (Error*)0;
}
if (minRGBBpc<=2) {
if (!WarningOK) return false;
Error::sev(Error::WARNING) << "SampleFormat: RGB2 would be better than Indexed8" << (Error*)0;
}
if (canGray) {
if (!WarningOK) return false;
Error::sev(Error::NOTICE) << "SampleFormat: Gray8 would be better than Indexed8" << (Error*)0;
}
assert(img->getTy()==img->TY_INDEXED);
{ iimg=PTS_dynamic_cast(Indexed*,img);
iimg->setBpc(8); /* should be a no-op */
if (iimg->wouldSetTranspc(Transparent)) return false;
iimg->setTranspcAndRepack(Transparent);
}
sf=SF_Indexed8; return true;
case SF_Rgb4:
// fprintf(stderr, "minrgbbpc=%d to=%d\n", minRGBBpc, TryOnly);
if (minRGBBpc>4 || zero || hasTransp || img->hasPixelRGB(Transparent)) return false;
if (TryOnly) return WarningOK || (nncols>256 && !canGray && !hasTransp);
if (hasTransp) {
if (!WarningOK) return false;
Error::sev(Error::WARNING) <<
(nncols<=255 ? "SampleFormat: Transparent8 may be better than "
: "SampleFormat: ignoring transparency for ")
<< "Rgb4" << (Error*)0;
}
if (canGray) {
if (!WarningOK) return false;
Error::sev(Error::WARNING) << "SampleFormat: Gray4 would be better than Rgb4" << (Error*)0;
}
if (nncols<=256) {
if (!WarningOK) return false;
Error::sev(Error::WARNING) << "SampleFormat: "
<< (hasTransp && nncols<=255 ? "Transparent8 may" : "Indexed8 would")
<< " be better than Rgb4" << (Error*)0;
}
img=img->toRGB(4); if (bak!=img) delete bak;
assert(img!=NULLP);
sf=SF_Rgb4; return true;
case SF_Rgb8:
if (zero || hasTransp || img->hasPixelRGB(Transparent)) return false;
if (TryOnly) return WarningOK || (nncols>256 && !canGray && minRGBBpc>=8 && !hasTransp);
if (hasTransp) {
if (!WarningOK) return false;
Error::sev(Error::WARNING) <<
(nncols<=255 ? "SampleFormat: Transparent8 may be better than "
: "SampleFormat: ignoring transparency for ")
<< "Rgb8" << (Error*)0;
}
if (canGray) {
if (!WarningOK) return false;
Error::sev(Error::WARNING) << "SampleFormat: Gray8 would be better than Rgb8" << (Error*)0;
}
if (minRGBBpc<=4) {
if (!WarningOK) return false;
Error::sev(Error::WARNING) << "SampleFormat: Rgb4 would be better than Rgb8" << (Error*)0;
}
if (nncols<=256) {
if (!WarningOK) return false;
Error::sev(Error::WARNING) << "SampleFormat: "
<< (hasTransp && nncols<=255 ? "Transparent8 may" : "Indexed8 would")
<< " be better than Rgb8" << (Error*)0;
}
img=img->toRGB(8); if (bak!=img) delete bak;
assert(img!=NULLP);
sf=SF_Rgb8; return true;
case SF_Asis:
if (img->getTy()!=img->TY_BLACKBOX) {
Error::sev(Error::WARNING) << "SampleFormat: cannot convert image to /Asis" << (Error*)0;
return false;
}
sf=SF_Asis; return true;
}
assert(0 && "unknown SampleFormat requested");
return false; /* NOTREACHED */
}

/* __END__ */