pi_savefile.m

Univ. of Calgary CREWS的免费地震研究软件

function SaveFile
% Saving files for Plot Image
global SCALE_OPT NUMBER_OF_COLORS GRAY_PCT CLIP COLOR_MAP NOBRIGHTEN PICKS PICKCOLOR XAXISTOP ZOOM_VALUE ZOOM_LOCKS
h=get(gcf,'userdata');
hmsg2=h(2);
hmsg=h(4);
hi=h(5);
hscale=h(6);
hclip=h(7);
hmaster=h(10);
hlimbox=h(14);
limboxflag=get(hlimbox,'value');
ampflag=get(hmaster,'value');
stuff=get(hmaster,'userdata');
mxs=stuff(1);
smean=stuff(2);
stddev=stuff(3);
%get the data
seis=get(hi,'cdata');
%determine old scaling
dat=get(hscale,'userdata');
oldscaleopt=dat(1);
mxsold=dat(2); mnsold=dat(3); smeanold=dat(4); stddevold=dat(5);
%new opt
newscaleopt=get(hscale,'value');
dat(1)=newscaleopt;
%get clip value
inewclip=get(hclip,'value');
ioldclip=get(hclip,'userdata');
set(hclip,'userdata',inewclip);
clips=get(hmsg,'userdata');
clipold=clips(ioldclip);
clipnew=clips(inewclip);

%get number of columns in colormap
clrmap=get(gcf,'colormap');
nkols=size(clrmap,1);

flag=computer;
flag='shit';

if(ampflag==1)
    % Figure is Independant
    col=[.8314 .8157 .7843];
    TTS='';
    stddev2=stddev;
    mxs2=mxs;
    smean2=smean;
elseif(ampflag==2)
    % Figure is Master
    col=[1 .5 0];
    TTS='';
    global SMEAN STDDEV MXS
    stddev2=stddev;
    mxs2=mxs;
    smean2=smean;
    SMEAN=smean2;
    STDDEV=stddev;
    MXS=mxs;
elseif(ampflag==3)
    % while figure is slaved, automatically Max Scaling
    col=[1 1 0];
    TTS='Scaling Automatically Forced Master Figure Scaling';
    % Figure is Slave
    limlns=findobj(gcf,'type','line','tag','LIMITLINE');
    limpts=findobj(gcf,'type','line','tag','LIMITPOINT');
    limdat=get(hlimbox,'userdata');
    if(~isempty(limlns))
        delete(limlns);
        delete(limpts);
        delete(limdat{3});
        set(hlimbox,'userdata',[]);
    end
    global SMEAN STDDEV MXS
    stddev2=STDDEV;
    mxs2=MXS;
    smean2=SMEAN;
end
set(hmaster,'backgroundcolor',col,'tooltipstring',TTS);
% 		%undo the old scaling
if( oldscaleopt == 1 ) %undo mean scaling
    mxsprime = min([smeanold+clipold*stddevold,mxsold]);
    mns=-mxsprime;
    smat = (seis-1)*(mxsprime-mns)/(nkols-1) + mns;
elseif( oldscaleopt == 2) %undo max scaling
    mns=-mxsold;
    smat = (seis-1)*(mxsold-mns)/(nkols-1) + mns;
end
seis=smat;
t=get(hi,'ydata');
x=get(hi,'xdata');
dat=get(gcbo,'userdata');
val=get(gcbo,'userdata');
switch val{1}
    case 1
        [file,path]=myuifile(gcf,'*.*','Saving File as SEGY','put');
        if(path==0)
            return
        elseif(isempty(findstr(file(end+4:end),'.sgy')))
            file=[file '.sgy'];
        end
        filename=[path file];
        dt=t(2)-t(1);
        altwritesegy(filename,seis,dt);
        AddMenu(filename);
    case 2
        [file,path]=myuifile(gcf,'*.mat','Saving File as .mat','put');
    if(path==0)
        return
    end
    if(isempty(findstr(file,'.mat')))
        file=[file '.mat'];
    end
    filename=[path file];
    SeisStruct.data=seis;
    SeisStruct.x=x;
    SeisStruct.t=t;
    SeisStruct.lastsave=['Plot Image.  ' date];
    cmt=['At this point, there is no way to tell whether ',...
            'the y (depth) is in time or distance.  The default ',...
            'is time.']; 
    SeisStruct.comment=cmt;
    save(filename,'SeisStruct');
    AddMenu(filename);
case 3
    save('plotimageproperties.mat','SCALE_OPT','NUMBER_OF_COLORS','GRAY_PCT','CLIP','COLOR_MAP','NOBRIGHTEN','PICKCOLOR','XAXISTOP');
    return
end
stringinfo=['Data has been save to "' filename '"'];
set(hmsg2,'string',stringinfo,'backgroundcolor',[1 1 1]);

function AddMenu(filename)
% Adding new Menu to quick menus
QuickMenu=findobj(gcf,'type','uimenu','tag','QUICK_OPENFILE');
udat=get(QuickMenu,'userdata');
CheckMenu=udat{3};
HoldNames=get(CheckMenu,'label');
set(CheckMenu(1),'label',filename,'visible','on');
jj=1;
ii=2;
while jj<=3
    if(strcmp(HoldNames{jj},filename))
        jj=jj+1;
    else
        if(isempty(deblank(HoldNames{jj})))
            vis='off';
        else
            vis='on';
        end
        set(CheckMenu(ii),'label',HoldNames{jj},'visible',vis);
        jj=jj+1;
        ii=ii+1;
    end
end
savenames=get(CheckMenu,'label');
save('pltimgfilemenu','savenames');