close all
clear all

net=RBFnet_creat;
theta=0:179;

while(1)
    num=input('Please input the images name :','s');   % 提示輸入待識別的圖片的名稱（輸入1～8的數字即可）
    if abs(num)<abs('1') || abs(num)>abs('8')          % 輸入不存在的圖片名稱時退出程序 
        close all
        clear all
        return
    else
        m=strcat('test',num,'.jpg');        % 將樣本圖片的文件名合并成字符串
        x=imread(m);                        % 讀入樣本圖片
        subplot(2,2,1);
        imshow(x);                          % 顯示原始圖
        title('原始圖像');

        level=graythresh(x);                % 自動計算灰度閾值
        bw=im2bw(x,level);                  % 將讀入的圖片變成二值圖像
        subplot(2,2,2);
        imshow(bw);                         % 顯示二值圖
        title('二值圖像');

        [R,xp]=radon(bw,theta,180);
        subplot(2,2,3);
        imagesc(theta,xp,R)                 % 顯示radon變換的圖
        title('radon變換圖像');
        
        [v,d]=eig(R);                       % 計算ht的特征值
        for k=1:180
            p(k,1)=real(d(k,k));            % 將特征值作為網絡的輸入
        end
       
        r=sim(net,p) ;                      % 對網絡進行仿真
        r=round(r);
        subplot(2,2,4);                     % 畫出識別出的圖
        if r(1,1)==0 && r(2,1)==0 && r(3,1)==1
            tr=imread('T.jpg');
            imshow(tr);
            title('識別結果');
        elseif  r(1,1)==0 && r(2,1)==1 && r(3,1)==0
            sq=imread('S.jpg');
            imshow(sq); 
            title('識別結果');
        else  
            error('wrong');
        end 
    end
end