#define get_channel(i, j, frame, channel) ((uchar *)(frame->imageData + j*frame->widthStep))[i*frame->nChannels + channel]

#define set_channel(i, j, frame, channel, value) ((uchar *)(frame->imageData + j*frame->widthStep))[i*frame->nChannels + channel] = value

int _tmain(int argc, _TCHAR* argv[])
{
    //Load the image
    IplImage *img = cvLoadImage(“beckham.bmp”);
    IplImage *imgGT = cvLoadImage(“beckhamgt.bmp”);
   
    //Show the image in a window
    cvNamedWindow(“Image:”, 1);
    cvShowImage(“Image:”, img);

   
   
    int TP = 0,TN = 0,FP = 0,FN = 0;

    pHSVImg = cvCreateImage(cvSize(img->width, img->height), img->depth, 3);
    pHSVgt = cvCreateImage(cvSize(imgGT->width, imgGT->height), imgGT->depth, 3);
            
    //Convert to HSV
    cvCvtColor(img,pHSVImg,CV_BGR2HSV);
   
    
    //Write the image in a different file
    if( !cvSaveImage(“picture_HSV.bmp”, pHSVImg))
    {
        fprintf(stderr, “fail to write image file\n”);
    }
    
    //Remember to release the memory after using it
    

   
    //convert to hsv, detect hue 
    IplImage *thresholdingFrame = cvCreateImage(cvSize(img->width, img->height), IPL_DEPTH_8U, 3);
   
   
    for(int i = 0; i < pHSVImg->width; i++)
    {
        for(int j = 0; j < pHSVImg->height; j++)
        {   
            int hue = get_channel(i, j, pHSVImg, 0) * 2; //openCV divides the hue by 2 to fit in the range [0,255]

            float saturation = get_channel(i, j, pHSVImg, 1) ;

            //float value = get_channel(i, j, hsvFrame, 2) / 255;
           
            /*if (hue < minHue)
            {
                minHue = hue;
                printf(“H min: %d\n”, minHue);
            }
           
            if (hue > maxHue)
            {
                maxHue = hue;
                printf(“H max: %d\n”, maxHue);
            }*/
           
            set_channel(i, j, pHSVImg, 1, hue / 2);
            set_channel(i, j, pHSVImg, 0, 0);
            set_channel(i, j, pHSVImg, 2, 0);
            set_channel(i, j, pHSVImg, 1, saturation );

            if (hue > 6 && hue < 38 )
            {
                set_channel(i, j, thresholdingFrame, 0, 255);
                set_channel(i, j, thresholdingFrame, 1, 255);
                set_channel(i, j, thresholdingFrame, 2, 255);
            }
            else
            {
                set_channel(i, j, thresholdingFrame, 0, 0);
                set_channel(i, j, thresholdingFrame, 1, 0);
                set_channel(i, j, thresholdingFrame, 2, 0);
            }
        }
    }
   
if( !cvSaveImage(“picture_thresholding.bmp”, thresholdingFrame))
    {
        fprintf(stderr, “fail to write image file\n”);
    }

//processing to skin detection
IplConvKernel* element = cvCreateStructuringElementEx( 5, 5, 2, 2, CV_SHAPE_ELLIPSE);
    cvDilate(thresholdingFrame, thresholdingFrame, element);
    cvErode(thresholdingFrame, thresholdingFrame, element);

cvSmooth(thresholdingFrame,thresholdingFrame,CV_MEDIAN, 3);
cvSmooth(thresholdingFrame,thresholdingFrame,CV_MEDIAN, 3);

//loop to read pixel per pixel and account if they are equal
for(int i = 0; i < thresholdingFrame->width; i++)
    {
        for(int j = 0; j < thresholdingFrame->height; j++)
        {   
           
                    int imgMy = get_channel(i,j,thresholdingFrame,0);
                    int imgGround = get_channel(i,j,imgGT,0);

                    if(imgMy == 255 && imgGround == 255)
                    { TP++;}

                    if(imgMy == 255 && imgGround == 0)
                    { FP++;}

                    if(imgMy == 0 && imgGround == 255)
                    { FN++;}

                    if(imgMy == 0 && imgGround == 0)
                    { TN++;}

           
        }
}
   

    printf(“TP is = %d\t”,TP);
    printf(“TN is = %d\t”,TN);
    printf(“FP is = %d\t”,FP);
    printf(“FN is = %d”,FN);

    cvNamedWindow(“threshold:”, 1);
    //cvNamedWindow(“ImageHSVgt:”, 1);
    cvShowImage(“threshold:”, thresholdingFrame);
    cvShowImage(“ImageHSV:”, pHSVImg);
    //cvShowImage(“ImageHSVgt:”, pHSVgt);
    cvWaitKey();
    cvDestroyWindow(“Image:”);
    cvDestroyWindow(“ImageHSV:”);
    cvReleaseImage(&img);
    cvReleaseImage(&pHSVImg);
    cvReleaseImage(&imgGT);
    cvReleaseImage(&thresholdingFrame);
    //cvReleaseImage(&pHSVgt);