Calculate the expected multitaper admittance and correlation spectra associated with the input global cross-power spectra of two functions.
call SHAdmitCorr (sgt, sgg, stt, lmax, tapers, lwin, k, admit, corr, mtdef, taper_wt)
sgt : input, real*8, dimension (lmax+1)
: The global cross-power spectrum of the functions G and T.
sgg : input, real*8, dimension (lmax+1)
: The global power spectrum of the function G.
stt : input, real*8, dimension (lmax+1)
: The global power spectrum of the function T.
lmax : input, integer
: The maximum spherical harmonic degree to use of the input spectra sgt, sgg, stt.
tapers : input, real*8, dimension (lwin+1, k)
: The spherical harmonic coefficients of the spherical cap localizing windows. Each column corresponds to the non-zero coefficients of a single angular order. Since all that is important is the power spectrum of each window, the exact angular order is not important. These are generated by a call to SHReturnTapers or SHReturnTapersM.
lwin : input, integer
: The spherical harmonic bandwidth of the localizing windows.
k : input, integer
: The number of localizing windows to use. Only the first k columns of tapers will be employed, which corresponds to the k best-concentrated localizing windows.
admit : output, real*8, dimension (lmax-lwin+1)
: The biased admittance spectrum obtained using the localized (cross-)power spectra of Sgt and Stt.
corr : output, real*8, dimension (lmax-lwin+1)
: The biased correlation spectrum obtained using the localized (cross-)power spectra of Sgt, Stt, and Sgg.
mtdef : optional, input, integer, default = 1
: (1) Calculate the biased multitaper cross-power spectra first, and use these to calculate the final admittance and correlation spectra. 2: Calculate an admittance and correlation spectrum using each individual taper, and then average these to get the final admittance and correlation spectra.
taper_wt : input, optional, real*8, dimension (k)
: The weights to apply to each individual windowed spectral estimate. The weights must sum to unity and are obtained from SHMTVarOpt.
Given the global cross-power spectra Sgt, Sgg and Stt of two functions G and T, SHBiasAdmitCorr will calculate the expected multitaper admittance and correlation spectra associated with the two global functions. This routine expects as input a matrix containing the spherical harmonic coefficients of the localizing windows, which can be generated by a call to SHReturnTapers or SHReturnTapersM. Only the k best-concentrated localization windows will be employed when calculating the biased cross-power spectra. The maximum calculated degree of the output biased admittance and correlation spectra corresponds to lmax-lwin, as it is assumed that the input cross-power spectra beyond lmax are unknown, and not zero.
Two manners of calculating the localized admittance and correlation spectra are possible according to the value of the optional parameter mtdef. In case 1, the multitaper cross-power spectra of Sgt, Sgg, and Stt are first calculated, and from these, the admittance and correlation spectra are formed. In case 2, the biased admittance and correlation spectra are calculated for each individual taper, and these are then averaged to obtain the biased multitaper admittance and correlation spectra.
The default is to apply equal weights to each individual windowed estimate of the spectrum, but this can be modified when mtdef is 1 by specifying the weights in the optional argument taper_wt. The weights must sum to unity and can be calculated by SHMTVarOpt.
Wieczorek, M. A. and F. J. Simons, Minimum-variance multitaper spectral estimation on the sphere, J. Fourier Anal. Appl., 13, doi:10.1007/s00041-006-6904-1, 665-692, 2007.
Simons, F. J., F. A. Dahlen and M. A. Wieczorek, Spatiospectral concentration on the sphere, SIAM Review, 48, 504-536, doi:10.1137/S0036144504445765, 2006.
Wieczorek, M. A. and F. J. Simons, Localized spectral analysis on the sphere, Geophys. J. Int., 162, 655-675.
shbias, shbiask, shreturntapers, shreturntapersm, shmtvaropt