144
Sk_s = cov(train_soy_clean );
Sk_w = cov(train_woods );
else
% JUNCAO DOS PIXELS ROTULADOS COM OS SEMI-ROTULADOS
RotSR_cn = [train_corn_notill; SemiRot{1}];
RotSR_cm = [train_corn_min ; SemiRot{2}];
RotSR_c = [train_corn ; SemiRot{3}];
RotSR_sn = [train_soy_notill ; SemiRot{4}];
RotSR_sm = [train_soy_min ; SemiRot{5}];
RotSR_s = [train_soy_clean ; SemiRot{6}];
RotSR_w = [train_woods ; SemiRot{7}];
% JUNCAO DOS PESOS DOS PIXELS ROTULADOS COM OS PESOS DOS SEMI-ROTULADOS
PesoRotSR_cn = [ones(size(train_corn_notill,1),1); Peso{1}];
PesoRotSR_cm = [ones(size(train_corn_min,1),1) ; Peso{2}];
PesoRotSR_sn = [ones(size(train_soy_notill,1),1) ; Peso{4}];
PesoRotSR_sm = [ones(size(train_soy_min,1),1) ; Peso{5}];
PesoRotSR_s = [ones(size(train_soy_clean,1),1) ; Peso{6}];
PesoRotSR_w = [ones(size(train_woods,1),1) ; Peso{7}];
% CALCULO DO VETOR DE MEDIAS COM VALORES PESOS ESTIMADOS
Mk_cn = mediaEst(train_corn_notill, PesoRotSR_cn);
Mk_cm = mediaEst(train_corn_min, PesoRotSR_cm);
Mk_c = mediaEst(train_corn, PesoRotSR_c );
Mk_sn = mediaEst(train_soy_notill, PesoRotSR_sn);
Mk_sm = mediaEst(train_soy_min, PesoRotSR_sm);
Mk_s = mediaEst(train_soy_clean, PesoRotSR_s );
Mk_w = mediaEst(train_woods, PesoRotSR_w );
% CALCULO DA MATRIZ COVARIANCIA COM VALORES ESTIMADOS
Sk_cn = covEst(RotSR_cn, Mk_cn, PesoRotSR_cn);
Sk_cm = covEst(RotSR_cm, Mk_cm, PesoRotSR_cm);
Sk_c = covEst(RotSR_c, Mk_c, PesoRotSR_c );
Sk_sn = covEst(RotSR_sn, Mk_sn, PesoRotSR_sn);
Sk_sm = covEst(RotSR_sm, Mk_sm, PesoRotSR_sm);
Sk_s = covEst(RotSR_s, Mk_s, PesoRotSR_s );
Sk_w = covEst(RotSR_w, Mk_w , PesoRotSR_w );
end
NumTROCA{ITERACAO} = 0; % Verifica o numero de trocas na iteracao
for Xpos=1 : size(Imagem,1)
X = Imagem(Xpos,:); % Pega o Pixel X na posicao Xpos na imagem geral
% FUNCOES DECISAO (Calcula Max. Ver. Classes)
MV_cn = fdecisao(X, Mk_cn, Sk_cn);
MV_cm = fdecisao(X, Mk_cm, Sk_cm);
MV_c = fdecisao(X, Mk_c, Sk_c );
MV_sn = fdecisao(X, Mk_sn, Sk_sn);
MV_sm = fdecisao(X, Mk_sm, Sk_sm);
MV_s = fdecisao(X, Mk_s, Sk_s );
MV_w = fdecisao(X, Mk_w, Sk_w );
% DEFINE CLASSIFICACAO(Classe vencedora = pixel > Max Ver)
XmaxverClasses = [MV_cn MV_cm MV_c MV_sn MV_sm MV_s MV_w];
[Xmaxver, Xclasse] = max(XmaxverClasses);
% DEFINE O PESO DO PIXEL X
Xpeso = PesoEst(XmaxverClasses, Xclasse);
% SALVANDO OS DADOS
vetClassif(Xpos,1)= Xclasse; % Salva classe vencedora