X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=src%2Fgdcmjpegls%2FDecoder%2Fmelcode.c;fp=src%2Fgdcmjpegls%2FDecoder%2Fmelcode.c;h=82e5cbc045374e2a605193016a228c73b7adb5ea;hb=8a8b4ae83ad19690f0cc51ebc9091d75ffc1cd8d;hp=0000000000000000000000000000000000000000;hpb=34eec700805c1b9795aac2bfdc6fc84268517df9;p=gdcm.git

diff --git a/src/gdcmjpegls/Decoder/melcode.c b/src/gdcmjpegls/Decoder/melcode.c
new file mode 100644
index 00000000..82e5cbc0
--- /dev/null
+++ b/src/gdcmjpegls/Decoder/melcode.c
@@ -0,0 +1,156 @@
+/* SPMG/JPEG-LS IMPLEMENTATION V.2.1
+   =====================================
+   These programs are Copyright (c) University of British Columbia. All rights reserved.
+   They may be freely redistributed in their entirety provided that this copyright
+   notice is not removed.  They may not be sold for profit or incorporated in
+   commercial programs without the written permission of the copyright holder.
+   Each program is provided as is, without any express or implied warranty,
+   without even the warranty of fitness for a particular purpose.
+
+   =========================================================
+   THIS SOFTWARE IS BASED ON HP's implementation of jpeg-ls:
+   =========================================================
+
+   LOCO-I/JPEG-LS IMPLEMENTATION V.0.90
+   -------------------------------------------------------------------------------
+   (c) COPYRIGHT HEWLETT-PACKARD COMPANY, 1995-1997.
+        HEWLETT-PACKARD COMPANY ("HP") DOES NOT WARRANT THE ACCURACY OR
+   COMPLETENESS OF THE INFORMATION GIVEN HERE.  ANY USE MADE OF, OR
+   RELIANCE ON, SUCH INFORMATION IS ENTIRELY AT USER'S OWN RISK.
+        BY DOWNLOADING THE LOCO-I/JPEG-LS COMPRESSORS/DECOMPRESSORS
+   ("THE SOFTWARE") YOU AGREE TO BE BOUND BY THE TERMS AND CONDITIONS
+   OF THIS LICENSING AGREEMENT.
+        YOU MAY DOWNLOAD AND USE THE SOFTWARE FOR NON-COMMERCIAL PURPOSES
+   FREE OF CHARGE OR FURTHER OBLIGATION.  YOU MAY NOT, DIRECTLY OR
+   INDIRECTLY, DISTRIBUTE THE SOFTWARE FOR A FEE, INCORPORATE THIS
+   SOFTWARE INTO ANY PRODUCT OFFERED FOR SALE, OR USE THE SOFTWARE
+   TO PROVIDE A SERVICE FOR WHICH A FEE IS CHARGED.
+        YOU MAY MAKE COPIES OF THE SOFTWARE AND DISTRIBUTE SUCH COPIES TO
+   OTHER PERSONS PROVIDED THAT SUCH COPIES ARE ACCOMPANIED BY
+   HEWLETT-PACKARD'S COPYRIGHT NOTICE AND THIS AGREEMENT AND THAT
+   SUCH OTHER PERSONS AGREE TO BE BOUND BY THE TERMS OF THIS AGREEMENT.
+        THE SOFTWARE IS NOT OF PRODUCT QUALITY AND MAY HAVE ERRORS OR DEFECTS.
+   THE JPEG-LS STANDARD IS STILL UNDER DEVELOPMENT. THE SOFTWARE IS NOT A
+   FINAL OR FULL IMPLEMENTATION OF THE STANDARD.  HP GIVES NO EXPRESS OR
+   IMPLIED WARRANTY OF ANY KIND AND ANY IMPLIED WARRANTIES OF
+   MERCHANTABILITY AND FITNESS FOR PURPOSE ARE DISCLAIMED.
+        HP SHALL NOT BE LIABLE FOR ANY DIRECT, INDIRECT, SPECIAL, INCIDENTAL,
+   OR CONSEQUENTIAL DAMAGES ARISING OUT OF ANY USE OF THE SOFTWARE.
+   -------------------------------------------------------------------------------
+*/
+
+/* melcode.c --- for processing in run mode
+ *
+ * Initial code by Alex Jakulin,  Aug. 1995
+ *
+ * Modified and optimized: Gadiel Seroussi, October 1995
+ *
+ * Color Enhancement: Guillermo Sapiro, August 1996
+ *
+ * Modified and added Restart marker and input tables by:
+ * David Cheng-Hsiu Chu, and Ismail R. Ismail march 1999
+ */
+
+#include <stdio.h>
+#include "global.h"
+#include "bitio.h"
+
+#define MELCSTATES  32  /* number of melcode states */
+
+static J[MELCSTATES] = {
+  0,0,0,0,1,1,1,1,2,2,2,2,3,3,3,3,4,4,5,5,6,6,7,
+  7,8,9,10,11,12,13,14,15 
+};
+
+
+static int  melcstate[MAX_COMPONENTS],    /* index to the state array */
+    melclen[MAX_COMPONENTS];      /* contents of the state array location
+                    indexed by melcstate: the "expected"
+             run length is 2^melclen, shorter runs are
+             encoded by a 1 followed by the run length
+             in binary representation, wit a fixed length
+             of melclen bits */
+
+static unsigned long melcorder[MAX_COMPONENTS];  /* 2^ melclen */
+
+
+void init_process_run(int maxrun)    /* maxrun is ignoreed when using MELCODE,
+          kept here for function compatibility */            
+{
+  int  n_c;
+
+  for (n_c=0;n_c<components;n_c++)
+  {
+    melcstate[n_c] = 0;
+    melclen[n_c] = J[0];
+    melcorder[n_c] = 1<<melclen[n_c];
+  }
+}
+
+
+
+
+/* decoding routine: reads bits from the input and returns a run length. */
+/* argument is the number of pixels left to  end-of-line (bound on run length) */
+int process_run_dec(int lineleft, int color)  
+{
+  int runlen = 0;
+    
+  do {
+    register temp, hits;
+    temp = zeroLUT[(byte)(~(reg >> 24))];   /* number of leading ones in the
+                 input stream, up to 8 */
+    for ( hits = 1; hits<=temp; hits++ ) 
+    {
+      runlen += melcorder[color];
+      if ( runlen >= lineleft )
+      { /* reached end-of-line */
+        if ( runlen==lineleft && melcstate[color] < MELCSTATES ) 
+        {
+          melclen[color] = J[++melcstate[color]];
+          melcorder[color] = (1L<<melclen[color]);
+        }
+        FILLBUFFER(hits); /* actual # of 1's consumed */
+        return lineleft; 
+      }
+      if ( melcstate[color] < MELCSTATES ) 
+      {
+        melclen[color] = J[++melcstate[color]];
+        melcorder[color] = (1L<<melclen[color]);
+      }
+    }  
+    if (temp != 8) 
+    {
+      FILLBUFFER(temp + 1);  /* consume the leading
+              0 of the remainder encoding */
+      break;
+        }
+        FILLBUFFER(8);
+  } while ( 1 );
+
+  /* read the length of the remainder */
+  if ( melclen[color] ) 
+  {
+    register temp;
+    GETBITS(temp, melclen[color]);  /*** GETBITS is a macro, not a function */
+    runlen += temp;
+  }
+  limit_reduce = melclen[color]+1;
+
+  /* adjust melcoder parameters */
+  if ( melcstate[color] ) 
+  {
+    melclen[color] = J[--melcstate[color]];
+    melcorder[color] = (1L<<melclen[color]);
+  }
+
+  return runlen;
+}
+
+
+
+void
+close_process_run()
+{
+/* retained for compatibility with ranked runs */
+}