Changeset 1157


Ignore:
Timestamp:
04.02.2016 22:48:32 (4 years ago)
Author:
klaus
Message:

determine slowness of phase picks

Location:
SHX/trunk/SeismicHandler
Files:
1 added
3 edited

Legend:

Unmodified
Added
Removed
  • SHX/trunk/SeismicHandler/basics/analysispar.py

    r1155 r1157  
    200200    def setDefaultWeight( self, weight ): 
    201201        self.default_weight = weight 
     202     
     203    def getPhaseListsByName( self, picktype=(PICKTYPE_MANUAL,PICKTYPE_AUTO), 
     204        minweight=None ): 
     205        "Re-sort phase lists into dictionary with phase name as key items." 
     206        pdict = {} 
     207        for station in self.getStations(): 
     208            for phase in self.getPhaseList(station): 
     209                if picktype and phase.picktype not in picktype: 
     210                    continue 
     211                if minweight != None and phase.weight < minweight: 
     212                    continue 
     213                if phase.name in pdict.keys(): 
     214                    pdict[phase.name].append( (station,phase.picktime) ) 
     215                else: 
     216                    pdict[phase.name] = [(station,phase.picktime)] 
     217        return pdict 
    202218 
    203219 
  • SHX/trunk/SeismicHandler/commands/beam.py

    r1074 r1157  
    152152     
    153153    def absToRelLocations( self, traces ): 
    154         "Returns relative locations of stations." 
     154        """ 
     155        Should use obspy utlGeoKm. !!! 
     156        Returns relative locations of stations. 
     157        """ 
    155158        print "dbg: compute rel from abs" 
    156159        stations = Stations() 
  • SHX/trunk/SeismicHandler/commands/phase.py

    r1150 r1157  
    99from SeismicHandler.basics.analysispar import PhaseList, AnalysisPar 
    1010from SeismicHandler.basics.error import ShxError 
     11from SeismicHandler.modules.seismics import getPhaseSlownesses 
    1112from SeismicHandler.basics import timeit 
    1213from SeismicHandler.modules.traces import addNetcode, number_of_visible_traces,\ 
     
    3233    PHASE DUMP 
    3334    PHASE PLOT [<outfile>] 
     35    PHASE SLOWNESS <phase> [<slovar>] [<azvar>] [<sloerrvar>] [<azerrvar>] 
    3436    """ 
    3537    numberOfParameters = [1,2,3,4,5,6] 
     
    175177                dspprog = 'display' 
    176178            os.system( "%s %s &" % (dspprog,outfile) ) 
     179        elif subcmd == 'slowness': 
     180            pdict = getPhaseSlownesses() 
     181            print "dbg:", pdict 
     182            pname = self.parameters[1] 
     183            if pname not in pdict.keys(): 
     184                return 
     185            if len(self.parameters) > 2 and self.parameters[2]: 
     186                varname = self.parameters[2][1:] 
     187                self.symbols.set( varname, pdict[pname][0] ) 
     188            if len(self.parameters) > 3 and self.parameters[3]: 
     189                varname = self.parameters[3][1:] 
     190                self.symbols.set( varname, pdict[pname][1] ) 
     191            if len(self.parameters) > 4 and self.parameters[4]: 
     192                varname = self.parameters[4][1:] 
     193                self.symbols.set( varname, pdict[pname][2] ) 
     194            if len(self.parameters) > 5 and self.parameters[5]: 
     195                varname = self.parameters[5][1:] 
     196                self.symbols.set( varname, pdict[pname][3] ) 
    177197        else: 
    178198            raise ShxError( "Command 'phase', unknown subcommand '%s'" % subcmd, 
Note: See TracChangeset for help on using the changeset viewer.