def rotate(self, rotation_angle): """ use mtpy.analysis.geometry to rotate a z array and recalculate res and phase """ import pek1dclasses as pek1dc if not hasattr(self, 'z'): self.read_respfile() new_z = np.zeros_like(self.z) zerr = np.zeros_like(self.z, dtype=float) for iz, zarray in enumerate(self.z): new_z[iz], ze = MTg.MTz.rotate_z(zarray, rotation_angle) self.z = new_z self.resistivity = np.zeros_like(self.z, dtype=float) self.phase = np.zeros_like(self.z, dtype=float) for iz in range(len(self.z)): r, re, p, pe = pek1dc._compute_res_phase(self.z[iz], zerr[iz], self.freq) self.resistivity[iz] = r # self.resistivity_err[iz] = re self.phase[iz] = p # self.phase_err[iz] = pe self.rotation_angle = rotation_angle
def rotate(self,rotation_angle): """ use mtpy.analysis.geometry to rotate a z array and recalculate res and phase """ import pek1dclasses as pek1dc if not hasattr(self,'z'): self.read_respfile() new_z = np.zeros_like(self.z) zerr = np.zeros_like(self.z,dtype=float) for iz,zarray in enumerate(self.z): new_z[iz],ze = MTg.MTz.rotate_z(zarray,rotation_angle) self.z = new_z self.resistivity = np.zeros_like(self.z,dtype=float) self.phase = np.zeros_like(self.z,dtype=float) for iz in range(len(self.z)): r,re,p,pe = pek1dc._compute_res_phase(self.z[iz],zerr[iz],self.freq) self.resistivity[iz] = r # self.resistivity_err[iz] = re self.phase[iz] = p # self.phase_err[iz] = pe self.rotation_angle = rotation_angle
def rotate(self, rotation_angle): """ use mtpy.analysis.geometry to rotate a z array and recalculate res and phase """ import pek1dclasses as pek1dc if not hasattr(self, 'z'): self.read_datafile() new_z = np.zeros_like(self.z) new_ze = np.zeros_like(self.zerr, dtype=float) # for iz,zarray in enumerate(self.z): new_z, new_ze = MTg.MTz.rotate_z(self.z, rotation_angle, zerr_array=self.zerr) self.z = new_z self.zerr = new_ze self.resistivity, self.resistivity_err, self.phase, self.phase_err = \ pek1dc._compute_res_phase(self.z,self.zerr,self.freq) self.rotation_angle = rotation_angle
def rotate(self,rotation_angle): """ use mtpy.analysis.geometry to rotate a z array and recalculate res and phase """ import pek1dclasses as pek1dc if not hasattr(self,'z'): self.read_datafile() new_z = np.zeros_like(self.z) new_ze = np.zeros_like(self.zerr,dtype=float) # for iz,zarray in enumerate(self.z): new_z,new_ze = MTg.MTz.rotate_z(self.z,rotation_angle,zerr_array = self.zerr) self.z = new_z self.zerr = new_ze self.resistivity, self.resistivity_err, self.phase, self.phase_err = \ pek1dc._compute_res_phase(self.z,self.zerr,self.freq) self.rotation_angle = rotation_angle