granicus.if.org Git - imagemagick/blob - coders/palm.c

   1 /*
   2 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
   3 %                                                                             %
   4 %                                                                             %
   5 %                                                                             %
   6 %                        PPPP    AAA   L      M   M                           %
   7 %                        P   P  A   A  L      MM MM                           %
   8 %                        PPPP   AAAAA  L      M M M                           %
   9 %                        P      A   A  L      M   M                           %
  10 %                        P      A   A  LLLLL  M   M                           %
  11 %                                                                             %
  12 %                                                                             %
  13 %                          Read/Write Palm Pixmap.                            %
  14 %                                                                             %
  15 %                                                                             %
  16 %                              Software Design                                %
  17 %                            Christopher R. Hawks                             %
  18 %                               December 2001                                 %
  19 %                                                                             %
  20 %  Copyright 1999-2004 ImageMagick Studio LLC, a non-profit organization      %
  21 %  dedicated to making software imaging solutions freely available.           %
  22 %                                                                             %
  23 %  You may not use this file except in compliance with the License.  You may  %
  24 %  obtain a copy of the License at                                            %
  25 %                                                                             %
  26 %    http://www.imagemagick.org/script/license.php                            %
  27 %                                                                             %
  28 %  Unless required by applicable law or agreed to in writing, software        %
  29 %  distributed under the License is distributed on an "AS IS" BASIS,          %
  30 %  WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.   %
  31 %  See the License for the specific language governing permissions and        %
  32 %  limitations under the License.                                             %
  33 %                                                                             %
  34 %                                                                             %
  35 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  36 %
  37 %  Based on pnmtopalm by Bill Janssen and ppmtobmp by Ian Goldberg.
  38 %
  39 */
  40 \f
  41 /*
  42   Include declarations.
  43 */
  44 #include "MagickCore/studio.h"
  45 #include "MagickCore/attribute.h"
  46 #include "MagickCore/blob.h"
  47 #include "MagickCore/blob-private.h"
  48 #include "MagickCore/cache.h"
  49 #include "MagickCore/color.h"
  50 #include "MagickCore/colormap.h"
  51 #include "MagickCore/colormap-private.h"
  52 #include "MagickCore/color-private.h"
  53 #include "MagickCore/colorspace.h"
  54 #include "MagickCore/colorspace-private.h"
  55 #include "MagickCore/constitute.h"
  56 #include "MagickCore/exception.h"
  57 #include "MagickCore/histogram.h"
  58 #include "MagickCore/image.h"
  59 #include "MagickCore/image-private.h"
  60 #include "MagickCore/list.h"
  61 #include "MagickCore/magick.h"
  62 #include "MagickCore/memory_.h"
  63 #include "MagickCore/monitor.h"
  64 #include "MagickCore/monitor-private.h"
  65 #include "MagickCore/paint.h"
  66 #include "MagickCore/pixel-accessor.h"
  67 #include "MagickCore/property.h"
  68 #include "MagickCore/quantize.h"
  69 #include "MagickCore/quantum-private.h"
  70 #include "MagickCore/static.h"
  71 #include "MagickCore/string_.h"
  72 #include "MagickCore/module.h"
  73 #include "MagickCore/utility.h"
  74 \f
  75 /*
  76   Define declarations.
  77 */
  78 #define PALM_IS_COMPRESSED_FLAG  0x8000
  79 #define PALM_HAS_COLORMAP_FLAG  0x4000
  80 #define PALM_HAS_FOUR_BYTE_FIELD  0x0200
  81 #define PALM_HAS_TRANSPARENCY_FLAG  0x2000
  82 #define PALM_IS_INDIRECT  0x1000
  83 #define PALM_IS_FOR_SCREEN  0x0800
  84 #define PALM_IS_DIRECT_COLOR  0x0400
  85 #define PALM_COMPRESSION_SCANLINE  0x00
  86 #define PALM_COMPRESSION_RLE  0x01
  87 #define PALM_COMPRESSION_NONE  0xFF
  88 \f
  89 /*
  90  The 256 color system palette for Palm Computing Devices.
  91 */
  92 static const unsigned char
  93   PalmPalette[256][3] =
  94   {
  95     {255, 255,255}, {255, 204,255}, {255, 153,255}, {255, 102,255},
  96     {255,  51,255}, {255,   0,255}, {255, 255,204}, {255, 204,204},
  97     {255, 153,204}, {255, 102,204}, {255,  51,204}, {255,   0,204},
  98     {255, 255,153}, {255, 204,153}, {255, 153,153}, {255, 102,153},
  99     {255,  51,153}, {255,   0,153}, {204, 255,255}, {204, 204,255},
 100     {204, 153,255}, {204, 102,255}, {204,  51,255}, {204,   0,255},
 101     {204, 255,204}, {204, 204,204}, {204, 153,204}, {204, 102,204},
 102     {204,  51,204}, {204,   0,204}, {204, 255,153}, {204, 204,153},
 103     {204, 153,153}, {204, 102,153}, {204,  51,153}, {204,   0,153},
 104     {153, 255,255}, {153, 204,255}, {153, 153,255}, {153, 102,255},
 105     {153,  51,255}, {153,   0,255}, {153, 255,204}, {153, 204,204},
 106     {153, 153,204}, {153, 102,204}, {153,  51,204}, {153,   0,204},
 107     {153, 255,153}, {153, 204,153}, {153, 153,153}, {153, 102,153},
 108     {153,  51,153}, {153,   0,153}, {102, 255,255}, {102, 204,255},
 109     {102, 153,255}, {102, 102,255}, {102,  51,255}, {102,   0,255},
 110     {102, 255,204}, {102, 204,204}, {102, 153,204}, {102, 102,204},
 111     {102,  51,204}, {102,   0,204}, {102, 255,153}, {102, 204,153},
 112     {102, 153,153}, {102, 102,153}, {102,  51,153}, {102,   0,153},
 113     { 51, 255,255}, { 51, 204,255}, { 51, 153,255}, { 51, 102,255},
 114     { 51,  51,255}, { 51,   0,255}, { 51, 255,204}, { 51, 204,204},
 115     { 51, 153,204}, { 51, 102,204}, { 51,  51,204}, { 51,   0,204},
 116     { 51, 255,153}, { 51, 204,153}, { 51, 153,153}, { 51, 102,153},
 117     { 51,  51,153}, { 51,   0,153}, {  0, 255,255}, {  0, 204,255},
 118     {  0, 153,255}, {  0, 102,255}, {  0,  51,255}, {  0,   0,255},
 119     {  0, 255,204}, {  0, 204,204}, {  0, 153,204}, {  0, 102,204},
 120     {  0,  51,204}, {  0,   0,204}, {  0, 255,153}, {  0, 204,153},
 121     {  0, 153,153}, {  0, 102,153}, {  0,  51,153}, {  0,   0,153},
 122     {255, 255,102}, {255, 204,102}, {255, 153,102}, {255, 102,102},
 123     {255,  51,102}, {255,   0,102}, {255, 255, 51}, {255, 204, 51},
 124     {255, 153, 51}, {255, 102, 51}, {255,  51, 51}, {255,   0, 51},
 125     {255, 255,  0}, {255, 204,  0}, {255, 153,  0}, {255, 102,  0},
 126     {255,  51,  0}, {255,   0,  0}, {204, 255,102}, {204, 204,102},
 127     {204, 153,102}, {204, 102,102}, {204,  51,102}, {204,   0,102},
 128     {204, 255, 51}, {204, 204, 51}, {204, 153, 51}, {204, 102, 51},
 129     {204,  51, 51}, {204,   0, 51}, {204, 255,  0}, {204, 204,  0},
 130     {204, 153,  0}, {204, 102,  0}, {204,  51,  0}, {204,   0,  0},
 131     {153, 255,102}, {153, 204,102}, {153, 153,102}, {153, 102,102},
 132     {153,  51,102}, {153,   0,102}, {153, 255, 51}, {153, 204, 51},
 133     {153, 153, 51}, {153, 102, 51}, {153,  51, 51}, {153,   0, 51},
 134     {153, 255,  0}, {153, 204,  0}, {153, 153,  0}, {153, 102,  0},
 135     {153,  51,  0}, {153,   0,  0}, {102, 255,102}, {102, 204,102},
 136     {102, 153,102}, {102, 102,102}, {102,  51,102}, {102,   0,102},
 137     {102, 255, 51}, {102, 204, 51}, {102, 153, 51}, {102, 102, 51},
 138     {102,  51, 51}, {102,   0, 51}, {102, 255,  0}, {102, 204,  0},
 139     {102, 153,  0}, {102, 102,  0}, {102,  51,  0}, {102,   0,  0},
 140     { 51, 255,102}, { 51, 204,102}, { 51, 153,102}, { 51, 102,102},
 141     { 51,  51,102}, { 51,   0,102}, { 51, 255, 51}, { 51, 204, 51},
 142     { 51, 153, 51}, { 51, 102, 51}, { 51,  51, 51}, { 51,   0, 51},
 143     { 51, 255,  0}, { 51, 204,  0}, { 51, 153,  0}, { 51, 102,  0},
 144     { 51,  51,  0}, { 51,   0,  0}, {  0, 255,102}, {  0, 204,102},
 145     {  0, 153,102}, {  0, 102,102}, {  0,  51,102}, {  0,   0,102},
 146     {  0, 255, 51}, {  0, 204, 51}, {  0, 153, 51}, {  0, 102, 51},
 147     {  0,  51, 51}, {  0,   0, 51}, {  0, 255,  0}, {  0, 204,  0},
 148     {  0, 153,  0}, {  0, 102,  0}, {  0,  51,  0}, { 17,  17, 17},
 149     { 34,  34, 34}, { 68,  68, 68}, { 85,  85, 85}, {119, 119,119},
 150     {136, 136,136}, {170, 170,170}, {187, 187,187}, {221, 221,221},
 151     {238, 238,238}, {192, 192,192}, {128,   0,  0}, {128,   0,128},
 152     {  0, 128,  0}, {  0, 128,128}, {  0,   0,  0}, {  0,   0,  0},
 153     {  0,   0,  0}, {  0,   0,  0}, {  0,   0,  0}, {  0,   0,  0},
 154     {  0,   0,  0}, {  0,   0,  0}, {  0,   0,  0}, {  0,   0,  0},
 155     {  0,   0,  0}, {  0,   0,  0}, {  0,   0,  0}, {  0,   0,  0},
 156     {  0,   0,  0}, {  0,   0,  0}, {  0,   0,  0}, {  0,   0,  0},
 157     {  0,   0,  0}, {  0,   0,  0}, {  0,   0,  0}, {  0,   0,  0},
 158     {  0,   0,  0}, {  0,   0,  0}, {  0,   0,  0}, {  0,   0,  0}
 159   };
 160 \f
 161 /*
 162   Forward declarations.
 163 */
 164 static MagickBooleanType
 165   WritePALMImage(const ImageInfo *,Image *,ExceptionInfo *);
 166 \f
 167 /*
 168 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 169 %                                                                             %
 170 %                                                                             %
 171 %                                                                             %
 172 %   F i n d C o l o r                                                         %
 173 %                                                                             %
 174 %                                                                             %
 175 %                                                                             %
 176 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 177 %
 178 %  FindColor() returns the index of the matching entry from PalmPalette for a
 179 %  given PixelInfo.
 180 %
 181 %  The format of the FindColor method is:
 182 %
 183 %      int FindColor(PixelInfo *pixel)
 184 %
 185 %  A description of each parameter follows:
 186 %
 187 %    o int: the index of the matching color or -1 if not found/
 188 %
 189 %    o pixel: a pointer to the PixelInfo to be matched.
 190 %
 191 */
 192 static ssize_t FindColor(PixelInfo *packet)
 193 {
 194   register ssize_t
 195     i;
 196
 197   for (i=0; i < 256; i++)
 198     if (ScaleQuantumToChar(ClampToQuantum(packet->red)) == PalmPalette[i][0] &&
 199         ScaleQuantumToChar(ClampToQuantum(packet->green)) == PalmPalette[i][1] &&
 200         ScaleQuantumToChar(ClampToQuantum(packet->blue)) == PalmPalette[i][2])
 201       return(i);
 202   return(-1);
 203 }
 204 \f
 205 /*
 206 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 207 %                                                                             %
 208 %                                                                             %
 209 %                                                                             %
 210 %   R e a d P A L M I m a g e                                                 %
 211 %                                                                             %
 212 %                                                                             %
 213 %                                                                             %
 214 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 215 %
 216 %  ReadPALMImage() reads an image of raw bites in LSB order and returns it.  It
 217 %  allocates the memory necessary for the new Image structure and returns a
 218 %  pointer to the new image.
 219 %
 220 %  The format of the ReadPALMImage method is:
 221 %
 222 %      Image *ReadPALMImage(const ImageInfo *image_info,
 223 %        ExceptionInfo *exception)
 224 %
 225 %  A description of each parameter follows:
 226 %
 227 %    o image_info: Specifies a pointer to an ImageInfo structure.
 228 %
 229 %    o exception: return any errors or warnings in this structure.
 230 %
 231 */
 232 static Image *ReadPALMImage(const ImageInfo *image_info,
 233   ExceptionInfo *exception)
 234 {
 235   Image
 236     *image;
 237
 238   MagickBooleanType
 239     status;
 240
 241   MagickOffsetType
 242     totalOffset,
 243     seekNextDepth;
 244
 245   PixelInfo
 246     transpix;
 247
 248   Quantum
 249     index;
 250
 251   register ssize_t
 252     i,
 253     x;
 254
 255   register Quantum
 256     *q;
 257
 258   size_t
 259     bytes_per_row,
 260     flags,
 261     bits_per_pixel,
 262     version,
 263     nextDepthOffset,
 264     transparentIndex,
 265     compressionType,
 266     byte,
 267     mask,
 268     redbits,
 269     greenbits,
 270     bluebits,
 271     one,
 272     pad,
 273     size,
 274     bit;
 275
 276   ssize_t
 277     count,
 278     y;
 279
 280   unsigned char
 281     *lastrow,
 282     *one_row,
 283     *ptr;
 284
 285   unsigned short
 286     color16;
 287
 288   /*
 289     Open image file.
 290   */
 291   assert(image_info != (const ImageInfo *) NULL);
 292   assert(image_info->signature == MagickCoreSignature);
 293   if (image_info->debug != MagickFalse)
 294     (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",
 295       image_info->filename);
 296   assert(exception != (ExceptionInfo *) NULL);
 297   assert(exception->signature == MagickCoreSignature);
 298   image=AcquireImage(image_info,exception);
 299   status=OpenBlob(image_info,image,ReadBinaryBlobMode,exception);
 300   if (status == MagickFalse)
 301     {
 302       (void) DestroyImageList(image);
 303       return((Image *) NULL);
 304     }
 305   totalOffset=0;
 306   do
 307   {
 308     image->columns=ReadBlobMSBShort(image);
 309     image->rows=ReadBlobMSBShort(image);
 310     if (EOFBlob(image) != MagickFalse)
 311       ThrowReaderException(CorruptImageError,"ImproperImageHeader");
 312     if ((image->columns == 0) || (image->rows == 0))
 313       ThrowReaderException(CorruptImageError,"NegativeOrZeroImageSize");
 314     status=SetImageExtent(image,image->columns,image->rows,exception);
 315     if (status == MagickFalse)
 316       return(DestroyImageList(image));
 317     bytes_per_row=ReadBlobMSBShort(image);
 318     flags=ReadBlobMSBShort(image);
 319     bits_per_pixel=(size_t) ReadBlobByte(image);
 320     if ((bits_per_pixel != 1) && (bits_per_pixel != 2) &&
 321         (bits_per_pixel != 4) && (bits_per_pixel != 8) &&
 322         (bits_per_pixel != 16))
 323       ThrowReaderException(CorruptImageError,"UnrecognizedBitsPerPixel");
 324     version=(size_t) ReadBlobByte(image);
 325     if ((version != 0) && (version != 1) && (version != 2))
 326       ThrowReaderException(CorruptImageError,"FileFormatVersionMismatch");
 327     nextDepthOffset=(size_t) ReadBlobMSBShort(image);
 328     transparentIndex=(size_t) ReadBlobByte(image);
 329     compressionType=(size_t) ReadBlobByte(image);
 330     if ((compressionType != PALM_COMPRESSION_NONE) &&
 331         (compressionType != PALM_COMPRESSION_SCANLINE ) &&
 332         (compressionType != PALM_COMPRESSION_RLE))
 333       ThrowReaderException(CorruptImageError,"UnrecognizedImageCompression");
 334     pad=ReadBlobMSBShort(image);
 335     (void) pad;
 336     /*
 337       Initialize image colormap.
 338     */
 339     one=1;
 340     if ((bits_per_pixel < 16) &&
 341         (AcquireImageColormap(image,one << bits_per_pixel,exception) == MagickFalse))
 342       ThrowReaderException(ResourceLimitError,"MemoryAllocationFailed");
 343     GetPixelInfo(image,&transpix);
 344     if (bits_per_pixel == 16)  /* Direct Color */
 345       {
 346         redbits=(size_t) ReadBlobByte(image);  /* # of bits of red */
 347         (void) redbits;
 348         greenbits=(size_t) ReadBlobByte(image);  /* # of bits of green */
 349         (void) greenbits;
 350         bluebits=(size_t) ReadBlobByte(image);  /* # of bits of blue */
 351         (void) bluebits;
 352         ReadBlobByte(image);  /* reserved by Palm */
 353         ReadBlobByte(image);  /* reserved by Palm */
 354         transpix.red=(double) (QuantumRange*ReadBlobByte(image)/31);
 355         transpix.green=(double) (QuantumRange*ReadBlobByte(image)/63);
 356         transpix.blue=(double) (QuantumRange*ReadBlobByte(image)/31);
 357       }
 358     if (bits_per_pixel == 8)
 359       {
 360         ssize_t
 361           index;
 362
 363         if (flags & PALM_HAS_COLORMAP_FLAG)
 364           {
 365             count=(ssize_t) ReadBlobMSBShort(image);
 366             for (i=0; i < (ssize_t) count; i++)
 367             {
 368               ReadBlobByte(image);
 369               index=ConstrainColormapIndex(image,255-i,exception);
 370               image->colormap[index].red=(MagickRealType)
 371                 ScaleCharToQuantum((unsigned char) ReadBlobByte(image));
 372               image->colormap[index].green=(MagickRealType)
 373                 ScaleCharToQuantum((unsigned char) ReadBlobByte(image));
 374               image->colormap[index].blue=(MagickRealType)
 375                 ScaleCharToQuantum((unsigned char) ReadBlobByte(image));
 376           }
 377         }
 378       else
 379         for (i=0; i < (ssize_t) (1L << bits_per_pixel); i++)
 380         {
 381           index=ConstrainColormapIndex(image,255-i,exception);
 382           image->colormap[index].red=(MagickRealType)
 383             ScaleCharToQuantum(PalmPalette[i][0]);
 384           image->colormap[index].green=(MagickRealType)
 385             ScaleCharToQuantum(PalmPalette[i][1]);
 386           image->colormap[index].blue=(MagickRealType)
 387             ScaleCharToQuantum(PalmPalette[i][2]);
 388         }
 389       }
 390     if (flags & PALM_IS_COMPRESSED_FLAG)
 391       size=ReadBlobMSBShort(image);
 392     (void) size;
 393     image->storage_class=DirectClass;
 394     if (bits_per_pixel < 16)
 395       {
 396         image->storage_class=PseudoClass;
 397         image->depth=8;
 398       }
 399     if (image_info->ping != MagickFalse)
 400       {
 401         (void) CloseBlob(image);
 402         return(image);
 403       }
 404     status=SetImageExtent(image,image->columns,image->rows,exception);
 405     if (status == MagickFalse)
 406       return(DestroyImageList(image));
 407     one_row=(unsigned char *) AcquireQuantumMemory(MagickMax(bytes_per_row,
 408       2*image->columns),sizeof(*one_row));
 409     if (one_row == (unsigned char *) NULL)
 410       ThrowReaderException(ResourceLimitError,"MemoryAllocationFailed");
 411     lastrow=(unsigned char *) NULL;
 412     if (compressionType == PALM_COMPRESSION_SCANLINE) {
 413       lastrow=(unsigned char *) AcquireQuantumMemory(MagickMax(bytes_per_row,
 414         2*image->columns),sizeof(*lastrow));
 415     if (lastrow == (unsigned char *) NULL)
 416       ThrowReaderException(ResourceLimitError,"MemoryAllocationFailed");
 417     }
 418     mask=(size_t) (1U << bits_per_pixel)-1;
 419     for (y=0; y < (ssize_t) image->rows; y++)
 420     {
 421       if ((flags & PALM_IS_COMPRESSED_FLAG) == 0)
 422         {
 423           /* TODO move out of loop! */
 424           image->compression=NoCompression;
 425           count=ReadBlob(image,bytes_per_row,one_row);
 426           if (count != (ssize_t) bytes_per_row)
 427             break;
 428         }
 429       else
 430         {
 431           if (compressionType == PALM_COMPRESSION_RLE)
 432             {
 433               /* TODO move out of loop! */
 434               image->compression=RLECompression;
 435               for (i=0; i < (ssize_t) bytes_per_row; )
 436               {
 437                 count=(ssize_t) ReadBlobByte(image);
 438                 if (count < 0)
 439                   break;
 440                 count=MagickMin(count,(ssize_t) bytes_per_row-i);
 441                 byte=(size_t) ReadBlobByte(image);
 442                 (void) ResetMagickMemory(one_row+i,(int) byte,(size_t) count);
 443                 i+=count;
 444               }
 445           }
 446         else
 447           if (compressionType == PALM_COMPRESSION_SCANLINE)
 448             {
 449               size_t
 450                 one;
 451
 452               /* TODO move out of loop! */
 453               one=1;
 454               image->compression=FaxCompression;
 455               for (i=0; i < (ssize_t) bytes_per_row; i+=8)
 456               {
 457                 count=(ssize_t) ReadBlobByte(image);
 458                 if (count < 0)
 459                   break;
 460                 byte=(size_t) MagickMin((ssize_t) bytes_per_row-i,8);
 461                 for (bit=0; bit < byte; bit++)
 462                 {
 463                   if ((y == 0) || (count & (one << (7 - bit))))
 464                     one_row[i+bit]=(unsigned char) ReadBlobByte(image);
 465                   else
 466                     one_row[i+bit]=lastrow[i+bit];
 467                 }
 468               }
 469               (void) CopyMagickMemory(lastrow, one_row, bytes_per_row);
 470             }
 471         }
 472       ptr=one_row;
 473       q=QueueAuthenticPixels(image,0,y,image->columns,1,exception);
 474       if (q == (Quantum *) NULL)
 475         break;
 476       if (bits_per_pixel == 16)
 477         {
 478           if (image->columns > (2*bytes_per_row))
 479             ThrowReaderException(CorruptImageError,"CorruptImage");
 480           for (x=0; x < (ssize_t) image->columns; x++)
 481           {
 482             color16=(*ptr++ << 8);
 483             color16|=(*ptr++);
 484             SetPixelRed(image,(Quantum) ((QuantumRange*((color16 >> 11) &
 485               0x1f))/0x1f),q);
 486             SetPixelGreen(image,(Quantum) ((QuantumRange*((color16 >> 5) &
 487               0x3f))/0x3f),q);
 488             SetPixelBlue(image,(Quantum) ((QuantumRange*((color16 >> 0) &
 489               0x1f))/0x1f),q);
 490             SetPixelAlpha(image,OpaqueAlpha,q);
 491             q+=GetPixelChannels(image);
 492           }
 493         }
 494       else
 495         {
 496           bit=8-bits_per_pixel;
 497           for (x=0; x < (ssize_t) image->columns; x++)
 498           {
 499             if ((size_t) (ptr-one_row) >= bytes_per_row)
 500               ThrowReaderException(CorruptImageError,"CorruptImage");
 501             index=(Quantum) (mask-(((*ptr) & (mask << bit)) >> bit));
 502             SetPixelIndex(image,index,q);
 503             SetPixelViaPixelInfo(image,image->colormap+(ssize_t) index,q);
 504             if (bit)
 505               bit-=bits_per_pixel;
 506             else
 507               {
 508                 ptr++;
 509                 bit=8-bits_per_pixel;
 510               }
 511             q+=GetPixelChannels(image);
 512           }
 513           if (SyncAuthenticPixels(image,exception) == MagickFalse)
 514             break;
 515         }
 516         if (image->previous == (Image *) NULL)
 517           {
 518             status=SetImageProgress(image,LoadImageTag,(MagickOffsetType) y,
 519               image->rows);
 520             if (status == MagickFalse)
 521               break;
 522           }
 523       }
 524     if (flags & PALM_HAS_TRANSPARENCY_FLAG)
 525       {
 526         ssize_t index=ConstrainColormapIndex(image,(ssize_t) (mask-
 527           transparentIndex),exception);
 528         if (bits_per_pixel != 16)
 529           transpix=image->colormap[index];
 530         (void) TransparentPaintImage(image,&transpix,(Quantum) TransparentAlpha,
 531           MagickFalse,exception);
 532       }
 533     one_row=(unsigned char *) RelinquishMagickMemory(one_row);
 534     if (compressionType == PALM_COMPRESSION_SCANLINE)
 535       lastrow=(unsigned char *) RelinquishMagickMemory(lastrow);
 536     if (EOFBlob(image) != MagickFalse)
 537       {
 538         ThrowFileException(exception,CorruptImageError,"UnexpectedEndOfFile",
 539           image->filename);
 540         break;
 541       }
 542     /*
 543       Proceed to next image. Copied from coders/pnm.c
 544     */
 545     if (image_info->number_scenes != 0)
 546       if (image->scene >= (image_info->scene+image_info->number_scenes-1))
 547         break;
 548     if (nextDepthOffset != 0)
 549       {
 550         /*
 551           Skip to next image.
 552         */
 553         totalOffset+=(MagickOffsetType) (nextDepthOffset*4);
 554         if (totalOffset >= (MagickOffsetType) GetBlobSize(image))
 555           ThrowReaderException(CorruptImageError,"ImproperImageHeader")
 556         else
 557           seekNextDepth=SeekBlob(image,totalOffset,SEEK_SET);
 558         if (seekNextDepth != totalOffset)
 559           ThrowReaderException(CorruptImageError,"ImproperImageHeader");
 560         /*
 561           Allocate next image structure. Copied from coders/pnm.c
 562         */
 563         AcquireNextImage(image_info,image,exception);
 564         if (GetNextImageInList(image) == (Image *) NULL)
 565           {
 566             (void) DestroyImageList(image);
 567             return((Image *) NULL);
 568           }
 569         image=SyncNextImageInList(image);
 570         status=SetImageProgress(image,LoadImagesTag,TellBlob(image),
 571           GetBlobSize(image));
 572         if (status == MagickFalse)
 573           break;
 574       }
 575   } while (nextDepthOffset != 0);
 576   (void) CloseBlob(image);
 577   return(GetFirstImageInList(image));
 578 }
 579 \f
 580 /*
 581 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 582 %                                                                             %
 583 %                                                                             %
 584 %                                                                             %
 585 %   R e g i s t e r P A L M I m a g e                                         %
 586 %                                                                             %
 587 %                                                                             %
 588 %                                                                             %
 589 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 590 %
 591 %  RegisterPALMImage() adds properties for the PALM image format to the list of
 592 %  supported formats.  The properties include the image format tag, a method to
 593 %  read and/or write the format, whether the format supports the saving of more
 594 %  than one frame to the same file or blob, whether the format supports native
 595 %  in-memory I/O, and a brief description of the format.
 596 %
 597 %  The format of the RegisterPALMImage method is:
 598 %
 599 %      size_t RegisterPALMImage(void)
 600 %
 601 */
 602 ModuleExport size_t RegisterPALMImage(void)
 603 {
 604   MagickInfo
 605     *entry;
 606
 607   entry=AcquireMagickInfo("PALM","PALM","Palm pixmap");
 608   entry->decoder=(DecodeImageHandler *) ReadPALMImage;
 609   entry->encoder=(EncodeImageHandler *) WritePALMImage;
 610   entry->flags|=CoderDecoderSeekableStreamFlag;
 611   entry->flags|=CoderEncoderSeekableStreamFlag;
 612   (void) RegisterMagickInfo(entry);
 613   return(MagickImageCoderSignature);
 614 }
 615 \f
 616 /*
 617 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 618 %                                                                             %
 619 %                                                                             %
 620 %                                                                             %
 621 %   U n r e g i s t e r P A L M I m a g e                                     %
 622 %                                                                             %
 623 %                                                                             %
 624 %                                                                             %
 625 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 626 %
 627 %  UnregisterPALMImage() removes format registrations made by the PALM
 628 %  module from the list of supported formats.
 629 %
 630 %  The format of the UnregisterPALMImage method is:
 631 %
 632 %      UnregisterPALMImage(void)
 633 %
 634 */
 635 ModuleExport void UnregisterPALMImage(void)
 636 {
 637   (void) UnregisterMagickInfo("PALM");
 638 }
 639 \f
 640 /*
 641 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 642 %                                                                             %
 643 %                                                                             %
 644 %                                                                             %
 645 %   W r i t e P A L M I m a g e                                               %
 646 %                                                                             %
 647 %                                                                             %
 648 %                                                                             %
 649 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 650 %
 651 %  WritePALMImage() writes an image of raw bits in LSB order to a file.
 652 %
 653 %  The format of the WritePALMImage method is:
 654 %
 655 %      MagickBooleanType WritePALMImage(const ImageInfo *image_info,
 656 %        Image *image,ExceptionInfo *exception)
 657 %
 658 %  A description of each parameter follows.
 659 %
 660 %    o image_info: Specifies a pointer to an ImageInfo structure.
 661 %
 662 %    o image:  A pointer to a Image structure.
 663 %
 664 %    o exception: return any errors or warnings in this structure.
 665 %
 666 */
 667 static MagickBooleanType WritePALMImage(const ImageInfo *image_info,
 668   Image *image,ExceptionInfo *exception)
 669 {
 670   MagickBooleanType
 671     status;
 672
 673   MagickOffsetType
 674     currentOffset,
 675     offset,
 676     scene;
 677
 678   MagickSizeType
 679     cc;
 680
 681   PixelInfo
 682     transpix;
 683
 684   QuantizeInfo
 685     *quantize_info;
 686
 687   register ssize_t
 688     x;
 689
 690   register const Quantum
 691     *p;
 692
 693   register Quantum
 694     *q;
 695
 696   ssize_t
 697     y;
 698
 699   size_t
 700     count,
 701     bits_per_pixel,
 702     bytes_per_row,
 703     nextDepthOffset,
 704     one;
 705
 706   unsigned char
 707     bit,
 708     byte,
 709     color,
 710     *lastrow,
 711     *one_row,
 712     *ptr,
 713     version;
 714
 715   unsigned int
 716     transparentIndex;
 717
 718   unsigned short
 719     color16,
 720     flags;
 721
 722   /*
 723     Open output image file.
 724   */
 725   assert(image_info != (const ImageInfo *) NULL);
 726   assert(image_info->signature == MagickCoreSignature);
 727   assert(image != (Image *) NULL);
 728   assert(image->signature == MagickCoreSignature);
 729   if (image->debug != MagickFalse)
 730     (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
 731   assert(exception != (ExceptionInfo *) NULL);
 732   assert(exception->signature == MagickCoreSignature);
 733   status=OpenBlob(image_info,image,WriteBinaryBlobMode,exception);
 734   if (status == MagickFalse)
 735     return(status);
 736   quantize_info=AcquireQuantizeInfo(image_info);
 737   flags=0;
 738   currentOffset=0;
 739   transparentIndex=0;
 740   transpix.red=0.0;
 741   transpix.green=0.0;
 742   transpix.blue=0.0;
 743   transpix.alpha=0.0;
 744   one=1;
 745   version=0;
 746   scene=0;
 747   do
 748   {
 749     (void) TransformImageColorspace(image,sRGBColorspace,exception);
 750     count=GetNumberColors(image,NULL,exception);
 751     for (bits_per_pixel=1;  (one << bits_per_pixel) < count; bits_per_pixel*=2) ;
 752     if (image_info->depth > 100)
 753       bits_per_pixel=image_info->depth-100;
 754     if (bits_per_pixel < 16)
 755       (void) TransformImageColorspace(image,image->colorspace,exception);
 756     if (bits_per_pixel < 8)
 757       {
 758         (void) TransformImageColorspace(image,GRAYColorspace,exception);
 759         (void) SetImageType(image,PaletteType,exception);
 760         (void) SortColormapByIntensity(image,exception);
 761       }
 762     if ((image->storage_class == PseudoClass) && (image->colors > 256))
 763       (void) SetImageStorageClass(image,DirectClass,exception);
 764     if (image->storage_class == PseudoClass)
 765       flags|=PALM_HAS_COLORMAP_FLAG;
 766     else
 767       flags|=PALM_IS_DIRECT_COLOR;
 768     (void) WriteBlobMSBShort(image,(unsigned short) image->columns); /* width */
 769     (void) WriteBlobMSBShort(image,(unsigned short) image->rows);  /* height */
 770     bytes_per_row=((image->columns+(16/bits_per_pixel-1))/(16/
 771       bits_per_pixel))*2;
 772     (void) WriteBlobMSBShort(image,(unsigned short) bytes_per_row);
 773     if ((image_info->compression == RLECompression) ||
 774         (image_info->compression == FaxCompression))
 775       flags|=PALM_IS_COMPRESSED_FLAG;
 776     (void) WriteBlobMSBShort(image, flags);
 777     (void) WriteBlobByte(image,(unsigned char) bits_per_pixel);
 778     if (bits_per_pixel > 1)
 779       version=1;
 780     if ((image_info->compression == RLECompression) ||
 781         (image_info->compression == FaxCompression))
 782       version=2;
 783     (void) WriteBlobByte(image,version);
 784     (void) WriteBlobMSBShort(image,0);  /* nextDepthOffset */
 785     (void) WriteBlobByte(image,(unsigned char) transparentIndex);
 786     if (image_info->compression == RLECompression)
 787       (void) WriteBlobByte(image,PALM_COMPRESSION_RLE);
 788     else
 789       if (image_info->compression == FaxCompression)
 790         (void) WriteBlobByte(image,PALM_COMPRESSION_SCANLINE);
 791       else
 792         (void) WriteBlobByte(image,PALM_COMPRESSION_NONE);
 793     (void) WriteBlobMSBShort(image,0);  /* reserved */
 794     offset=16;
 795     if (bits_per_pixel == 16)
 796       {
 797         (void) WriteBlobByte(image,5);  /* # of bits of red */
 798         (void) WriteBlobByte(image,6);  /* # of bits of green */
 799         (void) WriteBlobByte(image,5);  /* # of bits of blue */
 800         (void) WriteBlobByte(image,0);  /* reserved by Palm */
 801         (void) WriteBlobMSBLong(image,0);  /* no transparent color, YET */
 802         offset+=8;
 803       }
 804     if (bits_per_pixel == 8)
 805       {
 806         if (flags & PALM_HAS_COLORMAP_FLAG)  /* Write out colormap */
 807           {
 808             quantize_info->dither_method=IdentifyPaletteImage(image,exception)
 809               == MagickFalse ? RiemersmaDitherMethod : NoDitherMethod;
 810             quantize_info->number_colors=image->colors;
 811             (void) QuantizeImage(quantize_info,image,exception);
 812             (void) WriteBlobMSBShort(image,(unsigned short) image->colors);
 813             for (count = 0; count < image->colors; count++)
 814             {
 815               (void) WriteBlobByte(image,(unsigned char) count);
 816               (void) WriteBlobByte(image,ScaleQuantumToChar(ClampToQuantum(
 817                 image->colormap[count].red)));
 818               (void) WriteBlobByte(image,ScaleQuantumToChar(ClampToQuantum(
 819                 image->colormap[count].green)));
 820               (void) WriteBlobByte(image,ScaleQuantumToChar(ClampToQuantum(
 821                 image->colormap[count].blue)));
 822             }
 823             offset+=2+count*4;
 824           }
 825       else  /* Map colors to Palm standard colormap */
 826         {
 827           Image
 828             *affinity_image;
 829
 830           affinity_image=ConstituteImage(256,1,"RGB",CharPixel,&PalmPalette,
 831             exception);
 832           (void) TransformImageColorspace(affinity_image,
 833             affinity_image->colorspace,exception);
 834           (void) RemapImage(quantize_info,image,affinity_image,exception);
 835           for (y=0; y < (ssize_t) image->rows; y++)
 836           {
 837             q=GetAuthenticPixels(image,0,y,image->columns,1,exception);
 838             if (q == (Quantum *) NULL)
 839               break;
 840             for (x=0; x < (ssize_t) image->columns; x++)
 841             {
 842               SetPixelIndex(image,(Quantum) FindColor(&image->colormap[(ssize_t)
 843                 GetPixelIndex(image,q)]),q);
 844               q+=GetPixelChannels(image);
 845             }
 846           }
 847           affinity_image=DestroyImage(affinity_image);
 848         }
 849       }
 850     if (flags & PALM_IS_COMPRESSED_FLAG)
 851       (void) WriteBlobMSBShort(image,0);  /* fill in size later */
 852     lastrow=(unsigned char *) NULL;
 853     if (image_info->compression == FaxCompression)
 854       lastrow=(unsigned char *) AcquireQuantumMemory(bytes_per_row,
 855         sizeof(*lastrow));
 856       /* TODO check whether memory really was acquired? */
 857     one_row=(unsigned char *) AcquireQuantumMemory(bytes_per_row,
 858       sizeof(*one_row));
 859     if (one_row == (unsigned char *) NULL)
 860       ThrowWriterException(ResourceLimitError,"MemoryAllocationFailed");
 861     for (y=0; y < (ssize_t) image->rows; y++)
 862     {
 863       ptr=one_row;
 864       (void) ResetMagickMemory(ptr,0,bytes_per_row);
 865       p=GetVirtualPixels(image,0,y,image->columns,1,exception);
 866       if (p == (const Quantum *) NULL)
 867         break;
 868       if (bits_per_pixel == 16)
 869         {
 870           for (x=0; x < (ssize_t) image->columns; x++)
 871           {
 872             color16=(unsigned short) ((((31*(size_t) GetPixelRed(image,p))/
 873               (size_t) QuantumRange) << 11) | (((63*(size_t)
 874               GetPixelGreen(image,p))/(size_t) QuantumRange) << 5) |
 875               ((31*(size_t) GetPixelBlue(image,p))/(size_t) QuantumRange));
 876             if (GetPixelAlpha(image,p) == (Quantum) TransparentAlpha)
 877               {
 878                 transpix.red=(MagickRealType) GetPixelRed(image,p);
 879                 transpix.green=(MagickRealType) GetPixelGreen(image,p);
 880                 transpix.blue=(MagickRealType) GetPixelBlue(image,p);
 881                 transpix.alpha=(MagickRealType) GetPixelAlpha(image,p);
 882                 flags|=PALM_HAS_TRANSPARENCY_FLAG;
 883               }
 884             *ptr++=(unsigned char) ((color16 >> 8) & 0xff);
 885             *ptr++=(unsigned char) (color16 & 0xff);
 886             p+=GetPixelChannels(image);
 887           }
 888         }
 889       else
 890         {
 891           byte=0x00;
 892           bit=(unsigned char) (8-bits_per_pixel);
 893           for (x=0; x < (ssize_t) image->columns; x++)
 894           {
 895             if (bits_per_pixel >= 8)
 896               color=(unsigned char) GetPixelIndex(image,p);
 897             else
 898               color=(unsigned char) (GetPixelIndex(image,p)*
 899                 ((one << bits_per_pixel)-1)/MagickMax(1*image->colors-1,1));
 900             byte|=color << bit;
 901             if (bit != 0)
 902               bit-=(unsigned char) bits_per_pixel;
 903             else
 904               {
 905                 *ptr++=byte;
 906                 byte=0x00;
 907                 bit=(unsigned char) (8-bits_per_pixel);
 908               }
 909             p+=GetPixelChannels(image);
 910           }
 911           if ((image->columns % (8/bits_per_pixel)) != 0)
 912             *ptr++=byte;
 913         }
 914       if (image_info->compression == RLECompression)
 915         {
 916           x=0;
 917           while (x < (ssize_t) bytes_per_row)
 918           {
 919             byte=one_row[x];
 920             count=1;
 921             while ((one_row[++x] == byte) && (count < 255) &&
 922                    (x < (ssize_t) bytes_per_row))
 923               count++;
 924             (void) WriteBlobByte(image,(unsigned char) count);
 925             (void) WriteBlobByte(image,(unsigned char) byte);
 926           }
 927         }
 928       else
 929         if (image_info->compression == FaxCompression)
 930           {
 931             char
 932               tmpbuf[8],
 933               *tptr;
 934
 935             for (x = 0;  x < (ssize_t) bytes_per_row;  x += 8)
 936             {
 937               tptr = tmpbuf;
 938               for (bit=0, byte=0; bit < (unsigned char) MagickMin(8,(ssize_t) bytes_per_row-x); bit++)
 939               {
 940                 if ((y == 0) || (lastrow[x + bit] != one_row[x + bit]))
 941                   {
 942                     byte |= (1 << (7 - bit));
 943                     *tptr++ = (char) one_row[x + bit];
 944                   }
 945               }
 946               (void) WriteBlobByte(image, byte);
 947               (void) WriteBlob(image,tptr-tmpbuf,(unsigned char *) tmpbuf);
 948             }
 949             (void) CopyMagickMemory(lastrow,one_row,bytes_per_row);
 950           }
 951         else
 952           (void) WriteBlob(image,bytes_per_row,one_row);
 953       }
 954     if (flags & PALM_HAS_TRANSPARENCY_FLAG)
 955       {
 956         offset=SeekBlob(image,currentOffset+6,SEEK_SET);
 957         (void) WriteBlobMSBShort(image,flags);
 958         offset=SeekBlob(image,currentOffset+12,SEEK_SET);
 959         (void) WriteBlobByte(image,(unsigned char) transparentIndex);  /* trans index */
 960       }
 961     if (bits_per_pixel == 16)
 962       {
 963         offset=SeekBlob(image,currentOffset+20,SEEK_SET);
 964         (void) WriteBlobByte(image,0);  /* reserved by Palm */
 965         (void) WriteBlobByte(image,(unsigned char) ((31*transpix.red)/QuantumRange));
 966         (void) WriteBlobByte(image,(unsigned char) ((63*transpix.green)/QuantumRange));
 967         (void) WriteBlobByte(image,(unsigned char) ((31*transpix.blue)/QuantumRange));
 968       }
 969     if (flags & PALM_IS_COMPRESSED_FLAG)  /* fill in size now */
 970       {
 971         offset=SeekBlob(image,currentOffset+offset,SEEK_SET);
 972         (void) WriteBlobMSBShort(image,(unsigned short) (GetBlobSize(image)-
 973           currentOffset-offset));
 974       }
 975     if (one_row != (unsigned char *) NULL)
 976       one_row=(unsigned char *) RelinquishMagickMemory(one_row);
 977     if (lastrow != (unsigned char *) NULL)
 978       lastrow=(unsigned char *) RelinquishMagickMemory(lastrow);
 979     if (GetNextImageInList(image) == (Image *) NULL)
 980       break;
 981     /* padding to 4 byte word */
 982     for (cc=(GetBlobSize(image)) % 4; cc > 0; cc--)
 983       (void) WriteBlobByte(image,0);
 984     /* write nextDepthOffset and return to end of image */
 985     offset=SeekBlob(image,currentOffset+10,SEEK_SET);
 986     nextDepthOffset=(size_t) ((GetBlobSize(image)-currentOffset)/4);
 987     (void) WriteBlobMSBShort(image,(unsigned short) nextDepthOffset);
 988     currentOffset=(MagickOffsetType) GetBlobSize(image);
 989     offset=SeekBlob(image,currentOffset,SEEK_SET);
 990     image=SyncNextImageInList(image);
 991     status=SetImageProgress(image,SaveImagesTag,scene++,
 992       GetImageListLength(image));
 993     if (status == MagickFalse)
 994       break;
 995   } while (image_info->adjoin != MagickFalse);
 996   quantize_info=DestroyQuantizeInfo(quantize_info);
 997   (void) CloseBlob(image);
 998   return(MagickTrue);
 999 }