velocity_conversion.py 文件源码

def gaussian_filter(self,sigma_x=0.0,sigma_y=0.0):
      '''
      Applies a gaussian filter to the seismic velocity field to mimic
      the loss of spatial resolution introduced in tomographic imaging
      '''
      from scipy.ndimage.filters import gaussian_filter

      #filter absolute perturbations
      dvp_filtered  = gaussian_filter(self.dvp_abs,sigma=[sigma_x,sigma_y])
      dvs_filtered  = gaussian_filter(self.dvs_abs,sigma=[sigma_x,sigma_y])
      drho_filtered = gaussian_filter(self.drho_abs,sigma=[sigma_x,sigma_y])
      self.dvp_abs  = dvp_filtered
      self.dvs_abs  = dvs_filtered
      self.drho_abs = drho_filtered

      #filter relative perturbations
      dvp_filtered  = gaussian_filter(self.dvp_rel,sigma=[sigma_x,sigma_y])
      dvs_filtered  = gaussian_filter(self.dvs_rel,sigma=[sigma_x,sigma_y])
      drho_filtered = gaussian_filter(self.drho_rel,sigma=[sigma_x,sigma_y])
      self.dvp_rel  = dvp_filtered
      self.dvs_rel  = dvs_filtered
      self.drho_rel = drho_filtered