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

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,'optionname',optionvalue)
%      creation of an argus file newprofile
%
%   Available options:
%      - include: include list of existing ARGUS files
%      - color: line color (default='r')
%      - selectioncolor: line color of selected profiles (default='b')
%      - linestyle (default='-')
%      - linewidth (default=0.2)
%      - marker (default='+')
%      - markersize (default=7)
%      - markeredgecolor (default='r')
%
%   Usage:
%      expmaster(newfile,varargin)
%
%   Example:
%      expmaster('domain.exp','include',{'domain1.exp' 'domain2.exp'},'color','g','marker','+')
%
%   See also EXPDOC

%recover options
options=pairoptions(varargin{:});

%Some checks
if ~nargin | nargout
        error('expmaster usage: expmaster(newfile,varargin)')
elseif exist(newfile,'file'),
        %recursive call to expmaster if file already exists
        if ~exist(options,'include'),
                expmaster(newfile,'include',newfile,varargin{:});
                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

%Add default options
options=addfielddefault(options,'color','r');
options=addfielddefault(options,'selectioncolor','b');
options=addfielddefault(options,'LineStyle','-');
options=addfielddefault(options,'LineWidth',0.2);
options=addfielddefault(options,'Marker','+');
options=addfielddefault(options,'MarkerSize',7);
options=addfielddefault(options,'MarkerEdgeColor','r');

%put all the argus profiles given in input in one structure A
A=struct([]);
numprofiles=0;
numpoints=0;
closed=[];

%initialize the variables with files provided by 'include' option
if exist(options,'include'),
        files=getfieldvalue(options,'include');
        if ischar(files), files={files}; end
        for i=1:length(files),
                filename=files{i};
                if ~exist(filename,'file'),
                        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) & length(B(i).x)>1 )
                                        closed(numprofiles)=1;
                                else
                                        closed(numprofiles)=0;
                                end
                        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

%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)=[];

        %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
                        if length(A(i).x)==1,
                                plot(A(i).x,A(i).y,'color',getfieldvalue(options,'color'),'LineStyle',getfieldvalue(options,'LineStyle'),'LineWidth',getfieldvalue(options,'LineWidth'),...
                                        'MarkerEdgeColor',getfieldvalue(options,'MarkerEdgeColor'),'MarkerSize',getfieldvalue(options,'MarkerSize'),'Marker','o');
                        else
                                plot(A(i).x,A(i).y,'color',getfieldvalue(options,'color'),'LineStyle',getfieldvalue(options,'LineStyle'),'LineWidth',getfieldvalue(options,'LineWidth'));
                        end
                        prevplot2=prevplot2+1;
                end
        end

        %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,options);
                        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,options);
                        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,options);
                        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,options);
                        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,options);
                        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,options);
                        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,options);
                        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,options);
                        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,options);
                        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,options);
                        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,options);
                        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,options);
                        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

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;