OSTU阈值 图像边缘增强算法 MATLAB

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %传统边缘检测算子增强_laplacian of Gaussian，二阶梯度为0时为边界 %作者：穆为磊 %日期：2013年8月18日 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% clc; %清空command中内容 clear all; %清除workspace中的变量 close all; %关闭所有现有figure % A=imread('test3.bmp');%'原图1恢复-6.67x5.tif' A=imread('LENA.bmp'); %rgb转灰度 if isrgb(A)==1 I_gray=rgb2gray(A); else I_gray=A; end A=I_gray; figure,imshow(A); title('原图'); % Aa=medianforfun(A,3,3); % A=mat2gray(Aa); % figure,imshow(A); % title('中值滤波'); Scale_Enhance=1.24; [m,n]=size(A); Th=Otsu(A); for i=1:m for j=1:n if A(i,j)>=Th M_Mark(i,j)=255; else M_Mark(i,j)=0; end end end [BV,DV]=BVandDV(A,Th) figure,imshow(M_Mark),title('Otsu法'); %% LOG=EDGE(A,'log');%roberts,prewitt,sobel,log, figure,imshow(LOG,[]); title('log边缘'); %LOG算子 % log1=[-2 -4 -4 -4 -2; % -4 0 8 0 -4; % -4 8 24 8 4; % -4 0 8 0 -4; % -2 -4 -4 -4 -2]; % LOG1=imfilter(A,log1); % LOG1=conv2(A,log1); % LOG2=mat2gray(LOG1); % figure,imshow(LOG2,[]); % title('LOG'); % LOG2=LOG2*255; IE=A; for i=1:m for j=1:n if LOG(i,j)==1 IE(i,j)=IE(i,j)*Scale_Enhance; if IE(i,j)>255 IE(i,j)=255; end else IE(i,j)=IE(i,j); end end end [BV,DV]=BVandDV(IE,M_Mark) %BV为背景强度均值、DV为目标轻度均值 figure,imshow(IE),title('LOG增强');