Hola muy buenas, tengo un problema con el siguiente ejercicio:
Una laboratorio de investigación celular necesita desarrollar un software para la detección de
células madre no embrionarias procedentes de cultivos celulares en placas de Petri. Para ello, el
laboratorio toma imágenes digitales de las placas con objeto de realizar una detección de las células
obtenidas y determinar su número y sus características. Se dispone de imágenes consideradas como
una matriz bidimensional de píxeles cuyo valor numérico indica la intensidad de la correspondiente
región de la foto (0 indica negro, y 255 indica blanco). Una imagen puede tener hasta 512 filas y
hasta 512 columnas.
Para detectar una célula madre se considerará que:
• Una célula madre ocupa región de 2x2 pixels cuyo promedio de intensidad es igual o mayor
de 240.
• La detección de las células madre en la imagen se realiza comenzando por la esquina
superior izquierda y siguiendo un orden de izquierda a derecha y de arriba hacia abajo.
• Un píxel no puede pertenecer a dos células madre distintas. Es decir, una vez que se ha
marcado un píxel como perteneciente a una célula, ya no se tomará en cuenta para la
detección de otras células.
Se pide implementar un programa en C++ para la detección de células madre en una imagen. Debe
diseñarse la función DetectarCelulas que recibe como argumento una imagen y que se
encarga de mostrar por la pantalla:
• La intensidad promedio y posición de la esquina superior izquierda de cada célula detectada.
• El número total de células detectadas.
• La intensidad promedio de todas las células detectadas.
Debe aparecer por pantalla lo siguiente como resultado de la ejecución del programa:
Célula madre en pos (1,2) con promedio de intensidad 242.50
Célula madre en pos (1,4) con promedio de intensidad 246.25
Célula madre en pos (3,3) con promedio de intensidad 240.00
Célula madre en pos (4,6) con promedio de intensidad 250.00
Célula madre en pos (5,1) con promedio de intensidad 240.00
5 células madre con un promedio de de intensidad 243.75
Yo he dividido el programa en 3 partes:
- El algoritmo principal:
Código C++:
Ver original#include <iostream>
using namespace std;
// Constante global, que indica el numero maximo de de filas y de columnas que puede tener el array.
const int TAM=512;
// Se definen dos arrays con un tamaño de fila y de columna de 512 como maximo. Uno tipo float para almacenar los numeros de la celula, y otro tipo
// bool para almecenar si la casilla está o no ya seleccionada por otra celula.
typedef float Timagen[TAM][TAM];
typedef bool Tmarcar[TAM][TAM];
struct TCelula {
Tmarcar marcel;
Timagen imacel;
};
bool marcarcelula (TCelula &p);
void DetectarCelulas (TCelula &p);
int main() {
int ancho, alto, intensidadmaxima;
Timagen imagen;
cout << "Dime el ancho de la imagen: ";
cin >> ancho;
if (ancho>TAM) {
cout << "El ancho no puede ser mayor a 512. Vuelva a introducir el ancho: ";
cin >> ancho;
}
cout << "Dime el alto de la imagen: ";
cin >> alto;
if (alto>TAM) {
cout << "El alto de la imagen no puede ser mayor de 512. Vuelva a introducir el alto: ";
cin >> alto;
}
intensidadmaxima=255;
for (int j=0; j<alto; j++) {
for (int i=0; i<ancho; i++) {
cout << "Introduce el valor del pixel en la posicion (" << j << "," << i << "): " << endl;
cin >> imagen[j][i];
if (imagen[j][i] > intensidadmaxima) {
cout << "Introduce el valor del pixel en la posicion (" << j << "," << i << ") dentro del rango: " << endl;
cin >> imagen[j][i];
}
}
}
---Aqui irian luego las llamadas a las funciones.
return 0;
}
- Un subalgoritmo que busca una celula madre 2x2 con una intensidad promedio de 240 y la marca a true o false. Y aqui es donde tengo el problema: que no se como asociar el valor bool al valor de intensidad para poder operar con el.
Código C++:
Ver originalbool marcarcelula (TCelula &p) {
bool marca;
marca=false;
float marcado;
for (int y=0; y<TAM; y++) {
for (int x=0; x<TAM; x++) {
while (p.marcel[y][x] = false) {
for (int j=0; j<TAM; j++) {
for (int i=0; i<TAM; i++) {
marcado=p.imacel[j][i]+p.imacel[j][i+1]+p.imacel[j+1][i]+p.imacel[j+1][i+1];
if (marcado>240) {
p.marcel[y][x] = true;
p.marcel[y][x+1] = true;
p.marcel[y+1][x] = true;
p.marcel[y+1][x+1] =true;
}
}
}
}
}
}
return marca;
}
- Subalgoritmo detectar celula:
Código C++:
Ver originalvoid DetectarCelulas (TCelula &p) {
float intenpro, intentotal, sumaintensidad;
float celdet;
celdet=0;
intentotal=0;
for (int y=0; y<TAM; y++) {
for (int x=0; x<TAM; x++) {
if ((p.marcel[y][x] = true) && (p.marcel[y][x+1] = true) && (p.marcel[y+1][x] = true) && (p.marcel[y+1][x+1] = true)) {
for (int j=0; j<TAM; j++) {
for (int i=0; i<TAM; i++) {
intenpro = ((p.imacel[j][i] + p.imacel[j][i+1] + p.imacel[j+1][i] + p.imacel[j+1][i+1])/4);
cout << "Celula madre en posicion (" << j << "," << i << ") con promedio de intensidad " << intenpro << endl;
celdet++;
sumaintensidad=sumaintensidad+intenpro;
}
}
}
}
}
intentotal=sumaintensidad/celdet;
cout << "Se han encontrado " << celdet << " celulas madre, con un promedio de intensidad de " << intentotal << endl;
}
Espero que lo hayais entendido, si no es asi decirlo y os lo intentare explicar mejor lo que quiero hacer para ver si me podeis dar una solucion.
Muchas gracias.