def example2(): # generate random spherical harmonics coefficients with a given # power spectrum and plot its bandspectrum degrees = np.arange(201) scale = 10 power = 1. / (1. + (degrees / scale)**2)**2 coeffs = SHCoeffs.from_random(power) coeffs.plot_powerperband(show=False, fname='power.png') # expand coefficients into two different spatial grids and plot it grid1 = coeffs.expand(kind='GLQ') grid1.plot_rawdata(show=False, fname='GLQGrid.png') grid2 = coeffs.expand(kind='DH1') grid2.plot_rawdata(show=False, fname='DHGrid.png')
def example2(): # generate random spherical harmonics coefficients with a given # power spectrum and plot its bandspectrum degrees = np.arange(201) scale = 10 power = 1. / (1. + (degrees / scale)**2)**2 coeffs = SHCoeffs.from_random(power) coeffs.plot_powerperband(show=False, fname='power.png') # expand coefficients into two different spatial grids and plot it grid1 = coeffs.expand(kind='GLQ') grid1.plot_rawdata(show=False, fname='GLQGrid.png') grid2 = coeffs.expand(kind='DH1') grid2.plot_rawdata(show=False,fname='DHGrid.png')
def example1(): # generate random spherical harmonics coefficients with a given # power spectrum and plot its bandspectrum degrees = np.arange(201) scale = 10 power = 1. / (1. + (degrees / scale)**2)**2 coeffs = SHCoeffs.from_random(power) coeffs.plot_powerperband(show=False, fname='power.png') # expand coefficients into a spatial grid and plot it grid1 = coeffs.expand(kind='DH1') grid1.plot_rawdata(show=False, fname='DHGrid_unrotated.png') # rotate coefficients, expand to grid and plot again coeffs.rotate(40., 0., 0., degrees=True) grid2 = coeffs.expand(kind='DH1') grid2.plot_rawdata(show=False, fname='DHGrid_rotated.png')