source: issm/trunk/src/m/utils/Exp/expmaster.m@ 2652

function expmaster(newfile,varargin)
%EXPMASTER - allow to create, modify, add, cut, .. segments of domain outline together
%
%   this routine is used to create, modify, cut,... an Argus file (.exp)
%
%   expmaster(newprofile)
%      creation of an argus file newprofile
%
%   expmaster(newprofile,oldprofile)
%      the modification of the file oldprofile will be saved in newprofile
%
%   expmaster(newprofile,oldprofile1,oldprofile2,...)
%      the modification of the files oldprofile* will be saved in newprofile
%
%   Usage:
%      expmaster(newfile,varargin)
%
%   See also EXPDOC

%Some checks
if ~nargin | nargout
        error('expmaster usage: expmaster(newfile,varargin)')
elseif exist(newfile),
        %recursive call to expmaster if file already exists
        if nargin==1,
                expmaster(newfile,newfile);
                return;
        end

        %check modification
        choice=input(['A file ' newfile ' already exists, do you want to modify it? (y/n)'],'s');
        if ~strcmpi(choice,'y'),
                disp('no modification done ... exiting');
                return
        end
end

%put all the argus profiles given in input in one structure A
A=struct([]);
numprofiles=0;
numpoints=0;
closed=[];
for i=1:nargin-1
        filename=varargin{i};
        if ~exist(filename),
                error(['expmaster error message:, ' filename ' does not exist. Exiting...']);
        else
                %read file
                B=expread(filename,1);
                %go through all profiles of B
                for i=1:size(B,2)
                        %plug profile in A
                        if numprofiles
                                A(numprofiles+1)=B(i);
                        else
                                A=B(i);
                        end
                        %update numprofiles and numpoints
                        numpoints=numpoints+length(B(i).x);
                        numprofiles=numprofiles+1;
                        %figure out if the profile is closed or not
                        if (B(i).x(1)==B(i).x(end) & B(i).y(1)==B(i).y(end))
                                closed(numprofiles)=1;
                        else
                                closed(numprofiles)=0;
                        end
                end
        end
end

%Get root of newfile
[path root ext ver]=fileparts(newfile);

%get current figure
if ~isempty(get(0,'children')),%if there is already a figure (return the number of opened figures)
        set(gcf,'Renderer','zbuffer'); %fixes a bug on Mac OS X (not needed in future Matlab version)
        P=get(gcf,'position');
        F=getframe(gca);
        F=F.cdata;
        %get current axis
        xlim=get(gca,'Xlim');
        ylim=get(gca,'Ylim');
        %recreate x_m and y_m
        x_m=linspace(xlim(1),xlim(2),size(F,2));
        y_m=linspace(ylim(2),ylim(1),size(F,1)); %getframe reverse axis...
        %plot the data in another figure
        figure; set(gcf,'position',P);
        imagesc(x_m,y_m,F); set(gca,'Ydir','normal');
        prevplot=1;
        prevplot2=1;
else
        figure
        prevplot=0;
        prevplot2=0;
end

%plot existing profile if any

hold on
if numprofiles
        for i=1:numprofiles
                plot(A(i).x,A(i).y,'-r','MarkerSize',10);
        end
end

%Build backup structre for do and redo
backup=cell(1,3);
backup{1,1}=A;
backup{1,2}=numprofiles;
backup{1,3}=numpoints;
backup{1,4}=closed;

loop=1;
counter=1;
while loop

        %Go through A and rule out the empty profiles
        list=[];
        for i=1:size(A,2);
                if length(A(i).x)==0
                        list(end+1)=i;
                        numprofiles=numprofiles-1;
                end
        end
        A(list)=[];
        closed(list)=[];

        %display menu
        title('Main Menu','FontSize',14);
        button=menu('Menu','add a profile (open)',...%1
                'add a contour (closed)',...              %2
                'remove a profile',...                    %3
                'modify the position of a point',...      %4
                'add points inside a profile',...         %5
                'add points at the end of a profile',...  %6
                'remove points',...                       %7
                'remove several points',...               %8
                'cut a segment',...                       %9
                'cut a large area',...                    %10
                'merge profiles',...                      %11
                'close profile',...                       %12
                'undo',...                                %13
                'redo',...                                %14
                'quit');                                  %15


        %UNDO??
        if button==13;
                if counter==1
                        disp('Already at oldest change');
                else
                        counter=counter-1;
                        A=backup{counter,1};
                        numprofiles=backup{counter,2};
                        numpoints=backup{counter,3};
                        closed=backup{counter,4};
                end
        end

        %REDO??
        if button==14
                if counter==size(backup,1)
                        disp('Already at newest change');
                else
                        counter=counter+1;
                        A=backup{counter,1};
                        numprofiles=backup{counter,2};
                        numpoints=backup{counter,3};
                        closed=backup{counter,4};
                end
        end

        switch button

                case 1

                        [A,numprofiles,numpoints,closed]=addprofile(A,numprofiles,numpoints,closed,prevplot2,root);
                        counter=counter+1;
                        backup{counter,1}=A;
                        backup{counter,2}=numprofiles;
                        backup{counter,3}=numpoints;
                        backup{counter,4}=closed;

                case 2

                        [A,numprofiles,numpoints,closed]=addcontour(A,numprofiles,numpoints,closed,prevplot2,root);
                        counter=counter+1;
                        backup{counter,1}=A;
                        backup{counter,2}=numprofiles;
                        backup{counter,3}=numpoints;
                        backup{counter,4}=closed;

                case 3

                        [A,numprofiles,numpoints,closed]=removeprofile(A,numprofiles,numpoints,closed,prevplot2,root);
                        counter=counter+1;
                        backup{counter,1}=A;
                        backup{counter,2}=numprofiles;
                        backup{counter,3}=numpoints;
                        backup{counter,4}=closed;

                case 4

                        [A,numprofiles,numpoints,closed]=modifyposition(A,numprofiles,numpoints,closed,prevplot,root);
                        counter=counter+1;
                        backup{counter,1}=A;
                        backup{counter,2}=numprofiles;
                        backup{counter,3}=numpoints;
                        backup{counter,4}=closed;

                case 5

                        [A,numprofiles,numpoints,closed]=addinsideprofile(A,numprofiles,numpoints,closed,prevplot,root);
                        counter=counter+1;
                        backup{counter,1}=A;
                        backup{counter,2}=numprofiles;
                        backup{counter,3}=numpoints;
                        backup{counter,4}=closed;

                case 6

                        [A,numprofiles,numpoints,closed]=addendprofile(A,numprofiles,numpoints,closed,prevplot2,root);
                        counter=counter+1;
                        backup{counter,1}=A;
                        backup{counter,2}=numprofiles;
                        backup{counter,3}=numpoints;
                        backup{counter,4}=closed;

                case 7

                        [A,numprofiles,numpoints,closed]=removepoints(A,numprofiles,numpoints,closed,prevplot,root);
                        counter=counter+1;
                        backup{counter,1}=A;
                        backup{counter,2}=numprofiles;
                        backup{counter,3}=numpoints;
                        backup{counter,4}=closed;

                case 8

                        [A,numprofiles,numpoints,closed]=removeseveralpoints(A,numprofiles,numpoints,closed,prevplot,root);
                        counter=counter+1;
                        backup{counter,1}=A;
                        backup{counter,2}=numprofiles;
                        backup{counter,3}=numpoints;
                        backup{counter,4}=closed;

                case 9

                        [A,numprofiles,numpoints,closed]=cutprofile(A,numprofiles,numpoints,closed,prevplot,root);
                        counter=counter+1;
                        backup{counter,1}=A;
                        backup{counter,2}=numprofiles;
                        backup{counter,3}=numpoints;
                        backup{counter,4}=closed;

                case 10

                        [A,numprofiles,numpoints,closed]=cutarea(A,numprofiles,numpoints,closed,prevplot,root);
                        counter=counter+1;
                        backup{counter,1}=A;
                        backup{counter,2}=numprofiles;
                        backup{counter,3}=numpoints;
                        backup{counter,4}=closed;

                case 11

                        [A,numprofiles,numpoints,closed]=mergeprofiles(A,numprofiles,numpoints,closed,prevplot,root);
                        counter=counter+1;
                        backup{counter,1}=A;
                        backup{counter,2}=numprofiles;
                        backup{counter,3}=numpoints;
                        backup{counter,4}=closed;


                case 12

                        [A,numprofiles,numpoints,closed]=closeprofile(A,numprofiles,numpoints,closed,prevplot,root);
                        counter=counter+1;
                        backup{counter,1}=A;
                        backup{counter,2}=numprofiles;
                        backup{counter,3}=numpoints;
                        backup{counter,4}=closed;

                        %QUIT
                case 15

                        loop=0;

                otherwise

                        %do nothing

        end

        %Now erase all that have been done and plot the new structure A as it is
        undoplots(prevplot);
        if numprofiles
                prevplot2=1;
                for i=1:numprofiles
                        plot(A(i).x,A(i).y,'-r','MarkerSize',10);
                        prevplot2=prevplot2+1;
                end
        end
end

hold off

%write contour using expwrite
title('New file written, exiting...','FontSize',14);
if isempty(A)
        disp('Profile empty, no file written')
else
        expwrite(A,newfile);
end

%close window
close;