中心线提取.rar_条纹中心线_激光_激光提取_激光线提取_自适应激光条纹中心线提取算法

clear all;clc; close all; % I=imread('F:\1华科\D508\唐博\3.png'); % I=imread('F:\1华科\D508\唐博\x.bmp'); % I=imread('F:\1华科\D508\唐博\IMG_7290.JPG'); I=imread('F:\1华科\D508\唐博\双焊道.BMP'); G=rgb2gray(I); [h,w]=size(G); %获得图片尺寸 figure;imshow(G);title('灰度图'); gausFilter = fspecial('gaussian',[5,5],1.6); %高斯滤波 blurG=imfilter(G,gausFilter,'replicate'); figure;imshow(blurG); title('高斯滤波blurG'); [A,B]=zishiyinyuzhi1(blurG);title('提取中心线'); %----------用极大值法判断焊道最高点作为中心-------% [pks,locs] = findpeaks(A,'minpeakdistance',50); % [pks,locs] = findpeaks(A); figure; hold on; box on; plot(A); plot(locs,pks,'r^') hold off; m=length(A); acc=1; n=5;%5点滤波 for i=(n-1)/2+1:(m-(n-1)/2) p=i-(n-1)/2; q=i+(n-1)/2; A(acc)=sum(A(p:q))/n; acc=acc+1; end length_pks=length(pks); Max=zeros(length_pks,1);Max_x=zeros(length_pks,1);acc=1; mean_A=mean(A);%求平均值 for i=1:length(locs) if pks(i)>mean_A; Max(acc)= pks(i); Max_x(acc)=locs(i); acc=acc+1; end end Max_x(Max_x==0) = []; Max(Max==0) = []; figure; hold on; box on; plot(1:length(A),A); plot(Max_x,Max,'r^');axis equal hold off;