clear;close all;
fcol=['r','g','b']; %初始化色彩
for i = 1:2
    T=imread(strcat('IMG',num2str(i),'.jpg'));%讀入圖像，不顯示
    %巴特沃茲濾波，不顯示濾波后圖像
    for k=1:3
        A=T(:,:,k);
        for j=1:300
            B=A(j,:);
            C=double(B); %將整形轉換為雙精度型
            [M,N]=butter(8,0.4); %使用8階，截至頻率是0.4的巴特沃茲濾波器
            D=filter(M,N,C);
            A(k,:)=D;
            T(:,:,k)=A;
       end
    end
   %在濾波后圖像上截取目標圖
   %imshow(T);
%    [x,y]=ginput(4);
%    a=min(x);
%    b=min(y);
%    width=max(x)-min(x);
%    height=max(y)-min(y);
%    I=imcrop(T,[a,b,width,height]);% %在目標圖上求灰度（強度）值
   rect=imcrop(T,[880,680,300,110]);
   imshow(rect);
   %反應曲線
   x=1:300; y=100*ones(1,length(x)); %y=ones(a,b);表示形成a行b列全為1的數列。
                                 %詞句表示y=100的一條長為550的線段。
   Cur0=improfile(rect,x,y); %improfile用于沿著圖像中一條直線段路徑或直線路徑計算并繪制其強度（灰度）值。
                    %彩色圖像是將紅綠藍分離開了，個表達獨立的線條圖形。
   grid on; %加坐標格
   %三色線條顯示同一個圖上
   figure;
   for i=1:3
    Cur1(i,:)=Cur0(:,:,i)'; %將列向量轉換為行向量。
    plot(Cur1(i,:),fcol(i));
    hold on; grid on; 
    axis([0 300 110 190]);
    xlabel('水平位置');
    ylabel('光強分布');
    title('目標圖像的光強');
   end 
   figure,
   Cur2=Cur1(2,:);
   plot(Cur2,'g');grid on;

end