20 %  Copyright 1999-2019 ImageMagick Studio LLC, a non-profit organization      %
23 %  You may not use this file except in compliance with the License.  You may  %
45 #include "MagickCore/blob-private.h"
47 #include "MagickCore/color-private.h"
49 #include "MagickCore/colorspace-private.h"
52 #include "MagickCore/exception-private.h"
54 #include "MagickCore/image-private.h"
60 #include "MagickCore/monitor-private.h"
61 #include "MagickCore/pixel-accessor.h"
63 #include "MagickCore/quantum-private.h"
67 #include "MagickCore/string-private.h"
147   switch (image->depth >> 3)  in GetFITSPixel()
186   if (offset == -1)  in GetFITSPixelExtrema()
187     return(-1);  in GetFITSPixelExtrema()
188   number_pixels=(MagickSizeType) image->columns*image->rows;  in GetFITSPixelExtrema()
202 static inline double GetFITSPixelRange(const size_t depth)  in GetFITSPixelRange()  argument
204   return((double) ((MagickOffsetType) GetQuantumRange(depth)));  in GetFITSPixelRange()
214     pixels+=(bits_per_pixel >> 3)-1;  in SetFITSUnsignedPixels()
277   register Quantum  in ReadFITSImage()
289   assert(image_info->signature == MagickCoreSignature);  in ReadFITSImage()
290   if (image_info->debug != MagickFalse)  in ReadFITSImage()
292       image_info->filename);  in ReadFITSImage()
294   assert(exception->signature == MagickCoreSignature);  in ReadFITSImage()
396         (fits_info.bits_per_pixel != -32) && (fits_info.bits_per_pixel != -64))  in ReadFITSImage()
417       image->filename);  in ReadFITSImage()
424     image->columns=(size_t) fits_info.columns;  in ReadFITSImage()
425     image->rows=(size_t) fits_info.rows;  in ReadFITSImage()
426     image->depth=(size_t) (fits_info.bits_per_pixel < 0 ? -1 : 1)*  in ReadFITSImage()
428     image->endian=fits_info.endian;  in ReadFITSImage()
429     image->scene=(size_t) scene;  in ReadFITSImage()
430     if ((image_info->ping != MagickFalse) && (image_info->number_scenes != 0))  in ReadFITSImage()
431       if (image->scene >= (image_info->scene+image_info->number_scenes-1))  in ReadFITSImage()
433     status=SetImageExtent(image,image->columns,image->rows,exception);  in ReadFITSImage()
442         if ((fits_info.bits_per_pixel == -32) ||  in ReadFITSImage()
443             (fits_info.bits_per_pixel == -64))  in ReadFITSImage()
455     scale=QuantumRange*PerceptibleReciprocal(fits_info.max_data-  in ReadFITSImage()
457     for (y=(ssize_t) image->rows-1; y >= 0; y--)  in ReadFITSImage()
459       q=QueueAuthenticPixels(image,0,y,image->columns,1,exception);  in ReadFITSImage()
460       if (q == (Quantum *) NULL)  in ReadFITSImage()
462       for (x=0; x < (ssize_t) image->columns; x++)  in ReadFITSImage()
465         if ((image->depth == 16) || (image->depth == 32) ||  in ReadFITSImage()
466             (image->depth == 64))  in ReadFITSImage()
467           SetFITSUnsignedPixels(1,image->depth,image->endian,  in ReadFITSImage()
469         SetPixelGray(image,ClampToQuantum(scale*(fits_info.scale*(pixel-  in ReadFITSImage()
475       if (image->previous == (Image *) NULL)  in ReadFITSImage()
478             image->rows);  in ReadFITSImage()
486           image->filename);  in ReadFITSImage()
492     if (image_info->number_scenes != 0)  in ReadFITSImage()
493       if (image->scene >= (image_info->scene+image_info->number_scenes-1))  in ReadFITSImage()
495     if (scene < (ssize_t) (fits_info.number_planes-1))  in ReadFITSImage()
534 %  whether the format supports native in-memory I/O, and a brief
548   entry->decoder=(DecodeImageHandler *) ReadFITSImage;  in RegisterFITSImage()
549   entry->encoder=(EncodeImageHandler *) WriteFITSImage;  in RegisterFITSImage()
550   entry->magick=(IsImageFormatHandler *) IsFITS;  in RegisterFITSImage()
551   entry->flags^=CoderAdjoinFlag;  in RegisterFITSImage()
552   entry->flags|=CoderDecoderSeekableStreamFlag;  in RegisterFITSImage()
555   entry->decoder=(DecodeImageHandler *) ReadFITSImage;  in RegisterFITSImage()
556   entry->encoder=(EncodeImageHandler *) WriteFITSImage;  in RegisterFITSImage()
557   entry->magick=(IsImageFormatHandler *) IsFITS;  in RegisterFITSImage()
558   entry->flags^=CoderAdjoinFlag;  in RegisterFITSImage()
559   entry->flags|=CoderDecoderSeekableStreamFlag;  in RegisterFITSImage()
627   if (length > (size_t) (FITSBlocksize-offset))  in CopyFitsRecord()
628     length=FITSBlocksize-offset;  in CopyFitsRecord()
646   register const Quantum  in WriteFITSImage()
664   assert(image_info->signature == MagickCoreSignature);  in WriteFITSImage()
666   assert(image->signature == MagickCoreSignature);  in WriteFITSImage()
667   if (image->debug != MagickFalse)  in WriteFITSImage()
668     (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);  in WriteFITSImage()
670   assert(exception->signature == MagickCoreSignature);  in WriteFITSImage()
685   image->depth=GetImageQuantumDepth(image,MagickFalse);  in WriteFITSImage()
686   image->endian=MSBEndian;  in WriteFITSImage()
699     (long) ((quantum_info->format == FloatingPointQuantumFormat ? -1 : 1)*  in WriteFITSImage()
700     image->depth));  in WriteFITSImage()
708     (unsigned long) image->columns);  in WriteFITSImage()
712     (unsigned long) image->rows);  in WriteFITSImage()
726     image->depth > 8 ? (GetFITSPixelRange(image->depth)+1)/2.0 : 0.0);  in WriteFITSImage()
730     1.0*((MagickOffsetType) GetQuantumRange(image->depth)));  in WriteFITSImage()
736   if (image->endian == LSBEndian)  in WriteFITSImage()
758       for (y=(ssize_t) image->rows-1; y >= 0; y--)  in WriteFITSImage()
760         p=GetVirtualPixels(image,0,y,image->columns,1,exception);  in WriteFITSImage()
761         if (p == (const Quantum *) NULL)  in WriteFITSImage()
765         if (image->depth == 16)  in WriteFITSImage()
766           SetFITSUnsignedPixels(image->columns,image->depth,image->endian,  in WriteFITSImage()
768         if (((image->depth == 32) || (image->depth == 64)) &&  in WriteFITSImage()
769             (quantum_info->format != FloatingPointQuantumFormat))  in WriteFITSImage()
770           SetFITSUnsignedPixels(image->columns,image->depth,image->endian,  in WriteFITSImage()
776           image->rows);  in WriteFITSImage()
784       for (y=(ssize_t) image->rows-1; y >= 0; y--)  in WriteFITSImage()
786         p=GetVirtualPixels(image,0,y,image->columns,1,exception);  in WriteFITSImage()
787         if (p == (const Quantum *) NULL)  in WriteFITSImage()
791         if (image->depth == 16)  in WriteFITSImage()
792           SetFITSUnsignedPixels(image->columns,image->depth,image->endian,  in WriteFITSImage()
794         if (((image->depth == 32) || (image->depth == 64)) &&  in WriteFITSImage()
795             (quantum_info->format != FloatingPointQuantumFormat))  in WriteFITSImage()
796           SetFITSUnsignedPixels(image->columns,image->depth,image->endian,  in WriteFITSImage()
802           image->rows);  in WriteFITSImage()
807       for (y=(ssize_t) image->rows-1; y >= 0; y--)  in WriteFITSImage()
809         p=GetVirtualPixels(image,0,y,image->columns,1,exception);  in WriteFITSImage()
810         if (p == (const Quantum *) NULL)  in WriteFITSImage()
814         if (image->depth == 16)  in WriteFITSImage()
815           SetFITSUnsignedPixels(image->columns,image->depth,image->endian,  in WriteFITSImage()
817         if (((image->depth == 32) || (image->depth == 64)) &&  in WriteFITSImage()
818             (quantum_info->format != FloatingPointQuantumFormat))  in WriteFITSImage()
819           SetFITSUnsignedPixels(image->columns,image->depth,image->endian,  in WriteFITSImage()
825           image->rows);  in WriteFITSImage()
830       for (y=(ssize_t) image->rows-1; y >= 0; y--)  in WriteFITSImage()
832         p=GetVirtualPixels(image,0,y,image->columns,1,exception);  in WriteFITSImage()
833         if (p == (const Quantum *) NULL)  in WriteFITSImage()
837         if (image->depth == 16)  in WriteFITSImage()
838           SetFITSUnsignedPixels(image->columns,image->depth,image->endian,  in WriteFITSImage()
840         if (((image->depth == 32) || (image->depth == 64)) &&  in WriteFITSImage()
841             (quantum_info->format != FloatingPointQuantumFormat))  in WriteFITSImage()
842           SetFITSUnsignedPixels(image->columns,image->depth,image->endian,  in WriteFITSImage()
848           image->rows);  in WriteFITSImage()
854   length=(size_t) (FITSBlocksize-TellBlob(image) % FITSBlocksize);  in WriteFITSImage()