killint3.m,v

具有特色的地震数据处理源码

head	3.0;
access;
symbols;
locks; strict;
comment	@// @;


3.0
date	2000.06.13.19.22.48;	author gilles;	state Exp;
branches;
next	2.0;

2.0
date	99.05.21.18.47.11;	author mah;	state Exp;
branches;
next	1.1;

1.1
date	99.01.06.19.09.25;	author kay;	state Exp;
branches;
next	;


desc
@@


3.0
log
@Release 3
@
text
@function killint3(action,thm)

%interactive module for killing traces
%works with 'pickfb'
%saves killed trace data in matrix 'killflg'
%kills are only saved if 'save' is chosen on exiting program
%flagged traces are ignored by picking module and are diplayed as flat red lines
%thm is trace header matrix
%
%DSI customised processing software
%written by K.S. Beaty February, 1998

%$Id: killint3.m,v 2.0 1999/05/21 18:47:11 mah Exp gilles $
%$Log: killint3.m,v $
%Revision 2.0  1999/05/21 18:47:11  mah
%Release 2
%
%Revision 1.1  1999/01/06 19:09:25  kay
%Initial revision
%
%
%Copyright (C) 1998 Seismology and Electromagnetic Section/
%Continental Geosciences Division/Geological Survey of Canada
%
%This library is free software; you can redistribute it and/or
%modify it under the terms of the GNU Library General Public
%License as published by the Free Software Foundation; either
%version 2 of the License, or (at your option) any later version.
%
%This library is distributed in the hope that it will be useful,
%but WITHOUT ANY WARRANTY; without even the implied warranty of
%MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
%Library General Public License for more details.
%
%You should have received a copy of the GNU Library General Public
%License along with this library; if not, write to the
%Free Software Foundation, Inc., 59 Temple Place - Suite 330,
%Boston, MA  02111-1307, USA.
%
%DSI Consortium
%Continental Geosciences Division
%Geological Survey of Canada
%615 Booth St.
%Ottawa, Ontario
%K1A 0E9
%
%email: dsi@@cg.nrcan.gc.ca

global ptimes
global temp

data=get(gcf,'userdata');
axhand=data{5};
child=get(axhand,'children');
for n=1:3
 L{n}=flipud(findobj(child{n},'type','line'));
 P{n}=flipud(findobj(child{n},'type','patch'));
end %for
row=length(get(L{1}(1),'xdata'));
tr1=str2num(get(findobj('tag','EditTr1'),'string'));
tr2=str2num(get(findobj('tag','EditTr2'),'string'));
ntr=tr2-tr1+1;

switch action
 case 'firstcall'
  for n=1:3
   set(L{n}(1:ntr),'erasemode','background')
   set(P{n}(1:ntr),'erasemode','background')
   killflg{n}=zeros(row,ntr);
  end %for
  bad=find(thm(6,tr1:tr2)<0);

 case 'killbutton'
  wbdf=get(gcf,'windowbuttondownfcn');
  if ~isempty(wbdf) & length(wbdf)==11
   set(gcf,'windowbuttondownfcn','');
   set(findobj(gcf,'string','PICK'),'enable','on')
  else
   set(gcf,'windowbuttondownfcn','killint3 on')
   set(findobj(gcf,'string','PICK'),'enable','off')
   for n=1:3
    set(L{n}(1:ntr),'erasemode','background')
    set(P{n}(1:ntr),'erasemode','background')
   end %for
  end %if
  return;

 case 'on'
  pt=get(gca,'currentpoint');
  bad=round(pt(1,1))-tr1+1;
  killflg=get(findobj(gcf,'string','KILL'),'userdata');
  deadtr=find(killflg{1}(1,bad)~=0);
  if ~isempty(deadtr) %if trace has already been killed
   for n=1:3
    set(L{n}(bad),'xdata',killflg{n}(:,bad),'color','k')
    set(L{n}(bad),'linestyle','-')
    set(P{n}(bad),'visible','on')
    killflg{n}(:,bad)=0;
   end %for
   set(findobj(gcf,'string','KILL'),'userdata',killflg)
   return;
  end %if
end %case

for n=1:3
 for i=bad
  killflg{n}(:,i)=get(L{n}(i),'xdata')';
  flatline(1:row)=i+tr1-1;
  set(L{n}(i),'xdata',flatline,'color','r');
  set(P{n}(i),'visible','off');
 end %for
end %for
set(findobj(gcf,'string','KILL'),'userdata',killflg)

a=get(temp(3),'ydata');
a(bad)=NaN;
for n=3:5
 set(temp(n),'ydata',a)
end %for
a=find(ptimes(1,bad)~=0);
if ~isempty(a)
 b=find(ptimes(2,bad)~=0);
 if ~isempty(b)
  delete(ptimes(2:4,bad(b)))
 end %if
 ptimes(1:4,bad(a))=0;
end %if




@


2.0
log
@Release 2
@
text
@d13 1
a13 1
%$Id: killint3.m,v 1.1 1999/01/06 19:09:25 kay Exp mah $
d15 3
@


1.1
log
@Initial revision
@
text
@d13 4
a16 2
%$Id: $
%$Log:$
d18 1
@