Commit 5c8439191ab4347c270af53b77eb8c01d9520a6c

Authored by Arnaud Blanchard
1 parent 8547bdbe

Remove f_gradidient. use f_grad.sh instead

# Set the minimum version of cmake required to build this project
cmake_minimum_required(VERSION 2.6)
# Set the name of the project as the directory basename
project(f_gradient)
find_package(blc_channel)
find_package(blc_program)
add_definitions(${BL_DEFINITIONS})
include_directories(${BL_INCLUDE_DIRS})
add_executable(f_gradient f_gradient.cpp)
target_link_libraries(f_gradient ${BL_LIBRARIES})
# Copyright ETIS — ENSEA, Université de Cergy-Pontoise, CNRS (2011 - 2016)
# Author: Arnaud Blanchard (November 2016)
# This software is governed by the CeCILL v2.1 license under French law and abiding by the rules of distribution of free software.
# You can use, modify and/ or redistribute the software under the terms of the CeCILL v2.1 license as circulated by CEA, CNRS and INRIA at the following URL "http://www.cecill.info".
# As a counterpart to the access to the source code and rights to copy, modify and redistribute granted by the license,
# users are provided only with a limited warranty and the software's author, the holder of the economic rights, and the successive licensors have only limited liability.
# In this respect, the user's attention is drawn to the risks associated with loading, using, modifying and/or developing or reproducing the software by the user in light of its specific status of free software,
# that may mean that it is complicated to manipulate, and that also therefore means that it is reserved for developers and experienced professionals having in-depth computer knowledge.
# Users are therefore encouraged to load and test the software's suitability as regards their requirements in conditions enabling the security of their systems and/or data to be ensured  and, more generally, to use and operate it in the same conditions as regards security.
# The fact that you are presently reading this means that you have had knowledge of the CeCILL v2.1 license and that you accept its terms.
cd `dirname $0`
echo
echo "Checking `basename $PWD`"
echo "===================="
echo
echo "No check for now"
echo
//
// Created by Arnaud Blanchard on 22/12/15.
// Copyright ETIS 2015. All rights reserved.
//
#include "blc_channel.h"
#include "blc_program.h"
#include <unistd.h>
#include <math.h>
#define DEFAULT_OUTPUT_NAME ":gradient<pid>"
#define GRAD(x,y) grad.uchars[(x)+(y)*grad_width]
/** See how to use Filter of Scharr
*
* -3 0 3
* -10 0 10
* -3 0 3 / 32
*
*/
void grad_downscale(blc_array *grad, blc_array const *input, int scale){
int i, j, id, jd, value;
float sumx, sumy;
int in_width=input->dims[0].length;
int grad_width=grad->dims[0].length;
FOR(j, grad->dims[1].length){
FOR(i, grad_width){
sumx=sumy=0;
FOR(id, scale){
FOR(jd, scale){
value=input->uchars[i+id+(j+jd)*in_width];
sumx+=value/(float)(id-scale/2+0.5);
sumy+=value/(float)(jd-scale/2+0.5);
}
}
grad->uchars[i+j*grad_width]=CLIP_UCHAR((fabs(sumx)+fabs(sumy))/scale);
}
}
}
int main(int argc, char **argv){
blc_channel input, grad;
char const *output_name, *input_name, *downscale_str;
int in_width, in_height, grad_height, grad_width, downscale=1;
int i, j, id, jd;
blc_program_set_description("Compute grtadient of the image");
blc_program_add_option(&output_name, 'o', "output", "blc_channel-out", "channel name", DEFAULT_OUTPUT_NAME);
blc_program_add_option(&downscale_str, 'D', "downscale", "integer", "dowscale the resulting image", NULL);
blc_program_add_parameter(&input_name, "blc_channel-in", 1, "greyscale image to compute the gradient ", NULL);
blc_program_init(&argc, &argv, blc_quit);
if (downscale_str) SYSTEM_SUCCESS_CHECK(sscanf(downscale_str, "%d", &downscale), 1, "Parsing downscale '%s'", downscale_str);
if (strcmp(output_name, DEFAULT_OUTPUT_NAME)==0) SYSTEM_ERROR_CHECK(asprintf((char**)&output_name,":gradient%d", getpid()), -1, NULL);
input.open(input_name, BLC_CHANNEL_READ);
in_width=input.dims[0].length;
in_height=input.dims[1].length;
blc_loop_try_add_waiting_semaphore(input.sem_new_data);
blc_loop_try_add_posting_semaphore(input.sem_ack_data);
grad_width=in_width-downscale;
grad_height=in_height-downscale;
/*
grad_width=in_width/downscale-1;
grad_height=in_height/downscale-1;*/
grad.create_or_open(output_name, BLC_CHANNEL_WRITE, input.type, 'Y800', 2, grad_width, grad_height);
grad.publish();
blc_loop_try_add_waiting_semaphore(grad.sem_ack_data);
blc_loop_try_add_posting_semaphore(grad.sem_new_data);
if (downscale==1)
{
BLC_COMMAND_LOOP(0){
FOR(j, grad_height){
FOR(i, grad_width){
grad.uchars[i+j*grad_width]=(abs(input.uchars[i+j*in_width]-input.uchars[i+j*in_width+1])+abs(input.uchars[i+j*in_width]-input.uchars[i+(j+1)*in_width]))/2;
}
}
}
}
else {
BLC_COMMAND_LOOP(0){
grad_downscale(&grad, &input, downscale);
}
}
return EXIT_SUCCESS;
}