wiggleplot.m

这是matlab在地球物理数据处理方面的源码

% **********************************************
% *** Program to plot a seismic section      ***
% *** Written by Y. Bouzidi. Geophysics UofA ***
% *** January 1998                           ***
% **********************************************
%
      function [ht]=wigleplot(traces,offset,time,...
                  ntrpinch,nsecpinch,gain)
%
% *** get number of samples and traces ***
%
      [nt,nx]=size(traces);
      ampmax=max(max(abs(traces)));     
      dx=offset(2)-offset(1);
      dt=time(2)-time(1);
      for i=1:nx
        traces(:,i)=traces(:,i)./ampmax*dx*gain;
      end
      xticks=0:10:(nx-1)*dx+offset(1);
      yticks=0:nt*dt/10:nt*dt;
      hundreds=fix(nx/ntrpinch);
      xaxislength=hundreds+rem(nx-1,ntrpinch)/ntrpinch+...
                  offset(1)/ntrpinch/dx+2*dx*gain/ntrpinch;
      yaxislength=(nt-1)*dt/nsecpinch+.1;
      totwidth=xaxislength+1.25;
      totheight=yaxislength+1.25;
      xscale=1./totwidth;
      yscale=1./totheight;
      xaxispos=0.75*xscale;
      yaxispos=0.75*yscale;
      left=10;
      bottom=10;
      width=72*totwidth;
      height=72*totheight;
      if(height>600)
        bottom=-(height-600);
      end;
      rect1=[left,bottom,width,height];
      h=figure('position',rect1);
      rect2=[xaxispos yaxispos xaxislength*xscale yaxislength*yscale];
      axes1=axes('position',rect2,'box','on','drawmode','fast');
      hold on
      set(gca,'TickDir','out');
      set(gca,'YDir','reverse');
      set(gca,'xtick',xticks);
      set(gca,'ytick',yticks);
      set(gca,'xlim',[-ntrpinch*.1*dx*gain (nx-1)*dx+offset(1)+...
         ntrpinch*.1*dx*gain]);
      set(gca,'ylim',[-nsecpinch*.05 (nt+1)*dt+nsecpinch*.05]);
      xlabel('Distance (m)');
      ylabel('Two way time (sec)'); 
      for i=1:nx
        ht(i)=plot(traces(:,i)+offset(i),time,'k','erasemode','none');
      end
%
% *** end of program *****