guisegymat.m

这是用matlab对segy数据进行处理

zoom on

axes(handles.axBot)
bar(traces,data.HeaderInfo.cdp)
axBot=axis;axis([ax(1) ax(2) axBot(3) axBot(4)]); %% MAKE SURE X AXIS IS THE SAME
ylabel('CDP');
zoom on




    % Position axes
    axes(handles.axMain)
    set(handles.axMain,'XTickLabel','')
    axMainPOS=get(handles.axMain,'Position');

    axes(handles.axTop)
    set(handles.axTop,'XAxisLocation','top')
    set(handles.axTop,'YAxisLocation','Right')
    axTopPOS=get(handles.axTop,'Position');
    axTopPOS(2)=axMainPOS(2)+axMainPOS(4);
    axTopPOS(1)=axMainPOS(1);
    axTopPOS(3)=axMainPOS(3);
    set(handles.axTop,'Position',axTopPOS);

    axes(handles.axBot)
    set(handles.axBot,'YAxisLocation','Right')
    set(handles.axBot,'XAxisLocation','Bottom')
    axBotPOS=get(handles.axBot,'Position');
    axBotPOS(2)=axMainPOS(2)-axBotPOS(4);
    axBotPOS(1)=axMainPOS(1);
    axBotPOS(3)=axMainPOS(3);
    set(handles.axBot,'Position',axBotPOS);






actionSaveSettings(h,handles);
UpdatePlotSettings(h,handles);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function actionSaveSettings(h,handles)
data=guidata(h);
filename = ['.',mfilename,'.mat'];
Position=get(handles.figsegy,'Position');
SETTINGS=data.SETTINGS;
save(filename,'Position','SETTINGS');



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function actionUpdateSegyHeader(h,handles)
data=guidata(h);
if isfield(data.SegyHeader,'Line')==0, data.SegyHeader.Line=1001;end
if isfield(data.SegyHeader,'Job')==0, data.SegyHeader.Job=1002;end
if isfield(data.SegyHeader,'Reel')==0, data.SegyHeader.Reel=1003;end
set(handles.eLine,'String',num2str(data.SegyHeader.Line))
set(handles.eJob,'String',num2str(data.SegyHeader.Job))
set(handles.eReel,'String',num2str(data.SegyHeader.Reel))
set(handles.ens,'String',num2str(data.SegyHeader.ns))
set(handles.edt,'String',num2str(data.SegyHeader.dt))
guidata(h,data)

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function actionExit(h,handles);
uiresume

% --------------------------------------------------------------------
function actionSetColorMap(h, handles,cmap)
eval(['colormap(',cmap,')'])

cmap=colormap;
set(0,'DefaultFigureColormap',cmap)


%%% END ACTION LIST





% --------------------------------------------------------------------
function actionToolAmplitude(h, handles,action)
data=guidata(h);

if strcmp(action,'AGC');
  disp('agc')
  
  data.DataOld=data.Data;
  
  wlength=1;
  
  nt=length(data.SegyTraceHeaders);
  for it=1:nt
   disp([num2str(it),'/',num2str(nt)])
   ns=data.SegyTraceHeaders(it).ns;
   dt=data.SegyTraceHeaders(it).dt/1e+6;
   t=[1:1:ns].*dt;
   iwlength=wlength/dt;
   
   for is=1:iwlength:ns-iwlength;
     range=[is:1:(is+iwlength-1)];
     data.Data(range,it)=data.Data(range,it)./max(abs(data.Data(range,it)));
   end
   range=[is+1:1:ns];
   data.Data(range,it)=data.Data(range,it)./max(abs(data.Data(range,it)));
   
  end
  
  
end % END AGC

data.SETTINGS.cmin=min(data.Data(:));
data.SETTINGS.cmax=max(data.Data(:));

%data.SETTINGS.gain=1;


guidata(h,data);
actionPlotSegy(h,handles);
UpdatePlotSettings(h,handles);

% --------------------------------------------------------------------
function actionToolSort(h, handles,action)
data=guidata(h);
guidata(h,data);

%%% START OF ACTION LIST FOR TOOLS

% --------------------------------------------------------------------
function varargout = eSegyFile_Callback(h, eventdata, handles, varargin)




% --------------------------------------------------------------------
function varargout = pbGetSegyFile_Callback(h, eventdata, handles, varargin)
actionOpenSegy(h,handles)



% --------------------------------------------------------------------
function varargout = pbPlotSegyFile_Callback(h, eventdata, handles, varargin)
ActionPlotSegy(h,handles)


% --------------------------------------------------------------------
function varargout = pbPlotSegy_Callback(h, eventdata, handles, varargin)




% --------------------------------------------------------------------
function varargout = eLine_Callback(h, eventdata, handles, varargin)




% --------------------------------------------------------------------
function varargout = Untitled_1_Callback(h, eventdata, handles, varargin)




% --------------------------------------------------------------------
function varargout = Untitled_2_Callback(h, eventdata, handles, varargin)




% --------------------------------------------------------------------
function varargout = Untitled_3_Callback(h, eventdata, handles, varargin)




% --------------------------------------------------------------------
function varargout = pbRead_Callback(h, eventdata, handles, varargin)

actionReadSegy(h,handles);



% --------------------------------------------------------------------
function varargout = eCaxisMin_Callback(h, eventdata, handles, varargin)
data=guidata(h);
if ~isempty(str2num(get(h,'String')));
    data.SETTINGS.cmin = str2num(get(h,'String'));
end
guidata(h,data);
actionPlotSegy(h,handles)




% --------------------------------------------------------------------
function varargout = eCaxisMax_Callback(h, eventdata, handles, varargin)
data=guidata(h);
if ~isempty(str2num(get(h,'String')));
    data.SETTINGS.cmax = str2num(get(h,'String'));
end
guidata(h,data);
actionPlotSegy(h,handles)





% --------------------------------------------------------------------
function varargout = eCaxisGain_Callback(h, eventdata, handles, varargin)
data=guidata(h);
if ~isempty(str2num(get(h,'String')));
    gain=str2num(get(h,'String'));
    if gain<10^(get(handles.sCaxisGain,'min')), gain=19^(get(handles.sCaxisGain,'min')); end
    if gain>10^(get(handles.sCaxisGain,'max')), gain=19^(get(handles.sCaxisGain,'max')); end
    data.SETTINGS.gain = str2num(get(h,'String'));
end
guidata(h,data);
actionPlotSegy(h,handles)


% --------------------------------------------------------------------
function varargout = sCaxisGain_Callback(h, eventdata, handles, varargin)
data=guidata(h);
data.SETTINGS.gain = 10^(get(h,'value'));
guidata(h,data);
actionPlotSegy(h,handles)


% --------------------------------------------------------------------
function varargout = eJob_Callback(h, eventdata, handles, varargin)




% --------------------------------------------------------------------
function varargout = eReel_Callback(h, eventdata, handles, varargin)




% --------------------------------------------------------------------
function varargout = edt_Callback(h, eventdata, handles, varargin)




% --------------------------------------------------------------------
function varargout = ens_Callback(h, eventdata, handles, varargin)




% --------------------------------------------------------------------
function varargout = Untitled_6_Callback(h, eventdata, handles, varargin)




% --------------------------------------------------------------------
function varargout = Untitled_7_Callback(h, eventdata, handles, varargin)




% --------------------------------------------------------------------
function varargout = mColormaps_Callback(h, eventdata, handles, varargin)




% --------------------------------------------------------------------
function varargout = Untitled_9_Callback(h, eventdata, handles, varargin)





% --------------------------------------------------------------------
function varargout = mTools_Callback(h, eventdata, handles, varargin)




% --------------------------------------------------------------------
function varargout = Untitled_13_Callback(h, eventdata, handles, varargin)