~gerald-mwangi/+junk/ThesisMatlab

    clear all

    close all

    clc;

 

imstr={'Grove','Army','Evergreen','Dumptruck','Basketball','MiniCooper','Beanbags','Urban','Schefflera','Mequon'};

 

for s=1:length(imstr)

 

 

    load(strcat('./ResultsGeneralNewton/',char(imstr(s)),'.mat'));% load energies, prior energies and curvatures

 

    meanengen=[];

        meanenbas=[];

        vargen=[];

        varbas=[];

        for i=1:algo_iterations

            mgen=mean(mean(energies_gen(:,i)));

            mbas=mean(mean(energies_bas(:,i)));

            vgen=mean(mean((energies_gen(:,i)-mgen).^2));

            vbas=mean(mean((energies_bas(:,i)-mbas).^2));

            meanengen=[meanengen,mgen];

            meanenbas=[meanenbas,mbas];

            vargen=[vargen,vgen];

            varbas=[varbas,vbas];

        end

        

            meancurvgen=[];

        meancurvbas=[];

        varcurvgen=[];

        varcurvbas=[];

        for i=1:algo_iterations

            mgen=mean(mean(curvatures_gen(:,i)));

            mbas=mean(mean(curvatures_bas(:,i)));

            vgen=mean(mean((curvatures_gen(:,i)-mgen).^2));

            vbas=mean(mean((curvatures_bas(:,i)-mbas).^2));

            meancurvgen=[meancurvgen,mgen];

            meancurvbas=[meancurvbas,mbas];

            varcurvgen=[varcurvgen,vgen];

            varcurvbas=[varcurvbas,vbas];

        end

 

                figure(3);plot(meanengen,'-','LineWidth',2);

                hold on

            plot(meanenbas,'--','LineWidth',2);

 

                hold off

                title('mean $$E(\phi,\nabla\phi)$$','Interpreter','latex');

        set(gca,'FontSize',30)

        set( get( gca, 'Title' ), 'fontSize', 30 )

        legend('- = MOF','-- = BAS')

        print(strcat('./ResultsGeneralNewton/',char(imstr(s)),'-MeanEnergy'),'-dpng')

 

                figure(4);plot(sqrt(vargen),'-','LineWidth',2);

                hold on

            plot(sqrt(varbas),'--','LineWidth',2);

 

                hold off

                title('std $$E(\phi,\nabla\phi)$$','Interpreter','latex');

        set(gca,'FontSize',30)

        set( get( gca, 'Title' ), 'fontSize', 30 )

        legend('- = MOF','-- = BAS')

        print(strcat('./ResultsGeneralNewton/',char(imstr(s)),'-stdDevEnergy'),'-dpng')

        

        

                        figure(3);plot(meancurvgen,'-','LineWidth',2);

                hold on

            plot(meancurvbas,'--','LineWidth',2);

 

                hold off

                title('mean $$\left\vert \mathbf{K}\right\vert^2$$','Interpreter','latex');

        set(gca,'FontSize',30)

        set( get( gca, 'Title' ), 'fontSize', 30 )

        legend('- = MOF','-- = BAS')

        print(strcat('./ResultsGeneralNewton/',char(imstr(s)),'-MeanCurvature'),'-dpng')

 

                figure(4);plot(sqrt(varcurvgen),'-','LineWidth',2);

                hold on

            plot(sqrt(varcurvbas),'--','LineWidth',2);

 

                hold off

                title('std $$\left\vert \mathbf{K}\right\vert^2$$','Interpreter','latex');

        set(gca,'FontSize',30)

        set( get( gca, 'Title' ), 'fontSize', 30 )

        legend('- = MOF','-- = BAS')

        print(strcat('./ResultsGeneralNewton/',char(imstr(s)),'-stdDevCurvature'),'-dpng')

        

        imwrite(uint8(imgen),strcat('./ResultsGeneralNewton/',char(imstr(s)),'-imgen.png'));

        imwrite(uint8(imbasic),strcat('./ResultsGeneralNewton/',char(imstr(s)),'-imbas.png'));

        imwrite(uint8(im),strcat('./ResultsGeneralNewton/',char(imstr(s)),'-im.png'));

        imwrite(uint8(im+noise),strcat('./ResultsGeneralNewton/',char(imstr(s)),'-imnoise.png'));

end