matlab图像加权平均函数,matlab – 从电泳凝胶图像测量加权平均长度

这是我的尝试.它要求凝胶很好(即直通道和凝胶不应旋转),但应该相当普遍地工作.请注意,需要调整两个与图像大小相关的参数,以使其适用于不同大小的图像.

%# first size-dependent parameter: should be about 1/4th-1/5th

%# of the lane width in pixels.

minFilterWidth = 10;

%# second size-dependent parameter for filtering the

%# lane profiles

gaussWidth = 5;

%# read the image, normalize to 0...1

img = imread('http://img823.imageshack.us/img823/588/gele.png');

img = rgb2gray(img);

img = double(img)/255;

%# Otsu thresholding to (roughly) find lanes

thMsk = img < graythresh(img);

%# count the mask-pixels in each columns. Due to

%# lane separation, there will be fewer pixels

%# between lanes

cts = sum(thMsk,1);

%# widen the local minima, so that we get a nice

%# separation between lanes

ctsEroded = imerode(cts,ones(1,minFilterWidth));

%# use imregionalmin to identify the separation

%# between lanes. Invert to get a positive mask

laneMsk = ~repmat(imregionalmin(ctsEroded),size(img,1),1);

带有用于分析的车道的图像

%# for each lane, create an averaged profile

lblMsk = bwlabel(laneMsk);

nLanes = max(lblMsk(:));

profiles = zeros(size(img,1),nLanes);

midLane = zeros(1,nLanes);

for i = 1:nLanes

profiles(:,i) = mean(img.*(lblMsk==i),2);

midLane(:,i) = mean(find(lblMsk(1,:)==i));

end

%# Gauss-filter the profiles (each column is an

%# averaged intensity profile

G = exp(-(-gaussWidth*5:gaussWidth*5).^2/(2*gaussWidth^2));

G=G./sum(G);

profiles = imfilter(profiles,G','replicate'); %'

%# find the minima

[~,idx] = min(profiles,[],1);

%# plot

figure,imshow(img,[])

hold on, plot(midLane,idx,'.r')