convert color image to grayscale with pyopencl(example code)

grayscale.py

import pyopencl as cl
import sys
import Image
import numpy

def RoundUp(groupSize, globalSize):
	r = globalSize % groupSize;
	if r == 0:
		return globalSize;
	else:
		return globalSize + groupSize - r;

def main():

	imageObjects = [ 0, 0 ]

	if len(sys.argv) != 3:
		print "USAGE: " + sys.argv[0] + " <inputImageFile> <outputImageFile>"
		return 1

	# create context and command queue
	ctx = cl.create_some_context()
	queue = cl.CommandQueue(ctx)

	# load image
	im = Image.open(sys.argv[1])
	if im.mode != "RGBA":
		im = im.convert("RGBA")
	imgSize = im.size
	buffer = im.tostring()


	# Create ouput image object
	clImageFormat = cl.ImageFormat(cl.channel_order.RGBA, 
								cl.channel_type.UNSIGNED_INT8)
	imageObjects[0] = cl.Image(ctx,
								cl.mem_flags.READ_ONLY | cl.mem_flags.COPY_HOST_PTR,
								clImageFormat,
								imgSize,
								None,
								buffer)
	imageObjects[1] = cl.Image(ctx,
							cl.mem_flags.WRITE_ONLY,
							clImageFormat,
							imgSize)

	# load the kernel source code
	kernelFile = open("grayscale.cl", "r")
	kernelSrc = kernelFile.read()

	# Create OpenCL program
	program = cl.Program(ctx, kernelSrc).build()

	# Call the kernel directly
	localWorkSize = ( 16, 16 )
	globalWorkSize = ( RoundUp(localWorkSize[0], imgSize[0]),
					RoundUp(localWorkSize[1], imgSize[1]) )

	program.rgbaToGrayscale(queue,
							globalWorkSize,
							localWorkSize,
							imageObjects[0],
							imageObjects[1])

	# Read the output buffer back to the Host
	buffer = numpy.zeros(imgSize[0] * imgSize[1] * 4, numpy.uint8)
	origin = ( 0, 0, 0 )
	region = ( imgSize[0], imgSize[1], 1 )

	cl.enqueue_read_image(queue, imageObjects[1],
						origin, region, buffer).wait()

	print "Executed program succesfully."

	# Save the image to disk
	gsim = Image.fromstring("RGBA", imgSize, buffer.tostring())
	gsim.save(sys.argv[2])

main()

grayscale.cl

const sampler_t sampler = CLK_NORMALIZED_COORDS_FALSE | 
						  CLK_ADDRESS_CLAMP_TO_EDGE | 
						  CLK_FILTER_NEAREST;

__kernel void rgbaToGrayscale(__read_only image2d_t srcImg,
                              __write_only image2d_t dstImg)
{
	// Converts RGBA image to gray scale intensity using the following formula: 
	// I = 0.2126 * R + 0.7152 * G + 0.0722 * B

	int2 coord = (int2) (get_global_id(0), get_global_id(1));
	int width = get_image_width(srcImg);
	int height = get_image_height(srcImg);

	if (coord.x < width && coord.y < height)
	{
		uint4 color = read_imageui(srcImg, sampler, coord);
		float luminance = 0.2126f * color.x + 0.7152f * color.y + 0.0722f * color.z;
		color.x = color.y = color.z = (uint)luminance;

		// Write the output value to image
		write_imageui(dstImg, coord, color);
	}
}

[img]http://dl.iteye.com/upload/attachment/0081/3756/209a3012-150a-35c8-bf6b-2b7b93f7e612.jpg[/img]

[img]http://dl.iteye.com/upload/attachment/0081/3754/a0c07dd1-7338-3cfe-b2bb-9bc7d674e9c1.jpg[/img]