]> granicus.if.org Git - imagemagick/blob - coders/palm.c
(no commit message)
[imagemagick] / 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 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
233 static inline size_t MagickMax(const size_t x,const size_t y)
234 {
235   if (x > y)
236     return(x);
237   return(y);
238 }
239
240 static inline ssize_t MagickMin(const ssize_t x,const ssize_t y)
241 {
242   if (x < y)
243     return(x);
244   return(y);
245 }
246
247 static Image *ReadPALMImage(const ImageInfo *image_info,
248   ExceptionInfo *exception)
249 {
250   Image
251     *image;
252
253   MagickBooleanType
254     status;
255
256   MagickOffsetType
257     totalOffset,
258     seekNextDepth;
259
260   PixelInfo
261     transpix;
262
263   Quantum
264     index;
265
266   register ssize_t
267     i,
268     x;
269
270   register Quantum
271     *q;
272
273   ssize_t
274     count,
275     y;
276
277   size_t
278     bytes_per_row,
279     flags,
280     bits_per_pixel,
281     version,
282     nextDepthOffset,
283     transparentIndex,
284     compressionType,
285     byte,
286     mask,
287     redbits,
288     greenbits,
289     bluebits,
290     one,
291     pad,
292     size,
293     bit;
294
295   unsigned char
296     *lastrow,
297     *one_row,
298     *ptr;
299
300   unsigned short
301     color16;
302
303   /*
304     Open image file.
305   */
306   assert(image_info != (const ImageInfo *) NULL);
307   assert(image_info->signature == MagickSignature);
308   if (image_info->debug != MagickFalse)
309     (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",
310       image_info->filename);
311   assert(exception != (ExceptionInfo *) NULL);
312   assert(exception->signature == MagickSignature);
313   image=AcquireImage(image_info,exception);
314   status=OpenBlob(image_info,image,ReadBinaryBlobMode,exception);
315   if (status == MagickFalse)
316     {
317       (void) DestroyImageList(image);
318       return((Image *) NULL);
319     }
320   totalOffset=0;
321   do
322   {
323     image->columns=ReadBlobMSBShort(image);
324     image->rows=ReadBlobMSBShort(image);
325     if (EOFBlob(image) != MagickFalse)
326       ThrowReaderException(CorruptImageError,"ImproperImageHeader");
327     if ((image->columns == 0) || (image->rows == 0))
328       ThrowReaderException(CorruptImageError,"NegativeOrZeroImageSize");
329     bytes_per_row=ReadBlobMSBShort(image);
330     flags=ReadBlobMSBShort(image);
331     bits_per_pixel=(size_t) ReadBlobByte(image);
332     if ((bits_per_pixel == 0) || (bits_per_pixel > 16))
333       ThrowReaderException(CorruptImageError,"ImproperImageHeader");
334     version=(size_t) ReadBlobByte(image);
335     (void) version;
336     nextDepthOffset=(size_t) ReadBlobMSBShort(image);
337     transparentIndex=(size_t) ReadBlobByte(image);
338     compressionType=(size_t) ReadBlobByte(image);
339     pad=ReadBlobMSBShort(image);
340     (void) pad;
341     /*
342       Initialize image colormap.
343     */
344     one=1;
345     if ((bits_per_pixel < 16) &&
346         (AcquireImageColormap(image,one << bits_per_pixel,exception) == MagickFalse))
347       ThrowReaderException(ResourceLimitError,"MemoryAllocationFailed");
348     GetPixelInfo(image,&transpix);
349     if (bits_per_pixel == 16)  /* Direct Color */
350       {
351         redbits=(size_t) ReadBlobByte(image);  /* # of bits of red */
352         (void) redbits;
353         greenbits=(size_t) ReadBlobByte(image);  /* # of bits of green */
354         (void) greenbits;
355         bluebits=(size_t) ReadBlobByte(image);  /* # of bits of blue */
356         (void) bluebits;
357         ReadBlobByte(image);  /* reserved by Palm */
358         ReadBlobByte(image);  /* reserved by Palm */
359         transpix.red=(double) (QuantumRange*ReadBlobByte(image)/31);
360         transpix.green=(double) (QuantumRange*ReadBlobByte(image)/63);
361         transpix.blue=(double) (QuantumRange*ReadBlobByte(image)/31);
362       }
363     if (bits_per_pixel == 8)
364       {
365         ssize_t
366           index;
367
368         if (flags & PALM_HAS_COLORMAP_FLAG)
369           {
370             count=(ssize_t) ReadBlobMSBShort(image);
371             for (i=0; i < (ssize_t) count; i++)
372             {
373               ReadBlobByte(image);
374               index=ConstrainColormapIndex(image,255-i,exception);
375               image->colormap[index].red=(MagickRealType)
376                 ScaleCharToQuantum((unsigned char) ReadBlobByte(image));
377               image->colormap[index].green=(MagickRealType)
378                 ScaleCharToQuantum((unsigned char) ReadBlobByte(image));
379               image->colormap[index].blue=(MagickRealType)
380                 ScaleCharToQuantum((unsigned char) ReadBlobByte(image));
381           }
382         }
383       else
384         {
385           for (i=0; i < (ssize_t) (1L << bits_per_pixel); i++)
386           {
387             index=ConstrainColormapIndex(image,255-i,exception);
388             image->colormap[index].red=(MagickRealType)
389               ScaleCharToQuantum(PalmPalette[i][0]);
390             image->colormap[index].green=(MagickRealType)
391               ScaleCharToQuantum(PalmPalette[i][1]);
392             image->colormap[index].blue=(MagickRealType)
393               ScaleCharToQuantum(PalmPalette[i][2]);
394           }
395         }
396       }
397     if (flags & PALM_IS_COMPRESSED_FLAG)
398       size=ReadBlobMSBShort(image);
399     (void) size;
400     image->storage_class=DirectClass;
401     if (bits_per_pixel < 16)
402       {
403         image->storage_class=PseudoClass;
404         image->depth=8;
405       }
406     one_row=(unsigned char *) AcquireQuantumMemory(MagickMax(bytes_per_row,
407       2*image->columns),sizeof(*one_row));
408     if (one_row == (unsigned char *) NULL)
409       ThrowReaderException(ResourceLimitError,"MemoryAllocationFailed");
410     lastrow=(unsigned char *) NULL;
411     if (compressionType == PALM_COMPRESSION_SCANLINE) {
412       lastrow=(unsigned char *) AcquireQuantumMemory(MagickMax(bytes_per_row,
413         2*image->columns),sizeof(*lastrow));
414     if (lastrow == (unsigned char *) NULL)
415       ThrowReaderException(ResourceLimitError,"MemoryAllocationFailed");
416     }
417     mask=(size_t) (1U << bits_per_pixel)-1;
418     for (y = 0; y < (ssize_t) image->rows; y++)
419     {
420       if ((flags & PALM_IS_COMPRESSED_FLAG) == 0)
421         {
422           /* TODO move out of loop! */
423           image->compression=NoCompression;
424           count=ReadBlob(image,bytes_per_row,one_row);
425         }
426       else
427         {
428           if (compressionType == PALM_COMPRESSION_RLE)
429             { 
430               /* TODO move out of loop! */
431               image->compression=RLECompression;
432               for (i=0; i < (ssize_t) bytes_per_row; )
433               {
434                 count=(ssize_t) ReadBlobByte(image);
435                 count=MagickMin(count,(ssize_t) bytes_per_row-i);
436                 byte=(size_t) ReadBlobByte(image);
437                 (void) ResetMagickMemory(one_row+i,(int) byte,(size_t) count);
438                 i+=count;
439               }
440           }
441         else
442           if (compressionType == PALM_COMPRESSION_SCANLINE)
443             {
444               size_t
445                 one;
446
447               /* TODO move out of loop! */
448               one=1;
449               image->compression=FaxCompression;
450               for (i=0; i < (ssize_t) bytes_per_row; i+=8)
451               {
452                 count=(ssize_t) ReadBlobByte(image);
453                 byte=(size_t) MagickMin((ssize_t) bytes_per_row-i,8);
454                 for (bit=0; bit < byte; bit++)
455                 {
456                   if ((y == 0) || (count & (one << (7 - bit))))
457                     one_row[i+bit]=(unsigned char) ReadBlobByte(image);
458                   else
459                     one_row[i+bit]=lastrow[i+bit];
460                 }
461               }
462               (void) CopyMagickMemory(lastrow, one_row, bytes_per_row);
463             }
464         }
465       ptr=one_row;
466       q=QueueAuthenticPixels(image,0,y,image->columns,1,exception);
467       if (q == (Quantum *) NULL)
468         break;
469       if (bits_per_pixel == 16)
470         {
471           if (image->columns > (2*bytes_per_row))
472             ThrowReaderException(CorruptImageError,"CorruptImage");
473           for (x=0; x < (ssize_t) image->columns; x++)
474           {
475             color16=(*ptr++ << 8);
476             color16|=(*ptr++);
477             SetPixelRed(image,(Quantum) ((QuantumRange*((color16 >> 11) &
478               0x1f))/0x1f),q);
479             SetPixelGreen(image,(Quantum) ((QuantumRange*((color16 >> 5) &
480               0x3f))/0x3f),q);
481             SetPixelBlue(image,(Quantum) ((QuantumRange*((color16 >> 0) &
482               0x1f))/0x1f),q);
483             SetPixelAlpha(image,OpaqueAlpha,q);
484             q+=GetPixelChannels(image);
485           }
486         }
487       else
488         {
489           bit=8-bits_per_pixel;
490           for (x=0; x < (ssize_t) image->columns; x++)
491           {
492             if ((size_t) (ptr-one_row) >= bytes_per_row)
493               ThrowReaderException(CorruptImageError,"CorruptImage");
494             index=(Quantum) (mask-(((*ptr) & (mask << bit)) >> bit));
495             SetPixelIndex(image,index,q);
496             SetPixelInfoPixel(image,image->colormap+(ssize_t) index,q);
497             if (bit)
498               bit-=bits_per_pixel;
499             else
500               {
501                 ptr++;
502                 bit=8-bits_per_pixel;
503               }
504             q+=GetPixelChannels(image);
505           }
506           if (SyncAuthenticPixels(image,exception) == MagickFalse)
507             break;
508         }
509         if (image->previous == (Image *) NULL)
510           {
511             status=SetImageProgress(image,LoadImageTag,(MagickOffsetType) y,
512               image->rows);
513             if (status == MagickFalse)
514               break;
515           }
516       }
517     if (flags & PALM_HAS_TRANSPARENCY_FLAG)
518       {
519         if (bits_per_pixel != 16)
520           transpix=image->colormap[mask-transparentIndex];
521         (void) TransparentPaintImage(image,&transpix,(Quantum) TransparentAlpha,
522           MagickFalse,exception);
523       }
524     one_row=(unsigned char *) RelinquishMagickMemory(one_row);
525     if (compressionType == PALM_COMPRESSION_SCANLINE)
526       lastrow=(unsigned char *) RelinquishMagickMemory(lastrow);
527     /*
528       Proceed to next image. Copied from coders/pnm.c
529     */
530     if (image_info->number_scenes != 0)
531       if (image->scene >= (image_info->scene+image_info->number_scenes-1))
532         break;
533     if (nextDepthOffset != 0)
534       { 
535         /*
536           Skip to next image. 
537         */
538         totalOffset+=(MagickOffsetType) (nextDepthOffset*4);
539         if (totalOffset >= (MagickOffsetType) GetBlobSize(image))
540           {
541             ThrowReaderException(CorruptImageError,"ImproperImageHeader");
542           }
543         else
544           {
545             seekNextDepth=SeekBlob(image,totalOffset,SEEK_SET);
546           }
547         if (seekNextDepth != totalOffset)
548           ThrowReaderException(CorruptImageError,"ImproperImageHeader");
549         /*
550           Allocate next image structure. Copied from coders/pnm.c
551         */
552         AcquireNextImage(image_info,image,exception);
553         if (GetNextImageInList(image) == (Image *) NULL)
554           {
555             (void) DestroyImageList(image);
556             return((Image *) NULL);
557           }
558         image=SyncNextImageInList(image);
559         status=SetImageProgress(image,LoadImagesTag,TellBlob(image),
560           GetBlobSize(image));
561         if (status == MagickFalse)
562           break;
563       }
564   } while (nextDepthOffset != 0);
565   (void) CloseBlob(image);
566   return(GetFirstImageInList(image));
567 }
568 \f
569 /*
570 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
571 %                                                                             %
572 %                                                                             %
573 %                                                                             %
574 %   R e g i s t e r P A L M I m a g e                                         %
575 %                                                                             %
576 %                                                                             %
577 %                                                                             %
578 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
579 %
580 %  RegisterPALMImage() adds properties for the PALM image format to the list of
581 %  supported formats.  The properties include the image format tag, a method to
582 %  read and/or write the format, whether the format supports the saving of more
583 %  than one frame to the same file or blob, whether the format supports native
584 %  in-memory I/O, and a brief description of the format.
585 %
586 %  The format of the RegisterPALMImage method is:
587 %
588 %      size_t RegisterPALMImage(void)
589 %
590 */
591 ModuleExport size_t RegisterPALMImage(void)
592 {
593   MagickInfo
594     *entry;
595
596   entry=SetMagickInfo("PALM");
597   entry->decoder=(DecodeImageHandler *) ReadPALMImage;
598   entry->encoder=(EncodeImageHandler *) WritePALMImage;
599   entry->seekable_stream=MagickTrue;
600   entry->description=ConstantString("Palm pixmap");
601   entry->module=ConstantString("PALM");
602   (void) RegisterMagickInfo(entry);
603   return(MagickImageCoderSignature);
604 }
605 \f
606 /*
607 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
608 %                                                                             %
609 %                                                                             %
610 %                                                                             %
611 %   U n r e g i s t e r P A L M I m a g e                                     %
612 %                                                                             %
613 %                                                                             %
614 %                                                                             %
615 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
616 %
617 %  UnregisterPALMImage() removes format registrations made by the PALM
618 %  module from the list of supported formats.
619 %
620 %  The format of the UnregisterPALMImage method is:
621 %
622 %      UnregisterPALMImage(void)
623 %
624 */
625 ModuleExport void UnregisterPALMImage(void)
626 {
627   (void) UnregisterMagickInfo("PALM");
628 }
629 \f
630 /*
631 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
632 %                                                                             %
633 %                                                                             %
634 %                                                                             %
635 %   W r i t e P A L M I m a g e                                               %
636 %                                                                             %
637 %                                                                             %
638 %                                                                             %
639 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
640 %
641 %  WritePALMImage() writes an image of raw bits in LSB order to a file.
642 %
643 %  The format of the WritePALMImage method is:
644 %
645 %      MagickBooleanType WritePALMImage(const ImageInfo *image_info,
646 %        Image *image,ExceptionInfo *exception)
647 %
648 %  A description of each parameter follows.
649 %
650 %    o image_info: Specifies a pointer to an ImageInfo structure.
651 %
652 %    o image:  A pointer to a Image structure.
653 %
654 %    o exception: return any errors or warnings in this structure.
655 %
656 */
657 static MagickBooleanType WritePALMImage(const ImageInfo *image_info,
658   Image *image,ExceptionInfo *exception)
659 {
660   MagickBooleanType
661     status;
662
663   MagickOffsetType
664     currentOffset,
665     offset,
666     scene;
667
668   MagickSizeType
669     cc;
670
671   PixelInfo
672     transpix;
673
674   QuantizeInfo
675     *quantize_info;
676
677   register ssize_t
678     x;
679
680   register const Quantum
681     *p;
682
683   register Quantum
684     *q;
685
686   ssize_t
687     y;
688
689   size_t
690     count,
691     bits_per_pixel,
692     bytes_per_row,
693     nextDepthOffset,
694     one;
695
696   unsigned char
697     bit,
698     byte,
699     color,
700     *lastrow,
701     *one_row,
702     *ptr,
703     version;
704
705   unsigned int
706     transparentIndex;
707
708   unsigned short
709     color16,
710     flags;
711
712   /*
713     Open output image file.
714   */
715   assert(image_info != (const ImageInfo *) NULL);
716   assert(image_info->signature == MagickSignature);
717   assert(image != (Image *) NULL);
718   assert(image->signature == MagickSignature);
719   if (image->debug != MagickFalse)
720     (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
721   assert(exception != (ExceptionInfo *) NULL);
722   assert(exception->signature == MagickSignature);
723   status=OpenBlob(image_info,image,WriteBinaryBlobMode,exception);
724   if (status == MagickFalse)
725     return(status);
726   quantize_info=AcquireQuantizeInfo(image_info);
727   flags=0;
728   currentOffset=0;
729   transparentIndex=0;
730   transpix.red=0.0;
731   transpix.green=0.0;
732   transpix.blue=0.0;
733   transpix.alpha=0.0;
734   one=1;
735   version=0;
736   scene=0;
737   do
738   {
739     (void) TransformImageColorspace(image,sRGBColorspace,exception);
740     count=GetNumberColors(image,NULL,exception);
741     for (bits_per_pixel=1;  (one << bits_per_pixel) < count; bits_per_pixel*=2) ;
742     if (image_info->depth > 100)
743       bits_per_pixel=image_info->depth-100;
744     if (bits_per_pixel < 16)
745       (void) TransformImageColorspace(image,image->colorspace,exception);
746     if (bits_per_pixel < 8)
747       {
748         (void) TransformImageColorspace(image,GRAYColorspace,exception);
749         (void) SetImageType(image,PaletteType,exception);
750         (void) SortColormapByIntensity(image,exception);
751       }
752     if ((image->storage_class == PseudoClass) && (image->colors > 256))
753       (void) SetImageStorageClass(image,DirectClass,exception);
754     if (image->storage_class == PseudoClass)
755       flags|=PALM_HAS_COLORMAP_FLAG;
756     else
757       flags|=PALM_IS_DIRECT_COLOR;
758     (void) WriteBlobMSBShort(image,(unsigned short) image->columns); /* width */
759     (void) WriteBlobMSBShort(image,(unsigned short) image->rows);  /* height */
760     bytes_per_row=((image->columns+(16/bits_per_pixel-1))/(16/
761       bits_per_pixel))*2;
762     (void) WriteBlobMSBShort(image,(unsigned short) bytes_per_row);
763     if ((image_info->compression == RLECompression) ||
764         (image_info->compression == FaxCompression))
765       flags|=PALM_IS_COMPRESSED_FLAG;
766     (void) WriteBlobMSBShort(image, flags);
767     (void) WriteBlobByte(image,(unsigned char) bits_per_pixel);
768     if (bits_per_pixel > 1)
769       version=1;
770     if ((image_info->compression == RLECompression) ||
771         (image_info->compression == FaxCompression))
772       version=2;
773     (void) WriteBlobByte(image,version);
774     (void) WriteBlobMSBShort(image,0);  /* nextDepthOffset */
775     (void) WriteBlobByte(image,(unsigned char) transparentIndex);
776     if (image_info->compression == RLECompression)
777       (void) WriteBlobByte(image,PALM_COMPRESSION_RLE);
778     else
779       if (image_info->compression == FaxCompression)
780         (void) WriteBlobByte(image,PALM_COMPRESSION_SCANLINE);
781       else
782         (void) WriteBlobByte(image,PALM_COMPRESSION_NONE);
783     (void) WriteBlobMSBShort(image,0);  /* reserved */
784     offset=16;
785     if (bits_per_pixel == 16)
786       {
787         (void) WriteBlobByte(image,5);  /* # of bits of red */
788         (void) WriteBlobByte(image,6);  /* # of bits of green */
789         (void) WriteBlobByte(image,5);  /* # of bits of blue */
790         (void) WriteBlobByte(image,0);  /* reserved by Palm */
791         (void) WriteBlobMSBLong(image,0);  /* no transparent color, YET */
792         offset+=8;
793       }
794     if (bits_per_pixel == 8)
795       {
796         if (flags & PALM_HAS_COLORMAP_FLAG)  /* Write out colormap */
797           {
798             quantize_info->dither_method=IsPaletteImage(image,exception)
799               == MagickFalse ? RiemersmaDitherMethod : NoDitherMethod;
800             quantize_info->number_colors=image->colors;
801             (void) QuantizeImage(quantize_info,image,exception);
802             (void) WriteBlobMSBShort(image,(unsigned short) image->colors);
803             for (count = 0; count < image->colors; count++)
804             {
805               (void) WriteBlobByte(image,(unsigned char) count);
806               (void) WriteBlobByte(image,ScaleQuantumToChar(ClampToQuantum(
807                 image->colormap[count].red)));
808               (void) WriteBlobByte(image,ScaleQuantumToChar(ClampToQuantum(
809                 image->colormap[count].green)));
810               (void) WriteBlobByte(image,ScaleQuantumToChar(ClampToQuantum(
811                 image->colormap[count].blue)));
812             }
813             offset+=2+count*4;
814           }
815       else  /* Map colors to Palm standard colormap */
816         {
817           Image
818             *affinity_image;
819
820           affinity_image=ConstituteImage(256,1,"RGB",CharPixel,&PalmPalette,
821             exception);
822           (void) TransformImageColorspace(affinity_image,
823             affinity_image->colorspace,exception);
824           (void) RemapImage(quantize_info,image,affinity_image,exception);
825           for (y=0; y < (ssize_t) image->rows; y++)
826           {
827             q=GetAuthenticPixels(image,0,y,image->columns,1,exception);
828             for (x=0; x < (ssize_t) image->columns; x++)
829             {
830               SetPixelIndex(image,(Quantum) FindColor(&image->colormap[(ssize_t)
831                 GetPixelIndex(image,q)]),q);
832               q+=GetPixelChannels(image);
833             }
834           }
835           affinity_image=DestroyImage(affinity_image);
836         }
837       }
838     if (flags & PALM_IS_COMPRESSED_FLAG)
839       (void) WriteBlobMSBShort(image,0);  /* fill in size later */
840     lastrow=(unsigned char *) NULL;
841     if (image_info->compression == FaxCompression)
842       lastrow=(unsigned char *) AcquireQuantumMemory(bytes_per_row,
843         sizeof(*lastrow));
844       /* TODO check whether memory really was acquired? */
845     one_row=(unsigned char *) AcquireQuantumMemory(bytes_per_row,
846       sizeof(*one_row));
847     if (one_row == (unsigned char *) NULL)
848       ThrowWriterException(ResourceLimitError,"MemoryAllocationFailed");
849     for (y=0; y < (ssize_t) image->rows; y++)
850     {
851       ptr=one_row;
852       (void) ResetMagickMemory(ptr,0,bytes_per_row);
853       p=GetVirtualPixels(image,0,y,image->columns,1,exception);
854       if (p == (const Quantum *) NULL)
855         break;
856       if (bits_per_pixel == 16)
857         {
858           for (x=0; x < (ssize_t) image->columns; x++)
859           {
860             color16=(unsigned short) ((((31*(size_t) GetPixelRed(image,p))/
861               (size_t) QuantumRange) << 11) | (((63*(size_t)
862               GetPixelGreen(image,p))/(size_t) QuantumRange) << 5) |
863               ((31*(size_t) GetPixelBlue(image,p))/(size_t) QuantumRange));
864             if (GetPixelAlpha(image,p) == (Quantum) TransparentAlpha)
865               {
866                 transpix.red=(MagickRealType) GetPixelRed(image,p);
867                 transpix.green=(MagickRealType) GetPixelGreen(image,p);
868                 transpix.blue=(MagickRealType) GetPixelBlue(image,p);
869                 transpix.alpha=(MagickRealType) GetPixelAlpha(image,p);
870                 flags|=PALM_HAS_TRANSPARENCY_FLAG;
871               }
872             *ptr++=(unsigned char) ((color16 >> 8) & 0xff);
873             *ptr++=(unsigned char) (color16 & 0xff);
874             p+=GetPixelChannels(image);
875           }
876         }
877       else
878         {
879           byte=0x00;
880           bit=(unsigned char) (8-bits_per_pixel);
881           for (x=0; x < (ssize_t) image->columns; x++)
882           {
883             if (bits_per_pixel >= 8)
884               color=(unsigned char) GetPixelIndex(image,p);
885             else
886               color=(unsigned char) (GetPixelIndex(image,p)*
887                 ((one << bits_per_pixel)-1)/MagickMax(1*image->colors-1,1));
888             byte|=color << bit;
889             if (bit != 0)
890               bit-=(unsigned char) bits_per_pixel;
891             else
892               {
893                 *ptr++=byte;
894                 byte=0x00;
895                 bit=(unsigned char) (8-bits_per_pixel);
896               }
897             p+=GetPixelChannels(image);
898           }
899           if ((image->columns % (8/bits_per_pixel)) != 0)
900             *ptr++=byte;
901         }
902       if (image_info->compression == RLECompression)
903         {
904           x=0;
905           while (x < (ssize_t) bytes_per_row)
906           {
907             byte=one_row[x];
908             count=1;
909             while ((one_row[++x] == byte) && (count < 255) &&
910                    (x < (ssize_t) bytes_per_row))
911               count++;
912             (void) WriteBlobByte(image,(unsigned char) count);
913             (void) WriteBlobByte(image,(unsigned char) byte);
914           }
915         }
916       else
917         if (image_info->compression == FaxCompression)
918           {
919             char
920               tmpbuf[8],
921               *tptr;
922   
923             for (x = 0;  x < (ssize_t) bytes_per_row;  x += 8)
924             {
925               tptr = tmpbuf;
926               for (bit=0, byte=0; bit < (unsigned char) MagickMin(8,(ssize_t) bytes_per_row-x); bit++)
927               {
928                 if ((y == 0) || (lastrow[x + bit] != one_row[x + bit]))
929                   {
930                     byte |= (1 << (7 - bit));
931                     *tptr++ = (char) one_row[x + bit];
932                   }
933               }
934               (void) WriteBlobByte(image, byte);
935               (void) WriteBlob(image,tptr-tmpbuf,(unsigned char *) tmpbuf);
936             }
937             (void) CopyMagickMemory(lastrow,one_row,bytes_per_row);
938           }
939         else
940           (void) WriteBlob(image,bytes_per_row,one_row);
941       }
942     if (flags & PALM_HAS_TRANSPARENCY_FLAG)
943       {
944         offset=SeekBlob(image,currentOffset+6,SEEK_SET);
945         (void) WriteBlobMSBShort(image,flags);
946         offset=SeekBlob(image,currentOffset+12,SEEK_SET);
947         (void) WriteBlobByte(image,(unsigned char) transparentIndex);  /* trans index */
948       }
949     if (bits_per_pixel == 16)
950       {
951         offset=SeekBlob(image,currentOffset+20,SEEK_SET);
952         (void) WriteBlobByte(image,0);  /* reserved by Palm */
953         (void) WriteBlobByte(image,(unsigned char) ((31*transpix.red)/QuantumRange));
954         (void) WriteBlobByte(image,(unsigned char) ((63*transpix.green)/QuantumRange));
955         (void) WriteBlobByte(image,(unsigned char) ((31*transpix.blue)/QuantumRange));
956       }
957     if (flags & PALM_IS_COMPRESSED_FLAG)  /* fill in size now */
958       {
959         offset=SeekBlob(image,currentOffset+offset,SEEK_SET);
960         (void) WriteBlobMSBShort(image,(unsigned short) (GetBlobSize(image)-
961           currentOffset-offset));
962       }
963     if (one_row != (unsigned char *) NULL) 
964       one_row=(unsigned char *) RelinquishMagickMemory(one_row);
965     if (lastrow != (unsigned char *) NULL) 
966       lastrow=(unsigned char *) RelinquishMagickMemory(lastrow);
967     if (GetNextImageInList(image) == (Image *) NULL)
968       break;
969     /* padding to 4 byte word */
970     for (cc = (GetBlobSize(image))%4; cc > 0; cc--)
971     {
972       (void) WriteBlobByte(image,0);
973     } 
974     /* write nextDepthOffset and return to end of image */
975     offset=SeekBlob(image,currentOffset+10,SEEK_SET);
976     nextDepthOffset=(size_t) ((GetBlobSize(image)-currentOffset)/4);
977     (void) WriteBlobMSBShort(image,(unsigned short) nextDepthOffset);
978     currentOffset=(MagickOffsetType) GetBlobSize(image);
979     offset=SeekBlob(image,currentOffset,SEEK_SET);
980     image=SyncNextImageInList(image);
981     status=SetImageProgress(image,SaveImagesTag,scene++,
982       GetImageListLength(image));
983     if (status == MagickFalse)
984       break;
985   } while (image_info->adjoin != MagickFalse);
986   quantize_info=DestroyQuantizeInfo(quantize_info);
987   (void) CloseBlob(image);
988   return(MagickTrue);
989 }