Matrix<double> GetModules (const Matrix<int> & sobel_horizontal,
const Matrix<int> & sobel_vertical) {
    Matrix<double> modules(sobel_horizontal.n_rows, sobel_horizontal.n_cols);
        for (uint i = 0; i < modules.n_rows; i++) {
            for (uint j = 0; j < modules.n_cols; j++) {
            modules(i, j) = sqrt(sobel_horizontal(i, j) * sobel_horizontal(i, j) + sobel_vertical(i, j) * sobel_vertical(i, j));
        }
    }
    return modules;
}