matlab 图像提取中心线

clc clear all close all T=imread('simulate.jpg'); % 读入图像 T=rgb2gray(T); % 灰度转换 T=roicolor(T,0,30); BC1=imcomplement(T); BC2 = bwmorph(BC1,'remove'); %提取边缘 BC2=imcomplement(BC2); imshow(BC2) [r c]=find(BC2==0); A=[r,c]; e=zeros(2560,1); f=zeros(2560,1); for i=1:1:2560 d=find(c==i); f(i)=length(d)+0.5; e(i)=sum(r(d)); end h=e./f; s=round(h); B=ones(1440,2560); X=1:1:2560; for X=1:1:2560 B(s(X),X)=0; end figure,imshow(B)