Velocity Reviews > average of PIL images

# average of PIL images

vaneric
Guest
Posts: n/a

 02-18-2008
hi
i have a set of RGB images of diff faces (of people )as a 2 dim
numpyarray
...something like
threefaces=array([[xa1,xa2,xa3],
[xb1,xb2,xb3],
[xc1,xc2,xc3]])
where xa1,xa2,xa3 are tuples each representing rgb values of a pixel
of first image ..

i need to create the average face image and display it.problem is i
need to calculate (xa1+xb1+xc1)/3 etc to calculate avearge value of
each pixel.how can i do this calculation.do i need to convert the
r,g,b in to a single value for each pixel? the average value of a
pixel will be a float isn't it? how can i create a PIL image from
this?
any help,directive greatly appreciated
eric

Robert Kern
Guest
Posts: n/a

 02-18-2008
vaneric wrote:
> hi
> i have a set of RGB images of diff faces (of people )as a 2 dim
> numpyarray
> ..something like
> threefaces=array([[xa1,xa2,xa3],
> [xb1,xb2,xb3],
> [xc1,xc2,xc3]])
> where xa1,xa2,xa3 are tuples each representing rgb values of a pixel
> of first image ..
>
> i need to create the average face image and display it.problem is i
> need to calculate (xa1+xb1+xc1)/3 etc to calculate avearge value of
> each pixel.how can i do this calculation.

threefaces.mean(axis=0)

> do i need to convert the
> r,g,b in to a single value for each pixel?

It depends on the problem that you are trying to solve. If monochrome images are
acceptable for your problem, then it is probably best to convert all your images
to monochrome to do the averaging. At least for a first cut. Averaging color
images is tricky; you really shouldn't do it in the RGB colorspace. There are a
couple of colorspaces which you could choose; different problems require
different colorspaces.

> the average value of a
> pixel will be a float isn't it?

Yes.

> how can i create a PIL image from
> this?

# Cast the floating point data to bytes.
avgface = avgface.astype(numpy.uint

# Create the PIL image from the numpy data.
img = Image.fromstring('L', (width, height), avgface.tostring())

--
Robert Kern

"I have come to believe that the whole world is an enigma, a harmless enigma
an underlying truth."
-- Umberto Eco

Larry Bates
Guest
Posts: n/a

 02-18-2008
vaneric wrote:
> hi
> i have a set of RGB images of diff faces (of people )as a 2 dim
> numpyarray
> ..something like
> threefaces=array([[xa1,xa2,xa3],
> [xb1,xb2,xb3],
> [xc1,xc2,xc3]])
> where xa1,xa2,xa3 are tuples each representing rgb values of a pixel
> of first image ..
>
> i need to create the average face image and display it.problem is i
> need to calculate (xa1+xb1+xc1)/3 etc to calculate avearge value of
> each pixel.how can i do this calculation.do i need to convert the
> r,g,b in to a single value for each pixel? the average value of a
> pixel will be a float isn't it? how can i create a PIL image from
> this?
> any help,directive greatly appreciated
> eric

Take a look at ImageChops.difference. I've used it to calculate a
difference value as follows:

diff=ImageChops.difference(im1, im2)
totaldiff=sum(ImageStat.Stat(diff)._getmedian())

Maybe at least this will point you in the right direction.

-Larry Bates

7stud
Guest
Posts: n/a

 02-18-2008
On Feb 18, 10:18*am, vaneric <(E-Mail Removed)> wrote:
> hi
> i have a set of RGB images of diff faces (of people )as a 2 dim
> numpyarray
> ..something like
> threefaces=array([[xa1,xa2,xa3],
> * * * *[xb1,xb2,xb3],
> * * * *[xc1,xc2,xc3]])
> where xa1,xa2,xa3 are *tuples each representing rgb values of a pixel
> of first image ..
>
> i need to create the average face image and display it.problem is i
> need to calculate (xa1+xb1+xc1)/3 *etc to calculate avearge value of
> each pixel.how can i do this calculation.do i need to convert the
> r,g,b in to a single value for each pixel? the average value of a
> pixel will be a float isn't it? how can i create a PIL image from
> this?
> any help,directive greatly appreciated
> eric

import Numeric

arr = Numeric.array([
[1, 2, 3],
[4, 5, 6],
[7, 8, 9],
[10, 11, 12]
])

col1 = arr[0:,0]
print col1

sum = 0
count = 0
for num in col1:
sum += num
count += 1

avg = (sum * 1.0)/count #convert one term to a float to get float
result
print avg
print round(avg)
print int(round(avg))

print arr[0]

size = len(arr[0])
for i in range(size):
col = arr[0:, i]
sum = 0
count = 0

for num in col:
sum += num
count += 1

result = (sum * 1.0) /count
print result,
result = int(round(result))
print result

--output:--
[ 1 4 7 10]
5.5
6.0
6
[1 2 3]
5.5 6
6.5 7
7.5 8

7stud
Guest
Posts: n/a

 02-18-2008
On Feb 18, 1:58*pm, 7stud <(E-Mail Removed)> wrote:
> On Feb 18, 10:18*am, vaneric <(E-Mail Removed)> wrote:
>
>
>
> > hi
> > i have a set of RGB images of diff faces (of people )as a 2 dim
> > numpyarray
> > ..something like
> > threefaces=array([[xa1,xa2,xa3],
> > * * * *[xb1,xb2,xb3],
> > * * * *[xc1,xc2,xc3]])
> > where xa1,xa2,xa3 are *tuples each representing rgb values of a pixel
> > of first image ..

>
> > i need to create the average face image and display it.problem is i
> > need to calculate (xa1+xb1+xc1)/3 *etc to calculate avearge value of
> > each pixel.how can i do this calculation.do i need to convert the
> > r,g,b in to a single value for each pixel? the average value of a
> > pixel will be a float isn't it? how can i create a PIL image from
> > this?
> > any help,directive greatly appreciated
> > eric

>
> import Numeric
>
> arr = Numeric.array([
> * * [1, 2, 3],
> * * [4, 5, 6],
> * * [7, 8, 9],
> * * [10, 11, 12]
> ])
>
> col1 = arr[0:,0]
> print col1
>
> sum = 0
> count = 0
> for num in col1:
> * * sum += num
> * * count += 1
>
> avg = (sum * 1.0)/count * #convert one term to a float to get float
> result
> print avg
> print round(avg)
> print int(round(avg))
>
> print arr[0]
>
> size = len(arr[0])
> for i in range(size):
> * * col = arr[0:, i]
> * * sum = 0
> * * count = 0
>
> * * for num in col:
> * * * * sum += num
> * * * * count += 1
>
> * * result = (sum * 1.0) /count
> * * print result,
> * * result = int(round(result))
> * * print result
>
> --output:--
> [ 1 *4 *7 10]
> 5.5
> 6.0
> 6
> [1 2 3]
> 5.5 6
> 6.5 7
> 7.5 8

In this statement:

> col1 = arr[0:,0]

the first term is the row or row range, and the second term is the
column or column range. If you write this:

num = arr[0,0]

you get the element in row 0, column 0. But you can also specify
ranges for each col or row:

num = arr[1:, 2:]

That says to get all elements from row 1 to the bottom that are in
from column 2 to the end of the row.

7stud
Guest
Posts: n/a

 02-18-2008
On Feb 18, 2:05*pm, 7stud <(E-Mail Removed)> wrote:
> num = arr[1:, 2:]
>
> That says to get all elements from row 1 to the bottom that are in
> from column 2 to the end of the row.

err..

That says to get all elements from row 1 to the last row which are in
column 2, column 3, etc. to the end of the row.

vaneric
Guest
Posts: n/a

 02-19-2008
On Feb 19, 1:38 am, Robert Kern <(E-Mail Removed)> wrote:
>Averaging color
> images is tricky; you really shouldn't do it in the RGB colorspace.

hi,
thanx for the guidance and detailed replies..I tried to pack the
r,g,b into a single value like below(something a member posted in the
past)

def rgbTopixelvalue((r,g,b)):
alpha=255
return unpack("l", pack("BBBB", b, g, r, alpha))[0]

if i apply this for all images i can represent each image as an array
of longs instead of tuples.then for a pixel i can calculate the
average value
after this if i do the casting as you advised and create the avg
image
avgface = avgface.astype(numpy.uint
img = Image.fromstring('L', (width, height), avgface.tostring())

is there something wrong with my approach? I am a newbie in PIL/
imageprocessing ..so i would greately appreciate feedback
eric

vaneric
Guest
Posts: n/a

 02-19-2008
On Feb 19, 1:38 am, Robert Kern <(E-Mail Removed)> wrote:
>Averaging color
> images is tricky; you really shouldn't do it in the RGB colorspace.

hi,
thanx for the guidance and detailed replies..I tried to pack the
r,g,b into a single value like below(something a member posted in the
past)

def rgbTopixelvalue((r,g,b)):
alpha=255
return unpack("l", pack("BBBB", b, g, r, alpha))[0]

if i apply this for all images i can represent each image as an array
of longs instead of tuples.then for a pixel i can calculate the
average value
after this if i do the casting as you advised and create the avg
image
avgface = avgface.astype(numpy.uint

here if i use these pixelvalues to create an imag
img =Image.fromstring('RGB', (width, height), avgface.tostring())
it will fail because of -'not enough image data'..is there an
alternative to create average rgb color image ?(i want to keep the rgb
so can't convert to greyscale)

is there something wrong with my approach? I am a newbie in PIL/
imageprocessing ..so i would greately appreciate feedback
eric

Gabriel Genellina
Guest
Posts: n/a

 02-19-2008
En Tue, 19 Feb 2008 04:01:04 -0200, vaneric <(E-Mail Removed)>
escribió:
> On Feb 19, 1:38 am, Robert Kern <(E-Mail Removed)> wrote:

>> Averaging color
>> images is tricky; you really shouldn't do it in the RGB colorspace.

>
> hi,
> thanx for the guidance and detailed replies..I tried to pack the
> r,g,b into a single value like below(something a member posted in the
> past)
>
> def rgbTopixelvalue((r,g,b)):
> alpha=255
> return unpack("l", pack("BBBB", b, g, r, alpha))[0]

That's much worse than averaging the R,G,B components. First, you have to
omit the alfa value (or set it at the end). Then, consider this example:
(0,0,0)=black and (0,1,0)=almost black, average = (0,0,12 = dark blue, a
total nonsense.

As said above, try to compute using another color space, try HSL. The
colorsys module can transform from/to RGB.

--
Gabriel Genellina

vaneric
Guest
Posts: n/a

 02-19-2008

> > def rgbTopixelvalue((r,g,b)):
> > alpha=255
> > return unpack("l", pack("BBBB", b, g, r, alpha))[0]

>
> That's much worse than averaging the R,G,B components.

oops!
the intention was to pack r,g,b components into a single value sothat
calculations like finding covariant matrix of a set of images etc can
be done..(i need to represent each image as an array of values(either
long or float)..i can't work with an array of tuples of ints..

> As said above, try to compute using another color space, try HSL. The
> colorsys module can transform from/to RGB.

even if i convert from rgb to hsl i will have a tuple(h,s,l) for each
pixel and again i will have to convert it into a single value which i
can use in matrix multipln etc

is there a workaround sothat rgb color images can be worked on? any