el procesamiento de procesamiento de imagen morfológica Matlab- (con código) imagen digital

1.Objectives:

1. Estudio de MATLAB binarios algoritmos de procesamiento de imágenes morfológicas utilizan utilización;
2. El dominio de la imagen función de operaciones de MATLAB morfológica básica;
3. Aprender la morfología de las aplicaciones básicas.

Contenido 2.Experiment:

1. Programación de imagen binaria sustancialmente morfológica (erosión, dilatación, operación de apertura y cierre); la selección de una imagen diferente elemento de estructura de filtro del objeto.
2. Para lograr suavizar el ruido de imagen de escala de grises y el borde imagen extraída por la operación de la morfología.

3.Experiment Principio:

Ver el libro.

4.Experiment Pasos de Resultados y conlusion:

1, imagen binaria transformación morfológica de

Necesidad de imagen binaria función de transformación morfológica de escritura:

functionnewbuf = BwFilter (oldbuf, seleccionar) en la expansión de la llamada de función MATLAB, la función de correlación y la corrosión imagen cribado de algoritmos, el correspondiente procesamiento de la imagen binaria, el resultado final se almacena en la newbuf matriz. Programación BwFilter función de la función (). Los elementos estructurales también se pueden crear con unos y ceros función de función.
Aquí Insertar imagen Descripción

2, las operaciones morfológicas de imagen de entrada, es decir, la erosión, dilatación, operación de apertura y cierre, el cambio de la estructura de la forma de elemento, el tamaño, rehacer los experimentos anteriores, los resultados de experimentos comparativos, la influencia de la operación del análisis de elementos de la estructura;

Aquí Insertar imagen Descripción
La figura expansor Aquí Insertar imagen Descripción2 Corrosión Aquí Insertar imagen Descripciónfigura operación de apertura 3 Aquí Insertar imagen Descripciónen la Fig. 4 la operación de cierre

Conv: En comparación con la imagen original, como resultado de la erosión como antes para hacer que cada dólar relación se hace más pequeño, tan útil para la eliminación de pico de ruido. El resultado será el valor de gradación de tal manera que cada píxel se hace mayor que el elemento anterior se expandió, es adecuado para la eliminación de bajo nivel de ruido. Sin embargo, dado que este experimento demasiada pérdida de la imagen original, por lo que la expansión de las partículas más grandes, sólo las partículas de menos después del grabado.

= Funcionamiento Erosión a la operación de apertura, una operación de dilatación entonces (se ve bien ligada a dos de destino separada). Es decir, segmentado imagen.

= Operación de dilatación Primera y cierre, y luego grabado operación (bloque miradas juntos dos de cierre de conexión fino). Que hacer la imagen más satisfactoria.

forma elemento estructurante cambia, el tamaño del efecto:

(1) Expansión
Aquí Insertar imagen Descripción
(2) Corrosión Aquí Insertar imagen Descripciónoperación de apertura Aquí Insertar imagen Descripcióny cierre de operación Aquí Insertar imagen DescripciónConv: el cambio de la forma y tamaño del elemento de estructura, respectivamente, utilizando un tipo de 'disco', 'cuadrado' 'unos' y la matriz de encargo, etc., determinado por el parámetro de tamaño. Para ampliar, por ejemplo, puede verse bajo el mismo tipo de tamaño del elemento estructural se incrementa, manchas blancas más graves parcialmente expandido, y diferentes tipos y grado de expansión de diferentes direcciones. Para el grabado de la imagen a ser corroído, la imagen para la operación de apertura es parte de conexión delgada interrumpido, mayor es el tamaño de los parámetros estructurales, la porción de imagen menos blanco; Para la operación de cierre, el aumento de la conexión de la imagen, la porción blanca se hace grande y borrosa, el mismo tamaño aumenta, un nivel más profundo, diferentes tipos de efectos.

3 a la figura FigP0936 (bubbles_on_black_background) .tif Ejemplo funciones completas son los siguientes:

(1) Extracción y partículas límite de la imagen fundida
(2) se solapan entre sí para extraer las partículas
(3) no se superponen las partículas extraídas

Nota: se pueden utilizar algoritmos de llenado (1) de la zona. Como se muestra imagen de origen figura en, la imagen se puede convertir a una imagen binaria y, a continuación, sometido a invertida, de manera que los resultados para las partículas de relleno de área (operación de cierre) estarán conectados a la frontera, y después se compara con la imagen original, de imagen de partículas se puede derivar en el límite de la imagen fuente conectada.
(2) por ataque químico la operación de expansión de imagen. imagen Plantilla para la primera operación de ataque químico, debido a la inevitable superposición de partículas más grandes que el área de las partículas individuales, que quedan así parte después de la operación de grabado para las porciones superpuestas de las partículas, y luego ser expandido (operación abierta) con el que el el funcionamiento de la fuente de imagen se comparan, se puede concluir superposición de imagen de partículas.
(3) después de la superposición de las partículas obtenidas restando su imagen de origen, la imagen resultante es la distribución de partículas independiente
(1). Aquí Insertar imagen Descripción(2) (3)
Aquí Insertar imagen Descripción
Aquí Insertar imagen DescripciónConv: (. 1) en la imagen original se llena, puede ser significativamente hasta el punto blanca original observado aumentos en la imagen original para lograr el efecto de filtrado de ruido. (2) para la operación de la erosión en la imagen original, y luego una operación de dilatación (se ve bien vinculados a dos de destino separada). La imagen principal se divide es decir, operación de apertura. (3) Después de restar a cabo fue pequeñas partículas, es decir partículas dispersas figura.

[Apéndice] código de implementación

Un programa de

function newbuf=BwFilter(oldbuf,select) %子程序
if select==1                  %如果选择膨胀    
   A=strel('disk',5);        %构造结构元素    
   newbuf=imdilate(oldbuf,A);%进行膨胀处理
else
    if select==2              %如果选择腐蚀        
       B=strel('disk',5);        
       newbuf=imerode(oldbuf,B);
    else
        if select==3          %如果选择开运算           
          C=strel('square',3);           
          newbuf=imopen(oldbuf,C);
        else D=strel('square',3); %如果选择闭运算
            newbuf=imclose(oldbuf,D);        
        end
    end
end
end

procedimiento 2

clear all;clc;
addpath('E:\数字图像处理\程序与图像\dipum_toolbox_2.0.2(only run)');%添加相应的.p文件 
f=imread('lax.tif'); %读图像
m=input(''); %输入select种类,1代表膨胀,2代表腐蚀,3代表开运算,4代表闭运算

f1=im2bw(f);       %将图像转化为二值图像
f2=BwFilter(f1,m); %调用BwFilter子函数进行处理
subplot(121)       %显示原图(二值图像)
imshow(f1);
title('原图');
subplot(122)       %显示转化后图像imshow(f2);title('变换后图像');

procedimiento 3

clear all;clc;close all;                   %清空图像,释放变量
f=imread('FigP0936(bubbles_on_black_background).tif');%读图像
f1=im2bw(f);                               %转化为二值图像
f2=~f1;                                    %对图像取反
figure;imshow(f2);title('取反后图像');     %显示取反图像
f3=imfill(f2,'holes');                     %调用函数进行区域填充
figure;imshow(f3);title('区域填充后图像'); %显示区域填充后图像
g=f;                                       %使g等于f图像
[m,n]=size(g);        %取g图像的大小
for i=1:m             %将区域填充的结果与原图进行逐点比较,若该点在填充结果时为白,这点在g图像中便记为黑
    for j=1:n
        if f3(i,j)==1
            g(i,j)=0;
        end    
    end
end
figure;imshow(g);title('比较后图像');  %得到与边界相连的颗粒图像    

clear all;clc;
addpath('E:\数字图像处理\程序与图像\dipum_toolbox_2.0.2(only run)');%添加相应的.p文件 
f=imread('FigP0936(bubbles_on_black_background).tif');%读图像
mask1=strel(ones(15,15));  %用ones函数构造结构元素
mask1f2=imerode(f,mask1);       %对原图进行腐蚀处理
f3=filter2(fspecial('average',7),im2double(f2));%对腐蚀处理后结果进行滤波,滤波器形式由fspecial定义
f3=medfilt2(f3);           %进行中值滤波
f3=im2uint8(f3);           %将滤波后结果转化为无符号八位整型形式
figure;
imshow(f3);title('腐蚀并滤波后图像');%显示腐蚀并滤波后图像
mask2=strel('ball',12,12); %构造结构元素mask2
f4=imdilate(f3,mask2);     %再进行膨胀处理
figure;imshow(f4);title('膨胀后图像');%显示膨胀后图像
f5=f;                      %使f5等于f
[m,n]=size(f);             %得到原图的大小
for i=1:m                  %逐点比较处理后图像与原图,并做限幅处理,处理后图像中小于等于40的,一律省去,大于该值的取原图中的白点
    for j=1:n
        if f4(i,j)<=40
           f5(i,j)=0;
        else            f5(i,j)=f(i,j);        
        end    
    end
end

figure;imshow(f5);title('限幅比较后图像(交叠颗粒图像)');%显示交叠图像颗粒图
f6=f-f5;                   %将原图和交叠图像相减,得到独立分布颗粒图像
figure;imshow(f6);title('相减后图像(独立分布颗粒图像)');%显示独立分布颗粒图

PD: Soy una nota de investigación de Morfología Matemática

https://blog.csdn.net/weixin_42784535/article/details/105125411

Publicado 17 artículos originales · ganado elogios 12 · vistas 1657

Supongo que te gusta

Origin blog.csdn.net/weixin_42784535/article/details/105125124
Recomendado
Clasificación