void umbralizarImagen(Matriz &m, int coefInicial, int coefFinal, double umbral) { // El par (x0, y0) indica el primer elemento de la diagonal actual int x0 = 0; int y0 = 0; // El par (x, y) indica el coeficiente actual int x = x0; int y = y0; for(int i = 0; i < m.filas() * m.columnas(); i++) { // Umbralizo el coeficiente actual if(coefInicial <= i && i <= coefFinal) { if(abs(m.elem(y, x)) < umbral) m.elem(y, x) = 0; } // Avanzo al siguiente elemento de la diagonal x--; y++; // Verifico si llegué al final de la diagonal actual if(x < 0 || y >= m.filas()) { // Avanzo a la siguiente diagonal if(x0 < m.columnas() - 1) x0++; else y0++; // Avanzo al primer elemento de la nueva diagonal x = x0; y = y0; } } }
void agregarRuidoImpulsivo(Matriz &m, const double p) { int max, min; unsigned seed = chrono::system_clock::now().time_since_epoch().count(); std::default_random_engine generator(seed); std::uniform_real_distribution<double> distribution(0.0,1.0); if(m.columnas() == 1) { max = m.max(); min = m.min(); } else { max = 255; min = 0; } for(int i = 0;i < m.filas(); i++) { for(int j = 0;j < m.columnas(); j++) { double randNum = distribution(generator); if(randNum < p) { m.elem(i,j) = max; } else if(randNum >= 1 - p) { m.elem(i,j) = min; } } } }
void muestraMatriz(Matriz<T> &m){ //Aqui mostraremos la matriz for(int i=0; i<m.filas(); i++){ for(int j=0; j<m.columnas(); j++ ){ cout << m[i][j] << "\t"; } cout << endl; } }
void agregarRuidoAditivo(Matriz &m, const double mu, const double sigma) { unsigned seed = chrono::system_clock::now().time_since_epoch().count(); default_random_engine generator(seed); normal_distribution<double> distribution(mu,sigma); for(int i = 0; i < m.filas(); i++) { for(int j = 0; j < m.columnas(); j++) { //cout << "Que larga: " << distribution(generator) << endl; if(m.columnas() == 1) { m.elem(i,j) += distribution(generator); } else { m.elem(i,j) += floor(distribution(generator)); } } } }