6

I have a 12 bit packed image from a GigE camera. It is a little-endian file and each 3 bytes hold 2 12-bit pixels. I am trying to read this image using python and I tried something like this:

import bitstring
import numpy

with open('12bitpacked1.bin', 'rb') as f:
    data = f.read()
ii=numpy.zeros(2*len(data)/3)
ic = 0

for oo in range(0,len(data)/3):
    aa = bitstring.Bits(bytes=data[oo:oo+3], length=24)
    ii[ic],ii[ic+1] = aa.unpack('uint:12,uint:12')
    ic=ic+2

b = numpy.reshape(ii,(484,644))

In short: I read 3 bytes, convert them to bits and then unpack them as two 12-bit integers.

The result is, however, very different from what it should be. It looks like the image is separated into four quarters, each of them expanded to full image size and then overlapped.

What am I doing wrong here?

Update: Here are the test files:

12-bit packed

12-bit normal

They will not be identical, but they should show the same image. 12-bit normal has 12-bit pixel as uint16.

with open('12bit1.bin', 'rb') as f:
    a = numpy.fromfile(f, dtype=numpy.uint16)

b = numpy.reshape(a,(484,644))
Ivan
  • 65
  • 1
  • 4
  • Sounds like the image is interlaced, does it look correct if you combine the 4 quadrants into 1? See https://en.wikipedia.org/wiki/Interlacing_(bitmaps) – Dana the Sane Jun 13 '16 at 19:41
  • No, it is not interlaced. If it is it would at least horizontally look fine, but this is wrong in both directions. – Ivan Jun 13 '16 at 19:45

1 Answers1

6

With this code

for oo in range(0,len(data)/3):
  aa = bitstring.Bits(bytes=data[oo:oo+3], length=24)

you are reading bytes data[0:3], data[1:4], ... What you actually want is probably this:

for oo in range(0,len(data)/3):
  aa = bitstring.Bits(bytes=data[3*oo:3*oo+3], length=24)

[EDIT] Even more compact would be this:

for oo in range(0,len(data)-2,3):
  aa = bitstring.Bits(bytes=data[oo:oo+3], length=24)
coproc
  • 6,027
  • 2
  • 20
  • 31
  • Thanks, that was it. I had it first correct and then was testing something and changed it to this and forgot to put it back. – Ivan Jun 13 '16 at 20:30