示例#1
0
    def rpp_t(self, dt):

        if self.domain == 'time':
            if dt == self.dz:
                return self.rpp
            else:
                raise Exception('sampling mismatch')

        vpt = depth_to_time(self.vp, self.vp, self.dz, dt)
        times = np.arange(vpt.shape[0]) * dt
        vst = depth_to_time(self.vs, self.vp, self.dz, times)
        rhot = depth_to_time(self.rho, self.vp, self.dz, times)

        rpp = reflectivity_array(vpt, vst, rhot, self.theta)
        return rpp
示例#2
0
    def rpp_t(self, dt):

        if self.domain == 'time':
            if dt == self.dz:
                return self.rpp
            else:
                raise Exception('sampling mismatch')
            
        
        vpt = depth_to_time(self.vp, self.vp, self.dz, dt)
        times = np.arange(vpt.shape[0]) * dt
        vst = depth_to_time(self.vs, self.vp, self.dz, times)
        rhot = depth_to_time(self.rho, self.vp, self.dz, times)

        rpp = reflectivity_array(vpt, vst, rhot, self.theta)
        return rpp
示例#3
0
    def rpp(self):

        rpp = reflectivity_array(self.vp, self.vs, self.rho, self.theta)

        return rpp
示例#4
0
    def rpp(self):

        rpp = reflectivity_array(self.vp, self.vs, self.rho, self.theta)

        return rpp