#include <cv.h>
#include <highgui.h>
#include <stdio.h>

uchar Blue[640][480];
uchar Green[640][480];
uchar Red[640][480];
uchar Gray[640][480];

int main()
{
        IplImage *Image1;
        Image1 = cvLoadImage("rain.jpg",1);
        cvNamedWindow("Color",1);
        cvShowImage("Color",Image1);
        cvWaitKey(0);

        /* Load Image RGB Values */
        for(int i=0;i<Image1->height;i++)
        {
                for(int j=0;j<Image1->widthStep;j=j+3)
                {
                        Blue[i][(int)(j/3)]=Image1->imageData[i*Image1->widthStep+j];
                        Green[i][(int)(j/3)]=Image1->imageData[i*Image1->widthStep+j+1];
                        Red[i][(int)(j/3)]=Image1->imageData[i*Image1->widthStep+j+2];
                }
        }

        /* Implement Algorithms */
        for(int i=0;i<Image1->height;i++)
        {
                for(int j=0;j<Image1->width;j++)
                {
                        Gray[i][j]=(uchar)(0.299*Red[i][j] + 0.587*Green[i][j] + 0.114*Blue[i][j]);
                        Red[i][j]=Gray[i][j];
                        Green[i][j]=Gray[i][j];
                        Blue[i][j]=Gray[i][j];
                }
        }

        /* Save Image RGB Values */
        for(int i=0;i<Image1->height;i++)
        {
                for(int j=0;j<Image1->widthStep;j=j+3)
                {
                        Image1->imageData[i*Image1->widthStep+j]=Blue[i][(int)(j/3)];
                        Image1->imageData[i*Image1->widthStep+j+1]=Green[i][(int)(j/3)];
                        Image1->imageData[i*Image1->widthStep+j+2]=Red[i][(int)(j/3)];
                }
        }

        cvSaveImage("rain_gray.jpg",Image1);
        cvNamedWindow("Gray Level",1);
        cvShowImage("Gray Level",Image1);
        cvWaitKey(0);

        cvReleaseImage(&Image1);
        cvDestroyWindow("Gray Level");
}