fanalys.m

来自「峰值检测」· M 代码 · 共 53 行

function [MAX MIN]=fanalys(Name)
% 
% 只要在命令窗口输入fanalys('t20.m')就可以得出结果。但是对于太暗
% 的图像（18，24）会出现错误。如果使用彩色而不使用灰度可能会好，但是你的CCD将来是黑白的，不
% 能处理彩色的，献艺暂时有问题。你自己想想吧。


clc
I=rgb2gray(imread(Name));
imshow(I);
% U=imread('t20.png');
% I=rgb2gray(U);
x=1:400; y=120*ones(1,length(x));
C=improfile(I,x,y); grid on;%用impofile.m取得的光强，三维的。
A=C(:,:,1)';%取其中一维。 
subplot(211);
imshow(I);
subplot(212);
plot(A);
grid on;


B=[];
k=1;
for i=1:length(A)
    if abs(A(i)-min(A))<10
        B(k)=i;       
        k=k+1;
    else
        continue
    end
end
 MAX= mean(A(B(1):B(length(B))));      
disp('MAX：')
disp(MAX);

C=A(1:B(1)-40);

B=[];
k=1;
for i=1:length(C)
    if abs(C(i)-min(C))<10
        B(k)=i;       
        k=k+1;
    else
        continue
    end
end
MIN= mean(A(B(1):B(length(B))));       
disp('MIN：')
disp(MIN);