-
-#include <iostream>
-#include <sstream>
-
-#include <stdio.h>
-extern "C" {
-#include "gdcmjpeg/8/jconfig.h"
-#include "gdcmjpeg/8/jpeglib.h"
-#include "gdcmjpeg/8/jinclude.h"
-#include "gdcmjpeg/8/jerror.h"
-}
-
-#if defined(__BORLANDC__)
- #include <mem.h> // for memcpy
-#endif
-
-#include "gdcmJPEGFragment.h"
-#include <setjmp.h>
-#include <fstream>
-
-#include "jdatasrc.cxx"
-#include "jdatadst.cxx"
-
-typedef std::pair<size_t, uint32_t> JpegPair; //offset, jpeg size
-typedef std::vector<JpegPair> JpegVector;
-
-void WriteDICOMItems(std::ostream *fp, JpegVector &v)
-{
- // Item tag:
- uint16_t group = 0xfffe;
- uint16_t elem = 0xe000;
- gdcm::binary_write(*fp, group);
- gdcm::binary_write(*fp, elem);
- // Item Length
- uint32_t dummy = 0x12345678;
- size_t offset = fp->tellp();
- JpegPair jp;
- jp.first = offset;
- v.push_back(jp);
- gdcm::binary_write(*fp, dummy);
-}
-
-// PS 3.5, page 66
-void EncodeWithoutBasicOffsetTable(std::ostream *fp, int numFrag)// JpegVector& v) //, uint32_t length)
-{
- assert( numFrag == 1);
-
- // Item tag:
- uint16_t group = 0xfffe;
- uint16_t elem = 0xe000;
- gdcm::binary_write(*fp, group);
- gdcm::binary_write(*fp, elem);
- // Item Length
- uint32_t item_length = 0x0000;
- gdcm::binary_write(*fp, item_length);
-
-}
-
-// PS 3.5, page 67
-void EncodeWithBasicOffsetTable(std::ostream *fp, int numFrag, size_t &start)
-{
- // Item tag:
- uint16_t group = 0xfffe;
- uint16_t elem = 0xe000;
- gdcm::binary_write(*fp, group);
- gdcm::binary_write(*fp, elem);
- // Item Length
- uint32_t item_length = numFrag*4; // sizeof(uint32_t)
- gdcm::binary_write(*fp, item_length);
-
- // Just prepare the space
- start = fp->tellp(); //to be able to rewind
- for(int i=0; i<numFrag;++i)
- {
- uint32_t dummy = 0x0000;
- gdcm::binary_write(*fp, dummy);
- }
-}
-
-void UpdateBasicOffsetTable(std::ostream *fp, JpegVector const &v, size_t pos)
-{
- JpegVector::const_iterator i;
- fp->seekp( pos );
- const JpegPair &first = v[0];
- for(i=v.begin(); i!=v.end(); ++i)
- {
- const JpegPair &jp = *i;
- if(i == v.begin() ){ assert( jp.first - first.first == 0); }
- gdcm::binary_write(*fp, jp.first - first.first);
- std::cerr << "Updating Table:" << jp.first - first.first << std::endl;
- }
-}
-
-void UpdateJpegFragmentSize(std::ostream *fp, JpegVector const &v)
-{
- JpegVector::const_iterator i;
- for(i= v.begin(); i!=v.end(); ++i)
- {
- const JpegPair &jp = *i;
- fp->seekp( jp.first );
- gdcm::binary_write(*fp, jp.second );
- std::cerr << "Updating:" << jp.first << "," << jp.second << std::endl;
- }
-}
-
-void CloseJpeg(std::ostream *fp, JpegVector &v)
-{
- // sequence terminator
- uint16_t group = 0xfffe;
- uint16_t elem = 0xe000;
- gdcm::binary_write(*fp, group);
- gdcm::binary_write(*fp, elem);
-
- uint32_t length = 0x0;
- gdcm::binary_write(*fp, length);
-
- // Jpeg is done, now update the frag length
- UpdateJpegFragmentSize(fp, v);
-}
-
-bool InitializeJpeg(std::ostream *fp, int fragment_size, int image_width, int image_height,
- int sample_pixel, int quality, struct jpeg_compress_struct &cinfo, int &row_stride)
-{
-
- /* This struct contains the JPEG compression parameters and pointers to
- * working space (which is allocated as needed by the JPEG library).
- * It is possible to have several such structures, representing multiple
- * compression/decompression processes, in existence at once. We refer
- * to any one struct (and its associated working data) as a "JPEG object".
- */
- //struct jpeg_compress_struct cinfo;
- /* This struct represents a JPEG error handler. It is declared separately
- * because applications often want to supply a specialized error handler
- * (see the second half of this file for an example). But here we just
- * take the easy way out and use the standard error handler, which will
- * print a message on stderr and call exit() if compression fails.
- * Note that this struct must live as long as the main JPEG parameter
- * struct, to avoid dangling-pointer problems.
- */
- struct jpeg_error_mgr jerr;
- /* More stuff */
-
- /* Step 1: allocate and initialize JPEG compression object */
-
- /* We have to set up the error handler first, in case the initialization
- * step fails. (Unlikely, but it could happen if you are out of memory.)
- * This routine fills in the contents of struct jerr, and returns jerr's
- * address which we place into the link field in cinfo.
- */
- cinfo.err = jpeg_std_error(&jerr);
- /* Now we can initialize the JPEG compression object. */
- jpeg_create_compress(&cinfo);
-
- /* Step 2: specify data destination (eg, a file) */
- /* Note: steps 2 and 3 can be done in either order. */
-
- jpeg_stdio_dest(&cinfo, fp, fragment_size, 1);
-
- /* Step 3: set parameters for compression */
-
- /* First we supply a description of the input image.
- * Four fields of the cinfo struct must be filled in:
- */
- cinfo.image_width = image_width;/* image width and height, in pixels */
- cinfo.image_height = image_height;
- if ( sample_pixel == 3 )
- {
- cinfo.input_components = 3; /* # of color components per pixel */
- cinfo.in_color_space = JCS_RGB; /* colorspace of input image */
- }
- else
- {
- cinfo.input_components = 1; /* # of color components per pixel */
- cinfo.in_color_space = JCS_GRAYSCALE; /* colorspace of input image */
- }
- /* Now use the library's routine to set default compression parameters.
- * (You must set at least cinfo.in_color_space before calling this,
- * since the defaults depend on the source color space.)
- */
- jpeg_set_defaults(&cinfo);
- /*
- * http://www.koders.com/c/fid80DBBF1D49D004EF71CE7C493C34610C4F17D3D3.aspx
- * http://studio.imagemagick.org/pipermail/magick-users/2002-September/004685.html
- * You need to set -quality 101 or greater. If quality is 100 or less you
- * get regular JPEG output. This is not explained in the documentation, only
- * in the comments in coder/jpeg.c. When you have configured libjpeg with
- * lossless support, then
- *
- * quality=predictor*100 + point_transform
- *
- * If you don't know what these values should be, just use 101.
- * They only affect the compression ratio, not the image appearance,
- * which is lossless.
- */
- jpeg_simple_lossless (&cinfo, 1, 1);
- /* Now you can set any non-default parameters you wish to.
- * Here we just illustrate the use of quality (quantization table) scaling:
- */
- jpeg_set_quality(&cinfo, quality, TRUE /* limit to baseline-JPEG values */);
-
- /* Step 4: Start compressor */
-
- /* TRUE ensures that we will write a complete interchange-JPEG file.
- * Pass TRUE unless you are very sure of what you're doing.
- */
- jpeg_start_compress(&cinfo, TRUE);
-
- /* Step 5: while (scan lines remain to be written) */
- /* jpeg_write_scanlines(...); */
-
- /* Here we use the library's state variable cinfo.next_scanline as the
- * loop counter, so that we don't have to keep track ourselves.
- * To keep things simple, we pass one scanline per call; you can pass
- * more if you wish, though.
- */
- if (sample_pixel == 3)
- {
- row_stride = image_width * 3;/* JSAMPLEs per row in image_buffer */
- }
- else
- {
- row_stride = image_width * 1;/* JSAMPLEs per row in image_buffer */
- }
-
- /* everything was ok */
- return true;
-}
-
-bool FinalizeJpeg(struct jpeg_compress_struct &cinfo)
-{
- /* Step 6: Finish compression */
-
- jpeg_finish_compress(&cinfo);
-
- /* Step 7: release JPEG compression object */
-
- /* This is an important step since it will release a good deal of memory. */
- jpeg_destroy_compress(&cinfo);
-
- /* And we're done! */
- return true;
-}
-
-// If false then suspension return
-bool WriteScanlines(struct jpeg_compress_struct &cinfo, void *input_buffer, int row_stride)
-{
- JSAMPLE *image_buffer = (JSAMPLE*) input_buffer;
- JSAMPROW row_pointer[1]; /* pointer to JSAMPLE row[s] */
- row_pointer[0] = image_buffer;
-
- while (cinfo.next_scanline < cinfo.image_height) {
- /* jpeg_write_scanlines expects an array of pointers to scanlines.
- * Here the array is only one element long, but you could pass
- * more than one scanline at a time if that's more convenient.
- */
- //row_pointer[0] = & image_buffer[cinfo.next_scanline * row_stride];
-
- if( jpeg_write_scanlines(&cinfo, row_pointer, 1) != 1)
- {
- //entering suspension mode, basically we wrote the whole jpeg fragment
- // technically we could enforce that by checkig the value of row_pointer to
- // actually be at the end of the image...TODO
- return false;
- }
- row_pointer[0] += row_stride;
- }
-
- // Well looks like we are done writting the scanlines
- return true;
-}
-
-// input_buffer is ONE image
-// fragment_size is the size of this image (fragment)
-bool CreateOneFrame (std::ostream *fp, void *input_buffer, int fragment_size,
- int image_width, int image_height, int numZ, int sample_pixel, int quality, JpegVector &v)
-{
- struct jpeg_compress_struct cinfo;
- int row_stride; /* physical row width in image buffer */
- size_t beg = fp->tellp();
- bool r = InitializeJpeg(fp, fragment_size, image_width, image_height,
- sample_pixel, quality, cinfo, row_stride);
- assert( r );
- (void)numZ;
-
- uint8_t *pbuffer = (uint8_t*)input_buffer;
- //int i;
- //for(i=0; i<numZ; ++i)
-// {
- r = WriteScanlines(cinfo, pbuffer, row_stride);
- assert( r );
-// pbuffer+=fragment_size; //shift to next image
-
- //Upodate frag size
-// size_t end = fp->tellp();
-// std::cerr << "DIFF: " << end-beg << std::endl;
-
-// JpegPair &jp = v[i];
-// jp.second = end-beg;
- //beg = end; //
- // }
-
- r = FinalizeJpeg(cinfo);
- assert( r );
- size_t end = fp->tellp();
- static int i = 0;
- JpegPair &jp = v[i];
- jp.second = end-beg;
- std::cerr << "DIFF: " << i <<" -> " << end-beg << std::endl;
- ++i;
-
- //JpegPair &jp = v[0];
- //jp.second = 15328;
-
- return true;
-}
-
-//bool CreateMultipleFrames (std::ostream *fp, void *input_buffer, int fragment_size,
-// int image_width, int image_height, int sample_pixel, int quality, JpegVector &v)
-//{
-//}
-
-#define WITHOFFSETTABLE 1