%本程序用于對訓練樣本提取獨立主元，作為樣本特征，并送入SVM分類器中訓練
%圖像的預處理中不取對數，也無須做幅度歸一，由ICA的應用條件決定的。
%預處理后的圖像以向量的形式按行排列

%%%%%%%%%%%    圖像預處理部分    %%%%%%%%%%%%%%%%%%%%%
function main
clear all
disp('程序正在運行')
clear all
load T72_mag.mat
[T72]=subfunction1(mobandata);
load T72_fw.mat
fw_data=FW;
[fw_T72]=subfunction2(fw_data,T72);
clear mobandata;
clear FW;

load BMP_mag.mat
[BMP]=subfunction1(mobandata);
load BMP_fw.mat
fw_data=FW;
[fw_BMP]=subfunction2(fw_data,BMP);
clear mobandata;
clear FW;

load BTR_mag.mat
[BTR]=subfunction1(mobandata);
load BTR_fw.mat
fw_data=FW;
[fw_BTR]=subfunction2(fw_data,BTR);
clear mobandata;
clear FW;

%%%%%%%%%%%%  提取主元   %%%%%%%%%%%%%%%
m1=1;
m2=1;
m3=1;
for n=1:6
    L1=fw_T72(n);
    L2=fw_BMP(n); 
    L3=fw_BTR(n); 
    MIXimage_s=[T72(m1:m1+L1-1,:);BMP(m2:m2+L2-1,:);BTR(m3:m3+L3-1,:)];
    MIXimage_mean(n,:)=mean(MIXimage_s);
    X=MIXimage_s-ones(size(MIXimage_s,1),1)*MIXimage_mean(n,:);%?123*2113
    R=X*X';%123*123
    [v lambda]=eig(R);
    m=size(v,1);
    for i=1:m-40
        L(i)=lambda(m+1-i,m+1-i);
        u(:,i)=1/sqrt(L(i))*X'*v(:,m+1-i);
    end
    Proaxis{n}=u;%2113*83
    feature=X*u; %123*83
    clear L u;  
    %%%%%% SVM訓練 %%%%%%%%%%%%%%
    Parameters=[2,1,0.75,1,10,35,0.001,0,0.3,0.3];   %線性分類器 
    Samples=[feature'];
    Label=[1*ones(1,L1) 2*ones(1,L2) 3*ones(1,L3)];%row
    [Alphay,SVs,Bias,Para,Ns]=osuSVMTrain(Samples,Label,Parameters,0);
    svm_classpara{n}={Alphay SVs Bias Para Ns};
    m1=m1+L1;
    m2=m2+L2;
    m3=m3+L3; 
end
save traindata.mat Proaxis MIXimage_mean svm_classpara;
disp('程序結束')
    
%%%%%%%%%%%%%%%%%%%%%%     SUBFUNCTION 1       %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%    This function is to realize fft of  chips %%%%%%%%%%%%%%%%%%%%
 function [chips]=subfunction1(mobandata)
     for k=1:length(mobandata)
        chp=log10(mobandata{k}(33:97,33:97)+1);%取出65*65 不加一得出的是負值 
         fft_chip=fftshift(abs(fft2(chp)));%傅列變換，并把頻域中的圖像移到軸坐標中心
         fft_chip=fft_chip/max(max(fft_chip)); %幅度歸一
         cut_chip1=fft_chip(1:65,33:65);%取出一半幅頻信息
         X_chip1=reshape(cut_chip1,1,65*33);%形成列矢量
         chips(k,:)=X_chip1(33:65*33);%去掉一半直流    ？
     end       
%%%%%%%%%%%%%%%%%%%%%%%      SUBFUNCTION 2     %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%  This function is to separate the chip into subgroups   %%%%%%%%%%%%%
function [sub_chips]=subfunction2(fw_data,chips);
Az1=0;
Az2=60;         %*************參數方位組有關,10表示方位范圍

for n=1:6      %每10度方位內形成一個模板，共有4個。%********參數方位組有關
    fw1=find(fw_data>=Az1);
    fw2=find(fw_data>Az2);
    sub_chips(n)=abs(length(fw1)-length(fw2));%尋找10度方位內圖象的個數。
    Az1=Az1+60;  %********參數方位組有關
    Az2=Az2+60;  %********參數方位組有關
end