Reporte de actividades del proyecto de Investigación en ...
Reporte de actividades del proyecto de Investigación en ...
Reporte de actividades del proyecto de Investigación en ...
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
La realización <strong>de</strong>l programa #3 se implem<strong>en</strong>tó con máscaras <strong>de</strong> 3x3 (ver<br />
capítulo 4). Se com<strong>en</strong>zó con la segunda fila y segunda columna <strong>de</strong> la imag<strong>en</strong> para<br />
terminar con la p<strong>en</strong>última fila y p<strong>en</strong>última columna <strong>de</strong> la misma imag<strong>en</strong>; pero la<br />
imag<strong>en</strong> resultante para cada filtro ti<strong>en</strong>e dos columnas y dos filas m<strong>en</strong>os que la<br />
original.<br />
clear all; close all;<br />
im_color=imread('sacapuntas1.bmp');<br />
im_color = imadjust(im_color,stretchlim(im_color));<br />
rojo=double(im_color(:,:,1));<br />
ver<strong>de</strong>=double(im_color(:,:,2));<br />
azul=double(im_color(:,:,3));<br />
im_gris=round(rojo(:,:)*0.33)+round(ver<strong>de</strong>(:,:)*0.59)+round(azul(:,:)*0.11);<br />
im=double(im_gris);<br />
ruido_aleatorio=floor(rand(fil,col)*100);<br />
im_ruido=im+ruido_aleatorio;<br />
im=im_ruido;<br />
[fil col]=size(im);<br />
for x=2:1:fil-1<br />
for y=2:1:col-1<br />
a=im(x-1,y-1);b=im(x,y-1);c=im(x+1,y-1);d=im(x-1,y);e=im(x,y);f=im(x+1,y);<br />
g=im(x-1,y+1);h=im(x,y+1);i=im(x+1,y+1);<br />
im_fil_prom(x,y)=round((a+b+c+d+e+f+g+h+i)/9);<br />
vector(1:9)=sort([a b c d e f g h i]);<br />
im_fil_mediana(x,y)=round(sum(vector)/9);<br />
im_fil_minimo(x,y)=vector(1);<br />
im_fil_maximo(x,y)=vector(9);<br />
im_fil_punto_m(x,y)=round((vector(1)+vector(9))/2);<br />
<strong>en</strong>d<br />
<strong>en</strong>d<br />
figure<br />
hold on<br />
subplot(2,3,1); imshow(uint8(im_ruido));title('Imag<strong>en</strong> con ruido aleatorio');<br />
subplot(2,3,2); imshow(uint8(im_fil_prom));title('Filtro promedio');<br />
subplot(2,3,3); imshow(uint8(im_fil_mediana));title('Filtro mediana');<br />
subplot(2,3,4); imshow(uint8(im_fil_minimo));title('Filtro minimo');<br />
subplot(2,3,5); imshow(uint8(im_fil_maximo));title('Filtro maximo');<br />
subplot(2,3,6); imshow(uint8(im_fil_punto_m));title('Filtro punto medio');<br />
Programa #3: Filtros para imág<strong>en</strong>es <strong>en</strong> tonos <strong>de</strong> gris.<br />
89