2 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
6 % FFFFF IIIII TTTTT SSSSS %
13 % Read/Write Flexible Image Transport System Images. %
20 % Copyright 1999-2012 ImageMagick Studio LLC, a non-profit organization %
21 % dedicated to making software imaging solutions freely available. %
23 % You may not use this file except in compliance with the License. You may %
24 % obtain a copy of the License at %
26 % http://www.imagemagick.org/script/license.php %
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. %
34 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
42 #include "MagickCore/studio.h"
43 #include "MagickCore/attribute.h"
44 #include "MagickCore/blob.h"
45 #include "MagickCore/blob-private.h"
46 #include "MagickCore/cache.h"
47 #include "MagickCore/color-private.h"
48 #include "MagickCore/colorspace.h"
49 #include "MagickCore/colorspace-private.h"
50 #include "MagickCore/constitute.h"
51 #include "MagickCore/exception.h"
52 #include "MagickCore/exception-private.h"
53 #include "MagickCore/image.h"
54 #include "MagickCore/image-private.h"
55 #include "MagickCore/list.h"
56 #include "MagickCore/magick.h"
57 #include "MagickCore/memory_.h"
58 #include "MagickCore/module.h"
59 #include "MagickCore/monitor.h"
60 #include "MagickCore/monitor-private.h"
61 #include "MagickCore/pixel-accessor.h"
62 #include "MagickCore/property.h"
63 #include "MagickCore/quantum-private.h"
64 #include "MagickCore/static.h"
65 #include "MagickCore/statistic.h"
66 #include "MagickCore/string_.h"
67 #include "MagickCore/string-private.h"
68 #include "MagickCore/module.h"
73 #define FITSBlocksize 2880UL
78 static MagickBooleanType
79 WriteFITSImage(const ImageInfo *,Image *,ExceptionInfo *);
82 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
90 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
92 % IsFITS() returns MagickTrue if the image format type, identified by the
93 % magick string, is FITS.
95 % The format of the IsFITS method is:
97 % MagickBooleanType IsFITS(const unsigned char *magick,const size_t length)
99 % A description of each parameter follows:
101 % o magick: compare image format pattern against these bytes.
103 % o length: Specifies the length of the magick string.
106 static MagickBooleanType IsFITS(const unsigned char *magick,const size_t length)
110 if (LocaleNCompare((const char *) magick,"IT0",3) == 0)
112 if (LocaleNCompare((const char *) magick,"SIMPLE",6) == 0)
118 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
122 % R e a d F I T S I m a g e %
126 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
128 % ReadFITSImage() reads a FITS image file and returns it. It allocates the
129 % memory necessary for the new Image structure and returns a pointer to the
132 % The format of the ReadFITSImage method is:
134 % Image *ReadFITSImage(const ImageInfo *image_info,
135 % ExceptionInfo *exception)
137 % A description of each parameter follows:
139 % o image_info: the image info.
141 % o exception: return any errors or warnings in this structure.
145 static inline double GetFITSPixel(Image *image,int bits_per_pixel)
147 switch (image->depth >> 3)
150 return((double) ReadBlobByte(image));
152 return((double) ((short) ReadBlobShort(image)));
155 if (bits_per_pixel > 0)
156 return((double) ((int) ReadBlobLong(image)));
157 return((double) ReadBlobFloat(image));
161 if (bits_per_pixel > 0)
162 return((double) ((MagickOffsetType) ReadBlobLongLong(image)));
167 return(ReadBlobDouble(image));
170 static void GetFITSPixelExtrema(Image *image,const int bits_per_pixel,
171 double *minima,double *maxima)
182 register MagickOffsetType
185 offset=TellBlob(image);
186 number_pixels=(MagickSizeType) image->columns*image->rows;
187 *minima=GetFITSPixel(image,bits_per_pixel);
189 for (i=1; i < (MagickOffsetType) number_pixels; i++)
191 pixel=GetFITSPixel(image,bits_per_pixel);
197 (void) SeekBlob(image,offset,SEEK_SET);
200 static inline double GetFITSPixelRange(const size_t depth)
202 return((double) ((MagickOffsetType) GetQuantumRange(depth)));
205 static void SetFITSUnsignedPixels(const size_t length,
206 const size_t bits_per_pixel,unsigned char *pixels)
211 for (i=0; i < (ssize_t) length; i++)
214 pixels+=bits_per_pixel >> 3;
218 static Image *ReadFITSImage(const ImageInfo *image_info,
219 ExceptionInfo *exception)
221 typedef struct _FITSInfo
247 property[MaxTextExtent],
284 assert(image_info != (const ImageInfo *) NULL);
285 assert(image_info->signature == MagickSignature);
286 if (image_info->debug != MagickFalse)
287 (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",
288 image_info->filename);
289 assert(exception != (ExceptionInfo *) NULL);
290 assert(exception->signature == MagickSignature);
291 image=AcquireImage(image_info,exception);
292 status=OpenBlob(image_info,image,ReadBinaryBlobMode,exception);
293 if (status == MagickFalse)
295 image=DestroyImageList(image);
296 return((Image *) NULL);
299 Initialize image header.
301 (void) ResetMagickMemory(&fits_info,0,sizeof(fits_info));
302 fits_info.extend=MagickFalse;
303 fits_info.simple=MagickFalse;
304 fits_info.bits_per_pixel=8;
308 fits_info.number_planes=1;
309 fits_info.min_data=0.0;
310 fits_info.max_data=0.0;
313 fits_info.endian=MSBEndian;
317 for (comment=(char *) NULL; EOFBlob(image) == MagickFalse; )
319 for ( ; EOFBlob(image) == MagickFalse; )
324 count=ReadBlob(image,8,(unsigned char *) keyword);
327 for (i=0; i < 8; i++)
329 if (isspace((int) ((unsigned char) keyword[i])) != 0)
331 keyword[i]=tolower((int) ((unsigned char) keyword[i]));
334 count=ReadBlob(image,72,(unsigned char *) value);
342 while (isspace((int) ((unsigned char) *p)) != 0)
345 if (LocaleCompare(keyword,"end") == 0)
347 if (LocaleCompare(keyword,"extend") == 0)
348 fits_info.extend=(*p == 'T') || (*p == 't') ? MagickTrue : MagickFalse;
349 if (LocaleCompare(keyword,"simple") == 0)
350 fits_info.simple=(*p == 'T') || (*p == 't') ? MagickTrue : MagickFalse;
351 if (LocaleCompare(keyword,"bitpix") == 0)
352 fits_info.bits_per_pixel=StringToLong(p);
353 if (LocaleCompare(keyword,"naxis") == 0)
354 fits_info.number_axes=StringToLong(p);
355 if (LocaleCompare(keyword,"naxis1") == 0)
356 fits_info.columns=StringToLong(p);
357 if (LocaleCompare(keyword,"naxis2") == 0)
358 fits_info.rows=StringToLong(p);
359 if (LocaleCompare(keyword,"naxis3") == 0)
360 fits_info.number_planes=StringToLong(p);
361 if (LocaleCompare(keyword,"datamax") == 0)
362 fits_info.max_data=StringToDouble(p,(char **) NULL);
363 if (LocaleCompare(keyword,"datamin") == 0)
364 fits_info.min_data=StringToDouble(p,(char **) NULL);
365 if (LocaleCompare(keyword,"bzero") == 0)
366 fits_info.zero=StringToDouble(p,(char **) NULL);
367 if (LocaleCompare(keyword,"bscale") == 0)
368 fits_info.scale=StringToDouble(p,(char **) NULL);
369 if (LocaleCompare(keyword,"comment") == 0)
371 if (comment == (char *) NULL)
372 comment=ConstantString(p);
374 (void) ConcatenateString(&comment,p);
376 if (LocaleCompare(keyword,"xendian") == 0)
378 if (LocaleNCompare(p,"big",3) == 0)
379 fits_info.endian=MSBEndian;
381 fits_info.endian=LSBEndian;
383 (void) FormatLocaleString(property,MaxTextExtent,"fits:%s",keyword);
384 (void) SetImageProperty(image,property,p,exception);
387 while (((TellBlob(image) % FITSBlocksize) != 0) && (c != EOF))
388 c=ReadBlobByte(image);
389 if (fits_info.extend == MagickFalse)
391 number_pixels=(MagickSizeType) fits_info.columns*fits_info.rows;
392 if ((fits_info.simple != MagickFalse) && (fits_info.number_axes >= 1) &&
393 (fits_info.number_axes <= 4) && (number_pixels != 0))
397 Verify that required image information is defined.
399 if (comment != (char *) NULL)
401 (void) SetImageProperty(image,"comment",comment,exception);
402 comment=DestroyString(comment);
404 if (EOFBlob(image) != MagickFalse)
405 ThrowFileException(exception,CorruptImageError,"UnexpectedEndOfFile",
407 number_pixels=(MagickSizeType) fits_info.columns*fits_info.rows;
408 if ((fits_info.simple == MagickFalse) || (fits_info.number_axes < 1) ||
409 (fits_info.number_axes > 4) || (number_pixels == 0))
410 ThrowReaderException(CorruptImageError,"ImageTypeNotSupported");
411 for (scene=0; scene < (ssize_t) fits_info.number_planes; scene++)
413 image->columns=(size_t) fits_info.columns;
414 image->rows=(size_t) fits_info.rows;
415 image->depth=(size_t) (fits_info.bits_per_pixel < 0 ? -1 : 1)*
416 fits_info.bits_per_pixel;
417 image->endian=fits_info.endian;
418 image->scene=(size_t) scene;
419 if ((image_info->ping != MagickFalse) && (image_info->number_scenes != 0))
420 if (image->scene >= (image_info->scene+image_info->number_scenes-1))
423 Initialize image structure.
425 if ((fits_info.min_data != 0.0) || (fits_info.max_data != 0.0))
427 if ((fits_info.bits_per_pixel != 0) && (fits_info.max_data == 0.0))
428 fits_info.max_data=GetFITSPixelRange((size_t)
429 fits_info.bits_per_pixel);
432 GetFITSPixelExtrema(image,fits_info.bits_per_pixel,&fits_info.min_data,
433 &fits_info.max_data);
435 Convert FITS pixels to pixel packets.
437 scale=(double) QuantumRange/(fits_info.scale*(fits_info.max_data-
438 fits_info.min_data)+fits_info.zero);
439 for (y=(ssize_t) image->rows-1; y >= 0; y--)
441 q=QueueAuthenticPixels(image,0,y,image->columns,1,exception);
442 if (q == (Quantum *) NULL)
444 for (x=0; x < (ssize_t) image->columns; x++)
446 pixel=GetFITSPixel(image,fits_info.bits_per_pixel);
447 SetPixelRed(image,ClampToQuantum(scale*(fits_info.scale*
448 (pixel-fits_info.min_data)+fits_info.zero)),q);
449 SetPixelGreen(image,GetPixelRed(image,q),q);
450 SetPixelBlue(image,GetPixelRed(image,q),q);
451 q+=GetPixelChannels(image);
453 if (SyncAuthenticPixels(image,exception) == MagickFalse)
455 if (image->previous == (Image *) NULL)
457 status=SetImageProgress(image,LoadImageTag,(MagickOffsetType) y,
459 if (status == MagickFalse)
463 if (EOFBlob(image) != MagickFalse)
465 ThrowFileException(exception,CorruptImageError,"UnexpectedEndOfFile",
470 Proceed to next image.
472 if (image_info->number_scenes != 0)
473 if (image->scene >= (image_info->scene+image_info->number_scenes-1))
475 if (scene < (ssize_t) (fits_info.number_planes-1))
478 Allocate next image structure.
480 AcquireNextImage(image_info,image,exception);
481 if (GetNextImageInList(image) == (Image *) NULL)
483 image=DestroyImageList(image);
484 return((Image *) NULL);
486 image=SyncNextImageInList(image);
487 status=SetImageProgress(image,LoadImagesTag,TellBlob(image),
489 if (status == MagickFalse)
493 (void) CloseBlob(image);
494 return(GetFirstImageInList(image));
498 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
502 % R e g i s t e r F I T S I m a g e %
506 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
508 % RegisterFITSImage() adds attributes for the FITS image format to
509 % the list of supported formats. The attributes include the image format
510 % tag, a method to read and/or write the format, whether the format
511 % supports the saving of more than one frame to the same file or blob,
512 % whether the format supports native in-memory I/O, and a brief
513 % description of the format.
515 % The format of the RegisterFITSImage method is:
517 % size_t RegisterFITSImage(void)
520 ModuleExport size_t RegisterFITSImage(void)
525 entry=SetMagickInfo("FITS");
526 entry->decoder=(DecodeImageHandler *) ReadFITSImage;
527 entry->encoder=(EncodeImageHandler *) WriteFITSImage;
528 entry->magick=(IsImageFormatHandler *) IsFITS;
529 entry->adjoin=MagickFalse;
530 entry->seekable_stream=MagickTrue;
531 entry->description=ConstantString("Flexible Image Transport System");
532 entry->module=ConstantString("FITS");
533 (void) RegisterMagickInfo(entry);
534 entry=SetMagickInfo("FTS");
535 entry->decoder=(DecodeImageHandler *) ReadFITSImage;
536 entry->encoder=(EncodeImageHandler *) WriteFITSImage;
537 entry->magick=(IsImageFormatHandler *) IsFITS;
538 entry->adjoin=MagickFalse;
539 entry->seekable_stream=MagickTrue;
540 entry->description=ConstantString("Flexible Image Transport System");
541 entry->module=ConstantString("FTS");
542 (void) RegisterMagickInfo(entry);
543 return(MagickImageCoderSignature);
547 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
551 % U n r e g i s t e r F I T S I m a g e %
555 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
557 % UnregisterFITSImage() removes format registrations made by the
558 % FITS module from the list of supported formats.
560 % The format of the UnregisterFITSImage method is:
562 % UnregisterFITSImage(void)
565 ModuleExport void UnregisterFITSImage(void)
567 (void) UnregisterMagickInfo("FITS");
568 (void) UnregisterMagickInfo("FTS");
572 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
576 % W r i t e F I T S I m a g e %
580 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
582 % WriteFITSImage() writes a Flexible Image Transport System image to a
583 % file as gray scale intensities [0..255].
585 % The format of the WriteFITSImage method is:
587 % MagickBooleanType WriteFITSImage(const ImageInfo *image_info,
588 % Image *image,ExceptionInfo *exception)
590 % A description of each parameter follows.
592 % o image_info: the image info.
594 % o image: The image.
596 % o exception: return any errors or warnings in this structure.
599 static MagickBooleanType WriteFITSImage(const ImageInfo *image_info,
600 Image *image,ExceptionInfo *exception)
603 header[FITSBlocksize],
612 register const Quantum
627 Open output image file.
629 assert(image_info != (const ImageInfo *) NULL);
630 assert(image_info->signature == MagickSignature);
631 assert(image != (Image *) NULL);
632 assert(image->signature == MagickSignature);
633 if (image->debug != MagickFalse)
634 (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
635 assert(exception != (ExceptionInfo *) NULL);
636 assert(exception->signature == MagickSignature);
637 status=OpenBlob(image_info,image,WriteBinaryBlobMode,exception);
638 if (status == MagickFalse)
640 if (IsRGBColorspace(image->colorspace) == MagickFalse)
641 (void) TransformImageColorspace(image,RGBColorspace,exception);
643 Allocate image memory.
645 fits_info=(char *) AcquireQuantumMemory(FITSBlocksize,sizeof(*fits_info));
646 if (fits_info == (char *) NULL)
647 ThrowWriterException(ResourceLimitError,"MemoryAllocationFailed");
648 (void) ResetMagickMemory(fits_info,' ',FITSBlocksize*sizeof(*fits_info));
650 Initialize image header.
652 image->depth=GetImageQuantumDepth(image,MagickFalse);
653 quantum_info=AcquireQuantumInfo((const ImageInfo *) NULL,image);
654 if (quantum_info == (QuantumInfo *) NULL)
655 ThrowWriterException(ResourceLimitError,"MemoryAllocationFailed");
657 (void) FormatLocaleString(header,FITSBlocksize,
659 (void) strncpy(fits_info+offset,header,strlen(header));
661 (void) FormatLocaleString(header,FITSBlocksize,"BITPIX = %10ld",
662 (long) (quantum_info->format == FloatingPointQuantumFormat ? -1 : 1)*
664 (void) strncpy(fits_info+offset,header,strlen(header));
666 (void) FormatLocaleString(header,FITSBlocksize,"NAXIS = %10lu",
667 IsImageGray(image,exception) != MagickFalse ? 2UL : 3UL);
668 (void) strncpy(fits_info+offset,header,strlen(header));
670 (void) FormatLocaleString(header,FITSBlocksize,"NAXIS1 = %10lu",
671 (unsigned long) image->columns);
672 (void) strncpy(fits_info+offset,header,strlen(header));
674 (void) FormatLocaleString(header,FITSBlocksize,"NAXIS2 = %10lu",
675 (unsigned long) image->rows);
676 (void) strncpy(fits_info+offset,header,strlen(header));
678 if (IsImageGray(image,exception) == MagickFalse)
680 (void) FormatLocaleString(header,FITSBlocksize,
681 "NAXIS3 = %10lu",3UL);
682 (void) strncpy(fits_info+offset,header,strlen(header));
685 (void) FormatLocaleString(header,FITSBlocksize,"BSCALE = %E",1.0);
686 (void) strncpy(fits_info+offset,header,strlen(header));
688 (void) FormatLocaleString(header,FITSBlocksize,"BZERO = %E",
689 image->depth > 8 ? GetFITSPixelRange(image->depth) : 0.0);
690 (void) strncpy(fits_info+offset,header,strlen(header));
692 (void) FormatLocaleString(header,FITSBlocksize,"DATAMAX = %E",
693 1.0*((MagickOffsetType) GetQuantumRange(image->depth)));
694 (void) strncpy(fits_info+offset,header,strlen(header));
696 (void) FormatLocaleString(header,FITSBlocksize,"DATAMIN = %E",0.0);
697 (void) strncpy(fits_info+offset,header,strlen(header));
699 if (image->endian == LSBEndian)
701 (void) FormatLocaleString(header,FITSBlocksize,"XENDIAN = 'SMALL'");
702 (void) strncpy(fits_info+offset,header,strlen(header));
705 (void) FormatLocaleString(header,FITSBlocksize,"HISTORY %.72s",
706 GetMagickVersion((size_t *) NULL));
707 (void) strncpy(fits_info+offset,header,strlen(header));
709 (void) strncpy(header,"END",FITSBlocksize);
710 (void) strncpy(fits_info+offset,header,strlen(header));
712 (void) WriteBlob(image,FITSBlocksize,(unsigned char *) fits_info);
714 Convert image to fits scale PseudoColor class.
716 pixels=GetQuantumPixels(quantum_info);
717 if (IsImageGray(image,exception) != MagickFalse)
719 length=GetQuantumExtent(image,quantum_info,GrayQuantum);
720 for (y=(ssize_t) image->rows-1; y >= 0; y--)
722 p=GetVirtualPixels(image,0,y,image->columns,1,exception);
723 if (p == (const Quantum *) NULL)
725 length=ExportQuantumPixels(image,(CacheView *) NULL,quantum_info,
726 GrayQuantum,pixels,exception);
727 if (image->depth == 16)
728 SetFITSUnsignedPixels(image->columns,image->depth,pixels);
729 if (((image->depth == 32) || (image->depth == 64)) &&
730 (quantum_info->format != FloatingPointQuantumFormat))
731 SetFITSUnsignedPixels(image->columns,image->depth,pixels);
732 count=WriteBlob(image,length,pixels);
733 if (count != (ssize_t) length)
735 status=SetImageProgress(image,SaveImageTag,(MagickOffsetType) y,
737 if (status == MagickFalse)
743 length=GetQuantumExtent(image,quantum_info,RedQuantum);
744 for (y=(ssize_t) image->rows-1; y >= 0; y--)
746 p=GetVirtualPixels(image,0,y,image->columns,1,exception);
747 if (p == (const Quantum *) NULL)
749 length=ExportQuantumPixels(image,(CacheView *) NULL,quantum_info,
750 RedQuantum,pixels,exception);
751 if (image->depth == 16)
752 SetFITSUnsignedPixels(image->columns,image->depth,pixels);
753 if (((image->depth == 32) || (image->depth == 64)) &&
754 (quantum_info->format != FloatingPointQuantumFormat))
755 SetFITSUnsignedPixels(image->columns,image->depth,pixels);
756 count=WriteBlob(image,length,pixels);
757 if (count != (ssize_t) length)
759 status=SetImageProgress(image,SaveImageTag,(MagickOffsetType) y,
761 if (status == MagickFalse)
764 length=GetQuantumExtent(image,quantum_info,GreenQuantum);
765 for (y=(ssize_t) image->rows-1; y >= 0; y--)
767 p=GetVirtualPixels(image,0,y,image->columns,1,exception);
768 if (p == (const Quantum *) NULL)
770 length=ExportQuantumPixels(image,(CacheView *) NULL,quantum_info,
771 GreenQuantum,pixels,exception);
772 if (image->depth == 16)
773 SetFITSUnsignedPixels(image->columns,image->depth,pixels);
774 if (((image->depth == 32) || (image->depth == 64)) &&
775 (quantum_info->format != FloatingPointQuantumFormat))
776 SetFITSUnsignedPixels(image->columns,image->depth,pixels);
777 count=WriteBlob(image,length,pixels);
778 if (count != (ssize_t) length)
780 status=SetImageProgress(image,SaveImageTag,(MagickOffsetType) y,
782 if (status == MagickFalse)
785 length=GetQuantumExtent(image,quantum_info,BlueQuantum);
786 for (y=(ssize_t) image->rows-1; y >= 0; y--)
788 p=GetVirtualPixels(image,0,y,image->columns,1,exception);
789 if (p == (const Quantum *) NULL)
791 length=ExportQuantumPixels(image,(CacheView *) NULL,quantum_info,
792 BlueQuantum,pixels,exception);
793 if (image->depth == 16)
794 SetFITSUnsignedPixels(image->columns,image->depth,pixels);
795 if (((image->depth == 32) || (image->depth == 64)) &&
796 (quantum_info->format != FloatingPointQuantumFormat))
797 SetFITSUnsignedPixels(image->columns,image->depth,pixels);
798 count=WriteBlob(image,length,pixels);
799 if (count != (ssize_t) length)
801 status=SetImageProgress(image,SaveImageTag,(MagickOffsetType) y,
803 if (status == MagickFalse)
807 quantum_info=DestroyQuantumInfo(quantum_info);
808 length=(size_t) (FITSBlocksize-TellBlob(image) % FITSBlocksize);
811 (void) ResetMagickMemory(fits_info,0,length*sizeof(*fits_info));
812 (void) WriteBlob(image,length,(unsigned char *) fits_info);
814 fits_info=DestroyString(fits_info);
815 (void) CloseBlob(image);