Home | History | Annotate | Download | only in lzw
      1 // Copyright 2011 The Go Authors. All rights reserved.
      2 // Use of this source code is governed by a BSD-style
      3 // license that can be found in the LICENSE file.
      4 
      5 // Package lzw implements the Lempel-Ziv-Welch compressed data format,
      6 // described in T. A. Welch, ``A Technique for High-Performance Data
      7 // Compression'', Computer, 17(6) (June 1984), pp 8-19.
      8 //
      9 // In particular, it implements LZW as used by the GIF and PDF file
     10 // formats, which means variable-width codes up to 12 bits and the first
     11 // two non-literal codes are a clear code and an EOF code.
     12 //
     13 // The TIFF file format uses a similar but incompatible version of the LZW
     14 // algorithm. See the golang.org/x/image/tiff/lzw package for an
     15 // implementation.
     16 package lzw
     17 
     18 // TODO(nigeltao): check that PDF uses LZW in the same way as GIF,
     19 // modulo LSB/MSB packing order.
     20 
     21 import (
     22 	"bufio"
     23 	"errors"
     24 	"fmt"
     25 	"io"
     26 )
     27 
     28 // Order specifies the bit ordering in an LZW data stream.
     29 type Order int
     30 
     31 const (
     32 	// LSB means Least Significant Bits first, as used in the GIF file format.
     33 	LSB Order = iota
     34 	// MSB means Most Significant Bits first, as used in the TIFF and PDF
     35 	// file formats.
     36 	MSB
     37 )
     38 
     39 const (
     40 	maxWidth           = 12
     41 	decoderInvalidCode = 0xffff
     42 	flushBuffer        = 1 << maxWidth
     43 )
     44 
     45 // decoder is the state from which the readXxx method converts a byte
     46 // stream into a code stream.
     47 type decoder struct {
     48 	r        io.ByteReader
     49 	bits     uint32
     50 	nBits    uint
     51 	width    uint
     52 	read     func(*decoder) (uint16, error) // readLSB or readMSB
     53 	litWidth int                            // width in bits of literal codes
     54 	err      error
     55 
     56 	// The first 1<<litWidth codes are literal codes.
     57 	// The next two codes mean clear and EOF.
     58 	// Other valid codes are in the range [lo, hi] where lo := clear + 2,
     59 	// with the upper bound incrementing on each code seen.
     60 	// overflow is the code at which hi overflows the code width.
     61 	// last is the most recently seen code, or decoderInvalidCode.
     62 	clear, eof, hi, overflow, last uint16
     63 
     64 	// Each code c in [lo, hi] expands to two or more bytes. For c != hi:
     65 	//   suffix[c] is the last of these bytes.
     66 	//   prefix[c] is the code for all but the last byte.
     67 	//   This code can either be a literal code or another code in [lo, c).
     68 	// The c == hi case is a special case.
     69 	suffix [1 << maxWidth]uint8
     70 	prefix [1 << maxWidth]uint16
     71 
     72 	// output is the temporary output buffer.
     73 	// Literal codes are accumulated from the start of the buffer.
     74 	// Non-literal codes decode to a sequence of suffixes that are first
     75 	// written right-to-left from the end of the buffer before being copied
     76 	// to the start of the buffer.
     77 	// It is flushed when it contains >= 1<<maxWidth bytes,
     78 	// so that there is always room to decode an entire code.
     79 	output [2 * 1 << maxWidth]byte
     80 	o      int    // write index into output
     81 	toRead []byte // bytes to return from Read
     82 }
     83 
     84 // readLSB returns the next code for "Least Significant Bits first" data.
     85 func (d *decoder) readLSB() (uint16, error) {
     86 	for d.nBits < d.width {
     87 		x, err := d.r.ReadByte()
     88 		if err != nil {
     89 			return 0, err
     90 		}
     91 		d.bits |= uint32(x) << d.nBits
     92 		d.nBits += 8
     93 	}
     94 	code := uint16(d.bits & (1<<d.width - 1))
     95 	d.bits >>= d.width
     96 	d.nBits -= d.width
     97 	return code, nil
     98 }
     99 
    100 // readMSB returns the next code for "Most Significant Bits first" data.
    101 func (d *decoder) readMSB() (uint16, error) {
    102 	for d.nBits < d.width {
    103 		x, err := d.r.ReadByte()
    104 		if err != nil {
    105 			return 0, err
    106 		}
    107 		d.bits |= uint32(x) << (24 - d.nBits)
    108 		d.nBits += 8
    109 	}
    110 	code := uint16(d.bits >> (32 - d.width))
    111 	d.bits <<= d.width
    112 	d.nBits -= d.width
    113 	return code, nil
    114 }
    115 
    116 func (d *decoder) Read(b []byte) (int, error) {
    117 	for {
    118 		if len(d.toRead) > 0 {
    119 			n := copy(b, d.toRead)
    120 			d.toRead = d.toRead[n:]
    121 			return n, nil
    122 		}
    123 		if d.err != nil {
    124 			return 0, d.err
    125 		}
    126 		d.decode()
    127 	}
    128 }
    129 
    130 // decode decompresses bytes from r and leaves them in d.toRead.
    131 // read specifies how to decode bytes into codes.
    132 // litWidth is the width in bits of literal codes.
    133 func (d *decoder) decode() {
    134 	// Loop over the code stream, converting codes into decompressed bytes.
    135 loop:
    136 	for {
    137 		code, err := d.read(d)
    138 		if err != nil {
    139 			if err == io.EOF {
    140 				err = io.ErrUnexpectedEOF
    141 			}
    142 			d.err = err
    143 			break
    144 		}
    145 		switch {
    146 		case code < d.clear:
    147 			// We have a literal code.
    148 			d.output[d.o] = uint8(code)
    149 			d.o++
    150 			if d.last != decoderInvalidCode {
    151 				// Save what the hi code expands to.
    152 				d.suffix[d.hi] = uint8(code)
    153 				d.prefix[d.hi] = d.last
    154 			}
    155 		case code == d.clear:
    156 			d.width = 1 + uint(d.litWidth)
    157 			d.hi = d.eof
    158 			d.overflow = 1 << d.width
    159 			d.last = decoderInvalidCode
    160 			continue
    161 		case code == d.eof:
    162 			d.err = io.EOF
    163 			break loop
    164 		case code <= d.hi:
    165 			c, i := code, len(d.output)-1
    166 			if code == d.hi {
    167 				// code == hi is a special case which expands to the last expansion
    168 				// followed by the head of the last expansion. To find the head, we walk
    169 				// the prefix chain until we find a literal code.
    170 				c = d.last
    171 				for c >= d.clear {
    172 					c = d.prefix[c]
    173 				}
    174 				d.output[i] = uint8(c)
    175 				i--
    176 				c = d.last
    177 			}
    178 			// Copy the suffix chain into output and then write that to w.
    179 			for c >= d.clear {
    180 				d.output[i] = d.suffix[c]
    181 				i--
    182 				c = d.prefix[c]
    183 			}
    184 			d.output[i] = uint8(c)
    185 			d.o += copy(d.output[d.o:], d.output[i:])
    186 			if d.last != decoderInvalidCode {
    187 				// Save what the hi code expands to.
    188 				d.suffix[d.hi] = uint8(c)
    189 				d.prefix[d.hi] = d.last
    190 			}
    191 		default:
    192 			d.err = errors.New("lzw: invalid code")
    193 			break loop
    194 		}
    195 		d.last, d.hi = code, d.hi+1
    196 		if d.hi >= d.overflow {
    197 			if d.width == maxWidth {
    198 				d.last = decoderInvalidCode
    199 			} else {
    200 				d.width++
    201 				d.overflow <<= 1
    202 			}
    203 		}
    204 		if d.o >= flushBuffer {
    205 			break
    206 		}
    207 	}
    208 	// Flush pending output.
    209 	d.toRead = d.output[:d.o]
    210 	d.o = 0
    211 }
    212 
    213 var errClosed = errors.New("lzw: reader/writer is closed")
    214 
    215 func (d *decoder) Close() error {
    216 	d.err = errClosed // in case any Reads come along
    217 	return nil
    218 }
    219 
    220 // NewReader creates a new io.ReadCloser.
    221 // Reads from the returned io.ReadCloser read and decompress data from r.
    222 // If r does not also implement io.ByteReader,
    223 // the decompressor may read more data than necessary from r.
    224 // It is the caller's responsibility to call Close on the ReadCloser when
    225 // finished reading.
    226 // The number of bits to use for literal codes, litWidth, must be in the
    227 // range [2,8] and is typically 8. It must equal the litWidth
    228 // used during compression.
    229 func NewReader(r io.Reader, order Order, litWidth int) io.ReadCloser {
    230 	d := new(decoder)
    231 	switch order {
    232 	case LSB:
    233 		d.read = (*decoder).readLSB
    234 	case MSB:
    235 		d.read = (*decoder).readMSB
    236 	default:
    237 		d.err = errors.New("lzw: unknown order")
    238 		return d
    239 	}
    240 	if litWidth < 2 || 8 < litWidth {
    241 		d.err = fmt.Errorf("lzw: litWidth %d out of range", litWidth)
    242 		return d
    243 	}
    244 	if br, ok := r.(io.ByteReader); ok {
    245 		d.r = br
    246 	} else {
    247 		d.r = bufio.NewReader(r)
    248 	}
    249 	d.litWidth = litWidth
    250 	d.width = 1 + uint(litWidth)
    251 	d.clear = uint16(1) << uint(litWidth)
    252 	d.eof, d.hi = d.clear+1, d.clear+1
    253 	d.overflow = uint16(1) << d.width
    254 	d.last = decoderInvalidCode
    255 
    256 	return d
    257 }
    258