2 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
13 % Read/Write JPEG-2000 Image Format %
20 % Copyright 1999-2011 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/colorspace.h"
48 #include "MagickCore/colorspace-private.h"
49 #include "MagickCore/color.h"
50 #include "MagickCore/color-private.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/monitor.h"
59 #include "MagickCore/monitor-private.h"
60 #include "MagickCore/option.h"
61 #include "MagickCore/pixel-accessor.h"
62 #include "MagickCore/profile.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/module.h"
68 #if defined(MAGICKCORE_JP2_DELEGATE)
69 #ifndef JAS_IMAGE_CM_GRAY
70 #define JAS_IMAGE_CM_GRAY JAS_IMAGE_CS_GRAY
72 #ifndef JAS_IMAGE_CM_RGB
73 #define JAS_IMAGE_CM_RGB JAS_IMAGE_CS_RGB
76 #define uchar unsigned char
79 #define ushort unsigned short
82 #define uint unsigned int
84 #if !defined(ssize_tssize_t)
85 #define ssize_tssize_t long long
87 #if !defined(ussize_tssize_t)
88 #define ussize_tssize_t unsigned long long
93 #undef PACKAGE_TARNAME
94 #undef PACKAGE_VERSION
95 #include "jasper/jasper.h"
98 #undef PACKAGE_TARNAME
99 #undef PACKAGE_VERSION
104 Forward declarations.
106 #if defined(MAGICKCORE_JP2_DELEGATE)
107 static MagickBooleanType
108 WriteJP2Image(const ImageInfo *,Image *);
110 static volatile MagickBooleanType
111 instantiate_jp2 = MagickFalse;
115 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
123 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
125 % IsJP2() returns MagickTrue if the image format type, identified by the
126 % magick string, is JP2.
128 % The format of the IsJP2 method is:
130 % MagickBooleanType IsJP2(const unsigned char *magick,const size_t length)
132 % A description of each parameter follows:
134 % o magick: compare image format pattern against these bytes.
136 % o length: Specifies the length of the magick string.
139 static MagickBooleanType IsJP2(const unsigned char *magick,const size_t length)
143 if (memcmp(magick+4,"\152\120\040\040\015",5) == 0)
149 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
157 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
159 % IsJPC()() returns MagickTrue if the image format type, identified by the
160 % magick string, is JPC.
162 % The format of the IsJPC method is:
164 % MagickBooleanType IsJPC(const unsigned char *magick,const size_t length)
166 % A description of each parameter follows:
168 % o magick: compare image format pattern against these bytes.
170 % o length: Specifies the length of the magick string.
173 static MagickBooleanType IsJPC(const unsigned char *magick,const size_t length)
177 if (memcmp(magick,"\377\117",2) == 0)
183 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
187 % R e a d J P 2 I m a g e %
191 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
193 % ReadJP2Image() reads a JPEG 2000 Image file (JP2) or JPEG 2000
194 % codestream (JPC) image file and returns it. It allocates the memory
195 % necessary for the new Image structure and returns a pointer to the new
196 % image or set of images.
198 % JP2 support is originally written by Nathan Brown, nathanbrown@letu.edu.
200 % The format of the ReadJP2Image method is:
202 % Image *ReadJP2Image(const ImageInfo *image_info,
203 % ExceptionInfo *exception)
205 % A description of each parameter follows:
207 % o image_info: the image info.
209 % o exception: return any errors or warnings in this structure.
212 #if defined(MAGICKCORE_JP2_DELEGATE)
214 typedef struct _StreamManager
223 static int BlobRead(jas_stream_obj_t *object,char *buffer,const int length)
231 source=(StreamManager *) object;
232 count=ReadBlob(source->image,(size_t) length,(unsigned char *) buffer);
236 static int BlobWrite(jas_stream_obj_t *object,char *buffer,const int length)
244 source=(StreamManager *) object;
245 count=WriteBlob(source->image,(size_t) length,(unsigned char *) buffer);
249 static long BlobSeek(jas_stream_obj_t *object,long offset,int origin)
254 source=(StreamManager *) object;
255 return((long) SeekBlob(source->image,offset,origin));
258 static int BlobClose(jas_stream_obj_t *object)
263 source=(StreamManager *) object;
264 (void) CloseBlob(source->image);
266 source=(StreamManager *) NULL;
270 static inline size_t MagickMax(const size_t x,const size_t y)
277 static inline size_t MagickMin(const size_t x,const size_t y)
284 static jas_stream_t *JP2StreamManager(Image *image)
286 static jas_stream_ops_t
301 stream=(jas_stream_t *) jas_malloc(sizeof(*stream));
302 if (stream == (jas_stream_t *) NULL)
303 return((jas_stream_t *) NULL);
304 (void) ResetMagickMemory(stream,0,sizeof(*stream));
305 stream->rwlimit_=(-1);
306 stream->obj_=(jas_stream_obj_t *) jas_malloc(sizeof(StreamManager));
307 if (stream->obj_ == (jas_stream_obj_t *) NULL)
308 return((jas_stream_t *) NULL);
309 (void) ResetMagickMemory(stream->obj_,0,sizeof(StreamManager));
310 stream->ops_=(&StreamOperators);
311 stream->openmode_=JAS_STREAM_READ | JAS_STREAM_WRITE | JAS_STREAM_BINARY;
312 stream->bufbase_=(unsigned char *) jas_malloc(JAS_STREAM_BUFSIZE+
313 JAS_STREAM_MAXPUTBACK);
314 if (stream->bufbase_ == (void *) NULL)
316 stream->bufbase_=stream->tinybuf_;
321 stream->bufmode_=JAS_STREAM_FREEBUF | JAS_STREAM_BUFMODEMASK;
322 stream->bufsize_=JAS_STREAM_BUFSIZE;
324 stream->bufstart_=(&stream->bufbase_[JAS_STREAM_MAXPUTBACK]);
325 stream->ptr_=stream->bufstart_;
327 source=(StreamManager *) stream->obj_;
332 static Image *ReadJP2Image(const ImageInfo *image_info,ExceptionInfo *exception)
367 maximum_component_depth,
379 assert(image_info != (const ImageInfo *) NULL);
380 assert(image_info->signature == MagickSignature);
381 if (image_info->debug != MagickFalse)
382 (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",
383 image_info->filename);
384 assert(exception != (ExceptionInfo *) NULL);
385 assert(exception->signature == MagickSignature);
386 image=AcquireImage(image_info);
387 status=OpenBlob(image_info,image,ReadBinaryBlobMode,exception);
388 if (status == MagickFalse)
390 image=DestroyImageList(image);
391 return((Image *) NULL);
394 Initialize JPEG 2000 API.
396 jp2_stream=JP2StreamManager(image);
397 if (jp2_stream == (jas_stream_t *) NULL)
398 ThrowReaderException(DelegateError,"UnableToManageJP2Stream");
399 jp2_image=jas_image_decode(jp2_stream,-1,0);
400 if (jp2_image == (jas_image_t *) NULL)
402 (void) jas_stream_close(jp2_stream);
403 ThrowReaderException(DelegateError,"UnableToDecodeImageFile");
405 switch (jas_clrspc_fam(jas_image_clrspc(jp2_image)))
407 case JAS_CLRSPC_FAM_RGB:
409 components[0]=jas_image_getcmptbytype(jp2_image,JAS_IMAGE_CT_RGB_R);
410 components[1]=jas_image_getcmptbytype(jp2_image,JAS_IMAGE_CT_RGB_G);
411 components[2]=jas_image_getcmptbytype(jp2_image,JAS_IMAGE_CT_RGB_B);
412 if ((components[0] < 0) || (components[1] < 0) || (components[2] < 0))
414 (void) jas_stream_close(jp2_stream);
415 jas_image_destroy(jp2_image);
416 ThrowReaderException(CorruptImageError,"MissingImageChannel");
419 components[3]=jas_image_getcmptbytype(jp2_image,3);
420 if (components[3] > 0)
422 image->matte=MagickTrue;
427 case JAS_CLRSPC_FAM_GRAY:
429 components[0]=jas_image_getcmptbytype(jp2_image,JAS_IMAGE_CT_GRAY_Y);
430 if (components[0] < 0)
432 (void) jas_stream_close(jp2_stream);
433 jas_image_destroy(jp2_image);
434 ThrowReaderException(CorruptImageError,"MissingImageChannel");
439 case JAS_CLRSPC_FAM_YCBCR:
441 components[0]=jas_image_getcmptbytype(jp2_image,JAS_IMAGE_CT_YCBCR_Y);
442 components[1]=jas_image_getcmptbytype(jp2_image,JAS_IMAGE_CT_YCBCR_CB);
443 components[2]=jas_image_getcmptbytype(jp2_image,JAS_IMAGE_CT_YCBCR_CR);
444 if ((components[0] < 0) || (components[1] < 0) || (components[2] < 0))
446 (void) jas_stream_close(jp2_stream);
447 jas_image_destroy(jp2_image);
448 ThrowReaderException(CorruptImageError,"MissingImageChannel");
451 components[3]=jas_image_getcmptbytype(jp2_image,JAS_IMAGE_CT_UNKNOWN);
452 if (components[3] > 0)
454 image->matte=MagickTrue;
457 image->colorspace=YCbCrColorspace;
462 (void) jas_stream_close(jp2_stream);
463 jas_image_destroy(jp2_image);
464 ThrowReaderException(CoderError,"ColorspaceModelIsNotSupported");
467 image->columns=jas_image_width(jp2_image);
468 image->rows=jas_image_height(jp2_image);
469 image->compression=JPEG2000Compression;
470 for (i=0; i < (ssize_t) number_components; i++)
476 width=(size_t) (jas_image_cmptwidth(jp2_image,components[i])*
477 jas_image_cmpthstep(jp2_image,components[i]));
478 height=(size_t) (jas_image_cmptheight(jp2_image,components[i])*
479 jas_image_cmptvstep(jp2_image,components[i]));
480 x_step[i]=(unsigned int) jas_image_cmpthstep(jp2_image,components[i]);
481 y_step[i]=(unsigned int) jas_image_cmptvstep(jp2_image,components[i]);
482 if ((width != image->columns) || (height != image->rows) ||
483 (jas_image_cmpttlx(jp2_image,components[i]) != 0) ||
484 (jas_image_cmpttly(jp2_image,components[i]) != 0) ||
485 (x_step[i] != 1) || (y_step[i] != 1) ||
486 (jas_image_cmptsgnd(jp2_image,components[i]) != MagickFalse))
488 (void) jas_stream_close(jp2_stream);
489 jas_image_destroy(jp2_image);
490 ThrowReaderException(CoderError,"IrregularChannelGeometryNotSupported");
494 Convert JPEG 2000 pixels.
496 image->matte=number_components > 3 ? MagickTrue : MagickFalse;
497 maximum_component_depth=0;
498 for (i=0; i < (ssize_t) number_components; i++)
500 maximum_component_depth=(unsigned int) MagickMax((size_t)
501 jas_image_cmptprec(jp2_image,components[i]),(size_t)
502 maximum_component_depth);
503 pixels[i]=jas_matrix_create(1,(int) (image->columns/x_step[i]));
504 if (pixels[i] == (jas_matrix_t *) NULL)
506 for (--i; i >= 0; i--)
507 jas_matrix_destroy(pixels[i]);
508 jas_image_destroy(jp2_image);
509 ThrowReaderException(ResourceLimitError,"MemoryAllocationFailed");
512 image->depth=maximum_component_depth;
513 if (image_info->ping != MagickFalse)
515 (void) jas_stream_close(jp2_stream);
516 jas_image_destroy(jp2_image);
517 return(GetFirstImageInList(image));
519 for (i=0; i < (ssize_t) number_components; i++)
520 range[i]=GetQuantumRange((size_t) jas_image_cmptprec(jp2_image,
522 for (y=0; y < (ssize_t) image->rows; y++)
524 q=GetAuthenticPixels(image,0,y,image->columns,1,exception);
525 if (q == (const Quantum *) NULL)
527 for (i=0; i < (ssize_t) number_components; i++)
528 (void) jas_image_readcmpt(jp2_image,(short) components[i],0,
529 (jas_image_coord_t) (y/y_step[i]),(jas_image_coord_t) (image->columns/
530 x_step[i]),1,pixels[i]);
531 switch (number_components)
538 for (x=0; x < (ssize_t) image->columns; x++)
540 pixel=(QuantumAny) jas_matrix_getv(pixels[0],x/x_step[0]);
541 SetPixelRed(image,ScaleAnyToQuantum((QuantumAny) pixel,range[0]),q);
542 SetPixelGreen(image,GetPixelRed(image,q),q);
543 SetPixelBlue(image,GetPixelRed(image,q),q);
544 q+=GetPixelChannels(image);
553 for (x=0; x < (ssize_t) image->columns; x++)
555 pixel=(QuantumAny) jas_matrix_getv(pixels[0],x/x_step[0]);
556 SetPixelRed(image,ScaleAnyToQuantum((QuantumAny) pixel,range[0]),q);
557 pixel=(QuantumAny) jas_matrix_getv(pixels[1],x/x_step[1]);
558 SetPixelGreen(image,ScaleAnyToQuantum((QuantumAny) pixel,range[1]),q);
559 pixel=(QuantumAny) jas_matrix_getv(pixels[2],x/x_step[2]);
560 SetPixelBlue(image,ScaleAnyToQuantum((QuantumAny) pixel,range[2]),q);
561 q+=GetPixelChannels(image);
570 for (x=0; x < (ssize_t) image->columns; x++)
572 pixel=(QuantumAny) jas_matrix_getv(pixels[0],x/x_step[0]);
573 SetPixelRed(image,ScaleAnyToQuantum((QuantumAny) pixel,range[0]),q);
574 pixel=(QuantumAny) jas_matrix_getv(pixels[1],x/x_step[1]);
575 SetPixelGreen(image,ScaleAnyToQuantum((QuantumAny) pixel,range[1]),q);
576 pixel=(QuantumAny) jas_matrix_getv(pixels[2],x/x_step[2]);
577 SetPixelBlue(image,ScaleAnyToQuantum((QuantumAny) pixel,range[2]),q);
578 pixel=(QuantumAny) jas_matrix_getv(pixels[3],x/x_step[3]);
579 SetPixelAlpha(image,ScaleAnyToQuantum((QuantumAny) pixel,range[3]),q);
580 q+=GetPixelChannels(image);
585 if (SyncAuthenticPixels(image,exception) == MagickFalse)
587 status=SetImageProgress(image,LoadImageTag,(MagickOffsetType) y,
589 if (status == MagickFalse)
592 cm_profile=jas_image_cmprof(jp2_image);
593 icc_profile=(jas_iccprof_t *) NULL;
594 if (cm_profile != (jas_cmprof_t *) NULL)
595 icc_profile=jas_iccprof_createfromcmprof(cm_profile);
596 if (icc_profile != (jas_iccprof_t *) NULL)
601 icc_stream=jas_stream_memopen(NULL,0);
602 if ((icc_stream != (jas_stream_t *) NULL) &&
603 (jas_iccprof_save(icc_profile,icc_stream) == 0) &&
604 (jas_stream_flush(icc_stream) == 0))
614 Extract the icc profile, handle errors without much noise.
616 blob=(jas_stream_memobj_t *) icc_stream->obj_;
617 if (image->debug != MagickFalse)
618 (void) LogMagickEvent(CoderEvent,GetMagickModule(),
619 "Profile: ICC, %.20g bytes",(double) blob->len_);
620 profile=AcquireStringInfo(blob->len_);
621 SetStringInfoDatum(profile,blob->buf_);
622 icc_profile=(StringInfo *) GetImageProfile(image,"icc");
623 if (icc_profile == (StringInfo *) NULL)
624 (void) SetImageProfile(image,"icc",profile);
626 (void) ConcatenateStringInfo(icc_profile,profile);
627 profile=DestroyStringInfo(profile);
628 (void) jas_stream_close(icc_stream);
631 (void) jas_stream_close(jp2_stream);
632 jas_image_destroy(jp2_image);
633 for (i=0; i < (ssize_t) number_components; i++)
634 jas_matrix_destroy(pixels[i]);
635 return(GetFirstImageInList(image));
640 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
644 % R e g i s t e r J P 2 I m a g e %
648 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
650 % RegisterJP2Image() adds attributes for the JP2 image format to the list of
651 % supported formats. The attributes include the image format tag, a method
652 % method to read and/or write the format, whether the format supports the
653 % saving of more than one frame to the same file or blob, whether the format
654 % supports native in-memory I/O, and a brief description of the format.
656 % The format of the RegisterJP2Image method is:
658 % size_t RegisterJP2Image(void)
661 ModuleExport size_t RegisterJP2Image(void)
666 entry=SetMagickInfo("JP2");
667 entry->description=ConstantString("JPEG-2000 File Format Syntax");
668 entry->module=ConstantString("JP2");
669 entry->magick=(IsImageFormatHandler *) IsJP2;
670 entry->adjoin=MagickFalse;
671 entry->seekable_stream=MagickTrue;
672 entry->thread_support=NoThreadSupport;
673 #if defined(MAGICKCORE_JP2_DELEGATE)
674 entry->decoder=(DecodeImageHandler *) ReadJP2Image;
675 entry->encoder=(EncodeImageHandler *) WriteJP2Image;
677 (void) RegisterMagickInfo(entry);
678 entry=SetMagickInfo("JPC");
679 entry->description=ConstantString("JPEG-2000 Code Stream Syntax");
680 entry->module=ConstantString("JP2");
681 entry->magick=(IsImageFormatHandler *) IsJPC;
682 entry->adjoin=MagickFalse;
683 entry->seekable_stream=MagickTrue;
684 entry->thread_support=NoThreadSupport;
685 #if defined(MAGICKCORE_JP2_DELEGATE)
686 entry->decoder=(DecodeImageHandler *) ReadJP2Image;
687 entry->encoder=(EncodeImageHandler *) WriteJP2Image;
689 (void) RegisterMagickInfo(entry);
690 entry=SetMagickInfo("J2C");
691 entry->description=ConstantString("JPEG-2000 Code Stream Syntax");
692 entry->module=ConstantString("JP2");
693 entry->magick=(IsImageFormatHandler *) IsJPC;
694 entry->adjoin=MagickFalse;
695 entry->seekable_stream=MagickTrue;
696 entry->thread_support=NoThreadSupport;
697 #if defined(MAGICKCORE_JP2_DELEGATE)
698 entry->decoder=(DecodeImageHandler *) ReadJP2Image;
699 entry->encoder=(EncodeImageHandler *) WriteJP2Image;
701 (void) RegisterMagickInfo(entry);
702 entry=SetMagickInfo("JPX");
703 entry->description=ConstantString("JPEG-2000 File Format Syntax");
704 entry->module=ConstantString("JP2");
705 entry->magick=(IsImageFormatHandler *) IsJPC;
706 entry->adjoin=MagickFalse;
707 entry->seekable_stream=MagickTrue;
708 entry->thread_support=NoThreadSupport;
709 #if defined(MAGICKCORE_JP2_DELEGATE)
710 entry->decoder=(DecodeImageHandler *) ReadJP2Image;
711 entry->encoder=(EncodeImageHandler *) WriteJP2Image;
713 (void) RegisterMagickInfo(entry);
714 entry=SetMagickInfo("PGX");
715 entry->description=ConstantString("JPEG-2000 VM Format");
716 entry->module=ConstantString("JP2");
717 entry->magick=(IsImageFormatHandler *) IsJPC;
718 entry->adjoin=MagickFalse;
719 entry->seekable_stream=MagickTrue;
720 entry->thread_support=NoThreadSupport;
721 #if defined(MAGICKCORE_JP2_DELEGATE)
722 entry->decoder=(DecodeImageHandler *) ReadJP2Image;
724 (void) RegisterMagickInfo(entry);
725 #if defined(MAGICKCORE_JP2_DELEGATE)
726 if (instantiate_jp2 == MagickFalse)
729 instantiate_jp2=MagickTrue;
732 return(MagickImageCoderSignature);
736 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
740 % U n r e g i s t e r J P 2 I m a g e %
744 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
746 % UnregisterJP2Image() removes format registrations made by the JP2 module
747 % from the list of supported formats.
749 % The format of the UnregisterJP2Image method is:
751 % UnregisterJP2Image(void)
754 ModuleExport void UnregisterJP2Image(void)
756 (void) UnregisterMagickInfo("PGX");
757 (void) UnregisterMagickInfo("J2C");
758 (void) UnregisterMagickInfo("JPC");
759 (void) UnregisterMagickInfo("JP2");
760 #if defined(MAGICKCORE_JP2_DELEGATE)
761 if (instantiate_jp2 != MagickFalse)
764 instantiate_jp2=MagickFalse;
769 #if defined(MAGICKCORE_JP2_DELEGATE)
771 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
775 % W r i t e J P 2 I m a g e %
779 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
781 % WriteJP2Image() writes an image in the JPEG 2000 image format.
783 % JP2 support originally written by Nathan Brown, nathanbrown@letu.edu
785 % The format of the WriteJP2Image method is:
787 % MagickBooleanType WriteJP2Image(const ImageInfo *image_info,Image *image)
789 % A description of each parameter follows.
791 % o image_info: the image info.
793 % o image: The image.
796 static MagickBooleanType WriteJP2Image(const ImageInfo *image_info,Image *image)
800 magick[MaxTextExtent],
824 register const Quantum
841 assert(image_info != (const ImageInfo *) NULL);
842 assert(image_info->signature == MagickSignature);
843 assert(image != (Image *) NULL);
844 assert(image->signature == MagickSignature);
845 if (image->debug != MagickFalse)
846 (void) LogMagickEvent(TraceEvent,GetMagickModule(),"%s",image->filename);
847 status=OpenBlob(image_info,image,WriteBinaryBlobMode,&image->exception);
848 if (status == MagickFalse)
851 Initialize JPEG 2000 API.
853 if (IsRGBColorspace(image->colorspace) == MagickFalse)
854 (void) TransformImageColorspace(image,RGBColorspace);
855 jp2_stream=JP2StreamManager(image);
856 if (jp2_stream == (jas_stream_t *) NULL)
857 ThrowWriterException(DelegateError,"UnableToManageJP2Stream");
858 number_components=image->matte ? 4UL : 3UL;
859 if ((image_info->type != TrueColorType) &&
860 (IsImageGray(image,&image->exception) != MagickFalse))
862 if ((image->columns != (unsigned int) image->columns) ||
863 (image->rows != (unsigned int) image->rows))
864 ThrowWriterException(ImageError,"WidthOrHeightExceedsLimit");
865 (void) ResetMagickMemory(&component_info,0,sizeof(component_info));
866 for (i=0; i < (ssize_t) number_components; i++)
868 component_info[i].tlx=0;
869 component_info[i].tly=0;
870 component_info[i].hstep=1;
871 component_info[i].vstep=1;
872 component_info[i].width=(unsigned int) image->columns;
873 component_info[i].height=(unsigned int) image->rows;
874 component_info[i].prec=(int) MagickMax(MagickMin(image->depth,16),2);
875 component_info[i].sgnd=MagickFalse;
877 jp2_image=jas_image_create((int) number_components,component_info,
879 if (jp2_image == (jas_image_t *) NULL)
880 ThrowWriterException(DelegateError,"UnableToCreateImage");
881 if (number_components == 1)
886 jas_image_setclrspc(jp2_image,JAS_CLRSPC_SGRAY);
887 jas_image_setcmpttype(jp2_image,0,
888 JAS_IMAGE_CT_COLOR(JAS_CLRSPC_CHANIND_GRAY_Y));
895 jas_image_setclrspc(jp2_image,JAS_CLRSPC_SRGB);
896 jas_image_setcmpttype(jp2_image,0,
897 (jas_image_cmpttype_t) JAS_IMAGE_CT_COLOR(JAS_CLRSPC_CHANIND_RGB_R));
898 jas_image_setcmpttype(jp2_image,1,
899 (jas_image_cmpttype_t) JAS_IMAGE_CT_COLOR(JAS_CLRSPC_CHANIND_RGB_G));
900 jas_image_setcmpttype(jp2_image,2,
901 (jas_image_cmpttype_t) JAS_IMAGE_CT_COLOR(JAS_CLRSPC_CHANIND_RGB_B));
902 if (number_components == 4)
903 jas_image_setcmpttype(jp2_image,3,JAS_IMAGE_CT_OPACITY);
906 Convert to JPEG 2000 pixels.
908 for (i=0; i < (ssize_t) number_components; i++)
910 pixels[i]=jas_matrix_create(1,(int) image->columns);
911 if (pixels[i] == (jas_matrix_t *) NULL)
913 for (x=0; x < i; x++)
914 jas_matrix_destroy(pixels[x]);
915 jas_image_destroy(jp2_image);
916 ThrowWriterException(ResourceLimitError,"MemoryAllocationFailed");
919 range=GetQuantumRange((size_t) component_info[0].prec);
920 for (y=0; y < (ssize_t) image->rows; y++)
922 p=GetVirtualPixels(image,0,y,image->columns,1,&image->exception);
923 if (p == (const Quantum *) NULL)
925 for (x=0; x < (ssize_t) image->columns; x++)
927 if (number_components == 1)
928 jas_matrix_setv(pixels[0],x,(jas_seqent_t) ScaleQuantumToAny(
929 GetPixelIntensity(image,p),range));
932 jas_matrix_setv(pixels[0],x,(jas_seqent_t) ScaleQuantumToAny(
933 GetPixelRed(image,p),range));
934 jas_matrix_setv(pixels[1],x,(jas_seqent_t) ScaleQuantumToAny(
935 GetPixelGreen(image,p),range));
936 jas_matrix_setv(pixels[2],x,(jas_seqent_t) ScaleQuantumToAny(
937 GetPixelBlue(image,p),range));
938 if (number_components > 3)
939 jas_matrix_setv(pixels[3],x,(jas_seqent_t) ScaleQuantumToAny(
940 GetPixelAlpha(image,p),range));
942 p+=GetPixelChannels(image);
944 for (i=0; i < (ssize_t) number_components; i++)
945 (void) jas_image_writecmpt(jp2_image,(short) i,0,(unsigned int) y,
946 (unsigned int) image->columns,1,pixels[i]);
947 status=SetImageProgress(image,SaveImageTag,(MagickOffsetType) y,
949 if (status == MagickFalse)
952 (void) CopyMagickString(magick,image_info->magick,MaxTextExtent);
953 if (LocaleCompare(magick,"J2C") == 0)
954 (void) CopyMagickString(magick,"JPC",MaxTextExtent);
956 format=jas_image_strtofmt(magick);
957 options=(char *) NULL;
958 ResetImageOptionIterator(image_info);
959 key=GetNextImageOption(image_info);
960 for ( ; key != (char *) NULL; key=GetNextImageOption(image_info))
962 option=GetImageOption(image_info,key);
963 if (option == (const char *) NULL)
965 if (LocaleNCompare(key,"jp2:",4) == 0)
967 (void) ConcatenateString(&options,key+4);
970 (void) ConcatenateString(&options,"=");
971 (void) ConcatenateString(&options,option);
973 (void) ConcatenateString(&options," ");
976 option=GetImageOption(image_info,"jp2:rate");
977 if ((option == (const char *) NULL) &&
978 (image_info->compression != LosslessJPEGCompression) &&
979 (image->quality != UndefinedCompressionQuality) &&
980 ((double) image->quality <= 99.5) &&
981 ((image->rows*image->columns) > 2500))
984 option[MaxTextExtent];
993 alpha=115.0-image->quality;
994 rate=100.0/(alpha*alpha);
996 header_size+=(number_components-1)*142;
997 number_pixels=(double) image->rows*image->columns*number_components*
998 (GetImageQuantumDepth(image,MagickTrue)/8);
999 target_size=(number_pixels*rate)+header_size;
1000 rate=target_size/number_pixels;
1001 (void) FormatLocaleString(option,MaxTextExtent,"rate=%g",rate);
1002 (void) ConcatenateString(&options,option);
1004 status=jas_image_encode(jp2_image,jp2_stream,format,options) != 0 ?
1005 MagickTrue : MagickFalse;
1006 (void) jas_stream_close(jp2_stream);
1007 for (i=0; i < (ssize_t) number_components; i++)
1008 jas_matrix_destroy(pixels[i]);
1009 jas_image_destroy(jp2_image);
1010 if (status != MagickFalse)
1011 ThrowWriterException(DelegateError,"UnableToEncodeImageFile");