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
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
def rpp(self):
rpp = reflectivity_array(self.vp, self.vs, self.rho, self.theta)
return rpp
def rpp(self):
rpp = reflectivity_array(self.vp, self.vs, self.rho, self.theta)
return rpp