-
Notifications
You must be signed in to change notification settings - Fork 1
Blur Implementation
Jonathan edited this page Nov 28, 2017
·
4 revisions
As for secuential implementation we use opencv to read Image data
we first read the Image:
char* imageName = argv[1];
Mat image;
image = imread( imageName, IMREAD_COLOR );
if( !image.data )
{
std::cout<<" No image data \n ";
return -1;
}Create the Blurred Image as a copy, and assign the radius of kernel
Mat blur = image.clone();
int radius;
if(argv[2]==NULL || argc<3)
radius = 3;
else
radius = atoi(argv[2]);Create and offset to keep things safe and get the time to do analysis later:
//always radius less than protect
int protect = radius + 1;
//take the time
struct timespec start, finish;
double elapsed;
clock_gettime(CLOCK_MONOTONIC, &start);We then, Iterate through the structure and replace the pixels In the former Blurred Image created
for (int i = protect; i < image.rows - protect; ++i)
{
for (int j = protect; j < image.cols - protect; ++j)
{
//(i+1,j) right
//(i-1,j) left
//(i,j+1) down
//(i,j-1) top
//RED
blur.at<Vec3b>(i,j)[0] = (image.at<Vec3b>(i+radius,j)[0] + image.at<Vec3b>(i-radius,j)[0] + image.at<Vec3b>(i,j+radius)[0] + image.at<Vec3b>(i,j-radius)[0])/4;
//GREEN
blur.at<Vec3b>(i,j)[1] = (image.at<Vec3b>(i+radius,j)[1] + image.at<Vec3b>(i-radius,j)[1] + image.at<Vec3b>(i,j+radius)[1] + image.at<Vec3b>(i,j-radius)[1])/4;
//BLUE
blur.at<Vec3b>(i,j)[2] = (image.at<Vec3b>(i+radius,j)[2] + image.at<Vec3b>(i-radius,j)[2] + image.at<Vec3b>(i,j+radius)[2] + image.at<Vec3b>(i,j-radius)[2])/4;
}
}We finally take the time elapsed and write the output.
clock_gettime(CLOCK_MONOTONIC, &finish);
elapsed = (finish.tv_sec - start.tv_sec);
elapsed += (finish.tv_nsec - start.tv_nsec) / 1000000000.0;
std::cout<<"Time: "<< elapsed << std::endl;
//write in file
std::ofstream output;
output.open("mainResult.txt", std::ofstream::app | std::ofstream::out );
output<< imageName <<", "<< radius << ", "<< elapsed << std::endl;
output.close();Blur
POSIX Threads (pthreads)
Threads (openMP)
GPU (cuda)
How to run